Esope source files

   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 :