Télécharger la plaquette Cast3M


fr en

Télécharger source1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : source1.dgibi
  2. * section : thermique
  3. *----------------------------------------------------------------------*
  4. *                           SOURCE1.DGIBI                              *
  5. *----------------------------------------------------------------------*
  6. *
  7. * Objet :
  8. * -------
  9. *
  10. * Verfication / validation d'un modele de source de chaleur.
  11. * Formulation generale (THERMIQUE SOURCE).
  12. *
  13. * Description :
  14. * -------------
  15. * Comparaison des flux nodaux equivalents obtenus avec le modele a
  16. * ceux obtenus en fournissant directement la valeur de la densite de
  17. * chaleur a l'operateur SOURCE (ancienne syntaxe).
  18. *
  19. * Type de calcul : Aucun
  20. * Mode de calcul : 1D, 2D PLAN, AXIS et 3D
  21. * Type d'element : TRI3 TRI6 QUA4 QUA8 TET4 TE10 PYR5 PY13 PRI6 PR15
  22. *                  CUB8 CU20 BARR TUY2 COQ2 COQ3 COQ4 COQ6 COQ8
  23. * Objectif       : Ecart relatif des flux integres < 1.e-12
  24. *
  25. *----------------------------------------------------------------------*
  26.  
  27. *
  28. *----------------------------------------------------------------------*
  29. *
  30. *------------------------ 1D ELEMENTS LINEAIRES -----------------------*
  31. *
  32. opti dime 1 elem seg2 ;
  33.  
  34. lo1      = 50.e-3 ;
  35. l1       = (poin 0.) droi 10 (poin lo1) ;
  36.  
  37. * Massifs :
  38. * ---------
  39. *
  40. opti mode unid plan ;
  41.  
  42. mo1      = mode l1 thermique ;
  43. ma1      = mate mo1 k 25. ;
  44.  
  45. moq1     = mode l1 thermique source ;
  46. maq1     = mate moq1 'QVOL' 1. ;
  47.  
  48. sq1      = sour moq1 maq1 ;
  49. sqref1   = sour mo1 l1 1. ;
  50.  
  51. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  52. list err1 ;
  53. si (err1 > 1.e-12) ;
  54.   ERRE ' Probleme MODE SOURCE' ;
  55. fins ;
  56.  
  57. opti mode unid axis ;
  58.  
  59. mo1      = mode l1 thermique ;
  60. ma1      = mate mo1 k 25. ;
  61.  
  62. moq1     = mode l1 thermique source ;
  63. maq1     = mate moq1 'QVOL' 1. ;
  64.  
  65. sq1      = sour moq1 maq1 ;
  66. sqref1   = sour mo1 l1 1. ;
  67.  
  68. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  69. list err1 ;
  70. si (err1 > 1.e-12) ;
  71.   ERRE ' Probleme MODE SOURCE' ;
  72. fins ;
  73.  
  74. opti mode unid sphe ;
  75.  
  76. mo1      = mode l1 thermique ;
  77. ma1      = mate mo1 k 25. ;
  78.  
  79. moq1     = mode l1 thermique source ;
  80. maq1     = mate moq1 'QVOL' 1. ;
  81.  
  82. sq1      = sour moq1 maq1 ;
  83. sqref1   = sour mo1 l1 1. ;
  84.  
  85. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  86. list err1 ;
  87. si (err1 > 1.e-12) ;
  88.   ERRE ' Probleme MODE SOURCE' ;
  89. fins ;
  90.  
  91. *---------------------- 1D ELEMENTS QUADRATIQUES ---------------------*
  92. *
  93. opti dime 1 elem seg3 ;
  94.  
  95. lo1      = 50.e-3 ;
  96. l1       = (poin 0.) droi 10 (poin lo1) ;
  97.  
  98. * Massifs :
  99. * ---------
  100. *
  101. opti mode unid plan ;
  102.  
  103. mo1      = mode l1 thermique ;
  104. ma1      = mate mo1 k 25. ;
  105.  
  106. moq1     = mode l1 thermique source ;
  107. maq1     = mate moq1 'QVOL' 1. ;
  108.  
  109. sq1      = sour moq1 maq1 ;
  110. sqref1   = sour mo1 l1 1. ;
  111.  
  112. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  113. list err1 ;
  114. si (err1 > 1.e-12) ;
  115.   ERRE ' Probleme MODE SOURCE' ;
  116. fins ;
  117.  
  118. opti mode unid axis ;
  119.  
  120. mo1      = mode l1 thermique ;
  121. ma1      = mate mo1 k 25. ;
  122.  
  123. moq1     = mode l1 thermique source ;
  124. maq1     = mate moq1 'QVOL' 1. ;
  125.  
  126. sq1      = sour moq1 maq1 ;
  127. sqref1   = sour mo1 l1 1. ;
  128.  
  129. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  130. list err1 ;
  131. si (err1 > 1.e-12) ;
  132.   ERRE ' Probleme MODE SOURCE' ;
  133. fins ;
  134.  
  135. opti mode unid sphe ;
  136.  
  137. mo1      = mode l1 thermique ;
  138. ma1      = mate mo1 k 25. ;
  139.  
  140. moq1     = mode l1 thermique source ;
  141. maq1     = mate moq1 'QVOL' 1. ;
  142.  
  143. sq1      = sour moq1 maq1 ;
  144. sqref1   = sour mo1 l1 1. ;
  145.  
  146. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  147. list err1 ;
  148. si (err1 > 1.e-12) ;
  149.   ERRE ' Probleme MODE SOURCE' ;
  150. fins ;
  151.  
  152. *------------------------ 2D ELEMENTS LINEAIRES -----------------------*
  153. *
  154. opti dime 2 elem qua4 ;
  155.  
  156. ep1      = 10.e-3 ;
  157. l1       = (0 0) droi 10 (0 ep1) ;
  158. s1       = l1    tran 50 (lo1 0) ;
  159. l2       = s1 cote 3 ;
  160. s2       = l2    tran 50 (lo1 ep1) ;
  161. s2       = s2 chan tri3 ;
  162. s0       = s1 et s2 ;
  163.  
  164. *trac s0 titr ' Maillage test' ;
  165.  
  166. opti mode plan ;
  167. *
  168. * Massifs :
  169. * ---------
  170. *
  171. mo1      = mode s0 thermique ;
  172. ma1      = mate mo1 k 25. ;
  173.  
  174. moq1     = mode s0 thermique source ;
  175. maq1     = mate moq1 'QVOL' 1. ;
  176.  
  177. sq1      = sour moq1 maq1 ;
  178. sqref1   = sour mo1 s0 1. ;
  179.  
  180. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  181. list err1 ;
  182. si (err1 > 1.e-12) ;
  183.   ERRE ' Probleme MODE SOURCE' ;
  184. fins ;
  185.  
  186. * Barres :
  187. * --------
  188. *
  189. mo1      = mode l2 thermique barr ;
  190. ma1      = mate mo1 k 1 sect 0.01 ;
  191.  
  192. moq1     = mode l2 thermique source barr ;
  193. maq1     = mate moq1 'QVOL' 1. sect 0.01 ;
  194.  
  195. sq1      = sour moq1 maq1 ;
  196. sqref1   = sour mo1 l2 1. ma1 ;
  197.  
  198. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  199. list err1 ;
  200. si (err1 > 1.e-12) ;
  201.   ERRE ' Probleme MODE SOURCE' ;
  202.   ERRE ' Probleme MODE SOURCE' ;
  203.   ERRE ' Probleme MODE SOURCE' ;
  204. fins ;
  205.  
  206. * Tuyaux :
  207. * --------
  208. *
  209. mo1      = mode l2 thermique tuy2 ;
  210. ma1      = mate mo1 k 1 sect 0.01 ;
  211.  
  212. moq1     = mode l2 thermique source tuy2 ;
  213. maq1     = mate moq1 'QVOL' 1. sect 0.01 ;
  214.  
  215.  
  216. sq1      = sour moq1 maq1 ;
  217. sqref1   = sour mo1 l2 1. ma1 ;
  218.  
  219. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  220. list err1 ;
  221. si (err1 > 1.e-12) ;
  222.   ERRE ' Probleme MODE SOURCE' ;
  223. fins ;
  224.  
  225. * Coques :
  226. * --------
  227. *
  228. mo1      = mode l2 thermique coq2 ;
  229. ma1      = mate mo1 k 1 epai 0.01 ;
  230.  
  231. moq1     = mode l2 thermique source coq2 ;
  232. maq1     = mate moq1 'QVOL' 1. 'QINF' 1. 'QSUP' 1. epai 0.01 ;
  233.  
  234. sq1      = sour moq1 maq1 ;
  235. sqref1   = sour mo1 l2 1. ma1 ;
  236.  
  237. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  238. list err1 ;
  239. si (err1 > 1.e-12) ;
  240.   ERRE ' Probleme MODE SOURCE' ;
  241. fins ;
  242.  
  243. opti mode axis ;
  244. *
  245. * Massifs :
  246. * ---------
  247. *
  248. mo1      = mode s0 thermique ;
  249. ma1      = mate mo1 k 25. ;
  250.  
  251. moq1     = mode s0 thermique source ;
  252. maq1     = mate moq1 'QVOL' 1. ;
  253.  
  254. sq1      = sour moq1 maq1 ;
  255. sqref1   = sour mo1 s0 1. ;
  256.  
  257. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  258. list err1 ;
  259. si (err1 > 1.e-12) ;
  260.   ERRE ' Probleme MODE SOURCE' ;
  261. fins ;
  262.  
  263. * Barres :
  264. * --------
  265. *
  266. mo1      = mode l2 thermique barr ;
  267. ma1      = mate mo1 k 1 sect 0.01 ;
  268.  
  269. moq1     = mode l2 thermique source barr ;
  270. maq1     = mate moq1 'QVOL' 1. sect 0.01 ;
  271.  
  272. sq1      = sour moq1 maq1 ;
  273. sqref1   = sour mo1 l2 1. ma1 ;
  274.  
  275. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  276. list err1 ;
  277. si (err1 > 1.e-12) ;
  278.   ERRE ' Probleme MODE SOURCE' ;
  279. fins ;
  280.  
  281. * Tuyaux :
  282. * --------
  283. *
  284. mo1      = mode l2 thermique tuy2 ;
  285. ma1      = mate mo1 k 1 sect 0.01 ;
  286.  
  287. moq1     = mode l2 thermique source tuy2 ;
  288. maq1     = mate moq1 'QVOL' 1. sect 0.01 ;
  289.  
  290.  
  291. sq1      = sour moq1 maq1 ;
  292. sqref1   = sour mo1 l2 1. ma1 ;
  293.  
  294. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  295. list err1 ;
  296. si (err1 > 1.e-12) ;
  297.   ERRE ' Probleme MODE SOURCE' ;
  298. fins ;
  299.  
  300. * Coques :
  301. * --------
  302. *
  303. mo1      = mode l2 thermique coq2 ;
  304. ma1      = mate mo1 k 1 epai 0.01 ;
  305.  
  306. moq1     = mode l2 thermique source coq2 ;
  307. maq1     = mate moq1 'QVOL' 1. 'QINF' 1. 'QSUP' 1. epai 0.01 ;
  308.  
  309. sq1      = sour moq1 maq1 ;
  310. sqref1   = sour mo1 l2 1. ma1 ;
  311.  
  312. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  313. list err1 ;
  314. si (err1 > 1.e-12) ;
  315.   ERRE ' Probleme MODE SOURCE' ;
  316. fins ;
  317.  
  318. *---------------------- 2D ELEMENTS QUADRATIQUES ----------------------*
  319.  
  320. s0       = s0 chan quad ;
  321. l2       = l2 chan quad ;
  322.  
  323. opti mode plan ;
  324.  
  325. * Massifs :
  326. *----------
  327. *
  328. mo1      = mode s0 thermique ;
  329. ma1      = mate mo1 k 25. ;
  330.  
  331. moq1     = mode s0 thermique source ;
  332.  
  333. maq1     = mate moq1 'QVOL' 1. ;
  334.  
  335. sq1      = sour moq1 maq1 ;
  336. sqref1   = sour mo1 s0 1. ;
  337.  
  338. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  339. list err1 ;
  340. si (err1 > 1.e-12) ;
  341.   ERRE ' Probleme MODE SOURCE' ;
  342. fins ;
  343.  
  344. opti mode axis ;
  345.  
  346. * Massifs :
  347. *----------
  348. *
  349. mo1      = mode s0 thermique ;
  350. ma1      = mate mo1 k 25. ;
  351.  
  352. moq1     = mode s0 thermique source ;
  353.  
  354. maq1     = mate moq1 'QVOL' 1. ;
  355.  
  356. sq1      = sour moq1 maq1 ;
  357. sqref1   = sour mo1 s0 1. ;
  358.  
  359. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  360. list err1 ;
  361. si (err1 > 1.e-12) ;
  362.   ERRE ' Probleme MODE SOURCE' ;
  363. fins ;
  364.  
  365. *------------------------ 3D ELEMENTS LINEAIRES -----------------------*
  366. *
  367. opti dime 3 elem cub8 ;
  368. s0       = s0 chan line ;
  369. l2       = l2 chan line ;
  370.  
  371. v0       = s0 volu tran 5 (0 0 ep1) ;
  372. x1       = v0 coor 1 ;
  373. v3       = v0 elem appu stri (x1 poin supe (1.5 * lo1)) ;
  374. v0       = v0 diff v3 ;
  375. v0       = v0 et (v3 chan tet4) ;
  376.  
  377. el1      = (v0 elem cub8) elem 1 ;
  378. py1      = (el1 chan quaf) chan pyr5 ;
  379. v0       = (v0 diff el1) et py1 ;
  380.  
  381. * Massifs :
  382. *----------
  383. *
  384. mo1      = mode v0 thermique ;
  385. ma1      = mate mo1 k 25. ;
  386.  
  387. moq1     = mode v0 thermique source ;
  388. maq1     = mate moq1 'QVOL' 1. ;
  389.  
  390. sq1      = sour moq1 maq1 ;
  391. sqref1   = sour mo1 v0 1. ;
  392.  
  393. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  394. list err1 ;
  395. si (err1 > 1.e-12) ;
  396.   ERRE ' Probleme MODE SOURCE' ;
  397. fins ;
  398.  
  399. * Barres :
  400. * --------
  401. *
  402. mo1      = mode l2 thermique barr ;
  403. ma1      = mate mo1 k 1 sect 0.01 ;
  404.  
  405. moq1     = mode l2 thermique source barr ;
  406. maq1     = mate moq1 'QVOL' 1. sect 0.01 ;
  407.  
  408.  
  409. sq1      = sour moq1 maq1 ;
  410. sqref1   = sour mo1 l2 1. ma1 ;
  411.  
  412. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  413. list err1 ;
  414. si (err1 > 1.e-12) ;
  415.   ERRE ' Probleme MODE SOURCE' ;
  416. fins ;
  417.  
  418. * Tuyaux :
  419. * --------
  420. *
  421. mo1      = mode l2 thermique tuy2 ;
  422. ma1      = mate mo1 k 1 sect 0.01 ;
  423.  
  424. moq1     = mode l2 thermique source tuy2 ;
  425. maq1     = mate moq1 'QVOL' 1. sect 0.01 ;
  426.  
  427.  
  428. sq1      = sour moq1 maq1 ;
  429. sqref1   = sour mo1 l2 1. ma1 ;
  430.  
  431. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  432. list err1 ;
  433. si (err1 > 1.e-12) ;
  434.   ERRE ' Probleme MODE SOURCE' ;
  435. fins ;
  436.  
  437. * Coques :
  438. * --------
  439. *
  440. mo1      = mode s0 thermique coq3 coq4 ;
  441. ma1      = mate mo1 k 1 epai 0.01 ;
  442.  
  443. moq1     = mode s0 thermique source coq3 coq4 ;
  444. maq1     = mate moq1 'QVOL' 1. 'QINF' 1. 'QSUP' 1. epai 0.01 ;
  445.  
  446. sq1      = sour moq1 maq1 ;
  447. sqref1   = sour mo1 s0 1. ma1 ;
  448.  
  449. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  450. list err1 ;
  451. si (err1 > 1.e-12) ;
  452.   ERRE ' Probleme MODE SOURCE' ;
  453. fins ;
  454.  
  455. *---------------------- 3D ELEMENTS QUADRATIQUES ----------------------*
  456. *
  457. v0       = v0 chan quad ;
  458. s0       = s0 chan quad ;
  459.  
  460. * Massifs :
  461. * ---------
  462. *
  463. mo1      = mode v0 thermique ;
  464. ma1      = mate mo1 k 25. ;
  465.  
  466. moq1     = mode v0 thermique source ;
  467.  
  468. maq1     = mate moq1 'QVOL' 1. ;
  469.  
  470. sq1      = sour moq1 maq1 ;
  471. sqref1   = sour mo1 v0 1. ;
  472.  
  473. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  474. list err1 ;
  475. si (err1 > 1.e-12) ;
  476.   ERRE ' Probleme MODE SOURCE' ;
  477. fins ;
  478.  
  479. * Coques :
  480. * --------
  481. *
  482. mo1      = mode s0 thermique coq6 coq8 ;
  483. ma1      = mate mo1 k 1 epai 0.01 ;
  484.  
  485. moq1     = mode s0 thermique source coq6 coq8 ;
  486. maq1     = mate moq1 'QVOL' 1. 'QINF' 1. 'QSUP' 1. epai 0.01 ;
  487.  
  488. sq1      = sour moq1 maq1 ;
  489. sqref1   = sour mo1 s0 1. ma1 ;
  490.  
  491. err1     = (maxi abs (sq1 - sqref1)) / (maxi abs sqref1) ;
  492. list err1 ;
  493. si (err1 > 1.e-12) ;
  494.   ERRE ' Probleme MODE SOURCE' ;
  495. fins ;
  496.  
  497. *-------------------------------- FIN ---------------------------------*
  498. fin ;
  499.  
  500.  
  501.  
  502.  
  503.  

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