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

    智能桁架結(jié)構(gòu)最優(yōu)振動(dòng)控制與作動(dòng)器優(yōu)化配置

    2016-01-12 10:39:11曹玉巖,王志臣,付世欣
    振動(dòng)與沖擊 2015年5期
    關(guān)鍵詞:振動(dòng)控制模擬退火

    第一作者曹玉巖男,碩士,助理研究員,1986年2月生

    智能桁架結(jié)構(gòu)最優(yōu)振動(dòng)控制與作動(dòng)器優(yōu)化配置

    曹玉巖,王志臣,付世欣,吳慶林

    (中國(guó)科學(xué)院長(zhǎng)春光學(xué)精密機(jī)械與物理研究所,吉林130033)

    摘要:研究了智能桁架結(jié)構(gòu)最優(yōu)振動(dòng)控制和作動(dòng)器的優(yōu)化配置問(wèn)題。首先采用有限元方法,根據(jù)Hamilton原理推導(dǎo)了智能桁架結(jié)構(gòu)的機(jī)電耦合動(dòng)力學(xué)方程,根據(jù)線性二次型最優(yōu)控制理論,推導(dǎo)了結(jié)構(gòu)振動(dòng)控制的數(shù)學(xué)模型,通過(guò)最小化性能泛函,求解黎卡提矩陣代數(shù)方程確定了最優(yōu)控制輸入。然后通過(guò)對(duì)最優(yōu)控制性能指標(biāo)函數(shù)的修正,得到了與初始狀態(tài)無(wú)關(guān)的性能指標(biāo),以修正的性能指標(biāo)為目標(biāo)函數(shù),應(yīng)用模擬退火算法對(duì)作動(dòng)器位置進(jìn)行了優(yōu)化配置。最后給出了空間智能桁架結(jié)構(gòu)振動(dòng)控制算例驗(yàn)證建模過(guò)程和算法。算例結(jié)果表明,通過(guò)最優(yōu)振動(dòng)控制可以使結(jié)構(gòu)振動(dòng)快速衰減,達(dá)到振動(dòng)抑制的效果,而且通過(guò)模擬退火算法可以確定最佳的作動(dòng)器布置方式。

    關(guān)鍵詞:智能桁架;Hamilton原理;振動(dòng)控制;模擬退火

    收稿日期:2013-11-08修改稿收到日期:2014-02-12

    中圖分類號(hào):TP13文獻(xiàn)標(biāo)志碼:A

    Optimal vibration control and actuator placement for intelligent truss structures

    CAOYu-yan,WANGZhi-chen,FUShi-xin,WUQing-lin(Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences,Changchun 130033, China)

    Abstract:The optimal vibration control and actuator placement for intelligent truss structures was investigated. The mech-electric coupling dynamic formulation of piezoelectric intelligent truss structures based on Hamilton principle was presented. The mathematical model of vibration control for truss structures was derived using linear quadratic Gauss optimal control theory. The input of optimal control was determined by minimizing the performance function and solving the Riccati algebraic equation. Then, the performance criterion unrelated to initial state was obtained through correcting the performance criterion of optimal control. The optimization of actuator position was performed using simulated annealing algorithm so as to minimize the corrected performance criterion. The numerical example of vibration control for a spatial intelligent truss structure was given to verify the modeling process and optimization algorithm. The simulations results indicate that the vibration suppression effect is obvious and the vibration responses decrease quickly by the optimal vibration control.

    Key words:intelligent truss; Hamilton theory; vibration control; simulated annealing

    隨著航天事業(yè)的發(fā)展,航天器朝著輕量化、大型化方向發(fā)展,桁架結(jié)構(gòu)作為一種特殊的結(jié)構(gòu)得到了越來(lái)越普遍的應(yīng)用[1-2]。由于其自身低剛度、弱阻尼的特點(diǎn),振動(dòng)問(wèn)題很突出,在很長(zhǎng)一段時(shí)間內(nèi),桁架結(jié)構(gòu)振動(dòng)控制問(wèn)題一直是研究的熱點(diǎn),尤其是壓電陶瓷堆作動(dòng)器的出現(xiàn),使得復(fù)雜桁架結(jié)構(gòu)主動(dòng)振動(dòng)控制成為可能[3]。將作動(dòng)器和傳感器與桁架結(jié)構(gòu)結(jié)合一起形成自適應(yīng)結(jié)構(gòu),在振動(dòng)控制、噪聲抑制、形狀控制和精密定位領(lǐng)域具有廣泛的應(yīng)用前景。

    在國(guó)內(nèi)外,智能結(jié)構(gòu)的振動(dòng)主動(dòng)控制及優(yōu)化方法已經(jīng)取得了一定的進(jìn)展[1,4-6],結(jié)構(gòu)振動(dòng)控制和作動(dòng)器的最優(yōu)布置一直是該領(lǐng)域研究的主要問(wèn)題。Moita等[1]研究了壓電層合板結(jié)構(gòu)的振動(dòng)控制問(wèn)題,并采用模擬退火算法對(duì)貼片作動(dòng)器的位置進(jìn)行了優(yōu)化配置。Yan等[7]研究了空間桁架結(jié)構(gòu)的振動(dòng)控制問(wèn)題,作動(dòng)器的優(yōu)化配置采用遺傳算法實(shí)現(xiàn),但其優(yōu)化結(jié)果受到初始條件的影響很大,不同的初始條件將得到不同的優(yōu)化配置結(jié)果。徐亞蘭等[8]基于結(jié)構(gòu)模態(tài)空間范數(shù)研究了壓電柔性結(jié)構(gòu)振動(dòng)控制問(wèn)題。謝永等[9]對(duì)10跨桁架結(jié)構(gòu)的動(dòng)力特性和主動(dòng)振動(dòng)控制進(jìn)行實(shí)驗(yàn)研究。曾光等[10]研究了基于遺傳算法的空間智能桁架結(jié)構(gòu)振動(dòng)控制中的作動(dòng)器/傳感器位置優(yōu)化問(wèn)題。徐斌等[11]針對(duì)區(qū)間參數(shù)壓電智能結(jié)構(gòu)/控制一體化多目標(biāo)拓?fù)鋬?yōu)化問(wèn)題,提出了一種區(qū)間參數(shù)結(jié)構(gòu)/控制多目標(biāo)拓?fù)鋬?yōu)化方法。彭海軍等[12]針對(duì)桁架結(jié)構(gòu)的振動(dòng)抑制問(wèn)題,采用基于變形修正的物理降階方法對(duì)有限元模型降階,并設(shè)計(jì)線性二次高斯(LQG)最優(yōu)控制律。趙國(guó)偉等[13]對(duì)柔性自適應(yīng)桁架最優(yōu)振動(dòng)控制進(jìn)行了實(shí)驗(yàn)研究。

    本文研究了智能桁架結(jié)構(gòu)最優(yōu)振動(dòng)控制和作動(dòng)器的優(yōu)化配置問(wèn)題。首先采用有限元方法和最優(yōu)控制理論,建立智能桁架結(jié)構(gòu)最優(yōu)振動(dòng)控制模型,通過(guò)最小化性能泛函,并求解黎卡提(Riccati)矩陣代數(shù)方程確定了最優(yōu)控制輸入。然后應(yīng)用模擬退火算法對(duì)作動(dòng)器的位置進(jìn)行優(yōu)化。最后給出空間智能桁架結(jié)構(gòu)的仿真算例,驗(yàn)證建模過(guò)程和算法。

    1智能桁架結(jié)構(gòu)最優(yōu)振動(dòng)控制系統(tǒng)建模

    1.1智能桁架結(jié)構(gòu)動(dòng)力學(xué)建模

    智能桁架結(jié)構(gòu)的桿件包括壓電作動(dòng)器和主體桿件,其中壓電作動(dòng)器是由多層壓電片堆疊膠合而成的,李東旭等[14-16]詳細(xì)介紹了壓電材料單元的力學(xué)方程,這里僅作簡(jiǎn)要敘述。作動(dòng)器結(jié)構(gòu)如圖1所示,在電壓的作用下可產(chǎn)生軸向伸縮位移,其動(dòng)力學(xué)模型如圖2所示。

    圖1 壓電陶瓷堆作動(dòng)器結(jié)構(gòu) Fig.1 Structure of piezoelectric actuator

    圖2 壓電作動(dòng)器力學(xué)模型簡(jiǎn)化 Fig.2 Mechanic model of piezoelectric actuator

    根據(jù)Hamilton變分原理,軸向伸縮PZT壓電元件機(jī)電耦合動(dòng)力學(xué)方程為

    (1)

    式中:L=T-U為L(zhǎng)agrange函數(shù),T為系統(tǒng)動(dòng)能,U為系統(tǒng)勢(shì)能,δ′W為外載荷所做虛功(外載荷所做功的變分)。

    由式(1)可得壓電作動(dòng)器單元?jiǎng)恿W(xué)方程為[16]

    (2)

    (3)

    (4)

    壓電主動(dòng)桿是由壓電作動(dòng)器兩端與主體桿件連接而成的可軸向伸縮的桿件,結(jié)構(gòu)如圖3所示,其中單元1、3為主體材料單元,單元2為壓電作動(dòng)器單元。

    圖3 主動(dòng)桿示意圖 Fig.3 Diagram of active members

    單元1、3的有限元方程為

    (5)

    式中,上標(biāo)(h)表示主體桿件材料的相關(guān)矩陣,其含義與式(2)一致。

    對(duì)單元1或3,其方程的具體形式為

    (6)

    對(duì)單元1、2、3的質(zhì)量矩陣和剛度矩陣進(jìn)行組裝,得到總體坐標(biāo)系下結(jié)構(gòu)動(dòng)力學(xué)方程為

    (7)

    式中,上標(biāo)(c)表示主動(dòng)桿的相關(guān)矩陣。

    矩陣[H],[M]的表達(dá)式為

    式中,λ=k1k2+k2k3+k1k3,ki(i=1,2,3)為單元i的拉伸剛度,kuv為作動(dòng)器的耦合剛度。

    將式(8)代入式(7)得到主動(dòng)桿的機(jī)電耦合動(dòng)力學(xué)方程為

    [Kuv]e{δ}e+{Kvv}eφ=-{Q}e

    (10)

    通過(guò)縮聚變換,使壓電主動(dòng)桿的方程與主體桿單元具有相同形式。根據(jù)以上推導(dǎo)得到了壓電主動(dòng)桿件模型,通過(guò)坐標(biāo)變換將各單元的質(zhì)量矩陣和剛度矩陣變換為總體坐標(biāo)系下的特性矩陣,然后對(duì)其進(jìn)行組集,即得到壓電桁架結(jié)構(gòu)的無(wú)阻尼動(dòng)力學(xué)方程為

    (11)

    式中,

    1.2最優(yōu)振動(dòng)控制建模

    由于桁架結(jié)構(gòu)的特性矩陣如剛度矩陣、質(zhì)量矩陣等均是時(shí)不變的,結(jié)構(gòu)振動(dòng)控制系統(tǒng)屬于線性定常系統(tǒng),根據(jù)最優(yōu)控制理論,性能指標(biāo)是狀態(tài)變量的二次函數(shù),可以把這種動(dòng)態(tài)系統(tǒng)的最優(yōu)控制問(wèn)題視為線性二次型最優(yōu)控制問(wèn)題[17]。

    為了應(yīng)用線性二次型最優(yōu)控制,將式(11)變換為狀態(tài)空間方程表達(dá)式為

    (12)

    為了能在消耗最少控制能量的情況下使結(jié)構(gòu)振動(dòng)迅速衰減,根據(jù)線性二次型最優(yōu)控制理論,性能指標(biāo)函數(shù)表達(dá)為

    (13)

    式中,Q∈R2n×2n為半正定的狀態(tài)加權(quán)矩陣,R∈Rr×r為正定控制加權(quán)矩陣,Q,R為對(duì)陣矩陣且常取為對(duì)角陣。

    根據(jù)最優(yōu)控制理論,使性能指標(biāo)函數(shù)達(dá)到最小值的充分必要條件為

    u(t)=-Kx(t)=-R-1BTPx(t)

    (14)

    矩陣P可以通過(guò)求解黎卡提(Ricatii)矩陣代數(shù)方程得到,即

    ATP+PA+Q-PBR-1BTP=0

    (15)

    性能指標(biāo)函數(shù)的最小值為

    (16)

    將式(14)代入式(12)得

    (17)

    由式(17)可求得

    x=x0e(A-BR-1BTP)t

    (18)

    2壓電作動(dòng)器的優(yōu)化配置

    壓電作動(dòng)器的最優(yōu)配置問(wèn)題屬于離散規(guī)劃問(wèn)題,有很多優(yōu)化方法如遺傳算法、模擬退火算法可以求解該問(wèn)題,在這里采用模擬退火算法,與傳統(tǒng)的應(yīng)用梯度的優(yōu)化方法相比,模擬退火算法的優(yōu)點(diǎn)可以防止過(guò)早收斂與局部最優(yōu)解。

    2.1模擬退火算法

    模擬退火算法于1982年由Kirkpatrick等引入組合優(yōu)化領(lǐng)域,它源于對(duì)固體退火過(guò)程的模擬,采用Metropolis接受準(zhǔn)則,并適當(dāng)?shù)目刂茰囟鹊南陆颠^(guò)程實(shí)現(xiàn)模擬退火,從而達(dá)到求解全局優(yōu)化問(wèn)題的目的。模擬退火算法中按Metropolis接受準(zhǔn)則確定。

    算法中是否接受從當(dāng)前解x1到新解x2的轉(zhuǎn)移由Metropolis接受準(zhǔn)則決定,可以描述為

    2.2作動(dòng)器優(yōu)化配置

    根據(jù)線性二次型最優(yōu)控制理論,最小化性能指標(biāo)相當(dāng)于最小化系統(tǒng)能量與控制能量,將式(18)和式(14)代入式(13)得

    (20)

    從式(20)可以看出性能指標(biāo)函數(shù)與初始狀態(tài)有直接關(guān)系,為了消除初始狀態(tài)對(duì)優(yōu)化結(jié)果的影響,將對(duì)能量函數(shù)做如下修正[18],并將修正后的能量函數(shù)作為優(yōu)化目標(biāo)函數(shù)。假設(shè)系統(tǒng)的初始狀態(tài)x0在n維單位球面上隨機(jī)分布,則式(20)將變化為

    (21)

    式中,Aeq=A-BR-1BTP

    Qeq=Q+(R-1BTP)TR(R-1BTP)

    對(duì)式(21)計(jì)算關(guān)于初始狀態(tài)x0的數(shù)學(xué)期望為

    (22)

    (23)

    (24)

    修正后的系統(tǒng)性能指標(biāo)函數(shù)可以表達(dá)為

    (25)

    式(25)所表示的性能指標(biāo)的意義為:當(dāng)系統(tǒng)的初始狀態(tài)x0在n維單位球面上隨機(jī)分布時(shí),最小化修正性能指標(biāo)式(25)等價(jià)于最小化性能指標(biāo)式(25)關(guān)于初始狀態(tài)的數(shù)學(xué)期望。

    根據(jù)上述推導(dǎo),將修正的性能指標(biāo)函數(shù),即式(25)作為作動(dòng)器優(yōu)化配置的目標(biāo)函數(shù),可以建立作動(dòng)器優(yōu)化配置的數(shù)學(xué)模型為

    (26)

    式中,λj表示第j個(gè)桿件是否配置作動(dòng)器,0-不安裝,1-安裝,n,m分別表示桁架桿件總數(shù)和作動(dòng)器總數(shù)。

    2.3算法中初始溫度設(shè)置及新解構(gòu)造

    初始溫度T0及降溫準(zhǔn)則將對(duì)算法的性能有很大影響,必須設(shè)置足夠高,至少在初始時(shí)保證所有解都能被接受,這里采用多次隨機(jī)布置產(chǎn)生一組變量,計(jì)算相應(yīng)的目標(biāo)函數(shù)值,利用這組目標(biāo)函數(shù)的最大值與最小值之差來(lái)確定初始溫度狀態(tài),即

    (27)

    式中,pr為初始接受概率。

    (28)

    3仿真算例

    以如圖4所示20節(jié)點(diǎn)64桿空間智能桁架結(jié)構(gòu)為例,對(duì)其進(jìn)行相應(yīng)的仿真分析。桁架節(jié)點(diǎn)和單元編號(hào)如圖4所示,約束1~4號(hào)節(jié)點(diǎn)的全部自由度,桁架結(jié)構(gòu)共有16×3=48自由度。桁架結(jié)構(gòu)幾何參數(shù)和材料特性參數(shù)如表1所示。

    表1 桁架中桿件及作動(dòng)器相關(guān)參數(shù)

    圖4 空間桁架結(jié)構(gòu) Fig.4 Structure of truss

    仿真驗(yàn)證三種情況:首先在每個(gè)桿件上均布置作動(dòng)器,共64個(gè)作動(dòng)器,分析最優(yōu)控制前后結(jié)構(gòu)振動(dòng)衰減情況;然后通過(guò)分析布置10(1~64任取)個(gè)作動(dòng)器情況下優(yōu)化布置前后的控制效果,來(lái)驗(yàn)證優(yōu)化配置算法;最后分析布置不同數(shù)量作動(dòng)器情況,優(yōu)化性能指標(biāo)與作動(dòng)器數(shù)目的關(guān)系。其中最優(yōu)振動(dòng)控制加權(quán)矩陣選取Q=106I,R=I。

    3.1每個(gè)桿件均布置作動(dòng)器情況

    在每個(gè)桿件位置均布置作動(dòng)器的情況下,共有64

    個(gè)壓電作動(dòng)器,對(duì)三種不同的初始擾動(dòng)條件進(jìn)行了動(dòng)態(tài)響應(yīng)分析,并采用提出的方法進(jìn)行最優(yōu)振動(dòng)控制。三種不同的初始擾動(dòng)條件分別為:

    結(jié)構(gòu)自由振動(dòng)響應(yīng)和最優(yōu)振動(dòng)控制后的響應(yīng)如圖5所示。

    從圖5(a)~(c)中可以看出,通過(guò)最優(yōu)振動(dòng)控制可以使結(jié)構(gòu)振動(dòng)快速衰減,1s后衰減90%以上,1.5s后基本衰減為0。

    3.2配置10個(gè)壓電作動(dòng)器情況

    對(duì)個(gè)作動(dòng)器的位置進(jìn)行優(yōu)化配置,采用模擬退火算法對(duì)其進(jìn)行優(yōu)化配置,并對(duì)優(yōu)化配置前后最優(yōu)振動(dòng)控制的結(jié)果進(jìn)行對(duì)比。優(yōu)化求解參數(shù)如表2所示,優(yōu)化過(guò)程中,優(yōu)化前后個(gè)作動(dòng)器布置位置如圖6所示,目標(biāo)函數(shù)變化如圖7所示,新解接受與拒絕柱形圖如圖8所示,溫度下降曲線如圖9所示。

    圖5 布置64個(gè)壓電作動(dòng)器在三種工況下,智能桁架結(jié)構(gòu) 最優(yōu)振動(dòng)控制前后節(jié)點(diǎn)位移動(dòng)態(tài)響應(yīng),圖(a)(b)(c)分別對(duì)應(yīng)初始擾動(dòng)條件(Ⅰ)(Ⅱ)和(Ⅲ) Fig.5 Dynamic responses of the intelligent truss structure without control and with optimal control by actuators for different initial disturbances, figures (a), (b) and (c) show the response for initial disturbance case (Ⅰ), (Ⅱ) and (Ⅲ) respectively

    表2 優(yōu)化求解參數(shù)

    從圖10(a)和(b)可以看出,優(yōu)化配置前,20節(jié)點(diǎn)X向動(dòng)態(tài)位移衰減很慢,衰減時(shí)間較長(zhǎng),1.5 s時(shí)還存在較大幅度的振動(dòng)位移,而優(yōu)化配置后,動(dòng)態(tài)位移衰減很快,1 s內(nèi)位移得到了很好的抑制。

    圖6 優(yōu)化配置前后10個(gè)作動(dòng)器的位置 Fig.6 Initial and optimal locations for actuators

    3.3配置不同數(shù)量作動(dòng)器情況

    為了研究作動(dòng)器數(shù)目對(duì)控制性能的影響,對(duì)配置不同數(shù)目作動(dòng)器的情況進(jìn)行仿真分析,配置作動(dòng)器數(shù)目為5~62,性能指標(biāo)變化如圖11所示。

    從圖11可以看出,隨著作動(dòng)器數(shù)目的增加,性能指標(biāo)μ隨之減小,控制總能量不斷下降,即作動(dòng)器越多控制性能越好。

    從圖11還可以看出,在作動(dòng)器數(shù)目1~10階段,目標(biāo)函數(shù)下降速度較快,目標(biāo)函數(shù)值減小70%以上,而在作動(dòng)器數(shù)目17~62階段,目標(biāo)函數(shù)僅下降30%。由此可知,作動(dòng)器數(shù)目存在一個(gè)臨界值,在達(dá)到臨界數(shù)目前,作動(dòng)器數(shù)目的增加對(duì)控制性能影響較大,而超過(guò)臨界數(shù)目后,對(duì)控制性能影響不明顯。因此對(duì)桁架結(jié)構(gòu)的振動(dòng)控制,僅需配置很少一部分作動(dòng)器就能達(dá)到理想的控制效果。

    圖7 優(yōu)化配置10個(gè)作動(dòng)器目標(biāo)函數(shù)變化Fig.7Objectivefunctioncurveforoptimization圖8 優(yōu)化配置10個(gè)作動(dòng)器接受與拒絕圖Fig.8Acceptanceandrefusalcurveforoptimization圖9 優(yōu)化配置10個(gè)作動(dòng)器降溫曲線Fig.9Temperaturedescendcurveforoptimization

    圖10 優(yōu)化配置前后20節(jié)點(diǎn)X向動(dòng)態(tài)位移響應(yīng)Fig.10Thedynamicdisplacementresponseofnode2direction圖11 作動(dòng)器數(shù)目與性能指標(biāo)變化Fig.11Therelationsbetweenactuatornumbersandobjectionfunction

    4結(jié)論

    研究了智能桁架結(jié)構(gòu)最優(yōu)振動(dòng)控制與作動(dòng)器最優(yōu)配置問(wèn)題,首先采用有限元方法推導(dǎo)了智能桁架結(jié)構(gòu)振動(dòng)控制的數(shù)學(xué)模型,通過(guò)最小化性能泛函,求解黎卡提(Riccati)矩陣代數(shù)方程確定了最優(yōu)控制輸入。然后以修正的性能指標(biāo)為目標(biāo)函數(shù),應(yīng)用模擬退火算法對(duì)作動(dòng)器位置進(jìn)行了優(yōu)化配置。最后給出了桁架振動(dòng)控制算例驗(yàn)證建模過(guò)程和算法。

    算例結(jié)果表明,(1)通過(guò)最優(yōu)振動(dòng)控制可以使結(jié)構(gòu)振動(dòng)快速衰減,能夠達(dá)到很好的振動(dòng)抑制效果;(2)對(duì)于確定相同數(shù)量的作動(dòng)器,不同的布置方式控制效果差別很大,通過(guò)模擬退火算法優(yōu)化,可以確定出最佳布置方式,使控制性能得到明顯改善;(3)布置的作動(dòng)器數(shù)量越多,控制效果越好,但作動(dòng)器數(shù)目存在一個(gè)臨界值,小于該臨界值時(shí),作動(dòng)器數(shù)目的增加對(duì)控制性能的影響非常顯著,大于該臨界值時(shí),作動(dòng)器數(shù)目的增加對(duì)控制性能的影響不明顯。

    參考文獻(xiàn)

    [1]Moita J M S, Correia V M F, Martins P G, et al. Optimal design in vibration control of adaptive structures using a simulated annealing algorithm[J]. Composite Structures, 2006, 75(1):79-87.

    [2]Muc A, Kedziora P. Optimal design of smart laminated composite structures[J]. Materials And Manufacturing Processes, 2010, 25(4):272-280.

    [3]Irschik H. A review on static and dynamic shape control of structures by piezoelectric actuation[J]. Engineering Structures, 2002, 24(1):5-11.

    [4]Frecker M I. Recent advances in optimization of smart structures and actuators[J]. Journal of Intelligent Material Systems And Structures, 2003, 14(4-5):207-216.

    [5]Moita J S, Martins P G, Mota Soares C M, et al. Optimal dynamic control of laminated adaptive structures using a higher order model and a genetic algorithm[J]. Computers & Structures, 2008, 86(3-5):198-206.

    [6]聶潤(rùn)兔, 邵成勛, 鄒振祝. 智能桁架機(jī)電耦合動(dòng)力分析與振動(dòng)控制[J]. 振動(dòng)工程學(xué)報(bào),1997,10(2):3-8.

    NIE Run-tu, SHAO Cheng-xun, ZOU Zhen-zhu. Mech-electric coupling dynamic analysis and vibration control of intelligent truss structures[J]. Journal of Vibration Engineering, 1997, 10(2): 3-8.

    [7]Yan Y J, Yam L H. Optimal design of number and locations ofactuators in active vibration control of a space truss[J]. Smart Materials and Structures, 2002, 11(4):496-503.

    [8]徐亞蘭, 陳建軍. 空間范數(shù)在壓電柔性結(jié)構(gòu)振動(dòng)控制中的應(yīng)用研究[J]. 工程力學(xué),2013,30(10):1-4.

    XU Ya-lan,CHEN Jian-jun. Researcn on vibration control of piezoelectric flexible structures based on spatial norms[J]. Engineering Mechanics, 2013, 30(10): 1-4.

    [9]謝永, 趙童, 梁辰, 等. 桁架結(jié)構(gòu)振動(dòng)主動(dòng)控制的實(shí)驗(yàn)研究[J]. 力學(xué)季刊,2010,31(4):578-581.

    XIE Yong, ZHAO Tong, LIANG Chen, et al. Experimental Study of Active Control for a Trussed Structure[J]. Chinese Quarterly of Mechanics, 2010, 31(4): 578-581.

    [10]曾光, 李東旭. 空間智能桁架作動(dòng)器/傳感器位置優(yōu)化中的遺傳算法應(yīng)用[J]. 宇航學(xué)報(bào),2007,28(2):461-464.

    ZENG Guang, LI Dong-xu. Placement optimization of actuatorsPSensors of space intelligent truss using genetic algorithms[J]. Journal of Astronautics, 2007, 28(2): 461-464.

    [11]徐斌, 趙普猛, 黎瑩. 區(qū)間參數(shù)壓電智能桁架結(jié)構(gòu)/控制的多目標(biāo)非概率可靠性拓?fù)鋬?yōu)化[J]. 振動(dòng)工程學(xué)報(bào),2013,26(2):169-177.

    XU Bin, ZHAO Pu-meng, LI Ying. Multi-objective optimization for structural topology and control of piezoelectric smart truss with interval parameter considering non-probability reliability constraints[J].Journal of Vibration Engineering, 2013,26(2): 169-177.

    [12]彭海軍, 王文勝, 程耿東. 基于物理降階模型的桁架結(jié)構(gòu)振動(dòng)主動(dòng)控制[J]. 工程力學(xué),2013,30(12):1-7.

    PENG Hai-jun,WANG Wen-sheng, CHENG Geng-dong. Active vibration control of a trussed structure based on the physical reduced model[J]. Engineering Mechanics, 2013, 30(12): 1-7.

    [13]趙國(guó)偉, 黃海, 夏人偉. 柔性自適應(yīng)桁架及其振動(dòng)最優(yōu)控制實(shí)驗(yàn)[J]. 北京航空航天大學(xué)學(xué)報(bào),2005,31(4):434-438.

    ZHAO Guo-wei, HUANG Hai, XIA Ren-wei. Flexible adaptive truss testbed and its optimal control experiment for vibration[J]. Journal of Beijing University of Aeronautics and Astronautics, 2005, 31(4): 434-438.

    [14]李東旭. 大型撓性空間桁架結(jié)構(gòu)動(dòng)力學(xué)分析與模糊振動(dòng)控制[M]. 北京: 科學(xué)出版社, 2008.

    [15]司洪偉. 大撓性航天桁架結(jié)構(gòu)動(dòng)力學(xué)建模及其主動(dòng)模糊控制研究[D]. 長(zhǎng)沙: 國(guó)防科學(xué)技術(shù)大學(xué), 2006.

    [16]龍連春. 智能桁架結(jié)構(gòu)受力性態(tài)最優(yōu)控制的建模與分析[D]. 北京: 北京工業(yè)大學(xué), 2003.

    [17]劉豹, 唐萬(wàn)生. 現(xiàn)代控制理論[M]. 北京: 機(jī)械工業(yè)出版社, 2006.

    [18]Kleinman D L, Athans M. The Design of Suboptimal Linear Time-Varying Systems[J]. IEEE Transactions on Automatic Control, 1968, 13(2):150-159.

    猜你喜歡
    振動(dòng)控制模擬退火
    結(jié)合模擬退火和多分配策略的密度峰值聚類算法
    模擬退火遺傳算法在機(jī)械臂路徑規(guī)劃中的應(yīng)用
    某型發(fā)動(dòng)機(jī)傳動(dòng)裝配過(guò)程中的振動(dòng)控制技術(shù)研究
    阻尼器對(duì)懸索橋吊索扭轉(zhuǎn)振動(dòng)控制效果的數(shù)值研究
    剪切型轉(zhuǎn)動(dòng)粘彈性阻尼器在村鎮(zhèn)木結(jié)構(gòu)抗風(fēng)中的應(yīng)用
    應(yīng)急柴油機(jī)冷卻管道振動(dòng)控制優(yōu)化設(shè)計(jì)研究
    基于模擬退火剩余矩形算法的矩形件排樣
    軟件(2016年3期)2016-05-16 06:32:32
    基于模糊自適應(yīng)模擬退火遺傳算法的配電網(wǎng)故障定位
    SOA結(jié)合模擬退火算法優(yōu)化電容器配置研究
    基于遺傳-模擬退火算法的城市軌道交通快慢車停站方案
    热99re8久久精品国产| 免费日韩欧美在线观看| 国产男人的电影天堂91| 欧美国产精品va在线观看不卡| 一边摸一边抽搐一进一出视频| 老熟妇仑乱视频hdxx| 午夜免费观看性视频| 免费女性裸体啪啪无遮挡网站| 国产麻豆69| 国产av国产精品国产| 在线 av 中文字幕| 欧美av亚洲av综合av国产av| 欧美激情极品国产一区二区三区| 狠狠狠狠99中文字幕| 人成视频在线观看免费观看| 性色av乱码一区二区三区2| 亚洲三区欧美一区| 亚洲欧美激情在线| 亚洲精品国产区一区二| 老司机影院成人| 曰老女人黄片| 动漫黄色视频在线观看| 亚洲av成人不卡在线观看播放网 | 亚洲国产欧美一区二区综合| 亚洲avbb在线观看| 欧美精品啪啪一区二区三区 | 国产av又大| 精品一区在线观看国产| 亚洲熟女毛片儿| 久久毛片免费看一区二区三区| 在线看a的网站| 国产97色在线日韩免费| 一本久久精品| 国产日韩欧美在线精品| 国产有黄有色有爽视频| 一级片'在线观看视频| 中文字幕精品免费在线观看视频| 女人高潮潮喷娇喘18禁视频| 免费观看av网站的网址| 国产男女内射视频| 一区二区三区激情视频| 19禁男女啪啪无遮挡网站| 欧美日韩视频精品一区| 国产免费视频播放在线视频| 成人国产一区最新在线观看| 免费观看人在逋| av片东京热男人的天堂| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 亚洲性夜色夜夜综合| 国产精品秋霞免费鲁丝片| 在线观看免费高清a一片| 午夜视频精品福利| 亚洲九九香蕉| 亚洲综合色网址| 老司机影院成人| 日韩视频在线欧美| 欧美在线一区亚洲| 亚洲av电影在线观看一区二区三区| 少妇 在线观看| 美女午夜性视频免费| 巨乳人妻的诱惑在线观看| 成年女人毛片免费观看观看9 | 9191精品国产免费久久| 国产日韩欧美视频二区| 亚洲国产毛片av蜜桃av| 久久av网站| 视频区欧美日本亚洲| 国内毛片毛片毛片毛片毛片| 中国美女看黄片| 欧美日韩成人在线一区二区| 国产色视频综合| 成人国语在线视频| 脱女人内裤的视频| 曰老女人黄片| 亚洲欧美激情在线| 亚洲伊人久久精品综合| 免费久久久久久久精品成人欧美视频| 青青草视频在线视频观看| 亚洲va日本ⅴa欧美va伊人久久 | 97在线人人人人妻| kizo精华| 老司机影院毛片| av网站免费在线观看视频| 超色免费av| 91国产中文字幕| 久久这里只有精品19| av网站在线播放免费| 搡老乐熟女国产| 日日夜夜操网爽| 女人爽到高潮嗷嗷叫在线视频| 久久久久国产精品人妻一区二区| 天天躁夜夜躁狠狠躁躁| 波多野结衣av一区二区av| 人妻一区二区av| 久久人妻福利社区极品人妻图片| 麻豆av在线久日| 视频区图区小说| 国产高清videossex| av电影中文网址| 天天躁夜夜躁狠狠躁躁| 女人精品久久久久毛片| 精品卡一卡二卡四卡免费| 久久久国产欧美日韩av| 中国美女看黄片| 久久人人爽av亚洲精品天堂| av天堂在线播放| 俄罗斯特黄特色一大片| 久久ye,这里只有精品| 国产精品一二三区在线看| 十八禁网站网址无遮挡| 亚洲精品乱久久久久久| 秋霞在线观看毛片| 亚洲国产看品久久| 他把我摸到了高潮在线观看 | 老司机午夜福利在线观看视频 | 欧美成狂野欧美在线观看| 成人国产av品久久久| 国产成人免费观看mmmm| 成人影院久久| 亚洲欧美一区二区三区久久| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品av久久久久免费| av网站免费在线观看视频| 色婷婷av一区二区三区视频| 亚洲国产中文字幕在线视频| 亚洲国产精品成人久久小说| 精品一区二区三卡| 中文欧美无线码| 真人做人爱边吃奶动态| 日韩视频在线欧美| 午夜精品久久久久久毛片777| 精品卡一卡二卡四卡免费| 巨乳人妻的诱惑在线观看| 亚洲精品日韩在线中文字幕| 国精品久久久久久国模美| 日韩欧美国产一区二区入口| 亚洲精华国产精华精| 国产一区二区 视频在线| 国产高清视频在线播放一区 | 久久久久国产精品人妻一区二区| 视频区欧美日本亚洲| 可以免费在线观看a视频的电影网站| 欧美精品亚洲一区二区| 国产在线免费精品| 2018国产大陆天天弄谢| 麻豆乱淫一区二区| 国产精品av久久久久免费| 国产成人啪精品午夜网站| 精品国产超薄肉色丝袜足j| 视频区欧美日本亚洲| 777米奇影视久久| 国产亚洲av高清不卡| 亚洲国产日韩一区二区| 午夜久久久在线观看| 妹子高潮喷水视频| 肉色欧美久久久久久久蜜桃| 9色porny在线观看| 国产日韩欧美在线精品| 精品熟女少妇八av免费久了| 亚洲精品国产av成人精品| 一区福利在线观看| 操出白浆在线播放| 欧美日韩av久久| 亚洲成av片中文字幕在线观看| 99香蕉大伊视频| 91九色精品人成在线观看| 精品视频人人做人人爽| 久久香蕉激情| 久久久久国产精品人妻一区二区| 叶爱在线成人免费视频播放| 亚洲国产欧美在线一区| 在线观看舔阴道视频| 天天躁夜夜躁狠狠躁躁| 青春草亚洲视频在线观看| 精品久久蜜臀av无| 精品一区二区三区av网在线观看 | 国产欧美日韩综合在线一区二区| 久久99热这里只频精品6学生| 亚洲欧美精品综合一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 老汉色av国产亚洲站长工具| 1024视频免费在线观看| 国产精品 欧美亚洲| 97人妻天天添夜夜摸| 人人妻,人人澡人人爽秒播| videosex国产| 国产av又大| 看免费av毛片| 午夜激情av网站| 成年动漫av网址| a级毛片在线看网站| 婷婷丁香在线五月| 中文字幕精品免费在线观看视频| 男女无遮挡免费网站观看| 国产黄频视频在线观看| 成人18禁高潮啪啪吃奶动态图| 大片免费播放器 马上看| 久久久久久人人人人人| 亚洲黑人精品在线| av一本久久久久| 天堂俺去俺来也www色官网| 精品国产一区二区三区久久久樱花| 国产片内射在线| 亚洲伊人色综图| 中亚洲国语对白在线视频| 国产1区2区3区精品| 久久精品熟女亚洲av麻豆精品| 麻豆av在线久日| 热re99久久国产66热| 考比视频在线观看| 少妇人妻久久综合中文| 国产高清视频在线播放一区 | 蜜桃国产av成人99| 日本a在线网址| 国产片内射在线| 99国产极品粉嫩在线观看| a级毛片在线看网站| 丝袜美足系列| 一二三四社区在线视频社区8| 777米奇影视久久| 国产高清视频在线播放一区 | 黄色 视频免费看| 精品免费久久久久久久清纯 | 亚洲欧洲日产国产| 久久精品国产a三级三级三级| 免费看十八禁软件| 岛国在线观看网站| 精品国产一区二区三区久久久樱花| 制服人妻中文乱码| av欧美777| 国产精品秋霞免费鲁丝片| 久久热在线av| 日韩精品免费视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 亚洲中文字幕日韩| 老鸭窝网址在线观看| 久久国产亚洲av麻豆专区| 80岁老熟妇乱子伦牲交| 精品久久蜜臀av无| 新久久久久国产一级毛片| 两性夫妻黄色片| 欧美精品啪啪一区二区三区 | 久久综合国产亚洲精品| 国产高清视频在线播放一区 | av一本久久久久| 国产av又大| 亚洲欧洲精品一区二区精品久久久| 99国产精品一区二区三区| 老熟女久久久| 国产精品.久久久| 日韩有码中文字幕| 亚洲成人手机| 欧美日韩黄片免| 在线精品无人区一区二区三| 免费观看a级毛片全部| 日日爽夜夜爽网站| 亚洲国产欧美日韩在线播放| 欧美老熟妇乱子伦牲交| 免费观看a级毛片全部| 热re99久久国产66热| 精品亚洲乱码少妇综合久久| 大陆偷拍与自拍| 日韩大片免费观看网站| 在线 av 中文字幕| 黄色视频不卡| 人妻久久中文字幕网| 麻豆乱淫一区二区| 亚洲精品美女久久av网站| 久久久精品国产亚洲av高清涩受| 美女福利国产在线| 色视频在线一区二区三区| 精品国产超薄肉色丝袜足j| 国产精品一区二区在线观看99| 欧美日韩福利视频一区二区| 亚洲精品久久午夜乱码| 亚洲欧美色中文字幕在线| 视频在线观看一区二区三区| 亚洲成人手机| 99热全是精品| 欧美激情高清一区二区三区| 国产高清国产精品国产三级| 国产精品久久久久久精品古装| 国产精品.久久久| 操出白浆在线播放| 亚洲一区中文字幕在线| 亚洲中文av在线| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利视频精品| 欧美精品高潮呻吟av久久| 免费看十八禁软件| 国产av国产精品国产| 少妇精品久久久久久久| 免费看十八禁软件| 真人做人爱边吃奶动态| 国产真人三级小视频在线观看| 亚洲精品久久久久久婷婷小说| 黄色毛片三级朝国网站| 国产xxxxx性猛交| 国产精品偷伦视频观看了| 99久久综合免费| 精品一区二区三卡| 欧美成狂野欧美在线观看| 在线观看舔阴道视频| 久久久精品免费免费高清| 久久香蕉激情| 欧美人与性动交α欧美软件| 在线永久观看黄色视频| 久久精品久久久久久噜噜老黄| 在线十欧美十亚洲十日本专区| 一区二区日韩欧美中文字幕| 少妇裸体淫交视频免费看高清 | 国产精品九九99| 亚洲伊人久久精品综合| 两性午夜刺激爽爽歪歪视频在线观看 | 丝袜美足系列| 精品亚洲成a人片在线观看| 久久人妻熟女aⅴ| 伦理电影免费视频| 国产成人欧美| 精品久久久久久久毛片微露脸 | 亚洲一区中文字幕在线| 国产亚洲精品久久久久5区| 中文字幕高清在线视频| 国产成人精品在线电影| 91精品国产国语对白视频| 精品少妇黑人巨大在线播放| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女毛片av久久网站| 亚洲第一欧美日韩一区二区三区 | 成年人午夜在线观看视频| av天堂在线播放| 欧美在线一区亚洲| 久久热在线av| 中文字幕精品免费在线观看视频| 久久亚洲国产成人精品v| e午夜精品久久久久久久| 一区二区三区激情视频| 男女之事视频高清在线观看| 我的亚洲天堂| 12—13女人毛片做爰片一| 老司机影院毛片| 男女国产视频网站| 18禁观看日本| 亚洲精品自拍成人| 国产成人精品久久二区二区免费| 国产主播在线观看一区二区| 国产欧美日韩综合在线一区二区| 亚洲国产精品成人久久小说| www.av在线官网国产| 一级片'在线观看视频| 久久久久久久大尺度免费视频| 69av精品久久久久久 | 免费一级毛片在线播放高清视频 | 亚洲av日韩精品久久久久久密| 国产成人av激情在线播放| 欧美黄色片欧美黄色片| 国产精品一区二区免费欧美 | 在线观看免费高清a一片| 国产欧美亚洲国产| 中文精品一卡2卡3卡4更新| 十八禁高潮呻吟视频| 国产成人一区二区三区免费视频网站| 99国产极品粉嫩在线观看| 国产精品秋霞免费鲁丝片| 嫩草影视91久久| 亚洲国产av新网站| 无遮挡黄片免费观看| 日本vs欧美在线观看视频| 不卡一级毛片| 一进一出抽搐动态| 男女午夜视频在线观看| 91麻豆av在线| 两个人免费观看高清视频| 国产男人的电影天堂91| 国产野战对白在线观看| 国产精品偷伦视频观看了| 欧美人与性动交α欧美精品济南到| 色婷婷av一区二区三区视频| 黑人猛操日本美女一级片| 丝袜在线中文字幕| 电影成人av| 亚洲专区字幕在线| 欧美激情 高清一区二区三区| 日韩中文字幕视频在线看片| 亚洲专区中文字幕在线| 制服诱惑二区| 国产日韩欧美视频二区| 日本a在线网址| 99国产精品99久久久久| 亚洲精品一二三| 成年女人毛片免费观看观看9 | 日韩欧美免费精品| 久久99热这里只频精品6学生| 久久亚洲精品不卡| 国产精品秋霞免费鲁丝片| 啦啦啦啦在线视频资源| 欧美精品一区二区大全| 王馨瑶露胸无遮挡在线观看| 欧美老熟妇乱子伦牲交| 熟女少妇亚洲综合色aaa.| 极品少妇高潮喷水抽搐| 国产黄频视频在线观看| 下体分泌物呈黄色| 中文字幕另类日韩欧美亚洲嫩草| 精品卡一卡二卡四卡免费| 久久精品人人爽人人爽视色| 一区二区三区四区激情视频| 久久精品aⅴ一区二区三区四区| 日本撒尿小便嘘嘘汇集6| 欧美亚洲日本最大视频资源| 国产日韩欧美视频二区| 国产又色又爽无遮挡免| 菩萨蛮人人尽说江南好唐韦庄| 一边摸一边做爽爽视频免费| 成人免费观看视频高清| 亚洲欧洲精品一区二区精品久久久| 国产在线观看jvid| 久久久精品国产亚洲av高清涩受| 激情视频va一区二区三区| 久久精品国产亚洲av高清一级| 精品高清国产在线一区| 涩涩av久久男人的天堂| 国产欧美日韩综合在线一区二区| 激情视频va一区二区三区| 午夜福利在线观看吧| 免费在线观看视频国产中文字幕亚洲 | 久久人妻熟女aⅴ| 在线 av 中文字幕| 免费人妻精品一区二区三区视频| 亚洲精品国产av蜜桃| e午夜精品久久久久久久| 亚洲第一av免费看| av天堂久久9| 亚洲精品美女久久av网站| 丝袜美足系列| 大型av网站在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 97人妻天天添夜夜摸| 大香蕉久久成人网| 韩国精品一区二区三区| 老司机福利观看| 国产欧美日韩一区二区三 | 国产精品久久久人人做人人爽| 欧美另类一区| 在线天堂中文资源库| 王馨瑶露胸无遮挡在线观看| 国产精品秋霞免费鲁丝片| 麻豆国产av国片精品| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久免费高清国产稀缺| 99精品久久久久人妻精品| 人人妻,人人澡人人爽秒播| 可以免费在线观看a视频的电影网站| 日韩欧美一区视频在线观看| 国产成人av激情在线播放| 视频区欧美日本亚洲| 超碰成人久久| 黄片大片在线免费观看| av线在线观看网站| 大片电影免费在线观看免费| 欧美激情 高清一区二区三区| 亚洲色图综合在线观看| 啦啦啦 在线观看视频| 黄色怎么调成土黄色| 精品一区二区三卡| 亚洲男人天堂网一区| 精品一区在线观看国产| 91成人精品电影| 一进一出抽搐动态| 亚洲全国av大片| 免费在线观看黄色视频的| 黄网站色视频无遮挡免费观看| 亚洲人成电影观看| 亚洲熟女毛片儿| 欧美精品一区二区免费开放| 婷婷色av中文字幕| 别揉我奶头~嗯~啊~动态视频 | 最新的欧美精品一区二区| 美女国产高潮福利片在线看| 老鸭窝网址在线观看| 午夜两性在线视频| 亚洲五月婷婷丁香| 妹子高潮喷水视频| 久久99一区二区三区| 手机成人av网站| 亚洲激情五月婷婷啪啪| 老司机在亚洲福利影院| 交换朋友夫妻互换小说| 99热全是精品| 一区在线观看完整版| 三上悠亚av全集在线观看| 免费在线观看黄色视频的| 少妇粗大呻吟视频| 午夜老司机福利片| 少妇精品久久久久久久| 成人三级做爰电影| 午夜免费观看性视频| 亚洲精品国产色婷婷电影| 久久av网站| 免费在线观看日本一区| 亚洲av美国av| 亚洲第一欧美日韩一区二区三区 | 女人爽到高潮嗷嗷叫在线视频| 久久久国产精品麻豆| 精品卡一卡二卡四卡免费| videosex国产| 亚洲成人手机| 中文字幕另类日韩欧美亚洲嫩草| 午夜91福利影院| 成年美女黄网站色视频大全免费| 亚洲一区中文字幕在线| 最近最新免费中文字幕在线| 亚洲精品成人av观看孕妇| 成年人黄色毛片网站| 亚洲全国av大片| 最新在线观看一区二区三区| 久久狼人影院| 法律面前人人平等表现在哪些方面 | 啦啦啦 在线观看视频| 日本av手机在线免费观看| 黄色视频不卡| 中文字幕av电影在线播放| 老司机影院毛片| 久久精品亚洲熟妇少妇任你| 99精国产麻豆久久婷婷| 欧美精品啪啪一区二区三区 | 啦啦啦中文免费视频观看日本| 久久精品人人爽人人爽视色| 亚洲 欧美一区二区三区| 在线天堂中文资源库| 99久久99久久久精品蜜桃| 久久人妻熟女aⅴ| 手机成人av网站| 欧美在线一区亚洲| 久热爱精品视频在线9| 午夜福利在线观看吧| 国产一级毛片在线| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻1区二区| 狂野欧美激情性xxxx| 巨乳人妻的诱惑在线观看| 亚洲精品国产一区二区精华液| a级片在线免费高清观看视频| 一级毛片女人18水好多| 久久女婷五月综合色啪小说| 国产成+人综合+亚洲专区| 精品国产乱码久久久久久男人| 下体分泌物呈黄色| 桃红色精品国产亚洲av| 黄色 视频免费看| 99久久精品国产亚洲精品| 精品视频人人做人人爽| 久久久久久久久免费视频了| 日韩熟女老妇一区二区性免费视频| 亚洲欧美一区二区三区久久| 大型av网站在线播放| 天堂俺去俺来也www色官网| av福利片在线| videos熟女内射| 色综合欧美亚洲国产小说| 国产主播在线观看一区二区| 91大片在线观看| 亚洲精品美女久久av网站| netflix在线观看网站| 国产精品久久久久久人妻精品电影 | 精品一区在线观看国产| 91精品三级在线观看| netflix在线观看网站| 国产成人精品久久二区二区91| 亚洲欧美精品综合一区二区三区| 91国产中文字幕| 欧美日韩黄片免| 一区二区av电影网| 精品人妻一区二区三区麻豆| 成人三级做爰电影| 国产一卡二卡三卡精品| 亚洲欧美日韩高清在线视频 | 国产免费福利视频在线观看| 国产视频一区二区在线看| 欧美日韩视频精品一区| 伊人久久大香线蕉亚洲五| 菩萨蛮人人尽说江南好唐韦庄| av在线老鸭窝| 国产男女内射视频| 亚洲av美国av| 高潮久久久久久久久久久不卡| 久久久久精品人妻al黑| 亚洲中文日韩欧美视频| 性色av乱码一区二区三区2| 狂野欧美激情性bbbbbb| 两性夫妻黄色片| 日本五十路高清| 麻豆国产av国片精品| 久久精品国产a三级三级三级| 女警被强在线播放| 亚洲精品一区蜜桃| 人妻久久中文字幕网| 80岁老熟妇乱子伦牲交| 91老司机精品| 91成年电影在线观看| 不卡av一区二区三区| 久久久久国内视频| 久久久久久人人人人人| 亚洲精品乱久久久久久| 久久久欧美国产精品| 中文字幕av电影在线播放| 亚洲欧美精品自产自拍| 少妇 在线观看| 精品国产乱码久久久久久小说| 精品人妻1区二区| 欧美精品av麻豆av| 国产免费福利视频在线观看| 国产深夜福利视频在线观看| 国产精品一区二区免费欧美 |