! ANSYS Script for Stress Concentration

/PREP7  

!  Parameters for Geometry

r = 1.0
fac = 5
s = fac * r

! Element and Materials

ET,1,PLANE82
MP,EX,1,1.e7
MP,NUXY,1,0.3

!  Create Basic Geometry

rectng,0,s,0,s
cyl4,0,0,r
asba,1,2
adele,all

!  Create Two Areas for Meshing

ldiv,5
l,1,3
al,1,4,3,10
al,5,9,2,4

!  Parameters for Meshing

nd = 5
nt = 20
nr = 10

!  Meshing controls

lesize,1,,,nd
lesize,5,,,nd
lesize,10,,,nt,1/nr
lesize,9,,,nt,1/nr
lesize,4,,,nt,nr

!  Meshing commands

mshkey,1
amesh,all

!  Put in loads and boundary conditions

dl,9,,symm
dl,10,,symm
sfl,2,pres,-1000

finish

!  Solve the load step

/solu
solve
finish