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

間隔織物雙劍桿織機平行打緯機構箱體的振動分析

2024-07-20 00:00:00楊建成劉艷哲黃子趙全鵬劉健
天津工業大學學報 2024年3期

摘" " 要: 針對雙劍桿織機共軛凸輪與四連桿組合的平行打緯機構箱體隨著織機轉速的提高產生強烈振動,影響打緯精度和可靠性的問題,采用理論分析與虛擬仿真結合的方法對打緯機構及箱體進行動應力和振動特性研究。建立了平行打緯機構的動態靜力分析方程,運用ADAMS軟件對機構進行多體動力學仿真,求解箱體所受的支反力;構建箱體—彈簧阻尼—柔性基礎的系統模型,分析箱體在豎直方向的受迫振動,并在仿真過程中探索減小箱體振動的方法。研究結果表明:擺臂軸孔處箱體所受支反力峰值最大,打緯推程階段箱體的受力最為集中;當凸輪軸轉速為260 r/min時齒輪嚙合頻率與箱體一階固有頻率相等,會產生共振,應避免使用這一轉速工作;箱體在豎直方向上產生周期性振動,增加螺栓個數和附加動力減震器有助于改善箱體的振動特性,振幅分別減少38%和23%。

關鍵詞: 雙劍桿織機; 平行打緯機構; 支反力; 動態靜力分析; 彈簧阻尼; 減振

中圖分類號: TS103.135" " " " " "文獻標志碼: A" " " " " " " " 文章編號: 1671-024X(2024)03-0082-07

Vibration analysis of parallel beating up mechanism box of

double rapier loom for spacer fabric

YANG Jiancheng1,2,3, LIU Yanzhe1,2,3, HUANG Ziwen1,2,3, ZHAO Quanpeng1,2,3, LIU Jian1,2,3

(1. School of Mechanical Engineering, Tiangong University, Tianjin 300387, China; 2. Tianjin Key Laboratory of Advanced Mechatronics Equipment Technology, Tiangong University, Tianjin 300387, China; 3. Tianjin Virtual Simulation Experi-mental Teaching Center of Mechanical Fundamentals and Textile Equipment Design, Tiangong University, Tianjin 300387, China)

Abstract: Aiming at the problem that the parallel beating up mechanism box of double rapier loom with conjugate cam and four-bar linkage generates strong vibration with the increase of loom speed, which destroys the accuracy and reliability of beating up, the dynamic stress and vibration characteristics of the beating up mechanism and box are studied by combining theoretical analysis with virtual simulation. The dynamic static analysis equation of parallel beating up mechanism is established, and the multi-body dynamics simulation of the mechanism is carried out using ADAMS software to solve the reaction force on the box. The system model of box-spring damping-flexible foundation is built, the forced vibration of box in vertical direction is analyzed, and the method to reduce the vibration of box is explored in the simulation process. The results show that the peak value of the reaction force on the box body at the shaft hole of the swing arm is the largest, and the force on the box body is the most concentrated during the beating up and pushing process. When the cam shaft rotates at 260 r/min, the gear meshing frequency is equal to the first order natural frequency of the box, which will produce resonance. The speed should be avoided during operation. The box produces periodic vibration in the vertical direction. Increasing the number of bolts and adding dynamic shock absorbers will help to improve the vibration characteristics of the box, the amplitude decrease by 38% and 23%, respectively.

Key words: double rapier loom; parallel beating up mechanism; support reaction force; dynamic static analysis; spring damping; damping

共軛凸輪與四連桿組合的平行打緯機構主要應用于高速大間隔織物雙劍桿織機,這一打緯機構有效解決了織造間隔織物時要求每層緯紗所受打緯力應均勻的問題。目前裝備共軛凸輪平行打緯機構的國產雙劍桿織機轉速普遍不高,墻板、箱體的振動是影響打緯機構精度和壽命的重要因素之一。針對平行打緯機構進行動應力和振動分析對促進國內間隔織物雙劍桿織機的高速化進程具有重要意義。目前關于打緯機構振動特性方面的研究主要集中在四連桿打緯機構分析,金玉珍等[1]對四連桿打緯機構進行數值分析,確定了噴氣織機墻板的最大振動和最大應力位置。田建學等[2]通過仿真實驗總結出墻板所受沖擊力和變形隨曲軸轉速提高的變化規律。關于打緯機構動力學方面的研究主要集中在連桿打緯機構和共軛凸輪擺動打緯機構分析,沈丹峰等[3]針對墻板上曲軸孔與搖軸孔承受的支撐力,建立了動力學關系式,并對兩孔的安裝角度進行了優化設計。魏展等[4]應用凱恩方程建立了共軛凸輪擺動打緯機構的動力學模型。針對共軛凸輪平行打緯機構的動力學和振動研究尚處于起步階段。

本研究對平行打緯機構及箱體進行動態靜力分析、多體動力學仿真、模態分析和振動特性分析,并在仿真的過程中探究了減小箱體振動的改進方法。本研究旨在為共軛凸輪平行打緯機構的設計和優化提供參考。

1 機構原理及動應力分析

本研究以國外某先進雙劍桿天鵝絨織機中的共軛凸輪平行打緯機構為分析對象,其結構原理如圖1所示。

圖1中,打緯過程分為3個階段:推程階段、回程階段和靜止階段。共軛凸輪繞回轉軸做順時針勻速轉動,其主副凸輪分別對2個滾子施加正壓力,推桿與擺臂固結并繞回轉軸轉動,從而帶動鋼筘的往復擺動。靜止階段時滾子不受凸輪正壓力,鋼筘在打緯后死心位置靜止。

打緯機構通過回轉軸孔施加給箱體的作用力是箱體振動的主要激勵。對打緯機構進行動態靜力分析,假設凸輪按理想規律做等速回轉運動,不計入摩擦阻力、打緯阻力和凸輪滾子之間的碰撞,將推桿與滾子視為一個剛體。簡化后的打緯機構動態靜力分析模型如圖2所示。

分別對從動連桿機構和凸輪機構進行動力學分析。連桿部分各構件的受力如圖3所示。

圖3中:Fab、Mab分別為構件a對構件b的作用力和力矩;Fba = -Fab。Rca為鉸鏈c到構件a質心的矢量;Sa、Ga、φa分別為構件a的質心、重力和相對于整體坐標的角位移。

擺臂:

F32 x + F12 x = m2F32" y + F12 y = m2 + G2F32y R22x - F32x R22y + F12y R12 x - F12x R12y + M62 =J22(1)

連桿:

F43x - F32x = m3F43y - F32y = m3S3 + G3F43y R33x - F43x R33y - F32y R23x + F32x R23y = J33" (2)

搖桿:

F14x - F43 x = m4S4F14y - F43 y = m4S4+ G4F14y R44x - F14x R44y - F43y R34x + F43x R34y = J44" (3)

對凸輪機構進行動力學分析,各構件的受力如圖4所示。

凸輪:

F15x - FRx = 0F15y - FRy =" G5Md + FRy R5py- FRx R5py = J5 5" " "(4)

式中:FR為凸輪作用在滾子上的正壓力;Md為打緯機構的平衡力矩;p為凸輪與滾子的動態接觸點。推程階段主凸輪與滾子之間的接觸點受力,回程階段副凸輪與滾子之間的接觸點受力。

推桿:

F16x Rx = 0F16y - FRy =" G6FRy R1px - FRx R1px - M62 = J6 6" "(5)

聯立以上各式,即得到整個平行打緯機構的動態靜力分析方程,其是由15個線性方程組成的方程組,共有16個未知量,為了求解,可將平衡力矩Md用實際機器測量的凸輪軸輸入轉矩代入近似計算,以凸輪軸轉速300 r/min為例,測量凸輪軸驅動力矩一個周期內的數值記為打緯機構的平衡轉矩,這樣方程組未知量變為15個,可以進一步求解。將線性方程組寫成矩陣的形式:

AR = B(6)

式中:A為15 × 15的已知矩陣,其中元素與構件的質心位置有關;B為15維的已知列陣,其中包含了構件所受慣性力、慣性力矩和重力;R為15維的未知列陣,其中包含了各構件所受作用力和力矩.

R = A-1B (7)

矩陣A、B中元素大都是隨時間變化的,求解一次方程組只能得到某一位置的解,故還需對機構進行周期離散化,可以運用Matlab編程對m個離散位置分別求解。箱體所受支反力為箱體在回轉軸孔處對其他構件作用力的反力。

2 多體動力學仿真

通過對機構進行多體動力學仿真得到機構的動態響應和動應力。首先對機構進行參數化設計,使用作圖法,根據整機方案設計和工藝要求確定從動連桿機構的尺寸和擺臂動程,得到擺臂、連桿、搖桿的長度分別為0.29、0.11、0.2 m,擺臂角動程為33.49°。共軛凸輪從動件選用改進梯形運動規律,結合共軛凸輪機構的設計公式,運用Matlab編程得到共軛凸輪廓線數據和凸輪從動系統的幾何尺寸。在SolidWorks中建立機構的虛擬樣機模型,按照工程要求對機構進行結構設計,可提高仿真實驗的可靠性。將三維模型文件保存為x_t文件導入ADAMS軟件中,得到打緯機構的動力學模型如圖5所示。

對模型進行材料設置,其中擺臂、連桿、搖桿、推桿的材料均選用球墨鑄鐵,凸輪選用42GrMo,滾子選用NUP圓柱滾子向心滾動軸承,材料為軸承鋼。

添加約束和驅動后即可對系統進行動力學仿真和數據后處理。以凸輪軸設計轉速300 r/min進行仿真,筘座線性位移和速度響應分別如圖6和圖7所示。

凸輪推程角回程角均為90°,則0.05 s處為打緯前死心位置,0.1 s處為后死心位置。從圖6、圖7可以看出,鋼筘在前死心位置時x方向線性位移最大,x、y方向線性速度均為0。此時鋼筘處于豎直狀態并短時間內做水平運動,即實現鋼筘的平行打緯運動。從而驗證了理論設計的正確性和動力學模型的有效性。測量圖5所示的凸輪軸孔、擺臂軸孔、搖桿軸孔處回轉副的作用力的反力即為打緯箱體所受的支反力如圖8—圖10所示。

從圖8—圖10可以看出,打緯箱體回轉軸孔處所受支反力變化呈現周期性,且與凸輪旋轉周期一致。其中擺臂軸的支反力峰值大于凸輪軸的支反力峰值,搖桿軸支反力峰值最小。凸輪軸的最大支反力出現在回程階段,擺臂軸和搖桿軸的最大支反力出現在推程階段。這說明箱體在擺臂軸孔處應力沖擊較大,在打緯推程階段應力最為集中。在靜止階段3個軸孔的支反力值較小且無波動,此時各軸孔只承受機構靜平衡的支反力。

3 箱體的振動特性分析

運用ANAYS軟件中Workbench模塊對箱體進行模態分析,將箱體模型的x_t文件導入軟件,箱體材料為HT250,材料密度為7 280 kg/m3,彈性模量為1.38×1011 N/m2,泊松比0.156,劃分網格,包含了375 430個節點、214 744個四面體單元。添加底面固定約束后進行模態求解,得到打緯箱體的前10階固有頻率和各階振型描述如表1所示。

共軛凸輪通過主軸與凸輪軸的一級齒輪減速帶動轉動,齒輪嚙合頻率為:

f = (8)

式中:n為凸輪軸轉速;z為從動齒輪齒數。

由機構的理論設計可得主軸與凸輪軸間減速齒輪的齒數為97/19,將表1中模態分析計算的箱體前三階固有頻率代入式(8),當凸輪軸轉速為260 r/min、343 r/min、430 r/min時齒輪嚙合頻率分別與箱體一階、二階、三階固有頻率相等,引起箱體共振。目前國內雙劍桿織機凸輪轉速普遍在200 r/min以下,本課題來源項目的理想凸輪平行打緯目標轉速為300 r/min,對于本課題所涉織機,應避免使用凸輪260 r/min這一轉速工作。在后續的提速改進工程中,在不變更箱體結構時,還應避免凸輪343 r/min這一轉速,以減小箱體振動對織機產生的危害。

應用ANSYS中的2次開發語言APDL模塊求解柔性箱體受周期性支反力作用下的受迫振動[5-20],將第2節動力學仿真解得的箱體在凸輪軸、擺臂軸、搖桿軸所受周期性支反力離散成一系列沖擊載荷,逐個求解其動響應,在不同節點處拾取位移響應,可得最大變形出現在箱體下方節點處,最大變形為0.012 8 mm,振動變形較小,原因是此時離共振點較遠,沒有發生共振,在此不做深入研究。

忽略箱體的內部變形,將螺栓用彈簧阻尼單元描述[7],為了提高仿真的精確性,同時計入基礎的柔性變形。整個箱體結構可以等效成一個打緯箱體—彈簧阻尼—柔性基礎的耦合系統,模型如圖11所示。

圖11中,箱體在豎直方向上的運動為一單自由度系統受周期性激振力作用下的受迫振動。

打緯箱體所受支反力在y方向的分力引起箱體在豎直方向上的振動,設支反力y方向分力的合力為周期為T的函數F(t),將函數F(t)展開成傅里葉級數:

F(t)= a0 + (ak cos (kωt) + bk sin (kωt)) (9)

式中:a0、ak、bk為傅里葉系數。若系統的質量、剛度和阻尼分別為M、K和C,則此時受迫振動的穩態響應可以寫為:

x =βkcos(kω - arctan(bk /ak) - θk)(10)

式中:βk為第k 階響應的動力放大系數;θk為第k 階響應與第k 階激振力間的相位差。

由式(10)可以看出,在周期激振力作用下,系統的振動響應由無窮階簡諧振動而合成。 將箱體-基礎模型導入ADAMS,基礎為20 mm厚的鋼板焊接而成的中空橫梁底座,利用軟件中的創建柔性體功能將基礎柔性化,并對原有剛體模型進行覆蓋替換。在箱體底部與基礎間建立彈簧阻尼單元用以表述螺栓連接,根據結合面的螺栓影響錐理論用下式計算螺栓的等效剛度和等效阻尼:

K0 = π(?漬 + 2dtan γ)2K/4C0 = π(?漬 + 2d tanγ)2C/4 " "(11)

式中:φ為螺栓帽的外切圓直徑;d為螺釘長度;γ為影響錐的半錐角度數,一般取為 20 °;K、C分別為結合面單位面積剛度和阻尼。

本文中箱體底面固定螺栓選用 M16x50規格的內六角圓柱頭螺栓,連接方式為螺釘連接,基礎底座材料為Q235。將第2節動力學仿真所求箱體所受支反力的y方向的時變激勵數據以樣條曲線的方式導入模型,利用軟件中的AKISPL函數擬合曲線點集,在相應的軸孔添加作用力,約束箱體自由度為上下移動。運行交互仿真,忽略運動開始時載荷沖擊對系統造成的影響,得到箱體在周期性支反力激振下的穩態受迫振動位移響應如圖12所示。

從圖12可見,箱體在周期性激振力的作用下產生周期性振動,在一個周期內,箱體的質心位移響應不只有一個波峰。正方向上的振動幅值大于負方向上的振動幅值,最大振幅出現在打緯回程階段,為0.026 mm,且為向上振動。靜止階段箱體較為平穩,無劇烈振動。

以減小箱體正方向振幅為優化目標,提出2種減振方法:第1種,由公式(10)可以看出,單自由度系統受周期性激振力時,增大系統剛度可以減小穩態振幅,增大系統剛度最直接的方法就是增加螺栓個數,對底面筋板結構進行優化并在其上增加相應螺栓孔,將底面固定螺栓由4個增加至8個。第2種是采用附加動力減振器的方法,即在箱體附加1個輔助質量,用彈性元件使之與箱體連接如圖13所示。圖13所示為動力減振器的原理圖。

當主系統m1振動時,輔助質量m2也隨之振動。利用輔助質量的動力作用,使其加在主系統上的動力與激振力互相抵消,從而使主系統的振動得到抑制。這個附加系統就稱為動力減振器。本文在箱體上附加一個質量為40 kg、彈性元件剛度為2 ×105 N/mm構成的減振器。在仿真軟件中對上述2種方法進行驗證,得到箱體減振后穩態正方向振動位移如圖14所示。

由圖14可以看出,增加螺栓個數和附加動力減振器的措施對箱體正方向的振幅都有一定的抑制作用,相較不進行減振的原系統,增加螺栓個數后箱體振幅減小38%,附加動力減振器后振幅減小23%,但在實際的工程應用中,箱體結構不易變更,增加螺栓個數的方法具有局限性,且螺栓個數過多會造成螺栓功能的冗余,并對安裝造成不便。相較而言,附加動力減振器的方法可根據實際的減振效果調換附加質量和彈性元件,具有靈活性較強的優勢。

4 結 論

針對共軛凸輪平行打緯機構箱體的動應力和振動特性進行理論分析和仿真實驗分析,得出以下結論:

(1) 擺臂軸孔處箱體所受支反力峰值最大,打緯推程階段箱體的受力最為集中。

(2) 當凸輪軸轉速為260 r/min時齒輪嚙合頻率與箱體一階固有頻率相等,會引起箱體共振,應避免使用這一轉速進行工作。研究結果對實際的工程項目具有一定的指導意義。

(3) 箱體在周期性支反力的作用在豎直方向上產生周期性振動,最大振幅發生在打緯回程階段,為0.026 mm。僅從減小箱體正方向振幅的角度出發,增加螺紋個數和附加動力減振器均有利于改善箱體的振動特性,增加螺栓個數后箱體振幅減小38%,附加動力減振器后振幅減小23%。研究結果對共軛凸輪平行打緯機構箱體的設計和優化提供了參考。

參考文獻:

[1]" " 金玉珍, 胡小冬, 林培峰, 等. 噴氣織機打緯機構及墻板的振動特性[J]. 紡織學報, 2016, 37(7): 131-136, 141.

JIN Y Z, HU X D, LIN P F, et al. Vibration characteristics of wallboard and four-bar linkage beating-up mechanism of air-jet loom[J]. Journal of Textile Research, 2016, 37(7): 131-136, 141 (in Chinese).

[2]" " 田建學, 張玨, 屈曉淵, 等. 基于計算機仿真技術的織機打緯機構設計[J]. 西北大學學報(自然科學版), 2015, 45(6): 942-946.

TIAN J X, ZHANG J, QU X Y, et al. Design of loom′s beating-up mechanism based on computer simulation technology[J]. Journal of Northwest University(Natural Science Edition), 2015, 45(6): 942-946(in Chinese).

[3]" " 沈丹峰, 趙寧, 張華安, 等. 織機四連桿打緯機構對主墻板的動應力分析[J]. 紡織學報, 2013, 34(1): 110-115.

SHEN D F, ZHAO N, ZHANG H A, et al. Dynamic stress analysis of four-bar linkage beating-up mechanism on loom wallboard[J]. Journal of Textile Research, 2013, 34(1): 110-115(in Chinese).

[4]" " 魏展, 金國光, 袁汝旺, 等. 高速共軛凸輪打緯機構柔性動力學分析[J]. 機械工程學報, 2017, 53(3): 81-89.

WEI Z, JIN G G, YUAN R W, et al. Flexible dynamic analysis of high-speed conjugate cam beating-up mechanism[J]. Journal of Mechanical Engineering, 2017, 53(3): 81-89 (in Chinese).

[5]" " FU S P, LUO S M, HUANG H L. Analysis of transmission box vibration characteristics under random road torsional excitation[J]. Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, 2022, 236(12): 2582-2597.

[6]" " WU H, WU P B, XU K, et al. Research on vibration charac-teristics and stress analysis of gearbox housing in high-speed trains[J]. IEEE Access, 2019, 7: 102508-102518.

[7]" " 伍良生, 馬淑慧, 屈重年, 等. 彈簧-阻尼動力學單元螺栓連接結合面研究[J]. 機械設計與制造, 2014(1): 4-6.

WU L S, MA S H, QU C N, et al. The study of bolted joint surface of the spring-damping dynamics unit[J]. Machinery Design amp; Manufacture, 2014(1): 4-6 (in Chinese).

[8]" " LI M M, LIU Q W, DAI G H, et al. Vibration transfer path analysis of double-layer box for marine reducer[J]. Journal of Vibration Engineering amp; Technologies, 2021, 9(6): 1077-1089.

[9]" " 胡玉飛,張建超,陳湛,等.內部激勵下高速動車齒輪箱振動響應評估[J]. 北京交通大學學報,2022,46(4):148-156.

HU Y F, ZHANG J C, CHEN Z,et al. Evaluation of vibration response of high-speed train gearboxes under internal excita-tion[J]. Journal of Beijing Jiaotong University, 2022, 46(4):148-156(in Chinese).

[10]" YU X L, HUANGFU Y F, YANG Y, et al. Gear fault diagno-sis using gear meshing stiffness identified by gearbox housing vibration signals[J]. Frontiers of Mechanical Engineering, 2022, 17(4): 57.

[11]" 黃曉丹, 裴幫, 關鶴, 等. 基于剛柔耦合的地鐵齒輪箱故障仿真研究[J]. 機械傳動, 2021, 45(10): 137-143, 156.

HUANG X D, PEI B, GUAN H, et al. Fault simulation resea-rch of metro gearbox based on rigid-flexible coupling[J]. Jour-nal of Mechanical Transmission, 2021, 45(10): 137-143, 156(in Chinese).

[12]" 劉輝, 張喜清, 項昌樂. 多源動態激勵下變速箱箱體結構的動態響應分析[J]. 兵工學報, 2011, 32(2): 129-135.

LIU H, ZHANG X Q, XIANG C L. Dynamic response analysis of gearbox housing under multi-source dynamic excitation[J]. Acta Armamentarii, 2011, 32(2): 129-135(in Chinese).

[13]" SONG L M. Experimental research on vibration and noise re-duction of vehicular vacuum circuit breaker′s equipment box[J]. Science Discovery, 2018, 6(2): 85.

[14]" SEON S O, KIM J O. Vibration and acoustic characteristics depending on housing structures of ultrasonic sensors[J]. Tra-nsactions of the Korean Society for Noise and Vibration Engi-neering, 2019, 29(2): 157-165.

[15]" 任亞峰, 常山, 劉更, 等. 齒輪-箱體-基礎耦合系統的振動分析[J]. 華南理工大學學報(自然科學版), 2017, 45(5): 38-44.

REN Y F, CHANG S, LIU G, et al. Vibration analysis of gear-housing-foundation coupled system[J]. Journal of South China University of Technology(Natural Science Edition), 2017, 45(5): 38-44(in Chinese).

[16]" VIDHYA B, SRINIVAS K N. Vibration analysis including sta-tor, rotor, housing and dynamic response analysis of Flux Re-versal Generator[J]. Journal of Electrical Systems and Informa-tion Technology, 2018, 5(2): 144-157.

[17]" 黃冠華, 王興宇, 梅桂明, 等. 內外激勵下高速列車齒輪箱箱體動態響應分析[J]. 機械工程學報, 2015, 51(12): 95-100.

HUANG G H, WANG X Y, MEI G M, et al. Dynamic respo-nse analysis of gearbox housing system subjected to internal and external excitation in high-speed train[J]. Journal of Me-chanical Engineering, 2015, 51(12): 95-100 (in Chinese).

[18]" 王強龍, 李耀彬, 傘曉剛, 等. 螺栓連接條件下結構振動模態的時頻分析方法[J]. 機械設計與制造, 2016(7): 72-76.

WANG Q L, LI Y B, SAN X G, et al. Model analysis of struc-ture with bolt connections using the method of transient respo-nse and frequency transformation[J]. Machinery Design amp; Ma-nufacture, 2016(7): 72-76(in Chinese).

[19]" LEE T H. Muscular activity and acceleration of box vibration in manual holding tasks: Effects of load and height of the load′s center of gravity[J]. International Journal of Occupational Safe-ty and Ergonomics, 2020, 26(1):15-19.

[20]" 周香琴, 楊勝皎, 王琴龍, 等. 劍桿織機慣性載荷變化規律及提速試驗[J]. 紡織學報, 2017, 38(8):144-149.

ZHOU X Q, YANG S J, WANG Q L, et al. Inertial load change rule and speed increasing test of rapier loom[J]. Journal of Textile Research, 2017, 38(8): 144-149(in Chinese).

本文引文格式:

楊建成,劉艷哲,黃子文,等. 間隔織物雙劍桿織機平行打緯機構箱體的振動分析[J]. 天津工業大學學報,2024, 43(3): 82-88.

YANG J C, LIU Y Z, HUANG Z W, et al. Vibration analysis of parallel beating up mechanism box of double rapier loom for spacer fabric[J]. Journal of Tiangong University, 2024, 43(3): 82-88(in Chinese).

主站蜘蛛池模板: 精品亚洲麻豆1区2区3区| 又爽又大又黄a级毛片在线视频| 欧美区日韩区| 色亚洲成人| 五月婷婷激情四射| 国产人在线成免费视频| 高清码无在线看| 国产精品嫩草影院av| 亚洲综合色婷婷| 国产精品所毛片视频| 69国产精品视频免费| 国内99精品激情视频精品| 成人欧美在线观看| 欧美日本在线观看| 国产人妖视频一区在线观看| 国产精品任我爽爆在线播放6080| 色综合久久综合网| 国内精品免费| 亚洲免费毛片| 无码久看视频| 在线观看欧美国产| av在线人妻熟妇| 久久精品国产一区二区小说| 婷婷午夜影院| 91偷拍一区| 四虎永久在线精品国产免费| 91尤物国产尤物福利在线| 国产在线观看99| 国产va在线观看| 婷婷色丁香综合激情| 激情五月婷婷综合网| 国产经典免费播放视频| 成人免费黄色小视频| 欧美亚洲日韩不卡在线在线观看| 亚洲av日韩av制服丝袜| 色婷婷成人| 欧美国产精品不卡在线观看 | 69av免费视频| 国产成熟女人性满足视频| 欧美黄色网站在线看| 国产精品主播| 亚洲国产成熟视频在线多多 | 中国美女**毛片录像在线| 日韩精品免费在线视频| 毛片久久网站小视频| 九色视频一区| 91在线国内在线播放老师| 日韩精品久久久久久久电影蜜臀| 91久久偷偷做嫩草影院精品| 91在线激情在线观看| 亚洲人成网站18禁动漫无码 | 久久国产精品77777| 91精品专区国产盗摄| 欧美视频在线观看第一页| 欧美日韩导航| 午夜小视频在线| 国产成人喷潮在线观看| 一级做a爰片久久免费| 欧美精品成人一区二区视频一| 欧美69视频在线| 亚洲中文字幕av无码区| 19国产精品麻豆免费观看| 久久semm亚洲国产| 亚洲欧美日韩中文字幕在线一区| 日韩在线播放欧美字幕| 亚洲中字无码AV电影在线观看| 国语少妇高潮| 全免费a级毛片免费看不卡| 免费A级毛片无码无遮挡| 波多野衣结在线精品二区| 亚洲aaa视频| 欧美日韩国产在线观看一区二区三区| www.国产福利| 青草视频网站在线观看| 国产成人综合久久| 国产精品成人啪精品视频| 成人午夜在线播放| 97久久精品人人| 国产一级毛片网站| 日日噜噜夜夜狠狠视频| 天天操精品| 国产精品美人久久久久久AV|