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

    古箏琴板的振動(dòng)特性初探研究?

    2021-07-01 08:51:46嵇建書裴永超盧曉杰鄭梅生
    林產(chǎn)工業(yè) 2021年6期
    關(guān)鍵詞:古箏振型固有頻率

    嵇建書 裴永超 羅 敏 盧曉杰 鄭梅生

    (南京林業(yè)大學(xué)機(jī)械電子工程學(xué)院,江蘇 南京 210037)

    木材是木質(zhì)樂器的重要組成部分,樂器質(zhì)量等級(jí)主要由木材的材質(zhì)決定,因此材料的選擇對(duì)于樂器的制作至關(guān)重要。由于缺乏理論指導(dǎo)和科學(xué)手段,從業(yè)者在古箏等木質(zhì)樂器的實(shí)際制作中,主要通過看、聽、敲、掂等主觀判斷方法進(jìn)行篩選,因此不僅樂器質(zhì)量得不到保障,同時(shí)也造成木材的不合理利用與浪費(fèi)。據(jù)行業(yè)相關(guān)統(tǒng)計(jì),用于生產(chǎn)樂器共鳴板的原木出材率為10%~20%,出材率較低[1]。有研究表明,通過木材振動(dòng)試驗(yàn),分析各類木材的振動(dòng)特性,將木材振動(dòng)與聲學(xué)聯(lián)系起來,有助于發(fā)現(xiàn)木材振動(dòng)特性與木質(zhì)樂器質(zhì)量之間的關(guān)系,從而幫助生產(chǎn)廠家篩選、分類材料,提高我國的樂器生產(chǎn)制造水平和木材的利用率[2-5]。為此,本文以古箏琴板為例,通過理論與試驗(yàn)分析相結(jié)合對(duì)木制樂器板材的振動(dòng)特性進(jìn)行研究,以期為古箏琴板及其他樂器木質(zhì)板材的篩選與生產(chǎn)提供了一種有效的科學(xué)手段。

    1 古箏琴板振動(dòng)特性理論分析

    木制樂器構(gòu)造各異,但發(fā)聲原理基本相同。以古箏為例,其結(jié)構(gòu)如圖1所示。木制樂器的初始振動(dòng)由弦引起,琴弦振動(dòng)導(dǎo)致面板振動(dòng),弦聲通過面板、背板等形成的共鳴箱放大并傳播出去,從而使人們享受到美妙的音樂[6]。

    圖1 古箏結(jié)構(gòu)Fig. 1 Guzheng structure

    我國民族樂器共振板的材料主要為泡桐與云杉屬木材等[7-9]。共振面板是彈撥類樂器中的一個(gè)重要構(gòu)件。琴板在振動(dòng)時(shí),既有彎曲振動(dòng),又有扭轉(zhuǎn)振動(dòng)。本文首先對(duì)琴板振動(dòng)特性進(jìn)行理論分析,以求出振動(dòng)的固有頻率理論解,為后續(xù)琴板振動(dòng)試驗(yàn)提供理論依據(jù)。

    1.1 彈性體的彎曲振動(dòng)

    勻質(zhì)彈性體如矩形截面梁在xy平面內(nèi)作橫向彎曲振動(dòng),如圖2所示。設(shè)其長度為l,橫截面面積為A,ρ為每單位體積的質(zhì)量,EI為彎曲剛度且為常量[10-11]。

    圖2 板的受力分析Fig. 2 Force analysis of board

    從中取一個(gè)長度為dx的微段,上面有剪力Q(x,t)、彎矩M(x,t) ,在任何時(shí)間t,其橫向位移可表示為y=(x,t)。忽略轉(zhuǎn)動(dòng)慣量和剪切變形的影響,由圖3可知微段沿y方向力的平衡關(guān)系為:

    圖3 板的微段受力分析Fig. 3 Force analysis of micro segment of board

    解得其前5個(gè)值如表1所示。

    表1 頻率方程解Tab.1 Solution of frequency equation

    最后,得到彈性體的彎曲振動(dòng)固有頻率計(jì)算公式如下[7]:

    琴板屬于彈性體,因?yàn)槭歉飨虍愋圆牧?,如以板的振?dòng)研究會(huì)使問題復(fù)雜化。因此嘗試用公式(8)計(jì)算琴板的彎曲振動(dòng)固有頻率,并進(jìn)行試驗(yàn)驗(yàn)證。將琴板的基本參數(shù)代入式(8),得到其彎曲振動(dòng)固有頻率f(Hz)近似計(jì)算公式如下[12]:

    本試驗(yàn)用古箏琴板材質(zhì)為泡桐,長l=1.7 m,寬b=0.355 m,厚h=0.01 m,質(zhì)量m=1.62 kg,泊松比μ=0.4,彈性模量E=5.8×109Pa[9]。將具體數(shù)值代入公式(9),解得該琴板前五階彎曲固有頻率如表2所示。

    表2 琴板彎曲振動(dòng)固有頻率Tab.2 Natural frequency of bending vibration of the board

    1.2 彈性體的扭轉(zhuǎn)振動(dòng)

    勻質(zhì)彈性體如等直桿繞x軸作扭轉(zhuǎn)振動(dòng),如圖4所示。設(shè)桿的截面抗扭剛度、抗扭常數(shù)和剪切彈性模量分別為GJt(x)、Jt(x)和G[13]。

    圖4 桿的受力分析Fig. 4 Force analysis of the rod

    圖5 桿的微段受力分析Fig. 5 Force analysis of micro-segment of rod

    假設(shè)桿扭轉(zhuǎn)振動(dòng)時(shí),截面的翹曲可忽略。在平面內(nèi),截面將保持繞x軸作微擺動(dòng)[14],取微段dx,x截面的角位移表示為θ(x,t),微段對(duì)x軸的轉(zhuǎn)動(dòng)慣量Ip(x),圓截面的極慣性矩為Jp。根據(jù)力矩平衡關(guān)系得[15]:

    然后由邊界條件解得β值及振型函數(shù)H(x) ,最后得到板的扭轉(zhuǎn)振動(dòng)固有頻率計(jì)算公式如下:

    與琴板彎曲振動(dòng)同理,嘗試用公式(14)計(jì)算琴板的扭轉(zhuǎn)振動(dòng)固有頻率,并進(jìn)行試驗(yàn)驗(yàn)證。將琴板的基本參數(shù)代入式(14),得到其扭轉(zhuǎn)振動(dòng)固有頻率f(Hz)近似計(jì)算公式如下:

    本文試驗(yàn)用古箏琴板的剪切彈性模量G=2.07×109Pa,其他數(shù)值同上。將具體數(shù)值代入公式(15),解得該琴板的扭轉(zhuǎn)固有頻率如表3所示。

    表3 琴板扭轉(zhuǎn)振動(dòng)固有頻率Tab.3 Natural frequency of torsional vibration of the board

    1.3 有阻尼系統(tǒng)的衰減振動(dòng)

    自由振動(dòng)物體的定義為,振動(dòng)的振幅是恒定的,振動(dòng)將無限期地繼續(xù)下去。然而,實(shí)際的運(yùn)動(dòng)并非如此。自由振動(dòng)的物體不但受到恢復(fù)力的作用,還受到阻力的作用,振幅會(huì)隨著時(shí)間逐漸衰減,最終趨于零并停止振動(dòng),如圖6所示[17],其振動(dòng)函數(shù)如下:

    圖6 阻尼振動(dòng)運(yùn)動(dòng)圖Fig. 6 Damping vibration motion diagram

    通常把對(duì)振動(dòng)過程的阻力稱為阻尼,阻尼大小取決于物體的形狀、尺寸及介質(zhì)性質(zhì)。琴板振動(dòng)過程受阻尼的影響,呈現(xiàn)衰減振動(dòng)。該分析結(jié)果將為試驗(yàn)結(jié)果分析中的時(shí)域圖篩選提供一定的理論依據(jù)。

    1.4 轉(zhuǎn)動(dòng)慣量、剪切變形對(duì)板振動(dòng)的影響

    當(dāng)板的橫截面尺寸與長度之比不是很小,或者在分析高階振型時(shí),整個(gè)板被節(jié)點(diǎn)平面分成若干比較短的小段,這時(shí)需要考慮轉(zhuǎn)動(dòng)慣量和剪切變形的影響[18]。該板的受力分析圖如下:

    圖7 鐵木辛柯梁?jiǎn)卧w受力圖Fig. 7 Force diagram of Timoshenko beam element

    取微單元dx,繪制剪切力和彎矩引起的變形。在剪力為零的條件下,單元體dx的中心線垂直于橫截面。ψ為截面彎矩引起的夾角,β為剪力引起的剪切角,梁的軸線由彎矩和剪力引起的實(shí)際夾角為θ。通過建立動(dòng)力方程、運(yùn)動(dòng)方程并求解,可得固有頻率pi(Hz)計(jì)算公式如下:

    式(18)中p0(Hz)為不計(jì)轉(zhuǎn)動(dòng)慣量和剪切變形影響下的固有頻率:

    k為截面常數(shù),矩形截面k=1.2,圓形截面k=1.11;EI為彎曲剛度,N·m2[10-11];G為剪切彈性模量,Pa;ρ為每單位體積的質(zhì)量,kg/m3;A為橫截面面積,m2。

    由式(18)可知,考慮了剪切變形和轉(zhuǎn)動(dòng)慣量后,系統(tǒng)的固有頻率減小。由于系統(tǒng)的質(zhì)量和剛度可以決定固有頻率,剪切變形和轉(zhuǎn)動(dòng)慣量的影響不僅增加了系統(tǒng)的有效質(zhì)量,而且降低了系統(tǒng)的有效剛度,導(dǎo)致固有頻率下降。這種影響對(duì)高階頻率更為明顯。

    綜上所述,琴板的實(shí)際固有頻率要比理論計(jì)算的數(shù)值小。該分析結(jié)果將為試驗(yàn)結(jié)果分析中的時(shí)域圖篩選提供一定的理論依據(jù)。

    2 琴板振動(dòng)特性試驗(yàn)臺(tái)設(shè)計(jì)

    2.1 設(shè)計(jì)目的與技術(shù)要求

    設(shè)計(jì)并制作一種樂器木質(zhì)板材振動(dòng)特性測(cè)試裝置,用于測(cè)試古箏、吉他等樂器的木質(zhì)板材振動(dòng)頻率、彈性系數(shù)、阻尼系數(shù)等振動(dòng)特性參數(shù),為木質(zhì)樂器生產(chǎn)廠家提供科學(xué)篩選板材的方法。

    技術(shù)要求:1)板材懸吊后只限制其垂直向下自由移動(dòng),懸掛位置為手動(dòng)操作,有限位裝置; 2)支撐件位置可調(diào),可以測(cè)試各類不同規(guī)格樂器板材及其他物體的振動(dòng)特性;3)傳感器等測(cè)試元件位置、角度可調(diào),裝夾方便;4)試驗(yàn)臺(tái)易拆裝。

    2.2 試驗(yàn)臺(tái)整體結(jié)構(gòu)設(shè)計(jì)

    琴板測(cè)試臺(tái)由底板、立柱、掛棍、擋銷、橫梁、傳感器、傳感器夾、傳感器支撐架等構(gòu)成(圖8)。底板用于支撐其他構(gòu)件,其上開有橫、豎兩個(gè)方向的T型槽;兩根立柱及橫梁用于懸掛被測(cè)板材或其他物體,橫梁與立柱之間用G型夾固定;立柱可沿底板的T型槽作縱向移動(dòng)與固定,方便測(cè)量不同大小板材;兩個(gè)支撐架外側(cè)均設(shè)有掛棍與擋銷,擋銷設(shè)置在掛棍上方,同時(shí)掛棍上設(shè)有凹槽懸掛試驗(yàn)板材。傳感器支撐架可沿底板的T型槽作橫向移動(dòng)與固定,方便傳感器靠近被測(cè)物體以接收信號(hào);磁性表座安于傳感器支撐架上,并可快速固定;傳感器裝于磁性表座上,位置、角度均可調(diào),裝夾方便;整個(gè)試驗(yàn)臺(tái)易拆裝[19]。

    圖8 試驗(yàn)臺(tái)整體結(jié)構(gòu)三維圖Fig. 8 3D drawing of the overall structure of the test bench

    3 琴板振動(dòng)特性試驗(yàn)與結(jié)果分析

    3.1 琴板自由振動(dòng)頻率測(cè)試試驗(yàn)

    3.1.1 設(shè)備

    自由振動(dòng)試驗(yàn)設(shè)備均購于北京東方振動(dòng)噪聲研究所,包括振動(dòng)分析軟件DASP-V11、云智慧采集分析儀INV3062、聲壓傳感器INV9206、力錘、電荷放大器、電腦、試驗(yàn)臺(tái)等[20-22]。試驗(yàn)系統(tǒng)如圖9所示。

    圖9 自由振動(dòng)試驗(yàn)系統(tǒng)Fig.9 Free vibration test system

    3.1.2 試驗(yàn)方法

    利用試驗(yàn)臺(tái)將琴板懸吊,吊線選在節(jié)線位置;將聲壓傳感器的信號(hào)接收端固定于距離琴板中間3~5 cm處,接口端通過數(shù)據(jù)線與采集儀的2號(hào)端口連接;力錘通過電荷放大器與采集儀的1號(hào)端口連接;設(shè)置DASP軟件的初始環(huán)境;進(jìn)行錘擊測(cè)試,錘擊位置不固定,力錘對(duì)每個(gè)被測(cè)位置依次進(jìn)行3次敲擊,振動(dòng)信號(hào)通過聲壓傳感器傳至采集儀,由采集儀將數(shù)據(jù)傳輸至振動(dòng)分析系統(tǒng)進(jìn)行處理。所有試驗(yàn)室溫均為25 ℃,濕度為50%。

    3.1.3 結(jié)果與分析[23-25]

    經(jīng)過振動(dòng)分析系統(tǒng)處理后的時(shí)域圖如圖10所示,對(duì)應(yīng)的頻譜圖如圖11所示,數(shù)據(jù)列表如表4所示。

    圖10 時(shí)域圖Fig. 10 Time domain diagram

    圖11 頻譜圖Fig. 11 Spectrum diagram

    表4 數(shù)據(jù)列表(幅值譜Peak)Tab.4 Data list (Amplitude spectrum Peak)

    由圖10、11可見,時(shí)域信號(hào)呈現(xiàn)典型的有阻尼衰減振動(dòng)規(guī)律,頻譜圖的極值較為明顯,試驗(yàn)結(jié)果較好。幅值譜的頻率數(shù)據(jù),部分為琴板的彎曲振動(dòng)固有頻率,部分為琴板扭轉(zhuǎn)振動(dòng)固有頻率,其他為環(huán)境噪聲產(chǎn)生的干擾頻率。

    3.2 琴板振動(dòng)模態(tài)試驗(yàn)

    模態(tài)試驗(yàn)是對(duì)古箏琴板振型進(jìn)行試驗(yàn)分析。力錘錘擊琴板時(shí),3個(gè)加速度傳感器將振動(dòng)信號(hào)傳遞至采集儀,由其將信號(hào)傳遞至振動(dòng)分析系統(tǒng),采用模態(tài)與動(dòng)力學(xué)分析模塊進(jìn)行模態(tài)分析。

    3.2.1 彎曲振型

    1)主要試驗(yàn)設(shè)備同3.1.1所列。

    2)試驗(yàn)方法

    首先用2塊海綿支撐試驗(yàn)用琴板,與力錘連接的電荷放大器和3個(gè)加速度傳感器依次分別接入采集儀的1至4號(hào)通道,3個(gè)加速度傳感器分別置于選定的被測(cè)點(diǎn),并用雙面膠將其固定于琴板上;設(shè)置DASP軟件內(nèi)的信號(hào)采集模塊,對(duì)錘擊試驗(yàn)進(jìn)行預(yù)試驗(yàn),用力錘在每個(gè)被測(cè)點(diǎn)依次進(jìn)行3次敲擊,通過軟件內(nèi)的顯示情況,判斷預(yù)試驗(yàn)是否良好,若正常再進(jìn)行正式試驗(yàn)。

    3)結(jié)果與分析

    試驗(yàn)得到的數(shù)據(jù)為時(shí)域信號(hào)。首先從時(shí)域圖中選取一組信號(hào)良好的試驗(yàn)數(shù)據(jù),利用振動(dòng)分析系統(tǒng)中的模態(tài)與動(dòng)力學(xué)分析模塊對(duì)數(shù)據(jù)進(jìn)行分析計(jì)算,操作步驟依次為數(shù)據(jù)導(dǎo)入、參數(shù)設(shè)置、頻響函數(shù)計(jì)算、琴板結(jié)構(gòu)生成、輸入約束、脈沖響應(yīng)函數(shù)、特征系統(tǒng)實(shí)現(xiàn)算法等,最后得到古箏琴板彎曲振動(dòng)的前五階固有頻率、阻尼比(表5)及部分振型(圖12、13)。

    表5 彎曲振型固有頻率和阻尼比Tab.5 Bending mode natural frequency and damping

    圖12 琴板彎曲二階模態(tài)振型圖Fig. 12 Second order modal diagram of the board bending

    圖13 琴板彎曲三階模態(tài)振型圖Fig. 13 Third order modal diagram of the board bending

    二階模態(tài)振型主要表現(xiàn)為琴板凹陷與凸起的往復(fù)運(yùn)動(dòng)。三階模態(tài)振型主要表現(xiàn)為從左側(cè)起凸起、凹陷、凸起以及凹陷、凸起、凹陷的交替運(yùn)動(dòng)。通過對(duì)比觀察彎曲振型的固有頻率和自由振動(dòng)的頻率大小,可以在表4中找到與表5相近的頻率,再結(jié)合理論計(jì)算固有頻率值,整理獲得柱形圖14。

    圖14 彎曲振型固有頻率與自由振動(dòng)頻率以及理論計(jì)算固有頻率對(duì)比Fig.14 Comparison of natural frequency of bending mode and free vibration frequency and theoretical calculation natural frequency

    對(duì)比圖14中彎曲振型固有頻率和自由振動(dòng)頻率以及理論計(jì)算所得出的前五階頻率,發(fā)現(xiàn)各階頻率值均相近,平均誤差在5%以內(nèi),表明本文的研究結(jié)果可信度較高。

    3.2.2 扭轉(zhuǎn)振型

    1)主要試驗(yàn)設(shè)備同3.1.1所列。

    2)試驗(yàn)方法

    首先進(jìn)行DASP軟件的設(shè)置,根據(jù)采樣頻率約為試驗(yàn)結(jié)果最高頻率的5倍的試驗(yàn)要求,選用采樣頻率為2 048 Hz;試驗(yàn)采用4個(gè)通道,第1個(gè)為力錘通道,其他3個(gè)為加速度通道,設(shè)置每個(gè)通道的標(biāo)定值;試驗(yàn)為錘擊試驗(yàn),觸發(fā)次數(shù)為3次,采樣點(diǎn)數(shù)為采樣頻率數(shù)值的一半,即1 024個(gè),觸發(fā)力設(shè)置為10 N;進(jìn)行示波與零點(diǎn)平衡校準(zhǔn)。然后將琴板兩端放在備好的2塊海綿上,避開一階扭轉(zhuǎn)的節(jié)線位置,將加速度傳感器依次放在選定的被測(cè)點(diǎn),并用雙面膠將其固定;在琴板上依次選取7個(gè)激勵(lì)點(diǎn),力錘對(duì)每個(gè)激勵(lì)點(diǎn)進(jìn)行連續(xù)3次敲擊,并采集試驗(yàn)信號(hào)。

    3)結(jié)果與分析

    加速度傳感器將信號(hào)傳遞至分析儀進(jìn)行處理,在振動(dòng)分析系統(tǒng)中顯示為時(shí)域信號(hào),然后對(duì)時(shí)域信號(hào)進(jìn)行篩選,選用模態(tài)與動(dòng)力學(xué)分析模塊依次進(jìn)行分析,最后通過該模塊的脈沖響應(yīng)函數(shù)分析和特征系統(tǒng)實(shí)現(xiàn)算法得到琴板扭轉(zhuǎn)振動(dòng)的前五階固有頻率、阻尼比(表6)及部分振型。

    表6 扭轉(zhuǎn)振型固有頻率和阻尼比Tab.6 Torsional mode natural frequency and damping ratio

    圖15 琴板扭轉(zhuǎn)二階模態(tài)振型圖Fig. 15 Second-order mode shape of torsional vibration of the board

    圖16 琴板扭轉(zhuǎn)三階模態(tài)振型圖Fig. 16 Third-order mode shape of torsional vibration of the board

    通過對(duì)比觀察扭轉(zhuǎn)振型的固有頻率和自由振動(dòng)的頻率大小,可以在表4中找到與表6相近的頻率,再結(jié)合理論計(jì)算固有頻率值,整理得到柱形圖17。

    對(duì)比圖17中扭轉(zhuǎn)振型固有頻率和自由振動(dòng)頻率以及理論計(jì)算所得出的前五階頻率,發(fā)現(xiàn)各階頻率的值均相近,平均誤差在6%以內(nèi),表明本研究具有較高的可信度。

    圖17 扭轉(zhuǎn)振型固有頻率和自由振動(dòng)頻率以及理論計(jì)算固有頻率對(duì)比Fig.17 Comparison of natural frequency of torsional mode and free vibration frequency and theoretical calculation natural frequency

    4 結(jié)論

    本文以古箏為例研究了樂器木質(zhì)板材振動(dòng)特性,對(duì)古箏琴板彎曲振動(dòng)和扭轉(zhuǎn)振動(dòng)進(jìn)行了理論分析和有阻尼的衰減振動(dòng)分析,得到琴板的彎曲和扭轉(zhuǎn)振動(dòng)的前五階固有頻率理論計(jì)算值;設(shè)計(jì)并加工制作試驗(yàn)臺(tái);分別利用聲壓傳感器和加速度傳感器進(jìn)行了琴板的自由振動(dòng)頻率測(cè)試與彎曲和扭轉(zhuǎn)振動(dòng)模態(tài)試驗(yàn),獲得琴板自由振動(dòng)頻率試驗(yàn)值及其彎曲和扭轉(zhuǎn)振動(dòng)前五階固有頻率試驗(yàn)值及振型;琴板振動(dòng)的理論計(jì)算結(jié)果和2種試驗(yàn)結(jié)果對(duì)比表明:3個(gè)數(shù)值大小相近,平均誤差分別在5%及6%以內(nèi),驗(yàn)證了本文試驗(yàn)方法的可行性和試驗(yàn)數(shù)據(jù)的可信度。

    該研究成果對(duì)今后廠家篩選樂器木質(zhì)板材具有一定的指導(dǎo)意義。今后還需深入研究樂器聲學(xué)品質(zhì)與木材振動(dòng)性能的關(guān)系,并結(jié)合使用者對(duì)樂器成品的主觀評(píng)價(jià),綜合分析,才能客觀、科學(xué)地選取理想的樂器木質(zhì)板材。

    猜你喜歡
    古箏振型固有頻率
    我終于會(huì)彈古箏了
    關(guān)于模態(tài)綜合法的注記
    縱向激勵(lì)下大跨鋼桁拱橋高階振型效應(yīng)分析
    現(xiàn)場(chǎng)測(cè)定大型水輪發(fā)電機(jī)組軸系的固有頻率
    塔腿加過渡段輸電塔動(dòng)力特性分析
    第一次彈古箏
    小布老虎(2017年1期)2017-07-18 10:57:28
    古箏
    寶藏(2017年2期)2017-03-20 13:16:45
    我的古箏之旅
    結(jié)構(gòu)振型幾何辨識(shí)及應(yīng)用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    總溫總壓測(cè)頭模態(tài)振型變化規(guī)律研究
    av天堂在线播放| 国产区一区二久久| 久久人人爽av亚洲精品天堂| 国产亚洲av高清不卡| 免费看a级黄色片| 国产私拍福利视频在线观看| 午夜免费激情av| 校园春色视频在线观看| 欧美 亚洲 国产 日韩一| ponron亚洲| 亚洲男人天堂网一区| 日本 欧美在线| 精品久久蜜臀av无| 色婷婷久久久亚洲欧美| 午夜福利欧美成人| 国产精品久久久久久人妻精品电影| 性少妇av在线| 午夜日韩欧美国产| 欧美日本亚洲视频在线播放| 琪琪午夜伦伦电影理论片6080| 亚洲五月婷婷丁香| 天堂影院成人在线观看| 成人国产综合亚洲| 美女免费视频网站| 国产精品av久久久久免费| 免费在线观看黄色视频的| 最新美女视频免费是黄的| 午夜激情av网站| 在线视频色国产色| 69av精品久久久久久| 亚洲成国产人片在线观看| 欧美日本视频| 97人妻精品一区二区三区麻豆 | 亚洲视频免费观看视频| 久久久久国内视频| 精品一区二区三区av网在线观看| 亚洲精品国产区一区二| 亚洲精品国产精品久久久不卡| 国产色视频综合| 两个人看的免费小视频| 在线观看www视频免费| 91麻豆精品激情在线观看国产| 夜夜爽天天搞| 中文字幕最新亚洲高清| 亚洲五月色婷婷综合| 亚洲精品久久成人aⅴ小说| 欧美成人性av电影在线观看| 成人亚洲精品av一区二区| 一本大道久久a久久精品| 女人高潮潮喷娇喘18禁视频| 黄色视频,在线免费观看| 操美女的视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲第一电影网av| 国产一级毛片七仙女欲春2 | 老司机午夜十八禁免费视频| 亚洲国产欧美网| 亚洲 国产 在线| 视频区欧美日本亚洲| 别揉我奶头~嗯~啊~动态视频| 免费观看精品视频网站| 婷婷精品国产亚洲av在线| 国产一区二区三区在线臀色熟女| 久久精品国产99精品国产亚洲性色 | 亚洲少妇的诱惑av| 日本a在线网址| 自拍欧美九色日韩亚洲蝌蚪91| 欧美 亚洲 国产 日韩一| 啦啦啦免费观看视频1| 两个人视频免费观看高清| 国产成人欧美在线观看| 一边摸一边抽搐一进一出视频| 99re在线观看精品视频| 久久精品人人爽人人爽视色| 在线免费观看的www视频| 精品国产超薄肉色丝袜足j| 99久久久亚洲精品蜜臀av| av欧美777| 免费搜索国产男女视频| 9色porny在线观看| 波多野结衣高清无吗| 成在线人永久免费视频| 校园春色视频在线观看| 亚洲精品国产区一区二| 亚洲九九香蕉| 日韩欧美国产在线观看| 国产亚洲精品久久久久久毛片| 国产91精品成人一区二区三区| 国产成人影院久久av| 成人18禁高潮啪啪吃奶动态图| 999久久久精品免费观看国产| 天堂√8在线中文| 国产一区在线观看成人免费| 一区二区三区精品91| 黄片播放在线免费| 男人舔女人的私密视频| 亚洲欧美一区二区三区黑人| xxx96com| 99国产极品粉嫩在线观看| 国产精品av久久久久免费| 欧美一级毛片孕妇| 少妇 在线观看| 欧美老熟妇乱子伦牲交| 精品熟女少妇八av免费久了| 不卡av一区二区三区| 国产真人三级小视频在线观看| 电影成人av| 日韩欧美国产在线观看| 欧美一区二区精品小视频在线| 国产精品日韩av在线免费观看 | 日本免费一区二区三区高清不卡 | 少妇 在线观看| 99热只有精品国产| 黄色片一级片一级黄色片| 精品国产国语对白av| 美女大奶头视频| 亚洲国产精品成人综合色| 成人18禁在线播放| 亚洲美女黄片视频| 日韩成人在线观看一区二区三区| 亚洲国产欧美一区二区综合| 两个人视频免费观看高清| 国产欧美日韩综合在线一区二区| 亚洲国产欧美网| 97碰自拍视频| 国产高清有码在线观看视频 | 国产熟女午夜一区二区三区| av电影中文网址| 天天躁狠狠躁夜夜躁狠狠躁| 又大又爽又粗| 中文字幕av电影在线播放| 国产高清视频在线播放一区| 一夜夜www| 老司机福利观看| 欧美精品亚洲一区二区| 变态另类丝袜制服| aaaaa片日本免费| 亚洲一区二区三区不卡视频| 丰满的人妻完整版| 99精品久久久久人妻精品| 亚洲中文日韩欧美视频| 一a级毛片在线观看| 高清黄色对白视频在线免费看| 成年女人毛片免费观看观看9| 曰老女人黄片| 亚洲精品美女久久av网站| 午夜老司机福利片| 美女国产高潮福利片在线看| 国产欧美日韩精品亚洲av| bbb黄色大片| av电影中文网址| 国产主播在线观看一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲中文日韩欧美视频| 国产精品久久久人人做人人爽| 黄片大片在线免费观看| 亚洲av成人一区二区三| 久久精品国产亚洲av香蕉五月| 久久久久久久久久久久大奶| 成人欧美大片| 午夜福利欧美成人| 亚洲一卡2卡3卡4卡5卡精品中文| 天天一区二区日本电影三级 | 成人欧美大片| 一区福利在线观看| 黄色片一级片一级黄色片| 国产精品亚洲美女久久久| av视频免费观看在线观看| 欧美成狂野欧美在线观看| 99国产精品99久久久久| 欧美黑人精品巨大| 国产麻豆成人av免费视频| 国产成+人综合+亚洲专区| av欧美777| 午夜福利视频1000在线观看 | 日本a在线网址| 黄色丝袜av网址大全| 久久婷婷成人综合色麻豆| 99国产综合亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩综合在线一区二区| 九色亚洲精品在线播放| 久久精品人人爽人人爽视色| 成人18禁在线播放| 久久久久久人人人人人| 一个人观看的视频www高清免费观看 | 亚洲第一电影网av| 亚洲伊人色综图| 日韩av在线大香蕉| 国产成人免费无遮挡视频| 久久国产乱子伦精品免费另类| 一级片免费观看大全| 脱女人内裤的视频| 亚洲片人在线观看| 欧美丝袜亚洲另类 | 一级片免费观看大全| 亚洲 欧美一区二区三区| 美女大奶头视频| 国产精品 国内视频| 免费在线观看日本一区| 国产精品久久久久久亚洲av鲁大| 久久精品91无色码中文字幕| 91九色精品人成在线观看| 两个人免费观看高清视频| 国产欧美日韩一区二区精品| av免费在线观看网站| 欧美日本中文国产一区发布| 亚洲国产精品久久男人天堂| 日本五十路高清| 最近最新免费中文字幕在线| 国产熟女xx| 色综合亚洲欧美另类图片| 国产国语露脸激情在线看| 岛国在线观看网站| 午夜福利免费观看在线| 欧美一级毛片孕妇| 伊人久久大香线蕉亚洲五| 大型av网站在线播放| 成人18禁高潮啪啪吃奶动态图| 亚洲伊人色综图| 十八禁人妻一区二区| 亚洲午夜理论影院| 欧美黑人欧美精品刺激| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美国产一区二区入口| 给我免费播放毛片高清在线观看| 亚洲一区中文字幕在线| 露出奶头的视频| 久久伊人香网站| 人人妻,人人澡人人爽秒播| 午夜福利在线观看吧| 波多野结衣巨乳人妻| 久久婷婷人人爽人人干人人爱 | 天天躁狠狠躁夜夜躁狠狠躁| 三级毛片av免费| 成人国产综合亚洲| 亚洲第一av免费看| or卡值多少钱| 国产精华一区二区三区| 精品久久蜜臀av无| 91九色精品人成在线观看| 身体一侧抽搐| 欧美老熟妇乱子伦牲交| 大型av网站在线播放| 人妻久久中文字幕网| 熟妇人妻久久中文字幕3abv| 久久人人精品亚洲av| 成在线人永久免费视频| 国产极品粉嫩免费观看在线| 精品国产乱子伦一区二区三区| 动漫黄色视频在线观看| 他把我摸到了高潮在线观看| 亚洲精华国产精华精| 国产片内射在线| 欧美成人午夜精品| 久久精品国产清高在天天线| 亚洲一区高清亚洲精品| 午夜精品国产一区二区电影| 亚洲五月婷婷丁香| 国产精品99久久99久久久不卡| 国产精品自产拍在线观看55亚洲| 免费在线观看日本一区| 88av欧美| 丰满人妻熟妇乱又伦精品不卡| 久久人人精品亚洲av| 国产熟女午夜一区二区三区| 亚洲精品在线美女| 1024视频免费在线观看| 亚洲国产精品合色在线| 国产精品av久久久久免费| 国产精品一区二区免费欧美| 日韩欧美一区视频在线观看| 可以免费在线观看a视频的电影网站| 亚洲 欧美 日韩 在线 免费| 电影成人av| av电影中文网址| 在线天堂中文资源库| √禁漫天堂资源中文www| 男人舔女人的私密视频| 男人操女人黄网站| 亚洲中文字幕日韩| 一级,二级,三级黄色视频| 9191精品国产免费久久| 久久香蕉国产精品| 99精品在免费线老司机午夜| 免费在线观看完整版高清| 精品欧美一区二区三区在线| 女人高潮潮喷娇喘18禁视频| 一进一出抽搐动态| 国产成人欧美| 给我免费播放毛片高清在线观看| 亚洲三区欧美一区| 色精品久久人妻99蜜桃| 国语自产精品视频在线第100页| 亚洲狠狠婷婷综合久久图片| 午夜日韩欧美国产| 国产精品日韩av在线免费观看 | 亚洲av日韩精品久久久久久密| 制服人妻中文乱码| 黄色女人牲交| 美女扒开内裤让男人捅视频| 国产在线观看jvid| 成年女人毛片免费观看观看9| 好男人电影高清在线观看| 午夜免费激情av| 亚洲色图av天堂| 久久中文字幕人妻熟女| 色综合婷婷激情| 日本精品一区二区三区蜜桃| 91大片在线观看| 欧美黑人精品巨大| 日韩欧美在线二视频| 亚洲自拍偷在线| 操出白浆在线播放| 日本免费一区二区三区高清不卡 | 亚洲精品一区av在线观看| 99精品久久久久人妻精品| 日本黄色视频三级网站网址| 日韩大码丰满熟妇| 亚洲成人免费电影在线观看| 亚洲精品久久国产高清桃花| 国产精品一区二区精品视频观看| 日韩av在线大香蕉| 精品福利观看| 一级,二级,三级黄色视频| 天天一区二区日本电影三级 | 国产亚洲av高清不卡| 母亲3免费完整高清在线观看| 日韩国内少妇激情av| 国产成人精品在线电影| netflix在线观看网站| 亚洲熟妇熟女久久| 免费观看精品视频网站| 长腿黑丝高跟| 高潮久久久久久久久久久不卡| 18禁国产床啪视频网站| 久久久久精品国产欧美久久久| 精品一品国产午夜福利视频| 多毛熟女@视频| 国产又爽黄色视频| 大码成人一级视频| 一级a爱片免费观看的视频| 老司机午夜十八禁免费视频| 男女午夜视频在线观看| 桃色一区二区三区在线观看| 精品久久久久久,| 精品一区二区三区四区五区乱码| 自线自在国产av| 久久香蕉激情| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看 | 免费在线观看亚洲国产| 最好的美女福利视频网| 夜夜夜夜夜久久久久| 黑人巨大精品欧美一区二区蜜桃| 97超级碰碰碰精品色视频在线观看| 97碰自拍视频| 亚洲免费av在线视频| 久久香蕉精品热| 男女下面插进去视频免费观看| 欧美国产精品va在线观看不卡| 国产精品av久久久久免费| 91九色精品人成在线观看| 成在线人永久免费视频| 亚洲自偷自拍图片 自拍| 中文字幕另类日韩欧美亚洲嫩草| 一本大道久久a久久精品| 成在线人永久免费视频| 日日爽夜夜爽网站| 亚洲片人在线观看| 美女扒开内裤让男人捅视频| 国产主播在线观看一区二区| 桃色一区二区三区在线观看| 老司机靠b影院| 香蕉国产在线看| 久久精品91蜜桃| 国产成人影院久久av| 一级a爱视频在线免费观看| 色播亚洲综合网| 激情在线观看视频在线高清| 这个男人来自地球电影免费观看| 亚洲精品国产区一区二| 久久久国产成人免费| 老司机靠b影院| 在线观看免费午夜福利视频| av中文乱码字幕在线| 国内精品久久久久久久电影| 国产高清视频在线播放一区| 亚洲欧美激情综合另类| 国产av在哪里看| 欧美色欧美亚洲另类二区 | 成人18禁在线播放| 深夜精品福利| 一级毛片高清免费大全| 国产一区二区三区在线臀色熟女| 一本综合久久免费| 日日干狠狠操夜夜爽| 国产精品乱码一区二三区的特点 | 亚洲中文字幕一区二区三区有码在线看 | 日韩国内少妇激情av| 高清黄色对白视频在线免费看| 亚洲av五月六月丁香网| 又大又爽又粗| 99久久久亚洲精品蜜臀av| 午夜精品久久久久久毛片777| 女人被躁到高潮嗷嗷叫费观| 午夜老司机福利片| 国产精品久久久av美女十八| 免费久久久久久久精品成人欧美视频| 宅男免费午夜| 国产成人免费无遮挡视频| 精品久久久久久久毛片微露脸| 国产亚洲精品综合一区在线观看 | 日日摸夜夜添夜夜添小说| 欧美日本亚洲视频在线播放| 嫁个100分男人电影在线观看| 成人永久免费在线观看视频| 宅男免费午夜| 色综合欧美亚洲国产小说| 日韩欧美国产一区二区入口| 午夜福利视频1000在线观看 | 人成视频在线观看免费观看| 欧美日本视频| 99国产精品免费福利视频| 亚洲avbb在线观看| 亚洲精品在线观看二区| 好男人在线观看高清免费视频 | 国产欧美日韩一区二区三区在线| 亚洲精华国产精华精| 不卡av一区二区三区| 免费高清在线观看日韩| 淫妇啪啪啪对白视频| 亚洲激情在线av| 久久国产精品人妻蜜桃| 啦啦啦 在线观看视频| 国产1区2区3区精品| or卡值多少钱| 激情在线观看视频在线高清| 成人手机av| 午夜视频精品福利| 两人在一起打扑克的视频| 亚洲国产精品999在线| 91麻豆av在线| 久久天堂一区二区三区四区| 久久青草综合色| 成人亚洲精品一区在线观看| 少妇熟女aⅴ在线视频| 中国美女看黄片| 亚洲国产高清在线一区二区三 | 欧美日韩亚洲国产一区二区在线观看| 精品久久久久久成人av| 在线观看一区二区三区| 黄色a级毛片大全视频| 久久久久久大精品| 成人国产综合亚洲| 制服诱惑二区| 午夜福利18| 亚洲精品国产一区二区精华液| 欧美老熟妇乱子伦牲交| 国产成人精品久久二区二区免费| 亚洲aⅴ乱码一区二区在线播放 | 又紧又爽又黄一区二区| 麻豆国产av国片精品| 亚洲精品在线观看二区| www.精华液| www.999成人在线观看| 99久久99久久久精品蜜桃| 亚洲一码二码三码区别大吗| 欧美日韩黄片免| 国产成人av教育| 亚洲欧美日韩无卡精品| 在线播放国产精品三级| 国产亚洲欧美在线一区二区| 欧美绝顶高潮抽搐喷水| 久久久久国内视频| 日本免费一区二区三区高清不卡 | 亚洲av美国av| 久久久久久久久中文| 亚洲一码二码三码区别大吗| 成人特级黄色片久久久久久久| 99久久国产精品久久久| 午夜a级毛片| 午夜免费鲁丝| 久久人妻熟女aⅴ| 亚洲专区国产一区二区| bbb黄色大片| 中文字幕人成人乱码亚洲影| 久久热在线av| 亚洲人成77777在线视频| 成年版毛片免费区| 真人一进一出gif抽搐免费| 乱人伦中国视频| 国产精品久久电影中文字幕| 天堂动漫精品| 免费在线观看视频国产中文字幕亚洲| 一边摸一边抽搐一进一小说| 日韩免费av在线播放| 91麻豆av在线| 国产欧美日韩精品亚洲av| 久久精品亚洲精品国产色婷小说| 久久久久九九精品影院| 久久欧美精品欧美久久欧美| 精品久久久久久成人av| 女生性感内裤真人,穿戴方法视频| 无限看片的www在线观看| 18禁裸乳无遮挡免费网站照片 | 男女做爰动态图高潮gif福利片 | 免费在线观看完整版高清| 亚洲午夜理论影院| 一边摸一边抽搐一进一小说| 人人妻,人人澡人人爽秒播| 日本在线视频免费播放| 亚洲无线在线观看| 一级a爱视频在线免费观看| 国产极品粉嫩免费观看在线| 在线观看免费视频网站a站| 国产精品电影一区二区三区| 亚洲av第一区精品v没综合| 亚洲精品在线美女| 狂野欧美激情性xxxx| 亚洲激情在线av| 亚洲午夜理论影院| 日韩高清综合在线| 久久久久亚洲av毛片大全| 亚洲国产中文字幕在线视频| 亚洲熟妇中文字幕五十中出| 国产精品亚洲美女久久久| 嫁个100分男人电影在线观看| 女人高潮潮喷娇喘18禁视频| 成人18禁高潮啪啪吃奶动态图| 国产精品精品国产色婷婷| 国产精品1区2区在线观看.| 国产成人av教育| 国产av一区二区精品久久| 成人精品一区二区免费| 电影成人av| www.自偷自拍.com| 国产高清videossex| 日韩视频一区二区在线观看| 嫁个100分男人电影在线观看| 免费看a级黄色片| 黄色丝袜av网址大全| 久久久精品欧美日韩精品| 国产视频一区二区在线看| e午夜精品久久久久久久| 国产精品久久久人人做人人爽| 在线国产一区二区在线| 一级,二级,三级黄色视频| 亚洲伊人色综图| 成人18禁高潮啪啪吃奶动态图| 欧美激情高清一区二区三区| 一区二区三区激情视频| 亚洲免费av在线视频| 亚洲男人天堂网一区| 午夜福利影视在线免费观看| 如日韩欧美国产精品一区二区三区| 日本欧美视频一区| 国产99久久九九免费精品| 亚洲aⅴ乱码一区二区在线播放 | 欧美亚洲日本最大视频资源| 老司机午夜福利在线观看视频| 国产欧美日韩精品亚洲av| 一边摸一边抽搐一进一出视频| 1024视频免费在线观看| 亚洲成人久久性| 两人在一起打扑克的视频| 亚洲精品久久成人aⅴ小说| 一区在线观看完整版| 一级作爱视频免费观看| 久久欧美精品欧美久久欧美| 欧美中文综合在线视频| 午夜免费观看网址| 窝窝影院91人妻| 99国产综合亚洲精品| 亚洲最大成人中文| 91国产中文字幕| 婷婷精品国产亚洲av在线| 亚洲精品国产一区二区精华液| 国产精品二区激情视频| 韩国av一区二区三区四区| av网站免费在线观看视频| 亚洲五月色婷婷综合| 麻豆国产av国片精品| www.自偷自拍.com| 日本撒尿小便嘘嘘汇集6| 亚洲国产毛片av蜜桃av| 校园春色视频在线观看| 国产欧美日韩一区二区三| 黄片播放在线免费| 美女高潮到喷水免费观看| 亚洲国产日韩欧美精品在线观看 | av天堂在线播放| 亚洲欧美精品综合久久99| 亚洲视频免费观看视频| 一级毛片女人18水好多| 91成年电影在线观看| 又黄又爽又免费观看的视频| 午夜免费激情av| 嫩草影视91久久| 成年女人毛片免费观看观看9| 在线永久观看黄色视频| 精品国产乱子伦一区二区三区| 国产精品精品国产色婷婷| 国语自产精品视频在线第100页| tocl精华| 日本五十路高清| 日韩欧美在线二视频| 99国产精品免费福利视频| 99久久99久久久精品蜜桃| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利高清视频| 女警被强在线播放| 在线视频色国产色| 九色国产91popny在线| 午夜福利免费观看在线|