Télécharger precn1.eso

Retour à la liste

Numérotation des lignes :

  1. C PRECN1 SOURCE CHAT 05/10/27 21:17:35 5227
  2. SUBROUTINE PRECN1(PS1,EA,F1,F2,GANC,RMU0,FPRG,RH10,SLON,
  3. & FAI,SINI,SECT1,XLMBDA,SD,LAST,PSL,valstr,ICOMP)
  4. C======================================================================C
  5. C C
  6. C ENTREES : C
  7. C PS1 = EFFORT DE TENSION DU CABLE EN BOUT C
  8. C EA = MODULE d'YOUNG DE l'ACIER C
  9. C F1 = COEFFICIENT DE FROTTEMENT C
  10. C F2 = COEFFICIENT DE FROTTEMENT PAR METRE LINEAIRE DROIT C
  11. C GANC = GLISSEMENT A L'ANCRAGE C
  12. C RMU0 = COEFFICIENT EGAL A 0.43 POUR LES ACIERS TBR C
  13. C FPRG = CONTRAINTE DE RUPTURE GARANTIE C
  14. C RH10 = RELAXATION DES ACIERS A 1000 HEURES EN % C
  15. C SLON = LONGUEUR DU CABLE A PARTIR DE L'EXTREMITE C
  16. C FAI = VARIATION D'ANGLE TOTAL DE L ANCRAGE A L ELEMENT C
  17. C SINI = CONTRAINTE INITIALE AU POINT CONSIDERE C
  18. C sect1 = section du cable C
  19. C XLMBDA = DISTANCE D'INFLUENCE DU RECUL D'ANCRAGE
  20. C SD = VALEUR DE L'EXPOSANT POUR X=XLMBDA
  21. C LAST = 1 SI XLMDA DEPASSE LA LONGUEUR DU CABLE ( 0 SINON)
  22. C PSL = CONTRAINTE EN BOUT SI LAST = 1
  23. C SORTIES : C
  24. C valstr = EFFORT AU POINT DE GAUSS DE LELEMENT COURANT CABLE C
  25. C ICOMP = REPERE LA COMPATIBILITE DES DONNEES (=1 SI PB) C
  26. C C
  27. C----------------------------------------------------------------------C
  28. C Modif JMB NOV 97
  29. C======================================================================C
  30. IMPLICIT INTEGER(I-N)
  31. IMPLICIT REAL*8(A-H,O-Z)
  32. C
  33. -INC CCOPTIO
  34.  
  35.  
  36. icomp=0
  37. valstr = ps1
  38.  
  39. SP = F1*FAI + F2*SLON
  40.  
  41. * ... perte par frottement ...
  42.  
  43. valstr=PS1*EXP(-SP)
  44.  
  45. * ... perte éventuelle par glissement à l'ancrage ...
  46.  
  47. IF(XLMBDA.GT.1.e-20) THEN
  48. *
  49. IF(LAST.EQ.0) THEN
  50. * cas ou la longueur d'influence est inferieure
  51. * a la longueur du cable
  52.  
  53. IF (SLON.LT.XLMBDA) THEN
  54. valstr = valstr * EXP (2.D0*(SP-SD))
  55. ENDIF
  56. ELSE
  57. * cas ou la longueur d'influence est superieure
  58. * a la longueur du cable
  59.  
  60. VALSTR=PSL * EXP(SP-SD)
  61. ENDIF
  62. ENDIF
  63. *
  64.  
  65. IF ( valstr.LT.(0.D0) ) THEN
  66. ICOMP = 1
  67. RETURN
  68. ENDIF
  69.  
  70. ***** IF(FPRG.EQ.0.)FPRG=1.D-30
  71.  
  72. IF(FPRG.NE.0.D0) THEN
  73. VAL = (0.06D0*((valstr/(FPRG*SECT1) )-RMU0)*valstr*RH10)
  74. IF ( VAL.GT.(0.) ) THEN
  75. valstr = valstr - VAL
  76. ENDIF
  77. ENDIF
  78.  
  79. if ( valstr.lt.sini) valstr = sini
  80.  
  81. return
  82. END
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  

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