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

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

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

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

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

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

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

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

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

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

悬索桥颤振稳定性分析的精细时程积分法

研究精细时程积分法在悬索桥颤振稳定性分析中的应用,首先,将 气流整个作为一个系统,组集系统关于模态广义坐标的状态空间方程,然后,应用精细时程积分法计算状态的向量的时程响应,根据状态向量时程响应的对数衰减率判断系统的颤振稳定性,最后,以英国塞文悬索桥为数值算例,验证了本文方法的正确性。     

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

精细积分方法的评估与改进

详细分析了结构动力分析的精细积分方法的稳定性、计算精度,在此基础上提出了对现有精细积分方法的改进策略。算例证实了本文对精细积分方法改进的科学性与可行性。     

基于Householder方法的子域精细积分

有限元法生成的矩阵通常具有尺度和带宽很大的特点,应用子域精细积分面临着选择子域的困难。本文采用Householder方法将矩阵三对角化,从而使子域精细积分可用于大型有限元系统。通过在子域精细积分中引入预估-校正格式,可以在不需要迭代的情况下得到比蛙跳格式更高的精度。