船舶碰撞动力学过程数值仿真研究.pdfVIP

第21卷 第 1 期 爆 炸 与 冲 击 Vol.21, No. 1 2001年 1 月 EXPLOSION AND SHOCK WAVES Jan. , 2001 : 10012 1455(2001) 012 00292 0 王自力, 顾永宁 (上海交通大学船舶与海洋工程学院, 上海 200030) : 在分析显式非线性有 元基本理论和关键技术的基础上, 探讨了船舶动力学过程的数值仿真方 法, 并对碰撞过程中船体的力学性能进行了分析, 获得并讨论了碰撞力、能量吸收和损伤变形的时序结果, 所 得结论具有一般性。碰撞仿真分析中, 被撞船的舷侧受撞区域和撞击船的首部区域作为可变形结构处理, 而 其余区域则视为刚体。被撞船周围的水划分成Euler 有 体积网格, 船体结构则划分成Lagrange 有 元网格, 并采用耦合技术将它们联系在一起。撞击船周围的水的影响采用附连水质量进行处理。 : 船舶碰撞; 数值仿真; 结构耐撞性; 流体2 结构相互作用 : U 1143 : A X 1 船舶碰撞是船体结构在很短时间( 大约零点几秒到一点几秒) 内在巨大冲击载荷作用下的一种复杂 的非线性动态响应过程, 它具有非常明显的动力特性, 而且碰撞区构件一般都要迅速超越弹性阶段而进 [1] 入塑性流动状态, 并可能出现撕裂、屈曲等各种形式的破坏或失效 。船体是由钢板焊接而成的水中结 构物, 庞大的船体结构、品种繁多的构件以及流体介质的影响决定了船舶碰撞问题的分析要大大复杂于 陆上结构碰撞。就船舶碰撞问题而言, 如果想通过建立一个精确的数学模型而使之得到完全解析, 那几 乎是不可能的, 试验研究固然是一种可靠的方法, 但船舶碰撞试验却是一种极昂贵的破坏性试验, 而且 难以实施。因此, 随着计算机硬件技术的飞速发展, 特别是非线性有 元技术的日益进步和成熟, 显式 有 元数值仿真技术在船舶碰撞问题的研究中逐渐受到重视。 我们在分析显式非线性有 元基本理论和关键技术的基础上, 着力探讨了船舶碰撞动力学过程的 数值仿真方法。相撞船舶周围的水划分成Euler 有 体积网格, 船体结构部分则划分成Lagrange 有 元 网格, 并用流体2结构耦合算法将Euler 网格和Lagrange 网格结合在一起, 解决水和船体之间的相互作用 问题。通过仿真计算, 可以获得碰撞力、能量吸收和损伤变形的时序结果。 2 211 在总体坐标系中, 碰撞运动方程可表示为 Ma + Cv + Kd = F ex ( 1) 式中: M 为质量矩阵, C 为阻尼矩阵, K 为刚度矩阵, a 为加速度向量, v 为速度向量, d 为位移向量, F ex 为包括碰撞力在内的外力向量。 若令F in= Cv + Kd , 并设F re= F ex- F in, 则碰撞方程也可以写成 Ma = F re (2) 如果采用集中质量, 即质量矩阵M 变成对角矩阵, 那么各个自由度的方程将是相互独立的, 即 X : 200020422 ; : 20002 08214 : 王自力( 19 4) ) , 男, 博士研究生, 副教授。 30 爆 炸 与 冲 击 第21卷 re M a = F i = 1, 2, ,