Esope source files

   1 : C BET3D     SOURCE    KICH      09/09/11    21:15:00     6505           
   2 : C URGCST    SOURCE    PV        00/12/22    21:16:39     4064           
   3 :       SUBROUTINE BET3D(WRK52,WRK53,WRK54,WRK4,NSTRSS,
   4 :      1  IFOUR,DDT,IB,IGAU,HCAR)
   5 : C---------------------------------------------------------------------
   6 : C     MODELE BETON 3D
   7 : C---------------------------------------------------------------------
   8 : C
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8(A-H,O-Z)
  11 : -INC DECHE
  12 : C
  13 :       SEGMENT WRK4
  14 :         REAL*8 XE(3,NBNN)
  15 :       ENDSEGMENT
  16 : C
  17 :       SEGMENT BETINSA
  18 :        REAL*8 RT,RC,YOUN,XNU,GFT,GFC,CAR
  19 :        REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
  20 :        INTEGER IFIS,IPLA,IBB,IGAU1
  21 :       ENDSEGMENT
  22 : C      
  23 :       NCOMAT = nmatt
  24 :       NBNN=XE(/2)
  25 :       NVARI=VAR0(/1)
  26 :       SEGINI BETINSA
  27 : C
  28 :       CAR = HCAR
  29 : C
  30 :       YOUN = XMAT( 1)
  31 :       XNU  = XMAT( 2)
  32 :       RHO  = XMAT( 3)
  33 :       ALPH = XMAT( 4)
  34 :       RT   = XMAT( 5)
  35 :       RC   = XMAT( 6)
  36 :       GFT  = XMAT( 7)
  37 :       GFC  = XMAT( 8)
  38 : C
  39 : C---------------------------------------------------------------------
  40 : C
  41 : C     CORRESPONDANCE DES VARIABLES POUR URGCST
  42 : C
  43 : C
  44 :       IBB = IB
  45 :       IGAU1 = IGAU
  46 : C
  47 : C------------------------------------------------------
  48 : C   INITIALISATION DES VARIABLES INTERNES 
  49 : C------------------------------------------------------
  50 : C       A: MODELE PLASTIQUE
  51 : C
  52 :         IFIS  = INT(SNGL(VAR0(1)))
  53 :         IPLA   = INT(SNGL(VAR0(2)))
  54 :         SEQT   = VAR0(3)
  55 :         SEQC   = VAR0(4)
  56 :         DKT    = VAR0(5)
  57 :         DKC    = VAR0(6)
  58 :         ENDT    = VAR0(7)
  59 :         ENDC    = VAR0(8)
  60 :         IF (VAR0(9).EQ.0.D0) THEN
  61 : *******    COERAND = RAND()
  61 :            CALL TDRAND(COERAND)
  62 :          ELSE 
  63 :            COERAND = VAR0(9)
  64 :         ENDIF
  65 :         RT = (1. + ((2.*COERAND - 1.) * 0.1))*RT
  66 : C
  67 :         IF (SEQT.EQ.0.D0) THEN 
  68 :         SEQT = RT
  69 :         ENDIF
  70 :         IF (SEQC.EQ.0.D0) THEN 
  71 :         SEQC = RC/3.D0
  72 :         ENDIF
  73 : C
  74 : C---------------------------------------------------------------------
  75 : C
  76 :       GO TO (10,10,10,10,10,30),NSTRSS
  77 :    10 CONTINUE
  78 :       KERRE=437
  79 :       WRITE(*,*) '!! ATTENTION DANS BET3D  NSTRSS=',NSTRSS
  80 :       STOP
  81 :       GO TO 1000
  82 : C
  83 :    40 CONTINUE
  84 : C
  85 :    30 CONTINUE
  86 : C
  87 : C---------------------------------------------------------------------
  88 : C
  88 :          CALL BONE3D(SIG0,SIGF,DEPST,NSTRSS,BETINSA)
  89 : C
  90 : C---------------------------------------------------------------------
  91 : C              ACTUALISATION DES VARIABLES INTERNES
  92 : C---------------------------------------------------------------------
  93 : C       A: MODELE PLASTIQUE
  94 : C
  95 :       VARF( 1)=FLOAT(IFIS)
  96 :       VARF( 2)=FLOAT(IPLA)
  97 :       VARF( 3)=SEQT
  98 :       VARF( 4)=SEQC
  99 :       VARF( 5)=DKT     
 100 :       VARF( 6)=DKC
 101 :       VARF( 7)=ENDT     
 102 :       VARF( 8)=ENDC
 103 :       VARF( 9)=COERAND
 104 : 
 105 : C
 106 : C---------------------------------------------------------------------
 107 : C
 108 :  1000 CONTINUE
 109 : C
 110 :       SEGSUP BETINSA
 111 : C
 112 :       RETURN
 113 :       END
 114 : 
 115 : 
 116 :  
 117 :  
 118 :  
 119 :  
 120 :  
 121 :  
 122 :  
 123 :