• <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ù)值計算
    欧美成人一区二区免费高清观看| 有码 亚洲区| 精品熟女少妇av免费看| 亚洲人成网站高清观看| 丝袜美腿在线中文| 亚洲一级一片aⅴ在线观看| 婷婷精品国产亚洲av| 国产精品久久电影中文字幕| 国产精品久久视频播放| 久久中文看片网| 一进一出抽搐gif免费好疼| 免费人成在线观看视频色| 三级毛片av免费| 亚洲av一区综合| 国内揄拍国产精品人妻在线| 伦理电影大哥的女人| 日韩欧美国产在线观看| 好男人在线观看高清免费视频| 日韩欧美一区二区三区在线观看| 免费人成在线观看视频色| 久久精品夜色国产| 国产真实伦视频高清在线观看| 亚洲丝袜综合中文字幕| 国产成人a∨麻豆精品| 亚洲成人中文字幕在线播放| 国产白丝娇喘喷水9色精品| 级片在线观看| 麻豆成人午夜福利视频| 国产视频首页在线观看| 97人妻精品一区二区三区麻豆| 看十八女毛片水多多多| 又黄又爽又刺激的免费视频.| 一进一出抽搐动态| 青青草视频在线视频观看| 久久精品国产自在天天线| 人体艺术视频欧美日本| 精品国内亚洲2022精品成人| 人妻系列 视频| 亚洲国产精品成人综合色| 免费电影在线观看免费观看| 亚洲激情五月婷婷啪啪| 日韩成人伦理影院| 精品久久久久久成人av| 午夜福利视频1000在线观看| 色综合站精品国产| www日本黄色视频网| 九九在线视频观看精品| 国产成人午夜福利电影在线观看| 国产精品蜜桃在线观看 | 嫩草影院精品99| 3wmmmm亚洲av在线观看| 欧美区成人在线视频| 国产激情偷乱视频一区二区| 国产精品久久久久久久久免| 午夜老司机福利剧场| 亚洲无线在线观看| 久久精品久久久久久噜噜老黄 | 久久午夜福利片| 又粗又硬又长又爽又黄的视频 | av又黄又爽大尺度在线免费看 | 亚洲三级黄色毛片| 亚洲一区高清亚洲精品| 最后的刺客免费高清国语| 九九热线精品视视频播放| 少妇被粗大猛烈的视频| av在线蜜桃| 黄色配什么色好看| 成人美女网站在线观看视频| 国产69精品久久久久777片| 亚洲欧美清纯卡通| 舔av片在线| 亚洲综合色惰| 国产三级中文精品| 在线a可以看的网站| 少妇的逼水好多| 色尼玛亚洲综合影院| 国产精品美女特级片免费视频播放器| 日本免费a在线| 日韩三级伦理在线观看| 久久人人精品亚洲av| 欧美三级亚洲精品| 麻豆成人av视频| 黄色欧美视频在线观看| 老司机福利观看| 国产精品美女特级片免费视频播放器| 欧美在线一区亚洲| 色哟哟·www| 男女下面进入的视频免费午夜| 国产女主播在线喷水免费视频网站 | 国产精品99久久久久久久久| 高清在线视频一区二区三区 | 亚洲av男天堂| 国产亚洲5aaaaa淫片| 国产老妇女一区| 国产精品不卡视频一区二区| 天堂影院成人在线观看| 男人狂女人下面高潮的视频| 美女cb高潮喷水在线观看| 搡老妇女老女人老熟妇| 悠悠久久av| 91精品一卡2卡3卡4卡| 久久人人爽人人片av| 18禁裸乳无遮挡免费网站照片| 免费看美女性在线毛片视频| 久久久久久久久久黄片| 深夜a级毛片| 国产高潮美女av| 在线观看免费视频日本深夜| 国产午夜精品论理片| 尤物成人国产欧美一区二区三区| 国产日本99.免费观看| 亚洲四区av| 一个人观看的视频www高清免费观看| 国产一区二区三区在线臀色熟女| 日日撸夜夜添| 哪个播放器可以免费观看大片| 欧美精品一区二区大全| 成人亚洲欧美一区二区av| 欧美人与善性xxx| 小蜜桃在线观看免费完整版高清| 国产又黄又爽又无遮挡在线| 免费看日本二区| 一级黄片播放器| 亚洲国产精品成人久久小说 | 国产免费男女视频| 美女被艹到高潮喷水动态| 老女人水多毛片| 国产91av在线免费观看| 小蜜桃在线观看免费完整版高清| 中文字幕人妻熟人妻熟丝袜美| 女的被弄到高潮叫床怎么办| 人妻系列 视频| 乱系列少妇在线播放| 春色校园在线视频观看| 欧美成人a在线观看| 综合色丁香网| 一级毛片电影观看 | 一级黄色大片毛片| 看免费成人av毛片| 中国美白少妇内射xxxbb| 天堂影院成人在线观看| 久久精品久久久久久久性| 狂野欧美白嫩少妇大欣赏| 日韩亚洲欧美综合| 亚洲婷婷狠狠爱综合网| 两个人视频免费观看高清| 嫩草影院精品99| 嫩草影院新地址| a级毛片a级免费在线| 黑人高潮一二区| 日韩制服骚丝袜av| 国内揄拍国产精品人妻在线| 亚洲成人中文字幕在线播放| 黄片无遮挡物在线观看| 波多野结衣高清无吗| 不卡一级毛片| 天堂中文最新版在线下载 | 日本在线视频免费播放| 欧美成人精品欧美一级黄| 99热这里只有是精品在线观看| 国产成人aa在线观看| 久久久久九九精品影院| 在线观看午夜福利视频| 日本一本二区三区精品| 人人妻人人看人人澡| 国产三级中文精品| 在线免费观看不下载黄p国产| 久久精品综合一区二区三区| 午夜福利在线观看免费完整高清在 | 国产av一区在线观看免费| 毛片一级片免费看久久久久| 综合色av麻豆| 免费观看a级毛片全部| 久久久久久久久中文| 欧美bdsm另类| av在线蜜桃| 岛国毛片在线播放| 国产精品乱码一区二三区的特点| 国产精品久久电影中文字幕| 午夜a级毛片| 97人妻精品一区二区三区麻豆| 亚洲人成网站高清观看| 国产精品精品国产色婷婷| 变态另类成人亚洲欧美熟女| 日本免费a在线| 精品国产三级普通话版| 悠悠久久av| 伦理电影大哥的女人| 久久这里有精品视频免费| 又爽又黄a免费视频| 国产av一区在线观看免费| 国产伦理片在线播放av一区 | 午夜a级毛片| 激情 狠狠 欧美| videossex国产| 桃色一区二区三区在线观看| 免费观看精品视频网站| av在线亚洲专区| 白带黄色成豆腐渣| 国产高潮美女av| 久久中文看片网| 国产精品,欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 少妇被粗大猛烈的视频| 成人毛片a级毛片在线播放| 黄片无遮挡物在线观看| 99久久中文字幕三级久久日本| 美女黄网站色视频| 亚洲一级一片aⅴ在线观看| 身体一侧抽搐| 久久午夜福利片| 午夜精品在线福利| 欧美丝袜亚洲另类| 国产麻豆成人av免费视频| 国产日本99.免费观看| 可以在线观看的亚洲视频| 一卡2卡三卡四卡精品乱码亚洲| 成人一区二区视频在线观看| 人妻制服诱惑在线中文字幕| 国产极品天堂在线| 国产真实伦视频高清在线观看| 国内揄拍国产精品人妻在线| 亚洲av一区综合| av视频在线观看入口| 国产精品美女特级片免费视频播放器| 国产精华一区二区三区| 亚洲精品456在线播放app| 男人舔女人下体高潮全视频| 久久久精品94久久精品| 亚洲国产精品sss在线观看| 日韩欧美三级三区| 成人无遮挡网站| 亚洲国产精品国产精品| 免费av观看视频| 午夜福利视频1000在线观看| 日本在线视频免费播放| 国产精品99久久久久久久久| 三级国产精品欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人福利小说| 97在线视频观看| 狠狠狠狠99中文字幕| 亚洲精品影视一区二区三区av| 国产成人影院久久av| 人妻系列 视频| 免费黄网站久久成人精品| 亚洲18禁久久av| 91狼人影院| 人人妻人人看人人澡| 精品久久久久久久人妻蜜臀av| 亚洲,欧美,日韩| 日韩精品青青久久久久久| 亚洲精品色激情综合| 18禁裸乳无遮挡免费网站照片| 免费看日本二区| 亚洲成a人片在线一区二区| 熟女人妻精品中文字幕| 国产黄色视频一区二区在线观看 | 亚洲精品久久久久久婷婷小说 | 国产午夜精品一二区理论片| 国产黄色视频一区二区在线观看 | 久久九九热精品免费| 1024手机看黄色片| 男插女下体视频免费在线播放| 欧美极品一区二区三区四区| 久久久久免费精品人妻一区二区| 中文字幕av成人在线电影| 99久久精品热视频| 国产精品,欧美在线| 亚洲18禁久久av| 麻豆成人午夜福利视频| 99久久精品一区二区三区| 成人欧美大片| 国产毛片a区久久久久| 男人舔奶头视频| 久久久精品欧美日韩精品| 国产精品久久久久久久久免| 久久久久久久久中文| www.av在线官网国产| 成年女人看的毛片在线观看| 成人三级黄色视频| 麻豆国产97在线/欧美| 给我免费播放毛片高清在线观看| 亚洲在线观看片| av在线蜜桃| 国产精品三级大全| 久久国内精品自在自线图片| 成熟少妇高潮喷水视频| 黄色欧美视频在线观看| 亚洲性久久影院| 亚洲乱码一区二区免费版| 97在线视频观看| 黄色视频,在线免费观看| 91狼人影院| 中文亚洲av片在线观看爽| 久久精品国产亚洲av涩爱 | 日韩精品有码人妻一区| 日日摸夜夜添夜夜爱| 日韩 亚洲 欧美在线| 精品久久久久久久久久免费视频| 人人妻人人澡人人爽人人夜夜 | 观看免费一级毛片| 白带黄色成豆腐渣| 少妇的逼好多水| 夜夜夜夜夜久久久久| 亚洲中文字幕日韩| 深爱激情五月婷婷| 特大巨黑吊av在线直播| 精品不卡国产一区二区三区| 少妇的逼好多水| 在线免费观看不下载黄p国产| 成人永久免费在线观看视频| av免费在线看不卡| 可以在线观看的亚洲视频| 国产精品国产高清国产av| 国产精品,欧美在线| 国产日本99.免费观看| 波野结衣二区三区在线| 国产精品国产三级国产av玫瑰| 亚洲人与动物交配视频| 日本熟妇午夜| 国产伦精品一区二区三区视频9| 插阴视频在线观看视频| 99国产精品一区二区蜜桃av| 综合色av麻豆| 亚洲美女搞黄在线观看| 三级经典国产精品| 丰满乱子伦码专区| 国产色婷婷99| 亚洲av不卡在线观看| 99热6这里只有精品| 国产av一区在线观看免费| 波多野结衣巨乳人妻| 男女啪啪激烈高潮av片| 一夜夜www| 成人美女网站在线观看视频| 久久99蜜桃精品久久| 欧美xxxx黑人xx丫x性爽| 老师上课跳d突然被开到最大视频| 午夜福利高清视频| 91在线精品国自产拍蜜月| 免费观看精品视频网站| 美女大奶头视频| 国产免费一级a男人的天堂| 看十八女毛片水多多多| 男女下面进入的视频免费午夜| 国产亚洲精品久久久com| 麻豆成人av视频| 长腿黑丝高跟| 精品久久久久久久久久久久久| 国产亚洲5aaaaa淫片| 国产亚洲av嫩草精品影院| 伦精品一区二区三区| 日韩三级伦理在线观看| 婷婷色av中文字幕| av福利片在线观看| 嫩草影院精品99| 老师上课跳d突然被开到最大视频| 国产亚洲5aaaaa淫片| 亚洲aⅴ乱码一区二区在线播放| 国产成人a区在线观看| 亚洲中文字幕一区二区三区有码在线看| 99热网站在线观看| 午夜激情福利司机影院| 熟妇人妻久久中文字幕3abv| 亚洲精品成人久久久久久| 欧美色视频一区免费| 亚洲最大成人中文| 久久久久久久久久黄片| 最近的中文字幕免费完整| 国产高清不卡午夜福利| 欧美一区二区亚洲| 内射极品少妇av片p| 久久久久久国产a免费观看| 久久人妻av系列| 国产伦精品一区二区三区视频9| 日韩一区二区三区影片| ponron亚洲| 禁无遮挡网站| 久久鲁丝午夜福利片| 国产成人91sexporn| 99riav亚洲国产免费| 一级毛片电影观看 | 91久久精品国产一区二区三区| 成人鲁丝片一二三区免费| 亚洲色图av天堂| www.色视频.com| 国产精品久久久久久av不卡| 国产淫片久久久久久久久| a级毛色黄片| 久久久久免费精品人妻一区二区| 少妇被粗大猛烈的视频| 99在线视频只有这里精品首页| 久久人妻av系列| 天堂网av新在线| 12—13女人毛片做爰片一| 亚洲七黄色美女视频| 丝袜喷水一区| 22中文网久久字幕| 岛国在线免费视频观看| 一本精品99久久精品77| 3wmmmm亚洲av在线观看| 亚洲成人久久爱视频| 免费看光身美女| 国产成人福利小说| 天堂影院成人在线观看| 2021天堂中文幕一二区在线观| 男人和女人高潮做爰伦理| 日韩,欧美,国产一区二区三区 | 人人妻人人澡人人爽人人夜夜 | 欧美在线一区亚洲| 一本精品99久久精品77| 亚洲精品乱码久久久久久按摩| 狂野欧美激情性xxxx在线观看| 嘟嘟电影网在线观看| 夫妻性生交免费视频一级片| 国产视频首页在线观看| 亚洲av电影不卡..在线观看| 成人毛片a级毛片在线播放| 能在线免费看毛片的网站| 亚洲欧美成人精品一区二区| 国产高清视频在线观看网站| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 亚洲精品影视一区二区三区av| 天堂av国产一区二区熟女人妻| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 中文字幕人妻熟人妻熟丝袜美| 国产极品天堂在线| 18禁在线无遮挡免费观看视频| 尾随美女入室| 黄片wwwwww| 日本一二三区视频观看| 国产精品电影一区二区三区| 久久精品人妻少妇| 国产高清有码在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利在线观看免费完整高清在 | 97在线视频观看| 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 美女cb高潮喷水在线观看| 日韩国内少妇激情av| 成人特级av手机在线观看| 午夜久久久久精精品| 国产精品综合久久久久久久免费| 欧美+日韩+精品| 我的女老师完整版在线观看| 美女cb高潮喷水在线观看| 99视频精品全部免费 在线| 欧美成人一区二区免费高清观看| 久久鲁丝午夜福利片| 校园人妻丝袜中文字幕| 97超碰精品成人国产| 2022亚洲国产成人精品| 男人的好看免费观看在线视频| 日韩成人伦理影院| 亚洲精品色激情综合| 精品99又大又爽又粗少妇毛片| 97人妻精品一区二区三区麻豆| 久久精品国产亚洲av天美| 亚洲一级一片aⅴ在线观看| 天堂av国产一区二区熟女人妻| 午夜亚洲福利在线播放| 变态另类成人亚洲欧美熟女| 国产精品美女特级片免费视频播放器| 亚洲国产精品sss在线观看| 在线播放国产精品三级| 成年女人永久免费观看视频| 国产精品99久久久久久久久| 性色avwww在线观看| 熟女电影av网| 男人的好看免费观看在线视频| 在线天堂最新版资源| 给我免费播放毛片高清在线观看| av在线观看视频网站免费| 欧美成人免费av一区二区三区| 99热网站在线观看| 男女那种视频在线观看| 99热只有精品国产| 午夜精品一区二区三区免费看| 亚洲无线在线观看| 麻豆成人av视频| 99九九线精品视频在线观看视频| 在线播放无遮挡| 男女视频在线观看网站免费| 国产免费一级a男人的天堂| 夜夜夜夜夜久久久久| 国产伦精品一区二区三区视频9| 中文字幕久久专区| 99九九线精品视频在线观看视频| 一级黄色大片毛片| 亚洲av免费高清在线观看| av在线观看视频网站免费| 久久精品国产亚洲av天美| 天堂√8在线中文| 欧美日本视频| 在线观看午夜福利视频| 波多野结衣高清无吗| 国产精品久久电影中文字幕| 国产成人福利小说| 欧美激情久久久久久爽电影| 给我免费播放毛片高清在线观看| 久久人人精品亚洲av| 精品久久久久久久久久久久久| 国产精品日韩av在线免费观看| 国产精品乱码一区二三区的特点| 日本欧美国产在线视频| 自拍偷自拍亚洲精品老妇| 国产av一区在线观看免费| 一级毛片电影观看 | 久久久久国产网址| 1024手机看黄色片| 精华霜和精华液先用哪个| 亚洲av成人精品一区久久| 偷拍熟女少妇极品色| 日本成人三级电影网站| 非洲黑人性xxxx精品又粗又长| 91av网一区二区| 亚洲av成人精品一区久久| 久久综合国产亚洲精品| 国产熟女欧美一区二区| 日本黄色片子视频| 久久久久久伊人网av| 亚洲不卡免费看| 国产精品久久久久久精品电影| 中文字幕av在线有码专区| 校园人妻丝袜中文字幕| 身体一侧抽搐| 女的被弄到高潮叫床怎么办| 国产男人的电影天堂91| 日本熟妇午夜| 成人av在线播放网站| 久久精品久久久久久噜噜老黄 | 国内少妇人妻偷人精品xxx网站| 国产一级毛片在线| 又粗又硬又长又爽又黄的视频 | 久久久精品欧美日韩精品| 欧美最新免费一区二区三区| 能在线免费观看的黄片| 亚洲欧洲国产日韩| 久久这里有精品视频免费| 少妇猛男粗大的猛烈进出视频 | 老司机福利观看| 天堂√8在线中文| 岛国毛片在线播放| 国产美女午夜福利| www.色视频.com| 欧美高清成人免费视频www| 天天躁日日操中文字幕| 欧美+亚洲+日韩+国产| 久久精品国产鲁丝片午夜精品| 91久久精品国产一区二区三区| av福利片在线观看| 国产视频内射| 免费电影在线观看免费观看| 日日摸夜夜添夜夜添av毛片| 亚洲va在线va天堂va国产| 黄片无遮挡物在线观看| 精品久久久噜噜| 日本av手机在线免费观看| 国内精品美女久久久久久| 成人亚洲精品av一区二区| 男女边吃奶边做爰视频| 变态另类成人亚洲欧美熟女| 国产精品一二三区在线看| 免费av毛片视频| 亚洲欧美成人精品一区二区| 免费观看的影片在线观看| 黄片无遮挡物在线观看| 国产色婷婷99| 久久韩国三级中文字幕| 国产精品综合久久久久久久免费| 国产久久久一区二区三区| 哪里可以看免费的av片| 一个人看的www免费观看视频| 精品无人区乱码1区二区| 美女黄网站色视频| 日本熟妇午夜| 最近2019中文字幕mv第一页| 少妇人妻一区二区三区视频| 老女人水多毛片| 亚洲av第一区精品v没综合| 寂寞人妻少妇视频99o| 校园春色视频在线观看| 丰满乱子伦码专区| 亚洲中文字幕日韩| 老司机福利观看| 丰满乱子伦码专区| 国产精品久久久久久久电影| 欧美成人精品欧美一级黄| 亚洲精品色激情综合| 99久久九九国产精品国产免费| 国产成人aa在线观看| av免费在线看不卡| 少妇熟女aⅴ在线视频| 乱系列少妇在线播放| 国产精品久久久久久精品电影| 国产精品99久久久久久久久| 国产人妻一区二区三区在| 久久久成人免费电影| 久久久a久久爽久久v久久| 狂野欧美激情性xxxx在线观看| 亚洲电影在线观看av| 久久精品国产99精品国产亚洲性色| 麻豆国产97在线/欧美| 99视频精品全部免费 在线| 极品教师在线视频| 亚洲无线在线观看| 天堂中文最新版在线下载 | 亚洲精品自拍成人| 99热这里只有是精品50| 国产高清有码在线观看视频| 久久欧美精品欧美久久欧美|