• 
    

    
    

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

      ?

      復(fù)合式低頻隔振器的理論建模及其性能分析

      2017-11-06 02:29:41王林翔
      振動(dòng)與沖擊 2017年19期
      關(guān)鍵詞:振源固有頻率諧振

      吳 庭,王林翔

      (浙江大學(xué) 流體動(dòng)力及機(jī)電系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,杭州 310027)

      復(fù)合式低頻隔振器的理論建模及其性能分析

      吳 庭1,王林翔2

      (浙江大學(xué) 流體動(dòng)力及機(jī)電系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,杭州 310027)

      抑制低頻振動(dòng)是當(dāng)前振動(dòng)領(lǐng)域的一個(gè)熱點(diǎn)問(wèn)題,將慣容器與液壓蓄能器結(jié)合在一起,同時(shí)利用形狀記憶合金彈簧的超彈性,設(shè)計(jì)了一種將慣容器和擬零剛度隔振結(jié)合的復(fù)合式超低頻隔振器。說(shuō)明了超低頻隔振的原理,分析了隔振器的力傳遞比,研究了慣容器慣容值以及彈簧剛度對(duì)系統(tǒng)諧振頻率的影響。計(jì)算結(jié)果表明,所設(shè)計(jì)的超低頻隔振器的諧振頻率可以降低至1 Hz以內(nèi)。

      擬零剛度;慣容器;超彈性;超低頻隔振

      對(duì)低頻振動(dòng)激勵(lì)的隔離和抑制,一直是難以解決的問(wèn)題,因?yàn)橥ǔ2捎玫某R?guī)被動(dòng)隔振,無(wú)法適應(yīng)低頻隔振的要求[1]。普通的橡膠隔振器,對(duì)激勵(lì)頻率在5~6 Hz的振動(dòng)有較好的隔振效果。而對(duì)于2 Hz以下的振動(dòng),則多采用空氣彈簧,它是利用氣體壓縮的非線性恢復(fù)力和阻尼力來(lái)緩沖沖擊和振動(dòng)[2]??諝鈴椈筛粽衿鞯膭偠群妥枘峥蛇M(jìn)行調(diào)節(jié)與控制。然而主動(dòng)振動(dòng)控制的缺點(diǎn)就是成本較高,且需要時(shí)刻保證外加電源或氣源[3]。

      被動(dòng)振動(dòng)控制由于其結(jié)構(gòu)簡(jiǎn)單,經(jīng)濟(jì)性好,可靠性高,不附加能量仍然具有較寬廣的適用性[4]。如何對(duì)現(xiàn)有的被動(dòng)隔振裝置進(jìn)行改進(jìn)而使其具有低頻甚至超低頻隔振的特性,是許多學(xué)者的關(guān)注與研究的熱點(diǎn)。 為解決低頻隔振,一類方法是利用正負(fù)剛度并聯(lián)使系統(tǒng)獲得動(dòng)態(tài)零剛度。Carrella等[5]利用四根傾斜彈簧和一根豎直彈簧構(gòu)成了正負(fù)剛度并聯(lián)的機(jī)構(gòu),理論和實(shí)驗(yàn)驗(yàn)證了其具有零剛度的特性。2011年,Le等[6]將這一正負(fù)剛度聯(lián)合機(jī)構(gòu)設(shè)計(jì)到了汽車座椅的支撐結(jié)構(gòu)中去,有效地提高了系統(tǒng)的隔振效率。劉興天等[7]設(shè)計(jì)了類似的正負(fù)剛度并聯(lián)的機(jī)構(gòu),其中負(fù)剛度特性由屈曲歐拉梁產(chǎn)生。2008年,Mizuno等[8]曾提出過(guò)利用相互吸引的永磁鐵來(lái)給系統(tǒng)提供負(fù)剛度,在此基礎(chǔ)上,Carrella等[9]設(shè)計(jì)了將這一負(fù)剛度機(jī)構(gòu)與提供正剛度的機(jī)械彈簧聯(lián)合使用的機(jī)構(gòu)。

      第二類方法是增加被動(dòng)隔振系統(tǒng)的慣性值,2002年,Smith[10]提出了慣容器的概念,可以為被動(dòng)隔振系統(tǒng)提供高于自身質(zhì)量幾十倍的視在質(zhì)量,從而有效地降低系統(tǒng)的固有頻率,達(dá)到低頻隔振的效果?,F(xiàn)有的慣容器有多種形式,包括齒輪齒條式,滾珠絲杠式,液壓式慣容器[11-13]。慣容器會(huì)產(chǎn)生一個(gè)與兩個(gè)端點(diǎn)相對(duì)加速度成正比的力,這一比例系數(shù)又稱為慣容器的慣容值。國(guó)內(nèi)Xu等[14-16]嘗試將其應(yīng)用在汽車懸架上,從而提高汽車懸架的減振特性。

      本文提出了一種將慣容器與擬零剛度特性結(jié)合的復(fù)合式低頻隔振器,并且利用了形狀記憶合金彈簧(SMA(Shape Memory Alloy) spring)的超彈性[17]。對(duì)復(fù)合式低頻隔振器的動(dòng)力學(xué)性能進(jìn)行了建模,分析了隔振器的隔振效果,并分析了系統(tǒng)結(jié)構(gòu)參數(shù)對(duì)于系統(tǒng)諧振頻率的影響。為低頻隔振器的設(shè)計(jì)提供了理論支持。

      1 復(fù)合式低頻隔振器的工作原理

      本文所提出的復(fù)合式低頻隔振器結(jié)構(gòu)簡(jiǎn)圖如圖1所示。該復(fù)合式低頻隔振器由上部慣容器和底部的液壓緩沖與支撐裝置構(gòu)成。液壓緩沖部分的頂壁有一部分是由彈性膜與其他部分相連的,SMA彈簧通過(guò)這一活動(dòng)部分與內(nèi)部油壓平衡。

      靜止?fàn)顟B(tài)下,由于液壓油的初始油壓P0作用于活塞桿底部而支撐物體,使得系統(tǒng)具有給定的靜態(tài)剛度。慣容器中置有上下兩個(gè)很軟的對(duì)中彈簧。當(dāng)振源產(chǎn)生低頻振動(dòng)時(shí),由于與振源相連的慣容器會(huì)會(huì)產(chǎn)生一個(gè)較大的視在質(zhì)量,同時(shí)系統(tǒng)的振動(dòng)通過(guò)活塞桿傳入下方的液壓緩沖與支撐裝置中,形成了系統(tǒng)的動(dòng)態(tài)低剛度。兩者相結(jié)合使系統(tǒng)獲得很低的固有頻率,從而達(dá)到低頻隔振的目的。

      圖1 復(fù)合式低頻隔振器結(jié)構(gòu)簡(jiǎn)圖Fig.1 Structural diagram of hybrid isolator

      圖1中:活塞桿的面積為A1;活動(dòng)薄膜的作用面積為A2。

      2 復(fù)合式低頻隔振器的動(dòng)態(tài)模型

      如圖1所示,假設(shè)振動(dòng)位移為x,慣容器中設(shè)有對(duì)中彈簧,總剛度為k1,慣容器的慣容值為b,液壓油預(yù)置壓力為P0。記振動(dòng)質(zhì)量為m。當(dāng)振源無(wú)振動(dòng)時(shí),系統(tǒng)處于靜止?fàn)顟B(tài),整個(gè)機(jī)構(gòu)保持靜態(tài)平衡,可得

      P0A1=mg

      (1a)

      P0A2=f(y0)

      (1b)

      式中,f(y0)為SMA彈簧的預(yù)壓力。可以通過(guò)調(diào)整螺釘來(lái)調(diào)節(jié)彈簧的預(yù)壓力。

      當(dāng)外界輸入振動(dòng)F0eiωt時(shí),機(jī)構(gòu)狀態(tài)發(fā)生變化,整體保持動(dòng)態(tài)平衡,系統(tǒng)中的慣容器在系統(tǒng)的動(dòng)態(tài)過(guò)程中會(huì)產(chǎn)生與上下油口速度差成正比的力,這一比例系數(shù)即為慣容值b,系統(tǒng)的動(dòng)態(tài)平衡方程為

      (2)

      式中:ΔP為液壓油體積變化引起的壓力變化值;b為慣容器的慣容值;ε為系統(tǒng)的阻尼系數(shù)。

      設(shè)SMA彈簧位移變化為y,慣容器中的活塞桿向下運(yùn)動(dòng)時(shí),深入油液的長(zhǎng)度變化為x,所對(duì)應(yīng)的油液體積變化為ΔV1=A1x,形狀記憶合金彈簧壓縮量變化對(duì)應(yīng)的油液體積變化為ΔV2=A2y,總的體積變化為ΔV=A1x-A2y。

      ΔP與ΔV的關(guān)系可由液壓基礎(chǔ)知識(shí)得[18]

      (3)

      ΔP=C(A1x-A2y)

      (4)

      記SMA彈簧的動(dòng)態(tài)剛度為κ,彈簧力的變化與油液壓力的變化平衡,也可以記為動(dòng)態(tài)剛度和位移y的乘積

      ΔF=A2ΔP=κy

      (5)

      將式(4)代入式(5),可以得到x和y的關(guān)系,如下式所示

      (6)

      將式(6)代入式(5),即可得到壓力ΔP變化x的關(guān)系式,如式(7)所示

      (7)

      將式(7)代入式(2),即可得到系統(tǒng)的動(dòng)力學(xué)方程

      (8)

      將式(8)進(jìn)行無(wú)量綱化處理后得

      (9)

      3 系統(tǒng)的固有頻率分析

      由振動(dòng)理論可知,系統(tǒng)的固有頻率與外界的輸入無(wú)關(guān),為了分析系統(tǒng)的固有頻率,得到系統(tǒng)的間諧運(yùn)動(dòng)的微分方程

      (10)

      式中:ξ為系統(tǒng)的阻尼比;ωn為系統(tǒng)的固有頻率,其值為

      (11)

      由線性振動(dòng)理論可知,線性振動(dòng)系統(tǒng)的固有頻率為

      (12)

      對(duì)比兩者不難發(fā)現(xiàn),系統(tǒng)的等效質(zhì)量和等效剛度不再是原有的m和k,而變成了m′和k′,分別為

      m′ =m+b;

      (13)

      為了獲得很低的系統(tǒng)固有頻率,增加系統(tǒng)的等效質(zhì)量和獲得擬零剛度就成了兩個(gè)最主要的方法。由式(13)看出,慣容器的存在增加了系統(tǒng)的等效質(zhì)量。同時(shí)通過(guò)設(shè)計(jì),使得等效剛度接近于零,這樣就可以達(dá)到最終的低頻隔振的效果。

      式(13)中,κ為SMA彈簧的動(dòng)態(tài)剛度,不同于普通的機(jī)械彈簧,由于SMA彈簧中超彈性的存在,彈簧的力-位移曲線不再是線性曲線,而是呈現(xiàn)非線性的滯回曲線,本文選用文獻(xiàn)[19]中的SMA彈簧,力-位移曲線如圖2所示。

      圖2 SMA彈簧的力-位移曲線Fig.2 Force-displacement curve of SMA spring

      根據(jù)SMA彈簧的力-位移關(guān)系,為了讓SMA彈簧工作在超彈性區(qū)域,可以通過(guò)調(diào)整螺釘調(diào)整SMA彈簧的壓縮量使其工作在A點(diǎn)附近,此時(shí)SMA彈簧的動(dòng)態(tài)剛度最小,κ=10 N/mm。

      A1,A2的選取和初始條件有關(guān)。A1的選取與隔振目標(biāo)的重量及初始?jí)毫τ嘘P(guān),由式(1a)選定A1的面積。通過(guò)調(diào)節(jié)油壓可以為不同振源質(zhì)量的物體提供系統(tǒng)靜態(tài)剛度。根據(jù)圖2,將SMA彈簧的預(yù)壓力設(shè)為50 N。根據(jù)式(1b),同時(shí)可以確定薄膜的面積A2。

      由于系統(tǒng)中提供靜態(tài)剛度的不再是彈簧,因此對(duì)中彈簧的剛度系數(shù)可以設(shè)為小值。此處對(duì)中彈簧的作用是當(dāng)系統(tǒng)處于振動(dòng)平衡狀態(tài)時(shí),將活塞恢復(fù)至原位。只需克服活塞運(yùn)動(dòng)的靜摩擦力即可。

      本文中設(shè)定振源質(zhì)量為50 kg,選定活塞桿的面積A1=50 mm2,初始油壓調(diào)至1 MPa。當(dāng)振源質(zhì)量有改變時(shí),可以通過(guò)調(diào)節(jié)油壓來(lái)支撐不同的振源質(zhì)量。SMA彈簧工作在超彈性區(qū)間,調(diào)整其輸出力為50 N,選定薄膜面積A2=300 mm2。系統(tǒng)結(jié)構(gòu)參數(shù)如表1所示。

      表1 復(fù)合式隔振器的初始參數(shù)Tab.1 Initial parameter of hybrid isolator

      選定系統(tǒng)的慣容值為100 kg。

      將慣容器與擬零剛度結(jié)合在一起,從增加系統(tǒng)視在質(zhì)量與降低系統(tǒng)等效剛度兩個(gè)方面降低系統(tǒng)的固有頻率。 隔振系統(tǒng)的性能,可以通過(guò)傳遞到基礎(chǔ)的力和輸入力之比與輸入頻率之間的關(guān)系,即求得系統(tǒng)的傳遞函數(shù)。傳遞到基礎(chǔ)的力為

      (14)

      式中,F(xiàn)(s)和FT(s)分別為輸入力和傳遞力的拉式變換,通過(guò)聯(lián)立式(8),式(14),可以得到系統(tǒng)的傳遞函數(shù)

      (15)

      從幅頻曲線可以看出,系統(tǒng)的諧振頻率已經(jīng)降到0.04 Hz,慣容器和擬零剛度原理的結(jié)合增加了系統(tǒng)的隔振頻寬,可以實(shí)現(xiàn)真正的低頻隔振,提高隔振性能,如圖3所示。

      圖3 系統(tǒng)傳遞函數(shù)的幅頻曲線Fig.3 Amplitude-frequency curve of system transfer function

      4 系統(tǒng)結(jié)構(gòu)參數(shù)對(duì)響應(yīng)頻率的影響

      4.1k1,b值

      從系統(tǒng)的動(dòng)力學(xué)方程可以看出,對(duì)中彈簧的剛度k1與系統(tǒng)的等效動(dòng)態(tài)剛度有關(guān),k1越小,系統(tǒng)的等效剛度就越小,進(jìn)而就會(huì)降低系統(tǒng)的固有頻率。由于整個(gè)系統(tǒng)處于漂浮狀態(tài),對(duì)中彈簧的作用就是要讓活塞回到平衡位置,因此在保證彈簧的恢復(fù)力大于活塞的摩擦力的情況下,k1值越小,系統(tǒng)的諧振頻率就越低,隔振效果就越好。

      同理,增加系統(tǒng)的慣容值可以增加系統(tǒng)的視在質(zhì)量,從而降低系統(tǒng)的諧振頻率增加系統(tǒng)的隔振頻寬,達(dá)到良好的低頻隔振效果。常見的慣容器的慣容值都可以達(dá)到自身質(zhì)量的幾十倍甚至數(shù)百倍,但無(wú)限制的提高慣容值必定會(huì)帶來(lái)對(duì)慣容器外加負(fù)載的要求,因此選擇適當(dāng)?shù)膽T容器就已經(jīng)可以達(dá)到超低頻隔振的效果,無(wú)需單一的追求極大的慣容值。

      4.2不同A1/A2值的影響

      由式(13)的第二項(xiàng)可以看出設(shè)計(jì)參數(shù)中A1/A2的值對(duì)系統(tǒng)的等效剛度有著很大的影響,由于A1的選定與隔振物體重量和油壓的壓力范圍有關(guān)系,下面選定不同的A2值,作出系統(tǒng)的伯德圖,分析A1/A2的值對(duì)系統(tǒng)諧振頻率的影響。

      圖4 不同A1/A2下系統(tǒng)的伯德圖Fig.4 Bode diagrams of different ratios of A1/A2

      由圖4可以看出,增加活動(dòng)薄膜的面積A2,即使得A1/A2的比值越小,系統(tǒng)的等效動(dòng)態(tài)剛度就越小,從而可以獲得更寬的隔振頻寬,提高隔振性能。這一點(diǎn)很容易從物理上解釋,因?yàn)檫@一比值越小,說(shuō)明由于系統(tǒng)振動(dòng)引起的液壓緩沖裝置中的體積變化和壓力變化越小,從而使得等效的動(dòng)態(tài)剛度越小。

      5 結(jié) 論

      本文提出了一種將慣容器與擬零剛度原理相結(jié)合起的復(fù)合式低頻隔振器,可以隔離低頻甚至超低頻的振動(dòng)。通過(guò)分析,得到了系統(tǒng)的動(dòng)力學(xué)方程,分析了系統(tǒng)的固有頻率的影響因素,結(jié)果表明慣容器和擬零剛度原理的結(jié)合可以有效降低系統(tǒng)的諧振頻率,從而獲得超低頻的隔振性能。

      [1] RIVIN E I, RIVIN E I. Passive vibration isolation[M]. New York: Asme Press, 2003.

      [2] 張利國(guó), 張嘉鐘, 賈力萍, 等. 空氣彈簧的現(xiàn)狀及其發(fā)展[J]. 振動(dòng)與沖擊, 2007, 26(2): 146-151.

      ZHANG Liguo, ZHANG Jiazhong, JIA Liping. Current state and development of air spring[J]. Journal of Vibration and Shock, 2007, 26(2): 146-151.

      [3] 歐進(jìn)萍. 結(jié)構(gòu)振動(dòng)控制: 主動(dòng), 半主動(dòng)和智能控制[M]. 北京:科學(xué)出版社, 2003.

      [4] 杜華軍, 黃文虎, 鄒振祝. 航天支架結(jié)構(gòu)的被動(dòng)振動(dòng)控制[J]. 應(yīng)用力學(xué)學(xué)報(bào), 2002, 19(3): 10-13.

      DU Huajun, HUANG Wenhu, ZOU Zhenzhu. Passive vibration control of aerospace supporter[J]. Chinese Journal of Applied Mechanics, 2002, 19(3): 10-13.

      [5] CARRELLA A, BRENNAN M J, WATERS T P. Static analysis of a passive vibration isolator with quasi-zero-stiffness characteristic[J]. Journal of Sound and Vibration, 2007, 301(3): 678-689.

      [6] LE T D, AHN K K. A vibration isolation system in low frequency excitation region using negative stiffness structure for vehicle seat[J]. Journal of Sound and Vibration, 2011, 330(26): 6311-6335.

      [7] 劉興天, 孫靖雅, 肖鋒, 等. 準(zhǔn)零剛度微振動(dòng)隔振器的原理和性能研究[J]. 振動(dòng)與沖擊, 2013, 32(21): 69-73.

      LIU Xingtian,SUN Jingya,XIAO Feng. Principle and performance of a quasi-zero stiffness isolator for micro-vibration isolation[J]. Journal of Vibration and Shock, 2013, 32(21): 69-73.

      [8] MIZUNO T, TAKASAKI M, KISHITA D, et al. Vibration isolation system combining zero-power magnetic suspension with springs[J]. Control Engineering Practice, 2007, 15(2): 187-196.

      [9] CARRELLA A, BRENNAN M J, WATERS T P, et al. On the design of a high-static-low-dynamic stiffness isolator using linear mechanical springs and magnets[J]. Journal of Sound and Vibration, 2008, 315(3): 712-720.

      [10] SMITH M C. Synthesis of mechanical networks: the inerter[J]. Automatic Control, IEEE Transactions on, 2002, 47(10): 1648-1662.

      [11] 張孝良, 陳龍, 聶佳梅, 等. 2 級(jí)串聯(lián)型 ISD 懸架頻響特性分析與試驗(yàn)[J]. 江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 33(3): 255-258.

      ZHANG Xiaoliang,CHEN Long,NIE Jiamei. Analysis and experiment of frequency response characteristics of two-stage series-connected ISD suspension[J]. Journal of Jiangsu University(Natural Science), 2012, 33(3): 255-258.

      [12] 陳龍, 張孝良, 聶佳梅, 等. 基于半車模型的兩級(jí)串聯(lián)型 ISD 懸架性能分析[J]. 機(jī)械工程學(xué)報(bào), 2012, 48(6): 102-108.

      CHEN Long,ZHANG Xiaoliang,NIE Jiamei. Performance analysis of two-stage series-connected inerter-spring-damper suspension based on half-car model[J]. Journal of Mechanical Engineering, 2012, 48(6): 102-108.

      [13] 陳龍, 任皓, 汪若塵, 等. 液力式慣容器力學(xué)性能仿真與試驗(yàn)研究[J]. 振動(dòng)與沖擊, 2014, 33(12): 87-92.

      CHEN Long, REN Hao, WANG Ruochen. Simulations and tests for mechanical properties of a hydraulic inerter[J]. Journal of Vibration and Shock, 2014, 33(12): 87-92.

      [14] XU T. Design and analysis of a shock absorber with a variable moment of inertia flywheel for passive vehicle suspension[D]. Ottawa: University of Ottawa, 2013.

      [15] HU Y, LI C, CHEN M Z Q. Optimal control for semi-active suspension with inerter[C]//Control Conference (CCC), 2012 31st Chinese. Hefei: IEEE, 2012: 2301-2306.

      [16] ZHANG X J, AHMADIAN M, GUO K H. On the benefits of semi-active suspensions with inerters[J]. Shock and Vibration, 2012, 19(3): 257-272.

      [17] SAADAT S, SALICHS J, NOORI M, et al. An overview of vibration and seismic applications of NiTi shape memory alloy[J]. Smart Materials and Structures, 2002, 11(2): 218.

      [18] 王惠民. 流體力學(xué)基礎(chǔ)[M]. 北京:清華大學(xué)出版社, 2005.

      [19] LAGOUDAS D C, KHAN M M, MAYES J J, et al. Modelling of shape memory alloy springs for passive vibration isolation[C]//International Mechanical Engineering Congress and Exposition. New York: ASME, 2001.

      Designandperformanceanalysisofacompositelow-frequencyvibrationisolator

      WU Ting1, WANG Linxiang2

      (1. The State Key Lab for Fluid Power Transmission and Control, Zhejiang University, Hangzhou 310027,China)

      In current papers, a hybrid vibration isolator was proposed by combining the quasi-zero-stiffness isolator with an inerter. The inerter is employed to increase the effective mass of the isolator while the quasi-zero-stiffness part is to decrease the dynamic stiffness of the restoring elements. A superelastic shape memory alloy spring was integrated into the device. The low-frequency vibration isolation mechanism was specified. The force transmission ratio was analyzed and simulated. It is shown by simulation that the proposed vibration isolator could have a resonant frequency lower than 1 Hz. The influences of parameters on the system responses were discussed.

      quasi-zero-stiffness; inerter; superelasticity; low-frequency vibration

      TH212;TH213.3

      A

      10.13465/j.cnki.jvs.2017.19.035

      國(guó)家自然科學(xué)基金資助項(xiàng)目(51575478;61571007)

      2016-04-26 修改稿收到日期:2016-08-08

      吳庭 男,碩士生,1992年10月生

      王林翔 男,博士,教授,1971年1月生

      猜你喜歡
      振源固有頻率諧振
      Effects of Prandtl number in two-dimensional turbulent convection?
      現(xiàn)場(chǎng)測(cè)定大型水輪發(fā)電機(jī)組軸系的固有頻率
      基于諧振開關(guān)技術(shù)的低相噪LC VCO的設(shè)計(jì)
      考慮振源相互作用的船舶甲板減振方法
      船海工程(2019年1期)2019-03-04 01:58:04
      一種小型化低噪聲微波本振源設(shè)計(jì)
      電子制作(2017年9期)2017-04-17 03:00:36
      諧振式單開關(guān)多路輸出Boost LED驅(qū)動(dòng)電源
      基于CM6901 的LLC半橋諧振開關(guān)電源設(shè)計(jì)
      總溫總壓測(cè)頭模態(tài)振型變化規(guī)律研究
      A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
      轉(zhuǎn)向系統(tǒng)固有頻率設(shè)計(jì)研究
      汉中市| 丹巴县| 武义县| 驻马店市| 康保县| 雷波县| 南和县| 嘉荫县| 宁国市| 琼海市| 江津市| 乐都县| 新密市| 会同县| 丹巴县| 禄丰县| 聂荣县| 浦北县| 公安县| 兰溪市| 虹口区| 新巴尔虎右旗| 丹东市| 保亭| 略阳县| 汕头市| 东台市| 潞城市| 棋牌| 普宁市| 营山县| 甘南县| 盐池县| 监利县| 肃宁县| 铜川市| 东阳市| 即墨市| 扶余县| 大关县| 达拉特旗|