An exact solution of crack problems in piezoelectric materials

IntroductionInrecentyearscrackproblemsinpiezoelectricmaterialhavereceivedmuchattention.Manytheoreticalanalyseshavebeengivenby[1～16].Itshouldbe,however,notedthatalltheaboveanalysesarebasedonaso-calledimpermeablecrackassumphon,i.e.thecrackfacesareassumedtobeimpermeabletoelectricfield,sotheelectricdisplacementvanishesinsidethecrack.Usingthisassumption,onewillobtainthefollowingresultS[2'3'5,6'9'16]=whentheelectricloadsaresolelyaPPliedatinLfinity,theelectricdisplacementissquare-rootsingularatthe…     

Rate effects on toughness in elastic nonlinear viscous solids

A micromechanics-based constitutive relation for void growth in a nonlinear viscous solid is proposed to study rate effects on fracture toughness. This relation is incorporated into a microporous strip of cell elements embedded in a computational model for crack growth. The microporous strip is surrounded by an elastic nonlinear viscous solid referred to as the background material. Under steady-state crack growth, two dissipative processes contribute to the macroscopic fracture toughness—the work of separation in the strip of cell elements and energy dissipation by inelastic deformation in the background material. As the crack velocity increases, voids grow in the strain-rate strengthened microporous strip, thereby elevating the work of separation. In contrast, the energy dissipation in the background material decreases as the crack velocity increases. In the regime where the work of separation dominates energy dissipation, toughness increases with crack velocity. In the regime where energy dissipation is dominant, toughness decreases with crack velocity. Computational simulations show that the two regimes can exist in certain range of crack velocities for a given material. The existence of these regimes is greatly influenced by the rate dependence of the void growth mechanism (and the initial void size) as well as that of the bulk material. This competition between the two dissipative processes produces a U-shaped toughness-crack velocity curve. Our computational simulations predict trends that agree with fracture toughness vs. crack velocity data reported in several experimental studies for glassy polymers and rubber-modified epoxies.     

Angle cracks in two bonded functionally graded piezoelectric material under anti-plane shear

Consider two bonded functionally graded piezoelectric material (FGPM) with finite height. Each material contains an arbitrary oriented crack. The material properties are assumed in exponential forms in the direction normal to the interface. The crack surface condition is assumed to be electrically impermeable or permeable. Using the Fourier transform technique, the problem can be reduced to a system of singular integral equations, which are then solved numerically by applying the Gauss-Chebyshev integration formula to obtain the stress intensity factors at the crack tips. Numerical calculations are carried out to obtain the energy density factor S and the energy release rate G. In impermeable case, the energy release rate has been shown to be negative as the electric loads are applied. The positive definite characteristic of the energy density factor makes it possible for predicting the fracture behavior of the cracked structure. The influences of the non-homogeneous parameters and crack orientation on the energy density factors at the crack tips are discussed in detail. The results show that the energy density factor at the crack tip will be increased when the crack tip is located within the softer material.     

Multiple cracks in thermoelectroelastic bimaterials

The formulation for thermal stress and electric displacement in an infinite thermopiezoelectric plate with an interface and multiple cracks is presented. Using Green's function approach and the principle of superposition, a system of singular integral equations for the unknown temperature discontinuity defined on each crack face is developed and solved numerically. The formulation can then be used to calculate some fracture parameters such as the stress–electric displacement and strain energy density factor. The direction of crack growth for many cracks in thermopiezoelectric bimaterials is predicted by way of the strain energy density theory. Numerical results for stress–electric displacement factors and crack growth direction at a particular crack tip in two crack system of bimaterials are presented to illustrate the application of the proposed formulation.     

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.      

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.     

Fracture of piezoelectric materials: energy density criterion

In this paper, the concept of energy density factor S for piezoelectric materials is presented. In addition to the mechanical energy the electrical energy is included as well. The direction of crack initiation is assumed to occur when S_min reaches a critical value S_cr that can be used as an intrinsic materials parameter and is independent of the crack geometry and loading. The result agrees with empirical evidence qualitatively and explains rationally the effect of applied electric field on fracture strength: positive electric fields decrease the apparent fracture toughness of piezoelectric materials while negative electric fields increase it.     

压电复合材料中正n边形孔边裂纹的反平面问题研究

本文基于Cauchy积分理论和Schwarz-Christoffel(SC)变换技术,针对压电复合材料中带一条裂纹的正n边形孔口缺陷的反平面断裂力学进行了探究.假设满足电不可通边界条件,利用Cauchy积分公式和留数定理,获得了任意正n边形裂尖处应力和电位移两个场强度因子以及全能量释放率的封闭形式的显式解.当正n边形边...     

Antiplane internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane

This paper deals with the antiplane magnetoelectroelastic problem of an internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane. The properties of the material such as elastic modulus, piezoelectric constant, dielectric constant, piezomagnetic coefficient, magnetoelectric coefficient and magnetic permeability are assumed in exponential forms and vary along the crack direction. Fourier transforms are used to reduce the impermeable and permeable crack problems to a system of singular integral equations, which is solved numerically by using the Gauss-Chebyshev integration technique. The stress, electric displacement and magnetic induction intensity factors at the crack tips are determined numerically. The energy density theory is applied to study the effects of nonhomogeneous material parameter β, edge conditions, location of the crack and load ratios on the fracture behavior of the internal crack.     

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.     

Anti-plane fracture analysis for the weak-discontinuous interface in a non-homogeneous piezoelectric bi-material structure

The concept of weak discontinuity is extended to functionally graded piezoelectric bi-material interface, and fracture analysis for the weak discontinuous interface is performed by the methods of Fourier integral transform and Cauchy singular integral equation. Numerical results of the total energy release rate (TERR) and the mechanical strain energy release rate (MSERR) are obtained to show the effects of non-homogeneity parameters, geometrical parameters and loads. Parametric studies yield three conclusions: (1) To reduce the weak-discontinuity of the interface is beneficial to resisting interfacial fracture. The effect of the weak-discontinuity of the interface on TERR and MSERR still depends on the strip width. The wider the strip, the more sensitive the TERR and MSERR will be to the weak-discontinuity of the interface. (2) To predict the effect of electric load on crack propagation, MSERR is more appropriate than TERR to be used as a fracture parameter. To predict the effect of mechanical load on crack propagation, both of them could be used as fracture parameters, and MSERR is more conservative. (3) Mechanical load and negative electric displacement load would promote crack propagation, but positive electric displacement load would retard it. For the structure applied by combined mechanical and positive electric displacement loads, crack propagation may be impeded by appropriately selecting the strip width and the ratio of non-homogeneity parameters.     

Exact solutions of two semi-infinite collinear cracks in piezoelectric strip

Using the complex variable function method and the conformal mapping technique,the fracture problem of two semi-infinite collinear cracks in a piezoelectric strip is studied under the anti-plane shear stress and the in-plane electric load on the partial crack surface.Analytic solutions of the field intensity factors and the mechanical strain energy release rate are derived under the assumption that the surfaces of the crack are electrically impermeable.The results can be reduced to the well-known solutio...     

压电材料中的Bueckner功共轭积分及和J积分M积分的关系

研究横观各向同性压电材料中裂纹问题,提出了Bueckner功共轭积分在这类材料中的表达式：并通过引出两类辅助的应力-位移-电位移-电势场,证明功共轭积分和这类材料中的J积分和M积分仍然存在简单的两倍关系由此,各类在脆性材料断裂问题中已广泛应用的权函数方法可顺理成章地推广到压电材料的研究中来．这对独立地确定电位移强度因子和经典的I、II型应力强度因子提供了有力的数学上的工具．进而通过计算机械应变能释放率对压电材料中裂纹的稳定做出判断．     

Anti-plane analysis of semi-infinite crack in piezoelectric strip

Using the complex variable function method and the technique of the conformal mapping, the fracture problem of a semi-infinite crack in a piezoelectric strip is studied under the anti-plane shear stress and the in-plane electric load. The analytic solutions of the field intensity factors and the mechanical strain energy release rate are presented under the assumption that the surface of the crack is electrically impermeable. When the height of the strip tends to infinity, the analytic solutions of an infinitely large piezoelectric solid with a semi-infinite crack are obtained. Moreover, the present results can be reduced to the well-known solutions for a purely elastic material in the absence of the electric loading. In addition, numerical examples are given to show the influences of the loaded crack length, the height of the strip, and the applied mechanical/electric loads on the mechanical strain energy release rate.     

The mode III cracks originating from the edge of a circular hole in a piezoelectric solid

A mode III fracture problem of edge cracks originating from a circular hole in an infinite piezoelectric solid is studied based on complex variable method combined with the method of conformal mapping. Explicit and exact expressions for the complex potentials, field intensity factors and energy release rates are presented under the assumption that the surface of the cracks and hole is electrically impermeable. Numerical analysis is then conducted to discuss the influences of crack length and applied mechanical/electric loads on the field intensity factors and energy release rate for one and two edge cracks, respectively. It is found that for the case of a single edge crack, the field intensity factors are greater than those of double edge cracks, and moreover the electric loads can either promote or retard crack growth, depending on the magnitude and direction of the applied electric loads.     

利用焦耳效应提高含裂纹金属构件抗裂性能问题的研究

设一无限大金属薄板中含有一个线裂纹,对金属板施加恒定的电流场,在两个裂尖处产生的热量远远大于其余地方产生的热量,可简化成两个点热源．经求解得到了问题的解析解,包括裂纹尖端附近区域温度、应力、应变、应变能密度因子的解析表达式．计算结果表明,裂纹尖端处的材料发生熔化而形成一个焊点,裂纹尖端明显纯化,可抑制裂纹的进一步扩展,提高含裂纹金属构件的抗裂性能．     

各向异性压电材料平面裂纹的耦合场分析

用Stroh方法分析了各向异性压电材料电导通型裂纹问题的耦合场。结果表明,裂纹面上的切向电场强度和法向电位移均为常数,在裂纹尖端有由弹性场的耦事作用产生的奇异电导通裂纹模型中的静电场对裂纹尖端扩展的能量释放率不作贡献。     

裂纹面摩擦接触引起的断裂韧性增长的研究

采用弹黏塑性的材料本构关系, 建立了压、剪混合型裂纹常速准静 态扩展的力学模型, 求得了裂纹面摩擦接触条件下裂纹尖端场的数值解, 并基于数 值结果讨论了扩展裂纹的摩擦效应. 计算和分析表明, 裂纹面的摩擦效应主要表现 在两个方面. 第一方面是摩擦会导致裂纹尖端区材料的断裂韧性增高, 并且裂纹面间的摩擦作用越强, 增韧效果越显著. 摩擦增韧的机制可以解释为裂纹 面间的摩擦作用导致裂纹尖端塑性区尺寸变大, 使裂纹尖端场的塑性变形能增加, 从而使得裂纹尖端区材料增韧. 摩擦生热并不是导致材料断裂韧性增长的根本机制. 第二方面是摩擦会导致``断裂延缓'. 利用裂纹面的摩擦来提高构件的承载能力和延长构件的服役寿命具有较大的工程实用价值.     

Energy release rate and path-independent integral in dynamic fracture of magneto-electro-thermo-elastic solids

The energy release rate and associated energy flux integral in dynamic fracture of magneto-electro-thermo-elastic solids are formulated with the inclusion of multi-field fully coupled effects based on fundamental principles of thermodynamics. The difference between the global and local dynamic contour integrals is caused by unsteady state, mechanical body force, electricity conduction and thermal effect as the closed contour including crack faces is chosen. This formulation successfully captures the crack-tip singularity of coupled fields, offers the right expression for the crack driving force, and resolves the controversial issue on magneto-electro-thermo-elastic fracture criterion. Especially, for steady-state crack propagation in a magneto-electro-elastic solid, the path-independent dynamic contour integral is determined from the asymptotic near-tip field solution based on the Stroh-type formalism and the resulting dynamic energy release rate has an odd dependence on the dynamic magnetic induction intensity factor and the dynamic electric displacement intensity factor.