薄壁结构动力计算中一维有限元方法的应用-长江科学院院报.pdf

第2 3 卷第5 期 长　江　科　学　院　院　报 Vol.23 No.5 20 0 6 年10 月 Journal of Yangtze River Scientific Research Institute Oct.2 0 0 6 　　文章编号:1001 5485(2006)05 0041 04 薄壁结构动力计算中一维有限元方法的应用 1 1 2 汤建宏 , 聂孟喜, 梁应辰 (1.清华大学 水利水 工程系, 北京　100084;2.交通部三峡办公室, 北京　100084) 摘要:针对目前薄壁结构领域存在的建模与计算耗时太大的缺点, 将一维有限元的方法引入计算程序的编制, 可以 大幅度地降低单元与自由度的数目, 迅速得到薄壁结构的各项几何参数, 从而为后面的弯曲分析、扭转分析、乃至 动力学时程响应计算打下基础。 数值算例的检验结果表明, 这套程序计算所得到的结果与大型有限元软件 NAS- TRAN 计算的结果相一致, 但耗时却大大缩短, 这一特点可以便于结构设计师在结构的初始设计时对抗震动位移 迅速作出判断, 从而可以大大提高结构设计的效率。 关　键　词:薄壁结构;一维有限元方法;动力计算 中图分类号:TV37　　　文献标识码:A 　　由于薄壁结构高强、省材等优点, 这些年来在工 壁结构离散成等断面的单元, 得到沿高程方向一致 程领域内普遍产生了各种各样的高耸薄壁结构, 从 或线性变化的计算模型, 然后根据薄壁结构力学基 而也有越来越多的学者开始着眼于这种结构的力学 本理论建立单元刚度矩阵, 按有限元方法的一般步 [1, 2] 仿真分析 。 骤形成结构的总刚度矩阵, 并进行约束处理后即可 从原则上说, 一切结构都可以离散成有限元模 进行求解。 型, 用常规的结构分析有限元程序来计算。但在实 一维有限元方法有其应用范围的局限性, 即它 际操作中, 对于结构工程领域中的高层建筑、桥梁、 只能用于对薄壁结构进行宏观的扭转分析与动力计 船体等一些形状复杂的薄壁结构, 建立模型和准备 算等, 因为这种方法的一个重要假定是“刚周边假 输入数据需要大量的人力和时间, 而且由于截面复 定”, 这是由薄壁结构力学理论鼻祖———符拉索夫提 杂, 需要离散的单元数目较多, 从而对计算机的容量 出的, 即认为“薄壁杆件在外载荷作用下变形时, 断 和速度都有较高的要求, 计算费用也较大。所以目 面的形状不发生变化”, 因此这个假定被称为符拉索 [5] 前这一领域的发展受到了较大的制约, 常规的有限 夫假定。其中“薄壁杆件”的定义一般如下 :壁厚 元方法建模与计算的耗时让人们感到不甚麻烦, 故 与截面宽度之比t /b1/15, 壁厚与杆件长度之比 而需要发展一种高效的有限元方法。 t /L1/100, 此即为薄壁杆件。 值得一提的是, 目前已有部分学者对此领域进 在这个假定之下, 同一截面内各个质点的横、 [3, 4] 行研究 , 但大多并未形成完善的计算系统(比如 纵、垂向的平动位移以及速度、加速度、扭转角度都 前处理几何参数的计算等), 且没有涉及到动力学领 是一致的, 也就是一般有限元方法建模时有时会用 域, 而本文则在现有的研究基础之上对计算系统进 到的“刚化固结处理”。 行了完善与拓展。 有了这一假定, 建模时在截面内就可以不需细 分单元, 仅需沿高程方向划分即可, 从而可以大量减 少单元数目, 提高计算效率, 这也是一维有限