(57)【要約】【課題】ゴルファーにとってヘッドスピードが速くなるようなゴルフクラブシャフトを判定可能とし、ゴルファーに最適なゴルフクラブシャフトを選定可能とする。【解決手段】スイング中のシャフトの変形挙動を測定するシャフト挙動測定手段と、シャフトのＥＩ分布を算出するシャフトＥＩ算出手段と、スイング中のシャフトの変形形状を計算するシャフト形状計算手段とを有する第１解析システムと、ゴルファーのスイングを解析分類するスイング分類手段を有する第２解析システムとを備え、上記スイング中の上記シャフトの変形挙動を解析すると共に、上記ゴルファーのスイングを分類し、上記ゴルファーにとってヘッドスピードが最適なシャフトを選定する。【選択図】　図７
2004129687.tif 000002
【特許請求の範囲】
【請求項１】
スイング中のシャフトの変形挙動を測定するシャフト挙動測定手段と、
シャフトのＥＩ分布を算出するシャフトＥＩ算出手段と、
スイング中のシャフトの変形形状を計算するシャフト形状計算手段とを有する第１解析システムと、
ゴルファーのスイングを解析分類するスイング分類手段を有する第２解析システムとを備え、
上記スイング中の上記シャフトの変形挙動を解析すると共に、上記ゴルファーのスイングを分類し、上記ゴルファーにとってヘッドスピードが最適なシャフトを選定することを特徴とする最適ゴルフクラブシャフト選定システム。
【請求項２】
上記シャフト挙動測定手段によりスイング中のシャフトに生じるひずみ量を測定し、上記シャフト形状計算手段によりインパクト時の変形形状を計算すると共に、
上記シャフト分類手段によりゴルファーのスイングの特徴を抽出し、ヘッドスピードが速いシャフトを選定する請求項１に記載の最適ゴルフクラブシャフト選定システム。
【請求項３】
上記シャフトは繊維強化樹脂製とし、該シャフトにひずみゲージを４箇所以上取り付け、該各箇所に生じるひずみを測定し、上記ひずみの測定データから上記シャフトに生じる曲げモーメントを計算すると共に、複合材料に関する積層理論を用いて上記ＥＩ分布を算出し、
上記第１解析システムにより、撓み曲線の微分方程式、上記曲げモーメント、上記ＥＩ分布を用いてシャフト変形形状を計算し、シャフトの変形挙動を解析する請求項１または請求項２に記載の最適ゴルフクラブシャフト選定システム。
【請求項４】
上記スイング分類手段により、ゴルファーのスイングを、腕ターンあるいはボディターン、コック有あるいはコック無から選択される４つのパターンに分類している請求項１乃至請求項３のいずれか１項に記載の最適ゴルフクラブシャフト選定システム。
【請求項５】
上記スイング分類は、ＤＬＴ法を用いたスイング画像処理データ、あるいは、スイング時のシャフトに生じるひずみの時系列データを用いて行っている請求項４に記載の最適ゴルフクラブシャフト選定システム。
【請求項６】
上記ゴルファーがスイングした際の、三次元あるいは二次元座標系における手首関節角度の変化量によりコックの有無を判定し、
上記ゴルファーのスイング時のダウンスイング開始時刻からインパクトまでの三次元あるいは二次元座標系における頭・手首・左腰の変位量から求められる（腰部変位量）／（頭部変位量）、あるいは、（左腰変位量）の２つのパラメータによりボディターンと腕ターンを判定している請求項５に記載の最適ゴルフクラブシャフト選定システム。
【請求項７】
請求項１乃至請求項６のいずれか１項に記載の最適ゴルフクラブシャフト選定システムを用いて選定されたことを特徴とするゴルフクラブシャフト。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、最適ゴルフクラブシャフト選定システム及びゴルフクラブシャフトに関し、詳しくは、ゴルフクラブスイング中のシャフトの変形挙動を解析すると共に、ゴルファーのスイングを分類し、ゴルファーにとって最適なゴルフクラブ、特に、ヘッドスピードが速くなるゴルフクラブを選定するものである。
【０００２】
【従来の技術】
ゴルファーのスイングフォームは千差万別であるため、各ゴルファーにとって最適なゴルフクラブはそれぞれ異なっている。使用するゴルフクラブによりスコアに影響が出ることから、自分に合ったゴルフクラブを探し出すことはゴルファーにとって重要視されている。
【０００３】
ゴルファーにとって最適なゴルフクラブとしては、種々の要因が挙げられるが、特に大きな飛距離が得られることが要求され、よってスイング時のヘッドスピードが速くなるようなゴルフクラブが望まれている。従来、飛距離が得られるクラブを選定するためには、ゴルフ場及び練習場等で、実際にボールを打って飛距離、打ち出し方向等を計測することでクラブの適否を経験的、感覚的に判断している。
【０００４】
その他、スイング画像によりスイングフォームをチェックする、あるいは、シャフトにひずみゲージ等のセンサーを取り付け、シャフトのしなり具合を計測する等によりゴルフクラブの挙動等を計測し、ゴルファーに最適な条件を調査する種々の方法が提案されている。
【０００５】
例えば、特許第３２４３２１０号のスイング中のゴルフシャフトしなり挙動解析システムでは、ゴルフシャフトに取り付けた４ヶ所のひずみゲージにより測定したひずみの時系列と、シャフトに静的荷重を与えることによって求めたシャフトのＥＩ分布を用いて、シャフトの曲げ・曲率分布の変化率（時間微分）を含むしなり挙動を求めている。
【０００６】
また、特許第３２４３２０９号のゴルフクラブシャフトの調子の判定方法、その方法に表現される特性を特徴としたゴルフシャフトでは、特許第３２４３２１０号のしなり挙動解析システムを用いて、シャフトに発生する曲げモーメント分布を求め、多変量解析、クラスター分析によりゴルファーにスイングパターンを分類し、調子を判定している。
【０００７】
具体的には、以下の▲１▼、▲２▼に示す手法により調子を判定している。
▲１▼インパクト手前のスイング方向の曲げモーメント分布とシャフトのＥＩ分布からシャフトの曲率分布の変化率を解析し、シャフトを３等分あるいは等分した場合、全体の曲率分布の変化率に対するチップ・バット・中間の曲率分布の変化率の割合（動的なフレキシブル・レシオ：ＤＦＲ　Ｒａｔｉｏ）によって
（ｉ）ＤＦＲ　Ｒａｔｉｏ　Ｔｉｐ＞２０％　→先調子
ＤＦＲ　Ｒａｔｉｏ　Ｔｉｐ＜２０％　→手元・中調子
（ｉｉ）ＤＦＲ　Ｒａｔｉｏ　Ｂｕｔｔ／ＤＦＲ　ｍｉｄｄｌｅ＞１．２　手元調子
ＤＦＲ　Ｒａｔｉｏ　Ｂｕｔｔ／ＤＦＲ　ｍｉｄｄｌｅ＜１．２　中調子と判別する。また、曲率分布の変化率のピーク位置（動的なキックポイント）によって調子を判別する。
【０００８】
▲２▼インパクト手前のトウダウン方向の最大曲げモーメント分布、あるいは、トップオブスイングやコック解放時付近の最大曲げモーメント分布、シャフトのＥＩ分布から求まる曲率分布によってシャフトの特性を明確にする。
【０００９】
従来までの調子を判定するシステム（順式や逆式のフレックスの比・静的座屈による最大撓み位置・静的荷重を負荷することによる各部の曲率）では「調子」の性能が官能評価とは一致しない場合があったが、これらの手法によって、スイング中のシャフトの調子と官能評価が一致するようなゴルフシャフトの調子判定方法を開示している。
【００１０】
【特許文献１】
特許第３２４３２１０号
【００１１】
【特許文献２】
特許第３２４３２０９号
【００１２】
【発明が解決しようとする課題】
しかしながら、上記特許第３２４３２１０号では、スイング中のしなり挙動を解析することができるが、ヘッドスピードが速くなるようなゴルフクラブシャフトを判定するまでには至らず、ゴルファーに最適なゴルフクラブシャフトを選定できないという問題がある。
【００１３】
また、上記特許第３２４３２０９号では、動的なシャフト変形挙動を解明することで、シャフトの調子とフィーリングとの関係が官能評価と一致するとされているものの、具体的な官能評価と動的シャフト変形挙動（曲率分布の変化率）との関係が示されていない。よって、ヘッドスピードが速くなるようなゴルフクラブシャフトを判定するまでには至らず、ゴルファーに最適なゴルフクラブシャフトを選定できないという問題がある。
【００１４】
さらに、特許第３２４３２０９号、特許第３２４３２１０号で記載されているＥＩ分布の測定方法では、撓み量ｙを計測装置によって求め、シャフト軸方向の距離ｘとｙの最小自乗法によりｙとｘの近似式を算出し、シャフト先端に付加したモーメントと上記近似式からＥＩ分布を計算している。
【００１５】
しかし、この方法では、ＥＩ分布の測定対象が既知の鉄の棒の場合は、既知ＥＩ分布と計測ＥＩ分布は良く一致するが、ゴルフシャフト、特に、図１０に示すように、ＦＲＰ製のシャフト（Ａ〜Ｃの３種類）について、上記記載方法でＥＩ分布を測定すると、最小自乗法の近似の次数によって、チップ・バット部のＥＩ分布が発散してしまい、正確にＥＩ分布を求められないことがある。
【００１６】
本発明は上記した問題に鑑みてなされたものであり、ゴルファーにとってヘッドスピードが速くなるようなゴルフクラブシャフトを判定可能とし、ゴルファーに最適なゴルフクラブシャフトを選定できる最適ゴルフクラブシャフト選定システムを提供することを課題としている。
【００１７】
【課題を解決するための手段】
上記課題を解決するため、本発明は、スイング中のシャフトの変形挙動を測定するシャフト挙動測定手段と、
シャフトのＥＩ分布を算出するシャフトＥＩ算出手段と、
スイング中のシャフトの変形形状を計算するシャフト形状計算手段とを有する第１解析システムと、
ゴルファーのスイングを解析分類するスイング分類手段を有する第２解析システムとを備え、
上記スイング中の上記シャフトの変形挙動を解析すると共に、上記ゴルファーのスイングを分類し、上記ゴルファーにとってヘッドスピードが最適なシャフトを選定することを特徴とする最適ゴルフクラブシャフト選定システムを提供している。
【００１８】
シャフトの変形挙動及びシャフトのＥＩ分布を測定・算出し、シャフトの変形形状を計算すると共に、ゴルファーのスイングの特徴を把握した上で、ゴルファーに最適なシャフトを選定している。このため、各ゴルファーのスイング軌道に応じて、シャフトが適切に変形し、インパクト時にヘッドの重心深度相当分のシャフト先端にスイング方向に十分な撓みが生じるような、ヘッドスピードが速くなるゴルフクラブシャフトを選定することができる。
【００１９】
具体的には、本発明は、鋭意研究の結果、以下のことを見出したことに基づくものである。ヘッドスピードは、シャフトの撓み速度と手元の速度の２つが大きな要因となることを解明し、シャフトの撓み速度が最大になるのは、ヘッドの重心が最下点になるときである、即ち、重心深度相当分だけスイング方向に変形していればヘッドスピードが最も速くなることを解明した。スイングによってシャフト変形方向が異なるので、ヘッド重心がインパクトで最下点になるのもシャフトにより異なることとなるため、スイングパターンによって最適なシャフトが異なることとなる。
【００２０】
また、手元速度は、ゴルファーのスイングパターンに依存し、インパクト前に手元速度が急激に低下する場合や、逆に急激に上昇する場合がある。手元速度が低下する場合、ヘッドスピードのピークで打つためには柔らかいシャフトが最適シャフトとなる。一方、手元速度が上昇する場合、ヘッドスピードのピークで打つためには硬いシャフトが最適シャフトとなる。
【００２１】
また、本発明では、ひずみゲージによる測定ではなく、画像によるスイング分析から最適シャフトを選定することもできる。さらに、上記のようにシャフトの調子や硬さとヘッドスピードの関係を明確にしているため、ヘッドスピードが速くなるようなシャフトを選定可能となった。
【００２２】
即ち、シャフト挙動測定手段によりスイング中のシャフトに生じるひずみ量を測定し、上記シャフト形状計算手段によりインパクト時の変形形状を計算すると共に、上記シャフト分類手段によりゴルファーのスイングの特徴を抽出し、ヘッドスピードが速いシャフトを選定することが好ましい。即ち、撓み曲線の方程式等を用いて、シャフトの撓み量を計算することで、シャフトの変形形状を計算している。
【００２３】
シャフトに生じるひずみ量を測定するとシャフトの挙動を容易かつ正確に把握することができる。また、打球の飛距離や打球方向への影響が大きいインパクト時の変形形状を計算し、より精度の高い選定を可能にしている。
【００２４】
上記シャフトは繊維強化樹脂製とし、該シャフトにひずみゲージを４箇所以上取り付け、該各箇所に生じるひずみを測定し、上記ひずみの測定データから上記シャフトに生じる曲げモーメントを計算すると共に、複合材料に関する積層理論を用いて上記ＥＩ分布を算出し、
上記第１解析システムにより、撓み曲線の微分方程式、上記曲げモーメント、上記ＥＩ分布を用いてシャフト変形形状を計算し、シャフトの変形挙動を解析することが好ましい。
【００２５】
カーボン繊維等の強化繊維と熱硬化性樹脂等の樹脂中に含浸させたプリプレグの積層体等からなる繊維強化樹脂製のシャフトにおいても、複合材料に関する積層理論を用いてＥＩ分布を算出しているため、より正確なＥＩ分布を得ることができ、解析精度を向上することができる。ひずみゲージは複数個が良いが精度を高めるには４個以上が良く、シャフトの軸方向にできるだけ均等に、少なくともチップ側とバット側の両端に配置されるのが好ましく、有線あるいは無線により検知した電圧信号を送信することができる。ひずみゲージは、非常に軽量で小さいため、ゴルファーにとってスイング時に違和感もなく、測定から検知された信号を処理する時間も短い。なお、金属製等のシャフトの解析も可能であり、ドライバー、アイアン、パター用等、全ての種類のゴルフクラブシャフトの解析が可能である。
【００２６】
複合材料に関する積層理論とは、積層からなる材料（複合材料）のプライ（巻き付け）によって弾性定数は異なるので、応力は一様にならないことを考慮し、プライ間の平均応力を定義することによって多方向の積層材の等価弾性定数に置き換えた理論である。これによれば、シャフトを等分割（あるいは複数の要素）に分けてＥＩを求めることができ、チップ、バット部についても、シャフト軸方向のＥＩを補間せずに求めているので、他の部分（シャフト中央部）と同様の扱いで計算できる。よって、シャフトのチップ（ヘッド）側、バット（グリップ）側端部のＥＩが発散しないＥＩ分布を求めることができる。
【００２７】
さらに詳述すると積層理論とは、以下の２つの仮定の基に成り立つ理論である。
▲１▼積層板は対称である。従って、プライの配向角や弾性係数は積層板の中央面に対して対称である。
▲２▼ひずみは積層板の厚さを通じて一定である。これは、積層板の厚みが長さ・幅に対して小さいためである。
積層理論では、以上の仮定を基に、「積層板はプライによって弾性定数は異なるので応力分布は一様でない。従って、平均応力を定義すると容易になる。」という考えに基づき弾性定数を求めている。
【００２８】
積層方向の平均応力は下記の式１、積層１枚についての応力は斜向座標系での応力−ひずみの関係より下記の式２で求められる。また、図１（Ａ）に積層材料ＳＳの積層方向を示す。プリプレグをシャフトにした時のシャフト軸方向を方向１、方向１に垂直方向を方向２、せん断方向を方向６としている。厚み方向Ｚ、厚みはｈとしている。Ｑとは斜向軸系の工学的定数である。
【００２９】
【式１】
2004129687.tif 000003
【００３０】
【式２】
2004129687.tif 000004
【００３１】
式２を式１へ代入すると下記の式３が得られ、ひずみは仮定より定ひずみであるため以下のように変形され面内の合応力が得られる。
【００３２】
【式３】
2004129687.tif 000005
【００３３】
斜向軸系の工学的定数Ｑと弾性係数の関係は単軸引張試験を想定すると明確となる。図１（Ｂ）に示す積層材料において、等価弾性係数Ｅを求めるために下記の関係を得る。
【００３４】
【式４】
2004129687.tif 000006
【００３５】
式４を式３に代入し、また一様平均な応力が積層材に生じている（つまりプライ配向角＝０ｄｅｇの等方性材料）ので下記の関係が得られる。Ａ１１、Ａ１２はＱ１１、Ｑ１２の加算である。ＥＩを考える場合は、各プライの工学的定数の断面二次モーメントを掛け、加算することで積層材の等価ＥＩを計算できる。
【００３６】
【式５】
2004129687.tif 000007
【００３７】
上記スイング分類手段により、ゴルファーのスイングを、腕ターンあるいはボディターン、コック有あるいはコック無から選択される４つのパターンに分類していることが好ましい。ゴルファーのスイングをこのようなパターンに分類することで、各ゴルファーのスイングの特徴別に適したゴルフクラブシャフトの選定を行うことができる。
【００３８】
上記スイング分類は、ＤＬＴ法等を用いたスイング画像処理データ、あるいは、スイング時のシャフトに生じるひずみの時系列データを用いて行うことができる。
画像処理を用いてスイングパターンを分類するためには、人体各部（頭・首・肩・肘・手首・腰）等の複数箇所とシャフトの挙動を計測する必要がある。また、ひずみの時系列データを用いてスイングパターンを分類するためには、ひずみの波形データとスイング分類との相関関係を把握する必要がある。
【００３９】
上記ゴルファーがスイングした際の、三次元あるいは二次元座標系における手首関節角度の変化量によりコックの有無を判定し、上記ゴルファーのスイング時のダウンスイング開始時刻からインパクトまでの三次元あるいは二次元座標系における頭・手首・左腰の変位量から求められる（腰部変位量）／（頭部変位量）、あるいは、（左腰変位量）の２つのパラメータによりボディターンと腕ターンを判定することができる。
【００４０】
上記２つのパラメータは、予めレッスンプロ等のスイング診断の結果により、コックの有無の判定及びボディターンと腕ターンの判定との相関が得られているものである。このように事前に診断した相関関係と照合することで、判定精度をより高めている。
【００４１】
また、本発明は、本発明の最適ゴルフクラブシャフト選定システムを用いて選定されたことを特徴とするゴルフクラブシャフトを提供している。
本発明の最適ゴルフクラブシャフト選定システムでは、曲げ剛性分布等の特性が異なるシャフトの中からゴルファーに最適なゴルフクラブシャフトを選定することができる。このため、選定されたゴルフクラブシャフトは、各ゴルファーにとって使いやすく、容易に大きな飛距離が得られるため、非常に有用なものである。
【００４２】
【発明の実施の形態】
以下、本発明の最適ゴルフクラブシャフト選定システムの第１実施形態を図面を参照して説明する。
最適ゴルフクラブシャフト選定システムは、スイング中のシャフトの変形挙動を測定するシャフト挙動測定手段と、シャフトのＥＩ分布を算出するシャフトＥＩ算出手段と、スイング中のシャフトの変形形状を計算するシャフト形状計算手段とを有する第１解析システムと、ゴルファーのスイングを解析分類するスイング分類手段を有する第２解析システムとを備えており、スイング中のシャフトの変形挙動を解析すると共に、ゴルファーのスイングを分類し、ゴルファーにとってヘッドスピードが最適なシャフトを選定する構成としている。
【００４３】
具体的には、シャフト挙動測定手段によりスイング中のシャフトに生じるひずみ量を測定し、シャフト形状計算手段によりインパクト時の変形形状を計算すると共に、シャフト分類手段によりゴルファーのスイングの特徴を抽出し、ヘッドスピードが速いシャフトを選定するものである。
【００４４】
以下、最適シャフトの選定までの手順について説明する。
挙動解析を行うシャフトを準備する。図２に示すように、シャフト１は、カーボン繊維を強化繊維としたプリプレグの積層体からなり、小径側先端にヘッド２が取り付けられ、大径端側にグリップ３が取り付けられている。シャフト１には、ひずみゲージＧ（Ｇ１〜Ｇ４）が４箇所貼り付けられている。具体的には、ヘッド取付端から８５ｍｍの位置に第１ひずみゲージＧ１、さらに２００ｍｍずつ間隔をあけて第２ひずみゲージＧ２、第３ひずみゲージＧ３、第４ひずみゲージＧ４がシャフト１の長軸方向の同一軸線上の同面でスイング方向に対して垂直な面に貼り付けられ、スイング時にシャフト１に生じるひずみを測定可能な構成としている。
【００４５】
ひずみゲージＧで検出された電圧信号は、ブリッジボックスとアンプを介して動ひずみ計に送信され（図示せず）、コンピュータに取り込まれ、シャフト１の各ひずみ測定箇所の電圧信号を入力できる。スイング測定終了後、キャリブレーションを行うことによって、電圧をひずみに変換し、図３に示すように、各測定箇所でのひずみの時系列をシャフト挙動測定手段で計測している。図中、ひずみ（ｓｔｒａｉｎ）のグラフが０に近い側からＧ１、Ｇ２、Ｇ３、Ｇ４のデータとしている。
【００４６】
次に、複合材料に関する積層理論、即ち、複合材料の材料定数算出方法で用いられる積層理論を用いてシャフトのＥＩ分布を算出する。
具体的には、実測にて繊維軸系のヤング率・ポアソン比を求める。また、以下に繊維軸系の応力−ひずみの関係式を示す。これは、単軸縦方向、単軸横方向、縦方向せん断方向についてのひずみと応力の関係について、線形弾性材料に限定することで、重ね合わせを行うことで求めている。
【００４７】
【式６】
2004129687.tif 000008
【００４８】
式中、σは応力、Ｅは縦弾性係数、ｍは１／（１−νｘνｙ）で表される定数、εはひずみ、νはポアソン比を示す。ｘは引張方向、ｙは圧縮方向、Ｓはせん断方向を指す。
【００４９】
即ち、斜向軸系のひずみをひずみ変換則、斜向軸系の応力を応力逆変換則でそれぞれ斜向軸系と繊維軸系の関係を算出し、上記の応力−ひずみ関係式に代入することによって、繊維軸系の応力−ひすみ関係を求める。
また、繊維軸系の応力−ひずみ関係から斜向軸系の材料定数を決定し、積層されるプリプレグ１枚分の剛性を求める。
上記方法で繊維材であるプリプレグ毎の斜向軸系の材料定数を決定し、積層理論によってシャフトＥＩ算出手段によりＥＩ分布を算出する。
即ち、繊維軸系とは実測で求めた軸系に相当し、例えば、引張試験をした場合の軸系になる。斜向軸系とは、１枚のプリプレグを巻きつけた場合の軸系（シャフト長手方向を軸とする）になる。
【００５０】
そして、上記ひずみの測定データから測定各箇所においてシャフトに生じる曲げモーメントを計算する。撓み方向の変位をｙ、シャフト軸方向の距離をｘ、曲げモーメントをＭ、曲げ剛性をＥＩ、ひずみをε、シャフト内径をｒとすると、下記の撓み曲線の微分方程式が得られる。プリプレグを巻きつけると周方向に不均一な厚みになるのでｒを均一化して求めている。
【００５１】
【式７】
2004129687.tif 000009
【００５２】
曲げモーメントＭをｘを変数とする２次式で近似し、最小自乗近似によって下記式の係数ａ、ｂ、ｃを算出する。
【００５３】
【式８】
2004129687.tif 000010
【００５４】
算出された曲げモーメントＭの近似式を撓み曲線の微分方程式に代入し、１回積分することによって、撓み角度θ、２回積分することによって、撓み量ｙを算出する（ただし、グリップ部の撓み・撓み角度は０とする拘束条件を付加）。
【００５５】
【式９】
2004129687.tif 000011
【００５６】
撓み曲線の微分方程式、上記曲げモーメント、上記ＥＩ分布を用いてシャフト形状計算手段によりシャフト変形形状を計算し、第１解析システムにより、シャフトの変形挙動を解析する。このようにして、シャフトの変形形状を時々刻々算出している。
具体的には、図４において、スイング中のシャフトの変形形状の時間変化を示している。図中の右側から順にインパクトから１ｍｓ、５ｍｓ、１０ｍｓ、１５ｍｓ、２０ｍｓ、２５ｍｓ、３０ｍｓ、３５ｍｓ前としている。このように、シャフト各部の撓み量を計算することでシャフトの変形形状を把握している。
【００５７】
スイング分類手段を有する第２解析システムによりスイングパターンを分類する。具体的には、ゴルファーのスイングを、腕ターンあるいはボディターン、コック有あるいはコック無の４つのパターンに分類している。
【００５８】
スイングパターンを分類するために、ゴルファーの人体各部（頭・首・肩・肘・手首・腰等）とシャフトの動きを計測する。このような画像を計測する方法としてＤＬＴ法を用いる。具体的には、図５に示すように、２台の同期したカメラＣ１、Ｃ２によってスイングを測定し、各カメラ毎の人体各部の２次元座標を求め、既知の空間上の点（コントロールポイント）でカメラのキャリブレーションを行うことによって三次元座標を算出する手法である。図６（Ａ）（Ｂ）（Ｃ）に手首のポインティング結果を各軸（Ｘ軸、Ｙ軸、Ｚ軸）毎に示す。各プロットを５次式で近似しデータを補正した。
【００５９】
これらゴルファーがスイングした際の座標値から、三次元あるいは二次元座標系における手首関節角度の変化量によりコックの有無を判定する。
具体的には、まず、人体各部を正面（ＸＺ平面）に投影し、肘−手首ベクトル、シャフト−手首ベクトルから、その間の手首関節角度を求める。そして、コック解放からインパクト時までの手首関節角度の角度変化が６０ｄｅｇ以上であればコック有り、６０ｄｅｇ未満ではコック無しと判定する。この判定に６０ｄｅｇという数値を用いたのは事前にレッスンプロのスイング診断の結果によるコックの有無の判定と相関が高かったためである。
【００６０】
次に、スイングがボディターンか腕ターンかを判定する方法について説明する。
ゴルファーのスイング時のダウンスイング開始時刻からインパクトまでの三次元あるいは二次元座標系における頭・手首・左腰の変位量から求められる（腰部変位量）／（頭部変位量）あるいは、（左腰変位量）の２つのパラメータにより判定している。
【００６１】
例えば（腰部変位量）／（頭部変位量）の値が１２以上であれば腕ターンと判定し、（左腰変位量）の値が２５０以下であれば腕ターンと判定する。（腰部変位量）／（頭部変位量）と（左腰変位量）のどちらもが腕ターンと判断されない場合はボディターンと判定する。この判定に（腰部変位量）／（頭部変位量）の値が１２、（左腰変位量）の値が２５０という数値を用いたのは、事前にレッスンプロのスイング診断の結果によるボディターン・腕ターンの判定と相関が高かったためである。
【００６２】
これらの解析結果に基づき、ヘッドスピードが速くなるシャフトを選定する。第１解析システム、第２解析システムによって４つのスイング分類に分け、これら４つの分類と、それに適した曲げ剛性分布を有するシャフトとの関係を以下に示す。
【００６３】
シャフトのチップからの距離が、０〜２００ｍｍまでの部分をＡ（先）、２００〜４００ｍｍ（先よりの中央）までの部分をＢ、４００〜６００ｍｍ（中央）までの部分をＣ、６００〜８００ｍｍ（手元よりの中央）までの部分をＤ、８００〜１０００ｍｍ（手元）までの部分をＥとして、各部分のＥＩをそれぞれ以下のように設定している。
【００６４】
１．ボディターン・コック有りの場合
Ｅ部分の剛性を高くし、Ｄ部分の剛性を低くする。（後述する実施例のシャフト▲２▼のようなシャフト）
２．ボディターン・コック無しの場合
Ｅ部分の剛性を高くし、Ｂ・Ｃ・Ｄ部分の剛性をチップ先端から徐々に高くする。（後述する実施例のシャフト▲５▼のようなシャフト）
３．腕ターン・コック有りの場合
シャフトの剛性（フレックス）を低くする。例えば、シャフトの硬さをＳ→Ｒにする。（後述する実施例のシャフト▲４▼のようなシャフト）
４．腕ターン・コック無しの場合
Ｅ部分の剛性を低くする。（後述する実施例のシャフト▲３▼のようなシャフト）あるいは、Ｅ部分の剛性を高くし、Ｂ・Ｃ・Ｄ部分の剛性をチップ先端から徐々に高くする。（後述する実施例のシャフト▲５▼のようなシャフト）
【００６５】
このように各分類毎に選定されるヘッドスピードが最適なシャフトは、そのシャフトの変形形状がインパクト時にヘッドスピードが最適になる形状になっている。具体的には、ヘッドの重心深度相当分のシャフト先端に撓みが生じているので、ゴルファーにとってヘッドスピードが速くなる最適なゴルフシャフトを提供することができる。
【００６６】
このように、ゴルファーのスイングを画像処理、シャフトの変形をひずみゲージによって測定することで、ゴルファーのスイングを分類し、分類毎にヘッドスピードが最適なシャフトを選定することができる。
【００６７】
また、上記のように最適ゴルフクラブシャフト選定システムを用いて選定されたゴルフクラブシャフトは、シャフトの変形挙動に加え、ゴルファーのスイングを考慮しているため、ゴルファーのヘッドスピードが速くなり、非常に大きな飛距離を得ることができる。
【００６８】
また、スイング分類は、上記実施形態のようにＤＬＴ法を用いたスイング画像処理データを用いても良いが、以下に示すように、スイング時のシャフトに生じるひずみの時系列データを用いて行っても良い。
【００６９】
上記４つのスイング分類と、ひずみの波形データに基づくシャフトの変形挙動の関係を以下に示す。
図７、８に示すように、画像によるスイング分類とシャフトに生じるひずみの時系列（上記実施形態のように４箇所にひずみゲージを取り付けた場合）の関係から、ボディターンと腕ターン、コックの有無はそれぞれ１対１の関係になっていることがわかる。なお、ゲージ番号とＣＨの番号は一致している。図中、ひずみ（ｓｔｒａｉｎ）のグラフが０に近い側から１ｃｈ、２ｃｈ、３ｃｈ、４ｃｈのデータとしている。
【００７０】
１．ボディターン・コック有りの場合
ひずみのピークが２つで、最初のピークが大きい。
２．ボディターン・コック無しの場合
ひずみのピークが３つ。
３．腕ターン・コック有りの場合
ひずみのピークが２つで、最初のピークが小さい。
４．腕ターン・コック無しの場合
ひずみのピークが１つ、あるいは、インパクト前のピークがかなり小さい。
従って、この関係を利用すれば、画像処理ではなく、ひずみの時系列結果からもスイング分類を行うことができる。
【００７１】
以下、本発明の最適ゴルフクラブシャフト選定システムの実施例について詳述する。
図９に示すように５パターンのＥＩ（曲げ剛性）分布を示すシャフト（シャフト▲１▼〜シャフト▲５▼）を準備した。
【００７２】
まず、テスター９名（Ａ〜Ｉ）について、シャフト▲１▼の剛性分布を有する基準シャフトにヘッドを取り付けスイングし、各テスター毎にヘッドスピードを測定した。
【００７３】
次に、上記第１実施形態と同様に、画像によるスイングの測定、及びシャフトのひずみを測定し、スイング分類を行うことによってヘッドスピードが速くなるシャフトを選定した。また、選定結果を確認するために、シャフト▲２▼〜シャフト▲５▼についても各テスターがスイングし、ヘッドスピードを測定した。
各テスター毎のスイング分類の判定結果と各シャフトでのヘッドスピードの測定結果を下記の表１に示す。
【００７４】
【表１】
2004129687.tif 000012
【００７５】
表１に示すように、テスターＡ〜Ｃのスイングは、ボディターン・コック有であるため、本システムの解析結果によれば、シャフト▲２▼が最適となるはずであり、実際のヘッドスピードの測定結果は、この解析結果と一致していた。また、テスターＤ、Ｅのスイングは、腕ターン・コック有であるため、本システムの解析結果によれば、シャフト▲４▼が最適となるはずであり、実際のヘッドスピードの測定結果は、この解析結果と一致していた。
【００７６】
テスターＦ、Ｇのスイングは、ボディターン・コック無であるため、本システムの解析結果によれば、シャフト▲５▼が最適となるはずであり、実際のヘッドスピードの測定結果は、この解析結果と一致していた。また、テスターＨ、Ｉのスイングは、腕ターン・コック無であるため、本システムの解析結果によれば、シャフト▲３▼あるいはシャフト▲５▼が最適となるはずであり、実際のヘッドスピードの測定結果は、この解析結果と一致していた。
【００７７】
このように、各テスター全てにおいて、スイング分類によって分類毎に最適と判断されたシャフトが、５本のシャフトの中で最もヘッドスピードが速くなっているシャフトであることが確認できた。
【００７８】
【発明の効果】
以上の説明より明らかなように、本発明によれば、シャフトの変形挙動及びシャフトのＥＩ分布を算出し、シャフトの変形形状を計算すると共に、ゴルファーのスイングの特徴を把握した上で、ゴルファーに最適なシャフトを選定している。このため、各ゴルファーのスイング軌道に応じて、シャフトが適切に変形し、ヘッドの重心深度相当分のシャフト先端にスイング方向へ十分な撓みが生じるような、ヘッドスピードが速くなるゴルフクラブシャフトを選定することができる。また、複合材料に関する積層理論を用いてＥＩ分布を算出しているため、より正確なＥＩ分布を得ることができ、解析精度を向上することができる。
【００７９】
よって、各ゴルファーは、練習場や販売店等での試打の繰り返しを行うことなく、自分のスイングに適したゴルフクラブシャフトを的確に選定できるため、ゴルフ選びのための無駄なコストや時間を低減することができる。
【００８０】
さらに、本発明のゴルフクラブシャフトは、ゴルファーのヘッドスピードが速くなるようなシャフトであるため、容易に大きな飛距離を得ることができる。
【図面の簡単な説明】
【図１】（Ａ）（Ｂ）は積層理論の説明図である。
【図２】ゴルフクラブシャフトへのひずみゲージの取り付け状況を示す図である。
【図３】スイング中にシャフトに生じるひずみの時刻歴データを示す図である。
【図４】シャフトの変形形状の計算結果を示す図である。
【図５】ＤＬＴ法による測定状況を示す図である。
【図６】（Ａ）（Ｂ）（Ｃ）は、ＤＬＴ法による手首のポインティング結果を示す図である。
【図７】スイング分類（ボディターン）とひずみの波形データに基づくシャフト挙動の関係を示す図である。
【図８】スイング分類（腕ターン）とひずみの波形データに基づくシャフト挙動の関係を示す図である。
【図９】シャフトのＥＩ分布を示す図である。
【図１０】従来のＥＩ分布の計算結果を示す図である。
【符号の説明】
１　シャフト
２　ヘッド
３　グリップ
【図１】
2004129687.tif 000013
【図２】
2004129687.tif 000014
【図３】
2004129687.tif 000015
【図４】
2004129687.tif 000016
【図５】
2004129687.tif 000017
【図６】
2004129687.tif 000018
【図７】
2004129687.tif 000019
【図８】
2004129687.tif 000020
【図９】
2004129687.tif 000021
【図１０】
2004129687.tif 000022