Télécharger forecr.eso

Retour à la liste

Numérotation des lignes :

forecr
  1. C FORECR SOURCE CHAT 05/01/13 00:06:35 5004
  2. SUBROUTINE FORECR(DK,PAEC,IC,SEQ,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. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  13. * *
  14. * CALCUL DU PARAMETRE D'ECROUISSAGE *
  15. * _________________________________ *
  16. * *
  17. * ENTREES: *
  18. * *
  19. * DK VARIABLE D'ECROUISSAGE *
  20. * IC =1 : COURBE de TRACTION *
  21. * =2 : COURBE de COMPRESSION *
  22. * *
  23. * SORTIES: *
  24. * *
  25. * SEQ CONTRAINTE EQUIVALENTE *
  26. * PAEC PENTE *
  27. * *
  28. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  29. C
  30. BB=3.D0
  31. DT = 0.5
  32. DC = 0.25
  33. ENDL = 0.99
  34. C
  35. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  36. C Direction de traction
  37. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  38. C
  39. IF (IC.EQ.1) THEN
  40. F=RT
  41. A= -0.5D0
  42. B=F*(HC/GFT)*(1.D0+A/2.D0)
  43. C=B*LOG(1.D0-DT)/(LOG(((1+A)-SQRT(1+A*A))/(2.D0*A)))
  44. C
  45. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  46. C Direction de compression
  47. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  48. C
  49. ELSE
  50. F=RC/BB
  51. A=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
  52. B=F*(HC/GFC)*(1.D0+A/2.D0)
  53. C=B*LOG(1.D0-DC)/(LOG((1+A)/(2.D0*A)))
  54. ENDIF
  55. C
  56. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  57. C CALCUL LA CONTRAINTE UNIAXIALE ET DE LA PENTE
  58. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  59. C 1: CONTRAINTE UNIAXIALE
  60. C
  61. EXP1=EXP(-B*DK)
  62. DD=1.D0-C/B
  63. SEQ=F*((1.D0+A)*EXP1**(DD)-A*EXP1**(1.D0+DD))
  64. C
  65. C 2: CALCUL DE LA (+)*PENTE
  66. C
  67. PAEC=-B*F*(DD*(1.D0+A)*EXP1**(DD)-
  68. & (1.D0+DD)*A*EXP1**(1.D0+DD))
  69. C
  70. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  71. C
  72. IF(ABS(SEQ).LT.1.D-8)THEN
  73. SEQ=1.D-8
  74. PAEC=0.D0
  75. GOTO 10
  76. ENDIF
  77. C
  78. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  79. C
  80. 10 CONTINUE
  81. C
  82. RETURN
  83. END
  84.  
  85.  
  86.  

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