/filename,vem_11
/title, vem_11

!Data
	r=12.5
	a=400/2
	b=200/2
	t=10
	
	young=200E3
	nu=0.3
	p0=50
	
	m1=10
	m2=m1/2
	m3=m2/2
	pathsol=30
	pathdiv=20
	
/PNUM,KP,1  	
/PNUM,LINE,1
/PNUM,arEa,1

/PREP7
	
!Element type
	ET,1,PLANE182,,,3
	R,1,t
	
!Material data
	MPTEMP,1,0 
	MPDATA,EX,1,,young
	MPDATA,PRXY,1,,nu


*DO,I,0,3,1

	r=100/(2**I)

	/PREP7
	ACLEAR,ALL 
	ADELE,ALL,,,1   
	KDELE,ALL  
	
	*IF,I,eq,0,then
		!Geometry
		K,1,b,0
		K,2,b+r,0
		K,4,a,r
		K,5,a,a
		K,6,0,a
		K,7,0,0
		L,4,5
		L,5,6
		L,6,7
		L,7,1
		LARC,1,4,2,r
		AL,1,2,3,4,5
	*ELSE	
		!Geometry
		K,1,b,0
		K,2,b+r,0
		K,3,b+r,r
		K,4,a,r
		K,5,a,a
		K,6,0,a
		K,7,0,0
		L,4,5
		L,5,6
		L,6,7
		L,7,1
		L,3,4
		LARC,1,3,2,r
		AL,1,2,3,4,5,6
	*ENDIF
	
	!Meshing (free mesh)
	ASEL,All
	AATT,1,1
	AESIZE,ALL,m1
	AMESH,ALL
	!Loads
	DL,4,,UY,0
	DL,3,,UX,0
	SFL,2,PRES,-p0
			
	/SOLU
	Solve

	/POST1

	!Stress intensity (Mohr)
	PLNSOL,S,INT
	/IMAGE,SAVE,CHRVAL(r),BMP

*ENDDO