Télécharger la plaquette Cast3M


fr en

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