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.
  102. TEMP . 0 = T0;
  103. TEMP . 1 = T1;
  104.  
  105. TEMPER = TABL;
  106. TEMPER . 0 = MANU CHPO SUR1 1 T 0. ;
  107. TEMPER . 1 = MANU CHPO SUR1 1 T 500.;
  108. CHA2 = CHAR 'T' TEMP TEMPER;
  109.  
  110. * Chargement mecanique : chargement nul.
  111. LFORCE = PROG 0. 0.;
  112. FOR1 = FORCE FY 1. PA;
  113. CHPO1 = RESU FOR1;
  114. EVOL1 = EVOL MANU TEMPS LIS1 FORCE LFORCE;
  115. CHA1 = CHAR 'MECA' CHPO1 EVOL1;
  116.  
  117. * Parametres pour la procedure PASAPAS
  118. TABDEF = TABL ;
  119. TABDEF .'DEFTO' = VRAI;
  120. TABDEF .'DEFIN' = VRAI;
  121.  
  122. TAB1 = TABL ;
  123. TAB1.'CHARGEMENT' =(CHA1 ET CHA2);
  124. TAB1.'MODELE' = MO ;
  125. TAB1.'CARACTERISTIQUES' = MAT ;
  126. TAB1.'BLOCAGES_MECANIQUES' = CL ;
  127. TAB1.'TEMPS_CALCULES' = LIS1 ;
  128. TAB1.'TEMPS_SAUVES' = LIS1 ;
  129. tab1.'MES_SAUVEGARDES' = TABDEF;
  130.  
  131. * procedure PASAPAS
  132. PASAPAS TAB1;
  133.  
  134. *=========================================*
  135. * Extraction des resultats *
  136. *=========================================*
  137.  
  138. * Contraintes
  139. NN = 1;
  140. CHEL1 = TAB1.CONTRAINTES.NN;
  141.  
  142. * Interpolations aux noeuds :
  143. CHAM2 = CHAN NOEUD MO CHEL1;
  144. CHPO1 = CHAN CHPO MO CHAM2;
  145.  
  146. * Extractionsn des contraintes smxx et smyy
  147. SMXX1 = EXTR CHPO1 'SMXX' PA;
  148. SMYY1 = EXTR CHPO1 'SMYY' PA;
  149. SMXX1 = SMXX1 / 1.E6;
  150. SMYY1 = SMYY1 / 1.E6;
  151.  
  152. SMXX2 = -125.9;
  153. SMYY2 = -125.9;
  154.  
  155. * Calcul de l'erreur commise
  156. ERGSXX = 100.*(ABS((SMXX1-SMXX2)/ SMXX2 ));
  157. ERGSYY = 100.*(ABS((SMYY1-SMYY2)/ SMYY2 ));
  158.  
  159. * Deformations plastiques.
  160. CHEL10 = TAB1.DEFORMATIONS_INELASTIQUES.NN;
  161.  
  162. * Interpolations aux noeuds du maillage
  163. CHAM20 = CHAN NOEUD MO CHEL10;
  164. CHPO10 = CHAN CHPO MO CHAM20;
  165.  
  166. * Extractions des deformations plastiques
  167. EPXX1 = EXTR CHPO10 EIXX PA;
  168. EPYY1 = EXTR CHPO10 EIYY PA;
  169.  
  170. EPXX2 = -0.6399E-2;
  171. EPYY2 = -0.6399E-2;
  172.  
  173. * Calcul de l'erreur commise
  174. ERGEXX = 100.*(ABS((EPXX1-EPXX2)/ EPXX2 ));
  175. ERGEYY = 100.*(ABS((EPYY1-EPYY2)/ EPYY2 ));
  176.  
  177. * Nombre d'elements et noeuds du maillage
  178. NOSUR1 = NBNO SUR1 ;
  179. ELSUR1 = NBEL SUR1 ;
  180.  
  181. *=========================================*
  182. * Affichage des resultats *
  183. *=========================================*
  184. MESS ' RESULTATS ';
  185. MESS ' --------- ';
  186.  
  187. MESS ' Contraintes interpolees (MPa)';
  188. MESS ' -----------------------------';
  189. ' Contrainte theorique smxx :' SMXX2 'MPa';
  190. ' Contrainte calculee smxx :' SMXX1 'MPa';
  191. MESS ' Soit un ecart de :' ERGSXX '%';
  192.  
  193. ' Contrainte theorique smyy :' SMYY2 'MPa';
  194. ' Contrainte calculee smyy :' SMYY1 'MPa';
  195. MESS ' Soit un ecart de :' ERGSYY '%';
  196.  
  197. MESS 'Deformations plastiques interpolees';
  198. MESS '-----------------------------------';
  199.  
  200. ' Deformation plastique theorique epxx :'
  201. EPXX2;
  202. ' Deformation plastique calculee epxx :'
  203. EPXX1;
  204. MESS ' Soit un ecart de :' ERGEXX '%';
  205.  
  206. ' Deformation plastique theorique epyy :'
  207. EPYY2;
  208. ' Deformation plastique calculee epyy :'
  209. EPYY1;
  210. MESS ' Soit un ecart de :' ERGEYY '%';
  211.  
  212. MESS ' Maillage ';
  213. MESS ' -------- ';
  214. MESS ' Nombre de noeuds : ' NOSUR1 ;
  215. MESS ' Nombre d elements : ' ELSUR1 ;
  216.  
  217. SI (NEG GRAPH 'N');
  218. DEP1 = TAB1. 'DEPLACEMENTS' . NN;
  219. DEF0 = DEFO SUR1 DEP1 0.;
  220. DEF1 = DEFO SUR1 DEP1 10. BLEU;
  221. TITR 'THPL2 : DEFORMEE';
  222. TRAC (DEF0 ET DEF1);
  223.  
  224. TITR 'THPL2 : Deformation Plastique';
  225. TRAC CHEL10 MO;
  226. FINSI;
  227.  
  228. *=========================================*
  229. * Code fonctionnement et fin *
  230. *=========================================*
  231. * Maximum des erreurs relatives
  232. ERGMAX = MAXI
  233. (PROG ERGSXX ERGSYY ERGEXX ERGEYY);
  234.  
  235. SI (ERGMAX <EG 5.);
  236. ERRE 0 ;
  237. SINON;
  238. ERRE 5 ;
  239. FINSI;
  240. TEMPS ;
  241.  
  242. FIN;
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  

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