Télécharger la plaquette Cast3M


fr en

Télécharger umat03_cisXY_2122b.dgibi

Retour à la liste

Numérotation des lignes : 

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

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