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

    船用電子機(jī)柜力學(xué)性能仿真計(jì)算

    2024-12-03 00:00:00羅賢成
    指揮控制與仿真 2024年6期

    摘 要:為了指導(dǎo)和優(yōu)化結(jié)構(gòu)設(shè)計(jì),提出了一種快速預(yù)測(cè)船用電子機(jī)柜力學(xué)試驗(yàn)性能的數(shù)值模擬方法。根據(jù)船用電子設(shè)備結(jié)構(gòu)力學(xué)試驗(yàn)標(biāo)準(zhǔn),對(duì)某型船用電子機(jī)柜進(jìn)行力學(xué)試驗(yàn)仿真分析。首先,利用有限元軟件ANSYS建立電子機(jī)柜的有限元模型;其次,基于有限元數(shù)值求解方法,計(jì)算機(jī)柜的模態(tài)頻率、正弦掃頻響應(yīng)、顛震響應(yīng)和沖擊響應(yīng);最后,根據(jù)計(jì)算結(jié)果對(duì)機(jī)柜進(jìn)行強(qiáng)度校核,驗(yàn)證機(jī)柜力學(xué)設(shè)計(jì)可靠性。結(jié)果表明,機(jī)柜殼體頂部和底部直角、側(cè)面凸臺(tái)以及背部直角附近的應(yīng)力較高,但都低于機(jī)柜殼體鑄件的最大許用應(yīng)力;機(jī)柜在各個(gè)試驗(yàn)條件下具有較高的可靠性。船用機(jī)柜的力學(xué)性能仿真計(jì)算方法對(duì)同類型的其他電子設(shè)備力學(xué)分析具有較好的參考價(jià)值。

    關(guān)鍵詞:機(jī)柜;模態(tài)頻率;正弦掃頻;沖擊響應(yīng)

    中圖分類號(hào):U674.7 文獻(xiàn)標(biāo)志碼:A DOI:10.3969/j.issn.1673-3819.2024.06.018

    Simulation calculation on mechanical properties of marine electronic cabinet

    LUO Xiancheng

    (Jiangsu Automation Research Institute, Lianyungang 222061, China)

    Abstract:In order to guide and optimize structural design, this paper proposes a numerical simulation method for quickly predicting the mechanical test performance of marine electronic cabinets. According to the structural mechanics test standard of marine electronic equipment, the mechanical test simulation analysis of a type of marine electronic cabinet is carried out. Firstly, the finite element model of the electronic cabinet is established by ANSYS. Secondly, based on finite numerical solution method, the modal frequency, sine sweep vibration response, bumping response and shock response of the cabinet were obtained. Finally, according to the calculation results, the strength of the cabinet is checked, and the mechanical reliability of the cabinet is analyzed. The results show that the stresses near the right angles at the top and bottom, side boss and the right angles at the back are relatively high, but they are all lower than the maximum allowable stress of the cabinet, and the cabinet can work reliably under various test conditions. The modeling and analysis method presented in this paper can provide a good reference for the mechanical analysis of similar electronic equipment.

    Key words:cabinet; modal frequency; sine sweep; shock response

    收稿日期:2023-12-04 修回日期:2024-01-08

    作者簡(jiǎn)介:

    羅賢成(1993—),男,碩士,工程師,研究方向?yàn)榇半娮釉O(shè)備總體設(shè)計(jì)與結(jié)構(gòu)研發(fā)。

    由于船舶各種電子設(shè)備經(jīng)常受到船舶隨機(jī)振動(dòng)、顛震及沖擊等動(dòng)載荷作用,機(jī)柜等電子設(shè)備的結(jié)構(gòu)設(shè)計(jì)必須保證強(qiáng)度、剛度以及疲勞壽命滿足相應(yīng)要求。船用電子設(shè)備環(huán)境使用條件惡劣,由環(huán)境問題造成設(shè)備失效的比例高達(dá)50%,其中振動(dòng)因素占比達(dá)到27%。設(shè)備產(chǎn)生的共振問題可能會(huì)導(dǎo)致設(shè)備電性能差、零部件功能失效和設(shè)備疲勞損傷甚至破壞等問題[1]。同時(shí),相對(duì)于船舶有一定距離的水下非接觸爆炸作用將引起船舶突發(fā)的瞬態(tài)相對(duì)運(yùn)動(dòng),該運(yùn)動(dòng)對(duì)船舶而言就是沖擊載荷。水面船舶不可避免地遭受非接觸爆炸引起的沖擊問題,非接觸式水下爆炸產(chǎn)生氣泡脈動(dòng)和沖擊載荷[2],船舶設(shè)備在沖擊作用下極易失效。為提高船舶整體抗沖擊能力,必須在設(shè)計(jì)階段將船舶設(shè)備抗振動(dòng)和抗沖擊性能作為一項(xiàng)重要的設(shè)計(jì)因素[3],針對(duì)船舶電子設(shè)備剛強(qiáng)度沖擊特性進(jìn)行試驗(yàn)驗(yàn)證分析。

    國(guó)內(nèi)已有很多學(xué)者對(duì)電子設(shè)備抗振動(dòng)和沖擊性能展開了研究。高學(xué)斌等人[4]采用有限元分析軟件對(duì)電子機(jī)柜進(jìn)行了仿真分析,運(yùn)用整體法及質(zhì)量點(diǎn)法分別建立了機(jī)柜有限元模型,獲得了機(jī)柜在x、y、z方向的固有頻率值。朱曾輝等人[5]基于諧響應(yīng)理論分析,對(duì)氣液交換器的離心風(fēng)機(jī)安裝架進(jìn)行正弦振動(dòng)計(jì)算。金恒林[6]基于Workbench研究分析了某機(jī)載電子設(shè)備在惡劣振動(dòng)環(huán)境試驗(yàn)條件下的隨機(jī)振動(dòng)特性。崔健祿等人[7]采用Workbench有限元分析軟件探討電子設(shè)備在后坐力沖擊載荷下的靜力學(xué)和瞬態(tài)結(jié)構(gòu)響應(yīng)。郝云鵬等人[8]通過模擬振動(dòng)臺(tái)試驗(yàn),研究機(jī)柜的破壞情況。李健等人[9]根據(jù)艦載環(huán)境振動(dòng)和沖擊特性研究機(jī)柜減振器的特性。但這些研究仍然缺乏系統(tǒng)的力學(xué)試驗(yàn)仿真研究方法。本文從節(jié)省試驗(yàn)成本和快速預(yù)判結(jié)構(gòu)設(shè)計(jì)合理性角度出發(fā),以某型電子機(jī)柜作為研究對(duì)象,結(jié)合使用ANSYS軟件進(jìn)行了系統(tǒng)、全面的有限元仿真分析,參考相應(yīng)的國(guó)家標(biāo)準(zhǔn)[10-12],模擬機(jī)柜掃頻分析、半正弦波顛震和沖擊試驗(yàn)條件,對(duì)柜體進(jìn)行強(qiáng)度校核分析,并依據(jù)校核標(biāo)準(zhǔn)給出了校核結(jié)果。本文的研究方法對(duì)各種船用設(shè)備力學(xué)性能分析均具有參考價(jià)值。

    1 理論分析

    電子設(shè)備系統(tǒng)結(jié)構(gòu)復(fù)雜,在動(dòng)力學(xué)分析中常將其簡(jiǎn)化為多自由度系統(tǒng),而單自由度系統(tǒng)是多自由度系統(tǒng)中最簡(jiǎn)單的形式,多自由度系統(tǒng)振動(dòng)的許多基本概念將先在單自由度系統(tǒng)的討論中建立,進(jìn)而推廣到其他多自由度系統(tǒng)中[13]。圖1為單自由度彈簧阻尼系統(tǒng)。

    2 力學(xué)仿真分析

    基于以上理論分析,本文利用有限元分析軟件ANSYS對(duì)某典型電子機(jī)柜力學(xué)試驗(yàn)進(jìn)行了仿真計(jì)算分析,可以快速預(yù)測(cè)機(jī)柜在各試驗(yàn)工況下的力學(xué)性能,為機(jī)柜設(shè)計(jì)提供參考。

    2.1 力學(xué)模型

    機(jī)柜由柜體、各單元、減震器、轉(zhuǎn)接插板等組成,外形尺寸為540 mm(寬)×580 mm(深,不含減震器)×1 600 mm(高),機(jī)柜外觀見圖2。柜體具有安裝和保護(hù)各插箱、各模塊作用,是直接承受外部載荷的結(jié)構(gòu)基礎(chǔ),其結(jié)構(gòu)設(shè)計(jì)的好壞,將直接影響整體設(shè)備的工作穩(wěn)定性和可靠性。柜體主體結(jié)構(gòu)采用鋁合金(ZL101A)整體鑄造而成,其底部減震器可與船體連接。柜體內(nèi)的各插箱單元主要用于實(shí)現(xiàn)機(jī)柜顯示、計(jì)算、處理等功能,通過面板螺釘、導(dǎo)軌和定位銷與柜體相連。坐標(biāo)系滿足右手法則,規(guī)定為:柜前由右向左為X正向;由柜后向柜前為Y向正向;豎直向上為Z正向。

    各插箱單元通過面板螺釘、定位銷與柜體進(jìn)行連接。為了簡(jiǎn)化有限元計(jì)算,柜體內(nèi)單元均簡(jiǎn)化成塊,保證各單元的質(zhì)量配平,這樣既能考慮機(jī)柜和各插箱單元的剛度和質(zhì)量,又可以重點(diǎn)分析柜體主體框架的結(jié)構(gòu)設(shè)計(jì),結(jié)合振動(dòng)、顛震、沖擊等力學(xué)試驗(yàn)有限元分析,獲得柜體主體框架的強(qiáng)度校核結(jié)果。

    本文根據(jù)機(jī)柜的工作狀態(tài)和受力情況簡(jiǎn)化幾何模型,由于結(jié)構(gòu)較為復(fù)雜,為了減小計(jì)算規(guī)模,同時(shí)提高計(jì)算效率,將原幾何模型中對(duì)分析結(jié)果影響較小的幾何特征去除,如幾何模型中的倒角、圓角和非安裝連接孔等。本文后續(xù)力學(xué)分析均是基于簡(jiǎn)化模型進(jìn)行的。

    簡(jiǎn)化模型共分為兩個(gè)層次。第一層次為減震器,減震器可以看作帶有阻尼的彈簧,在實(shí)際分析過程中,結(jié)合工程經(jīng)驗(yàn),將減震器看作剛性體,將經(jīng)過減震器衰減后的載荷作為分析輸入。第二層次為空殼體,殼體和減震器、各個(gè)單元面板采用無摩擦接觸(frictionless接觸),就是說,一旦兩個(gè)物體之間出現(xiàn)了分離,則法向力為零。因此,當(dāng)外力發(fā)生改變時(shí),接觸面之間可能會(huì)分開,也可能會(huì)閉合,閉合情況下假設(shè)摩擦系數(shù)為零,即發(fā)生切向相對(duì)滑動(dòng)時(shí),沒有摩擦力,由于不考慮摩擦,計(jì)算所得的應(yīng)力值大于實(shí)際值,這實(shí)際是一種安全系數(shù)較高的算法。各個(gè)組件和螺釘頭及螺紋接觸的部分采用綁定接觸(bonded),使用綁定接觸以后,在接觸面或者接觸邊之間不存在切向的相對(duì)滑動(dòng)或者法向的相對(duì)分離,該種設(shè)置符合工程實(shí)際。

    機(jī)柜內(nèi)各個(gè)部件材料見表1。

    簡(jiǎn)化后的殼體網(wǎng)格劃分如圖3所示,網(wǎng)格總數(shù)為785 770,節(jié)點(diǎn)總數(shù)為1 305 467。

    通過安裝位置可知,機(jī)柜背部和底部均安裝減震器,共6個(gè)安裝孔,如圖4所示。

    2.2 模態(tài)分析

    機(jī)柜殼體的前六階固有頻率計(jì)算結(jié)果如表2所示。圖5為機(jī)柜殼體前六階各模態(tài)振型,第一階、第二階為殼體前后方向振動(dòng),第三階、第四階為殼體左右方向振動(dòng),第五階、第六階為殼體垂直方向振動(dòng)。掃頻試驗(yàn)及耐振試驗(yàn)可重點(diǎn)關(guān)注殼體整體共振點(diǎn)的應(yīng)力及變形情況。

    2.3 正弦掃頻分析

    振動(dòng)試驗(yàn)激勵(lì)條件見表3。對(duì)機(jī)柜殼體的有限元模型施加正弦振動(dòng)載荷譜作用,選取結(jié)構(gòu)的阻尼系數(shù)為0.05,由于載荷譜激勵(lì)的頻率范圍為1 Hz~160 Hz,施加載荷以保證覆蓋載荷譜激勵(lì)的頻率范圍。圖6為載荷輸入曲線。

    分別對(duì)機(jī)柜三個(gè)方向進(jìn)行正弦振動(dòng)掃頻分析,圖7~圖9為殼體在三個(gè)方向正弦振動(dòng)激勵(lì)作用下的加速度頻率響應(yīng)函數(shù),可以看出,當(dāng)外部激勵(lì)為64 Hz和128 Hz時(shí)曲線出現(xiàn)了峰值,此時(shí)外部激勵(lì)與殼體的低階和高階固有頻率均產(chǎn)生共振現(xiàn)象,殼體頻率響應(yīng)函數(shù)激勵(lì)曲線出現(xiàn)明顯激增,殼體將出現(xiàn)較強(qiáng)的應(yīng)力集中現(xiàn)象。

    圖10~圖12為殼體在三個(gè)方向正弦激勵(lì)作用下,激勵(lì)分別為64 Hz和128 Hz作用時(shí)殼體等效應(yīng)力云圖,統(tǒng)計(jì)結(jié)果如表4所示??梢钥闯觯?dāng)正弦激勵(lì)為y向64 Hz時(shí),殼體的等效應(yīng)力達(dá)到最大值,為95.524 MPa,但最大應(yīng)力仍在材料的許用范圍內(nèi),因此,殼體在承受三個(gè)方向正弦振動(dòng)載荷作用下的結(jié)構(gòu)設(shè)計(jì)是安全的。

    耐振試驗(yàn)是在指定頻率段內(nèi)若無危險(xiǎn)頻率段,即固有頻率不在頻率1 Hz~160 Hz區(qū)間,對(duì)機(jī)柜施加160 Hz的加速度載荷,加速度載荷取值25 m/s2;若有危險(xiǎn)頻率段,則在危險(xiǎn)頻率上施加加速度載荷。由此可見,耐振試驗(yàn)中的振動(dòng)應(yīng)力等同于正弦掃頻的最大應(yīng)力,殼體最大振動(dòng)應(yīng)力為95.524 MPa,滿足材料的強(qiáng)度要求。

    2.4 顛震分析

    對(duì)機(jī)柜沿豎直方向受到顛震試驗(yàn)下的結(jié)構(gòu)響應(yīng)進(jìn)行研究。激勵(lì)條件見表5,顛震波形為半正弦波。

    沖擊激勵(lì)函數(shù)為

    ag=70 000sin(196.35t),當(dāng)t≤0.016 s時(shí)

    0,當(dāng)t>0.016 s時(shí)單位:mm/s2。

    根據(jù)工程經(jīng)驗(yàn),激勵(lì)曲線和機(jī)柜底部響應(yīng)曲線見圖13(紅色實(shí)線為激勵(lì)曲線,藍(lán)色虛線為響應(yīng)曲線)。

    系統(tǒng)激勵(lì)的加速度值約為7g,響應(yīng)為5.8g,第一個(gè)響應(yīng)波時(shí)長(zhǎng)為71 ms,響應(yīng)總長(zhǎng)約為440 ms。外部激勵(lì)為7 g/16 ms,30次/min,脈沖間隔為2 s。在ANSYS分析中不考慮機(jī)柜減震器結(jié)構(gòu),所以直接施加5.8g沖擊載荷,其中脈沖時(shí)長(zhǎng)為70 ms。圖14為顛震載荷輸入曲線。

    在ANSYS中通過模態(tài)疊加法的瞬態(tài)分析進(jìn)行顛震仿真,模態(tài)疊加法又稱“振型疊加法”,它是以系統(tǒng)無阻尼的振型(模態(tài))為空間基底,通過坐標(biāo)變換,使原動(dòng)力方程解耦,求解n個(gè)相互獨(dú)立的方程獲得模態(tài)位移,進(jìn)而通過疊加各階模態(tài)的貢獻(xiàn)求得系統(tǒng)的響應(yīng)。

    表6給出了殼體在顛震條件下的應(yīng)力響應(yīng)結(jié)果。數(shù)據(jù)表明,殼體在承受豎直加速度載荷作用下的結(jié)構(gòu)設(shè)計(jì)是安全的。圖15給出了機(jī)柜殼體在豎直方向加速度顛震載荷作用下的應(yīng)力響應(yīng)云圖。

    2.5 沖擊分析

    根據(jù)GJB150.16-1986,機(jī)柜采用垂向打錘試驗(yàn),打錘輸入加速度量值為100g的半正弦波(半周期1 ms),機(jī)柜底部輸出響應(yīng)為30g(半周期2 ms)的逐漸衰減的類正弦波。

    機(jī)柜通過6個(gè)接口安裝在減震器上,在跟隨船舶進(jìn)行各種航行姿態(tài)時(shí)會(huì)受到不同方向的沖擊作用。沖擊工況也采用模態(tài)疊加法的瞬態(tài)分析進(jìn)行分析,以機(jī)柜底部實(shí)測(cè)響應(yīng)作為輸入,在ANSYS中輸入隨時(shí)間變化的沖擊載荷,圖16為垂向沖擊輸入。

    表7給出了殼體在沖擊試驗(yàn)下的應(yīng)力響應(yīng)結(jié)果。數(shù)據(jù)表明,殼體最大應(yīng)力發(fā)生在承受垂向沖擊時(shí),最大應(yīng)力為146.76 MPa,如圖17所示,殼體最大應(yīng)力小于許用應(yīng)力,各零部件在承受垂向沖擊是安全的。

    3 結(jié)束語

    本文基于有限元分析軟件ANSYS建立了船用電子機(jī)柜的有限元模型,根據(jù)相應(yīng)的試驗(yàn)標(biāo)準(zhǔn),確定了分析的邊界條件和載荷輸入條件,進(jìn)行了模態(tài)分析、正弦掃頻分析、顛震分析及沖擊響應(yīng)分析。其中,正弦掃頻分析并未考慮減震器的衰減作用,實(shí)際的安全系數(shù)應(yīng)該比文中的更高。顛震分析及沖擊響應(yīng)分析是以同類機(jī)柜實(shí)測(cè)的底部響應(yīng)作為輸入條件進(jìn)行分析的,該輸入條件考慮了減震器的衰減作用。計(jì)算結(jié)果表明,電子機(jī)柜在各種試驗(yàn)條件下工作具備較高的可靠性。本文的研究方法對(duì)其他船用設(shè)備力學(xué)分析具有較好的指導(dǎo)意義。

    參考文獻(xiàn):

    [1] 戈進(jìn)飛. 軍用電子設(shè)備結(jié)構(gòu)安全性設(shè)計(jì)概述[J]. 電子機(jī)械工程, 2010, 26(5): 1-6, 27.

    GE J F. Summary of military electronic equipment structure security design[J]. Electro-Mechanical Engineering, 2010, 26(5): 1-6, 27.

    [2] 董九亭, 劉建湖, 汪俊, 等. 水下爆炸下艦艇不同部位沖擊環(huán)境數(shù)值分析[J]. 中國(guó)艦船研究, 2018, 13(5): 32-38.

    DONG J T, LIU J H, WANG J, et al. Numerical analysis on shock environment in different ship regions subjected to underwater explosion[J]. Chinese Journal of Ship Research, 2018, 13(5): 32-38.

    [3] 羅思思.船舶電子設(shè)備的抗強(qiáng)沖擊試驗(yàn)分析[J]. 船舶科學(xué)技術(shù),2018,40(7):115-117.

    LUO S S. Test and analysis of anti strong impact of ship electronic equipment[J]. Ship Science and Technology, 2018,40(7):115-117.

    [4] 高學(xué)斌, 代廷振. 電子設(shè)備機(jī)柜模態(tài)分析研究及試驗(yàn)驗(yàn)證[J]. 電工技術(shù), 2020(23): 52-55, 59.

    GAO X B, DAI T Z. Modal analysis and test verification of electronic equipment cabinet[J]. Electric Engineering, 2020(23): 52-55, 59.

    [5] 朱曾輝, 李維忠. 艦船振動(dòng)環(huán)境下機(jī)柜氣液熱交換器振動(dòng)故障分析[J]. 雷達(dá)與對(duì)抗, 2023, 43(3): 57-61.

    ZHU Z H, LI W Z. Vibration fault analysis of gas-liquid heat exchanger in cabinet under ship vibration environment[J]. Radar & ECM, 2023, 43(3): 57-61.

    [6] 金恒林. 某機(jī)載電子設(shè)備的抗振能力分析[J]. 環(huán)境技術(shù), 2022, 40(4): 182-186, 191.

    JIN H L. Analysis of vibration resistance of airborne electronic equipment[J]. Environmental Technology, 2022, 40(4): 182-186, 191.

    [7] 崔健祿, 白鵬飛. 沖擊載荷下電子設(shè)備結(jié)構(gòu)強(qiáng)度響應(yīng)分析[J]. 光電技術(shù)應(yīng)用, 2022, 37(3): 74-80.

    CUI J L, BAI P F. Strength response analysis of electronic equipment structure under impact load[J]. Electro-Optic Technology Application, 2022, 37(3): 74-80.

    [8] 郝云鵬, 馮利飛, 毛晨曦. 典型通信設(shè)備機(jī)柜的振動(dòng)臺(tái)試驗(yàn)研究[J]. 地震研究, 2020, 43(3): 485-490.

    HAO Y P, FENG L F, MAO C X. Study on shaking table tests of a typical telecommunication cabinet[J]. Journal of Seismological Research, 2020, 43(3): 485-490.

    [9] 李健, 杭佳棋, 張恒, 等. 某型艦載電子機(jī)柜的減振器設(shè)計(jì)[J]. 噪聲與振動(dòng)控制, 2018, 38(1): 151-155.

    LI J, HANG J Q, ZHANG H, et al. Design of shock absorbers for a shipboard electronic cabinet[J]. Noise and Vibration Control, 2018, 38(1): 151-155.

    [10]海軍裝備論證研究中心標(biāo)準(zhǔn)規(guī)范研究室. GJB 150.16—86 軍用設(shè)備環(huán)境試驗(yàn)方法 振動(dòng)試驗(yàn)[S] . 北京:國(guó)防科學(xué)技術(shù)工業(yè)委員會(huì),1986.

    Naval Equipment Demonstration Research Center standard specification research office. GJB 150. 18-86 Environmental test methods for military equipment vibration test[S]. Beijing: Commission of Science, Technology and Industry for National Defense Standards Initial Issue Department, 1986.

    [11]海軍裝備論證研究中心標(biāo)準(zhǔn)規(guī)范研究室. GJB 150.18—86 軍用設(shè)備環(huán)境試驗(yàn)方法 沖擊試驗(yàn)[S] . 北京:國(guó)防科學(xué)技術(shù)工業(yè)委員會(huì),1986.

    Naval Equipment Demonstration Research Center standard specification research office. GJB 150. 18-86 Environmental test methods for military equipment impact test[S]. Beijing: Commission of Science, Technology and Industry for National Defense, 1986.

    [12]海軍裝備論證研究中心標(biāo)準(zhǔn)規(guī)范研究室. GJB 4.8-83 船舶電子設(shè)備環(huán)境試驗(yàn)顛震試驗(yàn)[S]. 北京:國(guó)防科學(xué)技術(shù)工業(yè)委員會(huì),1983.

    Naval Equipment Demonstration Research Center standard specification research office. GJB 4. 8-83 Environmental test for marine electronic equipment shock test[S]. Beijing: Commission of Science, Technology and Industry for National Defense, 1983.

    [13]邱成悌. 電子設(shè)備結(jié)構(gòu)設(shè)計(jì)原理[M]. 南京: 東南大學(xué)出版社, 2002.

    QIU C T. Structural design principle of electronic equipment[M]. Nanjing: Southeast University Press, 2002.

    [14]楊飛兵. 電子設(shè)備動(dòng)力學(xué)分析及其振動(dòng)控制[D]. 西安: 西安電子科技大學(xué),2012.

    YANG F B. Dynamic analysis of electronic equipment and vibration control[D]. Xi’an: Xidian University, 2012.

    (責(zé)任編輯:李楠)

    国产精品久久久久久久电影| 天堂中文最新版在线下载| 精品久久久噜噜| 欧美老熟妇乱子伦牲交| 91午夜精品亚洲一区二区三区| 亚洲欧美成人综合另类久久久| 亚洲欧洲日产国产| 亚洲国产色片| av在线老鸭窝| 亚洲av二区三区四区| 男人舔奶头视频| 亚洲自偷自拍三级| 婷婷色综合www| 国产欧美日韩一区二区三区在线 | 国产在线男女| 在线观看av片永久免费下载| 欧美精品国产亚洲| 午夜老司机福利剧场| 香蕉精品网在线| 欧美高清成人免费视频www| 亚洲第一区二区三区不卡| 日韩不卡一区二区三区视频在线| 啦啦啦视频在线资源免费观看| 99热6这里只有精品| 大又大粗又爽又黄少妇毛片口| 欧美激情国产日韩精品一区| 哪个播放器可以免费观看大片| 欧美少妇被猛烈插入视频| 久久国内精品自在自线图片| 高清在线视频一区二区三区| 99精国产麻豆久久婷婷| 精品国产一区二区久久| 嫩草影院入口| 国产伦精品一区二区三区四那| 人体艺术视频欧美日本| 亚洲国产av新网站| 日韩欧美精品免费久久| 国产亚洲av片在线观看秒播厂| 五月开心婷婷网| 国产精品无大码| 国产视频首页在线观看| 国产精品久久久久久精品电影小说| 久久久国产欧美日韩av| 日本-黄色视频高清免费观看| 国产精品秋霞免费鲁丝片| 男女边吃奶边做爰视频| 少妇的逼好多水| 七月丁香在线播放| 七月丁香在线播放| 中国三级夫妇交换| 99久久中文字幕三级久久日本| 自拍欧美九色日韩亚洲蝌蚪91 | 人妻少妇偷人精品九色| 久久精品久久久久久久性| 亚洲图色成人| 欧美变态另类bdsm刘玥| 边亲边吃奶的免费视频| 午夜激情久久久久久久| 久久精品夜色国产| 国产精品99久久久久久久久| 精品久久久久久久久亚洲| 在线观看美女被高潮喷水网站| 午夜激情福利司机影院| 精品少妇黑人巨大在线播放| 国产精品秋霞免费鲁丝片| 卡戴珊不雅视频在线播放| 特大巨黑吊av在线直播| 大陆偷拍与自拍| 国产精品秋霞免费鲁丝片| 午夜福利视频精品| 久久久a久久爽久久v久久| 精品国产一区二区三区久久久樱花| 菩萨蛮人人尽说江南好唐韦庄| 狂野欧美白嫩少妇大欣赏| 黄色欧美视频在线观看| 亚洲精品亚洲一区二区| 视频中文字幕在线观看| 久久久久久久国产电影| 一级毛片黄色毛片免费观看视频| 午夜激情久久久久久久| 久久精品国产亚洲网站| 黄色配什么色好看| 99热国产这里只有精品6| 免费观看性生交大片5| 另类精品久久| .国产精品久久| 亚洲美女视频黄频| av国产久精品久网站免费入址| 九九爱精品视频在线观看| 久久久久国产网址| 一级毛片 在线播放| 一区二区三区乱码不卡18| 久久精品国产鲁丝片午夜精品| 久久久久久久久久久免费av| 免费不卡的大黄色大毛片视频在线观看| 一区二区三区乱码不卡18| 丁香六月天网| 波野结衣二区三区在线| 777米奇影视久久| 老司机亚洲免费影院| 久久狼人影院| 在线观看免费日韩欧美大片 | 下体分泌物呈黄色| 人妻一区二区av| 街头女战士在线观看网站| 久久综合国产亚洲精品| 亚洲美女搞黄在线观看| 欧美xxⅹ黑人| 久久午夜综合久久蜜桃| 丝瓜视频免费看黄片| h日本视频在线播放| 精品国产一区二区久久| 久久久久精品性色| 熟女av电影| 日本av免费视频播放| 久久久久久久精品精品| 久久久国产欧美日韩av| 最新的欧美精品一区二区| 国产精品一区二区在线不卡| 国产91av在线免费观看| 国产熟女欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 午夜av观看不卡| 偷拍熟女少妇极品色| 91精品国产九色| 欧美激情极品国产一区二区三区 | 久久久亚洲精品成人影院| 男女国产视频网站| 99热这里只有精品一区| 欧美日韩综合久久久久久| 爱豆传媒免费全集在线观看| 老司机影院成人| av专区在线播放| 日韩av在线免费看完整版不卡| 国国产精品蜜臀av免费| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 亚洲国产精品专区欧美| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 免费观看a级毛片全部| 久久 成人 亚洲| 老熟女久久久| 亚洲第一区二区三区不卡| 日本猛色少妇xxxxx猛交久久| 久久久久国产精品人妻一区二区| 国产精品99久久久久久久久| 国产成人免费观看mmmm| 99热6这里只有精品| 久久av网站| av免费观看日本| 伦理电影免费视频| 欧美日韩av久久| 尾随美女入室| 国产成人免费观看mmmm| 国产日韩欧美视频二区| 女性生殖器流出的白浆| 五月玫瑰六月丁香| 国产无遮挡羞羞视频在线观看| 国产在线免费精品| 亚洲一级一片aⅴ在线观看| 乱系列少妇在线播放| 久久精品久久久久久噜噜老黄| 男男h啪啪无遮挡| 免费黄频网站在线观看国产| 只有这里有精品99| 性色avwww在线观看| 中文在线观看免费www的网站| 五月开心婷婷网| 国产精品成人在线| 97超视频在线观看视频| 日韩欧美精品免费久久| 99热这里只有是精品50| 国产女主播在线喷水免费视频网站| √禁漫天堂资源中文www| 99热6这里只有精品| 亚洲久久久国产精品| 观看av在线不卡| av卡一久久| 精品亚洲成国产av| 国产免费福利视频在线观看| 亚洲欧洲日产国产| 国产视频内射| 在线播放无遮挡| 国产成人精品无人区| 亚洲一区二区三区欧美精品| 看非洲黑人一级黄片| 亚洲美女视频黄频| 在线播放无遮挡| 伊人久久精品亚洲午夜| www.av在线官网国产| 久久99精品国语久久久| 欧美 日韩 精品 国产| 亚洲精品456在线播放app| 久久精品国产a三级三级三级| 午夜影院在线不卡| 人妻制服诱惑在线中文字幕| 纵有疾风起免费观看全集完整版| 如日韩欧美国产精品一区二区三区 | 人人妻人人澡人人看| 国产亚洲最大av| 精品国产国语对白av| 欧美丝袜亚洲另类| 国产欧美日韩精品一区二区| 欧美日韩综合久久久久久| 少妇被粗大猛烈的视频| 青春草亚洲视频在线观看| 成年人午夜在线观看视频| 久久久国产欧美日韩av| 波野结衣二区三区在线| 亚洲av电影在线观看一区二区三区| 久久精品国产a三级三级三级| 国产永久视频网站| 热re99久久国产66热| 成年女人在线观看亚洲视频| 99热网站在线观看| a级毛色黄片| 亚洲成人av在线免费| 激情五月婷婷亚洲| 亚洲不卡免费看| 天天躁夜夜躁狠狠久久av| 亚洲欧美日韩卡通动漫| 另类亚洲欧美激情| 国产一区亚洲一区在线观看| 极品教师在线视频| 亚洲激情五月婷婷啪啪| 国产亚洲精品久久久com| 国产av一区二区精品久久| 久久久久精品久久久久真实原创| 狂野欧美白嫩少妇大欣赏| 国产黄片美女视频| 偷拍熟女少妇极品色| 国产高清三级在线| 五月伊人婷婷丁香| 秋霞在线观看毛片| 大话2 男鬼变身卡| 美女内射精品一级片tv| 人人妻人人看人人澡| 国产一区二区在线观看av| 亚洲精品日本国产第一区| 中文字幕久久专区| 国产探花极品一区二区| 精品亚洲乱码少妇综合久久| 国产色婷婷99| 人人妻人人澡人人爽人人夜夜| 99久久精品热视频| 国产亚洲欧美精品永久| 免费观看的影片在线观看| 精品人妻一区二区三区麻豆| 人妻制服诱惑在线中文字幕| 久久99蜜桃精品久久| 少妇人妻久久综合中文| 人妻一区二区av| .国产精品久久| 久久久午夜欧美精品| 欧美精品国产亚洲| 嫩草影院入口| 国产无遮挡羞羞视频在线观看| 赤兔流量卡办理| 高清黄色对白视频在线免费看 | 久久国产亚洲av麻豆专区| 成人二区视频| 91精品一卡2卡3卡4卡| 日韩人妻高清精品专区| 精品久久国产蜜桃| 丰满少妇做爰视频| 亚洲美女搞黄在线观看| 波野结衣二区三区在线| 少妇人妻久久综合中文| 大香蕉97超碰在线| 成人美女网站在线观看视频| 国产黄色免费在线视频| 亚洲婷婷狠狠爱综合网| 久久久久久久亚洲中文字幕| 久久人妻熟女aⅴ| av国产精品久久久久影院| av专区在线播放| 久久久久久伊人网av| 成人二区视频| 在线观看av片永久免费下载| 国产老妇伦熟女老妇高清| 亚洲中文av在线| 欧美精品国产亚洲| 少妇精品久久久久久久| 日本wwww免费看| 另类精品久久| 在线观看免费日韩欧美大片 | 国产精品麻豆人妻色哟哟久久| 哪个播放器可以免费观看大片| 51国产日韩欧美| 欧美区成人在线视频| 校园人妻丝袜中文字幕| 亚洲国产精品999| 美女内射精品一级片tv| 午夜视频国产福利| 一边亲一边摸免费视频| 精品一区二区三卡| 日韩中文字幕视频在线看片| 99久久人妻综合| 亚洲天堂av无毛| 一级黄片播放器| 亚洲av日韩在线播放| 少妇猛男粗大的猛烈进出视频| 国产精品不卡视频一区二区| av一本久久久久| 丝袜喷水一区| 一区二区av电影网| 精品99又大又爽又粗少妇毛片| 欧美亚洲 丝袜 人妻 在线| 日韩亚洲欧美综合| 免费观看在线日韩| 能在线免费看毛片的网站| 黄色配什么色好看| 嫩草影院新地址| 亚洲欧美成人精品一区二区| 男女边吃奶边做爰视频| 秋霞伦理黄片| 国产国拍精品亚洲av在线观看| 国产乱人偷精品视频| 如日韩欧美国产精品一区二区三区 | 国产精品免费大片| 精品一区二区三卡| 大片电影免费在线观看免费| 免费少妇av软件| 亚洲av.av天堂| 中文字幕亚洲精品专区| 日韩中字成人| 欧美日韩亚洲高清精品| av免费观看日本| 国产精品一二三区在线看| 精品人妻偷拍中文字幕| 乱系列少妇在线播放| 寂寞人妻少妇视频99o| 久久午夜福利片| 九草在线视频观看| 精品少妇久久久久久888优播| 91精品一卡2卡3卡4卡| 国产精品欧美亚洲77777| 日韩视频在线欧美| 一级片'在线观看视频| 日本91视频免费播放| 亚洲av成人精品一二三区| 人妻少妇偷人精品九色| www.av在线官网国产| 久久人人爽人人爽人人片va| 欧美丝袜亚洲另类| 人妻夜夜爽99麻豆av| av在线播放精品| 丁香六月天网| 婷婷色麻豆天堂久久| 一区二区三区乱码不卡18| 欧美人与善性xxx| 亚洲av男天堂| 日本欧美国产在线视频| 18+在线观看网站| 日韩中字成人| 女性被躁到高潮视频| 国产精品99久久久久久久久| www.色视频.com| 欧美精品亚洲一区二区| 精品酒店卫生间| 午夜久久久在线观看| 欧美日韩亚洲高清精品| 在线看a的网站| 综合色丁香网| 丝袜喷水一区| 如日韩欧美国产精品一区二区三区 | 黑丝袜美女国产一区| 久久ye,这里只有精品| av在线app专区| 免费人妻精品一区二区三区视频| 大香蕉97超碰在线| 99re6热这里在线精品视频| 免费不卡的大黄色大毛片视频在线观看| 熟女人妻精品中文字幕| 欧美日韩综合久久久久久| 一级毛片久久久久久久久女| 日本欧美视频一区| 69精品国产乱码久久久| 在线 av 中文字幕| 成人免费观看视频高清| 搡老乐熟女国产| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 美女福利国产在线| 两个人免费观看高清视频 | 在线观看www视频免费| 一级毛片黄色毛片免费观看视频| 国产精品女同一区二区软件| 久久久国产欧美日韩av| 久久国产乱子免费精品| 黄色视频在线播放观看不卡| h日本视频在线播放| 我要看日韩黄色一级片| 国产日韩欧美亚洲二区| 亚洲欧美精品专区久久| 夜夜骑夜夜射夜夜干| 亚洲在久久综合| 一级a做视频免费观看| 中文字幕人妻熟人妻熟丝袜美| 国产午夜精品久久久久久一区二区三区| 精品人妻熟女av久视频| 五月开心婷婷网| 久久人人爽人人爽人人片va| 精品99又大又爽又粗少妇毛片| 中文天堂在线官网| 欧美性感艳星| 国产综合精华液| 久久 成人 亚洲| 看免费成人av毛片| 精品国产一区二区久久| 久久精品久久久久久久性| 97超碰精品成人国产| 久久久久视频综合| 久久精品国产亚洲网站| 国产在线免费精品| 免费不卡的大黄色大毛片视频在线观看| 韩国高清视频一区二区三区| 天天躁夜夜躁狠狠久久av| 国产日韩一区二区三区精品不卡 | 十分钟在线观看高清视频www | 国产一区二区三区综合在线观看 | 97超碰精品成人国产| 亚洲一级一片aⅴ在线观看| 一级av片app| 91午夜精品亚洲一区二区三区| 中文字幕精品免费在线观看视频 | 2021少妇久久久久久久久久久| 久久久久国产网址| 亚洲欧洲精品一区二区精品久久久 | 国产有黄有色有爽视频| 日本爱情动作片www.在线观看| 大又大粗又爽又黄少妇毛片口| av免费观看日本| 国产精品一区www在线观看| 五月开心婷婷网| 久久精品久久精品一区二区三区| 麻豆成人av视频| 欧美少妇被猛烈插入视频| 伊人久久国产一区二区| 男人舔奶头视频| 国产精品熟女久久久久浪| 黄色配什么色好看| 成人影院久久| 精品久久久精品久久久| 日本黄色日本黄色录像| 久久ye,这里只有精品| 国产免费又黄又爽又色| 日韩制服骚丝袜av| 特大巨黑吊av在线直播| 国产精品福利在线免费观看| 国产精品久久久久久精品电影小说| 黑丝袜美女国产一区| 亚洲,一卡二卡三卡| 亚洲国产精品999| 男女边吃奶边做爰视频| 大又大粗又爽又黄少妇毛片口| 精品国产一区二区久久| 狂野欧美激情性bbbbbb| 欧美日韩综合久久久久久| 日日摸夜夜添夜夜爱| 性高湖久久久久久久久免费观看| 免费观看的影片在线观看| 国产永久视频网站| 婷婷色av中文字幕| 久久久久久久久大av| 中文字幕久久专区| 青青草视频在线视频观看| 看十八女毛片水多多多| 在线 av 中文字幕| 欧美精品一区二区免费开放| 久久精品国产a三级三级三级| 亚洲综合色惰| 十八禁高潮呻吟视频 | av视频免费观看在线观看| 亚洲精品成人av观看孕妇| 大话2 男鬼变身卡| 建设人人有责人人尽责人人享有的| 少妇人妻一区二区三区视频| 国产日韩欧美亚洲二区| 欧美另类一区| 九草在线视频观看| 国产精品熟女久久久久浪| 2022亚洲国产成人精品| 亚洲不卡免费看| 中文精品一卡2卡3卡4更新| 国产成人午夜福利电影在线观看| 黄色毛片三级朝国网站 | 少妇熟女欧美另类| 人人妻人人澡人人看| 欧美一级a爱片免费观看看| 我要看日韩黄色一级片| 黄色日韩在线| 国产黄色视频一区二区在线观看| 亚洲熟女精品中文字幕| 国产高清三级在线| 99九九线精品视频在线观看视频| av天堂中文字幕网| 狂野欧美激情性bbbbbb| 精品久久久久久久久亚洲| 少妇精品久久久久久久| 中国国产av一级| 亚洲av二区三区四区| 亚洲国产色片| 精品久久久久久久久亚洲| 亚洲国产毛片av蜜桃av| 美女xxoo啪啪120秒动态图| 色94色欧美一区二区| 曰老女人黄片| 精品久久久久久久久亚洲| 午夜影院在线不卡| 欧美日韩视频高清一区二区三区二| 久久久久久久久久人人人人人人| 最近2019中文字幕mv第一页| 日韩精品有码人妻一区| 精品一区二区三卡| 久久99热这里只频精品6学生| 在线天堂最新版资源| 麻豆乱淫一区二区| 亚洲伊人久久精品综合| 国产在线男女| av卡一久久| 国内少妇人妻偷人精品xxx网站| 三级国产精品片| 少妇人妻久久综合中文| 蜜桃在线观看..| 好男人视频免费观看在线| 日韩人妻高清精品专区| 内地一区二区视频在线| 多毛熟女@视频| 亚洲中文av在线| 国产精品久久久久久精品电影小说| 免费观看a级毛片全部| 久久久久久久亚洲中文字幕| 观看美女的网站| 熟女电影av网| 久久精品久久久久久久性| 99久久中文字幕三级久久日本| 91精品一卡2卡3卡4卡| 欧美日韩国产mv在线观看视频| 亚洲精品视频女| 99久久人妻综合| 激情五月婷婷亚洲| av天堂中文字幕网| 人妻系列 视频| 一级毛片 在线播放| 欧美老熟妇乱子伦牲交| 日韩免费高清中文字幕av| 免费人妻精品一区二区三区视频| 日本91视频免费播放| 美女xxoo啪啪120秒动态图| 全区人妻精品视频| 免费看av在线观看网站| 777米奇影视久久| 国产在线视频一区二区| av女优亚洲男人天堂| 精品国产乱码久久久久久小说| 在线观看免费视频网站a站| 少妇猛男粗大的猛烈进出视频| 精品人妻偷拍中文字幕| 乱人伦中国视频| 久久久久人妻精品一区果冻| 色网站视频免费| 国产视频内射| 日韩中文字幕视频在线看片| 亚洲欧美日韩东京热| 下体分泌物呈黄色| 精品国产一区二区久久| 国产精品人妻久久久久久| 性色avwww在线观看| av播播在线观看一区| 成人美女网站在线观看视频| 国产高清有码在线观看视频| 成人毛片60女人毛片免费| 如何舔出高潮| 九色成人免费人妻av| 国产白丝娇喘喷水9色精品| 水蜜桃什么品种好| 国产成人91sexporn| 人妻一区二区av| 老司机影院成人| 日韩精品有码人妻一区| 人人妻人人看人人澡| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| 久久狼人影院| 亚洲婷婷狠狠爱综合网| 亚洲真实伦在线观看| 国产成人精品一,二区| 日日爽夜夜爽网站| 亚洲精品亚洲一区二区| 在线观看人妻少妇| 久久99热这里只频精品6学生| 伊人亚洲综合成人网| 久久精品国产a三级三级三级| 一级二级三级毛片免费看| 大话2 男鬼变身卡| 亚洲不卡免费看| 美女xxoo啪啪120秒动态图| 三级国产精品欧美在线观看| 中文字幕制服av| 亚洲欧洲精品一区二区精品久久久 | 久久婷婷青草| videossex国产| 国产真实伦视频高清在线观看| 日韩电影二区| 成人无遮挡网站| 如日韩欧美国产精品一区二区三区 | 免费看日本二区| 男女边吃奶边做爰视频| 九九在线视频观看精品| 极品教师在线视频| 伦理电影免费视频| 亚洲欧美精品自产自拍| 高清毛片免费看|