矩阵指数精细积分方法中参数的自适应选择

讨论了基于Pad逼近的矩阵指数精细积分方法中加权系数N和展开项数q的自适应选择问题.参数(N,q)的选择直接影响到矩阵指数计算的精度和效率.采用矩阵函数逼近理论,研究了参数Ⅳ和q的增加对精度的影响程度,据此,提出了参数(N,q)优化组合的递推自适应选择方法.该方法可以根据矩阵本身的性态选择合适的参数(N,q),而参数选择的计算量与矩阵指数的计算量相比几乎可以忽略,这对于增强矩阵指数精细积分方法的适应性和提高计算效率是很有益处的.算例验证了该方法的正确性和有效性.     

大规模动力系统改进的快速精细积分方法

提出一种针对大规模动力系统的改进的快速精细积分方法（FPIM）。以精细积分方法为基础,利用大规模动力系统矩阵的稀疏性和动力问题的物理特性,分析了矩阵指数的特殊结构,并基于此给出一种计算大规模动力系统矩阵指数及其动力响应的高效率方法。     

Adaptive selection of parameters for precise computation of matrix exponential based on padé approximation

讨论了基于Pad\'{e}逼近的矩阵指数精细积分方法中加权系数N 和展开项数q的自适应选择问题. 参数(N,q)的选择直接影响到矩阵指数计算的精度和效 率. 采用矩阵函数逼近理论，研究了参数N和q的增加对精度的影响程度，据此，提出了 参数(N,q)优化组合的递推自适应选择方法. 该方法可以根据矩阵本身的性态选择合适的参 数(N,q),而参数选择的计算量与矩阵指数的计算量相比几乎可以忽略，这对于增强矩阵指 数精细积分方法的适应性和提高计算效率是很有益处的. 算例验证了该方法的正确性和有效性.     

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

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

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

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

周期结构动力响应的高效数值方法

提出一种计算周期结构动力响应的高效率算法. 以精细积分方法为基础, 利用周期结构的对称性和动力问题的物理特性, 分析了周期结构对应矩阵指数的特殊结构, 并基于此给出一种计算周期结构对应矩阵指数的高效率方法. 在高效和精确计算周期结构对应矩阵指数的基础上, 得到了周期结构动力响应的高效率和高精度算法. 数值算例表明, 该方法效率高且节省存储要求.      

用等效力系变换矩阵求解静不定桁架极限载荷

采用等效力系变换矩阵研究了双模量静不定桁架极限载荷问题．首先证明了固体的等效力系变换矩阵与等效位移变换矩阵是互为转置的矩阵,采用等效力系变换矩阵求解双模量静不定桁架结构的内力,然后再利用静力方程确定双模量静不定桁架结构的极限载荷．当力的变换关系可以根据物理条件容易求得,而位移的变换关系不容易找出时,用等效力系变换矩阵求解静不定桁架极限载荷,就更能显示出其计算过程简洁、清晰等优点．用等效力系变换矩阵求解静不定桁架极限载荷不涉及材料的性质,对各向同性材料、双模量材料静不定桁架极限载荷的求解都适用．     

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

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

非线性动力方程的三次样条-增维精细算法

提出一种针对非线性动力方程的改进精细积分方法。该方法是在时间步长内采用分段的三次样条函数拟合非齐次项,保持高精度拟合的同时避免了求导运算和高次多项式插值带来的Runge现象。通过引入4×2个变量将动力方程增加四维转化为齐次方程,并建立相应的通解格式,避免了状态空间下系统矩阵求逆。将指数矩阵分为四个子模块,利用各模块的特点分别进行理论推导及基于精细积分法进行分步、分块计算得到相应的理论解和高精度数值解,无需反复计算整个指数矩阵,提高了解算效率。针对含未知状态量的非齐次项,引入预测-校正的方法进行迭代求解。数值计算结果表明了本文方法的有效性。     

结构动力方程的更新精细积分方法

将高斯积分方法与精细积分方法中的指数矩阵运算技巧结合起来,建立了精细积分法的更新形式及计算过程,对该更新精细积分方法的稳定性进行了论证与探讨。在实施精细积分过程中不必进行矩阵求逆,整个积分方法的精度取决于所选高斯积分点的数量。这种方法理论上可实现任意高精度,计算效率较高,其稳定性条件极易满足。数值例题也显示了这种方法的有效性。     

瞬态热传导方程的子结构精细积分方法

对线性定常结构动力系统提出的精细积分方法,在数值精度等方面表现出极大优越性,但是当矩阵尺度很大时在数值计算与存储中将产生困难,对此,本文对瞬态热传导方程,根据结构的概念,将结构分为若干个子结构,对各子结构分别进行指数矩阵运算并通过了结构间界面的物理量相联系,从而提高精细积分方法的计算效率。     

?????????????????????Krylov?????????

提出了一种新的精细时程积分法来求解大型动力系统. 结合Krylov子空间法、培德级数 近似以及一般载荷的维数扩展法,进一步提高精细时程积分法的计算效率. 利用维数扩展法 避免计算微分方程特解,并可处理任意载荷. 对于大型动力系统,通过Krylov子空间的降维 分析将问题转化到一个子空间,计算效率得到极大提高. 对于迭代次数$N$ 的选择作了详细讨论,进一步提高了计算效率.     

功能梯度材料瞬态热传导问题的降维精细积分法

针对功能梯度材料(FGMs)的二维瞬态热传导问题,提出了一种有效的降维精细积分法。通过空间坐标的分步离散,将原二维问题的精细积分转化为一系列一维问题的精细积分,解决了直接精细积分算法所面临的传递矩阵 规模过大的问题,从而大幅度降低了计算量和存储量。算例表明了本文方法的有效性。     

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

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

大型结构动力响应的状态方程的Krylov精细时程积分法

提出了一种新的精细时程积分法来求解大型动力系统.结合Krylov子空间法、培德级数近似以及一般载荷的维数扩展法,进一步提高精细时程积分法的计算效率.利用维数扩展法避免计算微分方程特解,并可处理任意载荷.对于大型动力系统,通过Krylov子空间的降维分析将问题转化到一个子空间,计算效率得到极大提高.对于迭代次数N的选择作了详细讨论,进一步提高了计算效率.     

无条件稳定的更新精细积分方法

提出将Pade逼近与精细积分方法中的指数矩阵运算技巧结合起来,建立了精细积分法的更新形式及计算过程,对该更新精细积分方法的稳定性进行了论证与探讨.结果表明,该更新精细积分方法是无条件稳定的,整个积分方法的精度取决于所取Pade逼近的阶数与高斯积分点的数量.数值例题也显示了该方法的高效率及其可行性.     

非齐次动力方程Duhamel项的精细积分

提出了不需要矩阵求逆运算的求解Duhamel积分项的精细积分方法.通过将精细积分法的关键思想--加法定理和增量存储--直接应用于Duhamel积分响应矩阵的求解,可给出当非齐次项分别为多项式、正弦/余弦以及指数函数等基本形式时Duhamel积分在计算机上的精确解.特别的,该算法不依赖于系统矩阵(或相关矩阵)的形态.当系统矩阵奇异或接近奇异时,其优越性更为显著.算例验证了该算法的有效性.     

一种广义精细积分法

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

Fast precise integration method for hyperbolic heat conduction problems

A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suit- able for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.