Sources Esope | Cast3M

ottest
C OTTEST    SOURCE    CHAT      05/01/13    02:06:57     5004
      SUBROUTINE OTTEST(XINVL,W,WMAX,SMAX,SIGMA,XLTR,PRECIZ,
     &                  BTR,WRUPT,
     &                  IB,IGAU,ICAS,KERRE)
*
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC PPARAM
-INC CCOPTIO
C
C
      DIMENSION SIGMA(6)
C
      DIMENSION WMAX(3),SMAX(7),W(3),
     &    XLTR(3),XINVL(3),WRUPT(3)
C
      DIMENSION WREOUV(3)
 
C
C     ICAS = 1   EN ENTREE
C     ICAS = 2   EN SORTIE
C
      XZER=0.D0
      KERRE=0
C
C------------------------------------------------------------------------
C            ON PASSE EN REVUE LES 3 DIRECTIONS POTENTIELLES
C                          DE FISSURATION
C------------------------------------------------------------------------
C
      DO 88 I=1,3
 
         IF(XINVL(I).NE.0.D0) THEN
 
           WREOUV(I)=BTR*MIN(WMAX(I),WRUPT(I))
*
*          TEST SUR OUVERTURE MINIMALE
*
           IF(W(I).LT.WREOUV(I)-1.D-10*WRUPT(I))  THEN
             KERRE=70
             WRITE(IOIMP,60000) ICAS
60000        FORMAT (/2X, '         PROBLEME   -  ICAS =',I3 //)
             WRITE(IOIMP,67675) I,W(I), WREOUV(I),WMAX(I)
67675        FORMAT(2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$' /
     &        4X,'ATTENTION   W &lt; WREOUV '/
     &        2X, ' I=',I4,2X,'W=',1PE12.5,2X,'WREOUV=',1PE12.5,
     &        2X,'WMAX=',1PE12.5/
     &        2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'/)
          ENDIF
*
*          TEST SUR DEPASSEMENT DE LA LIMITE EN TRACTION
*
           IF(SIGMA(I)-XLTR(I).GT.PRECIZ) THEN
             KERRE=70
             WRITE(IOIMP,60000) ICAS
             WRITE(IOIMP,77058) I,SIGMA(I)
77058        FORMAT(2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$' /
     &        4X,'ATTENTION   DEPASSEMENT DE RT ' /
     &        2X, ' I=',I4,2X,'SIGMA=',1PE12.5 /
     &        2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'/)
           ENDIF
 
*
*          TEST SUR SECANTE
*
          IF(BTR.NE.1.D0) THEN
            IF(WMAX(I).LT.WRUPT(I).AND.W(I).GT.WREOUV(I)) THEN
              PPP = SMAX(I)/( WMAX(I)-WREOUV(I))
              ZZOB = W(I) -WREOUV(I)
              SSS = ZZOB*PPP
              IF(SIGMA(I).GT.PRECIZ.AND.
     &           ABS(SIGMA(I)-SSS).GT.PRECIZ) THEN
 
             KERRE=70
             WRITE(IOIMP,60000) ICAS
                WRITE(IOIMP,79912) I,SIGMA(I),SSS,WMAX(I),
     &                       W(I),WREOUV(I),ZZOB,SMAX(I)
79912     FORMAT(2X,' ################################### '/
     &    2X,'I=',I4,2X,'SIGMA=',1PE12.5,
     &    2X,'SSS=',1PE12.5/2X,'WMAX=',1PE12.5/
     &    2X,'W=',1PE12.5,2X,'WREOUV=',1PE12.5,
     &    2X,'W-WREOUV = ',1PE12.5,2X,'SMAX= ',1PE12.5/
     &        2X,' ################################### '/)
           PRINT *,'SIGMA=',SIGMA(I)
           PRINT *,'SSS  =',SSS
           PRINT *,'WMAX =',WMAX(I)
           PRINT *,'W    =',W(I)
           PRINT *,'WREOUV=',WREOUV(I)
           PRINT *,'W - WREOUV=',ZZOB
           PRINT *,'SMAX    =',SMAX(I)
 
          ENDIF
        ENDIF
        ENDIF
 
      ENDIF
88    CONTINUE
 
      IF(KERRE.NE.0) THEN
          WRITE(IOIMP,70902) IB,IGAU
70902     FORMAT(2X,'SORTIE OTTEST  IB=',I6,2X,'IGAU=',I3////)
      ENDIF
C
      RETURN
      END