Sources Esope | Cast3M

princ
C PRINC     SOURCE    CB215821  16/04/21    21:18:06     8920
C PRINC     SOURCE    INSL       24/10/96
      SUBROUTINE PRINC(S,T,NSTRS)
C
C   CALCUL DES CONTRAINTES ET DIRECTIONS PRINCIPALES
C    L'ANGLE PHI=T(4) : EST EN DEGRE
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION S(4),T(4)
      PI=4.D0*ATAN(1.D0)
      A1=S(1)+S(2)
      A2=S(1)-S(2)
      A3=S(3)
      A4=S(3)
      IF(NSTRS.EQ.4.OR.NSTRS.EQ.6) A4=S(4)
      IF(ABS(A4) .GE. 1.D-8) GOTO 11
      IF(A2 .LT. 1.D-8) GOTO 121
      T(1)=S(1)
      T(2)=S(2)
      T(3)=S(3)
      T(4)=0.0D+00
      GO TO 10
  121 T(1)=S(2)
      T(2)=S(1)
      T(3)=S(3)
      T(4)=90.0D+00
      GO TO 10
   11 IF(ABS(A2) .GE. 1.D-8) GOTO 111
      T(1)=S(1)+ABS(A4)
      T(2)=S(2)-ABS(A4)
      T(3)=S(3)
      T(4)=45.0D+00*ABS(A4)/A4
      GO TO 10
  111 A5=SQRT(A2*A2+4.0D+00*A4*A4)
      T(1)=(A1+A5)/2.0D+00
      T(2)=(A1-A5)/2.0D+00
      T(3)=S(3)
      T(4)=ATAN(2.0D+00*A4/A2)*90.0D+00/PI
      IF(A2 .GE. 1.D-8) GOTO 10
      T(4)=T(4)+90.0D+00*ABS(A4)/A4
  10  CONTINUE
      IF(NSTRS.EQ.3.OR.NSTRS.EQ.5) T(3)=0.D0
      RETURN
      END