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

    滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響

    2018-08-02 01:32:16張文毅
    體育科學(xué) 2018年7期

    胡 齊,陳 騏,張文毅

    ?

    滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響

    胡 齊1,陳 騏1,張文毅2

    1.國家體育總局體育科學(xué)研究所, 北京 100061; 2.北京體育大學(xué), 北京 100084

    目的:分析研究滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響。方法:建立精細(xì)化三維幾何模型與網(wǎng)格模型,采用部分時(shí)均(partially averaged Navier-Stokes,PANS)湍流模型對氣動(dòng)特性進(jìn)行計(jì)算流體力學(xué)(computational fluid dynamics,CFD)數(shù)值模擬,獲取不同滑雪板夾角下(24°、28°、32°、36°和40°)力和力矩、壓力分布以及流場形態(tài)。結(jié)果:1)總升力與總阻力、升力系數(shù)與阻力系數(shù)、總力矩與滑雪板力矩均隨滑雪板夾角增大先減小后增大最后小幅減小,但總升阻比基本保持不變,總體呈略微下降趨勢;2)滑雪板升力隨滑雪板夾角增大呈現(xiàn)單調(diào)遞增的趨勢,而滑雪板阻力隨滑雪板夾角增大先減小后增大;3)滑雪板正面壓力分布規(guī)律基本相同且比較均勻,而滑雪板背面下半部分壓力隨滑雪板夾角變化比較明顯;4)運(yùn)動(dòng)員身后以及滑雪板身后的流場形態(tài)主要以回流渦形式出現(xiàn),滑雪板背面回流渦結(jié)構(gòu)明顯比運(yùn)動(dòng)員身后復(fù)雜得多,而且回流渦強(qiáng)度也明顯大得多。結(jié)論:1)滑雪板夾角變化對總升阻比影響很小,但對總升力、總阻力、升力系數(shù)以及阻力系數(shù)產(chǎn)生一定影響,同時(shí)對運(yùn)動(dòng)員氣動(dòng)特性也產(chǎn)生一定影響。2)滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中起著更為重要作用,不僅要關(guān)注運(yùn)動(dòng)員身體姿態(tài),而且應(yīng)更加關(guān)注滑雪板姿態(tài);3)總力矩大小和滑雪板力矩大小隨滑雪板夾角變化趨勢一致,而且滑雪板力矩在總力矩中始終占據(jù)主導(dǎo)地位,選取合適的滑雪板夾角對飛行階段穩(wěn)定性至關(guān)重要,建議優(yōu)選的滑雪板夾角范圍為24°~32°。

    氣動(dòng)特性;跳臺滑雪飛行階段;滑雪板夾角;計(jì)算流體力學(xué);V型

    1 前言

    整個(gè)跳臺滑雪過程通常依次分為4個(gè)不同的階段:助滑、起跳、飛行和著陸,空氣動(dòng)力學(xué)在這4個(gè)階段中均起著重要的作用。在起跳后,運(yùn)動(dòng)員應(yīng)在飛行階段早期盡快達(dá)到一個(gè)穩(wěn)定的飛行位置,并確保受到小阻力[13],同時(shí)應(yīng)完成向后旋轉(zhuǎn)和向前旋轉(zhuǎn)角動(dòng)量的完全平衡[16]。跳臺滑雪運(yùn)動(dòng)員在0.5s內(nèi)完成飛行姿態(tài)的穩(wěn)定[21],并且在這個(gè)短時(shí)間內(nèi)(約為15 m的飛行距離),很有可能運(yùn)動(dòng)員不能進(jìn)行任何有意的姿態(tài)調(diào)整,可以說起跳階段決定隨后的飛行姿態(tài)。同樣,這也意味著起跳時(shí)所發(fā)生的錯(cuò)誤不能在飛行階段得到有效糾正,但成功起跳的優(yōu)勢可能會(huì)因飛行過程中錯(cuò)誤而喪失。一般來說,飛行階段被認(rèn)為是整個(gè)跳臺滑雪過程中氣動(dòng)特性體現(xiàn)最為明顯的階段。諸多對跳臺滑雪飛行階段氣動(dòng)特性研究文獻(xiàn)[2,5,6,8-10,12,17,18,20]采用了風(fēng)洞試驗(yàn)測量或計(jì)算機(jī)模擬與現(xiàn)場測量相結(jié)合等方法。普遍認(rèn)為,滑雪器材(滑雪板和滑雪服)在跳臺滑雪飛行階段的作用不可忽略。與以往的滑雪板平行飛行方式相比,V型飛行方式具有更好的氣動(dòng)性能,使運(yùn)動(dòng)員能夠更多地向前傾斜[19,20]。近年來,為了使氣流經(jīng)過滑雪板更加平緩,在飛行過程中滑雪板姿態(tài)越來越受到人們的關(guān)注。Virmavirta等[22]通過風(fēng)洞試驗(yàn)研究單個(gè)滑雪板空氣動(dòng)力特性,將滑雪板懸掛在六分量風(fēng)洞平衡測量裝置上,這種試驗(yàn)研究方法或多或少偏理想化,顯然并不是所有的結(jié)果都適用于真實(shí)飛行階段。

    上述這些方法均不能直觀地顯示或分析運(yùn)動(dòng)員周圍的流場。能夠在運(yùn)動(dòng)過程中分析氣動(dòng)力、顯示壓力分布和詳細(xì)流場信息的另一種技術(shù)是計(jì)算流體力學(xué)(CFD),該技術(shù)最近被應(yīng)用于自行車和游泳等運(yùn)動(dòng)中。CFD不需要物理裝置,成本低,可以模擬大多數(shù)物理現(xiàn)象,只需一次模擬就可以收集大量的數(shù)據(jù)。盡管CFD已經(jīng)被認(rèn)為是未來跳臺滑雪空氣動(dòng)力學(xué)研究的重要工具,但迄今為止發(fā)表的CFD相關(guān)研究較少[1,3,4,7,14,15],尤其是針對飛行階段就更少[1,4,7,14]。Keizo等[3]人采用CFD技術(shù)研究了跳雪運(yùn)動(dòng)員在起跳過程中姿態(tài)對氣動(dòng)特性的影響,著重比較兩種不同跳躍者的起跳情況,結(jié)論認(rèn)為,起跳過程中氣動(dòng)力特性在較短的時(shí)間內(nèi)發(fā)生了動(dòng)態(tài)變化,兩種起跳方式的氣動(dòng)力特性差異較大。Gardan等[1]人采用CFD方法來探討迎風(fēng)角和速度對氣動(dòng)力的影響,數(shù)值結(jié)果表明,在早期飛行階段速度對升力系數(shù)和阻力系數(shù)的影響很小,相反,迎風(fēng)角變化對作用于運(yùn)動(dòng)員身體上的升力和阻力有很大的影響。

    綜上所述,國內(nèi)、外報(bào)道姿態(tài)與跳臺滑雪飛行階段氣動(dòng)特性的CFD研究主要集中于運(yùn)動(dòng)員身體姿態(tài),鮮有報(bào)道滑雪板姿態(tài)的CFD研究[7,14],即使有也僅是涉及而已。本研究擬建立運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的精細(xì)化三維幾何模型與網(wǎng)格模型,采用PANS湍流模型進(jìn)行CFD數(shù)值模擬,獲取不同滑雪板夾角下力和力矩、壓力分布以及流場形態(tài),分析研究V型飛行方式滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響。

    2 研究對象與方法

    2.1 研究對象

    研究對象為跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)。根據(jù)Müller等[9]人統(tǒng)計(jì)分析結(jié)果,跳臺滑雪運(yùn)動(dòng)員的身體形態(tài)特征選取平均值,即身高為177 cm,身體質(zhì)量指數(shù)(BMI)為19.5,軀干(坐高)與身高比值0.532,滑雪板長度258 cm,寬11.5 cm。

    2.2 研究方法

    2.2.1 控制方程

    本研究采用大渦模擬技術(shù),該技術(shù)對預(yù)測鈍體周圍的流動(dòng)分離有優(yōu)勢。為了獲得更精確的結(jié)果,采用了PANS湍流模型,其控制方程表達(dá)式如下:

    式中i為已分解流場速度,為時(shí)間,為流體密度,為介質(zhì)粘性系數(shù),為湍流粘性系數(shù),f為未分解湍動(dòng)能比率,為未分解湍動(dòng)能耗散率比率,k為未分解局部時(shí)均化湍動(dòng)能,為未分解局部時(shí)均化湍動(dòng)能耗散率。

    其中:

    采用有限體積法對控制方程進(jìn)行離散計(jì)算,壓力和速度的耦合采用SIMPLEC算法,時(shí)間采用二階差分格式進(jìn)行離散,湍動(dòng)能以及速度項(xiàng)采用二階迎風(fēng)格式,時(shí)間步長為0.000 1 s。

    2.2.2 模型建立與網(wǎng)格無關(guān)性驗(yàn)證

    根據(jù)研究對象特征,建立跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的三維實(shí)體模型,對運(yùn)動(dòng)員的身體特征進(jìn)行精細(xì)化建模,手指、耳朵、眼罩、面部特征可以從圖1中清晰分辨。

    圖1 運(yùn)動(dòng)員/滑雪板整體系統(tǒng)精細(xì)化實(shí)體模型圖

    Figure 1. The Refined Solid Model of the Whole Athlete/Skis System

    注:左圖為正視圖,右圖為側(cè)視圖。

    跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)計(jì)算區(qū)域的尺寸大小為長13.5 m、寬5 m、高6 m,如圖2所示。采用商用前處理軟件ICEM對計(jì)算區(qū)域進(jìn)行四面體網(wǎng)格離散,在計(jì)算區(qū)域中心對稱面上體積網(wǎng)格分布和運(yùn)動(dòng)員表面網(wǎng)格分布如圖3所示。運(yùn)動(dòng)員模型的表面網(wǎng)格分辨率最小達(dá)到0.5 mm。對于體網(wǎng)格,在運(yùn)動(dòng)員周圍區(qū)域配置了較細(xì)的網(wǎng)格單元,以高精度捕獲尾流結(jié)構(gòu)。另一方面,為了降低計(jì)算量,較粗的網(wǎng)格單元被分配到遠(yuǎn)離運(yùn)動(dòng)員的區(qū)域。

    圖2 計(jì)算區(qū)域示意圖

    Figure. 2. Computational Domain

    圖3 網(wǎng)格分布示意圖

    Figure 3. Mesh Distribution

    注:左圖為體網(wǎng)格,右圖為運(yùn)動(dòng)員表面網(wǎng)格。

    為了滿足PANS模型計(jì)算要求,運(yùn)動(dòng)員周邊采用網(wǎng)格加密,并保證y plus控制在100之內(nèi),此控制要求在先前的研究中獲得了證實(shí)[3]。針對上述網(wǎng)格模型,通過均勻加密選定了4種網(wǎng)格密度,網(wǎng)格點(diǎn)數(shù)從800萬~2 357萬,進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證,驗(yàn)證結(jié)果如表1所示。4種網(wǎng)格驗(yàn)證計(jì)算獲得的升阻比均為1.8。由此可見,即使選定800萬網(wǎng)格節(jié)點(diǎn)的計(jì)算域離散方案,也能夠準(zhǔn)確預(yù)測跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的氣動(dòng)力學(xué)性能。

    表1 網(wǎng)格無關(guān)性驗(yàn)證結(jié)果

    2.2.3 邊界條件與計(jì)算工況

    邊界條件的設(shè)置如下:

    1)進(jìn)口設(shè)置為速度進(jìn)口,進(jìn)口速度根據(jù)飛行速度設(shè)定;2)出口設(shè)置為壓力出口,壓力為大氣壓101 325 Pa;3)中間截面為周期性條件;4)其他壁面為無滑移邊界條件;5)氣體為不可壓縮空氣;6)在常重力下,重力加速度設(shè)置為g0=9.807 m/s2。

    跳臺滑雪運(yùn)動(dòng)員在飛行階段的姿態(tài)參數(shù)包括迎風(fēng)角、滑雪板與速度方向夾角α、滑雪板與身體夾角θ、上半身彎曲角度β、滑雪板夾角λ,如圖4所示。Müller等[10,17]人的研究成果給出這些姿態(tài)參數(shù)通常的取值范圍。同時(shí),圖5顯示了在飛行階段整個(gè)系統(tǒng)受力情況。阻力系數(shù)(C)和升力系數(shù)(C)的計(jì)算方法見公式(11)和(12),為模型的正面投影區(qū)域面積。

    在本研究中,開展了5種不同滑雪板夾角下跳臺滑雪飛行階段動(dòng)特性和流動(dòng)結(jié)構(gòu)特性CFD研究。其中,滑雪板夾角λ分別取值為24°、28°、32°、36°和40°,同時(shí)迎風(fēng)角φ=30°、滑雪板與速度方向夾角α=30°、滑雪板與身體夾角θ=16°、上半身彎曲角度β=18°,速度V=29 m/s。據(jù)此,為每個(gè)滑雪板夾角λ的運(yùn)動(dòng)員/滑雪板整體系統(tǒng)生成計(jì)算網(wǎng)格,分別進(jìn)行了CFD數(shù)值模擬。

    圖4 飛行階段姿態(tài)參數(shù)示意圖

    Figure 4. Attitude Parameters in Flight

    圖5 飛行階段整體系統(tǒng)受力情況圖

    Figure. 5. All the Forces Acting on the Athlete/Skis System during Flight Phase

    3 結(jié)果

    3.1 力和力矩

    作用于運(yùn)動(dòng)員/滑雪板整體系統(tǒng)上的氣動(dòng)力可分為升力和阻力,而且這些力的作用點(diǎn)絕大多數(shù)不會(huì)在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)質(zhì)心上,因此極有可能會(huì)產(chǎn)生相應(yīng)的力矩。表2列出了不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)以及滑雪板的力學(xué)特性結(jié)果。圖6~10顯示了運(yùn)動(dòng)員/滑雪板整體系統(tǒng)以及滑雪板的力學(xué)特性隨滑雪板夾角變化曲線。圖6和圖7中各項(xiàng)力均為作用在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)或滑雪板上同一性質(zhì)的合力。圖8中阻力系數(shù)與升力系數(shù)分別由公式(11)和(12)計(jì)算得到。圖9中力矩是相對于運(yùn)動(dòng)員/滑雪板整體系統(tǒng)質(zhì)心的力矩,“+”代表力矩作用后使整體系統(tǒng)逆時(shí)針旋轉(zhuǎn),“-”代表力矩作用后使整體系統(tǒng)順時(shí)針旋轉(zhuǎn)。圖10中的升阻比由升力除以阻力計(jì)算得到的。

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總升力以及滑雪板升力隨滑雪板夾角變化曲線如圖6所示??梢钥闯觯偵﹄S滑雪板夾角增大先減小后增大最后小幅減小,即當(dāng)夾角不大于28°時(shí),總升力隨夾角增大而減?。划?dāng)夾角大于28°且小于36°時(shí),總升力隨夾角增大而增大;當(dāng)夾角大于36°時(shí),總升力隨夾角增大小幅減小。滑雪板升力隨滑雪板夾角增大呈現(xiàn)單調(diào)遞增的趨勢,即當(dāng)夾角不大于28°時(shí),升力隨夾角增大而小幅增加;當(dāng)夾角大于28°且小于36°時(shí),升力隨夾角增大而大幅增大;當(dāng)夾角大于36°時(shí),升力隨夾角增大小幅增加。總升力的最大增幅是滑雪板升力的兩倍左右。

    表2 力學(xué)特性結(jié)果

    圖6 升力變化曲線圖

    Figure 6. Lift Force Changes with Ski Opening Angle

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總阻力以及滑雪板阻力隨滑雪板夾角變化曲線如圖7所示??梢钥闯?,總阻力隨滑雪板夾角增大變化趨勢與總升力一致,即先減小后增大最后小幅減小。滑雪板阻力隨滑雪板夾角增大先減小后增大,即當(dāng)夾角不大于28°時(shí),阻力隨夾角增大而減?。划?dāng)夾角大于28°且小于36°時(shí),阻力隨夾角增大而大幅增大;當(dāng)滑雪板夾角大于36°時(shí),阻力隨夾角增大小幅增加??傋枇Φ淖畲笤龇腔┌遄枇Φ膬杀蹲笥摇?/p>

    圖7 阻力變化曲線圖

    Figure 7. Drag Force Changes with Ski Opening Angle

    不同雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)阻力系數(shù)與升力系數(shù)隨滑雪板夾角變化曲線如圖8所示。不難看出,阻力系數(shù)隨滑雪板夾角增大變化趨勢與升力系數(shù)一致,即先減小后增大最后小幅減小,同樣與總升力以及總阻力隨滑雪板夾角增大變化趨勢一致。

    圖8 阻力系數(shù)、升力系數(shù)變化曲線圖

    Figure 8. Changes of Drag Coefficient and Lift Coefficient with Ski Opening Angle

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總力矩以及滑雪板力矩隨滑雪板夾角變化曲線如圖9所示??梢钥闯?,總力矩與滑雪板力矩均是產(chǎn)生逆時(shí)針旋轉(zhuǎn)的效果,其總力矩大小和滑雪板力矩大小隨滑雪板夾角增大變化趨勢一致,即先減小后增大最后小幅減小,同樣與總升力以及總阻力隨滑雪板夾角增大變化趨勢一致。

    圖9 力矩變化曲線圖

    Figure 9. Torque Changes with Ski Opening Angle

    圖10 升阻比變化曲線圖

    Figure 10. Lift-to-drag Ratio Changes with Ski Opening Angle

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總升阻比以及滑雪板升阻比隨滑雪板夾角變化曲線如圖10所示??梢钥闯?,總升阻比隨滑雪板夾角增大變化基本保持不變,總體呈略微下降趨勢。另一方面,滑雪板升阻比隨滑雪板夾角增大先大幅增大(在滑雪板夾角不大于28°之前),之后略微下降。

    3.2 壓力分布

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的正面壓力分布如圖11所示。對于運(yùn)動(dòng)員來說,正面壓力在不同夾角下分布規(guī)律基本一致,高壓區(qū)主要集中于頭部、上半身、胯關(guān)節(jié)、膝關(guān)節(jié)等部位,低壓區(qū)主要集中于身體側(cè)部。同時(shí)不難看出,夾角為28°和32°時(shí)正面高壓區(qū)面積偏小,尤其是滑雪板夾角為28°時(shí)正面高壓區(qū)面積明顯最小。對于滑雪板來說,正面壓力在不同夾角下分布規(guī)律基本相同,而且分布比較均勻。

    圖11 正面壓力分布圖

    Figure 11. Frontal Pressure Distribution

    不同滑雪板夾角λ下滑雪板背面壓力分布如圖12所示??梢钥闯觯哓?fù)壓區(qū)主要集中在滑雪板頂端和中部。進(jìn)一步可看出,滑雪板背面上半部分的壓力隨滑雪板夾角變化比較細(xì)微,然而滑雪板背面下半部分壓力變化比較明顯。

    圖12 滑雪板背面壓力分布圖

    Figure 12. Pressures Distribution of the Ski Back

    3.3 流場形態(tài)

    不同滑雪板夾角λ下渦流形態(tài)和氣流速度流線如圖13所示,左邊為背面正視圖,右邊為側(cè)視圖。運(yùn)動(dòng)員身后以及滑雪板后面的渦流形態(tài)均主要以回流渦結(jié)構(gòu)形式出現(xiàn),此結(jié)構(gòu)會(huì)加速能量耗散。對于運(yùn)動(dòng)員來說,其身后主要產(chǎn)生2對渦結(jié)構(gòu),分別在肩與上臂結(jié)合處以及腰與胯結(jié)合處的后方,當(dāng)夾角為28°時(shí),身后的回流渦結(jié)構(gòu)強(qiáng)度最低。對于滑雪板來說,無論角度是否改變,其后側(cè)有6對渦結(jié)構(gòu)產(chǎn)生,當(dāng)夾角為28°時(shí),后側(cè)的回流渦結(jié)構(gòu)強(qiáng)度也相對最低,之后隨著夾角增大,其后側(cè)的回流渦結(jié)果變得更加無序,尤其是滑雪板下半部分的后側(cè)。不難看出,當(dāng)夾角為24°時(shí),運(yùn)動(dòng)員腰與胯結(jié)合處后方的回流渦結(jié)構(gòu)強(qiáng)度較小,渦流尺度較大的區(qū)域主要集中在滑雪板后側(cè);當(dāng)夾角大于28°時(shí),運(yùn)動(dòng)員身后的回流渦強(qiáng)度高于滑雪板后側(cè),在腰與胯結(jié)合處后方形成尺度較大的回流渦。同時(shí),觀察到在負(fù)壓力梯度作用下,氣流速度在運(yùn)動(dòng)員以及滑雪板上均發(fā)生了流動(dòng)分離,尤其是在滑雪板下半部分。當(dāng)夾角為24°時(shí)流動(dòng)分離程度以及尾跡尺寸明顯最大,當(dāng)夾角為28°時(shí)流動(dòng)分離程度以及尾跡尺寸明顯最小,當(dāng)夾角大于28°時(shí)流動(dòng)分離程度以及尾跡尺寸逐漸變大。

    4 討論

    本研究中,所獲得的氣動(dòng)特性部分結(jié)果與Gardan等[1]人的部分研究結(jié)果相似。通過總結(jié)風(fēng)洞試驗(yàn)測量結(jié)果以及數(shù)值計(jì)算結(jié)果,他們公布了迎風(fēng)角為29°、滑雪板與身體夾角θ為17.4°、上半身彎曲角度β為16.6°、滑雪板夾角λ為37.9°的姿態(tài)條件下總升阻比為1.89。在本研究中,在迎風(fēng)角為30°、滑雪板與身體夾角θ為16°、上半身彎曲角度β為18°、滑雪板夾角λ為40°的姿態(tài)條件下總升阻比為1.818。兩者的結(jié)果非常接近,驗(yàn)證了本研究CFD計(jì)算結(jié)果的正確性。

    從圖6~8可以看出,總升力與總阻力、阻力系數(shù)與升力系數(shù)隨滑雪板夾角增大變化趨勢一致。同時(shí),圖10顯示出總升阻比隨滑雪板夾角變化基本保持不變,總體呈略微下降趨勢,也進(jìn)一步說明總升力和總阻力變化趨勢的一致性。因此,可以認(rèn)為,滑雪板夾角變化對運(yùn)動(dòng)員/滑雪板整體系統(tǒng)升阻比影響很小,但對阻力系數(shù)與升力系數(shù)有一定影響。

    圖13 流場形態(tài)示意圖

    Figure 13. Flow Field Form

    注:左圖為后視圖,右圖為側(cè)視圖。

    N?rstrud等[14]通過數(shù)值計(jì)算研究認(rèn)為,V型飛行方式代表了更為先進(jìn)的空氣動(dòng)力,滑雪板在產(chǎn)生良好的升力方面起著重要作用。從表2、圖6和圖7可以看出,不僅滑雪板升力一直占據(jù)總升力63%以上,而且滑雪板阻力也一直占據(jù)總阻力65%以上,可以認(rèn)為,滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中起著更為重要作用,滑雪板氣動(dòng)特性的改善將會(huì)獲得更為有利的效果,在分析飛行階段姿態(tài)氣動(dòng)特性時(shí),不僅要關(guān)注運(yùn)動(dòng)員身體姿態(tài),而且應(yīng)更加關(guān)注滑雪板姿態(tài)。

    另外,圖6展示了滑雪板升力隨滑雪板夾角增大呈現(xiàn)單調(diào)遞增的趨勢,同時(shí),圖8顯示了在滑雪板夾角大于28°后升力系數(shù)總體呈增大趨勢,即增大夾角能夠提高總氣動(dòng)升力系數(shù)。同時(shí),總力矩大小和滑雪板力矩大小隨滑雪板夾角變化趨勢一致,而且滑雪板力矩值是總力矩值的1.17倍以上,如圖9和表2所示。當(dāng)夾角大于28°時(shí),隨著夾角增大,總力矩值和滑雪板力矩值越大,產(chǎn)生向前傾的效果越明顯,可能會(huì)對運(yùn)動(dòng)員的穩(wěn)定飛行越不利。本研究結(jié)果與Virmavirta[22]的風(fēng)洞試驗(yàn)研究結(jié)論一致。他們重點(diǎn)研究了單個(gè)滑雪板的空氣動(dòng)力特性,結(jié)論認(rèn)為,大的滑雪板夾角改善了氣動(dòng)升力系數(shù),但是,張開如此寬度的滑雪板夾角的飛行姿態(tài)可能會(huì)使運(yùn)動(dòng)員飛行變得不穩(wěn)定。因此,運(yùn)動(dòng)員應(yīng)根據(jù)自身技術(shù)特點(diǎn)以及訓(xùn)練效果來選擇滑雪板夾角,切忌一味為了改善氣動(dòng)升力系數(shù)而選擇大的滑雪板夾角??偠灾?,滑雪板力矩在總力矩中始終占據(jù)主導(dǎo)地位,選取合適的滑雪板夾角對飛行穩(wěn)定性至關(guān)重要,優(yōu)選的滑雪板夾角范圍為24°~32°。

    隨著滑雪板夾角變化,滑雪板背面下半部分壓力變化比較明顯,同時(shí)滑雪板正面壓力基本不變且滑雪板背面上半部分的壓力變化也比較細(xì)微(圖11、12),這也是導(dǎo)致滑雪板升力、阻力以及力矩變化的主要原因。當(dāng)夾角為28°時(shí),滑雪板背面下半部分的次高負(fù)壓區(qū)域明顯偏少,低負(fù)壓區(qū)域明顯偏多,尤其在滑雪板底端部位,故此時(shí)滑雪板正面背面壓差偏小,從而滑雪板升力、阻力以及向前傾的力矩值偏小,尤其是力矩值,這也從表2以及圖6~9中不難看出。

    運(yùn)動(dòng)員身后以及滑雪板后側(cè)的流場形態(tài)主要以回流渦形式出現(xiàn),回流渦結(jié)構(gòu)會(huì)使能量耗散速度加快。如圖13所示,滑雪板后側(cè)回流渦結(jié)構(gòu)明顯比運(yùn)動(dòng)員身后復(fù)雜得多,而且綜合尺度及強(qiáng)度也明顯大得多,說明滑雪板背面能量損失更大,導(dǎo)致滑雪板的升力和阻力數(shù)值較大,并且在總升力和總阻力中占據(jù)較大的比重,從而進(jìn)一步驗(yàn)證了滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中的確起著更為重要作用。同時(shí),對于運(yùn)動(dòng)員來說,當(dāng)夾角為24°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度明顯偏大,即此時(shí)運(yùn)動(dòng)員的升力與阻力偏大;當(dāng)夾角為28°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度明顯最小,即此時(shí)運(yùn)動(dòng)員的升力與阻力最小,當(dāng)夾角大于28°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度逐漸變大,直至夾角為40°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度又略微變小一些,即運(yùn)動(dòng)員的升力與阻力先變大后略微變小。從表2以及圖6~9中不難看出,運(yùn)動(dòng)員的升力和阻力隨滑雪板夾角變化趨勢也與上述變化趨勢一致,這兩種結(jié)果之間也起到良好的相互驗(yàn)證效果。這些研究結(jié)果表明,滑雪板夾角變化對運(yùn)動(dòng)員氣動(dòng)特性產(chǎn)生一定影響,運(yùn)動(dòng)員的升力以及阻力的變化趨勢與總升力以及總阻力變化趨勢一致。

    5 結(jié)論

    1.滑雪板夾角變化對運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的升阻比影響很小,但對總升力、總阻力、升力系數(shù)以及阻力系數(shù)均產(chǎn)生一定影響,同時(shí)對運(yùn)動(dòng)員本身的氣動(dòng)特性也產(chǎn)生一定影響??偵涂傋枇﹄S滑雪板夾角的變化趨勢與運(yùn)動(dòng)員的升力和阻力變化趨勢一致。

    2.滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中起著更為重要作用,滑雪板氣動(dòng)特性的改善將會(huì)獲得更為有利的效果,在分析飛行階段姿態(tài)氣動(dòng)特性時(shí),不僅要關(guān)注運(yùn)動(dòng)員身體姿態(tài),更加要關(guān)注滑雪板姿態(tài)。

    3.總力矩大小和滑雪板力矩大小隨滑雪板夾角變化趨勢一致,而且滑雪板力矩在總力矩中始終占據(jù)主導(dǎo)地位,選取合適的滑雪板夾角對飛行階段穩(wěn)定性至關(guān)重要,運(yùn)動(dòng)員應(yīng)根據(jù)自身技術(shù)特點(diǎn)以及訓(xùn)練效果來選擇滑雪板夾角,切記一味為了改善氣動(dòng)升力系數(shù)而選擇大的滑雪板夾角,建議優(yōu)選的滑雪板夾角范圍為24°~32°。

    [1] GARDAN N, SCHNEIDER A, POLIDORI G,. Numerical investigation of the early flight phase in ski-jumping[J]. J Biomech, 2017, 59(1): 29-34.

    [2] JUNG A, STAAT M, MULLER W. Flight style optimization in ski jumping on normal, large, and ski flying hills[J]. J Biomech, 2014, 47(3): 716-722.

    [3] KEIZO Y, MAKOTO T, JUN I,. Effect of posture on the aerodynamic characteristics during take-off in ski jumping [J]. J Biomech, 2016, 49(15): 3688-3696.

    [4] LEE K D, PARK M J, KIM K Y. Optimization of ski jumper’s posture considering lift-to-drag ratio and stability[J]. J Biomech, 2012, 45(12): 2125–2132.

    [5] MARQUES-BRUNA P, GRIMSHAW P. Mechanics of flight in ski jumping: aero-dynamic stability in pitch[J]. Sports Technol, 2009, 2(1–2): 24–31.

    [6] MARQUES-BRUNA P, GRIMSHAW P. Mechanics of flight in ski jumping: aero-dynamic stability in roll and yaw[J]. Sports Technol, 2009, 2(3-4): 111–120.

    [7] MEILE W, REISENBERGER E, MAYER M,. Aerodynam-ics of ski jumping: experiments and CFD simulations[J]. Exp Fluids, 2006, 41(6): 949–964.

    [8] MULLER W.Performance factors in ski jumping[J]. J Biomech, 2006, 39(1): 192–213.

    [9] MULLER W, GROSCHL W, MULLER R,. Underweight in ski jumping: The solution of the problem[J]. Int J Sports Med, 2006, 27(11): 926-934.

    [10] MULLER W, PLATZER D, SCHMOLZER B. Dynamics of human flight on skis: improvements in safety and fairness in ski jumping[J]. J Biomech, 1996, 29(8): 1061–1068 .

    [11] MULLER W, PLATZER D, SCHMOLZER B. Scientific approach to ski safety[J]. Nature, 1995, 375(6531): 455-466.

    [12] MURAKAMI M, IWASE M, SEO K,.High-speed video image analysis of ski jumping flight posture[J]. Sports Eng, 2014, 17(4): 217–225.

    [13] MURAKAMI M, IWASE M, SEO K,. Ski jumping flight skill analysis based on high-speed video image[J]. Procedia Eng, 2010, 2(2): 2381–2386.

    [14] N?RSTRUD H, ?YE I J. On CFD simulation of ski jumping[M]//Computational Fluid Dynamics for Sport Simulation, 2009, 63–82.

    [15] RYU M, CHO L, CHO J. Aerodynamic analysis on postures of ski jumpers during flight using computational fluid dynamics[J]. TJp Soc Aeronaut S, 2015, 58(4): 204–212.

    [16] SCHWAMEDER H. Biomechanics research in ski jumping: 1991–2006[J]. Sport Biomech, 2008, 7(1): 114–136.

    [17] SCHMOLZER B, MULLER W. Individual flight styles in ski jumping: results obtained during Olympic games competitions[J]. J Biomech, 2005, 38(5): 1055–1065.

    [18] SCHMOLZER B, MULLER W. The importance of being light: aerodynamic forces and weight in ski jumping[J]. J Biomech, 2002, 35(8): 1059–1069.

    [19] SEO K, MURAKAMI M, YOSHIDA K. Optimal flight technique for V-style ski jumping[J]. Sports Eng, 2004, 7(2): 97–104.

    [20] SEO K, WATANABE I, MURAKAMI M. Aerodynamic force data for a V-style ski jumping flight[J]. Sports Eng, 2004, 7(1): 31–39.

    [21] VIRMAVIRTA M, ISOLEHTO J, KOMI P,. Characteristics of the early flight phase in the Olympic ski jumping competition[J]. J Biomech, 2005, 38(11): 2157–2163.

    [22] VIRMAVIRTA M, KIVEKAS J. Aerodynamics of isolated ski jumping ski[C]. XXV ISB Congress, Glasgow, 2015: 215-223.

    Effect of the Ski Opening Angle on the Aerodynamic Characteristics during Flight in Ski-jumping

    HU Qi1, CHEN Qi1, ZHANG Wen-yi2

    1. China Institute of Sport Science, Beijing 100061, China; 2. Beijing Sport University, Beijing 100084, China.

    Objective: This study aimed to investigate the effects of the ski opening angle on aerodynamic characteristics during the flight in ski-jumping. Methods: Refined 3D geometric models and mesh models were established, and aerodynamic characteristics were numerically simulated by using computational fluid dynamics (CFD) with partially averaged Navier-Stokes (PANS) turbulence model. The force and moment, pressure distribution and flow field form of different ski opening angles (24°, 28°, 32°, 36° and 40°) were obtained. Results: 1) As ski opening angle increases ,the total lift force and that drag force, the coefficients of lift and drag, the total torque and the skis’ torque all decrease firstly, then increase, and finally decrease slightly, but the total lift-to-drag ratio almost remains unchanged, and shows a slight downward trend as a whole. 2) The skis’ lift force increases monotonously with the increase of ski opening angle, while the skis’ drag force decreases firstly and then increases. 3) The front pressure distribution of skis is approximately the same and fairly uniform, but the back pressures distribution of the lower part of skis varies obviously with ski opening angle. 4) The flow field form behind athlete and skis mainly appears in the form of reflux vortex. The vortex structure behind skis is much more complicated than that behind the athlete, and the vortex intensity is also much greater. Conclusion: 1) The ski opening angle has little effect on the total lift-to-drag ratio, but has some influence on the total lift and total drag and the coefficients of lift and drag, as well as on the aerodynamic characteristics of the athlete. 2) The aerodynamic characteristics of skis play a more important role in the whole athlete/skis system, so more attention should be paid to skis’ attitude while paying attention to athletes' body posture. 3) The variation trend of total torque and skis’ torque with ski opening angle is consistent, and skis’ torque has always been dominant in the total torque, so it is very important to select the proper ski opening angle for the stability of flight phase, and it is suggested that the optimum range of ski opening angle is 24°~32°.

    1000-677X(2018)07-0042-08

    10.16469/j.css.201807017

    G804.6

    A

    2018-06-09;

    2018-07-12

    國家體育總局體育科學(xué)研究所基本科研業(yè)務(wù)費(fèi)資助項(xiàng)目(18-29)、中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助課題(2018RC016)。

    胡齊,男,助理研究員,碩士,主要研究方向?yàn)轶w育工程、流體與機(jī)械工程、運(yùn)動(dòng)生物力學(xué),E-mail:huqi@ciss.cn; 陳騏,男,副研究員,學(xué)士,主要研究方向?yàn)轶w育工程,E-mail:chenqi@ciss.cn; 張文毅,女,助理研究員,碩士,主要研究方向?yàn)檫\(yùn)動(dòng)生物力學(xué)、高等教育管理, E-mail: zhangwy1228@126.com。

    久久精品亚洲熟妇少妇任你| av欧美777| 免费不卡黄色视频| 欧美激情高清一区二区三区| 亚洲片人在线观看| 999精品在线视频| 欧美久久黑人一区二区| 男人的好看免费观看在线视频 | www.999成人在线观看| 精品国产超薄肉色丝袜足j| svipshipincom国产片| 视频在线观看一区二区三区| 久久婷婷人人爽人人干人人爱 | 黄色女人牲交| 久热爱精品视频在线9| 成人三级做爰电影| 满18在线观看网站| 免费av毛片视频| 中文字幕精品免费在线观看视频| 免费在线观看黄色视频的| 亚洲伊人色综图| 色婷婷久久久亚洲欧美| 久久午夜综合久久蜜桃| 两个人视频免费观看高清| www.自偷自拍.com| 久久青草综合色| 免费在线观看黄色视频的| 成人亚洲精品一区在线观看| 国产成人欧美| 黄色丝袜av网址大全| 一级,二级,三级黄色视频| 妹子高潮喷水视频| 法律面前人人平等表现在哪些方面| 一夜夜www| 日本 av在线| 亚洲色图av天堂| 欧美一级a爱片免费观看看 | 日韩一卡2卡3卡4卡2021年| 亚洲黑人精品在线| 日本精品一区二区三区蜜桃| 91精品国产国语对白视频| 久久九九热精品免费| 一级片免费观看大全| 久久精品国产亚洲av高清一级| 88av欧美| 天天一区二区日本电影三级 | 亚洲一区中文字幕在线| 欧美成人免费av一区二区三区| 女性生殖器流出的白浆| 欧美日韩黄片免| 狂野欧美激情性xxxx| 99在线人妻在线中文字幕| 在线观看午夜福利视频| 老司机深夜福利视频在线观看| 18禁裸乳无遮挡免费网站照片 | 亚洲精品粉嫩美女一区| 欧美精品啪啪一区二区三区| 日韩有码中文字幕| 国产精品久久久人人做人人爽| 欧美日韩精品网址| 日本 欧美在线| 日韩欧美国产一区二区入口| 国产成人啪精品午夜网站| 岛国视频午夜一区免费看| 男男h啪啪无遮挡| 亚洲欧美激情综合另类| 国产不卡一卡二| 久久久久久久久免费视频了| 午夜亚洲福利在线播放| 成人av一区二区三区在线看| 国产精品永久免费网站| 成人手机av| 国产一区在线观看成人免费| 色老头精品视频在线观看| 如日韩欧美国产精品一区二区三区| 老汉色av国产亚洲站长工具| 久久精品91蜜桃| 免费观看精品视频网站| 97超级碰碰碰精品色视频在线观看| 国产精品二区激情视频| 亚洲成人免费电影在线观看| 久久午夜综合久久蜜桃| 两个人视频免费观看高清| 亚洲av成人不卡在线观看播放网| 精品一区二区三区av网在线观看| 色老头精品视频在线观看| 亚洲伊人色综图| 国产亚洲欧美在线一区二区| 禁无遮挡网站| 美女扒开内裤让男人捅视频| netflix在线观看网站| 激情视频va一区二区三区| 精品国产一区二区三区四区第35| 婷婷六月久久综合丁香| 久久人妻熟女aⅴ| 国产97色在线日韩免费| 国产在线精品亚洲第一网站| 在线观看舔阴道视频| 悠悠久久av| av中文乱码字幕在线| 18禁国产床啪视频网站| 极品教师在线免费播放| 最近最新中文字幕大全免费视频| 一区二区三区精品91| 99久久精品国产亚洲精品| 最新美女视频免费是黄的| 一区二区三区国产精品乱码| 又黄又爽又免费观看的视频| 侵犯人妻中文字幕一二三四区| 美女大奶头视频| 午夜福利,免费看| 国产伦人伦偷精品视频| 亚洲国产精品成人综合色| 国产亚洲欧美98| 亚洲欧洲精品一区二区精品久久久| 国产精品一区二区免费欧美| 免费观看精品视频网站| 亚洲精品一卡2卡三卡4卡5卡| 97碰自拍视频| 亚洲第一青青草原| 琪琪午夜伦伦电影理论片6080| 中文字幕另类日韩欧美亚洲嫩草| 禁无遮挡网站| 亚洲国产精品久久男人天堂| 看免费av毛片| 黄色毛片三级朝国网站| av电影中文网址| 国产精品99久久99久久久不卡| 好男人电影高清在线观看| 丁香欧美五月| 日本精品一区二区三区蜜桃| 51午夜福利影视在线观看| 成年人黄色毛片网站| 搡老妇女老女人老熟妇| 亚洲aⅴ乱码一区二区在线播放 | 色老头精品视频在线观看| 亚洲人成网站在线播放欧美日韩| 久久久久久人人人人人| 两个人免费观看高清视频| 俄罗斯特黄特色一大片| 波多野结衣一区麻豆| 女人精品久久久久毛片| 成人18禁高潮啪啪吃奶动态图| 欧美乱码精品一区二区三区| 国产精品亚洲一级av第二区| 日韩中文字幕欧美一区二区| 亚洲五月色婷婷综合| 亚洲精品中文字幕在线视频| www.自偷自拍.com| 日本 av在线| 丝袜美腿诱惑在线| 别揉我奶头~嗯~啊~动态视频| 久久草成人影院| 午夜久久久久精精品| 一本大道久久a久久精品| 亚洲欧美日韩另类电影网站| 日韩免费av在线播放| 99re在线观看精品视频| 亚洲精品美女久久久久99蜜臀| 91精品三级在线观看| 免费看美女性在线毛片视频| 午夜免费激情av| 一卡2卡三卡四卡精品乱码亚洲| 女人被狂操c到高潮| 久久久久精品国产欧美久久久| 久久热在线av| 美国免费a级毛片| 亚洲欧美激情在线| 国产精品影院久久| 国产免费av片在线观看野外av| 在线天堂中文资源库| 大型av网站在线播放| 99在线视频只有这里精品首页| 日韩一卡2卡3卡4卡2021年| 两个人看的免费小视频| 国产成人精品久久二区二区免费| 又黄又粗又硬又大视频| 不卡一级毛片| 啦啦啦免费观看视频1| 久久精品国产清高在天天线| 1024香蕉在线观看| 亚洲熟妇熟女久久| 国产精品爽爽va在线观看网站 | 久久精品91蜜桃| 亚洲av成人不卡在线观看播放网| 涩涩av久久男人的天堂| 久久狼人影院| or卡值多少钱| 两性夫妻黄色片| 桃色一区二区三区在线观看| 久久午夜综合久久蜜桃| 热re99久久国产66热| 成年版毛片免费区| 自线自在国产av| 国产91精品成人一区二区三区| 久久中文字幕人妻熟女| 精品国产乱码久久久久久男人| 一本久久中文字幕| 午夜免费激情av| 十分钟在线观看高清视频www| 国产精品一区二区在线不卡| 精品卡一卡二卡四卡免费| 欧美一区二区精品小视频在线| 国产野战对白在线观看| 亚洲精品一区av在线观看| 亚洲激情在线av| 精品久久久久久久毛片微露脸| 亚洲成av人片免费观看| 欧美成狂野欧美在线观看| 美女 人体艺术 gogo| 老司机午夜福利在线观看视频| 国产精品一区二区免费欧美| 好看av亚洲va欧美ⅴa在| 午夜激情av网站| 91九色精品人成在线观看| 老司机深夜福利视频在线观看| 国产一区二区三区视频了| 色在线成人网| 黄色视频不卡| 视频区欧美日本亚洲| 大型av网站在线播放| 很黄的视频免费| svipshipincom国产片| 一区二区三区国产精品乱码| 侵犯人妻中文字幕一二三四区| 亚洲av成人不卡在线观看播放网| 国产精品久久久人人做人人爽| 精品久久久久久久毛片微露脸| 国产亚洲av嫩草精品影院| 久久草成人影院| 人妻丰满熟妇av一区二区三区| 一级黄色大片毛片| 欧美亚洲日本最大视频资源| av视频免费观看在线观看| 亚洲av日韩精品久久久久久密| 久久香蕉激情| 亚洲男人的天堂狠狠| 一本综合久久免费| 国产高清视频在线播放一区| 免费在线观看完整版高清| 午夜福利在线观看吧| 欧美人与性动交α欧美精品济南到| 波多野结衣一区麻豆| 黑人欧美特级aaaaaa片| 搡老熟女国产l中国老女人| 在线观看66精品国产| 男人舔女人的私密视频| 日本五十路高清| 色综合亚洲欧美另类图片| 黑丝袜美女国产一区| 国产精品 国内视频| 中文字幕人妻熟女乱码| 精品国产乱子伦一区二区三区| www.自偷自拍.com| 免费人成视频x8x8入口观看| 国产熟女午夜一区二区三区| 亚洲专区字幕在线| 悠悠久久av| 欧美日本中文国产一区发布| 成年女人毛片免费观看观看9| av免费在线观看网站| 在线永久观看黄色视频| 18禁黄网站禁片午夜丰满| 久久久久国产精品人妻aⅴ院| 国产欧美日韩一区二区三| 久久九九热精品免费| 丝袜美足系列| av视频免费观看在线观看| 精品欧美国产一区二区三| 中文字幕精品免费在线观看视频| 免费在线观看完整版高清| 日本一区二区免费在线视频| 搡老熟女国产l中国老女人| 又黄又爽又免费观看的视频| 不卡av一区二区三区| 在线观看www视频免费| av中文乱码字幕在线| 国产蜜桃级精品一区二区三区| 成在线人永久免费视频| 夜夜夜夜夜久久久久| 最近最新中文字幕大全电影3 | 国产一区二区激情短视频| 亚洲欧美激情综合另类| 老熟妇乱子伦视频在线观看| 欧美老熟妇乱子伦牲交| 可以在线观看的亚洲视频| 成人永久免费在线观看视频| 精品一品国产午夜福利视频| 99国产综合亚洲精品| 久久这里只有精品19| 国产成人免费无遮挡视频| 国产成人欧美| 国产欧美日韩精品亚洲av| 免费在线观看亚洲国产| 国产亚洲av高清不卡| 91成年电影在线观看| 国产av又大| 久久国产精品人妻蜜桃| 国产视频一区二区在线看| 欧美老熟妇乱子伦牲交| 久99久视频精品免费| 免费在线观看视频国产中文字幕亚洲| 免费av毛片视频| 操出白浆在线播放| 一区二区三区国产精品乱码| 97人妻精品一区二区三区麻豆 | 国产又爽黄色视频| 男女之事视频高清在线观看| 国产成人一区二区三区免费视频网站| 看免费av毛片| 久久久久久人人人人人| 国产伦人伦偷精品视频| 国产高清videossex| 色老头精品视频在线观看| 久久国产精品男人的天堂亚洲| 大陆偷拍与自拍| 精品无人区乱码1区二区| 亚洲精品av麻豆狂野| 韩国精品一区二区三区| 亚洲熟妇熟女久久| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆久久精品国产亚洲av| 午夜亚洲福利在线播放| av天堂久久9| 欧美+亚洲+日韩+国产| 制服丝袜大香蕉在线| 国产成人欧美在线观看| 曰老女人黄片| 变态另类丝袜制服| 十分钟在线观看高清视频www| 九色亚洲精品在线播放| 欧美午夜高清在线| 亚洲一卡2卡3卡4卡5卡精品中文| 一区二区三区国产精品乱码| 999久久久国产精品视频| 亚洲精品美女久久av网站| 麻豆av在线久日| 亚洲国产精品久久男人天堂| 欧美精品啪啪一区二区三区| 日本 av在线| 欧美人与性动交α欧美精品济南到| 国产亚洲精品综合一区在线观看 | 桃红色精品国产亚洲av| 高清在线国产一区| 亚洲欧美激情在线| 黄片小视频在线播放| 亚洲人成电影观看| 成人18禁高潮啪啪吃奶动态图| 在线观看免费视频日本深夜| 久久婷婷人人爽人人干人人爱 | 欧美日本视频| 婷婷精品国产亚洲av在线| 国产黄a三级三级三级人| 欧美人与性动交α欧美精品济南到| 91老司机精品| 纯流量卡能插随身wifi吗| 久久久水蜜桃国产精品网| 99在线视频只有这里精品首页| 成人18禁高潮啪啪吃奶动态图| 国产精品亚洲一级av第二区| 国产精品久久久久久精品电影 | 叶爱在线成人免费视频播放| 日本一区二区免费在线视频| 一区二区三区国产精品乱码| 午夜视频精品福利| 国产精品av久久久久免费| 成人亚洲精品av一区二区| 在线十欧美十亚洲十日本专区| 亚洲成人免费电影在线观看| 岛国视频午夜一区免费看| 国产亚洲av嫩草精品影院| 一级毛片女人18水好多| 亚洲国产欧美日韩在线播放| 99在线人妻在线中文字幕| 波多野结衣一区麻豆| 国产激情欧美一区二区| 12—13女人毛片做爰片一| videosex国产| 国产97色在线日韩免费| netflix在线观看网站| www.自偷自拍.com| 国产一区二区三区在线臀色熟女| 久久草成人影院| 老汉色∧v一级毛片| 亚洲精品国产色婷婷电影| 日韩国内少妇激情av| 亚洲精品粉嫩美女一区| 欧美大码av| 国产激情久久老熟女| 午夜激情av网站| 国产成人精品久久二区二区91| 久久人人爽av亚洲精品天堂| 丰满的人妻完整版| 成年女人毛片免费观看观看9| 老汉色∧v一级毛片| 日韩大码丰满熟妇| 亚洲色图av天堂| 欧美日韩亚洲综合一区二区三区_| 久久影院123| 午夜影院日韩av| bbb黄色大片| 亚洲精品中文字幕在线视频| 嫁个100分男人电影在线观看| 黄色成人免费大全| 国产精品 国内视频| 又黄又爽又免费观看的视频| 香蕉丝袜av| 精品卡一卡二卡四卡免费| 黄色视频不卡| 久久青草综合色| 亚洲第一青青草原| 午夜福利成人在线免费观看| 最好的美女福利视频网| 久久这里只有精品19| 色哟哟哟哟哟哟| 欧美成人一区二区免费高清观看 | 久久久水蜜桃国产精品网| 老司机午夜十八禁免费视频| 曰老女人黄片| 亚洲人成77777在线视频| 久久伊人香网站| 男人舔女人的私密视频| 黄频高清免费视频| 黑人欧美特级aaaaaa片| 涩涩av久久男人的天堂| 好男人在线观看高清免费视频 | 中文字幕av电影在线播放| 日本欧美视频一区| 这个男人来自地球电影免费观看| 咕卡用的链子| 岛国视频午夜一区免费看| 亚洲avbb在线观看| 国产视频一区二区在线看| 国产亚洲精品久久久久5区| 成人18禁在线播放| 亚洲伊人色综图| 黄片大片在线免费观看| 精品卡一卡二卡四卡免费| 少妇被粗大的猛进出69影院| 在线观看免费午夜福利视频| 日韩国内少妇激情av| 午夜福利,免费看| 久久精品国产99精品国产亚洲性色 | 成人免费观看视频高清| av视频在线观看入口| 国产成人av激情在线播放| 一a级毛片在线观看| 中文字幕人妻熟女乱码| 少妇熟女aⅴ在线视频| 成人亚洲精品一区在线观看| av天堂久久9| av中文乱码字幕在线| 大香蕉久久成人网| 两性午夜刺激爽爽歪歪视频在线观看 | 人成视频在线观看免费观看| 看免费av毛片| 日韩 欧美 亚洲 中文字幕| 亚洲精品中文字幕一二三四区| 黄色a级毛片大全视频| 他把我摸到了高潮在线观看| 91麻豆av在线| 亚洲第一av免费看| 亚洲av电影不卡..在线观看| 国产亚洲精品第一综合不卡| 波多野结衣一区麻豆| 国产一级毛片七仙女欲春2 | 在线观看免费午夜福利视频| 亚洲三区欧美一区| 国产成人精品无人区| 熟妇人妻久久中文字幕3abv| 国产精品日韩av在线免费观看 | 在线天堂中文资源库| 最新美女视频免费是黄的| 午夜福利一区二区在线看| 性欧美人与动物交配| 日韩精品青青久久久久久| 久久久久国内视频| 欧美日本视频| 在线免费观看的www视频| 人妻久久中文字幕网| ponron亚洲| 亚洲国产精品成人综合色| 黄色毛片三级朝国网站| 日本 av在线| 国产成人一区二区三区免费视频网站| 亚洲欧美精品综合一区二区三区| 大陆偷拍与自拍| 免费观看精品视频网站| 亚洲第一青青草原| 国产日韩一区二区三区精品不卡| 一区福利在线观看| 两人在一起打扑克的视频| 亚洲电影在线观看av| 免费看十八禁软件| 曰老女人黄片| 国产亚洲精品久久久久5区| 热re99久久国产66热| 美女高潮喷水抽搐中文字幕| 午夜福利成人在线免费观看| 乱人伦中国视频| 香蕉国产在线看| 曰老女人黄片| 天天添夜夜摸| 国内精品久久久久精免费| 日韩精品免费视频一区二区三区| 国产精品秋霞免费鲁丝片| 日韩中文字幕欧美一区二区| 亚洲欧美一区二区三区黑人| 亚洲精品av麻豆狂野| 亚洲激情在线av| 久久久久国产一级毛片高清牌| 纯流量卡能插随身wifi吗| 一级,二级,三级黄色视频| 亚洲专区字幕在线| 女生性感内裤真人,穿戴方法视频| 女人被躁到高潮嗷嗷叫费观| 亚洲黑人精品在线| 香蕉久久夜色| 琪琪午夜伦伦电影理论片6080| 久久人妻av系列| 在线观看一区二区三区| 老司机午夜十八禁免费视频| 免费在线观看影片大全网站| 99riav亚洲国产免费| 亚洲国产精品合色在线| 露出奶头的视频| 免费在线观看日本一区| 18禁国产床啪视频网站| 老司机福利观看| 人人妻人人澡欧美一区二区 | 免费在线观看日本一区| 亚洲激情在线av| 亚洲欧美激情综合另类| 国产一级毛片七仙女欲春2 | 69av精品久久久久久| 欧美乱码精品一区二区三区| 19禁男女啪啪无遮挡网站| 最新美女视频免费是黄的| 国产不卡一卡二| 在线av久久热| 自线自在国产av| 亚洲色图 男人天堂 中文字幕| 97人妻天天添夜夜摸| 久久久久久久久久久久大奶| 久久精品亚洲熟妇少妇任你| 亚洲性夜色夜夜综合| 国产成人免费无遮挡视频| 亚洲成av片中文字幕在线观看| 国产av一区在线观看免费| 一进一出好大好爽视频| 国产激情久久老熟女| 国产xxxxx性猛交| 性欧美人与动物交配| 神马国产精品三级电影在线观看 | 国产精品1区2区在线观看.| av视频在线观看入口| 国产免费男女视频| 国产精品自产拍在线观看55亚洲| 亚洲精品av麻豆狂野| 天堂√8在线中文| 精品一区二区三区av网在线观看| 国产欧美日韩一区二区三| 欧美中文综合在线视频| 99re在线观看精品视频| 精品无人区乱码1区二区| 日本免费一区二区三区高清不卡 | 欧美激情久久久久久爽电影 | 国产亚洲精品综合一区在线观看 | 50天的宝宝边吃奶边哭怎么回事| 国产91精品成人一区二区三区| www.精华液| 男女下面进入的视频免费午夜 | 亚洲熟妇熟女久久| 高潮久久久久久久久久久不卡| 99热只有精品国产| 欧美大码av| 欧美激情 高清一区二区三区| 亚洲欧洲精品一区二区精品久久久| 日韩欧美一区二区三区在线观看| 成人三级做爰电影| 波多野结衣av一区二区av| 在线av久久热| 亚洲欧美日韩无卡精品| 人人妻人人澡人人看| 国产欧美日韩一区二区三区在线| 香蕉国产在线看| 女人被躁到高潮嗷嗷叫费观| 欧美日韩福利视频一区二区| 伊人久久大香线蕉亚洲五| 三级毛片av免费| 国产精品影院久久| e午夜精品久久久久久久| 成人三级黄色视频| 中文字幕色久视频| 99re在线观看精品视频| 国产成人啪精品午夜网站| 国产高清激情床上av| 欧美国产日韩亚洲一区| 免费在线观看完整版高清| 亚洲欧美日韩高清在线视频| 久久香蕉国产精品| 国产亚洲精品av在线| 少妇裸体淫交视频免费看高清 | 日韩成人在线观看一区二区三区| 欧美丝袜亚洲另类 | 啪啪无遮挡十八禁网站| 国产精品98久久久久久宅男小说| 香蕉国产在线看| 国产欧美日韩一区二区三区在线| 亚洲成a人片在线一区二区| 啦啦啦免费观看视频1| 国产av一区在线观看免费| 99国产精品免费福利视频| 在线国产一区二区在线| 欧美成人性av电影在线观看|