Télécharger umat05.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : umat05.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *=======================================================================
  5. *
  6. * CAS TEST DE VALIDATION DE LA PRISE EN COMPTE D'UNE LOI DE
  7. *
  8. * COMPORTEMENT MECANIQUE NON LINEAIRE EXTERNE
  9. *
  10. *=======================================================================
  11. *
  12. * Modele 'NON_LINEAIRE' 'UTILISATEUR', integrateur specifique UMAT
  13. *
  14. * Test No5 : Modele PLASTIQUE ISOTROPE
  15. * Reprise du test isotro_cis.dgibi
  16. *
  17. *=======================================================================
  18.  
  19.  
  20. opti echo 0 ;
  21. * ajout de option epsilon lineaire pour la precision des test!
  22. OPTION epsilon lineaire;
  23.  
  24.  
  25. *=======================================================================
  26. *
  27. * PROCEDURES DEPENDANTES
  28. *
  29. *=======================================================================
  30.  
  31.  
  32.  
  33.  
  34. *=======================================================================
  35. *
  36. * JEU DE DONNEES PRINCIPAL
  37. *
  38. *=======================================================================
  39.  
  40. *=======================================================================
  41. * Resolution avec le modele 'PLASTIQUE' 'ISOTROPE' de CAST3M
  42. * Reprise du jeu de donnees isotro_cis.dgibi
  43. * Element CUB8
  44.  
  45. **************************************************
  46. * *
  47. * TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT *
  48. * -------------------------------------------- *
  49. * DE MATERIAU *
  50. * ----------- *
  51. * *
  52. * LOI DE COMPORTEMENT ELASTO-PLASTIQUE *
  53. * ISOTROPE *
  54. * *
  55. * COMPARAISON DE CALCULS SUR DES ELEMENTS: *
  56. * - COQUE EPAISSE ( MFR=5 ) *
  57. * - COQUE MINCE ( MFR=3 ) *
  58. * - MASSIF 3D ( MFR=1 ) *
  59. * *
  60. * STRUCTURE : UNE PLAQUE CARREE *
  61. * - COTE : .1 M *
  62. * - EPAISSEUR : 5.E-4 M *
  63. * *
  64. * CHARGEMENT : DEPLACEMENTS IMPOSES *
  65. * CISAILLEMENT PUR *
  66. * *
  67. **************************************************
  68. opti echo 0 dime 3 elem cub8 ;
  69. *
  70. *
  71. * Maillage
  72. *
  73. rayon0 = .05 ;
  74. epais0 = .0005 ;
  75. l1 = epais0 + rayon0 ;
  76. l2 = .1 ;
  77. n1 = 1 ;
  78. n2 = 1 ;
  79. p1 = rayon0 0. 0. ;
  80. p2 = l1 0. 0. ;
  81. p3 = l1 l2 0. ;
  82. p4 = rayon0 l2 0. ;
  83. p11 = rayon0 0. l2 ;
  84. p21 = l1 0. l2 ;
  85. p31 = l1 l2 l2 ;
  86. p41 = rayon0 l2 l2 ;
  87. d1 = d p1 n1 p2 ;
  88. d2 = d p2 n2 p3 ;
  89. d3 = d p3 n1 p4 ;
  90. d4 = d p4 n2 p1 ;
  91. d11 = d p11 n1 p21 ;
  92. d21 = d p21 n2 p31 ;
  93. d31 = d p31 n1 p41 ;
  94. d41 = d p41 n2 p11 ;
  95. *p01 = ( rayon0 + ( epais0 / 2. ) ) 0. l2 ;
  96. *p02 = ( rayon0 + ( epais0 / 2. ) ) l2 l2 ;
  97. sur1 = 'DALL' d1 d2 d3 d4 ;
  98. sur2 = 'DALL' d11 d21 d31 d41 ;
  99. vol1 = sur1 'VOLU' n2 sur2 ;
  100. *titr ' maillage ' ;
  101. *trace vol1 ;
  102. *
  103. * Deplacement maximal
  104. *
  105. depmax = .01 ;
  106. *
  107. * Courbe d'ecrouissage
  108. *
  109. beta0 = 1250. ;
  110. sig0 = 410.E6 ;
  111. xm = .5 ;
  112. youn0 = 7.34E10 ;
  113. epsp0 = 'PROG' 0. ;
  114. deb0 = 'PROG' 0. ;
  115. epsmax0 = depmax / l2 ;
  116. pas0 = epsmax0 / 50 ;
  117. epsp1 = 'PROG' 0. 'PAS' pas0 epsmax0 ;
  118. dim00 = 'DIME' epsp1 ;
  119. un = 'PROG' dim00 '*' 1. ;
  120. Y0 = sig0 ;
  121. epsp2 = beta0 * epsp1 ;
  122. epsp3 = epsp2 + un ;
  123. epsp31 = log epsp3 ;
  124. epsp32 = epsp31 * xm ;
  125. epsp4 = exp epsp32 ;
  126. Y1 = Y0 * epsp4 ;
  127. *
  128. epsp1 = epsp0 'ET' epsp1 ;
  129. Y1 = deb0 'ET' Y1 ;
  130. epsp2 = epsp1 + ( Y1 / youn0 ) ;
  131. trac00 = 'EVOL' 'MANU' 'defo' epsp1 'sigm' Y1 ;
  132. *
  133. * Modele de calcul
  134. *
  135. trac0 = 'EVOL' 'MANU' 'defo' epsp2 'sigm' Y1 ;
  136. *'DESS' trac0 ;
  137. mod0 = 'MODE' vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  138. 'PLASTIQUE' 'ISOTROPE' ;
  139. mat0 = 'MATE' mod0 'YOUN' youn0 'NU' 0.33
  140. 'TRAC' trac0 ;
  141. *
  142. * Conditions aux limites
  143. *
  144. cl1 = 'BLOQ' sur1 'UX' 'RY' ;
  145. cl2 = 'BLOQ' sur2 'UX' 'RY' ;
  146. cl3 = 'BLOQ' vol1 'UY' 'UZ' ;
  147. cl0 = cl1 'ET' cl2 'ET' cl3 ;
  148. *
  149. * Chargement
  150. *
  151. chp1 = 'COOR' 3 vol1 ;
  152. chp1 = chp1 * depmax / l2 ;
  153. chp1 = 'NOMC' chp1 'UY' ;
  154. chp2 = 'MANU' 'CHPO' vol1 1 'UZ' 0. ;
  155. *trace chp1 vol1 ;
  156. chp0 = chp1 + chp2 ;
  157. dep1 = 'DEPI' cl3 chp0 ;
  158. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1000. ) y ( 'PROG' 0. 1. ) ;
  159. cha0 = 'CHAR' 'DIMP' dep1 ev0 ;
  160. *vec1 = 'VECT' dep1 1.E-5 'FX' 'FY' 'FZ' 'ROUG' ;
  161. *titr ' chargement ' ;
  162. *trace vec1 vol1 ;
  163. *
  164. * Temps du calcul
  165. *
  166. dt0 = 50. ;
  167. tfin0 = 1000. ;
  168. *
  169. *
  170. * Resolution par PASAPAS
  171. *---------------------------------------
  172. *
  173. ta1 = 'TABLE' ;
  174. ta1 .'MODELE' = mod0 ;
  175. ta1 .'CARACTERISTIQUES' = mat0 ;
  176. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  177. ta1 .'CHARGEMENT' = cha0 ;
  178. *ta1 .'TEMPS_CALCULES' = ( ( 'PROG' 0. 'PAS' dt0 80. ) 'ET'
  179. *( 'PROG' 85. 'PAS' 5. 90. ) 'ET'
  180. *( 'PROG' 92.5 'PAS' 2.5 tfin0 ) ) ;
  181. ta1 .'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ;
  182. *
  183. * MODIF
  184. * MODIF_________________________________________________________________
  185. TEMPS 'ZERO' ;
  186. * MODIF_________________________________________________________________
  187. * MODIF
  188. TMASAU=table;
  189. ta1 . 'MES_SAUVEGARDES'=TMASAU;
  190. TMASAU .'DEFTO'=VRAI;
  191. TMASAU .'DEFIN'=VRAI;
  192. PASAPAS ta1 ;
  193. * MODIF
  194. * MODIF_________________________________________________________________
  195. TABTPS = TEMP 'NOEC';
  196. CPUint = TABTPS.'TEMPS_CPU'.'INITIAL';
  197. * MODIF_________________________________________________________________
  198. * MODIF
  199.  
  200.  
  201. *=======================================================================
  202. * Resolution avec le modele 'NON_LINEAIRE' 'UTILISATEUR'
  203.  
  204. * Definition du modele et affectation des proprietes de materiau--------
  205.  
  206. LCM05_0 = MOTS 'YOUN' 'NU ' ;
  207.  
  208. LCM05_1 = MOTS 'E001' 'S001' 'E002' 'S002' 'E003' 'S003' 'E004' 'S004'
  209. 'E005' 'S005' 'E006' 'S006' 'E007' 'S007' 'E008' 'S008'
  210. 'E009' 'S009' 'E010' 'S010' ;
  211.  
  212. LCM05_2 = MOTS 'E011' 'S011' 'E012' 'S012' 'E013' 'S013' 'E014' 'S014'
  213. 'E015' 'S015' 'E016' 'S016' 'E017' 'S017' 'E018' 'S018'
  214. 'E019' 'S019' 'E020' 'S020' ;
  215.  
  216. LCM05_3 = MOTS 'E021' 'S021' 'E022' 'S022' 'E023' 'S023' 'E024' 'S024'
  217. 'E025' 'S025' 'E026' 'S026' 'E027' 'S027' 'E028' 'S028'
  218. 'E029' 'S029' 'E030' 'S030' ;
  219.  
  220. LCM05_4 = MOTS 'E031' 'S031' 'E032' 'S032' 'E033' 'S033' 'E034' 'S034'
  221. 'E035' 'S035' 'E036' 'S036' 'E037' 'S037' 'E038' 'S038'
  222. 'E039' 'S039' 'E040' 'S040' ;
  223.  
  224. LCM05_5 = MOTS 'E041' 'S041' 'E042' 'S042' 'E043' 'S043' 'E044' 'S044'
  225. 'E045' 'S045' 'E046' 'S046' 'E047' 'S047' 'E048' 'S048'
  226. 'E049' 'S049' 'E050' 'S050' 'E051' 'S051' 'E052' 'S052' ;
  227.  
  228. LCMAT05 = LCM05_0 et LCM05_1 et LCM05_2 et LCM05_3 et
  229. LCM05_4 et LCM05_5 ;
  230.  
  231. LCVAR05 = MOTS 'EPSE' 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  232.  
  233. mo_util = MODE vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  234. 'NON_LINEAIRE' 'UTILISATEUR' 'NUME_LOI' 5
  235. 'C_MATERIAU' LCMAT05 'C_VARINTER' LCVAR05 ;
  236.  
  237. ma_u0 = MATE mo_util 'YOUN' youn0 'NU ' 0.33 ;
  238.  
  239. ma_u11 = MATE mo_util 'E001' (EXTR epsp2 1) 'S001' (EXTR Y1 1)
  240. 'E002' (EXTR epsp2 2) 'S002' (EXTR Y1 2)
  241. 'E003' (EXTR epsp2 3) 'S003' (EXTR Y1 3)
  242. 'E004' (EXTR epsp2 4) 'S004' (EXTR Y1 4)
  243. 'E005' (EXTR epsp2 5) 'S005' (EXTR Y1 5) ;
  244.  
  245. ma_u12 = MATE mo_util 'E006' (EXTR epsp2 6) 'S006' (EXTR Y1 6)
  246. 'E007' (EXTR epsp2 7) 'S007' (EXTR Y1 7)
  247. 'E008' (EXTR epsp2 8) 'S008' (EXTR Y1 8)
  248. 'E009' (EXTR epsp2 9) 'S009' (EXTR Y1 9)
  249. 'E010' (EXTR epsp2 10) 'S010' (EXTR Y1 10) ;
  250.  
  251. ma_u21 = MATE mo_util 'E011' (EXTR epsp2 11) 'S011' (EXTR Y1 11)
  252. 'E012' (EXTR epsp2 12) 'S012' (EXTR Y1 12)
  253. 'E013' (EXTR epsp2 13) 'S013' (EXTR Y1 13)
  254. 'E014' (EXTR epsp2 14) 'S014' (EXTR Y1 14)
  255. 'E015' (EXTR epsp2 15) 'S015' (EXTR Y1 15) ;
  256.  
  257. ma_u22 = MATE mo_util 'E016' (EXTR epsp2 16) 'S016' (EXTR Y1 16)
  258. 'E017' (EXTR epsp2 17) 'S017' (EXTR Y1 17)
  259. 'E018' (EXTR epsp2 18) 'S018' (EXTR Y1 18)
  260. 'E019' (EXTR epsp2 19) 'S019' (EXTR Y1 19)
  261. 'E020' (EXTR epsp2 20) 'S020' (EXTR Y1 20) ;
  262.  
  263. ma_u31 = MATE mo_util 'E021' (EXTR epsp2 21) 'S021' (EXTR Y1 21)
  264. 'E022' (EXTR epsp2 22) 'S022' (EXTR Y1 22)
  265. 'E023' (EXTR epsp2 23) 'S023' (EXTR Y1 23)
  266. 'E024' (EXTR epsp2 24) 'S024' (EXTR Y1 24)
  267. 'E025' (EXTR epsp2 25) 'S025' (EXTR Y1 25) ;
  268.  
  269. ma_u32 = MATE mo_util 'E026' (EXTR epsp2 26) 'S026' (EXTR Y1 26)
  270. 'E027' (EXTR epsp2 27) 'S027' (EXTR Y1 27)
  271. 'E028' (EXTR epsp2 28) 'S028' (EXTR Y1 28)
  272. 'E029' (EXTR epsp2 29) 'S029' (EXTR Y1 29)
  273. 'E030' (EXTR epsp2 30) 'S030' (EXTR Y1 30) ;
  274.  
  275. ma_u41 = MATE mo_util 'E031' (EXTR epsp2 31) 'S031' (EXTR Y1 31)
  276. 'E032' (EXTR epsp2 32) 'S032' (EXTR Y1 32)
  277. 'E033' (EXTR epsp2 33) 'S033' (EXTR Y1 33)
  278. 'E034' (EXTR epsp2 34) 'S034' (EXTR Y1 34)
  279. 'E035' (EXTR epsp2 35) 'S035' (EXTR Y1 35) ;
  280.  
  281. ma_u42 = MATE mo_util 'E036' (EXTR epsp2 36) 'S036' (EXTR Y1 36)
  282. 'E037' (EXTR epsp2 37) 'S037' (EXTR Y1 37)
  283. 'E038' (EXTR epsp2 38) 'S038' (EXTR Y1 38)
  284. 'E039' (EXTR epsp2 39) 'S039' (EXTR Y1 39)
  285. 'E040' (EXTR epsp2 40) 'S040' (EXTR Y1 40) ;
  286.  
  287. ma_u51 = MATE mo_util 'E041' (EXTR epsp2 41) 'S041' (EXTR Y1 41)
  288. 'E042' (EXTR epsp2 42) 'S042' (EXTR Y1 42)
  289. 'E043' (EXTR epsp2 43) 'S043' (EXTR Y1 43)
  290. 'E044' (EXTR epsp2 44) 'S044' (EXTR Y1 44)
  291. 'E045' (EXTR epsp2 45) 'S045' (EXTR Y1 45) ;
  292.  
  293. ma_u52 = MATE mo_util 'E046' (EXTR epsp2 46) 'S046' (EXTR Y1 46)
  294. 'E047' (EXTR epsp2 47) 'S047' (EXTR Y1 47)
  295. 'E048' (EXTR epsp2 48) 'S048' (EXTR Y1 48)
  296. 'E049' (EXTR epsp2 49) 'S049' (EXTR Y1 49)
  297. 'E050' (EXTR epsp2 50) 'S050' (EXTR Y1 50) ;
  298.  
  299. ma_u61 = MATE mo_util 'E051' (EXTR epsp2 51) 'S051' (EXTR Y1 51)
  300. 'E052' (EXTR epsp2 52) 'S052' (EXTR Y1 52) ;
  301.  
  302. ma_util = ma_u0 et ma_u11 et ma_u12 et ma_u21 et ma_u22 et
  303. ma_u31 et ma_u32 et ma_u41 et ma_u42 et
  304. ma_u51 et ma_u52 et ma_u61 ;
  305.  
  306. * Resolution par PASAPAS------------------------------------------------
  307.  
  308. TABU = TABLE ;
  309. TABU.'MODELE' = mo_util ;
  310. TABU.'CARACTERISTIQUES' = ma_util ;
  311. TABU.'BLOCAGES_MECANIQUES' = cl0 ;
  312. TABU.'CHARGEMENT' = cha0 ;
  313. TABU.'TEMPS_CALCULES' = 'PROG' 0. 'PAS' dt0 tfin0 ;
  314. TMASAU=table;
  315. tabu . 'MES_SAUVEGARDES'=TMASAU;
  316. TMASAU .'DEFTO'=VRAI;
  317. TMASAU .'DEFIN'=VRAI;
  318. TEMPS 'ZERO' ;
  319. PASAPAS TABU ;
  320. TABTPS = TEMP 'NOEC';
  321. CPUext = TABTPS.'TEMPS_CPU'.'INITIAL';
  322.  
  323.  
  324. *=======================================================================
  325. * Comparaison des resultats obtenus par les deux modeles
  326. *
  327. * N.B. Comme on simule un essai de cisaillement YZ :
  328. * - Critere strict (errrel < 1.d-12) sur GAYZ, GIYZ et SMYZ
  329. * - Pour les autres composantes on note les depassements du critere
  330. * si impression demandee
  331.  
  332. ***opti echo 1 ;
  333.  
  334. uxc = TABLE ; uyc = TABLE ; uzc = TABLE ;
  335. epxxc = TABLE ; epyyc = TABLE ; epzzc = TABLE ;
  336. gaxyc = TABLE ; gaxzc = TABLE ; gayzc = TABLE ;
  337. eixxc = TABLE ; eiyyc = TABLE ; eizzc = TABLE ;
  338. gixyc = TABLE ; gixzc = TABLE ; giyzc = TABLE ;
  339. smxxc = TABLE ; smyyc = TABLE ; smzzc = TABLE ;
  340. smxyc = TABLE ; smxzc = TABLE ; smyzc = TABLE ;
  341. epsec = TABLE ;
  342.  
  343. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  344. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  345. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  346. enxxu = TABLE ; enyyu = TABLE ; enzzu = TABLE ;
  347. gnxyu = TABLE ; gnxzu = TABLE ; gnyzu = TABLE ;
  348. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  349. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  350. epseu = TABLE ;
  351.  
  352. LDEPL = MOTS 'UX ' 'UY ' 'UZ ' ;
  353. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  354. LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  355. LDENN = MOTS 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  356. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  357. LVARI = MOTS 'EPSE' ;
  358.  
  359. idefo = 6 ;
  360. idein = 6 ;
  361. icont = 6 ;
  362.  
  363. ierr = 0 ;
  364. critere = 1.d-12 ;
  365. IMESS = 0 ;
  366. ***IMESS = 1 ;
  367.  
  368. nbno = 8 ;
  369. nbgau = 8 ;
  370.  
  371. netat = DIME (ta1.'TEMPS') ;
  372. ietat = -1 ;
  373.  
  374. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  375.  
  376. REPETER BECART netat ;
  377.  
  378. ietat = ietat + 1 ;
  379.  
  380. temps1 = ta1.'TEMPS'.ietat ;
  381. temps2 = TABU.'TEMPS'.ietat ;
  382. SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  383. ierr = 1 ;
  384. QUITTER BECART ;
  385. FINSI ;
  386.  
  387. * Extraction des resultats de calcul : Modele interne CAST3M-----------
  388.  
  389. * Deplacements
  390. uxc.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  391. uxc.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  392. uxc.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  393. uxc.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  394. uxc.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  395. uxc.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  396. uxc.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  397. uxc.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  398.  
  399. uyc.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  400. uyc.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  401. uyc.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  402. uyc.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  403. uyc.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  404. uyc.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  405. uyc.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  406. uyc.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  407.  
  408. uzc.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  409. uzc.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  410. uzc.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  411. uzc.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  412. uzc.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  413. uzc.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  414. uzc.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  415. uzc.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  416.  
  417. * Deformations totales
  418. igau = 0 ;
  419. REPETER BEPXXC nbgau ;
  420. igau = igau + 1 ;
  421. epxxc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  422. FIN BEPXXC ;
  423.  
  424. igau = 0 ;
  425. REPETER BEPYYC nbgau ;
  426. igau = igau + 1 ;
  427. epyyc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  428. FIN BEPYYC ;
  429.  
  430. igau = 0 ;
  431. REPETER BEPZZC nbgau ;
  432. igau = igau + 1 ;
  433. epzzc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  434. FIN BEPZZC ;
  435.  
  436. igau = 0 ;
  437. REPETER BGAXYC nbgau ;
  438. igau = igau + 1 ;
  439. gaxyc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  440. FIN BGAXYC ;
  441.  
  442. igau = 0 ;
  443. REPETER BGAXZC nbgau ;
  444. igau = igau + 1 ;
  445. gaxzc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  446. FIN BGAXZC ;
  447.  
  448. igau = 0 ;
  449. REPETER BGAYZC nbgau ;
  450. igau = igau + 1 ;
  451. gayzc.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  452. FIN BGAYZC ;
  453.  
  454. * Deformations inelastiques
  455. igau = 0 ;
  456. REPETER BEIXXC nbgau ;
  457. igau = igau + 1 ;
  458. eixxc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  459. 'EIXX' 1 1 igau ;
  460. FIN BEIXXC ;
  461.  
  462. igau = 0 ;
  463. REPETER BEIYYC nbgau ;
  464. igau = igau + 1 ;
  465. eiyyc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  466. 'EIYY' 1 1 igau ;
  467. FIN BEIYYC ;
  468.  
  469. igau = 0 ;
  470. REPETER BEIZZC nbgau ;
  471. igau = igau + 1 ;
  472. eizzc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  473. 'EIZZ' 1 1 igau ;
  474. FIN BEIZZC ;
  475.  
  476. igau = 0 ;
  477. REPETER BGIXYC nbgau ;
  478. igau = igau + 1 ;
  479. gixyc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  480. 'GIXY' 1 1 igau ;
  481. FIN BGIXYC ;
  482.  
  483. igau = 0 ;
  484. REPETER BGIXZC nbgau ;
  485. igau = igau + 1 ;
  486. gixzc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  487. 'GIXZ' 1 1 igau ;
  488. FIN BGIXZC ;
  489.  
  490. igau = 0 ;
  491. REPETER BGIYZC nbgau ;
  492. igau = igau + 1 ;
  493. giyzc.igau = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ietat)
  494. 'GIYZ' 1 1 igau ;
  495. FIN BGIYZC ;
  496.  
  497. * Contraintes
  498. igau = 0 ;
  499. REPETER BSMXXC nbgau ;
  500. igau = igau + 1 ;
  501. smxxc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  502. FIN BSMXXC ;
  503.  
  504. igau = 0 ;
  505. REPETER BSMYYC nbgau ;
  506. igau = igau + 1 ;
  507. smyyc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  508. FIN BSMYYC ;
  509.  
  510. igau = 0 ;
  511. REPETER BSMZZC nbgau ;
  512. igau = igau + 1 ;
  513. smzzc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  514. FIN BSMZZC ;
  515.  
  516. igau = 0 ;
  517. REPETER BSMXYC nbgau ;
  518. igau = igau + 1 ;
  519. smxyc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  520. FIN BSMXYC ;
  521.  
  522. igau = 0 ;
  523. REPETER BSMXZC nbgau ;
  524. igau = igau + 1 ;
  525. smxzc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  526. FIN BSMXZC ;
  527.  
  528. igau = 0 ;
  529. REPETER BSMYZC nbgau ;
  530. igau = igau + 1 ;
  531. smyzc.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  532. FIN BSMYZC ;
  533.  
  534. * Variables internes
  535. igau = 0 ;
  536. REPETER BEPSEC nbgau ;
  537. igau = igau + 1 ;
  538. epsec.igau = EXTR (ta1.'VARIABLES_INTERNES'.ietat)
  539. 'EPSE' 1 1 igau ;
  540. FIN BEPSEC ;
  541.  
  542. * Extraction des resultats de calcul : Modele externe utilisateur------
  543.  
  544. uxu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  545. uxu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  546. uxu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  547. uxu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  548. uxu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  549. uxu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  550. uxu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  551. uxu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  552.  
  553. uyu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  554. uyu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  555. uyu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  556. uyu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  557. uyu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  558. uyu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  559. uyu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  560. uyu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  561.  
  562. uzu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  563. uzu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  564. uzu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  565. uzu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  566. uzu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  567. uzu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  568. uzu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  569. uzu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  570.  
  571. * Deformations totales
  572. igau = 0 ;
  573. REPETER BEPXXU nbgau ;
  574. igau = igau + 1 ;
  575. epxxu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  576. FIN BEPXXU ;
  577.  
  578. igau = 0 ;
  579. REPETER BEPYYU nbgau ;
  580. igau = igau + 1 ;
  581. epyyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  582. FIN BEPYYU ;
  583.  
  584. igau = 0 ;
  585. REPETER BEPZZU nbgau ;
  586. igau = igau + 1 ;
  587. epzzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  588. FIN BEPZZU ;
  589.  
  590. igau = 0 ;
  591. REPETER BGAXYU nbgau ;
  592. igau = igau + 1 ;
  593. gaxyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  594. FIN BGAXYU ;
  595.  
  596. igau = 0 ;
  597. REPETER BGAXZU nbgau ;
  598. igau = igau + 1 ;
  599. gaxzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  600. FIN BGAXZU ;
  601.  
  602. igau = 0 ;
  603. REPETER BGAYZU nbgau ;
  604. igau = igau + 1 ;
  605. gayzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  606. FIN BGAYZU ;
  607.  
  608. * Deformations inelastiques (variables internes supplementaires)
  609. igau = 0 ;
  610. REPETER BENXXU nbgau ;
  611. igau = igau + 1 ;
  612. enxxu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  613. 'ENXX' 1 1 igau ;
  614. FIN BENXXU ;
  615.  
  616. igau = 0 ;
  617. REPETER BENYYU nbgau ;
  618. igau = igau + 1 ;
  619. enyyu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  620. 'ENYY' 1 1 igau ;
  621. FIN BENYYU ;
  622.  
  623. igau = 0 ;
  624. REPETER BENZZU nbgau ;
  625. igau = igau + 1 ;
  626. enzzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  627. 'ENZZ' 1 1 igau ;
  628. FIN BENZZU ;
  629.  
  630. igau = 0 ;
  631. REPETER BGNXYU nbgau ;
  632. igau = igau + 1 ;
  633. gnxyu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  634. 'GNXY' 1 1 igau ;
  635. FIN BGNXYU ;
  636.  
  637. igau = 0 ;
  638. REPETER BGNXZU nbgau ;
  639. igau = igau + 1 ;
  640. gnxzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  641. 'GNXZ' 1 1 igau ;
  642. FIN BGNXZU ;
  643.  
  644. igau = 0 ;
  645. REPETER BGNYZU nbgau ;
  646. igau = igau + 1 ;
  647. gnyzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  648. 'GNYZ' 1 1 igau ;
  649. FIN BGNYZU ;
  650.  
  651. * Contraintes
  652. igau = 0 ;
  653. REPETER BSMXXU nbgau ;
  654. igau = igau + 1 ;
  655. smxxu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  656. FIN BSMXXU ;
  657.  
  658. igau = 0 ;
  659. REPETER BSMYYU nbgau ;
  660. igau = igau + 1 ;
  661. smyyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  662. FIN BSMYYU ;
  663.  
  664. igau = 0 ;
  665. REPETER BSMZZU nbgau ;
  666. igau = igau + 1 ;
  667. smzzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  668. FIN BSMZZU ;
  669.  
  670. igau = 0 ;
  671. REPETER BSMXYU nbgau ;
  672. igau = igau + 1 ;
  673. smxyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  674. FIN BSMXYU ;
  675.  
  676. igau = 0 ;
  677. REPETER BSMXZU nbgau ;
  678. igau = igau + 1 ;
  679. smxzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  680. FIN BSMXZU ;
  681.  
  682. igau = 0 ;
  683. REPETER BSMYZU nbgau ;
  684. igau = igau + 1 ;
  685. smyzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  686. FIN BSMYZU ;
  687.  
  688. * Variables internes
  689. igau = 0 ;
  690. REPETER BEPSEU nbgau ;
  691. igau = igau + 1 ;
  692. epseu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  693. 'EPSE' 1 1 igau ;
  694. FIN BEPSEU ;
  695.  
  696. * Comparaison des deplacements-----------------------------------------
  697.  
  698. icomp = 1 ;
  699. ino = 0 ;
  700. REPETER BUX nbno ;
  701. ino = ino + 1 ;
  702. SI (EGA uxc.ino 0.0) ;
  703. SI (NON (EGA uxu.ino 0.0)) ;
  704. SI (EGA IMESS 1) ;
  705. Message 'Erreur absolue' uxu.ino
  706. 'Composante ' (EXTR LDEPL icomp)
  707. ' Noeud' ino 'Etat' ietat ;
  708. Message ;
  709. FINSI ;
  710. ************ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  711. FINSI ;
  712. SINON ;
  713. dux = ABS (uxc.ino - uxu.ino) ;
  714. errrel = ABS (dux / uxc.ino) ;
  715. SI (errrel > critere) ;
  716. SI (EGA IMESS 1) ;
  717. Message 'Erreur relative' errrel
  718. 'Composante ' (EXTR LDEPL icomp)
  719. ' Noeud' ino 'Etat' ietat ;
  720. Message 'uxc.ino =' uxc.ino 'uxu.ino =' uxu.ino ;
  721. Message ;
  722. FINSI ;
  723. ************ierr = 11 ; QUITTER BECART ;
  724. FINSI ;
  725. FINSI ;
  726. FIN BUX ;
  727.  
  728. icomp = 2 ;
  729. ino = 0 ;
  730. REPETER BUY nbno ;
  731. ino = ino + 1 ;
  732. SI (EGA uyc.ino 0.0) ;
  733. SI (NON (EGA uyu.ino 0.0)) ;
  734. SI (EGA IMESS 1) ;
  735. Message 'Erreur absolue' uyu.ino
  736. 'Composante ' (EXTR LDEPL icomp)
  737. ' Noeud' ino 'Etat' ietat ;
  738. Message ;
  739. FINSI ;
  740. ************ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  741. FINSI ;
  742. SINON ;
  743. duy = ABS (uyc.ino - uyu.ino) ;
  744. errrel = ABS (duy / uyc.ino) ;
  745. SI (errrel > critere) ;
  746. SI (EGA IMESS 1) ;
  747. Message 'Erreur relative' errrel
  748. 'Composante ' (EXTR LDEPL icomp)
  749. ' Noeud' ino 'Etat' ietat ;
  750. Message 'uyc.ino =' uyc.ino 'uyu.ino =' uyu.ino ;
  751. Message ;
  752. FINSI ;
  753. ************ierr = 11 ; QUITTER BECART ;
  754. FINSI ;
  755. FINSI ;
  756. FIN BUY ;
  757.  
  758. icomp = 3 ;
  759. ino = 0 ;
  760. REPETER BUZ nbno ;
  761. ino = ino + 1 ;
  762. SI (EGA uzc.ino 0.0) ;
  763. SI (NON (EGA uzu.ino 0.0)) ;
  764. SI (EGA IMESS 1) ;
  765. Message 'Erreur absolue' uzu.ino
  766. 'Composante ' (EXTR LDEPL icomp)
  767. ' Noeud' ino 'Etat' ietat ;
  768. Message ;
  769. FINSI ;
  770. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  771. FINSI ;
  772. SINON ;
  773. duz = ABS (uzc.ino - uzu.ino) ;
  774. errrel = ABS (duz / uzc.ino) ;
  775. SI (errrel > critere) ;
  776. SI (EGA IMESS 1) ;
  777. Message 'Erreur relative' errrel
  778. 'Composante ' (EXTR LDEPL icomp)
  779. ' Noeud' ino 'Etat' ietat ;
  780. Message 'uzc.ino =' uzc.ino 'uzu.ino =' uzu.ino ;
  781. Message ;
  782. FINSI ;
  783. ************ierr = 11 ; QUITTER BECART ;
  784. FINSI ;
  785. FINSI ;
  786. FIN BUZ ;
  787.  
  788. * Comparaison des deformations-----------------------------------------
  789.  
  790. icomp = 1 ;
  791. igau = 0 ;
  792. REPETER BEPXX nbgau ;
  793. igau = igau + 1 ;
  794. SI (EGA epxxc.igau 0.0) ;
  795. SI (NON (EGA epxxu.igau 0.0)) ;
  796. SI (EGA IMESS 1) ;
  797. Message 'Erreur absolue' epxxu.igau
  798. 'Composante ' (EXTR LDEFO icomp)
  799. ' Point Gauss' igau 'Etat' ietat ;
  800. Message ;
  801. FINSI ;
  802. ************ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  803. FINSI ;
  804. SINON ;
  805. depxx = ABS (epxxc.igau - epxxu.igau) ;
  806. errrel = ABS (depxx / epxxc.igau) ;
  807. SI (errrel > critere) ;
  808. SI (EGA IMESS 1) ;
  809. Message 'Erreur relative' errrel
  810. 'Composante ' (EXTR LDEFO icomp)
  811. ' Point Gauss' igau 'Etat' ietat ;
  812. Message 'epxxc.igau =' epxxc.igau 'epxxu.igau =' epxxu.igau;
  813. Message ;
  814. FINSI ;
  815. ************ierr = 21 ; QUITTER BECART ;
  816. FINSI ;
  817. FINSI ;
  818. FIN BEPXX ;
  819.  
  820. icomp = 2 ;
  821. igau = 0 ;
  822. REPETER BEPYY nbgau ;
  823. igau = igau + 1 ;
  824. SI (EGA epyyc.igau 0.0) ;
  825. SI (NON (EGA epyyu.igau 0.0)) ;
  826. SI (EGA IMESS 1) ;
  827. Message 'Erreur absolue' epyyu.igau
  828. 'Composante ' (EXTR LDEFO icomp)
  829. ' Point Gauss' igau 'Etat' ietat ;
  830. Message ;
  831. FINSI ;
  832. ************ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  833. FINSI ;
  834. SINON ;
  835. depyy = ABS (epyyc.igau - epyyu.igau) ;
  836. errrel = ABS (depyy / epyyc.igau) ;
  837. SI (errrel > critere) ;
  838. SI (EGA IMESS 1) ;
  839. Message 'Erreur relative' errrel
  840. 'Composante ' (EXTR LDEFO icomp)
  841. ' Point Gauss' igau 'Etat' ietat ;
  842. Message 'epyyc.igau =' epyyc.igau 'epyyu.igau =' epyyu.igau;
  843. Message ;
  844. FINSI ;
  845. ************ierr = 21 ; QUITTER BECART ;
  846. FINSI ;
  847. FINSI ;
  848. FIN BEPYY ;
  849.  
  850. icomp = 3 ;
  851. igau = 0 ;
  852. REPETER BEPZZ nbgau ;
  853. igau = igau + 1 ;
  854. SI (EGA epzzc.igau 0.0) ;
  855. SI (NON (EGA epzzu.igau 0.0)) ;
  856. SI (EGA IMESS 1) ;
  857. Message 'Erreur absolue' epzzu.igau
  858. 'Composante ' (EXTR LDEFO icomp)
  859. ' Point Gauss' igau 'Etat' ietat ;
  860. Message ;
  861. FINSI ;
  862. ************ierr = 20 ; errabs = epzzu.igau ; QUITTER BECART ;
  863. FINSI ;
  864. SINON ;
  865. depzz = ABS (epzzc.igau - epzzu.igau) ;
  866. errrel = ABS (depzz / epzzc.igau) ;
  867. SI (errrel > critere) ;
  868. SI (EGA IMESS 1) ;
  869. Message 'Erreur relative' errrel
  870. 'Composante ' (EXTR LDEFO icomp)
  871. ' Point Gauss' igau 'Etat' ietat ;
  872. Message 'epzzc.igau =' epzzc.igau 'epzzu.igau =' epzzu.igau;
  873. Message ;
  874. FINSI ;
  875. ************ierr = 21 ; QUITTER BECART ;
  876. FINSI ;
  877. FINSI ;
  878. FIN BEPZZ ;
  879.  
  880. icomp = 4 ;
  881. igau = 0 ;
  882. REPETER BGAXY nbgau ;
  883. igau = igau + 1 ;
  884. SI (EGA gaxyc.igau 0.0) ;
  885. SI (NON (EGA gaxyu.igau 0.0)) ;
  886. SI (EGA IMESS 1) ;
  887. Message 'Erreur absolue' gaxyu.igau
  888. 'Composante ' (EXTR LDEFO icomp)
  889. ' Point Gauss' igau 'Etat' ietat ;
  890. Message ;
  891. FINSI ;
  892. ************ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  893. FINSI ;
  894. SINON ;
  895. dgaxy = ABS (gaxyc.igau - gaxyu.igau) ;
  896. errrel = ABS (dgaxy / gaxyc.igau) ;
  897. SI (errrel > critere) ;
  898. SI (EGA IMESS 1) ;
  899. Message 'Erreur relative' errrel
  900. 'Composante ' (EXTR LDEFO icomp)
  901. ' Point Gauss' igau 'Etat' ietat ;
  902. Message 'gaxyc.igau =' gaxyc.igau 'gaxyu.igau =' gaxyu.igau;
  903. Message ;
  904. FINSI ;
  905. ************ierr = 21 ; QUITTER BECART ;
  906. FINSI ;
  907. FINSI ;
  908. FIN BGAXY ;
  909.  
  910. icomp = 5 ;
  911. igau = 0 ;
  912. REPETER BGAXZ nbgau ;
  913. igau = igau + 1 ;
  914. SI (EGA gaxzc.igau 0.0) ;
  915. SI (NON (EGA gaxzu.igau 0.0)) ;
  916. SI (EGA IMESS 1) ;
  917. Message 'Erreur absolue' gaxzu.igau
  918. 'Composante ' (EXTR LDEFO icomp)
  919. ' Point Gauss' igau 'Etat' ietat ;
  920. Message ;
  921. FINSI ;
  922. ************ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  923. FINSI ;
  924. SINON ;
  925. dgaxz = ABS (gaxzc.igau - gaxzu.igau) ;
  926. errrel = ABS (dgaxz / gaxzc.igau) ;
  927. SI (errrel > critere) ;
  928. SI (EGA IMESS 1) ;
  929. Message 'Erreur relative' errrel
  930. 'Composante ' (EXTR LDEFO icomp)
  931. ' Point Gauss' igau 'Etat' ietat ;
  932. Message 'gaxzc.igau =' gaxzc.igau 'gaxzu.igau =' gaxzu.igau;
  933. Message ;
  934. FINSI ;
  935. ************ierr = 21 ; QUITTER BECART ;
  936. FINSI ;
  937. FINSI ;
  938. FIN BGAXZ ;
  939.  
  940. icomp = 6 ;
  941. igau = 0 ;
  942. REPETER BGAYZ nbgau ;
  943. igau = igau + 1 ;
  944. SI (EGA gayzc.igau 0.0) ;
  945. SI (NON (EGA gayzu.igau 0.0)) ;
  946. ************SI (EGA IMESS 1) ;
  947. ************Message 'Erreur absolue' gayzu.igau
  948. ************ 'Composante ' (EXTR LDEFO icomp)
  949. ************ ' Point Gauss' igau 'Etat' ietat ;
  950. ************Message ;
  951. ************FINSI ;
  952. ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  953. FINSI ;
  954. SINON ;
  955. dgayz = ABS (gayzc.igau - gayzu.igau) ;
  956. errrel = ABS (dgayz / gayzc.igau) ;
  957. SI (errrel > critere) ;
  958. ************SI (EGA IMESS 1) ;
  959. ************Message 'Erreur relative' errrel
  960. ************ 'Composante ' (EXTR LDEFO icomp)
  961. ************ ' Point Gauss' igau 'Etat' ietat ;
  962. ************Message 'gayzc.igau =' gayzc.igau 'gayzu.igau =' gayzu.igau;
  963. ************Message ;
  964. ************FINSI ;
  965. ierr = 21 ; QUITTER BECART ;
  966. FINSI ;
  967. FINSI ;
  968. FIN BGAYZ ;
  969.  
  970. * Comparaison des deformations inelastiques----------------------------
  971.  
  972. icomp = 1 ;
  973. igau = 0 ;
  974. REPETER BEIXX nbgau ;
  975. igau = igau + 1 ;
  976. SI (EGA eixxc.igau 0.0) ;
  977. SI (NON (EGA enxxu.igau 0.0)) ;
  978. SI (EGA IMESS 1) ;
  979. Message 'Erreur absolue' enxxu.igau
  980. 'Composante ' (EXTR LDEIN icomp)
  981. ' Point Gauss' igau 'Etat' ietat ;
  982. Message ;
  983. FINSI ;
  984. ************ierr = 30 ; errabs = enxxu.igau ; QUITTER BECART ;
  985. FINSI ;
  986. SINON ;
  987. deixx = ABS (eixxc.igau - enxxu.igau) ;
  988. errrel = ABS (deixx / eixxc.igau) ;
  989. SI (errrel > critere) ;
  990. SI (EGA IMESS 1) ;
  991. Message 'Erreur relative' errrel
  992. 'Composante ' (EXTR LDEIN icomp)
  993. ' Point Gauss' igau 'Etat' ietat ;
  994. Message 'eixxc.igau =' eixxc.igau 'enxxu.igau =' enxxu.igau;
  995. Message ;
  996. FINSI ;
  997. ************ierr = 31 ; QUITTER BECART ;
  998. FINSI ;
  999. FINSI ;
  1000. FIN BEIXX ;
  1001.  
  1002. icomp = 2 ;
  1003. igau = 0 ;
  1004. REPETER BEIYY nbgau ;
  1005. igau = igau + 1 ;
  1006. SI (EGA eiyyc.igau 0.0) ;
  1007. SI (NON (EGA enyyu.igau 0.0)) ;
  1008. SI (EGA IMESS 1) ;
  1009. Message 'Erreur absolue' enyyu.igau
  1010. 'Composante ' (EXTR LDEIN icomp)
  1011. ' Point Gauss' igau 'Etat' ietat ;
  1012. Message ;
  1013. FINSI ;
  1014. ************ierr = 30 ; errabs = enyyu.igau ; QUITTER BECART ;
  1015. FINSI ;
  1016. SINON ;
  1017. deiyy = ABS (eiyyc.igau - enyyu.igau) ;
  1018. errrel = ABS (deiyy / eiyyc.igau) ;
  1019. SI (errrel > critere) ;
  1020. SI (EGA IMESS 1) ;
  1021. Message 'Erreur relative' errrel
  1022. 'Composante ' (EXTR LDEIN icomp)
  1023. ' Point Gauss' igau 'Etat' ietat ;
  1024. Message 'eiyyc.igau =' eiyyc.igau 'enyyu.igau =' enyyu.igau;
  1025. Message ;
  1026. FINSI ;
  1027. ************ierr = 31 ; QUITTER BECART ;
  1028. FINSI ;
  1029. FINSI ;
  1030. FIN BEIYY ;
  1031.  
  1032. icomp = 3 ;
  1033. igau = 0 ;
  1034. REPETER BEIZZ nbgau ;
  1035. igau = igau + 1 ;
  1036. SI (EGA eizzc.igau 0.0) ;
  1037. SI (NON (EGA enzzu.igau 0.0)) ;
  1038. SI (EGA IMESS 1) ;
  1039. Message 'Erreur absolue' enzzu.igau
  1040. 'Composante ' (EXTR LDEIN icomp)
  1041. ' Point Gauss' igau 'Etat' ietat ;
  1042. Message ;
  1043. FINSI ;
  1044. ************ierr = 30 ; errabs = enzzu.igau ; QUITTER BECART ;
  1045. FINSI ;
  1046. SINON ;
  1047. deizz = ABS (eizzc.igau - enzzu.igau) ;
  1048. errrel = ABS (deizz / eizzc.igau) ;
  1049. SI (errrel > critere) ;
  1050. SI (EGA IMESS 1) ;
  1051. Message 'Erreur relative' errrel
  1052. 'Composante ' (EXTR LDEIN icomp)
  1053. ' Point Gauss' igau 'Etat' ietat ;
  1054. Message 'eizzc.igau =' eizzc.igau 'enzzu.igau =' enzzu.igau;
  1055. Message ;
  1056. FINSI ;
  1057. ************ierr = 31 ; QUITTER BECART ;
  1058. FINSI ;
  1059. FINSI ;
  1060. FIN BEIZZ ;
  1061.  
  1062. icomp = 4 ;
  1063. igau = 0 ;
  1064. REPETER BGIXY nbgau ;
  1065. igau = igau + 1 ;
  1066. SI (EGA gixyc.igau 0.0) ;
  1067. SI (NON (EGA gnxyu.igau 0.0)) ;
  1068. SI (EGA IMESS 1) ;
  1069. Message 'Erreur absolue' gnxyu.igau
  1070. 'Composante ' (EXTR LDEIN icomp)
  1071. ' Point Gauss' igau 'Etat' ietat ;
  1072. Message ;
  1073. FINSI ;
  1074. ************ierr = 30 ; errabs = gnxyu.igau ; QUITTER BECART ;
  1075. FINSI ;
  1076. SINON ;
  1077. dgixy = ABS (gixyc.igau - gnxyu.igau) ;
  1078. errrel = ABS (dgixy / gixyc.igau) ;
  1079. SI (errrel > critere) ;
  1080. SI (EGA IMESS 1) ;
  1081. Message 'Erreur relative' errrel
  1082. 'Composante ' (EXTR LDEIN icomp)
  1083. ' Point Gauss' igau 'Etat' ietat ;
  1084. Message 'gixyc.igau =' gixyc.igau 'gnxyu.igau =' gnxyu.igau;
  1085. Message ;
  1086. FINSI ;
  1087. ************ierr = 31 ; QUITTER BECART ;
  1088. FINSI ;
  1089. FINSI ;
  1090. FIN BGIXY ;
  1091.  
  1092. icomp = 5 ;
  1093. igau = 0 ;
  1094. REPETER BGIXZ nbgau ;
  1095. igau = igau + 1 ;
  1096. SI (EGA gixzc.igau 0.0) ;
  1097. SI (NON (EGA gnxzu.igau 0.0)) ;
  1098. SI (EGA IMESS 1) ;
  1099. Message 'Erreur absolue' gnxzu.igau
  1100. 'Composante ' (EXTR LDEIN icomp)
  1101. ' Point Gauss' igau 'Etat' ietat ;
  1102. Message ;
  1103. FINSI ;
  1104. ************ierr = 30 ; errabs = gnxzu.igau ; QUITTER BECART ;
  1105. FINSI ;
  1106. SINON ;
  1107. dgixz = ABS (gixzc.igau - gnxzu.igau) ;
  1108. errrel = ABS (dgixz / gixzc.igau) ;
  1109. SI (errrel > critere) ;
  1110. SI (EGA IMESS 1) ;
  1111. Message 'Erreur relative' errrel
  1112. 'Composante ' (EXTR LDEIN icomp)
  1113. ' Point Gauss' igau 'Etat' ietat ;
  1114. Message 'gixzc.igau =' gixzc.igau 'gnxzu.igau =' gnxzu.igau;
  1115. Message ;
  1116. FINSI ;
  1117. ************ierr = 31 ; QUITTER BECART ;
  1118. FINSI ;
  1119. FINSI ;
  1120. FIN BGIXZ ;
  1121.  
  1122. icomp = 6 ;
  1123. igau = 0 ;
  1124. REPETER BGIYZ nbgau ;
  1125. igau = igau + 1 ;
  1126. SI (EGA giyzc.igau 0.0) ;
  1127. SI (NON (EGA gnyzu.igau 0.0)) ;
  1128. ************SI (EGA IMESS 1) ;
  1129. ************Message 'Erreur absolue' gnyzu.igau
  1130. ************ 'Composante ' (EXTR LDEIN icomp)
  1131. ************ ' Point Gauss' igau 'Etat' ietat ;
  1132. ************Message ;
  1133. ************FINSI ;
  1134. ierr = 30 ; errabs = gnyzu.igau ; QUITTER BECART ;
  1135. FINSI ;
  1136. SINON ;
  1137. dgiyz = ABS (giyzc.igau - gnyzu.igau) ;
  1138. errrel = ABS (dgiyz / giyzc.igau) ;
  1139. SI (errrel > critere) ;
  1140. ************SI (EGA IMESS 1) ;
  1141. ************Message 'Erreur relative' errrel
  1142. ************ 'Composante ' (EXTR LDEIN icomp)
  1143. ************ ' Point Gauss' igau 'Etat' ietat ;
  1144. ************Message 'giyzc.igau =' giyzc.igau 'gnyzu.igau =' gnyzu.igau;
  1145. ************Message ;
  1146. ************FINSI ;
  1147. ierr = 31 ; QUITTER BECART ;
  1148. FINSI ;
  1149. FINSI ;
  1150. FIN BGIYZ ;
  1151.  
  1152. * Comparaison des contraintes------------------------------------------
  1153.  
  1154. icomp = 1 ;
  1155. igau = 0 ;
  1156. REPETER BSMXX nbgau ;
  1157. igau = igau + 1 ;
  1158. SI (EGA smxxc.igau 0.0) ;
  1159. SI (NON (EGA smxxu.igau 0.0)) ;
  1160. SI (EGA IMESS 1) ;
  1161. Message 'Erreur absolue' smxxu.igau
  1162. 'Composante ' (EXTR LCONT icomp)
  1163. ' Point Gauss' igau 'Etat' ietat ;
  1164. Message ;
  1165. FINSI ;
  1166. ************ierr = 40 ; errabs = smxxu.igau ; QUITTER BECART ;
  1167. FINSI ;
  1168. SINON ;
  1169. dsmxx = ABS (smxxc.igau - smxxu.igau) ;
  1170. errrel = ABS (dsmxx / smxxc.igau) ;
  1171. SI (errrel > critere) ;
  1172. SI (EGA IMESS 1) ;
  1173. Message 'Erreur relative' errrel
  1174. 'Composante ' (EXTR LCONT icomp)
  1175. ' Point Gauss' igau 'Etat' ietat ;
  1176. Message 'smxxc.igau =' smxxc.igau 'smxxu.igau =' smxxu.igau;
  1177. Message ;
  1178. FINSI ;
  1179. ************ierr = 41 ; QUITTER BECART ;
  1180. FINSI ;
  1181. FINSI ;
  1182. FIN BSMXX ;
  1183.  
  1184. icomp = 2 ;
  1185. igau = 0 ;
  1186. REPETER BSMYY nbgau ;
  1187. igau = igau + 1 ;
  1188. SI (EGA smyyc.igau 0.0) ;
  1189. SI (NON (EGA smyyu.igau 0.0)) ;
  1190. SI (EGA IMESS 1) ;
  1191. Message 'Erreur absolue' smyyu.igau
  1192. 'Composante ' (EXTR LCONT icomp)
  1193. ' Point Gauss' igau 'Etat' ietat ;
  1194. Message ;
  1195. FINSI ;
  1196. ************ierr = 40 ; errabs = smyyu.igau ; QUITTER BECART ;
  1197. FINSI ;
  1198. SINON ;
  1199. dsmyy = ABS (smyyc.igau - smyyu.igau) ;
  1200. errrel = ABS (dsmyy / smyyc.igau) ;
  1201. SI (errrel > critere) ;
  1202. SI (EGA IMESS 1) ;
  1203. Message 'Erreur relative' errrel
  1204. 'Composante ' (EXTR LCONT icomp)
  1205. ' Point Gauss' igau 'Etat' ietat ;
  1206. Message 'smyyc.igau =' smyyc.igau 'smyyu.igau =' smyyu.igau;
  1207. Message ;
  1208. FINSI ;
  1209. ************ierr = 41 ; QUITTER BECART ;
  1210. FINSI ;
  1211. FINSI ;
  1212. FIN BSMYY ;
  1213.  
  1214. icomp = 3 ;
  1215. igau = 0 ;
  1216. REPETER BSMZZ nbgau ;
  1217. igau = igau + 1 ;
  1218. SI (EGA smzzc.igau 0.0) ;
  1219. SI (NON (EGA smzzu.igau 0.0)) ;
  1220. SI (EGA IMESS 1) ;
  1221. Message 'Erreur absolue' smzzu.igau
  1222. 'Composante ' (EXTR LCONT icomp)
  1223. ' Point Gauss' igau 'Etat' ietat ;
  1224. Message ;
  1225. FINSI ;
  1226. ************ierr = 40 ; errabs = smzzu.igau ; QUITTER BECART ;
  1227. FINSI ;
  1228. SINON ;
  1229. dsmzz = ABS (smzzc.igau - smzzu.igau) ;
  1230. errrel = ABS (dsmzz / smzzc.igau) ;
  1231. SI (errrel > critere) ;
  1232. SI (EGA IMESS 1) ;
  1233. Message 'Erreur relative' errrel
  1234. 'Composante ' (EXTR LCONT icomp)
  1235. ' Point Gauss' igau 'Etat' ietat ;
  1236. Message 'smzzc.igau =' smzzc.igau 'smzzu.igau =' smzzu.igau;
  1237. Message ;
  1238. FINSI ;
  1239. ************ierr = 41 ; QUITTER BECART ;
  1240. FINSI ;
  1241. FINSI ;
  1242. FIN BSMZZ ;
  1243.  
  1244. icomp = 4 ;
  1245. igau = 0 ;
  1246. REPETER BSMXY nbgau ;
  1247. igau = igau + 1 ;
  1248. SI (EGA smxyc.igau 0.0) ;
  1249. SI (NON (EGA smxyu.igau 0.0)) ;
  1250. SI (EGA IMESS 1) ;
  1251. Message 'Erreur absolue' smxyu.igau
  1252. 'Composante ' (EXTR LCONT icomp)
  1253. ' Point Gauss' igau 'Etat' ietat ;
  1254. Message ;
  1255. FINSI ;
  1256. ************ierr = 40 ; errabs = smxyu.igau ; QUITTER BECART ;
  1257. FINSI ;
  1258. SINON ;
  1259. dsmxy = ABS (smxyc.igau - smxyu.igau) ;
  1260. errrel = ABS (dsmxy / smxyc.igau) ;
  1261. SI (errrel > critere) ;
  1262. SI (EGA IMESS 1) ;
  1263. Message 'Erreur relative' errrel
  1264. 'Composante ' (EXTR LCONT icomp)
  1265. ' Point Gauss' igau 'Etat' ietat ;
  1266. Message 'smxyc.igau =' smxyc.igau 'smxyu.igau =' smxyu.igau;
  1267. Message ;
  1268. FINSI ;
  1269. ************ierr = 41 ; QUITTER BECART ;
  1270. FINSI ;
  1271. FINSI ;
  1272. FIN BSMXY ;
  1273.  
  1274. icomp = 5 ;
  1275. igau = 0 ;
  1276. REPETER BSMXZ nbgau ;
  1277. igau = igau + 1 ;
  1278. SI (EGA smxzc.igau 0.0) ;
  1279. SI (NON (EGA smxzu.igau 0.0)) ;
  1280. SI (EGA IMESS 1) ;
  1281. Message 'Erreur absolue' smxzu.igau
  1282. 'Composante ' (EXTR LCONT icomp)
  1283. ' Point Gauss' igau 'Etat' ietat ;
  1284. Message ;
  1285. FINSI ;
  1286. ************ierr = 40 ; errabs = smxzu.igau ; QUITTER BECART ;
  1287. FINSI ;
  1288. SINON ;
  1289. dsmxz = ABS (smxzc.igau - smxzu.igau) ;
  1290. errrel = ABS (dsmxz / smxzc.igau) ;
  1291. SI (errrel > critere) ;
  1292. SI (EGA IMESS 1) ;
  1293. Message 'Erreur relative' errrel
  1294. 'Composante ' (EXTR LCONT icomp)
  1295. ' Point Gauss' igau 'Etat' ietat ;
  1296. Message 'smxzc.igau =' smxzc.igau 'smxzu.igau =' smxzu.igau;
  1297. Message ;
  1298. FINSI ;
  1299. ************ierr = 41 ; QUITTER BECART ;
  1300. FINSI ;
  1301. FINSI ;
  1302. FIN BSMXZ ;
  1303.  
  1304. icomp = 6 ;
  1305. igau = 0 ;
  1306. REPETER BSMYZ nbgau ;
  1307. igau = igau + 1 ;
  1308. SI (EGA smyzc.igau 0.0) ;
  1309. SI (NON (EGA smyzu.igau 0.0)) ;
  1310. ************SI (EGA IMESS 1) ;
  1311. ************Message 'Erreur absolue' smyzu.igau
  1312. ************ 'Composante ' (EXTR LCONT icomp)
  1313. ************ ' Point Gauss' igau 'Etat' ietat ;
  1314. ************Message ;
  1315. ************FINSI ;
  1316. ierr = 40 ; errabs = smyzu.igau ; QUITTER BECART ;
  1317. FINSI ;
  1318. SINON ;
  1319. dsmyz = ABS (smyzc.igau - smyzu.igau) ;
  1320. errrel = ABS (dsmyz / smyzc.igau) ;
  1321. SI (errrel > critere) ;
  1322. ************SI (EGA IMESS 1) ;
  1323. ************Message 'Erreur relative' errrel
  1324. ************ 'Composante ' (EXTR LCONT icomp)
  1325. ************ ' Point Gauss' igau 'Etat' ietat ;
  1326. ************Message 'smyzc.igau =' smyzc.igau 'smyzu.igau =' smyzu.igau;
  1327. ************Message ;
  1328. ************FINSI ;
  1329. ierr = 41 ; QUITTER BECART ;
  1330. FINSI ;
  1331. FINSI ;
  1332. FIN BSMYZ ;
  1333.  
  1334. * Comparaison des variables internes-----------------------------------
  1335.  
  1336. icomp = 1 ;
  1337. igau = 0 ;
  1338. REPETER BEPSE nbgau ;
  1339. igau = igau + 1 ;
  1340. SI (EGA epsec.igau 0.0) ;
  1341. SI (NON (EGA epseu.igau 0.0)) ;
  1342. SI (EGA IMESS 1) ;
  1343. Message 'Erreur absolue' epseu.igau
  1344. 'Composante ' (EXTR LVARI icomp)
  1345. ' Point Gauss' igau 'Etat' ietat ;
  1346. Message ;
  1347. FINSI ;
  1348. ************ierr = 50 ; errabs = epseu.igau ; QUITTER BECART ;
  1349. FINSI ;
  1350. SINON ;
  1351. depse = ABS (epsec.igau - epseu.igau) ;
  1352. errrel = ABS (depse / epsec.igau) ;
  1353. SI (errrel > critere) ;
  1354. SI (EGA IMESS 1) ;
  1355. Message 'Erreur relative' errrel
  1356. 'Composante ' (EXTR LVARI icomp)
  1357. ' Point Gauss' igau 'Etat' ietat ;
  1358. Message 'epsec.igau =' epsec.igau 'epseu.igau =' epseu.igau;
  1359. Message ;
  1360. FINSI ;
  1361. ************ierr = 51 ; QUITTER BECART ;
  1362. FINSI ;
  1363. FINSI ;
  1364. FIN BEPSE ;
  1365.  
  1366. FIN BECART ;
  1367.  
  1368. * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  1369.  
  1370. * MESSAGES D'ERREUR_____________________________________________________
  1371.  
  1372. Message ;
  1373. Message 'CPU avec modele interne = ' CPUint 'centiemes de seconde' ;
  1374. Message 'CPU avec modele externe = ' CPUext 'centiemes de seconde' ;
  1375.  
  1376. ifin = netat-1 ;
  1377. giyzcmax = 0.0 ;
  1378. igau = 0 ;
  1379. REPETER BGIYZFC nbgau ;
  1380. igau = igau + 1 ;
  1381. giyz = EXTR (ta1.'DEFORMATIONS_INELASTIQUES'.ifin)
  1382. 'GIYZ' 1 1 igau ;
  1383. SI (giyz > giyzcmax) ; giyzcmax = giyz ; FINSI ;
  1384. FIN BGIYZFC ;
  1385. giyzumax = 0.0 ;
  1386. igau = 0 ;
  1387. REPETER BGIYZFU nbgau ;
  1388. igau = igau + 1 ;
  1389. giyz = EXTR (TABU.'VARIABLES_INTERNES'.ifin)
  1390. 'GNYZ' 1 1 igau ;
  1391. SI (giyz > giyzumax) ; giyzumax = giyz ; FINSI ;
  1392. FIN BGIYZFU ;
  1393. Message ;
  1394. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  1395. Message 'GIYZ max avec modele interne = ' giyzcmax ;
  1396. Message 'GIYZ max avec modele externe = ' giyzumax ;
  1397. Message ;
  1398.  
  1399. SI (ierr > 0) ;
  1400.  
  1401. SI (EGA ierr 1) ;
  1402. Message 'Instants de calcul non coincidants' ;
  1403. ERRE 5 ;
  1404. FINSI ;
  1405.  
  1406. SI (EGA ierr 10) ;
  1407. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1408. ' Noeud' ino 'Etat' ietat ;
  1409. ERRE 5 ;
  1410. FINSI ;
  1411.  
  1412. SI (EGA ierr 11) ;
  1413. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1414. ' Noeud' ino 'Etat' ietat ;
  1415. ERRE 5 ;
  1416. FINSI ;
  1417.  
  1418. SI (EGA ierr 20) ;
  1419. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  1420. ' Point Gauss' igau 'Etat' ietat ;
  1421. ERRE 5 ;
  1422. FINSI ;
  1423.  
  1424. SI (EGA ierr 21) ;
  1425. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  1426. ' Point Gauss' igau 'Etat' ietat ;
  1427. ERRE 5 ;
  1428. FINSI ;
  1429.  
  1430. SI (EGA ierr 30) ;
  1431. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  1432. ' Point Gauss' igau 'Etat' ietat ;
  1433. ERRE 5 ;
  1434. FINSI ;
  1435.  
  1436. SI (EGA ierr 31) ;
  1437. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  1438. ' Point Gauss' igau 'Etat' ietat ;
  1439. ERRE 5 ;
  1440. FINSI ;
  1441.  
  1442. SI (EGA ierr 40) ;
  1443. Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  1444. ' Point Gauss' igau 'Etat' ietat ;
  1445. ERRE 5 ;
  1446. FINSI ;
  1447.  
  1448. SI (EGA ierr 41) ;
  1449. Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  1450. ' Point Gauss' igau 'Etat' ietat ;
  1451. ERRE 5 ;
  1452. FINSI ;
  1453.  
  1454. SI (EGA ierr 50) ;
  1455. Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  1456. ' Point Gauss' igau 'Etat' ietat ;
  1457. ERRE 5 ;
  1458. FINSI ;
  1459.  
  1460. SI (EGA ierr 51) ;
  1461. Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  1462. ' Point Gauss' igau 'Etat' ietat ;
  1463. ERRE 5 ;
  1464. FINSI ;
  1465.  
  1466. SINON ;
  1467.  
  1468. Message 'Resultats concordants entre les deux modeles' ;
  1469. Message 'Aucune erreur relative superieure a ' critere
  1470. ' pour les composantes' ;
  1471. Message (EXTR LDEFO idefo) ' ' (EXTR LDEIN idein)
  1472. ' ' (EXTR LCONT icont) ;
  1473. ERRE 0 ;
  1474.  
  1475. FINSI ;
  1476.  
  1477. FIN ;
  1478.  
  1479.  
  1480.  
  1481.  
  1482.  
  1483.  
  1484.  
  1485.  

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