/TITLE,CINT_2mats,CENTERED CRACK PLATE UNDER THERMAL LOADING WITH ELASTIC-PERFECTLY PLASTIC MATERIAL
/COM, REFERENCE : E.M.REMZI, NAFEMS, TWO DIMENSIONAL TEST CASES IN POST YIELD FRACTURE 
/COM,             MECHANICS, REF: R00038
/COM,             CHAPTER 2.7 FROM REFERENCE ARTICLE
/COM,
/PREP7
! COEFFICIENT OF THERMAL EXPANISON 

ET,1,PLANE182			! PLANE 292 for thermal planar ELEMENT
! KEYOPT,1,1,1
KEYOPT,1,3,2	         		! PLANE STRAIN


k,1,0,0,0
k,2,76,0,0
K,3,76,10,0
K,4,0,10,0
K,5,13,5,0
K,6,63,5,0
K,7,63,6.5,0
k,8,13,6.5,0

l,1,2
l,2,3
l,3,4
l,4,1
l,5,6
l,6,7
l,7,8
l,8,5

Al,1,2,3,4,5,6,7,8
Al,5,6,7,8
AGLUE,1,2

!epoxy outer THERMAL MATERIAL PROPERTIES
MP,dens,2,1.2E0-6
MP, c,2, 1100 E03   !1100 j/kg C
MP,EX,2,2.5E03  !2.5Gpa
MP,NUXY,2,0.36
MP,kxx,2, 0.00017! w/mm celsius 0.17w/m c
MP,alpx,2, 50E-06 ! per degree celsius

! TB,	PRONY,	2, , ,	EXPERIMENTAL
! TB, EXPE,	2, 7, , EMODULUS
! !***************************at -10 degree***************
! TBTEMP, -10 ! temperature
! ! TBPT, DEFI, frequency, storage modulus, loss modulus, tan delta
	! TBPT, DEFI,1E-5 , 2538 ,266 , ,

	! TBPT, DEFI, 1E-4, 2638, 261, ,

	! TBPT, DEFI, 1E-3, 2743, 256, ,

	! TBPT, DEFI, 0.010, 2856, 253, ,

	! TBPT, DEFI, 0.100, 2979, 250, ,

	! TBPT, DEFI, 1.000, 3124, 248, ,

	! TBPT, DEFI, 10.000, 3296, 246, ,
	
! !***********************************at 5 degree***********************	
	! TBTEMP, 5 ! temperature
! ! TBPT, DEFI, frequency, storage modulus, loss modulus, tan delta
	! TBPT, DEFI,1E-5 , 1707 ,274 , ,

	! TBPT, DEFI, 1E-4, 1906, 274, ,

	! TBPT, DEFI, 1E-3, 2069, 274, ,

	! TBPT, DEFI, 0.010, 2209, 274, ,

	! TBPT, DEFI, 0.100, 2329, 274, ,

	! TBPT, DEFI, 1.000, 2435, 274, ,

	! TBPT, DEFI, 10.000, 2535, 274, ,
		
	! TBPT, DEFI, 100.0, 2635, 273, ,
	! TBPT, DEFI, 1000.00, 2742, 271, ,
	! TBPT, DEFI, 10000.00, 2857, 268, ,
	! TBPT, DEFI, 100000.00, 2983, 264, ,
	! TBPT, DEFI, 1E06, 3132, 259, ,
	! TBPT, DEFI, 1E07, 3306, 255, ,
	
	! ! TBPT, DEFI, 1E08 , , 251, ,
	! ! TBPT, DEFI, 1E09, , 249, ,
	
	
	! !++++++++++++++++++ AT 20 deg*********************
	! TBTEMP, 20 ! temperature
! ! TBPT, DEFI, frequency, storage modulus, loss modulus, tan delta
	! TBPT, DEFI,1E-5 , 372 ,249 , ,

	! TBPT, DEFI, 1E-4, 551, 264, ,

	! TBPT, DEFI, 1E-3, 799, 272, ,

	! TBPT, DEFI, 0.010, 1092, 277, ,

	! TBPT, DEFI, 0.100, 1384, 278, ,

	! TBPT, DEFI, 1.000, 1639, 277, ,

	! TBPT, DEFI, 10.000, 1851, 276, ,
		
	! TBPT, DEFI, 100.0, 2023, 275, ,
	! TBPT, DEFI, 1000.00, 2169, 274, ,
	! TBPT, DEFI, 10000.00, 2296, 274, ,
	! TBPT, DEFI, 100000.00, 2405, 274, ,
	! TBPT, DEFI, 1E06, 2507, 274, ,
	! TBPT, DEFI, 1E07, 2608, 274, ,
	
	! ! TBPT, DEFI, 1E08 , , 274, ,
	! ! TBPT, DEFI, 1E09, , 274, ,
	! ! TBPT, DEFI, 1E010, , 272, ,
	! ! TBPT, DEFI, 1E011, , 270, ,
		
	! !************************** at 35 degree***************
	! TBTEMP, 35 ! temperature
! ! TBPT, DEFI, frequency, storage modulus, loss modulus, tan delta
	! TBPT, DEFI,1E-5 , 142 ,64 , ,

	! TBPT, DEFI, 1E-4, 149, 85, ,

	! TBPT, DEFI, 1E-3, 166, 110, ,

	! TBPT, DEFI, 0.010, 204, 140, ,

	! TBPT, DEFI, 0.100, 276, 172, ,

	! TBPT, DEFI, 1.000, 402, 203, ,

	! TBPT, DEFI, 10.000, 596, 231, ,
		
	! TBPT, DEFI, 100.0, 858, 252, ,
	! TBPT, DEFI, 1000.00, 1156, 265, ,
	! TBPT, DEFI, 10000.00, 1442, 273, ,
	! TBPT, DEFI, 100000.00, 1689, 277, ,
	! TBPT, DEFI, 1E06, 1891, 278, ,
	! TBPT, DEFI, 1E07, 2057, 277, ,
	
	! TBPT, DEFI, 1E08 ,2199 , 276, ,
	! TBPT, DEFI, 1E09, 2321, 275, ,
	! TBPT, DEFI, 1E010, 2427, 274, ,
	
	! ! TBPT, DEFI, 1E011, , , ,
	
		! !************************** at 50  degree***************
	! TBTEMP, 50 ! temperature
! ! TBPT, DEFI, frequency, storage modulus, loss modulus, tan delta
	! TBPT, DEFI,1E-5 , 130 ,18 , ,

	! TBPT, DEFI, 1E-4, 133, 21, ,

	! TBPT, DEFI, 1E-3, 136, 24, ,

	! TBPT, DEFI, 0.010, 138, 29, ,

	! TBPT, DEFI, 0.100, 141, 37, ,

	! TBPT, DEFI, 1.000, 146, 48, ,

	! TBPT, DEFI, 10.000, 159, 63, ,
		
	! TBPT, DEFI, 100.0, 186, 84, ,
	! TBPT, DEFI, 1000.00, 239, 109, ,
	! TBPT, DEFI, 10000.00, 334, 138, ,
	! TBPT, DEFI, 100000.00, 487, 171, ,
	! TBPT, DEFI, 1E06, 706, 203, ,
	! TBPT, DEFI, 1E07, 982, 230, ,
	
	! TBPT, DEFI, 1E08 ,1271 , 251, ,
	! TBPT, DEFI, 1E09, 1536, 265, ,
	! TBPT, DEFI, 1E010, 1761, 273, ,
	
	! ! TBPT, DEFI, 1E011, 1945, , ,
	! ! TBPT, DEFI, 1E012, 2099, , ,


! !************************** at 65  degree***************
	! TBTEMP, 65 ! temperature
! ! TBPT, DEFI, frequency, storage modulus, loss modulus, tan delta
	! ! TBPT, DEFI,1E-5 ,  ,13 , ,

	! ! TBPT, DEFI, 1E-4, , 14, ,

	! TBPT, DEFI, 1E-3, 119, 13, ,

	! TBPT, DEFI, 0.010, 128, 13, ,

	! TBPT, DEFI, 0.100, 130, 14, ,

	! TBPT, DEFI, 1.000, 133, 17, ,

	! TBPT, DEFI, 10.000, 135, 20, ,
		
	! TBPT, DEFI, 100.0, 138, 23, ,
	! TBPT, DEFI, 1000.00, 141, 27, ,
	! TBPT, DEFI, 10000.00, 145, 33, ,
	! TBPT, DEFI, 100000.00, 156, 43, ,
	! TBPT, DEFI, 1E06, 181, 56, ,
	! TBPT, DEFI, 1E07, 231, 74, ,
	
	! TBPT, DEFI, 1E08 ,321 , 97, ,
	! TBPT, DEFI, 1E09, 468, 125, ,
	! TBPT, DEFI, 1E010, 650, 156, ,
	! TBPT, DEFI, 1E011, 952, 189, ,
	
	! ! TBPT, DEFI, 1E012, 1243, , ,
	
	! !************************** at 80  degree***************
	! TBTEMP, 80 ! temperature
! ! TBPT, DEFI, frequency, storage modulus, loss modulus, tan delta
	! ! TBPT, DEFI,1E-5 , 130 ,18 , ,

	! ! TBPT, DEFI, 1E-4, , 14, ,

	! ! TBPT, DEFI, 1E-3, 119, 13, ,

	! ! TBPT, DEFI, 0.010, 128, 13, ,

	! ! TBPT, DEFI, 0.100, 130, 14, ,

	! ! TBPT, DEFI, 1.000, 133, 17, ,

	! TBPT, DEFI, 10.000, 121, 14, ,
		
	! TBPT, DEFI, 100.0, 128, 13, ,
	! TBPT, DEFI, 1000.00, 130, 13, ,
	! TBPT, DEFI, 10000.00, 133, 15, ,
	! TBPT, DEFI, 100000.00, 136, 17, ,
	! TBPT, DEFI, 1E06, 138, 20, ,
	! TBPT, DEFI, 1E07, 141, 23, ,
	
	! TBPT, DEFI, 1E08 ,145 , 28, ,
	! TBPT, DEFI, 1E09, 156, 34, ,
	! TBPT, DEFI, 1E010, 180, 44, ,
	
	! ! TBPT, DEFI, 1E011, , 59, ,
	
	! ! TBPT, DEFI, 1E012, ,77 , ,
	! ! TBPT, DEFI, 1E012, ,101 , ,
	! TB,SHIFT,2,1,1,FICT
! TBDATA, 1,    130,      0.0,     1.0 ! SHIFT FUNCTION PARAMETERS
! TBDATA, 4,    130
! ,      1.0,     0.0 ! FICTIVE TEMP PARAMETERS, TFI, CFI, TAUFI
! TBDATA, 7,52.6E-7, 0.119E-7, -1.0E-11 ! GLASS CTE COEFFICIENTS
! TBDATA,12,52.6E-7, 0.119E-7, -1.0E-11 ! LIQUID CTE COEFFICIENTS


ESIZE,0.5	
! KSCON,9,0.1,1,8,,		
! kscon,5,0.08,1,8,
xnume = 601
ynume = 83
lsel,s,line, ,1,3,2
LESIZE,all, , ,xnume, , , , , 1
lsel,s,line, ,2,4,2
LESIZE,all, , ,ynume, , , , , 1
type,1
mat, 2
! MSHKEY,1 ! 
AMESH,1 !area 1 is the outer one
allsel,all

!stainless steel AISI201
MP, EX , 1, 197E03
MP, NUXY, 1, 0.265
MP, kxx, 1, 0.0605 !w/mm celsius
MP, alpx, 1, 18E-06 !/k
MP, dens, 1, 7.8E0-6
MP, c, 1, 500E03 j/kg celsius

ESIZE,0.2
! MSHKEY,0! ENFORCE MAPPED MESHING
! MSHAPE,0,2D ! ENFORCE QUADRILATERAL ELEMENT
type,1
MAT,1
AMESH,2 !area 2 is the inner one
ALLSEL,ALL

!*********************************************+XFENRICH*****************************************
! Element component required for XFENRICH command
nsel, s, loc, y, 0,10
nsel, r, loc, x,  0,38
esln
cm, testcmp, elem
allsel

! mesh the crack surface with mesh200 elements 
et, 2, 200, 0 ! keyopt(1) = 0 for mesh200 line elements
       
k, 10, 0.0, 5, 0.0
k, 11,76/601*102,5 , 0.0   
l,10, 11! define a line 

type, 2 
mat, 2  
lesize,9,,,12
lmesh,9
allsel  

! element component for mesh200 elements   
esel,s,type,, 2 
cm, m200el, elem   
allsel  
     
! mesh200 node component for crk front   
 nsel,s,loc,x,76/601*102
 nsel,r,loc,y,5.0
! nsel,a,loc,x,13.0
! nsel,r,loc,y,5.0
! cmsel,s,m200el,elem 
! nsle,r,all,
! nsel,s,loc,x,13
! nsel,r,loc,y,5
 nlist
cm, m200nd, node
allsel  

 ! Define enrichment identification
xfenrich, ENRICH1, testcmp, , SING,1.5,0.01
allsel

! Define LSM values for cut elements
xfcrkmesh, ENRICH1, m200el, m200nd
     
/com
/com *************************************************************  
/com
/com  		  INITIAL CRACK DATA  
/com
/com *************************************************************  
/com
xflist  
********************************************xfem end*********************

NSEL,S,LOC,X,13,63 !bottom side
NSEL,R,LOC,Y,5
D, ALL, All
allsel,all

NSEL,S,LOC,X,13,63 !top side
NSEL,R,LOC,Y,6.5
D, ALL, All
allsel,all

NSEL,S,LOC,X,13 !left side
NSEL,R,LOC,Y,5,6.5
D, ALL, All
allsel,all

NSEL,S,LOC,X,63
NSEL,R,LOC,Y,5,6.5
D, ALL, All
allsel,all

! condtraints set 2 on the epoxy
NSEL,S,LOC,X,0!bottom side
NSEL,R,LOC,Y,0
D, ALL, All
allsel,all

NSEL,S,LOC,X,0,76!bottom side
NSEL,R,LOC,Y,0
D, ALL, UY,0
allsel,all

finish

/SOLU

a= 130
b= -30
ANTYPE,static
!trnopt,all
!lumpm,0
OUTRES,CINT,All
NLGEOM,ON
autots,on
allsel,all
!NSUBS,20,1000,10
tref,30
*DO,I,1,2
 *IF,I,EQ,1,THEN 
 
		NSEL,S,LOC,x,0,76
		NSEL,R,LOC,y,0
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0
		NSEL,R,LOC,y,0,10
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0,76
		NSEL,R,LOC,y,10
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,76
		NSEL,r,LOC,y,0,10
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all

!ldread,temp,,,1200,0,'id2_ther_01','rth',' '
               	! auto time stepping
	allsel,all
	
	

	NSEL,S,LOC,Y,5
	NSEL,R,LOC,X,13 !node location of th node
	CM,CRACK1,NODE			! CRACK TIP NODE COMPONENT
	ALLSEL,ALL


	CINT,NEW,1
	CINT,NAME,CRACK1
	CINT,NCON,8		! NUMBER OF COUNTOURS AROUND CRACK TIP
	CINT,SYMM,Off			! SYMMETRY TURNED OFF
	CINT,NORM,0,
	ALLSEL,all
! **********************XFENRICH crackfront element
	!Fracture Parameter calculations  
	CINT, NEW, 2
	CINT, CXFE,  _XFCRKFREL1
	 CINT, TYPE, jint, 2 
	CINT, NCON, 8     
	CINT, NORM, 0, 2
		
	
	!CGROW calculations
	! cgrow, new, 2
	! cgrow, cid, 2
	! cgrow, method, xfem
	! cgrow, fcoption, mtab, 2
	
	ttime=1800*I			! Total simulated time, in seconds
	dtime=60				! Initial time increment
	dtimemax=60				! Max time increment
	kbc,0
	outres,all,-10			! Store results at 10 equal substeps
	time,ttime				! Total time 600 sec
	deltim,dtime,dtime,dtimemax
	inc = I
	lswrite,inc

	/solu
	! ANTYPE,static
	!trnopt,all
	!lumpm,0
	OUTRES,CINT,All
	!NLGEOM,ON
	autots,on
	allsel,all
	!NSUBS,20,1000,10
	
	!ldread,temp,,,1200,0,'id2_ther_01','rth',' '
					! auto time stepping
	allsel,all
	!NSEL,S,LOC,x,0,76
		NSEL,R,LOC,y,0
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0
		NSEL,R,LOC,y,0,10
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0,76
		NSEL,R,LOC,y,10
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,76
		NSEL,r,LOC,y,0,10
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all
	! !Fracture Parameter calculations  
		! CINT, NEW, 2
		CINT, CXFE,  _XFCRKFREL1
		 CINT, TYPE, jint, 2 
		CINT, NCON, 8     
		CINT, NORM, 0, 2
	
	ttime=3600*I		! Total simulated time, in seconds
	cycle_time=ttime !for inputing in the cycle 2 to n loop for time increement
	dtime=60				! Initial time increment
	dtimemax=60				! Max time increment
	kbc,0
	outres,all,-10			! Store results at 10 equal substeps
	time,ttime				! Total time 3600 sec
	deltim,dtime,dtime,dtimemax
	inc = inc + 1			!for incrementing the load step file number
	lswrite,inc

		NSEL,S,LOC,Y,5
		NSEL,R,LOC,X,13 !node location of th node
		CM,CRACK1,NODE			! CRACK TIP NODE COMPONENT
		ALLSEL,ALL

		! CINT,NEW,1
		CINT,NAME,CRACK1
		CINT,NCON,8	! NUMBER OF COUNTOURS AROUND CRACK TIP
		CINT,SYMM,Off		! SYMMETRY TURNED OFF
		CINT,NORM,0,
		ALLSEL,all
	! **********************XFENRICH crackfront element
		
			
		
		!CGROW calculations
		! cgrow, new, 2
		! cgrow, cid, 2
		! cgrow, method, xfem
		! cgrow, fcoption, mtab, 2
	
	*ELSE
	 allsel,all
	 /solu
	 OUTRES,CINT,All
	!NLGEOM,ON
	autots,on
	allsel,all
		NSEL,S,LOC,x,0,76
		NSEL,R,LOC,y,0
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0
		NSEL,R,LOC,y,0,10
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0,76
		NSEL,R,LOC,y,10
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,76
		NSEL,r,LOC,y,0,10
		bf,all,temp,a
		sf,all,conv,3.4e-6,a
		allsel,all
!ldread,temp,,,1200,0,'id2_ther_01','rth',' '
               	! auto time stepping
	allsel,all
	
	! !Fracture Parameter calculations  
	! ! CINT, NEW, 2
	CINT, CXFE,  _XFCRKFREL1
	 CINT, TYPE, jint, 2 
	CINT, NCON, 8     
	CINT, NORM, 0, 2
	
	ttime=cycle_time + 1800		! Total simulated time, in seconds
		dtime=60				! Initial time increment
		dtimemax=60				! Max time increment
		kbc,0
		outres,all,-10			! Store results at 10 equal substeps
		time,ttime				! Total time 600 sec
		deltim,dtime,dtime,dtimemax
		inc = inc +1 			!for increementing the load step file number
		lswrite,inc
		
	NSEL,S,LOC,Y,5
	NSEL,R,LOC,X,13 !node location of th node
	CM,CRACK1,NODE			! CRACK TIP NODE COMPONENT
	ALLSEL,ALL

	! CINT,NEW,1
	CINT,NAME,CRACK1
	CINT,NCON,8	! NUMBER OF COUNTOURS AROUND CRACK TIP
	CINT,SYMM,Off			! SYMMETRY TURNED OFF
	CINT,NORM,0,
	ALLSEL,all
		
	
	!CGROW calculations
	! cgrow, new, 2
	! cgrow, cid, 2
	! cgrow, method, xfem
	! cgrow, fcoption, mtab, 2
	
	/solu
	! ANTYPE,static
	!trnopt,all
	!lumpm,0
	OUTRES,CINT,All
	!NLGEOM,ON
	autots,on
	allsel,all
	!NSUBS,20,1000,10
		NSEL,R,LOC,y,0
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0
		NSEL,R,LOC,y,0,10
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,0,76
		NSEL,R,LOC,y,10
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all

		NSEL,S,LOC,x,76
		NSEL,r,LOC,y,0,10
		bf,all,temp,b
		sf,all,conv,3.4e-6,a
		allsel,all
	!ldread,temp,,,1200,0,'id2_ther_01','rth',' '
					! auto time stepping
	allsel,all
	
	
	!Fracture Parameter calculations  
		! CINT, NEW, 2
		CINT, CXFE,  _XFCRKFREL1
		 CINT, TYPE, jint, 2 
		CINT, NCON, 8     
		CINT, NORM, 0, 2
	!
	ttime=3600*I		! Total simulated time, in seconds
		cycle_time = ttime !for the first load step time in the next cycle
		dtime=60				! Initial time increment
		dtimemax=60				! Max time increment
		kbc,0
		outres,all,-10			! Store results at 10 equal substeps
		time,ttime				! Total time 3600 sec
		deltim,dtime,dtime,dtimemax
		inc = inc +1
		lswrite,inc

		NSEL,S,LOC,Y,5
		NSEL,R,LOC,X,13 !node location of th node
		CM,CRACK1,NODE			! CRACK TIP NODE COMPONENT
		ALLSEL,ALL

		! CINT,NEW,1
		CINT,NAME,CRACK1
		CINT,NCON,8	! NUMBER OF COUNTOURS AROUND CRACK TIP
		CINT,SYMM,Off		! SYMMETRY TURNED OFF
		CINT,NORM,0,
		ALLSEL,all
	! **********************XFENRICH crackfront element
		
			
		
		!CGROW calculations
		! cgrow, new, 2
		! cgrow, cid, 2
		! cgrow, method, xfem
		! cgrow, fcoption, mtab, 2
	*endif
*enddo
	
lssolve,1,inc,1

finish


/OUT,SCRATCH
/POST1
SET,LAST
! PRCINT,2
! *GET,MAXNODE,CINT,2,NNOD

! *GET,tipnode2,CINT,2,NODE,1

! *GET,J1,CINT,2,CTIP,tipnode2,,1,,
! *GET,J2,CINT,2,CTIP,tipnode2,,2,,
! *GET,J3,CINT,2,CTIP,tipnode2,,3,,
! *GET,J4,CINT,2,CTIP,tipnode2,,4,,
! *GET,J5,CINT,2,CTIP,tipnode2,,5,,
! *GET,J6,CINT,2,CTIP,tipnode2,,6,,
! *GET,J7,CINT,2,CTIP,tipnode2,,7,,
! *GET,J8,CINT,2,CTIP,tipnode2,,8,,

! JC2 = (ABS(J2)+ABS(J3)+ABS(J4)+ABS(J5)+ABS(J6)+ABS(J7)+ABS(J8))/7
! *STAT,JC2
! *DIM,LABEL1,CHAR,1
! *DIM,VALUE1,,1,3!assume 11
! LABEL1(1,1) = 'J-VALUE'
! *VFILL,VALUE1(1,1),DATA,55
! *VFILL,VALUE1(1,2),DATA,JC2
! *VFILL,VALUE1(1,3),DATA,ABS(JC2/55)

PRCINT,1
*GET,MAXNODE1,CINT,1,NNOD !****************for crack 1

*GET,tipnode1,CINT,1,NODE,1

*GET,J1,CINT,1,CTIP,tipnode1,,1,,
*GET,J2,CINT,1,CTIP,tipnode1,,2,,
*GET,J3,CINT,1,CTIP,tipnode1,,3,,
*GET,J4,CINT,1,CTIP,tipnode1,,4,,
*GET,J5,CINT,1,CTIP,tipnode1,,5,,
*GET,J6,CINT,1,CTIP,tipnode1,,6,,
*GET,J7,CINT,1,CTIP,tipnode1,,7,,
*GET,J8,CINT,1,CTIP,tipnode1,,8,,

JC2 = (ABS(J2)+ABS(J3)+ABS(J4)+ABS(J5)+ABS(J6)+ABS(J7)+ABS(J8))/7
*STAT,JC2
*DIM,LABEL1,CHAR,1
*DIM,VALUE1,,1,3!assume 11
LABEL1(1,1) = 'J-VALUE'
*VFILL,VALUE1(1,1),DATA,55
*VFILL,VALUE1(1,2),DATA,JC2
*VFILL,VALUE1(1,3),DATA,ABS(JC2/55)





/COM
PRCINT,1,,JINT
/COM
/COM >>> MODE 1 STRESS INTENSITY FACTOR
/COM
PRCINT,1,,K1
/COM
/COM >>> MODE 2 STRESS INTENSITY FACTOR
/COM
PRCINT,1,,K2
/COM






/OUT,vmr038-2g-182,vrt
/COM, -----------------vmr038-2g-182 RESULTS COMPARISON -----------------
/COM,
/COM,           |  TARGET  |   Mechanical APDL  |  RATIO
/COM,
/COM, PLAIN STRAIN WITH LOAD CONTROL
/COM, --------------------------------
/COM,
*VWRITE,LABEL1(1,1),VALUE1(1,1),VALUE1(1,2),VALUE1(1,3)
(1X,A8,'   ',F10.4,'  ',F14.4,'   ',F15.3)
/COM,
/COM,
/COM,------------------------------------------------------------------------
/OUT,
FINISH
*LIST,vmr038-2g-182,vrt