Influence of heat transfer and fluid flow on crack growth in multilayered porous/dense materials using XFEM: Application to Solid Oxide Fuel Cell like material design

An advanced numerical model is developed to investigate the influence of heat transfer and fluid flow on crack propagation in multi-layered porous materials. The fluid flow, governed by the Navier–Stokes and Darcy’s law, is discretized with the nonconforming Crouzeix–Raviart (CR) finite element method. A combination of Discontinuous Galerkin (DG) and Multi-Point Flux Approximation (MPFA) methods is used to solve the advection–diffusion heat transfer equation in the flow channel and in the fluid phase within the porous material. The crack is assumed to affect only the heat diffusion within the porous layer, therefore a time splitting technique is used to solve the heat transfer in the fluid and the solid phases separately. Thus, within the porous material, the crack induces a discontinuity of the temperature at the crack surfaces and a singularity of the flux at the crack tip. Conduction in the solid phase is solved using the eXtended Finite Element Method (XFEM) to better handle the discontinuities and singularities caused by the cracks. The XFEM is also used to solve the thermo-mechanical problem and to track the crack propagation. The multi-physics model is implemented then validated for the transient regime, this necessitated a post processing treatment in which, the stress intensity factors (SIF) are computed for each time step. The SIFs are then used in the crack propagation criterion and the crack orientation angle. The methodology seems to be robust accurate and the computational cost is reduced thanks to the XFEM.     

Asymptotic mechanical fields at the tip of a mode I crack in rubber-like solids

Asymptotic analyses of the mechanical fields in front of stationary and propagating cracks facilitate the understanding of the mechanical and physical state in front of crack tips, and they enable prediction of crack growth and failure. Furthermore, efficient modelling of arbitrary crack growth by use of XFEM (extended finite element method) requires accurate knowledge of the asymptotic crack tip fields. In the present work, we perform an asymptotic analysis of the mechanical fields in the vicinity of a propagating mode I crack in rubber. Plane deformation is assumed, and the material model is based on the Langevin function, which accounts for the finite extensibility of polymer chains. The Langevin function is approximated by a polynomial, and only the term of the highest order contributes to the asymptotic solution. The crack is predicted to adopt a wedge-like shape, i.e. the crack faces will be straight lines. The angle of the wedge and the order of the stress singularity depend on the hardening of the strain energy function. The present analysis shows that in materials with a significant hardening, the inertia term in the equations of motion becomes negligible in the asymptotic analysis. Hence, there is no upper theoretical limit to the crack speed.     

Finite strain analysis of crack tip fields in incompressible hyperelastic solids loaded in plane stress

A new method is developed to determine the dominant asymptotic stress and deformation fields near the tip of a Mode-I traction free plane stress crack. The analysis is based on the fully nonlinear equilibrium theory of incompressible hyperelastic solids. We show that the dominant singularity of the near tip stress field is governed by the asymptotic solution of a linear second order ordinary differential equation. Our method is applicable to any hyperelastic material with a smooth work function that depends only on the trace of the Cauchy-Green tensor and is particularly useful for materials that exhibit severe strain hardening. We apply this method to study two types of soft materials: generalized neo-Hookean solids and a solid that hardens exponentially. For the generalized neo-Hookean solids, our method is able to resolve a difficulty in the previous work by Geubelle and Knauss (1994a). Our theoretical results are compared with finite element simulations.     

Discussion on electromagnetic crack face boundary conditions for the fracture mechanics of magneto-electro-elastic materials

This paper discusses electromagnetic boundary conditions on crack faces in magneto- electroelastic materials, where piezoelectric, piezomagnetic and magnetoelectric effects are coupled. A notch of finite thickness in these materials is also addressed. Four idealized electromagnetic boundary conditions assumed for the crack-faces are separately investigated, i.e. (a) electrically and magnetically impermeable (crack-face), (b) electrically impermeable and magnetically permeable, (c) electrically permeable and magnetically impermeable, and (d) electrically and magnetically permeable. The influence of the notch thickness on important parameters, such as the field intensity factors, the energy release rate at the notch tips and the electromagnetic fields inside the notch, are studied and the results are obtained in closed-form. Results under different idealized electromagnetic boundary conditions on the crack-face are compared, and the applicability of these idealized assumptions is discussed.The project supported by the National Natural Science Foundation of China (10102004) The English text was polished by Yunming Chen.     

Non-local theory solution for in-plane shear of through crack

A non-local theory of elasticity is applied to obtain the plane strain stress and displacement field for a through crack under in-plane shear by using Schmidt's method. Unlike the classical elasticity solution, a lattice parameter enters into the problem that make the stresses finite at crack tip. Both the angular variations of the circumferential stress and strain energy density function are examined to associate their stationary value with locations of possible fracture initiation. The former criterion predicted a crack initiation angle of 54° from the plane of shear for the non-local solution as compared with about 75° for the classical elasticity solution. The latter criterion based on energy density yields a crack initiation angle of 80° for a Poisson's ratio of 0.28. This is much closer to the value that is predicted by the classical crack tips solution of elasticity.     

Experimental validation of the tip asymptotics for a fluid-driven crack

This paper provides experimental confirmation of the opening asymptotes that have been predicted to develop at the tip of fluid-driven cracks propagating in impermeable brittle elastic media. During propagation of such cracks, energy is dissipated not only by breaking of material bonds ahead of the tip but also by flow of viscous fluid. Theoretical analysis based on linear elastic fracture mechanics and lubrication theory predicts a complex multiscale asymptotic behavior of the opening in the tip region, which simplifies either as or as power law of the distance from the tip depending on whether the dominant mechanism of energy dissipation is bond breaking or viscous flow. The laboratory experiments entail the propagation of penny-shaped cracks by injection of glycerin or glucose based solutions in polymethyl methacrylate (PMMA) and glass specimens subjected to confining stresses. The full-field opening is measured from analysis of the loss of intensity as light passes through the dye-laden fluid that fills the crack. The experimental near-tip opening gives excellent agreement with theory and therefore confirms the predicted multi-scale tip asymptotics.     

Definition of a new predictor for multiaxial fatigue crack nucleation in rubber

From an engineering point of view, prediction of fatigue crack nucleation in automotive rubber parts is an essential prerequisite for the design of new components. We have derived a new predictor for fatigue crack nucleation in rubber. It is motivated by microscopic mechanisms induced by fatigue and developed in the framework of Configurational Mechanics. As the occurrence of macroscopic fatigue cracks is the consequence of the growth of pre-existing microscopic defects, the energy release rate of these flaws need to be quantified. It is shown that this microstructural evolution is governed by the smallest eigenvalue of the configurational (Eshelby) stress tensor. Indeed, this quantity appears to be a relevant multiaxial fatigue predictor under proportional loading conditions. Then, its generalization to non-proportional multiaxial fatigue problems is derived. Results show that the present predictor, which is related to the previously published predictors, is capable to unify multiaxial fatigue data.     

Impact response of a piezoelectric layered composite plate with a crack

The dynamic response of a central crack in a piezoelectric layered composite plate under normal impact is analyzed. The crack is oriented normally to the interfaces. The Laplace and Fourier transform techniques are used to formulate the problem in terms of a singular integral equation. The order of stress singularity around the tip of the terminated crack is also obtained. Numerical calculations are carried out, and the main results presented are the variations of the dynamic stress intensity factor and the dynamic energy density factor versus time as functions of the geometric parameters and the piezoelectric material properties of the layered composite plate.     

Analysis of a crack in a half-plane piezoelectric solid with traction-induction free boundary

In this paper, the problem of a crack embedded in a half-plane piezoelectric solid with traction-induction free boundary is analyzed. A system of singular integral equations is formulated for the materials with general anisotropic piezoelectric properties and for the crack with arbitrary orientation. The kernel functions developed are in complex form for general anisotropic piezoelectric materials and are then specialized to the case of transversely isotropic piezoelectric materials which are in real form. The obtained coupled mechanical and electric real kernel functions may be reduced to those kernel functions for purely elastic problems when the electric effects disappear. The system of singular integral equations is solved numerically and the coupling effects of the mechanical and electric phenomena are presented by the generalized stress intensity factors for transversely isotropic piezoelectric materials.     

A domain independent integral expression for the crack extension force of a curved crack in three dimensions

An integral expression that is domain independent in curvilinear coordinates and compatible with zero divergence of Eshelby's (Phil. Trans. Roy. Soc. (London) 244 (1951) 87.) energy momentum tensor was obtained from the principle of virtual work. By applying Eshelby's definition of the force on a material defect a general expression of the crack extension force for a curved crack in three dimensions, here called the F-integral, was derived from the domain independent integral expression. The F-integral is given explicitly for a number of curved cracks and found to be in agreement with previously known solutions, when available. The influence of crack surface and crack front curvature upon the various forms of the F-integral is discussed. The F-integral presented in this work is a generalisation of the J-integral (Rice, J. Appl. Mech. 35 (1968) 379.) to curved cracks in orthogonal curvilinear coordinates.     

Asymptotic Analysis of an Elasticity Equation for a Finger-Like Hydraulic Fracture

We derive a novel integral equation relating the fluid pressure in a finger-like hydraulic fracture to the fracture width. By means of an asymptotic analysis in the small height to length ratio limit we are able to establish the action of the integral operator for receiving points that lie within three distinct regions: (1) an outer expansion region in which the dimensionless pressure is shown to be equal to the dimensionless width plus a small correction term that involves the second derivative of the width, which accounts for the nonlocal effects of the integral operator. The leading order term in this expansion is the classic local elasticity equation in the PKN model that is widely used in the oil and gas industry; (2) an inner expansion region close to the fracture tip within which the action of the elastic integral operator is shown to be the same as that of a finite Hilbert transform associated with a state of plane strain. This result will enable pressure singularities and stress intensity factors to be incorporated into analytic models of these finger-like fractures in order to model the effect of material toughness; (3) an intermediate region within which the action of the Fredholm integral operator of the first kind is reduced to a second kind operator in which the integral term appears as a small perturbation which is associated with a convergent Neumann series. These results are important for deriving analytic models of finger-like hydraulic fractures that are consistent with linear elastic fracture mechanics. Submitted to Journal of Elasticity on February 5, 2007. Re-submitted with revisions on May 30, 2007.     

A weighted residual method for elastic-plastic analysis near a crack tip and the calculation of the plastic stress intensity factors

In this paper,a weighted residual method for the elastic-plastic analysis near a crack tip is systematically given by taking the model of power-law hardening under plane strain condition as a sample.The elastic-plastic solutions of the crack tip field and an approach based on the superposition of the nonlinear finite element method on the complete solution in the whole crack body field,to calculate the plastic stress intensity factors,are also developed.Therefore,a complete analysis based on the calculation both for the crack tip field and for the whole crack body field is provided.     

Analysis to the stress construction at the vicinity of crack tip for mode I fracture in 3-D state using finite element method

The stresses and strains are calculated for CT specimen of power hardening material in 3-D deformation state using ADINA finite element program, and the stress distribution at the vicinity of crack tip for Mode I fracture is analysed according to the results of calculation. It is found that the expression of stress can be written as the form of separation of variables of r and θ, then the function of r can be expanded in Laurant series. It is still found that the three normal stresses have the same order of magnitude. The conclusions offer two suppositions to obtain the theoretical expression of stress at the vinicity of crack tip for Mode I fracture with good ground, and the procedure of solution will be greatly simplified. This Work is Supported by the National Science Found.     

Observation of Stress Field Around an Oscillating Crack Tip in a Quenched Thin Glass Plate

The variation of stress field around an oscillating crack tip in a quenched thin glass plate is observed using instantaneous phase-stepping photoelasticity. The successive images around the propagating crack are recorded by a CCD camera that is equipped with a pixelated micro-retarder array. Then, the phase maps of the principal stress difference and the principal direction are easily obtained even though the photoelastic fringes cannot be visualized. The path of the crack growth as well as the stress intensity factors and the crack tip constraint are obtained from these phase distributions. Results show that the mode I stress intensity factor and the crack tip constraint vary remarkably with the crack growth. In addition, the results show that the mode-II stress intensity factor exists even though the crack propagates smoothly.