Télécharger prgurs.eso

Retour à la liste

Numérotation des lignes :

  1. C PRGURS SOURCE CHAT 05/01/13 02:26:34 5004
  2. SUBROUTINE PRGURS(SIG0,NSTRS,DEPST,VAR0,XMAT,NCOMAT,XCAR,ICARA,
  3. 1NVARI,SIGF,VARF,DEFP,MFR1,KERRE,wrkgur)
  4. C
  5. C ==================================================================
  6. C CE SOUS-PROGRAMME EST APPELE DANS "ECOUL2".
  7. C IL PREPARE L'INTEGRATION DE LA LOI DE GURSON
  8. C
  9. C ENTREES:
  10. C -------
  11. C NSTRS = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM.
  12. C SIG0(NSTRS) = CONTR. AU DEBUT DU PAS D'INTEGRATION
  13. C DEPST(NSTRS) = INCREMENT DES DEFORM. TOTALES
  14. C NVARI = NBR. DE VARIABLES INTERNES
  15. C VAR0(NVARI) = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  16. C CE TABLEAU CONTIENT DANS L'ORDRE:
  17. C - LA DEFORM. PLASTIQUE CUMULEE
  18. C - LA FRACTION DE VIDE INTERNE
  19. C - LA DENSITE
  20. C - LES DEFORM. PLASTIQUES
  21. C NCOMAT = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU
  22. C XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
  23. C MFR = INDICE DE LA FORMULATION MECANIQUE; SEULEMENT
  24. C MASSIF OU COQUE POUR LES MATERIAUX ENDOMMAGEABLES
  25. C ICARA = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  26. C XCAR(ICARA) = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  27. C
  28. C SORTIES:
  29. C -------
  30. C SIGF(NSTRS)= CONTR. A LA FIN DU PAS D'INTEGRATION
  31. C VARF(NVARI)= VARIABLES INTERNES A LA FIN DU PAS D'INTEGRATION
  32. C DEFP(NSTRS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS
  33. C D'INTEGRATION
  34. C KERRE = INDICE QUI REGIT LES ERREURS
  35. C = 77 SI LA DEFORM. PLAST. CUMULEE ENDOMMAGEE (2IEME VAR.
  36. C INT.) EST EN DEHORS DE LA COURBE DE TRACTION, DS.
  37. C LE CAS DE L'ECROUISSAGE ET DE L'ENDOMM. ISOTROPES.
  38. C CECI PEUT SE PRODUIRE SUITE A L'APPEL A "CRIDAM"
  39. C = 99 SI LA FORMULATION MECANIQUE N'EST PAS DISPONIBLE
  40. C POUR LE MODELE CONSIDERE OU S'IL Y A INCOMPATIBILITE
  41. C ENTRE MFR ET IFOUR
  42. C ==================================================================
  43. C ICI IL FAUT PROGRAMMER EN FORTRAN PUR
  44. C ===================================================================
  45. C
  46. IMPLICIT INTEGER(I-N)
  47. IMPLICIT REAL*8(A-H,O-Z)
  48. -INC PPARAM
  49. -INC CCOPTIO
  50. segment wrkgur
  51. real*8 sigbar, sy0,phi0,rho0,g,b,h
  52. real*8 epn,phin,sqrtj2,rho,sig(6)
  53. real*8 e(7),dt
  54. real*8 conv,tol1,tol2
  55. endsegment
  56. * common/prop/sigbar,sy0,phi0,rho0,g,b,h
  57. DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),SIGF(*),
  58. & VARF(*),DEFP(*)
  59. DIMENSION RSIG0(6),RDEFP(6),RSIGF(6) ,RDEPST(6)
  60. C
  61. C ADAPTATION DE L'OPTION DE CALCUL VERS LE 3D MASSIF DE SIG0 A RSIG0
  62. C
  63. IF (MFR1 .EQ. 1) THEN
  64. C MASSIF 3D
  65. IF (NSTRS .EQ. 6) THEN
  66. DO 10 I=1,NSTRS
  67. RSIG0(I)=SIG0(I)
  68. RDEPST(I)=DEPST(I)
  69. 10 CONTINUE
  70. ELSE IF ( NSTRS .EQ. 4 .AND. (IFOUR .EQ. 0 .OR.
  71. & IFOUR .EQ. -1)) THEN
  72. C calcul axisyemmtrique ou deformation plane
  73. DO 15 I=1,NSTRS
  74. RSIG0(I)=SIG0(I)
  75. RDEPST(I)=DEPST(I)
  76. 15 CONTINUE
  77. RSIG0(5)=0.D0
  78. RSIG0(6)=0.D0
  79. RDEPST(5)=0.D0
  80. RDEPST(6)=0.D0
  81. ENDIF
  82. ELSE
  83. KERRE = 99
  84. RETURN
  85. ENDIF
  86. C
  87. C DONNE MATERIAU
  88. C
  89. SIGBAR=XMAT(7)
  90. SY0=XMAT(5)
  91. phi0=XMAT(8)
  92. G=XMAT(1)/2.D0/(1.D0+XMAT(2))
  93. B=XMAT(1)/3.D0/(1.D0-2.D0*XMAT(2))
  94. H=XMAT(6)
  95. RHO0=1.D0
  96. C
  97. C METTRE LES CONTRAINTES SOUS FORME DEVIATORIC + PRESSION
  98. C
  99. DO 20 I=1,6
  100. RSIGF(I)=RSIG0(I)
  101. 20 CONTINUE
  102. C
  103. RP = -(RSIG0(1)+RSIG0(2)+RSIG0(3))/3.D0
  104. RSIG0(1)=RSIG0(1)+RP
  105. RSIG0(2)=RSIG0(2)+RP
  106. RSIG0(3)=RSIG0(3)+RP
  107.  
  108. * WRITE (*,*) '------------------------------'
  109. * WRITE (*,100) 'sig=' ,(RSIG0(I),I=1,6)
  110. * WRITE (*,*) 'P=',RP
  111. * WRITE (*,100) 'depst=' ,(DEPST(I),I=1,6)
  112. * WRITE (*,*) 'PHI0=',PHI0
  113. * WRITE (*,*) 'epse=',VAR0(1), ' phi=',VAR0(2),' dens=',VAR0(3)
  114.  
  115. C
  116. C NOUVEL ETAT APRES L'INCREMENT DE DEFORMATION
  117. C
  118. CALL GURSON(RSIG0,RP,RDEFP,VAR0(1),VAR0(2),VAR0(3),
  119. & RDEPST,1.D0,wrkgur)
  120.  
  121. C
  122. * WRITE (*,100) 'sig=' ,(RSIG0(I),I=1,6)
  123. * WRITE (*,*) 'P=',RP
  124. * WRITE (*,*) 'epse=',VAR0(1),' phi0=',VAR0(2),' dens=',VAR0(3)
  125.  
  126. RSIG0(1)=RSIG0(1)-RP
  127. RSIG0(2)=RSIG0(2)-RP
  128. RSIG0(3)=RSIG0(3)-RP
  129. C
  130. C MISE A JOUR DES VARIABLES
  131. C
  132. C VARIABLES INTERNES
  133. C
  134. DO 50 I=1,NVARI
  135. VARF(I)=VAR0(I)
  136. 50 CONTINUE
  137. C
  138. C CONTRAINTES
  139. C
  140. DO 60 I=1,6
  141. RSIGF(I)=RSIG0(I)
  142. 60 CONTINUE
  143. C
  144. C PASSAGE A L'OPTION DE CALCUL POUR LES CONTRAINTES
  145. C
  146. IF (MFR1 .EQ. 1) THEN
  147. IF (NSTRS .EQ. 6) THEN
  148. C massif 3D
  149. DO 70 I=1,NSTRS
  150. SIGF(I)=RSIGF(I)
  151. DEFP(I)=RDEFP(I)
  152. 70 CONTINUE
  153. ELSE IF ( NSTRS .EQ. 4 .AND. (IFOUR .EQ. 0
  154. & .OR. IFOUR .EQ. -1)) THEN
  155. C calcul axisymetrique ou deformations planes
  156. DO 80 I=1,NSTRS
  157. SIGF(I)=RSIGF(I)
  158. DEFP(I)=RDEFP(I)
  159. 80 CONTINUE
  160. IF ( ABS(RSIGF(5)) .GT. 1.E-6 .OR.
  161. & ABS(RSIGF(6)) .GT. 1.E-6) THEN
  162. * PRINT* ,'Incompatibilite de la valeur des contraintes'
  163. KERRE = 99
  164. * CALL ERREUR(5)
  165. ENDIF
  166. ENDIF
  167. ENDIF
  168. C
  169. C WRITE (*,100) 'sig=' ,(SIGF(I),I=1,6)
  170. C WRITE (*,100) 'var=' ,(VARF(I),I=1,3)
  171. C WRITE (*,100) 'dep=' ,(DEFP(I),I=1,6)
  172. C
  173. 100 FORMAT (A,6(1PE14.7,1X))
  174. RETURN
  175. END
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  

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