"表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する弾道シミュレーション方法であって、Aボールの打出される気流仮想空間内フィールドに、格子を生成し、B表面に複数のディンプルが形成された略球状のゴルフボールモデルを気流仮想空間内フィールドに設定し、Cその際、ゴルフボールモデルの近傍は、格子を仮想空間内フィールドの格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画し、Dゴルフボールモデルに重量を設定し、Eゴルフボールモデルに初期条件初期速度、打出し角度、回転数を与え、Fゴルフボールモデルが初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内フィールドの格子に戻し、G気流仮想空間内フィールドの各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内フィールドにおける飛行中のゴルフボールの揚力係数及び抗力係数を算出し、H飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定することを特徴とするゴルフボールの弾道シミュレーション方法。D-2上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定する請求項1記載のゴルフボールの弾道シミュレーション方法。格子形状が直交適合格子である請求項1又は2記載のゴルフボールの弾道シミュレーション方法。格子形状が非構造格子である請求項1又は2記載のゴルフボールの弾道シミュレーション方法。算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離を算出し、打出しから静止するまでの飛距離を算出する請求項1~4のいずれか1項に記載のゴルフボールの弾道シミュレーション方法。本発明は、このような問題に鑑みなされたもので、表面に複数のディンプルが形成されたゴルフボールの弾道を実物モデルによる試験評価によらずに評価して、効率的にゴルフボールを開発する方法を提供することを目的とする。本発明は、上記目的を達成するため、下記1~5のゴルフボールの弾道シミュレーション方法を提供する。1表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する弾道シミュレーション方法であって、Aボールの打出される気流仮想空間内フィールドに、格子を生成し、B表面に複数のディンプルが形成された略球状のゴルフボールモデルを気流仮想空間内フィールドに設定し、Cその際、ゴルフボールモデルの近傍は、格子を仮想空間内フィールドの格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画し、Dゴルフボールモデルに重量を設定し、Eゴルフボールモデルに初期条件初期速度、打出し角度、回転数を与え、Fゴルフボールモデルが初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内フィールドの格子に戻し、G気流仮想空間内フィールドの各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内フィールドにおける飛行中のゴルフボールの揚力係数及び抗力係数を算出し、H飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定することを特徴とするゴルフボールの弾道シミュレーション方法。2D-2上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定する1に記載のゴルフボールの弾道シミュレーション方法。3格子形状が直交適合格子である1又は2に記載のゴルフボールの弾道シミュレーション方法。4格子形状が非構造格子である1又は2に記載のゴルフボールの弾道シミュレーション方法。5算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離を算出し、打出しから静止するまでの飛距離を算出する1~4のいずれかに記載のゴルフボールの弾道シミュレーション方法。この弾道シミュレーション方法によれば、実物モデルによる試験評価を実施することなく、表面に複数のディンプルが形成されたゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときのゴルフボールの弾道を推定することができ、更にフォロー、アゲンスト、横風などの風の状態に応じた弾道を推定することも可能である。その結果、ゴルフボールの表面形状ディンプルの形状、配置、大きさなどの評価の時間が短縮され、また、評価の正確性、客観性も向上することから、ゴルフボールの種類やグレードに合わせてより高性能の製品を効率よく開発することができる。以下、図面を参照して、本発明につき更に詳しく説明する。本発明の弾道シミュレーション方法は、表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する方法である。本発明の弾道シミュレーション方法では、まず、コンピュータによりAボールの打ち出される気流仮想空間フィールドを設定し、図1Aに示したように、その気流仮想空間内フィールドに格子を生成する。なお、特に制限されるものではないが、この気流仮想空間内フィールドに形成される格子の形状としては、直交適合格子や非構造格子とすることができる。この場合、この気流仮想空間フィールドは、ゴルフボールが打出されてから飛行し、少なくとも着地するまでの領域全体を表すものであり、後述するゴルフボールモデルはこの気流仮想空間内フィールドを移動する。次に、図1に示されているように、Bこの気流仮想空間内フィールドに、表面に複数のディンプルが形成された略球状のゴルフボールモデルを設定する。この場合、このゴルフボールモデルは、3DCADなどにより作成することができ、例えば図2に示したゴルフボールモデルを例示することができる。更に、このゴルフボールモデルを気流仮想空間内フィールドに設定する際、Cゴルフボールモデルの近傍は、格子を仮想空間内フィールドの格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画する。このようにゴルフボールモデル近傍の格子を区画することにより、余分な計算を省くことができ、より効率的に計算を行うことができる。具体的には、まず、図2Aに示したように、ゴルフボールモデルの表面を例えば一辺を約0.002mmに区画して、三角形、四角形等の多辺形、又は略三角形、略四角形等の略多辺形の面区画を多数設定し、この各々の面区画を一面とするゴルフボールモデル表面に隣接した格子区画を設定する。ゴルフボールモデル表面に隣接する格子区画は、略四角柱形状等の略多角柱形状や略多角錐形状に設定される。そして、図1Bに示したように、このゴルフボールモデル表面に隣接する格子区画からゴルフボールから離間する方向に格子区画の体積が漸次増大するように、気流仮想空間フィールドの残部を格子状に区画して、気流仮想領域全域を格子区画で区画する。このゴルフボールモデル表面に隣接する格子区画以外の気流仮想空間フィールドの残部に形成される格子区画の形状としては、六面体、三角柱状五面体、四角錐状五面体、三角錐状四面体などの多面体が挙げられ、これらを適宜組み合わせて設定することができる。ゴルフボール周囲の気体流は、ゴルフボールに近いものの方がより大きく影響するため、図1Bに示されるように、格子区画をゴルフボールモデル近傍は細かく、気体流が与える影響が小さいゴルフボールモデルから遠方は粗く設定するのは上述の通りであるが、この場合ゴルフボールモデルのボール表面から離間する方向への格子区画の体積の増大は、連続的なものであっても、段階的なものであってもよい。次に、上記ゴルフボールモデルにD重量を設定し、E初期条件初期速度、打出し角度、回転数を与え、上記気流仮想空間フィールドでゴルフボールモデルを運動させる。このようにして、F上記Dの重量を有する上記ゴルフボールモデルが、上記気流仮想空間内フィールドで上記Eの初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内フィールドの格子に戻す操作を連続的に繰り返す。そして、ゴルフボールモデルが所定回転数で回転しながら、気流仮想空間内フィールドを飛行する際、この回転しているゴルフボールモデルに気流が当たることにより生じるゴルフボールモデル表面に作用する力によって発生する気体流の気流仮想空間内フィールドの運動要素の解析を格子区画毎に行う。ゴルフボールモデルが飛行することにより生じる気流がゴルフボールモデルに当たったときに生じる運動要素は、三次元空間座標系の各軸方向の気体流の速度、気体流の方向及びボールモデル表面に対する気体流の圧力であり、計算に当たって用いる基礎方程式、即ち、以下に示される質量保存則に相当する連続の式1~3及び物体の運動保存則に相当するナビエストークの方程式4~6に数値を代入して算出することができる。上記シミュレーションでは気流仮想空間の格子区画毎に空気の流れを演算により解析することができる。この演算に上記式1~6を用い、気流仮想空間を格子区画に区切ったことに対応して上記式1~6を離散化して演算を行うことができる。シミュレーションの方法は、シミュレーションの条件等を考慮して有限差分法、有限体積法、境界要素法、有限要素法等を適宜選択して行うことができる。ここで、場合によっては、D-2上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定し、風の状態に応じたシミュレーションを行うこともできる。即ち、気流仮想空間で横から、進行方向、又は逆方向の気流を作成し、フォローの風のとき、アゲンストの風のとき、横風のときの状態を考慮してシミュレーションすることもできる。具体的には、例えばアゲンストの風を考慮する場合は、気流仮想空間にゴルフボールモデル前方から所定速度の気流が流入する状態を設定し、気流仮想空間全体の各格子に上記気流が設定された状態で、ゴルフボールモデルを運動させる。アゲンストの風のとき以外の風の状態を考慮する場合は、上記流入する気流の方向を変える事でシミュレーションを行うことができる。次に、G気流仮想空間内フィールドの各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内フィールドにおける飛行中のゴルフボールの揚力係数及び抗力係数を算出する。この場合、気体流を散乱モデルとして乱流度を加味した速度値を代入して揚力係数CL及び抗力係数CDを算出することもできる。そして、上記揚力係数と抗力係数、上記Dで設定したボール重量、及びEの初期条件初期速度、打出し角度、回転数からボールモデルを運動させて、H飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定する。クラブにより打撃されて飛行するゴルフボールは、図4に示したように、重力Mg、空気による抵抗抗力D、更にボールがスピンを有するためにマグナス効果による揚力Lを受けることが知らされている。なお、図中、Vは飛行方向、cはボール中心を示し、この場合、ゴルフボールbは矢印R方向に回転している。この場合、ゴルフボールに働く力は下記弾道方程式7で表される。F=FLFDMg7F:ゴルフボールに働く力FL:揚力FD:抗力Mg:重力そして、以下の運動方程式により、ゴルフボールの運動を演算することができ、Fcd=0.5×CD×ρ×A×V28Fcl=0.5×CL×ρ×A×V29FcdFclmg=m×dV/dt10m:ボール重量g:重力加速度t:時間CL:揚力係数CD:抗力係数ρ:空気密度A:ゴルフボールの断面積V:ゴルフボールの対空気速度これから、微少時間飛行した後のゴルフボールの速度が算出され、更に、その結果を放物運動に関するニュートンNewtonの運動方程式に適用することにより、微少時間飛行した後のボールの飛距離及び高さ変化を算出することができる。上記微分方程式の演算は、オイラー法やルンゲクッタ法を使用することができ、これにより微少時間に対する速度の算出が可能である。一方、飛行するゴルフボールの回転数回転速度は時間の経過によって次第に減衰するが、微少時間経過後の回転数は次の式によって計算することができる。ω=ω0×Exp{-SRD1SRD2×V×t×β}11β=π×ρ×r4/l12ω0:ゴルフボールの初期回転角速度ω:ゴルフボールの回転角速度V:ゴルフボールの対空気速度t:時間ρ:空気密度r:ゴルフボールの半径l:ゴルフボールの慣性モーメントSRD1,SRD2:実験により求められたディンプル固有の減衰係数この結果、微少時間飛行した後のボールの飛距離、高さ変化が算出され、これを第1ステップとし、更に、算出された回転数及び速度を上記微少時間経過後の回転数及び速度として、上記操作を繰り返すことにより、ボール打出時から落下時までの各々の微少時間経過毎のボール飛距離及び高さを順次算出して弾道を推定することができる。なお、設定する微少時間は、シミュレーションの効率と精度の点から0.001~0.1秒とすることが好ましい。本発明の弾道シミュレーションのより具体的な実施例を示すと、例えば、図2に示したゴルフボールモデルにつき、ボール重量45.3g、直径42.7mm、慣性モーメント7.8×106kgm2と設定すると共に、打出し初期条件として、初速度67m/s、打出し角度10°、初期回転数2520rpmを設定する。このゴルフボールモデルを図1Aのように格子を生成させた気流仮想空間内フィールド内に設定し、かつ図2A及び図1Bに示した細かい格子をゴルフボールモデル近傍に生成させると共にボールモデルから離れる方向に格子のサイズが漸次増大するように区画する。そして、このボールモデルを上記初期条件で気流仮想空間内フィールド内を運動させ、揚力係数CL及び抗力係数CDを計算し、次いで上記式8,9及び10を用いて所定時間t例えば0.01秒飛行した後のゴルフボールの速度を算出し、更にその結果から、放物運動に関するニュートンNewtonの運動方程式を用いて、所定時間t飛行した後のボールの飛距離及び高さを算出し、また、一方で式11及び12を用いて所定時間t経過後の回転数を算出する。次に、この算出された所定時間t飛行した後のゴルフボールの速度及び回転数から、更に所定時間t飛行した後のボールの飛距離、左右のばらつき、高さ、回転数を上記した方法で同様に算出し、その後も同様にして、ゴルフボールが落下する位置まで、上記計算を繰り返し、図3に示したゴルフボールの弾道を得ることができる。また、必要に応じて、算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離ラン、転がり距離を算出し、打出しから静止するまでの飛距離を算出することもできる。ラン転がり距離は、下記式13により計算することができる。ラン転がり距離=A×LandXvelocity-B×LandYvelocityC×LandAngle-D×LandSpin13LandXvelocity水平落下速度成分m/sLandYvelocitySin2θ×垂直落下速度成分m/sLandAngleCosθLandSpin落下スピン量A,B,C,D:実際のゴルフボールの試験から逆算して得られる固有の係数である。これらの係数を変化させることで、地面の状態の違い例えば、芝とコンクリートもシミュレーションに反映させることができる。"