GENERAL SOLUTION OF PLANE PROBLEM OF PIEZOELECTRIC MEDIA EXPRESSED BY “HARMONIC FUNCTIONS”

First, based on the basic equations of two-dimensional piezoelectroelasticity, a displacement function is introduced and the general solution is then derived. Utilizing the generalized Almansi's theorem, the general solution is so simplified that all physical quantities can be expressed by three harmonic functions. Second, solutions of problems of a wedge loaded by point forces and point charge at the apex are also obtained in the paper. These solutions can be degenerated to those of problems of point forces and point charge acting on the line boundary of a piezoelectric half-plane.Project supported by the National Natural Science Foundation of China     

双材料应力分析中的镜像点方法

提出一种分析各类双材料中任一点受集中力作用问题的方法.通过将结合界面或其自由表面看作镜面,将应力函数或位移函数设定成固定于受载点及其镜像点上的局部坐标系下的形式,利用界面连续条件和Dirichlet的单值性原理,所有应力函数或位移函数就可由无限体中集中力的解或半无限体表面集中力的解的应力函数求得.这种方法不仅可适用于单一界面的情况,也可使用于多个界面并存的情况,并且也可适用于具有自由表面的结合材料.这一方法可应用于各类结合材料、涂层薄膜材料、板材等.     

条带域内二维各向异性压电介质的Green函数

以压电各向异性弹性介质广义平面变形的Stroh一般解为基础，采用复变函数方法（即保角变换技术），研究了条带域介质内物理场的封闭形式解，求得了介质内某一点同时存在广义线位错和广义线力作用时的简单明确解，它就是边界元法中的Green函数，还分析了极化介质表面的电荷分布情况，并进而讨论了线电荷q与边界分布电荷间的库仑力问题，文中结果不仅适用于平面或反平面变形问题，而且也适用于两者耦合的二维变形问题。     

弹性力学中Fredholm积分方程组解法的表达通式及其讨论

本文采用覆盖域的概念,给出受外力作用的任意形状弹性体采用Fredhikm积分方程组解法的统一表达式,而且就含洞的无限大体且远场应力不为零,边界上有集中力作用等特殊情况进行讨论并给出相应的表达式,对边界上作用力和覆盖域边界上作用力之间在合力和合力矩上的关系也作了讨论,文中给出的算例表明,对于一些简单情况可以求得解析解。而对于需要采用数值解的问题,本文工作具有精度高、收敛快的优点。     

An electric point charge moving along the poling direction of a transversely isotropic piezoelectric solid

The problem of an electric point charge moving constantly along the poling direction of a transversely isotropic piezoelectric solid is considered in a moving coordinate system, which moves together with the electric point charge. A general solution in the moving coordinate system is given, and all the field components, such as displacements, electric potential, stresses and electric displacements, can be concisely expressed in terms of four quasi-harmonic functions. We also present two examples to demonstrate the effect of the moving velocity on the values of _i. Once the general solution is given, the axisymmetric problem of a moving electric point charge can be easily solved. The explicit expressions of all the field components caused by the moving electric charge are presented, and the effect of the moving velocity on these field components is numerically investigated.     

Elastic Displacement in a Half-Space Under the Action of a Tensor Force. General Solution for the Half-Space with Point Forces

The elastic displacement in an isotropic elastic half-space with free surface is calculated for a point tensor force which may arise from the seismic moment of seismic sources concentrated at an inner point of the half-space. The starting point of the calculation is the decomposition of the displacement by means of the Helmholtz potentials and a simplified version of the Grodskii-Neuber-Papkovitch procedure. The calculations are carried out by using generalized Poisson equations and in-plane Fourier transforms, which are convenient for treating boundary conditions. As a general result we compute the displacement in the isotropic elastic half-space with free surface caused by point forces with arbitrary structure and orientation, localized either beneath the surface (generalized Mindlin problem) or on the surface (generalized Boussinesq-Cerruti problems). The inverse Fourier transforms are carried out by means of Sommerfeld-type integrals. For forces buried in the half-space explicit results are given for the surface displacement, which may exhibit finite values at the origin, or at distances on the surface of the order of the depth of the source. The problem presented here may be viewed as an addition to the well-known static problems of elastic equilibrium of a half-space under the action of concentrated loads. The application of the method to similar problems and another approach to the starting point of the general solution are discussed.     

Study on the 3D Green's functions of the fluid and piezoelectric bimaterials

Because most piezoelectric devices have interfaces with fluid in engineering,it is valuable to study the coupled field between fluid and piezoelectric media.As the fundamental problem,the 3 D Green's functions for point forces and point charge loaded in the fluid and piezoelectric bimaterials are studied in this paper.Based on the 3 D general solutions expressed by harmonic functions,we constructed the suitable harmonic functions with undetermined constants at first.Then,the couple field in the fluid and piezoelectric bimaterials can be derived by substitution of harmonic functions into general solutions.These constants can be obtained by virtue of the compatibility,boundary,and equilibrium conditions.At last,the characteristics of the electromechanical coupled fields are shown by numerical results.     

On the general expression of Fredholm Integral Equations Method in elasticity

In this paper, the concept of covering domain is introduced to develop a general expression for the Fredholm Integral Equations Method, by which elasticity problems of arbitrarily shaped bodies loaded by external forces can be solved. Some special expressions are given for a body with non-zero remote stresses, or subjected to some concentrated forces on its boundary. The relationship between the loading forces and solutions are also discussed. Some analytical solutions can be obtained for simple cases. When numerical computations are needed for the solution, the method proves to have high precision and fast convergency.     

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

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

正交异性体某些断裂动力学问题

首先对具有任意自相似指数的正交异性弹性动力学问题推导了解的一般表示,给出了一般解法,然后用这一方法对若干具体问题进行求解,利用本文方法可以迅速将所论问题化为半平面上的Ｒｉｍａｎｎ－Ｈｉｌｂｅｒｔ问题,并可以相当简单地得到问题的闭合解。     

An antisymmetric problem of a penny-shaped crack in a piezoelectric medium

Summary  An exact, three-dimensional analysis is developed for a penny-shaped crack in an infinite transversely isotropic piezoelectric medium. The crack is assumed to be parallel to the plane of isotropy, with its faces subjected to a couple of concentrated normal forces and a couple of point electric charges that are antisymmetric with respect to the crack plane. The fundamental solution of a concentrated force and a point charge acting on the surface of a piezoelectric half-space is employed to derive the integral equations for the general boundary value problem. For the above antisymmetric crack problem, complete expressions for the elastoelectric field are obtained. A numerical calculation is finally performed to show the piezoelectric effect in piezoelectric materials. It is noted here that the present analysis is an extension of Fabrikant's theory for elasticity. Received 30 August 1999; accepted for publication 1 March 2000     

边界积分方程中近奇异积分计算的一种变量替换法

准确估计近奇异边界积分是边界元分析中一项很重要的课题,其重要性仅次于对奇异积分的处理. 近年来已发展了许多方法,都取得了一定程度的成功,但这个问题至今仍未得到彻 底的解决. 基于一种新的变量变换的思想和观点,提交了一种通用的积分变换法, 它非常有效地改善了被积函数的震荡特性,从而消除了积分的近奇异性,在不增加计算量的情况 下, 极大地改进了近奇异积分计算的精度. 数值算例表明,其算法稳定,效率高, 并可达到很高的计算精度,即使区域内点非常地靠近边界,仍可取得很理想的结果.     

The general solution of spherically isotropic magnetoelectroelastic media and its applications

A general solution of the three-dimensional equilibrium problem of spherically isotropic magnetoelectroelastic media is presented. Base on the obtained general solution, exact and compact form solutions are obtained for (1) a spherically isotropic magnetoelectroelastic cone subjected to concentrated force, concentrated couple, a point charge and a point electric current at its apex; (2) a spherically isotropic magnetoelectroelastic space with a concentrated force at the origin; (3) a spherical shell under spherically symmetric deformation; and (4) stress concentration around a spherical cavity subjected to remote uniform tensile force, electric charge and electric current.     

Bending of isotropic thin square plates by concentrated edge couples and forces

Summary In this paper the complex variable method of Muskhelishvili is applied to problems of bending for small deflections of thin, isotropic, homogeneous plates by concentrated edge couples and forces. The functional equation involved in Muskhelishvili's method is solved by using function theory. The necessary conformal mapping function is found from the Schwarz-Christoffel formula or expansion of elliptic functions. A general solution is given for plates which can be mapped on a unit circle by polynomial type mapping functions. Three particular problems are worked out in detail; namely those of approximately square plates subjected, respectively, to two bending couples, to two twisting couples, both applied at the ends of a diagonal, and to four forces applied at the four corners. Numerical results are presented in the form of tables and graphs.     

Two-dimensional piezoelectricity. Part II: general solution,Green’s function and interface cracks

The complete solution space of a piezoelectric material is the direct sum of several orthogonal eigenspaces, one for each distinct eigenvalue. Each one of the 14 different classes of piezoelectric materials has a distinct form of the general solution, expressed in terms of the eigenvectors of the zeroth and higher orders and a kernel matrix containing analytic functions. When these functions are chosen to be logarithmic, one obtains, in a unified way, Green’s function of the infinite space as a single 8 × 8 matrix function G_∞ for the various load cases of concentrated line forces, dislocations, and a line charge. This expression of Green’s function is valid for all classes of nondegenerate and degenerate materials. With an appropriate choice of the parameters, it reduces to the solution of a half space with concentrated (line) forces at a boundary point, and with dislocations in the displacements. As another application, eigenvalues and eigensolutions are obtained for the bimaterial interface crack problem.     

POLYNOMIAL SOLUTIONS TO PIEZOELECTRIC BEAMS(Ⅱ)--ANALYTICAL SOLUTIONS TO TYPICAL PROBLEMS

For the orthotropic piezoelectric plane problem, a series of piezoelectric beams is solved and the corresponding analytical solutions are obtained with the trialand-error method on the basis of the general solution in the case of three distinct eigenvalues, in which all displacements, electrical potential, stresses and electrical displacements are expressed by three displacement functions in terms of harmonic polynomials. These problems are cantilever beam with cross force and point charge at free end, cantilever beam and simply-supported beam subjected to uniform loads on the upper and lower surfaces, and cantilever beam subjected to linear electrical potential.     

POLYNOMIAL SOLUTIONS TO PIEZOELECTRIC BEAMS (Ⅱ)——ANALYTICAL SOLUTIONS TO TYPICAL PROBLEMS

IntroductionThis paper is a continuation of Ref.[1],in which a series of orthotropic piezoelectricplane problems was solved and the corresponding exact solutions were obtained with the trial-and-error method,on the basis of the general solution expressed …     

Circular crack in a transversely isotropic piezoelectric space under point forces and point charges

In this paper, two kinds of circular crack including external circular crack and penny-shaped crack in a transversely isotropic piezoelectric space are considered. Firstly, we obtain the solution to the problem of an external circular crack in a transversely isotropic piezoelectric space subjected to antisymmetric normal point forces and point charges. Based on this, the solution of one-sided loading of an external circular crack is constructed. Secondly, the real shape of an external circular crack and the opening displacement of a penny-shaped crack under an arbitrary point force and point charge are further obtained. At last, the results are presented in a graphical form. The project supported by the National Natural Science Foundation of China (19872060 and 69982009) and the Postdoctoral Foundation of China     

On the levitation of a charged body in an electrostatic field

The possibility of levitation of a body bearing a point electric charge in the electric field of a fixed point charge of the same sign is demonstrated. Stabilization of the unstable equilibrium state of the body in which the body weight is balanced by the Coulomb force is provided by gyro forces due to body rotation. The conservative stability of the levitating body corresponds to a finite range of angular velocities. It is shown that dissipative and circulatory forces considered together can improve the stability of the system to asymptotic stability. The attraction domain of the stable equilibrium position is studied numerically as a function of system parameters.     

PARALLEL COMPUTING FOR STATIC RESPONSE ANALYSIS OF STRUCTURES WITH UNCERTAIN-BUT-BOUNDED PARAMETERS

The vertex solution for estimation on the static displacement bounds of structures with uncertain-but-bounded parameters is studied in this paper. For the linear static problem, when there are uncertain interval parameters in the stiffness matrix and the vector of applied forces, the static response may be an interval. Based on the interval operations, the interval solution obtained by the vertex solution is more accurate and more credible than other methods （such as the perturbation method）. However, the vertex solution method by traditional serial computing usually needs large computational efforts, especially for large structures. In order to avoid its disadvantages of large calculation and much runtime, its parallel computing which can be used in large-scale computing is presented in this paper. Two kinds of parallel computing algorithms are proposed based on the vertex solution. The parallel computing will solve many interval problems which cannot be resolved by traditional interval analysis methods.