关于creep子程序问题

zhoulan

写了一个蠕变子程序，不能收敛，望大佬解惑，仿真币太少了，见谅。              SUBROUTINE CREEP(DECRA,DESWA,STATEV,SERD,EC,ESW,P,QTILD,
     1TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME,CMNAME,LEXIMP,LEND,
     2COORDS,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME
C
      DIMENSION DECRA(5),DESWA(5),STATEV(*),PREDEF(*),DPRED(*),
     1TIME(3),EC(2),ESW(2),COORDS(*)
C
      A=6.2886E-6
      AN=6.0629
      C1=-13.68442+0.19957*TEMP-6.9875*0.0001*TEMP**2
      C2=-13.35202+0.11178*TEMP
      C3=37.23182-0.46776*TEMP+0.001488*TEMP**2
      C4=-6.07565+0.04292*TEMP
      C5=-49.85311+0.6351*TEMP-0.00198*TEMP**2
      Q=17439.84337-191.83477*TEMP+0.70838*TEMP**2
C
      DECRA(1)=A*QTILD**AN*(TIME(1)**C1+C2*TIME(1)**C3+C4*TIME(1)**C5-
     1         (TIME(1)-DTIME)**C1-C2*(TIME(1)-DTIME)**C3-
     2          C4*(TIME(1)-DTIME)**C5)*EXP(1.0)**(-Q/TEMP)
      write(*,*)DECRA(1)
      IF(LEXIMP.EQ.1) THEN
      DECRA(5)=A*AN*QTILD**(AN-1)*(TIME(1)**C1+C2*TIME(1)**C3+C4*TIME(1)
     1        **C5-(TIME(1)-DTIME)**C1-C2*(TIME(1)-DTIME)**C3-
     2        C4*(TIME(1)-DTIME)**C5)*EXP(1.0)**(-Q/TEMP)
      write(*,*)DECRA(5)
      END IF
C         
      RETURN
      END
下面为inp文件
*Heading
** Job name: CREEP Model name: Model-1
** Generated by: Abaqus/CAE 2016
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
      1,          10.,          4.5
      2,    8.9416666,          4.5
      3,   7.88333321,          4.5
      4,   6.82499981,          4.5
      5,   5.76666689,          4.5
      6,   4.70833349,          4.5
      7,    3.6500001,          4.5
      8,          10.,          3.5
      9,    8.9416666,          3.5
     10,   7.88333321,          3.5
     11,   6.82499981,          3.5
     12,   5.76666689,          3.5
     13,   4.70833349,          3.5
     14,    3.6500001,          3.5
     15,          10.,          2.5
     16,    8.9416666,          2.5
     17,   7.88333321,          2.5
     18,   6.82499981,          2.5
     19,   5.76666689,          2.5
     20,   4.70833349,          2.5
     21,    3.6500001,          2.5
     22,          10.,          1.5
     23,    8.9416666,          1.5
     24,   7.88333321,          1.5
     25,   6.82499981,          1.5
     26,   5.76666689,          1.5
     27,   4.70833349,          1.5
     28,    3.6500001,          1.5
     29,          10.,          0.5
     30,    8.9416666,          0.5
     31,   7.88333321,          0.5
     32,   6.82499981,          0.5
     33,   5.76666689,          0.5
     34,   4.70833349,          0.5
     35,    3.6500001,          0.5
     36,          10.,         -0.5
     37,    8.9416666,         -0.5
     38,   7.88333321,         -0.5
     39,   6.82499981,         -0.5
     40,   5.76666689,         -0.5
     41,   4.70833349,         -0.5
     42,    3.6500001,         -0.5
     43,          10.,         -1.5
     44,    8.9416666,         -1.5
     45,   7.88333321,         -1.5
     46,   6.82499981,         -1.5
     47,   5.76666689,         -1.5
     48,   4.70833349,         -1.5
     49,    3.6500001,         -1.5
     50,          10.,         -2.5
     51,    8.9416666,         -2.5
     52,   7.88333321,         -2.5
     53,   6.82499981,         -2.5
     54,   5.76666689,         -2.5
     55,   4.70833349,         -2.5
     56,    3.6500001,         -2.5
     57,          10.,         -3.5
     58,    8.9416666,         -3.5
     59,   7.88333321,         -3.5
     60,   6.82499981,         -3.5
     61,   5.76666689,         -3.5
     62,   4.70833349,         -3.5
     63,    3.6500001,         -3.5
     64,          10.,         -4.5
     65,    8.9416666,         -4.5
     66,   7.88333321,         -4.5
     67,   6.82499981,         -4.5
     68,   5.76666689,         -4.5
     69,   4.70833349,         -4.5
     70,    3.6500001,         -4.5
     71,          10.,         -5.5
     72,    8.9416666,         -5.5
     73,   7.88333321,         -5.5
     74,   6.82499981,         -5.5
     75,   5.76666689,         -5.5
     76,   4.70833349,         -5.5
     77,    3.6500001,         -5.5
     78,          10.,         -6.5
     79,    8.9416666,         -6.5
     80,   7.88333321,         -6.5
     81,   6.82499981,         -6.5
     82,   5.76666689,         -6.5
     83,   4.70833349,         -6.5
     84,    3.6500001,         -6.5
     85,          10.,         -7.5
     86,    8.9416666,         -7.5
     87,   7.88333321,         -7.5
     88,   6.82499981,         -7.5
     89,   5.76666689,         -7.5
     90,   4.70833349,         -7.5
     91,    3.6500001,         -7.5
     92,          10.,         -8.5
     93,    8.9416666,         -8.5
     94,   7.88333321,         -8.5
     95,   6.82499981,         -8.5
     96,   5.76666689,         -8.5
     97,   4.70833349,         -8.5
     98,    3.6500001,         -8.5
     99,          10.,         -9.5
    100,    8.9416666,         -9.5
    101,   7.88333321,         -9.5
    102,   6.82499981,         -9.5
    103,   5.76666689,         -9.5
    104,   4.70833349,         -9.5
    105,    3.6500001,         -9.5
    106,          10.,        -10.5
    107,    8.9416666,        -10.5
    108,   7.88333321,        -10.5
    109,   6.82499981,        -10.5
    110,   5.76666689,        -10.5
    111,   4.70833349,        -10.5
    112,    3.6500001,        -10.5
    113,          10.,        -11.5
    114,    8.9416666,        -11.5
    115,   7.88333321,        -11.5
    116,   6.82499981,        -11.5
    117,   5.76666689,        -11.5
    118,   4.70833349,        -11.5
    119,    3.6500001,        -11.5
    120,          10.,        -12.5
    121,    8.9416666,        -12.5
    122,   7.88333321,        -12.5
    123,   6.82499981,        -12.5
    124,   5.76666689,        -12.5
    125,   4.70833349,        -12.5
    126,    3.6500001,        -12.5
    127,          10.,        -13.5
    128,    8.9416666,        -13.5
    129,   7.88333321,        -13.5
    130,   6.82499981,        -13.5
    131,   5.76666689,        -13.5
    132,   4.70833349,        -13.5
    133,    3.6500001,        -13.5
    134,          10.,        -14.5
    135,    8.9416666,        -14.5
    136,   7.88333321,        -14.5
    137,   6.82499981,        -14.5
    138,   5.76666689,        -14.5
    139,   4.70833349,        -14.5
    140,    3.6500001,        -14.5
    141,          10.,        -15.5
    142,    8.9416666,        -15.5
    143,   7.88333321,        -15.5
    144,   6.82499981,        -15.5
    145,   5.76666689,        -15.5
    146,   4.70833349,        -15.5
    147,    3.6500001,        -15.5
    148,          10.,        -16.5
    149,    8.9416666,        -16.5
    150,   7.88333321,        -16.5
    151,   6.82499981,        -16.5
    152,   5.76666689,        -16.5
    153,   4.70833349,        -16.5
    154,    3.6500001,        -16.5
    155,          10.,        -17.5
    156,    8.9416666,        -17.5
    157,   7.88333321,        -17.5
    158,   6.82499981,        -17.5
    159,   5.76666689,        -17.5
    160,   4.70833349,        -17.5
    161,    3.6500001,        -17.5
    162,          10.,        -18.5
    163,    8.9416666,        -18.5
    164,   7.88333321,        -18.5
    165,   6.82499981,        -18.5
    166,   5.76666689,        -18.5
    167,   4.70833349,        -18.5
    168,    3.6500001,        -18.5
    169,          10.,        -19.5
    170,    8.9416666,        -19.5
    171,   7.88333321,        -19.5
    172,   6.82499981,        -19.5
    173,   5.76666689,        -19.5
    174,   4.70833349,        -19.5
    175,    3.6500001,        -19.5
    176,          10.,        -20.5
    177,    8.9416666,        -20.5
    178,   7.88333321,        -20.5
    179,   6.82499981,        -20.5
    180,   5.76666689,        -20.5
    181,   4.70833349,        -20.5
    182,    3.6500001,        -20.5
    183,          10.,        -21.5
    184,    8.9416666,        -21.5
    185,   7.88333321,        -21.5
    186,   6.82499981,        -21.5
    187,   5.76666689,        -21.5
    188,   4.70833349,        -21.5
    189,    3.6500001,        -21.5
    190,          10.,        -22.5
    191,    8.9416666,        -22.5
    192,   7.88333321,        -22.5
    193,   6.82499981,        -22.5
    194,   5.76666689,        -22.5
    195,   4.70833349,        -22.5
    196,    3.6500001,        -22.5
*Element, type=CPS4R
  1,   1,   2,   9,   8
  2,   2,   3,  10,   9
  3,   3,   4,  11,  10
  4,   4,   5,  12,  11
  5,   5,   6,  13,  12
  6,   6,   7,  14,  13
  7,   8,   9,  16,  15
  8,   9,  10,  17,  16
  9,  10,  11,  18,  17
10,  11,  12,  19,  18
11,  12,  13,  20,  19
12,  13,  14,  21,  20
13,  15,  16,  23,  22
14,  16,  17,  24,  23
15,  17,  18,  25,  24
16,  18,  19,  26,  25
17,  19,  20,  27,  26
18,  20,  21,  28,  27
19,  22,  23,  30,  29
20,  23,  24,  31,  30
21,  24,  25,  32,  31
22,  25,  26,  33,  32
23,  26,  27,  34,  33
24,  27,  28,  35,  34
25,  29,  30,  37,  36
26,  30,  31,  38,  37
27,  31,  32,  39,  38
28,  32,  33,  40,  39
29,  33,  34,  41,  40
30,  34,  35,  42,  41
31,  36,  37,  44,  43
32,  37,  38,  45,  44
33,  38,  39,  46,  45
34,  39,  40,  47,  46
35,  40,  41,  48,  47
36,  41,  42,  49,  48
37,  43,  44,  51,  50
38,  44,  45,  52,  51
39,  45,  46,  53,  52
40,  46,  47,  54,  53
41,  47,  48,  55,  54
42,  48,  49,  56,  55
43,  50,  51,  58,  57
44,  51,  52,  59,  58
45,  52,  53,  60,  59
46,  53,  54,  61,  60
47,  54,  55,  62,  61
48,  55,  56,  63,  62
49,  57,  58,  65,  64
50,  58,  59,  66,  65
51,  59,  60,  67,  66
52,  60,  61,  68,  67
53,  61,  62,  69,  68
54,  62,  63,  70,  69
55,  64,  65,  72,  71
56,  65,  66,  73,  72
57,  66,  67,  74,  73
58,  67,  68,  75,  74
59,  68,  69,  76,  75
60,  69,  70,  77,  76
61,  71,  72,  79,  78
62,  72,  73,  80,  79
63,  73,  74,  81,  80
64,  74,  75,  82,  81
65,  75,  76,  83,  82
66,  76,  77,  84,  83
67,  78,  79,  86,  85
68,  79,  80,  87,  86
69,  80,  81,  88,  87
70,  81,  82,  89,  88
71,  82,  83,  90,  89
72,  83,  84,  91,  90
73,  85,  86,  93,  92
74,  86,  87,  94,  93
75,  87,  88,  95,  94
76,  88,  89,  96,  95
77,  89,  90,  97,  96
78,  90,  91,  98,  97
79,  92,  93, 100,  99
80,  93,  94, 101, 100
81,  94,  95, 102, 101
82,  95,  96, 103, 102
83,  96,  97, 104, 103
84,  97,  98, 105, 104
85,  99, 100, 107, 106
86, 100, 101, 108, 107
87, 101, 102, 109, 108
88, 102, 103, 110, 109
89, 103, 104, 111, 110
90, 104, 105, 112, 111
91, 106, 107, 114, 113
92, 107, 108, 115, 114
93, 108, 109, 116, 115
94, 109, 110, 117, 116
95, 110, 111, 118, 117
96, 111, 112, 119, 118
97, 113, 114, 121, 120
98, 114, 115, 122, 121
99, 115, 116, 123, 122
100, 116, 117, 124, 123
101, 117, 118, 125, 124
102, 118, 119, 126, 125
103, 120, 121, 128, 127
104, 121, 122, 129, 128
105, 122, 123, 130, 129
106, 123, 124, 131, 130
107, 124, 125, 132, 131
108, 125, 126, 133, 132
109, 127, 128, 135, 134
110, 128, 129, 136, 135
111, 129, 130, 137, 136
112, 130, 131, 138, 137
113, 131, 132, 139, 138
114, 132, 133, 140, 139
115, 134, 135, 142, 141
116, 135, 136, 143, 142
117, 136, 137, 144, 143
118, 137, 138, 145, 144
119, 138, 139, 146, 145
120, 139, 140, 147, 146
121, 141, 142, 149, 148
122, 142, 143, 150, 149
123, 143, 144, 151, 150
124, 144, 145, 152, 151
125, 145, 146, 153, 152
126, 146, 147, 154, 153
127, 148, 149, 156, 155
128, 149, 150, 157, 156
129, 150, 151, 158, 157
130, 151, 152, 159, 158
131, 152, 153, 160, 159
132, 153, 154, 161, 160
133, 155, 156, 163, 162
134, 156, 157, 164, 163
135, 157, 158, 165, 164
136, 158, 159, 166, 165
137, 159, 160, 167, 166
138, 160, 161, 168, 167
139, 162, 163, 170, 169
140, 163, 164, 171, 170
141, 164, 165, 172, 171
142, 165, 166, 173, 172
143, 166, 167, 174, 173
144, 167, 168, 175, 174
145, 169, 170, 177, 176
146, 170, 171, 178, 177
147, 171, 172, 179, 178
148, 172, 173, 180, 179
149, 173, 174, 181, 180
150, 174, 175, 182, 181
151, 176, 177, 184, 183
152, 177, 178, 185, 184
153, 178, 179, 186, 185
154, 179, 180, 187, 186
155, 180, 181, 188, 187
156, 181, 182, 189, 188
157, 183, 184, 191, 190
158, 184, 185, 192, 191
159, 185, 186, 193, 192
160, 186, 187, 194, 193
161, 187, 188, 195, 194
162, 188, 189, 196, 195
*Nset, nset=_PickedSet2, internal, generate
   1,  196,    1
*Elset, elset=_PickedSet2, internal, generate
   1,  162,    1
** Section: Section-1
*Solid Section, elset=_PickedSet2, material=RUT300
,
*End Part
**  
**
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=Part-1-1, part=Part-1
*End Instance
**  
*Nset, nset=_PickedSet4, internal, instance=Part-1-1, generate
   1,  196,    1
*Elset, elset=_PickedSet4, internal, instance=Part-1-1, generate
   1,  162,    1
*Nset, nset=_PickedSet6, internal, instance=Part-1-1, generate
190,  196,    1
*Elset, elset=_PickedSet6, internal, instance=Part-1-1, generate
157,  162,    1
*Nset, nset=_PickedSet8, internal, instance=Part-1-1, generate
   1,  196,    1
*Elset, elset=_PickedSet8, internal, instance=Part-1-1, generate
   1,  162,    1
*Elset, elset=__PickedSurf7_S1, internal, instance=Part-1-1, generate
1,  6,  1
*Surface, type=ELEMENT, name=_PickedSurf7, internal
__PickedSurf7_S1, S1
*End Assembly
**
** MATERIALS
**
*Material, name=RUT300
*Conductivity
36.,
*Creep, law=USER
*Density
7.1e-09,
*Elastic
155000., 0.28
*Expansion
1.2e-05,
*Specific Heat
5.1e+08,
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: 对称/反对称/完全固定
*Boundary
_PickedSet6, ENCASTRE
** ----------------------------------------------------------------
**
** STEP: heat
**
*Step, name=heat, nlgeom=NO
*Static
1., 1., 1e-05, 1.
**
** LOADS
**
** Name: Load-1   Type: Pressure
*Dsload
_PickedSurf7, P, -200.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
U,
*Element Output, directions=YES
E, S
*Output, history, frequency=0
*End Step
** ----------------------------------------------------------------
**
** STEP: creep
**
*Step, name=creep, nlgeom=YES, inc=10000
*Visco, cetol=0.0001
0.1, 10., 1e-05, 10.
**
** PREDEFINED FIELDS
**
** Name: Predefined Field-3   Type: Temperature
*Temperature
_PickedSet8, 400.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-2
**
*Output, field
*Node Output
U,
*Element Output, directions=YES
CE, CEEQ, E, S
**
** FIELD OUTPUT: F-Output-1
**
*Node Output
U,
*Element Output, directions=YES
E, S
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step