发明名称 一种板坯连铸机基础辊缝工艺的制定方法
摘要 一种板坯连铸机基础辊缝工艺的制定方法,属于炼钢连铸技术领域。包括钢种热物性参数测量与计算,通过试验手段和运用商业软件JmatPro计算钢种在连铸工艺的冷速条件下,从液相线至室温范围内的瞬时热膨胀系数;凝固传热模型计算,二次开发子程序plotv.f中铸坯的收缩量计算;进行有限元分析计算,获得铸坯处于连铸机内不同位置时的铸坯厚度值,并换算为板坯连铸机各扇形段内凝固收缩量或扇形段锥度。优点在于,改善铸坯中心偏析或,消除中心疏松缺陷,方法简单,经济高效。
申请公布号 CN106001478A 申请公布日期 2016.10.12
申请号 CN201610538966.X 申请日期 2016.07.08
申请人 首钢总公司 发明人 刘珂;孙齐松;王立峰;刘洋;丁宁;徐士新;吕迺冰;佟倩;陈涛;王畅
分类号 B22D11/16(2006.01)I 主分类号 B22D11/16(2006.01)I
代理机构 北京华谊知识产权代理有限公司 11207 代理人 王普玉
主权项 一种板坯连铸机基础辊缝工艺的制定方法,其特征在于,具体步骤及参数如下:1)钢种热物性参数测量与计算:钢种热物性参数主要涉及瞬时热膨胀系数,需要测量、计算两部分工作;钢种在室温至1350℃范围内的瞬时热膨胀系数通过Gleeble热模拟试验机测试得到;但是1350℃至液相范围内的瞬时热膨胀系数受到试验设备能力的限制需要通过商业软件JmatPro根据钢种在凝固过程中的相变规律及各相密度随温度变化规律进行计算;2)凝固传热模型计算:建立适合特定板坯铸机的有限元分析模型,具体模型及热收缩计算原理如下:①根据铸坯几何及传热的对称性,取铸坯的四分之一断面作为凝固传热分析模型的计算域;忽略铸坯拉坯方向的热量损失,采用2D切片法计算铸坯在连铸机内的热状态,模型采用均匀网格划分;设板坯宽为A,厚为B,将1/2边长分别取m节点和n节点,取X,Y方向空间步长为Δx和Δy,于是有:<maths num="0001"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>x</mi><mo>=</mo><mfrac><mi>A</mi><mrow><mn>2</mn><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow>]]></math><img file="FDA0001045157150000011.GIF" wi="690" he="69" /></maths><maths num="0002"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>y</mi><mo>=</mo><mfrac><mi>B</mi><mrow><mn>2</mn><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow>]]></math><img file="FDA0001045157150000012.GIF" wi="664" he="71" /></maths>在垂直X‑Y平面的时间轴上,取时间步长为Δt,则Δt=t<sub>p+1</sub>‑t<sub>p</sub>=t<sub>p</sub>‑t<sub>p‑1</sub>;因此,研究区域<img file="FDA0001045157150000013.GIF" wi="641" he="63" />就转变为一系列离散点;记<img file="FDA0001045157150000014.GIF" wi="66" he="59" />为t<sub>p</sub>时刻结点(i,j)位置上点的温度;②二次开发子程序plotv.f中铸坯的收缩量计算方法如下:a.单位厚度的切片以一定的时间步长从结晶器弯月面处开始向下移动,依次通过结晶器冷却区、二次冷却区和空冷辐射区,直到出计算域为止;在每一时间点上计算该切片的温度分布;b.凝固收缩计算从弯月面初开始,因为粘滞性温度以上热膨胀系数被强制设定为零,这样就保证了从初始凝固开始具有收缩现象;c.在凝固计算过程中,每一个给定是时间点上的结点温度的改变是相对于该结点前一个时间点上的温度;每个结点的长度的变化是由该结点温度相对于前一时间的温度的改变而产生的;采用向前差分离散,宽度方向上的节点热应变为:<maths num="0003"><math><![CDATA[<mrow><msup><mi>&Delta;&epsiv;</mi><mrow><mi>t</mi><mi>h</mi></mrow></msup><mo>=</mo><mi>&alpha;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>&Delta;x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mrow><mrow><msubsup><mi>&Delta;x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup></mrow></mfrac><mo>-</mo><mn>1</mn></mrow>]]></math><img file="FDA0001045157150000021.GIF" wi="677" he="94" /></maths>其中,Δε<sup>th</sup>为(i,j)节点在p+1时刻的热收缩应变;<img file="FDA0001045157150000022.GIF" wi="203" he="63" />分别为(i,j)节点在p,p+1时刻的温度;α为<img file="FDA0001045157150000023.GIF" wi="184" he="67" />温度范围内的调和平均热膨胀系数;<img file="FDA0001045157150000024.GIF" wi="254" he="66" />为(i,j)节点在p、p+1时刻的宽度方向尺寸;该时间点p+1上宽度方向上结点尺寸为:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>&Delta;x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>&Delta;x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>+</mo><mi>&alpha;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mi>&Delta;x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup></mrow>]]></math><img file="FDA0001045157150000025.GIF" wi="745" he="67" /></maths>同理,厚度方向上的上结点尺寸为:<maths num="0005"><math><![CDATA[<mrow><msubsup><mi>&Delta;y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>&Delta;y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>+</mo><mi>&alpha;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mi>&Delta;y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup></mrow>]]></math><img file="FDA0001045157150000026.GIF" wi="738" he="63" /></maths>其中:<img file="FDA0001045157150000027.GIF" wi="382" he="110" /><img file="FDA0001045157150000028.GIF" wi="303" he="68" />分别为p、p+1时刻的瞬时热膨胀系数;d.每一个给定的时间点上铸坯在尺寸上由该方向上所有结点收缩叠加得到;时间点p+1上宽度方向上自然收缩量:<maths num="0006"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>A</mi></msub><mo>=</mo><msup><mi>A</mi><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>-</mo><msup><mi>A</mi><mn>0</mn></msup><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mrow><mo>(</mo><mi>p</mi><mo>+</mo><mn>1</mn><mo>)</mo><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi></mrow></msubsup><mrow><mo>(</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></msubsup><mi>&alpha;</mi><mo>&CenterDot;</mo><mo>(</mo><mrow><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup></mrow><mo>)</mo><mo>&CenterDot;</mo><msubsup><mi>&Delta;x</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001045157150000029.GIF" wi="1133" he="81" /></maths>其中,δ<sub>A</sub>为铸坯在宽度方向尺寸变形量;A<sup>p+1</sup>为铸坯宽度在p+1时刻尺寸;A<sup>0</sup>为铸坯宽度在初始状态下尺寸;时间点p+1上厚度方向上自然收缩量:<maths num="0007"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>B</mi></msub><mo>=</mo><msup><mi>B</mi><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>-</mo><msup><mi>B</mi><mn>0</mn></msup><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mrow><mo>(</mo><mi>p</mi><mo>+</mo><mn>1</mn><mo>)</mo><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>t</mi></mrow></msubsup><mrow><mo>(</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mi>&alpha;</mi><mo>&CenterDot;</mo><mo>(</mo><mrow><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>p</mi></msubsup></mrow><mo>)</mo><mo>&CenterDot;</mo><msubsup><mi>&Delta;y</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010451571500000210.GIF" wi="1143" he="86" /></maths>其中,δ<sub>B</sub>为铸坯在宽度方向尺寸变形量;B<sup>p+1</sup>为铸坯宽度在p+1时刻尺寸;B<sup>0</sup>为铸坯宽度在初始状态下尺寸;e.重复以上步骤a~d直到计算结束为止;3)进行有限元分析计算,获得铸坯处于连铸机内不同位置时的铸坯厚度值,并换算为板坯连铸机各扇形段内凝固收缩量或扇形段锥度。
地址 100041 北京市石景山区石景山路68号