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

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

一类新的具有无外力圆柱表面的杂交应力元

根据Hellinger-Reissner原理,建立了进一步改进的具有一个无外力圆柱表面三维杂交应力元.元内假定应力场满足以柱坐标表示的平衡方程,及圆柱面上的无外力边界条件.当退化为二维时,也满足协调方程.数值算例表明,当分析带圆弧的槽孔板、块时,在稀疏的有限元网格下,这类单元即可提供较以前各类特殊元、一般假定位移元及一般假定应力元远为准确的三维及二维应力分布.     

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

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

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

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

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

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

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

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

一种优化的中厚扁壳杂交元

自从产生第一个杂交应力元以来,杂交元理论不断发展完善,其应用范围也日益广泛。近年来,人们又采用广义变分原理作基础,在单元内以位移为Lagrange乘子放松单元应力平衡方程。由于应力不需要事先满足平衡条件,从而可采用自然坐标系,使得单元具有坐标不变性。     

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

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

偶应力理论的无网格法

根据偶应力理论的基本方程,在能量泛函中引入罚函数来满足偶应力理论的几何约束条件,从而得到偶应力理论的无网格法的控制方程.通过算例,分析了偶应力对带中心圆孔的无限平板在单轴拉伸及纯剪状态下的应力集中的影响,验证了该无网格方法的有效性和可行性.     

两端固支叠层圆柱厚壳轴对称问题的精确解析解

基于柱坐标系下的三维弹性力学基本方程,采用状态空间法得到两端固支单层与叠层圆柱厚壳轴对称问题的精确解析解。为严格满足固支端的边界条件,将固支端的边界位移函数作为状态变量引入状态方程,采用增维方法把非齐次状态方程变为齐次状态方程,并通过层合渐近技术将变系数状态矩阵转为常系数矩阵进行求解。所得到的解不仅严格满足三维弹性力学基本方程,而且严格满足固支边界条件,是真正意义上的三维精确解。算例表明,本研究解与有限元解吻合,具有很高的精度,且关于级数项数和分层数具有很好的收敛性。另外,通过圆柱厚壳各力学量沿径向和轴向的精确分布规律分析了厚径比和跨径比变化对位移和应力分布的影响。     

Coupling of assumed stress finite element and boundary element methods with stress-traction equilibrium

In this study, the stress based finite element method is coupled with the boundary element method in two different ways. In the first one, the ordinary distribution matrix is used for coupling. In the second one, the stress traction equilibrium is used at the interface line of both regions as a new coupling process. This new coupling procedure is presented without a distribution matrix. Several case studies are solved for the validation of the developed coupling procedure. The results of case studies are compared with the distribution matrix coupling, displacement based finite element method, assumed stress finite element method, boundary element method, ANSYS and analytical results whenever possible. It is shown that the coupling of the stress traction equilibrium with assumed stress finite elements gives as accurate results as those by the distribution matrix coupling.     

Free-edge interlaminar stress analysis of piezo-bonded composite laminates under symmetric electric excitation

A stress function-based approach is proposed to analyze the free-edge interlaminar stresses of piezo-bonded symmetric laminates. The proposed method satisfies the traction free boundary conditions, as well as surface free conditions. The symmetric laminated structure was excited under electric fields that can generate induced strain, resulting in pure extension in the laminated plate. The governing equations were obtained by taking the principle of complementary virtual work. To verify the proposed method, cross-ply, angle-ply and quasi-isotropic laminates were analyzed. The stress concentrations predicted by the present method were compared with those analyzed by the finite element method. The results show that the stress function-based analysis of piezo-bonded laminated composite structures is an efficient and accurate method for the initial design stage of piezo-bonded composite structures.     

THE STRESS SUBSPACE OF HYBRID STRESS ELEMENT AND THE DIAGONALIZATION METHOD FOR FLEXIBILITY MATRIX H

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓ,ｔｈｅｄｉｓｐｌａｃｅｍｅｎｔｅｌｅｍｅｎｔｓａｒｅｖｅｒｙｐｏｐｕｌａｒｂｅｃａｕｓｅｏｆｔｈｅｉｒｅａｓｙｃｏｎｓｔｒｕｃｔｉｏｎ .Ｈｏｗｅｖｅｒ,ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏａｓｓｕｍｅｔｈｅｃｏｎｔｉｎｕｏｕｓｄｉｓｐｌａｃｅｍｅｎｔｆｉｅｌｄｆｏｒｓａｔｉｓｆａｃｔｏｒｙｓｈｅｌｌｅｌｅｍｅｎｔｓｔｈａｔｉｓｆｒｅｅｆｒｏｍｓｈｅａｒｌｏｃｋｉｎｇ .Ｉｎ 1 964 ,ｔｈｅｈｙｂｒｉｄｓｔｒｅｓｓｅｌｅｍｅｎｔｗ…     

Orthogonal basic deformation mode method for zero-energy mode suppression of hybrid stress elements

A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The resulting orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. Moreover, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.     

Analysis of bonded anisotropic wedges with interface crack under anti-plane shear loading

The antiplane stress analysis of two anisotropic finite wedges with arbitrary radii and apex angles that are bonded together along a common edge is investigated. The wedge radial boundaries can be subjected to displacement-displacement boundary condi- tions, and the circular boundary of the wedge is free from any traction. The new finite complex transforms are employed to solve the problem. These finite complex transforms have complex analogies to both kinds of standard finite Mellin transforms. The traction free condition on the crack faces is expressed as a singular integral equation by using the exact analytical method. The explicit terms for the strength of singularity are extracted, showing the dependence of the order of the stress singularity on the wedge angle, material constants, and boundary conditions. A numerical method is used for solving the resul- tant singular integral equations. The displacement boundary condition may be a general term of the Taylor series expansion for the displacement prescribed on the radial edge of the wedge. Thus, the analysis of every kind of displacement boundary conditions can be obtained by the achieved results from the foregoing general displacement boundary condition. The obtained stress intensity factors （SIFs） at the crack tips are plotted and compared with those obtained by the finite element analysis （FEA）.     

杂交元本征应力模式和应力子空间的性质研究

详细讨论了有限元本征应力模式和应力子空间的性质，并着重讨论和进一步完善了与杂交应力有限元应力子空间有关的一些定理，为提出新方法提供了理论基础，主要包括：（1）证明了杂交元特征值不大于对应位移元的特征值；（2）证明了矩阵H非奇异的充分必要条件是假设应力模式线性无关；（3）证明了杂交元所对应位移元的本征应力模式形成的杂交元与该位移元相同；（4）证明了等价假设应力模式形成相同的杂交元；（5）证明了确定杂交元本征应力模式的充分必要条件是其范数平方等于所形成杂交元的变形模态特征值；（6）证明了杂交元假设应力模式与变形模态的能量一一对应的充分必要条件是假设应力模式彼此正交且与所对应位移元的本征应力模式除了一一对应者之外都正交。     

杂交元零能模式抑制的改进基本变形模式方法

基本变形模式方法是构造杂交元假设应力场的一种有效方法,但只给出了识别零能模式的必要条件.论文指出基本变形模式方法由于没有考虑耦合影响而存在应用上的局限,在此基础上利用柔度矩阵的正定性给出了判别零能机动模式的充要条件,并且通过内积向量无关性进一步完善了假设杂交元应力场的基本变形模式方法,可以适用于基本变形模式之间耦合的任...     

Extended Kantorovich method for local stresses in composite laminates upon polynomial stress functions

The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work, the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is estab-lished to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method (FEM). The convergent stresses have good agreements with those results obtained by three dimensional (3D) FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kan-torovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.     

The stress subspace of hybrid stress element and the diagonalization method for flexibility matrixH

The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix;2. ) The equivalent assumed stress modes lead to the identical hybrid element The Hilbert stress subspace of the assumed stress modes is established So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt’s method Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency is improved greatly The numerical examples show that the method is effective.