Package for Machine Design

Finite Element Analysis in Structural Mechanics

User Tools

Site Tools


en:example:dyn:1:1

Table of Contents

Hexahedra, BEAM56D1

Problem

Mesh

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

Boundary conditions

$u=v=w=0$ at nodes 1, 2, 3, 4, 21, 22, 23, 24

Solution

beam56d1.bat
rmd3 beam56d1.i1
rpd3 beam56d1.i2
srh3 beam56d1.i3
fefs beam56d1.i4
hmot beam56d1.iM
heig beam56d1.iE
hfrq beam56d1.iF

Input

beam56d1.i1
;  NELEM NNOD ITED
IP   4    56   56
;  CRIT SCALE
RP 1.01   1
XY N 1:20 X 4*0 4*0.25 4*0.5 4*0.75 4*1
          Y 5*(0 0.02 0.02 0)
          Z 5*(0 0 0.01 0.01)
EL E 1:4 N =A 1:8 =A =B 21:32 =B
           =4A =8B =8A =16B =12A =24B
EN
EN
beam56d1.i2
;  KREST
IP   1
         ; E    α ν   ρ
MP 1 T 1 V 2e11 0 0.3 7800
AS 1 /M 1 /B 0 N 1:4 21:24
EN
EN
beam56d1.i3
;  KREST
IP   1
EN
EN
beam56d1.i4
;  KREST
IP   1
EN
EN
beam56d1.iM
;  KDIAG KPRIN
IP   1     0
EN
EN
beam56d1.iE
;  KREST NROOT NITERX KTPR KEVP
IP   1     10     0     0    0
EN
EN
beam56d1.iF
;  KPRIN KOUT IVL IVU
IP   0     0   0   0
EN
EN

Output

beam56d1.oF
;  KPRIN KOUT IVL IVU                                                           
IP   0     0   0   0                                                            
EN                                                                              
EN                                                                              


                              FREQUENCIES [1/s]
  0.854545E+01  0.169337E+02  0.580379E+02  0.113843E+03  0.190357E+03
  0.365610E+03  0.491760E+03  0.636444E+03  0.940207E+03  0.127715E+04

                                * END OF HFRQ *
 TOTAL CPU: 00:00:00
en/example/dyn/1/1.txt · Last modified: 2022-03-15 08:28 by Petr Pařík