(57)【特許請求の範囲】
【請求項１】
粘弾性材料からなる製品の実使用状態を想定した測定条件下で、上記粘弾性材料に生じるひずみ、ひずみ速度、応力の時々刻々の値を測定し、
上記ひずみ、ひずみ速度、応力の時刻歴データと、上記粘弾性材料の粘性を考慮した粘弾性モデルとから、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで場合分けして上記粘弾性材料の粘性抵抗の時刻歴データを導出し、
上記粘弾性材料からなる製品を、解析対象の製品モデルとして設定し、該製品モデルを多数の要素に分割し、上記ひずみ、ひずみ速度、粘性抵抗の関係を入力し、上記ひずみの負荷状態時とひずみの除荷あるいは復元状態時とでの粘性抵抗の違いを考慮して有限要素法により解析を行い、上記粘弾性材料からなる製品モデルの性能を予測する際に、
上記要素毎に生じる全体座標系のひずみ及びひずみ速度を偏差成分と体積成分に分解し、上記偏差成分のひずみ及びひずみ速度を、全体座標系から主ひずみ座標系及び主ひずみ速度座標系に変換し、
上記変換した偏差主ひずみ及び偏差主ひずみ速度を用い、主ひずみ座標系及び主ひずみ速度座標系の各座標軸毎に、ひずみが同じ値でもひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定していることを特徴とする粘弾性材料からなる製品の性能予測のためのシミュレーション方法。
【請求項２】
上記要素毎に主ひずみの大きさであるノルムを計算し、シミュレーション時の前ステップのノルムと比較することによって、ノルムが増加している場合は負荷、ノルムが減少している場合は除荷と判断し、ひずみの負荷と除荷の状態を判定している請求項１に記載の粘弾性材料からなる製品の性能予測のためのシミュレーション方法。
【請求項３】
上記粘弾性材料は、材料物性が非線形性を示す材料である請求項１または請求項２に記載の粘弾性材料からなる製品の性能予測のためのシミュレーション方法。
【請求項４】
上記ひずみ、ひずみ速度、応力をスプリットホプキンソン棒試験機により測定している請求項１乃至請求項３のいずれか１項に記載の粘弾性材料からなる製品の性能予測のためのシミュレーション方法。
【請求項５】
上記測定時に粘弾性材料に生じる、上記ひずみの最大値が０．０５〜０．５０の範囲であり、あるいは／及び、上記ひずみ速度の最大値が５００／ｓ〜１００００／ｓの範囲である請求項１乃至請求項４のいずれか１項に記載の粘弾性材料からなる製品の性能予測のためのシミュレーション方法。
【請求項６】
上記粘弾性材料はゴルフボール用材料とし、上記製品モデルをゴルフボールとし、
上記ゴルフボールと、ゴルフクラブヘッドを想定した打撃物との衝突現象をシミュレーションし、該衝突時のゴルフボールの挙動を予測している請求項１乃至請求項５のいずれか１項に記載の粘弾性材料からなる製品の性能予測のためのシミュレーション方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、粘弾性材料からなる製品の性能予測のためのシミュレーション方法に関し、詳しくは、粘弾性材料からなる製品の性能をシミュレーションにより精度良く予測するシミュレーション方法に関するものである。
【０００２】
【従来の技術】
ゴムやエラストマーなどの高分子材料に代表される粘弾性材料は、タイヤ、スポーツ競技において使用される各種ボール、印刷機に用いられるロールなど各種製品に広く適用されている。
【０００３】
上記のような粘弾性材料や金属材料等の各種工業製品について、試作の費用と時間の節約等のためにシミュレーションを用いた製品開発が多方面で行われている。例えば、ゴルフボールの反発性能を予測するために、実際の打撃試験のシミュレーション方法が提案されている。
【０００４】
従来、上記のようなシミュレーションを行う場合、一般に材料の剛性、粘性などを測定する粘弾性スペクトルメーター、縦弾性係数（ヤング率）を測定する引張試験機などにより測定したボールの構成材料の物性値をシミュレーションへの入力値として用いている。特に、粘弾性スペクトルメーターでは、動的ひずみが与えられた試験片の物性値が測定されるため、粘弾性材料からなる製品のシミュレーションに有用である。
【０００５】
【発明が解決しようとする課題】
しかしながら、上記した粘弾性スペクトルメーター、縦弾性係数を測定する引張試験機などによる測定では試験片に大きな変形量を与えることができず、測定時に、粘弾性材料からなる試験片にかかる最大ひずみ速度は、０．００１／ｓから１．０／ｓ程度と小さく、また、最大圧縮ひずみも０．０００１から０．０２程度と小さな値となる。
【０００６】
一方、粘弾性材料からなる製品は、実使用時に外力の影響等で高速かつ大きな変形を起こす場合がある。例えば、ゴルフボールを構成する材料では、実際の打撃時に、最大ひずみ速度は５００／ｓ〜５０００／ｓ程度であり、最大圧縮ひずみは主に、０．０５〜０．５０程度であり、非常に大きな値となる。
【０００７】
このように、粘弾性スペクトルメーター、縦弾性係数を測定する引張試験機等による測定では、上記のような実使用条件と同様の高速かつ大きな変形条件下での粘弾性材料の物性値を測定することはできず、実使用時と、最大ひずみ速度及び最大圧縮ひずみが大きく異なる。このため、粘弾性スペクトルメーター、引張試験機等により得られた材料物性値を入力する従来のシミュレーションでは、粘弾性材料物性を考慮した正確なシミュレーションを行うことができないという問題がある。
【０００８】
即ち、衝撃荷重を受ける粘弾性材料の変形挙動は、静荷重を受ける場合とは異なり、変形量、あるいは変形速度によって著しく大きな影響を受けることが知られている。特にゴムやエラストマーなどに代表される高分子材料では、衝撃荷重により１万分あるいは１千分の数秒という非常に高速の変形を引き起こし、且つ変形量も数十％と非常に大きくなる。このような高速大変形の挙動を伴う粘弾性材料は多く、製品開発を効率的に行う上では精度の高いシミュレーションが要求されている。具体的には、ゴルフボール等のように使用時に衝撃を伴う製品では高速且つ大変形条件下での動的挙動および材料の特性が製品の性能を左右するため、実使用条件下での精度良いシミュレーションが製品開発に必要不可欠である。
【０００９】
また、粘弾性材料には、衝撃荷重等の外力が加えられた場合、損失係数等の物性値がひずみ、ひずみ速度の大きさにより変化する材料がある。即ち、粘弾性材料は、変形の速度、大きさが広範囲であるため、変形速度や大きさによって、物性値が一次直線的でなく、強い非線形的な変化を伴う性質がある。具体的には、外力により材料が変形し、ひずみが大きくなるに伴い、Ｓ−Ｓ（ひずみ−応力）曲線のループ面積が大きくなり、損失係数等の物性値が、材料の変形状態（変形の速度、大きさ）等により変化し非線形性を示す。粘弾性材料によっては、強い非線形性を持つ材料も多く、このような粘弾性材料を用いた製品についても精度の高いシミュレーションが要求されている。
【００１０】
しかしながら、粘弾性材料の非線形性を示す物性、例えば、損失係数が材料の変形の速度、大きさによって強く非線形的に変化する現象を精度良く表現できる手法がなかった。従来、粘弾性材料で構成されるゴルフボール等の製品については、粘弾性物性値がほとんど変化しないものとしてシミュレーションを行っていた。その結果、粘弾性材料からなる製品の実使用時の性能をシミュレーションにより正確に予測することができないという問題があり、実際には、試作品を製作することにより製品の性能を評価せざるを得ないという問題がある。
【００１１】
本発明は上記した問題に鑑みてなされたものであり、材料物性が非線形性を示す粘弾性材料からなる製品の実使用条件下での性能を、精度良くシミュレーションにより予測することを課題としている。
【００１２】
【課題を解決するための手段】
上記課題を解決するため、本発明は、粘弾性材料からなる製品の実使用状態を想定した測定条件下で、上記粘弾性材料に生じるひずみ、ひずみ速度、応力の時々刻々の値を測定し、
上記ひずみ、ひずみ速度、応力の時刻歴データと、上記粘弾性材料の粘性を考慮した粘弾性モデルとから、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで場合分けして上記粘弾性材料の粘性抵抗の時刻歴データを導出し、
上記粘弾性材料からなる製品を、解析対象の製品モデルとして設定し、該製品モデルを多数の要素に分割し、上記ひずみ、ひずみ速度、粘性抵抗の関係を入力し、上記ひずみの負荷状態時とひずみの除荷あるいは復元状態時とでの粘性抵抗の違いを考慮して有限要素法により解析を行い、上記粘弾性材料からなる製品モデルの性能を予測する際に、
上記要素毎に生じる全体座標系のひずみ及びひずみ速度を偏差成分と体積成分に分解し、上記偏差成分のひずみ及びひずみ速度を、全体座標系から主ひずみ座標系及び主ひずみ速度座標系に変換し、
上記変換した偏差主ひずみ及び偏差主ひずみ速度を用い、主ひずみ座標系及び主ひずみ速度座標系の各座標軸毎に、ひずみが同じ値でもひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定していることを特徴とする粘弾性材料からなる製品の性能予測のためのシミュレーション方法を提供している。
【００１３】
このように、実使用条件下で測定されたひずみ、ひずみ速度を体積成分と偏差成分に分解し、主ひずみ座標系及び主ひずみ速度座標系に変換した偏差成分の主ひずみ及び主ひずみ速度を用いて主ひずみ座標系及び主ひずみ速度座標系の各軸毎に、ひずみが同じ値の時でも、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定し、粘性抵抗を要素毎に各軸方向３成分について可変にしている。このため、材料の変形状態が負荷状態であるのか、除荷状態であるのか、実際の状況に対応した精度良い解析を行うことができると共に、３次元方向について正確に解析を行うことができ、シミュレーションの性能予測精度を向上することができる。
【００１４】
即ち、粘弾性材料で生じている変形状態は、加わるひずみが圧縮方向に増加する「負荷状態」と、圧縮量が徐々に減少する「復元状態」とに分けることができ、その状態によりひずみとひずみ速度の関係が異なってくる。従って、上記のように、粘弾性材料で生じているひずみ状態が負荷状態時と、除荷（あるいは復元）状態時との２つの状態を考慮し、負荷と除荷では異なったひずみ、ひずみ速度の関係を用いて粘性抵抗を決定し、ひずみがある同じ値のときでもひずみの負荷状態と除荷状態において異なる粘性抵抗値を導出することにより、実際の材料の変形状態を的確に捉えることができる。よって、材料が変形の速度、大きさ等により変化する現象を精度良く表現することができ、シミュレーションの精度をさらに向上させることができる。
【００１５】
また、上記各要素毎に主ひずみの大きさ（ノルム）を計算し、シミュレーション時の前ステップのノルムと比較することによって、ノルムが増加している場合は負荷、ノルムが減少している場合は除荷と判断し、ひずみの負荷と除荷の状態を判定していることが好ましい。
例えば、一方向の主ひずみだけでは、微小な変化で負荷と除荷が様々な変化を示すが、このように、多方向の大きさにすることで、ある程度の大きさの変化でないと負荷と除荷は変化しないため、計算の安定性を高めることができる。
【００１６】
本発明では、全体座標系でεｘ、εｙ、εｚ、εｘｙ、εｙｚ、εｚｘとする要素毎に生じるひずみ、ひずみ速度の偏差成分において粘弾性特性を考慮するために、上記ひずみ、ひずみ速度を体積成分と偏差成分に分解し、主ひずみ座標系及び主ひずみ速度座標系に変換した偏差成分の主ひずみ及び主ひずみ速度を用いて主ひずみ座標系及び主ひずみ速度座標系の各軸毎にひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定し、該粘性抵抗を要素毎に各軸方向３成分について可変にしている。このため、３次元方向について正確に粘弾性材料の粘弾性を考慮した解析を行うことができ、シミュレーションの精度を向上することができる。
ここで、体積成分とは、εｖ＝（εｘ＋εｙ＋εｚ）で表される成分を指し、偏差成分とは、εｘｙ’＝εｘｙ／２、εｙｚ’＝εｙｚ／２、εｚｘ’＝εｚｘ／２、εｘ’＝εｘ−εｖ／３、εｙ’＝εｙ−εｖ／３、εｚ’＝εｚ−εｖ／３で表される成分を指す。なお、せん断成分を１／２倍しているのは、物理ひずみを変換しているためである。
偏差成分のみを変換しているのは、体積成分の粘弾性成分は非常に小さく、また、現状では実測で求めることができないという理由による。即ち、体積成分の粘弾性成分はほとんどないので、偏差成分のみで判断し、変換は実測の結果から１軸方向の結果しかないために、εｘｙ、εｙｚ、εｚｘ成分＝０になるように変換し、ε１、ε２、ε３（主ひずみ３成分）のみの３成分について実測の結果を用いる。
なお、全体座標系とは、初期の形状、つまりモデルの形状を決定している座標系を指し、主ひずみ座標系及び主ひずみ速度座標系とは、ひずみ及びひずみ速度をせん断成分が０になるように変換した時、全体座標系を同じ角度変換した座標系を指す。
【００１７】
各要素毎に生じるひずみ（ひずみ速度）は、６成分であるが、高速大変形の条件下で測定される材料物性は、１軸方向の材料物性であるために、せん断方向の材料物性の考慮が不十分である。従って、各要素毎に生じるひずみ６成分を主ひずみ座標系に変換し、主軸方向のひずみ３成分を求める。
即ち、せん断方向のひずみを０にすることで全体座標系の６成分（εｘ、εｙ、εｚ、εｘｙ、εｙｚ、εｚｘ）を主軸方向のひずみ（主ひずみ座標系、ε１、ε２、ε３）３成分にする。これにより、せん断方向を考慮することなく、実測結果から粘弾性材料の材料特性を評価することができる。
【００１８】
具体的には、各要素毎に生じる全体座標系のひずみ及びひずみ速度を偏差成分と体積成分に分解し、偏差成分に粘弾性を考慮している。分解したひずみ、ひずみ速度の偏差６成分を、全体座標系から主ひずみ座標系及び主ひずみ速度座標系に変換し、主ひずみ座標系及び主ひずみ速度座標系での偏差主ひずみ、偏差主ひずみ速度の各３成分を用いている。これにより、主ひずみ座標系及び主ひずみ速度座標系で各軸毎に粘性抵抗を決定することができる。
【００１９】
また、上記粘性抵抗から粘弾性を考慮して、偏差成分の応力とひずみを計算し、該偏差成分の応力とひずみを元の全体座標系に再変換し、体積成分の応力とひずみ及び偏差成分の応力とひずみから、要素毎の応力とひずみを算出している。このように、要素毎の応力とひずみを算出することにより、時々刻々の応力、ひずみを算出することができる。
【００２０】
具体的には、下記の数式１、２より、偏差主ひずみ及び偏差主ひずみ速度を用い、各座標軸毎にひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定し、該粘性抵抗から粘弾性を考慮した偏差成分の応力とひずみを算出している。
（数式１）
0003910825.tif 000002
ここで、Ｅｉ：ｉ番目の要素の偏差主ひずみ、ｎ：解析時間ステップ、
βｉ：ｉ番目の要素のβ（ここでβ＝Ｅｔ／η、Ｅｔ＝Ｅｓｈｏｒｔ−Ｅｌｏｎｇ、η：粘性抵抗）、γｉ：ｉ番目の要素の偏差主ひずみ速度を指す。
（数式２）
0003910825.tif 000003
ここで、σｉ：ｉ番目の要素応力、Ｆｉ：剛性係数、Ｇｓｈｏｒｔ：短期せん断係数、Ｇｌｏｎｇ：長期せん断係数を指す。
数式１、２中のβ、剛性の係数Ｆを実験データより算出した後、数式１により、ひずみの偏差６成分の更新を行い、数式２により応力の計算を行う。この後、主ひずみ座標系、主ひずみ速度座標系から全体座標系に変換し、下記の数式３により応力の更新を行っている。
（数式３）
0003910825.tif 000004
ここで、Ｋ：体積弾性係数、εｖ：体積ひずみ、σｉｊ：全体座標系の応力を指す。
【００２１】
さらに、粘弾性材料の粘性を考慮した粘弾性モデルを用いひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を導出し、ひずみ、ひずみ速度、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なる粘性抵抗の関係を製品モデルに入力しシミュレーションを行っているため、粘弾性材料の非線形性を示す物性値が、材料の変形の速度、大きさによって非線形的に変化する現象を精度良く表現することができる。かつ、シミュレーションに用いられるひずみ、ひずみ速度、応力の値は、粘弾性材料からなる製品の実使用状態を想定した条件下で測定されているため、実際の粘弾性材料の種々の変形状態に対応したシミュレーションを行うことができる。従って、材料の変形状態により、ひずみ、ひずみ速度の関係が変化し、損失係数等の材料物性が非線形性を示すような粘弾性材料からなる製品においても、シミュレーションによる性能の精度良い予測を可能にしている。
【００２２】
本シミュレーション方法では、粘弾性材料からなる製品の実使用状態を想定した測定条件下で、上記粘弾性材料に生じるひずみ、ひずみ速度、応力の時々刻々の値を測定している。具体的には、実際の製品使用時に製品に外力が加わり、粘弾性材料が変形した状態を想定して測定条件を定めている。上記測定条件下で、粘弾性材料に生じるひずみ、ひずみ速度、応力の時々刻々の値を測定し、各時刻歴データを得ている。従って、上記ひずみ、ひずみ速度、応力の時刻歴データから、粘弾性材料に実使用状態を想定した外力を加わえた時の粘弾性材料の変形状態の情報を得ることができる。これにより、衝撃荷重よって高速大変形な挙動を伴う粘弾性材料についても正確に材料物性を予測することができる。
【００２３】
上記ひずみ、ひずみ速度、応力の時々刻々の値は、複数の異なる測定条件下で測定するのが好ましい。製品に加わる外力の大きさを変更し、複数のパターンの測定条件を設定することにより、ひずみ、ひずみ速度、応力について種々のパターンのデータを得ることができ、シミュレーションの入力値の精度を向上することができる。また、できる限り、多くのパターンのデータを得るために、上記ひずみ、ひずみ速度、応力の値は、粘弾性材料に外力を加えて、ひずみが生じてから、ひずみがほぼゼロになるまで測定するのが好ましく、シミュレーション精度の点より測定時間の間隔は短い方が良い。
【００２４】
粘弾性材料の粘性を考慮した粘弾性モデルとしては、バネとダッシュポットからなる粘弾性モデルが好ましい。このような粘弾性モデルにより、粘弾性材料の粘性を単純化することができるため、粘性が材料の変形状態に及ぼす影響を容易に考慮することができる。具体的には、マックスウェル（maxwell）モデル、フォークト（Voigt）モデルあるいは、さらに、複数のバネ、ダッシュポットの組み合わせが挙げられる。モデルの単純化という観点からは２要素モデルが好ましい。このような粘弾性材料モデルでは、ダッシュポットの粘性抵抗ηとバネの剛性（縦弾性係数Ｅあるいはせん断係数（横弾性係数））を可変にして用いている。せん断係数とは、縦弾性係数（ヤング率）Ｅとポアソン比とにより決定される物性値である。
【００２５】
上記ひずみ、ひずみ速度、応力の各時刻歴データと粘弾性材料の粘性を考慮した粘弾性モデルとから、粘弾性材料のひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗の時刻歴データを導出している。具体的には、上記粘弾性モデルにより、粘弾性材料に生じるひずみ、ひずみ速度、及び粘弾性材料の縦弾性係数、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗の各関係を定式化している。このように定式化することで、粘性抵抗をひずみ、ひずみ速度の関数として表している。そして、上記測定により得られた粘弾性材料に生じるひずみと応力により粘弾性材料の縦弾性係数を導出し、ひずみ、ひずみ速度と共に、上記関数に代入することでひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗の値を導出している。上記測定にて、ひずみ、ひずみ速度、応力についての時刻歴データを得ているため、粘性抵抗についても同様に、対応する時刻歴データを得ることができる。このように、粘性抵抗は、ひずみ、ひずみ速度の値によって決定され、その時間変化に対応して粘性抵抗も変化する。
【００２６】
上記粘弾性材料からなる製品を解析対象の製品モデルとして設定し、該製品モデルを含み、速度、拘束条件等を含んだ計算用入力データ（あるいはinput data）に上記ひずみ、ひずみ速度、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗の関係を入力し、上記ひずみ、ひずみ速度の違いによる粘性抵抗の変化を考慮したシミュレーションを行っている。
シミュレーション計算時に製品モデルに上記ひずみ、ひずみ速度、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗の関係を入力しているが、具体的には、ひずみとひずみ速度の関係、及び、ひずみと粘性抵抗の関係の各２次元データを入力し、演算することができる。また、ひずみとひずみ速度と粘性抵抗の関係を３次元データとし、粘性抵抗を関数化したデータとして入力し、演算することもできる。
【００２７】
ひずみとひずみ速度の関係、及び、ひずみと粘性抵抗の関係の各２次元データは、上記各関係を用いて、ひずみ、ひずみ速度とそれに対応するひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を入力データとして書き込むことになる。具体的には、複数の測定条件にて、ひずみ、ひずみ速度等の測定を行い、測定条件の異なる各パターンについて、ひずみ、ひずみ速度の時系列データから、各々ひずみとひずみ速度の対応関係を記録する。各々の曲線に対応する粘性抵抗の値も記録する。これらひずみ、ひずみ速度、粘性抵抗の３つの対応を整理することで、あるひずみ、ひずみ速度の時の粘性抵抗を精度良く導出し、計算を行っている。
【００２８】
異なるひずみ、ひずみ速度の条件下での測定数が多いほど精度良く粘弾性材料の物性を実現できるため、上記したように複数の異なった測定条件にてひずみ、ひずみ速度等の測定を行うのが好ましい。ただし、測定データが多いほど、シミュレーションの際、計算時間を要するため、定められた測定条件で測定した測定データと同一でないひずみ、ひずみ速度の場合には、補間を用いて粘性抵抗を算出するのが好ましい。上記補間は、例えば、近接するひずみ、ひずみ速度の条件下での粘性抵抗２値からの１次補間、あるいは測定した全ての条件下での個々の測定値を用いての補間等の種々の方法により行うことができる。このような補間作業を行うことにより、測定条件の違いによる材料に生じるひずみとひずみ速度の変化に対応して粘性抵抗を算出することができる。
【００２９】
上記粘弾性材料は、材料物性が非線形性を示す材料であっても本シミュレーション方法によれば、実使用状態を想定した材料物性及び変形挙動を正確にシミュレーションすることができる。粘弾性材料の材料物性が、一次直線的でなく、材料の変形速度や変形の大きさによって変化を伴う強い非線形性を示す材料においても、上記のように、粘性抵抗をひずみ、ひずみ速度の２つの値によって決定するような粘弾性モデルを用い粘性抵抗を算出することで、非線形性を考慮することができる。本シミュレーション方法では、特に、損失係数が強い非線形性を示す材料において、精度良くシミュレーションを行い、実使用状態を想定した製品の性能を予測することができる。
【００３０】
本発明のシミュレーション方法は、以下に示すように有限要素法により行っている。
上記シミュレーションを有限要素法解析により行う際は、上記製品モデルに多数の節点と要素を設定している。即ち、有限要素法で粘弾性材料で構成される製品をシミュレーションして材料物性を予測する際に、粘弾性モデルのバネの剛性（縦弾性係数あるいはせん断係数）を決定し、粘性を示す粘性抵抗を上記に示すように各要素ごとにその要素で生じているひずみ、ひずみ速度によって決定することで各要素ごとに適正なひずみ、ひずみ速度の条件下の材料物性を表すことができる。なお、縦弾性係数を入力する代わりに、ポアソン比との関係から、せん断係数で入力しても構わない。縦弾性係数で入力するか、せん断係数で入力するかは、有限要素法プログラムの仕様により変更できる。
【００３１】
上記ひずみ、ひずみ速度、応力をスプリットホプキンソン棒試験機により測定している。スプリットホプキンソン棒試験機を用いると、試験片に高速でかつ大変形のひずみを与えることができ、１万分あるいは１千分の数秒という高速で且つ変形量も数十％という、高速大変形条件下での粘弾性材料のひずみ、ひずみ速度、応力の時系列データを得ることができる。スプリットホプキンソン棒試験機の測定条件を、製品に衝撃荷重が加わった時の粘弾性材料に生じるひずみ、ひずみ速度と同等の条件とすれば、外力により生じた粘弾性材料の高速大変形状態等、種々の状態に対応した粘弾性材料の材料物性を得ることができる。従って、このスプリットホプキンソン棒試験機で測定された材料物性を用いることで、シミュレーションの精度を向上することができる。
【００３２】
また、スプリットホプキンソン棒試験機は、試験片に衝撃を加える打撃棒の衝突速度を変更するだけで、様々なひずみ、ひずみ速度の領域での試験片の材料物性を測定することができるため、容易に、種々のひずみ、ひずみ速度のパターンにおいて、材料物性を得ることができる。
【００３３】
スプリットホプキンソン棒試験機は、元来、金属材料の衝撃挙動の評価用に用いられていたが、本発明では、スプリットホプキンソン棒試験機を、粘性を持つ粘弾性材料の評価用に改良して用いている。スプリットホプキンソン棒試験機の測定方法等については後述する。
【００３４】
なお、試験片に高速でかつ大変形のひずみを与えることができ、製品が実際に使用される条件下でのひずみ、ひずみ速度の条件下で、材料物性を測定することができれば、上記スプリットホプキンソン棒以外の測定方法により材料物性を求めても良いことは言うまでもない。
【００３５】
上記粘弾性材料からなる製品の実使用状態を想定した測定条件下でのひずみ、ひずみ速度、応力の測定時に、粘弾性材料に生じる、上記ひずみの最大値は０．０５〜０．５０の範囲であり、あるいは／及び、上記ひずみ速度の最大値は５００／ｓ〜１００００／ｓ、好ましくは５００／ｓ〜５０００／ｓの範囲であるのが良い。上記ひずみの最大値と、ひずみ速度の最大値の範囲は、粘弾性材料の高速大変形時に生じるひずみ、ひずみ速度の条件であるため、高速大変形時の製品の性能を予測するには、この条件下のひずみ、ひずみ速度、応力の３つの時系列データを用いることが好ましい。
【００３６】
上記粘弾性材料はゴルフボール用材料とし、上記製品モデルをゴルフボールとしている。ゴルフボールは、粘弾性材料から構成され、実使用時に衝撃荷重等の外力を受け高速変形、あるいは大変形を強いられる製品であり、その高速、大変形状態がゴルフボール自体の性能に大きく影響を及ぼす。このため、本発明のシミュレーション方法による解析は、ゴルフボールの性能予測に非常に有用であり、試作をせずに、ゴルフボールの性能を精度良く予測することができるため、ゴルフボールの設計の効率化を図ることができる。
【００３７】
上記ゴルフボールと、ゴルフクラブヘッドを想定した打撃物との衝突現象をシミュレーションし、該衝突時のゴルフボールの挙動を予測している。実際にゴルフクラブヘッドにより、ゴルフボールが打撃されたときと同等のひずみ、ひずみ速度、応力の条件下で、ゴルフボールを構成する粘弾性材料の材料物性を予測することができるため、ゴルフボールの反発係数やゴルフボール打撃時のゴルフボールの変形状態等、ゴルフボールの挙動を予測することができる。
【００３８】
上記ゴルフボールは、架橋ゴム層等の単体からなる所謂１ピースゴルフボールであってもよく、架橋ゴム層等のコアにカバーが被覆された所謂２ピースゴルフボールでもよく、また、３層以上から構成されていてる所謂マルチピースゴルフボールであってもよく、粘弾性材料からなるあらゆる構造のゴルフボールに適応可能である。
【００３９】
上記粘弾性材料としては、粘弾性を有するあらゆる材料が含まれる。例えば、熱可塑性樹脂、熱硬化性樹脂、各種エラストマー、ゴム等が挙げられ、これらの単体、あるいは混合物を用いることができる。また、これらの単体、混合物には、着色剤、劣化防止剤、架橋剤等の各種添加剤が、必要に応じ配合されていてもよい。製品の実使用条件下で、ひずみ、ひずみ速度、応力を測定できる材料であれば、あらゆる材料に適応可能である。
【００４０】
具体的には、ゴルフボールの構成材料として用いられるアイオノマー樹脂等の合成樹脂、ポリブタジエン（ブタジエンゴム）、天然ゴム、ポリイソプレン、スチレン−ブタジエン共重合体、エチレン−プロピレン−ジエン共重合体（ＥＰＤＭ）、ウレタンゴム等が挙げられる。
【００４１】
粘弾性材料からなる製品としては、ゴルフボール以外に、プリンタ等の印刷機に用いられるゴムローラ、ゴムベルト、あるいはタイヤ、その他粘弾性材料からなるテニス用、ゴルフ用等のスポーツ用品等が挙げられる。粘弾性材料からなる製品としては、製品の少なくとも一部に粘弾性材料が使用されていれば良く、金属材料等の他の材料との複合成形品等であっても良く、製品の粘弾性材料に該当する部分の性能を予測することもできる。特に、使用時において、衝撃荷重を受け、高速で大きな変形量を示すような製品に好適であり、本シミュレーション方法により製品の性能や動的挙動を精度良く予測することができる。
【００４２】
また、本シミュレーション方法においては、上記製品モデルを含む計算用入力データに上記ひずみ、ひずみ速度、縦弾性係数あるいはせん断係数を入力し、上記ひずみ、ひずみ速度の違いによる縦弾性係数あるいはせん断係数の変化を考慮することにより、縦弾性係数が材料の変形状態（ひずみ、ひずみ速度の大きさ等）により変化する粘弾性材料についても正確な解析を行うことができ、シミュレーションの精度をさらに向上させることができる。
また、場合によっては、応力の負荷状態と除荷状態で場合分けすることもできる。
【００４３】
【発明の実施の形態】
以下、本発明の実施形態を図面を参照して説明する。
本発明の第一実施形態のシミュレーション方法では、非線形性を示す粘弾性材料として、ゴルフボールの構成材料であるブタジエンゴムを主成分とする材料を用いている。上記ブタジエンゴムを主成分とする材料を試験片として、粘弾性材料測定用に改良したスプリットホプキンソン棒試験機によって、ゴルフボールの実使用状態を想定した高速で大きな変形時の測定条件下で、上記ブタジエンゴムを主成分とする材料に生じるひずみ、ひずみ速度、応力の時々刻々の値を測定している。なお、スプリットホプキンソン棒試験機による測定方法については、後述する。
【００４４】
スプリットホプキンソン棒試験機は打撃棒の衝突速度を変更することで様々なひずみ、ひずみ速度の領域での材料物性を測定することができる。本実施形態では、４パターンの衝突速度（７ｍ／ｓ、１４ｍ／ｓ、２０ｍ／ｓ、２５ｍ／ｓ）を採用し、４つの測定条件下で測定を行い、それぞれの衝突パターンにてひずみ、ひずみ速度、応力の時刻歴データを得ている。上記ひずみε、ひずみ速度ε’、応力σの各時刻歴データを、それぞれ図１、図２、図３に示す。
【００４５】
スプリットホプキンソン棒試験機により打撃棒を衝突させた後、図１に示すように、ブタジエンゴムを主成分とする材料にはひずみが生じており、一定時間まではひずみの値は増加し、その後、なだらかに減少している。また、図２に示すように、ひずみ速度はひずみが最大値に達するまでは正の値を示し、ひずみの最大値以降は、負の値を示している。図３に示すように、応力の値も同様に、時間と共に変化している。
【００４６】
上記ひずみの時刻歴データと応力の時刻歴データからひずみ−応力線図を描くと図４のようになる。図４のひずみ−応力線図において、ひずみの値が最大時のひずみεｍａｘと、その時の応力σａの値を用い、下記の数式４により、各衝突パターン毎に試験片の縦弾性係数Ｅを算出している。
（数式４）
Ｅ＝σａ／εｍａｘ      −−−（４）
なお、本実施形態で用いたブタジエンゴムを主成分とする材料では、スプリットホプキンソン棒試験機において打撃棒の衝突速度を変更しても、各衝突パターンでの縦弾性係数はほぼ同じ値を示している。
【００４７】
粘弾性材料からなる製品について、粘性を考慮したシミュレーションを行うために、粘弾性材料の粘性を考慮した粘弾性モデルを設定している。具体的に、本実施形態では、図５に示すような、バネとダッシュポットからなる粘弾性モデルとして、基本的な２要素のVoigtモデルを用いている。即ち、上記粘弾性材料モデルでは、ダッシュポットの粘性抵抗ηとバネの剛性（縦弾性係数Ｅあるいはせん断係数）を可変にして用いている。
【００４８】
図５に示すように２要素のフォークト（Voigt）モデルとした粘弾性モデルにて、バネで生じる応力をσ１, ダッシュポットで生じる応力をσ２とすると全体で生じる応力σは以下のような数式５で表すことができる。
（数式５）
σ＝σ１＋σ２
＝Ｅε＋ηε’      −−−（５）
【００４９】
従って、上記数式５より、このような粘弾性モデルにて粘性抵抗ηは、以下の数式６で表すことができる。
（数式６）
η=(σ−Ｅε)/ε’      −−−（６）
【００５０】
従って、図１、図２、図３に示すひずみε、ひずみ速度ε’、応力σの時刻歴データと上記の数式３により、ひずみε、ひずみ速度ε’に対応する粘性抵抗ηを時々刻々算出できる。ここで、上記したように、上記ひずみ、ひずみ速度、応力の値は、時間の経過と共に変化しており、粘弾性材料で生じているひずみ状態は、加わるひずみが圧縮方向に増加する「負荷状態」と、圧縮量が徐々に減少する「復元状態」とに分けることができる。よって、粘性抵抗を算出する際にも、粘弾性材料で生じているひずみ状態が負荷状態時と、除荷（あるいは復元）状態時との２つに分けて算出し、図６に示すように、各衝突速度において粘性抵抗の時刻歴データを得ている。
【００５１】
このように、上記高速大変形条件下での測定により得られたひずみと応力より、粘弾性材料の縦弾性係数Ｅを決定した後、粘性を考慮した粘弾性モデルから導出した上記数式６により、負荷状態時と、除荷（あるいは復元）状態時とに分けて粘性抵抗ηを算出し、このひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗ηをシミュレーションに適用する。
【００５２】
上記ひずみの時刻歴データと粘性抵抗の時刻歴データより、ひずみと粘性抵抗との関係を図７に示す。図７に示すように、粘性抵抗の値は、ひずみの負荷状態時と除荷（あるいは復元）状態時とでは、ひずみが同じ値のときでも異なった値を示している。また、ひずみの時刻歴データとひずみ速度の時刻歴データより、ひずみとひずみ速度との関係を図８に示す。このように、粘性抵抗は、ひずみ、ひずみ速度によって変化し、かつ、ひずみの負荷状態時と除荷（あるいは復元）状態時とで異なる値を示すため、粘性抵抗はひずみ、ひずみ速度によって決定することができ、粘性抵抗をひずみの負荷状態時と除荷状態時とで場合分けして、ひずみ、ひずみ速度の関数としている。なお、図８において、４つの衝突速度の各データは、ひずみ速度の負の領域では、表示上一直線に重なっているように見えるが、実際の値は、各衝突速度毎にそれぞれ異なるひずみ、ひずみ速度の値となっている。
【００５３】
図７に示すひずみと粘性抵抗の関係と、図８に示すひずみとひずみ速度の関係とを、後述するように各衝突パターンにおいて、それぞれ、入力データとして書き込み、有限要素法解析により上記ひずみ、ひずみ速度の違いによる粘性抵抗の変化、即ち、ひずみの負荷状態時と除荷（あるいは復元）状態時との違いによる粘性抵抗の変化を考慮したシミュレーションを行っている。
【００５４】
本実施形態では、粘弾性材料からなる製品モデルとしてゴルフボールモデルを設定し、ゴルフクラブヘッド（打撃物）のゴルフボールへの衝突を想定したシミュレーションを行っている。図９に示すように、ゴルフボールモデル１０は、ブタジエンゴムを主成分とする１ピースボールを想定し、直径４２．８ｍｍの球状としている。
【００５５】
有限要素法解析を行うにあたり、製品モデルにおいて初期条件を設定している。即ち、ゴルフボールモデル１０の大きさ、形状、構造等の初期条件を設定すると共に、ゴルフボールモデル１０を多数の要素１１にメッシュ分割し、多数の節点１２を得ている。ゴルフボールモデルの全ての範囲の要素の数は、全体モデルで１０００〜１０００００が好ましく、さらに好ましくは、２５００〜２００００が良い。なお、上限の値は、現段階での計算機の能力を鑑みてのものであり、今後計算機の能力が向上するにつれ、解析の時間が短縮されるため上限の範囲は、変わることが容易に想像できる。
【００５６】
上記設定条件に基づいて、シミュレーション計算時に、上記ひずみ、ひずみ速度、粘性抵抗の関係とゴルフボールモデル１０のデータを入力データに書き込む。計算が実行されると時々刻々各要素に関して妥当な粘性抵抗を算出し、その粘性抵抗を用い式の計算を行っている。本シミュレーションでは、ひずみの負荷状態時と除荷（あるいは復元）状態時とで場合分けされたひずみとひずみ速度の関係、及び、ひずみと粘性抵抗の関係の各２次元データを入力し演算している。測定条件の異なる４つの各パターンについてのひずみとひずみ速度の時系列データから、各々ひずみとひずみ速度の対応関係を記録し、ひずみ、ひずみ速度と、それに対応する粘性抵抗を入力データとして書き込むことになる。４つの測定条件により測定されたひずみ、ひずみ速度の測定データと同一でないひずみ、ひずみ速度の場合には、近接するひずみ、ひずみ速度の条件下での粘性抵抗２値を用いて１次補間を行っている。
【００５７】
さらに具体的には、一つの要素に着目したとき、その要素に生じているひずみの大きさ（ノルム）を計算し、シミュレーション時の前ステップのノルムと比較することによって、ノルムが増加している場合は負荷、ノルムが減少している場合は除荷と判断し、ひずみの負荷と除荷の状態を判定し、材料状態を考慮したひずみとひずみ速度の情報を得ている。次に、測定から得られた各衝突ケース（測定ケース）のひずみ、ひずみ速度の関係から、各ケースでの同値のひずみが生じているときのひずみ速度の値を参照して、着目する要素のひずみ速度の値に近いひずみ速度の値を持つ２つのケースを探索する。この２つのケースに対応するひずみが同値のときの，ひずみ速度と粘性抵抗の値を用いて、補間により、その要素で生じるひずみ、ひずみ速度に対応した適切な粘性抵抗の値を算出することができる。本実施形態では、このように単純な一次補間による補間を行っている。このように、各要素ごとにその要素で生じているひずみ、ひずみ速度によって負荷と除荷では異なった粘性抵抗を決定することで、各要素ごとに適正なひずみ、ひずみ速度の条件下の材料物性を表している。
【００５８】
また、上記要素毎に生じる全体座標系のひずみ及びひずみ速度を偏差成分と体積成分に分解し、上記偏差成分のひずみ及びひずみ速度を、全体座標系から主ひずみ座標系及び主ひずみ速度座標系に変換し、変換した偏差主ひずみ及び偏差主ひずみ速度を用い、上述したように図７、８より各座標軸毎に、ひずみが同じ値でもひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定し、該数式１、２、３により粘性抵抗から粘弾性を考慮した偏差成分の応力とひずみを算出している。その後、上記偏差成分の応力とひずみを全体座標系に再変換し、上記体積成分の応力とひずみ、及び偏差成分の応力とひずみから、要素毎の応力とひずみを算出している。
【００５９】
即ち、各要素毎に生じるひずみ６成分を主ひずみ座標系に変換し、主軸方向のひずみ３成分を求めている。各要素毎に生じる全体座標系のひずみ及びひずみ速度を偏差成分と体積成分に分解し、偏差成分に粘弾性を考慮している。分解したひずみ、ひずみ速度の偏差６成分を、全体座標系から主ひずみ座標系及び主ひずみ速度座標系に変換し、主ひずみ座標系及び主ひずみ速度座標系での偏差主ひずみ、偏差主ひずみ速度の各３成分を用い、主ひずみ座標系及び主ひずみ速度座標系で各軸毎に、ひずみが同じ値でもひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定している。
【００６０】
本実施形態では、具体的に、図１０（Ａ）（Ｂ）（Ｃ）に示すように、ゴルフボールモデル１０に、打撃物として２００ｇ（ゴルフクラブヘッドの重さと同一）の円筒形状のアルミニウム製中空棒モデル２０を速度４５ｍ／ｓで衝突させた時のゴルフボールモデル１０の材料物性を有限要素法により解析し、シミュレーションを行っている。これにより、打撃物の衝突時から所定時間後までの所定時間におけるゴルフボールモデル１０の各要素１１に発生するひずみ量を演算している。
【００６１】
１要素に関して構成する各節点に質量を分配し、各節点を質点と置き換え、各節点の速度を質点の持つ速度として合計を割って、全体の速度としている。即ち、インパクト後のボールの速度Ｖｂｉ（ｉ＝ｘ、ｙ、ｚ）は、数式７により以下のように算出する。ボールの全運動量は全質点の運動量の和と考え、全運動量を総重量で割ったものをボールの飛び出し速度Ｖｂｉと定義する。
【００６２】
（数式７）
0003910825.tif 000005
ここで、Ｎ、Ｍは、ボールの節点数、総質量、Ｖｎｉ、Ｍｎはｎ番目の並進速度、節点を含む要素の質量をその要素に含まれる節点の数で割ったものである。
【００６３】
なお、アルミニウム製中空棒モデル２０の円形面２０ａを衝突面とし、衝突面はフラットであり、ゴルフボールモデル１０とフラットな正面衝突としている。また、アルミニウム製中空棒モデル２０の円形面２０ａの中心点２０ｂが、最初のゴルフボールモデル１０との衝突点となるように衝突させている。
【００６４】
上記した方法により、衝突前後の上記アルミニウム製中空棒モデル２０及びゴルフボール１０の速度を算出し、それぞれの速度及び重量からゴルフボール１０の反発係数を算出し予測している。
【００６５】
このように、上記ひずみ、ひずみ速度を体積成分と偏差成分に分解し、主ひずみ座標系及び主ひずみ速度座標系に変換した偏差成分の主ひずみ及び主ひずみ速度を用いて主ひずみ座標系及び主ひずみ速度座標系の各軸毎に、ひずみが同じ値でもひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗を決定し、上記粘性抵抗を要素毎に各軸方向３成分について可変にしている。このため、３次元方向について正確に粘弾性材料の粘弾性を考慮した解析を行うことができ、シミュレーションの精度を向上することができる。
【００６６】
また、ひずみの負荷状態時と除荷（あるいは復元）状態時とで場合分けされた粘性抵抗をひずみ、ひずみ速度の関数として決定し、上記ひずみ、ひずみ速度、粘性抵抗の３者の関係をゴルフボールモデルに入力して、有限要素法解析によりシミュレーションしている。このため、各要素の時々刻々のひずみ、ひずみ速度からそれに相当する粘性抵抗を算出でき、該粘性抵抗値はひずみの負荷状態時と除荷（あるいは復元）状態時との違いが考慮されており、高速大変形条件下においても強い非線形性を持つ粘弾性材料の特性を非常に精度良く予測することができる。
【００６７】
従って、実際のゴルフクラブヘッドによる打撃時と同等の条件下における、粘弾性材料からなるゴルフボールモデルの反発係数等の性能、動的挙動等を、試作をすることなく正確に把握することができる。
【００６８】
また、図１１示すように、応力−ひずみ曲線から、下記の数式８を用いて、位相角δが算出され、この位相角δより、損失係数（ｔａｎδ）を算出することもできる。
（数式８）
δ＝ｓｉｎ−１（（σａ−σｂ）／σｍａｘ）      −−−（８）
【００６９】
（スプリットホプキンソン棒試験機による材料物性の測定）
図１２は、粘弾性材料測定用に改良したスプリットホプキンソン棒試験機が示された模式的正面図である。
【００７０】
図１２に示すスプリットホプキンソン棒試験機は、打撃棒５１、入力棒５３及び出力棒５５を備えており、これらは直線上に配置されている。入力棒５３には、第一ひずみゲージ５７及び第二ひずみゲージ５９が取り付けられている。出力棒５５には、第三ひずみゲージ６１及び第四ひずみゲージ６３が取り付けられている。入力棒５３の後端５３ａと出力棒５５の前端５５ａとの間には、円柱状の試験片７０が挟持されている。
【００７１】
試験片７０は、測定すべき粘弾性材料を試験片の形状に成形したものでもよく、また、粘弾性材料を成形した製品状態から切り出されたものであってもよい。本実施形態では、試験片７０の長さ（すなわち入力棒５３の後端５３ａと出力棒５５の前端５５ａとの距離）は４ｍｍ、試験片７０の断面直径は１８ｍｍとしている。
【００７２】
打撃棒５１、入力棒５３及び出力棒５５はポリメチルメタアクリレート製の円柱であり、断面直径は２０ｍｍ、縦弾性係数は５３００ＭＰａ、比重は１．１９である。打撃棒５１の長さは、１００ｍｍである。入力棒５３及び出力棒５５（以下、この入力棒５３と出力棒５５とは、「応力棒」とも称される）の長さは、２０００ｍｍである。
第一ひずみゲージ５７は入力棒５３の後端５３ａから９００ｍｍの位置に取り付けられており、第二ひずみゲージ５９は入力棒５３の後端５３ａから６００ｍｍの位置に取り付けられている。また、第三ひずみゲージ６１は出力棒５５の前端５５ａから３００ｍｍの位置に取り付けられており、第四ひずみゲージ６３は出力棒５５の前端５５ａから６００ｍｍの位置に取り付けられている。
【００７３】
このように、図１２のスプリットホプキンソン棒試験機は、打撃棒５１、入力棒５３及び出力棒５５がポリメチルメタアクリレートからなる合成樹脂製であり、打撃棒５１及び入力棒５３が２０００ｍｍと大きく、しかも第一ひずみゲージ５７と入力棒５３の後端５３ａとの距離及び第二ひずみゲージ５９と入力棒５３の後端５３ａとの距離が大きいので、ゴルフボールに用いられるような架橋ゴム等の粘弾性を有する材料のひずみ、ひずみ速度、応力の測定に適している。
【００７４】
上記第一ひずみゲージ５７、第二ひずみゲージ５９、第三ひずみゲージ６１、第四ひずみゲージ６３として単軸プラスチック用ひずみゲージを用い、本実施形態では、株式会社共和電業製のＫＦＰ−５−３５０−Ｃ１−６５を用い、入力棒５３、出力棒５５の上記した位置に貼着している。これら第一ひずみゲージ５７乃至第４ひずみゲージ６３の入力棒５３及び出力棒５７への取付位置は、長さ方向において同一線上としている。
【００７５】
このスプリットホプキンソン棒試験機によるひずみ、ひずみ速度、応力の測定では、まず、打撃棒５１が入力棒５３の前端５３ａに、所定の速度で衝突する。これによって入射ひずみ波が生じ、この入射ひずみ波は入力棒５３の後端５３ａに向かって進む。この入射ひずみ波の一部は、入力棒５３の後端５３ａにおいて反射し、反射ひずみ波となって入力棒５３の前端５３ｂに向かって進む。入射ひずみ波の一部は、入力棒５３の後端５３ａから試験片７０を透過し、さらに出力棒５５に伝播して透過ひずみ波となり、出力棒５５の後端５５ｂに向かって進む。
【００７６】
入射ひずみ波は、第一ひずみゲージ５７及び第二ひずみゲージ５９によって実測される。実測された入射ひずみ波は、ローパスフィルターに通され、１０ｋＨｚ以上の高周波が除去される。さらに、入射ひずみ波の時刻歴は、そのベースライン値をゼロとするゼロ補正が施される。こうして得られた第一ひずみゲージ５７及び第二ひずみゲージ５９における時間軸ひずみのそれぞれがフーリエ変換され、周波数軸ひずみが求められる。この第一ひずみゲージ５７及び第二ひずみゲージ５９における周波数軸ひずみから、伝達関数が導出される。第一ひずみゲージ５７と入力棒５３の後端５３ａとの距離Ｘ１と、第二ひずみゲージ５９と入力棒５３の後端５３ａとの距離Ｘ２との比（Ｘ１：Ｘ２）が考慮されつつ、上記伝達関数に基づいて、入力棒５３の後端５３ａにおける周波数軸ひずみが推定される。この周波数軸ひずみがフーリエ逆変換されることにより、入力棒５３の後端５３ａにおける入射ひずみ波の時間軸ひずみ（ひずみの時刻歴）εｉが得られる。
【００７７】
同様に、入力棒５３の後端５３ａで反射して前端５３ｂに向かう反射ひずみ波が第二ひずみゲージ５９及び第一ひずみゲージ５７によって実測される。この実測された反射ひずみ波から、入力棒５３の後端５３ａにおける反射ひずみ波の時間軸ひずみ（ひずみの時刻歴）εｒが得られる。
【００７８】
また、出力棒５５の第三ひずみゲージ６１及び第四ひずみゲージ６３によって、試験片７０を経て出力棒５５に伝播される透過ひずみ波を実測し、この実測した透過ひずみ波から、出力棒５５の前端５５ａにおける透過ひずみ波の時間軸ひずみ（ひずみの時刻歴）εｔが得られる。
【００７９】
こうして得られたεｉ、εｒ及びεｔから、下記数式９によって、試験片７０のひずみ速度ε’が算出される。
（数式９）
ε’＝（Ｃ０／Ｌ）・（εｉ−εｒ−εｔ）
＝（（Ｅ／ρ）１／２／Ｌ）・（εｉ−εｒ−εｔ）    −−−（９）
（数式９において、Ｃ０は入力棒および出力棒中（応力棒）のひずみ波の伝播速度（ｍ／ｓ）を表し、Ｌは試験片の長さ（ｍ）を表し、Ｅは応力棒の縦弾性係数（Ｎ／ｍ２）を表し、ρは応力棒の密度（ｋｇ／ｍ３）を表す）
【００８０】
また、εｉ、εｒ及びεｔから、下記数式１０によって試験片７０のひずみεが算出される。
【００８１】
（数式１０）
0003910825.tif 000006
（数式１０において、Ｃ０は入力棒および出力棒中（応力棒）のひずみ波の伝播速度（ｍ／ｓ）を表し、Ｌは試験片の長さ（ｍ）を表し、Ｅは応力棒の縦弾性係数（Ｎ／ｍ２）を表し、ρは応力棒の密度（ｋｇ／ｍ３）を表す）
【００８２】
さらに、εｉ、εｒ及びεｔから、下記数式１１によって試験片７０の応力σが算出される。
（数式１１）
σ＝（Ｅ・Ａ／（２Ａｓ））・（εｉ＋εｒ＋εｔ）
＝（Ｅ・Ｄ２／（２（Ｄｓ）２））・（εｉ＋εｒ＋εｔ）−−−（１１）
（数式１１において、Ｅは入力棒および出力棒からなる応力棒の縦弾性係数（Ｎ／ｍ２）を表し、Ａは上記応力棒の断面積（ｍ２）を表し、Ａｓは試験片の断面積（ｍ２）を表し、Ｄは応力棒の直径（ｍ）を表し、Ｄｓは試験片の直径（ｍ）を表す）
【００８３】
こうして得られた試験片７０のひずみ時刻歴を、図１３のグラフに示す。図１３に示すように、曲線は、ピークＰ以降しばらくはなだらかであるが、その後、凹凸状となる。ピークＰ以降のなだらかな段階での点Ｓを選択し、この点Ｓにおける曲線に対する接線を画き、この接線と時間軸との交点から緩和時間λを導出し、下記数式１２によって求められる曲線を点Ｓ以降の曲線とすることによって、ひずみ時刻歴全体をなだらかな曲線（図１３中に点線で示す）とすることができる。これにより、最終的に得られる材料物性へのノイズの影響を除去することができる。
（数式１２）
ε（ｔ）＝ε０・ｅ−ｔ／λ      −−−（１２）
（数式１２において、ε０は接点におけるひずみを表す）
【００８４】
同様に、下記数式１３によって、応力時刻歴全体をなだらかな曲線とすることができ、これによって最終的に得られる材料物性へのノイズの影響を除去することができる。
（数式１３）
σ（ｔ）＝σ０・ｅ−ｔ／λ      −−−（１３）
（数式（１３）において、σ０は接点における応力を表す）
【００８５】
上記のようにして、試験片７０のひずみ時刻歴及び応力時刻歴については補正が行われている。
【００８６】
以上の方法により、スプリットホプキンソン棒試験機にて、高速大変形時のひずみの時刻歴データ、ひずみ速度の時刻歴データ、応力の時刻歴データを得ている。
【００８７】
以下、本発明の実施例及び比較例について詳述する。
まず、ウレタンゴムを主成分とする材料を用いてゴルフボールを作製した。ウレタンを主成分とする材料を１６０℃で３０分、直径４２．８ｍｍの金型で圧縮成型し、１ピースボールとした。
【００８８】
上記ウレタンゴムを主成分とする材料のひずみ、ひずみ速度、応力を上記したスプリットホプキンソン棒試験機により、打撃棒の衝突速度が７ｍ／ｓ、１４ｍ／ｓ、２０ｍ／ｓ、２５ｍ／ｓの４パターンについて測定した（室温２３℃、相対湿度５０％）。
測定時の最大ひずみと最大ひずみ速度を以下に示す。
衝突速度７ｍ／ｓ（最大ひずみ０．１２、最大ひずみ速度１３７８／ｓ）
衝突速度１４ｍ／ｓ（最大ひずみ０．２４、最大ひずみ速度２７０３／ｓ）
衝突速度２０ｍ／ｓ（最大ひずみ０．３５、最大ひずみ速度３８９８／ｓ）
衝突速度２５ｍ／ｓ（最大ひずみ０．４３、最大ひずみ速度４７１６／ｓ）
４パターンの衝突速度で測定しているため、各衝突速度における位相角δの値を下記の表１に示す。
【００８９】
【表１】
0003910825.tif 000007
【００９０】
（実施例）
ウレタンゴムを主成分とする材料を試験片としてスプリットホプキンソン棒試験機により測定した、ひずみ、ひずみ速度、応力の各時刻歴データと、第１実施形態と同様の粘弾性モデルとを用い、粘性抵抗を考慮してシミュレーションを行った。上記第１実施形態のシミュレーション方法において、さらに、上記製品モデルに上記ひずみ、ひずみ速度、縦弾性係数の関係を入力することで、上記ひずみ、ひずみ速度の違いによる縦弾性係数の変化を考慮した。使用ソフトは、日本総合研究所（株）製、汎用ソフトＬｓ−ｄｙｎａ９４０とした。
実施例のシミュレーションにより予測された位相角δを上記表１に示す。
【００９１】
（比較例）
上記実施例と同様のウレタンゴムを主成分とする材料を試験片としてスプリットホプキンソン棒試験機により測定した、ひずみ、ひずみ速度、応力の各時刻歴データを用いた。縦弾性係数は可変とし、負荷状態時のひずみ、ひずみ速度と負荷の粘性抵抗を利用したシミュレーションを行った。上記した点以外は実施例と同様とした。
比較例のシミュレーションにより予測された位相角δを上記表１に示す。
【００９２】
有限要素法解析により、上記材料ウレタンゴムを主成分とする材料よりなるゴルフボールに、２００ｇ（ヘッドの重さと同一）のアルミニウム製の中空棒モデルを初速度３５、４０、４５ｍ／ｓの３通りの速度で衝突させた時のゴルフボールの性能、材料の変形状態をシミュレートし、上記した方法でゴルフボールの解析での反発係数を算出した。ここで得られた上記ウレタンゴムを主成分とする材料よりなるゴルフボールの実施例及び比較例の解析での反発係数を、それぞれ下記の表２に示す。
【００９３】
【表２】
0003910825.tif 000008
【００９４】
また、上記ウレタンゴムを主成分とする材料で成型されたゴルフボールの実物を使用した実験により、下記の方法でゴルフボールの実物での反発係数の値を測定した。
上記実施例及び比較例の解析での反発係数と下記の実物を使用した実験での反発係数との差（％）を上記した表２に示す。
【００９５】
（ゴルフボールの実物を用いた実験による反発係数の測定）
ゴルフボールの反発係数を測定する方法として、室温２３℃の条件下、上記材料よりなるゴルフボールに、ゴルフクラブのヘッドの代用として、２００ｇ（ヘッドの重さと同一）のアルミニウム製の中空棒を速度３５、４０、４５ｍ／ｓの３通りの速度で衝突させた。衝突前後の上記中空棒及びゴルフボールの速度を測定し、それぞれの速度及び重量からゴルフボールの反発係数を算出した。
なお、アルミニウム製の中空棒の衝突面は、フラットであり、ゴルフボールとフラットな正面衝突とした。ゴルフクラブヘッドのように衝突面に角度がないため、衝突時にゴルフボールが回転しないので、ゴルフボールの反発性能のみを評価可能とした。
【００９６】
表１に示すように、実施例の解析での位相角δは、各衝突速度において、実験結果の位相角δとほぼ同一の値となった。一方、比較例の解析での位相角δは、負荷と負荷の曲線を用いてシミュレーションしているために、応力−ひずみ曲線においてエネルギーロスが実施例に比べて小さくなり、各衝突速度において比較例の位相角δの値は実施例より小さくなり、実験結果に比べその値も小さく、その差も大きかった。従って、実施例は精度良く実際の実験結果をシミュレーションできていることが確認できた。
【００９７】
また、表２に示すように、実施例の解析での反発係数の値は、実際にゴルフボールを試作して、実物により求めた実験での反発係数の値と、いずれの中空棒速度の場合においても、よく一致していた。実施例と実物実験結果との差は−４．８２％〜−５．９８％であり、実物の実験結果を正確にシミュレートしていることが確認できた。一方、比較例の解析での反発係数の値は、位相角δの値が小さいために、反発係数は実験値に対して大きくなり、実施例に比べてその差も大きかった。具体的には、比較例と実物実験結果との差は１６．１％〜２２．２％であり、実験結果と大きく異なっていた。
【００９８】
以上より、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで場合分けして、ひずみが同じ値でもひずみの負荷状態時とひずみの除荷あるいは復元状態時とでの粘性抵抗の違いを考慮した本発明のシミュレーションを行うことで、材料物性が非線形性を示す粘弾性材料からなる製品の実使用条件下での性能を、精度良くシミュレーションにより予測できることが確認できた。
【００９９】
【発明の効果】
以上の説明より明らかなように、本発明によれば、要素毎に生じるひずみ、ひずみ速度の偏差成分において粘弾性特性を考慮するために、上記ひずみ、ひずみ速度を体積成分と偏差成分に分解し、主ひずみ座標系及び主ひずみ速度座標系に変換した偏差成分の主ひずみ及び主ひずみ速度を用いて主ひずみ座標系及び主ひずみ速度座標系の各軸毎に粘性抵抗を決定し、粘性抵抗を要素毎に各軸方向３成分について可変にしている。このため、３次元方向について正確に粘弾性材料の粘弾性を考慮した解析を行うことができ、シミュレーションの精度を向上することができる。
【０１００】
また、粘弾性材料の粘性を考慮した粘弾性モデルを用い、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで場合分けして粘性抵抗を導出し、ひずみ、ひずみ速度、ひずみの負荷状態時とひずみの除荷あるいは復元状態時とで異なった粘性抵抗の関係を製品モデルに入力しシミュレーションを行っているため、実際の材料状態を正確に捉えることができ、粘弾性材料の非線形性を示す物性値が、材料の変形の速度、大きさによって非線形的に変化する現象を精度良く表現することができる。
【０１０１】
さらに、シミュレーションに用いられるひずみ、ひずみ速度、応力の値は、粘弾性材料からなる製品の実使用状態を想定した条件下で測定されているため、実際の粘弾性材料の種々の変形状態に対応したシミュレーションを行うことができる。
【０１０２】
従って、材料の変形状態により、ひずみ、ひずみ速度の関係が変化し、損失係数等の材料物性が非線形性を示すような粘弾性材料からなる製品においても、シミュレーションにより性能や製品の動的挙動を精度良く解析することができる。
【０１０３】
実際にゴルフクラブヘッドで、ゴルフボールが打撃されたときと同等のひずみ、ひずみ速度の条件下で、実際の打撃試験を正確にシミュレートすることができ、実際の打撃に近い状態でのゴルフボールの性能や変形挙動を精度良く予測することができる。これにより、ゴルフボールの性能の優劣を左右する材料物性の把握が容易となり、製品の性能向上に役立つだけでなく、ゴルフボールの設計段階において、実際のボール試作回数を減らし、試作に要する費用と時間を削減することができる。
【図面の簡単な説明】
【図１】  ひずみεの時刻歴データを示す図である。
【図２】  ひずみ速度ε’の時刻歴データを示す図である。
【図３】  応力σの時刻歴データを示す図である。
【図４】  ひずみ−応力線図と縦弾性係数の算出方法を示す図である。
【図５】  本実施形態の粘弾性モデルとして用いた２要素のVoigtモデルを示す図である。
【図６】  粘性抵抗の時刻歴データを示す図である。
【図７】  ひずみと粘性抵抗の関係を示す図である。
【図８】  ひずみとひずみ速度の関係を示す図である。
【図９】  ゴルフボールモデルのメッシュによる分割状況を示す図である。
【図１０】  アルミニウム製中空棒モデルのゴルフボールモデルへの衝突状況を示し、（Ａ）は衝突前、（Ｂ）は衝突中、（Ｃ）は衝突後の図である。
【図１１】  損失係数の算出方法を示す図である。
【図１２】  スプリットホプキンソン棒試験機が示された模式的正面図である。
【図１３】  試験片のひずみ時刻歴の状態が示された図である。
【符号の説明】
１０  ゴルフクラブモデル
１１  要素
１２  節点
２０  アルミニウム製中空棒モデル
【図１】
0003910825.tif 000009
【図２】
0003910825.tif 000010
【図３】
0003910825.tif 000011
【図４】
0003910825.tif 000012
【図５】
0003910825.tif 000013
【図６】
0003910825.tif 000014
【図７】
0003910825.tif 000015
【図８】
0003910825.tif 000016
【図９】
0003910825.tif 000017
【図１０】
0003910825.tif 000018
【図１１】
0003910825.tif 000019
【図１２】
0003910825.tif 000020
【図１３】
0003910825.tif 000021