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