Physiologically based boundary conditions in finite element modelling

Finite element analysis has been used extensively in the study of bone loading and implant performance, such as in the femur. The boundary conditions applied vary widely, generally producing excessive femoral deformation, and although it has been shown that the muscle forces influence femoral deflections and loading, little consideration has been given to the displacement constraints. It is hypothesised that careful application of physiologically based constraints can produce physiological deformation, and therefore straining, of the femur. Joint contact forces and a complete set of muscle forces were calculated based on the geometry of the Standardised Femur using previously validated musculoskeletal models. Five boundary condition cases were applied to a finite element model of the Standardised Femur: (A) diaphyseally constrained with hip contact and abductor forces; (B) case A plus vasti forces; (C) case A with complete set of muscle forces; (D) distally constrained with all muscle forces; (E) physiological constraints with all muscle forces. It was seen that only the physiological boundary conditions, case E, produced physiological deflections (< 2.0mm) of the femoral head in both the coronal and sagittal planes, which resulted in minimal reaction forces at the constrained nodes. Strains in the mid-diaphysis varied by up to 600 micro-strain under walking loads and 1000 micro-strain under stair climbing loads. The mode of loading, as indicated by the strain profiles on the cortex also varied substantially under these boundary conditions, which has important consequences for studies that examine localised bone loading such as fracture or bone remodelling simulations.     

Finite-element-analysis and experimental investigation of stresses in a femur

The state of stress of an intact femur was analysed using a three-dimensional finite element model. One of a pair of femora was used for determination of data for the 3-D model. The other was instrumented with 34 rosette strain gauges for experimental measurements. Good agreement of analytical and experimental data was achieved. For six modes of loading, the deflections and the principal and comparison stresses were determined and compared. The upper one third and the diaphysis of the femur are differently affected in their state of stress, by different modes of loading including simulation of the abductor muscles and the iliotibial tract. For qualitative stress studies of the diaphysis, loading by a single force on the femoral head parallel to the shaft axis may be adequate.     

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.     

Prediction of strength and strain of the proximal femur by a CT-based finite element method

Hip fractures are the most serious complication of osteoporosis and have been recognized as a major public health problem. In elderly persons, hip fractures occur as a result of increased fragility of the proximal femur due to osteoporosis. It is essential to precisely quantify the strength of the proximal femur in order to estimate the fracture risk and plan preventive interventions. CT-based finite element analysis could possibly achieve precise assessment of the strength of the proximal femur. The purpose of this study was to create a simulation model that could accurately predict the strength and surface strains of the proximal femur using a CT-based finite element method and to verify the accuracy of our model by load testing using fresh frozen cadaver specimens. Eleven right femora were collected. The axial CT scans of the proximal femora were obtained with a calibration phantom, from which the 3D finite element models were constructed. Materially nonlinear finite element analyses were performed. The yield and fracture loads were calculated, while the sites where elements failed and the distributions of the principal strains were determined. The strain gauges were attached to the proximal femoral surfaces. A quasi-static compression test of each femur was conducted. The yield loads, fracture loads and principal strains of the prediction significantly correlated with those measured (r=0.941, 0.979, 0.963). Finite element analysis showed that the solid elements and shell elements in undergoing compressive failure were at the same subcapital region as the experimental fracture site.     

Fracture prediction for the proximal femur using finite element models: Part I--Linear analysis.

Over 90 percent of the more than 250,000 hip fractures that occur annually in the United States are the result of falls from standing height. Despite this, the stresses associated with femoral fracture from a fall have not been investigated previously. Our objectives were to use three-dimensional finite element models of the proximal femur (with geometries and material properties based directly on quantitative computed tomography) to compare predicted stress distributions for one-legged stance and for a fall to the lateral greater trochanter. We also wished to test the correspondence between model predictions and in vitro strain gage data and failure loads for cadaveric femora subjected to these loading conditions. An additional goal was to use the model predictions to compare the sensitivity of several imaging sites in the proximal femur which are used for the in vivo prediction of hip fracture risk. In this first of two parts, linear finite element models of two unpaired human cadaveric femora were generated. In Part II, the models were extended to include nonlinear material properties for the cortical and trabecular bone. While there was poor correspondence between strain gage data and model predictions, there was excellent agreement between the in vitro failure data and the linear model, especially using a von Mises effective strain failure criterion. Both the onset of structural yielding (within 22 and 4 percent) and the load at fracture (within 8 and 5 percent) were predicted accurately for the two femora tested. For the simulation of one-legged stance, the peak stresses occurred in the primary compressive trabeculae of the subcapital region.(ABSTRACT TRUNCATED AT 250 WORDS)     

Simulation of physiological loading in total hip replacements

The determination of biomechanical force systems of implanted femurs to obtain adequate strain measurements has been neglected in many published studies. Due to geometric alterations induced by surgery and those inherent to the design of the prosthesis, the loading system changes because the lever arms are modified. This paper discusses the determination of adequate loading of the implanted femur based on the intact femur-loading configuration. Four reconstructions with Lubinus SPII, Charnley Roundback, Muller Straight and Stanmore prostheses were used in the study. Pseudophysiologic and nonphysiologic implanted system forces were generated and assessed with finite element analysis. Using an equilibrium system of forces composed by the Fx (medially direction) component of the hip contact force and the bending moments Mx (median plane) and My (coronal plane) allowed adequate, pseudo-physiological loading of the implanted femur. We suggest that at least the bending moment at the coronal plane must be restored in the implanted femur-loading configuration.     

Bone Load Estimation for the Proximal Femur Using Single Energy Quantitative CT Data

A density-based load estimation method was applied to determine femoral load patterns. Two-dimensional finite element models were constructed using single energy quantitative computed tomography (QCT) data from two femora. Basic load cases included parabolic pressure joint loads and constant tractions on the greater trochanter. An optimization procedure adjusted magnitudes of the basic load cases, such that the applied mechanical stimulus approached the ideal stimulus throughout each model. Dominant estimated load directions were generally consistent with published experimental data for gait. Other estimated loads suggested that loads at extreme joint orientations may be important to maintenance of bone structure. Remodeling simulations with the estimated loads produced density distributions qualitatively similar to the QCT data sets. Average nodal density errors between QCT data and predictions were 0.24 g/cm(3) and 0.28 g/cm(3). The results indicate that density-based load estimation could improve understanding of loading patterns on bones.     

Changes in the Architecture of Cancellous Bone in the Femora of Developing Mice as a Result of Short Duration Hypergravity

SYNOPSIS. Mice were used as a model system to demonstrate theeffects of increased loading on the trabecular bone in the developingproximal femur. A two-dimensional free body diagram revealedthat because the femur is normally held in a horizontal positionduring support, functional loading resulted in a net posteriorbending moment and a muscle generated axial compressive loadthat was 7.7 times the normal support weight for that limb.Experimental treatment consisted of four times normal gravityfor 10 min, six times per day for 30 days in addition to unrestrainednormal exercise. Unrestrained normal exercise served as thecontrol. An animal whose foot was amputated by his mother shortlyafter birth was used as a less than normal load example forcomparison with the hypergravity and control groups. All groupswere treated for 30 days. Two areas of trabecular bone wereexamined: a proximal central webbing area that may be underpredominantly tensile loading much of the time, and a more distalarea near the third trochanter that presumably experiences bendingstresses when the animal is exposed to normal activity or hypergravity.Results indicate that: 1) The central webbing area appears tohave a similar mass of bone regardless of loading conditions;however, with increased loading the trabeculae become orientedin directions that approximate directions of presumed principaltensile stresses. 2) There is no genetically predetermined amountof bony tissue in the third trochanter area; when the limb isrelatively unloaded no trabeculae are present, as the loadingincreases trabecular bone is formed. In the hypergravity group,trabeculae in the area of the third trochanter increased inlength and complexity and became more centrally positioned inthe lumen of the femur, a position that presumably would allowthe femur to better resist bending moments. These results demonstratethat short duration, high intensity functional loading can stimulategrowth of trabecular bone, and further, specific types of functionalloading can be correlated with specific architectural changes.     

Initiation and progression of mechanical damage in the intervertebral disc under cyclic loading using continuum damage mechanics methodology: A finite element study

It is difficult to study the breakdown of disc tissue over several years of exposure to bending and lifting by experimental methods. There is also no finite element model that elucidates the failure mechanism due to repetitive loading of the lumbar motion segment. The aim of this study was to refine an already validated poro-elastic finite element model of lumbar motion segment to investigate the initiation and progression of mechanical damage in the disc under simple and complex cyclic loading conditions. Continuum damage mechanics methodology was incorporated into the finite element model to track the damage accumulation in the annulus in response to the repetitive loading. The analyses showed that the damage initiated at the posterior inner annulus adjacent to the endplates and propagated outwards towards its periphery under all loading conditions simulated. The damage accumulated preferentially in the posterior region of the annulus. The analyses also showed that the disc failure is unlikely to happen with repetitive bending in the absence of compressive load. Compressive cyclic loading with low peak load magnitude also did not create the failure of the disc. The finite element model results were consistent with the experimental and clinical observations in terms of the region of failure, magnitude of applied loads and the number of load cycles survived.     

A finite element model of the human knee joint for the study of tibio-femoral contact

As a step towards developing a finite element model of the knee that can be used to study how the variables associated with a meniscal replacement affect tibio-femoral contact, the goals of this study were 1) to develop a geometrically accurate three-dimensional solid model of the knee joint with special attention given to the menisci and articular cartilage, 2) to determine to what extent bony deformations affect contact behavior, and 3) to determine whether constraining rotations other than flexion/extension affects the contact behavior of the joint during compressive loading. The model included both the cortical and trabecular bone of the femur and tibia, articular cartilage of the femoral condyles and tibial plateau, both the medial and lateral menisci with their horn attachments, the transverse ligament, the anterior cruciate ligament, and the medial collateral ligament. The solid models for the menisci and articular cartilage were created from surface scans provided by a noncontacting, laser-based, three-dimensional coordinate digitizing system with an root mean squared error (RMSE) of less than 8 microns. Solid models of both the tibia and femur were created from CT images, except for the most proximal surface of the tibia and most distal surface of the femur which were created with the three-dimensional coordinate digitizing system. The constitutive relation of the menisci treated the tissue as transversely isotropic and linearly elastic. Under the application of an 800 N compressive load at 0 degrees of flexion, six contact variables in each compartment (ie., medial and lateral) were computed including maximum pressure, mean pressure, contact area, total contact force, and coordinates of the center of pressure. Convergence of the finite element solution was studied using three mesh sizes ranging from an average element size of 5 mm by 5 mm to 1 mm by 1 mm. The solution was considered converged for an average element size of 2 mm by 2 mm. Using this mesh size, finite element solutions for rigid versus deformable bones indicated that none of the contact variables changed by more than 2% when the femur and tibia were treated as rigid. However, differences in contact variables as large as 19% occurred when rotations other than flexion/extension were constrained. The largest difference was in the maximum pressure. Among the principal conclusions of the study are that accurate finite element solutions of tibio-femoral contact behavior can be obtained by treating the bones as rigid. However, unrealistic constraints on rotations other than flexion/extension can result in relatively large errors in contact variables.     

Biomechanics of the hip and birth in early Homo

A complex of traits in the femur and pelvis of Homo ereclus and early “erectus-like” specimens has been described, but never satisfactorily explained. Here the functional relationships between pelvic and femoral structure in humans are explored using both theoretical biomechanical models and empirical tests within modern samples of diverse body form (Pecos Amerindians, East Africans). Results indicate that a long femoral neck increases mediolateral bending of the femoral diaphysis and decreases gluteal abductor and hip joint reaction forces. Increasing biacetabular breadth along with femoral neck length further increases M-L bending of the femoral shaft and maintains abductor and joint reaction forces at near “normal” levels. When compared to modern humans, Homo erectus and early “erectus-like” specimens are characterized by a long femoral neck and greatly increased M-L relative to A-P bending strength of the femoral shaft, coupled with no decrease in hip joint size and a probable increase in abductor force relative to body size. All of this strongly suggests that biacetabular breadth as well as femoral neck length was relatively large in early Homo. Several features preserved in early Homo partial hip bones also indicate that the true (lower) pelvis was very M-L broad, as well as A-P narrow. This is similar to the lower pelvic shape of australopithecines and suggests that nonrotational birth, in which the newborn's head is oriented transversely through the pelvic outlet, characterized early Homo as well as Australopithecus. Because M-L breadth of the pelvis is constrained by other factors, this may have limited increases in cranial capacity within Homo until rotational birth was established during the late Middle Pleistocene. During or after the transition to rotational birth biacetabular breadth decreased, reducing the body weight moment arm about the hip and allowing femoral neck length (abductor moment arm) to also decrease, both of which reduced M-L bending of the proximal femoral shaft. Variation in femoral structural properties within early Homo and other East African Early Pleistocene specimens has several taxonomic and phylogenetic implications. © 1995 Wiley-Liss, Inc.     

Effect of boundary conditions,impact loading and hydraulic stiffening on femoral fracture strength

Patient specific quantitative CT (QCT) imaging data together with the finite element (FE) method may provide an accurate prediction of a patient's femoral strength and fracture risk. Although numerous FE models investigating femoral fracture strength have been published, there is little consent on the effect of boundary conditions, dynamic loading and hydraulic strengthening due to intra-medullary pressure on the predicted fracture strength. We developed a QCT-derived FE model of a proximal femur that included node-specific modulus assigned based on the local bone density. The effect of three commonly used boundary conditions published in literature were investigated by comparing the resulting strain field due to an applied fracture load. The models were also augmented with viscoelastic material properties and subject to a realistic impact load profile to determine the effect of dynamic loads on the strain field. Finally, the effect of hydraulic strengthening was investigated by including node specific permeability and performing a coupled pore diffusion and stress analysis of the FE model. Results showed that all boundary conditions yield the same strain field patterns, but peak strains were 22% lower and fracture load was 18% higher when loaded at the greater trochanter than when loaded at the femoral head. Comparison of the dynamic models showed that material viscoelasticity was important, but inertial effects (vibration and shock) were not. Finally, pore pressure changes did not cause significant hydraulic strengthening of bone under fall impact loading.     

Mechanobiological prediction of proximal femoral deformities in children with cerebral palsy

Children with cerebral palsy (CP) walk with altered gait and frequently exhibit proximal femoral deformities, such as anteversion and coxa valga. The objective of this research was to investigate the effect of specific gait patterns on the femoral morphology in CP.

In this study, the mechanobiological principles were implemented on a 3D finite element (FE) model of the proximal femur in order to predict changes in morphology over time in healthy and CP children. This model relies on the assumption that cyclic octahedral shear stress promotes growth and cyclic hydrostatic compressive stress inhibits growth. Growth was simulated over 16 iterations, representing approximately 5 months of growth.

The FE model predicts an increase in the femoral anteversion and coxa valga for CP loading conditions when compared with healthy ones. Understanding the role of loading in skeletal morphogenesis may help prevent bone deformities and improve function in children with gait abnormalities.     

Bone Load Estimation for the Proximal Femur Using Single Energy Quantitative CT Data

A density-based load estimation method was applied to determine femoral load patterns. Two-dimensional finite element models were constructed using single energy quantitative computed tomography (QCT) data from two femora. basic load cases included parabolic pressure joint loads and constant tractions on the greater trochanter. An optimization procedure adjusted magnitudes of the basic load cases, such that the applied mechanical stimulus approached the ideal stimulus throughout each model. Dominant estimated load directions were generally consistent with published experimental data for gait. Other estimated loads suggested that loads at extreme joint orientations may be important to maintenance of bone structure. Remodeling simulations with the estimated loads produced density distributions qualitatively similar to the QCT data sets. Average nodal density errors between QCT data and predictions were 0·24 g/cm³ and 0·28 g/cm³. The results indicate that density-based load estimation could improve understanding of loading patterns on bones.     

Automated three-dimensional finite element modelling of bone: a new method

Three-dimensional finite element stress analysis of bone is a key to understanding bone remodelling, assessing fracture risk, and designing prostheses; however, the cost and complexity of predicting the stress field in bone with accuracy has precluded the routine use of this method. A new, automated method of generating patient-specific three-dimensional finite element models of bone is presented — it uses digital computed tomographic (CT) scan data to derive the geometry of the bone and to estimate its inhomogeneous material properties. Cubic elements of a user-specified size are automatically defined and then individually assigned the CT scan-derived material properties. The method is demonstrated by predicting the stress, strain, and strain energy in a human proximal femur in vivo. Three-dimensional loading conditions corresponding to the stance phase of gait were taken from the literature and applied to the model. Maximum principal compressive stresses of 8–23 MPa were computed for the medial femoral neck. Automated generation of additional finite element models with larger numbers of elements was used to verify convergence in strain energy.     

LISS钢板治疗股骨远端骨折有限元建模及分析

目的建立LISS-DF治疗股骨远端骨折近端螺钉不同单双皮质固定的三维有限元模型,并进行初步生物力学分析。方法提取CT图片相关数据,利用自行编写程序生成命令流文件,建立完整股骨以及16个不同LISS-DF治疗股骨远端AO分型33-A3型骨折的实体模型（钢板和股骨不接触、螺钉分别固定于钢板和股骨）,进行网格划分。分析不同载荷作用下完整股骨和LISS钢板近端螺钉全双皮质固定治疗骨折的模型受力状况。结果建立了相关的有限元模型。不同载荷作用下,LISS钢板近端螺钉全双皮质固定模型和完整股骨的应力集中均位于股骨颈内侧和股骨干外侧中下1/3处。相同载荷作用下,LISS钢板近端螺钉全双皮质固定模型的股骨颈部最大等效应力值略减小,股骨干最大等效应力值明显减小。结论研究建立的三维有限元模型,为应用LISS治疗股骨骨折的生物力学分析提供了良好的实验平台和基础。从生物力学角度而言,LISS-DF近端螺钉全双皮质固定为治疗股骨远端复杂骨折的有效方法。     

Proximal Femoral Density Patterns are Consistent with Bicentric Joint Loads

We developed an alternate method for density-based load estimation and applied it to estimate hip joint load distributions for two femora. Two-dimensional finite element models were constructed from single energy quantitative computed tomography (QCT) data. Load estimation was performed using five loading regions on the femoral head. Within each loading region, individual nodal loads, normal to the local surface, were supplied as input to the load estimation. An optimization procedure independently adjusted individual nodal load magnitudes in each region, and the magnitudes of muscle forces on the greater trochanter, such that the applied tissue stimulus approached the reference stimulus throughout the model. Dominant estimated load resultant directions were generally consistent with published experimental data for loads during gait. The estimated loads also suggested that loads near the extremes of the articulating surface may be important (even required) for development and maintenance of normal bone architecture. Estimated load distributions within nearly all regions predicted bicentric loading patterns, which are consistent with observations of hip joint incongruity. Remodeling simulations with the estimated loads predicted density distributions with features qualitatively similar to the QCT data sets. This study illustrates how applications of density-based bone load estimation can improve understanding of dominant loading patterns in other bones and joints. The prediction of bicentric loading suggests a very fine level of local adaptation to details of joint loading.     

Proximal Femoral Density Patterns are Consistent with Bicentric Joint Loads

We developed an alternate method for density-based load estimation and applied it to estimate hip joint load distributions for two femora. Two-dimensional finite element models were constructed from single energy quantitative computed tomography (QCT) data. Load estimation was performed using five loading regions on the femoral head. Within each loading region, individual nodal loads, normal to the local surface, were supplied as input to the load estimation. An optimization procedure independently adjusted individual nodal load magnitudes in each region, and the magnitudes of muscle forces on the greater trochanter, such that the applied tissue stimulus approached the reference stimulus throughout the model. Dominant estimated load resultant directions were generally consistent with published experimental data for loads during gait. The estimated loads also suggested that loads near the extremes of the articulating surface may be important (even required) for development and maintenance of normal bone architecture. Estimated load distributions within nearly all regions predicted bicentric loading patterns, which are consistent with observations of hip joint incongruity. Remodeling simulations with the estimated loads predicted density distributions with features qualitatively similar to the QCT data sets. This study illustrates how applications of density-based bone load estimation can improve understanding of dominant loading patterns in other bones and joints. The prediction of bicentric loading suggests a very fine level of local adaptation to details of joint loading.     

A novel methodology for generating 3D finite element models of the hip from 2D radiographs

Finite element (FE) modelling has been proposed as a tool for estimating fracture risk and patient-specific FE models are commonly based on computed tomography (CT). Here, we present a novel method to automatically create personalised 3D models from standard 2D hip radiographs. A set of geometrical parameters of the femur were determined from seven a–p hip radiographs and compared to the 3D femoral shape obtained from CT as training material; the error in reconstructing the 3D model from the 2D radiographs was assessed. Using the geometry parameters as the input, the 3D shape of another 21 femora was built and meshed, separating a cortical and trabecular compartment. The material properties were derived from the homogeneity index assessed by texture analysis of the radiographs, with focus on the principal tensile and compressive trabecular systems. The ability of these FE models to predict failure load as determined by experimental biomechanical testing was evaluated and compared to the predictive ability of DXA. The average reconstruction error of the 3D models was 1.77 mm (±1.17 mm), with the error being smallest in the femoral head and neck, and greatest in the trochanter. The correlation of the FE predicted failure load with the experimental failure load was r²=64% for the reconstruction FE model, which was significantly better (p<0.05) than that for DXA (r²=24%). This novel method for automatically constructing a patient-specific 3D finite element model from standard 2D radiographs shows encouraging results in estimating patient-specific failure loads.     

Mechanical significance of femoral head trabecular bone structure in Loris and Galago evaluated using micromechanical finite element models

Work on the interspecific and intraspecific variation of trabecular bone in the proximal femur of primates demonstrates important architectural variation between animals with different locomotor behaviors. This variation is thought to be related to the processes of bone adaptation whereby bone structure is optimized to the mechanical environment. Micromechanical finite element models were created for the proximal femur of the leaping Galago senegalensis and the climbing and quadrupedal Loris tardigradus by converting bone voxels from high-resolution X-ray computed tomography scans of the femoral head to eight-noded brick elements. The resulting models had approximately 1.8 million elements each. Loading conditions representing takeoff phase of a leap and more generalized load orientations were applied to the models, and the models were solved using the iterative "row-by-row" matrix-vector multiplication algorithm. The principal strain and Von Mises stress results for the leaping model were similar for both species at each load orientation. Similar hip joint reaction forces in the range of 4.9 x to 12 x body weight were calculated for both species under each loading condition, but the hip reaction values estimated for Loris were higher than predicted based on locomotor behavior. These results suggest that functional adaptation to hip joint loading may not fully explain the differences in femoral head trabecular bone structure in Galago and Loris. The finite element method represents a unique and useful tool for analyzing the functional adaptation of trabecular bone in a diversity of animals and for reconstructing locomotor behavior in extinct taxa.