/prep7
!change this line to et,1,82 for PLANE82
et,1,183
!sets the element for using the mixed formulation (u-p).
keyopt,1,6,1
mp,ex,1,1e3
mp,prxy,1,0.4999
k,1,0,0
k,2,2,0
k,3,2,3
k,4,0,2
k,5,0.4,0.4
k,6,1.4,0.6
k,7,1.5,2
k,8,0.3,1.6
!create lines and set their number of division to 1 for each line
l,1,2,1
l,2,3,1
l,3,4,1
l,4,1,1
l,5,6,1
l,6,7,1
l,7,8,1
l,8,5,1
l,1,5,1
l,2,6,1
l,3,7,1
l,4,8,1
a,1,2,6,5
a,3,4,8,7
a,1,5,8,4
a,5,6,7,8
amesh,all
!this area is meshed separately to force a quad mesh
a,2,3,7,6
AMAP,5,6,2,3,7
!apply displacement conditions to the surrounding nodes
d,1,all,0
d,2,ux,0.004
d,2,uy,0
d,9,ux,0.004
d,9,uy,-0.0018
d,10,ux,0
d,10,uy,-0.0012
d,3,ux,0.002
d,3,uy,0
d,20,ux,0.004
d,20,uy,-0.0006*1.5
d,11,ux,0.002
d,11,uy,-0.0006*2.5
d,18,ux,0
d,18,uy,-0.0006
fini