• 
    

    
    

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

      側(cè)風(fēng)對(duì)某發(fā)動(dòng)機(jī)工作參數(shù)穩(wěn)定性影響的監(jiān)控及分析

      2013-11-19 08:40:10朱彥偉曹高峰劉旭東
      航空發(fā)動(dòng)機(jī) 2013年2期
      關(guān)鍵詞:壓氣機(jī)風(fēng)向標(biāo)準(zhǔn)差

      朱彥偉,曹高峰,劉旭東

      (中航工業(yè)飛行試驗(yàn)研究院,西安710089)

      0 引言

      進(jìn)氣入口氣體流場(chǎng)的不均勻變化將影響發(fā)動(dòng)機(jī)的穩(wěn)定工作裕度[1],進(jìn)氣道與發(fā)動(dòng)機(jī)的相容性是國(guó)內(nèi)外一直在研究的重要課題,文獻(xiàn)[2]根據(jù)平行壓氣機(jī)理論計(jì)算分析了進(jìn)氣條件對(duì)單軸發(fā)動(dòng)機(jī)工作穩(wěn)定性的影響,結(jié)果表明進(jìn)氣條件對(duì)發(fā)動(dòng)機(jī)的工作穩(wěn)定性影響很顯著。對(duì)于雙軸發(fā)動(dòng)機(jī)來說,目前還沒有完善而實(shí)用的穩(wěn)定性分析理論方法,因此,對(duì)進(jìn)氣道與發(fā)動(dòng)機(jī)相容性評(píng)定的主要手段還是只能依靠試驗(yàn)[3-4]。在國(guó)內(nèi)一般是采用可移動(dòng)插板式畸變發(fā)生器逼喘的方法來模擬進(jìn)氣畸變,通過對(duì)發(fā)動(dòng)機(jī)進(jìn)口流場(chǎng)和工作參數(shù)的變化進(jìn)行分析的方法來評(píng)估進(jìn)氣畸變對(duì)發(fā)動(dòng)機(jī)工作穩(wěn)定性的影響[5]。進(jìn)氣畸變直接影響發(fā)動(dòng)機(jī)安全工作,為了有效地預(yù)測(cè)發(fā)動(dòng)機(jī)在畸變進(jìn)氣時(shí)對(duì)發(fā)動(dòng)機(jī)性能產(chǎn)生的影響,文獻(xiàn)[6-8]研究了進(jìn)氣畸變對(duì)渦扇發(fā)動(dòng)機(jī)工作穩(wěn)定的影響。而發(fā)動(dòng)機(jī)側(cè)風(fēng)地面試驗(yàn)是進(jìn)氣畸變項(xiàng)目中非常重要的適航驗(yàn)證科目。主要驗(yàn)證在規(guī)定的自然風(fēng)條件下,發(fā)動(dòng)機(jī)的工作穩(wěn)定性是否受到影響,發(fā)動(dòng)機(jī)的排氣溫度、高低壓轉(zhuǎn)速、振動(dòng)值等工作參數(shù)是否超限。通過地面?zhèn)蕊L(fēng)試驗(yàn)可以確定發(fā)動(dòng)機(jī)在地面?zhèn)蕊L(fēng)條件下的使用包線,為編寫飛行使用手冊(cè)提供支持。

      在國(guó)外已經(jīng)進(jìn)行了大量的工作研究地面風(fēng)速、風(fēng)向?qū)Πl(fā)動(dòng)機(jī)工作穩(wěn)定性的影響,俄羅斯相關(guān)研究表明,從風(fēng)對(duì)發(fā)動(dòng)機(jī)工作穩(wěn)定性影響的觀點(diǎn)來看,起作用的不是某個(gè)期限內(nèi)的風(fēng)速的平均值,而是風(fēng)速的現(xiàn)有值,即在接近地表層時(shí)空氣質(zhì)量紊流運(yùn)動(dòng)的速度[9]。美國(guó)GE公司在許多年前就進(jìn)行航空民用發(fā)動(dòng)機(jī)的固定風(fēng)源側(cè)風(fēng)試驗(yàn),并建立了1套比較完善的方法體系。但是,到目前為止,在自然風(fēng)條件下由FAA和中國(guó)適航當(dāng)局目擊的大涵道比民用發(fā)動(dòng)機(jī)地面?zhèn)蕊L(fēng)試驗(yàn),在世界發(fā)動(dòng)機(jī)試驗(yàn)史上是第1次,此次試驗(yàn)可以積累很多寶貴的經(jīng)驗(yàn),為今后其他型號(hào)發(fā)動(dòng)機(jī)側(cè)風(fēng)試驗(yàn)提供技術(shù)支持。本文根據(jù)側(cè)風(fēng)試驗(yàn)過程中的實(shí)時(shí)監(jiān)控情況對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行了分析,通過發(fā)動(dòng)機(jī)工作參數(shù)標(biāo)準(zhǔn)差對(duì)比方法和計(jì)算發(fā)動(dòng)機(jī)低壓壓氣機(jī)出口特征流量系數(shù)的方法,確定不同側(cè)風(fēng)風(fēng)向條件對(duì)發(fā)動(dòng)機(jī)工作穩(wěn)定性的影響程度;結(jié)合對(duì)比分析結(jié)果提出合理增加試驗(yàn)點(diǎn)的建議,為今后側(cè)風(fēng)試驗(yàn)的高效進(jìn)行提供方法支持。

      1 試驗(yàn)方法和條件

      發(fā)動(dòng)機(jī)側(cè)風(fēng)地面試驗(yàn)是進(jìn)氣畸變項(xiàng)目中非常重要的適航驗(yàn)證科目。主要驗(yàn)證在規(guī)定的自然風(fēng)條件下,發(fā)動(dòng)機(jī)的工作穩(wěn)定性是否受到影響,發(fā)動(dòng)機(jī)的排氣溫度、高低壓轉(zhuǎn)速等工作參數(shù)是否超限。

      本次側(cè)風(fēng)試驗(yàn)發(fā)動(dòng)機(jī)為美國(guó)GE公司生產(chǎn)的某型大涵道比發(fā)動(dòng)機(jī),發(fā)動(dòng)機(jī)的結(jié)構(gòu)外形如圖1所示,發(fā)動(dòng)機(jī)工作時(shí),空氣經(jīng)進(jìn)氣道進(jìn)入風(fēng)扇增壓的氣流分為2路,其中大部分的氣流進(jìn)入外涵道,其余部分空氣經(jīng)壓氣機(jī)進(jìn)口導(dǎo)流葉片進(jìn)入壓氣機(jī)、燃燒室和內(nèi)涵噴管。發(fā)動(dòng)機(jī)采用尾吊式安裝方式在飛機(jī)上進(jìn)行安裝,該發(fā)動(dòng)機(jī)在飛機(jī)上的安裝位置和側(cè)風(fēng)風(fēng)向的規(guī)定方法如圖2所示。

      圖1 發(fā)動(dòng)機(jī)外形結(jié)構(gòu)

      試驗(yàn)開始前利用經(jīng)過標(biāo)定后的氣象車測(cè)量距地面高10m處的風(fēng)速風(fēng)向,如果自然風(fēng)風(fēng)速滿足試驗(yàn)條件以后,將配裝試驗(yàn)發(fā)動(dòng)機(jī)的飛機(jī)拖至地面試驗(yàn)區(qū)并調(diào)整飛機(jī)的擺放方向。測(cè)試工程師檢查測(cè)試系統(tǒng)工作正常后,完成試驗(yàn)發(fā)動(dòng)機(jī)開車并暖機(jī)后,將油門桿推至最大起飛位置使發(fā)動(dòng)機(jī)在最大功率狀態(tài)工作3min,采用專用的設(shè)備對(duì)發(fā)動(dòng)機(jī)的高低壓轉(zhuǎn)速、高壓壓氣機(jī)出口靜壓、發(fā)動(dòng)機(jī)排氣溫度等參數(shù)進(jìn)行記錄,同時(shí)實(shí)時(shí)測(cè)量自然風(fēng)風(fēng)速和風(fēng)向的變化。第1個(gè)試驗(yàn)點(diǎn)完成以后,調(diào)整飛機(jī)的擺放位置進(jìn)行下一個(gè)試驗(yàn)點(diǎn)試驗(yàn)。

      圖2 側(cè)風(fēng)風(fēng)向規(guī)定

      2 數(shù)據(jù)分析方法

      2.1 發(fā)動(dòng)機(jī)工作參數(shù)計(jì)算

      發(fā)動(dòng)機(jī)工作參數(shù)的波動(dòng)是其不能穩(wěn)定工作的典型特征,發(fā)動(dòng)機(jī)工作參數(shù)的變化能夠直接反映其是否工作正常,通常監(jiān)控的發(fā)動(dòng)機(jī)工作參數(shù)有:低壓轉(zhuǎn)子轉(zhuǎn)速N1,高壓轉(zhuǎn)子轉(zhuǎn)速N2,高壓壓氣機(jī)出口靜壓PS3和發(fā)動(dòng)機(jī)排氣溫度TI[10]。通過計(jì)算對(duì)比不同風(fēng)向條件下工作參數(shù)的標(biāo)準(zhǔn)差,來反映發(fā)動(dòng)機(jī)工作參數(shù)的波動(dòng)情況。

      在給定風(fēng)速風(fēng)向條件下,當(dāng)發(fā)動(dòng)機(jī)達(dá)到給定功率狀態(tài)后,錄取一定時(shí)間(選取20s)的發(fā)動(dòng)機(jī)工作參數(shù),假定試驗(yàn)點(diǎn)樣本數(shù)為N,其統(tǒng)計(jì)方法如下:

      (1)將N個(gè)數(shù)據(jù)樣本點(diǎn)排列成數(shù)據(jù)序列;

      (2)求出數(shù)據(jù)序列的平均值和標(biāo)準(zhǔn)方差。

      其中,參數(shù)標(biāo)準(zhǔn)差S 的大小直接反映了數(shù)據(jù)參數(shù)離散程度,在給定發(fā)動(dòng)機(jī)狀態(tài)下(發(fā)動(dòng)機(jī)油門桿角度不變),工作參數(shù)的標(biāo)準(zhǔn)差S 越大,表明該參數(shù)的波動(dòng)量越大,則說明發(fā)動(dòng)機(jī)的工作穩(wěn)定性越差。

      2.2 特征流量系數(shù)計(jì)算

      根據(jù)氣體動(dòng)力學(xué)公式可知,可以用某一點(diǎn)的總壓和靜壓的比值,來表征壓氣機(jī)空氣流量系數(shù)的變化趨勢(shì),利用低壓壓氣機(jī)出口(25截面)測(cè)量耙上某點(diǎn)的總靜壓比(Pt25/Ps25)與該截面處發(fā)動(dòng)機(jī)的流量系數(shù)q(λ)存在正比函數(shù)關(guān)系[11]

      定義Q 為表征流量系數(shù),其關(guān)系式為

      Q值的變化大小可以反映低壓壓氣機(jī)出口(高壓壓氣機(jī)進(jìn)口)截面處氣流軸向速度的變化大小,因此,可以利用Q值的變化情況來反映該截面處氣流速度變化,進(jìn)而來反映發(fā)動(dòng)機(jī)的工作穩(wěn)定性。

      3 試驗(yàn)結(jié)果分析

      根據(jù)發(fā)動(dòng)機(jī)工作特性可知,隨著發(fā)動(dòng)機(jī)功率狀態(tài)的增加,發(fā)動(dòng)機(jī)各部件工作環(huán)境變得越惡劣,因此,以發(fā)動(dòng)機(jī)在最大功率狀態(tài)情況下的工作參數(shù)錄取數(shù)據(jù)為依據(jù),來分析側(cè)風(fēng)對(duì)發(fā)動(dòng)機(jī)工作穩(wěn)定性的影響。

      在不同風(fēng)速風(fēng)向情況下,利用式(1)和(2)對(duì)發(fā)動(dòng)機(jī)在最大工作狀態(tài)的20s工作參數(shù)統(tǒng)計(jì)見表1。從表中統(tǒng)計(jì)數(shù)據(jù)的“均值”數(shù)據(jù)列可見,風(fēng)向?qū)Πl(fā)動(dòng)機(jī)的壓氣機(jī)排放靜壓Ps3和發(fā)動(dòng)機(jī)排氣溫度TI的影響非常明顯,在180°(尾風(fēng))風(fēng)向?qū)Πl(fā)動(dòng)機(jī)排氣溫度的影響比較明顯,由于發(fā)動(dòng)機(jī)在尾風(fēng)條件下工作時(shí)將發(fā)動(dòng)機(jī)排氣吸進(jìn)進(jìn)氣道導(dǎo)致的。在230°風(fēng)向條件下發(fā)動(dòng)機(jī)工作參數(shù)的標(biāo)準(zhǔn)差要明顯的高于其他風(fēng)向條件下的工作參數(shù)的標(biāo)準(zhǔn)差,其中高壓壓氣機(jī)排放靜壓Ps3和發(fā)動(dòng)機(jī)排氣溫度TI的標(biāo)準(zhǔn)差明顯高于其他風(fēng)向條件,這2個(gè)參數(shù)的波動(dòng)最能夠反映發(fā)動(dòng)機(jī)工作情況,由此可知230°風(fēng)向?qū)Πl(fā)動(dòng)機(jī)工作穩(wěn)定性最大。

      表1 發(fā)動(dòng)機(jī)最大工作狀態(tài)的20s工作參數(shù)統(tǒng)計(jì)

      側(cè)風(fēng)風(fēng)向?qū)Πl(fā)動(dòng)機(jī)不同工作參數(shù)標(biāo)準(zhǔn)差的影響如圖3所示。從圖中明顯可見,反映發(fā)動(dòng)機(jī)工作參數(shù)穩(wěn)定性隨側(cè)風(fēng)風(fēng)向的變化趨勢(shì)。當(dāng)側(cè)風(fēng)風(fēng)向在0~90°之間變化時(shí),對(duì)發(fā)動(dòng)機(jī)工作穩(wěn)定性幾乎沒有造成影響;當(dāng)側(cè)風(fēng)風(fēng)向在90~230°之間變化時(shí),隨著風(fēng)向角度的增加,對(duì)發(fā)動(dòng)機(jī)穩(wěn)定工作影響越來越大,其中在風(fēng)向?yàn)?30°時(shí)其影響程度最大;當(dāng)側(cè)風(fēng)風(fēng)向在230~270°之間變化時(shí),隨著風(fēng)向角度的增加,對(duì)發(fā)動(dòng)機(jī)穩(wěn)定工作影響迅速變小。

      圖3 側(cè)風(fēng)風(fēng)向?qū)Πl(fā)動(dòng)機(jī)工作參數(shù)標(biāo)準(zhǔn)差的影響

      從發(fā)動(dòng)機(jī)工作參數(shù)的標(biāo)準(zhǔn)差可以很直觀地反映側(cè)風(fēng)風(fēng)向?qū)Πl(fā)動(dòng)機(jī)工作穩(wěn)定性的影響,而通過發(fā)動(dòng)機(jī)低壓壓氣機(jī)出口截面(25截面)的表征流量系數(shù),更能直觀地體現(xiàn)側(cè)風(fēng)風(fēng)向?qū)Πl(fā)動(dòng)機(jī)工作穩(wěn)定性的影響。不同風(fēng)向條件下發(fā)動(dòng)機(jī)25截面特征流量系數(shù)Q 的變化如圖4所示(圖中F、M、T分別代表不同測(cè)量耙的安裝位置)。從圖中可見,在90°風(fēng)向情況下特征流量系數(shù)Q 非常穩(wěn)定,趨于1條穩(wěn)定的直線,說明此時(shí)發(fā)動(dòng)機(jī)工作非常穩(wěn)定;在135°風(fēng)向情況下特征流量系數(shù)Q 就產(chǎn)生了輕微的波動(dòng),但整體趨勢(shì)還是1條直線,說明此時(shí)發(fā)動(dòng)機(jī)工作還能穩(wěn)定工作;在180°風(fēng)向情況下特征流量系數(shù)Q 就產(chǎn)生了較大的波動(dòng),其中QT值波動(dòng)開始變得比較不穩(wěn)定氣流,T位置區(qū)域的流速波動(dòng)要明顯高于F和M位置處的流速波動(dòng),說明此時(shí)發(fā)動(dòng)機(jī)工作開始出現(xiàn)不穩(wěn)定狀態(tài);在230°風(fēng)向情況下特征流量系數(shù)Q 就產(chǎn)生了嚴(yán)重的波動(dòng),在3個(gè)不同位置的表征流量系數(shù)(QF、QM和QT)都產(chǎn)生了非常嚴(yán)重的波動(dòng),說明此時(shí)發(fā)動(dòng)機(jī)工作已經(jīng)出現(xiàn)非常嚴(yán)重的不穩(wěn)定狀態(tài)。由于在進(jìn)行270°風(fēng)向試驗(yàn)過程中,測(cè)量耙數(shù)據(jù)沒有采集,沒有給出270°風(fēng)向條件下發(fā)動(dòng)機(jī)25截面特征流量系數(shù)Q 的變化。

      圖4 不同風(fēng)向條件下發(fā)動(dòng)機(jī)25截面特征流量系數(shù)Q 的變化

      從上面的分析可知,230°風(fēng)向側(cè)風(fēng)對(duì)發(fā)動(dòng)機(jī)的工作穩(wěn)定性影響最大,主要與發(fā)動(dòng)機(jī)的安裝形式有關(guān),對(duì)于尾吊式安裝的發(fā)動(dòng)機(jī)(試驗(yàn)發(fā)動(dòng)機(jī)ARJ21-700飛機(jī)安裝位置如圖5所示)來說,機(jī)身和垂尾就像1堵墻擋在氣流前面,當(dāng)風(fēng)速風(fēng)向時(shí)刻變化的側(cè)風(fēng)繞過機(jī)身和垂尾以后進(jìn)入發(fā)動(dòng)機(jī)進(jìn)氣道,將導(dǎo)致發(fā)動(dòng)機(jī)入口的氣流變得雜亂無章,帶有不同尺度的氣流團(tuán)進(jìn)氣發(fā)動(dòng)機(jī)壓氣機(jī),將導(dǎo)致發(fā)動(dòng)機(jī)不穩(wěn)定工作。試驗(yàn)過程中應(yīng)當(dāng)把180~270°風(fēng)向變化區(qū)域作為重點(diǎn)試驗(yàn)區(qū)域,在該區(qū)域的試驗(yàn)數(shù)據(jù)最能體現(xiàn)發(fā)動(dòng)機(jī)的抗側(cè)風(fēng)能力。

      圖5 發(fā)動(dòng)機(jī)安裝位置

      4 結(jié)論

      本文利用的發(fā)動(dòng)機(jī)工作參數(shù)標(biāo)準(zhǔn)差統(tǒng)計(jì)法和發(fā)動(dòng)機(jī)25截面表征流量系數(shù)法可以作為評(píng)定發(fā)動(dòng)機(jī)工作穩(wěn)定性的參考法,具有較高的工程使用價(jià)值,通過對(duì)不同側(cè)風(fēng)風(fēng)向條件下發(fā)動(dòng)機(jī)最大工作狀態(tài)工作參數(shù)的對(duì)比分析可以得到的結(jié)論如下:

      (1)側(cè)風(fēng)風(fēng)向在0~90°范圍內(nèi)變化時(shí),側(cè)風(fēng)對(duì)發(fā)動(dòng)機(jī)的工作穩(wěn)定性沒有影響。

      (2)側(cè)風(fēng)風(fēng)向在90~230°范圍內(nèi)變化時(shí),隨著風(fēng)向角度的增加,側(cè)風(fēng)對(duì)發(fā)動(dòng)機(jī)的工作穩(wěn)定性的影響越來越大。

      (3)側(cè)風(fēng)風(fēng)向在230~270°范圍內(nèi)變化時(shí),隨著風(fēng)向角度的增加,側(cè)風(fēng)對(duì)發(fā)動(dòng)機(jī)的工作穩(wěn)定性的影響急劇減小。

      (4)230°側(cè)風(fēng)風(fēng)向?qū)Πl(fā)動(dòng)機(jī)工作穩(wěn)定性影響最大,對(duì)于尾吊式安裝的發(fā)動(dòng)機(jī),可以將180~270°風(fēng)向變化區(qū)域作為重點(diǎn)試驗(yàn)區(qū)域,在該區(qū)域內(nèi)應(yīng)當(dāng)適當(dāng)增加試驗(yàn)點(diǎn),可以在節(jié)省試驗(yàn)成本的基礎(chǔ)上更快確定發(fā)動(dòng)機(jī)側(cè)風(fēng)條件下的工作邊界。

      [1]Society of Automotive Engineers.Gas turbine engine inlet flow distortion guidelines[R].SAE-ARP1420B.

      [2]Shen H.The prediction of performance of turbojet engine with distorted inlet flow and experimental studies [R].AIAA-83-7030.

      [3]Reid C.The response of axial flow comperssors to intake flow distortion[R].ASME 69-GT-29.

      [4]Chung K,Hosny W M,Steenken W G.Aerodynamic stability analysis of NASA J85-13 planar perssure pulse generator installation[R].NASA-CR-165141.

      [5]劉大響,葉培梁.航空發(fā)動(dòng)機(jī)手冊(cè):進(jìn)排氣裝置:第7冊(cè)[M].北京:航空工業(yè)出版社,2000:60-63.LIU Daxiang,YE Peiliang.Aviation engine manual:intake and exhaust device:senventh volume[M].Beijing:Aviation Industry Press,2000:60-63.(in Chinese)

      [6]Braithwaite W N,Dicus J H,Moss J E.Evaluation of a turofan engine of air jets as a steady-state inlet flow distortion device[R].NASA-TM-K-1955,1970.

      [7]Evans D G.Some comparisons of the flow characteristics of a turbofan compressor system with and without inlet rressure distortion[R].NASA-TM-K-71574,1974.

      [8]De Bogdan C E.Effect of a 180° extent inlet pressure distortion on the internal flow condition of a TF30-P-3 engine[R].NASA-TM-K-3267,1975.

      [9]A·П 列烏特, A·C 謝爾巴科夫.進(jìn)氣道與發(fā)動(dòng)機(jī)共同工作的特性及其飛行試驗(yàn)方法[M].劉兆林,譯.北京:國(guó)防工業(yè)出版社,1983:30-32.(in Chinese)ЛеУТA·П, Cщербаков·A.The compatibility of inlet/engine and flight test method[M].Liu zhaolin,translated.Beijing:National Defense Industry Press,1983:30-32.

      [10]莫翔宇,唐慶如,楊曉強(qiáng).探析民用大涵道比渦扇發(fā)動(dòng)機(jī)的狀態(tài)監(jiān)控技術(shù)[J].中國(guó)民航飛行學(xué)院學(xué)報(bào),2009,20(5):16-19.MO Xiangyu,TANG Qingru,YANG Xiaoqiang.Analysis of civil high bypass ratio turbofan engine condition monitoring technology[J].Journal of Civil Aviation Flight University of China,2009,20(5):16-19.(in Chinese)

      [11]王保國(guó),劉淑艷,劉艷明,等.空氣動(dòng)力學(xué)基礎(chǔ)[M].北京:國(guó)防工業(yè)出版社,2009:73-74.WANG Baoguo,LIU Shuyan,LIU Yanming,et al.Fundamentals of aerodynamics[M].Beijing:National Defense Industry Press,2009:73-74.(in Chinese)

      猜你喜歡
      壓氣機(jī)風(fēng)向標(biāo)準(zhǔn)差
      軸流壓氣機(jī)效率評(píng)定方法
      重型燃?xì)廨啓C(jī)壓氣機(jī)第一級(jí)轉(zhuǎn)子葉片斷裂分析
      用Pro-Kin Line平衡反饋訓(xùn)練儀對(duì)早期帕金森病患者進(jìn)行治療對(duì)其動(dòng)態(tài)平衡功能的影響
      壓氣機(jī)緊湊S形過渡段內(nèi)周向彎靜子性能數(shù)值計(jì)算
      自然與風(fēng)Feeling Nature
      行業(yè)統(tǒng)計(jì)帶來哪些風(fēng)向?
      對(duì)于平均差與標(biāo)準(zhǔn)差的數(shù)學(xué)關(guān)系和應(yīng)用價(jià)值比較研究
      風(fēng)向
      風(fēng)能(2015年8期)2015-02-27 10:15:11
      風(fēng)向
      風(fēng)能(2015年4期)2015-02-27 10:14:30
      高壓比離心壓氣機(jī)設(shè)計(jì)及試驗(yàn)驗(yàn)證
      南乐县| 万州区| 南和县| 海原县| 昭觉县| 子长县| 南投县| 屯昌县| 宁河县| 腾冲县| 霍林郭勒市| 白城市| 资兴市| 如皋市| 江油市| 额尔古纳市| 许昌市| 徐水县| 宁津县| 景德镇市| 仁怀市| 大同市| 临潭县| 于田县| 嵊泗县| 义乌市| 海盐县| 江阴市| 乐至县| 张家港市| 贡山| 丹巴县| 常宁市| 周至县| 靖州| 永安市| 庆云县| 平顶山市| 临夏县| 侯马市| 河东区|