en:example:heat:2:1
Table of Contents
Quadrilaterals, BEAM6T2
Problem
Mesh
| element | type | nodes |
|---|---|---|
| 1 | 6 | 1 3 4 2 12 14 13 11 |
| 2 | 6 | 3 5 6 4 15 17 16 14 |
| 3 | 6 | 5 7 8 6 18 20 19 17 |
| 4 | 6 | 7 9 10 8 21 23 22 20 |
Solution
- beam6t2.bat
xrm2 beam6t2.i1 xrpd beam6t2.iB xt2t beam6t2
Input
- beam6t2.i1
; NELEM NNOD ITED IP 4 23 6 ; 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:4 N =A 1 3 4 2 =A =B 12 14 13 11 =B =2A =3B =4A =6B =6A =9B EN EN
- beam6t2.iB
; 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 23*25 ; initial conditions AS 1 /M 1 /A 1 /G 1 /S 1 E 1 S4 /S 2 E 4 S2 /R 720000 36000 1 ; total time = 200 h, initial time step = 10 h EN EN
Output
- beam6t2.STR
# 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.8457E+01 0.0000E+00 0.0000E+00 15 -0.6898E+01 0.0000E+00 0.0000E+00 16 -0.6898E+01 0.0000E+00 0.0000E+00 17 -0.5347E+01 0.0000E+00 0.0000E+00 18 -0.3803E+01 0.0000E+00 0.0000E+00 19 -0.3803E+01 0.0000E+00 0.0000E+00 20 -0.2269E+01 0.0000E+00 0.0000E+00 21 -0.7422E+00 0.0000E+00 0.0000E+00 22 -0.7422E+00 0.0000E+00 0.0000E+00 23 0.7758E+00 0.0000E+00 0.0000E+00
en/example/heat/2/1.txt · Last modified: 2022-03-15 08:27 by Petr Pařík