Esope source files

   1 : C TRACTI    SOURCE    CHAT      05/01/13    03:44:22     5004           
   2 :       SUBROUTINE TRACTI (RESU,EPST,SIG,EPS,NCOURB,ICLE,IBI)                     
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       DIMENSION SIG(*),EPS(*)                                                   
   6 : C                                                                               
   7 : C  ICLE = 1  ON CALCULE LA PENTE A LA COURBE DE TRACTION                        
   8 : C  ICLE = 2  ON CALCULE LA CONTRAINTE                                           
   9 : C  IBI NON NUL  : ERREUR : ON DEPASSE LE BOUT DE LA COURBE
  10 : 
  11 : C                                                                               
  12 : C                ON INTERPOLE LINEAIREMENT                                      
  13 : C                                                                               
  14 :       IBI=0
  15 : *
  16 : *    TEST POUR VOIR SI ON DEPASSE EN EPSILON
  17 : *
  18 :       IF(EPST.GT.1.1D0*EPS(NCOURB)) THEN
  19 :          IBI=1
  20 :          RETURN
  21 :       ENDIF
  22 : *
  23 :       DO 11 I=2,NCOURB                                                          
  24 :       IP=I                                                                      
  25 :       IF(EPS(I).GT.EPST) GO TO 13                                               
  26 :   11  CONTINUE                                                                  
  27 :   13  CONTINUE                                                                  
  28 :       IP1=IP-1                                                                  
  29 :       PENTE=(SIG(IP)-SIG(IP1))/(EPS(IP)-EPS(IP1))                               
  30 :       GO TO (1,2),ICLE                                                          
  31 : 1     RESU=PENTE                                                                
  32 :       RETURN                                                                    
  33 : 2     RESU=SIG(IP1)+PENTE*(EPST-EPS(IP1))                                       
  34 :       RETURN                                                                    
  35 :       END                                                                       
  36 :  
  37 :