Télécharger thme2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : thme2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *----------------------------------------------------------------
  5. *
  6. * Réduction d'un jeu sous l'action d'une sollicitation
  7. * thermique en régime transitoire
  8. *
  9. *----------------------------------------------------------------
  10. *
  11. * Calcul thermo-mécanique ( mécanique et thermique linéaire )
  12. * Utilisation de la procédure PASAPAS
  13. * Utilisation de la procédure utilisateur perso1.
  14. *
  15. *
  16. * | x1 = 0.1 m x2 = 0.1002 m |
  17. * | T1 = 500°C T2 = ? |
  18. * |_______________________________ ____________________________|
  19. * | | | |
  20. * | T_initial = 20°C | | T_initial = 20°C |
  21. * |_______________________________| |____________________________|
  22. * | x0 = 0 m jeu x3 = 0.2002 m |
  23. * | T0 = 20°C T3 = 20°C |
  24. * | Encastrement Encastrement |
  25. *
  26. *
  27. * Le jeu sera considéré comblé quand la distance entre les
  28. * barreaux est inférieure à 0.1mm
  29. *
  30. * on étudie la diffusion thermique dans les barreaux et leurs
  31. * dilatations.
  32. *
  33. * caracteristiques des barreaux :
  34. *
  35. * l = 0.1m
  36. * K = 14.6
  37. * C = 460
  38. * RHO = 7800
  39. * ALPHA = 1E-5
  40. * YOUNG = 2.E11
  41. * NU = 0.00001
  42. *
  43. * Mailles de 0.01 m de cote
  44. *
  45. *----------------------------------------------------------------
  46. *
  47. OPTION ECHO 0;
  48. OPTION DIME 2 ELEM QUA8 MODE PLAN DEFO;
  49. *
  50. *-------------------- Creation du maillage ----------------------
  51. *
  52. DENS 1;
  53. P1 = 0. 0. ;
  54. P2 = 0.1 0. ;
  55. P3 = 0.1 0.01 ;
  56. P4 = 0. 0.01 ;
  57. P5 = 0.1002 0.;
  58. P6 = 0.2002 0.;
  59. P7 = 0.2002 0.01;
  60. P8 = 0.1002 0.01;
  61. P9 = 0.12 0.;
  62. D1 = DROI 10 P1 P2 ;
  63. D2 = DROI 1 P2 P3 ;
  64. D3 = DROI 10 P3 P4 ;
  65. D4 = DROI 1 P4 P1 ;
  66. D5 = DROI 10 P5 P6 ;
  67. D6 = DROI 1 P6 P7 ;
  68. D7 = DROI 10 P7 P8 ;
  69. D8 = DROI 1 P8 P5 ;
  70. L1 = D1 ET D2 ET D3 ET D4;
  71. L2 = D5 ET D6 ET D7 ET D8;
  72. SU1 = SURF L1 PLANE;
  73. SU2 = SURF L2 PLANE;
  74. P9 = SU2 POINT PROCHE P9;
  75. *
  76. *-------------------- Définition du modèle ----------------------
  77. *
  78. MOD1 = MODE (SU1 ET SU2) THERMIQUE ISOTROPE CONS UN;
  79. MOD3 = MODE (SU1 ET SU2) MECANIQUE ELASTIQUE ISOTROPE CONS DEUX ;
  80. MOD = MOD1 et MOD3;
  81. *
  82. *------------------- Définition des matériaux -------------------
  83. *
  84. MAT1 = MATE MOD1 K 14.6 'C' 460. RHO 7800. ;
  85. MAT3 = MATE MOD3 ALPHA 1E-5 YOUN 2E+11 NU 0.00001;
  86. MAT = MAT1 ET MAT3;
  87. *
  88. *------------------- blocages mécaniques ------------------------
  89. *
  90. CL1 = BLOQ UY (SU1 ET SU2);
  91. CL2 = BLOQ UX D4;
  92. CL3 = BLOQ UX D6;
  93. CL_ME = CL1 ET CL2 ET CL3;
  94. *
  95. *------------------- blocages thermiques ------------------------
  96. *
  97. CL4 = BLOQ T D4;
  98. CL5 = BLOQ T D2;
  99. CL6 = BLOQ T D6;
  100. TEM1 = DEPI (CL4 ET CL6) 20.;
  101. TEM2 = DEPI CL5 500.;
  102. CL_TH = CL4 ET CL5 ET CL6;
  103. *
  104. *------------------- températures imposées ----------------------
  105. *
  106. TEM1 = DEPI (CL4 ET CL6) 20.;
  107. TEM2 = DEPI CL5 500.;
  108. TIM = TEM1 ET TEM2;
  109. EV1 = EVOL MANU TEMPS (PROG 0. 200.) y (prog 1. 1.);
  110. CHA = CHAR 'TIMP' TIM EV1;
  111. *
  112. *------------------- température initiale -----------------------
  113. *
  114. TEM0 = MANU CHPO (SU1 ET SU2) 1 T 20.;
  115. *
  116. *--------------- creation de la table pour pasapas ------------
  117. *
  118. ta = table;
  119. ta.modele = mod;
  120. ta.caracteristiques = mat;
  121. ta.blocages_mecaniques = cl_me;
  122. ta.blocages_thermiques = cl_th;
  123. ta.temperatures = table;
  124. ta.temperatures . 0 = tem0;
  125. ta.chargement = cha;
  126. ta.procedure_thermique = DUPONT;
  127. ta.procedure_perso1 = vrai;
  128. ta.ok = faux;
  129. ta.temps_calcules = (prog 0. pas 5. 80.) et (prog 81. pas 1. 89.)
  130. et (prog 90. pas 5. 150.);
  131. *
  132. *------------ procédure utilisateur perso1 --------------------
  133. *
  134. debproc perso1 TTU*TABLE;
  135. dep = ttu.'ESTIMATION'.'DEPLACEMENTS';
  136. WTAB=TTU.'WTABLE';
  137. x = extr dep UX P2;
  138. list x;
  139. si ((WTAB.'TEMPS0' > 84.5) et ( WTAB.'TEMPS0' < 89.5));
  140. list ttu.'ESTIMATION'.'TEMPERATURES';
  141. finsi;
  142.  
  143. si (('>' x 0.0001) et (non ttu.ok));
  144. ttu.loup = WTAB.'TEMPS0';
  145. mess ' *************************';
  146. mess ' mise a jour du au contact pour le temps : ' ttu.loup ;
  147. mess ' *************************';
  148. CL7 = bloq T D8;
  149. tim = depi cl7 500.;
  150. WTAB.'BLOCAGES_THERMIQUES' = (WTAB.'BLOCAGES_THERMIQUES' ) et cl7;
  151. EV1 = EVOL MANU TEMPS (PROG 0. 200.) y (prog 1. 1.);
  152. CHA7 = CHAR 'TIMP' TIM EV1;
  153. WTAB.'CHARGEMENT_I' = WTAB.'CHARGEMENT_I' et cha7;
  154. WTAB.'CHARGEMENT' = WTAB.'CHARGEMENT' et cha7;
  155. ttu.ok = vrai;
  156. finsi;
  157. finproc;
  158. *
  159. *------------ Appel à la procédure PASAPAS --------------------
  160. *
  161. *
  162. *----------- Depouillement ------------------------------------
  163. *
  164. a = (14.6 / 7800.)/460.;
  165. tpr1 = prog;
  166. tpr3 = prog;
  167. tep_th = prog;
  168. x = 0.12;
  169. er1 = 0.;
  170. nb = 0;
  171. idim = (dime (ta.temps)) - 1;
  172. *
  173. repeter bou idim;
  174. ite = &bou;
  175. te = ta.temps.ite;
  176. tem = ta.temperatures.ite;
  177. tep1 = extr tem T p9;
  178. tpr1 = tpr1 et (prog tep1);
  179. tpr3 = tpr3 et (prog te);
  180. si ('>' te (ta.loup)) ;
  181. co1 = erf((x - .1002)/(2 * (( a * (te - (ta.loup))**0.5))));
  182. co2 = erf(0.1/(2.*(( a * (te - (ta.loup)))**0.5)));
  183. sol_te = 500. - ((500. - 20.)*co1/co2);
  184. er1 = er1 + (abs ((sol_te - tep1)/sol_te) * 100.);
  185. nb = nb + 1;
  186. sinon;
  187. sol_te = 20.;
  188. finsi;
  189. tep_th = tep_th et (prog sol_te);
  190. fin bou;
  191. *
  192. er1 = er1/nb;
  193. mess 'erreur sur la température :'er1'%';
  194. *
  195. *lvar = TEXT 'EVOLUTION DE LA TEMPERATURE AU POINT 0.12 m';
  196. *titre lvar;
  197. *ev_te = evol manu temps tpr3 temperature tep_th;
  198. *ev_tet = evol manu temps tpr3 temperature tpr1;
  199. *dess (ev_tet et ev_te) ;
  200. *
  201. si ( er1 > 5. ) ;
  202. mess 'Résultat incorrect' ;
  203. erreur 5 ;
  204. finsi ;
  205. *
  206. fin;
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  

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