/TITLE, Your name, ME 421, and date here, trussw.a
  ! Displacement and stress in a Warren truss

/COLOR,PBAK,off    ! black background
/PREP7             ! preprocessor section
p1 = -33e3         ! Load at nodes 2, 3, 4

ET,1,Link1         ! type 1 is a spar
R,1,3.0            ! area for material 1
MP,EX,1,29e6       ! Young's modulus, steel
MP,DENS,1,0.284    ! steel density, lbs/cu in
ACEL,,1            ! gravity

! bottom nodes, don't change the next three lines
N,1,0              ! lower-left node
N,5,40             ! lower-right node
FILL,1,5           ! remaining bottom nodes

! top nodes
N,6,5,5*sqrt(3)    ! upper-left node
N,7,(1*10+5),5*sqrt(3) ! upper-left center node
N,8,(2*10+5),5*sqrt(3) ! upper-right center node
N,9,(3*10+5),5*sqrt(3) ! upper-right node
/PNUM,node,1

TYPE,1             ! element type (ET) 1 is a spar
MAT,1              ! material property MP 1 (steel)
REAL,1             ! original area
! bottom elements
E,1,2              ! bottom left element  1
E,4,5              ! bottom right element 2
E,2,3              ! bottom left-center element 3
E,3,4              ! bottom right-center element 4

! left, top, and right elements
E,1,6              ! left 5
E,9,5              ! right 6
E,6,7              ! top-left 7
E,8,9              ! top-right 8
E,7,8              ! top-center 9

! inner elements
E,2,6              ! left pair 10
E,4,9              !
E,2,7              ! middle pair 12
E,4,8
E,3,7              ! center pair 14
E,3,8
APLOT,all           ! plot all areas
/PLOPTS,INFO,OFF    ! Removes entire legend column
/PLOPTS,FRAME,OFF   ! Removes frame

FINISH             ! preprocessor section

/SOLUTION          ! solution section
antype,static      ! analysis type is static
outpr,basic,1
outpr,nload,1
D,1,UX,0,,,,UY     ! pin at left end
D,5,UY,0           ! roller at right end
F,2,FY,p1          ! apply 3 loads
F,3,FY,p1
F,4,FY,p1

! /OUTPUT,Mass,txt,D:\Ansys\   !  mass
 /OUTPUT,Mass,txt,U:\win\Ansys\   !  mass

solve
/PBC,ALL,1
eplot

FINISH             ! solution section

/POST1             ! post processor
etable,STR,ls,1,3
pldisp,1           ! plot displacements
! /OUTPUT,Displ,txt,D:\Ansys\   !  displacements
 /OUTPUT,Displ,txt,U:\win\Ansys\   !  displacements

PRNSOL,U           ! show displacements
! /OUTPUT,Stress,txt,D:\Ansys\   !  forces and stresses
 /OUTPUT,Stress,txt,U:\win\Ansys\   !  forces and stresses
PRRSOL             ! show forces
PRETAB             ! stresses
! FINISH             ! post processor

!  TYPE    MASS
!     1   127.800

!    NODE      UX             UY
!       1  0.0000         0.0000
!       2  0.32886E-02   -0.38602E-01
!       3  0.10962E-01   -0.53158E-01
!       4  0.18635E-01   -0.38602E-01
!       5  0.21923E-01    0.0000
!       6  0.21922E-01   -0.20252E-01
!       7  0.15346E-01   -0.48096E-01
!       8  0.65773E-02   -0.48096E-01
!       9  0.98946E-06   -0.20252E-01

!    ELEM    STR
!       1   9536.9
!       2   9536.9
!       3   22252.
!       4   22252.
!       5  -19074.
!       6  -19074.
!       7  -19071.
!       8  -19071.
!       9  -25429.
!      10   19069.
!      11   19069.
!      12  -6360.7
!      13  -6360.7
!      14   6354.1
!      15   6354.1