bet3d
C BET3D SOURCE PV 17/12/08 21:15:10 9660 C URGCST SOURCE PV 00/12/22 21:16:39 4064 1 IFOUR,DDT,IB,IGAU,HCAR) C--------------------------------------------------------------------- C MODELE BETON 3D C--------------------------------------------------------------------- C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC DECHE C SEGMENT WRK4 REAL*8 XE(3,NBNN) ENDSEGMENT C SEGMENT BETINSA REAL*8 RT,RC,YOUN,XNU,GFT,GFC,CAR REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC INTEGER IFIS,IPLA,IBB,IGAU1 ENDSEGMENT C NCOMAT = nmatt NBNN=XE(/2) NVARI=VAR0(/1) SEGINI BETINSA C CAR = HCAR C YOUN = XMAT( 1) XNU = XMAT( 2) RHO = XMAT( 3) ALPH = XMAT( 4) RT = XMAT( 5) RC = XMAT( 6) GFT = XMAT( 7) GFC = XMAT( 8) C C--------------------------------------------------------------------- C C CORRESPONDANCE DES VARIABLES POUR URGCST C C IBB = IB IGAU1 = IGAU C C------------------------------------------------------ C INITIALISATION DES VARIABLES INTERNES C------------------------------------------------------ C A: MODELE PLASTIQUE C IFIS = INT(REAL(VAR0(1))) IPLA = INT(REAL(VAR0(2))) SEQT = VAR0(3) SEQC = VAR0(4) DKT = VAR0(5) DKC = VAR0(6) ENDT = VAR0(7) ENDC = VAR0(8) IF (VAR0(9).EQ.0.D0) THEN ******* COERAND = RAND() ELSE COERAND = VAR0(9) ENDIF RT = (1. + ((2.*COERAND - 1.) * 0.1))*RT C IF (SEQT.EQ.0.D0) THEN SEQT = RT ENDIF IF (SEQC.EQ.0.D0) THEN SEQC = RC/3.D0 ENDIF C C--------------------------------------------------------------------- C GO TO (10,10,10,10,10,30),NSTRSS 10 CONTINUE KERRE=437 WRITE(*,*) '!! ATTENTION DANS BET3D NSTRSS=',NSTRSS STOP GO TO 1000 C 40 CONTINUE C 30 CONTINUE C C--------------------------------------------------------------------- C C C--------------------------------------------------------------------- C ACTUALISATION DES VARIABLES INTERNES C--------------------------------------------------------------------- C A: MODELE PLASTIQUE C VARF( 1)=FLOAT(IFIS) VARF( 2)=FLOAT(IPLA) VARF( 3)=SEQT VARF( 4)=SEQC VARF( 5)=DKT VARF( 6)=DKC VARF( 7)=ENDT VARF( 8)=ENDC VARF( 9)=COERAND C C--------------------------------------------------------------------- C 1000 CONTINUE C SEGSUP BETINSA C RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales