Finite element analysis of thermo-elastic behavior of piezoelectric structures under finite deformation

It is noted that the behavior of most piezoelectric materials is temperature dependent and such piezo-thermo-elastic coupling phenomenon has become even more pronounced in the case of finite deformation. On the other hand, for the purpose of precise shape and vibration control of piezoelectric smart structures, their deformation under external excitation must be ideally modeled. This demands a thorough study of the coupled piezo-thermo-elastic response under finite deformation. In this study, the governing equations of piezoelectric structures are formulated through the theory of virtual displacement principle and a finite element method is developed. It should be emphasized that in the finite element method the fully coupled piezo-thermo-elastic behavior and the geometric non-linearity are considered. The method developed is then applied to simulate the dynamic and steady response of a clamped plate to heat flux acting on one side of the plate to mimic the behavior of a battery plate of satelite irradiated under the sun. The results obtained are compared against classical solutions, whereby the thermal conductivity is assumed to be independent of deformation. It is found that the full-coupled theory predicts less transient response of the temperature compared to the classic analysis. In the steady state limit, the predicted temperature distribution within the plate for small heat flux is almost the same for both analyses. However, it is noted that increasing the heat flux will increase the deviation between the predictions of the temperature distribution by the full coupled theory and by the classic analysis. It is concluded from the present study that, in order to precisely predict the deformation of smart structures, the piezo-thermo-elastic coupling, geometric non-linearity and the deformation dependent thermal conductivity should be taken into account. Project supported by the National Natural Science Foundation of China (Nos. 10132010 and 50135030) and the Foundation of In-service Doctors of Xi'an Jiaotong University.     

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.     

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

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

Periodic cracks in a functionally graded piezoelectric layer bonded to a piezoelectric half-plane

Studied is the problem of a periodic array of cracks in a functionally graded piezoelectric strip bonded to a homogeneous piezoelectric material. The properties of the functionally graded piezoelectric strip, such as elastic modulus, piezoelectric constant and dielectric constant, are assumed in exponential forms and vary along the crack direction. The crack surface condition is assumed to be electrically impermeable or permeable. Integral transform and dislocation density functions are employed to reduce the problem to the solution of a system of singular integral equations. The effects of the periodic crack spacing, material constants and the geometry parameters on the stress intensity factor, the energy release ratio and the energy density factor are studied.     

Asymptotic fields in the finite element analysis of electrically permeable interface cracks in piezoelectric bimaterials

Summary The interface crack problem for a piezoelectric bimaterial based on permeable conditions is studied numerically. To find the singular electromechanical field at the crack tip, an asymptotic solution is derived in connection with the conventional finite element method. For mechanical and electrical loads, the complex stress intensity factor for an interface crack is obtained. The influence of the applied loads on the electromechanical fields near the crack tip is also studied. For a particular case of a short crack with respect to the bimaterial size, the numerical results are compared with the exact analytical solutions, obtained for a piezoelectric bimaterial plane with an interface crack.One author (V.G.) gratefully acknowledges the support provided by the Alexander von Humboldt Foundation of Germany.accepted for publication 7 June 2004     

Finite element based form-finding algorithm for tensegrity structures

This paper presents a novel form-finding algorithm for tensegrity structures that is based on the finite element method. The required data for the form-finding is the topology of the structure, undeformed bar lengths, total cable length, prestress of cables and stiffness of bars. The form-finding is done by modifying the single cable lengths such that the total cable length is preserved and the potential energy of the system is minimized. Two- and three-dimensional examples are presented that demonstrate the excellent performance of the proposed algorithm.     

混凝土裂缝开展宽度计算的有限单元法

本文采用“最大拉应力单元开裂”准则，改进了非线性分析的迭代流程，使其最大可能地模拟实际结构的开裂情况，进而可以基于分布裂缝模式计算裂缝发生位置，裂缝间距，裂缝开展深度和宽度。算例表明，该法所计算的裂缝分布规律是正确的，计算裂缝宽度与试验实测值及规范公式计算值基本吻合，可以推广应用。     

热压电桁架结构屈曲稳定性有限元分析

梁、板、壳组合结构在航空航天和控制系统中得到了广泛应用。一般来说,桁架结构是作为整体结构的支撑来使用的,综合考虑热、机械和电场作用有助于充分发挥材料性能。本文从力的平衡方程出发,考虑温度场、机械场和电场多场耦合的情况,推导出了热压电桁架结构的静力分析有限元模型,并在此基础上给出了在多场荷载共同作用下的结构屈曲稳定性求解的有限元方程。数值算例表明,多场耦合效应能够对结构的屈曲稳定性影响较大。合理利用温度的影响,可以提高结构工作性能,并为梁、板、壳组合结构的优化和控制提供了理论基础。     

Interacting multiple cracks in piezoelectric materials

In this paper, a method is presented to calculate the plane electro–elastic fields in piezoelectric materials with multiple cracks. The cracks may be distributed randomly in locations, orientations and sizes. In the method, each crack is treated as a continuous distributed dislocations with the density function to be determined according to the conditions of external loads and crack surfaces. Some numerical examples are given to show the interacting effect among multiple cracks.     

Analysis of multiple axisymmetric annular cracks in a piezoelectric medium

An axisymmetric annular electric dislocation is defined. The solution of axisymmetric electric and Volterra climb and glide dislocations in an infinite transversely isotropic piezoelectric domain is obtained by means of Hankel transforms. The distributed dislocation technique is used to construct integral equations for a system of co-axial annular cracks with so-called permeable and impermeable electric boundary conditions on the crack faces where the domain is under axisymmetric electromechanical loading. These equations are solved numerically to obtain dislocation densities on the crack surfaces. The dislocation densities are employed to determine field intensity factors for a system of interacting annular and/or penny-shaped cracks.     

含压电片层合壳的有限元分析与控制仿真

本文首先针对含压电材料的一般结构，推导了按位移和电位移求解时的混合变分原理，在此基础上，通过对面电荷取变分，直接得到了含按电压驱动压电片层合壳的有限元方程，最后给出了利用压电片进行静变形和振动控制的仿真算例，并进行了分析。结果表明，本文建立的有限元方程准确、可靠。     

具有压电材料薄板弯曲控制的有限元法

考虑一均匀各向同性薄板结构,在其上下表面离散分布状压电执行元件,在横向外载和电压共同作用下,分析板的弯曲变形,并通过变化作用于执行元件上的电压对板的变形进行控制,根据逆压电效应将电压转换成作用于板上的等效作用量,用Ｈａｍｉｌｔｏｎ原理导出压电结构板弯曲变形的有限元公式,并对环形板情况进行数值计算。     

Finite element formulation for dynamics of delaminated composite beams with piezoelectric actuators

As a part of an effort to develop a model-supported method for detection of delaminations in composite beams with the use of time responses to external excitations, a finite element formulation for dynamics of a composite beam with delamination and attached piezoelectric actuators is developed. In this formulation account is taken of transverse shear deformation and nonlinear through-thickness variation of the longitudinal displacement. Parameters that characterize the delamination are incorporated into the formulation that makes the finite element model convenient for use in conjunction with damage identification (not discussed in the present paper). Computational predictions of frequencies show good agreement with experimental results.     

预应力混凝土平面杆系结构的有限元方法

建立了基于有限元方法的考虑材料和几何非线性的任意截面预应力混凝土平面杆系结构的数值分析模型,可用于模拟预应力混凝土大跨度梁、单向偏压细长柱等的非线性全过程结构响应。引入修正的Rodriguez截面模型确定截面切线刚度,其中混凝土的贡献通过截面边界顶点定义的梯形单元来实现;在此基础上利用传统的平面非线性杆单元导出了标准有限元公式。通过两个算例验证了该模型的可靠性和适用性。     

Finite element formulation for active functionally graded thin-walled structures

For the analysis and design process of smart structures with integrated piezoelectric patches, the finite element method provides an effective simulation approach. In this paper, an attempt on modeling and simulation of the behavior of hybrid active structures is carried out using developed Kirchhoff-type-four-node shell element.The finite element results are compared with reference solutions taking into account the electromechanical responses of smart structures with various geometries, and the results show very high agreement. The main aspect of the application of the proposed element is to predict the behavior of FGM shells containing piezoelectric layers. A set of numerical analyses is performed in order to highlight the applicability and effectiveness of the present finite element model, notably for smart FGM structures. A comprehensive parametric study is conducted to show the influence of material composition, the placement and the thickness of the piezoelectric layers on the deformation of the laminated structure.     

Closed-form solution for two collinear cracks in a piezoelectric strip

Under the permeable electric boundary condition, the problem of two collinear anti-plane shear cracks lying at the mid-plane of a piezoelectric strip is investigated. By using the Fourier transform, the associated problem is reduced to a singular integral equation. Solving the resulting equation analytically, the electro-elastic field intensity factors and energy release rates at the inner and outer crack tips can be determined in explicit form. Numerical results for PZT-5H piezoelectric ceramic are also presented graphically. The results reveal that the effect of electric field on crack growth in piezoelectric materials is dependent on applied elastic displacement.     

传统蒙古包结构在地震作用下的有限元分析

利用Solidworks及ABAQUS软件建立了传统蒙古包的有限元模型,探讨了传统蒙古包在地震作用下的受力机理。结果表明:在设防烈度为8度、最大基本加速度为0.30g的地震作用下,结构构件并未出现破坏,结构整体性较好;传统蒙古包结构上部为主要的耗能部位,乌尼与上部的套脑通过沿地震波方向的晃动释放能量;结构的易破坏处为乌尼杆、里围绳以及哈那与里围绳连接处、乌尼与套脑连接处。     

INTERACTION OF ELECTRIC CHARGES IN A PIEZOELECTRIC WITH RIGID EXTERNAL CRACKS

IntroductionRecently much attention has been paid to Green’s functions[1 -7]for anisotropic mediabecause of their important applications in boundary element methods used widely inengineering analysis. Generally speaking, it is relatively easy to obtain G…     

Out-of-plane interface cracks in dissimilar piezoelectric materials

Summary This paper investigates the problem of an anti-plane interfacial crack between two dissimilar piezoelectric material layers. A single crack is first considered. The effect of interaction of two collinear cracks in the medium on the field intensity factors is investigated. The solutions of several particular cases, including an infinite piezoelectric bi-material and a piezoelectric material bonded to an elastic medium, are given. The bi-material constants governing the behavior of the crack tip fields are identified. By considering the crack as a notch of finite thickness, it is shown that the thickness of the notch has a pronounced influence on the crack tip field. The results for the assumption of a permeable crack represent the limit case where the notch thickness is reduced to zero.BLW would like to thank the National Science Foundation of China (#10102004) and the City University of Hong Kong (DAG #7100219) for the support of this work. YGS also thanks the Multidiscipline Scientific Research Foundation Project (HIT. MD 2001. 39) of the Harbin Institute of Technology and the SRF for ROCS, SEM.accepted for publication 3 April 2003