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

    基于PFC3D的加權(quán)平均數(shù)建模法在三軸數(shù)值試驗中的應(yīng)用*

    2017-06-05 14:21:32張志華張謝東邱紅勝鄧雅思
    關(guān)鍵詞:粒土細(xì)觀半徑

    張志華 張謝東 邱紅勝 武 林 鄧雅思

    (武漢理工大學(xué)交通學(xué)院 武漢 430063)

    基于PFC3D的加權(quán)平均數(shù)建模法在三軸數(shù)值試驗中的應(yīng)用*

    張志華 張謝東 邱紅勝 武 林 鄧雅思

    (武漢理工大學(xué)交通學(xué)院 武漢 430063)

    基于室內(nèi)三軸試驗研究成果,以PFC3D為工具建立三維顆粒流模型,開展三軸數(shù)值試驗研究.針對大型三軸數(shù)值試驗大范圍粒徑顆粒生成困難問題,提出了加權(quán)平均數(shù)數(shù)值建模方法,并以轉(zhuǎn)動顆粒場分析手段為依據(jù)進行局部裂隙細(xì)觀分析.研究表明,加權(quán)平均數(shù)數(shù)值建模法得到的顆粒級配與室內(nèi)試驗顆粒級配及三軸試驗結(jié)果吻合度都較好,此法有利于提高計算效率;轉(zhuǎn)動顆粒場有利于三維立體模型剪切面的細(xì)觀分析,將其與顆粒位移場相結(jié)合進行細(xì)觀分析能更好的揭示粗粒土宏觀力學(xué)及變形特性的內(nèi)在機制.

    顆粒流模型;三軸試驗;加權(quán)平均數(shù)顆粒生成法;轉(zhuǎn)動顆粒場

    0 引 言

    由于粗粒土中存在大粒徑顆粒,一般的小型三軸試驗很難對粗粒土進行試驗研究,而大型三軸試驗試樣的體積較大,試驗研究耗時耗力,且對于大型三軸試驗的數(shù)值模型的建立又是一大難題.對于建立成功的模型,細(xì)觀分析手段的選用是揭示粗粒土宏觀特性成敗的關(guān)鍵.Cundall等[1-2]提出了離散元理論,而顆粒流(PFC)作為離散元方法的一種,解決顆粒體非連續(xù)介質(zhì)數(shù)值模擬問題,并將顆粒流的思想應(yīng)用到土體中.

    近幾年來,國內(nèi)外學(xué)者對顆粒流理論運用越來越多.Feng等[3-5]給出了二維區(qū)域內(nèi)部圓盤顆粒填充的前進邊法,而對于三維區(qū)域內(nèi)球體顆粒的前進面法并未給出具體的實現(xiàn)過程;耿慶東[6]提出了三維區(qū)域內(nèi)球體顆粒生成的前進面法,確定了前進面的選取順序,還通過改變半徑解決了漏填的情況,得到了三維區(qū)域內(nèi)部球體的高效填充方法;童朝霞等[7]采用截取法制樣,先重力沉積大量的密集顆粒體,再從中截取與室內(nèi)試樣尺寸相同的矩形試樣進行雙軸試驗;周劍等[8]對真實物體進行數(shù)字圖像二值處理,得到反映土石混合體結(jié)構(gòu)的圖像,將數(shù)字圖像轉(zhuǎn)換為PFC模型;Liu等[9]探討了不同黏結(jié)強度下的屈服應(yīng)力導(dǎo)致顆粒間裂隙的產(chǎn)生和發(fā)展不同;以細(xì)小顆粒簇模擬不規(guī)則非球狀的顆粒體進行數(shù)值模擬試驗,數(shù)值結(jié)果吻合程度較高,但此方法耗時耗力[10-12].

    基于此,文中提出了加權(quán)平均數(shù)顆粒生成法進行大型三軸數(shù)值試驗?zāi)M,并根據(jù)擬合室內(nèi)三軸試驗結(jié)果為依據(jù),以轉(zhuǎn)動顆粒場為細(xì)觀分析手段,并與顆粒位移場等分析手段相結(jié)合,進行試樣局部破壞過程分析,對今后粗粒土三軸顆粒流仿真試驗研究提供一定的參考價值.

    1 顆粒轉(zhuǎn)動模型

    在顆粒流理論中,顆粒間的相互作用看成一個動態(tài)的過程,在這個過程中顆粒內(nèi)部之間的力是平衡的,系統(tǒng)處于基本平衡狀態(tài)在運行.顆粒流理論運用的位移法與顆粒-顆粒間接觸點的相對位移和接觸力有關(guān),球與球的接觸,在平行黏結(jié)模型中,隨著黏結(jié)劑發(fā)生變形,接觸點處會產(chǎn)生額外的接觸力和力矩;在接觸黏結(jié)模型中,接觸點處只會產(chǎn)生接觸力,由于接觸黏結(jié)劑不能傳遞力矩而不會產(chǎn)生力矩.

    但是顆粒間的相互作用力大于黏結(jié)劑強度時黏結(jié)劑破裂,顆粒間形成裂隙.此時的顆粒間相互作用力就發(fā)生了改變,顆粒轉(zhuǎn)動模型見圖1.

    圖1 顆粒轉(zhuǎn)動模型

    (1)

    (2)

    (3)

    2 三軸試驗顆粒流模擬

    2.1 粗粒土三軸試驗宏觀力學(xué)參數(shù)

    試驗使用天水洪山公司制造的1 500kN微機控制電液伺服靜態(tài)三軸試驗機,試驗用料取自水布埡大壩工程當(dāng)年的施工料場,是水布埡茅口組灰?guī)r堆石料.參照《土工試驗規(guī)程》(SL237—1999)對土的飽和吸水率、飽和時極限抗壓強度、比重、軟化系數(shù)等基本物性指標(biāo)進行測試,得到其物理性能.結(jié)果見表1.

    表1 粗粒土的基本物性指標(biāo)

    試驗三軸試樣的尺寸直徑×高度為300 mm×600 mm,干密度為2 180 kg/m3,得到一個試樣的總重量為92.457 kg.采用篩分法和混合法[13]測試粗粒土的顆粒級配,得到粗粒土不同粒徑范圍內(nèi)顆粒的含量見表2.

    表2 粗粒土室內(nèi)試驗顆粒級配

    2.2 加權(quán)平均數(shù)顆粒生成法

    在進行顆粒生成時,如果按照試驗室內(nèi)的土料級配進行生成顆粒時,會發(fā)現(xiàn)一個很嚴(yán)重的問題,就是當(dāng)生成顆粒粒徑小于1 mm的顆粒時,雖然在該粒徑下的土料所占比例較小,但是由于單個顆粒的質(zhì)量很小,導(dǎo)致在顆粒流模型中會生成大量的粉末狀的小粒徑顆粒,嚴(yán)重影響了模型的運算速度.

    加權(quán)平均數(shù)法是根據(jù)室內(nèi)試驗材料的顆粒級配作為數(shù)據(jù)基礎(chǔ),利用若干個粒徑范圍內(nèi)的顆粒含量作為權(quán)數(shù)并與相應(yīng)的顆粒粒徑的乘積之和,計算得到的該加權(quán)算術(shù)平均數(shù).即

    (4)

    式中:Ri為顆粒平均半徑;Mi為某一粒徑下的顆粒含量(質(zhì)量分?jǐn)?shù));Wi為某一粒徑的加權(quán)數(shù),Wr為顆粒粒徑的加權(quán)平均數(shù);∑Mi為各顆粒粒徑占總質(zhì)量的含量之和,其值為1.

    根據(jù)室內(nèi)顆粒級配(見表2)和式(3),計算得到室內(nèi)試驗一個試樣內(nèi)的顆粒平均半徑為11.066 25mm.在顆粒流模型中,如果按照顆粒最小半徑及最大最小顆粒半徑比進行顆粒生成時,根據(jù)系統(tǒng)中的Random口令,顆粒半徑從最小到最大的生成概率是相同的.依據(jù)室內(nèi)試驗的加權(quán)平均顆粒半徑,如果固定顆粒的最小半徑,即可算出顆粒的最大半徑.計算式為

    (5)

    式中:Ri為表2中的顆粒半徑;ratioRmax/Rmin為顆粒最大半徑與最小半徑之比;Rmin為最小顆粒半徑.根據(jù)式(4)算得的Wr,即可算ratioRmax/Rmin的值.

    根據(jù)室內(nèi)試驗顆粒級配,選用模型內(nèi)最小顆粒半徑為0.5 mm,則根據(jù)式(4)得出,模型內(nèi)顆粒最大半徑與最小半徑之比為59.265,因此最大顆粒半徑為29.632 5 mm,接近室內(nèi)試驗試樣的最大顆粒尺寸.根據(jù)隨機生成概率相同的原理,得到顆粒流模型內(nèi)顆粒的級配,見表3,并進行室內(nèi)顆粒繼配及顆粒流模型中顆粒繼配曲線對比,見圖2.

    表3 顆粒流模型內(nèi)顆粒級配

    圖2 室內(nèi)試驗及數(shù)值試驗顆粒級配曲線

    由圖2可知,加權(quán)平均數(shù)建模法生成的顆粒級配與室內(nèi)試驗顆粒級配基本吻合.

    2.3 各細(xì)觀參數(shù)的確定

    根據(jù)室內(nèi)試驗結(jié)果的數(shù)據(jù)處理,并結(jié)合顆粒流模型內(nèi)考慮到的各細(xì)觀參數(shù)如初始彈性模量Ec,顆??辜魪姸萩,φ,顆粒剛度比kn/ks,顆粒最小半徑Rmin,顆粒最大、最小半徑比ratioRmax/Rmin;試樣的密度ρ,初始接觸點個數(shù)Nf,墻面的剛度Wstiff及墻面的摩擦系數(shù)Wfric等.文中根據(jù)室內(nèi)試驗結(jié)果中應(yīng)力應(yīng)變曲線的彈性階段(起始直線部分)計算初始彈性模量,根據(jù)摩爾應(yīng)力圓得到顆粒的抗剪強度指標(biāo),根據(jù)試驗試樣確定材料的密度,基于加權(quán)平均數(shù)計算方法計算得到數(shù)值試驗顆粒級配曲線中各顆粒粒徑的含量.得到一系列細(xì)觀參數(shù)值,見表4.

    表4 三軸試驗?zāi)P椭懈骷?xì)觀參數(shù)(圍壓0.8 MPa)

    2.4 數(shù)值結(jié)果與試驗結(jié)果對比

    根據(jù)表4中確定的模型內(nèi)的各細(xì)觀參數(shù)的值,并根據(jù)PFC3D手冊[14-15]中不同圍壓下模型內(nèi)顆粒的初始楊氏模量和顆粒剛度比不同.根據(jù)泊松比與顆粒剛度比之間的關(guān)系[16],在已知室內(nèi)試驗結(jié)束時測得試樣的泊松比,可計算各圍壓下的顆粒剛度比.圖3為不同圍壓下數(shù)值試驗應(yīng)力應(yīng)變曲線與室內(nèi)試驗結(jié)果的對比圖.

    圖3 不同圍壓下室內(nèi)試驗與數(shù)值試驗應(yīng)力應(yīng)變曲線

    3 以0.8 MPa圍壓為例的局部破壞過程分析

    在室內(nèi)三軸試驗中,低圍壓下試樣容易發(fā)生剪切破壞,而在高圍壓下試樣則不容易發(fā)生剪切破壞.但由圖3中應(yīng)力應(yīng)變曲線可知,在不同圍壓下,試樣在軸應(yīng)變5%~10%之間發(fā)生了從彈性變形到塑性變形的過渡.在這個過程中試樣內(nèi)部是否發(fā)生了剪切破壞,這就需要從細(xì)觀角度進行分析.

    為了分析數(shù)值試驗及室內(nèi)三軸試驗局部破壞過程,提出轉(zhuǎn)動顆粒場分析手段,并與顆粒位移場進行分析,對比模型剪切面的形成機理.

    3.1 不同的細(xì)觀分析方法

    3.1.1 轉(zhuǎn)動顆粒場

    依照位移法提出了轉(zhuǎn)動顆粒場細(xì)觀分析手段,根據(jù)局部破壞過程中的顆粒轉(zhuǎn)角大小分布情況,找出了二維顆粒流雙軸試驗剪切帶的分布情況,表明在模型加載過程中顆粒材料發(fā)生剪切破壞時局部顆粒伴隨有較大的轉(zhuǎn)動.采用X射線及光學(xué)測量手段分析了砂土的微觀結(jié)構(gòu)變形,發(fā)現(xiàn)在局部化砂土破壞形成過程中,局部區(qū)域的顆粒顯示出較大的顆粒旋轉(zhuǎn).

    圖4 不同軸應(yīng)變時的轉(zhuǎn)動顆粒場

    圖5給出了不同軸應(yīng)變時的顆粒旋轉(zhuǎn)速度場,當(dāng)顆粒的瞬時轉(zhuǎn)動速度達(dá)到110 rad/s時(以0.8 MPa為例,粘結(jié)顆粒發(fā)生破裂時,接觸力均達(dá)到兆帕,顆粒間接觸力很大,由式(3)可知,質(zhì)量很小的顆粒在較大力矩作用下產(chǎn)生的瞬時角加速度很大),顆粒被標(biāo)定為黑色,因此可以明顯的看出隨著軸應(yīng)變的增大顆粒在局部區(qū)域內(nèi)的旋轉(zhuǎn)速度情況,尤其是在剪切面附近的顆粒.因此,隨著剪切面的形成,剪切面附近的顆粒伴隨有較大的轉(zhuǎn)動.

    圖5 不同軸應(yīng)變時的顆粒位移場

    由圖4可知,軸應(yīng)變在5%~10%內(nèi)模型中轉(zhuǎn)動速度較大的顆粒數(shù)目較少,在軸應(yīng)變達(dá)到15%時,從外觀上可以很明顯的看到黑色顆粒的斜面分布規(guī)律,該斜面為剪切面.且隨著軸應(yīng)變不斷增大,分布在剪切面附近的顆粒數(shù)目越來越多.通過轉(zhuǎn)動速度較大的顆粒可以很好的找出模型中的剪切面的分布情況.比較圖3~4可知,試樣在壓縮剪切過程中達(dá)到應(yīng)力峰值之后的應(yīng)變軟化階段,剪切面才會完全形成.

    3.1.2 顆粒位移場

    觀察不同軸應(yīng)變時的顆粒位移場發(fā)現(xiàn),見圖5a),從壓縮剪切開始,由于上層加載板向下運動,導(dǎo)致模型上半部分顆粒集體向下運動,下半部分顆粒呈現(xiàn)兩側(cè)運動.在軸應(yīng)變5%左右時模型上半部分顆粒的位移分布比較均勻,上層顆粒均勻向下移動,基本沒有分塊移動;當(dāng)軸應(yīng)變達(dá)到10%時,上半部分顆粒開始出現(xiàn)一個三角區(qū)域,見圖5b),該三角區(qū)域是試樣剪切破壞之后的塊體,剪切面分布在三角區(qū)域的下邊緣.也就是說當(dāng)軸應(yīng)變在10%左右,試樣出現(xiàn)剪切面初始結(jié)構(gòu),試樣處于彈性到塑性階段的過渡過程,應(yīng)力應(yīng)變曲線出現(xiàn)應(yīng)力峰值,與圖3中室內(nèi)試驗結(jié)果吻合.在軸應(yīng)變?yōu)?5%~25%階段時,見圖5c)~圖5e),可以很明顯的看到模型右上角三角區(qū)域顆粒位移整體呈現(xiàn)右下方向,剪切面徹底形成,試樣處于塑性階段,應(yīng)力應(yīng)變曲線出現(xiàn)應(yīng)變軟化現(xiàn)象.因此剪切面的完整形成是在試樣應(yīng)力應(yīng)變曲線達(dá)到峰值強度之后的軟化階段.

    3.2 轉(zhuǎn)動顆粒與空間剪切面

    由于三軸試驗三維立體的內(nèi)部不可視性,導(dǎo)致形成的空間剪切面也無法看見.文中提出的轉(zhuǎn)動顆粒場以轉(zhuǎn)動顆粒跟蹤試樣空間剪切破壞面,比較圖5~6可以確定轉(zhuǎn)動顆粒場的可用性.基于此,文中在提出轉(zhuǎn)動顆粒場的基礎(chǔ)上,根據(jù)數(shù)值模型顆粒的空間處理,得到去除轉(zhuǎn)動顆粒的試樣空間剪切帶和轉(zhuǎn)動顆粒的空間分布圖,見圖7.

    圖6 去除轉(zhuǎn)動顆粒的空間剪切帶空隙

    圖7 不同軸應(yīng)變時的轉(zhuǎn)動顆粒空間分布

    圖6是去掉轉(zhuǎn)動顆粒之后的顆粒模型,可以很直觀的看出圍壓0.8 MPa下試樣的剪切破壞“面”,該破壞面并不是一個簡簡單單的單顆粒厚度的凸凹面,而是具有一定厚度,且該剪切面的厚度在模型邊緣較薄,在模型內(nèi)部較厚.

    圖7表示轉(zhuǎn)動顆粒在模型剪切壓縮過程中的生成與分布,由圖7a)~圖7b)可知,在初始壓縮階段,試樣內(nèi)部出現(xiàn)了大量的轉(zhuǎn)動顆粒,這是試樣內(nèi)產(chǎn)生局部裂隙導(dǎo)致的.圖7c)可知,轉(zhuǎn)動顆粒的分布呈斜面分布,與圖5c)比較得知,該斜面是剪切面.因此,在軸應(yīng)變10%~15%(含15%),某一時刻試樣被剪破,這與模型宏觀力學(xué)反應(yīng)-應(yīng)力應(yīng)變曲線是吻合的.

    將轉(zhuǎn)動顆粒與其位移場結(jié)合見圖8,可以看出轉(zhuǎn)動顆粒的位移場與顆粒的總位移場(圖6)是一致的.轉(zhuǎn)動顆粒在轉(zhuǎn)動過程中沿著剪切面的走勢滾動,滾動的顆粒更加促使了剪切面上相對摩擦力的減小,使得宏觀力學(xué)特性表現(xiàn)出抗剪強度下降的現(xiàn)象,從而形成了峰值過后應(yīng)變軟化的現(xiàn)象.

    圖9表示在三軸數(shù)值試驗過程中轉(zhuǎn)動顆粒的數(shù)量,得到轉(zhuǎn)動顆粒數(shù)量與軸應(yīng)變的關(guān)系式為一元二次方程,由于顆粒的轉(zhuǎn)動是顆粒間黏結(jié)劑破裂導(dǎo)致的,表明在剪切破壞初期,模型內(nèi)裂隙數(shù)量呈線性增加,當(dāng)試樣剪切破壞以后,模型內(nèi)部顆粒間裂隙的數(shù)量增加減緩,這是模型內(nèi)形成了剪切面造成的.剪切面的形成使得分布在剪切面上下的顆粒體會沿著最小抗剪強度的位置移動,也導(dǎo)致了模型整體抗剪強度的下降.

    圖8 軸應(yīng)變?yōu)?0%時的轉(zhuǎn)動顆粒位移場

    圖9 轉(zhuǎn)動顆粒的數(shù)量與軸應(yīng)變的關(guān)系

    4 結(jié) 論

    1) 以加權(quán)平均數(shù)顆粒生成法建立粗粒土三軸顆粒流模型,顆粒級配與室內(nèi)試樣級配大致吻合,且得到的應(yīng)力應(yīng)變曲線與室內(nèi)試驗結(jié)果基本吻合,驗證了該方法在三軸顆粒流數(shù)值試驗中的可用性.

    2) 轉(zhuǎn)動顆粒場分析手段能很好的揭示低圍壓下粗粒土應(yīng)變軟化現(xiàn)象機理,并將其與顆粒位移場相結(jié)合進行細(xì)觀分析,能更好的揭示粗粒土三軸試驗宏觀力學(xué)特性機理.尤其是以轉(zhuǎn)動顆粒為分析對象,充分解釋了低圍壓下剪切破壞初期模型內(nèi)部裂隙數(shù)量線性增加,裂隙貫穿形成剪切面的過程,而且根據(jù)轉(zhuǎn)動顆粒空間位置分布,得出剪切破壞面是一個厚度不均勻,形狀不規(guī)則的“斜面板”.

    [1]CUNDALL P A. A computer model for simulating progressive large scale movements in blocky rock systems[J]. Proceedings of the Symposium of the International Society of Rock Mechanics,1971(1):2-8.

    [2]CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Geotechnique,1979,29(1):47-65.

    [3]FENG Y T, HAN K, OWEN D R J. Filling domains with disks: an advancing front approach[J]. International Journal for Numerical Methods in Engineering,2003,56(3):699-713.

    [4]FENG Y T, HAN K, OWEN D R J. Filling domains with disks[D]. Glasgow: University of Glasgow,2001.

    [5]FENG Y T, HAN K, OWEN D R J. An advancing front packing of polygons, ellipses and spheres[C]. Proceedings of the 3rdinternational Conference on Diserete Element Methods: Numercial Modeling of Discontinua,Santa Fe, New Mexico, USA,2002.

    [6]耿慶東.公共面搜索算法改進及三維顆粒體生成算法[D].大連:大連理工大學(xué),2009.

    [7]童朝霞,周敏,張連衛(wèi),等.各項異性顆粒材料雙軸壓縮試驗的離散元數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報,2014,33(增刊2):4227-4232.

    [8]周劍,張路青,戴福初,等.基于黏結(jié)顆粒模型某滑坡土石混合體直剪試驗數(shù)值模擬[J].巖石力學(xué)與工程學(xué)報,2013,32(增刊1):2650-2659.

    [9]LIU S H, SUN D A, WANG Y S. Numerical study of soil collapse behavior by discrete element modeling[J]. Computers and Geotechnics,2003,30:399-408.

    [10]PRICE M, MURARIU V, MORRISON G. Sphere clump generation and trajectory comparison for real particles[J]. Proceedings of Discrete Element Modelling,2007,15:245-256.

    [11]MATSUSHIMA T. Discrete element simulation of an assembly of irregularly shaped grains: quantitative comparison with experiments[C]. 16th ASCE Engineering Mechanics Conference, University of Washington, Seattle,2003.

    [12]LEE Y. A packing algorithm for three-dimensional convex particles[J]. Granular Matter,2009,11(5):307-315.

    [13]李廣信.高等土力學(xué)[M].北京:清華大學(xué)出版社,2004.

    [14]Itasca Consulting Group Incorporation. PFC3Dtheory and background[R]. Minneapolis:Minnesota,2004.

    [15]Itasca Consulting Group Incorporation. PFC3Duser’s guide[R]. Minneapolis:Minnesota,2004.

    [16]張志華.基于PFC3D的粗粒土三軸試驗細(xì)觀數(shù)值模擬[D].宜昌:三峽大學(xué),2015.

    Weighted Average Assembling Method Using in Triaxial Numerical Experiment Based on PFC3D

    ZHANG Zhihua ZHANG Xiedong QIU Hongsheng WU Lin DENG Yasi

    (SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)

    Based on the results of indoor triaxial experimental research, a three-dimensional particle flow model is established based on PFC3Dto conduct triaxial numerical experiments. Based on the problem of a wide range of particle size in large-scale triaxial numerical test, numerical modeling method of weighted average assembling is proposed. Based on the analysis method of rotating particles field, local fracture microstructure analysis has been conducted. The results show that the particle gradations from laboratory agree well with the simulation results generated by weighted average generate method. In addition, the results of macro responses have great agreement. This method is beneficial to improving the calculation efficiency. Besides, the field of rotating particles in the shear mesoscopic analysis of three-dimensional model is advantageous. By combining the displacement field, the particle mesoscopic analysis can better reveal the internal mechanism of mechanics and deformation characteristics of coarse grained soil.

    particle flow model; triaxial test; weighted average generate method; field of rotating particles

    2016-02-06

    *國家自然科學(xué)基金項目資助(51308429)

    TU43

    10.3963/j.issn.2095-3844.2017.02.017

    張志華(1991—):男,博士生,主要研究領(lǐng)域為巖土體及隧道圍巖穩(wěn)定性離散元數(shù)值仿真

    猜你喜歡
    粒土細(xì)觀半徑
    粗粒土三軸試驗力學(xué)與長期變形特性研究綜述
    無膠結(jié)粗粒土初始結(jié)構(gòu)強度研究
    顆粒形狀對粗粒土破碎特性影響研究
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    連續(xù)展成磨削小半徑齒頂圓角的多刀逼近法
    一些圖的無符號拉普拉斯譜半徑
    加筋粗粒土大型直剪試驗及本構(gòu)模型適用性研究
    熱采水平井加熱半徑計算新模型
    基于四叉樹網(wǎng)格加密技術(shù)的混凝土細(xì)觀模型
    PBX炸藥的抗壓強度及抗拉強度細(xì)觀尺度的數(shù)值計算
    亚洲一区二区三区不卡视频| 亚洲一区高清亚洲精品| 天堂√8在线中文| 亚洲五月色婷婷综合| 免费久久久久久久精品成人欧美视频| 欧美日韩乱码在线| 国产野战对白在线观看| 香蕉久久夜色| 自拍欧美九色日韩亚洲蝌蚪91| 两性夫妻黄色片| 日日干狠狠操夜夜爽| 99热只有精品国产| 国产又色又爽无遮挡免费看| 国产欧美日韩一区二区三区在线| 国产免费男女视频| 亚洲情色 制服丝袜| 麻豆国产av国片精品| 亚洲精品中文字幕一二三四区| av在线播放免费不卡| 天天躁狠狠躁夜夜躁狠狠躁| 波多野结衣高清无吗| 午夜激情av网站| 亚洲一区二区三区不卡视频| 在线国产一区二区在线| 50天的宝宝边吃奶边哭怎么回事| 高清毛片免费观看视频网站| 淫秽高清视频在线观看| 国产欧美日韩一区二区精品| 搡老妇女老女人老熟妇| 亚洲狠狠婷婷综合久久图片| 国产成年人精品一区二区| av欧美777| 欧美av亚洲av综合av国产av| 国语自产精品视频在线第100页| 咕卡用的链子| 人人澡人人妻人| 日韩欧美国产在线观看| 多毛熟女@视频| 国产成人啪精品午夜网站| 国产成人一区二区三区免费视频网站| 久久天躁狠狠躁夜夜2o2o| 欧美日韩乱码在线| 少妇裸体淫交视频免费看高清 | 日韩欧美一区二区三区在线观看| 最好的美女福利视频网| 国产精品影院久久| 欧美亚洲日本最大视频资源| 亚洲中文字幕一区二区三区有码在线看 | 国产精品,欧美在线| 纯流量卡能插随身wifi吗| 国产精品综合久久久久久久免费 | 午夜福利欧美成人| 成人国产一区最新在线观看| 亚洲 欧美 日韩 在线 免费| 一边摸一边抽搐一进一出视频| 好看av亚洲va欧美ⅴa在| 午夜日韩欧美国产| 中文字幕高清在线视频| 亚洲在线自拍视频| 91老司机精品| 天堂√8在线中文| 日韩欧美一区二区三区在线观看| 久久精品国产亚洲av香蕉五月| 男女午夜视频在线观看| 18禁观看日本| 色哟哟哟哟哟哟| 一级a爱视频在线免费观看| 少妇熟女aⅴ在线视频| 咕卡用的链子| 亚洲av成人一区二区三| 欧美最黄视频在线播放免费| 精品国产乱子伦一区二区三区| 精品一区二区三区av网在线观看| 天堂影院成人在线观看| 午夜成年电影在线免费观看| 成人亚洲精品av一区二区| 操美女的视频在线观看| 久9热在线精品视频| 后天国语完整版免费观看| 一边摸一边抽搐一进一小说| 国产精品久久久久久人妻精品电影| 男女之事视频高清在线观看| 国产精品av久久久久免费| 黄色视频,在线免费观看| 99热只有精品国产| 人妻丰满熟妇av一区二区三区| 欧美黄色淫秽网站| 在线av久久热| 日韩有码中文字幕| 中文字幕人妻熟女乱码| 久久精品91无色码中文字幕| 十八禁网站免费在线| 亚洲片人在线观看| 村上凉子中文字幕在线| 夜夜看夜夜爽夜夜摸| 日本黄色视频三级网站网址| 亚洲中文字幕日韩| 色哟哟哟哟哟哟| 在线观看www视频免费| 777久久人妻少妇嫩草av网站| 亚洲最大成人中文| 中文字幕高清在线视频| 一本综合久久免费| 欧美国产精品va在线观看不卡| 天天一区二区日本电影三级 | 亚洲熟妇熟女久久| 女生性感内裤真人,穿戴方法视频| 色综合站精品国产| or卡值多少钱| 日韩视频一区二区在线观看| 亚洲人成77777在线视频| 亚洲久久久国产精品| 久久精品国产清高在天天线| 男女床上黄色一级片免费看| 长腿黑丝高跟| 黄片小视频在线播放| 巨乳人妻的诱惑在线观看| 久久久国产欧美日韩av| 亚洲精品在线观看二区| av在线播放免费不卡| 久热爱精品视频在线9| av天堂在线播放| 中文字幕人妻熟女乱码| 一边摸一边抽搐一进一小说| 男女下面插进去视频免费观看| 俄罗斯特黄特色一大片| 男女之事视频高清在线观看| 国产精品久久视频播放| 亚洲精华国产精华精| 亚洲一码二码三码区别大吗| 9191精品国产免费久久| 高清在线国产一区| 最新美女视频免费是黄的| 看免费av毛片| 91在线观看av| 精品国产国语对白av| 丝袜美足系列| 亚洲一区二区三区不卡视频| 少妇熟女aⅴ在线视频| 欧美 亚洲 国产 日韩一| 精品不卡国产一区二区三区| 国产午夜福利久久久久久| 老司机靠b影院| 国语自产精品视频在线第100页| 欧美激情久久久久久爽电影 | av免费在线观看网站| 国产亚洲精品综合一区在线观看 | 亚洲专区国产一区二区| 高清在线国产一区| 狂野欧美激情性xxxx| av片东京热男人的天堂| 欧美人与性动交α欧美精品济南到| 最近最新中文字幕大全免费视频| 免费观看人在逋| 黑丝袜美女国产一区| 久久久水蜜桃国产精品网| 亚洲视频免费观看视频| 91麻豆精品激情在线观看国产| 日韩欧美一区视频在线观看| 亚洲自拍偷在线| 亚洲美女黄片视频| 免费看a级黄色片| 日韩精品免费视频一区二区三区| 国产高清激情床上av| 亚洲自拍偷在线| 波多野结衣高清无吗| 亚洲全国av大片| 日韩免费av在线播放| www日本在线高清视频| 我的亚洲天堂| 可以免费在线观看a视频的电影网站| 好男人在线观看高清免费视频 | 19禁男女啪啪无遮挡网站| 亚洲熟妇中文字幕五十中出| 国产一区二区三区综合在线观看| 黑丝袜美女国产一区| 国产精品久久久久久人妻精品电影| 亚洲aⅴ乱码一区二区在线播放 | 午夜久久久久精精品| 成人亚洲精品av一区二区| 午夜久久久在线观看| 成熟少妇高潮喷水视频| 成年女人毛片免费观看观看9| 无人区码免费观看不卡| 十分钟在线观看高清视频www| 黄色片一级片一级黄色片| 美女 人体艺术 gogo| 老司机午夜福利在线观看视频| 校园春色视频在线观看| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲av一区麻豆| 可以免费在线观看a视频的电影网站| 免费在线观看亚洲国产| 国产欧美日韩综合在线一区二区| 精品熟女少妇八av免费久了| 国产麻豆成人av免费视频| 黄色片一级片一级黄色片| 女人精品久久久久毛片| 国产单亲对白刺激| 亚洲欧美精品综合久久99| 最新美女视频免费是黄的| 亚洲情色 制服丝袜| 男男h啪啪无遮挡| 国内毛片毛片毛片毛片毛片| 亚洲精品粉嫩美女一区| 19禁男女啪啪无遮挡网站| 99热只有精品国产| 1024香蕉在线观看| av欧美777| 色在线成人网| 91成人精品电影| 亚洲专区中文字幕在线| 97碰自拍视频| 人人妻,人人澡人人爽秒播| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品在线观看二区| 999精品在线视频| 亚洲伊人色综图| 天天躁夜夜躁狠狠躁躁| 天天添夜夜摸| 99re在线观看精品视频| 国产男靠女视频免费网站| 大型黄色视频在线免费观看| 国产精品综合久久久久久久免费 | 国产xxxxx性猛交| 啦啦啦 在线观看视频| 亚洲欧美一区二区三区黑人| 男女之事视频高清在线观看| 亚洲九九香蕉| 国产av精品麻豆| 91成年电影在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久天躁狠狠躁夜夜2o2o| 欧美成狂野欧美在线观看| 国产麻豆成人av免费视频| 亚洲免费av在线视频| 91麻豆av在线| 国产精品永久免费网站| 国产片内射在线| 麻豆一二三区av精品| 91精品国产国语对白视频| 免费不卡黄色视频| 免费在线观看亚洲国产| 久久午夜综合久久蜜桃| 成人精品一区二区免费| 中文字幕av电影在线播放| 欧美乱色亚洲激情| 亚洲欧美激情在线| 国产精品综合久久久久久久免费 | 制服丝袜大香蕉在线| 精品免费久久久久久久清纯| 男女之事视频高清在线观看| 免费看十八禁软件| 亚洲五月婷婷丁香| 精品无人区乱码1区二区| 国产激情久久老熟女| 最近最新中文字幕大全电影3 | av超薄肉色丝袜交足视频| 久久午夜亚洲精品久久| av视频在线观看入口| 高清毛片免费观看视频网站| 国产成人欧美在线观看| 日韩大尺度精品在线看网址 | 99精品欧美一区二区三区四区| 男女之事视频高清在线观看| 久久香蕉激情| 美女扒开内裤让男人捅视频| 9191精品国产免费久久| 久久香蕉精品热| 大型黄色视频在线免费观看| 免费av毛片视频| 黑丝袜美女国产一区| 久久午夜亚洲精品久久| 色av中文字幕| 麻豆一二三区av精品| 老熟妇乱子伦视频在线观看| 91字幕亚洲| 在线免费观看的www视频| 中文字幕另类日韩欧美亚洲嫩草| 99riav亚洲国产免费| 自拍欧美九色日韩亚洲蝌蚪91| 免费一级毛片在线播放高清视频 | 国产亚洲精品久久久久5区| 久久久国产成人免费| 中文字幕色久视频| 91大片在线观看| www日本在线高清视频| 97碰自拍视频| 国产91精品成人一区二区三区| 男女之事视频高清在线观看| 禁无遮挡网站| 三级毛片av免费| www.熟女人妻精品国产| 久久久久久久久中文| 国产亚洲精品久久久久久毛片| 日韩精品青青久久久久久| 久热这里只有精品99| 一二三四社区在线视频社区8| 亚洲人成77777在线视频| 91国产中文字幕| 国产97色在线日韩免费| www.熟女人妻精品国产| 国产精品亚洲美女久久久| 亚洲精品一卡2卡三卡4卡5卡| 18禁观看日本| 精品国产一区二区三区四区第35| 亚洲久久久国产精品| 免费搜索国产男女视频| 国产色视频综合| 久久这里只有精品19| 国产精品 欧美亚洲| 色播在线永久视频| 亚洲三区欧美一区| 亚洲自拍偷在线| 久久久水蜜桃国产精品网| av视频在线观看入口| 日韩一卡2卡3卡4卡2021年| 高潮久久久久久久久久久不卡| 老司机在亚洲福利影院| 悠悠久久av| 国产av一区二区精品久久| 一边摸一边做爽爽视频免费| 一级毛片精品| 搡老熟女国产l中国老女人| 女性被躁到高潮视频| 日韩欧美国产一区二区入口| 中文字幕人妻丝袜一区二区| 午夜亚洲福利在线播放| 亚洲人成电影免费在线| 亚洲精品国产色婷婷电影| 香蕉久久夜色| 国产午夜福利久久久久久| 国产真人三级小视频在线观看| 亚洲熟妇熟女久久| 嫩草影院精品99| 黄片大片在线免费观看| 99久久综合精品五月天人人| 欧美最黄视频在线播放免费| 亚洲中文字幕日韩| 97人妻精品一区二区三区麻豆 | 国产熟女午夜一区二区三区| 母亲3免费完整高清在线观看| 99国产精品一区二区蜜桃av| 欧美日韩福利视频一区二区| 级片在线观看| 老司机深夜福利视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 欧美日本中文国产一区发布| 电影成人av| 91麻豆av在线| xxx96com| 搡老熟女国产l中国老女人| 老汉色∧v一级毛片| 国产精品久久久av美女十八| 国产精品爽爽va在线观看网站 | 可以在线观看毛片的网站| 国产高清videossex| 日日摸夜夜添夜夜添小说| 久久婷婷成人综合色麻豆| 午夜福利欧美成人| 91精品三级在线观看| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色 | 丝袜美足系列| 欧美绝顶高潮抽搐喷水| 国产99白浆流出| 男女下面进入的视频免费午夜 | 村上凉子中文字幕在线| www.999成人在线观看| av超薄肉色丝袜交足视频| 99久久综合精品五月天人人| 香蕉久久夜色| 国产欧美日韩一区二区三区在线| av视频免费观看在线观看| 色播亚洲综合网| 夜夜夜夜夜久久久久| 日日摸夜夜添夜夜添小说| 99久久久亚洲精品蜜臀av| 精品国产美女av久久久久小说| 色播亚洲综合网| 国产亚洲欧美98| 无限看片的www在线观看| 国产精品综合久久久久久久免费 | a在线观看视频网站| 日韩欧美在线二视频| 欧美不卡视频在线免费观看 | 中文字幕色久视频| 韩国av一区二区三区四区| 人人妻人人澡欧美一区二区 | 免费在线观看日本一区| 亚洲中文字幕日韩| 美女高潮喷水抽搐中文字幕| 久久精品成人免费网站| 国产人伦9x9x在线观看| netflix在线观看网站| 成人av一区二区三区在线看| 国产精品久久久人人做人人爽| 叶爱在线成人免费视频播放| 久99久视频精品免费| 黄色a级毛片大全视频| 麻豆成人av在线观看| 日本欧美视频一区| videosex国产| 99精品久久久久人妻精品| 制服人妻中文乱码| 国语自产精品视频在线第100页| 午夜福利18| 51午夜福利影视在线观看| 欧美另类亚洲清纯唯美| 亚洲午夜精品一区,二区,三区| 村上凉子中文字幕在线| 亚洲七黄色美女视频| 亚洲一区高清亚洲精品| 中文字幕精品免费在线观看视频| 免费少妇av软件| 满18在线观看网站| 91麻豆av在线| 老熟妇乱子伦视频在线观看| 日韩欧美一区二区三区在线观看| 日本 欧美在线| 日韩大尺度精品在线看网址 | 在线观看免费日韩欧美大片| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费| 午夜福利,免费看| 国产熟女午夜一区二区三区| 国产激情欧美一区二区| 18禁观看日本| 久久久精品欧美日韩精品| 亚洲成国产人片在线观看| 女性被躁到高潮视频| 欧美在线黄色| 黄片大片在线免费观看| 看片在线看免费视频| 999久久久国产精品视频| 日韩欧美三级三区| 99在线视频只有这里精品首页| 啪啪无遮挡十八禁网站| 在线永久观看黄色视频| 在线观看午夜福利视频| av天堂久久9| 亚洲欧美日韩另类电影网站| 男女下面进入的视频免费午夜 | 999久久久精品免费观看国产| 午夜福利一区二区在线看| 香蕉丝袜av| 91麻豆av在线| 国产乱人伦免费视频| 亚洲专区中文字幕在线| 欧美黄色片欧美黄色片| 成人手机av| 色综合欧美亚洲国产小说| 97人妻天天添夜夜摸| 啦啦啦免费观看视频1| 免费观看人在逋| 精品熟女少妇八av免费久了| 精品国产乱子伦一区二区三区| 亚洲欧美激情在线| 在线永久观看黄色视频| 1024视频免费在线观看| 国产高清videossex| 不卡一级毛片| 亚洲精品美女久久av网站| 国产成人一区二区三区免费视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美色欧美亚洲另类二区 | 精品国产亚洲在线| 国产精品亚洲av一区麻豆| 99久久综合精品五月天人人| 亚洲人成电影免费在线| 制服丝袜大香蕉在线| 91麻豆精品激情在线观看国产| 在线观看66精品国产| 1024香蕉在线观看| or卡值多少钱| 日韩欧美在线二视频| 夜夜看夜夜爽夜夜摸| 久久精品aⅴ一区二区三区四区| 国产精品爽爽va在线观看网站 | 极品教师在线免费播放| 一卡2卡三卡四卡精品乱码亚洲| 身体一侧抽搐| 久久久国产精品麻豆| 国产片内射在线| 日韩欧美国产一区二区入口| 免费高清在线观看日韩| 女人被狂操c到高潮| 色哟哟哟哟哟哟| 国产精品av久久久久免费| www.熟女人妻精品国产| 99国产极品粉嫩在线观看| 满18在线观看网站| 午夜久久久久精精品| 欧美+亚洲+日韩+国产| 免费搜索国产男女视频| 精品久久久久久久人妻蜜臀av | 国产av精品麻豆| 午夜久久久久精精品| 亚洲成人久久性| 久久人妻熟女aⅴ| 91老司机精品| 免费少妇av软件| 一级作爱视频免费观看| 不卡av一区二区三区| 亚洲熟妇熟女久久| 午夜福利成人在线免费观看| 成人国产一区最新在线观看| a在线观看视频网站| 搡老熟女国产l中国老女人| 男女做爰动态图高潮gif福利片 | 欧美一级a爱片免费观看看 | 日韩一卡2卡3卡4卡2021年| 亚洲最大成人中文| 国产一区二区在线av高清观看| 欧美色欧美亚洲另类二区 | 成在线人永久免费视频| 精品日产1卡2卡| 午夜激情av网站| 亚洲成人国产一区在线观看| 香蕉国产在线看| 国产欧美日韩精品亚洲av| 日韩有码中文字幕| 九色亚洲精品在线播放| av欧美777| 一级片免费观看大全| 狠狠狠狠99中文字幕| 狂野欧美激情性xxxx| 好看av亚洲va欧美ⅴa在| 美女大奶头视频| a在线观看视频网站| 美国免费a级毛片| 男女做爰动态图高潮gif福利片 | 国产精品二区激情视频| 日韩av在线大香蕉| 国产国语露脸激情在线看| 欧美乱码精品一区二区三区| 美女大奶头视频| 一区二区三区精品91| 亚洲欧美日韩高清在线视频| 三级毛片av免费| 欧美日韩一级在线毛片| 午夜免费鲁丝| 黄网站色视频无遮挡免费观看| 亚洲七黄色美女视频| 成人18禁在线播放| 女同久久另类99精品国产91| 男女之事视频高清在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久青草综合色| 午夜两性在线视频| 深夜精品福利| 免费不卡黄色视频| 一区二区三区国产精品乱码| 亚洲全国av大片| 亚洲男人天堂网一区| 99精品久久久久人妻精品| 91在线观看av| 精品国产一区二区三区四区第35| 国产精品久久久久久人妻精品电影| 国产成人免费无遮挡视频| 男女午夜视频在线观看| 少妇裸体淫交视频免费看高清 | 国产一区在线观看成人免费| 国产高清有码在线观看视频 | 亚洲欧美激情综合另类| 窝窝影院91人妻| 欧美日韩黄片免| 国产成人精品在线电影| 国产亚洲欧美精品永久| 老司机靠b影院| 日本精品一区二区三区蜜桃| 欧美日韩福利视频一区二区| 女性被躁到高潮视频| 一区二区三区激情视频| 免费观看人在逋| 国产精品 欧美亚洲| 亚洲专区国产一区二区| 99国产精品免费福利视频| 狠狠狠狠99中文字幕| 精品国产乱子伦一区二区三区| 嫩草影视91久久| 免费高清在线观看日韩| 我的亚洲天堂| 亚洲精品一区av在线观看| а√天堂www在线а√下载| 黄色a级毛片大全视频| 国产精品久久久久久精品电影 | 999久久久精品免费观看国产| 母亲3免费完整高清在线观看| 大码成人一级视频| 亚洲va日本ⅴa欧美va伊人久久| 在线永久观看黄色视频| 嫁个100分男人电影在线观看| 一边摸一边抽搐一进一出视频| 国产精品久久久av美女十八| 欧美日韩乱码在线| 久久精品国产99精品国产亚洲性色 | 亚洲中文字幕日韩| 亚洲av日韩精品久久久久久密| 亚洲欧美激情综合另类| 真人做人爱边吃奶动态| 亚洲第一av免费看| 性色av乱码一区二区三区2| 女人被狂操c到高潮| 日日爽夜夜爽网站| 成人特级黄色片久久久久久久| 国产午夜精品久久久久久| 一区二区三区精品91| 淫妇啪啪啪对白视频| ponron亚洲| 一区福利在线观看| 啪啪无遮挡十八禁网站| 欧美一级a爱片免费观看看 | 黄色 视频免费看| 亚洲av第一区精品v没综合|