Télécharger sort_MAILLAGE.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : sort_MAILLAGE.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ***************************************************************
  6. * Test sort_MAILLAGE.dgibi: Jeux de données *
  7. * ----------------------------------------- *
  8. ***************************************************************
  9.  
  10. *___________________________________________________________________
  11. * TEST sort_MAILLAGE.dgibi *
  12. * *
  13. * Sortie d'un MAILLAGE avec SORT
  14. * Relecture avec LIRE (Option non testee a-priori) *
  15.  
  16. * Cas-Test de Verification & Validation
  17. *__________________________________________________________________*
  18. *
  19. OPTI ELEM CU20;
  20. GRAPH = MOT 'N';
  21. *
  22. *==================================================================*
  23. * geometrie : maillage *
  24. *==================================================================*
  25. *
  26.  
  27. K = 2;
  28.  
  29. PB = 0. 2.75 0.;
  30. PB1 = 0. 2.75 (((3.25 ** 2) - (2.75 ** 2)) ** 0.5);
  31. PC = 3.25 0. 0.;
  32. *
  33. C1 = PC CERC (6 * K) (0. 0. 0.) PB1;
  34. C2 = C1 PROJ DIRE (0. 0. 1) PLAN (0 0 0) (1 0 0) (0 1 0);
  35. *
  36. PA = 0. 1. 0.;
  37. PA1 = 0. 1. (((2. ** 2) - (1. ** 2)) ** 0.5);
  38. PD = 2. 0. 0.;
  39. *
  40. C3 = Pd CERC (6 * K) (0. 0. 0.) PA1;
  41. C4 = C3 PROJ DIRE (0. 0. 1) PLAN (0 0 0) (1 0 0) (0 1 0);
  42. *
  43. D1 = PA DROI (2 * K) (0. 1.583 0.) DROI (2 * K) PB;
  44. D3 = PC DROI (2 * K) (2.417 0. 0.) DROI (2 * K) PD;
  45. *
  46. ELIM (D1 ET C2 ET D3 ET C4) 0.0001;
  47. SUR1 = DALL D1 C2 D3 C4 PLAN;
  48. *
  49. VOL1 = SUR1 VOLU (2 * K) TRAN (0. 0. -0.6);
  50. *
  51. * Surfaces et arc d'ellipse pour conditions aux limites
  52. *
  53. FACE1 = D1 TRAN (2 * K) (0. 0. -0.6);
  54. FACE2 = C2 TRAN (2 * K) (0. 0. -0.6);
  55. FACE3 = D3 TRAN (2 * K) (0. 0. -0.6);
  56. C2CL = C2 PLUS (0. 0. -0.3);
  57. *
  58. ELIM (VOL1 ET FACE1 ET FACE2 ET FACE3 ET C2CL) 0.0001;
  59. *
  60. SI (NEG GRAPH 'N');
  61. TITR 'ELAS9 : MAILLAGE';
  62. TRAC (-1000 -1000 1000) FACE QUAL (COUL VOL1 BLAN);
  63. FINSI;
  64.  
  65. *
  66. *==================================================================*
  67. * sortie : maillage *
  68. *==================================================================*
  69. *
  70. OPTI SORT 'Sort_MAILLAGE.mesh';
  71. SORT VOL1;
  72. OPTI SORT 'BIDON';
  73. VOL2 = VOL1;
  74. OUBL VOL1;
  75. NBNO2 = NBNO VOL2;
  76. NBEL2 = NBEL VOL2;
  77. *
  78. *==================================================================*
  79. * Relecture : maillage *
  80. *==================================================================*
  81. *
  82. OPTI LECT 'Sort_MAILLAGE.mesh';
  83. NBNO1 = NBNO VOL1;
  84. NBEL1 = NBEL VOL1;
  85.  
  86. *
  87. *==================================================================*
  88. * Verification : Comparaison des MAILLAGES *
  89. *==================================================================*
  90. *
  91. ELIM VOL1 VOL2 1.E-5;
  92. UNIQ (VOL1 ET VOL2);
  93.  
  94. INT1 = INTE VOL1 VOL2;
  95. NBEI = NBEL INT1;
  96.  
  97. OPTI ECHO 0;
  98. MESS 'NOEUDS ' NBNO2 NBNO1;
  99. MESS 'ELEMENTS' NBEL2 NBEL1;
  100. MESS 'INTERSEC' NBEI ;
  101.  
  102. SI ((NEG NBNO2 NBNO1) OU (NEG NBEL2 NBEL1) OU (NEG NBEL1 NBEI));
  103. ERRE 5;
  104. FINS;
  105. OPTI ECHO 1;
  106. * OPTI DONN 5;
  107.  
  108. FIN;
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  

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