Télécharger unico1.eso

Retour à la liste

Numérotation des lignes :

  1. C UNICO1 SOURCE CHAT 05/01/13 03:58:30 5004
  2. C UNICO1 SOURCE NECH 99/05/26
  3. SUBROUTINE UNICO1(DK,PAEC,IC,SEQ,BETJEF,NECH0,NECH1)
  4. C
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. * COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,
  8. * & TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  9. * & ITER,ISIM,IBB,IGAU1,IZON
  10. * COMMON /DBETNECH0/DT,DC,ALFG,S0
  11. * COMMON /DBETNECH1/ENDL
  12. SEGMENT BETJEF
  13. REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,TDEF,
  14. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  15. INTEGER ICT,ICC,IMOD,IVISS,ITER,
  16. & ISIM,IBB,IGAU1,IZON
  17. ENDSEGMENT
  18. SEGMENT NECH0
  19. REAL*8 DT,DC,ALFG,S0,ENDO
  20. ENDSEGMENT
  21. SEGMENT NECH1
  22. REAL*8 ENDL
  23. ENDSEGMENT
  24. C
  25. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  26. * *
  27. * CALCUL DU PARAMETRE D'ECROUISSAGE *
  28. * _________________________________ *
  29. * *
  30. * ENTREES: *
  31. * *
  32. * DK VARIABLE D'ECROUISSAGE *
  33. * IC =1 : COURBE de TRACTION *
  34. * =2 : COURBE de COMPRESSION *
  35. * *
  36. * SORTIES: *
  37. * *
  38. * SEQ CONTRAINTE EQUIVALENTE *
  39. * PAEC PENTE *
  40. * *
  41. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  42. C
  43. BB=1.D0/AA
  44. C
  45. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  46. C Direction de traction
  47. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  48. C
  49. IF (IC.EQ.1) THEN
  50. F=ALFA*RB
  51. A= -0.5D0
  52. B=F*(HC/GFT)*(1.D0+A/2.D0)
  53. C=B*LOG(1.D0-DT)/(LOG(((1+A)-SQRT(1+A*A))/(2.D0*A)))
  54. C
  55. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  56. C Direction de compression
  57. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  58. C
  59. ELSE
  60. F=AA*RB
  61. A=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
  62. B=F*(HC/GFC)*(1.D0+A/2.D0)
  63. C=B*LOG(1.D0-DC)/(LOG((1+A)/(2.D0*A)))
  64. ENDIF
  65. C
  66. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  67. C CALCUL DE LA DEFORMATION PLASTIQUE LIMITE
  68. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  69. C
  70. DKL=LOG(1.D0-ENDL)/(-C)
  71. C
  72. IF(DK.GE.DKL)THEN
  73. DK=DKL
  74. ENDIF
  75. C
  76. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  77. C CALCUL LA CONTRAINTE UNIAXIALE ET DE LA PENTE
  78. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  79. C 1: CONTRAINTE UNIAXIALE
  80. C
  81. EXP1=EXP(-B*DK)
  82. DD=1.D0-C/B
  83. SEQ=F*((1.D0+A)*EXP1**(DD)-A*EXP1**(1.D0+DD))
  84. C
  85. C 2: CALCUL DE LA (+)*PENTE
  86. C
  87. PAEC=-B*F*(DD*(1.D0+A)*EXP1**(DD)-
  88. & (1.D0+DD)*A*EXP1**(1.D0+DD))
  89. C
  90. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  91. C
  92. IF(ABS(SEQ).LT.1.D-8)THEN
  93. SEQ=1.D-8
  94. PAEC=0.D0
  95. GOTO 10
  96. ENDIF
  97. C
  98. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  99. C
  100. 10 CONTINUE
  101. C
  102. RETURN
  103. END
  104.  
  105.  
  106.  

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