楊光輝， 張杰， 曹建國， 李洪波， 黃橋?qū)?/p>

(北京科技大學(xué) 機(jī)械工程學(xué)院， 北京 100083)

?

2 180 mm冷連軋機(jī)各機(jī)架對成品板形的影響

楊光輝， 張杰， 曹建國， 李洪波， 黃橋?qū)?/p>

(北京科技大學(xué) 機(jī)械工程學(xué)院， 北京 100083)

為了分析冷連軋機(jī)機(jī)架間的相互影響以及各機(jī)架對成品板形的影響，基于冷連軋凸度和平坦度模型，采用影響系數(shù)法建立了冷連軋機(jī)的板形調(diào)控功效模型. 針對某2 180 mm冷連軋機(jī)組，結(jié)合其4個典型軋制規(guī)程，對不同機(jī)架板形調(diào)控功效進(jìn)行了具體分析，可以看出，對于同一軋制規(guī)程，各機(jī)架對成品板形的影響差別很大，越靠近末機(jī)架，影響越大. 在冷連軋機(jī)的5個機(jī)架中，除了末機(jī)架即第5機(jī)架外，對成品板形影響最大的為第4機(jī)架，最大影響接近10 IU(1 IU=10-5)，且不同機(jī)架各板形調(diào)控手段對成品板形的影響系數(shù)在機(jī)架間相差10-3數(shù)量級. 最后，對板形數(shù)據(jù)進(jìn)行了統(tǒng)計分析及相關(guān)軋制試驗(yàn)來分析結(jié)論的可靠性.

冷軋機(jī)；帶鋼；成品板形；影響系數(shù)；機(jī)架

冷連軋的一個顯著特點(diǎn)就是機(jī)械、電氣、液壓控制系統(tǒng)和軋件變形之間存在緊密聯(lián)系，多機(jī)架連軋機(jī)由于帶鋼的聯(lián)系而形成一個系統(tǒng)，在各機(jī)架間帶鋼起著傳遞厚度、凸度、平坦度、溫度以及張力的作用. 前一個機(jī)架出口的板形影響著后一個機(jī)架的板形，各種擾動對前一機(jī)架出口板形的影響都將經(jīng)過一定的時間后，由被影響的那一段帶鋼帶到下一個機(jī)架而成為新的外擾源. 每一個外擾量產(chǎn)生的凸度波動、平坦度波動等都將成為新的外擾源. 因此，各種原始的和再生的外擾將在各個機(jī)架間相互影響[1-5]. 冷連軋機(jī)組中某一個機(jī)架參數(shù)發(fā)生變化后，除直接影響本機(jī)架的工作外，還將影響其他機(jī)架的工作，最終都會在成品的板形上表現(xiàn)出來. 但當(dāng)出現(xiàn)某些典型板形缺陷時，由于無法定量評價冷連軋機(jī)機(jī)架間的相互影響以及各機(jī)架對成品板形的影響，因此無法對軋制參數(shù)進(jìn)行更為精確合理的設(shè)置及補(bǔ)償，影響對板形的控制. 所以，需要建立冷連軋機(jī)組的各機(jī)架板形調(diào)控功效模型，研究冷連軋機(jī)機(jī)架間的相互影響以及對成品板形的影響.

本文通過建立冷連軋機(jī)各機(jī)架板形調(diào)控手段對成品板形的調(diào)控功效模型，針對某2 180 mm冷連軋機(jī)組，結(jié)合其4個典型的軋制規(guī)程，對冷連軋機(jī)各機(jī)架對成品機(jī)架板形調(diào)控影響進(jìn)行分析.

1 軋機(jī)出口板形的影響因素

某2 180 mm冷連軋機(jī)組由5機(jī)架6輥CVC軋機(jī)組成，如圖1所示. 該軋機(jī)裝備了自動厚度控制、自動板形控制、軋制速度控制等多組自動化控制系統(tǒng)，是中國目前最先進(jìn)的軋機(jī)之一. 其裝備了豐富的板形調(diào)控和檢測手段，有軋輥傾斜(S1～S5)、工作輥竄輥(S1～S5)、工作輥正負(fù)彎輥(S1～S5)、中間輥正負(fù)彎輥(S1～S5)、中間輥竄輥(S1～S5)和工作輥精細(xì)分段冷卻系統(tǒng)(S5). 在第1和第5機(jī)架后配備有先進(jìn)板形儀[6-7]. 可以生產(chǎn)最薄0.3 mm、寬至2 080 mm的帶鋼，設(shè)計年產(chǎn)最為2.15×106t.

2 180 mm冷連軋機(jī)組雖配備了先進(jìn)的軋機(jī)，具有豐富的板形控制手段，但生產(chǎn)的產(chǎn)品仍存在一定的板形質(zhì)量問題，有時板形不良甚至?xí)绊懞蠊ば虻捻樌a(chǎn)，造成跑偏事故. 現(xiàn)場生產(chǎn)中多出現(xiàn)中浪、邊中復(fù)合浪、復(fù)雜局部浪形等問題，如圖2所示. 在實(shí)際生產(chǎn)中，對板形的控制多依賴于控制系統(tǒng)的反饋控制，當(dāng)成品帶鋼出現(xiàn)板形缺陷時，無法在離線時評價各機(jī)架對成品板形的調(diào)控功效，從而不能對軋制參數(shù)進(jìn)行更合理的設(shè)置. 而目前應(yīng)用較為廣泛的模型又過于復(fù)雜，系數(shù)難以確定，實(shí)際運(yùn)用不太方便.

1.1 凸度和平坦度模型

為了研究冷連軋機(jī)的3種主要板形調(diào)控手段即工作輥彎輥、中間輥彎輥和中間輥竄輥對帶鋼板形的影響，同時考慮帶鋼入口板形(凸度、平坦度)的波動所引起出口凸度的改變，故需要建立含這5種變量發(fā)生變化時的板凸度方程. 帶鋼板凸度與影響參數(shù)的理論模型為

Ci=Cbi+KciCi-1+Kρiρi-1+

(1)

式中：下標(biāo)i為機(jī)架號；Ci為出口板凸度，mm；Cbi為基本板凸度，mm；Ci-1為入口板凸度，mm；ρi-1為入口平坦度，105IU；Sfi為中間輥竄輥量，mm；Fmi為中間輥彎輥力，kN；Fwi為工作輥彎輥力，kN；Kci為入口板凸度的影響系數(shù)，mm/mm；Kρi為入口平坦度的影響系數(shù)，mm/105IU；Ksi為中間輥竄輥的影響系數(shù)，mm/mm；Kmi為中間輥彎輥的影響系數(shù)，mm/kN；Kwi為工作輥彎輥的影響系數(shù)，mm/kN.

在忽略軋制區(qū)金屬橫向流動并假設(shè)軋機(jī)入口帶鋼被張力拉伸平直的條件下，以縱向纖維相對長度差表達(dá)的冷軋帶鋼平坦度模型為[8]

(2)

式中：ρi-1、ρi分別為第i機(jī)架帶鋼入口和出口的平坦度；Ci-1、Ci分別為第i機(jī)架帶鋼入口和出口的凸度；Hmi、hmi分別為第i機(jī)架帶鋼入口和出口的厚度.

1.2 板形調(diào)控功效模型

影響系數(shù)法是研究在擾動量或控制量作用下，冷連軋機(jī)的連軋過程從一個穩(wěn)態(tài)過渡到另一個穩(wěn)態(tài)后各參數(shù)變化量之間關(guān)系的一種行之有效的方法[9]. 假設(shè)在初始情況下，壓下量一定、前后張力一定、輥形(包括磨削輥形、磨損輥形、熱凸度輥形等)一定，工作輥彎輥力、中間輥彎輥力、中間輥竄輥量均為0，來料凸度為0、平坦度良好，此時軋制出來的板凸度稱為基本板凸度，用Cb表示. 上述任一參數(shù)的改變都可能會導(dǎo)致板形產(chǎn)生一個新的值，稱為再生板凸度. 各影響系數(shù)是指出口板凸度變化量同影響參數(shù)變化量的比值，其表征任一參數(shù)的單位改變量所引起的板凸度變化，即

(3)

式中:K為影響系數(shù)；δX為某個影響參數(shù)的變化量；δC為出口板凸度的變化量.

由以上分析可知，對于冷連軋機(jī)5個機(jī)架中的任何一個機(jī)架，工作輥彎輥、中間輥彎輥、中間輥竄輥、帶鋼入口凸度和入口平坦度等因素對帶鋼板形的影響系數(shù)可表示為如表1所示.

表1 各項(xiàng)因素對軋機(jī)出口板形的影響系數(shù)

Tab.1 Influence coefficient of various effect factors on strip shape of mill exit

參數(shù)無超聲300W超聲1000W超聲ZnMgCuZnMgCuZnMgCu晶內(nèi)質(zhì)量分?jǐn)?shù)3.6030.9050.6504.3091.2360.8014.9781.2800.911晶界質(zhì)量分?jǐn)?shù)21.8859.88012.33013.9706.21011.42011.4155.9509.580平均質(zhì)量分?jǐn)?shù)7.2252.7002.9866.2412.2312.9256.2652.2142.645

2 不同機(jī)架對成品板形調(diào)控功效模型的建立

在冷連軋機(jī)組中，某一機(jī)架的板形調(diào)控手段所造成的本機(jī)架帶鋼出口板形的變化，通過帶鋼的傳遞作用，必然會引起下一機(jī)架的入口板形波動，進(jìn)而造成下一機(jī)架的出口板形發(fā)生變化. 為了研究冷連軋機(jī)的某一機(jī)架的板形調(diào)控手段對成品板形的影響，需對各機(jī)架的板形調(diào)控功效模型進(jìn)行機(jī)架間迭代分析研究.

由于冷軋的板形主要是指平坦度，所以本文中以下所提到的各機(jī)架各板形調(diào)控手段對成品板形的調(diào)控功效，均指對成品平坦度的調(diào)控功效.

2.1 第5機(jī)架對成品板形調(diào)控功效

第5架的工作輥彎輥、中間輥彎輥、中間輥竄輥以及入口板形等不同影響因素對其軋機(jī)出口板形的影響系數(shù)及表達(dá)式如表2所示.

表2 第5架的各項(xiàng)因素對其軋機(jī)出口板形的影響系數(shù)及表達(dá)式

Tab.2 Influence coefficient of various effect factors of No.5 stand mill on its exit strip shape

影響因素對出口平坦度的影響系數(shù)對出口平坦度的影響系數(shù)表達(dá)式工作輥彎輥Kρ5w5-Kw5/hm5中間輥彎輥Kρ5m5-Km5/hm5中間輥竄輥Kρ5s5-Ks5/hm5入口凸度Kρ5c51/Hm5-Kc5/hm5入口平坦度Kρ5ρ51-Kρ5/hm5

2.2 第4機(jī)架對成品板形調(diào)控功效

2.2.1 中間輥竄輥

對于第4機(jī)架，中間輥竄輥會引起出口凸度和出口平坦度的變化，從而造成第5機(jī)架的入口凸度和入口平坦度的變化. 因此，首先需要分析第5機(jī)架入口板形對成品板形的影響. 在入口凸度和平坦度都有單位波動的情況下，所造成的成品平坦度波動為

(4)

第4機(jī)架中間輥竄輥對本機(jī)架出口凸度影響系數(shù)為Ks4，對本機(jī)架出口平坦度影響系數(shù)為-Ks4/hm4，將其代入式(4)，可得第4機(jī)架中間輥竄輥對成品板形的調(diào)控功效為

(5)

2.2.2 中間輥彎輥的影響

同理，可以求出中間輥彎輥對成品板形的調(diào)控功效，即

(6)

2.2.3 工作輥彎輥

同理，可以求出工作輥彎輥對成品板形的調(diào)控功效，即

(7)

2.2.4 第4機(jī)架入口板形對成品板形的影響

對于第4機(jī)架，入口凸度會引起出口凸度和出口平坦度的變化，從而造成第5機(jī)架的入口凸度和入口平坦度的變化. 所以，首先要分析第5機(jī)架入口板形對成品板形的影響. 對于第5機(jī)架，在入口凸度和平坦度都有單位波動的情況下，所造成的成品平坦度波動如式(4)所示.

其次要分析第4機(jī)架入口凸度對本機(jī)架出口板形的影響. 第4機(jī)架入口凸度對本機(jī)架出口凸度影響系數(shù)為Kc4，對本機(jī)架出口平坦度影響系數(shù)為1/Hm4-Kc4/hm4，將其代入式(4)，可得第4機(jī)架入口凸度對成品平坦度的影響系數(shù)為

(8)

同理，可得第4機(jī)架入口平坦度對成品平坦度的影響系數(shù)，即

(9)

同理可以推出第3、第2和第1機(jī)架的各板形調(diào)控手段對成品板形調(diào)控功效的影響以及入口板形對成品板形的影響.

3 典型軋制規(guī)程的實(shí)例分析

3.1 各軋制規(guī)程主要工藝參數(shù)

為了能夠更好地研究冷連軋機(jī)各機(jī)架板形調(diào)控功效，本文在參考某2 180 mm冷連軋機(jī)組軋制產(chǎn)品大綱的基礎(chǔ)上，選取了其中較為典型的4個軋制規(guī)程，其主要工藝參數(shù)如表3所示，各機(jī)架的帶鋼厚度和寬度分布如表4所示. 其中，軋制規(guī)程1和軋制規(guī)程2的寬度為1 000 mm，軋制規(guī)程3和軋制規(guī)程4的寬度為2 000 mm.

表3 軋制規(guī)程主要工藝參數(shù)

表4 各機(jī)架軋后帶鋼厚度

3.2 各機(jī)架不同板形調(diào)控手段對成品板形的影響

通過冷連軋機(jī)不同機(jī)架板形調(diào)控功效模型計算得到不同板形調(diào)控手段對成品板形的影響系數(shù)如表5所示，表中下標(biāo)KWb表示工作輥彎輥；KIb表示中間輥彎輥；KIs表示中間輥竄輥.

表5 各機(jī)架不同板形調(diào)控手段對成品板形的影響系數(shù)

可以看出，從首機(jī)架(S1)到成品機(jī)架(S5)，各機(jī)架板形調(diào)控手段對成品板形的影響系數(shù)越來越大. 因此，除了成品機(jī)架外，對成品板形影響最大的是第4機(jī)架. 如當(dāng)進(jìn)行滿量程調(diào)控時，4個軋制規(guī)程中的第4機(jī)架中間輥竄輥對成品板形的影響范圍分別為1.32、1.07、4.88和5.00 IU，可見第4機(jī)架中間輥竄輥對成品板形具有有一定的影響，經(jīng)計算得到第4機(jī)架3種板形控制手段對成品板形的最大影響約為10 IU.

而對于前3個機(jī)架，各板形調(diào)控手段對成品板形的影響系數(shù)的數(shù)量級很小，分別為10-12、10-9和10-6，也就是說即便4個軋制規(guī)程均進(jìn)行滿量程調(diào)控，其對成品板形的影響也很小，可以忽略. 這是因?yàn)楦鳈C(jī)架對成品板形的影響是通過帶鋼的傳遞作用實(shí)現(xiàn)的，即通過對本機(jī)架出口板形的影響，然后經(jīng)過機(jī)架間的迭代才能轉(zhuǎn)化為對成品板形的影響，而每經(jīng)過一個機(jī)架，都會有一定的削弱(不同機(jī)架各板形調(diào)控手段對成品板形的影響系數(shù)在機(jī)架間相差10-3數(shù)量級)，所以，越靠近成品機(jī)架，影響越大，越靠近首機(jī)架，影響越小.

3.3 各機(jī)架不同板形調(diào)控手段對本機(jī)架出口板形的影響

表6為冷連軋機(jī)各機(jī)架不同板形調(diào)控手段對本機(jī)架出口板形的影響系數(shù). 可以看出，各機(jī)架不同板形調(diào)控手段對本機(jī)架出口板形的影響系數(shù)基本處在同一數(shù)量級. 這是因?yàn)槔溥B軋機(jī)的5個機(jī)架有著相同或相近的工藝參數(shù)和調(diào)控手段，所以每個機(jī)架的板形調(diào)控手段對自身出口板形的影響類似.

表6 各機(jī)架不同板形調(diào)控手段對本機(jī)架出口板形的影響系數(shù)

Tab.6 Influence coefficients of different shape adjustment means of every stand on its own exit shape

軋制規(guī)程1234系數(shù)S1S2S3S4S5KWb/(IU·kN-1)1.201.601.902.402.40KIb/(IU·kN-1)0.590.740.921.101.10KIs/(IU·kN-1)0.650.770.911.101.10KWb/(IU·kN-1)1.001.301.501.901.80KIb/(IU·kN-1)0.490.60.720.870.86KIs/(IU·kN-1)0.560.650.740.860.82KWb/(IU·kN-1)2.503.103.904.904.90KIb/(IU·kN-1)1.181.501.802.302.30KIs/(IU·kN-1)2.573.103.604.304.20KWb/(IU·kN-1)2.002.503.003.703.70KIb/(IU·kN-1)0.981.201.401.801.70KIs/(IU·kN-1)2.242.603.003.503.30

4 板形數(shù)據(jù)的統(tǒng)計分析以及軋制試驗(yàn)

為明確超寬軋機(jī)產(chǎn)品板形質(zhì)量，將2 180 mm冷連軋機(jī)連續(xù)生產(chǎn)的500卷帶鋼的板形數(shù)據(jù)進(jìn)行分析和統(tǒng)計，表7和圖2為統(tǒng)計結(jié)果.

從表7中可以看出此批帶鋼綜合平坦度的在5IU以下的有85.7%，6IU以下的有98.1%，達(dá)到對板形質(zhì)量的基本要求. 圖3中的橫坐標(biāo)1～8分別代表OS(operating side)邊浪、DS(driving side)邊浪、中浪、雙邊浪、DS肋浪加OS單邊浪、OS肋浪加DS單邊浪、1/4浪和邊中復(fù)合浪，從圖3中可以看出此批帶鋼的板形缺陷以“邊中復(fù)合浪”為主占64.4%，“中浪”占16.7%和“OS單側(cè)肋浪加DS單邊浪”占12.9%，其他板形缺陷所占比例較小，均未超過5.0%.

表7 綜合平坦度分檔統(tǒng)計結(jié)果

為了研究冷連軋機(jī)各機(jī)架對成品板形的調(diào)控功效，從上述500卷帶鋼中挑出有代表性的2卷帶鋼進(jìn)行了現(xiàn)場在線軋制試驗(yàn)，一卷寬度為1 000 mm、主要板形缺陷為中浪，另一卷寬度為2 000 mm、主要板形缺陷為邊中復(fù)合浪. 對于寬度為1 000 mm的帶鋼，在其他各機(jī)架軋制參數(shù)保持不變的情況下，在帶鋼縱向位置1 700～1 900 m處將冷軋機(jī)的第4機(jī)架的工作輥彎輥力增大250 kN；對于寬度為2 000 mm的帶鋼，其他各機(jī)架的軋制參數(shù)保持不變，在帶鋼縱向位置700～900 m處將第4機(jī)架工作輥彎輥增大280 kN. 通過對軋制過程和最終產(chǎn)品的觀察和測量可以看到，第4機(jī)架工作輥彎輥力的變化對成品板形造成了一定的影響，如圖4所示. 可以看出，第4機(jī)架各板形控制手段對成品板形具有一定的影響，該影響是控制量還是干擾量應(yīng)該由板形控制策略以及軋機(jī)的實(shí)際板形控制能力等方面決定. 如果第5機(jī)架各板形控制手段的調(diào)控能力不足，那么第4機(jī)架可以起到一定的補(bǔ)充作用. 如果第5機(jī)架板形調(diào)控能力充足，那么該影響是一種干擾量，在板形調(diào)節(jié)過程中要對第5機(jī)架進(jìn)行一定的補(bǔ)償以抵消這種干擾.

5 結(jié) 論

對于同一軋制規(guī)程，各機(jī)架對成品板形的影響差別很大，越靠近成品機(jī)架，影響越大. 除了第5機(jī)架外，對成品板形影響最大的為第四機(jī)架，最大影響接近10 IU.

各機(jī)架對成品板形的影響是通過帶鋼的傳遞作用實(shí)現(xiàn)的，即通過對本機(jī)架出口板形的影響，然后經(jīng)過機(jī)架間的迭代才能轉(zhuǎn)化為對成品板形的影響，而每經(jīng)過一個機(jī)架，都會有一定的削弱，不同板形調(diào)控手段對成品板形的影響系數(shù)在機(jī)架間相差10-3數(shù)量級，所以，越靠近末機(jī)架，影響越大，越靠近首機(jī)架，影響越小.

[1] Pittner John， Simaan Marwan A. Control of a continuous tandem cold metal rolling process[J]. Control Engineering Practice， 2008，16(11):1379-1390.

[2] Venkata Reddy N， Suryanarayana G. A set-up model for tandem cold rolling mills[J]. Journal of Materials Processing Technology， 2001，116(2-3):269-277.

[3] 張云鵬，吳慶海，王長松.六輥CVC冷軋機(jī)板形控制性能研究[J].冶金設(shè)備，1998(6)：8-10.

Zhang Yunpeng， Wu Qinghai， Wang Changsong. Study of flatness control on 6-h CVC mill[J]. Metallurgical Equipment， 1998(6):8-10. (in Chinese)

[4] 王曉晨，楊荃.萬能凸度軋機(jī)非對稱板形調(diào)控特型與廣義整體板形調(diào)控策略研究[J].機(jī)械工程學(xué)報，2012，48(4)：58-65.

Wang Xiaochen， Yang Quan. Research on asymmetric shape control character and generalized overall shape control strategy for universal crown mill[J]. Journal of Mechanical Engineering， 2012，48(4):58-65. (in Chinese)

[5] 閆沁太，張杰，賈生暉，等.冷軋機(jī)板形調(diào)節(jié)能力分析方法的研究與應(yīng)用[J].機(jī)械工程學(xué)報，2011，47(4)：77-81.

Yan Qintai， Zhang Jie， Jia Shenghui， et al. Research on analyzing the flatness adjusting capacity for cold mill and its application[J]. Journal of Mechanical Engineering， 2011，47(4):77-81. (in Chinese)

[6] 楊光輝，張杰，曹建國，等.氣流激振及渦流測幅式板形儀帶鋼振幅與板形關(guān)系[J].北京理工大學(xué)學(xué)報，2015，35(7):671-676.

Yang Guanghui， Zhang Jie， Cao Jianguo， et al. Relationship between strip amplitude and shape for shapometer based on airflow excitation and eddy current[J]. Transactions of Beijing Institute of Technology， 2015，35(7):671-676. (in Chinese)

[7] 楊光輝，張杰，曹建國，等.SI-FLAT板形儀激振頻率設(shè)定[J].天津大學(xué)學(xué)報，2014，47(10):871-878.

Yang Guanghui， Zhang Jie， Cao Jianguo， et al. Setting of excitation frequency of SI-FLAT shapometer[J]. Journal of Tianjin University， 2014，47(10):871-878. (in Chinese)

[8] 劉曉英.板帶平直度與其他因素的關(guān)系[J].輕合金加工技術(shù)，1998，26(11)：16-17.

Liu Xiaoying. Relationship between strip flatness and other factors[J]. Light Alloy Fabrication Technology， 1998，26(11):16-17. (in Chinese)

[9] 孫一康.帶鋼冷連軋計算機(jī)控制[M].北京：冶金工業(yè)出版社，2002.

Sun Yikang. Computer control for tandem cold rolling of strip[M]. Beijing: Metallurgical Industry Press， 2002. (in Chinese)

(責(zé)任編輯：孫竹鳳)

Effect of Different Stands of 2 180 mm Tandem Cold Rolling Mill on Shape Control of Product

YANG Guang-hui， ZHANG Jie， CAO Jian-guo， LI Hong-bo， HUANG Qiao-bao

(School of Mechanical Engineering， University of Science and Technology Beijing， Beijing 100083， China)

In order to analyze the mutual influences among different stands of tandem cold rolling mill and the influence of each stand on the final product shape, a control efficacy model was built with the influence coefficient method based on the crown and flatness models of continuous cold rolling. Aiming at a 2 180 mm tandem cold rolling mill, combined with its four typical rolling schedules, the shape control efficacies of different stands were studied, and it was found that, for the same rolling schedule, the ability of each stand to affect the product shape varied a lot, and the closer to the last stand, the greater effect on product shape. Among the five stands of tandem cold rolling mill, except for the last stand, i.e, the fifth stand, the fourth stand had the largest effect on the product shape, the maximum value approached to the 10 IU(1 IU=10-5). The influence coefficient difference of the shape adjustment means of different stands on final product shape between the stands was about 10-3order of magnitude. Finally, the shape data were counted and analyzed and the rolling trial was carried out to prove the reliability of the conclusion.

cold rolling mill; strip; product shape; effect coefficient; stand

2015-08-24

北京高等學(xué)校青年英才計劃資助項(xiàng)目(YETP0369)；中央高?；究蒲袠I(yè)務(wù)費(fèi)基金資助項(xiàng)目(FRF-BR-15-047A)

楊光輝(1977—)，男，博士，副教授，E-mail：yanggh@ustb.edu.cn.

TG 333.71

A

1001-0645(2016)11-1111-06

10.15918/j.tbit1001-0645.2016.11.003