Package for Machine Design

Finite Element Analysis in Structural Mechanics

User Tools

Site Tools


en:example:dyn:1:3

Beam elements, BEAM53D1

Problem

Mesh

element type nodes
1 53 1 2
2 53 2 3
3 53 3 4
4 53 4 5

Boundary conditions

$u=v=w=\varphi_x=\varphi_y=\varphi_z=0$ at node 1

Solution

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

Input

beam53d1.i1
;  NELEM NNOD ITED
IP   4     5   53
;  CRIT SCALE ...
RP 1.01   1   8*0
 ; default values of cross-sectional characteristics
 ; A    Ik              Wk                Ieta
   2e-4 6.6666666666e-9 6.666666666666e-7 1.66666666666e-9
 ; Weta             Izeta               Wzeta              Px Py Pz
   3.33333333333e-7 6.66666666666666e-9 6.666666666666e-7  0  1  0
XY N 1:5 X 5*0:1 Y 5*0 Z 5*0
EL E 1:4 N 1 2  2 3  3 4  4 5
EN
EN
beam53d1.i2
;  KREST
IP   1
         ; E    α ν   ρ
MP 1 T 1 V 2e11 0 0.3 7800
AS 1 /M 1 /B 0 N 1
EN
EN
beam53d1.i3
;  KREST
IP   1
EN
EN
beam53d1.i4
;  KREST
IP   1
EN
EN
beam53d1.iM
;  KDIAG KPRIN
IP   1     0
EN
EN
beam53d1.iE
;  KREST NROOT NITERX KTPR KEVP
IP   1     10     0     0    0
EN
EN
beam53d1.iF
;  KPRIN KOUT IVL IVU
IP   0     0   0   0
EN
EN

Output

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


                              FREQUENCIES [1/s]
  0.818017E+01  0.163603E+02  0.513223E+02  0.102645E+03  0.144648E+03
  0.285359E+03  0.289296E+03  0.530754E+03  0.570718E+03  0.706729E+03

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