• <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| 一边摸一边抽搐一进一小说| 国产免费男女视频| 又爽又黄无遮挡网站| 亚洲精品影视一区二区三区av| www.色视频.com| 三级国产精品欧美在线观看| 天天一区二区日本电影三级| 久久久久久久久久久丰满| 男插女下体视频免费在线播放| 亚洲av熟女| 在线观看一区二区三区| 大香蕉久久网| 少妇熟女欧美另类| 男女边吃奶边做爰视频| 小说图片视频综合网站| 麻豆国产av国片精品| 一级黄色大片毛片| 国产91av在线免费观看| 一级毛片我不卡| av免费在线看不卡| 国产三级在线视频| 变态另类成人亚洲欧美熟女| 在线国产一区二区在线| 淫妇啪啪啪对白视频| 日韩亚洲欧美综合| 国产精品人妻久久久影院| 亚洲国产欧洲综合997久久,| 少妇猛男粗大的猛烈进出视频 | 亚洲性久久影院| 国语自产精品视频在线第100页| 欧美+亚洲+日韩+国产| 九九久久精品国产亚洲av麻豆| 国产精品精品国产色婷婷| 男插女下体视频免费在线播放| 99久国产av精品| av专区在线播放| 免费av不卡在线播放| 久久精品人妻少妇| 91久久精品国产一区二区三区| 色播亚洲综合网| 国产人妻一区二区三区在| 欧美成人精品欧美一级黄| 中文字幕久久专区| 欧美三级亚洲精品| 久久中文看片网| 综合色丁香网| 日本黄色片子视频| 国产伦在线观看视频一区| 91久久精品电影网| 长腿黑丝高跟| 老司机影院成人| 蜜桃久久精品国产亚洲av| 久久久欧美国产精品| 日本-黄色视频高清免费观看| 精品久久久久久久久亚洲| 色综合站精品国产| 久久久久国内视频| 色av中文字幕| 日韩,欧美,国产一区二区三区 | 精品久久久噜噜| 欧美激情久久久久久爽电影| 一级a爱片免费观看的视频| 激情 狠狠 欧美| 欧美色视频一区免费| 尤物成人国产欧美一区二区三区| 精品福利观看| 丰满的人妻完整版| 久久午夜福利片| 精品日产1卡2卡| 亚洲性久久影院| 波多野结衣高清无吗| 精品午夜福利在线看| 亚洲精品久久国产高清桃花| 在线播放国产精品三级| 成人亚洲欧美一区二区av| 国产精品福利在线免费观看| 成人漫画全彩无遮挡| 乱系列少妇在线播放| 亚洲精品乱码久久久v下载方式| 亚洲无线观看免费| 亚洲成人精品中文字幕电影| 亚洲18禁久久av| av视频在线观看入口| 男人的好看免费观看在线视频| av中文乱码字幕在线| 久久精品国产亚洲网站| 国内久久婷婷六月综合欲色啪| 干丝袜人妻中文字幕| 亚洲美女搞黄在线观看 | 蜜桃亚洲精品一区二区三区| 老司机福利观看| 成年免费大片在线观看| 99热这里只有精品一区| 中国美白少妇内射xxxbb| 成人亚洲欧美一区二区av| 精华霜和精华液先用哪个| 成年女人永久免费观看视频| 国产亚洲av嫩草精品影院| 国产精品久久久久久久久免| 在线播放国产精品三级| 内射极品少妇av片p| 亚洲欧美精品自产自拍| 美女免费视频网站| 久久中文看片网| 日本爱情动作片www.在线观看 | 亚洲欧美精品综合久久99| 九九爱精品视频在线观看| 日韩大尺度精品在线看网址| 青春草视频在线免费观看| 国产亚洲精品av在线| 亚洲av二区三区四区| 天天躁夜夜躁狠狠久久av| 天天躁夜夜躁狠狠久久av| 国产欧美日韩精品亚洲av| 久久精品人妻少妇| 久久国产乱子免费精品| 神马国产精品三级电影在线观看| 村上凉子中文字幕在线| 菩萨蛮人人尽说江南好唐韦庄 | 成年免费大片在线观看| 国产v大片淫在线免费观看| 久久午夜福利片| 一区二区三区四区激情视频 | 成年女人毛片免费观看观看9| 亚洲av免费高清在线观看| 卡戴珊不雅视频在线播放| 免费电影在线观看免费观看| 国产精品福利在线免费观看| 久久久久久九九精品二区国产| 亚洲经典国产精华液单| 国产精品野战在线观看| 国产亚洲精品综合一区在线观看| 国产欧美日韩精品一区二区| 国产成人freesex在线 | 在现免费观看毛片| 国产精品野战在线观看| 波多野结衣高清无吗| 狂野欧美激情性xxxx在线观看| 亚洲国产精品成人久久小说 | 欧美绝顶高潮抽搐喷水| 亚洲成人中文字幕在线播放| 美女内射精品一级片tv| 中文在线观看免费www的网站| 亚洲av二区三区四区| 久久婷婷人人爽人人干人人爱| 99久久无色码亚洲精品果冻| 中文字幕免费在线视频6| 亚洲欧美日韩东京热| 97碰自拍视频| 国产成人91sexporn| www.色视频.com| 午夜福利18| 日本五十路高清| 3wmmmm亚洲av在线观看| 日韩中字成人| 高清日韩中文字幕在线| 丰满乱子伦码专区| 欧美xxxx黑人xx丫x性爽| 夜夜夜夜夜久久久久| a级一级毛片免费在线观看| 伦理电影大哥的女人| 亚洲精品亚洲一区二区| 成人av一区二区三区在线看| 亚洲精品一区av在线观看| 婷婷色综合大香蕉| 欧美性感艳星| 黄色视频,在线免费观看| 免费av毛片视频| 精品久久久久久久久av| 亚洲七黄色美女视频| 久久午夜亚洲精品久久| 春色校园在线视频观看| 国产日本99.免费观看| 欧美中文日本在线观看视频| 亚洲在线观看片| 又黄又爽又免费观看的视频| 亚洲精品在线观看二区| 五月玫瑰六月丁香| 综合色丁香网| 国产精品一区www在线观看| 人妻久久中文字幕网| 成人综合一区亚洲| 欧美潮喷喷水| 亚洲无线观看免费| 日韩精品中文字幕看吧| 免费无遮挡裸体视频| 22中文网久久字幕| 波多野结衣高清无吗| 在线免费观看的www视频| 亚洲精品国产成人久久av| 成人特级av手机在线观看| 色播亚洲综合网| 在线免费十八禁| 日本黄色片子视频| 成人综合一区亚洲| 美女免费视频网站| 亚洲精品久久国产高清桃花| 看黄色毛片网站| 99国产极品粉嫩在线观看| 午夜亚洲福利在线播放| 真实男女啪啪啪动态图| 久久精品人妻少妇| 69人妻影院| 亚洲不卡免费看| 免费观看在线日韩| 丰满乱子伦码专区| 欧洲精品卡2卡3卡4卡5卡区| 成人亚洲欧美一区二区av| 菩萨蛮人人尽说江南好唐韦庄 | 特级一级黄色大片| 一区福利在线观看| 日韩亚洲欧美综合| 亚洲最大成人中文| 国产精品一及| av女优亚洲男人天堂| 国产在线男女| 久久精品久久久久久噜噜老黄 | 亚洲精品一区av在线观看| 亚洲七黄色美女视频| 高清毛片免费看| 国产精品一区二区性色av| 丰满乱子伦码专区| 成人二区视频| 国产午夜福利久久久久久| 亚洲无线在线观看| 亚洲图色成人| 亚洲真实伦在线观看| 国产高清有码在线观看视频| 特大巨黑吊av在线直播| 国产黄片美女视频| 亚洲精品粉嫩美女一区| 晚上一个人看的免费电影| 国产精品亚洲一级av第二区| 欧美一区二区国产精品久久精品| 成人亚洲精品av一区二区| 国产视频一区二区在线看| 女生性感内裤真人,穿戴方法视频| 最新在线观看一区二区三区| 国产伦一二天堂av在线观看| 国产精品1区2区在线观看.| 日韩欧美一区二区三区在线观看| 国产亚洲精品av在线| 网址你懂的国产日韩在线| 一进一出抽搐gif免费好疼| 精品国产三级普通话版| 国产成人福利小说| 最后的刺客免费高清国语| 国产精品嫩草影院av在线观看| 欧美3d第一页| 成人美女网站在线观看视频| 国产精品一区二区三区四区免费观看 | 日日摸夜夜添夜夜添小说| 欧美日韩综合久久久久久| 日本 av在线| 亚洲va在线va天堂va国产| 在线播放无遮挡| 免费电影在线观看免费观看| 国产久久久一区二区三区| 亚洲av中文av极速乱| 日韩欧美精品免费久久| 黄色日韩在线| 亚洲av第一区精品v没综合| 国产欧美日韩精品一区二区| 我的女老师完整版在线观看| 色尼玛亚洲综合影院| 亚洲自拍偷在线| 简卡轻食公司| 国产高清不卡午夜福利| 69av精品久久久久久| 亚洲中文日韩欧美视频| 亚洲国产欧美人成| 啦啦啦韩国在线观看视频| 丰满的人妻完整版| 99久国产av精品| 亚洲精品日韩在线中文字幕 | 亚洲三级黄色毛片| 国产一区亚洲一区在线观看| 国产乱人偷精品视频| 亚洲精品日韩av片在线观看| 日本撒尿小便嘘嘘汇集6| 人妻夜夜爽99麻豆av| 午夜福利在线观看吧| 内射极品少妇av片p| 国产精品日韩av在线免费观看| 丝袜喷水一区| 一进一出好大好爽视频| 高清毛片免费看| 国产一区二区三区av在线 | 18+在线观看网站| 亚洲av免费在线观看| 九九爱精品视频在线观看| 亚洲自拍偷在线| 午夜激情欧美在线| 亚洲激情五月婷婷啪啪| 国产片特级美女逼逼视频| 在线观看66精品国产| 亚洲av第一区精品v没综合| 久久久欧美国产精品| 亚洲av成人av| 青春草视频在线免费观看| 亚洲av.av天堂| 中文字幕免费在线视频6| 午夜福利18| 国产中年淑女户外野战色| 12—13女人毛片做爰片一| 久久久久九九精品影院| 国产黄片美女视频| 六月丁香七月| 亚洲精品一卡2卡三卡4卡5卡| 天美传媒精品一区二区| 久久精品国产亚洲av涩爱 | 亚洲成a人片在线一区二区| 国产高清视频在线播放一区| 国产精品嫩草影院av在线观看| 99视频精品全部免费 在线| av在线播放精品| 国产精品野战在线观看| 97在线视频观看| 国产色婷婷99| 精品午夜福利在线看| 国产 一区 欧美 日韩| 国产伦精品一区二区三区视频9| 99久国产av精品| 自拍偷自拍亚洲精品老妇| 可以在线观看的亚洲视频| 麻豆久久精品国产亚洲av| 美女cb高潮喷水在线观看| 成人亚洲欧美一区二区av| 久久韩国三级中文字幕| 老熟妇乱子伦视频在线观看| 性插视频无遮挡在线免费观看| 亚洲成av人片在线播放无| 国产一区二区亚洲精品在线观看| 国产在线男女| 中文字幕熟女人妻在线| 真人做人爱边吃奶动态| 波多野结衣高清作品| 乱码一卡2卡4卡精品| 国产高清激情床上av| 精品久久久久久成人av| 男女做爰动态图高潮gif福利片| 久久草成人影院| 国产一区二区三区av在线 | 日本熟妇午夜| 成人午夜高清在线视频| 久久久久国产精品人妻aⅴ院| 亚洲av一区综合| 91在线观看av| 偷拍熟女少妇极品色| 国产伦精品一区二区三区四那| 亚洲av成人av| 天堂av国产一区二区熟女人妻| 亚洲熟妇中文字幕五十中出| 成人欧美大片| 久久久久久大精品| 波多野结衣巨乳人妻| 精品久久久久久久久av| 久久久久久久久中文| 淫秽高清视频在线观看| 日韩成人伦理影院| 精品免费久久久久久久清纯| 国产成人freesex在线 | 欧美一区二区精品小视频在线| 午夜精品在线福利| 国产爱豆传媒在线观看| 韩国av在线不卡| 久久欧美精品欧美久久欧美| 午夜激情福利司机影院| 国产午夜精品久久久久久一区二区三区 | 精品久久久久久久久久免费视频| 国产男人的电影天堂91| 99热精品在线国产| 露出奶头的视频| 午夜视频国产福利| 免费电影在线观看免费观看| 亚洲欧美精品综合久久99| 亚洲国产精品sss在线观看| 在线免费十八禁| 欧美激情国产日韩精品一区| 偷拍熟女少妇极品色| 尾随美女入室| 91在线精品国自产拍蜜月| 精品久久久久久久人妻蜜臀av| 99久久中文字幕三级久久日本| 欧美极品一区二区三区四区| 国产高清激情床上av| 欧美国产日韩亚洲一区| 91麻豆精品激情在线观看国产| 久久人妻av系列| 免费不卡的大黄色大毛片视频在线观看 | 老司机影院成人| 国产高潮美女av| 一本一本综合久久| 日韩精品青青久久久久久| 热99re8久久精品国产| 神马国产精品三级电影在线观看| 国产精品伦人一区二区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产精品久久男人天堂| 99热只有精品国产| 亚洲在线自拍视频| 丰满乱子伦码专区| 亚洲美女搞黄在线观看 | 亚洲高清免费不卡视频| 啦啦啦韩国在线观看视频| 99国产极品粉嫩在线观看| 午夜免费激情av| 国产午夜精品论理片| eeuss影院久久| 成人av在线播放网站| 国产探花极品一区二区| 国产一区亚洲一区在线观看| 麻豆精品久久久久久蜜桃| 一个人看的www免费观看视频| 亚洲内射少妇av| 97超级碰碰碰精品色视频在线观看| 在线观看66精品国产| 亚洲18禁久久av| 国产精品久久久久久久久免| 国产爱豆传媒在线观看| 免费观看人在逋| 国产高清三级在线| 自拍偷自拍亚洲精品老妇| 国产欧美日韩一区二区精品| 波野结衣二区三区在线| 狠狠狠狠99中文字幕| 色综合色国产| 给我免费播放毛片高清在线观看| 亚洲成人久久性| 狂野欧美白嫩少妇大欣赏| 亚洲乱码一区二区免费版| 久久草成人影院| 国产伦精品一区二区三区视频9| 观看免费一级毛片| 丰满乱子伦码专区| 日韩,欧美,国产一区二区三区 | 搡老妇女老女人老熟妇| av天堂中文字幕网| 别揉我奶头~嗯~啊~动态视频| 午夜福利在线在线| 中文字幕免费在线视频6| 亚洲av.av天堂| 联通29元200g的流量卡| 国产精品女同一区二区软件| 日本熟妇午夜| 美女cb高潮喷水在线观看| 午夜福利高清视频| 最新在线观看一区二区三区| 婷婷色综合大香蕉| 久久人妻av系列| 国产一区二区激情短视频| 日韩欧美在线乱码| 免费人成视频x8x8入口观看| 不卡一级毛片| 三级毛片av免费| 日韩欧美精品v在线| 日韩成人伦理影院| 久久99热这里只有精品18| 晚上一个人看的免费电影| 日产精品乱码卡一卡2卡三| 国产人妻一区二区三区在| 国产极品精品免费视频能看的| 国产成人a∨麻豆精品| 少妇熟女aⅴ在线视频| 激情 狠狠 欧美| 看片在线看免费视频| 国产国拍精品亚洲av在线观看| 亚洲天堂国产精品一区在线| 菩萨蛮人人尽说江南好唐韦庄 | 国产一区二区三区av在线 | 在线天堂最新版资源| 国产成人aa在线观看| 无遮挡黄片免费观看| 日本三级黄在线观看| 国产欧美日韩精品亚洲av| 亚洲无线观看免费| 日韩欧美精品免费久久| 欧美一区二区精品小视频在线| 免费看av在线观看网站| 丰满的人妻完整版| 少妇的逼水好多| 一夜夜www| 在线天堂最新版资源| 51国产日韩欧美| 18禁黄网站禁片免费观看直播| 国产大屁股一区二区在线视频| 国产精品日韩av在线免费观看| 免费看光身美女| 国产成人影院久久av| 天天躁夜夜躁狠狠久久av| 如何舔出高潮| 亚洲精品色激情综合| 国产亚洲精品av在线| 美女内射精品一级片tv| 久久热精品热| 欧美成人免费av一区二区三区| 久久久久久久久大av| 看片在线看免费视频| 欧美激情在线99| 色哟哟哟哟哟哟| 性欧美人与动物交配| 久久久久久久久大av| 少妇的逼水好多| 一个人观看的视频www高清免费观看| 午夜福利在线观看免费完整高清在 | 精品午夜福利视频在线观看一区| 亚洲18禁久久av| 久久久久久久久久成人| 男插女下体视频免费在线播放| 国产av在哪里看| 又爽又黄a免费视频| 久久久久久大精品| av福利片在线观看| 久久人人爽人人片av| 国产高清激情床上av| 日日摸夜夜添夜夜添小说| 岛国在线免费视频观看| 99久久久亚洲精品蜜臀av| 极品教师在线视频| 长腿黑丝高跟| 亚洲国产精品合色在线| 成人av在线播放网站| 国产不卡一卡二| 精品国内亚洲2022精品成人| 一区二区三区四区激情视频 | 国产一区二区激情短视频| 国产毛片a区久久久久| 精品一区二区三区视频在线| eeuss影院久久| 青春草视频在线免费观看| 91在线精品国自产拍蜜月| 国产精品一区二区性色av| 99久国产av精品| 免费高清视频大片| 不卡视频在线观看欧美| 精品一区二区三区视频在线| 看非洲黑人一级黄片| 精品国内亚洲2022精品成人| 欧美激情在线99| 久久热精品热| 国产精品永久免费网站| 91久久精品国产一区二区三区| 波野结衣二区三区在线| 韩国av在线不卡| 色吧在线观看| 色5月婷婷丁香| 99久久精品国产国产毛片| 亚洲欧美日韩无卡精品| 人妻丰满熟妇av一区二区三区| 狠狠狠狠99中文字幕| 国产高潮美女av| 九九爱精品视频在线观看| 亚洲精品国产成人久久av| 成人亚洲欧美一区二区av| 国产精品一二三区在线看| 亚州av有码| 不卡视频在线观看欧美| 别揉我奶头~嗯~啊~动态视频| 久久精品久久久久久噜噜老黄 | 久久久精品大字幕| 欧美zozozo另类| 大型黄色视频在线免费观看| 久久精品国产亚洲av天美| 嫩草影院精品99| 免费观看在线日韩| 欧美成人免费av一区二区三区| 国产一区二区激情短视频| av福利片在线观看| 美女高潮的动态| 亚洲av成人av| 国产黄片美女视频| 在线看三级毛片| 国产精品人妻久久久影院| 午夜日韩欧美国产| 一本久久中文字幕| 五月伊人婷婷丁香| 免费av观看视频| 日韩欧美精品免费久久| 午夜福利18| 免费人成在线观看视频色| 欧美成人a在线观看| av在线观看视频网站免费| 天堂影院成人在线观看| 97在线视频观看| 九九在线视频观看精品| 亚洲精品456在线播放app| 国产一区二区激情短视频| 女同久久另类99精品国产91| 干丝袜人妻中文字幕| 九色成人免费人妻av| 网址你懂的国产日韩在线| 国产精品一二三区在线看| 国产男人的电影天堂91| h日本视频在线播放| 国产亚洲欧美98| 亚洲熟妇熟女久久| 小说图片视频综合网站| 亚洲精品456在线播放app| 久久人人爽人人爽人人片va| 伦理电影大哥的女人| 麻豆精品久久久久久蜜桃| 亚洲av成人av| 变态另类成人亚洲欧美熟女| 国产单亲对白刺激| 国产一区二区在线av高清观看| 国产精品亚洲一级av第二区| 国产v大片淫在线免费观看| 国产一区二区亚洲精品在线观看| 国产亚洲精品久久久久久毛片| 女人被狂操c到高潮| 有码 亚洲区| 全区人妻精品视频| 老司机福利观看| 最近视频中文字幕2019在线8|