发明名称 一种确定直缝焊管焊接温度的方法
摘要 本发明属于一种确定直缝焊管焊接温度的方法。其技术方案是:在焊接挤压辊和矫直辊之间的焊管表面上选取4~16个测温点进行实测;运用有限元法计算不同的传热边界及焊接热源密度条件下测温的温度值;建立温度计算值和实测值之差平方和与传热边界条件及焊接热源密度间的函数关系模型;以温度计算值和实测值之差平方和最小为目标,优化确定焊接时空气传热、接触传热、冷却水传热的传热边界条件以及焊接热源密度的实际热交换系数,据此确定直缝焊管焊接温度、焊缝附近温度及冷却速度场。本发明所确定的直缝焊管焊接温度、焊缝附近温度场和焊缝附近冷却速度场接近实际,能预测焊缝附近的微观组织和为生产过程的规范化管理提供基础数据。
申请公布号 CN103400037A 申请公布日期 2013.11.20
申请号 CN201310323133.8 申请日期 2013.07.29
申请人 武汉科技大学 发明人 李立新;刘升;李烨;朱少文;胡盛德;曾祥明;周绪昌
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 张火春
主权项 1.一种确定直缝焊管焊接温度的方法,其特征在于该方法的具体步骤是:第一步、温度实测在焊接挤压辊和矫直辊之间的焊管表面上选取4~16个测温点,实测所述测温点处直缝焊管表面的温度值;第二步、虚拟试验方案的设计直缝焊管对空气传热的待定热交换系数:<maths num="0001"><![CDATA[<math><mrow><msub><mi>h</mi><mi>r</mi></msub><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><mi>&sigma;&epsiv;</mi><mrow><mo>(</mo><msup><mi>T</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>T</mi><mi>f</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><msub><mi>T</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(1)中:x<sub>1</sub>为直缝焊管对空气传热的待定修正系数;σ为波尔兹曼常量,W/(m<sup>2</sup>·K<sup>4</sup>);ε为辐射率;T<sub>f</sub>为环境绝对温度,K;T为焊管表面的绝对温度,K;直缝焊管与轧辊接触传热的待定热交换系数:h<sub>c</sub>=2x<sub>2</sub>(kρct/π)<sup>1/2</sup>   (2)式(2)中:x<sub>2</sub>为直缝焊管与轧辊接触传热的待定修正系数;t为接触时间,s;ρ为材料密度,kg/m<sup>3</sup>;c为比热,J/(kg·℃);直缝焊管与冷却水传热的待定热交换系数:h<sub>w</sub>=125x<sub>3</sub>W<sup>0.663</sup>×10<sup>-0.00147T</sup>   (3)式(3)中:x<sub>3</sub>为直缝焊管与冷却水传热的待定修正系数;W为冷却水流量,L/s;待定的焊接热源密度:<maths num="0002"><![CDATA[<math><mrow><mi>q</mi><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mn>4</mn></msub><mi>P</mi><mo>&times;</mo><mn>86.5</mn><mo>%</mo></mrow><mi>V</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(4)中:x<sub>4</sub>为焊接热源密度的待定修正系数;P为焊接功率,W;V为涡流透入层体积,m<sup>3</sup>;以式(1)~式(4)中对应的直缝焊管对空气传热的待定修正系数x<sub>1</sub>、直缝焊管与轧辊接触传热的待定修正系数x<sub>2</sub>、直缝焊管与冷却水传热的待定修正系数x<sub>3</sub>和焊接热源密度的待定修正系数x<sub>4</sub>作为4个试验因素;先选取正交表,再以0.55和1.35为边界值,在0.55~1.35间均匀确定各试验因素的水平值,然后按正交表进行虚拟试验方案的设计;第三步、传热有限元模拟计算直缝焊管温度场的控制方程是:<maths num="0003"><![CDATA[<math><mrow><mi>&rho;c</mi><mfrac><mrow><mo>&PartialD;</mo><mi>T</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mi>&lambda;</mi><mrow><mo>(</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>T</mi></mrow><mrow><mo>&PartialD;</mo><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>T</mi></mrow><msup><mrow><mo>&PartialD;</mo><mi>y</mi></mrow><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><mrow><msup><mo>&PartialD;</mo><mn>2</mn></msup><mi>T</mi></mrow><msup><mrow><mo>&PartialD;</mo><mi>z</mi></mrow><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>+</mo><mi>q</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(5)中:ρ为材料密度,kg/m<sup>3</sup>;c为材料比热,J/(kg·℃);t为时间,s;λ为材料热导率,W/(m·℃);q=q(x,y,z,t)为材料内部热源密度,W/m<sup>3</sup>;确定直缝焊管焊接温度、焊缝附近温度场和焊缝附近冷却速度场所需的边界条件由式(1)~式(4)确定;则传热有限元模型由式(1)~式(5)组成;利用所述传热有限元模型,按正交表确定的虚拟试验方案进行传热有限元模拟计算,得到每个试验方案各测温点处的温度计算值。第四步、误差平方和模型的建立根据不同试验方案各测温点的温度计算值,结合各测温点的温度实测值,得到各测温点处的温度计算值与实测值之间的误差平方和,建立各测温点的温度计算值与实测值之间的误差平方和模型:<maths num="0004"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><msub><mi>a</mi><mi>k</mi></msub><msubsup><mi>x</mi><mi>i</mi><mi>m</mi></msubsup><msubsup><mi>x</mi><mi>j</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2,3,4</mn><mo>;</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>=</mo><mn>0,1,2</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(6)中:a<sub>k</sub>(k=1,2,…,41)为回归系数;x<sub>i</sub>(i=1,2,3,4)、x<sub>j</sub>(j=1,2,3,4)分别为直缝焊管对空气传热的修正系数x<sub>1</sub>、直缝焊管与轧辊接触传热的修正系数x<sub>2</sub>、直缝焊管与冷却水传热的修正系数x<sub>3</sub>和焊接热源密度的修正系数x<sub>4</sub>中的任一值;第五步、修正系数的优化以各测温点温度实测值与对应计算值之间误差平方和最小为追求目标,运用罚函数法,分别对直缝焊管对空气传热的待定修正系数x<sub>1</sub>、直缝焊管与轧辊接触传热的待定修正系数x<sub>2</sub>、直缝焊管与冷却水传热的待定修正系数x<sub>3</sub>和焊接热源密度的待定修正系数x<sub>4</sub>进行优化,得到直缝焊管对空气传热的实际修正系数<img file="FDA00003585161800032.GIF" wi="85" he="76" />直缝焊管与轧辊接触传热的实际修正系数<img file="FDA00003585161800033.GIF" wi="84" he="79" />直缝焊管与冷却水传热的实际修正系数<img file="FDA00003585161800034.GIF" wi="52" he="76" />和焊接热源密度的实际修正系数<img file="FDA00003585161800035.GIF" wi="80" he="77" />第六步、直缝焊管焊接温度的确定用直缝焊管对空气传热的实际修正系数<img file="FDA00003585161800036.GIF" wi="46" he="77" />代替式(1)中的直缝焊管对空气传热的待定修正系数x<sub>1</sub>,得到直缝焊管对空气传热的实际热交换系数:<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>h</mi><mi>r</mi><mo>&prime;</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mn>1</mn><mo>&prime;</mo></msubsup><mi>&sigma;&epsiv;</mi><mrow><mo>(</mo><msup><mi>T</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>T</mi><mi>f</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><msub><mi>T</mi><mi>f</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(7)中:<img file="FDA00003585161800038.GIF" wi="49" he="78" />为直缝焊管对空气传热的实际修正系数;σ为波尔兹曼常量,W/(m<sup>2</sup>·K<sup>4</sup>);ε为辐射率;T<sub>f</sub>为环境绝对温度,K;T为焊管表面的绝对温度,K;用直缝焊管与轧辊接触传热的实际修正系数<img file="FDA00003585161800041.GIF" wi="51" he="74" />代替式(2)中的直缝焊管与轧辊接触传热的待定修正系数x<sub>2</sub>,得到直缝焊管与轧辊接触传热的实际热交换系数:<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>h</mi><mi>c</mi><mo>&prime;</mo></msubsup><mo>=</mo><mn>2</mn><msubsup><mi>x</mi><mn>2</mn><mo>&prime;</mo></msubsup><msup><mrow><mo>(</mo><mi>k &rho;ct</mi><mo>/</mo><mi>&pi;</mi><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(8)中:<img file="FDA00003585161800043.GIF" wi="56" he="74" />为直缝焊管与轧辊接触传热的实际修正系数;t为接触时间,s;ρ为材料密度,kg/m<sup>3</sup>;c为比热,J/(kg·℃)。用直缝焊管与冷却水传热的实际修正系数<img file="FDA00003585161800044.GIF" wi="52" he="78" />替代式(3)中的直缝焊管与冷却水传热的待定修正系数x<sub>3</sub>,得到直缝焊管与冷却水传热的实际热交换系数:<maths num="0007"><![CDATA[<math><mrow><msub><mi>h</mi><mi>w</mi></msub><mo>=</mo><mn>125</mn><msubsup><mi>x</mi><mn>3</mn><mo>&prime;</mo></msubsup><msup><mi>W</mi><mn>0.663</mn></msup><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>0.00147</mn><mi>T</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(9)中:<img file="FDA00003585161800046.GIF" wi="52" he="76" />为直缝焊管与冷却水传热的实际修正系数;W为冷却水流量,L/s;用焊接热源密度的实际修正系数<img file="FDA00003585161800047.GIF" wi="52" he="77" />代替式(4)中的待定的焊接热源密度的待定修正系数x<sub>4</sub>,得到实际的焊接热源密度:<maths num="0008"><![CDATA[<math><mrow><mi>q</mi><mo>=</mo><mfrac><mrow><msubsup><mi>x</mi><mn>4</mn><mo>&prime;</mo></msubsup><mi>P</mi><mo>&times;</mo><mn>86.5</mn><mo>%</mo></mrow><mi>V</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(10)中:<img file="FDA00003585161800049.GIF" wi="56" he="78" />为焊接热源密度的实际修正系数;P为焊接功率,W;V为涡流透入层体积,m<sup>3</sup>;以式(7)~式(10)的计算结果为边界条件,运用传热有限元法,确定直缝焊管的焊接温度、焊缝附近的温度和冷却速度场。
地址 430081 湖北省武汉市青山区建设一路