Télécharger cridom.eso

Retour à la liste

Numérotation des lignes :

  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