Télécharger damag1.eso

Retour à la liste

Numérotation des lignes :

  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