Esope source files

   1 : C ENDAME    SOURCE    CHAT      05/01/12    23:38:48     5004           
   2 :       SUBROUTINE ENDAME(IC,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 : *--------------------------------------------------------------------*
  13 : *                                                                    *
  14 : *    CALCUL DE L'ENDOMMAGEMENT                                       *
  15 : *    _________________________                                       *
  16 : *                                                                    *
  17 : *    ENTREES:                                                        *
  18 : *                                                                    *
  19 : *        DKT       PARAMETRE D'ECROUISSAGE EN TRACTION               *
  20 : *        DKC       PARAMETRE D'ECROUISSAGE EN COMPRESSION            *
  21 : *        IC        =1 : COURBE de TRACTION                           *
  22 : *                  =2 : COURBE de COMPRESSION                        *  
  23 : *                                                                    *
  24 : *    SORTIES:                                                        *
  25 : *                                                                    *
  26 : *        ENDC      ENDOMMAGEMENT DE COMPRESSION                      *
  27 : *        ENDT      ENDOMMAGEMENT DE TRACTION                         *
  28 : *                                                                    *
  29 : *--------------------------------------------------------------------*
  30 : C
  31 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  32 : C                      Initialisation   
  33 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  34 : C
  35 :         BB=3.D0
  36 :         DT = 0.5D0
  37 :         DC = 0.25D0
  38 :         ENDL = 0.99D0
  39 : C
  40 : C      
  41 : C      
  42 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  43 : C                      Direction de traction   
  44 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  45 : C
  46 :            IF (IC.EQ.1) THEN
  47 :         FT=RT
  48 :              AT= -0.5D0 
  49 :              BT=FT*(HC/GFT)*(1.D0+AT/2.D0)
  50 :              CT=BT*LOG(1.D0-DT)/(LOG(((1+AT)-SQRT(1+AT*AT))/(2.D0*AT)))
  51 : C        
  52 :         ENDT=1- EXP(-CT*DKT)
  53 : C
  54 :         IF (ENDT.GE.ENDL) THEN
  55 :          IFIS=2
  56 :         ENDIF
  57 :         IF (ENDT.GT.0.D0) THEN
  58 :          IFIS=1
  59 :         ENDIF
  60 :            ENDIF
  61 : C      
  62 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  63 : C                      Direction de compression   
  64 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  65 : C
  66 :         IF (IC.EQ.2) THEN
  67 :          FC=RC
  68 :         AC=(2.D0*BB-1.D0)+2.D0* SQRT(BB**2-BB)
  69 :         BC=FC*(HC/GFC)*(1.D0+AC/2.D0) 
  70 :         CC=BC*LOG(1.D0-DC)/(LOG((1+AC)/(2.D0*AC)))
  71 : C
  72 :          ENDC=1- EXP(-CC*DKC)
  73 : C
  74 :          IF (ENDC.GE.ENDL) THEN
  75 :           IPLA=3
  76 :          ENDIF
  77 :          IF (ENDC.GE.DC) THEN
  78 :           IPLA=2
  79 :          ENDIF
  80 :          IF (ENDC.GT.0.D0) THEN
  81 :           IPLA=1
  82 :          ENDIF
  83 :         ENDIF
  84 : C
  85 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  86 : C                    CAS OU ENDO=ENDO_LIMITE
  87 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  88 : C
  89 :         IF(ENDT.GE.ENDL) ENDT= ENDL
  90 :         IF(ENDC.GE.ENDL) ENDC= ENDL     
  91 : C       
  92 :             IF ((ENDT.LT.0.D0).OR.(ENDC.LT.0.D0)) THEN
  93 :               WRITE(*,*)'PROBLEME DANS DAMAGE'
  94 :               WRITE(*,*)'ENDT=',ENDT,'ENDC=',ENDC
  95 : C             STOP
  96 :             ENDIF
  97 : C
  98 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  99 : C                       
 100 :         RETURN
 101 :         END 
 102 :  
 103 :