Télécharger ottoce.eso

Retour à la liste

Numérotation des lignes :

  1. C OTTOCE SOURCE CHAT 05/01/13 02:07:11 5004
  2. SUBROUTINE OTTOCE(MC,MM,SIGMA,DX,DXV1,W,WMAX,SMAX,
  3. & WRUPT,XLTR,XINVL,BTR,NFISSU,NVF,FCRIT,VF,YOUN,
  4. & PRECIZ,JRESU,XCOMP,XLAMC,DFF,DGG,KERRE)
  5. C=========================================================================
  6. C
  7. C 1 NOUVELLE FISSURATION SELON LA DIRECTION 1
  8. C 2 NOUVELLE FISSURATION SELON LA DIRECTION 2
  9. C 3 NOUVELLE FISSURATION SELON LA DIRECTION 3
  10. C 4 POURSUITE DE LA FISSURATION SELON LA DIRECTION 1
  11. C 5 POURSUITE DE LA FISSURATION SELON LA DIRECTION 2
  12. C 6 POURSUITE DE LA FISSURATION SELON LA DIRECTION 3
  13. C 7 RUPTURE SELON LA DIRECTION 1
  14. C 8 RUPTURE SELON LA DIRECTION 2
  15. C 9 RUPTURE SELON LA DIRECTION 3
  16. C 10 REOUVERTURE SELON LA DIRECTION 1
  17. C 11 REOUVERTURE SELON LA DIRECTION 2
  18. C 12 REOUVERTURE SELON LA DIRECTION 3
  19. C 13 REFERMETURE SELON LA DIRECTION 1
  20. C 14 REFERMETURE SELON LA DIRECTION 2
  21. C 15 REFERMETURE SELON LA DIRECTION 3
  22. C 16 COMPRESSION
  23. C
  24. C ENTREES :
  25. C SIGMA,W,WMAX,SMAX,WRUPT,XLTR,XINVL,NFISSU,NVF
  26. C
  27. C
  28. C SORTIES :
  29. C FCRIT,VF
  30. C
  31. C==========================================================================
  32. C
  33. IMPLICIT INTEGER(I-N)
  34. IMPLICIT REAL*8(A-H,O-Z)
  35. -INC CCOPTIO
  36. C
  37. PARAMETER (XZER=0.D0)
  38. C
  39. DIMENSION SIGMA(6),W(3),WMAX(3),WRUPT(3),XLTR(3),XINVL(3)
  40. DIMENSION FCRIT(*),MM(*),SMAX(*),DX(*),DXV1(*)
  41. DIMENSION FCF(3),VF(3,3),JRESU(3),XCOMP(*)
  42. DIMENSION DFF(*),DGG(*)
  43.  
  44. C
  45. C INITIALISATIONS
  46. C
  47. KERRE=0
  48. DO IC=1,16
  49. FCRIT(IC)=-1.D4*PRECIZ
  50. ENDDO
  51. CALL IANUL(JRESU,3)
  52. *
  53. *
  54. LAPPEL=0
  55. DO IC=1,MC
  56. JC=MM(IC)
  57. IF(JC.GE.1.AND.JC.LE.3) THEN
  58. LAPPEL=1
  59. GO TO 200
  60. ENDIF
  61. ENDDO
  62. *
  63. 200 CONTINUE
  64. *
  65. IF(XINVL(1)*XINVL(2)*XINVL(3).EQ.0.D0.AND.LAPPEL.EQ.1) THEN
  66. CALL OTTOFI(SIGMA,NFISSU,NVF,XLTR,XINVL,
  67. & VF,JRESU,FCF,PRECIZ,KERRE)
  68. IF(KERRE.NE.0) RETURN
  69. IF(IIMPI.EQ.42) THEN
  70. WRITE(IOIMP,76600)
  71. 76600 FORMAT(// 2X, ' OTTOCE - APRES OTTOFI ' /)
  72. WRITE(IOIMP,76601) (JRESU(IC),IC=1,3)
  73. 76601 FORMAT(2X, ' JRESU ',2X,3I3/)
  74. ENDIF
  75. ENDIF
  76. *
  77. *
  78. DO IC=1,MC
  79. JC=MM(IC)
  80. *
  81. GO TO (1,1,1,4,4,4,7,7,7,10,10,10,13,13,13,16),JC
  82. *
  83. KERRE=99
  84. RETURN
  85. *
  86. 1 CONTINUE
  87. *
  88. * nouvelle fissure selon une direction
  89. * ------------------------------------
  90. *
  91. IF(JRESU(JC).NE.0) THEN
  92. FCRIT(JC)=FCF(JC)
  93. ENDIF
  94. GO TO 100
  95. *
  96. 4 CONTINUE
  97. KC=JC-3
  98. *
  99. * poursuite de la fissuration selon la direction KC ( 1, 2 OU 3 )
  100. * ---------------------------------------------------------------
  101. *
  102. FCRIT(JC)= SIGMA(KC) - SMAX(KC)
  103. GO TO 100
  104. *
  105. 7 CONTINUE
  106. KC=JC-6
  107. *
  108. * rupture selon la direction KC ( 1, 2 OU 3 )
  109. * -------------------------------------------
  110. *
  111. FCRIT(JC)= -SIGMA(KC)
  112. GO TO 100
  113. *
  114. 10 CONTINUE
  115. KC=JC-9
  116. *
  117. * reouverture selon la direction KC ( 1, 2 OU 3 )
  118. * -----------------------------------------------
  119. *
  120. FCRIT(JC)= SIGMA(KC)
  121. GO TO 100
  122. *
  123. 13 CONTINUE
  124. KC=JC-12
  125. *
  126. * refermeture selon la direction KC ( 1, 2 OU 3 )
  127. * -----------------------------------------------
  128. *
  129. WREOUV=BTR*MIN(WMAX(KC),WRUPT(KC))
  130. WAUX=W(KC)+DXV1(KC)/XINVL(KC)
  131. FCRIT(JC)=YOUN*(WREOUV-WAUX)*XINVL(KC)
  132. GO TO 100
  133. *
  134. 16 CONTINUE
  135. *
  136. * compression
  137. * -----------
  138. *
  139. PRECIE=1.D-10
  140. CALL OTTOCP(SIGMA,FCR4,XLTR,DFF,DGG,H4,
  141. & PRECIE,PRECIZ,XCOMP,XLAMC,KERRE)
  142.  
  143.  
  144. FCRIT(JC) = FCR4
  145. GO TO 100
  146. *
  147. 100 CONTINUE
  148. *
  149. ENDDO
  150. *
  151. * IF(IIMPI.EQ.42) THEN
  152. * WRITE(IOIMP,77000)
  153. *77000 FORMAT(// 2X, ' EN SORTIE DE OTTOCE ' /)
  154. * DO IC=1,MC
  155. * JC=MM(IC)
  156. * WRITE(IOIMP,77001) IC,JC,FCRIT(JC)
  157. *77001 FORMAT( 2X, ' OTTOCE - IC=',I3,2X,'JC=',I3,2X,
  158. * & 'FCRIT=',1PE12.5/)
  159. * ENDDO
  160. * ENDIF
  161. *
  162. RETURN
  163. END
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  

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