(57)【要約】【解決手段】格子を生成させた気流仮想空間内にゴルフボールモデルを設定し、該ボールモデルに重量及び打ち出し初期条件を与えて飛行させ、各格子区画で算出した気体流の速度、方向及び圧力から、揚力係数及び抗力係数を算出し、飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を所定の微少時間毎に算出して、ゴルフボールの飛距離、左右へのばらつきを算出する場合に、レイノルズ数により飛行行程を２つ又は３つの領域に分け、ゴルフボールモデル表面と接する上記格子の高さを各領域毎に設定することにより、計算回数を少なくして計算時間の短縮を図る。【効果】実物モデルを用いる必要なく、ゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときの弾道を推定することができ、しかも推定を行なう計算時間を短縮することができる。【選択図】図１
2011224344.tif 000002
【特許請求の範囲】
【請求項１】
表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する弾道シミュレーション方法であって、
（Ａ）ボールの打出される気流仮想空間内（フィールド）に、格子を生成し、
（Ｂ）表面に複数のディンプルが形成された略球状のゴルフボールモデルを気流仮想空間内（フィールド）に設定し、
（Ｃ）その際、ゴルフボールモデルの近傍は、格子を仮想空間内（フィールド）の格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画し、
（Ｄ）ゴルフボールモデルに重量を設定し、
（Ｅ）ゴルフボールモデルに初期条件（初期速度、打出し角度、回転数）を与え、
（Ｆ）ゴルフボールモデルが初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内（フィールド）の格子に戻し、
（Ｇ）気流仮想空間内（フィールド）の各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内（フィールド）における飛行中のゴルフボールの揚力係数及び抗力係数を算出し、
（Ｈ）飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定する際に、
飛行中のゴルフボールモデル周囲のレイノルズ数（Ｒｅ）を算出すると共に、Ｒｅ＝２５００００〜１７００００を高速領域、Ｒｅ＝１７００００〜１０００００を中速領域、及びＲｅ＝１０００００以下を低速領域とする基準に従って、ゴルフボールモデルの飛行行程を２つ又は３つの領域に分け、上記高速領域におけるゴルフボールモデル表面と接する上記格子の高さ（格子高さＡ）、中速領域におけるゴルフボールモデル表面と接する上記格子の高さ（格子高さＢ）及び低速領域におけるゴルフボールモデル表面と接する上記格子の高さ（格子高さＣ）を、格子高さＡ＜格子高さＢ＜格子高さＣとなるように設定することを特徴とするゴルフボールの弾道シミュレーション方法。
【請求項２】
（Ｄ−２）上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定する請求項１記載のゴルフボールの弾道シミュレーション方法。
【請求項３】
上記高速領域における上記格子高さＡを０．００１〜０．０１５ｍｍの範囲、中速領域における上記格子高さＢを０．００５〜０．０２ｍｍの範囲、低速領域における上記格子高さＣを０．０１〜０．０２５ｍｍの範囲でそれぞれ設定する請求項１又は２記載のゴルフボールの弾道シミュレーション方法。
【請求項４】
上記気流仮想空間（フィールド）の下面に地面を想定した壁面を設定して、該気流仮想空間（フィールド）の下面を閉鎖面とする請求項１〜３のいずれか１項に記載のゴルフボールの弾道シミュレーション方法。
【請求項５】
上記（Ｅ）で設定する初期条件に、地面を想定した上記気流仮想空間（フィールド）の下面からゴルフボールモデル最下点までの距離を打撃高さとして加える請求項４記載のゴルフボールの弾道シミュレーション方法。
【請求項６】
格子形状が直交適合格子である請求項１〜５のいずれか１項に記載のゴルフボールの弾道シミュレーション方法。
【請求項７】
格子形状が非構造格子である請求項１〜５のいずれか１項に記載のゴルフボールの弾道シミュレーション方法。
【請求項８】
算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離を算出し、打出しから静止するまでの飛距離を算出する請求項１〜７のいずれか１項に記載のゴルフボールの弾道シミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、表面に複数のディンプルが形成されたゴルフボールの弾道を、コンピューター上にゴルフボールモデルを設定し、コンピューターによる演算を利用してゴルフボールモデルの運動要素を算出して推定する弾道シミュレーション方法に関する。
【背景技術】
【０００２】
ゴルフボールのような物体が大気中で飛行するときには、物体の周囲に気流の乱れが生じることが知られている。物体の表面が複雑な形状を形成していたり、回転しながら飛行していたりすると、飛行中の気流の乱れは複雑になり、物体の飛距離など飛び性能に大きな影響が及ぶこととなる。
【０００３】
ゴルフボールについては、多くの場合、平面視円形のディンプルが多数設けられているが、ディンプルの立体的な形状、配置、大きさなどの組合せが空力特性に及す影響が大きいことから、これらのディンプル要素と空力特性との因果関係を把握する必要がある。
【０００４】
通常、ディンプルの形状、構造及び配置等を変化させることによるゴルフボールの飛び性能に与える影響を調べるに当たっては、ボールの各種成形金型を作製してボールを種々試作し、ボールの打撃実験を行って、初速、スピン、弾道（飛距離、高さ）等を測定し、空力特性を評価することが多かった。
【０００５】
しかしながら、このような実物モデルによる試験評価は、多大な時間と費用を必要とすることに加えて、ディンプルの形状及び配置と空力特性との因果関係を明確に関連づけることができない。そのため、実験によって得られた評価結果に基づき新たに設計したゴルフボールが狙いどおりの性能を発揮するに至らないことがしばしば起こる。このような場合、その都度ボールの設計及び試作をやり直して、空力特性を確かめる必要があり、時間と費用が更に発生するため、効率的にボールの開発ができないという問題がある。
なお、本発明に関する先行技術文献情報としては、以下のものがある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００２−３５８４７３号公報
【特許文献２】特開２００６−２７５７２２号公報
【特許文献３】特開２００５−０３４３７８号公報
【特許文献４】特開２００２−３４０７３５号公報
【特許文献５】特開２００２−２５０７３９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明は、このような問題に鑑みなされたもので、表面に複数のディンプルが形成されたゴルフボールの弾道を実物モデルによる試験評価によらずに評価して、効率的にゴルフボールを開発する方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明は、上記目的を達成するため、下記〔１〕〜〔８〕のゴルフボールの弾道シミュレーション方法を提供する。
〔１〕表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する弾道シミュレーション方法であって、
（Ａ）ボールの打出される気流仮想空間内（フィールド）に、格子を生成し、
（Ｂ）表面に複数のディンプルが形成された略球状のゴルフボールモデルを気流仮想空間内（フィールド）に設定し、
（Ｃ）その際、ゴルフボールモデルの近傍は、格子を仮想空間内（フィールド）の格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画し、
（Ｄ）ゴルフボールモデルに重量を設定し、
（Ｅ）ゴルフボールモデルに初期条件（初期速度、打出し角度、回転数）を与え、
（Ｆ）ゴルフボールモデルが初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内（フィールド）の格子に戻し、
（Ｇ）気流仮想空間内（フィールド）の各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内（フィールド）における飛行中のゴルフボールの揚力係数及び抗力係数を算出し、
（Ｈ）飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定する際に、
飛行中のゴルフボールモデル周囲のレイノルズ数（Ｒｅ）を算出すると共に、Ｒｅ＝２５００００〜１７００００を高速領域、Ｒｅ＝１７００００〜１０００００を中速領域、及びＲｅ＝１０００００以下を低速領域とする基準に従って、ゴルフボールモデルの飛行行程を２つ又は３つの領域に分け、上記高速領域におけるゴルフボールモデル表面と接する上記格子の高さ（格子高さＡ）、中速領域におけるゴルフボールモデル表面と接する上記格子の高さ（格子高さＢ）及び低速領域におけるゴルフボールモデル表面と接する上記格子の高さ（格子高さＣ）を、格子高さＡ＜格子高さＢ＜格子高さＣとなるように設定することを特徴とするゴルフボールの弾道シミュレーション方法。
〔２〕（Ｄ−２）上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定する〔１〕記載のゴルフボールの弾道シミュレーション方法。
〔３〕上記高速領域における上記格子高さＡを０．００１〜０．０１５ｍｍの範囲、中速領域における上記格子高さＢを０．００５〜０．０２ｍｍの範囲、低速領域における上記格子高さＣを０．０１〜０．０２５ｍｍの範囲でそれぞれ設定する〔１〕又は〔２〕記載のゴルフボールの弾道シミュレーション方法。
〔４〕上記気流仮想空間（フィールド）の下面に地面を想定した壁面を設定して、該気流仮想空間（フィールド）の下面を閉鎖面とする〔１〕〜〔３〕のいずれかに記載のゴルフボールの弾道シミュレーション方法。
〔５〕上記（Ｅ）で設定する初期条件に、地面を想定した上記気流仮想空間（フィールド）の下面からゴルフボールモデル最下点までの距離を打撃高さとして加える〔４〕記載のゴルフボールの弾道シミュレーション方法。
〔６〕格子形状が直交適合格子である〔１〕〜〔５〕のいずれかに記載のゴルフボールの弾道シミュレーション方法。
〔７〕格子形状が非構造格子である〔１〕〜〔５〕のいずれかに記載のゴルフボールの弾道シミュレーション方法。
〔８〕算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離を算出し、打出しから静止するまでの飛距離を算出する〔１〕〜〔７〕のいずれかに記載のゴルフボールの弾道シミュレーション方法。
【０００９】
この弾道シミュレーション方法によれば、実物モデルによる試験評価を実施することなく、表面に複数のディンプルが形成されたゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときのゴルフボールの弾道を推定することができ、更にフォロー、アゲンスト、横風などの風の状態に応じた弾道を推定することも可能である。
【００１０】
この場合、本発明では、上記のように、飛行中のゴルフボールモデル周囲のレイノルズ数（Ｒｅ）を算出すると共に、Ｒｅ＝２５００００〜１７００００を高速領域、Ｒｅ＝１７００００〜１０００００を中速領域、及びＲｅ＝１０００００以下を低速領域とする基準に従って、ゴルフボールモデルの全飛行行程を２つ又は３つの領域に分け、それぞれの領域毎にゴルフボール表面と接する上記格子の高さ（格子高さＡ〜Ｃ）を設定することにより、十分な精度を担保しつつ計算時間の短縮を図るものである。
【００１１】
例えば、ボールが高速で打出されるドライバーショットなどでは、その飛行行程を、ボールが比較的高速で飛行する前域を高速領域（Ｒｅ＝２５００００〜１７００００）、中速で飛行する中盤を中速領域（Ｒｅ＝１７００００〜１０００００）、比較的低速で飛行する後域を低速領域（Ｒｅ＝１０００００以下）の３つの領域に分け、高速領域における上記格子高さ（格子高さＡ）、中速領域における上記格子高さ（格子高さＢ）及び低速領域における上記格子高さ（格子高さＣ）を格子高さＡ＜格子高さＢ＜格子高さＣとなるように設定することにより、計算時間の短縮を図るものである。つまり、レイノルズ数（Ｒｅ）の大きな高速領域ではボール表面近傍を流れる境界層が薄く、レイノルズ数（Ｒｅ）が小さくなる低速領域に移行するに従ってこの境界層は厚くなる。一方、本発明は上記のように、ボールモデル周囲に生成させた格子区画毎に気体流の速度、方向及び圧力を算出するものであるが、本発明者は、上記境界層が厚くなる上記レイノルズ数（Ｒｅ）の小さな領域では、ゴルフボールモデル周囲の格子高さを高くして格子区画の大きさを比較的大きくすることにより格子区画の数を少なくしても実用上問題のない精度を担保し得ることを見出し、更に検討を重ね、レイノルズ数（Ｒｅ）によって全飛行行程を上記高速領域、中速領域、低速領域の３つの領域に分けてＲｅ値が小さい領域ほどゴルフボールモデル表面と接する格子の高さを大きく設定して、精度を損なわない範囲で格子区画の数を少なくすることにより、良好な精度を維持しつつ計算時間を効果的に短縮し得ることを見出し、本発明を完成したものである。なお、飛行行程初期のＲｅ値が１７００００に満たないアイアンショットなどを想定する場合には、飛行行程は上記基準に従って中速領域（Ｒｅ＝１７００００〜１０００００）と低速領域（Ｒｅ＝１０００００以下）の２つの領域に分けられ、格子高さＢ＜格子高さＣとなるようにそれぞれの領域における格子高さが設定される。
【発明の効果】
【００１２】
このように、本発明によれば、実物モデルによる試験評価を実施することなく、ゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときのゴルフボールの弾道を推定することができる。しかも、その推定（シミュレーション）を行う際に生成させる格子区画の数を調整して効率よく計算を行なうことにより全計算時間を短縮することができ、その結果、ゴルフボールの表面形状（ディンプルの形状、配置、大きさなど）の評価の時間が短縮され、また評価の正確性、客観性も十分に担保されることから、ゴルフボールの種類やグレードに合わせてより高性能の製品を効率よく開発することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の方法おけるゴルフボールモデルと気流仮想空間を説明するための説明図であり、（Ａ）は気流仮想空間全体を示す概略図、（Ｂ）は（Ａ）のゴルフボールモデル近傍を示す拡大概略図である。
【図２】本発明におけるゴルフボールモデルの一例を示すもので、（Ａ）は表面に形成された面区画、（Ｂ）は表面に形成されたディンプルを示す図である。
【図３】本発明の弾道シミュレーション方法により推定したゴルフボールの弾道の一例を示す図である。
【図４】回転しながら飛行するゴルフボールに働く力を説明する説明図である。
【発明を実施するための形態】
【００１４】
以下、図面を参照して、本発明につき更に詳しく説明する。
本発明の弾道シミュレーション方法は、表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する方法である。
【００１５】
本発明の弾道シミュレーション方法では、まず、コンピュータにより（Ａ）ボールの打ち出される気流仮想空間（フィールド）を設定し、図１（Ａ）に示したように、その気流仮想空間内（フィールド）に格子を生成する。なお、特に制限されるものではないが、この気流仮想空間内（フィールド）に形成される格子の形状としては、直交適合格子や非構造格子とすることができる。
【００１６】
この場合、この気流仮想空間（フィールド）は、ゴルフボールが打出されてから飛行し、少なくとも着地するまでの領域全体を表すものであり、後述するゴルフボールモデルはこの気流仮想空間内（フィールド）を移動する。また、この気流仮想空間（フィールド）の下面は、開放された面として設定してもよいが、地面を想定した壁面を設定して閉鎖面とすることもできる。これにより、気流仮想空間（フィールド）の下面から空気が出て行く要素が計算から除かれ、より精度の高いシミュレーションを行うことができる。更に、実際にゴルフをプレーする場合には、ボールをティーに載せて打撃することも行われるが、この地面を想定した壁面を設けることにより、初期条件に打撃高さとしてこの壁面からボール最下点までの距離を設定することができ、より実際の使用時に近い条件でより精度の高いシミュレーションを行うことができるものである。
【００１７】
次に、図１に示されているように、（Ｂ）この気流仮想空間内（フィールド）に、表面に複数のディンプルが形成された略球状のゴルフボールモデルを設定する。この場合、このゴルフボールモデルは、３ＤＣＡＤなどにより作成することができ、例えば図２に示したゴルフボールモデルを例示することができる。
【００１８】
更に、このゴルフボールモデルを気流仮想空間内（フィールド）に設定する際、（Ｃ）ゴルフボールモデルの近傍は、格子を仮想空間内（フィールド）の格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画する。このようにゴルフボールモデル近傍の格子を区画することにより、余分な計算を省くことができ、より効率的に計算を行うことができる。
【００１９】
具体的には、まず、図２（Ａ）に示したように、ゴルフボールモデルの表面を例えば一辺を約０．００２ｍｍに区画して、三角形、四角形等の多辺形、又は略三角形、略四角形等の略多辺形の面区画を多数設定し、この各々の面区画を一面とするゴルフボールモデル表面に隣接した格子区画を設定する。ゴルフボールモデル表面に隣接する格子区画は、略四角柱形状等の略多角柱形状や略多角錐形状に設定される。そして、図１（Ｂ）に示したように、このゴルフボールモデル表面に隣接する格子区画からゴルフボールから離間する方向に格子区画の体積が漸次増大するように、気流仮想空間（フィールド）の残部を格子状に区画して、気流仮想領域全域を格子区画で区画する。
【００２０】
このゴルフボールモデル表面に隣接する格子区画以外の気流仮想空間（フィールド）の残部に形成される格子区画の形状としては、六面体、三角柱状五面体、四角錐状五面体、三角錐状四面体などの多面体が挙げられ、これらを適宜組み合わせて設定することができる。
【００２１】
ゴルフボール周囲の気体流は、ゴルフボールに近いものの方がより大きく影響するため、図１（Ｂ）に示されるように、格子区画をゴルフボールモデル近傍は細かく、気体流が与える影響が小さいゴルフボールモデルから遠方は粗く設定するのは上述の通りであるが、この場合ゴルフボールモデルのボール表面から離間する方向への格子区画の体積の増大は、連続的なものであっても、段階的なものであってもよい。
【００２２】
ここで、本発明では、飛行中のゴルフホールモデル周囲のレイノルズ数（Ｒｅ）を算出して、そのＲｅ値に応じて飛行行程を３つの領域に分割し、それぞれの領域毎にゴルフボールモデル表面に接する格子の高さを設定する。この飛行行程の領域区分及び格子高さの設定については、後で詳述する。
【００２３】
次に、上記ゴルフボールモデルに（Ｄ）重量を設定し、（Ｅ）初期条件（初期速度、打出し角度、回転数）を与え、上記気流仮想空間（フィールド）でゴルフボールモデルを運動させる。
【００２４】
ここで、上記のように、気流仮想空間（フィールド）の下面に地面を想定した壁面を設定して、該気流仮想空間（フィールド）の下面を閉鎖面とした場合には、この下面からゴルフボール最下点までの距離を、打撃高さとして初期条件の１つに加えることもできる。これにより、ボールをティーに載せて打撃した場合を想定したシミュレーションをより高精度に行うことができる。
【００２５】
このようにして、（Ｆ）上記（Ｄ）の重量を有する上記ゴルフボールモデルが、上記気流仮想空間内（フィールド）で上記（Ｅ）の初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内（フィールド）の格子に戻す操作を連続的に繰り返す。
【００２６】
そして、ゴルフボールモデルが所定回転数で回転しながら、気流仮想空間内（フィールド）を飛行する際、この回転しているゴルフボールモデルに気流が当たることにより生じるゴルフボールモデル表面に作用する力によって発生する気体流の気流仮想空間内（フィールド）の運動要素の解析を格子区画毎に行う。
【００２７】
ゴルフボールモデルが飛行することにより生じる気流がゴルフボールモデルに当たったときに生じる運動要素は、三次元空間座標系の各軸方向の気体流の速度、気体流の方向及びボールモデル表面に対する気体流の圧力であり、計算に当たって用いる基礎方程式、即ち、以下に示される質量保存則に相当する連続の式（１）〜（３）及び物体の運動保存則に相当するナビエストークの方程式（４）〜（６）に数値を代入して算出することができる。
【００２８】
【数１】
2011224344.tif 000003
【００２９】
上記シミュレーションでは気流仮想空間の格子区画毎に空気の流れを演算により解析することができる。この演算に上記式（１）〜（６）を用い、気流仮想空間を格子区画に区切ったことに対応して上記式（１）〜（６）を離散化して演算を行うことができる。シミュレーションの方法は、シミュレーションの条件等を考慮して有限差分法、有限体積法、境界要素法、有限要素法等を適宜選択して行うことができる。
【００３０】
ここで、場合によっては、（Ｄ−２）上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定し、風の状態に応じたシミュレーションを行うこともできる。
【００３１】
即ち、気流仮想空間で横から、進行方向、又は逆方向の気流を作成し、フォローの風のとき、アゲンストの風のとき、横風のときの状態を考慮してシミュレーションすることもできる。具体的には、例えばアゲンストの風を考慮する場合は、気流仮想空間にゴルフボールモデル前方から所定速度の気流が流入する状態を設定し、気流仮想空間全体の各格子に上記気流が設定された状態で、ゴルフボールモデルを運動させる。アゲンストの風のとき以外の風の状態を考慮する場合は、上記流入する気流の方向を変える事でシミュレーションを行うことができる。
【００３２】
次に、（Ｇ）気流仮想空間内（フィールド）の各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内（フィールド）における飛行中のゴルフボールの揚力係数及び抗力係数を算出する。
【００３３】
この場合、気体流を散乱モデルとして乱流度を加味した速度値を代入して揚力係数ＣＬ及び抗力係数ＣＤを算出することもできる。
【００３４】
そして、上記揚力係数と抗力係数、上記（Ｄ）で設定したボール重量、及び（Ｅ）の初期条件（初期速度、打出し角度、回転数、打撃高さ）からボールモデルを運動させて、（Ｈ）飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定する。
【００３５】
クラブにより打撃されて飛行するゴルフボールは、図４に示したように、重力Ｍｇ、空気による抵抗（抗力）Ｄ、更にボールがスピンを有するためにマグナス効果による揚力Ｌを受けることが知らされている。なお、図中、Ｖは飛行方向、ｃはボール中心を示し、この場合、ゴルフボールｂは矢印Ｒ方向に回転している。
【００３６】
この場合、ゴルフボールに働く力は下記弾道方程式（７）で表される。
Ｆ＝ＦＬ＋ＦＤ＋Ｍｇ  ・・・・（７）
Ｆ：ゴルフボールに働く力
ＦＬ：揚力
ＦＤ：抗力
Ｍｇ：重力
【００３７】
そして、以下の運動方程式により、ゴルフボールの運動を演算することができ、
Ｆcd＝０．５×ＣＤ×ρ×Ａ×Ｖ2    ・・・・（８）
Ｆcl＝０．５×ＣＬ×ρ×Ａ×Ｖ2    ・・・・（９）
Ｆcd＋Ｆcl＋ｍｇ＝ｍ×ｄＶ／ｄｔ   ・・・・（１０）
ｍ：ボール重量
ｇ：重力加速度
ｔ：時間
ＣＬ：揚力係数
ＣＤ：抗力係数
ρ：空気密度
Ａ：ゴルフボールの断面積
Ｖ：ゴルフボールの対空気速度
これから、微少時間飛行した後のゴルフボールの速度が算出され、更に、その結果を放物運動に関するニュートン（Ｎｅｗｔｏｎ）の運動方程式に適用することにより、微少時間飛行した後のボールの飛距離及び高さ変化を算出することができる。上記微分方程式の演算は、オイラー法やルンゲクッタ法を使用することができ、これにより微少時間に対する速度の算出が可能である。
【００３８】
一方、飛行するゴルフボールの回転数（回転速度）は時間の経過によって次第に減衰するが、微少時間経過後の回転数は次の式によって計算することができる。
ω＝ω0×Ｅxp｛−（ＳＲＤ1＋ＳＲＤ2×Ｖ）×ｔ×β｝  ・・・・（１１）
β＝（π×ρ×ｒ4）／ｌ    ・・・・（１２）
ω0：ゴルフボールの初期回転角速度
ω：ゴルフボールの回転角速度
Ｖ：ゴルフボールの対空気速度
ｔ：時間
ρ：空気密度
ｒ：ゴルフボールの半径
ｌ：ゴルフボールの慣性モーメント
ＳＲＤ1，ＳＲＤ2：実験により求められたディンプル固有の減衰係数
【００３９】
この結果、微少時間飛行した後のボールの飛距離、高さ変化が算出され、これを第１ステップとし、更に、算出された回転数及び速度を上記微少時間経過後の回転数及び速度として、上記操作を繰り返すことにより、ボール打出時から落下時までの各々の微少時間経過毎のボール飛距離及び高さを順次算出して弾道を推定することができる。
【００４０】
なお、設定する微少時間は、シミュレーションの効率と精度の点から０．００１〜０．５秒とすることが好ましい。
【００４１】
ここで、本発明では、上記操作により微少時間経過毎のボール飛距離及び高さを算出する際、上記ゴルフボールモデル周囲のレイノルズ数（Ｒｅ）に応じて区分した高速領域、中速領域及び低速領域の領域毎にゴルフボールモデル表面と接する格子の高さを設定する。
【００４２】
まず、飛行する上記ゴルフボールモデル周囲のレイノルズ数（Ｒｅ）を算出すると共に、そのレイノルズ数（Ｒｅ）が２５００００〜１７００００の飛行行程を高速領域、Ｒｅが１７００００〜１０００００の飛行行程を中速領域、Ｒｅが１０００００以下の飛行行程を低速領域とする基準に従って全飛行行程を２つ又は３つの領域に分ける。つまり、上記飛行中のゴルフボールモデル周囲のレイノルズ数（Ｒｅ）を計算して、打出し時からレイノルズ数（Ｒｅ）が１７００００に低下するまでの飛行行程を高速領域とし、Ｒｅ＝１７００００からＲｅ＝１０００００までの飛行行程を中速領域とし、Ｒｅ＝１０００００以下に低下してから着地までの飛行行程を低速領域とする。この場合、上記高速領域と中速領域との境界となるＲｅ＝１７００００は、高速領域内としても中速領域内としてもよいが、時間短縮を優先する場合には中速領域とすることが好ましく、精度を優先する場合は、高速領域とすることが好ましい。同様に、中速領域と低速領域との境界となるＲｅ＝１０００００も、中速領域内としても低速領域内としてもよいが、時間短縮を優先する場合には低速領域とすることが好ましく、精度を優先する場合は中速領域とすることが好ましい。なお、アイアンショットなどの打出し直後のレイノルズ数（Ｒｅ）が１７００００に満たない場合を想定するときには、上記の基準に従って全飛行行程は中速領域と低速領域の２つに分けられる。
【００４３】
上記レイノルズ数の算出は、下記式（１３）により算出される。
Ｒｅ＝ＵＬ／ｖ  ・・・（１３）
Ｒｅ：レイノルズ数
Ｕ：特性速度（ゴルフボールの対空気速度［m/s］）
Ｌ：特性長さ（ゴルフボールの直径［m］）
ｖ：動粘度［m2/s］
【００４４】
そして、本発明では、上記レイノルズ数（Ｒｅ）で区分した高速領域におけるゴルフボールモデル表面と接する格子の高さ（格子高さＡ）、中速領域におけるゴルフボールモデル表面と接する格子の高さ（格子高さＢ）及び低速領域におけるゴルフボールモデル表面と接する格子の高さ（格子高さＣ）を、格子高さＡ＜格子高さＢ＜格子高さＣとなるように設定する。勿論、打出し直後のレイノルズ数がＲｅ＝１７００００に達せず、飛行行程を中速領域と低速領域の２つの領域に分ける場合にも、同様に格子高さＢ＜格子高さＣと設定する。これにより、ゴルフボールモデル周囲の格子を比較的大きくとって格子区画の数を少なくしても実用上問題のない精度を担保し得る領域ほど格子区画の数を少なくすることになり、精度を損なわない範囲で計算回数を少なくし、計算時間の短縮を図ることができるものである。
【００４５】
ここで、上記格子高さＡ〜Ｃのそれぞれの高さは、想定するゴルフボールモデルの初速度や打出し角度、初期回転数などに応じて十分な有用性を有する精度と優位な時間短縮が得られるように適宜設定すればよく、特に制限されるものではないが、上記高速領域における上記格子高さＡは０．００１〜０．０１５ｍｍの範囲、中速領域における上記格子高さＢは０．００５〜０．０２ｍｍの範囲、低速領域における上記格子高さＣは０．０１〜０．０２５ｍｍの範囲でそれぞれ設定することが好ましい。この場合、上記格子高さ範囲の上限を超えると十分な精度が得られずに良好な評価を行うことができなくなる場合があり、一方下限を下回ると十分な時間短縮を図ることができない場合がある。
【００４６】
本発明の弾道シミュレーションを実施する際のより具体的手順を示せば、例えば次の通りである。まず、図１（Ａ）のように、ボールが打ち出される気流仮想空間をコンピュータに設定すると共に、この気流仮想空間内（フィルード）に格子を生成する。この場合、必要に応じて気流仮想空間の下面に地面を想定した壁面を設定する。一方、例えば図２に示したゴルフボールモデルにつき、ボール重量、直径、慣性モーメントを設定すると共に、打出し初期条件として、初速度、打出し角度、初期回転数を設定し、更に必要に応じて打撃高さを設定する。このゴルフボールモデルを図１（Ａ）のように格子を生成させた気流仮想空間内（フィールド）内に設定し、かつ図１（Ｂ）に示した細かい格子をゴルフボールモデル近傍に生成させると共にボールモデルから離れる方向に格子のサイズが漸次増大するように区画する。この場合、ゴルフボールモデル表面に接する格子の高さは、後述する高速領域に設定された格子高さＡとされる。そして、このボールモデルを上記初期条件で気流仮想空間内（フィールド）内を運動させ、揚力係数ＣＬ及び抗力係数ＣＤを計算し、次いで上記式（８），（９）及び（１０）を用いて所定の微少時間飛行した後のゴルフボールの速度を算出し、更にその結果から、放物運動に関するニュートン（Ｎｅｗｔｏｎ）の運動方程式を用いて、上記微少時間飛行した後のボールの飛距離及び高さを算出し、一方で式（１１）及び（１２）を用いて所定時間ｔ経過後の回転数を算出する。また、同時に上記（１３）を用いてボールモデル周囲のレイノルズ数（Ｒｅ）を算出する。
【００４７】
次に、この算出された微少時間飛行した後のゴルフボールの速度及び回転数から、更に微少時間飛行した後のボールの飛距離、左右のばらつき、高さ、回転数を上記した方法で同様に算出し、これを繰り返す。そして、上記レイノルズ数（Ｒｅ）が１７００００に低下した飛行点からボールが上記中速領域に達したと判断し、上記ゴルフボールモデル表面と接する格子の高さを中速領域に設定された所定の格子の高さ（格子高さＢ）へと高くする。更に、この格子高さＢで上記計算を繰り返し、上記レイノルズ数（Ｒｅ）が１０００００に低下した飛行点からボールが上記低速領域に達したと判断し、上記ゴルフボールモデル表面と接する格子の高さを低速領域に設定された所定の格子の高さ（格子高さＣ）へと更に高くし、ゴルフボールが落下する位置まで、この格子高さＣを用いて同様に上記計算を繰り返し、図３に示したゴルフボールの弾道を得ることができる。なお、上述したように、初期のレイノルズ数（Ｒｅ）が１７００００に満たない場合には、飛行行程を中速領域と低速領域との２領域に分けて同様に計算し、ゴルフボールの弾道を得る。
【００４８】
また、必要に応じて、算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離（ラン、転がり距離）を算出し、打出しから静止するまでの飛距離を算出することもできる。
【００４９】
ラン（転がり距離）は、下記式（１４）により計算することができる。
ラン（転がり距離）＝Ａ×Land X velocity −Ｂ×Land Y velocity ＋Ｃ×Land Angle−Ｄ×Land Spin      ・・・（１４）
Ｌand  X velocity= 水平落下速度成分（ｍ／ｓ）
Ｌand  Y velocity= Sin 2θ × 垂直落下速度成分（ｍ／ｓ）
Ｌand  Angle= Cosθ
Ｌand  Spin = 落下スピン量
Ａ，Ｂ，Ｃ，Ｄ：実際のゴルフボールの試験から逆算して得られる固有の係数である。これらの係数を変化させることで、地面の状態の違い（例えば、芝とコンクリート）もシミュレーションに反映させることができる。
【００５０】
本発明の弾道シミュレーション方法では、風洞試験のような実物モデルによる試験評価を実施することなく、ゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときのゴルフボールの弾道を推定することができ、更にフォロー、アゲンスト、横風などの風の状態に応じた弾道を推定することも可能である。しかも本発明によれば、レイノルズ数（Ｒｅ）によってボールの全飛行行程を２つ又は３つの領域に分け、それぞれの領域毎にゴルフボールモデル表面に接する格子の高さを設定して格子区画の大きさ及び数を調整することにより、十分な精度を担保しつつ計算時間の短縮を図ることができるものである。
【実施例】
【００５１】
以下、実施例及び比較例を示し、本発明を具体的に説明するが、本発明は下記実施例に制限されるものではない。
【００５２】
〔実施例、比較例〕
図２に示したゴルフボールモデルにつき、ボール重量４５．３ｇ、直径４２．７ｍｍ、慣性モーメント７．８×１０-6ｋｇ・ｍ2と設定すると共に、打出し初期条件として、初速度６５ｍ／ｓ、打出し角度１０．５°、初期回転数２８００ｒｐｍ、打撃高さ５０ｍｍを設定する。このゴルフボールモデルを図１（Ａ）のように格子を生成させた気流仮想空間内（フィールド）内に設定し、かつ図１（Ｂ）に示した細かい格子をゴルフボールモデル近傍に生成させると共にボールモデルから離れる方向に格子のサイズが漸次増大するように区画する。そして、このボールモデルを上記初期条件で気流仮想空間内（フィールド）内で運動させ、揚力係数ＣＬ及び抗力係数ＣＤを計算し、次いで上記式（８），（９）及び（１０）を用いて所定の微少時間（０．００１秒）飛行した後のゴルフボールの速度を算出し、更にその結果から、放物運動に関するニュートン（Ｎｅｗｔｏｎ）の運動方程式を用いて、上記微少時間（０．００１秒）飛行した後のボールの飛距離及び高さを算出し、また一方で式（１１）及び（１２）を用いて上記微少時間（０．００１秒）経過後の回転数を算出する。
【００５３】
次に、この算出された微少時間（０．００１秒）飛行した後のゴルフボールの速度及び回転数から、更に微少時間（０．００１秒）飛行した後のボールの飛距離、左右のばらつき、高さ、回転数を上記した方法で同様に算出し、その後も同様にして、ゴルフボールが落下する位置まで、上記計算を繰り返し、ゴルフボールの弾道をシミュレートして、飛距離（キャリー）を推定した。
【００５４】
このとき、実施例１及び実施例２のシミュレーションでは、飛行中のゴルフボール周囲のレイノルズ数（Ｒｅ）を上記式（１３）を用いて算出し、全飛行行程をＲｅ＝２５００００〜１７００００の高速領域、Ｒｅ＝１７００００〜１０００００の中速領域、及びＲｅ＝１０００００以下の低速領域の３つの領域に分け、各領域におけるゴルフボールモデル表面と接する格子の高さを表１の通りに設定して、上記シミュレーションを行った。一方、比較例にあっては領域の区分を行なわずに全飛行行程につき格子高さを０．００１ｍｍに設定してシミュレーションを行なった。推定された飛距離とシミュレートに要した時間（計算時間）を表１に示す。また、計算時間については、下記基準で時間短縮の効果を評価した。シミュレーションに用いたコンピュータの詳細は以下の通りである。なお、ボールモデル表面に接する格子区画は直交格子（立方体）とし、そのサイズは高さ（表１記載）×横０．５ｍｍ×厚さ０．５ｍｍとした。
【００５５】
〔計算時間についての評価基準〕
◎・・・１５０時間以内
○・・・２５０時間以内
×・・・２５０時間よりも長い
（なお、下記コンピュータを用いた今回のシミュレーションでは、このように評価基準を設定したが、計算時間の評価基準はコンピュータの性能により適宜変更する必要がある。）
〔コンピュータの詳細〕
ＯＳ：Windows（登録商標） XP
ＣＰＵ：intel Xeon X5482  3.20GHz/L2＝6MB×2  FSB＝1600MHz  2CPU/8Core
メモリ：64GB（4096MB×16）DDR2-800  ECC  REG  FB-DIMM
【００５６】
一方、上記シミュレーションで設定したゴルフボールモデルと同様のゴルフボールにつき、下記の方法で打撃試験を行なって飛距離を測定し、得られた実測値（２１３．２ヤード）と上記シミュレーションによる結果とを比較して、下記基準により精度の評価を行った。結果を表１に示す。
【００５７】
〔打撃試験〕
打撃ロボットにドライバー（Ｗ＃１）［ブリヂストンスポーツ社製「ＴｏｕｒＳｔａｇｅ  Ｘ−Ｄｒｉｖｅ  ７０１」（ロフト９°）］を装着し、実際に打撃して飛距離を測定して実測値とした。このとき、打撃条件は、ボール初速が６５ｍ／ｓ、打出し角度が１０．５°、初期バックスピン量が２８００ｒｐｍとなるように設定した。また、ボールをティーに載せて打撃を行い、地面からボール最下部までの高さを５０ｍｍとした。結果は、飛距離（キャリー）２１３．２ヤード、トータル距離（キャリー＋ラン）２２７．５ヤードであった。
〔精度の評価基準〕
◎・・・実測値の±５ヤード以内
○・・・実測値の±１０ヤード以内
×・・・実測値の±１０ヤードより大
【００５８】
【表１】
2011224344.tif 000004
【００５９】
表１の通り、全飛行行程を所定のレイノルズ数（Ｒｅ）により３つの領域に分け、Ｒｅ＝２５００００〜１７００００の高速領域におけるゴルフボールモデル表面と接する格子の高さ（格子高さＡ）、Ｒｅ＝１７００００〜１０００００の中速領域におけるゴルフボールモデル表面と接する格子の高さ（格子高さＢ）、及びＲｅ＝１０００００以下の低速領域におけるゴルフボールモデル表面と接する格子の高さ（格子高さＣ）を格子高さＡ＜格子高さＢ＜格子高さＣとなるように設定することにより、十分な精度を維持しつつ計算時間を効果的に短縮することができた。
【符号の説明】
【００６０】
Ｍｇ  重量
Ｄ  空気による抵抗（抗力）
Ｌ  揚力
Ｖ  飛行速度
Ｃ  ボールの中心
ｂ  ボール
Ｒ  回転方向
【図３】
2011224344.tif 000005
【図４】
2011224344.tif 000006
【図１】
2011224344.tif 000007
【図２】
2011224344.tif 000008