Télécharger la plaquette Cast3M


fr en

Télécharger sjonc2.eso

Retour à la liste

Numérotation des lignes : 

sjonc2

  1. C SJONC2    SOURCE    CHAT      05/01/13    03:20:33     5004
  2.       SUBROUTINE SJONC2(SIGI,DEPST,VARI,XMAT,TRAS,NTRAS,
  3.      >                                     SIGF,VARF,DEFP,KERRE)
  4. C-----------------------------------------------------------------------
  5. C
  6. C     ROUTINE DE CALCUL DE L'ECOULEMENT DANS LE JOINT
  7. C
  8. C     WORK ONLY IN 2-D PLANE STATE
  9. C
  10. C  INPUT
  11. C     SIGI    INITIAL STRESS
  12. C     DEPST   INCREMENT OF STRAIN
  13. C     XMAT    MATERIAL PROPERTIES
  14. C     TRAS    SHEAR CURVE
  15. C     VARI    INITIAL INTERNAL VARIABLES
  16. C
  17. C  OUTPUT
  18. C     SIGF    FINAL STRESS
  19. C     VARF    FINAL INTERNAL VARIABLES
  20. C     DEFP    INCREMENT OF PLASTIQUE STRAIN
  21. C
  22. C-----------------------------------------------------------------------
  23. C  P.PEGON   SEPTEMBER 2001
  24. C-----------------------------------------------------------------------
  25.       IMPLICIT INTEGER(I-N)
  26.       IMPLICIT REAL*8(A-H,O-Z)
  27. *
  28.       DIMENSION SIGI(*),DEPST(*),VARI(*),
  29.      .          SIGF(*),VARF(*) ,DEFP(*),
  30.      .          XMAT(*),TRAS(2,NTRAS)
  31. C
  32. C-----------------------------------------------------------------------
  33. C  MATERIAL PARAMETERS
  34. C  ===================
  35. C
  36. C      KS    = XMAT(1)   Shear modulus
  37. C      TRAS  = X/Y       PLASTIC Shear curve
  38. C      BETA  = XMAT(5)   Unloading Coefficient
  39. C-----------------------------------------------------------------------
  40. C  STRESS
  41. C  ======
  42. C
  43. C      SMS1.  = SIG.(1)
  44. C      SMS2   = SIG.(2)
  45. C-----------------------------------------------------------------------
  46. C  STRAIN
  47. C  ======
  48. C
  49. C      DRS1   = EPS(1)
  50. C      DRS2   = EPS(2)
  51. C-----------------------------------------------------------------------
  52. C  INTERNAL VARIABLES
  53. C  ==================
  54. C
  55. C
  56. C      EPSS   = VAR.(1) PLASTIC SHEAR COEFFICIENT
  57. C
  58. C
  59. C-----------------------------------------------------------------------
  60. C
  61. C        CURRENT POSITION ON THE LOADING CURVE
  62. C
  63.       EPSS=VARI(1)
  64.       CALL YOFXCU(EPSS,TRAS,NTRAS,SMSNP,DSMSNP,KERRE)
  65.       IF (KERRE.NE.0) RETURN
  66. C
  67. C        VALUE OF THE LOADING/UNLOADING MODULUS
  68. C
  69.       DMODU=SMSNP/(EPSS+(TRAS(2,1)/XMAT(1)))
  70.       DMODU=(1-XMAT(5))*XMAT(1)+DMODU*XMAT(5)
  71. C
  72. C        ELASTIC STRESS PREDICTION
  73. C
  74.       SIGF(1)=SIGI(1)+DMODU*DEPST(1)
  75.       SIGF(2)=SIGI(2)+DMODU*DEPST(2)
  76.       SMSN=SQRT(SIGF(1)**2+SIGF(2)**2)
  77. C
  78. C        TEST ON THE FLOW
  79. C
  80.       IF (SMSN.GT.SMSNP) THEN
  81. C
  82. C        -> YIELDING IN SHEAR ...
  83. C
  84.         EPSS= EPSS+(SMSN-SMSNP)/DMODU
  85.         CALL YOFXCU(EPSS,TRAS,NTRAS,SMSNP,DSMSNP,KERRE)
  86.         IF (KERRE.NE.0) RETURN
  87. C
  88.         SIGF(1)=SIGF(1)*SMSNP/SMSN
  89.         SIGF(2)=SIGF(2)*SMSNP/SMSN
  90. C
  91.       ENDIF
  92. C
  93.       VARF(1)= EPSS
  94.       DEFP(1) = VARF(1) - VARI(1)
  95. C
  96.       RETURN
  97.       END
  98.  
  99.  
  100.  

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