功能梯度压电材料圆环位移和电势边值问题的精确解

圆环形的压电材料器件在智能结构中得到了广泛的应用。本文推导了横观各向同性功能梯度压电材料圆环在内、外边界上给定位移和电势情况下的一般解。极化方向在圆环的半径方向,材料常数的梯度方向也设定在半径方向,并可表示为半径r的幂,本构关系为线性。然后推导了压电圆环外壁固定、接地,内壁沿垂向有一微小位移、电势分别为余弦分布和均匀分布的问题的精确解,并计算了该问题在这两种电势情况下产生的无量纲形式的径向和环向位移、电势、应力及电位移沿径向分布的数值结果。计算中考虑了不同的材料梯度,以及内壁的位移与电势的不同比例。     

轴对称任意梯度圆筒的弹性动力学解

基于轴对称平面应变问题的运动方程及弹性梯度材料的应力和位移关系,通过将圆筒分层使材料性质离散为分段常数函数,同时在时域内应用有限差分格式,求得了材料性质沿径向梯度变化的圆筒弹性动力学解。本文解不仅适合任意梯度的弹性圆筒,而且容易满足多种形式的初始条件和边界条件。通过对材料性质沿径向为连续函数分布和分段函数分布的梯度圆筒数值分析,并与已有文献结果比较,得出本文解与已有文献的解吻合较好,验证了本文解的正确性和有效性。对材料性质为分段函数的三层组合圆筒分析发现,中间功能梯度层的指数分布因子对圆筒的径向位移和应力随时间变化都会产生显著影响。     

功能梯度压电、压磁柱对称问题的振动分析

本文在不考虑体力、体电流和体电荷的情况下, 假定压电、压磁柱壳的材料参数沿圆柱厚度方向呈幂函数分布时, 研究了径向载荷作用下功能梯度压电、压磁空心柱壳的空间柱对称径向振动问题. 首先在柱坐标系下, 由功能梯度材料的参数、本构、梯度和平衡方程推导得出外激励作用下以Bessel函数表示圆柱壳的应力、电势、磁势等物理量的稳态解, 进而对空间柱对称的压电、压磁功能梯度材料动力控制问题进行了理论分析. 可以看出, 当梯度参数时, 即完全退化为横观各项同性压电、压磁柱对称的振动问题, 与文献[16]的基本方程为柱坐标下得出的结果一致. 最后给出数值算例, 结果表明材料不均匀性对沿径向振动各物理量的显著影响, 且用一个特定不均匀性参数值可以优化电磁力耦合的性能, 这在现代工程设计中尤为重要.     

功能梯度压电、压磁柱对称问题的振动分析

本文在不考虑体力、体电流和体电荷的情况下, 假定压电、压磁柱壳的材料参数沿圆柱厚度方向呈幂函数分布时, 研究了径向载荷作用下功能梯度压电、压磁空心柱壳的空间柱对称径向振动问题. 首先在柱坐标系下, 由功能梯度材料的参数、本构、梯度和平衡方程推导得出外激励作用下以Bessel函数表示圆柱壳的应力、电势、磁势等物理量的稳态解, 进而对空间柱对称的压电、压磁功能梯度材料动力控制问题进行了理论分析. 可以看出, 当梯度参数时, 即完全退化为横观各项同性压电、压磁柱对称的振动问题, 与文献[16]的基本方程为柱坐标下得出的结果一致. 最后给出数值算例, 结果表明材料不均匀性对沿径向振动各物理量的显著影响, 且用一个特定不均匀性参数值可以优化电磁力耦合的性能, 这在现代工程设计中尤为重要.     

压电材料中切口/接头端部平面电弹性场奇异性有限元分析

为了对平面载荷作用下压电材料中切口或接头端部附近电弹性场奇异性问题进行分析，首先以应力平衡方程、Maxwell方程和和边界条件为基础，得到一种求解压电材料特征问题的弱式方程；其次，假定楔形切口或接头端部附近单元内位移和电势沿径向分布为指数形式，而周向方向分布则采用泡函数插值，将其代入弱式方程，建立一种只需对楔形切口或接头端部附近周边进行离散的一维简单有限元方法．压电材料的极化轴可以是任意方向．利用该有限元模型讨论了楔形切口角度、极化轴方向和边界条件对奇性场的影响．通过和其它特定情况下的现有解相比，证实了该文有限元数值方法的有效性，而且精度很高．     

正交各向异性功能梯度圆柱壳的应力分析及材料剪裁

假设正交各向异性功能梯度材料的弹性系数沿圆柱壳的径向按照任意连续函数变化,采用应力函数法和加权残值法导出了圆柱壳在非轴对称载荷作用下应力分析的一种新的数值解.建立了圆柱壳内部应力状态给定时材料剪裁问题的基本方程,提出了实现圆柱壳内部一种特殊的应力分布时所需要的材料弹性系数沿径向变化的解析解.通过数值算例验证了本文所导出的应力分析的数值解的正确性和收敛性,分析了弹性系数沿径向的变化对圆柱壳内部应力分布的影响.数值算例还给出了实现圆柱壳内部环向应力和切应力沿径向均匀分布时功能梯度材料的弹性系数沿径向的三种不同变化形式.所得研究结果可为正交各向异性功能梯度材料圆柱壳的设计提供一定的参考,同时材料剪裁的解析结果也可作为其他数值方法计算结果验证的考题.     

功能梯度压电、压磁柱壳问题的振动分析

本文在不考虑体力、体电流和体电荷的情况下,假定压电、压磁柱壳的材料参数沿圆柱厚度方向呈幂函数分布,研究了径向载荷作用下功能梯度压电、压磁空心柱壳的空间柱对称径向振动问题.首先在柱坐标系下,由功能梯度压电、压磁空心柱壳的参数、本构、梯度和运动方程推导得出外激励作用下以Bessel函数表示圆柱壳的应力、电势、磁势等物理量的稳态解,进而对空间柱对称的功能梯度压电、压磁柱壳的动力控制问题进行了理论分析.进一步可以看出,当梯度参数β=0时,即完全退化为横观各向同性压电、压磁柱对称的振动问题,与文献[20]的基本方程为柱坐标下得出的结果完全一致.最后给出数值算例,数值结果表明,材料不均匀性对沿径向振动各物理量有显著影响,且用一个特定不均匀性参数β值可以优化力电磁耦合的性能,这在现代工程设计中尤为重要.     

一种特殊梯度分布的功能梯度圆柱壳的二维分析

功能梯度材料是一种新型材料,其结构分析已成为当今力学研究的热点。本文对一种特殊梯度分布的功能梯度材料圆柱壳进行了二维精确分析。从弹性力学平面应变问题的基本方程出发,引入应力函数,导出功能梯度材料圆柱壳受静载作用下的控制微分方程。假设材料的杨氏模量沿半径方向呈幂函数分布,泊松比为常数,利用分离变量法,导出了简支边界情况下功能梯度圆柱壳的精确解。通过算例分析了不同梯度变化时,功能梯度圆柱壳内的应力和位移变化规律。计算结果表明不同梯度分布的圆柱壳结构中的应力、位移沿厚度方向的变化规律是不同的,有时甚至差别很大。因此对于材料性质梯度变化的功能梯度材料圆柱壳,必须针对其自身特点,建立相应的理论分析模型。     

基于混合数值法的功能梯度材料板瞬态热响应分析

采用混合数值法分析了功能梯度材料板中瞬态热响应问题．功能梯度材料板的材料特性沿板厚方向连续变化,在梯度板中心作用温度变化的热源,此热源沿板厚方向连续分布．得出了温度变化对梯度板的影响,观测点处的位移响应与应力变化随时间分布的曲线以及瞬态的位移响应．研究表明,温度变化对板的作用效果与外载作用对板的作用效果是相当的．此研究是对混合数值法的在热领域应用的延伸与开拓,对热响应问题的研究具有一定的指导意义．     

内含弹性介质功能梯度压电球壳的径向振动调控

针对内含弹性介质的功能梯度压电球壳的径向振动问题,利用变量替换技术,求得了解析解。该解适用于材料参数沿厚度按幂律变化且密度的梯度指标可与其它物理量的梯度指标不同的功能梯度压电球壳,克服了以往分析中所有材料参数梯度指标均假设为相等的局限。数值结果表明,功能梯度压电球壳的径向振动特性可利用材料参数梯度指标和内部弹性介质的刚度进行有效调控。     

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.     

Isolated crack in three-dimensional piezoelectric solid. Part II: Stress intensity factors for circular crack

This is part II of the work concerned with finding the stress intensity factors for a circular crack in a solid with piezoelectric behavior. The method of solution involves reducing the problem to a system of hypersingular integral equations by application of the unit concentrated displacement discontinuity and the unit concentrated electric potential discontinuity derived in part I [1]. The near crack border elastic displacement, electric potential, stress and electric displacement are obtained. Stress and electric displacement intensity factors can be expressed in terms of the displacement and the potential discontinuity on the crack surface. Analogy is established between the boundary integral equations for arbitrary shaped cracks in a piezoelectric and elastic medium such that once the stress intensity factors in the piezoelectric medium can be determined directly from that of the elastic medium. Results for the penny-shaped crack are obtained as an example.     

压电介质二维边界积分方程中的基本解

由于压电介质的变形－电场耦合效应及压电响应的各向异性，使解析求解压电介质问题的工作变量十分复杂，若采用边界元数值方法求解，必须具备积分方程中的基本解，本文根据电磁场方程及连续介质力学的耦合性质论层出了二维无限域中分别在单位力及单位电荷载作用下的位移场，电势场、应力场和电位移场的解，从而确立了边界积分方程中所必需的八个基本解。     

Numerical solution of polarization saturation/dielectric breakdown model in 2D finite piezoelectric media

The polarization saturation (PS) model [Gao, H., Barnett, D.M., 1996. An invariance property of local energy release rates in a strip saturation model of piezoelectric fracture. Int. J. Fract. 79, R25–R29; Gao, H., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491–510], and the dielectric breakdown (DB) model [Zhang, T.Y., Zhao, M.H., Cao, C.F., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311–327] explain very well some experimental observations of fracture of piezoelectric ceramics. In this paper, the nonlinear hybrid extended displacement discontinuity-fundamental solution method (NLHEDD-FSM) is presented for numerical analysis of both the PS and DB models of two-dimensional (2D) finite piezoelectric media under impermeable and semi-permeable electric boundary conditions. In this NLHEDD-FSM, the solution is expressed approximately by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with sources placed at chosen points outside the domain of the problem under consideration, and the extended Crouch fundamental solutions with extended displacement discontinuities placed on the crack and the electric yielding zone. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy certain conditions on the boundary of the domain, on the crack face and the electric yielding zone. The zero electric displacement intensity factor in the PS model or the zero electric field strength intensity factor in the DB model at the outer tips of the electric yielding zone is used as a supplementary condition to determine the size of the electric yielding zone. Iteration approaches are adopted in the NLHEDD-FSM. The electric yielding zone is determined, and the extended intensity factors and the local J-integral are calculated for center cracks in piezoelectric strips. The effects of finite domain size, saturation property and different electric boundary conditions, as well as different models on the electric yielding zone and the local J-integral, are studied.     

压电、压磁材料球对称问题的通解

研究压电、压磁材料在球坐标系下,不计体力、体电荷和体电流的情况下,由平衡方程、梯度方程、压电和压磁的本构方程导出应力、应变、位移、电位移、电场强度、电位势、磁感强度和磁位势各未知量的通解．考虑不同的边界条件,将其通解应用到应力、电学短路以及位移、电学开路和磁场分布的边界条件中,得到不同边界条件下问题的解。     

Modified strip saturated models for two equal collinear cracks with coalesced zones in piezoelectric media

Two equal collinear cracks with coalesced interior electric saturation zones are analytically studied for two-dimensional(2D) arbitrary polarized semipermeable piezoelectric media based on a modified strip saturated model. The strip saturated model is modified here by varying the strip saturated constant electric displacement condition to the polynomially varying electric displacement conditions. Based on the linear, quadratic,and cubic electric displacement conditions on the inner and outer saturated zones, different modified strip saturated models are proposed and studied for two equal collinear cracks. With the Stroh formalism and the complex variable technique, these fracture problems are reduced into different types of non-homogeneous Riemann Hilbert problems in unknown generalized complex potential functions. These mathematical problems are then solved with the Riemann-Hilbert approach to obtain the stress and electric displacement components at any point of the domain. The explicit expressions for the outer saturated zone length, the crack opening potential(COP), the crack opening displacement(COD), and the local intensity factors(LIFs) are derived. A numerical study is presented for the modified strip saturated model in 2D arbitrary polarized semipermeable PZT-4 material. The obtained results are compared with those of the strip saturated model, and the effects of the polynomially varying saturation condition on the saturated zones and the applied electrical loadings are presented.     

压电弹性层合板静力机耦合特性的解析解

应用经典权壳理论和压电理论对层合板结构进行简化。对四边简支压电层合板在不同电学边界条件下,包括四边电学开路和电学短路。上下表面受外电压及无外加电压作用,进行分析,求得了电势和挠度的解析表达式,给出了压电层和基体的挠度、电势分布图。     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

Dielectric Breakdown Model for an Electrically Semi-Permeable Penny-Shaped Crack in Three-Dimensional Piezoelectric Media

The dielectric breakdown(DB) model for a penny-shaped crack under a semipermeable boundary condition in a three-dimensional piezoelectric medium is studied.An approximate analytical solution is derived by using the boundary integral equation with extended displacement discontinuity,and the corresponding boundary element method with double iterative approaches is developed to analyze the semi-permeable crack.The effect of electric boundary conditions on crack faces is discussed on the basis of DB model.By comparing the DB model with the polarization saturation(PS) model for different piezoelectric materials,some interesting phenomena related to the electric yielding zone and local J-integral are observed.