非线性最优控制系统的时程精细计算研究

针对非线性最优控制问题 ,通过一阶 Taylor级数展开 ,得到线性化的动力学方程 ,进而在方程原变量的基础上 ,引入对偶向量 (Lagrange乘子向量 ) ,将动力学方程从 Lagrange体系引入到了 Hamilton体系 ,在全状态下 ,从一个新的角度对非线性最优控制问题进行了描述 ,进一步基于时程精细积分理论 ,对其方程进行了有效的精细求解 ,并通过算例说明了文中方法的有效性     

基于高斯原理的非理想系统动力学建模

高斯原理给出了通过求函数极值、从可能运动中鉴别出真实运动的规则, 它可以使得多体系统动力学问题不需通过求解微分(代数)方程, 而是采用求解最小值的优化方法来解决, 从而提供了一种适用于优化算法的建模思路, 因此, 如何定义恰当的高斯拘束函数是动力学优化方法得以实现的前提. 对于理想系统而言, 约束对系统的作用可以通过约束方程来体现, 故高斯拘束可表达为系统质点加速度的函数, 系统的动力学问题因此可以描述为目标函数为高斯拘束函数、优化变量为质点加速度的约束最优化问题; 当系统中需要考虑干摩擦等非理想因素时, 部分相互作用不能被所定义的约束方程所涵盖而需要采用额外的物理规律来描述, 这种相互作用破坏了原有的针对理想系统的高斯拘束函数的极值特性. 基于变分类的高斯原理, 推导并证明了目标函数以理想约束力所表达的非理想系统的极值原理, 针对目前文献中用于非理想系统的高斯原理进行了讨论, 指出其实际为文中的极值原理在非理想约束力与理想约束力无明显关联时的一种特殊表达形式, 当非理想约束力与理想约束力有明显的函数关系(如库仑摩擦定律中滑动摩擦力与法向约束力间的线性关系)时, 该形式失效; 同时根据文中的极值原理, 得到了考虑库仑摩擦时非理想的多体系统动力学问题的优化模型. 例子中分析了优化模型及相应的线性互补性模型的关系, 分析发现在满足刚体滑动问题的唯一性条件下二者互为充分必要条件, 从而证明了文中优化模型的可靠性; 并采用优化计算方法进行了动力学模拟, 模拟结果显示了将高斯原理与优化算法相结合的可行性及有效性.      

DYNAMIC MODELING OF NONIDEAL SYSTEM BASED ON GAUSS'S PRINCIPLE$^{\bf 1)}$

高斯原理给出了通过求函数极值、从可能运动中鉴别出真实运动的规则, 它可以使得多体系统动力学问题不需通过求解微分(代数)方程, 而是采用求解最小值的优化方法来解决, 从而提供了一种适用于优化算法的建模思路, 因此, 如何定义恰当的高斯拘束函数是动力学优化方法得以实现的前提. 对于理想系统而言, 约束对系统的作用可以通过约束方程来体现, 故高斯拘束可表达为系统质点加速度的函数, 系统的动力学问题因此可以描述为目标函数为高斯拘束函数、优化变量为质点加速度的约束最优化问题; 当系统中需要考虑干摩擦等非理想因素时, 部分相互作用不能被所定义的约束方程所涵盖而需要采用额外的物理规律来描述, 这种相互作用破坏了原有的针对理想系统的高斯拘束函数的极值特性. 基于变分类的高斯原理, 推导并证明了目标函数以理想约束力所表达的非理想系统的极值原理, 针对目前文献中用于非理想系统的高斯原理进行了讨论, 指出其实际为文中的极值原理在非理想约束力与理想约束力无明显关联时的一种特殊表达形式, 当非理想约束力与理想约束力有明显的函数关系(如库仑摩擦定律中滑动摩擦力与法向约束力间的线性关系)时, 该形式失效; 同时根据文中的极值原理, 得到了考虑库仑摩擦时非理想的多体系统动力学问题的优化模型. 例子中分析了优化模型及相应的线性互补性模型的关系, 分析发现在满足刚体滑动问题的唯一性条件下二者互为充分必要条件, 从而证明了文中优化模型的可靠性; 并采用优化计算方法进行了动力学模拟, 模拟结果显示了将高斯原理与优化算法相结合的可行性及有效性.     

Optimal control of attitude for coupled-rigid-body spacecraft via Chebyshev-Gauss pseudospectral method

The attitude optimal control problem(OCP) of a two-rigid-body spacecraft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a discrete nonlinear programming(NLP) problem using the Chebyshev-Gauss pseudospectral method(CGPM). Solutions of the NLP problem can be obtained using the sequential quadratic programming(SQP) algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss(CG) points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform(FFT). To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and control variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.     

多体系统接触碰撞问题的牛顿-欧拉线性互补方法

基于非光滑动力学方法的多体系统接触碰撞分析是目前多体系统动力学的研究热点.本文采用牛顿-欧拉方法建立多体系统接触、碰撞问题的动力学模型,给出一种牛顿-欧拉型线性互补公式.该建模方法与目前一般采用的拉格朗日建模方法的不同之处是约束条件中除了库仑摩擦、单边约束之外还含有光滑等式约束.在建立系统动力学模型时,首先解除摩擦约束和单边约束得到原系统对应的基本系统.牛顿-欧拉方法采用最大数目坐标建立基本系统的动力学方程,由于坐标不相互独立,因此基本系统中带有等式约束,其数学模型为一组微分代数方程.借助约束雅可比矩阵,在基本系统微分代数方程中添加摩擦接触和单边约束对应的拉氏乘子,就可以得到系统全局运动的具有变拓扑结构特征的动力学方程,再结合非光滑约束互补条件便可构成完备的系统动力学模型.完备的动力学模型由动力学微分方程以及等式约束和不等式约束组成.线性互补公式采用分块矩阵形式进行推导,简化了推导过程.数值计算采用基于线性互补的时间步进算法.时间步进算法是目前流行的非光滑数值算法,其突出特点是可以免去数值积分中繁琐的事件检测过程,而数值积分过程中通过对线性互补问题的求解可以确定系统的触-离状态.通过对典型的曲柄滑块间隙机构进行数值分析,验证本文方法的有效性.     

The elastodynamic solution for a solid sphere and dynamic stress-focusing phenomenon

This paper presents an analytical method of solving the elastodynamic problem of a solid sphere.The basic solution of the elastodynamic problem is decomposed into a quasi-static solution satisfying the inhomogeneous compound boundary conditions and a dynamic solution satisfying the homogeneous compound boundary conditions.By utilizing the variable transform,the dynamic equation may be transformed into Bassel equation.By defining a finite Hankel transform,we can easily obtain the dynamic solution for the inhomogeneous dynamic equation.Thereby,the exact elastodynamic solution for a solid sphere can be obtained.From results carried out,we have observed that there exists the dynamic stress-focusing phenomenon at the center of a solid sphere under shock load and it results in very high dynamic stress-peak.     

Synthesis of a controller for stabilizing the motion of a rigid body about a fixed point

A method for the approximate design of an optimal controller for stabilizing the motion of a rigid body about a fixed point is considered. It is assumed that rigid body motion is nearly the motion in the classical Lagrange case. The method is based on the common use of the Bellman dynamic programming principle and the averagingmethod. The latter is used to solve theHamilton–Jacobi–Bellman equation approximately, which permits synthesizing the controller. The proposed method for controller design can be used in many problems close to the problem of motion of the Lagrange top (the motion of a rigid body in the atmosphere, the motion of a rigid body fastened to a cable in deployment of the orbital cable system, etc.).     

俯仰运动圆柱贮箱中液体的非线性晃动

首次对储仰运动圆柱贮箱中液体的有限幅值晃动问题进行了解析研究。首先建立了描述俯仰和／或偏航运动贮箱中液体晃动的非线性偏微分方程组,而后提出了相应的变分原理,建立了压力体积分形式的Ｌａｇｒａｎｇｅ函数,通过变分方程,最终得到措述俯仰和／或偏航运动圆柱贮箱中液体晃动的非线性动力学微分方程组,该动力学方程组自然满足液体自由表面的运动学和动力学办界条件。而后动用多尺度法求解了所得的动力学方程组,对非线性液     

多体动力学的几何积分方法研究进展

动力系统的几何积分研究是近20年来工程计算领域非常活跃的方向.多体动力学方程(微分方程, 微分代数方程)是一类典型的动力系统,将其从Lagrange体系向Hamilton系统过渡,目的在于从欧氏几何过渡到辛几何形态, 将对偶变量引入到力学研究中,然后利用辛几何的数学框架对多体系统动力学方程进行数值计算,可以预知多体动力学系统的一些定性信息,并在数值离散时能保持这些定性性质特征,尤其在表示关键的物理意义时需要强调保持这些几何性质.简要介绍多体系统(无约束多刚体系统、完整约束多刚体系统和柔性多体系统)的Hamilton正则方程的建立和几何积分方法的构造,着重介绍了在多体动力学计算中非常有应用前景的高阶辛算法(合成辛算法、分裂合成辛算法和辛精细积分法)、多辛算法,以及广义Hamilton 系统与Lie 群积分方法等计算几何力学方法, 并对Lie群积分的投影方法、流形局部坐标法等方法进行了阐述.      

带有末端集中质量的双连杆柔性机械臂主动控制

对带有末端集中质量的双连杆柔性机械臂的主动控制进行了研究,给出系统的动力学方程,采用非线性解耦反馈控制方法分别得出系统大范围运动方程和柔性臂的动力学方程,采用机械臂逆动力学方法和LQR方法分别设计大范围运动控制律和压电作动器控制律.仿真结果显示,本文控制方法能够有效地进行机械臂的轨迹跟踪,柔性臂的弹性振动可以得到有效抑制.     

用拉格朗日函数分析不同模量静不定桁架内力

在外载荷作用下的不同模量静不定桁架平衡问题,是任意有限多个自变量的多元函数在任意有限多个约束条件下的极值问题,对采用拉格朗日乘数法求解此类极值问题进行了数学证明．通过求解不同模量静不定桁架极限载荷的几个算例,阐述拉格朗日乘数法在计算不同模量静不定桁架极限载荷中的应用．研究结果表明：采用拉格朗日乘数法求解不同模量静不定桁架极限载荷的通用性较强,用拉格朗日乘数法求解不同模量静不定桁架极限载荷的方法不但克服了常规方法需利用几何关系建立协调方程的缺陷,且具有力学概念清晰直观、计算过程简便、便于工程设计人员在实际中掌握和应用．     

Scattering of harmonic antiplane shear waves by an arc-shaped interfacial crack in functionally graded annular bi-material guide rail

Abstract

The dynamic behavior of an arc-shaped interfacial crack in an orthotropic functionally graded annular bi-material structure is investigated. In order for the analysis to be executable, the material properties are assumed to vary with the power function of the radial coordinates. By applying the separation variable method, the boundary value problem of the partial differential equation describing the fracture problem of this article can be transformed into a Cauchy kernel singular integral equation with the unknown jump of displacements across the crack surfaces. The obtained integral equation is solved numerically by Lobatto–Chebyshev collocation method to show the effects of the geometric and physical parameters upon the dynamic stress field near the crack tips.

Communicated by Kuang-Hua Chang.     

非线性动力方程的增维精细积分法

对线性定常结构的动力系统提出的精细积分法，能得到在数值上逼近于精确解的结果。但是对于非齐次动力方程却涉及到矩阵求逆的困难，而且通常与时间有关的非齐次项不能进入精细积分的细化过程。采用增维的方法，将非齐次动力方程化为齐次方程，在实施精细积分的过程中不必进行矩阵求逆。这种处理方法对于程序实现和提高数值计算的稳定性十分有利，而且在大型问题中可明显提高计算效率，数值算例显示本文方法是有效的。     

考虑时滞的主动悬架系统控制策略对比研究全文替换

本文以二自由度四分之一汽车悬架系统为研究对象,采用不同控制策略对考虑时滞的悬架系统控制特性进行研究,并对控制效果进行对比分析．首先,采用第二类拉氏方程建立考虑时滞的二自由度悬架控制系统的动力学模型;然后分别基于状态变换法和H_∞控制理论,设计系统的时滞反馈控制律．其中状态变换法主要通过系统状态变量的转换,将系统时滞控制方程转换成不显含时滞的动力学方程,然后采用传统的二次型最优控制方法对系统进行控制．H_∞控制主要通过使用Lyapunov-Krasovskii泛函和自由权矩阵法,设计考虑时滞的H_∞控制律;最后在Matlab/Simulink平台上对控制系统进行仿真分析,并在相同时滞下对两种控制结果进行对比．研究表明,两种控制策略在考虑时滞的情况下均可保证系统的稳定性,且H_∞控制相比于状态变换法有更好的控制效果．     

Elastodynamic solution of a non-homogeneous orthotropic hollow cylinder

A theoretical method for analyzing the axisymmetric plane strain elastodynamic problem of a non-homogeneous orthotropic hollow cylinder is developed. Firstly, a new dependent variable is introduced to rewrite the governing equation, the boundary conditions and the initial conditions. Secondly, a special function is introduced to transform the inhomogeneous boundary conditions to homogeneous ones. By virtue of the orthogonal expansion technique, the equation with respect to the time variable is derived, of which the solution can be obtained. The displacement solution is finally obtained, which can be degenerated in a rather straightforward way into the solution for a homogeneous orthotropic hollow cylinder and isotropic solid cylinder as well as that for a non-homogeneous isotropic hollow cylinder. Using the present method, integral transform can be avoided and it can be used for hollow cylinders with arbitrary thickness and subjected to arbitrary dynamic loads. Numerical results are presented for a non-homogeneous orthotropic hollow cylinder subjected to dynamic internal pressure. The project supported by the National Natural Science Foundation of China (10172075 and 10002016)     

Dynamic modeling and simulation of deploying process for space solar power satellite receiver

To reveal some dynamic properties of the deploying process for the solar power satellite via an arbitrarily large phased array (SPS-ALPHA) solar receiver, the symplectic Runge-Kuttamethod is used to simulate the simplified model with the consideration of the Rayleigh damping effect. The system containing the Rayleigh damping can be separated and transformed into the equivalent nondamping system formally to insure the application condition of the symplectic Runge-Kutta method©First, the Lagrange equation with the Rayleigh damping governing the motion of the system is derived via the variational principle. Then, with some reasonable assumptions on the relations among the damping, mass, and stiffness matrices, the Rayleigh damping system is equivalently converted into the nondamping system formally, so that the symplectic Runge-Kutta method can be used to simulate the deploying process for the solar receiver. Finally, some numerical results of the symplectic Runge-Kutta method for the dynamic properties of the solar receiver are reported. The numerical results show that the proposed simplified model is valid for the deploying process for the SPS-ALPHA solar receiver, and the symplectic Runge-Kutta method can preserve the displacement constraints of the system well with excellent long-time numerical stability.     

Self-synchronization theory of dual motor driven vibration system with two-stage vibration isolation frame

This paper studies self-synchronization and stability of a dual-motor driven vibration system with a two-stage vibration isolation frame. Oscillation amplitude of the material box large enough can be ensured on the vibration system in order to screen materials. Reduction of the dynamic load transmitted to the foundation can also be achieved for the vibration system. A Lagrange equation is used to set up the motion differential equations of the system, and a dimensionless coupled equation of the eccentric rotors is obtained using a method of modified average small parameter. According to the existence condition of zero solution in the dimensionless coupled equation of the eccentric rotors, the precondition for commencing self-synchronization motion is achieved.The stability condition of self-synchronization is obtained based on the Routh-Hurwitz criterion. The theoretical analysis is validated by simulations and experiments.     

求解非线性结构动力方程的预估校正-辛时间子域法

提出了求解非线性结构动力方程的预估校正-辛时间子域法。首先,将结构非线性动力方程转换为状态空间方程,在任一时间子域内利用改进的欧拉法对各离散时刻的状态变量值进行预估和校正。然后,将离散的非线性项用Lagrange插值多项式展开并视为外荷载,结合辛时间子域法即可求解非线性动力系统的响应。这种方法不必对状态矩阵求逆,无需计算高阶导数,计算简单,格式统一,易于编程。算例结果表明,本文方法具有较高的计算精度、效率和稳定性,是一种求解非线性结构动力方程的有效方法。     

各向异性平板开孔动应力集中问题的研究

采用各向异性平板弯曲波动理论及摄动方法，对正交各向异性平板开孔弯曲波的散射及动应力集中问题进行了分析研究，得到了此种平板稳态弯曲波动问题的渐近形式的分析解。同时采用保角映射技术，为求解正交各向异性平板开孔弹性波的散射及动应力集中问题提供了一种统一规范的方法。     

Symplectic multi-level method for solving nonlinear optimal control problem

By converting an optimal control problem for nonlinear systems to a Hamiltonian system,a symplecitc-preserving method is proposed.The state and costate variables are approximated by the Lagrange polynomial.The state variables at two ends of the time interval are taken as independent variables.Based on the dual variable principle,nonlinear optimal control problems are replaced with nonlinear equations.Furthermore,in the implementation of the symplectic algorithm,based on the 2N algorithm,a multilevel method is proposed.When the time grid is refined from low level to high level,the initial state and costate variables of the nonlinear equations can be obtained from the Lagrange interpolation at the low level grid to improve efficiency.Numerical simulations show the precision and the efficiency of the proposed algorithm in this paper.