Télécharger clendo.eso

Retour à la liste

Numérotation des lignes :

  1. C CLENDO SOURCE CHAT 05/01/12 22:04:48 5004
  2. SUBROUTINE CLENDO (DDP,YUNG,XNU,ET,CISA,IENDO,KERRE)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. DIMENSION ET(*),CISA(*),DDP(6,6)
  6. C
  7. C***********************************************************************
  8. C***************** CALCUL DU MODULE D'ENDOMMAGEMENT ********************
  9. C***********************************************************************
  10. C
  11. C INITIALISATION
  12. C
  13. CALL ZDANUL(DDP,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. COEF=YUNG/(UPUN*UMDN)
  24. C
  25. C***********************************************************************
  26. C************************ PAS DE FISSURATION ***************************
  27. C***********************************************************************
  28. C
  29. DDP(1,1)=UMUN*COEF
  30. DDP(1,2)=XNU*COEF
  31. DDP(1,3)=DDP(1,2)
  32. DDP(2,1)=DDP(1,2)
  33. DDP(2,2)=DDP(1,1)
  34. DDP(2,3)=DDP(1,2)
  35. DDP(3,1)=DDP(1,2)
  36. DDP(3,2)=DDP(1,2)
  37. DDP(3,3)=DDP(1,1)
  38. DDP(4,4)=YUNG/(2.D0*UPUN)
  39. DDP(5,5)=DDP(4,4)
  40. DDP(6,6)=DDP(4,4)
  41. IF(IENDO.EQ.0) RETURN
  42. C
  43. GO TO (100,200,300,400,500,600,700),IENDO
  44. KERRE=468
  45. RETURN
  46. C
  47. C***********************************************************************
  48. C******************** FISSURATION EN X SEULEMENT ***********************
  49. C***********************************************************************
  50. C
  51. 100 DENO=(EMET1*UMUN)+(ET1*UPUN*UMDN)
  52. USDENO=1.D0/DENO
  53. DDP(1,1)=YUNG*ET1*UMUN*USDENO
  54. DDP(1,2)=XNU*YUNG*ET1*USDENO
  55. DDP(1,3)=DDP(1,2)
  56. DDP(2,1)=DDP(1,2)
  57. DDP(2,2)=YUNG*((EMET1/UPUN)+(ET1*UMUN))*USDENO
  58. DDP(2,3)=XNU*YUNG*((EMET1/UPUN)+ET1)*USDENO
  59. DDP(3,1)=DDP(1,3)
  60. DDP(3,2)=DDP(2,3)
  61. DDP(3,3)=DDP(2,2)
  62. DDP(4,4)=CISA(1)
  63. DDP(6,6)=CISA(3)
  64. RETURN
  65. C
  66. C***********************************************************************
  67. C******************** FISSURATION EN Y SEULEMENT ***********************
  68. C***********************************************************************
  69. C
  70. 200 DENO=(EMET2*UMUN)+(ET2*UPUN*UMDN)
  71. USDENO=1.D0/DENO
  72. DDP(1,1)=YUNG*((EMET2/UPUN)+(ET2*UMUN))*USDENO
  73. DDP(1,2)=XNU*YUNG*ET2*USDENO
  74. DDP(1,3)=XNU*YUNG*((EMET2/UPUN)+ET2)*USDENO
  75. DDP(2,1)=DDP(1,2)
  76. DDP(2,2)=YUNG*ET2*UMUN*USDENO
  77. DDP(2,3)=DDP(1,2)
  78. DDP(3,1)=DDP(1,3)
  79. DDP(3,2)=DDP(2,3)
  80. DDP(3,3)=DDP(1,1)
  81. DDP(4,4)=CISA(1)
  82. DDP(5,5)=CISA(2)
  83. RETURN
  84. C
  85. C***********************************************************************
  86. C******************** FISSURATION EN Z SEULEMENT ***********************
  87. C***********************************************************************
  88. C
  89. 400 DENO=(EMET3*UMUN)+(ET3*UPUN*UMDN)
  90. USDENO=1.D0/DENO
  91. DDP(1,1)=YUNG*((EMET3/UPUN)+(ET3*UMUN))*USDENO
  92. DDP(1,2)=XNU*YUNG*((EMET3/UPUN)+ET3)*USDENO
  93. DDP(1,3)=XNU*YUNG*ET3*USDENO
  94. DDP(2,1)=DDP(1,2)
  95. DDP(2,2)=DDP(1,1)
  96. DDP(2,3)=DDP(1,3)
  97. DDP(3,1)=DDP(1,3)
  98. DDP(3,2)=DDP(2,3)
  99. DDP(3,3)=YUNG*ET3*UMUN*USDENO
  100. DDP(5,5)=CISA(2)
  101. DDP(6,6)=CISA(3)
  102. RETURN
  103. C
  104. C***********************************************************************
  105. C******************** FISSURATION EN XY ENSEMBLE ***********************
  106. C***********************************************************************
  107. C
  108. 300 DENO=(EMET1*EMET2)+(ET1*EMET2*UPUN*UMUN)+(EMET1*ET2*UPUN*UMUN)+
  109. . (ET1*ET2*UPUN*UPUN*UMDN)
  110. USDENO=1.D0/DENO
  111. DDP(1,1)=YUNG*ET1*(EMET2+(ET2*UPUN*UMUN))*USDENO
  112. DDP(1,2)=XNU*YUNG*ET1*ET2*UPUN*USDENO
  113. DDP(1,3)=XNU*YUNG*ET1*(EMET2+(ET2*UPUN))*USDENO
  114. DDP(2,1)=DDP(1,2)
  115. DDP(2,2)=YUNG*ET2*(EMET1+(ET1*UPUN*UMUN))*USDENO
  116. DDP(2,3)=XNU*YUNG*ET2*(EMET1+(ET1*UPUN))*USDENO
  117. DDP(3,1)=DDP(1,3)
  118. DDP(3,2)=DDP(2,3)
  119. DDP(3,3)=YUNG*((EMET1*EMET2)+(ET1*EMET2)+(EMET1*ET2)+
  120. . (ET1*ET2*UPUN*UMUN))*USDENO
  121. DDP(4,4)=CISA(1)
  122. DDP(5,5)=CISA(2)
  123. DDP(6,6)=CISA(3)
  124. RETURN
  125. C
  126. C***********************************************************************
  127. C******************** FISSURATION EN YZ ENSEMBLE ***********************
  128. C***********************************************************************
  129. C
  130. 600 DENO=(EMET2*EMET3)+(ET2*EMET3*UPUN*UMUN)+(EMET2*ET3*UPUN*UMUN)+
  131. . (ET2*ET3*UPUN*UPUN*UMDN)
  132. USDENO=1.D0/DENO
  133. DDP(1,1)=YUNG*((EMET2*EMET3)+(ET2*EMET3)+(EMET2*ET3)+
  134. . (ET2*ET3*UPUN*UMUN))*USDENO
  135. DDP(1,2)=XNU*YUNG*ET2*(EMET3+(ET3*UPUN))*USDENO
  136. DDP(1,3)=XNU*YUNG*ET3*(EMET2+(ET2*UPUN))*USDENO
  137. DDP(2,1)=DDP(1,2)
  138. DDP(2,2)=YUNG*ET2*(EMET3+(ET3*UPUN*UMUN))*USDENO
  139. DDP(2,3)=XNU*YUNG*ET2*ET3*UPUN*USDENO
  140. DDP(3,1)=DDP(1,3)
  141. DDP(3,2)=DDP(2,3)
  142. DDP(3,3)=YUNG*ET3*(EMET2+(ET2*UPUN*UMUN))*USDENO
  143. DDP(4,4)=CISA(1)
  144. DDP(5,5)=CISA(2)
  145. DDP(6,6)=CISA(3)
  146. RETURN
  147. C
  148. C***********************************************************************
  149. C******************** FISSURATION EN ZX ENSEMBLE ***********************
  150. C***********************************************************************
  151. C
  152. 500 DENO=(EMET3*EMET1)+(ET3*EMET1*UPUN*UMUN)+(EMET3*ET1*UPUN*UMUN)+
  153. . (ET3*ET1*UPUN*UPUN*UMDN)
  154. USDENO=1.D0/DENO
  155. DDP(1,1)=YUNG*ET1*(EMET3+(ET3*UPUN*UMUN))*USDENO
  156. DDP(1,2)=XNU*YUNG*ET1*(EMET3+(ET3*UPUN))*USDENO
  157. DDP(1,3)=XNU*YUNG*ET3*ET1*UPUN*USDENO
  158. DDP(2,1)=DDP(1,2)
  159. DDP(2,2)=YUNG*((EMET3*EMET1)+(ET3*EMET1)+(EMET3*ET1)+
  160. . (ET3*ET1*UPUN*UMUN))*USDENO
  161. DDP(2,3)=XNU*YUNG*ET3*(EMET1+(ET1*UPUN))*USDENO
  162. DDP(3,1)=DDP(1,3)
  163. DDP(3,2)=DDP(2,3)
  164. DDP(3,3)=YUNG*ET3*(EMET1+(ET1*UPUN*UMUN))*USDENO
  165. DDP(4,4)=CISA(1)
  166. DDP(5,5)=CISA(2)
  167. DDP(6,6)=CISA(3)
  168. RETURN
  169. C
  170. C***********************************************************************
  171. C******************** FISSURATION EN XYZ ENSEMBLE **********************
  172. C***********************************************************************
  173. C
  174. 700 DENO=(EMET1*EMET2*EMET3)+
  175. . (ET1*EMET2*EMET3)+(EMET1*ET2*EMET3)+(EMET1*EMET2*ET3)+
  176. . (UPUN*UMUN*((EMET1*ET2*ET3)+(ET1*EMET2*ET3)+(ET1*ET2*EMET3)))+
  177. . (UPUN*UPUN*UMDN*(ET1*ET2*ET3))
  178. USDENO=1.D0/DENO
  179. DDP(1,1)=((EMET2*EMET3)+(ET2*EMET3)+(EMET2*ET3)+
  180. . (ET2*ET3*UPUN*UMUN))*YUNG*ET1*USDENO
  181. DDP(1,2)=XNU*YUNG*ET1*ET2*(EMET3+(ET3*UPUN))*USDENO
  182. DDP(1,3)=XNU*YUNG*ET1*ET3*(EMET2+(ET2*UPUN))*USDENO
  183. DDP(2,1)=DDP(1,2)
  184. DDP(2,2)=((EMET1*EMET3)+(ET1*EMET3)+(EMET1*ET3)+
  185. . (ET1*ET3*UPUN*UMUN))*YUNG*ET2*USDENO
  186. DDP(2,3)=XNU*YUNG*ET2*ET3*(EMET1+(ET1*UPUN))*USDENO
  187. DDP(3,1)=DDP(1,3)
  188. DDP(3,2)=DDP(2,3)
  189. DDP(3,3)=((EMET1*EMET2)+(ET1*EMET2)+(EMET1*ET2)+
  190. . (ET1*ET2*UPUN*UMUN))*YUNG*ET3*USDENO
  191. DDP(4,4)=CISA(1)
  192. DDP(5,5)=CISA(2)
  193. DDP(6,6)=CISA(3)
  194. RETURN
  195. C
  196. END
  197.  
  198.  

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