! Steuerdatei
fini
/clear,start
/TRIAD,LBOT 
/REPLOT 
/config,nproc,2
/FILNAME,VI
! Körperdemensionen in mm

/COM,ANSYS RELEASE  7.0    UP20021010       16:35:31    10/29/2004
!RESUME,mat,db,G:\users\nguyenc\FIB	!(Materialdatei laden)   
/COM,ANSYS RELEASE  7.0    UP20021010       16:35:31    10/29/2004
/prep7

!ANTYPE,STATIC                 ! STATIC ANALYSIS
!nlgeom,on

!deltim,0.1,0.001,0.2,


N_V = 1                       ! Version
JOBNAME =  'KS_%N_V%'         !
/TITLE, 21_07_05-2,hardness measuring  



!****************************************************************
	ET,1,PLANE182	
			
!Definition des elatischen Verhaltens

	MPTEMP,,,,,,,,  
	MPTEMP,1,0  
	MPDATA,EX,1,,80200
	MPDATA,PRXY,1,,0.3 
	TB,BISO,1                        
	TBDATA,1,185,1000
	eps=10e-8

!Elemanttypem definieren



!KOntaktelement
!	ET,3,CONTA171	

!KOntaktelement
!	ET,4,CONTA171	

!Targetelement
!	ET,5,TARGE169
!****************************************************************
!Definition der Varaiable

	Breite= 10e-3	!Globalmodel in Y-Richtung
	Höhe_1=300e-6	!Dicke der Schicht 1 in X-Richtung
	Höhe_2=10e-3	!Dicke der Schicht 2 in X-Richtung
	D_Höhe=600e-6	!Höhe des Diamantes
	distance =10e-6	!distance between substrate and diamond
	eps=1e-8	!tolerance
	Loch_B=4e-4	!Breite des Loches.


!***************************************************************

!generate the substrate volume 

	RECTNG,0,Breite,Höhe_2,Höhe_1+Höhe_2,	!Abmessung
	RECTNG,0,Breite,0,Höhe_1+Höhe_2,		!Abmessung
	RECTNG,0,Loch_B,0,Höhe_1+Höhe_2,

! area of diamond
	!RPR4,3,10e-3,10e-3,3e-3,270 


!Generate Keypoints for diamond area
	
	k,100,0,(Höhe_1+Höhe_2+distance),0
	k,102,0-(D_Höhe/tan(65.3)),(Höhe_1+Höhe_2+distance)+D_Höhe,0

!creat lines from Keypoints

	L,102,100

!choose all the created lines

	Lsel,s,loc,y,(Höhe_1+Höhe_2+10e-6),10

!create area from created lines
	!Al,all,
	!Alls
!overlap alls of areas
	aovlap,all  



!***************************************************************

!Material assignation

!choose the areas to reassign material data

	Asel,s,loc,y,(Höhe_1+Höhe_2+10e-6),10
	AATT,1,,,0 
	ALLS


	Asel,s,loc,y,(Höhe_1+Höhe_2+10e-6),Höhe_2
	AATT,1,,,0   
	ALLS


	Asel,s,loc,y,0,Höhe_2-eps
	AATT,1,,,0   
	ALLS

!***************************************************************

!meshing

	MSHAPE, 0, 2D
	mshkey,1


	lsel,s,loc,y,eps,Höhe_2-eps,
	lsel,r,loc,x,0+eps,Breite-eps,
	LESIZE,ALL,,,20,1/20

	lsel,s,loc,y,eps,Höhe_2-eps,
	lsel,r,loc,x,0+eps,0-eps,
	LESIZE,ALL,,,20,20

	lsel,s,loc,y,eps,Höhe_2-eps,
	lsel,r,loc,x,Breite-eps,Breite+eps,
	LESIZE,ALL,,,20,20

	ALLS

	lsel,s,loc,y,Höhe_2+eps,Höhe_1+Höhe_2-eps
	LESIZE,ALL,,,8
	ALLS

	
	lsel,s,loc,x,Loch_B+eps,Breite-eps,
	LESIZE,ALL,,,20,1/20
	ALLS
 	

	asel,s,loc,x,0,Loch_B+eps
	esize,5e-5
	AMESH,all  
	ALLS
	
		
	asel,s,loc,x,Loch_B+eps,Breite-eps,
	esize
	AMESH,all  
		
	ALLS

	n,400000,Breite/2,-1.5e-3,
!*************************************************************

/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,1,
MAT,1   
R,3 
REAL,3  
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,1   
R,3,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
LSEL,S,,,13 
CM,_TARGET,LINE 
TYPE,2  
LATT,-1,3,2,-1  
TYPE,2  
LMESH,ALL   
! Generate the contact surface  
LSEL,S,,,15 
LSEL,A,,,18 
CM,_CONTACT,LINE
TYPE,3  
NSLL,S,1
ESLN,S,0
ESURF   
*SET,_REALID,3  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
LSEL,S,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,2  
ESEL,A,TYPE,,3  
ESEL,R,REAL,,3  
LSEL,S,REAL,,3  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
/MREP,EPLOT 


ESEL,S,TYPE , , 2 
ESURF,,REVE 
/REPLOT 
ALLS
EPLOT   
/PSYMB,ESYS,1   
/REPLOT 

/PSYMB,ESYS,1
!**************************************************************


	


!**************************************************************


/COM, CONTACT PAIR PROPERTIES - START   
!Create a pilot node for contact pair ID 3  
REAL,3  
TYPE,2  
TSHAP,PILO  
E,400000
/COM, CONTACT PAIR PROPERTIES - END 
/MREP,EPLOT 



lsel,s,loc,y,0-eps,0+eps
DL,all, ,UY,0 
ALLS

nsel,s,loc,y,0-eps,0+eps
nsel,r,loc,x,0-eps,0+eps
D,all,UX,0 
alls

nsel,s,loc,x,0-eps,0+eps
D,all,UX,0 
alls

!**************************************************************
!Contact definition


R,3,0,1e-4,0.2, 0,1,0  
RMORE,0,0,0, 0,0,0  
RMORE,0,,, ,,   
RMORE,,,0, 0,0, 


!KOntakt Eigenschaft

!KEYOPT,1,1,0
!KEYOPT,1,3,1
!KEYOPT,1,6,0
!KEYOPT,1,10,0   
!KEYOPT,3,1,0
!KEYOPT,3,2,0
!KEYOPT,3,3,0
!KEYOPT,3,4,0
!KEYOPT,3,5,1
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,1   
KEYOPT,3,11,0   
KEYOPT,3,12,0   



!**************************************************************
save
!Lösung
/solu



	antype,0,new
	NLGEO,on
	SOLCON,ON
	NEQIT,50
	eqslv,
	EQSLV,FRont
	!NSUB,10,200
	deltim,0.1,0.02,0.1
	time,1
	OUTRES,ALL,all
	ERESX, NO
        kbc,0			! ramped load
	!NROPT, FULL , --,OFF
	Nsel,s,NOde,,400000
	
	D,400000,Ux,0
	
	D,400000,Uy,-477e-6
	
	allsel
	
	solve


	NLGEO,on
	SOLCON,ON
	NEQIT,50
	eqslv,
	EQSLV,FRont
	!NSUB,100,200
	deltim,0.005,0.001,0.2
	time,2
	OUTRES,ALL,all
	ERESX, NO
        kbc,0			! ramped load
	!NROPT, FULL , --,OFF
	Nsel,s,NOde,,400000
	
	D,400000,Ux,0
	D,400000,Uy,0
	
	allsel
	solve




/post26
!reaction Force : FY_3

RFORCE,2,400000,F,Y, FY_2   

!Displacement : UY_3

NSOL,3,400000,U,Y, UY_3 

!Defines an array parameter and its dimensions

*DIM,FY_2,ARRAY,50,1,1, , ,
*DIM,UY_3,ARRAY,50,1,1, , ,

VGET,FY_2,2,0,0
VGET,UY_3,3,0,0

*cfopen,RESULTS,dat

!*vwrite
!('rforce: displacment:')

*vwrite,	FY_2(1,1),-UY_3(1,1)
(E11.4,E11.4)

*cfclos



!****************************************************************
!3 D Darstellung
!*  
/EXPAND,18,AXIS,,,10
/REPLOT 




!****************************************************************