船舶船闸钢闸门碰撞的非线性数值仿真.pdfVIP

船舶-船闸钢闸门碰撞的非线性数值仿真 1 1 2 王政平 ，施泉 ，李林友 1 河海大学交通学院，南京（210098) 2 余姚市水利电力建筑勘测设计院，浙江余姚（315400 ） E-mail：sunwzhp@ 摘 要：研究船闸闸门的撞击性能，提高船闸的安全运作能力对内河航运业的稳定发展有 着重要的意义。本文基于显式有限单元法，利用Ls-dyna 计算平台，对船闸钢闸门遭遇船舶 撞击事件进行了数值仿真计算，并分析了闸门的撞击响应，克服了传统简单动力分析模式 的弊端，为船闸钢结构闸门的设计、制造和安装提供了参考依据。有限单元法数值仿真技 术为钢结构闸门动力研究提供了一个崭新的方法。 关键词：船舶，碰撞，钢闸门，有限元 中图分类号:U641.3+31 文献标识码:A 1. 前言 随着经济贸易的快速发展，水路运输讯速增长（江苏地区内河运量近年年均增长率达 8% ），船舶数量和航行速度都在不断增加[1] 。同时，由于通航船舶数量大，速度快，舵效 低，闸室壁、闸门遭遇船舶碰撞的事件越来越多，威胁着船闸的安全运作，给正常的水上航 运和工农业的输水灌溉带来安全隐患。经查，目前江苏省许多重要内河航运干线的闸门因碰 撞受损严重。这也是我国水运发展中的新问题。因此，研究船闸闸门的撞击性能，对高品质 闸门设计和船闸的安全运作有着十分重要的意义。 图1 船舶与闸门碰撞的实况照片 然而，船舶的撞击过程是复杂的非线性动态过程，与船舶排水量、撞击速度及撞击角度、 [2] 船舶结构及受撞物的变形特性等因素有关 。目前对船舶钢结构撞击的研究尚少，而且相关 研究也主要是通过简单冲量估算或物模实验来测算撞击力等，并不涉及结构及其整体的受撞 响应，现行规范也没有明确的相关规定。冲量估算无法表达整个撞击过程，且精度低；实验 测算精度很大程度上受初始条件的控制和测量设备精度限制，而且费用昂贵。 本文基于显式有限单元法，利用 Ls-dyna 计算平台[3][4]，对船舶撞击船闸钢闸门进行了 数值仿真计算，并分析了钢闸门的撞击响应。 2. 碰撞仿真有限元仿真模型的建立 江苏某船闸闸门宽 13m，闸门板采用对称三角钢架结构。闸门板高 6.4m，宽 6.54m， 单扇闸门重 24t 。闸门完全打开时，闸门主构架内侧上下两端及转轴固定。根据历年通航统 计资料，假设闸门在开闸通航时，侧壁遭遇重 1000t 的船舶的撞击，撞击速度为2m/s ，方向 - 1 - 0 与船闸轴线夹角 20 ，初始撞击点在闸门板的中央位置。 图2 船舶与闸门碰撞的平面示意图 2.1 几何模型 本文旨在研究钢闸门结构受船舶撞碰之响应，所以钢闸门构件及其连接方式（焊接、高 强螺栓铆接及与闸墩的转动）均按际情况建模，而船舶只考虑船体外壳刚度，对船舶内部的 设施和设备等进行了简化，并通过质量配重、调整，保持了船舶原有的惯性。 2.2 网格模型 闸门板、角钢、槽钢及船壳均采用四节点或三节点 Belytschko-Tsay 壳单元，面内单点 [3] 积分 ，单元的属性反映实际几何模型。 结构大位移大变形时出现网格畸形、材料破碎，使得简单连接区域变成多连接，甚至网 格重叠造成负质量，这将导致时间步长缩小、求解精度降低，甚至出错。本文采用 h-adaptive 方法进行自适应网格划分，使变化剧烈区的网格在迭代过程不断调节、细化，做到网格点分 布与物理解的耦合，