Télécharger enplas.eso

Retour à la liste

Numérotation des lignes :

  1. C ENPLAS SOURCE LJASON 07/09/24 21:15:08 5866
  2. SUBROUTINE ENPLAS(XMAT,NMATT,VAR0,VARF,NVARI,SIG0,
  3. & SIGF,DEPST,NSTRS,KERRE, ISTEP)
  4.  
  5. *******************************************************
  6. * Variables *
  7. * *
  8. * XMAT : Variables matériau *
  9. * *
  10. * VAR0 : Variables internes à l'état initial *
  11. * VARF : Variables internes après incrément *
  12. * *
  13. * SIG0 : Contrainte à l'incrément précédent *
  14. * SIGF : Contrainte à l'incrément suivant *
  15. * *
  16. * DEPST : Incrément de déplacement total imposé *
  17. *******************************************************
  18.  
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8(A-H,O-Z)
  21.  
  22. *****************************
  23. * Déclaration des variables *
  24. *****************************
  25.  
  26. INTEGER NMATT, NVARI, NSTRS, KERRE, ISTEP
  27.  
  28. * Entrees
  29. REAL*8 XMAT(NMATT), VAR0(NVARI), SIG0(NSTRS)
  30. REAL*8 DEPST(NSTRS) , SIG02(6)
  31.  
  32. * Sorties
  33. REAL*8 SIGF(NSTRS), VARF(NVARI)
  34.  
  35. * Matrice tangente
  36. REAL*8 DSDE(6,6)
  37.  
  38. * Variables materiau
  39. REAL*8 E,NU
  40. REAL*8 AC, AT, BC, BT, EPSD0, RC, RT
  41. REAL*8 P, AH, BH, CH, GAMMA, ALPHA
  42. REAL*8 A, K0, BETA
  43.  
  44. * Transport des deformations
  45. DO I = 4,NSTRS
  46. DEPST(I) = DEPST(I)/2.d0
  47. END DO
  48.  
  49. DO I=1,NSTRS
  50. SIG02(I) = SIG0(I)
  51. END DO
  52.  
  53. ********************
  54. * Partie plastique *
  55. ********************
  56.  
  57. CALL PPLAST(XMAT,NMATT,VAR0,VARF,NVARI,SIG0,SIGF,DEPST,
  58. & NSTRS,KERRE)
  59.  
  60. DO I = 1,NSTRS
  61. SIG0(I) = SIGF(I)
  62. END DO
  63.  
  64. ************************
  65. * Partie endommagement *
  66. ************************
  67.  
  68. * Appel a la routine MAZARS modifiee
  69.  
  70. CALL NMAZAR(XMAT,NMATT,VAR0,VARF,NVARI,SIG0,SIGF,
  71. & DEPST,NSTRS,ISTEP)
  72.  
  73. DO I = 4,NSTRS
  74. DEPST(I)=2.d0*DEPST(I)
  75. END DO
  76.  
  77. IF (ISTEP.EQ.1) THEN
  78. DO I=1,NSTRS
  79. SIGF(I)=SIG02(I)
  80. END DO
  81. END IF
  82.  
  83. IF (KERRE.NE.0.d0) THEN
  84. WRITE(*,*) KERRE
  85. END IF
  86.  
  87. RETURN
  88. END
  89.  
  90.  

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