Télécharger umat01.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : umat01.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 No1 : ELASTICITE ISOTROPE ISOTHERME
  15. *
  16. *=======================================================================
  17.  
  18.  
  19. opti echo 0 ;
  20.  
  21.  
  22. *=======================================================================
  23. *
  24. * PROCEDURES DEPENDANTES
  25. *
  26. *=======================================================================
  27.  
  28.  
  29.  
  30.  
  31. *=======================================================================
  32. *
  33. * JEU DE DONNEES PRINCIPAL
  34. *
  35. *=======================================================================
  36.  
  37.  
  38. *=======================================================================
  39. * Maillage
  40.  
  41. opti dime 3 elem cub8 ;
  42. option epsilon 'LINEAIRE';
  43.  
  44. long = 1.0 ;
  45. larg = 1.0 ;
  46. epai = 1.0 ;
  47.  
  48. p1 = 0.0 0.0 0.0 ;
  49. p2 = p1 PLUS (larg 0.0 0.0) ;
  50. p3 = p2 PLUS (0.0 long 0.0) ;
  51. p4 = p1 PLUS (0.0 long 0.0) ;
  52.  
  53. p11 = p1 PLUS (0.0 0.0 epai) ;
  54. p21 = p2 PLUS (0.0 0.0 epai) ;
  55. p31 = p3 PLUS (0.0 0.0 epai) ;
  56. p41 = p4 PLUS (0.0 0.0 epai) ;
  57.  
  58. p1p2 = DROIT 1 p1 p2 ;
  59. p4p3 = DROIT 1 p4 p3 ;
  60.  
  61. sbas = p1p2 REGL 1 p4p3 ;
  62.  
  63. p11p21 = DROIT 1 p11 p21 ;
  64. p41p31 = DROIT 1 p41 p31 ;
  65.  
  66. shaut = p11p21 REGL 1 p41p31 ;
  67.  
  68. cube = sbas VOLU 1 shaut ;
  69.  
  70.  
  71. *=======================================================================
  72. * Comportement elastique isotrope isotherme
  73. * Traction uniaxiale Z
  74.  
  75. *-----------------------------------------------------------------------
  76. * 1 - Calcul avec modele utilisateur
  77.  
  78. * 1.1 Definition du modele utilisateur
  79.  
  80. LCMAT = MOTS 'YOUN' 'NU ' ;
  81. mo_meca = MODE cube 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  82. 'NON_LINEAIRE' 'UTILISATEUR'
  83. 'NUME_LOI' 10 'C_MATERIAU' LCMAT ;
  84.  
  85. * 1.2 Affectation des proprietes de materiau
  86.  
  87. xyoun = 2.1e+11 ;
  88. xnu = 0.3 ;
  89. ma_meca = MATE mo_meca 'YOUN' xyoun 'NU' xnu ;
  90.  
  91. * 1.3 Evaluation de coefficients de la matrice de Hooke en vue du
  92. * calcul de la solution theorique
  93.  
  94. coef1 = xyoun / ((1.0+xnu) * (1.0-(2.0*xnu))) ;
  95. coezz = coef1 * (1.0-xnu) ;
  96. coexxyy = xnu * coef1 ;
  97.  
  98. * 1.4 Definition des chargements
  99.  
  100. * 1.4.0 Creation de progressions pour la ponderation des chargements,
  101. * definissant autant d'etats a calculer
  102.  
  103. L_temps = PROG 0.0 0.5 1.0 ;
  104. L_depl = PROG 0.0 0.5 1.0 ;
  105.  
  106. netat = DIME L_temps ;
  107.  
  108. * 1.4.1 Blocages mecaniques et deplacement axial impose
  109.  
  110. uz3 = 1.0e-2 ;
  111. clim1 = BLOQ sbas 'DEPL' ;
  112. clim2 = BLOQ shaut 'UX' 'UY' ;
  113. clim3 = BLOQ shaut 'UZ' ;
  114. depi3 = DEPI clim3 uz3 ;
  115. charmeca = CHAR 'DIMP' (EVOL MANU L_temps L_depl) depi3 ;
  116.  
  117. * 1.5 Resolution par PASAPAS
  118.  
  119. TABU = TABLE ;
  120. TABU.'MODELE' = mo_meca ;
  121. TABU.'CARACTERISTIQUES' = ma_meca ;
  122. TABU.'CHARGEMENT' = charmeca ;
  123. TABU.'BLOCAGES_MECANIQUES' = clim1 ET clim2 ET clim3 ;
  124. TABU.'TEMPS_CALCULES' = L_temps ;
  125. TMASAU=table;
  126. tabu . 'MES_SAUVEGARDES'=TMASAU;
  127. TMASAU .'DEFTO'=VRAI;
  128. TMASAU .'DEFIN'=VRAI;
  129. PASAPAS TABU ;
  130.  
  131. *-----------------------------------------------------------------------
  132. * 2 - Comparaison des resultats obtenus avec la solution theorique
  133.  
  134. * N.B. On simule un essai de traction ZZ
  135. * - Critere strict sur UZ , EPZZ, SMXX, SMYY et SMZZ :
  136. * errrel < 1.d-12
  137. * - Pour les autres composantes on note les depassements du critere
  138. * si impression demandee
  139.  
  140. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  141. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  142. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  143. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  144. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  145.  
  146. LDEPL = MOTS 'UX ' 'UY ' 'UZ ' ;
  147. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  148. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  149.  
  150. idepl = 3 ;
  151. idefo = 3 ;
  152. icon1 = 1 ;
  153. icon2 = 2 ;
  154. icon3 = 3 ;
  155.  
  156. ierr = 0 ;
  157. critere = 1.d-12 ;
  158. IMESS = 0 ;
  159. ***IMESS = 1 ;
  160.  
  161. nbno = 8 ;
  162. nbgau = 8 ;
  163.  
  164. ietat = 0 ;
  165.  
  166. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  167.  
  168. REPETER BECART (netat-1) ;
  169.  
  170. ietat = ietat + 1 ;
  171. temps = TABU.'TEMPS'.ietat ;
  172.  
  173. * Calcul de la solution theorique--------------------------------------
  174.  
  175. uz = temps * uz3 ;
  176. epzz = temps * (uz3 / epai) ;
  177. smxx = coexxyy * epzz ;
  178. smyy = coexxyy * epzz ;
  179. smzz = coezz * epzz ;
  180.  
  181. * Extraction des resultats de calcul : Modele externe utilisateur------
  182.  
  183. * Deplacements
  184. uxu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  185. uxu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  186. uxu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  187. uxu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  188. uxu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  189. uxu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  190. uxu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  191. uxu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  192.  
  193. uyu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  194. uyu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  195. uyu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  196. uyu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  197. uyu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  198. uyu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  199. uyu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  200. uyu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  201.  
  202. uzu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  203. uzu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  204. uzu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  205. uzu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  206. uzu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  207. uzu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  208. uzu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  209. uzu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  210.  
  211. * * Deformations totales
  212. igau = 0 ;
  213. REPETER BEPXXU nbgau ;
  214. igau = igau + 1 ;
  215. epxxu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  216. FIN BEPXXU ;
  217.  
  218. igau = 0 ;
  219. REPETER BEPYYU nbgau ;
  220. igau = igau + 1 ;
  221. epyyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  222. FIN BEPYYU ;
  223.  
  224. igau = 0 ;
  225. REPETER BEPZZU nbgau ;
  226. igau = igau + 1 ;
  227. epzzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  228. FIN BEPZZU ;
  229.  
  230. igau = 0 ;
  231. REPETER BGAXYU nbgau ;
  232. igau = igau + 1 ;
  233. gaxyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  234. FIN BGAXYU ;
  235.  
  236. igau = 0 ;
  237. REPETER BGAXZU nbgau ;
  238. igau = igau + 1 ;
  239. gaxzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  240. FIN BGAXZU ;
  241.  
  242. igau = 0 ;
  243. REPETER BGAYZU nbgau ;
  244. igau = igau + 1 ;
  245. gayzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  246. FIN BGAYZU ;
  247.  
  248. * Contraintes
  249. igau = 0 ;
  250. REPETER BSMXXU nbgau ;
  251. igau = igau + 1 ;
  252. smxxu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  253. FIN BSMXXU ;
  254.  
  255. igau = 0 ;
  256. REPETER BSMYYU nbgau ;
  257. igau = igau + 1 ;
  258. smyyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  259. FIN BSMYYU ;
  260.  
  261. igau = 0 ;
  262. REPETER BSMZZU nbgau ;
  263. igau = igau + 1 ;
  264. smzzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  265. FIN BSMZZU ;
  266.  
  267. igau = 0 ;
  268. REPETER BSMXYU nbgau ;
  269. igau = igau + 1 ;
  270. smxyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  271. FIN BSMXYU ;
  272.  
  273. igau = 0 ;
  274. REPETER BSMXZU nbgau ;
  275. igau = igau + 1 ;
  276. smxzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  277. FIN BSMXZU ;
  278.  
  279. igau = 0 ;
  280. REPETER BSMYZU nbgau ;
  281. igau = igau + 1 ;
  282. smyzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  283. FIN BSMYZU ;
  284.  
  285. * Comparaison des deplacements-----------------------------------------
  286.  
  287. icomp = 1 ;
  288. ino = 0 ;
  289. REPETER BUX nbno ;
  290. ino = ino + 1 ;
  291. SI (NON (EGA uxu.ino 0.0)) ;
  292. SI (EGA IMESS 1) ;
  293. Message 'Erreur absolue' uxu.ino
  294. 'Composante ' (EXTR LDEPL icomp)
  295. ' Noeud' ino 'Etat' ietat ;
  296. Message ;
  297. FINSI ;
  298. *********ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  299. FINSI ;
  300. FIN BUX ;
  301.  
  302. icomp = 2 ;
  303. ino = 0 ;
  304. REPETER BUY nbno ;
  305. ino = ino + 1 ;
  306. SI (NON (EGA uyu.ino 0.0)) ;
  307. SI (EGA IMESS 1) ;
  308. Message 'Erreur absolue' uyu.ino
  309. 'Composante ' (EXTR LDEPL icomp)
  310. ' Noeud' ino 'Etat' ietat ;
  311. Message ;
  312. FINSI ;
  313. *********ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  314. FINSI ;
  315. FIN BUY ;
  316.  
  317. icomp = 3 ;
  318. ino = 0 ;
  319. REPETER BUZ nbno ;
  320. ino = ino + 1 ;
  321. SI (ino < 5) ;
  322. SI (NON (EGA uzu.ino 0.0)) ;
  323. SI (EGA IMESS 1) ;
  324. Message 'Erreur absolue' uzu.ino
  325. 'Composante ' (EXTR LDEPL icomp)
  326. ' Noeud' ino 'Etat' ietat ;
  327. Message ;
  328. FINSI ;
  329. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  330. FINSI ;
  331. SINON ;
  332. duz = ABS (uz - uzu.ino) ;
  333. errrel = ABS (duz / uz) ;
  334. SI (errrel > critere) ;
  335. ************SI (EGA IMESS 1) ;
  336. ************Message 'Erreur relative' errrel
  337. ************ 'Composante ' (EXTR LDEPL icomp)
  338. ************ ' Noeud' ino 'Etat' ietat ;
  339. ************Message 'uz =' uz 'uzu.ino =' uzu.ino ;
  340. ************Message ;
  341. ************FINSI ;
  342. ierr = 11 ; QUITTER BECART ;
  343. FINSI ;
  344. FINSI ;
  345. FIN BUZ ;
  346.  
  347. * Comparaison des deformations totales---------------------------------
  348.  
  349. icomp = 1 ;
  350. igau = 0 ;
  351. REPETER BEPXX nbgau ;
  352. igau = igau + 1 ;
  353. SI (NON (EGA epxxu.igau 0.0)) ;
  354. SI (EGA IMESS 1) ;
  355. Message 'Erreur absolue' epxxu.igau
  356. 'Composante ' (EXTR LDEFO icomp)
  357. ' Point Gauss' igau 'Etat' ietat ;
  358. Message ;
  359. FINSI ;
  360. *********ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  361. FINSI ;
  362. FIN BEPXX ;
  363.  
  364. icomp = 2 ;
  365. igau = 0 ;
  366. REPETER BEPYY nbgau ;
  367. igau = igau + 1 ;
  368. SI (NON (EGA epyyu.igau 0.0)) ;
  369. SI (EGA IMESS 1) ;
  370. Message 'Erreur absolue' epyyu.igau
  371. 'Composante ' (EXTR LDEFO icomp)
  372. ' Point Gauss' igau 'Etat' ietat ;
  373. Message ;
  374. FINSI ;
  375. *********ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  376. FINSI ;
  377. FIN BEPYY ;
  378.  
  379. icomp = 3 ;
  380. igau = 0 ;
  381. REPETER BEPZZ nbgau ;
  382. igau = igau + 1 ;
  383. depzz = ABS (epzz - epzzu.igau) ;
  384. errrel = ABS (depzz / epzz) ;
  385. SI (errrel > critere) ;
  386. *********SI (EGA IMESS 1) ;
  387. *********Message 'Erreur relative' errrel
  388. ********* 'Composante ' (EXTR LDEFO icomp)
  389. ********* ' Point Gauss' igau 'Etat' ietat ;
  390. *********Message 'epzz =' epzz 'epzzu.igau =' epzzu.igau ;
  391. *********Message ;
  392. *********FINSI ;
  393. ierr = 21 ; QUITTER BECART ;
  394. FINSI ;
  395. FIN BEPZZ ;
  396.  
  397. icomp = 4 ;
  398. igau = 0 ;
  399. REPETER BGAXY nbgau ;
  400. igau = igau + 1 ;
  401. SI (NON (EGA gaxyu.igau 0.0)) ;
  402. SI (EGA IMESS 1) ;
  403. Message 'Erreur absolue' gaxyu.igau
  404. 'Composante ' (EXTR LDEFO icomp)
  405. ' Point Gauss' igau 'Etat' ietat ;
  406. Message ;
  407. FINSI ;
  408. *********ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  409. FINSI ;
  410. FIN BGAXY ;
  411.  
  412. icomp = 5 ;
  413. igau = 0 ;
  414. REPETER BGAXZ nbgau ;
  415. igau = igau + 1 ;
  416. SI (NON (EGA gaxzu.igau 0.0)) ;
  417. SI (EGA IMESS 1) ;
  418. Message 'Erreur absolue' gaxzu.igau
  419. 'Composante ' (EXTR LDEFO icomp)
  420. ' Point Gauss' igau 'Etat' ietat ;
  421. Message ;
  422. FINSI ;
  423. *********ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  424. FINSI ;
  425. FIN BGAXZ ;
  426.  
  427. icomp = 6 ;
  428. igau = 0 ;
  429. REPETER BGAYZ nbgau ;
  430. igau = igau + 1 ;
  431. SI (NON (EGA gayzu.igau 0.0)) ;
  432. SI (EGA IMESS 1) ;
  433. Message 'Erreur absolue' gayzu.igau
  434. 'Composante ' (EXTR LDEFO icomp)
  435. ' Point Gauss' igau 'Etat' ietat ;
  436. Message ;
  437. FINSI ;
  438. *********ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  439. FINSI ;
  440. FIN BGAYZ ;
  441.  
  442. * Comparaison des contraintes------------------------------------------
  443.  
  444. icomp = 1 ;
  445. igau = 0 ;
  446. REPETER BSMXX nbgau ;
  447. igau = igau + 1 ;
  448. dsmxx = ABS (smxx - smxxu.igau) ;
  449. errrel = ABS (dsmxx / smxx) ;
  450. SI (errrel > critere) ;
  451. *********SI (EGA IMESS 1) ;
  452. *********Message 'Erreur relative' errrel
  453. ********* 'Composante ' (EXTR LCONT icomp)
  454. ********* ' Point Gauss' igau 'Etat' ietat ;
  455. *********Message 'smxx =' smxx 'smxxu.igau =' smxxu.igau ;
  456. *********Message ;
  457. *********FINSI ;
  458. ierr = 31 ; QUITTER BECART ;
  459. FINSI ;
  460. FIN BSMXX ;
  461.  
  462. icomp = 2 ;
  463. igau = 0 ;
  464. REPETER BSMYY nbgau ;
  465. igau = igau + 1 ;
  466. dsmyy = ABS (smyy - smyyu.igau) ;
  467. errrel = ABS (dsmyy / smyy) ;
  468. SI (errrel > critere) ;
  469. *********SI (EGA IMESS 1) ;
  470. *********Message 'Erreur relative' errrel
  471. ********* 'Composante ' (EXTR LCONT icomp)
  472. ********* ' Point Gauss' igau 'Etat' ietat ;
  473. *********Message 'smyy =' smyy 'smyyu.igau =' smyyu.igau ;
  474. *********Message ;
  475. *********FINSI ;
  476. ierr = 31 ; QUITTER BECART ;
  477. FINSI ;
  478. FIN BSMYY ;
  479.  
  480. icomp = 3 ;
  481. igau = 0 ;
  482. REPETER BSMZZ nbgau ;
  483. igau = igau + 1 ;
  484. dsmzz = ABS (smzz - smzzu.igau) ;
  485. errrel = ABS (dsmzz / smzz) ;
  486. SI (errrel > critere) ;
  487. *********SI (EGA IMESS 1) ;
  488. *********Message 'Erreur relative' errrel
  489. ********* 'Composante ' (EXTR LCONT icomp)
  490. ********* ' Point Gauss' igau 'Etat' ietat ;
  491. *********Message 'smzz =' smzz 'smzzu.igau =' smzzu.igau ;
  492. *********Message ;
  493. *********FINSI ;
  494. ierr = 31 ; QUITTER BECART ;
  495. FINSI ;
  496. FIN BSMZZ ;
  497.  
  498. icomp = 4 ;
  499. igau = 0 ;
  500. REPETER BSMXY nbgau ;
  501. igau = igau + 1 ;
  502. SI (NON (EGA smxyu.igau 0.0)) ;
  503. SI (EGA IMESS 1) ;
  504. Message 'Erreur absolue' smxyu.igau
  505. 'Composante ' (EXTR LCONT icomp)
  506. ' Point Gauss' igau 'Etat' ietat ;
  507. Message ;
  508. FINSI ;
  509. *********ierr = 30 ; errabs = smxyu.igau ; QUITTER BECART ;
  510. FINSI ;
  511. FIN BSMXY ;
  512.  
  513. icomp = 5 ;
  514. igau = 0 ;
  515. REPETER BSMXZ nbgau ;
  516. igau = igau + 1 ;
  517. SI (NON (EGA smxzu.igau 0.0)) ;
  518. SI (EGA IMESS 1) ;
  519. Message 'Erreur absolue' smxzu.igau
  520. 'Composante ' (EXTR LCONT icomp)
  521. ' Point Gauss' igau 'Etat' ietat ;
  522. Message ;
  523. FINSI ;
  524. *********ierr = 30 ; errabs = smxzu.igau ; QUITTER BECART ;
  525. FINSI ;
  526. FIN BSMXZ ;
  527.  
  528. icomp = 6 ;
  529. igau = 0 ;
  530. REPETER BSMYZ nbgau ;
  531. igau = igau + 1 ;
  532. SI (NON (EGA smyzu.igau 0.0)) ;
  533. SI (EGA IMESS 1) ;
  534. Message 'Erreur absolue' smyzu.igau
  535. 'Composante ' (EXTR LCONT icomp)
  536. ' Point Gauss' igau 'Etat' ietat ;
  537. Message ;
  538. FINSI ;
  539. *********ierr = 30 ; errabs = smyzu.igau ; QUITTER BECART ;
  540. FINSI ;
  541. FIN BSMYZ ;
  542.  
  543. FIN BECART ;
  544.  
  545. * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  546.  
  547. * MESSAGES D'ERREUR_____________________________________________________
  548.  
  549. SI (ierr > 0) ;
  550.  
  551. SI (EGA ierr 10) ;
  552. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  553. ' Noeud' ino 'Etat' ietat ;
  554. ERRE 5 ;
  555. FINSI ;
  556.  
  557. SI (EGA ierr 11) ;
  558. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  559. ' Noeud' ino 'Etat' ietat ;
  560. ERRE 5 ;
  561. FINSI ;
  562.  
  563. SI (EGA ierr 20) ;
  564. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  565. ' Point Gauss' igau 'Etat' ietat ;
  566. ERRE 5 ;
  567. FINSI ;
  568.  
  569. SI (EGA ierr 21) ;
  570. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  571. ' Point Gauss' igau 'Etat' ietat ;
  572. ERRE 5 ;
  573. FINSI ;
  574.  
  575. SI (EGA ierr 30) ;
  576. Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  577. ' Point Gauss' igau 'Etat' ietat ;
  578. ERRE 5 ;
  579. FINSI ;
  580.  
  581. SI (EGA ierr 31) ;
  582. Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  583. ' Point Gauss' igau 'Etat' ietat ;
  584. ERRE 5 ;
  585. FINSI ;
  586.  
  587. SINON ;
  588.  
  589. Message 'Resultats concordants avec la solution theorique' ;
  590. Message 'Aucune erreur relative superieure a ' critere
  591. ' pour les composantes' ;
  592. Message (EXTR LDEPL idepl) ' ' (EXTR LDEFO idefo) ' '
  593. (EXTR LCONT icon1) ' ' (EXTR LCONT icon2) ' '
  594. (EXTR LCONT icon3) ;
  595. ERRE 0 ;
  596.  
  597. FINSI ;
  598.  
  599. FIN ;
  600.  
  601.  
  602.  
  603.  
  604.  
  605.  
  606.  
  607.  
  608.  
  609.  

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