(57)【要約】  【課題】ゴルフクラブシャフトの硬さの指標とされる、ゴルフクラブシャフトのグリップ部を固定端として振動させた時の振動数を用いて、ゴルフスウィング中のヘッドスピードを最大化させるゴルフクラブシャフトの振動数（最適振動数）を決定する。  【解決手段】ゴルファが行うゴルフスウィング中の、アドレスからバックスウィングを開始してトップの状態にいたる第１の時間Ｔ１（秒）と、トップの状態からゴルフボールの打撃にいたる第２の時間Ｔ２（秒）とを検出し、検出された第1の時間Ｔ１の第２の時間Ｔ２に対する比率であるスウィングリズム係数Ｒを用いて、下記式（１）によりゴルフクラブシャフトの最適振動数Ｆ（ｃｐｓ）を算出する。      Ｆ  ＝  ａ・Ｒ  ＋  ｂ／Ｔ２    （１）    （ａ，ｂは定数）【選択図】図１
2005034614.tif 000002
【特許請求の範囲】
【請求項１】
ゴルフクラブシャフトのグリップ側を固定端とし、このグリップ側と反対側を振動させたときのゴルフクラブシャフトの振動数であって、ゴルフクラブを把持して行うゴルフスウィングに対して、ゴルフボールの打撃時のヘッドスピードを略最大化させるときのゴルフクラブシャフトの特性を特徴付ける最適振動数を決定する決定方法であって、
ゴルファが行うゴルフスウィング中の、バックスウィングを開始してトップの状態にいたる第１の時間と、トップの状態からダウンスウィングを開始してゴルフボールの打撃にいたる第２の時間とを検出するステップと、
検出された前記第１および第２の時間の比率であるスウィングリズム係数を用いることにより、前記最適振動数を算出するステップと、を有することを特徴とするゴルフクラブシャフトの最適振動数の決定方法。
【請求項２】
前記最適振動数は、前記スウィングリズム係数の他に、前記第２の時間を用いて算出される請求項１に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項３】
前記最適振動数をＦ（ｃｐｓ）、前記スウィングリズム係数をＲ、前記第２の時間をＴ2（秒）とすると、前記最適振動数Ｆは下記式（１）で表される請求項２に記載のゴルフクラブシャフトの最適振動数の決定方法。
Ｆ  ＝  ａ・Ｒ  ＋  ｂ／Ｔ2                        （１）
但し、ａ、ｂは定数。
【請求項４】
前記定数ａおよびｂは、ゴルフクラブシャフトの長さに応じて変化するものである請求項３に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項５】
前記定数ａおよびｂの前記定数ａは０．６９〜０．９７の範囲にあり、かつ、前記定数ｂは０．３８〜０．７３の範囲にある請求項３または４に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項６】
前記定数ａおよびｂの前記定数ａは０．７０〜０．８４の範囲にあり、かつ、前記定数ｂは０．８４〜１．０２の範囲にある請求項３または４に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項７】
前記定数ａおよびｂの前記定数ａは１．０８〜１．２２の範囲にあり、かつ、前記定数ｂは−０．１８〜０．００の範囲にある請求項３または４に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項８】
前記最適振動数を算出するステップにおいて、前記定数ａおよびｂの値の組み合わせを複数設定して、前記定数ａおよびｂの組み合わせそれぞれに応じた前記最適振動数の候補値を複数求めることを特徴とする請求項３または４に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項９】
前記複数の定数ａおよびｂの組み合わせは、前記定数ａが０．６９〜０．９７の範囲かつ前記定数ｂが０．３８〜０．７３の範囲にある定数ａおよびｂの組み合わせ、および前記定数ａが０．７０〜０．８４の範囲かつ前記定数ｂが０．８４〜１．０２の範囲にある定数ａおよびｂの組み合わせ、および前記定数ａが１．０８〜１．２２の範囲かつ前記定数ｂが−０．１８〜０．００の範囲にある定数ａおよびｂの組み合わせ、の少なくとも２つの組み合わせであることを特徴とする請求項８に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項１０】
前記複数の最適振動数の候補値のうち、所定の数値範囲にある候補値を前記最適振動数の値として選択することを特徴とする請求項８または９に記載のゴルフクラブシャフトの最適振動数の決定方法。
【請求項１１】
請求項１〜１０のいずれか１項に記載のゴルフクラブシャフトの最適振動数の決定方法で算出された最適振動数を有するゴルフクラブをゴルファに最適なゴルフクラブとして選定することを特徴とするゴルフクラブの選定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ゴルファがゴルフクラブを把持して行うゴルフスウィングに対してゴルフボールの打撃時のヘッドスピードを略最大化させる最適なゴルフクラブシャフトの硬さを表す振動数を決定するゴルフクラブシャフトの最適振動数の決定方法および、この方法を用いたゴルフクラブの選定方法に関する。
【背景技術】
【０００２】
近年、ゴルフ用具の市場において、ゴルファ自身のゴルフスウィングの特徴に適したゴルフクラブヘッドがゴルファによって選択され購入されている。また、ゴルフ販売業者、製造業者等は、様々なゴルフスウィングに適する種々のゴルフクラブを提案し、ゴルファ各自のゴルフスウィングの特徴に適したゴルフクラブを選択させて販売させる方法を提案している。例えば、ゴルフボールの打ち出し時点のゴルフクラブヘッドの移動速度であるヘッドスピードの大小により、ゴルフクラブシャフトの材質を変え、また、ゴルフクラブヘッドの材質や形状を変えたゴルフクラブ等が提案されている。
一方、上記ゴルファに適したゴルフクラブの選択、販売をより効果的に行うために、実際に様々なゴルフクラブをゴルファに試打してもらいゴルフボールの初期弾道を計測し、好ましい初期弾道となるゴルフクラブを推奨することも提案されている。
【０００３】
このような初期弾道を計測する装置は、打ち出し直後のゴルフボールを撮影し、ゴルフボールの初速、スピン量等を計測するため、高速度カメラ等の高価な装置を必要とする。また、ゴルフクラブを購入しようとするゴルファが多数のゴルフクラブを試打する必要があり、試打したゴルフクラブの範囲内からしかゴルフクラブを選択することができない、といった不都合もある。
【０００４】
これに対して、下記特許文献１には、ゴルフスウィングのスタートからゴルフボールの打ち出しまでのスウィング時間を検出し、この検出したスウィング時間に基づいてゴルファに最適なフレックスのゴルフクラブシャフトを選定する選定装置が提案されている。
この選定装置では、ゴルフスウィングにおいて最もタイミングのとり易い、あるいは打ち易いとゴルファが感じるゴルフクラブシャフトのフレックスを選ぶために、スウィング時間とゴルフクラブシャフトの最適振動数との対応関係を関数で予め定めておき、これを利用する。
【特許文献１】特許第３０６１６４０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ところで、ゴルフボールを打ち出す時のゴルフクラブヘッドのヘッドスピードは、ゴルフクラブシャフトの種類によって種々変化することが知られている。これは、バックスウィングからトップの状態さらにはダウンスウィングに至る一連のゴルフスウィングがグリップ部を介してゴルフクラブに入力されるため、撓んだ状態にあるゴルフクラブシャフトをトップの状態から振り下ろす間（ダウンスウィング中）、ゴルフクラブヘッドに作用する遠心力によってゴルフクラブシャフトは過渡的に変形し、この時の変形速度がゴルフボールの打ち出し時点のゴルフクラブヘッドのヘッドスピードに加算されるからである。このため、大部分のゴルファが最も気にするゴルフボールの飛距離に大きな影響を与えるヘッドスピードを最大化させるためには、ゴルフクラブシャフトを最適に選定することは極めて重要である。
【０００６】
一方、ゴルフクラブシャフトは「Ｒ」、「Ｓ」、「Ｘ」等のような硬さ表示によってゴルフクラブシャフトのフレックスは区分けされている。また、これらの硬さは、ゴルフクラブシャフトを片持梁のようにゴルフクラブシャフトのグリップ部を固定端として反対側を振動させたときに計測される振動数、いわゆるゴルフクラブシャフトの硬さを表す振動数によって分類することもできる。しかし、これらのゴルフクラブシャフトは、ヘッドスピードの高いゴルファ程硬いゴルフクラブシャフトを用いればよいといったことが知られているだけで、ゴルフクラブシャフトの選定方法が定量的に詳細に知られているわけではない。
【０００７】
上記特許文献１に記載の選定装置は、ゴルフスウィング中のバックスウィングの開始時点からトップの状態さらにはダウンスウィング、そしてゴルフボールの打ち出しに至るスウィング時間を検出し、このスウィング時間に基づいてゴルファに最適なフレックスのゴルフクラブシャフトを選定するものであるが、この最適なフレックスのゴルフクラブシャフトは、ゴルフスウィングにおいて最もタイミングのとり易い、あるいは打ち易いとゴルファが感じるものに過ぎない。このため、ゴルフボールの飛距離に大きな影響を与えるヘッドスピードを最大化させることはできない。
【０００８】
そこで、本発明は、上記従来技術の問題点を解消するために、ゴルフクラブシャフトのグリップ側を固定端として反対側を振動させたときの振動数、すなわち、従来からゴルフクラブシャフトの硬さの指標として用いられているゴルフクラブシャフトの振動数であって、ゴルファがゴルフクラブを把持して行うゴルフスウィングに対して、ゴルフボールの打撃時のヘッドスピードを略最大化させるときのゴルフクラブシャフトの特性を特徴付ける最適振動数を決定する決定方法およびこれを用いたゴルフクラブの選定方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
そこで、本発明は、ゴルフクラブシャフトのグリップ側を固定端として反対側を振動させたときのゴルフクラブシャフトの振動数であって、ゴルフクラブを把持して行うゴルフスウィングに対して、ゴルフボールの打撃時のヘッドスピードを略最大化させるときのゴルフクラブシャフトの特性を特徴付ける最適振動数を決定する決定方法であって、ゴルファが行うゴルフスウィング中の、バックスウィングを開始してトップの状態にいたる第１の時間と、トップの状態からダウンスウィングを開始してゴルフボールの打撃にいたる第２の時間とを検出するステップと、検出された前記第１および第２の時間の比率であるスウィングリズム係数を用いることにより、前記最適振動数を算出するステップと、を有することを特徴とするゴルフクラブシャフトの最適振動数の決定方法を提供する。
【００１０】
ここで、前記最適振動数は、前記スウィングリズム係数の他に、前記第２の時間を用いて算出されるのが好ましい。
【００１１】
前記最適振動数をＦ（ｃｐｓ）、前記スウィングリズム係数をＲ、前記第２の時間をＴ2（秒）とすると、前記最適振動数Ｆは、例えば、下記式（１）で表される。
Ｆ  ＝  ａ・Ｒ  ＋  ｂ／Ｔ２                        （１）
但し、ａ、ｂは定数。
【００１２】
その際、前記定数ａおよびｂは、ゴルフクラブシャフトの長さに応じて変化するものであり、前記定数ａは０．６９〜０．９７の範囲にあり、定数ｂは０．３８〜０．７３の範囲にあることが好ましい。
【００１３】
前記定数ａおよびｂの前記定数ａは０．７０〜０．８４の範囲にあり、かつ、前記定数ｂは０．８４〜１．０２の範囲にあることもまた好ましい。
【００１４】
また、前記定数ａおよびｂの前記定数ａは１．０８〜１．２２の範囲にあり、かつ、前記定数ｂは−０．１８〜０．００の範囲にあることもまた好ましい。
【００１５】
また、前記最適振動数を算出するステップにおいて、前記定数ａおよびｂの組み合わせを複数設定して前記最適振動数の値の候補値を複数求め、求められた複数の候補値の中から１つの候補値を前記最適振動数の値として選択することが好ましい。
【００１６】
この際、前記複数の最適振動数の候補値のうち、所定の数値範囲にある候補値を前記最適振動数の値として選択することが好ましい。
【００１７】
また本発明は、前記ゴルフクラブシャフトの最適振動数の決定方法で算出された最適振動数を有するゴルフクラブをゴルファに最適なゴルフクラブとして選定することを特徴とするゴルフクラブの選定方法を提供する。
【発明の効果】
【００１８】
本発明によれば、バックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２の比率であるスウィングリズム係数を用いることで、従来より計測されてきたゴルフクラブシャフトの硬さの指標である振動数に関し、ゴルフスウィング中のヘッドスピードを最大化させるようにゴルフクラブシャフトの最適振動数を決定することができる。すなわち、ゴルフスウィングの可能性を最大限に引き出すようにゴルフクラブシャフトの最適振動数を決定することができる。この際、最適振動数を算出する式の係数を複数設定して、最適振動数の候補値を複数求め、これら最適振動数の候補値の中から最適振動数を絞り込むことで、初級者からプロに至る幅広い層のゴルファについて適正な最適振動数を決定することができる。さらに、最適振動数を算出する式の係数は、ゴルフクラブシャフトによって変化するので、ゴルフクラブ製造業者や販売業者は、ゴルフクラブシャフトの長さに応じて、最適振動数を詳細に設定することができる。これにより、ゴルフクラブのきめ細かい商品開発に有効に利用することができる。また、この最適振動数をゴルファが知ることで、ゴルフクラブシャフトが自分に適したものであるか否かを容易に判断することもできる。
【発明を実施するための最良の形態】
【００１９】
本発明のゴルフクラブシャフトの最適振動数の決定方法およびゴルフクラブの選定方法を添付の図面に示す好適実施形態に基づいて、以下に詳細に説明する。
図１は、ゴルフクラブシャフトの最適振動数を決定し、複数のゴルフクラブの中から最適なゴルフクラブを選定する装置１０の概略構成図である。図１では、ゴルファがゴルフクラブ３０を把持してゴルフスウィングする状態を示している。
【００２０】
装置１０は、計測装置２０と解析装置２４とを有する。解析装置２４において、計測装置２０から供給される３次元時系列データを用いて、従来よりゴルフクラブシャフトの硬さを特徴付ける指標の１つとされる、グリップ側を固定端としたときに発生するゴルフクラブシャフトの振動の振動数を最適に決定する装置である。ここで、３次元時系列データとはゴルフクラブのグリップ部３６の動きをグリップ部３６の位置と向きによって表したデータである。
【００２１】
計測装置２０は、図１に示すように、ゴルフスウィングするゴルファの背後に配置固定したトランスミッタ２０ａから３種類の所定の磁場を次々に発生させる。一方、移動および回転するグリップ部３６に固定されたレシーバ２０ｂが、トランスミッタ２０ａによって作られる３種類の磁場内の位置および向きに対応して磁気を感知して合計９つの出力電圧を出力する。この出力電圧からコントローラ２０ｃにおいてデータ処理がなされてレシーバ２０ｂの３次元位置と向き（オイラー角）のデータを得ることができる。
【００２２】
このような計測装置２０をより詳しく説明すると、図２に示されるように、所定の磁場を形成するトランスミッタ２０ａと、この磁場の強さや方向に応じて３軸方向の出力電圧を発生するレシーバ２０ｂと、コントローラ２０ｃとを有して構成される。
コントローラ２０ｃは、トランスミッタ２０ａに所定の３種類の磁場を順次発生する駆動信号を生成するとともに、レシーバ２０ｂから出力された信号を検出し、検出された信号よりデータ処理を行って、時系列データを演算して算出する。時系列データとは、所定の位置、例えばトランスミッタ２０ａの位置を基準位置とし、お互いに直交する３方向を基準方向とする３次元座標系における位置座標（ｘｍ，ｙｍ，ｚｍ）と、所定の基準方向、例えばＸ，Ｙ，Ｚ座標軸方向に対するレシーバ２０ｂの向きを表す姿勢角度、すなわちヨー角、ピッチ角およびロール角（以降では、オイラー角（θｙ，θｐ，θｒ）と表す）の時系列データである。
【００２３】
トランスミッタ２０ａおよびレシーバ２０ｂは、図２に示されるように、お互いに直交する３軸方向に各々ループ状に巻かれた３つのコイルによって構成され、トランスミッタ２０ａは、ゴルフスウィングするゴルファの後方に固定配置され、レシーバ２０ｂは、ゴルフクラブ３０のグリップ部３６の端部に固定される。
コントローラ２０ｃは、３種類の磁場を順次発生する駆動信号を生成する駆動回路２０ｄと、レシーバ２０ｂからの出力された信号を検出する検出回路２０ｅと、駆動回路２０ｄおよび検出回路２０ｅの制御を行い、検出回路２０ｅから送られる信号よりレシーバ２０ｂの３次元位置とオイラー角を求める制御ユニット２０ｆによって構成される。トランスミッタ２０ａは駆動回路２０ｄに、レシーバ２０ｂは検出回路２０ｅに、各々接続される。
計測装置２０は、以上のように構成される。
【００２４】
なお、グリップ部３６に固定されたレシーバ２０ｂの基準位置に対する３次元位置座標（ｘｍ，ｙｍ，ｚｍ）と基準方向に対するオイラー角（θｙ，θｐ，θｒ）の時系列データは以下のようにして得られる。
図２に示すように、駆動回路２０ｄは、制御ユニット２０ｆの指令信号にしたがって、周波数と位相が常時一定の同一信号を出力し、トランスミッタ２０ａの３軸方向に巻かれた３つのループ状コイルを順次励磁する。各ループ状コイルは、励磁のたびに各々異なる磁場を発生し、この磁場に基づいてレシーバ２０ｂの３軸方向に巻かれた３つのループ状コイルに各々独立な電圧を発生させる。この電圧は、トランスミッタ２０ａの３つのループ状コイルによって励磁される３つの磁場それぞれに応じて、レシーバ２０ｂの３つのループ状コイルに発生する３つの独立した電圧であるため、合計９個（３×３個）の電圧が得られる。
【００２５】
一方、磁場を形成させるトランスミッタ２０ａが所定の位置に固定設置されているので、発生する磁場の強さと方向に関する分布はトランスミッタ２０ａの設置された基準位置および、基準方向に対して既知となる。そこで、この形成された磁場によって生じる９つの電圧を用いることによって、上記基準位置に対するレシーバ２０ｂの３次元位置座標（ｘｍ，ｙｍ，ｚｍ）と上記基準方向に対するオイラー角（θｙ，θｐ，θｒ）の６つの未知数を求めることができる。
コントローラ２０ｃの制御ユニット２０ｆは、検出回路２０ｅから送られてきた９つの電圧を用いて、３次元位置座標（ｘｍ，ｙｍ，ｚｍ）とオイラー角（θｙ，θｐ，θｒ）のデータを演算して求める。
求められた３次元位置座標（ｘｍ，ｙｍ，ｚｍ）とオイラー角（θｙ，θｐ，θｒ）は、解析部２４に供給される。
このような計測装置２０として、例えば、3SPACE  FASTRAK（Polhemus社製）を挙げることができる。
【００２６】
解析装置２４は、計測装置２０から供給されたデータに基づいて、３次元時系列データを抽出するデータ抽出部２４ａと、３次元時系列データを用いて、バックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２を検出する時間検出部２４ｂと、検出したバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２を用いて最適振動数を算出する振動数算出部２４ｃと、算出した最適振動数に基づいて最適なゴルフクラブを選定する選定部２４ｄとを有する。解析装置２４には、モニタ２５および操作系（マウス、キーボード）２６が接続されている。
解析装置２４の上記各部分はモジュール化されたサブルーチンプログラムを実行することで機能が発揮するように構成されたコンピュータである。しかし、本発明においては、上記各部分が専用回路で組まれた専用装置であってもよい。
【００２７】
図１に示される解析装置２４には、解析部２４で行われる各処理が示されている。
まず、データ抽出部２４ａでは、前述した計測装置２０から供給された３次元位置座標（ｘｍ，ｙｍ，ｚｍ）とオイラー角（θｙ，θｐ，θｒ）の時系列データから、グリップ部３６の所定の位置を基準とした３次元位置と直交座標系（ＸＹＺ座標系）におけるグリップ部３６の３次元方向における向きについての３次元時系列データを演算して求める（ステップＳ１０）。例えば、３次元位置座標（ｘｍ，ｙｍ，ｚｍ）と、直交座標系における方位角と仰角とゴルフクラブシャフト３２の軸周りの回転の３次元時系列データを求める。グリップ部３６の３次元位置とこの直交座標系におけるグリップ部３６の３次元方向の向きの時系列データは、解析装置２４の図示されないメモリに記憶される。
さらに、算出されたグリップ部３６の３次元位置と３次元方向の向きからなる３次元時系列データは、設定された所定の方向から見たグリップ部３６の軌道、例えば、ゴルファと正対する正面からみたゴルフスウィング中のグリップ部３６の軌道のデータに変換されて、モニタ２５に供給される。また、３次元時系列データがモニタ２５にグラフとして表示される。
【００２８】
時間検出部２４ｂは、求められた３次元時系列データから、図３に示すようなバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２を検出する（ステップＳ２０）。バックスウィングの開始時点、トップの状態は、いずれもゴルフクラブの挙動は停止していることから、３次元時系列データの波形から、極値となっている箇所を読み取ることで特定することができる。また、ゴルフボールを打ち出すインパクト時点も例えば、ゴルフボールの打ち出しによる時系列データの変動から特定することができる。
【００２９】
図４（ａ）〜（ｃ）は、モニタ２５に表示されるＸＹＺ座標系で表されたグリップ部３６の挙動を表した３次元時系列データのグラフの一例である。
図４（ａ）〜（ｃ）からバックスウィング時間Ｔ１は約１．１（秒）、ダウンスウィング時間Ｔ２は約０．３５（秒）である。
【００３０】
振動数算出部２４ｃは、検出されたバックスウィング時間Ｔ１（秒）およびダウンスウィング時間Ｔ２（秒）から、ダウンスウィング時間Ｔ２に対するバックスウィング時間Ｔ１の比率であるスウィングリズム係数Ｒ（Ｔ１／Ｔ２）およびダウンスウィング時間Ｔ２を用いて最適振動数Ｆ（ｃｐｓ）を上記式（１）に従って算出する（ステップＳ３０）。ここで、振動数（ｃｐｓ）とは、ゴルフクラブシャフトのグリップ部をチャックにて固定把持し、ゴルフクラブヘッド側の先端部を自由端として振動させ、振動によってゴルフクラブシャフトの先端部分が光電管を横切る回数を計数し、１秒当たりの横切る回数から振動数を求めたものに相当し、最適振動数とは、ゴルフクラブシャフトがゴルフボールの打ち出し時のヘッドスピードを略最大化させるときの上記振動数である。あるいは、振動数（ｃｐｓ）とは、ゴルフクラブシャフト単体のグリップの装着される後端部を固定把持し、ゴルフクラブヘッドの装着される先端部に加速度ピックアップを貼り付けて、この加速度ピックアップで得られる加速度信号から求められる固有振動数に相当し、最適振動数とは、ゴルフクラブシャフトがゴルフボールの打ち出し時のヘッドスピードを略最大化させるときの上記固有振動数でもある。
【００３１】
ここで、上記式（１）における定数ａ，ｂは、ゴルフクラブの番手およびゴルフクラブシャフトの長さあるいはゴルファのスウィングのタイプに応じて変化するものである。
ドライバー（＃１）のゴルフクラブの場合であって、例えばゴルフクラブシャフトの長さが４５インチより長く４７インチ以下の長尺ゴルフクラブの場合、定数ａは０．７６を中心として±０．０７の範囲（０．６９〜０．８３）あり、定数ｂは、０．６４を中心として±０．０９の範囲（０．５５〜０．７３）にあることが好ましい。また、例えばゴルフクラブシャフトの長さが４３インチ以上４５インチ以下のゴルフクラブの場合、定数ａは０．９０を中心として±０．０７の範囲（０．８３〜０．９７）にあり、定数ｂは、０．４７を中心として±０．０９の範囲（０．３８〜０．５６）にあることが好ましい。
本発明で算出される最適振動数は、必ずしも１つの振動数の値によって定める必要はなく、定数ａおよびｂが上述のような範囲を持っていることから、最適振動数を、定数ａ、ｂのとり得る範囲によってとり得る振動数の範囲によって定めてもよい。この場合、例えばヘッドスピード別に最適振動数の範囲を定めてもよい。
【００３２】
選定部２４ｄでは、算出された最適振動数を有するゴルフクラブをゴルファに最適なゴルフクラブとして選定する（ステップＳ４０）。具体的には、選定部２４ｄに設けられた図示されないデータベースに、最適振動数とこの最適振動数を有するゴルフクラブシャフトの製品名や型式、さらには、このゴルフクラブシャフトを備えたゴルフクラブの型式や製品名をゴルフクラブシャフトと対応させて予め記憶させておく。そして、ステップＳ３０で算出された最適振動数からこの最適振動数に対応するゴルフクラブシャフトを備えるゴルフクラブの製品名、型式等の情報を検索して、検索結果をモニタ２５に表示させる。
このようにして、計測装置２０を用いて計測されたゴルファのゴルフクラブスウィングに最適なゴルフクラブ、すなわちゴルフボールの打ち出し時点のゴルフクラブヘッドのヘッドスピードを略最大化するゴルフクラブをゴルファに提供することができる。
【００３３】
本発明は、ゴルフスウィングの挙動を計測して得られるバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２の比率であるスウィングリズム係数Ｒ、さらにダウンスウィング時間Ｔ２を用いて、ゴルフスウィングに対してゴルフボールの打ち出し時点のヘッドスピードが略最大となるようなゴルフクラブシャフトの最適振動数を決定することを特徴とする。
このように、スウィングリズム係数Ｒおよびダウンスウィング時間Ｔ２を用いて最適振動数Ｆを決定することができるのは、ゴルフクラブヘッドの挙動を正確に再現することのできる有限要素法を利用したゴルフクラブのシミュレーションの解析結果によるものである。具体的には、ゴルフクラブの設計変更に応じて有限要素モデルを種々変更し、しかも、この設計変更を、ゴルフクラブヘッドのヘッドスピードが最大となるように目的関数を設定して、ゴルフクラブのシミュレーションの解析からゴルフクラブの最適設計仕様を探索するとともに、求められた最適設計仕様のゴルフクラブにおける振動数とバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２の関係を本願発明者が初めて見出したことによるものである。
【００３４】
具体的には、６人のゴルファに計測装置２０を用いてゴルフクラブシャフト長さが４６インチの＃１のウッド系のゴルフクラブ（ドライバー）でゴルフスウィングをしてもらい、この６人のゴルフスウィングのそれぞれに対して、ヘッドスピードが最大となる最適設計仕様のゴルフクラブモデルをそれぞれ探索し、このときのゴルフクラブモデルのグリップ側を固定し、その反対側を振動させたときの固有振動数を最適振動数Ｆとして、固有値解析によって算出した。なお、この固有値解析によって求まる固有振動数における振動形態は、後述するように図１１に示すゴルフスウィングのスウィングシミュレーションにおけるゴルフクラブの変形形態に略一致している。上記固有振動数はゴルフクラブシャフトの硬さを表す振動数に対応するものであるといえ、ヘッドスピードが最大となるときの上記固有振動数は最適振動数ということができる。
図５は、その時の結果である。図５中、横軸ｘにはバックスウィング時間Ｔ１（秒）、縦軸ｙには最適振動数Ｆ（ｃｐｓ）×ダウンスウィング時間Ｔ２（秒）をとっている。これによると、相関係数Ｒ２＝０．９６１５となって、ｙ＝０．７６８３・ｘ＋０．６４４２の直線に極めて良好に回帰していることがわかる。
これより、最適振動数Ｆとバックスウィング時間Ｔ１（秒）、ダウンスウィング時間Ｔ２との関係は、Ｆ＝０．７６８３・（Ｔ１／Ｔ２）＋０．６４４２／Ｔ２と表すことができる。
【００３５】
これに対し、図６は、上記特許文献１に開示される選定装置にて用いられる１分間の最適振動数Ｆ（ｃｐｍ）（＝Ｆ（ｃｐｓ）×６０）とスウィング時間Ｔ（＝バックスウィング時間Ｔ１＋ダウンスウィング時間Ｔ２）との関係を実線で示したものである。また、上述した６人のゴルフスウィングのスウィング時間と選定装置で決定される最適振動数（ｃｐｍ）との関係をプロットしている。これによると、バックスウィング時間Ｔ１＝１．２９５（秒）、ダウンスウィング時間Ｔ２（秒）＝０．４０５（秒）のゴルファＡと、バックスウィング時間Ｔ１＝１．２５０（秒）、ダウンスウィング時間Ｔ２秒）＝０．４５０（秒）のゴルファＢとは、スウィング時間Ｔ（＝Ｔ１＋Ｔ２）が同じであるものの最適振動数Ｆ（ｃｐｍ）が異なったものである。しかも、ゴルファＢは上記実線から大きく外れていることがわかる。
【００３６】
このように、本発明では、バックスウィング時間Ｔ１とダウンスウィングＴ２との比率を用いることで、ヘッドスピードを略最大化する最適振動数Ｆ（ｃｐｓ）を決定することができる。
ここで、ゴルフスウィングのリズムの速いといわれるゴルファ程スウィングリズム係数Ｒが大きく、その結果最適振動数Ｆ（ｃｐｓ）は高くなる。一方、ゴルフスウィングのリズムのゆったりとしているといわれるゴルファ程、スウィングリズム係数Ｒが小さく、その結果最適振動数Ｆ（ｃｐｓ）は低くなる。このことは、従来より定性的に言われてきた、ゴルフスウィングの速いゴルファは硬いゴルフクラブシャフトを、ゴルフスウィングのゆったりしているゴルファは軟らかいゴルフクラブシャフトを用いるとよい、といった事実と対応する。
【００３７】
このようなゴルフスウィングに対してゴルフボールの打ち出し時点のヘッドスピードを略最大化させるゴルフクラブシャフトを特徴付ける最適振動数Ｆ（ｃｐｓ）は、以下のようなゴルフクラブのゴルフスウィングシミュレーションの解析による最適設計仕様の探索結果により明らかとなったのである。
【００３８】
図７（ａ）〜（ｃ）は、ゴルフスウィングシミュレーションの解析に用いるゴルフクラブ３０の有限要素モデル化の説明図である。
図７（ａ）は、ゴルフクラブシャフト３２の一端にグリップ部３６を有するゴルフクラブ３０の一例を表しており、これに対して図７（ｂ）あるいは図７（ｃ）はゴルフクラブ３０のゴルフクラブモデルの一例を表している。
図７（ｂ）に示す例では、ゴルフクラブモデル４０は、各要素が梁要素のモデルで離散化された有限要素モデルであって、各要素の梁要素のモデルは断面積が一定の真直梁モデルである。ゴルフクラブシャフト３２に対応するゴルフクラブシャフトモデル４２は、要素１〜１２からなる真直梁モデルでモデル化され、グリップ側（バット側）と反対側（チップ側）の先端には、ゴルフクラブヘッド３４の質量をゴルフクラブヘッド３４の重心位置に対応する重心位置Ｇに集中質量を与えた質点モデル４４が付加され、さらに、バット側には、剛体要素である要素１３からなるグリップモデル４６が設けられている。質点モデル４４における重心位置Ｇおよび集中質量は、例えば、基準案におけるゴルフクラブヘッド３４の重心位置および質量が用いられる。
なお、図７（ｂ）および図７（ｃ）は、いずれもゴルフクラブヘッドを質点モデルで表したものであるが、ゴルフクラブヘッドの形状を忠実に再現した６面体ソリッド要素等の固体ソリッド要素モデルからなる有限要素モデルを作成してもよい。
【００３９】
なお、ゴルフクラブシャフトモデル４２における要素分割数は１２分割に限られず、５分割以上であればよいが、スウィングシミュレーションの解析において計算の効率化と計算精度の向上を同時に満足させる点から、１０〜１５分割が好ましく、より好ましくは１２分割であるのがよい。分割数が多いと設計変更する設計パラメータとしての変数が多くなり、計算を効率的に行うことはできない一方、分割数が少ないとヘッドスピードや振動数を精度良く再現することができない。
【００４０】
ここで、各要素１〜１２の材料定数は、各要素における真直梁モデルにおける曲げ剛性である。
例えば、補強層を巻き付けたゴルフクラブシャフトの積層構造について、構造力学および公知の古典積層理論が適用されて、ゴルフクラブシャフトの材料定数や寸法等の設計パラメータの値から、ヤング率および断面２次モーメントを算出して曲げ剛性ＥＩz （ヤング率Ｅ×断面２次モーメントＩz ）を中間パラメータとして算出する。より具体的には、補強層における補強繊維やマトリクス層の各種力学物性値（ヤング率、剪断弾性係数、ポアソン比）、補強層における補強繊維の配向角度、巻き付け位置、巻き付け枚数、マンドレルの内径および外径等の設計パラメータの値から、ヤング率および断面２次モーメントを算出し、真直梁モデルにおける曲げ剛性ＥＩz を中間パラメータとして算出する。
図７（ｃ）に示すゴルフクラブモデルは、図７（ｂ）に示すゴルフクラブモデルと異なり、各要素が梁要素で構成された要素１３、１４、１５からなるものである。
【００４１】
このようなゴルフクラブモデルのグリップモデル４６を表す要素１３や要素１３〜１５に、ゴルフスウィングを再現するために、図１に示す計測装置２０で計測されたゴルフスウィングのグリップ部の位置及び向きを表す３次元時系列データ（図４（ａ）〜（ｃ）参照）を、グリップモデル４６の境界条件として与えて、ゴルフクラブモデル４０のゴルフスウィングの挙動を算出する。この算出は、公知の有限要素ソルバーによって行われる。
こうして、ゴルフボールを打ち出す時のゴルフクラブヘッドのヘッドスピードを算出することができる。
以上のゴルフクラブモデルの作成とゴルフスウィングのシミュレーションの解析は、ゴルフクラブの設計変更に応じてゴルフクラブモデルを変更する度に行われ、ヘッドスピードが最大となるゴルフクラブの最適設計仕様を探索して求めることができる。この後、この最適設計仕様のゴルフクラブのゴルフクラブモデルのグリップモデル４６を表す要素１３や要素１３〜１５を固定して固有値解析を行い、従来より硬さの指標として計測されている振動数に対応する固有振動数を算出する。
【００４２】
図８は、ゴルフクラブモデルを用いてゴルフスウィングのシミュレーションの解析を行ったときの解析結果の一例を示す図である。より具体的には、ゴルフクラブモデルのグリップモデル４６に、計測された３次元時系列データを境界条件として入力した時のゴルフクラブモデルの変形の時間履歴を示した図である。図８からわかるように、ゴルフクラブシャフトモデルがゴルフスウィングにより変形していることがわかる。図９は、ゴルフクラブモデルを用いて固有値解析を行ったときの解析結果の一例を示す図であり、従来より硬さの指標として計測されている振動数に対応する固有振動数における振動形態を示している。
【００４３】
このような、ゴルフクラブヘッドのヘッドスピードを最大にするゴルフクラブシャフトの設計方法は、本願出願人が既に出願した特願平２００２−１６６８７２号に詳細に記載されている。
【００４４】
なお、上記装置１０では、計測装置２０を用いて得られた３次元時系列データからバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２を求めたが、これらの時間は計測装置２０によって求められるものに限定されない。例えば、高速度ビデオカメラや家庭用ビデオカメラで撮影された画像から、アドレスの状態、トップの状態およびゴルフボールの打ち出し（インパクト）の状態を表す画像の画像フレームを定め、アドレスの状態からトップの状態およびトップの状態からインパクトの状態にいたる画像フレーム数からバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２を求めてもよい。
本発明において最適振動数の算出に用いられる、ゴルフクラブにおける振動数（ゴルフクラブシャフトの硬さを表す振動数）とバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２の関係は、以上のようなゴルフクラブのゴルフスウィングシミュレーションの解析の結果によって本願発明者によって見出されたものである。
【００４５】
本発明では、前記定数ａおよびｂの値は１組の値に限定されず、より広範な様々なゴルファについてゴルフクラブシャフトの最適設計案を求める場合、ａおよびｂの値の組み合わせを複数設定して最適振動数の候補値を複数求め、複数の候補値から所定の条件を満たす適正な候補値を最適振動数の値として選択してもよい。
例えば、ヘッドスピードがおおよそ５０ｍ／ｓ以上のゴルフスウィング（上級者レベル、プロレベルのゴルファのゴルフスウィング）について、最適なゴルフクラブシャフトの硬さを表す振動数を求める場合、定数ａが０．６９．〜０．９７の範囲にあり、かつ定数ｂが０．３８〜０．７３の範囲にある定数ａおよびｂの組み合わせと、定数ａが０．７０〜０．８４の範囲にあり、かつ定数ｂが０．８４〜１．０２の範囲にある定数ａおよびｂの組み合わせと、定数ａが１．０８〜１．２２の範囲にあり、かつ定数ｂが−０．１８〜０．００の範囲にある定数ａおよびｂの組み合わせ、の３つの前記ａおよびｂの値の組み合わせを設定して最適振動数の候補値を複数求め、複数の候補値から所定の条件を満たす、例えば予め設定された振動数の範囲に含まれる適正な候補値を最適振動数として選択すればよい。この際、所定の条件を満たす適正な候補値は１つの振動数の値（候補値）である必要はなく、複数の候補値に絞られてもよい。
【００４６】
このように、バックスウィング時間Ｔ１とダウンスウィング時間Ｔ２との比率を用いた式（１）において、定数ａおよびｂの値の組み合わせを複数設定してゴルフスウィングに応じた最適振動数の候補値を複数求め、複数の候補値のなかから予め設定された数値範囲内にある候補値を最適振動数の値として選択することで、初級者、中級者、上級者やプロゴルファーなど幅広いゴルファについてゴルフクラブシャフトの最適振動数を求めることが出来る。
このように、定数ａおよびｂの値の組み合わせを複数設定することで幅広いゴルファに対して最適な振動数（ゴルフクラブシャフトの硬さを表す振動数）の値を求めることが出来るのは、上述のゴルフクラブのシミュレーション解析によるゴルフクラブの最適設計仕様の探索において、種々のゴルフクラブにおける振動数とヘッドスピードとの関係を基に、探索の際に算出した種々のゴルフクラブにおける振動数とバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２の関係を本願発明者が初めて見出したことによるものである。
【００４７】
具体的には、ゴルフクラブシャフト長が４６インチの所定のゴルフクラブを把持してゴルフスウィングを行った際の、Ｈ氏のゴルフスウィングとＮ氏のゴルフスウィングのそれぞれについて、上述のスウィング解析によるシミュレーション演算を行い、種々のゴルフクラブモデルについてゴルフクラブの振動数とヘッドスピードとの関係を求め、求められた関係に基づきヘッドスピードを略最大にする最適振動数の候補値（以降、最適振動数候補値とする）を各ゴルファ毎に複数求め、求められた最適振動数候補値とバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２の関係を見出したことによる。
ここで、Ｈ氏のゴルフスウィングはバックスウィング時間Ｔ１＝１．２９５（秒）、ダウンスウィング時間Ｔ２＝０．４０５（秒）のスウィングであり、Ｎ氏のゴルフスウィングはバックスウィング時間Ｔ１＝０．９０７（秒）、ダウンスウィング時間Ｔ２＝０．２９３（秒）のスウィングである。
【００４８】
図１０（ａ）および（ｂ）は、２人のゴルファ（Ｈ氏およびＮ氏）のゴルフスウィングについて、ゴルフクラブシャフトの硬さを表す振動数と、それを把持してゴルフスウィングを行った際のゴルフクラブヘッドのインパクト直前のヘッドスピードとの相関を表した散布図で、スウィングシミュレーションを行って得られたものである。図１０（ａ）はＨ氏についての散布図、図７（ｂ）はＮ氏についての散布図を示す。
図１０（ａ）および（ｂ）に示す散布図では、固有値解析によって求められた固有振動数の値をゴルフクラブシャフトの硬さを表す振動数の値として用いて、ゴルフクラブシャフトの硬さを表す振動数とゴルフクラブヘッドのインパクト直前のヘッドスピードとの相関を表している。上述のように、ゴルフクラブモデルのグリップ側を固定してその反対側を振動させたときの固有振動数は、このときの振動形態がゴルフスウィング時のゴルフクラブの変形形態に略一致しており、ゴルフクラブモデルの固有値解析によって求まる固有振動数の値は、従来からいわれているゴルフクラブシャフトの硬さを表す振動数に対応するものといえる。このことから、固有値解析によって求められた固有振動数の値をゴルフクラブシャフトの硬さを表す振動数の値として用いることで、ゴルフクラブシャフトの硬さを表す振動数とゴルフクラブヘッドのインパクト直前のヘッドスピードとの相関を表すことができるといえる。
【００４９】
図１０（ａ）および（ｂ）から判断されるように、ヘッドスピードはゴルフクラブシャフトの硬さを表す振動数の値に応じて一定周期で変動しており、複数の振動数の値において極大値をもっている。
ここでいう極大値とは、図１０（ａ）および（ｂ）に示す散布図をカーブフィッティング手法を用いて適当な関数で近似した際、この関数が極大値をもつ振動数の値を中心とした所定の振動数範囲におけるヘッドスピードの最大値のことをいう。
【００５０】
図１０（ａ）および（ｂ）から判断できるように、Ｈ氏とN氏はともに一般的に平均的な硬さのゴルフクラブシャフトの振動数として用いられる２４０ｃｐｍ〜２８０ｃｐｍの振動数範囲において、１つの基準極大値（図１０（ａ）中でＨ２、図１０（ｂ）中でＮ２として示す）を有している。
また、Ｈ氏およびＮ氏ともに、ヘッドスピードの値が基準極大値となる際の振動数よりも低い振動数で、上記基準極大値と隣り合う極大値（図１０（ａ）中でＨ１、図１０（ｂ）中でＮ１として示す）を有し、また、ヘッドスピードが基準極大値となる際の振動数よりも高い振動数で、基準となる極大値と隣り合う極大値（図１０（ａ）中でＨ３、図１０（ｂ）中でＮ３として示す）を有している。
ヘッドスピードがこれら極大値をとる際の各振動数は、各ゴルファ毎のゴルフスウィングにおいて比較的高いヘッドスピードとなる振動数であり最適振動数の候補値であるといえる。
【００５１】
図１１は、図１０（ａ）および（ｂ）に示す、Ｈ氏およびＮ氏のゴルフスウィングにおけるゴルフクラブの最適振動数候補値とバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２との関係を示すグラフである。
具体的には、ヘッドスピードが基準極大値（Ｈ２およびＮ２）をとる際の最適振動数候補値Ｆｘ２（Ｆｈ２＝２４５．９２５（ｃｐｍ）およびＦｎ２＝２６８．２５９（ｃｐｍ））と、Ｈ氏およびＮ氏のバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２との関係を、Ｈ氏およびＮ氏それぞれについてプロットした点（図１１中のＰｈ２およびＰｎ２）をそれぞれを結んだグラフ（図１１中式Ｌ２で表されるグラフ）と、最適振動数候補値Ｆｘ２（Ｆｈ２およびＦｎ２）よりも低い振動数で、基準極大値と隣り合う極大値（Ｈ１およびＮ１）をとる際の最適振動数候補値Ｆｘ１（Ｆｈ１＝２０８．８８８（ｃｐｍ）およびＦｎ１＝１９６．５８７（ｃｐｍ））の値と、Ｈ氏およびＮ氏のバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２との関係を、Ｈ氏およびＮ氏それぞれについてプロットした点を結んだグラフ（図１１中式Ｌ１で表されるグラフ）と、最適振動数候補値Ｆｘ２（Ｆｈ２およびＦｎ２）よりも高い振動数で、基準極大値と隣り合う極大値（Ｈ３およびＮ３）をとる際の最適振動数候補値Ｆｘ３（Ｆｈ３＝２８７．４０７（ｃｐｍ）およびＦｎ３＝３３５．８３６（ｃｐｍ））の値と、Ｈ氏およびＮ氏のバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２との関係を、Ｈ氏およびＮ氏それぞれについてプロットした点を結んだグラフ（図１１中式Ｌ３で表されるグラフ）である。
【００５２】
これら各グラフＬ１〜Ｌ３は、それぞれ、
Ｌ１；Ｆｘ１・Ｔ２＝１．１５９８Ｔ１−０．０９１９、および
Ｌ２；Ｆｘ２・Ｔ２＝０．８３４９Ｔ１＋０．５６０９、および
Ｌ３；Ｆｘ３・Ｔ２＝０．７７３２Ｔ１＋０．９３８７として表せる。
【００５３】
これらの式より、所定のゴルファのスウィングにおける複数の最適振動数候補値（Ｆｘ１〜Ｆｘ３）とバックスウィング時間Ｔ１（秒）、ダウンスウィング時間Ｔ２（秒）との関係は、上記式（１）において所定のａおよびｂの組み合わせの値が設定された下記式（２）〜（４）で表される。
Ｆｘ１＝１．１５９８・（Ｔ１／Ｔ２）−０．０９１９／Ｔ２  （２）
Ｆｘ２＝０．８３４９・（Ｔ１／Ｔ２）＋０．５６０９／Ｔ２  （３）
Ｆｘ３＝０．７７３２・（Ｔ１／Ｔ２）＋０．９３８７／Ｔ２  （４）
これら最適振動数候補値と、ゴルフスウィングにおけるバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２との関係を表す式（２）〜（４）を用い、Ｈ氏およびＮ氏以外のゴルファについても、ゴルファ毎に複数の最適振動数候補値が求められる。
【００５４】
表１に、異なるゴルファＡ、Ｂ、Ｃのバックスウィング時間Ｔ１およびダウンスウィング時間Ｔ２を用いて算出したゴルファＡ，Ｂ，Ｃ各々のゴルフスウィングにおける最適振動数候補値を示す。
表１中の最適振動数候補値Ｆｘ１は、上記式（２）を用いて算出したゴルファＡ，Ｂ，Ｃ各々の最適振動数候補値を示し、同様に、表１中の最適振動数候補値Ｆｘ２は上記の式（３）を用いて算出したゴルファＡ，Ｂ，Ｃ各々の最適振動数候補値を、表１中の最適振動数候補値Ｆｘ３は上記式（４）を用いて算出したゴルファＡ，Ｂ，Ｃ各々の最適振動数候補値である。
ここで、ゴルファＡ，Ｂ，Ｃともに上級者のゴルファであり、ゴルファＡのゴルフスウィングはバックスウィング時間Ｔ１＝１．２１（ｓ）、ダウンスウィング時間Ｔ２＝０．２７（ｓ）、ゴルファＢのゴルフスウィングはバックスウィング時間Ｔ１＝１．０４（ｓ）、ダウンスウィング時間Ｔ２＝０．３０（ｓ）、ゴルファＣのゴルフスウィングはバックスウィング時間Ｔ１＝１．１９（ｓ）、ダウンスウィング時間Ｔ２＝０．３８（ｓ）である。
【００５５】
【表１】
2005034614.tif 000003
【００５６】
ゴルファＡ，Ｂ，Ｃともに比較的速いスウィングスピードのスウィングを行うゴルファであり、一般的に実際の使用に適した実用的な振動数範囲とされるゴルフクラブシャフトの振動数の値の範囲は、おおよそ２４０ｃｐｍ〜３００ｃｐｍの範囲である。
この範囲から外れるものは、ゴルフスウィングの際、ゴルファのフィーリングで極端に硬い、柔らかいと感じられたり、ゴルフクラブのしなりが極端に小さく、または極端に大きく感じられ、ゴルファに不快感を与える。
上記表１に示す略最適振動数の値の中には、上述の実際の使用に適した実用的な振動数範囲から外れた範囲の振動数の値もある。このような、適正な振動数範囲から外れる振動数を有するゴルフクラブは、たとえヘッドスピードを略最大化させるものであっても、ゴルファにとって略最適なクラブとして適当ではない。
【００５７】
上記のように、３つの式（２）〜式（４）を用いることで、所定のゴルファのスウィングに対し最適振動数候補値を複数求め、この中から適正な振動数範囲にある最適振動数候補値をゴルファ毎の最適振動数候補値として抽出することが出来（表１中（＊）で示す振動数）、３人のゴルファＡ、Ｂ，Ｃ全てについて適正な振動数範囲の最適振動数を求めることが出来る。表１に示すゴルファＣの場合のように、本発明において選択される最適な振動数は１つの振動数のみであることに限らず、実際の使用に適した実用的な振動数範囲内で複数の最適振動数候補値を最適振動数として選択しても構わない。
本発明では、このようにａ，ｂの値の組み合わせを複数用いることで、１つの式のみを用いる場合に比べて幅広い範囲のゴルファに対して最適な振動数を絞りこむことが可能となる。
【００５８】
本発明では、このように、ゴルフスウィングのヘッドスピードを略最大にするゴルフクラブシャフトの最適振動数を選定することができる。また、複数の候補値（最適振動数の候補値）を絞り込むことができ、この最適振動数に基づいて最適なゴルフクラブを選定することができる。この際、ゴルフ販売店等などで、ゴルフクラブシャフトの振動数毎にグループ分けされた種々のゴルフクラブを用意しておき、求められた最適振動数に基づいてゴルフクラブシャフトが最適な振動数を有する複数のゴルフクラブからなるゴルフクラブ群を選択し、選択されたゴルフクラブ群の複数のゴルフクラブ中から、自分に合った最適なゴルフクラブを選択するようにしてもよい。
例えば、ゴルフクラブシャフトの振動数が全て同一であり（または所定の振動数範囲にあり）、ゴルフクラブシャフトの曲げ剛性分布がそれぞれ異なる複数のゴルフクラブからなるゴルフクラブ群を、所定の振動数（または所定の振動数範囲にあり）毎に複数用意しておき、求められた最適振動数に基づいて最適な振動数の（または所定の振動数範囲の）ゴルフクラブ群を選択する。このゴルフクラブ群のゴルフクラブを用いて実際にゴルフスウィングを試行することで、ゴルフクラブシャフトが全て最適な振動数（または所定の振動数範囲にあり）を有し、かつ曲げ剛性分布がそれぞれ異なる複数のゴルフクラブから、自分に合った最適な曲げ剛性分布をもったゴルフクラブを選択する。
このようにすることで、ゴルフ販売店等に大量に用意された種々のゴルフクラブの中から自分に合った最適なゴルフクラブを迅速に選択することができる。
【００５９】
以上、本発明のゴルフクラブシャフトの最適振動数の決定方法およびゴルフクラブの選定方法について説明したが、本発明のゴルフクラブシャフトの最適振動数の決定方法およびゴルフクラブの選定方法は上記実施例に限定されるものでなく、本発明の要旨を逸脱しない範囲において、各種の改良および変更を行ってもよいのはもちろんである。
【図面の簡単な説明】
【００６０】
【図１】本発明のゴルフクラブシャフトの最適振動数の決定方法およびゴルフクラブの選定方法を実施する装置の概略構成図である。
【図２】図１で示す計測装置の主要部の構成を説明する説明図である。
【図３】本発明のゴルフクラブシャフトの最適振動数の決定方法において用いるバックスウィング時間とダウンスウィング時間とを説明する説明図である。
【図４】（ａ）〜（ｃ）は、図１に示す装置で得られる３次元時系列データの一例を示す図である。
【図５】本発明のゴルフクラブシャフトの最適振動数の決定方法で用いられる算出式を説明する図である。
【図６】従来の最適振動数の決定方法で用いられるスウィング時間と最適振動数との関係を示す図である。
【図７】（ａ）〜（ｃ）は、ゴルフスウィングのシミュレーションの解析に用いるゴルフクラブの有限要素モデルの説明図である。
【図８】図７に示す有限要素モデルを用いてゴルフスウィングのシミュレーションの解析を行ったときの解析結果の一例を示す図である。
【図９】図７に示す有限要素モデルを用いて固有値解析を行ったときの解析結果の一例を示す図である。
【図１０】（ａ）および（ｂ）は、ゴルフクラブモデルが再現するゴルフクラブの、ゴルフクラブシャフトの硬さを表す振動数とゴルフクラブヘッドのインパクト直前のヘッドスピードとの相関を表す散布図である。
【図１１】本発明のゴルフクラブシャフトの最適振動数の決定方法で用いられる算出式を説明する図である。
【符号の説明】
【００６１】
１０  装置
２０  計測装置
２０ａ  トランスミッタ
２０ｂ  レシーバ
２０ｃ  コントローラ
２０ｄ  駆動回路
２０ｅ  検出回路
２０ｆ  制御ユニット
２４  解析装置
２４ａ  データ抽出部
２４ｂ  時間検出部
２４ｃ  振動数算出部
２４ｄ  選定部
２５  モニタ
２６  操作系
３０  ゴルフクラブ
３２  ゴルフクラブシャフト
３４  ゴルフクラブヘッド
３６  グリップ部
４０  ゴルフクラブモデル
４２  ゴルフクラブシャフトモデル
４４  質点モデル
４６  グリップモデル
【図１】
2005034614.tif 000004
【図２】
2005034614.tif 000005
【図３】
2005034614.tif 000006
【図４】
2005034614.tif 000007
【図５】
2005034614.tif 000008
【図６】
2005034614.tif 000009
【図７】
2005034614.tif 000010
【図８】
2005034614.tif 000011
【図９】
2005034614.tif 000012
【図１０】
2005034614.tif 000013
【図１１】
2005034614.tif 000014