500 DDELP=180.+( ( 180. /PI )*DELP) 



GO TO 80 

 600 DDELP=360.+ ( { 180. /PI >*DELP> 

 80 PRINT 6* Z , DDELA 



6 FORMAT ( 1H1 * 3H Z.= , F5 . 1 , 1 1 X , 1 3HDEL T A AC T I VE = , F 6 . 1 ) 

 PRINT 7 



7 FORMAT(1HO///,19X,4HSIG1,13X,2HT1,13X,5HS IG0A,11X,5HSIG2A*12X, 

 13HT2A) 



PRINT 21 i SIG1»T1,SIG0A,SIG2A»T2A 

 21 FORMAT (9X,5F16.1) 

 PRINT 23 



23 FORMAT(1HO*25X,5HALPHA,21X,5HSIG3A,12X,3HT3A,13X,3HP3A) 

 DO 24 I A L PHA = N»360,N 



AL PH AD= I ALPHA 



C ALPHAD=ANGLE OF PLANF FOR WHICH STRESS CONDITIONS ARE TO BE FOUND 

 IF ( ALPHAD-DDFLA ) 15*15*16 

 16 AL PHAD=ALPHAD+ 1 80 . 



I F ( ALPHAD-360 . ) 1 5 *15*92 

 9 2 GO TO 1O5 



15 ALPHA =(ALPHAD*PI )/ldO. 



T AN A=S INF ( ALPHA) AC OSF( ALPHA ) 



TANA2=TANA*TANA 



IF (ALPHAD-90.) 44*40,44 

 44 IF ( ALPHAD-2 70 . ) 48,40,48 



4 8 S IG3A= ( ( 2.*S IG0A-2 .*T2A*TANA + 2 .*S IG2A*TANA2 ) / ( 1 . + TANA2 ) ) ~ S I G 2 A 

 C S I G3A=NORMAL STRESS ON PLANF AT ANGLE ALPHA -(ACTIVF) 



T3A=T2A+{ S IG3A-SI G2A )*TANA 

 C T 3 A = ^HFAR STRESS ON PLANE AT ANGLE ALPHA -(ACTIVE) 



GO TO 41 



40 $ IG3A=S IG2A 

 T3A=-T2A 



41 P3A=SQRTF( S I G3A*S I G3 A+ Ti A* T3 A ) 



C P3A = RESULT ANT STRESS ON PLANE AT ANGLE ALPHA -(ACTIVE) 



24 PRINT 25 , AL PHAD , S I G 3 A , T 3 A , P 3 A 

 2 5 FORMAT ( 1 HO , 2 4 X , F 6 . 1 , 9 X , 3 F 1 6 . 1 ) 



10 5 PRINT 66*7. * DDE LP 

 66 FORMAT (1H1»3H Z = , F 5 . 1 , 1 1 X , 1 4HD E L T A PASS I VE= , F6. 1 ) 

 PRINT 70 



70 FORMAT ( 1 H^/// , 1 9X ,4HS I Gl , 1 3 X ,2HT 1 , 1 3* , 5H C I G^P, 1 1 X *5HS I G?P , 12X , 

 13HT2P) 



PRINT 21, SIG1 ,T1 ,SIG0P,SIG2P,T2P 

 PRINT 230 



230 FORMAT ( 1 HO * 2 5X * 5HALPHA * 2 1 X * 5HS I G3 P » 1 2 X , 3HT 3 P , 1 3 X * 3HP 3P ) 

 IF (SIG0P-SIG2P) tf 1,82, 82 

 82 DO 78 I ALPHA = N » 3 60 * N 

 AL PH AD= I ALPHA 



IF ( ALPHAD-DUFLP ) 115,115*116 

 116 AL PHAD=ALPHAD+ 1 80. 



.31 



