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

    組合金剛石圓鋸片模態(tài)和頻響分析

    2016-08-05 03:36:17葛健煜張進(jìn)生葉魯浩于家偉高麗君
    關(guān)鍵詞:振動(dòng)

    葛健煜, 黃 波, 王 志, 張進(jìn)生,, 葉魯浩, 于家偉, 高麗君

    (1. 山東大學(xué) 機(jī)械工程學(xué)院, 山東 濟(jì)南 250061; 2. 山東省石材工程技術(shù)研究中心, 山東 濟(jì)南 250061;3. 山東大學(xué) 高效潔凈機(jī)械制造教育部重點(diǎn)實(shí)驗(yàn)室, 山東 濟(jì)南 250061; 4. 山東大學(xué) 海易研究院, 山東 日照 276800)

    ?

    組合金剛石圓鋸片模態(tài)和頻響分析

    葛健煜1,2,3, 黃波1,2, 王志1,2,3, 張進(jìn)生1,2,3,4, 葉魯浩1,2,3, 于家偉1,2,3, 高麗君4

    (1. 山東大學(xué) 機(jī)械工程學(xué)院, 山東 濟(jì)南 250061; 2. 山東省石材工程技術(shù)研究中心, 山東 濟(jì)南 250061;3. 山東大學(xué) 高效潔凈機(jī)械制造教育部重點(diǎn)實(shí)驗(yàn)室, 山東 濟(jì)南 250061; 4. 山東大學(xué) 海易研究院, 山東 日照 276800)

    摘要:針對(duì)組合金剛石圓鋸片的振動(dòng)噪聲和振動(dòng)破壞問(wèn)題, 通過(guò)Virtual Lab中的聲學(xué)有限元模塊對(duì)組合金剛石圓鋸片進(jìn)行模態(tài)分析和單點(diǎn)激勵(lì)下的頻響分析, 同時(shí)對(duì)相同規(guī)格的單片鋸也進(jìn)行相應(yīng)分析, 分析了鋸片的固有頻率和主振型, 以及在軸向變頻載荷下鋸片的位移響應(yīng). 通過(guò)組合鋸與單片鋸的對(duì)比分析, 得出其相同點(diǎn)在于:為了避免在外界沖擊作用下產(chǎn)生的頻率與固有頻率相近而發(fā)生共振, 鋸片參數(shù)設(shè)置均要避開(kāi)共振轉(zhuǎn)速; 其不同點(diǎn)在于:在耦合作用下, 組合鋸出現(xiàn)兩個(gè)鋸片不同振型的情況, 其中一鋸片軸向振動(dòng)幅值為0, 振動(dòng)特性良好, 在實(shí)際工作中, 可以合理利用這些固有頻率振型. 研究結(jié)果為實(shí)際工作中組合鋸的振動(dòng)控制提供參考.

    關(guān)鍵詞:金剛石圓鋸片; 組合鋸; 振動(dòng); 模態(tài)分析; 頻響分析

    近年來(lái), 由于組合金剛石圓盤鋸能大幅提高生產(chǎn)效率和效益、 降低生產(chǎn)成本, 組合鋸開(kāi)始廣泛應(yīng)用于石材加工行業(yè). 當(dāng)前, 石材行業(yè)用組合金剛石圓盤鋸, 主要為不等直徑的大小片塔型組合鋸, 實(shí)現(xiàn)一次完全鋸解需要多次走刀, 效率較低, 使用等直徑鋸片組合可以極大地提高生產(chǎn)效率, 但等直徑組合鋸鋸片振動(dòng)噪聲問(wèn)題嚴(yán)重, 導(dǎo)致加工誤差大, 尤其體現(xiàn)在大板的加工中, 加工出板材厚度不一致性高, 資源浪費(fèi), 污染嚴(yán)重. 針對(duì)金剛石單片鋸的振動(dòng)噪聲問(wèn)題, 國(guó)內(nèi)外學(xué)者已經(jīng)進(jìn)行了廣泛而深入的研究,Beljo-Lucic教授研究了鋸切過(guò)程中的空氣動(dòng)力學(xué)噪聲[1],Vobolis介紹了一種能夠較為精確確定圓鋸片自由振動(dòng)的振型、 固有頻率以及振幅的方法[2],B.W.Huang研究了加工中鋸片的橫向振動(dòng)及其存在的穩(wěn)定性問(wèn)題[3],MingSongZhang研究發(fā)現(xiàn)鋸片開(kāi)槽影響噪聲衰減, 并研究了槽參數(shù)對(duì)鋸片振動(dòng)特性的影響[4],Beljo-Lucic通過(guò)實(shí)驗(yàn)驗(yàn)證了切割參數(shù)的優(yōu)化可以改善鋸片振動(dòng)性能, 降低噪聲[5],PohlM提出了在鋸片兩邊添加一種可以自主發(fā)電的壓電阻尼器達(dá)到降低噪聲的目的[6],G.S.Schajer認(rèn)為對(duì)圓鋸片采用導(dǎo)向裝置能減小圓鋸片的振動(dòng)噪聲[7]. 但針對(duì)組合金剛石圓盤鋸的振動(dòng)噪聲問(wèn)題, 尤其等直徑金剛石組合鋸的振動(dòng)噪聲問(wèn)題, 目前還沒(méi)有研究.

    本文對(duì)Φ1 650mm金剛石圓鋸片兩片等直徑組合進(jìn)行研究, 基于VirtualLab, 采用有限元方法, 計(jì)算分析組合金剛石圓鋸片固有振動(dòng)特性和外載荷作用下的振動(dòng)特性, 為低噪聲組合金剛石圓鋸片設(shè)計(jì)提供一定參考.

    1組合金剛石圓鋸片的模態(tài)分析

    1.1金剛石圓鋸片的振動(dòng)學(xué)方程

    采用彈性力學(xué)有限元法, 得到金剛石圓盤鋸基體的動(dòng)力學(xué)方程為[8]

    (1)

    (2)

    式中:ωi為第i階模態(tài)的固有頻率,i=1 ,2,…,n.

    將金剛石圓鋸片動(dòng)力學(xué)方程從物理坐標(biāo)系轉(zhuǎn)換到模態(tài)坐標(biāo)系, 得到振動(dòng)位移

    (3)

    則金剛石圓鋸片在物理坐標(biāo)系下速度為

    (4)

    第r階模態(tài)在結(jié)構(gòu)某點(diǎn)i處的振動(dòng)速度

    (5)

    式中:j為復(fù)數(shù)單位;φr,i是第r階模態(tài)在結(jié)構(gòu)i處的模態(tài)位移.

    1.2組合金剛石圓鋸片的有限元模型

    在VirtualLab中基于CATIA三維實(shí)體建模工具, 建立組合金剛石圓鋸片兩片組合的實(shí)體模型, 齒數(shù)為108, 建模結(jié)構(gòu)參數(shù)如表 1 所示, 有限元模型材料參數(shù)如表 2 所示.

    表 1 組合金剛石圓鋸片結(jié)構(gòu)參數(shù)

    表 2 組合金剛石圓鋸片的材料參數(shù)

    模型賦予2D屬性, 網(wǎng)格劃分時(shí)采用板殼Shell單元, 網(wǎng)格形狀為三角形, 劃分尺寸為20mm, 對(duì)組合圓鋸片進(jìn)行面網(wǎng)格劃分. 邊界條件:對(duì)組合圓鋸片主軸兩個(gè)端面處施加全約束. 組合圓鋸片有限元模型如圖 1 所示.

    圖 1 組合金剛石圓鋸片有限元網(wǎng)格Fig.1 Finite element mesh of combined diamond saw blades

    1.3組合金剛石圓鋸片模態(tài)結(jié)果分析

    實(shí)際工況下,Φ1 650組合鋸轉(zhuǎn)速多為n=300~360r/min, 因此鋸片工作時(shí)所受鋸切力的頻率[9]f=nz/60=540~648Hz, 其中鋸齒數(shù)z=108. 采用Lanczos法求得組合金剛石圓鋸片的固有頻率與振型, 計(jì)算范圍540~648Hz. 同時(shí), 采用相同方法對(duì)同一規(guī)格的單片鋸也進(jìn)行模態(tài)分析. 得到組合鋸與單片鋸模態(tài)分析的固有頻率對(duì)比如表 3 所示(“振型階次”從540Hz開(kāi)始計(jì)數(shù)為第一階次).

    由表 3, 可知:組合鋸固有頻率同單片鋸固有頻率表征的相同點(diǎn)在于:均隨著階次的增加, 固有頻率逐漸增加. 不同點(diǎn)在于:① 在同一頻率范圍內(nèi)(540~648Hz), 組合鋸的模態(tài)要復(fù)雜很多, 在模態(tài)耦合時(shí)出現(xiàn)了一些新的耦合頻率, 振型階次明顯多于單片鋸, 且有頻率相同的振型出現(xiàn), 如組合鋸第15階和第16階固有頻率相同, 第34階和第35階固有頻率相同, 這就導(dǎo)致當(dāng)組合鋸軸向激振力頻率達(dá)到其同時(shí)出現(xiàn)多個(gè)振型的固有頻率時(shí), 會(huì)激發(fā)鋸片更大的振動(dòng), 產(chǎn)生更大的噪聲, 而圓鋸片的主要激振力來(lái)自于它周期性的轉(zhuǎn)速, 這就指導(dǎo)組合圓鋸片鋸切時(shí)的轉(zhuǎn)速設(shè)置應(yīng)盡量遠(yuǎn)離各階固有頻率, 尤其是遠(yuǎn)離出現(xiàn)多個(gè)振型的固有頻率. ② 在頻率范圍540~648Hz內(nèi), 組合鋸固有頻率之間的差值相比于單片鋸更小, 模態(tài)更密集, 因此很容易發(fā)生共振, 故實(shí)際工況下, 考慮組合鋸鋸切參數(shù)設(shè)置時(shí), 要避開(kāi)的轉(zhuǎn)速也相應(yīng)增多.

    表 3 540~648 Hz內(nèi)組合鋸與單片鋸固有頻率

    通過(guò)總結(jié)分析, 提取了具有代表性固有模態(tài)的振型圖, 如圖 2 所示, 以兩片鋸組合為例, 振型查看取同一模態(tài)左右兩個(gè)鋸片振型圖, 同時(shí)對(duì)比相近固有頻率下單片鋸振型圖.

    圖 2 振型圖Fig.2 Vibration mode

    由圖 2 可以看出:① 組合鋸振型與單片鋸振型的相似點(diǎn)在于:在一些固有頻率點(diǎn), 組合鋸固有模態(tài)振型與單片鋸固有模態(tài)振型相近, 如圖2(c)中組合鋸在固有頻率578.007Hz時(shí)的振型圖與圖(f)單片鋸在固有頻率574.161Hz時(shí)振型圖相近, 即組合鋸保留了單片鋸的部分模態(tài)振型. ② 組合鋸振型與單片鋸振型的不同點(diǎn)在于:在耦合作用下, 組合鋸系統(tǒng)更多地呈現(xiàn)左右兩個(gè)鋸片不同振型的情況, 這些振型對(duì)組合鋸片間耦合作用敏感, 如圖2(a)中組合鋸564.877Hz振型圖及圖2(b)中組合鋸該固有頻率時(shí)左、 右兩邊鋸片振型圖, 左邊鋸片在564.877Hz時(shí)軸向振動(dòng)位移幅值為0, 而右邊鋸片軸向振動(dòng)主要集中在鋸片邊緣的周向圓上, 因此, 在實(shí)際鋸切時(shí), 可以利用固有頻率564.877Hz, 在該模態(tài)下, 只需對(duì)右邊鋸片的軸向振動(dòng)加以控制, 如可以在右邊鋸片邊緣周向圓上通過(guò)輥壓適張?zhí)幚砀纳其徠w應(yīng)力狀態(tài), 增強(qiáng)鋸片穩(wěn)定性, 也可以在鋸片邊緣周向圓處設(shè)置減小振動(dòng)的導(dǎo)向裝置, 而左邊鋸片的振動(dòng)特性良好, 可以不加處理.

    2組合金剛石圓鋸片的頻響分析

    2.1組合金剛石圓鋸片軸向載荷輸入

    軸向載荷是金剛石圓鋸片振動(dòng)產(chǎn)生的主要原因, 軸向載荷產(chǎn)生的機(jī)制有多個(gè)方面, 鋸片鋸齒的不對(duì)稱性、 被加工石材切口不規(guī)則、 鋸解過(guò)程中圓鋸片應(yīng)力狀態(tài)的變化、 圓鋸片或鋸機(jī)主軸的加工和裝配誤差等均能對(duì)鋸片產(chǎn)生軸向載荷[10].

    軸向載荷依據(jù)軸向力模型[11]

    (6)

    式中:Ft為鋸切力的切向分力; r0為鋸切力作用的半徑; x為鋸切力作用點(diǎn)處的橫向振動(dòng)位移; θ為鋸切力作用點(diǎn)的位置角.

    軸向載荷激勵(lì)力頻率范圍取500~700Hz, 激勵(lì)點(diǎn)位置如圖 3 所示.

    圖 3 軸向激勵(lì)力施加位置Fig.3 Applied position of axial excitation force

    2.2組合金剛石圓鋸片頻響函數(shù)求解參數(shù)設(shè)置

    在VirtualLab中, 通過(guò)模態(tài)疊加法計(jì)算組合金剛石圓鋸片的振動(dòng)頻響函數(shù), 以0~2 000Hz的組合金剛石圓鋸片振動(dòng)模態(tài)作為模態(tài)基建立組合金剛石圓鋸片模態(tài)空間[9], 計(jì)算組合金剛石圓鋸片在軸向激勵(lì)下的位移響應(yīng), 求解頻率范圍選540~650Hz, 計(jì)算步長(zhǎng)為1Hz. 此外計(jì)算過(guò)程中為防止共振時(shí)結(jié)構(gòu)的振動(dòng)響應(yīng)峰值趨向無(wú)限大, 需要添加模態(tài)阻尼. 金剛石圓鋸片的模態(tài)阻尼損失因子由試驗(yàn)測(cè)得, 根據(jù)已有的數(shù)據(jù), 計(jì)算取組合金剛石圓鋸片的各階模態(tài)損失因子為0.2%[12].

    2.3組合金剛石圓鋸片頻響結(jié)果分析

    取組合圓鋸片任一鋸片XOY平面上靠近激勵(lì)點(diǎn)一側(cè)的3個(gè)節(jié)點(diǎn)分析, 同時(shí)取相同載荷激勵(lì)下, 單片鋸相同位置節(jié)點(diǎn)分析, 節(jié)點(diǎn)位置如圖 4 所示, 位移響應(yīng)曲線如圖 5 所示.

    圖 4 分析節(jié)點(diǎn)選擇位置Fig.4 Analysis node selection

    圖 5 位移響應(yīng)曲線Fig.5 Displacement response curve

    由圖 5(a), (b)可知:組合鋸?fù)瑔纹徫灰祈憫?yīng)的相同之處是:① 在激振力下, 某些頻率下鋸片發(fā)生共振, 出現(xiàn)位移峰值, 且組合鋸與單片鋸主要位移峰值所在的頻率段基本對(duì)應(yīng), 如均在560~570Hz及590~600Hz時(shí)發(fā)生共振, 產(chǎn)生位移峰值, 因此, 在實(shí)際工作中, 組合鋸?fù)瑯右荛_(kāi)共振敏感區(qū)域使用鋸片, 避開(kāi)共振轉(zhuǎn)速. ② 距激振點(diǎn)越近的節(jié)點(diǎn), 引發(fā)共振的概數(shù)和程度均越大, 如節(jié)點(diǎn)C處位移峰明顯多于其他兩個(gè)節(jié)點(diǎn), 且位移峰值要高于其他兩個(gè)節(jié)點(diǎn).

    對(duì)比圖 5(c), (d)可知:組合鋸?fù)瑔纹徫灰祈憫?yīng)的不同之處是:在激勵(lì)力下, 組合鋸與單片鋸的位移最高峰值出現(xiàn)的頻段不同, 組合鋸位移主峰出現(xiàn)在590~600Hz之間, 而單片鋸位移主峰出現(xiàn)在560~570Hz之間, 且組合鋸的位移主峰值要高于單片鋸的位移主峰值, 這是因?yàn)椋?工況下組合鋸的每個(gè)鋸片均同時(shí)受軸向載荷, 鋸片間振動(dòng)耦合造成的, 故在實(shí)際工作中, 鋸片參數(shù)設(shè)置要特別參考組合鋸出現(xiàn)位移主峰的頻率段.

    3結(jié)論

    通過(guò)建立等直徑組合金剛石圓鋸片的有限元模型, 并對(duì)其進(jìn)行固有模態(tài)分析和單點(diǎn)軸向激勵(lì)作用下的頻響分析. 得出以下結(jié)論:

    1) 組合金剛石圓鋸片相比于金剛石單片鋸, 固有模態(tài)更密集, 且有頻率相同的振型出現(xiàn), 為了防止共振, 建議鋸切時(shí)鋸片的轉(zhuǎn)速設(shè)置應(yīng)盡量遠(yuǎn)離各階固有頻率, 尤其是遠(yuǎn)離出現(xiàn)多個(gè)振型的固有頻率.

    2) 在耦合作用下, 組合鋸保留了單片鋸的部分模態(tài)振型, 但更多地呈現(xiàn)兩個(gè)鋸片不同振型的情況:其一鋸片軸向振動(dòng)幅值為0, 振動(dòng)特性良好, 在實(shí)際工作中, 可以合理利用這些固有頻率振型, 對(duì)振動(dòng)明顯鋸片的軸向振動(dòng)可以通過(guò)輥壓適張?zhí)幚怼?設(shè)置導(dǎo)向裝置等加以控制.

    3) 在軸向激振力下, 組合鋸?fù)瑔纹徱粯樱?產(chǎn)生共振, 出現(xiàn)位移峰值, 但在耦合作用下, 其位移主峰所在的頻率段不同, 且峰值高于單片鋸, 故在實(shí)際工作中, 組合鋸?fù)瑯右荛_(kāi)共振敏感區(qū)域使用鋸片, 特別要考慮出現(xiàn)位移主峰的頻率段.

    參考文獻(xiàn):

    [1]Beljo-LucicR,VladoG.Suppressionofwhistlingnoiseinidlingcircularsawbydampingofclampingsystem[J].DrevarskyVyskum/WoodResearch, 2003(3):19-26.

    [2]UkvalberqieneK,VobolisJ.Experimentalstudiesofwoodcircularsawforms[J].WoodResearch, 2005, 50:47-58.

    [3]Bo-WunH,Jao-HwaK.Variationinthestabilityofarotatingbladediskwithalocalcrackdefect[J].JournalofSoundandVibration, 2006(3):486-502.

    [4]ZhangMingsong,ZhuPuxian.Diamondcircularsawbladesnoiseattenuationandfiniteelementmodalanalysis[J].AppliedMechanicsandMaterials. 2013(10):129-133.

    [5]Beljo-LucicR.Somepossibilitiesforreducingcircularsawidlingnoise[J].JournalofWoodScience, 2001, 47(5):389-393.

    [6]PohlM,RoseM.Piezoelectricshuntdampingofacircularsawbladewithautonomouspowersupplyfornoiseandvibrationreduction[J].JournalofSoundandVibration, 2015(1):20-31.

    [7]SchajerGS,WangSA.Effectofworkpieceinteractiononcircularsawcuttingstability[J].EuropeanJournalofWoodandWoodProducts, 2001(5): 388-393.

    [8]李慶華, 譚慶昌. 高速旋轉(zhuǎn)圓盤刀具薄刀片橫向振動(dòng)的有限元分析[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版), 2007, 37(4): 814-818.

    LiQinghua,TanQingchang.Finiteelementanalysisoflateralvibrationofhigh-speedrotatingdiskcuttersliceblade[J].JournalofJilinUniversity, 2007, 37(4): 814-818. (inChinese)

    [9]徐張明, 沈榮瀛. 利用FEM/IBEM計(jì)算流體介質(zhì)中的殼體的結(jié)構(gòu)聲耦合問(wèn)題[J].振動(dòng)工程學(xué)報(bào), 2002, 15(3):363-367.

    XuZhangming,ShenRongying.Structural-acousticcouplingproblemofanimmersedshellbyFEM/IBEM[J].JournalofVibrationEngineering, 2002, 15(3):363-367. (inChinese)

    [10]房懷英, 徐西鵬. 石材鋸切加工中圓鋸片的振動(dòng)研究[D]. 泉州:華僑大學(xué), 2012.

    [11]白碩瑋, 張進(jìn)生. 薄型鋸片鋸切硬脆石材橫向振動(dòng)模型[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2015(2):372-378.

    BaiShuowei,ZhangJinsheng.Analyticalmodeloftransversevibrationofthintypecircularsawinsawinghardandbrittlestone[J].TransactionsoftheChineseSocietyforAgriculturalMachinery, 2015(2):372-378. (inChinese)

    [12]VasquesCMA,CardosoLC.Viscoelasticdampingtechnologies:finiteelementmodelingandapplicationtocircularsawblades[M].SpringerNetherlands:VibrationandStructuralAcousticsAnalysis, 2011.

    文章編號(hào):1673-3193(2016)04-0375-06

    收稿日期:2015-12-17

    基金項(xiàng)目:泰山產(chǎn)業(yè)領(lǐng)軍人才工程專項(xiàng)經(jīng)費(fèi)資助項(xiàng)目; 山東省自主創(chuàng)新及成果轉(zhuǎn)化專項(xiàng)(2014CGZH0802); 山東省自然科學(xué)基金(ZR2012EEM032)

    作者簡(jiǎn)介:葛健煜(1991-), 男, 碩士生, 主要從事石材制品高效綠色加工技術(shù)與裝備的研究.

    通信作者:張進(jìn)生(1962-), 男, 教授, 博士研究生導(dǎo)師, 主要從事機(jī)械產(chǎn)品數(shù)字化設(shè)計(jì), 石材制品高效清潔生產(chǎn)技術(shù)與裝備的研究.

    中圖分類號(hào):TH113

    文獻(xiàn)標(biāo)識(shí)碼:A

    doi:10.3969/j.issn.1673-3193.2016.04.010

    Analysis on Modal and Frequency Response ofCombinedDiamondCircularSawBlades

    GE Jian-yu1,2,3, HUANG Bo1,2, WANG Zhi1,2,3, ZHANG Jin-sheng1,2,3,4,YELu-hao1,2,3,YUJia-wei1,2,3,GAOLi-jun4

    (1. School of Mechanical Engineering, Shandong University, Jinan 250061, China;2.StoneEngineeringTechnologyResearchCenterofShandongProvince,Jinan250061,China;3.KeyLaboratoryofHigh-EfficiencyandCleanMechanicalManufacture(ShandongUniversity),MinistryofEducation,Jinan250061,China; 4.HaiyiGroupInstituteofShandongUniversity,Rizhao276800,China)

    Abstract:Aiming at the vibration noise and vibration damage of combined diamond saw blades, through modal analysis and single point excitation frequency response analysis on combined diamond saw blades in acoustic finite element model of Virtual Lab,and the corresponding analysis for the same size of single chip,the natural frequency, the main vibration mode and the displacement response were analyzed. Comparing two kinds of analysis,come to the same points:In order to avoid the resonance when the natural frequency and the frequency of external shocks close,the resonance speed should be avoided when setting parameters; The difference lies in: under the coupling effect,combined saw shows two different vibration types,the axial vibration amplitude of one blade is 0, fine vibration characteristic, these vibration modes can be reasonably used in practical work.Simulation calculation practice references and recommendations for the vibration control in actual operation.

    Key words:diamond saw blades; combined saw; vibration; modal analysis; frequency response analysis

    猜你喜歡
    振動(dòng)
    振動(dòng)的思考
    某調(diào)相機(jī)振動(dòng)異常診斷分析與處理
    振動(dòng)與頻率
    This “Singing Highway”plays music
    具非線性中立項(xiàng)的廣義Emden-Fowler微分方程的振動(dòng)性
    中立型Emden-Fowler微分方程的振動(dòng)性
    基于ANSYS的高速艇艉軸架軸系振動(dòng)響應(yīng)分析
    船海工程(2015年4期)2016-01-05 15:53:26
    主回路泵致聲振動(dòng)分析
    UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
    帶有強(qiáng)迫項(xiàng)的高階差分方程解的振動(dòng)性
    国产淫语在线视频| 如何舔出高潮| 七月丁香在线播放| 欧美日韩亚洲综合一区二区三区_| 女性被躁到高潮视频| 成人漫画全彩无遮挡| 亚洲精品,欧美精品| 51午夜福利影视在线观看| 国产av精品麻豆| 国产欧美亚洲国产| 宅男免费午夜| 少妇精品久久久久久久| 日韩一区二区视频免费看| tube8黄色片| 十八禁高潮呻吟视频| 亚洲欧美色中文字幕在线| 秋霞伦理黄片| 国产又色又爽无遮挡免| 宅男免费午夜| 欧美av亚洲av综合av国产av | 少妇人妻精品综合一区二区| 亚洲人成77777在线视频| 免费观看人在逋| 国产精品.久久久| 女性生殖器流出的白浆| 亚洲久久久国产精品| 国产午夜精品一二区理论片| 午夜福利视频在线观看免费| 巨乳人妻的诱惑在线观看| 新久久久久国产一级毛片| 亚洲伊人色综图| 亚洲少妇的诱惑av| 无遮挡黄片免费观看| 99久久人妻综合| 国产野战对白在线观看| 巨乳人妻的诱惑在线观看| 国产又色又爽无遮挡免| 日韩,欧美,国产一区二区三区| 色94色欧美一区二区| 国产精品久久久久久人妻精品电影 | 18禁观看日本| 久久精品熟女亚洲av麻豆精品| 五月开心婷婷网| 日韩大码丰满熟妇| 最近最新中文字幕免费大全7| 91老司机精品| 国产伦人伦偷精品视频| 大香蕉久久网| 免费观看人在逋| 午夜av观看不卡| xxx大片免费视频| 99九九在线精品视频| 丝袜在线中文字幕| 青春草视频在线免费观看| 黑丝袜美女国产一区| 日本色播在线视频| 亚洲欧洲精品一区二区精品久久久 | 卡戴珊不雅视频在线播放| 亚洲欧美日韩另类电影网站| 欧美 日韩 精品 国产| 日韩欧美一区视频在线观看| 免费在线观看完整版高清| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区黑人| 电影成人av| 国产精品一国产av| 最近最新中文字幕免费大全7| 久久久国产一区二区| 丰满乱子伦码专区| 亚洲精品久久午夜乱码| 啦啦啦在线免费观看视频4| 日本爱情动作片www.在线观看| 久久 成人 亚洲| 一本久久精品| 水蜜桃什么品种好| 中文字幕亚洲精品专区| 一个人免费看片子| 亚洲,欧美,日韩| 亚洲国产成人一精品久久久| 久久人人爽av亚洲精品天堂| 亚洲国产最新在线播放| 免费日韩欧美在线观看| 亚洲熟女毛片儿| 国产精品蜜桃在线观看| 精品亚洲成国产av| 午夜福利一区二区在线看| 99re6热这里在线精品视频| 国产亚洲精品第一综合不卡| 人妻 亚洲 视频| 一本一本久久a久久精品综合妖精| 久久精品aⅴ一区二区三区四区| 99re6热这里在线精品视频| 午夜福利一区二区在线看| 在线天堂中文资源库| 国产高清不卡午夜福利| 满18在线观看网站| 国产精品国产三级专区第一集| 亚洲欧美成人综合另类久久久| 狂野欧美激情性xxxx| 国产成人欧美在线观看 | 超色免费av| 2018国产大陆天天弄谢| 久久久亚洲精品成人影院| 中文欧美无线码| 99香蕉大伊视频| 欧美另类一区| 热re99久久国产66热| 亚洲三区欧美一区| 在线观看免费视频网站a站| 黑丝袜美女国产一区| 久久99精品国语久久久| 国产一区二区在线观看av| 欧美av亚洲av综合av国产av | 高清视频免费观看一区二区| 肉色欧美久久久久久久蜜桃| 国产黄频视频在线观看| 在线观看免费午夜福利视频| 中文字幕色久视频| 欧美精品av麻豆av| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| tube8黄色片| 国产在线免费精品| 久久精品国产亚洲av高清一级| 一区二区三区精品91| 别揉我奶头~嗯~啊~动态视频 | av在线老鸭窝| 国产精品免费视频内射| 波多野结衣av一区二区av| 人人妻人人爽人人添夜夜欢视频| 黄色视频不卡| 亚洲国产精品一区二区三区在线| 亚洲色图 男人天堂 中文字幕| 99香蕉大伊视频| 男女边吃奶边做爰视频| 王馨瑶露胸无遮挡在线观看| 亚洲一区中文字幕在线| 九九爱精品视频在线观看| avwww免费| 亚洲成色77777| 久久精品亚洲av国产电影网| 又粗又硬又长又爽又黄的视频| 看免费av毛片| 亚洲精品日本国产第一区| 亚洲精品国产区一区二| 国产av一区二区精品久久| 国产伦人伦偷精品视频| 日韩成人av中文字幕在线观看| 丝瓜视频免费看黄片| 最近最新中文字幕大全免费视频 | 女人高潮潮喷娇喘18禁视频| 电影成人av| 日本av手机在线免费观看| 日韩制服丝袜自拍偷拍| 亚洲精品视频女| 亚洲av电影在线观看一区二区三区| 亚洲三区欧美一区| 九草在线视频观看| 18在线观看网站| 免费av中文字幕在线| 国产高清国产精品国产三级| 2021少妇久久久久久久久久久| 韩国av在线不卡| 亚洲少妇的诱惑av| 欧美老熟妇乱子伦牲交| 丁香六月欧美| 在线免费观看不下载黄p国产| 国产色婷婷99| 国产成人午夜福利电影在线观看| 黄频高清免费视频| 桃花免费在线播放| 亚洲欧美成人综合另类久久久| 高清不卡的av网站| 亚洲国产欧美日韩在线播放| av线在线观看网站| 国产伦理片在线播放av一区| 18在线观看网站| 少妇的丰满在线观看| av国产精品久久久久影院| 久久久久精品性色| 女人高潮潮喷娇喘18禁视频| 啦啦啦 在线观看视频| 一级片免费观看大全| 视频区图区小说| 欧美日韩视频高清一区二区三区二| av不卡在线播放| 免费黄网站久久成人精品| 亚洲,欧美精品.| av网站免费在线观看视频| 七月丁香在线播放| 中文字幕高清在线视频| 亚洲国产欧美日韩在线播放| 亚洲精华国产精华液的使用体验| 久久这里只有精品19| 国产精品久久久久久精品电影小说| 老汉色av国产亚洲站长工具| 久久久国产一区二区| 国产精品久久久久久久久免| 一级片免费观看大全| 国产1区2区3区精品| 一个人免费看片子| 久久精品久久久久久噜噜老黄| 国产精品嫩草影院av在线观看| 高清视频免费观看一区二区| 夫妻午夜视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看黄色视频的| 免费日韩欧美在线观看| 亚洲中文av在线| 日韩av免费高清视频| 如何舔出高潮| 色精品久久人妻99蜜桃| 日韩欧美精品免费久久| 18禁国产床啪视频网站| 免费久久久久久久精品成人欧美视频| 亚洲精品一区蜜桃| 波多野结衣一区麻豆| 国产老妇伦熟女老妇高清| 在线天堂中文资源库| 美女扒开内裤让男人捅视频| av线在线观看网站| 男人操女人黄网站| 欧美日韩精品网址| 波野结衣二区三区在线| 看免费av毛片| kizo精华| 最新在线观看一区二区三区 | 男人爽女人下面视频在线观看| 国产精品二区激情视频| 精品午夜福利在线看| 婷婷成人精品国产| 国产黄色免费在线视频| 99久国产av精品国产电影| 老汉色av国产亚洲站长工具| 制服人妻中文乱码| 极品少妇高潮喷水抽搐| 美女主播在线视频| 亚洲精品第二区| 美女脱内裤让男人舔精品视频| 黄片播放在线免费| 亚洲精品久久成人aⅴ小说| 99精品久久久久人妻精品| 女性被躁到高潮视频| 亚洲国产日韩一区二区| 欧美精品av麻豆av| 欧美黑人精品巨大| 中文字幕色久视频| 新久久久久国产一级毛片| 91aial.com中文字幕在线观看| 亚洲一区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| av国产精品久久久久影院| 久久久精品免费免费高清| 日日撸夜夜添| 18禁观看日本| 午夜av观看不卡| 美女国产高潮福利片在线看| 精品国产露脸久久av麻豆| 久久久久久久精品精品| 考比视频在线观看| 在线观看三级黄色| 精品亚洲成国产av| 色婷婷av一区二区三区视频| 中文字幕亚洲精品专区| 亚洲色图 男人天堂 中文字幕| 欧美日韩av久久| 波野结衣二区三区在线| 中文字幕av电影在线播放| av视频免费观看在线观看| 一本色道久久久久久精品综合| 国产成人欧美| 超色免费av| 中文乱码字字幕精品一区二区三区| 精品酒店卫生间| 欧美精品av麻豆av| 考比视频在线观看| 男女边吃奶边做爰视频| 丝袜喷水一区| 精品久久久久久电影网| 一区二区三区激情视频| 亚洲激情五月婷婷啪啪| 9热在线视频观看99| 国产av码专区亚洲av| 啦啦啦中文免费视频观看日本| 另类亚洲欧美激情| 亚洲综合色网址| av女优亚洲男人天堂| 波野结衣二区三区在线| 精品少妇内射三级| 美女高潮到喷水免费观看| 国产精品一区二区在线不卡| 久久毛片免费看一区二区三区| 爱豆传媒免费全集在线观看| 91成人精品电影| 亚洲欧洲国产日韩| 成人亚洲欧美一区二区av| 天美传媒精品一区二区| 亚洲国产精品999| 2018国产大陆天天弄谢| 国产淫语在线视频| 国产精品人妻久久久影院| 老司机靠b影院| 1024视频免费在线观看| 欧美日韩精品网址| 满18在线观看网站| 免费在线观看视频国产中文字幕亚洲 | 秋霞在线观看毛片| 亚洲av国产av综合av卡| 日韩 亚洲 欧美在线| 人体艺术视频欧美日本| 中文字幕高清在线视频| 成人亚洲精品一区在线观看| 欧美久久黑人一区二区| 久久天堂一区二区三区四区| 成人三级做爰电影| 国产色婷婷99| 日韩av免费高清视频| 国产精品久久久av美女十八| 国产一区亚洲一区在线观看| 亚洲国产精品国产精品| 一区二区av电影网| 国产爽快片一区二区三区| 青草久久国产| 在线观看www视频免费| 久久久亚洲精品成人影院| 成年人午夜在线观看视频| 五月开心婷婷网| 久久精品久久精品一区二区三区| 亚洲成人国产一区在线观看 | 亚洲精品第二区| 中文天堂在线官网| 亚洲国产欧美网| 美女中出高潮动态图| 大香蕉久久网| 国产日韩欧美在线精品| 欧美精品人与动牲交sv欧美| 免费av中文字幕在线| 在线免费观看不下载黄p国产| 久久精品国产亚洲av涩爱| 一区二区av电影网| 五月天丁香电影| 久久av网站| 日本午夜av视频| 人人澡人人妻人| 精品亚洲成a人片在线观看| 免费女性裸体啪啪无遮挡网站| 午夜久久久在线观看| 亚洲精品在线美女| 美女高潮到喷水免费观看| 中文字幕色久视频| 亚洲欧美日韩另类电影网站| 19禁男女啪啪无遮挡网站| 一本—道久久a久久精品蜜桃钙片| 在线 av 中文字幕| 亚洲久久久国产精品| 亚洲av男天堂| 亚洲精品成人av观看孕妇| 黑丝袜美女国产一区| 久久ye,这里只有精品| 精品一区二区三卡| 一级片'在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 成人三级做爰电影| 欧美变态另类bdsm刘玥| 亚洲精品国产色婷婷电影| 久热这里只有精品99| 午夜激情久久久久久久| 日韩电影二区| 久久久精品94久久精品| 1024香蕉在线观看| 久久影院123| 尾随美女入室| 亚洲 欧美一区二区三区| 高清视频免费观看一区二区| 中文字幕av电影在线播放| 国产成人精品福利久久| 少妇猛男粗大的猛烈进出视频| 性高湖久久久久久久久免费观看| av卡一久久| 久久精品人人爽人人爽视色| 国产成人a∨麻豆精品| 大香蕉久久网| 精品久久久精品久久久| 成人毛片60女人毛片免费| 亚洲欧美中文字幕日韩二区| 国产精品 国内视频| 免费黄色在线免费观看| 日本av免费视频播放| 久久精品久久久久久久性| 国产一区二区三区av在线| 国产亚洲一区二区精品| 国产精品一国产av| 久久久久久人人人人人| 操美女的视频在线观看| 国产精品麻豆人妻色哟哟久久| 国产成人欧美在线观看 | 亚洲在久久综合| 亚洲,一卡二卡三卡| 19禁男女啪啪无遮挡网站| 国产免费视频播放在线视频| 国产熟女欧美一区二区| 在现免费观看毛片| 亚洲一区二区三区欧美精品| 超碰97精品在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品国产a三级三级三级| 七月丁香在线播放| 亚洲国产看品久久| 亚洲精品国产色婷婷电影| 黄片小视频在线播放| 国产精品女同一区二区软件| 黄色视频不卡| 国产在视频线精品| 美女大奶头黄色视频| 一级黄片播放器| 亚洲美女黄色视频免费看| 欧美亚洲日本最大视频资源| 菩萨蛮人人尽说江南好唐韦庄| 免费黄频网站在线观看国产| 日韩av免费高清视频| 国产精品久久久人人做人人爽| 国产亚洲最大av| 亚洲精品一二三| 色视频在线一区二区三区| 九草在线视频观看| 日韩免费高清中文字幕av| 男女下面插进去视频免费观看| 最近的中文字幕免费完整| 成人手机av| 18禁动态无遮挡网站| av福利片在线| 美女视频免费永久观看网站| 免费人妻精品一区二区三区视频| 老司机靠b影院| 久久久精品国产亚洲av高清涩受| 在线观看人妻少妇| 欧美国产精品一级二级三级| 亚洲国产看品久久| 可以免费在线观看a视频的电影网站 | 久久天躁狠狠躁夜夜2o2o | 老司机靠b影院| 一区二区av电影网| 美女午夜性视频免费| 欧美中文综合在线视频| 欧美人与善性xxx| 卡戴珊不雅视频在线播放| 一级毛片黄色毛片免费观看视频| av女优亚洲男人天堂| 十八禁人妻一区二区| 婷婷成人精品国产| 亚洲国产av新网站| 国产淫语在线视频| 综合色丁香网| 99国产综合亚洲精品| 曰老女人黄片| 亚洲熟女精品中文字幕| 99久久人妻综合| 亚洲精品乱久久久久久| 久久天堂一区二区三区四区| 日日摸夜夜添夜夜爱| 99国产综合亚洲精品| 久久久精品区二区三区| 中文字幕色久视频| 悠悠久久av| 国产精品久久久人人做人人爽| 亚洲精品乱久久久久久| 伊人久久大香线蕉亚洲五| 久久久久久久久久久久大奶| 日韩免费高清中文字幕av| 黑丝袜美女国产一区| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| 大码成人一级视频| 97在线人人人人妻| 亚洲成av片中文字幕在线观看| 亚洲一码二码三码区别大吗| 午夜福利影视在线免费观看| 十八禁网站网址无遮挡| 男女国产视频网站| 美女午夜性视频免费| 中国三级夫妇交换| 国产精品久久久久久精品电影小说| 国产av码专区亚洲av| 国产高清国产精品国产三级| 日本欧美视频一区| 99热网站在线观看| 国产精品久久久久久久久免| 一区二区三区精品91| 国产 精品1| 在线看a的网站| 日日爽夜夜爽网站| 国产一区二区三区综合在线观看| 多毛熟女@视频| 搡老岳熟女国产| 日韩制服骚丝袜av| 国产成人欧美| 日本av免费视频播放| 韩国高清视频一区二区三区| 麻豆av在线久日| 2021少妇久久久久久久久久久| 国产精品av久久久久免费| 丝袜美足系列| 五月开心婷婷网| 一本大道久久a久久精品| 视频在线观看一区二区三区| 丝瓜视频免费看黄片| 亚洲国产中文字幕在线视频| 不卡视频在线观看欧美| 久久精品亚洲熟妇少妇任你| 只有这里有精品99| 国产亚洲av片在线观看秒播厂| 国产人伦9x9x在线观看| 性高湖久久久久久久久免费观看| 国产 精品1| 99热国产这里只有精品6| 欧美97在线视频| 最近中文字幕高清免费大全6| 免费女性裸体啪啪无遮挡网站| 青春草国产在线视频| 菩萨蛮人人尽说江南好唐韦庄| 成人午夜精彩视频在线观看| 亚洲,一卡二卡三卡| 满18在线观看网站| 18禁动态无遮挡网站| 啦啦啦在线观看免费高清www| 高清视频免费观看一区二区| 欧美人与性动交α欧美软件| 18在线观看网站| 人人妻,人人澡人人爽秒播 | 欧美日韩视频精品一区| tube8黄色片| 久久久久久久久免费视频了| 国产精品av久久久久免费| 秋霞在线观看毛片| 亚洲av男天堂| 久久av网站| 高清欧美精品videossex| 欧美日韩成人在线一区二区| 午夜福利视频在线观看免费| 国产成人av激情在线播放| 美女脱内裤让男人舔精品视频| 国产日韩欧美视频二区| 99久久精品国产亚洲精品| 亚洲精品一二三| 亚洲综合精品二区| 国产精品成人在线| videos熟女内射| 中文精品一卡2卡3卡4更新| 久久亚洲国产成人精品v| 好男人视频免费观看在线| 51午夜福利影视在线观看| 亚洲精品在线美女| 亚洲国产欧美在线一区| 九草在线视频观看| 日韩 亚洲 欧美在线| 精品久久久久久电影网| 另类亚洲欧美激情| 精品久久久精品久久久| 国产精品无大码| 少妇精品久久久久久久| 看免费av毛片| 亚洲国产精品一区二区三区在线| 国产欧美日韩一区二区三区在线| svipshipincom国产片| 日韩熟女老妇一区二区性免费视频| 国产精品免费视频内射| 国产成人免费观看mmmm| 国产精品二区激情视频| 一本大道久久a久久精品| 国产黄频视频在线观看| 国产xxxxx性猛交| 成人漫画全彩无遮挡| 老汉色av国产亚洲站长工具| 丰满少妇做爰视频| 亚洲美女黄色视频免费看| 美女午夜性视频免费| 午夜福利视频精品| 欧美亚洲 丝袜 人妻 在线| 欧美精品高潮呻吟av久久| 亚洲色图综合在线观看| 国产精品 国内视频| 国产日韩欧美亚洲二区| 欧美久久黑人一区二区| 成人毛片60女人毛片免费| 亚洲av成人不卡在线观看播放网 | 久久久久视频综合| 19禁男女啪啪无遮挡网站| 一级a爱视频在线免费观看| 十分钟在线观看高清视频www| 亚洲精品美女久久久久99蜜臀 | 久久毛片免费看一区二区三区| av电影中文网址| 婷婷色综合大香蕉| 国产午夜精品一二区理论片| 我要看黄色一级片免费的| 操出白浆在线播放| 国产精品国产av在线观看| 午夜福利影视在线免费观看| 别揉我奶头~嗯~啊~动态视频 | 国产免费又黄又爽又色| 中文字幕高清在线视频| 街头女战士在线观看网站| 久久精品熟女亚洲av麻豆精品| 99久久综合免费| 久久99精品国语久久久| 欧美97在线视频| netflix在线观看网站| 亚洲美女黄色视频免费看| 久热这里只有精品99| 又黄又粗又硬又大视频| 中文字幕人妻熟女乱码| av视频免费观看在线观看| 老司机深夜福利视频在线观看 | 午夜av观看不卡| 日韩欧美一区视频在线观看|