Fracture Prediction for the Proximal Femur Using Finite Element Models: Part II—Nonlinear Analysis

1991 ◽  
Vol 113 (4) ◽  
pp. 361-365 ◽  
Author(s):  
J. C. Lotz ◽  
E. J. Cheal ◽  
W. C. Hayes

In Part I we reported the results of linear finite element models of the proximal femur generated using geometric and constitutive data collected with quantitative computed tomography. These models demonstrated excellent agreement with in vitro studies when used to predict ultimate failure loads. In Part II, we report our extension of those finite element models to include nonlinear behavior of the trabecular and cortical bone. A highly nonlinear material law, originally designed for representing concrete, was used for trabecular bone, while a bilinear material law was used for cortical bone. We found excellent agreement between the model predictions and in vitro fracture data for both the onset of bone yielding and bone fracture. For bone yielding, the model predictions were within 2 percent for a load which simulated one-legged stance and 1 percent for a load which simulated a fall. For bone fracture, the model predictions were within 1 percent and 17 percent, respectively. The models also demonstrated different fracture mechanisms for the two different loading configurations. For one-legged stance, failure within the primary compressive trabeculae at the subcapital region occurred first, leading to load transfer and, ultimately, failure of the surrounding cortical shell. However, for a fall, failure of the cortical and trabecular bone occurred simultaneously within the intertrochanteric region. These results support our previous findings that the strength of the subcapital region is primarily due to trabecular bone whereas the strength of the intertrochanteric region is primarily due to cortical bone.

1991 ◽  
Vol 113 (4) ◽  
pp. 353-360 ◽  
Author(s):  
J. C. Lotz ◽  
E. J. Cheal ◽  
W. C. Hayes

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.


Bone ◽  
2012 ◽  
Vol 51 (4) ◽  
pp. 737-740 ◽  
Author(s):  
Janne E.M. Koivumäki ◽  
Jérôme Thevenot ◽  
Pasi Pulkkinen ◽  
Volker Kuhn ◽  
Thomas M. Link ◽  
...  

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Long Liu ◽  
Lifeng Wang ◽  
Ziwang Xiao

PurposeThe flexural reinforcement of bridges in-service has been an important research field for a long time. Anchoring steel plate at the bottom of beam is a simple and effective method to improve its bearing capacity. The purpose of this paper is to explore the influence of anchoring steel plates of different thicknesses on the bearing capacity of hollow slab beam and to judge its working status.Design/methodology/approachFirst, static load experiments are carried out on two in-service RC hollow slab beams; meanwhile, nonlinear finite element models are built to study the bearing capacity of them. The nonlinear material and shear slip effect of studs are considered in the models. Second, the finite element models are verified, and the numerical simulation results are in good agreement with the experimental results. Finally, the finite element models are adopted to carry out the research on the influence of different steel plate thicknesses on the flexural bearing capacity and ductility.FindingsWhen steel plates of different thicknesses are adopted to reinforce RC hollow slab beams, the bearing capacity increases with the increase of the steel plate thickness in a certain range. But when the steel plate thickness reaches a certain level, bearing capacity is no longer influenced. The displacement ductility coefficient decreases with the increase of steel plate thickness.Originality/valueBased on experimental study, this paper makes an extrapolation analysis of the bearing capacity of hollow slab beams reinforced with steel plates of different thicknesses through finite element simulation and discusses the influence on ductility. This method not only ensures the accuracy of bearing capacity evaluation but also does not need many samples, which is economical to a certain extent. The research results provide a basis for the reinforcement design of similar bridges.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Long Liu ◽  
Lifeng Wang ◽  
Ziwang Xiao

PurposeReinforcement of reinforced concrete (RC) beams in-service have always been an important research field, anchoring steel plate in the bottom of the beams is a kind of common reinforcement methods. In actual engineering, the contribution of pavement layer to the bearing capacity of RC beams is often ignored, which underestimates the bearing capacity and stiffness of RC beams to a certain extent. The purpose of this paper is to study the effect of pavement layer on the RC beams before and after reinforcement.Design/methodology/approachFirst, static load experiments are carried out on three in-service RC hollow slab beams, meanwhile, nonlinear finite element models are built to study the bearing capacity of them. The nonlinear material and shear slip effect of studs are considered in the models. Second, the finite element models are verified, and the numerical simulation results are in good agreement with the experimental results. Last, the finite element models are adopted to carry out the research on the influence of different steel plate thicknesses on the flexural bearing capacity and ductility.FindingsThe experimental results showed that pavement layers increase the flexural capacity of hollow slab beams by 16.7%, and contribute to increasing stiffness. Ductility ratio of SPRCB3 and PRCB2 was 30% and 24% lower than that of RCB1, respectively. The results showed that when the steel plate thickness was 1 mm–6 mm, the bearing capacity of the hollow slab beam increased gradually from 2158.0 kN.m to 2656.6 kN.m. As the steel plate thickness continuously increased to 8 mm, the ultimate bearing capacity increased to 2681.0 kN.m. The increased thickness did not cause difference to the bearing capacity, because of concrete crushing at the upper edge.Originality/valueIn this paper, based on the experimental study, the bearing capacity of hollow beam strengthened by steel plate with different thickness is extrapolated by finite element simulation, and its influence on ductility is discussed. This method not only guarantees the accuracy of the bearing capacity evaluation, but also does not require a large number of samples, and has certain economy. The research results provide a basis for the reinforcement design of similar bridges.


Author(s):  
João PO Freitas ◽  
Bruno Agostinho Hernandez ◽  
Paulo J Paupitz Gonçalves ◽  
Edmea C Baptista ◽  
Edson A Capello Sousa

Dental implants are widely used as a long-term treatment solution for missing teeth. A titanium implant is inserted into the jawbone, acting as a replacement for the lost tooth root and can then support a denture, crown or bridge. This allows discreet and high-quality aesthetic and functional improvement, boosting patient confidence. The use of implants also restores normal functions such as speech and mastication. Once an implant is placed, the surrounding bone will fuse to the titanium in a process known as osseointegration. The success of osseointegration is dependent on stress distribution within the surrounding bone and thus implant geometry plays an important role in it. Optimisation analyses are used to identify the geometry which results in the most favourable stress distribution, but the traditional methodology is inefficient, requiring analysis of numerous models and parameter combinations to identify the optimal solution. A proposed improvement to the traditional methodology includes the use of Design of Experiments (DOE) together with Response Surface Methodology (RSM). This would allow for a well-reasoned combination of parameters to be proposed. This study aims to use DOE, RSM and finite element models to develop a simplified optimisation analysis method for dental implant design. Drawing on data and results from previous studies, two-dimensional finite element models of a single Branemark implant, a multi-unit abutment, two prosthetic screws, a prosthetic crown and a region of mandibular bone were built. A small number of combinations of implant diameter and length were set based on the DOE method to analyse the influence of geometry on stress distribution at the bone-implant interface. The results agreed with previous studies and indicated that implant length is the critical parameter in reducing stress on cortical bone. The proposed method represents a more efficient analysis of multiple geometrical combinations with reduced time and computational cost, using fewer than a third of the models required by the traditional methods. Further work should include the application of this methodology to optimisation analyses using three-dimensional finite element models.


2006 ◽  
Vol 39 (13) ◽  
pp. 2457-2467 ◽  
Author(s):  
Fulvia Taddei ◽  
Luca Cristofolini ◽  
Saulo Martelli ◽  
H.S. Gill ◽  
Marco Viceconti

2004 ◽  
Vol 126 (1) ◽  
pp. 122-125 ◽  
Author(s):  
Xiang Wang, ◽  
Xiangyi Liu, and ◽  
Glen L. Niebur

The Orientation of trabecular bone specimens for mechanical testing must be carefully controlled. A method for accurately preparing on-axis cylindrical specimens using high-resolution micro-CT imaging was developed. Sixteen cylindrical specimens were prepared from eight bovine tibiae. High-resolution finite element models were generated from micro-CT images of parallelepipeds and used to determine the principal material coordinate system of each parallelepiped. A cylindrical specimen was then machined with a diamond coring bit. The resulting specimens were scanned again to evaluate the orientation. The average deviation between the principal fabric orientation and the longitudinal axis of the cylindrical specimen was only 4.70±3.11°.


Volume 2 ◽  
2004 ◽  
Author(s):  
Saiphon Charoenphan ◽  
Apiwon Polchai

The energy release rates in human cortical bone are investigated using a hybrid method of experimental and finite element modeling techniques. An explicit finite element analysis was implemented with an energy release rate calculation for evaluating this important fracture property of bones. Comparison of the critical value of the energy release rate, Gc, shows good agreement between the finite element models and analytical solutions. The Gc was found to be approximately 820–1150 J/m2 depending upon the samples. Specimen thickness appears to have little effect on the plane strain condition and pure mode I assumption. Therefore the energy release rate can be regarded as a material constant and geometry independent and can be determined with thinner specimens. In addition, the R curve resulting from the finite element models during slow crack growth shows slight ductility of the bone specimen that indicates an ability to resist crack propagation. Oscillations were found at the onset of the crack growth due to the nodal releasing application in the models. In this study light mass-proportional damping was used to suppress the noises. Although this techniques was found to be efficient for this slow crack growth simulation, other methods to continuously release nodes during the crack growth would be recommended for rapid crack propagation.


2011 ◽  
Vol 133 (4) ◽  
Author(s):  
Jonathan P. Vande Geest ◽  
B. R. Simon ◽  
Paul H. Rigby ◽  
Tyler P. Newberg

Finite element models (FEMs) including characteristic large deformations in highly nonlinear materials (hyperelasticity and coupled diffusive/convective transport of neutral mobile species) will allow quantitative study of in vivo tissues. Such FEMs will provide basic understanding of normal and pathological tissue responses and lead to optimization of local drug delivery strategies. We present a coupled porohyperelastic mass transport (PHEXPT) finite element approach developed using a commercially available ABAQUS finite element software. The PHEXPT transient simulations are based on sequential solution of the porohyperelastic (PHE) and mass transport (XPT) problems where an Eulerian PHE FEM is coupled to a Lagrangian XPT FEM using a custom-written FORTRAN program. The PHEXPT theoretical background is derived in the context of porous media transport theory and extended to ABAQUS finite element formulations. The essential assumptions needed in order to use ABAQUS are clearly identified in the derivation. Representative benchmark finite element simulations are provided along with analytical solutions (when appropriate). These simulations demonstrate the differences in transient and steady state responses including finite deformations, total stress, fluid pressure, relative fluid, and mobile species flux. A detailed description of important model considerations (e.g., material property functions and jump discontinuities at material interfaces) is also presented in the context of finite deformations. The ABAQUS-based PHEXPT approach enables the use of the available ABAQUS capabilities (interactive FEM mesh generation, finite element libraries, nonlinear material laws, pre- and postprocessing, etc.). PHEXPT FEMs can be used to simulate the transport of a relatively large neutral species (negligible osmotic fluid flux) in highly deformable hydrated soft tissues and tissue-engineered materials.


Sign in / Sign up

Export Citation Format

Share Document