next up previous contents
Next: OUTPUT Up: Using the penalty method Previous: Axisymmetric solution, TUBE7C1

INPUT

...rmd2.........tube7c1.I1
;  NELEM NNOD ITED 0 NGD 0 KPER NPER 0  KSS
IP   10   53    7  0  3  0   0    0  0   2
;  CRIT SCALE THDEF   XL   XU    YL   YU   ZL  ZU  ALPHA
RP 1.01 0.001   1    999  2001 -101    1    0   0    0
XY
   C      1  0.100000000000E+04  0.000000000000E+00
   C      2  0.100000000000E+04 -0.500000000000E+02
   C      3  0.100000000000E+04 -0.100000000000E+03
   C      4  0.105000000000E+04  0.000000000000E+00
   C      5  0.110000000000E+04  0.000000000000E+00
   C      6  0.105000000000E+04 -0.100000000000E+03
   C      7  0.110000000000E+04 -0.100000000000E+03
   C      8  0.110000000000E+04 -0.500000000000E+02
   C      9  0.115000000000E+04  0.000000000000E+00
   C     10  0.120000000000E+04  0.000000000000E+00
   C     11  0.115000000000E+04 -0.100000000000E+03
   C     12  0.120000000000E+04 -0.100000000000E+03
   C     13  0.120000000000E+04 -0.500000000000E+02
   C     14  0.125000000000E+04  0.000000000000E+00
   C     15  0.130000000000E+04  0.000000000000E+00
   C     16  0.125000000000E+04 -0.100000000000E+03
   C     17  0.130000000000E+04 -0.100000000000E+03
   C     18  0.130000000000E+04 -0.500000000000E+02
   C     19  0.135000000000E+04  0.000000000000E+00
   C     20  0.140000000000E+04  0.000000000000E+00
   C     21  0.135000000000E+04 -0.100000000000E+03
   C     22  0.140000000000E+04 -0.100000000000E+03
   C     23  0.140000000000E+04 -0.500000000000E+02
   C     24  0.145000000000E+04  0.000000000000E+00
   C     25  0.150000000000E+04  0.000000000000E+00
   C     26  0.145000000000E+04 -0.100000000000E+03
   C     27  0.150000000000E+04 -0.100000000000E+03
   C     28  0.150000000000E+04 -0.500000000000E+02
   C     29  0.155000000000E+04  0.000000000000E+00
   C     30  0.160000000000E+04  0.000000000000E+00
   C     31  0.155000000000E+04 -0.100000000000E+03
   C     32  0.160000000000E+04 -0.100000000000E+03
   C     33  0.160000000000E+04 -0.500000000000E+02
   C     34  0.165000000000E+04  0.000000000000E+00
   C     35  0.170000000000E+04  0.000000000000E+00
   C     36  0.165000000000E+04 -0.100000000000E+03
   C     37  0.170000000000E+04 -0.100000000000E+03
   C     38  0.170000000000E+04 -0.500000000000E+02
   C     39  0.175000000000E+04  0.000000000000E+00
   C     40  0.180000000000E+04  0.000000000000E+00
   C     41  0.175000000000E+04 -0.100000000000E+03
   C     42  0.180000000000E+04 -0.100000000000E+03
   C     43  0.180000000000E+04 -0.500000000000E+02
   C     44  0.185000000000E+04  0.000000000000E+00
   C     45  0.190000000000E+04  0.000000000000E+00
   C     46  0.185000000000E+04 -0.100000000000E+03
   C     47  0.190000000000E+04 -0.100000000000E+03
   C     48  0.190000000000E+04 -0.500000000000E+02
   C     49  0.195000000000E+04  0.000000000000E+00
   C     50  0.200000000000E+04  0.000000000000E+00
   C     51  0.195000000000E+04 -0.100000000000E+03
   C     52  0.200000000000E+04 -0.100000000000E+03
   C     53  0.200000000000E+04 -0.500000000000E+02
EL T  7 E      1 N      1      3      7      5      2      6
                        8      4
EL T  7 E      2 N      5      7     12     10      8     11
                       13      9
EL T  7 E      3 N     10     12     17     15     13     16
                       18     14
EL T  7 E      4 N     15     17     22     20     18     21
                       23     19
EL T  7 E      5 N     20     22     27     25     23     26
                       28     24
EL T  7 E      6 N     25     27     32     30     28     31
                       33     29
EL T  7 E      7 N     30     32     37     35     33     36
                       38     34
EL T  7 E      8 N     35     37     42     40     38     41
                       43     39
EL T  7 E      9 N     40     42     47     45     43     46
                       48     44
EL T  7 E     10 N     45     47     52     50     48     51
                       53     49
EN
EN
...rpd2.........tube7c1.I2
IP 1

;     polynomial dependence on effective stress (IV=9)
MP 1 T -1 I 9 

;    E     alpha    nu   rho  yield stress kinematic hard.
  V 2E11  V 1E-5  V 0.3  V 0       V 0         V 0 
  
;       effective strain rate              
            V 0 0 0 2E-28
            
SV 1 T 6 V -100E6
;***************************************************
AS 1 /M 1
     /B 0 C 2 N 1:53  ; plane strain
     /R 0 600 600     ; stiffness matrix for T=600 C
AS 2
    /S 1 E 1 S 1      ; internal pressure
    /R 0 600 600      ; reference temp.= current temp.=T=600 C
EN
EN
...srh2.........tube7c1.I3
IP 1
EN
EN
...fefs.........tube7c1.I4
IP 1
EN
EN
...hpp2.........tube7c1.IP
  ;KREST NLC NCYC KMOD KCRP KLARG KCNT KTPR 2*0 LOAD CASE SEQUENCE
IP   1    2    0    0    1    0     0    0  2*0 2  2
  ;           10*0              CREEP TIMES IN HOURS
RP            10*0              0 10
EN
EN
...hpls.........tube7c1.IL
  ;KMET KPOUT NSUBI NITER NINT KTPR
IP   0    1     0     0    0    0
  ; UTOL  RTOL  XTOL
RP   0     0     0
EN
EN
...str2.........tube7c1.I5
   ;KLC   0   KOUT ILC   0  IEL1 IEU1 IEL2 IEU2   0  KPROB KGRAF
IP   3    0    1    2    0    0    0    0    0    0    2    1
EN
EN