Esope source files

   1 : C UNICO1    SOURCE    CHAT      05/01/13    03:58:30     5004           
   2 : C     UNICO1    SOURCE   NECH   99/05/26
   3 :       SUBROUTINE UNICO1(DK,PAEC,IC,SEQ,BETJEF,NECH0,NECH1)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : *      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,      
   8 : *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
   9 : *     &                ITER,ISIM,IBB,IGAU1,IZON
  10 : *      COMMON /DBETNECH0/DT,DC,ALFG,S0
  11 : *      COMMON /DBETNECH1/ENDL
  12 :       SEGMENT BETJEF
  13 :        REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,TDEF,
  14 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  15 :        INTEGER ICT,ICC,IMOD,IVISS,ITER,
  16 :      &                ISIM,IBB,IGAU1,IZON
  17 :       ENDSEGMENT
  18 :       SEGMENT NECH0
  19 :        REAL*8 DT,DC,ALFG,S0,ENDO
  20 :       ENDSEGMENT
  21 :       SEGMENT NECH1
  22 :        REAL*8 ENDL
  23 :       ENDSEGMENT  
  24 : C 
  25 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  26 : *                                                                    *
  27 : *    CALCUL DU PARAMETRE D'ECROUISSAGE                               *
  28 : *    _________________________________                               *
  29 : *                                                                    *
  30 : *    ENTREES:                                                        *
  31 : *                                                                    *
  32 : *        DK        VARIABLE D'ECROUISSAGE                            *
  33 : *        IC        =1 : COURBE de TRACTION                           *
  34 : *                  =2 : COURBE de COMPRESSION                        * 
  35 : *                                                                    *
  36 : *    SORTIES:                                                        *
  37 : *                                                                    *
  38 : *        SEQ       CONTRAINTE EQUIVALENTE                            *
  39 : *        PAEC      PENTE                                             *
  40 : *                                                                    *
  41 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  42 : C
  43 :         BB=1.D0/AA      
  44 : C   
  45 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  46 : C                      Direction de traction   
  47 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  48 : C
  49 :            IF (IC.EQ.1) THEN 
  50 :        F=ALFA*RB
  51 :             A= -0.5D0   
  52 :             B=F*(HC/GFT)*(1.D0+A/2.D0)
  53 :             C=B*LOG(1.D0-DT)/(LOG(((1+A)-SQRT(1+A*A))/(2.D0*A)))
  54 : C
  55 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  56 : C                      Direction de compression  
  57 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  58 : C
  59 :            ELSE
  60 :        F=AA*RB
  61 :             A=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
  62 :             B=F*(HC/GFC)*(1.D0+A/2.D0)
  63 :             C=B*LOG(1.D0-DC)/(LOG((1+A)/(2.D0*A)))
  64 :            ENDIF
  65 : C
  66 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  67 : C       CALCUL DE LA DEFORMATION PLASTIQUE LIMITE
  68 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  69 : C
  70 :         DKL=LOG(1.D0-ENDL)/(-C)
  71 : C
  72 :            IF(DK.GE.DKL)THEN
  73 :             DK=DKL
  74 :            ENDIF        
  75 : C
  76 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~    
  77 : C       CALCUL LA CONTRAINTE UNIAXIALE ET DE LA PENTE
  78 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  79 : C       1: CONTRAINTE UNIAXIALE
  80 : C
  81 :             EXP1=EXP(-B*DK)
  82 :             DD=1.D0-C/B
  83 :             SEQ=F*((1.D0+A)*EXP1**(DD)-A*EXP1**(1.D0+DD))
  84 : C
  85 : C       2: CALCUL DE LA (+)*PENTE 
  86 : C
  87 :             PAEC=-B*F*(DD*(1.D0+A)*EXP1**(DD)-
  88 :      &                  (1.D0+DD)*A*EXP1**(1.D0+DD))
  89 : C
  90 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  91 : C
  92 :            IF(ABS(SEQ).LT.1.D-8)THEN
  93 :             SEQ=1.D-8
  94 :             PAEC=0.D0
  95 :             GOTO 10
  96 :            ENDIF        
  97 : C
  98 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  99 : C
 100 :    10   CONTINUE
 101 : C
 102 :         RETURN
 103 :         END 
 104 :  
 105 :