/TITLE, Your name, ME 421, and date here, torq3.a
  ! rotation and stress, shaft with hole

/PREP7              ! preprocessor section
torq2 = 120         ! moment at node 1
torq3 = -150        ! moment at node 2
d1 = 20e-3
d2 = 16e-3

ET,1,Pipe16         ! straight pipe
R,1,d1,d1/2         ! 20 mm diam
R,2,d1,(d1-d2)/2    ! with 16mm hole
MP,EX,1,200E9       ! elastic modulus
MP,GXY,1,77e9       ! shear modulus

N,1,0               ! node 1 at origin
N,2,0.25/2          ! use X axis
N,3,0.25            ! shaft is 250mm long

TYPE,1
MAT,1
REAL,1              ! shaft 1-2
E,1,2
REAL,2              ! shaft 2-3
E,2,3
FINISH              ! preprocessor section

/SOLUTION           ! solution section
antype,static
outpr,basic,1
outpr,nload,1

D,1,ALL             ! fix node 1
F,2,MX,torq2        ! moment at node 2
F,3,MX,torq3        ! moment at node 3
/OUTPUT,Stress,txt,D:\Ansys\   !  moments and stresses
! /OUTPUT,Stress,txt,U:\win\Ansys\   !  moments and stresses
solve
/PBC,ALL,1
FINISH             ! solution section

/POST1             ! post processor
etable,STR,LS,1,3
/OUTPUT,Rotate,txt,D:\Ansys\   !  rotation
! /OUTPUT,Rotate,txt,U:\win\Ansys\   !  rotation
! FINISH             ! post processor

!    NODE      ROTX            NODE       MX
!      1    0.                  1      30.0000
!      2   -0.310042E-02
!      3   -0.293574E-01

! EL= 1  NODES= 1  2  MAT= 1  S1MX= 0.19099E+08
! EL= 2  NODES= 2  3  MAT= 1  S1MX= 0.16174E+09 (max shear stress)