Télécharger la plaquette Cast3M


fr en

Télécharger precn1.eso

Retour à la liste

Numérotation des lignes : 

precn1

  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.  
  34. -INC PPARAM
  35. -INC CCOPTIO
  36.  
  37.  
  38.       icomp=0
  39.       valstr = ps1
  40.  
  41.       SP = F1*FAI + F2*SLON
  42.  
  43. *    ... perte par frottement ...
  44.  
  45.          valstr=PS1*EXP(-SP)
  46.  
  47. *    ... perte éventuelle par glissement à l'ancrage ...
  48.  
  49.         IF(XLMBDA.GT.1.e-20) THEN
  50. *
  51.           IF(LAST.EQ.0) THEN
  52. *          cas ou la longueur d'influence est inferieure
  53. *          a la longueur du cable
  54.  
  55.              IF  (SLON.LT.XLMBDA) THEN
  56.                valstr = valstr * EXP (2.D0*(SP-SD))
  57.              ENDIF
  58.           ELSE
  59. *          cas ou la longueur d'influence est superieure
  60. *          a la longueur du cable
  61.  
  62.              VALSTR=PSL * EXP(SP-SD)
  63.           ENDIF
  64.         ENDIF
  65. *
  66.  
  67.       IF ( valstr.LT.(0.D0) ) THEN
  68.          ICOMP = 1
  69.          RETURN
  70.       ENDIF
  71.  
  72. ***** IF(FPRG.EQ.0.)FPRG=1.D-30
  73.  
  74.       IF(FPRG.NE.0.D0) THEN
  75.         VAL = (0.06D0*((valstr/(FPRG*SECT1) )-RMU0)*valstr*RH10)
  76.         IF ( VAL.GT.(0.) ) THEN
  77.            valstr = valstr - VAL
  78.         ENDIF
  79.       ENDIF
  80.  
  81.       if ( valstr.lt.sini) valstr = sini
  82.  
  83.       return
  84.       END
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  

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