高尔夫球射程问题.docVIP

第３４卷　第４期武汉理工大学学报（交通科学与工程版）Ｖｏｌ．３４　Ｎｏ．４２０１０年８月Ｊｏｕｌ　ｏｆ　Ｗｕｈａｎ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　ＴｅｃｈｌｏｇｙＡｕｇ．２０１０（Ｔ(ａ) (ｎ) (ｒ) ｒａｎｓｐｏｒｔａｔｉｏｎ　Ｓｃｉｅｎｃｅ　＆Ｅｎｇｉｎｅｅｅｒｉｎｇ）(ｏ) (ｎ) 高尔夫球射程问题的动力学模拟研究＊ ｕ 高明向１）　陈定方２） （武汉理工大学理学院１）　武汉　４３００６３）　（武汉理工大学智能制造与控制研究所２）　武汉　４３００６３） 摘要：利用数值模拟方法研究了高尔夫球的射程问题．通过对初始速率、旋转速度及抛射角等多种 能够人为控制的因素对射程影响的模拟与分析，得到了这些因素影响射程的规律． 关键词：高尔夫球；射程；数值模拟 中图法分类号：Ｏ３１３．３ＤＯＩ：１０．３９６３／ｊ．ｉｓ ｓｓｎ．１００６－２８２３．２０１０．０４．０４８ ０　引　　言量，绕过质心的轴转动的角速度为ω；ｇ为重力加 速度；ρ为空气密度；Ｃｄ 为空气的阻力系数；μ为 高尔夫球是一种把享受大自然乐趣、体育锻空气的动态粘滞系数；Ａ为球的特征面积；Ｖ为球 炼和游戏集于一身的运动，是一种在优美环境中的体积；Ｆｗ 为风对球的作用力，也是位置和时间 进行的高尚娱乐活动．对高尔夫球类运动规律的的函数，可以根据实际需要建立不同的风力模型． 研究，将有助于高尔夫球运动的发展．在国外，对则球在空气中平动的常微分方程为 （ ｕ ｘ，ｔ）· 高尔夫球运动规律的研究方法除理论和实验研究 ¨ ＾ｕ （ｘ，ｔ）Ａ ｇ＋Ｃｄρ ＋ ｍｘ＝ｍ ２ 外，还采用数值模拟方法［１］．在国内则侧重于理论 ８πｒ ３ρω（ｘ，ｔ）×ｕ（ｘ，ｔ）＋ρＶｇ＋Ｆｗ（ｘ，ｔ） 分析［２－３］，而在实验研究和数值模拟上做得还不 ３ ＾ｕ 够．在这些文献中，为了分析和表示的方便，多将 一些对球的运动影响较小的因素忽略不计，并把 空气的阻力系数当成常量处理，因而得到的理论 分析结果和实际情况有一定差别． 尽可能多地考虑影响高尔夫球运动的因素， 建立高尔夫球运动微分方程，利用数值模拟方法， 对高尔夫球的运动进行动力学模拟，可以最大限 度地获得其真实的动力学行为，无论对于工程应 用还是游戏娱乐都具有很大的实际意义．本文针 对球员能够控制且对射程有明显影响的因素进行 了动力学模拟与分析． １　高尔夫球的运动微分方程 设高尔夫球的质量为ｍ；半径为ｒ；质心速度 为ｕ，是位置ｘ和时间ｔ的函数；为单位速度矢 （１）　　在空气中转动的常微分方程为 １θ(¨) ＝－８πｒ３μω（ｘ，ｔ）（２） 　　方程（１）右侧第１项为球的重力Ｇ；第２项为 空气阻力Ｆｄ；第３项为马格纳斯升力Ｆｌ；第４项 为空气浮力Ｆｆ；方程（２）右侧为空气对球的阻力 矩τ． 下面根据此２方程对高尔夫球的射程问题进 行数值模拟，并对结果进行分析．下文提到的粗糙 球是指现代高尔夫球（网纹高尔夫球）． ２　阻力与光滑程度和球速的关系 当高尔夫球在空气中飞行时，由于其速度比 较高，空气阻力所产生的影响就显得十分明显．对 于大小相同、初始速率相等的球来说，其所受空气 　　　　收稿日期：２０１０－０４－０８ 　　　　高明向（１９７０－）：男，博士，讲师，主要研究领域为虚拟现实技术的理论与实现 　　　＊湖北省科技攻关项目资助（批准号：２００４ＡＡ１０５Ｂ０４）  　第４期高明向，等：高尔夫球射程问题的动力学模拟研究·８３９· 阻力的差异表现在Ｃｄ 的不同，Ｃｄ 则与球表面粗 糙程度以及雷诺数Ｒｅ有关． 图１反映了Ｒｅ＞４×１０４时，粗糙球与光滑球 的Ｃｄ 与Ｒｅ间的关系．当Ｒｅ＜４×１０４时，由于光 滑球与粗糙球表面形成的边界层都是层流，边界 层的分离点大致相同，压差阻力也差不多一样．但 是，由于粗糙的表面会增加摩擦阻力，因而此时粗 糙球也会比光滑球受到更大的阻力． 图１　粗糙球与光滑球的Ｃｄ 与Ｒｅ间的关系 如图１所示，对于光滑球，当Ｒｅ超过３×１０５ 时Ｃｄ 突然减小，对于粗糙球，当Ｒｅ超过０．５×１０５时开始减少，而接近０．７５×１０５时Ｃｄ 达到最 小．优秀高尔夫球运动员的击球速度可达到６１ｍ／ｓ，可以计算，对于实际的高尔夫球，在１标准 大气压下，２０ ℃时，此速度下，Ｒｅ＝１．７３×１０５，根 据图１，在此雷诺数下空气对粗糙球的阻力系数 约为０．２６，而对光滑球的阻力系数约为０．４８（图 中已用虚线标出）．因此粗糙球受到的空气阻力要 比光滑球小得多．但速度为３１ｍ／ｓ时，Ｒｅ＝０．７５ ×１０５，根据图１，空气对二者的阻力系数相差 更多． 随着雷诺数的增加，当超过４×１０５时，从图 １中看到，空气对粗糙球的阻力系数反而超过了 光滑球．此时的球速大约为１４０ｍ／ｓ，一般情况 下，很难打出这么高的速度． ３　高尔夫球射程的动力学