A boundary integral equation formulation for the Reissner's plates resting on two-parameter foundation

In this paper a boundary integral equation formulation for the Reissner's plates resting on a two-parameter foundation is established. With the aid of the Hormander Operator method, the equations of the corresponding fundamental solutions are converted into a sixth order partial differential equation with a scale function as an unknown. In order to reduce the equation further, two auxiliary functions are introduced. They satisfy a second and a fourth order equation respectively. The expressions of the auxiliary functions can be derived easily. The fundamental solutions of the Reissnei's plates on the two-parameter foundation arc expressed by a linear combination of the auxiliary functions and their derivatives. The boundary integral equations are formulated by the use of the weighted residual procedure. The fundamental solutions obtained are taken as the kernel functions of the boundary integral equations. A few examples are studied. The numerical results show high accuacy and efficiency of the present formulation.This work was supported by the National Natural Science Foundation of China.     

An asymptotic analysis of composite beams with kinematically corrected end effects

A finite element-based beam analysis for anisotropic beams with arbitrary-shaped cross-sections is developed with the aid of a formal asymptotic expansion method. From the equilibrium equations of the linear three-dimensional (3D) elasticity, a set of the microscopic 2D and macroscopic 1D equations are systematically derived by introducing the virtual work concept. Displacements at each order are split into two parts, such as fundamental and warping solutions. First we seek the warping solutions via the microscopic 2D cross-sectional analyses that will be smeared into the macroscopic 1D beam equations. The variations of fundamental solutions enable us to formulate the macroscopic 1D beam problems. By introducing the orthogonality of asymptotic displacements to six beam fundamental solutions, the end effects of a clamped boundary are kinematically corrected without applying the sophisticated decay analysis method. The boundary conditions obtained herein are applied to composite beams with solid and thin-walled cross-sections in order to demonstrate the efficiency and accuracy of the formal asymptotic method-based beam analysis (FAMBA) presented in this paper. The numerical results are compared to those reported in literature as well as 3D FEM solutions.     

Closed form solutions for free vibrations of rectangular Mindlin plates

A new two-eigenfunctions theory, using the amplitude deflection and the generalized curvature as two fundamental eigenfunctions, is proposed for the free vibration solutions of a rectangular Mindlin plate. The three classical eigenvalue differential equations of a Mindlin plate are reformulated to arrive at two new eigenvalue differential equations for the proposed theory. The closed form eigensolutions, which are solved from the two differential equations by means of the method of separation of variables are identical with those via Kirchhoff plate theory for thin plate, and can be employed to predict frequencies for any combinations of simply supported and clamped edge conditions. The free edges can also be dealt with if the other pair of opposite edges are simply supported. Some of the solutions were not available before. The frequency parameters agree closely with the available ones through pb-2 Rayleigh-Ritz method for different aspect ratios and relative thickness of plate.     

On the half-infinite crack problem in thermo-electro-elasticity

This paper concerns itself with the fundamental solutions of the thermo-electro-elastic field in an infinite medium, weakened by a half-infinite plane crack with two identical point thermal loads applied on the crack surfaces. The corresponding mixed boundary value problem is solved by virtue of the potential theory method conjugated with the general solutions. The boundary governing equations are solved by using the results available in literature. Exact and complete three-dimensional (3D) fundamental solutions are presented in terms of elementary functions. The singularity at the crack tip is analyzed explicitly. The obtained solutions will be of high significance to the related BEM analysis.     

Variational principles for hydrodynamic impact problems

We first establish the rigorous field equations of the two continuous stages before and after entering water. Then correspondently, we obtain the specific variational principles, bounded theorems, and boundary integral equations of the second stage problems. The existence of solutions are proved and the scheme of solving the solutions are provided. Finally, as a numerical example, the ship's wave resistence problem is used to demonstrate the specific application of the second stage problems and its accuracy. Then we provide a rigorous and sound theoretical basis of variational finite element method and boundary element method for calculating the accurately fundamental equations.     

正交各向异性厚板的边界元解法

本文利用 Hormander 算子法和平面波分解法导出了计入剪切变形的正交各向异性厚板的基本解。建立了计入剪切变形的正交各向异性厚板的边界积分方程。文中详细地讨论了基本解的数值计算,并用边界元法分析了一些算例。     

Simulation of particle/fluid flows in vertical circular pipes

A two fluid continuum model is applied to the simulation of steady fully developed particle/fluid flow in a vertical circular pipe. Both closed form and numerical solutions are obtained to the associated governing equations. These solutions are compared to the predictions of another two fluid continuum model. These comparisons are used to illustrate the differences in the predictions of two widely used models, even in the fundamental problem of steady flow in a circular pipe.     

Boundary integral equations of unique solutions in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＢＥＭ）ｐｒｏｖｉｄｅｓａｎａｔｔｒａｃｔｉｖｅａｌｔｅｒｎａｔｉｖｅｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｉｔｓｍａｉｎａｄｖａｎｔａｇｅｓａｒｅｅｃｏｎｏｍｉｃａｌａｎｄｐａｒｔｉｃｕｌａｒｌｙｃｏｎｖｅｎｉｅｎｔｆｏｒｕｎｂｏｕｎｄｅｄｄｏｍａｉｎａｎｄｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｐｒｏｂｌｅｍｓ．Ｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎ（ＢＩＥ）ｉｓｔｈｅ…     

Mode I fracture analysis of a piezoelectric strip based on real fundamental solutions

In existing papers, mode I crack problems of piezoelectric ceramics are generally solved in complex domain because of the complex fundamental solutions of in-plane piezoelectric governing equations. In fact, these problems can alternatively be analyzed in real number field by recasting the solutions in real form instead. The main purpose of the present work is to develop such real fundamental solutions by detailed eigenvalue and eigenvector analyses. As an example of application, the widely studied fracture problem of a piezoelectric strip with a center-situated crack under mode I loading condition is then revisited based on the real fundamental solutions. Mixed boundary value conditions of the crack are transformed into Cauchy singular integral equations, which are then solved numerically to get fracture parameters including the energy release rate and intensity factors. Convergence behaviors of the kernel functions are surveyed. Theoretical derivation and computation are validated by the exact solution in a special case. The effect of a combined geometrical parameter on the crack is also investigated.     

Fundamental solutions of the streamfunction–vorticity formulation of the Navier–Stokes equations

A new boundary element procedure is developed for the solution of the streamfunction–vorticity formulation of the Navier–Stokes equations in two dimensions. The differential equations are stated in their transient version and then discretized via finite differences with respect to time. In this discretization, the non-linear inertial terms are evaluated in a previous time step, thus making the scheme explicit with respect to them. In the resulting discretized equations, fundamental solutions that take into account the coupling between the equations are developed by treating the non-linear terms as in homogeneities. The resulting boundary integral equations are solved by the regular boundary element method, in which the singular points are placed outside the solution domain.     

Elastic Fields for Point and Partial Line Loading in Transversely Isotropic Linear Elasticity

In this paper some fundamental concentrated loading solutions are derived for a transversely isotropic full space. As a starting point the potential function representation for the elastic field is re-examined in light of a recent result derived by the author. It is shown that expressions for two of the stress components need to be modified from what is given in some of the existing literature. The use of these new expressions is first demonstrated by considering two point loading cases. Subsequent analysis integrates these two point force solutions over finite line segments to obtain solutions for various cases of partial line loading. The ramifications of the two modified stress equations on the partial line loading solutions are also discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

On the fundamental equations of electromagnetoelastic media in variational form with an application to shell/laminae equations

The fundamental equations of elasticity with extensions to electromagnetic effects are expressed in differential form for a regular region of materials, and the uniqueness of solutions is examined. Alternatively, the fundamental equations are stated as the Euler–Lagrange equations of a unified variational principle, which operates on all the field variables. The variational principle is deduced from a general principle of physics by modifying it through an involutory transformation. Then, a system of two-dimensional shear deformation equations is derived in differential and fully variational forms for the high frequency waves and vibrations of a functionally graded shell. Also, a theorem is given, which states the conditions sufficient for the uniqueness in solutions of the shell equations. On the basis of a discrete layer modeling, the governing equations are obtained for the motions of a curved laminae made of any numbers of functionally graded distinct layers, whenever the displacements and the electric and magnetic potentials of a layer are taken to vary linearly across its thickness. The resulting equations in differential and fully variational, invariant forms account for various types of waves and coupled vibrations of one and two dimensional structural elements as well. The invariant form makes it possible to express the equations in a particular coordinate system most suitable to the geometry of shell (plate) or laminae. The results are shown to be compatible with and to recover some of earlier equations of plane and curved elements for special material, geometry and/or effects.     

Two-dimensional transient thermo-hydro-mechanical fundamental solutions of multiphase porous media in frequency and time domains

In this paper, the closed form two-dimensional fundamental solutions for a non-isothermal unsaturated deformable porous medium have been derived for a symmetric polar domain in both Laplace transform and time domains. The governing differential equations of the non-isothermal unsaturated soil consist of equilibrium, moisture, air and heat transfer equations including the suction effect, temperature effect and dissolved air in water. The derived fundamental solution has been verified mathematically by comparison with the previously presented corresponding fundamental solutions in three limiting cases including the steady-state thermo-hydro-mechanical, steady-state hydro-mechanical and elastostatic fundamental solutions. Also these 2D kernel functions are tested in comparison with a finite element method (FEM).     

轴对称环形片状界面裂纹问题分析

讨论受拉伸载荷作用的轴对称环形片状界而裂纹问题．该问题归结为求解一组超奇异积分-微分方程．方程中的未知位移间断近似表示为基本密度函数与多项式之积,其中基本密度函数考虑到问题的对称性用二维界面裂纹精确解表示．在圆形片状裂纹的情况下,数值结果与现有理论解作比较的结果表明,数值结果与相应界面圆形片状裂纹和均质体圆形片状裂纹的精确解均吻合得很好．文中以图表形式给出应力强度因子与材料组合和几何条件之间的关系．     

Characteristic equation solution strategy for deriving fundamental analytical solutions of 3D isotropic elasticity

A simple characteristic equation solution strategy for deriving the fundamental analytical solutions of 3D isotropic elasticity is proposed. By calculating the determinant of the differential operator matrix obtained from the governing equations of 3D elasticity, the characteristic equation which the characteristic general solution vectors must satisfy is established. Then, by substitution of the characteristic general solution vectors, which satisfy various reduced characteristic equations, into various reduced adjoint matrices of the differential operator matrix, the corresponding fundamental analytical solutions for isotropic 3D elasticity, including Boussinesq-Galerkin (B-G) solutions, modified Papkovich-Neuber solutions proposed by Min-zhong WANG (P-N-W), and quasi HU Hai-chang solutions, can be obtained. Furthermore, the independence characters of various fundamental solutions in polynomial form are also discussed in detail. These works provide a basis for constructing complete and independent analytical trial functions used in numerical methods.     

Singular boundary method for 2D dynamic poroelastic problems

This paper extends a strong-form meshless boundary collocation method, named the singular boundary method (SBM), for the solution of dynamic poroelastic problems in the frequency domain, which is governed by Biot equations in the form of mixed displacement–pressure formulation. The solutions to problems are represented by using the fundamental solutions of the governing equations in the SBM formulations. To isolate the singularities of the fundamental solutions, the SBM uses the concept of the origin intensity factors to allow the source points to be placed on the physical boundary coinciding with collocation points, which avoids the auxiliary boundary issue of the method of fundamental solutions (MFS). Combining with the origin intensity factors of Laplace and plane strain elastostatic problems, this study derives the SBM formulations for poroelastic problems. Five examples for 2D poroelastic problems are examined to demonstrate the efficiency and accuracy of the present method. In particular, we test the SBM to the multiply connected domain problem, the multilayer problem and the poroelastic problem with corner stress singularities, which are all under varied ranges of frequencies.     

Boundary element method for orthotropic thick plates

The fundamental solutions of the orthotropic thick plates taking into account the transverse shear deformation are derived by means of Hörmander's operator method and a plane-wave decomposition of the Dirac δ-function in this papey. The boundary integral equations of the thick plates have been formulated which are adapted to arbitrary boundary conditions and plane forms. The numerical calculation of the fundamental solutions is discussed in detail. Some numerical examples are analyzed with BEM.     

A non-linear rod theory for high-frequency vibrations of thermopiezoelectric materials

In relation to electroelastic media with thermopiezoelectric coupling, the system of one-dimensional equations is consistently derived so as to accommodate the high-frequency vibrations of a rod with temperature-dependent material. In the first part of the paper, a unified variational principle of differential type is presented which describes the fundamental equations of thermopiezoelectricity with second sound, including the physical and geometrical non-linearities. In the second part, the hierarchic system of rod equations is systematically deduced from the three-dimensional fundamental equations by use of Mindlin's method of reduction. The hierarchic system of equations which is derived in both differential and variational forms is capable of predicting the extensional, thickness-shear, flexural and torsional as well as coupled vibrations of the rod of uniform cross-section. All the higher-order effects are taken into account as deemed pertinent in any particular case. In the third part, attention is confined to certain cases involving special motions, materials and geometry. Besides, the uniqueness is investigated in solutions of the linearized system of rod equations and the sufficient conditions are enumerated for the uniqueness of solutions.     

New analytic solutions for static problems of rectangular thin plates point-supported at three corners

This paper deals with the bending of rectangular thin plates point-supported at three corners using an analytic symplectic superposition method. The problems are of fundamental importance in both civil and mechanical engineering, but there were no accurate analytic solutions reported in the literature. This is attributed to the difficulty in seeking the solutions that satisfy the governing fourth-order partial differential equation with the free boundary conditions at all the edges as well as the support conditions at the corners. In the following, the Hamiltonian system-based equation for plate bending is formulated, and two types of fundamental problems are analytically solved by the symplectic method. The analytic solutions of the plates point-supported at three corners are then obtained by superposition, where the constants are obtained by a set of linear equations. The solution procedure presented in this paper offers a rigorous way to yield analytic solutions of similar problems. Some numerical results, validated by the finite element method, are shown to provide useful benchmarks for comparison and validation of other solution methods.     

On the fundamental equations of piezoelasticity of quasicrystal media

The three-dimensional fundamental equations of elasticity of quasicrystals with extension to quasi-static electric effect are expresses in both differential and variational invariant forms for a regular region of quasicrystal material. The principle of conservation of energy is stated for the regular region and the constitutive relations are obtained for the piezoelasticity of material. A theorem is proved for the uniqueness in solutions of the fundamental equations by means of the energy argument. The sufficient boundary and initial conditions are enumerated for the uniqueness. Hamilton’s principle is stated for the regular region and a three-field variational principle is obtained under some constraint conditions. The constraint conditions, which are generally undesirable in computation, are removed by applying an involutory transformation. Then, a unified variational principle is obtained for the regular region, with one or more fixed internal surface of discontinuity. The variational principle operating on all the field variables generates all the fundamental equations of piezoelasticity of quasicrystals under the symmetry conditions of the phonon stress tensor and the initial conditions. The resulting equations, which are expressible in any system of coordinates and may be used through simultaneous approximation upon all the field variables in a direct method of solutions, pave the way to the study of important dislocation, fracture and interface problems of both elasticity and piezoelasticity of quasicrystal materials.