Télécharger la plaquette Cast3M


fr en

Télécharger thpl2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  thpl2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/06/13    PROVENANCE : TEST
  5.  
  6.  
  7. *******************************************
  8. *    Test thpl2.dgibi: jeux de données    *
  9. *    ---------------------------------    *
  10. *******************************************
  11.  
  12. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  13. * SINON SI GRAPH DIFFERENT DE N TOUS
  14. * LES GRAPHIQUES SONT AFFICHES
  15.  
  16. GRAPH = 'N' ;
  17.  
  18. SAUT PAGE;
  19. SI (NEG GRAPH 'N') ;
  20.   OPTI ECHO 1   ;
  21.   OPTI TRAC PSC   ;
  22. SINO ;
  23.   OPTI ECHO 0   ;
  24. FINSI ;
  25.  
  26. SAUT PAGE;
  27.  
  28. *******************************************
  29. *                  THPL2                  *
  30. *                                         *
  31. *   PLAQUE CARREE ENCASTREE SOUMIS A UN   *
  32. *           CHAMPS DE TEMPERTURE          *
  33. *                                         *
  34. *   Test de la commission V.P.C.S.        *
  35. *   Groupe Thermo-plasticite              *
  36. *   Statique non lineaire                 *
  37. *                                         *
  38. *   Une plaque carree de 10 cm de cote est*
  39. *   encastree et soumise a un champs      *
  40. *   constant de temperature de 500 degres *
  41. *   celsius.                              *
  42. *                                         *
  43. *  On se propose de calculer              *
  44. *   - contraintes smxx et smyy,           *
  45. *   - deformations plastiques epxx et epyy*
  46. *                                         *
  47. *      Les valeurs trouvees sont ensuites *
  48. *   comparees aux valeurs de reference    *
  49. *   obtenues analytiquement               *
  50. *******************************************
  51. OPTION DIME 2 ;
  52. OPTION MODE PLAN CONT;
  53. OPTION ELEM QUA4;
  54.  
  55. *=========================================*
  56. *          geometrie : maillage           *
  57. *=========================================*
  58.  
  59. * Dimensions en metres
  60. CO = 0.01;
  61.  
  62. PA = 0. 0.;
  63. PB = CO 0.;
  64. PC = CO CO;
  65. PD = 0. CO;
  66.  
  67. LAB = PA DROIT 8 PB;
  68. LBC = PB DROIT 8 PC;
  69. LCD = PC DROIT 8 PD;
  70. LDA = PD DROIT 8 PA;
  71.  
  72. LIG = LAB ET LBC ET LCD ET LDA;
  73.  
  74. SUR1 = DALL LAB LBC LCD LDA PLAN;
  75.  
  76. SI(NEG GRAPH 'N');
  77.   TITR 'THPL2 : MAILLAGE';
  78.   TRACE 'QUAL' SUR1;
  79. FINSI;
  80.  
  81. *=========================================*
  82. *        conditions aux limites           *
  83. *=========================================*
  84. CL = BLOQUE LIG DEPL ROTA;
  85.  
  86. *=========================================*
  87. *           MODE - materiau     
  88. *=========================================*
  89. MO  = MODE SUR1 MECANIQUE ELASTIQUE  PLASTIQUE CINEMATIQUE QUA4 ;
  90. MAT = MATE MO YOUN 163750.0E6 NU 0.3 ALPHA 13.875E-6 SIGY 100E6 H 2000E6 'TREF' 0. 'TALP' 0. ;
  91.  
  92. *=========================================*
  93. * champs relatifs a la procedure  PASAPAS *
  94. *=========================================*
  95. * Liste des pseudo-temps ou l'on desire des
  96. * resultats.
  97. T0 = 0. ;
  98. T1 = 1.;
  99. LIS1 =    PROG T0  T1;
  100.  
  101. * Champs points de temperature.
  103. TEMP . 0 = T0;
  104. TEMP . 1 = T1;
  105.  
  106. TEMPER = TABL;
  107. TEMPER . 0 = MANU CHPO SUR1 1 T 0.  ;
  108. TEMPER . 1 = MANU CHPO SUR1 1 T 500.;
  109. CHA2       = CHAR 'T' TEMP TEMPER;
  110.  
  111. * Chargement mecanique : chargement nul.
  112. LFORCE = PROG 0. 0.;
  113. FOR1   = FORCE FY 1. PA;
  114. CHPO1  = RESU FOR1;
  115. EVOL1  = EVOL MANU TEMPS LIS1 FORCE LFORCE;
  116. CHA1   = CHAR 'MECA' CHPO1 EVOL1;
  117.  
  118. * Parametres pour la procedure PASAPAS
  119. TABDEF          = TABL ;
  120. TABDEF .'DEFTO' = VRAI;
  121. TABDEF .'DEFIN' = VRAI;
  122.  
  123. TAB1                       = TABL  ;
  124. TAB1.'CHARGEMENT'          =(CHA1 ET CHA2);
  125. TAB1.'MODELE'              = MO   ;
  126. TAB1.'CARACTERISTIQUES'    = MAT   ;
  127. TAB1.'BLOCAGES_MECANIQUES' = CL    ;
  128. TAB1.'TEMPS_CALCULES'      = LIS1 ;
  129. TAB1.'TEMPS_SAUVES'        = LIS1 ;
  130. tab1.'MES_SAUVEGARDES'     = TABDEF;
  131.  
  132. * procedure PASAPAS
  133. PASAPAS TAB1;
  134.  
  135. *=========================================*
  136. *       Extraction des resultats          *
  137. *=========================================*
  138.  
  139. * Contraintes
  140. NN = 1;
  141. CHEL1 = TAB1.CONTRAINTES.NN;
  142.  
  143. * Interpolations aux noeuds :
  144. CHAM2 = CHAN NOEUD MO CHEL1;
  145. CHPO1 = CHAN CHPO  MO CHAM2;
  146.  
  147. * Extractionsn des contraintes smxx et smyy
  148. SMXX1 = EXTR CHPO1 'SMXX' PA;
  149. SMYY1 = EXTR CHPO1 'SMYY' PA;
  150. SMXX1 = SMXX1 / 1.E6;
  151. SMYY1 = SMYY1 / 1.E6;
  152.  
  153. SMXX2 = -125.9;
  154. SMYY2 = -125.9;
  155.  
  156. * Calcul de l'erreur commise
  157. ERGSXX = 100.*(ABS((SMXX1-SMXX2)/ SMXX2 ));
  158. ERGSYY = 100.*(ABS((SMYY1-SMYY2)/ SMYY2 ));
  159.  
  160. * Deformations plastiques.
  161. CHEL10 = TAB1.DEFORMATIONS_INELASTIQUES.NN;
  162.  
  163. * Interpolations aux noeuds du maillage
  164. CHAM20 = CHAN NOEUD MO CHEL10;
  165. CHPO10 = CHAN CHPO  MO CHAM20;
  166.  
  167. * Extractions des deformations plastiques
  168. EPXX1 = EXTR CHPO10 EIXX PA;
  169. EPYY1 = EXTR CHPO10 EIYY PA;
  170.  
  171. EPXX2 = -0.6399E-2;
  172. EPYY2 = -0.6399E-2;
  173.  
  174. * Calcul de l'erreur commise
  175. ERGEXX = 100.*(ABS((EPXX1-EPXX2)/ EPXX2 ));
  176. ERGEYY = 100.*(ABS((EPYY1-EPYY2)/ EPYY2 ));
  177.  
  178. * Nombre d'elements et noeuds du maillage
  179. NOSUR1 = NBNO SUR1 ;
  180. ELSUR1 = NBEL SUR1 ;
  181.  
  182. *=========================================*
  183. *        Affichage des resultats          *
  184. *=========================================*
  185. MESS ' RESULTATS ';
  186. MESS ' --------- ';
  188.  
  189. MESS '   Contraintes interpolees (MPa)';
  190. MESS '   -----------------------------';
  193. ' Contrainte theorique smxx :' SMXX2 'MPa';
  195. ' Contrainte calculee  smxx :' SMXX1 'MPa';
  196. MESS '      Soit un ecart de :' ERGSXX '%';
  198.  
  200. ' Contrainte theorique smyy :' SMYY2 'MPa';
  202. ' Contrainte calculee  smyy :' SMYY1 'MPa';
  203. MESS '      Soit un ecart de :' ERGSYY '%';
  205.  
  206. MESS 'Deformations plastiques interpolees';
  207. MESS '-----------------------------------';
  209.  
  211. ' Deformation plastique theorique epxx :' 
  212. EPXX2;
  214. ' Deformation plastique calculee  epxx :' 
  215. EPXX1;
  216. MESS '      Soit un ecart de :' ERGEXX '%';
  218.  
  220. ' Deformation plastique theorique epyy :'
  221. EPYY2;
  223. ' Deformation plastique calculee  epyy :' 
  224. EPYY1;
  225. MESS '      Soit un ecart de :' ERGEYY '%';
  227.  
  228. MESS '   Maillage ';
  229. MESS '   -------- ';
  230. MESS '     Nombre de noeuds  : ' NOSUR1 ;
  231. MESS '     Nombre d elements : ' ELSUR1 ;
  233.  
  234. SI (NEG GRAPH 'N');
  235.   DEP1 = TAB1. 'DEPLACEMENTS' . NN;
  236.   DEF0 = DEFO SUR1 DEP1 0.;
  237.   DEF1 = DEFO SUR1 DEP1 10. BLEU;
  238.   TITR 'THPL2 : DEFORMEE';
  239.   TRAC (DEF0 ET DEF1);
  240.  
  241.   TITR 'THPL2 : Deformation Plastique';
  242.   TRAC CHEL10 MO;
  243. FINSI;
  244.  
  245. *=========================================*
  246. *      Code fonctionnement et fin         *
  247. *=========================================*
  248. * Maximum des erreurs relatives
  249. ERGMAX = MAXI 
  250.         (PROG ERGSXX ERGSYY ERGEXX ERGEYY);
  251.  
  252. SI (ERGMAX <EG 5.);
  253.      ERRE  0 ;
  254. SINON;
  255.      ERRE  5 ;
  256. FINSI;
  258. TEMPS ;
  259.  
  260. FIN;
  261.  
  262.  
  263.  
  264.  
  265.  
  266.  
  267.  
  268.  
  269.  

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