1. 引言
流体的粘性在体育竞技比赛尤其是游泳比赛中具有显著的影响,张晓侠 [1] 综述了运用CFD模型建立运动员在游泳时的模型流体对运动员的作用力,肯定了CFD仿真技术能对游泳运动的技术指导及相关的流体科学研究进行了很好的解释,避免了高成本和复杂过程的试验;Gardano和Dabnichki [2] 在模拟自由泳手部划水时的动作,通过比较肘部与上臂不同夹角的阻力系数,发现肘部与上臂夹角为160˚时,阻力最大;Zaidi [3] 等人,利用CATIA软件建立二维测试图,采用GAMBIT划分三角非结构泳池结构,分析游泳时身体各处涡流的分布情况,发现表面剪切阻力随着人体与水流相对速度的增大而增大,人体表面几何复杂位置的剪切压力较其他部位更大;Marinho [4] 等人利用FLUENT软件模拟稳流划水运动中运动员手臂不同位置模型分析了行进阻力系数变化,发现当手臂与运动员身体更为贴近时,行进阻力更小,游泳的速度越快;Novais [5] 等人为分析人体在不同速度对滑行特性的影响,通过利用CFD模型仿真运动员在3种行进速度(1.5 m.s,2.0 m/s和2.5 m/s)、水深为0 m时的阻力情况,得出与水面正面接触越小,其形阻与摩阻越小的结论,另外,通过将运动员的运动深度放到水下,在做比较发现在水深0.25 m时阻力及阻力系数达到最大,在水深超过0.25 m后,阻力及阻力系数随水深逐渐下降,并得到0.75 m是综合考虑运动员的推进力和体力得到的综合结果;为研究多个运动员游泳时的水流特性,Sliva [6] 利用二维CFD模型,通过将运动员设置在不同间隔距离的情况下,对运动员周身的阻力系数做出比较,发现在两人间距0.50~6.00 m时,领先运动员的阻力随水流速度的增大而减小,而落后的运动员身边的阻力系数约为前者的56%,而当两人相距6.45~8.90 m,水流对两人的阻力系数不会受到干扰;另外Hanna [7] 为了获取更先进的游泳装备,同样可以采用CFD模型对身体的仿真来评价不同泳衣的材质和特性,减少了研发成本和研发时间。上述游泳运动的CFD研究,对体育运动中流体的研究起到了推进作用。
我们可以将这种方式运用到游泳池中,通过运用CFD仿真游泳赛道,分析仿真结果,得到更公平公正的赛道,或者在现有的赛道上得到更有利的游泳策略,为游泳比赛取得理论性的优势。
2. 流体粘性在游泳比赛上的理论分析
在游泳比赛的预选赛中,成绩拔尖的几名选手将会在决赛中分配在中央泳道,而成绩不太好的则会被分配到游泳池两侧靠近游泳池壁的泳道,这样做显然是有原因的。我们可以通过有限元的方法分析各个游泳者附近的水流速,并想办法改进游泳池的结构,使比赛更具公平性,下面我们首先简化物理模型,然后通过理论分析粗略计算粘性阻力。
2.1. 游泳者所受粘性阻力的前提假设
为了利用数学的方式来解决这个问题,我首先假设一些必须条件。
假设一,流体的各个物理量满足连续介质模型。
假设二,室温25˚,水的动力粘度μ = 0.894 × 10−3 Pa∙s。
假设三,游泳者的游动速度较低,且由于Re = ρvl/μ,从而使得雷诺数Re < 2000,可判断流体为层流。
假设四,共有7个泳道,且每个人在相同的游泳环境下游泳速度是相同的,均假设为v。
假设五,泳道的宽度为2y。
假设六,游泳者的游泳方式一样,即与侧身水流的接触面积一样,均为A。
假设七,由于隔离带的水深方向的长度较小,因而其对游泳者的影响较小。
2.2. 在中间的游泳者的粘性阻力分析
由假设可知,在游泳者游动过程中,由于隔离带在水深方向的尺寸不大,因而其对流体的传递流速的影响不大,可忽略隔离带的影响,于是游泳者附近的流体在粘性的作用下存在与游泳者游动方向相同的流速方向,即水流的向前速度在泳道宽度方向上是等梯度分布的,且由于水的动力粘度较大,因而可粗略地假设:在隔离带附近的流体流速为δV,δV > 0,由于速度是等梯度分布的,于是隔离带上水流与
游泳者侧身的水流的速度梯度
;
因而游泳者与左隔离带水流的粘性阻力为:
;
因而左右两隔离带对中央游泳者所产生的粘性阻力为:
游泳比赛粘性阻力
。
2.3. 在两侧泳道的游泳者受到的粘性阻力分析
由假设一可知,在两侧靠近游泳池壁的水的流速V1 = 0,游泳者与左侧隔离带间的流体的粘性阻力易知为:
;
游泳者与右侧游泳池壁间的水流的速度梯度为:
;
因而游泳者与右侧流体流层间的粘性阻力为:
;
因而左右水流对两侧游泳者所产生的粘性阻力为:
;
所以综上分析可知,在上述假设成立的情况下,即隔离带对水流的减速作用不明显时,显然中央游泳者所受到的粘性阻力近似等于两侧游泳者所受粘性阻力,但仍稍微小,即:
。
因而游泳比赛的不公平性得以体现。
3. 仿真模拟分析
3.1. 三维模型
为了进行实验分析,首先简化数学模型。因此在上述理论分析的建设的前提下,为了精确计算粘性阻力,应对某些模糊值量化。模型假设:为了凸显粘性阻力,有两种方法。一是增大粘性阻力,二是粘性阻力与其他阻力(手拍击水面的反力等)共同消耗水流的动能。因此可以减少其他阻力,以使粘性阻力所占比例增大,从而使得实验结果明显。于是游泳者在游泳池内潜泳,且假设其游泳面积为一圆形,半径为r = 0.2 m,并且增大游动速度。为了便于Autodesk CFD模拟计算,可利用直径为0.4 m的水泵在游泳池一端提供较大流速的水流,通过观察游泳池内流体分布,计算粘性阻力。并判断不同泳道的粘性阻力相对大小。量值假设:由上假设游泳者两侧与水流接触面积为半圆柱面,所以A = π × r × h = 1.131 m2 (人身高可设为1.8 m,即h = 1.8),游泳者假设圆心距离水面为0.45 m,水深为2 m,泳道宽度为2 m。共有七个泳道,即需七个水泵提供动能。由水泵及池壁组成的泳池网格化三维模型如下图1~4所示。
图1. 网格数为239,090的三维模型
Figure 1. 3D model with 239,090 grids
![](//html.hanspub.org/file/31-2570515x16_hanspub.png?20220524093020983)
Figure 2. 3D model with 187,534 grids
图2. 网格数为187,534的三维模型
![](//html.hanspub.org/file/31-2570515x17_hanspub.png?20220524093020983)
Figure 3. 3D model with 152,236 grids
图3. 网格数为152,236的三维模型
![](//html.hanspub.org/file/31-2570515x18_hanspub.png?20220524093020983)
Figure 4. 3D model with 136,260 grids
图4. 网格数为136,260的三维模型
如图5不同网格数量的粘性阻尼对比图所示,将上述4种方案进行计算泳池粘性阻尼的仿真中可以看出,后2种方案的结果相差是可以接受的,而前两种方案的结果要低得多。为了简化研究过程,在保证精度的前提下尽可能缩短研究时间,选择了第三种方案作为模型网格划分的依据,其他泳道在网格划分阶段也采用了相同的方案。
![](//html.hanspub.org/file/31-2570515x19_hanspub.png?20220524093020983)
Figure 5. Comparison of viscous damping with different mesh numbers
图5. 不同网格数量的粘性阻尼对比图
3.2. CFD模拟
进入Autodesk CFD,并导入Solidworks模型。且设置流体为水,温度为常温,水泵的流量为0.65 m3/min,转速为2000 r/min。然后自动网格划分,经过求解,得出以下分析和计算结果。
为了使立体的结果,平面化,即建立如下平面,过各水泵中心线,且与池底相垂直的各平面依照与左池壁的距离分别设置为A1,A2,A3,A4,A5,A6,A7,A8平面,过水泵中心线,且与池底平面相平行的平面,设为B平面,与A平面和B平面垂直,且距离泵芯距离50 cm的平面设置为C1平面,同时距离为100 cm,150 cm,200 cm,250 cm,300 cm的平面,依次设置为C2,C3,C4,C5,C6平面。
B平面流体流速如下图6所示,中间位置的流速比两边的流速较大。
![](//html.hanspub.org/file/31-2570515x20_hanspub.png?20220524093020983)
Figure 6. B plane absolute flow velocity distribution diagram
图6. B平面流体绝对流速分布图
下图7所示,由上至下分别为A1,A2,……,A7平面,且各平面上的流体绝对流速分布图如下所示。
![](//html.hanspub.org/file/31-2570515x21_hanspub.png?20220524093020983)
Figure 7. A plane absolute flow velocity distribution diagram
图7. A平面流体绝对流速分布图
为了计算水泵中心线上流体的绝对速度分布,我们已在上面建好了C1,C2,……,C7平面,再求出C平面上某点的绝对速度。如图8所示。
![](//html.hanspub.org/file/31-2570515x22_hanspub.png?20220524093020983)
Figure 8. Analysis diagram of taking point process
图8. 取点过程分析图
从而在C平面内得出各泵体中心点连线处流体绝对流速的分布情况,如下图9~14所示。
![](//html.hanspub.org/file/31-2570515x23_hanspub.png?20220524093020983)
Figure 9. Velocity analysis in C1 plane
图9. C1平面流速分析
将相关数据导入Excel中,并挑选出有用数据,如下表1和表2所示。
![](Images/Table_Tmp.jpg)
Table 1. Absolute velocity of water pump centerline at plane C
表1. 水泵中心线在C平面处的绝对速度
![](Images/Table_Tmp.jpg)
Table 2. Absolute velocity of the center point of the two water pumps at plane C
表2. 两水泵中心点在C平面处的绝对速度
注:C(1~6):表示数据点所在平面;m(1~7):各水泵中心线与C平面的交点;n(1~7):表示两水泵中心处与C平面的交点。
作出从左至右,1~7泳道,在其各自泳道水泵中心线方向的流体的绝对速度分布,如下图15所示。
![](//html.hanspub.org/file/31-2570515x29_hanspub.png?20220524093020983)
Figure 15. Absolute velocity distribution curve of water pump centerline
图15. 水泵中心线绝对速度分布曲线
作出从C1至C6平面,1~7泳道,在C平面的绝对速度分布如下图16所示。
![](//html.hanspub.org/file/31-2570515x30_hanspub.png?20220524093020983)
Figure 16. The absolute velocity distribution curve in the C plane
图16. C平面上绝对速度分布曲线
3.3. 数据分析
因为在里约奥运会男子200米比赛中,孙杨用时1分44秒夺冠,平均游泳速度为1.9231 m/s,因此可用上述C7平面进行分析。且由CFD的模拟结果,若游泳者的体质相同,即在游泳过程中可以产生相同的驱动力,则各泳道游泳者由于不同的粘性阻力,造成其所能达到的游泳速度不同。
若4泳道游泳者的速度为:v4 = 2.15725 m/s2。则1至7泳道游泳者的速度为:v1 = 1.47513 m/s2,v2 = 1.13232 m/s2,v3 = 1.12453 m/s2,v5 = 1.00036 m/s2,v6 = 0.810014 m/s2,v7 = 1.13051 m/s2。
且各泳道游泳者所受到的粘性阻力,由公式F = Aμdv/dy可求出,即从1至7泳道游泳者所受到的粘性阻力为:
。
同理可求,F2 = 0.00159 N,F3 = 0.00095 N,F4 = 0.00300 N,F5 = 0.00123 N,F6 = 0.00103 N,F7 = 0.00198 N。
其粘性阻力如下图17所示:
![](//html.hanspub.org/file/31-2570515x32_hanspub.png?20220524093020983)
Figure 17. Viscous resistance curves of each swimming lane
图17. 各泳道粘性阻力曲线
由上分析可知,中间泳道所受到的粘性阻力最大,但从CFD模拟求解结果可知,具有相同初动能的水流,在相同时刻,中间泳道的水流的绝对流速最大,即克服粘性阻力所做功小,但中间泳道的粘性阻力又是最大的,这两种条件显然是矛盾的。因而我猜测各泳道的粘性阻力与各泳道间的速度差值有关,即随着比赛时间的增大,各泳道速度差值渐增,从而导致中间泳道粘性阻力增大。① 比赛起始,各游泳者游动速度较低,受粘性阻力影响较低,因此各游泳者的速度近似相等,但由于两游泳者中间的流速不为0,从而致使中间游泳者的速度梯度小,粘性阻力较两侧的小。② 在游泳者加速过程中,由于中间游泳者所受到的粘性阻力小,其游泳速度升高较快,导致各游泳者的游动速度的差值增大,从而致使中间游泳者的速度梯度增大,其粘性阻力逐渐增大,即两侧的游泳者对中间的游泳者的游动起到阻碍的效果。③ 由于粘性阻力的增大,中间游泳者的游动速度将会下降。最后导致游泳者的游动速度之比趋于定值。即中间泳道的粘性阻力先逐渐增大,达到一定值时,粘性阻力在此值上下波动。
为了验证猜想,对实验数据加工,若下表3所示,求出各水泵中心线在C平面处的粘性切应力。
![](Images/Table_Tmp.jpg)
Table 3. Viscous shear stress of each swimming lane
表3. 各泳道粘性切应力
以上数据为基础,泳道速度相对差值从0开始增大,且假设7泳道粘性切应力为1,求出其他泳道的相对粘性切应力如下表4所示。
![](Images/Table_Tmp.jpg)
Table 4. Relative viscous shear stress of each swimming lane
表4. 各泳道相对粘性切应力
则中间泳道的相对粘性切应力如下图18所示;且下曲线符合假设内容,假设正确。
![](//html.hanspub.org/file/31-2570515x33_hanspub.png?20220524093020983)
Figure 18. Relative viscous shear stress curve of middle swimming lane
图18. 中间泳道相对粘性切应力曲线
4. 游泳比赛场地的改进
1) 由牛顿粘性应力公式,易知,最直接也是最可靠的减小粘性应力的方式是减小其流层间的速度梯度,但由上分析游泳者所受粘性阻力时,由于隔离带在水深方向的尺寸不大,因而其对流体的传递流速的影响不大,可忽略隔离带的影响,但若将隔离带换为隔离层的话,即隔离层的水深方向的尺寸大时,隔离带与池壁所起到的作用相同,从而致使隔离带附近流速为0,即将一个大游泳池在结构上划分为7个小游泳池,他们之间近似互不影响,且显然此时不论游泳者在哪个泳道,其所受到的粘性阻力近乎相同,粘性阻力对比赛成绩的影响很小,因而比赛更显公平。
2) 在游泳池壁附近增加传送带类装置,其装置如下图所示。且当游泳者向左游动时,隔离带随着传动系统向左移动,当游泳者向右游动时,隔离带随着传动系统向右移动,通过减小游泳者与游泳池壁间的速度梯度来减小各泳道粘性阻力差值。如图19所示为改进后的泳池结构。
![](//html.hanspub.org/file/31-2570515x34_hanspub.png?20220524093020983)
Figure 19. Conveyor belt type swimming pool wall
图19. 类传送带装置式泳池壁
5. 结论
总的来说,粘性阻力对中间游泳者起有利作用。中间泳道所受到的粘性阻力最大,但从CFD模拟求解结果可知,具有相同初动能的水流,在相同时刻,中间泳道的水流的绝对流速最大,即克服粘性阻力所做功小,但中间泳道的粘性阻力又是最大的,这两种条件显然是矛盾的。因而各泳道的粘性阻力与各泳道间的速度差值有关,即随着比赛时间的增大,各泳道速度差值渐增,从而导致中间泳道粘性阻力增大。① 比赛起始,各游泳者游动速度较低,受粘性阻力影响较低,因此各游泳者的速度近似相等,但由于两游泳者中间的流速不为0,从而致使中间游泳者的速度梯度小,粘性阻力较两侧的小。② 在游泳者加速过程中,由于中间游泳者所受到的粘性阻力小,其游泳速度升高较快,导致各游泳者的游动速度的差值增大,从而致使中间游泳者的速度梯度增大,其粘性阻力逐渐增大,即两侧的游泳者对中间的游泳者的游动起到阻碍的效果。③ 由于粘性阻力的增大,中间游泳者的游动速度将会下降。最后导致游泳者的游动速度之比趋于定值。即中间泳道的粘性阻力先逐渐增大,达到一定值时,粘性阻力在此值上下波动。
从以上分析中,我们可得出结论,比赛初始阶段,中间泳道粘性阻力较小,粘性阻力的影响对中间游泳者的起步加速是有利的。但若中间游泳者处于领先时,两侧的游泳者对其游动起阻碍作用,从而增大中间游泳者的粘性阻力,于是粘性阻力的影响对中间游泳者的稳定游动阶段是不利的。但当中间游泳者处于落后时,粘性阻力又起有利作用。