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 TRAC 'PSC' ;
  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 'TALP' 0. 'TREF' 20. 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. CL_TH = CL4 ET CL5 ET CL6;
  101. *
  102. *------------------- températures imposées ----------------------
  103. *
  104. TEM1 = DEPI (CL4 ET CL6) 20.;
  105. TEM2 = DEPI CL5 500.;
  106. TIM = TEM1 ET TEM2;
  107. EV1 = EVOL MANU TEMPS (PROG 0. 200.) y (prog 1. 1.);
  108. CHA = CHAR 'TIMP' TIM EV1;
  109. *
  110. *------------------- température initiale -----------------------
  111. *
  112. TEM0 = MANU CHPO (SU1 ET SU2) 1 T 20.;
  113. *
  114. *--------------- creation de la table pour pasapas ------------
  115. *
  116. ta = table;
  117. ta.modele = mod;
  118. ta.caracteristiques = mat;
  119. ta.blocages_mecaniques = cl_me;
  120. ta.blocages_thermiques = cl_th;
  121. ta.temperatures = table;
  122. ta.temperatures . 0 = tem0;
  123. ta.chargement = cha;
  124. ta.procedure_thermique = DUPONT;
  125. ta.procedure_perso1 = vrai;
  126. ta.ok = faux;
  127. ta.temps_calcules = prog 0. pas 5. 80. pas 0.5 100. pas 1. 150.;
  128. *
  129. *------------ procédure utilisateur perso1 --------------------
  130. *
  131. debproc perso1 TTU*TABLE;
  132. dep = ttu.'ESTIMATION'.'DEPLACEMENTS';
  133. WTAB=TTU.'WTABLE';
  134. x = extr dep UX P2;
  135. list x;
  136. si ((TTU.'ESTIMATION'.'TEMPS' > 84.5) et ( TTU.'ESTIMATION'.'TEMPS' < 89.5));
  137. list ttu.'ESTIMATION'.'TEMPERATURES';
  138. finsi;
  139.  
  140. si (('>' x 0.0001) et (non ttu.ok));
  141. ttu.loup = TTU.'ESTIMATION'.'TEMPS';
  142. mess ' *************************';
  143. mess ' mise a jour du au contact pour le temps : ' ttu.loup ;
  144. mess ' *************************';
  145. CL7 = bloq T D8;
  146. tim = depi cl7 500.;
  147. WTAB.'BLOCAGES_THERMIQUES' = (WTAB.'BLOCAGES_THERMIQUES' ) et cl7;
  148. EV1 = EVOL MANU TEMPS (PROG 0. 200.) y (prog 1. 1.);
  149. CHA7 = CHAR 'TIMP' TIM EV1;
  150. WTAB.'CHARGEMENT_I' = WTAB.'CHARGEMENT_I' et cha7;
  151. WTAB.'CHARGEMENT' = WTAB.'CHARGEMENT' et cha7;
  152. ttu.ok = vrai;
  153. finsi;
  154. finproc;
  155. *
  156. *------------ Appel à la procédure PASAPAS --------------------
  157. *
  158. *
  159. *----------- Depouillement ------------------------------------
  160. *
  161. a = (14.6 / 7800.)/460.;
  162. tpr1 = prog;
  163. tpr3 = prog;
  164. tep_th = prog;
  165. x = 0.12;
  166. er1 = 0.;
  167. nb = 0;
  168. idim = (dime (ta.temps)) - 1;
  169. *
  170. repeter bou idim;
  171. ite = &bou;
  172. te = ta.temps.ite;
  173. tem = ta.temperatures.ite;
  174. tep1 = extr tem T p9;
  175. tpr1 = tpr1 et (prog tep1);
  176. tpr3 = tpr3 et (prog te);
  177. si ('>' te (ta.loup)) ;
  178. co1 = erf((x - .1002)/(2 * (( a * (te - (ta.loup))**0.5))));
  179. co2 = erf(0.1/(2.*(( a * (te - (ta.loup)))**0.5)));
  180. sol_te = 500. - ((500. - 20.)*co1/co2);
  181. er1 = er1 + (abs ((sol_te - tep1)/sol_te) * 100.);
  182. nb = nb + 1;
  183. sinon;
  184. sol_te = 20.;
  185. finsi;
  186. tep_th = tep_th et (prog sol_te);
  187. fin bou;
  188. *
  189. er1 = er1/nb;
  190. mess 'erreur sur la température :'er1'%';
  191. *
  192. lvar = TEXT 'EVOLUTION DE LA TEMPERATURE AU POINT 0.12 m';
  193. titre lvar;
  194. ev_te = evol manu temps tpr3 temperature tep_th;
  195. ev_tet = evol manu temps tpr3 temperature tpr1;
  196. dess (ev_tet et ev_te) 'YBOR' 0. 250. ;
  197.  
  198. dess (EVOL 'TEMP' ta 'DEPLACEMENTS' 'UX' P2);
  199. *
  200. si ( er1 > 5. ) ;
  201. mess 'Resultat incorrect' ;
  202. erreur 5 ;
  203. finsi ;
  204. *
  205. fin;
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  

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