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

    不確定載荷下剛架體結(jié)構(gòu)穩(wěn)健優(yōu)化設(shè)計

    2015-08-23 09:37:06楊繼超戴宗妙
    關(guān)鍵詞:剛架穩(wěn)健性桿件

    楊繼超,戴宗妙

    (1.哈爾濱工程大學(xué)機電工程學(xué)院,黑龍江哈爾濱150001;2.中國船舶重工集團公司第七一三研究所,河南鄭州450015)

    優(yōu)化設(shè)計是指在特定約束條件下,獲得結(jié)構(gòu)某種性能指標(biāo)(如體積、質(zhì)量或擾度等)最優(yōu)。傳統(tǒng)的優(yōu)化設(shè)計,往往建立在確定性模型基礎(chǔ)上,沒有考慮變量參數(shù)實際值與名義值存在的差異,得到的優(yōu)化解位于可行域邊界。當(dāng)實際工程中變量參數(shù)(如載荷)發(fā)生擾動時,理論最優(yōu)解很可能落在不可行區(qū)域,給設(shè)計帶來潛在風(fēng)險。針對工程設(shè)計中不確定性,國內(nèi)外學(xué)者提出許多設(shè)計理論和方法。其中穩(wěn)健設(shè)計基本思想是在不消除不確定性源的前提下,提高產(chǎn)品對不確定因素抵抗力[1]。結(jié)構(gòu)穩(wěn)健優(yōu)化一個難點是對不確定性力的分析及對優(yōu)化過程影響,即建立不確定性優(yōu)化模型。

    目前,基于概率不確定模型應(yīng)用比較成熟。文獻[2-4]通過構(gòu)造約束和目標(biāo)函數(shù)靈敏度分析,進行穩(wěn)健優(yōu)化設(shè)計;文獻[5]采用Kriging模型,建立以目標(biāo)性能指標(biāo)均值和方差為優(yōu)化目的優(yōu)化模型;文獻[6]將穩(wěn)健優(yōu)化問題轉(zhuǎn)化成多目標(biāo)優(yōu)化問題,對結(jié)構(gòu)進行優(yōu)化設(shè)計;為了提高計算效率,文獻[7]基于隨機響應(yīng)量前四階矩建立穩(wěn)健優(yōu)化設(shè)計模型?;诟怕誓P头椒ㄐ枰罅繕颖緮?shù)據(jù)估算參數(shù)的不確定度或分布函數(shù),而這些數(shù)據(jù)在實際工程中往往很難實現(xiàn)。然而獲得不確定因素的邊界,是相對容易的。非概率方法就是根據(jù)不確定參數(shù)未知但有界的特點進行建模。文獻[8-10]提出了基于凸集理論的非概率不確定模型分析方法,對結(jié)構(gòu)進行穩(wěn)健性分析。文獻[11]在非概率凸模型可靠性約束下,通過超橢球域界定不確定參量,進行結(jié)構(gòu)優(yōu)化設(shè)計。文獻[12-13]運用區(qū)間模型對不確定性參數(shù)進行非概率性描述,建立結(jié)構(gòu)優(yōu)化和設(shè)計方法。目前,在集合理論凸模型中,大多方法都是求解結(jié)構(gòu)響應(yīng)上下界解析解的近似值[14];對于復(fù)雜的結(jié)構(gòu),區(qū)間模型計算量往往非常大,計算效率低下甚至無法實行。

    20世紀(jì)90年代,Ben-Haim提出非概率Information-Gap(I-G)理論描述不確定性模型,并以系統(tǒng)承受最大不確定性波動作為評價結(jié)構(gòu)穩(wěn)健性的指標(biāo)。文獻[15-16]給出了非概率的I-G模型處理不確定性的理論方法,文獻[17-18]基于I-G理論,對地下結(jié)構(gòu)在非確定因素影響下結(jié)構(gòu)的穩(wěn)健性進行分析。為此,本文針對剛架體結(jié)構(gòu)承受不確定載荷特點,引入I-G理論,建立結(jié)構(gòu)非概率穩(wěn)健性模型,并結(jié)合優(yōu)化設(shè)計方法,對剛架體結(jié)構(gòu)進行穩(wěn)健性優(yōu)化設(shè)計,提高對不確定因素的抵抗能力。

    1 不確定載荷I-G模型

    1.1 剛架體應(yīng)力約束方程

    在實際工程中,剛架體結(jié)構(gòu)通常會受到不確定載荷作用,且結(jié)構(gòu)需要滿足最基本應(yīng)力約束要求。令一個線彈性的三維剛架體結(jié)構(gòu)自由度為n,則系統(tǒng)平衡方程為

    式中:u∈Rn和f∈Rn分別表示節(jié)點的位移和外載荷矢量,K∈Sn為剛度矩陣。對于有m個桿件的剛架體結(jié)構(gòu),桿截面積a=(ai)∈Rm,則K可以寫成如下形式:

    式中:b=(bij)∈Rn是常向量。則由上述關(guān)系式,應(yīng)力約束方程為

    式中:E表示彈性模量,表示第i桿許用應(yīng)力。根據(jù)位移表示約束方程,則有

    式中:dl∈R,常量vl∈R。為了保證剛架體結(jié)構(gòu)的安全穩(wěn)定,桿件必須滿足式(3)或式(4)。當(dāng)上述約束關(guān)系式滿足時,不確定變量允許一定的變化幅度。這里通過穩(wěn)健函數(shù)來度量,它與影響剛度矩陣的桿截面積和不確定程度相關(guān)。

    1.2 不確定性因素I-G分析

    I-G模型是由凸集理論基礎(chǔ)發(fā)展而來。根據(jù)I-G理論,不確定度是采用已知信息與未知信息的偏離程度或斷缺狀態(tài)描述。由于外在環(huán)境的不可控性以及所需解決問題信息掌握貧乏,結(jié)合剛架體結(jié)構(gòu)特點,建立不確定性影響因素的I-G模型。不確定參數(shù)不確定性程度可表達為[17]

    式(6)表明,當(dāng)不考慮載荷不確定性(α=0)時,名義載荷即為實際載荷f,I-G模型退化成確定性模型。式(7)表明,不確定參數(shù)α決定不確定因素真實值和名義值之間的偏離程度,當(dāng)α小時,和f越接近;反之偏離程度大。相比較其他非概率模型,如固定區(qū)間模型或凸集模型,I-G模型U(α,)的嵌套凸集特性能更好反應(yīng)不確定性載荷的動態(tài)變化。

    1.3 基于I-G模型穩(wěn)健函數(shù)構(gòu)造

    穩(wěn)健性反映結(jié)構(gòu)對不確定性因素的抵抗程度,如果結(jié)構(gòu)能夠允許較大不確定性變化而不失效,就認(rèn)為該結(jié)構(gòu)穩(wěn)健強。上一節(jié),建立了不確定載荷的I-G模型,定義參數(shù)α表示不確定性變化幅度。

    在滿足約束條件式(3)或(4)的前提下,結(jié)合I-G理論,剛架體結(jié)構(gòu)系統(tǒng)的平衡方程(1)所有可能解集可表示成

    結(jié)合式(4)、(8),式(9)的最優(yōu)規(guī)劃問題可表示為

    式(10)即為穩(wěn)健函數(shù)。對于兩組不同剛架體截面積向量a1∈Rn和a2∈Rn,如果則a1比a2穩(wěn)健。

    2 結(jié)構(gòu)穩(wěn)健優(yōu)化模型

    2.1 一般優(yōu)化模型形式

    傳統(tǒng)的確定性問題數(shù)學(xué)優(yōu)化模型可以表示為

    式中:f(x)為優(yōu)化模型的目標(biāo)函數(shù),x優(yōu)化設(shè)計變量;gi(x)和n分別為約束條件以及約束條件數(shù)量;x1、x2為設(shè)計變量的取值范圍。

    當(dāng)存在不確定因素時,結(jié)構(gòu)穩(wěn)健優(yōu)化不但要考慮不確定因素對優(yōu)化模型的影響,還要考慮目標(biāo)函數(shù)的穩(wěn)健性,即能抵御不確定因素的影響。結(jié)合上一節(jié)的分析,注意到穩(wěn)健函數(shù)值大的結(jié)構(gòu)更為強健。基于優(yōu)化理論,建立穩(wěn)健函數(shù)最大化問題的優(yōu)化模型??紤]結(jié)構(gòu)體積(質(zhì)量)約束條件,橫截面向量a應(yīng)該滿足:

    式中:li>0 表示第i個桿長,為體積上限。對限定的位移uc和,應(yīng)力約束下魯棒函數(shù)最大化問題優(yōu)化模型可以寫為

    優(yōu)化目標(biāo)為穩(wěn)健性函數(shù)最優(yōu),即結(jié)構(gòu)抵抗不確定性載荷最大化。

    2.2 半定規(guī)劃模型

    上述過程中,基于I-G模型處理不確定載荷,則載荷集中包含無窮多載荷,常規(guī)方法很難對上式進行求解。因此通過半定規(guī)劃對穩(wěn)健優(yōu)化模型進行重新構(gòu)造。傳統(tǒng)半定規(guī)劃的數(shù)學(xué)形式為[14]

    式中 c,x∈Rn,F(xiàn)i(i=0,…,n)為n階常對稱矩陣,x 為設(shè)計變量,矩陣F(x)?0為正定。

    結(jié)合本文優(yōu)化模型式(13),令F?Rn表示滿足上述應(yīng)力約束(4)的力集合:

    進行如下定義:

    其中,p,t∈R,a∈Rm。將新參數(shù)(t,p)引入穩(wěn)健函數(shù)(a,uc)中,則可以得到:

    半定規(guī)劃算法有成熟的商業(yè)計算軟件包,可直接方便快捷調(diào)用,避免如牛頓法那樣推導(dǎo)雅可比矩陣和海森堡矩陣。

    3 數(shù)值仿真

    3.1 三桿系穩(wěn)健性優(yōu)化

    考慮如圖1所示的三桿件架體。其形式為對稱結(jié)構(gòu),桿件個數(shù)和自由度分別為m=3,n=2。為了結(jié)果分析方便,這里將單位歸一化處理。桿件3的長度l3=15且垂直于水平方向,θ=45°。桿件所受力狀態(tài)如圖1所示,名義載荷大小為100。考慮應(yīng)力約束式(3)=3.0(i=1,2,3)。設(shè)定結(jié)構(gòu)初始變量為a0=(20,20,20),以a=a0時結(jié)構(gòu)體積為基準(zhǔn),則結(jié)構(gòu)體積約束為≤1 148.53。

    圖1 三桿架體Fig.1 3-bar structure

    已知初始變量值a0,對式(18)進行求解,得到(a0,σc)=0.822。通過半定規(guī)劃模型式(19),對架結(jié)構(gòu)進行優(yōu)化設(shè)計,得到優(yōu)化迭代過程如表1。

    表1 優(yōu)化迭代過程Table 1 Iteration history of the optimization

    對結(jié)構(gòu)初始狀態(tài)a0進行分析。隨機產(chǎn)生一組力f,使其落在穩(wěn)健函數(shù)(a0,σc)限定的不確定載荷集之中。則結(jié)構(gòu)中桿件應(yīng)力σi(i=1,2,3)分布狀態(tài)如圖2所示。

    經(jīng)過觀察可以發(fā)現(xiàn),桿件3在應(yīng)力σ3=3.0附近有分布,而桿件1(或桿件2)應(yīng)力值始終沒有發(fā)揮最大功效。

    剛架經(jīng)過優(yōu)化后,最優(yōu)解為a=a*,經(jīng)過計算和分析,得到結(jié)構(gòu)桿件應(yīng)力σi(i=1,2,3)分布狀態(tài)如圖3所示。其中,應(yīng)力約束σi≤σc(i=1,2,3)在3個桿件中均產(chǎn)生作用,優(yōu)化后改善了結(jié)構(gòu)性能,最大程度發(fā)揮了剛架體結(jié)構(gòu)的潛力。

    圖2 當(dāng)a=a0時桿件應(yīng)力分布Fig.2 Stress states of the 3-barstructure with a=a0

    圖3 當(dāng)a=a*時桿件應(yīng)力分布Fig.3 Stress states of the 3-barstructure with a=a*

    3.2 多桿系剛架體結(jié)構(gòu)

    如圖4所示多桿架體結(jié)構(gòu),桿件個數(shù)和自由度分別為m=100,n=90。桿件在x軸和z軸方向的長度均為100,在y軸方向有桿長100和30(即中部連接處短桿)2種規(guī)格,其余為連接斜桿。結(jié)構(gòu)底座完全約束。載荷作用于如圖4所示位置,垂直方向載荷名義值為=(0,10,0);水平方向名義值=(2,0,0)。

    圖4 多桿架結(jié)構(gòu)Fig.4 Complex frame structure

    圖5 架體優(yōu)化設(shè)計Fig.5 Optimal design of the frame structure

    為了驗證優(yōu)化結(jié)果,觀察后標(biāo)記如圖6所示,對部分桿件進行分析。當(dāng)結(jié)構(gòu)中a=a0和a=a*時,產(chǎn)生一組隨機攝動載荷分別滿足應(yīng)力約束式(3)要求。

    圖6 架體平面結(jié)構(gòu)Fig.6 2-D structure of the frame structure

    圖7 隨機攝動載荷下最大應(yīng)力Fig.7 Fluctuation stress under random loading

    圖7表示優(yōu)化前后對比,即對應(yīng)a=a0、a*時,圖6中標(biāo)記對應(yīng)桿桿件最大應(yīng)力變化范圍。其中淺色為設(shè)計變量初始取值,深色為優(yōu)化后取值。通過觀察發(fā)現(xiàn),優(yōu)化后結(jié)構(gòu)在隨機攝動載荷下,最大應(yīng)力變化幅值總體上小于優(yōu)化前,結(jié)構(gòu)保持較大的安全富裕度,能更好抵抗不確定載荷的影響。

    4 結(jié)束語

    本文基于穩(wěn)健性函數(shù)構(gòu)建了剛架體結(jié)構(gòu)穩(wěn)健性優(yōu)化算法。該方法充分考慮了載荷的不確定因素,利用I-G理論處理不確定載荷,構(gòu)造應(yīng)力約束下的穩(wěn)健性函數(shù)最大值優(yōu)化模型,利用半定規(guī)劃松弛技術(shù)重新將優(yōu)化問題轉(zhuǎn)換成半定規(guī)劃問題,降低工程設(shè)計人員求解的難度,提高求解效率。通過數(shù)值仿真可以發(fā)現(xiàn)在優(yōu)化過程中約束條件發(fā)揮作用,在同等體積(質(zhì)量)條件下,雖結(jié)構(gòu)局部一些桿件應(yīng)力有所增加,但總體應(yīng)力狀態(tài)得到改善,提高了結(jié)構(gòu)抗擾動載荷的能力。

    [1]TAGUCHI G,CHOWDHURY S,TAGUCHI S.Robust engineering[M].New York:McGraw-Hill,2000:327-533.

    [2]張瑞軍,邱繼偉,賈慶軒,等.基于靈敏度附加項的多目標(biāo)可靠性穩(wěn)健優(yōu)化設(shè)計[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2013,21(4):777-785.ZHANG Ruijun,QIU Jiwei,JIA Qingxuan,et al.Multi-objective reliability-based robust optimization design based on the sensitivity additional items[J].Journal of Basic Science and Engineering,2013,21(4):777-785.

    [3]丁力平,馮毅雄,譚建榮,等.基于性能穩(wěn)健指數(shù)的產(chǎn)品族穩(wěn)健優(yōu)化設(shè)計[J].計算機集成制造系統(tǒng),2010,16(6):1121-1130.DING Liping,F(xiàn)ENG Yixiong,TAN Jianrong,et al.Robust optimization design of product family based on performance robustness index[J].Computer Integrated Manufacturing Systems,2010,16(6):1121-1130.

    [4]高一聰,馮毅雄,譚建榮,等.基于MSRE的機械產(chǎn)品質(zhì)量特性穩(wěn)健優(yōu)化設(shè)計方法[J].計算機集成制造系統(tǒng),2010,16(5):897-904.GAO Yicong,F(xiàn)ENG Yixiong,TAN Jianrong,et al.Product quality characteristics robust optimization design based on minimum sensitivity region estimation[J].Computer Integrated Manufacturing Systems,2010,16(5):897-904.

    [5]李小剛,程錦,劉振宇,等.基于雙層更新Kriging模型的機械結(jié)構(gòu)動態(tài)特性穩(wěn)健優(yōu)化設(shè)計[J].機械工程學(xué)報,2014,50(3):65-173.LI Xiaogang,CHENG Jin,LIU Zhenyu,et al.Robust optimization for dynamic characteristics of mechanical structures based on double renewal kriging model[J].Journal of Mechanical Engineering,2014,50(3):165-173.

    [6]于利磊,唐文勇,張圣坤,等.一種工程結(jié)構(gòu)的魯棒優(yōu)化設(shè)計方法[J].上海交通大學(xué)學(xué)報,2003,37(8):1189-1192.YU Lilei,TANG Wenyong,ZHANG Shengkun,et al.Robust optimization design method for engineering structures[J].Journal of Shanghai Jiaotong University,2003,37(8):1189-1192.

    [7]宋述芳,呂震宙.基于高階矩的穩(wěn)健優(yōu)化設(shè)計研究[J].力學(xué)學(xué)報,2012,44(4):735-744.SONG Shufang,LYU Zhenzhou.The robust optimization design based on moment estimation method[J].Chinese Journal of Theoretical and Applied Mechanics,2012,44(4):735-744.

    [8]BEN H Y,ELISHAKOFF I.Convex models of uncertainty in applied mechanics[M].Amsterdam:Elsevier,1990:175-230.

    [9]BEN-HAIM Y.Robust reliability of structures[J].Advances in Applied Mechanics,1997,33:1-41.

    [10]張成,吳慧,高博青,等.非概率不確定性結(jié)構(gòu)的魯棒性分析[J].計算力學(xué)學(xué)報,2013,30(1):51-56.ZHANG Cheng,WU Hui,GAO Boqing,et al.Analysis of structural robustness with non-probabilistic uncertainties[J].Chinese Journal of Computational Mechanics,2013,30(1):51-56.

    [11]王曉軍,王磊,賈曉,等.基于非概率凸模型可靠性的結(jié)構(gòu)優(yōu)化設(shè)計[J].北京航空航天大學(xué)學(xué)報,2012,38(5):630-635.WANG Xiaojun,WANG Lei,JIA Xiao,et al.Structural optimization design based on non-probabilistic convex modeling reliability[J].Journal of Beijing University of Aeronautics and Astronautics,2012,38(5)630-635.

    [12]劉曉,徐格寧,楊萍.區(qū)間干涉模型下的非概率可靠性應(yīng)用研究[J].中國礦業(yè)大學(xué)學(xué)報,2012,41(4):675-680.LIU Xiao,XU Gening,YANG Ping.Application of nonprobabilistic reliability based on interval interference model[J].Journal of China University of Minging & Technology,2012,41(4):675-680.

    [13]石海軍,錢林方,徐亞棟,等.具有參數(shù)區(qū)間不確定的機械臂結(jié)構(gòu)優(yōu)化[J].南京理工大學(xué)學(xué)報,2014,38(1):100-105.SHI Haijun,QIAN Linfang,XU Yadong,et al.Structural optimization for manipulator with interval uncertainties[J].Journal of Nanjing University of Science and Technology,2014,38(1):100-105.

    [14]WOLKOWICZ H,SAIGAL R,VANDENBERGHE L.Handbook of semidefinite programming:theory,algorithms,and applications[M].(s.l.):Kluwer Academic Publishers,2000:137-201.

    [15]YAKOV B H.Set-models of information-gap uncertainty:axioms and an inference scheme[J].Journal of the Franklin Institute,1999,336:1093-1117.

    [16]YAKOV B H.Uncertainty,probability and information-gaps[J].Reliability Engineering and System Safety,2004,85:249-266.

    [17]蘇永華,李翔.基于Info-Gap理論的地下結(jié)構(gòu)穩(wěn)健性分析方法[J].巖土工程學(xué)報,2011,33(2):227-233.SU Yonghua,LI Xiang.Robust reliability analysis for underground structures based on Info-Gap theory[J].Chinese Journal of Geotechnical Engineering,2011,33(2):227-233.

    [18]曹文貴,李媛,翟友成.基于Info-Gap理論的基樁下伏巖溶頂板穩(wěn)定性的主動分析方法[J].巖石力學(xué)與工程學(xué)報,2013,32(2):393-400.CAO Wengui,LI Yuan,ZHAI Youcheng.Active analysis method for stability of karst roof under foundation pile based on info-gap theory[J].Chinese Journal of Rock Mechanics and Engineering,2013,32(2):393-400.

    猜你喜歡
    剛架穩(wěn)健性桿件
    門式剛架結(jié)構(gòu)“借剛度”問題分析
    基于臨時支撐結(jié)構(gòu)的桿件初彎曲對其軸壓性能的影響
    四川建筑(2021年1期)2021-03-31 01:01:46
    塔式起重機拼裝式超長附著桿設(shè)計與應(yīng)用
    建筑機械化(2020年7期)2020-08-15 06:41:32
    會計穩(wěn)健性的定義和計量
    商情(2019年3期)2019-03-29 12:04:52
    會計穩(wěn)健性的文獻綜述
    財訊(2018年22期)2018-05-14 08:55:57
    剛架拱橋橫向整體性影響因素探討
    福建建筑(2018年3期)2018-03-29 01:14:13
    KD379:便攜折疊式衣架
    貨幣政策、會計穩(wěn)健性與銀行信貸關(guān)系探析
    某網(wǎng)架桿件彎曲的原因分析及處理
    平臺對門式剛架結(jié)構(gòu)穩(wěn)定性的影響分析
    netflix在线观看网站| 久久久精品国产亚洲av高清涩受| 狂野欧美激情性xxxx| 欧洲精品卡2卡3卡4卡5卡区| 国产激情欧美一区二区| 在线观看66精品国产| 亚洲av成人不卡在线观看播放网| 琪琪午夜伦伦电影理论片6080| 狠狠狠狠99中文字幕| 99久久国产精品久久久| 国产成人免费无遮挡视频| 精品少妇一区二区三区视频日本电影| 丁香欧美五月| 国产av精品麻豆| 久久精品成人免费网站| 欧美色视频一区免费| 一区二区三区高清视频在线| 亚洲国产欧美网| 欧美日本视频| 一边摸一边做爽爽视频免费| 欧美性长视频在线观看| 亚洲国产欧美网| 国产99久久九九免费精品| 国产麻豆69| 1024视频免费在线观看| 欧美精品亚洲一区二区| 性欧美人与动物交配| 久久久久久亚洲精品国产蜜桃av| bbb黄色大片| 青草久久国产| 色精品久久人妻99蜜桃| 这个男人来自地球电影免费观看| 亚洲国产精品成人综合色| 久久中文字幕人妻熟女| 搡老妇女老女人老熟妇| 亚洲欧美一区二区三区黑人| 国产精品 国内视频| 中亚洲国语对白在线视频| 日本 欧美在线| 亚洲黑人精品在线| 欧美成人性av电影在线观看| 精品一区二区三区四区五区乱码| 99精品在免费线老司机午夜| 91字幕亚洲| 亚洲av成人av| 一级作爱视频免费观看| 看片在线看免费视频| 亚洲五月色婷婷综合| av中文乱码字幕在线| 在线天堂中文资源库| 午夜福利视频1000在线观看 | 国产精品99久久99久久久不卡| 波多野结衣av一区二区av| 久久久国产精品麻豆| 真人一进一出gif抽搐免费| 级片在线观看| 国产99白浆流出| 国产欧美日韩一区二区三| 免费人成视频x8x8入口观看| 99久久精品国产亚洲精品| 窝窝影院91人妻| 亚洲中文字幕一区二区三区有码在线看 | 好男人电影高清在线观看| 亚洲男人天堂网一区| 国产精品久久视频播放| 欧美久久黑人一区二区| 国产成人欧美| 亚洲av成人av| 国产精品美女特级片免费视频播放器 | 国产精品综合久久久久久久免费 | 亚洲欧美精品综合久久99| 波多野结衣巨乳人妻| 成人av一区二区三区在线看| 欧美日本中文国产一区发布| 国产成人av教育| 一级毛片精品| 涩涩av久久男人的天堂| 最近最新中文字幕大全免费视频| 老司机午夜福利在线观看视频| 国产成+人综合+亚洲专区| 精品人妻1区二区| 亚洲人成77777在线视频| 美女 人体艺术 gogo| 亚洲 欧美一区二区三区| 在线观看66精品国产| 欧美一区二区精品小视频在线| 国产精品影院久久| 国产主播在线观看一区二区| 成人免费观看视频高清| 麻豆av在线久日| 操美女的视频在线观看| 欧美大码av| 国产99久久九九免费精品| 免费在线观看视频国产中文字幕亚洲| 亚洲一区中文字幕在线| 不卡av一区二区三区| 色哟哟哟哟哟哟| 午夜福利一区二区在线看| 90打野战视频偷拍视频| 男女下面插进去视频免费观看| 久久精品91无色码中文字幕| 18禁裸乳无遮挡免费网站照片 | 一个人免费在线观看的高清视频| 亚洲av第一区精品v没综合| 色精品久久人妻99蜜桃| 免费在线观看黄色视频的| 国产亚洲欧美98| 午夜成年电影在线免费观看| 精品一区二区三区av网在线观看| 精品高清国产在线一区| 一级a爱视频在线免费观看| 在线观看www视频免费| av在线播放免费不卡| 国产精品一区二区三区四区久久 | 人妻丰满熟妇av一区二区三区| 亚洲色图综合在线观看| 亚洲成a人片在线一区二区| 性少妇av在线| 91国产中文字幕| 免费在线观看亚洲国产| 亚洲avbb在线观看| 黑人欧美特级aaaaaa片| 国产在线精品亚洲第一网站| 99国产精品一区二区三区| 亚洲成人久久性| 国产主播在线观看一区二区| 色综合站精品国产| 亚洲av第一区精品v没综合| 1024视频免费在线观看| 精品熟女少妇八av免费久了| 嫩草影视91久久| 久久久久久人人人人人| 男女午夜视频在线观看| 亚洲视频免费观看视频| 成人亚洲精品一区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 老熟妇仑乱视频hdxx| 国产xxxxx性猛交| 变态另类丝袜制服| 午夜激情福利司机影院| 特级一级黄色大片| 日韩精品有码人妻一区| 在现免费观看毛片| 麻豆一二三区av精品| 人人妻,人人澡人人爽秒播| 我要看日韩黄色一级片| 成人鲁丝片一二三区免费| 亚洲,欧美,日韩| 精华霜和精华液先用哪个| 久久久久久久精品吃奶| 国内精品久久久久精免费| 欧美一区二区精品小视频在线| 久久久久久久午夜电影| 欧美中文日本在线观看视频| 亚洲中文字幕日韩| АⅤ资源中文在线天堂| 91精品国产九色| 国产白丝娇喘喷水9色精品| 免费av毛片视频| 国内精品宾馆在线| 欧美激情国产日韩精品一区| 九九久久精品国产亚洲av麻豆| 午夜免费激情av| av在线老鸭窝| 久久人人精品亚洲av| 最近中文字幕高清免费大全6 | 免费看a级黄色片| 一夜夜www| 蜜桃久久精品国产亚洲av| 小说图片视频综合网站| 亚洲狠狠婷婷综合久久图片| 婷婷丁香在线五月| 亚洲成人久久性| 听说在线观看完整版免费高清| 免费看av在线观看网站| 国产精品久久久久久久电影| 国产一区二区激情短视频| 两人在一起打扑克的视频| 亚洲欧美日韩东京热| 久久这里只有精品中国| 国产人妻一区二区三区在| 亚洲成av人片在线播放无| 美女免费视频网站| 美女cb高潮喷水在线观看| a级一级毛片免费在线观看| 午夜福利在线在线| 亚洲精品久久国产高清桃花| 18+在线观看网站| www日本黄色视频网| 性色avwww在线观看| 免费人成在线观看视频色| av在线观看视频网站免费| 女人十人毛片免费观看3o分钟| 别揉我奶头 嗯啊视频| 少妇高潮的动态图| 成人三级黄色视频| 国国产精品蜜臀av免费| 性插视频无遮挡在线免费观看| 午夜爱爱视频在线播放| 成人二区视频| 亚洲av中文字字幕乱码综合| 国产精品久久久久久久电影| 午夜免费成人在线视频| 国产精品精品国产色婷婷| 麻豆av噜噜一区二区三区| 最近最新中文字幕大全电影3| 91狼人影院| 日韩欧美 国产精品| 丰满的人妻完整版| 亚洲无线在线观看| 国产亚洲av嫩草精品影院| 亚洲欧美精品综合久久99| 中文字幕高清在线视频| 一区二区三区高清视频在线| 久久精品国产清高在天天线| 国产精品一区www在线观看 | 99久久精品国产国产毛片| 18+在线观看网站| 精品国产三级普通话版| 午夜激情福利司机影院| 亚洲专区国产一区二区| 我要看日韩黄色一级片| 久久久精品大字幕| 日韩精品有码人妻一区| 在线看三级毛片| 精品福利观看| 最近视频中文字幕2019在线8| 成人鲁丝片一二三区免费| 99视频精品全部免费 在线| 九色成人免费人妻av| 中出人妻视频一区二区| 亚洲电影在线观看av| 嫩草影院精品99| 人妻久久中文字幕网| 国产精品一区二区性色av| 久久久久久久久久黄片| 久99久视频精品免费| 日本熟妇午夜| 久久久久久国产a免费观看| 特大巨黑吊av在线直播| 日本黄大片高清| 男女之事视频高清在线观看| 日韩一区二区视频免费看| 日韩精品有码人妻一区| 免费在线观看成人毛片| 午夜福利在线观看免费完整高清在 | 亚洲电影在线观看av| 深爱激情五月婷婷| 免费在线观看日本一区| 久久精品国产自在天天线| 99精品久久久久人妻精品| 伊人久久精品亚洲午夜| 噜噜噜噜噜久久久久久91| 亚洲天堂国产精品一区在线| 午夜福利欧美成人| aaaaa片日本免费| 欧美最新免费一区二区三区| 日本免费a在线| 又黄又爽又刺激的免费视频.| 无遮挡黄片免费观看| 国产亚洲精品久久久com| 两性午夜刺激爽爽歪歪视频在线观看| 麻豆国产97在线/欧美| 中文字幕免费在线视频6| 91久久精品电影网| www.www免费av| 欧美人与善性xxx| 深夜精品福利| 直男gayav资源| 神马国产精品三级电影在线观看| 黄色女人牲交| 中文字幕精品亚洲无线码一区| 国产乱人伦免费视频| or卡值多少钱| 亚洲在线自拍视频| 免费av毛片视频| 女人被狂操c到高潮| 国内精品美女久久久久久| 女人十人毛片免费观看3o分钟| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 中文字幕免费在线视频6| 国产三级中文精品| 成人性生交大片免费视频hd| 亚洲中文字幕日韩| 久久午夜福利片| 国产精品亚洲美女久久久| 久久这里只有精品中国| 欧美黑人巨大hd| 又爽又黄a免费视频| 婷婷亚洲欧美| 3wmmmm亚洲av在线观看| 亚洲黑人精品在线| 亚洲男人的天堂狠狠| 看十八女毛片水多多多| 少妇丰满av| videossex国产| 日韩一区二区视频免费看| 午夜福利在线观看免费完整高清在 | 男人和女人高潮做爰伦理| 成人一区二区视频在线观看| 在线观看免费视频日本深夜| 国产精品av视频在线免费观看| 午夜福利视频1000在线观看| 成人欧美大片| 人妻夜夜爽99麻豆av| 久久精品影院6| 午夜精品在线福利| 99热这里只有是精品50| 全区人妻精品视频| 久久人人爽人人爽人人片va| 日韩 亚洲 欧美在线| 精品国内亚洲2022精品成人| 亚洲 国产 在线| 狂野欧美激情性xxxx在线观看| 好男人在线观看高清免费视频| 一本精品99久久精品77| 日韩中文字幕欧美一区二区| 哪里可以看免费的av片| 国产视频内射| 老熟妇乱子伦视频在线观看| 国产高清视频在线观看网站| 女人被狂操c到高潮| 国产白丝娇喘喷水9色精品| 国产精品一区二区三区四区久久| 国产麻豆成人av免费视频| 特级一级黄色大片| 国产91精品成人一区二区三区| 久99久视频精品免费| 在线观看舔阴道视频| 欧美bdsm另类| 禁无遮挡网站| 免费一级毛片在线播放高清视频| 国产 一区精品| 欧美黑人巨大hd| 美女被艹到高潮喷水动态| 亚洲精品在线观看二区| 观看美女的网站| 直男gayav资源| 一进一出抽搐gif免费好疼| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线播| 永久网站在线| 欧美一区二区亚洲| 级片在线观看| 欧美又色又爽又黄视频| 日本a在线网址| xxxwww97欧美| 极品教师在线免费播放| 亚洲天堂国产精品一区在线| 两个人的视频大全免费| 国产黄a三级三级三级人| 极品教师在线视频| 一区二区三区高清视频在线| 欧美zozozo另类| 美女被艹到高潮喷水动态| 香蕉av资源在线| 国产69精品久久久久777片| 成年女人看的毛片在线观看| 亚洲欧美日韩卡通动漫| 国产探花在线观看一区二区| 黄色女人牲交| 国产69精品久久久久777片| 三级男女做爰猛烈吃奶摸视频| 美女xxoo啪啪120秒动态图| 精品99又大又爽又粗少妇毛片 | 国内毛片毛片毛片毛片毛片| 真实男女啪啪啪动态图| 欧美日韩精品成人综合77777| av女优亚洲男人天堂| 波野结衣二区三区在线| 露出奶头的视频| 国产 一区 欧美 日韩| 久久国产乱子免费精品| 国产美女午夜福利| 精品福利观看| 国产一区二区三区在线臀色熟女| 别揉我奶头 嗯啊视频| 日韩欧美国产一区二区入口| 日本 欧美在线| 欧美高清成人免费视频www| 久久久久久久久久黄片| 两人在一起打扑克的视频| av.在线天堂| 中文字幕熟女人妻在线| 免费黄网站久久成人精品| 精品一区二区三区人妻视频| 少妇裸体淫交视频免费看高清| 欧美日韩亚洲国产一区二区在线观看| 看十八女毛片水多多多| 亚洲av日韩精品久久久久久密| 欧美国产日韩亚洲一区| 在线观看av片永久免费下载| 性插视频无遮挡在线免费观看| 免费一级毛片在线播放高清视频| 亚洲精品日韩av片在线观看| 99久久精品国产国产毛片| 亚洲,欧美,日韩| 日本精品一区二区三区蜜桃| 成人二区视频| 99久久久亚洲精品蜜臀av| 国产高清有码在线观看视频| 免费观看精品视频网站| 在线观看一区二区三区| 午夜爱爱视频在线播放| h日本视频在线播放| 嫩草影院精品99| 国产在线精品亚洲第一网站| 美女大奶头视频| 午夜精品久久久久久毛片777| 国产毛片a区久久久久| 久久久久久国产a免费观看| 日韩欧美国产在线观看| 深夜a级毛片| 别揉我奶头~嗯~啊~动态视频| 国产乱人伦免费视频| 国产精品自产拍在线观看55亚洲| 欧美高清性xxxxhd video| 国产爱豆传媒在线观看| 国产精品久久久久久av不卡| 色综合婷婷激情| 亚洲精品成人久久久久久| 一个人观看的视频www高清免费观看| 天堂av国产一区二区熟女人妻| 色综合亚洲欧美另类图片| 麻豆av噜噜一区二区三区| 精品一区二区三区人妻视频| 99久久无色码亚洲精品果冻| 国产精品不卡视频一区二区| 亚洲18禁久久av| 九九久久精品国产亚洲av麻豆| 一区二区三区四区激情视频 | 国产精品日韩av在线免费观看| 亚洲男人的天堂狠狠| 一级a爱片免费观看的视频| 99久国产av精品| 最新中文字幕久久久久| 久久人妻av系列| 长腿黑丝高跟| 国产毛片a区久久久久| 男女视频在线观看网站免费| 久久香蕉精品热| 久久午夜亚洲精品久久| 又粗又爽又猛毛片免费看| 欧美色欧美亚洲另类二区| 亚洲自拍偷在线| 老司机午夜福利在线观看视频| avwww免费| 老师上课跳d突然被开到最大视频| 天堂动漫精品| 在线国产一区二区在线| 国产高清有码在线观看视频| 亚洲av成人精品一区久久| 亚洲成人久久爱视频| 久久久久九九精品影院| 欧美日本视频| 国产亚洲91精品色在线| 日韩欧美国产在线观看| 无遮挡黄片免费观看| 男女之事视频高清在线观看| av在线老鸭窝| 九九久久精品国产亚洲av麻豆| 久9热在线精品视频| 日韩国内少妇激情av| 亚洲精品成人久久久久久| 嫁个100分男人电影在线观看| 亚洲无线观看免费| 亚洲国产精品sss在线观看| 久久久久九九精品影院| 欧美绝顶高潮抽搐喷水| 成人特级av手机在线观看| 久久久久久久亚洲中文字幕| 国产av麻豆久久久久久久| 欧美绝顶高潮抽搐喷水| 亚洲18禁久久av| 校园春色视频在线观看| 婷婷精品国产亚洲av在线| 国产精华一区二区三区| 网址你懂的国产日韩在线| 欧美性猛交╳xxx乱大交人| 日韩一区二区视频免费看| 亚洲av美国av| 人妻久久中文字幕网| 亚洲乱码一区二区免费版| 中文字幕精品亚洲无线码一区| 我的老师免费观看完整版| 嫁个100分男人电影在线观看| 亚洲专区国产一区二区| www.色视频.com| 亚洲精品影视一区二区三区av| 在线国产一区二区在线| 噜噜噜噜噜久久久久久91| 久久99热6这里只有精品| 色在线成人网| 亚洲黑人精品在线| 尾随美女入室| 舔av片在线| 如何舔出高潮| 午夜免费成人在线视频| 亚洲欧美清纯卡通| 小说图片视频综合网站| 一个人观看的视频www高清免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 我的女老师完整版在线观看| 高清在线国产一区| 久久精品夜夜夜夜夜久久蜜豆| xxxwww97欧美| 欧美又色又爽又黄视频| 99精品久久久久人妻精品| 国产精品久久久久久久电影| 午夜福利在线观看吧| 欧美不卡视频在线免费观看| 狠狠狠狠99中文字幕| 免费看光身美女| 观看美女的网站| 国产精品亚洲美女久久久| 深爱激情五月婷婷| 小说图片视频综合网站| 三级毛片av免费| 少妇丰满av| 人妻丰满熟妇av一区二区三区| 舔av片在线| .国产精品久久| 可以在线观看的亚洲视频| 97人妻精品一区二区三区麻豆| 国产在视频线在精品| 丰满乱子伦码专区| 白带黄色成豆腐渣| 麻豆国产av国片精品| 无遮挡黄片免费观看| 最近最新中文字幕大全电影3| 国产一区二区三区在线臀色熟女| 亚洲av二区三区四区| 国产人妻一区二区三区在| 日韩欧美精品v在线| 老女人水多毛片| 丰满乱子伦码专区| 国产精品1区2区在线观看.| 能在线免费观看的黄片| 男女下面进入的视频免费午夜| 人妻丰满熟妇av一区二区三区| 观看美女的网站| 日本-黄色视频高清免费观看| 最新在线观看一区二区三区| 美女黄网站色视频| 91在线观看av| 亚洲va在线va天堂va国产| 国产精品久久久久久亚洲av鲁大| 欧美高清性xxxxhd video| 有码 亚洲区| 少妇高潮的动态图| 亚洲成人久久爱视频| 男女下面进入的视频免费午夜| 日日夜夜操网爽| 久久6这里有精品| 国产精品一区二区性色av| 国产精品久久久久久久久免| 亚洲熟妇中文字幕五十中出| 少妇熟女aⅴ在线视频| 久久午夜亚洲精品久久| 美女免费视频网站| 大又大粗又爽又黄少妇毛片口| 岛国在线免费视频观看| 国产精品综合久久久久久久免费| 精品99又大又爽又粗少妇毛片 | 亚洲图色成人| 亚洲国产欧美人成| 国产精品98久久久久久宅男小说| 少妇人妻精品综合一区二区 | 国产精品99久久久久久久久| 午夜老司机福利剧场| 日本黄色视频三级网站网址| 久9热在线精品视频| 亚洲精品亚洲一区二区| 国产单亲对白刺激| 亚洲美女视频黄频| 国产精品一区二区免费欧美| 少妇的逼好多水| 亚洲久久久久久中文字幕| 亚洲精品在线观看二区| 亚洲av免费高清在线观看| 黄色丝袜av网址大全| 精品午夜福利视频在线观看一区| 校园春色视频在线观看| 国产伦在线观看视频一区| 午夜免费成人在线视频| 99热只有精品国产| 又爽又黄无遮挡网站| 1000部很黄的大片| 成人午夜高清在线视频| 欧美日韩亚洲国产一区二区在线观看| 999久久久精品免费观看国产| 欧美成人性av电影在线观看| 国产探花在线观看一区二区| 亚洲美女视频黄频| 亚洲精品456在线播放app | 国产av在哪里看| 真人一进一出gif抽搐免费| 欧美最新免费一区二区三区| 深爱激情五月婷婷| 国产精品野战在线观看| 我的女老师完整版在线观看| 国产精品久久电影中文字幕| 国产精品国产三级国产av玫瑰| 国产av麻豆久久久久久久| 深爱激情五月婷婷| 亚洲经典国产精华液单| 最近在线观看免费完整版| 国产一区二区在线观看日韩| 午夜视频国产福利| 高清毛片免费观看视频网站| 一级a爱片免费观看的视频| a级一级毛片免费在线观看| 在线观看午夜福利视频| 人人妻人人澡欧美一区二区|