Télécharger la plaquette Cast3M


fr en

Télécharger creep02_cisXY.dgibi

Retour à la liste

Numérotation des lignes : 

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

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