analysis job with a uel subroutine

94710 · Jun 5, 2007

Dear all

I try to run an abaqus job analysis with a uel subroutine written in c++ and I get

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

94710 · Jun 6, 2007

Thanks to all of you anyway.
I solve my problem and now my abaqus job completed.

94710 · Jun 7, 2007

Because of my unexperience I thought that my problems finished but it isn't so.

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.

94710 · Jun 12, 2007

My uel subroutine written in fortran, now finally give me correct
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:

extern "C" void
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:

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
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!

Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

analysis job with a uel subroutine

94710

Structural

Recommended for you

94710

Structural

94710

Structural

94710

Structural

Similar threads

Part and Inventory Search

Sponsor