Effects of electric/magnetic impact on the transient fracture of interface crack in piezoelectric-piezomagnetic sandwich structure: anti-plane case

Due to the incompatibility of the interlaminar deformations,the interface debonding or cracking usually happens in a layered magnetoelectric(ME)structure under an applied load.In this paper,the transient responses of the anti-plane interface cracks in piezoelectric(PE)-piezomagnetic(PM)sandwich structures are studied by the standard methods of the integral transform and singular integral equation.Discussion on the numerical examples indicates that the PE-PM-PE structure under electric impact is more likely to fracture than the PM-PE-PM structure under a magnetic impact.The dynamic stress intensity factors(DSIFs)are more sensitive to the variation of the active layer thickness.The effects of the material constants on the DSIFs are dependent on the roles played by PE and PM media during the deformation process.     

A domain-independent interaction integral for magneto-electro-elastic materials

Magneto-electro-elastic (MEE) materials usually consist of piezoelectric (PE) and piezomagnetic (PM) phases. Between different constituent phases, there exist lots of interfaces with discontinuous MEE properties. Complex interface distribution brings a great difficulty to the fracture analysis of MEE materials since the present fracture mechanics methods can hardly solve the fracture parameters efficiently of a crack surrounded by complex interfaces. This paper develops a new domain formulation of the interaction integral for the computation of the fracture parameters including stress intensity factors (SIFs), electric displacement intensity factor (EDIF) and magnetic induction intensity factor (MIIF) for linear MEE materials. The formulation derived here does not involve any derivatives of material properties and moreover, it can be proved that an arbitrary interface in the integral domain does not affect the validity and the value of the interaction integral. Namely, the interaction integral is domain-independent for material interfaces and thus, its application does not require material parameters to be continuous. Due to this advantage, the interaction integral becomes an effective approach for extracting the fracture parameters of MEE materials with complex interfaces. Combined with the extended finite element method (XFEM), the interaction integral is employed to solve several representative problems to verify its accuracy and domain-independence. Good results show the effectiveness of the present method in the fracture analysis of MEE materials with continuous and discontinuous properties. Finally, the particulate MEE composites composed of PE and PM phases are considered and four schemes of different property-homogenization level are proposed for comparing their effectiveness.     

DECAY RATE OF SAINT-VENANT END EFFECTS FOR PLANE DEFORMATIONS OF PIEZOELECTRIC-PIEZOMAGNETIC SANDWICH STRUCTURES

This paper is concerned with the decay of Saint-Venant end effects for plane deformations of piezoelectric （PE）-piezomagnetic （PM） sandwich structures, where a PM layer is located between two PE layers with the same material properties or reversely. The end of the sandwich structure is subjected to a set of self-equilibrated magneto-electro-elastic loads. The upper and lower surfaces of the sandwich structure axe mechanically free, electrically open or shorted as well as magnetically open or shorted. Firstly the constitutive equations of PE mate- rials and PM materials for plane strain are given and normalized. Secondly, the simplified state space approach is employed to arrange the constitutive equations into differential equations in a matrix form. Finally, by using the transfer matrix method, the characteristic equations for eigen- values or decay rates axe derived. Based on the obtained characteristic equations, the decay rates for the PE-PM-PE and PM-PE-PM sandwich structures are calculated. The influences of the electromagnetic boundary conditions, material properties of PE layers and volume fraction on the decay rates are discussed in detail.     

Comprehensive investigation of stress intensity factors in rotating disks containing three-dimensional semi-elliptical cracks

Initiation and propagation of cracks in rotating disks may cause catastrophic failures. Therefore, determination of fracture parameters under different working conditions is an essential issue. In this paper, a comprehensive study of stress intensity factors(SIFs) in rotating disks containing three-dimensional(3D) semi-elliptical cracks subjected to different working conditions is carried out. The effects of mechanical properties, rotational velocity, and orientation of cracks on SIFs in rotating disks under centrifugal loading are investigated. Also, the effects of using composite patches to reduce SIFs in rotating disks are studied. The effects of patching design variables such as mechanical properties, thickness, and ply angle are investigated separately. The modeling and analytical procedure are verified in comparison with previously reported results in the literature.     

Propagation behavior of elastic waves in one-dimensional piezoelectric/piezomagnetic phononic crystal with line defects

Using a stiffness matrix method, we in- vestigate the propagation behaviors of elastic waves in one-dimensional （1D） piezoelectric/piezomagnetic （PE/PM） phononic crystals （PCs） with line defects by calculating energy reflection/transmittion coefficients of quasi-pressure and quasi-shear waves. Line defects are created by the re- placement of PE or PM constituent layer. The defect modes existing in the first gap are considered and the influences on defect modes of the material properties and volume fraction of the defect layers, the type of incident waves, the location of defect layer and the number of structural layers are discussed in detail. Numerical results indicate that defect modes are the most obvious when the defect layers are inserted in the middle of the perfect PCs; the types of incidence wave and material properties of the defect layers have important effects on the numbers, the location of frequencies and the peaks of defect modes, and the defect modes are strongly de- pendent on volume fraction of the defect layers. We hope this paper will be found useful for the design of PE/PM acoustic filters or acoustic transducer with PCs structures.     

On the influence of electric boundary conditions on dynamic SIFs in piezoelectric materials

Dynamic stress intensity factors (SIFs) for a straight crack in a piezoelectric material under time-harmonic L- and SH-wave loading are determined for different electric boundary conditions. Impermeable, permeable and limited permeable cracks are compared. The problem is formulated and numerically solved using a nonhypersingular traction-based boundary integral equation method where the fundamental solution is obtained by Radon transform. A parametric study in the frequency domain shows the dependence of the SIFs on the choice of the electrical boundary conditions at the crack faces.     

Strip-coalesced interior zone model for two unequal collinear cracks weakening piezoelectric media

In this paper, a mathematical strip-saturation model is proposed for a poled transversely isotropic piezoelectric plate weakened by two impermeable unequal-collinear hairline straight cracks. Remotely applied in-plane unidirectional electromechanical loads open the cracks in mode-I such that the saturation zone developed at the interior tips of cracks gets coalesced. The developed saturation zones are arrested by distributing over their rims in-plane normal cohesive electrical displacement. The problem is solved using the Stroh formalism and the complex variable technique. The expressions are derived for the stress intensity factors （SIFs）, the lengths of the saturation zones developed, the crack opening displacement （COD）, and the energy release rate. An illustrative numerical case study is presented for the poled PZT-5H ceramic to investigate the effect of prescribed electromechanical loads on parameters affecting crack arrest. Also, the effect of different lengths of cracks on the SIFs and the local energy release rate （LERR） has been studied. The results obtained are graphically presented and analyzed.     

Finite element analyses of cracks in piezoelectric structures: a survey

Piezoelectric materials have widespread applications in modern technical areas such as mechatronics, smart structures or microsystem technology, where they serve as sensors or actuators. For the assessment of strength and reliability of piezoelectric structures under combined electrical and mechanical loading, the existence of cracklike defects plays an important role. Meanwhile, piezoelectric fracture mechanics has been established quite well, but its application to realistic crack configurations and loading situations in piezoelectric structures requires the use of numerical techniques as finite element methods (FEM) or boundary element methods (BEM). The aim of this paper is to review the state of the art of FEM to compute the coupled electromechanical boundary value problem of cracks in 2D and 3D piezoelectric structures under static and dynamic loading. In order to calculate the relevant fracture parameters very precisely and efficiently, the numerical treatment must account for the singularity of the mechanical and electrical fields at crack tips. The following specialized techniques are presented in detail (1) special singular crack tip elements, (2) determination of intensity factors K _I –K _IV from near tip fields, (3) modified crack closure integral, (4) computation of the electromechanical J-integral, and (5) exploitation of interaction integrals. Special emphasis is devoted to a realistic modeling of the dielectric medium inside the crack, leading to specific electric crack face boundary conditions. The accuracy, efficiency, and applicability of these techniques are examined by various example problems and discussed with respect to their advantages and drawbacks for practical applications.     

Fracture behaviors of piezoelectric materials

Theoretical analyses and experimental observations of the failure and fracture behaviors of piezoelectric materials are presented. The theoretical analyses are based on the Stroh formalism. A strip dielectric breakdown model is proposed to estimate the effect of electrical non-linearity on the piezoelectric fracture of electrically insulated cracks. The reviewed experiments include the indentation fracture test, the bending test on smooth samples, the fracture test on pre-notched or pre-cracked samples, the environment-assisted fracture test, etc. For electrically insulated cracks, the experimental results show a complicated fracture behavior under combined electrical and mechanical loading. Fracture data are greatly scattered when a static electric field is applied. For electrically conducting cracks, the experimental results demonstrate that static electric fields can fracture poled and depoled lead zirconate titanate (PZT) ceramics. A charge-free zone model is introduced to understand the failure behavior of conducting cracks in the depoled lead zirconate titanate ceramics under electrical and/or mechanical loading. These theoretical and experimental results indicate that fracture mechanics concepts are useful in the study of the failure behaviors of piezoelectric materials.     

带裂纹的椭圆孔口问题的应力分析

断裂现象与材料和结构中的孔洞、缺口或裂纹等缺陷密切相关,这是因为缺陷附近的应力集中明显.该文利用复变方法,通过保角映射研究了带裂纹的椭圆孔洞的平面弹性问题,给出了应力强度因子的解析解.并由此计算了两互相垂直的裂纹问题.     

裂纹体分析的权函数理论与应用: 回顾和展望

断裂力学是工程材料和结构的疲劳与断裂分析、损伤容限设计和结构完整性评定的理论基础. 应力强度因子作为线弹性裂纹尖端奇异场的单一表征参量和裂纹扩展驱动力, 在裂纹体的断裂力学分析中发挥着关键作用. 权函数法为复杂受载裂纹体的应力强度因子求解计算提供了强有力的解析工具, 不但具有远高于各类数值解法的计算效率, 而且精度可靠, 使用方便. 本文结合笔者团队在权函数法方面的长期研究工作, 对该方法自20世纪70年代初提出至今半个世纪以来, 国际断裂界在二维和三维权函数理论与应用方面的主要研究进展作了回顾和评述, 并对其未来发展提出了展望. 主要内容涵盖: 当前国际断裂界广泛应用的3种二维裂纹解析权函数法简介和以格林函数为基准的验证评价; 三维裂纹问题的片条合成权函数法和点载荷权函数法; 权函数法在复杂受载裂纹体的应力强度因子和裂纹张开位移等关键力学参量计算、内聚力/桥连等裂纹模型分析、共线多裂纹权函数理论及其在剩余强度预测等方面的应用, 以及复杂裂纹几何的工程化权函数分析和权函数法的反向应用问题.      

面向压电智能结构精确变形的协同优化设计方法

智能结构集智能材料与传统材料于一体,能够实现结构的主动控制,在航空航天等领域具有巨大的应用潜力.由于其系统复杂且具有多场耦合效应,智能结构的整体式优化设计方法成为结构控制技术研究的关键之一.为了提高压电智能结构的整体性能和变形精度,提出了同时考虑压电驱动器布局(分布位置及角度)和基体结构拓扑构型的协同优化设计新方法.采用多点约束方法 (multi-point constraints,MPC)建立压电驱动器和基体结构的连接,定义一种与测量点目标位移相关的权重函数,以实现结构的精确变形控制.通过协同优化设计,压电驱动器可以获得最优的分布位置及角度,同时基体结构获得最优的拓扑构型,从而提升了压电智能结构系统的整体驱动性能和变形精度.通过进一步分析,研究了精确变形、体分比约束与结构优化构型和整体刚度的关系,以及优化结果中可能存在的传力路径畸变现象.数值算例的设计结果表明,采用协同优化设计方法,能够扩大结构的寻优空间,有效减小变形误差,实现压电智能结构的精确变形控制.     

压电智能结构拓扑优化研究进展

压电材料因其变形精度高、反应速度快、易于制作成小型化元件已经被广泛应用于精密驱动、振动控制、精确定位等领域。改变压电智能结构中压电元器件的位置、大小、形状等参数能够有效地改善系统的力学性能,因而吸引了许多学者和工程师的关注和研究。拓扑优化作为有效的优化工具,已经成功应用于压电智能结构的优化设计中。本文首先阐述了压电智能结构拓扑优化的背景和意义,简要回顾了压电智能结构主动控制及分析方法,并综述了面向结构静变形控制的压电智能结构优化、面向振动控制的压电智能结构优化、压电俘能器的设计与优化等三个方面的研究进展。最后,简单归纳压电智能结构拓扑优化研究中值得关注的几个问题。     

压电智能结构拓扑优化研究进展

压电材料因其变形精度高、反应速度快、易于制作成小型化元件已经被广泛应用于精密驱动、振动控制、精确定位等领域。改变压电智能结构中压电元器件的位置、大小、形状等参数能够有效地改善系统的力学性能,因而吸引了许多学者和工程师的关注和研究。拓扑优化作为有效的优化工具,已经成功应用于压电智能结构的优化设计中。本文首先阐述了压电智能结构拓扑优化的背景和意义,简要回顾了压电智能结构主动控制及分析方法,并综述了面向结构静变形控制的压电智能结构优化、面向振动控制的压电智能结构优化、压电俘能器的设计与优化等三个方面的研究进展。最后,简单归纳压电智能结构拓扑优化研究中值得关注的几个问题。     

An analytical solution for the stress field and stress intensity factor in an infinite plane containing an elliptical hole with two unequal aligned cracks

The existing analytical solutions are extended to obtain the stress fields and the stress intensity factors (SIFs) of two unequal aligned cracks emanating from an elliptical hole in an infinite isotropic plane. A conformal mapping is proposed and combined with the complex variable method. Due to some difficulties in the calculation of the stress function, the mapping function is approximated and simplified via the applications of the series expansion. To validate the obtained solution, several examples are analyzed with the proposed method, the finite element method, etc. In addition, the effects of the lengths of the cracks and the ratio of the semi-axes of the elliptical hole (a/b) on the SIFs are studied. The results show that the present analytical solution is applicable to the SIFs for small cracks.     

Stress intensity factor measurement of cracks using a piezoelectric element

We present a stress intensity factor (SIF) measurement method of cracks using a piezoelectric element and electrostatic voltmeter. In the method, an isotropic piezoelectric element is first adhered near the crack tip. Then, the surface electrodes are attached to the three different positions on the piezoelectric element. The electric potentials of the surface electrodes, which are proportional to the strain sum (ɛ_x+ɛ_y) on the structural member, are measured by an electrostatic voltmeter during load cycling. Mode I and mode II SIFs of the crack are estimated using the relationship between the SIF and (σ_x+σ_y). The applicability of the proposed method is examined through experiments and numerical analysis.     

智能板结构振动特性的半解析法

为智能板结构振动特性的分析提出了一种半解析法.根据压电材料修正后的Hellinger-Re issner变分原理,推导了压电材料的Ham ilton ian等参元.智能板结构的基体板和压电块(压电传感器或驱动器)被看作独立的三维体,并用Ham ilton ian等参元分别离散板和压电块.考虑到板和压电块在连接界面上广义应力和广义位移的连续性,联立板和压电块的方程得到整个结构的方程组.数值实例的分析结果证明了方法的可靠性.     

A domain-independent interaction integral for fracture analysis of nonhomogeneous piezoelectric materials

Piezoelectric materials and structures contain more or less electromechanical interfaces in engineering applications. It is difficult to obtain the fracture parameters efficiently of the piezoelectric materials with complex interfaces. This paper presents a domain-independent interaction integral for material nonhomogeneity and discontinuity which can be used for solving the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) of piezoelectric materials with complex interfaces efficiently. The interaction integral is based on the J-integral by superimposition of two admissible states and the present formulation does not involve any derivatives of mechanical and electric properties. Moreover, it is proved that the interface in the integral domain does not affect the value of the interaction integral and thus, the present method does not require electromechanical parameters of piezoelectric materials to be continuous. The interaction integral method combined with the extended finite element method (XFEM) is used to investigate the influences of material continuity on the SIF and the EDIF and the results show that the material parameters and their first-order derivatives affect both the SIF and the EDIF greatly, while the higher-order derivatives affect both of them slightly.     

压电智能结构模糊自学习控制(FSLC)

将模糊逻辑与学习控制的基本思想相结合,根据控制系统的动态输出特性,采用模糊控制对学习控制律中的参数进行实时校正,实现系统的动态学习过程,提出了一种适用于压电智能结构振动控制的模糊自学控制方法FSLC（FuzzySelf-LearningContr01）。分别采用三维8节点实体单元（Solid45）和耦合单元模拟主结构和压电致动器／传感器,基于ANSYS参数化语言编写了压电智能结构振动控制分析的有限元程序。通过数值仿真证明了模糊自学习控制方法能有效控制压电结构的振动,并提高了自学习控制的收敛速度和获得了很好的控制效果。     

Propagation of a mode-III interfacial crack in a piezoelectric–piezomagnetic bi-material

The transient response of a semi-infinite mode-III interfacial crack propagating between piezoelectric (PE) and piezomagnetic (PM) half spaces is investigated in this paper. The integral transform method together with the Wiener–Hopf and Cagniard–de Hoop techniques is used to solve the mixed boundary value problem under consideration. The existence of generalized Maerfeld–Tournois interfacial wave is discussed and the solutions of the coupled fields are derived for four different cases of bulk shear wave velocity. The dynamic intensity factors of stress, electric displacement and magnetic induction as well as energy release rate (ERR) are obtained in explicit forms. The numerical results of the universal functions and dimensionless ERR for several different material combinations are presented and discussed in details. It is found that the Bleustein–Gulyaev (generalized Maerfeld–Tournois) waves dominate the dynamic characteristics of the interfacial crack propagation in PE–PM bi-material.