压电材料反平面裂纹自相似扩展分析

通过对耦合的波动方程和方程解耦,用自模拟方法研究了压电材料中反平面裂纹的自相似扩展问题。研究表明：对反平面问题,介质内的耦合场与裂纹扩展速度有关,在裂纹尖端有r^-1/2阶的奇异性;动态应力强度因子与电位称载荷有关,与静态结论不同;电位移强度因子与机械载 荷无关,与静态结果的表达形式一致。     

动态裂纹扩展中的分形效应

假设裂纹顶端沿着分形轨迹运动，建立了裂纹扩展的分形弯折（ｋｉｎｋｉｎｇ）模型来描述裂纹的动态扩展。根据这个模型，我们推导了分形裂纹扩展对劝态应力强度和裂纹速度的影响．动态应力强度因子与表观应力强度因子之比Ｋ（ｌ（ｔ），Ｖ）／Ｋ（Ｌ（ｔ），Ｏ）是表观裂纹速度Ｖ＿Ｏ，材料微结构参数（ｄ／Δａ），分维Ｄ和裂纹扩展路径的弯折角θ的函数。本文研究结果表明：在分形裂纹扩展中，表观（或量测）的裂纹速度Ｖ＿Ｏ很难接近Ｒａｙｌｅｉｇｈ波速Ｃ＿ｒ．动态断裂实验中Ｖ＿Ｏ明显低于Ｃ＿ｒ的原因可能是分形裂纹扩展效应所致。材料的微结构，裂纹扩展路径的分维和弯折角均很强地影响动态应力强度因子和裂纹扩展速度。     

功能梯度材料涂层平面运动裂纹分析

研究粘结于均匀材料基底上功能梯度材料涂层平面运动裂纹问题,假设功能梯度材料剪切模量和密度为坐标的指数函数,而泊松比为常数.采用Fourier变换和传递矩阵法将该混合边值问题转化为一对奇异积分方程,通过数值求解奇异积分方程组获得功能梯度材料涂层平面运动裂纹的应力强度因子.考察了结构几何尺寸、裂纹运动速度以及材料梯度参数对运动裂纹的应力强度因子的影响,发现材料梯度参数、结构几何尺寸、裂纹长度以及运动速度均对功能梯度材料动态断裂行为有显著影响.     

基于圆弧底试件的动态裂纹扩展及止裂规律研究

为了研究脆性材料的动态裂纹扩展及止裂规律，设计了一种带圆弧形底边的梯形开口边裂纹（trapezoidal opening crack with arc bottom，TOCAB）构型的试件。在落锤冲击设备加载下，对圆心角为0°、60°、90°和120°的TOCAB试件进行了冲击实验，并采用裂纹扩展计(crack propagation gauge，CPG)监测裂纹起裂和扩展时间，从而获得裂纹扩展速度。采用有限差分软件AUTODYN对落锤冲击设备和试件进行数值模拟，研究了裂纹的动态扩展过程及止裂规律。还基于实验和数值方法，计算了裂纹的临界动态应力强度因子。实验和数值结果均表明:3种弧度的TOCAB试件都可以实现运动裂纹止裂，该构型可用于研究动态裂纹止裂问题；数值计算的裂纹扩展路径与实验结果基本一致，验证了数值模型的有效性；裂纹起裂和止裂时刻的临界动态应力强度因子大于裂纹动态扩展过程中的临界动态应力强度因子。     

含节理岩体爆生裂纹扩展的动焦散模型实验研究

应用动态焦散线测试系统,模拟含节理岩体断裂爆破过程,进行了PMMA模型透射式动态焦散线实验,着重研究了爆炸初始裂纹与节理面不同夹角的情况下,裂纹尖端动态应力强度因子的变化规律,裂纹穿过节理面的扩展规律,以及裂纹扩展速度的变化规律。实验结果分析表明,爆生裂纹穿过节理面时,裂纹尖端的动态应力强度因子和裂纹扩展速度显著下降,穿过节理面后,应力强度因子又有所增强;裂纹穿过节理面时,裂纹会沿节理面偏离一段距离后沿初始裂纹方向继续扩展。研究结果可以为节理岩体的断裂爆破提供理论依据。     

爆炸载荷下板条边界斜裂纹的动态扩展行为

为了研究爆炸应力波作用下板条边界斜裂纹的动态扩展行为,首先分析了爆炸应力波在含边界斜裂纹板条中的传播,其次采用动态焦散线实验方法,进行了爆炸载荷下板条边界斜裂纹扩展规律的实验研究.研究结果表明,爆炸应力波作用下,板条试件边界斜裂纹的扩展过程中,裂纹扩展速度、扩展加速度和裂尖动态应力强度因子随时间波动变化,扩展速度最大值...     

含偏置裂纹三点弯曲梁的动态断裂行为研究

采用动态焦散线方法，对含偏置裂纹三点弯曲梁承受横向冲击的弯曲断裂行为进行了一系列动态断裂力学实验研究，分析了无量纲量ａ／ｌ的改变（ａ——初始裂纹偏离梁中心线的距离；ｌ——梁长度的一半）对于裂纹动态扩展行为（裂纹起始状态、裂纹尖端的复合应力强度因子、裂纹扩展速度、裂纹扩展轨迹）的影响，并借助动态光弹性应力分析，对应力波与扩展裂纹的相互作用以及应力波传播规律进行探讨．给出了裂纹尖端复合应力强度因子、裂纹扩展速度的变化、裂纹曲裂轨迹以及方向与梁中应力波传播的相互关系     

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

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

冲击载荷下扩展裂纹尖端动态能量释放率分布的焦散线分析

借助高速摄影捕捉裂纹瞬态扩展过程，利用动态焦散线研究了含有裂纹的三点弯曲梁在冲击载荷作用下扩展裂纹尖端的动态能量释放率分布规律；综合分析了裂纹扩展时间、长度、速度，以及扩展裂纹尖端动态应力强度因子与它的变化关系，表明了动态能量释放率在裂纹扩展过程中的驱动作用。     

粘弹性介质中扩展裂纹与J积分

本文提出了在线弹性及粘弹性介质中扩展裂纹与路径无关的J~*积分,并给出了严密的证明。文中证明了J~*积分与扩展裂纹尖端的张开位移(动态COD)之间有简单的关系,同时利用J~*积分求得了粘弹性介质中变速扩展裂纹尖端的奇异性。当裂纹以常速扩展时,J~*积分与能量释放率、动应力强度因子之间也有简单的关系。利用这些关系,我们给出了动态COD与动应力强度因子之间的关系式。     

Antiplane mechanical and inplane electric time-dependent load applied to two coplanar cracks in piezoelectric ceramic material

This work is concerned with the dynamic response of two coplanar cracks in a piezoelectric ceramic under antiplane mechanical and inplane electric time-dependent load. The cracks are assumed to act either as an insulator or as a conductor. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in Laplace transform domain. A numerical Laplace inversion algorithm is used to determine the dynamic stress and electric displacement factors that depend on time and geometry. A normalized equivalent parameter describing the ratio of the equivalent magnitude of electric load to that of mechanical load is introduced in the numerical computation of the dynamic stress intensity factor (DSIF) which has a similar trend as that for the pure elastic material. The results show that the dynamic electric field will impede or enhance crack propagation in a piezoelectric ceramic material at different stages of the dynamic electromechanical load. Moreover, the electromechanical response is greatly affected by the ratio of the crack length to the ligament between the cracks. The stress and electric displacement intensity factor can be combined by the energy density factor or function to address the fracture of piezoelectric materials under the combined influence of electromechanical loading.     

Griffith crack moving in a piezoelectric strip

Summary  The dynamic problem of an impermeable crack of constant length 2a propagating along a piezoelectric ceramic strip is considered under the action of uniform anti-plane shear stress and uniform electric field. The integral transform technique is employed to reduce the mixed-boundary-value problem to a singular integral equation. For the case of a crack moving in the mid-plane, explicit analytic expressions for the electroelastic field and the field intensity factors are obtained, while for an eccentric crack moving along a piezoelectric strip, numerical results are determined via the Lobatto–Chebyshev collocation method for solving a resulting singular integral equation. The results reveal that the electric-displacement intensity factor is independent of the crack velocity, while other field intensity factors depend on the crack velocity when referred to the moving coordinate system. If the crack velocity vanishes, the present results reduce to those for a stationary crack in a piezoelectric strip. In contrast to the results for a stationary crack, applied stress gives rise to a singular electric field and applied electric field results in a singular stress for a moving crack in a piezoelectric strip. Received 14 August 2001; accepted for publication 24 September 2002 The author is indebted to the AAM Reviewers for their helpful suggestions for improving this paper. The work was supported by the National Natural Science Foundation of China under Grant 70272043.     

Transient response of a cracked piezoelectric strip under arbitrary electro-mechanical impact

By using the well-developed integral transform methodology, the dynamic response of stress and electric displacement around a finite crack in an infinite piezoelectric strip are investigated under arbitrary dynamic anti-plane loads. The dynamic stress intensity factors and electric displacement are obtained analytically. It is shown that the dynamic crack-tip stress and electric field still have a square-root singularity. Numerical computations for the dynamic stress intensity factor show that the electric load has a significant influence on the dynamic response of stress field. The higher the ratio of the crack length to the width of the strip, the higher the peak value of the dynamic stress intensity factor is. On the other hand, the dynamic response of the electric field is determined solely by the applied electric load. The electric field will promote or retard the propagation of the crack depending on the time elapse since the application of the external electro-mechanical loads. The project supported by the National Natural Science Foundation of China and the Post-Doctor Science Foundation of China     

Multi-cracks problem for piezoelectric materials strip subjected to dynamic loading

This paper analyses and models the dynamic interaction among permeable multi-cracks in a piezoelectric strip under anti-plane shear waves by the Schmidt method. The Fourier transform is applied and then two pairs of triple integral equations can be solved using the Schmidt method. The results show that the stress and the electric displacement intensity factors of cracks depend on not only the crack length and the piezoelectric coefficient, but also the thickness of the piezoelectric strip, the distance between multi-cracks and the frequency of incident wave.     

Time-harmonic P-waves engulfing a rectangular limited-permeable crack in piezoelectric medium: Energy density and energy release

The dynamic behavior of a limited-permeable rectangular crack in a transversely isotropic piezoelectric material is impinged by to a P-wave. The generalized Almansi theorem and the Schmidt method are used to determine the stress intensity factor and energy density factor as the primary fracture criterion of failure. The mixed boundary value problem entails the evaluation of the appropriate crack edge stress singularities that are characteristics of the fundamental functions. The stress and electric displacement intensity factors are also used to find the energy release rate that can be computed numerically and compared with the results corresponding to those of the stress intensity factor, and energy density factor. Graphical presentation shows that the energy release rate is always negative for the boundary conditions considered while the energy density factors always remain positive. Under certain conditions, the stress and electric displacement intensity factors can be negative and subject to physical limitations. Piezoelectric material boundary value problem solutions should therefore be qualified by the application of failure criteria by fracture of otherwise, particularly when the mechanical and electrical energy can release by creating free surface at the macroscopic and microscopic scales. Negative energy release rate found for the piezoelectric medium in this work can be a case in point.Positive definiteness of the energy density factor can be applied to mutliscale fracture. This is not true for the stress intensity factor nor the energy release rate. Hence, crack initiation behavior for the permittivity of a rectangular crack due to the wave propagation effects may be studied. In particular, the initiation of micro-cracks may be identified with certain critical stress wave frequency band. Negative stress intensity factor may not enhance macrocracking but it does not exclude microcrack initiation.     

Effect of electric fields on fracture behavior of ferroelectric ceramics

The asymptotic problem of a semi-infinite crack perpendicular to the poling direction in a ferroelectric ceramic subjected to combined electric and mechanical loading is analyzed to investigate effect of electric fields on fracture behavior. Electromechanical coupling induced by the piezoelectric effect is neglected in this paper. The shape and size of the switching zone is shown to depend strongly on the relative magnitude between the applied electric field and stress field as well as on the ratio of the coercive electric field to the yield electric field. A universal relation between the crack tip stress intensity factor and the applied intensity factors of stress and electric field under small-scale conditions is obtained from the solution of the switching zone. It is found that the ratio of the coercive electric field to the yield electric field plays a significant role in determining the enhancement or reduction of the crack tip stress intensity factor. The fracture toughness variation of ferroelectrics under combined electric and mechanical loading is also discussed.     

功能梯度压电带中裂纹对SH波的散射

研究功能梯度压电带中裂纹对SH波的散射问题,为了便于分析,材料性质假定为指数模型,并假设裂纹面上的边界条件为电渗透型的.根据压电理论得到压电体的状态方程,利用Fourier积分变换,问题转化为对偶积分方程的求解.用Copson方法求解积分方程.求得了裂纹尖端动应力强度因子、电位移强度因子的解析表达式,最后数值结果显示了标准动应力强度因子与入射波数、材料参数、带宽、波数以及入射角之间的关系.     

The weight function for cracks in piezoelectrics

The weight function in fracture mechanics is the stress intensity factor at the tip of a crack in an elastic material due to a point load at an arbitrary location in the body containing the crack. For a piezoelectric material, this definition is extended to include the effect of point charges and the presence of an electric displacement intensity factor at the tip of the crack. Thus, the weight function permits the calculation of the crack tip intensity factors for an arbitrary distribution of applied loads and imposed electric charges. In this paper, the weight function for calculating the stress and electric displacement intensity factors for cracks in piezoelectric materials is formulated from Maxwell relationships among the energy release rate, the physical displacements and the electric potential as dependent variables and the applied loads and electric charges as independent variables. These Maxwell relationships arise as a result of an electric enthalpy for the body that can be formulated in terms of the applied loads and imposed electric charges. An electric enthalpy for a body containing an electrically impermeable crack can then be stated that accounts for the presence of loads and charges for a problem that has been solved previously plus the loads and charges associated with an unsolved problem for which the stress and electric displacement intensity factors are to be found. Differentiation of the electric enthalpy twice with respect to the applied loads (or imposed charges) and with respect to the crack length gives rise to Maxwell relationships for the derivative of the crack tip energy release rate with respect to the applied loads (or imposed charges) of the unsolved problem equal to the derivative of the physical displacements (or the electric potential) of the solved problem with respect to the crack length. The Irwin relationship for the crack tip energy release rate in terms of the crack tip intensity factors then allows the intensity factors for the unsolved problem to be formulated, thereby giving the desired weight function. The results are used to derive the weight function for an electrically impermeable Griffith crack in an infinite piezoelectric body, thereby giving the stress intensity factors and the electric displacement intensity factor due to a point load and a point charge anywhere in an infinite piezoelectric body. The use of the weight function to compute the electric displacement factor for an electrically permeable crack is then presented. Explicit results based on a previous analysis are given for a Griffith crack in an infinite body of PZT-5H poled orthogonally to the crack surfaces.     

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

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