共查询到17条相似文献,搜索用时 156 毫秒
1.
一维抛物型偏微分方程可以用精细积分方法精确求解。当精细积分中的矩阵指数函数用Pade逼近来代替时,可以得到一系列由简到繁、精度由低到高的差分格式,因而便于根据实际需要进行选取。常见的求解抛物型方程的差分格式如古典显式格式、隐式格式及六点差分格式为其中的特例。Pade逼近格式主要包括矩阵运算和线性方程组求解。本文利用Pade逼近格式对应的方程组系数矩阵为带状矩阵的特点,把原来在整个区域上求解的问题转化为分区域求解,在TRANSPUTER并行机上实现了该问题的并行算法,并对该并行算法的时间复杂度进行了分析。算例结果表明Pade逼近并行算法有很好的计算效果和并行效率。 相似文献
2.
结构动力方程的更新精细积分方法 总被引:26,自引:3,他引:26
将高斯积分方法与精细积分方法中的指数矩阵运算技巧结合起来,建立了精细积分法的更新形式及计算过程,对该更新精细积分方法的稳定性进行了论证与探讨。在实施精细积分过程中不必进行矩阵求逆,整个积分方法的精度取决于所选高斯积分点的数量。这种方法理论上可实现任意高精度,计算效率较高,其稳定性条件极易满足。数值例题也显示了这种方法的有效性。 相似文献
3.
矩阵指数精细积分方法中参数的自适应选择 总被引:2,自引:1,他引:1
讨论了基于Pad逼近的矩阵指数精细积分方法中加权系数N和展开项数q的自适应选择问题.参数(N,q)的选择直接影响到矩阵指数计算的精度和效率.采用矩阵函数逼近理论,研究了参数Ⅳ和q的增加对精度的影响程度,据此,提出了参数(N,q)优化组合的递推自适应选择方法.该方法可以根据矩阵本身的性态选择合适的参数(N,q),而参数选择的计算量与矩阵指数的计算量相比几乎可以忽略,这对于增强矩阵指数精细积分方法的适应性和提高计算效率是很有益处的.算例验证了该方法的正确性和有效性. 相似文献
4.
一维抛物型偏微分方程可以用精细积分方法精确求解.当精细积分中的矩阵指数函数用Padé逼近来代替时,可以得到一系列由简到繁、精度由低到高的差分格式,因而便于根据实际需要进行选取.常见的求解抛物型方程的差分格式如古典显式格式、隐式格式及六点差分格式为其中的特例.Padé逼近格式主要包括矩阵运算和线性方程组求解.本文利用Padé逼近格式对应的方程组系数矩阵为带状矩阵的特点,把原来在整个区域上求解的问题转化为分区域求解,在TRANSPUTER并行机上实现了该问题的并行算法,并对该并行算法的时间复杂度进行了分析.算例结果表明Padé逼近并行算法有很好的计算效果和并行效率. 相似文献
5.
Shujun Tan Zhigang Wu Wanxie Zhong 《力学学报》2009,41(6):961
讨论了基于Pad\'{e}逼近的矩阵指数精细积分方法中加权系数N
和展开项数q的自适应选择问题. 参数(N,q)的选择直接影响到矩阵指数计算的精度和效
率. 采用矩阵函数逼近理论,研究了参数N和q的增加对精度的影响程度,据此,提出了
参数(N,q)优化组合的递推自适应选择方法. 该方法可以根据矩阵本身的性态选择合适的参
数(N,q),而参数选择的计算量与矩阵指数的计算量相比几乎可以忽略,这对于增强矩阵指
数精细积分方法的适应性和提高计算效率是很有益处的. 算例验证了该方法的正确性和有效性. 相似文献
6.
一维抛物型偏微分方程可以用精细积分方法精确求解。当精细积分中的矩阵指数函数用Pad 逼近来代替时 ,可以得到一系列由简到繁、精度由低到高的差分格式 ,因而便于根据实际需要进行选取。常见的求解抛物型方程的差分格式如古典显式格式、隐式格式及六点差分格式为其中的特例。Pad 逼近格式主要包括矩阵运算和线性方程组求解。本文利用 Pad 逼近格式对应的方程组系数矩阵为带状矩阵的特点 ,把原来在整个区域上求解的问题转化为分区域求解 ,在 TRANSPUTER并行机上实现了该问题的并行算法 ,并对该并行算法的时间复杂度进行了分析。算例结果表明 Pad 逼近并行算法有很好的计算效果和并行效率。 相似文献
7.
8.
9.
一类指数矩阵函数及其应用 总被引:3,自引:0,他引:3
研究了一阶常微分方程组特解的精细积分方法. 针对非齐次项为多项式、指
数函数以及二者的乘积的情况,在Duhamel积分形式特解的基础上,引入了一类指数矩阵函
数. 通过该类函数的线性组合即可表达出非齐次方程的特解. 建立了该类指数矩阵函数的一
种高效递推算法,并在此基础上实现了特解的精细积分. 由于特解的积分过程能充分利用通
解精细积分过程的中间量,因此两个精细积分过程能有机地结合起来,形成了一种高效、统
一的广义精细积分法. 对上述递推算法做了进一步优化,并给出了通用的计算公式.
算例结果证明了该方法的有效性. 相似文献
10.
非齐次动力方程Duhamel项的精细积分 总被引:13,自引:1,他引:13
提出了不需要矩阵求逆运算的求解Duhamel积分项的精细积分方法.通过将精细积分法的关键思想--加法定理和增量存储--直接应用于Duhamel积分响应矩阵的求解,可给出当非齐次项分别为多项式、正弦/余弦以及指数函数等基本形式时Duhamel积分在计算机上的精确解.特别的,该算法不依赖于系统矩阵(或相关矩阵)的形态.当系统矩阵奇异或接近奇异时,其优越性更为显著.算例验证了该算法的有效性. 相似文献
11.
关于动力分析精细积分算法精度的讨论 总被引:9,自引:3,他引:6
对动力问题分析的精细积分算法的精度问题进行深入研究,并在此基础上提出对原有的算法的改进策略,改进后的算法可以较好地克服算法精度对积分时间步长的依赖性问题。 相似文献
12.
精细积分方法的评估与改进 总被引:8,自引:1,他引:8
详细分析了结构动力分析的精细积分方法的稳定性、计算精度,在此基础上提出了对现有精细积分方法的改进策略。算例证实了本文对精细积分方法改进的科学性与可行性。 相似文献
13.
本文在双连杆空间柔性机械臂系统非线性动力学方程的基础上,运用线性二次型(LQ)最优控制方法讨论了机械臂消除残余振动的控制问题。本文重点在于系统计算过程中,放弃了传统的差分类算法,对时变控制系统,引入时程精细积分方法,由于精细积分方法在有限的时间步长内又进行了更精细的划分,同时避免了差分法的许多计算障碍,使得该计算方法具有计算精度高及数值计算无条件稳定等特点。文中针对双连杆空间柔性机械臂系统这一典型结构,给出了其精细控制律,以说明精细积分法的优越性。 相似文献
14.
15.
16.
This paper describes an accurate and efficient method for calculating the first and second derivatives of dynamic response with respect to design variables for linear structural systems subjected to transient loads. An efficient algorithm to calculate the dynamic responses and their first and second derivatives is formulated based on Gauss precise time step integration method. The algorithm is achieved by direct differentiation and only a single dynamic analysis is required. Several numerical examples are comparatively demonstrated using the new developed method, analytical method, and central difference method. The results show that the new method is highly accurate compared with the analytical approach and is more efficient than the central difference method. 相似文献
17.
直接积分法是求解动力学方程的一种有效方法。应用一种预估-校正的Generalized-α法对结构大变形动力学问题进行分析求解,并与Newmark法和Bathe法进行对比研究。首先预估当前计算步的解,然后以预估值作为起始值进行非线性迭代计算,并对解不断校正,直到满足收敛条件,进入下一时间步的计算。在保证Generalized-α法性能的基础上,简化了非线性迭代公式,便于编程实现。通过壳和实体的大变形动力学算例,证明了本文方法具有较高的稳定性和精度。 相似文献