PROGRAM SMLDI A( INPUT ,OUTPUT .PUNCH ) 



DIMENSION X(ai),Y(4l),Z(41).YS(Al),VI(41I.FI(Al),FP(«l).F(41), 

 1FV(41),RES<41),G(I31),P(41),0(41),HI(4),HX(4) 

 SET TEST CONDITIONS 

 ANGLE=0. 

 D'Z.COO 

 PEAC IN DIGITIZED DATA 

 8 BEAD 1 ,UP,DN,UF ,0F ,0 I A,XC 

 tF(UP) 10,10,9 

 9 READ 2,CL ,T ,N, X«,XF,C,«a,FO, (HI < I ) , I =1 ,4> 

 READ 3, (FI ( I ) ,1 =1 ,40> 

 1 FORMAT { 2F3. 3, 2F j,0, 2F3,3 » 



2 FORMAT(F3.3,F3.2,I2,2F2,l,7F3.0) 



3 F0RMAT(24F3.0) 

 DETERMINE *»VE HEIGHT 



00 11 1=1,4 



IF(HI ( I )-«0)21 ,23,23 

 21 HX( I >=-(*0-HI ( I ) )•( ON/C)«( 1 ./XW> 



GO TO 11 

 23 ^iX(I)«(HI(I)-WO)*(UP/C)*(l./X*) 

 I 1 CONTI N IE 



H« (HX( 1 ) fMX(2)-MX{3)-HX(4) )/2. 

 CALCULATE CONSTANTS IN FORCE EQUATION 



PI=3. 1415926E36 



B«1.93B 



CALL «AV£L( T,0,XL) 



ZV»CL*(.5*0IA) 



A=EXP<(2.*PI*ZV)/XL) 



COSHA=At( 1. /A) 



SINHA=A-I1./A) 



B=EXP( (2.*PI*0)/XL) 



SINHa=B-( 1»/B) 



CS=C0SHA/SINHR 



SS=3INHA/SINHB 



CF 1=R*0I A»XC*H*PI*PI/(2.*T*T) 



FLV=CF1*H*CS*CS 



FDV«CF1»H»SS*SS 



FMV=CF1*DIA*PI*SS 



U«(H*CS*PI )/T 

 DETERMINE AVERAGE, MAXIMUM, AND MINIMUM FORCES 



FMAX=-.5 



FMINa.5 



SF=0. 



00 13 1=1,40 



IFIFK I )-FO)70,71,71 



70 FP< I)=(FI(I )-FO»*<DF/C)*.00 220 4/XF 

 GO TO 7 3 



71 FP(I»=(FI(I)-FO)*<UF/C)*.00 2204/XF 

 73 CONTINUE 



IF(FP( I ).GT,FMAX)902,90 1 



902 FMAX=FP(I) 



901 IF(FP( I ).LT.FMIN)903,12 



903 FMIN=FP( I ) 



12 SF=SF+FP(I) 



13 CONTINUE 



■ SF=SF/40. 



41 



