1、5应用壳单元 应用壳单元可以模拟结构，该结构一个方向的尺度（厚度）远小于其它方向的尺 度，并忽略沿厚度方向的应力。例如，压力容器结构的壁厚小于典型整体结构尺寸的 1/10，一般就可以用壳单元进行模拟。以下尺寸可以作为典型整体结构的尺寸： 支撑点之间的距离。 加强件之间的距离或截面厚度有很大变化部分之间的距离。 曲率半径。 所关注的最高阶振动模态的波长。 ABAQUS 壳单元假设垂直于壳面的横截面保持为平面。不要误解为在壳单元中也 要求厚度必须小于单元尺寸的1/10，高度精细的网格可能包含厚度尺寸大于平面内尺 寸的壳单元（尽管一般不推荐这样做），实体单元可能更适合这种情况。 5.1单元几何尺寸
2、在 ABAQUS中具有两种壳单元：常规的壳单元和基于连续体的壳单元。通过定 义单元的平面尺寸、表面法向和初始曲率，常规的壳单元对参考面进行离散。但是， 常规壳单元的节点不能定义壳的厚度；通过截面性质定义壳的厚度。另一方面，基于 连续体的壳单元类似于三维实体单元，它们对整个三维物体进行离散和建立数学描述， 其动力学和本构行为是类似于常规壳单元的。对于模拟接触问题，基于连续体的壳单 元与常规的壳单元相比更加精确，因为它可以在双面接触中考虑厚度的变化。然而， 对于薄壳问题，常规的壳单元提供更优良的性能。 在这本手册中，仅讨论常规的壳单元。因而，我们将常规的壳单元简单称为“壳 单元”。关于基于连续体的
3、壳单元的更多信息，请参阅ABAQUS分析用户手册的第 15.6.1 节“Shell elements：overview” 。 5.1.1壳体厚度和截面点（ section points） 需要用壳体的厚度来描述壳体的横截面，必须对它进行定义。除了定义壳体厚度 5-1 之外，无论是在分析过程中或者是在分析开始时，都可以选择横截面的刚度。 如果你选择在分析过程中计算刚度，ABAQUS 采用数值积分法沿厚度方向的每一 个截面点（ sectionpoints）（积分点）独立地计算应力和应变值，这样就允许了非线性 的材料行为。例如，弹塑性材料的壳在内部截面点还保持弹性时，其外部截面点可能 已经达到了屈服
4、。在S4R（4 节点、减缩积分）单元中唯一的积分点的位置和沿壳厚 度上截面点的分布如图5-1 所示。 图 5-1在数值积分壳中截面点的分布 当在分析过程中积分单元特性时，可指定壳厚度方向的截面点数目为任意奇数。 对性质均匀的壳单元，ABAQUS 默认在厚度方向上取5 个截面点， 对于大多数非线性 设计问题这是足够了。但是，对于一些复杂的模拟必须采用更多的截面点，尤其是当 预测会出现反向的塑性弯曲时（在这种情况下一般采用9个截面点是足够了） 。对于线 性问题， 3 个截面点已经提供了沿厚度方向的精确积分。当然，对于线弹性材料壳， 选择在分析开始时计算材料刚度更为有效。 如果选择仅在模拟开始时计算
5、横截面刚度，材料行为必须是线弹性的。在这种情 况下，所有的计算都是以整个横截面上的合力和合力矩的形式进行。如果需要输出应 力或应变，在壳底面、中面和顶面，ABAQUS 提供了默认的输出值。 5.1.2壳法线和壳面 壳单元的连接方式定义了它的正法线方向，如图5-2所示。 5-2 图 5-2壳的正法线 对于轴对称壳单元，从节点1 前进到节点2 的方向经逆时针旋转90 定义其正法 线方向。对于三维壳单元，根据出现在单元定义中的节点顺序，按右手法则围绕节点 前进给出其正法线方向。 壳体的顶表面是在正法线方向的表面，对于接触定义称其为SPOS 面；而底表面 是在沿着法线负方向的表面，对于接触定义称其为S
6、NEG 面。在相邻壳单元中的法线 必须是一致的。 正法线方向定义了基于单元的压力载荷（element-basedpressure load）应用的约定， 和随着壳厚度变化的量值的输出。施加于壳体单元上的正向压力载荷产生了作用在正 法线方向的载荷。 （基于单元的压力载荷的约定，对于壳单元是相反于对实体单元的约 定；基于表面的压力载荷的约定（surface-based pressure load），对于壳单元是相同于 对实体单元的约定。关于在基于单元的和基于表面的分布载荷之间的更多区别，请参 阅 ABAQUS 分析用户手册的第19.4.2 节“ Concentrated and distribut
7、ed loads” 。 ） 5.1.3壳的初始曲率 在 ABAQUS 中壳（除了单元类型S3/S3R、S3RS、S4R、S4RS、S4RSW 和 STRI3 之外）的公式是描述了真实的曲壳单元；真实的曲壳单元需要特别关注对初始壳面曲 5-3 率的精确计算。 在每一个壳单元的节点处，ABAQUS 自动地计算表面法线来估算壳的 初始曲率。应用相当精确的算法确定每一节点处的表面法线，在ABAQUS分析用户 手册 15.6.3 节“Defining the initial geometry of conventional shell elements”中详细地讨 论这种算法。 若采用图5-3 所示的粗
8、网格，在连接邻近单元的同一个节点上，ABAQUS 可能会 得到多个独立的表面法线。在单一节点上有多个法线的物理意义是在享用共同节点的 单元之间有一条折线。而你可能打算模拟这样一个结构，更希望它是一个拥有平滑曲 面的壳体； ABAQUS 将尝试在这种节点处创建一个平均的法线从而使得壳面平滑。 图 5-3网格细划对节点处表面法线的影响 所采用的基本平滑算法如下：如果与同一节点连接的所有壳单元在该节点处的法 线相互之间的夹角在20 以内，则这些法线将被平均化。平均法线将用作为所有与该 节点相连的单元在该节点的法线。如果ABAQUS未能光滑壳面，在数据文件中 （.dat）将发出一个警告信息。 有两种方
9、法可以改变默认的算法。为了在曲壳中引入折线或者用粗网格模拟曲壳， 或者是在节点坐标后面给出n2的分量，作为第4、第 5 和第 6 个数据值（这种方法需 要在文本编辑器中人工编辑由ABAQUS/CAE创建的输入文件） ；或者应用 *NORMAL 选项，直接规定法线方向（应用 ABAQUS/CAE的 Keywords Editor（关键词编辑器） 5-4 可以加入这个选项，见第6.1.2 节，“Cross-sectionorientation”）。如果应用两种方法， 后者优先。关于进一步详细的信息，请查阅ABAQUS分析用户手册的第15.6.3 节 “Defining the initial ge
10、ometry of conventional shell elements” 。 5.1.4参考面的偏置（ referance surface offset） 通过壳单元的节点和法线的定义来定义壳的参考面。当用壳单元建模时，典型的 参考面是重合于壳体的中面。然而在很多情况下，提出将参考面定义为中面的偏置更 为方便。例如，由CAD 软件包创建的面一般代表的或者是壳体的顶面或者是底面。 在这种情况下，定义参考面并与由CAD 创建的面一致是更容易的，因此，该参考面 偏置于壳体的中面。 对于接触问题，壳体的厚度是很重要的参数，壳体参考面的偏置也可以用于定义 更精确几何信息。另外一种情况是当模拟一个厚度
11、连续变化的壳体时，中面的偏置可 能是重要的，因为此时定义在壳体中面的节点可能是相当困难的。如果一个表面平滑 而另一个表面粗糙，比如在某些飞行器结构中，应用壳体参考面偏置定义在平滑表面 上的节点会是最容易的。 通过指定一个偏置量，可以引入偏置。定义偏置量作为从壳的中面到壳的参考表 面之间的壳体厚度的比值，如图5-4 所示。 图 5-4 对于偏置量为0.5的壳体偏置示意图 5-5 壳的自由度与其参考表面相关，在此处计算所有的动力学方程，包括计算单元的 面积。对于曲壳，大的偏置量可能导致面上积分的误差，会影响到壳截面的刚度、质 量和转动惯量。为了达到稳定性的目的，ABAQUS/Explicit也会按
12、偏置量平方的量级 自动地增大应用于壳单元的转动惯量，在过大的偏置量的动态分析中，这可能会导致 误差。当从壳中面的大偏置量是非常必要时，使用多点约束或刚体约束来代替偏置。 5.2壳体公式厚壳或薄壳 壳体问题一般可以归结为以下两类之一：薄壳问题和厚壳问题。厚壳问题假设横 向剪切变形对计算结果有重要的影响。另一方面，薄壳问题假设横向剪切变形是小到 足以忽略。图5-5(a) 描述了薄壳的横向剪切行为：初始垂直于壳面的材料线在整个变 形过程中保持直线和垂直。因此，横向剪切应变假设为零（0）。图 5-5(b) 描述了 厚壳的横向剪切行为：初始垂直于壳面的材料线在整个变形过程中并不要求保持垂直 于壳面，因此
13、，发生了横向剪切变形（0） 。 图 5-5在（a）薄壳和（ b）厚壳中的横截面行为 按照将壳单元应用于薄壳和厚壳问题来划分，ABAQUS 提供了多种壳单元。通用 目的的 （general-purpose）壳单元对于应用于薄壳和厚壳问题均有效。在某些特殊用途 的情况下，通过应用在ABAQUS/Standard 中的特殊用途壳单元可以获得增强的性能。 5-6 特殊用途的壳单元可归结为两类：仅为薄壳单元和仅为厚壳单元。所有特殊用途 的壳单元提供了可以有任意大的转动，但是限于小应变。薄壳单元强化了Kirchhoff 约束；即垂直于壳体中面的平截面保持垂直于壳中面，这样，或者是在单元公式的解 析解答（
14、STRI3 单元）或者是在通过罚函数约束的数值解答方面，Kirchhoff 约束得到 了强化。厚壳单元是二阶四边形单元，在小应变应用中，对于使解答沿壳的跨度方向 上平滑地变化的载荷，这种单元能产生比通用目的的壳单元更加精确的结果。 如何判断一个给定的应用是属于薄壳还是厚壳问题，我们可以提供几点指南。对 于厚壳，横向剪切变形是重要的，而对于薄壳它则可以忽略不计。通过厚度与跨度的 比值，可以评估在壳体中横向剪切的显著性。对于由单一各向同性材料组成的壳体， 当比值大于1/15 时可认为是厚壳；如果比值小于1/15，则可认为是薄壳。这些估计是 近似的；用户始终应当检验在模型中横向剪切的影响，以验证壳行
15、为的假设。在复合 材料层合壳结构中，由于横向剪切变形较为显著，对于应用薄壳理论，这个比值必须 是更小一些。采用高度柔软中间层的复合材料层合壳（即“三明治”复合）具有非常 低的横向剪切刚度，所以它们几乎总是要作为厚壳来模拟；如果平截面保持平面的假 设失效，则应采用实体单元。关于如何检验应用壳体理论的有效性的详细信息，请参 阅 ABAQUS 分析用户手册的第15.6.4 节“ Shell sectionbehavior” 。 通用目的壳单元和仅为厚壳单元考虑了横向剪力和剪切应变。对于三维单元，提 供了对于横向剪切应力的评估。这些应力的计算忽略了在弯曲和扭转变形之间的耦合 作用，并假设材料性质和弯矩
16、的空间梯度很小。 5.3壳的材料方向 与实体单元不同，每个壳体单元都使用局部材料方向。各向异型材料的数据（如 纤维增强复合材料）和单元输出变量（如应力和应变）都是以局部材料方向的形式定 义的。在大位移分析中，壳面上的局部材料坐标轴随着各积分点上材料的平均运动而 转动。 5.3.1默认的局部材料方向 局部材料的1 和 2 方向位于壳面内，默认的局部1 方向是整体坐标1 轴在壳面上 的投影。如果整体坐标1 轴是垂至于壳面，则局部1 方向则是整体坐标3 轴在壳面上 5-7 的投影。局部2 方向垂直于位于壳面中的局部1 方向，因此，局部1 方向、 2 方向和 壳体表面的正法线构成右手坐标系（如图5-6
17、 所示）。 图 5-6默认的壳体局部材料方向 局部材料方向的默认设置有时可能会产生问题；关于这方面的一个例子是圆柱形 壳体，如图5-7 所示。对于图中大多数单元，其局部1 方向就是环向。然而，有一行 单元垂直于整体1 轴，对于这些单元，局部1方向为整体3 轴在壳上的投影，使该处 的局部1 方向变为轴向，而不是环向。沿局部1 方向的应力 11的等值线图看起来就 会非常奇怪，由于大多数单元的 11为环向应力，而部分单元的11为轴向应力。在这 种情况下，对于模型需要定义更适合的局部方向，如在下一节中所讨论的。 图 5-7在圆柱形壳体中默认的局部材料1 方向 5.3.2建立可变的材料方向 应用局部的直
18、角、圆柱或者球坐标系，可以代替整体的笛卡尔坐标系，如图5-8 所示。定义一个局部(x, y,z) 坐标系的方向， 并使局部坐标轴的方向与材料方向一致。 为此，你必须先指定一个局部轴（1、 2 或 3），它最接近垂直于壳体的1 和 2材料方向， 并绕轴旋转。 ABAQUS 按照坐标轴的循环顺序（1,2,3）和按照你的选择将坐标轴投影 5-8 到壳体上，从而构成材料的1 方向。例如，如果你选择了x轴，ABAQUS 将 y 轴投 影到壳体上而构成材料的1 方向。由壳法线和材料1 方向的叉积来定义材料的2 方向。 一般情况下，最终的材料2 方向和其它局部坐标轴的投影，如本例中的z轴，对于曲 壳将是不一
19、致的。 图 5-8局部坐标系的定义 如果这些局部坐标轴没有建立理想的材料方向，你可以指定一个绕所选择轴的转 动，另外两个局部坐标轴在投影到壳面之前将按照该转动量进行转动，以得到最终的 材料方向。为了使投影容易实现，所选择的轴应尽可能地接近壳的法线。 例如，如果在图5-7 中的圆柱中心线与整体坐标3 轴一致，局部材料方向可以这 样定义，使局部材料1 方向总是沿着圆环方向，并使相应的局部材料2方向总是沿着 轴方向。其过程描述如下。 定义局部材料方向: 1.从 Property（特性）模块的主菜单栏中，选择Tools -Datum ，并定义一个圆 柱数据坐标系。 5-9 2.选择 Assign -M
20、aterial Orientation给部件赋于一个局部材料方向。当提示选 择坐标系时，选择在上一步中定义的数据坐标系，近似的壳体法线方向是 Axis-1 （1-轴）；不需要额外的转动。 5.4选择壳单元 对于需要考虑薄膜作用或含有弯曲模式沙漏的问题，以及具有平面弯曲的问题，当 希望得到更精确的解答时，可使用ABAQUS/Standard中的线性、有限薄膜应变、 完全积分的四边形壳单元（S4）。 线性、有限薄膜应变、减缩积分、四边形壳单元（S4R）是强健的，并适合应用于 广泛的问题。 线性、有限薄膜应变、三角形壳单元（S3/S3R）可作为通用目的的壳单元使用。因 为在单元中是常应变的近似场，求
21、解弯曲变形或者高应变梯度时可能需要精细的网 格划分。 在复合材料层合壳模型中，为了考虑剪切变形的影响，采用适合于模拟厚壳问题的 单元（ S4, S4R,S3/S3R, S8R）；并检验平截面保持平面的假定是否满足。 四边形或三角形的二次壳单元，对于应用于一般的小应变薄壳是很有效的，这些单 元对于剪力自锁或薄膜自锁都不敏感。 如果在接触模拟中一定要使用二阶单元，不要使用二阶三角形壳单元（STRI65）， 而要采用9节点的四边形壳单元（S9R5）。 对于规模非常大但仅经历几何线性行为的模型，使用线性、薄壳单元（S4R5）通常 比通用目的的壳单元更节约计算成本。 对于包含任意的大转动和小薄膜应变的显
22、式动态问题，小薄膜应变单元是有效的。 5.5例题：斜板 在这个例题中， 要求模拟如图5-9 所示的板， 该板与整体1 轴的夹角为30 ，一端 固支，另一端约束住了，但仅可沿平行于板轴的轨道运动。当板在均布载荷作用下， 需要确定跨中的挠度，也要评估线性分析对于该问题是否有效。计算分析将采用 ABAQUS/Standard 。 5-10 图 5-9斜板示意图 5.5.1前处理用 ABAQUS/CAE建立模型 为了进行这个模拟，应用ABAQUS/CAE创建整体模型。在本手册的在线文档第 A.3 节“Skewplate”提供了输入文件。当通过ABAQUS/CAE运行这个输入文件时， 将创建关于该问题的
23、完整的分析模型。根据下面给出的指导如果你遇到困难，或者如 果你希望检查你的工作，则可以运行这个输入文件。在附录A“ Example Files”中， 给出了如何提取和运行输入文件的指导。 如果你没有进入ABAQUS/CAE或者其它的前处理器，可以人工创建关于这个问 题的输入文件， 关于这方面的讨论， 见 GettingStarted withABAQUS/Standard ： Keywords Version，第 5.5 节“Example：Skewplate” 。 5-11 在开始建模前， 需要确定所使用的量纲系统。本例给出的尺寸单位为cm，但是给 出的载荷和材料属性的单位为MPa 和 GP
24、a。由于它们在量纲上不匹配，在模型应用中 必须选择一致性的量纲系统，并在输入数据中进行必要的换算。在以下的讨论中，采 用牛顿、米、千克和秒的国际量纲系统。 定义模型的几何形状 启动 ABAQUS/CAE ，进入 Part（部件） 模块， 并创建一个带有平面壳体基本 特征的三维变形体，命名该部件为Plate，并指定一个大致的部件尺寸为4.0 。 在下面的过程中，列出了一个所建议的创建部件几何形状的步骤： 绘制板几何： 1.在绘图区，用Create Lines（创建线） :Connected（连接）工具绘制一 条长度为0.4 m 的竖直线。 2.应用 Create Construction（创建辅
25、助项）: Line at anAngle（角度辅助线） 工具，通过直线的每个端点，创建一条与水平线成30 角方向的辅助线。 3.应用 Create Isolated Point（创建独立点）工具，在竖直线的右边画一个 与该线的水平距离为1.0 m 的独立点，然后通过这一点绘制一条竖直方向的辅 助线。 4.使用 Create Lines: Connected工具， 利用已建立的辅助线之间的交点定位各 个角点，绘制出斜矩形图。 最终的绘图如图5-10 所示。 5.在提示区，点击Done 完成绘制。 5-12 图 5-10绘制板的几何形状 定义材料、截面特性和局部材料方向 9Pa，泊松比 v = 板
26、的材料是各向同性的线弹性材料，其弹性模量E = 3010 0.3。进入 Property（特性）模块并创建材料定义；命名材料为Steel。 在整体坐标系下结构的方向如图5-9 所示。整体笛卡尔坐标系定义了默认的 材料方向，但是该板相对于这个坐标系统是倾斜的，如果使用默认的材料方向， 解释模拟的结果将是非常不容易的，因为沿着材料1 方向上的正应力 11将包含 来自板弯曲产生的轴向应力和与板轴线垂直的横向应力的贡献。如果材料的方向 是倾斜于板的轴线和与横向方向一致，则很容易解释模拟的结果。因此，需要定 义一个局部的直角坐标系，它的局部x方向沿着板的轴向（即与整体坐标系1 轴 的夹角为30 ） ，
27、y 轴方向也位于板的平面内。 定义壳的截面特性和局部材料方向： 1.定义一个均匀的壳体截面，命名为 PlateSection。赋值壳体厚度为0.8E-2， 并将 Steel材料定义到截面上。由于材料是线弹性，所以在分析前指定截面积 分的计算。 2.使用 Create Datum CSYS: 2Lines（两线法）工具，定义一个直角基准 5-13 坐标系，如图5-11 所示。 图 5 11用于定义局部材料方向的基准坐标系 3.从主菜单栏中，选择Assign -Material Orientation，并选择整个部件作为将 应用局部材料方向的区域。在视图中， 选择刚建立的基准坐标系，选择 Axis
28、 3 作为壳体法线近似的方向，无需绕此轴进行额外的旋转。 提示 ：为了检验已经赋值的局部材料方向是否正确，从主菜单栏中，选 择 Tools -Query ，进行关于材料方向的性能查询。 一旦在模型中的部件被网格剖分和创建成单元后，所有的单元变量将被定义在 这个局部坐标系下。 4最后，将截面的定义赋值于板。 创建装配件、定义分析步骤和指定输出要求 在 Assembly（装配件）模块中创建斜板部件的实体（instance）。在退出装配 件模块前定义几何集合体，以便于定义输出要求和边界条件。首先，需要在板的 跨中位置（ Midspan）将平板分割成两个区域以创建几何集合体。 分割板和定义几何集合:
29、1. 应用 Partition Face（切割面） : Shortest Path Between2Points（两点间最 5-14 短路径）工具，分割板为两个半区，采用板的斜边的中点创建切割分区， 如图 5-12 所示。 图 5 12应用分割在板的跨中定义几何集合体 2. 选择 Tools -Set-Create 为跨中创建一个几何集合, 命名为 MidSpan 。类似 地，为板的左和右边界各创建一个几何集合，分别命名为EndA 和 EndB 。 提示：通过选择主菜单栏中的Tools -Set -Manager ，可以查看已有的几何 集合。在Set Manager （集合管理器）对话框中，双
30、点击集合的名称。所选 择的集合在视图窗中以高亮度显示，若需要时可以对其定义进行编辑。 接下来，在Step模块中创建一个一般静态（static, general）分析步，命名为 Apply Pressure，并给出下面的步骤描述：Uniform pressure (20kPa) load 。接受所有的对分析步的默认设置。 在所需要的输出中，有节点位移、 约束反力和单元应力作为场变量输出数据， 这些数据将应用于在Visualization （可视化）模块中创建变形形状图、等值线图 和数据报表。 你也希望将跨中的位移写入历史数据，以便在可视化模块中创建X-Y 曲线图。 修改默认的输出请求: 1.编辑
31、场变量输出设置，因此只有整个模型的节点位移、约束反力和单元应力作 5-15 为场变量数据写入到.odb文件中。 2.编辑历史输出设置，因此只有MidSpan几何集合的节点位移作为历史变量数 据写入到 .odb文件中。 施加边界条件和荷载 如图 5-9 所示，板的左端是完全固支的；右端约束住了，但仅可沿平行于板 的轴向的轨道移动。由于右端边界条件的方向与整体坐标轴不一致，必须定义一 个局部坐标系，使一个轴与板的走向一致。你可以利用前面为定义局部材料方向 而创建的坐标系。在操作之前，切换到Load（载荷）模块。 在局部坐标系中定义边界条件： 1.选择BC-Create，在ApplyPressure
32、分析 步中 ， 定义 Displacement/Rotation（位移 /转动）力学边界条件，命名为Railboundary condition 。 在本例中，将边界条件定义在集合上，而不是直接在图形窗中选定区域。因此， 当提示选择施加边界条件的区域时，在图形窗中的提示区，点击 Sets（集合） 。 2. 从显示的RegionSelection（区域选择）对话框中，选择集合EndB 。选中 Highlight selectionsinviewport（高亮度显示选择区域）以确保选择了正确 的集合，此时板的右侧边会高亮度显示。点击Continue。 3. 在 EditBoundaryCondit
33、ion（编辑边界条件）对话框中，点击Edit （编辑） 以指定边界条件将要采用的局部坐标系。在图形窗中，选择前面为定义局部材 料方向而创建的坐标系，局部1 方向与板的轴向一致。 4. 在 Edit Boundary Condition对话框中，固定除了U1 以外的所有自由度。 现在，板的右边界被限制住了，仅能沿平行于板的轴向移动。一旦对板的模型 剖分了网格和生成了节点，所有打印的与这个区域相关的节点输出值（位移、 速度、约束反力等）都将被定义在这个局部坐标系中。 通过固定在平板左端（集合 EndA ）全部的自由度， 完成了边界条件的定义。 命名这个边界条件为Fix left end。对于这个边
34、界条件采用默认的整体方向。 最后， 在壳体的上部定义一均布压力载荷，命名为 Pressure，应用 Shift+ 点击， 选择部件的两个分区，并选定壳体的顶面（Magenta （紫红色）箭头）作 为施加压力载荷的面。为了更加清楚的区分板的顶面，可以旋转视图。指定载荷 5-16 值为 2.E4 Pa。 创建网格和定义作业 图 5-13 显示了对于该模拟所建议的网格剖分。 图 5-13对于斜板模拟所建议的网格设计 在选择单元类型前，必须回答以下问题：板是薄或者是厚？应变是小或者是 大？板是相当薄的，采用的厚度与最小跨度之比为0.02（厚度为0.8 cm，最小的 跨度为 40cm）。当我们现在不能轻
35、易地预测出在结构中的应变量级时，我们认为 应变是小的。基于这个信息，选用二次壳单元（S8R5），其原因是在小应变模拟 中，对于薄壳这类单元将给出精确的结果。关于壳单元选择的进一步详细内容， 请参阅 ABAQUS 分析用户手册的第15.6.2 节“Choosing ashell element” 。 进入 Mesh（网格）模块，采用整体单元的尺度为0.1 ，在部件上播撒种子。 从主菜单栏中， 选择 Mesh -Controls为本模型指定结构化网格剖分技术，应用 每节点有5个自由度的二次、减缩积分壳单元（S8R5），来创建一个四边形网格。 在操作前，重新命名模型为Linear。该模型构成了应用在
36、后面第8 章“非 线性”中所讨论的斜板例题的基础。 进入 Job（作业）模块，并定义一个命名为SkewPlate的作业，采用如下 的描述： 5-17 LinearElastic SkewPlate.20 kPa Load. 将模型保存到名为SkewPlate.cae的模型数据库文件中。 提交作业进行分析，并监控求解过程；修正由分析求解程序发现的任何模拟 错误，并调查引起任何警告的原因。 5.5.2后处理 这一节将讨论应用ABAQUS/CAE进行后处理。对于壳体分析结果的可视化，等 值线图（ contour）和矢量图（ symbol）都是很有用的。由于在第4章“应用实体单元” 中已经详细地讨论了
37、等值线图，这里将使用矢量图。 在位于工具栏的Module列表中，点击Visualization进入可视化模块，然后， 打开由该作业创建的.odb文件（ SkewPlate.odb） 。 ABAQUS/CAE默 认 绘 制 出 模 型 的 快 图 。 从 主 菜 单 栏 中 ， 通 过 选 择 Plot -Undeformed Shape或者点击在工具箱中的工具，绘制未变形模型图。 单元法线 应用未变形图检验模型的定义。对于斜板模型的单元法线，检验其定义是否 正确，并指向在正3方向。 显示单元法线: 1. 在提示区中，点击Undeformed Shape Plot Options（未变形图选项）
38、。 弹出 UndeformedShape Plot Options对话框。 2. 设置显示格式（renderstyle）为 Shaded （阴影图）。 3. 点击 Normals （法线）页。 4. 选中 Show normals（显示法线），并接受默认的Onelements设置。 5. 点击 OK 使设置生效并关闭对话框。 默认视角是等视图，可以应用视图（view）菜单中的选项或工具栏中的视图 工具（如）改变视图。 5-18 改变视图 : 1. 从主菜单栏中，选择View -Specify 。 弹出 SpecifyView（设置视图）对话框。 2. 从方法列表中，选择Viewpoint（视点
39、法）。 3. 键入视点矢量的X-、Y- 、Z- 坐标值分别为-0.2, -1、 0.8 ，和向上矢量的坐 标为 0,0,1。 4. 点击 OK。 ABAQUS/CAE显示出指定视图的模型，如图5-14 所示。 图 5-14斜板模型的壳单元法线 矢量图（ symbol plots） 矢量图用矢量显示了指定的变量，矢量的始端为节点或单元积分点。除了一 些主要是非力学输出变量和在节点上储存的单元结果（如节点力）之外，大多数 用张量和矢量表达的变量都可以绘制出矢量图。相应的箭头尺寸表明了结果值的 相应量级，矢量指向是沿着结果的整体方向。可以绘制出如位移(U)、约束反力 (RF) 等结果变量的矢量图；或
40、者可以绘制出这些变量的每个分量的矢量图。 5-19 生成位移矢量图： 1. 从主菜单栏中，选择Result -Field Output。 弹出 FieldOutput（场变量输出）对话框；默认选择Primary Variable（主要 变量）选项页。 2. 从输出变量列表中，选择位移（U）。 3. 从分量表中，选择3 方向位移（ U3）。 4. 点击 OK。 弹出 Select Plot Mode（选择绘图方式）对话框。 5. 选中 Symbol （矢量图），并点击OK 。 ABAQUS/CAE在模型变形图中显示出在3方向位移的矢量图。 6.为了修改矢量图的属性，在提示区，点击Symbol O
41、ptions（矢量图选项）。 弹出 Symbol Plot Options对话框；默认选择为Basic （基础）页。 7. 为了在未变形的模型中绘制矢量图，点击 Shape （形状） 页，选中 Undeformed shape （未变形形状图）。 8. 点击 OK 使设置生效并关闭对话框。 显示出关于未变形模型的矢量图，如图5-15 所示。 图 5-15位移矢量图 5-20 应用矢量图能够绘制出张量变量的主值，诸如应力。主应力值矢量图在每一 个积分点处用三个矢量来表示，每个矢量代表一个主值，其方向沿着相应的主方 向。压缩值的箭头指向积分点，而拉伸值的箭头背向积分点。也可以绘制出单独 的主值。
42、生成主应力的矢量图: 1. 从主菜单栏中，选择Result -FieldOutput。 弹出 FieldOutput对话框。 2. 从输出变量列表中，选择S（应力）；并从其不变量（invariants）列表中，选择 Max.Principal（最大主应力）。 3. 点击 OK 完成选择并关闭对话框。 ABAQUS/CAE显示出主应力的矢量图。 4. 为了改变箭头长度，在提示区，点击Symbol Options。 弹出 Symbol Plot Options对话框。 5. 点击 Color &Style（颜色与样式）页；然后点击Tensor （张量）子页。 6. 设置 Length （长度）选项
43、为Short （短）。 7. 点击 OK 确认设置，并关闭对话框。 显示的矢量图，如图5-16所示。 图 5 16平底面主应力的矢量图 5-21 8. 主应力默认在截面点1 处显示， 若要在其他非默认的截面点处显示应力，可从主 菜单栏中，选择Result -Section Points，弹出 SectionPoints（截面点）对 话框。 9. 选择理想的非默认的截面点绘图。 10. 在复杂模型中，单元网格线的存在可能遮掩了矢量图。为了消除所显示的单元 网格线，在 Symbol Plot Options对话框中的Basic 页中，选择 Feature edges （特征边界）。图5-17 显示
44、了默认截面点处主应力的矢量图，图中仅显示了板 的特征边界。 图 5-17特征边界图显示主应力的矢量 材料方向 ABAQUS/CAE也可以使单元材料方向可视化。这个功能是特别有助于确保 材料方向的正确性。 绘制材料方向: 1. 从主菜单栏中，选择Plot -MaterialOrientations；或在工具栏中点击工 具。 5-22 在变形后的模型中绘制出了材料的方向图，默认绘出了表示材料方向的三向坐 标，但是没有箭头。 2. 为了应用带箭头的坐标表示，在提示区， 点击 Material Orientation Options（材 料方向图选项）。 弹出 Material Orientation
45、Plot Options（材料方向图选项）对话框。 3. 点击 Color &Style页，再点击Triad （三向坐标系）页。 4. 设置 Arrowhead（箭头）选项，为三个方向增添上箭头。 5. 点击 OK，完成设置并关闭对话框。 6. 从主菜单栏中，选择View -Views Toolbox；或在工具箱中点击工具。 弹出 Views （视图）工具箱。 7. 使用在工具箱中预先设置的视角来显示板，如图5-18 所示，在这个图中，关闭 透视效果（ perspective）。点击在工具箱中的工具关闭透视效果。 默认的材料1 方向为蓝色，材料2 方向为黄色，如果显示的话，材料3 方向为 红色
46、。 图 5-18在板中的材料方位图 5-23 根据报表数据评估结果 通过检验打印的数据，进行其余的后处理工作。借助于显示组（display groups），对于整个模型的单元应力（分量S11、S22 和 S12），在约束点处的约 束反力（集合EndA 和 EndB ）和跨中节点的位移（集合MidSpan ），创建表格 数据报告。应力数据显示如下。 位置 Loc 1 和 Loc 2 标识了在单元中计算应力的截面点的位置。Loc 1（对 应于截面点1）位于壳表面的SNEG 面上，而 Loc 2 （对应于截面点3）位于 SPOS 面上。对于单元采用了局部材料方向：应力提供在局部坐标系中。 检查小应变
47、假设对于本模拟是否有效。对应于应力峰值的轴向应变为 5-24 110.0079。如果它是小于4% 或 5% ，可作为典型的小应变考虑，所以 0.0079 的应变属于应用单元S8R5 模拟的非常合适的范围。 在下面的表中，观察约束反力和反力矩： 在整体坐标系下输出了约束反力。检查约束反力和反力矩，对应于施加的外 载 荷 的 合 力 是 否 为 零 。 在3 方 向 的 非 零 约 束 反 力 等 于 竖 向 的 压 力 载 荷 （20kPa 1.0m0.4m）。除了约束反力，在转动自由度的约束处，压力载荷还 5-25 引起了自平衡的约束反力矩。 5.6相关的 ABAQUS例题 ABAQUS实例问
48、题手册（ABAQUSExample Problems Manual）的第2.1.6 节： “Pressurized fuel tank with variable shell thickness （变厚度压力燃料储罐）”。 ABAQUS 基准手册 （ABAQUS BenchmarksManual ）的第 1.1.2 节：“Analysis of an anisotropic layeredplate（各向异型层合板的分析）”。 ABAQUS 基准手册（ ABAQUSBenchmarks Manual ）的第 1.2.4 节：“Buckling of a simply supported squ
49、are plate（简支方板的屈曲）”。 ABAQUS 基准手册 （ABAQUS BenchmarksManual ）的第 2.3.1 节：“The barrel vault roof problem（圆柱形拱顶问题）”。 5.7建议阅读的文献 下面的参考文献提供了关于壳体理论的理论和计算方面更加深入的内容。 基本壳体理论 Timoshenko, S., Strength of Materials: Part II,KriegerPublishing Co., 1958. Timoshenko, S. and S. W. Krieger, Theory of Plates and Shells
50、, McGraw-Hill,Inc., 1959. Ugural, A. C., Stresses inPlates and Shells,McGraw-Hill,Inc.,1981. 基本计算壳体理论 Cook, R.D., D. S. Malkus, and M. E.Plesha, Concepts and Appliccations ofFinite Element Analysis, JohnWiley& Sons.1989. Hughes, T.J.R.,The Finite Element Method ,Prentice-Hill,Inc.,1987. 高等壳体理论 Budia
