Télécharger elas8.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : elas8.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************************
  6. * Test elas8.dgibi: Jeux de données *
  7. * --------------------------------- *
  8. *******************************************************
  9. * CAS TEST DU 91/06/13 PROVENANCE : TEST
  10.  
  11. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  12. * SINON SI GRAPH DIFFERENT DE N TOUS
  13. * LES GRAPHIQUES SONT AFFICHES
  14.  
  15. GRAPH = 'N' ;
  16.  
  17. SAUT PAGE;
  18. SI (NEG GRAPH 'N') ;
  19. OPTI ECHO 1 ;
  20. OPTI TRAC PSC ;
  21. SINO ;
  22. OPTI ECHO 0 ;
  23. FINSI ;
  24.  
  25. SAUT PAGE;
  26. *______________________________________________________
  27. * TEST ELAS8 *
  28. * *
  29. * MEMBRANE ELLIPTIQUE (Contraintes planes) *
  30. * *
  31. * cas-test NAFEMS : test numero LE1 *
  32. * *
  33. * *
  34. * Une membrane elliptique obtenue par projection *
  35. * d'arcs de cercles sur un plan, est soumis a une *
  36. * force lineique sur l'un de ces cotes. *
  37. * *
  38. * Vue les symetries le probleme est traite en *
  39. * contraintes planes. *
  40. * *
  41. * On se propose de calculer la contrainte smyy *
  42. * au point D. *
  43. * *
  44. *______________________________________________________
  45. *
  46. OPTI ELEM QUA8;
  47. *
  48. *=====================================================*
  49. * geometrie : maillage *
  50. *=====================================================*
  51. *
  52. * Dimension en metres
  53. *
  54. * Points
  55. *
  56. PB = 0. 2.75 0.;
  57. PB1 = 0. 2.75 (((3.25 ** 2) - (2.75 ** 2)) ** 0.5);
  58. PC = 3.25 0. 0.;
  59. *
  60. C1 = PC CERC 6 (0. 0. 0.) PB1;
  61. C2 = C1 PROJ CYLI (0. 0. 1) PLAN (0 0 0) (1 0 0)
  62. (0 1 0);
  63. *
  64. PA = 0. 1. 0.;
  65. PA1 = 0. 1. (((2. ** 2) - (1. ** 2)) ** 0.5);
  66. PD = 2. 0. 0.;
  67. *
  68. C3 = PD CERC 6 (0. 0. 0.) PA1;
  69. C4 = C3 PROJ CYLI (0. 0. 1) PLAN (0 0 0) (1 0 0)
  70. (0 1 0);
  71. *
  72. D1 = PA DROI 2 (0. 1.583 0.) DROI 2 PB;
  73. D3 = PC DROI 2 (2.417 0. 0.) DROI 2 PD;
  74. *
  75. ELIM (D1 ET C2 ET D3 ET C4) 0.0001;
  76. SUR1 = DALL D1 C2 D3 C4 PLAN;
  77. *
  78. SI (NEG GRAPH 'N');
  79. TITR 'ELAS8 : MAILLAGE';
  80. TRAC (0 0 1000) QUAL SUR1;
  81. FINSI;
  82. *
  83. OPTI MODE PLAN CONT;
  84. *
  85. *=====================================================*
  86. * modele - affecte - materiau - caracteristique *
  87. * rigidite - conditions aux limites *
  88. *=====================================================*
  89. *
  90. MODL1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE QUA8;
  91. MATR1 = MATE MODL1 YOUN 210E9 NU 0.3;
  92. RI1 = RIGI MODL1 MATR1;
  93. *
  94. * Conditions aux limites
  95. * - droite d1: symetrie par rapport a l'axe des ordo.
  96. * - droite d2: symetrie par rapport a l'axe des absc.
  97. *
  98. CL1 = SYMT DEPL ROTA (0 0) (0 1) SUR1 0.0001;
  99. CL2 = SYMT DEPL ROTA (0 0) (1 0) SUR1 0.0001;
  100. *
  101. RIG1 = RI1 ET CL1 ET CL2;
  102. *
  103. *=====================================================*
  104. * Chargement *
  105. *=====================================================*
  106. *
  107. * Pression externe uniforme de 10e6Pa.
  108. *
  109. CHA1 = PRES MASS MODL1 -10.E6 C2;
  110. *
  111. *=====================================================*
  112. * resolution : champs de deplacements *
  113. *=====================================================*
  114. *
  115. DEP1 = RESO RIG1 CHA1;
  116. *
  117. * Trace facultatif de la deformee et des reactions
  118. *
  119. SI (NEG GRAPH 'N');
  120. AMPVEC = 1.E-7;
  121. REA1 = REAC DEP1 RIG1;
  122. VECT1 = VECT CHA1 AMPVEC FX FY ROUGE;
  123. VECT2 = VECT REA1 AMPVEC FX FY BLEU;
  124. DEF0 = DEFO 0. DEP1 SUR1;
  125. DEF1 = DEFO 150. DEP1 SUR1 (VECT1 ET VECT2) VERT;
  126. TITR 'ELAS8 : DEFORMEE, CHARGEMENT, REACTIONS';
  127. TRAC SUR1 (DEF0 ET DEF1);
  128. FINSI;
  129. *
  130. *=====================================================*
  131. * champs de contraintes *
  132. *=====================================================*
  133. *
  134. CHAM1 = SIGM MODL1 MATR1 DEP1;
  135. *
  136. * Interpolation aux noeuds du maillage
  137. *
  138. CHAM2 = CHAN NOEUD CHAM1 MODL1;
  139. CHP2 = CHAN CHPO CHAM2 MODL1;
  140. *
  141. * Trace facultatif des contraintes aux noeuds.
  142. *
  143. SI (NEG GRAPH 'N');
  144. TITR 'ELAS8 : Contraintes Nodales';
  145. TRAC CACH SUR1 ((EXCO CHP2 SMYY)*1e-6);
  146. FINSI;
  147. *
  148. *=====================================================*
  149. * extraction et affichage d'une contrainte *
  150. *=====================================================*
  151. *
  152. * On cherche la contrainte smyy a point D.
  153. *
  154. SMYYD = EXTR CHP2 SMYY PD;SMYYD = SMYYD / 1.E6;
  155. *
  156. SMYYD1 = 92.7;
  157. ERG = 100 * (ABS ((SMYYD - SMYYD1) / SMYYD1));
  158. *
  159. MESS ' RESULTATS ';
  160. MESS ' --------- ';
  161. MESS ' CONTRAINTE SMYY THEORIQUE AU POINT D :' SMYYD1
  162. 'MPA';
  163. MESS ' CONTRAINTE SMYY CALCULEE AU POINT D :' SMYYD
  164. 'MPA';
  165. MESS ' SOIT UN ECART DE :' ERG '%';
  166. *
  167. NBNOSUR1 = NBNO SUR1;
  168. NBELSUR1 = NBEL SUR1;
  169. MESS ' NOMBRE D ELEMENTS :' NBELSUR1 'QUA8';
  170. MESS ' NOMBRE DE NOEUDS ;' NBNOSUR1 ;
  171. *
  172. *=====================================================*
  173. * code fonctionnement *
  174. *=====================================================*
  175. *
  176. SI (ERG <EG 1.);
  177. ERRE 0;
  178. SINON;
  179. ERRE 5;
  180. FINSI;
  181. *
  182. TEMPS;
  183. *
  184. FIN;
  185.  
  186.  
  187.  
  188.  
  189.  

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