Modeling of Dynamic Fracture and Damage in Two-Dimensional Trabecular Bone Microstructures Using the Cohesive Finite Element Method

2008 ◽  
Vol 130 (2) ◽  
Author(s):  
Vikas Tomar

Trabecular bone fracture is closely related to the trabecular architecture, microdamage accumulation, and bone tissue properties. Micro-finite-element models have been used to investigate the elastic and yield properties of trabecular bone but have only seen limited application in modeling the microstructure dependent fracture of trabecular bone. In this research, dynamic fracture in two-dimensional (2D) micrographs of ovine (sheep) trabecular bone is modeled using the cohesive finite element method. For this purpose, the bone tissue is modeled as an orthotropic material with the cohesive parameters calculated from the experimental fracture properties of the human cortical bone. Crack propagation analyses are carried out in two different 2D orthogonal sections cut from a three-dimensional 8mm diameter cylindrical trabecular bone sample. The two sections differ in microstructural features such as area fraction (ratio of the 2D space occupied by bone tissue to the total 2D space), mean trabecula thickness, and connectivity. Analyses focus on understanding the effect of the rate of loading as well as on how the rate variation interacts with the microstructural features to cause anisotropy in microdamage accumulation and in the fracture resistance. Results are analyzed in terms of the dependence of fracture energy dissipation on the microstructural features as well as in terms of the changes in damage and stresses associated with the bone architecture variation. Besides the obvious dependence of the fracture behavior on the rate of loading, it is found that the microstructure strongly influences the fracture properties. The orthogonal section with lesser area fraction, low connectivity, and higher mean trabecula thickness is more resistant to fracture than the section with high area fraction, high connectivity, and lower mean trabecula thickness. In addition, it is found that the trabecular architecture leads to inhomogeneous distribution of damage, irrespective of the symmetry in the applied loading with the fracture of the entire bone section rapidly progressing to bone fragmentation once the accumulated damage in any trabeculae reaches a critical limit.

Author(s):  
Vikas Tomar

Trabecular bone fracture is closely related to the trabecular architecture and microdamage accumulation. Micro-finite element models have been used to investigate the elastic and yield properties of trabecular bone but have only seen limited application in modeling the microstructure dependent fracture of trabecular bone, [1, 2]. In the presented research a cohesive finite element method (CFEM) based approach that can be used to model microstructure and loading rate dependent fracture in trabecular bone is developed for the first time. The emphasis is on understanding the effect of the rate of loading and its correlation with the bone microstructure on the microdamage accumulation and fracture behavior in the trabecular bone. Analyses focus on understanding the effect of the rate of loading, change in bone tissue properties with aging, and their correlation with the bone microstructure on the microdamage accumulation and the fracture behavior in the trabecular bone.


Author(s):  
G. Bobik ◽  
J. Żmudzki ◽  
K. Majewska

Purpose: Difference in the mechanical properties of bone and stiffer femoral implant causes bone tissue resorption, which may result in implant loosening and periprosthetic fractures. The introduction of porous material reduces the stiffness of the implant. The aim of the study was to analyse the influence of porous shell of femoral revision implant on bone tissue loading distribution with use the finite element method. Design/methodology/approach: Load transfer in the femur has been investigated using the finite element method (Ansys). Cementless implant models were placed in the anatomical femur model. Femur model included sponge bone and cortical bone. The solid implant was compared with the implant containing porous material in 70% in outer layer with a thickness of 2 mm. Load of 1500 N during gait was simulated. In addition, the forces of the ilio-tibial band and the abductor muscles were implemented, as well as the torque acting on the implant. Findings: Increase of stress for the porous model was found. The underload zones in bone have been reduced. Loading distribution was slightly more favourable, albeit rather in cortical bone. Stress value in cancellous bone around cementless implant margin has increased to a level that is dangerous for bone loss. Stress in the implant was not dangerous for damage. The stress distribution was different in the implant neck zone where the porous shell borne a little less load and high stress was shifted to the stiffer core. Research limitations/implications: Variable conditions for fitting the stem to the bone as well as the friction conditions were not investigated. Practical implications: Stress values in the spongy bone around the insertion edge of the cementless implant are consistent with long-term clinical results of the bone atrophy in 1 and 2 Gruen`s zones around the fully porous implants. Originality/value: The advantage of fully porous coated implant was the decrease of risk of trabecular bone tissue resorption around the implant tip and the increase of risk of trabecular bone tissue resorption around insertion edge of the implant.


2012 ◽  
Vol 135 (1) ◽  
Author(s):  
M. R. Hardisty ◽  
R. Zauel ◽  
S. M. Stover ◽  
D. P. Fyhrie

As the average age of the population has increased, the incidence of age-related bone fracture has also increased. While some of the increase of fracture incidence with age is related to loss of bone mass, a significant part of the risk is unexplained and may be caused by changes in intrinsic material properties of the hard tissue. This investigation focused on understanding how changes to the intrinsic damage properties affect bone fragility. We hypothesized that the intrinsic (μm) damage properties of bone tissue strongly and nonlinearly affect mechanical behavior at the apparent (whole tissue, cm) level. The importance of intrinsic properties on the apparent level behavior of trabecular bone tissue was investigated using voxel based finite element analysis. Trabecular bone cores from human T12 vertebrae were scanned using microcomputed tomography (μCT) and the images used to build nonlinear finite element models. Isotropic and initially homogenous material properties were used for all elements. The elastic modulus (Ei) of individual elements was reduced with a secant damage rule relating only principal tensile tissue strain to modulus damage. Apparent level resistance to fracture as a function of changes in the intrinsic damage properties was measured using the mechanical energy to failure per unit volume (apparent toughness modulus, Wa) and the apparent yield strength (σay, calculated using the 0.2% offset). Intrinsic damage properties had a profound nonlinear effect on the apparent tissue level mechanical response. Intrinsic level failure occurs prior to apparent yield strength (σay). Apparent yield strength (σay) and toughness vary strongly (1200% and 400%, respectively) with relatively small changes in the intrinsic damage behavior. The range of apparent maximum stresses predicted by the models was consistent with those measured experimentally for these trabecular bone cores from the experimental axial compressive loading (experimental: σmax = 3.0–4.3 MPa; modeling: σmax = 2–16 MPa). This finding differs significantly from previous studies based on nondamaging intrinsic material models. Further observations were that this intrinsic damage model reproduced important experimental apparent level behaviors including softening after peak load, microdamage accumulation before apparent yield (0.2% offset), unload softening, and sensitivity of the apparent level mechanical properties to variability of the intrinsic properties.


2012 ◽  
Vol 245 ◽  
pp. 57-62 ◽  
Author(s):  
Petr Marcián ◽  
Libor Borák ◽  
Ondřej Konečný ◽  
Petr Navrátil ◽  
Zdeněk Florian

This paper is focused on computational modeling of an interaction of dental implant with mandible bone. It describes creation of computational model including model of geometry, materials, loads and constraints. There is a comparative stress-strain analysis of the levels of cancellous bone model. Computations are performed with the use of finite element method. Results show differences between the model which includes trabecular architecture of cancellous bone tissue and the model with non-trabecular cancellous bone tissue. For better description of the processes in bone tissue and at the interface between bone tissue and implant, it is necessary to create the computational model on the highest possible level, i.e. with the trabecular bone tissue.


Author(s):  
Neander Berto Mendes ◽  
Lineu José Pedroso ◽  
Paulo Marcelo Vieira Ribeiro

ABSTRACT: This work presents the dynamic response of a lock subjected to the horizontal S0E component of the El Centro earthquake for empty and completely filled water chamber cases, by coupled fluid-structure analysis. Initially, the lock was studied by approximation, considering it similar to the case of a double piston coupled to a two-dimensional acoustic cavity (tank), representing a simplified analytical model of the fluid-structure problem. This analytical formulation can be compared with numerical results, in order to qualify the responses of the ultimate problem to be investigated. In all the analyses performed, modeling and numerical simulations were done using the finite element method (FEM), supported by the commercial software ANSYS.


Sign in / Sign up

Export Citation Format

Share Document