Growth of a penny-shaped crack with a nonsmall fracture process zone in a composite

The paper studies the stress rupture behavior of a reinforced viscoelastic composite through which a penny-shaped mode I crack propagates under a constant load. The composite has hexagonal symmetry and consists of elastic isotropic fibers and viscoelastic isotropic matrix. The material is modeled as a transversely isotropic homogeneous viscoelastic medium with effective characteristics. The crack is in the isotropy plane. The ring-shaped fracture process zone at the crack front is modeled by a modified Dugdale zone with time-dependent stresses. The viscoelastic properties of the matrix are characterized using a resolvent integral operator. Use is made of Volterra's principle, the method of operator continued fractions, and the theory of precritical crack growth in viscoelastic bodies. The problem is reduced to nonlinear integral equations. Numerical results are obtained for certain components of the composite, constant volume fractions, and different fracture strengths Translated from Prikladnaya Mekhanika, Vol. 44, No. 8, pp. 45–51, August 2008.     

基于结构化变形驱动的非局部宏-微观损伤模型的真Ⅱ型裂纹模拟

Ⅱ型载荷作用下裂纹变形模式也为Ⅱ型的破坏问题称为真Ⅱ型破坏.准确定量地把握真Ⅱ型破坏的全过程是具有挑战性的问题.本文采用结构化变形驱动的非局部宏-微观损伤模型对真Ⅱ型破坏问题进行了模拟.根据结构化变形理论将点偶的非局部应变分解为弹性应变与结构化应变两部分,进而利用Cauchy-Born准则与结构化应变计算点偶的结构化正伸长量.在本文中,结构化应变取为非局部应变的偏量部分.当点偶的结构化正伸长量超过临界伸长量时,微细观损伤开始在点偶层次发展.将微细观损伤在作用域中进行加权求和得到拓扑损伤,并通过能量退化函数将其嵌入到连续介质-损伤力学框架中进行数值求解.进一步地,本文采用Gauss-Lobatto积分格式计算点偶的非局部应变,将积分点数目降低到4个,显著降低了前处理和非线性分析的计算成本.通过对Ⅱ型加载下裂尖应变场的分析揭示了采用偏应变作为结构化应变的原因.对两个典型真Ⅱ型破坏问题的模拟结果表明,本文方法不仅可以把握Ⅱ型加载下的真Ⅱ型裂纹扩展模式,同时可以定量刻画加载过程中的载荷-变形曲线,且不具有网格敏感性.最后指出了需要进一步研究的问题.     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.     

层状磁电复合材料界面共线裂纹平面问题分析

本文研究了面内电磁势载荷作用下双层压电压磁复合材料中共线界面裂纹问题．考虑了压电材料的导磁性质和压磁材料的介电性质,引入了界面电位移和磁感强度的连续性条件．利用Fourier 变换得到一组第二类Cauchy 型奇异积分方程．进一步导出了相应问题的应力强度因子、电位移强度因子和磁感强度强度因子的表达式,给出了应力强度因子的数值结果．结果表明电磁载荷会导致界面裂纹尖端I、II 混合型应力奇异性,同时还伴随着电位移和磁感强度的奇异性．比较了双裂纹左右端的应力强度因子,发现在面内极化方向上施加面内磁势载荷时共线裂纹内侧尖端区域的两个法向应力场发生互相干涉增强．     

Transient response of an insulating crack between dissimilar piezoelectric layers under mechanical and electrical impacts

Summary  The dynamic response of an interface crack between two dissimilar piezoelectric layers subjected to mechanical and electrical impacts is investigated under the boundary condition of electrical insulation on the crack surface by using the integral transform and the Cauchy singular integral equation methods. The dynamic stress intensity factors, the dynamic electrical displacement intensity factor, and the dynamic energy release rate (DERR) are determined. The numerical calculation of the mode-I plane problem indicates that the DERR is more liable to be the token of the crack growth when an electrical load is applied. The dynamic response shows a significant dependence on the loading mode, the material combination parameters as well as the crack configuration. Under a given loading mode and a specified crack configuration, the DERR of an interface crack between piezoelectric media may be decreased or increased by adjusting the material combination parameters. It is also found that the intrinsic mechanical-electrical coupling plays a more significant role in the dynamic fracture response of in-plane problems than that in anti-plane problems. Received 4 September 2001; accepted for publication 23 July 2002 The work was supported by the National Natural Science Foundation under Grant Number 19891180, the Fundamental Research Foundation of Tsinghua University, and the Education Ministry of China.     

功能梯度板条断裂分析

现存文献关于功能梯度材料断裂问题的研究大都假设材料性质为坐标的指数函数或幂函数,而对其它函数形式较少采用。本文假设功能梯度材料剪切模量为坐标的双曲函数,而泊松比为常量,研究功能梯度板条的混合型裂纹问题。利用Fourier积分变换技术将混合边值问题转化为一对奇异积分方程,通过数值求解奇异积分方程获得含裂纹功能梯度板条分别在剪切和法向载荷作用下的I型和Ⅱ型应力强度因子,并讨论了材料的非均匀性和裂纹相对尺寸对裂纹尖端应力强度因子的影响。     

Crack energy density in arbitrary direction for piezoelectrics and its characteristic features in electromechanical mixed mode crack

The concepts of crack energy density (CED) and its derivatives in arbitrary direction were established for piezoelectric material and, keeping their application to mixed mode fracture in mind, the characteristic features of them as fracture parameters were investigated based on the approximate equations for CED and its derivatives. That is, CED and its derivatives in arbitrary direction are defined first and separation into their each mode contribution is made. Subsequently, path independent integral expressions of them are derived, and then using them, approximate equations of each mode contribution of CED are obtained concretely for the case where linear singular solution is known. The resulting equations are then used to investigate the effects of electric field and electrical boundary condition on CED and its derivatives. An infinite piezoelectric plane with a crack inclined with respect to the poling direction is considered as a numerical example. Mode I contribution of mechanical CED is mainly employed as a possible fracture parameter for the study and it was shown that applied electric field significantly influences on fracture parameters especially for the impermeable crack perpendicular to the poling direction. The effect of electric field has the tendency to decrease as crack inclination angle increases. It was also found that, even for the impermeable crack perpendicular to the poling direction, crack propagation could be deviated from self-similar direction under a strong negative electric field, and this fact is qualitatively consistent with an existing experimental observation. For the ideally sharp crack with no width, impermeable and Hao and Shen type boundary conditions are admissible showing qualitative agreement with experimental results, but exact boundary condition is not suitable and finally consistent with permeable boundary condition.     

Anti-plane analysis on a finite crack in a one-dimensional hexagonal quasicrystal strip

The anti-plane fracture problem for a finite crack in a one-dimensional hexagonal quasicrystal strip is analyzed. By using Fourier transforms, the mixed boundary value problems are reduced to the dual integral equations. The solution of the dual integral equations is then expressed by the complete elliptic integrals of the first and the third kinds. The expressions for stress, strains, displacements and field intensity factors of the phonon and phason fields near the crack tip are obtained exactly. The path-independent integral derived by a conservation law equals the energy release rate, which can be used as a fracture criterion for a mode III fracture problem.     

Mode-3 spontaneous crack propagation in unsymmetric functionally graded materials

The effects of the inhomogeneous material property variation on the mode-3 crack propagation characteristics have been analyzed. The spectral form of the elastodynamic boundary integral equations are derived for a functionally graded material, which is used as the numerical tool in our analysis. The material property gradient is assumed to vary unsymmetrically in the direction normal to the fracture plane. A parametric study has been performed by systematically varying the inhomogeneity length scale. In comparison with a homogeneous material, an unsymmetric functionally graded material offers lesser fracture resistance. The fracture resistance progressively decreases with increase in inhomogeneity, quantified by the increase in crack sliding displacement jumps, crack tip velocities and accelerations. The material property inhomogeneity affects only the transient crack propagation velocities, while the quasi-steady-state velocity remains unaltered.     

Modeling method for a crack problem of functionally graded materials with arbitrary properties—piecewise-exponential model

A new model, piecewise-exponential model (PE model), is developed to investigate the crack problem of the functionally graded materials (FGMs) with arbitrary properties. In the PE model, the functionally graded material is divided into some nonhomogeneous layers along the gradient direction of the properties, with each layer’s properties varying exponentially. By this way, the real material properties can be approached by a series of exponential functions. Since the real material properties are used on both surfaces of each nonhomogeneous layer, the nature of continuously varying properties of FGMs can be approached accurately. The influences of the local nonhomogeneity on the crack-tip fields can be fully considered. By using the new model, the fracture problem of a functionally graded strip with arbitrary properties and a crack vertical to the free surfaces is studied. The integral transform method, the theory of residues and the theory of singular integral equation are applied. Some representative samples with different kinds of nonhomogeneous properties are analyzed and the corresponding stress intensity factors (SIFs) are presented. It is shown that the PE mode is effective for investigating the crack problems of the FGMs with arbitrary properties.     

Penny shaped crack in a piezoelectric cylinder surrounded by an elastic medium subjected to combined in-plane mechanical and electrical loads

The fracture problem of a penny shaped crack in a piezoelectric ceramic cylinder surrounded by an infinite elastic medium under in-plane normal mechanical and electrical loads is considered with the electric continuous boundary conditions on the crack surface. By using the potential theory and Hankel transform, a system of dual integral equations is obtained, and expressed to a Fredholm integral equation of the second kind. The mechanical and electrical field equations and all sorts of field intensity factors of mode I are obtained, and the numerical values of various field intensity factors for PZT-6B piezoelectric ceramic surrounded by several different elastic media are graphically shown for a uniform load and a ring-shaped load, respectively. And the effects of the size of the piezoelectric cylinder and the elastic material properties on various field intensity factors are obtained.     

Mode III fracture problem of a cracked FGPM surface layer bonded to a cracked FGPM substrate

This work deals with the mode III fracture problem of a cracked functionally graded piezoelectric surface layer bonded to a cracked functionally graded piezoelectric substrate. The cracks are normal to the interface and the electro-elastic material properties are assumed to be varied along the crack direction. Potential and flux types of boundary condition are assigned on the edge of the surface layer. The problem under the assumptions of impermeable and permeable cracks can be formulated to the standard singular integral equations, which are solved by using the Gauss–Chebyshev technique. The effects of the boundary conditions, the material properties and crack interaction on the stress and electric displacement intensity factors are discussed.     

一种基于SHTB的Ⅱ型动态断裂实验技术

冲击剪切载荷作用下动态断裂韧性的测定是材料力学性能和断裂行为研究中重要组成部分.为了测定材料的Ⅱ型动态断裂韧性,许多学者采用不同的试样与实验方法进行了实验,但限于实验条件,裂纹断裂模式往往是I+Ⅱ复合型,而不是纯Ⅱ型,因而不能准确测得材料的Ⅱ型动态断裂韧性.鉴于此,本文基于分离式霍普金森拉杆(split Hopkinson tension bar,SHTB)实验技术,提出一种改进的紧凑拉伸剪切(modified compact tension shear,MCTS)试样,通过夹具对MCTS试样施加约束,从而保证试样按照纯Ⅱ型模式断裂.采用实验-数值方法对MCTS试样动态加载过程进行分析,将实验测得的波形输入有限元软件ANSYS-LSDYNA,得到了裂纹尖端应力强度因子-时间曲线,并与紧凑拉伸剪切(compact tension shear,CTS)试样进行了对比.同时采用数字图像相关法进行了实验,验证了有限元分析结果.结果表明,MCTS试样在整个加载过程中K_I K_Ⅱ,裂纹没有张开;而CTS试样在同样的加载过程中K_IK_Ⅱ,出现裂纹张开现象.这说明MCTS试样能够准确地测定材料的Ⅱ型动态断裂韧性,为材料动态力学测试提供了一种有效的实验技术.     

A UNIVERSAL CRACK EXTENSION CRITERION BASED ON THE EQUIVALENT STRESS GRADIENT:I.THEORY AND NUMERICAL VERIFICATION

In this contribution,the microscopic fracture mechanism and extension criterion for mixed type crack in ductile material under plane mixed mode loading are investigated in details.A universal extension...     

Basic solutions of multiple parallel symmetric mode-III cracks in functionally graded piezoelectric/piezomagnetic material plane

The Schmidt method is adopted to investigate the fracture problem of multiple parallel symmetric and permeable finite length mode-III cracks in a functionally graded piezoelectric/piezomagnetic material plane. This problem is formulated into dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. In order to obtain the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. The results show that the stress, the electric displacement, and the magnetic flux intensity factors of cracks depend on the crack length, the functionally graded parameter, and the distance among the multiple parallel cracks. The crack shielding effect is also obviously presented in a functionally graded piezoelectric/piezomagnetic material plane with mul- tiple parallel symmetric mode-III cracks.     

Growth of a mode II crack in an orthotropic plate made of a viscoelastic composite material

The growth of a straight mode II crack in a viscoelastic orthotropic plate is examined. The plate material is modeled by a viscoelastic anisotropic medium. The shear displacement in the fracture process zone is determined as a function of time using the corresponding elastic solution, the Volterra principle, and the method of operator continued functions. The time dependence of the crack length is constructed as integral equations of three phases of stable growth. The solution of these equations gives kinetic curves __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 9, pp. 89–97, September 2006.     

Anti-plane shear of an arbitrary oriented crack in a functionally graded strip bonded with two dissimilar half-planes

An internal crack located within a functionally graded material (FGM) strip bonded with two dissimilar half-planes and under an anti-plane load is considered. The crack is oriented in an arbitrary direction. The material properties of strip are assumed to vary exponentially in the thickness direction and two half-planes are assumed to be isotropic. Governing differential equations are derived and to reduce the difficulty of the problem dealing with solution of a system of singular integral equations Fourier integral transform is employed. Semi closed form solution for the stress distribution in the medium is obtained and mode III stress intensity factor (SIF), at the crack tip is calculated and its validity was verified. Finally, the effects of nonhomogeneous material parameter and crack orientation on the stress intensity factor are studied.     

Thermally induced fracture of a smart functionally graded composite structure

Discussed is the fracture behavior of a cracked smart actuator on a substrate under thermal load. The actuator is made of piezoelectric material with functionally graded material (FGM) properties. Integral transform method is used to reduce the problem to the solution of a set of singular integral equations and is solved numerically. This paper is completed by including graphical plots of the thermal flow, stress and electric displacement intensity factors around the crack for different crack positions and material gradients. Directions of crack initiation are also predicted by using the energy density criterion.     

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.     

LC4—M材料复合型韧断主要影响因素的分析

通过分析LC4－M铝合金材料在不同复合比载荷下的断裂试验结果,参照常规断裂现象,修正了一般断裂试验中认定裂纹启裂方向的方法,结合不同复合比下裂尖附近应力三维、主应力方向的计算分析,得到：在裂端的钝化变形区域、应力三维度的极大值处,对应于裂纹的启裂位置,即使在高韧性材料中发生剪切断裂的情况下也是如此;裂纹的启裂方向在拉伸断裂时与启裂点最大拉应力方向有关,在剪切断裂时启裂点最大剪应力方向有关,引起两种形式断裂的主要因数和破坏机理有很大不同。