analysis job with a uel subroutine
analysis job with a uel subroutine
(OP)
Dear all
I try to run an abaqus job analysis with a uel subroutine written in c++ and I get
[quote]
ABAQUS Error: The executable C:\ABAQUS\6.6-1\exec\standard.exe
aborted with system error code 142.
[\quote]
So I try to run it with a version of this uel subroutine written in fortran and I get the same error message.Has any one any idea to solve my problem? Thanks a lot in advance.
Here is my input file (cantBeamAbHC3_4x12x12.inp):
[code]
*HEADING
*PREPRINT, MODEL=YES
*NODE
1001001, 0.00000000, 0.00000000, 0.00000000
14001001, 0.00000000, 0.00000000, 30.00000000
1014001, 0.00000000, 30.00000000, 0.00000000
14014001, 0.00000000, 30.00000000, 30.00000000
1001006, 10.00000000, 0.00000000, 0.00000000
14001006, 10.00000000, 0.00000000, 30.00000000
1014006, 10.00000000, 30.00000000, 0.00000000
14014006, 10.00000000, 30.00000000, 30.00000000
1001002, 1.25000000, 0.00000000, 0.00000000
1001005, 8.75000000, 0.00000000, 0.00000000
1002001, 0.00000000, 1.25000000, 0.00000000
1013001, 0.00000000, 28.75000000, 0.00000000
*NGEN,NSET=LX1
1001002,1001005
*NGEN,NSET=LY1
1002001,1013001,1000
*NCOPY,OLDSET=LX1,NEWSET=LX2,CHANGENUMBER=13000,SHIFT
0.00000000, 30.00000000, 0.00000000
*NCOPY,OLDSET=LY1,NEWSET=LY2,CHANGENUMBER=5,SHIFT
10.00000000, 0.00000000, 0.00000000
*NCOPY,OLDSET=LX1,NEWSET=NSETA,CHANGENUMBER=1000,SHIFT
0.00000000, 1.25000000, 0.00000000
*NCOPY,OLDSET=LX1,NEWSET=NSETB,CHANGENUMBER=12000,SHIFT
0.00000000, 28.75000000, 0.00000000
*NFILL,NSET=SuperfInf
NSETA,NSETB,11,1000
*NSET,NSET=SuperfInf
LX1, LX2, LY1, LY2,1001001,1014001,1001006,1014006
*NCOPY,OLDSET=SuperfInf,CHANGENUMBER=13000000,SHIFT
0.00000000, 0.00000000, 30.00000000
*NCOPY,OLDSET=SuperfInf,NEWSET=NSETAA,CHANGENUMBER=1000000,SHIFT
0.00000000, 0.00000000, 1.25000000
*NCOPY,OLDSET=SuperfInf,NEWSET=NSETBB,CHANGENUMBER=12000000,SHIFT
0.00000000, 0.00000000, 28.75000000
*NFILL
NSETAA,NSETBB,11,1000000
**
*USER ELEMENT, NODES=27, TYPE=U111, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U111, ELSET=E111
1001001, 1001001, 2001001, 3001001, 1002001, 2002001, 3002001, 1003001, 2003001,
3003001, 1001002, 2001002, 3001002, 1002002, 2002002, 3002002, 1003002,
2003002, 3003002, 1001003, 2001003, 3001003, 1002003, 2002003, 3002003,
1003003, 2003003, 3003003
*UEL PROPERTY, ELSET=E111
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U112, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U112, ELSET=E112
2001001, 2001001, 3001001, 4001001, 2002001, 3002001, 4002001, 2003001, 3003001,
4003001, 2001002, 3001002, 4001002, 2002002, 3002002, 4002002, 2003002,
3003002, 4003002, 2001003, 3001003, 4001003, 2002003, 3002003, 4002003,
2003003, 3003003, 4003003
*ELGEN,ELSET=E112
2001001, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E112
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U113, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U113, ELSET=E113
12001001, 12001001, 13001001, 14001001, 12002001, 13002001, 14002001, 12003001, 13003001,
14003001, 12001002, 13001002, 14001002, 12002002, 13002002, 14002002, 12003002,
13003002, 14003002, 12001003, 13001003, 14001003, 12002003, 13002003, 14002003,
12003003, 13003003, 14003003
*UEL PROPERTY, ELSET=E113
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U121, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U121, ELSET=E121
1002001, 1002001, 2002001, 3002001, 1003001, 2003001, 3003001, 1004001, 2004001,
3004001, 1002002, 2002002, 3002002, 1003002, 2003002, 3003002, 1004002,
2004002, 3004002, 1002003, 2002003, 3002003, 1003003, 2003003, 3003003,
1004003, 2004003, 3004003
*ELGEN,ELSET=E121
1002001, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E121
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U122, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U122, ELSET=E122
2002001, 2002001, 3002001, 4002001, 2003001, 3003001, 4003001, 2004001, 3004001,
4004001, 2002002, 3002002, 4002002, 2003002, 3003002, 4003002, 2004002,
3004002, 4004002, 2002003, 3002003, 4002003, 2003003, 3003003, 4003003,
2004003, 3004003, 4004003
*ELGEN,ELSET=E122
2002001, 1,1,1, 10,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E122
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U123, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U123, ELSET=E123
12002001, 12002001, 13002001, 14002001, 12003001, 13003001, 14003001, 12004001, 13004001,
14004001, 12002002, 13002002, 14002002, 12003002, 13003002, 14003002, 12004002,
13004002, 14004002, 12002003, 13002003, 14002003, 12003003, 13003003, 14003003,
12004003, 13004003, 14004003
*ELGEN,ELSET=E123
12002001, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E123
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U131, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U131, ELSET=E131
1012001, 1012001, 2012001, 3012001, 1013001, 2013001, 3013001, 1014001, 2014001,
3014001, 1012002, 2012002, 3012002, 1013002, 2013002, 3013002, 1014002,
2014002, 3014002, 1012003, 2012003, 3012003, 1013003, 2013003, 3013003,
1014003, 2014003, 3014003
*UEL PROPERTY, ELSET=E131
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U132, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U132, ELSET=E132
2012001, 2012001, 3012001, 4012001, 2013001, 3013001, 4013001, 2014001, 3014001,
4014001, 2012002, 3012002, 4012002, 2013002, 3013002, 4013002, 2014002,
3014002, 4014002, 2012003, 3012003, 4012003, 2013003, 3013003, 4013003,
2014003, 3014003, 4014003
*ELGEN,ELSET=E132
2012001, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E132
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U133, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U133, ELSET=E133
12012001, 12012001, 13012001, 14012001, 12013001, 13013001, 14013001, 12014001, 13014001,
14014001, 12012002, 13012002, 14012002, 12013002, 13013002, 14013002, 12014002,
13014002, 14014002, 12012003, 13012003, 14012003, 12013003, 13013003, 14013003,
12014003, 13014003, 14014003
*UEL PROPERTY, ELSET=E133
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U211, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U211, ELSET=E211
1001002, 1001002, 2001002, 3001002, 1002002, 2002002, 3002002, 1003002, 2003002,
3003002, 1001003, 2001003, 3001003, 1002003, 2002003, 3002003, 1003003,
2003003, 3003003, 1001004, 2001004, 3001004, 1002004, 2002004, 3002004,
1003004, 2003004, 3003004
*ELGEN,ELSET=E211
1001002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E211
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U212, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U212, ELSET=E212
2001002, 2001002, 3001002, 4001002, 2002002, 3002002, 4002002, 2003002, 3003002,
4003002, 2001003, 3001003, 4001003, 2002003, 3002003, 4002003, 2003003,
3003003, 4003003, 2001004, 3001004, 4001004, 2002004, 3002004, 4002004,
2003004, 3003004, 4003004
*ELGEN,ELSET=E212
2001002, 2,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E212
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U213, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U213, ELSET=E213
12001002, 12001002, 13001002, 14001002, 12002002, 13002002, 14002002, 12003002, 13003002,
14003002, 12001003, 13001003, 14001003, 12002003, 13002003, 14002003, 12003003,
13003003, 14003003, 12001004, 13001004, 14001004, 12002004, 13002004, 14002004,
12003004, 13003004, 14003004
*ELGEN,ELSET=E213
12001002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E213
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U221, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U221, ELSET=E221
1002002, 1002002, 2002002, 3002002, 1003002, 2003002, 3003002, 1004002, 2004002,
3004002, 1002003, 2002003, 3002003, 1003003, 2003003, 3003003, 1004003,
2004003, 3004003, 1002004, 2002004, 3002004, 1003004, 2003004, 3003004,
1004004, 2004004, 3004004
*ELGEN,ELSET=E221
1002002, 2,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E221
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U222, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U222, ELSET=E222
2002002, 2002002, 3002002, 4002002, 2003002, 3003002, 4003002, 2004002, 3004002,
4004002, 2002003, 3002003, 4002003, 2003003, 3003003, 4003003, 2004003,
3004003, 4004003, 2002004, 3002004, 4002004, 2003004, 3003004, 4003004,
2004004, 3004004, 4004004
*ELGEN,ELSET=E222
2002002, 2,1,1, 10,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E222
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U223, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U223, ELSET=E223
12002002, 12002002, 13002002, 14002002, 12003002, 13003002, 14003002, 12004002, 13004002,
14004002, 12002003, 13002003, 14002003, 12003003, 13003003, 14003003, 12004003,
13004003, 14004003, 12002004, 13002004, 14002004, 12003004, 13003004, 14003004,
12004004, 13004004, 14004004
*ELGEN,ELSET=E223
12002002, 2,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E223
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U231, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U231, ELSET=E231
1012002, 1012002, 2012002, 3012002, 1013002, 2013002, 3013002, 1014002, 2014002,
3014002, 1012003, 2012003, 3012003, 1013003, 2013003, 3013003, 1014003,
2014003, 3014003, 1012004, 2012004, 3012004, 1013004, 2013004, 3013004,
1014004, 2014004, 3014004
*ELGEN,ELSET=E231
1012002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E231
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U232, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U232, ELSET=E232
2012002, 2012002, 3012002, 4012002, 2013002, 3013002, 4013002, 2014002, 3014002,
4014002, 2012003, 3012003, 4012003, 2013003, 3013003, 4013003, 2014003,
3014003, 4014003, 2012004, 3012004, 4012004, 2013004, 3013004, 4013004,
2014004, 3014004, 4014004
*ELGEN,ELSET=E232
2012002, 2,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E232
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U233, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U233, ELSET=E233
12012002, 12012002, 13012002, 14012002, 12013002, 13013002, 14013002, 12014002, 13014002,
14014002, 12012003, 13012003, 14012003, 12013003, 13013003, 14013003, 12014003,
13014003, 14014003, 12012004, 13012004, 14012004, 12013004, 13013004, 14013004,
12014004, 13014004, 14014004
*ELGEN,ELSET=E233
12012002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E233
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U311, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U311, ELSET=E311
1001004, 1001004, 2001004, 3001004, 1002004, 2002004, 3002004, 1003004, 2003004,
3003004, 1001005, 2001005, 3001005, 1002005, 2002005, 3002005, 1003005,
2003005, 3003005, 1001006, 2001006, 3001006, 1002006, 2002006, 3002006,
1003006, 2003006, 3003006
*UEL PROPERTY, ELSET=E311
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U312, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U312, ELSET=E312
2001004, 2001004, 3001004, 4001004, 2002004, 3002004, 4002004, 2003004, 3003004,
4003004, 2001005, 3001005, 4001005, 2002005, 3002005, 4002005, 2003005,
3003005, 4003005, 2001006, 3001006, 4001006, 2002006, 3002006, 4002006,
2003006, 3003006, 4003006
*ELGEN,ELSET=E312
2001004, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E312
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U313, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U313, ELSET=E313
12001004, 12001004, 13001004, 14001004, 12002004, 13002004, 14002004, 12003004, 13003004,
14003004, 12001005, 13001005, 14001005, 12002005, 13002005, 14002005, 12003005,
13003005, 14003005, 12001006, 13001006, 14001006, 12002006, 13002006, 14002006,
12003006, 13003006, 14003006
*UEL PROPERTY, ELSET=E313
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U321, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U321, ELSET=E321
1002004, 1002004, 2002004, 3002004, 1003004, 2003004, 3003004, 1004004, 2004004,
3004004, 1002005, 2002005, 3002005, 1003005, 2003005, 3003005, 1004005,
2004005, 3004005, 1002006, 2002006, 3002006, 1003006, 2003006, 3003006,
1004006, 2004006, 3004006
*ELGEN,ELSET=E321
1002004, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E321
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U322, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U322, ELSET=E322
2002004, 2002004, 3002004, 4002004, 2003004, 3003004, 4003004, 2004004, 3004004,
4004004, 2002005, 3002005, 4002005, 2003005, 3003005, 4003005, 2004005,
3004005, 4004005, 2002006, 3002006, 4002006, 2003006, 3003006, 4003006,
2004006, 3004006, 4004006
*ELGEN,ELSET=E322
2002004, 1,1,1, 10,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E322
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U323, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U323, ELSET=E323
12002004, 12002004, 13002004, 14002004, 12003004, 13003004, 14003004, 12004004, 13004004,
14004004, 12002005, 13002005, 14002005, 12003005, 13003005, 14003005, 12004005,
13004005, 14004005, 12002006, 13002006, 14002006, 12003006, 13003006, 14003006,
12004006, 13004006, 14004006
*ELGEN,ELSET=E323
12002004, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E323
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U331, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U331, ELSET=E331
1012004, 1012004, 2012004, 3012004, 1013004, 2013004, 3013004, 1014004, 2014004,
3014004, 1012005, 2012005, 3012005, 1013005, 2013005, 3013005, 1014005,
2014005, 3014005, 1012006, 2012006, 3012006, 1013006, 2013006, 3013006,
1014006, 2014006, 3014006
*UEL PROPERTY, ELSET=E331
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U332, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U332, ELSET=E332
2012004, 2012004, 3012004, 4012004, 2013004, 3013004, 4013004, 2014004, 3014004,
4014004, 2012005, 3012005, 4012005, 2013005, 3013005, 4013005, 2014005,
3014005, 4014005, 2012006, 3012006, 4012006, 2013006, 3013006, 4013006,
2014006, 3014006, 4014006
*ELGEN,ELSET=E332
2012004, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E332
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U333, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U333, ELSET=E333
12012004, 12012004, 13012004, 14012004, 12013004, 13013004, 14013004, 12014004, 13014004,
14014004, 12012005, 13012005, 14012005, 12013005, 13013005, 14013005, 12014005,
13014005, 14014005, 12012006, 13012006, 14012006, 12013006, 13013006, 14013006,
12014006, 13014006, 14014006
*UEL PROPERTY, ELSET=E333
200.000000,0.300000,2.500000,2.500000,2.500000
**
*BOUNDARY
SuperfInf,PINNED
**
*STEP,PERT
*STATIC
*INCLUDE, INPUT=Cloads_4x12x12.inp
*OUTPUT,FIELD
*NODE OUTPUT
UT,COORD
*NODE PRINT,F=0
*ENDSTEP
[\code]
Here is my fortran uel subroutine (AbHC3.for):
[code]
SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,
1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,
2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,
3 NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,
4 PERIOD)
C
INCLUDE 'ABA_PARAM.INC'
PARAMETER ( ZERO = 0.D0, ONE = 1.D0 , HALF = 0.5D0, TWO = 2.D0)
C
DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),
1 SVARS(NSVARS),ENERGY(8),PROPS(*),COORDS(MCRD,NNODE),
2 U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2),
3 PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),
4 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*),
5 JPROPS(*)
DIMENSION H(3), A00(3,3,3), A10(3,3,3), A01(3,3,3), A11(3,3,3),
1 XI(3)
C
INTEGER N, NR, NC, DR, DC
C
C UEL SUBROUTINE FOR A HORIZONTAL TRUSS ELEMENT
C
C SRESID - stores the static residual at time t+dt
C SVARS - In 1-6, contains the static residual at time t
C upon entering the routine. SRESID is copied to
C SVARS(1-6) after the dynamic residual has been
C calculated.
C - For half-step residual calculations: In 7-12,
C contains the static residual at the beginning
C of the previous increment. SVARS(1-6) are copied
C into SVARS(7-12) after the dynamic residual has
C been calculated.
C
YOUNG = PROPS(1)
POISSON = PROPS(2)
H(1) = PROPS(3)
H(2) = PROPS(4)
H(3) = PROPS(5)
C
LAMBDA = YOUNG*POISSON/(ONE+POISSON)/(ONE-TWO*POISSON)
MU = YOUNG/TWO/(ONE+POISSON)
C
A00(1,1,1) = 0.2000000000D0
A00(1,1,2) = 0.1166666667D0
A00(1,1,3) = 0.1666666667D-1
A00(1,2,1) = 0.1166666667D0
A00(1,2,2) = 0.2833333333D0
A00(1,2,3) = 0.1000000000D0
A00(1,3,1) = 0.1666666667D-1
A00(1,3,2) = 0.1000000000D0
A00(1,3,3) = 0.5000000000D-1
A00(2,1,1) = 0.5000000000D-1
A00(2,1,2) = 0.1083333333D0
A00(2,1,3) = 0.8333333333D-2
A00(2,2,1) = 0.1083333333D0
A00(2,2,2) = 0.4500000000D0
A00(2,2,3) = 0.1083333333D0
A00(2,3,1) = 0.8333333333D-2
A00(2,3,2) = 0.1083333333D0
A00(2,3,3) = 0.5000000000D-1
A00(3,1,1) = 0.5000000000D-1
A00(3,1,2) = 0.1000000000D0
A00(3,1,3) = 0.1666666667D-1
A00(3,2,1) = 0.1000000000D0
A00(3,2,2) = 0.2833333333D0
A00(3,2,3) = 0.1166666667D0
A00(3,3,1) = 0.1666666667D-1
A00(3,3,2) = 0.1166666667D0
A00(3,3,3) = 0.2000000000D0
A01(1,1,1) = -0.5000000000D0
A01(1,1,2) = 0.4166666667D0
A01(1,1,3) = 0.8333333333D-1
A01(1,2,1) = -0.4166666667D0
A01(1,2,2) = 0.1250000000D0
A01(1,2,3) = 0.2916666667D0
A01(1,3,1) = -0.8333333333D-1
A01(1,3,2) = -0.4166666667D-1
A01(1,3,3) = 0.1250000000D0
A01(2,1,1) = -0.1250000000D0
A01(2,1,2) = 0.8333333333D-1
A01(2,1,3) = 0.4166666667D-1
A01(2,2,1) = -0.3333333333D0
A01(2,2,2) = 0.0D0
A01(2,2,3) = 0.3333333333D0
A01(2,3,1) = -0.4166666667D-1
A01(2,3,2) = -0.8333333333D-1
A01(2,3,3) = 0.1250000000D0
A01(3,1,1) = -0.1250000000D0
A01(3,1,2) = 0.4166666667D-1
A01(3,1,3) = 0.8333333333D-1
A01(3,2,1) = -0.2916666667D0
A01(3,2,2) = -0.1250000000D0
A01(3,2,3) = 0.4166666667D0
A01(3,3,1) = -0.8333333333D-1
A01(3,3,2) = -0.4166666667D0
A01(3,3,3) = 0.5000000000D0
A10(1,1,1) = -0.5000000000D0
A10(1,1,2) = -0.4166666667D0
A10(1,1,3) = -0.8333333333D-1
A10(1,2,1) = 0.4166666667D0
A10(1,2,2) = 0.1250000000D0
A10(1,2,3) = -0.4166666667D-1
A10(1,3,1) = 0.8333333333D-1
A10(1,3,2) = 0.2916666667D0
A10(1,3,3) = 0.1250000000D0
A10(2,1,1) = -0.1250000000D0
A10(2,1,2) = -0.3333333333D0
A10(2,1,3) = -0.4166666667D-1
A10(2,2,1) = 0.8333333333D-1
A10(2,2,2) = 0.0D0
A10(2,2,3) = -0.8333333333D-1
A10(2,3,1) = 0.4166666667D-1
A10(2,3,2) = 0.3333333333D0
A10(2,3,3) = 0.1250000000D0
A10(3,1,1) = -0.1250000000D0
A10(3,1,2) = -0.2916666667D0
A10(3,1,3) = -0.8333333333D-1
A10(3,2,1) = 0.4166666667D-1
A10(3,2,2) = -0.1250000000D0
A10(3,2,3) = -0.4166666667D0
A10(3,3,1) = 0.8333333333D-1
A10(3,3,2) = 0.4166666667D0
A10(3,3,3) = 0.5000000000D0
A11(1,1,1) = 0.1333333333D1
A11(1,1,2) = -0.1000000000D1
A11(1,1,3) = -0.3333333333D0
A11(1,2,1) = -0.1000000000D1
A11(1,2,2) = 0.1D1
A11(1,2,3) = -0.6938893904D-17
A11(1,3,1) = -0.3333333333D0
A11(1,3,2) = -0.6938893904D-17
A11(1,3,3) = 0.3333333333D0
A11(2,1,1) = 0.3333333333D0
A11(2,1,2) = -0.1666666667D0
A11(2,1,3) = -0.1666666667D0
A11(2,2,1) = -0.1666666667D0
A11(2,2,2) = 0.3333333333D0
A11(2,2,3) = -0.1666666667D0
A11(2,3,1) = -0.1666666667D0
A11(2,3,2) = -0.1666666667D0
A11(2,3,3) = 0.3333333333D0
A11(3,1,1) = 0.3333333333D0
A11(3,1,2) = -0.6938893904D-17
A11(3,1,3) = -0.3333333333D0
A11(3,2,1) = -0.6938893904D-17
A11(3,2,2) = 0.1D1
A11(3,2,3) = -0.1000000000D1
A11(3,3,1) = -0.3333333333D0
A11(3,3,2) = -0.1000000000D1
A11(3,3,3) = 0.1333333333D1
C
IF(JTYPE.EQ.111)THEN
XI(1) = 1
XI(2) = 1
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.112)THEN
XI(1) = 1
XI(2) = 1
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.113)THEN
XI(1) = 1
XI(2) = 1
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.121)THEN
XI(1) = 1
XI(2) = 2
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.122)THEN
XI(1) = 1
XI(2) = 2
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.123)THEN
XI(1) = 1
XI(2) = 2
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.131)THEN
XI(1) = 1
XI(2) = 3
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.132)THEN
XI(1) = 1
XI(2) = 3
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.133)THEN
XI(1) = 1
XI(2) = 3
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.211)THEN
XI(1) = 2
XI(2) = 1
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.212)THEN
XI(1) = 2
XI(2) = 1
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.213)THEN
XI(1) = 2
XI(2) = 1
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.221)THEN
XI(1) = 2
XI(2) = 2
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.222)THEN
XI(1) = 2
XI(2) = 2
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.223)THEN
XI(1) = 2
XI(2) = 2
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.231)THEN
XI(1) = 2
XI(2) = 3
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.232)THEN
XI(1) = 2
XI(2) = 3
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.233)THEN
XI(1) = 2
XI(2) = 3
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.311)THEN
XI(1) = 3
XI(2) = 1
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.312)THEN
XI(1) = 3
XI(2) = 1
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.313)THEN
XI(1) = 3
XI(2) = 1
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.321)THEN
XI(1) = 3
XI(2) = 2
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.322)THEN
XI(1) = 3
XI(2) = 2
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.323)THEN
XI(1) = 3
XI(2) = 2
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.331)THEN
XI(1) = 3
XI(2) = 3
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.332)THEN
XI(1) = 3
XI(2) = 3
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.333)THEN
XI(1) = 3
XI(2) = 3
XI(3) = 3
C GOTO 999
END IF
C
C1111 = lambda + TWO*mu
C2222 = lambda + TWO*mu
C3333 = lambda + TWO*mu
C1122 = lambda
C1133 = lambda
C2233 = lambda
C1212 = mu
C1313 = mu
C2323 = mu
H123 = H(0)*H(1)/H(2)
H231 = H(1)*H(2)/H(0)
H312 = H(2)*H(0)/H(1)
N = 3
C
IF (LFLAGS(3).EQ.1) THEN
C Normal incrementation
IF (LFLAGS(1).EQ.1 .OR. LFLAGS(1).EQ.2) THEN
C *STATIC
IF (LFLAGS(4).NE.0) THEN
DO I = 1, N
DO J = 1, N
DO K = 1, N
DO P = 1, N
DO Q = 1, N
DO R = 1, N
NR = (I-1)*N*N + (J-1)*N + K
NC = (P-1)*N*N + (Q-1)*N + R
C
C contributo u1ijk u1pqr
C
DR = 1
DC = 1
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H231*C1111*A11(XI(0),I,P)*A00(XI(1),J,Q)*A00(XI(2),K,R) +
* H312*C1212*A00(XI(0),I,P)*A11(XI(1),J,Q)*A00(XI(2),K,R) +
* H123*C1313*A00(XI(0),I,P)*A00(XI(1),J,Q)*A11(XI(2),K,R)
C contributo u1ijk u2pqr
DR = 1
DC = 2
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(2)*C1122*A10(XI(0),I,P)*A01(XI(1),J,Q)*A00(XI(2),K,R) +
* H(2)*C1212*A01(XI(0),I,P)*A10(XI(1),J,Q)*A00(XI(2),K,R)
C contributo u1ijk u3pqr
DR = 1
DC = 3
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(1)*C1133*A10(XI(0),I,P)*A00(XI(1),J,Q)*A01(XI(2),K,R) +
* H(1)*C1313*A01(XI(0),I,P)*A00(XI(1),J,Q)*A10(XI(2),K,R)
C
C contributo u2ijk u2pqr
C
DR = 2
DC = 2
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H231*C1212*A11(XI(0),I,P)*A00(XI(1),J,Q)*A00(XI(2),K,R) +
* H312*C2222*A00(XI(0),I,P)*A11(XI(1),J,Q)*A00(XI(2),K,R) +
* H123*C2323*A00(XI(0),I,P)*A00(XI(1),J,Q)*A11(XI(2),K,R)
C contributo u2ijk u1pqr
DR = 2
DC = 1
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(2)*C1122*A01(XI(0),I,P)*A10(XI(1),J,Q)*A00(XI(2),K,R) +
* H(2)*C1313*A10(XI(0),I,P)*A01(XI(1),J,Q)*A00(XI(2),K,R)
C contributo u2ijk u3pqr
DR = 2
DC = 3
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(0)*C2233*A00(XI(0),I,P)*A10(XI(1),J,Q)*A01(XI(2),K,R) +
* H(0)*C2323*A00(XI(0),I,P)*A01(XI(1),J,Q)*A10(XI(2),K,R)
C
C contributo u3ijk u3pqr
C
DR = 3
DC = 3
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H231*C1313*A11(XI(0),I,P)*A00(XI(1),J,Q)*A00(XI(2),K,R) +
* H312*C2323*A00(XI(0),I,P)*A11(XI(1),J,Q)*A00(XI(2),K,R) +
* H123*C3333*A00(XI(0),I,P)*A00(XI(1),J,Q)*A11(XI(2),K,R)
C contributo u3ijk u1pqr
DR = 3
DC = 1
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(1)*C1133*A01(XI(0),I,P)*A00(XI(1),J,Q)*A10(XI(2),K,R) +
* H(1)*C1313*A10(XI(0),I,P)*A00(XI(1),J,Q)*A01(XI(2),K,R)
C contributo u3ijk u2pqr
DR = 2
DC = 2
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(0)*C2233*A00(XI(0),I,P)*A01(XI(1),J,Q)*A10(XI(2),K,R) +
* H(0)*C2323*A00(XI(0),I,P)*A10(XI(1),J,Q)*A01(XI(2),K,R)
END DO
END DO
END DO
END DO
END DO
END DO
C
DO I = 1, NDOFEL
RHS(I,1) = ZERO
DO J = 1, NDOFEL
RHS(I,1) = RHS(I,1) - AMATRX(I,J)*DU(J,1)
END DO
END DO
C
DO I = 1, NDOFEL
ENERGY(2) = ZERO
DO J = 1, NDOFEL
ENERGY(2) = ENERGY(2) + HALF*DU(I,1)*AMATRX(I,J)*U(J)
* + HALF*U(I)*AMATRX(I,J)*DU(J,1)
* + HALF*DU(I,1)*AMATRX(I,J)*DU(J,1)
END DO
END DO
C
END IF
END IF
ELSE IF(LFLAGS(3).EQ.100)THEN
C Output for perturbations
IF (LFLAGS(1).EQ.1 .OR. LFLAGS(1).EQ.2) THEN
C *STATIC
DO I = 1, NDOFEL
RHS(I,1) = ZERO
DO J = 1, NDOFEL
RHS(I,1) = RHS(I,1) - AMATRX(I,J)*DU(J,1)
END DO
END DO
C
DO I = 1, NDOFEL
ENERGY(2) = ZERO
DO J = 1, NDOFEL
ENERGY(2) = ENERGY(2) + HALF*DU(I,1)*AMATRX(I,J)*U(J)
* + HALF*U(I)*AMATRX(I,J)*DU(J,1)
* + HALF*DU(I,1)*AMATRX(I,J)*DU(J,1)
END DO
END DO
C
DO KVAR = 1, NSVARS
SVARS(KVAR) = RHS(KVAR,1)
END DO
C
END IF
END IF
C
RETURN
END
[\code]
Here is what I get in abaqus command window:
[quote]
C:\Documents and Settings\alorenzom\Documenti\HC3FEM\CantBeamAbHC3withFortSr>aba
qus make library=abhc3
ABAQUS JOB abhc3
Begin Compiling ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
End Compiling ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
Begin Linking ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
Creating library standardU.lib and object standardU.exp
End Linking ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
ABAQUS JOB abhc3 COMPLETED
C:\Documents and Settings\alorenzom\Documenti\HC3FEM\CantBeamAbHC3withFortSr>aba
qus job=cantbeamabhc3_4x12x12 user=abhc3-std.obj interactive
Old job files exist. Overwrite? (y/n): y
ABAQUS JOB cantbeamabhc3_4x12x12
ABAQUS Version 6.6-1
Begin Linking ABAQUS/Standard User Subroutines
6/5/2007 10:31:03 AM
Creating library standardU.lib and object standardU.exp
End Linking ABAQUS/Standard User Subroutines
6/5/2007 10:31:03 AM
Begin Analysis Input File Processor
6/5/2007 10:31:03 AM
Run pre.exe
ABAQUS License Manager checked out the following licenses:
ABAQUS/Standard checked out 5 tokens.
<13 out of 18 licenses remain available>.
6/5/2007 10:31:07 AM
End Analysis Input File Processor
Begin ABAQUS/Standard Analysis
6/5/2007 10:31:07 AM
Run standard.exe
ABAQUS License Manager checked out the following licenses:
ABAQUS/Standard checked out 5 tokens.
<13 out of 18 licenses remain available>.
6/5/2007 10:31:10 AM
ABAQUS Error: The executable C:\ABAQUS\6.6-1\exec\standard.exe
aborted with system error code 142.
Please check the .dat, .msg, and .sta files for error messages if the files
exist. If there are no error messages and you cannot resolve the problem,
please run the command "abaqus job=support information=support" to report and
save your system information. Use the same command to run ABAQUS that you
used when the problem occurred. Please contact your local ABAQUS support
office and send them the input file, the file support.log which you just
created, the executable name, and the error code.
ABAQUS/Analysis exited with errors
[\quote]
Note - In cantBeamAbHC3_4x12x12.dat and cantBeamAbHC3_4x12x12.msg there are no warning or error messages.Abaqus does not generate cantBeamAbHC3_4x12x12.sta
I try to run an abaqus job analysis with a uel subroutine written in c++ and I get
[quote]
ABAQUS Error: The executable C:\ABAQUS\6.6-1\exec\standard.exe
aborted with system error code 142.
[\quote]
So I try to run it with a version of this uel subroutine written in fortran and I get the same error message.Has any one any idea to solve my problem? Thanks a lot in advance.
Here is my input file (cantBeamAbHC3_4x12x12.inp):
[code]
*HEADING
*PREPRINT, MODEL=YES
*NODE
1001001, 0.00000000, 0.00000000, 0.00000000
14001001, 0.00000000, 0.00000000, 30.00000000
1014001, 0.00000000, 30.00000000, 0.00000000
14014001, 0.00000000, 30.00000000, 30.00000000
1001006, 10.00000000, 0.00000000, 0.00000000
14001006, 10.00000000, 0.00000000, 30.00000000
1014006, 10.00000000, 30.00000000, 0.00000000
14014006, 10.00000000, 30.00000000, 30.00000000
1001002, 1.25000000, 0.00000000, 0.00000000
1001005, 8.75000000, 0.00000000, 0.00000000
1002001, 0.00000000, 1.25000000, 0.00000000
1013001, 0.00000000, 28.75000000, 0.00000000
*NGEN,NSET=LX1
1001002,1001005
*NGEN,NSET=LY1
1002001,1013001,1000
*NCOPY,OLDSET=LX1,NEWSET=LX2,CHANGENUMBER=13000,SHIFT
0.00000000, 30.00000000, 0.00000000
*NCOPY,OLDSET=LY1,NEWSET=LY2,CHANGENUMBER=5,SHIFT
10.00000000, 0.00000000, 0.00000000
*NCOPY,OLDSET=LX1,NEWSET=NSETA,CHANGENUMBER=1000,SHIFT
0.00000000, 1.25000000, 0.00000000
*NCOPY,OLDSET=LX1,NEWSET=NSETB,CHANGENUMBER=12000,SHIFT
0.00000000, 28.75000000, 0.00000000
*NFILL,NSET=SuperfInf
NSETA,NSETB,11,1000
*NSET,NSET=SuperfInf
LX1, LX2, LY1, LY2,1001001,1014001,1001006,1014006
*NCOPY,OLDSET=SuperfInf,CHANGENUMBER=13000000,SHIFT
0.00000000, 0.00000000, 30.00000000
*NCOPY,OLDSET=SuperfInf,NEWSET=NSETAA,CHANGENUMBER=1000000,SHIFT
0.00000000, 0.00000000, 1.25000000
*NCOPY,OLDSET=SuperfInf,NEWSET=NSETBB,CHANGENUMBER=12000000,SHIFT
0.00000000, 0.00000000, 28.75000000
*NFILL
NSETAA,NSETBB,11,1000000
**
*USER ELEMENT, NODES=27, TYPE=U111, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U111, ELSET=E111
1001001, 1001001, 2001001, 3001001, 1002001, 2002001, 3002001, 1003001, 2003001,
3003001, 1001002, 2001002, 3001002, 1002002, 2002002, 3002002, 1003002,
2003002, 3003002, 1001003, 2001003, 3001003, 1002003, 2002003, 3002003,
1003003, 2003003, 3003003
*UEL PROPERTY, ELSET=E111
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U112, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U112, ELSET=E112
2001001, 2001001, 3001001, 4001001, 2002001, 3002001, 4002001, 2003001, 3003001,
4003001, 2001002, 3001002, 4001002, 2002002, 3002002, 4002002, 2003002,
3003002, 4003002, 2001003, 3001003, 4001003, 2002003, 3002003, 4002003,
2003003, 3003003, 4003003
*ELGEN,ELSET=E112
2001001, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E112
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U113, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U113, ELSET=E113
12001001, 12001001, 13001001, 14001001, 12002001, 13002001, 14002001, 12003001, 13003001,
14003001, 12001002, 13001002, 14001002, 12002002, 13002002, 14002002, 12003002,
13003002, 14003002, 12001003, 13001003, 14001003, 12002003, 13002003, 14002003,
12003003, 13003003, 14003003
*UEL PROPERTY, ELSET=E113
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U121, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U121, ELSET=E121
1002001, 1002001, 2002001, 3002001, 1003001, 2003001, 3003001, 1004001, 2004001,
3004001, 1002002, 2002002, 3002002, 1003002, 2003002, 3003002, 1004002,
2004002, 3004002, 1002003, 2002003, 3002003, 1003003, 2003003, 3003003,
1004003, 2004003, 3004003
*ELGEN,ELSET=E121
1002001, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E121
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U122, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U122, ELSET=E122
2002001, 2002001, 3002001, 4002001, 2003001, 3003001, 4003001, 2004001, 3004001,
4004001, 2002002, 3002002, 4002002, 2003002, 3003002, 4003002, 2004002,
3004002, 4004002, 2002003, 3002003, 4002003, 2003003, 3003003, 4003003,
2004003, 3004003, 4004003
*ELGEN,ELSET=E122
2002001, 1,1,1, 10,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E122
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U123, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U123, ELSET=E123
12002001, 12002001, 13002001, 14002001, 12003001, 13003001, 14003001, 12004001, 13004001,
14004001, 12002002, 13002002, 14002002, 12003002, 13003002, 14003002, 12004002,
13004002, 14004002, 12002003, 13002003, 14002003, 12003003, 13003003, 14003003,
12004003, 13004003, 14004003
*ELGEN,ELSET=E123
12002001, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E123
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U131, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U131, ELSET=E131
1012001, 1012001, 2012001, 3012001, 1013001, 2013001, 3013001, 1014001, 2014001,
3014001, 1012002, 2012002, 3012002, 1013002, 2013002, 3013002, 1014002,
2014002, 3014002, 1012003, 2012003, 3012003, 1013003, 2013003, 3013003,
1014003, 2014003, 3014003
*UEL PROPERTY, ELSET=E131
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U132, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U132, ELSET=E132
2012001, 2012001, 3012001, 4012001, 2013001, 3013001, 4013001, 2014001, 3014001,
4014001, 2012002, 3012002, 4012002, 2013002, 3013002, 4013002, 2014002,
3014002, 4014002, 2012003, 3012003, 4012003, 2013003, 3013003, 4013003,
2014003, 3014003, 4014003
*ELGEN,ELSET=E132
2012001, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E132
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U133, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U133, ELSET=E133
12012001, 12012001, 13012001, 14012001, 12013001, 13013001, 14013001, 12014001, 13014001,
14014001, 12012002, 13012002, 14012002, 12013002, 13013002, 14013002, 12014002,
13014002, 14014002, 12012003, 13012003, 14012003, 12013003, 13013003, 14013003,
12014003, 13014003, 14014003
*UEL PROPERTY, ELSET=E133
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U211, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U211, ELSET=E211
1001002, 1001002, 2001002, 3001002, 1002002, 2002002, 3002002, 1003002, 2003002,
3003002, 1001003, 2001003, 3001003, 1002003, 2002003, 3002003, 1003003,
2003003, 3003003, 1001004, 2001004, 3001004, 1002004, 2002004, 3002004,
1003004, 2003004, 3003004
*ELGEN,ELSET=E211
1001002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E211
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U212, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U212, ELSET=E212
2001002, 2001002, 3001002, 4001002, 2002002, 3002002, 4002002, 2003002, 3003002,
4003002, 2001003, 3001003, 4001003, 2002003, 3002003, 4002003, 2003003,
3003003, 4003003, 2001004, 3001004, 4001004, 2002004, 3002004, 4002004,
2003004, 3003004, 4003004
*ELGEN,ELSET=E212
2001002, 2,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E212
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U213, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U213, ELSET=E213
12001002, 12001002, 13001002, 14001002, 12002002, 13002002, 14002002, 12003002, 13003002,
14003002, 12001003, 13001003, 14001003, 12002003, 13002003, 14002003, 12003003,
13003003, 14003003, 12001004, 13001004, 14001004, 12002004, 13002004, 14002004,
12003004, 13003004, 14003004
*ELGEN,ELSET=E213
12001002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E213
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U221, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U221, ELSET=E221
1002002, 1002002, 2002002, 3002002, 1003002, 2003002, 3003002, 1004002, 2004002,
3004002, 1002003, 2002003, 3002003, 1003003, 2003003, 3003003, 1004003,
2004003, 3004003, 1002004, 2002004, 3002004, 1003004, 2003004, 3003004,
1004004, 2004004, 3004004
*ELGEN,ELSET=E221
1002002, 2,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E221
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U222, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U222, ELSET=E222
2002002, 2002002, 3002002, 4002002, 2003002, 3003002, 4003002, 2004002, 3004002,
4004002, 2002003, 3002003, 4002003, 2003003, 3003003, 4003003, 2004003,
3004003, 4004003, 2002004, 3002004, 4002004, 2003004, 3003004, 4003004,
2004004, 3004004, 4004004
*ELGEN,ELSET=E222
2002002, 2,1,1, 10,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E222
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U223, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U223, ELSET=E223
12002002, 12002002, 13002002, 14002002, 12003002, 13003002, 14003002, 12004002, 13004002,
14004002, 12002003, 13002003, 14002003, 12003003, 13003003, 14003003, 12004003,
13004003, 14004003, 12002004, 13002004, 14002004, 12003004, 13003004, 14003004,
12004004, 13004004, 14004004
*ELGEN,ELSET=E223
12002002, 2,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E223
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U231, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U231, ELSET=E231
1012002, 1012002, 2012002, 3012002, 1013002, 2013002, 3013002, 1014002, 2014002,
3014002, 1012003, 2012003, 3012003, 1013003, 2013003, 3013003, 1014003,
2014003, 3014003, 1012004, 2012004, 3012004, 1013004, 2013004, 3013004,
1014004, 2014004, 3014004
*ELGEN,ELSET=E231
1012002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E231
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U232, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U232, ELSET=E232
2012002, 2012002, 3012002, 4012002, 2013002, 3013002, 4013002, 2014002, 3014002,
4014002, 2012003, 3012003, 4012003, 2013003, 3013003, 4013003, 2014003,
3014003, 4014003, 2012004, 3012004, 4012004, 2013004, 3013004, 4013004,
2014004, 3014004, 4014004
*ELGEN,ELSET=E232
2012002, 2,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E232
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U233, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U233, ELSET=E233
12012002, 12012002, 13012002, 14012002, 12013002, 13013002, 14013002, 12014002, 13014002,
14014002, 12012003, 13012003, 14012003, 12013003, 13013003, 14013003, 12014003,
13014003, 14014003, 12012004, 13012004, 14012004, 12013004, 13013004, 14013004,
12014004, 13014004, 14014004
*ELGEN,ELSET=E233
12012002, 2,1,1, 1,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E233
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U311, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U311, ELSET=E311
1001004, 1001004, 2001004, 3001004, 1002004, 2002004, 3002004, 1003004, 2003004,
3003004, 1001005, 2001005, 3001005, 1002005, 2002005, 3002005, 1003005,
2003005, 3003005, 1001006, 2001006, 3001006, 1002006, 2002006, 3002006,
1003006, 2003006, 3003006
*UEL PROPERTY, ELSET=E311
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U312, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U312, ELSET=E312
2001004, 2001004, 3001004, 4001004, 2002004, 3002004, 4002004, 2003004, 3003004,
4003004, 2001005, 3001005, 4001005, 2002005, 3002005, 4002005, 2003005,
3003005, 4003005, 2001006, 3001006, 4001006, 2002006, 3002006, 4002006,
2003006, 3003006, 4003006
*ELGEN,ELSET=E312
2001004, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E312
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U313, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U313, ELSET=E313
12001004, 12001004, 13001004, 14001004, 12002004, 13002004, 14002004, 12003004, 13003004,
14003004, 12001005, 13001005, 14001005, 12002005, 13002005, 14002005, 12003005,
13003005, 14003005, 12001006, 13001006, 14001006, 12002006, 13002006, 14002006,
12003006, 13003006, 14003006
*UEL PROPERTY, ELSET=E313
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U321, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U321, ELSET=E321
1002004, 1002004, 2002004, 3002004, 1003004, 2003004, 3003004, 1004004, 2004004,
3004004, 1002005, 2002005, 3002005, 1003005, 2003005, 3003005, 1004005,
2004005, 3004005, 1002006, 2002006, 3002006, 1003006, 2003006, 3003006,
1004006, 2004006, 3004006
*ELGEN,ELSET=E321
1002004, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E321
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U322, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U322, ELSET=E322
2002004, 2002004, 3002004, 4002004, 2003004, 3003004, 4003004, 2004004, 3004004,
4004004, 2002005, 3002005, 4002005, 2003005, 3003005, 4003005, 2004005,
3004005, 4004005, 2002006, 3002006, 4002006, 2003006, 3003006, 4003006,
2004006, 3004006, 4004006
*ELGEN,ELSET=E322
2002004, 1,1,1, 10,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E322
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U323, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U323, ELSET=E323
12002004, 12002004, 13002004, 14002004, 12003004, 13003004, 14003004, 12004004, 13004004,
14004004, 12002005, 13002005, 14002005, 12003005, 13003005, 14003005, 12004005,
13004005, 14004005, 12002006, 13002006, 14002006, 12003006, 13003006, 14003006,
12004006, 13004006, 14004006
*ELGEN,ELSET=E323
12002004, 1,1,1, 10,1000,1000, 1,1000000,1000000
*UEL PROPERTY, ELSET=E323
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U331, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U331, ELSET=E331
1012004, 1012004, 2012004, 3012004, 1013004, 2013004, 3013004, 1014004, 2014004,
3014004, 1012005, 2012005, 3012005, 1013005, 2013005, 3013005, 1014005,
2014005, 3014005, 1012006, 2012006, 3012006, 1013006, 2013006, 3013006,
1014006, 2014006, 3014006
*UEL PROPERTY, ELSET=E331
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U332, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U332, ELSET=E332
2012004, 2012004, 3012004, 4012004, 2013004, 3013004, 4013004, 2014004, 3014004,
4014004, 2012005, 3012005, 4012005, 2013005, 3013005, 4013005, 2014005,
3014005, 4014005, 2012006, 3012006, 4012006, 2013006, 3013006, 4013006,
2014006, 3014006, 4014006
*ELGEN,ELSET=E332
2012004, 1,1,1, 1,1000,1000, 10,1000000,1000000
*UEL PROPERTY, ELSET=E332
200.000000,0.300000,2.500000,2.500000,2.500000
**
*USER ELEMENT, NODES=27, TYPE=U333, PROPERTIES=5, VARIABLES=81, COORDINATES=3
1
1,2
1,3
*ELEMENT,TYPE=U333, ELSET=E333
12012004, 12012004, 13012004, 14012004, 12013004, 13013004, 14013004, 12014004, 13014004,
14014004, 12012005, 13012005, 14012005, 12013005, 13013005, 14013005, 12014005,
13014005, 14014005, 12012006, 13012006, 14012006, 12013006, 13013006, 14013006,
12014006, 13014006, 14014006
*UEL PROPERTY, ELSET=E333
200.000000,0.300000,2.500000,2.500000,2.500000
**
*BOUNDARY
SuperfInf,PINNED
**
*STEP,PERT
*STATIC
*INCLUDE, INPUT=Cloads_4x12x12.inp
*OUTPUT,FIELD
*NODE OUTPUT
UT,COORD
*NODE PRINT,F=0
*ENDSTEP
[\code]
Here is my fortran uel subroutine (AbHC3.for):
[code]
SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,
1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,
2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,
3 NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,
4 PERIOD)
C
INCLUDE 'ABA_PARAM.INC'
PARAMETER ( ZERO = 0.D0, ONE = 1.D0 , HALF = 0.5D0, TWO = 2.D0)
C
DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),
1 SVARS(NSVARS),ENERGY(8),PROPS(*),COORDS(MCRD,NNODE),
2 U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2),
3 PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),
4 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*),
5 JPROPS(*)
DIMENSION H(3), A00(3,3,3), A10(3,3,3), A01(3,3,3), A11(3,3,3),
1 XI(3)
C
INTEGER N, NR, NC, DR, DC
C
C UEL SUBROUTINE FOR A HORIZONTAL TRUSS ELEMENT
C
C SRESID - stores the static residual at time t+dt
C SVARS - In 1-6, contains the static residual at time t
C upon entering the routine. SRESID is copied to
C SVARS(1-6) after the dynamic residual has been
C calculated.
C - For half-step residual calculations: In 7-12,
C contains the static residual at the beginning
C of the previous increment. SVARS(1-6) are copied
C into SVARS(7-12) after the dynamic residual has
C been calculated.
C
YOUNG = PROPS(1)
POISSON = PROPS(2)
H(1) = PROPS(3)
H(2) = PROPS(4)
H(3) = PROPS(5)
C
LAMBDA = YOUNG*POISSON/(ONE+POISSON)/(ONE-TWO*POISSON)
MU = YOUNG/TWO/(ONE+POISSON)
C
A00(1,1,1) = 0.2000000000D0
A00(1,1,2) = 0.1166666667D0
A00(1,1,3) = 0.1666666667D-1
A00(1,2,1) = 0.1166666667D0
A00(1,2,2) = 0.2833333333D0
A00(1,2,3) = 0.1000000000D0
A00(1,3,1) = 0.1666666667D-1
A00(1,3,2) = 0.1000000000D0
A00(1,3,3) = 0.5000000000D-1
A00(2,1,1) = 0.5000000000D-1
A00(2,1,2) = 0.1083333333D0
A00(2,1,3) = 0.8333333333D-2
A00(2,2,1) = 0.1083333333D0
A00(2,2,2) = 0.4500000000D0
A00(2,2,3) = 0.1083333333D0
A00(2,3,1) = 0.8333333333D-2
A00(2,3,2) = 0.1083333333D0
A00(2,3,3) = 0.5000000000D-1
A00(3,1,1) = 0.5000000000D-1
A00(3,1,2) = 0.1000000000D0
A00(3,1,3) = 0.1666666667D-1
A00(3,2,1) = 0.1000000000D0
A00(3,2,2) = 0.2833333333D0
A00(3,2,3) = 0.1166666667D0
A00(3,3,1) = 0.1666666667D-1
A00(3,3,2) = 0.1166666667D0
A00(3,3,3) = 0.2000000000D0
A01(1,1,1) = -0.5000000000D0
A01(1,1,2) = 0.4166666667D0
A01(1,1,3) = 0.8333333333D-1
A01(1,2,1) = -0.4166666667D0
A01(1,2,2) = 0.1250000000D0
A01(1,2,3) = 0.2916666667D0
A01(1,3,1) = -0.8333333333D-1
A01(1,3,2) = -0.4166666667D-1
A01(1,3,3) = 0.1250000000D0
A01(2,1,1) = -0.1250000000D0
A01(2,1,2) = 0.8333333333D-1
A01(2,1,3) = 0.4166666667D-1
A01(2,2,1) = -0.3333333333D0
A01(2,2,2) = 0.0D0
A01(2,2,3) = 0.3333333333D0
A01(2,3,1) = -0.4166666667D-1
A01(2,3,2) = -0.8333333333D-1
A01(2,3,3) = 0.1250000000D0
A01(3,1,1) = -0.1250000000D0
A01(3,1,2) = 0.4166666667D-1
A01(3,1,3) = 0.8333333333D-1
A01(3,2,1) = -0.2916666667D0
A01(3,2,2) = -0.1250000000D0
A01(3,2,3) = 0.4166666667D0
A01(3,3,1) = -0.8333333333D-1
A01(3,3,2) = -0.4166666667D0
A01(3,3,3) = 0.5000000000D0
A10(1,1,1) = -0.5000000000D0
A10(1,1,2) = -0.4166666667D0
A10(1,1,3) = -0.8333333333D-1
A10(1,2,1) = 0.4166666667D0
A10(1,2,2) = 0.1250000000D0
A10(1,2,3) = -0.4166666667D-1
A10(1,3,1) = 0.8333333333D-1
A10(1,3,2) = 0.2916666667D0
A10(1,3,3) = 0.1250000000D0
A10(2,1,1) = -0.1250000000D0
A10(2,1,2) = -0.3333333333D0
A10(2,1,3) = -0.4166666667D-1
A10(2,2,1) = 0.8333333333D-1
A10(2,2,2) = 0.0D0
A10(2,2,3) = -0.8333333333D-1
A10(2,3,1) = 0.4166666667D-1
A10(2,3,2) = 0.3333333333D0
A10(2,3,3) = 0.1250000000D0
A10(3,1,1) = -0.1250000000D0
A10(3,1,2) = -0.2916666667D0
A10(3,1,3) = -0.8333333333D-1
A10(3,2,1) = 0.4166666667D-1
A10(3,2,2) = -0.1250000000D0
A10(3,2,3) = -0.4166666667D0
A10(3,3,1) = 0.8333333333D-1
A10(3,3,2) = 0.4166666667D0
A10(3,3,3) = 0.5000000000D0
A11(1,1,1) = 0.1333333333D1
A11(1,1,2) = -0.1000000000D1
A11(1,1,3) = -0.3333333333D0
A11(1,2,1) = -0.1000000000D1
A11(1,2,2) = 0.1D1
A11(1,2,3) = -0.6938893904D-17
A11(1,3,1) = -0.3333333333D0
A11(1,3,2) = -0.6938893904D-17
A11(1,3,3) = 0.3333333333D0
A11(2,1,1) = 0.3333333333D0
A11(2,1,2) = -0.1666666667D0
A11(2,1,3) = -0.1666666667D0
A11(2,2,1) = -0.1666666667D0
A11(2,2,2) = 0.3333333333D0
A11(2,2,3) = -0.1666666667D0
A11(2,3,1) = -0.1666666667D0
A11(2,3,2) = -0.1666666667D0
A11(2,3,3) = 0.3333333333D0
A11(3,1,1) = 0.3333333333D0
A11(3,1,2) = -0.6938893904D-17
A11(3,1,3) = -0.3333333333D0
A11(3,2,1) = -0.6938893904D-17
A11(3,2,2) = 0.1D1
A11(3,2,3) = -0.1000000000D1
A11(3,3,1) = -0.3333333333D0
A11(3,3,2) = -0.1000000000D1
A11(3,3,3) = 0.1333333333D1
C
IF(JTYPE.EQ.111)THEN
XI(1) = 1
XI(2) = 1
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.112)THEN
XI(1) = 1
XI(2) = 1
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.113)THEN
XI(1) = 1
XI(2) = 1
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.121)THEN
XI(1) = 1
XI(2) = 2
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.122)THEN
XI(1) = 1
XI(2) = 2
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.123)THEN
XI(1) = 1
XI(2) = 2
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.131)THEN
XI(1) = 1
XI(2) = 3
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.132)THEN
XI(1) = 1
XI(2) = 3
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.133)THEN
XI(1) = 1
XI(2) = 3
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.211)THEN
XI(1) = 2
XI(2) = 1
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.212)THEN
XI(1) = 2
XI(2) = 1
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.213)THEN
XI(1) = 2
XI(2) = 1
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.221)THEN
XI(1) = 2
XI(2) = 2
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.222)THEN
XI(1) = 2
XI(2) = 2
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.223)THEN
XI(1) = 2
XI(2) = 2
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.231)THEN
XI(1) = 2
XI(2) = 3
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.232)THEN
XI(1) = 2
XI(2) = 3
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.233)THEN
XI(1) = 2
XI(2) = 3
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.311)THEN
XI(1) = 3
XI(2) = 1
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.312)THEN
XI(1) = 3
XI(2) = 1
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.313)THEN
XI(1) = 3
XI(2) = 1
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.321)THEN
XI(1) = 3
XI(2) = 2
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.322)THEN
XI(1) = 3
XI(2) = 2
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.323)THEN
XI(1) = 3
XI(2) = 2
XI(3) = 3
C GOTO 999
ELSE IF(JTYPE.EQ.331)THEN
XI(1) = 3
XI(2) = 3
XI(3) = 1
C GOTO 999
ELSE IF(JTYPE.EQ.332)THEN
XI(1) = 3
XI(2) = 3
XI(3) = 2
C GOTO 999
ELSE IF(JTYPE.EQ.333)THEN
XI(1) = 3
XI(2) = 3
XI(3) = 3
C GOTO 999
END IF
C
C1111 = lambda + TWO*mu
C2222 = lambda + TWO*mu
C3333 = lambda + TWO*mu
C1122 = lambda
C1133 = lambda
C2233 = lambda
C1212 = mu
C1313 = mu
C2323 = mu
H123 = H(0)*H(1)/H(2)
H231 = H(1)*H(2)/H(0)
H312 = H(2)*H(0)/H(1)
N = 3
C
IF (LFLAGS(3).EQ.1) THEN
C Normal incrementation
IF (LFLAGS(1).EQ.1 .OR. LFLAGS(1).EQ.2) THEN
C *STATIC
IF (LFLAGS(4).NE.0) THEN
DO I = 1, N
DO J = 1, N
DO K = 1, N
DO P = 1, N
DO Q = 1, N
DO R = 1, N
NR = (I-1)*N*N + (J-1)*N + K
NC = (P-1)*N*N + (Q-1)*N + R
C
C contributo u1ijk u1pqr
C
DR = 1
DC = 1
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H231*C1111*A11(XI(0),I,P)*A00(XI(1),J,Q)*A00(XI(2),K,R) +
* H312*C1212*A00(XI(0),I,P)*A11(XI(1),J,Q)*A00(XI(2),K,R) +
* H123*C1313*A00(XI(0),I,P)*A00(XI(1),J,Q)*A11(XI(2),K,R)
C contributo u1ijk u2pqr
DR = 1
DC = 2
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(2)*C1122*A10(XI(0),I,P)*A01(XI(1),J,Q)*A00(XI(2),K,R) +
* H(2)*C1212*A01(XI(0),I,P)*A10(XI(1),J,Q)*A00(XI(2),K,R)
C contributo u1ijk u3pqr
DR = 1
DC = 3
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(1)*C1133*A10(XI(0),I,P)*A00(XI(1),J,Q)*A01(XI(2),K,R) +
* H(1)*C1313*A01(XI(0),I,P)*A00(XI(1),J,Q)*A10(XI(2),K,R)
C
C contributo u2ijk u2pqr
C
DR = 2
DC = 2
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H231*C1212*A11(XI(0),I,P)*A00(XI(1),J,Q)*A00(XI(2),K,R) +
* H312*C2222*A00(XI(0),I,P)*A11(XI(1),J,Q)*A00(XI(2),K,R) +
* H123*C2323*A00(XI(0),I,P)*A00(XI(1),J,Q)*A11(XI(2),K,R)
C contributo u2ijk u1pqr
DR = 2
DC = 1
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(2)*C1122*A01(XI(0),I,P)*A10(XI(1),J,Q)*A00(XI(2),K,R) +
* H(2)*C1313*A10(XI(0),I,P)*A01(XI(1),J,Q)*A00(XI(2),K,R)
C contributo u2ijk u3pqr
DR = 2
DC = 3
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(0)*C2233*A00(XI(0),I,P)*A10(XI(1),J,Q)*A01(XI(2),K,R) +
* H(0)*C2323*A00(XI(0),I,P)*A01(XI(1),J,Q)*A10(XI(2),K,R)
C
C contributo u3ijk u3pqr
C
DR = 3
DC = 3
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H231*C1313*A11(XI(0),I,P)*A00(XI(1),J,Q)*A00(XI(2),K,R) +
* H312*C2323*A00(XI(0),I,P)*A11(XI(1),J,Q)*A00(XI(2),K,R) +
* H123*C3333*A00(XI(0),I,P)*A00(XI(1),J,Q)*A11(XI(2),K,R)
C contributo u3ijk u1pqr
DR = 3
DC = 1
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(1)*C1133*A01(XI(0),I,P)*A00(XI(1),J,Q)*A10(XI(2),K,R) +
* H(1)*C1313*A10(XI(0),I,P)*A00(XI(1),J,Q)*A01(XI(2),K,R)
C contributo u3ijk u2pqr
DR = 2
DC = 2
AMATRX( NR+(DR-1)*N*N*N , NC+(DC-1)*N*N*N ) =
* H(0)*C2233*A00(XI(0),I,P)*A01(XI(1),J,Q)*A10(XI(2),K,R) +
* H(0)*C2323*A00(XI(0),I,P)*A10(XI(1),J,Q)*A01(XI(2),K,R)
END DO
END DO
END DO
END DO
END DO
END DO
C
DO I = 1, NDOFEL
RHS(I,1) = ZERO
DO J = 1, NDOFEL
RHS(I,1) = RHS(I,1) - AMATRX(I,J)*DU(J,1)
END DO
END DO
C
DO I = 1, NDOFEL
ENERGY(2) = ZERO
DO J = 1, NDOFEL
ENERGY(2) = ENERGY(2) + HALF*DU(I,1)*AMATRX(I,J)*U(J)
* + HALF*U(I)*AMATRX(I,J)*DU(J,1)
* + HALF*DU(I,1)*AMATRX(I,J)*DU(J,1)
END DO
END DO
C
END IF
END IF
ELSE IF(LFLAGS(3).EQ.100)THEN
C Output for perturbations
IF (LFLAGS(1).EQ.1 .OR. LFLAGS(1).EQ.2) THEN
C *STATIC
DO I = 1, NDOFEL
RHS(I,1) = ZERO
DO J = 1, NDOFEL
RHS(I,1) = RHS(I,1) - AMATRX(I,J)*DU(J,1)
END DO
END DO
C
DO I = 1, NDOFEL
ENERGY(2) = ZERO
DO J = 1, NDOFEL
ENERGY(2) = ENERGY(2) + HALF*DU(I,1)*AMATRX(I,J)*U(J)
* + HALF*U(I)*AMATRX(I,J)*DU(J,1)
* + HALF*DU(I,1)*AMATRX(I,J)*DU(J,1)
END DO
END DO
C
DO KVAR = 1, NSVARS
SVARS(KVAR) = RHS(KVAR,1)
END DO
C
END IF
END IF
C
RETURN
END
[\code]
Here is what I get in abaqus command window:
[quote]
C:\Documents and Settings\alorenzom\Documenti\HC3FEM\CantBeamAbHC3withFortSr>aba
qus make library=abhc3
ABAQUS JOB abhc3
Begin Compiling ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
End Compiling ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
Begin Linking ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
Creating library standardU.lib and object standardU.exp
End Linking ABAQUS/Standard User Subroutines
6/5/2007 10:29:36 AM
ABAQUS JOB abhc3 COMPLETED
C:\Documents and Settings\alorenzom\Documenti\HC3FEM\CantBeamAbHC3withFortSr>aba
qus job=cantbeamabhc3_4x12x12 user=abhc3-std.obj interactive
Old job files exist. Overwrite? (y/n): y
ABAQUS JOB cantbeamabhc3_4x12x12
ABAQUS Version 6.6-1
Begin Linking ABAQUS/Standard User Subroutines
6/5/2007 10:31:03 AM
Creating library standardU.lib and object standardU.exp
End Linking ABAQUS/Standard User Subroutines
6/5/2007 10:31:03 AM
Begin Analysis Input File Processor
6/5/2007 10:31:03 AM
Run pre.exe
ABAQUS License Manager checked out the following licenses:
ABAQUS/Standard checked out 5 tokens.
<13 out of 18 licenses remain available>.
6/5/2007 10:31:07 AM
End Analysis Input File Processor
Begin ABAQUS/Standard Analysis
6/5/2007 10:31:07 AM
Run standard.exe
ABAQUS License Manager checked out the following licenses:
ABAQUS/Standard checked out 5 tokens.
<13 out of 18 licenses remain available>.
6/5/2007 10:31:10 AM
ABAQUS Error: The executable C:\ABAQUS\6.6-1\exec\standard.exe
aborted with system error code 142.
Please check the .dat, .msg, and .sta files for error messages if the files
exist. If there are no error messages and you cannot resolve the problem,
please run the command "abaqus job=support information=support" to report and
save your system information. Use the same command to run ABAQUS that you
used when the problem occurred. Please contact your local ABAQUS support
office and send them the input file, the file support.log which you just
created, the executable name, and the error code.
ABAQUS/Analysis exited with errors
[\quote]
Note - In cantBeamAbHC3_4x12x12.dat and cantBeamAbHC3_4x12x12.msg there are no warning or error messages.Abaqus does not generate cantBeamAbHC3_4x12x12.sta





RE: analysis job with a uel subroutine
I solve my problem and now my abaqus job completed.
RE: analysis job with a uel subroutine
My ABAQUS JOB COMPLETED but in msg file there are WARNING stating Solver problem-Zero pivot when processing 49 of my 78 loaded dofs and in odb there are warnings stating zero pivot for all my 78 loaded dofs.And results aren't those I was expecting, in particular displacements are all zero-valued.
RE: analysis job with a uel subroutine
results.
That isn't the same for the version in c++.I don't know why but I
have zero pivot for all my loaded degrees of freedom and then
all displacement zero-valued. I have printed for each of my elements
amatrx and they are correct.I have printed also for each of my
elements rhs and they have all components zero-valued. The only
statements about I'm doubtful if there may be difference between the
version in c++ of my subroutine and the one in fortran are the
following:
CODE
FOR_NAME(uel, UEL)(double** pRhs,double** pAmatrx,double*
pSvars,double* pEnergy,int* pNdofel,int* pNrhs,int* pNsvars,
double* pProps,int* pNprops,double** pCoords,int* pMcrd,int*
pNnode,double* pU,double** pDu,double* pV,double* pA,int*
pJtype,double* pTime,double* pDtime,
int* pKstep,int* pKinc,int* pJelem,double* pParams,int* pNdload,int**
pJdltyp,double** pAdlmag,double*** pPredef,
int* pNpredf,int* pLflags,int* pMlvarx,double** pDdlmag,int*
pMdload,double* pPnewdt,int* pJprops,int* pNjprop,
double* pPeriod){
pAmatrx = new double*[*pNdofel];
for(int i=0; i<*pNdofel; i++)
pAmatrx[i] = new double[*pNdofel];
pRhs = new double*[*pMlvarx];
for(int i=0; i<*pMlvarx; i++)
pRhs[i] = new double[1];
pDu = new double*[*pMlvarx];
for(int i=0; i<*pMlvarx; i++)
pDu[i] = new double[1];
pU = new double[*pNdofel];
CODE
1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,
2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,
3 NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,
4 PERIOD)
C
DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),
1 SVARS(NSVARS),ENERGY(8),PROPS(*),COORDS(MCRD,NNODE),
2 U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2),
3 PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),
4 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*),
5 JPROPS(*)
Please if anyone has implemented abaqus c++ user subroutine
HELP ME!