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

?

溫度和車速對半剛性基層瀝青路面模量場和疲勞壽命的影響

2021-07-01 20:54黃立葵宋偉姜正暉楊樹
湖南大學學報·自然科學版 2021年5期
關鍵詞:疲勞壽命道路工程瀝青路面

黃立葵 宋偉 姜正暉 楊樹

摘 ? 要:為分析半剛性基層瀝青路面中瀝青層模量在時間域和空間域上的分布特征,研究模量場年變化對結構層疲勞壽命的影響,將瀝青面層細分為1 cm厚的亞層,基于Odemark法經(jīng)迭代計算得到了瀝青層中的模量場分布,分析全年456種工況下路面結構的力學響應,基于Miner線性疲勞損傷模型計算得到了考慮模量場影響的疲勞壽命修正系數(shù). 結果表明:模量場隨時間顯著變化,車速為70 km/h時,瀝青層各深度處的模量年最大與最小值之比達20.1~34.9,且同一時間沿路面深度方向的模量不相等;車速對模量場有顯著影響,車速120 km/h與10 km/h時,瀝青層模量之比達1.12~2.28;路面結構的疲勞壽命分析中,應考慮模量場年變化的影響,在20 ℃、10 Hz標準工況下分析得到的上基層疲勞壽命偏保守,下基層疲勞壽命偏危險,建議采用考慮模量場影響的疲勞壽命修正系數(shù)調整標準工況下計算得到的結構層疲勞壽命.

關鍵詞:道路工程;瀝青路面;模量場;疲勞壽命;修正系數(shù)

中圖分類號:U416 ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼:A

Abstract:In order to analyze the distribution characteristics of the modulus field of semi-rigid asphalt pavement in the time domain and spatial domain and to study the effect of the annual change of the modulus field on the fatigue life of the structural layers, the asphalt layers were subdivided into several sub-layers with a thickness of 1 cm. The modulus field distribution in the asphalt layers was obtained by iterative calculations based on the Odemark method, and the fatigue life correction factors of structural layers were obtained based on the Miner linear fatigue damage model as well as the mechanical responses of the pavement structure under 456 conditions throughout the year. The results show that the modulus field varies significantly with time. At the speed of 70 km/h, the ratio of the annual maximum and minimum values of the modulus at each depth of the asphalt layer reaches 20.1~34.9, and the modulus with the pavement depths are not equal at the same time. The ratio of the modulus at the speed of 120km/h to that at 10 km/h is in the range of 1.12 ~ 2.28, which means that the vehicle speed has a significant effect on the modulus field. The influence of the annual change of the modulus field on the fatigue life of pavement structures should be considered. The fatigue life of the upper base calculated under the standard conditions at 20 ℃ and 10 Hz is conservative, while the fatigue life of the lower base is dangerous. The calculation results under the standard condition should be revised by the fatigue life correction factor when checking the fatigue life of the pavement structures.

Key words:road engineering;asphalt pavement;modulus field;fatigue life;correction factor

為獲得交通荷載作用下路面結構的真實力學響應,必須合理地選用路面材料的計算參數(shù). 作為一種典型的黏彈性材料,瀝青混合料的模量顯著依賴于溫度和荷載作用頻率. 受外界環(huán)境的影響,路面結構內部的溫度場變化顯著,而路面溫度場變化必將導致瀝青材料層沿路面深度方向的模量分布隨之變化. 此外,不同的行車速度對應不同的荷載作用時間,即使在相同的行車速度下,路面不同深度處的應力脈沖持續(xù)時間也存在較大差異,這等效于不同深度處的加載頻率不同. 因此,在確定路面材料計算參數(shù)時,必須考慮瀝青層模量的溫度和車速依賴性.

在有關路面力學響應的研究和我國現(xiàn)行路面設計規(guī)范中[1-9],大多將路面同一結構層的材料參數(shù)視為一致,這與真實情況相差很大,基于此計算得到的路面力學響應自然也不夠準確. 對于瀝青路面,荷載、環(huán)境、路面結構組成一個不可分離的系統(tǒng),在荷載、環(huán)境的綜合作用下,瀝青材料層模量以“場”的形式分布. 為了更合理地分析路面結構的力學響應,本文首先對長沙地區(qū)瀝青路面的模量場在時間域和空間域的分布特征展開研究,并分析了車速對模量場的影響. 在此基礎上,以模量場計算結果作為路面力學響應計算的材料輸入?yún)?shù),分析了模量場變化下半剛性基層瀝青路面的疲勞壽命,并以標準工況下路面結構的疲勞壽命為基準,計算得到了各工況下考慮模量場影響的疲勞壽命修正系數(shù).

1 ? 模量場計算方法

1.1 ? 路面溫度場

路面溫度場受多種因素共同影響,其中氣溫和太陽輻射是最主要的影響因素,其他因素的影響可以通過氣溫和太陽輻射體現(xiàn)出來[10-11]. 孫立軍等[11]在我國多個地區(qū)進行了大量實測,建立了式(1)所示的考慮地區(qū)差異的路面溫度場預估模型,可采用此模型計算路面不同時刻的溫度場.

1.2 ? 路面應力脈沖持續(xù)時間

Barksdale[13]、Brown[14]、Hu等[15]均提出了行車荷載下,路面結構中應力脈沖持續(xù)時間的計算方法,Odemark最早提出了等效厚度的概念,在美國最新版MEPDG中得到了應用[16],本文采用基于等效厚度的Odemark法計算路面結構中的應力脈沖持續(xù)時間.

通過式(5)將路面各結構層換算為與路基模量相同的等效層,并在換算后的等效路面結構中按45°的擴散角向下傳遞交通荷載,據(jù)此可經(jīng)迭代計算得到路面結構中任一深度處的應力影響長度Leff,應力影響長度的概念見圖1,該深度處的荷載作用時間由式(6)計算.

式中:he為待換算層的換算厚度,m;h1為換算層的實際厚度,m;E1為換算層的模量,MPa;μ1為換算層的泊松比;E2為路基的模量,MPa;μ2為路基的泊松比;t為應力脈沖持續(xù)時間,s;Leff為計算深度處的應力影響長度,m;v為行車速度,m/s.

1.3 ? 路面模量場

路面結構不同深度處的溫度和荷載作用頻率不同,受其影響,瀝青層不同深度處的模量也并不一致,為客觀反映瀝青層的模量沿路面深度方向的分布情況,應將其模量視為“場”,即“模量場”.

表2所示為我國常用的半剛性基層瀝青路面結構,選其作為研究對象具有典型性. 在計算路面結構力學響應時,瀝青層的材料參數(shù)由模量場計算結果賦值,其他結構層的材料參數(shù)按表中數(shù)據(jù)選取.

在瀝青層中,應力脈沖持續(xù)時間(荷載作用頻率)和模量兩者之間相互影響,且任一深度處的模量會受到該點上方各結構層模量的影響,因此,在求解模量場時需多次迭代計算直至結果收斂.

路面結構中模量場的具體計算步驟如下:

1)將瀝青面層細分為1 cm厚的亞層,取各亞層的中點為該層的計算點,由式(1)~(4)計算各亞層的溫度,并假設路面各深度處的荷載作用頻率均為10 Hz.

2)將路面各亞層的溫度和加載頻率數(shù)據(jù)代入式(8)~(10)中計算路面各亞層的模量.

3)基于步驟2)中計算得到的各亞層模量,由式(5)和式(6)計算各亞層的應力脈沖持續(xù)時間,并根據(jù)計算結果由式(7)計算各亞層的荷載作用頻率.

4)基于步驟3)中計算得到的荷載作用頻率數(shù)據(jù),按步驟2)重新計算路面各亞層的模量.

重復步驟2)~4),直至各亞層的模量在前后兩次計算結果中的差值小于1 MPa. 鑒于迭代計算過程繁瑣,編寫了相關程序,借助計算機自動完成迭代求解過程.

2 ? 模量場計算結果

如前所述,將瀝青面層細分為1 cm厚的若干亞層,取每個亞層的中點為該層的計算點,基于一整年的氣象數(shù)據(jù),以1 h為時間間隔,并考慮車速的影響,通過迭代計算得到了全年各工況下瀝青路面結構層逐小時模量場分布. 我國高速公路上行駛的重載貨車的車速大多為60~80 km/h,不失一般性,當分析模量場在時間域和空間域上的分布特征時,取基準車速為70 km/h.

2.1 ? 模量場時間域分布特征

瀝青結構層的模量場與外界環(huán)境緊密相關,在時間域上的波動性很大. 以E(z,v)表示車速v(km/h)時路面深度Z處的瀝青層模量. 圖2所示為路面不同深度處E (Z,70)的全年逐小時變化情況,1.5 cm、6.5 cm、13.5 cm、17.5 cm深度處的模量年最大值分別為20 224 MPa、23 246 MPa、23 865 MPa、24 982 MPa,年最小值分別為607 MPa、1 151 MPa、698 MPa、778 MPa,最大最小值之比分別達到33.3、20.2、34.2、32.1. 進一步分析了瀝青層其他深度處的模量年最大最小值之比,數(shù)據(jù)分析結果表明在70 km/h車速時,瀝青層各深度處的模量年最大最小值之比均位于20.1~34.9. 由此可見瀝青層各深度處的模量在全年的變化范圍均很大.

根據(jù)我國現(xiàn)行《公路瀝青路面設計規(guī)范》(JTG D50—2017),在分析路面結構的力學響應時,將20 ℃、10 Hz時的瀝青混合料模量作為瀝青面層的輸入?yún)?shù),由式(8)~(10)計算可得20 ℃、10 Hz時AC-13、AC-20、AC-25的模量分別為8 683 MPa、11 090 MPa、11 950 MPa,但這種材料參數(shù)取值方法并不能反映瀝青混合料的模量隨時間變化的情況. 為直觀反映瀝青層各深度處的模量隨時間的波動性,以1 000 MPa為模量區(qū)間,統(tǒng)計分析了全年8 760 h中路面各深度處的模量在各模量區(qū)間出現(xiàn)的小時數(shù),并將其與全年8760 h的比值匯總于表3中. 分析表3中數(shù)據(jù)可知,在路面0~4 cm深度范圍內(即AC-13結構層中),1.5 cm深度處的模量在全年8 760 h中處于8 000~9 000 MPa范圍內的小時數(shù)最多,但也只有456 h,僅占全年時間的5.21%;路面4~10 cm深度范圍內(即AC-20結構層中),7.5 cm深度處的模量在全年8 760 h中處于11 000~12 000 MPa范圍內的小時數(shù)最多,但也只有363 h,僅占全年時間的4.14%;路面10~18 cm深度范圍內(即AC-25結構層中),17.5 cm深度處的模量在全年8 760 h中處于11 000~12 000 MPa范圍內的小時數(shù)最多,但也只有397 h,僅占全年時間的4.53%. 由此可見,現(xiàn)行《公路瀝青路面設計規(guī)范》(JTG D50—2017)中有關瀝青混合料材料參數(shù)的取值方法并不能夠反映材料在實際工程中的情況. 因此,為保證計算結果的準確性,在進行路面結構驗算時,必須考慮瀝青層模量在時間域上的波動性對驗算結果的影響.

2.2 ? 模量場空間域分布特征

分別選取1月、4月、7月和10月作為冬、春、夏、秋季的代表月,并在每個代表月中選取日平均氣溫與月平均氣溫最為接近的一天作為該月的代表日. 圖3所示為車速70 km/h時1月28日、4月1日、7月19日和10月22日代表性時刻的路面模量場分布. 分析圖3可知,在各結構層的分層處,由于材料不同,模量具有顯著差異. 即使在同一結構層中,瀝青混合料的模量也會沿路面深度發(fā)生變化,如1月28日0:00,AC-13、AC-20、AC-25結構層中模量最大最小值之差分別為844 MPa、1 565 MPa、831 MPa,這自然也會影響路面結構的力學響應分析結果,因此在驗算路面結構的疲勞壽命時還應考慮瀝青層模量沿路面深度方向的分布情況,不應簡單地將同一結構層不同深度處的模量視為定值.

2.3 ? 瀝青層模量的速度依賴性

實際道路上行駛的車輛車速變化范圍很大,不同的行車速度對應路面結構不同的荷載作用時間,等效于荷載的加載頻率不同,進而對瀝青層各深度處的模量值和整個瀝青層中的模量場分布產(chǎn)生影響.

圖4、圖5分別為10月22日10:00瀝青層不同深度處的模量隨車速的變化情況. 數(shù)據(jù)分析表明:隨著車速的提高,各深度處的模量值均顯著增大,如當車速從10 km/h提高至120 km/h時,0.5cm深度處E (0.5,120)與E (0.5,10)之比達1.89,17.5 cm深度處E (17.5,120)與E (17.5,10)之比達1.72. 進一步分析發(fā)現(xiàn)瀝青層各深度處的模量在全年中任意時刻E (Z,120)與E (Z,10)之比均位于1.12~2.28,表明車速對瀝青路面的模量場具有顯著影響. 但在現(xiàn)行的 《公路瀝青路面設計規(guī)范》(JTG D50—2017)中,材料參數(shù)取為定值,這并不能反映不同車速時瀝青層模量的差異性,因此在路面結構驗算中還應考慮車速的影響.

3 ? 模量場在路面結構疲勞壽命驗算中的應用

不同的模量場分布特征下,路面結構的疲勞壽命具有顯著差異,因此在路面結構的疲勞壽命驗算中必須計入模量場的影響. 根據(jù)我國現(xiàn)行瀝青路面設計規(guī)范,驗算路面結構的疲勞壽命時,將20 ℃、10 Hz時的模量作為瀝青面層的輸入?yún)?shù),由前述分析可知,這不能反映真實情況,以此為基礎的計算結果自然也與實際狀況相差很大,因此很有必要在模量場計算結果的基礎上,對瀝青路面的疲勞壽命計算結果進行修正.

受外界環(huán)境因素的影響,一年中瀝青路面的模量場變化很大,需綜合分析全年中不同工況下,模量場對路面結構疲勞壽命的影響. 為此,在全年氣溫變化范圍內以1 ℃為步長,將氣溫劃分為若干個區(qū)段,在每個溫度區(qū)間內選擇一個代表時刻,將該時刻各行車速度下的模量場作為結構驗算時瀝青面層的材料輸入?yún)?shù),并基于路面結構層的力學響應計算結果和對應的疲勞方程計算對應工況下的路面結構疲勞壽命.

計算各工況下瀝青路面結構的力學響應時,使用ABAQUS軟件建立了橫向、縱向、豎向尺寸分別為8 m、8 m、6 m的路面結構有限元模型,采用C3D8單元將模型劃分為24萬個單元,模型底部和四周限制法向位移,假定各結構層間完全連續(xù). 為考慮模量場的影響,在模型中將瀝青面層細分為1 cm厚的亞層,每一亞層的材料參數(shù)根據(jù)模量場計算結果取值.

一天中不同時間段內交通量不同,因此在統(tǒng)計每個溫度區(qū)間的交通量時,需對每個時間段分別統(tǒng)計,表4所列為湖南省內某高速公路1 d中各時段內的交通量分布情況[18],在表4的基礎上按以下步驟統(tǒng)計各溫度區(qū)間的全年交通量:

1)將全年每一天的24 h按1 h的時間跨度劃分為24個區(qū)間,并假設每個時間段的交通量不隨日期變化,因此可將全年中每個時間區(qū)間的總計365 h劃入同一個統(tǒng)計時段.

2)統(tǒng)計全年8 760 h的逐小時氣溫數(shù)據(jù),統(tǒng)計結果表明長沙全年氣溫變化范圍為-7~+39 ℃,按1 ℃為溫度跨度將氣溫劃分為若干個溫度區(qū)間.

3)不同日期的同一時刻氣溫不同,基于1)中的假設,在每個統(tǒng)計時段的365 h中,處于各溫度區(qū)間的交通量占該統(tǒng)計時段365 h的總交通量的比例與所處該溫度區(qū)間的小時數(shù)占365 h的比例相同,據(jù)此可計算某統(tǒng)計時段位于某一溫度區(qū)間的交通量.

4)將不同時間段內同一溫度區(qū)間的交通量相加,即可得到全年中該溫度區(qū)間的累計交通量.

分析表6和表7中數(shù)據(jù)可知,受模量場年變化的影響,按標準工況計算得到的上基層疲勞壽命偏保守,而下基層疲勞壽命偏危險,且車速越低,對下基層疲勞壽命越不利. 因此,在路面結構驗算時應計入模量場年變化對驗算結果的影響. 為便于工程應用,可按本文方法計算我國常見氣候條件下常用路面結構的疲勞壽命修正系數(shù),路面設計時可仍取20 ℃、10 Hz時的模量作為瀝青面層的輸入?yún)?shù),計算出標準工況下基準路面結構的疲勞壽命后,乘以相應的考慮模量場影響的疲勞壽命修正系數(shù),即可得到更符合實際情況的路面結構層疲勞壽命.

4 ? 結 ? 論

本文在研究典型半剛性基層瀝青路面模量場的基礎上,基于Miner線性疲勞損傷模型和全年456種工況下的路面結構力學響應計算結果,分析了模量場年變化對半剛性基層疲勞壽命的影響,得到如下結論:

1)時間域上,瀝青層各深度處的模量均具有很大的波動性,車速為70 km/h時,瀝青層各深度處模量的年最大最小值之比達20.1~34.9,且同一時刻同一結構層中不同深度處的模量不相等.

2)車速變化時,瀝青層模量隨之顯著變化,路面各深度處的E (Z,120)與E (Z,10)之比達1.12~2.28.

3)從氣溫變化維度考慮模量場對半剛性基層疲勞壽命的影響,基于Miner線性疲勞損傷模型計算考慮模量場影響的路面結構全年損傷,定義疲勞壽命修正系數(shù)為20 ℃、10 Hz標準工況下基準路面結構的全年損傷與考慮模量場變化的路面結構全年損傷之比.

4)由于模量場年變化的影響,按標準工況計算得到的上基層疲勞壽命偏保守;而下基層疲勞壽命偏危險,且車速越低,對下基層疲勞壽命越不利.

5)為計入模量場年變化對瀝青路面疲勞壽命驗算結果的影響,可按照本文方法計算我國常見氣候條件下常用路面結構的疲勞壽命修正系數(shù),路面設計時可仍按標準工況計算路面結構的疲勞壽命,再乘以相應的疲勞壽命修正系數(shù)即可.

參考文獻

[1] ? ?YOU L Y,YAN K Z,HU Y B,et al. Spectral element method for dynamic response of transversely isotropic asphalt pavement under impact load[J]. Road Materials and Pavement Design,2018,19(1):223—238.

[2] ? ?董澤蛟,潘小康,邵顯智,等. 溫度場及非均布移動荷載作用下瀝青路面力學響應分析[J]. 同濟大學學報(自然科學版),2016,44(5):740—746.

DONG Z J,PAN X K,SHAO X Z,et al. Dynamic response analysis of asphalt pavement due to combined non-uniform moving load and temperature distribution[J]. Journal of Tongji University(Natural Science),2016,44(5):740—746. (In Chinese)

[3] ? ?劉能源,顏可珍,胡迎斌,等. 基于橫觀各向同性的瀝青路面加鋪層力學分析[J]. 湖南大學學報(自然科學版),2017,44(5):96—103.

LIU N Y,YAN K Z,HU Y B,et al. Mechanical analysis of asphalt concrete overlay over old asphalt pavements considering cross-anisotropy properties[J]. Journal of Hunan University(Natural Sciences),2017,44(5):96—103. (In Chinese)

[4] ? ?萬晨光,申愛琴,薛翠真,等. 層間接觸狀態(tài)對橋面鋪裝結構力學響應的影響[J]. 湖南大學學報(自然科學版),2016,43(9):113—119.

WAN C G,SHEN A Q,XUE C Z,et al. Effcet of interlayer contact state on mechanical response of bridge deck pavement structure[J]. Journal of Hunan University (Natural Sciences),2016,43(9):113—119. (In Chinese)

[5] ? ?郭寅川,王禮根,申愛琴,等. 動載作用下級配碎石基層瀝青路面力學分析[J]. 公路交通科技,2016,33(1):27—33.

GUO Y C,WANG L G,SHEN A Q,et al. Dynamic analysis of graded crushed stone as base of asphalt pavement[J]. Journal of Highway and Transportation Research and Development,2016,33(1):27—33. (In Chinese)

[6] ? ?SUN Y Z,GUO R,GAO L,et al. Study on dynamic response characteristics of saturated asphalt pavement under multi-field coupling[J]. Materials,2019,12(6):959.1—959.15.

[7] ? ?劉大維,戴宗宏,陳洋,等. 車輛多輪動載作用下柔性路面動應力響應[J]. 中國公路學報,2017,30(11):36—44.

LIU D W,DAI Z H,CHEN Y,et al. Dynamic stress response of flexible pavement under multiple wheel dynamic loads of vehicle[J]. China Journal of Highway and Transport,2017,30(11):36—44. (In Chinese)

[8] ? ?XIA R X,LI J H,HE J,et al. Effect analysis of vehicle system parameters on dynamic response of pavement[J]. Mathematical Problems in Engineering,2015,2015:1—8.

[9] ? ?公路瀝青路面設計規(guī)范:JTG D50—2017[S]. 北京:人民交通出版社,2017:12—20.

Specifications for design of highway asphalt pavement:JTG D50—2017[S]. Beijing:China Communications Press,2017:12—20. (In Chinese)

[10] ?宋小金,樊亮. 瀝青路面結構溫度隨深度變化規(guī)律研究[J]. 土木工程學報,2017,50(9):110—117.

SONG X J,F(xiàn)AN L. Study on the variation rules of temperature with depth for asphalt pavement structure[J]. China Civil Engineering Journal,2017,50(9):110—117. (In Chinese)

[11] ?孫立軍,秦健. 瀝青路面溫度場的預估模型[J]. 同濟大學學報(自然科學版),2006,34(4):480—483.

SUN L J,QIN J. Prediction model on temperature field in asphalt pavement[J]. Journal of Tongji University(Natural Science),2006,34(4):480—483. (In Chinese)

[12] ?曹其夢,于瑛,楊柳. 太陽逐時總輻射計算模型適用性分析:以我國部分地區(qū)為例[J]. 太陽能學報,2018,39(4):917—924.

CAO Q M,YU Y,YANG L. Applicability analysis of hourly total solar radiation calculation model-taking some regions of China as example[J]. Acta Energiae Solaris Sinica,2018,39(4):917—924. (In Chinese)

[13] ?BARKSDALE R D. Compressive stress pulse times in flexible pavements for use in dynamic testing[J]. Highway Research Record,1971,345(4):32—44.

[14] ?BROWN S F. Determination of youngs modulus for bituminous materials in pavement design[J]. Highway Research Record,1973,431:38—49.

[15] ?HU X D,ZHOU F J,HU S,et al. Proposed loading waveforms and loading time equations for mechanistic-empirical pavement design and analysis[J]. Journal of Transportation Engineering,2010,136(6):518—527.

[16] ?Applied Research Associates. Guide for mechanistic-empirical pavement design:NCHRP project 1-37A[R]. Washington DC:Transportation Research Board,2004:2—7.

[17] ?宋小金,曾夢瀾,樊亮. 瀝青路面反算模量與瀝青混合料動態(tài)模量的關系[J]. 土木工程學報,2019,52(4):121—128.

SONG X J,ZENG M L,F(xiàn)AN L. Relationship between back-calculated modulus of asphalt pavement and dynamic modulus of asphalt mixture[J]. China Civil Engineering Journal,2019,52(4):121—128. (In Chinese)

[18] ?李喜,王選倉,房娜仁,等. 基于溫度與荷載實際耦合的瀝青路面車轍預估[J]. 長安大學學報(自然科學版),2018,38(5):67—75.

LI X,WANG X C,F(xiàn)ANG N R,et al. Rutting prediction of asphalt pavement based on actual coupling of temperature and load[J]. Journal of Changan University (Natural Science Edition),2018,38(5):67—75. (In Chinese)

猜你喜歡
疲勞壽命道路工程瀝青路面
霧封層技術在瀝青路面養(yǎng)護中的應用
長壽命瀝青路面結構研究
透水瀝青路面淤塞問題分析
重載交通瀝青路面荷載圖式探討
激光沖擊強化對金屬材料疲勞壽命的影響及應用
淺析道路施工現(xiàn)場的組織與管理
隨機振動載荷下發(fā)射裝置尾罩疲勞壽命分析
探討道路工程路面水穩(wěn)層的施工與質量控制
芻議在道路工程瀝青混凝土施工中質量的控制對策
泛談飛機機體機構的強度設計
乌拉特后旗| 宁乡县| 巴东县| 嵩明县| 惠来县| 邳州市| 栾城县| 永靖县| 新丰县| 屏边| 凤阳县| 吉隆县| 新河县| 宁阳县| 平塘县| 科技| 宣武区| 军事| 栾城县| 罗城| 疏勒县| 华亭县| 铜山县| 沙湾县| 平原县| 阿尔山市| 洞头县| 贵溪市| 临汾市| 梅州市| 阳高县| 荥经县| 望奎县| 渝中区| 章丘市| 且末县| 朔州市| 敖汉旗| 呼玛县| 二连浩特市| 阿鲁科尔沁旗|