Télécharger endo.eso

Retour à la liste

Numérotation des lignes :

endo
  1. C ENDO SOURCE AM 09/12/08 21:15:13 6582
  2. SUBROUTINE ENDO(EPSILO,NSTRS,NVARI,CMAT,DEP,SIGF,VAR0,
  3. * VARF,NMAT,IFOUR)
  4. C
  5. C_________________________________________________________________
  6. C
  7. C CALCUL DE L ENDOMMAGEMENT ET DES CONTRAINTES APPARENTES FINALES
  8. C_________________________________________________________________
  9. C
  10. C declaration des variables
  11. C
  12. C
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. DIMENSION EPS33(3,3),EPSIPP(3),EPSILT(3),VALP33(3,3)
  16. DIMENSION SIGP(3),SIGPT(3),SIGPC(3)
  17. DIMENSION EPSILO(NSTRS),SIGF(NSTRS),CMAT(NMAT)
  18. DIMENSION DEP(NSTRS,NSTRS),VARF(NVARI),VAR0(NVARI)
  19. C
  20. XZERO=0.D0
  21. UN=1.D0
  22. XPETIT=1.D-12
  23. C
  24. C RECUPERATION DES PARAMETRES MATERIAUX
  25. C ET DES VARIABLES INTERNES
  26. C
  27. YOUN = CMAT(1)
  28. XNU = CMAT(2)
  29. EPSD0= CMAT(5)
  30. ATRA = CMAT(6)
  31. BTRA = CMAT(7)
  32. ACOM = CMAT(8)
  33. BCOM = CMAT(9)
  34. DINI = VAR0(2)
  35. C
  36. C CHANGEMENT EPSILO
  37. C CONTRAINTES PLANES
  38. C
  39. IF (IFOUR.EQ.-2) THEN
  40. temp = EPSILO(3)
  41. EPSILO(3) = EPSILO(4)
  42. EPSILO(4) = temp
  43. ENDIF
  44. C
  45. C REECRITURE DU VECTEUR DEFORMATION SOUS FORME MATRICIELLE
  46. C ENTREE : EPSILO -> VECTEUR DEFORMATION ELASTIQUE
  47. C SORTIE : EPS33 -> MATRICE DEFORMATION ELASTIQUE
  48. C Le 2 correspond au cisaillement
  49. C
  50. CALL ENDOCB (EPSILO,EPS33,2,IFOUR)
  51. C
  52. C DIAGONALISATION DE LA MATRICE DEFORMATION ELASTIQUE EPS33
  53. C SORTIE : EPSIPP(3) - VALEURS PROPRES
  54. C VALP33(3,3) - VECTEURS PROPRES
  55. C
  56. IF (IFOUR.EQ.2) THEN
  57. CALL JACOB3 (EPS33,3,EPSIPP,VALP33)
  58. ELSE
  59. CALL JACOB3 (EPS33,2,EPSIPP,VALP33)
  60. ENDIF
  61. C
  62. C CALCUL DES CONTRAINTES PRINCIPALES
  63. C
  64. DO 10 I=1,3
  65. SIGP(I)= XZERO
  66. DO 10 J=1,3
  67. SIGP(I)=SIGP(I)+DEP(I,J)*EPSIPP(J)
  68. 10 CONTINUE
  69. C
  70. C ON COMPLETE LA DEFORMATION DANS LE CAS DES
  71. C CONTRAINTES PLANES
  72. C
  73. IF (IFOUR.EQ. -2) THEN
  74. EPSIPP(3)= -(EPSIPP(1) + EPSIPP(2))*XNU / (UN-XNU)
  75. ENDIF
  76. C
  77. C CALCUL DE EPSILON TILDA (DEFORMATION POSITIVE)
  78. C
  79. EPSTIL=MAX( XZERO , EPSIPP(1) )**2 +
  80. * MAX( XZERO , EPSIPP(2) )**2 +
  81. * MAX( XZERO , EPSIPP(3) )**2
  82. EPSTIL=SQRT (EPSTIL)
  83. C
  84. C TEST EPSTIL>EPSD0
  85. C
  86. IF ( EPSTIL .GT. EPSD0) THEN
  87. C
  88. C CALCUL DE L ENDOMMAGEMENT
  89. C
  90. C CALCUL DE LA TRACE DES CONTRAINTES
  91. C NEGATIVE (COMPRESSION) : TRSIGC
  92. C POSITIVE (TRACTION) : TRSIGT
  93. C
  94. DO 20 I=1,3
  95. IF (SIGP(I).LT. XZERO) THEN
  96. SIGPC(I) = SIGP(I)
  97. SIGPT(I) = XZERO
  98. ELSE
  99. SIGPT(I) = SIGP(I)
  100. SIGPC(I) = XZERO
  101. ENDIF
  102. 20 CONTINUE
  103. C
  104. TRSIGT = SIGPT(1)+SIGPT(2)+SIGPT(3)
  105. TRSIGC = SIGPC(1)+SIGPC(2)+SIGPC(3)
  106. C
  107. C CALCUL DES DEFORMATIONS DUES AU CONTRAINTES POSITIVES
  108. C
  109. DO 30 I=1,3
  110. EPSILT(I) = (SIGPT(I)*(UN+XNU)-TRSIGT*XNU)/YOUN
  111. 30 CONTINUE
  112. C
  113. C CALCUL DE ALFAT ET ALFAC
  114. C
  115. ALFAT = MAX(XZERO,EPSIPP(1))*EPSILT(1) +
  116. * MAX(XZERO,EPSIPP(2))*EPSILT(2) +
  117. * MAX(XZERO,EPSIPP(3))*EPSILT(3)
  118. ALFAT = ALFAT/(EPSTIL*EPSTIL)
  119. ALFAC = UN - ALFAT
  120. C
  121. C AMELIORATION DE LA REPONSE EN CISAILLEMENT
  122. C POUR BETA > 1
  123. C
  124. BETA = UN
  125. IF (BETA .GT. UN) THEN
  126. IF ( ALFAT .GT. XPETIT ) THEN
  127. ALFAT=ALFAT**BETA
  128. ENDIF
  129. IF ( ALFAC .GT. XPETIT ) THEN
  130. ALFAC=ALFAC**BETA
  131. ENDIF
  132. ENDIF
  133. C
  134. c CALCUL DE DT, DC ET D
  135. C
  136. DT = UN - (EPSD0/EPSTIL)*(
  137. * (UN + ATRA)*EXP(-BTRA*(EPSTIL-EPSD0)) -
  138. * ATRA*EXP(-2.D0*BTRA*(EPSTIL-EPSD0)))
  139. DC = UN - (EPSD0/EPSTIL)*(
  140. * (UN + ACOM)*EXP(-BCOM*(EPSTIL-EPSD0)) -
  141. * ACOM*EXP(-2.D0*BCOM*(EPSTIL-EPSD0)))
  142. D = ALFAT*DT + ALFAC*DC
  143. C
  144. C ON BORNE LA VALEUR DE D a 0.99999
  145. C
  146. D = MIN ( D , UN-1.D-05 )
  147. C
  148. ELSE
  149. D = XZERO
  150. ENDIF
  151. C
  152. C ON TESTE LA CROISSANCE DE D
  153. C
  154. D=MAX ( D , DINI )
  155. C
  156. C STOCKAGE DES VARIABLES INTERNES FINALES
  157. C
  158. VARF(2) = D
  159. C
  160. C CONTRAINTES APPARENTES
  161. C
  162. DO 40 I=1,NSTRS
  163. SIGF(I) = SIGF(I)*(1.D0-D)
  164. 40 CONTINUE
  165. C
  166. C CHANGEMENT EPSILO
  167. C CONTRAINTES PLANES
  168. C
  169. IF (IFOUR.EQ.-2) THEN
  170. temp = EPSILO(3)
  171. EPSILO(3) = EPSILO(4)
  172. EPSILO(4) = temp
  173. ENDIF
  174. C
  175. END
  176.  
  177.  
  178.  

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