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

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

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

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

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

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

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

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

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

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

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

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

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

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

增维精细积分法的适用范围研究

对线性定常结构动力系统提出的增维精细积分法,能够将非齐次动力方程转化为齐次动力方程,不用对状态矩阵求逆就能方便高效地求解出结构的动力响应。本文在仔细分析增维精细积分法性质的基础上,提出了其适用条件,进一步拓宽了其应用范围,并给出了将荷载项展开成傅里叶级数时,相应增维精细积分法的表达式。同时,在一个时间步长内,通过对非齐次项作线性化假设,成功地将增维精细积分法应用到了非线性动力分析领域。本文方法计算格式统一,易于编程,具有很高的计算效率。数值算例证明了本文方法的有效性。     

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

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

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

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

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

基于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的总体方程,给出了求解此类时间相依方程的精细时间积分及空间积分策略,选取了两个典型算例对方法的有效性进行了验证,结果表明本文方法可以高效高精度地求解瞬态热传导问题。