/TITLE, Your name and date here, thermE01.a
  ! thermal load on stepped bar, both ends fixed
  ! Hibler prob 4-37, 38, page 146
/PREP7            ! preprocessor section
p2 = -2*75e3      ! Load at node 2
p3 =  2*100e3     ! Load at node 3
temper = 65       ! temperature
d1 = 50e-3        ! diam 1, aluminum
d2 = 20e-3        ! diam 2, steel

ET,1,Pipe16       ! Type 16 is a pipe
R,1,d1,d1/2       ! diam, wall thickness
R,2,d2,d2/2       ! diam, wall thickness
MP,EX,1, 70E9     ! Young-s modulus, aluminum
MP,EX,2,200E9     ! Young-s modulus, steel
MP,ALPX,1,23e-6   ! thermal expansion, aluminum
MP,ALPX,2,11.7e-6 ! thermal expansion, steel

N,1,0             ! Node 1 at 0
N,2,0.25          ! Node 2
N,3,(0.25+0.5)    ! Node 3
N,4,(2*0.25+0.5)  ! Node 3

TYPE,1
REAL,1            ! Area 1
MAT,1             ! aluminum
E,1,2             ! Element 1, nodes 1-2
E,3,4             ! Element 1, nodes 1-2
REAL,2            ! Area 2
MAT,2             ! steel
E,2,3             ! Element 2, nodes 2-3
FINISH            ! preprocessor section

/SOLUTION         ! solution section
antype,static
outpr,basic,1
outpr,nload,1
D,1,ALL,0         ! Fix node 1
D,4,ALL,0         ! Fix node 3
F,2,FX,p2         ! force at node 2
F,3,FX,p3         ! force at node 2
BFE,ALL,TEMP,1,temper  ! Temperature all elements
solve
/PBC,ALL,1
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

!    NODE      UX
!       1   0.0000
!       2   0.23873E-04
!       3   0.67073E-04
!       4   0.0000

!    NODE      FX
!       1   0.19235E+06
!       4  -0.24235E+06

!    ELEM    STR
!       1 -0.97966E+08
!       2 -0.12343E+09
!       3 -0.13482E+09