Télécharger la plaquette Cast3M


fr en

Télécharger unico1.eso

Retour à la liste

Numérotation des lignes : 

unico1

  1. C UNICO1    SOURCE    CHAT      05/01/13    03:58:30     5004
  2. C     UNICO1    SOURCE   NECH   99/05/26
  3.       SUBROUTINE UNICO1(DK,PAEC,IC,SEQ,BETJEF,NECH0,NECH1)
  4. C
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7. *      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,
  8. *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  9. *     &                ITER,ISIM,IBB,IGAU1,IZON
  10. *      COMMON /DBETNECH0/DT,DC,ALFG,S0
  11. *      COMMON /DBETNECH1/ENDL
  12.       SEGMENT BETJEF
  13.        REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,HC,ETA,TDEF,
  14.      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  15.        INTEGER ICT,ICC,IMOD,IVISS,ITER,
  16.      &                ISIM,IBB,IGAU1,IZON
  17.       ENDSEGMENT
  18.       SEGMENT NECH0
  19.        REAL*8 DT,DC,ALFG,S0,ENDO
  20.       ENDSEGMENT
  21.       SEGMENT NECH1
  22.        REAL*8 ENDL
  23.       ENDSEGMENT
  24. C
  25. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  26. *                                                                    *
  27. *    CALCUL DU PARAMETRE D'ECROUISSAGE                               *
  28. *    _________________________________                               *
  29. *                                                                    *
  30. *    ENTREES:                                                        *
  31. *                                                                    *
  32. *        DK        VARIABLE D'ECROUISSAGE                            *
  33. *        IC        =1 : COURBE de TRACTION                           *
  34. *                  =2 : COURBE de COMPRESSION                        *
  35. *                                                                    *
  36. *    SORTIES:                                                        *
  37. *                                                                    *
  38. *        SEQ       CONTRAINTE EQUIVALENTE                            *
  39. *        PAEC      PENTE                                             *
  40. *                                                                    *
  41. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*
  42. C
  43.         BB=1.D0/AA
  44. C
  45. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  46. C                      Direction de traction
  47. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  48. C
  49.            IF (IC.EQ.1) THEN
  50.        F=ALFA*RB
  51.             A= -0.5D0
  52.             B=F*(HC/GFT)*(1.D0+A/2.D0)
  53.             C=B*LOG(1.D0-DT)/(LOG(((1+A)-SQRT(1+A*A))/(2.D0*A)))
  54. C
  55. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  56. C                      Direction de compression
  57. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  58. C
  59.            ELSE
  60.        F=AA*RB
  61.             A=(2.D0*BB-1.D0)+2.D0*SQRT(BB**2-BB)
  62.             B=F*(HC/GFC)*(1.D0+A/2.D0)
  63.             C=B*LOG(1.D0-DC)/(LOG((1+A)/(2.D0*A)))
  64.            ENDIF
  65. C
  66. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  67. C       CALCUL DE LA DEFORMATION PLASTIQUE LIMITE
  68. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  69. C
  70.         DKL=LOG(1.D0-ENDL)/(-C)
  71. C
  72.            IF(DK.GE.DKL)THEN
  73.             DK=DKL
  74.            ENDIF
  75. C
  76. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  77. C       CALCUL LA CONTRAINTE UNIAXIALE ET DE LA PENTE
  78. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  79. C       1: CONTRAINTE UNIAXIALE
  80. C
  81.             EXP1=EXP(-B*DK)
  82.             DD=1.D0-C/B
  83.             SEQ=F*((1.D0+A)*EXP1**(DD)-A*EXP1**(1.D0+DD))
  84. C
  85. C       2: CALCUL DE LA (+)*PENTE
  86. C
  87.             PAEC=-B*F*(DD*(1.D0+A)*EXP1**(DD)-
  88.      &                  (1.D0+DD)*A*EXP1**(1.D0+DD))
  89. C
  90. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  91. C
  92.            IF(ABS(SEQ).LT.1.D-8)THEN
  93.             SEQ=1.D-8
  94.             PAEC=0.D0
  95.             GOTO 10
  96.            ENDIF
  97. C
  98. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  99. C
  100.    10   CONTINUE
  101. C
  102.         RETURN
  103.         END
  104.  
  105.  
  106.  

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