Télécharger g_decouplage_1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : g_decouplage_1.dgibi
  2. ****************************************************
  3. ****************************************************
  4. * *
  5. * VERIFICATION DE LA PROCEDURE G_THETA *
  6. * POUR LE CALCUL DES FIC POUR UNE FISSURE *
  7. * PENNY-SHAPED DANS UN CYLINDRE *
  8. * - *
  9. * VERIFICATION DU CALCUL DE KI VIA DECOUPLAGE *
  10. * EN 3D ET 2D AXISYMETRIQUE AVEC UNE *
  11. * SOLUTION ANALYTIQUE *
  12. * *
  13. * SOLUTION ANALYTIQUE TIREE DE : *
  14. * "PENNY-SHAPED CRACKS", GUIDERA ET LARDNER 1974 *
  15. * *
  16. ****************************************************
  17.  
  18. * I - CALCUL 3D COMPLET
  19. * ---------------------
  20.  
  21. OPTI 'DIME' 3 'ELEM' 'CUB8' ;
  22. BTRAC = FAUX ;
  23.  
  24. * MAILLAGE AUTOUR DE LA FISSURE
  25. L1 H1 LP = 10. 5. 5. ;
  26. L2 = L1 + 10. ;
  27. H2 = H1 + 5. ;
  28. DENS1 = 1. ;
  29. DENS DENS1 ;
  30.  
  31. P1 = DENS1 0. 0. ;
  32. P2 = DENS1 H1 0. ;
  33. P3 = L1 H1 0. ;
  34. P4 = L1 0. 0. ;
  35. D1 = DROI P1 P2 ;
  36. S1 = D1 TRAN ((0.-DENS1) 0. 0.) ;
  37. N1 = ((ENTI ((L1 - DENS1) / DENS1)) / 2) * 2 ;
  38. D2 = DROI N1 P2 P3 ;
  39. D3 = DROI P3 P4 ;
  40.  
  41. P5 = (LP - DENS1) 0. 0. ;
  42. P6 = (LP - DENS1) DENS1 0. ;
  43. P7 = (LP + DENS1) DENS1 0. ;
  44. P8 = (LP + DENS1) 0. 0. ;
  45. D4 = DROI (NBEL D1) P5 P6 ;
  46. D5 = DROI (NBEL D2) P6 P7 ;
  47. D6 = DROI (NBEL D3) P7 P8 ;
  48. DENS2 = DENS1 / H1 ;
  49. S2 = (D1 ET D2 ET D3) REGL (D4 ET D5 ET D6) 'DINI' DENS1 'DFIN' DENS2 ;
  50.  
  51. PP = LP 0. 0. ;
  52. D7 = DROI (N1/2) (DROI (N1/2) P8 PP) P5 ;
  53. S3 = DALL D4 D5 D6 D7 'PLAN' ;
  54. SURF1 = S1 ET S2 ET S3 ;
  55.  
  56. PSUP = (COOR 2 SURF1) POIN 'MAXI' ;
  57. ELSUP = (CONT SURF1) ELEM 'APPUYE' PSUP ;
  58. S4 = ELSUP TRAN (0. (H2 - H1) 0.) ;
  59. SURF1 = SURF1 ET S4 ;
  60. PDROI = (COOR 1 SURF1) POIN 'MAXI' ;
  61. ELDROI = (CONT SURF1) ELEM 'APPUYE' PDROI ;
  62. S5 = ELDROI TRAN ((L2 - L1) 0. 0.) ;
  63. SURF1 = SURF1 ET S5 ;
  64.  
  65. NFRON = 18 ;
  66. VOLU1 = SURF1 VOLU NFRON 'ROTA' 360. (0. 1. 0.) (0. 0. 0.) ;
  67. ELIM VOLU1 1.E-10 ;
  68. VOLU1 = REGE VOLU1 ;
  69. VOLU2 = VOLU1 SYME 'PLAN' (0. 0. 0.) (1. 0. 0.) (0. 0. 1.) ;
  70. OMEGA1 = VOLU1 ET VOLU2 ;
  71.  
  72. * DEFINITION DE LA FISSURE ET FUSION DES NOEUDS COINCIDENTS
  73. * PARTIE SUP
  74. ENVE1 = ENVE VOLU1 ;
  75. POIN1 = (COOR 2 ENVE1) POIN 'MINI' ;
  76. ELEM1 = ENVE1 ELEM 'APPUYE' POIN1 ;
  77. R1 = (((COOR 1 ELEM1)**2) + ((COOR 3 ELEM1)**2))**0.5 ;
  78. POIN1 = R1 POIN 'EGINFE' (LP + 0.02) ;
  79. LVSUP = ELEM1 ELEM 'APPUYE' POIN1 ;
  80. LIGSUP = ELEM1 DIFF LVSUP ;
  81. * PARTIE INF
  82. ENVE2 = ENVE VOLU2 ;
  83. POIN2 = (COOR 2 ENVE2) POIN 'MAXI' ;
  84. ELEM2 = ENVE2 ELEM 'APPUYE' POIN2 ;
  85. R2 = (((COOR 1 ELEM2)**2) + ((COOR 3 ELEM2)**2))**0.5 ;
  86. POIN2 = R2 POIN 'EGINFE' (LP + 0.02) ;
  87. LVINF = ELEM2 ELEM 'APPUYE' POIN2 ;
  88. LIGINF = ELEM2 DIFF LVINF ;
  89. * ELIMINATION
  90. ELIM LIGSUP LIGINF 1.E-10 ;
  91. * DEFINITION DU FRONT DE FISSURE
  92. FRONFISS = (CONT LIGSUP) INTE (CONT LVSUP) ;
  93. SI BTRAC ;
  94. TRAC SURF1 ;
  95. TRAC 'CACH' (VOLU1 ET (LVSUP COUL 'ROUG') ET (LIGSUP COUL 'VERT')) ;
  96. FINSI ;
  97.  
  98. * MODELE ET MATERIAU
  99. MOD1 = MODE OMEGA1 'MECANIQUE' 'ELASTIQUE' ;
  100. YOUN1 NU1 = 2.E5 0.3 ;
  101. MAT1 = MATE MOD1 'YOUN' YOUN1 'NU' NU1 ;
  102. RIG1 = RIGI MOD1 MAT1 ;
  103.  
  104. * CL EN EFFORTS
  105. SIG0 = ZERO MOD1 'CONTRAIN' ;
  106. SIG1 = MANU 'CHML' MOD1 'SMYY' 1. 'STRESSES' ;
  107. SIG2 = SIG0 + SIG1 ;
  108. DOMEGA1 = ENVE OMEGA1 ;
  109. DOMEGA1 = DIFF DOMEGA1 (LVSUP ET LVINF) ;
  110. T1 = REDU (BSIG MOD1 SIG2) DOMEGA1 ;
  111.  
  112. * CL EN DEPLACEMENTS
  113. PSUP = (COOR 2 DOMEGA1) POIN 'MAXI' ;
  114. P1 = PSUP POIN 'DROIT' (0. 0. 0.) (0. 1. 0.) ;
  115. PINF = (COOR 2 DOMEGA1) POIN 'MINI' ;
  116. P2 = PINF POIN 'DROIT' (0. 0. 0.) (0. 1. 0.) ;
  117. PMIL = DOMEGA1 POIN 'PLAN' (0. 0. 0.) (1. 0. 0.) (0. 0. 1.) ;
  118. P3 = (COOR 1 PMIL) POIN 'MAXI' ;
  119. P4 = (COOR 1 PMIL) POIN 'MINI' ;
  120. BLOQ1 = BLOQ (MOTS 'UX' 'UZ') (P1 ET P2) ;
  121. BLOQ3 = BLOQ 'UZ' P3 ;
  122. RELA3 = RELA 'UY' P3 + 'UY' P4 ;
  123. BLOQ0 = BLOQ1 ET BLOQ3 ET RELA3 ;
  124.  
  125. * RESO
  126. U1 = RESO (RIG1 ET BLOQ0) T1 ;
  127.  
  128. * G_THETA
  129. SUPTAB = TABL ;
  130. SUPTAB.'MODELE' = MOD1 ;
  131. SUPTAB.'CARACTERISTIQUES' = MAT1 ;
  132. SUPTAB.'BLOCAGES_MECANIQUES' = BLOQ0 ;
  133. SUPTAB.'CHARGEMENTS_MECANIQUES' = T1 ;
  134. SUPTAB.'OBJECTIF' = MOT 'DECOUPLAGE' ;
  135. SUPTAB.'LEVRE_SUPERIEURE' = LVSUP ;
  136. SUPTAB.'LEVRE_INFERIEURE' = LVINF ;
  137. SUPTAB.'COUCHE' = 4 ;
  138. SUPTAB.'FRONT_FISSURE' = FRONFISS ;
  139. SUPTAB.'SOLUTION_RESO' = U1 ;
  140.  
  141. G_THETA SUPTAB ;
  142.  
  143. KI3D = PROG ;
  144. KSUPTAB = SUPTAB.'RESULTATS'.'I' ;
  145. IND1 = INDE KSUPTAB ;
  146. REPE IFRON NFRON ;
  147. PFRON = IND1.&IFRON ;
  148. KI3D = KI3D ET (PROG KSUPTAB.PFRON) ;
  149. FIN IFRON ;
  150.  
  151.  
  152. * II - CALCUL 2D AXISYMETRIQUE
  153. * ----------------------------
  154.  
  155. OPTI 'DIME' 2 'MODE' 'AXIS' 'ELEM' 'QUA4' ;
  156.  
  157. L3 = L2 + 20. ;
  158. H3 = H2 + 10. ;
  159.  
  160. PSUP = (COOR 2 SURF1) POIN 'MAXI' ;
  161. ELSUP = (CONT SURF1) ELEM 'APPUYE' PSUP ;
  162. S4 = ELSUP TRAN (0. (H3 - H2)) ;
  163. SURF1 = SURF1 ET S4 ;
  164. PDROI = (COOR 1 SURF1) POIN 'MAXI' ;
  165. ELDROI = (CONT SURF1) ELEM 'APPUYE' PDROI ;
  166. S5 = ELDROI TRAN ((L3 - L2) 0.) ;
  167. SURF1 = SURF1 ET S5 ;
  168.  
  169. SURF2 = SURF1 SYME 'DROI' (0. 0.) (1. 0.) ;
  170. OMEGA2 = SURF1 ET SURF2 ;
  171.  
  172. * DEFINITION DE LA FISSURE ET FUSION DES NOEUDS COINCIDENTS
  173. C1 = CONT SURF1 ;
  174. POIN1 = (COOR 2 C1) POIN 'MINI' ;
  175. ELEM1 = C1 ELEM 'APPUYE' POIN1 ;
  176. POIN1 = (COOR 1 POIN1) POIN 'EGINFE' LP ;
  177. LVSUP = C1 ELEM 'APPUYE' POIN1 ;
  178. LIGSUP = ELEM1 DIFF LVSUP ;
  179. C2 = CONT SURF2 ;
  180. POIN2 = (COOR 2 C2) POIN 'MAXI' ;
  181. ELEM2 = C2 ELEM 'APPUYE' POIN2 ;
  182. POIN2 = (COOR 1 POIN2) POIN 'EGINFE' LP ;
  183. LVINF = C2 ELEM 'APPUYE' POIN2 ;
  184. LIGINF = ELEM2 DIFF LVINF ;
  185. ELIM LIGSUP LIGINF 1.E-10 ;
  186. SI BTRAC ;
  187. TRAC (OMEGA2 ET (LVSUP COUL 'ROUG')) ;
  188. FINSI ;
  189.  
  190. * DEFINITION DU FRONT DE FISSURE
  191. FRONFISS = ((COOR 1 LVSUP) POIN 'MAXI') POIN 1 ;
  192.  
  193. * MODELE ET MATERIAU
  194. MOD2 = MODE OMEGA2 'MECANIQUE' 'ELASTIQUE' ;
  195. MAT2 = MATE MOD2 'YOUN' YOUN1 'NU' NU1 ;
  196. RIG2 = RIGI MOD2 MAT2 ;
  197.  
  198. * CL EN EFFORTS
  199. SIG0 = ZERO MOD2 'CONTRAIN' ;
  200. SIG1 = MANU 'CHML' MOD2 'SMZZ' 1. 'STRESSES' ;
  201. SIG2 = SIG0 + SIG1 ;
  202. DOMEGA2 = CONT OMEGA2 ;
  203. DOMEGA2 = DIFF DOMEGA2 (LVSUP ET LVINF) ;
  204. T2 = REDU (BSIG MOD2 SIG2) DOMEGA2 ;
  205.  
  206. * CL EN DEPLACEMENTS
  207. ZOM = COOR 2 OMEGA2 ;
  208. PSUP = ZOM POIN 'MAXI' ;
  209. P1 = PSUP POIN 'DROIT' (0. 0.) (0. 1.) ;
  210. PINF = ZOM POIN 'MINI' ;
  211. P2 = PINF POIN 'DROIT' (0. 0.) (0. 1.) ;
  212. RELA1 = RELA 'UZ' P1 + 'UZ' P2 ;
  213. BLOQ0 = RELA1 ;
  214.  
  215. * RESO
  216. U2 = RESO (RIG2 ET BLOQ0) T2 ;
  217.  
  218. * G_THETA
  219. SUPTAB = TABL ;
  220. SUPTAB.'MODELE' = MOD2 ;
  221. SUPTAB.'CARACTERISTIQUES' = MAT2 ;
  222. SUPTAB.'BLOCAGES_MECANIQUES' = BLOQ0 ;
  223. SUPTAB.'CHARGEMENTS_MECANIQUES' = T2 ;
  224. SUPTAB.'OBJECTIF' = MOT 'DECOUPLAGE' ;
  225. SUPTAB.'LEVRE_SUPERIEURE' = LVSUP ;
  226. SUPTAB.'LEVRE_INFERIEURE' = LVINF ;
  227. SUPTAB.'COUCHE' = 5 ;
  228. SUPTAB.'FRONT_FISSURE' = FRONFISS ;
  229. SUPTAB.'SOLUTION_RESO' = U2 ;
  230.  
  231. G_THETA SUPTAB ;
  232. KI2D = PROG NFRON * (SUPTAB.'RESULTATS'.'I') ;
  233.  
  234.  
  235. * III - SOLUTION ANALYTIQUE
  236. * -------------------------
  237.  
  238. KIANA = 2. / PI * ((PI * LP)**0.5) ;
  239. MESS 'SOLUTION ANALYTIQUE :' KIANA ;
  240. KIANA = PROG NFRON * KIANA ;
  241.  
  242.  
  243. * IV - COMPARAISON
  244. * ----------------
  245.  
  246. KINUM = TABL ;
  247. KINUM.'2D' = KI2D ;
  248. KINUM.'3D' = KI3D ;
  249. ACCEPT = 2.5E-2 ;
  250. CMAX1 = MAXI KIANA 'ABS' ;
  251. DIM = 2 ;
  252. REPE I 2 ;
  253. MDIM = CHAI DIM 'D' ;
  254. CRIT1 = MAXI (KINUM.MDIM - KIANA) 'ABS' ;
  255. MESS 'ERREUR RELATIVE = ' (CRIT1 / CMAX1 * 100.) '%' ;
  256. SI (CRIT1 > (CMAX1*ACCEPT)) ;
  257. MESS 'LA SOLUTION' ' ' MDIM ' EST TROP ELOIGNEE DE LA SOLUTION' ;
  258. MESS 'ANALYTIQUE POUR LE CALCUL DU FIC EN MODE I' ;
  259. ERRE 5 ;
  260. FINSI ;
  261. DIM = DIM + 1 ;
  262. FIN I ;
  263.  
  264. FIN ;
  265.  
  266.  
  267.  

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