【JY】 ABAQUS子程序UEL的有限元原理与应用

不等待

即关注

【简述ABAQUS中UEL子程序】

ABAQUS作为成熟的商用有限元软件，可为高级用户提供特定的分析需求。ABAQUS常见的二次开发子程序包括：UMAT、VUMAT、UGENS、UEL和VUEL等。其中UEL/VUEL分别适用于ABAQUS的Standard/Explicit求解器。只有清楚有限元分析的基本原理，才能够较好地了解其分析的力学原理，才能对特定的分析需求编写合适的分析单元。

本文对ABAQUS/Standard中UEL子程序的有限元原理及编写原则进行初步梳理，并编写了平面三角形单元静力分析的UEL子程序，以加深对UEL的理解与认识。

【ABAQUS的UEL单元之有限元基本原理简述】

当用户需要用ABAQUS/Standard单元库中没有的单元进行分析时，可通过ABAQUS提供的UEL子程序接口进行二次开发，编写适用于特定分析的单元。下面以平面三角形单元为例，给出编写UEL单元的基本流程。

系统的势能Π由下式给出

其中为单元的应变能。

将应变-位移关系ε = Bq代入U_e中，由此对于势能方法形成的刚度矩阵如下：

其中t_e为单元厚度、A_e为单元面积。

有限元分析的第一步是将实体离散化为多个单元，此后构造出各个单元的单元刚度矩阵，在将单元刚度矩阵集成整体的刚度矩阵，最后通过整体刚度矩阵建立平衡方程从而求解各个节点的位移、应力、应变等响应。因此，根据有限元的分析原理，编写UEL的最终目的即是形成目标单元的单元刚度矩阵。以下对平面常应变三角形单元（CST）进行有限元分析，依据有限元分析思路编写UEL子程序，以初步了解有限单元的分析原理及UEL编写流程。

图2 平面三角形单元示意图

上图所示为平面三角形单元，由图示信息可得其节点坐标列阵u与位移列阵q分别为

1.形函数

设三角形单元节点1、2、3对应的形函数分别为N1、N2、N3，形函数满足条件：N1+N2+N3=1（在点i处，形函数Ni值为1，其余形函数值为0），形函数采用自然坐标ξ和η描述，有

2.高斯积分点及权重

此平面三角形单元所选取的高斯积分点及权重如下：

表1高斯积分点及权重

3.形函数与单元位移关系

得到形函数与单元节点位移后，根据等参元表示方法，单元内任意一点的位移都可用形状函数和未知结点位移场进行表示，有

可表示为u=Nq，其中N为