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

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

线形装药强爆炸地震反应谱与地震力计算

本文根据花岗岩、砾岩、黄土三种不同介质的近场爆炸地震加速度的实测资料和单层及多层建筑物的实测数据。计算出线形装药近区强爆炸地震加速度反应谱,并给出爆炸地震力的计算方法。     

柔性机械臂碰撞系统的一种建模方法及精细积分数值模拟

对一单连杆柔性机械臂与外界发生碰撞时的动力学特性进行了研究。在建立动力学模型时，利用有限元思想，用有限个刚性单元和抗扭连接弹簧来模拟机械臂及其柔性性质，得出了较为简单的非线性动力学方程，根据几何关系，建立了系统在发生碰撞时的约束条件。本文在系统的动力学方程中引入对偶变量，并将方程导向Hamilton体系，用精细积分法对系统的动力学特性进行了模拟。     

线形装药强爆炸地震反应谱与地震力计算

本文根据花岗岩、砾岩、黄土三种不同介质的近场爆炸地震加速度的实测资料和单层及多层建筑物的实测数据。计算出线形装药近区强爆炸地震加速度反应谱,并给出爆炸地震力的计算方法。     

矩阵黎卡提方程的精细积分法

选择恰当的参数，将２￣Ｎ类算法用于代数与微分黎卡提方程。证明了算得的解是如此精确，几乎是计算机上的精确解。数例验证了该结论。     

非线性动力方程精细积分法的自适应步长研究

基于Adams显式和隐式预估公式实现对时间步长的 自适应选择,利用当前时刻v(tk),采用预估公式的两种形式(显式与隐式),对v(tk+1)进行两次预估,利用两公式局部截断误差关系,得出误差估计值ξ(tk+1),并根据其大小 自适应调节时间步长.将该思想应用于预估型(求解过程需要用到预估公式)精细积分算法中,使精细积分...     

时滞抛物型方程的高精度精细积分法

对一类时滞抛物型方程初边值问题,提出了关于空间步长是四阶精度的高精度无条件稳定的精细积分法.数值算例表明,本文提出的精细积分法具有很高的精度,因而是一种有效的数值方法.     

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

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

隧洞开挖爆破地震反应谱特征研究

利用地震工程中应用比较成熟的反应谱理论对隧洞开挖爆破地震进行研究。针对反应谱计算中必须输入加速度数据的问题,对实测爆破振动速度数据进行数值微分计算,并采用小波分析方法消噪得到准确、清晰的加速度时程曲线,继而将其用于反应谱计算。结合不同工程中大量的爆破振动监测资料,对隧洞开挖爆破地震反应谱进行分析。结果表明:隧洞开挖爆破地震反应谱都有一些共同的特征,反应谱曲线一般存在一个明显的峰值点,峰值对应的周期集中在几毫秒到几十毫秒的区间内。另外,在归纳、整理大量地下工程开挖爆破地震反应谱的基础上,得到了隧洞开挖爆破地震的设计反应谱。     

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

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

桥梁受移动荷载动力响应的一种精细积分法

迄今为止,精细积分法都是应用于荷载作用位置固定不变的问题。本文将精细积分法推广到荷载作用点位置随时间而变化的情形,按有限元方法计算了桥梁受移动荷载作用时的动力响应,并与常规的Newmark方法、解析解做了比较。数值结果表明,精细积分法按本文的策略推广后,计算精度和效率均比通常的数值积分方法得到显著的提高。     

精细积分法在电报方程求解中的应用

将精细积分法应用到了二维的电报方程的数值计算之中。实例计算表明，该方法具有简单、计算精度高、无条件稳定、不需要进行复杂、费时的频域一时域转换及卷积积分，直接时域分析，处理非零初始值容易等优点。与传统的ＦＦＴ法及ＮＩＬＴ法相比，其效率更高，功能更强。     

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

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

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

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

ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

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

Characteristic Galerkin method for convection-diffusion equations and implicit algorithm using precise integration

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability. The project sponsored by the State Scientific and Technological Commission of China through “China State Key Project: the Theory and Methodology for Scientific and Engineering Computations with Large Scale”, the National Natural Science Foundation of China and the European Commission Research Project CI1^*CT94-0014.     

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

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

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

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

模拟相变问题的精细积分边界元法

将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。     

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

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