数值求解不可压粘性流体定常运动的格林函数方法

本文提出了一种数值求解不可压粘性流体定常运动的格林函数方法.在本文中利用Stokes方程的基本解作为格林函数将求解不可压粘性流体定常运动的边值问题化为求解速度场和边界应力的非线性积分方程组,在解出速度场和边界应力后可直接计算流场中各点的压力;用有限元近似将积分方程离散化而进行数值求解。对于小雷诺数流动,只归结为求解边界积分方程,使求解区域减少一个维度。对于非线性问题,可用迭代方法求解,在每次迭代中只须解出边界点上的速度或应力。通过几个简单的算例,表明本文所提出的方法具有精度高、处理边界条件简单、通用性强的优点,并具有求解各种复杂流动的潜力。     

压电介质平面问题的一般解和基本解

本文从压电介质平面问题基本方程出发，得到了含体积力的基本方程组的一般解。对状态方程进行Ｆｏｕｒｉｅｒ变换，由一般解得到Ｆｏｕｒｉｅｒ变换下状态方程的通解，对于单位集中力和单位点电荷情形，给出了压电介质平面问题各种情况下的有限形式的基本解。     

压电介质平面问题的基本解

应用复变函数的方法,对于压电介质平面问题,分析导出了无限介质或半无限介质受任意 集中载荷作用时的复势函数基本解;这些结果可作为边界元法的基本解,以求解具有复杂边界压电体的平面问题。     

压电陶瓷中圆币形裂纹在横向剪力下的机—电耦合行为

以弹性位移分量和电热函数基本未知量时，横观各向同性压电介质三维问题的场方程可化为四个联立的二阶线性偏微分方程组，本文导出了用四个调和函数表示位移分量及电势函数的表达式，即得到了该场方程的势函数能通解，作为通解的应用举例，文中求解了圆币形裂纹受横向剪切载荷下圆币形裂纹的尖端场及应力、电位移强度因子均具有明显的机－电耦合性质，而应力和电位移分量在裂尖仍具有－１／２的奇异性。     

变厚度板弯曲的边界元分析法

由于变厚度板弯曲问题的控制分方程复杂，直接求解其基本解推导边界积分方程建立边界元分析法较为困难，本文通过引入等效荷载，等效刚度，将此问题的控制微分方程化成与普通薄板弯曲问题基本方程相同的形式，利用求解通板弯曲问题的边界元迭代求解，建立了分析变厚度板弯曲问题的蛤法，算例表明本方法理正确，精度良好。     

关于压电介质的守恒方程和路径无关积分

借鉴分析动力学中的Ｊａｃｏｂｉ积分和循环积分概念，以及电磁场理论中的能量矩概念，导出了压电介质在静态场中的守恒方程形式，由这些守恒方程即可得在位错，断裂力学和其他缺陷理论中应用广泛的路径无关积分。     

复变形式的各向异性板弯曲问题的基本解

提出了求解各向异性板弯曲问题基本解的新方法。得到的基本解简捷明了，相应的法向弯矩和相当剪力的表达式易求，故便于应用在一般边界条件的各向异性板弯曲问题的边界积分方程。     

层状层电介质空间轴对称问题的状态空间解

从横观各向同性压电介质空间轴对称问题的基本方程出发,建立了压电介质空间轴对称问题的状态变量方程,对状态变量方程进行Hankel变换,得到以状态变量表示的单层压电介质在Hankel变换空间中的解,讨论了3种不同特征根的情况,利用提出的解得到了半无限压电体在垂直集中载荷和点电荷作出下的Boussinesq解。利用传递矩阵方法导出了多层压电介质空间轴对称问题解一般解析式。     

二维边界元奇异积分和多域缩聚法分析

基于基本解的一种新的表达式，对二维边界元分析中奇异积分的精确求解进行了讨论，从几何方面对基本解的奇异性进行了分析，给出了超参非连续元离散位势和弹性力学问题边界积分方程时奇异积分计算的精确式，从而为判断各种近似方法的优劣和间接方法的精度提供了依据，也为精确地分析了大规模问题提供了一条有效的途径。     

非线性大变形问题边界元分析的自校正方法

针对用增量法求解非线性方程解的漂移问题，在非线性问题边界元法计算中建立了自我校正方法，对在拖带坐标上建立的增量形式的基本方程，引入Ｌａｎｇｒａｎｇｅ校正因子，以全量形式的基本方程作为其辅助方程，在此基础上导出含校正项的边界积分方程，边界元自我校正方法的建立有效地保证了在非线性问题的计算中最终收敛在其解附近，提高了计算精度和运算效率。     

Isolated crack in three-dimensional piezoelectric solid: Part I - Solution by Hankel transform

The fundamental solutions are obtained for a unit concentrated electric potential discontinuity and unit concentrated displacement discontinuity in a three-dimensional piezoelectric medium. Displacements and stresses are derived by application of the boundary integral equation method. These expressions are used to obtain the stress intensity factors for a circular crack in Part II of the study.     

Penny-shaped crack in transversely isotropic piezoelectric materials

Using a method of potential functions introduced successively to integrate the field equations of three-dimensional problems for transversely isotropic piezoelectric materials, we obtain the so-called general solution in which the displacement components and electric potential functions are represented by a singular function satisfying some special partial differential equations of 6th order. In order to analyse the mechanical-electric coupling behaviour of penny-shaped crack for above materials, another form of the general solution is obtained under cylindrical coordinate system by introducing three quasi-harmonic functions into the general equations obtained above. It is shown that both the two forms of the general solutions are complete. Furthermore, the mechanical-electric coupling behaviour of penny-shaped crack in transversely isotropic piezoelectric media is analysed under axisymmetric tensile loading case, and the crack-tip stress field and electric displacement field are obtained. The results show that the stress and the electric displacement components near the crack tip have (r ^−1/2) singularity. The project supported by the Natural Science Foundation of Shaanxi Province, China     

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.     

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

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

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…     

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.     

RESPONSE ANALYSIS OF PIEZOELECTRIC SHELLS IN PLANE STRAIN UNDER RANDOM EXCITATIONS

The random response of a piezoelectric thick shell in plane strain state under boundary random excitations is studied and illustrated with a piezoelectric cylindrical shell. The differential equation for electric potential is integrated radially to obtain the electric potential as a function of displacement. The random stress boundary conditions are converted into homogeneous ones by transformation,which yields the electrical and mechanical coupling differential equation for displacement under random excitations. Then this partial differential equation is converted into ordinary differential equations using the Galerkin method and the Legendre polynomials,which represent a random multi-degree-of-freedom system with asymmetric stiffness matrix due to the electrical and mechanical coupling and the transformed boundary conditions. The frequency-response function matrix and response power spectral density matrix of the system are derived based on the theory of random vibration. The mean-square displacement and electric potential of the piezoelectric shell are finally obtained,and the frequency-response characteristics and the electrical and mechanical coupling properties are explored.     

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.     

压电材料平面应力状态的直线裂纹问题一般解

研究了含直线裂纹系的压电材料平面应力问题单个裂纹和双裂纹问题的封闭解答表明，在裂纹尖端，应力、电场强度和电位移有１／２阶的奇异性并与前人结果比较了产生电场奇异性的物理因素     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.