/TITLE, Your name, ME 421, and date here, torq2.a
  ! Ansys: find rotation and stress
  ! Beer 3.28, 3 gears, 2 shafts

/PREP7             ! preprocessor section
torq1 = 4e3         ! moment at node 1
torq2 = -10e3       ! moment at node 2
dia1 = 1.25
dia2 = 1.5

ET,1,Pipe16        ! type 16 is a pipe
R,1,dia1,dia1/2    ! diam, wall thickness
R,2,dia2,dia2/2
MP,EX,1,29E6       ! elastic modulus, steel
MP,GXY,1,11.2E6    ! shear modulus

N,1,(3+4)*12       ! use X axis
N,2,3*12
N,3,0              ! node 3 at origin

TYPE,1             ! element type (ET) 1
MAT,1              ! material property MP 1
REAL,1             ! area R 1
E,1,2
REAL,2             ! shaft 2-3
E,2,3
FINISH             ! preprocessor section

/SOLUTION          ! solution section
antype,static      ! analysis type is static
outpr,basic,1
outpr,nload,1
D,1,ALL            ! node 1 fully constrained
D,3,ALL            ! node 3 fully constrained
F,2,MX,torq2       ! moment at node 2
! /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
!   1     0.
!   2    -0.474943E-01
!   3     0.

!  NODE     MX
!   1     2656.18
!   3     7343.82

!  S1mx and S3mn are principal stresses (and max shear stress)
! EL= 1  NODES= 1  2  MAT= 1   S1MX=  6926.2
! EL= 2  NODES= 2  3  MAT= 1   S1MX= 11082.