Oblique edge crack in an anisotropic material under antiplane shear

An oblique edge crack in an anisotropic material under antiplane shear loadings is investigated. The antiplane problems are formulated based on a linear transformation method. An anisotropic solid containing an edge crack subjected to concentrated forces is first considered. The stress intensity factor for the edge crack with concentrated forces is obtained from the solution of the transformed edge crack in an isotropic material which is solved by using conformal mapping technique and complex function theory. The solution of the edge crack under concentrated loads is used to construct the stress intensity factor for the oblique edge crack in the anisotropic material subjected to antiplane distributed loads. Some numerical computations are carried out to calculate the stress intensity factors for the edge crack in inclined orthotropic materials subjected to point forces as well as distributed tractions.     

基于子模型法的带有表面裂纹钢丝应力强度因子研究

钢丝裂纹应力强度因子是进行钢丝疲劳断裂寿命评估、疲劳裂纹扩展分析和钢丝断裂强度评估等工作的重要参数。本文首先介绍了裂纹扩展分析软件FRANC3D,然后基于子模型法模拟研究了拉伸荷载作用下带有表面裂纹钢丝的应力强度因子,裂纹种类包括直线形裂纹和半圆形裂纹,最后拟合得到拉伸荷载作用下带表面裂纹钢丝的应力强度因子形状修正系数表达式,分析了利用该公式进行承载力评估时产生误差的原因。研究结果表明,利用子模型模拟分析拉伸荷载作用下带有表面裂纹的钢丝应力强度因子时计算精度高,计算速度快,对计算机硬件要求低;利用该方法得到的钢丝裂纹应力强度因子,在进行索承式桥梁吊索安全性能评估时,评估结果更精确。     

冲击载荷作用下裂纹动态响应的数值模拟

对垂直、剪切以及斜向等各种冲击载荷作用下裂纹的动态响应进行了数值模拟,得到了一系列随时间变化的动态应力场以及应变场图;根据其定义,计算出了相应的动态应力强度因子,进而分析了斜向载荷作用下裂纹起裂情况,并对最优断裂问题进行了阐述。     

A periodic array of cracks in a transversely isotropic magnetoelectroelastic material

Magnetoelectroelastic materials are inherently brittle and prone to fracture. Therefore, it is important to evaluate the fracture behavior of these advanced materials. In this paper, a periodic array of cracks in a transversely isotropic magnetoelectroelastic material is investigated. Hankel transform is applied to solve elastic displacements, electric potential and magnetic potential. The problem is reduced into a system of integral equations. Both impermeable and permeable crack-face electromagnetic boundary conditions assumptions are investigated. Quantities of the stress, electric displacement and magnetic induction and their intensity factor are obtained. Effect of the crack spacing on these quantities is investigated in details.     

平行于功能梯度材料夹层的币型裂纹起裂条件

分析了功能梯度材料中币型裂纹的扩展问题．裂纹平行于无限域中功能梯度材料夹层,受有与裂纹面成任意角度的拉应力．假定功能梯度材料夹层与两个半无限域均匀介质完全粘合,其弹性模量沿厚度方向变化．采用基于层状材料广义Kelvin基本解的边界元方法分析裂纹问题,给出了均布正应力和剪应力作用下裂纹的应力强度因子、将应力强度因子耦合于应变能密度断裂判据,讨论了裂纹体在拉伸应力作用下的起裂条件．     

The weight function for cracks in piezoelectrics

The weight function in fracture mechanics is the stress intensity factor at the tip of a crack in an elastic material due to a point load at an arbitrary location in the body containing the crack. For a piezoelectric material, this definition is extended to include the effect of point charges and the presence of an electric displacement intensity factor at the tip of the crack. Thus, the weight function permits the calculation of the crack tip intensity factors for an arbitrary distribution of applied loads and imposed electric charges. In this paper, the weight function for calculating the stress and electric displacement intensity factors for cracks in piezoelectric materials is formulated from Maxwell relationships among the energy release rate, the physical displacements and the electric potential as dependent variables and the applied loads and electric charges as independent variables. These Maxwell relationships arise as a result of an electric enthalpy for the body that can be formulated in terms of the applied loads and imposed electric charges. An electric enthalpy for a body containing an electrically impermeable crack can then be stated that accounts for the presence of loads and charges for a problem that has been solved previously plus the loads and charges associated with an unsolved problem for which the stress and electric displacement intensity factors are to be found. Differentiation of the electric enthalpy twice with respect to the applied loads (or imposed charges) and with respect to the crack length gives rise to Maxwell relationships for the derivative of the crack tip energy release rate with respect to the applied loads (or imposed charges) of the unsolved problem equal to the derivative of the physical displacements (or the electric potential) of the solved problem with respect to the crack length. The Irwin relationship for the crack tip energy release rate in terms of the crack tip intensity factors then allows the intensity factors for the unsolved problem to be formulated, thereby giving the desired weight function. The results are used to derive the weight function for an electrically impermeable Griffith crack in an infinite piezoelectric body, thereby giving the stress intensity factors and the electric displacement intensity factor due to a point load and a point charge anywhere in an infinite piezoelectric body. The use of the weight function to compute the electric displacement factor for an electrically permeable crack is then presented. Explicit results based on a previous analysis are given for a Griffith crack in an infinite body of PZT-5H poled orthogonally to the crack surfaces.     

Dynamic analysis of a penny-shaped dielectric crack in a magnetoelectroelastic solid under impacts

The transient response of a magneto-electro-elastic material with a penny-shaped dielectric crack subjected to in-plane magneto-electro-mechanical impacts is made. To simulate an opening crack with a dielectric interior, the crack-face electromagnetic boundary conditions are supposed to depend on the crack opening displacement and the jumps of electric and magnetic potentials across the crack. Four ideal crack-face electromagnetic boundary conditions involving a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions can be reduced. The Laplace and Hankel transform techniques are further utilized to solve the mixed initial-boundary-value problem. Three coupling Fredholm integral equations are obtained and solved by the composite Simpson's rule. Dynamic field intensity factors of stress, electric displacement, magnetic induction, crack opening displacement (COD), electric potential and magnetic potential are given in the Laplace transform domain. By means of a numerical inversion of the Laplace transform, numerical results are calculated to show the variations of the physical parameters of concern versus the normalized time in graphics. The effects of applied electric and magnetic loads on the dynamic intensity factors of stress and COD, and the dynamic energy release rate for a BaTiO₃-CoFe₂O₄ composite with a penny-shaped vacuum crack are discussed in detail.     

Eccentric crack in a piezoelectric strip bonded to half planes

We consider the problem of determining the singular stresses and electric fields in a piezoelectric ceramic strip containing an eccentric Griffith crack off the centre line bonded to two elastic half planes under anti-plane shear loading using the continuous crack-face condition. Fourier transforms are used to reduce the problem to the solution of two pairs of dual integral equations, which are then expressed to a Fredholm integral equation of the second kind. Numerical values on the stress intensity factor and energy release rate are obtained.     

Dielectric crack problem for a magnetoelectroelastic strip with functionally graded properties

The paper presents a fracture analysis for an electromagnetically dielectric crack in a functionally graded magnetoelectroelastic strip. It is considered that the material properties are varying exponentially along the width direction. Under the assumption of the in-plane magneto-electro-mechanical loadings, the dielectric crack is simulated by using the semi-permeable crack-face boundary conditions. The Fourier transform technique is applied to solve the boundary-value problem and four coupling singular integral equations are determined. A nonlinear system of algebraic equations is further derived and solved numerically to determine the electromagnetic field inside the crack. Then the field intensity factors of stress, electric displacement, and magnetic induction are given. Through the numerical computations, the effects of the material non-homogeneity and the permeability of crack interior on the electric displacement and the magnetic induction at the crack faces are studied. The variations of the intensity factors of stress, electric displacement, and magnetic induction versus the geometry of the crack, the strip width, and the material non-homogeneity are presented in graphics respectively.     

Dynamic crack propagation in a magneto-electro-elastic solid subjected to mixed loads: Transient Mode-III problem

The transient response of a Mode-III crack propagating in a magneto-electro-elastic solid subjected to mixed loads is investigated through solving the corresponding boundary-initial-value problem in both the cracked solid region and the interior fluid region with treatment of electro-magnetically permeable and impermeable crack face conditions in a unified way. The closed-form results for the dynamic field intensity factors are used to evaluate the dynamic energy release rate through the crack-tip dynamic contour integral. The permeability of the interior fluid region relative to the cracked solid region significantly affects the magneto-electro-mechanical coupling coefficient in the Bleustein–Gulyaev wave function and, consequently, the horizontal shear surface wave speed, the dynamic field intensity factors and the dynamic energy release rate. It is revealed from dynamic fracture mechanics analysis that the dynamic energy release rate thus obtained has an odd dependence on the dynamic electric displacement intensity factor and the dynamic magnetic induction intensity factor. It is also found that the horizontal shear surface wave speed provides the limiting velocity for the propagation of a Mode-III crack in a magneto-electro-elastic solid when there is only applied traction loading.     

An anti-plane propagating crack in a functionally graded piezoelectric strip bonded to a homogeneous piezoelectric strip

Summary A finite crack propagating at constant speed in a functionally graded piezoelectric strip (FGPS) bonded to a homogeneous piezoelectric strip is considered. It is assumed that the electroelastic material properties of the FGPS vary exponentially across the thickness of the strip, and that the bimaterial strip is under combined anti-plane mechanical shear and in-plane electrical loads. The analysis is conducted for the electrically unified crack boundary condition, which includes both the traditional permeable and the impermeable ones. By using the Fourier transform, the problem is reduced to the solution of Fredholm integral equations of the second kind. Numerical results for the stress intensity factor and the crack sliding displacement are presented to show the influences of the crack propagation speed, electric loads, FGPS gradation, crack length, electromechanical coupling coefficient, properties of the bonded homogeneous piezoelectric strip and crack location.     

Sudden deceleration or acceleration of an intersonic shear crack

Sudden jumps in the crack tip velocity were revealed by numerical simulation (in both continuum/cohesive element and molecular dynamics approaches) and experiments for rapid shear cracking. The cracking velocity may accelerate from a sub-Rayleigh speed to the intersonic range, or from an intersonic speed to a higher one, when the reflected impact wave reloads the crack tip. On the other hand, the cracking velocity may decelerate from an intersonic speed to a lower one or recede to the sub-Rayleigh range when the fracture driving force declines. The velocity change encountered during intersonic cracking plays a different role from that in the acceleration or deceleration of a subsonic crack. A crack propagating at an intersonic speed would leave a shear wave trailing behind. When the crack decelerates or accelerates, the effect of the trailing wave will lead to a transition period from one steady-state solution of crack tip singularity to another. This investigation aims at quantifying these processes. The full field solution of an intersonic mode II crack whose speed changed suddenly from one velocity (intersonic or subsonic) to another (intersonic or subsonic) is given in closed form. The solution is facilitated via superposing a series of propagating crack problems that are loaded by dislocations to seal the unwanted crack-face sliding or by concentrated forces moving at various speeds to negate the crack-face traction. In contrast to the subsonic solution, the results in the intersonic case indicate that the elastic fields around the crack tip depend on the deceleration or acceleration history that is traced back over a long time. Singularity matching dictates the jump that may actually take place.     

Torsional impact response of a penny-shaped crack in a functional graded strip

The torsional impact response of a penny-shaped crack in a nonhomogeneous strip is considered. The shear modulus is assumed to be functionally graded such that the mathematics is tractable. Laplace and Hankel transforms were used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Explicit expressions of both the dynamic stress intensity factor and the energy density factor were derived. And it is shown that, as crack driving force, they are equivalent for the present crack problem. Investigated are the effects of material nonhomogeneity and strip‘s highness on the dynamic fracture behavior.Numerical results reveal that the peak of the dynamic stress intensity factor can be suppressed by increasing the nonhomogeneity parameter of the shear modulus, and that the dynamic behavior varies little with the adjusting of the strip‘ s highness.     

冲击载荷作用下黑砂岩动态断裂参数的分形修正

基于分形理论研究了偏折裂纹扩展路径对动载荷作用下黑砂岩的动态断裂力学参数的测试误差影响作用,采用传统的分离式霍普金森压杆（split Hopkinson pressure bar, SHPB）实验装置对修正侧开单裂纹半孔板（improved single cleavage semi-circle specimen, ISCSC）试样进行动态冲击实验,随后采用裂纹扩展计进行裂纹起裂时间与裂纹扩展速度等动态断裂力学参数测试,采用分形理论对测试的裂纹扩展速度与动态应力强度因子进行修正,利用实验-数值法对黑砂岩的动态断裂韧度进行计算。研究结果表明,ISCSC构型构件能够有效应用于岩石材料动态裂纹扩展行为的研究,并发生了止裂现象,经分形修正的裂纹扩展速度与动态断裂韧度更接近实际裂纹动态扩展情况,修正前后得到黑砂岩材料的裂纹扩展速度误差为33.51%,动态断裂韧度最大误差为7.68%,说明利用分形理论对动态断裂韧度等动态断裂参数计算更合理。     

Moving crack for functionally grated material in an infinite length strip under antiplane shear

Considered is a Yoffe crack in an infinite strip of functionally grated material (FGM) subjected to antiplane shear. The shear moduli in two directions of FGM are assumed to be of exponential form. The dynamic stress intensity factor and strain energy density factor at the crack tip are obtained by using integral transforms and dual-integral equations. The numerical results show that the decrease of the strain energy density factor varies with the shear moduli gradient, and the increase of the strain energy density factor varies with the increase of the moving crack speed. The ratio of shear moduli in material vertical orientation has a great influence on the strain energy density factor.     

A penny-shaped interface crack between a functionally graded piezoelectric layer and a homogeneous piezoelectric layer

The problem of a penny-shaped interface crack between a functionally graded piezoelectric layer and a homogeneous piezoelectric layer is investigated. The surfaces of the composite structure are subjected to both mechanical and electrical loads. The crack surfaces are assumed to be electrically impermeable. Integral transform method is employed to reduce the problem to a Fredholm integral equation of the second kind. The stress intensity factor, electric displacement intensity factor and energy release rate are derived, some typical numerical results are plotted graphically. The effects of electrical loads, material nonhomogeneity and crack configuration on the fracture behaviors of the cracked composite structure are analyzed in detail.     

纤维复合材料复合型断裂正应力强度因子比准则

对纤维增强复合材料的裂纹起裂及开裂方向准则进行了研究，提出了复合型断裂的正应力强度因子比准则（Normal Stress Intensity Factor Raito Criterion).此准则是一种综合考虑了正应力强度因子和剪应力强度因子对裂纹起裂的推动的准则，并且不需要预先确定材料的特征尺度，使用较方便，且预测结果是很好的。     

Mode III stress intensity factors in an interfacial crack in dissimilar bonded materials

The antiplane shear deformation problem of two edge-bonded dissimilar isotropic wedges is considered. In the case when the sum of the two apex angles is equal to 2π, the problem reduces to that of two edge-bonded dissimilar materials with an interfacial crack subjected to concentrated antiplane shear tractions on the crack faces. An explicit expression is extracted for the stress intensity factor at the crack tip. In the special cases of different combinations of the apex angles, the obtained expression for the stress intensity factor may be simplified and relations of a simpler form are given for the stress intensity factor. It is shown that the stress intensity factor is dependent on the material properties as well as the geometry and loading. However, in special cases of equal apex angles as well as the case of similar materials the dependency of the stress intensity factor on the material properties disappears.     

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｇｒｅａｔａｔｔｅｎｔｉｏｎｓｈａｖｅｂｅｅｎｐａｉｄｔｏｔｈｅｒｅｓｅａｒｃｈｏｆＦｕｎｃｔｉｏｎａｌｌｙＧｒａｄｅｄＭａｔｅｒｉａｌｓ(ＦＧＭ).Ｆｒｏｍｔｈｅｖｉｅｗｐｏｉｎｔｓｏｆａｐｐｌｉｅｄｍｅｃｈａｎｉｃｓ,ＦＧＭａｒｅｎｏｎ＿ｈｏｍｏｇｅｎｅｏｕｓｓｏｌｉｄｓ.Ｔｈｅｎｏｎ＿ｈｏｍｏｇｅｎｅｉｔｙｏｆＦＧＭｈａｓａｇｒｅａｔｉｎｆｌｕｅｎｃｅｏｎｔｈｅｉｒｍｅｃｈａｎｉｃａｌｂｅｈａｖｉｏｒ,ｅｓｐｅｃｉａｌｌｙｗｈｅｎｔｈｅｃｏｍｐ…     

Crack shielding and material deterioration in damaged materials: an antiplane shear fracture problem

Summary A simple damage evolution model is proposed for a quasibrittle material in the framework of continuum damage mechanics. The model is used to obtain a closed form solution for a mode-III stationary crack under small scale damage conditions. It is found that the crack tip stress intensity factor is reduced, i.e., the crack is shielded by the damage. However, this shielding effect is completely offset by the material deterioration caused by the damage. It also holds for steady state crack growth. When the most effective shielding is reached for the stationary crack, the zone dominated by the stress intensity factor shrinks to the crack tip. The results for the antiplane shear problem should shed some light on the in- plane fracture problem. Received 4 August 1997; accepted for publication 7 October 1997