ANSYS在处理体壳连接问题中的应用.docxVIP

-1-

ANSYS在处理体壳连接问题中的应用

徐怀海

河海大学土木工程系，南京(210098)

摘要：实体单元和壳单元之间的连接由于本身自由度的不同而存在自由度不协调问题，这个问题是有限元计算中必须解决的问题。本文介绍了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在其文章中指出：在过渡单元中，如何确定应变和应力尚是一个需要研究的问题，在他的算例中，一般都规定材料的泊松比为零，可能正是为了回避上述问题。

另一种方法是建立实体单元节点与壳单元节点的多点约束方程，通过罚函数法或直接引入法使多点约束方程得到满足。

罚函数法是通过引入附加条件构造修正泛函，将场函数的问题转化为求修正泛函的驻值问题。其具体处理方法为：首先需要通过罚参数将约束方程引入到系统的能量泛函中，然后通过求解驻值条件可以得到满足约束方程，即满足接触面上位移协调条件的系统位移场。罚函数只能使多点约束方程近似得到满足，而且罚参数的选取是一件很难把握的事，选取不当

-2-

对结果影响很大，因此罚函数法不是一种好的方案。

建立实体单元节点与壳单元节点的多点约束方程的另一种方法是直接引入法。这种方法是通过接触面来定义。在接触面上，每个节点位移参数之间都存在着约束方程，取接触面上实体单元的节点位移参数作为独立变量，而接触面上壳单元的节点位移通过接触面的位移约束条件用实体单元中的节点位移参数来表示。直接引入法将与实体单元连接的壳单元的刚度矩阵和载荷向量进行修正，最终集合成系统的刚度矩阵和载荷向量，从而得到系统的求解方程组。直接引入法能使多点约束方程精确地得到满足，是值得采用的方法。

2．ANSYS中的处理方法

对于实体单元和壳单元的连接问题，在ANSYS中有不同的处理方式，并且各种方式应用的范围也不同。本文主要考虑了其中的三种不同处理方式，分别为用传统的耦合加约束方程处理接触面、用接触单元处理接触面以及用复合单元处理接触面，其中耦合加约束方程处理接触面和接触单元处理接触面属于直接引入法，复合单元处理接触面属于过渡单元法。

2.1耦合加约束方程处理接触面

当需要迫使两个或者多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起，如在两重复节点间形成铰链或滑移连接等。对于实体单元与壳单元的连接，如果只是用耦合处理接触面，由于实体单元只有平动自由度，则两种单元的节点表现为铰链，为了能够传递转动自由度，必须在接触面上加约束方程。

使用约束方程应注意的事项：所有的约束方程都以小转动理论为基础，因此，它应用在大转动分析中应当限制在约束方程所包含的自由度方向无重大变化的情况；约束方程的出现将产生不可预料的反力和节点力结果。

2.2用接触单元处理接