/filename,vem_11
/title, vem_11

!only the quarter problem has been modelled

!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
	
!Geometry
	!Keypoints
	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
	
	!Lines
	L,4,5
	L,5,6
	L,6,7
	L,7,1
	L,3,4
	
	!Arcs
	LARC,1,3,2,r
	
	!Areas
	AL,1,2,3,4,5,6
	
!Meshing (free mesh)
	!Attributes
	ASEL,All
	AATT,1,1
	AESIZE,ALL,m1
	AMESH,ALL
		
!Loads
	!Displacement
	DL,4,,UY,0
	DL,3,,UX,0
	
	!Pressure
	SFL,2,PRES,-p0
		
/SOLU
Solve

/POST1
!Plot the deformed shape
	PLDISP,1
	!Scale the displacement
	!/DSCALE,All,5
	
!Stresses
	!Y-component of stress
		PLNSOL,S,Y
	!Stress intensity (Mohr)
		PLNSOL,S,INT
	!Equivalent stress (HMH)
		PLNSOL,S,EQV
	
!Path operations
	!Y-component of stress on Path
		PATH,km,2,pathsol,pathdiv
		PPATH,1,,0,0,0
		PPATH,2,,b,0,0
		PDEF,sy,S,Y
		PLPATH,sy				!On graph
		!PLPAGM,sy				!On geometry
		PRPATH,sy

		
!----------------Mesh 2-------------------------------------

/PREP7

!Delet mesh
	ACLEAR,1
	
!Mesh again
	!Attributes
	ASEL,All
	AATT,1,1
	AESIZE,ALL,m2
	AMESH,ALL
	
/SOLU
Solve

/POST1
!Plot the deformed shape
	PLDISP,1
	!Scale the displacement
	!/DSCALE,All,5
	
!Stresses
	!Y-component of stress
		PLNSOL,S,Y
	!Stress intensity (Mohr)
		PLNSOL,S,INT
	!Equivalent stress (HMH)
		PLNSOL,S,EQV
	
!Path operations
	!Y-component of stress on Path
		PATH,km,2,pathsol,pathdiv
		PPATH,1,,0,0,0
		PPATH,2,,b,0,0
		PDEF,sy,S,Y
		PLPATH,sy				!On graph
		!PLPAGM,sy				!On geometry
		PRPATH,sy
		
!----------------Mesh 3-------------------------------------

/PREP7

!Delet mesh
	ACLEAR,1
	
!Mesh again
	!Attributes
	ASEL,All
	AATT,1,1
	AESIZE,ALL,m3
	AMESH,ALL
	
/SOLU
Solve

/POST1
!Plot the deformed shape
	PLDISP,1
	!Scale the displacement
	!/DSCALE,All,5
	
!Stresses
	!Y-component of stress
		PLNSOL,S,Y
	!Stress intensity (Mohr)
		PLNSOL,S,INT
	!Equivalent stress (HMH)
		PLNSOL,S,EQV
	
!Path operations
	!Y-component of stress on Path
		PATH,km,2,pathsol,pathdiv
		PPATH,1,,0,0,0
		PPATH,2,,b,0,0
		PDEF,sy,S,Y
		PLPATH,sy				!On graph
		!PLPAGM,sy				!On geometry
		PRPATH,sy