Télécharger elas5.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : elas5.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************************
  6. * Test elas5.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 ELAS5 *
  28. * *
  29. * CYLINDRE - SPHERE AXISYMETRIQUE *
  30. * *
  31. * cas-test NAFEMS : test numero LE7 *
  32. * *
  33. * La structure est constituée d'un cylindre surmonté d'une *
  34. * demi-sphère. Elle est soumise a une pression interne de 1MPa. *
  35. * *
  36. * Le problème présentant une symétrie axiale (géometrie et *
  37. * chargement), on se place en mode axisymétrique. *
  38. * *
  39. * On se propose de calculer la contrainte axiale Smyy sur la *
  40. * surface extérieure au point D. La valeur trouvée est ensuite *
  41. * comparée à une valeur théorique trouvée analytiquement. *
  42. * *
  43. *___________________________________________________________________
  44. OPTI MODE AXIS;
  45. OPTI ELEM SEG2;
  46. *
  47. *==================================================================*
  48. * geometrie : maillage *
  49. *==================================================================*
  50. * Dimension en metres
  51. *
  52. * Points
  53. *
  54. RO2 = 0.; ZO2 = 1.5; PO2 = RO2 ZO2;
  55. *
  56. RA = 0.; ZA = 2.5; PA = RA ZA;
  57. RB = 0.9814; ZB = 1.6920; PB = RB ZB;
  58. RC = 1.; ZC = 1.5; PC = RC ZC;
  59. RD = 1.; ZD = 1.4034; PD = RD ZD;
  60. RE = 1.; ZE = 1.1136; PE = RE ZE;
  61. RF = 1.; ZF = 0.; PF = RF ZF;
  62. *
  63. * Droites, cercles et lignes
  64. *
  65. * coefficient de finesse du maillage
  66. K = 2;
  67. *
  68. N1 = 6 * K;
  69. AB = PA CERC N1 PO2 PB;
  70. *
  71. N2 = 4 * K;
  72. BC = PB CERC N2 PO2 PC;
  73. *
  74. N3 = 2 * K;
  75. CD = PC DROI N3 PD;
  76. *
  77. N4 = 3 * K;
  78. DE = PD DROI N4 PE;
  79. *
  80. N5 = 4 * 2;
  81. EF = PE DROI N5 PF;
  82. *
  83. * Maillage
  84. *
  85. SUR1 = AB ET BC ET CD ET DE ET EF;
  86. *
  87. ELIM SUR1 0.001;
  88. *
  89. SI (NEG GRAPH 'N');
  90. TITR 'ELAS5 : MAILLAGE';
  91. TRAC QUAL SUR1;
  92. FINSI;
  93. *
  94. ELSUR1 = NBEL SUR1;
  95. NOSUR1 = NBNO SUR1;
  96. *
  97. *==================================================================*
  98. * modele - affecte - materiau - caracteristique *
  99. * rigidite - conditions aux limites *
  100. *==================================================================*
  101. MODL1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE COQ2;
  102. MATR1 = MATE MODL1 YOUN 210E9 NU 0.3 EPAI 0.025;
  103. *
  104. MOP = 'MODE' SUR1 'CHARGEMENT' 'PRESSION' COQ2 ;
  105. MAP = 'MATE' MOP 'PRES' 1.E6 ;
  106. *
  107. MODLT = MODL1 'ET' MOP ;
  108. MATRT = MATR1 'ET' MAP ;
  109. *
  110. RIG1 = RIGI MODLT MATRT;
  111. *
  112. * Conditions aux limites
  113. * - point pa : pas de deplacement radial et pas de rotation,
  114. * - point pf : pas de deplacement vertical.
  115. *
  116. CL1 = BLOQ ROTA PA;
  117. CL2 = BLOQ UR PA;
  118. CL3 = BLOQ UZ PF;
  119. *
  120. RIG1 = RIG1 ET CL1 ET CL2 ET CL3;
  121. *
  122. *==================================================================*
  123. * Chargement *
  124. *==================================================================*
  125. *
  126. * Pression interne uniforme de 1e6Pa.
  127. *
  128. CHA1 = 'PRES' MOP MAP ;
  129. *
  130. *==================================================================*
  131. * resolution : champs de deplacements *
  132. *==================================================================*
  133. *
  134. DEP1 = RESO RIG1 CHA1;
  135. *
  136. * Trace facultatif de la deformee et des reactions
  137. *
  138. SI (NEG GRAPH 'N');
  139. AMPVEC = .5E-6;
  140. REA1 = REAC DEP1 RIG1;
  141. VECT1 = VECT CHA1 (AMPVEC * 1) FR FZ ROUGE;
  142. VECT2 = VECT REA1 AMPVEC FR FZ BLEU ;
  143. DEF0 = DEFO 0. DEP1 SUR1;
  144. DEF1 = DEFO 100. DEP1 SUR1 (VECT1 ET VECT2) VERT;
  145. TITR 'ELAS5 : DEFORMEE, CHARGEMENT, REACTIONS';
  146. TRAC SUR1 (DEF0 ET DEF1);
  147. FINSI;
  148. *
  149. *==================================================================*
  150. * extraction d'une contrainte *
  151. *==================================================================*
  152. *
  153. * On cherche la contrainte axiale sigyy sur la surface exterieure
  154. * au point D.
  155. * Cette valeur est ensuite comparee a la valeur theorique obtenue
  156. * analytiquement.
  157. *
  158. * Calcul du champs de contraintes generalisees.
  159. *
  160. CHAM1 = SIGM MODL1 MATR1 DEP1;
  161. *
  162. * Calcul du champs de contraintes au sens des milieux continus.
  163. * Pour cela on utilise l'operateur CALP .
  164. *
  165. CHAM2 = CALP CHAM1 MATRT MODLT SUPE;
  166. *
  167. * Interpolation des contraintes au noeuds
  168. *
  169. CHAM3 = CHAN NOEUD MODL1 CHAM2;
  170. CHPO3 = CHAN CHPO CHAM3 MODL1;
  171. *
  172. * Extraction de la contrainte smyy au point D : Pour cela on extrait
  173. * la contrainte smxx du repere local des elements qui correspond a la
  174. * contrainte smyy dans le repere global.
  175. *
  176. SMYYD = EXTR CHPO3 SMXX PD;
  177. SMYYD = SMYYD / 1.E6;
  178. *
  179. SMYYD1 = 25.86;
  180. *
  181. * Calcul de l'ecart avec la valeur de reference
  182. *
  183. ERGYYD = 100 * (ABS ((SMYYD1 - SMYYD) / SMYYD1));
  184. *
  185. *==================================================================*
  186. * affichage des resultats *
  187. *==================================================================*
  188. *
  189. MESS ' Resultats ELAS5';
  190. MESS ' ---------------';
  191. *
  192. MESS ' La valeur cherchee est la contrainte axiale au point D, ';
  193. MESS ' sur la surface exterieure. Elle est comparee a la valeur ';
  194. MESS ' theorique obtenue analytiquement.';
  195. MESS ' Contrainte axiale theorique en D :' SMYYD1 'MPa';
  196. MESS ' Contrainte axiale calculee en D :' SMYYD 'MPa';
  197. MESS ' Soit un ecart de :' ERGYYD '%';
  198. *
  199. SI (K EGA 1);
  200. MESS ' Maillage grossier';
  201. SINON;
  202. MESS ' Maillage fin';
  203. FINSI;
  204. *
  205. MESS ' Nombre d elements :' ELSUR1;
  206. MESS ' Nombre de noeuds :' NOSUR1;
  207. MESS ' soit ' (NOSUR1 * 3) 'd.d.l.';
  208. *
  209. *==================================================================*
  210. * code fonctionnement *
  211. *==================================================================*
  212. *
  213. * L'ecart maximum entre valeur theorique et calculee doit etre
  214. * inferieure a 1%.
  215. *
  216. SAUT 2 LIGNE;
  217. SI (ERGYYD <EG 1.);
  218. ERRE 0;
  219. SINON;
  220. ERRE 5;
  221. FINSI;
  222. *
  223. *==================================================================*
  224. * Temps de calcul et fin *
  225. *==================================================================*
  226. *
  227. TEMPS;
  228. *
  229. FIN;
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  

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