Télécharger la plaquette Cast3M


fr en

Télécharger umat03_traXX.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  umat03_traXX.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. *            TRACTION XX
  16. *
  17. *=======================================================================
  18.  
  19.  
  20. opti echo 0 ;
  21.  
  22. *=======================================================================
  23. *
  24. *                       PROCEDURES DEPENDANTES
  25. *
  26. *=======================================================================
  27.  
  28.  
  29.  
  30.  
  31. *=======================================================================
  32. *
  33. *                      JEU DE DONNEES PRINCIPAL
  34. *
  35. *=======================================================================
  36.  
  37. *-----------------------------------------------------------------------
  38. * 0 - Maillage
  39. *     PLAQUE CARREE    COTE : 1 m    EPAISSEUR : 1 cm
  40.  
  41. opti dime 3 elem cub8 mode tridim ;
  42.  
  43. l1 = 1.0 ;
  44. l2 = 1.0 ;
  45. ep0 = 0.01 ;
  46.  
  47. p01 = 0.0 0.0 0.0 ;
  48. p02 = ep0 0.0 0.0 ;
  49. p03 = ep0 l2  0.0 ;
  50. p04 = 0.0 l2  0.0 ;
  51.  
  52. p10 = 0.0 0.0 l1 ;
  53. p20 = ep0 0.0 l1 ;
  54. p30 = ep0 l2  l1 ;
  55. p40 = 0.0 l2  l1 ;
  56.  
  57. li01 = 'DROIT' 1 p01 p02 ;
  58. li02 = 'DROIT' 1 p02 p03 ;
  59. li03 = 'DROIT' 1 p03 p04 ;
  60. li04 = 'DROIT' 1 p04 p01 ;
  61.  
  62. li10 = 'DROIT' 1 p10 p20 ;
  63. li20 = 'DROIT' 1 p20 p30 ;
  64. li30 = 'DROIT' 1 p30 p40 ;
  65. li40 = 'DROIT' 1 p40 p10 ;
  66.  
  67. sur01 = 'DALL' li01 li02 li03 li04 ;
  68. sur10 = 'DALL' li10 li20 li30 li40 ;
  69. mail1 = sur01 'VOLU' 1 sur10 ;
  70.  
  71. li2  = 'DROIT' 1 p01 p10 ;
  72. li3  = 'DROIT' 1 p10 p40 ;
  73. li4  = 'DROIT' 1 p40 p04 ;
  74. sbloq = 'DALL' li04 li2 li3 li4 ;
  75.  
  76. li5  = 'DROIT' 1 p03 p02 ;
  77. li6  = 'DROIT' 1 p02 p20 ;
  78. li7  = 'DROIT' 1 p20 p30 ;
  79. li8  = 'DROIT' 1 p30 p03 ;
  80. schar = 'DALL' li5 li6 li7 li8 ;
  81.  
  82. ELIM (mail1 ET sbloq ET schar) 1.e-5 ;
  83.  
  84. *-----------------------------------------------------------------------
  85. * Comportement 'FLUAGE' 'POLYNOMIAL'
  86. * Chargement en forces imposees, CISAILLEMENT PUR
  87. *-----------------------------------------------------------------------
  88. * 1 - Calcul avec le modele 'FLUAGE' 'POLYNOMIAL' de CAST3M
  89.  
  90. * 1.1 Definition du modele
  91.  
  92. mo_flupo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  93.                       'FLUAGE' 'POLYNOMIAL' ;
  94.  
  95. * 1.2 Affectation des constantes de materiau
  96.  
  97. densiu   = 0.95 ;
  98. dgrain   = 8.0E-6 ;
  99. tuo2     = 1500.0 ;
  100. tauxfisk = 1.5E19 ;
  101.  
  102. xyoun = (2.2693E11 - (1.5399E7*tuo2)) - (9.35973E3 * (tuo2**2.0)) ;
  103. xnu   = 0.34 ;
  104. xrho  = 10950.0 ;
  105. xalph = ( (9.828E-6 - (6.390E-10*tuo2)) + (1.330E-12 * (tuo2**2.0)) )
  106.         - (1.757E-17 * (tuo2**3.0)) ;
  107.  
  108. xsmax = xyoun / 1000.0 ;
  109. xaf0  = 0.0 ;
  110. xaf1  = (EXP(0.183*(100.0-densiu)))
  111.       * ( (0.2031E-13+(0.67861E-32*tauxfisk)) / (dgrain**2.0) )
  112.       * (EXP(-45294.4*(tuo2**-1.0))) ;
  113. xaf2  = 1.0 ;
  114. xaf3  = (EXP(0.241*(100.0-densiu)))
  115.       * (0.1524E-25+(0.50825E-44*tauxfisk))
  116.       * (EXP(-66431.8*(tuo2**-1.0))) ;
  117. xaf4  = 4.5 ;
  118. xaf5  = (1.447E-34*tauxfisk) * (EXP(-3624.0*(tuo2**-1.0))) ;
  119. xaf6  = 1.0 ;
  120.  
  121. ma_flupo = MATE mo_flupo 'YOUN' xyoun 'NU  ' xnu   'RHO ' xrho
  122.                          'ALPH' xalph
  123.                          'SMAX' xsmax 'AF0 ' xaf0  'AF1 ' xaf1
  124.                          'AF2 ' xaf2  'AF3 ' xaf3  'AF4 ' xaf4
  125.                          'AF5 ' xaf5  'AF6 ' xaf6 ;
  126.  
  127. * 1.3 Definition des conditions aux limites et chargements
  128.  
  129. cl1 = BLOQ mail1 'UY' 'UZ' ;
  130. cl2 = BLOQ sbloq 'UX' ;
  131. cl0 = cl1 ET cl2 ;
  132.  
  133. ********depmax = 5.0E5 ;
  134. depmax = 3.0E10 ;
  135. dep0 = FORCE schar 'FX' depmax 'FY' 0.0 'FZ' 0.0 ;
  136. ev0 = EVOL MANU temps (PROG 0.0 1.E-2 1000.0)
  137.                     y (PROG 0.0 1.0   1.0   ) ;
  138. cha0 = CHAR MECA dep0 ev0 ;
  139.  
  140. * 1.4 Definition des instants de calcul
  141.  
  142. pro0 = PROG 0.0 PAS 2.E-4 1.E-2 ;
  143. pro1 = PROG 0.0 PAS 2.E-4 1.E-2 ;
  144.  
  145. * 1.5 Resolution par PASAPAS
  146.  
  147. tabflup = TABLE ;
  148. tabflup.'MODELE' = mo_flupo ;
  149. tabflup.'CARACTERISTIQUES' = ma_flupo ;
  150. tabflup.'CHARGEMENT' = cha0 ;
  151. tabflup.'BLOCAGES_MECANIQUES' = cl0 ;
  152. tabflup.'TEMPS_CALCULES' = pro0 ;
  153. tabflup.'TEMPS_SAUVES' = pro1 ;
  154. TMASAU=table;
  155. tabflup . 'MES_SAUVEGARDES'=TMASAU;
  156. tabflup . 'PRECISION'=1E-8;
  157. TMASAU .'DEFTO'=VRAI;
  158. TMASAU .'DEFIN'=VRAI;
  159. TEMPS 'ZERO' ;
  160. tabflup.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE';
  161. PASAPAS tabflup ;
  162. TABTPS = TEMP 'NOEC';
  163. CPUint = TABTPS.'TEMPS_CPU'.'INITIAL';
  164.  
  165. *-----------------------------------------------------------------------
  166. * 2 - Calcul avec modele utilisateur
  167.  
  168. * 2.1 Definition du modele utilisateur
  169.  
  170. LCMAT21 = MOTS 'YOUN' 'NU  ' 'RHO ' 'ALPH'
  171.                'SMAX' 'AF0 ' 'AF1 ' 'AF2 ' 'AF3 ' 'AF4 ' 'AF5 ' 'AF6 ' ;
  172.  
  173. LCVAR21 = MOTS 'PXX ' 'PYY ' 'PZZ ' 'PXY ' 'PYZ ' 'PZX '
  174.                'VNXX' 'VNYY' 'VNZZ' 'VNXY' 'VNYZ' 'VNZX'
  175.                'EPCE' 'EPSH' 'LL  ' 'EPSE'
  176.                'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  177.  
  178. mo_util = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  179.                      'NON_LINEAIRE' 'UTILISATEUR'
  180.                      'NUME_LOI' 21
  181.                      'C_MATERIAU' LCMAT21 'C_VARINTER' LCVAR21 ;
  182.  
  183. * 2.2 Affectation des proprietes de materiau
  184.  
  185. ma_util = MATE mo_util 'YOUN' xyoun 'NU  ' xnu   'RHO ' xrho
  186.                        'ALPH' xalph
  187.                        'SMAX' xsmax 'AF0 ' xaf0  'AF1 ' xaf1
  188.                        'AF2 ' xaf2  'AF3 ' xaf3  'AF4 ' xaf4
  189.                        'AF5 ' xaf5  'AF6 ' xaf6 ;
  190.  
  191. * 2.3 Definition des conditions aux limites et chargements : IDEM 1.3
  192.  
  193. * 2.4 Definition des instants de calcul : IDEM 1.4
  194.  
  195. * 2.5 Initialisation des variables internes
  196.  
  197. chpxx0  = MANU CHML mo_util 'PXX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  198. chpyy0  = MANU CHML mo_util 'PYY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  199. chpzz0  = MANU CHML mo_util 'PZZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  200. chpxy0  = MANU CHML mo_util 'PXY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  201. chpyz0  = MANU CHML mo_util 'PYZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  202. chpzx0  = MANU CHML mo_util 'PZX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  203.  
  204. chvnxx0 = MANU CHML mo_util 'VNXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  205. chvnyy0 = MANU CHML mo_util 'VNYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  206. chvnzz0 = MANU CHML mo_util 'VNZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  207. chvnxy0 = MANU CHML mo_util 'VNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  208. chvnyz0 = MANU CHML mo_util 'VNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  209. chvnzx0 = MANU CHML mo_util 'VNZX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  210.  
  211. chepce0 = MANU CHML mo_util 'EPCE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  212. chepsh0 = MANU CHML mo_util 'EPSH' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  213. chll0   = MANU CHML mo_util 'LL  ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  214. chepse0 = MANU CHML mo_util 'EPSE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  215.  
  216. chenxx0 = MANU CHML mo_util 'ENXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  217. chenyy0 = MANU CHML mo_util 'ENYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  218. chenzz0 = MANU CHML mo_util 'ENZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  219. chgnxy0 = MANU CHML mo_util 'GNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  220. chgnxz0 = MANU CHML mo_util 'GNXZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  221. chgnyz0 = MANU CHML mo_util 'GNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  222.  
  223. var_ut0 = chpxx0  ET chpyy0  ET chpzz0  ET chpxy0  ET chpyz0  ET chpzx0
  224.  ET chvnxx0 ET chvnyy0 ET chvnzz0 ET chvnxy0 ET chvnyz0 ET chvnzx0
  225.  ET chepce0 ET chepsh0 ET chll0   ET chepse0
  226.  ET chenxx0 ET chenyy0 ET chenzz0 ET chgnxy0 ET chgnxz0 ET chgnyz0 ;
  227.  
  228. * 2.6 Resolution par PASAPAS
  229.  
  230. tabutil = TABLE ;
  231. tabutil.'MODELE' = mo_util ;
  232. tabutil.'CARACTERISTIQUES' = ma_util ;
  233. tabutil.'CHARGEMENT' = cha0 ;
  234. tabutil.'BLOCAGES_MECANIQUES' = cl0 ;
  235. tabutil.'VARIABLES_INTERNES' = TABLE ;
  236. tabutil.'VARIABLES_INTERNES' . 0 = var_ut0 ;
  237. tabutil.'TEMPS_CALCULES' = pro0 ;
  238. tabutil.'TEMPS_SAUVES' = pro1 ;
  239. TMASAU=table;
  240. tabutil . 'MES_SAUVEGARDES'=TMASAU;
  241. tabutil . 'PRECISION'=1E-8;
  242. TMASAU .'DEFTO'=VRAI;
  243. TMASAU .'DEFIN'=VRAI;
  244. TEMPS 'ZERO' ;
  245. tabutil.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE';
  246. PASAPAS tabutil ;
  247. TABTPS = TEMP 'NOEC';
  248. CPUext = TABTPS.'TEMPS_CPU'.'INITIAL';
  249.  
  250. *-----------------------------------------------------------------------
  251. * 3 - Post-traitement : comparaisons entre les deux modeles
  252. *
  253. *     N.B. Comme on simule un essai de traction XX
  254. *          - Critere strict (errrel < 1.d-10) sur EPXX, EIXX et SMXX
  255. *          - Pour les autres composantes on note les depassements du
  256. *            critere si impression demandee
  257.  
  258. uxc = TABLE ; uyc = TABLE ; uzc = TABLE ;
  259. epxxc = TABLE ; epyyc = TABLE ; epzzc = TABLE ;
  260. gaxyc = TABLE ; gaxzc = TABLE ; gayzc = TABLE ;
  261. eixxc = TABLE ; eiyyc = TABLE ; eizzc = TABLE ;
  262. gixyc = TABLE ; gixzc = TABLE ; giyzc = TABLE ;
  263. smxxc = TABLE ; smyyc = TABLE ; smzzc = TABLE ;
  264. smxyc = TABLE ; smxzc = TABLE ; smyzc = TABLE ;
  265. pxxc  = TABLE ; pyyc  = TABLE ; pzzc  = TABLE ;
  266. pxyc  = TABLE ; pyzc  = TABLE ; pzxc  = TABLE ;
  267. vnxxc = TABLE ; vnyyc = TABLE ; vnzzc = TABLE ;
  268. vnxyc = TABLE ; vnyzc = TABLE ; vnzxc = TABLE ;
  269. epcec = TABLE ; epshc = TABLE ; llc   = TABLE ; epsec = TABLE ;
  270.  
  271. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  272. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  273. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  274. enxxu = TABLE ; enyyu = TABLE ; enzzu = TABLE ;
  275. gnxyu = TABLE ; gnxzu = TABLE ; gnyzu = TABLE ;
  276. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  277. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  278. pxxu  = TABLE ; pyyu  = TABLE ; pzzu  = TABLE ;
  279. pxyu  = TABLE ; pyzu  = TABLE ; pzxu  = TABLE ;
  280. vnxxu = TABLE ; vnyyu = TABLE ; vnzzu = TABLE ;
  281. vnxyu = TABLE ; vnyzu = TABLE ; vnzxu = TABLE ;
  282. epceu = TABLE ; epshu = TABLE ; llu   = TABLE ; epseu = TABLE ;
  283.  
  284. LDEPL = MOTS 'UX  ' 'UY  ' 'UZ  ' ;
  285. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  286. LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  287. LDENN = MOTS 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  288. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  289. LVARI = MOTS 'PXX ' 'PYY ' 'PZZ ' 'PXY ' 'PYZ ' 'PZX '
  290.              'VNXX' 'VNYY' 'VNZZ' 'VNXY' 'VNYZ' 'VNZX'
  291.              'EPCE' 'EPSH' 'LL  ' 'EPSE' ;
  292.  
  293. idefo = 1 ;
  294. idein = 1 ;
  295. icont = 1 ;
  296.  
  297. ierr = 0 ;
  298. critere = 1.d-8 ;
  299. IMESS = 0 ;
  300. ***IMESS = 1 ;
  301.  
  302. nbno  = 8 ;
  303. nbgau = 8 ;
  304.  
  305. netat = DIME (tabflup.'TEMPS') ;
  306. ietat = -1 ;
  307.  
  308.  
  309. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  310.  
  311. REPETER BECART netat ;
  312.  
  313.    ietat = ietat + 1 ;
  314.  
  315.    temps1 = tabflup.'TEMPS'.ietat ;
  316.    temps2 = tabutil.'TEMPS'.ietat ;
  317.    SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  318.       ierr = 1 ;
  319.       QUITTER BECART ;
  320.    FINSI ;
  321.  
  322. *  Extraction des resultats de calcul : Modele interne CAST3M-----------
  323.  
  324. *  Deplacements
  325.    uxc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  326.    uxc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  327.    uxc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  328.    uxc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  329.    uxc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  330.    uxc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  331.    uxc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  332.    uxc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  333.  
  334.    uyc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  335.    uyc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  336.    uyc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  337.    uyc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  338.    uyc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  339.    uyc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  340.    uyc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  341.    uyc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  342.  
  343.    uzc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  344.    uzc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  345.    uzc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  346.    uzc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  347.    uzc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  348.    uzc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  349.    uzc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  350.    uzc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  351.  
  352. *  Deformations totales
  353.    igau = 0 ;
  354.    REPETER BEPXXC nbgau ;
  355.       igau = igau + 1 ;
  356.       epxxc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  357.    FIN BEPXXC ;
  358.  
  359.    igau = 0 ;
  360.    REPETER BEPYYC nbgau ;
  361.       igau = igau + 1 ;
  362.       epyyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  363.    FIN BEPYYC ;
  364.  
  365.    igau = 0 ;
  366.    REPETER BEPZZC nbgau ;
  367.       igau = igau + 1 ;
  368.       epzzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  369.    FIN BEPZZC ;
  370.  
  371.    igau = 0 ;
  372.    REPETER BGAXYC nbgau ;
  373.       igau = igau + 1 ;
  374.       gaxyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  375.    FIN BGAXYC ;
  376.  
  377.    igau = 0 ;
  378.    REPETER BGAXZC nbgau ;
  379.       igau = igau + 1 ;
  380.       gaxzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  381.    FIN BGAXZC ;
  382.  
  383.    igau = 0 ;
  384.    REPETER BGAYZC nbgau ;
  385.       igau = igau + 1 ;
  386.       gayzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  387.    FIN BGAYZC ;
  388.  
  389. *  Deformations inelastiques
  390.    igau = 0 ;
  391.    REPETER BEIXXC nbgau ;
  392.       igau = igau + 1 ;
  393.       eixxc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  394.                              'EIXX' 1 1 igau ;
  395.    FIN BEIXXC ;
  396.  
  397.    igau = 0 ;
  398.    REPETER BEIYYC nbgau ;
  399.       igau = igau + 1 ;
  400.       eiyyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  401.                              'EIYY' 1 1 igau ;
  402.    FIN BEIYYC ;
  403.  
  404.    igau = 0 ;
  405.    REPETER BEIZZC nbgau ;
  406.       igau = igau + 1 ;
  407.       eizzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  408.                              'EIZZ' 1 1 igau ;
  409.    FIN BEIZZC ;
  410.  
  411.    igau = 0 ;
  412.    REPETER BGIXYC nbgau ;
  413.       igau = igau + 1 ;
  414.       gixyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  415.                              'GIXY' 1 1 igau ;
  416.    FIN BGIXYC ;
  417.  
  418.    igau = 0 ;
  419.    REPETER BGIXZC nbgau ;
  420.       igau = igau + 1 ;
  421.       gixzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  422.                              'GIXZ' 1 1 igau ;
  423.    FIN BGIXZC ;
  424.  
  425.    igau = 0 ;
  426.    REPETER BGIYZC nbgau ;
  427.       igau = igau + 1 ;
  428.       giyzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  429.                              'GIYZ' 1 1 igau ;
  430.    FIN BGIYZC ;
  431.  
  432. *  Contraintes
  433.    igau = 0 ;
  434.    REPETER BSMXXC nbgau ;
  435.       igau = igau + 1 ;
  436.       smxxc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  437.    FIN BSMXXC ;
  438.  
  439.    igau = 0 ;
  440.    REPETER BSMYYC nbgau ;
  441.       igau = igau + 1 ;
  442.       smyyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  443.    FIN BSMYYC ;
  444.  
  445.    igau = 0 ;
  446.    REPETER BSMZZC nbgau ;
  447.       igau = igau + 1 ;
  448.       smzzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  449.    FIN BSMZZC ;
  450.  
  451.    igau = 0 ;
  452.    REPETER BSMXYC nbgau ;
  453.       igau = igau + 1 ;
  454.       smxyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  455.    FIN BSMXYC ;
  456.  
  457.    igau = 0 ;
  458.    REPETER BSMXZC nbgau ;
  459.       igau = igau + 1 ;
  460.       smxzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  461.    FIN BSMXZC ;
  462.  
  463.    igau = 0 ;
  464.    REPETER BSMYZC nbgau ;
  465.       igau = igau + 1 ;
  466.       smyzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  467.    FIN BSMYZC ;
  468.  
  469. *  Variables internes
  470.    igau = 0 ;
  471.    REPETER BPXXC nbgau ;
  472.       igau = igau + 1 ;
  473.       pxxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  474.                              'PXX ' 1 1 igau ;
  475.    FIN BPXXC ;
  476.  
  477.    igau = 0 ;
  478.    REPETER BPYYC nbgau ;
  479.       igau = igau + 1 ;
  480.       pyyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  481.                              'PYY ' 1 1 igau ;
  482.    FIN BPYYC ;
  483.  
  484.    igau = 0 ;
  485.    REPETER BPZZC nbgau ;
  486.       igau = igau + 1 ;
  487.       pzzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  488.                              'PZZ ' 1 1 igau ;
  489.    FIN BPZZC ;
  490.  
  491.    igau = 0 ;
  492.    REPETER BPXYC nbgau ;
  493.       igau = igau + 1 ;
  494.       pxyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  495.                              'PXY ' 1 1 igau ;
  496.    FIN BPXYC ;
  497.  
  498.    igau = 0 ;
  499.    REPETER BPYZC nbgau ;
  500.       igau = igau + 1 ;
  501.       pyzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  502.                              'PYZ ' 1 1 igau ;
  503.    FIN BPYZC ;
  504.  
  505.    igau = 0 ;
  506.    REPETER BPZXC nbgau ;
  507.       igau = igau + 1 ;
  508.       pzxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  509.                              'PZX ' 1 1 igau ;
  510.    FIN BPZXC ;
  511.  
  512.    igau = 0 ;
  513.    REPETER BVNXXC nbgau ;
  514.       igau = igau + 1 ;
  515.       vnxxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  516.                              'VNXX' 1 1 igau ;
  517.    FIN BVNXXC ;
  518.  
  519.    igau = 0 ;
  520.    REPETER BVNYYC nbgau ;
  521.       igau = igau + 1 ;
  522.       vnyyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  523.                              'VNYY' 1 1 igau ;
  524.    FIN BVNYYC ;
  525.  
  526.    igau = 0 ;
  527.    REPETER BVNZZC nbgau ;
  528.       igau = igau + 1 ;
  529.       vnzzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  530.                              'VNZZ' 1 1 igau ;
  531.    FIN BVNZZC ;
  532.  
  533.    igau = 0 ;
  534.    REPETER BVNXYC nbgau ;
  535.       igau = igau + 1 ;
  536.       vnxyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  537.                              'VNXY' 1 1 igau ;
  538.    FIN BVNXYC ;
  539.  
  540.    igau = 0 ;
  541.    REPETER BVNYZC nbgau ;
  542.       igau = igau + 1 ;
  543.       vnyzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  544.                              'VNYZ' 1 1 igau ;
  545.    FIN BVNYZC ;
  546.  
  547.    igau = 0 ;
  548.    REPETER BVNZXC nbgau ;
  549.       igau = igau + 1 ;
  550.       vnzxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  551.                              'VNZX' 1 1 igau ;
  552.    FIN BVNZXC ;
  553.  
  554.    igau = 0 ;
  555.    REPETER BEPCEC nbgau ;
  556.       igau = igau + 1 ;
  557.       epcec.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  558.                              'EPCE' 1 1 igau ;
  559.    FIN BEPCEC ;
  560.  
  561.    igau = 0 ;
  562.    REPETER BEPSHC nbgau ;
  563.       igau = igau + 1 ;
  564.       epshc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  565.                              'EPSH' 1 1 igau ;
  566.    FIN BEPSHC ;
  567.  
  568.    igau = 0 ;
  569.    REPETER BLLC nbgau ;
  570.       igau = igau + 1 ;
  571.       llc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  572.                              'LL  ' 1 1 igau ;
  573.    FIN BLLC ;
  574.  
  575.    igau = 0 ;
  576.    REPETER BEPSEC nbgau ;
  577.       igau = igau + 1 ;
  578.       epsec.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  579.                              'EPSE' 1 1 igau ;
  580.    FIN BEPSEC ;
  581.  
  582. *  Extraction des resultats de calcul : Modele externe utilisateur------
  583.  
  584. *  Deplacements
  585.    uxu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  586.    uxu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  587.    uxu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  588.    uxu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  589.    uxu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  590.    uxu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  591.    uxu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  592.    uxu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  593.  
  594.    uyu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  595.    uyu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  596.    uyu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  597.    uyu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  598.    uyu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  599.    uyu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  600.    uyu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  601.    uyu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  602.  
  603.    uzu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  604.    uzu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  605.    uzu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  606.    uzu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  607.    uzu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  608.    uzu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  609.    uzu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  610.    uzu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  611.  
  612. *  Deformations totales
  613.    igau = 0 ;
  614.    REPETER BEPXXU nbgau ;
  615.       igau = igau + 1 ;
  616.       epxxu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  617.    FIN BEPXXU ;
  618.  
  619.    igau = 0 ;
  620.    REPETER BEPYYU nbgau ;
  621.       igau = igau + 1 ;
  622.       epyyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  623.    FIN BEPYYU ;
  624.  
  625.    igau = 0 ;
  626.    REPETER BEPZZU nbgau ;
  627.       igau = igau + 1 ;
  628.       epzzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  629.    FIN BEPZZU ;
  630.  
  631.    igau = 0 ;
  632.    REPETER BGAXYU nbgau ;
  633.       igau = igau + 1 ;
  634.       gaxyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  635.    FIN BGAXYU ;
  636.  
  637.    igau = 0 ;
  638.    REPETER BGAXZU nbgau ;
  639.       igau = igau + 1 ;
  640.       gaxzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  641.    FIN BGAXZU ;
  642.  
  643.    igau = 0 ;
  644.    REPETER BGAYZU nbgau ;
  645.       igau = igau + 1 ;
  646.       gayzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  647.    FIN BGAYZU ;
  648.  
  649. *  Deformations inelastiques (variables internes supplementaires)
  650.    igau = 0 ;
  651.    REPETER BENXXU nbgau ;
  652.       igau = igau + 1 ;
  653.       enxxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  654.                               'ENXX' 1 1 igau ;
  655.    FIN BENXXU ;
  656.  
  657.    igau = 0 ;
  658.    REPETER BENYYU nbgau ;
  659.       igau = igau + 1 ;
  660.       enyyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  661.                               'ENYY' 1 1 igau ;
  662.    FIN BENYYU ;
  663.  
  664.    igau = 0 ;
  665.    REPETER BENZZU nbgau ;
  666.       igau = igau + 1 ;
  667.       enzzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  668.                               'ENZZ' 1 1 igau ;
  669.    FIN BENZZU ;
  670.  
  671.    igau = 0 ;
  672.    REPETER BGNXYU nbgau ;
  673.       igau = igau + 1 ;
  674.       gnxyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  675.                               'GNXY' 1 1 igau ;
  676.    FIN BGNXYU ;
  677.  
  678.    igau = 0 ;
  679.    REPETER BGNXZU nbgau ;
  680.       igau = igau + 1 ;
  681.       gnxzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  682.                               'GNXZ' 1 1 igau ;
  683.    FIN BGNXZU ;
  684.  
  685.    igau = 0 ;
  686.    REPETER BGNYZU nbgau ;
  687.       igau = igau + 1 ;
  688.       gnyzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  689.                               'GNYZ' 1 1 igau ;
  690.    FIN BGNYZU ;
  691.  
  692. *  Contraintes
  693.    igau = 0 ;
  694.    REPETER BSMXXU nbgau ;
  695.       igau = igau + 1 ;
  696.       smxxu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  697.    FIN BSMXXU ;
  698.  
  699.    igau = 0 ;
  700.    REPETER BSMYYU nbgau ;
  701.       igau = igau + 1 ;
  702.       smyyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  703.    FIN BSMYYU ;
  704.  
  705.    igau = 0 ;
  706.    REPETER BSMZZU nbgau ;
  707.       igau = igau + 1 ;
  708.       smzzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  709.    FIN BSMZZU ;
  710.  
  711.    igau = 0 ;
  712.    REPETER BSMXYU nbgau ;
  713.       igau = igau + 1 ;
  714.       smxyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  715.    FIN BSMXYU ;
  716.  
  717.    igau = 0 ;
  718.    REPETER BSMXZU nbgau ;
  719.       igau = igau + 1 ;
  720.       smxzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  721.    FIN BSMXZU ;
  722.  
  723.    igau = 0 ;
  724.    REPETER BSMYZU nbgau ;
  725.       igau = igau + 1 ;
  726.       smyzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  727.    FIN BSMYZU ;
  728.  
  729. *  Variables internes
  730.    igau = 0 ;
  731.    REPETER BPXXU nbgau ;
  732.       igau = igau + 1 ;
  733.       pxxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  734.                              'PXX ' 1 1 igau ;
  735.    FIN BPXXU ;
  736.  
  737.    igau = 0 ;
  738.    REPETER BPYYU nbgau ;
  739.       igau = igau + 1 ;
  740.       pyyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  741.                              'PYY ' 1 1 igau ;
  742.    FIN BPYYU ;
  743.  
  744.    igau = 0 ;
  745.    REPETER BPZZU nbgau ;
  746.       igau = igau + 1 ;
  747.       pzzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  748.                              'PZZ ' 1 1 igau ;
  749.    FIN BPZZU ;
  750.  
  751.    igau = 0 ;
  752.    REPETER BPXYU nbgau ;
  753.       igau = igau + 1 ;
  754.       pxyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  755.                              'PXY ' 1 1 igau ;
  756.    FIN BPXYU ;
  757.  
  758.    igau = 0 ;
  759.    REPETER BPYZU nbgau ;
  760.       igau = igau + 1 ;
  761.       pyzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  762.                              'PYZ ' 1 1 igau ;
  763.    FIN BPYZU ;
  764.  
  765.    igau = 0 ;
  766.    REPETER BPZXU nbgau ;
  767.       igau = igau + 1 ;
  768.       pzxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  769.                              'PZX ' 1 1 igau ;
  770.    FIN BPZXU ;
  771.  
  772.    igau = 0 ;
  773.    REPETER BVNXXU nbgau ;
  774.       igau = igau + 1 ;
  775.       vnxxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  776.                              'VNXX' 1 1 igau ;
  777.    FIN BVNXXU ;
  778.  
  779.    igau = 0 ;
  780.    REPETER BVNYYU nbgau ;
  781.       igau = igau + 1 ;
  782.       vnyyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  783.                              'VNYY' 1 1 igau ;
  784.    FIN BVNYYU ;
  785.  
  786.    igau = 0 ;
  787.    REPETER BVNZZU nbgau ;
  788.       igau = igau + 1 ;
  789.       vnzzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  790.                              'VNZZ' 1 1 igau ;
  791.    FIN BVNZZU ;
  792.  
  793.    igau = 0 ;
  794.    REPETER BVNXYU nbgau ;
  795.       igau = igau + 1 ;
  796.       vnxyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  797.                              'VNXY' 1 1 igau ;
  798.    FIN BVNXYU ;
  799.  
  800.    igau = 0 ;
  801.    REPETER BVNYZU nbgau ;
  802.       igau = igau + 1 ;
  803.       vnyzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  804.                              'VNYZ' 1 1 igau ;
  805.    FIN BVNYZU ;
  806.  
  807.    igau = 0 ;
  808.    REPETER BVNZXU nbgau ;
  809.       igau = igau + 1 ;
  810.       vnzxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  811.                              'VNZX' 1 1 igau ;
  812.    FIN BVNZXU ;
  813.  
  814.    igau = 0 ;
  815.    REPETER BEPCEU nbgau ;
  816.       igau = igau + 1 ;
  817.       epceu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  818.                              'EPCE' 1 1 igau ;
  819.    FIN BEPCEU ;
  820.  
  821.    igau = 0 ;
  822.    REPETER BEPSHU nbgau ;
  823.       igau = igau + 1 ;
  824.       epshu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  825.                              'EPSH' 1 1 igau ;
  826.    FIN BEPSHU ;
  827.  
  828.    igau = 0 ;
  829.    REPETER BLLU nbgau ;
  830.       igau = igau + 1 ;
  831.       llu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  832.                              'LL  ' 1 1 igau ;
  833.    FIN BLLU ;
  834.  
  835.    igau = 0 ;
  836.    REPETER BEPSEU nbgau ;
  837.       igau = igau + 1 ;
  838.       epseu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  839.                              'EPSE' 1 1 igau ;
  840.    FIN BEPSEU ;
  841.  
  842. *  Comparaison des deplacements-----------------------------------------
  843.  
  844.    icomp = 1 ;
  845.    ino = 0 ;
  846.    REPETER BUX nbno ;
  847.       ino = ino + 1 ;
  848.       SI (EGA uxc.ino 0.0) ;
  849.          SI (NON (EGA uxu.ino 0.0)) ;
  850.             SI (EGA IMESS 1) ;
  851.             Message 'Erreur absolue' uxu.ino
  852.                     'Composante ' (EXTR LDEPL icomp)
  853.                     ' Noeud' ino 'Etat' ietat ;
  854.             Message ;
  855.             FINSI ;
  856. ************ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  857.          FINSI ;
  858.       SINON ;
  859.          dux = ABS (uxc.ino - uxu.ino) ;
  860.          errrel = ABS (dux / uxc.ino) ;
  861.          SI (errrel > critere) ;
  862.             SI (EGA IMESS 1) ;
  863.             Message 'Erreur relative' errrel
  864.                     'Composante ' (EXTR LDEPL icomp)
  865.                     ' Noeud' ino 'Etat' ietat ;
  866.             Message 'uxc.ino =' uxc.ino 'uxu.ino =' uxu.ino ;
  867.             Message ;
  868.             FINSI ;
  869. ************ierr = 11 ; QUITTER BECART ;
  870.          FINSI ;
  871.       FINSI ;
  872.    FIN BUX ;
  873.  
  874.    icomp = 2 ;
  875.    ino = 0 ;
  876.    REPETER BUY nbno ;
  877.       ino = ino + 1 ;
  878.       SI (EGA uyc.ino 0.0) ;
  879.          SI (NON (EGA uyu.ino 0.0)) ;
  880.             SI (EGA IMESS 1) ;
  881.             Message 'Erreur absolue' uyu.ino
  882.                     'Composante ' (EXTR LDEPL icomp)
  883.                     ' Noeud' ino 'Etat' ietat ;
  884.             Message ;
  885.             FINSI ;
  886. ************ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  887.          FINSI ;
  888.       SINON ;
  889.          duy = ABS (uyc.ino - uyu.ino) ;
  890.          errrel = ABS (duy / uyc.ino) ;
  891.          SI (errrel > critere) ;
  892.             SI (EGA IMESS 1) ;
  893.             Message 'Erreur relative' errrel
  894.                     'Composante ' (EXTR LDEPL icomp)
  895.                     ' Noeud' ino 'Etat' ietat ;
  896.             Message 'uyc.ino =' uyc.ino 'uyu.ino =' uyu.ino ;
  897.             Message ;
  898.             FINSI ;
  899. ************ierr = 11 ; QUITTER BECART ;
  900.          FINSI ;
  901.       FINSI ;
  902.    FIN BUY ;
  903.  
  904.    icomp = 3 ;
  905.    ino = 0 ;
  906.    REPETER BUZ nbno ;
  907.       ino = ino + 1 ;
  908.       SI (EGA uzc.ino 0.0) ;
  909.          SI (NON (EGA uzu.ino 0.0)) ;
  910.             SI (EGA IMESS 1) ;
  911.             Message 'Erreur absolue' uzu.ino
  912.                     'Composante ' (EXTR LDEPL icomp)
  913.                     ' Noeud' ino 'Etat' ietat ;
  914.             Message ;
  915.             FINSI ;
  916. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  917.          FINSI ;
  918.       SINON ;
  919.          duz = ABS (uzc.ino - uzu.ino) ;
  920.          errrel = ABS (duz / uzc.ino) ;
  921.          SI (errrel > critere) ;
  922.             SI (EGA IMESS 1) ;
  923.             Message 'Erreur relative' errrel
  924.                     'Composante ' (EXTR LDEPL icomp)
  925.                     ' Noeud' ino 'Etat' ietat ;
  926.             Message 'uzc.ino =' uzc.ino 'uzu.ino =' uzu.ino ;
  927.             Message ;
  928.             FINSI ;
  929. ************ierr = 11 ; QUITTER BECART ;
  930.          FINSI ;
  931.       FINSI ;
  932.    FIN BUZ ;
  933.  
  934. *  Comparaison des deformations-----------------------------------------
  935.  
  936.    icomp = 1 ;
  937.    igau = 0 ;
  938.    REPETER BEPXX nbgau ;
  939.       igau = igau + 1 ;
  940.       SI (EGA epxxc.igau 0.0) ;
  941.          SI (NON (EGA epxxu.igau 0.0)) ;
  942. ************SI (EGA IMESS 1) ;
  943. ************Message 'Erreur absolue' epxxu.igau
  944. ************        'Composante ' (EXTR LDEFO icomp)
  945. ************        ' Point Gauss' igau 'Etat' ietat ;
  946. ************Message ;
  947. ************FINSI ;
  948.             ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  949.          FINSI ;
  950.       SINON ;
  951.          depxx = ABS (epxxc.igau - epxxu.igau) ;
  952.          errrel = ABS (depxx / epxxc.igau) ;
  953.          SI (errrel > critere) ;
  954. ************SI (EGA IMESS 1) ;
  955. ************Message 'Erreur relative' errrel
  956. ************        'Composante ' (EXTR LDEFO icomp)
  957. ************        ' Point Gauss' igau 'Etat' ietat ;
  958. ************Message 'epxxc.igau =' epxxc.igau 'epxxu.igau =' epxxu.igau;
  959. ************Message ;
  960. ************FINSI ;
  961.             ierr = 21 ; QUITTER BECART ;
  962.          FINSI ;
  963.       FINSI ;
  964.    FIN BEPXX ;
  965.  
  966.    icomp = 2 ;
  967.    igau = 0 ;
  968.    REPETER BEPYY nbgau ;
  969.       igau = igau + 1 ;
  970.       SI (EGA epyyc.igau 0.0) ;
  971.          SI (NON (EGA epyyu.igau 0.0)) ;
  972.             SI (EGA IMESS 1) ;
  973.             Message 'Erreur absolue' epyyu.igau
  974.                     'Composante ' (EXTR LDEFO icomp)
  975.                     ' Point Gauss' igau 'Etat' ietat ;
  976.             Message ;
  977.             FINSI ;
  978. ************ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  979.          FINSI ;
  980.       SINON ;
  981.          depyy = ABS (epyyc.igau - epyyu.igau) ;
  982.          errrel = ABS (depyy / epyyc.igau) ;
  983.          SI (errrel > critere) ;
  984.             SI (EGA IMESS 1) ;
  985.             Message 'Erreur relative' errrel
  986.                     'Composante ' (EXTR LDEFO icomp)
  987.                     ' Point Gauss' igau 'Etat' ietat ;
  988.             Message 'epyyc.igau =' epyyc.igau 'epyyu.igau =' epyyu.igau;
  989.             Message ;
  990.             FINSI ;
  991. ************ierr = 21 ; QUITTER BECART ;
  992.          FINSI ;
  993.       FINSI ;
  994.    FIN BEPYY ;
  995.  
  996.    icomp = 3 ;
  997.    igau = 0 ;
  998.    REPETER BEPZZ nbgau ;
  999.       igau = igau + 1 ;
  1000.       SI (EGA epzzc.igau 0.0) ;
  1001.          SI (NON (EGA epzzu.igau 0.0)) ;
  1002.             SI (EGA IMESS 1) ;
  1003.             Message 'Erreur absolue' epzzu.igau
  1004.                     'Composante ' (EXTR LDEFO icomp)
  1005.                     ' Point Gauss' igau 'Etat' ietat ;
  1006.             Message ;
  1007.             FINSI ;
  1008. ************ierr = 20 ; errabs = epzzu.igau ; QUITTER BECART ;
  1009.          FINSI ;
  1010.       SINON ;
  1011.          depzz = ABS (epzzc.igau - epzzu.igau) ;
  1012.          errrel = ABS (depzz / epzzc.igau) ;
  1013.          SI (errrel > critere) ;
  1014.             SI (EGA IMESS 1) ;
  1015.             Message 'Erreur relative' errrel
  1016.                     'Composante ' (EXTR LDEFO icomp)
  1017.                     ' Point Gauss' igau 'Etat' ietat ;
  1018.             Message 'epzzc.igau =' epzzc.igau 'epzzu.igau =' epzzu.igau;
  1019.             Message ;
  1020.             FINSI ;
  1021. ************ierr = 21 ; QUITTER BECART ;
  1022.          FINSI ;
  1023.       FINSI ;
  1024.    FIN BEPZZ ;
  1025.  
  1026.    icomp = 4 ;
  1027.    igau = 0 ;
  1028.    REPETER BGAXY nbgau ;
  1029.       igau = igau + 1 ;
  1030.       SI (EGA gaxyc.igau 0.0) ;
  1031.          SI (NON (EGA gaxyu.igau 0.0)) ;
  1032.             SI (EGA IMESS 1) ;
  1033.             Message 'Erreur absolue' gaxyu.igau
  1034.                     'Composante ' (EXTR LDEFO icomp)
  1035.                     ' Point Gauss' igau 'Etat' ietat ;
  1036.             Message ;
  1037.             FINSI ;
  1038. ************ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  1039.          FINSI ;
  1040.       SINON ;
  1041.          dgaxy = ABS (gaxyc.igau - gaxyu.igau) ;
  1042.          errrel = ABS (dgaxy / gaxyc.igau) ;
  1043.          SI (errrel > critere) ;
  1044.             SI (EGA IMESS 1) ;
  1045.             Message 'Erreur relative' errrel
  1046.                     'Composante ' (EXTR LDEFO icomp)
  1047.                     ' Point Gauss' igau 'Etat' ietat ;
  1048.             Message 'gaxyc.igau =' gaxyc.igau 'gaxyu.igau =' gaxyu.igau;
  1049.             Message ;
  1050.             FINSI ;
  1051. ************ierr = 21 ; QUITTER BECART ;
  1052.          FINSI ;
  1053.       FINSI ;
  1054.    FIN BGAXY ;
  1055.  
  1056.    icomp = 5 ;
  1057.    igau = 0 ;
  1058.    REPETER BGAXZ nbgau ;
  1059.       igau = igau + 1 ;
  1060.       SI (EGA gaxzc.igau 0.0) ;
  1061.          SI (NON (EGA gaxzu.igau 0.0)) ;
  1062.             SI (EGA IMESS 1) ;
  1063.             Message 'Erreur absolue' gaxzu.igau
  1064.                     'Composante ' (EXTR LDEFO icomp)
  1065.                     ' Point Gauss' igau 'Etat' ietat ;
  1066.             Message ;
  1067.             FINSI ;
  1068. ************ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  1069.          FINSI ;
  1070.       SINON ;
  1071.          dgaxz = ABS (gaxzc.igau - gaxzu.igau) ;
  1072.          errrel = ABS (dgaxz / gaxzc.igau) ;
  1073.          SI (errrel > critere) ;
  1074.             SI (EGA IMESS 1) ;
  1075.             Message 'Erreur relative' errrel
  1076.                     'Composante ' (EXTR LDEFO icomp)
  1077.                     ' Point Gauss' igau 'Etat' ietat ;
  1078.             Message 'gaxzc.igau =' gaxzc.igau 'gaxzu.igau =' gaxzu.igau;
  1079.             Message ;
  1080.             FINSI ;
  1081. ************ierr = 21 ; QUITTER BECART ;
  1082.          FINSI ;
  1083.       FINSI ;
  1084.    FIN BGAXZ ;
  1085.  
  1086.    icomp = 6 ;
  1087.    igau = 0 ;
  1088.    REPETER BGAYZ nbgau ;
  1089.       igau = igau + 1 ;
  1090.       SI (EGA gayzc.igau 0.0) ;
  1091.          SI (NON (EGA gayzu.igau 0.0)) ;
  1092.             SI (EGA IMESS 1) ;
  1093.             Message 'Erreur absolue' gayzu.igau
  1094.                     'Composante ' (EXTR LDEFO icomp)
  1095.                     ' Point Gauss' igau 'Etat' ietat ;
  1096.             Message ;
  1097.             FINSI ;
  1098. ************ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  1099.          FINSI ;
  1100.       SINON ;
  1101.          dgayz = ABS (gayzc.igau - gayzu.igau) ;
  1102.          errrel = ABS (dgayz / gayzc.igau) ;
  1103.          SI (errrel > critere) ;
  1104.             SI (EGA IMESS 1) ;
  1105.             Message 'Erreur relative' errrel
  1106.                     'Composante ' (EXTR LDEFO icomp)
  1107.                     ' Point Gauss' igau 'Etat' ietat ;
  1108.             Message 'gayzc.igau =' gayzc.igau 'gayzu.igau =' gayzu.igau;
  1109.             Message ;
  1110.             FINSI ;
  1111. ************ierr = 21 ; QUITTER BECART ;
  1112.          FINSI ;
  1113.       FINSI ;
  1114.    FIN BGAYZ ;
  1115.  
  1116. *  Comparaison des deformations inelastiques----------------------------
  1117.  
  1118.    icomp = 1 ;
  1119.    igau = 0 ;
  1120.    REPETER BEIXX nbgau ;
  1121.       igau = igau + 1 ;
  1122.       SI (EGA eixxc.igau 0.0) ;
  1123.          SI (NON (EGA enxxu.igau 0.0)) ;
  1124. ************SI (EGA IMESS 1) ;
  1125. ************Message 'Erreur absolue' enxxu.igau
  1126. ************        'Composante ' (EXTR LDEIN icomp)
  1127. ************        ' Point Gauss' igau 'Etat' ietat ;
  1128. ************Message ;
  1129. ************FINSI ;
  1130.             ierr = 30 ; errabs = enxxu.igau ; QUITTER BECART ;
  1131.          FINSI ;
  1132.       SINON ;
  1133.          deixx = ABS (eixxc.igau - enxxu.igau) ;
  1134.          errrel = ABS (deixx / eixxc.igau) ;
  1135.          SI (errrel > critere) ;
  1136. ************SI (EGA IMESS 1) ;
  1137. ************Message 'Erreur relative' errrel
  1138. ************        'Composante ' (EXTR LDEIN icomp)
  1139. ************        ' Point Gauss' igau 'Etat' ietat ;
  1140. ************Message 'eixxc.igau =' eixxc.igau 'enxxu.igau =' enxxu.igau;
  1141. ************Message ;
  1142. ************FINSI ;
  1143.             ierr = 31 ; QUITTER BECART ;
  1144.          FINSI ;
  1145.       FINSI ;
  1146.    FIN BEIXX ;
  1147.  
  1148.    icomp = 2 ;
  1149.    igau = 0 ;
  1150.    REPETER BEIYY nbgau ;
  1151.       igau = igau + 1 ;
  1152.       SI (EGA eiyyc.igau 0.0) ;
  1153.          SI (NON (EGA enyyu.igau 0.0)) ;
  1154.             SI (EGA IMESS 1) ;
  1155.             Message 'Erreur absolue' enyyu.igau
  1156.                     'Composante ' (EXTR LDEIN icomp)
  1157.                     ' Point Gauss' igau 'Etat' ietat ;
  1158.             Message ;
  1159.             FINSI ;
  1160. ************ierr = 30 ; errabs = enyyu.igau ; QUITTER BECART ;
  1161.          FINSI ;
  1162.       SINON ;
  1163.          deiyy = ABS (eiyyc.igau - enyyu.igau) ;
  1164.          errrel = ABS (deiyy / eiyyc.igau) ;
  1165.          SI (errrel > critere) ;
  1166.             SI (EGA IMESS 1) ;
  1167.             Message 'Erreur relative' errrel
  1168.                     'Composante ' (EXTR LDEIN icomp)
  1169.                     ' Point Gauss' igau 'Etat' ietat ;
  1170.             Message 'eiyyc.igau =' eiyyc.igau 'enyyu.igau =' enyyu.igau;
  1171.             Message ;
  1172.             FINSI ;
  1173. ************ierr = 31 ; QUITTER BECART ;
  1174.          FINSI ;
  1175.       FINSI ;
  1176.    FIN BEIYY ;
  1177.  
  1178.    icomp = 3 ;
  1179.    igau = 0 ;
  1180.    REPETER BEIZZ nbgau ;
  1181.       igau = igau + 1 ;
  1182.       SI (EGA eizzc.igau 0.0) ;
  1183.          SI (NON (EGA enzzu.igau 0.0)) ;
  1184.             SI (EGA IMESS 1) ;
  1185.             Message 'Erreur absolue' enzzu.igau
  1186.                     'Composante ' (EXTR LDEIN icomp)
  1187.                     ' Point Gauss' igau 'Etat' ietat ;
  1188.             Message ;
  1189.             FINSI ;
  1190. ************ierr = 30 ; errabs = enzzu.igau ; QUITTER BECART ;
  1191.          FINSI ;
  1192.       SINON ;
  1193.          deizz = ABS (eizzc.igau - enzzu.igau) ;
  1194.          errrel = ABS (deizz / eizzc.igau) ;
  1195.          SI (errrel > critere) ;
  1196.             SI (EGA IMESS 1) ;
  1197.             Message 'Erreur relative' errrel
  1198.                     'Composante ' (EXTR LDEIN icomp)
  1199.                     ' Point Gauss' igau 'Etat' ietat ;
  1200.             Message 'eizzc.igau =' eizzc.igau 'enzzu.igau =' enzzu.igau;
  1201.             Message ;
  1202.             FINSI ;
  1203. ************ierr = 31 ; QUITTER BECART ;
  1204.          FINSI ;
  1205.       FINSI ;
  1206.    FIN BEIZZ ;
  1207.  
  1208.    icomp = 4 ;
  1209.    igau = 0 ;
  1210.    REPETER BGIXY nbgau ;
  1211.       igau = igau + 1 ;
  1212.       SI (EGA gixyc.igau 0.0) ;
  1213.          SI (NON (EGA gnxyu.igau 0.0)) ;
  1214.             SI (EGA IMESS 1) ;
  1215.             Message 'Erreur absolue' gnxyu.igau
  1216.                     'Composante ' (EXTR LDEIN icomp)
  1217.                     ' Point Gauss' igau 'Etat' ietat ;
  1218.             Message ;
  1219.             FINSI ;
  1220. ************ierr = 30 ; errabs = gnxyu.igau ; QUITTER BECART ;
  1221.          FINSI ;
  1222.       SINON ;
  1223.          dgixy = ABS (gixyc.igau - gnxyu.igau) ;
  1224.          errrel = ABS (dgixy / gixyc.igau) ;
  1225.          SI (errrel > critere) ;
  1226.             SI (EGA IMESS 1) ;
  1227.             Message 'Erreur relative' errrel
  1228.                     'Composante ' (EXTR LDEIN icomp)
  1229.                     ' Point Gauss' igau 'Etat' ietat ;
  1230.             Message 'gixyc.igau =' gixyc.igau 'gnxyu.igau =' gnxyu.igau;
  1231.             Message ;
  1232.             FINSI ;
  1233. ************ierr = 31 ; QUITTER BECART ;
  1234.          FINSI ;
  1235.       FINSI ;
  1236.    FIN BGIXY ;
  1237.  
  1238.    icomp = 5 ;
  1239.    igau = 0 ;
  1240.    REPETER BGIXZ nbgau ;
  1241.       igau = igau + 1 ;
  1242.       SI (EGA gixzc.igau 0.0) ;
  1243.          SI (NON (EGA gnxzu.igau 0.0)) ;
  1244.             SI (EGA IMESS 1) ;
  1245.             Message 'Erreur absolue' gnxzu.igau
  1246.                     'Composante ' (EXTR LDEIN icomp)
  1247.                     ' Point Gauss' igau 'Etat' ietat ;
  1248.             Message ;
  1249.             FINSI ;
  1250. ************ierr = 30 ; errabs = gnxzu.igau ; QUITTER BECART ;
  1251.          FINSI ;
  1252.       SINON ;
  1253.          dgixz = ABS (gixzc.igau - gnxzu.igau) ;
  1254.          errrel = ABS (dgixz / gixzc.igau) ;
  1255.          SI (errrel > critere) ;
  1256.             SI (EGA IMESS 1) ;
  1257.             Message 'Erreur relative' errrel
  1258.                     'Composante ' (EXTR LDEIN icomp)
  1259.                     ' Point Gauss' igau 'Etat' ietat ;
  1260.             Message 'gixzc.igau =' gixzc.igau 'gnxzu.igau =' gnxzu.igau;
  1261.             Message ;
  1262.             FINSI ;
  1263. ************ierr = 31 ; QUITTER BECART ;
  1264.          FINSI ;
  1265.       FINSI ;
  1266.    FIN BGIXZ ;
  1267.  
  1268.    icomp = 6 ;
  1269.    igau = 0 ;
  1270.    REPETER BGIYZ nbgau ;
  1271.       igau = igau + 1 ;
  1272.       SI (EGA giyzc.igau 0.0) ;
  1273.          SI (NON (EGA gnyzu.igau 0.0)) ;
  1274.             SI (EGA IMESS 1) ;
  1275.             Message 'Erreur absolue' gnyzu.igau
  1276.                     'Composante ' (EXTR LDEIN icomp)
  1277.                     ' Point Gauss' igau 'Etat' ietat ;
  1278.             Message ;
  1279.             FINSI ;
  1280. ************ierr = 30 ; errabs = gnyzu.igau ; QUITTER BECART ;
  1281.          FINSI ;
  1282.       SINON ;
  1283.          dgiyz = ABS (giyzc.igau - gnyzu.igau) ;
  1284.          errrel = ABS (dgiyz / giyzc.igau) ;
  1285.          SI (errrel > critere) ;
  1286.             SI (EGA IMESS 1) ;
  1287.             Message 'Erreur relative' errrel
  1288.                     'Composante ' (EXTR LDEIN icomp)
  1289.                     ' Point Gauss' igau 'Etat' ietat ;
  1290.             Message 'giyzc.igau =' giyzc.igau 'gnyzu.igau =' gnyzu.igau;
  1291.             Message ;
  1292.             FINSI ;
  1293. ************ierr = 31 ; QUITTER BECART ;
  1294.          FINSI ;
  1295.       FINSI ;
  1296.    FIN BGIYZ ;
  1297.  
  1298. *  Comparaison des contraintes------------------------------------------
  1299.  
  1300.    icomp = 1 ;
  1301.    igau = 0 ;
  1302.    REPETER BSMXX nbgau ;
  1303.       igau = igau + 1 ;
  1304.       SI (EGA smxxc.igau 0.0) ;
  1305.          SI (NON (EGA smxxu.igau 0.0)) ;
  1306. ************SI (EGA IMESS 1) ;
  1307. ************Message 'Erreur absolue' smxxu.igau
  1308. ************        'Composante ' (EXTR LCONT icomp)
  1309. ************        ' Point Gauss' igau 'Etat' ietat ;
  1310. ************Message ;
  1311. ************FINSI ;
  1312.             ierr = 40 ; errabs = smxxu.igau ; QUITTER BECART ;
  1313.          FINSI ;
  1314.       SINON ;
  1315.          dsmxx = ABS (smxxc.igau - smxxu.igau) ;
  1316.          errrel = ABS (dsmxx / smxxc.igau) ;
  1317.          SI (errrel > critere) ;
  1318. ************SI (EGA IMESS 1) ;
  1319. ************Message 'Erreur relative' errrel
  1320. ************        'Composante ' (EXTR LCONT icomp)
  1321. ************        ' Point Gauss' igau 'Etat' ietat ;
  1322. ************Message 'smxxc.igau =' smxxc.igau 'smxxu.igau =' smxxu.igau;
  1323. ************Message ;
  1324. ************FINSI ;
  1325.             ierr = 41 ; QUITTER BECART ;
  1326.          FINSI ;
  1327.       FINSI ;
  1328.    FIN BSMXX ;
  1329.  
  1330.    icomp = 2 ;
  1331.    igau = 0 ;
  1332.    REPETER BSMYY nbgau ;
  1333.       igau = igau + 1 ;
  1334.       SI (EGA smyyc.igau 0.0) ;
  1335.          SI (NON (EGA smyyu.igau 0.0)) ;
  1336.             SI (EGA IMESS 1) ;
  1337.             Message 'Erreur absolue' smyyu.igau
  1338.                     'Composante ' (EXTR LCONT icomp)
  1339.                     ' Point Gauss' igau 'Etat' ietat ;
  1340.             Message ;
  1341.             FINSI ;
  1342. ************ierr = 40 ; errabs = smyyu.igau ; QUITTER BECART ;
  1343.          FINSI ;
  1344.       SINON ;
  1345.          dsmyy = ABS (smyyc.igau - smyyu.igau) ;
  1346.          errrel = ABS (dsmyy / smyyc.igau) ;
  1347.          SI (errrel > critere) ;
  1348.             SI (EGA IMESS 1) ;
  1349.             Message 'Erreur relative' errrel
  1350.                     'Composante ' (EXTR LCONT icomp)
  1351.                     ' Point Gauss' igau 'Etat' ietat ;
  1352.             Message 'smyyc.igau =' smyyc.igau 'smyyu.igau =' smyyu.igau;
  1353.             Message ;
  1354.             FINSI ;
  1355. ************ierr = 41 ; QUITTER BECART ;
  1356.          FINSI ;
  1357.       FINSI ;
  1358.    FIN BSMYY ;
  1359.  
  1360.    icomp = 3 ;
  1361.    igau = 0 ;
  1362.    REPETER BSMZZ nbgau ;
  1363.       igau = igau + 1 ;
  1364.       SI (EGA smzzc.igau 0.0) ;
  1365.          SI (NON (EGA smzzu.igau 0.0)) ;
  1366.             SI (EGA IMESS 1) ;
  1367.             Message 'Erreur absolue' smzzu.igau
  1368.                     'Composante ' (EXTR LCONT icomp)
  1369.                     ' Point Gauss' igau 'Etat' ietat ;
  1370.             Message ;
  1371.             FINSI ;
  1372. ************ierr = 40 ; errabs = smzzu.igau ; QUITTER BECART ;
  1373.          FINSI ;
  1374.       SINON ;
  1375.          dsmzz = ABS (smzzc.igau - smzzu.igau) ;
  1376.          errrel = ABS (dsmzz / smzzc.igau) ;
  1377.          SI (errrel > critere) ;
  1378.             SI (EGA IMESS 1) ;
  1379.             Message 'Erreur relative' errrel
  1380.                     'Composante ' (EXTR LCONT icomp)
  1381.                     ' Point Gauss' igau 'Etat' ietat ;
  1382.             Message 'smzzc.igau =' smzzc.igau 'smzzu.igau =' smzzu.igau;
  1383.             Message ;
  1384.             FINSI ;
  1385. ************ierr = 41 ; QUITTER BECART ;
  1386.          FINSI ;
  1387.       FINSI ;
  1388.    FIN BSMZZ ;
  1389.  
  1390.    icomp = 4 ;
  1391.    igau = 0 ;
  1392.    REPETER BSMXY nbgau ;
  1393.       igau = igau + 1 ;
  1394.       SI (EGA smxyc.igau 0.0) ;
  1395.          SI (NON (EGA smxyu.igau 0.0)) ;
  1396.             SI (EGA IMESS 1) ;
  1397.             Message 'Erreur absolue' smxyu.igau
  1398.                     'Composante ' (EXTR LCONT icomp)
  1399.                     ' Point Gauss' igau 'Etat' ietat ;
  1400.             Message ;
  1401.             FINSI ;
  1402. ************ierr = 40 ; errabs = smxyu.igau ; QUITTER BECART ;
  1403.          FINSI ;
  1404.       SINON ;
  1405.          dsmxy = ABS (smxyc.igau - smxyu.igau) ;
  1406.          errrel = ABS (dsmxy / smxyc.igau) ;
  1407.          SI (errrel > critere) ;
  1408.             SI (EGA IMESS 1) ;
  1409.             Message 'Erreur relative' errrel
  1410.                     'Composante ' (EXTR LCONT icomp)
  1411.                     ' Point Gauss' igau 'Etat' ietat ;
  1412.             Message 'smxyc.igau =' smxyc.igau 'smxyu.igau =' smxyu.igau;
  1413.             Message ;
  1414.             FINSI ;
  1415. ************ierr = 41 ; QUITTER BECART ;
  1416.          FINSI ;
  1417.       FINSI ;
  1418.    FIN BSMXY ;
  1419.  
  1420.    icomp = 5 ;
  1421.    igau = 0 ;
  1422.    REPETER BSMXZ nbgau ;
  1423.       igau = igau + 1 ;
  1424.       SI (EGA smxzc.igau 0.0) ;
  1425.          SI (NON (EGA smxzu.igau 0.0)) ;
  1426.             SI (EGA IMESS 1) ;
  1427.             Message 'Erreur absolue' smxzu.igau
  1428.                     'Composante ' (EXTR LCONT icomp)
  1429.                     ' Point Gauss' igau 'Etat' ietat ;
  1430.             Message ;
  1431.             FINSI ;
  1432. ************ierr = 40 ; errabs = smxzu.igau ; QUITTER BECART ;
  1433.          FINSI ;
  1434.       SINON ;
  1435.          dsmxz = ABS (smxzc.igau - smxzu.igau) ;
  1436.          errrel = ABS (dsmxz / smxzc.igau) ;
  1437.          SI (errrel > critere) ;
  1438.             SI (EGA IMESS 1) ;
  1439.             Message 'Erreur relative' errrel
  1440.                     'Composante ' (EXTR LCONT icomp)
  1441.                     ' Point Gauss' igau 'Etat' ietat ;
  1442.             Message 'smxzc.igau =' smxzc.igau 'smxzu.igau =' smxzu.igau;
  1443.             Message ;
  1444.             FINSI ;
  1445. ************ierr = 41 ; QUITTER BECART ;
  1446.          FINSI ;
  1447.       FINSI ;
  1448.    FIN BSMXZ ;
  1449.  
  1450.    icomp = 6 ;
  1451.    igau = 0 ;
  1452.    REPETER BSMYZ nbgau ;
  1453.       igau = igau + 1 ;
  1454.       SI (EGA smyzc.igau 0.0) ;
  1455.          SI (NON (EGA smyzu.igau 0.0)) ;
  1456.             SI (EGA IMESS 1) ;
  1457.             Message 'Erreur absolue' smyzu.igau
  1458.                     'Composante ' (EXTR LCONT icomp)
  1459.                     ' Point Gauss' igau 'Etat' ietat ;
  1460.             Message ;
  1461.             FINSI ;
  1462. ************ierr = 40 ; errabs = smyzu.igau ; QUITTER BECART ;
  1463.          FINSI ;
  1464.       SINON ;
  1465.          dsmyz = ABS (smyzc.igau - smyzu.igau) ;
  1466.          errrel = ABS (dsmyz / smyzc.igau) ;
  1467.          SI (errrel > critere) ;
  1468.             SI (EGA IMESS 1) ;
  1469.             Message 'Erreur relative' errrel
  1470.                     'Composante ' (EXTR LCONT icomp)
  1471.                     ' Point Gauss' igau 'Etat' ietat ;
  1472.             Message 'smyzc.igau =' smyzc.igau 'smyzu.igau =' smyzu.igau;
  1473.             Message ;
  1474.             FINSI ;
  1475. ************ierr = 41 ; QUITTER BECART ;
  1476.          FINSI ;
  1477.       FINSI ;
  1478.    FIN BSMYZ ;
  1479.  
  1480. *  Comparaison des variables internes-----------------------------------
  1481.  
  1482.    icomp = 1 ;
  1483.    igau = 0 ;
  1484.    REPETER BPXX nbgau ;
  1485.       igau = igau + 1 ;
  1486.       SI (EGA pxxc.igau 0.0) ;
  1487.          SI (NON (EGA pxxu.igau 0.0)) ;
  1488.             SI (EGA IMESS 1) ;
  1489.             Message 'Erreur absolue' pxxu.igau
  1490.                     'Composante ' (EXTR LVARI icomp)
  1491.                     ' Point Gauss' igau 'Etat' ietat ;
  1492.             Message ;
  1493.             FINSI ;
  1494. ************ierr = 50 ; errabs = pxxu.igau ; QUITTER BECART ;
  1495.          FINSI ;
  1496.       SINON ;
  1497.          dpxx = ABS (pxxc.igau - pxxu.igau) ;
  1498.          errrel = ABS (dpxx / pxxc.igau) ;
  1499.          SI (errrel > critere) ;
  1500.             SI (EGA IMESS 1) ;
  1501.             Message 'Erreur relative' errrel
  1502.                     'Composante ' (EXTR LVARI icomp)
  1503.                     ' Point Gauss' igau 'Etat' ietat ;
  1504.             Message 'pxxc.igau =' pxxc.igau 'pxxu.igau =' pxxu.igau;
  1505.             Message ;
  1506.             FINSI ;
  1507. ************ierr = 51 ; QUITTER BECART ;
  1508.          FINSI ;
  1509.       FINSI ;
  1510.    FIN BPXX ;
  1511.  
  1512.    icomp = 2 ;
  1513.    igau = 0 ;
  1514.    REPETER BPYY nbgau ;
  1515.       igau = igau + 1 ;
  1516.       SI (EGA pyyc.igau 0.0) ;
  1517.          SI (NON (EGA pyyu.igau 0.0)) ;
  1518.             SI (EGA IMESS 1) ;
  1519.             Message 'Erreur absolue' pyyu.igau
  1520.                     'Composante ' (EXTR LVARI icomp)
  1521.                     ' Point Gauss' igau 'Etat' ietat ;
  1522.             Message ;
  1523.             FINSI ;
  1524. ************ierr = 50 ; errabs = pyyu.igau ; QUITTER BECART ;
  1525.          FINSI ;
  1526.       SINON ;
  1527.          dpyy = ABS (pyyc.igau - pyyu.igau) ;
  1528.          errrel = ABS (dpyy / pyyc.igau) ;
  1529.          SI (errrel > critere) ;
  1530.             SI (EGA IMESS 1) ;
  1531.             Message 'Erreur relative' errrel
  1532.                     'Composante ' (EXTR LVARI icomp)
  1533.                     ' Point Gauss' igau 'Etat' ietat ;
  1534.             Message 'pyyc.igau =' pyyc.igau 'pyyu.igau =' pyyu.igau;
  1535.             Message ;
  1536.             FINSI ;
  1537. ************ierr = 51 ; QUITTER BECART ;
  1538.          FINSI ;
  1539.       FINSI ;
  1540.    FIN BPYY ;
  1541.  
  1542.    icomp = 3 ;
  1543.    igau = 0 ;
  1544.    REPETER BPZZ nbgau ;
  1545.       igau = igau + 1 ;
  1546.       SI (EGA pzzc.igau 0.0) ;
  1547.          SI (NON (EGA pzzu.igau 0.0)) ;
  1548.             SI (EGA IMESS 1) ;
  1549.             Message 'Erreur absolue' pzzu.igau
  1550.                     'Composante ' (EXTR LVARI icomp)
  1551.                     ' Point Gauss' igau 'Etat' ietat ;
  1552.             Message ;
  1553.             FINSI ;
  1554. ************ierr = 50 ; errabs = pzzu.igau ; QUITTER BECART ;
  1555.          FINSI ;
  1556.       SINON ;
  1557.          dpzz = ABS (pzzc.igau - pzzu.igau) ;
  1558.          errrel = ABS (dpzz / pzzc.igau) ;
  1559.          SI (errrel > critere) ;
  1560.             SI (EGA IMESS 1) ;
  1561.             Message 'Erreur relative' errrel
  1562.                     'Composante ' (EXTR LVARI icomp)
  1563.                     ' Point Gauss' igau 'Etat' ietat ;
  1564.             Message 'pzzc.igau =' pzzc.igau 'pzzu.igau =' pzzu.igau;
  1565.             Message ;
  1566.             FINSI ;
  1567. ************ierr = 51 ; QUITTER BECART ;
  1568.          FINSI ;
  1569.       FINSI ;
  1570.    FIN BPZZ ;
  1571.  
  1572.    icomp = 4 ;
  1573.    igau = 0 ;
  1574.    REPETER BPXY nbgau ;
  1575.       igau = igau + 1 ;
  1576.       SI (EGA pxyc.igau 0.0) ;
  1577.          SI (NON (EGA pxyu.igau 0.0)) ;
  1578.             SI (EGA IMESS 1) ;
  1579.             Message 'Erreur absolue' pxyu.igau
  1580.                     'Composante ' (EXTR LVARI icomp)
  1581.                     ' Point Gauss' igau 'Etat' ietat ;
  1582.             Message ;
  1583.             FINSI ;
  1584. ************ierr = 50 ; errabs = pxyu.igau ; QUITTER BECART ;
  1585.          FINSI ;
  1586.       SINON ;
  1587.          dpxy = ABS (pxyc.igau - pxyu.igau) ;
  1588.          errrel = ABS (dpxy / pxyc.igau) ;
  1589.          SI (errrel > critere) ;
  1590.             SI (EGA IMESS 1) ;
  1591.             Message 'Erreur relative' errrel
  1592.                     'Composante ' (EXTR LVARI icomp)
  1593.                     ' Point Gauss' igau 'Etat' ietat ;
  1594.             Message 'pxyc.igau =' pxyc.igau 'pxyu.igau =' pxyu.igau;
  1595.             Message ;
  1596.             FINSI ;
  1597. ************ierr = 51 ; QUITTER BECART ;
  1598.          FINSI ;
  1599.       FINSI ;
  1600.    FIN BPXY ;
  1601.  
  1602.    icomp = 5 ;
  1603.    igau = 0 ;
  1604.    REPETER BPYZ nbgau ;
  1605.       igau = igau + 1 ;
  1606.       SI (EGA pyzc.igau 0.0) ;
  1607.          SI (NON (EGA pyzu.igau 0.0)) ;
  1608.             SI (EGA IMESS 1) ;
  1609.             Message 'Erreur absolue' pyzu.igau
  1610.                     'Composante ' (EXTR LVARI icomp)
  1611.                     ' Point Gauss' igau 'Etat' ietat ;
  1612.             Message ;
  1613.             FINSI ;
  1614. ************ierr = 50 ; errabs = pyzu.igau ; QUITTER BECART ;
  1615.          FINSI ;
  1616.       SINON ;
  1617.          dpyz = ABS (pyzc.igau - pyzu.igau) ;
  1618.          errrel = ABS (dpyz / pyzc.igau) ;
  1619.          SI (errrel > critere) ;
  1620.             SI (EGA IMESS 1) ;
  1621.             Message 'Erreur relative' errrel
  1622.                     'Composante ' (EXTR LVARI icomp)
  1623.                     ' Point Gauss' igau 'Etat' ietat ;
  1624.             Message 'pyzc.igau =' pyzc.igau 'pyzu.igau =' pyzu.igau;
  1625.             Message ;
  1626.             FINSI ;
  1627. ************ierr = 51 ; QUITTER BECART ;
  1628.          FINSI ;
  1629.       FINSI ;
  1630.    FIN BPYZ ;
  1631.  
  1632.    icomp = 6 ;
  1633.    igau = 0 ;
  1634.    REPETER BPZX nbgau ;
  1635.       igau = igau + 1 ;
  1636.       SI (EGA pzxc.igau 0.0) ;
  1637.          SI (NON (EGA pzxu.igau 0.0)) ;
  1638.             SI (EGA IMESS 1) ;
  1639.             Message 'Erreur absolue' pzxu.igau
  1640.                     'Composante ' (EXTR LVARI icomp)
  1641.                     ' Point Gauss' igau 'Etat' ietat ;
  1642.             Message ;
  1643.             FINSI ;
  1644. ************ierr = 50 ; errabs = pzxu.igau ; QUITTER BECART ;
  1645.          FINSI ;
  1646.       SINON ;
  1647.          dpzx = ABS (pzxc.igau - pzxu.igau) ;
  1648.          errrel = ABS (dpzx / pzxc.igau) ;
  1649.          SI (errrel > critere) ;
  1650.             SI (EGA IMESS 1) ;
  1651.             Message 'Erreur relative' errrel
  1652.                     'Composante ' (EXTR LVARI icomp)
  1653.                     ' Point Gauss' igau 'Etat' ietat ;
  1654.             Message 'pzxc.igau =' pzxc.igau 'pzxu.igau =' pzxu.igau;
  1655.             Message ;
  1656.             FINSI ;
  1657. ************ierr = 51 ; QUITTER BECART ;
  1658.          FINSI ;
  1659.       FINSI ;
  1660.    FIN BPZX ;
  1661.  
  1662.    icomp = 7 ;
  1663.    igau = 0 ;
  1664.    REPETER BVNXX nbgau ;
  1665.       igau = igau + 1 ;
  1666.       SI (EGA vnxxc.igau 0.0) ;
  1667.          SI (NON (EGA vnxxu.igau 0.0)) ;
  1668.             SI (EGA IMESS 1) ;
  1669.             Message 'Erreur absolue' vnxxu.igau
  1670.                     'Composante ' (EXTR LVARI icomp)
  1671.                     ' Point Gauss' igau 'Etat' ietat ;
  1672.             Message ;
  1673.             FINSI ;
  1674. ************ierr = 50 ; errabs = vnxxu.igau ; QUITTER BECART ;
  1675.          FINSI ;
  1676.       SINON ;
  1677.          dvnxx = ABS (vnxxc.igau - vnxxu.igau) ;
  1678.          errrel = ABS (dvnxx / vnxxc.igau) ;
  1679.          SI (errrel > critere) ;
  1680.             SI (EGA IMESS 1) ;
  1681.             Message 'Erreur relative' errrel
  1682.                     'Composante ' (EXTR LVARI icomp)
  1683.                     ' Point Gauss' igau 'Etat' ietat ;
  1684.             Message 'vnxxc.igau =' vnxxc.igau 'vnxxu.igau =' vnxxu.igau;
  1685.             Message ;
  1686.             FINSI ;
  1687. ************ierr = 51 ; QUITTER BECART ;
  1688.          FINSI ;
  1689.       FINSI ;
  1690.    FIN BVNXX ;
  1691.  
  1692.    icomp = 8 ;
  1693.    igau = 0 ;
  1694.    REPETER BVNYY nbgau ;
  1695.       igau = igau + 1 ;
  1696.       SI (EGA vnyyc.igau 0.0) ;
  1697.          SI (NON (EGA vnyyu.igau 0.0)) ;
  1698.             SI (EGA IMESS 1) ;
  1699.             Message 'Erreur absolue' vnyyu.igau
  1700.                     'Composante ' (EXTR LVARI icomp)
  1701.                     ' Point Gauss' igau 'Etat' ietat ;
  1702.             Message ;
  1703.             FINSI ;
  1704. ************ierr = 50 ; errabs = vnyyu.igau ; QUITTER BECART ;
  1705.          FINSI ;
  1706.       SINON ;
  1707.          dvnyy = ABS (vnyyc.igau - vnyyu.igau) ;
  1708.          errrel = ABS (dvnyy / vnyyc.igau) ;
  1709.          SI (errrel > critere) ;
  1710.             SI (EGA IMESS 1) ;
  1711.             Message 'Erreur relative' errrel
  1712.                     'Composante ' (EXTR LVARI icomp)
  1713.                     ' Point Gauss' igau 'Etat' ietat ;
  1714.             Message 'vnyyc.igau =' vnyyc.igau 'vnyyu.igau =' vnyyu.igau;
  1715.             Message ;
  1716.             FINSI ;
  1717. ************ierr = 51 ; QUITTER BECART ;
  1718.          FINSI ;
  1719.       FINSI ;
  1720.    FIN BVNYY ;
  1721.  
  1722.    icomp = 9 ;
  1723.    igau = 0 ;
  1724.    REPETER BVNZZ nbgau ;
  1725.       igau = igau + 1 ;
  1726.       SI (EGA vnzzc.igau 0.0) ;
  1727.          SI (NON (EGA vnzzu.igau 0.0)) ;
  1728.             SI (EGA IMESS 1) ;
  1729.             Message 'Erreur absolue' vnzzu.igau
  1730.                     'Composante ' (EXTR LVARI icomp)
  1731.                     ' Point Gauss' igau 'Etat' ietat ;
  1732.             Message ;
  1733.             FINSI ;
  1734. ************ierr = 50 ; errabs = vnzzu.igau ; QUITTER BECART ;
  1735.          FINSI ;
  1736.       SINON ;
  1737.          dvnzz = ABS (vnzzc.igau - vnzzu.igau) ;
  1738.          errrel = ABS (dvnzz / vnzzc.igau) ;
  1739.          SI (errrel > critere) ;
  1740.             SI (EGA IMESS 1) ;
  1741.             Message 'Erreur relative' errrel
  1742.                     'Composante ' (EXTR LVARI icomp)
  1743.                     ' Point Gauss' igau 'Etat' ietat ;
  1744.             Message 'vnzzc.igau =' vnzzc.igau 'vnzzu.igau =' vnzzu.igau;
  1745.             Message ;
  1746.             FINSI ;
  1747. ************ierr = 51 ; QUITTER BECART ;
  1748.          FINSI ;
  1749.       FINSI ;
  1750.    FIN BVNZZ ;
  1751.  
  1752.    icomp = 10 ;
  1753.    igau = 0 ;
  1754.    REPETER BVNXY nbgau ;
  1755.       igau = igau + 1 ;
  1756.       SI (EGA vnxyc.igau 0.0) ;
  1757.          SI (NON (EGA vnxyu.igau 0.0)) ;
  1758.             SI (EGA IMESS 1) ;
  1759.             Message 'Erreur absolue' vnxyu.igau
  1760.                     'Composante ' (EXTR LVARI icomp)
  1761.                     ' Point Gauss' igau 'Etat' ietat ;
  1762.             Message ;
  1763.             FINSI ;
  1764. ************ierr = 50 ; errabs = vnxyu.igau ; QUITTER BECART ;
  1765.          FINSI ;
  1766.       SINON ;
  1767.          dvnxy = ABS (vnxyc.igau - vnxyu.igau) ;
  1768.          errrel = ABS (dvnxy / vnxyc.igau) ;
  1769.          SI (errrel > critere) ;
  1770.             SI (EGA IMESS 1) ;
  1771.             Message 'Erreur relative' errrel
  1772.                     'Composante ' (EXTR LVARI icomp)
  1773.                     ' Point Gauss' igau 'Etat' ietat ;
  1774.             Message 'vnxyc.igau =' vnxyc.igau 'vnxyu.igau =' vnxyu.igau;
  1775.             Message ;
  1776.             FINSI ;
  1777. ************ierr = 51 ; QUITTER BECART ;
  1778.          FINSI ;
  1779.       FINSI ;
  1780.    FIN BVNXY ;
  1781.  
  1782.    icomp = 11 ;
  1783.    igau = 0 ;
  1784.    REPETER BVNYZ nbgau ;
  1785.       igau = igau + 1 ;
  1786.       SI (EGA vnyzc.igau 0.0) ;
  1787.          SI (NON (EGA vnyzu.igau 0.0)) ;
  1788.             SI (EGA IMESS 1) ;
  1789.             Message 'Erreur absolue' vnyzu.igau
  1790.                     'Composante ' (EXTR LVARI icomp)
  1791.                     ' Point Gauss' igau 'Etat' ietat ;
  1792.             Message ;
  1793.             FINSI ;
  1794. ************ierr = 50 ; errabs = vnyzu.igau ; QUITTER BECART ;
  1795.          FINSI ;
  1796.       SINON ;
  1797.          dvnyz = ABS (vnyzc.igau - vnyzu.igau) ;
  1798.          errrel = ABS (dvnyz / vnyzc.igau) ;
  1799.          SI (errrel > critere) ;
  1800.             SI (EGA IMESS 1) ;
  1801.             Message 'Erreur relative' errrel
  1802.                     'Composante ' (EXTR LVARI icomp)
  1803.                     ' Point Gauss' igau 'Etat' ietat ;
  1804.             Message 'vnyzc.igau =' vnyzc.igau 'vnyzu.igau =' vnyzu.igau;
  1805.             Message ;
  1806.             FINSI ;
  1807. ************ierr = 51 ; QUITTER BECART ;
  1808.          FINSI ;
  1809.       FINSI ;
  1810.    FIN BVNYZ ;
  1811.  
  1812.    icomp = 12 ;
  1813.    igau = 0 ;
  1814.    REPETER BVNZX nbgau ;
  1815.       igau = igau + 1 ;
  1816.       SI (EGA vnzxc.igau 0.0) ;
  1817.          SI (NON (EGA vnzxu.igau 0.0)) ;
  1818.             SI (EGA IMESS 1) ;
  1819.             Message 'Erreur absolue' vnzxu.igau
  1820.                     'Composante ' (EXTR LVARI icomp)
  1821.                     ' Point Gauss' igau 'Etat' ietat ;
  1822.             Message ;
  1823.             FINSI ;
  1824. ************ierr = 50 ; errabs = vnzxu.igau ; QUITTER BECART ;
  1825.          FINSI ;
  1826.       SINON ;
  1827.          dvnzx = ABS (vnzxc.igau - vnzxu.igau) ;
  1828.          errrel = ABS (dvnzx / vnzxc.igau) ;
  1829.          SI (errrel > critere) ;
  1830.             SI (EGA IMESS 1) ;
  1831.             Message 'Erreur relative' errrel
  1832.                     'Composante ' (EXTR LVARI icomp)
  1833.                     ' Point Gauss' igau 'Etat' ietat ;
  1834.             Message 'vnzxc.igau =' vnzxc.igau 'vnzxu.igau =' vnzxu.igau;
  1835.             Message ;
  1836.             FINSI ;
  1837. ************ierr = 51 ; QUITTER BECART ;
  1838.          FINSI ;
  1839.       FINSI ;
  1840.    FIN BVNZX ;
  1841.  
  1842.    icomp = 13 ;
  1843.    igau = 0 ;
  1844.    REPETER BEPCE nbgau ;
  1845.       igau = igau + 1 ;
  1846.       SI (EGA epcec.igau 0.0) ;
  1847.          SI (NON (EGA epceu.igau 0.0)) ;
  1848.             SI (EGA IMESS 1) ;
  1849.             Message 'Erreur absolue' epceu.igau
  1850.                     'Composante ' (EXTR LVARI icomp)
  1851.                     ' Point Gauss' igau 'Etat' ietat ;
  1852.             Message ;
  1853.             FINSI ;
  1854. ************ierr = 50 ; errabs = epceu.igau ; QUITTER BECART ;
  1855.          FINSI ;
  1856.       SINON ;
  1857.          depce = ABS (epcec.igau - epceu.igau) ;
  1858.          errrel = ABS (depce / epcec.igau) ;
  1859.          SI (errrel > critere) ;
  1860.             SI (EGA IMESS 1) ;
  1861.             Message 'Erreur relative' errrel
  1862.                     'Composante ' (EXTR LVARI icomp)
  1863.                     ' Point Gauss' igau 'Etat' ietat ;
  1864.             Message 'epcec.igau =' epcec.igau 'epceu.igau =' epceu.igau;
  1865.             Message ;
  1866.             FINSI ;
  1867. ************ierr = 51 ; QUITTER BECART ;
  1868.          FINSI ;
  1869.       FINSI ;
  1870.    FIN BEPCE ;
  1871.  
  1872.    icomp = 14 ;
  1873.    igau = 0 ;
  1874.    REPETER BEPSH nbgau ;
  1875.       igau = igau + 1 ;
  1876.       SI (EGA epshc.igau 0.0) ;
  1877.          SI (NON (EGA epshu.igau 0.0)) ;
  1878.             SI (EGA IMESS 1) ;
  1879.             Message 'Erreur absolue' epshu.igau
  1880.                     'Composante ' (EXTR LVARI icomp)
  1881.                     ' Point Gauss' igau 'Etat' ietat ;
  1882.             Message ;
  1883.             FINSI ;
  1884. ************ierr = 50 ; errabs = epshu.igau ; QUITTER BECART ;
  1885.          FINSI ;
  1886.       SINON ;
  1887.          depsh = ABS (epshc.igau - epshu.igau) ;
  1888.          errrel = ABS (depsh / epshc.igau) ;
  1889.          SI (errrel > critere) ;
  1890.             SI (EGA IMESS 1) ;
  1891.             Message 'Erreur relative' errrel
  1892.                     'Composante ' (EXTR LVARI icomp)
  1893.                     ' Point Gauss' igau 'Etat' ietat ;
  1894.             Message 'epshc.igau =' epshc.igau 'epshu.igau =' epshu.igau;
  1895.             Message ;
  1896.             FINSI ;
  1897. ************ierr = 51 ; QUITTER BECART ;
  1898.          FINSI ;
  1899.       FINSI ;
  1900.    FIN BEPSH ;
  1901.  
  1902.    icomp = 15 ;
  1903.    igau = 0 ;
  1904.    REPETER BLL nbgau ;
  1905.       igau = igau + 1 ;
  1906.       SI (EGA llc.igau 0.0) ;
  1907.          SI (NON (EGA llu.igau 0.0)) ;
  1908.             SI (EGA IMESS 1) ;
  1909.             Message 'Erreur absolue' llu.igau
  1910.                     'Composante ' (EXTR LVARI icomp)
  1911.                     ' Point Gauss' igau 'Etat' ietat ;
  1912.             Message ;
  1913.             FINSI ;
  1914. ************ierr = 50 ; errabs = llu.igau ; QUITTER BECART ;
  1915.          FINSI ;
  1916.       SINON ;
  1917.          dll = ABS (llc.igau - llu.igau) ;
  1918.          errrel = ABS (dll / llc.igau) ;
  1919.          SI (errrel > critere) ;
  1920.             SI (EGA IMESS 1) ;
  1921.             Message 'Erreur relative' errrel
  1922.                     'Composante ' (EXTR LVARI icomp)
  1923.                     ' Point Gauss' igau 'Etat' ietat ;
  1924.             Message 'llc.igau =' llc.igau 'llu.igau =' llu.igau;
  1925.             Message ;
  1926.             FINSI ;
  1927. ************ierr = 51 ; QUITTER BECART ;
  1928.          FINSI ;
  1929.       FINSI ;
  1930.    FIN BLL ;
  1931.  
  1932.    icomp = 16 ;
  1933.    igau = 0 ;
  1934.    REPETER BEPSE nbgau ;
  1935.       igau = igau + 1 ;
  1936.       SI (EGA epsec.igau 0.0) ;
  1937.          SI (NON (EGA epseu.igau 0.0)) ;
  1938.             SI (EGA IMESS 1) ;
  1939.             Message 'Erreur absolue' epseu.igau
  1940.                     'Composante ' (EXTR LVARI icomp)
  1941.                     ' Point Gauss' igau 'Etat' ietat ;
  1942.             Message ;
  1943.             FINSI ;
  1944. ************ierr = 50 ; errabs = epseu.igau ; QUITTER BECART ;
  1945.          FINSI ;
  1946.       SINON ;
  1947.          depse = ABS (epsec.igau - epseu.igau) ;
  1948.          errrel = ABS (depse / epsec.igau) ;
  1949.          SI (errrel > critere) ;
  1950.             SI (EGA IMESS 1) ;
  1951.             Message 'Erreur relative' errrel
  1952.                     'Composante ' (EXTR LVARI icomp)
  1953.                     ' Point Gauss' igau 'Etat' ietat ;
  1954.             Message 'epsec.igau =' epsec.igau 'epseu.igau =' epseu.igau;
  1955.             Message ;
  1956.             FINSI ;
  1957. ************ierr = 51 ; QUITTER BECART ;
  1958.          FINSI ;
  1959.       FINSI ;
  1960.    FIN BEPSE ;
  1961.  
  1962. FIN BECART ;
  1963.  
  1964. * FIN   BOUCLE COMPARAISON DES RESULTATS________________________________
  1965.  
  1966. * MESSAGES D'ERREUR_____________________________________________________
  1967.  
  1968. Message ;
  1969. Message 'CPU avec modele interne = ' CPUint 'centiemes de seconde' ;
  1970. Message 'CPU avec modele externe = ' CPUext 'centiemes de seconde' ;
  1971.  
  1972. ifin = netat-1 ;
  1973. eixxcmax = 0.0 ;
  1974. igau = 0 ;
  1975. REPETER BEIXXFC nbgau ;
  1976.    igau = igau + 1 ;
  1977.    eixx = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ifin)
  1978.                     'EIXX' 1 1 igau ;
  1979.    SI (eixx > eixxcmax) ; eixxcmax = eixx ; FINSI ;
  1980. FIN BEIXXFC ;
  1981. eixxumax = 0.0 ;
  1982. igau = 0 ;
  1983. REPETER BEIXXFU nbgau ;
  1984.    igau = igau + 1 ;
  1985.    eixx = EXTR (tabutil.'VARIABLES_INTERNES'.ifin)
  1986.                     'ENXX' 1 1 igau ;
  1987.    SI (eixx > eixxumax) ; eixxumax = eixx ; FINSI ;
  1988. FIN BEIXXFU ;
  1989. Message ;
  1990. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  1991. Message 'EIXX max avec modele interne = ' eixxcmax ;
  1992. Message 'EIXX max avec modele externe = ' eixxumax ;
  1993. Message ;
  1994.  
  1995. SI (ierr > 0) ;
  1996.  
  1997.    SI (EGA ierr 1) ;
  1998.       Message 'Instants de calcul non coincidants' ;
  1999.       ERRE 5 ;
  2000.    FINSI ;
  2001.  
  2002.    SI (EGA ierr 10) ;
  2003.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  2004.               ' Noeud' ino 'Etat' ietat ;
  2005.       ERRE 5 ;
  2006.    FINSI ;
  2007.  
  2008.    SI (EGA ierr 11) ;
  2009.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  2010.               ' Noeud' ino 'Etat' ietat ;
  2011.       ERRE 5 ;
  2012.    FINSI ;
  2013.  
  2014.    SI (EGA ierr 20) ;
  2015.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  2016.               ' Point Gauss' igau 'Etat' ietat ;
  2017.       ERRE 5 ;
  2018.    FINSI ;
  2019.  
  2020.    SI (EGA ierr 21) ;
  2021.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  2022.               ' Point Gauss' igau 'Etat' ietat ;
  2023.       ERRE 5 ;
  2024.    FINSI ;
  2025.  
  2026.    SI (EGA ierr 30) ;
  2027.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  2028.               ' Point Gauss' igau 'Etat' ietat ;
  2029.       ERRE 5 ;
  2030.    FINSI ;
  2031.  
  2032.    SI (EGA ierr 31) ;
  2033.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  2034.               ' Point Gauss' igau 'Etat' ietat ;
  2035.       ERRE 5 ;
  2036.    FINSI ;
  2037.  
  2038.    SI (EGA ierr 40) ;
  2039.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  2040.               ' Point Gauss' igau 'Etat' ietat ;
  2041.       ERRE 5 ;
  2042.    FINSI ;
  2043.  
  2044.    SI (EGA ierr 41) ;
  2045.       Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  2046.               ' Point Gauss' igau 'Etat' ietat ;
  2047.       ERRE 5 ;
  2048.    FINSI ;
  2049.  
  2050.    SI (EGA ierr 50) ;
  2051.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  2052.               ' Point Gauss' igau 'Etat' ietat ;
  2053.       ERRE 5 ;
  2054.    FINSI ;
  2055.  
  2056.    SI (EGA ierr 51) ;
  2057.       Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  2058.               ' Point Gauss' igau 'Etat' ietat ;
  2059.       ERRE 5 ;
  2060.    FINSI ;
  2061.  
  2062. SINON ;
  2063.  
  2064.    Message 'Resultats concordants entre les deux modeles' ;
  2065.    Message 'Aucune erreur relative superieure a ' critere
  2066.            ' pour les composantes' ;
  2067.    Message (EXTR LDEFO idefo) '  ' (EXTR LDEIN idein)
  2068.                               '  ' (EXTR LCONT icont) ;
  2069.    ERRE 0 ;
  2070.  
  2071. FINSI ;
  2072.  
  2073. FIN ;
  2074.  
  2075.  
  2076.  
  2077.  
  2078.  
  2079.  
  2080.  
  2081.  
  2082.  
  2083.  

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