Télécharger cridam.eso

Retour à la liste

Numérotation des lignes :

  1. C CRIDAM SOURCE CB215821 16/04/21 21:16:09 8920
  2. SUBROUTINE CRIDAM(SIGG,NSTRS,VARR,XMAT,NCOMAT,XCAR,ICARA,
  3. 1INPLAS,NTRAC,TRACC,NC,NPNTER,NNVARI,MFR,XF,KERRE,SIGGD,
  4. 2YY,T,TRAC,NCOURB)
  5. C
  6. C ===================================================================
  7. C CE SOUS-PROGRAMME EST APPELE PAR "ENDOM".
  8. C IL CALCULE LA VALEUR DU CRITERE POUR CHAQUE MODELE DE MATERIAU
  9. C ENDOMMAGEABLE.
  10. C
  11. C ENTREES:
  12. C -------
  13. C NSTRS = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM.
  14. C SIGG(NSTRS) = CONTR.
  15. C NNVARI = NBR. DE VARIABLES INTERNES PILOTANT LES EQUATIONS
  16. C DU MODELE
  17. C VARR(2)... = VARIABLES INTERNES PILOTANT LES EQUATIONS DU MODELE
  18. C ...VARR(1+NNVARI)
  19. C NCOMAT = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU
  20. C XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
  21. C MFR = INDICE DE LA FORMULATION MECANIQUE; SEULEMENT
  22. C MASSIF OU COQUE POUR LES MATERIAUX ENDOMMAGEABLES
  23. C ICARA = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  24. C XCAR(ICARA) = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  25. C INPLAS = INDICE DU MODELE DECRIVANT LE MATERIAU ENDOMMAGEABLE
  26. C NTRAC = DIMENSION DE TRACC
  27. C TRACC(NTRAC) = COURBES DE TRACTION
  28. C NC = NBR. DE COURBES DE TRACTION
  29. C NPNTER(NC)= TABLEAU CONTENANT LE NOMBRE DE POINTS POUR CHAQUE
  30. C COURBE DE TRACTION (A T FIXEE)
  31. C T = TEMPERATURE A LAQUELLE ON CALCULE LE CRITERE
  32. C
  33. C SORTIES:
  34. C -------
  35. C XF= VALEUR DU CRITERE
  36. C KERRE= INDICE QUI REGIT LES ERREURS
  37. C = 77 SI LA DEFORM. PLAST. CUMULEE ENDOMMAGEE (VARIABLE VIR)
  38. C EST EN DEHORS DE LA COURBE DE TRACTION, DS. LE CAS DE
  39. C L'ECROUISSAGE ET DE L'ENDOMMAGEMENT ISOTROPES
  40. C = 99 SI LA FORMULATION MECANIQUE N'EST PAS DISPONIBLE POUR
  41. C LE MODELE CONSIDERE
  42. C
  43. C ==================================================================
  44. C
  45. IMPLICIT INTEGER(I-N)
  46. IMPLICIT REAL*8 (A-H,O-Z)
  47. DIMENSION SIGG(*),VARR(*),XMAT(*),XCAR(*),TRACC(*),SIGGD(*),
  48. 1 TRAC(*),NPNTER(*)
  49. C
  50. KERRE=0
  51. IPLAS=INPLAS-25
  52. C
  53. C -------------------------------------------------------------
  54. C LE CALCUL DU CRITERE DEPEND DU MODELE DONC DE INPLAS OU IPLAS
  55. C -------------------------------------------------------------
  56. GOTO (10,1000,1000) IPLAS
  57. C
  58. C
  59. C ==================================================
  60. C ECROUISSAGE ET ENDOMMAGEMENT ISOTROPES - VON MISES
  61. C ==================================================
  62. 10 VIR=VARR(2)
  63. XD=VARR(3)
  64. C
  65. IF (MFR.NE.1) GOTO 50
  66. C
  67. C -------
  68. C MASSIFS
  69. C -------
  70. C
  71. CALL TRACTT(TRACC,NPNTER,TRAC,NTRAC,NC,NCOURB,T)
  72. N2COUR=2*NCOURB
  73. CALL DERTRA(N2COUR,TRAC,VIR,YY,YPRIM,XINF,XSUP)
  74. C
  75. TRACE=SIGG(1)+SIGG(2)+SIGG(3)
  76. C
  77. C ..............................
  78. C CALCUL DU DEVIATEUR DES CONTR.
  79. C ..............................
  80. CALL ZDANUL(SIGGD,NSTRS)
  81. DO 30 I=1,NSTRS
  82. A=0.D0
  83. IF (I.LE.3) A=1.D0/3.D0
  84. SIGGD(I)=SIGG(I)-A*TRACE
  85. 30 CONTINUE
  86. C
  87. C ............................................
  88. C CALCUL DE LA CONTR. EQUIVALENTE DE VON MISES
  89. C ............................................
  90. SIGEQI=PROCON(SIGGD,SIGGD,NSTRS)
  91. SIGEQI=SQRT((3.D0/2.D0)*SIGEQI)
  92. C
  93. C .................
  94. C CALCUL DU CRITERE
  95. C .................
  96. XF=SIGEQI/(1.D0-XD)-YY
  97. GOTO 1000
  98. C
  99. 50 IF (MFR.NE.5) GOTO 998
  100. C
  101. C ---------------------------------------
  102. C COQUES EPAISSES - NON ENCORE DISPONIBLE
  103. C ---------------------------------------
  104. C
  105. 998 KERRE=99
  106. 1000 CONTINUE
  107. RETURN
  108. END
  109.  
  110.  
  111.  
  112.  
  113.  

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