HyperMesh 基础建模操作

     在结构进行受力分析过程中，可以采用Abaqus、Ansys、HyperWork等软件按照结构模型进行建模分析，不同的软件所转化的有限元模型是不一致的，以金属的壳单元为例，在Abaqus内采用的单元为S4、S4R、S3、S3R等，所赋予的属性为SHELLSECTION；但OptiStruct内则采用CQUAD4、CTRIA3等单元，赋予PSHELL属性。

     以下的建模方法、操作都是基于澳泰尔（Altair）公司旗下的前处理软件HyperMesh, 相关卡片介绍以求解器OptiStruct为主。

     备注：软件准备为HyperMesh2024；几何模型准备为梁模型（Short_Spar.CATPart）

     以下为软件界面：

    本文的主要内容预览：

一、几何模型导入处理

1. 导入方法

     方法一：打开HyperMesh软件，点击“下拉菜单栏”中的“File”按钮弹出对话框，选择“Import”下的“Geometry”，在软件左侧会弹出以下窗口，点击选择文件按钮，即可选择所需的几何模型（支持CATIA、Step、UG、IGES等文件格式）。

将所需的几何模型选中，呈高亮状态，点击“Import”按钮即可导入图形界面，如下图所示：

    同时，切换到“Model”窗口，可以看到刚刚导入的几何模型所携带的信息已被分别归纳到了“Components”、“Parts”、“Titles”的组别下，另外，可点击空白部分创建“Assembly” （装配体）、“Materials”（材料）、“Property”（属性）、“Load Collectors”（载荷集）、“Load Step”（工况）等其他组别，对模型进行分组管理，如下图所示：

     方法二：通过“工具栏”中的快捷菜单“Import Geometry”打开导入窗口，如下图所示：

2. 几何处理

     对于导入的实体模型，需先将模型的几何元素进行简化，而后根据实际情况将其转化为有限元模型，通常情况下，采用规则的二维单元（最好是矩形）来表征零部件；对于局部细节，可以采用三维实体单元表征。

     以导入的梁模型为例，首先要确定合适的建模面，可以选择中面建模、外表面建模或内表面建模，而后可通过“页面菜单”下的“Geom”和“Tool”功能将建模面处理为适合划分网格的几何面，具体处理方式与操作流程如下所示：

     （1）选择在外表面建模，需先将外表面从实体上抽取出来，然后将其作为独立部件进行几何处理。获取实体表面方法：选取完整的外表面，将其复制到一个新的文档中，类似于文件的复制粘贴过程。

          ① 在标签区的“Components”组别下创建一个新的组件，命名为Surface（按需即可），如下图所示：

           ② 在页面菜单使用“Geom”项调出几何处理面板，点击“Surface edit”命令调出面处理界面，如下图所示：

          找到“Surfs”键并点击，而后在图形区选择完整的外表面，选中后呈高亮状态。如下图所示：

          再次点击“Surfs”，在弹出的功能界面中点击“duplicate”，最后选择“current comp”完成表面复制；

          ③ 打开“Components”组别下的“surface”组件，右击选择“Isolate”，抽取出来的外表面则会单独显示在图形区，如下图所示：

          点击侧边“工具栏”的“Display fixed points”命令，图形表面的硬点（几何点）则会显示出来，如下图所示：

     注：硬点point是几何点，即图形上的线端点，参与有限元分析计算的是（网格）节点node而不是硬点point。

     （2）表面单独显示后先处理边线，应减少斜边和邻近边线的数量，尽可能使保留下来的边线将外表面划分为规则的四边形（矩形最佳）；面内的边线尽量以绿色的共享边为主，这有利于有限元网格中的节点处于共节点状态，以便于后续在其内部进行网格划分。

           ① 点击下侧“工具栏”的“Geometry color mode”命令，选择“By 2D Topo”切换为2D模式，如下图所示：

注：图形上的线有3种，一是自由边（红色），自由度最高，仅与一个曲面关联，另一侧是断开的；二是共享边（绿色或黄色），同时与两个或三个曲面相关联，曲面在共享边处几何连续；三是压缩边（蓝色），位于两曲面的缝合处，可将两个面合为一个面。因此，两条邻近的自由边可以合并为一条共享边，一条共享边又可压至面内形成一条压缩边。

         ② 通过几何Geom处理面板下的“quick edit”功能（或快捷键F11）调出快速处理界面，如下图所示：

          ③ 左击“toggle edge”按键降低边线自由度，对应曲面的合并（红变绿）、缝合（绿变蓝）；右击则升高，对应曲面的分解（蓝变绿、绿变红），需注意只有当两边间距低于所设容差“tolerance”时，才能完成合并或分解动作。按照上述操作，将表面的绿色的共享边处理为蓝色的压缩边，如下图所示：

     （3）简化边线后，需要对表面进行切割处理，绘制出一个容纳规整网格的表面。处理方式：将大尺寸方向的面作为被切割面，找到三个点组成一个切割平面，切割平面与被切割面的交线即为切割位置，达到一种面垂面或其相交面上的线延伸下来的“一刀切”式的效果。切割平面的搭建思路为：测量被切割面的相交线的两个端点的距离，取其中点沿其所在任意面的法向平移，不共线的三个点所在的面即为切割平面。

          ① 通过快捷键F8调取节点编辑界面，可以通过硬点生成节点。确认“points”键为高亮状态，选择硬点，最后点击“create”将表面硬点转化为可搭建网格的节点，如下图所示：

注：通过硬点生成的节点，并不代表硬点不存在于表面，表面依旧存在硬点。

          ② 通过几何面板“Geom”下的“temp nodes”键块（或快捷键F4）调出样点界面，测量两点间距，并在其中点位置创建一个新点，如下图所示：

          ③ 通过工具Tool面板下的“translate”功能（或快捷键Shift+F4）调出平移界面，选择刚刚创建的节点，确认“nodes”键为高亮状态，在所需移动矢量的一个垂面上依次选择绿-蓝-红三个点，新建节点将沿此面的法向平移，如下图所示：

注：关于+/-平移方向的确定，依据右手定则通过三个点的选取顺序来确定平面的法向，也就是，法向为+。

          ④ 通过几何处理Gemo面板下的“surface edit”键块调出表面编辑界面，选择“trim with surface/plane”功能，确认“surfs”键为高亮状态，先选取被切割面（一个或多个），然后选择确定好的三个点（两个边线节点和一个新建节点）组成切割平面，点击“trim”完成切割。如下图所示：

注：红色框内代表被切割面，黑色框内代表切割工具面。

          ⑤ 重复上述操作，一般在截面变化的地方进行切割，使硬点尽可能分布在一条直线上且得到的切割面内不存在硬点，整个外表面划分为连续的（近似）矩形的面。如下图所示：

二、建立网格模型

1. 常用单元类型介绍

     本次主要介绍的常用单元类型属于OptiStruct与Nastran求解器下所用的单元。

 CROD单元

     该类单元用于定义杆单元，仅具有轴向和扭转刚度。杆单元是所有与几何结构相关联的单元中最简单的单元，在结构仅受拉伸、压缩，以及扭转作用时，结构可以采用杆单元表征。下图为杆单元在求解文件中的书写格式：

    EID：单元的ID号，E表示element；

    PID：单元的属性ID号，P表示property;注意:不同的单元所采用的属性不一致；

    G1、G2：组成单元的节点ID号，同时在杆单元内G1点指向G2表示杆单元单元坐标系中的X轴；

CBAR单元

     该类单元用于定义等截面的梁单元，CBAR单元可以承受拉伸、压缩、扭转、两垂直平面的弯曲和剪切。CBAR单元只能表征剪切中心在中轴线上的梁。下图为CBAR单元在求解文件中的书写格式：

     EID：单元的ID号；

     PID：单元的属性ID号；

     GA、GB：表示组成单元的节点ID，并且CBAR单元内GA点指向GB表示CBAR单元单元坐标系中的X轴；

     X1、X2、X3：GA点与其表示一个向量，用于定义CBAR单元单元坐标系中的Y轴；

    OFFT：表示GA点与GB点在哪个坐标系下偏移；

    PA、PB：表示释放GA点与GB点在单元坐标系下的某个自由度；

    W1A、W2A、W3A、W1B、W2B、W3B:对GA点与GB点的偏移量；

注意：对于CBAR单元的单元坐标系在有限元分析中非常重要，CBAR单元的坐标规定如下图所示，首先通过GA点与GB点确定单元坐标系的X轴，其次通过单元坐标系的X轴与向量组成的平面确定单元坐标系的Y轴，最后通过单元坐标系的X轴、Y轴以及右手定则确定单元坐标系的Z轴。

 CQUAD4单元

     等参四边形单元，常用于表征板单元、壳单元以及膜单元，可以表示平面内的拉压剪、面内弯曲。下图为CQUAD4单元在求解文件中的书写格式；

     EID：单元的ID号；

     PID：单元的属性ID号；

     Gi：单元节点ID号；

     THETA：定义材料坐标系的方向角度；

     MCID：材料坐标系ID；

     ZOFFS：单元节点定义的平面到壳单元参考平面的偏移量；

     Ti：单元节点的厚度；如果Ti存在值，PSHELL中定义的厚度会被覆盖，单元的厚度为Ti的平均值；

     在有限元分析中，分析完成后得到的应力、应变是单元坐标系下的应力应变。CQUAD4单元中可以通过Gi与右手定则确定单元坐标系的Z轴，单元坐标系中的X轴可以通过G1G3线与G2G4线之间的二分之一夹角处确定，X方向为G1指向G2,最后再通过右手定则确定Y轴。CQUAD4的单元坐标系定义图如下所示：

     对于复合材料的分析，应对CQUAD4单元定义材料坐标系，材料坐标系的定义可以通过两种方法定义：法一.通过定义Theta中的角度，去确认材料坐标系的X轴，材料坐标系的Z轴与单元坐标系中的一致；法二.通过定义参考坐标系去定义单元的材料坐标。材料坐标系的定义图如下所示：

2. 建立网格模型操作方法

     完成表面处理之后，可以通过2d面板下的“automesh”打开划分2d网格的命令，也可以通过快捷键F12打开。如下图所示：

     在网格划分面板下，选择需要划分网格的表面，设置单元的尺寸、单元的类型（四边形单元、三角形单元）、划分的单元属于哪一个component、单元为高阶单元还是低阶单元、两块表面上的单元是否共节点，然后点击“mesh”，如下图所示：

      点击完“mesh”按钮后，图像区会出现网格，如下图所示，可以通过“density”调整上图确定的单元尺寸，点击“mesh”重构单元；“mesh style”可以改变上图所规定单元的单元类型，点击“mesh”可以重构单元；“biasing”可以改变单元沿着边线的分布规律（线性分布、指数分布、正态分布）；“check”可以检查单元的质量。

     重复上诉操作，完成表面网格的划分，如下图所示：

三、网格模型赋予属性

1. 常用材料属性类型介绍

各向同性材料(MAT1)

     所谓的各向同性材料即材料内任意一点在任意方向上，材料的性能都是相同的。对于该材料而言，只需要确定弹性模量（E），剪切模量（G）、泊松比（μ）中的两个参数就可以描述材料，三者存在如下关系：

     在MAT1卡片下，还可以定义材料的密度（RHO）、热膨胀系数（A）、参考温度（TREF）、应力极限（ST、SC、SS）,这些参数对于重力分析、动力学分析、模型检查都很重要。下图为MAT1在求解文件中的书写格式：

     通常对于静力分析而言，各向同性的材料参数可以只给弹性模量（E），剪切模量（G）、泊松比（μ ）中的两个参数。可以把材料的密度给上，方便测量模型的质量。

二维各向异性材料（MAT8）

     各向异性材料：在不同的方向上，材料的性能是不同的。对于该类材料而言，需要确定1方向的弹性模量、2方向的弹性模量、12平面的剪切弹性模量、泊松比；

     对于MAT8卡片，同样可以定义材料的密度（RHO）、纵向横向的热膨胀系数（Ai）、纵向的拉压许用应力应变（XT、XC）、横向的拉压许用应力应变（YT、YC）、面内的剪切许用应力应变、蔡吴失效准则所需要的应力空间强度参数（F12）。下图为MAT8在求解文件中的书写格式：

     通常可以采用MAT8卡片去定义复合材料，在进行定义时，若采用蔡吴准则去判定材料是否失效需要定义F12,需要定义强度许用值，许用值是应力还是应变可以通过STRN字段下是1.0还是空白判断，1.0表示许用值是应变许用值、空白表示许用值是应力许用值。（通常对于复合材料采用应变许用值去判断）

     在进行复合材料分析时，若只有G12，没有G1Z、G2Z的话，横向剪切刚度为无限刚，可以令G12 = G1Z = G2Z，如果采用MAT8表征高阶非线性单元的话，G1Z、G2Z必须大于0。

     不采用蔡吴失效准则分析时，可以只填写E1、E2、NU12、G12、G1Z、G2Z,进行静力分析。

杆单元的属性PROD

     PROD主要是用于定义CROD单元的截面积、材料、扭转常数；下图为PROD在求解文件中的书写格式；

     PID：属性ID；

     MID：引用的材料ID；如果进行静力学和动力学分析，该处所使用的材料编号只能是MAT1的材料编号；如果进行热分析，该处所使用的材料编号只能是MAT4或者MAT5的材料编号；

     A：CROD单元的截面；

     J：扭转常数；

    C：扭转应力系数，扭转应力可以通过以下公式计算：

梁单元的属性PBAR

     PBAR主要是用于定义CBAR单元的截面积、惯性矩、扭转常数。下图为PBAR在求解文件中的书写格式：

     MID: 引用的材料ID；如果进行静力学和动力学分析，该处所使用的材料编号只能是MAT1的材料编号；如果进行热分析，该处所使用的材料编号只能是MAT4或者MAT5的材料编号；

      A：梁的截面积；

     I1、I2、I12：惯性矩；

     J：扭转常数；

     NSM：单位长度下的非结构质量；

     Ci,Di,Ei,Fi：应力恢复系数；

     K1、K2：剪切刚度系数；默认情况下，K1与K2无限大，而横向剪切刚度为Ki*A*G,所以横向剪切刚度无限大，但当I12=0时，可以忽略K1、K2；当A=0.0时，K1与K2必须为空。

膜、板、壳单元的属性PSHELL

     PSHELL属性可以用于定义CQUAD4、CQUAD8、CTRIA3、CTRIA6单元，可以定义膜单元、板单元、壳单元，单元的厚度。下图为PSHELL在求解文件中的书写格式：

     MID1：薄膜材料属性的ID；

     T：CQUAD4、CQUAD8、CTRIA3、CTRIA6单元默认的厚度；

     MID2：弯曲材料属性的ID；

    12I/T**3：惯性弯矩比率，I（实际壳的惯性弯矩）与T3/12（均质壳的惯性弯矩）的比值，默认为1.0；

    MID3：横向剪切材料属性的ID;

    TS/T：横向剪切厚度比，横向剪切厚度TS与膜厚T的比值，默认为0.8333；

    NSM：单位面积内的非结构质量；

    MID4：薄膜—弯曲材料属性的ID；

    计算结构质量所采用的密度主要是MID1的密度，如果MID1为空的情况下，采用MID2的密度进行计算。

    如果横向剪切MID3与膜弯曲耦合项MID4引用MAT2，则G33必须为0；若MID3引用MAT8，则G1Z和G2Z不能为0。

复合材料壳单元的属性PCOMP

      PCOMP属性用于定义表征复合材料的CQUAD4、CQUAD8、CTRIA3、CTRIA6单元，可以定义单元的偏置，层间剪切许用应力、失效准则、参考温度、层合板的定义方式、层合板内每一层的材料和材料坐标系。下图为PCOMP在求解文件中的书写格式：

     Z0：参考平面到底面的距离，如果Z0≠-1/2单元厚度，并且指定PARAM,NOCOMPS,-1，均匀单元应力（homogeneous element stresses）是不正确的，而单元力和单元应力是正确的。要想获得均匀单元应力，需要使用ZOFFS。不建议使用Z0去指定非对称铺层。非对称铺层存在拉伸和弯曲之间的耦合，这会使扰度增加，也会导致层合板刚度的下降，降低屈曲载荷和震动的频率；

     NSM：单位面积内的非结构质量；

     SB：层间剪切许用应力；如果指定了失效准则，需要填写；

     FT：失效判断准则；若不填写，则不会进行失效分析计算，Nastran支持蔡希尔、蔡吴、霍夫曼、最大应力应变、哈辛准则；

     TREF：参考温度；

     GE：阻尼系数；PCOMP的GE会覆盖单层材料的GE，PCOMP的GE需要进行等效才能使用；在瞬态分析中，没有写入PARAM，W4，GE会被忽略；阻尼系数可以通过临界阻尼比（C/C0）乘2.0；

     LAM：层合板定义方式；“SYS”只指定一半板的铺层和厚度；“MEM”、“BEND”、“SMEAR”、“SMCORE”所有的铺层都需要定义；“SMCORE”指定夹心和面板的属性，在进行刚度计算时，将面板的厚度进行等分放到夹心的上下表面进行，刚度计算时丢层的顺序是被忽略的；

    MIDi：不同层所用的材料ID；

    Ti：不同层的厚度；

    THETAi：单层的材料的纵向与材料坐标系X轴之间的夹角；

    SOUTi：应力应变输出请求；

   当PARAM,NOCOMPS等于-1时,将不会输出单层的应力应变，输出的应力应变为均匀应变；

2. 赋予属性的方法

     建成网格模型后，需将网格赋予材料属性（弹性模量、剪切模量、泊松比等）和几何信息（主要是厚度、表征方式）以进行静力分析。材料属性的赋予可分为直接和间接两种方式，即将创建好的Property直接赋予到网格单元上，或将其赋予到特定的Component上，该Component内的单元就拥有了属性。如果同时对网格赋予直接属性和间接属性，在进行分析时，首选直接属性。

          ① 在标签区空白部分建立一个材料“Materials”分组，将其命名为“Al”；

          ② 在标签区的下方，根据需求选择材料卡片，金属等各向同性材料选择MAT1，复材等各向异性材料选择MAT8；

         ③ MAT1卡片下填实际的材料参数，如铝合金的弹性模量E填入71000、泊松比NU填入0.33、密度RHO填入2.796e-09。如下图所示：

          ④ 建好材料卡片后，点击标签区空白部分建立一个属性Property，并将属性名称改为“Surface”，Component与Property的名字保持一致，方便管理模型。

          ⑤ 在模型树下方，按照需求选择Card Image（一维单元中CROD单元选择PROD卡片、CBAR单元选择PBAR卡片；表征金属的二维单元可以选择PSHELL卡片、表征复合材料的二维单元可以选择PCOMP卡片；表征三维单元选择PSOLID卡片），以表征金属的二维单元为例：选择PSHELL卡片，点击“Materials”右边的“<Unspecified>”呈现高亮后左击鼠标选择材料“Al”,点击“OK”；在“T：”后面填写几何实体的厚度8.0，完成属性定义；（几何实体的厚度可以通过快捷键F4调出测量点之间距离的面板，通过“two points”测量两个硬点的距离）。属性定义如下图所示：

          ⑥将属性与网格关联，采用直接属性的方法如下：点击属性设置按钮“Properties”，调出属性设置界面，选择“assign”。 如下图所示：

            “element”按钮呈高亮后，选择需要赋予直接属性的单元,单元呈现高亮状态，如下图所示（选择单元，法一.可以直接在图形界面通过shift+鼠标左键选择所需单元；法二.可以通过鼠标左击“element”,弹出对话框通过by collector选择单元所在的Component——Surface）

            选择完单元之后，“Type”可以选择2D，也可以选择All,“property”选择“Surface”,点击“assign”，完成直接属性的赋予。如下图所示：

           可以通过模型树下的按属性显示检查单元是否有属性，图形区有的网格为属性Surface的颜色，证明属性赋予完成。如下图所示：

           ⑦将网格赋予间接属性，方法如下：在建立完属性之后，直接点击网格所在的component——“Surface”,在模型树下方将属性property与组件component关联，点击“Property”右边的“<Unspecified>”呈现高亮后左击鼠标选择“Surface”，点击“ok”,完成属性property与组件component的关联。如下图所示：

注意：第⑥步为直接属性的赋予，第⑦步为间接属性的赋予，建议最好使用间接属性区对网格赋予属性，通过Component和Assembly可以更好的管理模型。如果需要去除直接属性，可以属性设置界面下的“unassign”去掉单元的直接属性，操作如下：先选择需要去掉直接属性的单元，直接点击“unassign”即可，如下所示：

     谢谢大家阅览，有不足之处，请大家多多包涵、多多指教！