Télécharger djonl2.eso

Retour à la liste

Numérotation des lignes :

  1. C DJONL2 SOURCE CHAT 05/01/12 22:52:31 5004
  2. SUBROUTINE DJONL2(SIGI,DEPST,VARI,XMAT,SIGF,VARF,DEFP,KERRE)
  3. C-----------------------------------------------------------------------
  4. C
  5. C ROUTINE DE CALCUL DE L'ECOULEMENT DANS LE JOINT
  6. C
  7. C "DILATANT JOINT WITH COULOMBS FRICTION"
  8. C MODELE PROPOSED BY MF SNYMANN, WW BIRD AND JB MARTIN
  9. C
  10. C WORK ONLY IN 2-D PLANE (STRAIN) STATE
  11. C
  12. C INPUT
  13. C SIGI INITIAL STRESS
  14. C DEPST TRIAL INCREMENT OF ELASTIC STRAIN
  15. C XMAT MATERIAL PROPERTIES
  16. C VARI INITIAL INTERNAL VARIABLES
  17. C
  18. C OUTPUT
  19. C SIGF FINAL STRESS
  20. C VARF FINAL INTERNAL VARIABLES
  21. C DEFP INCREMENT OF PLASTIQUE STRAIN
  22. C
  23. C-----------------------------------------------------------------------
  24. C P.PEGON OCTOBER 93
  25. C-----------------------------------------------------------------------
  26. IMPLICIT INTEGER(I-N)
  27. IMPLICIT REAL*8(A-H,O-Z)
  28. -INC CCOPTIO
  29. -INC CCREEL
  30. *
  31. DIMENSION SIGI(*),DEPST(*),VARI(*),XMAT(*),
  32. . SIGF(*),VARF(*) ,DEFP(*)
  33. REAL*8 KS,KN,MU
  34. C
  35. C QUICK QUIT IN CASE OF DIMENSION ERROR
  36. C
  37. IF (IFOUR.NE.-3.AND.IFOUR.NE.-2.AND.IFOUR.NE.-1)THEN
  38. KERRE=99
  39. RETURN
  40. ELSE
  41. KERRE=0
  42. ENDIF
  43. C
  44. C-----------------------------------------------------------------------
  45. C MATERIAL PARAMETERS
  46. C ===================
  47. C
  48. C KS = XMAT(1) Shear modulus
  49. C KN = XMAT(2) Traction modulus
  50. C PHI = XMAT(5) Friction angle
  51. C MU = XMAT(6) Dilatation angle
  52. C FRTC = XMAT(7) Traction maximum (tan(PHI)*FRTC = cohesion)
  53. C-----------------------------------------------------------------------
  54. C STRESS
  55. C ======
  56. C
  57. C SMSN. = SIG.(1)
  58. C SMN. = SIG.(2)
  59. C-----------------------------------------------------------------------
  60. C STRAIN
  61. C ======
  62. C
  63. C DRSN = EPS(1)
  64. C DRN = EPS(2)
  65. C-----------------------------------------------------------------------
  66. C INTERNAL VARIABLES
  67. C ==================
  68. C
  69. C DS0. = VAR.(1) POSITION OF THE VERTEX (RELATIVE SHEAR)
  70. C DN0. = VAR.(2) POSITION OF THE VERTEX (EXTENSION)
  71. C DRSN. = VAR.(3) CURRENT DEFORMATION (RELATIVE SHEAR)
  72. C DRN. = VAR.(4) CURRENT DEFORMATION (EXTENSION)
  73. C DPSN. = VAR.(5) CURRENT PLASTIC DEFORMATION (RELATIVE SHEAR)
  74. C DPN. = VAR.(6) CURRENT PLASTIC DEFORMATION (EXTENSION)
  75. C
  76. C HC11 = VAR.(7)
  77. C HC12 = VAR.(8) CURRENT VALUES OF THE
  78. C HC21 = VAR.(9) CONSISTENT HOOK MATRIX
  79. C HC22 = VAR.(10)
  80. C-----------------------------------------------------------------------
  81. KS = XMAT(1)
  82. KN = XMAT(2)
  83. PSI = ATAN(KN/KS*TAN(XMAT(5)*XPI/180))
  84. MU = XMAT(6)*XPI/180
  85. FTRAC = XMAT(7)
  86. DTRAC = XMAT(7)/KN
  87. C
  88. TANPSI = TAN(PSI)
  89. TANMU = TAN(MU)
  90. C
  91. C TOTAL "STRAIN" FROM THE STRESS INCREMENT
  92. C
  93. DS0I = VARI(1)
  94. DN0I = VARI(2)
  95. DRSNF = VARI(3) + DEPST(1)
  96. DRNF = VARI(4) + DEPST(2)
  97. VARF(3)= DRSNF
  98. VARF(4)= DRNF
  99. C
  100. C FIRST REGION
  101. C
  102. IF ((DRNF-DTRAC).GE.(DN0I+ABS(DRSNF-DS0I)*TANMU))THEN
  103. VARF(5) = DRSNF
  104. VARF(6) = DRNF
  105. VARF(1) = DRSNF
  106. VARF(2) = DN0I+ABS(DRSNF-DS0I)*TANMU
  107. SIGF(1) = 0.D0
  108. SIGF(2) = FTRAC
  109. C
  110. VARF(7) = 0.D0
  111. VARF(8) = 0.D0
  112. VARF(9) = 0.D0
  113. VARF(10)= 0.D0
  114. C
  115. C SECOND REGION (ELASTIC CASE)
  116. C
  117. ELSEIF((DRNF-DTRAC).LE.(DN0I-ABS(DRSNF-DS0I)/TANPSI))THEN
  118. VARF(5) = DS0I
  119. VARF(6) = DN0I
  120. VARF(1) = DS0I
  121. VARF(2) = DN0I
  122. SIGF(1) = KS * (DRSNF-DS0I)
  123. SIGF(2) = KN * (DRNF -DN0I)
  124. C
  125. VARF(7) = KS
  126. VARF(8) = 0.D0
  127. VARF(9) = 0.D0
  128. VARF(10)= KN
  129. C
  130. C THIRD REGION
  131. C
  132. ELSE
  133. A=( (DRNF-DTRAC-DN0I)*SIN(PSI)+ABS(DRSNF-DS0I)*COS(PSI))
  134. > /COS(PSI-MU)
  135. B=(-(DRNF-DTRAC-DN0I)*COS(MU )+ABS(DRSNF-DS0I)*SIN(MU ))
  136. > /COS(PSI-MU)
  137. VARF(5) = DS0I + A * COS(MU) * SIGN(1.D0,DRSNF-DS0I)
  138. VARF(6) = DN0I + A * SIN(MU)
  139. VARF(1) = VARF(5)
  140. VARF(2) = VARF(6)
  141. SIGF(1) = KS * B * SIN(PSI) * SIGN(1.D0,DRSNF-DS0I)
  142. SIGF(2) = KN * B * COS(PSI) * (-1) + FTRAC
  143. C
  144. COHOOK = KN/(1+TANPSI*TANMU)
  145. TANPHI = TANPSI*KS/KN
  146. VARF(7) = COHOOK*TANPHI*TANMU
  147. VARF(8) =-COHOOK*SIGN(1.D0,DRSNF-DS0I)*TANPHI
  148. VARF(9) =-COHOOK*SIGN(1.D0,DRSNF-DS0I)*TANMU
  149. VARF(10)= COHOOK
  150. C
  151. ENDIF
  152. C
  153. C INCREMENT OF PLASTIC STRAIN
  154. C
  155. DEFP(1) = VARF(5) - VARI(5)
  156. DEFP(2) = VARF(6) - VARI(6)
  157. C
  158. RETURN
  159. END
  160.  
  161.  
  162.  
  163.  
  164.  

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