Télécharger la plaquette Cast3M


fr en

Télécharger camecthm.eso

Retour à la liste

Numérotation des lignes : 

camecthm

  1. C CAMECTHM  SOURCE    FD218221  26/02/13    21:15:07     12462          
  2. C CAMECTHM    SOURCE    AM        09/12/08    21:15:01     6582
  3.       SUBROUTINE CAMECTHM(NSTRS,DEPST,VAR0,NVARI,CMAT,NMAT,
  4.      *            IFOUR,IB,IGAU,SIGF,VARF,
  5.      * SIGI,DEP,EPSILO,EPSEQU,EPSFLU,DSIGT,SIGEF,AN,BN,CN,
  6.      * KERRE)
  7.  
  8. C
  9. C ROUTINE APPELE DANS BETONTHM
  10. C___________________________________________________________________
  11. C
  12. C             MODELE BETON FLUAGE-ENDOMMAGMENT
  13. C                   
  14. C       MODELE G. Sciumè (readatation M. Reviron-Benboudjema)
  15. C___________________________________________________________________
  16. C
  17.       IMPLICIT INTEGER(I-N)
  18.       IMPLICIT REAL*8(A-H,O-Z)
  19.       DIMENSION CMAT(NMAT)
  20.       DIMENSION VAR0(NVARI),VARF(NVARI)
  21. *
  22.       DIMENSION SIGI(NSTRS),DEP(NSTRS,NSTRS)
  23.       DIMENSION DEPST(NSTRS),EPSILO(NSTRS)
  24.       DIMENSION EPSEQU(NSTRS),EPSFLU(NSTRS)
  25.       DIMENSION DSIGT(NSTRS),SIGF(NSTRS),SIGEF(NSTRS)
  26.       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  27. C
  28. C______________________________________________________
  29. C Recuperation variable & intialisation
  30. C______________________________________________________
  31. C
  32.       CALL ZERO(AN,NSTRS,1)
  33.       CALL ZERO(BN,NSTRS,NSTRS)
  34.       CALL ZERO(CN,NSTRS,NSTRS)
  35.       CALL ZERO(EPSFLU,NSTRS,1)
  36. C
  37.       TDT  = CMAT(28)
  38.       BETA = CMAT(12)
  39. C
  40. C RECUPERATION DES VARIABLES INTERNES INITIALES
  41. C
  42.       SIGEF(1) = VAR0(3)
  43.       SIGEF(2) = VAR0(4)
  44. C
  45. C RECUPERATION DES CONTRAINTES INITIALES EFFECTIVES
  46. C
  47.       IF (IFOUR.EQ.-2) THEN
  48.         SIGEF(3) = VAR0(5)
  49.       ELSE
  50.         SIGEF(3) = VAR0(6)
  51.         SIGEF(4) = VAR0(5)
  52.         IF (IFOUR.EQ.2) THEN
  53.           SIGEF(5) = VAR0(7)
  54.           SIGEF(6) = VAR0(8)
  55.         ENDIF
  56.       ENDIF
  57. C
  58. C_______________________________________________________
  59. C
  60. C CONSTRUCTION DES MATRICES DE COMPLAISANCE ET ELASTIQUE
  61. C_______________________________________________________
  62. C
  63. C Matrice elastique
  64. C
  65.       CALL MATLAS(CMAT,NMAT,NSTRS,IFOUR,DEP,KERRE)
  66. C
  67. C Matrice de fluage
  68. C
  69. C   DFLUT = AN + BN * SIGMA_N + CN * SIGMA_N+1
  70. C
  71.       IF (TDT.NE.0.D0) THEN
  72.         CALL DEFCR(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  73.      $              AN,BN,CN,KERRE)
  74.       ENDIF
  75. C
  76. C_______________________________________________
  77. C
  78. C  CALCUL DE L INCREMENT DE CONTRAINTE EFFECTIVE
  79. C_______________________________________________
  80. C
  81.       CALL ELASF1(DEP,AN,BN,CN,DEPST,SIGEF,
  82.      $             NSTRS,DSIGT,KERRE)
  83. C
  84. C___________________________________________
  85. C
  86. C  CALCUL DES CONTRAINTES EFFECTIVES FINALES
  87. C___________________________________________
  88. C
  89.       DO 10 I=1,NSTRS
  90.         SIGF(I) = SIGEF(I) + DSIGT(I)
  91.    10 CONTINUE
  92. C__________________________________________________
  93. C
  94. C SAUVEGARDE DES DEFORMATIONS DE FLUAGE
  95. C  ET CALCUL DE LA DEFORMATION DE FLUAGE
  96. C   A LA FIN DU PAS DE TEMPS EPSFLU
  97. C__________________________________________________
  98. C
  99.       IF (TDT.NE.0.D0) THEN
  100.         CALL SAVCR(VARF,NVARI,CMAT,NMAT,SIGF,
  101.      $     SIGEF,NSTRS,IFOUR,EPSFLU,AN,BN,CN,KERRE)
  102.       ENDIF
  103. C__________________________________________________
  104. C
  105. C SAUVEGARDE DES CONTRAINTES EFFECTIVES
  106. C__________________________________________________
  107. C
  108.       VARF(3) = SIGF(1)
  109.       VARF(4) = SIGF(2)
  110.       IF (IFOUR.EQ.-2) THEN
  111.         VARF(5) = SIGF(3)
  112.       ELSE
  113.         VARF(6) = SIGF(3)
  114.         VARF(5) = SIGF(4)
  115.         IF (IFOUR.EQ.2) THEN
  116.           VARF(7) = SIGF(5)
  117.           VARF(8) = SIGF(6)
  118.         ENDIF
  119.       ENDIF
  120. C__________________________________________________
  121. C
  122. C  CALCUL DE LA DEFORMATION ELASTIQUE EN FIN DE PAS
  123. C__________________________________________________
  124. C
  125.       CALL DEFELA(DEP,EPSILO,NSTRS,SIGF,IFOUR,KERRE)
  126. C
  127. C__________________________________________________
  128. C
  129. C  CALCUL DE LA DEFORMATION EQUIVALENTE
  130. C__________________________________________________
  131. C
  132.       DO 20 I=1,NSTRS
  133.         EPSEQU(I) = EPSILO(I) + ( BETA * EPSFLU(I) )
  134.    20 CONTINUE
  135. C
  136. C____________________________________________________________
  137. C
  138. C  CALCUL ENDOMMAGEMENT ET DES CONTRAINTES APPARENTES FINALES
  139. C____________________________________________________________
  140. C
  141. C ENDOMMAGEMENT ET CONTRAINTES APPARENTES
  142. C
  143.       CALL ENDOTHM(EPSEQU,NSTRS,NVARI,CMAT,DEP,SIGF,VAR0,
  144.      $          VARF,NMAT,IFOUR)
  145. C
  146.       END
  147.  
  148.  

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