正交各向异性功能梯度材料Ⅲ型裂纹尖端动态应力场

研究了无限大正交各向异性功能梯度材料Griffith裂纹受反平面剪切冲击作用的问题.材料两个方向的剪切模量假定为成比例按特定梯度变化.通过采用积分变换-对偶积分方程方法,获得了裂纹尖端动态应力场.动态应力强度因子计算结果显示:增加剪切模量梯度或增加垂直于裂纹面方向的剪切模量可以抑制动态应力强度因子的幅值.     

各向异性材料动态裂纹扩展特性和动态应力强度因子

基于各向异性材料力学,研究了无限大各向异性材料中Ⅲ型裂纹的动态扩展问题．裂纹尖端的应力和位移被表示为解析函数的形式,解析函数可以表达为幂级数的形式,幂级数的系数由边界条件确定．确定了Ⅲ型裂纹的动态应力强度因子的表达式,得到了裂纹尖端的应力分量、应变分量和位移分量．裂纹扩展特性由裂纹扩展速度M和参数alpha反映,裂纹扩展越快,裂纹尖端的应力分量和位移分量越大;参数alpha对裂纹尖端的应力分量和位移分量有重要影响．     

基于裂纹尖端二阶弹性解的断裂过程区尺度

基于Westergaard应力函数裂纹尖端二阶弹性解,推导了裂纹尖端微裂区的轮廓线和特征尺寸的解析表达式;采用幂函数模型描述的拉应变软化模型,确定了在最大拉应力强度理论和最大拉应变强度理论下断裂过程区（FPZ）临界值的解析表达式;将基于Westergaard应力函数一阶弹性解及二阶弹性解、Muskhelishvili应力函数和Duan-Nakagawa模型确定的FPZ临界值进行了比较．结果表明裂纹尖端微裂区和FPZ临界值随着Poisson比的减小而增加并逐渐趋近于应用最大拉应力强度理论确定的结果;二阶弹性解确定的裂纹尖端微裂区和FPZ临界值大于一阶弹性解的值;FPZ临界值随着拉应变软化指数的增加而增加;二阶弹性解确定的FPZ临界值的精度远高于一阶弹性解确定的值.     

功能梯度条硬币型裂纹扭转冲击响应

研究非均匀条中硬币型裂纹的扭转冲击问题．材料的剪切模量假定按特定的梯度变化．采用Laplace 和Hankel 变换将问题化为求解Fredholm积分方程,通过将Bessel函数渐进展开获得裂纹尖端动态应力场．考查非均匀参数和功能梯度条高度对裂尖动态断裂行为的影响．动应力强度因子和能量密度因子的清晰表达式表明,作为裂纹扩展力,对于这里所研究的问题,二者是等价的．动应力强度因子的数值结果显示,增加剪切模量的非均匀参数可以抑制动应力强度因子的幅度,而条形域的高度对动态断裂特性的影响较小．     

一维六方压电准晶双材料界面共线裂纹问题北大核心CSCD

利用复变函数理论中的解析延拓、奇性主部分析和推广的Liouville定理,求解了一维六方压电准晶双材料在集中载荷作用下界面共线裂纹反平面弹性问题.导出了含有一条和两条有限长界面裂纹的封闭解,同时给出了裂纹尖端场强度因子(包含声子场和相位子场应力强度因子和电位移强度因子)的表达式.数值算例分析了外荷载与耦合系数之比对裂纹尖端场强度因子变化规律的影响.从数值结果中可以看出,当裂纹长度增加时,裂纹尖端场强度因子随之增加;应力强度因子随双材料耦合系数之比的增大而增大,电位移强度因子几乎不变;不同载荷作用下,裂纹尖端场强度因子随着裂纹长度改变时的变化趋势也不尽相同.研究结果可为压电准晶双材料的设计和制备提供一定的理论参考.     

压电压磁复合材料中一对平行裂纹对弹性波的散射

利用Schmidt方法对压电压磁复合材料中一对平行对称裂纹对反平面简谐波的散射问题进行了分析,借助富里叶变换得到了以裂纹面上的间断位移为未知变量的对偶积分方程．在求解对偶积分方程的过程中,裂纹面上的间断位移被展开成雅可比多项式的形式,最终获得了应力强度因子、电位移强度因子、磁通量强度因子三者之间的关系．结果表明,压电压磁复合材料中平行裂纹动态反平面断裂问题的应力奇异性与一般弹性材料中的动态反平面断裂问题的应力奇异性相同,同时讨论了裂纹间的屏蔽效应．     

楔型向错偶极子和裂纹的干涉效应

研究了晶体材料中一个楔型向错偶极子与裂纹的弹性干涉效应．运用复变函数方法获得了复势函数和应力场的封闭形式解答,导出了裂纹尖端应力强度因子和作用在向错偶极子中心点像力的解析表达式．获得了向错偶极子的位置、方向和偶臂长度对裂纹尖端应力强度因子的影响规律,并讨论了裂纹附近向错偶极子的平衡位置．结果表明向错偶极子靠近裂纹尖端时,对应力强度因子有明显的屏蔽或反屏蔽作用．     

压电材料中两平行不相等界面裂纹的动态特性研究

利用Schmidt方法,研究了压电材料中两个平行不相等的可导通界面裂纹对简谐反平面剪切波的散射问题．利用Fourier变换,使问题的求解转换为对两对以裂纹面张开位移为未知变量的对偶积分方程的求解．数值计算结果表明,动态应力强度因子及电位移强度因子受裂纹的几何参数、入射波频率的影响．在特殊情况下,与已有结果进行了比较分析．同时,电位移强度因子远小于不可导通电边界条件下相应问题的结果．     

位于两不同正交各向异性半平面间张开型界面裂纹的性能分析

利用Schmidt方法分析了位于正交各向异性材料中的张开型界面裂纹问题．经富立叶变换使问题的求解转换为求解两对对偶积分方程,其中对偶积分方程的变量为裂纹面张开位移．最终获得了应力强度因子的数值解．与以前有关界面裂纹问题的解相比,没遇到数学上难以处理的应力振荡奇异性,裂纹尖端应力场的奇异性与均匀材料中裂纹尖端应力场的奇异性相同．同时当上下半平面材料相同时,可以得到其精确解．     

垂直于双材料非完美界面Ⅱ型裂纹的断裂分析

研究了垂直于双材料非完美界面的Ⅱ型裂纹问题,采用线性弹簧模型模拟非完美界面．然后用Fourier积分变换方法把边值问题转化为求解具有Cauchy核的奇异积分方程,获得了裂纹两端应力强度因子的数值解．详细研究了问题的几种特例,并用数值实例分析了界面的非完美性对应力强度因子的影响．结果表明应力强度因子与界面参量有关并在完美界面和分离界面所对应的结果中变化．     

Transient response of two collinear dielectric cracks in a piezoelectric solid under inplane impacts

The paper is focused on the dynamic analysis of two collinear dielectric cracks in a piezoelectric material under the action of in-plane electromechanical impacts. Considering the dielectric permeability of crack interior, the electric displacements at the crack surfaces are governed by the jumps of electric potential and crack opening displacement across the cracks. The permeable and impermeable crack models are the limiting cases of the general one. The Laplace and Fourier transform techniques are further utilized to solve the mixed initial-boundary-value problem, and then to obtain the singular integral equations with Cauchy kernel, which are solved numerically. Dynamic intensity factors of stress, electric displacement and crack opening displacement are determined in time domain by means of a numerical inversion of the Laplace transform. Numerical results for PZT-5H are calculated to show the effects of the dielectric permeability inside the cracks, applied electric loadings and the geometry of the cracks on the fracture parameters in graphics. The observations reveal that based on the COD intensity factor, a positive electric field enhances the dynamic dielectric crack growth and a negative one impedes the dynamic dielectric crack growth in a piezoelectric solid.     

Transient thermoelastic response in a cracked strip of functionally graded materials via generalized fractional heat conduction

This work is devoted to analyzing a thermal shock problem of an elastic strip made of functionally graded materials containing a crack parallel to the free surface based on a generalized fractional heat conduction theory. The embedded crack is assumed to be insulated. The Fourier transform and the Laplace transform are employed to solve a mixed initial-boundary value problem associated with a time-fractional partial differential equation. Temperature and thermal stresses in the Laplace transform domain are evaluated by solving a system of singular integral equations. Numerical results of the thermoelastic fields in the time domain are given by applying a numerical inversion of the Laplace transform. The temperature jump between the upper and lower crack faces and the thermal stress intensity factors at the crack tips are illustrated graphically, and phase lags of heat flux, fractional orders, and gradient index play different roles in controlling heat transfer process. A comparison of the temperature jump and thermal stress intensity factors between the non-Fourier model and the classical Fourier model is made. Numerical results show that wave-like behavior and memory effects are two significant features of the fractional Cattaneo heat conduction, which does not occur for the classical Fourier heat conduction.     

On subcritical development of a shear crack in a composite with viscoelastic components

We present results of an investigation of the development of a transverse shear crack in a composite material with linearly viscoelastic components under external shear load. The solution is divided into the following two main stages: determination of the time dependence of the crack tip opening displacement and determination of the crack-growth kinetics as a result of the solution of integral equations. In the first stage, we use the solution of the corresponding elastic problem of determination of the crack opening displacement and the problem of determination of the effective moduli of the composite reinforced with unidirectional discrete fibers. Using the theoretically proved principle of elasto-viscoelastic analogy and the method of Laplace inverse transformation, we obtain a solution in a time domain. In the second stage, using the criterion of critical crack opening displacement for a transverse shear crack and an equation for the viscoelastic crack opening displacement of this crack, we construct an equation of crack growth. We present results of the numerical solution, which illustrate the influence of relations between the relaxation parameters of the materials of the components on the durability of the body with a crack.     

Transient response of a semi-infinite piezoelectric layer with a surface permeable crack

The transient response of a semi-infinite transversely isotropic piezoelectric layer containing a surface crack is analyzed for the case where anti-plane mechanical and in-plane electric impacts are suddenly exerted at the layer end. The integral transform techniques are used to reduce the associated mixed initial boundary value problem to a singular integral equation of the first kind, which can be solved numerically via the Lobatto–Chebyshev collocation technique. Dynamic field intensity factors are determined by employing a numerical inversion of the Laplace transform. The dynamic stress intensity factors are presented graphically and the effects of the material properties and geometric parameters are examined. Received: June 30, 2003     

An interface crack moving between magnetoelectroelastic and functionally graded elastic layers

A constant crack moving along the interface of magnetoelectroelastic and functionally graded elastic layers under anti-plane shear and in-plane electric and magnetic loading is investigated by the integral transform method. Fourier transforms are applied to reduce the mixed boundary value problem of the crack to dual integral equations, which are expressed in terms of Fredholm integral equations of the second kind. The singular stress, electric displacement and magnetic induction near the crack tip are obtained asymptotically and the corresponding field intensity factors are defined. Numerical results show that the stress intensity factors are influenced by the crack moving velocity, the material properties, the functionally graded parameter and the geometric size ratios. The propagation of the moving crack may bring about crack kinking, depending on the crack moving velocity and the material properties across the interface.     

条状功能梯度材料中偏心裂纹对反平面简谐波的散射问题

利用Schmidt方法研究了条状功能梯度材料中偏心裂纹对反平面简谐波的散射问题,裂纹垂直于条状功能梯度材料的边界．通过Fourier变换,问题可以转换为对一对未知变量是裂纹表面位移差的对偶积分方程求解．为了求解对偶积分方程,把裂纹表面的位移差展开为Jacobi多项式级数形式,进而得到了功能梯度参数、裂纹位置以及入射波频率对应力强度因子影响的规律．     

横观各向同性电磁弹性介质中裂纹对SH波的散射

研究横观各向同性电磁弹性介质中裂纹和反平面剪切波之间的相互作用．根据电磁弹性介质的平衡运动微分方程、电位移和磁感应强度微分方程,得到SH波传播的控制场方程．引入线性变换,将控制场方程简化为Helmholtz方程和两个Laplace方程A·D2通过Fourier变换,并采用非电磁渗透型裂面边界条件,得到了柯西奇异积分方程组．利用Chebyshev多项式求解积分方程,得到应力场、电场和磁场以及动应力强度因子的表达,并给出了数值算例．     

A singular perturbation approach to integral equations occurring in poroelasticity

Singular perturbation theory is used to solve the integral equationswhich occur when treating finite-length crack problems in porouselastic materials. The method provides the stress intensityfactors which characterize the near crack tip stress and displacementfields for small times. The method also gives the stress andpore pressure fields on the fracture plane for small times relativeto the diffusive time scale. In this paper, the authors treatcrack problems which are unmixed in the pore pressure boundarycondition on the fracture plane. The Abelian result that smalltimes correspond, in Laplace transform space, to large valuesof the transform variable is used to formulate the problemsin terms of a small parameter. Rescaling on this small parameterleads to inner problems which are eigensolutions of the semi-infiniteproblems treated earlier by the authors. The outer solutionsare given by elastic eigensolutions together with appropriatefluid dipole responses. These outer solutions give the completestress and pore pressure fields except in the neighbourhoodof the crack tips; in this region the outer solutions are asymptoticallymatched with inner solutions. The full outer solutions are givenhere as an asymptotic expansion for small times and enable thedevelopment of the outer fields to be followed in real time.A reciprocal theorem in Laplace transform space is used to checkthe small-time solutions. The inner problem is rescaled to asemi-infinite crack problem, so eigensolutions of this semi-infiniteproblem are used together with the known asymptotic behaviourof the real solution to identify the stress intensity factor.The stress intensity factor is then related to an integral involvingthe inner limit of the outer solution together with the eigensolutionof the semi-infinite problem. Using this integral, we recoverthe result for the stress intensity factor found using singularperturbation theory. A ‘nearly’ invariant integralanalogous to the invariant M integral used in elastostaticsis derived. Unfortunately, the poroelastic analogue is not invariant,although it is used to verify the small-time results.     

Magnetoelectroelastodynamic interaction of multiple arbitrarily oriented planar cracks

The problem of multiple arbitrarily oriented planar cracks in an infinite magnetoelectroelastic space under dynamic loadings is considered. An explicit solution to the problem is given in the Laplace transform domain in terms of suitable exponential Fourier integral representations. The unknown functions in the Fourier integrals are directly related to the Laplace transform of the jumps in the displacements, electric potential and magnetic potential across opposite crack faces and are to be determined by solving a system of hypersingular integral equations. Once the hypersingular integral equations are solved, the displacements, electric potential, magnetic potential and other quantities of interest such as the crack tip intensity factors may be easily computed in the Laplace transform domain and recovered in the physical space with the help of a suitable algorithm for inverting Laplace transforms.