Télécharger la plaquette Cast3M


fr en

Télécharger mobflu.eso

Retour à la liste

Numérotation des lignes : 

mobflu

  1. C MOBFLU    SOURCE    CB215821  16/04/21    21:17:46     8920
  2. C MOBFLU     SOURCE    INSL       24/10/96
  3.       SUBROUTINE MOBFLU(SIGR,STFL,D1,NSTRS,IFLU,RMOY,PHIF,TEMP0,
  4.      & DTEMP,TEMP1,POAR)
  5. C
  6.       IMPLICIT INTEGER(I-N)
  7.       IMPLICIT REAL*8(A-H,O-Z)
  8.       DIMENSION SIGR(NSTRS),STFL(NSTRS),D1(NSTRS,NSTRS)
  9. C----------------------------------------------------------------------
  10. *
  11.       IF(IFLU.EQ.0) GOTO 10
  12.       IF(IFLU.EQ.10.OR.IFLU.EQ.11.OR.IFLU.EQ.20.OR.IFLU.EQ.21) GOTO 20
  13.       IF(IFLU.EQ.30.OR.IFLU.EQ.31) GOTO 30
  14. *
  15.       WRITE(*,*) '!! ATTENTION CE MODELE DE FLUAGE N EST PAS DISPONIBLE'
  16.       WRITE(*,*) '      (DANS MOBFLU)  IFLU= ',IFLU
  17.       STOP
  18. c+mdj
  19. C      GOTO 100
  20. c+mdj
  21. *
  22.    10 CONTINUE
  23. *
  24. *    CE MATERIAUX NE FLUE PAS
  25. *
  26.         GOTO 100
  27. *
  28.    20 CONTINUE
  29. *
  30. *    MODELE DE FLUAGE BETON : LOI BPEL91
  31. *
  32. *       IFLU = 10 : FLUAGE (BPEL91_K(t)..) +  BETON ELASTIQUE
  33. *       IFLU = 11 : FLUAGE (BPEL91_K(t)..) +  BETON ELASTOPLASTIQUE
  34. *       IFLU = 20 : FLUAGE (BPEL91_K(t)=2) +  BETON ELASTIQUE
  35. *       IFLU = 21 : FLUAGE (BPEL91_K(t)=2) +  BETON ELASTOPLASTIQUE
  36. *
  37.         TEMPF= TEMP0 + DTEMP
  38.         TEMPR= TEMPF - TEMP1
  39.         IF(IFLU.EQ.10.OR.IFLU.EQ.11) THEN
  40.           DKS  = 1.D0/(1.D0+(20.D0*POAR))
  41.           DKE  = 0.4D0
  42.           DKC  = (120.D0-PHIF)/30.D0+.6667D0*(100.D0-PHIF)/(20.D0+RMOY)
  43.           DKT1 = 100.D0/(100.D0+TEMP1)
  44.           FACT = DKS*(DKE+DKC*DKT1)
  45.         ELSE
  46.           FACT = PHIF
  47.         ENDIF
  48. *
  49.         COEF0 = SQRT(TEMPR)
  50.         COEF1 = 5.D0*SQRT(RMOY)
  51. *
  52. *       CALCUL EN DEFORMATION DE FLUAGE TOTALE
  53. *
  54.         COEFT = FACT*COEF0/(COEF0+COEF1)
  55. *
  56.         CALL BST(D1,SIGR,NSTRS,NSTRS,STFL)
  57.         DO 2 I=1,NSTRS
  58.            STFL(I)=COEFT*STFL(I)
  59.    2    CONTINUE
  60.         TEMP0 = TEMPF
  61.         GOTO 100
  62. *
  63.    30 CONTINUE
  64. *
  65. *    MODELE DE FLUAGE BETON :LOI1 ----> EPSfl = A1 * T**A2 * Sig/E
  66. *
  67. *        IFLU = 31 : FLUAGE (LOI1) +  BETON ELASTOPLASTIQUE
  68. *        IFLU = 30 : FLUAGE (LOI1) +  BETON ELASTIQUE
  69. *
  70.         TEMPF = TEMP0 + DTEMP
  71.         COEFT = RMOY*(TEMPF**PHIF)
  72.         CALL BST(D1,SIGR,NSTRS,NSTRS,STFL)
  73.         DO 3 I=1,NSTRS
  74.            STFL(I)=COEFT*STFL(I)
  75.    3    CONTINUE
  76.         TEMP0 = TEMPF
  77.         GOTO 100
  78. C
  79.   100 CONTINUE
  80. C---------------------------------------------------------------------
  81. 1991  FORMAT(18(1X,E12.5))
  82.       RETURN
  83.       END
  84.  
  85.  
  86.  
  87.  
  88.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales