用有限元广义混合法分析不可压缩或几乎不可压缩弹性体

不可压缩或几乎不可压缩问题在数学上表现为最小 势能原理中的某些项趋于无穷大,使得有限元方程产生病态。本文给出了不可压缩或几乎不可压缩弹性分析的广义混合变分原理,以此为基础建立了该类问题的有限元广义混合法。该变分原理的泛函中不含有上面这种奇异项,故其有限元方程不会产生病态。算例表明该有限元法可以同时进行可压缩、不可压缩或几乎不可压缩弹性分析,且精度良好;有限元常规位移法及Hermann法是该法的特例。     

含压电片层合壳的有限元分析与控制仿真

本文首先针对含压电材料的一般结构，推导了按位移和电位移求解时的混合变分原理，在此基础上，通过对面电荷取变分，直接得到了含按电压驱动压电片层合壳的有限元方程，最后给出了利用压电片进行静变形和振动控制的仿真算例，并进行了分析。结果表明，本文建立的有限元方程准确、可靠。     

结构动态应力计算方法研究

复杂结构动态应力的准确计算是一个没有圆满解决的问题。本文以最小余能原理为基础，提出了计算结构动态应力的最小伤痛有法。该方法采用二次分析思想，首先采用常规有限元对结构进行适当离散，计算输出结构所需应力区域的有限元结点位移和加速度动力时程反应，再应用最小余能法计算所求部位的动态应力值。这种方法的优点是它可以与现有的有限元程序有机结合，方便使用；动应力在区域内的分布规律可以由计算者根据具体情况而确定，一般情况下，可以选用二次曲线来逼近动应力在区域内的实际分布，避免了常规有限元法计算结构动应力时必须对单元形函数求导的做法，从而提高了动应力计算精度。计算结果表明：本文方法计算结构动态应力结果较常规有限元法的计算结果有明显改进，特别是当结构变化剧烈时，改进效果更为明显。     

On parametric minimum complementary energy variational principle in soil mechanics

The parametric minimum complementary energy variational principle is given in this paper. To the problems of nonassociative flow rule in the theory of plasticity, the Drucker postulation can no longer be applied and the classical variational principles fail, the parametric variation principles can play the role instead. The parametric variational principle can be used for the finite element solution with sequential quadratic programming method to soil mechanics problems.     

Finite deformation plate theory and large rigid-body motions

A refined non-linear first-order theory of multilayered anisotropic plates undergoing finite deformations is elaborated. The effects of the transverse shear and transverse normal strains, and laminated anisotropic material response are included. On the basis of this theory, a simple and efficient finite element model in conjunction with the total Lagrangian formulation and Newton-Raphson method is developed. The precise representation of large rigid-body motions in the displacement patterns of the proposed plate elements is also considered. This consideration requires the development of the strain-displacement equations of the finite deformation plate theory with regard to their consistency with the arbitrarily large rigid-body motions. The fundamental unknowns consist of six displacements and 11 strains of the face planes of the plate, and 11 stress resultants. The element characteristic arrays are obtained by using the Hu-Washizu mixed variational principle. To demonstrate the accuracy and efficiency of this formulation and compare its performance with other non-linear finite element models reported in the literature, extensive numerical studies are presented.     

A variational solution of coupled vibration for moderately thick floating plate seated on the sea floor with variable water depth

Based on Reissner's plate theory, a localized variational principle for the coupled vibrations of moderately thick floating plates seated on the sea floor with variable water depth is presented in this paper. It will be a powerful tool to formulate various approximate methods of moderately thick floating plates. By means of this variational principle a formula for calculating the natural frequencies of moderately thick floating plates is also obtained. A numerical example shows that high accuracy and computational time saving may be achieved in calculation.     

变水深环境下中厚度浮板耦振问题的一个变分解

本文从Reissner 板的理论出发,提出并论证了变水深环境下中厚度浮板耦振问题的一个局部变分原理,它是建立中厚度浮板各种近似分析方法的一个有力工具,借助这个变分原理,文中还导得了一个计算中厚度浮板固有频率的一个变分式,数值算例表明本文方法具有精度高和省机时的优点。     

some general theorems and generalized and piecewise generalized variational principles for linear elastodynamics

From the concept of four-dimensional space and under the four kinds of time limit conditions, some general theorems for elastodynamics are developed, such as the principle of possible work action, the virtual displacement principle, the virtual stress-momentum principle, the reciprocal theorems and the related theorems of time terminal conditions derived from it. The variational principles of potential energy action and complementary energy action, the H-W principles, the H-R principles and the constitutive variational principles for elastodynamics are obtained. Hamilton's principle, Toupin's work and the formulations of Ref. [5], [17]-[24] may be regarded as some special cases of the general principles given in the paper. By considering three cases: piecewise space-time domain, piecewise space domain, piecewise time domain, the piecewise variational principles including the potential, the complementary and the mixed energy action fashions are given. Finally, the general formulation of piecewise variati     

含多椭圆孔(核)复合材料加筋层板的应力集中研究

基于各向异性体平面弹性理论中的复势方法,应用杂交变分原理建立了一种与常规有限元相协调的含任意椭圆核各向异性板杂交应力有限元,采用该杂交应力有限元来描述层板的椭圆核区域,采用杆单元来描述加强筋(杆单元的刚度取为层板沿筋条方向的刚度),其余区域采用常规8节点等参单元进行模拟,建立起分析含多椭圆核复合材料加筋壁板问题的力学分析方法,详细讨论了椭圆核大小、位置、筋条尺寸、相对位置、铺层比例等诸参数的影响规律,得到了一些有益的结论。     

一种高效率弹塑性有限元分析新方法：增量杂交／混合修正弦线模量法

本文在文献[1]给出的放松应力增量平衡约束的修正余能广义变分原理基础上,提出一种高效率的弹塑性有限元分析的新方法——增量杂交/混合修正弦线模量法。该法保持了文[1]方法的全部优点,而在迭代过程中,依据材料的单向拉伸应力—应变关系,不断改变过渡区和塑性区单元柔度矩阵和塑性矩阵中的弹性模量;并在体积压缩模量不变假设下,相应地改变过渡区单元矩阵中的泊松系数。从而大大降低了迭代收敛次数和单刚计算量,提高了多类变量弹塑性有限元分析的计算效率和收敛精度。     

On some general variational principles for creep with applications to thin shells

A constitutive relation for a viscous material subject to small strains but finite rotations is postulated and associated variational theorems are formulated. These are similar to the principle of minimum of potential energy and the Hellinger-Reissner theorem of an elastic solid. The derivation of strain-displacement relations for thin shells subject to small strains but moderately large rotations are given. On this basis a mixed variational principle for thin viscous shells is developed. For the problem of creep collapse of long cylindrical shells under external pressure it is demonstrated that the mixed variational principle may be advantageous compared to other variational theorems. A comparison with the creep collapse theory of Hoff et al. is given.     

The combined energy variational principles in linear elasticity

In this paper, a new kind of mixed energy variational principles in linear elasticity—the combined energy variational principles is presented. First, we define the conjugate body of an elastic body, which is obtained by changing the boundary conditions of the elastic body. Next, we decompose the conjugate body into two component-states, construct functionals of potential energy and complementary energy, respectively, for the component-states and define the additional hybrid energy between the component-states. Thus the functionals of combined energy can be constructed. Three typical combined energy variational principles are demonstrated and the other forms of combined energy variational principles are given. The application of the proposed principles to the calculation of thin plates with complicated boundaries is shown.     

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

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

Variational principles for Cosserat body

Two variational principles are derived for the mixed boundary value problem of Cosserat solid. These principles are a generalization of the stationary principle of potential energy and the stationary principle of complementary energy from non-linear theory of elasticity.     

厚板弹塑性混合状态Hamiltonian元的半解析法

采用初应力法并通过对增量形式Ｈｅｌｌｉｎｇｅｒ－Ｒｅｉｓｓｎｅｒ变分原理的修正，给出了厚板弹塑性分析的Ｈａｍｉｌｔｏｎ正则方程及其相应的半解析法。此法在厚板平面内采用通常的有限元离散，而沿板厚方向采用状态空间法给出半解析解答。     

Modified H-R mixed variational principle for magnetoelectroelastic bodies and state-vector equation

Based upon the Hellinger-Reissner (H-R) mixed variational principle for three-dimensional elastic bodies, the modified H-R mixed variational theorem for magnetoelectroelastic bodies was established. The state-vector equation of magnetoelectroelastic plates was derived from the proposed theorem by performing the variational operations. To lay a theoretical basis of the semi-analytical solution applied with the magnetoelectroelastic plates, the state-vector equation for the discrete element in plane was proposed through the use of the proposed principle. Finally, it is pointed out that the modified H-R mixed variational principle for pure elastic, single piezoelectric or single piezomagnetic bodies are the special cases of the present variational theorem.     

Canonical dual finite element method for solving post-buckling problems of a large deformation elastic beam

This paper presents a canonical dual mixed finite element method for the post-buckling analysis of planar beams with large elastic deformations. The mathematical beam model employed in the present work was introduced by Gao in 1996, and is governed by a fourth-order non-linear differential equation. The total potential energy associated with this model is a non-convex functional and can be used to study both the pre- and the post-buckling responses of the beams. Using the so-called canonical duality theory, this non-convex primal variational problem is transformed into a dual problem. In a proper feasible space, the dual variational problem corresponds to a globally concave maximization problem. A mixed finite element method involving both the transverse displacement field and the stress field as approximate element functions is derived from the dual variational problem and used to compute global optimal solutions. Numerical applications are illustrated by several problems with different boundary conditions.     

COMPUTATION OF STRESS INTENSITY FACTORS BY THE SUB-REGION MIXED FINITE ELEMENT METHOD OF LINES

Based on the sub-region generalized variational principle,a sub-region mixed ver- sion of the newly-developed semi-analytical‘finite element method of lines’(FEMOL)is pro- posed in this paper for accurate and efficient computation of stress intensity factors(SIFs)of two-dimensional notches/cracks.The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used,with the sought SIFs being among the unknown coefficients.The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements.A mixed system of ordinary differential equations(ODEs) and al- gebraic equations is derived via the sub-region generalized variational principle.A singularity removal technique that eliminates the stress parameters from the mixed equation system even- tually yields a standard FEMOL ODE system,the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver.A number of numerical examples,including bi-material notches/cracks in anti-plane and plane elasticity,are given to show the generally excellent performance of the proposed method.     

The method of transformation of the boundaries for structure the admissible displacements

In this paper the method of transformation of the boundaries for structure the admissible displacements with various boundary conditions is presented. What is called the method of transformation of the boundaries is that, first we transform the actual system into the basic system and additional boundary forces and displacements on the basis of the superposition principle, then apply variational principles to the basic system, finally find the permissible displacement of the actual system by means of the method of transformation of the series.In this paper, we also give mixed energy principles under variation of boundary conditions. The mixed energy principles as the potential and complementary energy principles under variation of boundary conditions are all the chief theoretical fundamental of the method of transformation of the boundaries.Applying the method of transformation of the boundaries, we form the permissible displacements of rectangular plates of plane stress and bending problems with various edge conditions.Because the method of transformation of the boundaries is in progress to follow the variational principles and definite program to form permissible displacements, the difficulty in supposing and piecing together permissible displacements in the Rayleigh-Ritz method will be overcome.     

基于余能原理的有限变形问题有限元列式

利用基面力概念,推导了一种基于余能原理的有限变形问题显式有限元列式,可应用于结构的大位移、大转动问题。以基面力为状态变量来表达单元的余能,将有限变形情况下的单元余能分解为变形余能部分和转动余能部分,利用Lagrange乘子法推导出余能原理有限元的控制方程,编制了相应的非线性有限元程序。通过算例分析,说明该列式和程序的可靠性和精确性。