Esope source files

   1 : C OTTEST    SOURCE    CHAT      05/01/13    02:06:57     5004           
   2 :       SUBROUTINE OTTEST(XINVL,W,WMAX,SMAX,SIGMA,XLTR,PRECIZ,
   3 :      &                  BTR,WRUPT,
   4 :      &                  IB,IGAU,ICAS,KERRE)
   5 : *
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)     
   8 : -INC CCOPTIO
   9 : C 
  10 : C
  11 :       DIMENSION SIGMA(6)
  12 : C
  13 :       DIMENSION WMAX(3),SMAX(7),W(3),
  14 :      &    XLTR(3),XINVL(3),WRUPT(3)
  15 : C
  16 :       DIMENSION WREOUV(3)
  17 : 
  18 : C
  19 : C     ICAS = 1   EN ENTREE
  20 : C     ICAS = 2   EN SORTIE
  21 : C
  22 :       XZER=0.D0
  23 :       KERRE=0
  24 : C
  25 : C------------------------------------------------------------------------
  26 : C            ON PASSE EN REVUE LES 3 DIRECTIONS POTENTIELLES
  27 : C                          DE FISSURATION
  28 : C------------------------------------------------------------------------
  29 : C
  30 :       DO 88 I=1,3
  31 : 
  32 :          IF(XINVL(I).NE.0.D0) THEN
  33 : 
  34 :            WREOUV(I)=BTR*MIN(WMAX(I),WRUPT(I))
  35 : *
  36 : *          TEST SUR OUVERTURE MINIMALE
  37 : *
  38 :            IF(W(I).LT.WREOUV(I)-1.D-10*WRUPT(I))  THEN
  39 :              KERRE=70
  40 :              WRITE(IOIMP,60000) ICAS
  41 : 60000        FORMAT (/2X, '         PROBLEME   -  ICAS =',I3 //)
  42 :              WRITE(IOIMP,67675) I,W(I), WREOUV(I),WMAX(I)
  43 : 67675        FORMAT(2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$' /
  44 :      &        4X,'ATTENTION   W < WREOUV '/
  45 :      &        2X, ' I=',I4,2X,'W=',1PE12.5,2X,'WREOUV=',1PE12.5,
  46 :      &        2X,'WMAX=',1PE12.5/
  47 :      &        2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'/)
  48 :           ENDIF
  49 : *
  50 : *          TEST SUR DEPASSEMENT DE LA LIMITE EN TRACTION
  51 : *
  52 :            IF(SIGMA(I)-XLTR(I).GT.PRECIZ) THEN
  53 :              KERRE=70
  54 :              WRITE(IOIMP,60000) ICAS
  55 :              WRITE(IOIMP,77058) I,SIGMA(I)
  56 : 77058        FORMAT(2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$' /
  57 :      &        4X,'ATTENTION   DEPASSEMENT DE RT ' /
  58 :      &        2X, ' I=',I4,2X,'SIGMA=',1PE12.5 /
  59 :      &        2X,'$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'/)
  60 :            ENDIF
  61 : 
  62 : *
  63 : *          TEST SUR SECANTE
  64 : *
  65 :           IF(BTR.NE.1.D0) THEN
  66 :             IF(WMAX(I).LT.WRUPT(I).AND.W(I).GT.WREOUV(I)) THEN
  67 :               PPP = SMAX(I)/( WMAX(I)-WREOUV(I))
  68 :               ZZOB = W(I) -WREOUV(I)
  69 :               SSS = ZZOB*PPP
  70 :               IF(SIGMA(I).GT.PRECIZ.AND.
  71 :      &           ABS(SIGMA(I)-SSS).GT.PRECIZ) THEN
  72 : 
  73 :              KERRE=70
  74 :              WRITE(IOIMP,60000) ICAS
  75 :                 WRITE(IOIMP,79912) I,SIGMA(I),SSS,WMAX(I),
  76 :      &                       W(I),WREOUV(I),ZZOB,SMAX(I)
  77 : 79912     FORMAT(2X,' ################################### '/
  78 :      &    2X,'I=',I4,2X,'SIGMA=',1PE12.5,
  79 :      &    2X,'SSS=',1PE12.5/2X,'WMAX=',1PE12.5/
  80 :      &    2X,'W=',1PE12.5,2X,'WREOUV=',1PE12.5,
  81 :      &    2X,'W-WREOUV = ',1PE12.5,2X,'SMAX= ',1PE12.5/
  82 :      &        2X,' ################################### '/)
  83 :            PRINT *,'SIGMA=',SIGMA(I)
  84 :            PRINT *,'SSS  =',SSS
  85 :            PRINT *,'WMAX =',WMAX(I)
  86 :            PRINT *,'W    =',W(I)
  87 :            PRINT *,'WREOUV=',WREOUV(I)
  88 :            PRINT *,'W - WREOUV=',ZZOB
  89 :            PRINT *,'SMAX    =',SMAX(I)
  90 : 
  91 :           ENDIF
  92 :         ENDIF
  93 :         ENDIF
  94 : 
  95 :       ENDIF
  96 : 88    CONTINUE
  97 : 
  98 :       IF(KERRE.NE.0) THEN
  99 :           WRITE(IOIMP,70902) IB,IGAU
 100 : 70902     FORMAT(2X,'SORTIE OTTEST  IB=',I6,2X,'IGAU=',I3////)
 101 :       ENDIF
 102 : C
 103 :       RETURN
 104 :       END     
 105 :  
 106 :  
 107 :  
 108 :  
 109 :  
 110 :  
 111 :  
 112 :  
 113 :  
 114 :  
 115 :  
 116 :  
 117 :