砰击载荷作用下船体结构载荷预报及结构响应探究.doc

砰击载荷作用下船体结构载荷预报及结构响应研究 1、相关定义 1.1、自定义模型介绍 为了研究加筋板入水砰击载荷与刚体入水砰击载荷的区别以及加筋板在砰 击载荷下的响应规律,现自行设计一个平底加筋板结构。 该结构为 1.5 m ×1.0m ×0.3m 平底方盒结构,底板上布置两根 T 型材,型 材的间距为 0.5m ,跨距为 1.0m 。T 型材的腹板尺寸为 10mm×120mm,面板 尺寸为 12 mm×100mm,结构的横剖面如图 5-7 所示。结构材料为钢材,密度为 8000Kg/ m 3 ,弹性模量2.1 E5 Pa,屈服极限 235MPa ,泊松比为 0.3。 图 5-7 结构的横剖面视图 图 5-8 为加筋板底面俯视图,现在底面定义几个点,P1 为底板中心点,P2 为加强筋中间点与底板相交点,P3 为右侧板的中心点,P4 点未在图上标出,它 的位置在 P2 点的上方,为 T 型材面板一侧的中心点。上述四点作为提取数据和 天津大学硕士学位论文 第五章 平底加筋板结构入水砰击的数值模拟 48 数据比较的参考点。 图 5-8 加筋板底面俯视图 根据上述数据建立有限元模型。如图 5-9 所示:流体区域的 XYZ 方向尺寸 为 6m ×10m ×6m ,下部绿色区域为水,高度 5m ,上部蓝色区域为空气,高 度 1m ,空气的一部分被隐去以便看到处于其中的结构(黄色区域)。水的密度 为 1000Kg/ m 3 ,体积模量为2.2E 9 。空气的密度为 1.2887Kg/ m 3 ,伽马值为 1.4 。 图 5-9 仿真模型的三维视图 图 5-10 平底结构的有限元模型 图 5-10 为平底结构的有限元模型,顶盖被隐去以便能看清内部的 T 型材。 结构表面被划分为若干个四边形面单元,用弹塑性材料填充。平底结构的外表面 被定义为耦合面,且耦合面的法线都指向外部。 天津大学硕士学位论文 第五章 平底加筋板结构入水砰击的数值模拟 49 1.2、摇臂与船体尾部结构接触定义 在摇臂与船尾结构发生碰撞的过程中,两者之间的接触区域、接触时间和相 互作用力是随时间变化的,并且是未知的。但是,在本文碰撞接触分析过程中, 摇臂与船尾结构的接触状态是接触研究的重要内容。因此,合理地定义接触类型 及计算方法是此碰撞分析中不可缺少的关键环节。 LS-DYNA 程序处理接触-碰撞界面主要采用三种不同的算法,即节点约束 法、对称罚函数法、分配参数法。其中,第一种方法现仅用于固连界面,第三种 仅用于滑动界面,而对称罚函数法具有对称性、动量守恒准确,不需要碰撞和释 放条件,很少引起 Hourglass 效应,因此,其已发展为一种非常用的接触界面算 法,在数值计算中被广泛应用。在计算过程中,每一时步先检查各从节点是否穿 透主表面,如果没有穿透,则对该节点不作任何处理。如果穿透,则在该从节点 与被穿透主表面之间引入一个较大的界面接触力,其大小与穿透深度、主片刚度 成正比,称为罚函数值。 另外,其主要有三种基本接触类型:单面接触,点-面接触,面-面接触。面 -面接触类型是最常用的,并且常用于任意形状且存在较大接触面积的物体接触, 另外,自动接触是使用最广泛的接触类型选项之一,它能通过自动接触算法自动 确定壳单元的接触表面方向。一般来说,平直的或凹面的、粗网格表面定义为 target 表面,凸面的、细网格表面定义为 contact 表面。 摇臂回转运动过程中,其前端首先与船体尾部的小垫板接触,继而与大垫板 接触,由于碰撞过程中不确定因素较多,其可能也与尾部甲板大面积接触,因此 , 本文采用自动面面接触(ASTS)计算,并定义摇臂表面为 target 表面,而船体 小垫板等结构为 contact 表面。 大部分船舶碰撞问题研究中,接触问题的摩擦现象被视为最简单的库仑摩 擦,而本文研究中,由于摇臂与驳船尾部的碰撞情况较特殊,属于绕定点转动碰 撞,碰撞过程中几乎不存在相对滑动,因此,本文忽略了摩擦系数的影响。 值得注意的是,定义接触之前需要先生成 PART,如下图所示: 34 天津大学硕士学位论文 第三章 大型下水驳船摇臂与船体尾部的碰撞过程分析 图 3-1 0 PART 生成、修改控制 图 3- 11 模型所有 PARTS 信息 生成 PART 之前,必须要进行网格划分,因为 PART 就是由相同单元类型, 实常数和材料号的组成单元集。 GUI: Main Menu: PreprocessorLS-DYNA optionsInput File/PARTSPARTS DATA,得到 3-10 图,可以直接生成、更新,删除 PARTS。另外,选择最后一 个选项,就能列出了模型 part 的全部信息,如 3-11 图所示,从这里也可以检测 材料属性的正确性。 定义接触类型和接