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

    基于固支直桿彎曲小試樣的P91/P92鋼蠕變性能對(duì)比評(píng)價(jià)研究

    2021-05-21 07:30:54余海洋周幗彥涂善東
    壓力容器 2021年4期
    關(guān)鍵詞:直桿固支壓頭

    韓 笑,余海洋,周幗彥,涂善東

    (華東理工大學(xué) 承壓系統(tǒng)與安全教育部重點(diǎn)實(shí)驗(yàn)室,上海 200237)

    0 引言

    近年來(lái),超(超)臨界機(jī)組已成為我國(guó)火力發(fā)電主流機(jī)組,隨著火電機(jī)組參數(shù)和效率的提高,機(jī)組所選用鋼材等級(jí)不斷提高。P91/P92鋼由于其較高的性價(jià)比及良好的性能,成為超(超)臨界機(jī)組關(guān)鍵部件的主要材料類型[1],包括主蒸汽管道、再熱蒸汽管道以及高溫聯(lián)箱等,如表1[2]所示。超(超)臨界機(jī)組關(guān)鍵部件的服役溫度范圍為600~700 ℃,服役壓力30.0 MPa,長(zhǎng)期服役于該高溫高壓工況下,這些結(jié)構(gòu)材料的性能劣化甚至斷裂,不但縮短機(jī)組壽命,還會(huì)對(duì)熱電廠的高效生產(chǎn)及安全性造成嚴(yán)重影響。其中蠕變斷裂為主要的失效形式。因此對(duì)比分析P91/P92材料的高溫蠕變性能及其影響因素,對(duì)火電機(jī)組關(guān)鍵部件的合理選材、高效設(shè)計(jì)以及安全運(yùn)行具有重要意義。

    表1 超(超)臨界鍋爐主要部件用鋼

    目前,國(guó)內(nèi)外已對(duì)P91/P92鋼材的蠕變性能進(jìn)行了廣泛的研究。BURTON等[3]提出了純扭轉(zhuǎn)的多軸螺旋彈簧蠕變?cè)囼?yàn)法,并用于研究低應(yīng)力下P91鋼的短程蠕變行為,分析了低應(yīng)力下組織退化對(duì)蠕變性能的影響。許樂(lè)等[4]對(duì)P92鋼焊接接頭蠕變損傷與裂紋擴(kuò)展進(jìn)行了有限元模擬,建立了裂紋萌生/擴(kuò)展行為與熱影響區(qū)寬度的關(guān)系。黃橋生等[5]對(duì)P92鋼焊接接頭進(jìn)行高溫持久試驗(yàn),從微觀角度分析了P92鋼焊接接頭蠕變損傷的形成機(jī)理。BENDICK等[6]總結(jié)了P92鋼的蠕變強(qiáng)度外推方法,擴(kuò)大了P92鋼蠕變強(qiáng)度的預(yù)測(cè)范圍。畢瑤等[7]通過(guò)蠕變持久和間斷試驗(yàn)制備了具有不同蠕變損傷狀態(tài)的P91鋼試樣,并利用3.5%NaCl電解液研究了P91高溫蠕變損傷過(guò)程中的電化學(xué)極化行為,得到了蠕變損傷速率與組織的關(guān)系。上述這些研究都是將P91/P92系列材料作為一類材料,并未對(duì)兩者進(jìn)行具體的對(duì)比分析。然而二者化學(xué)成分有所不同,所表現(xiàn)出的力學(xué)性能以及對(duì)溫度和壓力的敏感性也不盡相同,導(dǎo)致其構(gòu)件壽命和安全性的較大差異。因此,定量對(duì)比分析評(píng)估兩種材料的高溫力學(xué)性能,對(duì)超(超)臨界機(jī)組的精確選材和安全生產(chǎn)具有關(guān)鍵性的決定作用。

    為了盡可能減少對(duì)在役設(shè)備的損傷,國(guó)內(nèi)外也不乏利用小試樣方法對(duì)構(gòu)件材料P91/P92的高溫性能進(jìn)行研究。秦炎鋒等[8]采用P91材料,利用小沖桿法確定了該材料在一定溫度下的Norton 參數(shù);AL-ABEDY等[9]利用小沖桿試樣,采用Kocks-Mecking Estrin模型評(píng)價(jià)了P91高溫蠕變性能。SONG等[10]對(duì)比了小沖桿試樣與單軸試樣下的P92材料蠕變損傷特性,證明了小沖桿試樣評(píng)價(jià)P92鋼的蠕變性能相對(duì)于單軸試樣的優(yōu)越性。JEONG等[11]采用小沖桿蠕變的方法,對(duì)P92微觀析出相和空洞的變化進(jìn)行定量分析,得到了不同試驗(yàn)方式及加載方式對(duì)P92鋼析出相和空洞變化的影響規(guī)律。然而相關(guān)研究表明,小沖桿試驗(yàn)法具有一定的局限性,其試樣受力復(fù)雜,解析解較為困難,蠕變第三階段非常短暫,可采集數(shù)據(jù)有限。相比之下,固支直桿彎曲小試樣試驗(yàn)法[12-15]試樣結(jié)構(gòu)和受力簡(jiǎn)單,能夠獲得完整的蠕變?nèi)A段曲線,在評(píng)價(jià)在役和非在役設(shè)備材料性能方面具有潛在的優(yōu)越性[16]。因此,可以預(yù)見利用固支直桿小試樣方法對(duì)P91/P92材料進(jìn)行對(duì)比分析研究更方便,但其應(yīng)用的可行性還需進(jìn)一步驗(yàn)證。

    本文考慮到經(jīng)濟(jì)性、耗時(shí)長(zhǎng)短和試驗(yàn)精度等方面,主要基于固支直桿彎曲小試樣方法,再在確定P91/P92材料合適的試樣模型基礎(chǔ)上,采用有限元方法對(duì)P91/P92材料的蠕變性能進(jìn)行對(duì)比分析,并探討其對(duì)載荷及溫度的敏感性,為超(超)臨界機(jī)組用鋼正確經(jīng)濟(jì)的選取提供依據(jù)。

    1 P91/P92 鋼的理化性能

    1.1 P91/P92鋼化學(xué)成分

    P91是在9Cr-1Mo鋼的基礎(chǔ)上,以Cr,Mo為主要固溶強(qiáng)化合金元素,在限制C含量的同時(shí),嚴(yán)格控制P和S等殘余元素含量,并添加N,V和Nb等元素形成的低碳低合金耐熱鋼,現(xiàn)階段的服役溫度為600~650 ℃,屬于回火馬氏體鋼[17]。隨著高參數(shù)、大容量超(超)臨界機(jī)組的發(fā)展,P91鋼的性能已經(jīng)不能夠滿足日漸增長(zhǎng)的蒸汽溫度和壓力的要求,于是P92鋼應(yīng)運(yùn)而生。P92鋼是在P91鋼基礎(chǔ)上開發(fā)的新一代F/M鋼,主要采用增加W、減少M(fèi)o,同時(shí)加入少量的B及V,Nb的復(fù)合多元強(qiáng)化的手段開發(fā)出的[18],其高溫蠕變強(qiáng)度和服役溫度進(jìn)一步提高。根據(jù)ASTM A-335,P91/P92鋼的化學(xué)成分如表2[19]所示。

    表2 P91與P92鋼材料化學(xué)成分

    1.2 P91/P92鋼力學(xué)性能參數(shù)

    ASME規(guī)范中,P91/P92鋼的常溫力學(xué)性能如表3[20]所示。

    表3 P91與P92鋼的力學(xué)性能

    2 試樣模型參數(shù)的確定

    2.1 有限元模型的建立

    本文采用ABAQUS進(jìn)行數(shù)值模擬計(jì)算,對(duì)P91鋼與P92鋼的高溫蠕變性能進(jìn)行模擬對(duì)比分析。試樣的有限元模型如圖1所示,試樣為矩形截面直桿試樣,定義試樣厚度為2h、寬度b、總長(zhǎng)L、有效跨距l(xiāng)、壓頭半徑r。

    圖1 固支直桿彎曲小試樣有限元模型

    蠕變本構(gòu)采用Norton方程:

    (1)

    試樣兩端面直接固定(U1=U2=U3=0),恒定載荷P通過(guò)圓柱形壓頭施加在試樣中心,約束壓頭除y方向以外的全部自由度。有限元模擬過(guò)程中摩擦系數(shù)取f=0.3,試樣網(wǎng)格選用C3D8R類型,對(duì)試樣壓頭中心和固定端等變形較大位置的網(wǎng)格進(jìn)行加密,并且關(guān)閉幾何非線性。

    在本構(gòu)模型方面,文獻(xiàn)[12-13]中基于不同材料的試驗(yàn)結(jié)果,分析并修正了基于梁彎曲理論的固支直桿彎曲小試樣模型,與單軸蠕變?cè)囼?yàn)結(jié)果相比,其準(zhǔn)確度較高。故本文采用了修正后的固支直桿彎曲小試樣本構(gòu)模型進(jìn)行分析,其等效應(yīng)力和應(yīng)變公式[13]如下:

    (2)

    (3)

    為進(jìn)一步說(shuō)明固支直桿彎曲小試樣本構(gòu)模型對(duì)于本文所研究的P91/P92材料的有效性,采用文獻(xiàn)[21]中P91材料在650 ℃條件下的試驗(yàn)結(jié)果進(jìn)行驗(yàn)證。圖2為650 ℃下P91單軸拉伸試驗(yàn)與有限元模擬穩(wěn)態(tài)蠕變應(yīng)變速率對(duì)比圖。由有限元模擬反演得到650 ℃下P91材料的Norton蠕變參數(shù)為:n=8.41,B=1.17×10-20,與文獻(xiàn)[21]中試驗(yàn)結(jié)果n=8.234,B=2.472×10-20非常接近,蠕變指數(shù)n誤差僅為2.11%,蠕變常數(shù)B在相同數(shù)量級(jí)。說(shuō)明采用該本構(gòu)模型和有限元方法對(duì)P91材料進(jìn)行蠕變性能的模擬分析是可行且準(zhǔn)確的。

    圖2 650 ℃下P91單軸拉伸試驗(yàn)與有限元模擬關(guān)系對(duì)比Fig.2 Comparison diagram of relationship between uniaxial tensile test and finite element simulation of P91 at 650 ℃

    2.2 固支直桿彎曲小試樣尺寸參數(shù)的確定

    目前,固支直桿彎曲小試樣方法的試驗(yàn)參數(shù)尚未統(tǒng)一,各個(gè)研究者選用的試樣尺寸也不盡相同,導(dǎo)致得到的結(jié)果可比性較差。試樣的壓頭尺寸及對(duì)中度、試樣的尺寸、夾持端的結(jié)構(gòu)等的不同,也會(huì)影響固支直桿彎曲評(píng)價(jià)蠕變性能的準(zhǔn)確性。因此,為了排除試樣尺寸對(duì)模擬結(jié)果的影響,保證結(jié)果的有效性及準(zhǔn)確性,本文基于有限元法,運(yùn)用上述模型,精確地控制各個(gè)參數(shù)的變化,模擬分析不同參數(shù)對(duì)固支直桿三點(diǎn)彎曲試樣試驗(yàn)評(píng)價(jià)超(超)臨界機(jī)組用鋼蠕變性能的影響,以確定出最合適的模型參數(shù)。

    2.2.1 試樣有效跨距的確定

    由圖3所示的固支直桿三點(diǎn)彎曲試樣受力簡(jiǎn)圖可以看出,試樣受集中載荷發(fā)生橫力彎曲,試樣跨距同厚度滿足一定條件成為細(xì)長(zhǎng)梁時(shí),切應(yīng)力的影響可以忽略,此時(shí)橫力彎曲問(wèn)題簡(jiǎn)化為純彎曲問(wèn)題[22],三點(diǎn)彎梁橫截面上的最大切應(yīng)力與最大正應(yīng)力之比為:

    (4)

    當(dāng)2h∶l=1∶10時(shí),最大切應(yīng)力僅為正應(yīng)力的10%,可以忽略??紤]固支直桿三點(diǎn)彎曲小試樣固定端約束較大,選取最大比例為:2h∶l=1∶12。為了考慮有效跨距的影響作用,分別選取10,12,14,16四組尺寸,利用P91鋼分別在兩種最優(yōu)寬厚比(b∶2h=1∶1和b∶2h=2∶1)[15]下進(jìn)行模擬。P91鋼穩(wěn)態(tài)蠕變應(yīng)變速率與應(yīng)力的曲線見圖4。

    圖4 P91材料穩(wěn)態(tài)蠕變應(yīng)變速率與應(yīng)力σ關(guān)系曲線

    模擬結(jié)果與單軸拉伸試樣的穩(wěn)態(tài)位移速率誤差如表4所示。隨著有效跨距的增大,三點(diǎn)彎曲試樣的模擬結(jié)果與單軸拉伸試樣結(jié)果的偏差減小,但有效跨距對(duì)模擬結(jié)果的影響不大,偏差減小量均小于0.04,所以考慮到盡量減小試樣的尺寸,選取有效跨距為l=10 mm和l=12 mm兩種尺寸進(jìn)行后續(xù)模擬。

    表4 不同有效跨距的試樣與單軸試驗(yàn)對(duì)比誤差

    2.2.2 試樣具體寬厚尺寸的確定

    文獻(xiàn)[15]中經(jīng)過(guò)與單軸試驗(yàn)的穩(wěn)態(tài)應(yīng)變速率-應(yīng)力曲線對(duì)比可知,當(dāng)橫截面的寬厚比b∶2h=2∶1時(shí),選擇過(guò)渡段傾斜角(如圖5所示)為30°的試樣最優(yōu);當(dāng)橫截面的寬厚比b∶2h=1∶1時(shí),選擇過(guò)渡段傾斜角為45°的試樣最優(yōu)。

    圖5 過(guò)渡段結(jié)構(gòu)示意

    因此,選擇固定的寬厚比為b∶2h=2∶1和b∶2h=1∶1兩種比例,并改變具體幾何尺寸,利用P91材料進(jìn)行模擬,分析具體尺寸對(duì)模擬結(jié)果的影響。P91鋼穩(wěn)態(tài)蠕變應(yīng)變速率與應(yīng)力的關(guān)系曲線如圖6所示。

    由模擬結(jié)果可以看出,當(dāng)有效跨距l(xiāng)=10 mm時(shí),其整體趨勢(shì)均是隨著尺寸的增大,與單軸拉伸試驗(yàn)的擬合度越高,但誤差均在90%以上;當(dāng)有效跨距l(xiāng)=12 mm時(shí),擬合度更高,整體呈現(xiàn)出隨著尺寸增大,與單軸拉伸試驗(yàn)的擬合度先增大、后減小的趨勢(shì),并且在b=2h=1.1 mm,l=12 mm的尺寸時(shí)達(dá)到最高的擬合度。

    2.2.3 試樣壓頭尺寸的確定

    設(shè)計(jì)不同壓頭直徑與壓頭偏心距進(jìn)行有限元模擬,模擬結(jié)果如圖7,8所示。由圖7可以看出,曲線均較為平穩(wěn),壓頭直徑對(duì)模擬結(jié)果影響不大,所以沿用前人[12]所選取的壓頭直徑2.0 mm。由圖8可以看出,當(dāng)偏心距大于0.3 mm時(shí),曲線開始出現(xiàn)輕微的上下波動(dòng),偏心距對(duì)于模擬結(jié)果的影響開始顯現(xiàn)出來(lái)。為了控制試驗(yàn)誤差,偏心距應(yīng)控制在0.3 mm以內(nèi)。

    圖6 P91鋼穩(wěn)態(tài)蠕變應(yīng)變速率與應(yīng)力σ的關(guān)系曲線

    圖7 壓頭直徑D與穩(wěn)態(tài)位移速率關(guān)系

    圖8 壓頭偏心距e與穩(wěn)態(tài)位移速率關(guān)系

    綜合以上模擬結(jié)果,最終確定模擬所采用的固支直桿三點(diǎn)彎曲試樣的尺寸為:厚度2h=1.1 mm、寬度b=1.1 mm、總長(zhǎng)L=20 mm、有效跨距l(xiāng)=12 mm、壓頭直徑D=2.0 mm。利用此尺寸的固支直桿三點(diǎn)彎曲試樣進(jìn)行后續(xù)P91/P92材料的蠕變性能對(duì)比評(píng)估,可以有效地忽略試樣尺寸對(duì)模擬結(jié)果產(chǎn)生的影響。

    3 P91/P92 鋼蠕變性能對(duì)比分析

    3.1 P91/P92鋼蠕變性能對(duì)比

    利用上述確定的尺寸形式的固支直桿三點(diǎn)彎曲小試樣模型對(duì)P91/P92材料的蠕變性能進(jìn)行模擬對(duì)比。650 ℃下P92鋼的蠕變參數(shù)為:n=7.69,B=2.12×10-19[23]。

    圖9 P91和P92試樣中心位移d與時(shí)間t曲線對(duì)比

    在650 ℃下,當(dāng)載荷為30 N時(shí)對(duì)應(yīng)的兩種材料的中心蠕變位移d和時(shí)間t的曲線如圖9所示??梢钥闯?P91/P92材料中心位移-時(shí)間曲線相差較大,在相同載荷下,P91鋼試樣的中心位移迅速增加,在200 h左右達(dá)到穩(wěn)態(tài)蠕變狀態(tài);P92鋼試樣的中心位移速率遠(yuǎn)遠(yuǎn)小于P91鋼,且基本保持不變。P92鋼作為新型第三代超(超)臨界機(jī)組用鐵素體耐熱鋼,相比于P91鋼來(lái)說(shuō),其高溫持久強(qiáng)度和蠕變性能有了明顯提高。

    3.2 P91/P92鋼蠕變性能對(duì)載荷的敏感度對(duì)比

    兩者在相同應(yīng)力下的穩(wěn)態(tài)蠕變應(yīng)變速率的數(shù)值比較如表5所示。定義應(yīng)變速率增率υ,其計(jì)算公式如下:

    (5)

    圖10 P91和P92鋼穩(wěn)態(tài)蠕變應(yīng)變速率與應(yīng)力σ

    表5 P91/P92鋼穩(wěn)態(tài)蠕變應(yīng)變速率對(duì)比

    由表5可知,在小于400 MPa的應(yīng)力范圍內(nèi),P91鋼的應(yīng)變速率增率大于P92鋼的應(yīng)變速率增率,即P91材料的蠕變性能對(duì)于載荷的變化更加敏感。在相同應(yīng)力水平下,P91鋼的穩(wěn)態(tài)蠕變應(yīng)變速率比P92鋼高出約20~200倍,可見在P91鋼基礎(chǔ)上改進(jìn)后的P92鋼具有更優(yōu)良的蠕變性能。

    3.3 P91/P92鋼蠕變性能對(duì)溫度的敏感度對(duì)比

    查閱相關(guān)文獻(xiàn),得到P91鋼在580,600,620,625,650 ℃和P92鋼在595,610,625,640,650 ℃下的蠕變參數(shù)分別見表6,7[24-28]。

    表6 P91鋼蠕變參數(shù)

    表7 P92 鋼蠕變參數(shù)

    為了驗(yàn)證固支直桿三點(diǎn)彎曲小試樣本構(gòu)模型對(duì)所有溫度工況的適用性,分別對(duì)5個(gè)溫度工況下P91材料的蠕變性能(5個(gè)載荷工況)進(jìn)行了有限元模擬分析,并與對(duì)應(yīng)工況條件下的單軸蠕變?cè)囼?yàn)結(jié)果[24-25]進(jìn)行了對(duì)比,結(jié)果如表8所示??梢钥闯觯晒讨е睏U三點(diǎn)彎曲有限元模擬所得的蠕變參數(shù)n和B與同工況下由單軸蠕變?cè)囼?yàn)得到的結(jié)果較為接近,蠕變指數(shù)n最大偏差小于10%,蠕變常數(shù)B均在相鄰數(shù)量級(jí)。說(shuō)明該模型不受溫度工況條件的影響,且準(zhǔn)確性較高。

    表8 P91鋼單軸蠕變?cè)囼?yàn)與固支直桿三點(diǎn)彎曲有限元模擬所得的蠕變參數(shù)對(duì)比

    在此基礎(chǔ)上,在50 N的載荷下進(jìn)行模擬,對(duì)比P91/P92鋼蠕變性能對(duì)溫度的敏感度。模擬結(jié)果如圖11所示。

    (a)P91鋼

    圖12示出P91/P92鋼材料的穩(wěn)態(tài)位移速率-溫度變化情況??梢钥闯觯S著溫度的升高,P91和P92材料的穩(wěn)態(tài)位移速率均逐漸升高;在相同溫度條件下,P92鋼的穩(wěn)態(tài)蠕變位移速率明顯小于P91鋼,且兩者對(duì)于溫度變化都較敏感。

    為進(jìn)一步對(duì)比考察兩種材料對(duì)溫度的敏感性,定義蠕變速率溫度系數(shù)η(見公式(6))。計(jì)算得P91鋼、P92鋼在溫度低于625 ℃時(shí),蠕變應(yīng)變速率溫度系數(shù)分別為0.005 2 mm/(h·K)和0.003 8 mm/(h·K);在溫度高于625 ℃時(shí),P91鋼、P92鋼的蠕變應(yīng)變速率溫度系數(shù)分別為0.002 1 mm/(h·K)和0.004 5 mm/(h·K)。由此可見,小于625 ℃時(shí),P91材料的蠕變性能對(duì)于溫度更敏感;而高于625 ℃時(shí),P92材料對(duì)溫度更敏感。因此,綜合考慮選材的經(jīng)濟(jì)性和可靠性,在溫度低于625 ℃時(shí)可優(yōu)先選用P92鋼。

    (6)

    圖12 P91/P92鋼穩(wěn)態(tài)位移速率溫度(T)

    4 結(jié)論

    本文基于固支直桿三點(diǎn)彎曲小試樣梁彎曲本構(gòu)模型,采用有限元模擬的方法,確定了用于研究超(超)臨界機(jī)組材料高溫力學(xué)性能的固支直桿三點(diǎn)彎曲小試樣模型尺寸,并基于該尺寸模型對(duì)比分析了P91/P92鋼的蠕變性能。得到主要結(jié)論如下。

    (1)用于研究超(超)臨界機(jī)組材料高溫力學(xué)性能的固支直桿三點(diǎn)彎曲小試樣模型合理尺寸為b=2h=1.1 mm,l=12 mm,其中有效跨距l(xiāng)對(duì)固支直桿三點(diǎn)彎曲小試樣評(píng)價(jià)蠕變性能的影響不大,但當(dāng)寬厚比固定為b∶2h=2∶1和b∶2h=1∶1時(shí),有效跨距的影響較大。

    (2)在650 ℃,30 N載荷下,P92鋼中心位移比P91鋼有大幅下降,在600 h時(shí),P92鋼中心位移僅為P91鋼中心位移的10%,可見P92在高溫下具有更好的抗蠕變性能。

    (3)在100~400 MPa應(yīng)力范圍內(nèi),P91材料的蠕變性能對(duì)于載荷的變化更加敏感,在相同應(yīng)力水平下,P91鋼的穩(wěn)態(tài)蠕變應(yīng)變速率比P92鋼高出約20~200倍,故在高應(yīng)力條件下優(yōu)先選用P92鋼。

    (4)考慮溫度敏感性,在溫度低于625 ℃時(shí),P91材料的蠕變性能對(duì)溫度變化更敏感,可選用P92鋼;在溫度高于625 ℃時(shí),P92鋼的蠕變性能對(duì)溫度變化更敏感,但兩者差距不大,在保證設(shè)備安全性的前提下,考慮到選材的經(jīng)濟(jì)性,可選用P91鋼。

    猜你喜歡
    直桿固支壓頭
    固支結(jié)構(gòu)對(duì)壓阻式壓力傳感器的輸出影響研究
    基于ANSYS Workbench 的雜質(zhì)資源化處理設(shè)備壓頭靜力學(xué)及固有頻率分析*
    顫振環(huán)境軟金屬碰撞滑動(dòng)接觸摩擦的分子動(dòng)力學(xué)模擬
    一種可拆卸式厚板壓平機(jī)輔助壓頭
    中國(guó)科學(xué)技術(shù)館之“雙曲隧道”
    基于數(shù)據(jù)背景下采用“模擬退火算法”確定太陽(yáng)影子定位
    基于復(fù)合材料的山體滑坡固支模型及配套傳感器布置方法申請(qǐng)?zhí)? 202010514582.0
    傳感器世界(2021年2期)2021-03-27 19:36:46
    一角點(diǎn)支撐對(duì)面兩邊固支正交各向異性矩形薄板彎曲問(wèn)題的辛疊加解
    你的身高我知道
    直桿與圓環(huán)組合體落地反彈問(wèn)題的解法探討
    物理之友(2018年7期)2018-08-06 00:51:46
    午夜激情欧美在线| 国产免费av片在线观看野外av| 亚洲专区国产一区二区| 在线十欧美十亚洲十日本专区| 一个人看的www免费观看视频| 变态另类成人亚洲欧美熟女| 精品国产超薄肉色丝袜足j| 两个人的视频大全免费| 搡老岳熟女国产| 免费av毛片视频| 制服丝袜大香蕉在线| 18禁黄网站禁片午夜丰满| 亚洲第一电影网av| 一二三四在线观看免费中文在| 欧美av亚洲av综合av国产av| а√天堂www在线а√下载| 国产日本99.免费观看| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址| 日本精品一区二区三区蜜桃| 欧美成狂野欧美在线观看| 美女扒开内裤让男人捅视频| 久久精品亚洲精品国产色婷小说| 一进一出抽搐gif免费好疼| 伊人久久大香线蕉亚洲五| 97超视频在线观看视频| 成人亚洲精品av一区二区| 欧美午夜高清在线| 国产亚洲av嫩草精品影院| a在线观看视频网站| 免费看日本二区| 国产人伦9x9x在线观看| 淫妇啪啪啪对白视频| 天堂影院成人在线观看| 在线国产一区二区在线| 香蕉久久夜色| 国产亚洲精品久久久久久毛片| 精品电影一区二区在线| 久久欧美精品欧美久久欧美| 丰满的人妻完整版| 在线永久观看黄色视频| 亚洲精品国产精品久久久不卡| 最近最新中文字幕大全免费视频| 99久久国产精品久久久| 怎么达到女性高潮| 亚洲av片天天在线观看| 色哟哟哟哟哟哟| 久久久久久九九精品二区国产| 国产精品av视频在线免费观看| 丰满人妻一区二区三区视频av | 久久久久国产一级毛片高清牌| 蜜桃久久精品国产亚洲av| 中文字幕最新亚洲高清| 一区二区三区国产精品乱码| 怎么达到女性高潮| 成人精品一区二区免费| 丁香六月欧美| 国产精品野战在线观看| 欧美精品啪啪一区二区三区| 亚洲精品在线美女| 精华霜和精华液先用哪个| 国产成年人精品一区二区| 91av网站免费观看| 亚洲成av人片免费观看| 欧美成狂野欧美在线观看| 国产成人一区二区三区免费视频网站| 亚洲国产色片| 桃红色精品国产亚洲av| 欧美日韩一级在线毛片| 欧美中文日本在线观看视频| 欧美另类亚洲清纯唯美| 国产1区2区3区精品| 欧美乱色亚洲激情| 亚洲精品在线美女| 人人妻人人澡欧美一区二区| 99国产精品一区二区蜜桃av| 国产综合懂色| 成人av一区二区三区在线看| 亚洲九九香蕉| 国产视频内射| 欧美日韩精品网址| 久久天堂一区二区三区四区| 黄色片一级片一级黄色片| 真人做人爱边吃奶动态| 日韩成人在线观看一区二区三区| 国产成人精品无人区| 夜夜夜夜夜久久久久| 久久精品国产综合久久久| 国产三级中文精品| av国产免费在线观看| 久久精品影院6| 国产精品久久久久久人妻精品电影| 免费观看的影片在线观看| 精品久久蜜臀av无| 中文字幕高清在线视频| 国内精品美女久久久久久| 精品国产三级普通话版| 在线观看美女被高潮喷水网站 | 又黄又粗又硬又大视频| 欧美av亚洲av综合av国产av| 黄片大片在线免费观看| 国产探花在线观看一区二区| 午夜福利18| 51午夜福利影视在线观看| 老司机福利观看| 亚洲欧美精品综合一区二区三区| 天堂√8在线中文| 熟女人妻精品中文字幕| 久久久久国产精品人妻aⅴ院| 国产不卡一卡二| 欧美在线一区亚洲| 日本一本二区三区精品| 999精品在线视频| 国产av不卡久久| 一进一出抽搐动态| 少妇熟女aⅴ在线视频| 日韩欧美一区二区三区在线观看| 99在线视频只有这里精品首页| 成人欧美大片| 国产男靠女视频免费网站| 久久亚洲精品不卡| 亚洲激情在线av| 精品人妻1区二区| 一a级毛片在线观看| 99re在线观看精品视频| 国产私拍福利视频在线观看| 最好的美女福利视频网| 亚洲av成人av| 香蕉丝袜av| 99riav亚洲国产免费| 成人三级做爰电影| 看黄色毛片网站| 亚洲国产精品合色在线| 国内精品美女久久久久久| 一级黄色大片毛片| 哪里可以看免费的av片| e午夜精品久久久久久久| 日本熟妇午夜| 国产高清视频在线观看网站| 首页视频小说图片口味搜索| 久久国产精品影院| 国产一区二区三区在线臀色熟女| 国产久久久一区二区三区| 国产乱人伦免费视频| 国内揄拍国产精品人妻在线| 亚洲欧美日韩卡通动漫| 制服人妻中文乱码| 好男人在线观看高清免费视频| 一级作爱视频免费观看| 欧美成人免费av一区二区三区| 桃色一区二区三区在线观看| 国产在线精品亚洲第一网站| 亚洲最大成人中文| 黑人操中国人逼视频| 小说图片视频综合网站| 国产精品综合久久久久久久免费| 热99在线观看视频| 夜夜爽天天搞| 午夜精品在线福利| 久9热在线精品视频| 国产aⅴ精品一区二区三区波| 国产午夜精品论理片| 午夜福利在线观看免费完整高清在 | 国产真实乱freesex| 国产精品永久免费网站| 好看av亚洲va欧美ⅴa在| 黄色 视频免费看| 香蕉av资源在线| 69av精品久久久久久| 成人鲁丝片一二三区免费| 9191精品国产免费久久| av天堂中文字幕网| 亚洲 欧美 日韩 在线 免费| 婷婷精品国产亚洲av在线| 成年女人永久免费观看视频| av视频在线观看入口| 久99久视频精品免费| 一本精品99久久精品77| 搡老熟女国产l中国老女人| 成人特级黄色片久久久久久久| 国产又色又爽无遮挡免费看| а√天堂www在线а√下载| 国内精品久久久久久久电影| 怎么达到女性高潮| 国内精品一区二区在线观看| 亚洲av成人不卡在线观看播放网| 久久久水蜜桃国产精品网| 亚洲成人久久爱视频| 亚洲无线观看免费| 日韩欧美国产一区二区入口| 久久久久久大精品| 后天国语完整版免费观看| 午夜久久久久精精品| 亚洲av五月六月丁香网| 少妇丰满av| 免费大片18禁| 性色avwww在线观看| 一卡2卡三卡四卡精品乱码亚洲| 一本综合久久免费| 91字幕亚洲| 女生性感内裤真人,穿戴方法视频| 久久九九热精品免费| 亚洲国产精品久久男人天堂| 国产一区二区在线av高清观看| 长腿黑丝高跟| 又大又爽又粗| 99久久成人亚洲精品观看| 免费av不卡在线播放| 久久99热这里只有精品18| 国产精品久久久久久人妻精品电影| 熟女少妇亚洲综合色aaa.| 成人国产综合亚洲| 久久这里只有精品中国| 精品久久久久久久末码| 国产人伦9x9x在线观看| 成年女人看的毛片在线观看| 日本a在线网址| 制服人妻中文乱码| 不卡一级毛片| 国产av在哪里看| 日韩欧美三级三区| 一个人观看的视频www高清免费观看 | 免费观看的影片在线观看| 国产又黄又爽又无遮挡在线| 小说图片视频综合网站| 91在线精品国自产拍蜜月 | 国产一级毛片七仙女欲春2| 我要搜黄色片| 国产精品野战在线观看| 亚洲av免费在线观看| 两人在一起打扑克的视频| 九九在线视频观看精品| 美女高潮喷水抽搐中文字幕| 欧美色视频一区免费| 国产欧美日韩精品一区二区| 99在线视频只有这里精品首页| 亚洲国产看品久久| 97超视频在线观看视频| a在线观看视频网站| 脱女人内裤的视频| 2021天堂中文幕一二区在线观| 偷拍熟女少妇极品色| 亚洲男人的天堂狠狠| netflix在线观看网站| 色哟哟哟哟哟哟| 99久久精品热视频| 淫秽高清视频在线观看| 亚洲人成电影免费在线| 日日夜夜操网爽| 国产爱豆传媒在线观看| 在线永久观看黄色视频| 午夜福利视频1000在线观看| 色视频www国产| 日韩欧美在线乱码| 亚洲 国产 在线| 久久精品综合一区二区三区| 熟女人妻精品中文字幕| 9191精品国产免费久久| 人人妻人人澡欧美一区二区| 制服丝袜大香蕉在线| 欧美一级毛片孕妇| 激情在线观看视频在线高清| 久久99热这里只有精品18| 桃红色精品国产亚洲av| 中文字幕久久专区| 日韩人妻高清精品专区| 亚洲九九香蕉| 欧美一级毛片孕妇| 白带黄色成豆腐渣| 88av欧美| av女优亚洲男人天堂 | 国产野战对白在线观看| 国产美女午夜福利| 日韩国内少妇激情av| 一进一出好大好爽视频| aaaaa片日本免费| 亚洲午夜精品一区,二区,三区| 亚洲精品国产精品久久久不卡| 久久久久精品国产欧美久久久| 亚洲国产精品sss在线观看| 久久久国产成人免费| 国内精品一区二区在线观看| 观看美女的网站| 美女被艹到高潮喷水动态| 日韩欧美一区二区三区在线观看| 成人亚洲精品av一区二区| 女人高潮潮喷娇喘18禁视频| 久久九九热精品免费| 欧美不卡视频在线免费观看| 国产精品自产拍在线观看55亚洲| 人妻丰满熟妇av一区二区三区| 一区福利在线观看| 欧美成人免费av一区二区三区| 此物有八面人人有两片| 成人av在线播放网站| 欧美日韩精品网址| 亚洲av熟女| 9191精品国产免费久久| 欧美性猛交黑人性爽| 夜夜躁狠狠躁天天躁| 国内揄拍国产精品人妻在线| 精品99又大又爽又粗少妇毛片 | 成人鲁丝片一二三区免费| 午夜视频精品福利| 久久久色成人| 麻豆久久精品国产亚洲av| av在线天堂中文字幕| 国产精品电影一区二区三区| 99在线人妻在线中文字幕| 天堂影院成人在线观看| 国产私拍福利视频在线观看| 成熟少妇高潮喷水视频| 色噜噜av男人的天堂激情| 亚洲成a人片在线一区二区| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| 午夜两性在线视频| 美女cb高潮喷水在线观看 | 一级作爱视频免费观看| 人妻夜夜爽99麻豆av| 俺也久久电影网| 亚洲av五月六月丁香网| 欧美黑人巨大hd| 两性午夜刺激爽爽歪歪视频在线观看| 性色av乱码一区二区三区2| 国产伦精品一区二区三区四那| 国产精品免费一区二区三区在线| 一本精品99久久精品77| 久久精品aⅴ一区二区三区四区| 亚洲av美国av| 人妻久久中文字幕网| 男女午夜视频在线观看| av片东京热男人的天堂| 亚洲自拍偷在线| 亚洲人与动物交配视频| 日韩欧美国产一区二区入口| 热99re8久久精品国产| 1000部很黄的大片| 国产淫片久久久久久久久 | 中文字幕久久专区| 天堂√8在线中文| 九色国产91popny在线| 欧美午夜高清在线| svipshipincom国产片| 一区福利在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久人人做人人爽| 精品日产1卡2卡| 人人妻人人澡欧美一区二区| 美女高潮的动态| 人人妻,人人澡人人爽秒播| 三级毛片av免费| 国产精品综合久久久久久久免费| 黄色片一级片一级黄色片| 12—13女人毛片做爰片一| 给我免费播放毛片高清在线观看| 动漫黄色视频在线观看| 日本一二三区视频观看| 黄色 视频免费看| 国产亚洲精品久久久com| 亚洲午夜精品一区,二区,三区| 18禁美女被吸乳视频| 国产午夜精品论理片| 舔av片在线| 91在线观看av| 少妇裸体淫交视频免费看高清| 国产亚洲精品av在线| 亚洲国产欧美人成| 亚洲国产日韩欧美精品在线观看 | 99热6这里只有精品| 18禁观看日本| 久久久久性生活片| 亚洲自拍偷在线| 国产成人欧美在线观看| 啦啦啦免费观看视频1| 欧美成人免费av一区二区三区| 亚洲国产精品999在线| 国产v大片淫在线免费观看| 亚洲国产精品sss在线观看| 精品久久蜜臀av无| 少妇丰满av| 99久国产av精品| 999久久久国产精品视频| 亚洲国产色片| 又粗又爽又猛毛片免费看| 亚洲成人久久爱视频| 日韩免费av在线播放| 精品熟女少妇八av免费久了| www.自偷自拍.com| 成年免费大片在线观看| 国产成人系列免费观看| 欧美一级毛片孕妇| 日韩欧美一区二区三区在线观看| 国产综合懂色| 国产成人福利小说| 精品久久久久久久人妻蜜臀av| 久久久久久国产a免费观看| 亚洲av中文字字幕乱码综合| 变态另类丝袜制服| 亚洲av中文字字幕乱码综合| 国产精品免费一区二区三区在线| 脱女人内裤的视频| 90打野战视频偷拍视频| 不卡av一区二区三区| 制服丝袜大香蕉在线| 亚洲国产精品合色在线| 啪啪无遮挡十八禁网站| 日韩精品中文字幕看吧| 亚洲七黄色美女视频| 91在线观看av| 亚洲av第一区精品v没综合| 欧美激情久久久久久爽电影| 国产精品久久视频播放| 老汉色∧v一级毛片| 精品欧美国产一区二区三| 一区二区三区高清视频在线| 久久中文字幕人妻熟女| 亚洲成av人片在线播放无| 免费观看人在逋| 日本熟妇午夜| 99riav亚洲国产免费| 亚洲精品在线美女| 精品一区二区三区四区五区乱码| 欧美av亚洲av综合av国产av| 亚洲第一电影网av| 国产精品一区二区三区四区久久| 国产激情偷乱视频一区二区| 国产一区二区三区在线臀色熟女| 欧美激情久久久久久爽电影| 在线观看午夜福利视频| 国产精品综合久久久久久久免费| 99久久精品国产亚洲精品| 久久久国产欧美日韩av| 蜜桃久久精品国产亚洲av| 丰满的人妻完整版| 99精品在免费线老司机午夜| 欧美性猛交╳xxx乱大交人| 欧美精品啪啪一区二区三区| 国模一区二区三区四区视频 | 美女高潮的动态| 国产真人三级小视频在线观看| 午夜成年电影在线免费观看| 亚洲熟妇熟女久久| 日韩欧美国产在线观看| 哪里可以看免费的av片| 国产伦一二天堂av在线观看| 成人午夜高清在线视频| 久久精品91蜜桃| 99精品久久久久人妻精品| www.自偷自拍.com| 99久国产av精品| 久久婷婷人人爽人人干人人爱| 精品久久久久久久久久免费视频| 亚洲av日韩精品久久久久久密| 成人永久免费在线观看视频| cao死你这个sao货| 色哟哟哟哟哟哟| 美女 人体艺术 gogo| 久久久国产欧美日韩av| 国产精品 欧美亚洲| 日本精品一区二区三区蜜桃| 亚洲av片天天在线观看| 国产成人精品久久二区二区91| 黄色成人免费大全| 精品一区二区三区视频在线 | 嫩草影院精品99| 少妇熟女aⅴ在线视频| 中国美女看黄片| 午夜视频精品福利| 又爽又黄无遮挡网站| 国内揄拍国产精品人妻在线| 免费观看精品视频网站| 欧美性猛交╳xxx乱大交人| 欧美日韩精品网址| 1024手机看黄色片| 欧美日本亚洲视频在线播放| 午夜成年电影在线免费观看| 国产午夜精品久久久久久| 亚洲一区二区三区不卡视频| 丝袜人妻中文字幕| 好男人电影高清在线观看| 特大巨黑吊av在线直播| 欧美午夜高清在线| 真人做人爱边吃奶动态| 成人午夜高清在线视频| 我的老师免费观看完整版| 99久久精品国产亚洲精品| 国产精品av视频在线免费观看| 亚洲在线自拍视频| 精品久久久久久久毛片微露脸| 欧美日韩瑟瑟在线播放| 日本一二三区视频观看| 99国产极品粉嫩在线观看| 欧洲精品卡2卡3卡4卡5卡区| АⅤ资源中文在线天堂| 国产免费男女视频| 成人无遮挡网站| 欧美一区二区精品小视频在线| 97超级碰碰碰精品色视频在线观看| 午夜a级毛片| 精品久久久久久久久久免费视频| 老汉色av国产亚洲站长工具| 午夜日韩欧美国产| 国产精品 国内视频| 欧美一区二区精品小视频在线| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区| 制服人妻中文乱码| av在线天堂中文字幕| 一级毛片高清免费大全| 亚洲色图av天堂| 国产亚洲av嫩草精品影院| 亚洲精品一区av在线观看| 美女扒开内裤让男人捅视频| 亚洲精品粉嫩美女一区| 欧美绝顶高潮抽搐喷水| 男插女下体视频免费在线播放| 97碰自拍视频| 午夜福利成人在线免费观看| 久久婷婷人人爽人人干人人爱| 法律面前人人平等表现在哪些方面| 神马国产精品三级电影在线观看| 一级a爱片免费观看的视频| 欧美丝袜亚洲另类 | 亚洲欧美精品综合一区二区三区| 亚洲五月天丁香| 熟妇人妻久久中文字幕3abv| 99精品欧美一区二区三区四区| 高潮久久久久久久久久久不卡| 精品乱码久久久久久99久播| 熟女少妇亚洲综合色aaa.| 午夜福利18| 性欧美人与动物交配| 男女午夜视频在线观看| 99久久精品国产亚洲精品| 99国产极品粉嫩在线观看| 美女高潮的动态| 一级a爱片免费观看的视频| 久久久久国产一级毛片高清牌| 男女下面进入的视频免费午夜| av片东京热男人的天堂| 国产三级在线视频| 国产v大片淫在线免费观看| av国产免费在线观看| 蜜桃久久精品国产亚洲av| 国产精品av视频在线免费观看| 国内久久婷婷六月综合欲色啪| 国产精品,欧美在线| 欧美大码av| 国产精品一区二区三区四区久久| 制服人妻中文乱码| 久久这里只有精品中国| 亚洲人与动物交配视频| 国产高潮美女av| 无人区码免费观看不卡| 国产男靠女视频免费网站| 中出人妻视频一区二区| 久久久成人免费电影| 亚洲av成人精品一区久久| 琪琪午夜伦伦电影理论片6080| 欧美日韩国产亚洲二区| 日韩高清综合在线| 女人高潮潮喷娇喘18禁视频| 最新美女视频免费是黄的| 欧美日韩乱码在线| 色哟哟哟哟哟哟| 人妻久久中文字幕网| 久久久久久久午夜电影| av天堂中文字幕网| 欧美3d第一页| 亚洲,欧美精品.| 国产免费男女视频| 给我免费播放毛片高清在线观看| 国产野战对白在线观看| 亚洲在线自拍视频| 免费人成视频x8x8入口观看| 成人国产一区最新在线观看| 高清在线国产一区| 国产精品女同一区二区软件 | 免费看日本二区| 88av欧美| 岛国在线观看网站| 久久这里只有精品19| 国产精品一及| www日本在线高清视频| 成人性生交大片免费视频hd| 国产精品久久久久久精品电影| 亚洲国产日韩欧美精品在线观看 | 亚洲 欧美一区二区三区| 最近最新中文字幕大全免费视频| 伊人久久大香线蕉亚洲五| 黄色丝袜av网址大全| 国产伦人伦偷精品视频| 成人av一区二区三区在线看| 国产精品永久免费网站| 99国产综合亚洲精品| 不卡一级毛片| 成人亚洲精品av一区二区| 亚洲av电影在线进入| 国产一级毛片七仙女欲春2| 欧美日韩亚洲国产一区二区在线观看| 午夜福利免费观看在线| 亚洲欧美精品综合一区二区三区| 欧美3d第一页| 国产黄a三级三级三级人| 国产精品亚洲一级av第二区| 色视频www国产| 国产午夜精品久久久久久| 久久国产精品影院| 国产v大片淫在线免费观看| 19禁男女啪啪无遮挡网站| 99视频精品全部免费 在线 | 精品乱码久久久久久99久播| 国产精品国产高清国产av| 18禁黄网站禁片午夜丰满|