Télécharger la plaquette Cast3M


fr en

Télécharger darcy3_tetraedre_EFMH.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  darcy3_tetraedre_EFMH.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. ************************** CAS TEST : darcy3.dgibi ******************
  6. *
  7.  GRAPH = 'N'   ;
  8. 'SAUT' 'PAGE' ;
  9. *
  10. *-------------------------------------------------------------------
  11. * TEST DARCY3
  12. *                  CALCUL DARCY ORTHOTROPE 3D
  13. *    R�solution par une m�thode d'�l�ments finis mixtes hybrides.
  14. *
  15. * Ce test permet de v�rifier le bon fonctionnement des op�rateurs
  16. * utilis� afin de r�soudre les �quations de DARCY par une m�thode
  17. * d'�l�ments finis mixtes hybrides dans CASTEM2000.
  18. *
  19. * On effectue trois calculs sur un cube, maill� par des cubes
  20. * r�guliers.
  21. * Les conditions aux limites varient suivant le cas consid�r� :
  22. * On impose le flux ou la charge sur les cot�s du domaine.
  23. *
  24. * La solution analytique en charge est un polynome de degr� un,
  25. * la vitesse est constante et la conductivit� hydraulique orthotrope.
  26. *             __           __
  27. *            |               |
  28. *            |  1    0    0  |
  29. *        K = |  0   3/4   0  |
  30. *            |  0    0   1/2 |
  31. *            |__           __|
  32. *
  33. *        H(x,y,z) = -45 x -80 y -60z + 200.
  34. *        V(x,y,z) = ( 45 ; 60 ; 30 )
  35. *
  36. * On s'attend � une pr�cision de l'ordre de la pr�cision machine.
  37. *
  38. *-------------------------------------------------------------------
  39. *
  40. 'SAUT' 'PAGE' ;
  41. *
  42. *- Options g�n�rales de calcul.
  43. *
  44. 'TITR' 'EFMH DARCY ORTHOTROPE 3D Lineaire : darcy3.dgibi' ;
  45. 'OPTI' 'DIME' 3 'ELEM' 'TET4'                             ; 
  46.  'OPTI' 'ECHO' 1                                           ;
  47. *
  48. *                                                      ------------
  49. *                                                      = MAILLAGE =
  50. *                                                      ------------
  51. OEIL  = 5.D0  6.D0  7.D0 ;
  52. VECX  = 1.D0  0.D0  0.D0 ;
  53. VECY  = 0.D0  1.D0  0.D0 ;
  54. VECZ  = 0.D0  0.D0  1.D0 ;
  55. *
  56. ENX = 10          ;
  57. ENY = 10         ;
  58. ENZ = 10          ;
  59. DX  = 1.D0 / ENX ;
  60. DY  = 1.D0 / ENY ;
  61. DZ  = 1.D0 / ENZ ;
  62. *
  63. *- Cr�ation des points
  64. *
  65. A0 = 0.D0  0.D0  0.D0 ;
  66. B0 = 1.D0  0.D0  0.D0 ;
  67. C0 = 1.D0  1.D0  0.D0 ;
  68. D0 = 0.D0  1.D0  0.D0 ;
  69. E0 = 0.D0  0.D0  1.D0 ;
  70. F0 = 1.D0  0.D0  1.D0 ;
  71. G0 = 1.D0  1.D0  1.D0 ;
  72. H0 = 0.D0  1.D0  1.D0 ;
  73. *
  74. *- Cr�ation des droites
  75. *
  76. AB = 'DROI' ENX A0 B0 ;
  77. AD = 'DROI' ENY A0 D0 ;
  78. AE = 'DROI' ENZ A0 E0 ;
  79. BC = 'DROI' ENY B0 C0 ;
  80. BF = 'DROI' ENZ B0 F0 ;
  81. CD = 'DROI' ENX C0 D0 ;
  82. CG = 'DROI' ENZ C0 G0 ;
  83. DH = 'DROI' ENZ D0 H0 ;
  84. EF = 'DROI' ENX E0 F0 ;
  85. EH = 'DROI' ENY E0 H0 ;
  86. FG = 'DROI' ENY F0 G0 ;
  87. GH = 'DROI' ENX G0 H0 ;
  88. *
  89. FE = 'INVE' EF ;
  90. EA = 'INVE' AE ;
  91. DC = 'INVE' CD ;
  92. HD = 'INVE' DH ;
  93. DA = 'INVE' AD ;
  94. HE = 'INVE' EH ;
  95. GF = 'INVE' FG ;
  96. FB = 'INVE' BF ;
  97. *
  98. *- Creation des faces du cube
  99. *
  100. SDRO = 'DALL' AB BF FE EA 'PLAN' ;
  101. SGAU = 'DALL' DC CG GH HD 'PLAN' ;
  102. SBAS = 'DALL' AB BC CD DA 'PLAN' ;
  103. SHAU = 'DALL' EF FG GH HE 'PLAN' ;
  104. SDEV = 'DALL' BC CG GF FB 'PLAN' ;
  105. SDER = 'DALL' AD DH HE EA 'PLAN' ;
  106. *
  107. *- Cr�ation maillage g�om�trique
  108. *
  109. ENXM  = ENX + ENY + ENZ              ;
  110. ELI0  = 1.D0 / ENXM / 10.D0          ;
  111. 'SI' ('EGA' ('VALEUR' 'ELEM') 'CUB8') ;
  112.    CUBE1 = 'PAVE' SDER SDEV SBAS SHAU SGAU SDRO ;
  113. 'SINON' ;
  114.    CUBE1 = 'VOLU' (SDER 'ET' SDEV 'ET' SBAS 'ET' SHAU
  115.                    'ET' SGAU 'ET' SDRO) ;
  116. 'FINSI' ;
  117. QFTOT = CHANGE CUBE1 QUAF ;
  118. QFGAU = CHANGE SGAU QUAF ;
  119. QFDRO = CHANGE SDRO QUAF ;
  120. QFHAU = CHANGE SHAU QUAF ;
  121. QFBAS = CHANGE SBAS QUAF ;
  122. QFDEV = CHANGE SDEV QUAF ;
  123. QFDER = CHANGE SDER QUAF ;
  124.  ELIM ELI0 (QFTOT ET QFGAU ET QFDRO ET QFHAU ET QFBAS ET QFDEV ET
  125.               QFDER )  ;
  126. *
  127. *- Cr�ation maillage HYBRIDE et sous-objets (conditions aux limites)
  128. *
  129. MODHYB = MODE QFTOT 'DARCY' 'ANISOTROPE'     ;
  130. MODGAU = MODE QFGAU 'DARCY' 'ANISOTROPE'     ;
  131. MODDRO = MODE QFDRO 'DARCY' 'ANISOTROPE'     ;
  132. MODHAU = MODE QFHAU 'DARCY' 'ANISOTROPE'     ;
  133. MODBAS = MODE QFBAS 'DARCY' 'ANISOTROPE'     ;
  134. MODDEV = MODE QFDEV 'DARCY' 'ANISOTROPE'     ;
  135. MODDER = MODE QFDER 'DARCY' 'ANISOTROPE'     ;
  136. C11   = 'DOMA' MODHYB 'VOLUME'       ;
  137. CHYB1 = 'DOMA' MODHYB 'SURFACE'       ;
  138. CHYB2 = 'DOMA' MODHYB 'NORMALE'       ;
  139. CEGAU  = 'DOMA' MODGAU 'CENTRE' ;
  140. CEDRO  = 'DOMA' MODDRO 'CENTRE' ;
  141. CEHAU  = 'DOMA' MODHAU 'CENTRE' ;
  142. CEBAS  = 'DOMA' MODBAS 'CENTRE' ;
  143. CEDEV  = 'DOMA' MODDEV 'CENTRE' ;
  144. CEDER  = 'DOMA' MODDER 'CENTRE' ;
  145. *
  146. *- Solution analytique
  147. *
  148. XX  YY  ZZ  = 'COOR' (DOMA MODHYB 'FACE' )  ;
  149. XXC YYC ZZC = 'COOR' (DOMA MODHYB 'CENTRE') ;
  150. *
  151. VKX = 1.D0             ;
  152. VKY = 0.75D0           ;
  153. VKZ = 0.5D0            ;
  154.  
  155. VVKX = MANU 'CHPO'  (DOMA MODHYB 'CENTRE') 'K11' 1.D0 ; 
  156. VVKX = CHANGER ATTRIBUT VVKX  'NATU' DISCRET  ;
  157. VVKY = MANU 'CHPO'  (DOMA MODHYB 'CENTRE') 'K22' 0.75D0 ; 
  158. VVKY = CHANGER ATTRIBUT VVKY  'NATU' DISCRET  ;
  159. VVKZ = MANU 'CHPO'  (DOMA MODHYB 'CENTRE') 'K33' 0.5D0 ;
  160. VVKZ = CHANGER ATTRIBUT VVKZ  'NATU' DISCRET  ;
  161.  
  162. AA  = -45.D0           ;
  163. BB  = -80.D0           ;
  164. CC  = -60.D0           ;
  165. DD  = 200.D0           ;
  166. AAA = -1.D0 * VKX * AA ;
  167. BBB = -1.D0 * VKY * BB ;
  168. CCC = -1.D0 * VKZ * CC ;
  169. *
  170. PANAF = (AA * XX)  + (BB * YY)  + (CC * ZZ)  + DD                 ;
  171. PANAC = (AA * XXC) + (BB * YYC) + (CC * ZZC) + DD                 ;
  172. VANAC = 'MANU' 'CHPO' (DOMA MODHYB 'CENTRE') 3 'VX' AAA 
  173.          'VY' BBB 'VZ' CCC ;
  174. VANAF = 'MANU' 'CHPO' (DOMA MODHYB 'FACE')   3 'VX' AAA 
  175.           'VY' BBB 'VZ' CCC ;
  176. *
  177. *                                                  --------------
  178. *                                                  = RESOLUTION =
  179. *                                                  --------------
  180. *
  181. MATI3  = (NOMC 'K11' VVKX) et 
  182.          (NOMC 'K22' VVKY) et
  183.          (NOMC 'K33' VVKZ) et
  184.          (NOMC 'K21' (0.D0 * VVKX)) et
  185.          (NOMC 'K31' (0.D0 * VVKX)) et      
  186.          (NOMC 'K32' (0.D0 * VVKX))           ;
  187. *                                                        ;
  188. *- Conditions aux limites
  189. *
  190. BBGAU = 'BLOQ' CEGAU 'TH' ;
  191. BBDRO = 'BLOQ' CEDRO 'TH' ;
  192. BBHAU = 'BLOQ' CEHAU 'TH' ;
  193. BBBAS = 'BLOQ' CEBAS 'TH' ;
  194. BBDEV = 'BLOQ' CEDEV 'TH' ;
  195. BBDER = 'BLOQ' CEDER 'TH' ;
  196. *
  197. *- TH impos�e
  198. *
  199. TTIMP = 'REDU' PANAF CEGAU ;
  200. TTIM2 = 'EXCO' TTIMP 'SCAL' 'TH'    ;
  201. EEGAU = 'DEPI' BBGAU TTIM2          ;
  202. TTIMP = 'REDU' PANAF CEDRO ;
  203. TTIM2 = 'EXCO' TTIMP 'SCAL' 'TH'    ;
  204. EEDRO = 'DEPI' BBDRO TTIM2          ;
  205. TTIMP = 'REDU' PANAF CEBAS ;
  206. TTIM2 = 'EXCO' TTIMP 'SCAL' 'TH'    ;
  207. EEBAS = 'DEPI' BBBAS TTIM2          ;
  208. TTIMP = 'REDU' PANAF CEHAU ;
  209. TTIM2 = 'EXCO' TTIMP 'SCAL' 'TH'    ;
  210. EEHAU = 'DEPI' BBHAU TTIM2          ;
  211. TTIMP = 'REDU' PANAF CEDEV ;
  212. TTIM2 = 'EXCO' TTIMP 'SCAL' 'TH'    ;
  213. EEDEV = 'DEPI' BBDEV TTIM2          ;
  214. TTIMP = 'REDU' PANAF CEDER ;
  215. TTIM2 = 'EXCO' TTIMP 'SCAL' 'TH'    ;
  216. EEDER = 'DEPI' BBDER TTIM2          ;
  217. *
  218. *- Flux impos�
  219. *
  220. FLDRO = -60.D0 * (nomc 'FLUX' (doma moddro VOLUME));
  221. FLGAU  = 60.D0  * (nomc 'FLUX' (doma modgau VOLUME));
  222. FLHAU  =  30.D0 * (nomc 'FLUX' (doma modhau VOLUME));
  223. FLBAS  = -30.D0  * (nomc 'FLUX' (doma modbas VOLUME));
  224. FLDEV  =  45.D0 * (nomc 'FLUX' (doma moddev VOLUME)); 
  225. FLDER  = -45.D0  * (nomc 'FLUX' (doma modder VOLUME));
  226. *
  227. *
  228. *- Assemblage et r�solution en TH
  229. *
  230.  
  231. h_lim = 'RESOUD' ( BBHAU 'ET' BBBAS 'ET' BBDRO 'ET'
  232.                   BBGAU 'ET' BBDEV 'ET' BBDER)
  233.                  (EEHAU 'ET' EEBAS 'ET' EEDRO 'ET'
  234.                   EEGAU 'ET' EEDEV 'ET' EEDER) ;
  235. h_lim = 'EXCO' 'TH' h_lim;
  236. CHCLIM = TABLE;
  237. CHCLIM . 'DIRICHLET' = 'NOMC' 'I35' h_lim;
  238.  
  239.  
  240. GEOL1 = TABLE;
  241. GEOL1 . 'CONCENTRATION'     = 'NOMC' 'I35' (0.D0 * PANAC)  ;
  242. GEOL1 . 'LUMP'              = FAUX                  ;
  243. GEOL1 . 'TYPDISCRETISATION' = 'EFMH'                ;
  244. GEOL1 . 'THETA_DIFFUSION'   = 1.0D0                 ;
  245. GEOL1 . 'THETA_CONVECTION'  = 1.0D0                 ;
  246. GEOL1 . 'DECENTREMENT'      = FAUX                  ;
  247. GEOL1 . 'DELTAT'            = 1.D15                 ;
  248. GEOL1 . 'DIFFUSIVITE'       = mati3                 ;
  249. GEOL1 . 'SOLVEUR'           = 2                     ;
  250. GEOL1 . 'PRECONDITIONNEUR'  = 3                     ;
  251. GEOL1 . 'POROSITE'          = 0.D0 * PANAC          ;
  252. GEOL1 . 'CLIMITES'          = CHCLIM                ;
  253. GEOL1 . 'RECALCUL'          = VRAI                  ;
  254. *
  255. GEOL1 GEOL2 = TRANGEOL Modhyb GEOL1;
  256. CHTER1 = 'NOMC' 'TH' GEOL2 . 'TRACE_CONC';
  257. PCEN1 = 'NOMC' 'H' GEOL1 . 'CONCENTRATION' ;
  258. QFACE1 = 'NOMC' 'FLUX' GEOL1 . 'FLUXDIFF' ;
  259.  
  260.  
  261. * deuxieme probleme
  262.  
  263.  
  264.  
  265. h_lim = 'RESOUD' (BBHAU 'ET' BBBAS 'ET' BBGAU 'ET' BBDRO)
  266.                  (EEHAU 'ET' EEBAS 'ET' EEGAU 'ET' EEDRO) ;
  267. h_lim = 'EXCO' 'TH' h_lim;
  268. CHCLIM = TABLE;
  269. CHCLIM . 'NEUMANN' = ('NOMC' 'I35' (FLDEV 'ET' FLDER));
  270. CHCLIM . 'DIRICHLET' = 'NOMC' 'I35' h_lim;
  271.  
  272.  
  273. GEOL1 = TABLE;
  274. GEOL1 . 'CONCENTRATION'     = 'NOMC' 'I35' (0.D0 * PANAC)  ;
  275. GEOL1 . 'LUMP'              = FAUX                  ;
  276. GEOL1 . 'TYPDISCRETISATION' = 'EFMH'                ;
  277. GEOL1 . 'THETA_DIFFUSION'   = 1.0D0                 ;
  278. GEOL1 . 'THETA_CONVECTION'  = 1.0D0                 ;
  279. GEOL1 . 'DECENTREMENT'      = FAUX                  ;
  280. GEOL1 . 'DELTAT'            = 1.D15                 ;
  281. GEOL1 . 'DIFFUSIVITE'       = mati3                 ;
  282. GEOL1 . 'SOLVEUR'           = 2                     ;
  283. GEOL1 . 'PRECONDITIONNEUR'  = 3                     ;
  284. GEOL1 . 'POROSITE'          = 0.D0 * PANAC          ;
  285. GEOL1 . 'CLIMITES'          = CHCLIM                ;
  286. GEOL1 . 'RECALCUL'          = VRAI                  ;
  287. *
  288. GEOL1 GEOL2 = TRANGEOL Modhyb GEOL1;
  289. CHTER2 = 'NOMC' 'TH' GEOL2 . 'TRACE_CONC';
  290. PCEN2 = 'NOMC' 'H' GEOL1 . 'CONCENTRATION' ;
  291. QFACE2 = 'NOMC' 'FLUX' GEOL1 . 'FLUXDIFF' ;
  292.  
  293.  
  294. * troisieme probleme
  295.  
  296. h_lim = 'RESOUD' (BBGAU 'ET' BBDEV)
  297.                  (EEGAU 'ET' EEDEV) ;
  298. h_lim = 'EXCO' 'TH' h_lim;
  299. CHCLIM = TABLE;
  300. CHCLIM . 'NEUMANN' = ('NOMC' 'I35' ( FLHAU 'ET' FLBAS
  301.                                  'ET' FLDRO 'ET' FLDER));
  302. CHCLIM . 'DIRICHLET' = 'NOMC' 'I35' h_lim;
  303.  
  304.  
  305. GEOL1 = TABLE;
  306. GEOL1 . 'CONCENTRATION'     = 'NOMC' 'I35' (0.D0 * PANAC)  ;
  307. GEOL1 . 'LUMP'              = FAUX                  ;
  308. GEOL1 . 'TYPDISCRETISATION' = 'EFMH'                ;
  309. GEOL1 . 'THETA_DIFFUSION'   = 1.0D0                 ;
  310. GEOL1 . 'THETA_CONVECTION'  = 1.0D0                 ;
  311. GEOL1 . 'DECENTREMENT'      = FAUX                  ;
  312. GEOL1 . 'DELTAT'            = 1.D15                 ;
  313. GEOL1 . 'DIFFUSIVITE'       = mati3                 ;
  314. GEOL1 . 'SOLVEUR'           = 2                     ;
  315. GEOL1 . 'PRECONDITIONNEUR'  = 3                     ;
  316. GEOL1 . 'POROSITE'          = 0.D0 * PANAC          ;
  317. GEOL1 . 'CLIMITES'          = CHCLIM                ;
  318. GEOL1 . 'RECALCUL'          = VRAI                  ;
  319. *
  320. GEOL1 GEOL2 = TRANGEOL Modhyb GEOL1;
  321. GEOL1 GEOL2 = TRANGEOL Modhyb GEOL1 GEOL2;
  322.  
  323. CHTER3 = 'NOMC' 'TH' GEOL2 . 'TRACE_CONC';
  324. PCEN3 = 'NOMC' 'H' GEOL1 . 'CONCENTRATION' ;
  325. QFACE3 = 'NOMC' 'FLUX' GEOL1 . 'FLUXDIFF' ;
  326.  
  327.  
  328.  
  329. *
  330. *- Calcul de V
  331. *
  332. VCENT1 = 'HVIT' MODHYB  QFACE1         ;
  333. QFACE1 = 'EXCO' QFACE1 'FLUX' 'SCAL'             ;
  334. VFACE1 = QFACE1 * CHYB2 / CHYB1                  ;
  335. *
  336. VCENT2 = 'HVIT' MODHYB  QFACE2         ;
  337. QFACE2 = 'EXCO' QFACE2 'FLUX' 'SCAL'             ;
  338. VFACE2 = QFACE2 * CHYB2 / CHYB1                  ;
  339. *
  340. VCENT3 = 'HVIT' MODHYB  QFACE3         ;
  341. QFACE3 = 'EXCO' QFACE3 'FLUX' 'SCAL'             ;
  342. VFACE3 = QFACE3 * CHYB2 / CHYB1                  ;
  343. *
  344. *                                                   -----------------
  345. *                                                   = Calcul ERREUR =
  346. *                                                   -----------------
  347. * ERReur relative en Trace de charge TH aux faces des �l�ments
  348. * ERReur relative en charge H au centre des �l�ments
  349. * Erreur relative sur la vitesse au centre des �l�ments
  350. *
  351. ERRTP1 = 'EXCO' CHTER1 'TH' 'SCAL' ;
  352. ERRTP1 = (ERRTP1 - PANAF) ** 2     ; 
  353. ERRTP1 = 'RESULT' ERRTP1           ;
  354. ERRTP1 = ERRTP1 '/' ('RESULT' (PANAF ** 2));
  355. ERRTP1 = ('MAXIMUM' errtp1)**0.5D0 ;
  356.  
  357. ERRP1  = 'EXCO' PCEN1 'H' 'SCAL'   ;
  358. ERRP1  = (ERRP1 - PANAC) ** 2      ;
  359. ERRP1  = 'RESULT' ERRP1            ;
  360. ERRP1  = ERRP1 '/' ('RESULT' (PANAC ** 2));
  361. ERRP1  = ('MAXIMUM' errp1)**0.5D0  ;
  362. *
  363. ERRTP2 = 'EXCO' CHTER2 'TH' 'SCAL' ;
  364. ERRTP2 = (ERRTP2 - PANAF) ** 2     ; 
  365. ERRTP2 = 'RESULT' ERRTP2           ;
  366. ERRTP2 = ERRTP2 '/' ('RESULT' (PANAF ** 2));
  367. ERRTP2 = ('MAXIMUM' errtp2)**0.5D0 ;
  368.  
  369. ERRP2  = 'EXCO' PCEN2 'H' 'SCAL'   ;
  370. ERRP2  = (ERRP2 - PANAC) ** 2      ;
  371. ERRP2  = 'RESULT' ERRP2            ;
  372. ERRP2  = ERRP2 '/' ('RESULT' (PANAC ** 2));
  373. ERRP2  = ('MAXIMUM' errp2)**0.5D0  ;
  374. *
  375. *
  376. ERRTP3 = 'EXCO' CHTER3 'TH' 'SCAL' ;
  377. ERRTP3 = (ERRTP3 - PANAF) ** 2     ; 
  378. ERRTP3 = 'RESULT' ERRTP3           ;
  379. ERRTP3 = ERRTP3 '/' ('RESULT' (PANAF ** 2));
  380. ERRTP3 = ('MAXIMUM' errtp3)**0.5D0 ;
  381.  
  382.  
  383. ERRP3  = 'EXCO' PCEN3 'H' 'SCAL'   ;
  384. ERRP3  = (ERRP3 - PANAC) ** 2      ;
  385. ERRP3  = 'RESULT' ERRP3            ;
  386. ERRP3  = ERRP3 '/' ('RESULT' (PANAC ** 2));
  387. ERRP3  = ('MAXIMUM' errp3)**0.5D0  ;
  388. *
  389. *
  390. MOT1 = 'MOTS' 'VX' 'VY' 'VZ'        ;
  391. VDVD = 'PSCA' VANAC VANAC MOT1 MOT1 ;
  392. *
  393. VD1  = VANAC - VCENT1            ;
  394. VC1  = 'PSCA' VD1 VD1 MOT1 MOT1  ;
  395. SDC1 = 'MAXIMUM' ( ('RESULT'  VC1) / ('RESULT' VDVD )) ;
  396. SDC1 = SDC1 '**' 0.5             ;
  397. VD2  = VANAC - VCENT2            ;
  398. VC2  = 'PSCA' VD2 VD2 MOT1 MOT1  ;
  399. SDC2 = 'MAXIMUM' ( ('RESULT'  VC2) / ('RESULT' VDVD) ) ;
  400. SDC2 = SDC2 '**' 0.5             ;
  401. VD3  = VANAC - VCENT3            ;
  402. VC3  = 'PSCA' VD3 VD3 MOT1 MOT1  ;
  403. SDC3 = 'MAXIMUM' ( ('RESULT'  VC3) / ('RESULT' VDVD) ) ;
  404. SDC3 = SDC3 '**' 0.5             ;
  405. *
  406. *                                                 -------------------
  407. *                                                 = Trac� resultats =
  408. *                                                 -------------------
  409. 'SI' ('NEG' GRAPH 'N') ;
  410. *
  411. *- Transformation des quantit�s aux centres en MCHAML constant.
  412. *
  413. ERRP1 = 'KCHA' MODHYB 'CHAM' ERRP1 ;
  414. ERRP2 = 'KCHA' MODHYB 'CHAM' ERRP2 ;
  415. ERRP3 = 'KCHA' MODHYB 'CHAM' ERRP3 ;
  416. SDS1  = 'KCHA' MODHYB 'CHAM' SDC1  ;
  417. SDS2  = 'KCHA' MODHYB 'CHAM' SDC2  ;
  418. SDS3  = 'KCHA' MODHYB 'CHAM' SDC3  ;
  419. *
  420. * Dans chaque cas on trace
  421. *   L'erreur relative sur la charge au centre
  422. *   L'erreur relative sur la Vitesse au centre
  423. *
  424. 'TITR' 'darcy3/1 : Erreur relative sur la charge'  ;
  425. 'TRAC' MODHYB ERRP1                                ;
  426. 'TITR' 'darcy3/1 : Erreur relative sur la vitesse' ;
  427. 'TRAC' MODHYB SDS1                                 ;
  428. *
  429. 'TITR' 'darcy3/2 : Erreur relative sur la charge'  ;
  430. 'TRAC' MODHYB ERRP2                                ;
  431. 'TITR' 'darcy3/2 : Erreur relative sur la vitesse' ;
  432. 'TRAC' MODHYB SDS2                                 ;
  433. *
  434. 'TITR' 'darcy3/3 : Erreur relative sur la charge'  ;
  435. 'TRAC' MODHYB ERRP3                                ;
  436. 'TITR' 'darcy3/3 : Erreur relative sur la vitesse' ;
  437. 'TRAC' MODHYB SDS3                                 ;
  438. *
  439. 'FINSI' ;
  440. *
  441. *                                                 -------------------
  442. *                                                 = Gestion ERREURS =
  443. *                                                 -------------------
  444. MAXTP1 = ERRTP1 ;
  445. MAXTP2 = ERRTP2 ;
  446. MAXTP3 = ERRTP3 ;
  447. MAXP1  = ERRP1  ;
  448. MAXP2  = ERRP2  ;
  449. MAXP3  = ERRP3  ;
  450. MAXV1  = SDC1   ;
  451. MAXV2  = SDC2   ;
  452. MAXV3  = SDC3   ;
  453. *
  454. 'SAUT' 'PAGE'                                                       ;
  455. 'SAUT' 2 'LIGNE'                                                    ;
  456. 'MESS' '                             ERREURS RELATIVES            ' ;
  457. 'SAUT' 1 'LIGNE'                                                    ;
  458. 'MESS' ' cas test          TH                  H                 V' ;
  459. 'SAUT' 1 'LIGNE'                                                    ;
  460. 'MESS' ' numero1     ' maxtp1  '   ' maxp1 '   ' maxv1              ;
  461. 'SAUT' 1 'LIGNE'                                                    ;
  462. 'MESS' ' numero2     ' maxtp2  '   ' maxp2 '   ' maxv2              ;
  463. 'SAUT' 1 'LIGNE'                                                    ;
  464. 'MESS' ' numero3     ' maxtp3  '   ' maxp3 '   ' maxv3              ;
  465. 'SAUT' 2 'LIGNE'                                                    ;
  466. *
  467. EPS0 = 1.D-8 ;
  468. LOG1 = MAXTP1 > EPS0 ; LOG2 = MAXTP2 > EPS0 ; LOG3 = MAXTP3 > EPS0 ;
  469. LOG4 = MAXP1  > EPS0 ; LOG5 = MAXP2  > EPS0 ; LOG6 = MAXP3  > EPS0 ;
  470. LOG7 = MAXV1  > EPS0 ; LOG8 = MAXV2  > EPS0 ; LOG9 = MAXV3  > EPS0 ;
  471. LTP0 = LOG1 'OU' LOG2 'OU' LOG3 ;
  472. LP0  = LOG4 'OU' LOG5 'OU' LOG6 ;
  473. LV0  = LOG7 'OU' LOG8 'OU' LOG9 ;
  474. L0   = LTP0 'OU' LP0  'OU' LV0  ;
  475. 'SI' ( L0 ) ;
  476.    'ERRE' 5 ;
  477. 'SINO'      ;
  478.    'ERRE' 0 ;
  479. 'FINSI'     ;
  480. *
  481. 'FIN' ;
  482.  
  483.  
  484.  
  485.  
  486.  
  487.  
  488.  
  489.  
  490.  
  491.  
  492.  
  493.  
  494.  
  495.  
  496.  
  497.  
  498.  
  499.  
  500.  
  501.  
  502.  

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