!Macro used to retrieve deflection along the center of a elliptical diaphragm 
!for t = 150 um and L varying under 350 psi pressure loading.

/clear
/title, Deflection

! L is changing

ni = 2				!number of length parameter variatins
nj = 1				!number of thickness parameter variations
pp = 1				!number of pressure loadings simulated
				 
*dim,L_,array,ni		!nnx1 vector of possible diaphragm lengths
*dim,h_,array,nj		!nnx1 vector of possible thicknesses 
*dim,ccc,char,ni+2

!Configurations which yield desired strain for 350psi
L_(1)= 3e-3,4e-3	!fill length array (m)
h_(1)=150e-6	!fill thickness array (m)
ccc(1)='P1','P2','P3','P4'

*dim,disp_F,table,ni,nj,pp	!Create Table nnxnnxpp of displacement 
                            !results for fixed BC
*dim,strainx_F,table,ni,nj,pp	!create table nnxnnxpp of axial  
                              !strain at center point
*dim,strainy_F,table,ni,nj,pp	!create table nnxnnxpp of Y-direction 
                              !strain at center point

*dim,nnum,array,90

*do,i,1,ni
  *do,j,1,nj
	path,L_(i),h_(j),pp,i,j !calls macro path
	*enddo
 
*if,i,eq,ni,exit
  
/sol 
 lsclear,all		!clears loads and load step options from database
 finish		
 /prep7
 aclear,1		!Deletes nodes and area elements associated with area
 adele,all,1		!clear area and below
 ldele,all,1
 kdele,all		!deletes all unmeshed keypoints
 Finish

*enddo