Hexahedra, BEAM56G1

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 node 24
$u=w=0$ at node 22
$v=w=0$ at node 56
$w=0$ at node 54

Solution

The unit end force is replaced with the equivalent surface traction as $$q_x=\frac{-1}{0.01\times 0.02}=-5000\text{ Pa},\quad q_y=0\text{ Pa},\quad q_z=0\text{ Pa}.$$

The initial stress matrix is then computed from the stress field obtained by the standard linear analysis.

beam56g1.bat

rmd3 beam56g1.i1
rpd3 beam56g1.i2
srh3 beam56g1.i3
fefs beam56g1.i4
geo3 beam56g1.iG
heig beam56g1.iE
stab beam56g1.iS

Input

beam56g1.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

beam56g1.i2

;  KREST
IP   1
         ; E    α ν   ρ
MP 1 T 1 V 2e11 0 0.3 0
SV 1 T 9 V -5000 0 0 ; Fx = -1 N

AS 1 /M 1
  /B 0 N 24 /B 0 C 1 3 N 22     ; left end
  /B 0 C 2 3 N 56 /B 0 C 3 N 54 ; right end
  /S 1 E 4 S6

EN
EN

beam56g1.i3

;  KREST
IP   1
EN
EN

beam56g1.i4

;  KREST
IP   1
EN
EN

beam56g1.iG

;  ILC
IP  1
EN
EN

beam56g1.iE

;  KREST NROOT NITERX KTPR KEVP
IP   1     1      0     0    1
EN
EN

beam56g1.iS

IP 1
EN
EN

Output

beam56g1.oS

IP 1                                                                            
EN                                                                              
EN                                                                              


                              LOADING PARAMETERS
  0.348977E+04

                                * END OF STAB *
 TOTAL CPU: 00:00:00

Package for Machine Design

Table of Contents