精细时程积分的算法改进研究论文.pdfVIP

第十三届全国工程建设计算机应用学术会议论文集 广东佛山·2006 精细时程积分的算法改进研究 徐腾飞 向天宇赵人达 (西南交通大学桥梁工程系 四J1I成都610031) 【摘要】精细时程积分算法涉及到稠密矩阵，在计算大型结构中，计算规模大，需要大量的存储空间。 根据矩阵结构的特殊性以及指数函数的性质来计算指数矩阵，给出了精细时程积分的改进意见，使其更便 于在大型结构动力响应中的推广应用。算例证明了改进的正确性和可行性。 【关键词】结构工程；精细时程积分算法：动力响应；大型结构 1 引言， 钟万勰教授在文献【l】中提出了结构动力方程求解的精细时程积分法，结合指数矩阵的 精细算法能够获得高度精确结果，计算误差与计算机截断误差同量级，同时还具有无条件 稳定等优点【34J。精细时程法在将二阶微分方程降为一阶的同时，将r／阶方阵扩充为2咒阶方 阵，增大了计算的规模。在实际应用中在存储空间和运算时间两方面都要付出很高的代价， 这使精细时程法在大型结构动力响应分析中的应用受到了限制。 本文利用矩阵结构的特殊性与矩阵乘方的性质，改进算法，降低了精细时程积分的计 算规模，使其更便于在大型结构动力响应中的推广应用。 2精细时程积分方法 离散结构模型的动力学方程为 m肛)+[c肛)+【髟】{x)=妒(f)) (1) 。 ’{x(o)}=(x。)伍(o))=伍。) (2) 式中【M】，【C】和Ⅸ】分别表示质量矩阵、阻尼矩阵和刚度矩阵，{F(t))为给定的外部激励列 向量。则方程(1)可以表示为 p)=陋p)+{r) (3) 式中 (4) 式(3)称为结构动力响应的状态方程，其一般解为 {v(f))=etZlt{％)+fce。[ttlO-r){，．(力}df (5) 76 第十三届全国工程建设计算机应用学术会议论文集 广东佛山·2006 到递推公式 如(f。))=【丁()】{V(f。)}+r∥‰训{r(f)}df (6) 精细积分方法利用指数函数的加法定理，并做泰勒展开 r丁c△口，，=已cH，山=’[P【Ⅳ号]2“=c已c玎，r，2Ⅳ，f=多 (7) ．㈣，_--elXlr=Ⅲ嘻学 (8) 【Z】=2【巧一．】+眵。】·【Z_J“=1,2，…，Ⅳ) (9) 【T】=【，】+【瓦】 (10) (11) J=【丁】(fVI}+【Ⅳ】一1({％J+【Ⅳ】一1{‘1)卜【日r1【{％)+【日】-1{‘l+{‘l纠 ．{VM 这就是线性的非齐次项时的精细积分公式，被称为HPD．L格式‘l5】，此外针对非齐次项不同的 积分方法还有精细积分法的HPD．s格式‘41，高斯精细积分方法(GPTSM)‘51等等。 3精细时程积分算法的改进 在n个自由度的问题中，【明矩阵为2，z阶的方阵。而在文献【2】算法，要对【明进行多次乘 法运算。通过观察可以发现[明矩阵上半部分由O矩阵与单位矩阵构成。 记【明矩阵的i次方为 ．。