Télécharger la plaquette Cast3M


fr en

Télécharger damag1.eso

Retour à la liste

Numérotation des lignes : 

damag1

  1. C DAMAG1    SOURCE    CHAT      05/01/12    22:36:09     5004
  2.       SUBROUTINE DAMAG1(DKPT,DKPC,ENDC,ENDT,IFISU,IPLA,IC
  3.      &  ,BETJEF,NECH0,NECH1)
  4. C
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7. C
  8. *      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,
  9. *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  10. *     &                ITER,ISIM,IBB,IGAU1,IZON
  11. *      COMMON /DBETNECH0/DT,DC,ALFG,S0
  12. *      COMMON /DBETNECH1/ENDL
  13.       SEGMENT BETJEF
  14.        REAL*8 AA,BETA,RB,ALFA,EX,XNU,GFC,GFT,HC,ETA,TDEF,
  15.      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  16.        INTEGER ICT,ICC,IMOD,IVIS,ITER,
  17.      &                ISIM,IBB,IGAU,IZON
  18.       ENDSEGMENT
  19.       SEGMENT NECH0
  20.        REAL*8 DT,DC,ALFG,S0,ENDO
  21.       ENDSEGMENT
  22.       SEGMENT NECH1
  23.        REAL*8 ENDL
  24.       ENDSEGMENT
  25. C
  26. *--------------------------------------------------------------------*
  27. *                                                                    *
  28. *    CALCUL DE L'ENDOMMAGEMENT                                       *
  29. *    _________________________                                       *
  30. *                                                                    *
  31. *    ENTREES:                                                        *
  32. *                                                                    *
  33. *        DKT       PARAMETRE D'ECROUISSAGE EN TRACTION               *
  34. *        DKC       PARAMETRE D'ECROUISSAGE EN COMPRESSION            *
  35. *        IC        =1 : COURBE de TRACTION                           *
  36. *                  =2 : COURBE de COMPRESSION                        *
  37. *                       =3 : COURBE TRACTION COMPRESSION                  *
  38. *                                                                    *
  39. *    SORTIES:                                                        *
  40. *                                                                    *
  41. *        ENDC      ENDOMMAGEMENT DE COMPRESSION                      *
  42. *        ENDT      ENDOMMAGEMENT DE TRACTION                         *
  43. *                                                                    *
  44. *--------------------------------------------------------------------*
  45. C
  46. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  47. C                      Initialisation
  48. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  49. C
  50.       BB=1.D0/AA
  51. C
  52. C
  53. C
  54. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  55. C                      Direction de traction
  56. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  57. C
  58.            IF (IC.EQ.1.OR.IC.EQ.3) THEN
  59.          FT=ALFA*RB
  60.               AT= -0.5D0
  61.               BT=FT*(HC/GFT)*(1.D0+AT/2.D0)
  62.              CT=BT*LOG(1.D0-DT)/(LOG(((1+AT)-SQRT(1+AT*AT))/(2.D0*AT)))
  63. C
  64.         ENDT=1- EXP(-CT*DKPT)
  65.            ENDIF
  66. C
  67. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  68. C                      Direction de compression
  69. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  70. C
  71.         IF (IC.EQ.2.OR.IC.EQ.3) THEN
  72.          FC=AA*RB
  73.               AC=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
  74.               BC=FC*(HC/GFC)*(1.D0+AC/2.D0)
  75.               CC=BC*LOG(1.D0-DC)/(LOG((1+AC)/(2.D0*AC)))
  76. C
  77.          ENDC=1- EXP(-CC*DKPC)
  78.         ENDIF
  79. C
  80. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  81. C                    CAS OU ENDO=ENDO_LIMITE
  82. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  83. C
  84.         IF(ENDT.GE.ENDL) ENDT= ENDL
  85.         IF(ENDC.GE.ENDL) ENDC= ENDL
  86. C
  87.             IF ((ENDT.LT.0.D0).OR.(ENDC.LT.0.D0)) THEN
  88.               WRITE(*,*)'PROBLEME DANS DAMAGE'
  89.               WRITE(*,*)'ENDT=',ENDT,'ENDC=',ENDC
  90. C             STOP
  91.             ENDIF
  92. C
  93. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  94. C
  95.         RETURN
  96.         END
  97.  
  98.  
  99.  

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