/batch
/FILNAM,beam1
/title 2D Beam Sample Problem 1 - 10 elements

! Ten element 2D beam Sample Problem  10-1-98   KLL

/prep7
n,  1,  0.0,  0.0    * Node 1 is located at (0.0, 0.0)
n,  2,  1.0,  0.0
n,  3,  2.0,  0.0
n,  4,  3.0,  0.0
n,  5,  4.0,  0.0
n,  6,  5.0,  0.0
n,  7,  6.0,  0.0
n,  8,  7.0,  0.0
n,  9,  8.0,  0.0
n, 10,  9.0,  0.0
n, 11, 10.0,  0.0

et, 1, beam3         * Element type; no.1 is beam3

mp, ex, 1, 3.e7      * Material Property, Elastic modulus no. 1
mp, prxy, 1, 0.3

! Beam real constant set 1 for a 0.5 x 0.375 rectangular xsctn.
! Area, Izz (flexural Inertia), height 'h' as in sigma = Mc/I, c = h/2
! A = 0.1875, Izz = 0.0022, h = 0.375

r, 1, 0.1875, 0.0022, 0.375   

en,  1,  1,  2       * Element Number 1 connects nodes 1 & 2
en,  2,  2,  3
en,  3,  3,  4
en,  4,  4,  5
en,  5,  5,  6
en,  6,  6,  7
en,  7,  7,  8
en,  8,  8,  9
en,  9,  9, 10
en, 10, 10, 11

d, 1, ux, 0.         * Displacement at node 1 in x-dir is zero
d, 1, uy, 0.         * Displacement at node 1 in y-dir is zero
d, 1, rotz, 0.       * Rotation about z axis at node 1 is zero

f, 11, fy, -50.      * Force at node 11 in y-direction

/pnum, elem, 1       * Plot element numbers
eplot                * Plot the elements

finish

/solu                * Select static load solution
antype, static
solve
save
finish

/post1
etable, stress-x, s, max    * Create a table of element stress values
                            * P/A + Mc/Izz