Télécharger la plaquette Cast3M


fr en

Télécharger clouve.eso

Retour à la liste

Numérotation des lignes : 

clouve

  1. C CLOUVE    SOURCE    CHAT      05/01/12    22:08:47     5004
  2.       SUBROUTINE CLOUVE (AAA,YUNG,XNU,ET,CISA,IENDO,KERRE)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5.       DIMENSION ET(*),CISA(*),AAA(6,6),UMBETA(3)
  6. C
  7. C***********************************************************************
  8. C***************** CALCUL DU MODULE DE FISSURATION *********************
  9. C***********************************************************************
  10. C
  11. C     INITIALISATION
  12. C
  13.       CALL ZDANUL(AAA,36)
  14.       ET1=ET(1)
  15.       ET2=ET(2)
  16.       ET3=ET(3)
  17.       UPUN=1.D0+XNU
  18.       UMUN=1.D0-XNU
  19.       UMDN=1.D0-2.D0*XNU
  20.       EMET1=YUNG-ET1
  21.       EMET2=YUNG-ET2
  22.       EMET3=YUNG-ET3
  23.       CISEL=YUNG/(2.D0*UPUN)
  24.       UMBETA(1)=1.D0-(CISA(1)/CISEL)
  25.       UMBETA(2)=1.D0-(CISA(2)/CISEL)
  26.       UMBETA(3)=1.D0-(CISA(3)/CISEL)
  27. C
  28. C***********************************************************************
  29. C************************ PAS DE FISSURATION ***************************
  30. C***********************************************************************
  31. C
  32.       IF(IENDO.EQ.0) RETURN
  33. C
  34.       GO TO (100,200,300,400,500,600,700),IENDO
  35.       KERRE=468
  36.       RETURN
  37. C
  38. C***********************************************************************
  39. C******************** FISSURATION EN X SEULEMENT ***********************
  40. C***********************************************************************
  41. C
  42.   100 DENO=(EMET1*UMUN)+(ET1*UPUN*UMDN)
  43.       USDENO=1.D0/DENO
  44.       AAA(1,1)=EMET1*UMUN*USDENO
  45.       AAA(1,2)=EMET1*XNU*USDENO
  46.       AAA(1,3)=EMET1*XNU*USDENO
  47.       AAA(4,4)=UMBETA(1)
  48.       AAA(6,6)=UMBETA(3)
  49.       RETURN
  50. C
  51. C***********************************************************************
  52. C******************** FISSURATION EN Y SEULEMENT ***********************
  53. C***********************************************************************
  54. C
  55.   200 DENO=(EMET2*UMUN)+(ET2*UPUN*UMDN)
  56.       USDENO=1.D0/DENO
  57.       AAA(2,1)=EMET2*XNU*USDENO
  58.       AAA(2,2)=EMET2*UMUN*USDENO
  59.       AAA(2,3)=EMET2*XNU*USDENO
  60.       AAA(4,4)=UMBETA(1)
  61.       AAA(5,5)=UMBETA(2)
  62.       RETURN
  63. C
  64. C***********************************************************************
  65. C******************** FISSURATION EN Z SEULEMENT ***********************
  66. C***********************************************************************
  67. C
  68.   400 DENO=(EMET3*UMUN)+(ET3*UPUN*UMDN)
  69.       USDENO=1.D0/DENO
  70.       AAA(3,1)=EMET3*XNU*USDENO
  71.       AAA(3,2)=EMET3*XNU*USDENO
  72.       AAA(3,3)=EMET3*UMUN*USDENO
  73.       AAA(5,5)=UMBETA(2)
  74.       AAA(6,6)=UMBETA(3)
  75.       RETURN
  76. C
  77. C***********************************************************************
  78. C******************** FISSURATION EN XY ENSEMBLE ***********************
  79. C***********************************************************************
  80. C
  81.   300 DENO=(EMET1*EMET2)+(ET1*EMET2*UPUN*UMUN)+(EMET1*ET2*UPUN*UMUN)+
  82.      .     (ET1*ET2*UPUN*UPUN*UMDN)
  83.       USDENO=1.D0/DENO
  84.       AAA(1,1)=EMET1*(EMET2+(ET2*UPUN*UMUN))*USDENO
  85.       AAA(1,2)=EMET1*(XNU*ET2*UPUN)*USDENO
  86.       AAA(1,3)=XNU*EMET1*(EMET2+(ET2*UPUN))*USDENO
  87.       AAA(2,1)=EMET2*(XNU*ET1*UPUN)*USDENO
  88.       AAA(2,2)=EMET2*(EMET1+(ET1*UPUN*UMUN))*USDENO
  89.       AAA(2,3)=XNU*EMET2*(EMET1+(ET1*UPUN))*USDENO
  90.       AAA(4,4)=UMBETA(1)
  91.       AAA(5,5)=UMBETA(2)
  92.       AAA(6,6)=UMBETA(3)
  93.       RETURN
  94. C
  95. C***********************************************************************
  96. C******************** FISSURATION EN YZ ENSEMBLE ***********************
  97. C***********************************************************************
  98. C
  99.   600 DENO=(EMET2*EMET3)+(ET2*EMET3*UPUN*UMUN)+(EMET2*ET3*UPUN*UMUN)+
  100.      .     (ET2*ET3*UPUN*UPUN*UMDN)
  101.       USDENO=1.D0/DENO
  102.       AAA(2,1)=XNU*EMET2*(EMET3+(ET3*UPUN))*USDENO
  103.       AAA(2,2)=EMET2*(EMET3+(ET3*UPUN*UMUN))*USDENO
  104.       AAA(2,3)=EMET2*(XNU*ET3*UPUN)*USDENO
  105.       AAA(3,1)=XNU*EMET3*(EMET2+(ET2*UPUN))*USDENO
  106.       AAA(3,2)=EMET3*(XNU*ET2*UPUN)*USDENO
  107.       AAA(3,3)=EMET3*(EMET2+(ET2*UPUN*UMUN))*USDENO
  108.       AAA(4,4)=UMBETA(1)
  109.       AAA(5,5)=UMBETA(2)
  110.       AAA(6,6)=UMBETA(3)
  111.       RETURN
  112. C
  113. C***********************************************************************
  114. C******************** FISSURATION EN ZX ENSEMBLE ***********************
  115. C***********************************************************************
  116. C
  117.   500 DENO=(EMET3*EMET1)+(ET3*EMET1*UPUN*UMUN)+(EMET3*ET1*UPUN*UMUN)+
  118.      .     (ET3*ET1*UPUN*UPUN*UMDN)
  119.       USDENO=1.D0/DENO
  120.       AAA(1,1)=EMET1*(EMET3+(ET3*UPUN*UMUN))*USDENO
  121.       AAA(1,2)=XNU*EMET1*(EMET3+(ET3*UPUN))*USDENO
  122.       AAA(1,3)=EMET1*(XNU*ET3*UPUN)*USDENO
  123.       AAA(3,1)=EMET3*(XNU*ET1*UPUN)*USDENO
  124.       AAA(3,2)=XNU*EMET3*(EMET1+(ET1*UPUN))*USDENO
  125.       AAA(3,3)=EMET3*(EMET1+(ET1*UPUN*UMUN))*USDENO
  126.       AAA(4,4)=UMBETA(1)
  127.       AAA(5,5)=UMBETA(2)
  128.       AAA(6,6)=UMBETA(3)
  129.       RETURN
  130. C
  131. C***********************************************************************
  132. C******************** FISSURATION EN XYZ ENSEMBLE **********************
  133. C***********************************************************************
  134. C
  135.   700 DENO=(EMET1*EMET2*EMET3)+
  136.      .   (ET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+(EMET1*EMET2*ET3)+
  137.      .   (UPUN*UMUN*((EMET1*ET2*ET3)+(ET1*EMET2*ET3)+(ET1*ET2*EMET3)))+
  138.      .   (UPUN*UPUN*UMDN*(ET1*ET2*ET3))
  139.       USDENO=1.D0/DENO
  140.       AAA(1,1)=((EMET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+(EMET1*EMET2*ET3)+
  141.      .      (EMET1*ET2*ET3*UPUN*UMUN))*USDENO
  142.       AAA(1,2)=XNU*EMET1*ET2*(EMET3+(ET3*UPUN))*USDENO
  143.       AAA(1,3)=XNU*EMET1*ET3*(EMET2+(ET2*UPUN))*USDENO
  144.       AAA(2,1)=XNU*EMET2*ET1*(EMET3+(ET3*UPUN))*USDENO
  145.       AAA(2,2)=((EMET1*EMET2*EMET3)+(ET1*EMET2*EMET3)+(EMET1*EMET2*ET3)+
  146.      .      (ET1*EMET2*ET3*UPUN*UMUN))*USDENO
  147.       AAA(2,3)=XNU*EMET2*ET3*(EMET1+(ET1*UPUN))*USDENO
  148.       AAA(3,1)=XNU*EMET3*ET1*(EMET2+(ET2*UPUN))*USDENO
  149.       AAA(3,2)=XNU*EMET3*ET2*(EMET1+(ET1*UPUN))*USDENO
  150.       AAA(3,3)=((EMET1*EMET2*EMET3)+(ET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+
  151.      .      (ET1*ET2*EMET3*UPUN*UMUN))*USDENO
  152.       AAA(4,4)=UMBETA(1)
  153.       AAA(5,5)=UMBETA(2)
  154.       AAA(6,6)=UMBETA(3)
  155.       RETURN
  156. C
  157.       END
  158.  
  159.  

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