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

    基于Lanczos模態(tài)分析及密集度分析的發(fā)射裝置輕量化研究

    2016-07-16 03:26:47張英琦樂貴高馬大為馮國銅
    彈道學(xué)報 2016年2期
    關(guān)鍵詞:輕量化設(shè)計

    張英琦,樂貴高,馬大為,馮國銅

    (南京理工大學(xué) 機(jī)械工程學(xué)院,南京 210094)

    ?

    基于Lanczos模態(tài)分析及密集度分析的發(fā)射裝置輕量化研究

    張英琦,樂貴高,馬大為,馮國銅

    (南京理工大學(xué) 機(jī)械工程學(xué)院,南京 210094)

    摘要:針對某型號多管火箭發(fā)射裝置輕量化設(shè)計,采用有限元法建立其剛?cè)狁詈蟿恿W(xué)模型,對輕量化前后的模態(tài)響應(yīng)特性及其火箭彈發(fā)射密集度進(jìn)行對比分析,驗證了此輕量化設(shè)計的可行性。首先,對Lanczos迭代法在經(jīng)典的特征值問題中進(jìn)行了推導(dǎo)求解;其次,建立了某復(fù)合材料定向管的有限元模型,計算結(jié)果與此定向管樣機(jī)模態(tài)試驗結(jié)果進(jìn)行對比分析,結(jié)果表明基于Lanczos迭代法的有限元模型求解模態(tài)是正確的;最后,根據(jù)輕量化設(shè)計要求,建立了2種材料火箭炮發(fā)射裝置的整車有限元模型,采用Lanczos迭代法計算其模態(tài)特性,根據(jù)模態(tài)特性選擇合理的射擊間隔,并通過對比輕量化前后此發(fā)射裝置火箭彈的射擊密集度,驗證了此輕量化設(shè)計的可行性。

    關(guān)鍵詞:多管火箭炮;模態(tài)試驗;Lanczos迭代法;射擊密集度;輕量化設(shè)計

    多管火箭炮作為面殺傷武器,具有火力猛、威力大、機(jī)動靈活、反應(yīng)迅速等特點,它能夠在極短的時間內(nèi)發(fā)射多管火箭炮,給命中目標(biāo)以毀滅性打擊,給敵方心理上造成極大的震憾和威懾作用,在戰(zhàn)場上發(fā)揮著極其重要的作用[1-2]。輕量化設(shè)計是多管火箭炮發(fā)射裝置的重要發(fā)展趨勢。在確保不降低戰(zhàn)術(shù)技術(shù)性能的前提下對武器系統(tǒng)進(jìn)行輕量化設(shè)計,是提高作戰(zhàn)機(jī)動性的主要途徑。采用輕金屬、新型復(fù)合材料是實現(xiàn)多管火箭炮輕量化的主要措施之一。

    模態(tài)特征能夠從質(zhì)量和剛度特性上反映結(jié)構(gòu)的力學(xué)性能,故多管火箭炮發(fā)射裝置輕量化設(shè)計前后的模態(tài)特性是其動力學(xué)分析的重要基礎(chǔ),對控制和優(yōu)化火箭炮發(fā)射裝置的動態(tài)特性有著至關(guān)重要的作用。目前,我國在振動工程領(lǐng)域取得了很多研究成果,部分理論和實驗研究已達(dá)到國際先進(jìn)水平。但是,針對多管火箭炮結(jié)構(gòu)動力學(xué)的研究仍有待加強(qiáng)。

    本文通過模態(tài)分析反映多管火箭炮發(fā)射裝置的振動特性,從而依據(jù)振動特性調(diào)節(jié)多管火箭炮的合理發(fā)射間隔,在滿足發(fā)射密集度指標(biāo)的情況下,驗證輕量化設(shè)計的可行性。

    1模態(tài)求解基本理論

    大型結(jié)構(gòu)的特征值解法[3],主要有矩陣反迭代法、子空間迭代法、里茲向量直接疊加法和Lanczos迭代法等。Lanczos迭代法對n×n階矩陣進(jìn)行了數(shù)次迭代,可使矩陣三對角化,然后用Sturm定理和二分法求特征值,相應(yīng)的特征向量可用帶移位的反迭代求得。Lanczos迭代法可有效地克服因系統(tǒng)中存在等固有頻率或幾個固有頻率而引起的收斂過慢問題[4]。

    系統(tǒng)運動方程:

    (1)

    典型的無阻尼模態(tài)分析求解的基本方程就是經(jīng)典的特征值問題:

    Kφ-ω2Mφ=0

    (2)

    式中:φ為系統(tǒng)振型,ω為系統(tǒng)自振頻率。

    ω2=λ-δ

    (3)

    式中:δ>0,則λ即為正數(shù),這樣就將零特征值問題化解成非零特征值問題:

    Kφ=(λ-δ)Mφ

    (4)

    (5)

    M=LLT

    (6)

    (7)

    y=LTφ

    (8)

    則式(4)化為

    Ay=λy

    (9)

    尋找正交矩陣Q=(q1q2…qn),使得QTAQ=T,式中:n為正交基個數(shù),T為三對角陣,記為

    Lanczos算法遞推求解得到(當(dāng)j=1時,取βj-1=0):

    ①αj=(qjAqj)。

    ②殘余向量wj=Aqj-βj-1qj-1-αjqj。

    ③βj=‖wj‖2,qj+1=wj/βj。

    式中:j為遞推次數(shù)。上述遞推過程循環(huán)計算,可求得標(biāo)準(zhǔn)正交基(q1q2…qn)。采用Lanczos法求得正交變換矩陣Q和對角矩陣T后,可用二分法求解λ和ω。

    2復(fù)合材料定向管模態(tài)分析與試驗

    2.1復(fù)合材料定向管有限元模型

    根據(jù)某復(fù)合材料定向管的結(jié)構(gòu)和材料,模擬其鋪層形式,用殼單元建立有限元模型,如圖1所示。

    圖1復(fù)合材料定向管有限元模型

    復(fù)合材料定向管徑向由10層玻璃纖維增強(qiáng)樹脂基復(fù)合材料環(huán)向鋪層纏繞而成,相鄰兩鋪層之間相互垂直,以45°或-45°斜角鋪層。每一層的厚度為0.3mm,3mm厚的玻璃鋼定向管劃分后的單元共10層,如圖2所示。

    圖2 定向管材料鋪層

    材料參數(shù)利用文獻(xiàn)[5]中的玻璃纖維增強(qiáng)復(fù)合材料的性能參數(shù),如表1所示。表中,ρ為材料密度,E1為縱向彈性模量,E2為橫向彈性模量,ν為泊松比,G12為縱橫剪切彈性模量,Xt為縱向拉伸強(qiáng)度,Yt為橫向拉伸強(qiáng)度,Xc為縱向壓縮強(qiáng)度,Yc為橫向壓縮強(qiáng)度,S為剪切強(qiáng)度。

    表1 復(fù)合材料定向管材料參數(shù)

    2.2定向管樣機(jī)模態(tài)試驗驗證

    為驗證Lanczos迭代法在有限元模態(tài)分析中的準(zhǔn)確性,建立了復(fù)合材料定向管有限元模型,采用Lanczos迭代法進(jìn)行模態(tài)分析,同時制定了該防空火箭炮定向管樣機(jī)并進(jìn)行模態(tài)試驗。圖3和圖4為該復(fù)合材料定向管試驗?zāi)P秃湍B(tài)試驗測試系統(tǒng);圖5為錘擊法模態(tài)試驗原理框架圖;對該復(fù)合材料定向管仿真結(jié)果與試驗結(jié)果的前三階模態(tài)進(jìn)行對比如圖6所示。試驗過程中力錘在每個測試點上敲擊3次對其結(jié)果進(jìn)行平均,對每個點均觀察波形和試驗自譜,保證相干函數(shù)值不低于0.8。

    圖3 定向管試驗幾何模型

    圖4 模態(tài)試驗測試系統(tǒng)

    圖5 錘擊法模態(tài)試驗框圖

    如圖6和表2所示,d為模態(tài)階數(shù),fe為計算頻率,ft為試驗頻率,ε為結(jié)果偏差。對比定向管前三階自由模態(tài)的仿真結(jié)果與試驗結(jié)果可知,定向管振型一致,振動頻率誤差均在2%以內(nèi),兩者吻合度較高,考慮試驗測試存在不可控因素及數(shù)值解自身的誤差,兩者的吻合度是令人滿意的,表明了有限元模型的正確性??梢?基于Lanczos迭代法的有限元模型求解模態(tài)是可靠的,可以為后續(xù)的動力學(xué)分析提供研究基礎(chǔ)。

    圖6 仿真結(jié)果與試驗結(jié)果對比圖

    dfe/Hzft/Hzε/%191.61490.021.772270.94269.960.363607.47619.461.94

    3某多管火箭炮發(fā)射裝置輕量化設(shè)計

    3.1某車載多管火箭炮發(fā)射裝置數(shù)值模型

    發(fā)射裝置起豎仰角53°,回轉(zhuǎn)角度為0,采用輪胎和支腿的混合支撐方式,按給定射序發(fā)射6枚火箭彈。如圖7所示,建立了某型號車載多管火箭炮發(fā)射裝置的有限元模型及其射序。

    圖7 車載多管火箭炮有限元模型和射序示意圖

    仿真過程中的假設(shè)條件:

    ①發(fā)射車所處場坪為水平狀態(tài),忽略場坪的變化;

    ②忽略環(huán)境因素對發(fā)射裝置動力學(xué)的影響;

    ③忽略對發(fā)射平臺動力學(xué)影響較小的部件結(jié)構(gòu),如駕駛室、儀器艙、油箱、變速箱、發(fā)動機(jī)等部件,按照給定質(zhì)量參數(shù)對整車進(jìn)行配重;

    ④認(rèn)為發(fā)射系統(tǒng)各結(jié)構(gòu)部件在火箭彈發(fā)射過程中不會塑性變形和斷裂,變形均處于彈性變形范圍內(nèi)。

    該多管火箭炮發(fā)射裝置的射擊間隔初始設(shè)計值為0.7s,射擊間隔與系統(tǒng)的模態(tài)響應(yīng)密切相關(guān),應(yīng)盡量避開發(fā)射裝置的固有振動頻率[6-8],避免火箭炮發(fā)射過程中產(chǎn)生共振影響射擊密集度。因此,采用2.2節(jié)中驗證的Lanczos迭代法對發(fā)射裝置進(jìn)行模態(tài)仿真分析。給出該發(fā)射裝置的前四階自由模態(tài)如圖8所示,發(fā)射裝置前四階模態(tài)頻率及周期如表3所示,最大的第一階周期為0.619 6s,小于射擊間隔初始設(shè)計值,表明了該設(shè)計值的可行性。

    圖8 鋼結(jié)構(gòu)發(fā)射裝置前四階自由模態(tài)

    df/HzT1/sdf/HzT1/s11.61390.619634.04290.247323.64670.274244.69240.2131

    在0.7s發(fā)射間隔基礎(chǔ)上,通過仿真計算得到6枚火箭彈發(fā)射過程中的角速度、角位移響應(yīng)曲線,限于篇幅,給出6枚火箭彈的俯仰角速度響應(yīng)曲線如圖9所示,輕量化前火箭彈的起始擾動仿真結(jié)果如表4所示,表中,t為出筒時刻,ω1為俯仰角速度,φ1為俯仰角位移,ω2為偏航角速度,φ2為偏航角位移。

    圖9 輕量化前火箭炮俯仰角速度響應(yīng)

    t/sω1/(rad·s-1)φ1/radω2/(rad·s-1)φ2/rad0.20.001950-0.0001860.0061740.0002090.9-0.006202-0.0005970.0176180.0004131.6-0.0009190.001614-0.008548-0.0002262.3-0.010984-0.0007380.0111730.0002843.0-0.003252-0.000714-0.0003600.0001503.7-0.1117730.0002330.0058640.000183

    3.2車載多管火箭炮發(fā)射裝置輕量化設(shè)計

    3.1節(jié)中所示的車載多管火箭炮發(fā)射裝置,以其輕量化設(shè)計要求為背景,將發(fā)射箱、上架、下架和回轉(zhuǎn)盤等強(qiáng)度要求不高的部件均由鋼材料輕量化為鋁合金材料,即為發(fā)射裝置的輕量化初步設(shè)計方案。鋁合金材料密度為2 700kg/m3,彈性模量取70GPa,泊松比為0.33。以下對該輕量化設(shè)計方案從射擊密集度和射擊間隔兩方面進(jìn)行可行性驗證。

    采用2.2節(jié)中驗證的Lanczos迭代法對輕量化后的發(fā)射裝置進(jìn)行模態(tài)分析。給出輕量化后發(fā)射裝置的前四階自由模態(tài)如圖10所示,輕量化后發(fā)射裝置前四階模態(tài)頻率及周期如表5所示,輕量化后最大的第一階周期為0.811 0s,大于0.7s的射擊間隔設(shè)計值。研究發(fā)現(xiàn),發(fā)射裝置輕量化后,系統(tǒng)固有特性發(fā)生了變化,需要重新確定合理的發(fā)射間隔,如發(fā)射間隔沿用0.7s,連續(xù)發(fā)射時,前一發(fā)彈對后一發(fā)彈的影響明顯增大,使后一發(fā)火箭彈離軌時的起始擾動增大,使得密集度變差,達(dá)不到設(shè)計指標(biāo)要求。因此,需要重新確定發(fā)射時間間隔,結(jié)合結(jié)構(gòu)模態(tài)分析結(jié)果,并對比0.9s,1s,1.1s和1.2s的不同發(fā)射間隔對起始擾動的影響,確定了1.2s的發(fā)射時間間隔,并通過密集度計算對新發(fā)射時間間隔進(jìn)行了驗證。

    圖10 鋁合金材料發(fā)射裝置前四階自由模態(tài)

    df/HzT1/sdf/HzT1/s11.23300.811032.88620.346522.18120.458543.13720.3188

    將3.1節(jié)中數(shù)值模型的發(fā)射箱、上下架和回轉(zhuǎn)盤等部件賦予鋁合金材料屬性,仿真計算該裝置6枚火箭彈出筒的角速度、角位移響應(yīng)曲線,限于篇幅僅給出6枚火箭彈的俯仰角速度曲線如圖11所示,輕量化后火箭彈的起始擾動仿真結(jié)果如表6所示。

    圖11 輕量化后火箭炮俯仰角速度響應(yīng)

    t/sω1/(rad·s-1)φ1/radω2/(rad·s-1)φ2/rad0.20.012813-0.0006820.0133670.0004321.4-0.014839-0.0012890.0047530.0002982.60.021244-0.001189-0.0015990.0002723.8-0.010272-0.0013320.0269120.00066350.010787-0.0009730.0184110.0005886.20.007925-0.001304-0.0098160.000119

    3.3發(fā)射裝置輕量化前后密集度計算

    多管火箭炮的射擊密集度是對火箭彈散布的度量,是火箭炮武器系統(tǒng)的重要技術(shù)性能指標(biāo),是射擊精度的重要組成部分。故發(fā)射裝置的輕量化設(shè)計應(yīng)該在穩(wěn)定火箭炮武器發(fā)射密集度的基礎(chǔ)上進(jìn)行。

    運用3.1節(jié)和3.2節(jié)中輕量化前后火箭彈起始擾動仿真結(jié)果,并通過文獻(xiàn)[9],以六自由度剛體彈道模型為基礎(chǔ),綜合考慮質(zhì)量、質(zhì)量偏心等因素對火箭彈飛行過程及密集度的影響,建立了一般無控火箭彈外彈道模型。采用文獻(xiàn)[10~11]中密集度算法編程計算,獲得了該發(fā)射裝置輕量化前后的射擊密集度,輕量化前發(fā)射裝置偏航密集度為1/201rad,距離密集度為1/135rad,輕量化后發(fā)射裝置偏航密集度為1/190rad,距離密集度為1/110rad。

    輕量化未導(dǎo)致射擊密集度的降低,且有小幅度提升,滿足最大射程密集度設(shè)計指標(biāo)要求。通過對射擊密集度和射擊間隔指標(biāo)的分析,可知文中給出的火箭炮發(fā)射裝置的輕量化設(shè)計是可行的。改變發(fā)射間隔,在輕量化的基礎(chǔ)上密集度也有相應(yīng)的提高。

    4結(jié)束語

    定向管前三階自由模態(tài)的仿真結(jié)果與試驗吻合度較高,驗證了基于Lanczos迭代法的有限元模型求解模態(tài)的正確性。

    進(jìn)行輕量化設(shè)計后,通過有限元模態(tài)分析,發(fā)射裝置的最大的第一階周期為0.811 0s,大于0.7s的射擊間隔設(shè)計值,考慮發(fā)射過程中火箭彈的起始擾動,將輕量化后的發(fā)射間隔改為1.2s;進(jìn)行輕量化設(shè)計后,通過改變射擊時間間隔,發(fā)射裝置的密集度有一定的提高。

    本文提出的輕量化設(shè)計方案能夠有效地降低發(fā)射裝置的質(zhì)量,通過改變射擊時間間隔,減小火箭彈初始擾動,以此提高了射擊密集度。

    參考文獻(xiàn)

    [1]張敬.火箭炮的發(fā)展及性能提高的途徑[J].導(dǎo)彈與航天運載技術(shù),2001(2):58-61.

    ZHANG Jing.Development of rocket gun and methods of improving performance[J].Missiles and Space Vehicles,2001(2):58-61.(in Chinese)

    [2]CHRISTENSEN D E.Multiple rocket launcher characteristics and simulation technique,AD-A025379[R].1976.

    [3]王勖成.有限單元法[M].北京:清華大學(xué)出版社,2003.

    WANG Xu-cheng.Finite element method[M].Beijing:Tsinghua University Press,2003.(in Chinese)

    [4]徐稼軒,鄭鐵生.結(jié)構(gòu)動力分析的數(shù)值方法[M].西安:西安交通大學(xué)出版社,1993.

    XU Jia-xuan,ZHENG Tie-sheng.Numerical method for dynamic analysis of structures[M].Xi’an:Xi’an Jiaotong University Press,1993.(in Chinese)

    [5]蔡德詠,馬大為,于存貴,等.復(fù)合材料定向管輕量化設(shè)計[J].南京理工大學(xué)學(xué)報,2011,35(2):168-172.

    CAI De-yong,MA Da-wei,YU Cun-gui,et al.Light-weight design of composite materials direction tube[J].Journal of Nanjing University of Science and Technology,2011,35(2):168-172.(in Chinese)

    [6]潘宏俠,趙彥,陳國光,等.振動對火箭炮密集度影響研究[J].兵工學(xué)報,1993,11(4):17-22.

    PAN Hong-xia,ZHAO Yan,CHEN Guo-guang,et al.Effects of launcher vibration upon the hitting dispersion of rockets[J].Acta Armamentarii,1993,11(4):17-22.(in Chinese)

    [7]賀軍義,芮筱亭,王國平,等.提高多管火箭射擊密集度射擊技術(shù)研究[J].振動工程學(xué)報,2011,24(6):676-681.

    HE Jun-yi,RUI Xiao-ting,WANG Guo-ping,et al.Design technology of improving firing dispersion of MLRS[J].Journal of Vibration Engineering,2011,24(6):676-681.(in Chinese)

    [8]楊帆,王國平,芮筱亭,等.射擊間隔對多管火箭動態(tài)特性的影響[J].南京理工大學(xué)學(xué)報,2006,30(4):400-403.

    YANG Fan,WANG Guo-ping,RUI Xiao-ting,et al.Influence of firing interval on dynamics characteristic of MLRS[J].Journal of Nanjing University of Science and Technology,2006,30(4):400-403.(in Chinese)

    [9]徐明友.火箭外彈道學(xué)[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,2004.

    XU Ming-you.Exterior ballistics of rockets[M].Harbin:Harbin Institute of Technology Press,2004.(in Chinese)

    [10]王峰,于存貴,馬大為,等.艦載多功能火箭炮射擊密集度仿真分析[J].彈箭與制導(dǎo)學(xué)報,2006,26(2):737-739.

    WANG Feng,YU Cun-gui,MA Da-wei,et al.Simulation research on shooting dispersion of ship-borne multi-function launcher rocket[J].Journal of Projectiles,Rockets,Missiles and Guidance,2006,26(2):737-739.(in Chinese)

    [11]李魁武,王寶元.火炮射擊密集度研究方法[M].北京:國防工業(yè)出版社,2012:24-29.

    LI Kui-wu,WANG Bao-yuan.Research methods of gun firing dispersion[M].Beijing:National Defense Industry Press,2012:24-29.(in Chinese)

    Lightweight Research of MLRS Based on Modal Analysis With Lanczos-reduce Method and Dispersion Analysis

    ZHANG Ying-qi,LE Gui-gao,MA Da-wei,FENG Guo-tong

    (School of Mechanical Engineering,Nanjing University of Science and Technology,Nanjing 210094,China)

    Abstract:To research the lightweight of multiple launcher rocket system(MLRS),the rigid-flexible-coupling dynamic model with finite element method was applied.The modal response and firing dispersion after the lightweight was compared with that before the lightweight respectively,and the feasibility of lightweight was verified.The classical eigenvalue problem was solved by Lanczos-reduce method.The FEM model of composite materials direction tube was built,and the FEM results were compared with the modal test results.The result shows that the modal solved by Lanczos is correct. The FEM models of MLRS of two kinds of materials were built,and the modal characteristics were calculated by using Lanczos.The reasonable firing interval was selected according to modal characteristics.The firing dispersion of rocket after lightweight was compared with that before lightweight.The result shows that the design of lightweight is feasible.

    Key words:multiple launcher rocket system;modal test;Lanczos-reduce method;firing dispersion;lightweight design

    收稿日期:2015-12-18

    作者簡介:張英琦(1988- ),女,博士研究生,研究方向為結(jié)構(gòu)動力學(xué)。E-mail:zhangyingqi0627@163.com。

    中圖分類號:TJ303.4

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

    文章編號:1004-499X(2016)02-0068-06

    猜你喜歡
    輕量化設(shè)計
    FSC方程式賽車單體殼車架設(shè)計
    巴哈賽車車架仿真分析與優(yōu)化設(shè)計
    關(guān)于岸邊集裝箱起重機(jī)結(jié)構(gòu)輕量化設(shè)計及可靠性分析
    大型客機(jī)隔熱隔音棉的輕量化設(shè)計
    科技視界(2017年9期)2017-09-04 19:14:07
    伊朗設(shè)拉子地鐵車輛輕量化碳鋼車體設(shè)計
    科技視界(2017年2期)2017-04-18 09:07:27
    軌道扣件彈性墊板結(jié)構(gòu)優(yōu)化設(shè)計
    乘用車車身結(jié)構(gòu)輕量化設(shè)計技術(shù)研究與實踐
    FSAE賽車車架結(jié)構(gòu)拓?fù)鋬?yōu)化及輕量化設(shè)計研究
    帽形截面薄壁焊接結(jié)構(gòu)軸向耐撞性的優(yōu)化設(shè)計方法研究
    飛機(jī)平尾普通肋的優(yōu)化設(shè)計
    一区二区三区免费毛片| 少妇高潮的动态图| 91久久精品电影网| 国产精品人妻久久久影院| 十八禁高潮呻吟视频 | 国产综合精华液| 老司机影院成人| 国产乱人偷精品视频| 日本黄色日本黄色录像| 午夜av观看不卡| 九色成人免费人妻av| 亚洲精品国产av成人精品| 成年av动漫网址| 99re6热这里在线精品视频| 成人毛片60女人毛片免费| 桃花免费在线播放| 高清欧美精品videossex| 欧美成人午夜免费资源| 色网站视频免费| 亚洲欧洲国产日韩| 精品卡一卡二卡四卡免费| 另类精品久久| 色5月婷婷丁香| av.在线天堂| 精品久久久久久电影网| 王馨瑶露胸无遮挡在线观看| 精品亚洲成国产av| 亚洲欧美一区二区三区黑人 | 久久婷婷青草| 婷婷色麻豆天堂久久| 九色成人免费人妻av| 国产精品嫩草影院av在线观看| videos熟女内射| 夫妻性生交免费视频一级片| 亚洲精品色激情综合| 美女内射精品一级片tv| 久久久国产欧美日韩av| 在线看a的网站| 亚洲国产av新网站| 毛片一级片免费看久久久久| 91成人精品电影| 久久久久久久久久人人人人人人| 免费不卡的大黄色大毛片视频在线观看| 美女内射精品一级片tv| 久久久久久久亚洲中文字幕| 我的老师免费观看完整版| 日韩三级伦理在线观看| 99久久精品国产国产毛片| 在线观看人妻少妇| 亚洲精品aⅴ在线观看| 亚洲av福利一区| 夜夜骑夜夜射夜夜干| 又爽又黄a免费视频| 人人妻人人澡人人爽人人夜夜| 最黄视频免费看| 久久青草综合色| 丝袜在线中文字幕| 黄色配什么色好看| 亚洲av男天堂| 亚洲国产精品成人久久小说| 中国美白少妇内射xxxbb| 2022亚洲国产成人精品| av国产精品久久久久影院| 欧美3d第一页| 亚洲av二区三区四区| 日本午夜av视频| 91在线精品国自产拍蜜月| 如何舔出高潮| 99视频精品全部免费 在线| 美女国产视频在线观看| videossex国产| 十八禁高潮呻吟视频 | 亚洲欧美清纯卡通| 国产伦精品一区二区三区四那| 22中文网久久字幕| 成人漫画全彩无遮挡| 日本黄色片子视频| 久久av网站| 国产精品福利在线免费观看| 国产成人a∨麻豆精品| 亚洲图色成人| 在线天堂最新版资源| av又黄又爽大尺度在线免费看| 日本爱情动作片www.在线观看| 国内精品宾馆在线| a级毛色黄片| 狠狠精品人妻久久久久久综合| 国产精品蜜桃在线观看| 日日啪夜夜爽| 久久婷婷青草| 在线播放无遮挡| 日韩熟女老妇一区二区性免费视频| 人妻夜夜爽99麻豆av| 亚洲天堂av无毛| 自拍欧美九色日韩亚洲蝌蚪91 | 大又大粗又爽又黄少妇毛片口| 午夜影院在线不卡| 久久av网站| 精品国产一区二区三区久久久樱花| 婷婷色麻豆天堂久久| 最黄视频免费看| 久久国产乱子免费精品| 国产精品秋霞免费鲁丝片| 嫩草影院入口| 亚洲第一av免费看| 只有这里有精品99| 日韩欧美精品免费久久| 亚洲三级黄色毛片| 欧美日韩视频精品一区| 国产在视频线精品| 久久久欧美国产精品| 日本与韩国留学比较| 校园人妻丝袜中文字幕| 亚洲国产毛片av蜜桃av| 亚洲精品,欧美精品| 天美传媒精品一区二区| 久久久a久久爽久久v久久| 欧美精品一区二区大全| 国产男女内射视频| 国产在视频线精品| 欧美3d第一页| 精品亚洲乱码少妇综合久久| 观看免费一级毛片| 伊人久久精品亚洲午夜| 免费高清在线观看视频在线观看| 极品人妻少妇av视频| 搡老乐熟女国产| 国产免费一区二区三区四区乱码| 亚洲高清免费不卡视频| 一区二区三区免费毛片| 国产淫片久久久久久久久| 国产乱人偷精品视频| 午夜福利在线观看免费完整高清在| 大片电影免费在线观看免费| 国产一级毛片在线| 国产精品久久久久久久电影| 777米奇影视久久| 国产美女午夜福利| 少妇高潮的动态图| 多毛熟女@视频| 国产亚洲一区二区精品| 国产精品.久久久| 99热网站在线观看| 国产日韩欧美视频二区| 岛国毛片在线播放| 色哟哟·www| 国产69精品久久久久777片| 国产精品.久久久| 国产欧美另类精品又又久久亚洲欧美| 男人狂女人下面高潮的视频| 亚洲成色77777| av又黄又爽大尺度在线免费看| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 啦啦啦在线观看免费高清www| 99热6这里只有精品| 久久鲁丝午夜福利片| 亚洲av.av天堂| 国产精品一区www在线观看| 久久av网站| 一区二区av电影网| 少妇熟女欧美另类| 日本免费在线观看一区| 精品少妇久久久久久888优播| 日日啪夜夜爽| 黄片无遮挡物在线观看| 2021少妇久久久久久久久久久| 亚洲人成网站在线观看播放| 日日撸夜夜添| 99热这里只有是精品在线观看| 黄色一级大片看看| 97在线视频观看| 欧美日韩亚洲高清精品| 男的添女的下面高潮视频| 校园人妻丝袜中文字幕| 在线观看免费视频网站a站| 国产精品无大码| 成人无遮挡网站| 一级毛片aaaaaa免费看小| av播播在线观看一区| 国模一区二区三区四区视频| 人妻 亚洲 视频| 五月伊人婷婷丁香| 女性被躁到高潮视频| 中文资源天堂在线| 国产男人的电影天堂91| 自线自在国产av| 午夜福利网站1000一区二区三区| 美女福利国产在线| 一级毛片aaaaaa免费看小| 哪个播放器可以免费观看大片| 国产在视频线精品| 免费大片黄手机在线观看| 大码成人一级视频| 国产爽快片一区二区三区| 亚洲欧美一区二区三区黑人 | 欧美日本中文国产一区发布| 有码 亚洲区| 如何舔出高潮| 亚洲一区二区三区欧美精品| 一本大道久久a久久精品| 国产精品不卡视频一区二区| 精品午夜福利在线看| 丰满饥渴人妻一区二区三| 黑丝袜美女国产一区| 伦理电影免费视频| 午夜影院在线不卡| 搡女人真爽免费视频火全软件| 国产精品久久久久久精品电影小说| 十八禁网站网址无遮挡 | 国产伦理片在线播放av一区| 欧美精品国产亚洲| 99久久人妻综合| 男女边摸边吃奶| 大片免费播放器 马上看| 国产欧美日韩精品一区二区| 日日撸夜夜添| 亚洲人成网站在线播| 精品人妻偷拍中文字幕| 春色校园在线视频观看| 少妇高潮的动态图| 国产综合精华液| 国产成人精品一,二区| 91aial.com中文字幕在线观看| 我要看黄色一级片免费的| 欧美激情国产日韩精品一区| 久久久精品94久久精品| 99热这里只有是精品50| 少妇人妻 视频| 狂野欧美白嫩少妇大欣赏| 黄色日韩在线| 在线观看美女被高潮喷水网站| 亚洲国产毛片av蜜桃av| 亚洲国产av新网站| 晚上一个人看的免费电影| 亚洲精品aⅴ在线观看| 午夜福利视频精品| 免费看不卡的av| 国产成人91sexporn| 99热6这里只有精品| 插逼视频在线观看| 久久韩国三级中文字幕| 欧美+日韩+精品| 国产亚洲91精品色在线| 涩涩av久久男人的天堂| 少妇人妻 视频| 日韩一区二区三区影片| 久久热精品热| 国产老妇伦熟女老妇高清| 国产欧美亚洲国产| 9色porny在线观看| 成人午夜精彩视频在线观看| 亚洲自偷自拍三级| 精品亚洲成a人片在线观看| 中文欧美无线码| 女性被躁到高潮视频| 国产精品.久久久| 久久久久人妻精品一区果冻| 久久婷婷青草| 在线观看三级黄色| av在线app专区| 91aial.com中文字幕在线观看| 99久久综合免费| 女的被弄到高潮叫床怎么办| 人人妻人人爽人人添夜夜欢视频 | av一本久久久久| 国产在线一区二区三区精| av播播在线观看一区| av女优亚洲男人天堂| 国产男女超爽视频在线观看| 亚洲,一卡二卡三卡| 在线观看一区二区三区激情| 国产高清不卡午夜福利| 久久久久精品久久久久真实原创| 久久久国产一区二区| 在线 av 中文字幕| 国产深夜福利视频在线观看| 亚洲在久久综合| 只有这里有精品99| 综合色丁香网| 亚洲精品乱码久久久v下载方式| 女性生殖器流出的白浆| 久久久久久久久大av| 99热全是精品| 91精品一卡2卡3卡4卡| 午夜av观看不卡| 激情五月婷婷亚洲| 99九九在线精品视频 | 亚洲精品色激情综合| 欧美3d第一页| 国产有黄有色有爽视频| 日日撸夜夜添| 亚洲国产色片| 久久av网站| 搡女人真爽免费视频火全软件| 性色avwww在线观看| 欧美精品高潮呻吟av久久| 日本黄大片高清| 欧美国产精品一级二级三级 | 国产黄色免费在线视频| 91在线精品国自产拍蜜月| 精品亚洲成a人片在线观看| 高清午夜精品一区二区三区| 欧美bdsm另类| 午夜福利视频精品| 99热这里只有是精品50| 日韩大片免费观看网站| 91成人精品电影| 少妇人妻 视频| 91精品国产九色| 日日啪夜夜撸| 99热全是精品| 2022亚洲国产成人精品| 国产精品一区www在线观看| 一级毛片黄色毛片免费观看视频| 亚洲国产最新在线播放| 精品国产一区二区三区久久久樱花| 久久午夜综合久久蜜桃| 亚洲精品乱码久久久v下载方式| 亚洲国产精品一区二区三区在线| 国产欧美另类精品又又久久亚洲欧美| 男人和女人高潮做爰伦理| 成人午夜精彩视频在线观看| 2022亚洲国产成人精品| 在线观看免费高清a一片| 亚洲av成人精品一区久久| 久久99蜜桃精品久久| 性色avwww在线观看| a级一级毛片免费在线观看| 日本色播在线视频| 男男h啪啪无遮挡| 蜜桃久久精品国产亚洲av| 国产av精品麻豆| 少妇猛男粗大的猛烈进出视频| 99re6热这里在线精品视频| 久久影院123| 日韩电影二区| 国产精品免费大片| 日韩一区二区视频免费看| xxx大片免费视频| videos熟女内射| www.色视频.com| 男人狂女人下面高潮的视频| 国产有黄有色有爽视频| 亚洲中文av在线| 国产伦在线观看视频一区| www.色视频.com| 下体分泌物呈黄色| 久久女婷五月综合色啪小说| 国产精品女同一区二区软件| 久久久久久久亚洲中文字幕| 精品国产乱码久久久久久小说| 最近2019中文字幕mv第一页| 色网站视频免费| 日韩av免费高清视频| 亚洲不卡免费看| 日日啪夜夜撸| 成年人午夜在线观看视频| 九草在线视频观看| 丰满人妻一区二区三区视频av| 国产一级毛片在线| 亚洲精品乱码久久久久久按摩| 国产亚洲精品久久久com| 久久久久国产网址| 少妇的逼好多水| 中文在线观看免费www的网站| 最新的欧美精品一区二区| av国产久精品久网站免费入址| 亚洲国产精品999| 一级毛片 在线播放| 简卡轻食公司| 色94色欧美一区二区| 夫妻午夜视频| 婷婷色av中文字幕| 精品国产露脸久久av麻豆| 啦啦啦在线观看免费高清www| 91成人精品电影| 中国三级夫妇交换| 美女中出高潮动态图| 国产男女内射视频| 97超碰精品成人国产| 我要看日韩黄色一级片| 不卡视频在线观看欧美| 一二三四中文在线观看免费高清| 久久女婷五月综合色啪小说| 婷婷色综合www| av.在线天堂| 2021少妇久久久久久久久久久| 亚洲欧美一区二区三区黑人 | 九九在线视频观看精品| 欧美 亚洲 国产 日韩一| 精品熟女少妇av免费看| a 毛片基地| 日韩人妻高清精品专区| 色视频www国产| 成人黄色视频免费在线看| 极品教师在线视频| 国产一区二区在线观看av| 涩涩av久久男人的天堂| 久久久国产一区二区| 久热这里只有精品99| 噜噜噜噜噜久久久久久91| 免费看不卡的av| 国产精品女同一区二区软件| 亚洲中文av在线| 美女内射精品一级片tv| 永久网站在线| 在线免费观看不下载黄p国产| 最新的欧美精品一区二区| 在线观看www视频免费| av免费观看日本| 亚洲精品国产av蜜桃| 日韩制服骚丝袜av| 国产日韩欧美在线精品| 欧美国产精品一级二级三级 | 精品国产国语对白av| 人体艺术视频欧美日本| 国产淫片久久久久久久久| 亚洲精品自拍成人| 一级av片app| 久久久久精品久久久久真实原创| 免费黄频网站在线观看国产| 国产亚洲午夜精品一区二区久久| 国产伦理片在线播放av一区| 久久国产亚洲av麻豆专区| 最近中文字幕2019免费版| 国产亚洲精品久久久com| 中文天堂在线官网| 欧美丝袜亚洲另类| 日韩不卡一区二区三区视频在线| 一个人免费看片子| 免费久久久久久久精品成人欧美视频 | 特大巨黑吊av在线直播| 日本av免费视频播放| 三上悠亚av全集在线观看 | 午夜免费观看性视频| 成人影院久久| 少妇高潮的动态图| 99久国产av精品国产电影| 国产永久视频网站| 久久精品国产亚洲网站| 日韩精品有码人妻一区| 亚洲第一av免费看| 久久影院123| 看免费成人av毛片| 我要看黄色一级片免费的| 22中文网久久字幕| 美女视频免费永久观看网站| 国产成人91sexporn| 黄色欧美视频在线观看| 一级二级三级毛片免费看| av有码第一页| 国产精品嫩草影院av在线观看| 亚洲一区二区三区欧美精品| 99精国产麻豆久久婷婷| 国语对白做爰xxxⅹ性视频网站| 寂寞人妻少妇视频99o| 少妇人妻精品综合一区二区| 日韩伦理黄色片| 亚洲图色成人| 有码 亚洲区| av专区在线播放| 免费大片18禁| 欧美区成人在线视频| 欧美日韩av久久| 亚洲精品,欧美精品| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久影院| 免费看不卡的av| 高清不卡的av网站| 寂寞人妻少妇视频99o| 免费在线观看成人毛片| 日韩不卡一区二区三区视频在线| 精品一区二区三卡| 天美传媒精品一区二区| 亚洲av二区三区四区| 日本欧美国产在线视频| 国产精品福利在线免费观看| 在线亚洲精品国产二区图片欧美 | 亚洲欧美精品自产自拍| 亚洲av电影在线观看一区二区三区| av专区在线播放| 中文精品一卡2卡3卡4更新| av线在线观看网站| 视频区图区小说| 中文字幕制服av| freevideosex欧美| 韩国高清视频一区二区三区| 日本与韩国留学比较| 免费看av在线观看网站| 黄片无遮挡物在线观看| 建设人人有责人人尽责人人享有的| 亚洲精品亚洲一区二区| 51国产日韩欧美| 交换朋友夫妻互换小说| 色婷婷久久久亚洲欧美| 午夜福利在线观看免费完整高清在| 久久午夜福利片| 欧美激情国产日韩精品一区| 视频中文字幕在线观看| kizo精华| 97在线视频观看| 人人妻人人看人人澡| 啦啦啦视频在线资源免费观看| 国产成人精品久久久久久| 青春草国产在线视频| 少妇人妻 视频| 成年女人在线观看亚洲视频| 国产av国产精品国产| 寂寞人妻少妇视频99o| 欧美3d第一页| 久久久久视频综合| 午夜福利,免费看| 最后的刺客免费高清国语| 亚洲中文av在线| 久久青草综合色| 看免费成人av毛片| 在线观看免费高清a一片| 精品酒店卫生间| 最近中文字幕高清免费大全6| 国产成人一区二区在线| 亚洲精品久久午夜乱码| 国产免费福利视频在线观看| 国产一级毛片在线| 又粗又硬又长又爽又黄的视频| 久久久久久久久久成人| 欧美 亚洲 国产 日韩一| 黄片无遮挡物在线观看| 欧美精品高潮呻吟av久久| 免费大片18禁| 成年美女黄网站色视频大全免费 | 中文字幕久久专区| 久久久国产一区二区| 亚洲丝袜综合中文字幕| a级毛片免费高清观看在线播放| 两个人免费观看高清视频 | 国产一区有黄有色的免费视频| 夫妻午夜视频| 欧美最新免费一区二区三区| 久久人人爽人人片av| 美女国产视频在线观看| 蜜桃在线观看..| 国模一区二区三区四区视频| 只有这里有精品99| 国产一区有黄有色的免费视频| 天堂8中文在线网| 天堂中文最新版在线下载| 免费看av在线观看网站| 中文精品一卡2卡3卡4更新| 人妻制服诱惑在线中文字幕| 人妻少妇偷人精品九色| 亚洲欧美成人综合另类久久久| 午夜视频国产福利| 亚洲国产最新在线播放| 日本与韩国留学比较| 亚洲av成人精品一区久久| 日本爱情动作片www.在线观看| 人妻夜夜爽99麻豆av| 久久久久久久国产电影| 男女国产视频网站| 熟女电影av网| 免费看日本二区| 精品午夜福利在线看| 看免费成人av毛片| 亚洲人与动物交配视频| 亚洲国产精品专区欧美| 亚洲无线观看免费| 在现免费观看毛片| 国产欧美日韩综合在线一区二区 | 国产精品人妻久久久影院| 十八禁网站网址无遮挡 | 精华霜和精华液先用哪个| 精品一品国产午夜福利视频| 国产精品久久久久久av不卡| 免费高清在线观看视频在线观看| 极品少妇高潮喷水抽搐| 视频中文字幕在线观看| 男的添女的下面高潮视频| 乱人伦中国视频| 亚洲情色 制服丝袜| 亚洲国产精品一区二区三区在线| 久久久久久久国产电影| 噜噜噜噜噜久久久久久91| 亚洲电影在线观看av| 极品教师在线视频| 黑人猛操日本美女一级片| 精品国产一区二区三区久久久樱花| 欧美丝袜亚洲另类| 国产国拍精品亚洲av在线观看| 国产白丝娇喘喷水9色精品| 人妻系列 视频| 亚洲av福利一区| 伦理电影大哥的女人| 少妇熟女欧美另类| 精品午夜福利在线看| 国产黄片美女视频| 性色av一级| 九色成人免费人妻av| 色网站视频免费| 亚洲精品456在线播放app| 2018国产大陆天天弄谢| 久久国产精品大桥未久av | 久久久久视频综合| 伦理电影大哥的女人| 麻豆精品久久久久久蜜桃| 中文精品一卡2卡3卡4更新| 色94色欧美一区二区| 热re99久久精品国产66热6| 午夜久久久在线观看| 色网站视频免费| 欧美成人精品欧美一级黄| 狂野欧美激情性xxxx在线观看| 免费观看在线日韩| 自线自在国产av| 青青草视频在线视频观看|