999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

伺服驅(qū)動壓邊力控制系統(tǒng)建模與仿真

2015-12-30 05:15:14楊莉,王國佐,鄒鳳剛
中國機械工程 2015年19期
關鍵詞:系統(tǒng)仿真

伺服驅(qū)動壓邊力控制系統(tǒng)建模與仿真

楊莉王國佐鄒鳳剛秦泗吉

燕山大學,秦皇島,066004

摘要:將伺服驅(qū)動技術應用于壓邊力控制,可以更有效地控制拉深成形過程,提高成形效果。采用伺服電機驅(qū)動的壓邊力控制系統(tǒng)由機械和數(shù)控系統(tǒng)兩部分組成,系統(tǒng)參數(shù)匹配的合理性對壓邊力控制的動態(tài)性能影響較大,而系統(tǒng)的動態(tài)性能直接影響到壓邊力控制精度和效果。首先對雙滑塊六桿機構輸入輸出方程進行了線性化處理,建立了機械系統(tǒng)、數(shù)控伺服系統(tǒng)的數(shù)學模型,得到了壓邊力控制系統(tǒng)的數(shù)學模型和仿真模型,對控制系統(tǒng)進行了仿真分析,初步確定了控制系統(tǒng)的位置環(huán)PID參數(shù)。最后根據(jù)仿真分析結果進行了壓邊力控制實驗系統(tǒng)的參數(shù)設定,實驗表明,根據(jù)仿真值稍作調(diào)整就能達到實驗系統(tǒng)的控制要求,大大縮短了系統(tǒng)調(diào)試時間。

關鍵詞:壓邊力控制;伺服驅(qū)動;數(shù)學模型;PID控制;系統(tǒng)仿真

中圖分類號:TG385.9

收稿日期:2014-11-26

基金項目:國家自然科學基金資助項目(51175451)

作者簡介:楊莉,女,1962年生。燕山大學機械工程學院副教授、博士。研究方向為先進制造、數(shù)控重載驅(qū)動技術。王國佐,男,1989年生。燕山大學機械工程學院碩士研究生。鄒鳳剛,男,1988年生。燕山大學機械工程學院碩士研究生。秦泗吉,男,1963年生。燕山大學機械工程學院教授、博士。

Modeling and Simulation of Servo-driven BHF Control System

Yang LiWang GuozuoZou FenggangQin Siji

Yanshan University,Qinhuangdao,Hebei,066004

Abstract:Servo-driven control technology of BHF is more effective and can improve the forming effect in deep drawing process. BHF control system driven by servo motor was composed of the mechanical system and CNC system. The parameters of components in BHF system had a greater influence on the dynamic performance of the system, so it could have effects on accuracy and efficiency in BHF control. Firstly linearization of the input and output equations was made for the six-bar linkage and the double-slider mechanism. And mathematical models of mechanical system and CNC control system were obtained respectively. Then performance of the system was simulated and analyzed and the PID parameters of the position loop were adjusted. Finally, the parameters in the experiments of the BHF control system were given based on simulation results. The experimental results show that it can meet the requirements of the BHF control to adjust slightly the values of the parameters of the control components based on the modeling and simulation results, so the debug time of system control parameters can be shorten greatly.

Key words: blank holder force(BHF) control; servo drive; mathematical model; PID control; system simulation

0引言

拉深成形是板材成形中的典型工藝之一,在

生產(chǎn)中的應用非常廣泛。在成形過程中,通常將控制壓邊力作為控制成形過程的主要手段[1]。板坯在成形過程中的各個階段所需要的壓邊力是不同的,如果在各變形階段都施加相同的壓邊力,就可能出現(xiàn)在某個變形階段因壓邊力太小而使得起皺趨勢嚴重,而在另外的變形階段壓邊力過大導致某個部位出現(xiàn)破裂的危險[2]。在變壓邊力控制過程中,壓邊力參數(shù)的合理取值、調(diào)整和控制,成為控制金屬流動、防止出現(xiàn)起皺和破裂的最重要的手段之一[3-5]。

在理想情況下,拉深成形過程中的壓邊力是隨行程和位置同時變化的。但目前在生產(chǎn)實際中,采用通用壓力機或?qū)S脡毫C的成形,一般都不能滿足這樣的工藝要求,所以開發(fā)與拉深工藝要求相適應的、隨行程和位置同時變化的變壓邊力控制技術是非常必要的。

液壓傳動方式具有傳力大、易于控制等優(yōu)點,可用于大中型壓力機上的變壓邊力控制。變壓邊力控制系統(tǒng)是由比例溢流閥控制液壓缸內(nèi)的壓力而實現(xiàn)控制過程的[6]。這種控制系統(tǒng)結構復雜、響應時間慢、能耗高。

隨著控制技術的不斷發(fā)展,特別是計算機及其相關技術的發(fā)展,數(shù)控技術在成形領域得到了越來越多的應用。根據(jù)拉深工藝對壓邊力的要求,文獻[7]提出了以伺服電機作為驅(qū)動元件實現(xiàn)數(shù)控壓邊的方法。這種壓邊力控制方法采用伺服電機驅(qū)動,配以適當?shù)膲哼吜?zhí)行機構,對壓邊過程中壓料板的位置、速度、力等進行精確控制,使壓邊力按預先設定的規(guī)律隨拉深行程變化,實現(xiàn)不同的拉深工藝要求,達到控制目的。

筆者在前期工作[8-9]中,主要對應用于不同工藝形式的壓邊力執(zhí)行機構進行了研究,確定了機構構型和桿系尺寸。在此基礎上,還需要對整個控制系統(tǒng)進行建模和仿真,分析系統(tǒng)的性能,以便構建系統(tǒng)的實驗平臺。

伺服驅(qū)動壓邊力控制系統(tǒng)由數(shù)控系統(tǒng)和機械系統(tǒng)組成,這兩部分系統(tǒng)參數(shù)的合理匹配是整個控制系統(tǒng)得以良好運行的保證。本文以所研制的壓邊力控制系統(tǒng)為基礎,對壓邊力控制系統(tǒng)進行了建模、仿真,確定了系統(tǒng)控制器參數(shù),為實驗系統(tǒng)的PID參數(shù)設定提供了依據(jù),這種方法可顯著縮短系統(tǒng)調(diào)試時間。

1壓邊力控制系統(tǒng)的數(shù)學模型

1.1系統(tǒng)組成及原理

圖1為伺服驅(qū)動壓邊力控制系統(tǒng)組成框圖,其中數(shù)控系統(tǒng)包括上位機PC、下位機運動控制器、伺服電機、驅(qū)動器、位移和力檢測裝置等。機械系統(tǒng)由聯(lián)軸器、絲杠螺母副、六桿機構、壓料板、緩沖件和連接件等組成(圖2)。PC通過標準PCI總線連接具有實時控制功能的運動控制器,組成完整的數(shù)控裝置,數(shù)控裝置發(fā)出指令信號經(jīng)伺服驅(qū)動器轉換和放大后使電機轉動,伺服電機驅(qū)動絲杠螺母副和六桿機構帶動壓料板實現(xiàn)給定的壓邊運動,電機內(nèi)置脈沖編碼器用于檢測電機轉角和轉速,并實現(xiàn)反饋,力檢測裝置用于壓邊力的測量。

圖1 伺服驅(qū)動壓邊力控制系統(tǒng)組成框圖

1.伺服電機 2.聯(lián)軸器 3.下滑板 4.滾珠絲杠 5.六桿機構 6.螺母 7.凹模 8.測力環(huán) 9.壓料板組件 10.凸模 11.墊板 12.連接桿 13.緩沖件 14.連接板 圖2 伺服驅(qū)動壓邊力控制機械系統(tǒng)示意圖

1.2雙滑塊六桿機構數(shù)學模型

1.2.1雙滑塊六桿機構輸入和輸出關系

如圖2所示,六桿機構以螺母的勻速運動為輸入,以壓料板的變速直線運動為輸出,其輸入與輸出運動關系可按圖3所示的六桿機構簡圖以及給定相應參數(shù)導出。

圖3 六桿機構簡圖

如圖3所示,輸入滑塊為螺母,輸出滑塊為下滑板。可以認為,六桿機構由兩個四桿機構組成,上滑塊、桿AB、桿OB和機架組成一四桿機構,其中桿AB與上滑塊(輸入滑塊)相連。桿BC、桿CD、下滑塊和機架組成另一四桿機構,其中桿CD與下滑塊(輸出滑塊)相連,圖中桿OB和桿BC為剛性連接。

設桿OB、桿BA、桿CD與x軸夾角分別為θ1、θ2、θ3,OB、BA、BC、CD對應的各桿長度尺寸分別為l1、l2、l3、l4,輸入滑塊與固定鉸接點O的水平距離為a。圖3中x1、x2分別為輸入和輸出滑塊至點O的垂直距離,它們與θ1之間的關系為

(1)

(2)

六桿機構輸入與輸出關系可由以θ1為參數(shù)的式 (1)和式 (2)確定,其輸入與輸出關系為非線性關系,圖4所示為六桿機構輸入與輸出關系曲線。

1.2.2雙滑塊六桿機構的線性化處理

因六桿機構的輸入與輸出關系是非線性的,為了簡化建模和仿真,需進行線性化處理。

若輸入為x1、輸出為x2,設輸入輸出關系為

x2=f(x1)

(3)

若給定輸入點位置x10,對應輸出點位置x20,將式(3)展開成泰勒級數(shù):

(4)

(5)

(6)

對正裝模具結構,因壓邊力施加階段的工作行程很小,可在工作點附近,將六桿機構的輸入和輸出關系作為線性環(huán)節(jié)來處理。略去式(4)中(x1-x10)2項及其后面的高階項,得增量方程Δx2=kdΔx1,為簡化起見,省去增量符號Δ,增量方程改寫為[10]

x2=kdx1

(7)

式中,kd為比例系數(shù)。

其傳遞函數(shù)為

(8)

當壓邊工作點位置不同時,kd也不相同。

若壓邊行程在100~130mm范圍內(nèi),其對應的輸入與輸出二階導數(shù)關系曲線如圖5所示。隨著壓邊行程的增大,二階導數(shù)的絕對值越來越小,線性化程度也越來越好。若將工作點取x10=124.5mm,計算得出,該位置的輸入與輸出關系二階導數(shù)值小于0.005mm,說明該工作點附近線性化程度較好。

圖5 六桿執(zhí)行機構輸入與輸出二階導數(shù)關系曲線

1.3絲杠螺母副數(shù)學模型

絲杠螺母副以伺服電機的角位移為輸入,以螺母的直線位移為輸出,其動力平衡方程式為

(9)

式中,m為螺母等執(zhí)行部件的質(zhì)量;B1為導軌副上黏性阻尼系數(shù);FL為作用在螺母上的外載荷;Fp為螺母執(zhí)行部件的驅(qū)動力,F(xiàn)p=k(x0-x1);k為折算到絲杠軸上的傳動系統(tǒng)剛度;x1為螺母的直線位移;x0為電機對螺母執(zhí)行部件的輸入位移,x0=Lθ/(2π);θ為電機軸轉角;L為滾珠絲杠導程。

滾珠絲杠副傳動系統(tǒng)的剛度,可認為是與滾珠絲杠相關聯(lián)的零部件剛度的串聯(lián)總和,計算公式如下:

(10)

式中,ks1為絲杠軸拉壓剛度;kb為支承軸承的軸向接觸剛度;kc為滾珠與滾道的接觸剛度;kr為折合到滾珠絲杠副上的伺服剛度;kbr為軸承座的剛度;knr為螺母座的剛度;kt為聯(lián)軸器的剛度;ks2為滾珠絲杠副扭轉剛度折合到執(zhí)行部件的直線剛度。

式(10)中后幾項與前3項相比影響較小,實際計算時取前3項。

對式(9)進行Laplace變換得

(11)

以θ(s)為輸入,以x1(s)為輸出,不考慮外力FL時,其傳遞函數(shù)為

(12)

1.4伺服驅(qū)動系統(tǒng)的數(shù)學模型

壓邊力控制系統(tǒng)選用交流永磁同步伺服電機作為驅(qū)動元件,永磁同步電機在矢量控制條件下與直流電機是等價的,只是將電樞電流換成了與轉子磁鏈正交的交流電流分量iq,永磁同步電機的電壓平衡方程為

(13)

式中,Uq為電機電樞電壓;iq為電樞電流;L0為電樞電感;R為電樞電阻;E為電機反電勢,E=keω;ke為電機電勢系數(shù);ω為轉子角速度。

電機軸上的轉矩平衡方程為

(14)

式中,Mm為電機輸出轉矩,Mm=kTiq;kT為電機轉矩系數(shù);ML為負載轉矩;B2為等效到電機軸上的黏性阻尼系數(shù);J為折算到電機軸上的總轉動慣量。

若不考慮電機負載轉矩,將式(13)、式(14)整理后進行Laplace變換,得

以Uq(s)為輸入,ω(s)為輸出,其傳遞函數(shù)為

(15)

2壓邊力控制系統(tǒng)的仿真

2.1仿真模型

將伺服驅(qū)動系統(tǒng)和機械系統(tǒng)的數(shù)學模型綜合起來,并引入速度環(huán)和位置環(huán),得到圖6所示的壓邊力控制系統(tǒng)的數(shù)學模型,其中,Kv1、Ti分別為速度環(huán)增益和積分時間參數(shù),Kp2、Kv2、Ki2分別為比例、微分、積分增益。

利用仿真技術可以實現(xiàn)對系統(tǒng)動態(tài)性能的分析[11],并找到各參數(shù)間的匹配關系,因此借助動態(tài)系統(tǒng)仿真工具MATLAB/Simulink,對壓邊力控制系統(tǒng)建模,在空載且不計六桿機構與導桿摩擦情況下,系統(tǒng)仿真模型如圖7所示。

圖6 壓邊力控制系統(tǒng)模型

圖7 壓邊力控制系統(tǒng)仿真模型

壓邊力控制系統(tǒng)參數(shù)整定主要在速度環(huán)和位置環(huán)參數(shù)的匹配上,速度環(huán)在伺服驅(qū)動器內(nèi)實現(xiàn),通過電機上的脈沖編碼器反饋電機轉速信號構成,采用的是比例積分控制;位置環(huán)由運動控制器來完成,位置環(huán)是通過電機上的脈沖編碼器反饋電機轉角信號而構成的,采用典型的比例積分微分控制。

2.2仿真參數(shù)

根據(jù)所研制的壓邊力控制系統(tǒng),給定仿真模型的參數(shù)如表1所示,不考慮電機軸的阻尼力矩,選取工作點x10=100.7mm,計算kd=0.05,對壓邊力控制系統(tǒng)性能的影響進行仿真分析。

表1 仿真模型已知參數(shù)

由圖5和上文的分析可見,工作點坐標位置大于所選取的工作點x10時,其輸入輸出關系的線性化程度會更好。

2.3壓邊力控制系統(tǒng)的仿真分析

綜合考慮系統(tǒng)的快速性和穩(wěn)定性要求,根據(jù)表1參數(shù)值和壓邊力控制系統(tǒng)的仿真模型,對壓邊力控制系統(tǒng)進行PID參數(shù)仿真整定,得到整定參數(shù)如表2所示。當系統(tǒng)輸入為單位階躍信號r(t)時,得到系統(tǒng)的階躍響應信號c(t)(圖8),從響應曲線看,系統(tǒng)具有較好的動靜態(tài)性能。

表2 仿真整定參數(shù)

圖8 系統(tǒng)階躍響應信號

當改變工作點位置后,仍可按上述方法進行處理,仿真得到PID參數(shù)。

3壓邊力控制系統(tǒng)的實驗

圖9所示為所研制的數(shù)控系統(tǒng)和機械系統(tǒng)的部分組成,數(shù)控系統(tǒng)由PC、GE-200-SV-PCI運動控制器、SGMGH-20ACA61伺服電機、SGDM-20ADA驅(qū)動器、力檢測裝置及接口電路等組成,其中運動控制器附件端子板10連接運動控制器和伺服驅(qū)動器11,驅(qū)動器11與伺服電機1相連,電機內(nèi)置增量式光電脈沖編碼器。機械部分由六桿機構2和絲杠螺母副3等組成。

(a)機械系統(tǒng)(b)數(shù)控系統(tǒng) 1.伺服電機 2.六桿機構 3.絲杠螺母副 4.凹模 5.測力環(huán) 6.壓料板 7.凸模 8.開關電源 9.伺服變壓器 10.端子板 11.伺服驅(qū)動器 圖9 壓邊力控制實驗系統(tǒng)

將伺服電機的驅(qū)動器設置為速度控制方式,按仿真結果,設置速度環(huán)增益Kv1為40,積分時間常數(shù)Ti為20,仿真值與伺服驅(qū)動器出廠設置一致。

在電機調(diào)整好后,對運動控制器進行調(diào)試,設置±10V的模擬量輸出。不考慮前饋控制,速度前饋增益Kvff、加速度前饋增益Kaff設置為0,主要對比例Kp2、微分Kv2、積分增益Ki2三個PID基本參數(shù)進行設置。

根據(jù)表2仿真整定參數(shù)值,對實驗系統(tǒng)進行調(diào)試,設置壓邊力控制系統(tǒng)各參數(shù)的變化范圍為Kp2(2~8)、Kv2(0~9)、Ki2(12~26)。實驗時先確定Kp2,然后改變Kv2和Ki2,進行壓邊力加載實驗,結果表明Kp2對系統(tǒng)影響最大,當Kp2取8時電機產(chǎn)生振動;改變微分增益Kv2和積分增益Ki2對系統(tǒng)影響不顯著。表3為實驗系統(tǒng)中運動控制器調(diào)節(jié)后的PID參數(shù)值,該參數(shù)下系統(tǒng)的快速性和穩(wěn)定性最佳。

表3 運動控制器PID參數(shù)

仿真得到的PID參數(shù)值,為實驗系統(tǒng)提供了參考依據(jù),有效地縮短了調(diào)試時間。但是機械系統(tǒng)的摩擦、間隙和阻尼會影響系統(tǒng)性能,實驗系統(tǒng)PID參數(shù)與仿真結果相比稍有變化。

4結論

(1)在所選取的工作點,根據(jù)壓邊力的加載行程小和輸入輸出函數(shù)的特點,對六桿執(zhí)行機構的輸入和輸出關系進行了線性化處理。當工作點改變后,選取不同的比例系數(shù)kd,也可作類似的線性化處理。

(2)對壓邊力控制系統(tǒng)的數(shù)控系統(tǒng)、機械系統(tǒng)進行了數(shù)學建模,進一步建立了控制系統(tǒng)的數(shù)學模型和仿真模型,在保證系統(tǒng)具有良好的動態(tài)性能和穩(wěn)態(tài)精度的前提下,仿真得到了PID參數(shù)。

(3)利用仿真得到的PID參數(shù)對壓邊力控制實驗系統(tǒng)進行了參數(shù)設定,仿真為實驗系統(tǒng)PID參數(shù)設定提供了依據(jù),減少了系統(tǒng)調(diào)試時間。

參考文獻:

[1]DemirciHI,YaarM,DemirayK,etal.TheTheoreticalandExperimentalInvestigationofBlankHolderForcesPlateEffectinDeepDrawingProcessofAl1050Material[J].Materials&Design, 2008, 29(2): 526-532.

[2]ZhangWenfeng,ShivpuriR.ProbabilisticDesignofAluminumSheetDrawingforReducedRiskofWrinklingandFracture[J].ReliabilityEngineering&SystemSafety, 2009, 94(2): 152-161.

[3]SeniorBW.FlangeWrinklinginDeep-drawingOperations[J].JournaloftheMechanicsandPhysicsofSolids, 1956, 4(4): 235-246.

[4]SiegertK,DannenmannE,WagnerS,etal.Closed-loopControlSystemforBlankHolderForcesinDeepDrawing[J].CIRPAnnals-ManufacturingTechnology, 1995,44(1): 251-254.

[5]RodriguesDM,LeitaoC,MenezesLF.AMulti-stepAnalysisforDeterminingAdmissibleBlank-holderForcesinDeep-drawingOperations[J].MaterialsandDesign,2010,31:1475-1481.

[6]趙軍. 圓錐形件拉深成形智能化的研究[D]. 哈爾濱:哈爾濱工業(yè)大學,1997.

[7]秦泗吉. 壓邊力控制技術研究現(xiàn)狀及伺服數(shù)控壓邊方法可行性探討[J].中國機械工程,2007, 18(1):120-125.

QinSiji.State-of-the-artofBlankHoldingForceControlTechnologyandFeasibilityofNumericalServo-controlHolding[J].ChinaMechanicalEngineering,2007, 18(1):120-125.

[8]楊莉,秦泗吉,杜廣杰. 復合伺服驅(qū)動壓邊力控制方法及執(zhí)行機構設計[J].中國機械工程,2012, 23(6):1360-1363.

YangLi,QinSiji,DuGuangjie.ACompositeControlTechniqueofBHFDrivenbyServeoMotorandDesigmfortheActuator[J].ChinaMechanicalEngineering,2012, 23(6):1360-1363.

[9]楊莉,金振林,秦泗吉,等. 復合伺服驅(qū)動壓邊力控制執(zhí)行機構建模與優(yōu)化[J].機械工程學報,2015,51(3):139-145.

YangLi,JinZhenlin,QinSiji,etal.ModelingandOptimizationofCompositeandServo-drivenActuatorforBlank-holderForceControl[J].JournalofMechanicalEngineering,2015,51(3):139-145.

[10]鄒伯敏.自動控制理論[M].北京:機械工業(yè)出版社,1999.

[11]EbrahimiM,WhalleyR.Analysis,ModelingandSimulationofStiffnessinMachineToolDrives[J].Computers&IndustrialEngineering, 2000, 38(1): 93-105.

(編輯郭偉)

猜你喜歡
系統(tǒng)仿真
相控陣雷達系統(tǒng)仿真模型研究
基于Flexsim的集裝箱碼頭物流作業(yè)系統(tǒng)建模與仿真
物流科技(2016年12期)2017-04-01 03:06:34
油氣儲運專業(yè)仿真系統(tǒng)的開發(fā)和應用
科教導刊(2016年33期)2017-02-27 11:00:52
重大交通設施項目風險復雜動態(tài)交互演化機理與仿真分析
預測(2016年3期)2016-12-29 18:57:28
基于前饋補償?shù)闹鄙龣C控制律設計與仿真
航空兵器(2016年4期)2016-11-28 21:34:37
基于計算機系統(tǒng)仿真的電力機車模擬駕駛裝置研究
苗圃大棚自動灌溉控制系統(tǒng)的設計
大滯后系統(tǒng)中單神經(jīng)元PID控制器設計
雙饋感應發(fā)電機不脫網(wǎng)運行的系統(tǒng)仿真
小型民用風力發(fā)電系統(tǒng)最大功率跟蹤策略研究
主站蜘蛛池模板: 国产美女精品一区二区| 国产高潮流白浆视频| 毛片网站观看| 国产毛片久久国产| 免费a在线观看播放| 国产永久免费视频m3u8| 亚洲日韩精品无码专区| 美女毛片在线| 成人国产免费| 欧美激情伊人| 日韩欧美国产中文| 久久久久青草线综合超碰| 99视频国产精品| 国产一级毛片高清完整视频版| 不卡视频国产| 亚州AV秘 一区二区三区| 亚洲精品大秀视频| 国产欧美日韩专区发布| 99re热精品视频中文字幕不卡| 久久精品国产精品一区二区| 日韩国产亚洲一区二区在线观看| 国产精品第| 中国一级特黄视频| 干中文字幕| 这里只有精品在线播放| 九九视频在线免费观看| 中文字幕免费播放| 综合社区亚洲熟妇p| AV老司机AV天堂| 欧美第一页在线| 区国产精品搜索视频| 粗大猛烈进出高潮视频无码| 国产一级视频久久| 欧美另类精品一区二区三区| 国产香蕉国产精品偷在线观看| 中文字幕66页| 亚洲第一视频网| h网站在线播放| 国产肉感大码AV无码| 黄色网址手机国内免费在线观看| 丁香五月激情图片| 中文字幕 欧美日韩| 亚洲日韩精品无码专区| 亚洲第一综合天堂另类专| 国产屁屁影院| 国产人成乱码视频免费观看| 国产亚洲精品资源在线26u| 在线a网站| 亚洲精品无码av中文字幕| 露脸一二三区国语对白| 国产精品网址你懂的| 亚洲综合极品香蕉久久网| 亚洲成人77777| 色吊丝av中文字幕| 国产成熟女人性满足视频| 青青青视频免费一区二区| 高清色本在线www| 少妇极品熟妇人妻专区视频| 国产免费精彩视频| 自慰高潮喷白浆在线观看| 久久精品这里只有国产中文精品| 国产女人18水真多毛片18精品| 制服丝袜国产精品| 色播五月婷婷| www.91中文字幕| 波多野结衣视频网站| 亚洲精品国产日韩无码AV永久免费网 | 狠狠做深爱婷婷综合一区| 日韩最新中文字幕| 久久99这里精品8国产| 国产91在线|日本| 欧美综合成人| 91精品专区国产盗摄| 美女一级毛片无遮挡内谢| 国产成人高清精品免费5388| 伊人久久久久久久| 免费毛片全部不收费的| 国产综合精品一区二区| 色综合中文| 人人妻人人澡人人爽欧美一区| 免费看美女自慰的网站| 亚洲愉拍一区二区精品|