Predicting the yield of the proximal femur using high-order finite-element analysis with inhomogeneous orthotropic material properties

2010 ◽  
Vol 368 (1920) ◽  
pp. 2707-2723 ◽  
Author(s):  
Zohar Yosibash ◽  
David Tal ◽  
Nir Trabelsi
Keyword(s):  
Finite Element ◽  
Proximal Femur ◽  
Material Properties ◽  
Orthotropic Material ◽  
Isotropic Material ◽  
High Order ◽  
Principal Strain ◽  
Fe Model ◽  
Yield Criteria ◽  
Maximum Principal Strain

High-order finite-element (FE) analyses with inhomogeneous isotropic material properties have been shown to predict the strains and displacements on the surface of the proximal femur with high accuracy when compared with in vitro experiments. The same FE models with inhomogeneous orthotropic material properties produce results similar to those obtained with isotropic material properties. Herein, we investigate the yield prediction capabilities of these models using four different yield criteria, and the spread in the predicted load between the isotropic and orthotropic material models. Subject-specific high-order FE models of two human femurs were generated from CT scans with inhomogeneous orthotropic or isotropic material properties, and loaded by a simple compression force at the head. Computed strains and stresses by both the orthotropic and isotropic FE models were used to determine the load that predicts ‘yielding’ by four different ‘yield criteria’: von Mises, Drucker–Prager, maximum principal stress and maximum principal strain. One of the femurs was loaded by a simple load until fracture, and the force resulting in yielding was compared with the FE predicted force. The surface average of the ‘maximum principal strain’ criterion in conjunction with the orthotropic FE model best predicts both the yield force and fracture location compared with other criteria. There is a non-negligible influence on the predictions if orthotropic or isotropic material properties are applied to the FE model. All stress-based investigated ‘yield criteria’ have a small spread in the predicted failure. Because only one experiment was performed with a rather simplified loading configuration, the conclusions of this work cannot be claimed to be either reliable or sufficient, and future experiments should be performed to further substantiate the conclusions.

Patient-Specific Finite-Element Analyses of the Proximal Femur with Orthotropic Material Properties Validated by Experiments

2011 ◽  
Vol 133 (6) ◽  
Author(s):  
Nir Trabelsi ◽  
Zohar Yosibash
Keyword(s):  
Finite Element ◽  
Material Properties ◽  
Orthotropic Material ◽  
High Order ◽  
Patient Specific ◽  
Pronounced Difference ◽  
Quantitative Computer Tomography ◽  
Fresh Frozen ◽  
Orthotropic Properties

Patient-specific high order finite-element (FE) models of human femurs based on quantitative computer tomography (QCT) with inhomogeneous orthotropic and isotropic material properties are addressed. The point-wise orthotropic properties are determined by a micromechanics (MM) based approach in conjunction with experimental observations at the osteon level, and two methods for determining the material trajectories are proposed (along organs outer surface, or along principal strains). QCT scans on four fresh-frozen human femurs were performed and high-order FE models were generated with either inhomogeneous MM-based orthotropic or empirically determined isotropic properties. In vitro experiments were conducted on the femurs by applying a simple stance position load on their head, recording strains on femurs’ surface and head’s displacements. After verifying the FE linear elastic analyses that mimic the experimental setting for numerical accuracy, we compared the FE results to the experimental observations to identify the influence of material properties on models’ predictions. The strains and displacements computed by FE models having MM-based inhomogeneous orthotropic properties match the FE-results having empirically based isotropic properties well, and both are in close agreement with the experimental results. When only the strains in the femoral neck are being compared a more pronounced difference is noticed between the isotropic and orthotropic FE result. These results lay the foundation for applying more realistic inhomogeneous orthotropic material properties in FEA of femurs.

Age Dependent and Anisotropic Material Properties of Immature Porcine Sclera

2013 ◽  
Author(s):  
Jami M. Saffioti ◽  
Brittany Coats
Keyword(s):  
Finite Element ◽  
Material Properties ◽  
Anisotropic Material ◽  
Computational Models ◽  
Fe Model ◽  
Ocular Tissues ◽  
Load Bearing ◽  
Anisotropic Properties ◽  
Age Dependent ◽  
Testing Protocols

Current finite element (FE) models of the pediatric eye are based on adult material properties [2,3]. To date, there are no data characterizing the age dependent material properties of ocular tissues. The sclera is a major load bearing tissue and an essential component to most computational models of the eye. In preparation for the development of a pediatric FE model, age-dependent and anisotropic properties of sclera were evaluated in newborn (3–5 days) and toddler (4 weeks) pigs. Data from this study will guide future testing protocols for human pediatric specimens.

Fracture Prediction for the Proximal Femur Using Finite Element Models: Part I—Linear Analysis

1991 ◽  
Vol 113 (4) ◽  
pp. 353-360 ◽  
Author(s):  
J. C. Lotz ◽  
E. J. Cheal ◽  
W. C. Hayes
Keyword(s):  
Finite Element ◽  
Fracture Risk ◽  
Strain Gage ◽  
Proximal Femur ◽  
Material Properties ◽  
The United States ◽  
Nonlinear Material ◽  
Finite Element Models ◽  
Model Predictions

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. However, for a simulated fall, the peak stresses were in the intertrochanteric region. The Ward’s triangle (basicervical) site commonly used for the clinical assessment of osteoporosis was not heavily loaded in either situation. These findings suggest that the intertrochanteric region may be the most sensitive site for the assessment of fracture risk due to a fall and the subcapital region for fracture risk due to repetitive activities such as walking.

FEM-Based Simulation of Temperature in Speed Stroke Grinding with 3D Transient Moving Heat Sources

2011 ◽  
Vol 223 ◽  
pp. 733-742 ◽  
Author(s):  
Barbara Linke ◽  
Michael Duscha ◽  
Anh Tuan Vu ◽  
Fritz Klocke
Keyword(s):  
Finite Element ◽  
Material Properties ◽  
Material Removal ◽  
Structural Changes ◽  
Numerical Technique ◽  
Temperature Fields ◽  
Fe Model ◽  
Measured Temperature ◽  
Temperature Dependent ◽  
Removal Rates

The grinding process is one of the most important finishing processes to obtain high surface quality. Nowadays, grinding is also considered as a high performance process with high material removal rates. Nevertheless, to avoid thermally-induced structural changes poses a major challenge for this manufacturing technology. Until now, the Finite Element Method (FEM) has been widely applied as a proper numerical technique to predict workpiece properties in machining processes. However, actual models in grinding are limited to conventional grinding processes with simple workpiece profiles and low table speeds. In this paper, finite element simulations are expanded to 3-dimensional (3D) models with temperature-dependent material properties and heat source profiles derived from experimental results, i.e. tangential forces. Both temperature simulation and measurement were conducted for deep grinding, pendulum grinding and speed stroke grinding in the table speed range of vw= 12 m/min to 180 m/min and specific material removal rates of Q’w= 40 mm³/mms. Overall, the simulation results show a good agreement with the measured temperature and surface integrity after grinding. This research indicates that a 3D FE model with temperature dependent material properties can predict realistic temperature fields in speed stroke grinding. Therefore, the experiment and measurement costs and time can be reduced by FEM simulation.

scholarly journals Finite element modeling of proximal femur with quantifiable weight-bearing area in standing position

2020 ◽  
Author(s):  
Yang Peng ◽  
Tian-Ye Lin ◽  
Jing-Li Xu ◽  
Hui-Yu Zeng ◽  
Da Chen ◽  
...  
Keyword(s):  
Finite Element ◽  
Femoral Head ◽  
Proximal Femur ◽  
Physical Phenomenon ◽  
Weight Bearing ◽  
Positional Distribution ◽  
Standing Position ◽  
Fe Model ◽  
X Ray ◽  
Weight Bearing Area

Abstract BackgroundThe positional distribution and size of the weight-bearing area of femoral head in the standing position as well as the direct active surface of joint force can directly affect the result of finite element (FE) stress analysis, however in most studies related separate FE models of femur, the division of this area is vague, imprecise and un-individualized. The purpose of this study was to quantify the positional distribution and size of the weight-bearing area of femoral head in standing position by a set of simple methods, to realize individualized reconstruction of proximal femur FE model.MethodsFive adult volunteers were recruited for X-ray and CT examination in the same simulated bipedal standing position with a specialized patented device. We extracted these image data, calculated the 2D weight-bearing area on X-ray image, reconstructed the 3D model of proximal femur based on CT data, and registered them to realize the 2D weight-bearing area to 3D transformation as the quantified weight-bearing surface. One of the 3D models of proximal femur was randomly selected for finite element analysis (FEA), and we defined three different loading surfaces, and compared their FEA results.ResultsA total of 10 weight-bearing surfaces in 5 volunteers were constructed, they were mainly distributed on the dome and anterolateral of femoral head with crescent shape, in the range of 1,218.63mm2 - 1,871.06mm2. The results of FEA showed stress magnitude and distribution in proximal femur FE models among three different loading conditions were significant differences, the loading case with quantized weight-bearing area was more in accordance with the physical phenomenon of the hip.ConclusionThis study confirmed an effective FE modeling method of proximal femur, which can quantify weight-bearing area to define more reasonable load surface setting without increasing the actual modeling difficulty.

Automated Segmentation of Swine Skulls for Finite Element Model Creation Using High Resolution μ-CT Images

2019 ◽  
Vol 16 (03) ◽  
pp. 1842012 ◽  
Author(s):  
Zimo Zhu ◽  
Donna C. Jones ◽  
G. R. Liu ◽  
Sajjad Soleimani ◽  
Xu Huang ◽  
...  
Keyword(s):  
Finite Element ◽  
Material Properties ◽  
Complex Structure ◽  
Model Development ◽  
Element Model ◽  
Fe Model ◽  
Micro Computed Tomography ◽  
Mineral Density ◽  
Detailed Model ◽  
Bone Mineral Density Loss

Finite element (FE) analysis has been widely used to investigate bone responses to mechanical loading. Research in long bones has been straight forward because modeling of these bones requires only two material properties. Such an FE model may provide an adequate approximation of the anatomy for many cases. However, a more detailed model of skull bones is needed to accurately capture its complex structure of multiple bone pieces and the various mineral densities distributed throughout these bone pieces. Unfortunately, FE model development incorporating both complex geometries and anatomically accurate material properties is both computationally and labor intensive. In this study, a method is proposed to automatically segment micro-computed tomography ([Formula: see text]-CT) scan images of bone pieces to build an FE model of a full swine hemi-skull. Using the Digital Imaging and Communications in Medicine (DICOM) files from scanned bones, the complete geometry of each bone piece is recreated through seven customized processing algorithms. After assembling the bone pieces to form the skull, experimentally derived Young’s modulus values are correlated to grayscale values to produce a detailed FE model for accurate simulation. This detailed skull model can be used to predict strain/stress patterns in response to various loading regimes to facilitate research questions in fracture healing and growth, as well as bone tissue engineering and bone mineral density loss (e.g., osteoporosis).

scholarly journals Inverse Calculation of Timber-CFRP Composite Beams Using Finite Element Analysis

2020 ◽  
Author(s):  
Khaled Saad ◽  
András Lengyel
Keyword(s):  
Finite Element Analysis ◽  
Finite Element ◽  
Material Properties ◽  
Orthotropic Material ◽  
Material Model ◽  
Flexural Behavior ◽  
Fiber Reinforced Polymers ◽  
Element Analysis ◽  
Linear Elastic ◽  
Timber Beams

This study focuses on the flexural behavior of timber beams externally reinforced using carbon fiber-reinforced polymers (CFRP). Linear and non-linear finite element analysis were proposed and validated by experimental tests carried out on 44 timber beams to inversely determine the material properties of the timber and the CFRP. All the beams have the same geometrical properties and were loaded under four points bending. In this paper the general commercial software ANSYS was used, and three- and two-dimensional numerical models were evaluated for their ability to describe the behavior of the solid timber beams. The linear elastic orthotropic material model was assumed for the timber beams in the linear range and the 3D nonlinear rate-independent generalized anisotropic Hill potential model was assumed to describe the nonlinear behavior of the material. As for the CFRP, a linear elastic orthotropic material model was introduced for the fibers and a linear elastic isotropic model for the epoxy resin. No mechanical model was introduced to describe the interaction between the timber and the CFRP since failure occurred in the tensile zone of the wood. Simulated and measured load-mid-span deflection responses were compared and the material properties for timber-CFRP were numerically determined.

The Biomechanics of Human Femurs in Axial and Torsional Loading: Comparison of Finite Element Analysis, Human Cadaveric Femurs, and Synthetic Femurs

2006 ◽  
Vol 129 (1) ◽  
pp. 12-19 ◽  
Author(s):  
M. Papini ◽  
R. Zdero ◽  
E. H. Schemitsch ◽  
P. Zalzal
Keyword(s):  
Finite Element ◽  
Bone Quality ◽  
Material Properties ◽  
Element Model ◽  
Torsional Stiffness ◽  
Fe Model ◽  
Bone Stock ◽  
Torsional Loading ◽  
Input Material ◽  
Geometric Size

To assess the performance of femoral orthopedic implants, they are often attached to cadaveric femurs, and biomechanical testing is performed. To identify areas of high stress, stress shielding, and to facilitate implant redesign, these tests are often accompanied by finite element (FE) models of the bone/implant system. However, cadaveric bone suffers from wide specimen to specimen variability both in terms of bone geometry and mechanical properties, making it virtually impossible for experimental results to be reproduced. An alternative approach is to utilize synthetic femurs of standardized geometry, having material behavior approximating that of human bone, but with very small specimen to specimen variability. This approach allows for repeatable experimental results and a standard geometry for use in accompanying FE models. While the synthetic bones appear to be of appropriate geometry to simulate bone mechanical behavior, it has not, however, been established what bone quality they most resemble, i.e., osteoporotic or osteopenic versus healthy bone. Furthermore, it is also of interest to determine whether FE models of synthetic bones, with appropriate adjustments in input material properties or geometric size, could be used to simulate the mechanical behavior of a wider range of bone quality and size. To shed light on these questions, the axial and torsional stiffness of cadaveric femurs were compared to those measured on synthetic femurs. A FE model, previously validated by the authors to represent the geometry of a synthetic femur, was then used with a range of input material properties and change in geometric size, to establish whether cadaveric results could be simulated. Axial and torsional stiffnesses and rigidities were measured for 25 human cadaveric femurs (simulating poor bone stock) and three synthetic “third generation composite” femurs (3GCF) (simulating normal healthy bone stock) in the midstance orientation. The measured results were compared, under identical loading conditions, to those predicted by a previously validated three-dimensional finite element model of the 3GCF at a variety of Young’s modulus values. A smaller FE model of the 3GCF was also created to examine the effects of a simple change in bone size. The 3GCF was found to be significantly stiffer (2.3 times in torsional loading, 1.7 times in axial loading) than the presently utilized cadaveric samples. Nevertheless, the FE model was able to successfully simulate both the behavior of the 3GCF, and a wide range of cadaveric bone data scatter by an appropriate adjustment of Young’s modulus or geometric size. The synthetic femur had a significantly higher stiffness than the cadaveric bone samples. The finite element model provided a good estimate of upper and lower bounds for the axial and torsional stiffness of human femurs because it was effective at reproducing the geometric properties of a femur. Cadaveric bone experiments can be used to calibrate FE models’ input material properties so that bones of varying quality can be simulated.

Identification of printed circuit boards mechanical properties using response surface methods

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Mohammad Gharaibeh
Keyword(s):  
Finite Element ◽  
Response Surface ◽  
Material Properties ◽  
Printed Circuit Boards ◽  
Dynamic Testing ◽  
Mode Shapes ◽  
Fe Model ◽  
Circuit Boards ◽  
Content Type ◽  
Printed Circuit

Purpose This study aims to discuss the determination of the unknown in-plane mechanical material properties of printed circuit boards (PCBs) by correlating the results from dynamic testing and finite element (FE) models using the response surface method (RSM). Design/methodology/approach The first 10 resonant frequencies and vibratory mode shapes are measured using modal analysis with hammer testing experiment, and hence, systematically compared with finite element analysis (FEA) results. The RSM is consequently used to minimize the cumulative error between dynamic testing and FEA results by continuously modifying the FE model, to acquire material properties of PCBs. Findings Great agreement is shown when comparing FEA to measurements, the optimum in-plane material properties were identified, and hence, verified. Originality/value This paper used FEA and RSMs along with modal measurements to obtain in-plane material properties of PCBs. The methodology presented here can be easily generalized and repeated for different board designs and configurations.

A high-order control volume finite element method for thermoelastic analysis of functionally graded solids with mixed grids

2018 ◽  
Vol 233 (11) ◽  
pp. 3994-4013
Author(s):  
Qi Liu ◽  
Yan Yu ◽  
Pingjian Ming
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Material Properties ◽  
Control Volume ◽  
Functionally Graded ◽  
High Order ◽  
Shape Functions ◽  
Quadrilateral Element ◽  
Thermoelastic Analysis ◽  
Control Volume Finite Element

In this article, a new two-dimensional control volume finite element method has been developed for thermoelastic analysis in functionally graded materials. A nine-node quadrilateral element and a six-node triangular element are employed to deal with the mixed-grid problem. The unknown variables and material properties are defined at the node. The high-order shape functions of six-node triangular and nine-node quadrilateral element are employed to obtain the unknown variables and their derivatives. In addition, the material properties in functionally graded structure are also modeled by applying the high-order shape functions. The capabilities of the presented method to heat conduction problem, elastic problem, and thermoelastic problem have been validated. First, the defined location of material properties is found to be important for the accuracy of the numerical results. Second, the presented method is proven to be efficient and reliable for the elastic analysis in multi-phase materials. Third, the presented method is capable of high-order mixed grids. The memory and computational costs of the presented method are also compared with other numerical methods.

Sign in / Sign up

Export Citation Format

Share Document