集中载荷作用下不同边界条件的弯曲厚矩形板的受迫振动

运用边界积分法研究了四边简支、两对边固定另两对边简支、四边固定三种复杂边界条件下厚矩形板的受迫振动问题,求解过程清晰,从而给出了受迫振动控制方程和挠曲面方程。通过在Matlab平台上进行数值计算,得出了图表形式的计算结果,并与有限元模拟值进行对照。研究表明,边界积分法用于求解厚矩形板的受迫振动问题的准确性,本文推导的控制方程和挠曲面方程的正确性,进而对工程实际中的各种相关问题具有一定的现实意义,也为求解此类问题提供了一种新途径,可以直接运用到工程实际中。     

功能梯度材料动态断裂力学的径向积分边界元法

采用径向积分边界元法分析功能梯度材料动态断裂力学问题. 该方法使用与弹性模量无关的弹性静力学开尔文基本解作为问题的基本解,在导出的边界-域积分方程中含有由材料的非均质性和惯性项引起的域积分,通过径向积分法将域积分转化为等效的边界积分,得到只含边界积分的纯边界积分方程;从而建立只需边界离散的无内部网格边界元算法. 采用候博特方法求解关于时间二阶导数的系统离散的常微分方程组. 最后通过数值算例验证本文方法的精度和有效性.      

获得热传导问题“拟解析解”的精细积分算法

研究并讨论了用于热传导分析有限元解的精细积分算法，算法很好地克服了传统方法求解时的单调性问题，且对空间离散后所获得方程的解是解析的，因而算法的解将具有“拟解析解”的意义，论文证明了算法单调性．     

一期水管冷却效应的数值模拟新方法

采用有限单元离散空间域，视冷却水管所在单元为子结构，其它单元为常规单元，然后，把混凝土结构瞬态温度场的有限元控制方程表示成状态方程，最后，利用精细积分法求解．结果表明，子结构、状态方程和精细积分技术的引用降低了对计算机贮存量的要求和提高了冷却水管附近温度梯度的精度，显示出极大的优越性．     

Solution of the Incompressible Unsteady Navier–Stokes Equations Only in Terms of the Velocity Components

A method which uses only the velocity components as primitive variables is described for solution of the incompressible unsteady Navier–Stokes equations. The method involves the multiplication of the primitive variable-based Navier–Stokes equations with the unit normal vector of finite volume elements and the integration of the resulting equations along the boundaries of four-node quadrilateral finite volume elements. Therefore, the pressure term is eliminated from the governing equations and any difficulty associated with pressure or vorticity boundary conditions is avoided. The equations are discretized on four-node quadrilateral finite volume elements by using the second-order-accurate central finite differences with the mid-point integral rule in space and the first-order-accurate backward finite differences in time. The resulting system of algebraic equations is solved in coupled form using a direct solver. As a test case, an impulsively accelerated lid-driven cavity flow in a square enclosure is solved in order to verify the accuracy of the present method.     

THE MESHLESS VIRTUAL BOUNDARY METHOD AND ITS APPLICATIONS TO 2D ELASTICITY PROBLEMS

A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network (IRBFN) constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations (BIE) can be established while the fundamental solutions to the problems are introduced. The singular value decomposition (SVD) method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.     

Tractive rolling contact mechanics of graded coatings

In this study, the tractive rolling contact problem between a rigid cylinder and a graded coating is investigated. The main objective of this study is to investigate the effect of the stiffness ratio, the coefficient of friction and the coating thickness on the surface contact tractions, the surface in-plane stress, the stick zone length and the creep ratio parameter that may have a bearing on the fatigue life of the component. Assuming that the shear modulus varies exponentially through the thickness of the coating, the governing integral equations associated with the rolling contact problem are constructed. Furthermore, it is supposed that the contact patch is controlled by a central stick zone accompanied by two slip zones. The conventional Goodman approximation is employed in order to decouple the governing singular integral equations. Finally, the numerical solution of the integral equations is obtained by applying the Gauss–Chebyshev integration method.     

下卧刚性基岩的饱和地基上基础的动力分析

采用积分变换与积分方程研究了单层饱和地基上刚性基础的竖向振动问题，首先从饱和土的动力控制方程出发，考虑混合边值条件，获得了一组描述基础振动的对偶积分方程，进而给出了单层饱和地基表面的动力柔度系数。该文的工作可视为对饱和半空间地基上基础振动解答的推广。     

具有混合边界透水条件的多孔饱和半空间上刚性圆板的垂直振动

用积分变换和积分方程研究多孔饱和半空间上刚性圆板的垂直振动问题。首先应用逐次解耦方法求解多孔饱和固体的动力基本方程－Ｂｉｏｔ波动方程。然后考虑混合边界透水条件（半空间表面与圆板的接触面是不透水的,而其余表面是透水的）,建立子多孔饱和半空间上刚性圆板垂直振动的对偶积分方程,并化对偶积分方程为第二类Ｆｒｄｄｈｏｌｍ积分方程。     

SCATTERING OF ANTI-PLANE SHEAR WAVES BY A SINGLE CRACK IN AN UNBOUNDED TRANSVERSELY ISOTROPIC ELECTRO-MAGNETO-ELASTIC MEDIUM

A theoretical treatment of the scattering of anti-plane shear (SH) waves is provided by a single crack in an unbounded transversely isotropic electro-magneto-elastic medium. Based on the differential equations of equilibrium, electric displacement and magnetic induction intensity differential equations, the governing equations for SH waves were obtained. By means of a linear transform, the governing equations were reduced to one Helmholtz and two Laplace equations. The Cauchy singular integral equations were gained by making use of Fourier transform and adopting electro-magneto imperme ableboundary conditions. The closed form expression for the resulting stress intensity factor at the crack was achieved by solving the appropriate singular integral equations using Chebyshev polynomial. Typical examples are provided to show the loading frequency upon the local stress fields around the crack tips. The study reveals the importance of the electro-magneto-mechanical coupling terms upon the resulting dynamic stress intensity factor.     

多孔饱和半空间上刚体垂直振动的轴对称混合边值问题

研究圆柱形刚体在多孔饱和半空间上的垂直振动．首先应用Ｈａｎｋｅｌ变换求解多孔饱和固体的动力基本方程———Ｂｉｏｔ波动方程．然后按混合边值条件建立多孔饱和半空间上刚体垂直振动的对偶积分方程，用Ａｂｅｌ变换化对偶积分方程为第二类Ｆｒｅｄｈｏｌｍ积分方程．文末给出了多孔饱和半空间表面动力柔度系数的计算曲线．     

Dynamic behavior of piezoelectric/magnetostrictive composite hollow cylinder

The dynamic behavior of a multilayered, perfectly bonded piezoelectric/magnetostrictive composite hollow cylinder under radial deformation is investigated. The superposition method, the state space method as well as the separation of variables method are elegantly integrated in the solution approach. The governing equations are finally transformed into two Volterra integral equations of the second kind with respect to two functions of time. The elastic, electric and magnetic fields are finally obtained according to solving the integral equations. Free vibrations and transient responses are demonstrated by numerical experiments.     

A novel statistical linearization solution for randomly excited coupled bending-torsional beams resting on non-linear supports

A novel statistical linearization technique is developed for computing stationary response statistics of randomly excited coupled bending-torsional beams resting on non-linear elastic supports. The key point of the proposed technique consists in representing the non-linear coupled response in terms of constrained linear modes. The resulting set of non-linear equations governing the modal amplitudes is then replaced by an equivalent linear one via a classical statistical error minimization procedure, which provides algebraic non-linear equations for the second-order statistics of the beam response, readily solved by a simple iterative scheme. Data from Monte Carlo simulations, generated by a pertinent boundary integral method in conjunction with a Newmark numerical integration scheme, are used as benchmark solutions to check accuracy and reliability of the proposed statistical linearization technique.

     

上覆单相弹性层的饱和地基上刚性圆板的扭转振动分析

采用解析的方法研究了上覆单相弹层的饱和地基上刚性圆板的扭转振动。首先运用积分变换技术分别求解了单相弹性介质和饱和介质情况时的控制方程，然后按混合边值条件建立了上覆单相弹性层的饱和地基上刚性圆板扭转振动的对偶积分方程，并把对偶积分方程化为易于数值求解的第二类Fredholm积分方程，并给出了数值算例。     

三维N-S方程计算隐式有限体积法

利用NND有限差分格式，发展了一种新的完全隐式的有限体积数值方法，以求解与时间相关的N－S方程．对通过单元体界面的无粘流和粘性流通量均作隐式处理．对绕流钝锥体和不同攻角的气动辅助实验飞行器的高超声速粘性流和化学反应流获得了定常数值解．对流加热率和流场电子密度的计算值与实验数据进行了比较，符合较好，证实了本方法的精确性．     

The torsion of a non-homogeneous stratum with a hyperbolic variation of shear modulus

An exact formulation of the governing dual integral equations for the torsion of a non-homogeneous stratum due to a rigid circular body at its free surface is presented. The stratum varies in shear modulus according to the hyperbolic variation in a contemporary work [1]. It is shown that the unknown static stress distribution under the rigid body is governed by modified Bessel function of the first kind. By comparing the governing functions in the dual integral equations for five cases of elastic media: homogeneous half-space, and stratum, linearly non-homogeneous half-space and stratum and, finally, the present non-homogeneous stratum with hyperbolic variation, it is established that the surface shear modulus is the dominant parameter in the assessment of the stress and displacement fields in a non-homogeneous stratum where lateral variation of elastic properties is negligible.     

瞬态热传导问题的精细积分数值流形方法研究

数值流形方法（NMM）因其特有的双覆盖系统（数学覆盖和物理覆盖）在域离散方面具有独特的优势,而精细时间积分法则具有精度高、无条件稳定、无振荡以及计算结果不依赖于时间步长等特点。发展了用于研究二维瞬态热传导问题的精细积分NMM。结合待求问题的控制方程和边界条件,并基于修正变分原理导出了NMM的总体方程,给出了求解此类时间相依方程的精细时间积分及空间积分策略,选取了两个典型算例对方法的有效性进行了验证,结果表明本文方法可以高效高精度地求解瞬态热传导问题。     

Local boundary integral equations for orthotropic shallow shells

A meshless local Petrov–Galerkin (MLPG) formulation is presented for bending problems of shear deformable shallow shells with orthotropic material properties. Shear deformation of shells described by the Reissner theory is considered. Analyses of shells under static and dynamic loads are given here. For transient elastodynamic case the Laplace-transform is used to eliminate the time dependence of the field variables. A weak formulation with a unit test function transforms the set of governing equations into local integral equations on local subdomains in the plane domain of the shell. Nodal points are randomly spread in that domain and each node is surrounded by a circular subdomain to which local integral equations are applied. The meshless approximation based on the moving least-squares (MLS) method is employed for the implementation. Unknown Laplace-transformed quantities are computed from the local boundary integral equations. The time-dependent values are obtained by the Stehfest’s inversion technique.     

四边简支矩形中厚板的弯曲

本文采用Reissner中厚板理论求解了四边简支矩形中厚板的弯曲问题。文中首先对Reissner中厚板理论的控制方程进行了适当的变更,使之成为非耦联的二阶偏微分方程组,然后利用有限积分变换法求解所得新的控制方程,得到了四边简支矩形中厚板受均布载荷作用下的解析解。文中所述方法可用以求解具有其它边界条件和载荷的矩形中厚板的弯曲问题,同时还可移植应用于其它中厚板理论。     

Well-posedness of two-phase local/nonlocal integral polar models for consistent axisymmetric bending of circular microplates

Previous studies have shown that Eringen's differential nonlocal model would lead to the ill-posed mathematical formulation for axisymmetric bending of circular microplates. Based on the nonlocal integral models along the radial and circumferential directions, we propose nonlocal integral polar models in this work. The proposed strainand stress-driven two-phase nonlocal integral polar models are applied to model the axisymmetric bending of circular microplates. The governing differential equations and boundary conditions (BCs) as well as constitutive constraints are deduced. It is found that the purely strain-driven nonlocal integral polar model turns to a traditional nonlocal differential polar model if the constitutive constraints are neglected. Meanwhile, the purely strain-and stress-driven nonlocal integral polar models are ill-posed, because the total number of the differential orders of the governing equations is less than that of the BCs plus constitutive constraints. Several nominal variables are introduced to simplify the mathematical expression, and the general differential quadrature method (GDQM) is applied to obtain the numerical solutions. The results from the current models (CMs) are compared with the data in the literature. It is clearly established that the consistent softening and toughening effects can be obtained for the strain-and stress-driven local/nonlocal integral polar models, respectively. The proposed two-phase local/nonlocal integral polar models (TPNIPMs) may provide an e-cient method to design and optimize the plate-like structures for microelectro-mechanical systems.