Télécharger la plaquette Cast3M


fr en

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

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