finish
/clear

/title,Contact Elements
/filnam,conta3,1
/prep7

! Top Beam
X1=0
Y1=15
L1=100
H1=10

! Bottom Beam
X2=50
Y2=0
L2=100
H2=10

! Create Geometry
blc4,X1,Y1,L1,H1
blc4,X2,Y2,L2,H2

! define element type

ET,1,plane42	    ! element type 1
keyopt,1,3,3        ! plane stress w/thick 
type,1              ! activate element type 1
R, 1, 10            ! thickness 0.01   

! define material properties

MP,EX,  1, 200e3    ! Young's modulus
MP,NUXY,1, 0.3      ! Poisson's ratio

! meshing

esize,2             ! set meshing size
amesh,all           ! mesh area 1

ET,2,169
ET,4,172            

R,3,,,0.1,0.1,0,   
RMORE,,,1.0E20,0.0,1.0, 
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,4,3,0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,7,2
KEYOPT,4,8,0
KEYOPT,4,9,1
KEYOPT,4,10,2   
KEYOPT,4,11,0   
KEYOPT,4,12,0   
KEYOPT,4,2,0
REAL,3

LSEL,S,,,7  
 
TYPE,2  
NSLL,S,1
ESLN,S,0
ESURF 

! Generate the contact surface  
LSEL,S,,,1  

TYPE,4  
NSLL,S,1
ESLN,S,0
ESURF   
ALLSEL 

    

finish
/solut
antype,0

time,1			! Sets time at end of run to 1 sec
autots,on		! Auto time-stepping on
nsubst,100,1000,20	! Number of sub-steps
outres,all,last		! Write all output
neqit,100		! Max number of iterations
NLGEOM,1

nsel,s,loc,x,X1		! Constrain top beam
nsel,r,loc,y,Y1,(Y1+H1)
d,all,all
nsel,all

nsel,s,loc,x,(X2+L2)	! Constrain bottom beam
nsel,r,loc,y,Y2,(Y2+H2)
d,all,all
nsel,all

nsel,s,loc,x,(L1/2+4*H1)	! Apply load
nsel,r,loc,y,(Y1+H1)
f,all,fy,-10000
nsel,all

solve
finish

/post1
/dscale,1,1

PLNSOL,S,EQV,0,1