張釗楷，杜曉鐘

(太原科技大學(xué)重型機(jī)械教育部工程研究中心，太原 030024)

?

棒線材軋機(jī)計(jì)算機(jī)輔助孔型設(shè)計(jì)工藝軟件開發(fā)

張釗楷，杜曉鐘

(太原科技大學(xué)重型機(jī)械教育部工程研究中心，太原 030024)

計(jì)算機(jī)輔助設(shè)計(jì)發(fā)展以前，棒線材軋輥孔型形狀參數(shù)以及其相關(guān)力能參數(shù)計(jì)算非常耗時(shí)而且復(fù)雜。本軟件基于“兩圓夾一扁”的方法進(jìn)行孔型形狀參數(shù)的設(shè)計(jì)，通過烏薩托夫斯基方法等計(jì)算寬展，利用試驗(yàn)確定艾科隆德方法和西姆斯方法計(jì)算力能參數(shù)，依托C++語言進(jìn)行程序開發(fā)，其目的在于大幅度提高設(shè)計(jì)結(jié)果的精度，降低計(jì)算時(shí)間，準(zhǔn)確的做出孔型圖。通過對某廠22道次線材軋機(jī)的數(shù)據(jù)進(jìn)行采集和計(jì)算，對本軟件進(jìn)行驗(yàn)證。結(jié)果表明軟件具有良好的可靠性和嚴(yán)密性。

孔型設(shè)計(jì)；計(jì)算機(jī)輔助設(shè)計(jì)；棒線材軋制；工藝軟件

隨著科學(xué)技術(shù)的發(fā)展，計(jì)算機(jī)的廣泛應(yīng)用，工程設(shè)計(jì)由手工完成變?yōu)橛捎?jì)算機(jī)輔助完成，孔型設(shè)計(jì)也必然走同樣的發(fā)展道路[1]。計(jì)算機(jī)輔助孔型設(shè)計(jì)在這些年間的發(fā)展，隨著客戶市場的多元化，也要求軋鋼企業(yè)具有快速的響應(yīng)能力，也就要求計(jì)算機(jī)輔助孔型設(shè)計(jì)軟件具有滿足新產(chǎn)品要求和適應(yīng)實(shí)際生產(chǎn)環(huán)境變化的能力。這樣就需要更全面的孔型設(shè)計(jì)軟件來滿足現(xiàn)代棒線材軋鋼企業(yè)的柔性管理[2]。崔彥洲[3]等人開發(fā)的棒材CARD軟件，利用一維搜索法對棒材軋制孔型系統(tǒng)進(jìn)行優(yōu)化設(shè)計(jì)，賈麗娜[4]等人利用以等軸斷面插非等軸斷面的方法開發(fā)的棒材連續(xù)切分軋制計(jì)算機(jī)輔助孔型設(shè)計(jì)軟件，余萬華[5]等人以孔型設(shè)計(jì)實(shí)際經(jīng)驗(yàn)為基礎(chǔ)開發(fā)的計(jì)算機(jī)輔助孔型設(shè)計(jì)系統(tǒng)，呂立華[6]等人以軋制能耗最小為目標(biāo)開發(fā)的計(jì)算機(jī)輔助螺紋鋼孔型設(shè)計(jì)系統(tǒng)，唐文林[7]等人在簡單斷面型鋼計(jì)算機(jī)輔助優(yōu)化設(shè)計(jì)的研究中也開發(fā)了以軋制能耗最小為優(yōu)化目標(biāo)的系統(tǒng)。

此軟件使用Microsoft Visual Studio 2008作為編譯器。用C++語言編寫算法，大大提高軟件的計(jì)算速度，對軋制壓力的計(jì)算采用精準(zhǔn)的數(shù)學(xué)模型，使預(yù)報(bào)精度大大提高，并且針對不同的工況可以選擇不同的力能參數(shù)模型和孔型參數(shù)計(jì)算模型，增大了本系統(tǒng)的使用范圍。利用VC++專有的MFC類庫組成軟件的框架，其內(nèi)部的許多工具使軟件界面更人性化，方便使用者上手。

1 軟件數(shù)學(xué)模型的建立

1.1 斷面面積數(shù)學(xué)模型

在孔型參數(shù)計(jì)算過程中，孔型斷面面積的計(jì)算是重要的一環(huán)，本軟件通過孔型斷面面積的計(jì)算公式利用迭代法求出孔型的高度和寬度，由于遺傳效應(yīng)[7]的存在，使得斷面面積計(jì)算公式精度的高低對后續(xù)的計(jì)算有很大影響，只有降低每一道次的計(jì)算誤差，才能保證最后成品道次的精度。又由于當(dāng)前有限元法計(jì)算斷面面積耗時(shí)比較長，而且棒線材軋制過程條件復(fù)雜，變形不均勻的特點(diǎn)對有限元法的精度有一定影響，所以，現(xiàn)在多用數(shù)學(xué)解析模型計(jì)算軋件斷面面積。如圖1所示，本軟件的孔型斷面面積計(jì)算精度可以控制到不超過1%的高精度范圍內(nèi)，從而后續(xù)獲得的各個(gè)孔型參數(shù)會更加精確。

圖1 φ16.5 mm線材精軋段斷面面積的測量值與本軟件計(jì)算值比較Fig.1 φ16.5 mm wire rod finishing period of section area of the measured values comparing with the calculated values

孔型斷面面積數(shù)學(xué)模型如下：

箱型孔：

(1)

橢圓：

(2)

圓：

(3)

1.2 寬展模型

本軟件在眾多的軋制過程寬展模型中通過理論計(jì)算與實(shí)際數(shù)值的比較，選擇了在考慮不同工況下最接近實(shí)際數(shù)值的烏薩托夫斯基方法、斯米爾諾夫方法以及巴赫契諾夫方法[8]三種數(shù)學(xué)模型。其中烏薩托夫斯基方法主要考慮的影響軋件相對寬展系數(shù)的主要因素是壓下系數(shù)、軋件斷面形狀系數(shù)和輥徑系數(shù)等。斯米爾諾夫利用總功率最小的變分原理，得到的在簡單斷面孔型中軋制的相對寬展系數(shù)。巴赫奇諾夫公式為絕對寬展公式，這種方法在計(jì)算平輥軋制和箱形孔型軋制時(shí)的相對寬展可以得到接近實(shí)際值的數(shù)據(jù)。

1.2.1 烏薩托夫斯基方法

(4)

(5)

(6)

其中β為相對寬展系數(shù)，Wc為相對寬展指數(shù)，當(dāng)ηc在0.1～0.5之間時(shí)用式(5)，其他情況用式(6)，ηc為壓下系數(shù)的倒數(shù)，ηc=hc/Hc(hc、Hc分別為軋后和軋前軋件的平均高度)，δc為軋件斷面形狀系數(shù)，δc=B/Hc(B為軋前寬度)，εc為輥徑系數(shù)，εc=Hc/Dkc(Dkc為軋輥平均工作直徑)。

1.2.2 斯米爾諾夫方法

(7)

1.2.3 巴赫契諾夫方法

(8)

其中f為軋件與軋輥的接觸摩擦系數(shù)，△h為壓下量，H為來料高度，Rkc為軋輥平均工作半徑。

1.3 軋制溫度模型

△T=△Tf+△Tz+△Td-△Tb

(9)

其中△T為每一道的溫降，△Tf為由輻射引起的溫降，△Tz為由傳導(dǎo)引起的溫降，△Td為由對流引起的溫降，△Tb為由變形熱引起的溫升。

1.4 軋制力模型

此軟件在選擇壓力公式時(shí)，考慮軋制速度和軋制材料的屬性，并根據(jù)一些實(shí)測數(shù)據(jù)，通過分析比較，放棄一些誤差較大的公式，從而選擇了艾科隆德單位壓力公式和西姆斯單位壓力公式作為棒線材軋制壓力計(jì)算公式。

1.4.1 艾科隆德單位壓力公式

P=(1+m)(k+hm)

(10)

(11)

K=9.8(14-0.01t)(1.4+C+Mn)

(12)

h=0.1(14-0.01t)

(13)

(14)

其中△hc為軋輥平均壓下量，t為軋制溫度，υ為軋制速度，C、Mn為碳錳的百分含量。一般普碳鋼選用艾科隆德公式計(jì)算可得到準(zhǔn)確的計(jì)算結(jié)果，但是該公式不適用于高速軋制和合金鋼軋制。

1.4.2 西姆斯單位壓力公式

(15)

K=1.15σs

(16)

ε=△Hc/△hc

(17)

其中σs為變形抗力，在艾科隆德公式不適用時(shí)可按鋼種不同選用西姆斯公式。由于不同鋼種計(jì)算所得到的變形抗力不同，所以本軟件根據(jù)實(shí)際數(shù)據(jù)在軟件中包含了一些變形抗力的數(shù)據(jù)，可以作為備選方案，若設(shè)計(jì)人員所選鋼種不在這些備選方案內(nèi)，可通過輸入對應(yīng)鋼種的參數(shù)通過軟件自動計(jì)算變形抗力，如圖2為用兩種方法對φ16 mm普碳鋼計(jì)算后，得到軋制力分布折線圖。

圖2 φ16 mm普碳鋼22道次軋制力分布折線圖Fig.2 Broken line graph of 22 times rolling force distribution of φ16 mm Carbon-steel

2 軟件系統(tǒng)結(jié)構(gòu)

2.1 數(shù)據(jù)的輸入

原始數(shù)據(jù)模塊包括坯料輸入模塊，成品輸入模塊以及各道次參數(shù)輸入模塊三個(gè)部分[9]。如圖3所示為本軟件各個(gè)主要模塊關(guān)系示意圖，用戶根據(jù)提示輸入運(yùn)算參數(shù)及目標(biāo)函數(shù)以后即可進(jìn)行后續(xù)的一系列計(jì)算和輸出。

坯料輸入模塊中輸入坯料的形狀和尺寸，軋件鋼種的化學(xué)成分等。成品輸入模塊中，成品形狀的

圖3 軟件系統(tǒng)示意圖Fig.3 Schematic diagram of software system

選擇采用下拉菜單，選擇完成品形狀后輸入成品直徑以及偏差作為目標(biāo)函數(shù)使用。各道次的參數(shù)輸入模塊中，用戶根據(jù)需要通過下拉菜單選擇孔型系統(tǒng)及軋輥材質(zhì)，并輸入軋輥直徑，輥縫，允許壓力，轉(zhuǎn)速比，軋輥轉(zhuǎn)速等，對有現(xiàn)場實(shí)測溫度能力的軋機(jī)可以直接輸入實(shí)際溫度。

2.2 參數(shù)計(jì)算

圖4 等軸孔設(shè)計(jì)程序流程圖Fig.4 Program flow chart of lsometric hole design

參數(shù)計(jì)算包括軋輥轉(zhuǎn)速的計(jì)算，等軸孔參數(shù)計(jì)算以及扁孔計(jì)算三個(gè)部分。如圖4所示為等軸孔設(shè)計(jì)流程圖，以秒流量相等為限制條件，從而保證無張力或微張力軋制，用戶可以根據(jù)工況確定堆拉系數(shù)，軟件默認(rèn)堆拉系數(shù)為1，通過迭代計(jì)算滿足精度的孔型尺寸，圖5是本軟件的扁孔設(shè)計(jì)界面。

圖5 扁孔設(shè)計(jì)界面Fig.5 Interface of flat hole design

在軋輥轉(zhuǎn)速未知且傳動方式為組傳動的情況下在軋輥轉(zhuǎn)速計(jì)算模塊內(nèi)輸入想要計(jì)算的開始和終了道次數(shù)，軋制方向以及作為目標(biāo)函數(shù)的開始和終了道次的轉(zhuǎn)速來計(jì)算各道次軋輥轉(zhuǎn)速，下一步，進(jìn)入等軸孔計(jì)算模塊。首先，進(jìn)行成品孔計(jì)算，算出成品孔尺寸，其次，輸入需要計(jì)算的等軸孔起始和終了道次數(shù)，修改默認(rèn)的堆拉系數(shù)，最后，點(diǎn)擊等軸孔設(shè)計(jì)按鈕完成等軸孔的參數(shù)計(jì)算，各道次參數(shù)會顯示在對話框上的列表框內(nèi)。進(jìn)入扁孔設(shè)計(jì)之前，需要選擇各道次孔型的充滿度，孔型充滿度確定以后進(jìn)行扁孔設(shè)計(jì)，扁孔參數(shù)式根據(jù)兩圓夾一扁的原理完成設(shè)計(jì)的，在完成各道次的等軸孔設(shè)計(jì)以后，只需點(diǎn)擊計(jì)算按鈕就可以逐一顯示各個(gè)道次的多個(gè)符合條件的方案數(shù)，即可完成整個(gè)孔型設(shè)計(jì)的過程。

2.3 孔型圖輸出

孔型的各道次參數(shù)計(jì)算完成以后軟件可以根據(jù)所選擇的孔型調(diào)用相對應(yīng)的圖形代碼段，并根據(jù)計(jì)算所得的參數(shù)選擇對應(yīng)的圖紙及標(biāo)注尺寸，通過參數(shù)化設(shè)計(jì)的方法輸出孔型圖[10]。如圖6孔型圖都是直線和圓弧組成，每個(gè)道次的孔型都可以表示成直線和圓弧的有序排列。

圖6所示箱形孔型參數(shù)化繪圖示意圖，箱形孔型是對稱結(jié)構(gòu)，所以可以畫出其1/4的結(jié)構(gòu)以后通過鏡像獲得其整個(gè)圖形，圖1、3、5段為直線，2、4段為圓弧，通過計(jì)算獲得的孔型尺寸，在程序中將5段圖形代碼的首位坐標(biāo)值依次相接以后，按中心線鏡像即可獲得一個(gè)完整的箱形孔型圖。

圖6 箱形孔型參數(shù)化繪圖示意圖Fig.6 Schematic diagram of parametric drawing of box pass

2.4 軋制參數(shù)預(yù)報(bào)

此模塊包括溫度的預(yù)報(bào)和軋制力的預(yù)報(bào)。計(jì)算各道軋件軋制溫度，用以設(shè)定或校正初始輸入的軋制溫度。因?yàn)闇亟颠^程非常復(fù)雜，此系統(tǒng)充分考慮了軋制過程中，由熱輻射[11]，熱傳導(dǎo)，熱對流引起的溫降和變形熱引起的溫升對軋制溫度的影響。選擇了比較公認(rèn)的、考慮比較全面的溫降計(jì)算公式，計(jì)算各道次的軋制溫度。將通過計(jì)算得到的溫度稱為計(jì)算溫度，將通過鍵盤輸入的溫度成為設(shè)定溫度?？捎糜?jì)算溫度代替原設(shè)定溫度，也可不用此計(jì)算溫度，只做參考，設(shè)計(jì)人員可以靈活掌握。

圖7 參數(shù)輸入界面Fig.7 Interface of parameter input

3 計(jì)算實(shí)例

本文根據(jù)某廠數(shù)據(jù)輸入計(jì)算參數(shù)，如上圖7所示，根據(jù)軋機(jī)組的具體情況，計(jì)算得到如表1所示列表的右側(cè)數(shù)據(jù)為22道次軋制過程中的精軋段孔型設(shè)計(jì)的計(jì)算方案。左側(cè)數(shù)據(jù)為此22道次生產(chǎn)線的實(shí)際測量數(shù)據(jù)，對比可以發(fā)現(xiàn)本軟件的計(jì)算值和生產(chǎn)線的實(shí)測值非常接近，最大誤差可以控制在2%，平均誤差在1%以內(nèi)。

表1 φ16 mm線材精軋段孔型參數(shù)Tab.1 φ16 mm wire rod finishing period of groove parameters

此22道次軋制線采用4-6-6-6的軋機(jī)布置方式，其中粗軋機(jī)由4臺Φ840×750平-立交替布置，中軋機(jī)由6臺Φ730×750平-立交替布置，預(yù)精軋機(jī)由6臺Φ510×750平-立交替布置，精軋機(jī)由6臺Φ360×650平-立交替布置。如圖7所示，主要工藝參數(shù)包括來料尺寸，來料長度，總軋制道次數(shù)，成品尺寸及偏差，軋輥直徑，輥縫，軋輥轉(zhuǎn)速等。

4 結(jié)論

此套棒線材軋機(jī)計(jì)算機(jī)輔助孔型設(shè)計(jì)工藝軟件計(jì)算精度高，使用方便，充分利用了計(jì)算機(jī)的高速計(jì)算和準(zhǔn)確的推理判斷功能，使程序設(shè)計(jì)具有良好的可靠性，提高了孔型設(shè)計(jì)的效率和首次試軋的成功率，降低生產(chǎn)成本。同時(shí)，自動繪圖功能打破了傳統(tǒng)的人工軟件繪圖的局限性，符合生產(chǎn)實(shí)際要求。

[1] 趙松筠，唐文林，趙靜.棒線材軋機(jī)計(jì)算機(jī)輔助孔型設(shè)計(jì)[M].北京：冶金工業(yè)出版社，2011.

[2] 張?jiān)撇?面向敏捷制造的供應(yīng)鏈柔性管理[D].成都：西南交通大學(xué)，2004.

[3] 崔彥洲.棒材計(jì)算機(jī)輔助孔型設(shè)計(jì)及優(yōu)化[D].河北秦皇島：燕山大學(xué)，2005.

[4] 賈麗娜.棒材連續(xù)切分軋制計(jì)算機(jī)輔助孔型設(shè)計(jì)[J].鋼鐵，2006，41(3)：59-62.

[5] 余萬華，劉飛，高山.計(jì)算機(jī)輔助孔型設(shè)計(jì)系統(tǒng)開發(fā)[J].金屬制品，2014，40(1)：50-58.

[6] 呂立華，馬靳江，唐輝.計(jì)算機(jī)輔助螺紋鋼孔型設(shè)計(jì)系統(tǒng)的開發(fā)及應(yīng)用[J].重慶大學(xué)學(xué)報(bào)：自然科學(xué)版，2005，28(10)：27-30.

[7] 唐文林，顧卓，韓靜濤.簡單斷面型鋼計(jì)算機(jī)輔助優(yōu)化孔型設(shè)計(jì)的研究[J].西安冶金建筑學(xué)院學(xué)報(bào)，1990，22(3)：215-220.

[8] ESTIVALET，MARIE-CLAIRE，MALIMPENSA.Computer program for bar and rod pass design[J].Revue de Metallurgie.Cahiers D'Informations Techniques，1988，85(10)：795-808.

[9] 周志敏，吳迪，白光潤.棒線材孔型計(jì)算機(jī)優(yōu)化設(shè)計(jì)軟件[J].鋼鐵，1996，31(10)：26-29.

[10] 賈柱國，范銀平.計(jì)算機(jī)輔助孔型設(shè)計(jì)在型線材軋機(jī)上的應(yīng)用[J].軋鋼，2001，18(2)：17-19.

[11] 潘璐，帥美榮，黃慶學(xué).孔型設(shè)計(jì)對棒材連軋穩(wěn)定性的影響[J].太原科技大學(xué)學(xué)報(bào)，2012，33(2)：123-127.

Computer-aided Roll Pass Design for Wire Rod and Bar Rolling Mill Process of Software Development

ZHANG Zhao-kai，DU Xiao-zhong

(University of Science and TechnologyTaiyuan，The heavy mechanical engineering research center of Ministry of Education，Taiyuan 030024，China)

Taking into account the development of computer-aided design in the past，it is very time-consuming and complex in terms of calculation about shape parameters of rod and wire roll and its associated force parameters.Based on the approach of“two-round clip a flat”，the software has a design of shape parameters calculates spread by Rostov Chomsky method and force parameters with the Echo Lund and Sims by using experimental methods，and it has a program development relying on C++ language which greatly improves accuracy of design，and reduces the computation time with making an accurate pass figure.By verificating the software by data collection and calculation with the 22 passes of the wire mill factory，the results show that the software has good reliability and rigor.

pass design，computer-aided design，rod and wire rolling，process software

2015-03-09

973計(jì)劃前期研究專項(xiàng)(2012CB722801)

張釗楷(1989-)，男，研究生，主要研究方向?yàn)榘艟€材軋。

1673-2057(2015)06-0465-06

TG332.2

A

10.3969/j.issn.1673-2057.2015.06.011