Esope source files

   1 : C DAMAG1    SOURCE    CHAT      05/01/12    22:36:09     5004           
   2 :       SUBROUTINE DAMAG1(DKPT,DKPC,ENDC,ENDT,IFISU,IPLA,IC
   3 :      &  ,BETJEF,NECH0,NECH1)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : C
   8 : *      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,      
   9 : *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  10 : *     &                ITER,ISIM,IBB,IGAU1,IZON
  11 : *      COMMON /DBETNECH0/DT,DC,ALFG,S0
  12 : *      COMMON /DBETNECH1/ENDL
  13 :       SEGMENT BETJEF
  14 :        REAL*8 AA,BETA,RB,ALFA,EX,XNU,GFC,GFT,HC,ETA,TDEF,
  15 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  16 :        INTEGER ICT,ICC,IMOD,IVIS,ITER,
  17 :      &                ISIM,IBB,IGAU,IZON
  18 :       ENDSEGMENT
  19 :       SEGMENT NECH0
  20 :        REAL*8 DT,DC,ALFG,S0,ENDO
  21 :       ENDSEGMENT
  22 :       SEGMENT NECH1
  23 :        REAL*8 ENDL
  24 :       ENDSEGMENT
  25 : C
  26 : *--------------------------------------------------------------------*
  27 : *                                                                    *
  28 : *    CALCUL DE L'ENDOMMAGEMENT                                       *
  29 : *    _________________________                                       *
  30 : *                                                                    *
  31 : *    ENTREES:                                                        *
  32 : *                                                                    *
  33 : *        DKT       PARAMETRE D'ECROUISSAGE EN TRACTION               *
  34 : *        DKC       PARAMETRE D'ECROUISSAGE EN COMPRESSION            *
  35 : *        IC        =1 : COURBE de TRACTION                           *
  36 : *                  =2 : COURBE de COMPRESSION                        *
  37 : *                       =3 : COURBE TRACTION COMPRESSION                  *  
  38 : *                                                                    *
  39 : *    SORTIES:                                                        *
  40 : *                                                                    *
  41 : *        ENDC      ENDOMMAGEMENT DE COMPRESSION                      *
  42 : *        ENDT      ENDOMMAGEMENT DE TRACTION                         *
  43 : *                                                                    *
  44 : *--------------------------------------------------------------------*
  45 : C
  46 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  47 : C                      Initialisation   
  48 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  49 : C
  50 :       BB=1.D0/AA
  51 : C
  52 : C      
  53 : C      
  54 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  55 : C                      Direction de traction   
  56 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  57 : C
  58 :            IF (IC.EQ.1.OR.IC.EQ.3) THEN
  59 :          FT=ALFA*RB
  60 :               AT= -0.5D0        
  61 :               BT=FT*(HC/GFT)*(1.D0+AT/2.D0)
  62 :              CT=BT*LOG(1.D0-DT)/(LOG(((1+AT)-SQRT(1+AT*AT))/(2.D0*AT)))
  63 : C        
  64 :         ENDT=1- EXP(-CT*DKPT)
  65 :            ENDIF
  66 : C      
  67 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  68 : C                      Direction de compression   
  69 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  70 : C
  71 :         IF (IC.EQ.2.OR.IC.EQ.3) THEN
  72 :          FC=AA*RB
  73 :               AC=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)     
  74 :               BC=FC*(HC/GFC)*(1.D0+AC/2.D0) 
  75 :               CC=BC*LOG(1.D0-DC)/(LOG((1+AC)/(2.D0*AC)))
  76 : C
  77 :          ENDC=1- EXP(-CC*DKPC)
  78 :         ENDIF
  79 : C
  80 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  81 : C                    CAS OU ENDO=ENDO_LIMITE
  82 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  83 : C
  84 :         IF(ENDT.GE.ENDL) ENDT= ENDL
  85 :         IF(ENDC.GE.ENDL) ENDC= ENDL     
  86 : C       
  87 :             IF ((ENDT.LT.0.D0).OR.(ENDC.LT.0.D0)) THEN
  88 :               WRITE(*,*)'PROBLEME DANS DAMAGE'
  89 :               WRITE(*,*)'ENDT=',ENDT,'ENDC=',ENDC
  90 : C             STOP
  91 :             ENDIF
  92 : C
  93 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  94 : C                       
  95 :         RETURN
  96 :         END 
  97 :  
  98 :