• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    風(fēng)扇葉片凸肩工作面裂紋起始原因分析

    2015-11-19 08:43:02張海洋王繼業(yè)
    航空發(fā)動(dòng)機(jī) 2015年6期
    關(guān)鍵詞:磨痕倒角錯(cuò)位

    張海洋,曹 航,李 洋,王繼業(yè),邵 帥

    (中航工業(yè)沈陽(yáng)發(fā)動(dòng)機(jī)設(shè)計(jì)研究所,沈陽(yáng)110015)

    0 引言

    在航空發(fā)動(dòng)機(jī)風(fēng)扇葉片凸肩結(jié)構(gòu)設(shè)計(jì)中,主要從減小氣動(dòng)影響[1-2]和保證干摩擦減振效果[3-13]兩方面進(jìn)行設(shè)計(jì),而對(duì)其本身的強(qiáng)度設(shè)計(jì)研究較少。但在實(shí)際使用過(guò)程中,經(jīng)常出現(xiàn)凸肩工作面表面耐磨層掉渣等現(xiàn)象,甚至在工作面上發(fā)現(xiàn)了裂紋。對(duì)某型風(fēng)扇葉片工作后分解檢查發(fā)現(xiàn),風(fēng)扇葉片盆側(cè)凸肩工作面存在長(zhǎng)約3mm的裂紋,位于凸肩工作面尖部位置,并向上、下端面延伸一定距離。

    本文對(duì)凸肩工作面上出現(xiàn)的裂紋進(jìn)行了裂紋斷口分析。通過(guò)工作狀態(tài)下凸肩工作面接觸分析,確定了凸肩工作面擠壓應(yīng)力分布和變形錯(cuò)位規(guī)律,并分析了工作面初始裝配錯(cuò)位和工作面邊緣倒角等裝配與加工因素對(duì)擠壓應(yīng)力分布的影響,最后找到了初始裂紋產(chǎn)生原因。

    1 裂紋形貌和斷口分析檢查

    凸肩工作面裂紋形貌如圖1所示。裂紋位于葉片葉盆側(cè)凸肩耐磨層部位,靠近凸肩尖部,向上、下端面延伸一定距離。根據(jù)斷口分析,該裂紋為疲勞裂紋,起源于靠近凸肩下表面轉(zhuǎn)接區(qū)的凸肩工作面,如圖2中箭頭所指位置。

    圖1 凸肩工作面裂紋形貌

    圖2 裂紋斷口宏觀形貌

    在裂紋起始位置附近對(duì)凸肩工作面進(jìn)行電鏡檢查,耐磨涂層檢查結(jié)果如圖3所示。從圖中可見(jiàn),在裂紋起始位置,凸肩工作面表面耐磨層呈現(xiàn)多條微裂紋,表明在工作過(guò)程中受到較大的擠壓磨蹭應(yīng)力。而在遠(yuǎn)離裂紋起始位置的耐磨涂層檢查結(jié)果顯示表面無(wú)微裂紋,如圖4所示。

    圖3 裂紋起源附近凸肩工作面電鏡檢查

    圖4 遠(yuǎn)離裂紋起源附近凸肩工作面電鏡檢查

    對(duì)出現(xiàn)裂紋的凸肩工作面進(jìn)行表面磨痕檢查,盆側(cè)凸肩(存在裂紋)及與其配合的背側(cè)凸肩工作面表面磨痕分別如圖5、6所示。從圖中可見(jiàn),凸肩工作面尖部均為磨損較輕區(qū)域,葉盆、葉背側(cè)凸肩工作面分別靠近上、下表面轉(zhuǎn)接區(qū)域大部分都為磨損較輕區(qū)域。

    圖5 出現(xiàn)裂紋盆側(cè)凸肩工作面磨痕

    圖6 配合的背側(cè)凸肩工作面磨痕

    同時(shí),出現(xiàn)裂紋的盆側(cè)凸肩工作面磨痕表明存在初始裝配錯(cuò)位;并且出現(xiàn)裂紋的凸肩工作面邊緣倒角非常小(如圖2中工作面與上表面交接部位倒角),小于設(shè)計(jì)要求。

    綜上所述,根據(jù)裂紋形貌以及斷口分析檢查得出:

    (1)裂紋為疲勞裂紋,起源于靠近凸肩下表面轉(zhuǎn)接區(qū)的凸肩工作面,且此處無(wú)明顯的結(jié)構(gòu)和材質(zhì)缺陷;

    (2)在裂紋起源位置附近發(fā)現(xiàn)耐磨層表面存在較大的擠壓磨蹭應(yīng)力導(dǎo)致的微裂紋;

    (3)凸肩工作面存在工作中變形錯(cuò)位現(xiàn)象;

    (4)存在初始裝配錯(cuò)位;

    (5)邊緣倒角較小。

    根據(jù)檢查結(jié)果,推測(cè)裂紋產(chǎn)生原因?yàn)槿~片工作后,凸肩工作面發(fā)生變形錯(cuò)位,在裂紋起始位置出現(xiàn)較大的擠壓應(yīng)力集中[14-15];同時(shí)工作面的初始裝配錯(cuò)位和邊緣倒角較小等加工裝配因素加劇了擠壓應(yīng)力集中程度,使裂紋起源部位耐磨層出現(xiàn)微裂紋,擴(kuò)展形成裂紋。

    2 凸肩工作面接觸模型的建立

    推測(cè)裂紋產(chǎn)生與凸肩工作面擠壓應(yīng)力集中有關(guān),因此需要建立風(fēng)扇凸肩工作面接觸計(jì)算模型,分析工作狀態(tài)凸肩工作面擠壓應(yīng)力分布情況。

    為了節(jié)省計(jì)算時(shí)間,基于ANSYS平臺(tái)進(jìn)行凸肩工作面接觸計(jì)算,建立簡(jiǎn)化的凸肩工作面接觸模型[16-18]。首先,假設(shè)整圈風(fēng)扇葉片變形一致,建立2個(gè)相鄰風(fēng)扇葉片模型,將配合其凸肩工作面定義為接觸;建立1個(gè)經(jīng)過(guò)發(fā)動(dòng)機(jī)軸線的平面將未定義接觸的一側(cè)風(fēng)扇葉片凸肩切開(kāi),并要求不能切到工作面和根部倒角;然后,將切割平面繞發(fā)動(dòng)機(jī)軸線旋轉(zhuǎn)360°/N(N 為葉片數(shù))移動(dòng)到相鄰葉片的相同位置,切開(kāi)凸肩,刪掉兩側(cè)模型,如圖7所示;最后,由于假設(shè)整圈風(fēng)扇葉片變形一致,同時(shí)兩側(cè)的切開(kāi)面為相鄰葉片的同一位置,變形應(yīng)該一致,可以定義3向位移約束,如圖8所示。利用該方法定義的簡(jiǎn)化凸肩工作面接觸模型,可以考慮葉身可加載氣動(dòng)和離心載荷,模擬葉片實(shí)際工作情況;凸肩工作面可設(shè)置接觸參數(shù),如不同初始緊度、初始接觸狀態(tài)和摩擦系數(shù)等;可考慮不同凸肩工作面模型。

    圖7 凸肩工作面接觸計(jì)算模型

    圖8 被切開(kāi)面間的3向位移約束

    3 凸肩工作面擠壓應(yīng)力分布

    3.1 凸肩工作面擠壓應(yīng)力分析

    針對(duì)航空發(fā)動(dòng)機(jī)工作狀態(tài),考慮風(fēng)扇葉片工作載荷,采用凸肩工作面接觸計(jì)算模型,計(jì)算得到設(shè)計(jì)初始緊度時(shí)凸肩工作面擠壓應(yīng)力分布,如圖9所示,平均和最大擠壓應(yīng)力見(jiàn)表1。

    圖9 大狀態(tài)常用初始緊度時(shí)凸肩工作面擠壓應(yīng)力分布

    表1 凸肩工作面平均擠壓應(yīng)力和最大擠壓應(yīng)力

    從圖中可見(jiàn),擠壓應(yīng)力主要分布在凸肩工作面下側(cè),上邊緣擠壓應(yīng)力幾乎為0,表明工作中此處不接觸,可能出現(xiàn)了變形錯(cuò)位;在長(zhǎng)軸方向,靠近尖部存在一片大的擠壓應(yīng)力分布,而在根部工作面擠壓應(yīng)力較均勻且很小,因此,在長(zhǎng)軸方向也存在變形錯(cuò)位;從整個(gè)工作面擠壓應(yīng)力分布分析,整個(gè)工作面存在長(zhǎng)、短軸上的變形錯(cuò)位,導(dǎo)致整個(gè)工作面接觸不均勻,主要在凸肩尖部和下側(cè)接觸,在尖部下側(cè)區(qū)域形成大的擠壓應(yīng)力區(qū)域,在靠近裂紋起始位置的下側(cè)尖部位置形成擠壓應(yīng)力集中。根據(jù)工作面擠壓應(yīng)力分析,得出凸肩工作面呈現(xiàn)接觸不均勻現(xiàn)象,主要在盆側(cè)凸肩下側(cè)靠近尖部接觸,與凸肩工作面磨痕分析結(jié)果吻合;裂紋起始部位存在擠壓應(yīng)力集中。

    3.2 工作狀態(tài)凸肩工作面變形分析

    根據(jù)凸肩工作面接觸計(jì)算,得到工作狀態(tài)凸肩工作面相互位移變形,如圖10所示;采用凸肩局部坐標(biāo)系,如圖11所示;得到凸肩工作面盆側(cè)凸肩位移見(jiàn)表2。

    圖10 工作狀態(tài)凸肩相互變形錯(cuò)位(藍(lán)色為背側(cè)凸肩工作面,彩色為盆側(cè)裂紋凸肩工作面)

    根據(jù)計(jì)算結(jié)果,在工作狀態(tài)下,由于葉片受到的工作載荷和凸肩的初始緊度作用下,凸肩會(huì)在長(zhǎng)、短軸方向都發(fā)生錯(cuò)位變形,具體表現(xiàn)為:

    圖11 凸肩局部坐標(biāo)系

    表2 凸肩工作面盆側(cè)凸肩位移 mm

    (1)凸肩長(zhǎng)軸方向發(fā)生盆側(cè)凸肩尖部變形錯(cuò)位突出,會(huì)導(dǎo)致盆側(cè)凸肩尖部不發(fā)生接觸,磨痕較輕,與工作面磨痕檢查結(jié)果一致;

    (2)凸肩短軸方向發(fā)生盆側(cè)凸肩上側(cè)變形錯(cuò)位突出,會(huì)導(dǎo)致其工作面上側(cè)不發(fā)生接觸,磨痕較輕,與工作面磨痕檢查結(jié)果一致;

    (3)由于凸肩長(zhǎng)、短軸方向的位移,導(dǎo)致凸肩主要在盆側(cè)凸肩下側(cè)靠近尖部接觸,靠近裂紋起始位置出現(xiàn)擠壓應(yīng)力集中。

    4 凸肩初始裝配錯(cuò)位影響

    斷口分析檢查表明裂紋的凸肩工作面存在初始裝配錯(cuò)位。根據(jù)磨痕得到錯(cuò)位方式和尺寸,建立凸肩工作面接觸模型,如圖12所示。計(jì)算得到初始裝配錯(cuò)位情況下的凸肩工作面擠壓應(yīng)力分布,如圖13所示。錯(cuò)位和未錯(cuò)位時(shí)出現(xiàn)裂紋凸肩工作面擠壓應(yīng)力對(duì)比見(jiàn)表3。根據(jù)計(jì)算結(jié)果分析,由于初始裝配的錯(cuò)位,使工作面的整個(gè)軸向錯(cuò)位變形更加嚴(yán)重,導(dǎo)致工作面擠壓應(yīng)力發(fā)生變化,工作面上部分都不接觸;同時(shí)由于軸向錯(cuò)位加劇,擠壓應(yīng)力集中程度增加。

    圖12 錯(cuò)位模型的凸肩有限元計(jì)算模型

    圖13 錯(cuò)位模型的凸肩工作擠壓應(yīng)力分布

    表3 初始裝配錯(cuò)位和未錯(cuò)位時(shí)凸肩工作面擠壓應(yīng)力對(duì)比

    根據(jù)計(jì)算對(duì)比結(jié)果,凸肩工作面的初始裝配錯(cuò)位進(jìn)一步加劇了裂紋源位置的擠壓應(yīng)力集中,使擠壓應(yīng)力從139MPa增至267MPa,增幅達(dá)到92%。

    5 凸肩工作面邊緣倒角影響

    由裂紋斷口分析發(fā)現(xiàn),凸肩工作面與上、下表面的倒角非常小,小于設(shè)計(jì)要求;同時(shí)工作面擠壓應(yīng)力分析表明在工作面下側(cè)邊緣存在擠壓應(yīng)力集中,因此,工作面邊緣倒角可能影響擠壓應(yīng)力的集中程度。

    針對(duì)此情況,并考慮檢查發(fā)現(xiàn)的初始裝配錯(cuò)位情況,建立存在初始錯(cuò)位并且無(wú)工作面邊緣倒角情況下的凸肩接觸計(jì)算模型,如圖14所示。根據(jù)模型計(jì)算得到凸肩工作面擠壓應(yīng)力分布如圖15所示,與有倒角情況下凸肩工作面擠壓應(yīng)力對(duì)比見(jiàn)表4。從應(yīng)力分布和擠壓應(yīng)力對(duì)比可見(jiàn),當(dāng)凸肩工作面邊緣無(wú)倒角時(shí),最大擠壓應(yīng)力從267MPa增至336MPa,增幅為26%,擠壓應(yīng)力集中程度進(jìn)一步加劇。

    圖14 錯(cuò)位并無(wú)倒角模型的凸肩有限元計(jì)算模型

    圖15 錯(cuò)位并無(wú)倒角模型的凸肩工作擠壓應(yīng)力分布

    根據(jù)計(jì)算對(duì)比結(jié)果,存在裂紋的凸肩工作面邊緣倒角較小會(huì)導(dǎo)致裂紋起源位置的擠壓應(yīng)力集中進(jìn)一步加劇,最高增幅達(dá)26%。

    6 結(jié)論

    根據(jù)檢查和分析結(jié)果,得出凸肩工作面出現(xiàn)裂紋的原因如下:

    (1)凸肩工作面上的裂紋為疲勞裂紋,疲勞起源于靠近凸肩下表面轉(zhuǎn)接區(qū)的凸肩工作面表面;

    (2)在風(fēng)扇葉片工作載荷和凸肩初始緊度共同作用下,在凸肩工作面長(zhǎng)、短軸方向發(fā)生變形錯(cuò)位,導(dǎo)致凸肩工作面接觸不均勻,在裂紋源附近位置產(chǎn)生一定的擠壓應(yīng)力集中;

    (3)凸肩工作面的初始裝配錯(cuò)位導(dǎo)致裂紋源附近擠壓應(yīng)力集中程度增加;

    (4)凸肩工作面與下表面的倒角較小,進(jìn)一步加劇了裂紋源附近擠壓應(yīng)力集中程度;

    (5)較大的擠壓應(yīng)力集中導(dǎo)致裂紋源部位凸肩表面耐磨層出現(xiàn)微裂紋,擴(kuò)展形成初始裂紋。

    綜上所述,初始裂紋產(chǎn)生的主要原因?yàn)橥辜绻ぷ髅娴木植繑D壓應(yīng)力過(guò)大,而工作面擠壓應(yīng)力集中主要由于風(fēng)扇工作載荷和凸肩接觸面的初始緊度共同作用下產(chǎn)生的凸肩變形錯(cuò)位導(dǎo)致的,為了避免類似情況再次發(fā)生,可以進(jìn)行預(yù)變形設(shè)計(jì),增大接觸面積,減小擠壓應(yīng)力集中程度并使應(yīng)力集中位置位于工作面中部;同時(shí)保證凸肩裝配精度和凸肩工作面邊緣的倒角尺寸,防止擠壓應(yīng)力集中進(jìn)一步增加;并且提高耐磨層的抗擠壓能力。

    [1]黃致建,郝艷華,張峻峰,等.等強(qiáng)度低損失阻尼凸肩的設(shè)計(jì)與應(yīng)用[J].航空動(dòng)力學(xué)報(bào),2007,22(6):1002-1006.HUANGZhijian,HAOYanhua,ZHANG Junfeng,etal.Design and application of the damp shroud featuring equal strength and small loss[J].Journal of Aerospace Power,2007,22(6):1002-1006.(in Chinese)

    [2]黃致建.某種低損失阻尼凸肩的結(jié)構(gòu)設(shè)計(jì)[J].航空發(fā)動(dòng)機(jī),1995(4):25-32.HUANG Zhijian.Design of the damp shroud small loss[J].Aeroengine,1995(4):25-32.(in Chinese)

    [3]李宏新,趙開(kāi)寧,張連祥.航空發(fā)動(dòng)機(jī)風(fēng)扇葉片凸肩的結(jié)構(gòu)設(shè)計(jì)[J].航空發(fā)動(dòng)機(jī),2004,30(4):14-17.LI Hongxin,ZHAO Kaining,ZHANG Lianxiang.Structural design of mid-span shrouds for aircraft engine fan blades[J].Aeroengine,2004,30(4):14-17.(in Chinese)

    [4]李辛毅,季葆華,孟慶集.采用整體圍帶或翅翼形凸肩提供摩擦阻尼的葉片典型結(jié)構(gòu)[J].汽輪機(jī)技術(shù),1996,38(4):217-220.LI Xinyi,JI Baohua,MENG Qingji.Typical structure of blade using integral shroud or flapping-wing shroud to damping[J].1996,38(4):217-220.(in Chinese)

    [5]許鍔俊,田忠賢.論葉片摩擦阻尼減振[J].航空發(fā)動(dòng)機(jī),2000,26(1):36-40.XU Ejun,TIAN Zhongxian.Friction damper of blade[J].Aeroengine,2000,26(1):36-40.(in Chinese)

    [6]許鍔俊.凸肩(葉冠)接觸面摩擦機(jī)理及參數(shù)測(cè)定的實(shí)驗(yàn)研究[J].航空發(fā)動(dòng)機(jī),2002,28(1):26-32.XU Ejun.Experimental investigation of friction mechanism and parameters measurement in interfaces of blade shrouds[J].Aeroengine,2002,28(1):26-32.(in Chinese)

    [7]郝燕平,朱梓根.摩擦阻尼器參數(shù)對(duì)葉片振動(dòng)響應(yīng)的影響[J].航空發(fā)動(dòng)機(jī),2007,33(2):18-21.HAO Yanping,ZHU Zigen.Effects of friction damper parameters on blade vibration response[J].Aeroengine,2007,33(2):18-21.(in Chinese)

    [8]單穎春,朱鋅根,劉獻(xiàn)棟.凸肩結(jié)構(gòu)對(duì)葉片的干摩擦減振研究-規(guī)律分析[J].航空動(dòng)力學(xué)報(bào),2006,21(1):174-180.SHAN Yingchun,ZHU Zigen,LIU Xiandong.Investigation of the vibration control by frictional constraints between blade shrouds-theoretical method[J].Journal of Aerospace Power,2006,21(1):174-180.(in Chinese)

    [9]單穎春,朱鋅根,劉獻(xiàn)棟.凸肩結(jié)構(gòu)對(duì)葉片的干摩擦減振研究-理論方法[J].航空動(dòng)力學(xué)報(bào),2006,21(1):168-173.SHAN Yingchun,ZHU Zigen,LIU Xiandong.Investigation of the vibration control by frictional constraints between blade shrouds theoretical method[J].Journal of Aerospace Power,2006,21(1):168-173.(in Chinese)

    [10]郝燕平,朱梓根.摩擦阻尼器參數(shù)對(duì)葉片振動(dòng)響應(yīng)的影響[J].航空發(fā)動(dòng)機(jī),2007,33(2):18-21.HAO Yanping,ZHU Zigen.Effects of friction damper parameters on blade vibration response[J].Aeroengine,2007,33(2):18-21.(in Chinese)

    [11]單穎春,朱梓根,劉獻(xiàn)棟.某帶凸肩發(fā)動(dòng)機(jī)風(fēng)扇葉片非線性響應(yīng)及斷裂故障分析[J].振動(dòng)與沖擊,2007,26(1):95-99.SHAN Yingchun,ZHU Zigen,LIU Xiandong.Investigation of nonlinear responses and crack fault of an engine fan blade with shrouds[J].Journal of Vibration and Shock,2007,26(1):95-99.(in Chinese)

    [12]李宏新,黃致建,張力,等.一種排除帶凸肩風(fēng)扇葉片榫頭故障的新方法[J].航空發(fā)動(dòng)機(jī),2002,28(2):27-31.LI Hongxin,HUANG Zhijian,ZHANG Li,et al.New technique for eliminating fault of the tenon of fan blade with mid-span shroud[J].Aeroengine,2002,28(2):27-31.(in Chinese)

    [13]張海洋,蔚奪魁,王相平.凸肩摩擦減振效果影響因素研究[J].航空發(fā)動(dòng)機(jī),2013,39(5):22-25.ZHANG Haiyang,YU Duokui,WANG Xiangping.Influence factors of vibration control by frictional constraints between blade shrouds[J].Aeroengine,2013,39(5):22-25.(in Chinese)

    [14]李洋,佟文偉,韓振宇,等.發(fā)動(dòng)機(jī)引氣管卡箍斷裂原因分析[J].失效分析與預(yù)防,2013,8(3):167-172.LI Yang,TONG Wenwei,HAN Zhengyu,et al.Fracture analysis of aeroengine cited trachea clamp [J].Failure Analysis and Prevention,2013,8(3):167-172.(in Chinese)

    [15]李洋,佟文偉,欒旭,等.發(fā)動(dòng)機(jī)低壓渦輪工作葉片裂紋失效分析[J].沈陽(yáng)航空航天大學(xué)學(xué)報(bào),2013,30(3):1-5,48.LI Yang,TONG Wenwei,LUAN Xu,et al.Fracture analysis of aeroengine low pressure turbine blade[J].Journal of Shenyang Aerospace University,2013,30(3):1-5,48.(in Chinese)

    [16]秦飛,李英,張洪濤.帶鰭狀凸肩結(jié)構(gòu)汽輪機(jī)葉片的接觸分析[J].北京工業(yè)大學(xué)學(xué)報(bào),2005,31(6):566-569.QIN Fei,LI Ying,ZHANG Hongtao.Contact analysis of steam turbine blades with fin heaves[J].Journal of Beijing Polytechnic University,2005,31(6):566-569.(in Chinese)

    [17]王親猛,張錦,馬曉秋,等.干摩擦面上接觸應(yīng)力分布的混合分析法[J].中國(guó)航空學(xué)報(bào)(英文版),2002,15(2):77-81.WANG Qinmeng,ZHANG Jin,MA Xiaoqiu,et al.Hybrid method to analyze contact stress distribution on dry friction interfaces[J].Chinese Journal of Aeronautics,2002,15(2):77-81.(in Chinese)

    [18]徐凱,李建華,陳國(guó)定,等.航空發(fā)動(dòng)機(jī)帶凸肩葉片的強(qiáng)度計(jì)算方法研究[J].現(xiàn)代制造工程,2005(2):61-63.XU Kai,LI Jianhua,CHEN Guoding,et al.Analysis of strength calculation method of aeroengine shrouded blades[J].Machinery Manufacturing Engineer,2005(2):61-63.(in Chinese)

    猜你喜歡
    磨痕倒角錯(cuò)位
    關(guān)于J.Warburton方法的鋼絲交叉磨痕體積計(jì)算及誤差分析
    箱梁三面倒角的作圖方法
    有趣的錯(cuò)位攝影
    水介質(zhì)下打磨磨痕對(duì)鋼軌疲勞損傷的影響
    減少GDX2包裝機(jī)組“磨痕”煙包的數(shù)量
    關(guān)于機(jī)械設(shè)計(jì)與制造中的零件倒角問(wèn)題分析
    基于機(jī)械設(shè)計(jì)與制造中的零件倒角研究
    采用倒角坯改善鑄坯邊直裂缺陷的研究
    上海金屬(2016年2期)2016-11-23 05:34:40
    從輪胎磨痕準(zhǔn)確判斷裝備輪胎損傷
    避免“錯(cuò)位相減,一用就錯(cuò)”的錦囊妙計(jì)
    胶州市| 额敏县| 吐鲁番市| 浮梁县| 琼中| 景谷| 大港区| 格尔木市| 聂拉木县| 赫章县| 常熟市| 上饶市| 新邵县| 罗山县| 工布江达县| 内江市| 封丘县| 屯昌县| 磐安县| 柳林县| 屏东市| 周至县| 肇庆市| 古丈县| 津市市| 紫金县| 岳阳市| 聂拉木县| 海林市| 盖州市| 白水县| 长垣县| 富宁县| 长乐市| 彰化市| 永泰县| 英山县| 河北区| 永清县| 邵阳市| 原平市|