!		2D coupled field model, PZT, Adhesive, Steel - 2nodes support

!				Untersuchung der Bruchdehnung


fini
/cle

!**************** material property values *****************

Es=200e9			!steel
nus=0.28
rs=20e-3
hs=200e-6

lr=2.3e-3
load=0.1


!******************** element types ****************************
/prep7

et,1,plane183,,,1           !Steel
type,1              
mp,ex,1,200e6
mp,prxy,1,0.3
r,1

et,2,plane183,,,1
type,2
mp,ex,2,200e6
mp,prxy,2,0.3

et,3,targe169,,1	! Contact 
et,4,conta175,,1
keyo,4,7,2
r,2
r,3

!********************** geometry ***************************

BLC4,0,0,rs,hs
BLC4,0,0,lr,hs
AOVLAP,ALL 
    
BLC4,(rs-2e-3),-2e-3,4e-3,2e-3    
     
asel,all                 


!********************* mesh **********************************

mshkey,1
mshape,0,2d
ESIZE,(rs/200)

			
mat,1			!Steel Plate
type,1
real,1	
lesize,4,,,3
lesize,2,,,3
lesize,6,,,3
amesh,3
amesh,4

mat,2
type,2				!Support
real,1
amesh,1

allsel




!*****************  defining source and target nodes****************

asel,s,,,3
nsla,s,1
nsel,r,loc,x,(rs-2e-3),rs
nsel,r,loc,y,0
cm,source,node

asel,s,,,1
nsla,s,1
nsel,r,loc,x,(rs-2e-3),rs
nsel,r,loc,y,0
cm,target,node

!******************** contact elements***************************
allsel

type,3
mat,3
real,2
cmsel,s,target
esurf,,target

allsel
type,4
mat,3
cmsel,s,source
esurf,,source


type,3
mat,3
real,3
cmsel,s,source
esurf,,source

type,4
cmsel,s,target
esurf,,target

allsel

!********************** loads ***********************************


!nsel,s,loc,x,lr
!nsel,r,loc,y,hs			                
!*get,num,node,0,count         
!F,all,fy,-load   

nsel,s,loc,x,0
nsel,r,loc,y,0
d,all,uy,-0.42e-3

allsel    

!******************* boundary conditions ***********************

asel,s,,,1			!support
nsla,s,1
d,all,ux,0
d,all,uy,0

!nsel,all			!center constrain
!nsel,s,loc,x,0
!d,all,ux,0

lsel,s,loc,x,0        		!Symmetry constraints
dl,all,,symm
nsel,all   

fini


!********************** solve **********************************

/SOLU                          

antype,stat,new 
!nlgeom,on
allsel
time,1
nsubst,20,100
autots,on        ! auto time stepping
pred,on          ! predictor on
nropt,full,,on   ! Newton-Raphson on
solve   
fini


!********************* results ********************************

/post1


path,one,2,,100            	! Define path name and parameters
ppath,1,,0,0,0                    ! Define path along bimorph length
ppath,2,,rs,0,0
pdef,Uy,u,y
/axlab,x, Position (m)
/axlab,y, Deflection (m)
plpath,Uy


/replot