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

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

计算最优控制辛数值方法

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

非线性方程的保辛近似算法

保守体系的微分方程可用Hamilton体系的方法描述,其特点是保辛。两个辛矩阵之和不能保辛,两个辛矩阵的乘积仍是辛矩阵。最常用的小参数摄动法用的是加法,因此对辛矩阵不能保辛。从保辛的角度,要用正则变换。本文针对非线性微分方程,运用自变量坐标变换,对原系统进行变换。由此推导出变换后系统的变分原理。引入Hamilton对偶变量,通过数学变换,得到变系数非线性方程。针对该方程,本文提出了保辛摄动算法。通过数值算例,对不同步长下,保辛摄动法、多尺度摄动法、龙格库塔法和精确解的结果做了比较。数值例题表明,对于非线性方程,本文提出的保辛摄动算法有良好的精度。在步长增大的情况下,保辛摄动保持了良好的稳定性。     

Finite difference method of transinet nonlinear free surface wave problems

A finite difference method is developed for computing the two-dimensional transient potential flow generated by an impulse on the free surface. Both the dynamic and kinematic free surface conditions are considered in nonlinear version. the primary features of the present paper include the use of special coordinates transformations so that the geometry of the flow field is transformed into a time-invariant region, presents an iteration process, by which the velocity potential is computed as the solution of a Poisson equation, the application of fast Fourier transform (FFT) technique results in a tri-diagonal system of equations which can be readily solved by the Thomas algorithm, the computing time is significantly reduced. Thus an efficient technique for handling the transient potential problems is well justified. The feasibility of the present method has been verified by two examples including different initial disturbances respectively.     

一维离散结构能带结构与表面态的辛分析方法

本文基于辛几何方法推导了一维离散周期结构、半无穷周期结构和含杂质半无穷周期结构的本征方程,力求建立一个完整的辛分析体系。通过辛分析,将一维离散半无限周期结构转化到一个元胞上求解,大大简化了计算量。对于含杂质半无穷周期结构,结合辛分析和W-W算法,给出求解含杂质半无穷周期结构本征值问题的精确、稳定和高效算法。数值算例说明了本文算法的有效性。     

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.     

三维波纹型可延展结构振动特性的辛分析

基于三维组装技术的可延展结构具备优异的延展性和可调控性,使其成功应用于各类可延展电子器件的制备中。为了评估该类电子器件的稳定性,本文研究三维波纹型可延展结构的振动行为。首先,基于非线性的Euler-Bernoulli梁理论、Kelvin-Voigt粘弹性理论和考虑压电材料的表面压电效应,建立三维波纹结构的理论分析模型;其次,基于能量原理和扩展拉格朗日运动原理,推导出该结构的动力学控制方程;然后采用二级四阶辛Runge-Kutta求解该动力学方程。通过数值仿真实验验证了辛算法的优越性,同时,还发现随着三维波纹型可延展结构外界激励及其结构参数的变化,该结构的振动特性会从倍周期向分岔和混沌转化;本文结果为三维波纹型可延展结构的优化设计及应用提供理论基础。     

基于辛理论的载流碳纳米管能带分析

基于连续介质力学理论和辛弹性理论,将载流碳纳米管等效为铁木辛柯梁,采用哈密顿变分原理建立了载流碳纳米管的振动控制方程;引入对偶变量将振动控制方程从拉格朗日体系导入到哈密顿体系下;通过波传播方法分析了载流碳纳米管的能带结构;研究了流体密度、流速对载流碳纳米管能带结构的影响;同时计算了载流碳纳米管的散射矩阵. 研究发现:管内流速以及流体密度对剪切频率和弯曲频率有着非常重要的影响. 研究结果表明:载流碳纳米管的剪切频率和弯曲频率因流体的加入而减小,并随流速及流体密度的增大而减小;通过对数值结果的分析发现:载流碳纳米管由于管内流体、流速以及流体密度的作用,会使得载流碳纳米管变的更“软”. 其中,哈密顿体系下所得出的载流碳纳米管弯曲频率随管内流体密度的增加而变小,有别于在拉格朗日体系下非局部梁理论所得的结论. 同时,数值结果表明散射矩阵是酉矩阵,辛体系下的入射波功率流与反射波功率流相等,即功率流守恒,体现了辛弹性力学理论的优越性.      

Symplectic analysis for wave propagation in one-dimensional nonlinear periodic structures

The wave propagation problem in the nonlinear periodic mass-spring structure chain is analyzed using the symplectic mathematical method. The energy method is used to construct the dynamic equation, and the nonlinear dynamic equation is linearized using the small parameter perturbation method. Eigen-solutions of the symplectic matrix are used to analyze the wave propagation problem in nonlinear periodic lattices. Nonlinearity in the mass-spring chain, arising from the nonlinear spring stiffness effect, has profound effects on the overall transmission of the chain. The wave propagation characteristics are altered due to nonlinearity, and related to the incident wave intensity, which is a genuine nonlinear effect not present in the corresponding linear model. Numerical results show how the increase of nonlinearity or incident wave amplitude leads to closing of transmitting gaps. Comparison with the normal recursive approach shows effectiveness and superiority of the symplectic method for the wave propagation problem in nonlinear periodic structures.     

Bifurcation of coisotropic invariant tori under locally Hamiltonian perturbations of integrable systems and nondegenerate deformation of symplectic structure

We study the bifurcation problem for a Cantor set of coisotropic invariant tori in the case where a Liouville-integrable Hamiltonian system undergoes locally Hamiltonian perturbations and, simultaneously, a deformation of the symplectic structure of the phase space. We consider a new case where the deformed symplectic structure generates a nondegenerate matrix of the Poisson brackets of action variables. __________ Translated from Neliniini Kolyvannya, Vol. 9, No. 2, pp. 221–232, April–June, 2006.     

基于对偶变量变分原理的完整约束系统保辛算法

基于对偶变量变分原理,选择积分区间两端位移为独立变量,构造了求解完整约束哈密顿动力系统的高阶保辛算法。首先,利用拉格朗日多项式对作用量中的位移、动量及拉格朗日乘子进行近似;然后,对作用量中不包含约束的积分项采用Gauss积分近似,对作用量中包含约束的积分项采用Lobatto积分近似,从而得到近似作用量;最后,在此近似作用量的基础上,利用对偶变量变分原理,将求解完整约束哈密顿动力系统问题转化为一组非线性方程组的求解。算法具有保辛性和高阶收敛性,能够在位移的插值点处高精度地满足完整约束。算法的收敛阶数及数值性质通过数值算例验证。     

微分求积方法对于桩基大变形分析的应用

本文基于大变形的理论,采用弧坐标首先建立了具有初始位移的桩基的非线性数学模型,一组强非线性的微分-积分方程,其中,地基的抗力采用了Winkeler模型;其次,引入变数变换将微分-积分方程转化为一组非线性微分方程,并用微分求积方法离散了方程组,得到一组离散化的非线性代数方程;最后用Newton-Raphson迭代方法对离散化方程进行了求解,得到了桩基变形前后的构形、弯矩和剪力.计算中选取了两种不同类型的初始位移,并考察了它们对桩基大变形力学行为的影响.     

车辆-轨道系统垂向随机振动的辛方法分析

将轨道视为无限长的周期结构,建立车辆轨道垂向耦合模型. 使用虚拟激励法将随机的轨道不平顺激励转化为确定性的简谐激励,再用辛数学方法求解轨道结构的频率响应特性和耦合系统的响应功率谱. 整个计算模型只有26个自由度,求解过程快速而精确.数值算例中,将该方法与常规有限元方法进行了比较,验证了方法的高效性和正确性,讨论了车辆速度对系统随机响应的影响.      

Multiresolution Symplectic Scheme for wave propagation in complex media

A fast adaptive symplectic algorithm named Multiresolution Symplectic Scheme (MSS) was first presented to solve the problem of the wave propagation (WP) in complex media, using the symplectic scheme and Daubechies‘ compactly supported orthogonal wavelet transform to respectively discretise the time and space dimension of wave equation. The problem was solved in multiresolution symplectic geometry space under the conservative Hamiltonian system rather than the traditional Lagrange system. Due to the fascinating properties of the wavelets and symplectic scheme, MSS is a promising method because of little computational burden, robustness and reality of long-time simulation.     

An analytical symplectic approach to the vibration analysis of orthotropic graphene sheets

A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach.A Hamiltonian system is established by introduc-ing a total unknown vector consisting of the displacement amplitude,rotation angle,shear force,and bending moment. The high-order governing differential equation of the vibra-tion of SLGSs is transformed into a set of ordinary differential equations in symplectic space.Exact solutions for free vibra-tion are obtianed by the method of separation of variables without any trial shape functions and can be expanded in series of symplectic eigenfunctions. Analytical frequency equations are derived for all six possible boundary con-ditions. Vibration modes are expressed in terms of the symplectic eigenfunctions.In the numerical examples,com-parison is presented to verify the accuracy of the proposed method. Comprehensive numerical examples for graphene sheets with Levy-type boundary conditions are given.A para-metric study of the natural frequency is also included.     

High order symplectic conservative perturbation method for time-varying Hamiltonian system

This paper presents a high order symplectic conservative perturbation method for linear time-varying Hamiltonian system.Firstly,the dynamic equation of Hamiltonian system is gradually changed into a high order perturbation equation,which is solved approximately by resolving the Hamiltonian coefficient matrix into a "major component" and a "high order small quantity" and using perturbation transformation technique,then the solution to the original equation of Hamiltonian system is determined through a series of inverse transform.Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes,the transfer matrix is a symplectic matrix;furthermore,the exponential matrices can be calculated accurately by the precise time integration method,so the method presented in this paper has fine accuracy,efficiency and stability.The examples show that the proposed method can also give good results even though a large time step is selected,and with the increase of the perturbation order,the perturbation solutions tend to exact solutions rapidly.     

非线性Schrdinger方程的保辛数值求解

首先将非线性Schrdinger方程化为Hamilton正则方程形式,而后建立Hamilton体系下的变分原理。再用有限元法离散空间坐标,同时对时间坐标进行精细积分,最后运用混合能变分原理,提出非线性Schrdinger方程保辛数值解法。这种解法在保辛的同时,可以让能量和质量在积分格点上亦全部达到守恒。数值算例验证了该方法的有效性。     

温度场下岛-桥结构屈曲薄膜的动力学分析

基于岛-桥结构的柔性电子器件已被用于健康监测和皮肤电子等领域。但是柔性电子器件在工作中极易受工作温度变化等激励产生振动,进而影响器件的灵敏度与可靠性。因此本文研究在温度场作用下岛-桥结构屈曲薄膜的动力学问题。首先,基于Euler-Bernoulli梁理论,建立温度场作用下岛-桥结构屈曲薄膜的动力学控制方程。其次,通过引入新变量,将原动力学方程引入Hamilton体系中,得到相应的Hamilton正则方程。随后,采用辛Runge-Kutta方法求解该Hamilton正则方程,并与经典Runge-Kutta方法对比,数值结果显示了辛算法在求解非线性动力学方程时高精度、高数值稳定性的优势,进一步讨论了温度变化量、预应变、阻尼系数等对屈曲薄膜动力学响应的影响。本文研究为柔性电子器件动力学设计提供了理论参考。     

Structure-preserving properties of three differential schemes for oscillator system

A numerical method for the Hamiltonian system is required to preserve some structure-preserving properties. The current structure-preserving method satisfies the requirements that a symplectic method can preserve the symplectic structure of a finite dimension Hamiltonian system, and a multi-symplectic method can preserve the multi-symplectic structure of an infinite dimension Hamiltonian system. In this paper,the structure-preserving properties of three differential schemes for an oscillator system are investigated in detail. Both the theoretical results and the numerical results show that the results obtained by the standard forward Euler scheme lost all the three geometric properties of the oscillator system, i.e., periodicity, boundedness, and total energy,the symplectic scheme can preserve the first two geometric properties of the oscillator system, and the St¨ormer-Verlet scheme can preserve the three geometric properties of the oscillator system well. In addition, the relative errors for the Hamiltonian function of the symplectic scheme increase with the increase in the step length, suggesting that the symplectic scheme possesses good structure-preserving properties only if the step length is small enough.     

Saint-Venant problem of three-dimensional linear viscoelasticity in the Hamiltonian system

The traditional Saint-Venant problem of three-dimensional viscoelasticity is discussed under the Hamiltonia system with the use of the Laplace integral transformation, and the original problem is transformed into finding eigenvalues and eigenvectors of the Hamiltonia operator matrix. Since local effect near the boundary is usually neglected, all solutions of Saint-Venant problems can be obtained directly by the combinations of zero eigenvectors. Moreover, the adjoint relationships of the symplectic orthogonality of zero eigenvectors in the Laplace domain are generalized to the time domain. Therefore the problem can be discussed directly in the eigenvector space of the time domain, and the iterative application of Laplace transformation is not needed. Simply by applying the adjoint relationships of the symplectic orthogonality, an effective method for boundary condition is given. Based on this method, some typical examples are discussed, in which the whole character of total creep and relaxation of viscoelasticity is clearly revealed.