Télécharger la plaquette Cast3M


fr en

Télécharger cridam.eso

Retour à la liste

Numérotation des lignes : 

cridam

  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