Télécharger testkcha.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : testkcha.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************** CAS TEST : testkcha.dgibi ************************
  5. *
  6. * Cas-Test vérifiant le bon fonctionnement de l'opérateur 'KCHA' dans
  7. * les deux sens.
  8. *
  9. OPTION ECHO 1 ;
  10. CRIT1 = 1.D-14 ;
  11. 'SAUTER' 'PAGE';
  12.  
  13. 'TITRE' 'Validation du l opérateur KCHA';
  14.  
  15. 'OPTION' 'DIME' 2 'ELEM' 'QUA4' ;
  16.  
  17. *
  18. * Création maillage :
  19. *
  20. INX = 4 ;
  21. INY = 2 ;
  22. *
  23. X0 = 0. ;
  24. Y0 = 0. ;
  25. *
  26. LX = 1. ;
  27. LY = 1. ;
  28. *
  29. P1 = X0 Y0 ;
  30. P2 = P1 'PLUS' (LX 0.) ;
  31. P3 = P1 'PLUS' (LX LY) ;
  32. P4 = P1 'PLUS' (0. LY) ;
  33. *
  34. L1 = 'DROIT' P1 INX P2 ;
  35. L2 = 'DROIT' P2 INY P3 ;
  36. L3 = 'DROIT' P3 INX P4 ;
  37. L4 = 'DROIT' P4 3 P1 ;
  38. *
  39. LTOT = L1 'ET' L2 'ET' L3 'ET' L4 ;
  40. MTOT = 'SURFACE' LTOT;
  41. *
  42. MTOT1 = 'ELEM' MTOT 'TRI3' ;
  43. MTOT2 = 'ELEM' MTOT 'QUA4' ;
  44. *
  45. MTOTF = 'CHANGER' MTOT 'QUAF' ;
  46. MTOTF1 = 'CHANGER' MTOT1 'QUAF' ;
  47. MTOTF2 = 'CHANGER' MTOT2 'QUAF' ;
  48. ELIM 0.01 ( MTOTF 'ET' MTOTF1 'ET' MTOTF2) ;
  49. *TRAC MTOT;
  50. *
  51. * Construction modèles et compagnie
  52. *
  53. *HYTOT = 'DOMA' MTOT ;
  54. *HYDO1 = 'DOMA' MTOT1 'INCL' HYTOT 0.01;
  55. *HYDO2 = 'DOMA' MTOT2 'INCL' HYTOT 0.01;
  56. *
  57. MODHYB = 'MODELE' MTOTF 'DARCY' 'ISOTROPE' ;
  58. MODHY1 = 'MODELE' MTOTF1 'DARCY' 'ISOTROPE' ;
  59. MODHY2 = 'MODELE' MTOTF2 'DARCY' 'ISOTROPE' ;
  60.  
  61. 'LISTE' ('DOMA' modhyb 'TABLE');
  62. *
  63. * ==========
  64. * 1er test :
  65. * ==========
  66. *
  67. * Champ-point test, transformé en CHAMELEM puis ren CHPO :
  68. *
  69. CHP1 = 'COOR' 1 ('DOMA' MODHYB 'CENTRE' ) ;
  70. CHAM1 = 'KCHA' MODHYB 'CHAM' CHP1 ;
  71. CHP2 = 'KCHA' MODHYB 'CHPO' CHAM1 ;
  72.  
  73. * test :
  74. err1 = 'MAXIMUM' ('ABS' (CHP2 - CHP1));
  75. *
  76. * =========
  77. * 2e test :
  78. * =========
  79. *
  80. * Champ par élément test, transformé en champ-point puis ren CHAMELEM:
  81. * Dû à une erreur dans l'opérateur 'MANU', il faut -aujourd'hui-
  82. * imposer un nom de constituant à la main pour chaque chamelem créé.
  83. *
  84. CHAM31 = 'MANU' 'CHML' MODHY1 'COM1' 10. 'GRAVITE' ;
  85. CHAM31 = 'CHANGER' 'CONS' cham31 'HOP' ;
  86. CHAM32 = 'MANU' 'CHML' MODHY2 'COM2' 20. 'GRAVITE' ;
  87. CHAM32 = 'CHANGER' 'CONS' cham32 'HOP' ;
  88. CHAM3 = CHAM31 'ET' CHAM32 ;
  89. CHAM3 = 'REDU' cham3 ('DOMA' MODHYB 'MAILLAGE') ;
  90.  
  91. CHP3 = 'KCHA' MODHYB 'CHPO' CHAM3;
  92. CHAM4 = 'KCHA' MODHYB 'CHAM' CHP3;
  93. CHAM4 = 'CHANGER' 'CONS' cham4 'HOP' ;
  94.  
  95. * test :
  96. err2 = 'MAXIMUM' ('ABS' (CHAM3 - CHAM4));
  97.  
  98. *
  99. * =========
  100. * 3e test :
  101. * =========
  102. *
  103. * option 'QUAF' pour appuyer le chamelem sur le maillage quaf au lieu du
  104. * maillage de base
  105. CHAM4 = 'KCHA' MODHYB CHP3 'CHAM' 'QUAF' ;
  106. * il faut préciser que le chamelem s'appuie sur des quaf (sinon erreur)
  107. CHP4 = 'KCHA' modhyb cham4 'CHPO' 'QUAF' ;
  108.  
  109. * test :
  110. err3 = 'MAXIMUM' ('ABS' (CHP4 - CHP3));
  111.  
  112. *
  113. * =========
  114. * 4e test :
  115. * =========
  116. *
  117. * Création chpoint à partir d'un chamelem non aux centres
  118. * kcha fera alors la moyenne par élément sur tous les points
  119. xx = 'COOR' 1 MTOT ;
  120. * appuyé aux noeuds du quaf
  121. cham5 = 'CHANGER' xx modhyb 'CHAM' ;
  122. chpo5 = 'KCHA' cham5 modhyb 'CHPO' 'QUAF' ;
  123. * appuyé aux noeud du maillage
  124. cham6 = 'CHANGER' xx mtot 'CHAM' ;
  125. chpo6 = 'KCHA' cham6 modhyb 'CHPO' ;
  126.  
  127. * =======
  128. * Sortie
  129. * =======
  130.  
  131. 'SI' ((err1 > crit1) 'OU' (err2 > crit1) 'OU' (err3 > crit1)) ;
  132. 'ERREUR' 5;
  133. 'SINON';
  134. 'ERREUR' 0;
  135. 'FINSI';
  136.  
  137. *'TEMPS' 'SGAC' 'IMPR' ;
  138.  
  139. 'FIN' ;
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  

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