Télécharger pore2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : pore2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Pore2.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. complet = faux;
  8. *************************************************
  9. * pour calcul complet mettre complet à : vrai;
  10. ***************************************************
  11. *******************************************************
  12. * *
  13. * TEST PORE2 *
  14. * *
  15. * CYLINDRE EPAIS EN MILIEU POREUX *
  16. * AVEC EFFETS THERMIQUES *
  17. * *
  18. * REFERENCE : Benchmark INTERCLAY 1.4 Basic Model *
  19. * *
  20. * Le milieu est elastoplastique, *
  21. * modele Von-Mises associe *
  22. * *
  23. * On decharge la contrainte a l'interieur du *
  24. * trou puis on impose un champ de temperature *
  25. * *
  26. * On compare avec les valeurs de reference *
  27. * obtenues aussi par CASTEM2000 et en accord *
  28. * avec les resultats INTERCLAY *
  29. * *
  30. *******************************************************
  31. *
  32. OPTI ECHO 0;
  33. * ajout de option epsilon lineaire pour la precision des test!
  34. OPTION epsilon lineaire;
  35.  
  36. TITRE ' TEST PORE2 ' ;
  37. GRAPH = 'N';
  38. OPTIO DIME 2 MODEL AXIS ELEM QUA8 ;
  39. *
  40. *=====================================================*
  41. * geometrie : maillage *
  42. *=====================================================*
  43. *
  44. * Dimension en metres
  45. *
  46. DENS 0.5;
  47. P1 = 2.5 0. ; P2 = 25. 0. ; Q1 = 2.5 0.5 ;Q2 = 25. 0.5;
  48. PP1 = 5 0. ; QQ1 = 5 0.5 ;
  49. LI1 = P1 D P2 'DINI' 0.25 'DFIN' 2.5 ;
  50. LI2 = P2 D Q2 ;
  51. LI3 = INVE ( Q1 D Q2 'DINI' 0.25 'DFIN' 2.5 );
  52. LI4 = Q1 D P1 ;
  53. LI = LI1 ET LI2 ET LI3 ET LI4 ;
  54. CADRE = DALLER LI1 LI2 LI3 LI4 PLAN ;
  55. SI (NEG GRAPH 'N');
  56. TRAC CADRE ;
  57. FINSI;
  58. *
  59. *=====================================================*
  60. * modele - materiau - caracteristiques en Pa *
  61. * conditions aux limites *
  62. *=====================================================*
  63. *
  64. MO=MODE CADRE POREUX ELASTIQUE PLASTIQUE
  65. DRUCKER_PRAGER ;
  66. *
  67. MAT=MATE MO YOUNG 200.13E6 NU 0.25 MOB 7500E6 COB 1.
  68. PERM 4.E-12 VISCO 4.
  69. ALPH (5.E-5 / 3. ) KF 3000E6 RHOF 0.
  70. ALPF 40E-5 ALPM (19.E-5 / 3.) RHO 0.
  71. ALFA 0. BETA 0.57735027 K 1.E6
  72. GAMM 0. DELT 0.57735027
  73. ETA 0. MU 0.57735027 KL 1.E6 H 0. ;
  74. *
  75. * on bloque les faces superieure et inferieure
  76. *
  77. BL1=BLOQ UZ LI1 ;
  78. BL2=BLOQ UZ LI3 ;
  79. *
  80. BL = BL1 ET BL2 ;
  81. *
  82. *=====================================================*
  83. * 1-ere etape : chargement mecanique *
  84. * Initialisation de la table pour appel a PASAPAS *
  85. *=====================================================*
  86. *
  87. * contraintes initiales
  88. *
  89. ZOZO = MANU CHML MO SMRZ 0.
  90. SMRR -5.E6 SMTT -5.E6 SMZZ -5.E6 MSR0 0.
  91. TYPE CONTRAINTES STRESSES ;
  92. ZUZU = MANU CHML MO GARZ 0.
  93. EPRR 0. EPTT 0. EPZZ 0. PRES 2.5E6
  94. TYPE DEFORMATIONS STRESSES ;
  95. FF = BSIGMA MO ZOZO;
  96. EVT1= EVOL MANU 'TEMPS' (PROG 0 14 )
  97. 'CHARGE' (PROG 1 1 );
  98. CHA1 = CHAR 'MECA' FF EVT1 ;
  99. *
  100. * pression
  101. *
  102. FPI2 = PRESS MASS MO -1.E6 LI4 ;
  103. EVT2= EVOL MANU 'TEMPS' (PROG 0 PAS 1 14 )
  104. 'CHARGE' (PROG 0 PAS 0.25 3.5) ;
  105. CHA2 = CHAR 'MECA' FPI2 EVT2 ;
  106. CHATT = CHA1 ET CHA2 ;
  107. *
  108. TAB1 = TABLE;
  109. TAB1 . 'CONTRAINTES' = TABLE;
  110. TAB1 . 'DEFORMATIONS' = TABLE;
  111. TAB1 . 'DEPLACEMENTS' = TABLE;
  112. TAB1.'CONTRAINTES' . 0 = ZOZO;
  113. TAB1.'DEFORMATIONS'. 0 = ZUZU;
  114. TAB1.'DEPLACEMENTS' . 0 = MANU CHPO CADRE 1 P 2.5E6 ;
  115. TAB1.'TEMPS_CALCULES' = PROG 0 6 10 14;
  116. TAB1.'MODELE' = MO;
  117. TAB1.'CARACTERISTIQUES' = MAT;
  118. TAB1.'BLOCAGES_MECANIQUES' = BL;
  119. TAB1.'CHARGEMENT' = CHATT;
  120. TAB1.'PRECISION'=5.E-7;
  121. *
  122. TMASAU=table;
  123. tab1 . 'MES_SAUVEGARDES'=TMASAU;
  124. TMASAU .'DEFTO'=VRAI;
  125. TMASAU .'DEFIN'=VRAI;
  126. TAB1.'MONO'=monop;tab1.'NB_EL_PARA'=0.9;
  127. PASAPAS TAB1 ;
  128. *
  129. *=====================================================*
  130. * 2-eme etape : temperature *
  131. * Initialisation de la table pour appel a PASAPAS *
  132. *=====================================================*
  133. *
  134. TAB2 = TABLE;
  135. TAB2 . 'VARIABLES_INTERNES' = TABLE;
  136. TAB2 . 'DEPLACEMENTS' = TABLE;
  137. TAB2 . 'CONTRAINTES' = TABLE;
  138. TAB2 . 'DEFORMATIONS' = TABLE;
  139. SS0 = PECHE TAB1 CONTRAINTES ;
  140. EE0 = PECHE TAB1 DEFORMATIONS ;
  141. TAB2.'CONTRAINTES' . 0 = SS0;
  142. TAB2.'DEFORMATIONS' . 0 = EE0;
  143. TAB2.'PRECISION'=5.E-7;
  144. TAB2.'VARIABLES_INTERNES' . 0 =
  145. PECHE TAB1 VARIABLES_INTERNES ;
  146. TAB2.'DEPLACEMENTS' . 0 = PECHE TAB1 DEPLACEMENTS ;
  147. *
  148. FF = BSIGMA MO SS0 ;
  149. EVT1= EVOL MANU 'TEMPS' (PROG 0 100)
  150. 'CHARGE' (PROG 1 1 );
  151. CHA1 = CHAR 'MECA' FF EVT1 ;
  152. *
  153. LEFAC = 1. / ( LOG ( 2.5/25.) ) ;
  154. RAD = COOR 1 CADRE ;
  155. DT = NOMC ( LEFAC * ( LOG (RAD / 25. ) ) ) T ;
  156. *
  157. TABTER = TABLE ;TEMP = TABLE;
  158. TEMP . 0 = 0.;
  159. TEMP . 1 = 100.;
  160. TABTER . 0 = MANU CHPO CADRE 1 T 0. ;
  161. TABTER . 1 = (DT * 100.) ;
  162. CHA2 = CHAR 'T' TEMP TABTER;
  163. *
  164. CHATT = CHA1 ET CHA2;
  165. *
  166. si complet;
  167. tfin = 70.;
  168. sinon;
  169. tfin = 15.;
  170. finsi;
  171. TAB2.'TEMPS_CALCULES' = PROG 0 PAS 5. tfin; ;
  172. TAB2.'MODELE' = MO;
  173. TAB2.'CARACTERISTIQUES' = MAT;
  174. TAB2.'BLOCAGES_MECANIQUES' = BL;
  175. TAB2.'CHARGEMENT' = CHATT;
  176. *
  177. TMASAU=table;
  178. tab2 . 'MES_SAUVEGARDES'=TMASAU;
  179. TMASAU .'DEFTO'=VRAI;
  180. TMASAU .'DEFIN'=VRAI;
  181. TAB2.'PROCESSEURS'= MOT 'AUTOMATIQUE';
  182. PASAPAS TAB2 ;
  183. *
  184. *=====================================================*
  185. * affichage des resultats *
  186. *=====================================================*
  187. *
  188. MESS ' RESULTATS TEST PORE2 ';
  189. MESS ' -------------------- ';
  190. ' On cherche les contraintes et la pression ';
  191. ' Elles sont comparees a des valeurs de reference ';
  192. *
  193. *----------- r{sultats --------------------------------
  194. *
  195. * DISTRIBUTION DES CONTRAINTES ET PRESSION
  196. * POUR T = 70.
  197. *
  198. si complet;
  199. TT = 70. ;
  200. sinon;
  201. TT = 15.;
  202. finsi;
  203. SIGAA = PECHE TAB2 CONTRAINTES TT;
  204. EPSAA = PECHE TAB2 DEFORMATIONS TT;
  205. SIGAB = CHAN 'CHPO' MO SIGAA ;
  206. EPSAB = CHAN 'CHPO' MO EPSAA ;
  207. SIGARR= EVOL CHPO SIGAB SMRR LI1 ;
  208. SIGATT= EVOL CHPO SIGAB SMTT LI1 ;
  209. LAPP = EVOL CHPO EPSAB PRES LI1 ;
  210. SI (NEG GRAPH 'N');
  211. DESS SIGARR ; DESS SIGATT ; DESS LAPP ;
  212. FINSI;
  213. EVSIGAR = EXTR SIGARR ORDO 1;
  214. EVSIGAT = EXTR SIGATT ORDO 1;
  215. EVPP = EXTR LAPP ORDO 1;
  216. *
  217. * VALEURS DE REFERENCE calculees a la date du 1/2/93
  218. *
  219. si complet;
  220. * CONTRAINTE RADIALE
  221. SOLRR = PROG
  222. -1.59459E+06 -1.56051E+06 -1.78383E+06 -1.76850E+06
  223. -1.96931E+06 -1.97473E+06 -2.15711E+06 -2.17942E+06
  224. -2.34668E+06 -2.38270E+06 -2.53753E+06 -2.58467E+06
  225. -2.72927E+06 -2.78543E+06 -2.92159E+06 -2.98507E+06
  226. -3.11427E+06 -3.18371E+06 -3.30718E+06 -3.38150E+06
  227. -3.50030E+06 -3.57863E+06 -3.69361E+06 -3.77540E+06
  228. -3.88744E+06;
  229. SOLRR = SOLRR ET (PROG
  230. -3.97178E+06 -4.08114E+06 -4.16891E+06 -4.29345E+06
  231. -4.34768E+06 -4.48547E+06 -4.54414E+06 -4.62704E+06
  232. -4.68421E+06 -4.74775E+06 -4.79065E+06 -4.83863E+06
  233. -4.86963E+06 -4.90481E+06 -4.92612E+06 -4.95084E+06
  234. -4.96421E+06 -4.98036E+06 -4.98723E+06 -4.99637E+06
  235. -4.99792E+06 -5.00142E+06 );
  236. *
  237. * CONTRAINTE TANGENTIELLE
  238. SOLTT = PROG
  239. -3.59051E+06 -3.55538E+06 -3.77813E+06 -3.76143E+06
  240. -3.96142E+06 -3.96508E+06 -4.14632E+06 -4.16639E+06
  241. -4.33216E+06 -4.36546E+06 -4.51847E+06 -4.56247E+06
  242. -4.70501E+06 -4.75775E+06 -4.89181E+06 -4.95180E+06
  243. -5.07914E+06 -5.14532E+06 -5.26743E+06 -5.33906E+06
  244. -5.45722E+06 -5.53370E+06 -5.64905E+06 -5.72971E+06
  245. -5.84300E+06;
  246. SOLTT = SOLTT ET ( PROG
  247. -5.92761E+06 -6.03950E+06 -6.12667E+06 -6.24105E+06
  248. -6.30245E+06 -6.28413E+06 -6.13481E+06 -6.01293E+06
  249. -5.89978E+06 -5.80418E+06 -5.70545E+06 -5.62098E+06
  250. -5.53670E+06 -5.46384E+06 -5.39075E+06 -5.32680E+06
  251. -5.26302E+06 -5.20651E+06 -5.15031E+06 -5.09985E+06
  252. -5.04987E+06 -5.00444E+06 );
  253. *
  254. * PRESSION
  255. SOLPP = PROG
  256. 1.36747E+06 1.44034E+06 1.51321E+06 1.58674E+06
  257. 1.66027E+06 1.73387E+06 1.80747E+06 1.88113E+06
  258. 1.95479E+06 2.02829E+06 2.10179E+06 2.17490E+06
  259. 2.24800E+06 2.32049E+06 2.39297E+06 2.46467E+06
  260. 2.53636E+06 2.60709E+06 2.67783E+06 2.74773E+06
  261. 2.81762E+06 2.88621E+06 2.95480E+06 3.02386E+06
  262. 3.09291E+06;
  263. SOLPP = SOLPP ET (PROG
  264. 3.15682E+06 3.22073E+06 3.29927E+06 3.37782E+06
  265. 3.40824E+06 3.43867E+06 3.36939E+06 3.30011E+06
  266. 3.24558E+06 3.19105E+06 3.13296E+06 3.07487E+06
  267. 3.01765E+06 2.96042E+06 2.90299E+06 2.84556E+06
  268. 2.78819E+06 2.73081E+06 2.67343E+06 2.61604E+06
  269. 2.55866E+06 2.50128E+06 );
  270. sinon;
  271. SOLRR = prog
  272. -1.56503E+06 -1.57213E+06 -1.75778E+06 -1.77784E+06
  273. -1.94803E+06 -1.98244E+06 -2.14000E+06 -2.18611E+06
  274. -2.33334E+06 -2.38903E+06 -2.52782E+06 -2.59133E+06
  275. -2.72321E+06 -2.79313E+06 -2.91938E+06 -2.99454E+06
  276. -3.11618E+06 -3.19563E+06 -3.31352E+06 -3.39647E+06
  277. -3.51129E+06 -3.59714E+06 -3.70949E+06 -3.79760E+06
  278. -3.90787E+06;
  279. solrr = solrr et (prog
  280. -3.99816E+06 -4.12859E+06 -4.17682E+06 -4.32275E+06
  281. -4.37916E+06 -4.46316E+06 -4.52542E+06 -4.59310E+06
  282. -4.64335E+06 -4.69793E+06 -4.73815E+06 -4.78194E+06
  283. -4.81401E+06 -4.84898E+06 -4.87436E+06 -4.90211E+06
  284. -4.92201E+06 -4.94386E+06 -4.95928E+06 -4.97629E+06
  285. -4.98804E+06 -5.00116E+06);
  286. soltt = prog
  287. -3.56189E+06 -3.56846E+06 -3.75410E+06 -3.77364E+06
  288. -3.94386E+06 -3.97775E+06 -4.13540E+06 -4.18103E+06
  289. -4.32842E+06 -4.38367E+06 -4.52270E+06 -4.58583E+06
  290. -4.71802E+06 -4.78765E+06 -4.91426E+06 -4.98918E+06
  291. -5.11126E+06 -5.19054E+06 -5.30890E+06 -5.39174E+06
  292. -5.50708E+06 -5.59285E+06 -5.70571E+06 -5.79390E+06
  293. -5.90481E+06;
  294. soltt = soltt et ( prog
  295. -5.99480E+06 -6.11303E+06 -6.16665E+06 -6.10728E+06
  296. -6.00175E+06 -5.91312E+06 -5.83171E+06 -5.76581E+06
  297. -5.69639E+06 -5.63972E+06 -5.58251E+06 -5.53566E+06
  298. -5.48786E+06 -5.44850E+06 -5.40854E+06 -5.37545E+06
  299. -5.34188E+06 -5.31389E+06 -5.28556E+06 -5.26174E+06
  300. -5.23769E+06 -5.21735E+06);
  301.  
  302. solpp = prog
  303. 3.52566E+05 4.45884E+05 5.39202E+05 6.32610E+05
  304. 7.26018E+05 8.19101E+05 9.12184E+05 1.00502E+06
  305. 1.09786E+06 1.19042E+06 1.28298E+06 1.37529E+06
  306. 1.46759E+06 1.55963E+06 1.65167E+06 1.74345E+06
  307. 1.83523E+06 1.92692E+06 2.01860E+06 2.10969E+06
  308. 2.20078E+06 2.29305E+06 2.38533E+06 2.47274E+06
  309. 2.56016E+06;
  310. solpp = solpp et ( prog
  311. 2.66334E+06 2.76651E+06 2.80560E+06 2.84469E+06
  312. 2.82168E+06 2.79866E+06 2.78895E+06 2.77924E+06
  313. 2.76631E+06 2.75339E+06 2.74124E+06 2.72909E+06
  314. 2.71676E+06 2.70442E+06 2.69213E+06 2.67984E+06
  315. 2.66755E+06 2.65525E+06 2.64295E+06 2.63066E+06
  316. 2.61836E+06 2.60606E+06);
  317.  
  318. finsi;
  319.  
  320. *
  321. * ecarts
  322. *
  323.  
  324. ERSIGAR =(MAXI ABS ( EVSIGAR - SOLRR ))/
  325. (MAXI ABS SOLRR);
  326. ERSIGAT =(MAXI ABS ( EVSIGAT - SOLTT ))/
  327. (MAXI ABS SOLTT);
  328. ERPP =( MAXI ABS (EVPP - SOLPP))/
  329. (MAXI ABS SOLPP);
  330. *
  331. 'Ecart relatif en Contraintes radiales :' ERSIGAR ;
  332. 'Ecart relatif en Contraintes circonferentielles :' ERSIGAT ;
  333. 'Ecart relatif en Pressions :' ERPP ;
  334. *
  335. *
  336. *=====================================================*
  337. * code fonctionnement *
  338. *=====================================================*
  339. *
  340. * L'ecart maximum entre valeurs de reference et
  341. * calculee doit etre inferieur a 4.d-4
  342. *
  343. SAUT 2 LIGNE;
  344. SI (ERSIGAR <EG 0.001) ;
  345. SI (ERSIGAT <EG 0.001) ;
  346. SI (ERPP <EG 0.001) ;
  347. ERRE 0;
  348. SINON;
  349. ERRE 5;
  350. FINSI;
  351. SINON;
  352. ERRE 5;
  353. FINSI;
  354. SINON;
  355. ERRE 5;
  356. FINSI;
  357. *
  358. *
  359. FIN;
  360.  
  361.  
  362.  
  363.  
  364.  
  365.  
  366.  
  367.  
  368.  
  369.  
  370.  
  371.  
  372.  
  373.  
  374.  
  375.  
  376.  
  377.  
  378.  
  379.  
  380.  
  381.  
  382.  
  383.  
  384.  
  385.  
  386.  

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