ansys在处理体壳连接问题中的应用.pdfVIP

ANSYS 在处理体壳连接问题中的应用 徐怀海 河海大学土木工程系，南京 (210098) E-mail：haizhihun629611@163.com 摘 要：实体单元和壳单元之间的连接由于本身自由度的不同而存在自由度不协调问题，这 个问题是有限元计算中必须解决的问题。本文介绍了ANSYS 在处理实体单元和壳单元之间 的连接问题的几种方法，并通过实例分析比较，得出了几种方法的适用范围以及计算结果的 精确程度。 关键词：ANSYS ；实体单元；壳单元；连接问题 1．引言 在很多工程结构中，往往既有实体结构部分又有板壳结构部分，如苏通长江公路大桥主 航道索塔采用的钢锚箱与混凝土塔壁组合结构，索塔采用混凝土、钢锚箱采用钢材，这样充 分利用了两种材料，既满足索塔的受力要求，又降低了工程的造价。 当对这类结构进行数值分析时，为了使单元离散更切实际和有效，针对不同类型结构宜 采用相应合理的单元模式，因此，体壳组合结构模型是工程结构，特别是大型复杂组合结构 有限元分析中经常需要采用的模型，在体壳组合结构模型中，结构的一部分离散为实体单元， 另一部分离散为壳单元。 在壳单元中，每个节点具有六个自由度：沿 X 、Y 和 Z 方向的平动自由度 UX 、 UY 、 UZ 以及绕X 、Y 和 Z 轴的转动自由度 ROTX 、 ROTY 、ROTZ ，而实体单元的每个节点只 有 3 个沿着 X 、Y 和 Z 方向的平动自由度 UX 、UY 、UZ 。当实体单元和壳单元连接在一起 共同工作时，即存在自由度不协调问题。在使用 ANSYS 进行有限元分析时，如果实体单元 与壳单元相互之间不进行任何处理就直接进行运算，则系统会提示该运算不能进行或者会出 现节点位移超过最大值的警告，因此需要研究和解决两种不同类型单元的连接问题。 实体单元和壳单元的连接方法一般可以分为两大类： 一种是过渡单元法，Surana[1-4] 曾先后提出用于轴对称应力分析和三维应力分析的过渡 单元，以解决不同类型单元的连接。他所考虑的实体单元是等参单元，壳单元是退化壳单元， 过渡单元实际上是两种单元的结合。在过渡单元中，一部分位移插值函数取自实体单元，另 一部分位移插值函数取自壳单元，因而也将壳单元的法向约束假设（垂直于中面的法线在变 形后既不伸长也不缩短）引入了过渡单元位移模式，这一强制性法向约束会在过渡单元内产 生不合理的附加应力。所以 Surana 在其文章中指出：在过渡单元中，如何确定应变和应力 尚是一个需要研究的问题，在他的算例中，一般都规定材料的泊松比为零，可能正是为了回 避上述问题。 另一种方法是建立实体单元节点与壳单元节点的多点约束方程，通过罚函数法或直接引 入法使多点约束方程得到满足。 罚函数法是通过引入附加条件构造修正泛函，将场函数的问题转化为求修正泛函的驻值 问题。其具体处理方法为：首先需要通过罚参数将约束方程引入到系统的能量泛函中，然后 通过求解驻值条件可以得到满足约束方程，即满足接触面上位移协调条件的系统位移场。罚 函数只能使多点约束方程近似得到满足，而且罚参数的选取是一件很难把握的事，选取不当 -1- 对结果影响很大，因此罚函数法不是一种好的方案。 建立实体单元节点与壳单元节点的多点约束方程的另一种方法是直接引入法。这种方法 是通过接触面来定义。在接触面上，每个节点位移参数之间都存在着约束方程，取接触面上 实体单元的节点位移参数作为独立变量，而接触面上壳单元的节点位移通过接触面的位移约 束条件用实体单元中的节点位移参数来表示。直接引入法将与实体单元连接的壳单元的刚度 矩阵和载荷向量进行修正，最终集合成系统的刚度矩阵和载荷向量，从而得到系统的求解方 程组。直接引入法能使多点约束方程精确地得到满足，是值得采用的方法。 2 ．ANSYS 中的处理方法 对于实体单元和壳单元的连接问题，在 ANSYS 中有不同的处理方式，并且各种方式应 用的范围也不同。本文主要考虑了其中的三种不