Package for Machine Design

Finite Element Analysis in Structural Mechanics

User Tools

Site Tools


en:example:heat:2:2

Table of Contents

Triangles, BEAM4T2

Problem

Mesh

element type nodes
1 4 1 3 2 12 13 11
2 4 3 4 2 15 14 13 
3 4 3 5 4 16 17 15
4 4 5 6 4 19 18 17
5 4 5 7 6 20 21 19
6 4 7 8 6 23 22 21
7 4 7 9 8 24 25 23
8 4 9 10 8 27 26 25

Computation

The computation is executed with the following commands:

beam4t2.bat
xrm2 beam4t2.i1
xrpd beam4t2.iB
xt2t beam4t2

Input

beam4t2.i1
;  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
beam4t2.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 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

Output

beam4t2.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.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
en/example/heat/2/2.txt · Last modified: 2024-11-07 09:26 by Petr Pařík