The computation is executed with the following commands:
xrm2 beam4t2.i1 xrpd beam4t2.iB xt2t beam4t2
; NELEM NNOD ITED
IP 8 27 4
; CRIT SCALE THDEF
RP 1.01 1 0.01
XY N 1:10 X 2*0 2*0.25 2*0.5 2*0.75 2*1
Y 5*(0 0.02)
EL E 1:8 N =P 1 3 2 =P =Q 12 13 11 =Q =R 3 4 2 =R =S 15 14 13 =S
=2P =4Q =2R =4S =4P =8Q =4R =8S =6P =12Q =6R =12S
EN
EN
; KREST 0 KOUT INT3 NSAX NSTEPX KSU IP 1 0 2 0 10 1000 0 ; TIMS ERAL EDIF TOL DTRUN RP 0 1e-5 0.1 0.1 3600 MP 1 T 1 V 20 3e6 IV 1 T 5 V 0 0.1 SV 1 T 1 I 1 V 30 10 V 2*-20 ; x = 0 m, To = -20 °C SV 2 T 1 I 1 V 30 10 V 2*25 ; x = 1 m, To = 25 °C AV 1 T 6 N 1 1 1 V 4*0 ; KAPPR KAUTO KPRED GV 1 T 1 V 27*25 ; initial conditions AS 1 /M 1 /A 1 /G 1 /S 1 E 1 S3 /S 2 E 8 S1 /R 720000 36000 1 ; total time = 200 h, initial time step = 10 h EN EN
# 163 0.72000E+06
$TEMPERATURE
1 -0.1160E+02 0.0000E+00 0.0000E+00
2 -0.1160E+02 0.0000E+00 0.0000E+00
3 -0.8457E+01 0.0000E+00 0.0000E+00
4 -0.8457E+01 0.0000E+00 0.0000E+00
5 -0.5347E+01 0.0000E+00 0.0000E+00
6 -0.5347E+01 0.0000E+00 0.0000E+00
7 -0.2269E+01 0.0000E+00 0.0000E+00
8 -0.2269E+01 0.0000E+00 0.0000E+00
9 0.7758E+00 0.0000E+00 0.0000E+00
10 0.7758E+00 0.0000E+00 0.0000E+00
11 -0.1160E+02 0.0000E+00 0.0000E+00
12 -0.1002E+02 0.0000E+00 0.0000E+00
13 -0.1002E+02 0.0000E+00 0.0000E+00
14 -0.1002E+02 0.0000E+00 0.0000E+00
15 -0.8457E+01 0.0000E+00 0.0000E+00
16 -0.6898E+01 0.0000E+00 0.0000E+00
17 -0.6898E+01 0.0000E+00 0.0000E+00
18 -0.6898E+01 0.0000E+00 0.0000E+00
19 -0.5347E+01 0.0000E+00 0.0000E+00
20 -0.3803E+01 0.0000E+00 0.0000E+00
21 -0.3803E+01 0.0000E+00 0.0000E+00
22 -0.3803E+01 0.0000E+00 0.0000E+00
23 -0.2269E+01 0.0000E+00 0.0000E+00
24 -0.7423E+00 0.0000E+00 0.0000E+00
25 -0.7423E+00 0.0000E+00 0.0000E+00
26 -0.7423E+00 0.0000E+00 0.0000E+00
27 0.7758E+00 0.0000E+00 0.0000E+00