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)
  104. PLAN;
  105. *
  106. SI (NEG GRAPH 'N');
  107. SI (K EGA 1);
  108. TITR 'ELAS11 : MAILLAGE GROSSIER';
  109. SINON;
  110. TITR 'ELAS11 : MAILLAGE FIN';
  111. FINSI;
  112. TRAC QUAL SUR1;
  113. FINSI;
  114. *
  115. ELSUR1 = NBEL SUR1;
  116. NOSUR1 = NBNO SUR1;
  117. *
  118. *=====================================================*
  119. * modele - materiau *
  120. * caracteristique - rigidite *
  121. *=====================================================*
  122. *
  123. MODL1 = MODE SUR1 MECANIQUE ELASTIQUE ISOTROPE QUA8;
  124. MATR1 = MATE MODL1 YOUN 210E9 NU 0.3 ALPH 2.3E-4;
  125. RI1 = RIGI MODL1 MATR1;
  126. *
  127. *=====================================================*
  128. * Conditions aux limites *
  129. *=====================================================*
  130. *
  131. * - deplacement vertical nul pour les faces superieure
  132. * et inferieure.
  133. *
  134. CL1 = BLOQ UZ AB;
  135. CL2 = BLOQ UZ IH;
  136. *
  137. CLI1 = CL1 ET CL2;
  138. RIG1 = RI1 ET CLI1;
  139. *
  140. *=====================================================*
  141. * Chargement thermique *
  142. *=====================================================*
  143. *
  144. * Le solide est soumis a un gradient de temperature
  145. * lineaire, s'appliquant radialement et axialement.
  146. * La temperature est la suivante :
  147. *
  148. * T(r,z) = r + z
  149. *
  150. * Creation du champs de temperature
  151. *
  152. CHPR = COOR 1 SUR1;
  153. CHPZ = COOR 2 SUR1;
  154. *
  155. CHP4 = CHPR + CHPZ;
  156. *
  157. * Transformation du sous-type "scal" en sous-type "t"
  158. * pour chp4. Utilisation de l'operateur "vari".
  159. *
  160. LIS10 = PROG 0. 200.;
  161. LIS20 = PROG 0. 200.;
  162. EVOL1 = EVOL MANU SCAL LIS10 T LIS20;
  163. CHP5 = VARI CHP4 EVOL1 T;
  164. *
  165. * Creation du champ de contraintes equivalentes
  166. *
  167. CHAM10 = THET MODL1 MATR1 CHP5;
  168. *
  169. * Creation du champ de forces equivalentes
  170. *
  171. CHAR1 = BSIG MODL1 CHAM10;
  172. *
  173. *=====================================================*
  174. * Resolution : champs de deplacements *
  175. *=====================================================*
  176. *
  177. DEP1 = RESO RIG1 CHAR1;
  178. *
  179. * Trace facultatif de la deformee et des reactions
  180. *
  181. SI (NEG GRAPH 'N');
  182. AMPVEC = 1.E-8;
  183. REA1 = REAC DEP1 RIG1;
  184. VECT1 = VECT REA1 AMPVEC FR FZ BLEU;
  185. DEF0 = DEFO 0. DEP1 SUR1;
  186. DEF1 = DEFO 100. DEP1 SUR1 VECT1 VERT;
  187. TITR 'ELAS11 ; DEFORMEE, REACTIONS';
  188. TRAC SUR1 (DEF0 ET DEF1);
  189. FINSI;
  190. *
  191. *=====================================================*
  192. * extraction d'une contrainte *
  193. * et comparaison avec solution analytique *
  194. *=====================================================*
  195. *
  196. * On cherche la contrainte axiale sigzz au noeud A
  197. *
  198. * Calcul du champs de contraintes totales.
  199. *
  200. CHAM1 = SIGM MODL1 MATR1 DEP1;
  201. *
  202. * On retranche les contraintes d'origine thermique
  203. *
  204. CHAM1 = CHAM1 - CHAM10;
  205. *
  206. * Interpolation aux noeuds du maillage.
  207. *
  208. CHAM2 = CHAN NOEUD MODL1 CHAM1;
  209. CHPO10 = CHAN CHPO MODL1 CHAM2;
  210. *
  211. * Trace facultatif des contraintes aux noeuds.
  212. *
  213. SI (NEG GRAPH 'N');
  214. TITR 'ELAS9 : Contraintes Nodales';
  215. TRAC CACH SUR1 ((EXCO CHPO10 SMZZ)*1e-6);
  216. FINSI;
  217. *
  218. * Extraction de la contrainte recherchee.
  219. *
  220. SMZZA = EXTR CHPO10 SMZZ PA;
  221. SMZZA = SMZZA / 1.E6;
  222. *
  223. SMZZA1 = -105;
  224. *
  225. * Calcul de l'ecart avec la solution de reference.
  226. *
  227. ERGZZA = 100 * (ABS ((SMZZA - SMZZA1) / SMZZA1));
  228. *
  229. *=====================================================*
  230. * affichage des resultats *
  231. *=====================================================*
  232. *
  233. MESS ' RESULTATS ';
  234. MESS ' --------- ';
  235. *
  236. MESS ' La valeur cherchee est la contrainte axiale
  237. au point A,';
  238. ' elle est comparee a une valeur theorique obtenue
  239. analytiquement.';
  240. MESS ' Contrainte verticale theorique en A :' SMZZA1
  241. 'MPa';
  242. MESS ' Contrainte verticale calculee en A :' SMZZA
  243. 'MPa';
  244. MESS ' Soit un ecart de :' ERGZZA '%';
  245. *
  246. SI (K EGA 1);
  247. MESS ' Maillage grossier';
  248. SINON;
  249. MESS ' Maillage fin';
  250. FINSI;
  251. *
  252. MESS ' Nombre d elements :' ELSUR1;
  253. MESS ' Nombre de noeuds :' NOSUR1;
  254. MESS ' Soit :' (NOSUR1 * 3) 'd.d.l.';
  255. *
  256. *=====================================================*
  257. * code fonctionnement *
  258. *=====================================================*
  259. * L'ecart maximum entre valeur theorique et
  260. * calculee doit etre inferieure a 2%.
  261. *
  262. SAUT 1 LIGNE;
  263. SI (ERGZZA <EG 2);
  264. ERRE 0;
  265. SINON;
  266. ERRE 5;
  267. FINSI;
  268. *
  269. *=====================================================*
  270. * Temps de calcul et fin *
  271. *=====================================================*
  272. TEMPS;
  273. FIN;
  274.  
  275.  
  276.  
  277.  
  278.  

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