• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      飛鋸控制系統(tǒng)鋸車位置和速度的測量方法

      2016-12-06 08:02:43甄自源楊振強(qiáng)
      關(guān)鍵詞:誤差率管材編碼器

      甄自源, 楊振強(qiáng)

      (大連理工大學(xué) 電氣工程學(xué)院,遼寧 大連 116023)

      ?

      飛鋸控制系統(tǒng)鋸車位置和速度的測量方法

      甄自源, 楊振強(qiáng)

      (大連理工大學(xué) 電氣工程學(xué)院,遼寧 大連 116023)

      飛鋸控制系統(tǒng)作為一種高精度管材切割系統(tǒng),其精度和穩(wěn)定性很大程度上取決于鋸車位置和速度的測量精度和穩(wěn)定性。基于以上需求,使用TMS320F28069處理增量式光電編碼器發(fā)出的脈沖信號,提出了一種測量該系統(tǒng)鋸車位置與速度的方法。該方法引入脈沖長度簡化運(yùn)算和操作,提高了更換設(shè)備時的效率,同時給出了鋸車行進(jìn)長度的計算公式。在測量鋸車速度方面,優(yōu)化M/T加權(quán)平均過渡法以保證測速能夠平滑過渡,并通過誤差分析和仿真驗證了此方法的可行性,滿足了系統(tǒng)對精度和穩(wěn)定性的要求。

      飛鋸控制系統(tǒng);脈沖長度;M/T加權(quán)平均過渡法

      0 引言

      飛鋸控制系統(tǒng)可廣泛用于管材、型材等生產(chǎn)線的高精度連續(xù)自動追蹤切割,大幅提高生產(chǎn)效率。由于系統(tǒng)為動態(tài)在線切割,對位置與速度的精度要求較高[1],因此能夠準(zhǔn)確測量鋸車的位置和轉(zhuǎn)速顯得至關(guān)重要。該系統(tǒng)以DSP作為運(yùn)動控制器,其中的eQEP(增強(qiáng)正交編碼脈沖)模塊能夠高速采集光電編碼器的脈沖信號,再對其進(jìn)行處理和計算即可測量鋸車的位置和轉(zhuǎn)速,精確度與PLC控制相比大為提高。同時在計算轉(zhuǎn)速時使用更加合理的算法,可以使高轉(zhuǎn)速和低轉(zhuǎn)速之間能夠平滑過渡,進(jìn)一步提高系統(tǒng)的穩(wěn)定性。

      1 飛鋸控制系統(tǒng)概述

      1.1 飛鋸控制系統(tǒng)的組成

      本文所述飛鋸控制系統(tǒng)由以TMS320F28069處理器為核心的運(yùn)動控制器、VFD-VE伺服型變頻器、DOP-B05S100 觸摸屏、鋸車系統(tǒng)以及測速系統(tǒng)等組成,如圖1所示。

      圖1 飛鋸控制系統(tǒng)組成框圖

      運(yùn)動控制器對整個系統(tǒng)進(jìn)行監(jiān)視和控制。伺服型變頻器接收運(yùn)動控制器發(fā)出的電壓信號轉(zhuǎn)化為頻率來調(diào)整鋸車電機(jī)的轉(zhuǎn)速。觸摸屏可發(fā)出不同的動作指令,調(diào)整系統(tǒng)的部分參數(shù),并顯示系統(tǒng)的運(yùn)行狀態(tài)。鋸車系統(tǒng)由可前后移動的鋸車與固定軌道組成,鋸車上裝有鋸片,并連接繼電系統(tǒng)[2],其機(jī)械結(jié)構(gòu)示意圖如圖2所示[3-4]。繼電系統(tǒng)用來下達(dá)夾緊、松夾、落鋸、抬鋸等動作指令。測速系統(tǒng)可實(shí)時測量各電機(jī)的轉(zhuǎn)速,其中編碼器選用的規(guī)格是24V,1000線的增量式光電編碼器。

      圖2 鋸車系統(tǒng)機(jī)械結(jié)構(gòu)示意圖

      1.2 飛鋸控制系統(tǒng)鋸車的運(yùn)動狀態(tài)分析

      在本系統(tǒng)中,需要測量以下四部分組件的速度:鋸車行進(jìn)速度、管行進(jìn)速度、鋸片轉(zhuǎn)速和鋸片主軸移動速度。鋸車和管材的各階段運(yùn)動狀態(tài)圖如圖3所示,其中T1為鋸車加速階段,T2為夾緊前鋸車調(diào)整階段,T3、T4、T5、T6分別為夾緊、落鋸、抬鋸、松夾階段,T7為鋸車回零點(diǎn)階段,T8為鋸車零點(diǎn)調(diào)整階段[5]。

      圖3 鋸車各階段運(yùn)動狀態(tài)

      理想狀態(tài)下,管材勻速前進(jìn),鋸車與管行進(jìn)至兩者行程差為設(shè)定的切割管長時,兩者恰好同速,輸出夾緊管材的信號再落鋸切管;然而由于實(shí)際上管材并不總是勻速前進(jìn),且測量有誤差,此時鋸車與管可能不同速,但速度差最好盡量小。由此可見,鋸車行進(jìn)位置和速度的測量精度及穩(wěn)定性尤為重要。

      本系統(tǒng)管長設(shè)定范圍為2m~6m,定尺精度為3mm,即切割后的管材與設(shè)定管材的誤差最大值為±3mm。具體到各運(yùn)動階段來說,發(fā)送夾緊管材信號前定尺精度主要取決于管材與鋸車相對位置的測量精度;在發(fā)送夾緊管材信號后至管材被夾緊(由于繼電器有動作時間)時,定尺精度主要取決于管材與鋸車相對速度的測量精度,此部分定尺精度系統(tǒng)要求為1mm。本文以鋸車為例,提出其位置和速度的測量方法,又由于管材以勻速行進(jìn),則定尺精度以及系統(tǒng)穩(wěn)定性主要取決于鋸車位置和速度的精度和穩(wěn)定性。

      2 eQEP(增強(qiáng)正交編碼脈沖)模塊

      TMS320F28069處理器的eQEP模塊結(jié)合增量式光電編碼器可實(shí)現(xiàn)獲取運(yùn)動控制系統(tǒng)位置、速度以及方向信息的功能[6]。

      eQEP模塊的QEPA和QEPB引腳為正交時鐘模式輸入引腳,可接收編碼器輸出的相差90°的方波信號。如圖4所示,這一系列正交脈沖信號經(jīng)由正交解碼單元QDU解碼后得到4倍頻的位置脈沖信號和方向信號,再通過位置計數(shù)器進(jìn)行脈沖計數(shù)得到QPOSCNT和QDIR的值,即為電機(jī)的位置信息和方向信息[7]。

      圖4 eQEP輸入信號計數(shù)與方向解碼

      3 行進(jìn)長度的測量

      3.1 編碼器位置變化的脈沖計數(shù)測量

      本文主要關(guān)注的是鋸車的行進(jìn)長度。通過將光電編碼器與電機(jī)相連,根據(jù)1.1節(jié)和1.2節(jié)所述就可對鋸車位置進(jìn)行測量,由于行進(jìn)長度是和編碼器的累積行程相關(guān),而不用關(guān)注編碼器的旋轉(zhuǎn)周數(shù),因此只需記錄eQEP位置鎖存器QPOSSLAT的值來確定編碼器向前或向后行進(jìn)了多少個QCLK脈沖。

      3.2 脈沖長度的應(yīng)用

      若要計算鋸車行進(jìn)長度,除了記錄編碼器產(chǎn)生的脈沖計數(shù)值,還需要引入一個系數(shù)——脈沖長度,即編碼器的一個脈沖可換算為多少長度。由于鋸車系統(tǒng)結(jié)構(gòu)不同,因此鋸車脈沖長度Kc的計算方法也不同,如式(1)所示:

      (1)

      其中,Z為齒輪齒數(shù),M為齒輪模數(shù),X為編碼器計數(shù)的倍頻,P為編碼器的每轉(zhuǎn)脈沖數(shù),I為減速箱傳動比。本系統(tǒng)使用的76型飛鋸齒數(shù)Z為25,模數(shù)M為8,鋸車履帶傳動比I為7.1,鋸車編碼器為1000/rpm,4倍頻,則代入上面公式計算得鋸車脈沖長度Kc=0.0221mm。

      由于鋸車脈沖長度Kc均可通過觸摸屏輸入來隨時更改,因此當(dāng)系統(tǒng)的部分設(shè)備更換或有一定的磨損后,只需更改這個脈沖長度即可滿足設(shè)備變化的要求,簡化了計算和操作。

      3.3 行進(jìn)長度的計算

      根據(jù)3.1和3.2節(jié),可得到行進(jìn)長度的計算公式,如式(2)所示:

      L=CNT×K

      (2)

      其中,L為行進(jìn)長度(前進(jìn)方向為正),CNT為eQEP位置鎖存器QPOSSLAT的值,K為脈沖系數(shù)。式(2)也可應(yīng)用在其他組件計算行進(jìn)長度的情況。

      由于每次測量行進(jìn)長度時,脈沖計數(shù)值有±1個脈沖的最大誤差,因此本方法在測量行進(jìn)長度時的誤差為1個脈沖對應(yīng)的長度,即脈沖長度。同時為了避免系統(tǒng)長時間運(yùn)行導(dǎo)致誤差累積,鋸車在完成一定次數(shù)切管后,返回起始位置時應(yīng)對起始位置,即行進(jìn)長度為0的點(diǎn)進(jìn)行校正[8]。校正后的起始位置由限位開關(guān)確定。

      4 鋸車速度的測量

      4.1 M/T加權(quán)平均過渡法測速的原理

      M法和T法是電機(jī)轉(zhuǎn)速測量的常用方法。其中,M法是指通過測量在一定時間間隔T1內(nèi)編碼器產(chǎn)生的脈沖數(shù)M1來計算轉(zhuǎn)速的方法[9]。又由于本系統(tǒng)引入了脈沖長度,可將測量轉(zhuǎn)速轉(zhuǎn)化為測量脈沖頻率f1(每秒產(chǎn)生的脈沖數(shù)),故可得式(3):

      (3)

      若使用DSP的eQEP模塊進(jìn)行M法測速,則需使能eQEP的UTO(單位超時)中斷,通過設(shè)置單位周期寄存器QUPRD的值N1和系統(tǒng)時鐘頻率f0來確定T1,通過測量前后兩次進(jìn)中斷時的位置計數(shù)鎖存器QPOSLAT的差來確定M1。故式(3)可化為式(4):

      (4)

      在某些極端情況下M1會有±1的誤差,則在測量較高轉(zhuǎn)速時M1的相對誤差較小,而T1又為定值,因此M法在測量高速時誤差較小[10]。

      而T法是指通過測量編碼器產(chǎn)生一定數(shù)量的脈沖M2所需要的時間T2來計算轉(zhuǎn)速的方法。同M法,脈沖頻率f2的計算公式如式(5)[9]:

      (5)

      若使用DSP的eQEP模塊進(jìn)行T法測速,則需檢測eQEP的狀態(tài)寄存器QEPSTS的單位位置時間標(biāo)志位UPEVNT,通過設(shè)置捕捉控制寄存器QCAPCTL和單位位置事件預(yù)定標(biāo)UPPS來確定M2,而T2通過對分頻后的系統(tǒng)時鐘脈沖進(jìn)行計數(shù)的計數(shù)值N2來確定,其中時鐘頻率由系統(tǒng)時鐘頻率f0和捕捉定時器時鐘預(yù)定標(biāo)CCPS來確定,計數(shù)值由捕捉周期鎖存器QCPRDLAT來確定。則式(5)可化為式(6):

      (6)

      需要注意的是,eQEP的捕捉定時器發(fā)生上溢會導(dǎo)致鎖存器QCPRDLAT的偏差較大,因此需要監(jiān)測狀態(tài)寄存器QEPSTS捕捉上溢標(biāo)志位COEF,當(dāng)該標(biāo)志位為“1”時T2的值不能用于計算轉(zhuǎn)速。

      在某些極端情況下T2會有±1的誤差,則在測量較低轉(zhuǎn)速時T2的相對誤差較小,而M2又為定值,因此T法在測量低速時誤差較小[11]。

      M法和T法分別適用于測量高速和低速的情況,若將轉(zhuǎn)速以某一個值作為分界線來區(qū)分高速與低速,可能會出現(xiàn)M法和T法相互轉(zhuǎn)換時轉(zhuǎn)速突變的情況。根據(jù)這一情況,可采用M/T加權(quán)平均過渡法來進(jìn)行測速,即在高速和低速時分別使用M法和T法測速,在測速誤差接近需要切換時使用變系數(shù)加權(quán)平均法,使轉(zhuǎn)速平滑切換,測量誤差無波動,以提高系統(tǒng)的穩(wěn)定性[12]。

      M/T加權(quán)平均過渡法將電機(jī)轉(zhuǎn)速分為三個區(qū)域,即高速區(qū)、低速區(qū)和過渡區(qū),計算公式如式(7)所示:

      f=sf1+(1-s)f2

      (7)

      其中,f為加權(quán)計算后的脈沖頻率,與鋸車脈沖長度相乘即為鋸車行進(jìn)速度,故后文均將f的誤差率作為鋸車速度的誤差率;f1為M法測量的轉(zhuǎn)速,s為其權(quán)值;f2為T法測量的轉(zhuǎn)速,1-s為其權(quán)值。

      s的計算方法有多種,文獻(xiàn)[12]曾給出一種線性算法?;谠摲椒?,可應(yīng)用于飛鋸控制系統(tǒng)的s的計算公式如式(8)所示:

      (8)

      其中,fa是f1與f2的平均值,作為該時刻轉(zhuǎn)速處在何區(qū)域的判斷頻率值,即fa=(f1+f2)/2;fH是高速區(qū)與過渡區(qū)的分界限轉(zhuǎn)速值;fL是低速區(qū)與過渡區(qū)的分界限轉(zhuǎn)速值。

      該方法有一定的合理性,且計算簡單,但還可以進(jìn)一步優(yōu)化,具體分析將在4.2節(jié)論述。

      4.2 誤差分析與仿真

      由1.2節(jié)可知,在發(fā)送夾緊管材信號后至管材被夾緊時,定尺精度系統(tǒng)要求為1mm,定尺精度以及系統(tǒng)穩(wěn)定性主要取決于鋸車位置和速度的精度和穩(wěn)定性。設(shè)ΔL為鋸車行進(jìn)長度測量誤差最大值,δmax為測速誤差率最大值,fmax為脈沖頻率最大值,ΔT為夾緊繼電器最大動作時間,則其計算公式如式(9)所示:

      ΔL=δmaxfmaxKcΔT

      (9)

      其中,fmax為200kHz,ΔT為25ms。由3.2節(jié)Kc=0.0221可計算得,δmax=0.9%。故測速誤差率不能超過0.9%。

      為了使過渡區(qū)的轉(zhuǎn)速平滑切換,測速誤差不大于0.9%,則需保證過渡區(qū)的轉(zhuǎn)速誤差率函數(shù)連續(xù),且函數(shù)值不大于0.9%。M法和T法的測速誤差率分別如式(10)和式(11)所示[13]:

      (10)

      (11)

      而過渡區(qū)的測速誤差率如式(12)所示:

      (12)

      由于M法和T法的測速誤差率函數(shù)都是單調(diào)的,因此同一轉(zhuǎn)速下測速誤差率也相同時的脈沖頻率為測速誤差分界點(diǎn),即脈沖頻率大于此頻率時M法測速誤差較小,而小于此頻率時T法測速誤差較小,故這個測速誤差分界點(diǎn)應(yīng)該在過渡區(qū)內(nèi)。由式(13)即可算出該脈沖頻率:

      (13)

      在本文的飛鋸控制系統(tǒng)中,相關(guān)參數(shù)設(shè)置如下:f0=8MHz,N1=16000,UPPS=5,CCPS=4,計算可得測速誤差分界點(diǎn)時的脈沖頻率為89.2kHz。選取fH=100kHz,fL=80kHz,此時,若使用式(8)的線性算法s的取值計算公式如式(14)所示,相應(yīng)函數(shù)圖形如圖5所示。

      (14)

      圖5 權(quán)值線性算法函數(shù)圖

      高速區(qū)(M法)的誤差率在脈沖頻率為100kHz時為最大值0.5%,低速區(qū)(T法)的誤差率在脈沖頻率為80kHz時為最大值0.503%,而在過渡區(qū),令fa=f,將式(14)代入式(12)可算出在該區(qū)域內(nèi)以脈沖頻率f(單位為kHz)為自變量的線性算法誤差率函數(shù)如式(15)所示:

      (15)

      該函數(shù)的Matlab仿真圖如圖7虛線所示,該函數(shù)在(80,100)的區(qū)間內(nèi)連續(xù),且最大值約為0.56%。而本文提出M/T加權(quán)平均過渡法權(quán)值的另一種算法,s隨fa以二次函數(shù)規(guī)律變化,簡稱為二次算法。該函數(shù)如式(16)所示,相應(yīng)函數(shù)圖如圖6所示。

      (16)

      該方法在高速區(qū)和低速區(qū)的誤差率與線性算法相同。令fa=f,將式(14)代入式(12)可得過渡區(qū)內(nèi)以脈沖頻率f為自變量的二次算法誤差率函數(shù),如式(17)所示:

      (17)

      二次算法的誤差率函數(shù)仿真圖如圖6實(shí)線所示,圖中脈沖頻率處于(80,100)的區(qū)間為過渡區(qū),

      圖6 權(quán)值二次算法函數(shù)圖

      其他區(qū)域分別對應(yīng)低速區(qū)和高速區(qū),即分別應(yīng)用了T法和M法。由該圖易知,該函數(shù)在過渡區(qū)內(nèi)時連續(xù),且最大值與一次算法的最大值相同。兩種算法的結(jié)果最大的不同是,在除邊界值和最大值處之外,二次算法的轉(zhuǎn)速誤差率均比線性算法的誤差率低,并且在分界點(diǎn)80kHz和100kHz處,二次算法函數(shù)斜率與其他區(qū)域斜率連續(xù),即轉(zhuǎn)速誤差率變化更加穩(wěn)定,因此二次算法更加合理。

      圖7 過渡區(qū)內(nèi)誤差率仿真圖

      選用二次算法的M/T加權(quán)平均過渡法測速誤差率如表1所示。由表1可知,本方法測速誤差率最大值不超過0.6%,滿足誤差率連續(xù)且不大于0.9%的要求。

      表1 M/T加權(quán)平均過渡法的頻率分區(qū)和誤差率表

      5 結(jié)論

      本文介紹了飛鋸控制系統(tǒng),以及測量其鋸車位置和轉(zhuǎn)速的精度和穩(wěn)定度的重要性。同時介紹了DSP的eQEP(增強(qiáng)正交編碼脈沖)模塊測量該系統(tǒng)鋸車行進(jìn)長度和轉(zhuǎn)速的的方法和計算公式。其中,引入脈沖長度可優(yōu)化行進(jìn)長度和轉(zhuǎn)速的算法,簡化設(shè)備更換或磨損后的操作,提高了效率;使用M/T加權(quán)平均過渡法來測速保證了高轉(zhuǎn)速與低轉(zhuǎn)速能平滑轉(zhuǎn)換,且誤差率不大于0.9%,同時對權(quán)值的計算方法進(jìn)行了優(yōu)化,使用更加合理的二次算法,進(jìn)一步保證了系統(tǒng)的精度和穩(wěn)定性。

      [1] 王海鵬,陳和平,朱英韜,等. 工業(yè)機(jī)器人在板材切割生產(chǎn)線的應(yīng)用[J]. 組合機(jī)床與自動化加工技術(shù), 2015(5):141-144.

      [2] 李平,楊曉明,牛驍. 飛鋸機(jī)的控制方案簡析[J]. 重型機(jī)械科技,2006(2):36-39.

      [3] 于冬,顧培亮,高璇. 基于PLC技術(shù)的自動切管機(jī)控制系統(tǒng)的研究[J]. 組合機(jī)床與自動化加工技術(shù), 2004(10):7-9.

      [4] 任意,石小英,曾勝.全自動型材切割機(jī)設(shè)計與實(shí)現(xiàn)[J].電機(jī)與控制應(yīng)用,2009(9):30-33.

      [5] Yuhua Zhang, Xinhuan Song.The Motion Control of Profile Milling Flying Saw[C] // MATEC Web of Conferences 2015. 2015 7th International Conference on Mechanical and Electronics Engineering, 2015, 31:03014p1- 4.

      [6] Texas instruments. TMS320x2806x Piccolo Technical Reference Manual[S]. 2011,1:452-453.

      [7] 王子博. 編碼器四倍頻電路的單片機(jī)高速算法設(shè)計[J]. 組合機(jī)床與自動化加工技術(shù), 2007(11):73-74,78.

      [8] 楊紅波. 銅材自動定尺飛鋸機(jī)控制系統(tǒng)的研究[D].大連:大連交通大學(xué), 2010.

      [9] 鄧建,林樺. 基于DSP的絕對式光電編碼器的電機(jī)轉(zhuǎn)速測量[J].電機(jī)與控制應(yīng)用,2010(1):50-52,55.

      [10] 蘇指南. 仿形銑切飛鋸的運(yùn)動控制[D].天津:天津理工大學(xué), 2011.

      [11] 李昕奇,金勇. 基于DSP2812的電動機(jī)測速方法的研究[J]. 電氣自動化,2011(4):81-82.

      [12] 張濤,楊振強(qiáng),王曉旭. 應(yīng)用eQEP及編碼器測量電機(jī)位置與速度的方法[J]. 電氣傳動,2011(4):48-51.

      [13] 陳伯時. 電力拖動自動控制系統(tǒng)-運(yùn)動控制系統(tǒng)[M]. 3版.北京:機(jī)械工業(yè)出版社,2006.

      (編輯 李秀敏)

      Method of Measuring Saw Car’s Position and Speed in Flying Saw Control System

      ZHEN Zi-yuan,YANG Zhen-qiang

      (College of Electrical Engineering, Dalian University of Technology, Dalian Liaoning 116023, China)

      As a high precision pipe cutting system, accuracy and stability of flying saw control system depends largely on the accuracy and stability of the saw car’s position and speed. Based on the above requirements,introduced a method applied to this system to measure the saw car’s position and speed, which utilizes TMS320F28069 to process the pulse signals transmitted from incremental optical encoder. The method introduces pulse length to simplify the calculation and operation, which improves efficiency of changing instruments. The computational formula of moving length of the saw car is given. The method of M/T weighted average transition is optimized to guarantee smooth transition of measuring rotate speed, of which the practicability is proved by error analysis and simulation. The method meets the demand of accuracy and stability for the system.

      flying saw control system; pulse length; M/T weighted average transition

      1001-2265(2016)11-0100-05

      10.13462/j.cnki.mmtamt.2016.11.027

      2016-01-13

      甄自源(1990—),男,河北邢臺人,大連理工大學(xué)碩士研究生,研究方向為電機(jī)控制,(E-mail)zhen_ziyuan@163.com;通訊作者:楊振強(qiáng)(1964—),男,哈爾濱人,大連理工大學(xué)副教授,博士,研究方向為電機(jī)控制、模糊控制等,(E-mail)yangzq@dlut.edu.cn。

      TH162;TG506

      A

      猜你喜歡
      誤差率管材編碼器
      生化檢驗全程中質(zhì)量控制管理方式及應(yīng)用意義
      健康大視野(2020年1期)2020-03-02 11:33:53
      降低評吸人員單料煙感官評分誤差率探討
      基于FPGA的同步機(jī)軸角編碼器
      無線傳感器網(wǎng)絡(luò)定位算法在環(huán)境監(jiān)測中的應(yīng)用研究
      基于PRBS檢測的8B/IOB編碼器設(shè)計
      電工儀表測量中容易忽略的幾個問題
      JESD204B接口協(xié)議中的8B10B編碼器設(shè)計
      電子器件(2015年5期)2015-12-29 08:42:24
      大型管材在線切割設(shè)備的設(shè)計
      多總線式光電編碼器的設(shè)計與應(yīng)用
      SA-213T23管材焊后熱處理工藝
      阿拉善左旗| 郎溪县| 怀远县| 黄浦区| 达尔| 万安县| 万盛区| 铁力市| 桃源县| 泰顺县| 景泰县| 潞城市| 临安市| 宾川县| 岱山县| 深州市| 明水县| 司法| 丰宁| 嘉鱼县| 宁蒗| 永城市| 遂宁市| 渭南市| 思茅市| 保亭| 平武县| 永寿县| 大庆市| 报价| 桂林市| 信丰县| 开封县| 桦川县| 广宁县| 翁牛特旗| 巴林左旗| 南皮县| 托克托县| 霍邱县| 来凤县|