Télécharger la plaquette Cast3M


fr en

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

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