The strain energy density ratio criterion for predicting cracking direction in composite materials

The strain energy density ratio criterion for predicting cracking direction incomposite materials is proposed.The Tsai-Hill criterion and Norris criterion ofcomposite materials are extended to predict the cracking direction in composites.Thethree criteria are used to analyse the crack propagation problem of the unidirectionalfibre composite sheet with various fibre directions.The predicted results are comparedwith those of the existing normal stress ratio criterion and strain energy densitycriterion.     

Stress and energy density fields near a blunted crack front in power hardening material

Using Jaumann and Dienes rates of Euler stress in elastic-plastic constitutive equations of finite deformation, plane strain finite element analysis for a compact tension specimen with a blunted crack front is made. The Euler stress, Kirchhoff stress and volume strain energy density near a blunted crack tip are computed. Constitutive relations with different deformation rates affect the the near crack tip solution in a region within an order of magnitude of the crack opening displacement. The results differed from the corresponding solution of deformation plasticity (or nonlinear elasticity) with increasing deformation. They are smaller in a local region of about 2 to 10 times of the crack opening distance.The volume energy density near the crack tip is computed, the stationary values of which determine the locations of extensive yielding and possible sites of crack initiation. It remained nearly constant with increasing deformation. Such a character tends to support the volume energy density criterion as a means for quantifying the ductile fracture behavior of metals.     

Elastic interaction of partially debonded circular inclusions. II. Application to fibrous composite

A complete analytical solution has been obtained of the elasticity problem for a plane containing periodically distributed, partially debonded circular inclusions, regarded as the representative unit cell model of fibrous composite with interface damage. The displacement solution is written in terms of periodic complex potentials and extends the approach recently developed by Kushch et al. (2010) to the cell type models. By analytical averaging the local strain and stress fields, the exact formulas for the effective transverse elastic moduli have been derived. A series of the test problems have been solved to check an accuracy and numerical efficiency of the method. An effect of interface crack density on the effective elastic moduli of periodic and random structure FRC with interface damage has been evaluated. The developed approach provides a detailed analysis of the progressive debonding phenomenon including the interface cracks cluster formation, overall stiffness reduction and damage-induced anisotropy of the effective elastic moduli of composite.     

陶瓷基复合材料基体唇形裂纹失效规律研究

建立横向拉伸载荷下的唇形裂纹数学模型,采用复变函数的方法,通过保角映射,推导了唇形裂纹尖端应力场和位移场的解析解,建立了唇形裂纹的应力强度因子准则和最大能量释放率准则,结合算例分析陶瓷基复合材料基体唇形裂纹的几何参数、外载荷和纤维分布对失效准则的影响规律.结果 表明,(1)裂纹尖端应力场和位移场的解析解与有限元计算结果进行对比,验证了方法的有效性;(2)相较于Griffith裂纹和椭圆裂纹,基于唇形裂纹的失效准则对裂纹尖端的敏感性更高,适用于预测脆性陶瓷基体裂纹的扩展;(3)对于不同几何参数的唇形裂纹,采用最大能量释放率准则的基体裂纹的扩展速率要大于应力强度因子准则.     

A rigid triangular inclusion partially bonded in an elastic matrix

The two-dimensional problem of a rigid rounded-off angle triangular inclusion partially bonded in an infinite elastic plate is studied. The unbonded part of the inclusion boundary forms an interfacial crack. Based on the complex variable method for curvilinear boundaries, the problem is reduced to a non-homogeneous Hilbert problem and the stress and displacement fields in the plate are obtained in closed form. Special attention is paid in the investigation of the stress field in the vicinity of the crack tip. It is found that the stresses present an oscillatory singularity and the general equations for the local stresses are derived. The singular stress field is coupled with the maximum circumferential stress and the minimum strain energy density criteria to study the fracture characteristics of the composite plate. Results are given for the complex stress intensity factors, the local stresses, the crack extension angles and the critical applied loads for unstable crack growth from its more vulnerable tip or two types of interfacial cracks along the inclusion boundary.     

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.     

Interaction of general plane P wave and cylindrical inclusion partially debonded from its viscoelastic matrix

The interaction of a general plane P wave and an elastic cylindrical inclusion of infinite length partially debonded from its surrounding viscoelastic matrix of infinite extension is investigated. The debonded region is modeled as an arc-shaped interface crack between inclusion and matrix with non-contacting faces. With wave functions expansion and singular integral equation technique, the interaction problem is reduced to a set of simultaneous singular integral equations of crack dislocation density function. By analysis of the fundamental solution of the singular integral equation, it is found that dynamic stress field at the crack tip is oscillatory singular, which is related to the frequency of incident wave. The singular integral equations are solved numerically, and the crack open displacement and dynamic stress intensity factor are evaluated for various incident angles and frequencies. The project supported by the National Natural Science Foundation of China (19872002) and Climbing Foundation of Northern Jiaotong University     

Multiple cracks in thermoelectroelastic bimaterials

The formulation for thermal stress and electric displacement in an infinite thermopiezoelectric plate with an interface and multiple cracks is presented. Using Green's function approach and the principle of superposition, a system of singular integral equations for the unknown temperature discontinuity defined on each crack face is developed and solved numerically. The formulation can then be used to calculate some fracture parameters such as the stress–electric displacement and strain energy density factor. The direction of crack growth for many cracks in thermopiezoelectric bimaterials is predicted by way of the strain energy density theory. Numerical results for stress–electric displacement factors and crack growth direction at a particular crack tip in two crack system of bimaterials are presented to illustrate the application of the proposed formulation.     

Effect of orthotropy on crack propagation

The tendency of moving cracks to spread along the preferred directions of material anisotropy is treated. Depending on the velocity of crack propagation, the change of material properties in orthogonal planes is shown to affect the bifurcation characteristics. The problem is reduced to a system of dual integral equations that can be solved in a standard fashion. Of particular interest is the dynamic stress field near the tip of a moving crack in an orthotropic material. Although the 1√r stress singularity is preserved with r being the radial distance measured from the crack tip, the angular variations of the stresses are dependent on crack speed and material anisotropy. The possibility of crack bifurcation is examined by application of the strain energy density criterion for several composite systems. Crack branching is found to be enhanced by material anisotropy, a phenomenon that is not uncommon in composite materials.     

Inclined semi-elliptical crack for predicting crack growth direction based on apparent stress intensity factors

Linear elastic criterion of the inclined semi-elliptical crack growth direction is elaborated on the basis of the strain energy density theory. Stress and displacement fields are presented for higher order terms asymptotic expansion. Solutions for elastic stress intensity factors are accounting for the function describing of the crack tip fields near the free surface of plate. The mixed mode behavior of crack growth direction angle along the semi-elliptical crack front for different combination of biaxial loading, inclination crack angle and surface flaw geometry is determined.     

Magneto-electro-elastic antiplane analysis of a bimaterial BaTiO3–CoFe2O4 composite wedge with an interface crack

The antiplane analysis is made for a bimaterial BaTiO₃–CoFe₂O₄ composite wedge containing an interface crack. The coupled magneto-electro-elastic field is induced by the piezoelectric/piezomagnetic BaTiO₃–CoFe₂O₄ composite materials. For the crack problems, the intensity factors of stress, strain, electric displacement, electric field, magnetic induction and magnetic field at crack tips are derived analytically. Also, the energy density criterion is applied to predict the fracture behavior of the interface crack. The numerical results also show that the energy release rate for a crack in a single wedge is negative.     

Theoretical model of crack branching in magnetoelectric thermoelastic materials

Thermomagnetoelectroelastic crack branching of magnetoelectro thermoelastic materials is theoretically investigated based on Stroh formalism and continuous distribution of dislocation approach. The crack face boundary condition is assumed to be fully thermally, electrically and magnetically impermeable. Explicit Green’s functions for the interaction of a crack and a thermomagnetoelectroelastic dislocation (i.e., a thermal dislocation, a mechanical dislocation, an electric dipole and a magnetic dipole located at a same point) are presented. The problem is reduced to two sets of coupled singular integral equations with the thermal dislocation and magnetoelectroelastic dislocation densities along the branched crack line as the unknown variables. As a result, the formulations for the stress, electric displacement and magnetic induction intensity factors and energy release rate at the branched crack tip are expressed in terms of the dislocation density functions and the branch angle. Numerical results are presented to study the effect of applied thermal flux, electric field and magnetic field on the crack propagation path by using the maximum energy release rate criterion.     

Micro-crack tip fracture of commercial grade aluminum under mixed mode loading

In situ tensile tests were made in a scanning electron microscope (SEM) to investigate the deformation and micro-fracture in the immediate vicinity of a micro-crack tip in commercial pure aluminum with large-size crystal. Examined are the slip line field, stress intensity factor, strain energy density factor and crack tip opening displacement (CTOD) for mixed mode loading. Blunting and sharpening effects are observed. The latter is controlled by localized slip while the former by uniformed slip of the operating slip system with the highest crack tip Schmid factor. The operating slip system depends on the crystallographic orientation of crystal containing micro-cracks.The damage and fracture take place in the blunted region and depend on the coarsening and spacing of uniformed slip lines. The mixed mode micro-crack propagates along the direction where the voids grow and coalesce into the micro-crack. The direction also depends on the orientation of the applied loading. This suggests that the formation of macro-fracture mechanics could be applied. In particular, the minimum strain energy density criterion is suitable for determining the direction of micro-crack instability in the mixed mode. The in situ data were used to yield a nearly constant critical, minimum strain energy density factor for onset of micro-cracking.     

Thermo-electro-structural coupled analyses of crack arrest by Joule heating

Using the finite element method, thermo-electro-structural coupled analyses of the cracked conducting plate under high electric current have been solved. The crack contact condition and temperature-dependent material properties are considered in this analysis. The crack tip temperature, electric current density factor, stress intensity factor and strain energy density factor are obtained for discussions. Due to high electric current density and Joule heating at the crack tip, a circular melting area may exist around the tip. After cooling, a circular void or hole may occur at the crack tip and the crack arrest is achieved. The crack tip temperature decreases when the crack contact area increases. The proper tensile load is necessary for making the crack open enough and causing high current density at the crack tip and associated crack arrest. On the other hand, the crack tip temperature increases with time by the increasing external current and Joule heating. The values of mode-I stress intensity factor and strain energy density factor decrease with time due to the thermal deformation around the crack tip. Because of the temperature-dependent resistivity, the variation of the electric current density factor is complicated. In addition, it is not easy to create a crack-arrest condition when the crack length relative to the plate width is too small.     

ASYMPTOTIC SOLUTIONS OF MODE I STEADY GROWTH CRACK IN MATERIALS UNDER CREEP CONDITIONS

An asymptotic analysis is made on problems with a steady-state crack growth coupled with a creep law model under tensile loads. Asymptotic equations of crack tip fields in creep materials are derived and solved numerically under small scale conditions. Stress and strain functions are adopted under a polar coordinate system. The governing equations of asymptotic fields are obtained by inserting the stress field and strain field into the material law. The crack growth rate rather than fracture criterion plays an important role in the crack tip fields of materials with creep behavior.     

Pitting formation under rolling contact

The mechanism of pitting caused by rolling contact is analyzed using the fracture mechanics approach. The governing factors are the initial crack length, crack angle, contact force, friction, strain hardened layer, and the hydraulic pressure of trapped fluid acting on the crack surface. Mode I and II stress intensity and the strain energy density factors are calculated by application of the two-dimensional finite element method. The strain energy density criterion is applied to show that shallow angle crack under small rolling contact force and friction enhances the probability of pitting under the roller’s running surface. The presence of a strain hardened surface layer also tends to affect the fracture behavior. The analytical results agree well with the experimental observations.     

铁中Ⅰ型裂纹裂尖场与裂尖形变机制分析

应用有限元及连续介质力学计算了平面应交条件下体心立方铁中不同取向裂纹的裂尖应力场和弹性应变能密度分布,从宏微观相结合角度分析了裂尖场分布与裂尖微结构演化的相互关联.指出裂尖塑性变形的具体形式与裂尖滑移面上分切应力的大小密切相关,并从能量角度解释了裂尖相变产生的原因,最后结合相关分子动力学研究成果探讨了裂尖奇异性区域.     

断裂力学判据的评述

从Inglis 和Griffith 的著名论文到Irwin 和Rice 等的奠基性贡献,对断裂力学中的线弹性断裂力学的K判据,界面断裂力学的G判据,和弹塑性断裂力学的J 判据作了扼要的综述. 介绍了在界面断裂力学G判据的基础上提出的界面断裂力学的K判据,以说明断裂力学的判据存在改进的可能性. 在综述中归纳出断裂力学判据中目前还没有较好解决的几个问题. 在总结以往断裂力学研究经验的基础上,指出裂纹端应力奇异性的源是对断裂力学判据存在的问题作进一步研究的切入点. 探讨了裂纹端应变间断的奇点是裂纹端应力奇异性的源的问题,从而对裂纹端应力强度因子的物理意义进行了讨论. 最后,阐述了进行可靠的裂纹端应力场的弹塑性分析是改进弹塑性断裂力学判据的关键,而进行可靠的裂纹端应力场的弹塑性分析的前提是要通过裂纹端应力奇异性的源的研究来获得作用在裂纹端的造成裂纹端应变间断的有限值应力.      

Effect of strain rate on crack growth in aluminum alloy 1100-0

Stress and damage analysis are performed to analyze the Mode I crack growth behavior of a central crack panel made of aluminum alloy 1100-0. On account of the highly nonhomogeneous stress state, each material element would experience a different strain rate depending on the location and loading rate. A data bank of uniaxial stress and strain curves is provided to cover the range local strain rates depending on the load time history. Such a approach is referred to as the strain rate dependent model in contrast to plasticity that utilizes a single constitutive relation.The strain energy density criterion is applied to determine the onset of crack initiation, stable crack growth and final termination. A unique feature of the approach is that the same criterion could describe the foregoing three distinct events of fracture behavior. Results are obtained for applied loads with different strain rates and compared with those obtained from the classical theory of plasticity, which is unconservative.