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

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

瞬态热传导问题的精细积分-双重互易边界元法

采用双重互易边界元法结合精细积分法求解二维含热源的瞬态热传导问题。针对边界积分方程中热源项和温度关于时间导数项引起的域积分,采用双重互易法处理,将域积分转换为边界积分。采用边界元法将边界积分方程离散后,得到关于时间的微分方程组,并利用精细积分法处理其中的指数型矩阵;对于微分方程组中由边界条件和热源项引起的非齐次项,采用解析的方法计算。为了比较精细积分-双重互易边界元法的计算效果,同时使用有限差分法计算温度对时间的导数项。通过数值算例验证了本文方法的有效性和精确性。计算结果表明:时间步长对于精细积分-双重互易边界元法的结果影响较小,而有限差分法对时间步长比较敏感且只在时间步长选取较小时有效;当选取较大时间步长时,精细积分-双重互易边界元法依然具有良好的计算精度。     

一类非线性周期系统响应的精细积分法

对于一类非线性周期／变系数微分方程,提出基于精细积分法的数值解法,处理非线性周期／变系数微分方程系统的响应问题,其积分策略是：采用精细积分格式处理常系数部分;采用线性插值格式处理非线性周期／变系数部分,既继承精细积分方程高度准确的特点,又保证足够的精度与较小的计算量。通过数值算例,与以往与用的微分方程直接数值积分法（如预估－校正哈明法）求得的解加以比较表明,对于给定的精度要求,精细积分法更经济有效     

子域精细积分及偏微分方程数值解

对于偏微分方程半解析法的方程，精细时程积分虽然能求出高度准确的解，但往往面临矩阵尺度太大的困难，另一方面差分法虽然有带宽小的优点，但有稳定性及精度方面的问题，本文提出子域精细积分法，既可利用精细积分的数值优点，又有带宽小的好处，数值例题表明了子域精细积分法的效能。     

子域精细积分及偏微分方程数值解

对于偏微分方程半解析法的方程，精细时程积分虽然能求出高度准确的解，但往往面临矩阵尺度太大的困难；另一方面差分法虽然有带宽小的优点，但有稳定性及精度方面的问题．本文提出子域精细积分法，既可利用精细积分的数值优点，又有带宽小的好处．数值例题表明了子域精细积分法的效能．     

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

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

结构动力方程的精细与差分耦合时程积分法

提出一种将精细积分法与Newmark-β法耦合起来的结构动力学时程积分方法．该方法通过引入Newmark-β法的基本假设，将加速度分量从动力学方程中消去，动力学方程由二阶常微分方程组变为一阶常微分方程组，然后再用精细积分法进行逐步积分．与直接应用精细积分法相比，方程的个数可以减少一半．该文对这种方法进行了理论推导和算例验证，表明了该方法在结构动力分析中的有效性．     

一类指数矩阵函数及其应用

研究了一阶常微分方程组特解的精细积分方法. 针对非齐次项为多项式、指 数函数以及二者的乘积的情况，在Duhamel积分形式特解的基础上，引入了一类指数矩阵函 数. 通过该类函数的线性组合即可表达出非齐次方程的特解. 建立了该类指数矩阵函数的一 种高效递推算法，并在此基础上实现了特解的精细积分. 由于特解的积分过程能充分利用通 解精细积分过程的中间量，因此两个精细积分过程能有机地结合起来，形成了一种高效、统 一的广义精细积分法. 对上述递推算法做了进一步优化，并给出了通用的计算公式. 算例结果证明了该方法的有效性.     

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

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

非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

精细积分法在迷宫密封转子系统非线性动力学特性计算中的应用

基于Muszynska密封力模型，建立了迷宫密封转子系统的非线性动力学模型，将精细积分法推广应用于非线性情况，计算了迷宫密封不平衡转子系统的动力学特性，依据Floquet理论讨论其分岔特性。研究表明：在2^N类算法计算指数矩阵基础上提出的精细积分法和传统的数值计算方法相比，其精度高，在分析中通过取不同步长计算对比，表明该方法在某些情况下可以采取较大时间步长，有效提高了计算速度。     

瞬态热传导方程精细积分方法中对称性的利用

采用精细积分法求解瞬态热传导方程时,对指数矩阵进行变换后使其具有对称性,利用这一特性可使存贮量和计算量降低一半。变换后指数矩阵的带宽特性不变,采用子域精细积分可进一步提高算法的计算与存储效率。     

热传导问题灵敏度分析的伴随法

在热传导灵敏度分析的直接法的研究基础上，进一步探讨了稳态和瞬态热传导问题灵敏度分析的伴随法.推导了伴随法的计算列式，对于瞬态热传导问题，研究了瞬态约束处理的关键点方法，并提出伴随方程的精细积分解法。算例表明，稳态问题灵敏度计算，伴随法与直接法的结果是一致的；瞬态问题灵敏度计算，两种方法的精度相当。     

Timoshenko梁谐响应求解新方法探索

在空间域上采用只与结点有关的无网格方法离散，在时间域上采用精细积分方法求 解. 无网格离散过程中，利用伽辽金积分等效弱形式代替微分形式的控制方程，并 用修正变分原理满足位移边界条件，采用移动最小二乘法求解离散的形函数，把形 函数代入等效积分弱形式得到离散的二阶方程；精细积分过程中非齐次项采 用Romberg积分. 同时给出了两种不同边界条件的谐响 应求解的两个数值算例，得到了精确的数值结果.     

一种广义精细积分法

提出了求解非齐次动力方程特解的一种精细数值积分法，该方法与通解 精细积分法具有相同精度. 首先选取一个积分形式的非齐次方程特解，将积分区域划分为 2$^{N}$份，并对之进行精细的数值积分；然后针对载荷为多项式、指数函数及三角函数的情 况，将积分求和转化为一个递推过程，按此只需$n$次矩阵乘法就能计算出积分和，从而得到 非齐次方程的特解. 该方法的优点是能与通解的精细积分过程有机地结合起来，具有极高的 精度和效率，同时还具有较广泛的适用范围. 算例结果证明了该方法的有效性.     

ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

IntroductionThepreciseintegrationmethod(PIM) [1],whichwasproposedforsolvingstructuraldynamicequations.Thismethodissimplerandpossesseshigherprecision .Forlinearsteadystructuraldynamicsystems,itsnumericalresultsattheintegrationpointsarealmostequaltothatoftheexactsolutioninmachineaccuracy .InthepreciseintegrationmethodforsolvingPDEs,theequationsshouldbediscretizedinthephysicalspaceforobtainingthesystemofODEsintime ,whichisoftenexecutedbythefinitedifferencemethodorthefiniteelementmethod .Inrec…     

结构动力方程的增维精细积分法

对线性定常结构动力系统提出的精细积分方法,能够得到在数值上逼近于精确解的结果,但对于非齐次动力方程涉及到矩阵求逆的困难。提出采用增维的办法,将非齐次动力方程转化为齐次动力方程,在实施精细积分过程中不必进行矩阵求逆,这种方法对于程序实现和提高数值稳定性十分有利,而且在大型问题中计算效率较高,从而改进了精细积分方法的应用,数值例题显示了本文方法的有效性。     

Precise integration method for LQG optimal measurement feedback control problem

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

瞬态热传导问题的精细积分数值流形方法研究

数值流形方法（NMM）因其特有的双覆盖系统（数学覆盖和物理覆盖）在域离散方面具有独特的优势，而精细时间积分法则具有精度高、无条件稳定、无振荡以及计算结果不依赖于时间步长等特点。发展了用于研究二维瞬态热传导问题的精细积分NMM。结合待求问题的控制方程和边界条件，并基于修正变分原理导出了NMM的总体方程，给出了求解此类时间相依方程的精细时间积分及空间积分策略，选取了两个典型算例对方法的有效性进行了验证，结果表明本文方法可以高效高精度地求解瞬态热传导问题。