Télécharger la plaquette Cast3M


fr en

Télécharger creep01_cisXZ.dgibi

Retour à la liste

Numérotation des lignes : 

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

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