具有压电材料弹性环形板的屈曲

对在不同位置粘有任意多对压电片的弹性环形板在电压和径向压力作用下的轴对称屈曲进行了研究．根据压电效应的等效作用量，初参数法以及传递矩阵方法得到了环板屈曲的特征方程．对不同位置及不同宽度的压电片，计算了结构的屈曲载荷．并得到了稳定性边界曲线．     

分区样条等参元方法分析加肋轴对称组合壳

本文提出分区样条等参元方法分析加肋轴对称组合壳在静水外压作用下的应力及稳定性。该方法根据母线几何形状、壳板厚度变化及边界约束状况将加肋轴对称组合壳划分为若干个区域，在不同区域采用不同步长的样条等参元进行分析，因而便于对复杂的组合壳进行分析，并可大大减少计算所需自由度数，减少计算工作量，同时计算精度很高，并且易于进行边界约束处理。应用该方法对加肋锥－柱、加肋锥－环－柱组合壳进行了计算，其结果与其它方法计算值及实验值吻合良好。     

分区样条等参元方法分析加助轴对称组合壳

本文提出分区样条行动等参元方法分析加肋轴对称组合壳在静水外压作用下的应力及稳定性，该方法根据母线几何形状，壳板厚度变化及边界约束状况将加肋轴对称组合壳划分为若干个区域，在不同区域采用不同步长的样条等参元进行分析，因而便于对复杂的组合壳进行分析，并可大大减少计算所需自由度数，减少计算工作量，同时计算精度很高，并且易于进行边界约束处理，应用该方法对加肋锥－柱，加肋锥－环－柱组合壳进行了计算，其结果与其     

层状饱和土Biot固结问题状态空间法

针对饱和多孔介质空间非轴对Biot固结问题，引入状态变量，构造了两组相比独立的状态变量方程，利用Fourier级数和Laplace-Hankel变换，将状态变量方程转换为两组一阶常微分方程组，提出了均质饱和多孔介质空间非轴对称Biot固结问题的传递矩阵，得到以状态变量和传递矩阵乘积的形式表示的均质饱和多孔介质空间非轴对称Biot固结问题的解，利用层间完全接触的条件，可得到N层饱和多孔介质空间非轴对称Biot固结问题的一般解析表达式，文中考虑几种不同的边界条件，分析了两个算例，数值结果表明该方法具有较高的计算精度和良好的计算稳定性。     

葵花型索穹顶结构等效平面桁架静力计算方法

选取n次轴对称葵花型索穹顶结构的1/n,将空间轴对称结构等效为一次超静定平面桁架结构.给出了葵花型索穹顶结构撑杆不在同一径向上的内力和位移处理方法,推导了力法求解约束力的计算公式,提出了小变形时各点位移计算的位移传递矩阵和计算公式.最后,通过一个算例介绍了该方法的计算过程并与有限元计算结果进行了比较,内力和位移误差均在工程设计要求内,验证了方法的正确性和准确性.     

用积分变换及边界积方法求解多层地基的静力问题

本文利用积分变换及矩阵递推方法得到了任意ｎ层弹性体平面应变及轴对称问题的ｍｉｎｄｌｉｎ解。再把此解作为基本解，利用Ｓｏｍｉｇｌｉａｎａ关系式，得到计算多层弹性体内部任意点位移的简便方法，利用此法很容易编制程序，且具有较高的计算精度与速度。     

各向异性展状介质中弹性波的传播矩阵解法

基于广义胡克定律及混和变量弹性波方程，解析求得各层介质内位移、应力传递矩阵，给出了直角坐标系下各向异性层状介质中弹性波的传播矩阵解法．该方法适用于非轴对称各向异性和点源作用，较好地解决了数值计算中有效数字精度损失问题．数值结果表明，计算效率、准确性及稳定性均较好．     

层状地基轴对称问题的Mindlin解

本文在文献〔１〕的基础上，对层状地基应力和位移的计算方法做了一些创新的工作，首先作者提出求解传矩阵的逆矩阵的简便公式。并在应力、位移的计算上，对荷载作用面上下分别采用二种不同的传递矩阵。加快了计算速度，避免了层状体系计算中常见的溢出问题及大数减大数造成的计算困难。     

求解多层弹性半空间轴对称问题的精确刚度矩阵法

本文首先从弹性力学的基本方程出发，利用Hankel积分变换等数学手段，推导出了单层弹性半空问轴对称问题的刚度矩阵，然后按传统的有限元方法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程和Hankel积分逆变换就可解出静荷载作用下多层弹性半空间轴对称问题的精确解。由于刚度矩阵的元素中不含有正指数项，计算时不会出现溢出的现象，从而克服了传递矩阵法的缺点。由于在推导过程中摒弃了应力函数的选择，使得问题的求解更加理论化和合理化，同时也为进一步研究这类问题如温度场，动力学等方向奠定了理论基础。最后，文中还给出了计算实例来证明推导结果的准确性。     

求解动荷载作用下多层粘弹性半空间轴对称问题的精确刚度矩阵法

利用Laplace—Hankel联合积分变换，推导出了单层粘弹性半空间轴对称问题在动荷载作用下，层间完全接触情况的刚度矩阵，然后按传统的有限元方法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程就可解出动荷载作用下多层粘弹性半空间轴对称问题的矩阵。由于刚度矩阵的元素中只含有负指数项，计算时不会出现溢出的现象。本文还成功地应用了Durbin的Laplace逆变换的数值方法，求解出了多层粘弹性体的时域解。最后，文中还给出了路面动弯沉的计算结果与实测结果的对比来证明推导结果的准确性。     

Application of an improved boundary element method on the problem of elastic torsion

An improved boundary clement method has been used in analyzing and calculating the problems of the torsion of a prismatic bar with elliptical cross-section. In this paper the calculated results correspond with the values of boundary element method. However, the quantity of data required by the improved boundary element method is much less than that required by boundary element method, and the calculating time will be greatly reduced. Therefore, the procedure of this paper is an economical and efficient numerical computational way for solving Poisson equation problem.     

A 2D boundary element method for simulating the deformation of axisymmetric compound non‐Newtonian drops

The boundary integral formulation of the solution to the Stokes equations is used to describe the deformation of small compound non‐Newtonian axisymmetric drops suspended in a Newtonian fluid that is subjected to an axisymmetric flow field. The non‐Newtonian stress is treated as a source term in the Stokes equations, which yields an extra integral over the domains containing non‐Newtonian material. By transforming the integral representation for the velocity to cylindrical co‐ordinates and performing the integration over the azimuthal direction analytically, the dimension of the problem can be reduced from three to two. A boundary element method for the remaining two‐dimensional problem aimed at the simulation of the deformation of such axisymmetric compound non‐Newtonian drops is developed. Apart from a numerical validation of the method, simulation results for a drop consisting of an Oldroyd‐B fluid and a viscoelastic material are presented. Moreover, the method is extended to compound drops that are composed of a viscous inner core encapsulated by a viscoelastic material. The simulation results for these drops are verified against theoretical results from literature. Moreover, it is shown that the method can be used to identify the dominant break‐up mechanism of compound drops in relation to the specific non‐Newtonian character of the membrane. Copyright © 1999 John Wiley & Sons, Ltd.     

Simulation of three-dimensional polymer mould filling processes using a pseudo-concentration method

Mould filling processes, in which a material flow front advances through a mould, are typical examples of moving boundary problems. The moving boundary is accompanied by a moving contact line at the mould walls causing, from a macroscopic modelling viewpoint, a stress singularity. In order to be able to simulate such processes, the moving boundary and moving contact line problem must be overcome. A numerical model for both two- and three-dimensional mould filling simulations has been developed. It employs a pseudo-concentration method in order to avoid elaborate three-dimensional remeshing, and has been implemented in a finite element program. The moving contact line problem has been overcome by employing a Robin boundary condition at the mould walls, which can be turned into a Dirichlet (no-slip) or a Neumann (free-slip) boundary condition depending on the local pseudo-concentration. Simulation results for two-dimensional test cases demonstrate the model's ability to deal with flow phenomena such as fountain flow and flow in bifurcations. The method is by no means limited to two-dimensional flows, as is shown by a pilot simulation for a simple three-dimensional mould. The reverse problem of mould filling is the displacement of a viscous fluid in a tube by a less viscous fluid, which has had considerable attention since the 1960's. Simulation results for this problem are in good agreement with results from the literature. © 1998 John Wiley & Sons, Ltd.     

A practical method for the evaluation of eigenfunctions from compound matrix variables in finite elastic bifurcation problems

We show how the compound matrix method can be extended to give eigenfunctions as well as generalised eigenvalues to bifurcation problems in non-linear elasticity. When the incremental problem is formulated in terms of displacements only there are significant difficulties that arise from the non-trivial boundary conditions. In order to avoid these problems we adopt a Stroh formulation of the incremental problem. This then produces trivial boundary conditions for the compound matrix eigenvalue problem and more importantly known initial conditions for the compound matrix eigenfunction problem. This results in a straightforward and robust calculation for the eigenfunctions.     

求解二点边值问题打靶法的一种改进方法

本文讨论两点边值问题的数值解法,利用最小二乘法修正打靶法所需初始参数,将边值问题转化为相应的初值问题与左梯度控制方程合并,引入精细积分法进行求解。通过数值算例表明,该方法是一种求解两点边值问题的有效方法。     

Characteristics‐based boundary conditions for the Euler adjoint problem

Over the last decade, the adjoint method has been consolidated as one of the most versatile and successful tools for aerodynamic design. It has become a research area on its own, spawning a large variety of applications and a prolific literature. Yet, some relevant aspects of the method remain relatively less explored in the literature. Such is the case with the adjoint boundary problem. In particular for Euler flows, both fluid dynamic and adjoint equations entail complementary Riemann problems, and these yield boundary conditions that are fully consistent with well‐posedness. In the literature, this approach has been pursued solely in terms of Riemann variables. This work formulates the adjoint boundary problem so as to correspond precisely to that imposed on the flow, as it is given in terms of primitive variables. Test results have shown to be in agreement with the traditional approach for external flow problems. Copyright © 2012 John Wiley & Sons, Ltd.     

二维边界层方程的迭代求解

针对二维边界层方程,提出了分析积分迭代法.首先将该方法应用于Blasius方程和Falkner-Skan方程的求解,数值计算结果稳定,计算精度高;然后对外部有势流不能达到自相似要求、必须二维求解的二维层流边界层问题,在分析积分迭代法中加上计算力学的松弛迭代法,形成了一套有效的算法.数值结果表明该方法用于层流边界层的计算是很有效的.     

有限长轴承非稳态油膜力自由-移动边界问题的有限条解法

讨论了不可压缩流体润滑的动载径向滑动轴承油膜压力分布的自由移动边界问题的有限条解法．将自由边界问题转化为全域(矩形域)的具有不等式约束的微分方程的边值问题，进一步化为具有不等式约束的泛函极值问题。借助有限条法在矩形域上离散这个泛函，得到了一个特殊的二次泛函的规化问题。通过变量平移变换，使其化为标准的二次规划问题。然后借助于牛顿非光滑算法，迭代求解非线性的互补方程。给出了有限长轴承真实的油膜应力分布。对于所求解方程的系数矩阵的高度稀疏性。给出了紧缩存储算法。节省了存储空间和减少了计算量。算例表明该方法是有效的。     

一种改进格林元方法及在渗流问题中的应用

采用格林公式和基本解推导出直接边界积分方程来求解渗流问题.边界积分方程数值离散基于格林元方法(Green element methond),改进了原方法中压力和压力导数的求解方法,命名为混合边界元方法(Mixed boundary element method).相较于格林元类方法,该方法显式考虑了求解节点的外法向流量值和压力值,并使求得的数值解在求解区域上能够连续,符合实际的物理过程,在不增加额外未知数的情况下提高了计算精度.分析了不同网格类型对模拟计算结果的影响,并对稳定渗流问题、非稳定(瞬态)渗流问题和非稳态问题进行了实例计算,结果显示改进方法提高了计算精度,并对各类渗流问题有较好的适应性.     

插值型维数分裂无网格法中本质边界条件施加方法研究

本文以求解三维波动方程为例,介绍了改进的插值型维数分裂无单元Galerkin方法,推导了方程的弱形式,构造了具有插值特性的逼近函数,建立了可直接施加本质边界条件的离散方程组,研究不同本质边界条件施加方法对计算结果的影响．本文列举了三种常用的处理本质边界条件的方法：直接配点法、对角元素置大数法和对角元素化一法．选取了三个数值算例,分别采用不同的本质边界条件施加方法,分析计算结果,证明了三种施加方法的有效性,讨论了每种施加方法的优缺点,并针对问题需求选出合适的施加本质边界条件的方法．与改进的无单元Galerkin方法相比,改进的插值型维数分裂无单元Galerkin方法具有更高的计算精度和更快的计算速度．