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

    基于仿真的動(dòng)態(tài)環(huán)境試驗(yàn)設(shè)備系統(tǒng)設(shè)計(jì)

    2017-10-13 16:30:14牛寶良
    裝備環(huán)境工程 2017年8期
    關(guān)鍵詞:振動(dòng)臺(tái)仿真技術(shù)加速度

    牛寶良

    基于仿真的動(dòng)態(tài)環(huán)境試驗(yàn)設(shè)備系統(tǒng)設(shè)計(jì)

    牛寶良

    (中國(guó)工程物理研究院總體工程研究所,四川 綿陽(yáng) 621900)

    目的 介紹計(jì)算機(jī)仿真技術(shù)的特點(diǎn)及其在動(dòng)態(tài)設(shè)備系統(tǒng)設(shè)計(jì)中的應(yīng)用。方法 首先簡(jiǎn)要介紹仿真技術(shù)的發(fā)展和特點(diǎn),其次介紹一些仿真經(jīng)驗(yàn),最后通過(guò)若干實(shí)例介紹Simulink軟件在動(dòng)態(tài)設(shè)備系統(tǒng)中的應(yīng)用。結(jié)果 基于Simulink仿真,可以快速獲得系統(tǒng)運(yùn)行過(guò)程中的各種所關(guān)注的參數(shù)運(yùn)行結(jié)果,為系統(tǒng)性能評(píng)估和優(yōu)化提供可靠的依據(jù)。結(jié)論 計(jì)算機(jī)仿真技術(shù)已經(jīng)日趨成熟,應(yīng)用于動(dòng)態(tài)環(huán)境試驗(yàn)系統(tǒng)設(shè)計(jì)中,可提升設(shè)計(jì)水平。

    仿真;動(dòng)態(tài)系統(tǒng);代數(shù)環(huán);環(huán)境試驗(yàn)

    環(huán)境試驗(yàn)設(shè)備中,振動(dòng)、沖擊、加速度、搖擺試驗(yàn)設(shè)備等都屬于動(dòng)態(tài)試驗(yàn)設(shè)備。動(dòng)態(tài)試驗(yàn)設(shè)備系統(tǒng)的特性一般比較復(fù)雜,多個(gè)參數(shù)相互耦合,且往往是頻率的函數(shù),這給系統(tǒng)的設(shè)計(jì)帶來(lái)困難。早期的設(shè)計(jì)技術(shù)只能采用簡(jiǎn)化的基于解析公式的設(shè)計(jì)計(jì)算,以液壓振動(dòng)臺(tái)的伺服閥為例,計(jì)算輸出流量時(shí)不考慮負(fù)載影響、或者負(fù)載按某個(gè)固定值來(lái)處理,而實(shí)際上,負(fù)載是時(shí)時(shí)刻刻都在變動(dòng)中。伺服閥隨著頻率的增加,流量會(huì)衰減,相位滯后也會(huì)不斷增加,相位滯后對(duì)系統(tǒng)的穩(wěn)定性也會(huì)有顯著影響。這樣一種處處簡(jiǎn)化的設(shè)計(jì)手段,對(duì)系統(tǒng)動(dòng)態(tài)特性的預(yù)估和了解是粗糙的、不全面的。因此也增加系統(tǒng)設(shè)計(jì)的風(fēng)險(xiǎn)。

    仿真技術(shù)是分析動(dòng)態(tài)系統(tǒng)的有效手段,它把各種相關(guān)的動(dòng)態(tài)變量放到一組微分方程中,用模型(電路或者軟件)表達(dá)系統(tǒng)的動(dòng)態(tài)特性,運(yùn)行模型得到各量隨時(shí)間的變化曲線。后期隨著全數(shù)字仿真技術(shù)的發(fā)展,不僅僅是微分方程,各種非線性環(huán)節(jié)或者模塊也都可以進(jìn)到系統(tǒng)模型中,因此逼真程度進(jìn)一步提高。隨著微型計(jì)算機(jī)和仿真軟件的普及,仿真技術(shù)在各行各業(yè)應(yīng)用越來(lái)越多。文中通過(guò)若干實(shí)例介紹仿真技術(shù)在動(dòng)態(tài)環(huán)境試驗(yàn)設(shè)備設(shè)計(jì)中的應(yīng)用。

    1 計(jì)算機(jī)仿真技術(shù)歷程與現(xiàn)狀簡(jiǎn)介

    仿真技術(shù)經(jīng)歷了模擬計(jì)算機(jī)、數(shù)字模擬混合計(jì)算機(jī)、全數(shù)字計(jì)算機(jī)三個(gè)大的階段[1]。1946年,一臺(tái)能夠模擬導(dǎo)彈、飛機(jī)飛行軌跡等物理現(xiàn)象的模擬計(jì)算機(jī)在美國(guó)問(wèn)世。標(biāo)志著模擬仿真技術(shù)的開(kāi)始。此后,涌現(xiàn)了大量的模擬機(jī),20世紀(jì)50年代至60年代初是模擬計(jì)算機(jī)的黃金時(shí)期。盡管模擬仿真計(jì)算機(jī)取得了巨大的成功,但也暴露出數(shù)據(jù)精度不夠、實(shí)現(xiàn)插值函數(shù)困難、不能滿足數(shù)字控制系統(tǒng)的需要等問(wèn)題。20世紀(jì)50年代末,由于導(dǎo)彈技術(shù)發(fā)展的需要,促使了混合仿真計(jì)算機(jī)的誕生。20世紀(jì)60年代至70年代中葉是混合仿真計(jì)算機(jī)的鼎盛時(shí)期。美國(guó)ADI公司1978年研制出全數(shù)字仿真計(jì)算機(jī)AD10,它的誕生標(biāo)志著在那場(chǎng)混合仿真計(jì)算機(jī)與數(shù)字仿真計(jì)算機(jī)的激烈競(jìng)爭(zhēng)中,全數(shù)字仿真計(jì)算機(jī)終領(lǐng)風(fēng)騷,也標(biāo)志著全數(shù)字仿真新時(shí)代的來(lái)臨。1973年美國(guó)研制成功專(zhuān)用液壓仿真軟件HYDSIM,之后歐美各國(guó)不斷發(fā)展完善液壓仿真軟件,1980年代先后推出了DSH(德國(guó))、HASP(英國(guó))等軟件。

    1985年10月15日,我國(guó)第一臺(tái)全數(shù)字仿真計(jì)算機(jī)“銀河仿真I型計(jì)算機(jī)(YH-F1)”研制成功,標(biāo)志著中國(guó)進(jìn)入全數(shù)字仿真時(shí)代。此后又有多個(gè)新型的仿真計(jì)算機(jī)問(wèn)世。

    航空、航天是最先開(kāi)始應(yīng)用仿真技術(shù)的幾個(gè)傳統(tǒng)領(lǐng)域,是仿真計(jì)算機(jī)誕生的搖籃。也只有這樣的一些重要的、復(fù)雜的動(dòng)態(tài)系統(tǒng)設(shè)計(jì)才用得起。

    仿真技術(shù)的普遍應(yīng)用則是在微型計(jì)算機(jī)普及和仿真軟件的圖形化編程之后,標(biāo)志性的節(jié)點(diǎn)是MATLAB/Simulink仿真軟件的推出。MATLAB/Sim- ulink在仿真軟件中具有非常重要的地位??梢哉J(rèn)為Simulink是一個(gè)分水嶺,之前建模是基于語(yǔ)句描述,模塊之間的關(guān)系靠語(yǔ)句定義,不能直觀看到系統(tǒng)的模塊組成和模塊之間的關(guān)系。Simulink提供了圖形化的建模方式,只要進(jìn)行鼠標(biāo)的簡(jiǎn)單拖拉操作就可構(gòu)造出復(fù)雜的仿真模型。它外表以方塊圖形式呈現(xiàn),且采用分層結(jié)構(gòu)。從建模角度講,這既適于自上而下(Top-down)的設(shè)計(jì)流程(概念、功能、系統(tǒng)、子系統(tǒng)、直至器件),又適于自下而上(Bottom-up)逆向程序設(shè)計(jì)。從分析研究角度講,這種Simulink模型不僅能讓用戶知道具體環(huán)節(jié)的動(dòng)態(tài)細(xì)節(jié),而且能讓用戶清晰地了解各器件、各子系統(tǒng)、各系統(tǒng)間的信息交換,掌握各部分之間的交互影響。這給建模帶來(lái)極大的方便,也有利于模型的查錯(cuò)和修改。除了建??旖莘奖阃?,其他功能也很強(qiáng)大,它可以處理的系統(tǒng)包括:線性、非線性系統(tǒng);離散、連續(xù)及混合系統(tǒng);單任務(wù)、多任務(wù)離散事件系統(tǒng)。 Simulink 6.5是里程碑版本,它的加速模式能使運(yùn)算速度提高幾十上百倍。

    除了MATLAB/Simulink外,還有一些仿真軟件,比如AMESim 是當(dāng)今領(lǐng)先的傳動(dòng)系統(tǒng)和液壓/機(jī)械系統(tǒng)建模、仿真及動(dòng)力學(xué)分析軟件。再比如ADAMS,是應(yīng)用較多的一種動(dòng)力學(xué)仿真軟件。比如Modelica,是所謂的多領(lǐng)域建模軟件。它們各有特色,或者對(duì)某領(lǐng)域有特別的優(yōu)勢(shì),或者多學(xué)科綜合建模方面有優(yōu)勢(shì)。

    實(shí)時(shí)仿真系統(tǒng)也是仿真系統(tǒng)的一個(gè)重要分支,不僅要求獲得系統(tǒng)響應(yīng),而且要求響應(yīng)與仿真對(duì)象實(shí)際的響應(yīng)時(shí)間是嚴(yán)格的1:1。當(dāng)具有這樣一種時(shí)間上的1:1特性時(shí),就可以與真實(shí)的硬件系統(tǒng)連接,構(gòu)成所謂的硬件在環(huán)仿真系統(tǒng)。這類(lèi)系統(tǒng)通常都不是普通PC所能勝任的,需要專(zhuān)門(mén)的比較高端的硬件系統(tǒng)。國(guó)內(nèi)外也都有相應(yīng)的產(chǎn)品,國(guó)外產(chǎn)品有dSPACE,RT-LAB等,國(guó)內(nèi)的有國(guó)防科技大學(xué)的KDRTS、西安二炮工程學(xué)院的CHY-Ⅲ、航天三院的海鷹實(shí)時(shí)仿真工作站、西北工業(yè)大學(xué)的并行計(jì)算機(jī)PD-100、北京華力創(chuàng)通科技有限公司的HRT1000、北京廣思特公司的半實(shí)物仿真系統(tǒng)、北京經(jīng)緯恒潤(rùn)科技有限公司的Higale半實(shí)物仿真系統(tǒng)等。

    在環(huán)境試驗(yàn)設(shè)備領(lǐng)域的仿真,應(yīng)該說(shuō)起步不算晚,李坦在1985年就在模擬計(jì)算機(jī)上開(kāi)展了液壓振動(dòng)臺(tái)的仿真[2],但是由于仿真軟硬件等不成熟不方便等原因,并未大量普及。大約在2000年之后,仿真應(yīng)用逐漸多了起來(lái)。2010年之后,特別是近幾年,進(jìn)入了普遍應(yīng)用階段[3—7,9—16]。

    如今MATLAB/Simulink仍然是工程技術(shù)人員喜歡的仿真軟件之一。文中基于MATLAB/Simulink來(lái)介紹若干仿真實(shí)例,展示仿真給環(huán)境試驗(yàn)設(shè)備設(shè)計(jì)帶來(lái)的便利。

    2 仿真技術(shù)應(yīng)用實(shí)例

    2.1 仿真的幾個(gè)基本概念

    2.1.1 連續(xù)與離散

    以單自由度振動(dòng)系統(tǒng)為例,它的運(yùn)動(dòng)方程是

    式中:為質(zhì)量塊的質(zhì)量,kg;為運(yùn)動(dòng)阻尼,N/(m·s-1);為彈簧剛度,N/m;為質(zhì)量塊的位移,m;為的二階導(dǎo)數(shù),表示加速度,m/s2;為的一階導(dǎo)數(shù),表示速度,m/s;為作用在質(zhì)量塊上的力,N。

    對(duì)應(yīng)的拉氏變換的表達(dá)式為:

    式中:是位移的拉氏變換;是力的拉氏變換。

    式(2)對(duì)應(yīng)的模型就是連續(xù)模型,它的Simulink模型可以用連續(xù)積分器1/構(gòu)建,如圖1中的上面一個(gè)模型。把式(2)離散化(有多種離散辦法,應(yīng)用較多的是Tustin變換),或者用連續(xù)積分器替代連續(xù)積分器1/s。這里用替換法,用梯形積分(Simulink還提供了前向積分、后向積分,都屬于矩形積分,精度比梯形積分差)替換連續(xù)積分,所得模型如圖1中下面一個(gè)模型。從模型形式上看,連續(xù)系統(tǒng)用的變量,離散系統(tǒng)用的變量。在內(nèi)部解算上,連續(xù)系統(tǒng)有多種解算方法可選,而且這些算法一般都是變步長(zhǎng)的多步算法,離散系統(tǒng)的算法則是定步長(zhǎng)單步法。

    采用離散系統(tǒng)建模有以下好處:離散系統(tǒng)支持Simulink的加速度模式,仿真速度快;單步法,運(yùn)算量是固定的,仿真耗時(shí)也是固定的;半實(shí)物仿真系統(tǒng)必須是離散的,采用離散建模方便轉(zhuǎn)換到半實(shí)物仿真系統(tǒng);記錄的響應(yīng)數(shù)據(jù)是等間隔的,便于后續(xù)處理。文獻(xiàn)[11]給出了詳細(xì)實(shí)現(xiàn)方法。

    當(dāng)用離散積分器(梯形)替換完連續(xù)積分器1/之后,點(diǎn)擊運(yùn)行,會(huì)提示“The model 'xxx' has an algebraic loop. Accelerator mode is not supported for models with algebraic loops.”,意思是模型有一個(gè)代數(shù)環(huán),加速模式不支持代數(shù)環(huán)。

    2.1.2 代數(shù)環(huán)

    盡管離散系統(tǒng)有上述的好處,但是代數(shù)環(huán)是個(gè)問(wèn)題。所謂的代數(shù)環(huán),就是模型有環(huán)路(有反饋必然有環(huán)路),求解當(dāng)前某個(gè)變量時(shí),又要用到這個(gè)變量。對(duì)于動(dòng)態(tài)系統(tǒng),這幾乎是必然要遇到的。解決代數(shù)環(huán)有幾個(gè)辦法,最簡(jiǎn)單的辦法是加延遲環(huán)節(jié),或者積分器采用前向歐拉法(Forward Euler)。這些方法精度不高,有時(shí)還會(huì)發(fā)散。筆者提出了一個(gè)預(yù)測(cè)法,較好地解決了代數(shù)環(huán)。所謂的預(yù)測(cè)法,就是用前面若干點(diǎn)預(yù)測(cè)當(dāng)前點(diǎn)的值。相應(yīng)的Simulink模型見(jiàn)圖2和圖3。

    圖2 預(yù)測(cè)環(huán)節(jié)

    注意:預(yù)測(cè)環(huán)節(jié)中有三個(gè)延遲器,如果知道初值,比如圖1所示的例子,加速度初值是一個(gè)重力加速度,就可以把它們的初值都置為9.806 65 m/s2;如果不知道初值,就置為0。帶來(lái)的后果就是仿真開(kāi)始的幾個(gè)數(shù)據(jù)點(diǎn)不正確,之后就正確了,基本不影響仿真。

    把這個(gè)模塊(模塊名pred)插入在mass1后面,見(jiàn)圖1下面的一個(gè)模型,代數(shù)環(huán)的問(wèn)題就解決了。

    對(duì)圖1所示的仿真模型,給定一組參數(shù),=10 kg,=39478 N/m,=125.66 N/( m·s-1),可以仿真得到連續(xù)系統(tǒng)、離散系統(tǒng)的仿真結(jié)果,見(jiàn)圖4和圖5??梢钥吹剑介L(zhǎng)合適的情況下,連續(xù)系統(tǒng)、離散系統(tǒng)結(jié)果在采樣時(shí)刻是一致的,離散系統(tǒng)仿真結(jié)果波形是階梯狀的。這個(gè)模型同時(shí)也顯示了Simulink其實(shí)是可以在一個(gè)模型里同時(shí)有連續(xù)系統(tǒng)、離散系統(tǒng)。

    2.1.3 采樣間隔

    當(dāng)采用離散建模時(shí),采樣間隔是一個(gè)很關(guān)鍵的參數(shù),采樣間隔的倒數(shù)是采樣頻率。一般來(lái)說(shuō),采樣頻率與系統(tǒng)中振蕩環(huán)節(jié)(也稱二階環(huán)節(jié))的共振頻率相關(guān)。首先,為了能看清振蕩的響應(yīng),采樣頻率應(yīng)該是振蕩環(huán)節(jié)的共振頻率的10倍以上。其次,為了仿真系統(tǒng)運(yùn)行穩(wěn)定(不發(fā)散),采樣間隔也應(yīng)減小。越小的采樣間隔,一般精度越高,但同時(shí)仿真一次的耗時(shí)越長(zhǎng)。因此,仿真間隔也是要從這兩方面取舍。通??梢韵热?×10-4s或者1×10-5s,不行再調(diào)整。

    2.1.4 上采樣、下采樣、抽取

    一個(gè)物理系統(tǒng),有時(shí)由快慢相差很大的模塊構(gòu)成,這時(shí)如果采用統(tǒng)一的采樣間隔,那就只能確??斓哪K,即采樣間隔取小值。這樣一來(lái),對(duì)于慢模塊來(lái)說(shuō),數(shù)據(jù)量就很大,許多數(shù)據(jù)是多余的,而且容易導(dǎo)致計(jì)算機(jī)內(nèi)存溢出。

    一種辦法就是采用抽取,即在示波器、數(shù)據(jù)保存模塊中,設(shè)置Sampling參數(shù)的Decimation,比如設(shè)置Decimation值為100,表示每運(yùn)算100步,記錄1個(gè)數(shù)據(jù)到示波器、數(shù)組、或者文件中。這種辦法數(shù)據(jù)量可以減下來(lái),但是模型計(jì)算量一點(diǎn)不減少。

    另一種辦法是采用上采樣、下采樣的辦法。即慢速模塊用一個(gè)比較大的步長(zhǎng),而快的模塊用一個(gè)小的步長(zhǎng),慢模塊的輸出進(jìn)到快模塊前用上采樣(Upsample),快模塊的輸出進(jìn)到慢模塊前用下采樣(Downsample)。Upsample、Downsample在Signal Processing Blockset下的Signal Operations目錄下可以找到。下采樣模塊直接使用沒(méi)有問(wèn)題。上采樣模塊不能直接使用,因?yàn)樗遣捎眉?的辦法提高采樣頻率,把一個(gè)連續(xù)的信號(hào)變成了脈沖信號(hào),如圖5的上圖所示。筆者給出了修正的上采樣模塊,如圖3所示,得到的上采樣后的波形如圖5的下圖所示。這樣,慢模塊的運(yùn)算量的減小了很多,整個(gè)系統(tǒng)的仿真速度就能大幅提高。

    總之,適當(dāng)設(shè)置抽取參數(shù)、采用上采樣、下采樣有利于避免內(nèi)存溢出,有利于減小數(shù)據(jù)、方便后續(xù)處理和顯示,有利于提高仿真速度。抽取的原則是信號(hào)的每個(gè)周期有10~100個(gè)數(shù)據(jù)點(diǎn)為佳。點(diǎn)數(shù)再少波形就太粗糙,點(diǎn)數(shù)再多也沒(méi)有太大意義。對(duì)于掃頻信號(hào),只能按最高頻率來(lái)設(shè)置。

    避免內(nèi)存溢出的另一個(gè)辦法是在示波器上勾選(limit data points to last)選項(xiàng),填寫(xiě)數(shù)據(jù)點(diǎn)數(shù),只保留仿真最后的若干數(shù)據(jù)點(diǎn)。把數(shù)據(jù)寫(xiě)入文件既能把全程數(shù)據(jù)保留下來(lái),也不會(huì)導(dǎo)致內(nèi)存溢出。

    2.1.5 單位

    單位其實(shí)很重要。由于歷史原因,各種非國(guó)際單位制的單位仍沒(méi)有完全淘汰,再者,有些單位比較直觀方便,仍在使用。筆者的觀點(diǎn)是方便的或者習(xí)慣的單位可以用,但是進(jìn)入仿真程序時(shí)一律換算為國(guó)際單位制,仿真結(jié)果用于顯示時(shí),可以再換算為日常習(xí)慣的單位。比如加速度,通常習(xí)慣用作單位,仿真時(shí)必須換算為m/s2,流量用L/min,仿真時(shí)必須換算為m3/s。有些弄簡(jiǎn)成繁的單位,比如質(zhì)量用kg·s2/cm,則應(yīng)堅(jiān)決淘汰。單位正確是仿真結(jié)果正確的重要保證。

    2.2 液壓振動(dòng)臺(tái)

    液壓振動(dòng)臺(tái)在裝備環(huán)境試驗(yàn)、汽車(chē)道路模擬試驗(yàn)、地震模擬試驗(yàn)等方面廣泛應(yīng)用。雖然技術(shù)指標(biāo)有所不同,但是主要構(gòu)成、控制原理是相同的。下面以某10 t振動(dòng)臺(tái)為例介紹仿真的應(yīng)用,結(jié)構(gòu)如圖6所示。主要參數(shù):位移為±120 mm,速度為1 m/s,加速度為10,負(fù)載為1 t。

    液壓振動(dòng)臺(tái)設(shè)計(jì)最主要的參數(shù)就是活塞有效面積,傳統(tǒng)的設(shè)計(jì)是基于最佳功率匹配原則[8],或者簡(jiǎn)單的負(fù)載壓強(qiáng)=2/3油源壓強(qiáng)?;诜抡娴脑O(shè)計(jì),可以以2/3原則計(jì)算初值,通過(guò)仿真全面了解在整個(gè)頻帶上的效果,然后再進(jìn)行迭代。液壓臺(tái)的流量也是系統(tǒng)的主要參數(shù)之一。有文獻(xiàn)分析過(guò),共振頻率處最低,超過(guò)共振頻率后,逐步上升。通過(guò)仿真可以了解全頻帶的流量需求。建立的液壓振動(dòng)臺(tái)Simulink模型如圖7所示。

    圖6 液壓振動(dòng)臺(tái)結(jié)構(gòu)

    利用該仿真模型可以仿真研究[9]:理想情況下的性能;活塞摩擦對(duì)波形的影響;基礎(chǔ)質(zhì)量對(duì)波形的影響;控制算法對(duì)性能的影響;伺服閥遮蓋對(duì)波形的影響;負(fù)載振動(dòng)頻率特性;基礎(chǔ)的振動(dòng)響應(yīng)。

    這里僅介紹兩種三狀態(tài)控制仿真情況。一種是電壓輸入,一種是加速度輸入。

    電壓輸入時(shí),希望最大電壓在不同頻率對(duì)應(yīng)不同的量的最大值,低頻10 V對(duì)應(yīng)最大位移,中頻段,10 V對(duì)應(yīng)最大速度,高頻段,10 V對(duì)應(yīng)最大加速度。這樣既能發(fā)揮振動(dòng)臺(tái)的最大能力,也不會(huì)因?yàn)轱柡拖鞑▽?dǎo)致波形失真。優(yōu)化后的效果見(jiàn)圖8。

    另外一種是加速度輸入,給定的輸入信號(hào)與目標(biāo)加速度信號(hào)成比例,控制器優(yōu)化的目標(biāo)是讓響應(yīng)加速度盡可能逼近給定加速度。從傳函來(lái)說(shuō),希望傳遞函數(shù)接近一根平直的直線。優(yōu)化后的某地震波控制如圖9所示。由于低頻段受位移限制,加速度很小,所以采用加速度輸入進(jìn)行地震波試驗(yàn)時(shí),必須先做一次三參量生成,看看位移、速度是否超振動(dòng)臺(tái)極限,確認(rèn)不超極限時(shí),方可進(jìn)行試驗(yàn)。

    2.3 二維搖擺試驗(yàn)機(jī)

    二維搖擺試驗(yàn)機(jī)是中國(guó)工程物理研究院總體工程研究所研制的一套試驗(yàn)系統(tǒng)。在研制階段進(jìn)行了仿真分析,其結(jié)構(gòu)和仿真模型如圖10和11所示。

    系統(tǒng)主要參數(shù):運(yùn)動(dòng)自由度為兩自由度,即可同時(shí)進(jìn)行橫搖和縱搖;最大負(fù)載質(zhì)量不小于1 000 kg;負(fù)載高度極限不低于2 200 mm;臺(tái)面尺寸不小于 2 200 mm×2 200 mm;負(fù)載質(zhì)心高度極限不低于1 000 mm;傾斜試驗(yàn)?zāi)芰?,最大橫傾角為±45°,最大縱傾角為±30°;搖擺試驗(yàn)?zāi)芰?,最大橫搖角為±60°、橫搖周期為3~14 s,最大縱搖角為±30°、縱搖周期為4~10 s;傾斜搖擺復(fù)合試驗(yàn)?zāi)芰Γㄅ?0°傾斜夾具),最大橫搖角為±60°、橫搖周期為3~14 s,最大縱搖角為±15°、縱搖周期為4~8 s;

    通過(guò)仿真,得到如圖12所示的閥流量、角位移,得到圖13所示的負(fù)載角速度、驅(qū)動(dòng)力矩、機(jī)械功率,以及其他關(guān)注量,比如跟蹤誤差、齒隙對(duì)搖擺角度影響等。對(duì)閥選型、供電功率配置等提供了有力支撐。

    2.4 離心機(jī)

    中國(guó)工程物理研究院總體工程研究所2號(hào)離心機(jī),容量為200t,最大離心加速度為200。系統(tǒng)結(jié)構(gòu)主體見(jiàn)圖14。系統(tǒng)的仿真模型見(jiàn)圖15。細(xì)節(jié)可參考文獻(xiàn)[10—12]。

    通過(guò)仿真,可以確定多大的上升速率可行,或者給定的一個(gè)上升速率,需要的多大電壓、電流。如圖16所示。也可以為新離心機(jī)設(shè)計(jì)提供支撐。

    2.5 游戲機(jī)

    某升降游戲機(jī)的結(jié)構(gòu)如圖17所示。負(fù)載為4 000 kg,最大升程為20 m,通過(guò)滑輪組實(shí)現(xiàn)3倍的位移放大。負(fù)載通過(guò)兩根鋼絲繩提升,相當(dāng)于每根鋼絲繩承受0.5倍載荷。

    圖16 加速度、電流時(shí)域曲線

    系統(tǒng)仿真模型如圖18所示。包括信號(hào)源、控制器、伺服閥、振動(dòng)系統(tǒng)等幾個(gè)模塊。控制器采用PIDF控制算法。伺服閥采用雙出口伺服閥模型,液壓缸采用單出桿模型,可真實(shí)反映上下兩腔的油液壓強(qiáng)及液壓缸的輸出力。振動(dòng)部分用鋼絲繩把液壓缸活塞和負(fù)載連接,鋼絲繩用單邊剛度彈簧表示,只承受拉,不承受壓,壓剛度為0?,F(xiàn)通過(guò)仿真設(shè)計(jì)系統(tǒng)的主要參數(shù),比如可實(shí)現(xiàn)的加速度、液壓伺服閥的流量、液壓缸的缸徑桿徑、油源的壓力、流量、一次升降需要的油液等。

    通過(guò)多次迭代優(yōu)化,確定伺服閥額定流量為800 L/min,油源壓力為21 MPa、流量為600 L/min。負(fù)載運(yùn)動(dòng)的位移、加速度見(jiàn)圖19。

    可以看到,負(fù)載較好跟蹤了給定位移,有三個(gè)時(shí)刻達(dá)到失重(-1加速度),鋼絲繩峰值載荷為81.26 kN。作為游戲機(jī),主要是給游客加速度感覺(jué),這樣的結(jié)果是滿意的。如果是作為工業(yè)控制,需要精確跟隨給定位移,有若干辦法來(lái)提高精度,比如增加活塞面積、改變鋼絲繩剛度、嘗試不同的控制策略和控制參數(shù)等。

    3 對(duì)仿真的再認(rèn)識(shí)

    仿真技術(shù)的作用已經(jīng)無(wú)需質(zhì)疑,對(duì)于復(fù)雜系統(tǒng),通過(guò)仿真可以提高設(shè)計(jì)效率和質(zhì)量,節(jié)約研制經(jīng)費(fèi)。在此,簡(jiǎn)要總結(jié)一下對(duì)仿真的幾點(diǎn)認(rèn)識(shí)。

    1)深化對(duì)系統(tǒng)的認(rèn)識(shí)。仿真模型的建立依賴一定的科學(xué)規(guī)律,而這些科學(xué)規(guī)律往往是理想條件下對(duì)單一問(wèn)題的科學(xué)表述,比如牛頓第二定律=,黏性阻尼力=·,彈性定律=。一個(gè)單自由度系統(tǒng),包含了上述三個(gè)規(guī)律,理想情況下可以得到解析式。如果阻尼力是=·v,或彈簧力=kx,或彈簧是鋼絲繩類(lèi)(受拉有剛度,受壓剛度為0),或很多的自由度,解析表達(dá)式就可能得不到,因而也就不能得到其響應(yīng)和認(rèn)識(shí)其規(guī)律。這些對(duì)于仿真來(lái)說(shuō),則是非常容易的。因此,科學(xué)規(guī)律是建模的基礎(chǔ),而建模與仿真可以獲得對(duì)系統(tǒng)更深入全面的認(rèn)識(shí)。

    2)系統(tǒng)參數(shù)設(shè)計(jì)與優(yōu)化。當(dāng)給定一組設(shè)計(jì)目標(biāo)參數(shù)后,基于建模仿真,就可以得出系統(tǒng)相關(guān)的許多參數(shù),比如液壓振動(dòng)臺(tái)設(shè)計(jì),通過(guò)仿真可以得到基礎(chǔ)質(zhì)量、油源壓力、電源功率等一系列外界供給參數(shù),由此分析系統(tǒng)的建造規(guī)模、經(jīng)費(fèi)等。通過(guò)多輪次迭代,可以進(jìn)行優(yōu)化,以最低的成本獲得最大收益。

    3)控制器設(shè)計(jì)與優(yōu)化。含有控制系統(tǒng)的試驗(yàn)系統(tǒng),其控制系統(tǒng)往往成為系統(tǒng)的核心子系統(tǒng),它的算法決定了系統(tǒng)的穩(wěn)定性(穩(wěn)定余度)、控制精度等。在仿真系統(tǒng)不用擔(dān)心系統(tǒng)不穩(wěn)定導(dǎo)致的設(shè)備損壞;可以施加多種不同的載荷,檢查控制算法的適宜性;算法、控制器參數(shù)、被控對(duì)象等的修改只需修改仿真模型,操作快速。整個(gè)過(guò)程安全、便捷、高效。

    4)故障與異常分析。對(duì)于真實(shí)的系統(tǒng),人為制造某種故障檢查系統(tǒng)的反映是有很大風(fēng)險(xiǎn)的,在仿真系統(tǒng)中,則不存在任何風(fēng)險(xiǎn)??梢园踩⒎奖愕赜^察系統(tǒng)的反應(yīng)、故障聯(lián)鎖機(jī)制是否起效。

    5)人員培訓(xùn)。比較操作復(fù)雜的重要系統(tǒng),比如飛機(jī)駕駛、高鐵駕駛、核電站的操作,采用仿真機(jī)進(jìn)行培訓(xùn),是安全、高效、節(jié)約的培訓(xùn)方式。

    6)分析工程因素的影響。在做系統(tǒng)設(shè)計(jì)時(shí),往往用的理想模型。比如伺服閥,理想情況是零遮蓋。實(shí)際四個(gè)節(jié)流邊的遮蓋量肯定是非零的,且不相同??梢越㈦p出口四邊滑閥模型,人為給定不同的正負(fù)遮蓋量,通過(guò)仿真看看它們對(duì)系統(tǒng)波形失真、控制精度等的影響。再比如摩擦力,分析原理時(shí)往往不考慮摩擦或者給一個(gè)固定的摩擦力值,實(shí)際上與加工間隙或過(guò)盈量、密封件的材質(zhì)、壓縮量等都有關(guān)系,可以人為給一些不同的摩擦力,看看對(duì)系統(tǒng)性能的影響,或找到摩擦力允許的上限值。

    7)仿真的“真”與“不真”。在正確建模的情況下,可以說(shuō)仿真得到的規(guī)律是“真”的,仿真的值是“不真”的。因?yàn)槿绻P驼_,仿真結(jié)果就能反映系統(tǒng)主要運(yùn)行規(guī)律。如果模型不正確,規(guī)律也就不正確了。

    以液壓振動(dòng)臺(tái)為例,用二階振蕩環(huán)節(jié)表示一個(gè)伺服閥時(shí),振動(dòng)波形不失真,與實(shí)際系統(tǒng)相差甚遠(yuǎn)。因?yàn)槟P鸵呀?jīng)忽略了閥壓降這樣一個(gè)顯著的非線性因素。加上閥壓降這個(gè)非線性因素,則振動(dòng)波形失真特性就能顯示出來(lái)。再以摩擦為例,不考慮摩擦,則波形光滑不失真。加上摩擦這個(gè)實(shí)際存在的因素,結(jié)果就與實(shí)際相近,如圖20所示。

    實(shí)際的系統(tǒng)是非常復(fù)雜的,常常不愿或者不能建立反映每個(gè)參數(shù)細(xì)節(jié)的模型。以液壓缸的摩擦為例,具體每個(gè)液壓缸的摩擦力都是不一樣的(每個(gè)液壓缸的密封圈槽尺寸有加工誤差、每個(gè)密封圈有加工誤差、運(yùn)行溫度也會(huì)影響密封圈的壓縮量)。再以油液粘度為例,它是溫度、壓力的函數(shù),建模時(shí)也常常把它當(dāng)作常數(shù),也會(huì)帶來(lái)誤差。伺服閥的頻率特性每一臺(tái)各有不同,即使是同一臺(tái),它的頻率特性與流量幅值還有關(guān)系,小流量頻響范圍寬,大流量頻響范圍小等。這些因素導(dǎo)致仿真輸入?yún)?shù)與實(shí)際存在差異,因此仿真的值是“不真”的,但是這個(gè)所謂的不真,是有限度的,不能因此而否定仿真。

    總之仿真可以反映系統(tǒng)的主要、基本的規(guī)律,得到的結(jié)果具體值與實(shí)際系統(tǒng)會(huì)有差異,認(rèn)識(shí)這點(diǎn),才能正確看待仿真技術(shù)。

    4 結(jié)語(yǔ)

    文中簡(jiǎn)略回顧了仿真技術(shù)的仿真歷程,給出了一些經(jīng)驗(yàn)和若干仿真在環(huán)境試驗(yàn)設(shè)備方面的應(yīng)用實(shí)例,總結(jié)了仿真技術(shù)的認(rèn)識(shí)。從中可以看出,計(jì)算機(jī)仿真技術(shù)在環(huán)境試驗(yàn)設(shè)備研制中大有作為,將來(lái)伴隨著仿真技術(shù)的發(fā)展,環(huán)境試驗(yàn)設(shè)備研制也會(huì)越來(lái)越快速、準(zhǔn)確。

    [1] 李興瑋, 曹娟. 仿真計(jì)算機(jī)的過(guò)去?現(xiàn)在和未來(lái)[J]. 系統(tǒng)仿真學(xué)報(bào), 2009,21(S2): 106-111.

    [2] 李坦. 位移及加速度控制的液壓伺服振動(dòng)合的動(dòng)態(tài)優(yōu)化設(shè)計(jì)[J]. 試驗(yàn)機(jī)與材料試驗(yàn), 1985(2): 15-22.

    [3] 李長(zhǎng)安, 王長(zhǎng)江, 何長(zhǎng)安. 不對(duì)稱結(jié)構(gòu)的防滑剎車(chē)壓力伺服閥的仿真研究[C]// 97中國(guó)航空學(xué)會(huì)液壓與氣動(dòng)專(zhuān)業(yè)學(xué)術(shù)討論會(huì). 1997.

    [4] 吳志清. 四通閥控非對(duì)稱液壓缸伺服系統(tǒng)的動(dòng)態(tài)特性[J]. 煤礦機(jī)械, 2004(7): 32-34.

    [5] 江桂云, 王勇勤, 嚴(yán)興春. 液壓伺服閥控缸動(dòng)態(tài)特性數(shù)學(xué)建模及仿真分析[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版), 2008(5): 195-198.

    [6] 任杉, 李瑋, 韓青. 基于MATLAB 的閥控缸伺服系統(tǒng)仿真分析[J]. 裝備制造技術(shù), 2009(9): 32-33.

    [7] 徐敏, 薛敏東, 王曦. 噴嘴擋板式力反饋電液伺服閥非線性動(dòng)態(tài)模型[C]// 2010年中國(guó)國(guó)防科技會(huì)議. 西安, 2010.

    [8] 關(guān)景泰. 機(jī)電液控制技術(shù)[M]. 上海: 同濟(jì)大學(xué)出版社, 2003.

    [9] 牛寶良, 王玨. 液壓振動(dòng)臺(tái)時(shí)域建模與仿真[J]. 裝備環(huán)境工程, 2010, 7(6): 247-251.

    [10] 鈕王杰, 張斌超, 劉春艷. 基于MATLAB/Simulink轉(zhuǎn)速閉環(huán)直流調(diào)速控制系統(tǒng)的仿真研究[J]. 山西大同大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 27(5): 33-35.

    [11] 牛寶良. 一種計(jì)算結(jié)構(gòu)時(shí)程響應(yīng)的簡(jiǎn)單直接積分法[C]// 第23屆全國(guó)結(jié)構(gòu)工程學(xué)術(shù)會(huì)議. 蘭州,2014.

    [12] 牛寶良. 某離心機(jī)模型辨識(shí)與仿真[J]. 裝備環(huán)境工程, 2015, 12(5): 61-66.

    [13] 劉增光, 岳大靈, 安林超, 等. 基于MATLAB 的力反饋兩級(jí)電液伺服閥建模與仿真[J]. 液壓與氣動(dòng), 2015(5): 83-85.

    [14] 張宏, 高銘, 郭凱. 伺服閥開(kāi)口度對(duì)電液伺服系統(tǒng)動(dòng)態(tài)特性的影響[J]. 液壓氣動(dòng)與密封, 2015(10): 27-30.

    [15] 張軍, 左哲清, 王書(shū)銘. 數(shù)字伺服閥卡滯問(wèn)題研究[J]. 液壓與氣動(dòng), 2014(6): 105-107.

    [16] 沈偉, 崔霞. 比例伺服閥靜動(dòng)態(tài)特性仿真建模[J]. 計(jì)算機(jī)仿真, 2015, 32(2): 236-240.

    Dynamic Environment Test Equipment System Design Based on Simulating

    NIU Bao-liang

    (Institute of System Engineering, China Academy of Engineering Physics, Mianyang 621900, China)

    ObjectiveTo introduce characteristics of simulation technology and its application in designing dynamic equipment system. Methods Development, characteristics and experience of simulation technology were introduced firstly. Application of Simulink software in dynamic equipment system was introduced with several examples. Results Main characteristic parameters concerned could be obtained quickly with Simulink simulation. It gave dependable data to evaluate or optimize the equipment system. Conclusion The computer simulation technology is quite mature and can be used in the dynamic environment test system design to advance the design level.

    simulation; equipment; dynamic system; algebraic loops; environmental test

    10.7643/ issn.1672-9242.2017.08.001

    TJ01;TB114

    A

    1672-9242(2017)08-0001-10

    2017-04-18;

    2017-05-15

    牛寶良 (1963—),男,陜西人,碩士,研究員,主要從事振動(dòng)、離心試驗(yàn)及相關(guān)設(shè)備研發(fā)。

    猜你喜歡
    振動(dòng)臺(tái)仿真技術(shù)加速度
    “鱉”不住了!從26元/斤飆至38元/斤,2022年甲魚(yú)能否再跑出“加速度”?
    基于振動(dòng)臺(tái)試驗(yàn)的通信機(jī)柜地震易損性分析
    虛擬仿真技術(shù)支持情景構(gòu)建
    天際加速度
    創(chuàng)新,動(dòng)能轉(zhuǎn)換的“加速度”
    金橋(2018年4期)2018-09-26 02:24:46
    死亡加速度
    AMESim仿真技術(shù)在汽車(chē)空調(diào)制冷系統(tǒng)中的應(yīng)用
    大型液壓離心振動(dòng)臺(tái)控制策略的仿真研究
    應(yīng)用三維仿真技術(shù)模擬彈子鎖異常開(kāi)啟
    420 kV避雷器振動(dòng)臺(tái)抗震試驗(yàn)
    女警被强在线播放| 狂野欧美激情性bbbbbb| 午夜福利视频在线观看免费| 国产精品九九99| 老汉色av国产亚洲站长工具| 成年av动漫网址| 国产高清不卡午夜福利| 成人亚洲精品一区在线观看| 国产三级黄色录像| 性高湖久久久久久久久免费观看| 久久中文字幕一级| 黑丝袜美女国产一区| 欧美xxⅹ黑人| 亚洲熟女精品中文字幕| 国产精品一国产av| 亚洲色图 男人天堂 中文字幕| 亚洲av国产av综合av卡| 国产成人精品无人区| 免费观看av网站的网址| 满18在线观看网站| 亚洲av国产av综合av卡| 免费观看a级毛片全部| 女警被强在线播放| 一本—道久久a久久精品蜜桃钙片| 精品熟女少妇八av免费久了| av片东京热男人的天堂| 国产激情久久老熟女| 在线观看国产h片| 中文字幕av电影在线播放| 咕卡用的链子| 日韩 亚洲 欧美在线| 午夜视频精品福利| tube8黄色片| 成人午夜精彩视频在线观看| 久久久久久久久久久久大奶| 黑人猛操日本美女一级片| 老司机在亚洲福利影院| 91麻豆精品激情在线观看国产 | 免费看十八禁软件| 人人澡人人妻人| 又紧又爽又黄一区二区| 国产欧美日韩一区二区三 | 久9热在线精品视频| 亚洲av电影在线进入| 丝袜人妻中文字幕| 两个人看的免费小视频| 免费在线观看影片大全网站 | 午夜久久久在线观看| 人妻人人澡人人爽人人| 国产淫语在线视频| 一区二区三区精品91| 久久久久久久久免费视频了| 亚洲精品国产区一区二| 男女床上黄色一级片免费看| 视频区图区小说| 91精品伊人久久大香线蕉| 亚洲专区中文字幕在线| 你懂的网址亚洲精品在线观看| 操美女的视频在线观看| 国产片特级美女逼逼视频| 日韩精品免费视频一区二区三区| 十八禁人妻一区二区| 国产精品一区二区在线不卡| 日韩中文字幕欧美一区二区 | 亚洲国产最新在线播放| 极品少妇高潮喷水抽搐| av网站免费在线观看视频| 免费久久久久久久精品成人欧美视频| 国产精品国产三级国产专区5o| 欧美老熟妇乱子伦牲交| 18禁黄网站禁片午夜丰满| 少妇猛男粗大的猛烈进出视频| 国产爽快片一区二区三区| 日韩熟女老妇一区二区性免费视频| 国产极品粉嫩免费观看在线| 亚洲av片天天在线观看| 久久精品熟女亚洲av麻豆精品| 国产成人精品久久二区二区91| av网站在线播放免费| 欧美日韩精品网址| 尾随美女入室| 日本vs欧美在线观看视频| 欧美日韩国产mv在线观看视频| 亚洲成人免费av在线播放| 中文字幕人妻丝袜制服| 国产精品三级大全| 悠悠久久av| 女人久久www免费人成看片| 老司机午夜十八禁免费视频| 久久久久网色| 国产成人一区二区在线| 亚洲国产av影院在线观看| 精品国产国语对白av| 丝袜美腿诱惑在线| 欧美中文综合在线视频| 韩国高清视频一区二区三区| 人人妻人人澡人人爽人人夜夜| 日本vs欧美在线观看视频| 夜夜骑夜夜射夜夜干| 巨乳人妻的诱惑在线观看| 大型av网站在线播放| 亚洲欧洲日产国产| 久久久久久免费高清国产稀缺| 日本黄色日本黄色录像| 狂野欧美激情性xxxx| 曰老女人黄片| 捣出白浆h1v1| 中文字幕色久视频| 欧美成人午夜精品| 国产男女超爽视频在线观看| 国产激情久久老熟女| 日韩一区二区三区影片| 看十八女毛片水多多多| 精品欧美一区二区三区在线| 亚洲av综合色区一区| 大香蕉久久网| 桃花免费在线播放| 久热爱精品视频在线9| 亚洲第一青青草原| 性少妇av在线| 9热在线视频观看99| 国产一卡二卡三卡精品| 亚洲人成电影免费在线| 中文乱码字字幕精品一区二区三区| 成人午夜精彩视频在线观看| 国产成人av教育| 国产在视频线精品| 欧美激情 高清一区二区三区| 99精品久久久久人妻精品| 亚洲欧美一区二区三区久久| 国产又爽黄色视频| 成人黄色视频免费在线看| 国产片内射在线| 色视频在线一区二区三区| 亚洲一码二码三码区别大吗| 久久99精品国语久久久| 麻豆国产av国片精品| 成人国产av品久久久| 韩国高清视频一区二区三区| 成年动漫av网址| 午夜福利在线免费观看网站| 国产国语露脸激情在线看| 日日夜夜操网爽| 丰满饥渴人妻一区二区三| 黄色 视频免费看| 国产亚洲av片在线观看秒播厂| 大香蕉久久网| 99热国产这里只有精品6| 日日夜夜操网爽| 婷婷色综合www| 亚洲国产精品国产精品| 国产欧美亚洲国产| 亚洲视频免费观看视频| 女性生殖器流出的白浆| 亚洲黑人精品在线| 精品久久久久久电影网| 欧美精品av麻豆av| 欧美精品啪啪一区二区三区 | 青春草亚洲视频在线观看| 国产成人精品在线电影| 在线观看免费日韩欧美大片| 校园人妻丝袜中文字幕| 久久精品久久久久久久性| 午夜免费观看性视频| 国产精品av久久久久免费| 欧美日韩一级在线毛片| av天堂久久9| 手机成人av网站| 晚上一个人看的免费电影| 国产精品香港三级国产av潘金莲 | 99热网站在线观看| 啦啦啦中文免费视频观看日本| 侵犯人妻中文字幕一二三四区| av在线老鸭窝| 国产精品.久久久| 免费观看a级毛片全部| 亚洲欧洲日产国产| 不卡av一区二区三区| 人人妻,人人澡人人爽秒播 | 午夜av观看不卡| 亚洲精品av麻豆狂野| 国产亚洲精品第一综合不卡| 亚洲激情五月婷婷啪啪| 午夜日韩欧美国产| 欧美日韩亚洲高清精品| 国产高清视频在线播放一区 | 在线观看免费午夜福利视频| 国产男女内射视频| 制服人妻中文乱码| 在现免费观看毛片| 国产野战对白在线观看| 亚洲人成77777在线视频| 别揉我奶头~嗯~啊~动态视频 | 丁香六月欧美| 少妇粗大呻吟视频| 欧美黄色片欧美黄色片| 日韩大码丰满熟妇| √禁漫天堂资源中文www| 老司机影院成人| 国产亚洲精品第一综合不卡| 成人亚洲欧美一区二区av| 极品人妻少妇av视频| 亚洲精品久久午夜乱码| 日韩欧美一区视频在线观看| 两人在一起打扑克的视频| 亚洲五月色婷婷综合| 国产免费现黄频在线看| 黄片播放在线免费| 精品一品国产午夜福利视频| 亚洲av欧美aⅴ国产| 又粗又硬又长又爽又黄的视频| 妹子高潮喷水视频| 亚洲国产精品一区二区三区在线| 精品一区在线观看国产| 18禁黄网站禁片午夜丰满| 啦啦啦 在线观看视频| 99精国产麻豆久久婷婷| 国产精品国产三级专区第一集| 黄色片一级片一级黄色片| 在线观看国产h片| av在线播放精品| 汤姆久久久久久久影院中文字幕| 久久久国产一区二区| 精品卡一卡二卡四卡免费| 亚洲欧美日韩另类电影网站| 国语对白做爰xxxⅹ性视频网站| 亚洲国产毛片av蜜桃av| 日韩av在线免费看完整版不卡| 两个人看的免费小视频| 婷婷色av中文字幕| 性高湖久久久久久久久免费观看| 国产极品粉嫩免费观看在线| netflix在线观看网站| 亚洲精品中文字幕在线视频| 男女午夜视频在线观看| 欧美少妇被猛烈插入视频| 无限看片的www在线观看| avwww免费| 亚洲一码二码三码区别大吗| 精品一区二区三区四区五区乱码 | 欧美日韩视频高清一区二区三区二| 另类精品久久| 精品人妻一区二区三区麻豆| 国产精品.久久久| 两性夫妻黄色片| 亚洲国产欧美在线一区| 亚洲av男天堂| 欧美黄色淫秽网站| 欧美精品av麻豆av| 手机成人av网站| 两个人免费观看高清视频| 欧美激情高清一区二区三区| 欧美黄色淫秽网站| 免费在线观看日本一区| 免费久久久久久久精品成人欧美视频| 久热爱精品视频在线9| 熟女av电影| 视频在线观看一区二区三区| 久久久国产一区二区| 狂野欧美激情性bbbbbb| 高清视频免费观看一区二区| 亚洲欧美一区二区三区黑人| 人人妻人人澡人人爽人人夜夜| 久久精品久久久久久噜噜老黄| 国产一区二区 视频在线| 国产成人欧美| 建设人人有责人人尽责人人享有的| 成年人免费黄色播放视频| 亚洲免费av在线视频| 亚洲精品第二区| 成人亚洲精品一区在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲中文字幕日韩| 久久久国产精品麻豆| 777米奇影视久久| 久久人人爽人人片av| 性高湖久久久久久久久免费观看| 国产精品亚洲av一区麻豆| 亚洲熟女毛片儿| 欧美日韩综合久久久久久| 国产成人啪精品午夜网站| 久久精品久久久久久久性| 国产伦人伦偷精品视频| 久久ye,这里只有精品| 午夜两性在线视频| 午夜福利,免费看| 色精品久久人妻99蜜桃| av天堂久久9| 黄网站色视频无遮挡免费观看| 深夜精品福利| 在线观看免费日韩欧美大片| 激情五月婷婷亚洲| 女性被躁到高潮视频| 80岁老熟妇乱子伦牲交| 18在线观看网站| 男女高潮啪啪啪动态图| 婷婷色综合www| 日韩制服骚丝袜av| 校园人妻丝袜中文字幕| 欧美日韩亚洲高清精品| 日本a在线网址| 91老司机精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人影院久久av| 大片免费播放器 马上看| 1024香蕉在线观看| 一级毛片 在线播放| 亚洲五月色婷婷综合| 久久中文字幕一级| 多毛熟女@视频| 亚洲精品日本国产第一区| 国产精品久久久av美女十八| 久久久久久久久久久久大奶| 日本猛色少妇xxxxx猛交久久| 国语对白做爰xxxⅹ性视频网站| 嫩草影视91久久| 人妻一区二区av| 免费看十八禁软件| 青春草亚洲视频在线观看| 久久久久久久大尺度免费视频| 亚洲成色77777| 可以免费在线观看a视频的电影网站| 下体分泌物呈黄色| 久久精品国产亚洲av涩爱| 国产成人91sexporn| 欧美性长视频在线观看| 久久av网站| 成人国产av品久久久| 国产精品久久久久久精品电影小说| 婷婷色av中文字幕| 19禁男女啪啪无遮挡网站| 欧美黄色淫秽网站| 国产精品国产三级国产专区5o| 五月开心婷婷网| svipshipincom国产片| 国产亚洲av片在线观看秒播厂| 两性夫妻黄色片| 2021少妇久久久久久久久久久| 中文字幕人妻丝袜一区二区| 日韩一卡2卡3卡4卡2021年| 啦啦啦在线观看免费高清www| 亚洲精品久久午夜乱码| 国产男人的电影天堂91| 黄网站色视频无遮挡免费观看| 免费女性裸体啪啪无遮挡网站| 亚洲欧美日韩高清在线视频 | 国产精品免费大片| 男女午夜视频在线观看| 免费高清在线观看视频在线观看| 夜夜骑夜夜射夜夜干| 十八禁网站网址无遮挡| www.精华液| 久久久久久人人人人人| 成年美女黄网站色视频大全免费| 日本猛色少妇xxxxx猛交久久| 欧美日韩黄片免| 成人免费观看视频高清| 亚洲成国产人片在线观看| 国产精品国产三级国产专区5o| 久久鲁丝午夜福利片| 亚洲成国产人片在线观看| 国产亚洲一区二区精品| 国产精品.久久久| 又粗又硬又长又爽又黄的视频| 久久久久久人人人人人| 电影成人av| 色视频在线一区二区三区| 男女午夜视频在线观看| 日韩av在线免费看完整版不卡| 交换朋友夫妻互换小说| 9色porny在线观看| 精品亚洲成国产av| 十八禁人妻一区二区| 无限看片的www在线观看| 久久久久精品人妻al黑| 国产亚洲欧美精品永久| 女人被躁到高潮嗷嗷叫费观| 国产成人精品久久二区二区91| 精品第一国产精品| 国产精品 欧美亚洲| 中文字幕最新亚洲高清| 久久久久久久大尺度免费视频| 香蕉丝袜av| 精品一区在线观看国产| 好男人电影高清在线观看| 成人亚洲欧美一区二区av| 大型av网站在线播放| 99久久综合免费| 热99久久久久精品小说推荐| 七月丁香在线播放| 一级黄片播放器| 夫妻性生交免费视频一级片| www.999成人在线观看| 成人国产av品久久久| 国产亚洲精品第一综合不卡| 久久精品成人免费网站| 亚洲精品日本国产第一区| 美女脱内裤让男人舔精品视频| 日日摸夜夜添夜夜爱| 欧美激情 高清一区二区三区| 国产激情久久老熟女| 国产免费又黄又爽又色| 各种免费的搞黄视频| 老熟女久久久| 精品国产乱码久久久久久小说| 欧美亚洲日本最大视频资源| 一区二区三区精品91| 秋霞在线观看毛片| 欧美精品人与动牲交sv欧美| 最近手机中文字幕大全| 视频区欧美日本亚洲| 校园人妻丝袜中文字幕| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 国产成人系列免费观看| 国产97色在线日韩免费| 午夜精品国产一区二区电影| 18在线观看网站| 一级黄色大片毛片| 香蕉丝袜av| 多毛熟女@视频| 黑人巨大精品欧美一区二区蜜桃| 国产黄色免费在线视频| 人妻人人澡人人爽人人| 一二三四社区在线视频社区8| 又粗又硬又长又爽又黄的视频| 乱人伦中国视频| 亚洲人成电影观看| 老司机亚洲免费影院| 国产欧美亚洲国产| 丰满迷人的少妇在线观看| 国产精品久久久人人做人人爽| 美女扒开内裤让男人捅视频| 99久久人妻综合| 久久99精品国语久久久| 久久人人爽人人片av| 1024视频免费在线观看| 国产精品一区二区精品视频观看| 国产成人精品久久久久久| 中文字幕高清在线视频| 1024香蕉在线观看| 欧美黄色淫秽网站| 国产亚洲精品第一综合不卡| 99精品久久久久人妻精品| 久久久精品区二区三区| 亚洲av欧美aⅴ国产| 又黄又粗又硬又大视频| 精品久久久精品久久久| 免费久久久久久久精品成人欧美视频| 欧美成人午夜精品| 一区二区三区乱码不卡18| 人妻一区二区av| 欧美精品av麻豆av| 好男人电影高清在线观看| bbb黄色大片| 中文字幕制服av| 成人黄色视频免费在线看| 欧美成人午夜精品| 国产精品久久久久久人妻精品电影 | 老司机亚洲免费影院| 香蕉丝袜av| 国精品久久久久久国模美| 久久久久精品人妻al黑| 亚洲欧美日韩高清在线视频 | 老司机亚洲免费影院| 日韩,欧美,国产一区二区三区| 久久人人爽av亚洲精品天堂| av视频免费观看在线观看| 色综合欧美亚洲国产小说| 国产日韩欧美亚洲二区| 婷婷色综合大香蕉| 免费高清在线观看视频在线观看| 亚洲精品一二三| 伊人亚洲综合成人网| 两个人看的免费小视频| 国产欧美日韩一区二区三区在线| 国产一级毛片在线| av电影中文网址| videos熟女内射| 亚洲一区二区三区欧美精品| 天堂8中文在线网| 国产无遮挡羞羞视频在线观看| 亚洲av日韩精品久久久久久密 | 欧美精品一区二区免费开放| 伦理电影免费视频| 丝瓜视频免费看黄片| 亚洲国产欧美一区二区综合| 久久青草综合色| 国产一区有黄有色的免费视频| 老鸭窝网址在线观看| 脱女人内裤的视频| 一个人免费看片子| 1024视频免费在线观看| 性色av乱码一区二区三区2| 国产成人影院久久av| 啦啦啦在线观看免费高清www| 精品少妇黑人巨大在线播放| 久久av网站| 日本猛色少妇xxxxx猛交久久| 悠悠久久av| 又大又爽又粗| 欧美日韩精品网址| 国产成人精品久久久久久| av国产精品久久久久影院| a级毛片黄视频| 久久99一区二区三区| 老司机亚洲免费影院| 永久免费av网站大全| 久久久久网色| 51午夜福利影视在线观看| 夜夜骑夜夜射夜夜干| 欧美日韩亚洲综合一区二区三区_| 久久国产亚洲av麻豆专区| 高清视频免费观看一区二区| 成人影院久久| 亚洲国产毛片av蜜桃av| 女性生殖器流出的白浆| 国产成人免费无遮挡视频| 成年美女黄网站色视频大全免费| 狂野欧美激情性xxxx| 97在线人人人人妻| 国产一区亚洲一区在线观看| 亚洲欧美色中文字幕在线| 午夜激情久久久久久久| 在线观看免费视频网站a站| 欧美国产精品一级二级三级| 赤兔流量卡办理| 中国国产av一级| 蜜桃在线观看..| 亚洲精品国产色婷婷电影| 国产主播在线观看一区二区 | 最近手机中文字幕大全| 国产一区二区三区综合在线观看| 激情五月婷婷亚洲| 人妻 亚洲 视频| 多毛熟女@视频| 久久国产精品男人的天堂亚洲| 一级毛片 在线播放| 国产精品麻豆人妻色哟哟久久| 91麻豆av在线| av欧美777| 欧美久久黑人一区二区| 女人久久www免费人成看片| 国精品久久久久久国模美| 日本黄色日本黄色录像| 久久 成人 亚洲| 欧美日本中文国产一区发布| 国产成人精品久久二区二区免费| 性高湖久久久久久久久免费观看| 国精品久久久久久国模美| 欧美日韩综合久久久久久| 9热在线视频观看99| 国产色视频综合| 另类精品久久| 少妇人妻久久综合中文| 国产又爽黄色视频| 搡老乐熟女国产| 欧美黑人欧美精品刺激| 一级毛片女人18水好多 | 91老司机精品| 捣出白浆h1v1| 日日摸夜夜添夜夜爱| videosex国产| 精品久久久精品久久久| 国产免费福利视频在线观看| 在现免费观看毛片| 久久免费观看电影| 久久精品久久久久久久性| 国产熟女午夜一区二区三区| 国产一级毛片在线| 国产精品秋霞免费鲁丝片| 无遮挡黄片免费观看| 欧美中文综合在线视频| 久久亚洲国产成人精品v| 天堂8中文在线网| 大香蕉久久网| 国产成人91sexporn| 国产精品一区二区免费欧美 | 欧美精品高潮呻吟av久久| 国产男人的电影天堂91| 一区福利在线观看| 黑人猛操日本美女一级片| 亚洲,欧美精品.| 国产av国产精品国产| 精品国产一区二区三区久久久樱花| 无限看片的www在线观看| 纵有疾风起免费观看全集完整版| 欧美日韩精品网址| 国产精品国产三级专区第一集| 久久天躁狠狠躁夜夜2o2o | 国产野战对白在线观看| 啦啦啦在线观看免费高清www| 一区二区日韩欧美中文字幕| 欧美大码av| 日本黄色日本黄色录像| 免费一级毛片在线播放高清视频 | 久久久国产一区二区| 80岁老熟妇乱子伦牲交| 母亲3免费完整高清在线观看| 91麻豆av在线| 日韩欧美一区视频在线观看| av国产精品久久久久影院| 黄频高清免费视频| 日韩一区二区三区影片| 久久久精品国产亚洲av高清涩受| 国产xxxxx性猛交| 亚洲欧美日韩高清在线视频 | 丝袜美足系列| 久久久久久久久久久久大奶| 女人被躁到高潮嗷嗷叫费观| 精品欧美一区二区三区在线| 亚洲,欧美精品.| 国产免费一区二区三区四区乱码| 欧美精品人与动牲交sv欧美| 国产一级毛片在线|