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

    船舶肘板拓?fù)鋬?yōu)化設(shè)計

    2015-02-07 02:55:28程遠(yuǎn)勝劉甜甜劉均
    中國艦船研究 2015年5期
    關(guān)鍵詞:肘板板結(jié)構(gòu)甲板

    程遠(yuǎn)勝,劉甜甜,劉均

    華中科技大學(xué)船舶與海洋工程學(xué)院,湖北武漢430074

    0 引 言

    組成船體縱向或橫向框架的構(gòu)件常采用三角形肘板連接。現(xiàn)代船舶的建造運營表明,肘板的損傷與裂紋的出現(xiàn)大量存在著,其中肘板損傷數(shù)約占整個結(jié)構(gòu)損傷數(shù)的26.85%[1]。這些損傷主要由肘板處的應(yīng)力集中引起,因此,提出新的肘板結(jié)構(gòu)型式,對有效降低節(jié)點結(jié)構(gòu)的應(yīng)力集中程度具有較大的實用價值。

    Kim 等[2]對船舶肘板的材料進(jìn)行了研究和探索,設(shè)計出了一種高強度的低碳素鑄鋼,其屈服極限和極限拉伸強度分別為480 和600 MPa,在保證應(yīng)力集中不惡化的情況下,肘板尺寸和重量分別下降了30%和50%。郭信川等[3]提出了基于修改的粒子群算法的船舶肘板優(yōu)化方法,對肘板節(jié)點結(jié)構(gòu)進(jìn)行了形狀優(yōu)化,最終優(yōu)化后的新型肘板節(jié)點結(jié)構(gòu)較原始肘板節(jié)點結(jié)構(gòu)的應(yīng)力集中程度降低了10.03%。Lim 等[4]研究了肘板對框架接頭處應(yīng)力分布和極限強度的影響,其根據(jù)節(jié)點的力學(xué)特性,重新設(shè)計了肘板的形狀,有效改善了其結(jié)構(gòu)應(yīng)力狀態(tài)。姜以威等[5]研究了梁肘板的應(yīng)力集中,認(rèn)為它是由梁肘板節(jié)點的內(nèi)部特性和外部輸入特性所引起,并分析了梁肘板的焊接工藝性、節(jié)點形式和尺寸對應(yīng)力集中的影響。王波和楊平[6]針對船體的梁連接節(jié)點,在ANSYS 中建立了梁連接節(jié)點的殼單元模型,探討了肘板尺寸變化對節(jié)點承載能力的影響規(guī)律,并對比了幾種常見節(jié)點的強度和屈曲性能。史戰(zhàn)新[7]也針對肘板的優(yōu)化做了一定的工作,其針對水下結(jié)構(gòu)物的肘板結(jié)構(gòu),基于子模型法分析了節(jié)點的結(jié)構(gòu)應(yīng)力,并在此基礎(chǔ)上應(yīng)用ANSYS 與Matlab 的聯(lián)合雙目標(biāo)遺傳算法對其進(jìn)行了形狀優(yōu)化,優(yōu)化后的肘板結(jié)構(gòu)疲勞壽命得到了很大的改善。田旭軍等[8]通過ANSYS軟件和改進(jìn)的遺傳算法,對水下結(jié)構(gòu)物的肘板結(jié)構(gòu)進(jìn)行形狀優(yōu)化設(shè)計,降低了肘板節(jié)點的應(yīng)力集中。

    本文將提出船舶典型節(jié)點肘板的拓?fù)鋬?yōu)化設(shè)計數(shù)學(xué)模型,基于Optistruct 軟件進(jìn)行肘板拓?fù)鋬?yōu)化設(shè)計,并將獲得的肘板拓?fù)鋬?yōu)化設(shè)計結(jié)果進(jìn)行適當(dāng)?shù)墓こ袒幚?,以有效降低?jié)點結(jié)構(gòu)的應(yīng)力集中程度。

    1 船舶肘板節(jié)點結(jié)構(gòu)應(yīng)力分析

    1.1 整體模型強度計算

    選取某船舶三艙段結(jié)構(gòu)作為分析對象,其典型特征是中間艙段第1 層、第2 層的甲板設(shè)有大開口,其橫截面如圖1所示。三艙段結(jié)構(gòu)總長L=30 m,寬B=18.36 m,型深D=11.3 m,大開口寬度9.12 m,第1 層(頂層)甲板結(jié)構(gòu)布置如圖2 所示。結(jié)構(gòu)材料彈性模量E=200 GPa,泊松比μ=0.3,材料密度ρ=7 800 kg/m3。

    圖1 艙段結(jié)構(gòu)橫截面示意圖Fig.1 The cross-section of the middle cabin structure

    圖2 第1 層甲板結(jié)構(gòu)布置圖Fig.2 The structural layout plan of the first deck

    子模型區(qū)域為連接船舶艙壁垂直桁材與第2層甲板縱桁的肘板節(jié)點結(jié)構(gòu)。大開口艙段的有限元模型不含肘板結(jié)構(gòu)。甲板板、舷側(cè)外板、外底板、肋板、龍骨以及甲板桁材腹板采用板殼單元Shell 181 模擬,子模型區(qū)域外的甲板縱骨、甲板桁材面板和舷側(cè)縱骨等采用梁單元Beam 188 模擬,而子模型區(qū)域內(nèi)的這些結(jié)構(gòu)則全部采用板殼單元Shell 181 模擬。有限元整體模型如圖3 所示。共劃分有188 599 個單元,其中板殼單元135 659 個。

    圖3 三艙段結(jié)構(gòu)有限元整體模型Fig.3 The FEM model of a three-cabin structure

    有限元模型全局坐標(biāo)系為直角坐標(biāo)系,船長方向向艏為X 軸正方向,船寬方向向左為Y 軸正方向,型深方向向上為Z 軸正方向。為了模擬大開口艙段真實的受力狀態(tài),在艙段兩端建立剛性域并在兩端形心處加載463 094.4 kN·m 的彎矩,在第1 層甲板表面加載4.9 kPa 的均布壓力,在第2,3,4 層甲板上加載9.8 kPa 的均布壓力。邊界條件為:在艙段一端的主節(jié)點約束其所有平動自由度和X,Z 方向的轉(zhuǎn)動自由度,釋放Y 方向的轉(zhuǎn)動自由度;另一端的主節(jié)點約束其Y,Z 方向的平動自由度以及X,Z 方向的轉(zhuǎn)動自由度,釋放X 方向的平動和Y 方向的轉(zhuǎn)動。整體模型計算結(jié)果如圖4 和圖5 所示。

    圖4 整體模型應(yīng)力云圖Fig.4 Mises stress contours of full model

    圖5 整體模型變形云圖Fig.5 Displacement contours of full model

    1.2 子模型應(yīng)力計算

    1.2.1 強度計算的有限元子模型

    本文采用子模型法對肘板節(jié)點結(jié)構(gòu)進(jìn)行應(yīng)力分析。子模型法基于圣維南原理,即如果實際分布載荷被等效載荷代替,應(yīng)力和應(yīng)變只在載荷施加位置附近有改變。因此,只要子模型切割邊界避開載荷集中及應(yīng)力集中位置,子模型內(nèi)部就可以得到較精確的解[9]。在應(yīng)力的精細(xì)化分析中,子模型法得到了廣泛應(yīng)用[10-11]。

    本文在ANSYS 和Hyperworks 軟件中建立有限元子模型的步驟如下:

    1)在ANSYS 中生成并分析較粗糙網(wǎng)格的整體模型,保留整體模型.db 文件及結(jié)果.rst文件;

    2)在ANSYS 中生成子模型的幾何模型,在關(guān)注區(qū)域添加肘板結(jié)構(gòu),將幾何模型文件導(dǎo)入Hyperworks 并劃分細(xì)網(wǎng)格,生成節(jié)點;

    3)將在Hyperworks 軟件中生成的肘板節(jié)點結(jié)構(gòu)的有限元子模型導(dǎo)入ANSYS 中,提取并保存子模型邊界節(jié)點,以及后綴為.node 的文件,默認(rèn)文件名與子模型文件名一致;

    4)在ANSYS 子模型工作環(huán)境下恢復(fù)整體模型.db 文件,在后處理模塊讀入整體模型結(jié)果文件,利用命令“cbdof”形成后綴為.cbdo 的子模型邊界插值文件;

    5)在ANSYS 軟件中恢復(fù)子模型.db 文件,在求解模塊讀入上一步生成的.cbdo 文件,并根據(jù)子模型范圍內(nèi)的實際位移邊界條件及外載荷進(jìn)行加載求解。

    子模型的切割邊界需遠(yuǎn)離關(guān)注區(qū)域,邊界的選取也要經(jīng)過嘗試驗證后才能確定。根據(jù)整體模型的應(yīng)力分布,初步選取Z 方向范圍為(6 580~11 300),X 方向范圍為(-6 000~7 250),Y 方向范圍為(-950~950),其在全局坐標(biāo)系下的位置、材料參數(shù)以及單元類型等均與整體模型中一致。

    在子模型的其他部分,采用20 mm 的網(wǎng)格進(jìn)行劃分,并在肘板關(guān)注區(qū)域以5 mm 的網(wǎng)格進(jìn)行劃分。關(guān)注區(qū)域為艙壁垂直桁材與第2 層甲板縱桁連接處,如圖6 所示。加載插值邊界及子模型區(qū)域內(nèi)其他載荷后的子模型有限元模型如圖7所示。

    完成邊界條件的切割、插值以及對子模型的分析后,應(yīng)首先驗證切割的邊界是否離應(yīng)力集中區(qū)域足夠遠(yuǎn)。本文對比了子模型的切割邊界處和整體模型對應(yīng)位置的Mises 應(yīng)力,發(fā)現(xiàn)兩者的應(yīng)力分布與大小基本一致,認(rèn)為子模型的邊界選取較合適。

    圖6 肘板節(jié)點結(jié)構(gòu)子模型Fig.6 FE model of sub-model

    圖7 子模型邊界插值及其加載Fig.7 Boundary interpolation and loading of sub-model

    1.2.2 典型節(jié)點結(jié)構(gòu)子模型法計算結(jié)果分析

    本文關(guān)注的節(jié)點為連接船舶艙壁垂直桁材與第2 層甲板縱桁的肘板節(jié)點,如圖8 所示。為便于表述,圖中桁材面板的淺色區(qū)域稱為區(qū)域2,桁材腹板的淺色區(qū)域稱為區(qū)域3。

    圖8 拓?fù)鋬?yōu)化對象初始形狀Fig.8 Original structure of bracket

    應(yīng)用子模型法計算該處有肘板和無肘板情況下的應(yīng)力分布,計算結(jié)果如圖9 所示。在有肘板的情況下,區(qū)域2 的最大Mises 應(yīng)力為145.0 MPa,出現(xiàn)在肘板與艙壁桁材面板的連接處;區(qū)域3 的最大Mises 應(yīng)力為72.39 MPa,出現(xiàn)在區(qū)域3 的中心處;肘板區(qū)域的最大Mises 應(yīng)力為260.30 MPa,出現(xiàn)在肘板與第2 層甲板縱桁面板的交匯處。在無肘板的情況下,區(qū)域2 的最大Mises 應(yīng)力為207.20 MPa,出現(xiàn)在兩桁材面板的連接處;區(qū)域3的最大Mises 應(yīng)力為83.52 MPa,出現(xiàn)在區(qū)域3 的中心處。由此可知,肘板的存在極大地改善了區(qū)域2 的應(yīng)力集中,但是在肘板的角趾處產(chǎn)生了新的應(yīng)力集中點。

    圖9 有、無肘板應(yīng)力分布結(jié)果Fig.9 Mises stress contours of sub-model with and without bracket

    2 船舶肘板節(jié)點結(jié)構(gòu)的拓?fù)鋬?yōu)化

    2.1 拓?fù)鋬?yōu)化數(shù)學(xué)模型

    本文基于Optistruct 軟件,采用變密度法對船舶肘板節(jié)點結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化。所謂變密度法[12],其實是一種偽密度法,即人為假設(shè)的一種材料物理參數(shù)(如彈性模量)與密度之間的一種線性或者非線性關(guān)系,經(jīng)過拓?fù)鋬?yōu)化后的單元密度一般按照0~1 分布在給定的初始拓?fù)鋬?yōu)化區(qū)域上。通過控制單元的密度向0 或者1 兩端收斂,可以使得拓?fù)鋬?yōu)化結(jié)果抽象成可以加工的結(jié)構(gòu)。

    本文的優(yōu)化對象為圖8 中所示的肘板,設(shè)計變量為肘板結(jié)構(gòu)的單元密度。肘板優(yōu)化設(shè)計的首要目標(biāo)是降低區(qū)域2 的應(yīng)力集中。由于Optistruct軟件在進(jìn)行優(yōu)化求解的過程中不能提取設(shè)置為設(shè)計變量區(qū)域結(jié)構(gòu)的應(yīng)力,但鑒于區(qū)域2 的應(yīng)力水平與肘板自身的應(yīng)力水平具有一定的對應(yīng)關(guān)系,所以將區(qū)域2 的Mises 應(yīng)力最小化來作為優(yōu)化目標(biāo)。為保證桁材的抗彎強度,約束函數(shù)設(shè)定為區(qū)域2 的Mises 應(yīng)力。在有肘板的情況下,區(qū)域3 的最大Mises 應(yīng)力為72.39 MPa,在無肘板的情況下,區(qū)域3 的最大Mises 應(yīng)力為83.52 MPa,因此,為保證桁材的抗彎強度,本文在進(jìn)行肘板優(yōu)化設(shè)計時限定區(qū)域3 的Mises 應(yīng)力不得大于73 MPa。

    2.2 拓?fù)鋬?yōu)化結(jié)果及工程化處理

    求解上述數(shù)學(xué)模型,至優(yōu)化迭代20 步時停止計算。為便于觀察,隱藏非設(shè)計區(qū)域單元,只顯示設(shè)計區(qū)域單元,優(yōu)化后,節(jié)點結(jié)構(gòu)的密度云圖如圖10 所示。目標(biāo)函數(shù)隨迭代步數(shù)的變化歷程曲線如圖11 所示,圖中,橫坐標(biāo)為迭代步數(shù),縱坐標(biāo)為目標(biāo)函數(shù)的數(shù)值。由圖11 可以看出,最后的目標(biāo)函數(shù)已經(jīng)收斂,說明肘板的優(yōu)化在該數(shù)學(xué)模型下已經(jīng)達(dá)到最優(yōu)。

    圖10 肘板拓?fù)鋬?yōu)化密度云圖Fig.10 Density contours of topology optimization

    圖11 目標(biāo)函數(shù)迭代歷程曲線Fig.11 Objective function--time process of topology optimization

    肘板最后一步迭代的應(yīng)力云圖如圖12 所示。肘板自身最大Mises 應(yīng)力隨迭代步數(shù)的變化歷程曲線如圖13 所示,圖中,橫坐標(biāo)為迭代步數(shù),縱坐標(biāo)為肘板的最大Mises 應(yīng)力。

    對應(yīng)于優(yōu)化的最終迭代步數(shù),肘板區(qū)域的最大Mises 應(yīng)力為194.20 MPa,區(qū)域3 的最大Mises應(yīng)力為72.66 MPa,區(qū)域2 的最大Mises 應(yīng)力為102.98 MPa。密度云圖中,越接近紅色(深色部分),就表示這部分的材料越重要,越接近藍(lán)色(淺色部分),就表示這部分的材料越可以舍棄。綜合以上結(jié)果,可以對肘板的拓?fù)鋬?yōu)化結(jié)果進(jìn)行定性分析,認(rèn)為肘板的3 個尖角是引起桁材面板和肘板自身應(yīng)力集中的主要因素。

    圖12 拓?fù)鋬?yōu)化應(yīng)力云圖Fig.12 Mises stress contours of topology optimization

    圖13 肘板應(yīng)力迭代歷程曲線Fig.13 Bracket stress--time process of topology optimization

    根據(jù)圖10 所示的肘板拓?fù)鋬?yōu)化結(jié)果及分析,對肘板結(jié)構(gòu)進(jìn)行了如下工程化處理,并在肘板的3 個尖角處設(shè)置了相同半徑的圓弧開口,如圖14所示。采用前述的子模型分析法,分別對圓弧半徑取10~46 mm 的方案進(jìn)行肘板節(jié)點結(jié)構(gòu)應(yīng)力計算。

    圖14 修改后的肘板結(jié)構(gòu)型式Fig.14 New bracket structure type

    應(yīng)用上述子模型法對不同圓弧半徑的肘板進(jìn)行計算,結(jié)果如表1 所示。在計算過程中,為避免網(wǎng)格劃分對結(jié)果的影響,模型的網(wǎng)格密度、網(wǎng)格質(zhì)量以及網(wǎng)格參數(shù)等均保持一致。

    表1 不同圓弧半徑時應(yīng)力結(jié)果Tab.1 Mises stress results of model with different radii

    表1 的結(jié)果表明,當(dāng)圓弧半徑較大時,在保證區(qū)域2 和區(qū)域3 的應(yīng)力水平變化不大的情況下,肘板自身的應(yīng)力水平與傳統(tǒng)的三角形肘板相比有所下降。當(dāng)圓弧半徑為43 mm 時,對肘板自身的應(yīng)力降低效果最好,和傳統(tǒng)的三角形肘板相比,應(yīng)力降低了16.56%。

    3 結(jié) 論

    本文通過引入子模型技術(shù),對船舶肘板節(jié)點結(jié)構(gòu)進(jìn)行了應(yīng)力精細(xì)化分析,并在此基礎(chǔ)上,以肘板結(jié)構(gòu)的密度為優(yōu)化變量,以與肘板相連的桁材面板的Mises 應(yīng)力為目標(biāo),以與肘板相連的桁材腹板的Mises 應(yīng)力為約束條件,建立了合理的拓?fù)鋬?yōu)化數(shù)學(xué)模型并進(jìn)行了肘板優(yōu)化設(shè)計,在根據(jù)肘板拓?fù)鋬?yōu)化結(jié)果進(jìn)行適當(dāng)?shù)墓こ袒幚砗?,提出了一種新型的肘板結(jié)構(gòu)型式,得到如下結(jié)論:

    1)肘板的存在較好地改善了桁材面板交匯處的應(yīng)力集中。對于本文的算例,桁材面板交匯處的最大Mises 應(yīng)力由207.2 MPa 降低到了145.0 MPa,降低了30.02%。但肘板的存在又引起了新的應(yīng)力集中點,肘板與桁材面板交匯處的應(yīng)力達(dá)到了260.30 MPa。

    2)基于肘板拓?fù)鋬?yōu)化結(jié)果提出的新型肘板結(jié)構(gòu)型式是一種在傳統(tǒng)三角形肘板的3 個尖角處設(shè)置1/4 圓弧的結(jié)構(gòu)。該種肘板在本文所采用的典型載荷下有效降低了節(jié)點結(jié)構(gòu)的應(yīng)力集中。對于本文的算例,肘板與桁材面板交匯處的應(yīng)力由260.30 MPa 降低到了217.20 MPa,降低了16.56%。所提出的肘板結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計方法為獲得創(chuàng)新性的肘板結(jié)構(gòu)型式提供了一種新的手段與方法。

    [1]JORDAN C R,COCHRAN C S. In service performance of structural details:SSC-272[R].Ship Structure Com?mittee,US,1978:154-159.

    [2]KIM T W,LIM S S,SEOK H H,et al. Concurrent engi?neering solution for the design of ship and offshore bracket parts and fabrication process[J]. International Journal of Naval Architecture and Ocean Engineering,2013,5(3):376-391.

    [3]郭信川,吳俊仁,陳慶忠. 艙口圍緣之鍵槽端肘板構(gòu)件最佳化之研究[J]. 中國造船暨輪機工程學(xué)刊,2009,8(3):131-142.KUO H C,WU J L,CHEN C J. Optimization design of a match-stick-hole end-bracket device of a hatch[J].Journal of Taiwan Society of Naval Architects and Ma?rine Engineers,2009,8(3):131-142.

    [4]LIM J B P,NETHERCOT D A. Design and develop?ment of a general cold-formed steel portal framing sys?tem[J]. The Structural Engineer,2002,80(21):31-40.

    [5]姜以威,王廣戈.梁肘板節(jié)點的應(yīng)力集中研究[J].鎮(zhèn)江船舶學(xué)院學(xué)報,1990,4(4):47-56.

    [6]王波,楊平.船舶結(jié)構(gòu)節(jié)點的承載力分析[J].船海工程,2010,39(2):18-21.WANG Bo,YANG Ping. Analysis of the load bearing capacity of ship structural details[J]. Ship and Ocean Engineering,2010,39(2):18-21.

    [7]史戰(zhàn)新.基于Ansys 子模型法的肘板結(jié)構(gòu)優(yōu)化[J].艦船科學(xué)技術(shù),2014,36(8):19-26.SHI Zhanxin. The bracket structure optimization based on Ansys sub-model method[J]. Ship Science and Technology,2014,36(8):19-26.

    [8]田旭軍,閆國強,胡剛義. 基于遺傳算法的肘板結(jié)構(gòu)型式優(yōu)化[J].中國艦船研究,2007,2(2):23-26.TIAN Xujun,YAN Guoqiang,HU Gangyi. Optimiza?tion of bracket plate using genetic algorithm[J]. Chi?nese Journal of Ship Research,2007,2(2):23-26.

    [9]夏偉,胡成,瞿爾仁.Ansys 子模型分析技術(shù)在處理應(yīng)力集中時的應(yīng)用[J]. 工程與建設(shè),2006,20(2):92-94.XIA Wei,HU Cheng,QU Erren. Applications of Ansys submodeling technology in dealing with stress concen?tration[J]. Engineering and Construction,2006,20(2):92-94.

    [10]朱揚,程遠(yuǎn)勝,劉均.激光焊接夾層甲板板格強度計算的子模型方法[J]. 船舶力學(xué),2014,18(10):1228-1236.ZHU Yang,CHENG Yuansheng,LIU Jun. Sub-model method for strength calculation of a laser-welded steel sandwich panel structure[J].Journal of Ship Me?chanics,2014,18(10):1228-1236.

    [11]王虎,程遠(yuǎn)勝,劉均. I 型金屬夾層結(jié)構(gòu)連接構(gòu)件強度數(shù)值計算方法[J]. 中國艦船研究,2012,7(3):51-56.WANG Hu,CHENG Yuansheng,LIU Jun. Strength analysis on model I-core steel sandwich panel joints[J]. Chinese Journal of Ship Research,2012,7(3):51-56.

    [12]朱劍鋒,林逸,陳瀟凱,等. 汽車變速箱殼體結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計[J]. 吉林大學(xué)學(xué)報(工學(xué)版),2013,43(3):584-589.ZHU Jianfeng,LIN Yi,CHEN Xiaokai,et al. Struc?tural topology optimization based design of automotive transmission housing structure[J]. Journal of Jilin University(Engineering and Technology Edition),2013,43(3):584-589.

    猜你喜歡
    肘板板結(jié)構(gòu)甲板
    環(huán)氧樹脂基復(fù)合材料加筋板結(jié)構(gòu)吸濕行為研究
    礦用顎式破碎機肘板結(jié)構(gòu)的改進(jìn)分析
    科考船木甲板安裝工藝
    基于子模型細(xì)化分析的VLCC底邊艙上折角半檔肘板選型
    殲-15艦載戰(zhàn)斗機在遼寧艦飛行甲板準(zhǔn)備起飛
    國防(2017年12期)2017-04-09 07:52:09
    空氣管及通風(fēng)管的應(yīng)力分析
    替代甲板結(jié)合負(fù)壓吸引治療修復(fù)甲床缺損
    甲板前置治療甲床部分缺損
    板結(jié)構(gòu)-聲場耦合分析的FE-LSPIM/FE法
    彈性邊界板結(jié)構(gòu)拓?fù)鋬?yōu)化聲輻射
    国产午夜福利久久久久久| 啦啦啦 在线观看视频| 亚洲无线在线观看| 一本久久中文字幕| 老司机午夜福利在线观看视频| 中文字幕久久专区| 精品电影一区二区在线| 亚洲精品美女久久av网站| 麻豆av在线久日| 热99re8久久精品国产| 久久影院123| 老司机福利观看| 中文字幕高清在线视频| 国产成人欧美| 国产一区二区在线av高清观看| 国产精华一区二区三区| 一级毛片女人18水好多| 亚洲专区中文字幕在线| 久久久精品欧美日韩精品| a在线观看视频网站| √禁漫天堂资源中文www| 香蕉久久夜色| 国产成人系列免费观看| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 色老头精品视频在线观看| 亚洲av成人一区二区三| 欧美激情高清一区二区三区| 久久久久国产一级毛片高清牌| netflix在线观看网站| ponron亚洲| 色哟哟哟哟哟哟| 国产精品,欧美在线| 最近最新免费中文字幕在线| 国产精品永久免费网站| 欧美不卡视频在线免费观看 | av免费在线观看网站| 亚洲情色 制服丝袜| 亚洲精品av麻豆狂野| 欧美不卡视频在线免费观看 | 99国产综合亚洲精品| 男女之事视频高清在线观看| 18禁裸乳无遮挡免费网站照片 | 高清黄色对白视频在线免费看| 久久久久九九精品影院| 热99re8久久精品国产| 美女扒开内裤让男人捅视频| 亚洲情色 制服丝袜| 美国免费a级毛片| 国产亚洲av高清不卡| 亚洲精品一区av在线观看| 91精品三级在线观看| 国产精品久久视频播放| 制服诱惑二区| 免费观看精品视频网站| 国产精品久久久久久亚洲av鲁大| 搡老妇女老女人老熟妇| www.熟女人妻精品国产| 身体一侧抽搐| 国产一区二区激情短视频| 欧美日韩一级在线毛片| 国产亚洲精品av在线| 亚洲精品在线观看二区| 满18在线观看网站| 欧美久久黑人一区二区| 黄色 视频免费看| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 免费久久久久久久精品成人欧美视频| 日韩精品中文字幕看吧| 精品国产美女av久久久久小说| 麻豆av在线久日| 国产日韩一区二区三区精品不卡| 亚洲精品在线美女| 久久久久久大精品| 老鸭窝网址在线观看| 一边摸一边抽搐一进一小说| 久久久久久久精品吃奶| 亚洲精华国产精华精| 色精品久久人妻99蜜桃| 国产国语露脸激情在线看| 国产精品久久电影中文字幕| 国产成人免费无遮挡视频| 777久久人妻少妇嫩草av网站| 天天躁夜夜躁狠狠躁躁| 少妇被粗大的猛进出69影院| 天天添夜夜摸| 一区二区三区精品91| 女人精品久久久久毛片| 久久久久久久久免费视频了| 变态另类成人亚洲欧美熟女 | 一进一出抽搐动态| 色综合婷婷激情| 亚洲美女黄片视频| a级毛片在线看网站| 亚洲天堂国产精品一区在线| 咕卡用的链子| 黄网站色视频无遮挡免费观看| 波多野结衣av一区二区av| 亚洲中文av在线| 欧美人与性动交α欧美精品济南到| aaaaa片日本免费| 国产激情欧美一区二区| 欧美激情久久久久久爽电影 | 欧美黄色淫秽网站| 日日摸夜夜添夜夜添小说| 国产精品亚洲一级av第二区| 亚洲伊人色综图| 亚洲av成人av| 操出白浆在线播放| av电影中文网址| 日本欧美视频一区| 电影成人av| 亚洲第一电影网av| 99久久久亚洲精品蜜臀av| av有码第一页| 狂野欧美激情性xxxx| 一边摸一边抽搐一进一出视频| 国产av一区二区精品久久| 嫁个100分男人电影在线观看| 黄色丝袜av网址大全| 欧美乱色亚洲激情| 丝袜在线中文字幕| 欧美+亚洲+日韩+国产| 国产欧美日韩一区二区三| av视频在线观看入口| 久99久视频精品免费| 制服人妻中文乱码| 超碰成人久久| 久久婷婷成人综合色麻豆| 国产欧美日韩一区二区三| 久久久久久久精品吃奶| 亚洲成a人片在线一区二区| 啦啦啦观看免费观看视频高清 | 亚洲精品中文字幕一二三四区| 久久这里只有精品19| 淫妇啪啪啪对白视频| 老司机在亚洲福利影院| 精品久久久久久,| 日韩视频一区二区在线观看| 久99久视频精品免费| 免费在线观看日本一区| 一边摸一边做爽爽视频免费| 国产欧美日韩精品亚洲av| x7x7x7水蜜桃| 亚洲一区中文字幕在线| 天堂√8在线中文| 99久久精品国产亚洲精品| 国产精品乱码一区二三区的特点 | 黄色毛片三级朝国网站| 国产xxxxx性猛交| 18禁裸乳无遮挡免费网站照片 | 国产精品国产高清国产av| 国产欧美日韩一区二区三| 别揉我奶头~嗯~啊~动态视频| 国产熟女xx| av在线播放免费不卡| 嫁个100分男人电影在线观看| 看片在线看免费视频| 久久人人精品亚洲av| 亚洲黑人精品在线| 久久国产亚洲av麻豆专区| 午夜福利18| 电影成人av| 久久青草综合色| 色在线成人网| 日韩大码丰满熟妇| 国产成人免费无遮挡视频| 国产又爽黄色视频| 麻豆国产av国片精品| 香蕉久久夜色| 久久精品国产清高在天天线| 欧美丝袜亚洲另类 | 久久午夜亚洲精品久久| 无限看片的www在线观看| 欧美午夜高清在线| 狂野欧美激情性xxxx| 国产视频一区二区在线看| av有码第一页| 欧美大码av| 嫩草影视91久久| 69av精品久久久久久| 夜夜躁狠狠躁天天躁| 国产野战对白在线观看| 十八禁人妻一区二区| 午夜日韩欧美国产| 人人妻人人爽人人添夜夜欢视频| 国产单亲对白刺激| 亚洲中文字幕一区二区三区有码在线看 | 久久香蕉激情| 久久精品亚洲熟妇少妇任你| 琪琪午夜伦伦电影理论片6080| 亚洲欧美精品综合久久99| 激情视频va一区二区三区| 黑丝袜美女国产一区| 亚洲国产精品久久男人天堂| 亚洲专区中文字幕在线| 一级片免费观看大全| 精品久久久久久久人妻蜜臀av | 久久这里只有精品19| 777久久人妻少妇嫩草av网站| 亚洲精品国产精品久久久不卡| 大码成人一级视频| 国产成人欧美| 男女下面插进去视频免费观看| 亚洲欧美日韩高清在线视频| 婷婷六月久久综合丁香| 十分钟在线观看高清视频www| 国产精品亚洲一级av第二区| 亚洲精品中文字幕一二三四区| 我的亚洲天堂| 一区二区日韩欧美中文字幕| 可以在线观看的亚洲视频| 亚洲av五月六月丁香网| 日韩欧美免费精品| 久久久久久久久免费视频了| 日韩欧美一区视频在线观看| 久久久久久亚洲精品国产蜜桃av| 麻豆久久精品国产亚洲av| 性色av乱码一区二区三区2| 91老司机精品| 美女大奶头视频| 精品欧美国产一区二区三| 欧美绝顶高潮抽搐喷水| 人人妻人人爽人人添夜夜欢视频| 大码成人一级视频| 老熟妇仑乱视频hdxx| www日本在线高清视频| 88av欧美| √禁漫天堂资源中文www| 嫩草影视91久久| 视频在线观看一区二区三区| 欧美日韩瑟瑟在线播放| av在线天堂中文字幕| 精品第一国产精品| av视频免费观看在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品亚洲一级av第二区| 日韩欧美免费精品| 在线国产一区二区在线| 一二三四社区在线视频社区8| 精品久久蜜臀av无| 女人高潮潮喷娇喘18禁视频| tocl精华| 一边摸一边做爽爽视频免费| 老鸭窝网址在线观看| 色播在线永久视频| 欧美黄色片欧美黄色片| 亚洲精品国产色婷婷电影| av中文乱码字幕在线| 成人手机av| 亚洲最大成人中文| 中文字幕另类日韩欧美亚洲嫩草| 女人被狂操c到高潮| 色综合欧美亚洲国产小说| 黑人巨大精品欧美一区二区蜜桃| 19禁男女啪啪无遮挡网站| 精品无人区乱码1区二区| 波多野结衣巨乳人妻| 欧美av亚洲av综合av国产av| 色播亚洲综合网| 欧美久久黑人一区二区| 村上凉子中文字幕在线| 色精品久久人妻99蜜桃| 久久久久久久久中文| 中文字幕人成人乱码亚洲影| www.www免费av| 亚洲国产欧美网| 婷婷六月久久综合丁香| x7x7x7水蜜桃| 亚洲人成网站在线播放欧美日韩| 成人亚洲精品一区在线观看| 久久伊人香网站| 一区二区三区高清视频在线| 午夜亚洲福利在线播放| 最新美女视频免费是黄的| 日韩欧美国产一区二区入口| 中文字幕另类日韩欧美亚洲嫩草| 多毛熟女@视频| 色哟哟哟哟哟哟| 免费看美女性在线毛片视频| 人妻久久中文字幕网| 亚洲精品国产区一区二| 99国产综合亚洲精品| 成人国语在线视频| 老司机在亚洲福利影院| 精品乱码久久久久久99久播| 国产精品美女特级片免费视频播放器 | 热99re8久久精品国产| 中文字幕精品免费在线观看视频| 一a级毛片在线观看| 91国产中文字幕| 中文字幕久久专区| 男女午夜视频在线观看| a在线观看视频网站| 波多野结衣巨乳人妻| 搡老岳熟女国产| 给我免费播放毛片高清在线观看| 咕卡用的链子| 一进一出好大好爽视频| 国产麻豆成人av免费视频| 日韩国内少妇激情av| 日韩欧美一区视频在线观看| 制服诱惑二区| 9191精品国产免费久久| 久久狼人影院| 日韩大码丰满熟妇| 人人澡人人妻人| 亚洲成人精品中文字幕电影| 两人在一起打扑克的视频| 悠悠久久av| 午夜免费激情av| 国产精品国产高清国产av| 日韩欧美免费精品| 国产免费男女视频| 国产亚洲精品久久久久5区| 中文字幕高清在线视频| 91成年电影在线观看| 欧美最黄视频在线播放免费| 久久久水蜜桃国产精品网| 91大片在线观看| 婷婷六月久久综合丁香| 国内久久婷婷六月综合欲色啪| 久久国产精品男人的天堂亚洲| aaaaa片日本免费| 亚洲全国av大片| 久久久国产成人免费| 久久精品91蜜桃| 麻豆av在线久日| 精品福利观看| 18美女黄网站色大片免费观看| 色老头精品视频在线观看| 午夜亚洲福利在线播放| 亚洲va日本ⅴa欧美va伊人久久| 国产精品精品国产色婷婷| 一级毛片高清免费大全| 久久久水蜜桃国产精品网| 我的亚洲天堂| 亚洲精品美女久久av网站| 久久午夜综合久久蜜桃| 88av欧美| 国产99久久九九免费精品| www.999成人在线观看| av天堂久久9| 少妇被粗大的猛进出69影院| 亚洲天堂国产精品一区在线| 在线永久观看黄色视频| 国产人伦9x9x在线观看| 色综合亚洲欧美另类图片| 日韩欧美在线二视频| 国产不卡一卡二| 色综合欧美亚洲国产小说| 美女免费视频网站| 久久久久国产一级毛片高清牌| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩亚洲国产一区二区在线观看| 一级a爱片免费观看的视频| 亚洲成国产人片在线观看| 村上凉子中文字幕在线| 欧美一级毛片孕妇| avwww免费| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品男人的天堂亚洲| 精品久久蜜臀av无| 亚洲片人在线观看| 欧美黄色淫秽网站| 亚洲国产欧美日韩在线播放| 亚洲av熟女| 精品国内亚洲2022精品成人| 日本五十路高清| 日本 欧美在线| 国产精品电影一区二区三区| 国产男靠女视频免费网站| 国产高清有码在线观看视频 | 欧美国产精品va在线观看不卡| 黄色视频不卡| 亚洲精品国产一区二区精华液| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久av网站| 一卡2卡三卡四卡精品乱码亚洲| 国产精品综合久久久久久久免费 | 色尼玛亚洲综合影院| 好男人电影高清在线观看| 啪啪无遮挡十八禁网站| 在线国产一区二区在线| 91麻豆精品激情在线观看国产| 免费不卡黄色视频| 在线观看www视频免费| 狠狠狠狠99中文字幕| 麻豆一二三区av精品| 在线永久观看黄色视频| 国产97色在线日韩免费| 久久久久久免费高清国产稀缺| 免费人成视频x8x8入口观看| 国产1区2区3区精品| 久久婷婷人人爽人人干人人爱 | 波多野结衣av一区二区av| 亚洲色图综合在线观看| 免费av毛片视频| 狂野欧美激情性xxxx| 日本欧美视频一区| 美女扒开内裤让男人捅视频| 亚洲成人国产一区在线观看| 啦啦啦韩国在线观看视频| 久久香蕉国产精品| 欧美乱码精品一区二区三区| 久久精品国产亚洲av高清一级| 激情在线观看视频在线高清| 日韩高清综合在线| 国产高清视频在线播放一区| 操美女的视频在线观看| 91成人精品电影| 国产单亲对白刺激| 婷婷六月久久综合丁香| 最新美女视频免费是黄的| 国产精品二区激情视频| 亚洲中文日韩欧美视频| 18禁美女被吸乳视频| 亚洲一区中文字幕在线| 大香蕉久久成人网| 欧美一区二区精品小视频在线| 人人妻人人澡人人看| 亚洲少妇的诱惑av| 成人亚洲精品av一区二区| 午夜福利一区二区在线看| 精品第一国产精品| 视频区欧美日本亚洲| 香蕉久久夜色| 电影成人av| 一区二区三区激情视频| 一进一出抽搐动态| 国产成人系列免费观看| 天堂影院成人在线观看| 久久香蕉国产精品| 老司机深夜福利视频在线观看| av视频在线观看入口| 又紧又爽又黄一区二区| 一区在线观看完整版| 成人国产一区最新在线观看| 中出人妻视频一区二区| 亚洲视频免费观看视频| 黄色毛片三级朝国网站| 精品免费久久久久久久清纯| 国产视频一区二区在线看| 国产欧美日韩一区二区精品| 亚洲少妇的诱惑av| 久久伊人香网站| 中文字幕人妻丝袜一区二区| 成人欧美大片| 少妇裸体淫交视频免费看高清 | 国产高清videossex| 99精品欧美一区二区三区四区| 精品久久久久久久人妻蜜臀av | 在线永久观看黄色视频| 欧美国产日韩亚洲一区| 老汉色av国产亚洲站长工具| 国产成+人综合+亚洲专区| 啪啪无遮挡十八禁网站| 黄色 视频免费看| 性少妇av在线| 日韩精品免费视频一区二区三区| 国产在线精品亚洲第一网站| 国产一区二区激情短视频| 成人永久免费在线观看视频| 久久精品91无色码中文字幕| ponron亚洲| 少妇的丰满在线观看| 日韩欧美一区二区三区在线观看| 窝窝影院91人妻| 少妇裸体淫交视频免费看高清 | 九色国产91popny在线| 久久香蕉国产精品| www日本在线高清视频| 18禁裸乳无遮挡免费网站照片 | 日韩精品免费视频一区二区三区| 淫妇啪啪啪对白视频| 精品久久久久久久毛片微露脸| 色婷婷久久久亚洲欧美| 黄片大片在线免费观看| 人人妻,人人澡人人爽秒播| 色av中文字幕| 天堂√8在线中文| 午夜福利影视在线免费观看| 99riav亚洲国产免费| 天天躁夜夜躁狠狠躁躁| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av美国av| 咕卡用的链子| 又大又爽又粗| 亚洲伊人色综图| 国产成人啪精品午夜网站| 一二三四社区在线视频社区8| 美国免费a级毛片| 欧美成人一区二区免费高清观看 | 午夜福利成人在线免费观看| 亚洲精品美女久久久久99蜜臀| 亚洲av成人一区二区三| 色精品久久人妻99蜜桃| 嫩草影视91久久| 亚洲精品国产一区二区精华液| 十八禁网站免费在线| 又黄又粗又硬又大视频| 91国产中文字幕| 国产亚洲精品第一综合不卡| 亚洲天堂国产精品一区在线| 嫩草影院精品99| 久热爱精品视频在线9| 国产麻豆成人av免费视频| 在线观看一区二区三区| 免费在线观看日本一区| 最近最新中文字幕大全免费视频| 999精品在线视频| 久久久国产成人精品二区| 亚洲第一青青草原| 国产熟女xx| 又大又爽又粗| 国产男靠女视频免费网站| 国产片内射在线| 久久久国产精品麻豆| 精品福利观看| 无人区码免费观看不卡| 两个人免费观看高清视频| 看片在线看免费视频| 国产精品久久久久久人妻精品电影| 不卡一级毛片| 精品日产1卡2卡| 国产精品一区二区三区四区久久 | 亚洲欧美激情在线| aaaaa片日本免费| 国产精品久久久久久精品电影 | 亚洲精品一区av在线观看| 黄色丝袜av网址大全| or卡值多少钱| 中文亚洲av片在线观看爽| 两人在一起打扑克的视频| 18禁国产床啪视频网站| 午夜福利影视在线免费观看| 成人手机av| 国产精品久久久av美女十八| 中文字幕久久专区| 人成视频在线观看免费观看| 午夜精品国产一区二区电影| 熟女少妇亚洲综合色aaa.| 亚洲精品美女久久av网站| 久久久久国产精品人妻aⅴ院| 别揉我奶头~嗯~啊~动态视频| 久久狼人影院| 一区二区三区国产精品乱码| 久久精品91无色码中文字幕| 日本五十路高清| 国产日韩一区二区三区精品不卡| 免费观看人在逋| 欧美亚洲日本最大视频资源| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品粉嫩美女一区| 九色亚洲精品在线播放| 日韩一卡2卡3卡4卡2021年| 精品国产美女av久久久久小说| 久久精品国产综合久久久| 男人的好看免费观看在线视频 | 神马国产精品三级电影在线观看 | 国产精品久久电影中文字幕| 久久久久久大精品| 一边摸一边做爽爽视频免费| 性欧美人与动物交配| 久久精品影院6| 老司机在亚洲福利影院| 极品人妻少妇av视频| 精品乱码久久久久久99久播| 亚洲va日本ⅴa欧美va伊人久久| 午夜精品久久久久久毛片777| av欧美777| 国产成人精品久久二区二区免费| 亚洲人成电影观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲色图综合在线观看| 午夜福利在线观看吧| 午夜影院日韩av| 精品国产乱子伦一区二区三区| 两个人视频免费观看高清| 男人舔女人的私密视频| 久久午夜综合久久蜜桃| 亚洲七黄色美女视频| 男人舔女人的私密视频| 国产精品一区二区在线不卡| 18禁美女被吸乳视频| 999久久久国产精品视频| 日韩欧美国产在线观看| 少妇裸体淫交视频免费看高清 | 久久久久久久久中文| 一本久久中文字幕| 啦啦啦免费观看视频1| 一区二区三区国产精品乱码| 黄片大片在线免费观看| 国产又色又爽无遮挡免费看| 亚洲av成人一区二区三| tocl精华| 国产精品久久久人人做人人爽| 亚洲国产日韩欧美精品在线观看 | 看片在线看免费视频| 日韩高清综合在线| 免费少妇av软件| 丝袜美足系列| tocl精华| 国产精品久久久人人做人人爽| 无限看片的www在线观看| 免费在线观看黄色视频的| 亚洲成人国产一区在线观看| 久久久久国产精品人妻aⅴ院| 精品无人区乱码1区二区| 国产精品二区激情视频| 国产精品一区二区精品视频观看| 欧美黄色淫秽网站| 美女高潮到喷水免费观看| 亚洲 国产 在线|