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

    排水通道分布式布設(shè)下雙層地基平面應(yīng)變固結(jié)分析

    2020-01-17 01:38:32陳益峰胡少偉梅國雄
    工程力學(xué) 2020年1期
    關(guān)鍵詞:邊界條件雙層分布式

    陳 征,張 峰,陳益峰,胡少偉,梅國雄

    (1.武漢大學(xué)水資源與水電工程科學(xué)國家重點實驗室,湖北,武漢 430072;2.武漢大學(xué)水工巖石力學(xué)教育部重點實驗室,湖北,武漢 430072;3.中建海峽建設(shè)發(fā)展有限公司科技與設(shè)計管理部, 福建,福州 350015;4.南京水利科學(xué)研究院材料結(jié)構(gòu)研究所,江蘇,南京 210024;5.廣西大學(xué)土木建筑工程學(xué)院,廣西,南寧 530004)

    工程中,排水固結(jié)法作為處理軟土地基的有效方法之一,在我國沿海地區(qū)被廣泛使用[1]。而水平排水體作為固結(jié)排水通道,其布設(shè)方式對軟土地基土固結(jié)效率具有重要的影響[2]。因此,研究如何合理布設(shè)軟土地基排水通道具有十分重要的工程應(yīng)用價值。

    目前,針對地基表面排水通道的研究主要分為排水通道的數(shù)學(xué)表征及優(yōu)化兩方面。其中,排水通道在固結(jié)理論中反映的是地基土的邊界條件,其數(shù)學(xué)表征主要分為三類:Terzaghi等[3]提出的完全透水邊界、Gray[4]提出的豎向半透水邊界、Gibson和Shefford[5]提出的水平向透水邊界。不同學(xué)者針對上述排水通道的不同數(shù)學(xué)表征方式開展了大量研究[6-13]。排水通道的優(yōu)化研究則主要體現(xiàn)在不同水平排水材料(如砂墊層、碎石墊層、復(fù)合土工布等)排水效果檢驗[14-16]和水平排水體厚度、寬度及滲透性優(yōu)化等方面[17-19]。

    總體來說,上述研究工作皆基于地基表面排水通道全面布設(shè)的前提下開展的。然而由Terzaghi固結(jié)理論推演出的典型地基無量綱出水速率-時間因數(shù)曲線關(guān)系[20](如圖1所示)可知,地基土排水速率在固結(jié)初期急劇下降(時間因數(shù)10-4增至10-1時,無量綱排水速率由56.42降至1.91)。固結(jié)中后期,地基排水速率相對較低,最終趨于不排水。因此,過量的排水通道被極大地浪費。工程中最為關(guān)心的通常不是固結(jié)前期,而是后期(如固結(jié)度達(dá)85%左右)。在達(dá)到設(shè)計固結(jié)度所額外增長的時間在工程中能被接受的前提下,為盡可能地節(jié)約資源,范鶴濱和梅國雄[21]提出了排水通道等間距條形分布式布設(shè)的新方法,變排水通道全面布設(shè)為分布式布設(shè),使得排水通道總面積大幅降低,并通過數(shù)值模擬驗證了分布式布設(shè)排水通道對工程中所關(guān)心的固結(jié)后期的固結(jié)度影響很小。隨后,Yao等[22]基于Biot固結(jié)理論對成層地基下排水通道分布式布設(shè)問題進(jìn)行了數(shù)值分析,得到了影響分布式排水通道布設(shè)方式適用性的各種影響因素。上述研究皆基于數(shù)值模擬,難以將分布式排水通道應(yīng)用于工程設(shè)計。為此,筆者在前人研究的基礎(chǔ)上,給出了排水通道條形分布式布設(shè)下單層地基平面應(yīng)變固結(jié)半解析解,并繪出了設(shè)計圖,為指導(dǎo)分布式排水通道高效且經(jīng)濟(jì)的布設(shè)提供了理論依據(jù)[20]。

    圖1 Terzaghi一維固結(jié)地基土無量綱排水速率及平均固結(jié)度隨時間變化曲線[20]Fig.1 Variations of dimensionless drainage velocity and average degree of consolidation with time factor under one dimensional Terzaghi’s consolidation[20]

    本文在已有研究工作的基礎(chǔ)上,對條形分布式排水通道下雙層地基平面應(yīng)變固結(jié)問題展開深入研究,并對不同水平排水通道幾何參數(shù)及地基土物理參數(shù)進(jìn)行固結(jié)敏感性分析以探究雙層地基條件下影響分布式排水邊界適用性的主要參數(shù),為工程設(shè)計提供理論參考。

    1 數(shù)學(xué)模型

    1.1 問題描述

    如圖2(a)所示,雙層飽和粘土上下層厚度分別為H1和H2,且H1+H2=H;頂面透水,底面不透水;kh1、kv1、mv1分別為上層土的水平滲透系數(shù)、豎向滲透系數(shù)和體積壓縮系數(shù);kh2、kv2、mv2分別為下層土的水平滲透系數(shù)、豎向滲透系數(shù)和體積壓縮系數(shù);寬度為2L的條形水平排水通道以間距2D在地基表面等間距布設(shè)。參考Barron[23]研究砂井地基固結(jié)理論的方法,本文根據(jù)條形排水通道等寬、等間距布設(shè)的對稱特性,選取圖2(a)中特征單元為研究對象,具體參數(shù)及邊界如圖2(b)所示。

    圖2 固結(jié)模型簡圖Fig.2 Sketch of the consolidation model

    1.2 控制方程

    在滿足平面應(yīng)變假設(shè)的基礎(chǔ)上,雙層地基二維自由應(yīng)變固結(jié)控制方程為:

    式中:γw為水的重度;u1和u2分別為上、下土層超靜孔壓;t為固結(jié)時間。

    1.3 求解條件

    如圖2(b)所示,地基表面分為兩部分區(qū)域:水平導(dǎo)水層布設(shè)區(qū)與未布設(shè)區(qū)。其中,未布設(shè)砂墊層區(qū)域地基邊界條件介于完全排水與隔水之間,而未布設(shè)砂墊層區(qū)域排水性越好,排水通道分布式布設(shè)對地基固結(jié)影響越小。因此,為保守起見,選取對排水通道分布式布設(shè)最不利工況進(jìn)行研究,即假定排水通道未布設(shè)區(qū)為隔水邊界條件,則地基表面邊界條件為:

    地基底部為隔水邊界,滿足:

    層間孔壓及滲流連續(xù)條件為:

    模型側(cè)向為對稱邊界條件,即:

    地基表面瞬時施加恒定荷載u0,初始條件為:

    1.4 無量綱化

    為便于方程求解及敏感性分析,定義如下無量綱參數(shù)及變量:

    應(yīng)用上述無量綱參數(shù)及變量,模型簡圖如圖2(c)所示,模型控制方程及求解條件轉(zhuǎn)換為:

    控制方程:

    豎向邊界及滲流連續(xù)條件:

    側(cè)向邊界條件:

    初始條件:

    由控制方程及求解條件可知,控制排水通道分布式布設(shè)雙層地基平面應(yīng)變固結(jié)的主要參數(shù)包括時間變量Tv1;幾何參數(shù)λ、η、δ;地基土物理參數(shù)κ1、κ2、ζ、μ。

    2 求解過程

    2.1 邊界轉(zhuǎn)換

    常規(guī)邊界條件下,固結(jié)模型的求解在空間上常采用積分變換法[24-26]。本文模型由于水平排水通道的分布式布設(shè),地基表面邊界為混合邊界條件(第一類邊界條件和第二類邊界條件同時存在于同一邊界上),而混合邊界條件下的固結(jié)問題難以直接采用積分變換法求解。為此,本文根據(jù)Chen等[20]提出的邊界轉(zhuǎn)換法,通過中間函數(shù)將地基表面混合邊界條件轉(zhuǎn)換為同類邊界條件。為便于中間函數(shù)的求解,將z=0處邊界式(10)及式(11)統(tǒng)一轉(zhuǎn)換為第二類邊界條件,即:

    式中,q(X,Tv1)為中間函數(shù)。

    2.2 超靜孔壓及平均固結(jié)度解答

    首先,根據(jù)初始條件對固結(jié)控制方程及邊界條件進(jìn)行Laplace變換;然后,根據(jù)側(cè)向邊界條件對Laplace域內(nèi)方程及豎向邊界條件進(jìn)行有限Fourier余弦變換及逆變換,從而得到無量綱超靜孔壓在Laplace域內(nèi)的解答;最后,通過局部離散法求得中間函數(shù)q在Laplace域內(nèi)節(jié)點值,并代入原解答,具體求解過程參見附錄A。結(jié)果羅列如下:

    式中:

    其余各參數(shù)參見附錄A。

    雙層地基土基于孔壓計算的平均固結(jié)度在Laplace域內(nèi)為:

    式中:

    在式(21)的矩陣計算中,Dij內(nèi)級數(shù)求和的收斂速度直接影響的求解速度。本文采用Shanks法加速級數(shù)的收斂[27]。將求得的代入式(17)、式(20)再采用數(shù)值Laplace逆變換法即可得到無量綱超靜孔隙水壓力及平均固結(jié)度[28]。

    3 解答特例及數(shù)值驗證

    3.1 解答特例

    3.1.1λ=100%

    當(dāng)λ= 1時,根據(jù)三角函數(shù)的性質(zhì),有:

    故函數(shù)gmi(X)滿足:

    使得

    將式(24)代入式(21),并通過高斯消元法可求得:

    從式(25)可以看出地基表面各位置均相等,故當(dāng)λ= 100%時,全面布設(shè)排水通道下二維固結(jié)模型等價于雙層地基一維固結(jié)模型[7]。

    3.1.2κ1η2→ ∞

    當(dāng)κ1η2→ ∞時,根據(jù)式(A11)矩陣運算可求得:

    故有

    將式(27)代入式(21)可得:

    式(28)代入式(18)與λ= 1時式(25)代入式(18)所求得的相同,故通過式(20)所求的平均固結(jié)度亦相同,即κ1η2→ ∞時,分布式排水邊界下二維固結(jié)模型可轉(zhuǎn)換為全面布設(shè)排水通道下二維固結(jié)模型,亦即是雙層地基一維固結(jié)模型[7]。

    3.2 數(shù)值驗證

    為深入驗證本文解答的正確性,采用數(shù)值有限元法與本文半解析解進(jìn)行對比驗證。為簡化數(shù)值計算,本文所建立的模型是與式(9)~式(15)相對應(yīng)的無量綱化數(shù)學(xué)模型,模型參數(shù)如表1所示。時間間隔的選取隨時間因數(shù)的增大而逐漸增大,時間在每個相鄰冪次段皆被等分為90段(例如,時間因數(shù)在100~101之間的時間間隔為10-1)。計算模型的網(wǎng)格均勻劃分,網(wǎng)格大小為0.5×0.5(如圖3所示)。圖3中標(biāo)注了A、B、C、D和E五個點用以進(jìn)行數(shù)值解與本文半解析解的對比。

    表1 數(shù)值模型計算參數(shù)Table 1 Parameters of the numerical model

    圖3 模型的網(wǎng)格及邊界條件Fig.3 Mesh and boundary conditions of the model

    圖4給出了表1工況下數(shù)值解與本文半解析解所計算的超靜孔壓對比曲線。從圖4明顯可以看出本文解與數(shù)值解幾乎完全重合,從而驗證了本文解答的正確性,也證實了本文解答所采用的局部離散法及Laplace數(shù)值逆變換法的有效性與準(zhǔn)確性。

    圖4 本文半解析解與數(shù)值解所計算的超靜孔壓對比曲線Fig.4 Comparison of excess pore water pressure calculated from the proposed semi-analytical solution and numerical solution

    4 固結(jié)敏感性分析

    本文分別從排水通道物理參數(shù)和地基土物理及力學(xué)參數(shù)兩個角度來研究排水通道分布式布設(shè)對地基固結(jié)特性的影響。進(jìn)行敏感性分析的參數(shù)包括水平排水通道參數(shù):λ、η;地基土參數(shù):上層地基土κ1和下層地基土κ2、ζ、μ。為便于分析,提出了固結(jié)時間增長率Et(U)來表征排水通道的分布式布設(shè)相比于全面布設(shè)達(dá)到特定固結(jié)度U所用時間的增長率。Et(U)可表示為:

    式中,td(U)和tf(U)分別表示水平排水通道分布式布設(shè)和全面布設(shè)達(dá)到固結(jié)度U所用時間。

    4.1 水平排水通道幾何參數(shù)

    參數(shù)λ=L/(L+D)表征條形水平排水通道布設(shè)面積與地基土表面積之比,即布設(shè)率。當(dāng)λ=100%時,水平排水通道分布式布設(shè)退化為全面布設(shè)。圖5反映了不同布設(shè)率對固結(jié)度的影響,從該圖可以看出,當(dāng)λ=60%時,固結(jié)度曲線已幾乎與全面布設(shè)下的固結(jié)度曲線重疊。相對于水平排水通道全面布設(shè),分布式布設(shè)情況下地基平均固結(jié)度與全面布設(shè)間的差值(ΔU(Tv1) =UT(Tv1) ?U(Tv1),其中UT(Tv1)為全面布設(shè)砂墊層地基平均固結(jié)度)由初始時刻0逐漸增大至極值,隨后逐漸減小至0。

    表2給出了固結(jié)后期不同布設(shè)率下時間增長率隨平均固結(jié)度的變化關(guān)系。值得一提的是,當(dāng)固結(jié)度越高,采用分布式布設(shè)水平排水通道與全面布設(shè)相比對固結(jié)效率的影響越小。例如在布設(shè)率為60%的情況下,地基土平均固結(jié)度達(dá)到85%所需固結(jié)時間相比于全面布設(shè)增加了4.54%,而達(dá)到90%所需的固結(jié)時間則僅增加了4.37%。

    圖5 不同布設(shè)率λ下地基平均固結(jié)度演化規(guī)律Fig.5 Evolution of average degree of consolidation for different pave rate λ

    表2 不同布設(shè)率λ下時間增長率隨平均固結(jié)度的變化Table 2 Increasing rate of consolidation time for a certain average degree of consolidation with different λ

    結(jié)合圖5及表2分析表明,相比排水通道全面布設(shè),分布式布設(shè)雖然不能完全滿足固結(jié)初期的地基排水需求,但能很好地滿足固結(jié)中后期的排水要求,同時對于工程中較為關(guān)注的地基固結(jié)后期的固結(jié)效率影響較小。

    參數(shù)η=H/(L+D)為特征單元的厚-寬比,表征條形水平排水通道布設(shè)條數(shù)。圖6給出了定布設(shè)率λ下厚-寬比η對地基固結(jié)的影響。隨著厚-寬比的增大,水平排水通道布設(shè)條數(shù)逐漸增多,水平滲流路徑逐漸減小,從而加快了水平向滲流速度使得平均固結(jié)度亦隨之增大。當(dāng)η趨于∞時,水平滲流路徑逐漸減小直至趨于0,徑向滲流趨于瞬時完成,土體僅產(chǎn)生一維豎向滲流,使得地基土平均固結(jié)度曲線逐漸靠近水平排水通道全面布設(shè)情況下的平均固結(jié)度曲線,這也與3.1.2節(jié)中η →∞結(jié)果一致。表3為不同厚-寬比η下不同平均固結(jié)度對應(yīng)的時間增長率。由表3可見,厚-寬比越大,固結(jié)時間增長率越低。如地基土平均固結(jié)度達(dá)到90%時,厚-寬比η= 4對應(yīng)的固結(jié)時間增長率為8.99%,而厚-寬比η= 8對應(yīng)的固結(jié)時間增長率僅為4.37%,可見厚-寬比對固結(jié)效率影響較為顯著。

    圖6 不同厚-寬比η下地基平均固結(jié)度隨時間變化曲線Fig.6 Evolution of average degree of consolidation for different thickness-width ratios η

    表3 不同厚-寬比η下時間增長率隨平均固結(jié)度的變化Table 3 Increasing rate of consolidation time for a certain average degree of consolidation with different thickness-width ratios η

    綜合表2及表3可以發(fā)現(xiàn),增加厚-寬比(排水通道的布設(shè)條數(shù))雖然可以顯著地降低時間增長率,但并不利于施工的開展。為降低過量排水通道造成的極大浪費,又對工程所關(guān)心的地基固結(jié)度影響較小,僅僅通過增加排水通道條數(shù)顯然是不可取的,需要聯(lián)合布設(shè)率進(jìn)行整體優(yōu)化設(shè)計。

    4.2 地基土物理參數(shù)

    4.2.1 上層土

    為研究上層土滲透各向異性系數(shù)κ1對固結(jié)效率的影響,圖7給出了四種不同κ1值下地基土平均固結(jié)度隨時間的演化規(guī)律??梢园l(fā)現(xiàn),排水通道分布式布設(shè)情況下,隨著κ1值的增大,平均固結(jié)度曲線逐漸向全面布設(shè)情況下的平均固結(jié)度曲線靠近。這是由于隨著κ1值的增大,上層土水平向滲透系數(shù)隨之增大,使得土體水平向滲流速度加快,從而增加了土體固結(jié)速度。當(dāng)κ1→ ∞時,土體固結(jié)過程中水平向滲流趨于瞬時完成,土體二維固結(jié)等價于一維豎向固結(jié),這與4.1.2節(jié)中特例的結(jié)論相符。

    結(jié)合表4及圖7可以看出:當(dāng)κ1> 100時,分布式排水通道的布設(shè)對固結(jié)效率的不利影響幾乎可以忽略;當(dāng)10-1<κ1<100時,分布式排水通道的布設(shè)對固結(jié)已經(jīng)產(chǎn)生了較為明顯的影響,但依然可以通過增加布設(shè)率λ及厚-寬比η來降低固結(jié)時間增長率;當(dāng)κ1<10-1時,分布式排水對固結(jié)的影響較大,即使通過增加布設(shè)率λ及厚-寬比η,也難以應(yīng)用于實際工程。

    圖7 不同滲透各向異性系數(shù)κ1對地基平均固結(jié)度的影響Fig.7 Influence of anisotropy coefficients κ1 on the average degree of consolidation

    表4 不同參數(shù)κ1下時間增長率隨平均固結(jié)度的變化Table 4 Increasing rate of consolidation time for a certain average degree of consolidation with different κ1

    4.2.2 下層土

    下層土滲透各向異性系數(shù)κ2值的變化體現(xiàn)在雙層地基下層土水平向滲透系數(shù)的變化。圖8反映了排水通道分布式布設(shè)和全面布設(shè)兩種情況下κ2值對固結(jié)過程的影響。從圖8可以發(fā)現(xiàn),4組不同κ2值下分布式排水平均固結(jié)度曲線與全面排水平均固結(jié)度曲線幾乎完全重合,即κ2值的變化對地基整體固結(jié)影響相對較小,表5中的數(shù)據(jù)也證實了這一點。這是因為排水通道在上層土表面,土體固結(jié)自上而下,下層土水平向滲透系數(shù)的變化對土體整體固結(jié)并不敏感。

    參數(shù)ζ的變化表征了地基下層土滲透性的變化。圖9及表6給出了分布式排水條件下參數(shù)ζ對地基土固結(jié)效率的影響。分析圖9及表6發(fā)現(xiàn):雙層土地基固結(jié)中,當(dāng)上層土滲透性較下層土強時(即ζ<100),排水通道分布式布設(shè)對地基土固結(jié)的不利影響較?。环粗?,不利影響則相對較大。

    圖8 不同滲透各向異性系數(shù)κ2對地基平均固結(jié)度的影響Fig.8 Influence of anisotropy coefficients κ2 on the average degree of consolidation

    表5 不同參數(shù)κ2下時間增長率隨平均固結(jié)度的變化Table 5 Increasing rate of consolidation time for a certain average degree of consolidation with different κ2

    圖9 不同參數(shù)ζ下地基平均固結(jié)度隨時間變化曲線Fig.9 Evolution of average degree of consolidation for different ζ

    參數(shù)μ的變化反映了下層土體積壓縮系數(shù)的變化。圖10及表7給出了分布式排水條件下參數(shù)μ對地基土固結(jié)效率的影響。分布式排水通道對固結(jié)不利影響隨著參數(shù)μ的增大而逐漸減小。由固結(jié)方程式(9)可知,參數(shù)ζ越大或μ越小,下層土固結(jié)系數(shù)越大,間接相當(dāng)于雙層地基上層土體的固結(jié)系數(shù)越小,而固結(jié)從靠近排水通道的上層土體開始,使得分布式排水通道對固結(jié)效率造成的不利影響亦隨之增大。

    表6 不同參數(shù)ζ下時間增長率隨平均固結(jié)度的變化Table 6 Increasing rate of consolidation time for a certain average degree of consolidation with different ζ

    圖10 不同參數(shù)μ下地度隨時變化曲Fig.10 Evolution of average degree of consolidation for different μ

    表7 不同參數(shù)μ下時間增長率隨平均固結(jié)度的變化Table 7 Increasing rate of consolidation time for a certain average degree of consolidation with different μ

    通過總結(jié)雙層地基上、下層土體物理參數(shù)對分布式排水通道下地基土固結(jié)效率的影響,可以發(fā)現(xiàn):總體來說,下層土體物理參數(shù)相較于上層土體對分布式排水通道下地基土固結(jié)效率地影響較小。在工程設(shè)計中,分布式排水通道的布設(shè)應(yīng)優(yōu)先考慮上層土體的物理力學(xué)性質(zhì),通過合理的選擇分布式排水通道的幾何參數(shù)達(dá)到與全面布設(shè)排水通道幾乎相同的固結(jié)效率,從而在不延誤工期的情況下盡可能地提高資源利用率、降低資源損耗。

    5 結(jié)論

    本文基于平面應(yīng)變假設(shè),建立了條形分布式布設(shè)排水通道下雙層地基平面應(yīng)變固結(jié)模型,并通過邊界轉(zhuǎn)換法求得了固結(jié)問題的半解析解。在通過數(shù)值法驗證了解答正確性的基礎(chǔ)上,開展了固結(jié)參數(shù)敏感性分析,具體結(jié)論羅列如下:

    (1) 相比于排水通道全面布設(shè),分布式布設(shè)率λ越高,分布式排水通道對固結(jié)的抑制作用越小,且固結(jié)時間增長率越低。固結(jié)后期,相比于全面布設(shè)排水通道,分布式布設(shè)對地基整體固結(jié)效率的降低幾乎可以忽略,尤其當(dāng)布設(shè)率λ較高時。

    (2) 隨著厚-寬比η的增大,固結(jié)時間增長率顯著降低,但并不利于施工的開展。為減少排水通道的布設(shè),又達(dá)到與全面布設(shè)排水通道幾乎相同的固結(jié)效率,需要聯(lián)合布設(shè)率λ及厚-寬比η進(jìn)行整體優(yōu)化設(shè)計。

    (3) 上層土滲透各向異性系數(shù)κ1值越大,排水通道分布式布設(shè)對地基整體固結(jié)效率的影響越低,且越利于分布式排水通道的應(yīng)用。

    (4) 地基下層土參數(shù)中,滲透各向異性系數(shù)κ2的變化對固結(jié)過程影響極其微?。欢鴿B透參數(shù)ζ越大或土體積壓縮參數(shù)μ越小,分布式排水通道對固結(jié)效率造成的不利影響則越大。

    (5) 總體上,雙層地基下層土體物理參數(shù)相較于上層土體對分布式排水通道下地基土固結(jié)效率的影響較小。在工程設(shè)計中,應(yīng)優(yōu)先考慮上層土體的物理力學(xué)性質(zhì)來合理地選擇分布式排水通道的幾何參數(shù)。

    附錄A:

    根據(jù)模型初始條件式(15)先對控制方程及邊界條件進(jìn)行Laplace變換,可得:

    由模型邊界條件(A5),對控制方程(A1)及邊界條件(A2)~(A4)進(jìn)行有限Fourier余弦變換,得:

    方程組(A6)為二階常微分方程,根據(jù)邊界條件及連續(xù)條件求得其解答為:

    式中,n=1,2;待定系數(shù)A11、A12、A21、A22可由邊界條件(A7)~(A9)求出,待定系數(shù)滿足:

    對解答(A1)進(jìn)行有限Fourier余弦逆變換,無量綱超靜孔隙水壓力可寫為:

    圖A1 水平排水通道局部離散Fig.A1 Discretization of horizontal drainage body

    猜你喜歡
    邊界條件雙層分布式
    一類帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動邊值問題的漸近解
    墨爾本Fitzroy雙層住宅
    分布式光伏熱錢洶涌
    能源(2017年10期)2017-12-20 05:54:07
    分布式光伏:爆發(fā)還是徘徊
    能源(2017年5期)2017-07-06 09:25:54
    次級通道在線辨識的雙層隔振系統(tǒng)振動主動控制
    基于DDS的分布式三維協(xié)同仿真研究
    傳統(tǒng)Halbach列和雙層Halbach列的比較
    帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
    一種雙層寬頻微帶天線的設(shè)計
    久久影院123| 国产欧美日韩精品一区二区| 国产黄片视频在线免费观看| 日本色播在线视频| 另类亚洲欧美激情| 99热这里只有是精品在线观看| 精品卡一卡二卡四卡免费| 国国产精品蜜臀av免费| 26uuu在线亚洲综合色| 亚洲久久久国产精品| 极品少妇高潮喷水抽搐| 亚洲欧美日韩另类电影网站| 亚洲国产精品专区欧美| 国产免费一级a男人的天堂| 黄色配什么色好看| 国产精品不卡视频一区二区| 女的被弄到高潮叫床怎么办| 午夜91福利影院| 大片电影免费在线观看免费| 美女脱内裤让男人舔精品视频| 国产成人freesex在线| 亚洲真实伦在线观看| 男人爽女人下面视频在线观看| 久久国内精品自在自线图片| 最后的刺客免费高清国语| 亚洲内射少妇av| 国产免费一级a男人的天堂| 久久精品国产亚洲av天美| 在现免费观看毛片| 十分钟在线观看高清视频www | 欧美一级a爱片免费观看看| 成年人午夜在线观看视频| 大片电影免费在线观看免费| 日韩人妻高清精品专区| 国产av一区二区精品久久| 亚洲不卡免费看| 毛片一级片免费看久久久久| 亚洲怡红院男人天堂| 久久久亚洲精品成人影院| 少妇被粗大的猛进出69影院 | 久久久久久久久久人人人人人人| 日日啪夜夜爽| 大片免费播放器 马上看| 日日撸夜夜添| 男女免费视频国产| 少妇精品久久久久久久| 国产精品免费大片| 久久久久国产网址| 在线亚洲精品国产二区图片欧美 | 久久6这里有精品| 女人精品久久久久毛片| 在线天堂最新版资源| 男人添女人高潮全过程视频| 高清在线视频一区二区三区| 亚洲精品国产av成人精品| 国产av国产精品国产| 成人特级av手机在线观看| 又粗又硬又长又爽又黄的视频| 97超视频在线观看视频| 在线观看免费日韩欧美大片 | 国产伦精品一区二区三区四那| 狠狠精品人妻久久久久久综合| 精品少妇内射三级| 狂野欧美激情性xxxx在线观看| 亚洲内射少妇av| 久久6这里有精品| 亚洲欧美中文字幕日韩二区| 精品国产国语对白av| 国产乱来视频区| 黑人巨大精品欧美一区二区蜜桃 | 精品久久久久久电影网| 我的女老师完整版在线观看| 亚洲国产欧美在线一区| 中国国产av一级| 中国国产av一级| 国产精品福利在线免费观看| 精品久久久噜噜| 国产日韩欧美在线精品| 一级毛片aaaaaa免费看小| 最近中文字幕2019免费版| 少妇猛男粗大的猛烈进出视频| 高清黄色对白视频在线免费看 | a级毛片免费高清观看在线播放| 永久网站在线| 永久网站在线| 亚洲第一区二区三区不卡| 男女免费视频国产| 亚洲精品亚洲一区二区| 99热网站在线观看| 国产精品三级大全| 毛片一级片免费看久久久久| 国产精品一二三区在线看| 国产日韩欧美在线精品| 亚洲无线观看免费| 欧美日本中文国产一区发布| 人妻夜夜爽99麻豆av| 欧美亚洲 丝袜 人妻 在线| 亚洲国产精品专区欧美| 2018国产大陆天天弄谢| 男人舔奶头视频| 亚洲欧美精品专区久久| 我的老师免费观看完整版| 午夜精品国产一区二区电影| 成年人午夜在线观看视频| av福利片在线观看| 亚洲国产日韩一区二区| 不卡视频在线观看欧美| 内地一区二区视频在线| 欧美 日韩 精品 国产| 久久精品国产亚洲av天美| 亚洲第一av免费看| 亚洲丝袜综合中文字幕| 97超碰精品成人国产| 日日摸夜夜添夜夜添av毛片| 日韩电影二区| 80岁老熟妇乱子伦牲交| 99久久精品一区二区三区| 精品一品国产午夜福利视频| 精品一品国产午夜福利视频| 日韩视频在线欧美| 在线观看免费视频网站a站| 成年女人在线观看亚洲视频| 精品人妻熟女av久视频| 亚洲综合精品二区| 久久精品夜色国产| 亚洲av免费高清在线观看| 中文乱码字字幕精品一区二区三区| 丰满少妇做爰视频| 嫩草影院入口| 欧美精品国产亚洲| 激情五月婷婷亚洲| 18禁在线播放成人免费| 精品久久国产蜜桃| 啦啦啦在线观看免费高清www| 少妇的逼好多水| 男女边吃奶边做爰视频| 免费观看无遮挡的男女| 久久国产亚洲av麻豆专区| 丰满少妇做爰视频| 日韩亚洲欧美综合| 亚洲av免费高清在线观看| 建设人人有责人人尽责人人享有的| 国产视频内射| 啦啦啦啦在线视频资源| 欧美 亚洲 国产 日韩一| 最近手机中文字幕大全| 大片免费播放器 马上看| 午夜免费鲁丝| 午夜免费鲁丝| 亚洲人与动物交配视频| 亚洲精品成人av观看孕妇| 国产精品国产三级国产专区5o| 一级av片app| 日本欧美国产在线视频| 成人18禁高潮啪啪吃奶动态图 | 国产精品麻豆人妻色哟哟久久| 99热这里只有是精品在线观看| 在线免费观看不下载黄p国产| 国产亚洲精品久久久com| 国产男女内射视频| 欧美高清成人免费视频www| 日韩电影二区| 91精品伊人久久大香线蕉| 日韩成人av中文字幕在线观看| 亚洲国产精品一区三区| 国产亚洲最大av| 久久99精品国语久久久| 国产av一区二区精品久久| 久久久久视频综合| 日韩制服骚丝袜av| 在线看a的网站| 久久精品国产亚洲av涩爱| 人人妻人人澡人人爽人人夜夜| 性色av一级| 人人妻人人爽人人添夜夜欢视频 | 亚洲av免费高清在线观看| 五月玫瑰六月丁香| 亚洲精品色激情综合| 国产精品人妻久久久影院| 欧美日韩视频精品一区| 22中文网久久字幕| 欧美人与善性xxx| av女优亚洲男人天堂| 日韩欧美精品免费久久| 王馨瑶露胸无遮挡在线观看| 少妇猛男粗大的猛烈进出视频| 91久久精品国产一区二区成人| 午夜免费鲁丝| 亚洲欧美清纯卡通| 午夜免费鲁丝| 在线观看免费高清a一片| av卡一久久| 插阴视频在线观看视频| 国产片特级美女逼逼视频| 日韩 亚洲 欧美在线| 国产中年淑女户外野战色| 99热全是精品| 婷婷色av中文字幕| 建设人人有责人人尽责人人享有的| 日韩一区二区三区影片| 水蜜桃什么品种好| 妹子高潮喷水视频| 九九爱精品视频在线观看| 欧美97在线视频| 一级毛片aaaaaa免费看小| 国产日韩欧美视频二区| 丰满饥渴人妻一区二区三| 蜜桃久久精品国产亚洲av| 国产免费又黄又爽又色| 久久精品国产鲁丝片午夜精品| 伊人久久国产一区二区| 亚洲第一av免费看| 亚洲欧美成人精品一区二区| 看免费成人av毛片| 伦理电影大哥的女人| 一级a做视频免费观看| 国产欧美亚洲国产| 性高湖久久久久久久久免费观看| 久久久久久久久久久久大奶| 一本—道久久a久久精品蜜桃钙片| 久久99一区二区三区| 久久青草综合色| 国语对白做爰xxxⅹ性视频网站| 欧美另类一区| 一级av片app| 日本免费在线观看一区| 男女国产视频网站| 日韩一区二区三区影片| 女人久久www免费人成看片| 草草在线视频免费看| 国产亚洲最大av| 中文字幕免费在线视频6| 欧美xxⅹ黑人| av天堂久久9| 婷婷色av中文字幕| 最近中文字幕2019免费版| 最近最新中文字幕免费大全7| a 毛片基地| 国产在线视频一区二区| 最近最新中文字幕免费大全7| 在线 av 中文字幕| 免费看不卡的av| 99热国产这里只有精品6| 一本久久精品| 少妇的逼水好多| 午夜av观看不卡| 成人无遮挡网站| h视频一区二区三区| 午夜av观看不卡| 国产精品久久久久成人av| 成人美女网站在线观看视频| av网站免费在线观看视频| 亚洲av电影在线观看一区二区三区| 韩国高清视频一区二区三区| 亚洲在久久综合| 午夜av观看不卡| 男男h啪啪无遮挡| 亚洲欧美中文字幕日韩二区| 免费黄频网站在线观看国产| 尾随美女入室| 五月伊人婷婷丁香| 国产成人a∨麻豆精品| 欧美日韩av久久| 国产美女午夜福利| 精品一品国产午夜福利视频| 日日啪夜夜撸| 久久午夜综合久久蜜桃| 久久精品久久久久久久性| 国产黄片美女视频| 成人国产av品久久久| 女的被弄到高潮叫床怎么办| 99国产精品免费福利视频| 国产极品粉嫩免费观看在线 | 久久人人爽人人爽人人片va| 国产精品熟女久久久久浪| 美女中出高潮动态图| 大话2 男鬼变身卡| 日韩av免费高清视频| 蜜臀久久99精品久久宅男| 我要看日韩黄色一级片| 男人狂女人下面高潮的视频| 国产毛片在线视频| 十分钟在线观看高清视频www | 蜜桃久久精品国产亚洲av| 人体艺术视频欧美日本| 国产精品偷伦视频观看了| 久久久久久伊人网av| 国产综合精华液| 亚洲国产色片| 男女无遮挡免费网站观看| 在线观看免费日韩欧美大片 | 免费大片18禁| 少妇 在线观看| 欧美成人午夜免费资源| av在线老鸭窝| 日日啪夜夜爽| 免费大片黄手机在线观看| 人人妻人人澡人人看| 国产av国产精品国产| 中文在线观看免费www的网站| 国产熟女欧美一区二区| 亚洲伊人久久精品综合| 国语对白做爰xxxⅹ性视频网站| 久久这里有精品视频免费| 22中文网久久字幕| 久久久久人妻精品一区果冻| 黄色视频在线播放观看不卡| 十八禁网站网址无遮挡 | 丰满迷人的少妇在线观看| 能在线免费看毛片的网站| 日本欧美视频一区| freevideosex欧美| 久久综合国产亚洲精品| 秋霞在线观看毛片| 亚洲,欧美,日韩| 国产极品天堂在线| 国产精品秋霞免费鲁丝片| 亚洲,一卡二卡三卡| 春色校园在线视频观看| 内地一区二区视频在线| 成人亚洲精品一区在线观看| 免费人妻精品一区二区三区视频| 亚洲内射少妇av| 久久久久精品性色| 中文字幕久久专区| 天天躁夜夜躁狠狠久久av| 伦理电影免费视频| 国产黄色视频一区二区在线观看| h视频一区二区三区| freevideosex欧美| 国产高清有码在线观看视频| 亚洲精品456在线播放app| 黄色日韩在线| 亚洲欧美成人综合另类久久久| 激情五月婷婷亚洲| 乱系列少妇在线播放| 亚洲人成网站在线观看播放| 日韩av不卡免费在线播放| 国产精品免费大片| 国产成人aa在线观看| 久久久精品免费免费高清| 日日摸夜夜添夜夜添av毛片| 亚洲精品日韩av片在线观看| 久久6这里有精品| av在线app专区| av天堂中文字幕网| 女人久久www免费人成看片| 久久久久久久久久久久大奶| 亚洲欧洲国产日韩| 久久人人爽av亚洲精品天堂| 成人二区视频| 国产伦在线观看视频一区| 一区二区三区精品91| av天堂久久9| 91精品伊人久久大香线蕉| 九九爱精品视频在线观看| 欧美xxⅹ黑人| 黑人猛操日本美女一级片| 91精品国产九色| 亚洲三级黄色毛片| 两个人的视频大全免费| 久久影院123| 亚洲av综合色区一区| 成人免费观看视频高清| 亚洲国产av新网站| 在现免费观看毛片| 久久毛片免费看一区二区三区| 国产熟女午夜一区二区三区 | 六月丁香七月| 亚洲精品成人av观看孕妇| 在线天堂最新版资源| 亚洲精品国产色婷婷电影| 三上悠亚av全集在线观看 | 国产成人免费观看mmmm| 亚洲精品日本国产第一区| 国产精品一区二区三区四区免费观看| 亚洲国产av新网站| 最后的刺客免费高清国语| 欧美 日韩 精品 国产| 欧美区成人在线视频| 日韩大片免费观看网站| 亚洲国产色片| 人妻 亚洲 视频| 国产一区二区三区av在线| 久久久久久久久大av| 亚洲va在线va天堂va国产| 韩国高清视频一区二区三区| 国产精品国产av在线观看| 国产精品成人在线| 一级二级三级毛片免费看| 精品少妇内射三级| a级毛色黄片| 国产伦精品一区二区三区四那| av福利片在线观看| 亚洲欧美精品专区久久| 狂野欧美激情性xxxx在线观看| 久久久久网色| 精品一区二区三区视频在线| 观看免费一级毛片| 美女脱内裤让男人舔精品视频| 日韩av在线免费看完整版不卡| www.色视频.com| 狂野欧美白嫩少妇大欣赏| 中文天堂在线官网| 嫩草影院新地址| 国产深夜福利视频在线观看| 精品国产一区二区三区久久久樱花| 在线观看三级黄色| 性色avwww在线观看| av视频免费观看在线观看| 日本与韩国留学比较| 观看美女的网站| 热re99久久精品国产66热6| 国产精品一区二区在线观看99| 亚洲国产av新网站| 国产日韩欧美视频二区| 亚洲欧美成人精品一区二区| 日日撸夜夜添| 国产淫片久久久久久久久| 欧美区成人在线视频| 国产亚洲5aaaaa淫片| 一个人看视频在线观看www免费| 少妇被粗大猛烈的视频| 亚洲综合精品二区| 欧美高清成人免费视频www| 看非洲黑人一级黄片| 一本一本综合久久| 日韩欧美一区视频在线观看 | 色哟哟·www| 国产精品国产三级国产专区5o| 熟女av电影| 哪个播放器可以免费观看大片| 亚洲国产最新在线播放| 亚洲精品一二三| 只有这里有精品99| 日日爽夜夜爽网站| 日韩视频在线欧美| 国产 一区精品| 人人澡人人妻人| 成年av动漫网址| 欧美最新免费一区二区三区| 国产精品免费大片| 嘟嘟电影网在线观看| 寂寞人妻少妇视频99o| 黑人猛操日本美女一级片| 极品教师在线视频| 久久久国产精品麻豆| xxx大片免费视频| 久久精品久久久久久久性| 99视频精品全部免费 在线| 男女国产视频网站| 亚洲国产成人一精品久久久| 色网站视频免费| 久久影院123| 内地一区二区视频在线| 天天操日日干夜夜撸| 日本av手机在线免费观看| 在线观看www视频免费| 一个人看视频在线观看www免费| 一本—道久久a久久精品蜜桃钙片| 亚洲国产精品成人久久小说| a级片在线免费高清观看视频| 最黄视频免费看| 好男人视频免费观看在线| 午夜福利,免费看| 大片电影免费在线观看免费| 在线观看免费高清a一片| 国产爽快片一区二区三区| 久久久久久久久久久久大奶| 久久人人爽人人片av| 国产精品熟女久久久久浪| 尾随美女入室| 久久久久久伊人网av| 国产日韩欧美视频二区| 97超碰精品成人国产| 春色校园在线视频观看| 色5月婷婷丁香| 最新中文字幕久久久久| 国产精品久久久久久久久免| 国产一区二区三区综合在线观看 | 热re99久久精品国产66热6| 亚洲av欧美aⅴ国产| 国产精品熟女久久久久浪| 人人妻人人澡人人爽人人夜夜| 亚洲精品久久午夜乱码| 欧美成人午夜免费资源| 99热国产这里只有精品6| 男人爽女人下面视频在线观看| 久久久精品免费免费高清| 欧美亚洲 丝袜 人妻 在线| 六月丁香七月| 男女无遮挡免费网站观看| 大香蕉久久网| √禁漫天堂资源中文www| av女优亚洲男人天堂| 五月开心婷婷网| 美女视频免费永久观看网站| 一级毛片黄色毛片免费观看视频| 中文字幕人妻丝袜制服| 久久99热这里只频精品6学生| 成人综合一区亚洲| 自拍偷自拍亚洲精品老妇| 久久久久人妻精品一区果冻| 亚洲人成网站在线观看播放| 看免费成人av毛片| 国产男人的电影天堂91| 丝瓜视频免费看黄片| 亚洲中文av在线| 免费大片18禁| 美女视频免费永久观看网站| 亚洲成色77777| 亚洲色图综合在线观看| 国产欧美日韩精品一区二区| www.色视频.com| 视频区图区小说| 日本爱情动作片www.在线观看| 久久这里有精品视频免费| 国产 一区精品| 欧美丝袜亚洲另类| 久久免费观看电影| 亚洲国产精品成人久久小说| 老司机影院毛片| 国产男女内射视频| 精品一区二区三卡| 91aial.com中文字幕在线观看| 亚洲人与动物交配视频| 人体艺术视频欧美日本| 国产精品人妻久久久影院| 黑人巨大精品欧美一区二区蜜桃 | 日本-黄色视频高清免费观看| 一本色道久久久久久精品综合| 亚洲欧美日韩另类电影网站| 国产成人精品一,二区| 日本wwww免费看| av不卡在线播放| 国产精品一区二区三区四区免费观看| 亚洲自偷自拍三级| 自线自在国产av| 日韩在线高清观看一区二区三区| 久久久久久久久大av| 在线观看一区二区三区激情| 免费不卡的大黄色大毛片视频在线观看| 性色avwww在线观看| 伦精品一区二区三区| 亚洲一级一片aⅴ在线观看| 日韩欧美精品免费久久| 亚洲av成人精品一区久久| 精品一区二区免费观看| 国语对白做爰xxxⅹ性视频网站| 黄色一级大片看看| 精品久久久久久久久av| 免费在线观看成人毛片| 国产成人精品无人区| 日日摸夜夜添夜夜爱| 国产精品不卡视频一区二区| 欧美精品人与动牲交sv欧美| 一区二区三区乱码不卡18| 久久这里有精品视频免费| 99久久中文字幕三级久久日本| 简卡轻食公司| 国产高清不卡午夜福利| 精品少妇久久久久久888优播| 国产日韩欧美亚洲二区| 中文欧美无线码| 久久精品国产亚洲网站| 全区人妻精品视频| 午夜老司机福利剧场| h视频一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品一二三| 亚洲欧美成人综合另类久久久| 亚洲成人一二三区av| 亚洲欧美成人综合另类久久久| 欧美成人午夜免费资源| 丰满人妻一区二区三区视频av| 国产一区有黄有色的免费视频| 一本色道久久久久久精品综合| 男女无遮挡免费网站观看| 国产精品国产三级国产av玫瑰| 涩涩av久久男人的天堂| 岛国毛片在线播放| 中文字幕精品免费在线观看视频 | 成人无遮挡网站| 日本黄大片高清| 建设人人有责人人尽责人人享有的| 成人亚洲精品一区在线观看| 日韩中字成人| 下体分泌物呈黄色| 国产精品福利在线免费观看| 99热6这里只有精品| 日韩伦理黄色片| 国产免费又黄又爽又色| 九色成人免费人妻av| 亚洲国产日韩一区二区| 国产成人a∨麻豆精品| 国产探花极品一区二区| 一本久久精品| freevideosex欧美| 乱系列少妇在线播放| 成人国产麻豆网| 日韩电影二区| 一级毛片 在线播放| 99热这里只有是精品50| 又粗又硬又长又爽又黄的视频| 亚洲av不卡在线观看| 成人黄色视频免费在线看| 亚洲欧美中文字幕日韩二区| 国产午夜精品久久久久久一区二区三区| 在线免费观看不下载黄p国产| 看非洲黑人一级黄片| 新久久久久国产一级毛片| av在线观看视频网站免费| 99热这里只有是精品在线观看| 亚洲av欧美aⅴ国产| 嫩草影院入口|