Télécharger incvar.eso

Retour à la liste

Numérotation des lignes :

  1. C INCVAR SOURCE PV 15/04/13 21:15:10 8474
  2. C INCREP SOURCE ROUX G-M 10/11/04
  3. SUBROUTINE INCVAR(SIG,VAR,EPSVPT,VARPT,XMAT,NSTRS0,
  4. & MFR,NVARI,NCOMAT,INPLAS)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION SIG(*),VAR(*),EPSVPT(*),VARPT(*),XMAT(*)
  7. DIMENSION XX(6),ST(3,3),D(3),XP(3,3)
  8.  
  9. C***********************************************
  10. C CALCUL MECANIQUE
  11. C***********************************************
  12.  
  13. C------------------------------
  14. C Initialisation des variables
  15. C------------------------------
  16. IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN
  17. IF (INPLAS.EQ.137) THEN
  18. C1=XMAT(11)
  19. gam1=XMAT(12)
  20. C2=XMAT(13)
  21. ELSE
  22. C1=XMAT(10)
  23. gam1=XMAT(11)
  24. C2=XMAT(12)
  25. ENDIF
  26. ENDIF
  27.  
  28. C-----------------
  29. C Calcul de J2(ss)
  30. C-----------------
  31. C calcul de ss, déviateur des contraintes
  32. TRACE =(SIG(1)+SIG(2)+SIG(3))/3.0D0
  33. DO 70 I=1,NSTRS0
  34. A = 0.0D0
  35. IF (I.LE.3) A=1.0D0
  36. IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN
  37. XX(I) = SIG(I)-A*TRACE
  38. ENDIF
  39. IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN
  40. XXX = SIG(I)-A*TRACE-VAR(I)*C1
  41. XX(I) = XXX-VAR(I+NSTRS0)*C2
  42. ENDIF
  43. 70 CONTINUE
  44. AJ2 = PROCON (XX,XX,NSTRS0)
  45. AJ2 = SQRT(1.5D0*AJ2)
  46.  
  47. C------------------
  48. C test d'élasticité
  49. C------------------
  50. IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN
  51. RR=VAR(1)*XMAT(6)*XMAT(7)
  52. ENDIF
  53. IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN
  54. RR=VAR((2*NSTRS0)+1)*XMAT(6)*XMAT(7)
  55. ENDIF
  56. SY=XMAT(5)
  57. EE=AJ2-SY-RR
  58. IF (EE.GT.0D0) THEN
  59. TESTE=0.0D0
  60. ELSE
  61. TESTE=1.0D0
  62. ENDIF
  63.  
  64. C-------------------------------------------------
  65. C vitesse de déformation plastique cumulée ppt=rpt
  66. C-------------------------------------------------
  67. IF (TESTE.EQ.1.0D0) THEN
  68. RPT=0.0D0
  69. ELSE
  70. IF (INPLAS.EQ.136.OR.INPLAS.EQ.137) THEN
  71. RPT1=EE/XMAT(8)
  72. rpt1=max(-335d0,min(335d0,rpt1))
  73. RPT2=(exp(RPT1)-exp(-1.D0*RPT1))/2.d0
  74. RPT=XMAT(10)*(RPT2)**XMAT(9)
  75. ENDIF
  76. IF (INPLAS.EQ.138.OR.INPLAS.EQ.139) THEN
  77. RPT=(EE/XMAT(8))**XMAT(9)
  78. ENDIF
  79. ENDIF
  80. IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN
  81. VARPT(2)=RPT
  82. ENDIF
  83. IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN
  84. VARPT((2*NSTRS0)+2)=RPT
  85. ENDIF
  86.  
  87. C----------------------------------
  88. C vitesses de déformation plastique
  89. C----------------------------------
  90. IF (TESTE.EQ.1.0D0) THEN
  91. DO 80 I=1,NSTRS0
  92. EPSVPT(I)=0.0D0
  93. 80 CONTINUE
  94. ELSE
  95. DO 90 I=1,NSTRS0
  96. EPSVPT(I)=(1.5D0*RPT*XX(I))/AJ2
  97. 90 CONTINUE
  98. ENDIF
  99.  
  100. C----------------------------------------------
  101. C vitesse d'évolution de l'écrouissage isotrope
  102. C----------------------------------------------
  103. IF (INPLAS.EQ.136.OR.INPLAS.EQ.138) THEN
  104. Dr=(1-XMAT(6)*VAR(1))*RPT
  105. VARPT(1)=Dr
  106. ENDIF
  107. IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN
  108. Dr=(1-XMAT(6)*VAR((2*NSTRS0)+1))*RPT
  109. VARPT((2*NSTRS0)+1)=Dr
  110. ENDIF
  111.  
  112. C-------------------------------------------------
  113. C vitesse d'évolution de l'écrouissage cinematique
  114. C-------------------------------------------------
  115. IF (INPLAS.EQ.137.OR.INPLAS.EQ.139) THEN
  116. C======NON-LINEAIRE
  117. DO 96 I=1,NSTRS0
  118. VARPT(I)=EPSVPT(I)-gam1*RPT*VAR(I)
  119. 96 CONTINUE
  120. C======NON-LINEAIRE
  121. DO 97 I=1,NSTRS0
  122. VARPT(I+NSTRS0)=EPSVPT(I)
  123. 97 CONTINUE
  124. ENDIF
  125.  
  126.  
  127.  
  128.  
  129. RETURN
  130. END
  131.  
  132.  
  133.  

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