双材料界面裂纹的改进广义有限差分法

采用数值方法进行断裂力学分析时,裂纹尖端奇异区域处理的好坏直接关系到最终断裂力学参数的求解精度。与传统均匀介质不同,复合材料界面裂纹渐近位移和应力场表现出剧烈的振荡特性,许多用于表征经典的平方根和负平方根物理场渐近性的传统方法也因此失效。论文提出了一种改进的广义有限差分法,该方法基于多元函数泰勒级数展开和移动最小二乘法的思想,将节点变量的各阶导数由相邻点集函数的加权线性累加来近似,具有无网格、无数值积分、数据准备简单、稀疏矩阵快速求解等优点。为提高该方法求解断裂力学问题的计算精度和数值稳定性,论文引入了裂尖奇异区域局部点簇的自动创建技术和一种基于局部点簇几何尺寸的矩阵正则化算法。数值算例表明,所提算法稳定,效率高,在不增加计算量的前提下,显著提高了裂尖近场力学参量和断裂力学参数的求解精度和数值稳定性。     

应用半权函数法计算界面裂纹的应力强度因子

应用半权函数法求解双材料界面裂纹的应力强度因子，得到以半权函数对参考位移与应力加权积分的形式表示的应力强度因子。针对特征值为复数λ的双材料界面裂纹裂尖应力和位移场，设置与之对应特征值为-λ的位移函数，即半权函数。半权函数的应力函数满足平衡方程，应力应变关系，界面的连续条件以及在裂纹面上面力为0；半权函数与裂纹体的几何尺寸无关，对边界条件没有要求。由功的互等定理得到应力强度因子KⅠ和KⅡ的积分形式表达式。本文计算了多种情况下界面裂纹应力强度因子的算例，与文献结果符合得很好。由于裂尖应力的振荡奇异性已经在积分中避免，只需考虑绕裂尖远场的任意路径上位移和应力，即使采用该路径上较粗糙的参考解也可以得到较精确的结果。     

求解界面裂纹应力强度因子的围线积分法

本文基于Ｂｅｔｔｉ功互等定理和双材料界面裂纹辅助场，提出了一种求解界面裂纹应力强度因子的方法，即远场围线积分法。此方法与积分径的选择无关，用有元元法计算出远离裂纹尖端的位移场和应力场，应可通过计算绕裂尖围线的积分，精确地给出界面裂纹应力强度因子ＫＩ和ＫⅡ。     

各向异性双材料界面裂纹扩展裂尖场

应用界面断裂力学理论和Stroh方法,研究了广义平面变形下动态裂纹沿着各向异性双材料界面扩展时的裂尖奇异应力及动态应力强度因子.双材料界面的动态裂尖区域特性主要由两个实矩阵W和D确定,且裂尖奇异应力和动态应力强度因子可以由包含这两个矩阵的柯西奇异积分方程确定,同时给出了动态应力强度因子和能量释放率的显示表达式.算例得出当裂纹以小速度扩展时,裂尖振荡因子ε与静态时几乎相同,当界面裂纹扩展速度接近瑞利波速时,ε趋于无穷大;同时得出应力强度因子及能量释放率随裂纹扩展速度的变化关系.     

夹杂对两相材料界面裂纹的影响

根据界面上应力和位移的连续条件,得到了单向拉伸状态下,含有椭圆夹杂的无限大双材料组合板的复势解。进一步通过求解Hilbert问题,得到了含有夹杂和半无限界面裂纹的无限大板的应力场,并由此给出了裂尖的应力强度因子K。计算了夹杂的形状、夹杂的位置、夹杂的材料选取以及上、下半平面材料与夹杂材料的不同组合对裂尖应力强度的影响。计算结果表明夹杂到裂尖的距离和夹杂材料的性质对K影响较大,对于不同材料组合,该影响有较大差异。夹杂距裂尖较近时,会对K产生明显屏蔽作用,随着夹杂远离裂尖,对K的影响也逐渐减小。另外,软夹杂对K有屏蔽作用,硬夹杂对K有反屏蔽作用,而夹杂形状对K几乎没有影响。     

双压电体界面上的电偶极子和裂纹5

求得了压电体双材料界面上的孤立二维电偶极子的解析解，结果表明某点电偶极子激发的应力－电位移场与该点到电偶极子的距离的平方成反比。研究了压电体双材料界面上的电偶极子对裂纹的作用，得到了问题的闭合解。在电偶极子的作用下，界面裂纹裂尖近区应力－电位移仍具有r^-1/2 iεα的振荡奇异性，文中求得裂尖应力强度因子，当电偶极子距裂尖距离ρ很近时，裂尖应力强度因子与ρ^-3/2-iεα成比例。     

黏弹性材料界面裂纹应力场奇性分析

研究两半无限大黏弹性体间Griffith界面裂纹在简谐载荷作用下裂纹尖端动应力场的奇异特性.通过引入裂纹张开位移和裂纹位错密度函数,相应的混合边值问题归结为一组耦合的奇异积分方程.渐近分析表明裂尖动应力场的奇异特征完全包含在奇异积分方程的基本解中.通过对基本解的深入分析发现黏弹性材料界面裂纹裂尖动应力场具有与材料参数和外载荷频率相关的振荡奇异特性.以标准线性固体黏弹材料为例讨论了材料参数和载荷频率对奇性指数和振荡指数的影响.     

裂纹面无摩擦接触的弹塑性界面裂尖应力场

本文通过对弹塑性幂硬化双材料界面裂纹尖端应力场的高阶渐近分析，获得了裂纹面无摩擦接触的裂尖一阶和二阶应力场解答，位移场在界面处呈现交叉匹配是本文解答的一个重要特点．最后结果表明，当界面上下材料的硬化指数之差大于１时，（即ｎ１－ｎ２＞１时），二阶应力场角分布为一常数解，而当０＜ｎ１－ｎ２≤１时，二阶应力场角分布函数则随θ变化而变化。     

裂纹面无磨擦接触的弹塑性界面裂尖应力场

本文通过对弹塑性幂硬化双材料界面裂纹尖端应力场的高阶渐近分析，获得了裂纹面无磨擦接触的裂尖一阶和二阶应力场解答，位移场在界面处呈现交叉匹配是本文解答的一个重要特点，最后结果表明，当界面上下材料的硬化指数之差大于１时（即ｎ１－ｎ２＞１时），二阶应力场角分布为一常数解；而当０＜ｎ１－ｎ２≤１时，二阶应力场角分布函数则随θ变化而变化。     

热释电材料问题的通解与界面裂纹

该文讨论了热释电材料中的热弹性问题的一般解，进而求解了共线界面裂纹问题．利用Ｓｔｒｏｈ方法，把热释电材料的热弹性界面裂纹问题化为一向量形式的Ｈｉｌｂｅｒｔ问题，求出这一Ｈｉｌｂｅｒｔ问题的通解，进而求得了热释电材料热弹性界面裂纹的闭合解，得到了温度、热流、位移、电势、应力和电位移的全场解，得到了裂纹张开位移及电势差的精确表达式．在此基础上，还求得了均匀热释电体中单个热弹性裂纹裂尖场，单个界面裂纹裂尖场以及点热源与界面裂纹的作用．此外，该文还对界面裂纹顶点附近的端部场作了渐近分析．     

Analysis method of planar interface cracks of arbitrary shape in three-dimensional transversely isotropic magnetoelectroelastic bimaterials

Using the fundamental solutions for three-dimensional transversely isotropic magnetoelectroelastic bimaterials, the extended displacements at any point for an internal crack parallel to the interface in a magnetoelectroelastic bimaterial are expressed in terms of the extended displacement discontinuities across the crack surfaces. The hyper-singular boundary integral–differential equations of the extended displacement discontinuities are obtained for planar interface cracks of arbitrary shape under impermeable and permeable boundary conditions in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. An analysis method is proposed based on the analogy between the obtained boundary integral–differential equations and those for interface cracks in purely elastic media. The singular indexes and the singular behaviors of near crack-tip fields are studied. Three new extended stress intensity factors at crack tip related to the extended stresses are defined for interface cracks in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. A penny-shaped interface crack in magnetoelectroelastic bimaterials is studied by using the proposed method.The results show that the extended stresses near the border of an impermeable interface crack possess the well-known oscillating singularity r^?1/2±iε or the non-oscillating singularity r^?1/2±κ. Three-dimensional transversely isotropic magnetoelectroelastic bimaterials are categorized into two groups, i.e., ε-group with non-zero value of ε and κ-group with non-zero value of κ. The two indexes ε and κ do not coexist for one bimaterial. However, the extended stresses near the border of a permeable interface crack have only oscillating singularity and depend only on the mechanical loadings.     

用子域边界元法研究各向异性材料中的界面裂纹

用子域边界元法研究各向异性材料中的界面裂纹，在边界元公式中，采用了带特征根的基本解，以增量形式的边界积分方程为基础，通过二次等参元及国分之一面力奇异远离散化处理，可以得到各子域的代数方程组，依据凝集技术，可得到仅含有子域公共边界及裂纹边界未知量的求解方程组，通过迭代法，可以寻求到每种载荷作用下的裂纹所处的真实状态，然后，由文献「２」中的方法求解界面裂纹的应力强度因子。结果表明，子域边界方法是正确的     

On the crack-tip stress singularity of interfacial cracks in transversely isotropic piezoelectric bimaterials

In this paper, characteristics of the interface crack-tip stress and electric displacement fields in transversely isotropic piezoelectric bimaterials are studied. The authors have proven, within the framework of the generalized Stroh formalism for piezoelectric bimaterials, that there is no coexistence of the parameters (oscillating) and κ (non-oscillating) in the interface crack-tip generalized stress field for all transversely isotropic piezoelectric bimaterials. This leads to the classification of piezoelectric bimaterials into one group that exhibits the oscillating property in the interface crack-tip generalized stress field and the other that does not. Fifteen (15) pair-combinations of six (6) piezoelectric materials PZT-4, PZT-5H, PZT-6B, PZT-7A, P-7, and BaTiO₃, which are commonly used in practice, are numerically analyzed in this study, and the results backup the above theoretical conclusions. Moreover, the associated eigenvectors for such material systems (with either =0 or κ=0) are also obtained numerically, and the result show that there still exist four linear independent associate eigenvectors for each bimaterial.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

Interface crack-tip generalized stress field and stress intensity factors in transversely isotropic piezoelectric bimaterials

Transversely isotropic piezoelectric (TIP) bimaterials with an impermeable interface crack have been classified [Int. J. Frac. 119 (2003) L41] into two classes corresponding to the vanishing of the two singularity parameters or κ. It is shown in the present paper that the related eigenvalue problems for either =0 or κ=0 are not degenerate. The crack-tip generalized stress fields are obtained subsequently. A new definition of crack-tip intensity factors is presented for interface cracks in practical TIP bimaterial of practical interest. Such defined intensity factors are real numbers, which dominate the maximum crack-tip stress singularity and do not generate any phase angle change under any dimension system transformation for physical quantities.     

轴对称环形片状界面裂纹问题分析

讨论受拉伸载荷作用的轴对称环形片状界而裂纹问题．该问题归结为求解一组超奇异积分-微分方程．方程中的未知位移间断近似表示为基本密度函数与多项式之积,其中基本密度函数考虑到问题的对称性用二维界面裂纹精确解表示．在圆形片状裂纹的情况下,数值结果与现有理论解作比较的结果表明,数值结果与相应界面圆形片状裂纹和均质体圆形片状裂纹的精确解均吻合得很好．文中以图表形式给出应力强度因子与材料组合和几何条件之间的关系．     

双材料界面Ⅰ型扩展裂纹尖端的弹黏塑性场

双材料界面中存在材料黏性效应, 对界面裂纹尖端场的分布和界面本身性能 的变化起着重要的影响. 考虑裂纹尖端的奇异性, 建立了双材料界面扩展裂纹尖端的弹黏塑 性控制方程. 引入界面裂纹尖端的位移势函数和边界条件, 对刚性-弹黏塑性界面I型界面 裂纹进行了数值分析, 求得了界面裂纹尖端应力应变场, 并讨论了界面裂纹尖端场随各影响 参数的变化规律. 计算结果表明, 黏性效应是研究界面扩展裂纹尖端场时的一个主要因素, 界面裂纹尖端为弹黏塑性场, 其场受材料的黏性系数、马赫数和奇异性指数控制.     

Crack problems of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in piezoelectrics or on the interfaces of piezoelectric bimaterials. A class of boundary problems involving many cracks is also solved. For homogeneous materials it is found that the normal electric displacementD ₂ induced by the crack is constant along the crack faces which depends only on the applied remote stress field. Within the crack slit, the electric fields induced by the crack are also constant and not affected by the applied electric field. For the bimaterials with realH, the normal electric displacementD ₂ is constant along the crack faces and electric fieldE ₂ has the singularity ahead of the crack tip and a jump across the interface. The project is supported by the National Natural Science Foundation of China(No. 19704100) and the Natural Science Foundation of Chinese Academy of Sciences(No. KJ951-1-201).     

FATIGUE GROWTH MODELING OF MIXED-MODE CRACK IN PLANE ELASTIC MEDIA

This paper presents an extension of a displacement discontinuity method with cracktip elements （a boundary element method） proposed by the author for fatigue crack growth analysis in plane elastic media under mixed-mode conditions. The boundary element method consists of the non-singular displacement discontinuity elements presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or right crack-tip element is placed locally at the corresponding left or right crack tip on top of the non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. Crack growth is simulated with an incremental crack extension analysis based on the maximum circumferential stress criterion. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not required because of an intrinsic feature of the numerical approach. Crack growth is modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characteristics of some related elements are adjusted according to the manner in which the boundary element method is implemented. As an example, the fatigue growth process of cracks emanating from a circular hole in a plane elastic plate is simulated using the numerical simulation approach.