專(zhuān)利名稱(chēng)：平整機(jī)成品板面粗糙度預(yù)報(bào)與控制技術(shù)的制作方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種平整軋制工藝技術(shù)，特別涉及一種平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)。
背景技術(shù)：
表面粗糙度作為平整成品帶鋼重要的特性之一，它不僅影響帶鋼沖壓時(shí)的變形行為和涂鍍后的外觀面貌，而且可以改變材料的耐蝕性。從目前的質(zhì)量發(fā)展趨勢(shì)來(lái)看，帶鋼表面的粗糙度控制已經(jīng)成為重要的質(zhì)量指標(biāo)，尤其對(duì)汽車(chē)板至關(guān)重要。如何建立平整機(jī)機(jī)板面粗糙度預(yù)報(bào)模型，對(duì)要生產(chǎn)的帶鋼表面粗糙度值進(jìn)行預(yù)報(bào)，最終達(dá)到事先調(diào)整，避免粗糙度超差已經(jīng)成為現(xiàn)場(chǎng)技術(shù)攻關(guān)的重點(diǎn)與難點(diǎn)。對(duì)于這方面的研究，國(guó)外已經(jīng)開(kāi)始起步，并作了一些相關(guān)研究工作，可見(jiàn)的文獻(xiàn)有劍持一仁的“US430鋼板冷軋過(guò)程中表面光潔度影響因素的研究”、“軋制條件對(duì)SUS430鋼板表面粗糙度的影響”[“剣持一仁的“US430鋼板の冷間壓延における表面光沢おょびヒ一トストリ一クに及ぼすロ一ルコ一テイングの影響”(Journal of the JSTP，1999，40(457)129-133)、“SUS430鋼板の表面光沢に及ぼす調(diào)質(zhì)壓延條件の影響”(Journal of the JSTP，2000，41(468)54-58)]，但尚未檢索到有任何國(guó)外文獻(xiàn)提到關(guān)于平整機(jī)機(jī)板面粗糙度預(yù)報(bào)與控制模型方面的內(nèi)容。國(guó)內(nèi)對(duì)于板面粗糙度方面的研究還停留在對(duì)冷連軋機(jī)進(jìn)行定性研究方面，而且偏重于其它方面如板形的研究，主要文獻(xiàn)見(jiàn)張清東等人的“冷軋帶鋼板形和表面粗糙度協(xié)同控制”(北京科技大學(xué)學(xué)報(bào)，2005，27(2)232-234)。而國(guó)內(nèi)對(duì)于平整機(jī)成品帶鋼板面粗糙度預(yù)報(bào)模型方面的研究屬于空白，沒(méi)有發(fā)現(xiàn)任何相關(guān)的文獻(xiàn)。

發(fā)明內(nèi)容
本發(fā)明的目的是提供一種適合平整軋制工藝的板面粗糙度在線(xiàn)預(yù)報(bào)與控制模型，使現(xiàn)場(chǎng)能夠?qū)崿F(xiàn)根據(jù)用戶(hù)對(duì)帶材表面粗糙度的要求適時(shí)調(diào)整軋制工藝參數(shù)以控制成品帶鋼的表面粗糙度。為了實(shí)現(xiàn)上述目的，本發(fā)明采用了以下技術(shù)方案這種平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，包括以下步驟a、建立反映平整機(jī)工作輥粗糙度與工作輥原始粗糙度、換輥后軋制公里數(shù)之間關(guān)系的數(shù)學(xué)模型；按下述步驟建立反映平整機(jī)工作輥粗糙度與工作輥原始粗糙度、換輥后軋制公里數(shù)之間關(guān)系的數(shù)學(xué)模型
1)以Rar=Rar0·e-BL·L]]>的形式構(gòu)造出平整機(jī)工作輥粗糙度與軋制公里數(shù)之間關(guān)系的基本模型框架；式中Rar-平整機(jī)工作輥粗糙度；Rar0-平整機(jī)工作輥原始粗糙度；L-工作輥換輥后的軋制公里數(shù)；BL.-工作輥粗糙度衰減系數(shù)；2)在現(xiàn)場(chǎng)收集一組關(guān)于不同軋制公里數(shù)Li下所對(duì)應(yīng)的工作輥粗糙度Rari以及相應(yīng)的工作輥原始粗糙度Rar0i的數(shù)據(jù)；3)給定工作輥粗糙度衰減系數(shù)初始設(shè)定值X0＝{BL′}與計(jì)算精度∑＝0.00001；4)將BL′、Li、Rar0i代入公式Rar=Rar0·e-BL·L,]]>求出相應(yīng)的工作輥粗糙度理論計(jì)算值Rari′；5)計(jì)算出相應(yīng)的目標(biāo)函數(shù)值F(X)；6)判斷F(X)≤∑是否成立，若不成立改變初始值X0重復(fù)上述步驟4)、5)，直到判斷條件成立，結(jié)束計(jì)算，得出特定平整機(jī)的工作輥粗糙度衰減系數(shù)BL以及相應(yīng)的反映平整機(jī)工作輥粗糙度與工作輥原始粗糙度、換輥后軋制公里數(shù)之間關(guān)系的數(shù)學(xué)模型；b、構(gòu)造壓印率Krs與遺傳率Kss的函數(shù)表達(dá)式；c、構(gòu)造平整后成品板面粗糙度Rastrip的初始模型；d、根據(jù)特定平整機(jī)的現(xiàn)場(chǎng)實(shí)驗(yàn)數(shù)據(jù)，建立特定平整機(jī)的板面粗糙度預(yù)報(bào)模型，按下述步驟建立特定平整機(jī)的板面粗糙度模型1)、現(xiàn)場(chǎng)采集一組試驗(yàn)數(shù)據(jù)，滿(mǎn)足以下要求帶材的來(lái)料厚度都為hε，材料強(qiáng)度都為kε，而壓下量共有ε1，ε2，…，εi，…，εn共n個(gè)，而每個(gè)壓下量εi對(duì)應(yīng)兩個(gè)不同的平整前入口帶材表面粗糙度與工作輥軋輥粗糙度的值(關(guān)于這點(diǎn)，可以通過(guò)選擇不同的工作輥軋制公里數(shù)來(lái)實(shí)現(xiàn))，設(shè)其為Rage1i，Rage2i，Rar1i，Rar2i。這樣，與之對(duì)應(yīng)，通過(guò)現(xiàn)場(chǎng)實(shí)際測(cè)量就可以得出Rastrip11，Rastrip21，Rastrip12，…，Rastrip1i，Rastriop2i，…，Rastrip1n，Rastrip2n共計(jì)2n個(gè)帶鋼表面粗糙度的值；2)定義中間函數(shù)G(ε)、F(ε)；3)根據(jù)實(shí)驗(yàn)數(shù)據(jù)求解出相應(yīng)的G(εi)、F(εi)的表達(dá)式；4)根據(jù)G(εi)、F(εi)的表達(dá)式通過(guò)曲線(xiàn)擬合的方法給出G(ε)、F(ε)的具體表達(dá)式；5)采集材料強(qiáng)度都為kh、壓下量都為εh，而帶材來(lái)料厚度分別為h1，h2，…，hi，…，hn的第二組實(shí)驗(yàn)數(shù)據(jù)與帶材來(lái)料厚度都為hk、壓下量都為εk，而材料強(qiáng)度分別為k1，k2，…，ki，…，kn的第三組實(shí)驗(yàn)數(shù)據(jù)；6)定義中間函數(shù)G1(h)、F1(h)、G2(h)、F2(h)；7)根據(jù)實(shí)驗(yàn)數(shù)據(jù)通過(guò)曲線(xiàn)擬合的方法給出G1(h)、F1(h)、G2(h)、F2(h)的具體表達(dá)式；8)根據(jù)G(ε)、F(ε)、G1(h)、F1(h)、G2(h)、F2(h)求出壓印率Krs與遺傳率Kss的具體函數(shù)表達(dá)式；9)輸出特定平整機(jī)的成品板面粗糙度預(yù)報(bào)模型。
e、根據(jù)板面粗糙度模型在線(xiàn)預(yù)報(bào)成品板面粗糙度。
按下述步驟實(shí)現(xiàn)成品板面粗糙度的預(yù)報(bào)1)收集待平整帶鋼的主要特征參數(shù)；2)收集平整機(jī)工作輥主要特征參數(shù)；3)收集主要平整軋制工藝參數(shù)；4)計(jì)算平整機(jī)工作輥粗糙度；5)將平整帶鋼參數(shù)、工作輥特征參數(shù)、平整軋制工藝參數(shù)等代入平整機(jī)成品板面粗糙度模型，計(jì)算相應(yīng)成品板面粗糙度。
f、成品板面粗糙度的在線(xiàn)控制，按下述步驟實(shí)現(xiàn)成品板面粗糙度的在線(xiàn)控制1)收集待平整帶鋼平整前的板面粗糙度Rage、帶鋼強(qiáng)度k等主要參數(shù)；2)收集工作輥換輥后軋制公里數(shù)L、工藝技術(shù)規(guī)程所規(guī)定的工作輥軋制公里數(shù)(即換輥周期)Lmax0等軋輥主要特征參數(shù)；3)收集用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip以及用戶(hù)允許最大延伸率εmax、最小延伸率εmin等需求參數(shù)；4)按照現(xiàn)場(chǎng)經(jīng)驗(yàn)設(shè)定工作輥的原始粗糙度的初始值Rar00；5)引入相關(guān)系數(shù)α、β，并給定初始值α＝1.0、β＝0；6)令Rar0＝Rar00，而軋制公里數(shù)分別取L1＝α·Lmax0、L2＝β·Lmax0(代表最小軋制公里數(shù)，新輥時(shí)最小軋制公里數(shù)為零)代入步驟a)所述模型計(jì)算出相應(yīng)的工作輥實(shí)際可能最小粗糙度Rarmin與最大粗糙度Rarmax；7)將平整前的板面粗糙度Rage、帶鋼強(qiáng)度k、工作輥實(shí)際可能最大粗糙度Rarmax、實(shí)際可能最小粗糙度Rarmin以及用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip代入成品板面粗糙度模型，反算出相應(yīng)的最大延伸率εmax′、最小延伸率εmin′；8)判斷不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax是否同時(shí)成立，如果成立，則繼續(xù)步驟10)；如果不成立，調(diào)整工作輥的原始粗糙度的初始設(shè)定值Rar0，重復(fù)步驟4)、5)以及6)判斷能否找到合適的原始粗糙度設(shè)定值Rar0y，使得上述兩個(gè)不等式同時(shí)成立，如果可以找到則將工作輥原始粗糙度設(shè)定為Rar0y，繼續(xù)步驟10)；如無(wú)法合適的原始粗糙度值，則進(jìn)入步驟9)；9)判斷εmax′≥εmin是否成立，如果成立，則減小工作輥原始粗糙度的設(shè)定值至等式εmin′＝εmin成立，以一定的步長(zhǎng)逐步增大β的值，重復(fù)步驟5)、6)、7)，使得不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax同時(shí)成立，然后進(jìn)入步驟10)；如果不成立，則增大工作輥的原始粗糙度的初始設(shè)定值至等式εmax′＝εmax成立，以一定的步長(zhǎng)逐步減小α的值，重復(fù)步驟5)、6)、7)，使得不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax同時(shí)成立，然后進(jìn)入步驟10)；10)將平整前的板面粗糙度Rage、帶鋼強(qiáng)度k、工作輥實(shí)際粗糙度Rar以及用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip代入成品板面粗糙度模型，反算出延伸率ε，進(jìn)行延伸率設(shè)定。同時(shí)注意，當(dāng)α<1.0]]>時(shí)，則在工作輥軋制公里數(shù)達(dá)到α·Lmax0時(shí)，及時(shí)轉(zhuǎn)換規(guī)格或換輥；而當(dāng)β>0]]>時(shí)，則應(yīng)該先安排生產(chǎn)其它規(guī)格的產(chǎn)品，直到工作輥軋制公里數(shù)已經(jīng)達(dá)到β·Lmax0時(shí)開(kāi)始生產(chǎn)該計(jì)劃規(guī)格產(chǎn)品。
步驟a中所述目標(biāo)函數(shù)定義為F(X)=Σi=1n(Rari-Rari′)2]]>式中n-樣本數(shù)步驟b中所述壓印率Krs與遺傳率Kss的函數(shù)表達(dá)式為Krs=f1(h)·f2(k)·f(ϵ)Kss=g1(h)·g2(k)·g(ϵ)]]>式中h-所平整帶材的來(lái)料厚度k-所平整帶材的強(qiáng)度ε-平整壓下率步驟c中所述平整板面粗糙度的初始模型表達(dá)式如下Rastrip=Kss·Rage+Krs·Rar0·e-BL·L]]>式中Rastrip-平整后成品帶鋼表面粗糙度Rage-平整前來(lái)料帶材表面粗糙度步驟d中所述中間函數(shù)G(ε)、F(ε)的初始表達(dá)式如下
G(ε)＝g1(hε)·g2(Kε)·g(ε)F(ε)＝f1(hε)·f2(kε)·f(ε)步驟d中所述求解出的G(εl)、F(εi)的表達(dá)式如下G(ϵi)=Rar2i·astrip1i-Rar1i·Rastrip2iRar2i·Rage1i-Rar1i·Rage2iF(ϵi)=Rage2i·astrip1i-Rage1i·Rastrip2iRage2i·Rar1i-Rage1i·Rar2i]]>步驟d中所述中間函數(shù)G1(h)、F1(h)、G2(h)、F2(h)的表達(dá)式如下G1(h)=g2(kh)·g(ϵh)·g1(h)F1(h)=f2(kh)·f(ϵh)·f1(h)G2(k)=g1(hk)·g(ϵk)·g2(k)F2(k)=f1(hk)·f(ϵk)·f2(k)]]>步驟d中所述印率Krs與遺傳率Kss的具體函數(shù)表達(dá)式如下Krs=f1(h)·f2(k)·f(ϵ)=F(ϵ)·F1(h)·F2(k)F(ϵh)·F2(kh)]]>Kss=g1(h)·g2(k)·g(ϵ)=G(ϵ)·G1(h)·G2(k)G(ϵh)·G2(kh)]]>步驟d中所述成品板面粗糙度模型如下Rastrip=G(ϵ)·G1(h)·G2(k)G(ϵh)·G2(kh)·Rage+F(ϵ)·F1(h)·F2(k)F(ϵh)·F2(kh)·Rar0·e-BL·L]]>步驟e中所述平整帶鋼的主要特征參數(shù)包括平整前帶鋼表面粗糙度、帶鋼強(qiáng)度。
步驟e中所述平整機(jī)工作輥主要特征參數(shù)包括工作輥原始粗糙度、工作輥換輥后軋制公里數(shù)。
步驟e中所述平整軋制工藝參數(shù)包括平整軋制延伸率。
步驟f中所述α、β平整機(jī)軋制計(jì)劃系數(shù)，當(dāng)α<1.0]]>時(shí)，應(yīng)將特定規(guī)格產(chǎn)品安排在換輥初期軋制；而當(dāng)β>0]]>時(shí)，則應(yīng)該安排換輥后期軋制。


通過(guò)以下結(jié)合附圖對(duì)本發(fā)明較佳實(shí)施例的描述，可以進(jìn)一步理解本發(fā)明的目的、特征和優(yōu)點(diǎn)。
圖1是本發(fā)明中平整機(jī)工作輥粗糙度模型計(jì)算流程圖；圖2是本發(fā)明一個(gè)較佳實(shí)施例的特定平整機(jī)的板面粗糙度模型求解計(jì)算流程圖；
圖3是按照本發(fā)明一個(gè)較佳實(shí)施例的平整機(jī)成品板面粗糙度在線(xiàn)預(yù)報(bào)計(jì)算流程圖；圖4是按照本發(fā)明一個(gè)較佳實(shí)施例的平整機(jī)成品板面粗糙度在線(xiàn)控制流程圖；圖5是本發(fā)明中平整機(jī)成品板面粗糙度預(yù)報(bào)與控制總流程圖具體實(shí)施方式
實(shí)施例1根據(jù)軋制理論與現(xiàn)場(chǎng)經(jīng)驗(yàn)可以知道，對(duì)于平整軋制過(guò)程而言，成品帶鋼的板面粗糙度主要是由工作輥復(fù)印到帶鋼上面的粗糙度與平整前帶鋼粗糙度的遺傳部分疊加而成。其中，工作輥復(fù)印到帶鋼表面的粗糙度與工作輥本身的粗糙度以及壓印率(即軋輥粗糙度傳遞到帶鋼上的能力)有關(guān)；平整前帶鋼遺傳部分的粗糙度則取于平整前帶鋼的原始粗糙度與相應(yīng)的遺傳率。由此可見(jiàn)，要實(shí)現(xiàn)粗糙度的在線(xiàn)預(yù)報(bào)，求出壓印率與遺傳率的表達(dá)式非常關(guān)鍵，本發(fā)明采用下式來(lái)表達(dá)平整后成品帶鋼的表面粗糙度Rastrip＝Kss·Rage+Krs·Rar(28)式中，Rastrip為成品帶鋼表面粗糙度，Rage為平整前帶材表面粗糙度，Rar為平整過(guò)程中工作輥的輥面粗糙度，Kss為平整前帶鋼表面粗糙度的遺傳率，Krs為軋輥粗糙度傳遞到帶鋼上能力即壓印率。顯然，在成品板面粗糙度的預(yù)報(bào)過(guò)程中，平整前帶材的表面粗糙度是可以實(shí)際測(cè)量出來(lái)的，是一個(gè)已知值，而平整過(guò)程中工作輥的表面粗糙度Rar則可以通過(guò)圖1所述工作輥輥面粗糙度計(jì)算模型由工作輥軋制公里數(shù)L以及工作輥原始粗糙度Ra0求解得到，也認(rèn)為是已經(jīng)知道的。所以只要將Kss、Krs的表達(dá)式求解出來(lái)，那么特定平整機(jī)的粗糙度模型就建立起來(lái)了。
圖2為按照本發(fā)明一個(gè)較佳實(shí)施例的特定平整機(jī)成品板面粗糙度模型求解計(jì)算的流程圖。該實(shí)施例被用于特定平整機(jī)成品板面粗糙度模型的建立，并且具有大量的關(guān)于該平整機(jī)在平整前入口帶材表面粗糙度與工作輥軋輥粗糙度的值的實(shí)驗(yàn)數(shù)據(jù)。
在步驟21中，首先構(gòu)造出遺傳率Kss與復(fù)印率Krs的影響因素表達(dá)式，并給出相應(yīng)的成品帶鋼板面粗糙度影響因素模型，分別如下式所示Krs=f1(h)·f2(k)·f(ϵ)Kss=g1(h)·g2(k)·g(ϵ)---(29)]]>Rastrip=g1(h)·g2(k)·g(ϵ)·Rage+f1(h)·f2(k)·f(ϵ)·Rar0·e-BL·L---(30)]]>式(30)中，Rar0·e-BL·L是平整過(guò)程中工作輥輥面粗糙度模型，可由圖1求出，其中Rar0表示工作輥原始粗糙度，L表示軋制公里數(shù)，BL表示平整機(jī)工作輥粗糙度隨著軋制公里數(shù)的衰減率。
隨后在步驟22中，選擇一組試驗(yàn)數(shù)據(jù)，滿(mǎn)足以下要求帶材的來(lái)料厚度都為hε，材料強(qiáng)度都為kε，而壓下量共有ε1，ε2，…，εi，…，εn共n個(gè)，而每個(gè)壓下量εi對(duì)應(yīng)兩個(gè)不同的平整前帶材表面粗糙度與工作輥粗糙度的值(關(guān)于這點(diǎn)，可以通過(guò)選擇不同的工作輥軋制公里數(shù)來(lái)實(shí)現(xiàn))，設(shè)其為Rage1i，Rage2i，Rar1i，Rar2i。這樣，與之對(duì)應(yīng)，通過(guò)現(xiàn)場(chǎng)實(shí)際測(cè)量就可以得出Rastrip11，Rastrip21，Rastrip12，Rastrip22，…，Rastrip1i，Rastrip2i，…，Rastrip1n，Rastrip2n共計(jì)2n個(gè)帶鋼表面粗糙度的值。以壓下量為εi為例，將相關(guān)實(shí)驗(yàn)數(shù)據(jù)代入公式(30)得到方程組式(31)Rastrip1i=g1(hϵ)·g2(kϵ)·g(ϵi)·Rage1i+f1(hϵ)·f2(kϵ)·f(ϵi)·Rar1iRastrip2i=g1(hϵ)·g2(kϵ)·g(ϵi)·Rage2i+f1(hϵ)·f2(kϵ)·f(ϵi)·Rar2i---(31)]]>顯然，在方程組(31)中，由于帶材的來(lái)料厚度與材料強(qiáng)度不變，因此g1(hε)、g2(kε)、f1(hε)、f2(kε)等函數(shù)實(shí)際上是常數(shù)，相應(yīng)的g1(hε)·g2(kε)、f1(hε)·f2(kε)也就是常數(shù)，如果令C1=g1(hϵ)·g2(kϵ)C2=f1(hϵ)·f2(kϵ)---(32)]]>則將式(32)代入式(31)并整理得Rastrip1i=Rage1i·[C1·g(ϵi)]+Rar1i·[C2·f(ϵi)]Rastrip2i=Rage2i·[C1·g(ϵi)]+Rar2i·[C2·f(ϵi)]---(33)]]>由于在實(shí)驗(yàn)數(shù)據(jù)采集過(guò)程中，不同組實(shí)驗(yàn)數(shù)據(jù)在軋制時(shí)工作輥的原始粗糙度以及換輥后的軋制公里數(shù)是已知的，因此相應(yīng)的平整前帶材表面粗糙度與工作軋輥粗糙度Rage1i，Rage2i，Rar1i，Rar2i的值是已知的，而帶鋼表面粗糙度Rastrip1i，Rastrip2i通過(guò)實(shí)測(cè)得來(lái)，也是已知的。這樣，如果以C1·g(εi)、C2·f(εi)為未知數(shù)，求解方程式(33)得C1·g(ϵi)=Rar2i·astrip1i-Rar1i·Rastrip2iRar2i·Rage1i-Rar1i·Rage2iC2·f(ϵi)=Rage2i·astrip1i-Rage1i·Rastrip2iRage2i·Rar1i-Rage1i·Rar2i---(34)]]>隨后，在步驟23中，如果令G(εi)＝C1·g(εi)、F(εi)＝C2·f(εi)，將其代入式(34)得G(ϵi)=Rar2i·astrip1i-Rar1i·Rastrip2iRar2i·Rage1i-Rar1i·Rage2iF(ϵi)=Rage2i·astrip1i-Rage1i·Rastrip2iRage2i·Rar1i-Rage1i·Rar2i---(35)]]>接著，在步驟24中，將相應(yīng)的Rage1i，Rage2i，Rar1i，Rar2i以及Rastrip1i，Rastrip2i的值代入式(35)，并以εi為橫坐標(biāo)，對(duì)應(yīng)的G(εi)、F(εi)為縱坐標(biāo)，作出相應(yīng)的坐標(biāo)圖，然后進(jìn)行曲線(xiàn)擬合，求出函數(shù)G(ε)、F(ε)的表達(dá)式。
接著，在步驟25中，選擇強(qiáng)度為kh、壓下量為εh，厚度為h1，h2，…，hi，…，hn的第二組實(shí)驗(yàn)數(shù)據(jù)。
隨后，在步驟26中，選擇厚度為hk、壓下量為εk，強(qiáng)度為k1，k2，…，ki，…，kn的第三組實(shí)驗(yàn)數(shù)據(jù)。
接著，在步驟27中，令D1=g2(kh)·g(ϵh)D2=f2(kh)·f(ϵh)G1(hi)=D1·g1(hi)F1(hi)=D2·f1(hi)---(36)]]>E1=g1(hk)·g(ϵk)E2=f1(hk)·f(ϵk)G2(ki)=E1·g2(ki)F2(ki)=E2·f2(ki)---(47)]]>從而定義出相應(yīng)的中間函數(shù)G1(h)、F1(h)、G2(h)、F2(h)。
隨后，在步驟28中，采用與步驟24相同的方法，根據(jù)實(shí)驗(yàn)數(shù)據(jù)求出給出G1(h)、F1(h)、G2(h)、F2(h)的具體表達(dá)式。
接著，在步驟29中根據(jù)G(ε)、F(ε)、G1(h)、F1(h)、G2(h)、F2(h)求出Krs、Kss的具體表達(dá)式，方法如下根據(jù)相關(guān)函數(shù)定義有G(ε)＝g1(hε)·g2(kε)·g(ε) (48)G1(h)＝g2(kh)·g(εh)·g1(h) (49)G2(k)＝g1(hk)·g(εk)·g2(k) (50)如果將式(48)、式(49)、式(50)左右兩邊相乘得G(ε)·G1(h)·G2(k)＝g1(hε)·g2(kε)·g(ε)·g2(kh)·g(εh)(51)·g1(h)·g1(hk)·g(εk)·g2(k)整理式(51)得G(ε)·G1(h)·G2(k)＝[g1(hε)·g2(kε)·g(εh)]· (52)[g1(hk)·g(εk)·g2(kh)]·g1(h)·g2(k)·g(ε)顯然，由式(48)、式(50)得g1(hε)·g2(kε)·g(εh)＝G(εh) (53)g1(hk)·g(εk)·g2(kh)＝G2(kh) (54)將式(53)、式(54)代入式(52)，并整理得到遺傳率Kss的表達(dá)式Kss=g1(h)·g2(k)·g(ϵ)=G(ϵ)·G1(h)·G2(k)G(ϵh)·G2(kh)---(55)]]>同樣，采用類(lèi)似的方法，可以得到壓印率的Krs的表達(dá)式
Krs=f1(h)·f2(k)·f(ϵ)=F(ϵ)·F1(h)·F2(k)F(ϵh)·F2(kh)---(56)]]>隨后，在步驟30中，將式(55)、式(56)代入式(30)就可以得到相應(yīng)的實(shí)用特定的平整機(jī)組成品板面粗糙度模型，如式(57)所示Rastrip=G(ϵ)·G1(h)·G2(k)G(ϵh)·G2(kh)·Rage+F(ϵ)·F1(h)·F2(k)F(ϵh)·F2(kh)·Rar0·e-BL·L---(57)]]>實(shí)施例2一般而言，對(duì)于特定平整機(jī)而言，建立成品板面粗糙度模型的目的之一是為了實(shí)現(xiàn)根據(jù)已知軋制工藝參數(shù)對(duì)成品的板面粗糙度進(jìn)行在線(xiàn)預(yù)報(bào)。
圖3是按照本發(fā)明一個(gè)較佳實(shí)施例的平整機(jī)成品板面粗糙度在線(xiàn)預(yù)報(bào)計(jì)算流程圖。該實(shí)施例被用于特定平整機(jī)成品板面粗糙度的在線(xiàn)預(yù)報(bào)，并且關(guān)于該平整機(jī)在平整前入口帶材表面粗糙度與工作輥原始粗糙度及換輥后的軋制公里數(shù)已知。
首先，在步驟31中，收集待平整帶鋼的主要特征參數(shù)，包括帶鋼的厚度h、強(qiáng)度k、原始表面粗糙度Rage等。
隨后，在步驟32中，收集平整機(jī)工作輥主要特征參數(shù)包括原始粗糙度Rar0與軋制公里數(shù)L。
接著，在步驟33中，收集主要平整軋制工藝參數(shù)，主要是指平整的延伸率。
隨后，在步驟34中，調(diào)用相應(yīng)的工作輥粗糙度模型，利用下式求出相應(yīng)的工作輥表面粗糙度Rar=Rar0·e-BL·L---(58)]]>式(58)中，Rar為平整機(jī)工作輥表面粗糙度，Rar0為平整機(jī)工作輥原始粗糙度，L為工作輥換輥后的軋制公里數(shù)，BL為工作輥粗糙度衰減系數(shù)。
接著，在步驟35中，調(diào)平整機(jī)成品板面粗糙度預(yù)報(bào)模型。
最后，在步驟36中，將平整帶鋼參數(shù)、工作輥特征參數(shù)、平整軋制工藝參數(shù)等代入平整機(jī)成品板面粗糙度模型，計(jì)算相應(yīng)成品板面粗糙度，實(shí)現(xiàn)成品板面粗糙度的在線(xiàn)預(yù)報(bào)。
實(shí)施例3對(duì)于特定平整機(jī)而言，建立成品板面粗糙度模型的除了能夠?qū)崿F(xiàn)根據(jù)已知軋制工藝參數(shù)對(duì)成品的板面粗糙度進(jìn)行在線(xiàn)預(yù)報(bào)之外，還有一個(gè)重要的目的就是對(duì)成品板面粗糙度進(jìn)行控制，這也是本發(fā)明的意義所在。
圖4是按照本發(fā)明一個(gè)較佳實(shí)施例的平整機(jī)成品板面粗糙度在線(xiàn)控制流程圖。該實(shí)施例被用于特定平整機(jī)成品板面粗糙度的在線(xiàn)控制，并且關(guān)于該平整機(jī)在平整前入口帶材表面粗糙度與工作輥原始粗糙度及換輥后的軋制公里數(shù)、換輥周期等已知。
首先，在步驟41中，收集待平整帶鋼平整前的板面粗糙度Rage、被平整帶鋼的強(qiáng)度k、厚度等主要參數(shù)；隨后，在步驟42中，收集工作輥換輥后軋制公里數(shù)L、工藝技術(shù)規(guī)程所規(guī)定的工作輥軋制公里數(shù)(即換輥周期)Lmax0等軋輥主要特征參數(shù)；接著，在步驟43中，收集用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip以及用戶(hù)允許最大延伸率εmax、最小延伸率εmin等需求參數(shù)；隨后，在步驟44中，按照現(xiàn)場(chǎng)經(jīng)驗(yàn)設(shè)定工作輥的原始粗糙度的初始值Rar00；接著，在步驟45中，引入描述軋制公里數(shù)的相關(guān)系數(shù)α、β，并給定初始值α＝1.0、β＝0；隨后，在步驟46中，令Rar0＝Rar00，而軋制公里數(shù)分別取L1＝α·Lmax0、L2＝β·Lmax0(代表最小軋制公里數(shù)，新輥時(shí)最小軋制公里數(shù)為零)。
接著，在步驟47中，依靠工作輥粗糙度模型計(jì)算出相應(yīng)的工作輥實(shí)際可能最小粗糙度Rarmin與最大粗糙度Rarmax；隨后，在步驟48中，將平整前的板面粗糙度Rage、帶鋼強(qiáng)度k、工作輥實(shí)際可能最大粗糙度Rarmax、實(shí)際可能最小粗糙度Rarmin以及用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip代入成品板面粗糙度模型，反算出相應(yīng)的最大延伸率εmax′、最小延伸率εmin′；并判斷不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax是否同時(shí)成立，如果成立，則繼續(xù)步驟52)；如果不成立，調(diào)整工作輥的原始粗糙度的初始設(shè)定值Rar0，重復(fù)步驟44)、45)以及46)判斷能否找到合適的原始粗糙度設(shè)定值Rar0y，使得上述兩個(gè)不等式同時(shí)成立，如果可以找到則將工作輥原始粗糙度設(shè)定為Rar0y，繼續(xù)步驟52)；如無(wú)法合適的原始粗糙度值，則進(jìn)入步驟49)；接著，在步驟49中，判斷εmax′≥εmin是否成立，如果成立，則進(jìn)入步驟50，減小工作輥原始粗糙度的設(shè)定值至等式εmin′＝εmin成立，以一定的步長(zhǎng)逐步增大β的值，重復(fù)步驟45)、46)、47)，使得不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax同時(shí)成立，然后進(jìn)入步驟52)；如果不成立，則進(jìn)入步驟51，增大工作輥的原始粗糙度的初始設(shè)定值至等式εmax′＝εmax成立，以一定的步長(zhǎng)逐步減小α的值，重復(fù)步驟5)、6)、7)，使得不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax同時(shí)成立，然后進(jìn)入步驟52)；隨后，在步驟52中，將平整前的板面粗糙度Rage、帶鋼強(qiáng)度k、工作輥實(shí)際粗糙度Rar以及用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip代入成品板面粗糙度模型，反算出延伸率ε，進(jìn)行延伸率設(shè)定。同時(shí)注意，當(dāng)α<1.0]]>時(shí)，則在工作輥軋制公里數(shù)達(dá)到α·Lmax0時(shí)，及時(shí)轉(zhuǎn)換規(guī)格或換輥；而當(dāng)β>0]]>時(shí)，則應(yīng)該先安排生產(chǎn)其它規(guī)格的產(chǎn)品，直到工作輥軋制公里數(shù)已經(jīng)達(dá)到β·Lmax0時(shí)開(kāi)始生產(chǎn)該計(jì)劃規(guī)格產(chǎn)品。
最后，在步驟53中，輸出Rar0、ε、α、β，并按此值進(jìn)行相關(guān)參數(shù)的設(shè)定，完成成品板面粗糙度的設(shè)定計(jì)算。
以上借助實(shí)施例描述了本發(fā)明的具體實(shí)施方式
，但是應(yīng)該理解的是，這里具體的描述不應(yīng)理解為對(duì)本發(fā)明的范圍的限定，本領(lǐng)域內(nèi)的技術(shù)人員在閱讀本說(shuō)明書(shū)后對(duì)上述實(shí)施例作出的各種修改，這些都屬于本發(fā)明的范圍。
本發(fā)明的平整過(guò)程成品板面粗糙度控制技術(shù)，是通過(guò)對(duì)工作輥輥面粗糙度、軋制參數(shù)和所軋帶材表面粗糙度對(duì)應(yīng)關(guān)系的深入分析，在引入壓印率與遺傳率概念的基礎(chǔ)上，針對(duì)平整工藝的生產(chǎn)特點(diǎn)，經(jīng)過(guò)大量的現(xiàn)場(chǎng)試驗(yàn)與理論研究而建立的，在實(shí)踐中能夠?qū)崿F(xiàn)根據(jù)用戶(hù)對(duì)帶材表面粗糙度的要求適時(shí)調(diào)整軋制工藝參數(shù)以控制成品帶鋼的表面粗糙度。
權(quán)利要求
1.一種平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，包括以下步驟a、建立反映平整機(jī)工作輥粗糙度與工作輥原始粗糙度、換輥后軋制公里數(shù)之間關(guān)系的數(shù)學(xué)模型；按下述步驟建立反映平整機(jī)工作輥粗糙度與工作輥原始粗糙度、換輥后軋制公里數(shù)之間關(guān)系的數(shù)學(xué)模型1)以Rar=Rar0·e-BL·L]]>的形式構(gòu)造出平整機(jī)工作輥粗糙度與軋制公里數(shù)之間關(guān)系的基本模型框架；式中Rar-平整機(jī)工作輥粗糙度；Rar0-平整機(jī)工作輥原始粗糙度；L-工作輥換輥后的軋制公里數(shù)；LBL-工作輥粗糙度衰減系數(shù)；2)在現(xiàn)場(chǎng)收集一組關(guān)于不同軋制公里數(shù)Ll下所對(duì)應(yīng)的工作輥粗糙度Rari以及相應(yīng)的工作輥原始粗糙度Rar0i的數(shù)據(jù)；3)給定工作輥粗糙度衰減系數(shù)初始設(shè)定值X0＝{BL′}與計(jì)算精度∑＝0.00001；4)將BL′、Li、Rar0i代入公式Rar=Rar0·e-BL·L,]]>求出相應(yīng)的工作輥粗糙度理論計(jì)算值Rart′；5)計(jì)算出相應(yīng)的目標(biāo)函數(shù)值F(X)；6)判斷F(X)≤∑是否成立，若不成立改變初始值X0重復(fù)上述步驟4)、5)，直到判斷條件成立，結(jié)束計(jì)算，得出特定平整機(jī)的工作輥粗糙度衰減系數(shù)BL以及相應(yīng)的反映平整機(jī)工作輥粗糙度與工作輥原始粗糙度、換輥后軋制公里數(shù)之間關(guān)系的數(shù)學(xué)模型；b、構(gòu)造壓印率Krs與遺傳率Kss的函數(shù)表達(dá)式；c、構(gòu)造平整后成品板面粗糙度Rastrip的初始模型；d、根據(jù)特定平整機(jī)的現(xiàn)場(chǎng)實(shí)驗(yàn)數(shù)據(jù)，建立特定平整機(jī)的板面粗糙度預(yù)報(bào)模型，按下述步驟建立特定平整機(jī)的板面粗糙度模型1)、現(xiàn)場(chǎng)采集一組試驗(yàn)數(shù)據(jù)，滿(mǎn)足以下要求帶材的來(lái)料厚度都為hε，材料強(qiáng)度都為kε，而壓下量共有ε1，ε2，…，εi，…，εn共n個(gè)，而每個(gè)壓下量εi對(duì)應(yīng)兩個(gè)不同的平整前入口帶材表面粗糙度與工作輥軋輥粗糙度的值，關(guān)于這點(diǎn)，可以通過(guò)選擇不同的工作輥軋制公里數(shù)來(lái)實(shí)現(xiàn)，設(shè)其為Rage1i，Rage2i，Rar1i，Rar2i。這樣，與之對(duì)應(yīng)，通過(guò)現(xiàn)場(chǎng)實(shí)際測(cè)量就可以得出Rastrip11，Rastrip21，Rastrip12，Rastrip22，…，Rastrip1i，Rastrip2i，…，Rastrip1n，Rastrip2n共計(jì)2n個(gè)帶鋼表面粗糙度的值；2)定義中間函數(shù)G(ε)、F(ε)；3)根據(jù)實(shí)驗(yàn)數(shù)據(jù)求解出相應(yīng)的G(εi)、F(εi)的表達(dá)式；4)根據(jù)G(εi)、F(εi)的表達(dá)式通過(guò)曲線(xiàn)擬合的方法給出G(ε)、F(ε)的具體表達(dá)式；5)采集材料強(qiáng)度都為kh、壓下量都為εh，而帶材來(lái)料厚度分別為h1，h2，…，hi，…，hn的第二組實(shí)驗(yàn)數(shù)據(jù)與帶材來(lái)料厚度都為hk、壓下量都為εk，而材料強(qiáng)度分別為k1，k2，…，ki，…，kn的第三組實(shí)驗(yàn)數(shù)據(jù)；6)定義中間函數(shù)G1(h)、F1(h)、G2(h)、F2(h)；7)根據(jù)實(shí)驗(yàn)數(shù)據(jù)通過(guò)曲線(xiàn)擬合的方法給出G1(h)、F1(h)、G2(h)、F2(h)的具體表達(dá)式；8)根據(jù)G(ε)、F(ε)、G1(h)、F1(h)、G2(h)、F2(h)求出壓印率Krs與遺傳率Kss的具體函數(shù)表達(dá)式；9)輸出特定平整機(jī)的成品板面粗糙度預(yù)報(bào)模型；e、根據(jù)板面粗糙度模型在線(xiàn)預(yù)報(bào)成品板面粗糙度，按下述步驟實(shí)現(xiàn)成品板面粗糙度的預(yù)報(bào)1)收集待平整帶鋼的主要特征參數(shù)，所述平整帶鋼的主要特征參數(shù)包括平整前帶鋼表面粗糙度、帶鋼強(qiáng)度；2)收集平整機(jī)工作輥主要特征參數(shù)，所述平整機(jī)工作輥主要特征參數(shù)包括工作輥原始粗糙度、工作輥換輥后軋制公里數(shù)；3)收集主要平整軋制工藝參數(shù)，所述平整軋制工藝參數(shù)包括平整軋制延伸率；4)計(jì)算平整機(jī)工作輥粗糙度；5)將平整帶鋼參數(shù)、工作輥特征參數(shù)、平整軋制工藝參數(shù)等代入平整機(jī)成品板面粗糙度模型，計(jì)算相應(yīng)成品板面粗糙度；f、成品板面粗糙度的在線(xiàn)控制，按下述步驟實(shí)現(xiàn)成品板面粗糙度的在線(xiàn)控制1)收集待平整帶鋼平整前的板面粗糙度Rage、帶鋼強(qiáng)度k等主要參數(shù)；2)收集工作輥換輥后軋制公里數(shù)L、工藝技術(shù)規(guī)程所規(guī)定的工作輥軋制公里數(shù)(即換輥周期)Lmax0等軋輥主要特征參數(shù)；3)收集用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip以及用戶(hù)允許最大延伸率εmaz、最小延伸率εmin等需求參數(shù)；4)按照現(xiàn)場(chǎng)經(jīng)驗(yàn)設(shè)定工作輥的原始粗糙度的初始值Rar00；5)引入相關(guān)系數(shù)α、β，并給定初始值α＝1.0、β＝0；6)令Rar0＝Rar00，而軋制公里數(shù)分別取L1＝α·Lmax0、L2＝β·Lmax0代入步驟a)所述模型計(jì)算出相應(yīng)的工作輥實(shí)際可能最小粗糙度Rarmin與最大粗糙度Rarmax；代表最小軋制公里數(shù)，新輥時(shí)最小軋制公里數(shù)為零；7)將平整前的板面粗糙度Rage、帶鋼強(qiáng)度k、工作輥實(shí)際可能最大粗糙度Rarmax、實(shí)際可能最小粗糙度Rarmin以及用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip代入成品板面粗糙度模型，反算出相應(yīng)的最大延伸率εmax′、最小延伸率εmin′；8)判斷不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax是否同時(shí)成立，如果成立，則繼續(xù)步驟10)；如果不成立，調(diào)整工作輥的原始粗糙度的初始設(shè)定值Rar0，重復(fù)步驟4)、5)以及6)判斷能否找到合適的原始粗糙度設(shè)定值Rar0y，使得上述兩個(gè)不等式同時(shí)成立，如果可以找到則將工作輥原始粗糙度設(shè)定為Rar0y，繼續(xù)步驟10)；如無(wú)法合適的原始粗糙度值，則進(jìn)入步驟9)；9)判斷εmax′≥εmin是否成立，如果成立，則減小工作輥原始粗糙度的設(shè)定值至等式εmin′＝εmin成立，以一定的步長(zhǎng)逐步增大β的值，重復(fù)步驟5)、6)、7)，使得不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax同時(shí)成立，然后進(jìn)入步驟10)；如果不成立，則增大工作輥的原始粗糙度的初始設(shè)定值至等式εmax′＝εmax成立，以一定的步長(zhǎng)逐步減小α的值，重復(fù)步驟5)、6)、7)，使得不等式εmin≤εmax′≤εmax、εmin≤εmin′≤εmax同時(shí)成立，然后進(jìn)入步驟10)；10)將平整前的板面粗糙度Rage、帶鋼強(qiáng)度k、工作輥實(shí)際粗糙度Rar以及用戶(hù)對(duì)成品帶鋼板面粗糙度的需求Rastrip代入成品板面粗糙度模型，反算出延伸率ε，進(jìn)行延伸率設(shè)定。同時(shí)注意，當(dāng)α＜1.0時(shí)，則在工作輥軋制公里數(shù)達(dá)到α·Lmax0時(shí)，及時(shí)轉(zhuǎn)換規(guī)格或換輥；而當(dāng)β＞0時(shí)，則應(yīng)該先安排生產(chǎn)其它規(guī)格的產(chǎn)品，直到工作輥軋制公里數(shù)已經(jīng)達(dá)到β·Lmax0時(shí)開(kāi)始生產(chǎn)該計(jì)劃規(guī)格產(chǎn)品。
2.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟a中所述目標(biāo)函數(shù)定義為F(X)=Σi=1n(Rari-Rari′)2]]>式中n-樣本數(shù)。
3.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟b中所述壓印率Krs與遺傳率Kss的函數(shù)表達(dá)式為Krs=f1(h)·f2(k)·f(ϵ)Kss=g1(h)·g2(k)·g(ϵ)]]>式中h-所平整帶材的來(lái)料厚度k-所平整帶材的強(qiáng)度ε-平整壓下率。
4.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟c中所述平整板面粗糙度的初始模型表達(dá)式如下Rastrip=Kss·Rage+Krs·Rar0·e-BL·L]]>式中Rastrip-平整后成品帶鋼表面粗糙度Rage-平整前來(lái)料帶材表面粗糙度。
5.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟d中所述中間函數(shù)G(ε)、F(ε)的初始表達(dá)式如下G(ε)＝g1(hε)·g2(kε)·g(ε)f(ε)＝f1(hε)·f2(kε)·f(ε)。
6.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟d中所述求解出的G(εi)、F(εi)的表達(dá)式如下G(ϵi)=Rar2i·astrip1i-Rar1i·Rastrip2iRar2iRage1i-Rar1i·Rage2iF(ϵi)=Rage2i·Rastrip1i-Rage1i·Rastrip2iRage2i·Rar1i-Rage1i·Rar2i.]]>
7.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟d中所述中間函數(shù)G1(h)、F1(h)、G2(h)、F2(h)的表達(dá)式如下G1(h)=g2(kh)·g(ϵh)·g1(h)F1(h)=f2(kh)·f(ϵh)·f1(h)G2(k)=g1(hk)·g(ϵk)·g2(k)F2(k)=f1(hk)·f(ϵk)·f2(k).]]>
8.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟d中所述印率Krs與遺傳率Kss的具體函數(shù)表達(dá)式如下Krs=f1(h)·f2(k)·f(ϵ)=F(ϵ)·F1(h)·F2(k)F(ϵh)·F2(kh)]]>Kss=g1(h)·g2(k)·g(ϵ)=G(ϵ)·G1(h)·G2(k)G(ϵh)·G2(kh).]]>
9.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟d中所述成品板面粗糙度數(shù)學(xué)模型如下Rastrip=G(ϵ)·G1(h)·G2(k)G(ϵh)·G2(kh)·Rage+F(ϵ)·F1(h)·F2(k)F(ϵh)·F2(kh)·Rar0·e-BL·L.]]>
10.根據(jù)權(quán)利要求1所述的平整機(jī)的成品板面粗糙度預(yù)報(bào)與控制技術(shù)，其特征在于步驟f中所述α、β平整機(jī)軋制計(jì)劃系數(shù)，當(dāng)α＜1.0時(shí)，應(yīng)將特定規(guī)格產(chǎn)品安排在換輥初期軋制；而當(dāng)β＞0時(shí)，則應(yīng)該安排換輥后期軋制。
全文摘要
一種平整機(jī)成品板面粗糙度預(yù)報(bào)與控制技術(shù)，包括以下步驟a.建立反映平整機(jī)工作輥粗糙度與工作輥原始粗糙度、換輥后軋制公里數(shù)之間關(guān)系的數(shù)學(xué)模型；b.構(gòu)造壓印率K
文檔編號(hào)G05B13/04GK1804739SQ20051004818
公開(kāi)日2006年7月19日 申請(qǐng)日期2005年12月12日 優(yōu)先權(quán)日2005年12月12日
發(fā)明者白振華, 周慶田 申請(qǐng)人:燕山大學(xué)