Development of a three-dimensional finite element model of a human tibia using experimental modal analysis.

The modal analysis of a human tibia consisted of characterizing its dynamic behavior by determining natural frequency, damping ratio and mode shapes. Two methods were used to perform the modal analysis: (1) a finite element method (structural model); (2) an experimental modal analysis (modal model). The experimental modal model was used to optimize the structural model. After optimization, differences in results between the two models were found to be due only to mechanical properties and mass distribution. The influences of boundary conditions and geometric properties (such as inertia and length) were eliminated by the finite element model itself. The percent relative error between the two methods was approximately 3%, corresponding to the standard deviation of the measured frequencies. For the frequency range considered, the mode shapes were bending modes in two different vibration planes (latero-medial and sagittal), with a slight torsion effect due to the twisted geometry of the tibia.     

An accurate finite element model of the cervical spine under quasi-static loading

Cervical disc injury due to impact has been observed in clinical and biomechanical investigations; however, there is a lack of data that helps to elucidate the mechanisms of disc injury during these collisions. Therefore, it is necessary to understand the behavior of the cervical spine under different types of loading situations. A three dimensional finite element (FE) model for the multi-level cervical spine segment (C0-C7) was developed using computed tomography (CT) data and applied to study the internal stresses and strains of the intervertebral discs under quasi-static loading conditions. The intervertebral discs were treated as nonlinear, anisotropic and incompressible subjected to large deformations. The model accuracy was validated by comparing it with previously published experimental and numerical results for different movements. It was shown that the use of a fiber reinforced model to describe the behavior of the annulus of the discs would predict higher maximum shear strains than an isotropic one, being therefore important the use of complex constitutive models in order to be able to detect the appearance of injured zones, since those strains and stresses are supposed to be related with damage to soft tissues. Several movements were analyzed: flexion, extension and axial rotation, obtaining that the maximum shear stresses in the disc were higher for a flexo-extension movement.     

Biomechanical study of pediatric human cervical spine: a finite element approach

Although considerable effort has been made to understand the biomechanical behavior of the adult cervical spine, relatively little information is available on the response of the pediatric cervical spine to external forces. Since significant anatomical differences exist between the adult and pediatric cervical spines, distinct biomechanical responses are expected. The present study quantified the biomechanical responses of human pediatric spines by incorporating their unique developmental anatomical features. One-, three-, and six-year-old cervical spines were simulated using the finite element modeling technique, and their responses computed and compared with the adult spine response. The effects of pure overall structural scaling of the adult spine, local component developmental anatomy variations that occur to the actual pediatric spines, and structural scaling combined with local component anatomy variations on the responses of the pediatric spines were studied. Age- and component-related developmental anatomical features included variations in the ossification centers, cartilages, growth plates, vertebral centrum, facet joints, and annular fibers and nucleus pulposus of the intervertebral discs. The flexibility responses of the models were determined under pure compression, pure flexion, pure extension, and varying degrees of combined compression-flexion and compression-extension. The pediatric spine responses obtained with the pure overall (only geometric) scaling of the adult spine indicated that the flexibilities consistently increase in a uniform manner from six- to one-year-old spines under all loading cases. In contrast, incorporation of local anatomic changes specific to the pediatric spines of the three age groups (maintaining the same adult size) not only resulted in considerable increases in flexibilities, but the responses also varied as a function of the age of the pediatric spine and type of external loading. When the geometric scaling effects were added to these spines, the increases in flexibilities were slightly higher; however, the pattern of the responses remained the same as found in the previous approach. These results indicate that inclusion of developmental anatomical changes characteristic of the pediatric spines has more of a predominant effect on biomechanical responses than extrapolating responses of the adult spine based on pure overall geometric scaling.     

Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis

Magnetic resonance (MR) imaging has been widely used to evaluate the thickness and volume of articular cartilage both in vivo and in vitro. While morphological information on the cartilage can be obtained using MR images, image processing for extracting geometric boundaries of the cartilage may introduce variations in the thickness of the cartilage. To evaluate the variability of using MR images to construct finite element (FE) knee cartilage models, five investigators independently digitized the same set of MR images of a human knee. The topology of cartilage thickness was determined using a minimal distance algorithm. Less than 8 percent variation in cartilage thickness was observed from the digitized data. The effect of changes in cartilage thickness on contact stress analysis was then investigated using five FE models of the knee. One FE model (average FE model) was constructed using the mean values of the digitized contours of the cartilage, and the other four were constructed by varying the thickness of the average FE model by +/- 5 percent and +/- 10 percent, respectively. The results demonstrated that under axial tibial compressive loading (up to 1,400 N), variations of cartilage thickness caused by digitization of MR images may result in a difference of approximately 10 percent in peak contact stresses (surface pressure, von Mises stress, and hydrostatic pressure) in the cartilage. A reduction of cartilage thickness caused increases of contact stresses, while an increase of cartilage thickness reduced contact stresses. Furthermore, the effect of variation of material properties of the cartilage on contact stress analysis was investigated. The peak contact stress increased almost linearly with the Young's modulus of the cartilage. The peak von Mises stress was dramatically reduced when the Poisson,s ratio was increased from 0.05 to 0.49 under an axial compressive load of 1,400 N, while peak hydrostatic pressure was dramatically increased. Peak surface pressure was also increased with the Poisson's ratio, but with a lower magnitude compared to von Mises stress and hydrostatic pressure. In conclusion, the imaging process may cause 10 percent variations in peak contact stress, and the predicted stress distribution is sensitive to the accuracy of the material properties of the cartilage model, especially to the variation of Poisson's ratio.     

An EMG-driven biomechanical model of the canine cervical spine

Due to the frequency of cervical spine injuries in canines, the purpose of this effort was to develop an EMG-driven dynamic model of the canine cervical spine to assess a biomechanical understanding that enables one to investigate the risk of neck disorders. A canine subject was recruited in this investigation in order to collect subject specific data. Reflective markers and a motion capture system were used for kinematic measurement; surface electrodes were used to record electromyography signals, and with the aid of force plate kinetics were recorded. A 3D model of the canine subject was reconstructed from an MRI dataset. Muscles lines of action were defined through a new technique with the aid of 3D white light scanner. The model performed well with a 0.73 weighted R² value in all three planes. The weighted average absolute error of the predicted moment was less than 10% of the external moment. The proposed model is a canine specific forward-dynamics model that precisely tracks the canine subject head and neck motion, calculates the muscle force generated from the twelve major moment producing muscles, and estimates resulting loads on specific spinal tissues.     

A three-dimensional finite element model for arterial clamping

Clamp induced injuries of the arterial wall may determine the outcome of surgical procedures. Thus, it is important to investigate the underlying mechanical effects. We present a three-dimensional finite element model, which allows the study of the mechanical response of an artery-treated as a two-layer tube-during arterial clamping. The important residual stresses, which are associated with the load-free configuration of the artery, are also considered. In particular, the finite element analysis of the deformation process of a clamped artery and the associated stress distribution is presented. Within the clamping area a zone of axial tensile peak-stresses was identified, which (may) cause intimal and medial injury. This is an additional injury mechanism, which clearly differs from the commonly assumed wall damage occurring due to compression between the jaws of the clamp. The proposed numerical model provides essential insights into the mechanics of the clamping procedure and the associated injury mechanisms. It allows detailed parameter studies on a virtual clamped artery, which can not be performed with other methodologies. This approach has the potential to identify the most appropriate clamps for certain types of arteries and to guide optimal clamp design.     

Development of an accurate three-dimensional finite element knee model

This paper presents the development of a detailed articulating three-dimensional finite-element model of the human knee, derived from MRI scan images. The model utilises precise material models and many contact interfaces in order to produce a realistic kinematic response. The behaviour of the model was examined within two fields of biomechanical simulations: general life and car-crash. These simulations were performed with the non-linear explicit dynamic code PAM-SAFE trade mark. The knee model produced results that compared favourably with existing literature. Such a model (together with other joint models that could be constructed using the same techniques) would be a valuable tool for examining new designs of prosthesis and mechanisms of injury.     

Effect of mechanical loading on heterotopic ossification in cervical total disc replacement: a three-dimensional finite element analysis

The development of heterotopic ossification (HO) is considered one of the major complications following cervical total disc replacement (TDR). Even though previous studies have identified clinical and biomechanical conditions that may stimulate HO, the mechanism of HO formation has not been fully elucidated. The objective of this study is to investigate whether mechanical loading is a biomechanical condition that plays a substantial role to decide the HO formation. A finite element model of TDR on the C5–C6 was developed, and HO formation was predicted by simulating a bone adaptation process under various physiological mechanical loadings. The distributions of strain energy on vertebrae were assessed after HO formation. For the compressive force, most of the HO formation occurred on the vertebral endplates uncovered by the implant footplate which was similar to the Type 1 HO. For the anteriorly directed shear force, the HO was predominantly formed in the anterior parts of both the upper and lower vertebrae as the Type 2 HO. For both the flexion and extension moments, the HO shapes were similar to those for the shear force. The total strain energy was reduced after HO formation for all loading conditions. Two distinct types of HO were predicted based on mechanically induced bone adaptation processes, and our findings were consistent with those of previous clinical studies. HO formation might have a role in compensating for the non-uniform strain energy distribution which is one of the mechanical parameters related to the bone remodeling after cervical TDR.     

Investigation of effective intrusion and extrusion force for maxillary canine using finite element analysis

Abstract

Orthodontic tooth movement is mainly regulated by the biomechanical responses of loaded periodontal ligament (PDL). We investigated the effective intervals of orthodontic force in pure maxillary canine intrusion and extrusion referring to PDL hydrostatic stress and logarithmic strain. Finite element analysis (FEA) models, including a maxillary canine, PDL and alveolar bone, were constructed based on computed tomography (CT) images of a patient. The material properties of alveolar bone were non-uniformly defined using HU values of CT images; PDL was assumed to be a hyperelastic–viscoelastic material. The compressive stress and tensile stress ranging from 0.47 to 12.8?kPa and 18.8 to 51.2?kPa, respectively, were identified as effective for tooth movement; a strain 0.24% was identified as the lower limit of effective strain. The stress/strain distributions within PDL were acquired in canine intrusion and extrusion using FEA; root apex was the main force-bearing area in intrusion–extrusion movements and was more prone to resorption. Owing to the distinction of PDL biomechanical responses to compression and tension, the effective interval of orthodontic force was substantially lower in canine intrusion (80–90?g) than in canine extrusion (230–260?g). A larger magnitude of force remained applicable in canine extrusion. This study revised and complemented orthodontic biomechanical behaviours of tooth movement with intrusive–extrusive force and could further help optimize orthodontic treatment.     

A three-dimensional finite element analysis of the upper tibia

A three-dimensional finite element model of the proximal tibia has been developed to provide a base line for further modeling of prosthetic resurfaced tibiae. The geometry for the model was developed by digitizing coronal and transverse sections made with the milling machine, from one fresh tibia of average size. The load is equally distributed between the medial and lateral compartments over contact areas that were reported in the literature. An indentation test has been used to measure the stiffness and the ultimate strength of cancellous bone in four cadaver tibiae. These values provided the statistical basis for characterising the inhomogeneous distribution of the cancellous bone properties in the proximal tibia. All materials in the model were assumed to be linearly elastic and isotropic. Mechanical properties for the cortical bone and cartilage have been taken from the literature. Results have been compared with strain gage tests and with a two-dimensional axisymmetric finite element model both from the literature. Qualitative comparison between trabecular alignment, and the direction of the principal compressive stresses in the cancellous bone, showed a good relationship. Maximum stresses in the cancellous bone and cortical bone, under a load which occurs near stance phase during normal gait, show safety factors of approximately eight and twelve, respectively. The load sharing between the cancellous bone and the cortical bone has been plotted for the first 40 mm distally from the tibial eminence.     

A three-dimensional model of the human cervical spine for impact simulation

A three-dimensional analytical model of the cervical spine is described. The cervical vertebrae and the head are modeled as rigid bodies which are interconnected by deformable elements representing the intervertebral disks, facet joints, ligaments and muscles. A special pentahedral continuum element for representing the articular facets is described which effectively maintains stability of the cervical spine in both lateral and frontal plane accelerations, which is very difficult with multi-spring models of the facets. A simplified representation is used for the spine and body below the level of T1. The neck musculature is modeled by over 100 muscle elements representing 22 major muscle groups in the neck. The model has been validated for frontal and sideways impact accelerations by simulating published experimental data. Results are also presented to show the effects of the stretch reflex response on the dynamics of the head and neck under moderate acceleration.     

Finite element analysis of moment-rotation relationships for human cervical spine

A comprehensive, geometrically accurate, nonlinear C0-C7 FE model of head and cervical spine based on the actual geometry of a human cadaver specimen was developed. The motions of each cervical vertebral level under pure moment loading of 1.0 Nm applied incrementally on the skull to simulate the movements of the head and cervical spine under flexion, tension, axial rotation and lateral bending with the inferior surface of the C7 vertebral body fully constrained were analysed. The predicted range of motion (ROM) for each motion segment were computed and compared with published experimental data. The model predicted the nonlinear moment-rotation relationship of human cervical spine. Under the same loading magnitude, the model predicted the largest rotation in extension, followed by flexion and axial rotation, and least ROM in lateral bending. The upper cervical spines are more flexible than the lower cervical levels. The motions of the two uppermost motion segments account for half (or even higher) of the whole cervical spine motion under rotational loadings. The differences in the ROMs among the lower cervical spines (C3-C7) were relatively small. The FE predicted segmental motions effectively reflect the behavior of human cervical spine and were in agreement with the experimental data. The C0-C7 FE model offers potentials for biomedical and injury studies.     

Non-homogeneous analysis of three-dimensional transmural finite deformation in canine ventricular myocardium

A new method has been developed for analyzing transmural distributions of finite deformation in canine ventricular myocardium without the need to assume that the strain in a finite volume of the wall is homogeneous. The three-dimensional nodal geometric parameters of bilinear-cubic or bilinear-quadratic finite elements are fitted by least squares to the measured coordinates of 12-18 radiopaque markers implanted in the left ventricular free wall. For six dog hearts, root-mean-squared errors in the fitted in-plane coordinates ranged from 0.079-0.556 mm in the end-diastolic reference state and 0.142-0.622 mm at end-systole. The corresponding error ranges in the radial coordinate were 0.042-0.264 mm at end-diastole and 0.106-0.279 mm at end-systole. Smoothly continuous transmural profiles of wall strain computed as the element deformed during the cardiac cycle from end-diastole to end-systole showed good agreement with the discrete results of conventional homogeneous analysis. Using the kinematics of a thick-walled incompressible cylinder, overall absolute errors due to the non-homogeneity of myocardial deformation were found to be reduced in the new analysis by 30-35% for typical experimental parameters. Overall relative errors were also reduced (from 23 to 20%). Since measurement errors in the reconstructed marker coordinates were spatially smoothed by the fitting procedure, noise in the computed deformations was also substantially attenuated, and transmural gradients of three-dimensional strain components could be obtained with improved accuracy. Hence physiological factors affected by transmural stress and strain distributions, such as myocardial blood flow, ischemia and hypertrophy, may be better understood.     

A three-dimensional inverse finite element analysis of the heel pad

Quantification of plantar tissue behavior of the heel pad is essential in developing computational models for predictive analysis of preventive treatment options such as footwear for patients with diabetes. Simulation based studies in the past have generally adopted heel pad properties from the literature, in return using heel-specific geometry with material properties of a different heel. In exceptional cases, patient-specific material characterization was performed with simplified two-dimensional models, without further evaluation of a heel-specific response under different loading conditions. The aim of this study was to conduct an inverse finite element analysis of the heel in order to calculate heel-specific material properties in situ. Multidimensional experimental data available from a previous cadaver study by Erdemir et al. ("An Elaborate Data Set Characterizing the Mechanical Response of the Foot," ASME J. Biomech. Eng., 131(9), pp. 094502) was used for model development, optimization, and evaluation of material properties. A specimen-specific three-dimensional finite element representation was developed. Heel pad material properties were determined using inverse finite element analysis by fitting the model behavior to the experimental data. Compression dominant loading, applied using a spherical indenter, was used for optimization of the material properties. The optimized material properties were evaluated through simulations representative of a combined loading scenario (compression and anterior-posterior shear) with a spherical indenter and also of a compression dominant loading applied using an elevated platform. Optimized heel pad material coefficients were 0.001084 MPa (μ), 9.780 (α) (with an effective Poisson's ratio (ν) of 0.475), for a first-order nearly incompressible Ogden material model. The model predicted structural response of the heel pad was in good agreement for both the optimization (<1.05% maximum tool force, 0.9% maximum tool displacement) and validation cases (6.5% maximum tool force, 15% maximum tool displacement). The inverse analysis successfully predicted the material properties for the given specimen-specific heel pad using the experimental data for the specimen. The modeling framework and results can be used for accurate predictions of the three-dimensional interaction of the heel pad with its surroundings.     

A three-dimensional finite element model of the human anterior cruciate ligament: a computational analysis with experimental validation

In this study, the force and stress distribution within the anteromedial (AM) and posterolateral (PL) bundles of the anterior cruciate ligament (ACL) in response to an anterior tibial load with the knee at full extension was calculated using a validated three-dimensional finite element model (FEM) of a human ACL. The interaction between the AM and PL bundles, as well as the contact and friction caused by the ACL wrapping around the bone during knee motion, were included in the model. The AM and PL bundles of the ACL were simulated as incompressible homogeneous and isotropic hyperelastic materials. The multiple-degrees-of-freedom (DOF) knee kinematics of a cadaveric knee were first obtained using a robotic/universal force-moment sensor testing system. These data were used as the boundary conditions for the FEM of the ACL to calculate the forces in the ACL. The calculated forces were compared to the in situ force in the ACL, determined experimentally, to validate the model. The validated FEM was then used to calculate the force and stress distribution within the ACL under an anterior tibial load at full extension. The AM and PL bundles shared the force, and the stress distribution was non-uniform within both bundles with the highest stress localized near the femoral insertion site. The contact and friction caused by the ACL wrapping around the bone during knee motion played the role of transferring the force from the ACL to the bone, and had a direct effect on the force and stress distribution of the ACL. This validated model will enable the analysis of force and stress distribution in the ACL in response to more complex loading conditions and has the potential to help design improved surgical procedures following ACL injuries.     

A three-dimensional finite element model for the mechanics of cell-cell interactions

Technical challenges, including significant ones associated with cell rearrangement, have hampered the development of three-dimensional finite element models for the mechanics of embryonic cells. These challenges have been overcome by a new formulation in which the contents of each cell, assumed to have a viscosity mu, are modeled using a system of orthogonal dashpots. This approach overcomes a stiffening artifact that affects more traditional models, in which space-filling viscous elements are used to model the cytoplasm. Cells are assumed to be polyhedral in geometry, and each n-sided polygonal face is subdivided into n triangles with a common node at the face center so that it needs not remain flat. A constant tension gamma is assumed to act along each cell-cell interface, and cell rearrangements occur through one of two complementary topological transformations. The formulation predicts mechanical interactions between pairs of similar or dissimilar cells that are consistent with experiments, two-dimensional simulations, contact angle theory, and intracellular pressure calculations. Simulations of the partial engulfment of one tissue type by another show that the formulation is able to model aggregates of several hundred cells without difficulty. Simulations carried out using this formulation suggest new experimental approaches for measuring cell surface tensions and interfacial tensions. The formulation holds promise as a tool for gaining insight into the mechanics of isolated or aggregated embryonic cells and for the design and interpretation of experiments that involve them.     

Biomechanical consequences of subtalar joint arthroereisis in treating posterior tibial tendon dysfunction: a theoretical analysis using finite element analysis

Subtalar joint arthroereisis (SJA) has been introduced to control the hyperpronation in cases of flatfoot. The objective of this study is to evaluate the biomechanical consequence of SJA to restore the internal stress and load transfer to the intact state from the attenuated biomechanical condition induced by posterior tibial tendon dysfunction (PTTD). A three-dimensional finite element model of the foot and ankle complex was constructed based on clinical images of a healthy female (age 28 years, height 165 cm, body mass 54 kg). The boundary and loading condition during walking was acquired from the gait experiment of the model subject. Five sets of simulations (conditions) were completed: intact condition, mild PTTD, severe PTTD, mild PTTD with SJA, severe PTTD with SJA. The maximum von Mises stress of the metatarsal shafts and the load transfer along the midfoot during stance were analyzed. Generally, SJA deteriorated the joint force of the medial cuneonavicular and calcaneocuboid joints during late stance, while that of the metatarsocuneiform joints during early stance were over-corrected. Only the calcaneocuboid joint force at 45% stance demonstrated a trend of improvement. Besides, SJA exaggerated the increased stress of the metatarsals compared to the PTTD conditions, except that of the first metatarsal. Our study did not support the hypothesis that SJA can restore the internal load transfer and midfoot stress. SJA cannot compensate the salvage of midfoot stability attributed by PTTD and could be biomechanically insufficient to restore the biomechanical environment. Additional procedures such as orthotic intervention may be necessary.     

Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1–L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.