Télécharger la plaquette Cast3M


fr en

Télécharger creep01_traZZ.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  creep01_traZZ.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 No1 : Modele de FLUAGE POLYNOMIAL
  15. *            TRACTION ZZ
  16. *
  17. *=======================================================================
  18.  
  19.  
  20. opti echo 0 ;
  21.  
  22.  
  23. *=======================================================================
  24. *
  25. *                       PROCEDURES DEPENDANTES
  26. *
  27. *=======================================================================
  28.  
  29.  
  30.  
  31.  
  32. *=======================================================================
  33. *
  34. *                      JEU DE DONNEES PRINCIPAL
  35. *
  36. *=======================================================================
  37.  
  38. *-----------------------------------------------------------------------
  39. * 0 - Maillage
  40. *     PLAQUE CARREE    COTE : 1 m    EPAISSEUR : 1 cm
  41.  
  42. opti dime 3 elem cub8 mode tridim ;
  43.  
  44. l1 = 1.0 ;
  45. l2 = 1.0 ;
  46. ep0 = 0.01 ;
  47.  
  48. p01 = 0.0 0.0 0.0 ;
  49. p02 = ep0 0.0 0.0 ;
  50. p03 = ep0 l2  0.0 ;
  51. p04 = 0.0 l2  0.0 ;
  52.  
  53. p10 = 0.0 0.0 l1 ;
  54. p20 = ep0 0.0 l1 ;
  55. p30 = ep0 l2  l1 ;
  56. p40 = 0.0 l2  l1 ;
  57.  
  58. li01 = 'DROIT' 1 p01 p02 ;
  59. li02 = 'DROIT' 1 p02 p03 ;
  60. li03 = 'DROIT' 1 p03 p04 ;
  61. li04 = 'DROIT' 1 p04 p01 ;
  62.  
  63. li10 = 'DROIT' 1 p10 p20 ;
  64. li20 = 'DROIT' 1 p20 p30 ;
  65. li30 = 'DROIT' 1 p30 p40 ;
  66. li40 = 'DROIT' 1 p40 p10 ;
  67.  
  68. sur01 = 'DALL' li01 li02 li03 li04 ;
  69. sur10 = 'DALL' li10 li20 li30 li40 ;
  70. mail1 = sur01 'VOLU' 1 sur10 ;
  71.  
  72. *-----------------------------------------------------------------------
  73. * 1 - Calcul avec le modele 'FLUAGE' 'POLYNOMIAL' interne CAST3M
  74.  
  75. * 1.1 Definition du modele
  76.  
  77. mo_flupo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  78.                       'FLUAGE' 'POLYNOMIAL' ;
  79.  
  80. * 1.2 Affectation des proprietes de materiau
  81.  
  82. densiu   = 0.95 ;
  83. dgrain   = 8.0E-6 ;
  84. tuo2     = 1500.0 ;
  85. tauxfisk = 1.5E19 ;
  86.  
  87. xyoun = (2.2693E11 - (1.5399E7*tuo2)) - (9.35973E3 * (tuo2**2.0)) ;
  88. xcis  = (8.583E10 - (5.157E6*tuo2)) - (3.747E3 * (tuo2**2.0)) ;
  89. xnu   = (xyoun/(2.0*xcis)) - 1.0 ;
  90. xrho  = 10950.0 ;
  91. xalph = ( (1.1833E-5 - (5.013E-9*tuo2)) + (3.756E-12 * (tuo2**2.0)) )
  92.         - (6.125E-17 * (tuo2**3.0)) ;
  93.  
  94. xsmax = xyoun / 1000.0 ;
  95. xaf0  = 0.0 ;
  96. xaf1  = (EXP(0.183*(100.0-densiu)))
  97.       * ( (0.2031E-13+(0.67861E-32*tauxfisk)) / (dgrain**2.0) )
  98.       * (EXP(-45294.4*(tuo2**-1.0))) ;
  99. xaf2  = 1.0 ;
  100. xaf3  = (EXP(0.241*(100.0-densiu)))
  101.       * (0.1524E-25+(0.50825E-44*tauxfisk))
  102.       * (EXP(-66431.8*(tuo2**-1.0))) ;
  103. xaf4  = 4.5 ;
  104. xaf5  = (1.447E-34*tauxfisk) * (EXP(-3624.0*(tuo2**-1.0))) ;
  105. xaf6  = 1.0 ;
  106.  
  107. ma_flupo = MATE mo_flupo 'YOUN' xyoun 'NU  ' xnu   'RHO ' xrho
  108.                          'ALPH' xalph
  109.                          'SMAX' xsmax 'AF0 ' xaf0  'AF1 ' xaf1
  110.                          'AF2 ' xaf2  'AF3 ' xaf3  'AF4 ' xaf4
  111.                          'AF5 ' xaf5  'AF6 ' xaf6 ;
  112.  
  113. * 1.3 Definition des conditions aux limites et chargements
  114.  
  115. cl1 = BLOQ mail1 'UX' 'UY' ;
  116. cl2 = BLOQ sur01 'UZ' ;
  117. cl0 = cl1 ET cl2 ;
  118.  
  119. depmax = 5.0E8 ;
  120. dep0 = FORCE sur10 'FX' 0.0 'FY' 0.0 'FZ' depmax ;
  121. ev0 = EVOL MANU temps (PROG 0.0 1.E-2 1000.0)
  122.                     y (PROG 0.0 1.0   1.0   ) ;
  123. cha0 = CHAR MECA dep0 ev0 ;
  124.  
  125. * 1.4 Definition des instants de calcul
  126.  
  127. pro0 = PROG 0.0 PAS 2.E-4 6.E-3 ;
  128. pro1 = PROG 0.0 PAS 2.E-4 6.E-3 ;
  129.  
  130. * 1.5 Resolution par PASAPAS
  131.  
  132. tabflup = TABLE ;
  133. tabflup.'MODELE' = mo_flupo ;
  134. tabflup.'CARACTERISTIQUES' = ma_flupo ;
  135. tabflup.'CHARGEMENT' = cha0 ;
  136. tabflup.'BLOCAGES_MECANIQUES' = cl0 ;
  137. tabflup.'TEMPS_CALCULES' = pro0 ;
  138. tabflup.'TEMPS_SAUVES' = pro1 ;
  139. TMASAU=table;
  140. tabflup . 'MES_SAUVEGARDES'=TMASAU;
  141. TMASAU .'DEFTO'=VRAI;
  142. TMASAU .'DEFIN'=VRAI;
  143. TEMPS 'ZERO' ;
  144. PASAPAS tabflup ;
  145. TABTPS = TEMP 'NOEC';
  146. CPUint = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  147.  
  148. *-----------------------------------------------------------------------
  149. * 2 - Calcul avec modele 'VISCO_EXTERNE' 'GENERAL' no21
  150. *     Modele 'FLUAGE' 'POLYNOMIAL' externalise dans CREEP
  151.  
  152. * 2.1 Definition du modele utilisateur externalise dans CREEP
  153.  
  154. mo_util = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  155.                      'VISCO_EXTERNE' 'GENERAL' 'NUME_LOI' 21 ;
  156.  
  157. * 2.2 Affectation des proprietes de materiau
  158.  
  159. ma_util = MATE mo_util 'YOUN' xyoun 'NU  ' xnu   'RHO ' xrho
  160.                        'ALPH' xalph ;
  161.  
  162. * 2.3 Definition des conditions aux limites et chargements : idem 1.3
  163.  
  164. * 2.4 Definition des instants de calcul : idem 1.4
  165.  
  166. * 2.5 Resolution par PASAPAS
  167.  
  168. tabutil = TABLE ;
  169. tabutil.'MODELE' = mo_util ;
  170. tabutil.'CARACTERISTIQUES' = ma_util ;
  171. tabutil.'CHARGEMENT' = cha0 ;
  172. tabutil.'BLOCAGES_MECANIQUES' = cl0 ;
  173. tabutil.'TEMPS_CALCULES' = pro0 ;
  174. tabutil.'TEMPS_SAUVES' = pro1 ;
  175. TMASAU=table;
  176. tabutil . 'MES_SAUVEGARDES'=TMASAU;
  177. TMASAU .'DEFTO'=VRAI;
  178. TMASAU .'DEFIN'=VRAI;
  179. TEMPS 'ZERO' ;
  180. PASAPAS tabutil ;
  181. TABTPS = TEMP 'NOEC';
  182. CPUext = TABTPS.'TEMPS_CPU'.'INITIAL' ;
  183.  
  184. *-----------------------------------------------------------------------
  185. * 3 - Post-traitement : comparaisons entre les deux modeles
  186. *
  187. *     N.B. Comme on simule un essai de traction ZZ :
  188. *          - Critere strict (errrel < 1.d-6) sur EPZZ, EIZZ et SMZZ
  189. *          - Pour les autres composantes on note les depassements du
  190. *            critere si impression demandee
  191.  
  192. uxc = TABLE ; uyc = TABLE ; uzc = TABLE ;
  193. epxxc = TABLE ; epyyc = TABLE ; epzzc = TABLE ;
  194. gaxyc = TABLE ; gaxzc = TABLE ; gayzc = TABLE ;
  195. eixxc = TABLE ; eiyyc = TABLE ; eizzc = TABLE ;
  196. gixyc = TABLE ; gixzc = TABLE ; giyzc = TABLE ;
  197. smxxc = TABLE ; smyyc = TABLE ; smzzc = TABLE ;
  198. smxyc = TABLE ; smxzc = TABLE ; smyzc = TABLE ;
  199.  
  200. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  201. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  202. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  203. eixxu = TABLE ; eiyyu = TABLE ; eizzu = TABLE ;
  204. gixyu = TABLE ; gixzu = TABLE ; giyzu = TABLE ;
  205. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  206. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  207.  
  208. LDEPL = MOTS 'UX  ' 'UY  ' 'UZ  ' ;
  209. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  210. LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  211. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  212.  
  213. idefo = 3 ;
  214. idein = 3 ;
  215. icont = 3 ;
  216.  
  217. ierr = 0 ;
  218. critere = 1.d-6 ;
  219. IMESS = 0 ;
  220. ***IMESS = 1 ;
  221.  
  222. nbno  = 8 ;
  223. nbgau = 8 ;
  224.  
  225. netat = DIME (tabflup.'TEMPS') ;
  226. ietat = -1 ;
  227.  
  228.  
  229. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  230.  
  231. REPETER BECART netat ;
  232.  
  233.    ietat = ietat + 1 ;
  234.  
  235.    temps1 = tabflup.'TEMPS'.ietat ;
  236.    temps2 = tabutil.'TEMPS'.ietat ;
  237.    SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  238.       ierr = 1 ;
  239.       QUITTER BECART ;
  240.    FINSI ;
  241.  
  242. *  Extraction des resultats de calcul : Modele interne CAST3M-----------
  243.  
  244. *  Deplacements
  245.    uxc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  246.    uxc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  247.    uxc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  248.    uxc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  249.    uxc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  250.    uxc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  251.    uxc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  252.    uxc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  253.  
  254.    uyc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  255.    uyc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  256.    uyc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  257.    uyc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  258.    uyc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  259.    uyc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  260.    uyc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  261.    uyc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  262.  
  263.    uzc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  264.    uzc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  265.    uzc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  266.    uzc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  267.    uzc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  268.    uzc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  269.    uzc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  270.    uzc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  271.  
  272. *  Deformations totales
  273.    igau = 0 ;
  274.    REPETER BEPXXC nbgau ;
  275.       igau = igau + 1 ;
  276.       epxxc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  277.    FIN BEPXXC ;
  278.  
  279.    igau = 0 ;
  280.    REPETER BEPYYC nbgau ;
  281.       igau = igau + 1 ;
  282.       epyyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  283.    FIN BEPYYC ;
  284.  
  285.    igau = 0 ;
  286.    REPETER BEPZZC nbgau ;
  287.       igau = igau + 1 ;
  288.       epzzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  289.    FIN BEPZZC ;
  290.  
  291.    igau = 0 ;
  292.    REPETER BGAXYC nbgau ;
  293.       igau = igau + 1 ;
  294.       gaxyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  295.    FIN BGAXYC ;
  296.  
  297.    igau = 0 ;
  298.    REPETER BGAXZC nbgau ;
  299.       igau = igau + 1 ;
  300.       gaxzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  301.    FIN BGAXZC ;
  302.  
  303.    igau = 0 ;
  304.    REPETER BGAYZC nbgau ;
  305.       igau = igau + 1 ;
  306.       gayzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  307.    FIN BGAYZC ;
  308.  
  309. *  Deformations inelastiques
  310.    igau = 0 ;
  311.    REPETER BEIXXC nbgau ;
  312.       igau = igau + 1 ;
  313.       eixxc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  314.                              'EIXX' 1 1 igau ;
  315.    FIN BEIXXC ;
  316.  
  317.    igau = 0 ;
  318.    REPETER BEIYYC nbgau ;
  319.       igau = igau + 1 ;
  320.       eiyyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  321.                              'EIYY' 1 1 igau ;
  322.    FIN BEIYYC ;
  323.  
  324.    igau = 0 ;
  325.    REPETER BEIZZC nbgau ;
  326.       igau = igau + 1 ;
  327.       eizzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  328.                              'EIZZ' 1 1 igau ;
  329.    FIN BEIZZC ;
  330.  
  331.    igau = 0 ;
  332.    REPETER BGIXYC nbgau ;
  333.       igau = igau + 1 ;
  334.       gixyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  335.                              'GIXY' 1 1 igau ;
  336.    FIN BGIXYC ;
  337.  
  338.    igau = 0 ;
  339.    REPETER BGIXZC nbgau ;
  340.       igau = igau + 1 ;
  341.       gixzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  342.                              'GIXZ' 1 1 igau ;
  343.    FIN BGIXZC ;
  344.  
  345.    igau = 0 ;
  346.    REPETER BGIYZC nbgau ;
  347.       igau = igau + 1 ;
  348.       giyzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  349.                              'GIYZ' 1 1 igau ;
  350.    FIN BGIYZC ;
  351.  
  352. *  Contraintes
  353.    igau = 0 ;
  354.    REPETER BSMXXC nbgau ;
  355.       igau = igau + 1 ;
  356.       smxxc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  357.    FIN BSMXXC ;
  358.  
  359.    igau = 0 ;
  360.    REPETER BSMYYC nbgau ;
  361.       igau = igau + 1 ;
  362.       smyyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  363.    FIN BSMYYC ;
  364.  
  365.    igau = 0 ;
  366.    REPETER BSMZZC nbgau ;
  367.       igau = igau + 1 ;
  368.       smzzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  369.    FIN BSMZZC ;
  370.  
  371.    igau = 0 ;
  372.    REPETER BSMXYC nbgau ;
  373.       igau = igau + 1 ;
  374.       smxyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  375.    FIN BSMXYC ;
  376.  
  377.    igau = 0 ;
  378.    REPETER BSMXZC nbgau ;
  379.       igau = igau + 1 ;
  380.       smxzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  381.    FIN BSMXZC ;
  382.  
  383.    igau = 0 ;
  384.    REPETER BSMYZC nbgau ;
  385.       igau = igau + 1 ;
  386.       smyzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  387.    FIN BSMYZC ;
  388.  
  389. *  Extraction des resultats de calcul : Modele externe utilisateur------
  390.  
  391. *  Deplacements
  392.    uxu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  393.    uxu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  394.    uxu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  395.    uxu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  396.    uxu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  397.    uxu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  398.    uxu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  399.    uxu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  400.  
  401.    uyu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  402.    uyu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  403.    uyu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  404.    uyu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  405.    uyu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  406.    uyu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  407.    uyu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  408.    uyu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  409.  
  410.    uzu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  411.    uzu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  412.    uzu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  413.    uzu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  414.    uzu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  415.    uzu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  416.    uzu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  417.    uzu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  418.  
  419. *  Deformations totales
  420.    igau = 0 ;
  421.    REPETER BEPXXU nbgau ;
  422.       igau = igau + 1 ;
  423.       epxxu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  424.    FIN BEPXXU ;
  425.  
  426.    igau = 0 ;
  427.    REPETER BEPYYU nbgau ;
  428.       igau = igau + 1 ;
  429.       epyyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  430.    FIN BEPYYU ;
  431.  
  432.    igau = 0 ;
  433.    REPETER BEPZZU nbgau ;
  434.       igau = igau + 1 ;
  435.       epzzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  436.    FIN BEPZZU ;
  437.  
  438.    igau = 0 ;
  439.    REPETER BGAXYU nbgau ;
  440.       igau = igau + 1 ;
  441.       gaxyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  442.    FIN BGAXYU ;
  443.  
  444.    igau = 0 ;
  445.    REPETER BGAXZU nbgau ;
  446.       igau = igau + 1 ;
  447.       gaxzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  448.    FIN BGAXZU ;
  449.  
  450.    igau = 0 ;
  451.    REPETER BGAYZU nbgau ;
  452.       igau = igau + 1 ;
  453.       gayzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  454.    FIN BGAYZU ;
  455.  
  456. *  Deformations inelastiques
  457.    igau = 0 ;
  458.    REPETER BEIXXU nbgau ;
  459.       igau = igau + 1 ;
  460.       eixxu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  461.                               'EIXX' 1 1 igau ;
  462.    FIN BEIXXU ;
  463.  
  464.    igau = 0 ;
  465.    REPETER BEIYYU nbgau ;
  466.       igau = igau + 1 ;
  467.       eiyyu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  468.                               'EIYY' 1 1 igau ;
  469.    FIN BEIYYU ;
  470.  
  471.    igau = 0 ;
  472.    REPETER BEIZZU nbgau ;
  473.       igau = igau + 1 ;
  474.       eizzu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  475.                               'EIZZ' 1 1 igau ;
  476.    FIN BEIZZU ;
  477.  
  478.    igau = 0 ;
  479.    REPETER BGIXYU nbgau ;
  480.       igau = igau + 1 ;
  481.       gixyu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  482.                               'GIXY' 1 1 igau ;
  483.    FIN BGIXYU ;
  484.  
  485.    igau = 0 ;
  486.    REPETER BGIXZU nbgau ;
  487.       igau = igau + 1 ;
  488.       gixzu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  489.                               'GIXZ' 1 1 igau ;
  490.    FIN BGIXZU ;
  491.  
  492.    igau = 0 ;
  493.    REPETER BGIYZU nbgau ;
  494.       igau = igau + 1 ;
  495.       giyzu.igau = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ietat)
  496.                               'GIYZ' 1 1 igau ;
  497.    FIN BGIYZU ;
  498.  
  499. *  Contraintes
  500.    igau = 0 ;
  501.    REPETER BSMXXU nbgau ;
  502.       igau = igau + 1 ;
  503.       smxxu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  504.    FIN BSMXXU ;
  505.  
  506.    igau = 0 ;
  507.    REPETER BSMYYU nbgau ;
  508.       igau = igau + 1 ;
  509.       smyyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  510.    FIN BSMYYU ;
  511.  
  512.    igau = 0 ;
  513.    REPETER BSMZZU nbgau ;
  514.       igau = igau + 1 ;
  515.       smzzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  516.    FIN BSMZZU ;
  517.  
  518.    igau = 0 ;
  519.    REPETER BSMXYU nbgau ;
  520.       igau = igau + 1 ;
  521.       smxyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  522.    FIN BSMXYU ;
  523.  
  524.    igau = 0 ;
  525.    REPETER BSMXZU nbgau ;
  526.       igau = igau + 1 ;
  527.       smxzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  528.    FIN BSMXZU ;
  529.  
  530.    igau = 0 ;
  531.    REPETER BSMYZU nbgau ;
  532.       igau = igau + 1 ;
  533.       smyzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  534.    FIN BSMYZU ;
  535.  
  536. *  Comparaison des deplacements-----------------------------------------
  537.  
  538.    icomp = 1 ;
  539.    ino = 0 ;
  540.    REPETER BUX nbno ;
  541.       ino = ino + 1 ;
  542.       SI (EGA uxc.ino 0.0) ;
  543.          SI (NON (EGA uxu.ino 0.0)) ;
  544.             SI (EGA IMESS 1) ;
  545.             Message 'Erreur absolue' uxu.ino
  546.                     'Composante ' (EXTR LDEPL icomp)
  547.                     ' Noeud' ino 'Etat' ietat ;
  548.             Message ;
  549.             FINSI ;
  550. ************ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  551.          FINSI ;
  552.       SINON ;
  553.          dux = ABS (uxc.ino - uxu.ino) ;
  554.          errrel = ABS (dux / uxc.ino) ;
  555.          SI (errrel > critere) ;
  556.             SI (EGA IMESS 1) ;
  557.             Message 'Erreur relative' errrel
  558.                     'Composante ' (EXTR LDEPL icomp)
  559.                     ' Noeud' ino 'Etat' ietat ;
  560.             Message 'uxc.ino =' uxc.ino 'uxu.ino =' uxu.ino ;
  561.             Message ;
  562.             FINSI ;
  563. ************ierr = 11 ; QUITTER BECART ;
  564.          FINSI ;
  565.       FINSI ;
  566.    FIN BUX ;
  567.  
  568.    icomp = 2 ;
  569.    ino = 0 ;
  570.    REPETER BUY nbno ;
  571.       ino = ino + 1 ;
  572.       SI (EGA uyc.ino 0.0) ;
  573.          SI (NON (EGA uyu.ino 0.0)) ;
  574.             SI (EGA IMESS 1) ;
  575.             Message 'Erreur absolue' uyu.ino
  576.                     'Composante ' (EXTR LDEPL icomp)
  577.                     ' Noeud' ino 'Etat' ietat ;
  578.             Message ;
  579.             FINSI ;
  580. ************ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  581.          FINSI ;
  582.       SINON ;
  583.          duy = ABS (uyc.ino - uyu.ino) ;
  584.          errrel = ABS (duy / uyc.ino) ;
  585.          SI (errrel > critere) ;
  586.             SI (EGA IMESS 1) ;
  587.             Message 'Erreur relative' errrel
  588.                     'Composante ' (EXTR LDEPL icomp)
  589.                     ' Noeud' ino 'Etat' ietat ;
  590.             Message 'uyc.ino =' uyc.ino 'uyu.ino =' uyu.ino ;
  591.             Message ;
  592.             FINSI ;
  593. ************ierr = 11 ; QUITTER BECART ;
  594.          FINSI ;
  595.       FINSI ;
  596.    FIN BUY ;
  597.  
  598.    icomp = 3 ;
  599.    ino = 0 ;
  600.    REPETER BUZ nbno ;
  601.       ino = ino + 1 ;
  602.       SI (EGA uzc.ino 0.0) ;
  603.          SI (NON (EGA uzu.ino 0.0)) ;
  604.             SI (EGA IMESS 1) ;
  605.             Message 'Erreur absolue' uzu.ino
  606.                     'Composante ' (EXTR LDEPL icomp)
  607.                     ' Noeud' ino 'Etat' ietat ;
  608.             Message ;
  609.             FINSI ;
  610. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  611.          FINSI ;
  612.       SINON ;
  613.          duz = ABS (uzc.ino - uzu.ino) ;
  614.          errrel = ABS (duz / uzc.ino) ;
  615.          SI (errrel > critere) ;
  616.             SI (EGA IMESS 1) ;
  617.             Message 'Erreur relative' errrel
  618.                     'Composante ' (EXTR LDEPL icomp)
  619.                     ' Noeud' ino 'Etat' ietat ;
  620.             Message 'uzc.ino =' uzc.ino 'uzu.ino =' uzu.ino ;
  621.             Message ;
  622.             FINSI ;
  623. ************ierr = 11 ; QUITTER BECART ;
  624.          FINSI ;
  625.       FINSI ;
  626.    FIN BUZ ;
  627.  
  628. *  Comparaison des deformations-----------------------------------------
  629.  
  630.    icomp = 1 ;
  631.    igau = 0 ;
  632.    REPETER BEPXX nbgau ;
  633.       igau = igau + 1 ;
  634.       SI (EGA epxxc.igau 0.0) ;
  635.          SI (NON (EGA epxxu.igau 0.0)) ;
  636.             SI (EGA IMESS 1) ;
  637.             Message 'Erreur absolue' epxxu.igau
  638.                     'Composante ' (EXTR LDEFO icomp)
  639.                     ' Point Gauss' igau 'Etat' ietat ;
  640.             Message ;
  641.             FINSI ;
  642. ************ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  643.          FINSI ;
  644.       SINON ;
  645.          depxx = ABS (epxxc.igau - epxxu.igau) ;
  646.          errrel = ABS (depxx / epxxc.igau) ;
  647.          SI (errrel > critere) ;
  648.             SI (EGA IMESS 1) ;
  649.             Message 'Erreur relative' errrel
  650.                     'Composante ' (EXTR LDEFO icomp)
  651.                     ' Point Gauss' igau 'Etat' ietat ;
  652.             Message 'epxxc.igau =' epxxc.igau 'epxxu.igau =' epxxu.igau;
  653.             Message ;
  654.             FINSI ;
  655. ************ierr = 21 ; QUITTER BECART ;
  656.          FINSI ;
  657.       FINSI ;
  658.    FIN BEPXX ;
  659.  
  660.    icomp = 2 ;
  661.    igau = 0 ;
  662.    REPETER BEPYY nbgau ;
  663.       igau = igau + 1 ;
  664.       SI (EGA epyyc.igau 0.0) ;
  665.          SI (NON (EGA epyyu.igau 0.0)) ;
  666.             SI (EGA IMESS 1) ;
  667.             Message 'Erreur absolue' epyyu.igau
  668.                     'Composante ' (EXTR LDEFO icomp)
  669.                     ' Point Gauss' igau 'Etat' ietat ;
  670.             Message ;
  671.             FINSI ;
  672. ************ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  673.          FINSI ;
  674.       SINON ;
  675.          depyy = ABS (epyyc.igau - epyyu.igau) ;
  676.          errrel = ABS (depyy / epyyc.igau) ;
  677.          SI (errrel > critere) ;
  678.             SI (EGA IMESS 1) ;
  679.             Message 'Erreur relative' errrel
  680.                     'Composante ' (EXTR LDEFO icomp)
  681.                     ' Point Gauss' igau 'Etat' ietat ;
  682.             Message 'epyyc.igau =' epyyc.igau 'epyyu.igau =' epyyu.igau;
  683.             Message ;
  684.             FINSI ;
  685. ************ierr = 21 ; QUITTER BECART ;
  686.          FINSI ;
  687.       FINSI ;
  688.    FIN BEPYY ;
  689.  
  690.    icomp = 3 ;
  691.    igau = 0 ;
  692.    REPETER BEPZZ nbgau ;
  693.       igau = igau + 1 ;
  694.       SI (EGA epzzc.igau 0.0) ;
  695.          SI (NON (EGA epzzu.igau 0.0)) ;
  696. ************SI (EGA IMESS 1) ;
  697. ************Message 'Erreur absolue' epzzu.igau
  698. ************        'Composante ' (EXTR LDEFO icomp)
  699. ************        ' Point Gauss' igau 'Etat' ietat ;
  700. ************Message ;
  701. ************FINSI ;
  702.             ierr = 20 ; errabs = epzzu.igau ; QUITTER BECART ;
  703.          FINSI ;
  704.       SINON ;
  705.          depzz = ABS (epzzc.igau - epzzu.igau) ;
  706.          errrel = ABS (depzz / epzzc.igau) ;
  707.          SI (errrel > critere) ;
  708. ************SI (EGA IMESS 1) ;
  709. ************Message 'Erreur relative' errrel
  710. ************        'Composante ' (EXTR LDEFO icomp)
  711. ************        ' Point Gauss' igau 'Etat' ietat ;
  712. ************Message 'epzzc.igau =' epzzc.igau 'epzzu.igau =' epzzu.igau;
  713. ************Message ;
  714. ************FINSI ;
  715.             ierr = 21 ; QUITTER BECART ;
  716.          FINSI ;
  717.       FINSI ;
  718.    FIN BEPZZ ;
  719.  
  720.    icomp = 4 ;
  721.    igau = 0 ;
  722.    REPETER BGAXY nbgau ;
  723.       igau = igau + 1 ;
  724.       SI (EGA gaxyc.igau 0.0) ;
  725.          SI (NON (EGA gaxyu.igau 0.0)) ;
  726.             SI (EGA IMESS 1) ;
  727.             Message 'Erreur absolue' gaxyu.igau
  728.                     'Composante ' (EXTR LDEFO icomp)
  729.                     ' Point Gauss' igau 'Etat' ietat ;
  730.             Message ;
  731.             FINSI ;
  732. ************ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  733.          FINSI ;
  734.       SINON ;
  735.          dgaxy = ABS (gaxyc.igau - gaxyu.igau) ;
  736.          errrel = ABS (dgaxy / gaxyc.igau) ;
  737.          SI (errrel > critere) ;
  738.             SI (EGA IMESS 1) ;
  739.             Message 'Erreur relative' errrel
  740.                     'Composante ' (EXTR LDEFO icomp)
  741.                     ' Point Gauss' igau 'Etat' ietat ;
  742.             Message 'gaxyc.igau =' gaxyc.igau 'gaxyu.igau =' gaxyu.igau;
  743.             Message ;
  744.             FINSI ;
  745. ************ierr = 21 ; QUITTER BECART ;
  746.          FINSI ;
  747.       FINSI ;
  748.    FIN BGAXY ;
  749.  
  750.    icomp = 5 ;
  751.    igau = 0 ;
  752.    REPETER BGAXZ nbgau ;
  753.       igau = igau + 1 ;
  754.       SI (EGA gaxzc.igau 0.0) ;
  755.          SI (NON (EGA gaxzu.igau 0.0)) ;
  756.             SI (EGA IMESS 1) ;
  757.             Message 'Erreur absolue' gaxzu.igau
  758.                     'Composante ' (EXTR LDEFO icomp)
  759.                     ' Point Gauss' igau 'Etat' ietat ;
  760.             Message ;
  761.             FINSI ;
  762. ************ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  763.          FINSI ;
  764.       SINON ;
  765.          dgaxz = ABS (gaxzc.igau - gaxzu.igau) ;
  766.          errrel = ABS (dgaxz / gaxzc.igau) ;
  767.          SI (errrel > critere) ;
  768.             SI (EGA IMESS 1) ;
  769.             Message 'Erreur relative' errrel
  770.                     'Composante ' (EXTR LDEFO icomp)
  771.                     ' Point Gauss' igau 'Etat' ietat ;
  772.             Message 'gaxzc.igau =' gaxzc.igau 'gaxzu.igau =' gaxzu.igau;
  773.             Message ;
  774.             FINSI ;
  775. ************ierr = 21 ; QUITTER BECART ;
  776.          FINSI ;
  777.       FINSI ;
  778.    FIN BGAXZ ;
  779.  
  780.    icomp = 6 ;
  781.    igau = 0 ;
  782.    REPETER BGAYZ nbgau ;
  783.       igau = igau + 1 ;
  784.       SI (EGA gayzc.igau 0.0) ;
  785.          SI (NON (EGA gayzu.igau 0.0)) ;
  786.             SI (EGA IMESS 1) ;
  787.             Message 'Erreur absolue' gayzu.igau
  788.                     'Composante ' (EXTR LDEFO icomp)
  789.                     ' Point Gauss' igau 'Etat' ietat ;
  790.             Message ;
  791.             FINSI ;
  792. ************ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  793.          FINSI ;
  794.       SINON ;
  795.          dgayz = ABS (gayzc.igau - gayzu.igau) ;
  796.          errrel = ABS (dgayz / gayzc.igau) ;
  797.          SI (errrel > critere) ;
  798.             SI (EGA IMESS 1) ;
  799.             Message 'Erreur relative' errrel
  800.                     'Composante ' (EXTR LDEFO icomp)
  801.                     ' Point Gauss' igau 'Etat' ietat ;
  802.             Message 'gayzc.igau =' gayzc.igau 'gayzu.igau =' gayzu.igau;
  803.             Message ;
  804.             FINSI ;
  805. ************ierr = 21 ; QUITTER BECART ;
  806.          FINSI ;
  807.       FINSI ;
  808.    FIN BGAYZ ;
  809.  
  810. *  Comparaison des deformations inelastiques----------------------------
  811.  
  812.    icomp = 1 ;
  813.    igau = 0 ;
  814.    REPETER BEIXX nbgau ;
  815.       igau = igau + 1 ;
  816.       SI (EGA eixxc.igau 0.0) ;
  817.          SI (NON (EGA eixxu.igau 0.0)) ;
  818.             SI (EGA IMESS 1) ;
  819.             Message 'Erreur absolue' eixxu.igau
  820.                     'Composante ' (EXTR LDEIN icomp)
  821.                     ' Point Gauss' igau 'Etat' ietat ;
  822.             Message ;
  823.             FINSI ;
  824. ************ierr = 30 ; errabs = eixxu.igau ; QUITTER BECART ;
  825.          FINSI ;
  826.       SINON ;
  827.          deixx = ABS (eixxc.igau - eixxu.igau) ;
  828.          errrel = ABS (deixx / eixxc.igau) ;
  829.          SI (errrel > critere) ;
  830.             SI (EGA IMESS 1) ;
  831.             Message 'Erreur relative' errrel
  832.                     'Composante ' (EXTR LDEIN icomp)
  833.                     ' Point Gauss' igau 'Etat' ietat ;
  834.             Message 'eixxc.igau =' eixxc.igau 'eixxu.igau =' eixxu.igau;
  835.             Message ;
  836.             FINSI ;
  837. ************ierr = 31 ; QUITTER BECART ;
  838.          FINSI ;
  839.       FINSI ;
  840.    FIN BEIXX ;
  841.  
  842.    icomp = 2 ;
  843.    igau = 0 ;
  844.    REPETER BEIYY nbgau ;
  845.       igau = igau + 1 ;
  846.       SI (EGA eiyyc.igau 0.0) ;
  847.          SI (NON (EGA eiyyu.igau 0.0)) ;
  848.             SI (EGA IMESS 1) ;
  849.             Message 'Erreur absolue' eiyyu.igau
  850.                     'Composante ' (EXTR LDEIN icomp)
  851.                     ' Point Gauss' igau 'Etat' ietat ;
  852.             Message ;
  853.             FINSI ;
  854. ************ierr = 30 ; errabs = eiyyu.igau ; QUITTER BECART ;
  855.          FINSI ;
  856.       SINON ;
  857.          deiyy = ABS (eiyyc.igau - eiyyu.igau) ;
  858.          errrel = ABS (deiyy / eiyyc.igau) ;
  859.          SI (errrel > critere) ;
  860.             SI (EGA IMESS 1) ;
  861.             Message 'Erreur relative' errrel
  862.                     'Composante ' (EXTR LDEIN icomp)
  863.                     ' Point Gauss' igau 'Etat' ietat ;
  864.             Message 'eiyyc.igau =' eiyyc.igau 'eiyyu.igau =' eiyyu.igau;
  865.             Message ;
  866.             FINSI ;
  867. ************ierr = 31 ; QUITTER BECART ;
  868.          FINSI ;
  869.       FINSI ;
  870.    FIN BEIYY ;
  871.  
  872.    icomp = 3 ;
  873.    igau = 0 ;
  874.    REPETER BEIZZ nbgau ;
  875.       igau = igau + 1 ;
  876.       SI (EGA eizzc.igau 0.0) ;
  877.          SI (NON (EGA eizzu.igau 0.0)) ;
  878. ************SI (EGA IMESS 1) ;
  879. ************Message 'Erreur absolue' eizzu.igau
  880. ************        'Composante ' (EXTR LDEIN icomp)
  881. ************        ' Point Gauss' igau 'Etat' ietat ;
  882. ************Message ;
  883. ************FINSI ;
  884.             ierr = 30 ; errabs = eizzu.igau ; QUITTER BECART ;
  885.          FINSI ;
  886.       SINON ;
  887.          deizz = ABS (eizzc.igau - eizzu.igau) ;
  888.          errrel = ABS (deizz / eizzc.igau) ;
  889.          SI (errrel > critere) ;
  890. ************SI (EGA IMESS 1) ;
  891. ************Message 'Erreur relative' errrel
  892. ************        'Composante ' (EXTR LDEIN icomp)
  893. ************        ' Point Gauss' igau 'Etat' ietat ;
  894. ************Message 'eizzc.igau =' eizzc.igau 'eizzu.igau =' eizzu.igau;
  895. ************Message ;
  896. ************FINSI ;
  897.             ierr = 31 ; QUITTER BECART ;
  898.          FINSI ;
  899.       FINSI ;
  900.    FIN BEIZZ ;
  901.  
  902.    icomp = 4 ;
  903.    igau = 0 ;
  904.    REPETER BGIXY nbgau ;
  905.       igau = igau + 1 ;
  906.       SI (EGA gixyc.igau 0.0) ;
  907.          SI (NON (EGA gixyu.igau 0.0)) ;
  908.             SI (EGA IMESS 1) ;
  909.             Message 'Erreur absolue' gixyu.igau
  910.                     'Composante ' (EXTR LDEIN icomp)
  911.                     ' Point Gauss' igau 'Etat' ietat ;
  912.             Message ;
  913.             FINSI ;
  914. ************ierr = 30 ; errabs = gixyu.igau ; QUITTER BECART ;
  915.          FINSI ;
  916.       SINON ;
  917.          dgixy = ABS (gixyc.igau - gixyu.igau) ;
  918.          errrel = ABS (dgixy / gixyc.igau) ;
  919.          SI (errrel > critere) ;
  920.             SI (EGA IMESS 1) ;
  921.             Message 'Erreur relative' errrel
  922.                     'Composante ' (EXTR LDEIN icomp)
  923.                     ' Point Gauss' igau 'Etat' ietat ;
  924.             Message 'gixyc.igau =' gixyc.igau 'gixyu.igau =' gixyu.igau;
  925.             Message ;
  926.             FINSI ;
  927. ************ierr = 31 ; QUITTER BECART ;
  928.          FINSI ;
  929.       FINSI ;
  930.    FIN BGIXY ;
  931.  
  932.    icomp = 5 ;
  933.    igau = 0 ;
  934.    REPETER BGIXZ nbgau ;
  935.       igau = igau + 1 ;
  936.       SI (EGA gixzc.igau 0.0) ;
  937.          SI (NON (EGA gixzu.igau 0.0)) ;
  938.             SI (EGA IMESS 1) ;
  939.             Message 'Erreur absolue' gixzu.igau
  940.                     'Composante ' (EXTR LDEIN icomp)
  941.                     ' Point Gauss' igau 'Etat' ietat ;
  942.             Message ;
  943.             FINSI ;
  944. ************ierr = 30 ; errabs = gixzu.igau ; QUITTER BECART ;
  945.          FINSI ;
  946.       SINON ;
  947.          dgixz = ABS (gixzc.igau - gixzu.igau) ;
  948.          errrel = ABS (dgixz / gixzc.igau) ;
  949.          SI (errrel > critere) ;
  950.             SI (EGA IMESS 1) ;
  951.             Message 'Erreur relative' errrel
  952.                     'Composante ' (EXTR LDEIN icomp)
  953.                     ' Point Gauss' igau 'Etat' ietat ;
  954.             Message 'gixzc.igau =' gixzc.igau 'gixzu.igau =' gixzu.igau;
  955.             Message ;
  956.             FINSI ;
  957. ************ierr = 31 ; QUITTER BECART ;
  958.          FINSI ;
  959.       FINSI ;
  960.    FIN BGIXZ ;
  961.  
  962.    icomp = 6 ;
  963.    igau = 0 ;
  964.    REPETER BGIYZ nbgau ;
  965.       igau = igau + 1 ;
  966.       SI (EGA giyzc.igau 0.0) ;
  967.          SI (NON (EGA giyzu.igau 0.0)) ;
  968.             SI (EGA IMESS 1) ;
  969.             Message 'Erreur absolue' giyzu.igau
  970.                     'Composante ' (EXTR LDEIN icomp)
  971.                     ' Point Gauss' igau 'Etat' ietat ;
  972.             Message ;
  973.             FINSI ;
  974. ************ierr = 30 ; errabs = giyzu.igau ; QUITTER BECART ;
  975.          FINSI ;
  976.       SINON ;
  977.          dgiyz = ABS (giyzc.igau - giyzu.igau) ;
  978.          errrel = ABS (dgiyz / giyzc.igau) ;
  979.          SI (errrel > critere) ;
  980.             SI (EGA IMESS 1) ;
  981.             Message 'Erreur relative' errrel
  982.                     'Composante ' (EXTR LDEIN icomp)
  983.                     ' Point Gauss' igau 'Etat' ietat ;
  984.             Message 'giyzc.igau =' giyzc.igau 'giyzu.igau =' giyzu.igau;
  985.             Message ;
  986.             FINSI ;
  987. ************ierr = 31 ; QUITTER BECART ;
  988.          FINSI ;
  989.       FINSI ;
  990.    FIN BGIYZ ;
  991.  
  992. *  Comparaison des contraintes------------------------------------------
  993.  
  994.    icomp = 1 ;
  995.    igau = 0 ;
  996.    REPETER BSMXX nbgau ;
  997.       igau = igau + 1 ;
  998.       SI (EGA smxxc.igau 0.0) ;
  999.          SI (NON (EGA smxxu.igau 0.0)) ;
  1000.             SI (EGA IMESS 1) ;
  1001.             Message 'Erreur absolue' smxxu.igau
  1002.                     'Composante ' (EXTR LCONT icomp)
  1003.                     ' Point Gauss' igau 'Etat' ietat ;
  1004.             Message ;
  1005.             FINSI ;
  1006. ************ierr = 40 ; errabs = smxxu.igau ; QUITTER BECART ;
  1007.          FINSI ;
  1008.       SINON ;
  1009.          dsmxx = ABS (smxxc.igau - smxxu.igau) ;
  1010.          errrel = ABS (dsmxx / smxxc.igau) ;
  1011.          SI (errrel > critere) ;
  1012.             SI (EGA IMESS 1) ;
  1013.             Message 'Erreur relative' errrel
  1014.                     'Composante ' (EXTR LCONT icomp)
  1015.                     ' Point Gauss' igau 'Etat' ietat ;
  1016.             Message 'smxxc.igau =' smxxc.igau 'smxxu.igau =' smxxu.igau;
  1017.             Message ;
  1018.             FINSI ;
  1019. ************ierr = 41 ; QUITTER BECART ;
  1020.          FINSI ;
  1021.       FINSI ;
  1022.    FIN BSMXX ;
  1023.  
  1024.    icomp = 2 ;
  1025.    igau = 0 ;
  1026.    REPETER BSMYY nbgau ;
  1027.       igau = igau + 1 ;
  1028.       SI (EGA smyyc.igau 0.0) ;
  1029.          SI (NON (EGA smyyu.igau 0.0)) ;
  1030.             SI (EGA IMESS 1) ;
  1031.             Message 'Erreur absolue' smyyu.igau
  1032.                     'Composante ' (EXTR LCONT icomp)
  1033.                     ' Point Gauss' igau 'Etat' ietat ;
  1034.             Message ;
  1035.             FINSI ;
  1036. ************ierr = 40 ; errabs = smyyu.igau ; QUITTER BECART ;
  1037.          FINSI ;
  1038.       SINON ;
  1039.          dsmyy = ABS (smyyc.igau - smyyu.igau) ;
  1040.          errrel = ABS (dsmyy / smyyc.igau) ;
  1041.          SI (errrel > critere) ;
  1042.             SI (EGA IMESS 1) ;
  1043.             Message 'Erreur relative' errrel
  1044.                     'Composante ' (EXTR LCONT icomp)
  1045.                     ' Point Gauss' igau 'Etat' ietat ;
  1046.             Message 'smyyc.igau =' smyyc.igau 'smyyu.igau =' smyyu.igau;
  1047.             Message ;
  1048.             FINSI ;
  1049. ************ierr = 41 ; QUITTER BECART ;
  1050.          FINSI ;
  1051.       FINSI ;
  1052.    FIN BSMYY ;
  1053.  
  1054.    icomp = 3 ;
  1055.    igau = 0 ;
  1056.    REPETER BSMZZ nbgau ;
  1057.       igau = igau + 1 ;
  1058.       SI (EGA smzzc.igau 0.0) ;
  1059.          SI (NON (EGA smzzu.igau 0.0)) ;
  1060. ************SI (EGA IMESS 1) ;
  1061. ************Message 'Erreur absolue' smzzu.igau
  1062. ************        'Composante ' (EXTR LCONT icomp)
  1063. ************        ' Point Gauss' igau 'Etat' ietat ;
  1064. ************Message ;
  1065. ************FINSI ;
  1066.             ierr = 40 ; errabs = smzzu.igau ; QUITTER BECART ;
  1067.          FINSI ;
  1068.       SINON ;
  1069.          dsmzz = ABS (smzzc.igau - smzzu.igau) ;
  1070.          errrel = ABS (dsmzz / smzzc.igau) ;
  1071.          SI (errrel > critere) ;
  1072. ************SI (EGA IMESS 1) ;
  1073. ************Message 'Erreur relative' errrel
  1074. ************        'Composante ' (EXTR LCONT icomp)
  1075. ************        ' Point Gauss' igau 'Etat' ietat ;
  1076. ************Message 'smzzc.igau =' smzzc.igau 'smzzu.igau =' smzzu.igau;
  1077. ************Message ;
  1078. ************FINSI ;
  1079.             ierr = 41 ; QUITTER BECART ;
  1080.          FINSI ;
  1081.       FINSI ;
  1082.    FIN BSMZZ ;
  1083.  
  1084.    icomp = 4 ;
  1085.    igau = 0 ;
  1086.    REPETER BSMXY nbgau ;
  1087.       igau = igau + 1 ;
  1088.       SI (EGA smxyc.igau 0.0) ;
  1089.          SI (NON (EGA smxyu.igau 0.0)) ;
  1090.             SI (EGA IMESS 1) ;
  1091.             Message 'Erreur absolue' smxyu.igau
  1092.                     'Composante ' (EXTR LCONT icomp)
  1093.                     ' Point Gauss' igau 'Etat' ietat ;
  1094.             Message ;
  1095.             FINSI ;
  1096. ************ierr = 40 ; errabs = smxyu.igau ; QUITTER BECART ;
  1097.          FINSI ;
  1098.       SINON ;
  1099.          dsmxy = ABS (smxyc.igau - smxyu.igau) ;
  1100.          errrel = ABS (dsmxy / smxyc.igau) ;
  1101.          SI (errrel > critere) ;
  1102.             SI (EGA IMESS 1) ;
  1103.             Message 'Erreur relative' errrel
  1104.                     'Composante ' (EXTR LCONT icomp)
  1105.                     ' Point Gauss' igau 'Etat' ietat ;
  1106.             Message 'smxyc.igau =' smxyc.igau 'smxyu.igau =' smxyu.igau;
  1107.             Message ;
  1108.             FINSI ;
  1109. ************ierr = 41 ; QUITTER BECART ;
  1110.          FINSI ;
  1111.       FINSI ;
  1112.    FIN BSMXY ;
  1113.  
  1114.    icomp = 5 ;
  1115.    igau = 0 ;
  1116.    REPETER BSMXZ nbgau ;
  1117.       igau = igau + 1 ;
  1118.       SI (EGA smxzc.igau 0.0) ;
  1119.          SI (NON (EGA smxzu.igau 0.0)) ;
  1120.             SI (EGA IMESS 1) ;
  1121.             Message 'Erreur absolue' smxzu.igau
  1122.                     'Composante ' (EXTR LCONT icomp)
  1123.                     ' Point Gauss' igau 'Etat' ietat ;
  1124.             Message ;
  1125.             FINSI ;
  1126. ************ierr = 40 ; errabs = smxzu.igau ; QUITTER BECART ;
  1127.          FINSI ;
  1128.       SINON ;
  1129.          dsmxz = ABS (smxzc.igau - smxzu.igau) ;
  1130.          errrel = ABS (dsmxz / smxzc.igau) ;
  1131.          SI (errrel > critere) ;
  1132.             SI (EGA IMESS 1) ;
  1133.             Message 'Erreur relative' errrel
  1134.                     'Composante ' (EXTR LCONT icomp)
  1135.                     ' Point Gauss' igau 'Etat' ietat ;
  1136.             Message 'smxzc.igau =' smxzc.igau 'smxzu.igau =' smxzu.igau;
  1137.             Message ;
  1138.             FINSI ;
  1139. ************ierr = 41 ; QUITTER BECART ;
  1140.          FINSI ;
  1141.       FINSI ;
  1142.    FIN BSMXZ ;
  1143.  
  1144.    icomp = 6 ;
  1145.    igau = 0 ;
  1146.    REPETER BSMYZ nbgau ;
  1147.       igau = igau + 1 ;
  1148.       SI (EGA smyzc.igau 0.0) ;
  1149.          SI (NON (EGA smyzu.igau 0.0)) ;
  1150.             SI (EGA IMESS 1) ;
  1151.             Message 'Erreur absolue' smyzu.igau
  1152.                     'Composante ' (EXTR LCONT icomp)
  1153.                     ' Point Gauss' igau 'Etat' ietat ;
  1154.             Message ;
  1155.             FINSI ;
  1156. ************ierr = 40 ; errabs = smyzu.igau ; QUITTER BECART ;
  1157.          FINSI ;
  1158.       SINON ;
  1159.          dsmyz = ABS (smyzc.igau - smyzu.igau) ;
  1160.          errrel = ABS (dsmyz / smyzc.igau) ;
  1161.          SI (errrel > critere) ;
  1162.             SI (EGA IMESS 1) ;
  1163.             Message 'Erreur relative' errrel
  1164.                     'Composante ' (EXTR LCONT icomp)
  1165.                     ' Point Gauss' igau 'Etat' ietat ;
  1166.             Message 'smyzc.igau =' smyzc.igau 'smyzu.igau =' smyzu.igau;
  1167.             Message ;
  1168.             FINSI ;
  1169. ************ierr = 41 ; QUITTER BECART ;
  1170.          FINSI ;
  1171.       FINSI ;
  1172.    FIN BSMYZ ;
  1173.  
  1174. FIN BECART ;
  1175.  
  1176. * FIN   BOUCLE COMPARAISON DES RESULTATS________________________________
  1177.  
  1178. * MESSAGES D'ERREUR_____________________________________________________
  1179.  
  1180. Message ;
  1181. Message 'CPU avec modele interne = ' CPUint 'centiemes de seconde' ;
  1182. Message 'CPU avec modele externe = ' CPUext 'centiemes de seconde' ;
  1183. Message ;
  1184.  
  1185. ifin = netat-1 ;
  1186.  
  1187. epzzcmax = 0.0 ;
  1188. igau = 0 ;
  1189. REPETER BEPZZFC nbgau ;
  1190.    igau = igau + 1 ;
  1191.    epzz = EXTR (tabflup.'DEFORMATIONS'.ifin)
  1192.                     'EPZZ' 1 1 igau ;
  1193.    SI (epzz > epzzcmax) ; epzzcmax = epzz ; FINSI ;
  1194. FIN BEPZZFC ;
  1195. epzzumax = 0.0 ;
  1196. igau = 0 ;
  1197. REPETER BEPZZFU nbgau ;
  1198.    igau = igau + 1 ;
  1199.    epzz = EXTR (tabutil.'DEFORMATIONS'.ifin)
  1200.                     'EPZZ' 1 1 igau ;
  1201.    SI (epzz > epzzumax) ; epzzumax = epzz ; FINSI ;
  1202. FIN BEPZZFU ;
  1203. Message 'Deformation totale finale (maxi sur points de Gauss)' ;
  1204. Message 'EPZZ max avec modele interne = ' epzzcmax ;
  1205. Message 'EPZZ max avec modele externe = ' epzzumax ;
  1206. Message ;
  1207.  
  1208. eizzcmax = 0.0 ;
  1209. igau = 0 ;
  1210. REPETER BEIZZFC nbgau ;
  1211.    igau = igau + 1 ;
  1212.    eizz = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ifin)
  1213.                     'EIZZ' 1 1 igau ;
  1214.    SI (eizz > eizzcmax) ; eizzcmax = eizz ; FINSI ;
  1215. FIN BEIZZFC ;
  1216. eizzumax = 0.0 ;
  1217. igau = 0 ;
  1218. REPETER BEIZZFU nbgau ;
  1219.    igau = igau + 1 ;
  1220.    eizz = EXTR (tabutil.'DEFORMATIONS_INELASTIQUES'.ifin)
  1221.                     'EIZZ' 1 1 igau ;
  1222.    SI (eizz > eizzumax) ; eizzumax = eizz ; FINSI ;
  1223. FIN BEIZZFU ;
  1224. Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  1225. Message 'EIZZ max avec modele interne = ' eizzcmax ;
  1226. Message 'EIZZ max avec modele externe = ' eizzumax ;
  1227. Message ;
  1228.  
  1229. smzzcmax = 0.0 ;
  1230. igau = 0 ;
  1231. REPETER BSMZZFC nbgau ;
  1232.    igau = igau + 1 ;
  1233.    smzz = EXTR (tabflup.'CONTRAINTES'.ifin)
  1234.                     'SMZZ' 1 1 igau ;
  1235.    SI (smzz > smzzcmax) ; smzzcmax = smzz ; FINSI ;
  1236. FIN BSMZZFC ;
  1237. smzzumax = 0.0 ;
  1238. igau = 0 ;
  1239. REPETER BSMZZFU nbgau ;
  1240.    igau = igau + 1 ;
  1241.    smzz = EXTR (tabutil.'CONTRAINTES'.ifin)
  1242.                     'SMZZ' 1 1 igau ;
  1243.    SI (smzz > smzzumax) ; smzzumax = smzz ; FINSI ;
  1244. FIN BSMZZFU ;
  1245. Message 'Contrainte finale (maxi sur points de Gauss)' ;
  1246. Message 'SMZZ max avec modele interne = ' smzzcmax ;
  1247. Message 'SMZZ max avec modele externe = ' smzzumax ;
  1248. Message ;
  1249.  
  1250. SI (ierr > 0) ;
  1251.  
  1252.    SI (EGA ierr 1) ;
  1253.       Message 'Instants de calcul non coincidants' ;
  1254.       ERRE 5 ;
  1255.    FINSI ;
  1256.  
  1257.    SI (EGA ierr 10) ;
  1258.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  1259.               ' Noeud' ino 'Etat' ietat ;
  1260.       ERRE 5 ;
  1261.    FINSI ;
  1262.  
  1263.    SI (EGA ierr 11) ;
  1264.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  1265.               ' Noeud' ino 'Etat' ietat ;
  1266.       ERRE 5 ;
  1267.    FINSI ;
  1268.  
  1269.    SI (EGA ierr 20) ;
  1270.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  1271.               ' Point Gauss' igau 'Etat' ietat ;
  1272.       ERRE 5 ;
  1273.    FINSI ;
  1274.  
  1275.    SI (EGA ierr 21) ;
  1276.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  1277.               ' Point Gauss' igau 'Etat' ietat ;
  1278.       ERRE 5 ;
  1279.    FINSI ;
  1280.  
  1281.    SI (EGA ierr 30) ;
  1282.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  1283.               ' Point Gauss' igau 'Etat' ietat ;
  1284.       ERRE 5 ;
  1285.    FINSI ;
  1286.  
  1287.    SI (EGA ierr 31) ;
  1288.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  1289.               ' Point Gauss' igau 'Etat' ietat ;
  1290.       ERRE 5 ;
  1291.    FINSI ;
  1292.  
  1293.    SI (EGA ierr 40) ;
  1294.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  1295.               ' Point Gauss' igau 'Etat' ietat ;
  1296.       ERRE 5 ;
  1297.    FINSI ;
  1298.  
  1299.    SI (EGA ierr 41) ;
  1300.       Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  1301.               ' Point Gauss' igau 'Etat' ietat ;
  1302.       ERRE 5 ;
  1303.    FINSI ;
  1304.  
  1305.    SI (EGA ierr 50) ;
  1306.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  1307.               ' Point Gauss' igau 'Etat' ietat ;
  1308.       ERRE 5 ;
  1309.    FINSI ;
  1310.  
  1311.    SI (EGA ierr 51) ;
  1312.       Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  1313.               ' Point Gauss' igau 'Etat' ietat ;
  1314.       ERRE 5 ;
  1315.    FINSI ;
  1316.  
  1317. SINON ;
  1318.  
  1319.    Message 'Resultats concordants entre les deux modeles' ;
  1320.    Message 'Aucune erreur relative superieure a ' critere
  1321.            ' pour les composantes' ;
  1322.    Message (EXTR LDEFO idefo) '  ' (EXTR LDEIN idein)
  1323.                               '  ' (EXTR LCONT icont) ;
  1324.    ERRE 0 ;
  1325.  
  1326. FINSI ;
  1327.  
  1328. FIN ;
  1329.  
  1330.  
  1331.  
  1332.  
  1333.  
  1334.  
  1335.  
  1336.  
  1337.  
  1338.  

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