Esope source files

   1 : C TUFK1A    SOURCE    CHAT      05/01/13    03:56:36     5004           
   2 :       SUBROUTINE TUFK1A(XM,XP,YOU,RAYMO,EPAI,TETA,XK1,AIRE)
   3 : C=======================================================================
   4 : C   CALCULE LE FACTEUR D INTENSITE DES CONTRAINTES
   5 : C                ET L AIRE DE BRECHE
   6 : C    ENTREES
   7 : C     XM    = MOMENT
   8 : C     XP    = EFFORT NORMAL
   9 : C     YOU   = MODULE D YOUNG
  10 : C     RAYMO = RAYON MOYEN
  11 : C     EPAI  = EPAISSEUR
  12 : C     TETA  = 1/2 ANGLE DE LA FISSURE EN RADIANS
  13 : C    SORTIES
  14 : C     XK1   = FACTEUR D INTENSITE DES CONTRAINTES
  15 : C     AIRE  = AIRE DE BRECHE
  16 : C     
  17 : C    APPELE PAR TUFPLA
  18 : C=======================================================================
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)
  21 : C    Include contenant quelques constantes dont XPI :
  22 : -INC CCREEL
  23 :       PARAMETER(DEUX=2.D0)
  24 : C     
  25 : C     CALCUL DE A   COEFIICIENT ZAHOR
  26 : C     
  27 :       RSURT=RAYMO / EPAI
  28 :       IF(RSURT.LE.10.D0) THEN
  29 :          AXX = ( .125D0*RSURT - .25D0 ) **.25D0
  30 :       ELSE
  31 :          AXX = ( .4D0*RSURT - 3.D0 ) **.25D0
  32 :       ENDIF
  33 : C     
  34 :       TETA1= (TETA * 180.D0) / XPI
  35 :       TESPI = TETA/XPI
  35 :       CALL TUFIFP(TESPI,AXX,FP,FM,FMP,FOP,FOM)
  36 : C     
  37 : C     FACTEUR D INTENSITE DES CONTRAINTES
  38 : C     
  39 :       IF(TETA1.LE.(0.5D0))THEN
  40 :          XK1=XZERO
  41 :          GOTO 666
  42 :       ENDIF
  43 :       SQQ= XPI * RAYMO * TETA
  44 :       SQQ= SQRT(SQQ)
  45 :       XEX= SQQ * FOP/(DEUX * XPI * RAYMO *EPAI)
  46 :       XFL= SQQ * FOM/(XPI * RAYMO * RAYMO *EPAI)
  47 :       XK1=XEX * XP + XFL * XM
  48 : C     
  49 : C     CALCUL DES AIRES DE BRECHE NOTE TECHNIQUE DRE/STRE/LMA 85/695
  50 : C     
  51 :       SIGM = XP /( DEUX * XPI * RAYMO * EPAI )
  52 :       SIGF = XM /( XPI * RAYMO * RAYMO * EPAI )
  53 :       XIM =  XPI * RAYMO * RAYMO * DEUX * TETA * TETA * FP /YOU
  54 :       XIF =  XIM * ( .75D0 +(.25D0 * COS ( TETA )))
  55 :       AIRE=XIM * SIGM  + XIF * SIGF
  56 : C     
  57 :  666  CONTINUE
  58 :       RETURN
  59 :       END
  60 : 
  61 : 
  62 : 
  63 :  
  64 :  
  65 :