单点子域积分与差分

通过稳定性分析、显式与隐式积分，表明了单点子域积分相对于差分法的优越性．     

单点子域积分与多点子域积分

基于单点子域精细积分的思想,针对抛物线型热传导方程初边值问题,提出了多点子域积分的概念,推出了一种多点子域积分的FTCS格式。该格式为显格式,并证明其为无条件稳定。数值算例表明,多点子域积分的FTCS格式具有比单点子域积分的FTCS格式收敛速度快的特点。     

非线性动力学问题的一个显式精细积分算法

对非线性动力学问题,给出了一个显式的精细积分算法,适用于多自由度,强非线性,非保守系统,算例表明本方法精度高,计算量较少。     

二维扩散方程的单点子域精细积分法

建立了二维扩散方程的单点子域精细积分法,并通过稳定性分析,表明了单点子域精细积分法相对于差分法的优越性。     

泊松方程中域积分化为边界积分的方法

本文讨论了二维和三维泊松方程中域积分化为边界积分的方法。对于形如x~ig_x(y,z)、y~ig_x(x,z)和z~ig_z(x,y)的荷载给出了域积分转化为边界积分的正确公式。而对于复杂荷载,利用泰勒展开将域积分近似地转化为边界积分并给出了误差估计。计算结果表明利用本文方法可大大节省计算时间。因此,本文方法是一种十分有效的方法。     

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

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

结构动力响应分析的三阶显隐式时程积分方法

基于泰勒级数展开式提出了一种用于结构动力响应分析的高精度时程积分方法,该方法假设t时刻的速度和加速度由t?Δt时刻、t时刻、t+Δt时刻的速度和加速度加权表示,并可根据求解需要调节权值,将积分算法构造成隐式格式或显式格式.通过理论分析和数值算例,计算讨论了该算法的稳定性和精度,确定了最佳的权值和允许的时间步长.结果表明:本文算法最高具有三阶精度,且具有振幅衰减率低、周期延长率极小等优点.最后结合一个铁道工程实例,表明本文算法适用于大型非线性动态响应的精确快速求解.     

有限元方法中的光滑积分伪弱形式

提出了一种光滑积分伪弱形式,将光滑积分拓展至被积函数非偏导项求解。结合光滑应变技术和伪弱形式,可实现有限元系统方程统一光滑积分求解,即对刚度矩阵和质量矩阵中的应变矩阵和形函数矩阵均可进行光滑积分处理,并转化为光滑子域的边界积分。光滑积分伪弱形式与光滑应变技术比较,增加了形函数矩阵不定积分处理过程,且没有降低有限元求解对形函数连续性的要求。不过,伪弱形式改变了单元积分的求解形式,连续质量矩阵求解也无需坐标映射和雅可比矩阵计算。以轴对称二维问题为研究对象,结果表明极度不规则三角形和四边形单元光滑积分伪弱形式在静态和动态有限元方程求解中也具有很好的精度。     

一族新的高精度显式差分格式

对求解三维抛物型偏微分方程,利用待定参数法构造出一族新的高精度的三层显式差分格式,其精度为O(Δt3+Δx4+Δy4+Δz4),并论证了其稳定性,通过数值实例可见其精度较文[1]提高2位有效数字.     

色散方程的高稳定性两层四点显格式的单点精细积分法

基于单点精细积分的思想,对色散方程Ut=aUxxx构造了一类高稳定性的两层四点显式差分格式,其局部截断误差为O(τ+h),稳定性条件为|R|=|aτ/h3|≤f(β),其中f(β)对任意正实数β为单调递增函数.它们不仅显著地改善了同类格式的稳定性条件|R|≤0.25,而且也优于众多三层多点(5点或5点以上)显格式的稳定性条件.     

一种广义精细积分法

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

色散方程的高稳定性三层五点显格式的广义单点精细积分法

基于文[1]中的单点精细积分方法,对色散方程Ut=aUxxx提出了一种构造高稳定性三层五点(蛙跳)显格式的广义单点精细积分法.文中格式的局部截断误差为O(τ2+h2),而稳定性条件为|R|≤g(β)(其中g对任意正实数是单调递增函数),是同类格式中最好的[2].     

一种高效的局部径向基点插值无网格方法

提出了一种弹性动力分析的高效局部径向基点插值无网格方法(MLRPI).该方法采用径向基点插值形函数近似解变量,运用局部Petrov-Galerkin法推导出了相应的离散方程,并根据波动模拟的精度要求,得到某一结点的动力方程.然后采用Newmark常平均加速度法和中心差分法相结合的显式积分格式进行时域积分,得到每个自由度的一种解耦递推格式.最后,对一平面应变问题进行了求解,比较了该文提出的解耦MI.RPI方法、常规MLRPI方法和ANSYS有限元方法的精度和计算时间,结果表明解耦MLRPI方法与常规MLRPI方法的精度相当,但计算效率大大提高.     

基于微分求积法求解非粘滞阻尼系统的时程响应

研究一种含有指数型非粘滞阻尼线性多自由度振动系统的时程分析问题。该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表述为质点速度与核函数的卷积。由于阻尼模型的改变,常用的数值积分方法（如Newmark-β法、Wilson-θ法）不能直接应用于这种非粘滞阻尼系统。基于一种无条件稳定的微分求积方法,给出了这种非粘滞阻...     

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

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

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

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

精细辛算法的高效格式和简化计算

对精细辛几何算法设计了高效的迭代过程，减少了精细积分的计算量，同时提出了精细辛算法的简化形式，避免了复杂的矩阵求逆运算，并给出了相应的误差估计，最后编制了程序进行验证，证明了所采取的方法能够使计算快捷，精度高，稳定性好．     

获得热传导问题“拟解析解”的精细积分算法

研究并讨论了用于热传导分析有限元解的精细积分算法，算法很好地克服了传统方法求解时的单调性问题，且对空间离散后所获得方程的解是解析的，因而算法的解将具有“拟解析解”的意义，论文证明了算法单调性．