插值矩阵法分析正交各向异性板切口应力奇异性

基于切口尖端附近区域位移场渐近展开,提出了分析正交各向异性复合材料板切口奇异性的新方法.将位移场的渐近展开式的典型项代入弹性板的基本方程,得到关于正交各向异性板切口奇异性指数的一组非线性常微分方程的特征值问题;再采用变量代换法,将非线性特征问题转化为线性特征问题,用插值矩阵法求解获得的正交各向异性板切口若干阶应力奇异性指数和相应特征函数.该法可由相应的特征角函数对板切口的平面应力和反平面奇异特征值加以区分,并将计算结果与现有结果对照,表明了该文方法的有效性.     

微分求积法分析平面接头应力奇异性北大核心CSCD

对于双材料平面接头问题提出了一个分析应力奇性指数的新方法:微分求积法(DQM).首先,将平面接头连接点处位移场的径向渐近展开格式代入平面弹性力学控制方程,获得了关于应力奇性指数的常微分方程组(ODEs)特征值问题.然后,基于DQM理论,将ODEs的特征值问题转化为标准型广义代数方程组特征值问题,求解之可一次性地计算出双材料平面接头连接点处应力奇性指数,同时,一并求出了接头连接点处相应的位移和应力特征函数.数值计算结果说明该文DQM计算平面接头连接点处应力奇性指数的结果是正确的.     

预制V型裂纹尖端应力强度因子的研究

本文提出了爆炸载荷作用下预制V型裂纹的复变应力函数,并用Westergaard方法推导了预制V型裂纹尖端的应力场和位移场,从而得到了V型裂纹尖端的应力强度因子.爆破试验结果表明了公式的正确性.     

求解双材料裂纹结构全域应力场的扩展边界元法

在线弹性理论中,复合材料裂纹尖端具有多重应力奇异性,常规数值方法不易求解.该文建立的扩展边界元法(XBEM)对围绕尖端区域位移函数采用自尖端径向距离r的渐近级数展开式表达,其幅值系数作为基本未知量,而尖端外部区域采用常规边界元法离散方程.两方程联立求解可获得裂纹结构完整的位移和应力场.对两相材料裂纹结构尖端的两个材料域分别采用合理的应力特征对,然后对其进行计算,通过计算结果的对比分析,表明了扩展边界元法求解两相材料裂纹结构全域应力场的准确性和有效性.     

弹性-粘塑性材料Ⅰ型动态扩展裂纹尖端场的渐近解

提出了一种新的弹性-粘塑性模型用于分析Ⅰ型动态扩展裂纹尖端的应力应变场.给出了适当的位移模式,推导了渐近方程并且给出了数值解.分析和计算表明:对于低粘性情况,裂纹尖端场具有对数奇异性;对于高粘性情况,渐近方程无解.分析比较表明该结果具有高玉臣提出的单参数解的所有优点,并且消除了粘性区随裂纹扩展而移动的不足.     

弹性—粘塑性材料I型动态扩展裂纹尖端场的渐近解

提出了一种新的弹性－粘塑性模型用于分析I型动态扩展裂纹尖端的应力应变场。给出了适当的位移模式，推导了渐近方程并且给出了数值解。分析和计算表明：对于低粘性情况，裂纹尖端场具有对数奇异性；对于高粘性情况，渐近方程无解。分析比较表明该结果具有高压臣提出的单参数解的所有优点，并且消除了粘性区随裂纹扩展而移动的不足。     

正交异性双材料的Ⅱ型界面裂纹尖端场

通过引入含16个待定实系数和两个实应力奇异指数的应力函数,再借助边界条件,得到了两个八元非齐次线性方程组.求解该方程组,在双材料工程参数满足适当条件下,确定了两个实应力奇异指数.根据极限唯一性定理,求出了全部系数,得到了应力函数的表示式.代入相应的力学公式,推出了当特征方程组两个判别式都小于0时,每种材料的裂纹尖端应力强度因子、应力场和位移场的理论解.裂纹尖端附近的应力和位移有混合型断裂特征,但没有振荡奇异性和裂纹面相互嵌入现象作为特例,当两种正交异性材料相同时,可以推出正交异性单材料Ⅱ型断裂的应力奇异指数、应力强度因子公式、应力场、位移场表示式.     

相变增韧陶瓷平面应变Ⅰ型定常扩展裂纹的渐近分析(Ⅱ)*

本文采用一种考虑相交剪切变形的陶瓷材料本构关系,对平面应变Ⅰ型定常扩展裂纹尖端场进行渐近分析.给出了裂纹尖端附近环形域内的应力、速率分布以及应力奇异性指数.对不同材料参数下的变化规律进行了详细的分析和讨论.     

反平面剪切动态扩展裂纹尖端的弹-粘塑性场

采用Bingham弹性-粘塑性模型对反平面剪切动态扩展裂纹尖端的应力应变场进行了渐近分析．给出了适当的位移模式、推导了渐近方程并且给出了数值解．分析和计算表明对于低粘性情况,裂纹尖端场具有对数奇异性．对于高粘性情况,裂纹尖场具有幂奇异性A·D2对于临界情况,两种奇异性可以相互转换．揭示了粘性在裂纹尖端场研究中的重要作用．     

相变增韧陶瓷平面应变I型定常扩展裂纹的渐近分析（Ⅱ）

本文采用一种考虑相变剪切变形的陶瓷材料本构关系，对平面应Ｉ型定常扩展裂纹尖端场进行渐近分析。给出了裂纹尖端附近环形域内的应力，速率分布以及应力奇异性指数，对不同材料参数下的变化规律进行了详细的分析和讨论。     

Complete elastic–plastic stress asymptotic solutions near general plane V-notch tips

An efficient method is developed to determine the multiple term eigen-solutions of the elastic–plastic stress fields at the plane V-notch tip in power-law hardening materials. By introducing the asymptotic expansions of stress and displacement fields around the V-notch tip into the fundamental equations of elastic–plastic theory, the governing ordinary differential equations (ODEs) with the stress and displacement eigen-functions are established. Then the interpolating matrix method is employed to solve the resulting nonlinear and linear ODEs. Consequently, the first four and even more terms of the stress exponents and the associated eigen-solutions are obtained. The present method has the advantages of greater versatility and high accuracy, which is capable of dealing with the V-notches with arbitrary opening angle under plane strain and plane stress. In the present analysis, both the elastic and the plastic deformations are considered, thus the complete elastic and plastic stress asymptotic solutions are evaluated. Numerical examples are shown to demonstrate the accuracy and effectiveness of the present method.     

Evaluation of the stress singularities of plane V-notches in bonded dissimilar materials

According to the linear theory of elasticity, there exists a combination of different orders of stress singularity at a V-notch tip of bonded dissimilar materials. The singularity reflects a strong stress concentration near the sharp V-notches. In this paper, a new way is proposed in order to determine the orders of singularity for two-dimensional V-notch problems. Firstly, on the basis of an asymptotic stress field in terms of radial coordinates at the V-notch tip, the governing equations of the elastic theory are transformed into an eigenvalue problem of ordinary differential equations (ODEs) with respect to the circumferential coordinate θ around the notch tip. Then the interpolating matrix method established by the first author is further developed to solve the general eigenvalue problem. Hence, the singularity orders of the V-notch problem are determined through solving the corresponding ODEs by means of the interpolating matrix method. Meanwhile, the associated eigenvectors of the displacement and stress fields near the V-notches are also obtained. These functions are essential in calculating the amplitude of the stress field described as generalized stress intensity factors of the V-notches. The present method is also available to deal with the plane V-notch problems in bonded orthotropic multi-material. Finally, numerical examples are presented to illustrate the accuracy and the effectiveness of the method.     

压电压磁复合材料中界面裂纹对弹性波的散射

利用Schmidt方法分析了压电压磁复合材料中可导通界面裂纹对反平面简谐波的散射问题．经过富里叶变换得到了以裂纹面上的间断位移为未知变量的对偶积分方程A·D2在求解对偶积分方程的过程中,裂纹面上的间断位移被展开成雅可比多项式的形式．数值模拟分析了裂纹长度、波速和入射波频率对应力强度因子、电位移强度因子、磁通量强度因子的影响A·D2从结果中可以看出,压电压磁复合材料中可导通界面裂纹的反平面问题的应力奇异性形式与一般弹性材料中的反平面问题应力奇异性形式相同．     

不同功能梯度压电压磁层状介质中共线界面裂纹的动态性能分挝

对不同功能梯度压电压磁层状介质中,共线界面裂纹对简谐应力波作用下的动态问题,进行了分析.经Fourier变换,使问题的求解转换为求解以裂纹面上位移间断为未知量的三重对偶积分方程,三重对偶方程可以采用Schmidt方法来求解,进而分析了功能梯度参数、入射波频率和层状介质厚度对应力、电位移和磁通量强度因子的影响.     

Analysis of the stress singularity for a bi-material V-notch by the boundary element method

A new algorithm coupling the boundary element technique with the characteristic expansion method is proposed for the computation of the singular stress field in the V-notched bi-material structure. After the stress asymptotic expansions are introduced into the linear elasticity equilibrium equations, the governing equations at the small sector dug out from the bi-material V-notch tip region are transformed into the ordinary differential eigen-equations. All the parameters in the asymptotic expansions except the combination coefficients can be achieved by solving the established eigen-equations with the interpolating matrix method. Furthermore, the conventional boundary element method is applied to modeling the remaining structure without the notch tip region. The combination coefficients in the asymptotic expansion forms can be computed by the discretized boundary integral equations. Thus, the singular stress field at the V-notch tip and the generalized stress intensity factors of the bi-material notch are successfully calculated. The accurate singular stress field obtained here is very useful in the evaluation of the fracture property and the fatigue life of the V-notched bi-material structure.     

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

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

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.     

A general method for constructing Timoshenko-type theories

The equations of the plane theory of for the elasticity bending of a long strip are reduced by the method of simple iterations to the solution of a system of two equations for the displacement of the axis of the strip and the shear stress. If the transverse load varies slowly along the strip, the resolvent equations reduce to a single equation that is identical to the classical equation for the bend of a beam. When a local load is applied, the resolvent equation acquires an additional singular term that is the solution of the equation for the shear stresses under the assumption that the displacement (deflection) is a function of small variability. The convergence of the solution in an asymptotic sense is demonstrated. The application of the method of simple iterations to the dynamic equations for the bending of a strip also leads to a system of two resolvent equations in the displacement of the axis of the strip and the shear stress. These equations reduce to a single equation that is identical with the well-known Timoshenko equation. Hence, the procedure for using the method of simple iterations that has been developed can be classified as a general method for obtaining Timoshenko-type theories. An equation is derived for the bending of a strip on an elastic base with an isolated functional singular part with two bed coefficients, corresponding to the transverse and longitudinal springiness of the base.     

Analysis of anti-plane interface cracks in one-dimensional hexagonal quasicrystal coating

The displacement discontinuity method is extended to study the fracture behavior of interface cracks in one-dimensional hexagonal quasicrystal coating subjected to anti-plane loading. The Fredholm integral equation of the first kind is established in terms of displacement discontinuities. The fundamental solution for anti-plane displacement discontinuity is derived by the Fourier transform method. The singularity of stress near the crack front is analyzed, and Chebyshev polynomials of the second kind are numerically adopted to solve the integral equations. The displacement discontinuities across crack faces, the stress intensity factors, and the energy release rate are calculated from the coefficients of Chebyshev polynomials. In combination with numerical simulations, a comprehensive study of influencing factors on the fracture behavior is conducted.