暂态历程的精细计算方法.pdfVIP

维普资讯 第12卷第 】期 计算结构力学及其应用 Vo1．12No．1 1995皇2月COMPUTATIONALSTRUCTURALMECHANICS PLIc皇! Fe—brua—ry1995 I_ ： } f．』 ，／z( 暂态历程的精细计算方法 要 准确的条件 关键词 堑生互 ／时程积分；2类算法 和 红， ，晔缓 分娄号 0241．4 鲫 婀 1 关于暂态计算的方程 暂态历程的计算方法是广泛关心的课题 ～ 虽然暂态分析往往要处理非线性或时变系 统 ，但是其基础仍是首先应精细地将线性时不变系统的分析做好，然后就可 用不同的方法进 一 步处理各类较复杂的课题 。 大量的工程 问题要求解热传导，扩散对流问题，以及结构动力学问题的暂态过程 。通常都 采用半解析法将空间坐标先行离散而得到常微分方程组来求解 半离散的热传导方程是 H n4-Ca十 ， 一 0 ． (1．1) 算 分 一 其 中C是热容量矩阵，对称正定，H 是热传导矩阵，通常为对称非负 × 阵}，是热源 向量， 暂 析 以上 皆为已知 n是待求向量函数，其初始值n。为已知。 在流体的热输运中，在热传导同时还有对流问题 此时热窖矩阵虽然仍为对称正定，但热 传导 H 不对称了。因此对于一般问题 ，H 不能作为对称看待 。 D 对于固体通过流体来强迫对流散热的问题 ，得到的热容矩阵c往往数值上相差较大 ，固 为气体的热容量比固体的热容量小很多 此时得到的方程虽然仍是(1．1)这种形式 ，但其齐次 2 鹊 蚧 方程的本征值分布大小相差也加大了。这对于通常的积分法带来许多问题 。 结构动力学对于暂态历程的需求是显然的结构动力学方程为 主+ G 4-K — r0)， 维 (1．2) 其 中 为对称正定质量阵，G为阻尼阵(对称非负)及陀螺力阵(反对称)之和，K为结构刚度 阵(对称非负) 当处理流体共同作用时，它可以是一般的矩阵，r是外力向量，皆为已知 。 是 待求位移 向量，其初态为已知 。， 0 为 已知 (1．3) 以上方程是用位移有限元法推导来的离散形式 如果采用一般变分原理来推导，则将得到2 个一阶常数微分方程组。 精细积分法宜于处理一阶方程 ，因此采用哈密顿体系常用的方法 ，见文 [6]，化成为 率文于 1993年 9月 2日收到 *国家 自然科学基金资助 维普资讯 计算结构力学及其应用 -0= Hv+ (1．4) 的形式 ，其中。是 2 维的待求向量，有 以下公式 一 户 户 ： Mic+ Gx／2， H — A D 日 ，1]，A：一M一G／2． B —GM—G／4一K， C：一GM一 ／2， D =M- ， f一 (0，r (1．5) 当 为对称正定 ，G为反对称，K为对称，且 ，=0时，即成为 自由陀螺系统，体系成为保守；此 时 占，D为对称矩阵，A=一 ，而H 成为哈密顿矩阵 。