Télécharger la plaquette Cast3M


fr en

Télécharger enplas.eso

Retour à la liste

Numérotation des lignes : 

enplas

  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