/TITLE, simply supported, triangular load, triangmat.a
     ! load from support 0 to 60 to 0 lb/in, 
     ! two elements, W8 x 31
     ! convert to concentrated load for Matlab
     ! rotation and deflection at nodes are correct
/PREP7          ! preprocessor section
f2 = -4200      ! force
m3 = 10e3       ! moment

ET,1,Beam3      ! Type 3 is beam
MP,EX,1,30E6    ! Young-s modulus, steel
MP,GXY,1,11.2E6 ! shear modulus
R,1,9.1,110,8   ! Area, Iz, height
N,1,0           ! Node 1 at 0
N,2,100          ! Node 2 at 50
N,3,200         ! Node 2 at 100
E,1,2           ! Element 1
E,2,3           ! Element 2
FINISH          ! preprocessor section

/SOLUTION       ! solution section
antype,static
outpr,basic,1
outpr,nload,1
D,1,UX,0,,,,UY  ! Pin node 1
D,3,UY,0        ! Fix Y at node 3
F,2,FY,f2       ! force at node 2
F,3,MZ,m3       ! Moment at node 1
F,1,MZ,-m3      ! Moment at node 3
/output,Displ,txt,D:\Ansys\   !  displacements
! /output,Displ,txt,U:\win\Ansys\   !  displacements
solve
/PBC,ALL,1
FINISH          ! solution section

/POST1          ! post processor
etable,STR,ls,1,3
/output,Force,txt,D:\Ansys\   !  forces and stresses
! /output,Force,txt,U:\win\Ansys\   !  forces and stresses
pretab
prrsol
FINISH           ! post processor

!    NODE      UX              UY              ROTZ
!       1     0.00000         0.00000       -0.378788E-02
!       2     0.00000        -0.242424       0.00000    
!       3     0.00000         0.00000        0.378788E-02

! Stress and forces are incorrect
!EL=       1  NODES=       1       2  MAT=  1    BEAM3
! LOCATION   SMAX        SMIN     
!  1 (I)    363.64     -363.64    
!  2 (J)    8000.0     -8000.0    

! EL=       2  NODES=       2       3  MAT=  1   BEAM3
! LOCATION   SMAX        SMIN     
!  1 (I)    8000.0     -8000.0    
!  2 (J)    363.64     -363.64    

!    NODE      FX              FY              MZ  
!       1     0.00000         2100.00    
!       3                     2100.00