Patran船舶结构有限元建模的

何勇

目前，市面上主流的通用结构有限元分析前处理器种类众多，包括 Patran、Ansys、Abaqus、HyperMesh和 Femap等。这五款软件在国内船舶行业均有应用，其中 Patran的用户数量很多。本文结合笔者多年运用 Patran进行船舶结构有限元建模积累的经验，为广大同行介绍一些建模方面的技巧，以期对读者有一些帮助。

一、网格划分与编辑技巧

1.根据现有的节点撒种子

对于具有曲线边界的面，例如首尾区域垂向加强的横舱壁，横舱壁与外板交线上的网格种子并不是等间距的，此时可先在横舱壁与甲板的交线上生成等间距的节点，然后将这些节点投射到横舱壁与外板的交线上，最后用Create＼Mesh Seed＼Tabular方法（Coordinate Type选Node and Point）根据投射得到的节点生成网格种子，这样得到的网格种子可以保证横向对齐。

2.径向扫略

由低阶单元变换得到高阶单元时，Sweep＼Element＼ Extrude是一种经常使用的方法，此方法为线性扫略。但径向扫略的作用也不容忽视，图 1所示的立柱环筋的创建就是一个典型的例子。

径向扫略通过 Sweep＼Element＼Radial Cyl.方法实现。扫略时，Axis填写扫略平面的法向，Radial Distance填写扫略距离，为正时背离原点扫略，为负时指向原点扫略。

3.跨越边界的单元的修改

通过扫略创建甲板边板壳单元时，在非平行中体区域，会出现最外侧的一行单元跨越甲板边线的情况，此时需移动节点将其拉回到甲板边线上。有些用户习惯使用 Modify＼Node＼Move方法实现这个目标，但这种方法要求输入节点要移动到的新位置，而这个位置的坐标通常是未知的，因此需先通过求交点的方法找到这个位置。

更快捷的做法是使用 Modify＼Node＼Project方法，Project onto选 Surface，Direction选 Define Vector，即将节点沿指定的方向投射到指定的面上。为此，先用 Create＼Surface＼Extrude方法将甲板边线沿垂向扫略成一个面，然后将节点沿横向投射到这个面上，这样可一次性完成所有节点的移动。

4.单元切割

手工修改网格时，经常需要对单元进行切割，可通过Modify＼Tria|Quad＼Split方法进行，该方法提供了多种切割模式。例如，对于四边形单元，有图 2（a）所示的 6种模式。可以看到，这 6种切割模式还是不够丰富的，作为补充，Utilities>FEM-Elements>Break Elements菜单命令提供了更全面的切割模式，如图 2（b）所示，用户可按需选用。

（a）E lements模块切割模式

（b）Ut ilities切割模式

二、单元赋属性技巧

在 Patran中，单元赋属性是一项很繁重的工作，因为与其他前处理器不同，Patran中没有当前材料、板厚和梁截面的概念，因此所有单元必须指定属性。此外，船舶结构包含大量的加强筋，用梁单元模拟这些加强筋时，赋属性需指定方向矢量和偏置矢量，比较耗时。因此掌握一些单元赋属性的技巧非常有必要。

1.自动赋属性

使用 Create＼Element＼Edit方法手工创建单元时，可单击操作面板上的 Select Existing Prop按钮，在弹出的 Property Sets面板中选择适当的单元属性集，将新创建的单元添加到所选择属性集的作用区域中，从而实现自动赋属性的功能。

2.在单元坐标系中填写梁单元的偏置矢量

默认情况下，梁单元的偏置矢量以梁端节点的分析坐标系（默认情况下为全局坐标系）为参考进行填写。这种填法的不便之处在于用户需要随时判断偏置矢量的非零分量及其正负。例如，甲板纵骨的偏置矢量是 Z分量非零且为负；左舷舷侧纵骨和右舷舷侧纵骨的偏置矢量虽然都是 Y分量非零，但前者为负，而后者为正。当加强筋腹板不与某个主坐标平面平行时，偏置矢量的三个分量可能均为非零，填写尤其麻烦。那么是否有一种更高效且不易出错的填写方法，使偏置矢量不随加强筋的类型和方向而变呢？答案是肯定的，那就是在单元坐标系中填写偏置矢量。

图 3为底边舱斜板（与水平面成 45°）纵骨（左舷，规格为 L125x80x8）的偏置矢量在全局坐标系和单元坐标系中的填写方法，两种方法产生的实际效果都如图4所示，可以看出在单元坐标系中填写更为方便。

（a）全局坐标系中

（b）单元坐标系中

将参考坐标系切换为单元坐标系的方法为向右拖动 Input Properties面板中的水平滚动条，将Offset @ Node1/2文本框最右侧的下拉列表选项切换为 Element（默认为Analysis）。切换后便可在单元坐标系中填写偏置矢量。为方便读者参考，图 5列出了三种常用加强筋截面（单元类型为 CBEAM）的单元坐标系。

3.使用局部坐标系作为梁单元的方向矢量

默认情况下，梁单元的方向矢量以第一个节点的分析坐标系（默认为全局坐标系）为参考进行填写。在某些情况下，这种填法非常麻烦，图 6所示的起重机立柱垂向加强筋就是一个典型的例子。

图 6中各垂向加强筋每隔 45°均匀分布。很显然，在全局坐标系中，各加强筋的方向矢量各不相同，这就意味着要为每根加强筋分别定义一个属性集。事实上，Patran支持在用户自己定义的局部坐标系中填写方向矢量，方法为在填写 Bar Orientation时，先输入方向矢量，然后空一格再输入“Coord n”（n为局部坐标系的编号），这样就表示方向矢量以编号为 n的局部坐标系为参考。对于上例，可以以立柱的中心为原点创建一个柱坐标系（假设为Coord 2），这样加强筋的方向矢量均沿该柱坐标系的径向并指向原点，即方向矢量可填为“<-1 0 0> Coord 2”。

三、其他常用技巧

1.组的导出与导入

计算书送审或者在不同版本的 Patran之间传递有限元模型时，通常会用到 bdf文件，但是 bdf文件中并无模型分组的信息，当导入 bdf文件重新构建模型时，读入的节点和单元将位于一个组中。如果模型比较复杂，不分组时模型操作是非常不便的，因此有必要重新生成原数据库文件中的组，步骤如下。

（1）执行 Utilities>Group>Groups To Session File菜单命令，将创建组的命令写入一个对话文件。

（2）导入 bdf文件后，执行 File>Session>Play菜单命令播放上一步生成的对话文件。

2.Region工具的使用

有时我们希望把某些节点或单元归到一起，供后续的某些操作使用，但是为避免分组过多，且不希望把节点或单元放到一个组里，此时可以借助 Regions工具。创建 Region的菜单命令为 Tools>Regions，创建好的 Region可以作为载荷 /边界条件和单元属性集的作用区域，使用方法为在 Select Application Region面板中将Select切换为Region。

3.PCL表达式的应用

有时我们在 Patran中输入的数据可能需要通过计算才能得到，此时大多数用户的做法是先用计算器等外部工具算好再填入 Patran，利用 PCL表达式可以省略这个中间过程。以图a）偏置矢量的填写为例，用户可以输入“<0 `121*cosd（45）``-121*sind（45）`>”，单击 Apply按钮后 Patran将自动将表达式替换为其结果。用户需特别注意表达式两侧的撇号，此符号位于键盘左上角，与“～”占据同一个键。