The Effect of Varying the Density-Modulus Relationship Used to Apply Material Properties in a Finite Element Model of the Distal Ulna

Author(s):  
Rebecca L. Austman ◽  
Jaques S. Milner ◽  
David W. Holdsworth ◽  
Cynthia A. Dunning

In many areas of orthopaedic biomechanics, such as implant design, properly developed Finite Element (FE) models can be a great companion to in-vitro studies, as they may allow a wider range of experimental variables to be explored in a cost-effective and timely manner. One challenge in developing these models is the assignment of accurate material properties to bone. Through the use of computed tomography (CT), many recent studies have developed subject-specific FE models, where material properties of bone are assigned based on density information derived from the scans. This involves the use of an equation to relate density and elastic modulus. There are several such relationships from which to choose in the literature. Most FE studies tend to use one of these multiple equations without justification or investigation into its appropriateness for the model.

Author(s):  
Emiliano P Ravera ◽  
Marcos J Crespo ◽  
Paola A Catalfamo Formento

Analysis of the human locomotor system using rigid-body musculoskeletal models has increased in the biomechanical community with the objective of studying muscle activations of different movements. Simultaneously, the finite element method has emerged as a complementary approach for analyzing the mechanical behavior of tissues. This study presents an integrative biomechanical framework for gait analysis by linking a musculoskeletal model and a subject-specific finite element model of the pelvis. To investigate its performance, a convergence study was performed and its sensitivity to the use of non-subject-specific material properties was studied. The total hip joint force estimated by the rigid musculoskeletal model and by the finite element model showed good agreement, suggesting that the integrative approach estimates adequately (in shape and magnitude) the hip total contact force. Previous studies found movements of up to 1.4 mm in the anterior–posterior direction, for single leg stance. These results are comparable with the displacement values found in this study: 0–0.5 mm in the sagittal axis. Maximum von Mises stress values of approximately 17 MPa were found in the pelvic bone. Comparing this results with a previous study of our group, the new findings show that the introduction of muscular boundary conditions and the flexion–extension movement of the hip reduce the regions of high stress and distributes more uniformly the stress across the pelvic bone. Thus, it is thought that muscle force has a relevant impact in reducing stresses in pelvic bone during walking of the finite element model proposed in this study. Future work will focus on including other deformable structures, such as the femur and the tibia, and subject-specific material properties.


Author(s):  
R L Austman ◽  
J S Milner ◽  
D W Holdsworth ◽  
C E Dunning

Assigning an appropriate density—modulus relationship is an important factor when applying inhomogeneous material properties to finite element models of bone. The purpose of this study was to develop a customized density—modulus equation for the distal ulna, using beam theory combined with experimental results. Five custom equations of the form E = aρb were used to apply material properties to models of eight ulnae. All equations passed through a point (1.85, Ec), where ρ = 1.85 g/cm3 represents the average density of cortical bone. For custom equations (1) to (3), Ec was predicted using beam theory, and the value of b was varied within the range reported in the literature. Custom equations (4) and (5) used other values of Ec from the literature, while keeping b constant. Results obtained from the custom equations were compared with those from other equations in the literature, and with experimental results. The beam theory analysis predicted Ec = 21 ± 1.6 GPa, and the three custom equations using this value tended to have the lowest errors. The power of the equations did not affect the results as much as the value used for Ec. Overall, a customized density—modulus relationship for the ulna was generated, which provided improved results over using previously reported density—modulus equations.


2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Sandeep Panwar Jogi ◽  
Rafeek Thaha ◽  
Sriram Rajan ◽  
Vidur Mahajan ◽  
Vasantha Kumar Venugopal ◽  
...  

Abstract Background Appropriate structural and material properties are essential for finite-element-modeling (FEM). In knee FEM, structural information could extract through 3D-imaging, but the individual subject’s tissue material properties are inaccessible. Purpose The current study's purpose was to develop a methodology to estimate the subject-specific stiffness of the tibiofemoral joint using finite-element-analysis (FEA) and MRI data of knee joint with and without load. Methods In this study, six Magnetic Resonance Imaging (MRI) datasets were acquired from 3 healthy volunteers with axially loaded and unloaded knee joint. The strain was computed from the tibiofemoral bone gap difference (ΔmBGFT) using the knee MR images with and without load. The knee FEM study was conducted using a subject-specific knee joint 3D-model and various soft-tissue stiffness values (1 to 50 MPa) to develop subject-specific stiffness versus strain models. Results Less than 1.02% absolute convergence error was observed during the simulation. Subject-specific combined stiffness of weight-bearing tibiofemoral soft-tissue was estimated with mean values as 2.40 ± 0.17 MPa. Intra-subject variability has been observed during the repeat scan in 3 subjects as 0.27, 0.12, and 0.15 MPa, respectively. All subject-specific stiffness-strain relationship data was fitted well with power function (R2 = 0.997). Conclusion The current study proposed a generalized mathematical model and a methodology to estimate subject-specific stiffness of the tibiofemoral joint for FEM analysis. Such a method might enhance the efficacy of FEM in implant design optimization and biomechanics for subject-specific studies. Trial registration The institutional ethics committee (IEC), Indian Institute of Technology, Delhi, India, approved the study on 20th September 2017, with reference number P-019; it was a pilot study, no clinical trail registration was recommended.


2011 ◽  
Vol 9 (67) ◽  
pp. 351-361 ◽  
Author(s):  
O. Panagiotopoulou ◽  
S. D. Wilshin ◽  
E. J. Rayfield ◽  
S. J. Shefelbine ◽  
J. R. Hutchinson

Finite element modelling is well entrenched in comparative vertebrate biomechanics as a tool to assess the mechanical design of skeletal structures and to better comprehend the complex interaction of their form–function relationships. But what makes a reliable subject-specific finite element model? To approach this question, we here present a set of convergence and sensitivity analyses and a validation study as an example, for finite element analysis (FEA) in general, of ways to ensure a reliable model. We detail how choices of element size, type and material properties in FEA influence the results of simulations. We also present an empirical model for estimating heterogeneous material properties throughout an elephant femur (but of broad applicability to FEA). We then use an ex vivo experimental validation test of a cadaveric femur to check our FEA results and find that the heterogeneous model matches the experimental results extremely well, and far better than the homogeneous model. We emphasize how considering heterogeneous material properties in FEA may be critical, so this should become standard practice in comparative FEA studies along with convergence analyses, consideration of element size, type and experimental validation. These steps may be required to obtain accurate models and derive reliable conclusions from them.


2008 ◽  
Vol 295 (3) ◽  
pp. H1141-H1149 ◽  
Author(s):  
Gaurav Krishnamurthy ◽  
Daniel B. Ennis ◽  
Akinobu Itoh ◽  
Wolfgang Bothe ◽  
Julia C. Swanson ◽  
...  

We measured leaflet displacements and used inverse finite-element analysis to define, for the first time, the material properties of mitral valve (MV) leaflets in vivo. Sixteen miniature radiopaque markers were sewn to the MV annulus, 16 to the anterior MV leaflet, and 1 on each papillary muscle tip in 17 sheep. Four-dimensional coordinates were obtained from biplane videofluoroscopic marker images (60 frames/s) during three complete cardiac cycles. A finite-element model of the anterior MV leaflet was developed using marker coordinates at the end of isovolumic relaxation (IVR; when the pressure difference across the valve is ∼0), as the minimum stress reference state. Leaflet displacements were simulated during IVR using measured left ventricular and atrial pressures. The leaflet shear modulus ( Gcirc-rad) and elastic moduli in both the commisure-commisure ( Ecirc) and radial ( Erad) directions were obtained using the method of feasible directions to minimize the difference between simulated and measured displacements. Group mean (±SD) values (17 animals, 3 heartbeats each, i.e., 51 cardiac cycles) were as follows: Gcirc-rad= 121 ± 22 N/mm2, Ecirc= 43 ± 18 N/mm2, and Erad= 11 ± 3 N/mm2( Ecirc> Erad, P < 0.01). These values, much greater than those previously reported from in vitro studies, may result from activated neurally controlled contractile tissue within the leaflet that is inactive in excised tissues. This could have important implications, not only to our understanding of mitral valve physiology in the beating heart but for providing additional information to aid the development of more durable tissue-engineered bioprosthetic valves.


Author(s):  
William J. Newman ◽  
Richard E. Debski ◽  
Susan M. Moore ◽  
Jeffrey A. Weiss

The shoulder is one of the most complex and often injured joints in the human body. The inferior glenohumeral ligament (IGHL), composed of the anterior band (AB), posterior band (PB) and the axillary pouch, has been shown to be an important contributor to anterior shoulder stability (Turkel, 1981). Injuries to the IGHL of the glenohumeral capsule are especially difficult to diagnose and treat effectively. The objective of this research was to develop a methodology for subject-specific finite element (FE) modeling of the ligamentous structures of the glenohumeral joint, specifically the IGHL, and to determine how changes in material properties affect predicted strains in the IGHL at 60° of external rotation. Using the techniques developed in this research, an improved understanding of the contribution of the IGHL to shoulder stability can be acquired.


Author(s):  
Sean M. Finley ◽  
J. Harley Astin ◽  
Evan Joyce ◽  
Andrew T. Dailey ◽  
Douglas L. Brockmeyer ◽  
...  

OBJECTIVE The underlying biomechanical differences between the pediatric and adult cervical spine are incompletely understood. Computational spine modeling can address that knowledge gap. Using a computational method known as finite element modeling, the authors describe the creation and evaluation of a complete pediatric cervical spine model. METHODS Using a thin-slice CT scan of the cervical spine from a 5-year-old boy, a 3D model was created for finite element analysis. The material properties and boundary and loading conditions were created and model analysis performed using open-source software. Because the precise material properties of the pediatric cervical spine are not known, a published parametric approach of scaling adult properties by 50%, 25%, and 10% was used. Each scaled finite element model (FEM) underwent two types of simulations for pediatric cadaver testing (axial tension and cardinal ranges of motion [ROMs]) to assess axial stiffness, ROM, and facet joint force (FJF). The authors evaluated the axial stiffness and flexion-extension ROM predicted by the model using previously published experimental measurements obtained from pediatric cadaveric tissues. RESULTS In the axial tension simulation, the model with 50% adult ligamentous and annulus material properties predicted an axial stiffness of 49 N/mm, which corresponded with previously published data from similarly aged cadavers (46.1 ± 9.6 N/mm). In the flexion-extension simulation, the same 50% model predicted an ROM that was within the range of the similarly aged cohort of cadavers. The subaxial FJFs predicted by the model in extension, lateral bending, and axial rotation were in the range of 1–4 N and, as expected, tended to increase as the ligament and disc material properties decreased. CONCLUSIONS A pediatric cervical spine FEM was created that accurately predicts axial tension and flexion-extension ROM when ligamentous and annulus material properties are reduced to 50% of published adult properties. This model shows promise for use in surgical simulation procedures and as a normal comparison for disease-specific FEMs.


2021 ◽  
Author(s):  
Sinan Yıldırım ◽  
Ufuk Çoban ◽  
Mehmet Çevik

Suspension linkages are one of the fundamental structural elements in each vehicle since they connect the wheel carriers i.e. axles to the body of the vehicle. Moreover, the characteristics of suspension linkages within a suspension system can directly affect driving safety, comfort and economics. Beyond these, all these design criteria are bounded to the package space of the vehicle. In last decades, suspension linkages have been focused on in terms of design development and cost reduction. In this study, a control arm of a diesel public bus was taken into account in order to get the most cost-effective design while improving the strength within specified boundary conditions. Due to the change of the supplier, the control arm of a rigid axle was redesigned to find an economical and more durable solution. The new design was analyzed first by the finite element analysis software Ansys and the finite element model of the control arm was validated by physical tensile tests. The outputs of the study demonstrate that the new design geometry reduces the maximum Von Mises stress 15% while being within the elastic region of the material in use and having found an economical solution in terms of supplier’s criteria.


Sign in / Sign up

Export Citation Format

Share Document