Elastic properties and masticatory bone stress in the macaque mandible

One important limitation of mechanical analyses with strain gages is the difficulty in directly estimating patterns of stress or loading in skeletal elements from strain measurements. Because of the inherent anisotropy in cortical bone, orientation of principal strains and stresses do not necessarily coincide, and it has been demonstrated theoretically that such differences may be as great as 45 degrees (Cowin and Hart, 1990). Likewise, relative proportions of stress and strain magnitudes may differ. This investigation measured the elastic properties of a region of cortical bone on both the buccal and lingual surfaces of the lower border of the macaque mandible. The elastic property data was then combined with macaque mandibular strain data from published and a new in vivo strain gage experiment to determine directions and magnitudes of maximum and minimum principal stresses. The goal was to compare the stresses and strains and assess the differences in orientation and relative magnitude between them. The main question was whether these differences might lead to different interpretations of mandibular function. Elastic and shear moduli, and Poisson's ratios were measured using an ultrasonic technique from buccal and lingual cortical surfaces in 12 macaque mandibles. Mandibular strain gage data were taken from a published set of experiments (Hylander, 1979), and from a new experiment in which rosette strain gauges were fixed to the buccal and lingual cortices of the mandibular corpus of an adult female Macaca fascicularis, after which bone strain was recorded during mastication. Averaged elastic properties were combined with strain data to calculate an estimate of stresses in the mandibular corpus. The elastic properties were similar to those of the human mandibular cortex. Near its lower border, the macaque mandible was most stiff in a longitudinal direction, less stiff in an inferosuperior direction, and least stiff in a direction normal to the bone's surface. The lingual aspect of the mandible was slightly stiffer than the buccal aspect. Magnitudes of stresses calculated from average strains ranged from a compressive stress of -16.00 GPa to a tensile stress of 8.84 GPa. The orientation of the principal stresses depended on whether the strain gage site was on the working or balancing side. On the balancing side of the mandibles, maximum principal stresses were oriented nearly perpendicular to the lower border of the mandible. On the working side of the mandibles, the orientation of the maximum principal stresses was more variable than on the balancing side, indicating a larger range of possible mechanisms of loading. Near the lower border of the mandible, differences between the orientation of stresses and strains were 12 degrees or less. Compared to ratios between maximum and minimum strains, ratios between maximum and minimum stresses were more divergent from a ratio of 1.0. Results did not provide any major reinterpretations of mandibular function in macaques, but rather confirmed and extended existing work. The differences between stresses and strains on the balancing side of the mandible generally supported the view that during the power stroke the mandible was bent and slightly twisted both during mastication and transducer biting. The calculated stresses served to de-emphasize the relative importance of torsion. On the working side, the greater range of variability in the stress analysis compared to the strain analysis suggested that a more detailed examination of loadings and stress patterns in each individual experiment would be useful to interpret the results. Torsion was evident on the working side; but in a number of experiments, further information was needed to interpret other superimposed regional loading patterns, which may have included parasagittal bending and reverse parasagittal bending.     

Effect of superficial collagen patterns and fibrillation of femoral articular cartilage on knee joint mechanics-a 3D finite element analysis

Collagen fibrils of articular cartilage have specific depth-dependent orientations and the fibrils bend in the cartilage surface to exhibit split-lines. Fibrillation of superficial collagen takes place in osteoarthritis. We aimed to investigate the effect of superficial collagen fibril patterns and collagen fibrillation of cartilage on stresses and strains within a knee joint. A 3D finite element model of a knee joint with cartilage and menisci was constructed based on magnetic resonance imaging. The fibril-reinforced poroviscoelastic material properties with depth-dependent collagen orientations and split-line patterns were included in the model. The effects of joint loading on stresses and strains in cartilage with various split-line patterns and medial collagen fibrillation were simulated under axial impact loading of 1000 N. In the model, the collagen fibrils resisted strains along the split-line directions. This increased also stresses along the split-lines. On the contrary, contact and pore pressures were not affected by split-line patterns. Simulated medial osteoarthritis increased tissue strains in both medial and lateral femoral condyles, and contact and pore pressures in the lateral femoral condyle. This study highlights the importance of the collagen fibril organization, especially that indicated by split-line patterns, for the weight-bearing properties of articular cartilage. Osteoarthritic changes of cartilage in the medial femoral condyle created a possible failure point in the lateral femoral condyle. This study provides further evidence on the importance of the collagen fibril organization for the optimal function of articular cartilage.     

Masticatory stress, orbital orientation and the evolution of the primate postorbital bar

A postorbital bar is one of a suite of derived features which distinguishes basal primates from their putative sister taxon, plesiadapiforms. Two hypotheses have been put forward to explain postorbital bar development and variation in circumorbital form: the facial torsion model and visual predation hypothesis. To test the facial torsion model, we employ strain data on circumorbital and mandibular loading patterns in representative primates with a postorbital bar and masticatory apparatus similar to basal primates. To examine the visual predation hypothesis, we employ metric data on orbit orientation in Paleocene and Eocene primates, as well as several clades of visual predators and foragers that vary interspecifically in postorbital bar formation.A comparison of galago circumorbital and mandibular peak strains during powerful mastication demonstrates that circumorbital strains are quite low. This indicates that, as in anthropoids, the strepsirhine circumorbital region is excessively overbuilt for countering routine masticatory loads. The fact that circumorbital peak-strain levels are uniformly low in both primate suborders undermines any model which posits that masticatory stresses are determinants of circumorbital form, function and evolution. This is interpreted to mean that sufficient cortical bone must exist to prevent structural failure due to non-masticatory traumatic forces. Preliminary data also indicate that the difference between circumorbital and mandibular strains is greater in larger taxa.Comparative analyses of several extant analogs suggest that the postorbital bar apparently provides rigidity to the lateral orbital margins to ensure a high level of visual acuity during chewing and biting. The origin of the primate postorbital bar is linked to changes in orbital convergence and frontation at smaller sizes due to nocturnal visual predation and increased encephalization. By incorporating in vivo and fossil data, we reformulate the visual predation hypothesis of primate origins and thus offer new insights into major adaptive transformations in the primate skull.     

Split-line patterns of the mandible following masseterectomy in adult and growing monkeys.

Architectural changes of the mandibular ramus following masseterectomy were analyzed by means of the split-line technique in both adult and growing monkeys. The right masseter muscle was resected in ten tufted capuchin monkeys (Cebus apella) and they were sacrificed following a 12-month period. Along with modifications in the shape of the mandible, the masseterectomy induced structural changes as revealed by alterations in the split-line patterns. These changes were much more evident in the adult animals. In growing monkeys, intrinsic forces determine an architecture of mandibular growth, represented by the direction of the splits. This is not influenced by the lack of mechanical factors like muscle action. Conversely, the final fibrous architecture of the bone in the adult monkey may be modified by extrinsic mechanical factors.     

Main patterns and individual differences in baboon skull split-lines and theories of causes of split-line orientation in bone

Split-line patterns are reported in skulls of five adult male baboons. While variations in pattern occur in all parts of the skull, these variations are relatively minor in the following regions: supraorbital, lateral orbital, medial orbital, nasal bones, zygomatico-alveolar crest, nasal opening, alveolar process of maxilla and mandible. Wide differences in pattern occur in these regions: infraorbital, zygomatic bone, body of maxilla, and frontal bone posterior to the supraorbital area. The major variability in split-line orientation indicates that oversimplified interpretations of the patterns in terms of (1) conformity to gross structure, or (2) direction of bone growth, are untenable. The variations do not contradict a functional interpretation in which mechanical forces and skull form interact to different degrees in different individuals, however. Skulls of a variety of primates are useful for functional analysis, because they have similar structural plans, but the differences are well outside the normal range of variation for a single species.     

Three-dimensional finite element stress analysis of the dentate human mandible.

The biomechanical events which accompany functional loading of the human mandible are not fully understood. The techniques normally used to record them are highly invasive. Computer modelling offers a promising alternative approach in this regard, with the additional ability to predict regional stresses and strains in inaccessible locations. In this study, we built two three-dimensional finite element (FE) models of a human mandible reconstructed from tomographs of a dry dentate jaw. The first model was used for a complete mechanical characterization of physical events. It also provided comparative data for the second model, which had an increased vertical corpus depth. In both cases, boundary conditions included rigid restraints at the first right molar and endosteal cortical surfaces of the articular eminences of temporal bones. Groups of parallel multiple vectors simulated individual masticatory muscle loads. The models were solved for displacements, stresses, strains, and forces. The simulated muscle loads in the first model deformed the mandible helically upward and toward its right (working) side. The highest principal stresses occurred at the bite point, anterior aspects of the coronoid processes, symphyseal region, and right and left sides of the mandibular corpus. In general, the observed principal stresses and strains were highest on the periosteal cortical surface and alveolar bone. At the symphyseal region, maximum principal stresses and strains were highest on the lower lingual mandibular aspect, whereas minimum principal stresses and strains were highest on its upper labial side. Subcondylar principal strains and condylar forces were higher on the left (balancing or nonbiting) side than on the right mandibular side, with condylar forces more concentrated on the anteromedial aspect of the working-side condyle and on the central and lateral aspects of the left. When compared with in vivo strain data from macaques during comparable biting events, the predictive strain values from the first model were qualitatively similar. In the second model, the reduced tensile stress on the working-side, and decreased shear stress bilaterally, confirmed that lower stresses occurred on the lower mandibular border with increased jaw depth. Our results suggested that although the mandible behaved in a beam-like manner, its corpus acted more like a combination of open and closed cross sections due to the presence of tooth sockets, at least for the task modelled.(ABSTRACT TRUNCATED AT 400 WORDS)     

Uniform stress state in bone structure with residual stress.

Residual stress and strain in living tissues have been investigated from the viewpoint of mechanical optimality maintained by adaptive remodeling. In this study, the residual stresses in the cortical-cancellous bone complex of bovine coccygeal vertebrae were examined. Biaxial strain gages were bonded onto the cortical surface, so that the gage axes were aligned in the cephalocaudal and circumferential directions. Strains induced by removal of the end-plate and the cancellous bone were recorded sequentially. The results revealed the existence of compressive residual stress in the cortical bone and tensile residual stress in the cancellous bone in both the cephalocaudal and the circumferential direction. The observed strains were examined on the basis of the uniform stress hypothesis using a three-bar model for the cephalocaudal direction and a three-layered cylinder model for the circumferential direction. In this model study, the magnitude of effective stresses, which is defined as the macroscopic stress divided by the area fraction of bone material, was found not to differ significantly between cephalocaudal and circumferential directions, although they were evaluated using independent models. These results demonstrate that the uniform stress state of the cortical-cancellous bone structure is consistent with results obtained in the cutting experiment when the existence of residual stress is taken into account.     

The effects of misalignment during in vivo loading of bone: Techniques to detect the proximity of objects in three-dimensional models

Theories of mechanical adaptation of bone suggest that mechanical loading causes bone formation at discrete locations within bone microstructure experiencing the greatest mechanical stress/strain. Experimental testing of such theories requires in vivo loading experiments and high-resolution finite element models to determine the distribution of mechanical stresses. Finite element models of in vivo loading experiments typically assume idealized boundary conditions with applied load perfectly oriented on the bone, however small misalignments in load orientation during an in vivo experiment are unavoidable, and potentially confound the ability of finite element models to predict locations of bone formation at the scale of micrometers. Here we demonstrate two different three-dimensional spatial correlation methods to determine the effects of misalignment in load orientation on the locations of high mechanical stress/strain in the rodent tail loading model. We find that, in cancellous bone, the locations of tissue with high stress are maintained under reasonable misalignments in load orientation (p<0.01). In cortical bone, however, angular misalignments in the dorsal direction can alter the locations of high mechanical stress, but the locations of tissue with high stress are maintained under other misalignments (p<0.01). We conclude that, when using finite element models of the rodent tail loading model, small misalignments in loading orientation do not affect the predicted locations of high mechanical stress within cancellous bone.     

Variations in cortical material properties throughout the human dentate mandible

Material properties and their variations in individual bone organs are important for understanding bone adaptation and quality at a tissue level, and are essential for accurate mechanical models. Yet material property variations have received little systematic study. Like all other material property studies in individual bone organs, studies of the human mandible are limited by a low number of both specimens and sampled regions. The aims of this study were to determine: 1) regional variability in mandibular material properties, 2) the effect of this variability on the modeling of mandibular function, and 3) the relationship of this variability to mandibular structure and function. We removed 31 samples on both facial and lingual cortices of 10 fresh adult dentate mandibles, measured cortical thickness and density, determined the directions of maximum stiffness with a pulse transmission ultrasonic technique, and calculated elastic properties from measured ultrasonic velocities. Results showed that each of these elastic properties in the dentate human mandible demonstrates unique regional variation. The direction of maximum stiffness was near parallel to the occlusal plane within the corpus. On the facial ramus, the direction of maximum stiffness was more vertically oriented. Several sites in the mandible did not show a consistent direction of maximum stiffness among specimens, although all specimens exhibited significant orthotropy. Mandibular cortical thickness varied significantly (P < 0.001) between sites, and decreased from 3.7 mm (SD = 0.9) anteriorly to 1.4 mm posteriorly (SD = 0.1). The cortical plate was also significantly thicker (P < 0.003) on the facial side than on the lingual side. Bone was 50-100% stiffer in the longitudinal direction (E(3), 20-30 GPa) than in the circumferential or tangential directions (E(2) or E(1); P < 0.001). The results suggest that material properties and directional variations have an important impact on mandibular mechanics. The accuracy of stresses calculated from strains and average material properties varies regionally, depending on variations in the direction of maximum stiffness and anisotropy. Stresses in some parts of the mandible can be more accurately calculated than in other regions. Limited evidence suggests that the orientations and anisotropies of cortical elastic properties correspond with features of cortical bone microstructure, although the relationship with functional stresses and strains is not clear.     

Anisotropy and inter-condyle heterogeneity of cartilage under large-strain shear

We were the first to examine the mechanical responses of skeletally mature bovine femoral cartilage under large-strain simple shear (up to ±20%) using a multiaxial shear testing device. Since shear loading is critical in both tissue failure and chondrocyte responses, we aimed to probe (1) anisotropy with respect to the split-line direction (principal alignment of the collagen fibers near the articulating surface), (2) heterogeneity between femoral condyles, and (3) the influence of local cartilage thickness. We harvested a total of 48 cuboid cartilage specimens from four bovine knee joints. With each specimen we applied shear strains both parallel and perpendicular to the local split-line direction at a rate of 75 μm/min and calculated the peak-to-peak shear stresses, shear strain–energy dissipation densities, and peak effective shear moduli. The Wilcoxon signed rank test revealed that the medial condyle was anisotropic in some mechanical measures at applied shear strains above 5%, while the lateral condyle was mechanically isotropic at all applied shear strains. The Kruskal–Wallis test revealed no significant differences in the median mechanical behavior of the lateral and medial condyles. Spearman׳s rank correlations revealed statistically significant negative monotonic correlations among thickness and most of our mechanical measures for both lateral and medial condyles at most applied strains and directions of applied shear. These results suggest that large-strain analyses account for nonlinear, anisotropic and location-dependent effects not fully realized at small strains. Our findings may inspire new experiments and models that consider anisotropy and heterogeneity of cartilage in ways previously ignored.     

Bone structure and the patterns of force transmission in the cat skull (Felis catus)

Projection microradiography was used to determine the density and orientation of the force transmitting structures, i.e., trabeculae and bone lying between approximately parallel vascular canals, within the bones of cat skulls. The organisation in the skulls was confirmed statistically for a total of ten cats. The results of the observations showed that within specific areas of the skull a high degree of structural orientation and an increased density of osseous structures was present. The distribution of these characters corresponded in contiguous bones such that a continuum of structural organisation was established between the alveolar region and the site of attachment of the temporalis and masseter muscles and the glenoid region. The patterns of force transmission during jaw closure were determined when a resistance was placed initially between the canines and then the carnassials. An analysis was first carried out on dry skulls using colophonium resin to determine the direction of the force distribution. The nature and the approximate magnitude of the forces were ascertained by replacing the resin with strain gauges. The basic similarities in the strain patterns recorded from the dry skulls and those from the ten anaesthetised cats in which strain gauges had been intra-vitally implanted, substantiated the recordings made on the dry skulls. Combination of the results from the three sets of experiments defined the patterns of force distribution in the cat skull during the closure of the mandible against a resistance. The results showed that: (1) the combined action of the temporalis and masseter muscles tended to reduce the overall strain in the skull bones, and that the deformations produced by the action of the masseter were greater than that exerted by the temporalis muscles; (2) during biting, whether the resistance was placed between the canines or carnassials, compressive forces predominated in the facial bones; (3) small movements observed between facial bones indicated the presence of a flexible component within the skull, thus allowing large forces to be exerted during biting without overstressing the facial bones; (4) the glenoid fossa is part of a force bearing joint; (5) forces generated during biting were resisted within the skull by forces of an opposite nature generated within the system, the incompressible nature of bone and by the effect of the soft tissues; (6) the nature and the magnitude of the strain altered when a resistance was placed at the canines and then at the carnassials; however, the pattern of force distribution within the skull remained the same; (7) there was a direct correspondence between the detailed structural organisation of the bones and the patterns of force distribution. This conclusion would appear to apply in general to mammalian skulls. The study also emphasises the importance, neglected hitherto, of carrying out a variety of experiments to determine the patterns of force distribution in bones. The Trajectorial Theory of bone organisation is discussed and, on the basis of the results obtained, a modified theory is proposed. This states that: the structural continuum is common to the compact and cancellous bone and comprises bony bars which are aligned in the optimum direction for the transmission of force to a region in the bone or bones where it is effectively resisted.     

Biaxial constitutive relations for the passive canine diaphragm.

Samples of the muscular sheet excised from the midcostal region of dog diaphragms were subjected to biaxial loading. That is, stresses in the direction of the muscle fibers and in the direction perpendicular to the fibers in the plane of the sheet were measured at different combinations of strains in the two directions. Stress-strain relations were obtained by fitting equations to these data. In the direction of the muscle fibers, for strains up to 0.7, stress is a modestly nonlinear function of strain and ranges up to approximately 60 g/cm. In the direction perpendicular to the fibers, the sheet is stiffer and more strongly nonlinear. At a strain in the perpendicular direction of approximately 0.35, stress increases abruptly. The stress-strain relation in the muscle direction is consistent with observations of passive muscle shortening in vivo. However, the stiffness in the perpendicular direction is not high enough to explain the observation that strains in the perpendicular direction in vivo are nearly zero. We conclude that, in the passive diaphragm in vivo, stress in the direction perpendicular to the muscle fibers is small.     

A comparison of stress distributions for different surgical procedures, screw dimensions and orientations for a Temporomandibular joint implant

Finite element analysis is a useful analytical tool for the design of biomedical implants. The aim of this study was to investigate the behavior of temporomandibular joint implants with multiple design variables of the screws used for fixation of the implant. A commercially available implant with full mandible was analyzed using a finite element software package. The effects of different design variables such as orientation, diameter and stem length of the screws on the stress distribution in bone for two different surgical procedures were investigated. Considering the microstrain in bone as a principal factor, the acceptable ranges for screw diameter and length were determined. Parallel orientation of the screws performed better from a stress point of view when compared to the zig-zag orientation. Sufficient contact between the implant collar and mandibular condyle was shown to reduce the peak stresses which may lead to long term success. The distance between screw holes in the parallel orientation was much closer when compared to the zig-zag orientation. However, the stresses in bone near the screw hole area for the parallel orientation were within acceptable limits.     

Internal strains and stresses measured in cortical bone via high-energy X-ray diffraction

High-energy synchrotron X-ray diffraction was used to study internal stresses in bone under in situ compressive loading. A transverse cross-section of a 12-14 year old beagle fibula was studied with 80.7 keV radiation, and the transmission geometry was used to quantify internal strains and corresponding stresses in the mineral phase, carbonated hydroxyapatite. The diffraction patterns agreed with tabulated patterns, and the distribution of diffracted intensity around 00.2/00.4 and 22.2 diffraction rings was consistent with the imperfect 00.1 fiber texture expected along the axis of a long bone. Residual compressive stress along the bone's longitudinal axis was observed in the specimen prior to testing: for 22.2 this stress equaled -95 MPa and for 00.2/00.4 was between -160 and -240 MPa. Diffraction patterns were collected for applied compressive stresses up to -110 MPa, and, up to about -100 MPa, internal stresses rose proportionally with applied stress but at a higher rate, corresponding to stress concentration in the mineral of 2.8 times the stress applied. The widths of the 00.2 and 00.4 diffraction peaks indicated that crystallite size perpendicular to the 00.1 planes increased from t=41 nm before stress was applied to t=44 nm at -118 MPa applied stress and that rms strain epsilon(rms) rose from 2200 muepsilon before loading to 4600 muepsilon at the maximum applied stress. Small angle X-ray scattering of the unloaded sample, recorded after deformation was complete, showed a collagen D-period of 66.4 nm (along the bone axis).     

Regulation of bone stress and strain in the immature and mature rat femur

Bone in vivo stresses and moments were determined from rosette strain recordings obtained from the mid-diaphysis of growing exercising rats. Two activity groups were examined beginning at 3 weeks of age: 2 min day-1 and 45 min day-1 at 0.2 ms-1 in an exercise wheel. In vitro moment-strain curves were obtained during mechanical calibration tests on intact femora, and area inertial properties were determined from the mid-diaphysis cross-sections. The mechanical calibration and histomorphometry procedures were then used to compute functional stresses and moments based on the in vivo rosette strain recordings. During the period 6-30 weeks of age the rats increased in body weight over threefold, but no significant changes in principal strain and stress magnitude or orientation were found. Peak in vivo compressive and tensile moments increased during growth in proportion to the animal mass squared, but the ratio of these moments to animal body weight times bone length (BWBL) remained constant throughout growth and in the adult. The parameter BWBL appears, therefore, to be a useful predictor of long bone functional strength. Peak torsional moments remained a constant 8.1 +/- 3.0% of the ultimate torsional strength, providing a safety factor of approximately 12. Differences in the in vivo moments between the two activity groups were found, which were due primarily to adaptive, but not significant, changes in bone geometry. These findings support the hypothesis that long bones model and model during growth and altered activity in order to regulate the functional strains at a predefined level.     

The effects of the periodontal ligament on mandibular stiffness: a study combining finite element analysis and geometric morphometrics

It is generally accepted that the periodontal ligament (PDL) plays a crucial role in transferring occlusal forces from the teeth to the alveolar bone. Studies using finite element analysis (FEA) have helped to better understand this role and show that the stresses and strains in the alveolar bone are influenced by whether and how PDL is included in FE models. However, when the overall distribution of stresses and strains in crania and mandibles are of interest, PDL is often not included in FE models, although little is known about how this affects the results. Here we study the effect of representing PDL as a layer of solid material with isotropic homogeneous properties in an FE model of a human mandible using a novel application of geometric morphometrics. The results show that the modelling of the PDL affects the deformation and thus strain magnitudes not only of the alveolar bone around the biting tooth, but that the whole mandible deforms differently under load. As a result, the strain in the mandibular corpus is significantly increased when PDL is included, while the strain in the bone beneath the biting tooth is reduced. These results indicate the importance of the PDL in FE studies. Thus we recommend that the PDL should be included in FE models of the masticatory apparatus, with tests to assess the sensitivity of the results to changes in the Young's modulus of the PDL material.     

Innovative Approaches to the Relationship Between Diet and Mandibular Morphology in Primates

Attempts to establish relationships between mandibular morphology and either traditional dietary categories or geometric and material properties of primate diets have not been particularly successful. Using our conceptual framework of the feeding factors impacting mandibular morphology, we argue that this is because dietary categories and food geometric and material properties affect mandibular morphology only through intervening variables that are currently poorly understood, i.e., feeding behavior, mandibular loading, and stress and strain regimes. Our studies of 3-dimensional jaw kinematics in macaques and capuchins show that, although jaw movement profiles during chewing are affected by food material properties and species-level effects, patterns of jaw movements in these two species are broadly similar. However, because mandibular loading, stress, and strain regimes are determined by interactions between feeding behavior (such as jaw kinematics) and mandibular morphology, it is difficult to say whether these similarities in chewing kinematics also mean similarities in loading, stress, and strain. Comparative analyses of the scaling of daily feeding time and chew cycle duration reveal only weak support for the hypothesis that larger primates chew more than smaller primates. Consideration of these results suggests that better data are needed on the relationship between dietary categories, food material and geometric properties, the amount of time/cycles associated with different feeding behaviors (ingestion, premolar biting, mastication), and mandible stress and strain patterns if we are to understand fully relationships between mandibular morphology and diet in primates.     

In vivo tibial stiffness is maintained by whole bone morphology and cross-sectional geometry in growing female mice

Whole bone morphology, cortical geometry, and tissue material properties modulate skeletal stresses and strains that in turn influence skeletal physiology and remodeling. Understanding how bone stiffness, the relationship between applied load and tissue strain, is regulated by developmental changes in bone structure and tissue material properties is important in implementing biophysical strategies for promoting healthy bone growth and preventing bone loss. The goal of this study was to relate developmental patterns of in vivo whole bone stiffness to whole bone morphology, cross-sectional geometry, and tissue properties using a mouse axial loading model. We measured in vivo tibial stiffness in three age groups (6, 10, 16 wk old) of female C57Bl/6 mice during cyclic tibial compression. Tibial stiffness was then related to cortical geometry, longitudinal bone curvature, and tissue mineral density using microcomputed tomography (microCT). Tibial stiffness and the stresses induced by axial compression were generally maintained from 6 to 16 wks of age. Growth-related increases in cortical cross-sectional geometry and longitudinal bone curvature had counteracting effects on induced bone stresses and, therefore, maintained tibial stiffness similarly with growth. Tissue mineral density increased slightly from 6 to 16 wks of age, and although the effects of this increase on tibial stiffness were not directly measured, its role in the modulation of whole bone stiffness was likely minor over the age range examined. Thus, whole bone morphology, as characterized by longitudinal curvature, along with cortical geometry, plays an important role in modulating bone stiffness during development and should be considered when evaluating and designing in vivo loading studies and biophysical skeletal therapies.     

Osteocyte lacunae tissue strain in cortical bone

Current theories suggest that bone modeling and remodeling are controlled at the cellular level through signals mediated by osteocytes. However, the specific signals to which bone cells respond are still unknown. Two primary theories are: (1) osteocytes are stimulated via the mechanical deformation of the perilacunar bone matrix and (2) osteocytes are stimulated via fluid flow generated shear stresses acting on osteocyte cell processes within canaliculi. Recently, much focus has been placed on fluid flow theories since in vitro experiments have shown that bone cells are more responsive to analytically estimated levels of fluid shear stress than to direct mechanical stretching using macroscopic strain levels measured on bone in vivo. However, due to the complex microstructural organization of bone, local perilacunar bone tissue strains potentially acting on osteocytes cannot be reliably estimated from macroscopic bone strain measurements. Thus, the objective of this study was to quantify local perilacunar bone matrix strains due to macroscopically applied bone strains similar in magnitude to those that occur in vivo. Using a digital image correlation strain measurement technique, experimentally measured bone matrix strains around osteocyte lacunae resulting from macroscopic strains of approximately 2000 microstrain are significantly greater than macroscopic strain on average and can reach peak levels of over 30,000 microstrain locally. Average strain concentration factors ranged from 1.1 to 3.8, which is consistent with analytical and numerical estimates. This information should lead to a better understanding of how bone cells are affected by whole bone functional loading.     

Evaluation of threshold stress for bone resorption around screws based on in vivo strain measurement of miniplate

The purpose of this study is to investigate the critical threshold stress causing bone resorption evaluated from strain measurement in vivo, comparing the various finite element models. In this study strains of miniplates used for mandibular fractures were measured once a week until the strains reduced. The maximum bite force for each patient was applied in the incisal, right molar and left molar region. The strains increased and reached a peak level at 2-4 weeks, whereas the bite forces increased during the period of measurements. A 3-D osteosynthesis model using finite element method showed that the compressive stresses of the bone surrounding screws ranged within approximately -40 MPa under the condition generating the same amounts of strains measured in the miniplates. Furthermore, various finite element models simulating mandibular reconstruction using the fibular graft were constructed. The models for reconstruction using single strut fibula showed distinct stress concentration in the cortical bone surrounding screws, and the peak stress levels were 2 to 3 times as strong as that of the fracture model. We conclude that critical threshold for bone resorption should be approximately -50 MPa (3600 micro strain).