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

    基于中心組合設計的顆粒流平直節(jié)理模型宏-細觀參數(shù)相關性研究

    2021-04-16 06:37:26
    中南大學學報(自然科學版) 2021年3期
    關鍵詞:細觀節(jié)理單軸

    (中南大學土木工程學院,湖南長沙,410075)

    在礦產(chǎn)資源開采、邊坡、隧道、石油開采等巖體工程中,天然巖石材料存在的細觀缺陷(如裂紋或孔隙)極大地影響著巖石的宏觀力學性能。在外載作用下,裂紋或孔隙的演變發(fā)展形成宏觀缺陷會引起巖體工程災害隱患,因此,研究巖石細觀缺陷發(fā)展至宏觀破壞過程及其機理,對巖體工程的強度計算、安全評估和防災減災等具有重要理論和實際意義。由于現(xiàn)有試驗測試手段的局限性和巖石內(nèi)部缺陷的復雜性,人們難以在巖體工程中實時監(jiān)測到巖石細觀損傷積累至宏觀破壞的全過程,故多借助于數(shù)值模擬方法。目前,研究巖石裂紋擴展的數(shù)值模擬方法主要有離散元法、有限元法、邊界元法等,其中,基于非連續(xù)介質(zhì)力學理論的離散元法在模擬巖石裂隙擴展時具有獨特的優(yōu)勢,因為該方法既能避免選取復雜的本構(gòu)關系和裂尖奇異性,又能夠從細觀層面上模擬巖石宏觀力學特性并揭示其破壞機理。離散元顆粒流軟件PFC(particle flow code)在巖石類材料的數(shù)值模擬中應用較廣泛[1]。PFC 軟件是將巖石離散為剛性顆粒,通過定義顆粒接觸本構(gòu)關系來描述顆粒之間的相互作用,以此模擬巖石的宏觀力學特性。

    目前,PFC 顆粒黏結(jié)模型大多為接觸黏結(jié)(linear contact bond)、平行黏結(jié)(parallel bond)和平直節(jié)理模型(flat jointed model)。在接觸黏結(jié)模型中,黏結(jié)斷裂時顆粒接觸剛度仍然存在,這與巖石材料不相符[2];在平行黏結(jié)模型中,模擬得到的巖石抗拉強度與抗壓強度的比值偏大,與巖石材料的力學性能也不相符[3-4],通常采用簇平行和新的膠接模型來解決該問題[5-6]。在平直節(jié)理模型中[7],圓形顆粒被構(gòu)造成多邊形從而抑制了顆粒的旋轉(zhuǎn),使得該模型計算得到的單軸抗拉強度和抗壓強度之比減小,更符合巖石力學特性。采用顆粒流PFC2D軟件模擬巖石的宏觀力學特性時,首先需要確定顆粒黏結(jié)模型的細觀參數(shù)。目前,采用較多的是試錯法,如文獻[8-9]基于試錯法標定了平行黏結(jié)的細觀參數(shù);文獻[10]基于試錯法標定了平直節(jié)理接觸模型的細觀參數(shù),并與平行黏結(jié)模型的細觀參數(shù)進行對比,分析了平直節(jié)理模型的優(yōu)勢,但未建立宏-細觀參數(shù)之間的關系式。由于試錯法具有一定的盲目性,細觀參數(shù)較多且標定過程較繁瑣,從而導致試驗次數(shù)很多,因此,研究者尋求更有效的細觀參數(shù)定量確定方法,如:文獻[11]采用單因素和多因素分析方法,分析了平行黏結(jié)模型的細觀參數(shù)對宏觀參數(shù)的影響,建立了平行黏結(jié)模型的宏-細觀參數(shù)關系式;文獻[12]采用BP神經(jīng)網(wǎng)絡法,建立了簇平行黏結(jié)模型的宏-細觀參數(shù)的非線性模型,但BP 神經(jīng)網(wǎng)絡需要大量的隨機組合樣本才能具有較高的精度;文獻[13]采用球面對稱設計法建立了平行黏結(jié)宏-細觀參數(shù)的關系式,但考慮的宏觀參數(shù)有限;文獻[14-15]通過對平直節(jié)理模型細觀參數(shù)進行正交分析,得到了細觀參數(shù)對宏觀參數(shù)的影響規(guī)律,并利用回歸分析建立了宏-細觀參數(shù)關系式,但因試驗次數(shù)較少導致精度較低??傊壳瓣P于PFC2D平直節(jié)理接觸模型的細觀參數(shù)定量確定方法偏少。為提高精度、減少試驗次數(shù),需要尋求一種更合適的方法來建立PFC2D平直節(jié)理接觸模型的宏-細觀參數(shù)關系式。中心組合設計(又稱星點設計,central composite design)[16-17]是一種多因素五水平的試驗設計方法,相對于正交試驗和均勻化試驗,它在提高精度的同時不需要大量增加實驗次數(shù),特別適用于多因素、多水平試驗情形。

    本文作者基于PFC2D平直節(jié)理模型,通過控制變量法得到細觀參數(shù)的取值范圍和對宏觀參數(shù)的影響規(guī)律,采用中心組合設計和回歸分析,建立巖石PFC2D宏-細觀參數(shù)關系式以定量確定細觀參數(shù),模擬計算巖石在單軸和雙軸壓縮下完整巖石試件的應力-應變曲線和含多裂紋巖石試件的斷裂軌跡,并與試驗結(jié)果進行對比,以驗證所建立的巖石PFC2D宏-細觀參數(shù)關系式的有效性。

    1 PFC平直節(jié)理模型

    平直節(jié)理模型是針對硬巖提出的一種模型,如圖1所示。該模型將圓形顆粒構(gòu)造成多邊形,從而限制了顆粒的旋轉(zhuǎn),致使模型單軸抗壓強度和抗拉強度的比值增大,更加符合巖石類材料的力學性能。平直節(jié)理模型采用扁平面將2 個顆粒黏結(jié),黏結(jié)面由片段的單元組成。每個單元具有黏結(jié)和未黏結(jié)2種模式,當黏結(jié)界面上的荷載單元達到黏接強度時,黏結(jié)單元便會破壞成非黏結(jié)單元,其本構(gòu)關系如圖2所示。表1所示為平直節(jié)理模型的細觀參數(shù),包括顆粒細觀參數(shù)(6個)和平直節(jié)理參數(shù)(7個)。

    圖1 PFC2D平直節(jié)理模型[7]Fig.1 Flat joint model of PFC2D[7]

    圖2 PFC2D平直節(jié)理模型的本構(gòu)模型[7]Fig.2 Constitutive relation of PFC2D flat joint model[7]

    2 PFC細觀參數(shù)確定方法

    2.1 細觀參數(shù)取值范圍的確定

    為了模擬巖石類材料宏觀力學特性,首先需要確定PFC2D平直節(jié)理模型的細觀參數(shù)。由于平直節(jié)理模型的細觀參數(shù)較多(共13個),為了減少標定次數(shù),根據(jù)前人的研究成果[13],進行以下假設:

    表1 PFC2D平直節(jié)理模型的細觀參數(shù)[7]Table 1 Microscopic parameters of PFC2D flat joint model[7]

    1)顆粒接觸模量Eb等于平直節(jié)理接觸模量Ec;

    2)顆粒剛度比(kn/ks)b等于平直節(jié)理接觸剛度比kn/ks;

    3)最大、最小顆粒半徑之比Rmax/Rmin=1.66;

    4)顆粒密度ρ=2 300 kg/m3;

    5)平直節(jié)理半徑系數(shù)-λ=1;

    6)交界面段數(shù)N=4;

    7)平直節(jié)理抗拉強度小于抗剪強度。

    為了確定PFC2D細觀參數(shù),需要分析細觀參數(shù)對巖石材料宏觀力學參數(shù)(含彈性模量E、泊松比ν、抗壓強度σf、抗拉強度σt、黏聚力C2、摩擦角φ)的影響,以此確定細觀參數(shù)的范圍。針對剩余的7個細觀參數(shù)采用控制變量法分析,選取適當?shù)某跏技氂^參數(shù),即Ec=40 GPa,kn/ks=1,C1=30 MPa,σc=3 MPa,μb=0.1,μ=0.1,R=0.2 mm[18],在改變單一細觀參數(shù)的情況下開展巖石材料的單軸壓縮(E,ν,σf)、單軸拉伸(σt)和雙軸壓縮(C2,φ)數(shù)值模擬試驗,計算得到的單一細觀參數(shù)對宏觀力學參數(shù)的影響結(jié)果如圖3~9 所示,其中橫坐標為PFC2D細觀參數(shù),縱坐標為歸一化處理后的巖石宏觀參數(shù),即細觀參數(shù)改變后通過模擬計算得到的宏觀參數(shù)與由初始細觀參數(shù)模擬計算得到的宏觀參數(shù)之比。

    從圖3可知:接觸模量Ec對宏觀彈性模量E影響最大,對其他宏觀參數(shù)幾乎沒有影響(在圖中表現(xiàn)為其他參數(shù)曲線基本重合)。由圖4可知:剛度比kn/ks主要影響泊松比ν,隨著kn/ks的增加,宏觀彈性模量E、抗拉強度σt略減少(在圖中表現(xiàn)為E,σt和C2曲線基本重合),其他宏觀參數(shù)基本保持不變。由圖5可知:平直節(jié)理黏聚力C1主要影響σf和黏聚力C2,對泊松比ν和摩擦角φ的影響較小。由圖6可知:平直節(jié)理抗拉強度σc主要影響σt,隨著σc增加,σf略增加。由圖7和圖8可知:顆粒摩擦因素μb和平直節(jié)理內(nèi)摩擦因數(shù)μ主要影響摩擦角φ,且μb對φ的影響更顯著;隨著μ增加,σf略微增加,C2略微減少,對其他宏觀參數(shù)影響較小(在圖中表現(xiàn)為E、σt曲線基本重合)。由于μb和μ對宏觀參數(shù)的影響具有相同的規(guī)律,可取μ=μb[7]。由圖9可知,顆粒半徑R對宏觀參數(shù)影響的幅度較小,這與文獻[19]中的結(jié)論一致??紤]到計算效率和精度,本文選取最小半徑Rmin=0.3 mm。因此,待定的獨立細觀參數(shù)只有5個:Ec,kn/ks,C1,σc和μb。

    圖3 細觀參數(shù)Ec對宏觀參數(shù)的影響Fig.3 Effect of microscopic parameters Ec on macroscopic parameters

    圖4 細觀參數(shù)kn/ks對宏觀參數(shù)的影響Fig.4 Effect of microscopic parameters kn/ks on macroscopic parameters

    圖5 細觀參數(shù)C1對宏觀參數(shù)的影響Fig.5 Effect of microscopic parameters C1 on macroscopic parameters

    圖6 細觀參數(shù)σc對宏觀參數(shù)的影響Fig.6 Effect of microscopic parameters σc on macroscopic parameters

    圖7 細觀參數(shù)μb對宏觀參數(shù)的影響Fig.7 Effect of microscopic parameters μb on macroscopic parameters

    圖8 細觀參數(shù)μ對宏觀參數(shù)的影響Fig.8 Effect of microscopic parameters μ on macroscopic parameters

    圖9 細觀參數(shù)R對宏觀參數(shù)的影響Fig.9 Effect of microscopic parameters R on macroscopic parameters

    由上述分析可知,細觀參數(shù)Ec,kn/ks,σc和μb與宏觀參數(shù)關聯(lián)度最大的分別是E,v,σt和φ,可由對應的宏觀力學參數(shù)范圍來確定該細觀參數(shù)范圍。細觀參數(shù)C1與宏觀參數(shù)σf和C2關聯(lián)度均較大,考慮到巖石類材料的黏聚力C2差異較大,變化范圍難以確定,故C1的取值范圍由σf的取值范圍來確定。對于大多數(shù)巖石類材料,宏觀力學參數(shù)的大致范圍為[20-21]:E=30~120 GPa,v=0.10~0.35,σf=50~300 MPa,σt=3~25 MPa,φ=30°~65°。故由圖3~9 可大致確定細觀參數(shù)的取值范圍為:Ec=35~110 GPa,kn/ks=1.6~5.3,C1=25~150 MPa,σc=5.3~30 MPa,μb=0.1~0.7。

    2.2 基于中心組合設計樣本的細觀參數(shù)計算公式

    中心組合設計又稱星點設計,是一種多因素多水平的試驗設計方法,它是在二水平析因設計方法的基礎之上加上極值點和中心點組成的。可見,中心組合設計相對于正交、均勻設計具有較高的實驗精度,并且也不過多地增加實驗次數(shù),更適合于多水平多因素的試驗設計。為了定量確定PFC2D細觀參數(shù),采用中心組合設計法建立細觀參數(shù)與宏觀參數(shù)之間的關系式。考慮到待定的細觀參數(shù)有5個(Ec,kn/ks,C1,σc和μb),利用Design-Expert 軟件生成5 因素5 水平的中心組合設計代碼表(表2),X1,X2,X3,X4和X5分別代表Ec,kn/ks,C1,σc和μb。每個元素均有5 個水平代碼:-2.38,-1,0,+1 和+2.38,其中,代碼-2.38,0 和+2.38對應的值分別為該因素范圍的最小值、平均值、最大值,代碼-1和+1對應的值可通過“任意2個水平代碼之差的比值與該水平代碼對應值之差的比值相等”原則計算求得。以Ec為例,其范圍為35~110 GPa,則Ec的水平代碼-2.38,0和+2.38所對應的值分別為35,72.5和110。設其水平代碼-1和+1對應的值分別為x和y,根據(jù)上述原則,有

    經(jīng)求解得到x=56.74,y=88.26。同理,可求得其他因素水平代碼對應的值,見表3。將表3中各因素水平代碼對應的值代入表2,可得到中心組合設計試驗表,如表4所示。

    根據(jù)上述控制變量法分析(圖3~9)、中心組合設計試驗表(表4)可知:E受Ec和kn/ks的影響,v主要受kn/ks的影響,σf受C1,σc和μb的影響,σt受σc和kn/ks的影響,C2主要受C1和μb的影響,φ主要受C1和μb的影響。對6個數(shù)值模型宏觀力學性能(E,v,σf,σt,C2,φ)和5 個細觀參數(shù)(Ec,kn/ks,C1,σc,μb),利用SPSS軟件進行非線性回歸,得到模型宏觀力學參數(shù)與PFC2D細觀參數(shù)之間的關系式:

    表2 中心組合設計代碼表Table 2 Central composite design code table

    表3 細觀參數(shù)的中心組合設計代碼值Table 3 CCD code value of microscopic parameters

    其中,R2為相關系數(shù)。可見,相關系數(shù)R2均較高,表明擬合精度較高,能夠準確地反映巖石宏-細觀參數(shù)的關聯(lián)性。若已知巖石的宏觀力學參數(shù),則由式(2)可計算出PFC2D細觀參數(shù),其中Ec,kn/ks和σc由E,v和σt的關系式聯(lián)立求解,C1和μb通過相關系數(shù)R2較高的σf和φ關系式求得。

    3 試驗驗證

    3.1 巖石單軸和雙軸壓縮試驗

    為了驗證宏-細觀參數(shù)關系式(2),選取產(chǎn)自于云南的紅砂巖開展試驗研究。表5所示為本試驗測得的紅砂巖以及其他4 種巖石材料[22-23]的宏觀力學參數(shù)。在TRW-30001 軸伺服試驗機上開展巖石的單軸(直徑×高度為50 mm×100 mm的圓柱體試件)、雙軸(長×寬×高為100 mm×100 mm×40 mm 的立方體試件)壓縮試驗,其中雙軸壓縮中側(cè)向壓力為6 MPa,加載速率為500 N/s。在試驗過程中自動記錄應力-應變曲線。

    根據(jù)不同巖石的宏觀力學參數(shù)(表5)和模型宏觀力學參數(shù)與PFC2D細觀參數(shù)之間的關系式,可分別計算得到相應的PFC2D細觀參數(shù)(Ec,kn/ks,C1,σc和μb),如表6所示。采用該細觀參數(shù),分別開展不同巖石的單軸壓縮(E,ν,σf)、單軸拉伸(σt)和雙軸壓縮(C2,φ)模擬計算,得到的巖石宏觀力學參數(shù)(E,v,σf,σt,C2和φ)見表7,并與試驗結(jié)果進行對比。從表7可見:各參數(shù)的實驗結(jié)果與模擬結(jié)果相對誤差均小于10%,模擬精度較高,從而驗證了采用中心組合設計建立的模型宏觀力學參數(shù)與PFC2D細觀參數(shù)之間的關系式是合理、有效的。

    表4 中心組合設計試驗表Table 4 Central composite design table

    表5 不同巖石的宏觀力學參數(shù)Table 5 Macroscopic mechanical parameters of different rocks

    表6 不同巖石的PFC2D細觀參數(shù)計算結(jié)果Table 6 Calculated results of PFC2D microscopic parameters of different rocks

    圖10所示為紅砂巖在單軸、雙軸壓縮下數(shù)值模擬計算得到的應力-應變曲線,并與試驗結(jié)果進行對比。從圖10可見:在單軸、雙軸壓縮下,實測的紅砂巖應力-應變曲線(OABCD)可分為壓密階段(OA)、彈性變形階段(AB)、損傷階段(BC)和破壞階段(CD)。模擬結(jié)果(OB′C′D曲線)中沒有出現(xiàn)試驗曲線的壓密階段,這是因為PFC2D建模時顆粒為連續(xù)體,而實際巖石存在有天然的孔隙、裂紋等缺陷,在壓縮時會有微孔、微裂縫等閉合壓密階段(OA)。計算得到的紅砂巖單軸壓縮彈性模量、抗壓強度分別為11.38 GPa 和59.45 MPa,實測的單軸彈性模量、抗壓強度分別為10.69 GPa 和57.20 MPa,相對誤差分別為6.45%和3.93%;計算得到的雙軸壓縮彈性模量和抗壓強度分別為10.83 GPa和78.20 MPa,實測的雙軸彈性模量和抗壓強度分別為11.82 GPa 和82.40 MPa,相對誤差分別為8.38%和5.10%。可見,彈性模量和抗壓強度模擬結(jié)果(未考慮初始缺陷)均略高于實驗結(jié)果,雙軸壓縮比單軸壓縮的抗壓強度更高。忽略初始壓密階段,單、雙軸壓縮試驗曲線與模擬曲線基本重合,彈性模量、抗壓強度均較吻合,從而驗證了采用宏-細觀參數(shù)關系式(2)計算PFC2D 細觀參數(shù)方法的合理性和有效性。

    表7 不同巖石的宏觀力學參數(shù)(模擬與試驗結(jié)果對比)Table 7 Macroscopic mechanical parameters of different rocks(Comparison of simulation and test results)

    圖10 紅砂巖單、雙軸壓縮應力-應變曲線Fig.10 Stress-strain curves of sandstone

    3.2 巖石多裂紋斷裂試驗

    圖11 含3條平行裂紋的紅砂巖試件模型Fig.11 Model of sandstone specimen with three parallel cracks

    為了進一步驗證宏-細觀參數(shù)關系式(2),對含3條平行裂紋(裂紋長度為20 mm,寬度為1 mm)的紅砂巖試件進行斷裂試驗。采用與上述壓縮試驗相同的紅砂巖,宏觀力學參數(shù)、細觀參數(shù)與表5和表6所示的相同。試件是長×寬×高為100 mm×100 mm×20 mm的立方體,預制3條等長平行裂紋AB,CD和EF,裂紋長度為20 mm,寬度為1 mm,傾角α=45°,間距為15 mm,如圖11所示。在DNS100 伺服加載系統(tǒng)上進行單軸壓縮斷裂試驗,最大荷載為1 000 kN,以位移加載速率為0.1 mm/s進行加載,采用攝像機記錄試件斷裂破壞過程。

    圖12和圖13所示分別為單軸壓縮下紅砂巖試件細觀裂紋演變和宏觀裂紋擴展軌跡的數(shù)值模擬結(jié)果,其中,黑色、紅色、白色分別代表顆粒、細觀裂紋、宏觀裂紋。從圖12和圖13可見:在初始加載時,3 條裂紋的6 個裂紋尖端都會產(chǎn)生初始損傷而出現(xiàn)微裂紋(圖12(a)),此時沒有宏觀裂紋(圖13(a));隨著進一步加載,各裂紋尖端的損傷不斷積累會產(chǎn)生越來越多的微裂紋并集聚(圖12(b)),由于3條裂紋之間的損傷區(qū)相互作用,裂紋之間的微裂紋首先貫通而形成了最初的2 條宏觀裂紋BC和DE(圖13(b))。當外載不斷增加到破壞荷載時,上部裂紋AB的尖端A和下部裂紋EF的尖端F附近微裂紋急劇增多(圖12(c)),最后形成了2條沿軸向加載方向的宏觀主裂紋(尖端A向上延伸、尖端F向下延伸)以及衍生的次生裂紋,破壞軌跡如圖13(c)所示。

    圖12 單軸壓縮下含3條平行裂紋的紅砂巖試件細觀裂紋演變過程(模擬結(jié)果)Fig.12 Micro-crack evolutions of sandstone specimen with three parallel cracks under uniaxial compression(simulation results)

    圖13 單軸壓縮下含3條平行裂紋的紅砂巖試件宏觀裂紋擴展軌跡(模擬結(jié)果)Fig.13 Macro-crack trajectories of sandstone specimen with three parallel cracks under uniaxial compression(simulation results)

    圖14所示為攝像機記錄的單軸壓縮下紅砂巖試件裂紋擴展軌跡試驗結(jié)果。從圖14可見:在加載初期,各裂紋尖端會有局部損傷產(chǎn)生微裂紋,但沒有出現(xiàn)宏觀裂紋(圖14(a));隨著外載增加,3條裂紋之間的微裂紋聚集、演變、貫通而產(chǎn)生了2條宏觀裂紋;當外載不斷增加到破壞荷載時,裂紋尖端A和尖端F的微裂紋不斷發(fā)育,最后形成了沿軸向方向延伸的2 條主裂紋(尖端A向上延伸、尖端F向下延伸)和衍生的次生裂紋,破壞軌跡如圖14(c)所示,其中,試件的上部碎片是試驗機卸載時能量釋放沖壓所致。

    比較圖13和圖14可見:除了最后破壞階段的次生裂紋有細微差別外(這是實際巖石材料的非均質(zhì)性所致),單軸壓縮下含3 條平行裂紋的紅砂巖試件裂紋擴展軌跡試驗結(jié)果與PFC2D模擬結(jié)果較吻合,從而進一步驗證了采用宏-細觀參數(shù)關系式(2)計算PFC2D細觀參數(shù)方法的合理性和可行性。

    圖14 單軸壓縮下含3條平行裂紋的紅砂巖試件裂紋擴展軌跡(試驗結(jié)果)Fig.14 Propagation trajectories of sandstone specimen with three parallel cracks under uniaxial compression(test results)

    4 結(jié)論

    1)采用中心組合設計法和回歸分析,建立了巖石的PFC2D平直節(jié)理模型宏-細觀參數(shù)關系式,可定量計算巖石的細觀參數(shù)。與現(xiàn)有的試錯法和正交設計方法相比,該方法試驗次數(shù)較少,精度較高。

    2)PFC2D細觀參數(shù)對巖石的宏觀力學參數(shù)的影響規(guī)律為:Ec對宏觀彈性模量E影響最大,kn/ks對泊松比ν影響最大,C1對σf和C2影響最大,μb對φ影響最大,它們各自呈正相關關系。

    3)通過開展單、雙軸壓縮試驗和斷裂試驗,測得了完整巖石試件的應力-應變曲線、宏觀力學參數(shù)和含多裂紋巖石試件的斷裂軌跡,該試驗結(jié)果與PFC2D模擬計算結(jié)果較吻合,從而驗證了所建立的PFC2D平直節(jié)理模型宏-細觀參數(shù)關系式的有效性。

    4)采用中心組合設計來確定巖石PFC2D平直節(jié)理模型宏-細觀參數(shù)關系的方法,可推廣到其他脆性材料的細觀參數(shù)確定,為PFC 定量分析多裂紋多場耦合斷裂的細觀機理提供可靠依據(jù)。

    猜你喜歡
    細觀節(jié)理單軸
    基于細觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    單軸壓縮條件下巖石峰后第Ⅱ種類型應力——應變曲線的新解釋
    新疆阜康白楊河礦區(qū)古構(gòu)造應力場特征
    CFRP-鋼復合板的單軸拉伸力學性能
    工程與建設(2019年1期)2019-09-03 01:12:38
    新疆阜康白楊河礦區(qū)構(gòu)造節(jié)理發(fā)育特征
    中國煤層氣(2018年3期)2018-07-27 11:41:42
    單軸應變Si NMOS電流模型研究
    電子測試(2017年12期)2017-12-18 06:35:42
    Effect of Magnetic Field on Forced Convection between Two Nanofluid Laminar Flows in a Channel
    斜單軸跟蹤式光伏組件的安裝傾角優(yōu)化設計
    基于四叉樹網(wǎng)格加密技術(shù)的混凝土細觀模型
    PBX炸藥的抗壓強度及抗拉強度細觀尺度的數(shù)值計算
    火炸藥學報(2014年1期)2014-03-20 13:17:25
    久久国产精品大桥未久av| 男女边吃奶边做爰视频| 日韩,欧美,国产一区二区三区| 一级黄片播放器| 老司机影院成人| 国产淫语在线视频| 国产精品久久久久久精品电影小说| 91老司机精品| 十八禁人妻一区二区| 在线看a的网站| 90打野战视频偷拍视频| 亚洲中文字幕日韩| 女人爽到高潮嗷嗷叫在线视频| 女警被强在线播放| 美女主播在线视频| 国产精品秋霞免费鲁丝片| 国产精品久久久久久人妻精品电影 | 免费在线观看黄色视频的| 中文欧美无线码| 日韩人妻精品一区2区三区| 国产男女超爽视频在线观看| 国产三级黄色录像| 亚洲,欧美,日韩| 午夜福利乱码中文字幕| 日韩精品免费视频一区二区三区| 亚洲欧美色中文字幕在线| 日本91视频免费播放| 乱人伦中国视频| 纯流量卡能插随身wifi吗| 久热爱精品视频在线9| 国产精品一区二区精品视频观看| 曰老女人黄片| 曰老女人黄片| 美女国产高潮福利片在线看| 国产在视频线精品| 真人做人爱边吃奶动态| 高潮久久久久久久久久久不卡| 久久影院123| 成年女人毛片免费观看观看9 | 高清视频免费观看一区二区| 国产欧美日韩一区二区三 | 国语对白做爰xxxⅹ性视频网站| 日日摸夜夜添夜夜爱| netflix在线观看网站| 亚洲中文日韩欧美视频| 99国产精品99久久久久| 国产精品久久久久久人妻精品电影 | 国产91精品成人一区二区三区 | 天堂8中文在线网| 91九色精品人成在线观看| 亚洲欧美一区二区三区国产| 亚洲国产精品国产精品| 天堂俺去俺来也www色官网| 丝袜在线中文字幕| 国产野战对白在线观看| 久久久国产欧美日韩av| videos熟女内射| 精品人妻在线不人妻| www.999成人在线观看| 男女之事视频高清在线观看 | 午夜激情久久久久久久| 黑人欧美特级aaaaaa片| 一本综合久久免费| 悠悠久久av| 美女午夜性视频免费| 久久精品熟女亚洲av麻豆精品| 国产亚洲欧美在线一区二区| 狠狠婷婷综合久久久久久88av| 国产精品 国内视频| 大话2 男鬼变身卡| 午夜av观看不卡| 久久这里只有精品19| 大片电影免费在线观看免费| 久久99热这里只频精品6学生| 日韩免费高清中文字幕av| 亚洲三区欧美一区| 男女午夜视频在线观看| 免费看不卡的av| 国产av精品麻豆| 亚洲精品久久久久久婷婷小说| av福利片在线| 七月丁香在线播放| 男女边吃奶边做爰视频| 国产一区二区 视频在线| 国产av国产精品国产| 制服诱惑二区| 十八禁网站网址无遮挡| 人人妻人人爽人人添夜夜欢视频| 在线观看免费视频网站a站| 欧美激情 高清一区二区三区| 国精品久久久久久国模美| 欧美大码av| 久久综合国产亚洲精品| 久久 成人 亚洲| av电影中文网址| 午夜两性在线视频| 少妇粗大呻吟视频| 女人被躁到高潮嗷嗷叫费观| 丝瓜视频免费看黄片| 久久精品aⅴ一区二区三区四区| 成年人免费黄色播放视频| 国产精品av久久久久免费| 国产成人精品久久二区二区91| 久久性视频一级片| 操出白浆在线播放| 熟女少妇亚洲综合色aaa.| 高清av免费在线| 日本午夜av视频| 精品久久久久久电影网| 熟女少妇亚洲综合色aaa.| 国产亚洲欧美精品永久| 日本色播在线视频| 色精品久久人妻99蜜桃| 人成视频在线观看免费观看| 天天躁夜夜躁狠狠久久av| 如日韩欧美国产精品一区二区三区| 亚洲天堂av无毛| 日韩大片免费观看网站| 男女高潮啪啪啪动态图| 亚洲欧美日韩高清在线视频 | 精品亚洲成国产av| 新久久久久国产一级毛片| 日韩视频在线欧美| 国产成人一区二区在线| 叶爱在线成人免费视频播放| 丝袜美腿诱惑在线| 亚洲av日韩在线播放| 日韩av在线免费看完整版不卡| 成人黄色视频免费在线看| 欧美黄色淫秽网站| 国产视频首页在线观看| 久久人人爽av亚洲精品天堂| 久久国产亚洲av麻豆专区| 99国产精品一区二区三区| 久久国产精品影院| 好男人电影高清在线观看| 日本欧美视频一区| 天天添夜夜摸| 国产一卡二卡三卡精品| 亚洲av电影在线进入| 欧美在线一区亚洲| 视频区图区小说| 久久久久国产一级毛片高清牌| 久久影院123| 99九九在线精品视频| 一区二区三区激情视频| 高清黄色对白视频在线免费看| 美女脱内裤让男人舔精品视频| 亚洲欧美精品自产自拍| 国产精品国产三级专区第一集| 欧美日韩一级在线毛片| 色综合欧美亚洲国产小说| 久久午夜综合久久蜜桃| 高清不卡的av网站| 成年人午夜在线观看视频| 国产老妇伦熟女老妇高清| 另类精品久久| 亚洲av在线观看美女高潮| 在线观看人妻少妇| 91国产中文字幕| 最黄视频免费看| avwww免费| 久久久精品免费免费高清| 久久毛片免费看一区二区三区| 色94色欧美一区二区| 欧美在线一区亚洲| 亚洲熟女精品中文字幕| 成人国语在线视频| 青春草亚洲视频在线观看| 日日摸夜夜添夜夜爱| netflix在线观看网站| 少妇 在线观看| 伊人久久大香线蕉亚洲五| 久久精品熟女亚洲av麻豆精品| 国产精品二区激情视频| 精品免费久久久久久久清纯 | 亚洲av美国av| 亚洲国产日韩一区二区| 免费观看av网站的网址| 欧美日韩视频精品一区| 久久久国产一区二区| 日日摸夜夜添夜夜爱| 日韩一区二区三区影片| 一个人免费看片子| 日本黄色日本黄色录像| 久久久久久久久久久久大奶| 在线观看免费午夜福利视频| av又黄又爽大尺度在线免费看| 在线观看国产h片| 天天躁日日躁夜夜躁夜夜| 久久久精品区二区三区| 欧美av亚洲av综合av国产av| 国产一区二区三区综合在线观看| 最近最新中文字幕大全免费视频 | 色婷婷久久久亚洲欧美| 久久久精品区二区三区| 国产精品麻豆人妻色哟哟久久| 午夜福利影视在线免费观看| 飞空精品影院首页| 蜜桃在线观看..| 高清黄色对白视频在线免费看| 精品一区二区三区四区五区乱码 | 女人爽到高潮嗷嗷叫在线视频| 侵犯人妻中文字幕一二三四区| 在线观看人妻少妇| www.av在线官网国产| 亚洲精品日韩在线中文字幕| 在线观看免费午夜福利视频| 丰满人妻熟妇乱又伦精品不卡| 美女福利国产在线| 亚洲欧美中文字幕日韩二区| 纯流量卡能插随身wifi吗| 成年动漫av网址| 免费在线观看完整版高清| 免费观看av网站的网址| 亚洲色图 男人天堂 中文字幕| 天堂中文最新版在线下载| 满18在线观看网站| 欧美成狂野欧美在线观看| 一区二区av电影网| 一级片'在线观看视频| 亚洲国产最新在线播放| 中文乱码字字幕精品一区二区三区| 国产男人的电影天堂91| 亚洲中文日韩欧美视频| 久热爱精品视频在线9| 免费在线观看视频国产中文字幕亚洲 | 精品视频人人做人人爽| 午夜福利,免费看| 欧美精品一区二区大全| 日本一区二区免费在线视频| 桃花免费在线播放| 一本综合久久免费| 午夜影院在线不卡| 最近中文字幕2019免费版| 性少妇av在线| 人人澡人人妻人| 老司机影院成人| 在线亚洲精品国产二区图片欧美| 欧美97在线视频| 免费观看人在逋| 亚洲美女黄色视频免费看| 人人澡人人妻人| 五月天丁香电影| 国产精品麻豆人妻色哟哟久久| 美女脱内裤让男人舔精品视频| 免费人妻精品一区二区三区视频| 少妇的丰满在线观看| 黑人欧美特级aaaaaa片| netflix在线观看网站| 一区二区三区激情视频| 两个人免费观看高清视频| 国产成人精品久久二区二区91| 99热网站在线观看| 青春草亚洲视频在线观看| 久久99热这里只频精品6学生| 久久精品国产综合久久久| 欧美性长视频在线观看| 90打野战视频偷拍视频| 国产精品九九99| 欧美xxⅹ黑人| 免费观看a级毛片全部| 啦啦啦在线观看免费高清www| 久久国产亚洲av麻豆专区| 男男h啪啪无遮挡| 久久久久网色| 婷婷色av中文字幕| 国产真人三级小视频在线观看| 日韩视频在线欧美| 欧美日韩亚洲综合一区二区三区_| 波多野结衣一区麻豆| 搡老乐熟女国产| 超碰成人久久| 亚洲精品国产av蜜桃| 久久鲁丝午夜福利片| 成人手机av| 91老司机精品| 久久久久国产精品人妻一区二区| 男人添女人高潮全过程视频| 亚洲国产精品一区二区三区在线| 午夜老司机福利片| 国产精品欧美亚洲77777| 亚洲中文日韩欧美视频| 欧美人与性动交α欧美软件| 亚洲熟女毛片儿| 女人被躁到高潮嗷嗷叫费观| 美女中出高潮动态图| 99国产精品免费福利视频| kizo精华| 肉色欧美久久久久久久蜜桃| 七月丁香在线播放| 亚洲av电影在线观看一区二区三区| 精品国产国语对白av| 亚洲专区中文字幕在线| videosex国产| 午夜av观看不卡| 熟女少妇亚洲综合色aaa.| 亚洲,一卡二卡三卡| 国产精品熟女久久久久浪| 国产成人欧美| 亚洲国产毛片av蜜桃av| 亚洲成人国产一区在线观看 | 十分钟在线观看高清视频www| 高清欧美精品videossex| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久久毛片微露脸 | 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 美女高潮到喷水免费观看| 日本午夜av视频| 美女福利国产在线| 日韩 欧美 亚洲 中文字幕| 男女之事视频高清在线观看 | 精品久久久精品久久久| 亚洲精品第二区| 在线 av 中文字幕| h视频一区二区三区| 亚洲av欧美aⅴ国产| 久久青草综合色| 精品亚洲成a人片在线观看| 亚洲一码二码三码区别大吗| 青春草亚洲视频在线观看| 飞空精品影院首页| 国产精品二区激情视频| 观看av在线不卡| 人人澡人人妻人| 亚洲精品一区蜜桃| 亚洲欧美一区二区三区黑人| 夫妻午夜视频| 日韩制服丝袜自拍偷拍| 免费观看av网站的网址| 丰满人妻熟妇乱又伦精品不卡| 老司机亚洲免费影院| 日日夜夜操网爽| 精品亚洲成a人片在线观看| 电影成人av| 国产人伦9x9x在线观看| 精品少妇一区二区三区视频日本电影| 午夜福利一区二区在线看| 欧美日韩成人在线一区二区| 热re99久久国产66热| 丝瓜视频免费看黄片| 岛国毛片在线播放| 免费在线观看黄色视频的| 高清不卡的av网站| 亚洲 国产 在线| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 一级毛片我不卡| 欧美精品高潮呻吟av久久| 高清黄色对白视频在线免费看| 久久久国产一区二区| 国产精品99久久99久久久不卡| 人妻 亚洲 视频| 久久精品亚洲av国产电影网| 久久99热这里只频精品6学生| 国产爽快片一区二区三区| 国产免费现黄频在线看| 少妇人妻久久综合中文| av不卡在线播放| 免费少妇av软件| 久久久亚洲精品成人影院| 国产日韩一区二区三区精品不卡| 久久精品久久久久久久性| 久久久久精品人妻al黑| 一区二区av电影网| 欧美人与性动交α欧美软件| 一区福利在线观看| 亚洲人成电影观看| 亚洲精品中文字幕在线视频| 十八禁网站网址无遮挡| 亚洲欧美成人综合另类久久久| 人妻人人澡人人爽人人| 久久国产精品男人的天堂亚洲| 伊人亚洲综合成人网| 亚洲av电影在线观看一区二区三区| 老司机靠b影院| 夫妻性生交免费视频一级片| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 国产91精品成人一区二区三区 | 老司机在亚洲福利影院| avwww免费| 久久性视频一级片| 啦啦啦视频在线资源免费观看| 精品人妻1区二区| 亚洲七黄色美女视频| www.熟女人妻精品国产| 免费在线观看影片大全网站 | 女警被强在线播放| 亚洲欧美清纯卡通| 赤兔流量卡办理| 麻豆国产av国片精品| 精品国产国语对白av| 日韩电影二区| 欧美成狂野欧美在线观看| 丝袜美足系列| 欧美在线黄色| 久久性视频一级片| 99久久99久久久精品蜜桃| 国产免费又黄又爽又色| 久久99精品国语久久久| 久久国产亚洲av麻豆专区| 国产在线观看jvid| 18在线观看网站| 真人做人爱边吃奶动态| 久久国产精品男人的天堂亚洲| 国产成人影院久久av| 日日摸夜夜添夜夜爱| av在线老鸭窝| 汤姆久久久久久久影院中文字幕| netflix在线观看网站| 男女午夜视频在线观看| 色精品久久人妻99蜜桃| 真人做人爱边吃奶动态| 国产xxxxx性猛交| 一区二区三区精品91| 99久久综合免费| 在线观看免费视频网站a站| 欧美日韩视频精品一区| 日韩免费高清中文字幕av| 亚洲人成网站在线观看播放| 香蕉国产在线看| 日韩免费高清中文字幕av| 18禁国产床啪视频网站| 国产欧美日韩综合在线一区二区| 热99国产精品久久久久久7| 蜜桃在线观看..| a级毛片在线看网站| 久久久久久久精品精品| 国产在线观看jvid| 欧美日韩视频精品一区| 亚洲自偷自拍图片 自拍| 国产精品国产av在线观看| 男女边摸边吃奶| 真人做人爱边吃奶动态| 一区二区日韩欧美中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机深夜福利视频在线观看 | 国产亚洲精品第一综合不卡| 亚洲精品久久久久久婷婷小说| av国产久精品久网站免费入址| 久久毛片免费看一区二区三区| 国产国语露脸激情在线看| 成人手机av| 我要看黄色一级片免费的| 欧美人与善性xxx| 中文字幕色久视频| 国产无遮挡羞羞视频在线观看| 一级毛片黄色毛片免费观看视频| 国产精品人妻久久久影院| 色精品久久人妻99蜜桃| 免费女性裸体啪啪无遮挡网站| 七月丁香在线播放| 女性被躁到高潮视频| 国产成人欧美| 国产精品九九99| 七月丁香在线播放| 日韩 亚洲 欧美在线| 国产伦理片在线播放av一区| 亚洲五月婷婷丁香| 亚洲av片天天在线观看| 久久毛片免费看一区二区三区| 满18在线观看网站| 精品第一国产精品| 韩国精品一区二区三区| 国产成人精品无人区| 青青草视频在线视频观看| 人人妻,人人澡人人爽秒播 | 一级片免费观看大全| 91老司机精品| 精品亚洲成国产av| 久久99一区二区三区| 我的亚洲天堂| av在线老鸭窝| 国产黄频视频在线观看| 国产1区2区3区精品| 十分钟在线观看高清视频www| 久久女婷五月综合色啪小说| 国产有黄有色有爽视频| 国产一区二区三区综合在线观看| 亚洲久久久国产精品| 亚洲激情五月婷婷啪啪| 免费久久久久久久精品成人欧美视频| 亚洲精品日本国产第一区| 亚洲国产精品999| 国产一区二区三区av在线| 亚洲自偷自拍图片 自拍| 成年动漫av网址| 免费高清在线观看视频在线观看| 精品高清国产在线一区| 精品亚洲成a人片在线观看| 高清黄色对白视频在线免费看| 亚洲av电影在线观看一区二区三区| 国产精品99久久99久久久不卡| 人妻 亚洲 视频| 香蕉国产在线看| 亚洲,一卡二卡三卡| tube8黄色片| 中国美女看黄片| 脱女人内裤的视频| 天天影视国产精品| 电影成人av| av国产精品久久久久影院| 亚洲成人免费av在线播放| 男的添女的下面高潮视频| 激情五月婷婷亚洲| 中文字幕亚洲精品专区| 亚洲一码二码三码区别大吗| 久久影院123| 女人高潮潮喷娇喘18禁视频| 亚洲精品日韩在线中文字幕| 啦啦啦视频在线资源免费观看| 国产野战对白在线观看| 国产色视频综合| 无限看片的www在线观看| 亚洲成av片中文字幕在线观看| 国产精品人妻久久久影院| 亚洲 国产 在线| 久久久亚洲精品成人影院| 亚洲国产欧美网| 午夜福利,免费看| 日韩制服骚丝袜av| 建设人人有责人人尽责人人享有的| 国产成人精品在线电影| 亚洲国产中文字幕在线视频| 亚洲少妇的诱惑av| 国产黄色免费在线视频| 色婷婷久久久亚洲欧美| 国产在线视频一区二区| 国产精品三级大全| 国产爽快片一区二区三区| 一级片免费观看大全| 国产av国产精品国产| 一区二区日韩欧美中文字幕| 亚洲国产av新网站| 亚洲成人免费电影在线观看 | 亚洲av综合色区一区| 国产精品久久久久久人妻精品电影 | 日韩伦理黄色片| 亚洲精品第二区| 最近中文字幕2019免费版| 精品一区二区三区av网在线观看 | 亚洲av美国av| 又粗又硬又长又爽又黄的视频| 中文字幕亚洲精品专区| 汤姆久久久久久久影院中文字幕| 亚洲av片天天在线观看| 国产精品国产av在线观看| 国产成人av教育| 午夜两性在线视频| 亚洲黑人精品在线| 免费一级毛片在线播放高清视频 | 男女边摸边吃奶| 亚洲欧洲精品一区二区精品久久久| 蜜桃国产av成人99| 美女脱内裤让男人舔精品视频| 男女高潮啪啪啪动态图| 狠狠精品人妻久久久久久综合| 欧美精品人与动牲交sv欧美| 高清不卡的av网站| 亚洲人成网站在线观看播放| 国产成人91sexporn| 亚洲精品第二区| 51午夜福利影视在线观看| 国产99久久九九免费精品| kizo精华| 国产精品免费大片| 97在线人人人人妻| 久久99精品国语久久久| 热99国产精品久久久久久7| 看免费成人av毛片| 一级毛片电影观看| 搡老岳熟女国产| 午夜老司机福利片| 欧美成人午夜精品| 久久免费观看电影| 国产免费现黄频在线看| 国产高清videossex| 国产日韩欧美亚洲二区| 两性夫妻黄色片| av不卡在线播放| 波野结衣二区三区在线| 亚洲av电影在线观看一区二区三区| 欧美av亚洲av综合av国产av| 久久九九热精品免费| 亚洲欧美精品综合一区二区三区| 免费高清在线观看日韩| 免费在线观看日本一区| 免费在线观看完整版高清| 校园人妻丝袜中文字幕| 看免费成人av毛片| 亚洲成人手机| 欧美日韩综合久久久久久| 咕卡用的链子| 中文字幕另类日韩欧美亚洲嫩草| 久久精品亚洲熟妇少妇任你| 美女大奶头黄色视频| 亚洲精品第二区| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩亚洲高清精品| 最黄视频免费看| 十八禁高潮呻吟视频| 18禁裸乳无遮挡动漫免费视频| 久久国产精品男人的天堂亚洲| 各种免费的搞黄视频| 欧美在线黄色| 精品一区在线观看国产| 亚洲第一青青草原| 看免费av毛片| 99国产综合亚洲精品| 一边摸一边抽搐一进一出视频| 男女免费视频国产| 九色亚洲精品在线播放| 91国产中文字幕|