有限厚度板穿透裂纹前缘附近三维弹性应力场分析

通过三维有限元计算来研究有限宽度、有限厚度含有穿透裂纹板的裂纹前缘应力场，从中找出应力强度因子与板的厚度、裂纹长度之间的关系，同时还分析了裂尖的三维约束程度和三维约束区的大小。分析结果表明：应力强度因子沿厚度的分布是不均匀的，应力强度因子的最大值及其位置与厚度有关；有限厚度板中面应力强度因子(KI)m-p及最大应力强度因子(KI)max均大于平面应力或平面应变的应力强度因子。对有限厚度裂纹问题，按平面应力或平面应变来考虑是不安全的；板中面的应力强度因子(KI)m-p及最大应力强度因子(KI)max是厚度B／a的函数；板的中面离面约束系数Tx最大，自由面(z=B)Tx=0。沿厚度方向裂尖附近的离面约束系数Tx也是z／B和B／a的函数，随着厚度的增加离面约束系数Tx增大，离中面越近离面约束系数Tx越大。Tx随着x的增大急剧减小，三维约束影响区域大小大约为板厚的一半，且裂纹长度a／W对应力强度因子沿厚度变化规律及Tx影响区域大小影响较小。     

三点弯曲试样动态应力强度因子计算研究

利用Hopkinson压杆对三点弯曲试样进行冲击加载,采集了垂直裂纹面距裂尖2mm和与裂纹面成60°距裂尖5mm处的应变信号。根据裂尖附近测试的应变信号计算试样的动态应力强度因子,并与有限元计算结果进行比较,结果表明由于裂尖有一段疲劳裂纹区,通过裂尖附近应变信号来计算动态应力强度因子时,如果裂尖位置确定不准及粘贴应变片位置不够准确对计算结果将带来很大影响。因此利用应变片法计算动态应力强度因子时,为了获得更准确的计算结果,在实验后应对试件裂纹面进行分析测量,重新确定裂尖位置,必要时需对应变片至裂尖距离进行修正后再计算动态应力强度因子值。     

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

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

V型缺口裂端的三维应力状态及约束分析

利用三维有限元方法,研究了有限厚度板中V型缺口根部穿透裂纹前沿的三维弹性应力场。对不同厚度、不同缺口张开角和裂纹长度对应力强度因子及裂尖附近的三维约束程度的影响进行了分析,同时还讨论了三维约束区的大小。研究结果显示：当缺口张开角小于60度时,不同缺口的应力强度因子和离面约束因子的分布基本一致,角度的影响不明显。应力强度因子是厚度的函数,板中面的应力强度因子随厚度的增加逐渐减小趋近干平面值,最大为1．08倍的平面值。当板厚超过15倍的缺口深度时,应力强度因子最大值将从中面转移至接近自由表面位置,距中面约0．4倍板厚。三维约束非常明显的区域在裂尖前沿约0.45倍厚度的范围内．二维到三维的过渡区在裂尖前沿1.5倍厚度的区域内;在中面上三维效应影响区最大,随着离中面距离的增加逐渐减小,在自由表面上降为0。     

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

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

正交异性钢桥面板肋-面板焊缝疲劳分析

基于焊缝的局部三维断裂力学模型和超重多轴货车的载荷谱,进行正交异性钢桥面板的肋-面板焊缝表面裂纹的疲劳寿命分析。采用Schwartz-Neuman交替法计算肋-面板焊缝处半椭圆表面裂纹的应力强度因子,基于裂尖反向塑性区模型考虑循环载荷中压应力对疲劳裂纹扩展的作用。正交异性钢桥面板的肋-面板焊缝的应力计算结果表明:超载货车作用下肋-面板焊缝处的横向应力峰值和应力幅都有明显增加;相比于标准疲劳荷载车,超载货车作用下肋-面板焊缝处半椭圆表面裂纹的裂纹扩展率增大了6.1倍;对应于正交异性钢桥面板的肋-面板焊缝处的拉-压循环应力,平面应变状态下的裂尖反向塑性区使裂纹扩展率增加了3.7倍;基于所得裂纹扩展速率,本研究给出仅在严重超载的五轴和六轴货车作用下正交异性钢桥面板的肋-面板焊缝疲劳寿命不足20年,远远低于桥梁的设计寿命。因此,考虑超载多轴货车的载荷谱和循环载荷中的压应力对肋-面板焊缝疲劳裂纹扩展的影响十分重要。     

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

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

复合材料断裂力学参量的实验测定

本文采用三维云纹干涉仪同步测量含裂纹复合材料弯曲板的面内位移场及离面位移场，并对离面位移场云纹图进行计算机自动采集和处理，编制了云纹数字图像处理软件。文中推导了正交异性复合材料弯曲问题的离面位移场表达式，得到了弯曲板裂纹尖端位移w与应力场强度因子K1的关系式。通过实测位移，结合裂尖位移场推导出应力强度因子K1值。     

三相压电复合本构模型中的弧形界面裂纹

深入研究了三相同心圆柱压电复合本构模型中的弧形绝缘界面裂纹问题。采用复势方法获得了该问题的级数形式的解答,并给出了应力、应变、电位移和电场强度等物理量在全场及界面上的分布,同时推导了裂尖处广义强度因子及裂面张开位移和裂面上电势差的表达式。具体计算表明该级数解答收敛迅速,同时显示出第三相混杂区的影响是不能忽略的。由于裂尖处应力奇异性为－1／2,则这种解答不会出现平面应变状态下界面裂纹裂尖处的振荡奇异性,从而不会产生违反物理实际的裂面相互嵌入现象,则该弹性解答也是建立了坚实的物理基础之上。     

裂尖大变形条件下的约束研究

通过面内应力比Tx和离面约束Tz在塑性区分布规律的深入分析，发现约束在塑性区不同区域的影响是不同的，在裂尖的钝化区，面内应力比Tx起重要作用，而在钝化区外的塑性区内，离面约束Tz起主要作用，在J－Q理论的基础上引入离面约束Tz可以很好地描述钝化区以外的应力应变分布，在钝化区内引入面内应力比Tx可以对裂尖前方的应力应变进行大变形分析和描述，研究表明综合两个约束的共同作用，可以很好地描述大变形条件下裂尖前方的应力分布。     

Tuning of non-uniform switch toughening in ferroelectric composites by an electric field

This paper deals with a mode III interfacial crack subject to anti-plane stress and in-plane electric fields. The analysis concentrates on the tuning of fracture toughness from non-uniform ferroelectric-ferroelastic domain switch-ing by an electric field. The electric loading changes the size of the asymmetric switching zone. Employing the weight function method, we obtain the electrically-dependent switch toughening for stationary and quasi-static growing interfacial cracks, respectively. Multi-domain solutions are derived for non-poled and fully-poled ferroelectric composites. Numer-ical results are presented on the electric field tuning of the critical applied stress intensity factor. The research provides ways to optimize fracture properties of ferroelectric compos-ites by altering the electric field.     

An energy-consistent fracture model for ferroelectrics

The fracture behavior of ferroelectrics has been intensively studied in recent decades, though currently a widely accepted fracture mechanism is still lacking. In this work, enlightened by previous experimental observations that crack propagation in ferroelectrics is always accompa-nied by domain switching, we propose a micromechanical model in which both crack propagation and domain switch-ing are controlled by energy-based criteria. Both electric energy and mechanical energy can induce domain switching, while only mechanical energy can drive crack propagation. Furthermore, constrained domain switching is considered in this model, leading to the gradient domain switching zone near the crack tip. Analysis results show that stress-induced ferroelastic switching always has a toughening effect as the mechanical energy release rate serves as the driving force for both fracture and domain switching. In compari-son, the electric-field-induced switching may have either a toughening or detoughening effect. The proposed model can qualitatively agree with the existing experimental results.     

PHASE FIELD SIMULATION OF DOMAIN SWITCHING IN FERROELECTRIC SINGLE CRYSTAL WITH ELECTRICALLY PERMEABLE AND IMPERMEABLE CRACKS

Domain switching around electrically permeable and impermeable cracks in ferroelectric single crystals subjected to a mechanical load is investigated by using a phase field model.It is found that the e...     

An I-integral method for crack-tip intensity factor variation due to domain switching in ferroelectric single-crystals

In the present study, an I-integral method is established for solving the crack-tip intensity factors of ferroelectric single-crystals. The I-integral combined with the phase field model is successfully used to investigate crack-tip intensity factor variations due to domain switching in ferroelectricity subjected to electromechanical loadings, which exhibits several advantages over previous methods based on small-scale switching. First, the shape of the switching zone around a crack tip is predicted by the time-dependent Ginzburg–Landau equation, which does not require preset energy-based switching criterion. Second, the I-integral can directly solve the crack-tip intensity factors and decouple the crack-tip intensity factors of different modes based on superimposing an auxiliary state onto an actual state. Third, the I-integral is area-independent, namely, the I-integral is not affected by the integral area size, the polarization distributions, or domain walls. This makes the I-integral applicable to large-scale domain switching. To this end, the electro-elastic field intensity factors of an impermeable crack in PbTiO₃ ferroelectric single crystals are evaluated under electrical, mechanical, and combined loading. The intensity factors obtained by the I-integral agree well with those obtained by the extrapolation technique. From numerical results, the following conclusions can be drawn with respect to fracture behavior of ferroelectrics under large-scale switching. Under displacement controlled mechanical loading, the stress intensity factors (SIFs) decrease monotonically due to the domain switching process, which means a crack tip shielding or effective switching-induced toughening occurs. If an external electric field is applied, the electric displacement intensity factor (EDIF) increases in all cases, i.e., the formed domain patterns enhance the electric crack tip loading. The energy release rate, expressed by the crack-tip J-integral, is reduced by the domain switching in all examples, which underlines the switching-induced-toughening effect. In contrast, under stress controlled load, the SIF evolves due to large-scale switching to a stable value, which is higher than the non-switching initial value, i.e., fracture is promoted in this case.     

Toughening under non-uniform ferro-elastic domain switching

The existing models of switch-toughening seldom consider the effect of non-uniform ferro-elastic domain switching in the vicinity of a crack. To explore this issue, an evolution law for the volume fraction of the switched portion under applied electromechanical loading is established from the minimum energy principle. Based on this law, a switching model capable of dealing with the non-uniform distribution of switching strain is developed. The domain switching zone is divided into a saturated inner core and an active surrounding annulus. Mono-domain solution of ferro-elastic toughening is obtained under the model of small scale domain switching. Toughening for ferroelectrics with different poling states is estimated via Reuss type approximation. Two sets of solutions are obtained according to spherical and cylindrical inclusions. The interval of toughening defined by these two models covers the range of experimental data. The same conclusion is reached for the size of the switching zone.     

Ⅱ-Ⅲ混合型裂纹相变增韧分析

本文采用压力敏感准则和松函数法对相变增韧陶瓷Ⅱ-Ⅲ混合型裂纹的增韧效应进行了理论预测。分别给出了静止裂纹和定常扩展裂纹相变塑性屏蔽的理论表达式。结果表明：相变对静止裂纹有甚微的负屏蔽效应，并随KⅢ/KⅡ的比值而波动；相变对定常扩展裂纹的增韧结果除与材料弹性模量、相变尾区高度和相变体积分数有关外，还与KⅢ/KⅡ的比值有关。     

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u _i , electric displacement D _i and volume fraction ρ _I of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction ρ _I of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch. The project supported by the National Natural Science Foundation of China (10572138).     

Toughening analysis of mode III crack in transformation toughened ceramics

The fracture toughness enhancement of mode III crack in transformation toughened ceramics such as Ce-TZP and Mg-PSZ is predicted by using a pressure sensitive transformation criterion and the weight function method. The theoretical analyses of the toughening effect for both stationary and steady-state growing crack are given, respectively. The results show that the toughening is only associated with crack growth and there is no toughening effect for stationary crack.The project was supported by the National Science Foundation of China     

陶瓷相变对Ⅰ—Ⅱ混合型裂纹增韧的理论分析

采用压力敏感准则和权函数法对相变增韧陶瓷I－II混合型裂纹的增韧效应进行了理论预测,分别给出了静止裂纹和定常扩展裂纹相变生屏蔽的理论表达式,并通过计算机进行了数值计算,结果表明：相变对静止裂纹有负屏蔽效应,并随KII/KI的比值增大而增大,对定常扩展裂纹的增韧结果除与材料弹性模量,相变尾区高度和相变体积分数有关外,还随KII／KI的比值增大而增大,说明相变对II型裂纹的增韧作用比对I型裂纹的增韧作用更显著。     

相变增韧陶瓷Ⅱ型裂纹增韧分析

本文采用压力敏感和权函数法，对相变增韧陶瓷Ⅱ型裂纹的增韧效应进行了理论预测。分别给出了静止裂纹和定常扩展裂纹相变性屏蔽的理论表达式，结果表明；相变对静止裂纹无增韧效应，纹尖端屏蔽来自于裂纹扩展尾区的贡献。