On solutions of the problem in stresses with the use of Maxwell stress functions

The spatial problems of elasticity are mainly solved in displacements [1, 2], i.e., the Lamé equations are taken as the initial equations. This is related to the lack of general solutions for the system of basic equations of elasticity expressed in stresses. In this connection, a new variational statement of the problem in stresses was developed in [3, 4]; this statement consists in solving six generalized equations of compatibility for six independent components of the stress tensor, while the three equilibrium equations are transferred to the set of boundary conditions. This method is more convenient for the numerical solution of problems in stresses and has been tested when solving various boundary value problems. In the present paper, analyzing the completeness of the Saint-Venant identities and using the Maxwell stress functions, we obtain a new resolving system of three differential equations of strain compatibility for the three desired stress functions φ, ξ, and ψ. This system is an alternative to the three Lamé equilibrium equations for three desired displacement components u, v, w and is simpler in structure. Moreover, both of these systems of resolving equations can be solved by the new recursive-operator method [5, 6]. In contrast to well-known methods for constructing general solutions of linear differential equations and their systems, the solutions obtained by the recursive-operator method are constructed as operator-power series acting on arbitrary analytic functions of real variables (not necessarily harmonic), and the series coefficients are determined from recursive relations (matrix in the case of systems of equations). The arbitrary functions contained in the general solution can be determined directly either from the boundary conditions (the obtained system of inhomogeneous equations with a right-hand side can also be solved by the recursive-operator method [6]) or by choosing them from various classes of analytic functions (elementary, special); a complete set of particular solutions can be obtained in the same function classes, and the coefficients of linear combinations of particular solutions can be determined by the Trefftz method, the least-squares method, and the collocation method.     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

无网格方法在平面粘弹性力学问题中的应用

本文综合应用无网格方法(EFGM)、线性粘弹性与弹性力学之间的对应原理,Laplace变换和逆变换等方法求解了拟静态平面弹性和粘弹性力学问题。首先,利用Laplace变换和逆变换推导了平面问题的粘弹性本构关系,建立了拟静态粘弹性平面问题的边值问题;其次,利用粘弹性与弹性力学之间的对应原理得到了Laplace变换域中平面问题的基本方程,在Laplace变换域中建立了相应的泛函,并得到了用无网格方法离散的控制方程;同时,求解了几个拟静态弹性和粘弹性平面问题,给出了它们的表达式和数值结果;最后,采用Laplace逆变换和数值逆变换,得到了粘弹性力学平面问题在物理空间中的解,并比较了由解析解和无网格数值方法所得到的数值结果,可以看到它们是非常吻合的。说明本文方法的正确性和有效性。     

Asymptotic solutions of coupled equations of supercritically axially moving beam

In supercritical regime, the coupled model equations for the axially moving beam with simple support boundary conditions are considered. The critical speed is determined by linear bifurcation analysis, which is in agreement with the results in the literature. For the corresponding static equilibrium state, the second-order asymptotic nontrivial solutions are obtained through the multiple scales method. Meantime, the numerical solutions are also obtained based on the finite difference method. Comparisons among the analytical solutions, numerical solutions and solutions of integro-partial-differential equation of transverse which is deduced from coupled model equations are made. We find that the second-order asymptotic analytical solutions can well capture the nontrivial equilibrium state regardless of the amplitude of transverse displacement. However, the integro-partial-differential equation is only valid for the weak small-amplitude vibration axially moving slender beams.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

Dynamic instability of truncated conical shells under periodic axial load

Based on the Dynamic version of Donnell type basic equations neglecting bending deformations before instability, the parametric instability of truncated conical shells subjected to periodic axial load is studied under four different boundary conditions. Applying Galerkin's method, the basic equations are reduced to a system of coupled Mathieu equations, from which the instability regions are determined by using Hsu's results. As a numerical example, instability regions for a completely clamped shell are determined for relatively wide range of frequencies. The effects of static axial load as well as damping force on the instability regions are also examined.     

功能梯度材料结构分析的半解析数值方法研究

一种典型的半解析数值方法——线法被引入功能梯度材料的结构分析。首先推导了功能梯度材料位移形式的平衡方程和边界条件,然后阐述了线法功能梯度材料结构分析的基本步骤和数值原理。该方法的基本思想是通过有限差分将问题的控制方程半离散为定义在沿梯度方向离散节线上的常微分方程组,然后应用B样条函数Gauss配点法求解该常微分方程组得到问题的解答。为演示线法在功能梯度材料结构分析中的应用,给出了线性梯度和指数梯度功能梯度材料板分别受恒定位移、均匀拉伸载荷和弯曲载荷作用的数值算例。与相应问题解析解和其他数值方法的比较表明,线法的计算结果具有很高的精度,而且不需要任何特殊的考虑就能够有效模拟材料内部物性参数的连续变化,也无需事先选取满足特定条件的待定场函数,是一种非常适合功能梯度材料结构形式和材料特点的半解析数值方法。     

复合材料层合板的三维研究

本文提出了固支复合材料各向异性层合圆板受均布横向载荷作用下的满足三维弹性力学基本微分方程和边界条件的解析解答。文中采用一种发展的摄动方法进行求解，板中的每个应力和位移都展开为无量纲厚度参数ε的摄动级数，并采用二维板理论解答作为其相应三维摄动解答的一个基本解的形式，通过摄动方法逐级求解而获得完整的三维解答。文中以解析形式和数值形式给出了高精确度的三维应力和位移结果，结果表明，本文求解三维问题的解析方法是合理有效的。     

多体系统接触碰撞问题的牛顿-欧拉线性互补方法

基于非光滑动力学方法的多体系统接触碰撞分析是目前多体系统动力学的研究热点.本文采用牛顿-欧拉方法建立多体系统接触、碰撞问题的动力学模型,给出一种牛顿-欧拉型线性互补公式.该建模方法与目前一般采用的拉格朗日建模方法的不同之处是约束条件中除了库仑摩擦、单边约束之外还含有光滑等式约束.在建立系统动力学模型时,首先解除摩擦约束和单边约束得到原系统对应的基本系统.牛顿-欧拉方法采用最大数目坐标建立基本系统的动力学方程,由于坐标不相互独立,因此基本系统中带有等式约束,其数学模型为一组微分代数方程.借助约束雅可比矩阵,在基本系统微分代数方程中添加摩擦接触和单边约束对应的拉氏乘子,就可以得到系统全局运动的具有变拓扑结构特征的动力学方程,再结合非光滑约束互补条件便可构成完备的系统动力学模型.完备的动力学模型由动力学微分方程以及等式约束和不等式约束组成.线性互补公式采用分块矩阵形式进行推导,简化了推导过程.数值计算采用基于线性互补的时间步进算法.时间步进算法是目前流行的非光滑数值算法,其突出特点是可以免去数值积分中繁琐的事件检测过程,而数值积分过程中通过对线性互补问题的求解可以确定系统的触-离状态.通过对典型的曲柄滑块间隙机构进行数值分析,验证本文方法的有效性.     

Implementation of high‐order compact finite‐difference method to parabolized Navier–Stokes schemes

The numerical solution to the parabolized Navier–Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth‐order compact finite‐difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite‐difference algorithm of Beam and Warming type with a high‐order compact accuracy. A shock‐fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously computed with the fourth‐order accuracy. The computations are carried out for a benchmark case: hypersonic axisymmetric flow over a blunt cone at Mach 8. A sensitivity study is performed for the basic flowfield, including profiles and their derivatives obtained from the fourth‐order compact PNS and IPNS solutions, and the effects of grid size and numerical dissipation term used are discussed. The present results for the flowfield variables and also their derivatives are compared with those of other basic flow models to demonstrate the accuracy and efficiency of the proposed method. The present work represents the first known application of a high‐order compact finite‐difference method to the PNS schemes, which are computationally more efficient than Navier–Stokes solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

Nonlinear three-dimension analysis for axially symmetrical circular plates and multilayered plates

Analytic nonlinear three-dimension solutions are presented for axially symmetrical homogeneous isotropic circular plates and multilayered plates with rigidly clamped boundary conditions and under transverse load.The geometric nonlinearily from a moderately large deflection is considered.A developmental perturbation method is used to solve the complicated nonlinear three-dimension differential equations of equilibrium.The basic idea of this perturbation method is using the two-dimension solutions as a basic form of the corresponding three-dimension solutions,and then processing the perturbation procedure to obtain the three-dimension perturbation solutions.The nonlinear three-dimension results in analytic expressions and in numerical forms for ordinary plates and multilayered plates are presented.All of the plate stresses are shown in figures.The results show that this perturbation method used to analyse nonlinear three-dimension problems of plates is effective.     

层体边界元法在二维层体断裂分析中的应用

在刚度矩阵法的基础上建立了用于进行二维多层体结构断裂分析的边界单元法（ＢＥＭＬＭ）由于ＢＥＭＬＭ的基本方程中已经包含了层体表面和裂纹缝面的边界条件，因而不需要对这些边界进行单元离散，从而其断裂分析可望有较好的精度通过与柯西积分方程法进行结合，算例表明ＢＥ ＭＬＭ是可靠并有效的     

求解不连续中厚板自由振动的微分容积单元法

基于区域叠加原理和微分容积法，发展了一种新型的数值方法——微分容积单元法，用以分析具有不连续几何特征的中厚板的自由振动。根据板的不连续情况将其划分为若干单元，在每个单元内用微分容积法将控制微分方程离散成为一组线性代数方程．在相邻的单元连接处应用位移连续条件和平衡条件，引入边界约束条件后得到一套关于各配点位移的齐次线性代数方程，由此可导出求解系统固有频率的特征方程。本文用子空间迭代法求解特征方程，并以开孔板、混合边界条件板和突变厚度板为例研究了方法的收敛性和计算精度。     

A mixture differential quadrature method for solving two-dimensional incompressible Navier-Stokes equations

Ｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ（ＤＱＭ）ｐｒｏｐｏｓｅｄｂｙＲ．Ｂｅｌｌｍａｎ［１,２］ｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｅｍｐｌｏｙｅｄｉｎｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｓｏｆｐｒｏｂｌｅｍｓｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｈｙｓｉｃａｌｓｃｉｅｎｃｅ．ＢｅｃａｕｓｅｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｎａｌｌｇｒｉｄｐｏｉｎｔｓｉｓｕｓｅｄｔｏｆｉｔｔｈｅｄｅｒｉｖａｔｉｖｅｓａｔｇｒｉｄｐｏｉｎｔｓｉｎｔｈｅＤＱＭ,ｉｔｉｓｅｎｏｕｇｈｔｏｏｂｔａ…     

The study of quasi wavelets based numerical method applied to Burgers' equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｏｎｃｅｐｔｏｆｗａｖｅｌｅｔｗａｓｆｉｒｓｔｉｎｔｒｏｄｕｃｅｄｂｙｇｅｏｐｈｙｓｉｃｉｓｔＭｏｒｌｅｔ[1].Ａｗａｖｅｌｅｔｆｕｎｃｔｉｏｎｔｈａｔｉｓｆｉｎｉｔｅ＿ｅｎｅｒｇｙｆｕｎｃｔｉｏｎｓｗｉｔｈｗｅｌｌｌｏｃａｌｉｚａｔｉｏｎｐｒｏｐｅｒｔｉｅｓｃａｎｂｅｄｅｓｉｇｎａｔｅｄｂｙｄｉｌａｔｉｏｎａｎｄｔｒａｎｓｌａｔｉｏｎ .Ｔｈｅａｂｉｌｉｔｙｏｆｗａｖｅｌｅｔｂａｓｅｄｍｅｔｈｏｄｔｏａｎａｌｙｓｅｆｕｎｃｔｉｏｎｓｉｎｔｅｒｍｏｆｔｈｅｉｒｌｏｃａ…     

Non-linear and linear evolution of perturbation in stochastic basic flows

In this paper, we examine the non-linear and linear evolutions of perturbation in stochastic basic flows with two-dimensional quasi-geostrophic equations on a sphere. As the analytic solutions for the considered quasi-geostrophic equations are not available, the Fourier finite volume element method is used to perform numerical simulation. It is found that, the non-linear and linear evolutions of perturbation in stochastic basic flow will be consistent for a short period of time and small stochastic fluctuations when they are consistent in the deterministic basic flow. However, the tangent linear model will fail to approximate the original non-linear model when the time period is considerably long and stochastic fluctuation becomes large. Moreover, the global energy decays faster for stochastic basic flow with stronger fluctuations.     

矩形悬臂薄板的解析解

首先把弹性薄板弯曲问题的控制方程表示成为Hamilton正则方程,然后利用辛几何方法对全状态相变量进行分离变量,求出其本征值后,再按本征函数展开的方法求出矩形悬臂薄板的解析解。由于在求解过程中不需要事先人为地选取挠度函数,而是从薄板弯曲的基本方程出发,直接利用数学的方法求出可以满足其边界条件的这类问题的解析解,使得问题的求解更加理论化和合理化。文中的最后还给出了计算实例来验证本文所采用的方法以及所推导出的公式的正确性。     

计入横向剪切变形时层合板位移差分解法

在考虑横向剪切变形对层合板弹性解的影响时，本文提出一种数值计算方法。由边界条件给出边界结点位移的表达式，根据薄板的经典理论和一阶横向剪切变形理论导出位移增量所满足的平衡微分方程，引用经典理论计算的横向剪力修正了荷载列阵。致使在较粗的网格划分时、宽广的层合板长厚比范围内，仍能得到与解析解颇为一致的数值解。     

Analysis of a laminar boundary layer flow over a flat plate with injection or suction

An analysis is performed to study a laminar boundary layer flow over a porous flat plate with injection or suction imposed at the wall. The basic equations of this problem are reduced to a system of nonlinear ordinary differential equations by means of appropriate transformations. These equations are solved analytically by the optimal homotopy asymptotic method (OHAM), and the solutions are compared with the numerical solution (NS). The effect of uniform suction/injection on the heat transfer and velocity profile is discussed. A constant surface temperature in thermal boundary conditions is used for the horizontal flat plate.     

任意四边形薄板大挠度弯曲问题研究

本文用pb-2 Ritz能量法求解任意四边形薄板的几何非线性弯曲问题.首先通过坐标变换将任意四边形区域转换到一个2×2单位正方形求解区域,并建立求解域内的能量泛函,然后取一个完备的二元多项式级数(p-2)与描述边界形状的一个基本函数(b)的乘积作为Ritz函数,由能量最小原理建立结构的刚度方程,非线性代数方程采用拟牛顿法求解.文中给出了详细的数学公式,并对四边简支方板和四边固定菱形板进行了数值分析,算例表明,本方法具有公式简单、精度较高的优点,用以分析大挠度问题是非常有效的.