Development and numerical validation of a finite element model of the muscle standardized femur

Author(s):  
K Polgar ◽  
H S Gill ◽  
M Viceconti ◽  
D W Murray ◽  
J J O'Connor

The human femur is one of the parts of the musculo-skeletal system most frequently analysed by means of the finite element (FE) method. Most FE studies of the human femur are based on computed tomography data sets of a particular femur. Since the geometry of the chosen sample anatomy influences the computed results, direct comparison across various models is often difficult or impossible. The aim of the present work was to develop and validate a novel three-dimensional FE model of the human femur based on the muscle standardized femur (MuscleSF) geometry. In the new MuscleSF FE model, the femoral attachment of each muscle was meshed separately on the external bone surface. The model was tested under simple load configurations and the results showed good agreement with the converged solution of a former study. In the future, using the validated MuscleSF FE model for numerical studies of the human femur will provide the following benefits: (a) the numerical accuracy of the model is known; (b) muscle attachment areas are incorporated in the model, therefore physiological loading conditions can be easily defined; (c) analyses of the femur under physiological load cases will be replicable; (d) results based on different load configurations could be compared across various studies.

2015 ◽  
Vol 23 (9) ◽  
pp. 1480-1494 ◽  
Author(s):  
Zhaoye Qin ◽  
Delin Cui ◽  
Shaoze Yan ◽  
Fulei Chu

Due to frictional slippage between the joint components, clamp band joints may generate nonlinear stiffness and friction damping, which will affect the dynamics of the joint structures. Accurate modeling of the frictional behavior in clamp band joints is crucial for reliable estimation of the joint structure dynamics. While the finite element (FE) method is a powerful tool to analyze structures assembled with joints, it is computationally expensive and inefficient to perform transient analyses with three-dimensional (3D) FE models involving contact nonlinearity. In this paper, a two-dimensional (2D) FE model of much more efficiency is applied to investigate the dynamics of a clamp band jointed structure subjected to longitudinal base excitations. Prior to dynamic analyses, the sources of the model inaccuracy are determined, upon which a two-step model updating technique is proposed to improve the accuracy of the 2D model in accordance with the quasi-static test data. Then, based on the updated 2D model, the nonlinear influence of the clamp band joint on the dynamic response of the joint structure is investigated. Sine-sweep tests are carried out to validate the updated 2D FE model. The FE modeling and updating techniques proposed here can be applied to other types of structures of cyclic symmetry to develop accurate model with high computational efficiency.


2000 ◽  
Author(s):  
Subramanya Uppala ◽  
Robert X. Gao ◽  
Scott Cowan ◽  
K. Francis Lee

Abstract The strength and stability of the lumbar spine are determined not only by the bone and muscles, but also by the visco-elastic structures and the interplay between the different components of the spine, such as ligaments, capsules, annulus fibrosis, and articular cartilage. In this paper we present a non-linear three-dimensional Finite Element model of the lumbar spine. Specifically, a three-dimensional FE model of the L4-5 one-motion segment/2 vertebrae was developed. The cortical shell and the cancellous bone of the vertebral body were modeled as 3D isoparametric eight-nodal elements. Finite element models of spinal injuries with fixation devices are also developed. The deformations across the different sections of the spine are observed under the application of axial compression, flexion/extension, and lateral bending. The developed FE models provided input to both the fixture design and experimental studies.


2019 ◽  
Vol 281 ◽  
pp. 01006 ◽  
Author(s):  
Majid M.A. Kadhim ◽  
Mohammed J Altaee ◽  
Ali Hadi Adheem ◽  
Akram R. Jawdhari

Fibre reinforced cementitious matric (FRCM) is a recent application of fibre reinforced polymer (FRP) reinforcement, developed to overcome several limitations associated with the use of organic adhesive [e.g. epoxies] in FRPs. It consists of two dimensional FRP mesh saturated with a cement mortar, which is inorganic in nature and compatible with concrete and masonry substrates. In this study, a robust three-dimensional (3D) finite element (FE) model has been developed to study the behaviour of slender reinforced concrete columns confined by FRCM jackets, and loaded concentrically and eccentrically. The model accounts for material nonlinearities in column core and cement mortar, composite failure of FRP mesh, and global buckling. The model response was validated against several laboratory tests from literature, comparing the ultimate load, load-lateral deflection and failure mode. Maximum divergence between numerical and experimental results was 12%. Following the validation, the model will be used later in a comprehensive parametric analysis to gain a profound knowledge of the strengthening system, and examine the effects of several factors expected to influence the behaviour of confined member.


2019 ◽  
Vol 43 (4) ◽  
pp. 443-453
Author(s):  
Stephen M. Handrigan ◽  
Sam Nakhla

An investigation to determine the effect of porosity concentration and location on elastic modulus is performed. Due to advancements in testing methods, the manufacturing and testing of microbeams to obtain mechanical response is possible through the use of focused ion beam technology. Meanwhile, rigorous analysis is required to enable accurate extraction of the elastic modulus from test data. First, a one-dimensional investigation with beam theory, Euler–Bernoulli and Timoshenko, was performed to estimate the modulus based on load-deflection curve. Second, a three-dimensional finite element (FE) model in Abaqus was developed to identify the effect of porosity concentration. Furthermore, the current work provided an accurate procedure to enable accurate extraction of the elastic modulus from load-deflection data. The use of macromodels such as beam theory and three-dimensional FE model enabled enhanced understanding of the effect of porosity on modulus.


2014 ◽  
Vol 2014 ◽  
pp. 1-9 ◽  
Author(s):  
Namkeun Kim ◽  
You Chang ◽  
Stefan Stenfelt

A three-dimensional finite-element (FE) model of a human dry skull was devised for simulation of human bone-conduction (BC) hearing. Although a dry skull is a simplification of the real complex human skull, such model is valuable for understanding basic BC hearing processes. For validation of the model, the mechanical point impedance of the skull as well as the acceleration of the ipsilateral and contralateral cochlear bone was computed and compared to experimental results. Simulation results showed reasonable consistency between the mechanical point impedance and the experimental measurements when Young’s modulus for skull and polyurethane was set to be 7.3 GPa and 1 MPa with 0.01 and 0.1 loss factors at 1 kHz, respectively. Moreover, the acceleration in the medial-lateral direction showed the best correspondence with the published experimental data, whereas the acceleration in the inferior-superior direction showed the largest discrepancy. However, the results were reasonable considering that different geometries were used for the 3D FE skull and the skull used in the published experimental study. The dry skull model is a first step for understanding BC hearing mechanism in a human head and simulation results can be used to predict vibration pattern of the bone surrounding the middle and inner ear during BC stimulation.


2017 ◽  
Vol 36 (2) ◽  
pp. 35-45
Author(s):  
Henry M. Kiwelu

Experiments were performed on scaled glue laminated bending specimens to observetime dependent development of deformations during drying and wetting. Measurementsdetermined changes in the average moisture content and external shape and dimensionsbetween when specimens were placed into constant or variable climates. Alterations inthe external shape and dimensions reflected changes in the average value anddistribution of moisture and mechanosorptive creep in the glulam. The results are beingused to develop a sequentially-coupled three-dimensional hygrothermal Finite Element(FE) model for predicting temporally varying internal strains and external deformationsof drying or wetting solid wood structural components. The model implies temporallyvarying, and eventual steady, state internal stress distributions in members based onelastic and creep compliances that represent wood within glulam as a continuousorthotropic homogenised material. Thus, predictions are consistent with smearedengineering stress analysis methods rather than being a physically correct analogue ofhow solid wood behaves. This paper discusses limitations of and intended improvementsto the FE modelling. Complementary investigations are underway to address otheraspects of the hygrothermal behaviour of structural members of wood and othermaterials (e.g. reinforced concrete) embedded within superstructure frameworks ofmulti-storey hybrid buildings.


Author(s):  
Jingwen Hu ◽  
Zhigang Li ◽  
Jinhuan Zhang

Head injury is the leading cause of pediatric fatality and disability in the United States (1). Although finite element (FE) method has been widely used for investigating head injury under impact, there are only a few 3D pediatric head FE models available in the literature, including a 6-month-old child head model developed by Klinich et al (2), a newborn, a 6-month-old and a 3-year-old child head model developed by Roth et al. (3, 4, 5), and a 1.5-month-old infant head model developed by Coats et al (6). Each of these models only represents a head at a single age with single head geometry. Nowadays, population-based simulations are getting more and more attention. In population-based injury simulations, impact responses for not only an individual but also a group of people can be predicted, which takes into account variations among people thus providing more realistic predictions. However, a parametric pediatric head model capable of simulating head responses for different children at different ages is currently not available. Therefore, the objective of this study is to develop a fast and efficient method to build pediatric head FE models with different head geometries and skull thickness distributions. The method was demonstrated by morphing a 6-month-old infant head FE model into three newborn infant head FE models and by validating three morphed head models against limited cadaveric test data.


1995 ◽  
Vol 117 (3) ◽  
pp. 272-278 ◽  
Author(s):  
M. Dalstra ◽  
R. Huiskes ◽  
L. van Erning

Due to both its shape and its structural architecture, the mechanics of the pelvic bone are complex. In Finite Element (FE) models, these aspects have often been (over) simplified, sometimes leading to conclusions which did not bear out in reality. The purpose of this study was to develop a more realistic FE model of the pelvic bone. This not only implies that the model has to be three-dimensional, but also that the thickness of the cortical shell and the density distribution of the trabecular bone throughout the pelvic bone have to be incorporated in the model in a realistic way. For this purpose, quantitative measurements were performed on computer tomography scans of several pelvic bones, after which the measured quantities were allocated to each element of the mesh individually. To validate this FE model, two fresh pelvic bones were fitted with strain gages and loaded in a testing machine. Stresses calculated from the strain data of this experiment were compared to the results of a simulation with the developed pelvic FE model.


2005 ◽  
Vol 288 (4) ◽  
pp. H1844-H1850 ◽  
Author(s):  
Alan B. C. Dang ◽  
Julius M. Guccione ◽  
Jacob M. Mishell ◽  
Peng Zhang ◽  
Arthur W. Wallace ◽  
...  

Infarcted segments of myocardium demonstrate functional impairment ranging in severity from hypokinesis to dyskinesis. We sought to better define the contributions of passive material properties (stiffness) and active properties (contracting myocytes) to infarct thickening. Using a finite-element (FE) model, we tested the hypothesis that infarcted myocardium must contain contracting myocytes to be akinetic and not dyskinetic. A three-dimensional FE mesh of the left ventricle was developed with echocardiographs from a reperfused ovine anteroapical infarct. The nonlinear stress-strain relationship for the diastolic myocardium was anisotropic with respect to the local muscle fiber direction, and an elastance model for active fiber stress was incorporated. The diastolic stiffness ( C) and systolic material property (isometric tension at longest sarcomere length and peak intracellular calcium concentration, Tmax) of the uninfarcted remote myocardium were assumed to be normal ( C = 0.876 kPa, Tmax = 135.7 kPa). Diastolic and systolic properties of the infarct necessary to produce akinesis, defined as an average radial strain between −0.01 and 0.01, were determined by assigning a range of diastolic stiffnesses and scaling infarct Tmax to represent the percentage of contracting myocytes between 0% and 100%. As C was increased to 11 times normal ( C = 10 kPa) the percentage of Tmax necessary for akinesis increased from 20% to 50%. Without contracting myocytes, C = 250 kPa was necessary to achieve akinesis. If infarct stiffness is <285 times normal, contracting myocytes are required to prevent dyskinetic infarct wall motion.


2017 ◽  
Vol 2017 ◽  
pp. 1-8 ◽  
Author(s):  
Siti Nurfaezah Zahari ◽  
Mohd Juzaila Abd Latif ◽  
Nor Raihanah Abdull Rahim ◽  
Mohammed Rafiq Abdul Kadir ◽  
Tunku Kamarul

The present study was conducted to examine the effects of body weight on intradiscal pressure (IDP) and annulus stress of intervertebral discs at lumbar spine. Three-dimensional finite element model of osseoligamentous lumbar spine was developed subjected to follower load of 500 N, 800 N, and 1200 N which represent the loads for individuals who are normal and overweight with the pure moments at 7.5 Nm in flexion and extension motions. It was observed that the maximum IDP was 1.26 MPa at L1-L2 vertebral segment. However, the highest increment of IDP was found at L4-L5 segment where the IDP was increased to 30% in flexion and it was more severe at extension motion reaching to 80%. Furthermore, the maximum annulus stress also occurred at the L1-L2 segment with 3.9 MPa in extension motion. However, the highest increment was also found at L4-L5 where the annulus stress increased to 17% in extension motion. Based on these results, the increase of physiological loading could be an important factor to the increment of intradiscal pressure and annulus fibrosis stress at all intervertebral discs at the lumbar spine which may lead to early intervertebral disc damage.


Sign in / Sign up

Export Citation Format

Share Document