Télécharger clouve.eso

Retour à la liste

Numérotation des lignes :

  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