Télécharger la plaquette Cast3M


fr en

Télécharger creep03_cisXY.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  creep03_cisXY.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 'VISCO_EXTERNE' 'GENERAL', integre par CCREEP (schema CCONST)
  13. *
  14. * Test No3 : Modele de FLUAGE POLYNOMIAL
  15. *            CISAILLEMENT XY
  16. *            Modele externalise non instanciable (ILOI=21) compare
  17. *            au modele externalise non instanciable avec variables
  18. *            internes supplementaires (ILOI=23)
  19. *
  20. *=======================================================================
  21.  
  22.  
  23. opti echo 1 ;
  24.  
  25.  
  26. *=======================================================================
  27. *
  28. *                       PROCEDURES DEPENDANTES
  29. *
  30. *=======================================================================
  31.  
  32.  
  33.  
  34.  
  35. *=======================================================================
  36. *
  37. *                      JEU DE DONNEES PRINCIPAL
  38. *
  39. *=======================================================================
  40.  
  41. *-----------------------------------------------------------------------
  42. * 0 - Maillage
  43. *     PLAQUE CARREE    COTE : 1 m    EPAISSEUR : 1 cm
  44.  
  45. opti dime 3 elem cub8 mode tridim ;
  46.  
  47. l1 = 1.0 ; 
  48. l2 = 1.0 ; 
  49. ep0 = 0.01 ; 
  50.  
  51. p01 = 0.0 0.0 0.0 ; 
  52. p02 = ep0 0.0 0.0 ; 
  53. p03 = ep0 l2  0.0 ; 
  54. p04 = 0.0 l2  0.0 ; 
  55.  
  56. p10 = 0.0 0.0 l1 ;
  57. p20 = ep0 0.0 l1 ;
  58. p30 = ep0 l2  l1 ;
  59. p40 = 0.0 l2  l1 ; 
  60.  
  61. li01 = 'DROIT' 1 p01 p02 ; 
  62. li02 = 'DROIT' 1 p02 p03 ;
  63. li03 = 'DROIT' 1 p03 p04 ;
  64. li04 = 'DROIT' 1 p04 p01 ;
  65.  
  66. li10 = 'DROIT' 1 p10 p20 ;
  67. li20 = 'DROIT' 1 p20 p30 ;
  68. li30 = 'DROIT' 1 p30 p40 ;
  69. li40 = 'DROIT' 1 p40 p10 ;
  70.  
  71. sur01 = 'DALL' li01 li02 li03 li04 ;
  72. sur10 = 'DALL' li10 li20 li30 li40 ;
  73. mail1 = sur01 'VOLU' 1 sur10 ;
  74.  
  75. li2  = 'DROIT' 1 p02 p20 ;
  76. li3  = 'DROIT' 1 p20 p10 ;
  77. li4  = 'DROIT' 1 p10 p01 ;
  78. sbloq = 'DALL' li01 li2 li3 li4 ;
  79.  
  80. li5  = 'DROIT' 1 p04 p03 ;
  81. li6  = 'DROIT' 1 p03 p30 ;
  82. li7  = 'DROIT' 1 p30 p40 ;
  83. li8  = 'DROIT' 1 p40 p04 ;
  84. schar = 'DALL' li5 li6 li7 li8 ;
  85.  
  86. ELIM (mail1 ET sbloq ET schar) 1.e-5 ;
  87.  
  88. *-----------------------------------------------------------------------
  89. * 1 - Calcul avec modele 'VISCO_EXTERNE' 'GENERAL' no21
  90. *     Modele 'FLUAGE' 'POLYNOMIAL' externalise dans CREEP
  91. *     Modele non instanciable au niveau des iterations internes
  92.  
  93. * 1.1 Definition du modele utilisateur externalise dans CREEP
  94.  
  95. mo_uti1 = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  96.                      'VISCO_EXTERNE' 'GENERAL' 'NUME_LOI' 21 ;
  97.  
  98. * 1.2 Affectation des proprietes de materiau
  99.  
  100. densiu   = 0.95 ;
  101. dgrain   = 8.0E-6 ;
  102. tuo2     = 1500.0 ;
  103. tauxfisk = 1.5E19 ;
  104.  
  105. xyoun = (2.2693E11 - (1.5399E7*tuo2)) - (9.35973E3 * (tuo2**2.0)) ;
  106. xcis  = (8.583E10 - (5.157E6*tuo2)) - (3.747E3 * (tuo2**2.0)) ;
  107. xnu   = (xyoun/(2.0*xcis)) - 1.0 ;
  108. xrho  = 10950.0 ;
  109. xalph = ( (1.1833E-5 - (5.013E-9*tuo2)) + (3.756E-12 * (tuo2**2.0)) )
  110.         - (6.125E-17 * (tuo2**3.0)) ;
  111.  
  112. ma_uti1 = MATE mo_uti1 'YOUN' xyoun 'NU  ' xnu   'RHO ' xrho
  113.                        'ALPH' xalph ;
  114.  
  115. * 1.3 Definition des conditions aux limites et chargements
  116.  
  117. cl1 = BLOQ mail1 'UY' 'UZ' ;
  118. cl2 = BLOQ sbloq 'UX' ;
  119. cl0 = cl1 ET cl2 ;
  120.  
  121. depmax = 5.0E8 ;
  122. dep0 = FORCE schar 'FX' depmax 'FY' 0.0 'FZ' 0.0 ;
  123. ev0 = EVOL MANU temps (PROG 0.0 1.E-2 1000.0)
  124.                     y (PROG 0.0 1.0   1.0   ) ;
  125. cha0 = CHAR MECA dep0 ev0 ;
  126.  
  127. * 1.4 Definition des instants de calcul
  128.  
  129. pro0 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  130. pro1 = PROG 0.0 PAS 1.E-6 3.0E-5 ;
  131.  
  132. * 1.5 Resolution par PASAPAS
  133.  
  134. tabuti1 = TABLE ;
  135. tabuti1.'MODELE' = mo_uti1 ;
  136. tabuti1.'CARACTERISTIQUES' = ma_uti1 ;
  137. tabuti1.'CHARGEMENT' = cha0 ;
  138. tabuti1.'BLOCAGES_MECANIQUES' = cl0 ;
  139. tabuti1.'TEMPS_CALCULES' = pro0 ;
  140. tabuti1.'TEMPS_SAUVES' = pro1 ;
  141. TMASAU=table;
  142. tabuti1 . 'MES_SAUVEGARDES'=TMASAU;
  143. TMASAU .'DEFTO'=VRAI;
  144. TMASAU .'DEFIN'=VRAI;
  145. TEMPS 'ZERO' ;
  146. PASAPAS tabuti1 ;
  147. TABTPS = TEMP 'NOEC';
  148. CPUuti1 = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  149.  
  150. *-----------------------------------------------------------------------
  151. * 2 - Calcul avec modele 'VISCO_EXTERNE' 'GENERAL' no23
  152. *     Modele 'FLUAGE' 'POLYNOMIAL' externalise dans CREEP
  153. *     Modele non instanciable au niveau des iterations internes,
  154. *     comprenant deux variables internes supplementaires
  155.  
  156. * 2.1 Definition du modele utilisateur externalise dans CREEP
  157.  
  158. LCVAR23 = MOTS 'PSUP' 'QSUP' ;
  159.  
  160. mo_uti2 = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  161.                      'VISCO_EXTERNE' 'GENERAL'
  162.                      'NUME_LOI' 23 'C_VARINTER' LCVAR23 ;
  163.  
  164. * 2.2 Affectation des proprietes de materiau
  165.  
  166. ma_uti2 = MATE mo_uti2 'YOUN' xyoun 'NU  ' xnu   'RHO ' xrho
  167.                        'ALPH' xalph ;
  168.  
  169. * 2.3 Definition des conditions aux limites et chargements : idem 1.3
  170.  
  171. * 2.4 Definition des instants de calcul : idem 1.4
  172.  
  173. * 2.5 Resolution par PASAPAS
  174.  
  175. tabuti2 = TABLE ;
  176. tabuti2.'MODELE' = mo_uti2 ;
  177. tabuti2.'CARACTERISTIQUES' = ma_uti2 ;
  178. tabuti2.'CHARGEMENT' = cha0 ;
  179. tabuti2.'BLOCAGES_MECANIQUES' = cl0 ;
  180. tabuti2.'TEMPS_CALCULES' = pro0 ;
  181. tabuti2.'TEMPS_SAUVES' = pro1 ;
  182. TMASAU=table;
  183. tabuti2 . 'MES_SAUVEGARDES'=TMASAU;
  184. TMASAU .'DEFTO'=VRAI;
  185. TMASAU .'DEFIN'=VRAI;
  186. TEMPS 'ZERO' ;
  187. PASAPAS tabuti2 ;
  188. TABTPS = TEMP 'NOEC';
  189. CPUuti2 = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  190.  
  191. *-----------------------------------------------------------------------
  192. * 3 - Post-traitement : comparaisons entre les deux modeles
  193. *
  194. *     N.B. Comme on simule un essai de cisaillement XY :
  195. *          - Critere strict (errrel < 1.d-8) :
  196. *            * sur GAXY, GIXY et SMXY,
  197. *            * sur les variables internes EC0, ESW0, P, QTILD,
  198. *            * pour le modele externe no23 : sur les variables internes
  199. *              supplementaires PSUP et QSUP (comparees a P et QTILD).
  200. *          - Pour les autres composantes on note les depassements du
  201. *            critere si impression demandee
  202.  
  203. uxu1 = TABLE ; uyu1 = TABLE ; uzu1 = TABLE ;
  204. epxxu1 = TABLE ; epyyu1 = TABLE ; epzzu1 = TABLE ;
  205. gaxyu1 = TABLE ; gaxzu1 = TABLE ; gayzu1 = TABLE ;
  206. eixxu1 = TABLE ; eiyyu1 = TABLE ; eizzu1 = TABLE ;
  207. gixyu1 = TABLE ; gixzu1 = TABLE ; giyzu1 = TABLE ;
  208. smxxu1 = TABLE ; smyyu1 = TABLE ; smzzu1 = TABLE ;
  209. smxyu1 = TABLE ; smxzu1 = TABLE ; smyzu1 = TABLE ;
  210. ec0u1  = TABLE ; esw0u1 = TABLE ; pu1    = TABLE ; qtldu1 = TABLE ;
  211.  
  212. uxu2 = TABLE ; uyu2 = TABLE ; uzu2 = TABLE ;
  213. epxxu2 = TABLE ; epyyu2 = TABLE ; epzzu2 = TABLE ;
  214. gaxyu2 = TABLE ; gaxzu2 = TABLE ; gayzu2 = TABLE ;
  215. eixxu2 = TABLE ; eiyyu2 = TABLE ; eizzu2 = TABLE ;
  216. gixyu2 = TABLE ; gixzu2 = TABLE ; giyzu2 = TABLE ;
  217. smxxu2 = TABLE ; smyyu2 = TABLE ; smzzu2 = TABLE ;
  218. smxyu2 = TABLE ; smxzu2 = TABLE ; smyzu2 = TABLE ;
  219. ec0u2  = TABLE ; esw0u2 = TABLE ; pu2    = TABLE ; qtldu2 = TABLE ;
  220. psupu2 = TABLE ; qsupu2 = TABLE ;
  221.  
  222. LDEPL = MOTS 'UX  ' 'UY  ' 'UZ  ' ;
  223. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  224. LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  225. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  226. LVARI = MOTS 'EC0 ' 'ESW0' 'P   ' 'QTLD' ;
  227. LVSUP = MOTS 'PSUP' 'QSUP' ;
  228.  
  229. idefo = 4 ;
  230. idein = 4 ;
  231. icont = 4 ;
  232. ivar1 = 1 ;
  233. ivar2 = 2 ;
  234. ivar3 = 3 ;
  235. ivar4 = 4 ;
  236. ivas1 = 1 ;
  237. ivas2 = 2 ;
  238.  
  239. ierr = 0 ;
  240. critere = 1.d-8 ;
  241. IMESS = 0 ;
  242. ***IMESS = 1 ;
  243.  
  244. nbno  = 8 ;
  245. nbgau = 8 ;
  246.  
  247. netat = DIME (tabuti1.'TEMPS') ;
  248. ietat = -1 ;
  249.  
  250.  
  251. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  252.  
  253. REPETER BECART netat ;
  254.  
  255.    ietat = ietat + 1 ;
  256.  
  257.    temps1 = tabuti1.'TEMPS'.ietat ;
  258.    temps2 = tabuti2.'TEMPS'.ietat ;
  259.    SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  260.       ierr = 1 ;
  261.       QUITTER BECART ;
  262.    FINSI ;
  263.  
  264. *  Extraction des resultats de calcul : Modele utilisateur no21---------
  265.  
  266. *  Deplacements
  267.    uxu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  268.    uxu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  269.    uxu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  270.    uxu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  271.    uxu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  272.    uxu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  273.    uxu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  274.    uxu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  275.  
  276.    uyu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  277.    uyu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  278.    uyu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  279.    uyu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  280.    uyu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  281.    uyu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  282.    uyu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  283.    uyu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  284.  
  285.    uzu1.1 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  286.    uzu1.2 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  287.    uzu1.3 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  288.    uzu1.4 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  289.    uzu1.5 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  290.    uzu1.6 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  291.    uzu1.7 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  292.    uzu1.8 = EXTR (tabuti1.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  293.  
  294. *  Deformations totales
  295.    igau = 0 ;
  296.    REPETER BEPXXU1 nbgau ;
  297.       igau = igau + 1 ;
  298.       epxxu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau;
  299.    FIN BEPXXU1 ;
  300.  
  301.    igau = 0 ;
  302.    REPETER BEPYYU1 nbgau ;
  303.       igau = igau + 1 ;
  304.       epyyu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau;
  305.    FIN BEPYYU1 ;
  306.  
  307.    igau = 0 ;
  308.    REPETER BEPZZU1 nbgau ;
  309.       igau = igau + 1 ;
  310.       epzzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau;
  311.    FIN BEPZZU1 ;
  312.  
  313.    igau = 0 ;
  314.    REPETER BGAXYU1 nbgau ;
  315.       igau = igau + 1 ;
  316.       gaxyu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau;
  317.    FIN BGAXYU1 ;
  318.  
  319.    igau = 0 ;
  320.    REPETER BGAXZU1 nbgau ;
  321.       igau = igau + 1 ;
  322.       gaxzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau;
  323.    FIN BGAXZU1 ;
  324.  
  325.    igau = 0 ;
  326.    REPETER BGAYZU1 nbgau ;
  327.       igau = igau + 1 ;
  328.       gayzu1.igau = EXTR (tabuti1.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau;
  329.    FIN BGAYZU1 ;
  330.  
  331. *  Deformations inelastiques
  332.    igau = 0 ;
  333.    REPETER BEIXXU1 nbgau ;
  334.       igau = igau + 1 ;
  335.       eixxu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  336.                              'EIXX' 1 1 igau ;
  337.    FIN BEIXXU1 ;
  338.  
  339.    igau = 0 ;
  340.    REPETER BEIYYU1 nbgau ;
  341.       igau = igau + 1 ;
  342.       eiyyu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  343.                              'EIYY' 1 1 igau ;
  344.    FIN BEIYYU1 ;
  345.  
  346.    igau = 0 ;
  347.    REPETER BEIZZU1 nbgau ;
  348.       igau = igau + 1 ;
  349.       eizzu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  350.                              'EIZZ' 1 1 igau ;
  351.    FIN BEIZZU1 ;
  352.  
  353.    igau = 0 ;
  354.    REPETER BGIXYU1 nbgau ;
  355.       igau = igau + 1 ;
  356.       gixyu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  357.                              'GIXY' 1 1 igau ;
  358.    FIN BGIXYU1 ;
  359.  
  360.    igau = 0 ;
  361.    REPETER BGIXZU1 nbgau ;
  362.       igau = igau + 1 ;
  363.       gixzu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  364.                              'GIXZ' 1 1 igau ;
  365.    FIN BGIXZU1 ;
  366.  
  367.    igau = 0 ;
  368.    REPETER BGIYZU1 nbgau ;
  369.       igau = igau + 1 ;
  370.       giyzu1.igau = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ietat)
  371.                              'GIYZ' 1 1 igau ;
  372.    FIN BGIYZU1 ;
  373.  
  374. *  Contraintes
  375.    igau = 0 ;
  376.    REPETER BSMXXU1 nbgau ;
  377.       igau = igau + 1 ;
  378.       smxxu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  379.    FIN BSMXXU1 ;
  380.  
  381.    igau = 0 ;
  382.    REPETER BSMYYU1 nbgau ;
  383.       igau = igau + 1 ;
  384.       smyyu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  385.    FIN BSMYYU1 ;
  386.  
  387.    igau = 0 ;
  388.    REPETER BSMZZU1 nbgau ;
  389.       igau = igau + 1 ;
  390.       smzzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  391.    FIN BSMZZU1 ;
  392.  
  393.    igau = 0 ;
  394.    REPETER BSMXYU1 nbgau ;
  395.       igau = igau + 1 ;
  396.       smxyu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  397.    FIN BSMXYU1 ;
  398.  
  399.    igau = 0 ;
  400.    REPETER BSMXZU1 nbgau ;
  401.       igau = igau + 1 ;
  402.       smxzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  403.    FIN BSMXZU1 ;
  404.  
  405.    igau = 0 ;
  406.    REPETER BSMYZU1 nbgau ;
  407.       igau = igau + 1 ;
  408.       smyzu1.igau = EXTR (tabuti1.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  409.    FIN BSMYZU1 ;
  410.  
  411. *  Variables internes
  412.    igau = 0 ;
  413.    REPETER BEC0U1 nbgau ;
  414.       igau = igau + 1 ;
  415.       ec0u1.igau =  EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  416.                          'EC0 ' 1 1 igau ;
  417.    FIN BEC0U1 ;
  418.  
  419.    igau = 0 ;
  420.    REPETER BESW0U1 nbgau ;
  421.       igau = igau + 1 ;
  422.       esw0u1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  423.                          'ESW0' 1 1 igau ;
  424.    FIN BESW0U1 ;
  425.  
  426.    igau = 0 ;
  427.    REPETER BPU1 nbgau ;
  428.       igau = igau + 1 ;
  429.       pu1.igau =    EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  430.                          'P   ' 1 1 igau ;
  431.    FIN BPU1 ;
  432.  
  433.    igau = 0 ;
  434.    REPETER BQTLDU1 nbgau ;
  435.       igau = igau + 1 ;
  436.       qtldu1.igau = EXTR (tabuti1.'VARIABLES_INTERNES'.ietat)
  437.                          'QTLD' 1 1 igau ;
  438.    FIN BQTLDU1 ;
  439.  
  440. *  Extraction des resultats de calcul : Modele utilisateur no23---------
  441.  
  442. *  Deplacements
  443.    uxu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  444.    uxu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  445.    uxu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  446.    uxu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  447.    uxu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  448.    uxu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  449.    uxu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  450.    uxu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  451.  
  452.    uyu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  453.    uyu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  454.    uyu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  455.    uyu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  456.    uyu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  457.    uyu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  458.    uyu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  459.    uyu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  460.  
  461.    uzu2.1 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  462.    uzu2.2 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  463.    uzu2.3 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  464.    uzu2.4 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  465.    uzu2.5 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  466.    uzu2.6 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  467.    uzu2.7 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  468.    uzu2.8 = EXTR (tabuti2.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  469.  
  470. *  Deformations totales
  471.    igau = 0 ;
  472.    REPETER BEPXXU2 nbgau ;
  473.       igau = igau + 1 ;
  474.       epxxu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau;
  475.    FIN BEPXXU2 ;
  476.  
  477.    igau = 0 ;
  478.    REPETER BEPYYU2 nbgau ;
  479.       igau = igau + 1 ;
  480.       epyyu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau;
  481.    FIN BEPYYU2 ;
  482.  
  483.    igau = 0 ;
  484.    REPETER BEPZZU2 nbgau ;
  485.       igau = igau + 1 ;
  486.       epzzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau;
  487.    FIN BEPZZU2 ;
  488.  
  489.    igau = 0 ;
  490.    REPETER BGAXYU2 nbgau ;
  491.       igau = igau + 1 ;
  492.       gaxyu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau;
  493.    FIN BGAXYU2 ;
  494.  
  495.    igau = 0 ;
  496.    REPETER BGAXZU2 nbgau ;
  497.       igau = igau + 1 ;
  498.       gaxzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau;
  499.    FIN BGAXZU2 ;
  500.  
  501.    igau = 0 ;
  502.    REPETER BGAYZU2 nbgau ;
  503.       igau = igau + 1 ;
  504.       gayzu2.igau = EXTR (tabuti2.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau;
  505.    FIN BGAYZU2 ;
  506.  
  507. *  Deformations inelastiques
  508.    igau = 0 ;
  509.    REPETER BEIXXU2 nbgau ;
  510.       igau = igau + 1 ;
  511.       eixxu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  512.                               'EIXX' 1 1 igau ;
  513.    FIN BEIXXU2 ;
  514.  
  515.    igau = 0 ;
  516.    REPETER BEIYYU2 nbgau ;
  517.       igau = igau + 1 ;
  518.       eiyyu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  519.                               'EIYY' 1 1 igau ;
  520.    FIN BEIYYU2 ;
  521.  
  522.    igau = 0 ;
  523.    REPETER BEIZZU2 nbgau ;
  524.       igau = igau + 1 ;
  525.       eizzu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  526.                               'EIZZ' 1 1 igau ;
  527.    FIN BEIZZU2 ;
  528.  
  529.    igau = 0 ;
  530.    REPETER BGIXYU2 nbgau ;
  531.       igau = igau + 1 ;
  532.       gixyu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  533.                               'GIXY' 1 1 igau ;
  534.    FIN BGIXYU2 ;
  535.  
  536.    igau = 0 ;
  537.    REPETER BGIXZU2 nbgau ;
  538.       igau = igau + 1 ;
  539.       gixzu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  540.                               'GIXZ' 1 1 igau ;
  541.    FIN BGIXZU2 ;
  542.  
  543.    igau = 0 ;
  544.    REPETER BGIYZU2 nbgau ;
  545.       igau = igau + 1 ;
  546.       giyzu2.igau = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ietat)
  547.                               'GIYZ' 1 1 igau ;
  548.    FIN BGIYZU2 ;
  549.  
  550. *  Contraintes
  551.    igau = 0 ;
  552.    REPETER BSMXXU2 nbgau ;
  553.       igau = igau + 1 ;
  554.       smxxu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  555.    FIN BSMXXU2 ;
  556.  
  557.    igau = 0 ;
  558.    REPETER BSMYYU2 nbgau ;
  559.       igau = igau + 1 ;
  560.       smyyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  561.    FIN BSMYYU2 ;
  562.  
  563.    igau = 0 ;
  564.    REPETER BSMZZU2 nbgau ;
  565.       igau = igau + 1 ;
  566.       smzzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  567.    FIN BSMZZU2 ;
  568.  
  569.    igau = 0 ;
  570.    REPETER BSMXYU2 nbgau ;
  571.       igau = igau + 1 ;
  572.       smxyu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  573.    FIN BSMXYU2 ;
  574.  
  575.    igau = 0 ;
  576.    REPETER BSMXZU2 nbgau ;
  577.       igau = igau + 1 ;
  578.       smxzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  579.    FIN BSMXZU2 ;
  580.  
  581.    igau = 0 ;
  582.    REPETER BSMYZU2 nbgau ;
  583.       igau = igau + 1 ;
  584.       smyzu2.igau = EXTR (tabuti2.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  585.    FIN BSMYZU2 ;
  586.  
  587. *  Variables internes
  588.    igau = 0 ;
  589.    REPETER BEC0U2 nbgau ;
  590.       igau = igau + 1 ;
  591.       ec0u2.igau =  EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  592.                          'EC0 ' 1 1 igau ;
  593.    FIN BEC0U2 ;
  594.  
  595.    igau = 0 ;
  596.    REPETER BESW0U2 nbgau ;
  597.       igau = igau + 1 ;
  598.       esw0u2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  599.                          'ESW0' 1 1 igau ;
  600.    FIN BESW0U2 ;
  601.  
  602.    igau = 0 ;
  603.    REPETER BPU2 nbgau ;
  604.       igau = igau + 1 ;
  605.       pu2.igau =    EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  606.                          'P   ' 1 1 igau ;
  607.    FIN BPU2 ;
  608.  
  609.    igau = 0 ;
  610.    REPETER BQTLDU2 nbgau ;
  611.       igau = igau + 1 ;
  612.       qtldu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  613.                          'QTLD' 1 1 igau ;
  614.    FIN BQTLDU2 ;
  615.  
  616. *  Variables internes supplementaires
  617.    igau = 0 ;
  618.    REPETER BPSUPU2 nbgau ;
  619.       igau = igau + 1 ;
  620.       psupu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  621.                          'PSUP' 1 1 igau ;
  622.    FIN BPSUPU2 ;
  623.  
  624.    igau = 0 ;
  625.    REPETER BQSUPU2 nbgau ;
  626.       igau = igau + 1 ;
  627.       qsupu2.igau = EXTR (tabuti2.'VARIABLES_INTERNES'.ietat)
  628.                          'QSUP' 1 1 igau ;
  629.    FIN BQSUPU2 ;
  630.  
  631. *  Comparaison des deplacements-----------------------------------------
  632.  
  633.    icomp = 1 ;
  634.    ino = 0 ;
  635.    REPETER BUX nbno ;
  636.       ino = ino + 1 ;
  637.       SI (EGA uxu1.ino 0.0) ;
  638.          SI (NON (EGA uxu2.ino 0.0)) ;
  639.             SI (EGA IMESS 1) ;
  640.             Message 'Erreur absolue' uxu2.ino
  641.                     'Composante ' (EXTR LDEPL icomp)
  642.                     ' Noeud' ino 'Etat' ietat ;
  643.             Message ;
  644.             FINSI ;
  645. ************ierr = 10 ; errabs = uxu2.ino ; QUITTER BECART ;
  646.          FINSI ;
  647.       SINON ;
  648.          dux = ABS (uxu1.ino - uxu2.ino) ;
  649.          errrel = ABS (dux / uxu1.ino) ;
  650.          SI (errrel > critere) ;
  651.             SI (EGA IMESS 1) ;
  652.             Message 'Erreur relative' errrel
  653.                     'Composante ' (EXTR LDEPL icomp)
  654.                     ' Noeud' ino 'Etat' ietat ;
  655.             Message 'uxu1.ino =' uxu1.ino 'uxu2.ino =' uxu2.ino ;
  656.             Message ;
  657.             FINSI ;
  658. ************ierr = 11 ; QUITTER BECART ;
  659.          FINSI ;
  660.       FINSI ;
  661.    FIN BUX ;
  662.  
  663.    icomp = 2 ;
  664.    ino = 0 ;
  665.    REPETER BUY nbno ;
  666.       ino = ino + 1 ;
  667.       SI (EGA uyu1.ino 0.0) ;
  668.          SI (NON (EGA uyu2.ino 0.0)) ;
  669.             SI (EGA IMESS 1) ;
  670.             Message 'Erreur absolue' uyu2.ino
  671.                     'Composante ' (EXTR LDEPL icomp)
  672.                     ' Noeud' ino 'Etat' ietat ;
  673.             Message ;
  674.             FINSI ;
  675. ************ierr = 10 ; errabs = uyu2.ino ; QUITTER BECART ;
  676.          FINSI ;
  677.       SINON ;
  678.          duy = ABS (uyu1.ino - uyu2.ino) ;
  679.          errrel = ABS (duy / uyu1.ino) ;
  680.          SI (errrel > critere) ;
  681.             SI (EGA IMESS 1) ;
  682.             Message 'Erreur relative' errrel
  683.                     'Composante ' (EXTR LDEPL icomp)
  684.                     ' Noeud' ino 'Etat' ietat ;
  685.             Message 'uyu1.ino =' uyu1.ino 'uyu2.ino =' uyu2.ino ;
  686.             Message ;
  687.             FINSI ;
  688. ************ierr = 11 ; QUITTER BECART ;
  689.          FINSI ;
  690.       FINSI ;
  691.    FIN BUY ;
  692.  
  693.    icomp = 3 ;
  694.    ino = 0 ;
  695.    REPETER BUZ nbno ;
  696.       ino = ino + 1 ;
  697.       SI (EGA uzu1.ino 0.0) ;
  698.          SI (NON (EGA uzu2.ino 0.0)) ;
  699.             SI (EGA IMESS 1) ;
  700.             Message 'Erreur absolue' uzu2.ino
  701.                     'Composante ' (EXTR LDEPL icomp)
  702.                     ' Noeud' ino 'Etat' ietat ;
  703.             Message ;
  704.             FINSI ;
  705. ************ierr = 10 ; errabs = uzu2.ino ; QUITTER BECART ;
  706.          FINSI ;
  707.       SINON ;
  708.          duz = ABS (uzu1.ino - uzu2.ino) ;
  709.          errrel = ABS (duz / uzu1.ino) ;
  710.          SI (errrel > critere) ;
  711.             SI (EGA IMESS 1) ;
  712.             Message 'Erreur relative' errrel
  713.                     'Composante ' (EXTR LDEPL icomp)
  714.                     ' Noeud' ino 'Etat' ietat ;
  715.             Message 'uzu1.ino =' uzu1.ino 'uzu2.ino =' uzu2.ino ;
  716.             Message ;
  717.             FINSI ;
  718. ************ierr = 11 ; QUITTER BECART ;
  719.          FINSI ;
  720.       FINSI ;
  721.    FIN BUZ ;
  722.  
  723. *  Comparaison des deformations-----------------------------------------
  724.  
  725.    icomp = 1 ;
  726.    igau = 0 ;
  727.    REPETER BEPXX nbgau ;
  728.       igau = igau + 1 ;
  729.       SI (EGA epxxu1.igau 0.0) ;
  730.          SI (NON (EGA epxxu2.igau 0.0)) ;
  731.             SI (EGA IMESS 1) ;
  732.             Message 'Erreur absolue' epxxu2.igau
  733.                     'Composante ' (EXTR LDEFO icomp)
  734.                     ' Point Gauss' igau 'Etat' ietat ;
  735.             Message ;
  736.             FINSI ;
  737. ************ierr = 20 ; errabs = epxxu2.igau ; QUITTER BECART ;
  738.          FINSI ;
  739.       SINON ;
  740.          depxx = ABS (epxxu1.igau - epxxu2.igau) ;
  741.          errrel = ABS (depxx / epxxu1.igau) ;
  742.          SI (errrel > critere) ;
  743.             SI (EGA IMESS 1) ;
  744.             Message 'Erreur relative' errrel
  745.                     'Composante ' (EXTR LDEFO icomp)
  746.                     ' Point Gauss' igau 'Etat' ietat ;
  747.             Message 'epxxu1.igau =' epxxu1.igau
  748.                     'epxxu2.igau =' epxxu2.igau;
  749.             Message ;
  750.             FINSI ;
  751. ************ierr = 21 ; QUITTER BECART ;
  752.          FINSI ;
  753.       FINSI ;
  754.    FIN BEPXX ;
  755.  
  756.    icomp = 2 ;
  757.    igau = 0 ;
  758.    REPETER BEPYY nbgau ;
  759.       igau = igau + 1 ;
  760.       SI (EGA epyyu1.igau 0.0) ;
  761.          SI (NON (EGA epyyu2.igau 0.0)) ;
  762.             SI (EGA IMESS 1) ;
  763.             Message 'Erreur absolue' epyyu2.igau
  764.                     'Composante ' (EXTR LDEFO icomp)
  765.                     ' Point Gauss' igau 'Etat' ietat ;
  766.             Message ;
  767.             FINSI ;
  768. ************ierr = 20 ; errabs = epyyu2.igau ; QUITTER BECART ;
  769.          FINSI ;
  770.       SINON ;
  771.          depyy = ABS (epyyu1.igau - epyyu2.igau) ;
  772.          errrel = ABS (depyy / epyyu1.igau) ;
  773.          SI (errrel > critere) ;
  774.             SI (EGA IMESS 1) ;
  775.             Message 'Erreur relative' errrel
  776.                     'Composante ' (EXTR LDEFO icomp)
  777.                     ' Point Gauss' igau 'Etat' ietat ;
  778.             Message 'epyyu1.igau =' epyyu1.igau
  779.                     'epyyu2.igau =' epyyu2.igau;
  780.             Message ;
  781.             FINSI ;
  782. ************ierr = 21 ; QUITTER BECART ;
  783.          FINSI ;
  784.       FINSI ;
  785.    FIN BEPYY ;
  786.  
  787.    icomp = 3 ;
  788.    igau = 0 ;
  789.    REPETER BEPZZ nbgau ;
  790.       igau = igau + 1 ;
  791.       SI (EGA epzzu1.igau 0.0) ;
  792.          SI (NON (EGA epzzu2.igau 0.0)) ;
  793.             SI (EGA IMESS 1) ;
  794.             Message 'Erreur absolue' epzzu2.igau
  795.                     'Composante ' (EXTR LDEFO icomp)
  796.                     ' Point Gauss' igau 'Etat' ietat ;
  797.             Message ;
  798.             FINSI ;
  799. ************ierr = 20 ; errabs = epzzu2.igau ; QUITTER BECART ;
  800.          FINSI ;
  801.       SINON ;
  802.          depzz = ABS (epzzu1.igau - epzzu2.igau) ;
  803.          errrel = ABS (depzz / epzzu1.igau) ;
  804.          SI (errrel > critere) ;
  805.             SI (EGA IMESS 1) ;
  806.             Message 'Erreur relative' errrel
  807.                     'Composante ' (EXTR LDEFO icomp)
  808.                     ' Point Gauss' igau 'Etat' ietat ;
  809.             Message 'epzzu1.igau =' epzzu1.igau
  810.                     'epzzu2.igau =' epzzu2.igau;
  811.             Message ;
  812.             FINSI ;
  813. ************ierr = 21 ; QUITTER BECART ;
  814.          FINSI ;
  815.       FINSI ;
  816.    FIN BEPZZ ;
  817.  
  818.    icomp = 4 ;
  819.    igau = 0 ;
  820.    REPETER BGAXY nbgau ;
  821.       igau = igau + 1 ;
  822.       SI (EGA gaxyu1.igau 0.0) ;
  823.          SI (NON (EGA gaxyu2.igau 0.0)) ;
  824. ************SI (EGA IMESS 1) ;
  825. ************Message 'Erreur absolue' gaxyu2.igau
  826. ************        'Composante ' (EXTR LDEFO icomp)
  827. ************        ' Point Gauss' igau 'Etat' ietat ;
  828. ************Message ;
  829. ************FINSI ;
  830.             ierr = 20 ; errabs = gaxyu2.igau ; QUITTER BECART ;
  831.          FINSI ;
  832.       SINON ;
  833.          dgaxy = ABS (gaxyu1.igau - gaxyu2.igau) ;
  834.          errrel = ABS (dgaxy / gaxyu1.igau) ;
  835.          SI (errrel > critere) ;
  836. ************SI (EGA IMESS 1) ;
  837. ************Message 'Erreur relative' errrel
  838. ************        'Composante ' (EXTR LDEFO icomp)
  839. ************        ' Point Gauss' igau 'Etat' ietat ;
  840. ************Message 'gaxyu1.igau =' gaxyu1.igau
  841. ************        'gaxyu2.igau =' gaxyu2.igau;
  842. ************Message ;
  843. ************FINSI ;
  844.             ierr = 21 ; QUITTER BECART ;
  845.          FINSI ;
  846.       FINSI ;
  847.    FIN BGAXY ;
  848.  
  849.    icomp = 5 ;
  850.    igau = 0 ;
  851.    REPETER BGAXZ nbgau ;
  852.       igau = igau + 1 ;
  853.       SI (EGA gaxzu1.igau 0.0) ;
  854.          SI (NON (EGA gaxzu2.igau 0.0)) ;
  855.             SI (EGA IMESS 1) ;
  856.             Message 'Erreur absolue' gaxzu2.igau
  857.                     'Composante ' (EXTR LDEFO icomp)
  858.                     ' Point Gauss' igau 'Etat' ietat ;
  859.             Message ;
  860.             FINSI ;
  861. ************ierr = 20 ; errabs = gaxzu2.igau ; QUITTER BECART ;
  862.          FINSI ;
  863.       SINON ;
  864.          dgaxz = ABS (gaxzu1.igau - gaxzu2.igau) ;
  865.          errrel = ABS (dgaxz / gaxzu1.igau) ;
  866.          SI (errrel > critere) ;
  867.             SI (EGA IMESS 1) ;
  868.             Message 'Erreur relative' errrel
  869.                     'Composante ' (EXTR LDEFO icomp)
  870.                     ' Point Gauss' igau 'Etat' ietat ;
  871.             Message 'gaxzu1.igau =' gaxzu1.igau
  872.                     'gaxzu2.igau =' gaxzu2.igau;
  873.             Message ;
  874.             FINSI ;
  875. ************ierr = 21 ; QUITTER BECART ;
  876.          FINSI ;
  877.       FINSI ;
  878.    FIN BGAXZ ;
  879.  
  880.    icomp = 6 ;
  881.    igau = 0 ;
  882.    REPETER BGAYZ nbgau ;
  883.       igau = igau + 1 ;
  884.       SI (EGA gayzu1.igau 0.0) ;
  885.          SI (NON (EGA gayzu2.igau 0.0)) ;
  886.             SI (EGA IMESS 1) ;
  887.             Message 'Erreur absolue' gayzu2.igau
  888.                     'Composante ' (EXTR LDEFO icomp)
  889.                     ' Point Gauss' igau 'Etat' ietat ;
  890.             Message ;
  891.             FINSI ;
  892. ************ierr = 20 ; errabs = gayzu2.igau ; QUITTER BECART ;
  893.          FINSI ;
  894.       SINON ;
  895.          dgayz = ABS (gayzu1.igau - gayzu2.igau) ;
  896.          errrel = ABS (dgayz / gayzu1.igau) ;
  897.          SI (errrel > critere) ;
  898.             SI (EGA IMESS 1) ;
  899.             Message 'Erreur relative' errrel
  900.                     'Composante ' (EXTR LDEFO icomp)
  901.                     ' Point Gauss' igau 'Etat' ietat ;
  902.             Message 'gayzu1.igau =' gayzu1.igau
  903.                     'gayzu2.igau =' gayzu2.igau;
  904.             Message ;
  905.             FINSI ;
  906. ************ierr = 21 ; QUITTER BECART ;
  907.          FINSI ;
  908.       FINSI ;
  909.    FIN BGAYZ ;
  910.  
  911. *  Comparaison des deformations inelastiques----------------------------
  912.  
  913.    icomp = 1 ;
  914.    igau = 0 ;
  915.    REPETER BEIXX nbgau ;
  916.       igau = igau + 1 ;
  917.       SI (EGA eixxu1.igau 0.0) ;
  918.          SI (NON (EGA eixxu2.igau 0.0)) ;
  919.             SI (EGA IMESS 1) ;
  920.             Message 'Erreur absolue' eixxu2.igau
  921.                     'Composante ' (EXTR LDEIN icomp)
  922.                     ' Point Gauss' igau 'Etat' ietat ;
  923.             Message ;
  924.             FINSI ;
  925. ************ierr = 30 ; errabs = eixxu2.igau ; QUITTER BECART ;
  926.          FINSI ;
  927.       SINON ;
  928.          deixx = ABS (eixxu1.igau - eixxu2.igau) ;
  929.          errrel = ABS (deixx / eixxu1.igau) ;
  930.          SI (errrel > critere) ;
  931.             SI (EGA IMESS 1) ;
  932.             Message 'Erreur relative' errrel
  933.                     'Composante ' (EXTR LDEIN icomp)
  934.                     ' Point Gauss' igau 'Etat' ietat ;
  935.             Message 'eixxu1.igau =' eixxu1.igau
  936.                     'eixxu2.igau =' eixxu2.igau;
  937.             Message ;
  938.             FINSI ;
  939. ************ierr = 31 ; QUITTER BECART ;
  940.          FINSI ;
  941.       FINSI ;
  942.    FIN BEIXX ;
  943.  
  944.    icomp = 2 ;
  945.    igau = 0 ;
  946.    REPETER BEIYY nbgau ;
  947.       igau = igau + 1 ;
  948.       SI (EGA eiyyu1.igau 0.0) ;
  949.          SI (NON (EGA eiyyu2.igau 0.0)) ;
  950.             SI (EGA IMESS 1) ;
  951.             Message 'Erreur absolue' eiyyu2.igau
  952.                     'Composante ' (EXTR LDEIN icomp)
  953.                     ' Point Gauss' igau 'Etat' ietat ;
  954.             Message ;
  955.             FINSI ;
  956. ************ierr = 30 ; errabs = eiyyu2.igau ; QUITTER BECART ;
  957.          FINSI ;
  958.       SINON ;
  959.          deiyy = ABS (eiyyu1.igau - eiyyu2.igau) ;
  960.          errrel = ABS (deiyy / eiyyu1.igau) ;
  961.          SI (errrel > critere) ;
  962.             SI (EGA IMESS 1) ;
  963.             Message 'Erreur relative' errrel
  964.                     'Composante ' (EXTR LDEIN icomp)
  965.                     ' Point Gauss' igau 'Etat' ietat ;
  966.             Message 'eiyyu1.igau =' eiyyu1.igau
  967.                     'eiyyu2.igau =' eiyyu2.igau;
  968.             Message ;
  969.             FINSI ;
  970. ************ierr = 31 ; QUITTER BECART ;
  971.          FINSI ;
  972.       FINSI ;
  973.    FIN BEIYY ;
  974.  
  975.    icomp = 3 ;
  976.    igau = 0 ;
  977.    REPETER BEIZZ nbgau ;
  978.       igau = igau + 1 ;
  979.       SI (EGA eizzu1.igau 0.0) ;
  980.          SI (NON (EGA eizzu2.igau 0.0)) ;
  981.             SI (EGA IMESS 1) ;
  982.             Message 'Erreur absolue' eizzu2.igau
  983.                     'Composante ' (EXTR LDEIN icomp)
  984.                     ' Point Gauss' igau 'Etat' ietat ;
  985.             Message ;
  986.             FINSI ;
  987. ************ierr = 30 ; errabs = eizzu2.igau ; QUITTER BECART ;
  988.          FINSI ;
  989.       SINON ;
  990.          deizz = ABS (eizzu1.igau - eizzu2.igau) ;
  991.          errrel = ABS (deizz / eizzu1.igau) ;
  992.          SI (errrel > critere) ;
  993.             SI (EGA IMESS 1) ;
  994.             Message 'Erreur relative' errrel
  995.                     'Composante ' (EXTR LDEIN icomp)
  996.                     ' Point Gauss' igau 'Etat' ietat ;
  997.             Message 'eizzu1.igau =' eizzu1.igau
  998.                     'eizzu2.igau =' eizzu2.igau;
  999.             Message ;
  1000.             FINSI ;
  1001. ************ierr = 31 ; QUITTER BECART ;
  1002.          FINSI ;
  1003.       FINSI ;
  1004.    FIN BEIZZ ;
  1005.  
  1006.    icomp = 4 ;
  1007.    igau = 0 ;
  1008.    REPETER BGIXY nbgau ;
  1009.       igau = igau + 1 ;
  1010.       SI (EGA gixyu1.igau 0.0) ;
  1011.          SI (NON (EGA gixyu2.igau 0.0)) ;
  1012. ************SI (EGA IMESS 1) ;
  1013. ************Message 'Erreur absolue' gixyu2.igau
  1014. ************        'Composante ' (EXTR LDEIN icomp)
  1015. ************        ' Point Gauss' igau 'Etat' ietat ;
  1016. ************Message ;
  1017. ************FINSI ;
  1018.             ierr = 30 ; errabs = gixyu2.igau ; QUITTER BECART ;
  1019.          FINSI ;
  1020.       SINON ;
  1021.          dgixy = ABS (gixyu1.igau - gixyu2.igau) ;
  1022.          errrel = ABS (dgixy / gixyu1.igau) ;
  1023.          SI (errrel > critere) ;
  1024. ************SI (EGA IMESS 1) ;
  1025. ************Message 'Erreur relative' errrel
  1026. ************        'Composante ' (EXTR LDEIN icomp)
  1027. ************        ' Point Gauss' igau 'Etat' ietat ;
  1028. ************Message 'gixyu1.igau =' gixyu1.igau
  1029. ************        'gixyu2.igau =' gixyu2.igau;
  1030. ************Message ;
  1031. ************FINSI ;
  1032.             ierr = 31 ; QUITTER BECART ;
  1033.          FINSI ;
  1034.       FINSI ;
  1035.    FIN BGIXY ;
  1036.  
  1037.    icomp = 5 ;
  1038.    igau = 0 ;
  1039.    REPETER BGIXZ nbgau ;
  1040.       igau = igau + 1 ;
  1041.       SI (EGA gixzu1.igau 0.0) ;
  1042.          SI (NON (EGA gixzu2.igau 0.0)) ;
  1043.             SI (EGA IMESS 1) ;
  1044.             Message 'Erreur absolue' gixzu2.igau
  1045.                     'Composante ' (EXTR LDEIN icomp)
  1046.                     ' Point Gauss' igau 'Etat' ietat ;
  1047.             Message ;
  1048.             FINSI ;
  1049. ************ierr = 30 ; errabs = gixzu2.igau ; QUITTER BECART ;
  1050.          FINSI ;
  1051.       SINON ;
  1052.          dgixz = ABS (gixzu1.igau - gixzu2.igau) ;
  1053.          errrel = ABS (dgixz / gixzu1.igau) ;
  1054.          SI (errrel > critere) ;
  1055.             SI (EGA IMESS 1) ;
  1056.             Message 'Erreur relative' errrel
  1057.                     'Composante ' (EXTR LDEIN icomp)
  1058.                     ' Point Gauss' igau 'Etat' ietat ;
  1059.             Message 'gixzu1.igau =' gixzu1.igau
  1060.                     'gixzu2.igau =' gixzu2.igau;
  1061.             Message ;
  1062.             FINSI ;
  1063. ************ierr = 31 ; QUITTER BECART ;
  1064.          FINSI ;
  1065.       FINSI ;
  1066.    FIN BGIXZ ;
  1067.  
  1068.    icomp = 6 ;
  1069.    igau = 0 ;
  1070.    REPETER BGIYZ nbgau ;
  1071.       igau = igau + 1 ;
  1072.       SI (EGA giyzu1.igau 0.0) ;
  1073.          SI (NON (EGA giyzu2.igau 0.0)) ;
  1074.             SI (EGA IMESS 1) ;
  1075.             Message 'Erreur absolue' giyzu2.igau
  1076.                     'Composante ' (EXTR LDEIN icomp)
  1077.                     ' Point Gauss' igau 'Etat' ietat ;
  1078.             Message ;
  1079.             FINSI ;
  1080. ************ierr = 30 ; errabs = giyzu2.igau ; QUITTER BECART ;
  1081.          FINSI ;
  1082.       SINON ;
  1083.          dgiyz = ABS (giyzu1.igau - giyzu2.igau) ;
  1084.          errrel = ABS (dgiyz / giyzu1.igau) ;
  1085.          SI (errrel > critere) ;
  1086.             SI (EGA IMESS 1) ;
  1087.             Message 'Erreur relative' errrel
  1088.                     'Composante ' (EXTR LDEIN icomp)
  1089.                     ' Point Gauss' igau 'Etat' ietat ;
  1090.             Message 'giyzu1.igau =' giyzu1.igau
  1091.                     'giyzu2.igau =' giyzu2.igau;
  1092.             Message ;
  1093.             FINSI ;
  1094. ************ierr = 31 ; QUITTER BECART ;
  1095.          FINSI ;
  1096.       FINSI ;
  1097.    FIN BGIYZ ;
  1098.  
  1099. *  Comparaison des contraintes------------------------------------------
  1100.  
  1101.    icomp = 1 ;
  1102.    igau = 0 ;
  1103.    REPETER BSMXX nbgau ;
  1104.       igau = igau + 1 ;
  1105.       SI (EGA smxxu1.igau 0.0) ;
  1106.          SI (NON (EGA smxxu2.igau 0.0)) ;
  1107.             SI (EGA IMESS 1) ;
  1108.             Message 'Erreur absolue' smxxu2.igau
  1109.                     'Composante ' (EXTR LCONT icomp)
  1110.                     ' Point Gauss' igau 'Etat' ietat ;
  1111.             Message ;
  1112.             FINSI ;
  1113. ************ierr = 40 ; errabs = smxxu2.igau ; QUITTER BECART ;
  1114.          FINSI ;
  1115.       SINON ;
  1116.          dsmxx = ABS (smxxu1.igau - smxxu2.igau) ;
  1117.          errrel = ABS (dsmxx / smxxu1.igau) ;
  1118.          SI (errrel > critere) ;
  1119.             SI (EGA IMESS 1) ;
  1120.             Message 'Erreur relative' errrel
  1121.                     'Composante ' (EXTR LCONT icomp)
  1122.                     ' Point Gauss' igau 'Etat' ietat ;
  1123.             Message 'smxxu1.igau =' smxxu1.igau
  1124.                     'smxxu2.igau =' smxxu2.igau;
  1125.             Message ;
  1126.             FINSI ;
  1127. ************ierr = 41 ; QUITTER BECART ;
  1128.          FINSI ;
  1129.       FINSI ;
  1130.    FIN BSMXX ;
  1131.  
  1132.    icomp = 2 ;
  1133.    igau = 0 ;
  1134.    REPETER BSMYY nbgau ;
  1135.       igau = igau + 1 ;
  1136.       SI (EGA smyyu1.igau 0.0) ;
  1137.          SI (NON (EGA smyyu2.igau 0.0)) ;
  1138.             SI (EGA IMESS 1) ;
  1139.             Message 'Erreur absolue' smyyu2.igau
  1140.                     'Composante ' (EXTR LCONT icomp)
  1141.                     ' Point Gauss' igau 'Etat' ietat ;
  1142.             Message ;
  1143.             FINSI ;
  1144. ************ierr = 40 ; errabs = smyyu2.igau ; QUITTER BECART ;
  1145.          FINSI ;
  1146.       SINON ;
  1147.          dsmyy = ABS (smyyu1.igau - smyyu2.igau) ;
  1148.          errrel = ABS (dsmyy / smyyu1.igau) ;
  1149.          SI (errrel > critere) ;
  1150.             SI (EGA IMESS 1) ;
  1151.             Message 'Erreur relative' errrel
  1152.                     'Composante ' (EXTR LCONT icomp)
  1153.                     ' Point Gauss' igau 'Etat' ietat ;
  1154.             Message 'smyyu1.igau =' smyyu1.igau
  1155.                     'smyyu2.igau =' smyyu2.igau;
  1156.             Message ;
  1157.             FINSI ;
  1158. ************ierr = 41 ; QUITTER BECART ;
  1159.          FINSI ;
  1160.       FINSI ;
  1161.    FIN BSMYY ;
  1162.  
  1163.    icomp = 3 ;
  1164.    igau = 0 ;
  1165.    REPETER BSMZZ nbgau ;
  1166.       igau = igau + 1 ;
  1167.       SI (EGA smzzu1.igau 0.0) ;
  1168.          SI (NON (EGA smzzu2.igau 0.0)) ;
  1169.             SI (EGA IMESS 1) ;
  1170.             Message 'Erreur absolue' smzzu2.igau
  1171.                     'Composante ' (EXTR LCONT icomp)
  1172.                     ' Point Gauss' igau 'Etat' ietat ;
  1173.             Message ;
  1174.             FINSI ;
  1175. ************ierr = 40 ; errabs = smzzu2.igau ; QUITTER BECART ;
  1176.          FINSI ;
  1177.       SINON ;
  1178.          dsmzz = ABS (smzzu1.igau - smzzu2.igau) ;
  1179.          errrel = ABS (dsmzz / smzzu1.igau) ;
  1180.          SI (errrel > critere) ;
  1181.             SI (EGA IMESS 1) ;
  1182.             Message 'Erreur relative' errrel
  1183.                     'Composante ' (EXTR LCONT icomp)
  1184.                     ' Point Gauss' igau 'Etat' ietat ;
  1185.             Message 'smzzu1.igau =' smzzu1.igau
  1186.                     'smzzu2.igau =' smzzu2.igau;
  1187.             Message ;
  1188.             FINSI ;
  1189. ************ierr = 41 ; QUITTER BECART ;
  1190.          FINSI ;
  1191.       FINSI ;
  1192.    FIN BSMZZ ;
  1193.  
  1194.    icomp = 4 ;
  1195.    igau = 0 ;
  1196.    REPETER BSMXY nbgau ;
  1197.       igau = igau + 1 ;
  1198.       SI (EGA smxyu1.igau 0.0) ;
  1199.          SI (NON (EGA smxyu2.igau 0.0)) ;
  1200. ************SI (EGA IMESS 1) ;
  1201. ************Message 'Erreur absolue' smxyu2.igau
  1202. ************        'Composante ' (EXTR LCONT icomp)
  1203. ************        ' Point Gauss' igau 'Etat' ietat ;
  1204. ************Message ;
  1205. ************FINSI ;
  1206.             ierr = 40 ; errabs = smxyu2.igau ; QUITTER BECART ;
  1207.          FINSI ;
  1208.       SINON ;
  1209.          dsmxy = ABS (smxyu1.igau - smxyu2.igau) ;
  1210.          errrel = ABS (dsmxy / smxyu1.igau) ;
  1211.          SI (errrel > critere) ;
  1212. ************SI (EGA IMESS 1) ;
  1213. ************Message 'Erreur relative' errrel
  1214. ************        'Composante ' (EXTR LCONT icomp)
  1215. ************        ' Point Gauss' igau 'Etat' ietat ;
  1216. ************Message 'smxyu1.igau =' smxyu1.igau
  1217. ************        'smxyu2.igau =' smxyu2.igau;
  1218. ************Message ;
  1219. ************FINSI ;
  1220.             ierr = 41 ; QUITTER BECART ;
  1221.          FINSI ;
  1222.       FINSI ;
  1223.    FIN BSMXY ;
  1224.  
  1225.    icomp = 5 ;
  1226.    igau = 0 ;
  1227.    REPETER BSMXZ nbgau ;
  1228.       igau = igau + 1 ;
  1229.       SI (EGA smxzu1.igau 0.0) ;
  1230.          SI (NON (EGA smxzu2.igau 0.0)) ;
  1231.             SI (EGA IMESS 1) ;
  1232.             Message 'Erreur absolue' smxzu2.igau
  1233.                     'Composante ' (EXTR LCONT icomp)
  1234.                     ' Point Gauss' igau 'Etat' ietat ;
  1235.             Message ;
  1236.             FINSI ;
  1237. ************ierr = 40 ; errabs = smxzu2.igau ; QUITTER BECART ;
  1238.          FINSI ;
  1239.       SINON ;
  1240.          dsmxz = ABS (smxzu1.igau - smxzu2.igau) ;
  1241.          errrel = ABS (dsmxz / smxzu1.igau) ;
  1242.          SI (errrel > critere) ;
  1243.             SI (EGA IMESS 1) ;
  1244.             Message 'Erreur relative' errrel
  1245.                     'Composante ' (EXTR LCONT icomp)
  1246.                     ' Point Gauss' igau 'Etat' ietat ;
  1247.             Message 'smxzu1.igau =' smxzu1.igau
  1248.                     'smxzu2.igau =' smxzu2.igau;
  1249.             Message ;
  1250.             FINSI ;
  1251. ************ierr = 41 ; QUITTER BECART ;
  1252.          FINSI ;
  1253.       FINSI ;
  1254.    FIN BSMXZ ;
  1255.  
  1256.    icomp = 6 ;
  1257.    igau = 0 ;
  1258.    REPETER BSMYZ nbgau ;
  1259.       igau = igau + 1 ;
  1260.       SI (EGA smyzu1.igau 0.0) ;
  1261.          SI (NON (EGA smyzu2.igau 0.0)) ;
  1262.             SI (EGA IMESS 1) ;
  1263.             Message 'Erreur absolue' smyzu2.igau
  1264.                     'Composante ' (EXTR LCONT icomp)
  1265.                     ' Point Gauss' igau 'Etat' ietat ;
  1266.             Message ;
  1267.             FINSI ;
  1268. ************ierr = 40 ; errabs = smyzu2.igau ; QUITTER BECART ;
  1269.          FINSI ;
  1270.       SINON ;
  1271.          dsmyz = ABS (smyzu1.igau - smyzu2.igau) ;
  1272.          errrel = ABS (dsmyz / smyzu1.igau) ;
  1273.          SI (errrel > critere) ;
  1274.             SI (EGA IMESS 1) ;
  1275.             Message 'Erreur relative' errrel
  1276.                     'Composante ' (EXTR LCONT icomp)
  1277.                     ' Point Gauss' igau 'Etat' ietat ;
  1278.             Message 'smyzu1.igau =' smyzu1.igau
  1279.                     'smyzu2.igau =' smyzu2.igau;
  1280.             Message ;
  1281.             FINSI ;
  1282. ************ierr = 41 ; QUITTER BECART ;
  1283.          FINSI ;
  1284.       FINSI ;
  1285.    FIN BSMYZ ;
  1286.  
  1287. *  Comparaison des variables internes-----------------------------------
  1288.  
  1289.    icomp = 1 ;
  1290.    igau = 0 ;
  1291.    REPETER BEC0 nbgau ;
  1292.       igau = igau + 1 ;
  1293.       SI (EGA ec0u1.igau 0.0) ;
  1294.          SI (NON (EGA ec0u2.igau 0.0)) ;
  1295. ************SI (EGA IMESS 1) ;
  1296. ************Message 'Erreur absolue' ec0u2.igau
  1297. ************        'Composante ' (EXTR LVARI icomp)
  1298. ************        ' Point Gauss' igau 'Etat' ietat ;
  1299. ************Message ;
  1300. ************FINSI ;
  1301.             ierr = 50 ; errabs = ec0u2.igau ; QUITTER BECART ;
  1302.          FINSI ;
  1303.       SINON ;
  1304.          dec0 = ABS (ec0u1.igau - ec0u2.igau) ;
  1305.          errrel = ABS (dec0 / ec0u1.igau) ;
  1306.          SI (errrel > critere) ;
  1307. ************SI (EGA IMESS 1) ;
  1308. ************Message 'Erreur relative' errrel
  1309. ************        'Composante ' (EXTR LVARI icomp)
  1310. ************        ' Point Gauss' igau 'Etat' ietat ;
  1311. ************Message 'ec0u1.igau =' ec0u1.igau
  1312. ************        'ec0u2.igau =' ec0u2.igau;
  1313. ************Message ;
  1314. ************FINSI ;
  1315.             ierr = 51 ; QUITTER BECART ;
  1316.          FINSI ;
  1317.       FINSI ;
  1318.    FIN BEC0 ;
  1319.  
  1320.    icomp = 2 ;
  1321.    igau = 0 ;
  1322.    REPETER BESW0 nbgau ;
  1323.       igau = igau + 1 ;
  1324.       SI (EGA esw0u1.igau 0.0) ;
  1325.          SI (NON (EGA esw0u2.igau 0.0)) ;
  1326. ************SI (EGA IMESS 1) ;
  1327. ************Message 'Erreur absolue' esw0u2.igau
  1328. ************        'Composante ' (EXTR LVARI icomp)
  1329. ************        ' Point Gauss' igau 'Etat' ietat ;
  1330. ************Message ;
  1331. ************FINSI ;
  1332.             ierr = 50 ; errabs = esw0u2.igau ; QUITTER BECART ;
  1333.          FINSI ;
  1334.       SINON ;
  1335.          desw0 = ABS (esw0u1.igau - esw0u2.igau) ;
  1336.          errrel = ABS (desw0 / esw0u1.igau) ;
  1337.          SI (errrel > critere) ;
  1338. ************SI (EGA IMESS 1) ;
  1339. ************Message 'Erreur relative' errrel
  1340. ************        'Composante ' (EXTR LVARI icomp)
  1341. ************        ' Point Gauss' igau 'Etat' ietat ;
  1342. ************Message 'esw0u1.igau =' esw0u1.igau
  1343. ************        'esw0u2.igau =' esw0u2.igau;
  1344. ************Message ;
  1345. ************FINSI ;
  1346.             ierr = 51 ; QUITTER BECART ;
  1347.          FINSI ;
  1348.       FINSI ;
  1349.    FIN BESW0 ;
  1350.  
  1351.    icomp = 3 ;
  1352.    igau = 0 ;
  1353.    REPETER BP nbgau ;
  1354.       igau = igau + 1 ;
  1355.       SI (EGA pu1.igau 0.0) ;
  1356.          SI (NON (EGA pu2.igau 0.0)) ;
  1357. ************SI (EGA IMESS 1) ;
  1358. ************Message 'Erreur absolue' pu2.igau
  1359. ************        'Composante ' (EXTR LVARI icomp)
  1360. ************        ' Point Gauss' igau 'Etat' ietat ;
  1361. ************Message ;
  1362. ************FINSI ;
  1363.             ierr = 50 ; errabs = pu2.igau ; QUITTER BECART ;
  1364.          FINSI ;
  1365.       SINON ;
  1366.          dp = ABS (pu1.igau - pu2.igau) ;
  1367.          errrel = ABS (dp / pu1.igau) ;
  1368.          SI (errrel > critere) ;
  1369. ************SI (EGA IMESS 1) ;
  1370. ************Message 'Erreur relative' errrel
  1371. ************        'Composante ' (EXTR LVARI icomp)
  1372. ************        ' Point Gauss' igau 'Etat' ietat ;
  1373. ************Message 'pu1.igau =' pu1.igau
  1374. ************        'pu2.igau =' pu2.igau;
  1375. ************Message ;
  1376. ************FINSI ;
  1377.             ierr = 51 ; QUITTER BECART ;
  1378.          FINSI ;
  1379.       FINSI ;
  1380.    FIN BP ;
  1381.  
  1382.    icomp = 4 ;
  1383.    igau = 0 ;
  1384.    REPETER BQTLD nbgau ;
  1385.       igau = igau + 1 ;
  1386.       SI (EGA qtldu1.igau 0.0) ;
  1387.          SI (NON (EGA qtldu2.igau 0.0)) ;
  1388. ************SI (EGA IMESS 1) ;
  1389. ************Message 'Erreur absolue' qtldu2.igau
  1390. ************        'Composante ' (EXTR LVARI icomp)
  1391. ************        ' Point Gauss' igau 'Etat' ietat ;
  1392. ************Message ;
  1393. ************FINSI ;
  1394.             ierr = 50 ; errabs = qtldu2.igau ; QUITTER BECART ;
  1395.          FINSI ;
  1396.       SINON ;
  1397.          dqtld = ABS (qtldu1.igau - qtldu2.igau) ;
  1398.          errrel = ABS (dqtld / qtldu1.igau) ;
  1399.          SI (errrel > critere) ;
  1400. ************SI (EGA IMESS 1) ;
  1401. ************Message 'Erreur relative' errrel
  1402. ************        'Composante ' (EXTR LVARI icomp)
  1403. ************        ' Point Gauss' igau 'Etat' ietat ;
  1404. ************Message 'qtldu1.igau =' qtldu1.igau
  1405. ************        'qtldu2.igau =' qtldu2.igau;
  1406. ************Message ;
  1407. ************FINSI ;
  1408.             ierr = 51 ; QUITTER BECART ;
  1409.          FINSI ;
  1410.       FINSI ;
  1411.    FIN BQTLD ;
  1412.  
  1413. *  Comparaison des variables internes supplementaires (modele no23)-----
  1414.  
  1415.    icomp = 1 ;
  1416.    igau = 0 ;
  1417.    REPETER BPSUP nbgau ;
  1418.       igau = igau + 1 ;
  1419.       SI (EGA pu2.igau 0.0) ;
  1420.          SI (NON (EGA psupu2.igau 0.0)) ;
  1421. ************SI (EGA IMESS 1) ;
  1422. ************Message 'Erreur absolue' psupu2.igau
  1423. ************        'Composante ' (EXTR LVSUP icomp)
  1424. ************        ' Point Gauss' igau 'Etat' ietat ;
  1425. ************Message ;
  1426. ************FINSI ;
  1427.             ierr = 60 ; errabs = psupu2.igau ; QUITTER BECART ;
  1428.          FINSI ;
  1429.       SINON ;
  1430.          dpsup = ABS (pu2.igau - psupu2.igau) ;
  1431.          errrel = ABS (dpsup / pu2.igau) ;
  1432.          SI (errrel > critere) ;
  1433. ************SI (EGA IMESS 1) ;
  1434. ************Message 'Erreur relative' errrel
  1435. ************        'Composante ' (EXTR LVSUP icomp)
  1436. ************        ' Point Gauss' igau 'Etat' ietat ;
  1437. ************Message 'pu2.igau =' pu2.igau
  1438. ************        'psupu2.igau =' psupu2.igau;
  1439. ************Message ;
  1440. ************FINSI ;
  1441.             ierr = 61 ; QUITTER BECART ;
  1442.          FINSI ;
  1443.       FINSI ;
  1444.    FIN BPSUP ;
  1445.  
  1446.    icomp = 2 ;
  1447.    igau = 0 ;
  1448.    REPETER BQSUP nbgau ;
  1449.       igau = igau + 1 ;
  1450.       SI (EGA qtldu2.igau 0.0) ;
  1451.          SI (NON (EGA qsupu2.igau 0.0)) ;
  1452. ************SI (EGA IMESS 1) ;
  1453. ************Message 'Erreur absolue' qsupu2.igau
  1454. ************        'Composante ' (EXTR LVSUP icomp)
  1455. ************        ' Point Gauss' igau 'Etat' ietat ;
  1456. ************Message ;
  1457. ************FINSI ;
  1458.             ierr = 60 ; errabs = qsupu2.igau ; QUITTER BECART ;
  1459.          FINSI ;
  1460.       SINON ;
  1461.          dqsup = ABS (qtldu2.igau - qsupu2.igau) ;
  1462.          errrel = ABS (dqsup / qtldu2.igau) ;
  1463.          SI (errrel > critere) ;
  1464. ************SI (EGA IMESS 1) ;
  1465. ************Message 'Erreur relative' errrel
  1466. ************        'Composante ' (EXTR LVSUP icomp)
  1467. ************        ' Point Gauss' igau 'Etat' ietat ;
  1468. ************Message 'qtldu2.igau =' qtldu2.igau
  1469. ************        'qsupu2.igau =' qsupu2.igau;
  1470. ************Message ;
  1471. ************FINSI ;
  1472.             ierr = 61 ; QUITTER BECART ;
  1473.          FINSI ;
  1474.       FINSI ;
  1475.    FIN BQSUP ;
  1476.  
  1477. FIN BECART ;
  1478.  
  1479. * FIN   BOUCLE COMPARAISON DES RESULTATS________________________________
  1480.  
  1481. * MESSAGES D'ERREUR_____________________________________________________
  1482.  
  1483. Message ;
  1484. Message 'CPU avec modele externe 21 = ' CPUuti1 'centiemes de seconde' ;
  1485. Message 'CPU avec modele externe 23 = ' CPUuti2 'centiemes de seconde' ;
  1486. Message ;
  1487.  
  1488. ifin = netat-1 ;
  1489.  
  1490. gaxyu1ma = 0.0 ;
  1491. igau = 0 ;
  1492. REPETER BGAXYFU1 nbgau ;
  1493.    igau = igau + 1 ;
  1494.    gaxy = EXTR (tabuti1.'DEFORMATIONS'.ifin)
  1495.                     'GAXY' 1 1 igau ;
  1496.    SI (gaxy > gaxyu1ma) ; gaxyu1ma = gaxy ; FINSI ;
  1497. FIN BGAXYFU1 ;
  1498. gaxyu2ma = 0.0 ;
  1499. igau = 0 ;
  1500. REPETER BGAXYFU2 nbgau ;
  1501.    igau = igau + 1 ;
  1502.    gaxy = EXTR (tabuti2.'DEFORMATIONS'.ifin)
  1503.                     'GAXY' 1 1 igau ;
  1504.    SI (gaxy > gaxyu2ma) ; gaxyu2ma = gaxy ; FINSI ;
  1505. FIN BGAXYFU2 ;
  1506. Message 'Deformation totale finale (maxi sur points de Gauss)' ;
  1507. Message 'GAXY max avec modele externe 21 = ' gaxyu1ma ;
  1508. Message 'GAXY max avec modele externe 23 = ' gaxyu2ma ;
  1509. Message ;
  1510.  
  1511. gixyu1ma = 0.0 ;
  1512. igau = 0 ;
  1513. REPETER BGIXYFU1 nbgau ;
  1514.    igau = igau + 1 ;
  1515.    gixy = EXTR (tabuti1.'DEFORMATIONS_INELASTIQUES'.ifin)
  1516.                     'GIXY' 1 1 igau ;
  1517.    SI (gixy > gixyu1ma) ; gixyu1ma = gixy ; FINSI ;
  1518. FIN BGIXYFU1 ;
  1519. gixyu2ma = 0.0 ;
  1520. igau = 0 ;
  1521. REPETER BGIXYFU2 nbgau ;
  1522.    igau = igau + 1 ;
  1523.    gixy = EXTR (tabuti2.'DEFORMATIONS_INELASTIQUES'.ifin)
  1524.                     'GIXY' 1 1 igau ;
  1525.    SI (gixy > gixyu2ma) ; gixyu2ma = gixy ; FINSI ;
  1526. FIN BGIXYFU2 ;
  1527. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  1528. Message 'GIXY max avec modele externe 21 = ' gixyu1ma ;
  1529. Message 'GIXY max avec modele externe 23 = ' gixyu2ma ;
  1530. Message ;
  1531.  
  1532. smxyu1ma = 0.0 ;
  1533. igau = 0 ;
  1534. REPETER BSMXYFU1 nbgau ;
  1535.    igau = igau + 1 ;
  1536.    smxy = EXTR (tabuti1.'CONTRAINTES'.ifin)
  1537.                     'SMXY' 1 1 igau ;
  1538.    SI (smxy > smxyu1ma) ; smxyu1ma = smxy ; FINSI ;
  1539. FIN BSMXYFU1 ;
  1540. smxyu2ma = 0.0 ;
  1541. igau = 0 ;
  1542. REPETER BSMXYFU2 nbgau ;
  1543.    igau = igau + 1 ;
  1544.    smxy = EXTR (tabuti2.'CONTRAINTES'.ifin)
  1545.                     'SMXY' 1 1 igau ;
  1546.    SI (smxy > smxyu2ma) ; smxyu2ma = smxy ; FINSI ;
  1547. FIN BSMXYFU2 ;
  1548. Message 'Contrainte finale (maxi sur points de Gauss)' ;
  1549. Message 'SMXY max avec modele externe 21 = ' smxyu1ma ;
  1550. Message 'SMXY max avec modele externe 23 = ' smxyu2ma ;
  1551. Message ;
  1552.  
  1553. SI (ierr > 0) ;
  1554.  
  1555.    SI (EGA ierr 1) ;
  1556.       Message 'Instants de calcul non coincidants' ;
  1557.       ERRE 5 ;
  1558.    FINSI ;
  1559.  
  1560.    SI (EGA ierr 10) ;
  1561.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1562.               ' Noeud' ino 'Etat' ietat ;
  1563.       ERRE 5 ;
  1564.    FINSI ;
  1565.  
  1566.    SI (EGA ierr 11) ;
  1567.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1568.               ' Noeud' ino 'Etat' ietat ;
  1569.       ERRE 5 ;
  1570.    FINSI ;
  1571.  
  1572.    SI (EGA ierr 20) ;
  1573.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  1574.               ' Point Gauss' igau 'Etat' ietat ;
  1575.       ERRE 5 ;
  1576.    FINSI ;
  1577.  
  1578.    SI (EGA ierr 21) ;
  1579.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  1580.               ' Point Gauss' igau 'Etat' ietat ;
  1581.       ERRE 5 ;
  1582.    FINSI ;
  1583.  
  1584.    SI (EGA ierr 30) ;
  1585.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  1586.               ' Point Gauss' igau 'Etat' ietat ;
  1587.       ERRE 5 ;
  1588.    FINSI ;
  1589.  
  1590.    SI (EGA ierr 31) ;
  1591.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  1592.               ' Point Gauss' igau 'Etat' ietat ;
  1593.       ERRE 5 ;
  1594.    FINSI ;
  1595.  
  1596.    SI (EGA ierr 40) ;
  1597.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  1598.               ' Point Gauss' igau 'Etat' ietat ;
  1599.       ERRE 5 ;
  1600.    FINSI ;
  1601.  
  1602.    SI (EGA ierr 41) ;
  1603.       Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  1604.               ' Point Gauss' igau 'Etat' ietat ;
  1605.       ERRE 5 ;
  1606.    FINSI ;
  1607.  
  1608.    SI (EGA ierr 50) ;
  1609.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  1610.               ' Point Gauss' igau 'Etat' ietat ;
  1611.       ERRE 5 ;
  1612.    FINSI ;
  1613.  
  1614.    SI (EGA ierr 51) ;
  1615.       Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  1616.               ' Point Gauss' igau 'Etat' ietat ;
  1617.       ERRE 5 ;
  1618.    FINSI ;
  1619.  
  1620.    SI (EGA ierr 60) ;
  1621.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LVSUP icomp)
  1622.               ' Point Gauss' igau 'Etat' ietat ;
  1623.       ERRE 5 ;
  1624.    FINSI ;
  1625.  
  1626.    SI (EGA ierr 61) ;
  1627.       Message 'Erreur relative' errrel 'Composante ' (EXTR LVSUP icomp)
  1628.               ' Point Gauss' igau 'Etat' ietat ;
  1629.       ERRE 5 ;
  1630.    FINSI ;
  1631.  
  1632. SINON ;
  1633.  
  1634.    Message 'Resultats concordants entre les deux modeles' ;
  1635.    Message 'Aucune erreur relative superieure a ' critere
  1636.            ' pour les composantes' ;
  1637.    Message (EXTR LDEFO idefo) '  ' (EXTR LDEIN idein) '  '
  1638.            (EXTR LCONT icont) '  '
  1639.            (EXTR LVARI ivar1) '  ' (EXTR LVARI ivar2) '  '
  1640.            (EXTR LVARI ivar3) '  ' (EXTR LVARI ivar4) ;
  1641.    Message ;
  1642.    Message 'Modele externe 23 : aucune erreur relative superieure a '
  1643.             critere ;
  1644.    Message 'entre les composantes  '
  1645.            (EXTR LVARI ivar3) '  ' (EXTR LVARI ivar4) '  et  '
  1646.            (EXTR LVSUP ivas1) '  ' (EXTR LVSUP ivas2) ;
  1647.    ERRE 0 ;
  1648.  
  1649. FINSI ;
  1650.  
  1651. FIN ;
  1652.  
  1653.  
  1654.  
  1655.  
  1656.  
  1657.  
  1658.  
  1659.  
  1660.  
  1661.  
  1662.  
  1663.  

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