Télécharger la plaquette Cast3M


fr en

Télécharger creep01_traXX.dgibi

Retour à la liste

Numérotation des lignes : 

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

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