Télécharger thm1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : thm1.dgibi
  2. *---------------------------------------------------------------------*
  3. * *
  4. * Cas-test du modele THERMOHYDRIQUE SCHREFLER *
  5. * *
  6. * Description : Simulation d'une structure composite constituee d'un *
  7. * ------------- modele THM et d'un modele THERMIQUE classique en 2D *
  8. * Deformation Plane. De plus, CL de CONVECTION appli- *
  9. * -quee a une extremite du milieu THM. *
  10. * *
  11. * Par defaut, comportement mecanique endommageable pour *
  12. * verifier possibilite couplage avec endommagement. *
  13. * *
  14. * Test de non-regression sur resultats du 7/10/2008. *
  15. * *
  16. * *
  17. *---------------------------------------------------------------------*
  18. *
  19. *
  20.  
  21.  
  22.  
  23. ************tc on enleve ce cas test pour l'instant !
  24. fin;
  25.  
  26.  
  27. *
  28. 'OPTI' 'DIME' 2 'ELEM' 'QUA8' 'MODE' 'PLAN' 'DEFO' ;
  29. * *
  30. * Calcul sur tous les pas de temps : ICOMPLET = VRAI ;
  31. ICOMPLET = FAUX ;
  32. * *
  33. * Simulation avec un modele d'endommagement : IMAZARS = VRAI ;
  34. IMAZARS = VRAI ;
  35. * *
  36. * Traces : ITRAC1 = VRAI ;
  37. ITRAC1 = FAUX ;
  38. * *
  39. 'OPTI' 'ECHO' 0 ;
  40. * *
  41. *---------------------------------------------------------------------*
  42. * *
  43. * I. Geometrie et Donnees Initiales *
  44. * *
  45. *---------------------------------------------------------------------*
  46. * *
  47. *--------------------------- I.1. Geometrie --------------------------*
  48. * *
  49. PA1 = 0. 0. ;
  50. PB1 = 0. 0.01 ;
  51. PC1 = 0.3 0.01 ;
  52. PD1 = 0.3 0. ;
  53. * *
  54. * Nombre d'elements suivant Ox et Oy : *
  55. NBEX1 = 60 ;
  56. NBEY1 = 1 ;
  57. * Droites contour de la geometrie : *
  58. DAB1 = PA1 'DROI' NBEY1 PB1 ;
  59. DCD1 = PC1 'DROI' NBEY1 PD1 ;
  60. DBC1 = PB1 'DROI' NBEX1 PC1 ;
  61. DDA1 = PD1 'DROI' NBEX1 PA1 ;
  62. DFE1 = DCD1 'PLUS' (0.01 0.) ;
  63. * *
  64. SURF1 = 'DALL' DAB1 DBC1 DCD1 DDA1 ;
  65. SURF2 = DCD1 'REGL' 1 DFE1 ;
  66. * *
  67. *-------------------- I.2. Discretisation en Temps -------------------*
  68. * *
  69. 'SI' ICOMPLET ;
  70. LTPSCAL1 = PROG 0. PAS 100. 36000. ;
  71. 'SINO' ;
  72. LTPSCAL1 = PROG 0. PAS 100. 300. ;
  73. 'FINS' ;
  74. * *
  75. *--------------------- I.3. Conditions Initiales ---------------------*
  76. * *
  77. * Temperature (TK0), Pression de Gaz (PG0), Humidite Relative (HR0) *
  78. * et Pression Capillaire (PC0) Initiales : *
  79. TK0 = 298.15 ;
  80. PG0 = 101325. ;
  81. HR0 = 0.60895 ;
  82. PC0 = 68.3596E6 ;
  83. * *
  84. *---------------------------------------------------------------------*
  85. * *
  86. * II. Modelisation *
  87. * *
  88. *---------------------------------------------------------------------*
  89. * *
  90. *---------------- II.1. Modeles et Caracteristiques ------------------*
  91. * *
  92. MODTH1 = 'MODE' SURF1 'THERMOHYDRIQUE' ;
  93. MODTH2 = 'MODE' SURF2 'THERMIQUE' ;
  94. MODCO1 = 'MODE' DAB1 'THERMIQUE' 'CONVECTION' ;
  95. 'SI' IMAZARS ;
  96. MODM1 = 'MODE' SURF1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  97. 'ENDOMMAGEMENT' 'MAZARS' ;
  98. 'SINO' ;
  99. MODM1 = 'MODE' SURF1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  100. 'FINS' ;
  101. MODM2 = 'MODE' SURF2 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' ;
  102. MOD0 = MODTH1 'ET' MODM1 'ET' MODTH2 'ET' MODM2 ;
  103. MOD0 = MODTH1 'ET' MODM1 'ET' MODTH2 'ET' MODM2 'ET' MODCO1 ;
  104. * *
  105. * Module d'Young (MY1), Coefficient de Poisson (NU1) et Dilatation *
  106. * Thermique (ALPH1) : *
  107. MY1 = 'EVOL' 'MANU'
  108. 'T ' ('PROG' 20. 60. 110. 250.)
  109. 'YOUN' ('PROG' 38.5E9 33.E9 29.8E9 20.4E9) ;
  110. NU1 = 0.25 ;
  111. ALPH1 = 11.E-6 ;
  112. * *
  113. * Parametres pour modele MAZARS : *
  114. EPSD1 = 2.E-04 ;
  115. AT1 = 1. ;
  116. BT1 = 1.E4 ;
  117. AC1 = 1.2 ;
  118. BC1 = 1.5E3 ;
  119. BETA1 = 1.06 ;
  120. * *
  121. CH0 = 'MANU' 'CHPO' SURF1 3 'PG' PG0 'PC' PC0 'T' TK0 ;
  122. MATTH1 F0 = @MATETHM MODTH1 CH0 ;
  123. MATTH2 = 'MATE' MODTH2 'RHO' 1. 'C' 0. 'K' 1.E-6 ;
  124. MATCO1 = 'MATE' MODCO1 'H' 1.E6 ;
  125. 'SI' IMAZARS ;
  126. MATM1 = 'MATE' MODM1 'YOUN' MY1 'NU' NU1 'ALPH' ALPH1 'RHO' 2.E3
  127. 'KTR0' EPSD1 'ATRA' AT1 'BTRA' BT1 'ACOM' AC1 'BCOM' BC1 'BETA' BETA1 ;
  128. 'SINO' ;
  129. MATM1 = 'MATE' MODM1 'YOUN' MY1 'NU' NU1 'ALPH' ALPH1 ;
  130. 'FINS' ;
  131. MATM2 = 'MATE' MODM2 'YOUN' 1.E9 'NU' 0.3 'ALPH' 1.E-6 ;
  132. MAT0 = MATTH1 'ET' MATM1 'ET' MATTH2 'ET' MATM2 ;
  133. MAT0 = MATTH1 'ET' MATM1 'ET' MATTH2 'ET' MATM2 'ET' MATCO1 ;
  134. * *
  135. *--------------------- II.2. C.L. et Chargements ---------------------*
  136. * *
  137. * II.2.1. ThermoHydriques : *
  138. CLPG1 = 'BLOQ' 'PG' DAB1 ;
  139. CLPC1 = 'BLOQ' 'PC' DAB1 ;
  140. CLTH0 = CLPG1 'ET' CLPC1 ;
  141. * *
  142. CHPG1 = 'DEPI' CLPG1 1. ;
  143. CHPC1 = 'DEPI' CLPC1 1. ;
  144. CHTK1 = 'MANU' 'CHPO' DAB1 1 'TE ' 1. ;
  145. EVPG1 = 'EVOL' 'MANU' 'TEMPS' ('PROG' 0. 300. 36000.)
  146. ('PROG' 101325. 431325. 758925.) ;
  147. EVPC1 = 'EVOL' 'MANU' 'TEMPS' ('PROG' 0. 300. 36000.)
  148. ('PROG' 6.8362E7 4.39004E7 4.08575E7) ;
  149. EVTK1 = 'EVOL' 'MANU' 'TEMPS' ('PROG' 0. 300. 36000.)
  150. ('PROG' 298.15 418. 440.) ;
  151. CGPG1 = 'CHAR' 'TIMP' CHPG1 EVPG1 ;
  152. CGPC1 = 'CHAR' 'TIMP' CHPC1 EVPC1 ;
  153. CGTK1 = 'CHAR' 'TECO' CHTK1 EVTK1 ;
  154. CG0 = CGPG1 'ET' CGPC1 'ET' CGTK1 ;
  155. * *
  156. * Mecaniques : *
  157. CLME1 = 'BLOQ' 'UX' DAB1 ;
  158. CLME2 = 'BLOQ' 'UY' DBC1 ;
  159. CLME3 = 'BLOQ' 'UY' DDA1 ;
  160. CLME0 = CLME1 'ET' CLME2 'ET' CLME3 ;
  161. * *
  162. *---------------------------------------------------------------------*
  163. * *
  164. * III. Resolution PASAPAS *
  165. * *
  166. *---------------------------------------------------------------------*
  167. * *
  168. *---------------------- III.1. Table de Donnees ----------------------*
  169. * *
  170. TPAS1 = 'TABL' ;
  171. TPAS1 . 'MODELE' = MOD0 ;
  172. TPAS1 . 'CARACTERISTIQUES' = MAT0 ;
  173. TPAS1 . 'BLOCAGES_THERMIQUES' = CLTH0 ;
  174. TPAS1 . 'BLOCAGES_MECANIQUES' = CLME0 ;
  175. TPAS1 . 'CHARGEMENT' = CG0 ;
  176. TPAS1 . 'TEMPS_CALCULES' = LTPSCAL1 ;
  177. TPAS1 . 'PROCEDURE_THERMIQUE' = 'NONLINEAIRE' ;
  178. TPAS1 . 'RELAXATION_THETA' = 1.0 ;
  179. TPAS1 . 'MOVA' = 'D' ;
  180. * *
  181. TPAS1 . 'TEMPERATURES' = 'TABL' ;
  182. TPAS1 . 'TEMPERATURES' . 0 = CH0 ;
  183. * *
  184.  
  185. *------------------------- III.2. Resolution -------------------------*
  186. * *
  187. PASAPAS TPAS1 ;
  188. * *
  189. *---------------------------------------------------------------------*
  190. * *
  191. * IV. Post-traitement *
  192. * *
  193. *---------------------------------------------------------------------*
  194. * *
  195. *------------------- IV.1. Solution de comparaison -------------------*
  196. * *
  197. EVPGR1 = 'EVOL' 'TURQ' 'MANU'
  198. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  199. 'PG' ('PROG' 1.01325E+05 1.01325E+05 1.01325E+05) ;
  200. EVPCR1 = 'EVOL' 'TURQ' 'MANU'
  201. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  202. 'PC' ('PROG' 6.83596E+07 6.83596E+07 6.83596E+07) ;
  203. EVTKR1 = 'EVOL' 'TURQ' 'MANU'
  204. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  205. 'T ' ('PROG' 2.98150E+02 2.98150E+02 2.98150E+02) ;
  206. EVPGR2 = 'EVOL' 'VERT' 'MANU'
  207. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  208. 'PG' ('PROG' 4.31325E+05 1.03086E+05 1.01345E+05) ;
  209. EVPCR2 = 'EVOL' 'VERT' 'MANU'
  210. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  211. 'PC' ('PROG' 4.39004E+07 6.75367E+07 6.83477E+07) ;
  212. EVTKR2 = 'EVOL' 'VERT' 'MANU'
  213. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  214. 'T ' ('PROG' 4.17982E+02 3.01067E+02 2.98192E+02) ;
  215. 'SI' ICOMPLET ;
  216. EVPGR3 = 'EVOL' 'JAUN' 'MANU'
  217. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  218. 'PG' ('PROG' 5.11160E+05 2.70857E+05 1.21151E+05) ;
  219. EVPCR3 = 'EVOL' 'JAUN' 'MANU'
  220. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  221. 'PC' ('PROG' 4.31589E+07 4.31737E+07 5.23776E+07) ;
  222. EVTKR3 = 'EVOL' 'JAUN' 'MANU'
  223. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  224. 'T ' ('PROG' 4.23359E+02 3.83641E+02 3.52063E+02) ;
  225. EVPGR4 = 'EVOL' 'ORAN' 'MANU'
  226. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  227. 'PG' ('PROG' 6.76337E+05 5.09370E+05 3.86782E+05) ;
  228. EVPCR4 = 'EVOL' 'ORAN' 'MANU'
  229. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  230. 'PC' ('PROG' 4.16246E+07 3.94484E+07 4.15511E+07) ;
  231. EVTKR4 = 'EVOL' 'ORAN' 'MANU'
  232. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  233. 'T ' ('PROG' 4.34453E+02 4.08195E+02 3.86481E+02) ;
  234. EVPGR5 = 'EVOL' 'ROUG' 'MANU'
  235. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  236. 'PG' ('PROG' 7.58925E+05 6.08924E+05 5.02148E+05) ;
  237. EVPCR5 = 'EVOL' 'ROUG' 'MANU'
  238. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  239. 'PC' ('PROG' 4.08575E+07 3.79207E+07 3.89898E+07) ;
  240. EVTKR5 = 'EVOL' 'ROUG' 'MANU'
  241. 'X ' ('PROG' 0.00000E+00 5.00000E-02 1.00000E-01)
  242. 'T ' ('PROG' 4.39999E+02 4.16870E+02 3.97931E+02) ;
  243. EVPGR0 = EVPGR1 'ET' EVPGR2 'ET' EVPGR3 'ET' EVPGR4 'ET' EVPGR5 ;
  244. EVPCR0 = EVPCR1 'ET' EVPCR2 'ET' EVPCR3 'ET' EVPCR4 'ET' EVPCR5 ;
  245. EVTKR0 = EVTKR1 'ET' EVTKR2 'ET' EVTKR3 'ET' EVTKR4 'ET' EVTKR5 ;
  246. 'SINO' ;
  247. EVPGR0 = EVPGR1 'ET' EVPGR2 ;
  248. EVPCR0 = EVPCR1 'ET' EVPCR2 ;
  249. EVTKR0 = EVTKR1 'ET' EVTKR2 ;
  250. 'FINS' ;
  251. * *
  252. *------------ IV.2. Depouillement de la solution calculee ------------*
  253. * *
  254. TTEMPS1 = TPAS1 . 'TEMPS' ;
  255. IDIM1 = 'DIME' TTEMPS1 ;
  256. TINDEX1 = 'TABL' ;
  257. 'REPE' BOU1 IDIM1 ;
  258. IM1 = &BOU1 - 1 ;
  259. TPSI1 = TTEMPS1 . IM1 ;
  260. TINDEX1 . TPSI1 = IM1 ;
  261. 'FIN' BOU1 ;
  262. * *
  263. TSOL1 = TPAS1 . 'TEMPERATURES' ;
  264. * *
  265. PX1 = DBC1 'POIN' 'PROC' (0.05 0.) ;
  266. PX2 = DBC1 'POIN' 'PROC' (0.10 0.) ;
  267. DX1 = (PB1 'DROI' 1 PX1) 'ET' (PX1 'DROI' PX2) ;
  268. * *
  269. SOLPG1 = 'EVOL' 'TURQ' 'CHPO' (TSOL1 . (TINDEX1 . 0.)) 'PG' DX1 ;
  270. SOLPC1 = 'EVOL' 'TURQ' 'CHPO' (TSOL1 . (TINDEX1 . 0.)) 'PC' DX1 ;
  271. SOLTK1 = 'EVOL' 'TURQ' 'CHPO' (TSOL1 . (TINDEX1 . 0.)) 'T ' DX1 ;
  272. SOLPG2 = 'EVOL' 'VERT' 'CHPO' (TSOL1 . (TINDEX1 . 300.)) 'PG' DX1 ;
  273. SOLPC2 = 'EVOL' 'VERT' 'CHPO' (TSOL1 . (TINDEX1 . 300.)) 'PC' DX1 ;
  274. SOLTK2 = 'EVOL' 'VERT' 'CHPO' (TSOL1 . (TINDEX1 . 300.)) 'T ' DX1 ;
  275. 'SI' ICOMPLET ;
  276. SOLPG3 = 'EVOL' 'JAUN' 'CHPO' (TSOL1 . (TINDEX1 . 9.E3)) 'PG' DX1 ;
  277. SOLPC3 = 'EVOL' 'JAUN' 'CHPO' (TSOL1 . (TINDEX1 . 9.E3)) 'PC' DX1 ;
  278. SOLTK3 = 'EVOL' 'JAUN' 'CHPO' (TSOL1 . (TINDEX1 . 9.E3)) 'T ' DX1 ;
  279. SOLPG4 = 'EVOL' 'ORAN' 'CHPO' (TSOL1 . (TINDEX1 . 27.E3)) 'PG' DX1 ;
  280. SOLPC4 = 'EVOL' 'ORAN' 'CHPO' (TSOL1 . (TINDEX1 . 27.E3)) 'PC' DX1 ;
  281. SOLTK4 = 'EVOL' 'ORAN' 'CHPO' (TSOL1 . (TINDEX1 . 27.E3)) 'T ' DX1 ;
  282. SOLPG5 = 'EVOL' 'ROUG' 'CHPO' (TSOL1 . (TINDEX1 . 36.E3)) 'PG' DX1 ;
  283. SOLPC5 = 'EVOL' 'ROUG' 'CHPO' (TSOL1 . (TINDEX1 . 36.E3)) 'PC' DX1 ;
  284. SOLTK5 = 'EVOL' 'ROUG' 'CHPO' (TSOL1 . (TINDEX1 . 36.E3)) 'T ' DX1 ;
  285. 'FINS' ;
  286. * *
  287. 'SI' ICOMPLET ;
  288. EVPGTOT1 = SOLPG1 'ET' SOLPG2 'ET' SOLPG3 'ET' SOLPG4 'ET' SOLPG5 ;
  289. EVPCTOT1 = SOLPC1 'ET' SOLPC2 'ET' SOLPC3 'ET' SOLPC4 'ET' SOLPC5 ;
  290. EVTPTOT1 = SOLTK1 'ET' SOLTK2 'ET' SOLTK3 'ET' SOLTK4 'ET' SOLTK5 ;
  291. 'SINO' ;
  292. EVPGTOT1 = SOLPG1 'ET' SOLPG2 ;
  293. EVPCTOT1 = SOLPC1 'ET' SOLPC2 ;
  294. EVTPTOT1 = SOLTK1 'ET' SOLTK2 ;
  295. 'FINS' ;
  296. * *
  297. *--------------------------- IV.3. Traces ----------------------------*
  298. * *
  299. 'SI' ITRAC1 ;
  300. TDESS1 = 'TABL' ;
  301. TDESS1 . 1 = 'MOT' 'MARQ SS CROI NOLI' ;
  302. TDESS1 . 2 = 'MOT' 'MARQ SS CROI NOLI' ;
  303. TDESS1 . 'TITRE' = 'TABL' ;
  304. TDESS1 . 'TITRE' . 1 = 'MOT' 'PAS DE LEGENDE' ;
  305. TDESS1 . 'TITRE' . 2 = 'MOT' 'PAS DE LEGENDE' ;
  306. TDESS1 . 'TITRE' . 6 = 'MOT' 'Initial' ;
  307. TDESS1 . 'TITRE' . 7 = 'MOT' 't = 5 mn' ;
  308. 'SI' ICOMPLET ;
  309. TDESS1 . 3 = 'MOT' 'MARQ SS CROI NOLI' ;
  310. TDESS1 . 4 = 'MOT' 'MARQ SS CROI NOLI' ;
  311. TDESS1 . 5 = 'MOT' 'MARQ SS CROI NOLI' ;
  312. TDESS1 . 'TITRE' . 3 = 'MOT' 'PAS DE LEGENDE' ;
  313. TDESS1 . 'TITRE' . 4 = 'MOT' 'PAS DE LEGENDE' ;
  314. TDESS1 . 'TITRE' . 5 = 'MOT' 'PAS DE LEGENDE' ;
  315. TDESS1 . 'TITRE' . 8 = 'MOT' 't = 2 h 30' ;
  316. TDESS1 . 'TITRE' . 9 = 'MOT' 't = 7 h 30 mn' ;
  317. TDESS1 . 'TITRE' . 10 = 'MOT' 't = 10 h' ;
  318. 'FINS' ;
  319. * *
  320. TIT1 = 'MOT'
  321. 'Simulation (lignes) Vs. Solution de Refence (points).' ;
  322. 'DESS' TDESS1 (EVPGR0 'ET' EVPGTOT1) 'TITR' TIT1 ;
  323. 'DESS' TDESS1 (EVPCR0 'ET' EVPCTOT1) 'TITR' TIT1 ;
  324. 'DESS' TDESS1 (EVTKR0 'ET' EVTPTOT1) 'TITR' TIT1 ;
  325. 'FINS' ;
  326. * *
  327. *------------------------ IV.4. Calcul erreur ------------------------*
  328. * *
  329. ERRPG1 = ('MAXI' 'ABS' ('EXTR' (EVPGR1 - SOLPG1) 'ORDO'))
  330. / ('MAXI' ('EXTR' EVPGR1 'ORDO')) ;
  331. ERRPC1 = ('MAXI' 'ABS' ('EXTR' (EVPCR1 - SOLPC1) 'ORDO'))
  332. / ('MAXI' ('EXTR' EVPCR1 'ORDO')) ;
  333. ERRTK1 = ('MAXI' 'ABS' ('EXTR' (EVTKR1 - SOLTK1) 'ORDO'))
  334. / ('MAXI' ('EXTR' EVTKR1 'ORDO')) ;
  335. ERRPG2 = ('MAXI' 'ABS' ('EXTR' (EVPGR2 - SOLPG2) 'ORDO'))
  336. / ('MAXI' ('EXTR' EVPGR2 'ORDO')) ;
  337. ERRPC2 = ('MAXI' 'ABS' ('EXTR' (EVPCR2 - SOLPC2) 'ORDO'))
  338. / ('MAXI' ('EXTR' EVPCR2 'ORDO')) ;
  339. ERRTK2 = ('MAXI' 'ABS' ('EXTR' (EVTKR2 - SOLTK2) 'ORDO'))
  340. / ('MAXI' ('EXTR' EVTKR2 'ORDO')) ;
  341. 'SI' ICOMPLET ;
  342. ERRPG3 = ('MAXI' 'ABS' ('EXTR' (EVPGR3 - SOLPG3) 'ORDO'))
  343. / ('MAXI' ('EXTR' EVPGR3 'ORDO')) ;
  344. ERRPC3 = ('MAXI' 'ABS' ('EXTR' (EVPCR3 - SOLPC3) 'ORDO'))
  345. / ('MAXI' ('EXTR' EVPCR3 'ORDO')) ;
  346. ERRTK3 = ('MAXI' 'ABS' ('EXTR' (EVTKR3 - SOLTK3) 'ORDO'))
  347. / ('MAXI' ('EXTR' EVTKR3 'ORDO')) ;
  348. ERRPG4 = ('MAXI' 'ABS' ('EXTR' (EVPGR4 - SOLPG4) 'ORDO'))
  349. / ('MAXI' ('EXTR' EVPGR4 'ORDO')) ;
  350. ERRPC4 = ('MAXI' 'ABS' ('EXTR' (EVPCR4 - SOLPC4) 'ORDO'))
  351. / ('MAXI' ('EXTR' EVPCR4 'ORDO')) ;
  352. ERRTK4 = ('MAXI' 'ABS' ('EXTR' (EVTKR4 - SOLTK4) 'ORDO'))
  353. / ('MAXI' ('EXTR' EVTKR4 'ORDO')) ;
  354. ERRPG5 = ('MAXI' 'ABS' ('EXTR' (EVPGR5 - SOLPG5) 'ORDO'))
  355. / ('MAXI' ('EXTR' EVPGR5 'ORDO')) ;
  356. ERRPC5 = ('MAXI' 'ABS' ('EXTR' (EVPCR5 - SOLPC5) 'ORDO'))
  357. / ('MAXI' ('EXTR' EVPCR5 'ORDO')) ;
  358. ERRTK5 = ('MAXI' 'ABS' ('EXTR' (EVTKR5 - SOLTK5) 'ORDO'))
  359. / ('MAXI' ('EXTR' EVTKR5 'ORDO')) ;
  360. 'FINS' ;
  361. * *
  362. 'SI' ICOMPLET ;
  363. ERRPG0 = 'MAXI' ('PROG' ERRPG1 ERRPG2 ERRPG3 ERRPG4 ERRPG5) ;
  364. ERRPC0 = 'MAXI' ('PROG' ERRPC1 ERRPC2 ERRPC3 ERRPC4 ERRPC5) ;
  365. ERRTK0 = 'MAXI' ('PROG' ERRTK1 ERRTK2 ERRTK3 ERRTK4 ERRTK5) ;
  366. 'SINO' ;
  367. ERRPG0 = 'MAXI' ('PROG' ERRPG1 ERRPG2) ;
  368. ERRPC0 = 'MAXI' ('PROG' ERRPC1 ERRPC2) ;
  369. ERRTK0 = 'MAXI' ('PROG' ERRTK1 ERRTK2) ;
  370. 'FINS' ;
  371. ERR0 = 'MAXI' ('PROG' ERRPG0 ERRPC0 ERRTK0) ;
  372. * *
  373. 'SI' (ERR0 '>' 1.E-3) ;
  374. 'ERRE' 5 ;
  375. 'FIN' ;
  376. 'SINO' ;
  377. 'MESS' ' ****** Ecart relatif max. solution de reference =' ERR0 ;
  378. 'SAUT' 1 'LIGN' ;
  379. 'MESS'
  380. '----------------------- CAS-TEST THM1 REUSSI ! -----------------------'
  381. ;
  382. 'SAUT' 1 'LIGN' ;
  383. 'FINS' ;
  384. * *
  385. 'FIN' ;
  386.  
  387.  
  388.  
  389.  
  390.  
  391.  
  392.  

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