Télécharger la plaquette Cast3M


fr en

Télécharger creep01_cisXY.dgibi

Retour à la liste

Numérotation des lignes : 

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

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