Télécharger la plaquette Cast3M


fr en

Télécharger creep04_cisXY.dgibi

Retour à la liste

Numérotation des lignes : 

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

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