Télécharger la plaquette Cast3M


fr en

Télécharger tufk1a.eso

Retour à la liste

Numérotation des lignes : 

tufk1a

  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
  36.       CALL TUFIFP(TESPI,AXX,FP,FM,FMP,FOP,FOM)
  37. C
  38. C     FACTEUR D INTENSITE DES CONTRAINTES
  39. C
  40.       IF(TETA1.LE.(0.5D0))THEN
  41.          XK1=XZERO
  42.          GOTO 666
  43.       ENDIF
  44.       SQQ= XPI * RAYMO * TETA
  45.       SQQ= SQRT(SQQ)
  46.       XEX= SQQ * FOP/(DEUX * XPI * RAYMO *EPAI)
  47.       XFL= SQQ * FOM/(XPI * RAYMO * RAYMO *EPAI)
  48.       XK1=XEX * XP + XFL * XM
  49. C
  50. C     CALCUL DES AIRES DE BRECHE NOTE TECHNIQUE DRE/STRE/LMA 85/695
  51. C
  52.       SIGM = XP /( DEUX * XPI * RAYMO * EPAI )
  53.       SIGF = XM /( XPI * RAYMO * RAYMO * EPAI )
  54.       XIM =  XPI * RAYMO * RAYMO * DEUX * TETA * TETA * FP /YOU
  55.       XIF =  XIM * ( .75D0 +(.25D0 * COS ( TETA )))
  56.       AIRE=XIM * SIGM  + XIF * SIGF
  57. C
  58.  666  CONTINUE
  59.       RETURN
  60.       END
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  
  67.  

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