基于Lie群的刚体动力学建模及数值计算方法研究

基于Lie群和Lie代数之间的指数映射等价关系,推导了基于Lie群的自由刚体连续动力学方程.结合离散变分原理,推导了其Lie群离散变分积分子.通过证明可知连续和离散动力学系统都具有动量守恒性.对连续动力学方程进行同维化处理,使其变为常规非线性方程组的形式,利用Runge-Kutta法进行求解;基于Runge-Kutta基本理论,推导了直接用于Lie群的Runge-Kutta法,从而使Runge-Kutta法可用于求解变维非线性方程组;通过Lie代数变换,利用Kelly变换和Newton迭代对Lie群离散变分积分子进行求解.仿真对比结果表明,3种算法下的计算结果高度吻合,且能高精度地保持系统的结构守恒和动量守恒性.     

求解周期性分段线性系统动态响应的高效数值方法

基于参变量变分原理,提出了一种求解具有大量间隙弹簧的周期性分段线性系统动态响应的高效率数值方法.通过参变量变分原理来描述间隙弹簧,将复杂的非线性动力问题转化为线性互补问题求解,避免了求解过程中的迭代和刚度阵更新,该算法能准确判断间隙弹簧的压缩和松弛状态.基于结构的周期性和能量传播速度的有限性,提出了一种求解系统动态响应的高效率精细积分方法.该算法指出周期结构的矩阵指数中存在大量的相同元素和零元素,从而不需要重复计算和存储这部分元素,节省了计算量并降低了计算机存储要求.分析了一个五自由度分段线性系统在简谐荷载作用下的动力学行为,包括稳定的周期运动、准周期运动和混沌运动.通过与Runge-Kutta方法的比较,该文方法的正确性和高效率得到了验证.     

一类中立型延迟积分微分方程Runge-Kutta方法的稳定性

讨论了一类非线性中立型延迟积分微分方程Runge-Kutta方法的稳定性.在适当的条件下证明了运用Runge-Kutta方法求解这类方程既是数值稳定的也是渐近稳定的.     

偏微分方程的区间小波自适应精细积分法

利用插值小波理论构造了拟Shannon区间小波,并结合外推法给出了一种求解非线性常微分方程组的时间步长自适应精细积分法,在此基础上构造了求解非线性偏微分方程的区间小波自适应精细积分法（AIWPIM）．数值结果表明,该方法在计算精度上优于将小波和四阶Runge-Kutta法组合得到的偏微分方程的数值求解方法,而计算量则相差不大．该文方法通过Burgers方程给出,但适用于一般情形．     

非线性中立型泛函微分方程Runge-Kutta 法的稳定性和收敛性

本文涉及Runge-Kutta 法变步长求解非线性中立型泛函微分方程(NFDEs) 的稳定性和收敛性.为此, 基于Volterra 泛函微分方程Runge-Kutta 方法的B- 理论, 引入了中立型泛函微分方程Runge-Kutta 方法的EB (expanded B-theory)-稳定性和EB-收敛性概念. 之后获得了Runge-Kutta 方法变步长求解此类方程的EB - 稳定性和EB- 收敛性. 这些结果对中立型延迟微分方程和中立型延迟积分微分方程也是新的.     

常微分方程的半隐Runge-Kutta法

针对某些非线性常微分方程,提出一种算子分裂半隐Runge-Kutta方法,对于非线性部分采用显式计算,对于刚性强的线性部分采用隐式处理.给出了格式的推导,分析了绝对稳定性,并证明了半隐二阶格式的收敛性.相比于显式Runge-Kutta法,半隐格式计算量相近,但改进了稳定性,数值结果显示了方法的合理性和有效性.最后,将算子分裂半隐Runge-Kutta方法应用于数值求解Zakharov偏微分方程组.     

求解非线性Black-Scholes模型的自适应小波精细积分法

针对非线性Black-Scholes方程,基于quasi-Shannon小波函数给出了一种求解非线性偏微分方程的自适应多尺度小波精细积分法.该方法首先利用插值小波理论构造了用于逼近连续函数的多尺度小波插值算子,利用该算子可以将非线性Black-Scholes方程自适应离散为非线性常微分方程组;然后将用于求解常微分方程组的精细积分法和小波变换的动态过程相结合,并利用非线性处理技术(如同伦分析技术)可有效求解非线性Black-Scholes方程.数值结果表明了该方法在数值精度和计算效率方面的优越性.     

辛Runge-Kutta方法在卫星交会对接中的非线性动力学应用研究

卫星交会对接问题是实现太空平台等空间系统的关键问题之一.考虑了由于地球引力作用而引起的卫星交会对接中的非线性动力学问题.首先,采用能量方法给出Lagrange函数;然后,通过引入广义坐标和广义动量,以及Legendre变换,得到Hamilton方程;随后,采用辛Runge-Kutta方法求解该Hamilton方程,并与传统的四阶Runge-Kutta方法对比.数值结果表明:辛Runge-Kutta方法能够在积分过程中长时间保持系统的固有特性,为天体动力学问题的研究提供了良好的数值方法.     

刚性多滞量积分微分方程的Runge-Kutta 方法

本文研究了求解刚性多滞量积分微分方程的Runge-Kutta方法的非线性稳定性和计算有效性.经典Runge—Kutta方法连同复合求积公式和Pouzet求积公式被改造用于求解一类刚性多滞量Volterra型积分微分方程.其分析导出了:在适当条件下,扩展的Runge-Kutta方法是渐近稳定和整体稳定的.此外,数值试验表明所给出的方法是高度有效的.     

基于非结构自适应网格的二维Euler方程数值求解方法研究

提出了一种基于非结构自适应网格的二维Euler方程的数值解法.采用有限体积法进行空间离散,通量计算采用Jamson中心格式,使得它适用于任意多边形计算单元.为了得到定常解,采用一种显式的四步Runge-Kutta迭代方法对时间进行积分.根据流场参数的变化梯度确定加密边,由加密准则进行自适应网格剖分,然后得到分布合理的加密过后的网格.求解二维Euler方程,对NACA0012翼型进行了数值模拟,通过对自适应前后的数值解的对比,说明所建立的方法是正确的.     

小垂度粘弹性索非线性响应及振动主动控制

研究具有初始应力的小垂度粘弹性索的非线性动态响应及振动主动控制。在假定索材料的本构关系为一般微分本构类型的基础上,建立小垂度粘弹性索的运动微分方程;应用Galerkin方法将其转化为可用Runge-Kutta数值积分方法求解的一系列三阶非线性常微分方程。在仅考虑面内的横向振动及忽略非线性的情况下得到了连续状态空间中的状态方程,将状态方程离散为差分方程形式,并用矩阵指数来逐步近似状态转移矩阵;基于二次性能指标的最小化得到了最优的控制力与状态向量。最后通过数值仿真研究说明了粘性参数对索动态响应的影响。     

Diffraction of nonstationary waves at curvilinear cavities in orthotropic bedrock

We propose a method of solving the two-dimensional dynamic problem for an elastic orthotropic body. On the basis of the fundamental solutions constructed for the equations of motion in displacements in the Laplace transform space with respect to time and the boundary element method, solving the boundary problem is reduced to solving a finite system of linear algebraic equations. We carry out an investigation of the stressed state of the bedrock near a tunnel cavity of circular cross-section. Two figures. Bibliography: littes. Translated fromTeoreticheskaya i Prikladnaya Mekhanika, No. 22, pp. 96–102 1991.     

Approximation, integration and differentiation of time functions using a set of orthogonal hybrid functions (HF) and their application to solution of first order differential equations

Differential equations of different types and orders are of utmost importance for mathematical modeling of control system problems. State variable method uses the concept of expressing n number of first order differential equations in vector matrix form to model and analyze/synthesize control systems.The present work proposes a new set of orthogonal hybrid functions (HF) which evolved from synthesis of sample-and-hold functions (SHF) and triangular functions (TF). This HF set is used to approximate a time function in a piecewise linear manner with the mean integral square error (MISE) much less than block pulse function based approximation which always provides staircase solutions.The operational matrices for integration and differentiation in HF domain are also derived and employed for solving non-homogeneous and homogeneous differential equations of the first order as well as state equations. The results are compared with exact solutions, the 4th order Runge-Kutta method and its further improved versions proposed by Simos [6]. The presented HF domain theory is well supported by a few illustrations.     

Time integration of inelastic material models exhibits an order reduction for higher-order methods – can this be avoided?

For the numerical treatment of inelastic material behavior within the finite element method a partitioned ansatz is standard in most of the software frameworks; the weak form of equilibrium is discretized in space and solved on a global level, whereas the initial value problem for the evolution equations of internal state variables is separately solved on a local, i.e. Gauss-point level, where strains, derived from global displacements, serve as input, [1]. Applying higher order methods (p > 2) to the time integration of plasticity models an order reduction is reported where Runge-Kutta schemes have shown hardly more than order two at best [2, 3]. In the present contribution, we analyze the reason for order reduction and in doing so, introduce an improved strain approximation and switching point detection which play a crucial role for the convergence order of multi-stage methods used in this context. We apply Runge-Kutta methods of Radau IIa class to the evolution equations of viscoelastic and elastoplastic material models and show ther improved performence in numerical examples. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An improved approximate Newton method for implicit Runge-Kutta formulas

Implicit Runge-Kutta (IRK) methods (such as the s-stage Radau IIA method with s=3,5, or 7) for solving stiff ordinary differential equation systems have excellent stability properties and high solution accuracy orders, but their high computing costs in solving their nonlinear stage equations have seriously limited their applications to large scale problems. To reduce such a cost, several approximate Newton algorithms were developed, including a commonly used one called the simplified Newton method. In this paper, a new approximate Jacobian matrix and two new test rules for controlling the updating of approximate Jacobian matrices are proposed, yielding an improved approximate Newton method. Theoretical and numerical analysis show that the improved approximate Newton method can significantly improve the convergence and performance of the simplified Newton method.     

多自由度强非线性颤振分析的增量谐波平衡法

对多个自由度上含有强非线性项系统的颤振问题,推广应用增量谐波平衡法进行分析．考虑带有强非线性立方平移和俯仰刚度项的二元机翼颤振方程,首先将方程用矩阵形式表示,然后把振动过程分解成为振动瞬态的持续增量过程,再采用振幅作为控制参数应用谐波平衡法,以这种推广的增量谐波平衡法求得方程解的表达式,并由此分析系统的分岔现象、极限环颤振现象和谐波项数的取值问题,最后用龙格-库塔数值方法进行验算,结果表明:分析多个自由度的强非线性颤振,增量谐波平衡法是精确有效的．     

Runge-Kutta collocation methods for differential-algebraic equations of indices 2 and 3

Stiffly accurate Runge-Kutta collocation methods with explicit first stage are examined. The parameters of these methods are chosen so as to minimize the errors in the solutions to differential-algebraic equations of indices 2 and 3. This construction results in methods for solving such equations that are superior to the available Runge-Kutta methods.     

On Time Discretizations for Spectral Methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.     

求解二维浅水波方程的旋转混合格式北大核心CSCD

针对二维浅水波方程数值求解问题,构造了一种旋转通量混合格式.空间方向上,该算法利用浅水波方程通量函数的旋转不变性,在单元界面法线方向及单元界面切线方向上采用可消除红斑现象的HLL与满足热力学第二定律的熵稳定加权混合数值通量函数,时间方向上采用三阶强稳定Runge-Kutta法.数值结果表明,该混合格式对于二维浅水波方程数值求解具有分辨率高的良好特性.