Télécharger la plaquette Cast3M


fr en

Télécharger forecr.eso

Retour à la liste

Numérotation des lignes : 

forecr

  1. C FORECR    SOURCE    CHAT      05/01/13    00:06:35     5004
  2.       SUBROUTINE FORECR(DK,PAEC,IC,SEQ,BETINSA)
  3. C
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6.       SEGMENT BETINSA
  7.        REAL*8 RT,RC,YOUN,XNU,GFT,GFC,HC
  8.        REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
  9.        INTEGER IFIS,IPLA,IBB,IGAU
  10.       ENDSEGMENT
  11. C
  12. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  13. *                                                                    *
  14. *    CALCUL DU PARAMETRE D'ECROUISSAGE                               *
  15. *    _________________________________                               *
  16. *                                                                    *
  17. *    ENTREES:                                                        *
  18. *                                                                    *
  19. *        DK        VARIABLE D'ECROUISSAGE                            *
  20. *        IC        =1 : COURBE de TRACTION                           *
  21. *                  =2 : COURBE de COMPRESSION                        *
  22. *                                                                    *
  23. *    SORTIES:                                                        *
  24. *                                                                    *
  25. *        SEQ       CONTRAINTE EQUIVALENTE                            *
  26. *        PAEC      PENTE                                             *
  27. *                                                                    *
  28. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  29. C
  30.         BB=3.D0
  31.         DT = 0.5
  32.         DC = 0.25
  33.         ENDL = 0.99
  34. C
  35. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  36. C                      Direction de traction
  37. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  38. C
  39.            IF (IC.EQ.1) THEN
  40.        F=RT
  41.             A= -0.5D0
  42.             B=F*(HC/GFT)*(1.D0+A/2.D0)
  43.             C=B*LOG(1.D0-DT)/(LOG(((1+A)-SQRT(1+A*A))/(2.D0*A)))
  44. C
  45. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  46. C                      Direction de compression
  47. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  48. C
  49.            ELSE
  50.        F=RC/BB
  51.             A=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
  52.             B=F*(HC/GFC)*(1.D0+A/2.D0)
  53.        C=B*LOG(1.D0-DC)/(LOG((1+A)/(2.D0*A)))
  54.            ENDIF
  55. C
  56. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  57. C       CALCUL LA CONTRAINTE UNIAXIALE ET DE LA PENTE
  58. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  59. C       1: CONTRAINTE UNIAXIALE
  60. C
  61.             EXP1=EXP(-B*DK)
  62.             DD=1.D0-C/B
  63.             SEQ=F*((1.D0+A)*EXP1**(DD)-A*EXP1**(1.D0+DD))
  64. C
  65. C       2: CALCUL DE LA (+)*PENTE
  66. C
  67.             PAEC=-B*F*(DD*(1.D0+A)*EXP1**(DD)-
  68.      &                 (1.D0+DD)*A*EXP1**(1.D0+DD))
  69. C
  70. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  71. C
  72.            IF(ABS(SEQ).LT.1.D-8)THEN
  73.             SEQ=1.D-8
  74.             PAEC=0.D0
  75.             GOTO 10
  76.            ENDIF
  77. C
  78. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  79. C
  80.    10   CONTINUE
  81. C
  82.         RETURN
  83.         END
  84.  
  85.  
  86.  

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