Télécharger la plaquette Cast3M


fr en

Télécharger creep01_traYY.dgibi

Retour à la liste

Numérotation des lignes : 

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

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