考虑厚度方向应力的弹塑性大变形壳单元.PDF

张华：考虑厚度方向应力的弹塑性大变形壳单元 279 考虑厚度方向应力的弹塑性大变形壳单元 张华 重庆大学工程力学系．重庆400044 摘要庄板料成形数值模拟中通常使用壳单元，然而，传统的壳单元在计算“回弹”(弹性恢复)变形时 精度很差，其原因之一是由于这些壳单元忽略了厚度方向的应力而代之为平面应力条件。为了改善这个缺 陷．本文提出的壳单元，通过构造厚度方向的二次形函数引入厚度方向的变形．使得三维本构关系可以不 经任何改变被直接使用，并且壳上F表面的载荷对厚度方向应力的影响也能较好地表现出来。通过几个数 值算例证明该壳单元对弹塑性大变形问题分析是有效的。 关键词壳单元；弹塑性；大变形：厚向应力；板料成形：回弹 1 引言 近年来，板料成形数值模拟在汽车冲压件产品和模具的设计中得到了广泛的应用，并且在对板 料成形缺陷一破裂、起皱的预测方面取得了很好的成效。然而，对板料成形缺陷之一的“回弹”的 预测至今还未达到令人满意得效果。其原因之一是由于目前使用的壳单元通常都忽略了厚度方向的 应力而代之为平面应力条件。当板料成形的弯曲半径很小的时候以及板料的上、下表面同时都与模 具发生接触的时候，忽略厚度方向的应力会给“回弹”计算带来较大的误差。因此，为了克服这个 缺陷，需要开发一种能直接使用三维本构关系、正确反映三维应力状态的壳单元。为了能直接使用 个位移和横截面的2个转动)为基础，加入壳厚向纤维的变形，由于弯曲应力沿厚度方向至少呈线 性分布，为了与之协调需要把厚度方向的位移提高到二次插值，使厚度方向的应变为线性分布，增 加了一个非位移变量：另一种方法i4刮是以块体单元为基础，通过把位移转化成上下表面节点的相对 位移．成为相对自由度壳单元，由于在壳的厚度方向的位移插值只是线性关系，其应变是常量，为 了}，J、补这个缺陷，在厚度方向需要加入强制假定应变(EAs)。本文的方法是在壳的中面和上、下表 面都分布节点，使之在壳的厚度方向上形成位移的二次插值关系，这样就可以很自然地得到沿厚度 方向呈线性变化的应变，然后利用壳单元的直线假设，即上、中、下面的节点在变形前后始终处在 同一直线上．将原有的9个自由度变成中面上3个位移自由度、横截面的2个转动自由度和上下表 面节点沿壳厚度方向的2个相对位移自由度，共7个自由度。本文分为5个部分，第一部分是开发 考虑厚度方向应力壳单元的目的；第二部分讨论这种新壳单元的坐标和位移的插值形函数；第三部 分推导该单元的有限元方程．包括单元刚度矩阵和上、下表面载荷的转换方式；第四部分给出几个 数值算例，以证明该壳单元的有效性：第五部分是对本文的总结。 2壳单元的基本公式 为了得到能表现三维应力状态的壳单元t首先从块体单元的概念出发，构造一个壳单元(如F嘻 !塑 堡兰!耋垒生垦主鱼壁垒竺垡望兰苎至兰查兰翌L—————————————一 1所示)。陔壳单元分别由中面、下表面和上表面的节点组成，因此，壳单元内任意一点的坐标可以 用插值关系表示为 【1) x=∑N“(#，口)·[Ⅳ：(f)x：+．v；(f)x：+Ⅳ；(f)x；】 式中：x：、』；和x；分别表示中面、下表面和上表面的节点Ⅱ的坐标，Ⅳ“(f川)表示壳单元面内的 插值形函数，Ⅳ：(f)、．v；(f)和Ⅳ?(f)分别是对应中面、下表面和上表面的节点的沿厚度方向的形 函数．它们是二次插值关系，可以写成 (2) 2：Ⅳ；(f)=号f(f—1)；Ⅳ；(f)=÷f(f+1) Ⅳ二(O=1一f Fig．1中，壳单元的局部坐标系^可以定义为‘71 铲南’铲尚’岛。，”l 。’ 这里 gi：荐0x (