一类带单侧约束力学问题的变分方法

力学中许多问题如弹塑性问题，弹塑性接触问题、塑性耦合问题，渗流问题等，其约束条件或者边界形态是无法预知的，只能用一组不等式表示，属于非线性问题，经典的变人原理在求解此类问题时只能采取繁琐的迭代过程。     

关于平面问题分类的再探讨

关于平面问题分类的再探讨张雷顺（郑州工学院，郑州４５０００２）由于平面问题比较简单，所以在弹性力学教材中，对平面问题给出了较多的篇幅，也列举了较多的例题．在工程中，相当多的问题也都近似为平面问题处理．但在认定平面问题和进行平面问题分类方面，作者发现时...     

关于平面问题分类的再探讨

关于平面问题分类的再探讨张雷顺（郑州工学院，郑州４５０００２）由于平面问题比较简单，所以在弹性力学教材中，对平面问题给出了较多的篇幅，也列举了较多的例题．在工程中，相当多的问题也都近似为平面问题处理．但在认定平面问题和进行平面问题分类方面，作者发现时...     

动态设计变量优化方法及其在工程静力学通用程序设计中的应用

本文提出动态设计变量优化新方法并解决理论和工程问题，该方法首先构建目标函数框架，然后根据具体问题的输入条件，动态地进行设计变量分配和合理排序，形成实际问题的动态目标函数。基于动态设计变量优化方法，编制出一个能解决单刚体、刚体系的平面和空间问题、摩擦问题和桁架所有工程静力学平衡问题的通用程序，并通过实例分析验证。为解决更多工程领域问题提供有效新观点。     

变厚度板弯曲的边界元分析法

由于变厚度板弯曲问题的控制分方程复杂，直接求解其基本解推导边界积分方程建立边界元分析法较为困难，本文通过引入等效荷载，等效刚度，将此问题的控制微分方程化成与普通薄板弯曲问题基本方程相同的形式，利用求解通板弯曲问题的边界元迭代求解，建立了分析变厚度板弯曲问题的蛤法，算例表明本方法理正确，精度良好。     

套管-水泥环-地层应力分布的理论解

利用弹性力学理论研究地应力场中套管-水泥环-地层系统应力分布的理论解，在求解过程中，将原问题分解为两个相对简单的子问题．考虑到套管井系统是层状结构的特点，且各层受力情况具有相似性，进一步地，基于两个弹性力学平面应变基本问题的解，采用结构力学的求解思想，由位移连续条件分别求得了两个子问题的解，最后由叠加原理得到了原问题的理论解．     

求解线性方程组的一种新方法

将线性方程组的一般系数矩阵转化为对称正定矩阵，从而把原线性方程组的求解问题转化为一个等价变分问题的极少值点寻优问题，借助对分寻优法进行求解。算例结果表明，本文方法不仅对于良态线性方程组的求解问题是有效的，而且对于病态线性方程组的求解问题同样是有效的。     

一维杆与结构接触冲击问题的计算

一维杆与结构接触冲击问题的计算黄剑敏，任文敏（清华大学工程力学系，北京１０００８４）１引言工程中大量存在着杆与结构的接触冲击问题，对这类问题的研究极为重要，但这方面的论著甚少，理论所能解决的问题更少，仅在极个别的情况下有精确解。Ｔ．Ｔ．Ｒ．Ｈｕｇｈｅ...     

一维杆与结构接触冲击问题的计算

一维杆与结构接触冲击问题的计算黄剑敏，任文敏（清华大学工程力学系，北京１０００８４）１引言工程中大量存在着杆与结构的接触冲击问题，对这类问题的研究极为重要，但这方面的论著甚少，理论所能解决的问题更少，仅在极个别的情况下有精确解。Ｔ．Ｔ．Ｒ．Ｈｕｇｈｅ...     

复合材料的裂纹动力学问题的模型

复合材料产生裂纹后，其纤维处形成“桥连”，这是一个不可避免的现象。由于桥连问题很复杂，在数学方法的处理上有很大困难，至今人们研究的大多是桥连的静力学问题，而对其动力学问题研究得很少。只有建立复合材料的桥连动力学模型，才能更好地研究复合材料的断裂动力学问题。为了便于分析复合材料的问题，将桥连处用载荷代替，当裂纹高速扩展时，其纤维也连续地断裂。通过复变函数论的方法，将所讨论的问题转化为Riemann-Hilbet问题。利用建立的动态模型和自相似方法，得到了正交异性体中扩展裂纹受运动的集中力Px/t及均布载荷作用下位移、应力和动态应力强度因子的解析解，并通过迭加原理，最终求得了该模型的解。     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题，推导了相应的离散方程. 与 其它基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不 需要划分网格，用配点法离散控制方程，可以直接施加边界条件，不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

A stability study of Navier-Stokes equations (III)

In this paper, the necessary conditions of the existence of C² solutions in some initial problems of Navier-Stokes equations are given, and examples of instability of initial value (at t=0) problems are also given. The initial value problem of Navier-Stokes equation is one of the most fundamental problem for this equation various authors studies this problem and contributed a number of results. J. Lerav, a French professor, proved the existence of Navier-Stokes equation under certain defined initial and boundary value conditions. In this paper, with certain rigorously defined key concepts, based upon the basic theory of J. Hadamard partial differential equations¹, gives a fundamental theory of instability of Navier-Stokes equations. Finally, many examples are given, proofs referring to Ref. [4].     

基于能量等效原理的应变局部化分析:Ⅱ.有限元解法

以非局部塑性理论为基础,应用状态空间理论,通过局部和非局部两个状态空间的塑性能量耗散率等效原理,提出了一种求解应变局部化问题的新方法,以得到与网格无关的数值解.针对二维问题的屈服函数和流动法则导出了求解非局部内变量的一般方程,并提出了在有限元环境中求解应变局部化问题的应力更新算法.为了验证所提出的方法,对1个一维拉杆和3个二维平面应变加载试件进行了有限元分析.数值结果表明,塑性应变的分布和载荷-位移曲线都随着网格的变小而稳定地收敛,应变局部化区域的尺寸只与材料内尺度有关,而对有限元网格的大小不敏感.对于一维问题,当有限元网格尺寸减小时,数值解收敛于解析解.对于二维剪切带局部化问题,数值解随着网格尺寸的减小而稳定地向唯一解收敛.当网格尺寸减小时,剪切带的宽度和方向基本上没有变化.而且得到的塑性应变分布和网格变形是平滑的.这说明,所提方法可以克服经典连续介质力学模型导致的网格相关性问题,从而获得具有物理意义的客观解.此模型只需要单元之间的位移插值函数具有C~0连续性,因而容易在现有的有限元程序中实现而无需对程序作大的修改.     

Analogy solutions of axisymmetrical twist problems by axisymmetrical finite element program

It used to be considered that an axisymmetrical problem and a twist problem of an axisymmetrical body cannot be simulated by each other, because the number of unknown variables in an axisymmetrical problem is greater than that in a twist problem, and the governing equations are not the same. This paper proposes a degenerated analogy method, by which the twist problems of axisymmetrical bodies can be simulated by axisymmetrical problems with finite element programs.An ordinary structural analysis method can be used to analyze an axisymmetrical problem, but a twist problem of axisymmetrical bodies is treated as a 3-dimensional problem usually. According to the method proposed in this paper, the analysis of a twist problem can be simulated by the analysis of an axisymmetrical body with a structural analysis problem. The example of analysis computation is also given. Thecomputed result is in agreement with the theoretical result.In this paper, the constitutive relation of the degenerated analogy problem is given.The authors suggest that a twist problem of a body made of any materials is simulated by an axisymmetrical problem of a body made of orthotropic material. If you have to use some program for the axisymmetrical problem to be limited to isotropic materials the penalty coefficient method can be used to solve the problem.     

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.     

Studies of quality of geometrically nonlinear elasticity theory for small strains and arbitrary displacements

Earlier it was shown in [1, 2] that the equations of classical nonlinear elasticity constructed for the case of small strains and arbitrary displacements are ill posed, because their use in specific problems may result in the appearance of “spurious” bifurcation points. A detailed analysis of these equations and the construction, in their stead, of consistent equations of geometrically nonlinear theory of elasticity can be found in [3]. Certain steps in this direction were also made in [4, 5]. In [3], it was also stated that the methods and applied program packages (APPs) based on the use of the classical relations of nonlinear elasticity require some revision and correction. In the present paper, this conclusion is justified and confirmed by numerical finite-element solutions of several three-dimensional geometrically nonlinear deformation problems and linearized problems on the stability of equilibrium of a rectilinear beam. These solutions were obtained by using two APPs developed by the authors and the well-known APP “ANSYS.” It is shown that the classical equations of the geometrically nonlinear theory of elasticity, which underly the first of the developed APP and the well-known APP “ANSYS,” often lead to overestimated buckling loads for structural members as compared with the consistent equations proposed in [1–3].     

变边界粘弹性轴对称问题的复变函数法

对边界几何形状、位置随时间变化的变边界结构,给出了用复变函数求解粘弹问题的解析方法。文中用拉普拉斯变换结合平面弹性复变方法,对内外边界变化时粘弹性轴对称问题进行求解。引入两个与时间、空间相关的解析函数,给出了变边界情况下应力、位移以及边界条件与解析函数的关系。当解析函数形式部分确定,则可用边界条件求解其中与时间相关的待定函数。求解待定函数的方程一般情况下为一系列积分方程,特殊情况可求得解析解。对轴对称问题中应力边值问题、位移边值问题以及混合边值问题,分别利用边界条件求得相关系数,从而得到了应力与位移的解析表达。当取Boltzmann粘弹模型时,进行不同边值问题的分析。分析显示,应力、位移的形态与大小均与边界变化过程相关,与固定边界粘弹性问题有较大不同。本文解答可用于粘弹性轴对称问题内外边界任意变化及各种边值问题的力学分析。此外,该法可进一步进行荷载非对称、复杂孔型变边界问题的求解。     

The Schrödinger equation in theory of plates and shells with orthorhombic anisotropy

This work is the continuation of the discussion of Refs. [1-5]. In this paper:[A] The Love-Kirchhoff equations of vibration problem with small deflection for orthorhombic misotropic thin shells or orthorhombic anisotropic thin plates on Winkler’s base are classified as several of the same solutions of Schrodmger equation, and we can obtain the general solutions for the two above-mentioned problems by the method in Refs. [1] and [3-5].[B] The. von Karman-Vlasov equations of large deflection problem for shallow shells with orthorhombic anisotropy (their special cases are the von Harmon equations of large deflection problem for thin plates with orthorhombic anisotropy) are classified as the solutions of AKNS equation or Dirac equation, and we can obtain the exact solutions for the two abovementioned problems by the inverse scattering method in Refs. [4-5].The general solution of small deflection problem or the exact solution of large deflection problem for the corrugated or rib-reinforced plates and shells as special cases is included in this paper.     

Construction of approximate solutions of boundary value impact strain problems

The method for constructing approximate solutions of boundary value problems of impact strain dynamics in the form of ray expansions behind the strain discontinuity fronts is generalized to the case of curvilinear and diverging rays. This proposed generalization is illustrated by an example of dynamics of an antiplane motion of an elastic medium. The ray method is one of the methods for constructing approximate solutions of nonstationary boundary value problems of strain dynamics. It was proposed in [1, 2] and then widely used in nonstationary problems of mathematical physics involving surfaces on which the desired function or its derivatives have discontinuities [3–7]. A complete, qualified survey of papers in this direction can be found in [8]. This method is based on the expansion of the solution in a Taylor-type series behind the moving discontinuity surface rather than in a neighborhood of a stationary point. The coefficients of this series are the jumps of the derivatives of the unknown functions, for which, as a consequence of the compatibility conditions, one can obtain ordinary differential equations, i.e., discontinuity damping equations. In the case where the problem with velocity discontinuity surfaces is considered in a nonlinear medium, this method cannot be used directly, because one cannot obtain the damping equation. A modification of this method for the purpose of using it to solve problems of that type was proposed in [9–11], where, as an example, the solutions of several one-dimensional problems were considered. In the present paper, we show how this method can be transferred to the case of multidimensional impact strain problems in which the geometry of the ray is not known in advance and the rays become curvilinear and diverging. By way of example, we consider a simple problem on the antiplane motion of a nonlinearly elastic incompressible medium.     

饱和地基上弹性圆板的动力响应

研究弹性圆板在饱和地基上的垂直振动特性,即首先应用Hankel变换方法求解饱和土波动方程,然后按混合边值条件建立饱和地基上圆板垂直振动的对偶积分方程,用一种简便的方法,对偶积分方程可化为易于数值计算的第二类Fredholm积分方程。文末的数值分析得出了板振动的一些规律性,由此表明当板的挠曲刚度D趋于无穷大且不计板的质量时,其结果和无质量刚性圆盘在饱和地基上的振动特性完全一致。