Télécharger t2.eso

Retour à la liste

Numérotation des lignes :

  1. C T2 SOURCE CHAT 05/01/13 03:30:00 5004
  2. SUBROUTINE T2(SIG,DSIG,YOUN,ANU,RT1,RT2,H2,RDP,ADP,
  3. 1 XLAM1,XLAM2,IDAM,ANG)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. -INC CCOPTIO
  7. DIMENSION SIG(3),DSIG(3),SFC(3),DSFC(3),TSFC(3),SFG(3),DSFG(3),
  8. 1 EPC(3),A(4),B(2),IDAM(3),SG(3),DSG(3)
  9. Y=YOUN/(1.D0-ANU*ANU)
  10. IF(RT2.LT.1.E-10) THEN
  11. H2=0.D0
  12. RT2=0.D0
  13. ENDIF
  14. C
  15. C--------------------------------------------
  16. C BOITE D ECOULEMENT DE TRACTION 2
  17. C--------------------------------------------
  18. C
  19. C-----------------------------------------------
  20. C ON SE PLACE DANS LE REPERE DE FISSURATION
  21. C-----------------------------------------------
  22. C
  23. CALL CHREP(SIG,SFG,ANG)
  24. CALL CHREP(DSIG,DSFG,ANG)
  25. IF(IIMPI.EQ.9) WRITE(IOIMP,9999)
  26. 9999 FORMAT(1X,'T2 ECOULEMENT SUIVANT TRACTION 2',/)
  27. 1 DL2=(SFG(2)+DSFG(2)-RT2)/(Y-H2)
  28. IF(DL2.LT.0.D0) THEN
  29. IDAM(2)=0
  30. RETURN
  31. ENDIF
  32. C
  33. C----------------------------------------------------------------
  34. C CAS OU LA LIMITE EN TRACTION RESTE SUPERIEUR A 0 (1000)
  35. C----------------------------------------------------------------
  36. C
  37. IF((RT2-H2*DL2).GE.-1.E-10) GOTO 1000
  38. C
  39. C--------------------------------------------------------
  40. C CAS OU LA LIMITE EN TRACTION DEVIENT NULLE
  41. C--------------------------------------------------------
  42. C
  43. DL2=RT2/H2
  44. EPC(2)=DL2
  45. EPC(1)=0.D0
  46. EPC(3)=0.D0
  47. CALL CPHOOB(EPC,SFC,YOUN,ANU)
  48. X=0.D0
  49. IF(DSFG(2).NE.0.D0) X=(SFC(2)-SFG(2))/DSFG(2)
  50. DO 10 ITYP=1,3
  51. SFC(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-SFC(ITYP)
  52. 10 DSFC(ITYP)=(1.D0-X)*DSFG(ITYP)
  53. CALL CDP(SFC,ADP,RDP,VCDP)
  54. CALL CTRAF(SFC(1),RT1,VCTR1)
  55. CALL CCOAF(SFC(1),XLAM1,VCCO1)
  56. IF(VCTR1.GT.0.D0) GOTO 2000
  57. IF(VCCO1.GT.0.D0) GOTO 3000
  58. IF(VCDP.GT.0.D0) GOTO 4000
  59. RT2=0.D0
  60. H2=0.D0
  61. XLAM2=XLAM2+DL2
  62. DO 20 ITYP=1,3
  63. SFG(ITYP)=SFC(ITYP)
  64. 20 DSFG(ITYP)=DSFC(ITYP)
  65. GOTO 1
  66. 1000 EPC(2)=DL2
  67. EPC(1)=0.D0
  68. EPC(3)=0.D0
  69. CALL CPHOOB(EPC,SFC,YOUN,ANU)
  70. DO 1010 ITYP=1,3
  71. 1010 SFC(ITYP)=SFG(ITYP)+DSFG(ITYP)-SFC(ITYP)
  72. C
  73. C-----------------------------------------------
  74. C ON REGARDE SI LES AUTRES CRITERES NE
  75. C SONT PAS ENDOMMAGES PENDANT L ECOULEMENT
  76. C-----------------------------------------------
  77. C
  78. CALL CDP(SFC,ADP,RDP,VCDP)
  79. CALL CTRAF(SFC(1),RT1,VCTR1)
  80. CALL CCOAF(SFC(1),XLAM1,VCCO1)
  81. C
  82. C------------------------------------
  83. C CRITERE DE TRACTION ENDOMMAGE
  84. C------------------------------------
  85. C
  86. IF(VCTR1.GT.0.D0) GOTO 2000
  87. C
  88. C----------------------------------------
  89. C CRITERE DE COMPRESSION ENDOMMAGE
  90. C----------------------------------------
  91. C
  92. IF(VCCO1.GT.0.D0) GOTO 3000
  93. C
  94. C------------------------------------
  95. C CRITERE DE DRUCKER ENDOMMAGE
  96. C------------------------------------
  97. C
  98. IF(VCDP.GT.0.D0) GOTO 4000
  99. C
  100. C--------------------------------------------------------
  101. C AUCUN AUTRE CRITERE ENDOMMAGE PENDANT L ECOULEMENT
  102. C--------------------------------------------------------
  103. C
  104. CALL CHREP(SFC,SIG,-ANG)
  105. DO 1020 ITYP=1,3
  106. 1020 DSIG(ITYP)=0.D0
  107. RT2=RT2-H2*DL2
  108. IF(RT2.LE.1.E-10) THEN
  109. H2=0.D0
  110. RT2=0.D0
  111. ENDIF
  112. IDAM(1)=0
  113. IDAM(2)=0
  114. IDAM(3)=0
  115. XLAM2=XLAM2+DL2
  116. RETURN
  117. C
  118. C----------------------------------------
  119. C CAS OU EN ENDOMMAGE LA TRACTION1
  120. C----------------------------------------
  121. C
  122. 2000 A(1)=DSFG(1)
  123. A(3)=DSFG(2)
  124. A(2)=-Y*ANU
  125. A(4)=-Y+H2
  126. DET=A(1)*A(4)-A(2)*A(3)
  127. IF(DET.EQ.0.D0) THEN
  128. IDAM(1)=1
  129. IDAM(2)=1
  130. IDAM(3)=0
  131. RETURN
  132. ENDIF
  133. B(1)=RT1-SFG(1)
  134. B(2)=RT2-SFG(2)
  135. CALL SYLIN2(A,B,X,DL2)
  136. EPC(1)=0.D0
  137. EPC(2)=DL2
  138. EPC(3)=0.D0
  139. CALL CPHOOB(EPC,SFC,YOUN,ANU)
  140. DO 2010 ITYP=1,3
  141. DSFC(ITYP)=(1.D0-X)*DSFG(ITYP)
  142. 2010 SFC(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-SFC(ITYP)
  143. CALL CDP(SFC,ADP,RDP,VCDP)
  144. IF(VCDP.GT.0.D0) GOTO 4000
  145. RT2=RT2-H2*DL2
  146. IF(RT2.LE.1.E-10) THEN
  147. H2=0.D0
  148. RT2=0.D0
  149. ENDIF
  150. IDAM(1)=1
  151. IDAM(2)=1
  152. IDAM(3)=0
  153. XLAM2=XLAM2+DL2
  154. CALL CHREP(DSFC,DSIG,-ANG)
  155. CALL CHREP(SFC,SIG,-ANG)
  156. RETURN
  157. C
  158. C----------------------------------------
  159. C CAS OU EN ENDOMMAGE LA COMPRESSION 1
  160. C----------------------------------------
  161. C
  162. 3000 A(1)=DSFG(1)
  163. A(3)=DSFG(2)
  164. A(2)=-Y*ANU
  165. A(4)=-Y+H2
  166. DET=A(1)*A(4)-A(2)*A(3)
  167. IF(DET.EQ.0.D0) THEN
  168. IDAM(1)=-1
  169. IDAM(2)=1
  170. IDAM(3)=0
  171. RETURN
  172. ENDIF
  173. B(1)=-SFG(1)
  174. B(2)=RT2-SFG(2)
  175. CALL SYLIN2(A,B,X,DL2)
  176. EPC(1)=0.D0
  177. EPC(2)=DL2
  178. EPC(3)=0.D0
  179. CALL CPHOOB(EPC,SFC,YOUN,ANU)
  180. DO 3010 ITYP=1,3
  181. DSFC(ITYP)=(1.D0-X)*DSFG(ITYP)
  182. 3010 SFC(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-SFC(ITYP)
  183. CALL CDP(SFC,ADP,RDP,VCDP)
  184. IF(VCDP.GT.0.D0) GOTO 4000
  185. RT2=RT2-H2*DL2
  186. IF(RT2.LE.1.E-10) THEN
  187. H2=0.D0
  188. RT2=0.D0
  189. ENDIF
  190. IDAM(1)=-1
  191. IDAM(2)=1
  192. IDAM(3)=0
  193. XLAM2=XLAM2+DL2
  194. CALL CHREP(DSFC,DSIG,-ANG)
  195. CALL CHREP(SFC,SIG,-ANG)
  196. RETURN
  197. C
  198. C----------------------------------------
  199. C CAS OU EN ENDOMMAGE LE DP
  200. C----------------------------------------
  201. C
  202. 4000 CONTINUE
  203. EPC(1)=0.D0
  204. EPC(2)=1.D0
  205. EPC(3)=0.D0
  206. CALL CPHOOB(EPC,SFC,YOUN,ANU)
  207. E=Y-H2
  208. F=(SFG(2)-RT2)/E
  209. G=DSFG(2)/E
  210. DO 4010 ITYP=1,3
  211. SG(ITYP)=SFG(ITYP)-F*SFC(ITYP)
  212. 4010 DSG(ITYP)=DSFG(ITYP)-G*SFC(ITYP)
  213. CALL XDP(SG,DSG,RDP,ADP,X,ITEST)
  214. IF(ITEST.EQ.1) THEN
  215. IDAM(1)=0
  216. IDAM(2)=1
  217. IDAM(3)=1
  218. RETURN
  219. ENDIF
  220. DL2=(SFG(2)+X*DSFG(2)-RT2)/(Y-H2)
  221. EPC(1)=0.D0
  222. EPC(2)=DL2
  223. EPC(3)=0.D0
  224. RT2=RT2-H2*DL2
  225. IF(RT2.LE.1.E-10) THEN
  226. H2=0.D0
  227. RT2=0.D0
  228. ENDIF
  229. CALL CPHOOB(EPC,SFC,YOUN,ANU)
  230. DO 4020 ITYP=1,3
  231. SFC(ITYP)=X*DSFG(ITYP)-SFC(ITYP)+SFG(ITYP)
  232. 4020 DSFC(ITYP)=DSFG(ITYP)*(1.D0-X)
  233. XLAM2=XLAM2+DL2
  234. CALL CHREP(DSFC,DSIG,-ANG)
  235. CALL CHREP(SFC,SIG,-ANG)
  236. IDAM(1)=0
  237. IDAM(2)=1
  238. IDAM(3)=1
  239. RETURN
  240. END
  241.  
  242.  

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