Télécharger la plaquette Cast3M


fr en

Télécharger umat03_cisXY_2122.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  umat03_cisXY_2122.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 No3 : Modele de FLUAGE POLYNOMIAL
  15. *            CISAILLEMENT XY
  16. *
  17. *=======================================================================
  18.  
  19.  
  20. opti echo 0 ;
  21.  
  22.  
  23. *=======================================================================
  24. *
  25. *                       PROCEDURES DEPENDANTES
  26. *
  27. *=======================================================================
  28.  
  29.  
  30.  
  31.  
  32. *=======================================================================
  33. *
  34. *                      JEU DE DONNEES PRINCIPAL
  35. *
  36. *=======================================================================
  37.  
  38. *-----------------------------------------------------------------------
  39. * 0 - Maillage
  40. *     PLAQUE CARREE    COTE : 1 m    EPAISSEUR : 1 cm
  41.  
  42. opti dime 3 elem cub8 mode tridim ;
  43.  
  44. l1 = 1.0 ; 
  45. l2 = 1.0 ; 
  46. ep0 = 0.01 ; 
  47.  
  48. p01 = 0.0 0.0 0.0 ; 
  49. p02 = ep0 0.0 0.0 ; 
  50. p03 = ep0 l2  0.0 ; 
  51. p04 = 0.0 l2  0.0 ; 
  52.  
  53. p10 = 0.0 0.0 l1 ;
  54. p20 = ep0 0.0 l1 ;
  55. p30 = ep0 l2  l1 ;
  56. p40 = 0.0 l2  l1 ; 
  57.  
  58. li01 = 'DROIT' 1 p01 p02 ; 
  59. li02 = 'DROIT' 1 p02 p03 ;
  60. li03 = 'DROIT' 1 p03 p04 ;
  61. li04 = 'DROIT' 1 p04 p01 ;
  62.  
  63. li10 = 'DROIT' 1 p10 p20 ;
  64. li20 = 'DROIT' 1 p20 p30 ;
  65. li30 = 'DROIT' 1 p30 p40 ;
  66. li40 = 'DROIT' 1 p40 p10 ;
  67.  
  68. sur01 = 'DALL' li01 li02 li03 li04 ;
  69. sur10 = 'DALL' li10 li20 li30 li40 ;
  70. mail1 = sur01 'VOLU' 1 sur10 ;
  71.  
  72. li2  = 'DROIT' 1 p02 p20 ;
  73. li3  = 'DROIT' 1 p20 p10 ;
  74. li4  = 'DROIT' 1 p10 p01 ;
  75. sbloq = 'DALL' li01 li2 li3 li4 ;
  76.  
  77. li5  = 'DROIT' 1 p04 p03 ;
  78. li6  = 'DROIT' 1 p03 p30 ;
  79. li7  = 'DROIT' 1 p30 p40 ;
  80. li8  = 'DROIT' 1 p40 p04 ;
  81. schar = 'DALL' li5 li6 li7 li8 ;
  82.  
  83. ELIM (mail1 ET sbloq ET schar) 1.e-5 ;
  84.  
  85. *-----------------------------------------------------------------------
  86. * 1 - Calcul avec modele utilisateur no21
  87. *     Modele 'FLUAGE' 'POLYNOMIAL' externalise
  88.  
  89. * 1.1 Definition du modele utilisateur
  90.  
  91. LCMAT21 = MOTS 'YOUN' 'NU  ' 'RHO ' 'ALPH'
  92.                'SMAX' 'AF0 ' 'AF1 ' 'AF2 ' 'AF3 ' 'AF4 ' 'AF5 ' 'AF6 ';
  93.  
  94. LCVAR21 = MOTS 'PXX ' 'PYY ' 'PZZ ' 'PXY ' 'PYZ ' 'PZX '
  95.                'VNXX' 'VNYY' 'VNZZ' 'VNXY' 'VNYZ' 'VNZX'
  96.                'EPCE' 'EPSH' 'LL  ' 'EPSE'
  97.                'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  98.  
  99. mo_uti1 = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  100.                      'NON_LINEAIRE' 'UTILISATEUR'
  101.                      'NUME_LOI' 21
  102.                      'C_MATERIAU' LCMAT21 'C_VARINTER' LCVAR21 ;
  103.  
  104. * 1.2 Affectation des proprietes de materiau
  105.  
  106. densiu   = 0.95 ;
  107. dgrain   = 8.0E-6 ;
  108. tuo2     = 1500.0 ;
  109. tauxfisk = 1.5E19 ;
  110.  
  111. xyoun = (2.2693E11 - (1.5399E7*tuo2)) - (9.35973E3 * (tuo2**2.0)) ;
  112. xcis  = (8.583E10 - (5.157E6*tuo2)) - (3.747E3 * (tuo2**2.0)) ;
  113. xnu   = (xyoun/(2.0*xcis)) - 1.0 ;
  114. xrho  = 10950.0 ;
  115. xalph = ( (1.1833E-5 - (5.013E-9*tuo2)) + (3.756E-12 * (tuo2**2.0)) )
  116.         - (6.125E-17 * (tuo2**3.0)) ;
  117.  
  118. xsmax = xyoun / 1000.0 ;
  119. xaf0  = 0.0 ;
  120. xaf1  = (EXP(0.183*(100.0-densiu)))
  121.       * ( (0.2031E-13+(0.67861E-32*tauxfisk)) / (dgrain**2.0) )
  122.       * (EXP(-45294.4*(tuo2**-1.0))) ;
  123. xaf2  = 1.0 ;
  124. xaf3  = (EXP(0.241*(100.0-densiu)))
  125.       * (0.1524E-25+(0.50825E-44*tauxfisk))
  126.       * (EXP(-66431.8*(tuo2**-1.0))) ;
  127. xaf4  = 4.5 ;
  128. xaf5  = (1.447E-34*tauxfisk) * (EXP(-3624.0*(tuo2**-1.0))) ;
  129. xaf6  = 1.0 ;
  130.  
  131. ma_uti1 = MATE mo_uti1 'YOUN' xyoun 'NU  ' xnu   'RHO ' xrho
  132.                        'ALPH' xalph
  133.                        'SMAX' xsmax 'AF0 ' xaf0  'AF1 ' xaf1
  134.                        'AF2 ' xaf2  'AF3 ' xaf3  'AF4 ' xaf4
  135.                        'AF5 ' xaf5  'AF6 ' xaf6  ;
  136.  
  137. * 1.3 Definition des conditions aux limites et chargements
  138.  
  139. cl1 = BLOQ mail1 'UY' 'UZ' ;
  140. cl2 = BLOQ sbloq 'UX' ;
  141. cl0 = cl1 ET cl2 ;
  142.  
  143. ********depmax = 5.0E5 ;
  144. depmax = 5.0E8 ;
  145. dep0 = FORCE schar 'FX' depmax 'FY' 0.0 'FZ' 0.0 ;
  146. ev0 = EVOL MANU temps (PROG 0.0 1.E-2 1000.0)
  147.                     y (PROG 0.0 1.0   1.0   ) ;
  148. cha0 = CHAR MECA dep0 ev0 ;
  149.  
  150. * 1.4 Definition des instants de calcul
  151.  
  152. pro0 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  153. pro1 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  154.  
  155. * 1.5 Initialisation des variables internes
  156.  
  157. chpxx01  = MANU CHML mo_uti1 'PXX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  158. chpyy01  = MANU CHML mo_uti1 'PYY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  159. chpzz01  = MANU CHML mo_uti1 'PZZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  160. chpxy01  = MANU CHML mo_uti1 'PXY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  161. chpyz01  = MANU CHML mo_uti1 'PYZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  162. chpzx01  = MANU CHML mo_uti1 'PZX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  163.  
  164. chvnxx01 = MANU CHML mo_uti1 'VNXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  165. chvnyy01 = MANU CHML mo_uti1 'VNYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  166. chvnzz01 = MANU CHML mo_uti1 'VNZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  167. chvnxy01 = MANU CHML mo_uti1 'VNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  168. chvnyz01 = MANU CHML mo_uti1 'VNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  169. chvnzx01 = MANU CHML mo_uti1 'VNZX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  170.  
  171. chepce01 = MANU CHML mo_uti1 'EPCE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  172. chepsh01 = MANU CHML mo_uti1 'EPSH' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  173. chll01   = MANU CHML mo_uti1 'LL  ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  174. chepse01 = MANU CHML mo_uti1 'EPSE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  175.  
  176. chenxx01 = MANU CHML mo_uti1 'ENXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  177. chenyy01 = MANU CHML mo_uti1 'ENYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  178. chenzz01 = MANU CHML mo_uti1 'ENZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  179. chgnxy01 = MANU CHML mo_uti1 'GNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  180. chgnxz01 = MANU CHML mo_uti1 'GNXZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  181. chgnyz01 = MANU CHML mo_uti1 'GNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  182.  
  183. var_ut01 = chpxx01  ET chpyy01  ET chpzz01  ET
  184.            chpxy01  ET chpyz01  ET chpzx01  ET
  185.            chvnxx01 ET chvnyy01 ET chvnzz01 ET
  186.            chvnxy01 ET chvnyz01 ET chvnzx01 ET
  187.            chepce01 ET chepsh01 ET chll01   ET chepse01 ET
  188.            chenxx01 ET chenyy01 ET chenzz01 ET
  189.            chgnxy01 ET chgnxz01 ET chgnyz01 ;
  190.  
  191. * 1.6 Resolution par PASAPAS
  192.  
  193. tabuti1 = TABLE ;
  194. tabuti1.'MODELE' = mo_uti1 ;
  195. tabuti1.'CARACTERISTIQUES' = ma_uti1 ;
  196. tabuti1.'CHARGEMENT' = cha0 ;
  197. tabuti1.'BLOCAGES_MECANIQUES' = cl0 ;
  198. tabuti1.'VARIABLES_INTERNES' = TABLE ;
  199. tabuti1.'VARIABLES_INTERNES' . 0 = var_ut01 ;
  200. tabuti1.'TEMPS_CALCULES' = pro0 ;
  201. tabuti1.'TEMPS_SAUVES' = pro1 ;
  202. tabuti1.'PRECISION' = 1E-10;
  203. TMASAU=table;
  204. tabuti1 . 'MES_SAUVEGARDES'=TMASAU;
  205. TMASAU .'DEFTO'=VRAI;
  206. TMASAU .'DEFIN'=VRAI;
  207. TEMPS 'ZERO' ;
  208. PASAPAS tabuti1 ;
  209. TABTPS = TEMP 'NOEC';
  210. CPUext1 =  TABTPS.'TEMPS_CPU'.'INITIAL';
  211.  
  212. *-----------------------------------------------------------------------
  213. * 2 - Calcul avec modele utilisateur no22
  214. *     Modele 'FLUAGE' 'POLYNOMIAL' instanciable au niveau
  215. *     des iterations internes
  216.  
  217. * 2.1 Definition du modele utilisateur
  218.  
  219. LCMAT22 = MOTS 'YOUN' 'NU  ' 'RHO ' 'ALPH'
  220.                'SMAX' 'AF0 ' 'AF1 ' 'AF2 ' 'AF3 ' 'AF4 ' 'AF5 ' 'AF6 ' 'TALP' 'TREF' ;
  221.  
  222. LCVAR22 = MOTS 'PXX ' 'PYY ' 'PZZ ' 'PXY ' 'PYZ ' 'PZX '
  223.                'VNXX' 'VNYY' 'VNZZ' 'VNXY' 'VNYZ' 'VNZX'
  224.                'EPCE' 'EPSH' 'LL  ' 'EPSE'
  225.                'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  226.  
  227. LCPAR22 = MOTS 'T   ' 'PORO' 'YOGC' 'ALPC'
  228.                'TFIS' 'TUO2' 'FACF' 'DSIU' 'DGRA' ;
  229.  
  230. mo_uti2 = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  231.                      'NON_LINEAIRE' 'UTILISATEUR'
  232.                      'NUME_LOI' 22 'PARA_LOI' LCPAR22
  233.                      'C_MATERIAU' LCMAT22 'C_VARINTER' LCVAR22 ;
  234.  
  235. * 2.2 Affectation des proprietes de materiau
  236.  
  237. lpyoun = MOTS 'T   ' 'PORO' 'YOGC' ;
  238. lpnu   = MOTS 'T   ' 'PORO' 'YOGC' ;
  239. lpalph = MOTS 'T   ' 'ALPC' ;
  240. lpsmax = MOTS 'T   ' 'PORO' 'YOGC' ;
  241. lpaf0  = MOTS 'TFIS' 'TUO2' 'FACF' ;
  242. lpaf1  = MOTS 'TFIS' 'TUO2' 'FACF' 'DSIU' 'DGRA' ;
  243. lpaf2  = MOTS 'TFIS' 'TUO2' 'FACF' ;
  244. lpaf3  = MOTS 'TFIS' 'TUO2' 'FACF' 'DSIU' ;
  245. lpaf4  = MOTS 'TFIS' 'TUO2' 'FACF' ;
  246. lpaf5  = MOTS 'TFIS' 'TUO2' 'FACF' ;
  247. lpaf6  = MOTS 'TFIS' 'TUO2' 'FACF' ;
  248.  
  249. ma_uti2 = MATE mo_uti2 'YOUN' lpyoun 'NU  ' lpnu   'RHO ' xrho
  250.                        'ALPH' lpalph
  251.                        'SMAX' lpsmax 'AF0 ' lpaf0  'AF1 ' lpaf1
  252.                        'AF2 ' lpaf2  'AF3 ' lpaf3  'AF4 ' lpaf4 
  253.                        'AF5 ' lpaf5  'AF6 ' lpaf6  'TALP' 0. 'TREF' tuo2;
  254.  
  255. * 2.3 Definition des conditions aux limites et chargements
  256. *     IDEM 1.3 pour les CL et le chargement mecaniques
  257. *     Chargements supplementaires correspondant aux parametres externes
  258.  
  259. chT    = MANU CHPO mail1 'T   ' tuo2     ;
  260. chPORO = MANU CHPO mail1 'PORO' 0.0      ;
  261. chYOGC = MANU CHPO mail1 'YOGC' 1.0      ;
  262. chALPC = MANU CHPO mail1 'ALPC' 1.0      ;
  263. chTFIS = MANU CHPO mail1 'TFIS' tauxfisk ;
  264. chTUO2 = MANU CHPO mail1 'TUO2' tuo2     ;
  265. chFACF = MANU CHPO mail1 'FACF' 1.0      ;
  266. chDSIU = MANU CHPO mail1 'DSIU' densiu   ;
  267. chDGRA = MANU CHPO mail1 'DGRA' dgrain   ;
  268.  
  269. chT    = CHAN ATTRIBUT chT    'NATURE' 'DISCRET' ;
  270. chPORO = CHAN ATTRIBUT chPORO 'NATURE' 'DISCRET' ;
  271. chYOGC = CHAN ATTRIBUT chYOGC 'NATURE' 'DISCRET' ;
  272. chALPC = CHAN ATTRIBUT chALPC 'NATURE' 'DISCRET' ;
  273. chTFIS = CHAN ATTRIBUT chTFIS 'NATURE' 'DISCRET' ;
  274. chTUO2 = CHAN ATTRIBUT chTUO2 'NATURE' 'DISCRET' ;
  275. chFACF = CHAN ATTRIBUT chFACF 'NATURE' 'DISCRET' ;
  276. chDSIU = CHAN ATTRIBUT chDSIU 'NATURE' 'DISCRET' ;
  277. chDGRA = CHAN ATTRIBUT chDGRA 'NATURE' 'DISCRET' ;
  278.  
  279. Heavysid = EVOL MANU 'temps'  (PROG 0.0 1000.0)
  280.                      'signal' (PROG 1.0 1.0) ;
  281.  
  282. charT    = CHAR 'T   ' chT    Heavysid ;
  283. charPORO = CHAR 'PORO' chPORO Heavysid ;
  284. charYOGC = CHAR 'YOGC' chYOGC Heavysid ;
  285. charALPC = CHAR 'ALPC' chALPC Heavysid ;
  286. charTFIS = CHAR 'TFIS' chTFIS Heavysid ;
  287. charTUO2 = CHAR 'TUO2' chTUO2 Heavysid ;
  288. charFACF = CHAR 'FACF' chFACF Heavysid ;
  289. charDSIU = CHAR 'DSIU' chDSIU Heavysid ;
  290. charDGRA = CHAR 'DGRA' chDGRA Heavysid ;
  291.  
  292. charuti2 = cha0     ET charT    ET charPORO ET charYOGC ET charALPC ET
  293.            charTFIS ET charTUO2 ET charFACF ET charDSIU ET charDGRA ;
  294.  
  295. * 2.4 Definition des instants de calcul : IDEM 1.4
  296.  
  297. * 2.5 Initialisation des variables internes
  298.  
  299. chpxx02  = MANU CHML mo_uti2 'PXX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  300. chpyy02  = MANU CHML mo_uti2 'PYY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  301. chpzz02  = MANU CHML mo_uti2 'PZZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  302. chpxy02  = MANU CHML mo_uti2 'PXY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  303. chpyz02  = MANU CHML mo_uti2 'PYZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  304. chpzx02  = MANU CHML mo_uti2 'PZX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  305.  
  306. chvnxx02 = MANU CHML mo_uti2 'VNXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  307. chvnyy02 = MANU CHML mo_uti2 'VNYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  308. chvnzz02 = MANU CHML mo_uti2 'VNZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  309. chvnxy02 = MANU CHML mo_uti2 'VNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  310. chvnyz02 = MANU CHML mo_uti2 'VNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  311. chvnzx02 = MANU CHML mo_uti2 'VNZX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  312.  
  313. chepce02 = MANU CHML mo_uti2 'EPCE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  314. chepsh02 = MANU CHML mo_uti2 'EPSH' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  315. chll02   = MANU CHML mo_uti2 'LL  ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  316. chepse02 = MANU CHML mo_uti2 'EPSE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  317.  
  318. chenxx02 = MANU CHML mo_uti2 'ENXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  319. chenyy02 = MANU CHML mo_uti2 'ENYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  320. chenzz02 = MANU CHML mo_uti2 'ENZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  321. chgnxy02 = MANU CHML mo_uti2 'GNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  322. chgnxz02 = MANU CHML mo_uti2 'GNXZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  323. chgnyz02 = MANU CHML mo_uti2 'GNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  324.  
  325. var_ut02 = chpxx02  ET chpyy02  ET chpzz02  ET
  326.            chpxy02  ET chpyz02  ET chpzx02  ET
  327.            chvnxx02 ET chvnyy02 ET chvnzz02 ET
  328.            chvnxy02 ET chvnyz02 ET chvnzx02 ET
  329.            chepce02 ET chepsh02 ET chll02   ET chepse02 ET
  330.            chenxx02 ET chenyy02 ET chenzz02 ET
  331.            chgnxy02 ET chgnxz02 ET chgnyz02 ;
  332.  
  333. * 2.6 Resolution par PASAPAS
  334.  
  335. tabuti2 = TABLE ;
  336. tabuti2.'MODELE' = mo_uti2 ;
  337. tabuti2.'CARACTERISTIQUES' = ma_uti2 ;
  338. tabuti2.'CHARGEMENT' = charuti2 ;
  339. tabuti2.'BLOCAGES_MECANIQUES' = cl0 ;
  340. tabuti2.'VARIABLES_INTERNES' = TABLE ;
  341. tabuti2.'VARIABLES_INTERNES' . 0 = var_ut02 ;
  342. tabuti2.'TEMPS_CALCULES' = pro0 ;
  343. tabuti2.'TEMPS_SAUVES' = pro1 ;
  344. TMASAU=table;
  345. tabuti2 . 'MES_SAUVEGARDES'=TMASAU;
  346. tabuti2.'PRECISION' = 1E-10;
  347. TMASAU .'DEFTO'=VRAI;
  348. TMASAU .'DEFIN'=VRAI;
  349. TEMPS 'ZERO' ;
  350. PASAPAS tabuti2 ;
  351. TABTPS = TEMP 'NOEC';
  352. CPUext2 = TABTPS.'TEMPS_CPU'.'INITIAL';
  353.  
  354. *-----------------------------------------------------------------------
  355. * 3 - Post-traitement : comparaisons entre les deux modeles
  356. *
  357. *     N.B. Comme on simule un essai de cisaillement XY :
  358. *          - Critere strict (errrel < 1.d-8) sur GAXY, GNXY et SMXY
  359. *          - Pour les autres composantes on note les depassements du
  360. *            critere si impression demandee
  361.  
  362. uxu1 = TABLE ; uyu1 = TABLE ; uzu1 = TABLE ;
  363. epxxu1 = TABLE ; epyyu1 = TABLE ; epzzu1 = TABLE ;
  364. gaxyu1 = TABLE ; gaxzu1 = TABLE ; gayzu1 = TABLE ;
  365. enxxu1 = TABLE ; enyyu1 = TABLE ; enzzu1 = TABLE ;
  366. gnxyu1 = TABLE ; gnxzu1 = TABLE ; gnyzu1 = TABLE ;
  367. smxxu1 = TABLE ; smyyu1 = TABLE ; smzzu1 = TABLE ;
  368. smxyu1 = TABLE ; smxzu1 = TABLE ; smyzu1 = TABLE ;
  369. pxxu1  = TABLE ; pyyu1  = TABLE ; pzzu1  = TABLE ;
  370. pxyu1  = TABLE ; pyzu1  = TABLE ; pzxu1  = TABLE ;
  371. vnxxu1 = TABLE ; vnyyu1 = TABLE ; vnzzu1 = TABLE ;
  372. vnxyu1 = TABLE ; vnyzu1 = TABLE ; vnzxu1 = TABLE ;
  373. epceu1 = TABLE ; epshu1 = TABLE ; llu1   = TABLE ; epseu1 = TABLE ;
  374.  
  375. uxu2 = TABLE ; uyu2 = TABLE ; uzu2 = TABLE ;
  376. epxxu2 = TABLE ; epyyu2 = TABLE ; epzzu2 = TABLE ;
  377. gaxyu2 = TABLE ; gaxzu2 = TABLE ; gayzu2 = TABLE ;
  378. enxxu2 = TABLE ; enyyu2 = TABLE ; enzzu2 = TABLE ;
  379. gnxyu2 = TABLE ; gnxzu2 = TABLE ; gnyzu2 = TABLE ;
  380. smxxu2 = TABLE ; smyyu2 = TABLE ; smzzu2 = TABLE ;
  381. smxyu2 = TABLE ; smxzu2 = TABLE ; smyzu2 = TABLE ;
  382. pxxu2  = TABLE ; pyyu2  = TABLE ; pzzu2  = TABLE ;
  383. pxyu2  = TABLE ; pyzu2  = TABLE ; pzxu2  = TABLE ;
  384. vnxxu2 = TABLE ; vnyyu2 = TABLE ; vnzzu2 = TABLE ;
  385. vnxyu2 = TABLE ; vnyzu2 = TABLE ; vnzxu2 = TABLE ;
  386. epceu2 = TABLE ; epshu2 = TABLE ; llu2   = TABLE ; epseu2 = TABLE ;
  387.  
  388. LDEPL = MOTS 'UX  ' 'UY  ' 'UZ  ' ;
  389. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  390. LDENN = MOTS 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  391. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  392. LVARI = MOTS 'PXX ' 'PYY ' 'PZZ ' 'PXY ' 'PYZ ' 'PZX '
  393.              'VNXX' 'VNYY' 'VNZZ' 'VNXY' 'VNYZ' 'VNZX'
  394.              'EPCE' 'EPSH' 'LL  ' 'EPSE' ;
  395.  
  396. idefo = 4 ;
  397. idenn = 4 ;
  398. icont = 4 ;
  399.  
  400. ierr = 0 ;
  401. critere = 1.5d-8 ;
  402. IMESS = 0 ;
  403. ***IMESS = 1 ;
  404.  
  405. nbno  = 8 ;
  406. nbgau = 8 ;
  407.  
  408. netat = DIME (tabuti1.'TEMPS') ;
  409. ietat = -1 ;
  410.  
  411.  
  412. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  413.  
  414. REPETER BECART netat ;
  415.  
  416.    ietat = ietat + 1 ;
  417.  
  418.    temps1 = tabuti1.'TEMPS'.ietat ;
  419.    temps2 = tabuti2.'TEMPS'.ietat ;
  420.    SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  421.       ierr = 1 ;
  422.       QUITTER BECART ;
  423.    FINSI ;
  424.  
  425. *  Extraction des resultats de calcul : Modele utilisateur no21---------
  426.  
  427. *  Deplacements
  428.    uxu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  429.    uxu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  430.    uxu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  431.    uxu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  432.    uxu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  433.    uxu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  434.    uxu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  435.    uxu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  436.  
  437.    uyu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  438.    uyu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  439.    uyu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  440.    uyu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  441.    uyu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  442.    uyu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  443.    uyu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  444.    uyu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  445.  
  446.    uzu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  447.    uzu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  448.    uzu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  449.    uzu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  450.    uzu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  451.    uzu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  452.    uzu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  453.    uzu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  454.  
  455. *  Deformations totales
  456.    igau = 0 ;
  457.    REPETER BEPXXU1 nbgau ;
  458.       igau = igau + 1 ;
  459.       epxxu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau;
  460.    FIN BEPXXU1 ;
  461.  
  462.    igau = 0 ;
  463.    REPETER BEPYYU1 nbgau ;
  464.       igau = igau + 1 ;
  465.       epyyu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau;
  466.    FIN BEPYYU1 ;
  467.  
  468.    igau = 0 ;
  469.    REPETER BEPZZU1 nbgau ;
  470.       igau = igau + 1 ;
  471.       epzzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau;
  472.    FIN BEPZZU1 ;
  473.  
  474.    igau = 0 ;
  475.    REPETER BGAXYU1 nbgau ;
  476.       igau = igau + 1 ;
  477.       gaxyu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau;
  478.    FIN BGAXYU1 ;
  479.  
  480.    igau = 0 ;
  481.    REPETER BGAXZU1 nbgau ;
  482.       igau = igau + 1 ;
  483.       gaxzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau;
  484.    FIN BGAXZU1 ;
  485.  
  486.    igau = 0 ;
  487.    REPETER BGAYZU1 nbgau ;
  488.       igau = igau + 1 ;
  489.       gayzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau;
  490.    FIN BGAYZU1 ;
  491.  
  492. *  Deformations inelastiques (variables internes supplementaires)
  493.    igau = 0 ;
  494.    REPETER BENXXU1 nbgau ;
  495.       igau = igau + 1 ;
  496.       enxxu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  497.                               'ENXX' 1 1 igau ;
  498.    FIN BENXXU1 ;
  499.  
  500.    igau = 0 ;
  501.    REPETER BENYYU1 nbgau ;
  502.       igau = igau + 1 ;
  503.       enyyu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  504.                               'ENYY' 1 1 igau ;
  505.    FIN BENYYU1 ;
  506.  
  507.    igau = 0 ;
  508.    REPETER BENZZU1 nbgau ;
  509.       igau = igau + 1 ;
  510.       enzzu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  511.                               'ENZZ' 1 1 igau ;
  512.    FIN BENZZU1 ;
  513.  
  514.    igau = 0 ;
  515.    REPETER BGNXYU1 nbgau ;
  516.       igau = igau + 1 ;
  517.       gnxyu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  518.                               'GNXY' 1 1 igau ;
  519.    FIN BGNXYU1 ;
  520.  
  521.    igau = 0 ;
  522.    REPETER BGNXZU1 nbgau ;
  523.       igau = igau + 1 ;
  524.       gnxzu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  525.                               'GNXZ' 1 1 igau ;
  526.    FIN BGNXZU1 ;
  527.  
  528.    igau = 0 ;
  529.    REPETER BGNYZU1 nbgau ;
  530.       igau = igau + 1 ;
  531.       gnyzu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  532.                               'GNYZ' 1 1 igau ;
  533.    FIN BGNYZU1 ;
  534.  
  535. *  Contraintes
  536.    igau = 0 ;
  537.    REPETER BSMXXU1 nbgau ;
  538.       igau = igau + 1 ;
  539.       smxxu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  540.    FIN BSMXXU1 ;
  541.  
  542.    igau = 0 ;
  543.    REPETER BSMYYU1 nbgau ;
  544.       igau = igau + 1 ;
  545.       smyyu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  546.    FIN BSMYYU1 ;
  547.  
  548.    igau = 0 ;
  549.    REPETER BSMZZU1 nbgau ;
  550.       igau = igau + 1 ;
  551.       smzzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  552.    FIN BSMZZU1 ;
  553.  
  554.    igau = 0 ;
  555.    REPETER BSMXYU1 nbgau ;
  556.       igau = igau + 1 ;
  557.       smxyu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  558.    FIN BSMXYU1 ;
  559.  
  560.    igau = 0 ;
  561.    REPETER BSMXZU1 nbgau ;
  562.       igau = igau + 1 ;
  563.       smxzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  564.    FIN BSMXZU1 ;
  565.  
  566.    igau = 0 ;
  567.    REPETER BSMYZU1 nbgau ;
  568.       igau = igau + 1 ;
  569.       smyzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  570.    FIN BSMYZU1 ;
  571.  
  572. *  Variables internes
  573.    igau = 0 ;
  574.    REPETER BPXXU1 nbgau ;
  575.       igau = igau + 1 ;
  576.       pxxu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  577.                              'PXX ' 1 1 igau ;
  578.    FIN BPXXU1 ;
  579.  
  580.    igau = 0 ;
  581.    REPETER BPYYU1 nbgau ;
  582.       igau = igau + 1 ;
  583.       pyyu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  584.                              'PYY ' 1 1 igau ;
  585.    FIN BPYYU1 ;
  586.  
  587.    igau = 0 ;
  588.    REPETER BPZZU1 nbgau ;
  589.       igau = igau + 1 ;
  590.       pzzu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  591.                              'PZZ ' 1 1 igau ;
  592.    FIN BPZZU1 ;
  593.  
  594.    igau = 0 ;
  595.    REPETER BPXYU1 nbgau ;
  596.       igau = igau + 1 ;
  597.       pxyu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  598.                              'PXY ' 1 1 igau ;
  599.    FIN BPXYU1 ;
  600.  
  601.    igau = 0 ;
  602.    REPETER BPYZU1 nbgau ;
  603.       igau = igau + 1 ;
  604.       pyzu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  605.                              'PYZ ' 1 1 igau ;
  606.    FIN BPYZU1 ;
  607.  
  608.    igau = 0 ;
  609.    REPETER BPZXU1 nbgau ;
  610.       igau = igau + 1 ;
  611.       pzxu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  612.                              'PZX ' 1 1 igau ;
  613.    FIN BPZXU1 ;
  614.  
  615.    igau = 0 ;
  616.    REPETER BVNXXU1 nbgau ;
  617.       igau = igau + 1 ;
  618.       vnxxu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  619.                              'VNXX' 1 1 igau ;
  620.    FIN BVNXXU1 ;
  621.  
  622.    igau = 0 ;
  623.    REPETER BVNYYU1 nbgau ;
  624.       igau = igau + 1 ;
  625.       vnyyu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  626.                              'VNYY' 1 1 igau ;
  627.    FIN BVNYYU1 ;
  628.  
  629.    igau = 0 ;
  630.    REPETER BVNZZU1 nbgau ;
  631.       igau = igau + 1 ;
  632.       vnzzu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  633.                              'VNZZ' 1 1 igau ;
  634.    FIN BVNZZU1 ;
  635.  
  636.    igau = 0 ;
  637.    REPETER BVNXYU1 nbgau ;
  638.       igau = igau + 1 ;
  639.       vnxyu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  640.                              'VNXY' 1 1 igau ;
  641.    FIN BVNXYU1 ;
  642.  
  643.    igau = 0 ;
  644.    REPETER BVNYZU1 nbgau ;
  645.       igau = igau + 1 ;
  646.       vnyzu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  647.                              'VNYZ' 1 1 igau ;
  648.    FIN BVNYZU1 ;
  649.  
  650.    igau = 0 ;
  651.    REPETER BVNZXU1 nbgau ;
  652.       igau = igau + 1 ;
  653.       vnzxu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  654.                              'VNZX' 1 1 igau ;
  655.    FIN BVNZXU1 ;
  656.  
  657.    igau = 0 ;
  658.    REPETER BEPCEU1 nbgau ;
  659.       igau = igau + 1 ;
  660.       epceu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  661.                              'EPCE' 1 1 igau ;
  662.    FIN BEPCEU1 ;
  663.  
  664.    igau = 0 ;
  665.    REPETER BEPSHU1 nbgau ;
  666.       igau = igau + 1 ;
  667.       epshu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  668.                              'EPSH' 1 1 igau ;
  669.    FIN BEPSHU1 ;
  670.  
  671.    igau = 0 ;
  672.    REPETER BLLU1 nbgau ;
  673.       igau = igau + 1 ;
  674.       llu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  675.                              'LL  ' 1 1 igau ;
  676.    FIN BLLU1 ;
  677.  
  678.    igau = 0 ;
  679.    REPETER BEPSEU1 nbgau ;
  680.       igau = igau + 1 ;
  681.       epseu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  682.                              'EPSE' 1 1 igau ;
  683.    FIN BEPSEU1 ;
  684.  
  685. *  Extraction des resultats de calcul : Modele utilisateur no22---------
  686.  
  687. *  Deplacements
  688.    uxu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  689.    uxu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  690.    uxu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  691.    uxu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  692.    uxu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  693.    uxu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  694.    uxu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  695.    uxu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  696.  
  697.    uyu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  698.    uyu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  699.    uyu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  700.    uyu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  701.    uyu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  702.    uyu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  703.    uyu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  704.    uyu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  705.  
  706.    uzu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  707.    uzu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  708.    uzu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  709.    uzu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  710.    uzu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  711.    uzu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  712.    uzu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  713.    uzu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  714.  
  715. *  Deformations totales
  716.    igau = 0 ;
  717.    REPETER BEPXXU2 nbgau ;
  718.       igau = igau + 1 ;
  719.       epxxu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau;
  720.    FIN BEPXXU2 ;
  721.  
  722.    igau = 0 ;
  723.    REPETER BEPYYU2 nbgau ;
  724.       igau = igau + 1 ;
  725.       epyyu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau;
  726.    FIN BEPYYU2 ;
  727.  
  728.    igau = 0 ;
  729.    REPETER BEPZZU2 nbgau ;
  730.       igau = igau + 1 ;
  731.       epzzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau;
  732.    FIN BEPZZU2 ;
  733.  
  734.    igau = 0 ;
  735.    REPETER BGAXYU2 nbgau ;
  736.       igau = igau + 1 ;
  737.       gaxyu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau;
  738.    FIN BGAXYU2 ;
  739.  
  740.    igau = 0 ;
  741.    REPETER BGAXZU2 nbgau ;
  742.       igau = igau + 1 ;
  743.       gaxzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau;
  744.    FIN BGAXZU2 ;
  745.  
  746.    igau = 0 ;
  747.    REPETER BGAYZU2 nbgau ;
  748.       igau = igau + 1 ;
  749.       gayzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau;
  750.    FIN BGAYZU2 ;
  751.  
  752. *  Deformations inelastiques (variables internes supplementaires)
  753.    igau = 0 ;
  754.    REPETER BENXXU2 nbgau ;
  755.       igau = igau + 1 ;
  756.       enxxu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  757.                               'ENXX' 1 1 igau ;
  758.    FIN BENXXU2 ;
  759.  
  760.    igau = 0 ;
  761.    REPETER BENYYU2 nbgau ;
  762.       igau = igau + 1 ;
  763.       enyyu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  764.                               'ENYY' 1 1 igau ;
  765.    FIN BENYYU2 ;
  766.  
  767.    igau = 0 ;
  768.    REPETER BENZZU2 nbgau ;
  769.       igau = igau + 1 ;
  770.       enzzu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  771.                               'ENZZ' 1 1 igau ;
  772.    FIN BENZZU2 ;
  773.  
  774.    igau = 0 ;
  775.    REPETER BGNXYU2 nbgau ;
  776.       igau = igau + 1 ;
  777.       gnxyu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  778.                               'GNXY' 1 1 igau ;
  779.    FIN BGNXYU2 ;
  780.  
  781.    igau = 0 ;
  782.    REPETER BGNXZU2 nbgau ;
  783.       igau = igau + 1 ;
  784.       gnxzu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  785.                               'GNXZ' 1 1 igau ;
  786.    FIN BGNXZU2 ;
  787.  
  788.    igau = 0 ;
  789.    REPETER BGNYZU2 nbgau ;
  790.       igau = igau + 1 ;
  791.       gnyzu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  792.                               'GNYZ' 1 1 igau ;
  793.    FIN BGNYZU2 ;
  794.  
  795. *  Contraintes
  796.    igau = 0 ;
  797.    REPETER BSMXXU2 nbgau ;
  798.       igau = igau + 1 ;
  799.       smxxu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  800.    FIN BSMXXU2 ;
  801.  
  802.    igau = 0 ;
  803.    REPETER BSMYYU2 nbgau ;
  804.       igau = igau + 1 ;
  805.       smyyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  806.    FIN BSMYYU2 ;
  807.  
  808.    igau = 0 ;
  809.    REPETER BSMZZU2 nbgau ;
  810.       igau = igau + 1 ;
  811.       smzzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  812.    FIN BSMZZU2 ;
  813.  
  814.    igau = 0 ;
  815.    REPETER BSMXYU2 nbgau ;
  816.       igau = igau + 1 ;
  817.       smxyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  818.    FIN BSMXYU2 ;
  819.  
  820.    igau = 0 ;
  821.    REPETER BSMXZU2 nbgau ;
  822.       igau = igau + 1 ;
  823.       smxzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  824.    FIN BSMXZU2 ;
  825.  
  826.    igau = 0 ;
  827.    REPETER BSMYZU2 nbgau ;
  828.       igau = igau + 1 ;
  829.       smyzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  830.    FIN BSMYZU2 ;
  831.  
  832. *  Variables internes
  833.    igau = 0 ;
  834.    REPETER BPXXU2 nbgau ;
  835.       igau = igau + 1 ;
  836.       pxxu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  837.                              'PXX ' 1 1 igau ;
  838.    FIN BPXXU2 ;
  839.  
  840.    igau = 0 ;
  841.    REPETER BPYYU2 nbgau ;
  842.       igau = igau + 1 ;
  843.       pyyu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  844.                              'PYY ' 1 1 igau ;
  845.    FIN BPYYU2 ;
  846.  
  847.    igau = 0 ;
  848.    REPETER BPZZU2 nbgau ;
  849.       igau = igau + 1 ;
  850.       pzzu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  851.                              'PZZ ' 1 1 igau ;
  852.    FIN BPZZU2 ;
  853.  
  854.    igau = 0 ;
  855.    REPETER BPXYU2 nbgau ;
  856.       igau = igau + 1 ;
  857.       pxyu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  858.                              'PXY ' 1 1 igau ;
  859.    FIN BPXYU2 ;
  860.  
  861.    igau = 0 ;
  862.    REPETER BPYZU2 nbgau ;
  863.       igau = igau + 1 ;
  864.       pyzu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  865.                              'PYZ ' 1 1 igau ;
  866.    FIN BPYZU2 ;
  867.  
  868.    igau = 0 ;
  869.    REPETER BPZXU2 nbgau ;
  870.       igau = igau + 1 ;
  871.       pzxu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  872.                              'PZX ' 1 1 igau ;
  873.    FIN BPZXU2 ;
  874.  
  875.    igau = 0 ;
  876.    REPETER BVNXXU2 nbgau ;
  877.       igau = igau + 1 ;
  878.       vnxxu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  879.                              'VNXX' 1 1 igau ;
  880.    FIN BVNXXU2 ;
  881.  
  882.    igau = 0 ;
  883.    REPETER BVNYYU2 nbgau ;
  884.       igau = igau + 1 ;
  885.       vnyyu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  886.                              'VNYY' 1 1 igau ;
  887.    FIN BVNYYU2 ;
  888.  
  889.    igau = 0 ;
  890.    REPETER BVNZZU2 nbgau ;
  891.       igau = igau + 1 ;
  892.       vnzzu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  893.                              'VNZZ' 1 1 igau ;
  894.    FIN BVNZZU2 ;
  895.  
  896.    igau = 0 ;
  897.    REPETER BVNXYU2 nbgau ;
  898.       igau = igau + 1 ;
  899.       vnxyu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  900.                              'VNXY' 1 1 igau ;
  901.    FIN BVNXYU2 ;
  902.  
  903.    igau = 0 ;
  904.    REPETER BVNYZU2 nbgau ;
  905.       igau = igau + 1 ;
  906.       vnyzu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  907.                              'VNYZ' 1 1 igau ;
  908.    FIN BVNYZU2 ;
  909.  
  910.    igau = 0 ;
  911.    REPETER BVNZXU2 nbgau ;
  912.       igau = igau + 1 ;
  913.       vnzxu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  914.                              'VNZX' 1 1 igau ;
  915.    FIN BVNZXU2 ;
  916.  
  917.    igau = 0 ;
  918.    REPETER BEPCEU2 nbgau ;
  919.       igau = igau + 1 ;
  920.       epceu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  921.                              'EPCE' 1 1 igau ;
  922.    FIN BEPCEU2 ;
  923.  
  924.    igau = 0 ;
  925.    REPETER BEPSHU2 nbgau ;
  926.       igau = igau + 1 ;
  927.       epshu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  928.                              'EPSH' 1 1 igau ;
  929.    FIN BEPSHU2 ;
  930.  
  931.    igau = 0 ;
  932.    REPETER BLLU2 nbgau ;
  933.       igau = igau + 1 ;
  934.       llu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  935.                              'LL  ' 1 1 igau ;
  936.    FIN BLLU2 ;
  937.  
  938.    igau = 0 ;
  939.    REPETER BEPSEU2 nbgau ;
  940.       igau = igau + 1 ;
  941.       epseu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  942.                              'EPSE' 1 1 igau ;
  943.    FIN BEPSEU2 ;
  944.  
  945. *  Comparaison des deplacements-----------------------------------------
  946.  
  947.    icomp = 1 ;
  948.    ino = 0 ;
  949.    REPETER BUX nbno ;
  950.       ino = ino + 1 ;
  951.       SI (EGA uxu1.ino 0.0) ;
  952.          SI (NON (EGA uxu2.ino 0.0)) ;
  953.             SI (EGA IMESS 1) ;
  954.             Message 'Erreur absolue' uxu2.ino
  955.                     'Composante ' (EXTR LDEPL icomp)
  956.                     ' Noeud' ino 'Etat' ietat ;
  957.             Message ;
  958.             FINSI ;
  959. ************ierr = 10 ; errabs = uxu2.ino ; QUITTER BECART ;
  960.          FINSI ;
  961.       SINON ;
  962.          dux = ABS (uxu1.ino - uxu2.ino) ;
  963.          errrel = ABS (dux / uxu1.ino) ;
  964.          SI (errrel > critere) ;
  965.             SI (EGA IMESS 1) ;
  966.             Message 'Erreur relative' errrel
  967.                     'Composante ' (EXTR LDEPL icomp)
  968.                     ' Noeud' ino 'Etat' ietat ;
  969.             Message 'uxu1.ino =' uxu1.ino 'uxu2.ino =' uxu2.ino ;
  970.             Message ;
  971.             FINSI ;
  972. ************ierr = 11 ; QUITTER BECART ;
  973.          FINSI ;
  974.       FINSI ;
  975.    FIN BUX ;
  976.  
  977.    icomp = 2 ;
  978.    ino = 0 ;
  979.    REPETER BUY nbno ;
  980.       ino = ino + 1 ;
  981.       SI (EGA uyu1.ino 0.0) ;
  982.          SI (NON (EGA uyu2.ino 0.0)) ;
  983.             SI (EGA IMESS 1) ;
  984.             Message 'Erreur absolue' uyu2.ino
  985.                     'Composante ' (EXTR LDEPL icomp)
  986.                     ' Noeud' ino 'Etat' ietat ;
  987.             Message ;
  988.             FINSI ;
  989. ************ierr = 10 ; errabs = uyu2.ino ; QUITTER BECART ;
  990.          FINSI ;
  991.       SINON ;
  992.          duy = ABS (uyu1.ino - uyu2.ino) ;
  993.          errrel = ABS (duy / uyu1.ino) ;
  994.          SI (errrel > critere) ;
  995.             SI (EGA IMESS 1) ;
  996.             Message 'Erreur relative' errrel
  997.                     'Composante ' (EXTR LDEPL icomp)
  998.                     ' Noeud' ino 'Etat' ietat ;
  999.             Message 'uyu1.ino =' uyu1.ino 'uyu2.ino =' uyu2.ino ;
  1000.             Message ;
  1001.             FINSI ;
  1002. ************ierr = 11 ; QUITTER BECART ;
  1003.          FINSI ;
  1004.       FINSI ;
  1005.    FIN BUY ;
  1006.  
  1007.    icomp = 3 ;
  1008.    ino = 0 ;
  1009.    REPETER BUZ nbno ;
  1010.       ino = ino + 1 ;
  1011.       SI (EGA uzu1.ino 0.0) ;
  1012.          SI (NON (EGA uzu2.ino 0.0)) ;
  1013.             SI (EGA IMESS 1) ;
  1014.             Message 'Erreur absolue' uzu2.ino
  1015.                     'Composante ' (EXTR LDEPL icomp)
  1016.                     ' Noeud' ino 'Etat' ietat ;
  1017.             Message ;
  1018.             FINSI ;
  1019. ************ierr = 10 ; errabs = uzu2.ino ; QUITTER BECART ;
  1020.          FINSI ;
  1021.       SINON ;
  1022.          duz = ABS (uzu1.ino - uzu2.ino) ;
  1023.          errrel = ABS (duz / uzu1.ino) ;
  1024.          SI (errrel > critere) ;
  1025.             SI (EGA IMESS 1) ;
  1026.             Message 'Erreur relative' errrel
  1027.                     'Composante ' (EXTR LDEPL icomp)
  1028.                     ' Noeud' ino 'Etat' ietat ;
  1029.             Message 'uzu1.ino =' uzu1.ino 'uzu2.ino =' uzu2.ino ;
  1030.             Message ;
  1031.             FINSI ;
  1032. ************ierr = 11 ; QUITTER BECART ;
  1033.          FINSI ;
  1034.       FINSI ;
  1035.    FIN BUZ ;
  1036.  
  1037. *  Comparaison des deformations-----------------------------------------
  1038.  
  1039.    icomp = 1 ;
  1040.    igau = 0 ;
  1041.    REPETER BEPXX nbgau ;
  1042.       igau = igau + 1 ;
  1043.       SI (EGA epxxu1.igau 0.0) ;
  1044.          SI (NON (EGA epxxu2.igau 0.0)) ;
  1045.             SI (EGA IMESS 1) ;
  1046.             Message 'Erreur absolue' epxxu2.igau
  1047.                     'Composante ' (EXTR LDEFO icomp)
  1048.                     ' Point Gauss' igau 'Etat' ietat ;
  1049.             Message ;
  1050.             FINSI ;
  1051. ************ierr = 20 ; errabs = epxxu2.igau ; QUITTER BECART ;
  1052.          FINSI ;
  1053.       SINON ;
  1054.          depxx = ABS (epxxu1.igau - epxxu2.igau) ;
  1055.          errrel = ABS (depxx / epxxu1.igau) ;
  1056.          SI (errrel > critere) ;
  1057.             SI (EGA IMESS 1) ;
  1058.             Message 'Erreur relative' errrel
  1059.                     'Composante ' (EXTR LDEFO icomp)
  1060.                     ' Point Gauss' igau 'Etat' ietat ;
  1061.             Message 'epxxu1.igau =' epxxu1.igau
  1062.                     'epxxu2.igau =' epxxu2.igau;
  1063.             Message ;
  1064.             FINSI ;
  1065. ************ierr = 21 ; QUITTER BECART ;
  1066.          FINSI ;
  1067.       FINSI ;
  1068.    FIN BEPXX ;
  1069.  
  1070.    icomp = 2 ;
  1071.    igau = 0 ;
  1072.    REPETER BEPYY nbgau ;
  1073.       igau = igau + 1 ;
  1074.       SI (EGA epyyu1.igau 0.0) ;
  1075.          SI (NON (EGA epyyu2.igau 0.0)) ;
  1076.             SI (EGA IMESS 1) ;
  1077.             Message 'Erreur absolue' epyyu2.igau
  1078.                     'Composante ' (EXTR LDEFO icomp)
  1079.                     ' Point Gauss' igau 'Etat' ietat ;
  1080.             Message ;
  1081.             FINSI ;
  1082. ************ierr = 20 ; errabs = epyyu2.igau ; QUITTER BECART ;
  1083.          FINSI ;
  1084.       SINON ;
  1085.          depyy = ABS (epyyu1.igau - epyyu2.igau) ;
  1086.          errrel = ABS (depyy / epyyu1.igau) ;
  1087.          SI (errrel > critere) ;
  1088.             SI (EGA IMESS 1) ;
  1089.             Message 'Erreur relative' errrel
  1090.                     'Composante ' (EXTR LDEFO icomp)
  1091.                     ' Point Gauss' igau 'Etat' ietat ;
  1092.             Message 'epyyu1.igau =' epyyu1.igau
  1093.                     'epyyu2.igau =' epyyu2.igau;
  1094.             Message ;
  1095.             FINSI ;
  1096. ************ierr = 21 ; QUITTER BECART ;
  1097.          FINSI ;
  1098.       FINSI ;
  1099.    FIN BEPYY ;
  1100.  
  1101.    icomp = 3 ;
  1102.    igau = 0 ;
  1103.    REPETER BEPZZ nbgau ;
  1104.       igau = igau + 1 ;
  1105.       SI (EGA epzzu1.igau 0.0) ;
  1106.          SI (NON (EGA epzzu2.igau 0.0)) ;
  1107.             SI (EGA IMESS 1) ;
  1108.             Message 'Erreur absolue' epzzu2.igau
  1109.                     'Composante ' (EXTR LDEFO icomp)
  1110.                     ' Point Gauss' igau 'Etat' ietat ;
  1111.             Message ;
  1112.             FINSI ;
  1113. ************ierr = 20 ; errabs = epzzu2.igau ; QUITTER BECART ;
  1114.          FINSI ;
  1115.       SINON ;
  1116.          depzz = ABS (epzzu1.igau - epzzu2.igau) ;
  1117.          errrel = ABS (depzz / epzzu1.igau) ;
  1118.          SI (errrel > critere) ;
  1119.             SI (EGA IMESS 1) ;
  1120.             Message 'Erreur relative' errrel
  1121.                     'Composante ' (EXTR LDEFO icomp)
  1122.                     ' Point Gauss' igau 'Etat' ietat ;
  1123.             Message 'epzzu1.igau =' epzzu1.igau
  1124.                     'epzzu2.igau =' epzzu2.igau;
  1125.             Message ;
  1126.             FINSI ;
  1127. ************ierr = 21 ; QUITTER BECART ;
  1128.          FINSI ;
  1129.       FINSI ;
  1130.    FIN BEPZZ ;
  1131.  
  1132.    icomp = 4 ;
  1133.    igau = 0 ;
  1134.    REPETER BGAXY nbgau ;
  1135.       igau = igau + 1 ;
  1136.       SI (EGA gaxyu1.igau 0.0) ;
  1137.          SI (NON (EGA gaxyu2.igau 0.0)) ;
  1138. ************SI (EGA IMESS 1) ;
  1139. ************Message 'Erreur absolue' gaxyu2.igau
  1140. ************        'Composante ' (EXTR LDEFO icomp)
  1141. ************        ' Point Gauss' igau 'Etat' ietat ;
  1142. ************Message ;
  1143. ************FINSI ;
  1144.             ierr = 20 ; errabs = gaxyu2.igau ; QUITTER BECART ;
  1145.          FINSI ;
  1146.       SINON ;
  1147.          dgaxy = ABS (gaxyu1.igau - gaxyu2.igau) ;
  1148.          errrel = ABS (dgaxy / gaxyu1.igau) ;
  1149.          SI (errrel > critere) ;
  1150. ************SI (EGA IMESS 1) ;
  1151. ************Message 'Erreur relative' errrel
  1152. ************        'Composante ' (EXTR LDEFO icomp)
  1153. ************        ' Point Gauss' igau 'Etat' ietat ;
  1154. ************Message 'gaxyu1.igau =' gaxyu1.igau
  1155. ************        'gaxyu2.igau =' gaxyu2.igau;
  1156. ************Message ;
  1157. ************FINSI ;
  1158.             ierr = 21 ; QUITTER BECART ;
  1159.          FINSI ;
  1160.       FINSI ;
  1161.    FIN BGAXY ;
  1162.  
  1163.    icomp = 5 ;
  1164.    igau = 0 ;
  1165.    REPETER BGAXZ nbgau ;
  1166.       igau = igau + 1 ;
  1167.       SI (EGA gaxzu1.igau 0.0) ;
  1168.          SI (NON (EGA gaxzu2.igau 0.0)) ;
  1169.             SI (EGA IMESS 1) ;
  1170.             Message 'Erreur absolue' gaxzu2.igau
  1171.                     'Composante ' (EXTR LDEFO icomp)
  1172.                     ' Point Gauss' igau 'Etat' ietat ;
  1173.             Message ;
  1174.             FINSI ;
  1175. ************ierr = 20 ; errabs = gaxzu2.igau ; QUITTER BECART ;
  1176.          FINSI ;
  1177.       SINON ;
  1178.          dgaxz = ABS (gaxzu1.igau - gaxzu2.igau) ;
  1179.          errrel = ABS (dgaxz / gaxzu1.igau) ;
  1180.          SI (errrel > critere) ;
  1181.             SI (EGA IMESS 1) ;
  1182.             Message 'Erreur relative' errrel
  1183.                     'Composante ' (EXTR LDEFO icomp)
  1184.                     ' Point Gauss' igau 'Etat' ietat ;
  1185.             Message 'gaxzu1.igau =' gaxzu1.igau
  1186.                     'gaxzu2.igau =' gaxzu2.igau;
  1187.             Message ;
  1188.             FINSI ;
  1189. ************ierr = 21 ; QUITTER BECART ;
  1190.          FINSI ;
  1191.       FINSI ;
  1192.    FIN BGAXZ ;
  1193.  
  1194.    icomp = 6 ;
  1195.    igau = 0 ;
  1196.    REPETER BGAYZ nbgau ;
  1197.       igau = igau + 1 ;
  1198.       SI (EGA gayzu1.igau 0.0) ;
  1199.          SI (NON (EGA gayzu2.igau 0.0)) ;
  1200.             SI (EGA IMESS 1) ;
  1201.             Message 'Erreur absolue' gayzu2.igau
  1202.                     'Composante ' (EXTR LDEFO icomp)
  1203.                     ' Point Gauss' igau 'Etat' ietat ;
  1204.             Message ;
  1205.             FINSI ;
  1206. ************ierr = 20 ; errabs = gayzu2.igau ; QUITTER BECART ;
  1207.          FINSI ;
  1208.       SINON ;
  1209.          dgayz = ABS (gayzu1.igau - gayzu2.igau) ;
  1210.          errrel = ABS (dgayz / gayzu1.igau) ;
  1211.          SI (errrel > critere) ;
  1212.             SI (EGA IMESS 1) ;
  1213.             Message 'Erreur relative' errrel
  1214.                     'Composante ' (EXTR LDEFO icomp)
  1215.                     ' Point Gauss' igau 'Etat' ietat ;
  1216.             Message 'gayzu1.igau =' gayzu1.igau
  1217.                     'gayzu2.igau =' gayzu2.igau;
  1218.             Message ;
  1219.             FINSI ;
  1220. ************ierr = 21 ; QUITTER BECART ;
  1221.          FINSI ;
  1222.       FINSI ;
  1223.    FIN BGAYZ ;
  1224.  
  1225. *  Comparaison des deformations inelastiques----------------------------
  1226.  
  1227.    icomp = 1 ;
  1228.    igau = 0 ;
  1229.    REPETER BENXX nbgau ;
  1230.       igau = igau + 1 ;
  1231.       SI (EGA enxxu1.igau 0.0) ;
  1232.          SI (NON (EGA enxxu2.igau 0.0)) ;
  1233.             SI (EGA IMESS 1) ;
  1234.             Message 'Erreur absolue' enxxu2.igau
  1235.                     'Composante ' (EXTR LDENN icomp)
  1236.                     ' Point Gauss' igau 'Etat' ietat ;
  1237.             Message ;
  1238.             FINSI ;
  1239. ************ierr = 30 ; errabs = enxxu2.igau ; QUITTER BECART ;
  1240.          FINSI ;
  1241.       SINON ;
  1242.          deixx = ABS (enxxu1.igau - enxxu2.igau) ;
  1243.          errrel = ABS (deixx / enxxu1.igau) ;
  1244.          SI (errrel > critere) ;
  1245.             SI (EGA IMESS 1) ;
  1246.             Message 'Erreur relative' errrel
  1247.                     'Composante ' (EXTR LDENN icomp)
  1248.                     ' Point Gauss' igau 'Etat' ietat ;
  1249.             Message 'enxxu1.igau =' enxxu1.igau
  1250.                     'enxxu2.igau =' enxxu2.igau;
  1251.             Message ;
  1252.             FINSI ;
  1253. ************ierr = 31 ; QUITTER BECART ;
  1254.          FINSI ;
  1255.       FINSI ;
  1256.    FIN BENXX ;
  1257.  
  1258.    icomp = 2 ;
  1259.    igau = 0 ;
  1260.    REPETER BENYY nbgau ;
  1261.       igau = igau + 1 ;
  1262.       SI (EGA enyyu1.igau 0.0) ;
  1263.          SI (NON (EGA enyyu2.igau 0.0)) ;
  1264.             SI (EGA IMESS 1) ;
  1265.             Message 'Erreur absolue' enyyu2.igau
  1266.                     'Composante ' (EXTR LDENN icomp)
  1267.                     ' Point Gauss' igau 'Etat' ietat ;
  1268.             Message ;
  1269.             FINSI ;
  1270. ************ierr = 30 ; errabs = enyyu2.igau ; QUITTER BECART ;
  1271.          FINSI ;
  1272.       SINON ;
  1273.          deiyy = ABS (enyyu1.igau - enyyu2.igau) ;
  1274.          errrel = ABS (deiyy / enyyu1.igau) ;
  1275.          SI (errrel > critere) ;
  1276.             SI (EGA IMESS 1) ;
  1277.             Message 'Erreur relative' errrel
  1278.                     'Composante ' (EXTR LDENN icomp)
  1279.                     ' Point Gauss' igau 'Etat' ietat ;
  1280.             Message 'enyyu1.igau =' enyyu1.igau
  1281.                     'enyyu2.igau =' enyyu2.igau;
  1282.             Message ;
  1283.             FINSI ;
  1284. ************ierr = 31 ; QUITTER BECART ;
  1285.          FINSI ;
  1286.       FINSI ;
  1287.    FIN BENYY ;
  1288.  
  1289.    icomp = 3 ;
  1290.    igau = 0 ;
  1291.    REPETER BENZZ nbgau ;
  1292.       igau = igau + 1 ;
  1293.       SI (EGA enzzu1.igau 0.0) ;
  1294.          SI (NON (EGA enzzu2.igau 0.0)) ;
  1295.             SI (EGA IMESS 1) ;
  1296.             Message 'Erreur absolue' enzzu2.igau
  1297.                     'Composante ' (EXTR LDENN icomp)
  1298.                     ' Point Gauss' igau 'Etat' ietat ;
  1299.             Message ;
  1300.             FINSI ;
  1301. ************ierr = 30 ; errabs = enzzu2.igau ; QUITTER BECART ;
  1302.          FINSI ;
  1303.       SINON ;
  1304.          deizz = ABS (enzzu1.igau - enzzu2.igau) ;
  1305.          errrel = ABS (deizz / enzzu1.igau) ;
  1306.          SI (errrel > critere) ;
  1307.             SI (EGA IMESS 1) ;
  1308.             Message 'Erreur relative' errrel
  1309.                     'Composante ' (EXTR LDENN icomp)
  1310.                     ' Point Gauss' igau 'Etat' ietat ;
  1311.             Message 'enzzu1.igau =' enzzu1.igau
  1312.                     'enzzu2.igau =' enzzu2.igau;
  1313.             Message ;
  1314.             FINSI ;
  1315. ************ierr = 31 ; QUITTER BECART ;
  1316.          FINSI ;
  1317.       FINSI ;
  1318.    FIN BENZZ ;
  1319.  
  1320.    icomp = 4 ;
  1321.    igau = 0 ;
  1322.    REPETER BGNXY nbgau ;
  1323.       igau = igau + 1 ;
  1324.       SI (EGA gnxyu1.igau 0.0) ;
  1325.          SI (NON (EGA gnxyu2.igau 0.0)) ;
  1326. ************SI (EGA IMESS 1) ;
  1327. ************Message 'Erreur absolue' gnxyu2.igau
  1328. ************        'Composante ' (EXTR LDENN icomp)
  1329. ************        ' Point Gauss' igau 'Etat' ietat ;
  1330. ************Message ;
  1331. ************FINSI ;
  1332.             ierr = 30 ; errabs = gnxyu2.igau ; QUITTER BECART ;
  1333.          FINSI ;
  1334.       SINON ;
  1335.          dgixy = ABS (gnxyu1.igau - gnxyu2.igau) ;
  1336.          errrel = ABS (dgixy / gnxyu1.igau) ;
  1337.          SI (errrel > critere) ;
  1338. ************SI (EGA IMESS 1) ;
  1339. ************Message 'Erreur relative' errrel
  1340. ************        'Composante ' (EXTR LDENN icomp)
  1341. ************        ' Point Gauss' igau 'Etat' ietat ;
  1342. ************Message 'gnxyu1.igau =' gnxyu1.igau
  1343. ************        'gnxyu2.igau =' gnxyu2.igau;
  1344. ************Message ;
  1345. ************FINSI ;
  1346.             ierr = 31 ; QUITTER BECART ;
  1347.          FINSI ;
  1348.       FINSI ;
  1349.    FIN BGNXY ;
  1350.  
  1351.    icomp = 5 ;
  1352.    igau = 0 ;
  1353.    REPETER BGNXZ nbgau ;
  1354.       igau = igau + 1 ;
  1355.       SI (EGA gnxzu1.igau 0.0) ;
  1356.          SI (NON (EGA gnxzu2.igau 0.0)) ;
  1357.             SI (EGA IMESS 1) ;
  1358.             Message 'Erreur absolue' gnxzu2.igau
  1359.                     'Composante ' (EXTR LDENN icomp)
  1360.                     ' Point Gauss' igau 'Etat' ietat ;
  1361.             Message ;
  1362.             FINSI ;
  1363. ************ierr = 30 ; errabs = gnxzu2.igau ; QUITTER BECART ;
  1364.          FINSI ;
  1365.       SINON ;
  1366.          dgixz = ABS (gnxzu1.igau - gnxzu2.igau) ;
  1367.          errrel = ABS (dgixz / gnxzu1.igau) ;
  1368.          SI (errrel > critere) ;
  1369.             SI (EGA IMESS 1) ;
  1370.             Message 'Erreur relative' errrel
  1371.                     'Composante ' (EXTR LDENN icomp)
  1372.                     ' Point Gauss' igau 'Etat' ietat ;
  1373.             Message 'gnxzu1.igau =' gnxzu1.igau
  1374.                     'gnxzu2.igau =' gnxzu2.igau;
  1375.             Message ;
  1376.             FINSI ;
  1377. ************ierr = 31 ; QUITTER BECART ;
  1378.          FINSI ;
  1379.       FINSI ;
  1380.    FIN BGNXZ ;
  1381.  
  1382.    icomp = 6 ;
  1383.    igau = 0 ;
  1384.    REPETER BGNYZ nbgau ;
  1385.       igau = igau + 1 ;
  1386.       SI (EGA gnyzu1.igau 0.0) ;
  1387.          SI (NON (EGA gnyzu2.igau 0.0)) ;
  1388.             SI (EGA IMESS 1) ;
  1389.             Message 'Erreur absolue' gnyzu2.igau
  1390.                     'Composante ' (EXTR LDENN icomp)
  1391.                     ' Point Gauss' igau 'Etat' ietat ;
  1392.             Message ;
  1393.             FINSI ;
  1394. ************ierr = 30 ; errabs = gnyzu2.igau ; QUITTER BECART ;
  1395.          FINSI ;
  1396.       SINON ;
  1397.          dgiyz = ABS (gnyzu1.igau - gnyzu2.igau) ;
  1398.          errrel = ABS (dgiyz / gnyzu1.igau) ;
  1399.          SI (errrel > critere) ;
  1400.             SI (EGA IMESS 1) ;
  1401.             Message 'Erreur relative' errrel
  1402.                     'Composante ' (EXTR LDENN icomp)
  1403.                     ' Point Gauss' igau 'Etat' ietat ;
  1404.             Message 'gnyzu1.igau =' gnyzu1.igau
  1405.                     'gnyzu2.igau =' gnyzu2.igau;
  1406.             Message ;
  1407.             FINSI ;
  1408. ************ierr = 31 ; QUITTER BECART ;
  1409.          FINSI ;
  1410.       FINSI ;
  1411.    FIN BGNYZ ;
  1412.  
  1413. *  Comparaison des contraintes------------------------------------------
  1414.  
  1415.    icomp = 1 ;
  1416.    igau = 0 ;
  1417.    REPETER BSMXX nbgau ;
  1418.       igau = igau + 1 ;
  1419.       SI (EGA smxxu1.igau 0.0) ;
  1420.          SI (NON (EGA smxxu2.igau 0.0)) ;
  1421.             SI (EGA IMESS 1) ;
  1422.             Message 'Erreur absolue' smxxu2.igau
  1423.                     'Composante ' (EXTR LCONT icomp)
  1424.                     ' Point Gauss' igau 'Etat' ietat ;
  1425.             Message ;
  1426.             FINSI ;
  1427. ************ierr = 40 ; errabs = smxxu2.igau ; QUITTER BECART ;
  1428.          FINSI ;
  1429.       SINON ;
  1430.          dsmxx = ABS (smxxu1.igau - smxxu2.igau) ;
  1431.          errrel = ABS (dsmxx / smxxu1.igau) ;
  1432.          SI (errrel > critere) ;
  1433.             SI (EGA IMESS 1) ;
  1434.             Message 'Erreur relative' errrel
  1435.                     'Composante ' (EXTR LCONT icomp)
  1436.                     ' Point Gauss' igau 'Etat' ietat ;
  1437.             Message 'smxxu1.igau =' smxxu1.igau
  1438.                     'smxxu2.igau =' smxxu2.igau;
  1439.             Message ;
  1440.             FINSI ;
  1441. ************ierr = 41 ; QUITTER BECART ;
  1442.          FINSI ;
  1443.       FINSI ;
  1444.    FIN BSMXX ;
  1445.  
  1446.    icomp = 2 ;
  1447.    igau = 0 ;
  1448.    REPETER BSMYY nbgau ;
  1449.       igau = igau + 1 ;
  1450.       SI (EGA smyyu1.igau 0.0) ;
  1451.          SI (NON (EGA smyyu2.igau 0.0)) ;
  1452.             SI (EGA IMESS 1) ;
  1453.             Message 'Erreur absolue' smyyu2.igau
  1454.                     'Composante ' (EXTR LCONT icomp)
  1455.                     ' Point Gauss' igau 'Etat' ietat ;
  1456.             Message ;
  1457.             FINSI ;
  1458. ************ierr = 40 ; errabs = smyyu2.igau ; QUITTER BECART ;
  1459.          FINSI ;
  1460.       SINON ;
  1461.          dsmyy = ABS (smyyu1.igau - smyyu2.igau) ;
  1462.          errrel = ABS (dsmyy / smyyu1.igau) ;
  1463.          SI (errrel > critere) ;
  1464.             SI (EGA IMESS 1) ;
  1465.             Message 'Erreur relative' errrel
  1466.                     'Composante ' (EXTR LCONT icomp)
  1467.                     ' Point Gauss' igau 'Etat' ietat ;
  1468.             Message 'smyyu1.igau =' smyyu1.igau
  1469.                     'smyyu2.igau =' smyyu2.igau;
  1470.             Message ;
  1471.             FINSI ;
  1472. ************ierr = 41 ; QUITTER BECART ;
  1473.          FINSI ;
  1474.       FINSI ;
  1475.    FIN BSMYY ;
  1476.  
  1477.    icomp = 3 ;
  1478.    igau = 0 ;
  1479.    REPETER BSMZZ nbgau ;
  1480.       igau = igau + 1 ;
  1481.       SI (EGA smzzu1.igau 0.0) ;
  1482.          SI (NON (EGA smzzu2.igau 0.0)) ;
  1483.             SI (EGA IMESS 1) ;
  1484.             Message 'Erreur absolue' smzzu2.igau
  1485.                     'Composante ' (EXTR LCONT icomp)
  1486.                     ' Point Gauss' igau 'Etat' ietat ;
  1487.             Message ;
  1488.             FINSI ;
  1489. ************ierr = 40 ; errabs = smzzu2.igau ; QUITTER BECART ;
  1490.          FINSI ;
  1491.       SINON ;
  1492.          dsmzz = ABS (smzzu1.igau - smzzu2.igau) ;
  1493.          errrel = ABS (dsmzz / smzzu1.igau) ;
  1494.          SI (errrel > critere) ;
  1495.             SI (EGA IMESS 1) ;
  1496.             Message 'Erreur relative' errrel
  1497.                     'Composante ' (EXTR LCONT icomp)
  1498.                     ' Point Gauss' igau 'Etat' ietat ;
  1499.             Message 'smzzu1.igau =' smzzu1.igau
  1500.                     'smzzu2.igau =' smzzu2.igau;
  1501.             Message ;
  1502.             FINSI ;
  1503. ************ierr = 41 ; QUITTER BECART ;
  1504.          FINSI ;
  1505.       FINSI ;
  1506.    FIN BSMZZ ;
  1507.  
  1508.    icomp = 4 ;
  1509.    igau = 0 ;
  1510.    REPETER BSMXY nbgau ;
  1511.       igau = igau + 1 ;
  1512.       SI (EGA smxyu1.igau 0.0) ;
  1513.          SI (NON (EGA smxyu2.igau 0.0)) ;
  1514. ************SI (EGA IMESS 1) ;
  1515. ************Message 'Erreur absolue' smxyu2.igau
  1516. ************        'Composante ' (EXTR LCONT icomp)
  1517. ************        ' Point Gauss' igau 'Etat' ietat ;
  1518. ************Message ;
  1519. ************FINSI ;
  1520.             ierr = 40 ; errabs = smxyu2.igau ; QUITTER BECART ;
  1521.          FINSI ;
  1522.       SINON ;
  1523.          dsmxy = ABS (smxyu1.igau - smxyu2.igau) ;
  1524.          errrel = ABS (dsmxy / smxyu1.igau) ;
  1525.          SI (errrel > critere) ;
  1526. ************SI (EGA IMESS 1) ;
  1527. ************Message 'Erreur relative' errrel
  1528. ************        'Composante ' (EXTR LCONT icomp)
  1529. ************        ' Point Gauss' igau 'Etat' ietat ;
  1530. ************Message 'smxyu1.igau =' smxyu1.igau
  1531. ************        'smxyu2.igau =' smxyu2.igau;
  1532. ************Message ;
  1533. ************FINSI ;
  1534.             ierr = 41 ; QUITTER BECART ;
  1535.          FINSI ;
  1536.       FINSI ;
  1537.    FIN BSMXY ;
  1538.  
  1539.    icomp = 5 ;
  1540.    igau = 0 ;
  1541.    REPETER BSMXZ nbgau ;
  1542.       igau = igau + 1 ;
  1543.       SI (EGA smxzu1.igau 0.0) ;
  1544.          SI (NON (EGA smxzu2.igau 0.0)) ;
  1545.             SI (EGA IMESS 1) ;
  1546.             Message 'Erreur absolue' smxzu2.igau
  1547.                     'Composante ' (EXTR LCONT icomp)
  1548.                     ' Point Gauss' igau 'Etat' ietat ;
  1549.             Message ;
  1550.             FINSI ;
  1551. ************ierr = 40 ; errabs = smxzu2.igau ; QUITTER BECART ;
  1552.          FINSI ;
  1553.       SINON ;
  1554.          dsmxz = ABS (smxzu1.igau - smxzu2.igau) ;
  1555.          errrel = ABS (dsmxz / smxzu1.igau) ;
  1556.          SI (errrel > critere) ;
  1557.             SI (EGA IMESS 1) ;
  1558.             Message 'Erreur relative' errrel
  1559.                     'Composante ' (EXTR LCONT icomp)
  1560.                     ' Point Gauss' igau 'Etat' ietat ;
  1561.             Message 'smxzu1.igau =' smxzu1.igau
  1562.                     'smxzu2.igau =' smxzu2.igau;
  1563.             Message ;
  1564.             FINSI ;
  1565. ************ierr = 41 ; QUITTER BECART ;
  1566.          FINSI ;
  1567.       FINSI ;
  1568.    FIN BSMXZ ;
  1569.  
  1570.    icomp = 6 ;
  1571.    igau = 0 ;
  1572.    REPETER BSMYZ nbgau ;
  1573.       igau = igau + 1 ;
  1574.       SI (EGA smyzu1.igau 0.0) ;
  1575.          SI (NON (EGA smyzu2.igau 0.0)) ;
  1576.             SI (EGA IMESS 1) ;
  1577.             Message 'Erreur absolue' smyzu2.igau
  1578.                     'Composante ' (EXTR LCONT icomp)
  1579.                     ' Point Gauss' igau 'Etat' ietat ;
  1580.             Message ;
  1581.             FINSI ;
  1582. ************ierr = 40 ; errabs = smyzu2.igau ; QUITTER BECART ;
  1583.          FINSI ;
  1584.       SINON ;
  1585.          dsmyz = ABS (smyzu1.igau - smyzu2.igau) ;
  1586.          errrel = ABS (dsmyz / smyzu1.igau) ;
  1587.          SI (errrel > critere) ;
  1588.             SI (EGA IMESS 1) ;
  1589.             Message 'Erreur relative' errrel
  1590.                     'Composante ' (EXTR LCONT icomp)
  1591.                     ' Point Gauss' igau 'Etat' ietat ;
  1592.             Message 'smyzu1.igau =' smyzu1.igau
  1593.                     'smyzu2.igau =' smyzu2.igau;
  1594.             Message ;
  1595.             FINSI ;
  1596. ************ierr = 41 ; QUITTER BECART ;
  1597.          FINSI ;
  1598.       FINSI ;
  1599.    FIN BSMYZ ;
  1600.  
  1601. *  Comparaison des variables internes-----------------------------------
  1602.  
  1603.    icomp = 1 ;
  1604.    igau = 0 ;
  1605.    REPETER BPXX nbgau ;
  1606.       igau = igau + 1 ;
  1607.       SI (EGA pxxu1.igau 0.0) ;
  1608.          SI (NON (EGA pxxu2.igau 0.0)) ;
  1609.             SI (EGA IMESS 1) ;
  1610.             Message 'Erreur absolue' pxxu2.igau
  1611.                     'Composante ' (EXTR LVARI icomp)
  1612.                     ' Point Gauss' igau 'Etat' ietat ;
  1613.             Message ;
  1614.             FINSI ;
  1615. ************ierr = 50 ; errabs = pxxu2.igau ; QUITTER BECART ;
  1616.          FINSI ;
  1617.       SINON ;
  1618.          dpxx = ABS (pxxu1.igau - pxxu2.igau) ;
  1619.          errrel = ABS (dpxx / pxxu1.igau) ;
  1620.          SI (errrel > critere) ;
  1621.             SI (EGA IMESS 1) ;
  1622.             Message 'Erreur relative' errrel
  1623.                     'Composante ' (EXTR LVARI icomp)
  1624.                     ' Point Gauss' igau 'Etat' ietat ;
  1625.             Message 'pxxu1.igau =' pxxu1.igau 'pxxu2.igau =' pxxu2.igau;
  1626.             Message ;
  1627.             FINSI ;
  1628. ************ierr = 51 ; QUITTER BECART ;
  1629.          FINSI ;
  1630.       FINSI ;
  1631.    FIN BPXX ;
  1632.  
  1633.    icomp = 2 ;
  1634.    igau = 0 ;
  1635.    REPETER BPYY nbgau ;
  1636.       igau = igau + 1 ;
  1637.       SI (EGA pyyu1.igau 0.0) ;
  1638.          SI (NON (EGA pyyu2.igau 0.0)) ;
  1639.             SI (EGA IMESS 1) ;
  1640.             Message 'Erreur absolue' pyyu2.igau
  1641.                     'Composante ' (EXTR LVARI icomp)
  1642.                     ' Point Gauss' igau 'Etat' ietat ;
  1643.             Message ;
  1644.             FINSI ;
  1645. ************ierr = 50 ; errabs = pyyu2.igau ; QUITTER BECART ;
  1646.          FINSI ;
  1647.       SINON ;
  1648.          dpyy = ABS (pyyu1.igau - pyyu2.igau) ;
  1649.          errrel = ABS (dpyy / pyyu1.igau) ;
  1650.          SI (errrel > critere) ;
  1651.             SI (EGA IMESS 1) ;
  1652.             Message 'Erreur relative' errrel
  1653.                     'Composante ' (EXTR LVARI icomp)
  1654.                     ' Point Gauss' igau 'Etat' ietat ;
  1655.             Message 'pyyu1.igau =' pyyu1.igau 'pyyu2.igau =' pyyu2.igau;
  1656.             Message ;
  1657.             FINSI ;
  1658. ************ierr = 51 ; QUITTER BECART ;
  1659.          FINSI ;
  1660.       FINSI ;
  1661.    FIN BPYY ;
  1662.  
  1663.    icomp = 3 ;
  1664.    igau = 0 ;
  1665.    REPETER BPZZ nbgau ;
  1666.       igau = igau + 1 ;
  1667.       SI (EGA pzzu1.igau 0.0) ;
  1668.          SI (NON (EGA pzzu2.igau 0.0)) ;
  1669.             SI (EGA IMESS 1) ;
  1670.             Message 'Erreur absolue' pzzu2.igau
  1671.                     'Composante ' (EXTR LVARI icomp)
  1672.                     ' Point Gauss' igau 'Etat' ietat ;
  1673.             Message ;
  1674.             FINSI ;
  1675. ************ierr = 50 ; errabs = pzzu2.igau ; QUITTER BECART ;
  1676.          FINSI ;
  1677.       SINON ;
  1678.          dpzz = ABS (pzzu1.igau - pzzu2.igau) ;
  1679.          errrel = ABS (dpzz / pzzu1.igau) ;
  1680.          SI (errrel > critere) ;
  1681.             SI (EGA IMESS 1) ;
  1682.             Message 'Erreur relative' errrel
  1683.                     'Composante ' (EXTR LVARI icomp)
  1684.                     ' Point Gauss' igau 'Etat' ietat ;
  1685.             Message 'pzzu1.igau =' pzzu1.igau 'pzzu2.igau =' pzzu2.igau;
  1686.             Message ;
  1687.             FINSI ;
  1688. ************ierr = 51 ; QUITTER BECART ;
  1689.          FINSI ;
  1690.       FINSI ;
  1691.    FIN BPZZ ;
  1692.  
  1693.    icomp = 4 ;
  1694.    igau = 0 ;
  1695.    REPETER BPXY nbgau ;
  1696.       igau = igau + 1 ;
  1697.       SI (EGA pxyu1.igau 0.0) ;
  1698.          SI (NON (EGA pxyu2.igau 0.0)) ;
  1699.             SI (EGA IMESS 1) ;
  1700.             Message 'Erreur absolue' pxyu2.igau
  1701.                     'Composante ' (EXTR LVARI icomp)
  1702.                     ' Point Gauss' igau 'Etat' ietat ;
  1703.             Message ;
  1704.             FINSI ;
  1705. ************ierr = 50 ; errabs = pxyu2.igau ; QUITTER BECART ;
  1706.          FINSI ;
  1707.       SINON ;
  1708.          dpxy = ABS (pxyu1.igau - pxyu2.igau) ;
  1709.          errrel = ABS (dpxy / pxyu1.igau) ;
  1710.          SI (errrel > critere) ;
  1711.             SI (EGA IMESS 1) ;
  1712.             Message 'Erreur relative' errrel
  1713.                     'Composante ' (EXTR LVARI icomp)
  1714.                     ' Point Gauss' igau 'Etat' ietat ;
  1715.             Message 'pxyu1.igau =' pxyu1.igau 'pxyu2.igau =' pxyu2.igau;
  1716.             Message ;
  1717.             FINSI ;
  1718. ************ierr = 51 ; QUITTER BECART ;
  1719.          FINSI ;
  1720.       FINSI ;
  1721.    FIN BPXY ;
  1722.  
  1723.    icomp = 5 ;
  1724.    igau = 0 ;
  1725.    REPETER BPYZ nbgau ;
  1726.       igau = igau + 1 ;
  1727.       SI (EGA pyzu1.igau 0.0) ;
  1728.          SI (NON (EGA pyzu2.igau 0.0)) ;
  1729.             SI (EGA IMESS 1) ;
  1730.             Message 'Erreur absolue' pyzu2.igau
  1731.                     'Composante ' (EXTR LVARI icomp)
  1732.                     ' Point Gauss' igau 'Etat' ietat ;
  1733.             Message ;
  1734.             FINSI ;
  1735. ************ierr = 50 ; errabs = pyzu2.igau ; QUITTER BECART ;
  1736.          FINSI ;
  1737.       SINON ;
  1738.          dpyz = ABS (pyzu1.igau - pyzu2.igau) ;
  1739.          errrel = ABS (dpyz / pyzu1.igau) ;
  1740.          SI (errrel > critere) ;
  1741.             SI (EGA IMESS 1) ;
  1742.             Message 'Erreur relative' errrel
  1743.                     'Composante ' (EXTR LVARI icomp)
  1744.                     ' Point Gauss' igau 'Etat' ietat ;
  1745.             Message 'pyzu1.igau =' pyzu1.igau 'pyzu2.igau =' pyzu2.igau;
  1746.             Message ;
  1747.             FINSI ;
  1748. ************ierr = 51 ; QUITTER BECART ;
  1749.          FINSI ;
  1750.       FINSI ;
  1751.    FIN BPYZ ;
  1752.  
  1753.    icomp = 6 ;
  1754.    igau = 0 ;
  1755.    REPETER BPZX nbgau ;
  1756.       igau = igau + 1 ;
  1757.       SI (EGA pzxu1.igau 0.0) ;
  1758.          SI (NON (EGA pzxu2.igau 0.0)) ;
  1759.             SI (EGA IMESS 1) ;
  1760.             Message 'Erreur absolue' pzxu2.igau
  1761.                     'Composante ' (EXTR LVARI icomp)
  1762.                     ' Point Gauss' igau 'Etat' ietat ;
  1763.             Message ;
  1764.             FINSI ;
  1765. ************ierr = 50 ; errabs = pzxu2.igau ; QUITTER BECART ;
  1766.          FINSI ;
  1767.       SINON ;
  1768.          dpzx = ABS (pzxu1.igau - pzxu2.igau) ;
  1769.          errrel = ABS (dpzx / pzxu1.igau) ;
  1770.          SI (errrel > critere) ;
  1771.             SI (EGA IMESS 1) ;
  1772.             Message 'Erreur relative' errrel
  1773.                     'Composante ' (EXTR LVARI icomp)
  1774.                     ' Point Gauss' igau 'Etat' ietat ;
  1775.             Message 'pzxu1.igau =' pzxu1.igau 'pzxu2.igau =' pzxu2.igau;
  1776.             Message ;
  1777.             FINSI ;
  1778. ************ierr = 51 ; QUITTER BECART ;
  1779.          FINSI ;
  1780.       FINSI ;
  1781.    FIN BPZX ;
  1782.  
  1783.    icomp = 7 ;
  1784.    igau = 0 ;
  1785.    REPETER BVNXX nbgau ;
  1786.       igau = igau + 1 ;
  1787.       SI (EGA vnxxu1.igau 0.0) ;
  1788.          SI (NON (EGA vnxxu2.igau 0.0)) ;
  1789.             SI (EGA IMESS 1) ;
  1790.             Message 'Erreur absolue' vnxxu2.igau
  1791.                     'Composante ' (EXTR LVARI icomp)
  1792.                     ' Point Gauss' igau 'Etat' ietat ;
  1793.             Message ;
  1794.             FINSI ;
  1795. ************ierr = 50 ; errabs = vnxxu2.igau ; QUITTER BECART ;
  1796.          FINSI ;
  1797.       SINON ;
  1798.          dvnxx = ABS (vnxxu1.igau - vnxxu2.igau) ;
  1799.          errrel = ABS (dvnxx / vnxxu1.igau) ;
  1800.          SI (errrel > critere) ;
  1801.             SI (EGA IMESS 1) ;
  1802.             Message 'Erreur relative' errrel
  1803.                     'Composante ' (EXTR LVARI icomp)
  1804.                     ' Point Gauss' igau 'Etat' ietat ;
  1805.             Message 'vnxxu1.igau =' vnxxu1.igau
  1806.                     'vnxxu2.igau =' vnxxu2.igau;
  1807.             Message ;
  1808.             FINSI ;
  1809. ************ierr = 51 ; QUITTER BECART ;
  1810.          FINSI ;
  1811.       FINSI ;
  1812.    FIN BVNXX ;
  1813.  
  1814.    icomp = 8 ;
  1815.    igau = 0 ;
  1816.    REPETER BVNYY nbgau ;
  1817.       igau = igau + 1 ;
  1818.       SI (EGA vnyyu1.igau 0.0) ;
  1819.          SI (NON (EGA vnyyu2.igau 0.0)) ;
  1820.             SI (EGA IMESS 1) ;
  1821.             Message 'Erreur absolue' vnyyu2.igau
  1822.                     'Composante ' (EXTR LVARI icomp)
  1823.                     ' Point Gauss' igau 'Etat' ietat ;
  1824.             Message ;
  1825.             FINSI ;
  1826. ************ierr = 50 ; errabs = vnyyu2.igau ; QUITTER BECART ;
  1827.          FINSI ;
  1828.       SINON ;
  1829.          dvnyy = ABS (vnyyu1.igau - vnyyu2.igau) ;
  1830.          errrel = ABS (dvnyy / vnyyu1.igau) ;
  1831.          SI (errrel > critere) ;
  1832.             SI (EGA IMESS 1) ;
  1833.             Message 'Erreur relative' errrel
  1834.                     'Composante ' (EXTR LVARI icomp)
  1835.                     ' Point Gauss' igau 'Etat' ietat ;
  1836.             Message 'vnyyu1.igau =' vnyyu1.igau
  1837.                     'vnyyu2.igau =' vnyyu2.igau;
  1838.             Message ;
  1839.             FINSI ;
  1840. ************ierr = 51 ; QUITTER BECART ;
  1841.          FINSI ;
  1842.       FINSI ;
  1843.    FIN BVNYY ;
  1844.  
  1845.    icomp = 9 ;
  1846.    igau = 0 ;
  1847.    REPETER BVNZZ nbgau ;
  1848.       igau = igau + 1 ;
  1849.       SI (EGA vnzzu1.igau 0.0) ;
  1850.          SI (NON (EGA vnzzu2.igau 0.0)) ;
  1851.             SI (EGA IMESS 1) ;
  1852.             Message 'Erreur absolue' vnzzu2.igau
  1853.                     'Composante ' (EXTR LVARI icomp)
  1854.                     ' Point Gauss' igau 'Etat' ietat ;
  1855.             Message ;
  1856.             FINSI ;
  1857. ************ierr = 50 ; errabs = vnzzu2.igau ; QUITTER BECART ;
  1858.          FINSI ;
  1859.       SINON ;
  1860.          dvnzz = ABS (vnzzu1.igau - vnzzu2.igau) ;
  1861.          errrel = ABS (dvnzz / vnzzu1.igau) ;
  1862.          SI (errrel > critere) ;
  1863.             SI (EGA IMESS 1) ;
  1864.             Message 'Erreur relative' errrel
  1865.                     'Composante ' (EXTR LVARI icomp)
  1866.                     ' Point Gauss' igau 'Etat' ietat ;
  1867.             Message 'vnzzu1.igau =' vnzzu1.igau
  1868.                     'vnzzu2.igau =' vnzzu2.igau;
  1869.             Message ;
  1870.             FINSI ;
  1871. ************ierr = 51 ; QUITTER BECART ;
  1872.          FINSI ;
  1873.       FINSI ;
  1874.    FIN BVNZZ ;
  1875.  
  1876.    icomp = 10 ;
  1877.    igau = 0 ;
  1878.    REPETER BVNXY nbgau ;
  1879.       igau = igau + 1 ;
  1880.       SI (EGA vnxyu1.igau 0.0) ;
  1881.          SI (NON (EGA vnxyu2.igau 0.0)) ;
  1882.             SI (EGA IMESS 1) ;
  1883.             Message 'Erreur absolue' vnxyu2.igau
  1884.                     'Composante ' (EXTR LVARI icomp)
  1885.                     ' Point Gauss' igau 'Etat' ietat ;
  1886.             Message ;
  1887.             FINSI ;
  1888. ************ierr = 50 ; errabs = vnxyu2.igau ; QUITTER BECART ;
  1889.          FINSI ;
  1890.       SINON ;
  1891.          dvnxy = ABS (vnxyu1.igau - vnxyu2.igau) ;
  1892.          errrel = ABS (dvnxy / vnxyu1.igau) ;
  1893.          SI (errrel > critere) ;
  1894.             SI (EGA IMESS 1) ;
  1895.             Message 'Erreur relative' errrel
  1896.                     'Composante ' (EXTR LVARI icomp)
  1897.                     ' Point Gauss' igau 'Etat' ietat ;
  1898.             Message 'vnxyu1.igau =' vnxyu1.igau
  1899.                     'vnxyu2.igau =' vnxyu2.igau;
  1900.             Message ;
  1901.             FINSI ;
  1902. ************ierr = 51 ; QUITTER BECART ;
  1903.          FINSI ;
  1904.       FINSI ;
  1905.    FIN BVNXY ;
  1906.  
  1907.    icomp = 11 ;
  1908.    igau = 0 ;
  1909.    REPETER BVNYZ nbgau ;
  1910.       igau = igau + 1 ;
  1911.       SI (EGA vnyzu1.igau 0.0) ;
  1912.          SI (NON (EGA vnyzu2.igau 0.0)) ;
  1913.             SI (EGA IMESS 1) ;
  1914.             Message 'Erreur absolue' vnyzu2.igau
  1915.                     'Composante ' (EXTR LVARI icomp)
  1916.                     ' Point Gauss' igau 'Etat' ietat ;
  1917.             Message ;
  1918.             FINSI ;
  1919. ************ierr = 50 ; errabs = vnyzu2.igau ; QUITTER BECART ;
  1920.          FINSI ;
  1921.       SINON ;
  1922.          dvnyz = ABS (vnyzu1.igau - vnyzu2.igau) ;
  1923.          errrel = ABS (dvnyz / vnyzu1.igau) ;
  1924.          SI (errrel > critere) ;
  1925.             SI (EGA IMESS 1) ;
  1926.             Message 'Erreur relative' errrel
  1927.                     'Composante ' (EXTR LVARI icomp)
  1928.                     ' Point Gauss' igau 'Etat' ietat ;
  1929.             Message 'vnyzu1.igau =' vnyzu1.igau
  1930.                     'vnyzu2.igau =' vnyzu2.igau;
  1931.             Message ;
  1932.             FINSI ;
  1933. ************ierr = 51 ; QUITTER BECART ;
  1934.          FINSI ;
  1935.       FINSI ;
  1936.    FIN BVNYZ ;
  1937.  
  1938.    icomp = 12 ;
  1939.    igau = 0 ;
  1940.    REPETER BVNZX nbgau ;
  1941.       igau = igau + 1 ;
  1942.       SI (EGA vnzxu1.igau 0.0) ;
  1943.          SI (NON (EGA vnzxu2.igau 0.0)) ;
  1944.             SI (EGA IMESS 1) ;
  1945.             Message 'Erreur absolue' vnzxu2.igau
  1946.                     'Composante ' (EXTR LVARI icomp)
  1947.                     ' Point Gauss' igau 'Etat' ietat ;
  1948.             Message ;
  1949.             FINSI ;
  1950. ************ierr = 50 ; errabs = vnzxu2.igau ; QUITTER BECART ;
  1951.          FINSI ;
  1952.       SINON ;
  1953.          dvnzx = ABS (vnzxu1.igau - vnzxu2.igau) ;
  1954.          errrel = ABS (dvnzx / vnzxu1.igau) ;
  1955.          SI (errrel > critere) ;
  1956.             SI (EGA IMESS 1) ;
  1957.             Message 'Erreur relative' errrel
  1958.                     'Composante ' (EXTR LVARI icomp)
  1959.                     ' Point Gauss' igau 'Etat' ietat ;
  1960.             Message 'vnzxu1.igau =' vnzxu1.igau
  1961.                     'vnzxu2.igau =' vnzxu2.igau;
  1962.             Message ;
  1963.             FINSI ;
  1964. ************ierr = 51 ; QUITTER BECART ;
  1965.          FINSI ;
  1966.       FINSI ;
  1967.    FIN BVNZX ;
  1968.  
  1969.    icomp = 13 ;
  1970.    igau = 0 ;
  1971.    REPETER BEPCE nbgau ;
  1972.       igau = igau + 1 ;
  1973.       SI (EGA epceu1.igau 0.0) ;
  1974.          SI (NON (EGA epceu2.igau 0.0)) ;
  1975.             SI (EGA IMESS 1) ;
  1976.             Message 'Erreur absolue' epceu2.igau
  1977.                     'Composante ' (EXTR LVARI icomp)
  1978.                     ' Point Gauss' igau 'Etat' ietat ;
  1979.             Message ;
  1980.             FINSI ;
  1981. ************ierr = 50 ; errabs = epceu2.igau ; QUITTER BECART ;
  1982.          FINSI ;
  1983.       SINON ;
  1984.          depce = ABS (epceu1.igau - epceu2.igau) ;
  1985.          errrel = ABS (depce / epceu1.igau) ;
  1986.          SI (errrel > critere) ;
  1987.             SI (EGA IMESS 1) ;
  1988.             Message 'Erreur relative' errrel
  1989.                     'Composante ' (EXTR LVARI icomp)
  1990.                     ' Point Gauss' igau 'Etat' ietat ;
  1991.             Message 'epceu1.igau =' epceu1.igau
  1992.                     'epceu2.igau =' epceu2.igau;
  1993.             Message ;
  1994.             FINSI ;
  1995. ************ierr = 51 ; QUITTER BECART ;
  1996.          FINSI ;
  1997.       FINSI ;
  1998.    FIN BEPCE ;
  1999.  
  2000.    icomp = 14 ;
  2001.    igau = 0 ;
  2002.    REPETER BEPSH nbgau ;
  2003.       igau = igau + 1 ;
  2004.       SI (EGA epshu1.igau 0.0) ;
  2005.          SI (NON (EGA epshu2.igau 0.0)) ;
  2006.             SI (EGA IMESS 1) ;
  2007.             Message 'Erreur absolue' epshu2.igau
  2008.                     'Composante ' (EXTR LVARI icomp)
  2009.                     ' Point Gauss' igau 'Etat' ietat ;
  2010.             Message ;
  2011.             FINSI ;
  2012. ************ierr = 50 ; errabs = epshu2.igau ; QUITTER BECART ;
  2013.          FINSI ;
  2014.       SINON ;
  2015.          depsh = ABS (epshu1.igau - epshu2.igau) ;
  2016.          errrel = ABS (depsh / epshu1.igau) ;
  2017.          SI (errrel > critere) ;
  2018.             SI (EGA IMESS 1) ;
  2019.             Message 'Erreur relative' errrel
  2020.                     'Composante ' (EXTR LVARI icomp)
  2021.                     ' Point Gauss' igau 'Etat' ietat ;
  2022.             Message 'epshu1.igau =' epshu1.igau
  2023.                     'epshu2.igau =' epshu2.igau;
  2024.             Message ;
  2025.             FINSI ;
  2026. ************ierr = 51 ; QUITTER BECART ;
  2027.          FINSI ;
  2028.       FINSI ;
  2029.    FIN BEPSH ;
  2030.  
  2031.    icomp = 15 ;
  2032.    igau = 0 ;
  2033.    REPETER BLL nbgau ;
  2034.       igau = igau + 1 ;
  2035.       SI (EGA llu1.igau 0.0) ;
  2036.          SI (NON (EGA llu2.igau 0.0)) ;
  2037.             SI (EGA IMESS 1) ;
  2038.             Message 'Erreur absolue' llu2.igau
  2039.                     'Composante ' (EXTR LVARI icomp)
  2040.                     ' Point Gauss' igau 'Etat' ietat ;
  2041.             Message ;
  2042.             FINSI ;
  2043. ************ierr = 50 ; errabs = llu2.igau ; QUITTER BECART ;
  2044.          FINSI ;
  2045.       SINON ;
  2046.          dll = ABS (llu1.igau - llu2.igau) ;
  2047.          errrel = ABS (dll / llu1.igau) ;
  2048.          SI (errrel > critere) ;
  2049.             SI (EGA IMESS 1) ;
  2050.             Message 'Erreur relative' errrel
  2051.                     'Composante ' (EXTR LVARI icomp)
  2052.                     ' Point Gauss' igau 'Etat' ietat ;
  2053.             Message 'llu1.igau =' llu1.igau 'llu2.igau =' llu2.igau;
  2054.             Message ;
  2055.             FINSI ;
  2056. ************ierr = 51 ; QUITTER BECART ;
  2057.          FINSI ;
  2058.       FINSI ;
  2059.    FIN BLL ;
  2060.  
  2061.    icomp = 16 ;
  2062.    igau = 0 ;
  2063.    REPETER BEPSE nbgau ;
  2064.       igau = igau + 1 ;
  2065.       SI (EGA epseu1.igau 0.0) ;
  2066.          SI (NON (EGA epseu2.igau 0.0)) ;
  2067.             SI (EGA IMESS 1) ;
  2068.             Message 'Erreur absolue' epseu2.igau
  2069.                     'Composante ' (EXTR LVARI icomp)
  2070.                     ' Point Gauss' igau 'Etat' ietat ;
  2071.             Message ;
  2072.             FINSI ;
  2073. ************ierr = 50 ; errabs = epseu2.igau ; QUITTER BECART ;
  2074.          FINSI ;
  2075.       SINON ;
  2076.          depse = ABS (epseu1.igau - epseu2.igau) ;
  2077.          errrel = ABS (depse / epseu1.igau) ;
  2078.          SI (errrel > critere) ;
  2079.             SI (EGA IMESS 1) ;
  2080.             Message 'Erreur relative' errrel
  2081.                     'Composante ' (EXTR LVARI icomp)
  2082.                     ' Point Gauss' igau 'Etat' ietat ;
  2083.             Message 'epseu1.igau =' epseu1.igau
  2084.                     'epseu2.igau =' epseu2.igau;
  2085.             Message ;
  2086.             FINSI ;
  2087. ************ierr = 51 ; QUITTER BECART ;
  2088.          FINSI ;
  2089.       FINSI ;
  2090.    FIN BEPSE ;
  2091.  
  2092. FIN BECART ;
  2093.  
  2094. * FIN   BOUCLE COMPARAISON DES RESULTATS________________________________
  2095.  
  2096. * MESSAGES D'ERREUR_____________________________________________________
  2097.  
  2098. Message ;
  2099. Message 'CPU avec modele externe 21 = ' CPUext1 'centiemes de seconde' ;
  2100. Message 'CPU avec modele externe 22 = ' CPUext2 'centiemes de seconde' ;
  2101.  
  2102. ifin = netat-1 ;
  2103. gnxy1max = 0.0 ;
  2104. igau = 0 ;
  2105. REPETER BGNXYF1 nbgau ;
  2106.    igau = igau + 1 ;
  2107.    gnxy = EXTR (tabuti1.'VARIABLES_INTERNES'.ifin)
  2108.                     'GNXY' 1 1 igau ;
  2109.    SI (gnxy > gnxy1max) ; gnxy1max = gnxy ; FINSI ;
  2110. FIN BGNXYF1 ;
  2111. gnxy2max = 0.0 ;
  2112. igau = 0 ;
  2113. REPETER BGNXYF2 nbgau ;
  2114.    igau = igau + 1 ;
  2115.    gnxy = EXTR (tabuti2.'VARIABLES_INTERNES'.ifin)
  2116.                     'GNXY' 1 1 igau ;
  2117.    SI (gnxy > gnxy2max) ; gnxy2max = gnxy ; FINSI ;
  2118. FIN BGNXYF2 ;
  2119. Message ;
  2120. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  2121. Message 'GNXY max avec modele externe 21 = ' gnxy1max ;
  2122. Message 'GNXY max avec modele externe 22 = ' gnxy2max ;
  2123. Message ;
  2124.  
  2125. SI (ierr > 0) ;
  2126.  
  2127.    SI (EGA ierr 1) ;
  2128.       Message 'Instants de calcul non coincidants' ;
  2129.       ERRE 5 ;
  2130.    FINSI ;
  2131.  
  2132.    SI (EGA ierr 10) ;
  2133.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  2134.               ' Noeud' ino 'Etat' ietat ;
  2135.       ERRE 5 ;
  2136.    FINSI ;
  2137.  
  2138.    SI (EGA ierr 11) ;
  2139.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  2140.               ' Noeud' ino 'Etat' ietat ;
  2141.       ERRE 5 ;
  2142.    FINSI ;
  2143.  
  2144.    SI (EGA ierr 20) ;
  2145.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  2146.               ' Point Gauss' igau 'Etat' ietat ;
  2147.       ERRE 5 ;
  2148.    FINSI ;
  2149.  
  2150.    SI (EGA ierr 21) ;
  2151.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  2152.               ' Point Gauss' igau 'Etat' ietat ;
  2153.       ERRE 5 ;
  2154.    FINSI ;
  2155.  
  2156.    SI (EGA ierr 30) ;
  2157.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDENN icomp)
  2158.               ' Point Gauss' igau 'Etat' ietat ;
  2159.       ERRE 5 ;
  2160.    FINSI ;
  2161.  
  2162.    SI (EGA ierr 31) ;
  2163.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDENN icomp)
  2164.               ' Point Gauss' igau 'Etat' ietat ;
  2165.       ERRE 5 ;
  2166.    FINSI ;
  2167.  
  2168.    SI (EGA ierr 40) ;
  2169.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  2170.               ' Point Gauss' igau 'Etat' ietat ;
  2171.       ERRE 5 ;
  2172.    FINSI ;
  2173.  
  2174.    SI (EGA ierr 41) ;
  2175.       Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  2176.               ' Point Gauss' igau 'Etat' ietat ;
  2177.       ERRE 5 ;
  2178.    FINSI ;
  2179.  
  2180.    SI (EGA ierr 50) ;
  2181.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  2182.               ' Point Gauss' igau 'Etat' ietat ;
  2183.       ERRE 5 ;
  2184.    FINSI ;
  2185.  
  2186.    SI (EGA ierr 51) ;
  2187.       Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  2188.               ' Point Gauss' igau 'Etat' ietat ;
  2189.       ERRE 5 ;
  2190.    FINSI ;
  2191.  
  2192. SINON ;
  2193.  
  2194.    Message 'Resultats concordants entre les deux modeles' ;
  2195.    Message 'Aucune erreur relative superieure a ' critere
  2196.            ' pour les composantes' ;
  2197.    Message (EXTR LDEFO idefo) '  ' (EXTR LDENN idenn)
  2198.                               '  ' (EXTR LCONT icont) ;
  2199.    ERRE 0 ;
  2200.  
  2201. FINSI ;
  2202.  
  2203. *-----------------------------------------------------------------------
  2204. * 4 - Test supplementaire de sauvegarde
  2205. *     BUG wrmodl.eso + limodl.eso corrige le 06/02/2003
  2206.  
  2207. ***opti sauv 'umat03_cisXY_2122.sauv' ;
  2208. ***sauv tabuti1 tabuti2 ;
  2209.  
  2210. ***Message ;
  2211. ***Message 'Sauvegarde des deux tables produites par PASAPAS' ;
  2212. ***ERRE 0 ;
  2213.  
  2214. FIN ;
  2215.  
  2216.  
  2217.  
  2218.  
  2219.  
  2220.  
  2221.  
  2222.  
  2223.  
  2224.  
  2225.  
  2226.  
  2227.  
  2228.  
  2229.  
  2230.  

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