(57)【特許請求の範囲】
【請求項１】
ゴルフボールモデルを節点からなる多数の要素にメッシュ分割し、ゴルフボールの材料物性を入力し、有限要素法解析により、ゴルフクラブヘッドとゴルフボールの衝突を想定したシミュレーションを行い、上記衝突時にゴルフボールモデルに発生するひずみ量を演算し、
上記ゴルフボールモデルの各要素の応力、ひずみ成分、節点座標値を出力し、各要素の時々刻々の６成分の応力、ひずみの値を算出し、
上記６成分の応力、ひずみの値より、各要素の各成分における応力とひずみの関係を求め、各要素のエネルギーロス値を算出し、
上記エネルギーロス値をゴルフボールモデルに含まれる全ての要素に関して算出し、可視化ソフトに読み込ませ、ゴルフボールモデルのエネルギーロスを可視化していることを特徴とするゴルフボールのエネルギーロス評価方法。
【請求項２】
上記エネルギーロスを可視化する際に、上記各要素の節点の節点座標値からゴルフボールモデルの変形形状を表示している請求項１に記載のゴルフボールのエネルギーロス評価方法。
【請求項３】
上記材料物性を、実際のゴルフクラブヘッドによるゴルフボール打撃時と同等のひずみ及びひずみ速度の条件下で測定している請求項１または請求項２に記載のゴルフボールのエネルギーロス評価方法。
【請求項４】
上記材料物性を、スプリットホプキンソン棒試験機により、実際のゴルフクラブヘッドによるゴルフボール打撃時と同等のひずみ及びひずみ速度の条件下で測定している請求項１または請求項２に記載のゴルフボールのエネルギーロス評価方法。
【請求項５】
上記材料物性として、ヤング率、せん断係数、粘性係数、あるいは、損失係数等の材料の剛性と粘性に関わる材料物性を測定している請求項１乃至請求項４のいずれか１項に記載のゴルフボールのエネルギーロス評価方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ゴルフボールのエネルギーロス評価方法に関し、特に、ゴルフボール内部に生じるエネルギーロスを、実際にボールを試作することなくシミュレーションで可視化できるようにするものである。
【０００２】
【従来の技術】
ゴルファーがゴルフボールに求める重要な性能の１つに、ボールの飛距離が挙げられる。飛距離の大きなゴルフボールはゴルファーに爽快感を与え、また、スコアメイクにも寄与する。この飛距離向上のためには、ゴルフボールの反発性能の向上が必要である。
【０００３】
従って、ゴルフボールの設計において、ボールが打撃されたときの反発性能は、大きな評価項目の一つであり、従来、ゴルフボールの反発性能を向上させるための種々の提案が挙げられている。
【０００４】
その一つとして、ゴルフボールの打撃時にボール内部で生じるエネルギーのロス（損失）に着目することが考えられている。即ち、ゴルフボールは打撃物によりインパクト（衝突）されると大変形が生じた後、打撃物から分離して反発を生じる。このときのゴルフボールの反発係数は衝突中にボール内部で生じるエネルギーのロス（損失）に大きく影響されることが判っている。
【０００５】
また、ゴルフボールを打撃した際、スイートスポットを外すと、衝突のエネルギーがボールの外に逃げてしまい、エネルギーがロスすることが判っている。このように、打撃時にゴルフボール内部で生じるエネルギーロスは、ゴルフボールの飛距離に大きな影響を与えるため、ゴルフボールの反発性能を向上させるために、上記エネルギーロスは有益な情報である。
【０００６】
【発明が解決しようとする課題】
しかしながら、打撃物による衝突（インパクト）中のゴルフボール内部を、実際に観察することは極めて困難であるため、現実には打撃物のインパクトによって生じるゴルフボール内部のエネルギーロスを評価できていないのが実状である。従って、ゴルフボール内部のどの部分（例えば、ボール表面のカバー部分、あるいは、ボールの中心付近等）で、どの程度エネルギーロスが生じているかを予測、観察できていないため、エネルギーロスをボールの反発性能向上に十分に活かせていないという問題がある。
【０００７】
また、従来のゴルフボールの設計では、ボールを構成する材料を決定し、その後、実際にボールの試作を行い、試作したボールを実際に打撃することで、設計したボールの反発性能に関しての評価を行っていた。この設計方法では、ボールの試作において非常に費用と時間を要するという問題がある。
【０００８】
本発明は上記した問題に鑑みてなされたものであり、実際の打撃物による衝突（インパクト）時にゴルフボール内部に生じるエネルギーロスの評価を可能とし、ゴルフボールの反発性能の向上を図ると共に、ゴルフボールの設計を容易にすることを課題としている。
【０００９】
【課題を解決するための手段】
上記課題を解決するため、本発明は、ゴルフボールモデルを節点からなる多数の要素にメッシュ分割し、ゴルフボールの材料物性を入力し、有限要素法解析により、ゴルフクラブヘッドとゴルフボールの衝突を想定したシミュレーションを行い、上記衝突時にゴルフボールモデルに発生するひずみ量を演算し、
上記ゴルフボールモデルの各要素の応力、ひずみ成分、節点座標値を出力し、各要素の時々刻々の６成分の応力、ひずみの値を算出し、
上記６成分の応力、ひずみの値より、各要素の各成分における応力とひずみの関係を求め、各要素のエネルギーロス値を算出し、
上記エネルギーロス値をゴルフボールモデルに含まれる全ての要素に関して算出し、可視化ソフトに読み込ませ、ゴルフボールモデルのエネルギーロスを可視化していることを特徴とするゴルフボールのエネルギーロス評価方法を提供している。
【００１０】
このように、打撃物とゴルフボールの衝突時に、ゴルフボール内部で発生するエネルギーロスの状態を可視化し観察することにより、ゴルフボールの反発性能に大きく影響するエネルギーロスの評価が可能となる。これにより、ゴルフボールの設計上、反発性能の向上に有益な情報を得ることができるため、反発性能に優れたゴルフボールを設計することができる。
【００１１】
上記の算出法で算出した各要素のエネルギーのロスを３次元的に可視化、表示することで、ゴルフボール内部のどの部分で、エネルギーのロス（損失）がどの程度生じているかを予測、観察することが可能となる。なお、エネルギーロスとは、各要素の各成分の変形（ヒステリシス）のロスである。
【００１２】
また、有限要素法を用いてシミュレーションを行っているため、試作の費用と時間の節約が実現でき、様々な材料を用いたあらゆる構造のゴルフボールの設計を短期間に達成することができる。
【００１３】
ゴルフボールモデルを節点からなる多数の要素にメッシュ分割し、ゴルフボールの材料物性を入力し、有限要素法解析により、ゴルフクラブヘッドとゴルフボールの衝突を想定したシミュレーションを行い、上記衝突時にゴルフボールモデルに発生するひずみ量を演算している。
【００１４】
解析モデルにおいて、有限要素法解析をするにあたり、初期条件を設定している。設定条件としては、ゴルフボールモデルの大きさ、形状、構造及び材料物性、及び、衝突する物体（例えばゴルフクラブヘッドモデル）の、速度、形状、構造、材料物性等が挙げられる。
【００１５】
上記ゴルフボールモデルの各要素の応力、ひずみ成分、及び節点座標値を出力し、各要素の時々刻々の６成分の応力、ひずみの値を算出している。
【００１６】
ひずみ、応力は３次元空間では垂直方向３成分、せん断方向３成分の６つの成分から構成される。従って、有限要素法によるシミュレーションにおいて、ゴルフボールを構成する各要素に関して、それぞれ上記６成分でエネルギーロスの算出を行っている。これにより、実際の打撃物とゴルフボールの衝突時とほぼ同じ状態での解析が可能となり、シミュレーションの精度が向上する。
【００１７】
また、ゴルフボールモデルの各要素のひずみ、応力の３次元空間での垂直方向３成分、せん断方向３成分の基準となる座標系は有限要素法での要素座標系としても良い。即ち、各要素の要素座標変換に必要な節点座標値を出力し、各要素の時々刻々の要素座標系での６成分の応力、ひずみの値を算出しても良い。衝突時のゴルフボールの変形に伴い各要素に生じる回転移動によって座標軸が移動するため、要素座標系にすることで座標軸を回転移動に合わせて移動させることができ、回転による影響を除去し変形のみを考慮することができる。
【００１８】
要素座標系でなく、全体座標系で６成分の有限要素法解析を行いエネルギーロスを計算すると、衝突によりゴルフボールモデルの要素が回転の移動を伴う場合に、要素の回転移動の影響が含まれる。回転移動の影響を除去し、各方向のエネルギーロスの成分を算出する場合には、要素座標系の垂直、せん断方向のそれぞれ３成分を算出し、エネルギーロスを算出することが好ましい。
【００１９】
６成分の応力、ひずみの値より、各要素の各成分における応力とひずみの関係を求め、各要素のエネルギーロス値を算出している。
具体的には、各要素の各成分において、時々刻々の応力とひずみの関係から応力−ひずみ曲線を描く。上記応力−ひずみ曲線より、ひずみの絶対値が最大となる点を求め、衝突後ひずみ量が増加し、上記ひずみの絶対値が最大となる点までの仕事量を算出する。同様に、ひずみが除荷の方向、つまりひずみが減少しひずみ０になるまでの仕事量を算出する。上記２つの仕事量よりエネルギーロス値を算出することができる。
上記エネルギーロス値は上記６成分それぞれで存在するため、６成分の合計値がその要素のエネルギーロス値となる。
【００２０】
上記エネルギーロス値をゴルフボールモデルに含まれる全ての要素に関して算出し、可視化ソフトに読み込ませ、ゴルフボールモデルのエネルギーロスを可視化している。
具体的な可視化の実現方法としては、例えば可視化汎用ソフトｅｎｓｉｇｈｔ等が挙げられる。これにより、ボールを構成する各要素で生じるエネルギーロス値を例えば色により区別した図（ｃｏｎｔｏｕｒ図）等で表示することができる。
【００２１】
上記エネルギーロスを可視化する際に、上記各要素の節点の節点座標値からゴルフボールモデルの変形形状を表示している。これにより、ゴルフボールのエネルギーロスと衝突時の変形形状とを併せて可視化することができ、ゴルフボールの設計に有効である。
【００２２】
また、エネルギーロスの可視化において、各要素のエネルギーロスの値の大小により色分けし、コンター図においてエネルギーロスの違いを色の違いで（エネルギーロスの大小により暖色系から寒色系に色を徐々に変化させる）表示することで、可視化を実現すると、一目で、より明確にエネルギーロスの評価を行うことができる。
【００２３】
上記スプリットホプキンソン棒試験機による測定では、試験片に高速でかつ大変形のひずみを与えることができる。従って、実際にゴルフクラブヘッドにより、ゴルフボールが打撃されたときと同等のひずみ、ひずみ速度、応力の条件下で、ゴルフボールを構成する材料物性を測定することができる。このスプリットホプキンソン棒試験機で測定された材料物性を用いることで、有限要素法解析により、実際の打撃時にボールを構成する材料が受けているひずみ、応力の状態を正確にシミュレーションすることができるため、エネルギーロスの可視化の精度を向上させることができる。
【００２４】
実際のゴルフクラブヘッドによる打撃時に、ゴルフボールに発生する最大圧縮ひずみ、最大ひずみ速度を想定して、最大圧縮ひずみが０．０５〜０．５０であり、最大ひずみ速度が５００／ｓ〜１００００／ｓ、好ましくは５００／ｓ〜５０００／ｓであるような条件下で測定された材料物性を用いて、有限要素法により解析することが好ましい。高速大変形時に生じるひずみ、ひずみ速度の条件であるため、評価の精度を向上することができる。
【００２５】
なお、試験片に高速でかつ大変形のひずみを与えることができ、実際にゴルフクラブヘッドにより、ゴルフボールが打撃されたときと同等の最大圧縮ひずみ、最大ひずみ速度の条件下で、材料物性を測定することができれば、上記スプリットホプキンソン棒以外の測定方法により材料物性を求めても良いことは言うまでもない。
【００２６】
有限要素法解析に用いられる材料物性としては、ヤング率、せん断係数、粘性係数、あるいは、損失係数等の材料の剛性と粘性に関わる材料物性が好ましい。上記のような材料物性のいずれか１つ、あるいは、いくつかを組み合わせて有限要素法解析に用いることにより、ゴルフボールのエネルギーロスの評価の精度を高めることができる。
【００２７】
有限要素法では解析対象（ゴルフボールモデル）を多数の要素にメッシュ分割しているが、ゴルフボールモデル全範囲の要素の数は、全体モデルで１０００〜１０００００が好ましく、さらに好ましくは、２５００〜６００００、より好ましくは２５００〜２００００が良い。なお、上限の値は、現段階での計算機の能力を鑑みてのものであり、今後計算機の能力が向上するにつれ、解析の時間が短縮されるため上限の範囲は、変わることが容易に想像できる。
【００２８】
要素の形状は６面体が好ましく、要素のサイズは、解析精度の観点より細かければ細かいほど好ましいが、計算時間は細かいほど、その分増加する。
【００２９】
要素座標系を定義する方法（要素座標系の決定方法）としては、例えば、以下のような方法が挙げられる。要素に含まれる２つの節点より要素座標の第１軸を決定し、さらに、もう１つの節点で平面を定義する。この時、第１軸に垂直で、かつ上記平面の法線方向になる方向を第２軸と定義する。さらに、第１軸及び第２軸に垂直な方向を第３軸と定義する。上記第１軸、第２軸、第３軸から要素座標を決定する。
【００３０】
本発明によりシミュレーションされるゴルフボールは、架橋ゴム層等の単体からなる所謂１ピースゴルフボールであってもよく、架橋ゴム層等のコアにカバーが被覆された所謂２ピースゴルフボールでもよく、また、３層以上から構成されていてる所謂マルチピースゴルフボールであってもよく、あらゆる構造のゴルフボールに適応可能である。
【００３１】
また、本発明によりシミュレーションされるゴルフボールの材料としては、一般にゴルフボールに用いられる材料であればよく、また、各種添加剤が配合されていてもよい。即ち、有限要素法解析に必要な材料物性を測定できる材料であれば、あらゆる材料からなるゴルフボールに適応可能である。
【００３２】
【発明の実施の形態】
以下、本発明の実施形態を図面を参照して説明する。
図１は、本発明のシミュレーション方法のフローチャートを示す。
（ステップ１）ゴルフボールモデルを節点からなる多数の要素にメッシュ分割し、ゴルフボールの材料物性を入力し、有限要素法解析により、ゴルフクラブヘッドとゴルフボールの衝突を想定したシミュレーションを行い、上記衝突時にゴルフボールモデルに発生するひずみ量を演算する。
（ステップ２）上記ゴルフボールモデルの各要素の応力、ひずみ成分、節点座標値を出力し、各要素の時々刻々の６成分の応力、ひずみの値を算出する。
（ステップ３）上記６成分の応力、ひずみの値より、各要素の各成分における応力とひずみの関係を求め、各要素のエネルギーロス値を算出する。
（ステップ４）上記エネルギーロス値をゴルフボールモデルに含まれる全ての要素に関して算出し、可視化ソフトに読み込ませ、ゴルフボールモデルのエネルギーロスを可視化する。また、ゴルフボールモデルの各節点座標値より、変形形状も合わせて可視化する。
【００３３】
これにより、ゴルフボールの反発性能に大きく影響するエネルギーロスの評価が可能となり、ゴルフボールの設計上、反発性能の向上に有益な情報を得ることができる。
【００３４】
ステップ１では、ゴルフボールモデルを節点からなる多数の要素にメッシュ分割し、ゴルフボールの材料物性を入力し、有限要素法解析により、ゴルフクラブヘッドとゴルフボールの衝突を想定したシミュレーションを行い、上記衝突時にゴルフボールモデルに発生するひずみ量を演算する。
【００３５】
本発明のシミュレーション方法により、シミュレーションを行うゴルフボールモデル１０を示す。図２に示すように、ゴルフボールモデル１０は、ハイシスポリブタジエンゴムを主成分とする１ピースボールを想定し、直径４２．８ｍｍの球状としている。
【００３６】
有限要素法解析を行うにあたり、解析モデルにおいて初期条件を設定する。即ち、上記ゴルフボールモデル１０の大きさ、形状、構造及び材料物性等の初期条件を設定すると共に、ゴルフボールモデル１０を多数の要素１１にメッシュ分割し、多数の節点１２を得ている。本実施形態では、３２００個の要素１１に分割し、要素１１の形状は６面体、要素１１のサイズは約１．２５ｍｍ〜５．６２ｍｍ、要素１１の体積は約１．９５ｍｍ３〜６３．１ｍｍ３である。なお、図中の格子状の線は、メッシュを示す線である。
上記設定条件に基づいて、打撃物の衝突時から所定時間後のゴルフボールモデル１０に発生するひずみ量を演算する。
【００３７】
また、上記ハイシスポリブタジエンゴムを主成分とする材料の物性を、後述するスプリットホプキンソン棒試験機により測定する。これにより、実際のゴルフクラブヘッドによる打撃時と同等のひずみ、ひずみ速度の条件下における材料の物性値を得ている。
上記材料の物性値を入力し、以下に示すように、ゴルフボールモデル１０について、有限要素法によりシミュレーションを行っている。
【００３８】
ゴルフクラブヘッド（打撃物）のゴルフボールへの衝突を想定し、図３（Ａ）（Ｂ）（Ｃ）に示すように、シミュレーションを行っている。即ち、上記ゴルフボールモデル１０に、打撃物として２００ｇ（ゴルフクラブヘッドの重さと同一）の円筒形状のアルミニウム製中空棒モデル２０を速度４５ｍ／ｓで衝突させた時のゴルフボールモデル１０の状態を有限要素法により解析し、シミュレーションを行っている。これにより、各要素１１における所定時間におけるひずみ量を演算している。
なお、アルミニウム製中空棒モデル２０の円形面２０ａを衝突面とし、衝突面はフラットであり、ゴルフボールモデル１０とフラットな正面衝突としている。また、アルミニウム製中空棒モデル２０の円形面２０ａの中心点２０ｂが、最初のゴルフボールモデル１０との衝突点となるように衝突させている。
【００３９】
ステップ２では、上記ゴルフボールモデルの各要素の応力、ひずみ成分、及び節点座標値を出力し、各要素の時々刻々の６成分の応力、ひずみの値を算出する。
【００４０】
まず、上記有限要素法解析によるシミュレーション結果から、ゴルフボールモデル１０の各要素１１が、打撃物（アルミニウム製中空棒モデル２０）との衝突時から受けている時々刻々の応力、ひずみ成分と、上記各要素１１の節点１２の座標データ（節点座標値）とを算出する。
【００４１】
次に、上記算出値より、各要素１１の時々刻々の応力、ひずみの値を下記の６成分について算出する。各要素１１について、ひずみ及び応力は、垂直方向３成分、せん断方向３成分の計６成分から構成されるため、上記６成分について算出している。
【００４２】
なお、要素座標変換に必要な節点座標を出力し、各要素の時々刻々の要素座標系での６成分の応力、ひずみの値を算出することもできる。ここで、要素座標系を定義する方法として、具体的に、６面体形状の要素３０における要素座標の決定方法を図４（Ａ）（Ｂ）（Ｃ）に示す。
６面体形状の要素３０を構成する８つの節点ａ〜節点ｆのうち、まず、２つの節点を選ぶ（本具体例では節点ａと節点ｂ）。この２つの節点ａと節点ｂとを通る直線を要素座標の第１軸とする。
さらに、もう１つの節点ｃと第１軸とで平面αを定義する（平面αは面ａｂｃｄとなる）。節点ａを通る平面αの法線を第２軸とする。
次に、第１軸と第２軸の両方に垂直な軸を第３軸とする。このように、第１軸、第２軸、第３軸により要素座標を決定する。
【００４３】
本実施形態においても、上記のような方法で要素座標を決定しても良い。要素座標系とすることで、要素が回転の移動を伴う場合でも、回転移動の影響を除去した各成分のエネルギーロスを算出することができる。
【００４４】
ステップ３では、上記６成分の応力、ひずみの値より、各要素の各成分における応力とひずみの関係を求め、各要素のエネルギーロス値を算出する。
以下、具体的なエネルギーロス算出方法の一例を示す。以下の方法以外にも、応力−ひずみ曲線のループ内面積を算出すれば良い。
【００４５】
まず、ステップ２で算出した各要素の時々刻々の６成分の応力、ひずみの値を用い、上記応力とひずみの関係をグラフに示す。
即ち、上記有限要素法解析によるシミュレーションでは、ゴルフボールモデル１０を構成する各要素１１の各成分において、打撃物（アルミニウム製中空棒モデル２０）が衝突し材料が応力を受け、ひずみが生じた後、材料が元に戻る過程において、図５のようなひずみと応力の関係（応力−ひずみ曲線）を示す。図５において、横軸はひずみ、縦軸は応力を表し、ゴルフボールモデル１０中の１要素の１成分の応力、ひずみの関係を示している。
【００４６】
図５中のひずみと応力の関係を示す曲線において、材料に衝撃が加えられるのに伴い、上記曲線は矢印に示すように時々刻々、点Ｏから点Ａ（ひずみの絶対値が最大値εｍａｘとなる点）へ行き、再び点Ｏに戻る。
その際、点Ｏから点Ａへ行くまでの仕事量と、再び点Ｏに戻るまでの仕事量との差（図中縦縞部分）が生じる。この仕事量の差（縦縞部分の面積）の値により、打撃物のゴルフボールへの衝突時に生じるエネルギーロスを求めることができる。上記仕事量の差は、ゴルフボールモデル１０を構成する各要素１１の各成分で、それぞれ生じている。
【００４７】
上記仕事量の差を求めるため、次に、各要素１１の各成分の応力、ひずみに関して、ひずみの値が０から、その絶対値εｍａｘが最大値に達するまでの仕事量を算出する。
【００４８】
具体的には、図６に示すように、応力−ひずみ曲線が点Ｏから点Ａへ行く（ひずみ量が増加している）過程において、ひずみがε１（応力：σ１）からε２（応力：σ２）にｄεだけ変化したときの仕事量（図６中の斜め縞の面積）ｄＬ１２を下記数式（１）により算出する。
（数式１）
ｄＬ１２＝１／２（σ１＋σ２）・ｄε−−−（１）
上記同様の方法により、点Ｏから点Ａまでｄεずつ変化したときの仕事量をそれぞれ求め、それら仕事量の値を合計し、図７に示すように、点Ｏから点Ａへ行くまでの仕事量の総和Ｌ１（縞線部分に示す仕事量）を算出する。
【００４９】
さらに、各要素１１の各成分の応力、ひずみに関して、点Ａ（ひずみの絶対値の最大値εｍａｘ）からひずみが除荷の方向、つまりひずみが減少し、ひずみ０近傍になるまでの仕事量Ｌ２を算出する。
【００５０】
具体的には、図８に示すように、応力−ひずみ曲線が点Ａから再び点Ｏへ戻る（ひずみ量が減少している）までの過程において、ひずみがε３（応力：σ３）からε４（応力：σ４）にｄεだけ変化したときの仕事量（図８中の斜め縞の面積）ｄＬ３４を下記数式（２）により算出する。
（数式２）
ｄＬ３４＝１／２（σ３＋σ４）・  │  ｄε  │  −−−（２）
上記同様の方法により、点Ａから点Ｏまでｄεずつ変化したときの仕事量をそれぞれ求め、それらの値を合計し、図９に示すように、点Ａから再び点Ｏに戻るまでの仕事量の総和Ｌ２（縞線部分に示す仕事量）を算出する。
【００５１】
上記の方法により算出された点Ｏから点Ａへ行くまでの仕事量の総和Ｌ１と、点Ａから再び点Ｏに戻るまでの仕事量の総和Ｌ２の値から、上記仕事量の差Ｌ（図５に示す縦縞部分の面積）を下記数式（３）により算出する。
（数式３）
Ｌ＝Ｌ１―Ｌ２−−−（３）
【００５２】
上記仕事量の差Ｌの値と、各要素１１の初期形状での体積Ｖとの積により、下記数式（４）に示すように、各要素１１の各成分におけるエネルギーロスを算出できる。
（数式４）
エネルギーロス＝Ｌ・Ｖ−−−（４）
エネルギーロスの値は、６成分それぞれで存在するため、６成分各々において、上記同様の方法でエネルギーロスの値を算出し、６成分のエネルギーロスの合計値がその要素のエネルギーロス値となる。
【００５３】
ステップ４では、上記エネルギーロス値をゴルフボールモデルに含まれる全ての要素に関して単位体積当たりのエネルギーロス値を算出し、可視化ソフトに読み込ませ、ゴルフボールモデルのエネルギーロスを可視化する。
以下、具体的なエネルギーロスの可視化方法を示す。
【００５４】
まず、ステップ３で算出された各要素１１のエネルギーロス値を各要素１１がもつエネルギーロス値として可視化ソフトに読み込ませる。
【００５５】
次に、図１０に示すように、ゴルフボールモデル１０の各要素におけるエネルギーロスをコンター（ｃｏｎｔｏｕｒ）図（断面図）で示し、ゴルフボールモデル１０の内部に生じるエネルギーロスを可視化する。
具体的には、ゴルフボールモデル１０において、エネルギーロスの値により表示領域を定め、ゴルフボール内部の状態を図示している。即ち、最もエネルギーロスが大きい領域を領域（赤）、最もエネルギーロスが小さい領域を領域（青）とし、エネルギーロスの値により順次領域を定め（領域赤＞…黄＞…緑＞…青）、ゴルフボールの断面に表示している。
【００５６】
また、上記解析により得られた各要素１１の節点１２の座標データ（節点座標値）からゴルフボールモデル１０の変形形状も併せて表示し、打撃物の衝突時のゴルフボールモデル１０の形状も可視化する。
【００５７】
以上より、打撃物の衝突時にゴルフボール内部のどの部分で、どの程度エネルギーロスが生じているかを予測、観察すると共に、ゴルフボールの変形形状をも予測することができる。
【００５８】
エネルギーロスの可視化においては、各要素のエネルギーロスの値の大小により色分けし、コンター図においてエネルギーロスの違いを色で表示することで、可視化を実現するのが最適である。例えば、図１０に示すように、各要素のエネルギーをそれぞれ色分けして表示すると、一目でエネルギーロスの大小を評価することができる。上記実施形態ではエネルギーロスの値により色の表示を分けて、ゴルフボール内部の状態を表示しているが、表示の際の領域の数や各領域におけるエネルギーロスの値の範囲は、適宜設定することができる。
【００５９】
上記実施形態で示した方法以外のエネルギーロスの算出方法として、以下のような手法が挙げられる。
ゴルフボールモデルの要素のうち、ある１要素における要素座標系の１方向の応力−ひずみ曲線を図１１に示す。まず、図１１の曲線より、ひずみが最大となるεｍａｘを決定する。次ぎに、図１２に示すように、例えば、ひずみが０．１からｄεだけ進んだときのエネルギーロスを、図示するように作図し図中斜め縞の範囲として算出する。よって、ひずみが０から（０．１＋ｄε）まで進んだときのエネルギーロスを、図１３に示す斜め縞の範囲として算出することができる。このようにして、応力−ひずみ曲線が０近傍に戻る前での途中時点でのエネルギーロスを算出することができる。この手法によれば、変形の途中段階で、その時点までで生じるエネルギーロスを算出することができる。また、ひずみが０に戻りきらず０近傍で終了しループが閉になりきらないこともある。この場合でも、ひずみの最大値の絶対値に対して０％〜２０％、好ましくは０％〜１０％、さらに好ましくは０％〜５％のひずみ値程度まで小さくなり、０近傍に達すれば、ループが閉じなくてもヒステリシスロスを計算し、エネルギーロスとして算出することができる。
【００６０】
本実施形態では、１ピースゴルフボールモデルをシミュレーションしているが、本発明によりシミュレーションされるゴルフボールは、２ピースゴルフボール、あるいは、３層以上から構成されていてるマルチピースゴルフボール等あらゆる構造のゴルフボールにおいて、有限要素法によりシミュレーションし、エネルギーロスを可視化することができる。
【００６１】
また、ゴルフボールに用いる材料においても、材料物性の測定が可能な材料であれば、あらゆる材料を用いた場合のシミュレーションを行うことができる。
【００６２】
従って、あらゆる構造のゴルフボールにおいて、あらゆる材料を用いた場合のエネルギーロスを可視化できるため、試作をすることなく様々なゴルフボールを容易に設計することができる。
【００６３】
なお、本実施形態では、スプリットホプキンソン棒試験機により、材料物性を測定しているが、試験片に高速でかつ大変形のひずみが与えられ、実際にゴルフクラブヘッドにより、ゴルフボールが打撃されたときと同等の最大圧縮ひずみ、最大ひずみ速度の条件下で、材料物性（ヤング率や損失係数等）を測定することができれば、上記スプリットホプキンソン棒以外の測定方法により材料物性（ヤング率や損失係数等）を求めても良い。
【００６４】
（実施例１）
ゴルフボールモデルは、ハイシスポリブタジエンゴムを主成分とする１ピースボール（ハイシスポリブタジエンゴム：ＢＲ０１    ＪＳＲ製    １００ｐｈｒ、アクリル酸亜鉛（ＺＤＡ）３１ｐｈｒ、酸化亜鉛（ＺｎＯ）２０ｐｈｒ、ジクミルパーオキサイド（ＤＣＰ）０．４ｐｈｒからなる混合物を１６０℃で３０分、直径４２．８ｍｍの金型で圧縮成型）を想定し、直径４２．８ｍｍの球状とした。
【００６５】
上記ハイシスポリブタジエンゴムを主成分とする材料の材料物性（ヤング率と損失係数）を、後述するスプリットホプキンソン棒試験機により測定した（室温２３℃、相対湿度５０％、打撃棒の衝突速度２５ｍ／ｓ）。測定により得られたヤング率の値は７７．８ＭＰａ、損失係数の値は０．５４９３であった。また、測定時の最大圧縮ひずみの値は０．２８９７、最大ひずみ速度の値は３９０８．１［／ｓ］であった。
【００６６】
上記実施形態と同様の方法により、上記実施例１のゴルフボールモデル１０’について、図１４に示すように、初速４５ｍ／ｓでアルミニウム製中空棒モデル２０’を衝突させた時のエネルギーロスを算出した。ゴルフボールモデルは、１２７９７個の要素に分割し、要素の形状は６面体、サイズは０．５ｍｍ〜１．６９ｍｍ、体積は約０．１２５ｍｍ３〜２．７３ｍｍ３とした。図１５に示すように、ゴルフボールモデルの要素のうち、１要素におけるの１方向の応力−ひずみ曲線を得た。
上記仕事量の差Ｌの算出式に従い、Ｌの値を計算すると、
Ｌ＝２．２４ＭＰａであった。
同様に、他の５方向に関しても仕事量の差Ｌの算出を行い、６方向の仕事量の差の合計を求めた。６方向の仕事量の差の合計と要素の体積１．１５×１０−９ｍ３を掛けて、１要素のエネルギーロスを算出すると、２．６０２×１０−３Ｎ・ｍであった。
【００６７】
ゴルフボールモデルを構成する全ての要素に関して同様の計算を行い、単位体積当たりのエネルギーの値を、図１６に示すコンター図で可視化した。図中、ＣＯＮＴＯＵＲ−ＬＥＶＥＬにて、エネルギーの値［ＭＰａ］と、それに対応する色の関係を示した。要素のエネルギーロスの値が約６．０ＭＰａより大きい部分は赤系の色で表示している一方、エネルギーロスの値が小さい部分は青系の色で表示した。ゴルフボールモデルにおいてアルミニウム製中空棒が衝突した点付近が、特にエネルギーのロスが大きいことが確認できた。これにより、ゴルフボールのどの部分でエネルギーロスが大きいかを一目で判断可能であることが確認できた。なお、このコンター図において、１つの要素に関して複数の色が存在するのは、要素のエネルギーロスの値を節点に振り分け、隣接する要素の値も加えられるために、そのようになった。
【００６８】
（スプリットホプキンソン棒試験機による材料物性の測定）
図１７は、本発明のゴルフボールのエネルギーロス評価方法において、有限要素法解析に必要な材料物性（ヤング率や損失係数等）を測定するスプリットホプキンソン棒試験機が示された模式的正面図である。
【００６９】
図１７に示すスプリットホプキンソン棒試験機は、打撃棒５１、入力棒５３及び出力棒５５を備えており、これらは直線上に配置されている。入力棒５３には、第一ひずみゲージ５７及び第二ひずみゲージ５９が取り付けられている。出力棒５５には、第三ひずみゲージ６１及び第四ひずみゲージ６３が取り付けられている。入力棒５３の後端５３ａと出力棒５５の前端５５ａとの間には、円柱状の試験片７０が挟持されている。
【００７０】
試験片７０は、ゴルフボールに用いられるゴム組成物等の材料を試験片の形状に成形したものでもよく、また、球状に成形されたゴルフボールから切り出されたものであってもよい。本実施形態では、試験片７０の長さ（すなわち入力棒５３の後端５３ａと出力棒５５の前端５５ａとの距離）は４ｍｍ、試験片７０の断面直径は１８ｍｍとしている。
【００７１】
打撃棒５１、入力棒５３及び出力棒５５はポリメチルメタアクリレート製の円柱であり、断面直径は２０ｍｍ、ヤング率は５３００ＭＰａ、比重は１．１９である。打撃棒５１の長さは、１００ｍｍである。入力棒５３及び出力棒５５（以下、この入力棒５３と出力棒５５とは、「応力棒」とも称される）の長さは、２０００ｍｍである。
第一ひずみゲージ５７は入力棒５３の後端５３ａから９００ｍｍの位置に取り付けられており、第二ひずみゲージ５９は入力棒５３の後端５３ａから６００ｍｍの位置に取り付けられている。また、第三ひずみゲージ６１は出力棒５５の前端５５ａから３００ｍｍの位置に取り付けられており、第四ひずみゲージ６３は出力棒５５の前端５５ａから６００ｍｍの位置に取り付けられている。
【００７２】
このように、図１７のスプリットホプキンソン棒試験機は、打撃棒５１、入力棒５３及び出力棒５５がポリメチルメタアクリレートからなる合成樹脂製であり、打撃棒５１及び入力棒５３が２０００ｍｍと大きく、しかも第一ひずみゲージ５７と入力棒５３の後端５３ａとの距離及び第二ひずみゲージ５９と入力棒５３の後端５３ａとの距離が大きいので、ゴルフボールに用いられるような架橋ゴム等の粘弾性を有する材料の材料物性（ヤング率や損失係数等）の測定に適している。
【００７３】
上記第一ひずみゲージ５７、第二ひずみゲージ５９、第三ひずみゲージ６１、第四ひずみゲージ６３として単軸プラスチック用ひずみゲージを用い、本実施形態では、株式会社共和電業製のＫＦＰ−５−３５０−Ｃ１−６５を用い、入力棒５３、出力棒５５の上記した位置に貼着している。これら第一ひずみゲージ５７乃至第４ひずみゲージ６３の入力棒５３及び出力棒５７への取付位置は、長さ方向において同一線上としている。
【００７４】
このスプリットホプキンソン棒試験機によるヤング率や損失係数等の材料物性の測定では、まず、打撃棒５１が入力棒５３の前端５３ａに、２５ｍ／ｓの速度で衝突する。これによって入射ひずみ波が生じ、この入射ひずみ波は入力棒５３の後端５３ａに向かって進む。この入射ひずみ波の一部は、入力棒５３の後端５３ａにおいて反射し、反射ひずみ波となって入力棒５３の前端５３ｂに向かって進む。入射ひずみ波の一部は、入力棒５３の後端５３ａから試験片７０を透過し、さらに出力棒５５に伝播して透過ひずみ波となり、出力棒５５の後端５５ｂに向かって進む。
【００７５】
入射ひずみ波は、第一ひずみゲージ５７及び第二ひずみゲージ５９によって実測される。実測された入射ひずみ波は、ローパスフィルターに通され、１０ｋＨｚ以上の高周波が除去される。さらに、入射ひずみ波の時刻歴は、そのベースライン値をゼロとするゼロ補正が施される。こうして得られた第一ひずみゲージ５７及び第二ひずみゲージ５９における時間軸ひずみのそれぞれがフーリエ変換され、周波数軸ひずみが求められる。この第一ひずみゲージ５７及び第二ひずみゲージ５９における周波数軸ひずみから、伝達関数が導出される。第一ひずみゲージ５７と入力棒５３の後端５３ａとの距離Ｘ１と、第二ひずみゲージ５９と入力棒５３の後端５３ａとの距離Ｘ２との比（Ｘ１：Ｘ２）が考慮されつつ、上記伝達関数に基づいて、入力棒５３の後端５３ａにおける周波数軸ひずみが推定される。この周波数軸ひずみがフーリエ逆変換されることにより、入力棒５３の後端５３ａにおける入射ひずみ波の時間軸ひずみ（ひずみの時刻歴）εｉが得られる。
【００７６】
同様に、入力棒５３の後端５３ａで反射して前端５３ｂに向かう反射ひずみ波が第二ひずみゲージ５９及び第一ひずみゲージ５７によって実測される。この実測された反射ひずみ波から、入力棒５３の後端５３ａにおける反射ひずみ波の時間軸ひずみ（ひずみの時刻歴）εｒが得られる。
【００７７】
また、出力棒５５の第三ひずみゲージ６１及び第四ひずみゲージ６３によって、試験片７０を経て出力棒５５に伝播される透過ひずみ波を実測し、この実測した透過ひずみ波から、出力棒５５の前端５５ａにおける透過ひずみ波の時間軸ひずみ（ひずみの時刻歴）εｔが得られる。
【００７８】
こうして得られたεｉ、εｒ及びεｔから、下記数式（５）によって、試験片７０のひずみ速度ε’が算出される。
（数式５）
0003597517.tif 000002
（数式（５）において、Ｃ０は入力棒および出力棒中（応力棒）のひずみ波の伝播速度（ｍ／ｓ）を表し、Ｌは試験片の長さ（ｍ）を表し、Ｅは応力棒のヤング率（Ｎ／ｍ２）を表し、ρは応力棒の密度（ｋｇ／ｍ３）を表す）
【００７９】
また、εｉ、εｒ及びεｔから、下記数式（６）によって試験片７０のひずみεが算出される。
【００８０】
（数式６）
0003597517.tif 000003
（数式（６）において、Ｃ０は入力棒および出力棒中（応力棒）のひずみ波の伝播速度（ｍ／ｓ）を表し、Ｌは試験片の長さ（ｍ）を表し、Ｅは応力棒のヤング率（Ｎ／ｍ２）を表し、ρは応力棒の密度（ｋｇ／ｍ３）を表す）
【００８１】
さらに、εｉ、εｒ及びεｔから、下記数式（７）によって試験片７０の応力σが算出される。
（数式７）
0003597517.tif 000004
（数式（７）において、Ｅは入力棒および出力棒からなる応力棒のヤング率（Ｎ／ｍ２）を表し、Ａは上記応力棒の断面積（ｍ２）を表し、Ａｓは試験片の断面積（ｍ２）を表し、Ｄは応力棒の直径（ｍ）を表し、Ｄｓは試験片の直径（ｍ）を表す）
【００８２】
こうして得られた試験片７０のひずみ時刻歴を、図１８のグラフに示す。図１８に示すように、曲線は、ピークＰ以降しばらくはなだらかであるが、その後、凹凸状となる。ピークＰ以降のなだらかな段階での点Ｓを選択し、この点Ｓにおける曲線に対する接線を画き、この接線と時間軸との交点から緩和時間λを導出し、下記数式（８）によって求められる曲線を点Ｓ以降の曲線とすることによって、ひずみ時刻歴全体をなだらかな曲線（図１８中に点線で示す）とすることができる。これにより、最終的に得られる材料物性（ヤング率や損失係数等）へのノイズの影響を除去することができる。
（数式８）
ε（ｔ）＝ε０・ｅ−ｔ／λ            −−−（８）
（数式（８）において、ε０は接点におけるひずみを表す）
【００８３】
同様に、下記数式（９）によって、応力時刻歴全体をなだらかな曲線とすることができ、これによって最終的に得られる材料物性（ヤング率や損失係数等）へのノイズの影響を除去することができる。
（数式９）
σ（ｔ）＝σ０・ｅ−ｔ／λ            −−−（９）
（数式（９）において、σ０は接点における応力を表す）
【００８４】
このような補正が行われた試験片７０のひずみ時刻歴及び応力時刻歴から、応力−ひずみ曲線が決定される。図１９は、典型的な応力−ひずみ曲線が示されたグラフである。この応力−ひずみ曲線から、下記の数式（１０）を用いて、試験片７０のヤング率Ｅｓが算出される。
（数式１０）
Ｅｓ＝σｍａｘ／εｍａｘ            −−−（１０）
【００８５】
また、図１９の応力−ひずみ曲線から、下記の数式（１１）を用いて、位相角δが算出される。
（数式１１）
δ＝ｓｉｎ−１（（σａ−σｂ）／σｍａｘ）            −−−（１１）
そして、この位相角δより、損失係数（ｔａｎδ）が算出される。
【００８６】
上記した打撃棒が入力棒の前端に２５ｍ／ｓの速度で衝突したときの試験片の最大ひずみ速度は２０００／ｓから２５００／ｓ程度であり、また、最大圧縮ひずみは０．１５から０．２５程度である。このひずみの変形挙動は、ゴルフボールが打撃された場合の変形挙動に近いものである。すなわち、このスプリットホプキンソン棒試験機によって測定されるヤング率や損失係数等の材料物性は、ゴルフボールが実際に打撃されたときに近い状態における材料物性である。
なお、打撃棒の衝突速度は、２５ｍ／ｓに限らず、打撃棒が衝突したときの試験片のひずみの変形挙動が、ゴルフボールが打撃された場合の変形挙動に近くなるような速度であればよい。
【００８７】
【発明の効果】
以上の説明より明らかなように、本発明によれば、打撃物の衝突（インパクト）によって生じるゴルフボール内部のエネルギーロスを実際に観察し評価することができる。従って、ボール内部のどの部分（例えば、ボール表面のカバー部分、あるいは、ボールの中心付近等）で、どの程度エネルギーロスが生じているかを予測、観察することができる。これにより、エネルギーロスをゴルフボールの反発性能向上に有用なデータとして用いることができる。
【００８８】
また、有限要素法を用いたシミュレーションにより評価を行っているため、実際にボールの試作を行う必要がなく、あらゆるタイプのゴルフボールの設計を容易に行うことができる。これにより、ゴルフボールの試作に要していた、費用や時間を削減することができる。
【００８９】
以上より、打撃物の衝突時に生じるエネルギーロスの少ない、つまり反発性能に関して最適なゴルフボールを短期間に提供することができる。
【図面の簡単な説明】
【図１】本発明のエネルギーロス評価方法を示すフローチャートである。
【図２】本発明のエネルギーロス評価方法のシミュレーションを行ったゴルフボールモデルの概略図である。
【図３】アルミニウム製中空棒モデルのゴルフボールモデルへの衝突を想定したシミュレーションにおける各モデルの状態を示し、（Ａ）は衝突前、（Ｂ）は衝突中、（Ｃ）は衝突後の状態図である。
【図４】（Ａ）（Ｂ）（Ｃ）は要素座標の決定方法を示す概略図である。
【図５】ゴルフボールモデル中の１要素の１成分における応力、ひずみの関係を示す図である。
【図６】要素中に発生するひずみ量が増加している際の仕事量の算出方法を示す図である。
【図７】要素中に発生するひずみ量が増加している際の仕事量の総和を示す図である。
【図８】要素中に発生するひずみ量が減少している際の仕事量の算出方法を示す図である。
【図９】要素中に発生するひずみ量が減少している際の仕事量の総和を示す図である。
【図１０】ゴルフボールモデルと打撃物との衝突時のゴルフボールモデルの各要素のエネルギーロスを可視化した状態を示す図である。
【図１１】ゴルフボールモデルのある１要素における要素座標系の１方向の応力−ひずみ曲線を示す。
【図１２】ひずみが０．１からｄεだけ進んだときのエネルギーロスの算出方法を示す。
【図１３】ひずみが０から（０．１＋ｄε）まで進んだときのエネルギーロスの算出方法を示す。
【図１４】実施例１のゴルフボールモデルとアルミニウム製中空棒の衝突時の状態を示す図である。
【図１５】実施例１におけるゴルフボールモデル中の１要素の１成分における応力、ひずみの関係を示す図である。
【図１６】実施例１におけるゴルフボールモデルと打撃物との衝突時のゴルフボールモデルの各要素のエネルギーロスを可視化した状態を示す図である。
【図１７】ヤング率や損失係数等の材料物性が測定されるスプリットホプキンソン棒試験機が示された模式的正面図である。
【図１８】試験片のひずみ時刻歴の状態が示されたグラフである。
【図１９】スプリットホプキンソン棒試験機による材料物性の測定時の応力−ひずみ曲線が示されたグラフである。
【符号の説明】
１０    ゴルフボールモデル
１１    要素
１２    節点
２０    アルミニウム製中空棒モデル
【図１】
0003597517.tif 000005
【図２】
0003597517.tif 000006
【図３】
0003597517.tif 000007
【図４】
0003597517.tif 000008
【図５】
0003597517.tif 000009
【図６】
0003597517.tif 000010
【図７】
0003597517.tif 000011
【図８】
0003597517.tif 000012
【図９】
0003597517.tif 000013
【図１０】
0003597517.tif 000014
【図１１】
0003597517.tif 000015
【図１２】
0003597517.tif 000016
【図１３】
0003597517.tif 000017
【図１４】
0003597517.tif 000018
【図１５】
0003597517.tif 000019
【図１６】
0003597517.tif 000020
【図１７】
0003597517.tif 000021
【図１８】
0003597517.tif 000022
【図１９】
0003597517.tif 000023