力电耦合扩展无网格伽辽金法在压电柔性臂断裂分析中的研究

为提高含裂纹压电柔性臂在断裂分析中的求解精度,基于压电材料断裂力学理论,建立了压电柔性臂的力学分析模型。将描述裂纹尖端奇异性的位移场函数和电场函数引入到传统无网格伽辽金法中,提出了力电耦合扩展无网格伽辽金法。与传统无网格伽辽金法相比,本方法只需要较小的影响域来描述裂尖奇异场,并且节点影响域不会被裂纹线影响,不要可视准则和衍射准则,提高了计算精度。由虚拟裂纹闭合法推导了压电材料能量释放率,讨论了不同高斯点密度对强计算结果的影响。与解析解、有限元法的计算结果进行了比较,在高斯点个数为6×6时,扩展无网格伽辽金法的计算误差为1.88%,明显好于有限元的计算误差2.5%。数值算例结果表明本方法具有较高的计算精度。     

CTS试样Ⅰ/Ⅱ混合型断裂特性数值计算研究

采用CTS试样研究Ⅰ/Ⅱ混合型断裂特性计算裂纹前缘应力强度因子时可采用解析公式,一旦裂纹发生扩展,解析公式便不再适用.文中采用有限元法研究紧凑拉伸剪切(CTS)试样在Ⅰ/Ⅱ平面混合型加载下的裂纹扩展行为.采用ANSYS建立CTS试样Ⅰ/Ⅱ混合型测试系统有限元模型,为模拟真实受力状态,在CTS试样-销-扇型夹具以及扇型夹具-销-U型夹具之间分别建立接触对进行接触力学分析.通过与解析公式结果进行对比验证了该数值方法的可靠性.采用最大环向应力准则(MTS),模拟了CTS试样不同加载角度下的裂纹扩展路径,获得了裂纹扩展路径中应力强度因子随裂纹长度的变化曲线,解释了裂纹扩展路径不与外载荷方向垂直的原因.结合文中计算结果,在CTS试样Ⅰ/Ⅱ混合型裂纹扩展速率实验测得裂纹长度与寿命的关系曲线a-N的基础上,便可得到材料Ⅰ/Ⅱ型混合型裂纹扩展速率曲线.     

裂纹起裂方向及扩展路径的数值模拟

采用FRANC2D断裂力学分析软件分析了含裂纹构件的断裂问题,通过在裂纹尖端附近采用1/4奇异等参元,很好的刻画了裂纹尖端场的1/√r奇异性,通过对若干典型断裂问题的计算,得到了裂纹尖端的断裂参数,模拟了有限板中心斜裂纹和边斜裂纹的扩展过程,并对裂纹的起裂方向进行了预测,与实验结果进行对比具有较好的一致性.分析结果表明,FRANC2D软件对线弹性材料的断裂问题计算结果精度很高,可以模拟裂纹的扩展路径.利用该软件的计算结果对进一步深入分析断裂问题有很大帮助,可以为工程断裂分析提供参考.     

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

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

求解I型裂纹构元J积分的半解析方法

表征裂纹尖端应力应变场程度的J积分是一个定义明确、理论严密的弹塑性断裂力学基础参量. 目前J积分的计算主要是依靠塑性因子法和有限元法,但对各类裂纹构元获得J积分以及载荷-位移关系的解析公式以实现材料断裂韧性理论预测和材料测试是断裂力学的重要和困难的任务. 以J积分为参量的材料断裂测试中应用最广的是I型裂纹试样的断裂韧性测试. 本文在平面应变条件下,针对断裂韧性测试中使用的6种I型裂纹构元,基于能量等效假设,提出了J积分-载荷和载荷-位移的工程半解析统一表征方法,进而结合有限元分析的少量计算获得J积分-载荷和载荷-位移关系的半解析公式待定参数. 分析表明,6种I型裂纹构元的J积分-载荷和载荷-位移统一公式的预测结果与有限元结果吻合良好. 新提出的J积分-载荷工程半解析公式包含了材料的弹性模量、应力强度系数和应变硬化指数,能够广泛适应不同的材料,且运用该公式能够方便获取任意载荷点对应的J积分值. 应用新方法可便于获得各类I型裂纹构元的J积分-载荷和载荷-位移工程半解析公式.      

CTS试样I/II混合型断裂特性数值计算研究

采用CTS试样研究I/II混合型断裂特性计算裂纹前缘应力强度因子时可采用解析公式,一旦裂纹发生扩展,解析公式便不再适用。文中采用有限元法研究紧凑拉伸剪切（CTS）试样在I/II平面混合型加载下的裂纹扩展行为。采用ANSYS建立CTS试样I/II混合型测试系统有限元模型,为模拟真实受力状态,在CTS试样-销-扇型夹具以及扇型夹具-销-U型夹具之间分别建立接触对进行接触力学分析。通过与解析公式结果进行对比验证了该数值方法的可靠性。采用最大环向应力准则（MTS）,模拟了CTS试样不同加载角度下的裂纹扩展路径,获得了裂纹扩展路径中应力强度因子随裂纹长度的变化曲线,解释了裂纹扩展路径不与外载荷方向垂直的原因。结合文中计算结果,在CTS试样I/II混合型裂纹扩展速率实验测得裂纹长度与寿命的关系曲线a-N的基础上,便可得到材料I/II型混合型裂纹扩展速率曲线。     

平面幂律塑性Ⅰ型裂纹尖端应力场全解

在航空航天、船舶、石油管道和核电等领域,服役结构或部件在长期极端条件下运行,不可避免地会产生裂纹,因此,为研究含裂纹结构的准静态断裂行为,必须了解裂纹尖端附近区域的应力应变场特点.对于幂律材料裂纹构元,研究平面应变和平面应力条件下Ⅰ型裂纹尖端应力场的解析分布.基于能量密度等效和量纲分析,推导了能量密度中值点代表性体积单元(representative volume element, RVE)的等效应力解析方程,并定义其为应力因子,进而针对有限平面应变和平面应力紧凑拉伸(compact tension, CT)试样和单边裂纹弯曲(single edge bend, SEB)试样,以应力因子作为应力特征量,并构造用于表征裂尖等效应力等值线的蝶翅轮廓式和扇贝轮廓式三角特殊函数,提出描述幂律塑性条件下平面I型裂纹尖端应力场的半解析模型.该半解析模型形式简单,对CT和SEB试样的裂尖应力场的预测结果与有限元分析的结果比较表明,两者之间均密切吻合,模型公式可直接用于预测Ⅰ型裂纹尖端应力分布,方便于断裂安全评价和理论发展.     

功能梯度压电板条中电绝缘型运动裂纹的电弹性场

基于三维弹性理论和压电理论，对材料系数按指数函数规律分布的功能梯度压电板条中的反平面运动裂纹问题进行了求解。利用Fourier积分变换方法将电绝缘型运动裂纹问题化为对偶积分方程，并进一步归结为易于求解的第二类Fredholm积分方程。通过渐近分析，获得了裂纹尖端应力、应变、电位移和电场的解析解，给出了裂纹尖端场各个变量的角分布函数，并求得了裂纹尖端场的强度因子，分析了压电材料物性梯度参数、几何尺寸及裂纹运动速度对它们的影响。结果表明，对于电绝缘型裂纹，功能梯度压电板条中运动裂纹尖端附近的各个场变量都具有-1／2阶的奇异性；当裂纹运动速度增大时，裂纹扩展的方向会偏离裂纹面。     

层状压电陶瓷致动器中力电耦合场奇异性的数值分析

首先推导了不同压电材料界面裂纹尖端处的扇形区域内包含基本方程、裂纹面D-P边界条件和交界面处边界条件的弱形式。通过假设力电耦合位移场(位移和电势)与到裂纹尖端距离的(λ 1)次方成正比,可以得到一个分析压电材料裂纹尖端处力电耦合场奇异性的特殊的一维有限元列式。该一维有限元列式只需对扇形区域在角度方向上离散,最后的总体方程为一个关于λ的二次特征根方程。探讨了层状压电陶瓷致动器中可能出现奇异力电耦合场的部位的裂纹面边界条件及交界面处边界条件,进而将该一维有限元法进行推广,用于研究了这些部位的力电耦合场的奇异性。通过数值算例与相应的精确解的比较表明该方法是正确的,而且仅用很少单元就可以得到非常精确的结果。     

求解Ⅰ型裂纹构元J积分的半解析方法

表征裂纹尖端应力应变场程度的J积分是一个定义明确、理论严密的弹塑性断裂力学基础参量.目前J积分的计算主要是依靠塑性因子法和有限元法,但对各类裂纹构元获得J积分以及载荷-位移关系的解析公式以实现材料断裂韧性理论预测和材料测试是断裂力学的重要和困难的任务.以J积分为参量的材料断裂测试中应用最广的是Ⅰ型裂纹试样的断裂韧性测试.本文在平面应变条件下,针对断裂韧性测试中使用的6种Ⅰ型裂纹构元,基于能量等效假设,提出了J积分-载荷和载荷-位移的工程半解析统一表征方法,进而结合有限元分析的少量计算获得J积分-载荷和载荷-位移关系的半解析公式待定参数.分析表明,6种Ⅰ型裂纹构元的J积分-载荷和载荷-位移统一公式的预测结果与有限元结果吻合良好.新提出的J积分-载荷工程半解析公式包含了材料的弹性模量、应力强度系数和应变硬化指数,能够广泛适应不同的材料,且运用该公式能够方便获取任意载荷点对应的J积分值.应用新方法可便于获得各类Ⅰ型裂纹构元的J积分-载荷和载荷-位移工程半解析公式.     

Rectangular tensile sheet with a center notch root crack

This paper deals with the rectangular tensile sheet with a center notch crack. Such a crack problem is called a center notch crack problem for short. By using a hybrid displacement discontinuity method (a boundary element method) proposed recently by Yan, two center notch models are analyzed in detail. By changing the geometrical forms and parameters of the center notch, and by comparing the SIFs of the center notch crack problem with those of the center cracked plate tension specimen (CCT), which is a model frequently used in fracture mechanics, the effect of the geometrical forms and parameters of the center notch on the stress intensity factors (SIFs) of the center cracked plate tension specimen, is revealed. Some geometric characterestic parameters are introduced here, which are used to formulate the notch length and the branch crack length, which are to be determined in mechanical machining of the center cracked plate tension specimen. So we can say that the geometric characterestic parameters and the formulae used to determine the notch length and the branch crack length presented in this paper perhaps have some guidance role for mechanical machining of the center cracked plate tension specimen. In addition, the numerical investigation proves that the conventional angular notched specimen is much less sensitive to the size of notch than is the circular notched specimen.     

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.     

Closed form solution for two unequal collinear semi-permeable straight cracks in a piezoelectric media

The problem of two unequal collinear straight cracks weakening a poled transversely isotropic piezoelectric ceramic is addressed under semi-permeable electric boundary conditions on the crack faces. The plate has been subjected to combined in-plane normal(to the faces of the cracks) mechanical and electric loads. Problem is formulated employing Stroh formalism and solved using complex variable technique. The elastic field, electric field and energy release rate are obtained in closed analytic form. A case study is presented for poled PZT-5H cracked plate to study the effect of prescribed mechanical load, electric load, inter-crack distance and crack lengths on crack arrest parameters stress intensity factor (SIF), electric displacement intensity factor (EDIF) and mechanical and total energy release rates (ERR). Moreover a comparative study is done of impermeable and semi-permeable crack face boundary conditions on SIF, EDIF and ERR, and results obtained is presented graphically. It is observed that the effect of dielectric medium in the crack gap cannot be ignored.     

Analysis of electric boundary condition effects on crack propagation in piezoelectric ceramics

There are three types of cracks: impermeable crack, permeable crack and conducting crack, with different electric boundary conditions on faces of cracks in piezoelectric ceramics, which poses difficulties in the analysis of piezoelectric fracture problems. In this paper, in contrast to our previous FEM formulation, the numerical analysis is based on the used of exact electric boundary conditions at the crack faces, thus the common assumption of electric impermeability in the FEM analysis is avoided. The crack behavior and elasto-electric fields near a crack tip in a PZT-5 piezoelectric ceramic under mechanical, electrical and coupled mechanical-electrical loads with different electric boundary conditions on crack faces are investigated. It is found that the dielectric medium between the crack faces will reduce the singularity of stress and electric displacement. Furthermore, when the permittivity of the dielectric medium in the crack gap is of the same order as that of the piezoelectric ceramic, the crack becomes a conducting crack, the applied electric field has no effect on the crack propagation. The project supported by the National Natural Science Foundation of China (19672026, 19891180)     

Crack energy density of a piezoelectric material under general electromechanical loading

Crack energy density is considered and used as a possible fracture parameter in piezoelectricity under arbitrary electromechanical remote loads. The closed-form solution of a crack in a piezoelectric infinite plate subjected to general static electromechanical loading is obtained through a method alternative to the more common Stroh’s formalism. This analytical method, which is based on the spectral theorem of linear algebra, involves a transformation of similarity induced by the fundamental matrix in order to express the equations governing the problem in terms of complex potentials. The application of the mechanical boundary condition of stress-free crack and of one of the three considered electric boundary conditions (impermeable, permeable or semipermeable) leads then to the formulation of a Hilbert problem whose solution yields the stress and displacement fields. The crack energy density factors for mixed mode are then calculated under different mechanical and electrical loadings, as well as under different electric boundary conditions. The non-singular terms of the stress expressions are retained as well. The definition of the minimum energy density fracture criterion, as proposed by Sih, is given, and the influence of load biaxiality and positive or negative applied electric field on the criterion results is analyzed. The prediction of the incipient branching angle as from the energy density approach is also compared to that arising from the maximum circumferential stress theory for a mixed mode loading condition. Numerical results and graphs are presented and discussed for a PZT-4 piezoelectric ceramic.     

Applicability of the crack faces thermoelectric boundary conditions for thermopiezoelectric materials

The applicability and effect of the crack surfaces thermoelectric boundary conditions in thermopiezoelectric fracture mechanics problem are discussed by using the finite thickness notch approach. The stress and electric displacement intensity factors at the notch tips, and thermal flux and electric displacement inside the notch are derived in closed-form. The numerical results are compared with the ideal crack solutions. It is found that the electrically impermeable crack boundary condition assumption is reasonable if the flaw in the material is a notch with finite width, and the thermal conductivity of air or vacuum inside the crack must be considered.     

Effects of electric field on crack growth for a penny-shaped dielectric crack in a piezoelectric layer

The assumptions of impermeable and permeable cracks give rise to significant errors in determining electro-elastic behavior of a cracked piezoelectric material. The former simply imposes that the permittivity or electric displacement of the crack interior vanishes, and the latter neglects also the effects of the dielectric of an opening crack interior. Considering the presence of the dielectric of an opening crack interior and the permeability of the crack surfaces for electric field, this paper analyzes electro-elastic behavior induced by a penny-shaped dielectric crack in a piezoelectric ceramic layer. In the cases of prescribed displacement or prescribed stress at the layer surfaces, the Hankel transform technique is employed to reduce the problem to Fredholm integral equations with a parameter dependent nonlinearly on the unknown functions. For an infinite piezoelectric space, a closed-form solution can be derived explicitly, while for a piezoelectric layer, an iterative technique is suggested to solve the resulting nonlinear equations. Field intensity factors are obtained in terms of the solution of the equations. Numerical results of the crack opening displacement intensity factors are presented for a cracked PZT-5H layer and the effect of applied electric field on crack growth are examined for both cases. The results indicate that the fracture toughness of a piezoelectric ceramic is affected by the direction of applied electric fields, dependent on the elastic boundary conditions.     

Applicability of the crack-face electromagnetic boundary conditions for fracture of magnetoelectroelastic materials

This paper discusses the different electromagnetic boundary conditions on the crack-faces in magnetoelectroelastic materials, which possess coupled piezoelectric, piezomagnetic and magnetoelectric effects. A notch of finite thickness in these materials containing air (or vacuum) is also addressed. Four ideal crack-face electromagnetic boundary condition assumptions, that is, (a) electrically and magnetically impermeable crack, (b) electrically impermeable and magnetically permeable crack, (c) electrically permeable and magnetically impermeable crack and (d) electrically and magnetically permeable crack, are investigated separately. The influence of notch thickness on the field intensity factors at notch tips and the electromagnetic field inside the notch are obtained in closed-form. The results are compared with the ideal crack solutions. Applicability of crack-face electromagnetic boundary condition assumptions is discussed.     

基于新型裂尖杂交元的压电材料断裂力学研究

提出了一种裂尖邻域杂交元模型，将其与标准杂交应力元结合来求解压电材料裂纹尖 端的奇性电弹场和断裂参数的数值解．裂纹尖端杂交元的建立步骤为：1) 利用高次内插有限元特征法求解特征问题，得到反映裂尖奇异性电弹场状况的特 征值和特征角分布函数；2) 利用广义Hellinger-Reissner变分泛函以及特征问题的解来建立裂尖邻域杂交元模型．该 方法求解电弹场时，摒弃了传统有限元方法中裂尖奇异性场需要借助解析解的做法，也避免 了单纯有限元方法中需要在裂尖端部进行高密度单元划分．采用PZT5板中心裂纹问题 作为考核例，数值结果显示了良好的精确性．作为进一步应用，求解了含中心界面裂纹 的PZT4-PZT5两相压电材料的应力强度因子和电位移强度因子．所有的算例都考虑 了3种裂纹面电边界条件．     

Crack energy density in arbitrary direction for piezoelectrics and its characteristic features in electromechanical mixed mode crack

The concepts of crack energy density (CED) and its derivatives in arbitrary direction were established for piezoelectric material and, keeping their application to mixed mode fracture in mind, the characteristic features of them as fracture parameters were investigated based on the approximate equations for CED and its derivatives. That is, CED and its derivatives in arbitrary direction are defined first and separation into their each mode contribution is made. Subsequently, path independent integral expressions of them are derived, and then using them, approximate equations of each mode contribution of CED are obtained concretely for the case where linear singular solution is known. The resulting equations are then used to investigate the effects of electric field and electrical boundary condition on CED and its derivatives. An infinite piezoelectric plane with a crack inclined with respect to the poling direction is considered as a numerical example. Mode I contribution of mechanical CED is mainly employed as a possible fracture parameter for the study and it was shown that applied electric field significantly influences on fracture parameters especially for the impermeable crack perpendicular to the poling direction. The effect of electric field has the tendency to decrease as crack inclination angle increases. It was also found that, even for the impermeable crack perpendicular to the poling direction, crack propagation could be deviated from self-similar direction under a strong negative electric field, and this fact is qualitatively consistent with an existing experimental observation. For the ideally sharp crack with no width, impermeable and Hao and Shen type boundary conditions are admissible showing qualitative agreement with experimental results, but exact boundary condition is not suitable and finally consistent with permeable boundary condition.