国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

葉片吸力面不同結(jié)構(gòu)對(duì)離心泵空化初生的影響

2022-05-04 08:17:28趙偉國(guó)亢艷東李清華薛子陽(yáng)
振動(dòng)與沖擊 2022年7期
關(guān)鍵詞:空泡揚(yáng)程離心泵

趙偉國(guó), 亢艷東, 李清華, 薛子陽(yáng)

(1.蘭州理工大學(xué) 能源與動(dòng)力工程學(xué)院,蘭州 730050;2.蘭州理工大學(xué) 甘肅省流體機(jī)械及系統(tǒng)重點(diǎn)實(shí)驗(yàn)室,蘭州 730050)

離心泵作為應(yīng)用最為廣泛的通用機(jī)械之一,在生產(chǎn)生活中發(fā)揮著十分重要的作用。然而,在離心泵運(yùn)行過(guò)程中難免會(huì)出現(xiàn)空化現(xiàn)象,從而影響離心泵的水力性能,嚴(yán)重時(shí)引發(fā)的振動(dòng)和噪聲會(huì)影響泵的穩(wěn)定運(yùn)行[1-2]。

影響離心泵流動(dòng)穩(wěn)定的主要原因是由流量變化造成的流動(dòng)不穩(wěn)定性和由壓力變化造成的空化不穩(wěn)定性[3]。離心泵在小流量下運(yùn)行時(shí),葉輪進(jìn)口會(huì)出現(xiàn)回流現(xiàn)象,使流動(dòng)變得復(fù)雜[4];在離心泵中流動(dòng)分離多發(fā)生在葉輪流道的入口,造成較大的流動(dòng)損失;在葉輪流道內(nèi),由于葉片背面與工作面存在壓差,工作面的壓力高,速度低,而背面正好相反,此壓力梯度將使流體向葉片背面流動(dòng),形成二次流并加劇流動(dòng)損失[5];陳紅勛等[6-7]設(shè)計(jì)開(kāi)發(fā)了縫隙引流離心泵葉輪并對(duì)其進(jìn)行了性能研究,結(jié)果發(fā)現(xiàn)縫隙引流葉輪可以改善低比轉(zhuǎn)速離心泵內(nèi)部流動(dòng)結(jié)構(gòu)。本文布置的結(jié)構(gòu)尺寸較小,對(duì)葉片進(jìn)口的排擠較小,對(duì)離心泵流動(dòng)不穩(wěn)定性影響較??;而且布置三種結(jié)構(gòu)后,降低了葉片吸力面的負(fù)壓梯度,對(duì)離心泵的流動(dòng)有一定的改善。

離心泵內(nèi)空化主要發(fā)生在葉片前緣且前緣空化是造成揚(yáng)程下降的起源[8]。由于前緣空化是附著在葉片上的而且在靠近壁面處會(huì)形成空泡破裂,因此具有很強(qiáng)的空化破壞能力[9]。離心泵空化的發(fā)展過(guò)程可以分為空化初生階段、準(zhǔn)穩(wěn)態(tài)空化階段和不穩(wěn)定空化階段;在空化初生階段從泵外特性來(lái)看泵的性能受空化影響較小,然而如果泵長(zhǎng)期在這個(gè)階段運(yùn)行,會(huì)對(duì)葉片進(jìn)口邊造成空蝕破壞,因此,對(duì)空化初生階段進(jìn)行深入研究有非常重要的意義[10-11]。水力機(jī)械中抑制空化的方法主要包括兩種:第一種是優(yōu)化葉輪進(jìn)口幾何參數(shù),來(lái)提高離心泵的空化性能;劉宜等[12]通過(guò)改變?nèi)~輪進(jìn)口直徑與葉片進(jìn)口安放角,對(duì)離心泵進(jìn)行了全流場(chǎng)數(shù)值模擬,結(jié)果表明,在一定范圍內(nèi),隨著離心泵葉輪進(jìn)口直徑與葉片進(jìn)口安放角的增加,離心泵效率有所提高,空化性能有所改善;第二種是利用附加結(jié)構(gòu)增大葉輪進(jìn)口處的壓力達(dá)到抑制空化的目的;在離心葉輪前加裝誘導(dǎo)輪,對(duì)流體預(yù)做功,可以增大流體在葉輪進(jìn)口處的壓力,從而提高離心泵的空化性能[13],程效銳等[14]通過(guò)改變等螺距誘導(dǎo)輪的螺距,探究其對(duì)離心泵空化性能的影響,結(jié)果表明,在一定范圍內(nèi)增大誘導(dǎo)輪螺距,可以提高泵的抗汽蝕性能;在空化發(fā)生區(qū)域注入液體,聚合物溶液或氣體,通過(guò)施加外部干擾迫使邊界層向湍流過(guò)渡來(lái)達(dá)到抑制空化的目的[15],王巍等[16]提出了水翼表面主動(dòng)射流對(duì)繞水翼周?chē)鲃?dòng)加以控制的方法,以探究回射流與流場(chǎng)空化特性的關(guān)系,結(jié)果表明,射流使水翼吸力面上的邊界層減薄,增強(qiáng)了流動(dòng)的抗逆壓梯度能力,一定程度上阻擋了回射流向水翼前緣的流動(dòng);在空化發(fā)生區(qū)域設(shè)置規(guī)則或不規(guī)則圖形的障礙物,通過(guò)提高低壓區(qū)壓力來(lái)達(dá)到抑制空化的目的[17],趙偉國(guó)等[18]提出了一種在葉片表面加障礙物的方法來(lái)抑制空化初生及發(fā)展的方法,結(jié)果表明,障礙物可以有效增大葉片近壁湍動(dòng)能,改變壓力分布,對(duì)離心泵內(nèi)各個(gè)階段空化均有抑制作用。

本文通過(guò)在離心泵葉片吸力面前緣處布置凹槽、橫向障礙物及不連續(xù)障礙物三種不同結(jié)構(gòu),基于修正的SSTk-ω湍流模型和Zwart-Gerber-Belamri空化模型,采用瞬態(tài)數(shù)值模擬方法,對(duì)比分析這三種結(jié)構(gòu)對(duì)空化初生時(shí)的離心泵葉輪內(nèi)壓力分布、空泡體積及速度分布等的影響與原因,以期得到對(duì)離心泵空化初生抑制效果最好的結(jié)構(gòu),為離心泵空化抑制提供新思路。

1 模型參數(shù)與網(wǎng)格劃分

1.1 模型設(shè)計(jì)參數(shù)

選取一臺(tái)比轉(zhuǎn)速為32的低比轉(zhuǎn)速離心泵作為研究對(duì)象,其基本參數(shù)如表1所示。

表1 泵的基本參數(shù)

通過(guò)Pro/E軟件對(duì)離心泵進(jìn)行三維全流道建模,整個(gè)計(jì)算域包括進(jìn)口延長(zhǎng)段、葉輪、前后腔體、蝸殼和出口延長(zhǎng)段五部分,如圖1所示。

圖1 計(jì)算域

根據(jù)課題組前期對(duì)在葉片表面布置障礙物位置的研究[19],將三種結(jié)構(gòu)布置于葉輪半徑的28%處,截面形狀均為1 mm×1 mm的矩形,葉輪流道寬度為b,凹槽和橫向障礙物寬度取b,三個(gè)不連續(xù)障礙物寬度與中間部分空隙寬度都取0.2b。其單個(gè)葉片模型分別如圖2(a)、2(b)、2(c)所示。

(a) 凹槽結(jié)構(gòu)模型

1.2 網(wǎng)格劃分情況

通過(guò)ICEM軟件對(duì)計(jì)算域各部分進(jìn)行網(wǎng)格劃分,其中蝸殼部分采用自適應(yīng)性很好的四面體網(wǎng)格,其他部分均采用六面體網(wǎng)格,并對(duì)葉輪近壁面邊界層進(jìn)行了加密處理,首層網(wǎng)格高度設(shè)置為0.05 mm,變化率為1.2,在邊界層內(nèi)至少設(shè)置了10個(gè)網(wǎng)格節(jié)點(diǎn)。葉輪網(wǎng)格與布置三種結(jié)構(gòu)附近網(wǎng)格如圖3所示。

(a) 葉輪網(wǎng)格

為了提高數(shù)值模擬的精度,更好地捕捉近壁面區(qū)域邊界層內(nèi)的流動(dòng)情況,須要保證該區(qū)域內(nèi)有足夠的節(jié)點(diǎn)數(shù)。通常用y+值來(lái)檢驗(yàn)距離壁面最近節(jié)點(diǎn)的位置[20],y+值表示距離壁面最近的網(wǎng)格節(jié)點(diǎn)到壁面的距離,是無(wú)量綱參數(shù),定義如下

(1)

式中:τω為壁面切應(yīng)力;ρ為流體密度;Δn為與壁面最近兩個(gè)網(wǎng)格節(jié)點(diǎn)間距離;v為運(yùn)動(dòng)黏度。

本次模擬選用的SSTk-ω模型近壁區(qū)采用k-ω模型,y+在100以?xún)?nèi)基本滿(mǎn)足k-ω湍流模型對(duì)近壁區(qū)網(wǎng)格質(zhì)量要求[21]。葉輪葉片的y+值范圍分布如圖4所示。

圖4 葉輪葉片的y+分布

為了確定數(shù)值計(jì)算所需要的合適網(wǎng)格,對(duì)原型泵進(jìn)行了網(wǎng)格無(wú)關(guān)性驗(yàn)證,如表2所示。方案1計(jì)算得到的揚(yáng)程與效率值與試驗(yàn)值偏差較大,而隨著網(wǎng)格數(shù)的增大,揚(yáng)程與效率值趨于穩(wěn)定,葉輪內(nèi)空泡體積也趨于平穩(wěn),滿(mǎn)足網(wǎng)格無(wú)關(guān)性要求[22],綜合考慮模擬精度與計(jì)算資源選擇方案3網(wǎng)格進(jìn)行計(jì)算。

表2 網(wǎng)格無(wú)關(guān)性檢查

2 數(shù)值計(jì)算模型

2.1 多相流模型

本次模擬汽液兩相流動(dòng)選擇ANSYS-CFX歐拉多相流模型中的均相模型,其控制方程如下

連續(xù)性方程

(2)

動(dòng)量方程

(3)

混合相密度為

ρm=αvρv+ρl(1-αv)

(4)

式中:ρm、ρv、ρl分別為混合相密度、汽相密度與液相密度;ui、uj、uk為速度分量;μ、μt分別為層流黏度與湍流黏度;αv為汽相體積分?jǐn)?shù);δij為克羅內(nèi)克常數(shù)。

2.2 湍流模型

湍流模型選用修正的SSTk-ω模型,此模型能更好地捕捉空泡的流動(dòng)特性[23-26]。其表達(dá)式為

(5)

Gω-Yω+Dω+Sω

(6)

(7)

(8)

式中:Gk、Gω分別為湍動(dòng)能k和耗散率ω的生成項(xiàng);Γk、Γω分別為k和ω的有效擴(kuò)散系數(shù);Yk、Yω分別為k和ω的耗散;Dω為交叉擴(kuò)散項(xiàng);Sk、Sω為用戶(hù)自定義源相;Cω為經(jīng)驗(yàn)常數(shù),取1。

2.3 空化模型

Zwart-Gerber-Belamri空化模型[27]是基于簡(jiǎn)化的Rayleigh-Plesset方程發(fā)展而來(lái)的,適用于模擬離心泵內(nèi)的非定??栈匦?,其表達(dá)式為

(9)

(10)

(11)

式中:fv為汽相質(zhì)量分?jǐn)?shù);Re、Rc分別為液相汽化速率與汽相凝結(jié)速率;Rb為空泡成核點(diǎn)半徑;Pv為飽和蒸汽壓力;αnuc為空泡成核點(diǎn)體積分?jǐn)?shù);Ce、Cc分別為汽化和凝結(jié)經(jīng)驗(yàn)系數(shù)。

3 數(shù)值模擬與試驗(yàn)驗(yàn)證

3.1 數(shù)值模擬設(shè)置

通過(guò)ANSYS CFX軟件對(duì)離心泵進(jìn)行數(shù)值計(jì)算,設(shè)置進(jìn)、出口邊界為壓力進(jìn)口與質(zhì)量流量出口,通過(guò)逐步降低離心泵進(jìn)口總壓的方式來(lái)控制泵內(nèi)空化的劇烈程度,進(jìn)行定常和非定常模擬。非定常計(jì)算一般要計(jì)算5~8個(gè)葉輪旋轉(zhuǎn)周期才能得到可靠的解[28],設(shè)置葉輪旋轉(zhuǎn)圈數(shù)為8圈,葉輪每旋轉(zhuǎn)0.5°作為一個(gè)時(shí)間步長(zhǎng),計(jì)算總步數(shù)為5 760步,收斂精度為10-5。為保證葉輪內(nèi)流體進(jìn)入穩(wěn)定流動(dòng)狀態(tài),選取最后一圈的最后一個(gè)時(shí)間步長(zhǎng)的計(jì)算結(jié)果進(jìn)行非定常分析,計(jì)算結(jié)果的葉輪與蝸殼隔舌的相對(duì)位置是相同的。

3.2 試驗(yàn)裝置

為了驗(yàn)證數(shù)值模擬的可靠性,在離心泵閉式試驗(yàn)臺(tái)上進(jìn)行了外特性試驗(yàn)。試驗(yàn)裝置如圖5所示。

1-電動(dòng)機(jī);2-聯(lián)軸器;3-扭矩儀;4-變頻器;5-穩(wěn)流罐;6-壓力脈動(dòng)監(jiān)測(cè)裝置;7-出口壓力表;8-進(jìn)口壓力表;9-可視化離心泵;10-電磁流量計(jì);11-出口調(diào)節(jié)閥門(mén);12-真空表;13-汽蝕罐;14-液環(huán)真空泵

3.3 外特性試驗(yàn)驗(yàn)證

圖6(a)和圖6(b)分別為原型泵和布置三種結(jié)構(gòu)離心泵的流量揚(yáng)程對(duì)比曲線(xiàn)與流量效率對(duì)比曲線(xiàn)。從圖6中可以看出,原型泵的揚(yáng)程模擬曲線(xiàn)與試驗(yàn)曲線(xiàn)有相同變化趨勢(shì),揚(yáng)程模擬值降低在3.9%以?xún)?nèi);效率模擬曲線(xiàn)與試驗(yàn)曲線(xiàn)變化趨勢(shì)也基本一致,模擬效率值增大在3.4%以?xún)?nèi),模擬值與試驗(yàn)值吻合較好,從而驗(yàn)證了數(shù)值模擬的可靠性;布置三種結(jié)構(gòu)離心泵的揚(yáng)程和效率模擬曲線(xiàn)與原型泵的模擬曲線(xiàn)也有相同變化趨勢(shì),在設(shè)計(jì)工況點(diǎn),布置凹槽后揚(yáng)程下降1.7%,效率增大2.4%,布置橫向障礙物和不連續(xù)障礙物后揚(yáng)程分別增大2.2%、1.6%,效率分別下降2.6%、2.3%,說(shuō)明布置三種結(jié)構(gòu)對(duì)離心泵的外特性性能影響較小。

(a) 流量-揚(yáng)程曲線(xiàn)

3.4 空化試驗(yàn)驗(yàn)證

在流體機(jī)械計(jì)算中通常采用無(wú)量綱空化數(shù)來(lái)描述空化發(fā)生的嚴(yán)重程度,公式定義如下

(12)

(13)

式中:P1為參考?jí)毫Γ捎帽玫倪M(jìn)口壓力;U為基準(zhǔn)速度,采用葉片進(jìn)口邊和前蓋板相交點(diǎn)的圓周速度;n為葉輪的轉(zhuǎn)速;D0為葉片進(jìn)口邊和前蓋板交界處直徑。

圖7為原型泵和布置三種結(jié)構(gòu)離心泵在設(shè)計(jì)工況1.0Q下的空化性能曲線(xiàn)。從圖中可以看出,原型泵的空化性能模擬曲線(xiàn)與試驗(yàn)曲線(xiàn)變化趨勢(shì)基本一致,與試驗(yàn)值相比,模擬值增大在4.1%以?xún)?nèi),驗(yàn)證了空化模型的可行性;布置三種結(jié)構(gòu)離心泵的空化性能模擬曲線(xiàn)與原型泵的模擬曲線(xiàn)變化趨勢(shì)也相同,布置凹槽后斷裂揚(yáng)程提高9.2%,布置橫向障礙物后斷裂揚(yáng)程提高18.1%,布置不連續(xù)障礙物后斷裂揚(yáng)程提高13.1%,說(shuō)明布置三種結(jié)構(gòu)均能提高離心泵的空化性能。

圖7 空化性能曲線(xiàn)

4 計(jì)算結(jié)果及分析

4.1 三種結(jié)構(gòu)對(duì)湍動(dòng)能分布的影響

湍動(dòng)能描述的是湍流的脈動(dòng)程度,隨著湍動(dòng)能的增大,流動(dòng)損失也會(huì)增大。圖8為不同空化數(shù)下葉輪中間截面的湍動(dòng)能分布。從圖中可以看出:隨著空化數(shù)的減小,在原模型和布置三種結(jié)構(gòu)的葉輪中,湍動(dòng)能分布區(qū)域從葉輪進(jìn)口向葉輪出口方向延伸,在葉輪出口處由于葉輪與蝸殼之間的動(dòng)靜干涉作用引起了葉輪內(nèi)流動(dòng)的不穩(wěn)定使得在葉輪出口處的湍動(dòng)能分布區(qū)域最大;而在同一空化數(shù)下,與原模型葉輪相比,布置三種結(jié)構(gòu)的葉輪中湍動(dòng)能強(qiáng)度明顯降低,分布區(qū)域顯著縮小,說(shuō)明能量耗散減少,流動(dòng)變得相對(duì)穩(wěn)定,對(duì)流場(chǎng)結(jié)構(gòu)有明顯改善,其中布置不連續(xù)障礙物的葉輪改善效果最佳。

(a) σ=1.16

4.2 三種結(jié)構(gòu)對(duì)壓力分布的影響

圖9為原模型和布置三種結(jié)構(gòu)的葉輪流道中間切面(0.5 span)的絕對(duì)壓力分布圖,其中span定義為葉輪前蓋板與后蓋板之間的無(wú)量綱距離。從圖9中可以看出:隨著空化數(shù)的減小,在原模型和布置三種結(jié)構(gòu)的葉輪進(jìn)口處空化發(fā)生區(qū)域增大,葉輪內(nèi)空泡體積也增大,壓力分布也有明顯變化,表現(xiàn)為葉輪進(jìn)口處低壓區(qū)增大并向出口處延伸,這也是離心泵內(nèi)空化發(fā)展的主要原因;而在同一空化數(shù)下,與原模型葉輪相比,布置三種結(jié)構(gòu)的葉輪中低壓分布范圍明顯減小,這是由于布置三種結(jié)構(gòu)后均在不同程度上誘發(fā)了葉片進(jìn)口處相對(duì)高壓區(qū)域的形成,而壓力分布的變化會(huì)導(dǎo)致空泡形態(tài)的改變,因此相對(duì)高壓區(qū)的形成會(huì)對(duì)空泡的產(chǎn)生有明顯抑制作用。其中布置橫向障礙物的葉輪中低壓分布范圍減小的最為明顯,對(duì)空泡的抑制效果最好。

(a) σ=1.16

4.3 三種結(jié)構(gòu)對(duì)空泡體積的影響

葉輪內(nèi)空泡體積為

(14)

式中:N為葉輪中控制單元總數(shù)量;αv,i為每個(gè)控制單元中的蒸汽體積分?jǐn)?shù);Vi為每個(gè)控制單元體積。

圖10為不同空化數(shù)下周期內(nèi)葉輪中空泡總體積的平均值。從圖中可以看出:在原模型和布置三種結(jié)構(gòu)的葉輪中,隨著空化數(shù)的減小,葉輪內(nèi)空泡體積顯著增大,葉片進(jìn)口邊空化區(qū)域有一定的擴(kuò)張;而在同一空化數(shù)下,與原模型葉輪相比,布置三種結(jié)構(gòu)的葉輪中空泡體積明顯減少,當(dāng)空化數(shù)σ=1.16時(shí),布置橫向障礙物很大程度地減少了葉輪內(nèi)空泡體積,減少量為82%,對(duì)空化抑制效果最為明顯,布置凹槽和不連續(xù)障礙物減少量分別為41%、67%,對(duì)空化抑制效果也很明顯;空化數(shù)減小至σ=0.82時(shí),布置橫向障礙物對(duì)空化抑制效果最佳,葉輪內(nèi)空泡體積減少量為75%,布置凹槽和不連續(xù)障礙物對(duì)空化抑制效果也很明顯,減少量分別為46%、52%。

(a) σ=1.16

4.4 三種結(jié)構(gòu)對(duì)速度分布的影響

圖11為原模型和布置三種結(jié)構(gòu)的葉輪中間截面的速度矢量圖。從圖中可以看出:隨著空化數(shù)的減小,在原模型和布置三種結(jié)構(gòu)的葉輪中,葉輪流速都有所增大,而且流速分布較為一致,說(shuō)明在空化初生階段,空泡對(duì)葉輪流道的阻塞程度較小,葉輪內(nèi)流速較為平穩(wěn);而在同一空化數(shù)下,與原模型葉輪相比,布置三種結(jié)構(gòu)的葉輪中流速都有不同程度的減小,布置凹槽后,部分流體進(jìn)入凹槽旋轉(zhuǎn)后匯入主流減小了主流流速,從而誘發(fā)了凹槽附近局部壓力增大,抑制了空泡的產(chǎn)生;布置橫向障礙物與不連續(xù)障礙物后,由于障礙物阻擋了部分流體使得主流流速降低,從而誘發(fā)了障礙物附近局部壓力增大,對(duì)空泡的產(chǎn)生起到了抑制作用。

(a) σ=1.16

5 結(jié) 論

(1) 離心泵外特性試驗(yàn)與數(shù)值模擬結(jié)果吻合較好,而且布置三種結(jié)構(gòu)后,在設(shè)計(jì)工況點(diǎn),布置凹槽后揚(yáng)程下降1.7%,效率增大2.4%,布置橫向障礙物和不連續(xù)障礙物后揚(yáng)程分別增大2.2%、1.6%,效率分別下降2.6%、2.3%,說(shuō)明布置三種結(jié)構(gòu)對(duì)離心泵的外特性性能影響較小。

(2) 在空化初生階段,布置三種結(jié)構(gòu)的葉輪中湍動(dòng)能強(qiáng)度均有不同程度的下降,對(duì)流場(chǎng)結(jié)構(gòu)有明顯改善。其中,布置不連續(xù)障礙物葉輪改善效果最佳是由于中間留有間隙能夠更好的改善近壁面流場(chǎng)分布。

(3) 在空化初生階段,布置三種結(jié)構(gòu)的葉輪中進(jìn)口低壓分布區(qū)域均有不同程度的收縮,而相對(duì)高壓區(qū)的形成對(duì)空泡的產(chǎn)生起到了很大的抑制作用,使空泡體積均有很大程度的減少。其中,布置橫向障礙物葉輪中低壓分布范圍減小的最為明顯,空泡體積減少的最多,抑制空化效果最好。

(4) 在空化初生階段,布置三種結(jié)構(gòu)后,葉輪流道內(nèi)流速都有不同程度的減小,使葉輪內(nèi)局部壓力增大,對(duì)空泡的產(chǎn)生起到了抑制作用。其中,布置凹槽后由于部分流體進(jìn)入凹槽旋轉(zhuǎn)后匯入主流從而減小了主流流速,布置橫向障礙物與不連續(xù)障礙物后由于障礙物阻擋了部分流體使得主流流速降低。

猜你喜歡
空泡揚(yáng)程離心泵
一種改善離心泵運(yùn)行狀態(tài)的方法
大型立式單級(jí)引黃離心泵軸向力平衡的研究
離心泵流場(chǎng)外特性及激勵(lì)數(shù)值計(jì)算
管路受力誘發(fā)高揚(yáng)程離心泵振動(dòng)加劇原因分析
高比速離心泵的高效設(shè)計(jì)研究
水下航行體雙空泡相互作用數(shù)值模擬研究
基于LPV的超空泡航行體H∞抗飽和控制
基于CFD的對(duì)轉(zhuǎn)槳無(wú)空泡噪聲的仿真預(yù)報(bào)
船海工程(2015年4期)2016-01-05 15:53:28
新型中高揚(yáng)程大流量水錘泵結(jié)構(gòu)技術(shù)改進(jìn)研究
SPH在水下高速物體空泡發(fā)展模擬中的應(yīng)用
二手房| 利川市| 体育| 庆元县| 江阴市| 疏勒县| 禄丰县| 农安县| 长宁县| 邢台县| 漳平市| 洛南县| 手游| 恩施市| 册亨县| 锡林郭勒盟| 惠东县| 开封市| 岑溪市| 巫溪县| 儋州市| 资兴市| 安吉县| 大庆市| 班玛县| 改则县| 麻栗坡县| 云龙县| 吉水县| 唐山市| 元阳县| 平阴县| 运城市| 吉隆县| 萍乡市| 酒泉市| 盱眙县| 治县。| 新密市| 桐梓县| 乐东|