/TITLE, Your name, ME 421, and date here, therm6.a
     ! thermal load stepped bar, fixed right end
     ! Beer prob 1.2, page 13

/PREP7             ! preprocessor section
p1 = -12e3         ! Load at node 1
p2 = (2*25e3)      ! Load at node 1
temper=35          ! temperature change
d1 = 1.5           ! diam 1
d2 = 2.5           ! diam 2

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,11e6       ! Young-s modulus, aluminum
MP,EX,2,6.5e6      ! Young-s modulus, magnesium
MP,ALPX,1,13e-6    ! thermal expansion, aluminum
MP,ALPX,2,14e-6    ! thermal expansion, magnesium
N,1,0              ! Node 1 at 0
N,2,3*12           ! Node 2
N,3,(3+4)*12       ! Node 3

TYPE,1             ! element type1 (pipe)
REAL,1             ! area R 1
MAT,1              ! material property MP 1
E,1,2              ! Element 1, nodes 1-2
REAL,2             ! Area 2
MAT,2              ! material property MP 2
E,2,3              ! Element 2, nodes 2-3
FINISH             ! preprocessor section

/SOLUTION          ! solution section
antype,static
outpr,basic,1
outpr,nload,1
D,3,ALL            ! Fix node 3
F,1,FX,p1          ! force at node 1
F,2,FX,p2          ! 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.49573E-02
!       2  0.33646E-01
!       3  0.

!    NODE     FX
!       3   -38000.

!    ELEM    STR
!       1   6790.6
!       2  -7741.3