Télécharger umat03_traXX.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : umat03_traXX.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *=======================================================================
  5. *
  6. * CAS TEST DE VALIDATION DE LA PRISE EN COMPTE D'UNE LOI DE
  7. *
  8. * COMPORTEMENT MECANIQUE NON LINEAIRE EXTERNE
  9. *
  10. *=======================================================================
  11. *
  12. * Modele 'NON_LINEAIRE' 'UTILISATEUR', integrateur specifique UMAT
  13. *
  14. * Test No3 : Modele de FLUAGE POLYNOMIAL
  15. * TRACTION XX
  16. *
  17. *=======================================================================
  18.  
  19.  
  20. opti echo 0 ;
  21. * ajout de option epsilon lineaire pour la precision des test!
  22. OPTION epsilon lineaire;
  23.  
  24. *=======================================================================
  25. *
  26. * PROCEDURES DEPENDANTES
  27. *
  28. *=======================================================================
  29.  
  30.  
  31.  
  32.  
  33. *=======================================================================
  34. *
  35. * JEU DE DONNEES PRINCIPAL
  36. *
  37. *=======================================================================
  38.  
  39. *-----------------------------------------------------------------------
  40. * 0 - Maillage
  41. * PLAQUE CARREE COTE : 1 m EPAISSEUR : 1 cm
  42.  
  43. opti dime 3 elem cub8 mode tridim ;
  44.  
  45. l1 = 1.0 ;
  46. l2 = 1.0 ;
  47. ep0 = 0.01 ;
  48.  
  49. p01 = 0.0 0.0 0.0 ;
  50. p02 = ep0 0.0 0.0 ;
  51. p03 = ep0 l2 0.0 ;
  52. p04 = 0.0 l2 0.0 ;
  53.  
  54. p10 = 0.0 0.0 l1 ;
  55. p20 = ep0 0.0 l1 ;
  56. p30 = ep0 l2 l1 ;
  57. p40 = 0.0 l2 l1 ;
  58.  
  59. li01 = 'DROIT' 1 p01 p02 ;
  60. li02 = 'DROIT' 1 p02 p03 ;
  61. li03 = 'DROIT' 1 p03 p04 ;
  62. li04 = 'DROIT' 1 p04 p01 ;
  63.  
  64. li10 = 'DROIT' 1 p10 p20 ;
  65. li20 = 'DROIT' 1 p20 p30 ;
  66. li30 = 'DROIT' 1 p30 p40 ;
  67. li40 = 'DROIT' 1 p40 p10 ;
  68.  
  69. sur01 = 'DALL' li01 li02 li03 li04 ;
  70. sur10 = 'DALL' li10 li20 li30 li40 ;
  71. mail1 = sur01 'VOLU' 1 sur10 ;
  72.  
  73. li2 = 'DROIT' 1 p01 p10 ;
  74. li3 = 'DROIT' 1 p10 p40 ;
  75. li4 = 'DROIT' 1 p40 p04 ;
  76. sbloq = 'DALL' li04 li2 li3 li4 ;
  77.  
  78. li5 = 'DROIT' 1 p03 p02 ;
  79. li6 = 'DROIT' 1 p02 p20 ;
  80. li7 = 'DROIT' 1 p20 p30 ;
  81. li8 = 'DROIT' 1 p30 p03 ;
  82. schar = 'DALL' li5 li6 li7 li8 ;
  83.  
  84. ELIM (mail1 ET sbloq ET schar) 1.e-5 ;
  85.  
  86. *-----------------------------------------------------------------------
  87. * Comportement 'FLUAGE' 'POLYNOMIAL'
  88. * Chargement en forces imposees, CISAILLEMENT PUR
  89. *-----------------------------------------------------------------------
  90. * 1 - Calcul avec le modele 'FLUAGE' 'POLYNOMIAL' de CAST3M
  91.  
  92. * 1.1 Definition du modele
  93.  
  94. mo_flupo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  95. 'FLUAGE' 'POLYNOMIAL' ;
  96.  
  97. * 1.2 Affectation des constantes de materiau
  98.  
  99. densiu = 0.95 ;
  100. dgrain = 8.0E-6 ;
  101. tuo2 = 1500.0 ;
  102. tauxfisk = 1.5E19 ;
  103.  
  104. xyoun = (2.2693E11 - (1.5399E7*tuo2)) - (9.35973E3 * (tuo2**2.0)) ;
  105. xnu = 0.34 ;
  106. xrho = 10950.0 ;
  107. xalph = ( (9.828E-6 - (6.390E-10*tuo2)) + (1.330E-12 * (tuo2**2.0)) )
  108. - (1.757E-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.0E5 ;
  136. depmax = 3.0E10 ;
  137. dep0 = FORCE schar 'FX' depmax 'FY' 0.0 'FZ' 0.0 ;
  138. ev0 = EVOL MANU temps (PROG 0.0 1.E-2 1000.0)
  139. y (PROG 0.0 1.0 1.0 ) ;
  140. cha0 = CHAR MECA dep0 ev0 ;
  141.  
  142. * 1.4 Definition des instants de calcul
  143.  
  144. pro0 = PROG 0.0 PAS 2.E-4 1.E-2 ;
  145. pro1 = PROG 0.0 PAS 2.E-4 1.E-2 ;
  146.  
  147. * 1.5 Resolution par PASAPAS
  148.  
  149. tabflup = TABLE ;
  150. tabflup.'MODELE' = mo_flupo ;
  151. tabflup.'CARACTERISTIQUES' = ma_flupo ;
  152. tabflup.'CHARGEMENT' = cha0 ;
  153. tabflup.'BLOCAGES_MECANIQUES' = cl0 ;
  154. tabflup.'TEMPS_CALCULES' = pro0 ;
  155. tabflup.'TEMPS_SAUVES' = pro1 ;
  156. TMASAU=table;
  157. tabflup . 'MES_SAUVEGARDES'=TMASAU;
  158. tabflup . 'PRECISION'=1E-8;
  159. TMASAU .'DEFTO'=VRAI;
  160. TMASAU .'DEFIN'=VRAI;
  161. TEMPS 'ZERO' ;
  162. PASAPAS tabflup ;
  163. TABTPS = TEMP 'NOEC';
  164. CPUint = TABTPS.'TEMPS_CPU'.'INITIAL';
  165.  
  166. *-----------------------------------------------------------------------
  167. * 2 - Calcul avec modele utilisateur
  168.  
  169. * 2.1 Definition du modele utilisateur
  170.  
  171. LCMAT21 = MOTS 'YOUN' 'NU ' 'RHO ' 'ALPH'
  172. 'SMAX' 'AF0 ' 'AF1 ' 'AF2 ' 'AF3 ' 'AF4 ' 'AF5 ' 'AF6 ' ;
  173.  
  174. LCVAR21 = MOTS 'PXX ' 'PYY ' 'PZZ ' 'PXY ' 'PYZ ' 'PZX '
  175. 'VNXX' 'VNYY' 'VNZZ' 'VNXY' 'VNYZ' 'VNZX'
  176. 'EPCE' 'EPSH' 'LL ' 'EPSE'
  177. 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  178.  
  179. mo_util = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  180. 'NON_LINEAIRE' 'UTILISATEUR'
  181. 'NUME_LOI' 21
  182. 'C_MATERIAU' LCMAT21 'C_VARINTER' LCVAR21 ;
  183.  
  184. * 2.2 Affectation des proprietes de materiau
  185.  
  186. ma_util = MATE mo_util 'YOUN' xyoun 'NU ' xnu 'RHO ' xrho
  187. 'ALPH' xalph
  188. 'SMAX' xsmax 'AF0 ' xaf0 'AF1 ' xaf1
  189. 'AF2 ' xaf2 'AF3 ' xaf3 'AF4 ' xaf4
  190. 'AF5 ' xaf5 'AF6 ' xaf6 ;
  191.  
  192. * 2.3 Definition des conditions aux limites et chargements : IDEM 1.3
  193.  
  194. * 2.4 Definition des instants de calcul : IDEM 1.4
  195.  
  196. * 2.5 Initialisation des variables internes
  197.  
  198. chpxx0 = MANU CHML mo_util 'PXX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  199. chpyy0 = MANU CHML mo_util 'PYY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  200. chpzz0 = MANU CHML mo_util 'PZZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  201. chpxy0 = MANU CHML mo_util 'PXY ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  202. chpyz0 = MANU CHML mo_util 'PYZ ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  203. chpzx0 = MANU CHML mo_util 'PZX ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  204.  
  205. chvnxx0 = MANU CHML mo_util 'VNXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  206. chvnyy0 = MANU CHML mo_util 'VNYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  207. chvnzz0 = MANU CHML mo_util 'VNZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  208. chvnxy0 = MANU CHML mo_util 'VNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  209. chvnyz0 = MANU CHML mo_util 'VNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  210. chvnzx0 = MANU CHML mo_util 'VNZX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  211.  
  • chepce0 = MANU CHML mo_util 'EPCE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chepsh0 = MANU CHML mo_util 'EPSH' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chll0 = MANU CHML mo_util 'LL ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chepse0 = MANU CHML mo_util 'EPSE' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  •  
  • chenxx0 = MANU CHML mo_util 'ENXX' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chenyy0 = MANU CHML mo_util 'ENYY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chenzz0 = MANU CHML mo_util 'ENZZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chgnxy0 = MANU CHML mo_util 'GNXY' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chgnxz0 = MANU CHML mo_util 'GNXZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  • chgnyz0 = MANU CHML mo_util 'GNYZ' 0.0 'STRESSES' 'TYPE' 'VARINTER' ;
  •  
  • var_ut0 = chpxx0 ET chpyy0 ET chpzz0 ET chpxy0 ET chpyz0 ET chpzx0
  • ET chvnxx0 ET chvnyy0 ET chvnzz0 ET chvnxy0 ET chvnyz0 ET chvnzx0
  • ET chepce0 ET chepsh0 ET chll0 ET chepse0
  • ET chenxx0 ET chenyy0 ET chenzz0 ET chgnxy0 ET chgnxz0 ET chgnyz0 ;
  •  
  • * 2.6 Resolution par PASAPAS
  •  
  • tabutil = TABLE ;
  • tabutil.'MODELE' = mo_util ;
  • tabutil.'CARACTERISTIQUES' = ma_util ;
  • tabutil.'CHARGEMENT' = cha0 ;
  • tabutil.'BLOCAGES_MECANIQUES' = cl0 ;
  • tabutil.'VARIABLES_INTERNES' = TABLE ;
  • tabutil.'VARIABLES_INTERNES' . 0 = var_ut0 ;
  • tabutil.'TEMPS_CALCULES' = pro0 ;
  • tabutil.'TEMPS_SAUVES' = pro1 ;
  • TMASAU=table;
  • tabutil . 'MES_SAUVEGARDES'=TMASAU;
  • tabutil . 'PRECISION'=1E-8;
  • TMASAU .'DEFTO'=VRAI;
  • TMASAU .'DEFIN'=VRAI;
  • TEMPS 'ZERO' ;
  • PASAPAS tabutil ;
  • TABTPS = TEMP 'NOEC';
  • CPUext = TABTPS.'TEMPS_CPU'.'INITIAL';
  •  
  • *-----------------------------------------------------------------------
  • * 3 - Post-traitement : comparaisons entre les deux modeles
  • *
  • * N.B. Comme on simule un essai de traction XX
  • * - Critere strict (errrel < 1.d-10) sur EPXX, EIXX et SMXX
  • * - Pour les autres composantes on note les depassements du
  • * critere si impression demandee
  •  
  • uxc = TABLE ; uyc = TABLE ; uzc = TABLE ;
  • epxxc = TABLE ; epyyc = TABLE ; epzzc = TABLE ;
  • gaxyc = TABLE ; gaxzc = TABLE ; gayzc = TABLE ;
  • eixxc = TABLE ; eiyyc = TABLE ; eizzc = TABLE ;
  • gixyc = TABLE ; gixzc = TABLE ; giyzc = TABLE ;
  • smxxc = TABLE ; smyyc = TABLE ; smzzc = TABLE ;
  • smxyc = TABLE ; smxzc = TABLE ; smyzc = TABLE ;
  • pxxc = TABLE ; pyyc = TABLE ; pzzc = TABLE ;
  • pxyc = TABLE ; pyzc = TABLE ; pzxc = TABLE ;
  • vnxxc = TABLE ; vnyyc = TABLE ; vnzzc = TABLE ;
  • vnxyc = TABLE ; vnyzc = TABLE ; vnzxc = TABLE ;
  • epcec = TABLE ; epshc = TABLE ; llc = TABLE ; epsec = TABLE ;
  •  
  • uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  • epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  • gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  • enxxu = TABLE ; enyyu = TABLE ; enzzu = TABLE ;
  • gnxyu = TABLE ; gnxzu = TABLE ; gnyzu = TABLE ;
  • smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  • smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  • pxxu = TABLE ; pyyu = TABLE ; pzzu = TABLE ;
  • pxyu = TABLE ; pyzu = TABLE ; pzxu = TABLE ;
  • vnxxu = TABLE ; vnyyu = TABLE ; vnzzu = TABLE ;
  • vnxyu = TABLE ; vnyzu = TABLE ; vnzxu = TABLE ;
  • epceu = TABLE ; epshu = TABLE ; llu = TABLE ; epseu = TABLE ;
  •  
  • LDEPL = MOTS 'UX ' 'UY ' 'UZ ' ;
  • LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  • LDEIN = MOTS 'EIXX' 'EIYY' 'EIZZ' 'GIXY' 'GIXZ' 'GIYZ' ;
  • LDENN = MOTS 'ENXX' 'ENYY' 'ENZZ' 'GNXY' 'GNXZ' 'GNYZ' ;
  • LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  • LVARI = MOTS 'PXX ' 'PYY ' 'PZZ ' 'PXY ' 'PYZ ' 'PZX '
  • 'VNXX' 'VNYY' 'VNZZ' 'VNXY' 'VNYZ' 'VNZX'
  • 'EPCE' 'EPSH' 'LL ' 'EPSE' ;
  •  
  • idefo = 1 ;
  • idein = 1 ;
  • icont = 1 ;
  •  
  • ierr = 0 ;
  • critere = 1.d-8 ;
  • IMESS = 0 ;
  • ***IMESS = 1 ;
  •  
  • nbno = 8 ;
  • nbgau = 8 ;
  •  
  • netat = DIME (tabflup.'TEMPS') ;
  • ietat = -1 ;
  •  
  •  
  • * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  •  
  • REPETER BECART netat ;
  •  
  • ietat = ietat + 1 ;
  •  
  • temps1 = tabflup.'TEMPS'.ietat ;
  • temps2 = tabutil.'TEMPS'.ietat ;
  • SI ((ABS (temps1-temps2)) > 1.0E-9) ;
  • ierr = 1 ;
  • QUITTER BECART ;
  • FINSI ;
  •  
  • * Extraction des resultats de calcul : Modele interne CAST3M-----------
  •  
  • * Deplacements
  • uxc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  • uxc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  • uxc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  • uxc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  • uxc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  • uxc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  • uxc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  • uxc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  •  
  • uyc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  • uyc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  • uyc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  • uyc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  • uyc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  • uyc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  • uyc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  • uyc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  •  
  • uzc.1 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  • uzc.2 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  • uzc.3 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  • uzc.4 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  • uzc.5 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  • uzc.6 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  • uzc.7 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  • uzc.8 = EXTR (tabflup.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  •  
  • * Deformations totales
  • igau = 0 ;
  • REPETER BEPXXC nbgau ;
  • igau = igau + 1 ;
  • epxxc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  • FIN BEPXXC ;
  •  
  • igau = 0 ;
  • REPETER BEPYYC nbgau ;
  • igau = igau + 1 ;
  • epyyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  • FIN BEPYYC ;
  •  
  • igau = 0 ;
  • REPETER BEPZZC nbgau ;
  • igau = igau + 1 ;
  • epzzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  • FIN BEPZZC ;
  •  
  • igau = 0 ;
  • REPETER BGAXYC nbgau ;
  • igau = igau + 1 ;
  • gaxyc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  • FIN BGAXYC ;
  •  
  • igau = 0 ;
  • REPETER BGAXZC nbgau ;
  • igau = igau + 1 ;
  • gaxzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  • FIN BGAXZC ;
  •  
  • igau = 0 ;
  • REPETER BGAYZC nbgau ;
  • igau = igau + 1 ;
  • gayzc.igau = EXTR (tabflup.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  • FIN BGAYZC ;
  •  
  • * Deformations inelastiques
  • igau = 0 ;
  • REPETER BEIXXC nbgau ;
  • igau = igau + 1 ;
  • eixxc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  • 'EIXX' 1 1 igau ;
  • FIN BEIXXC ;
  •  
  • igau = 0 ;
  • REPETER BEIYYC nbgau ;
  • igau = igau + 1 ;
  • eiyyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  • 'EIYY' 1 1 igau ;
  • FIN BEIYYC ;
  •  
  • igau = 0 ;
  • REPETER BEIZZC nbgau ;
  • igau = igau + 1 ;
  • eizzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  • 'EIZZ' 1 1 igau ;
  • FIN BEIZZC ;
  •  
  • igau = 0 ;
  • REPETER BGIXYC nbgau ;
  • igau = igau + 1 ;
  • gixyc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  • 'GIXY' 1 1 igau ;
  • FIN BGIXYC ;
  •  
  • igau = 0 ;
  • REPETER BGIXZC nbgau ;
  • igau = igau + 1 ;
  • gixzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  • 'GIXZ' 1 1 igau ;
  • FIN BGIXZC ;
  •  
  • igau = 0 ;
  • REPETER BGIYZC nbgau ;
  • igau = igau + 1 ;
  • giyzc.igau = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ietat)
  • 'GIYZ' 1 1 igau ;
  • FIN BGIYZC ;
  •  
  • * Contraintes
  • igau = 0 ;
  • REPETER BSMXXC nbgau ;
  • igau = igau + 1 ;
  • smxxc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  • FIN BSMXXC ;
  •  
  • igau = 0 ;
  • REPETER BSMYYC nbgau ;
  • igau = igau + 1 ;
  • smyyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  • FIN BSMYYC ;
  •  
  • igau = 0 ;
  • REPETER BSMZZC nbgau ;
  • igau = igau + 1 ;
  • smzzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  • FIN BSMZZC ;
  •  
  • igau = 0 ;
  • REPETER BSMXYC nbgau ;
  • igau = igau + 1 ;
  • smxyc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  • FIN BSMXYC ;
  •  
  • igau = 0 ;
  • REPETER BSMXZC nbgau ;
  • igau = igau + 1 ;
  • smxzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  • FIN BSMXZC ;
  •  
  • igau = 0 ;
  • REPETER BSMYZC nbgau ;
  • igau = igau + 1 ;
  • smyzc.igau = EXTR (tabflup.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  • FIN BSMYZC ;
  •  
  • * Variables internes
  • igau = 0 ;
  • REPETER BPXXC nbgau ;
  • igau = igau + 1 ;
  • pxxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'PXX ' 1 1 igau ;
  • FIN BPXXC ;
  •  
  • igau = 0 ;
  • REPETER BPYYC nbgau ;
  • igau = igau + 1 ;
  • pyyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'PYY ' 1 1 igau ;
  • FIN BPYYC ;
  •  
  • igau = 0 ;
  • REPETER BPZZC nbgau ;
  • igau = igau + 1 ;
  • pzzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'PZZ ' 1 1 igau ;
  • FIN BPZZC ;
  •  
  • igau = 0 ;
  • REPETER BPXYC nbgau ;
  • igau = igau + 1 ;
  • pxyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'PXY ' 1 1 igau ;
  • FIN BPXYC ;
  •  
  • igau = 0 ;
  • REPETER BPYZC nbgau ;
  • igau = igau + 1 ;
  • pyzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'PYZ ' 1 1 igau ;
  • FIN BPYZC ;
  •  
  • igau = 0 ;
  • REPETER BPZXC nbgau ;
  • igau = igau + 1 ;
  • pzxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'PZX ' 1 1 igau ;
  • FIN BPZXC ;
  •  
  • igau = 0 ;
  • REPETER BVNXXC nbgau ;
  • igau = igau + 1 ;
  • vnxxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'VNXX' 1 1 igau ;
  • FIN BVNXXC ;
  •  
  • igau = 0 ;
  • REPETER BVNYYC nbgau ;
  • igau = igau + 1 ;
  • vnyyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'VNYY' 1 1 igau ;
  • FIN BVNYYC ;
  •  
  • igau = 0 ;
  • REPETER BVNZZC nbgau ;
  • igau = igau + 1 ;
  • vnzzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'VNZZ' 1 1 igau ;
  • FIN BVNZZC ;
  •  
  • igau = 0 ;
  • REPETER BVNXYC nbgau ;
  • igau = igau + 1 ;
  • vnxyc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'VNXY' 1 1 igau ;
  • FIN BVNXYC ;
  •  
  • igau = 0 ;
  • REPETER BVNYZC nbgau ;
  • igau = igau + 1 ;
  • vnyzc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'VNYZ' 1 1 igau ;
  • FIN BVNYZC ;
  •  
  • igau = 0 ;
  • REPETER BVNZXC nbgau ;
  • igau = igau + 1 ;
  • vnzxc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'VNZX' 1 1 igau ;
  • FIN BVNZXC ;
  •  
  • igau = 0 ;
  • REPETER BEPCEC nbgau ;
  • igau = igau + 1 ;
  • epcec.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'EPCE' 1 1 igau ;
  • FIN BEPCEC ;
  •  
  • igau = 0 ;
  • REPETER BEPSHC nbgau ;
  • igau = igau + 1 ;
  • epshc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'EPSH' 1 1 igau ;
  • FIN BEPSHC ;
  •  
  • igau = 0 ;
  • REPETER BLLC nbgau ;
  • igau = igau + 1 ;
  • llc.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'LL ' 1 1 igau ;
  • FIN BLLC ;
  •  
  • igau = 0 ;
  • REPETER BEPSEC nbgau ;
  • igau = igau + 1 ;
  • epsec.igau = EXTR (tabflup.'VARIABLES_INTERNES'.ietat)
  • 'EPSE' 1 1 igau ;
  • FIN BEPSEC ;
  •  
  • * Extraction des resultats de calcul : Modele externe utilisateur------
  •  
  • * Deplacements
  • uxu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p01 ;
  • uxu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p02 ;
  • uxu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p03 ;
  • uxu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p04 ;
  • uxu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p10 ;
  • uxu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p20 ;
  • uxu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p30 ;
  • uxu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UX' p40 ;
  •  
  • uyu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p01 ;
  • uyu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p02 ;
  • uyu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p03 ;
  • uyu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p04 ;
  • uyu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p10 ;
  • uyu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p20 ;
  • uyu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p30 ;
  • uyu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UY' p40 ;
  •  
  • uzu.1 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p01 ;
  • uzu.2 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p02 ;
  • uzu.3 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p03 ;
  • uzu.4 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p04 ;
  • uzu.5 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p10 ;
  • uzu.6 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p20 ;
  • uzu.7 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p30 ;
  • uzu.8 = EXTR (tabutil.'DEPLACEMENTS'.ietat) 'UZ' p40 ;
  •  
  • * Deformations totales
  • igau = 0 ;
  • REPETER BEPXXU nbgau ;
  • igau = igau + 1 ;
  • epxxu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  • FIN BEPXXU ;
  •  
  • igau = 0 ;
  • REPETER BEPYYU nbgau ;
  • igau = igau + 1 ;
  • epyyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  • FIN BEPYYU ;
  •  
  • igau = 0 ;
  • REPETER BEPZZU nbgau ;
  • igau = igau + 1 ;
  • epzzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  • FIN BEPZZU ;
  •  
  • igau = 0 ;
  • REPETER BGAXYU nbgau ;
  • igau = igau + 1 ;
  • gaxyu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  • FIN BGAXYU ;
  •  
  • igau = 0 ;
  • REPETER BGAXZU nbgau ;
  • igau = igau + 1 ;
  • gaxzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  • FIN BGAXZU ;
  •  
  • igau = 0 ;
  • REPETER BGAYZU nbgau ;
  • igau = igau + 1 ;
  • gayzu.igau = EXTR (tabutil.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  • FIN BGAYZU ;
  •  
  • * Deformations inelastiques (variables internes supplementaires)
  • igau = 0 ;
  • REPETER BENXXU nbgau ;
  • igau = igau + 1 ;
  • enxxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'ENXX' 1 1 igau ;
  • FIN BENXXU ;
  •  
  • igau = 0 ;
  • REPETER BENYYU nbgau ;
  • igau = igau + 1 ;
  • enyyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'ENYY' 1 1 igau ;
  • FIN BENYYU ;
  •  
  • igau = 0 ;
  • REPETER BENZZU nbgau ;
  • igau = igau + 1 ;
  • enzzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'ENZZ' 1 1 igau ;
  • FIN BENZZU ;
  •  
  • igau = 0 ;
  • REPETER BGNXYU nbgau ;
  • igau = igau + 1 ;
  • gnxyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'GNXY' 1 1 igau ;
  • FIN BGNXYU ;
  •  
  • igau = 0 ;
  • REPETER BGNXZU nbgau ;
  • igau = igau + 1 ;
  • gnxzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'GNXZ' 1 1 igau ;
  • FIN BGNXZU ;
  •  
  • igau = 0 ;
  • REPETER BGNYZU nbgau ;
  • igau = igau + 1 ;
  • gnyzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'GNYZ' 1 1 igau ;
  • FIN BGNYZU ;
  •  
  • * Contraintes
  • igau = 0 ;
  • REPETER BSMXXU nbgau ;
  • igau = igau + 1 ;
  • smxxu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  • FIN BSMXXU ;
  •  
  • igau = 0 ;
  • REPETER BSMYYU nbgau ;
  • igau = igau + 1 ;
  • smyyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  • FIN BSMYYU ;
  •  
  • igau = 0 ;
  • REPETER BSMZZU nbgau ;
  • igau = igau + 1 ;
  • smzzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  • FIN BSMZZU ;
  •  
  • igau = 0 ;
  • REPETER BSMXYU nbgau ;
  • igau = igau + 1 ;
  • smxyu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  • FIN BSMXYU ;
  •  
  • igau = 0 ;
  • REPETER BSMXZU nbgau ;
  • igau = igau + 1 ;
  • smxzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  • FIN BSMXZU ;
  •  
  • igau = 0 ;
  • REPETER BSMYZU nbgau ;
  • igau = igau + 1 ;
  • smyzu.igau = EXTR (tabutil.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  • FIN BSMYZU ;
  •  
  • * Variables internes
  • igau = 0 ;
  • REPETER BPXXU nbgau ;
  • igau = igau + 1 ;
  • pxxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'PXX ' 1 1 igau ;
  • FIN BPXXU ;
  •  
  • igau = 0 ;
  • REPETER BPYYU nbgau ;
  • igau = igau + 1 ;
  • pyyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'PYY ' 1 1 igau ;
  • FIN BPYYU ;
  •  
  • igau = 0 ;
  • REPETER BPZZU nbgau ;
  • igau = igau + 1 ;
  • pzzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'PZZ ' 1 1 igau ;
  • FIN BPZZU ;
  •  
  • igau = 0 ;
  • REPETER BPXYU nbgau ;
  • igau = igau + 1 ;
  • pxyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'PXY ' 1 1 igau ;
  • FIN BPXYU ;
  •  
  • igau = 0 ;
  • REPETER BPYZU nbgau ;
  • igau = igau + 1 ;
  • pyzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'PYZ ' 1 1 igau ;
  • FIN BPYZU ;
  •  
  • igau = 0 ;
  • REPETER BPZXU nbgau ;
  • igau = igau + 1 ;
  • pzxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'PZX ' 1 1 igau ;
  • FIN BPZXU ;
  •  
  • igau = 0 ;
  • REPETER BVNXXU nbgau ;
  • igau = igau + 1 ;
  • vnxxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'VNXX' 1 1 igau ;
  • FIN BVNXXU ;
  •  
  • igau = 0 ;
  • REPETER BVNYYU nbgau ;
  • igau = igau + 1 ;
  • vnyyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'VNYY' 1 1 igau ;
  • FIN BVNYYU ;
  •  
  • igau = 0 ;
  • REPETER BVNZZU nbgau ;
  • igau = igau + 1 ;
  • vnzzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'VNZZ' 1 1 igau ;
  • FIN BVNZZU ;
  •  
  • igau = 0 ;
  • REPETER BVNXYU nbgau ;
  • igau = igau + 1 ;
  • vnxyu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'VNXY' 1 1 igau ;
  • FIN BVNXYU ;
  •  
  • igau = 0 ;
  • REPETER BVNYZU nbgau ;
  • igau = igau + 1 ;
  • vnyzu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'VNYZ' 1 1 igau ;
  • FIN BVNYZU ;
  •  
  • igau = 0 ;
  • REPETER BVNZXU nbgau ;
  • igau = igau + 1 ;
  • vnzxu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'VNZX' 1 1 igau ;
  • FIN BVNZXU ;
  •  
  • igau = 0 ;
  • REPETER BEPCEU nbgau ;
  • igau = igau + 1 ;
  • epceu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'EPCE' 1 1 igau ;
  • FIN BEPCEU ;
  •  
  • igau = 0 ;
  • REPETER BEPSHU nbgau ;
  • igau = igau + 1 ;
  • epshu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'EPSH' 1 1 igau ;
  • FIN BEPSHU ;
  •  
  • igau = 0 ;
  • REPETER BLLU nbgau ;
  • igau = igau + 1 ;
  • llu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'LL ' 1 1 igau ;
  • FIN BLLU ;
  •  
  • igau = 0 ;
  • REPETER BEPSEU nbgau ;
  • igau = igau + 1 ;
  • epseu.igau = EXTR (tabutil.'VARIABLES_INTERNES'.ietat)
  • 'EPSE' 1 1 igau ;
  • FIN BEPSEU ;
  •  
  • * Comparaison des deplacements-----------------------------------------
  •  
  • icomp = 1 ;
  • ino = 0 ;
  • REPETER BUX nbno ;
  • ino = ino + 1 ;
  • SI (EGA uxc.ino 0.0) ;
  • SI (NON (EGA uxu.ino 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' uxu.ino
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dux = ABS (uxc.ino - uxu.ino) ;
  • errrel = ABS (dux / uxc.ino) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message 'uxc.ino =' uxc.ino 'uxu.ino =' uxu.ino ;
  • Message ;
  • FINSI ;
  • ************ierr = 11 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BUX ;
  •  
  • icomp = 2 ;
  • ino = 0 ;
  • REPETER BUY nbno ;
  • ino = ino + 1 ;
  • SI (EGA uyc.ino 0.0) ;
  • SI (NON (EGA uyu.ino 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' uyu.ino
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • duy = ABS (uyc.ino - uyu.ino) ;
  • errrel = ABS (duy / uyc.ino) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message 'uyc.ino =' uyc.ino 'uyu.ino =' uyu.ino ;
  • Message ;
  • FINSI ;
  • ************ierr = 11 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BUY ;
  •  
  • icomp = 3 ;
  • ino = 0 ;
  • REPETER BUZ nbno ;
  • ino = ino + 1 ;
  • SI (EGA uzc.ino 0.0) ;
  • SI (NON (EGA uzu.ino 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' uzu.ino
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • duz = ABS (uzc.ino - uzu.ino) ;
  • errrel = ABS (duz / uzc.ino) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • Message 'uzc.ino =' uzc.ino 'uzu.ino =' uzu.ino ;
  • Message ;
  • FINSI ;
  • ************ierr = 11 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BUZ ;
  •  
  • * Comparaison des deformations-----------------------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BEPXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epxxc.igau 0.0) ;
  • SI (NON (EGA epxxu.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' epxxu.igau
  • ************ 'Composante ' (EXTR LDEFO icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depxx = ABS (epxxc.igau - epxxu.igau) ;
  • errrel = ABS (depxx / epxxc.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LDEFO icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'epxxc.igau =' epxxc.igau 'epxxu.igau =' epxxu.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPXX ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BEPYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epyyc.igau 0.0) ;
  • SI (NON (EGA epyyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epyyu.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depyy = ABS (epyyc.igau - epyyu.igau) ;
  • errrel = ABS (depyy / epyyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epyyc.igau =' epyyc.igau 'epyyu.igau =' epyyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPYY ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BEPZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epzzc.igau 0.0) ;
  • SI (NON (EGA epzzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epzzu.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = epzzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depzz = ABS (epzzc.igau - epzzu.igau) ;
  • errrel = ABS (depzz / epzzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epzzc.igau =' epzzc.igau 'epzzu.igau =' epzzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPZZ ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BGAXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gaxyc.igau 0.0) ;
  • SI (NON (EGA gaxyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gaxyu.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgaxy = ABS (gaxyc.igau - gaxyu.igau) ;
  • errrel = ABS (dgaxy / gaxyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gaxyc.igau =' gaxyc.igau 'gaxyu.igau =' gaxyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGAXY ;
  •  
  • icomp = 5 ;
  • igau = 0 ;
  • REPETER BGAXZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gaxzc.igau 0.0) ;
  • SI (NON (EGA gaxzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gaxzu.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgaxz = ABS (gaxzc.igau - gaxzu.igau) ;
  • errrel = ABS (dgaxz / gaxzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gaxzc.igau =' gaxzc.igau 'gaxzu.igau =' gaxzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGAXZ ;
  •  
  • icomp = 6 ;
  • igau = 0 ;
  • REPETER BGAYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gayzc.igau 0.0) ;
  • SI (NON (EGA gayzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gayzu.igau
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgayz = ABS (gayzc.igau - gayzu.igau) ;
  • errrel = ABS (dgayz / gayzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gayzc.igau =' gayzc.igau 'gayzu.igau =' gayzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 21 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGAYZ ;
  •  
  • * Comparaison des deformations inelastiques----------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BEIXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA eixxc.igau 0.0) ;
  • SI (NN(EGA enxxu.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' enxxu.igau
  • ************ 'Composante ' (EXTR LDEIN icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 30 ; errabs = enxxu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • deixx = ABS (eixxc.igau - enxxu.igau) ;
  • errrel = ABS (deixx / eixxc.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LDEIN icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'eixxc.igau =' eixxc.igau 'enxxu.igau =' enxxu.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEIXX ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BEIYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA eiyyc.igau 0.0) ;
  • SI (NON (EGA enyyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' enyyu.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = enyyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • deiyy = ABS (eiyyc.igau - enyyu.igau) ;
  • errrel = ABS (deiyy / eiyyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'eiyyc.igau =' eiyyc.igau 'enyyu.igau =' enyyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEIYY ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BEIZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA eizzc.igau 0.0) ;
  • SI (NON (EGA enzzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' enzzu.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = enzzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • deizz = ABS (eizzc.igau - enzzu.igau) ;
  • errrel = ABS (deizz / eizzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'eizzc.igau =' eizzc.igau 'enzzu.igau =' enzzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEIZZ ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BGIXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gixyc.igau 0.0) ;
  • SI (NON (EGA gnxyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gnxyu.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = gnxyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgixy = ABS (gixyc.igau - gnxyu.igau) ;
  • errrel = ABS (dgixy / gixyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gixyc.igau =' gixyc.igau 'gnxyu.igau =' gnxyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGIXY ;
  •  
  • icomp = 5 ;
  • igau = 0 ;
  • REPETER BGIXZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA gixzc.igau 0.0) ;
  • SI (NON (EGA gnxzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gnxzu.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = gnxzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgixz = ABS (gixzc.igau - gnxzu.igau) ;
  • errrel = ABS (dgixz / gixzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'gixzc.igau =' gixzc.igau 'gnxzu.igau =' gnxzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGIXZ ;
  •  
  • icomp = 6 ;
  • igau = 0 ;
  • REPETER BGIYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA giyzc.igau 0.0) ;
  • SI (NON (EGA gnyzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' gnyzu.igau
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 30 ; errabs = gnyzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dgiyz = ABS (giyzc.igau - gnyzu.igau) ;
  • errrel = ABS (dgiyz / giyzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'giyzc.igau =' giyzc.igau 'gnyzu.igau =' gnyzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 31 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BGIYZ ;
  •  
  • * Comparaison des contraintes------------------------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BSMXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smxxc.igau 0.0) ;
  • SI (NON (EGA smxxu.igau 0.0)) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur absolue' smxxu.igau
  • ************ 'Composante ' (EXTR LCONT icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message ;
  • ************FINSI ;
  • ierr = 40 ; errabs = smxxu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmxx = ABS (smxxc.igau - smxxu.igau) ;
  • errrel = ABS (dsmxx / smxxc.igau) ;
  • SI (errrel > critere) ;
  • ************SI (EGA IMESS 1) ;
  • ************Message 'Erreur relative' errrel
  • ************ 'Composante ' (EXTR LCONT icomp)
  • ************ ' Point Gauss' igau 'Etat' ietat ;
  • ************Message 'smxxc.igau =' smxxc.igau 'smxxu.igau =' smxxu.igau;
  • ************Message ;
  • ************FINSI ;
  • ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMXX ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BSMYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smyyc.igau 0.0) ;
  • SI (NON (EGA smyyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smyyu.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smyyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmyy = ABS (smyyc.igau - smyyu.igau) ;
  • errrel = ABS (dsmyy / smyyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smyyc.igau =' smyyc.igau 'smyyu.igau =' smyyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMYY ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BSMZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smzzc.igau 0.0) ;
  • SI (NON (EGA smzzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smzzu.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smzzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmzz = ABS (smzzc.igau - smzzu.igau) ;
  • errrel = ABS (dsmzz / smzzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smzzc.igau =' smzzc.igau 'smzzu.igau =' smzzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMZZ ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BSMXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smxyc.igau 0.0) ;
  • SI (NON (EGA smxyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smxyu.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smxyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmxy = ABS (smxyc.igau - smxyu.igau) ;
  • errrel = ABS (dsmxy / smxyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smxyc.igau =' smxyc.igau 'smxyu.igau =' smxyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMXY ;
  •  
  • icomp = 5 ;
  • igau = 0 ;
  • REPETER BSMXZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smxzc.igau 0.0) ;
  • SI (NON (EGA smxzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smxzu.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smxzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmxz = ABS (smxzc.igau - smxzu.igau) ;
  • errrel = ABS (dsmxz / smxzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smxzc.igau =' smxzc.igau 'smxzu.igau =' smxzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMXZ ;
  •  
  • icomp = 6 ;
  • igau = 0 ;
  • REPETER BSMYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA smyzc.igau 0.0) ;
  • SI (NON (EGA smyzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' smyzu.igau
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 40 ; errabs = smyzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dsmyz = ABS (smyzc.igau - smyzu.igau) ;
  • errrel = ABS (dsmyz / smyzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'smyzc.igau =' smyzc.igau 'smyzu.igau =' smyzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 41 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BSMYZ ;
  •  
  • * Comparaison des variables internes-----------------------------------
  •  
  • icomp = 1 ;
  • igau = 0 ;
  • REPETER BPXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pxxc.igau 0.0) ;
  • SI (NON (EGA pxxu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' pxxu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = pxxu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dpxx = ABS (pxxc.igau - pxxu.igau) ;
  • errrel = ABS (dpxx / pxxc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'pxxc.igau =' pxxc.igau 'pxxu.igau =' pxxu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BPXX ;
  •  
  • icomp = 2 ;
  • igau = 0 ;
  • REPETER BPYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pyyc.igau 0.0) ;
  • SI (NON (EGA pyyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' pyyu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = pyyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dpyy = ABS (pyyc.igau - pyyu.igau) ;
  • errrel = ABS (dpyy / pyyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'pyyc.igau =' pyyc.igau 'pyyu.igau =' pyyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BPYY ;
  •  
  • icomp = 3 ;
  • igau = 0 ;
  • REPETER BPZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pzzc.igau 0.0) ;
  • SI (NON (EGA pzzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' pzzu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = pzzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dpzz = ABS (pzzc.igau - pzzu.igau) ;
  • errrel = ABS (dpzz / pzzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'pzzc.igau =' pzzc.igau 'pzzu.igau =' pzzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BPZZ ;
  •  
  • icomp = 4 ;
  • igau = 0 ;
  • REPETER BPXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pxyc.igau 0.0) ;
  • SI (NON (EGA pxyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' pxyu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = pxyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dpxy = ABS (pxyc.igau - pxyu.igau) ;
  • errrel = ABS (dpxy / pxyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'pxyc.igau =' pxyc.igau 'pxyu.igau =' pxyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BPXY ;
  •  
  • icomp = 5 ;
  • igau = 0 ;
  • REPETER BPYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pyzc.igau 0.0) ;
  • SI (NON (EGA pyzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' pyzu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = pyzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dpyz = ABS (pyzc.igau - pyzu.igau) ;
  • errrel = ABS (dpyz / pyzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'pyzc.igau =' pyzc.igau 'pyzu.igau =' pyzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BPYZ ;
  •  
  • icomp = 6 ;
  • igau = 0 ;
  • REPETER BPZX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA pzxc.igau 0.0) ;
  • SI (NON (EGA pzxu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' pzxu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = pzxu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dpzx = ABS (pzxc.igau - pzxu.igau) ;
  • errrel = ABS (dpzx / pzxc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'pzxc.igau =' pzxc.igau 'pzxu.igau =' pzxu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BPZX ;
  •  
  • icomp = 7 ;
  • igau = 0 ;
  • REPETER BVNXX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA vnxxc.igau 0.0) ;
  • SI (NON (EGA vnxxu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' vnxxu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = vnxxu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dvnxx = ABS (vnxxc.igau - vnxxu.igau) ;
  • errrel = ABS (dvnxx / vnxxc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'vnxxc.igau =' vnxxc.igau 'vnxxu.igau =' vnxxu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BVNXX ;
  •  
  • icomp = 8 ;
  • igau = 0 ;
  • REPETER BVNYY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA vnyyc.igau 0.0) ;
  • SI (NON (EGA vnyyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' vnyyu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = vnyyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dvnyy = ABS (vnyyc.igau - vnyyu.igau) ;
  • errrel = ABS (dvnyy / vnyyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'vnyyc.igau =' vnyyc.igau 'vnyyu.igau =' vnyyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BVNYY ;
  •  
  • icomp = 9 ;
  • igau = 0 ;
  • REPETER BVNZZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA vnzzc.igau 0.0) ;
  • SI (NON (EGA vnzzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' vnzzu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = vnzzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dvnzz = ABS (vnzzc.igau - vnzzu.igau) ;
  • errrel = ABS (dvnzz / vnzzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'vnzzc.igau =' vnzzc.igau 'vnzzu.igau =' vnzzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BVNZZ ;
  •  
  • icomp = 10 ;
  • igau = 0 ;
  • REPETER BVNXY nbgau ;
  • igau = igau + 1 ;
  • SI (EGA vnxyc.igau 0.0) ;
  • SI (NON (EGA vnxyu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' vnxyu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = vnxyu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dvnxy = ABS (vnxyc.igau - vnxyu.igau) ;
  • errrel = ABS (dvnxy / vnxyc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'vnxyc.igau =' vnxyc.igau 'vnxyu.igau =' vnxyu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BVNXY ;
  •  
  • icomp = 11 ;
  • igau = 0 ;
  • REPETER BVNYZ nbgau ;
  • igau = igau + 1 ;
  • SI (EGA vnyzc.igau 0.0) ;
  • SI (NON (EGA vnyzu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' vnyzu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = vnyzu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dvnyz = ABS (vnyzc.igau - vnyzu.igau) ;
  • errrel = ABS (dvnyz / vnyzc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'vnyzc.igau =' vnyzc.igau 'vnyzu.igau =' vnyzu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BVNYZ ;
  •  
  • icomp = 12 ;
  • igau = 0 ;
  • REPETER BVNZX nbgau ;
  • igau = igau + 1 ;
  • SI (EGA vnzxc.igau 0.0) ;
  • SI (NON (EGA vnzxu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' vnzxu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = vnzxu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dvnzx = ABS (vnzxc.igau - vnzxu.igau) ;
  • errrel = ABS (dvnzx / vnzxc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'vnzxc.igau =' vnzxc.igau 'vnzxu.igau =' vnzxu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BVNZX ;
  •  
  • icomp = 13 ;
  • igau = 0 ;
  • REPETER BEPCE nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epcec.igau 0.0) ;
  • SI (NON (EGA epceu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epceu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = epceu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depce = ABS (epcec.igau - epceu.igau) ;
  • errrel = ABS (depce / epcec.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epcec.igau =' epcec.igau 'epceu.igau =' epceu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPCE ;
  •  
  • icomp = 14 ;
  • igau = 0 ;
  • REPETER BEPSH nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epshc.igau 0.0) ;
  • SI (NON (EGA epshu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epshu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = epshu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depsh = ABS (epshc.igau - epshu.igau) ;
  • errrel = ABS (depsh / epshc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epshc.igau =' epshc.igau 'epshu.igau =' epshu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPSH ;
  •  
  • icomp = 15 ;
  • igau = 0 ;
  • REPETER BLL nbgau ;
  • igau = igau + 1 ;
  • SI (EGA llc.igau 0.0) ;
  • SI (NON (EGA llu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' llu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = llu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • dll = ABS (llc.igau - llu.igau) ;
  • errrel = ABS (dll / llc.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'llc.igau =' llc.igau 'llu.igau =' llu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BLL ;
  •  
  • icomp = 16 ;
  • igau = 0 ;
  • REPETER BEPSE nbgau ;
  • igau = igau + 1 ;
  • SI (EGA epsec.igau 0.0) ;
  • SI (NON (EGA epseu.igau 0.0)) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur absolue' epseu.igau
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message ;
  • FINSI ;
  • ************ierr = 50 ; errabs = epseu.igau ; QUITTER BECART ;
  • FINSI ;
  • SINON ;
  • depse = ABS (epsec.igau - epseu.igau) ;
  • errrel = ABS (depse / epsec.igau) ;
  • SI (errrel > critere) ;
  • SI (EGA IMESS 1) ;
  • Message 'Erreur relative' errrel
  • 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • Message 'epsec.igau =' epsec.igau 'epseu.igau =' epseu.igau;
  • Message ;
  • FINSI ;
  • ************ierr = 51 ; QUITTER BECART ;
  • FINSI ;
  • FINSI ;
  • FIN BEPSE ;
  •  
  • FIN BECART ;
  •  
  • * FIN BOUCLE COMPARAISON DES RESULTATS________________________________
  •  
  • * MESSAGES D'ERREUR_____________________________________________________
  •  
  • Message ;
  • Message 'CPU avec modele interne = ' CPUint 'centiemes de seconde' ;
  • Message 'CPU avec modele externe = ' CPUext 'centiemes de seconde' ;
  •  
  • ifin = netat-1 ;
  • eixxcmax = 0.0 ;
  • igau = 0 ;
  • REPETER BEIXXFC nbgau ;
  • igau = igau + 1 ;
  • eixx = EXTR (tabflup.'DEFORMATIONS_INELASTIQUES'.ifin)
  • 'EIXX' 1 1 igau ;
  • SI (eixx > eixxcmax) ; eixxcmax = eixx ; FINSI ;
  • FIN BEIXXFC ;
  • eixxumax = 0.0 ;
  • igau = 0 ;
  • REPETER BEIXXFU nbgau ;
  • igau = igau + 1 ;
  • eixx = EXTR (tabutil.'VARIABLES_INTERNES'.ifin)
  • 'ENXX' 1 1 igau ;
  • SI (eixx > eixxumax) ; eixxumax = eixx ; FINSI ;
  • FIN BEIXXFU ;
  • Message ;
  • Message 'Deformation de fluage finale (maxi sur points de Gauss)' ;
  • Message 'EIXX max avec modele interne = ' eixxcmax ;
  • Message 'EIXX max avec modele externe = ' eixxumax ;
  • Message ;
  •  
  • SI (ierr > 0) ;
  •  
  • SI (EGA ierr 1) ;
  • Message 'Instants de calcul non coincidants' ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 10) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 11) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  • ' Noeud' ino 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 20) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 21) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 30) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 31) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LDEIN icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 40) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 41) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 50) ;
  • Message 'Erreur absolue' errabs 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SI (EGA ierr 51) ;
  • Message 'Erreur relative' errrel 'Composante ' (EXTR LVARI icomp)
  • ' Point Gauss' igau 'Etat' ietat ;
  • ERRE 5 ;
  • FINSI ;
  •  
  • SINON ;
  •  
  • Message 'Resultats concordants entre les deux modeles' ;
  • Message 'Aucune erreur relative superieure a ' critere
  • ' pour les composantes' ;
  • Message (EXTR LDEFO idefo) ' ' (EXTR LDEIN idein)
  • ' ' (EXTR LCONT icont) ;
  • ERRE 0 ;
  •  
  • FINSI ;
  •  
  • FIN ;
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  •  
  • © Cast3M 2003 - Tous droits réservés.
    Mentions légales