Patterns of strain in the macaque tibia during functional activity.

The strain environment of the tibial midshaft of two female macaques was evaluated through in vivo bone strain experiments using three rosette gauges around the circumference of the bones. Strains were collected for a total of 123 walking and galloping steps as well as several climbing cycles. Principal strains and the angle of the maximum (tensile) principal strain with the long axis of the bone were calculated for each gauge site. In addition, the normal strain distribution throughout the cross section was determined from the longitudinal normal strains (strains in the direction of the long axis of the bone) at each of the three gauge sites, and at the corresponding cross-sectional geometry of the bone. This strain distribution was compared with the cross-sectional properties (area moments) of the midshaft. For both animals, the predominant loading regime was found to be bending about an oblique axis running from anterolateral to posteromedial. The anterior and part of the medial cortex are in tension; the posterior and part of the lateral cortex are in compression. The axis of bending does not coincide with the maximum principal axis of the cross section, which runs mediolaterally. The bones are not especially buttressed in the plane of bending, but offer the greatest strength anteroposteriorly. The cross-sectional geometry therefore does not minimize strain or bone tissue. Peak tibial strains are slightly higher than the peak ulnar strains reported earlier for the same animals (Demes et al. [1998] Am J Phys Anthropol 106:87-100). Peak strains for both the tibia and the ulna are moderate in comparison to strains recorded during walking and galloping activities in nonprimate mammals.     

In vivo bone strain and bone functional adaptation

Mechanistic interpretations of bone cross-sectional shapes are based on the paradigm of shape optimization such that bone offers maximum mechanical resistance with a minimum of material. Recent in vivo strain studies (Demes et al., Am J Phys Anthropol 106 (1998) 87-100, Am J Phys Anthropol 116 (2001) 257-265; Lieberman et al., Am J Phys Anthropol 123 (2004) 156-171) have questioned these interpretations by demonstrating that long bones diaphyses are not necessarily bent in planes in which they offer maximum resistance to bending. Potential limitations of these in vivo studies have been pointed out by Ruff et al. (Am J Phys Anthropol 129 (2006) 484-498). It is demonstrated here that two loading scenarios, asymmetric bending and buckling, would indeed not lead to correct predictions of loads from strain. It is also shown that buckling is of limited relevance for many primate long bones. This challenges a widely held view that circular bone cross sections make loading directions unpredictable for bones which is based on a buckling load model. Asymmetric bending is a potentially confounding factor for bones with directional differences in principal area moments (I(max) > I(min)). Mathematical corrections are available and should be applied to determine the bending axis in such cases. It is concluded that loads can be reliably extrapolated from strains. More strain studies are needed to improve our understanding of the relationships between activities, bone loading regimes associated with them, and the cross-sectional geometry of bones.     

Major and minor centroidal axes serve as objective, automatable reference points to test mechanobiological hypotheses using histomorphometry

Recent studies show that minor and major centroidal axes (CA) of long bone cross sections provide valuable predictions of prevailing loading patterns in age and treatment matched cohorts of animals. Furthermore, using CA, we recently showed that the degree of mineralization and area of woven bone laid down in the first two weeks after creation of a critical sized bone defect relate inversely and correlate significantly to loading patterns. Here, we aim to determine how closely independent measures of the spatial distribution of bone apposition determined using the major and minor CA as reference points correlate to those using anatomically defined axes as reference points. In histological sections from the previous study, we found no statistically significant difference between the anatomical and centroidal axes with respect to the centroid location or axis rotation, but there is a significant albeit small difference in the average distance between centroids. Outcome measures calculated in areas of bone defined by 15°, 30°, 45°, 60°, or 90° sectors when using the CA differ less than 5% from those calculated using anatomical axes as reference points. Hence, the major and minor CA provide objective reference points for comparison of mechanobiological outcome measures between animals in matched cohorts. Calculation of major and minor CA is automated, which reduces the potential for observer bias. A major advantage of using the major and minor CA as reference points is that it allows for direct relation of outcome measures to loading patterns in age and treatment matched cohorts, ultimately providing a tool to test mechanobiological hypotheses on histological cross sections of bone.     

A comparison of primate, carnivoran and rodent limb bone cross-sectional properties: are primates really unique?

The cross-sectional properties of mammalian limb bones provide an important source of information about their loading history and locomotor adaptations. It has been suggested, for instance, that the cross-sectional strength of primate limb bones differs from that of other mammals as a consequence of living in a complex arboreal environment (Kimura, 1991, 1995). In order to test this hypothesis more rigorously, we have investigated cross-sectional properties in samples of humeri and femora of 71 primate species, 30 carnivorans and 59 rodents. Primates differ from carnivorans and rodents in having limb bones with greater cross-sectional strength than mammals of similar mass. This might imply that primates have stronger bones than carnivorans and rodents. However, primates also have longer proximal limb bones than other mammals. When cross-sectional dimensions are regressed against bone length, primates appear to have more gracile bones than other mammals. These two seemingly contradictory findings can be reconciled by recognizing that most limb bones experience bending as a predominant loading regime. After regressing cross-sectional strength against the product of body mass and bone length, a product which should be proportional to the bending moments applied to the limb, primates are found to overlap considerably with carnivorans and rodents. Consequently, primate humeri and femora are similar to those of nonprimates in their resistance to bending. Comparisons between arboreal and terrestrial species within the orders show that the bones of arboreal carnivorans have greater cross-sectional properties than those of terrestrial carnivorans, thus supporting Kimura's general notion. However, no differences were found between arboreal and terrestrial rodents. Among primates, the only significant difference was in humeral bending rigidity, which is higher in the terrestrial species. In summary, arboreal and terrestrial species do not show consistent differences in long bone reinforcement, and Kimura's conclusions must be modified to take into account the interaction of bone length and cross-sectional geometry.     

Cross-sectional structural parameters from densitometry

Bone densitometry has previously been used to obtain cross-sectional properties of bone from a single X-ray projection across the bone width. Using three unique projections, we have extended the method to obtain the principal area moments of inertia and orientations of the principal axes at each scan cross-section along the length of the scan. Various aluminum phantoms were used to examine scanner characteristics to develop the highest accuracy possible for in vitro non-invasive analysis of cross-sectional properties. Factors considered included X-ray photon energy, initial scan orientation, the angle spanned by the three scans (included angle), and I(min)/I(max) ratios. Principal moments of inertia were accurate to within +/-3.1% and principal angles were within +/-1 degrees of the expected value for phantoms scanned with included angles of 60 degrees and 90 degrees at the higher X-ray photon energy (140 kVp). Low standard deviations in the error (0.68-1.84%) also indicate high precision of calculated measurements with these included angles. Accuracy and precision decreased slightly when the included angle was reduced to 30 degrees. The method was then successfully applied to a pair of excised cadaveric tibiae. The accuracy and insensitivity of the algorithms to cross-sectional shape and changing isotropy (I(min)/I(max)) values when various included angles are used make this technique viable for future in vivo studies.     

Patterns of strain in the macaque ulna during functional activity

In vivo bone strain experiments were performed on the ulnae of three female rhesus macaques to test how the bone deforms during locomotion. The null hypothesis was that, in an animal moving its limbs predominantly in sagittal planes, the ulna experiences anteroposterior bending. Three rosette strain gauges were attached around the circumference of the bone slightly distal to midshaft. They permit a complete characterization of the ulna's loading environment. Strains were recorded during walking and galloping activities. Principal strains and strain directions relative to the long axis of the bone were calculated for each gauge site. In all three animals, the lateral cortex experienced higher tensile than compressive principal strains during the stance phase of walking. Compressive strains predominated at the medial cortex of two animals (the gauge on this cortex of the third animal did not function). The posterior cortex was subject to lower strains; the nature of the strain was highly dependent on precise gauge position. The greater principal strains were aligned closely with the long axis of the bone in two animals, whereas they deviated up to 45° from the long axis in the third animal. A gait change from walk to gallop was recorded for one animal. It was not accompanied by an incremental change in strain magnitudes. Strains are at the low end of the range of strain magnitudes recorded for walking gaits of nonprimate mammals. The measured distribution of strains in the rhesus monkey ulna indicates that mediolateral bending, rather than anteroposterior bending, is the predominant loading regime, with the neutral axis of bending running from anterior and slightly medial to posterior and slightly lateral. A variable degree of torsion was superimposed over this bending regime. Ulnar mediolateral bending is apparently caused by a ground reaction force vector that passes medial to the forearm. The macaque ulna is not reinforced in the plane of bending. The lack of buttressing in the loaded plane and the somewhat counterintuitive bending direction recommend caution with regard to conventional interpretations of long bone cross-sectional geometry. Am J Phys Anthropol 106:87–100, 1998. © 1998 Wiley-Liss, Inc.     

Determining effective centroid position in biomechanical testing: a technique for simplifying whole bone analysis

BACKGROUND: Whole bone in vitro biomechanical compressive testing can be complicated by three factors: sample asymmetry, heterogeneous material properties, and unknown effective centroid location. METHOD OF APPROACH: The technique presented here facilitates the calculation of effective centroid position, modulus of elasticity and equivalent uniform strain magnitude for a cross section of bone from a simple whole bone compressive test. Simplification of section response to load is achieved through a combination of linear beam and simple planer geometry theory. The technique requires three longitudinal strain gauges be affixed around the test specimen cross section of interest, gauge position need not be determined. Sample loading is then accomplished using a simple four point loading jig. RESULTS: Results of the technique are presented on an object with known elasticity and geometry, an aluminium tube, and seven pairs of equine third metacarpal whole bones. CONCLUSIONS: Average cross section modulus of elasticity, equivalent uniform cross section strain, and effective centroid locations were all predicted to within the range of published values. Employing the testing setup and analysis technique presented in this paper resulted in a significant savings in both implementation complexity and cost over previously available techniques.     

Coordinate systems for the carpal bones of the wrist

The eight small and complexly shaped carpal bones of the wrist articulate in six degrees of freedom with each other and to some extent with the radius and the metacarpals. With the increasing number and sophistication of studies of the carpus, a standardized definition for a coordinate system for each the carpal bones would aid in the reporting and comparison of findings. This paper presents a method for defining and constructing a coordinate system specific to each of the eight carpal bones based upon the inertial properties of the bone, derived from surface models constructed from three-dimensional (3-D) medical image volumes. Surface models from both wrists of 5 male and 5 female subjects were generated from CT image volumes in two neutral wrist positions (functional and clinical). An automated algorithm found the principal inertial axes and oriented them according to preset conditions in 85% of the bones, the remaining bones were corrected manually. Six of the eight carpal bones were significantly more extended in the functional neutral position than in the clinical neutral position. Gender had no significant effect on carpal bone posture in either wrist position. Correlations between the 3-D carpal posture and the commonly used 2-D clinical radiographic carpal angles are established. 3-D coordinate systems defined by the anatomy of the carpal bone, such as the ones presented here, are necessary to completely describe 3-D changes in the posture of the carpal bones.     

Computed-tomography-based finite-element models of long bones can accurately capture strain response to bending and torsion

Finite element (FE) models of long bones constructed from computed-tomography (CT) data are emerging as an invaluable tool in the field of bone biomechanics. However, the performance of such FE models is highly dependent on the accurate capture of geometry and appropriate assignment of material properties. In this study, a combined numerical-experimental study is performed comparing FE-predicted surface strains with strain-gauge measurements. Thirty-six major, cadaveric, long bones (humerus, radius, femur and tibia), which cover a wide range of bone sizes, were tested under three-point bending and torsion. The FE models were constructed from trans-axial volumetric CT scans, and the segmented bone images were corrected for partial-volume effects. The material properties (Young's modulus for cortex, density-modulus relationship for trabecular bone and Poisson's ratio) were calibrated by minimizing the error between experiments and simulations among all bones. The R(2) values of the measured strains versus load under three-point bending and torsion were 0.96-0.99 and 0.61-0.99, respectively, for all bones in our dataset. The errors of the calculated FE strains in comparison to those measured using strain gauges in the mechanical tests ranged from -6% to 7% under bending and from -37% to 19% under torsion. The observation of comparatively low errors and high correlations between the FE-predicted strains and the experimental strains, across the various types of bones and loading conditions (bending and torsion), validates our approach to bone segmentation and our choice of material properties.     

Locomotor variation and bending regimes of capuchin limb bones

Primates are very versatile in their modes of progression, yet laboratory studies typically capture only a small segment of this variation. In vivo bone strain studies in particular have been commonly constrained to linear locomotion on flat substrates, conveying the potentially biased impression of stereotypic long bone loading patterns. We here present substrate reaction forces (SRF) and limb postures for capuchin monkeys moving on a flat substrate (“terrestrial”), on an elevated pole (“arboreal”), and performing turns. The angle between the SRF vector and longitudinal axes of the forearm or leg is taken as a proxy for the bending moment experienced by these limb segments. In both frontal and sagittal planes, SRF vectors and distal limb segments are not aligned, but form discrepant angles; that is, forces act on lever arms and exert bending moments. The positions of the SRF vectors suggest bending around oblique axes of these limb segments. Overall, the leg is exposed to greater moments than the forearm. Simulated arboreal locomotion and turns introduce variation in the discrepancy angles, thus confirming that expanding the range of locomotor behaviors studied will reveal variation in long bone loading patterns that is likely characteristic of natural locomotor repertoires. “Arboreal” locomotion, even on a linear noncompliant branch, is characterized by greater variability of force directions and discrepancy angles than “terrestrial” locomotion (significant for the forearm only), partially confirming the notion that life in trees is associated with greater variation in long bone loading. Directional changes broaden the range of external bending moments even further. Am J Phys Anthropol, 2009. © 2009 Wiley‐Liss, Inc.     

Computer modelling of bone’s adaptation: the role of normal strain,shear strain and fluid flow

Bone loss is a serious health problem. In vivo studies have found that mechanical stimulation may inhibit bone loss as elevated strain in bone induces osteogenesis, i.e. new bone formation. However, the exact relationship between mechanical environment and osteogenesis is less clear. Normal strain is considered as a prime stimulus of osteogenic activity; however, there are some instances in the literature where osteogenesis is observed in the vicinity of minimal normal strain, specifically near the neutral axis of bending in long bones. It suggests that osteogenesis may also be induced by other or secondary components of mechanical environment such as shear strain or canalicular fluid flow. As it is evident from the literature, shear strain and fluid flow can be potent stimuli of osteogenesis. This study presents a computational model to investigate the roles of these stimuli in bone adaptation. The model assumes that bone formation rate is roughly proportional to the normal, shear and fluid shear strain energy density above their osteogenic thresholds. In vivo osteogenesis due to cyclic cantilever bending of a murine tibia has been simulated. The model predicts results close to experimental findings when normal strain, and shear strain or fluid shear were combined. This study also gives a new perspective on the relation between osteogenic potential of micro-level fluid shear and that of macro-level bending shear. Attempts to establish such relations among the components of mechanical environment and corresponding osteogenesis may ultimately aid in the development of effective approaches to mitigating bone loss.     

Allometry between length and cross-sectional dimensions of the femur and tibia in Homo sapiens sapiens

Allometric equations relating length and cross-sectional geometric properties of the femur and tibia are generated using skeletal remains from three recent human population samples. Approximate isometry, or geometric similarity, is found both within and between samples. Cross-sectional areas scale to approximately length2, while second moments of area scale to approximately length4. It is shown that this is consistent with the maintenance of equivalent mechanical stress in long bones of different length under dynamic loadings in vivo. Other evidence indicates that bending and torsional loadings are more critical than axial loadings in the determination of lower limb bone cross-sectional dimensions.     

Determination of mechanical stiffness of bone by pQCT measurements: correlation with non-destructive mechanical four-point bending test data

Mechanical tests of bone provide valuable information about material and structural properties important for understanding bone pathology in both clinical and research settings, but no previous studies have produced applicable non-invasive, quantitative estimates of bending stiffness. The goal of this study was to evaluate the effectiveness of using peripheral quantitative computed tomography (pQCT) data to accurately compute the bending stiffness of bone. Normal rabbit humeri (N=8) were scanned at their mid-diaphyses using pQCT. The average bone mineral densities and the cross-sectional moments of inertia were computed from the pQCT cross-sections. Bending stiffness was determined as a function of the elastic modulus of compact bone (based on the local bone mineral density), cross-sectional moment of inertia, and simulated quasistatic strain rate. The actual bending stiffness of the bones was determined using four-point bending tests. Comparison of the bending stiffness estimated from the pQCT data and the mechanical bending stiffness revealed excellent correlation (R2=0.96). The bending stiffness from the pQCT data was on average 103% of that obtained from the four-point bending tests. The results indicate that pQCT data can be used to accurately determine the bending stiffness of normal bone. Possible applications include temporal quantification of fracture healing and risk management of osteoporosis or other bone pathologies.     

How do trabeculae affect the calculation of structural properties of bone?

One difficulty that arises in an analysis of the cross-sectional properties of bone is whether to include cancellous bone in the analysis. The purpose of this paper is to determine how different amounts of cancellous bone affect the measurement of structural properties of bone cross-sections. Thirty-two tibial and femoral cross-sections were chosen at random from a series of cross-sectioned nonhuman primate bones. Geometrical properties were calculated for the cross-sections, and torsional and bending stress analyses were performed. The results suggest that the effect of including cancellous bone in the analysis is closely related to the amount of bone, where it lies within the cross-section, and the type of analysis performed. Including cancellous bone in calculations of structural properties of bone cross-sections may cause the strength and stiffness of the bone to be exaggerated.     

In vivo strains in the femur of the nine‐banded armadillo (Dasypus novemcinctus)

The capacity of limb bones to resist the locomotor loads they encounter depends on both the pattern of those loads and the material properties of the skeletal elements. Among mammals, understanding of the interplay between these two factors has been based primarily on evidence from locomotor behaviors in upright placentals, which show limb bones that are loaded predominantly in anteroposterior bending with minimal amounts of torsion. However, loading patterns from the femora of opossums, marsupials using crouched limb posture, show appreciable torsion while the bone experiences mediolateral (ML) bending. These data indicated greater loading diversity in mammals than was previously recognized, and suggested the possibility that ancestral loading patterns found in sprawling lineages (e.g., reptilian sauropsids) might have been retained among basal mammals. To further test this hypothesis, we recorded in vivo locomotor strains from the femur of the nine‐banded armadillo (Dasypus novemcinctus), a member of the basal xenarthran clade of placental mammals that also uses crouched limb posture. Orientations of principal strains and magnitudes of shear strains indicate that armadillo femora are exposed to only limited torsion; however, bending is essentially ML, placing the medial aspect of the femur in compression and the lateral aspect in tension. This orientation of bending is similar to that found in opossums, but planar strain analyses indicate much more of the armadillo femur experiences tension during bending, potentially due to muscles pulling on the large, laterally positioned third trochanter. Limb bone safety factors were estimated between 3.3 and 4.3 in bending, similar to other placental mammals, but lower than opossums and most sprawling taxa. Thus, femoral loading patterns in armadillos show a mixture of similarities to both opossums (ML bending) and other placentals (limited torsion and low safety factors), along with unique features (high axial tension) that likely relate to their distinctive hindlimb anatomy. J. Morphol. 26:889–899, 2015. © 2015 Wiley Periodicals, Inc.     

Comparison of the hindfoot axes of a multi-segment foot model to the underlying bony anatomy

Musculoskeletal models used in gait analysis require coordinate systems to be identified for the body segments of interest. It is not obvious how hindfoot (or rearfoot) axes defined by skin-mounted markers relate to the anatomy of the underlying bones. The aim of this study was to compare the marker-based axes of the hindfoot in a multi-segment foot model to the orientations of the talus and calcaneus as characterized by their principal axes of inertia. Twenty adult females with no known foot deformities had radio-opaque markers placed on their feet and ankles at the foot model marker locations. CT images of the feet were acquired as the participants lay supine with their feet in a semi-weight bearing posture. The spatial coordinates of the markers were obtained from the images and used to define the foot model axes. Segmented masks of the tali and calcanei were used to create 3D bone models, from which the principal axes of the bones were obtained. The orientations of the principal axes were either within the range of typical values reported in the imaging literature or differed in ways that could be explained by variations in how the angles were defined. The model hindfoot axis orientations relative to the principal axes of the bones had little bias but were highly variable. Consideration of coronal plane hindfoot alignment as measured clinically and radiographically suggested that the model hindfoot coordinate system represents the posterior calcaneal tuberosity, rather than the calcaneus as a whole.     

Relationships of loading history and structural and material characteristics of bone: development of the mule deer calcaneus

If a bone's morphologic organization exhibits the accumulated effects of its strain history, then the relative contributions of a given strain stimulus to a bone's development may be inferred from a bone's hierarchical organization. The artiodactyl calcaneus is a short cantilever, loaded habitually in bending, with prevalent compression in the cranial (Cr) cortex, tension in the caudal (Cd) cortex, and shear in the medial and lateral cortices (i.e., neutral axis). Artiodactyl calcanei demonstrate unusually heterogeneous structural and material organization between these cortices. This study examines potential relationships between developmental morphologic variations and the functional strain distribution of the deer calcaneus. One calcaneus was obtained from each of 36 (fetus to adult) wild deer. Predominant collagen fiber orientation (CFO), microstructural characteristics, mineral content (% ash), and geometric parameters were determined from transversely cut segments. Radiographs were examined for arched trabeculae, which may reflect tension/compression stress trajectories. Results showed that cross-sectional shape changes with age from quasi-circular to quasi-elliptical, with the long axis in the cranial-caudal direction of habitual bending. Cranial ("compression") cortical thickness increased at a greater rate than the Cd ("tension") cortex. Fetal bones exhibited arched trabeculae. Percent ash was not uniform (Cr > Cd), and this disparity increased with age (absolute differences: 2.5% fetuses, 4.3% adults). Subadult bones showed progressively more secondary osteons and osteocyte lacunae in the Cr cortex, but the Cd cortex tended to have more active remodeling in the subadult and adult bones. Nonuniform Cr:Cd CFO patterns first consistently appear in the subadults, and are correlated with secondary bone formation and habitual strain mode. Medial and lateral cortices in these groups exhibited elongated secondary osteons. These variations may represent "strain-mode-specific" (i.e., tension, compression, shear) adaptations. The heterogeneous organization may also be influenced by variations in longitudinal strain magnitude (highest in the Cr cortex) and principal strain direction-oblique in medial-lateral cortices (where shear strains also predominate). Other factors such as local reductions in longitudinal strain may influence the increased remodeling activity of the Cd cortex. Some structural variations, such as arched trabeculae, that are established early in ontogeny may be strongly influenced by genetic- or epigenetic-derived processes. Material variations, such as secondary osteon population densities and CFO, which appear later, may be products of extragenetic factors, including microdamage.     

Cortical Structure of Hallucal Metatarsals and Locomotor Adaptations in Hominoids

Diaphyseal morphology of long bones, in part, reflects in vivo loads experienced during the lifetime of an individual. The first metatarsal, as a cornerstone structure of the foot, presumably expresses diaphyseal morphology that reflects loading history of the foot during stance phase of gait. Human feet differ substantially from those of other apes in terms of loading histories when comparing the path of the center of pressure during stance phase, which reflects different weight transfer mechanisms. Here we use a novel approach for quantifying continuous thickness and cross-sectional geometric properties of long bones in order to test explicit hypotheses about loading histories and diaphyseal structure of adult chimpanzee, gorilla, and human first metatarsals. For each hallucal metatarsal, 17 cross sections were extracted at regularly-spaced intervals (2.5% length) between 25% and 65% length. Cortical thickness in cross sections was measured in one degree radially-arranged increments, while second moments of area were measured about neutral axes also in one degree radially-arranged increments. Standardized thicknesses and second moments of area were visualized using false color maps, while penalized discriminant analyses were used to evaluate quantitative species differences. Humans systematically exhibit the thinnest diaphyseal cortices, yet the greatest diaphyseal rigidities, particularly in dorsoplantar regions. Shifts in orientation of maximum second moments of area along the diaphysis also distinguish human hallucal metatarsals from those of chimpanzees and gorillas. Diaphyseal structure reflects different loading regimes, often in predictable ways, with human versus non-human differences probably resulting both from the use of arboreal substrates by non-human apes and by differing spatial relationships between hallux position and orientation of the substrate reaction resultant during stance. The novel morphological approach employed in this study offers the potential for transformative insights into form-function relationships in additional long bones, including those of extinct organisms (e.g., fossils).     

Characterizing gait induced normal strains in a murine tibia cortical bone defect model

The critical role that mechanical stimuli serve in mediating bone repair is recognized but incompletely understood. Further, previous attempts to understand this role have utilized application of externally applied mechanical loads to study the tissue’s response. In this project, we have therefore endeavored to capitalize on bone’s own consistently diverse loading environment to develop a novel model that would enable assessment of the influence of physiologically engendered mechanical stimuli on cortical defect repair. We used an inverse dynamics approach with finite element analysis (FEA) to first quantify normal strain distributions generated in mouse tibia during locomotion. The strain environment of the tibia, as previously reported for other long bones, was found to arise primarily due to bending and was consistent in orientation through the stance phase of gait. Based on these data, we identified three regions within a transverse cross-section of the mid-diaphysis as uniform locations of either peak tension, peak compression, or the neutral axis of bending (i.e. minimal strain magnitude). We then used FEA to quantify the altered strain environment that would be produced by a 0.6 mm diameter cylindrical cortical bone defect at each diaphyseal site and, in an in situ study confirmed our ability to accurately place defects at the desired diaphyseal locations. The resulting model will enable the exploration of cortical bone healing within the context of physiologically engendered mechanical strain.