变宽度板在拉伸和弯曲作用下的三维应力集中

根据一种修正的余能原理，建立了具有一个无外力圆柱面的三维杂交应力元，元内假定应力场满足三维柱坐标表示的平衡方程及无外力圆柱面上的外力边界条件；当元退化为二维时，也满足协调条件。单元位移场选择与相邻单元协调。数值算例表明这种特殊杂交应力元在相当粗的网格下，能十分有效地分析变宽度薄／厚板在拉伸与弯曲作用下的三维(及二维)应力集中。     

两个相等圆孔板弯曲时的三维应力集中

根据Hellinger-Reissner原理建立了具有一个无外力圆柱表面的三维八节点杂交应力元,其假设应力场严格满足柱坐标下三维平衡方程及圆柱面上无外力边界条件;当元退化为二维时也满足协调方程。数值算例表明,这种特殊杂交应力元可高效地分析具有两个圆孔薄板和厚板的应力集中,特别是三维应力集中。     

含一个无外力圆柱面带转动的三维杂交应力元

根据一种修正的余能原理,建立了一类具有一个无外力圆柱表面及结点含转动自由度的8 结点新型三维杂交应力元. 单元边界位移场选择二次位移插值函数,且与相邻元协调;单元内假定应力场满足以柱坐标表示的平衡方程及圆柱面上无外力边界条件. 数值算例表明,这种特殊杂交应力元在相当粗的网格下即能十分准确地分析圆弧形槽口附近及曲梁的三维（及二维）的孔边应力分布.      

含一个无外力圆柱面带转动的三维杂交应力元简

根据一种修正的余能原理,建立了一类具有一个无外力圆柱表面及结点含转动自由度的8 结点新型三维杂交应力元. 单元边界位移场选择二次位移插值函数,且与相邻元协调;单元内假定应力场满足以柱坐标表示的平衡方程及圆柱面上无外力边界条件. 数值算例表明,这种特殊杂交应力元在相当粗的网格下即能十分准确地分析圆弧形槽口附近及曲梁的三维（及二维）的孔边应力分布.     

具有均匀受压圆柱面的三维杂交应力元

1.前言 应用具有圆柱表面的特殊杂交应力元分析具有圆柱面的物体时,将给出较一般的假定位移元、及一般的杂交应力元远为准确的结果。本文研究圆柱面上承受均匀法向力时的三维杂交应力元,它将用于更准确地分析一些弹性力学不易解决的压力管道等问题。     

扩展的Hellinger-Reissner原理及特殊层合杂交应力元

建立了一种非匀质材料新的、扩展的Hellinger-Reissner原理,发展了当一个单元域划分为不同材料特性子域、其元内应力场沿子域表面不连续、且位移场在子域表面也急剧变化时,一个非匀质有限元刚度列式便利方法。这种列式亦可用于对每层横向剪应变均独立处置的厚层板。基于此变分原理建立了新的具有一个无外力圆柱表面的层合杂交应力元,单元各层独立假定的应力场通过以自然坐标表示的非协调位移为权函数使齐次平衡方程变分满足的理性方法及严格满足给定圆柱面上无外力条件得到,位移场在元间及层间连续条件则分别通过Lagrange乘子进行了松弛。数值算例表明：这类新型元可有效地分析具有多类圆柱形槽孔的厚、中厚及薄层板自由孔边应力分布。     

岩土的一种强度准则及其变换应力法

基于岩土摩擦性,假设岩土破坏是由其物理空间内特征面上的应力比决定,提出了等效应力比的概念,即物理空间特征面上的剪应力合力与正应力合力的比值.在二维条件下,等效应力比可表示为σ-τ坐标系下与摩尔圆相切的直线扣除截距正切值;在三维条件下,假设在XYZ空间内存在一三维物理空间平面,此三维空间特征平面的等效应力比为影响材料强度特性的决定性因素,基于上述三维空间特征面建立了强度准则并称之为a准则.SMP准则以及广义Mises准则都是a准则的特例,当二维坐标中的截距为零时,则强度准则退化为SMP(spatially mobilized plane)强度准则,而当正切角为零时,则强度准则退化为广义Mises准则.而当截距与外切角均不为零时,则强度准则为介于上述两者之间的一种强度准则,在偏平面上为介于SMP曲边三角形与广义Mises圆形之间的曲边三角形.在子午面上,采用考虑岩土压剪耦合的屈服准则,破坏准则采用幂函数表达式.在偏平面上提出了基于a准则的形状函数,并采用真三维应力状态表示的破坏强度准则表示在三轴压缩路径下以p,q二维应力变量表达的准则公式,推导得到了基于a准则的变换应力公式,可简单地将一般以p,q为基本变量的二维模型转变为三维应力模型.通过强度以及多种应力路径的测试对比,验证了a准则及基于该准则的变换应力公式的合理性.     

新的一个求解带孔薄板弹性问题的二维杂交应力单元

建立了一个新的求解带圆孔薄板弹性问题的二维杂交应力单元,该单元为四节点四边形平面单元,名为P-HS4-8β。由极坐标系下的物理方程和几何方程求解出了一个极坐标方向的应力,通过将这个应力带入由Hellinger-Reissner原理推导的极坐标系下平面应力问题的能量方程中,得到了消除了该应力的能量方程,基于这个能量方程建立了杂交应力单元列式。根据圆孔边无外力条件和相容方程,推导了适用于求解带圆孔薄板问题的极坐标系下的二应力插值矩阵,并将此矩阵应用于新的有限单元列式中。数值算例表明新单元在求解孔边附近的应力时具有较高的精度。     

考虑膨胀力的非饱和介质热-水-应力耦合二维有限元分析

从建立应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,使用Galerk in方法,将各控制方程分别在空间域和时间域进行离散,开发出了一个可考虑膨胀力的用于分析非饱和介质中热-水-应力耦合弹塑性问题的二维有限元程序.通过对一个假定的核废料地下处置库的热-水-应力耦合问题的数值计算,比较了无、有膨胀力时的情况,在定性上验证了该程序的正确性.     

多层复合壳体三维振动分析的谱——微分求积混合法

对于较厚的多层复合壳体,其振动位移沿厚度方向呈锯齿形变化且层间剪切和拉、压应力呈三维耦合状态,采用传统的等效单层理论分析已不能满足精度要求. 建立不受结构厚度、铺层材料性质和铺层方式限制的三维分析方法具有重要的研究价值. 本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种适用一般边界条件和铺层方式的多层复合壳体三维分析新方法——谱--微分求积混合法. 该方法应用三维弹性理论对独立铺层进行精确建模,有效克服了二维简化理论对横向变形以及层间应力估计不确切的缺点;引入微分求积技术对铺层进行数值离散,将三维偏微分问题转化为二维偏微分问题,降低了求解维度和难度;应用广义谱方法近似地表述离散计算面上的场变量,将获取的二维偏微分方程转化为以场变量谱展开系数为未知量的线性代数方程组,避免了对超越方程的求解. 数值验证结果表明该方法收敛性好,计算精度高.      

A quadrilateral membrane hybrid stress element with drilling degrees of freedom

A new kind of quadrilateral assumed stress hybrid membrane element with drilling degrees of freedom and a traction-free inclined side has been developed based on an extended Hellinger-Reissner principle which is established by expanding the essential terms of the assumed stress field as polynomials in the natural coordinates of the element. The homogeneous equilibrium equations are imposed in a variational sense through the internal displacements which are also expanded in the natural coordinates, while the tractionfree conditions along the inclined side are satisfied exactly. The use of such special element in the finite element solution is shown to be highly accurate when only a very coarse element mesh is used for plates with V-shaped rounded notches or inclined sides.     

DEFORMATION RIGIDITY OF ASSUMED STRESS MODES IN HYBRID ELEMENTS

The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And the formula of the additional element deformation rigidity due to additional mode into the assumed stress field is derived. Based on, it is concluded in theory that the zero-energy stress mode cannot suppress the zero-energy deformation modes but increase the extra rigidity to the nonzero-energy deformation modes of the element instead. So they should not be employed to assume the stress field. In addition, the parasitic stress modes will produce the spurious parasitic energy and result the element behaving over rigidity. Thus, they should not be used into the assumed stress field even though they can suppress the zero-energy deformation modes of the element. The numerical examples show the performance of the elements including the zero-energy stress modes or the parasitic stress modes.     

含任意椭圆核各向异性板杂交应力有限元

基于含椭圆核有限大各向异性板弹性问题的复变函数级数解,应用杂交变分原理建立了一种与常规有限元相协调的含任意椭圆核各向异性板杂交应力有限元.单元内的应力场和位移场采用满足平衡方程、几何方程与物理方程的复变函数级数解,假设的复变函数级数解精确满足椭圆核边界处的位移协调条件和应力连续条件,单元外边界上的位移场按常规有限元位移场假设,单元内椭圆核的长轴可以与材料主轴不重合.单元刚度矩阵采用Gauss积分求得,并给出了建立刚度矩阵的主要公式和推倒过程.数值计算结果表明该单元具有计算精度高、计算工作量小等优点.     

Axisymmetric solid-of-revolution elements based on the assumed stress hybrid model

A series of axisymmetric solid-of-revolution elements with 8-node and quadrilateral cross section have been developed based on the assumed stress hybrid model. A quadratic boundary displacement assumption is employed for each element and a variety of interior stress assumptions have been made. Two different kinds of procedure used for developing stress field have been studied. Example problems of a thick-walled cylinder under internal pressure and a thick-walled sphere under internal pressure are utilized to evaluate the various elements and a desirable stress assumption has been identified. Comparisons of present results with those obtained by the use of 8-node element based on the assumed displacement model indicate that this hybrid stress element is far superior in predicting the stress distribution.     

圆柱型各向异性弹性力学平面问题

本文对圆柱型各向异性弹性力学平面问题的基本方程进行了改写。在此基础上，导出了应力函数Ｇ和位移函数φ，它们满足相同的控制方程，比文〔１〕的应力函数Ｆ的控制方程要简单，便于求得特解，并有Ｆ＝ｒＧ的关系。还对若干经典问题进行了求解。