求解非线性动力方程的一种齐次扩容精细积分法

提出了求解非线性动力方程的一种齐次扩容精细积分法.首先利用泰勒公式将动力方程的非线性部分在tk时刻展开至二阶或更高阶级数,然后将1,(t-tk)和(t-tk)2/2等扩充到状态方程中,建立了便于用精细积分法计算的齐次方程形式.该方法能有效避免系统矩阵的求逆问题,且在保证具有较高计算精度的前提下,能使积分步长有效拓宽,提高了计算效率.为适应实际计算,还提出了一种通过迭代修正间接计算导数的方法.计算结果表明所提出的方法具有较好的计算精度和可靠性,是一种求解非线性动力方程的有效方法.     

用精细积分法求解非线性薛定谔方程

文章将精细积分法从求解线性定常结构动力系统推广应用于求解非线性薛定谔方程上.首先将非线性薛定谔方程变形为齐次方程的形式,然后用精细积分法模拟其随时间的演化过程.具体模拟了变系数非线性薛定谔方程的解,给出了两周期性孤子的相互作用情况及两个光脉冲的相互干涉情况.通过具体算例,说明该方法是可以对非线性薛定谔方程所反应的问题进行模拟计算的,并且有切实的效果.     

多自由度变弹性系统动力响应的多步积分法

采用微分方程的多步数值积分方法——亚当姆斯法来计算多自由度弹性系统的强迫动力响应。     

精细积分法在非线性动力学问题中的应用

 针对非齐次结构动力方程Duhamel形式的特解,建立了一种高效的特解精细积分法,对于非齐次项为幂函数和指数函数的情况,该方法能给出计算机上最高精度的解答。上述特解精细积分过程能与通解精细积分过程有机地结合起来,并形成一种高效的广义精细积分法。在此基础上,建立了非线性动力学方程的一种迭代算法。该方法具有很高的精度和效率以及较大的适用范围。算例结果证明了该方法的有效性。     

二维抛物型方程精细积分法与差分法比较

可用单内点子域精细积分法,求解二维抛物型方程初值问题。当单内点精细积分中的传递函数即指数函数用Taylor展开式的一阶近似来替代时,精细积分转化为差分方程。研究这一对应关系,使各种常见差分格式均找到对应的单点精细积分格式,并在单点精细积分的一般公式中获得统一表达式。     

结构非线性动力学方程的自适应精细积分算法

将定常结构动力方程的精细时程积分算法推广应用于非线性动力学问题时,对非线性项的线性化处理使该方法的计算精度对时间步长非常敏感.为此本研究中将龙贝格积分法引入该方法,提出了由此而产生的指数矩阵的快速精细算法,从而使时间步长的选择具有了自适应性,且计算精度和效率均得到提高.文中的数值算例给出了该方法的计算精度和效率.     

非线性动力学方程的自适应精细积分

将定常结构动力方程的精细积分算法推广应用于非线性动力学问题的求解．对非线性项的线性化处理使该方法的计算精度对时间步长非常敏感,为此将龙贝格积分法引入该方法,提出了由此而产生的指数矩阵的快速精细算法,从而使时间步长的选择具有了自适应性,计算精度和效率均得到提高。     

一种提高增维精细积分法计算精度的方法

提出了一种提高常系数非齐次常微分方程组增维精细积分法计算精度的方法。对于常系数非齐次常微分方程组,一般增维精细积分方法在每一个时间步内把非齐次项当成常数,并取其值为该时间步的初始值。在每一个时间步长内,仍然将非齐次项当成常数,但是该常数的值取为该时间段内不同时刻值的平均值,或者取为中间时刻的值,计算精度得到了很大的改善。数值算例显示了方法的有效性。     

非线性动力方程精细积分级数解的并行算法

非线性动力方程通过变量变换可以转化为一阶微分方程,该方程的解由表示初值影响的齐次方程解和反映荷载作用的积分之和组成.其中:第一项用指数矩阵计算;第二项在文中采用级数解计算(设计了3种相应的并行算法),算法1对级数解的每一项先做若干个向量的线性组合,再做矩阵向量乘1次;算法2与算法1原理相同,只是将矩阵的幂运算转换成乘积;算法3先做若干个矩阵向量乘,再做若干个向量的线性组合.算法1的并行效率最好,但存储空间需求大,不利于大型结构的求解.算法2、3利用动力方程的稀疏变换改善了算法1的不足,算法3中级数解每一项计算均在其前一项基础上进行,一般能比算法2节省时间.最后,给出了算例验证,三种算法都获得了较好的加速比.     

结构动力方程的精细积分-FFT方法

 离散结构动力方程为差分方程,并假设始、末时刻位移已知,将初值问题形式上转化为边值问题,然后利用快速傅立叶变换(FFT)进行求解,从而得到非齐次结构动力方程的一个数值特解。将该数值特解与通过精细积分法求得的齐次方程通解相结合,建立了求解结构动力方程的一种新方法。该方法具有较高的精度和计算效率,算例的数值结果证明了本文方法的有效性。     

一种改进的精细-龙格库塔法

 提出了求解非线性动力学方程的一种改进的精细龙格库塔法。首先对于线性问题,利用等步长的Newton-Cotes积分公式计算非齐次方程Duhamel积分形式的特解。由于在此过程中提出了一种简便的算法,与常规的同精度数值积分法相比,能较大程度地降低计算量和存储量。然后将上述方法推广到非线性问题,对于各积分点上未知的状态参量,参照龙格库塔法的几何解释进行一次预估。与已有的精细龙格库塔法相比,在精度和效率上均有较大程度的改善。算例结果充分证明了该方法的有效性。     

扩散方程的子域精细积分恒稳定的隐式差分格式

基于子域精细积分的思想,提出了一个求解扩散方程初边值问题的含参数a＞0恒稳定的隐式差分格式.它的局部截断误差为0(ατ2 h2),其中α＞0,τ和h分别为时间步长和空间步长.文末的数值实验表明,该方法有较高的精度和良好的实用性.     

对流方程的样条子域精细积分(SSPI)格式

针对对流方程第一类初边值问题,基于子域精细积分的思想,结合三次样条函数逼近,提出一个含参数α(α>0)无条件稳定的样条子域精细积分(SSPI)格式,并进行数值实验.SSPI格式求解对流方程有效,而且局部截断误差为O(ατ2 τ2 h4).SSPI格式不仅能够求解对流方程的第一类边值问题,而且能够求解第二类、第三类初边值问题,是一种有效的算法.     

对流-扩散方程的样条子域精细积分隐格式

针对对流-扩散方程的初边值问题, 利用子域精细积分的思想, 结合三次样条函数逼近, 提出含参数(α>0)的一族无条件稳定的隐格式,其局部截断误差阶为O(ατ+τ2+h2).当参数0<α≤τ时,其精度相当于O(τ2+h2), 且可用三对角线追赶法容易地求解. 数值计算表明,理论分析与实际例子相符合.     

流固耦合系统动力响应分析的精细时程积分法

通过对结构与理想流体耦合问题的分析 ,利用有限元方法对流固耦合系统动力响应进行了研究 .采用精细时程积分法、威尔逊θ法和纽马克法进行计算 .算例表明 ,精细时程积分方法具有精度高、不受时间步长的严格限制和计算工作量小等优点 ,适合于流固耦合系统的动力响应分析     

高维非线性Cahn—Hilliard方程的谱方法

考察一类非线性Ｃａｈｎ－Ｈｉｌｉａｒｄ方程的谱方法,构靠了一类有条件稳定的半离散和全离散格式,采用先验估计和Ｓｏｂｏｌｅｖ不等式,证明有了其格式的收敛性与稳定性。     

Precise Integration Method of C-integral for Linear Viscoelastic Solids with Crack

Generally, viscoelastic solid constitutive equation can be written into differential form and integral form. From differential form of constitutive equation the mathematical representation of the state space equation can be derived. Due to the state space equation, general constitutive equation can be solved by precise integration method that can be used in many fields with the advantages of highly precision and convenience. For linear viscoelastic solids with crack, the finite elements program of the precise integration method for viscoelastic solid is developed, which appears to be efficient and precise. C-integral is used to be the characterising parameter.