Esope source files

   1 : C REELLE    SOURCE    PV        11/03/07    21:18:04     6885           
   2 : C     REELLE    SOURCE   NECHNECH   99/06/08
   3 :       SUBROUTINE REELLE(S1,DEP,DEPSITF1,DEPSITF2,DEPE,S2,ENDM,
   4 :      &            NSTRS,IFIS,IPLA,BETJEF,NECH0,NECH1)
   5 : C
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 :       DIMENSION S1(4), S2(4), V1(4)
   9 :       DIMENSION DEP(4,4),DEPE(4,4)
  10 : *      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,      
  11 : *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  12 : *     &                ITER,ISIM,IBB,IGAU1,IZON      
  13 : *      COMMON /DBETNECH0/DT,DC,ALFG,S0
  14 : *      COMMON /DBETNECH1/ENDL
  15 :       SEGMENT BETJEF
  16 :        REAL*8 AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  17 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  18 :        INTEGER ICT,ICC,IMOD,IVIS,ITER,
  19 :      &                ISIM,IBB,IGAU,IZON
  20 :       ENDSEGMENT
  21 :       SEGMENT NECH0
  22 :        REAL*8 DT,DC,ALFG,S0,ENDO
  23 :       ENDSEGMENT
  24 :       SEGMENT NECH1
  25 :        REAL*8 ENDL
  26 :       ENDSEGMENT
  27 : C
  28 : *--------------------------------------------------------------------*
  29 : *    CALCUL  LA MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE ET            *
  30 : *    DE LA CONTRAINTE REELLE                                         *
  31 : *                                                                    *
  32 : *                                                                    *
  33 : *    ENTREES:                                                        *
  34 : *                                                                    *
  35 : *        S1        CONTRAINTE EFFECTIVE                              *
  36 : *        NSTRS     NOMBRE DE CONTRAINTES                             *
  37 : *        DPSTF1    DEFORMATION EQUIVALENTE EN TRACTION               *
  38 : *        DPSTF2    DEFORMATION EQUIVALENTE EN COMPRESSION            *
  39 : *        DEP       MATRICE ELSATO-PLASTIQUE                          *
  40 : *                                                                    *
  41 : *    SORTIES:                                                        *
  42 : *                                                                    *
  43 : *        S2        CONTRAINTE REELLE                                 *
  44 : *        DEPE      MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE            *
  45 : *        ENDO      ENDOMMAGEMENT                                     *
  46 : *--------------------------------------------------------------------*
  47 : C       INITIALISATION
  48 : C
  49 : C      S0 = 0.25
  50 :       ENDT = 0.D0
  51 :       ENDC = 0.D0
  52 :       ENDM = 0.D0
  53 :            NN=NSTRS-1
  53 :         CALL ZERO(DEPE,4,4)
  53 :         CALL ZERO(S2,4,1)
  54 :         SOM1 = 0.D0
  55 :         SOM2 = 0.D0
  56 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  57 : C       CALCUL DE L'ENDOMMAGEMENT
  58 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  59 : C       WRITE(*,*)'K1=',DEPSITF1
  60 : C       WRITE(*,*)'K2=',DEPSITF2
  61 : C
  61 :            CALL DAMAG1(DEPSITF1,DEPSITF2,ENDC,
  62 :      &          ENDT,IFIS,IPLA,3,BETJEF,NECH0,NECH1)
  63 : C                        
  64 : C
  65 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  66 : C       CALCUL DES VALEURS PROPRES DE LA MATRICE DE VECTEUR CONTRAINTE
  67 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  68 : C
  68 :            CALL PRINC(S1,V1,NSTRS)
  69 : C
  70 :            DO 10 I=1, NN
  71 :             SOM1=SOM1+(V1(I)+DABS(V1(I)))/2.D0
  72 :             SOM2=SOM2+DABS(V1(I))
  73 :    10 CONTINUE
  74 :         IF(SOM2.EQ.0.D0) THEN
  75 :               SOM2=SOM2+1.D-10
  76 :              ENDIF
  77 :            RSIG=SOM1/SOM2
  78 :            RSIG=ABS(RSIG)
  79 :            S=S0+(1-S0)*RSIG
  80 : C      S = 1.D0
  81 : C
  82 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  83 : C       CALCUL DE L'ENDOMMAGEMENT TOTAL
  84 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  85 : C
  86 :            ENDM=1.D0-((1.D0-ENDC)*(1.D0-S*ENDT))
  87 : C
  88 :            IF((ENDM.LT.0.D0).OR.(ENDM.GT.1D0))THEN
  89 :              WRITE(*,*)'ENDO=',ENDM
  90 : C            STOP
  91 :            ENDIF
  92 : C
  93 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  94 : C
  95 :         DO 20 I=1, NSTRS
  96 :            S2(I)=(1.D0-ENDM)*S1(I)
  97 :    20 CONTINUE
  98 : C       
  99 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
 100 : C       CALCUL DE LA MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE
 101 : C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 102 : C
 103 :       DO 30 I=1, NSTRS
 104 :       DO 40 J=1, NSTRS
 105 :        DEPE(I,J)=(1.D0-ENDM)*DEP(I,J)
 106 :    40 CONTINUE
 107 :    30 CONTINUE
 108 : C
 109 :       RETURN
 110 :       END 
 111 :  
 112 :  
 113 :