• 
    

    
    

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

      ?

      鋼軌無損檢測中的超聲導波技術(shù)

      2015-05-25 10:20:36周建民徐清瑤廖曉蘇
      儀表技術(shù)與傳感器 2015年6期
      關(guān)鍵詞:群速度軌頭導波

      周建民,徐清瑤,李 鵬,萬 琪,廖曉蘇

      (華東交通大學機電工程學院,江西南昌 330013)

      0 引言

      鋼軌在使用過程中,由于自然因素以及列車載荷的作用,致使其表面和內(nèi)部易產(chǎn)生各類損傷和缺陷,嚴重時甚至會造成鋼軌斷裂、列車脫軌等重大事故[1]。因此,開展鋼軌無損檢測技術(shù)的研究具有重要的意義。

      超聲波檢測是鋼軌無損檢測的重要手段之一。目前,常規(guī)超聲波檢測方法有脈沖反射法、穿透法和共振法等[2]。這些方法在技術(shù)方面比較成熟,但其檢測過程采用逐點掃描方式,檢測速度慢,已較難滿足現(xiàn)階段我國鐵路軌道高速發(fā)展的需求。與常規(guī)超聲波檢測方法相比,新型的超聲導波技術(shù)具有沿傳播路徑衰減小、傳播距離遠、檢測速度快等優(yōu)點,已被廣泛應(yīng)用于管道等大型構(gòu)件的無損檢測中[3],并且該技術(shù)在鋼軌大范圍缺陷的完整性檢測上具有良好的應(yīng)用前景。因此,鋼軌超聲導波檢測技術(shù)已成為近年來國內(nèi)外無損檢測領(lǐng)域研究的熱點之一。

      1 導波檢測原理

      機械振動在彈性介質(zhì)中的波稱為彈性波(即聲波),而頻率大于20 kHz的聲波則稱為超聲波。當超聲波被局限在具有邊界的介質(zhì)(如平板、管道等)內(nèi)傳播時,超聲波將會在邊界處不斷地反射,從而沿著介質(zhì)的方向傳播形成超聲導波[4-5]。產(chǎn)生超聲導波的這種具有邊界的介質(zhì)稱為波導,軌頭、軌腰和軌底是鋼軌中三種典型的波導。由于這些波導的結(jié)構(gòu)不同,在其中傳播的導波頻率和速度也相應(yīng)發(fā)生變化。

      用超聲導波對鋼軌進行無損檢測時,可以通過信號發(fā)生裝置產(chǎn)生激勵信號,經(jīng)功率放大器放大后由導波傳感器在鋼軌的一端激發(fā)超聲導波。如果導波沿著沒有損傷的軌頭、軌腰和軌底傳播,那么導波的群速度和相速度就基本保持一致。如果導波在傳播過程中遇到界面不連續(xù)處(如軌底腐蝕),則可能發(fā)生反射、散射和模式轉(zhuǎn)換,這樣便會產(chǎn)生攜帶局部缺陷特征的回波[6],如圖1所示。通過對回波信號進行分析就可以確定缺陷的位置,回波幅值還能夠用于鋼軌損傷程度的評定。

      圖1 超聲導波檢測原理(以腐蝕為例)

      2 導波的頻散及多模態(tài)特性

      相對于體波(常見的超聲波),鋼軌中超聲導波的顯著特點就是具有頻散特性。導波在鋼軌軌頭、軌腰和軌底中傳播時,由于受波導截面形狀和尺寸的影響,使得在鋼軌中傳播的超聲波速度依賴于波的頻率,從而導致超聲波的幾何彌散,即導波的相速度隨頻率的不同而發(fā)生變化,這種現(xiàn)象就稱為鋼軌中超聲導波的頻散現(xiàn)象[7]。頻散對超聲導波的傳播有很大影響,主要表現(xiàn)為導波信號在傳播一定距離后時域波包的寬度會拉長,幅度明顯減小,這不僅降低了檢測的靈敏度,而且會給后續(xù)信號處理帶來很大的困難。因此,應(yīng)盡量避免將頻散嚴重的導波模態(tài)用于鋼軌的無損檢測。

      多模態(tài)是鋼軌中超聲導波的另一重要特性。波導結(jié)構(gòu)一般在某一頻率下至少存在兩個或兩個以上的導波模態(tài),它們在波導中以不同的群速度傳播。鋼軌中的超聲導波存在很多種不同的模態(tài),多種模態(tài)疊加在一起并分別以不同的速度傳播可能會使數(shù)據(jù)信號解析變得困難[8]。因此,需要合理選取合適的模態(tài)來檢測鋼軌不同部位、不同類型的典型缺陷。

      3 鋼軌中超聲導波傳播特性的數(shù)值計算

      超聲導波的傳播特性(頻散、波結(jié)構(gòu)等)在很大程度上影響著其檢測效果[9]。因此,應(yīng)用導波技術(shù)對鋼軌進行無損檢測時,需要根據(jù)導波的傳播特性來選擇合適的導波模態(tài)和頻率范圍。規(guī)則截面波導(如板、管)中導波傳播特性的數(shù)值計算方法已經(jīng)很成熟[10-12],但由于鋼軌是異型截面波導,一般的商用軟件無法完成其傳播特性的數(shù)值計算,因此,鋼軌中超聲導波傳播特性的數(shù)值計算一直是導波檢測研究的熱點。

      國外很多學者應(yīng)用半解析有限元(Semi-analy-tical finite element,SAFE)法對鋼軌中超聲導波的傳播特性進行了相關(guān)研究。該方法由Taweel等[13]提出,可用于求解任意截面波導中超聲導波的傳播特性。將該方法應(yīng)用于鋼軌時,只需對鋼軌橫截面進行離散,沿波傳播方向上的位移用諧波指數(shù)函數(shù)eikx(x為鋼軌長度,k為導波的波數(shù))來描述,從總體運動方程特征值系統(tǒng)中解得特征值ki(即,第i個共振模式的波數(shù))后,便可獲得頻散關(guān)系。

      Hayashi等[15]采用半解析有限元法計算日制JIS50型鋼軌中超聲導波的頻散曲線,得到了導波在0~100 kHz范圍內(nèi)的相速度、群速度頻散曲線,分別如圖2和圖3所示(其中,橫坐標均為激勵頻率,縱坐標分別為鋼軌中導波隨激勵頻率變化的相速度和群速度)。從圖中可以看出,鋼軌在該頻率范圍內(nèi)存在眾多傳播模式,有些模式的導波頻散曲線在一定的頻率范圍內(nèi)較為平緩,展現(xiàn)了較好的傳播特性,可用于檢測。為了驗證數(shù)值計算結(jié)果,Hayashi等分別對軌頭踏面和軌頭側(cè)面進行導波檢測實驗,實驗得到的相速度、群速度頻散曲線與理論得到的結(jié)果非常吻合。

      圖2 鋼軌相速度頻散曲線

      圖3 鋼軌群速度頻散曲線

      同樣采用半解析有限元法對鋼軌中超聲導波的傳播特性進行 數(shù) 值 計 算 的 還 有 Coccia[16]、Bartoli[17]、Loveday[18-19]、Lee[20]、Hesse[21-22]、盧超[23-24]以及??〗埽?5]等人,Hayashi[26]還結(jié)合半解析有限元法開發(fā)了一款用于計算鋼軌中超聲導波頻散曲線的軟件Rail Dispersion。

      除了半解析有限元法,一些學者也采用其他方法來解決鋼軌中導波傳播特性的數(shù)值計算問題。Sanderson等[27]根據(jù)固有頻率分析可以提供任意頻率下結(jié)構(gòu)所有傳播模式的位移模態(tài),采用三維固有頻率分析方法求出了英制BS113A型鋼軌在0~60 kHz范圍內(nèi)的頻散曲線。國內(nèi)何存富等[28]根據(jù)結(jié)構(gòu)振動解與波動解的互化原理,采用振動模態(tài)分析方法提取出鋼軌中導波的相速度和群速度頻散曲線,并由鋼軌振動分析中獲得的變形信息求解出導波的波結(jié)構(gòu)。

      以上研究進展有利地推動了鋼軌超聲導波檢測技術(shù)的發(fā)展,計算出的導波傳播特性能為鋼軌檢測過程中激勵頻率及導波模態(tài)的選擇提供依據(jù)。

      4 激勵頻率及導波模態(tài)的選擇

      通過對鋼軌中超聲導波傳播特性的數(shù)值計算表明,鋼軌中存在很多不同種類的導波模態(tài),且隨著頻率增大鋼軌中的導波模態(tài)也會增多,因此,在用導波技術(shù)對鋼軌進行無損檢測時應(yīng)選取合適的激勵頻率及導波模態(tài)。

      Gharaibeh等[29]通過ABAQUS建立三維有限元模型進行數(shù)值模擬,得到了英制BS113A型鋼軌中超聲導波的8種模態(tài),如圖4所示。對這些模態(tài)進行分析,發(fā)現(xiàn)F3、T2、F2模態(tài)的導波能量都分別集中于單一部位(軌頭、軌腰或軌底),且在70 kHz時這三種模態(tài)的導波相速度基本一致。因此,分別選用70 kHz時F3、T2、F2模態(tài)的導波對鋼軌缺陷進行檢測,結(jié)果成功檢測出離激勵點9 m遠處2 mm深的軌頭橫向缺陷和8 m遠處5 mm深的軌底橫向缺陷,從而驗證了激勵頻率及導波模態(tài)選擇的正確性。

      Hayashi等[30]指出可將鋼軌軌底主要的導波模態(tài)分為縱向振動模態(tài)、垂直振動模態(tài)和橫向振動模態(tài),縱向振動模態(tài)的導波群速度頻散曲線總體上很雜亂且無規(guī)律,不適合用于檢測,而橫向振動模態(tài)和垂直振動模態(tài)的導波在60~200 kHz頻率范圍內(nèi)群速度頻散曲線較為平緩,可用于檢測。

      Lee等[21]指出對鋼軌軌頭剝離不敏感,而對軌頭剝離下方的橫向裂紋敏感的頻率和模態(tài)能用于檢測軌頭剝離下方的橫向裂紋,并得出結(jié)論:低頻導波可以更好地檢測位于鋼軌軌頭的橫向裂紋,且能最大限度地減小剝離帶來的影響。

      圖4 鋼軌不同部位的導波模態(tài)

      Rose等[31]認為適合鋼軌檢測的頻率范圍為40~80 kHz,并用電磁超聲傳感器激勵60 kHz的導波進行實驗驗證。實驗表明,該方法可以有效地檢測鋼軌的橫向裂紋,但信噪比較低。

      以上研究成果表明,選擇激勵頻率和導波模態(tài)時,應(yīng)保證在所選頻率范圍內(nèi)所選模態(tài)的導波傳播能力強,對缺陷具有較高的靈敏度,同時保證該導波的頻散曲線較為平緩,以避開頻散嚴重的檢測頻帶,并且盡可能選擇較低的激勵頻率,以減少鋼軌中超聲導波模態(tài)的數(shù)量。

      5 導波信號處理方法

      應(yīng)用超聲導波對鋼軌進行檢測時,對于給定的激勵頻率至少會包含2種模態(tài),并且隨著頻率的增大模態(tài)還會增多,只有通過信號處理方法來區(qū)分接收信號的不同模態(tài)才能獲得導波檢測所需的模態(tài)。此外,實驗獲得的檢測信號還可能受到實驗環(huán)境引入的噪聲等因素的干擾,這些干擾會影響甚至淹沒檢測中的有用信號,因此,需要用信號處理方法來增強目標信號[32]。

      目前,常用的信號處理方法主要有時域分析、頻域分析、時頻分析以及小波變換等。由于鋼軌中的導波信號同時受時間和頻率的限制,具有典型的非平穩(wěn)隨機信號特征,若用單一的時域方法或頻域方法則難以達到良好的處理效果,而時頻分析是瞬時頻譜分析,因此,用時頻分析方法處理可以獲得更豐富、更準確的信息。

      Scale等[33]采用基于小波變換的時頻分析方法來表征鋼軌中具有多模態(tài)及頻散特征的導波的傳播。對長為4.6 m的美制115-lb AREMA型鋼軌進行導波檢測實驗后,通過時頻分析方法識別出了1~7 kHz范圍內(nèi)的3種縱向振動模態(tài)、一種橫向振動模態(tài)以及一種垂直振動模態(tài)的導波,并采用Gabor小波變換從檢測信號的小波尺度圖中提取出了導波的群速度頻散曲線和衰減系數(shù)曲線。同樣采用時頻分析方法對鋼軌中導波檢測信號處理的還有 McNamara[34]、Scalea[35]以及 Bartolia[36]等人。

      針對超聲導波檢測時檢測靈敏度和分辨率受鋼軌中頻散波模態(tài)影響較大的問題,Kappatos等[37]提出了一種基于離散小波變換的信號處理方法,并通過添加信噪比為-20~20 dB的頻散波形來對這種信號處理方法進行評估,結(jié)果表明,該方法能在很大程度上減少頻散波帶來的影響。

      6 導波技術(shù)在鋼軌無損檢測中的工程應(yīng)用

      目前,國外已有一些研究機構(gòu)研制出了鋼軌導波檢測設(shè)備[38-42],并將其應(yīng)用于實際線路的鋼軌無損檢測。Gurvich等[38]研制了型號為AKR1224的便攜式低頻超聲導波檢測儀,如圖5所示。該儀器將由12個超聲波傳感器組成的無線陣列式低頻導波檢測探頭作為其發(fā)射和接收單元,并采用脈沖回波法進行檢測,無需掃描整個鋼軌表面就能快速尋找并定位鋼軌中一定距離的缺陷,且能分別對軌頭、軌腰和軌底缺陷進行檢測。此外,該導波檢測儀還內(nèi)置存儲器,用于存放測量結(jié)果。大量的現(xiàn)場測試表明,該導波檢測儀最長檢測距離為30 m,缺陷定位精度能達到0.1 m,且鋼軌不同部位的缺陷檢測不會相互干擾,檢測性能良好。

      圖5 便攜式低頻超聲導波檢測儀

      Wilcox等[39]研制了主要用于檢測鋼軌垂直裂紋和鋁熱焊縫缺陷的G-Scan鋼軌超聲導波檢測裝置,該裝置通過一個機械/氣動控制系統(tǒng)將傳感器陣列干耦合到鋼軌上來進行檢測,且能控制檢測模式和傳播方向,滿足鐵路軌道快速檢測的技術(shù)指標。

      Scalea等[40]研制了基于超聲導波的非接觸式軌檢車,該設(shè)備采用激光激勵超聲導波,并用空氣耦合式換能器來接收,能快速檢測出水平裂縫下的垂直裂紋。此外,該設(shè)備還通過離散小波變換來提高非接觸檢測的信噪比。

      由以上研究成果得知,超聲導波技術(shù)應(yīng)用于鋼軌時檢測速度快且能保證靈敏度,檢測性能良好,滿足現(xiàn)代鐵路軌道高速發(fā)展的需求。

      7 結(jié)束語

      超聲導波技術(shù)作為一種新型的無損檢測技術(shù),檢測速度快,可以避免逐點檢測的費時和高成本,在鋼軌的長距離檢測中有著明顯的優(yōu)勢。目前,鋼軌超聲導波檢測技術(shù)在以下幾個方面仍有待進一步研究。

      (1)近年來國內(nèi)外多采用半解析有限元法來求解鋼軌中超聲導波的頻散曲線,探索新的頻散曲線求解方法對于激勵頻率及導波模態(tài)的選擇將會更有幫助。

      (2)導波傳播過程中存在很多種不同的模態(tài),不排除其他模態(tài)的導波對鋼軌無損檢測會有更佳的效果。此外,還可以選擇不同的導波模態(tài)對鋼軌同一部位缺陷進行檢測,通過特定數(shù)據(jù)融合方法,或許能更為全面地反映實際缺陷的特征。

      (3)鋼軌缺陷有很多種不同的類型,而目前鋼軌超聲導波檢測技術(shù)在檢測信號智能分析方面還鮮有文獻報道,因此,很有必要選擇合適的信號處理方法對檢測信號進行準確的特征提取并實現(xiàn)檢測信號類型的智能識別。

      [1]任遠,張友鵬,史宏章,等.基于LabVIEW的超聲導波斷軌檢測系統(tǒng)的設(shè)計.鐵道運營技術(shù),2010,16(2):16-19.

      [2]何學科.鋼軌探傷.北京:中國鐵道出版社,2009:34-35.

      [3]王軍陣,王建斌,張軒碩.基于ARM和Linux的超聲導波管道檢測系統(tǒng).儀表技術(shù)與傳感器,2011(4):51-52.

      [4]ROSE J L.Ultrasonic waves in solid media.Cambridge,England:Cambridge University Press,1999:82 -92.

      [5]丁輝.計算超聲學—聲場分析及應(yīng)用.北京:科學出版社,2010:102-132.

      [6]徐柳娟,王秋萍.導波技術(shù)在管道損傷檢測中的應(yīng)用研究評述.人民黃河,2013,35(3):124 -126.

      [7]金建華,金紀東.基于磁致伸縮效應(yīng)的管道縱向超聲導波檢測傳感器.儀表技術(shù)與傳感器,2004(11):3-4.

      [8]王雪梅,倪文波,王平.高速鐵路軌道無損探傷技術(shù)的研究現(xiàn)狀和發(fā)展趨勢.無損檢測,2013,35(2):10-17.

      [9]焦敬品,鐘茜,王炯耿,等.窄板中超聲導波傳播特性試驗研究.工程力學,2013,30(7):255 -261.

      [10]閻石,張海鳳,蒙彥宇.Lamb波頻散曲線的數(shù)值計算及試驗驗證.華中科技大學學報(城市科學版),2010,27(1):1-4.

      [11]何存富,李隆濤,吳斌.超聲導波在管道中傳播的數(shù)值模擬.北京工業(yè)大學學報,2004,30(2):129 -133.

      [12]鄭祥明,趙玉珍,史耀武.蘭姆波頻散曲線的計算.無損檢測,2003,25(2):66 -68.

      [13]TAWEEL H,DONG S B,KAZIC M.Wave reflection from the free end of a cylinder with an arbitrary cross-section.International Journal of Solids and Structures,2000,37(12):1701 -1726.

      [14]HAYASHI T,ROSE J L.Guided wave simulation and visualization by a Semi- Analytical Finite Element Method.Materials Evaluation,2003,61(1):75 -79.

      [15]HAYASHI T,SONG W J,ROSE J L.Guided wave dispersion curves for a bar with an arbitrary cross-section,a rod and rail example.Ultrasonics,2003,41:175 -183.

      [16]COCCIA S,BARTOLI I,MARZANI A,et al.Numerical and Experimental Study of Guided Waves for Detection of Defects in the Rail Head.NDT&E International,2011,44:93 -100.

      [17]BARTOLI I,MARZANI A,SCALEAA F L,et al.Modeling wave propagation in damped waveguides of arbitrary cross-section.Journal of Sound and Vibration ,2006,295:685 -707.

      [18]LOVEDAY P W.Measurement of modal amplitudes of guided waves in rails.The International Society for Optical Engineering,2008,6935(69351J):1-8.

      [19]LOVEDAY P W.Guided wave inspection and monitoring of railway track.J Nondestruct Eval,2012,31:303 -309.

      [20]LEE C M,ROSE J L,CHO Y.A guided wave approach to defect detection under shelling in rail.NDT&E International,2009,42(3):174-180.

      [21]HESSE D,CAWLEY P.Surface wave modes in rails.Journal of the Acoustical Society of America,2006,120(2):733 -740.

      [22]HESSE D,CAWLEY P.A single probe spatial averaging technique for guided waves and its applica- tion to surface wave rail inspection.IEEE Transactions on Ultrasonics,2007,54(11):2344 -2356.

      [23]LU C,MEN P,LI L L.An experimential study EMAT ultrasonic surface waves modes in rail head.International Joural of Applied E-lectromagnetics and Mechanics,2010,33:1127 -1133.

      [24]盧超,李誠,??〗?鋼軌軌底垂直振動模式導波檢測技術(shù)的實驗研究.實驗力學,2012,27(5):593 -600.

      [25]CHANG J J.Nondestructive Evaluation for Bottom Edges of a Rail using Guided Wave.2011年全球華人無損檢測技術(shù)高峰論壇,廈門,2011.

      [26]HAYASHI T,MIYAZAKI Y,MURASE M.Guided wave inspection for bottom edge of rails.Review of Progress in Quantitative in NDE,2007,26:169 -175.

      [27]SANDERSON R M,SMITH S D.The Application of finite element modelling to guided wave testing systems.Insight:Non-Destructive Testing and Condition Monitoring,2002,44(6):359 -363.

      [28]何存富,劉青青,焦敬品,等.基于振動模態(tài)分析的鋼軌中超聲導波傳播特性數(shù)值計算方法.振動與沖擊,2014,33(3):9-13.

      [29]GHARAIBEH Y,SANDERSON R,MUDGE P,et al.Investigation of the behavior of selected ultrasonic guided wave modes to inspect rails for long - range testing and monitoring.Proc.ImechE Vol.225 Part F:Rail and Rapid Transit,2010,44(2):311 -324.

      [30]HAYASHI T.Guided wave dispersion curves drived with a semianalytical finite element method and its application to nondestructive inspection.Japanese Journal of Applied Physics ,2008,47(5):3865-3870.

      [31]ROSE J L,AVIOLI M J,MUDGE P,et al.Guided wave inspection potential of defects in rail.NDT&E International,2004,37:153 -161.

      [32]吳斌,鄧菲,何存富.超聲導波無損檢測中的信號處理研究進展.北京工業(yè)大學學報,2007,33(4):342-348.

      [33]SCALE F L,MCNAMARA J.Measuring high - frequency wave propagation in railroad tracks by joint time-frequency analysis.Journal of Sound and Vibration,2004,273:637 -651.

      [34]MCNAMARA J.Health monitoring of railroad tracks by elastic -wave based on non-destructive testing.PhD Dissertation,San Diego:Unive-rsity of California,2003.

      [35]SCALEA L,MCNAMARA J.Ultrasonic NDE of railroad tracks:air-coupled cross-sectional inspection and long-range inspection.Insight,2003,45(6):394 –401.

      [36]BARTOLI I,SCALEAA F L,F(xiàn)ATEH M,et al.Modeling guided wave propagation with application to the long-range defect detection in rail- road tracks.NDT&E International,2005,38:325 -334.

      [37]KAPPATOS V,SELUCK C,GAN T H,et al.Decreasing the influence of dispersive wave modes in long-range ultras- onic rail testing using wavelets.5th IET Conference on Railway Condition Monitoring and Non-Destructive Testing,Derby,United kingdom,2011.

      [38]GURVICH A K,SAMOKRUTOV A A,SHEVALDYKIN V G.Guided wave ultrasonic flaw detection in rails Application Experience[EB/OL].(2011-08-10)[2014 -04 -27].http://www.ndt.net/article/ecndt 2006/doc/We.4.5.4.pdf.

      [39]WILCOX P D,PAVLAKOVIC B N,EVANS M J,et al.Long range inspection of rail using guided waves.Review of Progress in Quantitative in NDE,2003,22:236 -243.

      [40]SCALEA L,RIZZO F,COCCIA P,et al.Non - contact ultrasonic inspection of rails and signal processing for automatic defect detection and classification.Insight,2005,47(6):1–8.

      [41]COCCIA S.Ultrasonic guided waves for structural health monitoring and application to rail inspection prototype for the Federal Railroad Administration .San Diego:University of California,2007.

      [42]MARIANI S,NGUYEN T,PHILLIPS R R,et al.Noncontact ultrasonic guided wave inspection of rails.Structural Health Monitoring,2013(12):539-548.

      猜你喜歡
      群速度軌頭導波
      VTI介質(zhì)地震波群速度三維擴展各向異性線性近似表征
      超聲導波技術(shù)在長輸管道跨越段腐蝕檢測中的應(yīng)用
      重載鐵路直線段鋼軌斷裂原因分析
      鐵道建筑(2021年11期)2021-03-14 10:01:50
      激光調(diào)制摻鉺光纖中光速及其在高靈敏光纖傳感領(lǐng)域的潛在應(yīng)用研究
      卷簧缺陷檢測的超聲導波傳感器研制
      電子制作(2019年9期)2019-05-30 09:42:00
      物質(zhì)波的波速與頻率公式中的能量
      小議超光速
      科教導刊(2018年11期)2018-07-25 11:34:48
      基于Gocator視覺傳感器的軌頭參數(shù)計算
      淬火引軌早期失效分析和預防對策
      磁致伸縮導波激勵傳感器模型及輸出特性
      泾源县| 琼中| 西乌珠穆沁旗| 大连市| 正安县| 九龙县| 长海县| 夏津县| 天长市| 宕昌县| 青岛市| 宜春市| 湘潭市| 清水县| 微山县| 徐汇区| 蓬莱市| 宕昌县| 盈江县| 高平市| 青冈县| 当涂县| 翁牛特旗| 拜泉县| 同心县| 怀远县| 托克托县| 乡宁县| 巴中市| 许昌市| 仙游县| 左云县| 秦皇岛市| 西乌| 伊春市| 开化县| 双城市| 习水县| 凤山县| 襄城县| 陆川县|