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

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

NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS

A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.     

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

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

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

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

CrMnFeCoNi纳米晶温度相关的拉伸行为研究

摘要：高熵合金是一种由多种主元元素组成的新型合金。实验研究表明等原子比CrMnFeCoNi高熵合金在低温下具有比室温更高的拉伸强度和断裂韧性。本文针对这一现象,利用分子动力学模拟对平均晶粒尺寸为6 nm的CrMnFeCoNi纳米晶在300、200和77 K下分别进行拉伸模拟。模拟研究揭示了纳米尺度CrMnFeCoNi高熵合金力学行为的温度效应和强韧机理。微结构演化分析表明：低温下,塑性变形阶段,滑移系开动的较少,位错滑移所受的阻力越大,屈服强度和抗拉强度越大;模型破坏时,孔洞缺陷形核较慢,更多孔洞缺陷演化成断口,更多的断口分摊拉伸应变,使得高熵合金纳米晶的低温韧性更好。     

计算最优控制辛数值方法

针对最优控制问题（OCP）的辛数值方法研究及应用进行综述。主要涉及内容包括,动力学系统为常微分方程描述的一般无约束、含不等式约束和状态时滞的最优控制问题,微分代数方程描述的一般无约束、含不等式约束和含切换系统的最优控制问题,以及闭环最优控制问题。从间接法和直接法两个求解框架出发,重点介绍本课题组在保辛算法方面的研究工作。在间接法框架下,首先基于生成函数和变分原理,将OCP保辛离散为非线性方程组,再数值求解方程组。在直接法框架下,将OCP保辛离散为有限维的非线性规划问题（NLP）,再数值求解。针对闭环最优控制问题,提出了保辛模型预测控制、滚动时域估计和瞬时最优控制算法。研究表明,保辛算法具有高精度和高效率的特点,在航空航天和机器人等领域有着广泛应用前景和价值。     

Precise integration method for LQG optimal measurement feedback control problem

ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｓｉｄｅｒｔｈｅｌｉｎｅａｒｓｙｓｔｅｍｏｆｔｈｅｍｅａｓｕｒｅｍｅｎｔｆｅｅｄｂａｃｋｃｏｎｔｒｏｌ ｘ=Ａｘ Ｂｗ Ｂ2 ｕ , ( 1 )ｙ =Ｃｘ ｖ ,( 2 )ｗｈｅｒｅｘｉｓｔｈｅｎ＿ｄｉｍｅｎｓｉｏｎａｌｓｔａｔｅｖｅｃｔｏｒ,ｙｉｓａｑ＿ｖｅｃｔｏｒｏｆｍｅａｓｕｒｅｍｅｎｔｓ,ｕｉｓａｎｍ＿ｖｅｃｔｏｒｏｆｃｏｎｔｒｏｌｉｎｐｕｔｓ,ｗａｎｄｖａｒｅｌ＿ｖｅｃｔｏｒ,ｑ＿ｖｅｃｔｏｒｏｆｗｈｉｔｅ＿ｎｏｉｓｅｐｒｏｃｅｓｓｗｉｔｈｋｎｏｗｎｓｔａｔｉｓｔｉｃａｌｐｒｏｐｅ…     

A stabilized complementarity formulation for nonlinear analysis of 3D bimodular materials

Bi-modulus materials with different mechanical responses in tension and compression are often found in civil, composite, and biological engineering. Numerical analysis of bimodular materials is strongly nonlinear and convergence is usually a problem for traditional iterative schemes. This paper aims to develop a stabilized computational method for nonlinear analysis of 3D bimodular materials. Based on the parametric variational principle, a unified constitutive equa-tion of 3D bimodular materials is proposed, which allows the eight principal stress states to be indicated by three para-metric variables introduced in the principal stress directions. The original problem is transformed into a standard linear complementarity problem (LCP) by the parametric virtual work principle and a quadratic programming algorithm is developed by solving the LCP with the classic Lemke’s algo-rithm. Update of elasticity and stiffness matrices is avoided and, thus, the proposed algorithm shows an excellent conver-gence behavior compared with traditional iterative schemes. Numerical examples show that the proposed method is valid and can accurately analyze mechanical responses of 3D bimodular materials. Also, stability of the algorithm is greatly improved.     

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

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

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

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

Numerical method for dynamics of multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints

Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here,the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem(NCP). An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a constraint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.     

A combined parametric quadratic programming and precise integration method based dynamic analysis of elastic-plastic hardening/softening problems

The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (19872016, 50178016, 19832010) and the Foundation for University Key Teacher by the Ministry of Education of China     

Numerical solutions of linear quadratic control for time-varying systems via symplectic conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented. It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation. The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

Numerical solutions of linear quadratic control for time-varying systems via symplectic par conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented.It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation.The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

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

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

Discrete optimal control for Birkhoffian systems

In this paper, we investigate a discrete variational optimal control for mechanical systems that admit a Birkhoffian representation. Instead of discretizing the original equations of motion, our research is based on a direct discretization of the Pfaff–Birkhoff–d’Alembert principle. The resulting discrete forced Birkhoffian equations then serve as constraints for the minimization of the objective functional. In this way, the optimal control problem is transformed into a finite-dimensional optimization problem, which can be solved by standard methods. This approach yields discrete dynamics, which is more faithful to the continuous equations of motion and consequently yields more accurate solutions to the optimal control problem which is to be approximated. We illustrate the method numerically by optimizing the control for the damped oscillator.     

Solving frictional contact problems by two aggregate-function-based algorithms

Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-function-based algorithms for solving complementarity problems are proposed. One is called the self-adjusting interior point algorithm, the other is called the aggregate function smoothing algorithm. Numerical experiment shows the efficiency of the proposed two algorithms. The project supported by the National Natural Science Foundation of China(10225212, 50178016, 10302007), the National Key Basic Research Special Foundation and the Ministry of Education of China The English text was polished by Ron Marshall.     

The optimal control variational principle and finite elements analysis for viscoplastic dynamics

This paper presents the optimal control variational principle for Perzyna model which is one of the main constitutive relation of viscoplasticity in dynamics. And it could also be transformed to solve the parametric quadratic programming problem. The FEM form of this problem and its implementation have also been discussed in the paper.     

The Method of Integrodifferential Relations for Linear Elasticity Problems

Some possible modifications of the governing equations of the linear theory of elasticity are considered. The stress–strain relation is specified by an integral equality instead of the local Hooke’s law. The modified integrodifferential boundary value problem is reduced to the minimization of a nonnegative functional under differential constraints. A numerical algorithm based on polynomial approximations of unknown functions (stresses and displacements) is developed and applied to linear elasticity problems. The bilateral estimation criteria of solution errors are proposed in order to analyze the algorithm convergence rate. The numerical results obtained by applying the integrodifferential relation method and the conventional variational method are compared and discussed.