A Three-Dimensional Crack Tip Element for Energy Release Rate Determination in Layered Elastic Structures

2001 ◽  
Vol 35 (6) ◽  
pp. 457-488 ◽  
Author(s):  
Lijie Yu ◽  
Barry D. Davidson
2011 ◽  
Vol 488-489 ◽  
pp. 363-366
Author(s):  
Karsten Wippler ◽  
Meinhard Kuna

A general purpose direct BEM code has been developed for three-dimensional crack problems in piezoelectric structures. Special 3D non-continuous crack tip elements and several techniques for determining crack tip parameters were implemented. To calculate the electromechanical energy release rate for a virtual crack extension, the θ-method is employed, which was originally suggested by BONNET for linear elastic materials. The paper presents the generalization and numerical realization of the θ-method to 3D piezoelectric cracks. The great advantage of the θ-method is the direct computation of energy release rate, whereas the way via K-factors and the IRWIN matrix is more complicated. The efficiency and accuracy of the technique are shown for various example problems by comparing with analytical solutions.


1995 ◽  
Vol 62 (2) ◽  
pp. 294-305 ◽  
Author(s):  
B. D. Davidson ◽  
Hurang Hu ◽  
R. A. Schapery

A previously developed linear elastic crack-tip element analysis is reviewed briefly, and then extended and refined for practical applications. The element provides analytical expressions for total energy release rate and mode mix in terms of plate theory force and moment resultants near the crack tip. The element may be used for cracks within or between homogeneous isotropic or orthotropic layers, as well as for delamination of laminated composites. Classical plate theory is used to derive the equations for total energy release rate and mode mix; a “mode mix parameter,” Ω, as obtained from a separate continuum analysis is necessary to complete the mode mix decomposition. This parameter depends upon the elastic and geometrical properties of the materials above and below the crack plane, but not on the loading. A relatively simple finite element technique for determining the mode-mix parameter is presented and convergence in terms of mesh refinement is studied. Specific values of Ω are also presented for a large number of cases. For those interfaces where a linear elastic solution predicts an oscillatory singularity, an approach is described which allows a unique, physically meaningful value of fracture mode ratio to be defined. This approach is shown to provide predictions of crack growth between dissimilar homogeneous materials that are equivalent to those obtained from the oscillatory field solution. Application of the approach to delamination in fiber-reinforced laminated composites is also discussed.


2003 ◽  
Vol 18 (10) ◽  
pp. 2379-2386 ◽  
Author(s):  
Dov Sherman ◽  
Ilan Be'ery

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.


2020 ◽  
Vol 10 (12) ◽  
pp. 4227
Author(s):  
Yali Yang ◽  
Seok Jae Chu ◽  
Wei song Huang ◽  
Hao Chen

The evaluation of energy release rate with angle is still a challenging task in metal crack propagation analysis, especially for the mixed Mode I-II-III loading situation. In this paper, the energy release rate associated with stress intensity factors at an arbitrary angle under mixed mode loadings has been investigated using both a numerical method and theoretical derivation. A relatively simple and precise numerical method was established through a series of spatial-inclined ellipses in Mode I-II and ellipsoids in Mode I-II-III, with different propagation angles computed from simulation. Meanwhile, a theoretical expression of the energy release rate with angle for a crack tip under a I-II-III mixed mode crack was deduced based on the propagation mechanism of the crack tip under the influence of a stress field. It is confirmed that the theoretical expression deduced could provide results as accurately as the present numerical method. The present results were confirmed to be effective and accurate by comparison with experimental data and other literature.


2014 ◽  
Vol 81 (10) ◽  
Author(s):  
Theocharis Baxevanis ◽  
Chad M. Landis ◽  
Dimitris C. Lagoudas

A finite element analysis of steady-state crack growth in pseudoelastic shape memory alloys under the assumption of adiabatic conditions is carried out for plane strain, mode I loading. The crack is assumed to propagate at a critical level of the crack-tip energy release rate and the fracture toughness is obtained as the ratio of the far-field applied energy release rate to the crack-tip critical value. Results related to the influence of latent heat on the near-tip stress field and fracture toughness are presented for a range of parameters related to thermomechanical coupling. The levels of fracture toughness enhancement, associated with the energy dissipated by the transformed material in the wake of the growing crack, are found to be lower under adiabatic conditions than under isothermal conditions [Baxevanis et al., 2014, J. Appl. Mech., 81, 041005]. Given that in real applications of shape memory alloy (SMA) components the processes are usually not adiabatic, which is the case with the lowest energy dissipation during a cyclic loading–unloading process (hysteresis), it is expected that the actual level of transformation toughening would be higher than the one corresponding to the adiabatic case.


Sign in / Sign up

Export Citation Format

Share Document