Télécharger la plaquette Cast3M


fr en

Télécharger testfer.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  testfer.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************************
  6. *         Test testfer.dgibi: jeux de données         *
  7. *         -----------------------------------         *
  8. *                                                     *
  9. * Toit cylindrique soumis à son propre poids          *
  10. *******************************************************
  11.  
  12. *OPTI ECHO 0;
  13.  
  14. *** affichage
  15. GRAPH = 'N' ;
  16.  
  17. *** titre du maillage
  18. NOM1 = 'TOIT';
  19.  
  20. *** données du problème ------------------------------------------------
  21. *** caractéristiques du matériau
  22. YOUN1 = 3e+10;
  23. NU1   = 0.;
  24. RHO1  = 2.0833e+4; 
  25.  
  26. *** dimensions
  27. LON1  = 3.;
  28. EPAI1 = LON1 / 100.;
  29.  
  30. *** options de calcul --------------------------------------------------
  32.  
  33. *** maillage -----------------------------------------------------------
  34. OPTI ELEM SEG3;
  35.  
  36. *** coordonnées des points
  37. P1   =  0.                  LON1               0.;
  38. P2   =  (LON1 * (sin 40.)) (LON1 * (cos 40.))  0.;
  39. CEN1 =  0.                  0.                 0.;
  40.  
  41. *** nb éléments
  42. N1 = 20;
  43.  
  44. *** lignes
  45. L1 = CERC N1 P1 CEN1 P2;
  46.  
  47. *** création de la surface
  48. OPTI ELEM TRI6;
  49.  
  50. VEC1 = 0. 0. LON1;
  51. S200 = L1 TRAN  N1  VEC1;
  52.  
  53. *** affichage
  54. SI(NEG GRAPH 'N');
  55.    TRAC FACE S200; 
  56.  
  57. *'LISTE' (noeud 84) ;
  58. *'OPTION' donn 5 ;
  59.    FINSI;
  60.  
  61. *** modèle mécanique ---------------------------------------------------
  62. MOD1 = MODE S200 MECANIQUE ELASTIQUE ISOTROPE COQ6;
  63.  
  64. *** matériau et caractéristiques ---------------------------------------
  65. MAT1 = MATER MOD1 YOUNG YOUN1 NU NU1 RHO RHO1;
  66. CAR1 = CARA MOD1 EPAI EPAI1;
  67. MAT1 = MAT1 ET CAR1;
  68.  
  69. *** construction de la matrice de rigidité -----------------------------
  70. RIG1 = RIGI MOD1 MAT1;
  71.  
  72. *** définition des conditions aux limites ------------------------------
  73. *** identification des côtés pour application des conditions 
  74. *** aux limites 
  75. LC10 LC20 LC30 LC40 = COTE S200;
  76.  
  77. SI(NEG GRAPH 'N');
  78.    TRAC QUALI (LC10 ET LC20 ET LC30 ET LC40);
  79. FINSI;
  80.  
  81. BLOQ30 = BLOQ UX UY RZ LC30;
  82. BLOQ40 = BLOQ UX RY RZ LC40;
  83. BLOQ10 = BLOQ UZ RX RY LC10;
  84.  
  85. RIG1   = RIG1 ET BLOQ10 ET BLOQ40 ET BLOQ30;
  86.  
  87. *** poids propre -------------------------------------------------------
  88. MAS1 = MASS MOD1 MAT1;
  89. CHA1 = MANU CHPO S200 1 UY (-1.*10.);
  90. CHA1 = MAS1 * CHA1;
  91.  
  92. SI(NEG GRAPH 'N');
  93.    VEC1 = VECT CHA1 .01 FX FY FZ ROUGE;
  94.    TRAC VEC1 S200;
  95. FINSI;
  96.  
  97. *** résolution ---------------------------------------------------------
  98. TEMPS ZERO;
  99. RES1 = RESO RIG1 CHA1;
  100. TEMPS;
  101. SAUT 1 LIGNE;
  102.  
  103. *** post-traitement ----------------------------------------------------
  104.  
  105. SI(NEG GRAPH 'N');
  106. ****** maillage initial et déformée
  107.    DEF0 = DEFO S200 RES1 0.;
  108.    DEF1 = DEFO S200 RES1 VERT;
  109.    TITR 'TOIT : DEFORMEE';
  110.    TRAC S200 (DEF0 ET DEF1);
  111. FINSI;
  112.  
  113. *** déplacement aux points P1 et P2
  114. MESS 'Deplacement Uy au point P1  :' (EXTR RES1 UY P1);
  115. MESS 'Deplacement Uy au point P2  :' (EXTR RES1 UY P2);
  116.  
  117. *** CHPO pour tester SORTIR option FERVIEW -----------------------------
  118.  
  119. baba = 'NOMC' 'XX' ('COORDONNEE' 1 S200) ;
  120. bobo = 'NOMC' 'YY' ('COORDONNEE' 2 S200) ;
  121. boba = 'NOMC' 'ZZ' ('COORDONNEE' 3 S200) ;
  122. bubu = baba 'ET' bobo 'ET' boba ;
  123. byby =(1. '-' baba) 'ET' (1. '-' bobo) 'ET' (1. '-' boba) ; 
  124.  
  125. 'OPTION' SORT 'testfer_MAILLAGE.asc' ;
  126. 'SORTIR'  'FER' S200  ;
  127.  
  128. 'OPTION' SORT 'testfer_CHPOINT.asc' ;
  129. 'SORTIR' 'FER' S200  bubu  RES1 'TITR' 'Chpoints test' 'TEMP' 1. ;
  130.  
  131. 'SORTIR' 'FER'  bubu RES1 S200 'TITR' 'Chpoints test' 'TEMP' 1. ;
  132. 'SORTIR' 'FER' S200  byby RES1 'SUIT' 'TEMP' 2.  ;
  133. 'SORTIR' 'FER' S200 byby  RES1 'SUIT' 2 'TEMP' 2. ;
  134.  
  135.  
  136. FIN;
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  

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