Télécharger la plaquette Cast3M


fr en

Télécharger cridom.eso

Retour à la liste

Numérotation des lignes : 

cridom

  1. C CRIDOM    SOURCE    CB215821  16/04/21    21:16:09     8920
  2.       SUBROUTINE CRIDOM(SIGG,NSTRS,VARR,XMAT,NCOMAT,XCAR,ICARA,
  3.      1INPLAS,NNVARI,MFR,XG,KERRE,SIGGD,Y)
  4. C
  5. C ===================================================================
  6. C CE SOUS-PROGRAMME EST APPELE PAR "ENDOM".
  7. C IL CALCULE LA VALEUR DU CRITERE D'ENDOMMAGEMENT POUR CHAQUE
  8. C MODELE DE MATERIAU ENDOMMAGEABLE.
  9. C
  10. C ENTREES:
  11. C -------
  12. C NSTRS        = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM.
  13. C SIGG(NSTRS)  = CONTR.
  14. C NNVARI       = NBR. DE VARIABLES INTERNES PILOTANT LES EQUATIONS
  15. C                DU MODELE
  16. C VARR(2)...   = VARIABLES INTERNES PILOTANT LES EQUATIONS DU MODELE
  17. C ...VARR(1+NNVARI)
  18. C NCOMAT       = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU
  19. C XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
  20. C MFR          = INDICE DE LA FORMULATION MECANIQUE; SEULEMENT
  21. C                MASSIF OU COQUE POUR LES MATERIAUX ENDOMMAGEABLES
  22. C ICARA        = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  23. C XCAR(ICARA)  = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  24. C INPLAS       = INDICE DU MODELE DECRIVANT LE MATERIAU ENDOMMAGEABLE
  25. C
  26. C SORTIES:
  27. C -------
  28. C XG= VALEUR DU CRITERE
  29. C SIGGD(NSTRS)= DEVIATEUR DES CONTR.
  30. C Y= SEUIL D'ENDOMMAGEMENT
  31. C KERRE= INDICE QUI REGIT LES ERREURS
  32. C      = 99 SI LA FORMULATION MECANIQUE N'EST PAS DISPONIBLE POUR
  33. C           LE MODELE CONSIDERE
  34. C
  35. C ==================================================================
  36. C
  37.       IMPLICIT INTEGER(I-N)
  38.       IMPLICIT REAL*8 (A-H,O-Z)
  39.       DIMENSION SIGG(*),VARR(*),XMAT(*),XCAR(*),SIGGD(*)
  40. C
  41.       KERRE=0
  42.       IPLAS=INPLAS-25
  43. C
  44. C -------------------------------------------------------------
  45. C LE CALCUL DU CRITERE DEPEND DU MODELE DONC DE INPLAS OU IPLAS
  46. C -------------------------------------------------------------
  47.       GOTO (10,1000,1000) IPLAS
  48. C
  49. C
  50. C ==================================================
  51. C ECROUISSAGE ET ENDOMMAGEMENT ISOTROPES - VON MISES
  52. C ==================================================
  53.    10 XP=VARR(1)
  54.       XNU=XMAT(2)
  55.       Y=XMAT(6)
  56. C
  57.       IF (MFR.NE.1) GOTO 50
  58. C
  59. C -------
  60. C MASSIFS
  61. C -------
  62. C
  63.       TRACE=SIGG(1)+SIGG(2)+SIGG(3)
  64. C
  65. C ..............................
  66. C CALCUL DU DEVIATEUR DES CONTR.
  67. C ..............................
  68.       CALL ZDANUL(SIGGD,NSTRS)
  69.       DO 30 I=1,NSTRS
  70.                        A=0.D0
  71.                        IF (I.LE.3) A=1.D0/3.D0
  72.                        SIGGD(I)=SIGG(I)-A*TRACE
  73.    30 CONTINUE
  74. C
  75. C ............................................
  76. C CALCUL DE LA CONTR. EQUIVALENTE DE VON MISES
  77. C ............................................
  78.       SIGEQI=PROCON(SIGGD,SIGGD,NSTRS)
  79.       SIGEQI=SQRT((3.D0/2.D0)*SIGEQI)
  80. C
  81. C .............................................................
  82. C                     CALCUL DU CRITERE
  83. C LORSQUE LE PREMIER PAS NON NUL DE L'UTILISATEUR EST PLASTIQUE
  84. C LE VON MISES RISQUE D'ETRE NUL POUR LE PAS PRECEDENT : POUR
  85. C EVITER LA DIVISION PAR ZERO ON RAJOUTE 1.D-10
  86. C .............................................................
  87.       SIGEQI=SIGEQI+1.D-10
  88.       TRIAXI=TRACE/(3.D0*SIGEQI)
  89.       TRIAXI=TRIAXI*TRIAXI
  90.       TRIAXI=3.D0*(1.D0-2.D0*XNU)*TRIAXI+(2.D0*(1.D0+XNU))/3.D0
  91.       XG=(XP*TRIAXI)-Y
  92.       GOTO 1000
  93. C
  94.    50 IF (MFR.NE.5) GOTO 998
  95. C
  96. C ---------------------------------------
  97. C COQUES EPAISSES - NON ENCORE DISPONIBLE
  98. C ---------------------------------------
  99. C
  100.   998 KERRE=99
  101.  1000 CONTINUE
  102.       RETURN
  103.       END
  104.  
  105.  
  106.  
  107.  

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