Abaqus复合材料壳单元建模—姊妹篇1：常规建模step-by-step

采用商业有限元软件Abaqus进行复合材料结构建模时，一般有两种建模方法：常规建模方法和Composite layup快速建模方法，主要差异在创建属性、赋属性和指定铺层坐标系方面，常规建模方法和一般商业软件类似，将创建材料、创建属性、赋属性和指定铺层坐标系四个步骤分离，通用性较强，尤其是对于包含UMAT/VUMAT子程序开发的复合材料分析模型或者是三维实体单元显式动力学分析模型，仅支持该类建模方法；Composite layup快速建模方法将创建属性、赋属性和指定铺层坐标系三部分内容集成在一起，可一次性完成设置，效率较高。本文先从最基本的常规建模方法讲起。

一般对于大尺寸复合材料结构，跨厚度比例大，满足板壳理论的假设，采用壳单元就能获得高的求解精度。壳单元计算效率高，结合二维损伤起始判据判据（Hashin, Tsai-W, Maxe, Maxs等）可以预测结构的危险区域和危险程度，另外，Abaqus自身还内嵌了二维Hashin的渐进损伤分析模型，采用Hashin失效判据去判断损伤起始，损伤起始以后采用基于能量演化的连续退化准则对材料刚度进行退化。

Abaqus中常用的壳单元类型有S4、S4R、S8R等。以下介绍复合材料开孔板壳单元模型的建模步骤。

第1步：绘制几何

在Part模块下绘制几何，几何类型为3D-Deformable- Shell，草图如下：

绘制完草图后，退出草图，得到开孔板的几何模型，如下：

第2步：创建材料

与复合材料壳单元对应的是2D材料模型Lamina，将视图切换至Property模块，点击创建材料按钮，在跳出窗口中选择Mechanical→Elasticity→Elastic选项，在材料类型下拉框中选择Lamina，如下图所示。

表格中的6个数据分别为纵向（沿纤维方向）弹性模量、横向（垂直于纤维方向）弹性模量、面内泊松比以及三个方向的剪切模量。与其他商业有限元软件不同的是，即使是对于二维材料模型，仍然需要输入面外的剪切模量G13和G23，这两项数据是用于定义壳的横向剪切行为。

一个复合材料分析模型中可以包含多种材料模型，例如不同的铺层采用不同的材料。

第3步：创建截面属性

点击创建属性按钮，选择Shell→Composite类型，如下图所示。

点击“Continue”，进入属性设置界面。

该界面中首选要定义的是截面属性积分方式，有两种积分方式：During analysis和Before analysis，默认的是第一种。During analysis积分方式会在计算过程中，每个增量步重新计算截面属性，计算量较大，准确度高，当材料有明显非线性（塑性、损伤后刚度退化等）时，建议使用During analysis；Before analysis是在计算之前对截面属性进行预积分，分析过程中截面属性不再变化，计算量小。

该界面中第二个要定义的就是铺层信息表，共5列，分别是铺层材料、铺层绝对厚度、铺层角度、单层积分点数量、铺层名称。表格行数根据实际铺层数来定义。

当结构中不同区域有不同的铺层材料或者铺层顺序时，可以创建多个截面属性。

第4步：赋属性

点击赋属性按钮，按照主窗口下方提示区的提示信息去选择要赋予属性的区域，当有多个不同的铺层区域时，分别选择对应的截面属性逐个设置。

当几何结构被赋予属性时，颜色会由默认的灰色变成绿色，如下图所示。

第5步：指定铺层坐标系

复合材料分析模型与金属结构分析模型不同，由于材料具有各向异性，因此还需要额外指定铺层的坐标系。

点击创建坐标系按钮，有三种坐标系可以选择：直角坐标系、圆柱坐标系和球坐标系。可根据自己的结构特征选择合适的坐标系，针对本模型，选择直角坐标系，然后指定坐标原点、X轴方向以及XY平面即可。

坐标系创建完成以后，点击指定材料方位按钮，然后选择创建好的坐标系，并指定法向方向和附加的旋转角度，一般壳单元默认设置即可，在设定材料方位时，视图中几何模型上会显示相应的局部坐标系，可以将Additional Rotation附件旋转角度在0°、45°、90°及-45°之间不断切换来检查模型坐标系设置的对错与否。（注：在Abaqus复合材料分析中，一般有总体坐标系、铺层坐标系、材料坐标系以及断面坐标系等，将在后续文章中做深入讲解）。

第6步：创建装配

将视图切换到Assembly模块，点击创建实例按钮，进入以下界面：

创建实例有两种方式，一种是基于自身模型的Part创建，一种是基于其他模型创建，会将所选择的模型中的实例导入到当前模型下。

创建完实例之后，视图中会出现蓝色装配体。

第7步：划分网格

将视图切换到Mesh模块，进行网格划分，一般建议创建完几何之后，优先划分网格再去设置边界、载荷、接触等信息，对于复杂结构，在网格划分过程中需要对结构进行切分，切分会破坏原始构型，导致先前设置的边界、载荷、接触等信息丢失，造成工作反复，因此将划分网格这一步提前。

首先是设置网格密度，共有两个按钮，，左侧的为全局网格密度设置，右侧为局部网格密度设置，可以对特别的边线或者面设置不同的网格密度.

在设置网格时，如果出现以下提示，是因为创建实例时选择了Dependent类型，此类实例仅能在零件下划分网格或设置网格属性。

解决办法是将网格划分的对象由Assembly切换到Part，如下图所示。

在全局网格密度设置界面下设定近似全局网格尺寸为2mm，其他默认设置。

然后点击网格划分按钮去划分网格，划分网格以后视图如下：

可以通过网格控制按钮来设置单元形状和网格生成算法，

还可以通过单元类型设置按钮设置单元类型，默认的是四节点减缩积分单元S4R。

第8步：设置分析步

将视图切换至Step模块，点击按钮创建分析步，常用的分析步类型有Static,General和Dynamic,Expicit两种，本案例选择Static,General类型，创建一个静力分析步。

分析步总时间默认为1秒，增量设置可以按下图进行修改。

第9步：设置输出

Abaqus中的输出有两类，场输出和历史输出，场输出适合针对大量的对象进行低频次输出（如查案应力、应变及变形云图），历史输出适合针对少量的对象进行高频次输出（如监控系统能量或者绘制载荷位移曲线）

另外需要注意的是，复合材料结构在设置场输出变量时与各向同性材料也有不同，默认针对复合材料结构，场输出至输出上下表层，不会输出中间层，因此需要进行特殊设置。如下图所示，在Output at layered section points中应输入所有的截面点编号，比如壳单元有四层铺层，每一层有3个截面点，则一共需要输入12个截面点才能保证后处理中能查看所有铺层的场变量数据。

第10步：设置边界和载荷

视图切换到Load模块，设置边界条件和载荷。与常规各向同性材料模型一样，一端设置固支，一端给定拉伸位移1mm，此处略去。

需要注意的是，边界和载荷在整个建模过程中属于最不起眼、最容易被忽视的步骤，但同时也是对计算结果影响最大的因素之一，万不可小觑。

第11步：创建任务

切换到Job模块，创建任务，提交计算。与传统模型无异，此处略去。

第12步：查看结果

计算完成后，在Job模块下点击Results按钮，可直接打开结果文件并跳转至后处理界面。

查看应力时，默认显示的是第一层应力水平。

当需要查看每一层应力水平时，点击菜单Result→Section Points，如下图所示，切换至Plies选项，将会出现在截面属性中定义的所有铺层，每一层的输出位置又分为上中下三个位置，如下图所示。

第1层应力水平

第2层应力水平

【完结】

欢迎关注“复合材料力学”微信公众平台。同时，欢迎参加由复合材料力学公众平台与技术邻共同举办的Abaqus复合材料分析培训班，为期三天，白天上课，晚上练习指导、独家讲义、内容全面细致，由浅入深，理论与实际操作结合，带你一次掌握Abaqus复合材料分析。

培训大纲如下：

基础班

高级班

【独家讲义】

【联系人】

微信jm19961996，可扫码添加微信。