A PRECISE INTEGRATION APPROACH FOR THE DYNAMIC-STIFFNESS MATRIX OF STRIP FOOTINGS ON A LAYERED MEDIUM 1）

提出应用精细积分算法计算多层地基的动力刚度问题.精细积分是计算层状介质中波传播的高效而精确的数值方法.利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组,运用精细积分方法求解格林函数,最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内,进而得到刚性条带基础频率域的动力柔度或刚度矩阵.所建议的精细积分算法,可以避免一般传递矩阵计算中的指数溢出问题,对各种情况有广泛的适应性,计算稳定,在高频段可以保障收敛性,并能达到较高的计算精度.     

多层地基条带基础动力刚度矩阵的精细积分算法

提出应用精细积分算法计算多层地基的动力刚度问题. 精细积分是计算层状介质中波传播的高效而精确的数值方法. 利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组, 运用精细积分方法求解格林函数, 最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内, 进而得到刚性条带基础频率域的动力柔度或刚度矩阵. 所建议的精细积分算法, 可以避免一般传递矩阵计算中的指数溢出问题, 对各种情况有广泛的适应性, 计算稳定, 在高频段可以保障收敛性, 并能达到较高的计算精度.      

非线性动力学积分方程分块积分解法

对于非线性动力学方程组分块地应用精细积分算法,使其化成积分方程表达式,求解的表达式中具有相对低阶的转换矩阵,从而使精细积分更适用于多自由度、强非线性、变系数、非保守系统,针对积分方程提出了一个显示预测-校正的单步四阶精度自起步的精细积分算法。算例表明本方法是有效的。     

基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

精细时程积分法及其数值衍生格式应用评估

旋翼气动弹性耦合动力学方程本质上是一组刚性比较大的非线性偏微分方程。在有限元结构离散后,可改写为非齐次微分方程组,其中非齐次项是桨叶运动量(位移与速度)和气动载荷的函数。针对这类方程,本文尝试引入精细积分法及其衍生格式,借助数值方法计算Duhamel积分项。从积分精度与数值稳定性方面比较研究具有代表性的精细库塔法和高精度直接积分法。结合隐式积分算法,评估精细积分法应用于旋翼动力学方程的可行性。算例表明,精细积分法对矩形直桨叶动力学方程具有足够的求解精度。     

精细时程积分法及其数值衍生格式应用评估

旋翼气动弹性耦合动力学方程本质上是一组刚性比较大的非线性偏微分方程。在有限元结构离散后,可改写为非齐次微分方程组,其中非齐次项是桨叶运动量（位移与速度）和气动载荷的函数。针对这类方程,本文尝试引入精细积分法及其衍生格式,借助数值方法计算Duhamel积分项。从积分精度与数值稳定性方面比较研究具有代表性的精细库塔法和高精度直接积分法。结合隐式积分算法,评估精细积分法应用于旋翼动力学方程的可行性。算例表明,精细积分法对矩形直桨叶动力学方程具有足够的求解精度。     

求解非线性动力学方程的预测-校正数值算法

提出了一种求解非线性动力学方程的预测-校正数值算法.由于校正过程消除了线性化带来的误差,使计算精度得到较大提高.数值算例表明该方法正确、有效.     

饱和两相介质动力固结问题时域求解的精细时程积分方法

针对u-p形式的饱和两相介质波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例,对该方法的计算精度进行了校核。开展了该方法相关算法特性的研究,对采用不同数值积分方法计算非齐次波动方程特解项计算精度的差异进行了对比研究,并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明,(1)该方法具有良好的计算精度。(2)计算非齐次波动方程特解项的数值积分方法中,梯形积分法的计算精度最差,高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度。(3)增加高斯积分点数目对于提高计算精度的作用并不显著。     

层状地基任意形状刚性基础动力响应求解

提出了基于积分变换、对偶方程与精细积分算法求解多层地基任意形状刚性基础的动力刚度问题. 首先在频率波数域内圆柱坐标体系中利用圆形微元的对称与反对称特性建立多层地基中格林影响函数的波动方程,然后将应力和位移关系表示成对偶形式进行精细积分求解以提高计算精度和稳定性. 再将任意形状刚性基础与地基的交界面离散化为一系列圆形微元,利用格林影响函数建立其平动与转动动力刚度的矩阵方程. 该求解方法高效、准确并且计算稳定,适于任意复杂多层地基任意形状基础动力刚度的计算.      

Duhamel项的精细积分方法在非线性微分方程数值求解中的应用

基于Duhamel项的精细积分方法,构造了几种求解非线性微分方程的数值算法。首先将非线性微分方程在形式上划分为线性部分和非线性部分,对非线性部分进行多项式近似,利用Duhamel积分矩阵,导出了非线性方程求解的一般格式。然后结合传统的数值积分技术,例如Adams线性多步法等,构造了基于精细积分方法的相应算法。本文算法利用了精细积分方法对线性部分求解高度精确的优点,大大提高了传统算法的数值精度和稳定性,尤其是对于刚性问题。本文构造的算法不需要对线性系统矩阵求逆,可以方便的考察不同的线性系统矩阵对算法性能的影响。数值算例验证了本文算法的有效性,并表明非线性系统的线性化矩阵作为线性部分是比较合理的选择。     

精细积分时域平均法和随机扩阶系统法

讨论含随机参数结构的动力响应的计算问题,发展了精细积分时域平均法（ＴＡＰＩＭ）,它可以用来计算确定性系统受到随机激励时的动力响应;结合随机扩阶系统方法与随机有限元法,将ＴＡＰＩＭ方法应用于计算随机参数结构的动力响应,取得了较好的结果。结出了数值算例,结果表明随机扩阶系统法,随机有限元法与精细积分时域平均法的结合是计算 随机参数结构动力响应的有效方法。     

精细积分方法的发展与扩展应用

钟万勰院士于1991年首先提出计算矩阵指数的精细积分方法,其要点是2^N类算法和增量存储。精细积分方法可给出矩阵指数在计算机意义上的精确解,为常微分方程的数值计算提供了高精度、高稳定性的算法,现已成功应用于结构动力响应、随机振动、热传导以及最优控制等众多领域。本文首先介绍矩阵指数精细积分方法的提出、基本思想和发展;然后依次介绍在时不变/时变线性微分方程、非线性微分方程以及大规模问题求解中发展起来的各种精细积分方法,分析了其优缺点和适用范围;最后介绍了精细积分方法的基本思想在两点边值问题、椭圆函数和病态代数方程等问题的扩展应用,进一步展示了该思想的特色。     

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

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

一个求解浅水波方程的时域分段展开算法

发展了一种时域分段展开自适应方法求解一维非线性浅水波方程。通过时域分段展开,将一个非线性的时空耦合初边值问题转化为一系列的线性空间边值问题,并采用有限元方法递推求解;通过展开阶数的递进,实现了分段时域的自适应计算,当不同步长时可保持稳定的计算精度。研究结果表明,当步长较大而Heun’s法、四阶Runge-Kutta法不能得到合理结果时,本文算法仍能保证足够的计算精度。     

Direct Integration Methods with Integral Model for Dynamic Systems

ＩｎｔｒｏｄｕｃｔｉｏｎＴｒａｎｓｉｅｎｔｓｔａｔｅａｎａｌｙｓｉｓｈａｓｂｅｅｎａｎａｃｔｉｖｅｒｅｓｅａｒｃｈａｒｅａｉｎｍａｎｙｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ.Ｉｎｐｒａｃｔｉｃａｌｓｉｔｕａｔｉｏｎ ,ｌｉｋｅｓｔｒｕｃｔｕｒａｌｍｅｃｈａｎｉｃｓ,ｔｈｅｓｙｓｔｅｍｓｂｅｉｎｇｓｔｕｄｉｅｄａｒｅｕｓｕａｌｌｙｎｏｎｌｉｎｅａｒａｎｄｔｉｍｅ＿ｄｅｐｅｎｄｅｎｔ.Ｔｈｅａｎａｌｙｔｉｃａｌｍｅｔｈｏｄｓｔｏｓｏｌｖｅｔｈｅｓｅｐｒｏｂｌｅｍｓａｒｅｉｍｐｏｓｓｉｂｌｅ .Ｔｈｅｎｕｍ…     

基于黎曼几何的非线性振子动力学分析递推解析算法

基于黎曼几何和变分原理,推导了黎曼流形上非线性耗散动力系统的二阶微分动力学方程,并运用流形收缩的概念将动力学方程离散化,进而建立了相应的递推求解格式。选取3个自治非线性阻尼振子系统,分别采用递推解析算法和龙格库塔法求解微分动力学方程,并比较分析了不同的时间步长下两种算法的计算耗时。结果表明,与龙格库塔法相比,基于黎曼几何的递推算法不仅能得到每一时步的解析表达式,而且计算耗时短,计算效率高。基于黎曼流形的动力学方程递推算法为非线性动力学系统的解析求解提供了新思路。     

基于非线性状态空间辨识的气动弹性模型降阶

高维、非线性气动弹性系统的模型降阶是当前气动弹性力学与控制领域的研究热点之一.然而国内外现有的非线性模型降阶方法仍存在辨识算法复杂、精度有待提高等问题.本研究提出了一种基于非线性状态空间辨识的跨音速气动弹性模型降阶方法. 首先,该方法基于非定常空气动力的单位脉冲响应数据,采用特征系统实现算法对非线性状态空间模型的线性动力学部分进行系统辨识. 其次,引入状态和控制输入的非线性函数, 采用优化算法对非线性函数的系数矩阵进行优化,进而得到考虑非线性效应的空气动力降阶模型.为了验证该降阶模型在预测跨音速气动弹性力学行为的精确性,本文以三维机翼为研究对象,分别从基于非线性降阶模型的气动力辨识、跨声速颤振边界计算和极限环振荡预测三方面进行了算例验证,并与现有的模型降阶方法进行了对比, 进一步说明本文所提出方法的有效性.研究结果表明, 该降阶模型对上述三类问题的计算精度与直接流-固耦合方法相吻合,可用于高效预测飞行器跨声速气动弹性力学行为.      

非线性阻尼动力方程的复合积分法

非线性动力方程直接积分法的基础是构造$t$时刻与t+\Delta t时刻状态量间的关系, 由此形成基本量的非线性方程组, 再在每个时间步内采用 Newton-Raphson或BFGS等迭代方法求解. 该文基于Bathe复合积分法(composite implicit time integration), 提出了非线性阻尼系统基于速度变量的复合时间 积分迭代格式. 以非线性黏滞阻尼Sdof系统为例, 按上述方法以及基于BFGS迭代 的Newmark-\beta法编制Fortran程序, 结果与Adina软件对比, 验证了该文方法的有效性.     

On the limitations of linear beams for the problems of moving mass-beam interaction using a meshfree method

This paper deals with the capabilities of linear and nonlinear beam theories in predicting the dynamic response of an elastically supported thin beam traversed by a moving mass. To this end, the discrete equations of motion are developed based on Lagrange’s equations via reproducing kernel particle method (RKPM). For a particular case of a simply supported beam, Galerkin method is also employed to verify the results obtained by RKPM, and a reasonably good agreement is achieved. Variations of the maximum dynamic deflection and bending moment associated with the linear and nonlinear beam theories are investigated in terms of moving mass weight and velocity for various beam boundary conditions. It is demonstrated that for majority of the moving mass velocities, the differences between the results of linear and nonlinear analyses become remarkable as the moving mass weight increases, particularly for high levels of moving mass velocity. Except for the cantilever beam, the nonlinear beam theory predicts higher possibility of moving mass separation from the base beam compared to the linear one. Furthermore, the accuracy levels of the linear beam theory are determined for thin beams under large deflections and small rotations as a function of moving mass weight and velocity in various boundary conditions.     

多柔体系统响应计算的子循环计算方法研究

过去近30年中,柔性多体系统动力学研究取得了巨大的进展,人们的兴趣集中在柔性多体系统建模、计算及实验研究等3个方面. Belytschko等于1979年提出的子循环算法已经成功地应用于结构动力响应的有限元计算中,然而有关柔性多体动力学的子循环算法研究尚未见报道. 该文提出了一种适合于柔性多体系统响应计算的中心差分类子循环算法,在将非线性微分-代数混合方程组(DAEs)缩并为纯微分方程组(ODE)的基础上,推导出快、慢变分量的同步更新公式和子步更新公式;在变量的数值积分过程中,采用能量平衡计算检查算法的稳定性;算例结果表明该算法可以在保持合适的精度要求下,有效地提高响应的计算效率;对积分步长进行摄动修正可以保持算法的稳定性.