Télécharger la plaquette Cast3M


fr en

Télécharger umat01.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  umat01.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 No1 : ELASTICITE ISOTROPE ISOTHERME
  15. *
  16. *=======================================================================
  17.  
  18.  
  19. opti echo 0 ;
  20.  
  21.  
  22. *=======================================================================
  23. *
  24. *                       PROCEDURES DEPENDANTES
  25. *
  26. *=======================================================================
  27.  
  28.  
  29.  
  30.  
  31. *=======================================================================
  32. *
  33. *                      JEU DE DONNEES PRINCIPAL
  34. *
  35. *=======================================================================
  36.  
  37.  
  38. *=======================================================================
  39. * Maillage
  40.  
  41. opti dime 3 elem cub8 ;
  42.  
  43. long = 1.0 ;
  44. larg = 1.0 ;
  45. epai = 1.0 ;
  46.  
  47. p1  = 0.0 0.0 0.0 ;
  48. p2  = p1 PLUS (larg 0.0 0.0) ;
  49. p3  = p2 PLUS (0.0 long 0.0) ;
  50. p4  = p1 PLUS (0.0 long 0.0) ;
  51.  
  52. p11 = p1 PLUS (0.0 0.0 epai) ;
  53. p21 = p2 PLUS (0.0 0.0 epai) ;
  54. p31 = p3 PLUS (0.0 0.0 epai) ;
  55. p41 = p4 PLUS (0.0 0.0 epai) ;
  56.  
  57. p1p2 = DROIT 1 p1 p2 ;
  58. p4p3 = DROIT 1 p4 p3 ;
  59.  
  60. sbas = p1p2 REGL 1 p4p3 ;
  61.  
  62. p11p21 = DROIT 1 p11 p21 ;
  63. p41p31 = DROIT 1 p41 p31 ;
  64.  
  65. shaut = p11p21 REGL 1 p41p31 ;
  66.  
  67. cube = sbas VOLU 1 shaut ;
  68.  
  69.  
  70. *=======================================================================
  71. * Comportement elastique isotrope isotherme
  72. * Traction uniaxiale Z
  73.  
  74. *-----------------------------------------------------------------------
  75. * 1 - Calcul avec modele utilisateur
  76.  
  77. * 1.1 Definition du modele utilisateur
  78.  
  79. LCMAT = MOTS 'YOUN' 'NU  ' ;
  80. mo_meca = MODE cube 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  81.                     'NON_LINEAIRE' 'UTILISATEUR'
  82.                     'NUME_LOI' 10 'C_MATERIAU' LCMAT ;
  83.  
  84. * 1.2 Affectation des proprietes de materiau
  85.  
  86. xyoun = 2.1e+11 ;
  87. xnu = 0.3 ;
  88. ma_meca = MATE mo_meca 'YOUN' xyoun 'NU' xnu ;
  89.  
  90. * 1.3 Evaluation de coefficients de la matrice de Hooke en vue du
  91. *     calcul de la solution theorique
  92.  
  93. coef1 = xyoun / ((1.0+xnu) * (1.0-(2.0*xnu))) ;
  94. coezz = coef1 * (1.0-xnu) ;
  95. coexxyy = xnu * coef1 ;
  96.  
  97. * 1.4 Definition des chargements
  98.  
  99. * 1.4.0 Creation de progressions pour la ponderation des chargements,
  100. *       definissant autant d'etats a calculer
  101.  
  102. L_temps = PROG 0.0 0.5 1.0 ;
  103. L_depl  = PROG 0.0 0.5 1.0 ;
  104.  
  105. netat = DIME L_temps ;
  106.  
  107. * 1.4.1 Blocages mecaniques et deplacement axial impose
  108.  
  109. uz3 = 1.0e-2 ;
  110. clim1 = BLOQ sbas 'DEPL' ;
  111. clim2 = BLOQ shaut 'UX' 'UY' ;
  112. clim3 = BLOQ shaut 'UZ' ;
  113. depi3 = DEPI clim3 uz3 ;
  114. charmeca = CHAR 'DIMP' (EVOL MANU L_temps L_depl) depi3 ;
  115.  
  116. * 1.5 Resolution par PASAPAS
  117.  
  118. TABU = TABLE ;
  119. TABU.'MODELE' = mo_meca ;
  120. TABU.'CARACTERISTIQUES' = ma_meca ;
  121. TABU.'CHARGEMENT' = charmeca ;
  122. TABU.'BLOCAGES_MECANIQUES' = clim1 ET clim2 ET clim3 ;
  123. TABU.'TEMPS_CALCULES' = L_temps ;
  124. TMASAU=table;
  125. tabu . 'MES_SAUVEGARDES'=TMASAU;
  126. TMASAU .'DEFTO'=VRAI;
  127. TMASAU .'DEFIN'=VRAI;
  128. TABU.'HYPOTHESE_DEFORMATIONS'='LINEAIRE';
  129. PASAPAS TABU ;
  130.  
  131. *-----------------------------------------------------------------------
  132. * 2 - Comparaison des resultats obtenus avec la solution theorique
  133.  
  134. * N.B. On simule un essai de traction ZZ
  135. *      - Critere strict sur UZ  , EPZZ, SMXX, SMYY et SMZZ :
  136. *        errrel < 1.d-12
  137. *      - Pour les autres composantes on note les depassements du critere
  138. *        si impression demandee
  139.  
  140. uxu = TABLE ; uyu = TABLE ; uzu = TABLE ;
  141. epxxu = TABLE ; epyyu = TABLE ; epzzu = TABLE ;
  142. gaxyu = TABLE ; gaxzu = TABLE ; gayzu = TABLE ;
  143. smxxu = TABLE ; smyyu = TABLE ; smzzu = TABLE ;
  144. smxyu = TABLE ; smxzu = TABLE ; smyzu = TABLE ;
  145.  
  146. LDEPL = MOTS 'UX  ' 'UY  ' 'UZ  ' ;
  147. LDEFO = MOTS 'EPXX' 'EPYY' 'EPZZ' 'GAXY' 'GAXZ' 'GAYZ' ;
  148. LCONT = MOTS 'SMXX' 'SMYY' 'SMZZ' 'SMXY' 'SMXZ' 'SMYZ' ;
  149.  
  150. idepl = 3 ;
  151. idefo = 3 ;
  152. icon1 = 1 ;
  153. icon2 = 2 ;
  154. icon3 = 3 ;
  155.  
  156. ierr = 0 ;
  157. critere = 1.d-12 ;
  158. IMESS = 0 ;
  159. ***IMESS = 1 ;
  160.  
  161. nbno  = 8 ;
  162. nbgau = 8 ;
  163.  
  164. ietat = 0 ;
  165.  
  166. * DEBUT BOUCLE COMPARAISON DES RESULTATS________________________________
  167.  
  168. REPETER BECART (netat-1) ;
  169.  
  170.    ietat = ietat + 1 ;
  171.    temps = TABU.'TEMPS'.ietat ;
  172.  
  173. *  Calcul de la solution theorique--------------------------------------
  174.  
  175.    uz = temps * uz3 ;
  176.    epzz = temps * (uz3 / epai) ;
  177.    smxx = coexxyy * epzz ;
  178.    smyy = coexxyy * epzz ;
  179.    smzz = coezz * epzz ;
  180.  
  181. *  Extraction des resultats de calcul : Modele externe utilisateur------
  182.  
  183. *  Deplacements
  184.    uxu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p1 ;
  185.    uxu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p2 ;
  186.    uxu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p3 ;
  187.    uxu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p4 ;
  188.    uxu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p11 ;
  189.    uxu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p21 ;
  190.    uxu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p31 ;
  191.    uxu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UX' p41 ;
  192.  
  193.    uyu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p1 ;
  194.    uyu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p2 ;
  195.    uyu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p3 ;
  196.    uyu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p4 ;
  197.    uyu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p11 ;
  198.    uyu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p21 ;
  199.    uyu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p31 ;
  200.    uyu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UY' p41 ;
  201.  
  202.    uzu.1 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p1 ;
  203.    uzu.2 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p2 ;
  204.    uzu.3 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p3 ;
  205.    uzu.4 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p4 ;
  206.    uzu.5 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p11 ;
  207.    uzu.6 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p21 ;
  208.    uzu.7 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p31 ;
  209.    uzu.8 = EXTR (TABU.'DEPLACEMENTS'.ietat) 'UZ' p41 ;
  210.  
  211. *  *  Deformations totales
  212.    igau = 0 ;
  213.    REPETER BEPXXU nbgau ;
  214.       igau = igau + 1 ;
  215.       epxxu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPXX' 1 1 igau ;
  216.    FIN BEPXXU ;
  217.  
  218.    igau = 0 ;
  219.    REPETER BEPYYU nbgau ;
  220.       igau = igau + 1 ;
  221.       epyyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPYY' 1 1 igau ;
  222.    FIN BEPYYU ;
  223.  
  224.    igau = 0 ;
  225.    REPETER BEPZZU nbgau ;
  226.       igau = igau + 1 ;
  227.       epzzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'EPZZ' 1 1 igau ;
  228.    FIN BEPZZU ;
  229.  
  230.    igau = 0 ;
  231.    REPETER BGAXYU nbgau ;
  232.       igau = igau + 1 ;
  233.       gaxyu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXY' 1 1 igau ;
  234.    FIN BGAXYU ;
  235.  
  236.    igau = 0 ;
  237.    REPETER BGAXZU nbgau ;
  238.       igau = igau + 1 ;
  239.       gaxzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAXZ' 1 1 igau ;
  240.    FIN BGAXZU ;
  241.  
  242.    igau = 0 ;
  243.    REPETER BGAYZU nbgau ;
  244.       igau = igau + 1 ;
  245.       gayzu.igau = EXTR (TABU.'DEFORMATIONS'.ietat) 'GAYZ' 1 1 igau ;
  246.    FIN BGAYZU ;
  247.  
  248. *  Contraintes
  249.    igau = 0 ;
  250.    REPETER BSMXXU nbgau ;
  251.       igau = igau + 1 ;
  252.       smxxu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXX' 1 1 igau ;
  253.    FIN BSMXXU ;
  254.  
  255.    igau = 0 ;
  256.    REPETER BSMYYU nbgau ;
  257.       igau = igau + 1 ;
  258.       smyyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYY' 1 1 igau ;
  259.    FIN BSMYYU ;
  260.  
  261.    igau = 0 ;
  262.    REPETER BSMZZU nbgau ;
  263.       igau = igau + 1 ;
  264.       smzzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMZZ' 1 1 igau ;
  265.    FIN BSMZZU ;
  266.  
  267.    igau = 0 ;
  268.    REPETER BSMXYU nbgau ;
  269.       igau = igau + 1 ;
  270.       smxyu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXY' 1 1 igau ;
  271.    FIN BSMXYU ;
  272.  
  273.    igau = 0 ;
  274.    REPETER BSMXZU nbgau ;
  275.       igau = igau + 1 ;
  276.       smxzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMXZ' 1 1 igau ;
  277.    FIN BSMXZU ;
  278.  
  279.    igau = 0 ;
  280.    REPETER BSMYZU nbgau ;
  281.       igau = igau + 1 ;
  282.       smyzu.igau = EXTR (TABU.'CONTRAINTES'.ietat) 'SMYZ' 1 1 igau ;
  283.    FIN BSMYZU ;
  284.  
  285. *  Comparaison des deplacements-----------------------------------------
  286.  
  287.    icomp = 1 ;
  288.    ino = 0 ;
  289.    REPETER BUX nbno ;
  290.       ino = ino + 1 ;
  291.       SI (NON (EGA uxu.ino 0.0)) ;
  292.          SI (EGA IMESS 1) ;
  293.          Message 'Erreur absolue' uxu.ino
  294.                  'Composante ' (EXTR LDEPL icomp)
  295.                  ' Noeud' ino 'Etat' ietat ;
  296.          Message ;
  297.          FINSI ;
  298. *********ierr = 10 ; errabs = uxu.ino ; QUITTER BECART ;
  299.       FINSI ;
  300.    FIN BUX ;
  301.  
  302.    icomp = 2 ;
  303.    ino = 0 ;
  304.    REPETER BUY nbno ;
  305.       ino = ino + 1 ;
  306.       SI (NON (EGA uyu.ino 0.0)) ;
  307.          SI (EGA IMESS 1) ;
  308.          Message 'Erreur absolue' uyu.ino
  309.                  'Composante ' (EXTR LDEPL icomp)
  310.                  ' Noeud' ino 'Etat' ietat ;
  311.          Message ;
  312.          FINSI ;
  313. *********ierr = 10 ; errabs = uyu.ino ; QUITTER BECART ;
  314.       FINSI ;
  315.    FIN BUY ;
  316.  
  317.    icomp = 3 ;
  318.    ino = 0 ;
  319.    REPETER BUZ nbno ;
  320.       ino = ino + 1 ;
  321.       SI (ino < 5) ;
  322.          SI (NON (EGA uzu.ino 0.0)) ;
  323.             SI (EGA IMESS 1) ;
  324.             Message 'Erreur absolue' uzu.ino
  325.                     'Composante ' (EXTR LDEPL icomp)
  326.                     ' Noeud' ino 'Etat' ietat ;
  327.             Message ;
  328.             FINSI ;
  329. ************ierr = 10 ; errabs = uzu.ino ; QUITTER BECART ;
  330.          FINSI ;
  331.       SINON ;
  332.          duz = ABS (uz - uzu.ino) ;
  333.          errrel = ABS (duz / uz) ;
  334.          SI (errrel > critere) ;
  335. ************SI (EGA IMESS 1) ;
  336. ************Message 'Erreur relative' errrel
  337. ************        'Composante ' (EXTR LDEPL icomp)
  338. ************        ' Noeud' ino 'Etat' ietat ;
  339. ************Message 'uz =' uz 'uzu.ino =' uzu.ino ;
  340. ************Message ;
  341. ************FINSI ;
  342.             ierr = 11 ; QUITTER BECART ;
  343.          FINSI ;
  344.       FINSI ;
  345.    FIN BUZ ;
  346.  
  347. *  Comparaison des deformations totales---------------------------------
  348.  
  349.    icomp = 1 ;
  350.    igau = 0 ;
  351.    REPETER BEPXX nbgau ;
  352.       igau = igau + 1 ;
  353.       SI (NON (EGA epxxu.igau 0.0)) ;
  354.          SI (EGA IMESS 1) ;
  355.          Message 'Erreur absolue' epxxu.igau
  356.                  'Composante ' (EXTR LDEFO icomp)
  357.                  ' Point Gauss' igau 'Etat' ietat ;
  358.          Message ;
  359.          FINSI ;
  360. *********ierr = 20 ; errabs = epxxu.igau ; QUITTER BECART ;
  361.       FINSI ;
  362.    FIN BEPXX ;
  363.  
  364.    icomp = 2 ;
  365.    igau = 0 ;
  366.    REPETER BEPYY nbgau ;
  367.       igau = igau + 1 ;
  368.       SI (NON (EGA epyyu.igau 0.0)) ;
  369.          SI (EGA IMESS 1) ;
  370.          Message 'Erreur absolue' epyyu.igau
  371.                  'Composante ' (EXTR LDEFO icomp)
  372.                  ' Point Gauss' igau 'Etat' ietat ;
  373.          Message ;
  374.          FINSI ;
  375. *********ierr = 20 ; errabs = epyyu.igau ; QUITTER BECART ;
  376.       FINSI ;
  377.    FIN BEPYY ;
  378.  
  379.    icomp = 3 ;
  380.    igau = 0 ;
  381.    REPETER BEPZZ nbgau ;
  382.       igau = igau + 1 ;
  383.       depzz = ABS (epzz - epzzu.igau) ;
  384.       errrel = ABS (depzz / epzz) ;
  385.       SI (errrel > critere) ;
  386. *********SI (EGA IMESS 1) ;
  387. *********Message 'Erreur relative' errrel
  388. *********        'Composante ' (EXTR LDEFO icomp)
  389. *********        ' Point Gauss' igau 'Etat' ietat ;
  390. *********Message 'epzz =' epzz 'epzzu.igau =' epzzu.igau ;
  391. *********Message ;
  392. *********FINSI ;
  393.          ierr = 21 ; QUITTER BECART ;
  394.       FINSI ;
  395.    FIN BEPZZ ;
  396.  
  397.    icomp = 4 ;
  398.    igau = 0 ;
  399.    REPETER BGAXY nbgau ;
  400.       igau = igau + 1 ;
  401.       SI (NON (EGA gaxyu.igau 0.0)) ;
  402.          SI (EGA IMESS 1) ;
  403.          Message 'Erreur absolue' gaxyu.igau
  404.                  'Composante ' (EXTR LDEFO icomp)
  405.                  ' Point Gauss' igau 'Etat' ietat ;
  406.          Message ;
  407.          FINSI ;
  408. *********ierr = 20 ; errabs = gaxyu.igau ; QUITTER BECART ;
  409.       FINSI ;
  410.    FIN BGAXY ;
  411.  
  412.    icomp = 5 ;
  413.    igau = 0 ;
  414.    REPETER BGAXZ nbgau ;
  415.       igau = igau + 1 ;
  416.       SI (NON (EGA gaxzu.igau 0.0)) ;
  417.          SI (EGA IMESS 1) ;
  418.          Message 'Erreur absolue' gaxzu.igau
  419.                  'Composante ' (EXTR LDEFO icomp)
  420.                  ' Point Gauss' igau 'Etat' ietat ;
  421.          Message ;
  422.          FINSI ;
  423. *********ierr = 20 ; errabs = gaxzu.igau ; QUITTER BECART ;
  424.       FINSI ;
  425.    FIN BGAXZ ;
  426.  
  427.    icomp = 6 ;
  428.    igau = 0 ;
  429.    REPETER BGAYZ nbgau ;
  430.       igau = igau + 1 ;
  431.       SI (NON (EGA gayzu.igau 0.0)) ;
  432.          SI (EGA IMESS 1) ;
  433.          Message 'Erreur absolue' gayzu.igau
  434.                  'Composante ' (EXTR LDEFO icomp)
  435.                  ' Point Gauss' igau 'Etat' ietat ;
  436.          Message ;
  437.          FINSI ;
  438. *********ierr = 20 ; errabs = gayzu.igau ; QUITTER BECART ;
  439.       FINSI ;
  440.    FIN BGAYZ ;
  441.  
  442. *  Comparaison des contraintes------------------------------------------
  443.  
  444.    icomp = 1 ;
  445.    igau = 0 ;
  446.    REPETER BSMXX nbgau ;
  447.       igau = igau + 1 ;
  448.       dsmxx = ABS (smxx - smxxu.igau) ;
  449.       errrel = ABS (dsmxx / smxx) ;
  450.       SI (errrel > critere) ;
  451. *********SI (EGA IMESS 1) ;
  452. *********Message 'Erreur relative' errrel
  453. *********        'Composante ' (EXTR LCONT icomp)
  454. *********        ' Point Gauss' igau 'Etat' ietat ;
  455. *********Message 'smxx =' smxx 'smxxu.igau =' smxxu.igau ;
  456. *********Message ;
  457. *********FINSI ;
  458.          ierr = 31 ; QUITTER BECART ;
  459.       FINSI ;
  460.    FIN BSMXX ;
  461.  
  462.    icomp = 2 ;
  463.    igau = 0 ;
  464.    REPETER BSMYY nbgau ;
  465.       igau = igau + 1 ;
  466.       dsmyy = ABS (smyy - smyyu.igau) ;
  467.       errrel = ABS (dsmyy / smyy) ;
  468.       SI (errrel > critere) ;
  469. *********SI (EGA IMESS 1) ;
  470. *********Message 'Erreur relative' errrel
  471. *********        'Composante ' (EXTR LCONT icomp)
  472. *********        ' Point Gauss' igau 'Etat' ietat ;
  473. *********Message 'smyy =' smyy 'smyyu.igau =' smyyu.igau ;
  474. *********Message ;
  475. *********FINSI ;
  476.          ierr = 31 ; QUITTER BECART ;
  477.       FINSI ;
  478.    FIN BSMYY ;
  479.  
  480.    icomp = 3 ;
  481.    igau = 0 ;
  482.    REPETER BSMZZ nbgau ;
  483.       igau = igau + 1 ;
  484.       dsmzz = ABS (smzz - smzzu.igau) ;
  485.       errrel = ABS (dsmzz / smzz) ;
  486.       SI (errrel > critere) ;
  487. *********SI (EGA IMESS 1) ;
  488. *********Message 'Erreur relative' errrel
  489. *********        'Composante ' (EXTR LCONT icomp)
  490. *********        ' Point Gauss' igau 'Etat' ietat ;
  491. *********Message 'smzz =' smzz 'smzzu.igau =' smzzu.igau ;
  492. *********Message ;
  493. *********FINSI ;
  494.          ierr = 31 ; QUITTER BECART ;
  495.       FINSI ;
  496.    FIN BSMZZ ;
  497.  
  498.    icomp = 4 ;
  499.    igau = 0 ;
  500.    REPETER BSMXY nbgau ;
  501.       igau = igau + 1 ;
  502.       SI (NON (EGA smxyu.igau 0.0)) ;
  503.          SI (EGA IMESS 1) ;
  504.          Message 'Erreur absolue' smxyu.igau
  505.                  'Composante ' (EXTR LCONT icomp)
  506.                  ' Point Gauss' igau 'Etat' ietat ;
  507.          Message ;
  508.          FINSI ;
  509. *********ierr = 30 ; errabs = smxyu.igau ; QUITTER BECART ;
  510.       FINSI ;
  511.    FIN BSMXY ;
  512.  
  513.    icomp = 5 ;
  514.    igau = 0 ;
  515.    REPETER BSMXZ nbgau ;
  516.       igau = igau + 1 ;
  517.       SI (NON (EGA smxzu.igau 0.0)) ;
  518.          SI (EGA IMESS 1) ;
  519.          Message 'Erreur absolue' smxzu.igau
  520.                  'Composante ' (EXTR LCONT icomp)
  521.                  ' Point Gauss' igau 'Etat' ietat ;
  522.          Message ;
  523.          FINSI ;
  524. *********ierr = 30 ; errabs = smxzu.igau ; QUITTER BECART ;
  525.       FINSI ;
  526.    FIN BSMXZ ;
  527.  
  528.    icomp = 6 ;
  529.    igau = 0 ;
  530.    REPETER BSMYZ nbgau ;
  531.       igau = igau + 1 ;
  532.       SI (NON (EGA smyzu.igau 0.0)) ;
  533.          SI (EGA IMESS 1) ;
  534.          Message 'Erreur absolue' smyzu.igau
  535.                  'Composante ' (EXTR LCONT icomp)
  536.                  ' Point Gauss' igau 'Etat' ietat ;
  537.          Message ;
  538.          FINSI ;
  539. *********ierr = 30 ; errabs = smyzu.igau ; QUITTER BECART ;
  540.       FINSI ;
  541.    FIN BSMYZ ;
  542.  
  543. FIN BECART ;
  544.  
  545. * FIN   BOUCLE COMPARAISON DES RESULTATS________________________________
  546.  
  547. * MESSAGES D'ERREUR_____________________________________________________
  548.  
  549. SI (ierr > 0) ;
  550.  
  551.    SI (EGA ierr 10) ;
  552.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEPL icomp)
  553.               ' Noeud' ino 'Etat' ietat ;
  554.       ERRE 5 ;
  555.    FINSI ;
  556.  
  557.    SI (EGA ierr 11) ;
  558.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEPL icomp)
  559.               ' Noeud' ino 'Etat' ietat ;
  560.       ERRE 5 ;
  561.    FINSI ;
  562.  
  563.    SI (EGA ierr 20) ;
  564.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LDEFO icomp)
  565.               ' Point Gauss' igau 'Etat' ietat ;
  566.       ERRE 5 ;
  567.    FINSI ;
  568.  
  569.    SI (EGA ierr 21) ;
  570.       Message 'Erreur relative' errrel 'Composante ' (EXTR LDEFO icomp)
  571.               ' Point Gauss' igau 'Etat' ietat ;
  572.       ERRE 5 ;
  573.    FINSI ;
  574.  
  575.    SI (EGA ierr 30) ;
  576.       Message 'Erreur absolue' errabs 'Composante ' (EXTR LCONT icomp)
  577.               ' Point Gauss' igau 'Etat' ietat ;
  578.       ERRE 5 ;
  579.    FINSI ;
  580.  
  581.    SI (EGA ierr 31) ;
  582.       Message 'Erreur relative' errrel 'Composante ' (EXTR LCONT icomp)
  583.               ' Point Gauss' igau 'Etat' ietat ;
  584.       ERRE 5 ;
  585.    FINSI ;
  586.  
  587. SINON ;
  588.  
  589.    Message 'Resultats concordants avec la solution theorique' ;
  590.    Message 'Aucune erreur relative superieure a ' critere
  591.            ' pour les composantes' ;
  592.    Message (EXTR LDEPL idepl) '  ' (EXTR LDEFO idefo) '  '
  593.            (EXTR LCONT icon1) '  ' (EXTR LCONT icon2) '  '
  594.            (EXTR LCONT icon3) ;
  595.    ERRE 0 ;
  596.  
  597. FINSI ;
  598.  
  599. FIN ;
  600.  
  601.  
  602.  
  603.  
  604.  
  605.  
  606.  
  607.  
  608.  
  609.  
  610.  
  611.  

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