Télécharger umat02.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : umat02.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 No2 : ELASTICITE ISOTROPE NON ISOTHERME
  15. *
  16. *=======================================================================
  17.  
  18.  
  19. opti echo 1 ;
  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.  
  43. long = 1.0 ;
  44. larg = 1.0 ;
  45. epai = 1.0 ;
  46.  
  47. p1 = 0.0 0.0 0.0 ;
  48. p2 = p1 PLUS (larg 0.0 0.0) ;
  49. p3 = p2 PLUS (0.0 long 0.0) ;
  50. p4 = p1 PLUS (0.0 long 0.0) ;
  51.  
  52. p11 = p1 PLUS (0.0 0.0 epai) ;
  53. p21 = p2 PLUS (0.0 0.0 epai) ;
  54. p31 = p3 PLUS (0.0 0.0 epai) ;
  55. p41 = p4 PLUS (0.0 0.0 epai) ;
  56.  
  57. p1p2 = DROIT 1 p1 p2 ;
  58. p4p3 = DROIT 1 p4 p3 ;
  59.  
  60. sbas = p1p2 REGL 1 p4p3 ;
  61.  
  62. p11p21 = DROIT 1 p11 p21 ;
  63. p41p31 = DROIT 1 p41 p31 ;
  64.  
  65. shaut = p11p21 REGL 1 p41p31 ;
  66.  
  67. cube = sbas VOLU 1 shaut ;
  68.  
  69.  
  70. *=======================================================================
  71. * Comportement elastique isotrope non isotherme
  72. * Traction uniaxiale Z
  73.  
  74. *-----------------------------------------------------------------------
  75. * 1 - Calcul avec le modele 'ELASTIQUE' 'ISOTROPE' de CAST3M
  76. * Resolution quasi-statique
  77.  
  78. * 1.1 Definition du modele elastique
  79.  
  80. mo_elas = MODE cube 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  81. epsilon lineaire ;
  82.  
  83. * 1.2 Affectation des proprietes de materiau
  84.  
  85. xyoun = 2.1e+11 ; xnu = 0.3 ; xalph = 1.0e-6 ; xtref = 0.0 ;
  86. xdydt = -2.1e+07 ; xdndt = -3.0e-05 ;
  87.  
  88. xyoun1 = xyoun - (xdydt * 10000.0) ;
  89. xnu1 = xnu - (xdndt * 10000.0) ;
  90. xyoun2 = xyoun + (xdydt * 10000.0) ;
  91. xnu2 = xnu + (xdndt * 10000.0) ;
  92.  
  93. L_T = PROG (xtref-10000.0) (xtref+10000.0) ;
  94. L_Y = PROG xyoun1 xyoun2 ;
  95. L_N = PROG xnu1 xnu2 ;
  96.  
  97. evoyoun = EVOL MANU 'T ' L_T 'YOUN' L_Y ;
  98. evonu = EVOL MANU 'T ' L_T 'NU' L_N ;
  99.  
  100. ma_elas = MATE mo_elas 'YOUN' evoyoun 'NU' evonu 'ALPH' xalph ;
  101.  
  102. * 1.3 Definition des chargements
  103.  
  104. * 1.3.0 Creation de progressions pour la ponderation des chargements,
  105. * definissant autant d'etats a calculer
  106.  
  107. L_temps = PROG 0.0 PAS 1.0 3.0 ;
  108. L_depl = PROG 0.0 PAS 1.0 3.0 ;
  109. L_teta = PROG 0.0 PAS 1.0 3.0 ;
  110.  
  111. netat = DIME L_temps ;
  112.  
  113. * 1.3.1 Blocages mecaniques et deplacement axial impose
  114.  
  115. uzimpo = 1.0e-2 ;
  116. clim1 = BLOQ sbas 'DEPL' ;
  117. clim2 = BLOQ shaut 'UX' 'UY' ;
  118. clim3 = BLOQ shaut 'UZ' ;
  119. depi3 = DEPI clim3 uzimpo ;
  120. chardimp = CHAR 'DIMP' (EVOL MANU L_temps L_depl) depi3 ;
  121.  
  122. * 1.3.2 Profil spatial nominal de temperature
  123.  
  124. tetaz0 = 0.0 ; tetaz1 = 300.0 ; dtetadz = (tetaz1-tetaz0)/epai ;
  125. ZZ = COOR 3 cube ;
  126. chteta = dtetadz * ZZ ;
  127. chteta = NOMC 'TETA' chteta ;
  128. chteta = CHAN ATTRIBUT chteta 'NATURE' 'DISCRET' ;
  129. charteta = CHAR 'TETA' chteta (EVOL MANU L_temps L_teta) ;
  130.  
  131. chT = NOMC 'T ' chteta ;
  132. chT = CHAN ATTRIBUT chT 'NATURE' 'DISCRET' ;
  133. charT = CHAR 'T ' chT (EVOL MANU L_temps L_teta) ;
  134.  
  135. * 1.4 Definition d'un champ de temperature de reference pour la
  136. * dilatation thermique
  137.  
  138. chpxtref = MANU CHPO cube 'SCAL' xtref ;
  139. chTref = NOMC 'T ' chpxtref ;
  140.  
  141. * 1.5 Resolution des etats successifs
  142.  
  143. * 1.5.1 Creation et initialisation de la table de resultats
  144.  
  145. ta1 = TABLE ;
  146.  
  147. ta1.'MATERIAU' = TABLE ;
  148. ta1.'MATERIAU' . 0 = VARI 'NUAG' mo_elas ma_elas chTref 'RIGIDITE' ;
  149.  
  150. ta1.'DEPLACEMENTS' = TABLE ;
  151. UX0 = MANU CHPO cube 'UX' 0.0 ;
  152. UY0 = MANU CHPO cube 'UY' 0.0 ;
  153. UZ0 = MANU CHPO cube 'UZ' 0.0 ;
  154. ta1.'DEPLACEMENTS' . 0 = COLI (COLI UX0 1.0 UY0 1.0) 1.0 UZ0 1.0 ;
  155.  
  156. ta1.'DEFORMATIONS' = TABLE ;
  157. ta1.'DEFORMATIONS' . 0 = ZERO mo_elas 'DEFORMAT' ;
  158.  
  159. ta1.'DEFORMATIONS_ELASTIQUES' = TABLE ;
  160. ta1.'DEFORMATIONS_ELASTIQUES' . 0 = ZERO mo_elas 'DEFORMAT' ;
  161.  
  162. ta1.'DEFORMATIONS_THERMIQUES' = TABLE ;
  163. ta1.'DEFORMATIONS_THERMIQUES' . 0 = ZERO mo_elas 'DEFORMAT' ;
  164.  
  165. ta1.'CONTRAINTES' = TABLE ;
  166. ta1.'CONTRAINTES' . 0 = ZERO mo_elas 'CONTRAIN' ;
  167.  
  168. * 1.5.2 Resolution quasi-statique
  169.  
  170. ietat = 0 ;
  171.  
  172. REPETER BETAT (netat-1) ;
  173.  
  174. ietat = ietat + 1 ;
  175. ietatm = ietat - 1 ;
  176.  
  177. coe_depl = EXTR L_depl (ietat+1) ;
  178. coe_teta = EXTR L_teta (ietat+1) ;
  179.  
  180. * Calcul des chargements imposes a l'instant courant :
  181. * - Champ de temperature et ecart par rapport a la temperature de
  182. * reference pour la dilatation thermique ;
  183. * - Deplacement impose.
  184.  
  185. chpoT = coe_teta * chT ;
  186. chpoT = NOMC 'SCAL' chpoT ;
  187. chpodT = COLI chpoT 1.0 chpxtref -1.0 ;
  188. chpoT = NOMC 'T ' chpoT ;
  189. chpodT = NOMC 'T ' chpodT ;
  190.  
  191. depietat = coe_depl * depi3 ;
  192.  
  193. * Champ de materiau a l'instant courant
  194.  
  195. ta1.'MATERIAU'.ietat = VARI 'NUAG' mo_elas ma_elas chpoT 'RIGIDITE' ;
  196.  
  197. * Rigidite structurelle a l'instant courant
  198. * Rigidite totale
  199.  
  200. rig_stru = RIGI mo_elas (ta1.'MATERIAU'.ietat) ;
  201. rig_etat = rig_stru ET clim1 ET clim2 ET clim3 ;
  202.  
  203. * Resolution du champ de deplacement a l'instant courant
  204.  
  205. ta1.'DEPLACEMENTS'.ietat = RESO rig_etat depietat ;
  206.  
  207. * Calcul des deformations totales a l'instant courant
  208.  
  209. ta1.'DEFORMATIONS'.ietat = EPSI mo_elas (ta1.'DEPLACEMENTS'.ietat) ;
  210.  
  211. * Calcul des deformations thermiques a l'instant courant
  212.  
  213. ta1.'DEFORMATIONS_THERMIQUES'.ietat = EPTH mo_elas
  214. (ta1.'MATERIAU'.ietat) chpodT ;
  215.  
  216. * Calcul des deformations elastiques a l'instant courant
  217.  
  218. ta1.'DEFORMATIONS_ELASTIQUES'.ietat = (ta1.'DEFORMATIONS'.ietat)
  219. - (ta1.'DEFORMATIONS_THERMIQUES'.ietat) ;
  220.  
  221. * Calcul des contraintes a l'instant courant
  222.  
  223. ta1.'CONTRAINTES'.ietat = ELAS mo_elas
  224. (ta1.'DEFORMATIONS_ELASTIQUES'.ietat)
  225. (ta1.'MATERIAU'.ietat) ;
  226.  
  227. FIN BETAT ;
  228.  
  229. *-----------------------------------------------------------------------
  230. * 2 - Calcul avec le modele 'ELASTIQUE' 'ISOTROPE' de CAST3M
  231. * Resolution par PASAPAS
  232.  
  233. ta2 = TABLE ;
  234. ta2.'MODELE' = mo_elas ;
  235. ta2.'CARACTERISTIQUES' = ma_elas ;
  236. ta2.'CHARGEMENT' = chardimp ET charT ;
  237. ta2.'TALPHA_REFERENCE' = xtref ;
  238. ta2.'BLOCAGES_MECANIQUES' = clim1 ET clim2 ET clim3 ;
  239. ta2.'TEMPS_CALCULES' = L_temps ;
  240. TMASAU=table;
  241. ta2 . 'MES_SAUVEGARDES'=TMASAU;
  242. TMASAU .'DEFTO'=VRAI;
  243. TMASAU .'DEFIN'=VRAI;
  244. PASAPAS ta2 ;
  245.  
  246. *-----------------------------------------------------------------------
  247. * 3 - Calcul avec modele utilisateur
  248. * Temperature 'TETA' definie en tant que parametre externe
  249. * => Instanciation du materiau et prise en charge de la dilatation
  250. * thermique dans UMAT
  251.  
  252. * 3.1 Definition du modele utilisateur
  253.  
  254. LCMAT11 = MOTS 'YOUN' 'NU ' 'ALPH' 'TREF' 'DYDT' 'DNDT' ;
  255. LCVAR11 = MOTS 'EEXX' 'EEYY' 'EEZZ' 'EEXY' 'EEXZ' 'EEYZ'
  256. 'ETXX' 'ETYY' 'ETZZ' 'ETXY' 'ETXZ' 'ETYZ' ;
  257. LCPAR11 = MOTS 'TETA' ;
  258. opti epsilon lineaire;
  259. mo_util = MODE cube 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  260. 'NON_LINEAIRE' 'UTILISATEUR'
  261. 'NUME_LOI' 11 'PARA_LOI' LCPAR11
  262. 'C_MATERIAU' LCMAT11 'C_VARINTER' LCVAR11;
  263.  
  264. * 3.2 Affectation des proprietes de materiau
  265.  
  266. ma_util = MATE mo_util 'YOUN' xyoun 'NU' xnu 'ALPH' xalph 'TREF' xtref
  267. 'DYDT' xdydt 'DNDT' xdndt ;
  268.  
  269. * 3.3 Resolution par PASAPAS
  270.  
  271. TABU = TABLE ;
  272. TABU.'MODELE' = mo_util ;
  273. TABU.'CARACTERISTIQUES' = ma_util ;
  274. TABU.'CHARGEMENT' = chardimp ET charteta ;
  275. TABU.'BLOCAGES_MECANIQUES' = clim1 ET clim2 ET clim3 ;
  276. TABU.'TEMPS_CALCULES' = L_temps ;
  277. TMASAU=table;
  278. tabu . 'MES_SAUVEGARDES'=TMASAU;
  279. TMASAU .'DEFTO'=VRAI;
  280. TMASAU .'DEFIN'=VRAI;
  281. tabU.'MOVA'='EEXX';
  282. PASAPAS TABU ;
  283.  
  284. opti echo 0;
  285. *=======================================================================
  286. * Comparaison des resultats obtenus par les trois modeles
  287. *
  288. * N.B. Comme on simule un essai de traction ZZ
  289. * - Critere strict (errrel < 1.d-12) sur EPZZ, EEZZ, ETZZ, SMXX,
  290. * SMYY et SMZZ
  291. * - Pour les autres composantes on note les depassements du critere
  292. * si impression demandee
  293.  
  294. uxc1 = TABLE ; uyc1 = TABLE ; uzc1 = TABLE ;
  295. epxxc1 = TABLE ; epyyc1 = TABLE ; epzzc1 = TABLE ;
  296. gaxyc1 = TABLE ; gaxzc1 = TABLE ; gayzc1 = TABLE ;
  297. eexxc1 = TABLE ; eeyyc1 = TABLE ; eezzc1 = TABLE ;
  298. gexyc1 = TABLE ; gexzc1 = TABLE ; geyzc1 = TABLE ;
  299. etxxc1 = TABLE ; etyyc1 = TABLE ; etzzc1 = TABLE ;
  300. gtxyc1 = TABLE ; gtxzc1 = TABLE ; gtyzc1 = TABLE ;
  301. smxxc1 = TABLE ; smyyc1 = TABLE ; smzzc1 = TABLE ;
  302. smxyc1 = TABLE ; smxzc1 = TABLE ; smyzc1 = TABLE ;
  303.  
  304. uxc2 = TABLE ; uyc2 = TABLE ; uzc2 = TABLE ;
  305. smxxc2 = TABLE ; smyyc2 = TABLE ; smzzc2 = TABLE ;
  306. smxyc2 = TABLE ; smxzc2 = TABLE ; smyzc2 = TABLE ;
  307.  
  308. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  309. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  310. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  311. eexxu = TABLE ; eeyyu = TABLE ; eezzu = TABLE ;
  312. gexyu = TABLE ; gexzu = TABLE ; geyzu = TABLE ;
  313. etxxu = TABLE ; etyyu = TABLE ; etzzu = TABLE ;
  314. gtxyu = TABLE ; gtxzu = TABLE ; gtyzu = TABLE ;
  315. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  316. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  317.  
  318. LDEPL = MOTS 'UX ' 'UY ' 'UZ ' ;
  319. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  320. LDEFE = MOTS 'EEXX' 'EEYY' 'EEZZ' 'GEXY' 'GEXZ' 'GEYZ' ;
  321. LDEFT = MOTS 'ETXX' 'ETYY' 'ETZZ' 'GTXY' 'GTXZ' 'GTYZ' ;
  322. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  323.  
  324. idefo = 3 ;
  325. idefe = 3 ;
  326. ideft = 3 ;
  327. icon1 = 1 ;
  328. icon2 = 2 ;
  329. icon3 = 3 ;
  330.  
  331. ierr = 0 ;
  332. critere = 1.d-12 ;
  333. IMESS = 0 ;
  334. ***IMESS = 1 ;
  335.  
  336. nbno = 8 ;
  337. nbgau = 8 ;
  338.  
  339. ietat = -1 ;
  340.  
  341. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  342.  
  343. REPETER BECART netat ;
  344.  
  345. ietat = ietat + 1 ;
  346.  
  347. temps1 = EXTR L_temps (ietat+1) ;
  348. temps2 = ta2.'TEMPS'.ietat ;
  349. temps3 = TABU.'TEMPS'.ietat ;
  350. SI (OU ((ABS (temps1-temps2)) > 1.0E-9)
  351. ((ABS (temps2-temps3)) > 1.0E-9) ) ;
  352. ierr = 1 ;
  353. QUITTER BECART ;
  354. FINSI ;
  355.  
  356. * Extraction des resultats de calcul : Modele interne CAST3M-----------
  357. * Resultats issus de la resolution quasi-statique
  358.  
  359. * Deplacements
  360. uxc1.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  361. uxc1.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  362. uxc1.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  363. uxc1.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  364. uxc1.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  365. uxc1.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  366. uxc1.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  367. uxc1.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  368.  
  369. uyc1.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  370. uyc1.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  371. uyc1.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  372. uyc1.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  373. uyc1.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  374. uyc1.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  375. uyc1.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  376. uyc1.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  377.  
  378. uzc1.1 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  379. uzc1.2 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  380. uzc1.3 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  381. uzc1.4 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  382. uzc1.5 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  383. uzc1.6 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  384. uzc1.7 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  385. uzc1.8 = EXTR (ta1.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  386.  
  387. * Deformations totales
  388. igau = 0 ;
  389. REPETER BEPXXC1 nbgau ;
  390. igau = igau + 1 ;
  391. epxxc1.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  392. FIN BEPXXC1 ;
  393.  
  394. igau = 0 ;
  395. REPETER BEPYYC1 nbgau ;
  396. igau = igau + 1 ;
  397. epyyc1.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  398. FIN BEPYYC1 ;
  399.  
  400. igau = 0 ;
  401. REPETER BEPZZC1 nbgau ;
  402. igau = igau + 1 ;
  403. epzzc1.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  404. FIN BEPZZC1 ;
  405.  
  406. igau = 0 ;
  407. REPETER BGAXYC1 nbgau ;
  408. igau = igau + 1 ;
  409. gaxyc1.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  410. FIN BGAXYC1 ;
  411.  
  412. igau = 0 ;
  413. REPETER BGAXZC1 nbgau ;
  414. igau = igau + 1 ;
  415. gaxzc1.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  416. FIN BGAXZC1 ;
  417.  
  418. igau = 0 ;
  419. REPETER BGAYZC1 nbgau ;
  420. igau = igau + 1 ;
  421. gayzc1.igau = EXTR (ta1.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  422. FIN BGAYZC1 ;
  423.  
  424. * Deformations elastiques
  425. igau = 0 ;
  426. REPETER BEEXXC1 nbgau ;
  427. igau = igau + 1 ;
  428. eexxc1.igau = EXTR (ta1.'DEFORMATIONS_ELASTIQUES'.ietat)
  429. 'EPXX' 1 1 igau ;
  430. FIN BEEXXC1 ;
  431.  
  432. igau = 0 ;
  433. REPETER BEEYYC1 nbgau ;
  434. igau = igau + 1 ;
  435. eeyyc1.igau = EXTR (ta1.'DEFORMATIONS_ELASTIQUES'.ietat)
  436. 'EPYY' 1 1 igau ;
  437. FIN BEEYYC1 ;
  438.  
  439. igau = 0 ;
  440. REPETER BEEZZC1 nbgau ;
  441. igau = igau + 1 ;
  442. eezzc1.igau = EXTR (ta1.'DEFORMATIONS_ELASTIQUES'.ietat)
  443. 'EPZZ' 1 1 igau ;
  444. FIN BEEZZC1 ;
  445.  
  446. igau = 0 ;
  447. REPETER BGEXYC1 nbgau ;
  448. igau = igau + 1 ;
  449. gexyc1.igau = EXTR (ta1.'DEFORMATIONS_ELASTIQUES'.ietat)
  450. 'GAXY' 1 1 igau ;
  451. FIN BGEXYC1 ;
  452.  
  453. igau = 0 ;
  454. REPETER BGEXZC1 nbgau ;
  455. igau = igau + 1 ;
  456. gexzc1.igau = EXTR (ta1.'DEFORMATIONS_ELASTIQUES'.ietat)
  457. 'GAXZ' 1 1 igau ;
  458. FIN BGEXZC1 ;
  459.  
  460. igau = 0 ;
  461. REPETER BGEYZC1 nbgau ;
  462. igau = igau + 1 ;
  463. geyzc1.igau = EXTR (ta1.'DEFORMATIONS_ELASTIQUES'.ietat)
  464. 'GAYZ' 1 1 igau ;
  465. FIN BGEYZC1 ;
  466.  
  467. * Deformations thermiques
  468. igau = 0 ;
  469. REPETER BETXXC1 nbgau ;
  470. igau = igau + 1 ;
  471. etxxc1.igau = EXTR (ta1.'DEFORMATIONS_THERMIQUES'.ietat)
  472. 'EPXX' 1 1 igau ;
  473. FIN BETXXC1 ;
  474.  
  475. igau = 0 ;
  476. REPETER BETYYC1 nbgau ;
  477. igau = igau + 1 ;
  478. etyyc1.igau = EXTR (ta1.'DEFORMATIONS_THERMIQUES'.ietat)
  479. 'EPYY' 1 1 igau ;
  480. FIN BETYYC1 ;
  481.  
  482. igau = 0 ;
  483. REPETER BETZZC1 nbgau ;
  484. igau = igau + 1 ;
  485. etzzc1.igau = EXTR (ta1.'DEFORMATIONS_THERMIQUES'.ietat)
  486. 'EPZZ' 1 1 igau ;
  487. FIN BETZZC1 ;
  488.  
  489. igau = 0 ;
  490. REPETER BGTXYC1 nbgau ;
  491. igau = igau + 1 ;
  492. gtxyc1.igau = EXTR (ta1.'DEFORMATIONS_THERMIQUES'.ietat)
  493. 'GAXY' 1 1 igau ;
  494. FIN BGTXYC1 ;
  495.  
  496. igau = 0 ;
  497. REPETER BGTXZC1 nbgau ;
  498. igau = igau + 1 ;
  499. gtxzc1.igau = EXTR (ta1.'DEFORMATIONS_THERMIQUES'.ietat)
  500. 'GAXZ' 1 1 igau ;
  501. FIN BGTXZC1 ;
  502.  
  503. igau = 0 ;
  504. REPETER BGTYZC1 nbgau ;
  505. igau = igau + 1 ;
  506. gtyzc1.igau = EXTR (ta1.'DEFORMATIONS_THERMIQUES'.ietat)
  507. 'GAYZ' 1 1 igau ;
  508. FIN BGTYZC1 ;
  509.  
  510. * Contraintes
  511. igau = 0 ;
  512. REPETER BSMXXC1 nbgau ;
  513. igau = igau + 1 ;
  514. smxxc1.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  515. FIN BSMXXC1 ;
  516.  
  517. igau = 0 ;
  518. REPETER BSMYYC1 nbgau ;
  519. igau = igau + 1 ;
  520. smyyc1.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  521. FIN BSMYYC1 ;
  522.  
  523. igau = 0 ;
  524. REPETER BSMZZC1 nbgau ;
  525. igau = igau + 1 ;
  526. smzzc1.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  527. FIN BSMZZC1 ;
  528.  
  529. igau = 0 ;
  530. REPETER BSMXYC1 nbgau ;
  531. igau = igau + 1 ;
  532. smxyc1.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  533. FIN BSMXYC1 ;
  534.  
  535. igau = 0 ;
  536. REPETER BSMXZC1 nbgau ;
  537. igau = igau + 1 ;
  538. smxzc1.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  539. FIN BSMXZC1 ;
  540.  
  541. igau = 0 ;
  542. REPETER BSMYZC1 nbgau ;
  543. igau = igau + 1 ;
  544. smyzc1.igau = EXTR (ta1.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  545. FIN BSMYZC1 ;
  546.  
  547. * Extraction des resultats de calcul : Modele interne CAST3M-----------
  548. * Resultats issus de la resolution par PASAPAS
  549.  
  550. * Deplacements
  551. uxc2.1 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  552. uxc2.2 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  553. uxc2.3 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  554. uxc2.4 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  555. uxc2.5 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  556. uxc2.6 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  557. uxc2.7 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  558. uxc2.8 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  559.  
  560. uyc2.1 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  561. uyc2.2 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  562. uyc2.3 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  563. uyc2.4 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  564. uyc2.5 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  565. uyc2.6 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  566. uyc2.7 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  567. uyc2.8 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  568.  
  569. uzc2.1 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  570. uzc2.2 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  571. uzc2.3 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  572. uzc2.4 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  573. uzc2.5 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  574. uzc2.6 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  575. uzc2.7 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  576. uzc2.8 = EXTR (ta2.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  577.  
  578. * Contraintes
  579. igau = 0 ;
  580. REPETER BSMXXC2 nbgau ;
  581. igau = igau + 1 ;
  582. smxxc2.igau = EXTR (ta2.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  583. FIN BSMXXC2 ;
  584.  
  585. igau = 0 ;
  586. REPETER BSMYYC2 nbgau ;
  587. igau = igau + 1 ;
  588. smyyc2.igau = EXTR (ta2.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  589. FIN BSMYYC2 ;
  590.  
  591. igau = 0 ;
  592. REPETER BSMZZC2 nbgau ;
  593. igau = igau + 1 ;
  594. smzzc2.igau = EXTR (ta2.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  595. FIN BSMZZC2 ;
  596.  
  597. igau = 0 ;
  598. REPETER BSMXYC2 nbgau ;
  599. igau = igau + 1 ;
  600. smxyc2.igau = EXTR (ta2.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  601. FIN BSMXYC2 ;
  602.  
  603. igau = 0 ;
  604. REPETER BSMXZC2 nbgau ;
  605. igau = igau + 1 ;
  606. smxzc2.igau = EXTR (ta2.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  607. FIN BSMXZC2 ;
  608.  
  609. igau = 0 ;
  610. REPETER BSMYZC2 nbgau ;
  611. igau = igau + 1 ;
  612. smyzc2.igau = EXTR (ta2.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  613. FIN BSMYZC2 ;
  614.  
  615. * Extraction des resultats de calcul : Modele externe utilisateur------
  616.  
  617. * Deplacements
  618. uxu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  619. uxu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  620. uxu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  621. uxu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  622. uxu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  623. uxu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  624. uxu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  625. uxu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  626.  
  627. uyu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  628. uyu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  629. uyu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  630. uyu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  631. uyu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  632. uyu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  633. uyu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  634. uyu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  635.  
  636. uzu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  637. uzu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  638. uzu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  639. uzu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  640. uzu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  641. uzu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  642. uzu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  643. uzu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  644.  
  645. * Deformations totales
  646. igau = 0 ;
  647. REPETER BEPXXU nbgau ;
  648. igau = igau + 1 ;
  649. epxxu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  650. FIN BEPXXU ;
  651.  
  652. igau = 0 ;
  653. REPETER BEPYYU nbgau ;
  654. igau = igau + 1 ;
  655. epyyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  656. FIN BEPYYU ;
  657.  
  658. igau = 0 ;
  659. REPETER BEPZZU nbgau ;
  660. igau = igau + 1 ;
  661. epzzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  662. FIN BEPZZU ;
  663.  
  664. igau = 0 ;
  665. REPETER BGAXYU nbgau ;
  666. igau = igau + 1 ;
  667. gaxyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  668. FIN BGAXYU ;
  669.  
  670. igau = 0 ;
  671. REPETER BGAXZU nbgau ;
  672. igau = igau + 1 ;
  673. gaxzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  674. FIN BGAXZU ;
  675.  
  676. igau = 0 ;
  677. REPETER BGAYZU nbgau ;
  678. igau = igau + 1 ;
  679. gayzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  680. FIN BGAYZU ;
  681.  
  682. * Deformations elastiques
  683. igau = 0 ;
  684. REPETER BEEXXU nbgau ;
  685. igau = igau + 1 ;
  686. eexxu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  687. 'EEXX' 1 1 igau ;
  688. FIN BEEXXU ;
  689.  
  690. igau = 0 ;
  691. REPETER BEEYYU nbgau ;
  692. igau = igau + 1 ;
  693. eeyyu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  694. 'EEYY' 1 1 igau ;
  695. FIN BEEYYU ;
  696.  
  697. igau = 0 ;
  698. REPETER BEEZZU nbgau ;
  699. igau = igau + 1 ;
  700. eezzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  701. 'EEZZ' 1 1 igau ;
  702. FIN BEEZZU ;
  703.  
  704. igau = 0 ;
  705. REPETER BGEXYU nbgau ;
  706. igau = igau + 1 ;
  707. gexyu.igau = (EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  708. 'EEXY' 1 1 igau) * 2.d0 ;
  709. FIN BGEXYU ;
  710.  
  711. igau = 0 ;
  712. REPETER BGEXZU nbgau ;
  713. igau = igau + 1 ;
  714. gexzu.igau = (EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  715. 'EEXZ' 1 1 igau) * 2.d0 ;
  716. FIN BGEXZU ;
  717.  
  718. igau = 0 ;
  719. REPETER BGEYZU nbgau ;
  720. igau = igau + 1 ;
  721. geyzu.igau = (EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  722. 'EEYZ' 1 1 igau) * 2.d0 ;
  723. FIN BGEYZU ;
  724.  
  725. * Deformations thermiques
  726. igau = 0 ;
  727. REPETER BETXXU nbgau ;
  728. igau = igau + 1 ;
  729. etxxu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  730. 'ETXX' 1 1 igau ;
  731. FIN BETXXU ;
  732.  
  733. igau = 0 ;
  734. REPETER BETYYU nbgau ;
  735. igau = igau + 1 ;
  736. etyyu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  737. 'ETYY' 1 1 igau ;
  738. FIN BETYYU ;
  739.  
  740. igau = 0 ;
  741. REPETER BETZZU nbgau ;
  742. igau = igau + 1 ;
  743. etzzu.igau = EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  744. 'ETZZ' 1 1 igau ;
  745. FIN BETZZU ;
  746.  
  747. igau = 0 ;
  748. REPETER BGTXYU nbgau ;
  749. igau = igau + 1 ;
  750. gtxyu.igau = (EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  751. 'ETXY' 1 1 igau) * 2.d0 ;
  752. FIN BGTXYU ;
  753.  
  754. igau = 0 ;
  755. REPETER BGTXZU nbgau ;
  756. igau = igau + 1 ;
  757. gtxzu.igau = (EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  758. 'ETXZ' 1 1 igau) * 2.d0 ;
  759. FIN BGTXZU ;
  760.  
  761. igau = 0 ;
  762. REPETER BGTYZU nbgau ;
  763. igau = igau + 1 ;
  764. gtyzu.igau = (EXTR (TABU.'VARIABLES_INTERNES'.ietat)
  765. 'ETYZ' 1 1 igau) * 2.d0 ;
  766. FIN BGTYZU ;
  767.  
  768. * Contraintes
  769. igau = 0 ;
  770. REPETER BSMXXU nbgau ;
  771. igau = igau + 1 ;
  772. smxxu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  773. FIN BSMXXU ;
  774.  
  775. igau = 0 ;
  776. REPETER BSMYYU nbgau ;
  777. igau = igau + 1 ;
  778. smyyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  779. FIN BSMYYU ;
  780.  
  781. igau = 0 ;
  782. REPETER BSMZZU nbgau ;
  783. igau = igau + 1 ;
  784. smzzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  785. FIN BSMZZU ;
  786.  
  787. igau = 0 ;
  788. REPETER BSMXYU nbgau ;
  789. igau = igau + 1 ;
  790. smxyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  791. FIN BSMXYU ;
  792.  
  793. igau = 0 ;
  794. REPETER BSMXZU nbgau ;
  795. igau = igau + 1 ;
  796. smxzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  797. FIN BSMXZU ;
  798.  
  799. igau = 0 ;
  800. REPETER BSMYZU nbgau ;
  801. igau = igau + 1 ;
  802. smyzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  803. FIN BSMYZU ;
  804.  
  805. * Comparaison des deplacements-----------------------------------------
  806.  
  807. icomp = 1 ;
  808. ino = 0 ;
  809. REPETER BUX1 nbno ;
  810. ino = ino + 1 ;
  811. SI (EGA uxc1.ino 0.0) ;
  812. SI (NON (EGA uxu.ino 0.0)) ;
  813. SI (EGA IMESS 1) ;
  814. Message 'Erreur absolue' uxu.ino
  815. 'par rapport a la resolution quasi-statique' ;
  816. Message 'Composante ' (EXTR LDEPL icomp)
  817. ' Noeud' ino 'Etat' ietat ;
  818. Message ;
  819. FINSI ;
  820. ************ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  821. FINSI ;
  822. SINON ;
  823. dux = ABS (uxc1.ino - uxu.ino) ;
  824. errrel = ABS (dux / uxc1.ino) ;
  825. SI (errrel > critere) ;
  826. SI (EGA IMESS 1) ;
  827. Message 'Erreur relative' errrel
  828. 'par rapport a la resolution quasi-statique' ;
  829. Message 'Composante ' (EXTR LDEPL icomp)
  830. ' Noeud' ino 'Etat' ietat ;
  831. Message 'uxc1.ino =' uxc1.ino 'uxu.ino =' uxu.ino ;
  832. Message ;
  833. FINSI ;
  834. ************ierr = 11 ; QUITTER BECART ;
  835. FINSI ;
  836. FINSI ;
  837. FIN BUX1 ;
  838. ino = 0 ;
  839. REPETER BUX2 nbno ;
  840. ino = ino + 1 ;
  841. SI (EGA uxc2.ino 0.0) ;
  842. SI (NON (EGA uxu.ino 0.0)) ;
  843. SI (EGA IMESS 1) ;
  844. Message 'Erreur absolue' uxu.ino
  845. 'par rapport a la resolution par PASAPAS' ;
  846. Message 'Composante ' (EXTR LDEPL icomp)
  847. ' Noeud' ino 'Etat' ietat ;
  848. Message ;
  849. FINSI ;
  850. ************ierr = 12 ; errabs = uxu.ino ; QUITTER BECART ;
  851. FINSI ;
  852. SINON ;
  853. dux = ABS (uxc2.ino - uxu.ino) ;
  854. errrel = ABS (dux / uxc2.ino) ;
  855. SI (errrel > critere) ;
  856. SI (EGA IMESS 1) ;
  857. Message 'Erreur relative' errrel
  858. 'par rapport a la resolution par PASAPAS' ;
  859. Message 'Composante ' (EXTR LDEPL icomp)
  860. ' Noeud' ino 'Etat' ietat ;
  861. Message 'uxc2.ino =' uxc2.ino 'uxu.ino =' uxu.ino ;
  862. Message ;
  863. FINSI ;
  864. ************ierr = 13 ; QUITTER BECART ;
  865. FINSI ;
  866. FINSI ;
  867. FIN BUX2 ;
  868.  
  869. icomp = 2 ;
  870. ino = 0 ;
  871. REPETER BUY1 nbno ;
  872. ino = ino + 1 ;
  873. SI (EGA uyc1.ino 0.0) ;
  874. SI (NON (EGA uyu.ino 0.0)) ;
  875. SI (EGA IMESS 1) ;
  876. Message 'Erreur absolue' uyu.ino
  877. 'par rapport a la resolution quasi-statique' ;
  878. Message 'Composante ' (EXTR LDEPL icomp)
  879. ' Noeud' ino 'Etat' ietat ;
  880. Message ;
  881. FINSI ;
  882. ************ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  883. FINSI ;
  884. SINON ;
  885. duy = ABS (uyc1.ino - uyu.ino) ;
  886. errrel = ABS (duy / uyc1.ino) ;
  887. SI (errrel > critere) ;
  888. SI (EGA IMESS 1) ;
  889. Message 'Erreur relative' errrel
  890. 'par rapport a la resolution quasi-statique' ;
  891. Message 'Composante ' (EXTR LDEPL icomp)
  892. ' Noeud' ino 'Etat' ietat ;
  893. Message 'uyc1.ino =' uyc1.ino 'uyu.ino =' uyu.ino ;
  894. Message ;
  895. FINSI ;
  896. ************ierr = 11 ; QUITTER BECART ;
  897. FINSI ;
  898. FINSI ;
  899. FIN BUY1 ;
  900. ino = 0 ;
  901. REPETER BUY2 nbno ;
  902. ino = ino + 1 ;
  903. SI (EGA uyc2.ino 0.0) ;
  904. SI (NON (EGA uyu.ino 0.0)) ;
  905. SI (EGA IMESS 1) ;
  906. Message 'Erreur absolue' uyu.ino
  907. 'par rapport a la resolution par PASAPAS' ;
  908. Message 'Composante ' (EXTR LDEPL icomp)
  909. ' Noeud' ino 'Etat' ietat ;
  910. Message ;
  911. FINSI ;
  912. ************ierr = 12 ; errabs = uyu.ino ; QUITTER BECART ;
  913. FINSI ;
  914. SINON ;
  915. duy = ABS (uyc2.ino - uyu.ino) ;
  916. errrel = ABS (duy / uyc2.ino) ;
  917. SI (errrel > critere) ;
  918. SI (EGA IMESS 1) ;
  919. Message 'Erreur relative' errrel
  920. 'par rapport a la resolution par PASAPAS' ;
  921. Message 'Composante ' (EXTR LDEPL icomp)
  922. ' Noeud' ino 'Etat' ietat ;
  923. Message 'uyc2.ino =' uyc2.ino 'uyu.ino =' uyu.ino ;
  924. Message ;
  925. FINSI ;
  926. ************ierr = 13 ; QUITTER BECART ;
  927. FINSI ;
  928. FINSI ;
  929. FIN BUY2 ;
  930.  
  931. icomp = 3 ;
  932. ino = 0 ;
  933. REPETER BUZ1 nbno ;
  934. ino = ino + 1 ;
  935. SI (EGA uzc1.ino 0.0) ;
  936. SI (NON (EGA uzu.ino 0.0)) ;
  937. SI (EGA IMESS 1) ;
  938. Message 'Erreur absolue' uzu.ino
  939. 'par rapport a la resolution quasi-statique' ;
  940. Message 'Composante ' (EXTR LDEPL icomp)
  941. ' Noeud' ino 'Etat' ietat ;
  942. Message ;
  943. FINSI ;
  944. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  945. FINSI ;
  946. SINON ;
  947. duz = ABS (uzc1.ino - uzu.ino) ;
  948. errrel = ABS (duz / uzc1.ino) ;
  949. SI (errrel > critere) ;
  950. SI (EGA IMESS 1) ;
  951. Message 'Erreur relative' errrel
  952. 'par rapport a la resolution quasi-statique' ;
  953. Message 'Composante ' (EXTR LDEPL icomp)
  954. ' Noeud' ino 'Etat' ietat ;
  955. Message 'uzc1.ino =' uzc1.ino 'uzu.ino =' uzu.ino ;
  956. Message ;
  957. FINSI ;
  958. ************ierr = 11 ; QUITTER BECART ;
  959. FINSI ;
  960. FINSI ;
  961. FIN BUZ1 ;
  962. ino = 0 ;
  963. REPETER BUZ2 nbno ;
  964. ino = ino + 1 ;
  965. SI (EGA uzc2.ino 0.0) ;
  966. SI (NON (EGA uzu.ino 0.0)) ;
  967. SI (EGA IMESS 1) ;
  968. Message 'Erreur absolue' uzu.ino
  969. 'par rapport a la resolution par PASAPAS' ;
  970. Message 'Composante ' (EXTR LDEPL icomp)
  971. ' Noeud' ino 'Etat' ietat ;
  972. Message ;
  973. FINSI ;
  974. ************ierr = 12 ; errabs = uzu.ino ; QUITTER BECART ;
  975. FINSI ;
  976. SINON ;
  977. duz = ABS (uzc2.ino - uzu.ino) ;
  978. errrel = ABS (duz / uzc2.ino) ;
  979. SI (errrel > critere) ;
  980. SI (EGA IMESS 1) ;
  981. Message 'Erreur relative' errrel
  982. 'par rapport a la resolution par PASAPAS' ;
  983. Message 'Composante ' (EXTR LDEPL icomp)
  984. ' Noeud' ino 'Etat' ietat ;
  985. Message 'uzc2.ino =' uzc2.ino 'uzu.ino =' uzu.ino ;
  986. Message ;
  987. FINSI ;
  988. ************ierr = 13 ; QUITTER BECART ;
  989. FINSI ;
  990. FINSI ;
  991. FIN BUZ2 ;
  992.  
  993. * Comparaison des deformations totales---------------------------------
  994.  
  995. icomp = 1 ;
  996. igau = 0 ;
  997. REPETER BEPXX nbgau ;
  998. igau = igau + 1 ;
  999. SI (EGA epxxc1.igau 0.0) ;
  1000. SI (NON (EGA epxxu.igau 0.0)) ;
  1001. SI (EGA IMESS 1) ;
  1002. Message 'Erreur absolue' epxxu.igau
  1003. 'Composante ' (EXTR LDEFO icomp)
  1004. ' Point Gauss' igau 'Etat' ietat ;
  1005. Message ;
  1006. FINSI ;
  1007. ************ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  1008. FINSI ;
  1009. SINON ;
  1010. depxx = ABS (epxxc1.igau - epxxu.igau) ;
  1011. errrel = ABS (depxx / epxxc1.igau) ;
  1012. SI (errrel > critere) ;
  1013. SI (EGA IMESS 1) ;
  1014. Message 'Erreur relative' errrel
  1015. 'Composante ' (EXTR LDEFO icomp)
  1016. ' Point Gauss' igau 'Etat' ietat ;
  1017. Message 'epxxc1.igau =' epxxc1.igau 'epxxu.igau ='
  1018. epxxu.igau;
  1019. Message ;
  1020. FINSI ;
  1021. ************ierr = 21 ; QUITTER BECART ;
  1022. FINSI ;
  1023. FINSI ;
  1024. FIN BEPXX ;
  1025.  
  1026. icomp = 2 ;
  1027. igau = 0 ;
  1028. REPETER BEPYY nbgau ;
  1029. igau = igau + 1 ;
  1030. SI (EGA epyyc1.igau 0.0) ;
  1031. SI (NON (EGA epyyu.igau 0.0)) ;
  1032. SI (EGA IMESS 1) ;
  1033. Message 'Erreur absolue' epyyu.igau
  1034. 'Composante ' (EXTR LDEFO icomp)
  1035. ' Point Gauss' igau 'Etat' ietat ;
  1036. Message ;
  1037. FINSI ;
  1038. ************ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  1039. FINSI ;
  1040. SINON ;
  1041. depyy = ABS (epyyc1.igau - epyyu.igau) ;
  1042. errrel = ABS (depyy / epyyc1.igau) ;
  1043. SI (errrel > critere) ;
  1044. SI (EGA IMESS 1) ;
  1045. Message 'Erreur relative' errrel
  1046. 'Composante ' (EXTR LDEFO icomp)
  1047. ' Point Gauss' igau 'Etat' ietat ;
  1048. Message 'epyyc1.igau =' epyyc1.igau 'epyyu.igau ='
  1049. epyyu.igau;
  1050. Message ;
  1051. FINSI ;
  1052. ************ierr = 21 ; QUITTER BECART ;
  1053. FINSI ;
  1054. FINSI ;
  1055. FIN BEPYY ;
  1056.  
  1057. icomp = 3 ;
  1058. igau = 0 ;
  1059. REPETER BEPZZ nbgau ;
  1060. igau = igau + 1 ;
  1061. SI (EGA epzzc1.igau 0.0) ;
  1062. SI (NON (EGA epzzu.igau 0.0)) ;
  1063. ************SI (EGA IMESS 1) ;
  1064. ************Message 'Erreur absolue' epzzu.igau
  1065. ************ 'Composante ' (EXTR LDEFO icomp)
  1066. ************ ' Point Gauss' igau 'Etat' ietat ;
  1067. ************Message ;
  1068. ************FINSI ;
  1069. ierr = 20 ; errabs = epzzu.igau ; QUITTER BECART ;
  1070. FINSI ;
  1071. SINON ;
  1072. depzz = ABS (epzzc1.igau - epzzu.igau) ;
  1073. errrel = ABS (depzz / epzzc1.igau) ;
  1074. SI (errrel > critere) ;
  1075. ************SI (EGA IMESS 1) ;
  1076. ************Message 'Erreur relative' errrel
  1077. ************ 'Composante ' (EXTR LDEFO icomp)
  1078. ************ ' Point Gauss' igau 'Etat' ietat ;
  1079. ************Message 'epzzc1.igau =' epzzc1.igau 'epzzu.igau ='
  1080. ************ epzzu.igau;
  1081. ************Message ;
  1082. ************FINSI ;
  1083. ierr = 21 ; QUITTER BECART ;
  1084. FINSI ;
  1085. FINSI ;
  1086. FIN BEPZZ ;
  1087.  
  1088. icomp = 4 ;
  1089. igau = 0 ;
  1090. REPETER BGAXY nbgau ;
  1091. igau = igau + 1 ;
  1092. SI (EGA gaxyc1.igau 0.0) ;
  1093. SI (NON (EGA gaxyu.igau 0.0)) ;
  1094. SI (EGA IMESS 1) ;
  1095. Message 'Erreur absolue' gaxyu.igau
  1096. 'Composante ' (EXTR LDEFO icomp)
  1097. ' Point Gauss' igau 'Etat' ietat ;
  1098. Message ;
  1099. FINSI ;
  1100. ************ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  1101. FINSI ;
  1102. SINON ;
  1103. dgaxy = ABS (gaxyc1.igau - gaxyu.igau) ;
  1104. errrel = ABS (dgaxy / gaxyc1.igau) ;
  1105. SI (errrel > critere) ;
  1106. SI (EGA IMESS 1) ;
  1107. Message 'Erreur relative' errrel
  1108. 'Composante ' (EXTR LDEFO icomp)
  1109. ' Point Gauss' igau 'Etat' ietat ;
  1110. Message 'gaxyc1.igau =' gaxyc1.igau 'gaxyu.igau ='
  1111. gaxyu.igau;
  1112. Message ;
  1113. FINSI ;
  1114. ************ierr = 21 ; QUITTER BECART ;
  1115. FINSI ;
  1116. FINSI ;
  1117. FIN BGAXY ;
  1118.  
  1119. icomp = 5 ;
  1120. igau = 0 ;
  1121. REPETER BGAXZ nbgau ;
  1122. igau = igau + 1 ;
  1123. SI (EGA gaxzc1.igau 0.0) ;
  1124. SI (NON (EGA gaxzu.igau 0.0)) ;
  1125. SI (EGA IMESS 1) ;
  1126. Message 'Erreur absolue' gaxzu.igau
  1127. 'Composante ' (EXTR LDEFO icomp)
  1128. ' Point Gauss' igau 'Etat' ietat ;
  1129. Message ;
  1130. FINSI ;
  1131. ************ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  1132. FINSI ;
  1133. SINON ;
  1134. dgaxz = ABS (gaxzc1.igau - gaxzu.igau) ;
  1135. errrel = ABS (dgaxz / gaxzc1.igau) ;
  1136. SI (errrel > critere) ;
  1137. SI (EGA IMESS 1) ;
  1138. Message 'Erreur relative' errrel
  1139. 'Composante ' (EXTR LDEFO icomp)
  1140. ' Point Gauss' igau 'Etat' ietat ;
  1141. Message 'gaxzc1.igau =' gaxzc1.igau 'gaxzu.igau ='
  1142. gaxzu.igau;
  1143. Message ;
  1144. FINSI ;
  1145. ************ierr = 21 ; QUITTER BECART ;
  1146. FINSI ;
  1147. FINSI ;
  1148. FIN BGAXZ ;
  1149.  
  1150. icomp = 6 ;
  1151. igau = 0 ;
  1152. REPETER BGAYZ nbgau ;
  1153. igau = igau + 1 ;
  1154. SI (EGA gayzc1.igau 0.0) ;
  1155. SI (NON (EGA gayzu.igau 0.0)) ;
  1156. SI (EGA IMESS 1) ;
  1157. Message 'Erreur absolue' gayzu.igau
  1158. 'Composante ' (EXTR LDEFO icomp)
  1159. ' Point Gauss' igau 'Etat' ietat ;
  1160. Message ;
  1161. FINSI ;
  1162. ************ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  1163. FINSI ;
  1164. SINON ;
  1165. dgayz = ABS (gayzc1.igau - gayzu.igau) ;
  1166. errrel = ABS (dgayz / gayzc1.igau) ;
  1167. SI (errrel > critere) ;
  1168. SI (EGA IMESS 1) ;
  1169. Message 'Erreur relative' errrel
  1170. 'Composante ' (EXTR LDEFO icomp)
  1171. ' Point Gauss' igau 'Etat' ietat ;
  1172. Message 'gayzc1.igau =' gayzc1.igau 'gayzu.igau ='
  1173. gayzu.igau;
  1174. Message ;
  1175. FINSI ;
  1176. ************ierr = 21 ; QUITTER BECART ;
  1177. FINSI ;
  1178. FINSI ;
  1179. FIN BGAYZ ;
  1180.  
  1181. * Comparaison des deformations elastiques------------------------------
  1182.  
  1183. icomp = 1 ;
  1184. igau = 0 ;
  1185. REPETER BEEXX nbgau ;
  1186. igau = igau + 1 ;
  1187. SI (EGA eexxc1.igau 0.0) ;
  1188. SI (NON (EGA eexxu.igau 0.0)) ;
  1189. SI (EGA IMESS 1) ;
  1190. Message 'Erreur absolue' eexxu.igau
  1191. 'Composante ' (EXTR LDEFE icomp)
  1192. ' Point Gauss' igau 'Etat' ietat ;
  1193. Message ;
  1194. FINSI ;
  1195. ************ierr = 30 ; errabs = eexxu.igau ; QUITTER BECART ;
  1196. FINSI ;
  1197. SINON ;
  1198. deexx = ABS (eexxc1.igau - eexxu.igau) ;
  1199. errrel = ABS (deexx / eexxc1.igau) ;
  1200. SI (errrel > critere) ;
  1201. SI (EGA IMESS 1) ;
  1202. Message 'Erreur relative' errrel
  1203. 'Composante ' (EXTR LDEFE icomp)
  1204. ' Point Gauss' igau 'Etat' ietat ;
  1205. Message 'eexxc1.igau =' eexxc1.igau 'eexxu.igau ='
  1206. eexxu.igau;
  1207. Message ;
  1208. FINSI ;
  1209. ************ierr = 31 ; QUITTER BECART ;
  1210. FINSI ;
  1211. FINSI ;
  1212. FIN BEEXX ;
  1213.  
  1214. icomp = 2 ;
  1215. igau = 0 ;
  1216. REPETER BEEYY nbgau ;
  1217. igau = igau + 1 ;
  1218. SI (EGA eeyyc1.igau 0.0) ;
  1219. SI (NON (EGA eeyyu.igau 0.0)) ;
  1220. SI (EGA IMESS 1) ;
  1221. Message 'Erreur absolue' eeyyu.igau
  1222. 'Composante ' (EXTR LDEFE icomp)
  1223. ' Point Gauss' igau 'Etat' ietat ;
  1224. Message ;
  1225. FINSI ;
  1226. ************ierr = 30 ; errabs = eeyyu.igau ; QUITTER BECART ;
  1227. FINSI ;
  1228. SINON ;
  1229. deeyy = ABS (eeyyc1.igau - eeyyu.igau) ;
  1230. errrel = ABS (deeyy / eeyyc1.igau) ;
  1231. SI (errrel > critere) ;
  1232. SI (EGA IMESS 1) ;
  1233. Message 'Erreur relative' errrel
  1234. 'Composante ' (EXTR LDEFE icomp)
  1235. ' Point Gauss' igau 'Etat' ietat ;
  1236. Message 'eeyyc1.igau =' eeyyc1.igau 'eeyyu.igau ='
  1237. eeyyu.igau;
  1238. Message ;
  1239. FINSI ;
  1240. ************ierr = 31 ; QUITTER BECART ;
  1241. FINSI ;
  1242. FINSI ;
  1243. FIN BEEYY ;
  1244.  
  1245. icomp = 3 ;
  1246. igau = 0 ;
  1247. REPETER BEEZZ nbgau ;
  1248. igau = igau + 1 ;
  1249. SI (EGA eezzc1.igau 0.0) ;
  1250. SI (NON (EGA eezzu.igau 0.0)) ;
  1251. ************SI (EGA IMESS 1) ;
  1252. ************Message 'Erreur absolue' eezzu.igau
  1253. ************ 'Composante ' (EXTR LDEFE icomp)
  1254. ************ ' Point Gauss' igau 'Etat' ietat ;
  1255. ************Message ;
  1256. ************FINSI ;
  1257. ierr = 30 ; errabs = eezzu.igau ; QUITTER BECART ;
  1258. FINSI ;
  1259. SINON ;
  1260. deezz = ABS (eezzc1.igau - eezzu.igau) ;
  1261. errrel = ABS (deezz / eezzc1.igau) ;
  1262. SI (errrel > critere) ;
  1263. ************SI (EGA IMESS 1) ;
  1264. ************Message 'Erreur relative' errrel
  1265. ************ 'Composante ' (EXTR LDEFE icomp)
  1266. ************ ' Point Gauss' igau 'Etat' ietat ;
  1267. ************Message 'eezzc1.igau =' eezzc1.igau 'eezzu.igau ='
  1268. ************ eezzu.igau;
  1269. ************Message ;
  1270. ************FINSI ;
  1271. ierr = 31 ; QUITTER BECART ;
  1272. FINSI ;
  1273. FINSI ;
  1274. FIN BEEZZ ;
  1275.  
  1276. icomp = 4 ;
  1277. igau = 0 ;
  1278. REPETER BGEXY nbgau ;
  1279. igau = igau + 1 ;
  1280. SI (EGA gexyc1.igau 0.0) ;
  1281. SI (NON (EGA gexyu.igau 0.0)) ;
  1282. SI (EGA IMESS 1) ;
  1283. Message 'Erreur absolue' gexyu.igau
  1284. 'Composante ' (EXTR LDEFE icomp)
  1285. ' Point Gauss' igau 'Etat' ietat ;
  1286. Message ;
  1287. FINSI ;
  1288. ************ierr = 30 ; errabs = gexyu.igau ; QUITTER BECART ;
  1289. FINSI ;
  1290. SINON ;
  1291. dgexy = ABS (gexyc1.igau - gexyu.igau) ;
  1292. errrel = ABS (dgexy / gexyc1.igau) ;
  1293. SI (errrel > critere) ;
  1294. SI (EGA IMESS 1) ;
  1295. Message 'Erreur relative' errrel
  1296. 'Composante ' (EXTR LDEFE icomp)
  1297. ' Point Gauss' igau 'Etat' ietat ;
  1298. Message 'gexyc1.igau =' gexyc1.igau 'gexyu.igau ='
  1299. gexyu.igau;
  1300. Message ;
  1301. FINSI ;
  1302. ************ierr = 31 ; QUITTER BECART ;
  1303. FINSI ;
  1304. FINSI ;
  1305. FIN BGEXY ;
  1306.  
  1307. icomp = 5 ;
  1308. igau = 0 ;
  1309. REPETER BGEXZ nbgau ;
  1310. igau = igau + 1 ;
  1311. SI (EGA gexzc1.igau 0.0) ;
  1312. SI (NON (EGA gexzu.igau 0.0)) ;
  1313. SI (EGA IMESS 1) ;
  1314. Message 'Erreur absolue' gexzu.igau
  1315. 'Composante ' (EXTR LDEFE icomp)
  1316. ' Point Gauss' igau 'Etat' ietat ;
  1317. Message ;
  1318. FINSI ;
  1319. ************ierr = 30 ; errabs = gexzu.igau ; QUITTER BECART ;
  1320. FINSI ;
  1321. SINON ;
  1322. dgexz = ABS (gexzc1.igau - gexzu.igau) ;
  1323. errrel = ABS (dgexz / gexzc1.igau) ;
  1324. SI (errrel > critere) ;
  1325. SI (EGA IMESS 1) ;
  1326. Message 'Erreur relative' errrel
  1327. 'Composante ' (EXTR LDEFE icomp)
  1328. ' Point Gauss' igau 'Etat' ietat ;
  1329. Message 'gexzc1.igau =' gexzc1.igau 'gexzu.igau ='
  1330. gexzu.igau;
  1331. Message ;
  1332. FINSI ;
  1333. ************ierr = 31 ; QUITTER BECART ;
  1334. FINSI ;
  1335. FINSI ;
  1336. FIN BGEXZ ;
  1337.  
  1338. icomp = 6 ;
  1339. igau = 0 ;
  1340. REPETER BGEYZ nbgau ;
  1341. igau = igau + 1 ;
  1342. SI (EGA geyzc1.igau 0.0) ;
  1343. SI (NON (EGA geyzu.igau 0.0)) ;
  1344. SI (EGA IMESS 1) ;
  1345. Message 'Erreur absolue' geyzu.igau
  1346. 'Composante ' (EXTR LDEFE icomp)
  1347. ' Point Gauss' igau 'Etat' ietat ;
  1348. Message ;
  1349. FINSI ;
  1350. ************ierr = 30 ; errabs = geyzu.igau ; QUITTER BECART ;
  1351. FINSI ;
  1352. SINON ;
  1353. dgeyz = ABS (geyzc1.igau - geyzu.igau) ;
  1354. errrel = ABS (dgeyz / geyzc1.igau) ;
  1355. SI (errrel > critere) ;
  1356. SI (EGA IMESS 1) ;
  1357. Message 'Erreur relative' errrel
  1358. 'Composante ' (EXTR LDEFE icomp)
  1359. ' Point Gauss' igau 'Etat' ietat ;
  1360. Message 'geyzc1.igau =' geyzc1.igau 'geyzu.igau ='
  1361. geyzu.igau;
  1362. Message ;
  1363. FINSI ;
  1364. ************ierr = 31 ; QUITTER BECART ;
  1365. FINSI ;
  1366. FINSI ;
  1367. FIN BGEYZ ;
  1368.  
  1369. * Comparaison des deformations thermiques------------------------------
  1370.  
  1371. icomp = 1 ;
  1372. igau = 0 ;
  1373. REPETER BETXX nbgau ;
  1374. igau = igau + 1 ;
  1375. SI (EGA etxxc1.igau 0.0) ;
  1376. SI (NON (EGA etxxu.igau 0.0)) ;
  1377. SI (EGA IMESS 1) ;
  1378. Message 'Erreur absolue' etxxu.igau
  1379. 'Composante ' (EXTR LDEFT icomp)
  1380. ' Point Gauss' igau 'Etat' ietat ;
  1381. Message ;
  1382. FINSI ;
  1383. ************ierr = 40 ; errabs = etxxu.igau ; QUITTER BECART ;
  1384. FINSI ;
  1385. SINON ;
  1386. detxx = ABS (etxxc1.igau - etxxu.igau) ;
  1387. errrel = ABS (detxx / etxxc1.igau) ;
  1388. SI (errrel > critere) ;
  1389. SI (EGA IMESS 1) ;
  1390. Message 'Erreur relative' errrel
  1391. 'Composante ' (EXTR LDEFT icomp)
  1392. ' Point Gauss' igau 'Etat' ietat ;
  1393. Message 'etxxc1.igau =' etxxc1.igau 'etxxu.igau ='
  1394. etxxu.igau;
  1395. Message ;
  1396. FINSI ;
  1397. ************ierr = 41 ; QUITTER BECART ;
  1398. FINSI ;
  1399. FINSI ;
  1400. FIN BETXX ;
  1401.  
  1402. icomp = 2 ;
  1403. igau = 0 ;
  1404. REPETER BETYY nbgau ;
  1405. igau = igau + 1 ;
  1406. SI (EGA etyyc1.igau 0.0) ;
  1407. SI (NON (EGA etyyu.igau 0.0)) ;
  1408. SI (EGA IMESS 1) ;
  1409. Message 'Erreur absolue' etyyu.igau
  1410. 'Composante ' (EXTR LDEFT icomp)
  1411. ' Point Gauss' igau 'Etat' ietat ;
  1412. Message ;
  1413. FINSI ;
  1414. ************ierr = 40 ; errabs = etyyu.igau ; QUITTER BECART ;
  1415. FINSI ;
  1416. SINON ;
  1417. detyy = ABS (etyyc1.igau - etyyu.igau) ;
  1418. errrel = ABS (detyy / etyyc1.igau) ;
  1419. SI (errrel > critere) ;
  1420. SI (EGA IMESS 1) ;
  1421. Message 'Erreur relative' errrel
  1422. 'Composante ' (EXTR LDEFT icomp)
  1423. ' Point Gauss' igau 'Etat' ietat ;
  1424. Message 'etyyc1.igau =' etyyc1.igau 'etyyu.igau ='
  1425. etyyu.igau;
  1426. Message ;
  1427. FINSI ;
  1428. ************ierr = 41 ; QUITTER BECART ;
  1429. FINSI ;
  1430. FINSI ;
  1431. FIN BETYY ;
  1432.  
  1433. icomp = 3 ;
  1434. igau = 0 ;
  1435. REPETER BETZZ nbgau ;
  1436. igau = igau + 1 ;
  1437. SI (EGA etzzc1.igau 0.0) ;
  1438. SI (NON (EGA etzzu.igau 0.0)) ;
  1439. ************SI (EGA IMESS 1) ;
  1440. ************Message 'Erreur absolue' etzzu.igau
  1441. ************ 'Composante ' (EXTR LDEFT icomp)
  1442. ************ ' Point Gauss' igau 'Etat' ietat ;
  1443. ************Message ;
  1444. ************FINSI ;
  1445. ierr = 40 ; errabs = etzzu.igau ; QUITTER BECART ;
  1446. FINSI ;
  1447. SINON ;
  1448. detzz = ABS (etzzc1.igau - etzzu.igau) ;
  1449. errrel = ABS (detzz / etzzc1.igau) ;
  1450. SI (errrel > critere) ;
  1451. ************SI (EGA IMESS 1) ;
  1452. ************Message 'Erreur relative' errrel
  1453. ************ 'Composante ' (EXTR LDEFT icomp)
  1454. ************ ' Point Gauss' igau 'Etat' ietat ;
  1455. ************Message 'etzzc1.igau =' etzzc1.igau 'etzzu.igau ='
  1456. ************ etzzu.igau;
  1457. ************Message ;
  1458. ************FINSI ;
  1459. ierr = 41 ; QUITTER BECART ;
  1460. FINSI ;
  1461. FINSI ;
  1462. FIN BETZZ ;
  1463.  
  1464. icomp = 4 ;
  1465. igau = 0 ;
  1466. REPETER BGTXY nbgau ;
  1467. igau = igau + 1 ;
  1468. SI (EGA gtxyc1.igau 0.0) ;
  1469. SI (NON (EGA gtxyu.igau 0.0)) ;
  1470. SI (EGA IMESS 1) ;
  1471. Message 'Erreur absolue' gtxyu.igau
  1472. 'Composante ' (EXTR LDEFT icomp)
  1473. ' Point Gauss' igau 'Etat' ietat ;
  1474. Message ;
  1475. FINSI ;
  1476. ************ierr = 40 ; errabs = gtxyu.igau ; QUITTER BECART ;
  1477. FINSI ;
  1478. SINON ;
  1479. dgtxy = ABS (gtxyc1.igau - gtxyu.igau) ;
  1480. errrel = ABS (dgtxy / gtxyc1.igau) ;
  1481. SI (errrel > critere) ;
  1482. SI (EGA IMESS 1) ;
  1483. Message 'Erreur relative' errrel
  1484. 'Composante ' (EXTR LDEFT icomp)
  1485. ' Point Gauss' igau 'Etat' ietat ;
  1486. Message 'gtxyc1.igau =' gtxyc1.igau 'gtxyu.igau ='
  1487. gtxyu.igau;
  1488. Message ;
  1489. FINSI ;
  1490. ************ierr = 41 ; QUITTER BECART ;
  1491. FINSI ;
  1492. FINSI ;
  1493. FIN BGTXY ;
  1494.  
  1495. icomp = 5 ;
  1496. igau = 0 ;
  1497. REPETER BGTXZ nbgau ;
  1498. igau = igau + 1 ;
  1499. SI (EGA gtxzc1.igau 0.0) ;
  1500. SI (NON (EGA gtxzu.igau 0.0)) ;
  1501. SI (EGA IMESS 1) ;
  1502. Message 'Erreur absolue' gtxzu.igau
  1503. 'Composante ' (EXTR LDEFT icomp)
  1504. ' Point Gauss' igau 'Etat' ietat ;
  1505. Message ;
  1506. FINSI ;
  1507. ************ierr = 40 ; errabs = gtxzu.igau ; QUITTER BECART ;
  1508. FINSI ;
  1509. SINON ;
  1510. dgtxz = ABS (gtxzc1.igau - gtxzu.igau) ;
  1511. errrel = ABS (dgtxz / gtxzc1.igau) ;
  1512. SI (errrel > critere) ;
  1513. SI (EGA IMESS 1) ;
  1514. Message 'Erreur relative' errrel
  1515. 'Composante ' (EXTR LDEFT icomp)
  1516. ' Point Gauss' igau 'Etat' ietat ;
  1517. Message 'gtxzc1.igau =' gtxzc1.igau 'gtxzu.igau ='
  1518. gtxzu.igau;
  1519. Message ;
  1520. FINSI ;
  1521. ************ierr = 41 ; QUITTER BECART ;
  1522. FINSI ;
  1523. FINSI ;
  1524. FIN BGTXZ ;
  1525.  
  1526. icomp = 6 ;
  1527. igau = 0 ;
  1528. REPETER BGTYZ nbgau ;
  1529. igau = igau + 1 ;
  1530. SI (EGA gtyzc1.igau 0.0) ;
  1531. SI (NON (EGA gtyzu.igau 0.0)) ;
  1532. SI (EGA IMESS 1) ;
  1533. Message 'Erreur absolue' gtyzu.igau
  1534. 'Composante ' (EXTR LDEFT icomp)
  1535. ' Point Gauss' igau 'Etat' ietat ;
  1536. Message ;
  1537. FINSI ;
  1538. ************ierr = 40 ; errabs = gtyzu.igau ; QUITTER BECART ;
  1539. FINSI ;
  1540. SINON ;
  1541. dgtyz = ABS (gtyzc1.igau - gtyzu.igau) ;
  1542. errrel = ABS (dgtyz / gtyzc1.igau) ;
  1543. SI (errrel > critere) ;
  1544. SI (EGA IMESS 1) ;
  1545. Message 'Erreur relative' errrel
  1546. 'Composante ' (EXTR LDEFT icomp)
  1547. ' Point Gauss' igau 'Etat' ietat ;
  1548. Message 'gtyzc1.igau =' gtyzc1.igau 'gtyzu.igau ='
  1549. gtyzu.igau;
  1550. Message ;
  1551. FINSI ;
  1552. ************ierr = 41 ; QUITTER BECART ;
  1553. FINSI ;
  1554. FINSI ;
  1555. FIN BGTYZ ;
  1556.  
  1557. * Comparaison des contraintes------------------------------------------
  1558.  
  1559. icomp = 1 ;
  1560. igau = 0 ;
  1561. REPETER BSMXX1 nbgau ;
  1562. igau = igau + 1 ;
  1563. SI (EGA smxxc1.igau 0.0) ;
  1564. SI (NON (EGA smxxu.igau 0.0)) ;
  1565. ************SI (EGA IMESS 1) ;
  1566. ************Message 'Erreur absolue' smxxu.igau
  1567. ************ 'par rapport a la resolution quasi-statique' ;
  1568. ************Message 'Composante ' (EXTR LCONT icomp)
  1569. ************ ' Point Gauss' igau 'Etat' ietat ;
  1570. ************Message ;
  1571. ************FINSI ;
  1572. ierr = 50 ; errabs = smxxu.igau ; QUITTER BECART ;
  1573. FINSI ;
  1574. SINON ;
  1575. dsmxx = ABS (smxxc1.igau - smxxu.igau) ;
  1576. errrel = ABS (dsmxx / smxxc1.igau) ;
  1577. SI (errrel > critere) ;
  1578. ************SI (EGA IMESS 1) ;
  1579. ************Message 'Erreur relative' errrel
  1580. ************ 'par rapport a la resolution quasi-statique' ;
  1581. ************Message 'Composante ' (EXTR LCONT icomp)
  1582. ************ ' Point Gauss' igau 'Etat' ietat ;
  1583. ************Message 'smxxc1.igau =' smxxc1.igau 'smxxu.igau ='
  1584. ************ smxxu.igau;
  1585. ************Message ;
  1586. ************FINSI ;
  1587. ierr = 51 ; QUITTER BECART ;
  1588. FINSI ;
  1589. FINSI ;
  1590. FIN BSMXX1 ;
  1591. igau = 0 ;
  1592. REPETER BSMXX2 nbgau ;
  1593. igau = igau + 1 ;
  1594. SI (EGA smxxc2.igau 0.0) ;
  1595. SI (NON (EGA smxxu.igau 0.0)) ;
  1596. ************SI (EGA IMESS 1) ;
  1597. ************Message 'Erreur absolue' smxxu.igau
  1598. ************ 'par rapport a la resolution par PASAPAS' ;
  1599. ************Message 'Composante ' (EXTR LCONT icomp)
  1600. ************ ' Point Gauss' igau 'Etat' ietat ;
  1601. ************Message ;
  1602. ************FINSI ;
  1603. ierr = 52 ; errabs = smxxu.igau ; QUITTER BECART ;
  1604. FINSI ;
  1605. SINON ;
  1606. dsmxx = ABS (smxxc2.igau - smxxu.igau) ;
  1607. errrel = ABS (dsmxx / smxxc2.igau) ;
  1608. SI (errrel > critere) ;
  1609. ************SI (EGA IMESS 1) ;
  1610. ************Message 'Erreur relative' errrel
  1611. ************ 'par rapport a la resolution par PASAPAS' ;
  1612. ************Message 'Composante ' (EXTR LCONT icomp)
  1613. ************ ' Point Gauss' igau 'Etat' ietat ;
  1614. ************Message 'smxxc2.igau =' smxxc2.igau 'smxxu.igau ='
  1615. ************ smxxu.igau;
  1616. ************Message ;
  1617. ************FINSI ;
  1618. ierr = 53 ; QUITTER BECART ;
  1619. FINSI ;
  1620. FINSI ;
  1621. FIN BSMXX2 ;
  1622.  
  1623.  
  1624. icomp = 2 ;
  1625. igau = 0 ;
  1626. REPETER BSMYY1 nbgau ;
  1627. igau = igau + 1 ;
  1628. SI (EGA smyyc1.igau 0.0) ;
  1629. SI (NON (EGA smyyu.igau 0.0)) ;
  1630. ************SI (EGA IMESS 1) ;
  1631. ************Message 'Erreur absolue' smyyu.igau
  1632. ************ 'par rapport a la resolution quasi-statique' ;
  1633. ************Message 'Composante ' (EXTR LCONT icomp)
  1634. ************ ' Point Gauss' igau 'Etat' ietat ;
  1635. ************Message ;
  1636. ************FINSI ;
  1637. ierr = 50 ; errabs = smyyu.igau ; QUITTER BECART ;
  1638. FINSI ;
  1639. SINON ;
  1640. dsmyy = ABS (smyyc1.igau - smyyu.igau) ;
  1641. errrel = ABS (dsmyy / smyyc1.igau) ;
  1642. SI (errrel > critere) ;
  1643. ************SI (EGA IMESS 1) ;
  1644. ************Message 'Erreur relative' errrel
  1645. ************ 'par rapport a la resolution quasi-statique' ;
  1646. ************Message 'Composante ' (EXTR LCONT icomp)
  1647. ************ ' Point Gauss' igau 'Etat' ietat ;
  1648. ************Message 'smyyc1.igau =' smyyc1.igau 'smyyu.igau ='
  1649. ************ smyyu.igau;
  1650. ************Message ;
  1651. ************FINSI ;
  1652. ierr = 51 ; QUITTER BECART ;
  1653. FINSI ;
  1654. FINSI ;
  1655. FIN BSMYY1 ;
  1656. igau = 0 ;
  1657. REPETER BSMYY2 nbgau ;
  1658. igau = igau + 1 ;
  1659. SI (EGA smyyc2.igau 0.0) ;
  1660. SI (NON (EGA smyyu.igau 0.0)) ;
  1661. ************SI (EGA IMESS 1) ;
  1662. ************Message 'Erreur absolue' smyyu.igau
  1663. ************ 'par rapport a la resolution par PASAPAS' ;
  1664. ************Message 'Composante ' (EXTR LCONT icomp)
  1665. ************ ' Point Gauss' igau 'Etat' ietat ;
  1666. ************Message ;
  1667. ************FINSI ;
  1668. ierr = 52 ; errabs = smyyu.igau ; QUITTER BECART ;
  1669. FINSI ;
  1670. SINON ;
  1671. dsmyy = ABS (smyyc2.igau - smyyu.igau) ;
  1672. errrel = ABS (dsmyy / smyyc2.igau) ;
  1673. SI (errrel > critere) ;
  1674. ************SI (EGA IMESS 1) ;
  1675. ************Message 'Erreur relative' errrel
  1676. ************ 'par rapport a la resolution par PASAPAS' ;
  1677. ************Message 'Composante ' (EXTR LCONT icomp)
  1678. ************ ' Point Gauss' igau 'Etat' ietat ;
  1679. ************Message 'smyyc2.igau =' smyyc2.igau 'smyyu.igau ='
  1680. ************ smyyu.igau;
  1681. ************Message ;
  1682. ************FINSI ;
  1683. ierr = 53 ; QUITTER BECART ;
  1684. FINSI ;
  1685. FINSI ;
  1686. FIN BSMYY2 ;
  1687.  
  1688.  
  1689. icomp = 3 ;
  1690. igau = 0 ;
  1691. REPETER BSMZZ1 nbgau ;
  1692. igau = igau + 1 ;
  1693. SI (EGA smzzc1.igau 0.0) ;
  1694. SI (NON (EGA smzzu.igau 0.0)) ;
  1695. ************SI (EGA IMESS 1) ;
  1696. ************Message 'Erreur absolue' smzzu.igau
  1697. ************ 'par rapport a la resolution quasi-statique' ;
  1698. ************Message 'Composante ' (EXTR LCONT icomp)
  1699. ************ ' Point Gauss' igau 'Etat' ietat ;
  1700. ************Message ;
  1701. ************FINSI ;
  1702. ierr = 50 ; errabs = smzzu.igau ; QUITTER BECART ;
  1703. FINSI ;
  1704. SINON ;
  1705. dsmzz = ABS (smzzc1.igau - smzzu.igau) ;
  1706. errrel = ABS (dsmzz / smzzc1.igau) ;
  1707. SI (errrel > critere) ;
  1708. ************SI (EGA IMESS 1) ;
  1709. ************Message 'Erreur relative' errrel
  1710. ************ 'par rapport a la resolution quasi-statique' ;
  1711. ************Message 'Composante ' (EXTR LCONT icomp)
  1712. ************ ' Point Gauss' igau 'Etat' ietat ;
  1713. ************Message 'smzzc1.igau =' smzzc1.igau 'smzzu.igau ='
  1714. ************ smzzu.igau;
  1715. ************Message ;
  1716. ************FINSI ;
  1717. ierr = 51 ; QUITTER BECART ;
  1718. FINSI ;
  1719. FINSI ;
  1720. FIN BSMZZ1 ;
  1721. igau = 0 ;
  1722. REPETER BSMZZ2 nbgau ;
  1723. igau = igau + 1 ;
  1724. SI (EGA smzzc2.igau 0.0) ;
  1725. SI (NON (EGA smzzu.igau 0.0)) ;
  1726. ************SI (EGA IMESS 1) ;
  1727. ************Message 'Erreur absolue' smzzu.igau
  1728. ************ 'par rapport a la resolution par PASAPAS' ;
  1729. ************Message 'Composante ' (EXTR LCONT icomp)
  1730. ************ ' Point Gauss' igau 'Etat' ietat ;
  1731. ************Message ;
  1732. ************FINSI ;
  1733. ierr = 52 ; errabs = smzzu.igau ; QUITTER BECART ;
  1734. FINSI ;
  1735. SINON ;
  1736. dsmzz = ABS (smzzc2.igau - smzzu.igau) ;
  1737. errrel = ABS (dsmzz / smzzc2.igau) ;
  1738. SI (errrel > critere) ;
  1739. ************SI (EGA IMESS 1) ;
  1740. ************Message 'Erreur relative' errrel
  1741. ************ 'par rapport a la resolution par PASAPAS' ;
  1742. ************Message 'Composante ' (EXTR LCONT icomp)
  1743. ************ ' Point Gauss' igau 'Etat' ietat ;
  1744. ************Message 'smzzc2.igau =' smzzc2.igau 'smzzu.igau ='
  1745. ************ smzzu.igau;
  1746. ************Message ;
  1747. ************FINSI ;
  1748. ierr = 53 ; QUITTER BECART ;
  1749. FINSI ;
  1750. FINSI ;
  1751. FIN BSMZZ2 ;
  1752.  
  1753.  
  1754. icomp = 4 ;
  1755. igau = 0 ;
  1756. REPETER BSMXY1 nbgau ;
  1757. igau = igau + 1 ;
  1758. SI (EGA smxyc1.igau 0.0) ;
  1759. SI (NON (EGA smxyu.igau 0.0)) ;
  1760. SI (EGA IMESS 1) ;
  1761. Message 'Erreur absolue' smxyu.igau
  1762. 'par rapport a la resolution quasi-statique' ;
  1763. Message 'Composante ' (EXTR LCONT icomp)
  1764. ' Point Gauss' igau 'Etat' ietat ;
  1765. Message ;
  1766. FINSI ;
  1767. ************ierr = 50 ; errabs = smxyu.igau ; QUITTER BECART ;
  1768. FINSI ;
  1769. SINON ;
  1770. dsmxy = ABS (smxyc1.igau - smxyu.igau) ;
  1771. errrel = ABS (dsmxy / smxyc1.igau) ;
  1772. SI (errrel > critere) ;
  1773. SI (EGA IMESS 1) ;
  1774. Message 'Erreur relative' errrel
  1775. 'par rapport a la resolution quasi-statique' ;
  1776. Message 'Composante ' (EXTR LCONT icomp)
  1777. ' Point Gauss' igau 'Etat' ietat ;
  1778. Message 'smxyc1.igau =' smxyc1.igau 'smxyu.igau ='
  1779. smxyu.igau;
  1780. Message ;
  1781. FINSI ;
  1782. ************ierr = 51 ; QUITTER BECART ;
  1783. FINSI ;
  1784. FINSI ;
  1785. FIN BSMXY1 ;
  1786. igau = 0 ;
  1787. REPETER BSMXY2 nbgau ;
  1788. igau = igau + 1 ;
  1789. SI (EGA smxyc2.igau 0.0) ;
  1790. SI (NON (EGA smxyu.igau 0.0)) ;
  1791. SI (EGA IMESS 1) ;
  1792. Message 'Erreur absolue' smxyu.igau
  1793. 'par rapport a la resolution par PASAPAS' ;
  1794. Message 'Composante ' (EXTR LCONT icomp)
  1795. ' Point Gauss' igau 'Etat' ietat ;
  1796. Message ;
  1797. FINSI ;
  1798. ************ierr = 52 ; errabs = smxyu.igau ; QUITTER BECART ;
  1799. FINSI ;
  1800. SINON ;
  1801. dsmxy = ABS (smxyc2.igau - smxyu.igau) ;
  1802. errrel = ABS (dsmxy / smxyc2.igau) ;
  1803. SI (errrel > critere) ;
  1804. SI (EGA IMESS 1) ;
  1805. Message 'Erreur relative' errrel
  1806. 'par rapport a la resolution par PASAPAS' ;
  1807. Message 'Composante ' (EXTR LCONT icomp)
  1808. ' Point Gauss' igau 'Etat' ietat ;
  1809. Message 'smxyc2.igau =' smxyc2.igau 'smxyu.igau ='
  1810. smxyu.igau;
  1811. Message ;
  1812. FINSI ;
  1813. ************ierr = 53 ; QUITTER BECART ;
  1814. FINSI ;
  1815. FINSI ;
  1816. FIN BSMXY2 ;
  1817.  
  1818.  
  1819. icomp = 5 ;
  1820. igau = 0 ;
  1821. REPETER BSMXZ1 nbgau ;
  1822. igau = igau + 1 ;
  1823. SI (EGA smxzc1.igau 0.0) ;
  1824. SI (NON (EGA smxzu.igau 0.0)) ;
  1825. SI (EGA IMESS 1) ;
  1826. Message 'Erreur absolue' smxzu.igau
  1827. 'par rapport a la resolution quasi-statique' ;
  1828. Message 'Composante ' (EXTR LCONT icomp)
  1829. ' Point Gauss' igau 'Etat' ietat ;
  1830. Message ;
  1831. FINSI ;
  1832. ************ierr = 50 ; errabs = smxzu.igau ; QUITTER BECART ;
  1833. FINSI ;
  1834. SINON ;
  1835. dsmxz = ABS (smxzc1.igau - smxzu.igau) ;
  1836. errrel = ABS (dsmxz / smxzc1.igau) ;
  1837. SI (errrel > critere) ;
  1838. SI (EGA IMESS 1) ;
  1839. Message 'Erreur relative' errrel
  1840. 'par rapport a la resolution quasi-statique' ;
  1841. Message 'Composante ' (EXTR LCONT icomp)
  1842. ' Point Gauss' igau 'Etat' ietat ;
  1843. Message 'smxzc1.igau =' smxzc1.igau 'smxzu.igau ='
  1844. smxzu.igau;
  1845. Message ;
  1846. FINSI ;
  1847. ************ierr = 51 ; QUITTER BECART ;
  1848. FINSI ;
  1849. FINSI ;
  1850. FIN BSMXZ1 ;
  1851. igau = 0 ;
  1852. REPETER BSMXZ2 nbgau ;
  1853. igau = igau + 1 ;
  1854. SI (EGA smxzc2.igau 0.0) ;
  1855. SI (NON (EGA smxzu.igau 0.0)) ;
  1856. SI (EGA IMESS 1) ;
  1857. Message 'Erreur absolue' smxzu.igau
  1858. 'par rapport a la resolution par PASAPAS' ;
  1859. Message 'Composante ' (EXTR LCONT icomp)
  1860. ' Point Gauss' igau 'Etat' ietat ;
  1861. Message ;
  1862. FINSI ;
  1863. ************ierr = 52 ; errabs = smxzu.igau ; QUITTER BECART ;
  1864. FINSI ;
  1865. SINON ;
  1866. dsmxz = ABS (smxzc2.igau - smxzu.igau) ;
  1867. errrel = ABS (dsmxz / smxzc2.igau) ;
  1868. SI (errrel > critere) ;
  1869. SI (EGA IMESS 1) ;
  1870. Message 'Erreur relative' errrel
  1871. 'par rapport a la resolution par PASAPAS' ;
  1872. Message 'Composante ' (EXTR LCONT icomp)
  1873. ' Point Gauss' igau 'Etat' ietat ;
  1874. Message 'smxzc2.igau =' smxzc2.igau 'smxzu.igau ='
  1875. smxzu.igau;
  1876. Message ;
  1877. FINSI ;
  1878. ************ierr = 53 ; QUITTER BECART ;
  1879. FINSI ;
  1880. FINSI ;
  1881. FIN BSMXZ2 ;
  1882.  
  1883.  
  1884. icomp = 6 ;
  1885. igau = 0 ;
  1886. REPETER BSMYZ1 nbgau ;
  1887. igau = igau + 1 ;
  1888. SI (EGA smyzc1.igau 0.0) ;
  1889. SI (NON (EGA smyzu.igau 0.0)) ;
  1890. SI (EGA IMESS 1) ;
  1891. Message 'Erreur absolue' smyzu.igau
  1892. 'par rapport a la resolution quasi-statique' ;
  1893. Message 'Composante ' (EXTR LCONT icomp)
  1894. ' Point Gauss' igau 'Etat' ietat ;
  1895. Message ;
  1896. FINSI ;
  1897. ************ierr = 50 ; errabs = smyzu.igau ; QUITTER BECART ;
  1898. FINSI ;
  1899. SINON ;
  1900. dsmyz = ABS (smyzc1.igau - smyzu.igau) ;
  1901. errrel = ABS (dsmyz / smyzc1.igau) ;
  1902. SI (errrel > critere) ;
  1903. SI (EGA IMESS 1) ;
  1904. Message 'Erreur relative' errrel
  1905. 'par rapport a la resolution quasi-statique' ;
  1906. Message 'Composante ' (EXTR LCONT icomp)
  1907. ' Point Gauss' igau 'Etat' ietat ;
  1908. Message 'smyzc1.igau =' smyzc1.igau 'smyzu.igau ='
  1909. smyzu.igau;
  1910. Message ;
  1911. FINSI ;
  1912. ************ierr = 51 ; QUITTER BECART ;
  1913. FINSI ;
  1914. FINSI ;
  1915. FIN BSMYZ1 ;
  1916. igau = 0 ;
  1917. REPETER BSMYZ2 nbgau ;
  1918. igau = igau + 1 ;
  1919. SI (EGA smyzc2.igau 0.0) ;
  1920. SI (NON (EGA smyzu.igau 0.0)) ;
  1921. SI (EGA IMESS 1) ;
  1922. Message 'Erreur absolue' smyzu.igau
  1923. 'par rapport a la resolution par PASAPAS' ;
  1924. Message 'Composante ' (EXTR LCONT icomp)
  1925. ' Point Gauss' igau 'Etat' ietat ;
  1926. Message ;
  1927. FINSI ;
  1928. ************ierr = 52 ; errabs = smyzu.igau ; QUITTER BECART ;
  1929. FINSI ;
  1930. SINON ;
  1931. dsmyz = ABS (smyzc2.igau - smyzu.igau) ;
  1932. errrel = ABS (dsmyz / smyzc2.igau) ;
  1933. SI (errrel > critere) ;
  1934. SI (EGA IMESS 1) ;
  1935. Message 'Erreur relative' errrel
  1936. 'par rapport a la resolution par PASAPAS' ;
  1937. Message 'Composante ' (EXTR LCONT icomp)
  1938. ' Point Gauss' igau 'Etat' ietat ;
  1939. Message 'smyzc2.igau =' smyzc2.igau 'smyzu.igau ='
  1940. smyzu.igau;
  1941. Message ;
  1942. FINSI ;
  1943. ************ierr = 53 ; QUITTER BECART ;
  1944. FINSI ;
  1945. FINSI ;
  1946. FIN BSMYZ2 ;
  1947.  
  1948.  
  1949. FIN BECART ;
  1950.  
  1951. * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  1952.  
  1953. * MESSAGES D'ERREUR_____________________________________________________
  1954.  
  1955. SI (ierr > 0) ;
  1956.  
  1957. SI (EGA ierr 1) ;
  1958. Message 'Instants de calcul non coincidants' ;
  1959. ERRE 5 ;
  1960. FINSI ;
  1961.  
  1962. SI (EGA ierr 10) ;
  1963. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1964. ' Noeud' ino 'Etat' ietat ;
  1965. Message 'Erreur par rapport a la resolution quasi-statique' ;
  1966. ERRE 5 ;
  1967. FINSI ;
  1968.  
  1969. SI (EGA ierr 11) ;
  1970. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1971. ' Noeud' ino 'Etat' ietat ;
  1972. Message 'Erreur par rapport a la resolution quasi-statique' ;
  1973. ERRE 5 ;
  1974. FINSI ;
  1975.  
  1976. SI (EGA ierr 12) ;
  1977. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1978. ' Noeud' ino 'Etat' ietat ;
  1979. Message 'Erreur par rapport a la resolution par PASAPAS' ;
  1980. ERRE 5 ;
  1981. FINSI ;
  1982.  
  1983. SI (EGA ierr 13) ;
  1984. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1985. ' Noeud' ino 'Etat' ietat ;
  1986. Message 'Erreur par rapport a la resolution par PASAPAS' ;
  1987. ERRE 5 ;
  1988. FINSI ;
  1989.  
  1990. SI (EGA ierr 20) ;
  1991. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  1992. ' Point Gauss' igau 'Etat' ietat ;
  1993. ERRE 5 ;
  1994. FINSI ;
  1995.  
  1996. SI (EGA ierr 21) ;
  1997. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  1998. ' Point Gauss' igau 'Etat' ietat ;
  1999. ERRE 5 ;
  2000. FINSI ;
  2001.  
  2002. SI (EGA ierr 30) ;
  2003. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFE icomp)
  2004. ' Point Gauss' igau 'Etat' ietat ;
  2005. ERRE 5 ;
  2006. FINSI ;
  2007.  
  2008. SI (EGA ierr 31) ;
  2009. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFE icomp)
  2010. ' Point Gauss' igau 'Etat' ietat ;
  2011. ERRE 5 ;
  2012. FINSI ;
  2013.  
  2014. SI (EGA ierr 40) ;
  2015. Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFT icomp)
  2016. ' Point Gauss' igau 'Etat' ietat ;
  2017. ERRE 5 ;
  2018. FINSI ;
  2019.  
  2020. SI (EGA ierr 41) ;
  2021. Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFT icomp)
  2022. ' Point Gauss' igau 'Etat' ietat ;
  2023. ERRE 5 ;
  2024. FINSI ;
  2025.  
  2026. SI (EGA ierr 50) ;
  2027. Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  2028. ' Point Gauss' igau 'Etat' ietat ;
  2029. Message 'Erreur par rapport a la resolution quasi-statique' ;
  2030. ERRE 5 ;
  2031. FINSI ;
  2032.  
  2033. SI (EGA ierr 51) ;
  2034. Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  2035. ' Point Gauss' igau 'Etat' ietat ;
  2036. Message 'Erreur par rapport a la resolution quasi-statique' ;
  2037. ERRE 5 ;
  2038. FINSI ;
  2039.  
  2040. SI (EGA ierr 52) ;
  2041. Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  2042. ' Point Gauss' igau 'Etat' ietat ;
  2043. Message 'Erreur par rapport a la resolution par PASAPAS' ;
  2044. ERRE 5 ;
  2045. FINSI ;
  2046.  
  2047. SI (EGA ierr 53) ;
  2048. Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  2049. ' Point Gauss' igau 'Etat' ietat ;
  2050. Message 'Erreur par rapport a la resolution par PASAPAS' ;
  2051. ERRE 5 ;
  2052. FINSI ;
  2053.  
  2054. SINON ;
  2055.  
  2056. Message 'Resultats concordants entre les deux modeles' ;
  2057. Message 'Aucune erreur relative superieure a ' critere
  2058. ' pour les composantes' ;
  2059. Message (EXTR LDEFO idefo) ' ' (EXTR LDEFE idefe) ' '
  2060. (EXTR LDEFT ideft) ' ' (EXTR LCONT icon1) ' '
  2061. (EXTR LCONT icon2) ' ' (EXTR LCONT icon3) ;
  2062. ERRE 0 ;
  2063.  
  2064. FINSI ;
  2065.  
  2066. FIN ;
  2067.  
  2068.  
  2069.  
  2070.  
  2071.  
  2072.  
  2073.  
  2074.  
  2075.  
  2076.  

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