Télécharger la plaquette Cast3M


fr en

Télécharger endame.eso

Retour à la liste

Numérotation des lignes : 

endame

  1. C ENDAME    SOURCE    CHAT      05/01/12    23:38:48     5004
  2.       SUBROUTINE ENDAME(IC,BETINSA)
  3. C
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6.       SEGMENT BETINSA
  7.        REAL*8 RT,RC,YOUN,XNU,GFT,GFC,HC
  8.        REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
  9.        INTEGER IFIS,IPLA,IBB,IGAU
  10.       ENDSEGMENT
  11. C
  12. *--------------------------------------------------------------------*
  13. *                                                                    *
  14. *    CALCUL DE L'ENDOMMAGEMENT                                       *
  15. *    _________________________                                       *
  16. *                                                                    *
  17. *    ENTREES:                                                        *
  18. *                                                                    *
  19. *        DKT       PARAMETRE D'ECROUISSAGE EN TRACTION               *
  20. *        DKC       PARAMETRE D'ECROUISSAGE EN COMPRESSION            *
  21. *        IC        =1 : COURBE de TRACTION                           *
  22. *                  =2 : COURBE de COMPRESSION                        *
  23. *                                                                    *
  24. *    SORTIES:                                                        *
  25. *                                                                    *
  26. *        ENDC      ENDOMMAGEMENT DE COMPRESSION                      *
  27. *        ENDT      ENDOMMAGEMENT DE TRACTION                         *
  28. *                                                                    *
  29. *--------------------------------------------------------------------*
  30. C
  31. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  32. C                      Initialisation
  33. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  34. C
  35.         BB=3.D0
  36.         DT = 0.5D0
  37.         DC = 0.25D0
  38.         ENDL = 0.99D0
  39. C
  40. C
  41. C
  42. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  43. C                      Direction de traction
  44. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  45. C
  46.            IF (IC.EQ.1) THEN
  47.         FT=RT
  48.              AT= -0.5D0
  49.              BT=FT*(HC/GFT)*(1.D0+AT/2.D0)
  50.              CT=BT*LOG(1.D0-DT)/(LOG(((1+AT)-SQRT(1+AT*AT))/(2.D0*AT)))
  51. C
  52.         ENDT=1- EXP(-CT*DKT)
  53. C
  54.         IF (ENDT.GE.ENDL) THEN
  55.          IFIS=2
  56.         ENDIF
  57.         IF (ENDT.GT.0.D0) THEN
  58.          IFIS=1
  59.         ENDIF
  60.            ENDIF
  61. C
  62. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  63. C                      Direction de compression
  64. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  65. C
  66.         IF (IC.EQ.2) THEN
  67.          FC=RC
  68.         AC=(2.D0*BB-1.D0)+2.D0* SQRT(BB**2-BB)
  69.         BC=FC*(HC/GFC)*(1.D0+AC/2.D0)
  70.         CC=BC*LOG(1.D0-DC)/(LOG((1+AC)/(2.D0*AC)))
  71. C
  72.          ENDC=1- EXP(-CC*DKC)
  73. C
  74.          IF (ENDC.GE.ENDL) THEN
  75.           IPLA=3
  76.          ENDIF
  77.          IF (ENDC.GE.DC) THEN
  78.           IPLA=2
  79.          ENDIF
  80.          IF (ENDC.GT.0.D0) THEN
  81.           IPLA=1
  82.          ENDIF
  83.         ENDIF
  84. C
  85. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  86. C                    CAS OU ENDO=ENDO_LIMITE
  87. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  88. C
  89.         IF(ENDT.GE.ENDL) ENDT= ENDL
  90.         IF(ENDC.GE.ENDL) ENDC= ENDL
  91. C
  92.             IF ((ENDT.LT.0.D0).OR.(ENDC.LT.0.D0)) THEN
  93.               WRITE(*,*)'PROBLEME DANS DAMAGE'
  94.               WRITE(*,*)'ENDT=',ENDT,'ENDC=',ENDC
  95. C             STOP
  96.             ENDIF
  97. C
  98. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  99. C
  100.         RETURN
  101.         END
  102.  
  103.  
  104.  

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