Télécharger cameca.eso

Retour à la liste

Numérotation des lignes :

  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