• <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ū)消除方法
    国产成人av教育| av欧美777| 久久久精品国产亚洲av高清涩受| 伦理电影免费视频| 国产男人的电影天堂91| 亚洲avbb在线观看| 国产男人的电影天堂91| 亚洲情色 制服丝袜| 亚洲精品一区蜜桃| 亚洲一区中文字幕在线| 搡老熟女国产l中国老女人| 国产精品熟女久久久久浪| 9191精品国产免费久久| 久久精品国产综合久久久| 国产1区2区3区精品| 丝袜美足系列| 热re99久久国产66热| 亚洲精品国产av成人精品| 五月天丁香电影| 亚洲精品第二区| 欧美亚洲日本最大视频资源| 精品人妻一区二区三区麻豆| 美女高潮喷水抽搐中文字幕| avwww免费| 久久 成人 亚洲| 久9热在线精品视频| 亚洲第一欧美日韩一区二区三区 | 免费在线观看视频国产中文字幕亚洲 | 亚洲精品国产一区二区精华液| 亚洲情色 制服丝袜| av国产精品久久久久影院| 亚洲精品国产av成人精品| 美女视频免费永久观看网站| 亚洲国产毛片av蜜桃av| 久久久久精品人妻al黑| 欧美精品一区二区大全| 夜夜夜夜夜久久久久| 色视频在线一区二区三区| av不卡在线播放| 精品福利观看| 亚洲 国产 在线| 宅男免费午夜| 免费观看人在逋| 视频在线观看一区二区三区| 精品国产一区二区久久| 十分钟在线观看高清视频www| 久久国产精品人妻蜜桃| av在线播放精品| 美女脱内裤让男人舔精品视频| 亚洲伊人色综图| 涩涩av久久男人的天堂| 久久久精品免费免费高清| 亚洲中文av在线| 十八禁人妻一区二区| 爱豆传媒免费全集在线观看| 亚洲精品久久午夜乱码| 亚洲av成人一区二区三| 亚洲中文日韩欧美视频| 美女国产高潮福利片在线看| 啦啦啦啦在线视频资源| 欧美在线黄色| 亚洲中文日韩欧美视频| 久久国产亚洲av麻豆专区| 中文字幕高清在线视频| 老司机深夜福利视频在线观看 | 一个人免费在线观看的高清视频 | 中文字幕色久视频| √禁漫天堂资源中文www| 青青草视频在线视频观看| 在线观看舔阴道视频| 99国产精品一区二区三区| 亚洲精品第二区| 18禁观看日本| 亚洲国产看品久久| 91麻豆精品激情在线观看国产 | 国产成人免费观看mmmm| 69av精品久久久久久 | 成人手机av| 久9热在线精品视频| 熟女少妇亚洲综合色aaa.| 国产1区2区3区精品| 色老头精品视频在线观看| 女人被躁到高潮嗷嗷叫费观| 久久久欧美国产精品| 我的亚洲天堂| 交换朋友夫妻互换小说| 国产精品偷伦视频观看了| 国产成人免费观看mmmm| 极品人妻少妇av视频| 午夜免费成人在线视频| 新久久久久国产一级毛片| 一个人免费看片子| 欧美日韩国产mv在线观看视频| 亚洲 欧美一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 黑人巨大精品欧美一区二区mp4| 午夜精品国产一区二区电影| 亚洲精品一区蜜桃| 国产免费av片在线观看野外av| av有码第一页| 69av精品久久久久久 | 纵有疾风起免费观看全集完整版| 操出白浆在线播放| 丰满人妻熟妇乱又伦精品不卡| 啦啦啦中文免费视频观看日本| 黄色 视频免费看| 国产精品二区激情视频| 老司机午夜十八禁免费视频| 亚洲精品国产av成人精品| 少妇精品久久久久久久| 麻豆av在线久日| 亚洲精品国产色婷婷电影| 亚洲成人免费电影在线观看| 亚洲国产日韩一区二区| 十分钟在线观看高清视频www| 欧美 亚洲 国产 日韩一| 一个人免费看片子| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产色婷婷电影| 一区二区三区激情视频| 90打野战视频偷拍视频| 老司机福利观看| 曰老女人黄片| 一二三四社区在线视频社区8| 久9热在线精品视频| 中文字幕高清在线视频| 一区二区日韩欧美中文字幕| 最黄视频免费看| 男男h啪啪无遮挡| 涩涩av久久男人的天堂| 亚洲av国产av综合av卡| 大片免费播放器 马上看| 男女午夜视频在线观看| 69av精品久久久久久 | 777久久人妻少妇嫩草av网站| 考比视频在线观看| bbb黄色大片| 咕卡用的链子| 国产精品久久久人人做人人爽| 国产精品成人在线| 在线观看舔阴道视频| 久久国产精品男人的天堂亚洲| av网站免费在线观看视频| 中文字幕最新亚洲高清| 亚洲少妇的诱惑av| 亚洲一码二码三码区别大吗| 亚洲欧洲精品一区二区精品久久久| 在线观看www视频免费| 99国产精品免费福利视频| 一进一出抽搐动态| 日本撒尿小便嘘嘘汇集6| 老司机亚洲免费影院| 久久久欧美国产精品| 午夜精品久久久久久毛片777| 99久久综合免费| 精品亚洲成国产av| 纵有疾风起免费观看全集完整版| 人人澡人人妻人| 欧美+亚洲+日韩+国产| 精品国产乱码久久久久久男人| 午夜久久久在线观看| 亚洲熟女毛片儿| 亚洲中文av在线| 精品久久蜜臀av无| 欧美亚洲日本最大视频资源| 久久99一区二区三区| 国产精品一区二区在线不卡| av在线app专区| 亚洲性夜色夜夜综合| 午夜激情久久久久久久| 少妇粗大呻吟视频| 丝袜美腿诱惑在线| 中国国产av一级| 伊人久久大香线蕉亚洲五| 91av网站免费观看| 黄色怎么调成土黄色| 精品亚洲成a人片在线观看| 色精品久久人妻99蜜桃| 国产成人a∨麻豆精品| 男男h啪啪无遮挡| 欧美精品亚洲一区二区| 日韩欧美国产一区二区入口| 汤姆久久久久久久影院中文字幕| 亚洲国产av新网站| 久久99热这里只频精品6学生| 国产精品国产三级国产专区5o| 久热这里只有精品99| 国产免费福利视频在线观看| 99香蕉大伊视频| 青春草亚洲视频在线观看| 在线观看舔阴道视频| 一本一本久久a久久精品综合妖精| 亚洲精品在线美女| 亚洲,欧美精品.| 亚洲美女黄色视频免费看| 欧美国产精品一级二级三级| 亚洲精品在线美女| 免费在线观看完整版高清| 又大又爽又粗| 一级毛片精品| 老司机亚洲免费影院| 国产精品香港三级国产av潘金莲| 99国产综合亚洲精品| 成人三级做爰电影| 一本久久精品| 国产精品麻豆人妻色哟哟久久| 伊人亚洲综合成人网| 午夜成年电影在线免费观看| 一区二区三区激情视频| 啦啦啦 在线观看视频| 日本wwww免费看| 午夜日韩欧美国产| 一个人免费看片子| 日韩,欧美,国产一区二区三区| 国产国语露脸激情在线看| 窝窝影院91人妻| 在线观看人妻少妇| 交换朋友夫妻互换小说| 多毛熟女@视频| 美女视频免费永久观看网站| 国产精品.久久久| 欧美亚洲 丝袜 人妻 在线| 午夜免费鲁丝| av免费在线观看网站| 美女中出高潮动态图| 69av精品久久久久久 | 国产成人影院久久av| 精品一区在线观看国产| 水蜜桃什么品种好| 97在线人人人人妻| 人人妻人人澡人人爽人人夜夜| 一区二区av电影网| 亚洲色图综合在线观看| 9热在线视频观看99| 亚洲少妇的诱惑av| 男女免费视频国产| 久久av网站| 美女高潮喷水抽搐中文字幕| 久久久国产精品麻豆| 桃花免费在线播放| 国产成人欧美在线观看 | 日本av免费视频播放| 2018国产大陆天天弄谢| √禁漫天堂资源中文www| 免费一级毛片在线播放高清视频 | 久久精品国产亚洲av香蕉五月 | 久久国产精品大桥未久av| 午夜福利免费观看在线| 人妻一区二区av| 肉色欧美久久久久久久蜜桃| 亚洲欧美日韩另类电影网站| 亚洲精品一区蜜桃| 国产一卡二卡三卡精品| 亚洲精品日韩在线中文字幕| av片东京热男人的天堂| 男女国产视频网站| 欧美中文综合在线视频| 国产精品久久久久久精品电影小说| www.自偷自拍.com| 国产在线观看jvid| 亚洲男人天堂网一区| cao死你这个sao货| 中文字幕最新亚洲高清| 悠悠久久av| 国产精品一区二区在线观看99| 别揉我奶头~嗯~啊~动态视频 | 欧美激情 高清一区二区三区| 欧美精品一区二区免费开放| 高潮久久久久久久久久久不卡| 成人av一区二区三区在线看 | 在线观看舔阴道视频| 亚洲第一av免费看| 大型av网站在线播放| 久久午夜综合久久蜜桃| 啦啦啦啦在线视频资源| 成人免费观看视频高清| 成人国语在线视频| 久久久久网色| 亚洲视频免费观看视频| 免费观看a级毛片全部| 国产成人av激情在线播放| 青春草亚洲视频在线观看| 日韩大码丰满熟妇| 亚洲综合色网址| 欧美乱码精品一区二区三区| 18禁观看日本| 人人妻,人人澡人人爽秒播| 久久中文字幕一级| 色综合欧美亚洲国产小说| 曰老女人黄片| 久久中文看片网| 免费在线观看黄色视频的| 午夜福利影视在线免费观看| 一区福利在线观看| 中文欧美无线码| 在线观看免费高清a一片| 国产成人系列免费观看| 亚洲国产精品一区二区三区在线| xxxhd国产人妻xxx| 亚洲成av片中文字幕在线观看| 午夜成年电影在线免费观看| 久久精品aⅴ一区二区三区四区| 最新在线观看一区二区三区| 午夜免费成人在线视频| 日本撒尿小便嘘嘘汇集6| 亚洲国产成人一精品久久久| 午夜精品久久久久久毛片777| 亚洲成人免费电影在线观看| 十八禁人妻一区二区| a 毛片基地| 成年美女黄网站色视频大全免费| 超碰成人久久| 亚洲精品中文字幕一二三四区 | 女性被躁到高潮视频| 国产精品自产拍在线观看55亚洲 | 久久午夜综合久久蜜桃| 欧美午夜高清在线| 欧美另类亚洲清纯唯美| 一区福利在线观看| 欧美中文综合在线视频| 80岁老熟妇乱子伦牲交| 国产高清videossex| av欧美777| 精品久久蜜臀av无| 国产xxxxx性猛交| 国产97色在线日韩免费| 精品亚洲成国产av| 夜夜夜夜夜久久久久| 亚洲av成人不卡在线观看播放网 | 亚洲国产毛片av蜜桃av| 成人国语在线视频| 免费少妇av软件| 成人国语在线视频| 免费少妇av软件| 首页视频小说图片口味搜索| 久久久国产欧美日韩av| 国产97色在线日韩免费| 欧美av亚洲av综合av国产av| 少妇的丰满在线观看| 人人妻,人人澡人人爽秒播| 在线观看免费午夜福利视频| 久久人人97超碰香蕉20202| 国产欧美日韩一区二区精品| 韩国高清视频一区二区三区| 亚洲熟女精品中文字幕| 激情视频va一区二区三区| 这个男人来自地球电影免费观看| 九色亚洲精品在线播放| 久久免费观看电影| 国产视频一区二区在线看| a级毛片在线看网站| 老司机在亚洲福利影院| 男女边摸边吃奶| 下体分泌物呈黄色| 亚洲精品自拍成人| 日韩一区二区三区影片| 亚洲美女黄色视频免费看| 一区二区三区激情视频| 亚洲天堂av无毛| 国产无遮挡羞羞视频在线观看| 欧美亚洲 丝袜 人妻 在线| 国产亚洲av高清不卡| 日韩人妻精品一区2区三区| 欧美日韩精品网址| 1024视频免费在线观看| 亚洲三区欧美一区| tocl精华| 久久人人97超碰香蕉20202| 18禁黄网站禁片午夜丰满| 纯流量卡能插随身wifi吗| 69精品国产乱码久久久| 欧美人与性动交α欧美软件| 久久久久久久精品精品| 免费在线观看影片大全网站| 天天操日日干夜夜撸| 亚洲精品久久午夜乱码| 99国产极品粉嫩在线观看| 又紧又爽又黄一区二区| 亚洲国产欧美日韩在线播放| 51午夜福利影视在线观看| 久久青草综合色| 午夜免费成人在线视频| 嫁个100分男人电影在线观看| 欧美中文综合在线视频| 午夜影院在线不卡| 久久这里只有精品19| av免费在线观看网站| 啦啦啦视频在线资源免费观看| 国产精品国产三级国产专区5o| 99精国产麻豆久久婷婷| 男女国产视频网站| 纯流量卡能插随身wifi吗| 岛国毛片在线播放| 一本色道久久久久久精品综合| 美国免费a级毛片| 美女中出高潮动态图| 欧美日韩成人在线一区二区| 亚洲成人手机| 男人添女人高潮全过程视频| 超色免费av| 我要看黄色一级片免费的| www.av在线官网国产| 青春草亚洲视频在线观看| 一区福利在线观看| 一进一出抽搐动态| 九色亚洲精品在线播放| 亚洲七黄色美女视频| 女人爽到高潮嗷嗷叫在线视频| 成年人免费黄色播放视频| 久久午夜综合久久蜜桃| 热99久久久久精品小说推荐| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久国产电影| 亚洲精品日韩在线中文字幕| 一本综合久久免费| 国产三级黄色录像| 淫妇啪啪啪对白视频 | 色老头精品视频在线观看| 国产主播在线观看一区二区| 夜夜夜夜夜久久久久| 麻豆乱淫一区二区| 少妇裸体淫交视频免费看高清 | 少妇粗大呻吟视频| 国产精品久久久久久人妻精品电影 | 99精品久久久久人妻精品| 97精品久久久久久久久久精品| 亚洲欧美一区二区三区黑人| 精品国产乱子伦一区二区三区 | 人人妻人人添人人爽欧美一区卜| 国产不卡av网站在线观看| www.精华液| 美女中出高潮动态图| 十分钟在线观看高清视频www| 人人妻,人人澡人人爽秒播| 别揉我奶头~嗯~啊~动态视频 | 精品熟女少妇八av免费久了| 国产亚洲一区二区精品| 最近最新中文字幕大全免费视频| 国产在视频线精品| 性高湖久久久久久久久免费观看| 精品少妇久久久久久888优播| 18禁国产床啪视频网站| 亚洲欧洲日产国产| 一个人免费看片子| 免费不卡黄色视频| 成年av动漫网址| 中文字幕人妻熟女乱码| 超碰97精品在线观看| 桃红色精品国产亚洲av| 少妇 在线观看| 亚洲精品自拍成人| 午夜激情av网站| 国产福利在线免费观看视频| 国产免费一区二区三区四区乱码| av有码第一页| 精品国内亚洲2022精品成人 | 亚洲综合色网址| 纯流量卡能插随身wifi吗| 国产欧美日韩综合在线一区二区| 亚洲成人免费av在线播放| 国产成人啪精品午夜网站| a级片在线免费高清观看视频| av在线播放精品| 国产男女超爽视频在线观看| 一本综合久久免费| 久久精品成人免费网站| 亚洲全国av大片| 各种免费的搞黄视频| h视频一区二区三区| 黑人欧美特级aaaaaa片| 黑丝袜美女国产一区| www日本在线高清视频| 久久久久国内视频| 成人国产一区最新在线观看| 久久久久视频综合| 亚洲,欧美精品.| 国产免费视频播放在线视频| 高清黄色对白视频在线免费看| 国产精品一二三区在线看| 国产成人av教育| 国产xxxxx性猛交| 黑人欧美特级aaaaaa片| 俄罗斯特黄特色一大片| 色老头精品视频在线观看| 成人手机av| 日韩欧美国产一区二区入口| 国产男女内射视频| 啦啦啦 在线观看视频| 我要看黄色一级片免费的| 又大又爽又粗| 久久人人97超碰香蕉20202| 国产男人的电影天堂91| www.999成人在线观看| 国产免费一区二区三区四区乱码| av超薄肉色丝袜交足视频| 飞空精品影院首页| 老司机午夜福利在线观看视频 | 久久国产精品大桥未久av| 老司机在亚洲福利影院| 桃红色精品国产亚洲av| 女人被躁到高潮嗷嗷叫费观| 中文字幕色久视频| 国产精品二区激情视频| 日韩免费高清中文字幕av| 男女边摸边吃奶| 男人舔女人的私密视频| 高清av免费在线| 高清视频免费观看一区二区| 又黄又粗又硬又大视频| 一区二区三区精品91| 99国产极品粉嫩在线观看| 一级毛片电影观看| 国产欧美亚洲国产| 曰老女人黄片| 一区在线观看完整版| 国产在线观看jvid| 精品少妇久久久久久888优播| 天天添夜夜摸| 国产精品一区二区精品视频观看| 最黄视频免费看| 美女午夜性视频免费| 亚洲精品久久午夜乱码| 肉色欧美久久久久久久蜜桃| 国产无遮挡羞羞视频在线观看| 首页视频小说图片口味搜索| 深夜精品福利| 两性夫妻黄色片| 日韩三级视频一区二区三区| 一区二区日韩欧美中文字幕| 亚洲精品乱久久久久久| 黄色怎么调成土黄色| 欧美老熟妇乱子伦牲交| 色视频在线一区二区三区| 日韩三级视频一区二区三区| 高清欧美精品videossex| 久久天躁狠狠躁夜夜2o2o| 国产精品1区2区在线观看. | 免费一级毛片在线播放高清视频 | 亚洲综合色网址| 99热全是精品| 麻豆国产av国片精品| 精品人妻一区二区三区麻豆| √禁漫天堂资源中文www| 一二三四社区在线视频社区8| 9191精品国产免费久久| 欧美精品一区二区大全| 欧美日韩亚洲高清精品| 老司机午夜福利在线观看视频 | 国产av精品麻豆| 日韩 亚洲 欧美在线| 免费高清在线观看视频在线观看| 男人舔女人的私密视频| √禁漫天堂资源中文www| 岛国毛片在线播放| a级毛片黄视频| 日本av免费视频播放| 国产深夜福利视频在线观看| √禁漫天堂资源中文www| 美女高潮到喷水免费观看| 免费日韩欧美在线观看| 视频在线观看一区二区三区| 国产精品影院久久| 正在播放国产对白刺激| 国产男女超爽视频在线观看| 欧美成人午夜精品| 午夜老司机福利片| 亚洲第一av免费看| 黄片播放在线免费| 日本a在线网址| 国产无遮挡羞羞视频在线观看| 亚洲精品乱久久久久久| 捣出白浆h1v1| www日本在线高清视频| 亚洲天堂av无毛| 男人添女人高潮全过程视频| 日日爽夜夜爽网站| 日韩人妻精品一区2区三区| 欧美+亚洲+日韩+国产| 国产精品欧美亚洲77777| 国产男女超爽视频在线观看| 国产亚洲一区二区精品| 国产精品亚洲av一区麻豆| 欧美午夜高清在线| 在线观看免费高清a一片| 1024视频免费在线观看| 精品人妻一区二区三区麻豆| 淫妇啪啪啪对白视频 | 丝袜美足系列| 18在线观看网站| 久久久国产欧美日韩av| 淫妇啪啪啪对白视频 | 无限看片的www在线观看| 黄色视频,在线免费观看| 一本大道久久a久久精品| 青草久久国产| 午夜激情av网站| 老司机午夜十八禁免费视频| av在线app专区| 久久99一区二区三区| 黄色 视频免费看| 欧美中文综合在线视频| av有码第一页| 又紧又爽又黄一区二区| 亚洲免费av在线视频| 91精品三级在线观看| 9热在线视频观看99| 别揉我奶头~嗯~啊~动态视频 | 天堂8中文在线网| 99国产精品一区二区蜜桃av | 999久久久精品免费观看国产| 高清av免费在线| 国产免费现黄频在线看| 亚洲三区欧美一区| 欧美一级毛片孕妇| 欧美日韩中文字幕国产精品一区二区三区 |