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

    基于改進(jìn)Bingham模型的磁流變阻尼器力學(xué)建模及試驗(yàn)研究

    2021-05-12 04:34:28祝世興楊麗昆祝恒佳
    關(guān)鍵詞:剪切應(yīng)力建模振動(dòng)

    祝世興,楊麗昆,魏 戩,祝恒佳

    (中國民航大學(xué)航空工程學(xué)院,天津 300300)

    磁流變阻尼器(magnetorheological damper,MRD)是一種阻尼可調(diào)的半主動(dòng)智能控制器件,其傳動(dòng)介質(zhì)磁流變液(magnetorheological fluid,MRF)是一種受磁場(chǎng)影響的特殊材料[1]。MRF的黏度可隨外加磁場(chǎng)強(qiáng)度不同而迅速發(fā)生變化,進(jìn)而改變輸出阻尼力的大小。因此MRD的輸出阻尼可以在外加磁場(chǎng)(一般通過電流大小來控制)作用下實(shí)現(xiàn)動(dòng)態(tài)調(diào)節(jié),并以其結(jié)構(gòu)簡單、可調(diào)范圍寬、響應(yīng)快、可靠性和適應(yīng)性強(qiáng)等優(yōu)點(diǎn)被廣泛應(yīng)用于建筑[2-4]、汽車[5-8]、航空[9-10]、航海[11-12]等領(lǐng)域。

    一般阻尼器多應(yīng)用于結(jié)構(gòu)的減震、隔振及減擺等維持系統(tǒng)的穩(wěn)定性方面。在削弱振動(dòng)過程中往往由于非線性因素[13](本文MRD的非線性因素主要是MRF)的存在,使系統(tǒng)的運(yùn)動(dòng)具有非線性特性,非線性特性對(duì)多體系統(tǒng)動(dòng)力學(xué)建模[14-15]的影響是不可避免的。黎崛珉等[16]研究了剛度非線性與阻尼非線性對(duì)隔振系統(tǒng)的影響。姚文莉等[17]將非理想系統(tǒng)的動(dòng)力學(xué)建模問題描述為目標(biāo)函數(shù)為高斯拘束函數(shù)、優(yōu)化變量為質(zhì)點(diǎn)加速度的約束最優(yōu)化問題,并驗(yàn)證其方法的有效性。Zhu[18,19]等人研究了帶有阻尼器的空氣彈簧非線性動(dòng)力學(xué)模型,用來描述車輛減震系統(tǒng)以提高乘坐舒適性。

    為了準(zhǔn)確模擬MRD復(fù)雜的非線性阻尼力特性,實(shí)現(xiàn)對(duì)MRD精準(zhǔn)控制,研究人員基于多體動(dòng)力學(xué)[20]及MRF的材料特性對(duì)MRD的力學(xué)建模方法進(jìn)行了深入探究[21-22]。Shou[23]、李趙春[24]等人研究了MRD在沖擊載荷作用下動(dòng)力學(xué)響應(yīng)特性模型。Ma[25]、Gurubasavaraju[26]、Yang[27]等研究了MRF剪切應(yīng)力模型,對(duì)MRF電導(dǎo)率特性進(jìn)行了 實(shí) 驗(yàn) 研 究。Zolfagharian[28]、Peng[29]、陳 世嵬[30]等人對(duì)MRD進(jìn)行參數(shù)化建模并對(duì)模型的參數(shù)辨識(shí)方法進(jìn)行了研究。龔微[31]等提出了一種可以求解阻尼器逆模型的3段線性變阻尼恢復(fù)力模型,適用于MRD控制器的設(shè)計(jì)。梅真等[32]對(duì)MRD進(jìn)行動(dòng)力性能測(cè)試,建立了參數(shù)化模型和非參數(shù)化模型,通過對(duì)比2種建模方式的結(jié)果,發(fā)現(xiàn)非參數(shù)化模型的精度高于參數(shù)化模型。

    以上研究結(jié)果表明,由于MRF的流變機(jī)理復(fù)雜,目前尚無統(tǒng)一、精確的力學(xué)模型描述,使得參數(shù)化模型的建模方式復(fù)雜化。非參數(shù)化模型是基于實(shí)驗(yàn)測(cè)試結(jié)果建立的,利用數(shù)值計(jì)算方法對(duì)模型參數(shù)進(jìn)行辨識(shí),得到較為精確化的模型,可簡化建模方式和提高模型精度??紫闁|等[33]提出了Bingham-多項(xiàng)式力學(xué)模型,并與試驗(yàn)測(cè)試結(jié)果進(jìn)行對(duì)比,解決了Bingham模型精度低、多項(xiàng)式模型高階不穩(wěn)定的問題。祝世興等[34]提出用傅立葉級(jí)數(shù)來表達(dá)MRD的非線性關(guān)系,改進(jìn)了線性模型精度差的問題。劉永強(qiáng)等[35]提出一種對(duì)MRD模型進(jìn)行參數(shù)辨識(shí)的新方法,分析模型參數(shù)對(duì)MRD輸出阻尼的影響規(guī)律,獲取調(diào)控阻尼力的最佳參數(shù),并擬合參數(shù)值和電流之間的函數(shù)關(guān)系。胡海剛等[36]提出一種將遺傳算法和最小二乘法相結(jié)合的方法,對(duì)MRD的力學(xué)模型進(jìn)行參數(shù)辨識(shí),并用Simulink工具箱進(jìn)行仿真驗(yàn)證,所得結(jié)果精度較高,可靠性滿足要求。周勇等[37]提出一種改進(jìn)型RBF神經(jīng)網(wǎng)絡(luò)模型,基于實(shí)驗(yàn)實(shí)測(cè)結(jié)果,利用線性插值法建立模型算法,并與改進(jìn)前的模型進(jìn)行對(duì)比,發(fā)現(xiàn)改進(jìn)后模型精度更高,可以更加可靠地描述MRD的動(dòng)力學(xué)特性。

    現(xiàn)有的非參數(shù)化模型參數(shù)較多,建模方式較為復(fù)雜,分段模型存在分段方式模糊、拐點(diǎn)不明確等問題。為了精確描述磁流變阻尼器在不同電流下力與速度之間的非線性關(guān)系、減少力學(xué)模型參數(shù)、簡化分段模型的分段方式,本文作者通過對(duì)臺(tái)架試驗(yàn)的實(shí)驗(yàn)結(jié)果進(jìn)行分析,建立Rational-Bingham模型,采用最小二乘法對(duì)模型參數(shù)進(jìn)行辨識(shí),基于臺(tái)架試驗(yàn)結(jié)果對(duì)Rational-Bingham模型和多項(xiàng)式-Bingham模型的精度進(jìn)行對(duì)比研究。

    1 實(shí)驗(yàn)

    實(shí)驗(yàn)所采用的磁流變阻尼器結(jié)構(gòu)原理圖如圖1,其中缸筒內(nèi)填充MRF。

    圖1 磁流變阻尼器結(jié)構(gòu)原理圖

    MRD實(shí)物圖如圖2,通過活塞桿引出的導(dǎo)線可以為MRD施加不同大小的電流。阻尼特性測(cè)試平臺(tái)如圖3所示,將MRD的上端固定,下端活塞端連接激振臺(tái)中心,激振臺(tái)通過液壓作動(dòng)為阻尼器提供正弦位移激勵(lì)。

    實(shí)驗(yàn)采用正弦位移激勵(lì),以振幅為2 mm、頻率為3 Hz為例,通過圖2引出的導(dǎo)線分別對(duì)MRD施加0、0.5、1.0、1.5、2.0、2.5 A的電流,研究在不同磁場(chǎng)強(qiáng)度下MRD的阻尼特性。將不同電流下力-位移關(guān)系轉(zhuǎn)換為力-速度關(guān)系,結(jié)果如圖4所示,其中橫坐標(biāo)表示MRD的輸入速度v,縱坐標(biāo)表示MRD的輸出阻尼力F。

    圖2 實(shí)物圖

    圖3 阻尼特性測(cè)試平臺(tái)

    圖4 不同電流下MRD力-速度關(guān)系

    從不同電流下MRD的動(dòng)態(tài)響應(yīng)曲線可以看出,將速度換向處(A處和C處)和加速度換向處(B處和D處)作為分段點(diǎn)進(jìn)行分段建??蓪⒁粋€(gè)運(yùn)動(dòng)周期分為2段線性曲線和2段為非線性曲線。

    2 Bingham模型中MRF力學(xué)特性

    MRF在無外加磁場(chǎng)的影響下可近似為牛頓流體,此時(shí)MRF的剪切應(yīng)力與剪切應(yīng)變率之間的關(guān)系為:

    式中:τ為剪切應(yīng)力;η為黏度系數(shù);γ為剪切應(yīng)變率。

    MRF在外加磁場(chǎng)影響下具有黏塑性流體特性,如圖5所示,在不同的磁場(chǎng)強(qiáng)度下,流動(dòng)模式可分為屈服前區(qū)、屈服區(qū)、屈服后區(qū)等3種變化過程。其中MRF在屈服前區(qū)剪切應(yīng)力隨剪切應(yīng)變率的變化呈正相關(guān);在屈服區(qū)剪切應(yīng)力隨剪切應(yīng)變率的變化呈不穩(wěn)定狀態(tài);在屈服后區(qū)呈穩(wěn)態(tài),剪切應(yīng)力不隨剪切應(yīng)變率變化。

    圖5 MRF的剪切應(yīng)力與剪切應(yīng)變率關(guān)系

    出于本文的實(shí)際工況及MRD的工作狀態(tài)考慮,選取穩(wěn)態(tài)工作狀態(tài)下MRD中MRF的所處狀態(tài)進(jìn)行研究,此狀態(tài)下MRF穩(wěn)定流動(dòng)處于屈服后區(qū),剪切應(yīng)力為定值,MRF黏度為定值。

    Bingham模型將MRF的剪切應(yīng)力描述為:

    式中:τy為MRF的剪切應(yīng)力;τ0為MRF的臨界剪切應(yīng)力;H為MRF所處磁場(chǎng)強(qiáng)度;γ為MRF的剪切應(yīng)變率;η為MRF達(dá)到屈服后區(qū)的黏度系數(shù)。

    Bingham模型中MRF在不同磁場(chǎng)強(qiáng)度影響下剪切應(yīng)力與剪切應(yīng)變率的量綱為一化關(guān)系如圖6所示。

    圖6 Bingham模型中MRF本構(gòu)關(guān)系

    由圖6可知,磁場(chǎng)強(qiáng)度的變化會(huì)改變Bingham模型中的臨界剪切應(yīng)力,而MRF在剪切屈服前后的黏度η在Bingham模型中被視為定值。

    描述MRD輸出力與激勵(lì)速度之間線性關(guān)系的常用 Bingham力學(xué)模型,是由 Stanway等[38]提出的由庫倫摩擦力和黏滯阻尼器構(gòu)成的理想化模型,即:

    式中:fk為庫倫阻尼力,與MRF的屈服應(yīng)力相關(guān);x為MRD的輸出位移;˙x為MRD的輸出速度;c0為黏滯阻尼系數(shù);f0為補(bǔ)償元件存在而產(chǎn)生的力,結(jié)合本文MRD實(shí)際結(jié)構(gòu)及實(shí)際工況此項(xiàng)可忽略。

    將圖4的實(shí)驗(yàn)結(jié)果與Bingham力學(xué)模型進(jìn)行對(duì)比分析,可以發(fā)現(xiàn):圖4的線性部分在理論上近似表示MRF達(dá)到屈服后期且黏度不變時(shí)MRD的運(yùn)動(dòng)狀態(tài),可通過Bingham模型模擬,但其非線性區(qū)域的滯回特性無法用Bingham模型進(jìn)行精確描述。本文針對(duì)Bingham模型對(duì)描述MRD非線性區(qū)域力學(xué)特性存在的局限性問題,對(duì)傳統(tǒng)的Bingham模型進(jìn)行改進(jìn),從而獲得更準(zhǔn)確的磁流變阻尼器非線性力學(xué)模型。

    3 力學(xué)建模與參數(shù)辨識(shí)

    本文在Bingham模型基礎(chǔ)上進(jìn)行分段建模,對(duì)圖4的線性部分采用Bingham模型,非線性部分分別采用多項(xiàng)式函數(shù)和Rational函數(shù)進(jìn)行描述。

    3.1 多項(xiàng)式-Bingham模型建立與參數(shù)辨識(shí)

    對(duì)圖4進(jìn)行分析可知,B-C段和D-A段所表示的MRD輸出力與速度之間的線性關(guān)系可通過Bingham模型進(jìn)行描述,為了在形式上進(jìn)一步簡化分段模型,本文將式(3)中fksgn(˙x)視為一個(gè)整體 f,c0視為 c,則參數(shù) f和 c為待識(shí)別參數(shù)。A-B段和C-D段為非線性部分,首先采用多項(xiàng)式函數(shù)進(jìn)行描述。為了防止高階振蕩失真和更好地描述滯回特性,擬采用3次多項(xiàng)式建模,可得多項(xiàng)式-Bingham模型(以下簡稱為多項(xiàng)式函數(shù)模型),其表達(dá)式為:

    式中:v為MRD的速度;˙v為MRD的加速度;c1、f1、p11、p12…為待識(shí)別參數(shù);式(4)中函數(shù)分別依次對(duì)應(yīng)表示圖4中D-A段、A-B段、B-C段、CD段。

    本文采用最小二程曲線擬合法對(duì)模型參數(shù)進(jìn)行辨識(shí),最小二成曲線擬合法的基本原理是對(duì)給定數(shù)據(jù)(xi,yi)(i=0,1,…,n),在函數(shù)類 φ中,求得g(x)∈φ,使誤差的平方和R最小,可表示為:

    從幾何意義上來說,就是找到與給定點(diǎn)(xi,yi)(i=0,1,…,n)的距離平方和為最小的函數(shù)曲線 y=g(x),函數(shù) g(x)稱為最小二乘解或擬合函數(shù),求擬合函數(shù)g(x)的方法稱為最小二乘曲線擬合法。

    在振幅為2 mm、頻率為3 Hz的振動(dòng)工況下,將加載電流分別為 0、0.5、1.0、1.5、2.0、2.5 A的實(shí)驗(yàn)結(jié)果分別對(duì)多項(xiàng)式函數(shù)模型進(jìn)行參數(shù)辨識(shí),得到不同電流下對(duì)應(yīng)的參數(shù)值,并建立不同參數(shù)與電流之間的關(guān)系。具體過程以加載電流為1.5 A的參數(shù)辨識(shí)結(jié)果為例,利用Matlab擬合工具箱對(duì) CD段系數(shù)為 p21、p22、p23、p24的多項(xiàng)式函數(shù)進(jìn)行曲線擬合,擬合結(jié)果如圖7所示,可得CD段的各參數(shù) p21、p22、p23、p24分別為:0.000 016、0.001 25、0.035 35、0.063 13。

    根據(jù)擬合結(jié)果可以看出,多項(xiàng)式函數(shù)與實(shí)際散點(diǎn)之間的擬合偏差較小,因此三階多項(xiàng)式具有一定的有效性。通過對(duì)多項(xiàng)式函數(shù)模型進(jìn)行參數(shù)辨識(shí),可得到不同電流下多項(xiàng)式函數(shù)模型的各參數(shù)值,部分結(jié)果如表1所示。

    圖7 1.5 A電流下C-D段多項(xiàng)式函數(shù)的擬合結(jié)果

    表1 多項(xiàng)式-Bingham模型參數(shù)辨識(shí)結(jié)果

    通過對(duì)參數(shù)辨識(shí)結(jié)果進(jìn)行分析可知,多項(xiàng)式函數(shù)模型各參數(shù)與電流之間為非線性關(guān)系,作者擬采用三階多項(xiàng)式來建立各個(gè)參數(shù)與電流之間的關(guān)系,可得到各參數(shù)與電流的函數(shù)關(guān)系,部分結(jié)果如下:

    以上為多項(xiàng)式函數(shù)模型中各參數(shù)與電流之間建立的映射關(guān)系,結(jié)合式(4),可以得到完整的多項(xiàng)式函數(shù)模型,其中模型的輸入分別為電流和速度,輸出為對(duì)應(yīng)輸入電流和速度下的輸出阻尼力。

    3.2 Rational-Bingham模型建立與參數(shù)辨識(shí)

    為了更好地描述同一振動(dòng)工況不同電流下MRD的動(dòng)態(tài)響應(yīng),本文建立了有理函數(shù)模型與多項(xiàng)式函數(shù)模型表進(jìn)行對(duì)比研究。本質(zhì)上有理函數(shù)(rational函數(shù))是由2個(gè)多項(xiàng)式的商表示的函數(shù),表示如下:

    式中:P(x)為分子多項(xiàng)式函數(shù);Q(x)為分母多項(xiàng)式函數(shù)。

    當(dāng)多項(xiàng)式擬合效果不夠好時(shí),Rational函數(shù)擬合在精度上有一定的優(yōu)越性[39]。為了減少模型參數(shù),將Rational函數(shù)的分子分母擬采用一次線性函數(shù)進(jìn)行建模,可得到Rational-Bingham模型(以下簡稱為Rational函數(shù)模型)為:

    式中:v為MRD的速度;˙v為MRD的加速度;a1、b1、d1、a2、b2、d2為待識(shí)別參數(shù);式中函數(shù)分別對(duì)應(yīng)表示圖4中 D-A段、A-B段、B-C段、CD段。

    同樣在振幅為2 mm、頻率為3 Hz的振動(dòng)工況下,將加載電流分別為 0、0.5、1.0、1.5、2.0、2.5 A實(shí)驗(yàn)結(jié)果分別對(duì)Rational函數(shù)模型進(jìn)行參數(shù)辨識(shí),得到不同電流下對(duì)應(yīng)的參數(shù)值,并建立不同參數(shù)與電流之間的關(guān)系。具體過程仍以加載電流為1.5 A的參數(shù)辨識(shí)結(jié)果為例,利用Matlab中的擬合工具箱對(duì)CD段系數(shù)為a2、b2、d2的Rational函數(shù)進(jìn)行曲線擬合,擬合結(jié)果如圖8所示,可得CD段的各參數(shù) a2、b2、d2分別為:-0.591 7、0.677 1、-11.78。

    圖8 1.5 A電流下C-D段Rational函數(shù)的擬合結(jié)果

    根據(jù)擬合結(jié)果可以看出,與多項(xiàng)式函數(shù)模型相比,Rational函數(shù)模型與實(shí)際散點(diǎn)的擬合效果好一些,通過對(duì)Rational函數(shù)模型進(jìn)行參數(shù)辨識(shí),得到不同電流下Rational函數(shù)模型中各參數(shù)值,部分結(jié)果如表2所示。

    表2 Rational-Bingham模型參數(shù)辨識(shí)結(jié)果

    同樣,對(duì)Rational函數(shù)模型參數(shù)辨識(shí)結(jié)果進(jìn)行分析可知,各參數(shù)與電流為非線性關(guān)系,選擇采用三階多項(xiàng)式來建立各參數(shù)與電流之間的關(guān)系,得到各參數(shù)與電流之間的函數(shù)關(guān)系,部分結(jié)果如下:

    以上為Rational函數(shù)模型中各參數(shù)與電流之間建立的映射關(guān)系,結(jié)合式(8)可得到完整的Rational函數(shù)模型,與多項(xiàng)式函數(shù)模型相似,Rational函數(shù)模型的輸入分別為電流和速度,輸出為對(duì)應(yīng)輸入電流和速度下的輸出阻尼力。

    4 模型仿真分析

    本文在Matlab/Simulink環(huán)境中對(duì)模型在不同工況下力與速度之間的非線性關(guān)系進(jìn)行仿真分析。如圖9所示,模型輸入分別為可以自定義的輸入電流及正弦位移激勵(lì),模型輸出為MRD輸出阻尼力與速度之間的關(guān)系。其中SubSystem模塊中分別為多項(xiàng)式函數(shù)模型和Rational函數(shù)模型,是仿真系統(tǒng)的子系統(tǒng)。系統(tǒng)可分別輸出在不同工況下多項(xiàng)式函數(shù)模型和Rational函數(shù)模型的F-v仿真結(jié)果。

    圖9 模型的Simulink仿真框圖

    4.1 辨識(shí)工況(電流不同)

    在振幅A為2 mm、頻率f為3 Hz的振動(dòng)工況下,將不同電流分別輸入多項(xiàng)式函數(shù)模型仿真系統(tǒng)和Rational函數(shù)模型仿真系統(tǒng),可分別得到輸出F-v曲線并與實(shí)驗(yàn)結(jié)果生成對(duì)比,受篇幅限制,本節(jié)以施加電流大小分別為0、0.5、1.0、1.5 A為例,所得結(jié)果如圖10所示。

    圖10 辨識(shí)工況下兩模型仿真結(jié)果與實(shí)驗(yàn)結(jié)果

    通過分析圖10中模型仿真結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比,可以得出以下結(jié)果:

    1)從整體上看,在同一振動(dòng)工況不同電流下2種模型的擬合效果都比較好,都可以不同程度較為準(zhǔn)確地描述經(jīng)典Bingham力學(xué)模型無法模擬的非線性滯回區(qū)域。

    2)從局部上看,多項(xiàng)式函數(shù)模型在低于0.5 A時(shí)高速段存在局部不穩(wěn)定性、在速度換向處存在擬合效果不佳的缺陷,而Rational函數(shù)模型可以修正多項(xiàng)式函數(shù)這一缺陷,在穩(wěn)定性上更加符合實(shí)際。

    4.2 驗(yàn)證工況(激勵(lì)頻率、振幅不同)

    為了分別檢驗(yàn)2種模型在同一電流不同振動(dòng)工況下的適用性,將圖9中輸入電流設(shè)為I=1.0 A固定不變,首先改變輸入正弦位移激勵(lì)信號(hào)的頻率,以振幅A為2 mm、頻率f為5 Hz為例,所得多項(xiàng)式函數(shù)模型和Rational函數(shù)模型的仿真結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比如圖11所示。

    圖11 I=1.0 A,A=2 mm,f=5 Hz工況下,兩模型仿真結(jié)果與實(shí)驗(yàn)結(jié)果

    電流I=1.0 A固定不變,改變輸入位移激勵(lì)信號(hào)的振幅,以振幅A為4 mm、頻率f為3 Hz的振動(dòng)工況為例,所得多項(xiàng)式函數(shù)模型和Rational函數(shù)模型的仿真結(jié)果與實(shí)驗(yàn)結(jié)果如圖12所示。

    通過分析改變正弦激勵(lì)輸入信號(hào)的頻率或振幅后,2種模型的仿真結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比(見圖11、圖12),可以得出以下結(jié)果:

    1)多項(xiàng)式函數(shù)模型存在的缺陷更加顯著,改變頻率或振幅后,多項(xiàng)式函數(shù)模型在一個(gè)行程周期的高速段均存在不收斂的現(xiàn)象。線性區(qū)的仿真結(jié)果趨近于實(shí)驗(yàn)結(jié)果,非線性低速區(qū)的仿真結(jié)果存在一定程度的誤差。

    圖12 I=1.0 A,A=4 mm,f=3 Hz工況下,兩模型仿真結(jié)果與實(shí)驗(yàn)結(jié)果

    2)為解決多項(xiàng)式函數(shù)模型在改變振動(dòng)工況后高速段存在不收斂問題,可將多項(xiàng)式函數(shù)與線性函數(shù)相交之后的高速段取線性函數(shù),低速非線性區(qū)用多項(xiàng)式函數(shù)來表達(dá),這將改變分段模型的分段方式,因此多項(xiàng)式函數(shù)模型在本文的分段方式下對(duì)同一電流不同振動(dòng)工況的適用性較弱。

    3)在改變頻率或振幅后,Rational函數(shù)模型在非線性區(qū)的表達(dá)存在不同程度的誤差,但總體仿真結(jié)果與實(shí)驗(yàn)結(jié)果相符,改變振動(dòng)工況后仍然可以較好地模擬MRD輸出阻尼力的最大值。此外Rational函數(shù)模型可以修正多項(xiàng)式函數(shù)模型高速段不收斂的震蕩現(xiàn)象,對(duì)同一輸入電流不同振動(dòng)工況下的適用性較強(qiáng)。

    2種模型的非線性區(qū)均存在一定程度誤差的原因是:不同振動(dòng)工況(正弦輸入信號(hào))直接影響模型的輸入速度,所以簡化模型中各參數(shù)與正弦輸入信號(hào)的關(guān)系,在同一振動(dòng)工況下建立各參數(shù)與電流之間的絕對(duì)關(guān)系,使得改變振動(dòng)工況后非線性低速區(qū)的仿真結(jié)果存在一定程度的誤差。

    5 模型誤差對(duì)比評(píng)估

    2種非線性模型與實(shí)驗(yàn)結(jié)果之間的誤差可通過均方根誤差RMSE來進(jìn)行評(píng)價(jià),計(jì)算公式如下:

    式中:RMSE為均方根誤差;yi為擬合數(shù)據(jù)為實(shí)驗(yàn)數(shù)據(jù)。

    5.1 辨識(shí)工況

    對(duì)同一振動(dòng)工況(A=2 mm,f=3 Hz)不同電流下,2種模型的仿真結(jié)果與實(shí)驗(yàn)結(jié)果各段的均方根誤差值及Rational函數(shù)與多項(xiàng)式函數(shù)相比Rational函數(shù)的誤差降低率,如表3所示。

    CRMSE(累積均方根誤差)計(jì)算公式表示為:

    ERR(誤差降低率)計(jì)算公式表示為:

    表3 2種模型在不同電流下的均方根誤差

    通過對(duì)表3進(jìn)行分析,可以得出以下結(jié)果:

    1)2種模型在同一振動(dòng)工況不同電流下各段的擬合精度平均可達(dá)0.99以上;

    2)不同電流下,Rational函數(shù)模型的累積擬合誤差與多項(xiàng)式函數(shù)模型的累積擬合誤差相比均有一定程度的降低,與2種模型各自的仿真結(jié)果一致;

    3)在電流低于0.5 A時(shí),Rational函數(shù)模型的誤差降低率隨著電流的增大而增大,電流高于0.5 A后Rational函數(shù)模型的誤差降低率隨著電流的增大而減小。

    5.2 驗(yàn)證工況

    對(duì)同一電流(I=1.0 A)不同振動(dòng)工況(分別改變正弦位移激勵(lì)輸入信號(hào)的頻率和振幅)下,2種模型的仿真結(jié)果與實(shí)驗(yàn)結(jié)果各段的均方根誤差值及Rational函數(shù)模型與多項(xiàng)式函數(shù)模型相比Rational函數(shù)模型的誤差降低率,如表4所示。

    通過對(duì)表4進(jìn)行分析,可以得出以下結(jié)果:

    1)2種模型線性區(qū)的誤差結(jié)果較為接近,與模型結(jié)構(gòu)相符合;非線性區(qū)多項(xiàng)式函數(shù)模型的誤差與Rational函數(shù)模型相比均有不同程度的增大,與仿真結(jié)果一致;

    2)同一電流不同振動(dòng)工況下Rational函數(shù)模型與多項(xiàng)式函數(shù)模型相比Rational函數(shù)模型的誤差降低率均在20%以上,結(jié)合圖11、圖12的仿真結(jié)果可知,Rational函數(shù)模型的適用性強(qiáng)于多項(xiàng)式函數(shù)模型。

    表4 2種模型在不同振動(dòng)工況下的均方根誤差

    6 結(jié)論

    1)相比于多項(xiàng)式函數(shù)模型,Rational函數(shù)模型的分子、分母項(xiàng)均為一次線性函數(shù),模型參數(shù)較少,且模型結(jié)構(gòu)簡單,分段方式明確。

    2)在相同振動(dòng)工況下,Rational函數(shù)模型可改進(jìn)多項(xiàng)式函數(shù)模型在電流低于0.5 A時(shí)高速段存在的不穩(wěn)定性和在速度換向拐點(diǎn)處擬合效果不佳的現(xiàn)象,之后隨著電流增大,2種模型擬合誤差的差距逐漸縮小。Rational函數(shù)模型整體上可更加精確地描述MRD在不同電流下的非線性動(dòng)態(tài)特性。

    3)在相同電流下,改變輸入正弦位移激勵(lì)的頻率或振幅,通過不同頻率、不同振幅激勵(lì)工況的驗(yàn)證表明,Rational函數(shù)模型與多項(xiàng)式函數(shù)模型相比,具有更高的精度和更好的適用性。

    猜你喜歡
    剪切應(yīng)力建模振動(dòng)
    振動(dòng)的思考
    聯(lián)想等效,拓展建?!浴皫щ娦∏蛟诘刃?chǎng)中做圓周運(yùn)動(dòng)”為例
    振動(dòng)與頻率
    心瓣瓣膜區(qū)流場(chǎng)中湍流剪切應(yīng)力對(duì)瓣膜損害的研究進(jìn)展
    基于PSS/E的風(fēng)電場(chǎng)建模與動(dòng)態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    不對(duì)稱半橋變換器的建模與仿真
    中立型Emden-Fowler微分方程的振動(dòng)性
    剪切應(yīng)力對(duì)聚乳酸結(jié)晶性能的影響
    中國塑料(2016年6期)2016-06-27 06:34:24
    動(dòng)脈粥樣硬化病變進(jìn)程中血管細(xì)胞自噬的改變及低剪切應(yīng)力對(duì)血管內(nèi)皮細(xì)胞自噬的影響*
    硫化氫在低剪切應(yīng)力導(dǎo)致內(nèi)皮細(xì)胞自噬障礙中的作用
    国产精品麻豆人妻色哟哟久久 | 国产成人精品婷婷| 伦理电影大哥的女人| 午夜精品国产一区二区电影 | 亚洲中文字幕日韩| 大香蕉久久网| 日韩欧美在线乱码| 九草在线视频观看| 国产精品久久电影中文字幕| 日韩视频在线欧美| 99热这里只有是精品在线观看| 嘟嘟电影网在线观看| 精品久久久久久成人av| 中文字幕亚洲精品专区| 亚洲欧美中文字幕日韩二区| 最后的刺客免费高清国语| 久久精品久久精品一区二区三区| 99九九线精品视频在线观看视频| 亚洲av成人av| 99久久中文字幕三级久久日本| 日本色播在线视频| 日韩欧美精品v在线| 亚洲真实伦在线观看| 久久这里只有精品中国| 亚洲五月天丁香| 嘟嘟电影网在线观看| 尾随美女入室| 日本色播在线视频| 日日干狠狠操夜夜爽| 午夜激情福利司机影院| 国产精品野战在线观看| 成人毛片a级毛片在线播放| 99久久中文字幕三级久久日本| av卡一久久| 国产单亲对白刺激| 日本猛色少妇xxxxx猛交久久| 国产精品无大码| 男女国产视频网站| 午夜激情福利司机影院| 麻豆成人av视频| 亚洲精品亚洲一区二区| 最近视频中文字幕2019在线8| 亚洲欧美中文字幕日韩二区| 中文字幕av在线有码专区| 人妻夜夜爽99麻豆av| 亚州av有码| 国产精品久久视频播放| 久久精品国产亚洲av天美| 九草在线视频观看| 99热网站在线观看| 日本欧美国产在线视频| 99热精品在线国产| 亚洲精品久久久久久婷婷小说 | 晚上一个人看的免费电影| 在线播放无遮挡| 日韩精品青青久久久久久| 日本猛色少妇xxxxx猛交久久| 久久久久网色| 99久久精品一区二区三区| 国产老妇女一区| 亚洲综合色惰| 日本欧美国产在线视频| 小蜜桃在线观看免费完整版高清| 久久久午夜欧美精品| 男女国产视频网站| 桃色一区二区三区在线观看| 国产探花极品一区二区| 欧美日韩在线观看h| 久久久久久国产a免费观看| 亚洲国产色片| 久久亚洲国产成人精品v| 波多野结衣高清无吗| 色综合站精品国产| 成人鲁丝片一二三区免费| 欧美日韩国产亚洲二区| 91aial.com中文字幕在线观看| 久久久精品欧美日韩精品| 三级国产精品欧美在线观看| 久久精品国产自在天天线| 日本色播在线视频| 国内揄拍国产精品人妻在线| 国产av在哪里看| 久久精品人妻少妇| 亚洲国产精品sss在线观看| 亚洲成人中文字幕在线播放| 亚洲欧洲国产日韩| 99久久中文字幕三级久久日本| 97超视频在线观看视频| 久久久久久大精品| 一级二级三级毛片免费看| 亚洲av熟女| 午夜老司机福利剧场| 久99久视频精品免费| 欧美3d第一页| 国产日韩欧美在线精品| 美女脱内裤让男人舔精品视频| 久久久久久伊人网av| 日本一二三区视频观看| 国产探花在线观看一区二区| 日韩欧美精品v在线| av在线播放精品| 成人亚洲精品av一区二区| 九草在线视频观看| 婷婷六月久久综合丁香| 午夜精品在线福利| 亚洲欧美一区二区三区国产| 成人无遮挡网站| 精品99又大又爽又粗少妇毛片| 免费黄色在线免费观看| 亚洲精品自拍成人| 久久精品人妻少妇| 黄片无遮挡物在线观看| 欧美日本亚洲视频在线播放| 26uuu在线亚洲综合色| 在线免费观看的www视频| 国产高清视频在线观看网站| 男女边吃奶边做爰视频| 99国产精品一区二区蜜桃av| 能在线免费看毛片的网站| 免费在线观看成人毛片| 国产单亲对白刺激| 大又大粗又爽又黄少妇毛片口| 欧美区成人在线视频| 亚洲综合色惰| 亚洲欧美日韩高清专用| 国产成人一区二区在线| www.色视频.com| 春色校园在线视频观看| 欧美潮喷喷水| 午夜福利在线在线| 亚洲高清免费不卡视频| 亚洲av成人精品一区久久| 国产女主播在线喷水免费视频网站 | 久久久久久久国产电影| 亚洲成人av在线免费| 97超视频在线观看视频| 韩国高清视频一区二区三区| h日本视频在线播放| 亚洲一级一片aⅴ在线观看| 精品一区二区三区视频在线| www.av在线官网国产| 国产伦精品一区二区三区视频9| 国模一区二区三区四区视频| 波多野结衣巨乳人妻| 国产一区二区三区av在线| 狠狠狠狠99中文字幕| 国产久久久一区二区三区| 国产伦精品一区二区三区四那| 国产精品电影一区二区三区| 69人妻影院| 色噜噜av男人的天堂激情| 色综合亚洲欧美另类图片| 国产一区二区亚洲精品在线观看| 男女视频在线观看网站免费| 久热久热在线精品观看| 99热精品在线国产| 国语自产精品视频在线第100页| 国产精品国产三级专区第一集| 久久久久久久久大av| 国产男人的电影天堂91| 一边亲一边摸免费视频| 一个人免费在线观看电影| 校园人妻丝袜中文字幕| 亚洲va在线va天堂va国产| 久久99热6这里只有精品| 晚上一个人看的免费电影| 亚洲av中文字字幕乱码综合| 国产精品国产高清国产av| 国产精品一区二区性色av| 亚洲真实伦在线观看| 黑人高潮一二区| 青春草亚洲视频在线观看| 免费搜索国产男女视频| 99热全是精品| 精品国产一区二区三区久久久樱花 | 亚洲自偷自拍三级| 99久国产av精品国产电影| 亚洲三级黄色毛片| 黄片无遮挡物在线观看| 亚洲精品国产成人久久av| 一区二区三区乱码不卡18| 女人被狂操c到高潮| 国产淫语在线视频| 在线观看66精品国产| 国产精品国产三级专区第一集| 亚洲精品乱久久久久久| 一区二区三区乱码不卡18| 免费看av在线观看网站| 尾随美女入室| 国产精品国产三级专区第一集| 日韩三级伦理在线观看| 国产成人a∨麻豆精品| 三级毛片av免费| 色尼玛亚洲综合影院| 国产伦理片在线播放av一区| 久久99精品国语久久久| 成人漫画全彩无遮挡| 国产亚洲精品久久久com| 久久精品久久久久久久性| 在线播放国产精品三级| 日韩一区二区三区影片| 色5月婷婷丁香| 色尼玛亚洲综合影院| 国产真实乱freesex| 国产精品美女特级片免费视频播放器| 中文在线观看免费www的网站| 国产精品1区2区在线观看.| 少妇人妻一区二区三区视频| 中文字幕亚洲精品专区| 亚洲熟妇中文字幕五十中出| a级毛片免费高清观看在线播放| 狠狠狠狠99中文字幕| 乱系列少妇在线播放| 欧美丝袜亚洲另类| 国产乱人偷精品视频| 欧美变态另类bdsm刘玥| 国产成人freesex在线| 国内精品美女久久久久久| 黄片wwwwww| 九九热线精品视视频播放| 插阴视频在线观看视频| 国产一级毛片在线| 亚洲国产成人一精品久久久| 真实男女啪啪啪动态图| 亚洲人成网站在线观看播放| 久久精品国产亚洲av涩爱| 日韩精品青青久久久久久| 久久精品影院6| 最近手机中文字幕大全| videos熟女内射| 美女脱内裤让男人舔精品视频| 99久久中文字幕三级久久日本| 成人漫画全彩无遮挡| 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| 校园人妻丝袜中文字幕| 一个人观看的视频www高清免费观看| 亚洲人与动物交配视频| 观看免费一级毛片| 一个人看视频在线观看www免费| 中文天堂在线官网| 亚洲欧美精品专区久久| 久久99精品国语久久久| 黄片无遮挡物在线观看| 日本五十路高清| 国产高潮美女av| 国产精品一区二区性色av| 国内少妇人妻偷人精品xxx网站| 亚洲精品亚洲一区二区| 日韩欧美精品v在线| 少妇被粗大猛烈的视频| 国产三级中文精品| 成人性生交大片免费视频hd| 成人午夜高清在线视频| 亚洲最大成人av| 国产 一区精品| 亚洲中文字幕一区二区三区有码在线看| 亚洲电影在线观看av| 少妇的逼好多水| 欧美极品一区二区三区四区| 男的添女的下面高潮视频| 嫩草影院入口| 亚洲,欧美,日韩| 国产精品人妻久久久影院| 午夜激情欧美在线| 日韩视频在线欧美| 日韩欧美三级三区| 少妇人妻精品综合一区二区| 亚洲人成网站高清观看| 波野结衣二区三区在线| 久久欧美精品欧美久久欧美| 国产国拍精品亚洲av在线观看| 淫秽高清视频在线观看| 国产亚洲精品久久久com| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产精品国产精品| 日韩大片免费观看网站 | 久久99热这里只有精品18| 丰满乱子伦码专区| 亚洲国产色片| videos熟女内射| 亚洲第一区二区三区不卡| 美女国产视频在线观看| 亚洲高清免费不卡视频| 中文字幕免费在线视频6| 网址你懂的国产日韩在线| 国产精品福利在线免费观看| 精品一区二区免费观看| 春色校园在线视频观看| 欧美变态另类bdsm刘玥| 欧美日韩在线观看h| 色播亚洲综合网| 青春草视频在线免费观看| 成人毛片a级毛片在线播放| 国国产精品蜜臀av免费| 超碰97精品在线观看| 国产成人精品久久久久久| 国产成人a∨麻豆精品| 国产伦理片在线播放av一区| 麻豆国产97在线/欧美| 亚洲天堂国产精品一区在线| 亚洲精品成人久久久久久| 午夜激情欧美在线| 久久久亚洲精品成人影院| 春色校园在线视频观看| 最近中文字幕2019免费版| 久久99热这里只频精品6学生 | 日韩大片免费观看网站 | 久久99精品国语久久久| 国产中年淑女户外野战色| 九九在线视频观看精品| 亚洲国产最新在线播放| www.av在线官网国产| 国产精品精品国产色婷婷| 久久人妻av系列| 国产午夜精品久久久久久一区二区三区| 禁无遮挡网站| 精品国产一区二区三区久久久樱花 | 日本三级黄在线观看| 亚洲欧美日韩高清专用| 午夜精品国产一区二区电影 | av播播在线观看一区| 国产成人a∨麻豆精品| 国产片特级美女逼逼视频| 99热网站在线观看| 我要搜黄色片| 精品国产三级普通话版| 午夜福利在线观看吧| 国产成人a区在线观看| 亚洲精品久久久久久婷婷小说 | 性色avwww在线观看| a级一级毛片免费在线观看| 2021天堂中文幕一二区在线观| av视频在线观看入口| 日本一本二区三区精品| 人人妻人人澡欧美一区二区| 美女内射精品一级片tv| 99热这里只有是精品50| 只有这里有精品99| 91精品伊人久久大香线蕉| 国产淫语在线视频| 精品久久久久久久久亚洲| 午夜老司机福利剧场| 免费黄网站久久成人精品| 麻豆成人午夜福利视频| 久久久久久久国产电影| 九草在线视频观看| 国产成人免费观看mmmm| 精品人妻一区二区三区麻豆| 高清日韩中文字幕在线| 精品久久久久久久久久久久久| 国产真实伦视频高清在线观看| 三级国产精品片| 久久久亚洲精品成人影院| 中文乱码字字幕精品一区二区三区 | 日韩欧美国产在线观看| 我的女老师完整版在线观看| 老师上课跳d突然被开到最大视频| 少妇丰满av| 嘟嘟电影网在线观看| 村上凉子中文字幕在线| 成年免费大片在线观看| 只有这里有精品99| 成人二区视频| 3wmmmm亚洲av在线观看| 蜜桃久久精品国产亚洲av| 欧美bdsm另类| 国产v大片淫在线免费观看| 亚洲精品乱码久久久久久按摩| 一本一本综合久久| 国产精品一区二区性色av| 日本欧美国产在线视频| 亚洲欧美精品自产自拍| av女优亚洲男人天堂| 天堂中文最新版在线下载 | 日韩大片免费观看网站 | 一级毛片久久久久久久久女| 国产黄色视频一区二区在线观看 | 少妇丰满av| 亚洲欧美日韩卡通动漫| 久久久久性生活片| 麻豆成人av视频| 丰满人妻一区二区三区视频av| 一级黄片播放器| h日本视频在线播放| 久久精品综合一区二区三区| 中文字幕免费在线视频6| 久久热精品热| 性插视频无遮挡在线免费观看| 久久精品熟女亚洲av麻豆精品 | 亚洲av日韩在线播放| 深爱激情五月婷婷| 国产综合懂色| 国产 一区精品| 亚洲欧美中文字幕日韩二区| 亚洲av二区三区四区| 超碰av人人做人人爽久久| 又爽又黄a免费视频| 精品人妻偷拍中文字幕| 中国美白少妇内射xxxbb| 日韩精品青青久久久久久| 能在线免费看毛片的网站| 欧美日韩国产亚洲二区| 亚洲精品色激情综合| 国产精品熟女久久久久浪| 免费av观看视频| 国产一区亚洲一区在线观看| 99久久成人亚洲精品观看| 精品国产三级普通话版| 国产亚洲一区二区精品| 午夜福利高清视频| 欧美成人一区二区免费高清观看| 熟妇人妻久久中文字幕3abv| 少妇丰满av| 毛片一级片免费看久久久久| 欧美丝袜亚洲另类| 国产一区有黄有色的免费视频 | 欧美成人精品欧美一级黄| 夜夜看夜夜爽夜夜摸| 国模一区二区三区四区视频| 18禁在线无遮挡免费观看视频| 成年免费大片在线观看| 亚洲国产日韩欧美精品在线观看| 久久精品久久久久久噜噜老黄 | 午夜福利在线观看吧| 99久国产av精品| 亚洲丝袜综合中文字幕| 国产极品精品免费视频能看的| 国产免费一级a男人的天堂| 免费av毛片视频| 日韩高清综合在线| 国产成人精品一,二区| 99久久人妻综合| 22中文网久久字幕| 性插视频无遮挡在线免费观看| 国产亚洲午夜精品一区二区久久 | 美女xxoo啪啪120秒动态图| 国产色婷婷99| 久久婷婷人人爽人人干人人爱| 精品欧美国产一区二区三| 不卡视频在线观看欧美| 国产精品不卡视频一区二区| 亚洲国产精品sss在线观看| 国产黄色视频一区二区在线观看 | 两个人视频免费观看高清| 久久这里有精品视频免费| 亚洲av中文av极速乱| 日韩一区二区三区影片| 国产高清不卡午夜福利| 国产成人福利小说| 亚洲欧美精品综合久久99| 大话2 男鬼变身卡| 一边摸一边抽搐一进一小说| 欧美高清成人免费视频www| 色噜噜av男人的天堂激情| 日本欧美国产在线视频| 久久亚洲国产成人精品v| 亚洲av男天堂| 亚洲av不卡在线观看| 男人舔奶头视频| 国产精品久久久久久精品电影小说 | 免费av不卡在线播放| 性插视频无遮挡在线免费观看| 日本黄色视频三级网站网址| 国产真实乱freesex| 自拍偷自拍亚洲精品老妇| 女的被弄到高潮叫床怎么办| 亚洲美女视频黄频| 最近最新中文字幕大全电影3| 久久人人爽人人爽人人片va| 亚洲国产日韩欧美精品在线观看| 欧美成人a在线观看| 国内精品美女久久久久久| 久久精品夜色国产| 亚洲精品成人久久久久久| 亚洲久久久久久中文字幕| av国产免费在线观看| 成人av在线播放网站| 精品一区二区三区视频在线| 伦精品一区二区三区| 天堂√8在线中文| 男女那种视频在线观看| 久久99热这里只有精品18| 婷婷色av中文字幕| 国产精品久久久久久久电影| 国产成人a∨麻豆精品| 久久久久免费精品人妻一区二区| 亚洲精品乱久久久久久| 亚洲图色成人| 欧美日本亚洲视频在线播放| 简卡轻食公司| 欧美精品一区二区大全| 99在线视频只有这里精品首页| 你懂的网址亚洲精品在线观看 | av播播在线观看一区| 国产人妻一区二区三区在| 高清午夜精品一区二区三区| 欧美成人a在线观看| 伦理电影大哥的女人| 岛国在线免费视频观看| 尾随美女入室| 91久久精品国产一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 99久久成人亚洲精品观看| 国产三级在线视频| 麻豆av噜噜一区二区三区| 91aial.com中文字幕在线观看| 国产麻豆成人av免费视频| 日本黄色视频三级网站网址| 丝袜美腿在线中文| 搞女人的毛片| 亚洲av一区综合| 尾随美女入室| 免费看美女性在线毛片视频| 午夜免费男女啪啪视频观看| 亚洲精品国产成人久久av| 免费av观看视频| 卡戴珊不雅视频在线播放| 国产午夜福利久久久久久| 国产亚洲一区二区精品| 亚洲怡红院男人天堂| 成年av动漫网址| 国产黄a三级三级三级人| 搞女人的毛片| 国产精品嫩草影院av在线观看| 亚洲av中文av极速乱| 亚洲国产精品久久男人天堂| 亚洲国产欧美人成| 亚洲,欧美,日韩| 精品久久久久久久末码| 欧美又色又爽又黄视频| 国产精品一二三区在线看| 欧美日本亚洲视频在线播放| 亚洲国产精品成人久久小说| 久久久久性生活片| 日本黄色视频三级网站网址| 成人无遮挡网站| 欧美日韩国产亚洲二区| 国产亚洲午夜精品一区二区久久 | 亚州av有码| 观看免费一级毛片| 亚洲精品色激情综合| 精品人妻一区二区三区麻豆| 黄片wwwwww| 亚洲一区高清亚洲精品| 婷婷色麻豆天堂久久 | 大话2 男鬼变身卡| 欧美成人一区二区免费高清观看| 欧美三级亚洲精品| 最近手机中文字幕大全| 91精品伊人久久大香线蕉| 九九在线视频观看精品| 少妇猛男粗大的猛烈进出视频 | 久热久热在线精品观看| 婷婷六月久久综合丁香| av线在线观看网站| 99热网站在线观看| 欧美一级a爱片免费观看看| 国产精品日韩av在线免费观看| 国产伦理片在线播放av一区| 免费观看人在逋| 日韩,欧美,国产一区二区三区 | 美女被艹到高潮喷水动态| 亚洲av男天堂| 别揉我奶头 嗯啊视频| 日韩制服骚丝袜av| 九九爱精品视频在线观看| 欧美bdsm另类| 亚洲av电影不卡..在线观看| 国产一区二区三区av在线| 久热久热在线精品观看| 国产伦精品一区二区三区视频9| 丝袜喷水一区| 国产成人精品婷婷| 日韩在线高清观看一区二区三区| 高清av免费在线| 乱码一卡2卡4卡精品| 国产精品av视频在线免费观看| 亚洲图色成人| 日韩欧美三级三区| 欧美97在线视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产淫语在线视频| 久久久久久大精品| 中国美白少妇内射xxxbb| 国产免费视频播放在线视频 | 青春草视频在线免费观看| 精品午夜福利在线看| 国产色婷婷99| 亚洲av.av天堂| 亚洲在线自拍视频| 九色成人免费人妻av| 干丝袜人妻中文字幕| 人人妻人人澡欧美一区二区| 一本一本综合久久| 欧美又色又爽又黄视频| 亚洲欧美成人综合另类久久久 | 亚洲av日韩在线播放| 99热6这里只有精品| 汤姆久久久久久久影院中文字幕 | 少妇丰满av| 高清毛片免费看| 久久久久久伊人网av| 国产精品久久电影中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 舔av片在线| 亚洲国产精品合色在线| 欧美xxxx黑人xx丫x性爽| 国产免费一级a男人的天堂| 三级国产精品片| 亚洲成人中文字幕在线播放| 蜜桃久久精品国产亚洲av| 免费不卡的大黄色大毛片视频在线观看 | 六月丁香七月| 免费观看在线日韩|