Esope source files

   1 : C FORECR    SOURCE    CHAT      05/01/13    00:06:35     5004           
   2 :       SUBROUTINE FORECR(DK,PAEC,IC,SEQ,BETINSA)
   3 : C
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 :       SEGMENT BETINSA
   7 :        REAL*8 RT,RC,YOUN,XNU,GFT,GFC,HC
   8 :        REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
   9 :        INTEGER IFIS,IPLA,IBB,IGAU
  10 :       ENDSEGMENT
  11 : C 
  12 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  13 : *                                                                    *
  14 : *    CALCUL DU PARAMETRE D'ECROUISSAGE                               *
  15 : *    _________________________________                               *
  16 : *                                                                    *
  17 : *    ENTREES:                                                        *
  18 : *                                                                    *
  19 : *        DK        VARIABLE D'ECROUISSAGE                            *
  20 : *        IC        =1 : COURBE de TRACTION                           *
  21 : *                  =2 : COURBE de COMPRESSION                        * 
  22 : *                                                                    *
  23 : *    SORTIES:                                                        *
  24 : *                                                                    *
  25 : *        SEQ       CONTRAINTE EQUIVALENTE                            *
  26 : *        PAEC      PENTE                                             *
  27 : *                                                                    *
  28 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  29 : C
  30 :         BB=3.D0
  31 :         DT = 0.5
  32 :         DC = 0.25
  33 :         ENDL = 0.99      
  34 : C   
  35 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  36 : C                      Direction de traction   
  37 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  38 : C
  39 :            IF (IC.EQ.1) THEN 
  40 :        F=RT
  41 :             A= -0.5D0   
  42 :             B=F*(HC/GFT)*(1.D0+A/2.D0)
  43 :             C=B*LOG(1.D0-DT)/(LOG(((1+A)-SQRT(1+A*A))/(2.D0*A)))
  44 : C
  45 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  46 : C                      Direction de compression  
  47 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  48 : C
  49 :            ELSE
  50 :        F=RC/BB
  51 :             A=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
  52 :             B=F*(HC/GFC)*(1.D0+A/2.D0)
  53 :        C=B*LOG(1.D0-DC)/(LOG((1+A)/(2.D0*A)))
  54 :            ENDIF
  55 : C
  56 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    
  57 : C       CALCUL LA CONTRAINTE UNIAXIALE ET DE LA PENTE
  58 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  59 : C       1: CONTRAINTE UNIAXIALE
  60 : C
  61 :             EXP1=EXP(-B*DK)
  62 :             DD=1.D0-C/B
  63 :             SEQ=F*((1.D0+A)*EXP1**(DD)-A*EXP1**(1.D0+DD))
  64 : C
  65 : C       2: CALCUL DE LA (+)*PENTE 
  66 : C
  67 :             PAEC=-B*F*(DD*(1.D0+A)*EXP1**(DD)-
  68 :      &                 (1.D0+DD)*A*EXP1**(1.D0+DD))
  69 : C
  70 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  71 : C
  72 :            IF(ABS(SEQ).LT.1.D-8)THEN
  73 :             SEQ=1.D-8
  74 :             PAEC=0.D0
  75 :             GOTO 10
  76 :            ENDIF        
  77 : C
  78 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  79 : C
  80 :    10   CONTINUE
  81 : C
  82 :         RETURN
  83 :         END 
  84 :  
  85 :