发明名称 一种大型铸件铸造过程表面温度变化间接测量的方法
摘要 本发明公开了属于砂型铸造过程中铸件温度变化测量技术领域的一种大型铸件铸造过程表面温度变化间接测量的方法。它包括以下步骤:1)在砂型中距离铸件表面按不同距离处布置第1热电偶~第n热电偶;2)在时刻t,由第1热电偶~第n热电偶分别测得温度<img file="DDA0000113829350000011.GIF" wi="156" he="52" />3)计算时刻t经过各个热电偶的热通量;4)计算时刻t通过铸件与砂型的界面的热通量;5)在时刻t计算砂型内表面温度<img file="DDA0000113829350000012.GIF" wi="64" he="52" />6)计算时刻t铸件表面温度<img file="DDA0000113829350000013.GIF" wi="137" he="58" />本发明的有益效果是:K型热电偶只测量多点砂型温度,不用接触铸件,稳定可靠,不影响铸件质量,经济方便。可以在测得砂型温度变化的同时,计算出铸件温度变化,方便工厂掌握铸件温度变化规律,进而调整一些工艺方案,例如准确地把握落砂时间。
申请公布号 CN102407295B 申请公布日期 2013.07.03
申请号 CN201110388031.5 申请日期 2011.11.29
申请人 清华大学;宁夏共享铸钢有限公司 发明人 康进武;游锐;聂刚;郝小坤;龙海敏;黄天佑;付龙;马斌
分类号 B22D2/00(2006.01)I;B22D46/00(2006.01)I 主分类号 B22D2/00(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 朱琨
主权项 1.一种大型铸件铸造过程表面温度变化间接测量的方法,其特征在于,分为以下步骤:1)在砂型中距离铸件表面按不同距离依次布置第1热电偶(TC<sub>1</sub>)~第n热电偶(TC<sub>n</sub>),它们距离铸件表面的距离分别为d<sub>1</sub>~d<sub>n</sub>,第1热电偶(TC<sub>1</sub>)~第n热电偶(TC<sub>n</sub>)安装在垂直于铸件表面的同一水平线上;3≤n≤10;2)在时刻t,由第1热电偶(TC<sub>1</sub>)~第n热电偶(TC<sub>n</sub>)分别测得温度<img file="FDA00002817211000011.GIF" wi="163" he="57" />并由温度记录仪记录以上时刻与温度;3)计算通过第i热电偶(TC<sub>i</sub>)到第i+1热电偶(TC<sub>i+1</sub>)的热通量<img file="FDA00002817211000012.GIF" wi="108" he="61" />其计算公式如下:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>q</mi><mrow><mi>i</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>t</mi></msubsup><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><msubsup><mi>T</mi><mi>i</mi><mi>t</mi></msubsup><mo>-</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>t</mi></msubsup><mo>)</mo></mrow><msub><mi>k</mi><mi>sand</mi></msub></mrow><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>d</mi><mi>i</mi></msub><mo>)</mo></mrow></mfrac></mrow></math>]]></maths>上式中,<img file="FDA00002817211000014.GIF" wi="178" he="57" />分别为第i热电偶(TC<sub>i</sub>)和第i+1热电偶(TC<sub>i+1</sub>)在时刻t测得的温度;K<sub>sand</sub>为砂型的导热系数;d<sub>i</sub>,d<sub>i+1</sub>分别为第i热电偶(TC<sub>i</sub>)和第i+1热电偶(TC<sub>i+1</sub>)到铸件表面的距离,i的取值为1至n-1;4)计算时刻t通过铸件与砂型的界面(1)的热通量:对<maths num="0002"><![CDATA[<math><mrow><mo>[</mo><mfrac><mrow><msub><mi>d</mi><mn>1</mn></msub><mo>+</mo><msub><mi>d</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>q</mi><mn>1,2</mn><mi>t</mi></msubsup><mo>]</mo><mo>,</mo><mo>[</mo><mfrac><mrow><msub><mi>d</mi><mn>2</mn></msub><mo>+</mo><msub><mi>d</mi><mn>3</mn></msub></mrow><mn>2</mn></mfrac><mo>,</mo><msubsup><mi>q</mi><mn>2,3</mn><mi>t</mi></msubsup><mo>]</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>[</mo><mfrac><mrow><msub><mi>d</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>d</mi><mi>n</mi></msub></mrow><mn>2</mn></mfrac><mo>,</mo><msubsup><mi>q</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi></mrow><mi>t</mi></msubsup><mo>]</mo></mrow></math>]]></maths>这n-1组数据采用最小二乘法拟合成线性模型曲线,并得到如下关系式:<maths num="0004"><![CDATA[<math><mrow><msup><mi>y</mi><mi>t</mi></msup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>b</mi><mi>k</mi><mi>t</mi></msubsup><msup><mi>x</mi><mrow><mi>n</mi><mo>-</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup></mrow></math>]]></maths>上式中,x为距离铸件表面的距离,y<sup>t</sup>为时刻t在距离铸件表面x处通过的热通量,当x=0时,<img file="FDA00002817211000021.GIF" wi="162" he="57" />为时刻t通过铸件与砂型的界面(1)的热通量<img file="FDA00002817211000022.GIF" wi="67" he="57" />采用最小二乘法得到关于<img file="FDA00002817211000023.GIF" wi="41" he="57" />的如下方程组,解方程组可得到<img file="FDA00002817211000024.GIF" wi="68" he="57" /><maths num="0005"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>s</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><msubsup><mi>b</mi><mi>k</mi><mi>t</mi></msubsup><mo>=</mo><msub><mi>s</mi><mrow><mi>j</mi><mo>,</mo><mi>y</mi></mrow></msub></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>w</mi><mi>m</mi></msub><msubsup><mi>x</mi><mi>m</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>-</mo><mi>j</mi></mrow></msubsup><msubsup><mi>x</mi><mi>m</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>-</mo><mi>k</mi></mrow></msubsup></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>,</mo><mi>y</mi></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>w</mi><mi>m</mi></msub><msubsup><mi>x</mi><mi>m</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn><mo>-</mo><mi>j</mi></mrow></msubsup><msubsup><mi>y</mi><mi>m</mi><mi>t</mi></msubsup></mrow></math>]]></maths>上述方程组中,w<sub>1</sub>设置为<img file="FDA00002817211000028.GIF" wi="156" he="99" />w<sub>2</sub>至w<sub>n-1</sub>均设置为<img file="FDA00002817211000029.GIF" wi="156" he="109" />j的取值为1至n-1,x<sub>m</sub>为<img file="FDA000028172110000210.GIF" wi="219" he="111" /><img file="FDA000028172110000211.GIF" wi="51" he="58" />为<img file="FDA000028172110000212.GIF" wi="133" he="62" />5)在时刻t砂型内表面温度<img file="FDA000028172110000213.GIF" wi="44" he="57" />的计算:对<img file="FDA000028172110000214.GIF" wi="612" he="58" /><img file="FDA000028172110000215.GIF" wi="70" he="58" />这n组数据也采用最小二乘法拟合成线性模型曲线,并得到如下关系式:<maths num="0008"><![CDATA[<math><mrow><msup><mi>T</mi><mi>t</mi></msup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>a</mi><mi>k</mi><mi>t</mi></msubsup><msup><mi>d</mi><mrow><mi>n</mi><mo>-</mo><mi>k</mi></mrow></msup></mrow></math>]]></maths>上式中,d为距离铸件表面的距离,T<sup>t</sup>为时刻t在距离铸件表面d处砂型温度,当d=0时,<img file="FDA000028172110000217.GIF" wi="139" he="56" />为时刻t砂型内表面温度<img file="FDA000028172110000218.GIF" wi="72" he="56" />采用最小二乘法得到关于<img file="FDA000028172110000219.GIF" wi="43" he="56" />的如下方程组,解方程组可得到<img file="FDA000028172110000220.GIF" wi="70" he="56" /><maths num="0009"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>u</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow></msub><msubsup><mi>a</mi><mi>k</mi><mi>t</mi></msubsup><mo>=</mo><msub><mi>u</mi><mrow><mi>l</mi><mo>,</mo><mi>T</mi></mrow></msub></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><msub><mi>u</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>m</mi></msub><msubsup><mi>d</mi><mi>m</mi><mrow><mi>n</mi><mo>-</mo><mi>l</mi></mrow></msubsup><msubsup><mi>d</mi><mi>m</mi><mrow><mi>n</mi><mo>-</mo><mi>k</mi></mrow></msubsup></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><msub><mi>u</mi><mrow><mi>l</mi><mo>,</mo><mi>T</mi></mrow></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>m</mi></msub><msubsup><mi>d</mi><mi>m</mi><mrow><mi>n</mi><mo>-</mo><mi>l</mi></mrow></msubsup><msubsup><mi>T</mi><mi>m</mi><mi>t</mi></msubsup></mrow></math>]]></maths>上述方程组中,w为权重,w<sub>1</sub>设置为<img file="FDA00002817211000031.GIF" wi="122" he="108" />w<sub>2</sub>至w<sub>n</sub>均设置为<img file="FDA00002817211000032.GIF" wi="89" he="108" />l的取值为1至n;6)时刻t铸件表面温度<img file="FDA00002817211000033.GIF" wi="110" he="62" />的计算公式如下:<maths num="0012"><![CDATA[<math><mrow><msubsup><mi>T</mi><mi>casting</mi><mi>t</mi></msubsup><mo>=</mo><mfrac><msubsup><mi>q</mi><mn>0</mn><mi>t</mi></msubsup><mi>h</mi></mfrac><mo>+</mo><msubsup><mi>T</mi><mn>0</mn><mi>t</mi></msubsup></mrow></math>]]></maths>上式中,<img file="FDA00002817211000035.GIF" wi="42" he="58" />为时刻t通过铸件与砂型的界面(1)热通量;<img file="FDA00002817211000036.GIF" wi="45" he="58" />为时刻t拟合得到的砂型内表面温度;h为铸件与砂型的界面(1)换热系数;计算不同时刻的铸件表面温度后,得到铸件表面温度与时间的关系曲线。
地址 100084 北京市海淀区北京市100084-82信箱