Package for Machine Design

Finite Element Analysis in Structural Mechanics

User Tools

Site Tools


en:example:heat:1:5

Transition elements, BEAM71T1

Problem

Mesh

element type nodes
1 56 1 2 3 4 5 6 7 8 25 26 27 28 29 30 31 32 33 34 35 36
2 61 9 10 12 11 37 38 39 40
3 56 13 14 15 16 17 18 19 20 41 42 43 44 45 46 47 48 49 50 51 52
4 61 21 22 24 23 53 54 55 56
5 71 element 1 face S6 ↔ element 2 edge L1 
6 71 element 2 edge L3 ↔ element 3 face S1
7 71 element 3 face S6 ↔ element 4 edge L1

Solution

beam71t1.bat
xrm3 beam71t1.i1
xrpd beam71t1.iB
xt3s beam71t1

Input

beam71t1.i1
;  NELEM NNOD ITED
IP   7    56   56
;  CRIT SCALE THDEF
RP 1.01   1    0.01
XY N 1:24 X 4*0 6*0.25 6*0.5 6*0.75 2*1
          Y =A 2*(0 0.02 0.02 0) =A
            =B 2*(0 0.02) =B =A =B
          Z =A 2*(0 0 0.01 0.01) =A
            =B 4*0.005 =B =A =B
EL E 1 3 N =A 1:8 =A =B 25:36 =B =12A =16B
EL T 61 E 2 4 N =A 9 10 12 11 =A =B 37:40 =B =12A =16B
CN T -3 E 5  1 6  2 1
        E 6  2 3  3 1
        E 7  3 6  4 1
EN
EN
beam71t1.iB
;  KREST 0 KOUT INT3 NSAX NSTEPX KSU
IP   1   0   2    0   10     0    1
;  TIMS ERAL EDIF
RP   0  1e-5  0.1
MP 1 T 1 V 20 0
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
LV 1 T 1 I 1 V 30 10 V 2*25  ; x = 1 m, To =  25 °C
AV 1 T 6 N 1 0 0 V 4*0       ; KAPPR
GV 1 T 1 V 72*25             ; initial estimate

AS 1 /M 1 /A 1 /G 1
  /S 1 E 1 S1
  /L 1 E 4 L3

EN
EN

Output

beam71t1.STR
#    1    0.00000E+00
$TEMPERATURE
      1 -0.1183E+02  0.0000E+00  0.0000E+00
      2 -0.1183E+02  0.0000E+00  0.0000E+00
      3 -0.1183E+02  0.0000E+00  0.0000E+00
      4 -0.1183E+02  0.0000E+00  0.0000E+00
      5 -0.8750E+01  0.0000E+00  0.0000E+00
      6 -0.8750E+01  0.0000E+00  0.0000E+00
      7 -0.8750E+01  0.0000E+00  0.0000E+00
      8 -0.8750E+01  0.0000E+00  0.0000E+00
      9 -0.8750E+01 -0.9660E-13  0.0000E+00
     10 -0.8750E+01  0.8158E-13  0.0000E+00
     11 -0.5684E+01 -0.1416E-13  0.0000E+00
     12 -0.5684E+01  0.2719E-14  0.0000E+00
     13 -0.5684E+01  0.0000E+00  0.0000E+00
     14 -0.5684E+01  0.0000E+00  0.0000E+00
     15 -0.5684E+01  0.0000E+00  0.0000E+00
     16 -0.5684E+01  0.0000E+00  0.0000E+00
     17 -0.2593E+01  0.0000E+00  0.0000E+00
     18 -0.2593E+01  0.0000E+00  0.0000E+00
     19 -0.2593E+01  0.0000E+00  0.0000E+00
     20 -0.2593E+01  0.0000E+00  0.0000E+00
     21 -0.2593E+01 -0.4532E-14  0.0000E+00
     22 -0.2593E+01  0.1889E-14  0.0000E+00
     23  0.4728E+00 -0.1573E-14  0.0000E+00
     24  0.4728E+00  0.5332E-15  0.0000E+00
     25 -0.1182E+02  0.0000E+00  0.0000E+00
     26 -0.1183E+02  0.0000E+00  0.0000E+00
     27 -0.1182E+02  0.0000E+00  0.0000E+00
     28 -0.1183E+02  0.0000E+00  0.0000E+00
     29 -0.1029E+02  0.0000E+00  0.0000E+00
     30 -0.1029E+02  0.0000E+00  0.0000E+00
     31 -0.1029E+02  0.0000E+00  0.0000E+00
     32 -0.1029E+02  0.0000E+00  0.0000E+00
     33 -0.8760E+01  0.0000E+00  0.0000E+00
     34 -0.8753E+01  0.0000E+00  0.0000E+00
     35 -0.8760E+01  0.0000E+00  0.0000E+00
     36 -0.8753E+01  0.0000E+00  0.0000E+00
     37 -0.8750E+01 -0.7464E-14  0.0000E+00
     38 -0.7217E+01 -0.1114E-13  0.0000E+00
     39 -0.5684E+01 -0.5446E-14  0.0000E+00
     40 -0.7217E+01  0.1416E-13  0.0000E+00
     41 -0.5669E+01  0.0000E+00  0.0000E+00
     42 -0.5680E+01  0.0000E+00  0.0000E+00
     43 -0.5669E+01  0.0000E+00  0.0000E+00
     44 -0.5680E+01  0.0000E+00  0.0000E+00
     45 -0.4139E+01  0.0000E+00  0.0000E+00
     46 -0.4139E+01  0.0000E+00  0.0000E+00
     47 -0.4139E+01  0.0000E+00  0.0000E+00
     48 -0.4139E+01  0.0000E+00  0.0000E+00
     49 -0.2608E+01  0.0000E+00  0.0000E+00
     50 -0.2597E+01  0.0000E+00  0.0000E+00
     51 -0.2608E+01  0.0000E+00  0.0000E+00
     52 -0.2597E+01  0.0000E+00  0.0000E+00
     53 -0.2593E+01 -0.1471E-14  0.0000E+00
     54 -0.1060E+01 -0.3158E-15  0.0000E+00
     55  0.4727E+00 -0.4319E-15  0.0000E+00
     56 -0.1060E+01  0.7517E-15  0.0000E+00
en/example/heat/1/5.txt · Last modified: 2022-03-15 08:26 by Petr Pařík