Package for Machine Design

Finite Element Analysis in Structural Mechanics

User Tools

Site Tools


en:example:stat:3:3

Beam elements, BEAM53S1

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

The end force $F_x,$ $F_y,$ $F_z$ is directly assigned to node 5 as a concentrated (nodal) force. The distributed loading is replaced with the equivalent body force $$\begin{align} Q_x&=0, \\ Q_y&=\frac{-100}{0.01\times0.02}=-500\times10^3\text{ N/m}^3, \\ Q_z&=1000\times10^3\text{ N/m}^3, \end{align}$$ which is assigned to elements 1 to 4.

beam53s1.bat
rmd3 beam53s1.i1
rpd3 beam53s1.i2
srh3 beam53s1.i3
fefs beam53s1.i4
str3 beam53s1.i5

Input

beam53s1.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
beam53s1.i2
;  KREST
IP   1
         ; E    α ν   ρ
MP 1 T 1 V 2e11 0 0.3 0
NV 1 T 6 V 20000 -20 0     ; Fx = 20000 N, Fy = -20 N, Fz =  0 N
NV 2 T 6 V 0 0 10          ; Fx =     0 N, Fy =   0 N, Fz = 10 N
VV 1 T 6 V 0 -500e3 0      ; ly = -100 N/m
VV 2 T 6 V 0 -500e3 1000e3 ; ly = -100 N/m, lz = 200 N/m

AS 1 /M 1
  /B 0 N 1
  ; in-plane loading ly, Fx, Fy
  /N 1 N 5
  /V 1 E 1:4

AS 2
  ; three-dimensional loading ly, lz, Fx, Fy, Fz
  /N 1 N 5
  /N 2 N 5
  /V 2 E 1:4

EN
EN
beam53s1.i3
;  KREST
IP   1
EN
EN
beam53s1.i4
;  KREST
IP   1
EN
EN
beam53s1.i5
;  KLC 0 KOUT ILC  ...  KPROB KGRAF
IP  1  0   1   0   6*0    0     2
EN
EN

Output

beam53s1.o5
;  KLC 0 KOUT ILC  ...  KPROB KGRAF                                             
IP  1  0   1   0   6*0    0     2                                               
EN                                                                              
EN                                                                              
                                                      1TH LOAD CASE

                                                    DISPLACEMENTS
           NODE             U                V                W             ALPHA            BETA             GAMMA
                           [m]              [m]              [m]           [rad]            [rad]            [rad]

              1       0.0000000E+00    0.0000000E+00    0.0000000E+00    0.0000000E+00    0.0000000E+00    0.0000000E+00
              2       0.1250000E-03   -0.1418457E-02    0.0000000E+00    0.0000000E+00    0.0000000E+00   -0.1050781E-01
              3       0.2500000E-03   -0.4882813E-02    0.0000000E+00    0.0000000E+00    0.0000000E+00   -0.1656250E-01
              4       0.3750000E-03   -0.9426270E-02    0.0000000E+00    0.0000000E+00    0.0000000E+00   -0.1933594E-01
              5       0.5000000E-03   -0.1437500E-01    0.0000000E+00    0.0000000E+00    0.0000000E+00   -0.2000000E-01

                                                       REACTIONS
           NODE             RX               RY               RZ            RALPHA           RBETA            RGAMMA
                           [N]              [N]              [N]            [Nm]             [Nm]             [Nm]

              1      -0.2000000E+05    0.1200000E+03    0.0000000E+00    0.0000000E+00    0.0000000E+00    0.7000000E+02



                                                 EQUILIBRIUM OF FORCES AND REACTIONS

                                                 SUM OF FORCES [N]             SUM OF REACTIONS [N]


                         DIRECTION X:                       20000.00                     -20000.00
                         DIRECTION Y:                        -120.00                        120.00
                         DIRECTION Z:                           0.00                          0.00




                                                   STRESSES AT GAUSSPOINTS



    IE   KSI_1      KSI_2    SIG-KSI SIG-ETA_1 SIG-ZETA_1 SIG-ETA_2 SIG-ZETA_2   TAU-K    HMH_1    HMH_2
          [m]        [m]       [MPa]     [MPa]      [MPa]     [MPa]      [MPa]   [MPa]    [MPa]    [MPa]

     1 0.0000E+00 0.0000E+00   100.0     105.0        0.0     105.0        0.0     0.0    205.0    205.0
     2 0.0000E+00 0.0000E+00   100.0      64.7        0.0      64.7        0.0     0.0    164.7    164.7
     3 0.0000E+00 0.0000E+00   100.0      33.7        0.0      33.7        0.0     0.0    133.7    133.7
     4 0.0000E+00 0.0000E+00   100.0      12.2        0.0      12.2        0.0     0.0    112.2    112.2
                                                      2TH LOAD CASE

                                                    DISPLACEMENTS
           NODE             U                V                W             ALPHA            BETA             GAMMA
                           [m]              [m]              [m]           [rad]            [rad]            [rad]

              1       0.0000000E+00    0.0000000E+00    0.0000000E+00    0.0000000E+00    0.0000000E+00    0.0000000E+00
              2       0.1250000E-03   -0.1418457E-02    0.8769531E-02    0.0000000E+00   -0.6437500E-01   -0.1050781E-01
              3       0.2500000E-03   -0.4882813E-02    0.2968750E-01    0.0000000E+00   -0.9875000E-01   -0.1656250E-01
              4       0.3750000E-03   -0.9426270E-02    0.5642578E-01    0.0000000E+00   -0.1125000E+00   -0.1933594E-01
              5       0.5000000E-03   -0.1437500E-01    0.8500000E-01    0.0000000E+00   -0.1150000E+00   -0.2000000E-01

                                                       REACTIONS
           NODE             RX               RY               RZ            RALPHA           RBETA            RGAMMA
                           [N]              [N]              [N]            [Nm]             [Nm]             [Nm]

              1      -0.2000000E+05    0.1200000E+03   -0.2100000E+03    0.0000000E+00    0.1100000E+03    0.7000000E+02



                                                 EQUILIBRIUM OF FORCES AND REACTIONS

                                                 SUM OF FORCES [N]             SUM OF REACTIONS [N]


                         DIRECTION X:                       20000.00                     -20000.00
                         DIRECTION Y:                        -120.00                        120.00
                         DIRECTION Z:                         210.00                       -210.00




                                                   STRESSES AT GAUSSPOINTS



    IE   KSI_1      KSI_2    SIG-KSI SIG-ETA_1 SIG-ZETA_1 SIG-ETA_2 SIG-ZETA_2   TAU-K    HMH_1    HMH_2
          [m]        [m]       [MPa]     [MPa]      [MPa]     [MPa]      [MPa]   [MPa]    [MPa]    [MPa]

     1 0.0000E+00 0.0000E+00   100.0     105.0      330.0     105.0      330.0     0.0    535.0    535.0
     2 0.0000E+00 0.0000E+00   100.0      64.7      191.2      64.7      191.2     0.0    355.9    355.9
     3 0.0000E+00 0.0000E+00   100.0      33.7       90.0      33.7       90.0     0.0    223.7    223.7
     4 0.0000E+00 0.0000E+00   100.0      12.2       26.2      12.2       26.2     0.0    138.4    138.4

                                * END OF STR3 *
 TOTAL CPU: 00:00:00
en/example/stat/3/3.txt · Last modified: 2022-03-15 08:24 by Petr Pařík