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

    改進(jìn)型負(fù)輸入整形與最優(yōu)控制結(jié)合的振動(dòng)抑制方法

    2014-08-11 14:50:07趙志剛游斌弟
    振動(dòng)與沖擊 2014年8期
    關(guān)鍵詞:最優(yōu)控制改進(jìn)型柔性

    趙志剛, 游斌弟, 趙 陽(yáng)

    (1.中國(guó)民航大學(xué) 空中交通管理學(xué)院,天津 300300;2.哈爾濱工業(yè)大學(xué) 航天學(xué)院,哈爾濱 150001;3.哈爾濱工業(yè)大學(xué) (威海)船舶與海洋工程學(xué)院,山東 威海 264209)

    改進(jìn)型負(fù)輸入整形與最優(yōu)控制結(jié)合的振動(dòng)抑制方法

    趙志剛1,2, 游斌弟3, 趙 陽(yáng)2

    (1.中國(guó)民航大學(xué) 空中交通管理學(xué)院,天津 300300;2.哈爾濱工業(yè)大學(xué) 航天學(xué)院,哈爾濱 150001;3.哈爾濱工業(yè)大學(xué) (威海)船舶與海洋工程學(xué)院,山東 威海 264209)

    針對(duì)諸如柔性機(jī)械臂這類柔性系統(tǒng)的主動(dòng)振動(dòng)控制問題,提出了基于改進(jìn)型負(fù)輸入整形和最優(yōu)控制結(jié)合的振動(dòng)抑制方法。以柔性機(jī)械臂為研究對(duì)象,設(shè)計(jì)最優(yōu)狀態(tài)反饋控制實(shí)現(xiàn)其旋轉(zhuǎn)機(jī)動(dòng)任務(wù)及柔性振動(dòng)抑制。為增強(qiáng)柔性振動(dòng)抑制效果和改善系統(tǒng)的運(yùn)動(dòng)時(shí)間,根據(jù)引入線性二次型調(diào)節(jié)器(LQR)反饋后整個(gè)閉環(huán)系統(tǒng)的振動(dòng)頻率和阻尼比設(shè)計(jì)改進(jìn)型負(fù)輸入整形器作為前饋控制器。將前饋控制與反饋控制相結(jié)合能發(fā)揮其各自的優(yōu)點(diǎn),提高系統(tǒng)的性能。仿真分析結(jié)果表明,所設(shè)計(jì)的混合控制策略可以有效地抑制柔性振動(dòng),且可以減少系統(tǒng)響應(yīng)時(shí)間的延遲和加快系統(tǒng)響應(yīng)速度。

    主動(dòng)振動(dòng)控制;改進(jìn)型負(fù)輸入整形;最優(yōu)控制;線性二次型調(diào)節(jié)器;柔性機(jī)械臂;振動(dòng)抑制

    現(xiàn)代工業(yè)和航空航天等領(lǐng)域廣泛應(yīng)用柔性結(jié)構(gòu)來滿足輕型化、柔性化的設(shè)計(jì)要求。其中,柔性機(jī)械臂(如高空作業(yè)車的伸縮臂,空間機(jī)械臂等)作為一種典型代表常用來完成各種作業(yè)任務(wù)。然而,柔性機(jī)械臂在執(zhí)行操作任務(wù)時(shí)容易產(chǎn)生振動(dòng),影響其定位精度和工作效率,嚴(yán)重時(shí)甚至危害機(jī)械臂系統(tǒng)結(jié)構(gòu)的安全。解決柔性機(jī)械臂的柔性振動(dòng)抑制問題可以提高其工作性能,增強(qiáng)系統(tǒng)的安全性和可靠性,并且可以降低系統(tǒng)的操作成本。因此,這一課題具有重要的工程意義,已引起了國(guó)內(nèi)外許多研究者的興趣。

    針對(duì)柔性結(jié)構(gòu)的振動(dòng)抑制問題,國(guó)內(nèi)外學(xué)者提出了許多種研究方法[1-2],如最優(yōu)控制、自適應(yīng)控制、魯棒控制、奇異攝動(dòng)控制、模糊控制、滑??刂频?,上述方法及其復(fù)合控制方法等獲得了較深入的研究。其中,最優(yōu)控制是通過求解最優(yōu)控制力來使系統(tǒng)振動(dòng)達(dá)到理想的抑制效果,且兼顧響應(yīng)與控制兩方面相互矛盾的要求,是滿足一定最優(yōu)條件的反饋控制。對(duì)于具有二次型性能指標(biāo)的線性系統(tǒng),因其最優(yōu)反饋控制律能用解析的形式表示,易于工程實(shí)現(xiàn),歷來受到人們的重視。戈新生等[3]針對(duì)雙連桿空間剛?cè)嵝詸C(jī)械臂的殘余振動(dòng)問題,采用LQ最優(yōu)控制方法來消除殘余振動(dòng),通過數(shù)值仿真表明所提出方法的有效性。

    除了上面列舉的控制方法之外,還有一種越來越吸引人們興趣的主動(dòng)振動(dòng)控制方法——輸入整形(Input Shaping,IS)(又稱為輸入成形或輸入成型)振動(dòng)抑制方法。該方法作為一種典型的前饋控制方法,自從Singer等[4]正式提出以來就得到了廣泛而深入的研究。近年來,國(guó)內(nèi)學(xué)者在此領(lǐng)域也取得了豐富的研究成果。董明曉等[5]設(shè)計(jì)最優(yōu)輸入整形器(Input Shaper),應(yīng)用于變參數(shù)橋式起重機(jī)等的防擺控制中。張鵬等[6]通過引入偏值點(diǎn),提出了基于零點(diǎn)配置技術(shù)的EI成型器設(shè)計(jì)方法。輸入整形作為一種簡(jiǎn)單、高效的前饋補(bǔ)償器與反饋控制相結(jié)合時(shí),不僅可以提高控制性能,甚至還可以降低反饋控制器的設(shè)計(jì)難度。胡慶雷等[7-10]將輸入成形方法與滑模變結(jié)構(gòu)控制相結(jié)合應(yīng)用于撓性航天器大角度姿態(tài)機(jī)動(dòng)時(shí)撓性結(jié)構(gòu)的振動(dòng)控制。Pai[11]提出了輸入整形結(jié)合自適應(yīng)滑模控制的方案解決帶有參數(shù)不確定性和外擾的柔性結(jié)構(gòu)振動(dòng)控制問題??讘椚实萚12]將輸入成形和PD反饋控制器相結(jié)合來解決柔性航天器的振動(dòng)抑制問題。李瑞川等[13]將PD控制結(jié)合輸入整形控制策略應(yīng)用到含有剛性模態(tài)的彈性機(jī)構(gòu)中,研究了輸入整形參數(shù)對(duì)系統(tǒng)控制力的作用規(guī)律以及降低系統(tǒng)能耗的機(jī)理。王曉軍等[14]研究了輸入整形技術(shù)結(jié)合線形二次型調(diào)節(jié)器(LQR)的混合控制方案來解決橋式起重機(jī)這一類不完全控制系統(tǒng)的負(fù)載殘留振蕩抑制問題。上述研究所取得的成果大多都是基于傳統(tǒng)的輸入整形設(shè)計(jì)方法,陜晉軍等[15]通過放寬設(shè)計(jì)整形器需滿足包含最少脈沖個(gè)數(shù)的約束,提出了改進(jìn)型的輸入成形方法(Modified Input Shaping, MIS)并應(yīng)用該方法結(jié)合PD控制用于柔性機(jī)械臂的振動(dòng)抑制,理論分析表明在一些方面比傳統(tǒng)的輸入成形方法有更好的性能。輸入整形法能夠很好地抑制柔性系統(tǒng)振動(dòng)所付出的代價(jià)是系統(tǒng)運(yùn)動(dòng)響應(yīng)時(shí)間的延遲,為了獲得更快的運(yùn)動(dòng)響應(yīng),可以通過在輸入整形器脈沖序列中引入負(fù)脈沖來提高系統(tǒng)的上升時(shí)間。

    綜上所述,為進(jìn)一步提高控制效果、改善控制系統(tǒng)性能,尤其是提高系統(tǒng)上升時(shí)間,本文提出改進(jìn)型負(fù)輸入整形控制和最優(yōu)狀態(tài)反饋控制相結(jié)合的控制方法。首先設(shè)計(jì)LQR反饋控制器,然后設(shè)計(jì)改進(jìn)型負(fù)輸入整形器作為前饋控制器,構(gòu)成前饋控制和反饋控制的混合控制方案來解決柔性結(jié)構(gòu)的振動(dòng)抑制問題。

    1 輸入整形

    輸入整形是指由一系列脈沖序列與一定的期望輸入相卷積,所形成的指令作為控制系統(tǒng)運(yùn)動(dòng)的輸入。其中,根據(jù)振動(dòng)系統(tǒng)頻率和阻尼比所設(shè)計(jì)的脈沖序列稱之為輸入整形器[2]。其成形過程如圖1所示。

    圖1 輸入整形器對(duì)輸入信號(hào)進(jìn)行成形Fig.1 The input shaping process

    為研究輸入整形方法,考慮單自由度有阻尼二階振動(dòng)系統(tǒng):

    (1)

    式中:ω為系統(tǒng)自然振動(dòng)頻率,ζ為系統(tǒng)阻尼比。

    (2)

    其中:

    (3)

    (4)

    為使成形后的指令可以使系統(tǒng)達(dá)到與原控制指令相同的運(yùn)動(dòng)位置,還需要脈沖幅值約束條件,即

    (5)

    此外,整形器的時(shí)間長(zhǎng)度越短,給系統(tǒng)帶來的時(shí)間滯后就會(huì)越短,有利于提高系統(tǒng)的響應(yīng)速度。因此,時(shí)間約束條件為

    min(tn)

    (6)

    根據(jù)式(2)至式(6),對(duì)于確定的系統(tǒng),已知其自然頻率ω和阻尼比ζ,就可根據(jù)殘余振動(dòng)V(ω,ζ)的值為零或低于某一期望水平來確定脈沖幅值A(chǔ)i和作用時(shí)間ti。

    1.1 典型輸入整形器

    傳統(tǒng)最簡(jiǎn)單的正輸入整形器是包含兩脈沖的整形器,稱之為ZV(Zero Vibration)整形器[4],其數(shù)學(xué)表示形式為

    (7)

    其中KT、Td表示如下式

    (8)

    滿足約束條件式(2)-式(5),并在同一振動(dòng)周期內(nèi)選擇多個(gè)脈沖可以獲得MIS ZV整形器[15],對(duì)于3脈沖MIS ZV整形器可表示如下:

    (9)

    其中KM表示如下式

    (10)

    1.2 改進(jìn)型負(fù)輸入整形

    設(shè)計(jì)輸入整形器的關(guān)鍵在于確定脈沖序列中各脈沖的幅值和作用時(shí)間。除此之外,確定整形器中所包含的脈沖個(gè)數(shù)也至關(guān)重要。本文給出一種改進(jìn)型負(fù)輸入整形(Modified Negative Input Shaping,MNIS)方法,該方法放寬整形器包含最少脈沖個(gè)數(shù)的約束,且在同一振動(dòng)周期內(nèi)選擇包含正負(fù)交替的脈沖。

    滿足約束條件式(2)為零,即保證系統(tǒng)殘余振動(dòng)為零,且同時(shí)滿足式(5)的MNIS ZV整形器的脈沖作用時(shí)間選擇為

    (11)

    式中:n為整形器脈沖序列所包括脈沖個(gè)數(shù)。式(11)同時(shí)表明:隨著脈沖個(gè)數(shù)的增加,每個(gè)脈沖可以獲得更短的延遲時(shí)間。

    與脈沖作用時(shí)間式(11)相對(duì)應(yīng)的脈沖幅值為

    (12)

    式中:

    (13)

    M=-K+…+(-1)i-1Ki-1+…+(-1)n-1Kn-1=

    (14)

    式(11)-式(14)給出了改進(jìn)型負(fù)輸入整形器的封閉解,根據(jù)以上內(nèi)容,只要知道系統(tǒng)的振動(dòng)頻率與阻尼比以及所確定的脈沖數(shù)目,就可以很方便地設(shè)計(jì)MNIS ZV整形器。通過該方法所設(shè)計(jì)的整形器可以得到更短的延遲時(shí)間,且脈沖序列中所包含的脈沖個(gè)數(shù)可以靈活選擇。這給研究人員在實(shí)際工程設(shè)計(jì)中留有更多的選擇空間。例如,對(duì)于3脈沖MNIS ZV整形器,可表示如下:

    (15)

    式中:

    (16)

    2 線性二次型最優(yōu)控制

    LQ最優(yōu)控制因其產(chǎn)生狀態(tài)線性反饋的最優(yōu)控制律,易于構(gòu)成閉環(huán)最優(yōu)控制,對(duì)從事主動(dòng)振動(dòng)控制的研究人員具有很大吸引力,在工程中得到了廣泛應(yīng)用。

    LQR最優(yōu)控制是要求在滿足動(dòng)力學(xué)方程

    (17)

    且能使系統(tǒng)在消耗較少控制能量的情況下,確定最優(yōu)控制輸入u*,使下面二次型性能指標(biāo)Jm最小。

    (18)

    式中:Q是對(duì)狀態(tài)變量加權(quán)矩陣,為半正定對(duì)稱矩陣;R是對(duì)控制變量加權(quán)矩陣,為正定對(duì)稱矩陣。

    LQR最優(yōu)控制問題就是對(duì)線性時(shí)不變系統(tǒng)式(18)確定最優(yōu)控制的輸入規(guī)律,即

    其實(shí),我和阿花也是老朋友了,阿花當(dāng)年從韶關(guān)來深圳時(shí),我們就認(rèn)識(shí)了。林強(qiáng)信將雪茄掐在煙灰缸里,說,我們?cè)羞^多次合作,而且很愉快。所以你去她廠里,我贊成。你幫她,就等于我在幫她嘛。哈哈哈……

    u=-Kx

    (19)

    式中:K為狀態(tài)反饋增益矩陣。

    利用極值原理,使得二次型性能指標(biāo)式(19)最小的K表示為

    K=R-1BTP

    (20)

    其中:P為正定對(duì)稱常數(shù)矩陣,滿足下列Riccati代數(shù)方程

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

    (21)

    由式(19)和式(20),則可得LQR最優(yōu)控制為

    u*=-R-1BTPx

    (22)

    3 柔性機(jī)械臂系統(tǒng)

    以柔性機(jī)械臂作為被控對(duì)象進(jìn)行應(yīng)用研究,仿真模型如圖2所示。柔性機(jī)械臂系統(tǒng)由中心剛體、柔性臂桿及末端剛性負(fù)載mp組成。圖2中OXY和Oxy分別是慣性坐標(biāo)系和柔性臂體固聯(lián)坐標(biāo)系,Ox與未變形時(shí)柔性臂的軸線重合,其原點(diǎn)位于臂與中心剛體的連接處,θ(t)表示兩坐標(biāo)系間的相對(duì)轉(zhuǎn)動(dòng)關(guān)系。柔性臂相對(duì)坐標(biāo)系Oxy的柔性變形用δ表示。l,ρ,E,φ,Φ分別表示柔性臂長(zhǎng)度、密度、彈性模量、內(nèi)外圓直徑。Jo表示中心剛體轉(zhuǎn)動(dòng)慣量。

    圖2 柔性機(jī)械臂的簡(jiǎn)化模型Fig.2 Sketch of the flexible manipulator model

    應(yīng)用有限元法描述系統(tǒng)的彈性變形,采用拉格朗日方法建立系統(tǒng)的動(dòng)力學(xué)方程,可得由系統(tǒng)廣義自由度θ和節(jié)點(diǎn)位移列陣d表征的動(dòng)力學(xué)方程:

    (23)

    式中:J為系統(tǒng)總轉(zhuǎn)動(dòng)慣量;I為剛?cè)狁詈暇仃?;M為柔性臂質(zhì)量陣;C為柔性臂阻尼陣;K為柔性臂剛度陣;T為驅(qū)動(dòng)力矩。

    利用約束模態(tài)法展開式(23),可得柔性臂運(yùn)動(dòng)與振動(dòng)耦合的動(dòng)力學(xué)方程如下:

    (24)

    式中:ωc=diag(ωc1,…,ωcN)和ζ分別表示約束模態(tài)頻率矩陣和阻尼比矩陣;ηc=[ηc1,…,ηcN]T為約束模態(tài)坐標(biāo)向量;Hc為1×N維的模態(tài)坐標(biāo)耦合矩陣,EN為N×N維單位陣,N表示約束模態(tài)截?cái)嚯A數(shù)。

    4 改進(jìn)型負(fù)輸入整形與LQR聯(lián)合設(shè)計(jì)

    則可以將式(24)寫為

    (25)

    式中:E為相應(yīng)維數(shù)的單位陣。

    即可得狀態(tài)空間的形式方程

    y=CX

    (26)

    其中:

    Bc=[0 0 1 0]T,C=[1 0 0 0]

    進(jìn)一步改寫為如下狀態(tài)方程的形式

    y=CX

    (27)

    本文提出的改進(jìn)型負(fù)輸入整形控制和最優(yōu)狀態(tài)反饋控制相結(jié)合方法主要分為以下兩個(gè)步驟:

    第一步,設(shè)計(jì)最優(yōu)狀態(tài)反饋控制律,實(shí)現(xiàn)柔性機(jī)械臂旋轉(zhuǎn)機(jī)動(dòng)任務(wù)的同時(shí)對(duì)柔性振動(dòng)進(jìn)行抑制。

    將式(22)代入式(27),可得閉環(huán)系統(tǒng)的狀態(tài)方程為

    y=CX

    (28)

    式中:As=A-BR-1BTP為閉環(huán)系統(tǒng)矩陣。

    第二步,根據(jù)整個(gè)閉環(huán)系統(tǒng)的特征值設(shè)計(jì)改進(jìn)型負(fù)輸入整形器,輸入形成器的控制作用在LQR反饋回路的前端,目的是在保證被控對(duì)象完成旋轉(zhuǎn)機(jī)動(dòng)任務(wù)的同時(shí)進(jìn)一步抑制其柔性振動(dòng)。

    閉環(huán)系統(tǒng)矩陣As的特征值,可求出系統(tǒng)模態(tài)的振動(dòng)頻率和阻尼比

    (29)

    其中:λsi為矩陣系統(tǒng)矩陣As的第i個(gè)特征值,ωsi和ζsi分別為系統(tǒng)模態(tài)的振動(dòng)頻率和阻尼比。

    基于改進(jìn)型負(fù)輸入整形與LQR聯(lián)合設(shè)計(jì)實(shí)現(xiàn)機(jī)械臂振動(dòng)抑制的原理如圖3所示。

    圖3 實(shí)現(xiàn)振動(dòng)抑制的控制系統(tǒng)原理框圖Fig.3 Block diagram of control system for vibration suppression

    其中:r為期望輸入,K1取決于狀態(tài)反饋增益K=[K1K2…K2(N+1)]中對(duì)應(yīng)的角位移反饋系數(shù)。

    根據(jù)式(29)中閉環(huán)系統(tǒng)振動(dòng)模態(tài)的振動(dòng)頻率ωsi和阻尼比ζsi可以靈活設(shè)計(jì)一系列包含不同脈沖個(gè)數(shù)的MNIS整形器,使機(jī)械臂在改進(jìn)型負(fù)輸入整形器和LQR控制器聯(lián)合控制下完成指定旋轉(zhuǎn)任務(wù),同時(shí)實(shí)現(xiàn)振動(dòng)抑制。

    當(dāng)然,實(shí)際情況中受控對(duì)象一般包含多階振動(dòng)模態(tài),這就需要設(shè)計(jì)多模態(tài)整形器。多模態(tài)系統(tǒng)的輸入整形器設(shè)計(jì)方法一般有聯(lián)立方程法[16]和級(jí)聯(lián)法[4]??紤]到前者的計(jì)算相當(dāng)復(fù)雜,往往需要借助專門的計(jì)算軟件進(jìn)行求解和設(shè)計(jì),而后者設(shè)計(jì)簡(jiǎn)單,且對(duì)高模態(tài)的不靈敏性要比前者更高,所以本文應(yīng)用級(jí)聯(lián)法實(shí)現(xiàn)多模態(tài)整形器的設(shè)計(jì),即采用如下的卷積形式來獲得:

    Smj=S1*S2*…*Si*…*SN

    j=1,2,…,N

    (30)

    其中“*”表示進(jìn)行卷積運(yùn)算處理,Si代表對(duì)應(yīng)系統(tǒng)第i個(gè)振動(dòng)模態(tài)的輸入整形器。

    5 仿真算例

    為了進(jìn)一步驗(yàn)證本文所提出的MNIS整形器的有效性,系統(tǒng)主要仿真參數(shù):轉(zhuǎn)動(dòng)慣量Jo=40 kg·m2,臂桿長(zhǎng)度l=4 m,材料密度ρ=2 700 kg/m3,彈性模量E=70 GPa,外圓直徑Ф=0.03 m,內(nèi)圓直徑φ=0.028 m,負(fù)載質(zhì)量mp=2 kg。有限單元體個(gè)數(shù)n=10,初始狀態(tài)為零。為了對(duì)比分析仿真效果,應(yīng)用LQR控制器結(jié)合不同輸入整形器設(shè)計(jì)了四種控制方案完成參考輸入為60°機(jī)械臂旋轉(zhuǎn)機(jī)動(dòng)任務(wù),詳見表1。

    表1 四種控制方案

    考慮到系統(tǒng)存在剛體運(yùn)動(dòng)振動(dòng)和柔性模態(tài)振動(dòng),這里抑制剛體運(yùn)動(dòng)模態(tài)(ωr=0.337 2 rad/s,ζr=0.706 0)和一階柔性系統(tǒng)模態(tài)(ωf=5.517 7 rad/s,ζf=0.023 5),采用級(jí)聯(lián)法進(jìn)行設(shè)計(jì)。首先為剛體運(yùn)動(dòng)的振動(dòng)模態(tài)設(shè)計(jì)相應(yīng)的整形器Sr,然后為柔性模態(tài)振動(dòng)設(shè)計(jì)相應(yīng)的整形器Sf。根據(jù)式(30),將獲得的整形器進(jìn)行級(jí)聯(lián),可得一個(gè)新的整形器Sm=Sr*Sf,以實(shí)現(xiàn)同時(shí)對(duì)剛體運(yùn)動(dòng)的振動(dòng)和柔性模態(tài)振動(dòng)的振動(dòng)抑制。為便于比較,四種控制方案中LQR均取相同反饋增益,K=[10.00 0.025 3 41.96 1.749],K1=10.00(取LQR反饋中x1的增益);所用各種整形器參數(shù)如下:

    i=1,2

    (31)

    i=1,2,3

    (32)

    i=1,2,3

    (33)

    四種方案的仿真結(jié)果如圖4-圖9所示,分析其控制效果如下:

    四種方案均有效地完成了柔性機(jī)械臂的旋轉(zhuǎn)60°機(jī)動(dòng)任務(wù)。由圖4-圖9可以看出,單獨(dú)采用LQR控制的方案1能夠?qū)崿F(xiàn)機(jī)動(dòng)任務(wù),所用上升時(shí)間約為9.87 s,但是存在超調(diào)量,調(diào)整時(shí)間大約為23 s;而且存在殘余振動(dòng),一階模態(tài)振動(dòng)持續(xù)時(shí)間約30 s,機(jī)械臂末端彈性變形也較大,難于迅速精確定位于期望的位置。而引入整形器的方案2和方案3剛體運(yùn)動(dòng)都非常平穩(wěn),沒有超調(diào)量,所用上升時(shí)間分別約為13.7 s和18.3 s。這說明LQR通過和輸入整形技術(shù)結(jié)合,產(chǎn)生的混合方案能夠有效地減小機(jī)械臂的振動(dòng),但由于整形器時(shí)滯時(shí)間的存在而增加了系統(tǒng)的上升時(shí)間。為了降低這種不利影響,方案4基于改進(jìn)型負(fù)輸入整形與LQR聯(lián)合控制使剛體運(yùn)動(dòng)非常平穩(wěn),沒有超調(diào)量,且柔性振動(dòng)模態(tài)抑制效果也很好。所用上升時(shí)間約為9.15 s。

    方案2、方案3和方案4與方案1相比,由于輸入整形器的作用使得角速度曲線變得比較光滑,角速度變化更平穩(wěn),如圖6所示。方案4與方案2相比雖然所用整形器包含脈沖個(gè)數(shù)多,但是由于負(fù)脈沖的引入,仍然有更短的上升時(shí)間,減少了4.55 s;方案4的上升時(shí)間僅為方案3的一半,大大減少了正輸入整形器導(dǎo)致的系統(tǒng)響應(yīng)的時(shí)間延遲,這說明了方案4在改善機(jī)械臂運(yùn)動(dòng)時(shí)間上比方案2和方案3具有較大優(yōu)越性。

    綜上所述,仿真實(shí)驗(yàn)結(jié)果表明改進(jìn)型負(fù)輸入整形器和LQR控制器聯(lián)合控制方案使機(jī)械臂在下完成指定的旋轉(zhuǎn)任務(wù)的同時(shí)有效抑制了柔性振動(dòng),且在減少系統(tǒng)響應(yīng)延遲時(shí)間、提高了工作效率上優(yōu)于LQR控制器和其它輸入整形器組成的方案。

    圖4 控制輸入指令Fig.4Controlinputcommand圖5 角位移Fig.5Angulardisplacement圖6 角速度Fig.6Angularvelocity

    圖7 一階模態(tài)坐標(biāo)Fig.7Thefirstmodalcoordinate圖8 末端彈性變形Fig.8Deformationangleoftip圖9 末端彈性變形速度Fig.9Deformationvelocityoftip

    6 結(jié) 論

    以柔性機(jī)械臂為研究對(duì)象,提出了基于改進(jìn)型負(fù)輸入整形主動(dòng)振動(dòng)控制和最優(yōu)控制相結(jié)合的控制器設(shè)計(jì)方法來解決諸如柔性機(jī)械臂這類柔性系統(tǒng)的振動(dòng)抑制問題,具體工作及結(jié)論如下:

    (1) 為完成柔性機(jī)械臂旋轉(zhuǎn)機(jī)動(dòng)任務(wù),以消耗較少控制能量為性能指標(biāo)設(shè)計(jì)LQR反饋控制器,并同時(shí)進(jìn)行振動(dòng)抑制;

    (2) 為進(jìn)一步改善控制效果、提高系統(tǒng)響應(yīng)速度,根據(jù)LQR反饋后的整個(gè)閉環(huán)系統(tǒng)的振動(dòng)頻率和阻尼比設(shè)計(jì)改進(jìn)型負(fù)輸入整形器,應(yīng)用級(jí)聯(lián)法設(shè)計(jì)剛體運(yùn)動(dòng)振動(dòng)和柔性模態(tài)振動(dòng)的多模態(tài)整形器作為前饋控制器,進(jìn)行振動(dòng)抑制;

    (3) 對(duì)比文中設(shè)計(jì)的四種控制方案:改進(jìn)型負(fù)輸入整形器和LQR控制器聯(lián)合控制的方案增強(qiáng)了振動(dòng)抑制效果,有效減小系統(tǒng)響應(yīng)超調(diào)量,在改善機(jī)械臂的運(yùn)動(dòng)時(shí)間,提高工作效率上優(yōu)于LQR控制器和其它輸入整形器組成的方案;

    (4) 進(jìn)一步工作重點(diǎn)是將本文方法應(yīng)用于柔性結(jié)構(gòu)振動(dòng)抑制試驗(yàn)研究中。此外,如果系統(tǒng)的轉(zhuǎn)動(dòng)慣量等參數(shù)的實(shí)際值和名義值不同,即模型參數(shù)不精確,這將影響系統(tǒng)的頻率和阻尼比信息的獲得。根據(jù)不精確信息設(shè)計(jì)的整形器可能會(huì)導(dǎo)致系統(tǒng)殘余振動(dòng)不能被完全消除,這就需要設(shè)計(jì)具有魯棒性的改進(jìn)型負(fù)輸入整形器來解決此問題,這也是進(jìn)一步要研究的問題。

    [ 1 ] Benosman M, Vey G L. Control of flexible manipulators: A survey[J]. Robotica, 2004, 22(5), 533-545.

    [ 2 ] 王樹新, 員今天, 石菊榮, 等. 柔性機(jī)械臂建模理論與控制方法研究綜述[J]. 機(jī)器人, 2002, 24(1), 86-91, 96. WANG Shu-xin, YUAN Jin-tian, SHI ju-rong, et al. A roadmap of research on modeling and control strategy for flexible manipulators[J]. Robot, 2002, 24(1), 86-91, 96.

    [ 3 ] 戈新生, 姜兵利, 劉延柱.空間剛?cè)嵝詸C(jī)械臂振動(dòng)抑制的LQ最優(yōu)控制方法[J]. 振動(dòng)與沖擊, 1999, 18(1):49-52, 56. GE Xin-sheng, JIANG Bing-li, LIU Yan-zhu. The LQ optimal control method for vibration suppression of rigid flexible space manipulator[J]. Journal of Vibration and Shock, 1999, 18(1): 49-52, 56.

    [ 4 ] Singer N C, Seering W P. Preshaping command inputs to reduce system vibration[J]. Journal of Dynamics Systems, Measurement and Control, 1990, 112(3): 76-82.

    [ 5 ] 董明曉, 脫建智, 任意翔, 等. 最優(yōu)輸入整形抑制變參數(shù)橋式起重機(jī)載荷擺動(dòng)[J]. 振動(dòng)與沖擊, 2009, 28(10): 207-209. DONG Ming-xiao, TUO Jian-zhi, REN Yi-xiang, et al. Optimal input-shaper to suppress payload oscillation of bridge cranes with varying parameters[J]. Journal of Vibration and Shock, 2009, 28(10): 207-209.

    [ 6 ] 張鵬, 李元春. 基于混合EI成型器的多模態(tài)柔性結(jié)構(gòu)振動(dòng)控制[J]. 力學(xué)學(xué)報(bào), 2010, 42(4): 774-781. ZHANG Peng, LI Yuan-chun. Vibration control of multi-mode flexible structure based on hybrid extra-insensitive shaper[J]. Chinese Journal of Theoretical and Applied Mechanics, 2010, 42(4): 774-781.

    [ 7 ] 胡慶雷, 馬廣富.基于滑模輸出反饋與輸入成形控制相結(jié)合的撓性航天器主動(dòng)振動(dòng)抑制方法[J]. 振動(dòng)與沖擊, 2007, 26 (6):133-138, 189. HU Qing-lei, MA Guang-fu. Active vibration suppression in flexible spacecraft with mismatched uncertainty via sliding mode/shaped input control[J]. Journal of Vibration and Shock, 2007, 26 (6):133-138, 189.

    [ 8 ] HU Qing-lei, SHI Peng, GAO Hui-jun. Adaptive variable structure and commanding shaped vibration control of flexible spacecraft[J]. Journal of Guidance, Control, and Dynamics, 2007, 30(3): 804-815.

    [ 9 ] HU Qing-lei, MA Guang-fu. Active vibration suppression in flexible spacecraft attitude maneuver using variable structure control and input shaping technique[C]. AIAA Guidance, Navigation, and Control Conference and Exhibit, San Francisco, California, Aug. 15-18, AIAA, 2005:2005-6246.

    [10] HU Qing-lei. Input shaping and variable structure control for simultaneous precision positioning and vibration reduction of flexible spacecraft with saturation compensation[J]. Journal of Sound and Vibration, 2008, 318 (1-2): 18-35.

    [11] Pai M C. Closed-loop input shaping control of vibration in flexible structures via adaptive sliding mode control[J]. Shock and Vibration, 2012, 19(2): 217-229.

    [12] 孔憲仁, 楊正賢, 葉東, 等.基于輸入成形的柔性航天器振動(dòng)閉環(huán)抑制方法研究[J]. 振動(dòng)與沖擊, 2010, 29(3): 72-76.

    KONG Xian-ren, YANG Zheng-xian,Ye Dong, et al. Feedback control in conjunction with input shaping for flexible spacecraft vibration suppression[J]. Journal of Vibration and Shock, 2010, 29(3): 72-76.

    [13] 李瑞川, 徐濤, 左文杰, 等. 基于PD結(jié)合輸入整形控制的含有剛性模態(tài)彈性機(jī)構(gòu)的力及能耗分析[J]. 振動(dòng)與沖擊, 2010, 29(12): 158-161. LI Rui-chuan, XU Tao, ZUO Wen-jie, et al. Analysis on force and energy consumption for elastic mechanisms with a rigid mode motion based on PD combined with input-shaping control[J]. Journal of Vibration and Shock, 2010, 29(12): 158-161.

    [14] 王曉軍, 邵惠鶴. 基于輸入整形的線性二次型調(diào)節(jié)器及其應(yīng)用[J]. 上海交通大學(xué)學(xué)報(bào), 2006, 40 (5): 848-851, 855. WANG Xiao-jun, SHAO Hui-he. The linear quadratic regulator based on input shaping technique and its application[J]. Journal of Shanghai Jiaotong University, 2006, 40 (5): 848-851, 855.

    [15] Shan J J, Liu H T, Sun D. Modified input shaping for a rotating single-link flexible manipulator[J]. Journal of Sound and Vibration, 2005, 285(1-2): 187-207.

    [16] Hyde J M, Seering W P. Inhibiting multiple mode vibration in controlled flexible systems[C]. American Control Conference. Boston: IEEE, 1991: 2449-2454.

    Combination of modified negative input shaping and optimal control for vibration suppression

    ZHAO Zhi-gang1,2, YOU Bin-di3, ZHAO Yang2

    (1.College of Air Traffic Management, Civil Aviation University of China, Tianjin 300300, China;2. School of Astronautics,Harbin Institute of Technology, Harbin 150001,China; 3. School of Naval Architecture and Ocean Engineering, Harbin Institute of Technology at Weihai, Weihai 264209, China)

    A new vibration reduction control strategy was presented. With it, the modified negative input shaping (MNIS) technique was combined with the optimal control method for active vibration control of flexible systems, such as, a flexible manipulator. An optimal state feedback controller was designed for the flexible manipulator maneuver and its vibration eliminating. To improve the performance of vibration control and reduce the system response time duration, a modified negative input shaper was used as a feed-forward controller, it was designed to utilize the vibration frequency and damping ratio of the whole closed loop system with a linear quadratic regulator (LQR). Combining the feed-forward control and feedback control and utilizing the respective advantages, the performance of the control system can be improved. Both analytical and numerical results were presented to demonstrate the effectiveness of the hybrid control strategy. It was shown that the proposed method can be used to reduce the time delay of the system’s response and increase the response speed of the system.

    active vibration control; modified negative input shaping (MNIS); optimal control; linear quadratic regulator (LQR); flexible manipulator; vibration suppression

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目(2013CB733000);國(guó)家自然科學(xué)基金資助項(xiàng)目(50975056,11072066);國(guó)家青年科學(xué)基金資助項(xiàng)目(51205079)

    2012-12-06 修改稿收到日期:2013-03-23

    趙志剛 男,博士生, 1982年生

    V414; TP271

    A

    10.13465/j.cnki.jvs.2014.08.035

    猜你喜歡
    最優(yōu)控制改進(jìn)型柔性
    一種柔性拋光打磨頭設(shè)計(jì)
    Cr5改進(jìn)型支承輥探傷無底波原因分析
    灌注式半柔性路面研究進(jìn)展(1)——半柔性混合料組成設(shè)計(jì)
    石油瀝青(2021年5期)2021-12-02 03:21:18
    條件平均場(chǎng)隨機(jī)微分方程的最優(yōu)控制問題
    高校學(xué)生管理工作中柔性管理模式應(yīng)用探索
    帶跳躍平均場(chǎng)倒向隨機(jī)微分方程的線性二次最優(yōu)控制
    改進(jìn)型CKF算法及其在GNSS/INS中的應(yīng)用
    Timoshenko梁的邊界最優(yōu)控制
    采用最優(yōu)控制無功STATCOM 功率流的解決方案
    改進(jìn)型逆變器無效開關(guān)死區(qū)消除方法
    久久国产精品影院| 1024视频免费在线观看| 亚洲免费av在线视频| 精品国内亚洲2022精品成人 | 午夜福利在线观看吧| 一本—道久久a久久精品蜜桃钙片| 久久 成人 亚洲| 热re99久久精品国产66热6| 国产在线观看jvid| 狠狠婷婷综合久久久久久88av| 精品国产一区二区久久| 少妇猛男粗大的猛烈进出视频| 国产亚洲av高清不卡| 午夜福利乱码中文字幕| 啪啪无遮挡十八禁网站| 伦理电影免费视频| 亚洲五月色婷婷综合| 国产真人三级小视频在线观看| 美女脱内裤让男人舔精品视频| 国产国语露脸激情在线看| 人人妻人人添人人爽欧美一区卜| 成人av一区二区三区在线看 | 狂野欧美激情性xxxx| 99久久综合免费| 欧美 日韩 精品 国产| 久久av网站| 美女大奶头黄色视频| 色综合欧美亚洲国产小说| 丰满少妇做爰视频| 精品人妻一区二区三区麻豆| 黄频高清免费视频| 国产成人欧美在线观看 | 高潮久久久久久久久久久不卡| 欧美成人午夜精品| 波多野结衣av一区二区av| 亚洲精品国产精品久久久不卡| 电影成人av| 亚洲第一欧美日韩一区二区三区 | 亚洲专区国产一区二区| 在线永久观看黄色视频| 十八禁人妻一区二区| 男女高潮啪啪啪动态图| 亚洲五月色婷婷综合| 热99久久久久精品小说推荐| 亚洲精品久久久久久婷婷小说| 国产精品成人在线| 亚洲av电影在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 建设人人有责人人尽责人人享有的| 国产又色又爽无遮挡免| 亚洲熟女精品中文字幕| 午夜福利视频精品| 亚洲国产毛片av蜜桃av| 少妇裸体淫交视频免费看高清 | 成人国产一区最新在线观看| 91av网站免费观看| 国产激情久久老熟女| 五月天丁香电影| 亚洲欧美色中文字幕在线| 狠狠婷婷综合久久久久久88av| 久久久久视频综合| 老司机影院成人| 日本欧美视频一区| 又大又爽又粗| 50天的宝宝边吃奶边哭怎么回事| 亚洲伊人色综图| 久久女婷五月综合色啪小说| 亚洲第一欧美日韩一区二区三区 | 国产精品久久久久久精品古装| 午夜福利,免费看| 久久av网站| 成人手机av| 两个人看的免费小视频| 精品国产超薄肉色丝袜足j| 亚洲男人天堂网一区| 两性夫妻黄色片| 91老司机精品| 少妇粗大呻吟视频| 91成年电影在线观看| 性少妇av在线| 久久精品亚洲熟妇少妇任你| 另类亚洲欧美激情| 最新的欧美精品一区二区| 正在播放国产对白刺激| 午夜视频精品福利| 岛国毛片在线播放| 99久久精品国产亚洲精品| 咕卡用的链子| 免费少妇av软件| 满18在线观看网站| 91大片在线观看| 精品国产超薄肉色丝袜足j| 老汉色av国产亚洲站长工具| 成人国产av品久久久| 另类精品久久| 久久久国产精品麻豆| 热99re8久久精品国产| 国产免费一区二区三区四区乱码| h视频一区二区三区| 中亚洲国语对白在线视频| 成人影院久久| 18禁观看日本| 日本a在线网址| 一级片'在线观看视频| 亚洲欧美精品综合一区二区三区| 精品少妇黑人巨大在线播放| 黄色片一级片一级黄色片| 王馨瑶露胸无遮挡在线观看| 在线观看人妻少妇| 国产精品 国内视频| 久久热在线av| 三上悠亚av全集在线观看| 成人亚洲精品一区在线观看| 老司机福利观看| 精品福利观看| 亚洲精品美女久久av网站| 人人妻人人澡人人看| 亚洲欧美一区二区三区久久| 国产av一区二区精品久久| 欧美日韩成人在线一区二区| 国产成人精品久久二区二区免费| 亚洲国产看品久久| 亚洲美女黄色视频免费看| 窝窝影院91人妻| 乱人伦中国视频| 少妇裸体淫交视频免费看高清 | 国产精品久久久av美女十八| a在线观看视频网站| 50天的宝宝边吃奶边哭怎么回事| av在线老鸭窝| 国产一区二区三区av在线| 一区二区三区四区激情视频| 亚洲国产毛片av蜜桃av| 在线观看免费日韩欧美大片| 国产成人精品久久二区二区免费| 80岁老熟妇乱子伦牲交| 三级毛片av免费| 无限看片的www在线观看| 美女中出高潮动态图| 亚洲av国产av综合av卡| 亚洲中文日韩欧美视频| av欧美777| 桃红色精品国产亚洲av| 咕卡用的链子| 午夜福利,免费看| 久久精品久久久久久噜噜老黄| 午夜91福利影院| 日韩电影二区| 欧美日韩福利视频一区二区| 亚洲国产精品成人久久小说| 国产三级黄色录像| 多毛熟女@视频| 999久久久精品免费观看国产| 成人免费观看视频高清| 国产激情久久老熟女| 国产无遮挡羞羞视频在线观看| 亚洲av电影在线进入| 大码成人一级视频| 黄色 视频免费看| 亚洲性夜色夜夜综合| 高清av免费在线| 亚洲三区欧美一区| 国内毛片毛片毛片毛片毛片| 777米奇影视久久| 精品少妇内射三级| 成人国产一区最新在线观看| 欧美xxⅹ黑人| 脱女人内裤的视频| 久久狼人影院| 青草久久国产| 亚洲va日本ⅴa欧美va伊人久久 | 欧美少妇被猛烈插入视频| 亚洲国产欧美日韩在线播放| 亚洲一码二码三码区别大吗| 国产av国产精品国产| 少妇精品久久久久久久| 亚洲天堂av无毛| 欧美人与性动交α欧美软件| 一二三四社区在线视频社区8| 老司机午夜十八禁免费视频| 国产日韩欧美视频二区| 国产av精品麻豆| 国产片内射在线| 纯流量卡能插随身wifi吗| 精品免费久久久久久久清纯 | 免费久久久久久久精品成人欧美视频| 欧美日韩黄片免| 黄色视频,在线免费观看| 欧美中文综合在线视频| 精品人妻1区二区| 亚洲一区二区三区欧美精品| 亚洲精品国产av成人精品| 少妇被粗大的猛进出69影院| 日韩有码中文字幕| 日韩一区二区三区影片| 国产一区二区三区在线臀色熟女 | 国产一区二区三区综合在线观看| 一级片'在线观看视频| 国产一级毛片在线| 国产免费福利视频在线观看| 精品少妇久久久久久888优播| 精品人妻1区二区| 丝袜人妻中文字幕| 黑人操中国人逼视频| av视频免费观看在线观看| 亚洲欧美成人综合另类久久久| www日本在线高清视频| 欧美黑人欧美精品刺激| 欧美另类亚洲清纯唯美| 午夜免费成人在线视频| 国产熟女午夜一区二区三区| 国产精品免费视频内射| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 国产精品99久久99久久久不卡| 国产亚洲欧美在线一区二区| 国产一区二区 视频在线| tube8黄色片| 麻豆av在线久日| 王馨瑶露胸无遮挡在线观看| 国产一区二区三区综合在线观看| 美女扒开内裤让男人捅视频| 精品人妻一区二区三区麻豆| av福利片在线| 亚洲三区欧美一区| 伊人亚洲综合成人网| 黄色毛片三级朝国网站| 国产91精品成人一区二区三区 | 色婷婷久久久亚洲欧美| 日韩熟女老妇一区二区性免费视频| 免费av中文字幕在线| 自线自在国产av| 免费看十八禁软件| 国产亚洲精品一区二区www | 国产亚洲精品一区二区www | 国产伦理片在线播放av一区| 久久久国产一区二区| 亚洲 国产 在线| 中亚洲国语对白在线视频| 色播在线永久视频| 精品人妻熟女毛片av久久网站| 中文字幕色久视频| 99久久人妻综合| 在线观看免费日韩欧美大片| 不卡av一区二区三区| 日韩,欧美,国产一区二区三区| 中文精品一卡2卡3卡4更新| 飞空精品影院首页| 人成视频在线观看免费观看| 91九色精品人成在线观看| 免费在线观看完整版高清| 首页视频小说图片口味搜索| 国产精品影院久久| 一进一出抽搐动态| 婷婷成人精品国产| 亚洲一区二区三区欧美精品| 色老头精品视频在线观看| videosex国产| av网站在线播放免费| 水蜜桃什么品种好| 大型av网站在线播放| 欧美另类亚洲清纯唯美| 国产一级毛片在线| 久久天躁狠狠躁夜夜2o2o| 91九色精品人成在线观看| 一二三四社区在线视频社区8| 亚洲欧洲日产国产| 日韩电影二区| 色综合欧美亚洲国产小说| www.自偷自拍.com| 在线观看免费午夜福利视频| 999久久久精品免费观看国产| 欧美精品人与动牲交sv欧美| 一本综合久久免费| 欧美日本中文国产一区发布| 国产精品一区二区免费欧美 | 日日夜夜操网爽| 亚洲视频免费观看视频| 少妇的丰满在线观看| 国产高清视频在线播放一区 | 999久久久国产精品视频| 欧美亚洲 丝袜 人妻 在线| 99国产精品一区二区三区| 性高湖久久久久久久久免费观看| 日日夜夜操网爽| 久久中文看片网| 777米奇影视久久| 亚洲伊人色综图| 久9热在线精品视频| 又紧又爽又黄一区二区| 啪啪无遮挡十八禁网站| 国产精品久久久久久精品电影小说| 国产精品自产拍在线观看55亚洲 | 国产片内射在线| 国产xxxxx性猛交| 国产一区二区在线观看av| 亚洲国产欧美日韩在线播放| 欧美亚洲 丝袜 人妻 在线| 曰老女人黄片| 99国产精品一区二区蜜桃av | 一区二区三区四区激情视频| 亚洲精品美女久久av网站| 丰满迷人的少妇在线观看| 精品一品国产午夜福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 极品少妇高潮喷水抽搐| 精品一区二区三区四区五区乱码| av线在线观看网站| 国产一区二区在线观看av| 少妇粗大呻吟视频| 99精品久久久久人妻精品| 日韩一区二区三区影片| 嫁个100分男人电影在线观看| 美女国产高潮福利片在线看| 欧美激情高清一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产精品自产拍在线观看55亚洲 | 久久久久久久久免费视频了| 久久中文看片网| 久久久精品94久久精品| 精品人妻在线不人妻| 老司机福利观看| av在线老鸭窝| 久久久久国产一级毛片高清牌| 亚洲精品一区蜜桃| 99热网站在线观看| 女人高潮潮喷娇喘18禁视频| 最新在线观看一区二区三区| 久久久国产欧美日韩av| 日韩中文字幕视频在线看片| 国产精品麻豆人妻色哟哟久久| e午夜精品久久久久久久| 好男人电影高清在线观看| 男人操女人黄网站| 美女视频免费永久观看网站| 99香蕉大伊视频| 老司机靠b影院| 在线观看免费日韩欧美大片| 不卡一级毛片| 精品熟女少妇八av免费久了| 老熟妇仑乱视频hdxx| 后天国语完整版免费观看| 热re99久久精品国产66热6| a在线观看视频网站| 亚洲五月色婷婷综合| 亚洲精华国产精华精| 色老头精品视频在线观看| 最近最新中文字幕大全免费视频| 亚洲国产精品成人久久小说| 久久天躁狠狠躁夜夜2o2o| 水蜜桃什么品种好| 免费看十八禁软件| 久久久精品区二区三区| 欧美日韩视频精品一区| 久久久久久久国产电影| 极品人妻少妇av视频| 亚洲av成人不卡在线观看播放网 | 成年动漫av网址| 欧美日韩亚洲国产一区二区在线观看 | 99久久国产精品久久久| 亚洲av成人不卡在线观看播放网 | 一区福利在线观看| 亚洲欧洲日产国产| 国产一卡二卡三卡精品| 黄色a级毛片大全视频| 亚洲,欧美精品.| 久久ye,这里只有精品| 日韩视频在线欧美| 午夜久久久在线观看| avwww免费| 人妻久久中文字幕网| 免费少妇av软件| 欧美一级毛片孕妇| 精品视频人人做人人爽| 欧美国产精品一级二级三级| 91麻豆精品激情在线观看国产 | 国产在线免费精品| 丝袜美足系列| videosex国产| 女人爽到高潮嗷嗷叫在线视频| 美女脱内裤让男人舔精品视频| 一区二区三区乱码不卡18| 日韩欧美免费精品| 在线十欧美十亚洲十日本专区| 成在线人永久免费视频| 国产成+人综合+亚洲专区| 免费女性裸体啪啪无遮挡网站| 无限看片的www在线观看| 欧美xxⅹ黑人| 亚洲欧美精品综合一区二区三区| 91国产中文字幕| 亚洲精品中文字幕在线视频| 欧美日韩一级在线毛片| 97在线人人人人妻| 可以免费在线观看a视频的电影网站| tube8黄色片| 老司机福利观看| 欧美黑人精品巨大| 久久综合国产亚洲精品| 深夜精品福利| 中文字幕av电影在线播放| 岛国毛片在线播放| 一区二区三区乱码不卡18| 黄频高清免费视频| av免费在线观看网站| 中国美女看黄片| 日韩一区二区三区影片| 精品熟女少妇八av免费久了| 涩涩av久久男人的天堂| 美女扒开内裤让男人捅视频| a级片在线免费高清观看视频| 9色porny在线观看| 国产精品.久久久| 丁香六月欧美| 国产亚洲精品久久久久5区| 久久久久网色| 国产成人影院久久av| 韩国高清视频一区二区三区| bbb黄色大片| 久久精品国产a三级三级三级| 在线av久久热| 考比视频在线观看| 久久99热这里只频精品6学生| 国产欧美亚洲国产| 日韩免费高清中文字幕av| 久久99一区二区三区| 97精品久久久久久久久久精品| 久久性视频一级片| 久久久久久久国产电影| 亚洲精品乱久久久久久| 欧美黑人精品巨大| 亚洲欧美日韩高清在线视频 | 国产麻豆69| 国产精品麻豆人妻色哟哟久久| 亚洲国产中文字幕在线视频| 欧美黄色片欧美黄色片| 国产免费视频播放在线视频| 香蕉丝袜av| 亚洲一区中文字幕在线| 成人手机av| 天天躁夜夜躁狠狠躁躁| 亚洲第一欧美日韩一区二区三区 | 五月开心婷婷网| 亚洲欧美一区二区三区久久| 亚洲av欧美aⅴ国产| 国产日韩欧美视频二区| 18禁黄网站禁片午夜丰满| 丝瓜视频免费看黄片| 亚洲精品乱久久久久久| 亚洲视频免费观看视频| 欧美性长视频在线观看| 中亚洲国语对白在线视频| 国产无遮挡羞羞视频在线观看| 国产国语露脸激情在线看| 波多野结衣av一区二区av| 亚洲欧洲日产国产| 日韩,欧美,国产一区二区三区| netflix在线观看网站| 亚洲一码二码三码区别大吗| 动漫黄色视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 巨乳人妻的诱惑在线观看| 久久ye,这里只有精品| 黄色视频不卡| 夜夜骑夜夜射夜夜干| 精品福利永久在线观看| 丝袜在线中文字幕| 亚洲精品乱久久久久久| 1024香蕉在线观看| 亚洲成人手机| 99国产精品一区二区蜜桃av | 国产精品免费视频内射| 最黄视频免费看| 女人高潮潮喷娇喘18禁视频| 国产精品 欧美亚洲| 十八禁网站网址无遮挡| 我要看黄色一级片免费的| 日韩人妻精品一区2区三区| 国产成人欧美| 亚洲欧美成人综合另类久久久| 天天影视国产精品| 久久综合国产亚洲精品| 国产精品av久久久久免费| 手机成人av网站| 欧美日韩视频精品一区| 天天躁狠狠躁夜夜躁狠狠躁| 欧美精品av麻豆av| 午夜福利在线免费观看网站| 国产精品一区二区精品视频观看| 亚洲熟女精品中文字幕| 国产一区有黄有色的免费视频| 欧美精品av麻豆av| 午夜福利在线免费观看网站| av免费在线观看网站| 99精品久久久久人妻精品| 少妇精品久久久久久久| 亚洲情色 制服丝袜| 亚洲性夜色夜夜综合| 亚洲精品一卡2卡三卡4卡5卡 | 黄网站色视频无遮挡免费观看| 国产成人一区二区三区免费视频网站| 亚洲精品久久成人aⅴ小说| 久久青草综合色| 少妇裸体淫交视频免费看高清 | 天天操日日干夜夜撸| 国产精品一二三区在线看| 精品一区二区三区四区五区乱码| 两个人免费观看高清视频| 中文字幕av电影在线播放| 亚洲国产精品一区二区三区在线| 精品国产超薄肉色丝袜足j| 操美女的视频在线观看| 国产一卡二卡三卡精品| 午夜91福利影院| 人成视频在线观看免费观看| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲欧美清纯卡通| 一级,二级,三级黄色视频| av免费在线观看网站| 免费看十八禁软件| 国产精品影院久久| 欧美精品av麻豆av| 国产av精品麻豆| 欧美成人午夜精品| 国产区一区二久久| 黄片小视频在线播放| 搡老熟女国产l中国老女人| 久久天躁狠狠躁夜夜2o2o| 韩国精品一区二区三区| av一本久久久久| 国产精品一区二区在线观看99| 狂野欧美激情性bbbbbb| 国产一区二区在线观看av| 久久毛片免费看一区二区三区| 亚洲成人手机| 最新的欧美精品一区二区| 97在线人人人人妻| 精品第一国产精品| 色94色欧美一区二区| 久久久久国产精品人妻一区二区| bbb黄色大片| 两个人免费观看高清视频| 精品少妇内射三级| 最新的欧美精品一区二区| 精品人妻一区二区三区麻豆| 久久国产精品人妻蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 动漫黄色视频在线观看| 亚洲国产欧美在线一区| 日本黄色日本黄色录像| 国产精品久久久久久精品电影小说| 在线观看免费午夜福利视频| 精品一品国产午夜福利视频| 亚洲中文日韩欧美视频| 亚洲精品久久久久久婷婷小说| 精品亚洲成a人片在线观看| 国产高清国产精品国产三级| 精品国产乱码久久久久久男人| 999久久久精品免费观看国产| 欧美一级毛片孕妇| 国产精品九九99| 97人妻天天添夜夜摸| 日本一区二区免费在线视频| 国产精品国产av在线观看| 咕卡用的链子| 制服诱惑二区| 在线观看免费视频网站a站| 亚洲中文日韩欧美视频| 欧美亚洲 丝袜 人妻 在线| 少妇精品久久久久久久| 一区二区三区四区激情视频| a级毛片在线看网站| 精品人妻1区二区| 久久综合国产亚洲精品| 国产成人精品久久二区二区91| 亚洲 国产 在线| 国产高清国产精品国产三级| 免费女性裸体啪啪无遮挡网站| 飞空精品影院首页| 欧美国产精品一级二级三级| 久久女婷五月综合色啪小说| 久久久精品区二区三区| 久久久国产一区二区| 久久久国产精品麻豆| 午夜两性在线视频| 国产男人的电影天堂91| 成人免费观看视频高清| 亚洲三区欧美一区| 成人三级做爰电影| 精品久久久久久久毛片微露脸 | 国产精品 国内视频| 极品少妇高潮喷水抽搐| 久久99热这里只频精品6学生| 亚洲国产欧美一区二区综合| 国产野战对白在线观看| 黄色怎么调成土黄色| 欧美精品人与动牲交sv欧美| av天堂在线播放| 黄色怎么调成土黄色| av超薄肉色丝袜交足视频| 91成人精品电影| 亚洲色图综合在线观看| 亚洲精品日韩在线中文字幕| 不卡一级毛片| 夫妻午夜视频| 少妇裸体淫交视频免费看高清 | 一区二区三区精品91| 男女午夜视频在线观看| 色视频在线一区二区三区| 天天操日日干夜夜撸| 新久久久久国产一级毛片| 亚洲 欧美一区二区三区| 这个男人来自地球电影免费观看| 国产在线观看jvid| 大片免费播放器 马上看|