Télécharger la plaquette Cast3M


fr en

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