Télécharger la plaquette Cast3M


fr en

Télécharger thpl1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  thpl1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *   Test Thpl1.dgibi: Jeux de données     *
  5. *   ---------------------------------     *
  6. *                                         *
  7. * CAS TEST DU 92/03/20     PROVENANCE : TC1
  8.  
  9.  
  10. *******************************************
  11. *    Test thpl1.dgibi: jeux de données    *
  12. *    ---------------------------------    *
  13. *******************************************
  14.  
  15. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  16. * SINON SI GRAPH DIFFERENT DE N TOUS
  17. * LES GRAPHIQUES SONT AFFICHES
  18.  
  19. GRAPH = 'N' ;
  20.  
  21. SAUT PAGE;
  22. SI (NEG GRAPH 'N') ;
  23.   OPTI ECHO 1   ;
  24.   OPTI TRAC PSC   ;
  25. SINO ;
  26.   OPTI ECHO 0   ;
  27. FINSI ;
  28.  
  29. SAUT PAGE;
  30.  
  31. *******************************************
  32. *                   THPL1                 *
  33. *                                         *
  34. *  BARREAU CYLINDRIQUE SOUMIS A UN CHAMPS *
  35. *            DE TEMPERATURE               *
  36. *                                         *
  37. *    Test de la Commission V.P.C.S.       *
  38. *    Groupe Thermo-plasticite  Statique   *
  39. *    non-lineaire                         *
  40. *    Test numero SSNA01/90                *
  41. *                                         *
  42. *        Un barreau cylindrique de 100 mm *
  43. *    de hauteur et de 20 mm de diametre   *
  44. *    est bloque axialement a ces deux     *
  45. *    extremites. Il est soumis a un champ *
  46. *    constant de temperature, de 500°C.   *
  47. *                                         *
  48. *        On traite le probleme en mode    *
  49. *    axisymetrique avec la procedure      *
  50. *    PASAPAS car le materiau est plastique*
  51. *                                         *
  52. *        On se propose de calculer        *
  53. *     - la contrainte axiale              *
  54. *     - la deformation plastique axiale   *
  55. *******************************************
  56. OPTION DIME 2;
  57. OPTION MODE AXIS ELEM QUA4;
  58.  
  59. *=========================================*
  60. *         geometrie : maillage            *
  61. *=========================================*
  62. * Dimensions en millimetres.
  63. RA = 0.0 ;  ZA = 0. ;     PA = RA ZA ;
  64. RB = 10. ;  ZB = 0. ;     PB = RB ZB ;
  65. RC = 10. ;  ZC = 100.0 ;  PC = RC ZC ;
  66. RD = 0.0 ;  ZD = 100.0 ;  PD = RD ZD ;
  67.  
  68. DAB = DROI 1  PA PB;
  69. DBC = DROI 10 PB PC;
  70. DCD = DROI 1  PC PD;
  71. DDA = DROI 10 PD PA;
  72.  
  73. SUR = DALL DAB DBC DCD DDA PLAN;
  74.  
  75. ELIM SUR 0.001;
  76.  
  77. SI (NEG GRAPH 'N');
  78.    TITR 'THPL1 : MAILLAGE';
  79.    TRAC QUAL SUR;
  80. FINSI;
  81.  
  82. ELSUR = NBEL SUR;
  83. NOSUR = NBNO SUR;
  84.  
  85. *=========================================*
  86. *        conditions aux limites           *
  87. *=========================================*
  88. CL1 = BLOQ UZ DAB ;
  89. CL2 = BLOQ UZ DCD ;
  90. CL  = CL1 ET CL2  ;
  91.  
  92. *=========================================*
  93. *         MODE - materiau                 *
  94. *=========================================*
  95. MOD = MODE SUR MECANIQUE ELASTIQUE 
  96.                PLASTIQUE CINEMATIQUE QUA4;
  97.  
  98. MAT = MATE MOD YOUN 200E3 NU 0.3 SIGY 300
  99.                H 2000     ALPH 1.2E-5 'TALP' 0. 'TREF' 0.;
  100.  
  101. *=========================================*
  102. * champs relatifs a la procedure PASAPAS  *
  103. *=========================================*
  104. * Liste des pseudo-temps ou l'on desire des
  105. * resultats
  106.  
  107. t0 = 0. ; 
  108. t1 = 1. ;
  109. LISTE = PROG T0 T1 ;
  110.  
  111. * Champs-point de temperature aux instants
  112. * t0 et t1.
  114. TEMP . 0 = t0;
  115. TEMP . 1 = t1;
  116.  
  117. TEMPER = TABL;
  118. TEMPER . 0 = MANU CHPO SUR 1 T 0.  ;
  119. TEMPER . 1 = MANU CHPO SUR 1 T 500.;
  120. CHA2       = CHAR 'T' TEMP TEMPER;
  121.  
  122. * Chargement mecanique de la structure:
  123. * chargement nul.
  124. LFORCE= PROG 0. 0.;
  125. FOR   = FORCE FR 1000. PC;
  126. EVOL1 = EVOL MANU TEMPS LISTE FORCE LFORCE;
  127. CHA1  = CHAR 'MECA' FOR EVOL1;
  128.  
  129. * parametres pour la procedure PASAPAS
  130. TABDEF          = TABL ;
  131. TABDEF .'DEFTO' = VRAI;
  132. TABDEF .'DEFIN' = VRAI;
  133.  
  134. TAB1                       = TABL  ;
  135. TAB1.'CHARGEMENT'          =(CHA1 ET CHA2);
  136. TAB1.'MODELE'              = MOD   ;
  137. TAB1.'CARACTERISTIQUES'    = MAT   ;
  138. TAB1.'BLOCAGES_MECANIQUES' = CL    ;
  139. TAB1.'TEMPS_CALCULES'      = LISTE ;
  140. TAB1.'TEMPS_SAUVES'        = LISTE ;
  141. tab1.'MES_SAUVEGARDES'     = TABDEF;
  142.  
  143. * procedure PASAPAS
  144. PASAPAS TAB1;
  145.  
  146. *=========================================*
  147. *       extractions des resultats         *
  148. *=========================================*
  149.  
  150. * Contraintes
  151. nn = 1;
  152. CHEL1 = TAB1. 'CONTRAINTES'.nn;
  153.  
  154. CHAM2 = CHAN NOEUD MOD CHEL1;
  155. CHPO1 = CHAN CHPO  MOD CHAM2;
  156.  
  157. * Extraction du la contrainte SMZZ
  158. SMZZ1 = EXTR CHPO1 SMZZ PA;
  159. SMZZ2 = -309.0 ;
  160.  
  161. * Deformations plastiques
  162. CHEL10 = TAB1. 
  163.            'DEFORMATIONS_INELASTIQUES' .nn;
  164.  
  165. * Interpolations aux noeuds du maillage
  166. * Utilisation des nouveaux chamelems.
  167. CHAM20 = CHAN NOEUD MOD CHEL10;
  168. CHPO10 = CHAN CHPO  MOD CHAM20;
  169.  
  170. * Extraction de la deformation plastique
  171. * EPZZ
  172. EPZZ1 = EXTR CHPO10 EIZZ PA;
  173. EPZZ2 = -.4455E-2 ;
  174.  
  175. *=========================================*
  176. *       affichage des resultats           *
  177. *=========================================*
  179. MESS ' RESULTATS ';
  180. MESS ' --------- ';
  182.  
  183. MESS '   1) Contraintes ' ;
  184. MESS '      ----------- ' ;
  186.  
  187. MESS ' Contrainte theorique smzz:'
  188.                               SMZZ2 'Mpa' ;
  189. MESS ' Contrainte calculee smzz :'
  190.                               SMZZ1 'Mpa' ;
  191. ERG1 = (ABS((SMZZ2-SMZZ1)/SMZZ2)) * 100. ;
  192. MESS '       Soit un ecart de  :' ERG1 '%';
  194. MESS '   2) Deformations plastiques ';
  195. MESS '      ----------------------- ' ;
  197.  
  198. MESS ' Deformation plastique theorique:'
  199.                                      EPZZ2;
  200. MESS ' Deformation plastique calculee :'
  201.                                      EPZZ1;
  202. ERG2 = (ABS((EPZZ2-EPZZ1)/EPZZ2)) * 100. ;
  204. '             Soit un ecart de :' ERG2 '%';
  205.  
  207. MESS '   Maillage ';
  208. MESS '   -------- ';
  209. MESS '   Nombre de noeuds        : ' NOSUR;
  210. MESS '   Nombre d elements (qua4): ' ELSUR;
  212.  
  213. SI (NEG GRAPH 'N');
  214.   DEP1 = TAB1. 'DEPLACEMENTS' . nn;
  215.   DEF0 = DEFO SUR DEP1 0.;
  216.   DEF1 = DEFO SUR DEP1 10. BLEU;
  217.   TITR 'THPL1 : DEFORMEE';
  218.   TRAC (DEF0 ET DEF1);
  219. FINSI;
  220.  
  221. *=========================================*
  222. *       code fonctionnement et fin        *
  223. *=========================================*
  224.  
  225. * Maximum des erreurs relatives
  226. ERGMAX = MAXI (PROG ERG1 ERG2);
  227.  
  228. * Il y a erreur si un des pourcentage erg*
  229. * est superieur à 5%.
  230.  
  231. SI (ERGMAX <EG 0.03);
  232.      ERRE  0 ;
  233. SINON;
  234.      ERRE  5 ;
  235. FINSI;
  236.  
  238. TEMPS ;
  240.  
  241. FIN;
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  

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