单元与网格选择技术在有限元分析中应用.docVIP

单元与网格选择技术在有限元分析中应用 　　本文结合SolidWorks Simulation软件中单元的技术特点与网格划分的规则，阐述了结构有限元分析中常见的单元类型和力学模型，并在此基础上讨论了仿真分析中单元的选取技术和网格的离散规则，为得到准确的有限元仿真结果提供了技术参考。 　　一、引言 　　有限元法解决工程实际问题一般包括前处理、计算和后处理三部分。其中前处理部分的工作量占整个分析过程工作量的40%～50%，后处理过程占45%～55%，而计算只占5%左右。计算和后处理过程是建立在前处理过程的基础上的，因此有限元的前处理过程相当重要。在有限元的前处理过程中，进行的主要工作是划分网格。网格划分的好坏直接关系到计算与分析的结果，是有限元分析的关键。SolidWorksSimulation用四面体实体单元、三角形壳单元分别划分实体与面几何信息，只有这两种网格单元才能对几乎任何几何实体或面进行可靠的网格划分。 　　二、SolidWorks Simulation网格特点 　　SolidWorksSimulation是一款基于三维设计软件SolidWorks界面的有限元仿真分析软件。根据模型几何特点并考虑设计工程师对操作速度的要求，Simulation中对网格分为五种基本单元类型：实体单元、壳单元、梁（杆）单元和2D单元。其他如接触单元、刚性杆单元等存在于软件内核中，使用者不需要直接操作此类单元。 　　1.实体单元 　　Simulation中的3D实体单元分为两种类型：一阶四面体单元也称为【草稿品质】，二阶四面体单元称为【高品质】。 　　一阶四面体单元在实体内沿着面和边缘模拟一阶线性位移场，即一阶单元。每一个一阶四面体单元共有四个节点，分别对应四面体的四个角点。每个节点有三个自由度，节点位移可完全由三个位移分量来表示。一阶单元的边是直线，面是平面。在单元加载变形后，这些边和面仍保持直线和平面，不能很好模拟曲面型几何模型。一阶单元变形如图1所示。 　　二阶四面体单元模拟了二阶位移场以及相应的一阶应力场，每个二阶四面体单元有十个节点，并且每个节点有3个自由度。当单元因加载而变形时，如果单元需要模拟曲线型几何模型，二阶单元的边和面就可以是曲线型形状。二阶单元变形如图2所示。 　　2.壳单元 　　Simulation中的壳单元分为一阶三角形壳单元和二阶三角形壳单元。 　　一阶三角形壳单元沿其面和边模拟线性位移场和常数应变、应力。当单元变形时，一阶壳单元的直边仍保持为直线，如图3所示。每个一阶壳单元上有3个节点，每个节点上有6个自由度，它的位移可完全由三个平移分量和三个转动分量描述。 　　二阶三角形壳单元模拟二阶位移场和一阶线性应力场，每个二阶壳单元有六个节点，三个角点和三个中间点。在划分网格过程中，当单元因加载而变形时，如果需要模拟曲线型几何模型，二阶壳单元的边和面就可以模拟曲线形状，如图4所示。 　　3.梁（杆）单元 　　相对于一阶实体和壳单元，两个节点的梁单元把两个面挠度模拟为三次函数并把轴向平移和扭转模拟为线性，两节点梁单元在每个端点处都有六个自由度：三个平移自由度和三个旋转自由度。 　　上述的一些网格信息在商用的分析软件中理论、技术都十分成熟，在现实的工程分析应用时往往需要选择合适的网格或者使用混合网格离散模型。所以要想运用好网格还需要对网格有深度的了解。 　　三、网格划分与单元选择规则 　　Simulation提供“标准网格”、“基于曲率的网格”两种网格处理器。基于曲率的网格处理器可以在网格化之前检查实体之间的干涉状态。如果检查到干涉，则网格化会中止。同时基于曲率的网格始终兼容钣金实体和表面实体的接触或部分接触边缘。基于曲率的网格处理器还支持相触实体面之间的兼容网格化。当分析模型曲面较多、不规则，建议尝试“基于曲率的网格”。可以提高生成网格的成功率。在【选项】→【默认选项】可以设置默认网格，如图5所示。 　　当模型在一个尺寸很小另外两个尺寸很大的情况下，为保证计算精度需要在小的尺寸方向至少有两层网格，所以实体划分网格会很不经济。针对有此特点的模型，我们需要使用壳单元。壳网格属于二维网格信息，即面网格。将实体模型转换为曲面模型，模型中面是首选的面。定义壳的厚度信息默认也是以中面两侧加厚处理的，当模型比较复杂时所有实体的中面不能得到也属正常，往往选取实体的内表面、外表面作为壳面对结果影响不大。软件还提供【偏移】设置，可以保证壳面正确的模拟实际模型，如图6所示。 　　壳单元分为【薄壳单元】、【厚壳单元】。薄壳单元：技术假定横截面垂直于中面且保持平直，同时保证变形之后仍然垂直于中面。薄壳单元忽略厚度方向的剪切变形及应力。当宽比厚值大于20时应定义为薄壳单元。厚壳单元：技术假定横截面垂直于中面并在变