更一般的杂交广义变分原理及有限元模型

本文根据钱伟长教授提出的更一般的广义变分原理,给出了适用于有限元法中的更广义杂交变分原理,并由此建立了新的广义杂交模理. 进一步以变厚度薄板弯曲单元为例,对基于各种不同的广义变分原理建立的各种杂交元做了比较.     

基于参数变分原理的非均质材料弹塑性有限元分析的Voronoi单元法

在非均质材料的有限元数值模拟中,采用了Voronoi单元(VCFEM)以克服经典位移元的局限性．基于参数变分原理和二次规划法进行了Voronoi单元的二维弹塑性分析A·D2推导了有限元列式并形成最终的二次规划求解模型．研究了非均质材料微观夹杂对整体力学性能的影响．数值算例证明了该方法的正确和可行性．     

三角形单元协调与非协调位移的能量正交关系

基于组合稳定化变分原理,周天孝提出的组合杂交法是绝对收敛和稳定的,它给出了一种系统化的增强应力/应变方法,并建立了一簇低阶仿射等价的n-cube(n=2,3)单元。本文论证了单元上应力插值为线性,位移插值为协调线性部分和非协调二次部分之和的三角形组合杂交单元其协调部分与非协调部分的能量正交关系,进而得到此三角形单元刚度矩阵等同于协调的三角形线性元刚度矩阵,即非协调部分无应变增强特性。     

Stokes流动的罚-杂交有限元分析*

本文建议了一种用于分析Stokes流动的罚-杂交变分原理,其中,偏应力张量和静水压力事先满足线动量平衡.建立了相应的有限元模型.由此,压力可在列式过程中消去,使得有限元矩阵方程仅以节点速度作为唯一的求解未知量.推导了几种4-节点和8-节点四边形单元.通过数值算例,显示了单元性能.     

含多个任意参数的广义变分原理及换元乘子法

弹性力学变分原理的泛函变换可分为三种格式:Ⅰ、放松格式,Ⅱ、增广格式,Ⅲ、等价格式. 根据格式Ⅲ,提出含多个任意参数的广义变分原理及其泛函表示式,其中包括:以位移u为一类泛函变量的多参数广义变分原理;以位移u和应力σ为二类泛函变量的多参数广义变分原理;以位移u和应变ε为二类泛函变量的多参数广义变分原理;以位移u应变ε和应力σ为三类泛函变量的多参数广义变分原理.由这些原理可得出等价泛函一系列新形式,此外,通过参数的合理选择,可构造出一系列有限元模型. 本文还讨论了拉氏乘子法“失效”问题,指出“失效”现象产生的原因,提出乘子法“恢复有效”的作法——换元乘子法.     

基于广义变分原理的矩形薄板单元

本文应用广义变分原理,构造了适合正交各向异性薄板静动力分析的矩形单元MR—12.计算结果表明,基于广义变分原理的非协调元具有很好的收敛性和计算精度.证明了广义变分原理在建立非协调单元中的有效性和优越性.MR—12单元的计算格式和普通矩形板元无原则性的差别,极易推广使用.     

基于20节点辛元的复合材料层合板应力分析

通常情况下,常规位移有限元法获得的应力结果比位移精度低一阶次,且面外应力难以满足连续性要求.联合最小势能原理和H-R变分原理,构造出包含位移和3个面外应力两类变量的20节点六面体辛元.由于两类变量采用高阶插值函数近似,无需引入单元内部的非协调位移项,因此相关理论的推导过程非常简单.与Hamilton部分混合元不同,该辛元涉及的变量沿3个坐标方向均做离散处理,不受单元厚度和结构几何形状的限制.数值实例表明20节点辛元的数值结果收敛稳定.在粗糙网格的情况下,与20节点位移元相比,该文单元的面外应力更接近精确解.     

非线性弹性体的弹性动力学变分原理

本文根据文献[1],对非线性应力应变关系的弹性体,导出了弹性动力学问题的变分原理和广义变分原理,提出了混合位移协调元和混合应力协调元的瞬时广义变分原理.     

弹性理论中各种变分原理的分类

本文按弹性理论中各种变分原理的约束条件的不同,对所有变分原理进行分类.我们在前文中业已指出,应力应变关系这样的约束条件是不能用拉氏乘子法解除的.剩下的可能约束条件共有四种:(1)平衡方程,(2)应变位移关系,(3)边界外力已知的边界条件,和(4)边界位移已知的边界条件.弹性理论的各种变分原理中,有的只有一种约束条件,有的有两种或三种,最多只能有四种约束条件.这样一共可能有15种变分原理,但是每种变分原理既可以用应变能A表示,又可以用余能B表示.这样,我们一共应有30种形式完全不同的变分原理,我们全部列出了这三十种形式的变分原理.     

纯位移线弹性方程Locking-Free非协调三棱柱单元的构造分析

主要构造了三维空间中线弹性问题纯位移变分形式下无闭锁三棱柱单元.此单元是具有18个自由度的非协调元.单元的形函数满足位移的散度属于零次多项式空间,通过分析得到有限元解和真解误差的能量模具有一阶收敛性,L~2模具有二阶收敛性.     

电动力学电磁场边值问题的广义变分原理

给出了线性各项异性电磁场边值问题的广义虚功原理表达式,运用钱伟长教授提出的方法建立了该问题的广义变分原理,可直接反映该问题的全部特征,即4个Maxwell方程、2个场强-位势方程、2个本构方程和8个边界条件．继而导出了一族有先决条件的广义变分原理．作为例证,导出了两个退化形式的广义变分原理,和已知的广义变分原理等价．此外还导出了两个修正的广义变分原理,可为该问题提供杂交有限元模型．建立的各广义变分原理可为电磁场边值问题的有限元应用提供更为完善的理论基础．     

粘性流体动力学的混合协调元和混合杂交非协调元变分法

本文提出并论证了粘性正压流体流动的混合协调有限元变分原理,在论证中发现应力协调条件会自然满足。同时提出论证了有关的混合杂交非协调有限元广义变分原理,这能够简化粘性流动的非协调元计算。     

Variational formulation of the material failure process in solids by embedded discontinuities model

This work presents a variational formulation of the material failure process, idealized as strain or displacement discontinuities, by weak, strong, or discrete embedded discontinuities into a continuum. It is shown that the solution of the proposed variational formulation may be approximated by different types of finite elements with embedded discontinuities. The developed displacement approximation of a finite element split by the discontinuity leads to a symmetric stiffness matrix, which considers not only the continuity of tractions but also the rigid body relative motions of the portions in which the element is split. The variational formulation of a continuum with more than one discontinuity in its interior is developed. It is shown that this formulation may lead to finite elements with embedded discontinuities that can be classified as displacement, force, mixed, and hybrid models. To show the effectiveness of the proposed formulation, the classical example of a bar under tension is solved using one and 2D finite element approximations. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

饱和多孔介质耦合系统的变分原理

本文采用变积方法,建立了等温准静态下饱和多孔介质的六类变量的广义变分原理.在此基础上,通过引入约束条件得到各级变分原理,其中包括五类变量,四类变量,三类变量和二类变量的变分原理.除得到文献中已有的变分原理外,本文给出了许多新的变分原理,为建立饱和多孔介质的有限元模型提供了基础.     

Advanced Shell Element Formulations for Coupled Electromechanical Systems

With the significantly increasing applications of smart structures, piezoelectric material is widely used in branches of engineering sciences. Normally, the Finite Element Method is employed in the numerical analysis of these structures [2]. In this contribution, in order to avoid the locking effects and zero energy modes, the Assumed Natural Strain (ANS) Method [4] is implemented into four‐node piezoelectric shallow shell elements, by using the two‐field variational formulation in which displacements and electric potentials serve as independent variables and the three‐field variational formulation in which the dielectric displacement is taken as an independent variable additionally [3]. Moreover, a quadratic variation of the electric potential through the thickness direction is applied in the two‐field formulation. Numerical examples of piezoelectric sensors and actuators are presented, showing the behaviour of the shell elements by using different hybrid finite element formulations. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

弹性厚板的分区广义变分原理

本文提出弹性厚板分区广义变分原理,其要点如下:1.各分区可任意定为势能区或余能区.分区势能、分区余能、分区混合变分原理是它的三种特殊形式.2.每个分区中独立变分变量的个数可任意规定.每个分区可定为单类变量区、二类变量区或三类变量区.3.每个交界线上的位移和力的连接条件可以放宽.这个原理为非协调元的厚板有限元法提供理论基础.各种厚板有限元模型可看作这个原理的特殊应用.特别是弹性厚板分区混合变分原理的提出为分区混合有限元法应用于厚板问题打下了基础.     

Stabilized hybrid finite element methods based on the combination of saddle point principles of elasticity problems

How, in a discretized model, to utilize the duality and complementarity of two saddle point variational principles is considered in the paper. A homology family of optimality conditions, different from the conventional saddle point conditions of the domain-decomposed Hellinger-Reissner principle, is derived to enhance stability of hybrid finite element schemes. Based on this, a stabilized hybrid method is presented by associating element-interior displacement with an element-boundary one in a nonconforming manner. In addition, energy compatibility of strain-enriched displacements with respect to stress terms is introduced to circumvent Poisson-locking.

     

Shape-Aware Matching of Implicit Surfaces Based on Thin Shell Energies

A shape sensitive, variational approach for the matching of surfaces considered as thin elastic shells is investigated. The elasticity functional to be minimized takes into account two different types of nonlinear energies: a membrane energy measuring the rate of tangential distortion when deforming the reference shell into the template shell, and a bending energy measuring the bending under the deformation in terms of the change of the shape operators from the undeformed into the deformed configuration. The variational method applies to surfaces described as level sets. It is mathematically well-posed, and an existence proof of an optimal matching deformation is given. The variational model is implemented using a finite element discretization combined with a narrow band approach on an efficient hierarchical grid structure. For the optimization, a regularized nonlinear conjugate gradient scheme and a cascadic multilevel strategy are used. The features of the proposed approach are studied for synthetic test cases and a collection of geometry processing examples.     

Structure-preserving space-time discretization of nonlinear structural dynamics based on a mixed variational formulation

The present work deals with the design of structure-preserving numerical methods in the field of nonlinear elastodynamics and structural dynamics. Structure-preserving schemes such as energy-momentum consistent (EMC) methods are known to exhibit superior numerical stability and robustness. Most of the previously developed schemes are relying on a displacement-based variational formulation of the underlying mechanical model. In contrast to that we present a mixed variational framework for the systematic design of EMC schemes. The newly proposed mixed approach accomodates high-performance mixed finite elements such as the shell element due to Wagner & Gruttmann [1] and the brick element due to Kasper & Taylor [2]. Accordingly, the proposed approach makes possible the structure-preserving extension to the dynamic regime of those high-performance elements. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)