Télécharger elas11.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : elas11.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *******************************************************
  6. * Test elas11.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 ELAS11 *
  28. * *
  29. * SOLIDE (CYLINDRE - CONE - SPHERE) *
  30. * SOUMIS A UN GRADIENT *
  31. * DE TEMPERATURE (Calcul axisymetrique 2D). *
  32. * *
  33. * Test NAFEMS : numero LE11 *
  34. * Nom : solid cylinder / taper / sphere temperature *
  35. * *
  36. * La structure est une enceinte epaisse. Elle est *
  37. * constituee d'une partie spherique surmontee d'une *
  38. * partie conique, elle -meme surmontee d'une partie *
  39. * cylindrique. *
  40. * *
  41. * Elle est soumis a un gradient lineaire de *
  42. * temperature s'exercant radialement et axialement. *
  43. * Le champs de temperature s'exprime comme suit : *
  44. * *
  45. * T(x,y,z) = (x**2 + y**2)**.5 + z *
  46. * ou *
  47. * T(r,z) = r + z *
  48. * *
  49. * L'ensemble est maille avec des elements massifs *
  50. * CU20. Le maillage est obtenue par rotation de 90 *
  51. * degres d'une coupe verticale. *
  52. * *
  53. * On se propose de calculer la contrainte smyy au *
  54. * point A. Cette valeur est ensuite comparee avec *
  55. * la valeur de reference obtenue a l'aide d'un *
  56. * calcul axisymetrique sur un maillage fin. *
  57. * *
  58. *______________________________________________________
  59. *
  60. OPTI MODE AXIS;
  61. OPTI ELEM QUA8;
  62. *
  63. *=====================================================*
  64. * geometrie : maillage *
  65. *=====================================================*
  66. *
  67. * Dimension en metres
  68. *
  69. * Points
  70. *
  71. PO = 0. 0.;
  72. P1 = 0. 1.;
  73. PA = 1. 0.;
  74. PB = 1.4 0.;
  75. PE = (1. * (COS 30)) (1. * (SIN 30));
  76. PD = (1.4 * (COS 30)) 0.7;
  77. PC = (COS 30) (SIN 30);
  78. PE = (COS 45) (SIN 45);
  79. PH = 0.7071 1.79;
  80. PI = 1. 1.79;
  81. PG = 1. 1.39;
  82. *
  83. * Droites, cercles et lignes
  84. *
  85. * coefficient de finesse du maillage
  86. *
  87. K = 2;
  88. *
  89. N1 = 1 * K;
  90. AB = PA DROI N1 PB;
  91. GI = PG DROI N1 PI;
  92. EC = PE CERC N1 PO PC;
  93. IH = PI DROI N1 PH;
  94. *
  95. N2 = 2 * K;
  96. BD = PB CERC N2 PO PD;
  97. DG = PD DROI N2 PG;
  98. CA = PC CERC N2 PO PA;
  99. HE = PH DROI N2 PE;
  100. *
  101. * maillage par elements massifs du volume
  102. *
  103. SUR1 = DALL AB (BD ET DG ET GI) IH (HE ET EC ET CA) PLAN;
  104. *
  105. SI (NEG GRAPH 'N');
  106. SI (K EGA 1);
  107. TITR 'ELAS11 : MAILLAGE GROSSIER';
  108. SINON;
  109. TITR 'ELAS11 : MAILLAGE FIN';
  110. FINSI;
  111. TRAC QUAL SUR1;
  112. FINSI;
  113. *
  114. ELSUR1 = NBEL SUR1;
  115. NOSUR1 = NBNO SUR1;
  116. *
  117. *=====================================================*
  118. * modele - materiau *
  119. * caracteristique - rigidite *
  120. *=====================================================*
  121. *
  122. MODL1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE QUA8;
  123. MATR1 = MATE MODL1 YOUN 210E9 NU 0.3 ALPH 2.3E-4;
  124. RI1 = RIGI MODL1 MATR1;
  125. *
  126. *=====================================================*
  127. * Conditions aux limites *
  128. *=====================================================*
  129. *
  130. * - deplacement vertical nul pour les faces superieure
  131. * et inferieure.
  132. *
  133. CL1 = BLOQ UZ AB;
  134. CL2 = BLOQ UZ IH;
  135. *
  136. CLI1 = CL1 ET CL2;
  137. RIG1 = RI1 ET CLI1;
  138. *
  139. *=====================================================*
  140. * Chargement thermique *
  141. *=====================================================*
  142. *
  143. * Le solide est soumis a un gradient de temperature
  144. * lineaire, s'appliquant radialement et axialement.
  145. * La temperature est la suivante :
  146. *
  147. * T(r,z) = r + z
  148. *
  149. * Creation du champs de temperature
  150. *
  151. CHPR = COOR 1 SUR1;
  152. CHPZ = COOR 2 SUR1;
  153. *
  154. CHP4 = CHPR + CHPZ;
  155. *
  156. * Transformation du sous-type "scal" en sous-type "t"
  157. * pour chp4. Utilisation de l'operateur "vari".
  158. *
  159. LIS10 = PROG 0. 200.;
  160. LIS20 = PROG 0. 200.;
  161. EVOL1 = EVOL MANU SCAL LIS10 T LIS20;
  162. CHP5 = VARI CHP4 EVOL1 T;
  163. *
  164. * Creation du champ de contraintes equivalentes
  165. *
  166. CHAM10 = THET MODL1 MATR1 CHP5;
  167. *
  168. * Creation du champ de forces equivalentes
  169. *
  170. CHAR1 = BSIG MODL1 CHAM10;
  171. *
  172. *=====================================================*
  173. * Resolution : champs de deplacements *
  174. *=====================================================*
  175. *
  176. DEP1 = RESO RIG1 CHAR1;
  177. *
  178. * Trace facultatif de la deformee et des reactions
  179. *
  180. SI (NEG GRAPH 'N');
  181. AMPVEC = 1.E-8;
  182. REA1 = REAC DEP1 RIG1;
  183. VECT1 = VECT REA1 AMPVEC FR FZ BLEU;
  184. DEF0 = DEFO 0. DEP1 SUR1;
  185. DEF1 = DEFO 100. DEP1 SUR1 VECT1 VERT;
  186. TITR 'ELAS11 ; DEFORMEE, REACTIONS';
  187. TRAC SUR1 (DEF0 ET DEF1);
  188. FINSI;
  189. *
  190. *=====================================================*
  191. * extraction d'une contrainte *
  192. * et comparaison avec solution analytique *
  193. *=====================================================*
  194. *
  195. * On cherche la contrainte axiale sigzz au noeud A
  196. *
  197. * Calcul du champs de contraintes totales.
  198. *
  199. CHAM1 = SIGM MODL1 MATR1 DEP1;
  200. *
  201. * On retranche les contraintes d'origine thermique
  202. *
  203. CHAM1 = CHAM1 - CHAM10;
  204. *
  205. * Interpolation aux noeuds du maillage.
  206. *
  207. CHAM2 = CHAN NOEUD MODL1 CHAM1;
  208. CHPO10 = CHAN CHPO MODL1 CHAM2;
  209. *
  210. * Trace facultatif des contraintes aux noeuds.
  211. *
  212. SI (NEG GRAPH 'N');
  213. TITR 'ELAS9 : Contraintes Nodales';
  214. TRAC CACH SUR1 ((EXCO CHPO10 SMZZ)*1e-6);
  215. FINSI;
  216. *
  217. * Extraction de la contrainte recherchee.
  218. *
  219. SMZZA = EXTR CHPO10 SMZZ PA;
  220. SMZZA = SMZZA / 1.E6;
  221. *
  222. SMZZA1 = -105;
  223. *
  224. * Calcul de l'ecart avec la solution de reference.
  225. *
  226. ERGZZA = 100 * (ABS ((SMZZA - SMZZA1) / SMZZA1));
  227. *
  228. *=====================================================*
  229. * affichage des resultats *
  230. *=====================================================*
  231. *
  232. MESS ' RESULTATS ';
  233. MESS ' --------- ';
  234. *
  235. MESS ' La valeur cherchee est la contrainte axiale au point A,';
  236. MESS ' elle est comparee a une valeur theorique obtenue analytiquement.';
  237. MESS ' Contrainte verticale theorique en A :' SMZZA1 'MPa';
  238. MESS ' Contrainte verticale calculee en A :' SMZZA 'MPa';
  239. MESS ' Soit un ecart de :' ERGZZA '%';
  240. *
  241. SI (K EGA 1);
  242. MESS ' Maillage grossier';
  243. SINON;
  244. MESS ' Maillage fin';
  245. FINSI;
  246. *
  247. MESS ' Nombre d elements :' ELSUR1;
  248. MESS ' Nombre de noeuds :' NOSUR1;
  249. MESS ' Soit :' (NOSUR1 * 3) 'd.d.l.';
  250. *
  251. *=====================================================*
  252. * code fonctionnement *
  253. *=====================================================*
  254. * L'ecart maximum entre valeur theorique et
  255. * calculee doit etre inferieure a 2%.
  256. *
  257. SAUT 1 LIGNE;
  258. SI (ERGZZA <EG 2);
  259. ERRE 0;
  260. SINON;
  261. ERRE 5;
  262. FINSI;
  263. *
  264. *=====================================================*
  265. * Temps de calcul et fin *
  266. *=====================================================*
  267. TEMPS;
  268. FIN;
  269.  
  270.  
  271.  
  272.  
  273.  
  274.  
  275.  

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