Télécharger elas7.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : elas7.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************************
  6. * Test elas7.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 1 ;
  23. FINSI ;
  24.  
  25. SAUT PAGE
  26. *______________________________________________________
  27. * TEST ELAS7 *
  28. * *
  29. * COQUE AVEC BRANCHEMENT *
  30. * *
  31. * TEST ELEMENTAIRE NAFEMS : ANALYSE LINEAIRE *
  32. * ELASTIQUE *
  33. * numero LE9 AXISYMETRIQUE *
  34. * *
  35. * L'enceinte est constituee par un cylindre et *
  36. * par une sphere (demie), reposant sur un cylindre. *
  37. * L'enceinte est soumise a une pression interne *
  38. * uniforme de 1e6Pa. *
  39. * On cherche la contrainte axiale sur la surface *
  40. * exterieure du cylindre superieure au point C. La *
  41. * contrainte obtenue est comparee au resultat *
  42. * theorique. *
  43. * *
  44. *______________________________________________________
  45. *
  46. OPTI MODE AXIS;
  47. OPTI ELEM SEG2;
  48. *
  49. *=====================================================*
  50. * geometrie : maillage *
  51. *=====================================================*
  52. *
  53. * Dimension en metres
  54. *
  55. * Points
  56. *
  57. COS75 = COS 7.5; COS375 = COS 37.5; COS45 = COS 45;
  58. SIN75 = SIN 7.5; SIN375 = SIN 37.5; SIN45 = SIN 45;
  59. R = 2 ** -.5;
  60. *
  61. RA = R; ZA = 0.; PA = RA ZA;
  62. RB = 0; ZB = COS45; PB = RB ZB;
  63. RC = R; ZC = 1.; PC = RC ZC;
  64. RD = R; ZD = 2.; PD = RD ZD;
  65. RO = 0; ZO = 1 + COS45; PO = RO ZO;
  66. PA1 = PC PLUS (0 -0.125);
  67. PD1 = PC PLUS (0 0.125);
  68. PB1 = PB PLUS (SIN375 (1 - COS375));
  69. *
  70. * Droites, cercles et lignes
  71. *
  72. * coefficient de finesse du maillage
  73. *
  74. K = 1;
  75. *
  76. SI (K EGA 1);
  77. TITR ' LE9 : MAILLAGE GROSSIER';
  78. SINON;
  79. TITR ' LE9 : MAILLAGE FIN';
  80. FINSI;
  81. *
  82. N1 = 4 * K;
  83. AA1 = PA DROI N1 PA1;
  84. A1C = PA1 DROI N1 PC;
  85. CD1 = PC DROI N1 PD1;
  86. D1D = PD1 DROI N1 PD;
  87. B1C = PB1 CERC N1 PO PC;
  88. *
  89. N2 = 5 * K;
  90. BB1 = PB CERC N2 PO PB1;
  91. *
  92. * Maillage
  93. *
  94. SUR1 = BB1 ET B1C;
  95. SUR2 = CD1 ET D1D;
  96. SUR3 = AA1 ET A1C;
  97. SURF1 = SUR1 ET SUR2 ET SUR3;
  98. *
  99. SI (NEG GRAPH 'N');
  100. TITR 'ELAS7 : MAILLAGE';
  101. TRAC QUAL SURF1;
  102. FINSI;
  103. *
  104. ELSURF1 = NBEL SURF1;
  105. NOSURF1 = NBNO SURF1;
  106. *
  107. *=====================================================*
  108. * modele - affecte - materiau - caracteristique *
  109. * rigidite - conditions aux limites *
  110. *=====================================================*
  111. *
  112. MODL1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE COQ2;
  113. MODL2 = MODE SUR2 MECANIQUE ELASTIQUE ISOTROPE COQ2;
  114. MODL3 = MODE SUR3 MECANIQUE ELASTIQUE ISOTROPE COQ2;
  115. *
  116. MATR1 = MATE MODL1 YOUN 210E9 NU 0.3 EPAI 0.01;
  117. MATR2 = MATE MODL2 YOUN 210E9 NU 0.3 EPAI 0.01;
  118. MATR3 = MATE MODL3 YOUN 210E9 NU 0.3 EPAI 0.01;
  119. *
  120. RI1 = RIGI MODL1 MATR1;
  121. RI2 = RIGI MODL2 MATR2;
  122. RI3 = RIGI MODL3 MATR3;
  123. *
  124. * Conditions aux limites
  125. * - point pa : encastrement complet.
  126. *
  127. CL1 = BLOQ ROTA DEPL PA;
  128. *
  129. RIG1 = RI1 ET RI2 ET RI3 ET CL1;
  130. *
  131. *=====================================================*
  132. * Chargement *
  133. *=====================================================*
  134. *
  135. * Pression interne uniforme de 1e6Pa.
  136. *
  137. MOP1 = 'MODE' SUR1 'CHARGEMENT' 'PRESSION' 'COQ2' ;
  138. MAP1 = 'PRES' MOP1 'PRES' -1.E6 ;
  139. CH1 = 'BSIG' MOP1 MAP1 ;
  140. MOP2 = 'MODE' SUR2 'CHARGEMENT' 'PRESSION' 'COQ2' ;
  141. MAP2 = 'PRES' MOP2 'PRES' -1.E6 ;
  142. CH2 = 'BSIG' MOP2 MAP2 ;
  143. CHA1 = CH1 ET CH2;
  144. *
  145. *=====================================================*
  146. * resolution : champs de deplacements *
  147. *=====================================================*
  148. *
  149. DEP1 = RESO RIG1 CHA1;
  150. *
  151. * Trace facultatif de la deformee et des reactions
  152. *
  153. SI (NEG GRAPH 'N');
  154. AMPVEC = .4E-6;
  155. REA1 = REAC DEP1 RIG1;
  156. VECT1 = VECT CHA1 (AMPVEC * 1) FR FZ ROUGE;
  157. VECT2 = VECT REA1 (AMPVEC * 1) FR FZ BLEU ;
  158. DEF0 = DEFO 0. DEP1 SURF1;
  159. DEF1 = DEFO 40. DEP1 SURF1 (VECT1 ET VECT2) VERT;
  160. TITR 'ELAS7 : DEFORMEE, CHARGEMENT, REACTIONS';
  161. TRAC SURF1 (DEF0 ET DEF1);
  162. FINSI;
  163. *
  164. *=====================================================*
  165. * extraction d'une contrainte *
  166. * et comparaison avec solution analytique *
  167. *=====================================================*
  168. *
  169. * On cherche la contrainte axiale sur la surface
  170. * exterieure du cylindre superieure au point C.
  171. *
  172. * Calcul du champs de contraintes generalisees
  173. * uniquement sur le cylindre superieur (sur2).
  174. *
  175. * Reduction du chpoint des deplacements au cylindre
  176. * superieur.
  177. *
  178. DEP2 = REDU DEP1 SUR2;
  179. CHAM1 = SIGM MODL2 MATR2 DEP2;
  180. *
  181. * Calcul du champs de contraintes au sens des milieux
  182. * continus. Pour cela on utilise l'operateur CALP .
  183. *
  184. CHAM2 = CALP CHAM1 MATR2 MODL2 INFE;
  185. *
  186. * Calcul des contraintes aux noeuds du maillage
  187. *
  188. CHAM3 = CHAN NOEUD MODL2 CHAM2;
  189. CHPO3 = CHAN CHPO MODL2 CHAM3;
  190. *
  191. * Extraction de la contrainte recherchee
  192. *
  193. SMZZC = EXTR CHPO3 SMXX PC;
  194. SMZZC = SMZZC / 1.E6;
  195. *
  196. SMZZC1 = -319.9;
  197. *
  198. * Calcul de l'ecart avec la solution de reference
  199. *
  200. ERGZZC = 100 * (ABS ((SMZZC - SMZZC1) / SMZZC1));
  201. *
  202. *=====================================================*
  203. * affichage des resultats *
  204. *=====================================================*
  205. *
  206. MESS ' RESULTATS : TEST ELAS7 ';
  207. MESS ' ---------------------- ';
  208. *
  209. MESS ' On cherche la contrainte axiale sur la surface
  210. exterieure';
  211. MESS ' du cylindre superieure au point C.';
  212. ' Elle est comparee a une valeur theorique obtenue
  213. analytiquement.';
  214. MESS ' Contrainte axiale theorique en C:' SMZZC1 'MPa';
  215. MESS ' Contrainte axiale calculee en C :' SMZZC 'MPa';
  216. MESS ' Soit un ecart de :' ERGZZC '%';
  217. *
  218. SI (K EGA 1);
  219. MESS ' MAILLAGE GROSSIER';
  220. SINON;
  221. MESS ' MAILLAGE FIN';
  222. FINSI;
  223. *
  224. MESS ' Nombre d elements :' ELSURF1;
  225. MESS ' Nombre de noeuds :' NOSURF1;
  226. MESS ' soit ' (NOSURF1 * 3) 'd.d.l.';
  227. *
  228. *=====================================================*
  229. * code fonctionnement *
  230. *=====================================================*
  231. *
  232. * L'ecart maximum entre valeur theorique et calculee
  233. * doit etre inferieure a 2.5%.
  234. *
  235. SI (ERGZZC <EG 2.5);
  236. ERRE 0;
  237. SINON;
  238. ERRE 5;
  239. FINSI;
  240. *
  241. *=====================================================*
  242. * Temps de calcul et fin *
  243. *=====================================================*
  244. *
  245. TEMPS;
  246. *
  247. FIN;
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  

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