压电材料反平面界面裂纹的杂交元分析

基于奇异性电弹场数值特征解开发了一种新型反平面界面裂纹尖端单元。将新型单元与四节点压电P-S单元组装,求解从绝缘到导通的任意电边界条件下,压电结构反平面界面裂纹尖端电弹场的数值解。考察了层厚、载荷类型和裂纹面间电边界条件等对反平面界面裂纹尖端断裂参数的影响。算例证明新型单元能使P-S单元数显著降低,计算结果更为精确。     

双材料自由边界面端部问题的数值分析

对各向异性双材料自由边界面端部奇异性场问题进行了研究,利用有限元分析法所得到的各向异性双材料自由边界面端部的应力奇异性指数以及角分布函数,构造了一个自由边界面端部单元,据此建立了自由边界面端部奇异性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型.与四节点单元相结合,提出一种求解自由边界面端部广义应力强度因子的杂交元法.考核例结果表明:本文方法的数值解精度高,可应用于各向异性材料双材料自由边界面端部问题.     

基于有限元方法的航空发动机叶片应力强度因子计算

为研究叶片裂纹尖端的应力奇异性,以某型航空发动机压气机叶片为例,利用有限元方法研究了叶片裂纹尖端应力强度因子的计算方法,并研究了旋转叶片振动状态下裂尖应力强度因子随裂纹长度的变化规律。建立计算模型时,在裂纹尖端划分了三维奇异单元,在裂尖外围划分了过渡单元。计算结果表明：研究旋转叶片振动状态下的裂尖应力奇异性,仅利用I型应力强度因子就具有足够的精度;对于同一裂纹,绝大多数情况下叶盆面应力强度因子大于叶背面应力强度因子,故研究叶片应力强度因子时只需研究叶盆应力强度因子即可;随着裂纹扩展,叶盆面I型应力强度因子不断增大。本文的研究方法及结论为进一步研究叶片的裂纹扩展规律及损伤容限奠定了基础。     

矩形压电体中反平面裂纹的电弹性场

基于线性压电理论,本文获得了含有中心反平面裂纹的矩形压电体中的奇异应力和电场。利用Fourier积分变换和Fourier正弦级数将电绝缘型裂纹问题化为对偶积分方程,并进一步归结为易于求解的第二类Fred-holm积分方程。获得了裂纹尖端应力、应变、电位移和电场的解析解,求得了裂纹尖端场的强度因子及能量释放率。分析了压电矩形体的几何尺寸对它们的影响。结果表明,对于电绝缘型裂纹,裂纹尖端附近的各个场变量都具有-1/2阶的奇异性,能量释放率与电荷载的方向及大小有关,并且有可能为负值。     

压电陶瓷中圆币形裂纹在横向剪力下的机—电耦合行为

以弹性位移分量和电热函数基本未知量时，横观各向同性压电介质三维问题的场方程可化为四个联立的二阶线性偏微分方程组，本文导出了用四个调和函数表示位移分量及电势函数的表达式，即得到了该场方程的势函数能通解，作为通解的应用举例，文中求解了圆币形裂纹受横向剪切载荷下圆币形裂纹的尖端场及应力、电位移强度因子均具有明显的机－电耦合性质，而应力和电位移分量在裂尖仍具有－１／２的奇异性。     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

基于数字散斑相关方法测定Ⅰ型裂纹应力强度因子

提出了一种通过数字散斑相关方法测定金属材料Ⅰ型裂纹尖端位置和应力强度因子的实验方法.实验采用疲劳试验机对含Ⅰ型缺口的Cr12MoV钢试件预制裂纹,通过数字散斑相关方法测试试件在三点弯曲加载条件下裂纹的扩展过程及裂尖区域的位移场.将位移场数据代入裂尖位移场方程组,采用牛顿-拉普森方法求解含未知参量的裂尖非线性位移场方程组,计算裂尖位置和应力强度因子.实验结果表明,采用该方法可以准确地测定金属材料Ⅰ型裂纹应力强度因子、裂尖位置及裂纹扩展长度,解决了以往研究中因不能准确测定裂纹尖端位置,而无法准确计算Ⅰ型裂纹裂尖断裂参数的难题,揭示了金属材料裂纹扩展过程中应力强度因子演化特征.     

压电薄膜表面贯穿裂纹的有限元分析

基于有限元软件ANSYS数值模拟,计算了激光作用下的压电薄膜表面贯穿裂纹外场应力强度因子和电位移强度因子,并且研究了90°畴变所诱致的畴变增韧行为。首先,求解无裂纹压电薄膜在激光作用下的热-力-电响应,将求得的应力和电位移场反向作用于裂纹面,求解裂纹尖端处的外场应力和电位移强度因子,然后基于小范围畴变理论求解了90°畴变所致的屏蔽应力强度因子。讨论了薄膜表面裂纹的外场应力强度因子、电位移强度因子及屏蔽应力强度因子随激光作用时间和裂纹位置的变化关系,从而预测压电薄膜体系在加热工作状况下的裂纹扩展和断裂行为。     

用围线积分法计算压电材料面外剪切问题的应力强度因子

本文把Ｂｅｔｉ功互等原理推广到压电材料的面外剪切问题中，并且根据Ｐａｋ的压电材料Ⅲ型裂纹问题复势解，给出了其裂端位移、电势、应力和电位移的渐近解及相应的辅助场具体形式。然后，把有限元数值解作为真实平衡状态，把推导出的辅助场作为辅助平衡状态，利用围线积分法计算出了压电材料Ⅲ型裂纹问题的应力强度因子ＫⅢ和电强度因子ＫⅣ。算例表明，计算结果与理论解符合得很好     

用围线积分法计算压电材料面外剪切问题的应力强度因子

本文把Ｂｅｔｔｉ功互等原理推广到压电材料的面外剪切问题中，并且根据Ｐａｋ的压电材料Ⅲ型裂纹问题复势解，给出了其裂端位移，电势，应力和电位移的渐近解及相应的辅助场具体，然后，把有限元数值作为真实平衡状态，把推导出的辅助场作为辅助平衡状态，利用围线积分法出了压电材料Ⅲ型裂纹问题的应力强度因子ＫⅡ和电强度因子ＫⅣ。算例表明，计算结果与理论解符合得很好。     

???仯??????????????????????о?

采用有限元方法对表面含有两个不同大小半椭圆裂纹的有限厚矩形板在拉伸载荷作用下 进行参数化求解，得到两裂纹取不同尺寸和位置坐标时附属裂纹(尺寸较小裂纹)对主 导裂纹(尺寸较大裂纹)前沿参考点($\theta= 0, \pi/2, \pi $)处的应力强 度因子的影响系数$\beta$, 然后结合神经网络技术建立了相邻裂纹尺寸和位置参 数到主导裂纹前沿点处$\beta$的多变量非线性映射关系.     

Permeable cracks between two dissimilar piezoelectric materials

An interfacial crack with electrically permeable surfaces between two dissimilar piezoelectric ceramics under electromechanical loading is investigated. An exact expression for singular stress and electric fields near the tip of a permeable crack between two dissimilar anisotropic piezoelectric media are obtained. The interfacial crack-tip fields are shown to consist of both an inverse square root singularity and a pair of oscillatory singularities. It is found that the singular fields near the permeable interfacial crack tip are uniquely characterized by the real valued stress intensity factors proposed in this paper. The energy release rate is obtained in terms of the stress intensity factors. The exact solution of stress and electric fields for a finite interfacial crack problem is also derived.     

压电圆柱形夹杂对邻近裂纹三维应力场的影响

研究了压电复合材料薄板中压电圆柱形夹杂与邻近宏观钝裂纹间的相互作用。重点分析了外加电场，裂尖与压电圆柱形夹杂间韧带长度对裂尖三维应力场的影响。计算结果表明：在不同的外加电场作用下，压电体不仅能改变裂尖张开应力的大小，还能改变其分布。所得结果对进一步探讨线弹性介质中裂纹的启裂控制有参考价值。     

Mixed mode crack initiation in piezoelectric ceramic strip

When piezoelectric ceramics are subjected to mechanical and electrical load, they can fracture prematurely due to their brittle behavior. Hence, it is important to know the electro–elastic interaction and fracture behavior of piezoelectric materials. The problem of a through crack in a piezoelectric strip of finite thickness is studied in this paper. Fourier transforms are used to reduce the problem to the solution of singular integral equations. The model technique can solve for polarization in an arbitrary direction and material anisotropy. Numerical values of the crack-tip field amplification for a piezoelectric strip under in-plane electromechanical loading are obtained. Energy density factor criterion is applied to obtain the maximum of the minimum energy density and direction of crack initiation. The influence of crack length and crack position on stress intensity and energy density factors is discussed.     

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.     

Combination of fracture and damage mechanics for numerical failure analysis

A new approach for the analysis of crack propagation in brittle materials is proposed, which is based on a combination of fracture mechanics and continuum damage mechanics within the context of the finite element method. The approach combines the accuracy of singular crack-tip elements from fracture mechanics theories with the flexibility of crack representation by softening zones in damage mechanics formulations. A super element is constructed in which the typical elements are joined together. The crack propagation is decided on either of two fracture criteria; one criterion is based on the energy release rate or the J-integral, the other on the largest principal stress in the crack-tip region. Contrary to many damage mechanics methods, the combined fracture⧹damage approach is not sensitive to variations in the finite element division. Applications to situations of mixed-mode crack propagation in both two- and three-dimensional problems reveal that the calculated crack paths are independent of the element size and the element orientation and are accurate within one element from the theoretical (curvilinear) crack paths.     

A conducting arc crack between a circular piezoelectric inclusion and an unbounded matrix

The present paper investigates the problem of a conducting arc crack between a circular piezoelectric inclusion and an unbounded piezoelectric matrix. The original boundary value problem is reduced to a standard Riemann–Hilbert problem of vector form by means of analytical continuation. Explicit solutions for the stress singularities δ=−(1/2)±iε are obtained, closed form solutions for the field potentials are then derived through adopting a decoupling procedure. In addition, explicit expressions for the field component distributions in the whole field and along the circular interface are also obtained. Different from the interface insulating crack, stresses, strains, electric displacements and electric fields at the crack tips all exhibit oscillatory singularities. We also define a complex electro-elastic field concentration vector to characterize the singular fields near the crack tips and derive a simple expression for the energy release rate, which is always positive, in terms of the field concentration vector. The condition for the disappearance of the index ε is also discussed. When the index ε is zero, we obtain conventionally defined electro-elastic intensity factors. The examples demonstrate the physical behavior and the correctness of the obtained solution.     

Closed-form solution for two collinear cracks in a piezoelectric strip

Under the permeable electric boundary condition, the problem of two collinear anti-plane shear cracks lying at the mid-plane of a piezoelectric strip is investigated. By using the Fourier transform, the associated problem is reduced to a singular integral equation. Solving the resulting equation analytically, the electro-elastic field intensity factors and energy release rates at the inner and outer crack tips can be determined in explicit form. Numerical results for PZT-5H piezoelectric ceramic are also presented graphically. The results reveal that the effect of electric field on crack growth in piezoelectric materials is dependent on applied elastic displacement.     

Plastic yielding on inclined slip-planes at a crack tip

Plastic yield at crack tips on singular slip-planes, inclined to the crack plane, has been studied under plane-strain conditions for combined tension, hydrostatic stress, and in-plane shear. The singular integral equation, which represents the equilibrium condition of edge dislocations on the slip-planes, is transformed into a Fredholm integral equation in order to avoid difficulties that occur with its numerical solution. Results are presented for the slip-band length, the plastic crack-tip opening displacement, stress fields, and crack-opening contours. A series expansion of the results obtained numerically confirms approximate analytical expressions given by J.R. Rice (1974), up to the third-order in the applied stresses. The results of finite element methods agree with values of the crack-tip opening displacement obtained here to within 10 per cent. Ahead of the crack tip, the principal tensile stresses exceed the principal shear stresses by a factor of 10, approximately.     

改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.