Finite Element Modeling for Energy Release Rate in Human Cortical Bone

Volume 2 ◽  
2004 ◽  
Author(s):  
Saiphon Charoenphan ◽  
Apiwon Polchai
Keyword(s):  
Finite Element ◽  
Crack Propagation ◽  
Finite Element Modeling ◽  
Energy Release Rate ◽  
Cortical Bone ◽  
Energy Release ◽  
Crack Growth ◽  
Release Rate ◽  
Slow Crack Growth ◽  
Finite Element Models

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.

The Characterization of Mixed-Mode Energy Release Rates in Orthotropic Materials

1996 ◽  
Vol 210 (1) ◽  
pp. 33-42 ◽  
Author(s):  
C A Walker ◽  
Jamasri
Keyword(s):  
Finite Element ◽  
Energy Release Rate ◽  
Energy Release ◽  
Crack Growth ◽  
Release Rate ◽  
Mode I ◽  
Mode Ii ◽  
Finite Element Models ◽  
Release Rates ◽  
Energy Release Rates

The aim of this work was to predict, from the material constants, mixed-mode energy release rates in orthotropic materials, in particular the general cases in which the crack is aligned at a random angle to the principal material direction, normal to the plane of orthotropy. Two-dimensional finite element models with various fibre orientations were generated. The finite element models were validated by comparing two sets of contour plots of deformation, one resulting from the finite element analysis and the other from moiré interferograms of the experimental work. On comparison there was shown to be a strict similarity between experimentally determined and computational deformation fields. Variations of the energy release rates were investigated for both rapid and stable crack growth. This was accomplished by generating two-dimensional stable crack growth finite element models. In general, energy release rates were found to be strongly affected by the fibre orientation. An increase of the angle of the crack growth direction caused a decrease of the mode I energy release rate and, by contrast, an increase of the mode II energy release rate, but the mode II energy release rate was always a small fraction of the mode I value. Crack extension caused a gradual increase of the mode I energy release rate both for coplanar and non-coplanar crack growth. However, there was no significant effect found on the mode II energy release rate.

Efficient Computation of Strain Energy Release Rate in Crack Growth Simulation

1999 ◽  
Author(s):  
D. J. Chen
Keyword(s):  
Finite Element ◽  
Crack Propagation ◽  
Energy Release Rate ◽  
Energy Release ◽  
Release Rate ◽  
Strain Energy ◽  
Strain Energy Release Rate ◽  
Strain Energy Release ◽  
Error Estimator ◽  
Efficient Computation

Abstract This paper utilizes an automated process to simplify the calculation of the strain energy release rate (SERR) during the crack propagation. The convergence of a finite element solution is achieved by adaptive re-meshing scheme with an error estimator of the linear strain triangular (LST) elements. As the desired mesh density is achieved, computation of the SERR using virtual crack closure technique (VCCT) can be obtained by using the static condensation scheme without re-analyzing the finite element models. Thus, the amount of computational and modeling time can be significantly reduced in the analysis of the crack propagation.

An Application of Finite Element-Based Fracture Mechanics Analysis to Cord-Rubber Structures

1996 ◽  
Vol 24 (3) ◽  
pp. 220-235 ◽  
Author(s):  
T. G. Ebbott
Keyword(s):  
Finite Element ◽  
Energy Release Rate ◽  
Energy Release ◽  
Crack Growth ◽  
Release Rate ◽  
Fatigue Behavior ◽  
Internal Crack ◽  
Materials Testing ◽  
Element Analysis ◽  
Rate Evaluation

Abstract A finite element-based method to analyze the severity of internal cracks in cord-rubber structures is presented. The method includes materials testing to characterize rubber fatigue behavior, a global-local finite element analysis to provide the detail necessary to model explicitly an internal crack, and use of the J-integral and virtual crack closure techniques for energy release rate evaluation. Analysis of the multiaxial and cyclic fracture situation is carried out by considering the cycle of each mode of fracture separately and then combining the effect of each mode to determine the total effect. Crack growth rates in the structure are assumed to be the same as the crack growth rate in a laboratory specimen at the same level of cyclic energy release rate. Results are presented for a material change in a critical tire region.

Crack Propagation Modeling in a PWR Under PTS Using XFEM

2020 ◽  
Author(s):  
Diego F. Mora ◽  
Markus Niffenegger ◽  
Roman Mukin
Keyword(s):  
Fracture Toughness ◽  
Finite Element ◽  
Crack Propagation ◽  
Energy Release Rate ◽  
Energy Release ◽  
Release Rate ◽  
Damage Evolution ◽  
Evolution Model ◽  
Critical Energy Release Rate ◽  
Damage Initiation

Abstract The finite element simulation of fracture propagation of BCC metals is challenging, as it needs to incorporate the brittle, ductile-brittle transition and ductile behavior presented by the fracture toughness. In this contribution, we restrict ourselves to the use of XFEM method to simulate the cleavage fracture due to initial flaws in the reactor pressure vessel of a reference design of the two-loop PWR nuclear power plant. A hypothetical large break loss of coolant accident is selected as accident scenario to obtain the loading conditions under which the crack is subjected. The thermal-mechanical calculation is performed using a finite element model of the whole RPV and the initial and boundary conditions are determined from the thermal-hydraulic simulation of the transient in TRACE. The method proposed in this contribution is based on the cohesive segment approach implemented in ABAQUS, which requires the definition of the damage properties of the material. The segment approach does not use the fracture toughness as failure criterion. Instead, it uses a traction separation law that is able to capture the brittle fracture behavior of ferritic steel. The crack propagation in XFEM uses a propagation criterion based on a cohesive damage initiation criterion and a damage evolution model. In order to implement the damage evolution model, the fracture energy release rate is directly related to the fracture toughness. The postulated crack is inserted in a submodel to reduce the computational cost of the calculation. The location of such submodel corresponds to the region of the core that suffers maximum irradiation and is subjected to high tensile stresses due to the cooling plume generated during the transient PTS cooling. The crack propagation analysis of postulated axial crack showed that its propagation happens in axial direction in those finite elements close to the inner surface because the energy release rate GI is larger than the critical energy release rate GIC. At the deepest point of the crack, the stresses in the finite element fulfil the damage initiation criterion but the crack does not propagate in radial direction (GI < GIC).

Calculating Energy Release Rate as a Function of Crack Length Using a Multiple-Step Crack Closure Technique in Tire Finite Element Models

2018 ◽  
Vol 46 (3) ◽  
pp. 130-152
Author(s):  
Dennis S. Kelliher
Keyword(s):  
Finite Element ◽  
Energy Release Rate ◽  
Crack Length ◽  
Energy Release ◽  
Crack Closure ◽  
Release Rate ◽  
Fatigue Behavior ◽  
Three Dimensions ◽  
Quantitative Prediction ◽  
Closure Technique

ABSTRACT When performing predictive durability analyses on tires using finite element methods, it is generally recognized that energy release rate (ERR) is the best measure by which to characterize the fatigue behavior of rubber. By addressing actual cracks in a simulation geometry, ERR provides a more appropriate durability criterion than the strain energy density (SED) of geometries without cracks. If determined as a function of crack length and loading history, and augmented with material crack growth properties, ERR allows for a quantitative prediction of fatigue life. Complications arise, however, from extra steps required to implement the calculation of ERR within the analysis process. This article presents an overview and some details of a method to perform such analyses. The method involves a preprocessing step that automates the creation of a ribbon crack within an axisymmetric-geometry finite element model at a predetermined location. After inflating and expanding to three dimensions to fully load the tire against a surface, full ribbon sections of the crack are then incrementally closed through multiple solution steps, finally achieving complete closure. A postprocessing step is developed to determine ERR as a function of crack length from this enforced crack closure technique. This includes an innovative approach to calculating ERR as the crack length approaches zero.

scholarly journals Effect of the Characteristic Size and Content of Graphene on the Crack Propagation Path of Alumina/Graphene Composite Ceramics

Materials ◽  
2021 ◽  
Vol 14 (3) ◽  
pp. 611
Author(s):  
Benshuai Chen ◽  
Guangchun Xiao ◽  
Mingdong Yi ◽  
Jingjie Zhang ◽  
Tingting Zhou ◽  
...  
Keyword(s):  
Crack Propagation ◽  
Energy Release Rate ◽  
Energy Release ◽  
Release Rate ◽  
Volume Content ◽  
Phase Interface ◽  
Average Energy ◽  
Characteristic Size ◽  
Composite Ceramics ◽  
Graphene Composite

In this paper, the Voronoimosaic model and the cohesive element method were used to simulate crack propagation in the microstructure of alumina/graphene composite ceramic tool materials. The effects of graphene characteristic size and volume content on the crack propagation behavior of microstructure model of alumina/graphene composite ceramics under different interfacial bonding strength were studied. When the phase interface is weak, the average energy release rate is the highest as the short diameter of graphene is 10–50 nm and the long diameter is 1600–2000 nm. When the phase interface is strong, the average energy release rate is the highest as the short diameter of graphene is 50–100 nm and the long diameter is 800–1200 nm. When the volume content of graphene is 0.50 vol.%, the average energy release rate reaches the maximum. When the velocity load is 0.005 m s−1, the simulation result is convergent. It is proven that the simulation results are in good agreement with the experimental phenomena.

Shape and energies of a dynamically propagating crack under bending

2003 ◽  
Vol 18 (10) ◽  
pp. 2379-2386 ◽  
Author(s):  
Dov Sherman ◽  
Ilan Be'ery
Keyword(s):  
Crack Propagation ◽  
Energy Release Rate ◽  
Energy Release ◽  
Release Rate ◽  
Crack Front ◽  
Three Dimensional ◽  
Strain Energy Release ◽  
Governing Parameter ◽  
Propagating Crack ◽  
Front Shape

We report on the exact shape of a propagating crack in a plate with a high width/thickness ratio and subjected to bending deformation. Fracture tests were carried out with brittle solids—single crystal, polycrystalline, and amorphous. The shape of the propagating crack was determined from direct temporal crack length measurements and from the surface perturbations generated during rapid crack propagation. The shape of the crack profile was shown to be quarter-elliptical with a straight, long tail; the governing parameter of the ellipse axes is the specimen's thickness at most length of crack propagation. Universality of the crack front shape is demonstrated. The continuum mechanics approach applicable to two-dimensional problems was used in this three-dimensional problem to calculate the quasistatic strain energy release rate of the propagating crack using the formulations of the dynamic energy release rate along the crack loci. Knowledge of the crack front shape in the current geometry and loading configuration is important for practical and scientific aspects.

Sign in / Sign up

Export Citation Format

Share Document