Télécharger la plaquette Cast3M


fr en

Télécharger cameca.eso

Retour à la liste

Numérotation des lignes : 

cameca

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

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