壳单元VUEL子程序，inp文件出错，请帮忙看看？

图样图乃咿唔

最近在写学VUEL的编写，已经走通了一个二维平面应力的例子，但是想写壳单元的VUEL出现了下面错误，（）inp文件没有写对，都没法进行子程序调试):
***ERROR: THE SPECIFIED DEGREES OF FREEDOM FOR THE USER ELEMENT 1 ARE
           INCORRECT. IF ONE TRANSLATIONAL/ROTATIONAL DEGREE OF FREEDOM IS
           INCLUDED AT A NODE, ALL TRANSLATIONAL/ROTATIOANAL DEGREES OF FREEDOM
           MUST BE INCLUDED AT THAT NODE. PLEASE SEE THE USERS MANUAL FOR THE
           REQUIREMENTS IN ACTIVATING DEGREES OF FREEDOM ON USER ELEMENTS IN
           Abaqus/Explicit.
请问这个是什么原因？（似乎是自由度定义的问题）可查看二楼.inp文件,自由度定义应该是没有错的，不知道是什么原因，谢谢了。

图样图乃咿唔

*Heading
** Job name: Job-1 Model name: Model-1
** Generated by: Abaqus/CAE 6.13-1
*Node
      1,           0.,           0.,           0.
      2,          0.5,           0.,           0.
      3,           1.,           0.,           0.
      4,          1.5,           0.,           0.
      5,           2.,           0.,           0.
      6,           0.,          0.5,           0.
      7,          0.5,          0.5,           0.
      8,           1.,          0.5,           0.
      9,          1.5,          0.5,           0.
      10,           2.,         0.5,           0.
      11,           0.,          1.,           0.
      12,          0.5,          1.,           0.
      13,           1.,          1.,           0.
      14,          1.5,          1.,           0.
      15,           2.,          1.,           0.
*Nset, Nset=allnodes
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15
*NSET, NSET=RIGHTNODE
5, 10, 15
*NSET, NSET=LEFTNODE
1, 6, 11
*User Element,Nodes=3,Type=vu2,Properties=10,Coordinates=3,Variables=28
1,2,3,4,5,6
*Element, Type=VU2, Elset=AllElements
1,  1,  2,  7
2,  7,  6,  1
3,  2,  3,  8
4,  8,  7,  2
5,  3,  4,  9
6,  9,  8,  3
7,  4,  5, 10
8, 10,  9,  4
9,  6,  7, 12
10, 12, 11,  6
11,  7,  8, 13
12, 13, 12,  7
13,  8,  9, 14
14, 14, 13,  8
15,  9, 10, 15
16, 15, 14,  9
*Uel Property, Elset=AllElements
80.9e9,0.3,2750,50.0,118.0,0.00193,0.20,6500,0.25,1.27e-3
** ----------------------------------------------------------------
*Boundary
LeftNode, 1, 1
** LeftNode, 2, 2
** LeftNode, 3, 3
** LeftNode, 4, 4
** LeftNode, 5, 5
** LeftNode, 6, 6
**
** STEP: Step-1
**
*Amplitude, Name=LoadAmp, Definition=Smooth Step
0.0, 0.0, 1.0, 1.0
*Step, nlgeom=yes
*Dynamic, Explicit, element
, 1.0
*Variable Mass Scaling, type=below min, Dt=1.0e-3, freq=2
*Bulk Viscosity
0.00, 0.0
*Cload, Amplitude=LoadAmp
RightNode, 2, 1.0e2
**
** OUTPUT REQUESTS
** FIELD OUTPUT: F-Output-1
**
**Output, field, variable=PRESELECT
**Output, History, Freq=1
**Node Output, Nset=ALLNodeS
** U1, V1, A1
*End Step

旧日的轨迹

图样图乃咿唔 发表于 2016-8-9 21:09
*Heading
** Job name: Job-1 Model name: Model-1
** Generated by: Abaqus/CAE 6.13-1

*User Element,Nodes=3,Type=vu2,Properties=10,Coordinates=3,Variables=2
8
1,2,3,4,5,6
第一行的8是指所有节点包含自由度8？一个节点能有8个自由度吗？
第二行表示单元第一个节点除了自由度8以外，还有自由度2，3，4，5，6？

图样图乃咿唔

那个8是上面variables=28

图样图乃咿唔

旧日的轨迹 发表于 2016-8-10 10:51
*User Element,Nodes=3,Type=vu2,Properties=10,Coordinates=3,Variables=2
8
1,2,3,4,5,6

那个8是上面variables=28，可能复制的时候换行了。 就是想表示所有的点自由度都是1,2,3,4,5,6

wu330324

能否传分平面应力的子程序学习一下，必定酬谢啊！

felixyu623

你好楼主，不知道你的问题解决没有？我也想编写一个壳单元的VUEL，但是一直不太理解如何将第三个自由度加入壳单元的有限元列式。我的理解是壳单元的两个转动自由度是根据由壳单元法向确定的节点坐标系来确定的，但是ABAQUS中却有三个转动自由度，那么这样6个自由度的单元是根据什么理论实现的呢？还请楼主帮忙解惑，或者分享一些文献资料？谢谢了！

海枝丹文

想问一下，楼主知不知道怎么在VUEL中控制单元删除的，求回复，谢谢

a56789

楼主能方便分享一下壳单元的子程序吗？QQ：2625007064