Télécharger la plaquette Cast3M


fr en

Télécharger ancrel.eso

Retour à la liste

Numérotation des lignes : 

ancrel

  1. C ANCREL    SOURCE    CHAT      05/01/12    21:22:18     5004
  2. C
  3.       SUBROUTINE ANCREL(SIGI,DEPST,VARI,XMAT,
  4.      .                       SIGF,VARF,DEFP,KERRE)
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8(A-H,O-Z)
  7.       DIMENSION XMAT(*),VARI(*),VARF(*),SIGI(*)
  8.      . ,DEPST(*),DEFP(*),SIGF(*)
  9. C----------------------------------------------------------------------
  10. C
  11. C     LOI D'ANCRAGE ELIGEHAUSEN SUR UN ELEMENT DE JOINT 2D
  12. C----------------------------------------------------------------------
  13. C     Didier COMBESCURE - EMSI - MAI 2000
  14. C----------------------------------------------------------------------
  15. C
  16. C      XMAT( 1) XKS
  17. C      XMAT( 2) XKN
  18. C      XMAT( 3) RHO
  19. C      XMAT( 4) ALPHA
  20. C      XMAT( 5) S1T
  21. C      XMAT( 6) S2T
  22. C      XMAT( 7) S3T
  23. C      XMAT( 8) T1T
  24. C      XMAT( 9) T3T
  25. C      XMAT(10) ALFA
  26. C      XMAT(11) PERI
  27. C
  28. C  STRESS
  29. C  ======
  30. C
  31. C      SMSN.  = SIG.(1)
  32. C      SMN.   = SIG.(2)
  33. C-----------------------------------------------------------------------
  34. C  STRAIN
  35. C  ======
  36. C
  37. C      DRSN   = EPS(1)
  38. C      DRN    = EPS(2)
  39. C-----------------------------------------------------------------------
  40. C  INTERNAL VARIABLES
  41. C  ==================
  42. C
  43. C      EPPP = VAR0(1) VARF(1) EPPP GLISSEMENT PLASTIQUE POSITIF
  44. C      EPPM = VAR0(2) VARF(2) EPPM GLISSEMENT PLASTIQUE NEGATIF
  45. C      EPSO = VAR0(3) VARF(3) EPSO GLISSEMENT TOTAL
  46. C
  47. C  MATRICE TANGENTE
  48. C
  49. C      KSTAN   = VARF(4)
  50. C----------------------------------------------------------------------
  51. C
  52.        XKS  = XMAT(1)
  53.        XKN  = XMAT(2)
  54.        S1T = XMAT(5)
  55.        S2T = XMAT(6)
  56.        S3T = XMAT(7)
  57.        T1T = XMAT(8)
  58.        T3T = XMAT(9)
  59.        ALFA = XMAT(10)
  60.        PERI = XMAT(11)
  61. C
  62.        G12 = XKS/PERI
  63. C
  64. C  EFFORT NORMAL (ELASTIQUE)
  65. C
  66.         SIGF(2) = SIGI(2)+(XKN*DEPST(2))
  67. C
  68. C  CISAILLEMENT - GLISSEMENT
  69. C
  70.        EPPP   = VARI(1)
  71.        EPPM   = VARI(2)
  72. C
  73.        EPSO = VARI(3) + DEPST(1)
  74. C
  75.        CALL ELIGEH(EPSO,G12,S1T,S2T,S3T,T1T,T3T,ALFA,
  76.      .            EPPP,EPPM,EPPPF,EPPMF,XKTAN,TAUF)
  77. C
  78.        SIGF(1) = PERI * TAUF
  79. C
  80.        VARF(1) = EPPPF
  81.        VARF(2) = EPPMF
  82.        VARF(3) = EPSO
  83. C
  84.        VARF(4) = PERI * XKTAN
  85.        VARF(5) = XKN
  86. C
  87.        RETURN
  88.        END
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  

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