I F ( AL PHAD- 3 6°.) 115,115,71 



71 GO TO 51 

 115 ALPHA =(ALPHAD*PI 1/180. 



TANA=S I NF( ALPHA) /COS F( ALPHA ) 



TANA2=T ANA* TANA 



IF (ALPHAD-^0. ) 76,74,76 

 76 IF ( ALPHAD-270. ) 73,74,73 



73 SIG3P=( ( 2.*SI GOP-2 .*T2P*TANA+2.*S I G2 P*TANA2 ) / ( 1 . + TANA 2 ) )-S I G2P 

 C SIG3P=NORMAL STRESS ON PLANE AT ANGLE ALPHA -(PASSIVE) 



T3P=T2P+(SIG3P-SIG2P>*TANA 

 C T 3P = SHFAR STRESS ON PLANE AT ANGLE ALPHA -(PASSIVF) 

 GO TO 75 



74 SIG3P=SIG2P 

 T3P=-T2P 



75 P3P = SQRTF( S I G3P*S I G3P+T3P*T3 P J 



C P3P=RESULTANT STRESS ON PLANE AT ANGLE ALPHA -(PASSIVE' 

 78 PRINT 25 , AL PHAD , S I G 3P , T 3 P , P 3 P 



GO TO 51 

 81 DO 84 I ALPHA = N » 3 60 * N 

 87 AL PHAD= I ALPHA 



IF ( ALPHAD+180.-DDELP) 84,85,85 



85 IF (ALP HAD- DDE LP) 86,51,51 



86 AL PHA= ( ALPHAD*P I ) / 1 80. 

 TANA=S I NF( ALPHA ) /COS F( ALPHA) 

 T A N A ? = T A N A * T A N A 



I F ( AL PHAD-90 • ) 9 1 , 99 , 9 1 

 91 IF ( ALPHAD-270. ) 93 ,99 ,93 



9 3 S I G 3 P = ( ( 2 .*S IGOP-2 ,*T2 P*TANA + 2 .*S I G2P*TANA2 ) / ( 1 -+TANA2 > )-S IG2P 



T3P=T2P+ ( SIG3P-SI G2P )*TANA 



GO TO 94 

 99 SIG3P=SIG2P 



T3P=-T2P 



94 P3P=SORTF(SlG3P- ,fC, IG3P+T3P*T3P) 



95 PRINT 25 » ALPHAD* S IG3P * T3P » P3P 

 84 CONT INUE 



51 CONTINUF 

 100 5 CALL EXIT 

 END 



32 



