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

    缺陷/應力交互對碳鋼Lcr波聲彈性系數(shù)的影響

    2017-07-25 11:20:11石常亮繆文炳董世運
    材料工程 2017年7期
    關鍵詞:盲孔時間差傳播速度

    劉 彬,石常亮,繆文炳,董世運

    (1 江蘇科技大學 材料科學與工程學院,江蘇 鎮(zhèn)江 212003; 2 廣東省工業(yè)分析檢測中心,廣州 510650;3 裝甲兵工程學院裝備再制造技術國防科技重點實驗室,北京 100072)

    缺陷/應力交互對碳鋼Lcr波聲彈性系數(shù)的影響

    劉 彬1,石常亮2,繆文炳1,董世運3

    (1 江蘇科技大學 材料科學與工程學院,江蘇 鎮(zhèn)江 212003; 2 廣東省工業(yè)分析檢測中心,廣州 510650;3 裝甲兵工程學院裝備再制造技術國防科技重點實驗室,北京 100072)

    基于Lcr波聲彈性理論,探討缺陷及其尺寸對Lcr波評價應力的影響機理。結合“當量法”預制不同直徑盲孔,采用互相關系數(shù)函數(shù)計算Lcr波時間差,通過線性擬合得到Lcr波聲彈性系數(shù),基于彈塑性變形和圓孔應力集中理論澄清盲孔直徑對Lcr波聲彈性系數(shù)的影響機理。結果表明:各直徑盲孔Lcr波時間差隨應力增大基本呈線性增加,但其非線性特征亦逐漸明顯,線性階段的最大應力值小于試樣屈服強度;Lcr波聲彈性系數(shù)隨盲孔直徑增大逐漸減小,并趨于平穩(wěn)。分析認為,盲孔應力集中是導致上述結果的主要原因,試樣各向異性組織及盲孔深度也是其重要因素。

    Lcr波聲彈性理論;應力集中;無損評價;彈塑性變形;交互影響

    應力(殘余應力或工作應力)是影響甚至決定產品性能與服役安全的重要因素,因而進行應力評價對保證產品服役安全與可靠性就顯得極為重要。目前,應力評價方法分為無損評價與有損評價兩類。有損評價法[1,2]雖可實現(xiàn)應力的評價,但是在破壞產品完整性基礎上進行,且評價結果為抽樣檢測結果,因而給產品安全帶來隱患。無損評價法,如超聲波法、光學法等[3-6],是在保證產品完整性基礎上實現(xiàn)應力的評價,可為產品服役安全與可靠性的在線評價和監(jiān)測提供技術支撐。其中,超聲波具有方便、安全、適于現(xiàn)場檢測及可實現(xiàn)在線監(jiān)測等優(yōu)點,因此在應力評價領域引起眾多學者的廣泛關注。

    超聲波法是以超聲波聲彈性理論為基礎而實現(xiàn)應力評價的一種無損評價方法,目前在石油管道、航空航天、壓力容器等領域得到推廣應用,但其評價機理仍待澄清。為此,國內外學者開展了系列相關實驗研究與理論探討,國外Yashar等[7,8],Seyedali等[9]對Lcr波(臨界折射縱波)評價焊接殘余應力開展了系列研究,并與數(shù)值模擬結果及小孔法結果進行了對比分析。國內也有學者采用不同波形超聲波對焊接殘余應力進行了評價[10-13]。Hu等[14]采用超聲波對碳鋼結構應力進行了評價,并探討了影響評價結果的因素。魏智等[15]基于材料本征參數(shù)測量對超聲表面波在低碳鋼與鋁合金中的傳播速度和應力間關系進行了理論研究,并與實驗結果進行了對比分析。上述研究雖取得了一定成果,但均未考慮缺陷(尤其是內部不可見缺陷)對超聲波評價應力結果的影響,而如焊接缺陷(氣孔、裂紋、夾雜)以及與凝固相關的缺陷等均是難以避免的,且相關研究表明缺陷不僅會改變應力的分布狀態(tài),而且會影響應力的數(shù)值;因而,探討缺陷及其尺寸對超聲波評價應力結果的影響與機理對超聲波技術的推廣以及結果保證就顯得極為重要。

    為避免其他因素的影響,本工作選取Lcr波對低碳鋼應力進行評價,基于“當量法”,以深度恒定與不同直徑的規(guī)則盲孔模擬內部缺陷,利用互相關系數(shù)函數(shù)計算Lcr波時間差,基于Lcr波聲彈性理論,提取并建立Lcr波聲彈性系數(shù)與盲孔直徑間關系,結合彈塑性變形及圓孔應力集中理論,澄清盲孔直徑對Lcr波聲彈性系數(shù)的影響機理。

    1 實驗材料與方法

    實驗材料選擇綜合性能良好的Q235鋼,其力學性能如表1所示。

    表1 Q235鋼力學性能Table 1 Mechanical properties of Q235 steel

    依據(jù)本工作Lcr波采集方式和實驗要求加工預制盲孔缺陷的非標準室溫靜載拉伸試樣,如圖1所示。結合“當量法”,采用電火花打孔方法沿試樣寬度方向中心線每間隔30mm加工深度為3.0mm,直徑分別為0.5,1.0,1.5,2.0,2.5,3.0mm的盲孔。前期研究表明,由機加工等導致的試樣初始應力是影響超聲波評價應力結果的重要因素,為此采用WZS-20雙室真空燒結爐對Q235鋼試樣進行去應力退火處理,加熱溫度為(600±10)℃,保溫時間為30min,爐冷至100℃,再空冷至室溫。

    圖1 預置盲孔靜載拉伸試樣示意圖Fig.1 Schematic diagram of preset blind-hole static tensile specimen

    超聲波聲彈性效應是一種弱效應,一般而言應力改變100MPa,超聲波在鋼中的傳播速度改變量約為0.01%[16],因而超聲波傳播速度的精確測量對超聲波評價應力結果至關重要。鑒于此,探討可替代Lcr波傳播速度的應力評價特征參數(shù)的計算方法就亟待解決?;诖?,本工作固定Lcr波在試樣中的傳播距離恒定不變,以Lcr波傳播時間差(Lcr波時間差)表征Lcr波傳播速度改變量。目前,超聲波傳播時間差計算方法眾多,基于簡單、方便等優(yōu)點,互相關函數(shù)法在該領域得到廣泛應用。結合本工作Lcr波信號特點,對互相關函數(shù)進行歸一化離散處理,得到適用于離散數(shù)據(jù)分析的歸一化離散互相關函數(shù),并定義為互相關系數(shù)函數(shù),可表示為:

    (1)

    式中:x(i)和y(i)為不同應力時Q235鋼Lcr波信號;n為計算點數(shù)(本工作定義互相關步長)。

    基于超聲波聲彈性理論和各向同性材料假設,Lcr波在薄板中的傳播速度與應力間關系可表示為[17]:

    (2)

    式中:V1和V0分別為Lcr波在應力與無應力介質中的傳播速度;σ11和σ22分別為平行于加載方向和垂直于加載方向的應力;k1和k2分別為平行于加載方向和垂直于加載方向的聲彈性系數(shù)。

    相關研究表明,薄板中Lcr波聲彈性系數(shù)滿足k1?k2,即式(2)可簡化為:

    (3)

    結合本工作研究方法可知,Lcr波在Q235鋼中的傳播距離恒定不變,因此式(3)可簡化為:

    Δt=t0·k1·σ11=k·σ

    (4)

    式中:k為Lcr波聲彈性系數(shù);Δt為Lcr波時間差。

    鑒于此,搭建Lcr波應力評價系統(tǒng),主要包括超聲波脈沖發(fā)射接收儀5072PR、熒光數(shù)字示波器DPO3034、自主研發(fā)的中心頻率為2.5MHz的雙Lcr波換能器和自主開發(fā)的基于Matlab的Lcr波時間差計算系統(tǒng)。

    2 結果與分析

    2.1 分析討論

    將預制盲孔缺陷的非標準室溫靜載拉伸Q235鋼試樣豎直夾持于SANS-CMT5205電子萬能試驗機,以0.5kN/s緩慢加載至預定載荷不卸載,將Lcr波換能器沿平行于加載方向放置于盲孔正上方的檢測面采集Lcr波,如圖1所示,直至屈服載荷停止加載。為避免Lcr波換能器與Q235鋼間耦合對應力評價結果的影響,采用自制Lcr波換能器夾具保持換能器與Q235鋼間耦合壓力恒定不變,并對各載荷Q235鋼Lcr波重復采集3次。對比可知,載荷相同時,Q235鋼Lcr波基本重合。鑒于此,以第1次采集的Q235鋼Lcr波為分析對象,探討不同直徑盲孔時Q235鋼Lcr波沿時間軸隨應力的變化規(guī)律。不同應力時,直徑分別為0.5,3.0mm盲孔處Q235鋼Lcr波(接收信號)如圖2所示。

    由圖2可知,隨應力增大,直徑為0.5,3.0mm盲孔處Q235鋼Lcr波沿時間軸的變化規(guī)律基本相同,即隨應力增大,盲孔處Lcr波沿時間軸逐漸向右移動,這表明Lcr波通過相同距離所需傳播時間逐漸延長,因而Lcr波在Q235鋼中的傳播速度隨應力增大逐漸減小,這與Lcr波聲彈性理論是相符的;當應力分別達到200MPa和230MPa(小于試樣屈服強度)時,再隨應力增大,盲孔處Lcr波沿時間軸呈不規(guī)律變化,即Lcr波傳播速度隨應力增大呈不規(guī)律變化,這與Lcr波聲彈性理論略有不符。分析認為,盲孔處Q235鋼的變形狀態(tài)是決定Lcr波在其中的傳播速度隨應力變化的關鍵,盲孔應力集中效應引起的局部區(qū)域應力遠大于試樣平均應力是導致上述結果的主要原因。

    圖2 不同直徑盲孔處Q235鋼Lcr波 (a)0.5mm;(b)3.0mmFig.2 Lcr wave of Q235 steel corresponding to different diameters of blind-hole (a)0.5mm;(b)3.0mm

    為獲得Q235鋼Lcr波聲彈性系數(shù),通過式(1)計算各應力時Lcr波時間差。前期相關研究[18]表明,互相關步長是影響式(1)計算結果的重要因素,隨互相關步長增大,互相關系數(shù)函數(shù)計算結果精度越高,Lcr波時間差計算結果精度亦越高。鑒于此,本工作選擇1周期長度為互相關步長計算盲孔處Q235鋼Lcr波時間差,并建立Lcr波時間差與應力間關系,如圖3所示。

    圖3 不同直徑盲孔處Lcr波時間差與應力曲線Fig.3 Curves of the difference in time of flight between Lcr wave and stress corresponding to different diameters of blind-hole

    圖3表明,盲孔直徑不同時,盲孔處Lcr波時間差隨應力變化的規(guī)律基本相同。低應力階段,盲孔處Lcr波時間差基本呈線性規(guī)律隨應力增大而增加,與理論變化規(guī)律差別不明顯;高應力階段,盲孔處Lcr波時間差隨應力增大而增加,但Lcr波時間差波動較明顯,且隨應力增大,實驗規(guī)律與理論規(guī)律間差別亦越明顯;當應力分別達到200,208,215,220,225,230MPa時,盲孔處Lcr波時間差隨應力增大呈非線性規(guī)律變化。對比相同應力時盲孔處Lcr波時間差可知,隨盲孔直徑增大,Lcr波時間差逐漸減小,且應力越大,Lcr波時間差的變化也越大。為驗證上述實驗結果,采用ANSYS有限元軟件對預制盲孔的Q235鋼應力進行數(shù)值模擬計算,其應力云圖及沿盲孔中心的應力曲線如圖4和圖5所示。

    圖4 不同載荷時Q235鋼應力分布 (a)15kN;(b)25kN;(c)50kN;(d)62.5kNFig.4 Stress distributions of Q235 steel with different loads (a)15kN;(b)25kN;(c)50kN;(d)62.5kN

    圖5 沿盲孔中心的應力曲線Fig.5 Stress curves along the center of blind-hole

    圖4給出了載荷分別為15,25,50,62.5kN時,Q235鋼的應力云圖。由圖4可知,盲孔處應力值遠高于試樣平均應力,且盲孔處應力呈明顯啞鈴狀分布,這與無限大薄板圓孔應力集中理論是相符的。圖5給出了載荷分別為15,25,50,62.5kN時,沿盲孔中心分布的應力曲線。由于本工作Lcr波時間差是其傳播距離內平均應力的反映,因而對盲孔處應力平均值進行分析。對比可知,載荷相同時,盲孔處應力平均值隨盲孔直徑增大逐漸減??;盲孔直徑相同時,與載荷由15kN增加至25kN相比,載荷由50kN增加至62.5kN時盲孔處的應力變化更為明顯。這與圖3中不同應力階段時盲孔處Lcr波時間差隨應力變化的規(guī)律是相符的。

    基于公式(4),采用線性函數(shù)對線性階段Q235鋼Lcr波時間差與應力進行擬合,得到其聲彈性系數(shù)k,建立k與盲孔直徑間關系,為獲得盲孔直徑對Lcr波聲彈性系數(shù)影響的分析,采用冪函數(shù)對其進行擬合,結果如式(5)和圖6所示。

    圖6 Lcr波聲彈性系數(shù)與盲孔直徑關系Fig.6 Relationship between Lcr wave acoustoelastic coefficient and diameter of blind-holek=0.2154·D-0.2327

    (5)

    式中:k為Lcr波聲彈性系數(shù);D為盲孔直徑。

    由圖6可知,隨盲孔直徑增大,Q235鋼盲孔處Lcr波聲彈性系數(shù)k逐漸減小,并趨于平穩(wěn),這與文獻[19]中數(shù)值模擬結果的變化規(guī)律是相符的;因而,通過式(5)即可獲得任意直徑盲孔時Q235鋼Lcr波聲彈性系數(shù),這也就表明采用本工作方法可實現(xiàn)缺陷/應力交互作用對Lcr波聲彈性系數(shù)影響的研究。

    2.2 結果分析

    介質彈性變形是Lcr波評價應力的前提,也是Lcr波聲彈性效應的理論基礎,該理論表明介質處于彈性變形狀態(tài)時,Lcr波在介質中的傳播速度與應力呈嚴格線性關系,這與圖3中Q235鋼盲孔處Lcr波時間差隨應力的變化規(guī)律基本相符,但略有不同。分析認為,盲孔引起的Q235鋼局部應力集中是導致上述結果的主要原因。應力集中理論表明,受單向均勻載荷作用的無限大薄板中圓孔會導致其周圍應力分布發(fā)生改變,即圓孔邊緣出現(xiàn)應力峰值,且應力呈非線性規(guī)律隨距圓孔距離增大而減小,直至達到均勻應力值,如圖7所示。

    圖7 無限大薄板中圓孔應力集中示意圖Fig.7 Schematic diagram of stress concentration caused by hole in infinite plate

    結合本工作Lcr波信號采集方式可知,圖2中Lcr波沿時間軸的變化規(guī)律是其傳播距離內應力平均值的反映,而盲孔應力集中不僅改變了應力的分布狀態(tài),而且影響應力的數(shù)值,因此圖3中實驗規(guī)律與理論規(guī)律并不完全相同。對圖3中不同應力階段Lcr波時間差與應力關系分析認為,低應力階段時,由盲孔引起的應力集中程度較小,如25kN時Q235鋼應力云圖與應力曲線所示,盲孔周圍最高應力值約為83MPa,遠小于Q235鋼屈服強度,因而試樣處于彈性變形狀態(tài),滿足Lcr波聲彈性理論要求,Lcr波時間差與應力間線性關系良好;高應力階段時,由盲孔引起的應力集中程度逐漸變大,試樣整體雖仍處于彈性變形狀態(tài),但其最高應力值處Q235鋼已進入塑性變形狀態(tài),即Lcr波傳播距離內試樣變形狀態(tài)不再滿足Lcr波聲彈性理論要求。結合本工作Lcr波評價應力方法可知,Lcr波傳播速度的變化是其傳播距離內Q235鋼變形狀態(tài)的“平均反映”,而應力越大,早于試樣整體進入塑性變形狀態(tài)的區(qū)域越大,即塑性變形面積越大,因而高應力階段Lcr波時間差與應力間關系具有明顯非線性特征,且隨應力增大,該非線性特征越明顯;當Lcr波傳播距離內盲孔應力集中引起的試樣塑性變形對Lcr波傳播速度的影響起主導作用時,Q235鋼變形狀態(tài)不再滿足Lcr波聲彈性理論要求,因此Lcr波時間差與應力不再呈線性規(guī)律變化,這也是應力小于Q235鋼屈服強度時盲孔處Lcr波時間差隨應力呈非線性規(guī)律變化的主要原因。

    2.3 影響因素分析

    2.3.1 盲孔深度的影響

    結合本工作預制盲孔方法可知,電火花打孔法的火花放電時間決定了盲孔深度,但很難保證盲孔深度恒定不變。相關研究表明,盲孔深度不僅會影響Q235鋼試樣的受力截面積,而且會影響盲孔周圍應力的分布和數(shù)值,從而影響Lcr波聲彈性系數(shù)與盲孔直徑間關系。

    2.3.2 Lcr波形畸變的影響

    互相關系數(shù)函數(shù)是基于信號波形相似程度分析而實現(xiàn)超聲波信號間時間差計算的一種方法,波形相似程度越高,互相關系數(shù)函數(shù)計算結果精度越高。前期相關研究表明[20],超聲波信號波形畸變是影響互相關系數(shù)函數(shù)計算結果精度的重要因素,波形畸變程度越高,互相關系數(shù)函數(shù)計算結果精度越低,即信號間時間差結果精度越低。結合本工作研究方法認為,Q235鋼組織結構、Lcr波應力評價系統(tǒng)等引起的“噪聲”是導致Q235鋼Lcr波發(fā)生畸變的主要原因。為此,對Q235鋼組織結構進行分析,如圖8所示。由圖8可知該組織呈明顯方向性(與軋制方向同向),滲碳體(黑色物質)沿界面分布,規(guī)律性明顯。由超聲波與介質相互作用原理可知,超聲波在異質界面(聲阻抗不同)處發(fā)生反射、透射、折射等現(xiàn)象,并形成相應的模式轉換波,進而與初始超聲波相互作用,導致超聲波波形發(fā)生畸變,且其波形畸變程度與異質界面面積、組織各向異性程度等成正比。這也是Q235鋼Lcr波形發(fā)生畸變的主要原因。

    圖8 Q235鋼顯微組織Fig.8 Microstructure of Q235 steel

    3 結論

    (1)盲孔直徑不同,Lcr波時間差隨應力變化的規(guī)律基本相同,即Lcr波時間差隨應力增大基本呈線性規(guī)律增加,但隨應力增大,Lcr波時間差的波動程度亦逐漸變大。

    (2)盲孔引起的應力集中是導致Lcr波時間差與應力呈線性變化階段最大應力值小于試樣屈服強度的主要原因,Lcr波傳播距離內試樣的平均變形狀態(tài)是影響Lcr波時間差隨應力變化規(guī)律的關鍵因素,預制盲孔深度及Lcr波形畸變也是影響應力評價結果的重要因素。

    (3)隨盲孔直徑增大,Lcr波聲彈性系數(shù)逐漸減小,并趨于平穩(wěn),其擬合結果表明采用本工作方法可實現(xiàn)任意直徑盲孔對Lcr波評價Q235鋼應力影響的研究。

    [1] 孫鶴立,張騰. 車體鋁合金氬弧焊與攪拌摩擦焊殘余應力研究[J]. 鐵道機車車輛, 2016, 36(4): 99-102.

    SUN H L, ZHANG T. Research on residual stress of argon arc welding and friction stir welding for aluminium alloy of carbody[J]. Railway Locomotive & Car, 2016, 36(4): 99-102.

    [2] 劉曉佳, 林健, 雷永平, 等. 鈦合金電子束焊接表面殘余應力的測試和有限元分析[J]. 航空材料學報, 2016, 36(4): 35-40.

    LIU X J, LIN J, LEI Y P, et al. Residual stress test and finite element analysis of titanium alloy surface obtained by electron beam welding[J]. Journal of Aeronautical Materials, 2016, 36(4): 35-40.

    [3] 趙翠華. 殘余應力超聲波測量方法研究[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2008.

    ZHAO C H. Study on ultrasonic measurement of residual stress[D]. Harbin: Harbin Institute of Technology, 2008.

    [4] 王曉, 史亦韋, 梁菁, 等. 聲彈性法測量鋁合金預拉伸板中的應力[J]. 材料工程, 2015, 43(12): 95-100.

    WANG X, SHI Y W, LIANG J, et al. Stress in pre-stretched aluminum alloy plate by acoustic elasticity[J]. Journal of Materials Engineering, 2015, 43(12): 95-100.

    [5] 齊紅宇, 劉金龍, 楊曉光, 等. 基于云紋干涉法與鉆孔法的等離子熱障涂層殘余應力實驗[J]. 航空動力學報, 2011, 26(3): 617-621.

    QI H Y, LIU J L, YANG X G, et al. Investigation on residual stress of plasma sprayed thermal barrier coatings based on moire interferometry and hole drilling method[J]. Journal of Aerospace Power, 2011, 26(3): 617-621.

    [6] 吳劍劍. 淬火硬化層殘余應力分布磁記憶檢測[D]. 南昌: 南昌航空大學, 2015.

    WU J J. Magnetic memory detection of the distribution of residual stress on hardened layer[D]. Nanchang: Nangchang Hangkong University, 2015.

    [7] YASHAR J, MEHDI A, MEHDI A N. Using finite element and ultrasonic method to evaluate welding longitudinal residual stress through the thickness in austenitic stainless steel plates[J]. Materials & Design, 2013, 45: 628-642.

    [8] YASHAR J, SEYEDALI S, MEHDI A N. Taguchi optimization and ultrasonic measurement of residual stresses in the friction stir welding[J]. Materials & Design, 2014, 55: 27-34.

    [9] SEYEDALI S, MEHDI A N, YASHAR J, et al. Using ultrasonic waves and finite element method to evaluate through-thickness residual stresses distribution in the friction stir welding of aluminum plates[J]. Materials & Design, 2013, 52: 870-880.

    [10] 路浩. 基于全包絡權重算法的超聲波法殘余應力無損測量系統(tǒng)[J]. 焊接學報, 2015, 36(3): 101-104.

    LU H. Ultrasonic residual stress measurement system based on entire envelope weighting algorithm[J]. Transactions of the China Welding Institution, 2015, 36(3): 101-104.

    [11] 路浩, 劉雪松, 孟立春, 等. 高速列車車體服役狀態(tài)殘余應力超聲波法無損測量及驗證[J]. 焊接學報, 2009, 30(4): 81-84.

    LU H, LIU X S, MENG L C, et al. Residual stress evaluation of high-speed train body structure by ultrasonic method and verification[J]. Transactions of the China Welding Institution, 2009, 30(4): 81-84.

    [12] 石一飛, 沈中華, 倪曉武, 等. 激光激發(fā)瑞利波測量鋁合金焊接殘余應力[J]. 中國激光, 2008, 35(10): 1627-1631.

    SHI Y F, SHEN Z H, NI X W, et al. Welding stress measurement using laser-generated Rayleigh wave in aluminum alloys[J]. Chinese Journal of Lasers, 2008, 35(10): 1627-1631.

    [13] 馬子奇, 劉雪松, 張世平, 等. 超聲波法曲面工件殘余應力測量[J]. 焊接學報, 2011, 32(11): 25-28.

    MA Z Q, LIU X S, ZHANG S P, et al. Curved convex face residual stress measurement by ultrasonic method[J]. Transactions of the China Welding Institution, 2011, 32(11): 25-28.

    [14] HU E Y, HE Y M, CHEN Y M. Experimental study on the surface stress measurement with Rayleigh wave detection technique[J]. Applied Acoustics, 2009, 70(2): 356-360.

    [15] 魏智, 徐蔚, 曲云霞, 等. 表面波聲彈性應力檢測基礎[J]. 河北工業(yè)大學學報, 2002, 31(1): 10-14.

    WEI Z, XU W, QU Y X, et al. On acoustoelasticity of surface acoustic wave for ultrasonic NDT[J]. Journal of Hebei University of Technology, 2002, 31(1): 10-14.

    [16] AHARON Z, KENNETH J. Measurement of acoustoelastic coefficients of Rayleigh waves in steel alloys[J]. Journal of Nondestructive Evaluation, 1982, 3(2):115-124.

    [17] HUGUES D S, KELLY J L. Second-order elastic deformation of solid[J]. Physical Review, 1951, 92(5): 1145-1149.

    [18] 劉彬, 董世運, 徐濱士, 等. 互相關函數(shù)步長影響超聲波評價涂層應力的實驗研究[J]. 材料工程, 2011, (4): 54-57.

    LIU B, DONG S Y, XU B S, et al. Study on influence of calculating length on evaluation result of stress by surface acoustic wave for coating[J]. Journal of Materials Engineering, 2011, (4): 54-57.

    [19] 張寧鋒. 基于ANSYS 的有限寬板孔邊應力集中分析[J]. 蘭州工業(yè)高等??茖W校學報, 2007, 14(1): 34-38.

    ZHANG N F. The analysis of stress concentration based on ANSYS[J]. Journal of Lanzhou Polytechnic College, 2007, 14(1): 34-38.

    [20] 劉彬. 再制造金屬涂層質量的超聲波/磁記憶綜合無損評價研究[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2012.

    LIU B. Ultrasonic and metal magnetic memory testing method for quality nondestructive evaluation of remanufacturing coating[D]. Harbin: Harbin Institute of Technology, 2012.

    (本文責編:寇鳳梅)

    Interaction Influence of Flaw and Stress on LcrWave Acoustoelastic Coefficient of Carbon Steel

    LIU Bin1,SHI Chang-liang2,MIAO Wen-bing1,DONG Shi-yun3

    (1 School of Materials Science and Engineering,Jiangsu University of Science and Technology,Zhenjiang 212003,Jiangsu,China;2 Guangdong Industrial Analysis and Testing Center,Guangzhou 510650,China; 3 National Key Laboratory for Remanufacturing,Academy of Armored Forces Engineering,Beijing 100072,China)

    Based on Lcrwave acoustoelastic theory, the influence mechanism of blind-hole on stress evaluation with Lcrwave was discussed. Combined with equivalent method, the blind-holes with different diameters were preset, the difference in time of flight between Lcrwave was calculated by cross correlation coefficient function, the Lcrwave acoustoelastic coefficient was determined with linear fitting function, the influence mechanism of blind-hole diameter on Lcrwave acoustoelastic coefficient was clarified based on the elastic plastic deformation theory and stress concentration theory of hole. The results show that when the diameter of blind-hole is different, the change rule of the difference in time of flight between Lcrwave is basically the same as stress increases, the difference in time of flight between Lcrwave increases linearly as stress increases, while the nonlinear feature is also gradually obvious, the maximum stress in linear stage is less than the yield strength of specimen, the Lcrwave acoustoelastic coefficient decreases gradually as blind-hole diameter increases, and it tends to be constant. The analysis shows that the stress concentration is seen as the main reason for above results, anisotropic structure of specimen and blind-hole depth are also the important reasons for that.

    Lcrwave acoustoelastic theory;stress concentration;nondestructive evaluation;elastic plastic deformation;interaction influence

    10.11868/j.issn.1001-4381.2015.001437

    TG115.2

    A

    1001-4381(2017)07-0097-06

    國家自然科學基金資助項目(51305172);中國博士后科學基金資助項目(2016M591795);廣東省公益研究與能力建設專項(2015-A030401073)

    2015-11-24;

    2017-02-07

    劉彬(1983-),男,博士,講師,研究方向:表面工程、質量無損評價及壽命預測,聯(lián)系地址:江蘇省鎮(zhèn)江市夢溪路2號江蘇科技大學材料科學與工程學院(212003),E-mail:liubindely@163.com

    猜你喜歡
    盲孔時間差傳播速度
    代謝綜合征患者臂踝脈搏波傳播速度與頸動脈粥樣硬化的關系
    盲孔內反倒角的加工工藝研究
    量子定位系統(tǒng)中符合計數(shù)與到達時間差的獲取
    射孔槍耐壓性能的影響因素分析
    測井技術(2019年4期)2019-12-25 03:28:20
    新雷
    一類廣義canmassa—Holm方程的無限傳播速度與漸近行為
    基于BP網(wǎng)絡的GIS局部放電聲電聯(lián)合檢測故障定位方法
    立體聲音樂節(jié)目后期制作中聲像定位的探討
    演藝科技(2017年8期)2017-09-25 16:08:33
    一種跨層盲孔制作及對位方式研究
    一種快速填盲孔的工藝及原理研究
    久久精品成人免费网站| 人成视频在线观看免费观看| 一区二区三区国产精品乱码| 电影成人av| 日韩成人在线观看一区二区三区| 桃色一区二区三区在线观看| 欧美激情 高清一区二区三区| 久久久国产一区二区| 午夜福利免费观看在线| 国产精品乱码一区二三区的特点 | 欧美久久黑人一区二区| 久久狼人影院| xxx96com| 天天添夜夜摸| 亚洲国产欧美日韩在线播放| 成年人黄色毛片网站| 午夜福利,免费看| 在线十欧美十亚洲十日本专区| 在线观看日韩欧美| 日韩欧美国产一区二区入口| 中出人妻视频一区二区| 久久天堂一区二区三区四区| 日本一区二区免费在线视频| 另类亚洲欧美激情| 中文字幕人妻熟女乱码| 久久久国产欧美日韩av| 日本三级黄在线观看| 1024香蕉在线观看| 老司机深夜福利视频在线观看| 久9热在线精品视频| 日韩成人在线观看一区二区三区| 视频区图区小说| 亚洲av日韩精品久久久久久密| 黑人巨大精品欧美一区二区mp4| 色综合欧美亚洲国产小说| 视频区欧美日本亚洲| 亚洲一区二区三区色噜噜 | 国产欧美日韩精品亚洲av| 成人精品一区二区免费| 久久香蕉激情| 亚洲精品久久成人aⅴ小说| 女人爽到高潮嗷嗷叫在线视频| 成人亚洲精品一区在线观看| 激情视频va一区二区三区| 极品人妻少妇av视频| 黄色视频不卡| 亚洲精品av麻豆狂野| 欧美激情极品国产一区二区三区| 国产视频一区二区在线看| 国产精品久久电影中文字幕| 少妇裸体淫交视频免费看高清 | 久久伊人香网站| 亚洲人成电影免费在线| 国产一区二区三区在线臀色熟女 | 免费在线观看日本一区| 欧美日韩一级在线毛片| 欧美日韩中文字幕国产精品一区二区三区 | 免费人成视频x8x8入口观看| avwww免费| 久久精品91蜜桃| 咕卡用的链子| 国产极品粉嫩免费观看在线| 男女之事视频高清在线观看| 美女午夜性视频免费| 国产无遮挡羞羞视频在线观看| 亚洲色图综合在线观看| 亚洲人成电影免费在线| 黄色视频不卡| 日本五十路高清| 757午夜福利合集在线观看| 成人国语在线视频| 夜夜爽天天搞| 1024香蕉在线观看| 九色亚洲精品在线播放| www.www免费av| 久久中文看片网| 免费在线观看影片大全网站| 成人影院久久| 电影成人av| www.自偷自拍.com| 在线观看免费日韩欧美大片| 午夜久久久在线观看| 国产有黄有色有爽视频| 成人国产一区最新在线观看| 久久国产精品人妻蜜桃| 日日干狠狠操夜夜爽| 18禁国产床啪视频网站| 午夜福利影视在线免费观看| 一区二区三区精品91| 男女做爰动态图高潮gif福利片 | 777久久人妻少妇嫩草av网站| 91老司机精品| 亚洲一区高清亚洲精品| 久久精品国产亚洲av高清一级| 亚洲国产毛片av蜜桃av| 在线观看66精品国产| 日韩免费高清中文字幕av| 亚洲一区高清亚洲精品| 亚洲av第一区精品v没综合| 91麻豆精品激情在线观看国产 | 波多野结衣av一区二区av| 两性午夜刺激爽爽歪歪视频在线观看 | 宅男免费午夜| 久久久久久久久中文| 电影成人av| 亚洲午夜理论影院| 国产免费av片在线观看野外av| 黄片播放在线免费| 十八禁网站免费在线| 久久热在线av| 91成人精品电影| 成人18禁在线播放| 乱人伦中国视频| 成人三级做爰电影| 国产欧美日韩一区二区三| 亚洲国产欧美网| 成年人黄色毛片网站| 国产精品久久久av美女十八| 亚洲色图综合在线观看| 亚洲五月天丁香| 国产精品免费一区二区三区在线| 女同久久另类99精品国产91| 国产成人av激情在线播放| 最近最新中文字幕大全电影3 | 欧美一级毛片孕妇| 亚洲熟女毛片儿| 可以免费在线观看a视频的电影网站| 国产午夜精品久久久久久| 亚洲欧美一区二区三区黑人| 热99re8久久精品国产| 伦理电影免费视频| 老熟妇乱子伦视频在线观看| 亚洲成人精品中文字幕电影 | 国产精品久久久久久人妻精品电影| 精品熟女少妇八av免费久了| 国产精品一区二区在线不卡| 黄色成人免费大全| 99久久人妻综合| 免费高清视频大片| 国产av又大| 一级片'在线观看视频| 亚洲 欧美一区二区三区| 老司机深夜福利视频在线观看| 亚洲国产欧美日韩在线播放| 搡老熟女国产l中国老女人| 波多野结衣高清无吗| 电影成人av| 成年人免费黄色播放视频| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜精品久久久久久| 国产精品一区二区精品视频观看| 天堂俺去俺来也www色官网| 国产av在哪里看| 嫩草影视91久久| 波多野结衣高清无吗| 欧美午夜高清在线| 久久久精品欧美日韩精品| 欧美成狂野欧美在线观看| 性欧美人与动物交配| 免费在线观看日本一区| 亚洲av五月六月丁香网| 久久久国产欧美日韩av| 叶爱在线成人免费视频播放| 欧美乱妇无乱码| 老司机福利观看| 咕卡用的链子| 欧美+亚洲+日韩+国产| 18禁国产床啪视频网站| 国产精品久久久久久人妻精品电影| 很黄的视频免费| 搡老乐熟女国产| 亚洲欧美一区二区三区久久| 久久人人97超碰香蕉20202| 午夜福利一区二区在线看| 亚洲九九香蕉| aaaaa片日本免费| 国产精品一区二区免费欧美| 老熟妇乱子伦视频在线观看| 美女午夜性视频免费| 国产精品偷伦视频观看了| 久久人人爽av亚洲精品天堂| 亚洲精品中文字幕在线视频| 天堂√8在线中文| 又大又爽又粗| 高清毛片免费观看视频网站 | 午夜成年电影在线免费观看| 一边摸一边抽搐一进一出视频| 国产男靠女视频免费网站| 久久久水蜜桃国产精品网| 久久久久国内视频| 国产精品一区二区在线不卡| 身体一侧抽搐| a级片在线免费高清观看视频| 免费高清在线观看日韩| 亚洲第一欧美日韩一区二区三区| netflix在线观看网站| 午夜视频精品福利| 国产亚洲欧美在线一区二区| 无人区码免费观看不卡| 亚洲五月色婷婷综合| 日本黄色日本黄色录像| 欧美另类亚洲清纯唯美| 亚洲全国av大片| 欧美成人性av电影在线观看| 国产日韩一区二区三区精品不卡| 中文亚洲av片在线观看爽| 少妇被粗大的猛进出69影院| www.www免费av| 国产片内射在线| 久久精品亚洲精品国产色婷小说| 亚洲精品一区av在线观看| 久久99一区二区三区| 国产午夜精品久久久久久| 亚洲全国av大片| www.自偷自拍.com| 岛国在线观看网站| 久久婷婷成人综合色麻豆| 午夜福利欧美成人| 在线观看一区二区三区激情| 天堂√8在线中文| 深夜精品福利| 他把我摸到了高潮在线观看| 亚洲国产精品sss在线观看 | 97碰自拍视频| 女人被躁到高潮嗷嗷叫费观| 精品国产乱子伦一区二区三区| 亚洲精品久久午夜乱码| 一本大道久久a久久精品| 久久影院123| 国产精品一区二区精品视频观看| 国产精品一区二区三区四区久久 | 在线播放国产精品三级| 久久草成人影院| 欧美精品一区二区免费开放| 一区二区日韩欧美中文字幕| 欧美不卡视频在线免费观看 | 999久久久国产精品视频| 欧洲精品卡2卡3卡4卡5卡区| 两个人看的免费小视频| 性色av乱码一区二区三区2| 一夜夜www| 亚洲精品国产区一区二| 亚洲欧美精品综合一区二区三区| 欧美日韩av久久| 欧美精品啪啪一区二区三区| 欧美激情久久久久久爽电影 | 亚洲九九香蕉| 一夜夜www| 少妇裸体淫交视频免费看高清 | 国产成+人综合+亚洲专区| 两性午夜刺激爽爽歪歪视频在线观看 | a级毛片在线看网站| 亚洲狠狠婷婷综合久久图片| 视频在线观看一区二区三区| 热99re8久久精品国产| 欧美成人免费av一区二区三区| 久久人人爽av亚洲精品天堂| 99久久99久久久精品蜜桃| 美女国产高潮福利片在线看| 在线观看66精品国产| 新久久久久国产一级毛片| 99久久人妻综合| 天天躁夜夜躁狠狠躁躁| 欧美在线黄色| 啦啦啦免费观看视频1| 国产一卡二卡三卡精品| 女生性感内裤真人,穿戴方法视频| 天天添夜夜摸| 午夜福利在线观看吧| 国产熟女午夜一区二区三区| 久久草成人影院| 精品一区二区三区av网在线观看| 国产伦一二天堂av在线观看| 搡老熟女国产l中国老女人| 少妇 在线观看| 国产精品亚洲一级av第二区| 人人妻,人人澡人人爽秒播| 精品一区二区三区视频在线观看免费 | 亚洲精品国产精品久久久不卡| 天堂√8在线中文| 国产亚洲精品久久久久5区| 午夜福利在线观看吧| 国产欧美日韩精品亚洲av| 男女之事视频高清在线观看| 久久久久久久久久久久大奶| 三级毛片av免费| 国产精品久久电影中文字幕| 在线av久久热| 国产人伦9x9x在线观看| 一级毛片高清免费大全| 亚洲欧美精品综合一区二区三区| 天堂俺去俺来也www色官网| 免费看a级黄色片| 亚洲成人免费av在线播放| 亚洲国产欧美日韩在线播放| 日韩精品免费视频一区二区三区| av视频免费观看在线观看| 国产精品免费一区二区三区在线| 大陆偷拍与自拍| a在线观看视频网站| 99久久精品国产亚洲精品| 丰满迷人的少妇在线观看| 久久人妻av系列| 香蕉久久夜色| 亚洲色图av天堂| 亚洲全国av大片| 他把我摸到了高潮在线观看| 婷婷六月久久综合丁香| 久久狼人影院| 亚洲精品美女久久av网站| 午夜亚洲福利在线播放| 后天国语完整版免费观看| 欧美日韩福利视频一区二区| 成熟少妇高潮喷水视频| 美女福利国产在线| 亚洲一区中文字幕在线| 国产视频一区二区在线看| 亚洲成人免费av在线播放| 国产精品乱码一区二三区的特点 | 亚洲成a人片在线一区二区| 岛国在线观看网站| 亚洲自拍偷在线| 久久精品成人免费网站| 欧美大码av| 国产成人欧美| 国产精品av久久久久免费| 美女扒开内裤让男人捅视频| 一边摸一边做爽爽视频免费| 十八禁网站免费在线| 国产区一区二久久| 老司机亚洲免费影院| 18禁国产床啪视频网站| 日韩视频一区二区在线观看| 99在线视频只有这里精品首页| 黑人欧美特级aaaaaa片| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 中亚洲国语对白在线视频| 久久人人爽av亚洲精品天堂| av中文乱码字幕在线| 高清黄色对白视频在线免费看| 成人影院久久| 身体一侧抽搐| 久久精品国产清高在天天线| 欧美乱妇无乱码| 亚洲av片天天在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲黑人精品在线| xxx96com| 欧美日本亚洲视频在线播放| 大型av网站在线播放| 人人澡人人妻人| 亚洲精品在线观看二区| 精品免费久久久久久久清纯| 国产高清videossex| 欧美日韩精品网址| 亚洲激情在线av| 又大又爽又粗| 日日摸夜夜添夜夜添小说| 91成人精品电影| 色精品久久人妻99蜜桃| 国产激情欧美一区二区| 身体一侧抽搐| 日本免费a在线| 丰满迷人的少妇在线观看| 亚洲专区中文字幕在线| a级毛片在线看网站| 久久人人爽av亚洲精品天堂| 亚洲美女黄片视频| 日韩av在线大香蕉| 欧美最黄视频在线播放免费 | 亚洲九九香蕉| 亚洲欧洲精品一区二区精品久久久| 香蕉丝袜av| 国产精品一区二区免费欧美| 精品国产一区二区三区四区第35| 久久精品亚洲av国产电影网| 不卡av一区二区三区| 中文字幕人妻丝袜一区二区| 99久久综合精品五月天人人| 老熟妇仑乱视频hdxx| 免费看十八禁软件| 欧美乱码精品一区二区三区| 亚洲熟妇中文字幕五十中出 | 一区福利在线观看| 日本wwww免费看| 超碰97精品在线观看| 琪琪午夜伦伦电影理论片6080| 久久久国产成人精品二区 | 男女之事视频高清在线观看| 黄色丝袜av网址大全| 亚洲全国av大片| 日韩欧美免费精品| 国产精品一区二区在线不卡| 韩国精品一区二区三区| 满18在线观看网站| 国产aⅴ精品一区二区三区波| 日韩精品青青久久久久久| 国产精品 国内视频| 亚洲精品av麻豆狂野| 亚洲国产精品999在线| 久久久久九九精品影院| 精品高清国产在线一区| 一级a爱片免费观看的视频| 男人的好看免费观看在线视频 | 亚洲五月色婷婷综合| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美一区视频在线观看| 国产国语露脸激情在线看| 亚洲欧美精品综合一区二区三区| 国产成人精品在线电影| 级片在线观看| 亚洲国产毛片av蜜桃av| 高清毛片免费观看视频网站 | 久热这里只有精品99| 麻豆国产av国片精品| 欧美日韩av久久| 长腿黑丝高跟| 在线播放国产精品三级| 久久人人精品亚洲av| 亚洲人成电影观看| 日韩国内少妇激情av| 亚洲中文av在线| www.精华液| 老熟妇仑乱视频hdxx| 最近最新免费中文字幕在线| 精品第一国产精品| 国产在线观看jvid| 日韩有码中文字幕| 9191精品国产免费久久| 久久精品aⅴ一区二区三区四区| 中文字幕人妻熟女乱码| 一二三四在线观看免费中文在| 琪琪午夜伦伦电影理论片6080| 大码成人一级视频| 亚洲欧美激情在线| 嫩草影院精品99| 欧美久久黑人一区二区| 在线视频色国产色| 一区二区三区精品91| 成年版毛片免费区| 亚洲狠狠婷婷综合久久图片| 我的亚洲天堂| 亚洲精品久久成人aⅴ小说| 免费在线观看黄色视频的| 色老头精品视频在线观看| 琪琪午夜伦伦电影理论片6080| 在线播放国产精品三级| 最好的美女福利视频网| 成人手机av| 国产精品98久久久久久宅男小说| av在线天堂中文字幕 | 欧美日韩亚洲高清精品| 国产伦人伦偷精品视频| 黄片播放在线免费| 国产精品久久久人人做人人爽| e午夜精品久久久久久久| 亚洲中文字幕日韩| 国产成人精品无人区| 免费在线观看日本一区| 老司机亚洲免费影院| 国产精品综合久久久久久久免费 | 亚洲精品一二三| 久久久国产一区二区| 大型黄色视频在线免费观看| www.www免费av| 国产熟女xx| 老汉色av国产亚洲站长工具| 中文字幕色久视频| 成人av一区二区三区在线看| 久久九九热精品免费| 色老头精品视频在线观看| 91成年电影在线观看| 国产精品一区二区在线不卡| 在线观看免费视频网站a站| 女人被狂操c到高潮| 乱人伦中国视频| 久久香蕉激情| 久热爱精品视频在线9| 大型av网站在线播放| 欧美人与性动交α欧美软件| 欧美激情极品国产一区二区三区| 人人澡人人妻人| 欧美午夜高清在线| av网站在线播放免费| 亚洲国产精品999在线| 亚洲欧美激情综合另类| 亚洲全国av大片| 欧美日韩福利视频一区二区| 不卡一级毛片| 视频在线观看一区二区三区| 老熟妇乱子伦视频在线观看| 亚洲欧美一区二区三区黑人| 一本综合久久免费| 国产精品自产拍在线观看55亚洲| 99精品在免费线老司机午夜| 亚洲久久久国产精品| 国产在线精品亚洲第一网站| a级片在线免费高清观看视频| 露出奶头的视频| 久久人妻av系列| 少妇被粗大的猛进出69影院| 美女扒开内裤让男人捅视频| 亚洲第一av免费看| 女人被狂操c到高潮| 婷婷精品国产亚洲av在线| 日韩有码中文字幕| 亚洲精品中文字幕在线视频| 这个男人来自地球电影免费观看| 不卡av一区二区三区| 国产有黄有色有爽视频| 搡老熟女国产l中国老女人| 欧美日韩中文字幕国产精品一区二区三区 | 在线永久观看黄色视频| 国产精品 欧美亚洲| 国产野战对白在线观看| 国产成人一区二区三区免费视频网站| 欧美日韩亚洲高清精品| 国产av又大| 丰满饥渴人妻一区二区三| 99riav亚洲国产免费| 九色亚洲精品在线播放| 国产精品久久久久久人妻精品电影| 欧美人与性动交α欧美软件| 少妇裸体淫交视频免费看高清 | 看黄色毛片网站| www.999成人在线观看| 91老司机精品| 少妇被粗大的猛进出69影院| 多毛熟女@视频| 另类亚洲欧美激情| 黄色丝袜av网址大全| 国产欧美日韩一区二区三| av天堂久久9| 亚洲,欧美精品.| 久久国产精品人妻蜜桃| 亚洲av美国av| 亚洲欧洲精品一区二区精品久久久| 在线永久观看黄色视频| 搡老乐熟女国产| 黄片小视频在线播放| av天堂久久9| 免费在线观看完整版高清| 久久国产精品人妻蜜桃| 91麻豆精品激情在线观看国产 | 久久狼人影院| 美女扒开内裤让男人捅视频| 黄色女人牲交| 高清av免费在线| 日韩 欧美 亚洲 中文字幕| 国产激情久久老熟女| av视频免费观看在线观看| 成人三级做爰电影| 国产成人系列免费观看| 校园春色视频在线观看| 亚洲精品久久成人aⅴ小说| 亚洲国产中文字幕在线视频| 日韩欧美一区视频在线观看| 亚洲第一欧美日韩一区二区三区| 在线观看一区二区三区| 午夜久久久在线观看| 在线永久观看黄色视频| 久久国产精品人妻蜜桃| 亚洲,欧美精品.| 级片在线观看| 国产主播在线观看一区二区| 中文字幕最新亚洲高清| 岛国视频午夜一区免费看| 精品一区二区三区av网在线观看| 国产又色又爽无遮挡免费看| 黑人猛操日本美女一级片| 久久精品国产亚洲av香蕉五月| 国产精品永久免费网站| 18禁国产床啪视频网站| 久久久久国产一级毛片高清牌| 欧美在线黄色| 国内久久婷婷六月综合欲色啪| 国产高清激情床上av| 久久国产精品人妻蜜桃| 亚洲久久久国产精品| 成人精品一区二区免费| 天堂中文最新版在线下载| 成熟少妇高潮喷水视频| 成人特级黄色片久久久久久久| 桃红色精品国产亚洲av| 欧美一区二区精品小视频在线| 日本撒尿小便嘘嘘汇集6| 久久欧美精品欧美久久欧美| 欧美日韩一级在线毛片| 999久久久精品免费观看国产| 国产精品爽爽va在线观看网站 | 99久久久亚洲精品蜜臀av| 国产一区二区在线av高清观看| 午夜免费激情av| 国产一区二区三区视频了| 亚洲一卡2卡3卡4卡5卡精品中文| 久久婷婷成人综合色麻豆| 美女 人体艺术 gogo| 欧美av亚洲av综合av国产av| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利,免费看| 日日爽夜夜爽网站| 精品国产乱码久久久久久男人| 欧美日韩av久久| tocl精华| 波多野结衣av一区二区av| 大码成人一级视频| 国产蜜桃级精品一区二区三区| 夜夜躁狠狠躁天天躁| 午夜福利免费观看在线| 中文字幕精品免费在线观看视频| 久久久久久久久久久久大奶| 国产无遮挡羞羞视频在线观看| 视频区图区小说| 欧美在线黄色|