Télécharger la plaquette Cast3M


fr en

Télécharger creep01_cisYZ.dgibi

Retour à la liste

Numérotation des lignes : 

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

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