/TITLE, cantilever, tapered load, canttmat.a
     ! load from support 0 to 60 lb/in, W8 x 31
     ! one element
     ! convert to concentrated load for Matlab

/PREP7          ! preprocessor section
w2 = 60         ! uniform load
f2 = -2100      ! force
m2 = 30e3       ! 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 100
E,1,2           ! Element 1, nodes 1-2
FINISH           ! preprocessor section

/SOLUTION        ! solution section
antype,static
outpr,basic,1
outpr,nload,1
D,1,ALL               ! Fix node 1
F,2,FY,f2       ! force at node 2
F,2,MZ,m2       ! Moment at node 1

! SFBEAM,1,1,PRES,0,w2  ! Tapered load 0-60 lb/in
/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.00000    
!     2     0.00000    -0.166667   -0.227273E-02

! Incorrect stress
! EL=       1  NODES=       1       2  MAT=  1  BEAM3
! LOCATION   SMAX        SMIN     
!  1 (I)    6545.5     -6545.5    
!  2 (J)    1090.9     -1090.9    

! should be the following
! EL=       1  NODES=       1      2  MAT=  1  BEAM3
! PRES LOAD KEY = 1  FACE NODES =       1           2
!                PRESSURES(F/L) =   0.0000      60.000    
! LOCATION   SMAX        SMIN     
!  1 (I)    7272.7     -7272.7    
!  2 (J)    0.0000      0.0000    

! incorrect force, moment
!    NODE      FX              FY              MZ  
!       1     0.00000         2100.00         180000.    

! should be the following
!    NODE      FX          FY          MZ    
!       1   0.0000      3000.0     0.20000E+06