参变量变分原理的提出、发展与应用

参变量变分原理及其参数二次规划算法是由钟万勰院士1985年针对弹性接触边界非线性问题首次提出来的,经过将近40年的不断发展,目前参变量变分原理已经成功应用于各个领域,其中包括弹塑性分析、接触问题、润滑力学、岩土力学、变刚度杆系结构、先进材料性能分析、材料的蠕变与损伤、柔性结构力学和LQ最优控制等各个工程领域。本文首先回顾了参变量变分原理的起源,介绍了参变量变分原理的基本概念,然后以弹塑性分析问题为例,阐明建立参变量变分原理的理论模型以及实现数值参数二次规划求解原理,最后详细回顾了参变量变分原理的基本理论与相应数值算法在各个领域的发展及其工程应用,展示了参变量变分原理在求解各类非线性问题的特色与优势。     

材料弹塑性性能数值模拟的多尺度方法

将均匀化方法和渐近分析(Asymptotic Analysis)与参变量变分原理相结合提出了一种模拟复合材料非线性性能的多尺度数值方法.该方法用渐近分析建立宏-细观变量之间的联系,利用参变量变分原理计算非线性响应,求解过程采用迭代算法.为提高计算精度,针对Von-Mises准则和Tsai-Hill准则,提出了一个基于参变量变分原理的改进算法,算例表明该方法可以显著消除传统方法采用线性展开式构造线性互补条件所带来的误差.     

基于参变量变分原理的非均质多边形微结构材料弹塑性分析

为了更好地模拟复合材料及含夹杂非均质材料等的宏观弹塑性力学性能,简化有限元建模时间和减少有限元模拟计算量。本文基于参变量变分原理,提出了一种采用任意多边形弹塑性单元进行结构非线性分析的参数二次规划算法,给出了参变量最小势能原理以及最终的二次规划模型,并在有限元分析与优化设计软件系统JIFEX上进行了程序实现。数值算例证明了本文方法的正确与可行性。     

解Drucker-Prager塑性问题的二阶锥互补法

基于经典弹塑性理论中多数屈服准则具有凸锥数学结构的事实,将在大规模计算中更具潜力的锥规划法引入弹塑性分析。考虑到弹塑性流动理论有关联与非关联之分,本文提出利用锥型互补法求解弹塑性问题。具体以Drucker-Prager弹塑性模型为例,首先利用最大塑性功耗散原理和圆锥对偶理论等工具,建立了弹塑性本构方程的等价二阶锥互补模型;然后,基于参变量变分原理和有限元技术,建立了弹塑性增量分析的二阶锥线性互补模型;最后,利用一类半光滑Newton算法求解。数值算例表明了本文方法的有效性。     

正交各向异性弹塑性摩擦接触问题的数值求解

采用正交各向异性摩擦定律对三维弹塑性摩擦接触问题进行分析,基于参变量变分原理,经过有限元离散,将问题化为线性互补问题模型,之后给出一个求解互补问题的非内点光滑化算法．对三维接触问题,滑动方向的确定一直是个难点,为此,该文采用作者提出的组合规划法和迭代法对各向异性摩擦本构模型进行分析,数值结果说明了模型与算法的正确性。     

索网天线的参变量变分及非线性有限元方法

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.      

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.     

基于参变量变分原理的直升机系留载荷高性能计算方法

为了保证直升机在舰船上的安全性,必须使用系留设备将直升机系留在舰船上。直升机的系留问题可简化为由机身刚体、索具和起落架组成的杆件系统,索具只承受拉力而不承受压力,起落架只承受压力而不承受拉力。因此,直升机系留问题为典型的强非线性问题,需要发展有效的求解算法。在考虑大变形的情况下,基于参变量变分原理建立了求解直升机系留载荷的高性能计算方法。 该方法利用参变量变分原理能够准确判断索具和起落架的拉压状态,并将材料非线性静力问题转换为线性静力互补问题求解,极大地提高了结果的收敛性。数值算例中,通过与有限元通用软件NASTRAN和ABAQUS计算结果比较,证实了该方法的精确性、收敛性及高效性。     

解三维摩擦接触问题的一个二阶锥线性互补法

针对三维摩擦接触问题的求解,给出了一种基于参变量变分原理的二阶锥线性互补法. 首先,基于三维Coulomb摩擦锥在数学表述上属于二阶锥的事实,利用二阶锥规划对偶理论,建立了三维Coulomb摩擦接触条件的参变量二阶锥线性互补模型,它是二维Coulomb摩擦接触条件参变量线性互补模型在三维情形下的自然推广;随后,利用参变量变分原理与有限元方法,建立了求解三维摩擦接触问题的二阶锥线性互补法. 较之于将三维Coulomb摩擦锥进行显式线性化的线性互补法,该方法无需对三维Coulomb摩擦锥进行线性化,因而在保证精度的前提下所解问题的规模要小很多. 最后通过算例展示了该方法的特点.      

输入受限LQ控制的参变量变分原理和算法

在最优控制理论中根据模拟理论思想发展了塑性力学和接触力学中的参变量变分原理, 并建立了控制输入受限的线性二次(linear quadratic, LQ)最优控制问题的求解新方程---耦合的Hamilton正则方程与线性互补方程. 通过将连续时间离散成一系列等间距时间区段, 在离散时域内采用参数二次规划方法给出数值求解输入受限的LQ最优控制问题的新算法. 数值仿真验证了该算法在求解控制输入受限的LQ最优控制问题中的有效性, 并且该算法具有较快的收敛性, 在大步长下具有较高的计算精度.      

Rigorous bounds on the effective moduli of composites and inhomogeneous bodies with negative-stiffness phases

We review the theoretical bounds on the effective properties of linear elastic inhomogeneous solids (including composite materials) in the presence of constituents having non-positive-definite elastic moduli (so-called negative-stiffness phases). Using arguments of Hill and Koiter, we show that for statically stable bodies the classical displacement-based variational principles for Dirichlet and Neumann boundary problems hold but that the dual variational principle for traction boundary problems does not apply. We illustrate our findings by the example of a coated spherical inclusion whose stability conditions are obtained from the variational principles. We further show that the classical Voigt upper bound on the linear elastic moduli in multi-phase inhomogeneous bodies and composites applies and that it imposes a stability condition: overall stability requires that the effective moduli do not surpass the Voigt upper bound. This particularly implies that, while the geometric constraints among constituents in a composite can stabilize negative-stiffness phases, the stabilization is insufficient to allow for extreme overall static elastic moduli (exceeding those of the constituents). Stronger bounds on the effective elastic moduli of isotropic composites can be obtained from the Hashin–Shtrikman variational inequalities, which are also shown to hold in the presence of negative stiffness.     

On a regular convex solver potential for an elastic-damage constitutive model: a theoretical analysis

This work is related with the proposition of a so-called regular or convex solver potential to be used in numerical simulations involving a certain class of constitutive elastic-damage models. All the mathematical aspects involved are based on convex analysis, which is employed aiming a consistent variational formulation of the potential and its conjugate one. It is shown that the constitutive relations for the class of damage models here considered can be derived from the solver potentials by means of sub-differentials sets. The optimality conditions of the resulting minimisation problem represent in particular a linear complementarity problem. Finally, a simple example is present in order to illustrate the possible integration errors that can be generated when finite step analysis is performed.     

Using the Euler–Lagrange variational principle to obtain flow relations for generalized Newtonian fluids

The Euler–Lagrange variational principle is used to obtain analytical and numerical flow relations in cylindrical tubes. The method is based on minimizing the total stress in theflow duct using the fluid constitutive relation between stress and rate of strain. Newtonian and non-Newtonian fluid models, which include power law, Bingham, Herschel–Bulkley, Carreau, and Cross, are used for demonstration.     

Towards variational constitutive updates for non-associative plasticity models at finite strain: Models based on a volumetric-deviatoric split

In this paper, an enhanced variational constitutive update suitable for a class of non-associative plasticity theories at finite strain is proposed. In line with classical numerical formulations for plasticity models, such as the by now established return-mapping algorithm, variational constitutive updates represent a numerical method for computing the unknown state variables. However, in contrast to conventional algorithms, variational constitutive updates are fully variational, i.e., all unknown variables follow jointly from minimizing a certain potential. In addition to the physical and mathematical elegance of these variational schemes, they show several practical advantages as well. For instance, numerically efficient and robust optimization schemes can be directly employed for solving the resulting minimization problem. Since mathematically, plasticity is a non-smooth problem and often, it leads to highly singular systems of equations as known from single crystal plasticity, a robust implementation is of utmost importance. So far, variational constitutive updates have been developed for different classes of standard dissipative solids, i.e., solids characterized by associative evolution equations and flow rules. In the present paper, this framework is extended to a certain class of non-associative plasticity models at finite strain. All models falling into this class show a volumetric-deviatoric split of the Helmholtz energy and the yield function. Typical prototypes are Drucker-Prager or Mohr-Coulomb models playing an important role in soil mechanics. The efficiency and robustness of the resulting algorithmic formulation is demonstrated by means of selected numerical examples.     

Analysis of material instabilities in inelastic solids by incremental energy minimization and relaxation methods: evolving deformation microstructures in finite plasticity

We propose an approach to the definition and analysis of material instabilities in rate-independent standard dissipative solids at finite strains based on finite-step-sized incremental energy minimization principles. The point of departure is a recently developed constitutive minimization principle for standard dissipative materials that optimizes a generalized incremental work function with respect to the internal variables. In an incremental setting at finite time steps this variational problem defines a quasi-hyperelastic stress potential. The existence of this potential allows to be recast a typical incremental boundary-value problem of quasi-static inelasticity into a principle of minimum incremental energy for standard dissipative solids. Mathematical existence theorems for sufficiently regular minimizers then induce a definition of the material stability of the inelastic material response in terms of the sequentially weakly lower semicontinuity of the incremental variational functional. As a consequence, the incremental material stability of standard dissipative solids may be defined in terms of the quasi-convexity or the rank-one convexity of the incremental stress potential. This global definition includes the classical local Hadamard condition but is more general. Furthermore, the variational setting opens up the possibility to analyze the post-critical development of deformation microstructures in non-stable inelastic materials based on energy relaxation methods. We outline minimization principles of quasi- and rank-one convexifications of incremental non-convex stress potentials for standard dissipative solids. The general concepts are applied to the analysis of evolving deformation microstructures in single-slip plasticity. For this canonical model problem, we outline details of the constitutive variational formulation and develop numerical and semi-analytical solution methods for a first-level rank-one convexification. A set of representative numerical investigations analyze the development of deformation microstructures in the form of rank-one laminates in single slip plasticity for homogeneous macro-deformation modes as well as inhomogeneous macroscopic boundary-value problems. The well-posedness of the relaxed variational formulation is indicated by an independence of typical finite element solutions on the mesh-size.     

Generalized variational principles for boundary value problem of electromagnetic field in electrodynamics

An expression of the generalized principle of virtual work for the boundary value problem of the linear and anisotropic electromagnetic field is given. Using Chien＇s method, a pair of generalized variational principles （GVPs） are established, which directly leads to all four Maxwell＇s equations, two intensity-potential equations, two constitutive equations, and eight boundary conditions. A family of constrained variational principles is derived sequentially. As additional verifications, two degenerated forms are obtained, equivalent to two known variational principles. Two modified GVPs are given to provide the hybrid finite element models for the present problem.     

PHASE FIELD SIMULATION OF DOMAIN SWITCHING IN FERROELECTRIC SINGLE CRYSTAL WITH ELECTRICALLY PERMEABLE AND IMPERMEABLE CRACKS

Domain switching around electrically permeable and impermeable cracks in ferroelectric single crystals subjected to a mechanical load is investigated by using a phase field model.It is found that the e...     

不可压饱和粘弹性多孔介质固结问题的Gurtin型变分原理

基于多孔介质理论,在固相骨架和孔隙流体微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,利用卷积积分的性质,本文首先建立了以固相骨架位移、孔隙流体相对速度和孔隙流体压力为宗量的流体饱和粘弹性多孔介质固结问题的一个Gurtin型变分原理.其次,利用Lagrange乘子法解除相关的变分约束条件,建立了流体饱和粘弹性多孔介质固结问题的若干广义Gurtin型变分原理,包括第三类的Hu-Washizu型变分原理.最后,简单讨论了等价初边值问题的相应变分原理.这些Gurtin型变分原理的建立不仅丰富了饱和粘弹性多孔介质的相关理论,而且为相关数值模拟方法,如有限元法、无网格法等的建立奠定了理论基础.     

橡胶类材料大变形本构关系及其有限元方法

讨论大变形拟不可压缩模胶类材料的本构关系及有限元分析方法，采用乘法分解，将变形梯度表示成等容和体积变形两部分，在此基础上，推导了克希荷夫应力和格林应力表示的Ｙｅｏｈ形式应变能模胶类材料的本构关系及数值处理方法，为处理不可压缩问题，采用三场变分原理，其中静水压力，体积膨胀，以及位移均作为独立变量进行处理，并指出该变分原理同胡－鹫津广义变分原理的联系，变形采用相容等参插值，压力及体积膨胀采用低阶插值，