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

?

鑄態(tài)50Cr5MoV鋼的動(dòng)態(tài)再結(jié)晶行為分析

2017-02-05 08:59:24王耀琨李釗庫(kù)
大型鑄鍛件 2017年1期
關(guān)鍵詞:鑄態(tài)再結(jié)晶軋輥

王耀琨 杜 佳 陳 陽(yáng) 李釗庫(kù)

(中冶陜壓重工設(shè)備有限公司,陜西711711)

鑄態(tài)50Cr5MoV鋼的動(dòng)態(tài)再結(jié)晶行為分析

王耀琨 杜 佳 陳 陽(yáng) 李釗庫(kù)

(中冶陜壓重工設(shè)備有限公司,陜西711711)

采用Gleeble-3800熱模擬實(shí)驗(yàn)機(jī)對(duì)鑄態(tài)50Cr5MoV鋼進(jìn)行熱壓縮實(shí)驗(yàn),研究其在變形溫度為900℃、950℃、1 000℃、1 050℃、1 100℃和應(yīng)變速率為0.01 s-1、0.1 s-1、1 s-1條件下的動(dòng)態(tài)再結(jié)晶行為。結(jié)果表明:鑄態(tài)50Cr5MoV鋼在熱變形過(guò)程中發(fā)生了動(dòng)態(tài)再結(jié)晶,且隨著應(yīng)變速率的降低以及變形溫度的升高,動(dòng)態(tài)再結(jié)晶越容易發(fā)生。

鑄態(tài)50Cr5MoV鋼;動(dòng)態(tài)再結(jié)晶;本構(gòu)方程;臨界應(yīng)變模型

軋輥是軋機(jī)的重要設(shè)備,由于軋機(jī)特殊的工作環(huán)境,要求軋輥必須具有良好的耐磨性、高屈服強(qiáng)度和抗剝落性能,輥頸具備高強(qiáng)韌性和抗斷裂性能。50Cr5MoV鋼是繼第一代2%Cr材質(zhì)、第二代3%Cr材質(zhì)之后研發(fā)出的第三代鍛輥用鋼,其淬透性、疲勞強(qiáng)度以及抗事故性能均比第一代、第二代材質(zhì)有明顯提高,已逐漸成為現(xiàn)代化板帶軋機(jī)軋輥的主要用鋼,并將在較長(zhǎng)一段時(shí)間內(nèi)占據(jù)主導(dǎo)地位[1-2]。

動(dòng)態(tài)再結(jié)晶是金屬熱加工生產(chǎn)過(guò)程中控制組織與性能的有效手段之一[3-4]。動(dòng)態(tài)再結(jié)晶使得50Cr5MoV鋼在熱鍛成型過(guò)程中的微觀(guān)組織發(fā)生轉(zhuǎn)變,通過(guò)調(diào)節(jié)熱鍛工藝,可以細(xì)化奧氏體組織,優(yōu)化材料的力學(xué)性能,從而提高50Cr5MoV軋輥的使用壽命。本文利用Gleeble-3800熱模擬機(jī)對(duì)鑄態(tài)50Cr5MoV軋輥鋼進(jìn)行高溫壓縮試驗(yàn),研究其動(dòng)態(tài)再結(jié)晶行為,并建立了動(dòng)態(tài)再結(jié)晶模型,為制定合理的鍛造工藝提供理論依據(jù),從而優(yōu)化鍛鋼軋輥的綜合性能。

1 實(shí)驗(yàn)材料及方法

本實(shí)驗(yàn)所用設(shè)備為Gleeble-3800熱力模擬實(shí)驗(yàn)機(jī),實(shí)驗(yàn)材料為50Cr5MoV軋輥鋼(鑄錠),其化學(xué)成分如表1。

表1 50Cr5MoV軋輥鋼的化學(xué)成分(質(zhì)量分?jǐn)?shù),%)Table 1 Chemical compositions of 50Cr5MoV roll steel (mass fraction, %)

從鑄錠上截取樣品,制成?10 mm×15 mm的圓柱型試樣,在Gleeble-3800型熱模擬機(jī)上進(jìn)行壓縮實(shí)驗(yàn)。以10℃/s速度將試樣加熱至1 150℃以便奧氏體化,進(jìn)行3 min保溫,然后將試樣以10℃/s的冷速分別冷卻至900℃、950℃、1 000℃、1 050℃和1 100℃,為了消除試樣內(nèi)部的溫度差,樣品保溫20 s,進(jìn)行壓縮變形,變形量按70%控制,變形速率分別按照0.01 s-1、0.1 s-1和1 s-1進(jìn)行控制。實(shí)驗(yàn)過(guò)程中記錄應(yīng)力-應(yīng)變數(shù)據(jù),變形后的樣品立即淬火,將奧氏體組織固定到室溫。實(shí)驗(yàn)過(guò)程如圖1所示。

2 實(shí)驗(yàn)結(jié)果與分析

通過(guò)壓縮變形過(guò)程所記錄的數(shù)據(jù),采用Origin繪圖軟件對(duì)數(shù)據(jù)進(jìn)行處理,可以得到50Cr5MoV軋輥鋼在不同實(shí)驗(yàn)條件下的真應(yīng)力-真應(yīng)變曲線(xiàn),如圖2所示。

圖1 實(shí)驗(yàn)鋼熱模擬實(shí)驗(yàn)工藝Figure 1 The thermal simulation process of tested steel

(a)

(b)

(c)圖2 50Cr5MoV軋輥鋼的真應(yīng)力-應(yīng)變曲線(xiàn)Figure 2 The true stress-strain curve of 50Cr5MoV roll steel

由圖2可知,50Cr5MoV鋼流變應(yīng)力受溫度及應(yīng)變速率的影響顯著。在變形開(kāi)始階段,隨著形變量的不斷增加,流變應(yīng)力以較快的速度增大。分析認(rèn)為,在變形開(kāi)始階段,形變量增加,位錯(cuò)發(fā)生了滑移等運(yùn)動(dòng),位錯(cuò)運(yùn)動(dòng)引起的軟化作用不足以抵消位錯(cuò)密度增加所導(dǎo)致的硬化作用,因此,流變應(yīng)力迅速增加,形變強(qiáng)化作用顯著[5]。后續(xù)變形的曲線(xiàn)呈現(xiàn)兩種不同的主導(dǎo)軟化態(tài)勢(shì),即動(dòng)態(tài)再結(jié)晶(DRX)和動(dòng)態(tài)回復(fù)(DRV)。

在0.01 s-1~0.1 s-1的低應(yīng)變速率條件下,如圖2(b)所示,流變曲線(xiàn)具有明顯的應(yīng)力峰值,說(shuō)明壓縮變形過(guò)程中發(fā)生了動(dòng)態(tài)再結(jié)晶,隨著動(dòng)態(tài)再結(jié)晶程度的不斷增大直至發(fā)生完全動(dòng)態(tài)再結(jié)晶,位錯(cuò)增值所引起的加工硬化與位錯(cuò)運(yùn)動(dòng)、位錯(cuò)脫釘?shù)纫鸬膭?dòng)態(tài)軟化達(dá)到動(dòng)態(tài)平衡,變形進(jìn)入了穩(wěn)態(tài)流動(dòng)階段[5]。

在1 s-1的較高的應(yīng)變速率條件下,如圖2(c)中溫度為900℃的曲線(xiàn),應(yīng)力到達(dá)一定程度后,隨著應(yīng)變量的增加,流變應(yīng)力趨于平穩(wěn),為動(dòng)態(tài)回復(fù)型軟化。這是因?yàn)榧庸び不什皇请S著變形程度的增加而一直增大,在材料發(fā)生加工硬化的同時(shí),變形產(chǎn)生的位錯(cuò)能夠通過(guò)交滑移和攀移等方式運(yùn)動(dòng),使部分位錯(cuò)消失,部分位錯(cuò)重排,從而造成金屬材料的回復(fù)軟化。這種情況下,整個(gè)軟化過(guò)程由動(dòng)態(tài)回復(fù)控制[6]。

當(dāng)應(yīng)變速率一定時(shí),動(dòng)態(tài)再結(jié)晶行為隨著變形溫度的升高越明顯,而最大流變應(yīng)力相應(yīng)降低。因?yàn)樽冃螠囟壬咴鰪?qiáng)了變形的熱激活作用,晶界可動(dòng)性也增強(qiáng),有利于動(dòng)態(tài)再結(jié)晶形核和晶粒長(zhǎng)大以及位錯(cuò)湮滅,這將會(huì)降低流變應(yīng)力[7]。當(dāng)變形溫度一定時(shí),再結(jié)晶晶粒形核和長(zhǎng)大時(shí)間隨著應(yīng)變速率的降低變得更為充分,動(dòng)態(tài)再結(jié)晶進(jìn)行得越充分。同樣位錯(cuò)之間的交互作用導(dǎo)致位錯(cuò)密度急劇增加,增加了變形儲(chǔ)存能,從而引起再結(jié)晶的臨界溫度的降低,動(dòng)態(tài)再結(jié)晶程度也相應(yīng)增加。

3 動(dòng)態(tài)再結(jié)晶動(dòng)力學(xué)模型

3.1 熱變形本構(gòu)關(guān)系模型

研究表明,金屬材料的高溫流變應(yīng)力既與成分有關(guān),又與形變溫度和應(yīng)變速率有關(guān),在材料的成分不變時(shí),采用Arrhenius雙曲正弦函數(shù)描述50Cr5MoV鋼的本構(gòu)關(guān)系[8]:

C.M.Sellars和W.J.M.Tegart提出并實(shí)驗(yàn)驗(yàn)證了熱變形條件,通常可以用溫度補(bǔ)償?shù)膽?yīng)變速率因子Zener-Hollomon參數(shù)Z來(lái)表示,Z參數(shù)可用下式表示[9]:

根據(jù)應(yīng)力條件高低的不同,式(1)可以簡(jiǎn)化為以下兩種方式:

其中,系數(shù)α、β及n1之間滿(mǎn)足α=β/n1。

因此,可以通過(guò)式(3)和式(4)確定α值。根據(jù)文獻(xiàn)的求解方法[10],可得α=0.007 7。

當(dāng)對(duì)式(1)進(jìn)行Tayler級(jí)數(shù)展開(kāi)時(shí)發(fā)現(xiàn),在低應(yīng)力條件下較為接近式(3),在高應(yīng)力條件下接近式(4),因而采用式(1)基本能適用于描述整個(gè)應(yīng)力范圍的熱變形過(guò)程,且眾多研究結(jié)果[11-12]表明,式(1)能夠較好地描述金屬材料的常規(guī)熱變形過(guò)程。

對(duì)式(1)進(jìn)行變換處理可得:

對(duì)式(5)兩邊取對(duì)數(shù)并求偏微分,可得熱變形激活能Qdef的表達(dá)式:

經(jīng)過(guò)計(jì)算,可得50Cr5MoV鋼熱變形激活能Qdef=333.47kJ/mol,b=5 689.33,n=7.05。因此,50Cr5MoV鋼的動(dòng)態(tài)再結(jié)晶數(shù)學(xué)模型為:

(7)

熱變形過(guò)程中,通過(guò)對(duì)lnZ和ln[sinh(ασp)]進(jìn)行線(xiàn)性擬合(見(jiàn)圖4),可求得lnA=28.11。結(jié)合圖4、式(7),50Cr5MoV鋼的高溫?zé)嶙冃伪緲?gòu)方程可表示為:

(8)

圖4 ln[sinh(ασp)]與lnZ的關(guān)系圖Figure 4 The relationship between ln[sinh(ασp)] and lnZ

3.2 動(dòng)態(tài)再結(jié)晶臨界應(yīng)變模型

動(dòng)態(tài)再結(jié)晶是否發(fā)生,通常用臨界應(yīng)變量(εc)來(lái)判斷,只有形變量大于εc時(shí),才會(huì)發(fā)生動(dòng)態(tài)再結(jié)晶[13]。因此,準(zhǔn)確地確定臨界應(yīng)變量εc,對(duì)研究熱變形工藝參數(shù)尤為重要。在以往的科研工作中,研究人員提出了許多數(shù)學(xué)模型用于預(yù)測(cè)動(dòng)態(tài)再結(jié)晶開(kāi)始的臨界條件[5,14]。Poliak和Jone等人提出了計(jì)算動(dòng)態(tài)再結(jié)晶臨界應(yīng)變的方法,可由下式表示:

式中,θ表示加工硬化率,θ=dσ/dε。

動(dòng)態(tài)再結(jié)晶發(fā)生時(shí)的臨界應(yīng)變點(diǎn)由(-dθ/dε)-ε關(guān)系曲線(xiàn)確定。而Najafizadeh-Jonas對(duì)以上方法做了簡(jiǎn)化,通過(guò)lnθ-ε曲線(xiàn)分析法確定臨界應(yīng)變。在本研究工作中,采用類(lèi)似Najafizadeh-Jonas的方法,通過(guò)對(duì)lnθ-ε曲線(xiàn)的三次多項(xiàng)式擬合來(lái)確定拐點(diǎn)的方法來(lái)確定動(dòng)態(tài)再結(jié)晶臨界應(yīng)變(εc)[15]。該多項(xiàng)式如下:

式中,A1、A2、A3、A4為相應(yīng)變形條件下的常數(shù)。

在動(dòng)態(tài)再結(jié)晶開(kāi)始時(shí)的臨界應(yīng)變處,有下式成立:

可得:

由于實(shí)驗(yàn)鋼的實(shí)際流變應(yīng)力曲線(xiàn)并不是數(shù)學(xué)意義上的光滑曲線(xiàn),很難直接測(cè)量或計(jì)算出加工硬化率θ,因此,先將應(yīng)力-應(yīng)變曲線(xiàn)進(jìn)行九次多項(xiàng)式函數(shù)進(jìn)行擬合,再進(jìn)行微分運(yùn)算。變形溫度1 000℃,變形速率分別為0.01s-1、0.1s-1、1s-1時(shí)的lnθ-ε曲線(xiàn)及其三次多項(xiàng)式擬合值如圖5,圖中箭頭所指拐點(diǎn)為開(kāi)始發(fā)生動(dòng)態(tài)再結(jié)晶的臨界點(diǎn)(εc)。從圖中能夠看出,實(shí)驗(yàn)值和計(jì)算值吻合度較高,說(shuō)明三次多項(xiàng)式擬合可以描述加工硬化率與應(yīng)變之間的關(guān)系。

通過(guò)式(12)可以計(jì)算變形條件為1 000℃、應(yīng)變速率0.01s-1、0.1s-1、1s-1時(shí)的εc分別為0.095、0.119和0.164。根據(jù)峰值應(yīng)變,可求得對(duì)應(yīng)的εc/εp值分別為0.511、0.511和0.524。通過(guò)以上方法,統(tǒng)計(jì)不同應(yīng)變條件時(shí)的εc和εp(表2)??汕蟮忙與/εp的平均值如下:

根據(jù)臨界應(yīng)變、峰值應(yīng)變數(shù)據(jù)及Z參數(shù)的計(jì)算結(jié)果,繪制50Cr5MoV鋼動(dòng)態(tài)再結(jié)晶lnεp-lnZ和lnεc-lnZ的關(guān)系曲線(xiàn)圖(圖6)??梢杂^(guān)察到,Z參數(shù)可以反映臨界應(yīng)變與變形條件之間的規(guī)律,且lnε與lnZ呈現(xiàn)較好的線(xiàn)性吻合。通過(guò)線(xiàn)性擬合得到動(dòng)態(tài)再結(jié)晶臨界模型如下:

圖5 加工硬化率與應(yīng)變關(guān)系曲線(xiàn)及其三次多項(xiàng)式擬合結(jié)果Figure 5 The relation curve between work hardening ratesand strains and its three times polynomial fitting results

溫度℃0.01s-10.1s-11s-1εcεpεc/εpεcεpεc/εpεcεpεc/εp9009501000105011000.1070.0980.0950.0820.0760.2050.1910.1860.1550.1510.5220.5130.5110.5290.5030.1350.1260.1190.1080.0950.2650.2450.2330.2120.1910.5090.5140.5110.5090.4970.2110.1720.1640.1420.1280.4250.3430.3130.2770.2550.4960.5010.5240.5130.502

圖6 lnεp-lnZ和lnεc-lnZ的關(guān)系Figure 6 The relationship between lnεp-lnZ and lnεc-lnZ

3.3 動(dòng)態(tài)再結(jié)晶晶粒尺寸模型

動(dòng)態(tài)再結(jié)晶過(guò)程中,變形工藝對(duì)再結(jié)晶晶粒尺寸影響很大。當(dāng)溫度升高時(shí),位錯(cuò)的運(yùn)動(dòng)比低溫時(shí)更容易發(fā)生,這有利于動(dòng)態(tài)再結(jié)晶的形核與晶粒長(zhǎng)大[16]。動(dòng)態(tài)再結(jié)晶晶粒尺寸Ddrex正比于Z參數(shù)的指函數(shù)[17],其表達(dá)式為:

其中,C、a為與材料有關(guān)的常數(shù)。

利用Image tool軟件,統(tǒng)計(jì)了完全動(dòng)態(tài)再結(jié)晶時(shí)的晶粒尺寸,計(jì)算其平均值(表3),將不同變形條件下測(cè)得的晶粒尺寸d代入式(16),進(jìn)行回歸處理(圖7),得到不同變形參數(shù)下再結(jié)晶晶粒尺寸lnDdrex與lnZ呈線(xiàn)性關(guān)系,從而求得C=2.64×103,a=-0.159。得出再結(jié)晶晶粒尺寸與變形參數(shù)之間的模型為:

Ddrex=2.64×103Z-0.159 (17)

圖7 lnDdrex與lnZ的關(guān)系Figure 7 The relationshop between lnDdrex and lnZ

4 結(jié)論

(1)通過(guò)熱模擬壓縮變形實(shí)驗(yàn),可知50Cr5MoV鋼存在動(dòng)態(tài)再結(jié)晶,并且隨著應(yīng)變速率的降低,形變溫度的升高,動(dòng)態(tài)再結(jié)晶臨界應(yīng)變量越小,峰值應(yīng)力越大,再結(jié)晶造成的軟化越明顯。

(2)對(duì)熱壓縮實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析,得到50Cr5MoV軋輥鋼的動(dòng)態(tài)再結(jié)晶激活能為333.47 kJ/mol。

(3)通過(guò)對(duì)鑄態(tài)50Cr5MoV鋼的動(dòng)態(tài)再結(jié)晶行為研究,建立了動(dòng)態(tài)再結(jié)晶本構(gòu)方程、動(dòng)態(tài)再結(jié)晶臨界應(yīng)變模型以及動(dòng)態(tài)再結(jié)晶晶粒尺寸模型,為大型鍛鋼軋輥的成型工藝奠定理論基礎(chǔ)。

[1] 俞誓達(dá),陳菡.我圍軋輥業(yè)現(xiàn)狀及發(fā)展中應(yīng)注意的問(wèn)題[J]. 鋼鐵.2007,42(7):1.

[2] 武建平,賈建平,王自東.V-N微合金化含鎢改進(jìn)型50Cr5MoVW冷軋工作輥用鋼[J].特殊鋼,2007,28(5):40.

[3] Alexander T, Piotr W, Mihaela A. Modelling the dynamic recrystallization in C-Mn micro-alloyed steel during thermo-mechanical treatment using cellular automata[J] .Computational Materials Science,2014,94:85-94.

[4] 向嵩,譚智林. Q690低碳微合金鋼熱變形微觀(guān)組織演變及加工圖[J]. 湖南大學(xué)學(xué)報(bào)(自然科學(xué)版),2013, 40(4):80-86.

[5] Cao Yu,Di Hong-shuang,Zhang Jin-qi. Research on hot deformation behavior and hot work ability of alloy 800H[J].Acta Metallurgica Sinica, 2013, 7 (49):811-821.

[6] 易嘯. 2Cr13鋼熱變形行為研究[D]. 沈陽(yáng): 東北大學(xué)材料與冶金學(xué)院, 2012 :21-23.

[7] Mcqueen H J,Ryan N D.Constitutive analysis in hot working[J]. Materials Science and Engineering A, 2002, 322:43-46.

[8] Naderi M, Durrenberger L, Bleck W. Constitutive relationships for 22MnB5 boron steel deformed isothermally at high temperatures[J].Materials Science and Engineering, 2008; A478:130.

[9] Zener C, Hollomon J H. Effect of strain-rate upon the plastic flow of steel[J]. J Appl Phys, 1944, 15(1): 22-27.

[10] Mahmoud R, Ghandehari F, Davood N. Modeling the high temperature flow behavior and dynamic recrystallization kinetics of a medium carbon microalloyed steel[J]. ASM International, 2014, 23:1077-1087.

[11] 馬博, 彭艷,劉云飛. 低合金鋼Q345B動(dòng)態(tài)再結(jié)晶動(dòng)力學(xué)模型[J].材料熱處理學(xué)報(bào), 2010, 31(4): 141-147.

[12] Imbert C A C, Mcqueen H J. Peak strength strain hardening and dynamic restoration of A2 and M2 tool steels in hot deformation[J]. Materials Science and Engineering A, 2001, 313:88-103.

[13] Quan Guo-zheng, Zhao Lei, Shi Yu. An identification and characterization for the dynamic recrystallization critical conditions of Ti-6Al-2Zr-1Mo-1V alloy[J]. Journal of Functional Materials, 2012, 43(2):222-226.

[14] Fen Di,Zhang Xin-ming,Liu Sheng-dan. Rate controlling mechanisms in hot deformation of 7A55 aluminum alloy[J].Transactions of Nonferrous Metals Society of China, 2014, 24:28-35.

[15] Jonas J J,Qujelennec X,Jiang L. Theavrami kinetics of dynamic recrystallization[J]. Acta Materialia,2009,57( 9) : 2748-2756.

[16] 余永寧.金屬學(xué)原理[M].北京:冶金工業(yè)出版社,2003.

[17] Wahabi M E,Cabrera J M,Prado J M.Hot working of two AISI 304 steels: a comparative study[J]. Materials Science and Engineering A,2003,343: 116-125.

編輯 杜青泉

Analysis of Dynamic Recrystallization Behaviour of As-cast 50Cr5MoV Steel

Wang Yaokun, Du Jia, Chen Yang, Li Zhaoku

By using Gleeble-3800 thermomechanical simulator, the hot compression test for as-cast 50Cr5MoV steel has been carried out. Under conditions of the deformation temperature of 900℃, 950℃, 1 000℃, 1 050℃ and 1 100℃ and the strain rate of 0.01 s-1, 0.1 s-1, 1 s-1, the dynamic recrystallization behaviour of steel has been studied. The results show that the dynamic recrystallization appears during the hot deformation process of as-cast 50Cr5MoV steel. And the dynamic recrystallization is more likely to appear along with the reducing of the strain rate and the increasing of the deformation temperature.

as-cast 50Cr5MoV steel; dynamic recrystallization; constitutive equation; critical strain model

2016—09—01

王耀琨(1964—)男,本科,高級(jí)工程師。

TG111.7

B

猜你喜歡
鑄態(tài)再結(jié)晶軋輥
梅鋼1780產(chǎn)線(xiàn)高速鋼軋輥使用維護(hù)
鑄態(tài)QTRSi4Mo1材料的研制
?;に噷?duì)低溫Hi-B鋼初次及二次再結(jié)晶的影響
上海金屬(2016年3期)2016-11-23 05:19:38
數(shù)控軋輥磨床主輔系統(tǒng)分離的探討
鑄態(tài)30Cr2Ni4MoV鋼動(dòng)態(tài)再結(jié)晶行為研究
大型鑄鍛件(2015年1期)2016-01-12 06:32:58
Cu元素對(duì)7XXX 系列鋁合金再結(jié)晶的影響
上海金屬(2014年3期)2014-12-19 13:09:04
Q460GJE鋼形變奧氏體的動(dòng)態(tài)再結(jié)晶行為研究
上海金屬(2014年3期)2014-12-19 13:09:03
軋輥用012Al鋼熱處理試驗(yàn)研究
上海金屬(2014年1期)2014-12-18 06:52:03
一種鑄態(tài)含氮M2高速鋼的熱變形行為研究
成分、工藝和組織對(duì)鑄態(tài)球墨鑄鐵低溫韌性的影響①
那坡县| 本溪市| 罗江县| 乐山市| 隆回县| 临海市| 南宫市| 南阳市| 广宗县| 长丰县| 宣汉县| 乌拉特中旗| 五家渠市| 从化市| 桂平市| 江安县| 左贡县| 麻江县| 陆丰市| 通山县| 梁山县| 南靖县| 莱阳市| 金乡县| 万载县| 苏尼特左旗| 集安市| 策勒县| 太康县| 柳林县| 高州市| 峨边| 汤阴县| 青河县| 江山市| 中方县| 大关县| 巴里| 新营市| 广宗县| 双城市|