(57)【要約】【課題】たわみ分布を効率良く、高い精度で求めることができるゴルフクラブヘッドのフェース面のたわみ分布測定方法を提供する。【解決手段】本発明は、ゴルフクラブヘッドのフェース面に対して、複数の測定点を予め設定しておき、各測定点に既知の加振力を付与し、このときのゴルフクラブヘッドにおいて各測定点から外れる場所における応答信号を測定する工程と、既知の加振力と応答信号から周波数応答関数を求め、この周波数応答関数から虚部を抽出し、周波数応答関数の虚部と周波数応答関数の周波数との応答関数を作成する工程と、応答関数に基づいて１次共振周波数を求める工程と、この応答関数を１次共振周波数を中心周波数として所定の幅の周波数範囲に区画する工程と、応答関数において区画された周波数範囲における正の値を合計して、この合計値を測定点におけるたわみ量とする工程とを有する。【選択図】図２
2011050440.tif 000002
【特許請求の範囲】
【請求項１】
ゴルフクラブヘッドのフェース面に対して、複数の測定点を予め設定しておき、前記各測定点に、既知の加振力を付与し、このときの前記ゴルフクラブヘッドにおいて前記各測定点から外れる場所における応答信号を測定する工程と、
前記既知の加振力と前記応答信号から周波数応答関数を求め、前記周波数応答関数から虚部を抽出し、前記周波数応答関数の虚部と前記周波数応答関数の周波数との応答関数を作成する工程と、
前記応答関数に基づいて１次共振周波数を求める工程と、
前記応答関数を、前記１次共振周波数を中心周波数として所定の幅の周波数範囲に区画する工程と、
前記応答関数において区画された周波数範囲における正の値を合計して、この合計値を前記測定点におけるたわみ量とする工程とを有することを特徴とするゴルフクラブヘッドのフェース面のたわみ分布測定方法。
【請求項２】
さらに、前記各測定点のたわみ量を用いて、前記フェース面上におけるたわみ分布の中心を求める工程を有する請求項１に記載のゴルフクラブヘッドのフェース面のたわみ分布測定方法。
【請求項３】
前記たわみ分布の中心は、フェース面に設けられた基準点から各測定点までの距離の重み付けによる加重平均により求められるものである請求項２に記載のゴルフクラブヘッドのフェース面のたわみ分布測定方法。
【請求項４】
前記所定の幅の周波数範囲は、前記１次共振周波数を中心周波数とした１／３オクターブバンドの範囲である請求項１〜３のいずれか１項に記載のゴルフクラブヘッドのフェース面のたわみ分布測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、飛距離、スイートスポットの大きさ等のゴルフクラブヘッドの性能の評価に用いられるフェース面のたわみ分布測定方法に関し、特に、たわみ分布を効率良く、高い精度で求めることができると共に、たわみ分布の中心も高い精度で求めることができるゴルフクラブヘッドのフェース面のたわみ分布測定方法に関するものである。
【背景技術】
【０００２】
現在、ゴルフクラブヘッドについて飛距離、スイートスポットの大きさ等の性能の評価をするために、ゴルフクラブヘッドのフェース面の反発係数の分布が測定されている。この反発係数以外に、フェース面のたわみ分布も測定されている（特許文献１等）。
【０００３】
特許文献１には、例えば、フェース面の縁部に加速度検出器を取り付けた状態で、インパルスハンマにてフェース面上の多数の測定点を打撃し、インパルスハンマの加振力（Ｆ）に対する振動加速度（Ａ）の比（Ａ／Ｆ）が、最も大きい点が最大撓み位置に相当することが記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００３−７７２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上述の特許文献１のような最大撓み位置の測定方法では、インパルスハンマの加振力（Ｆ）に対する振動加速度（Ａ）の比（Ａ／Ｆ）を求めて撓みを算出している。しかしながら、特許文献１の測定方法では、加振力（Ｆ）に対する振動加速度（Ａ）の比（Ａ／Ｆ）を求めているだけなので、得られる撓みでは、ゴルフクラブヘッドの飛距離、スイートスポットの大きさ等の性能をより高い精度で評価するには十分ではない。
【０００６】
また、特許文献１の測定方法でたわみ分布の中心を求める場合、インパルスハンマで打撃した点にたわみ分布の中心があるとは限らず、測定精度が十分ではない。特許文献１において、たわみ分布の中心の測定精度を高くするには、インパルスハンマでの打撃点を多くする必要があり、測定に多くの時間を要するという問題点がある。
このように、ゴルフクラブヘッドの飛距離、スイートスポットの大きさ等の性能をより高い精度で評価するためには、フェースのたわみ分布を高い精度で、しかも効率良く求める方法が望まれている。さらには、たわみ分布の中心を高い精度で求めることも望まれている。
【０００７】
本発明の目的は、前記従来技術に基づく問題点を解消し、たわみ分布を効率良く、高い精度で求めることができるゴルフクラブヘッドのフェース面のたわみ分布測定方法を提供することにある。
本発明の他の目的は、たわみ分布の中心を高い精度で求めることができるゴルフクラブヘッドのフェース面のたわみ分布測定方法を提供することにある。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、本発明は、ゴルフクラブヘッドのフェース面に対して、複数の測定点を予め設定しておき、前記各測定点に、既知の加振力を付与し、このときの前記ゴルフクラブヘッドにおいて前記各測定点から外れる場所における応答信号を測定する工程と、前記既知の加振力と前記応答信号から周波数応答関数を求め、前記周波数応答関数から虚部を抽出し、前記周波数応答関数の虚部と前記周波数応答関数の周波数との応答関数を作成する工程と、前記応答関数に基づいて１次共振周波数を求める工程と、前記応答関数を、前記１次共振周波数を中心周波数として所定の幅の周波数範囲に区画する工程と、前記応答関数において区画された周波数範囲における正の値を合計して、この合計値を前記測定点におけるたわみ量とする工程とを有することを特徴とするゴルフクラブヘッドのフェース面のたわみ分布測定方法を提供するものである。
【０００９】
本発明において、さらに、前記各測定点のたわみ量を用いて、前記フェース面上におけるたわみ分布の中心を求める工程を有することが好ましい。
また、本発明において、前記たわみ分布の中心は、フェース面に設けられた基準点から各測定点までの距離の重み付けによる加重平均により求められるものであることが好ましい。
また、本発明において、前記所定の幅の周波数範囲は、例えば、前記１次共振周波数を中心周波数とした１／３オクターブバンドの範囲である。
【００１０】
また、本発明において、例えば、前記応答信号は、加速度センサを用いて測定された加速度信号、歪ゲージを用いて測定された歪量を表す信号（歪信号）、または音圧計を用いて測定された音圧レベルを表す信号（音圧信号）である。
また、本発明において、例えば、前記既知の加振力は、インパルスハンマにより与えられる。
【発明の効果】
【００１１】
本発明によれば、周波数応答関数の虚部と周波数応答関数の周波数との応答関数から１次共振周波数を求め、この応答関数において１次共振周波数を中心周波数として区画された周波数範囲における正の値を合計して、この合計値を測定点におけるたわみ量とすることにより、ゴルフクラブヘッドのフェース面のたわみ分布を高い精度で求めることができる。これにより、飛距離、スイートスポットの大きさ等のゴルフクラブヘッドの性能をより高い精度で評価することができる。
更には、本発明においては、各測定点のたわみ量を用いて、フェース面上におけるたわみ分布の中心を求めることができる。この場合、たわみ量を高い精度で求めているため、たわみ分布の中心も高い精度で求めることができる。
【図面の簡単な説明】
【００１２】
【図１】（ａ）は、本発明の実施形態のゴルフクラブヘッドのフェース面のたわみ分布測定方法に用いられる測定装置を示す模式図であり、（ｂ）は、ゴルフクラブヘッドのフェース面の測定点を示す模式図である。
【図２】（ａ）は、縦軸に周波数応答関数の虚部をとり、横軸に周波数をとって、周波数応答関数の結果を示すグラフであり、（ｂ）は、縦軸に周波数応答関数の虚部をとり、横軸に周波数をとって、たわみ量の求め方を説明するためのグラフである。
【図３】（ａ）〜（ｅ）は、それぞれ、本発明の実施形態のゴルフクラブヘッドのフェース面のたわみ分布測定方法による測定結果の例を示す模式図である。
【発明を実施するための形態】
【００１３】
以下、添付の図面に示す実施形態に基づいて、本発明のゴルフクラブヘッドのフェース面のたわみ分布測定方法を詳細に説明する。
図１（ａ）は、本発明の実施形態のゴルフクラブヘッドのフェース面のたわみ分布測定方法に用いられる測定装置を示す模式図であり、（ｂ）は、ゴルフクラブヘッドのフェース面の測定点を示す模式図である。
【００１４】
図１（ａ）に示す測定装置１０は、例えば、ゴルフクラブヘッド２０のフェース面２２のたわみ量を測定するとともに、たわみ分布を求め、更には、たわみ分布の中心を求めるものである。
図１（ａ）に示す測定装置１０は、所定の加振力をゴルフクラブヘッド２０に加えるインパルスハンマ１２と、ＦＦＴアナライザ１４と、加速度センサ１６と、アンプ１８と、コンピュータ（以下、ＰＣという）１９とを有する。
加速度センサ１６はアンプ１８を介してＦＦＴアナライザ１４に接続されており、インパルスハンマ１２はＦＦＴアナライザ１４に接続されている。
【００１５】
インパルスハンマ１２は、構造物の固有振動数測定やモーダル解析を行うための力センサを内蔵した加振ハンマである。このインパルスハンマ１２で、例えば、ゴルフクラブヘッド２０のフェース面２２上に予め設定された測定点ｐを打撃した場合、内蔵されている力センサから加振力信号がＦＦＴアナライザ１４に出力される。このインパルスハンマ１２を用いることにより、ゴルフクラブヘッド２０のフェース面２２の各測定点ｐに既知の加振力を付与することができる。
【００１６】
加速度センサ１６は、ゴルフクラブヘッド２０のフェース面２２の測定点ｐがインパルスハンマで打撃されたときに、ゴルフクラブヘッド２０に生じる加速度を測定するものである。すなわち、加速度センサ１６は、ゴルフクラブヘッド２０のフェース面２２の各測定点ｐに加振力が印加された場合、ゴルフクラブヘッド２０に生じる加速度を測定するものである。
この加速度センサ１６は、たわみを測定するフェース面２２における測定点ｐから離れた位置に取り付けられるものである。加速度センサ１６は、例えば、ゴルフクラブヘッド２０のホーゼル部２４の下方に取り付けられる。この加速度センサ１６が取り付けられる位置を取付位置Ｄという。
加速度センサ１６には、例えば、軽量加速度ピックアップが用いられる。この軽量加速度ピックアップは、加速度信号を電荷の形態で出力する。なお、加速度ピックアップは超軽量のものが好ましい。例えば、質量が０．６ｇのものを用いるとよい。これは、質量が大きいとフェース面２２の振動に対する影響が大きくなるためである。
【００１７】
アンプ１８は、加速度センサ１６に接続されている。このアンプ１８には、加速度センサ１６で得られた加速度信号が、例えば、電荷の形態で入力される。アンプ１８は、その電荷の形態で入力された各加速度信号を電圧に変換し、増幅してＦＦＴアナライザ１４に加速度信号として出力する。
【００１８】
ＦＦＴアナライザ１４は、アンプ１８から電圧の形態で出力される加速度センサ１６による加速度信号およびインパクトハンマ１２から出力される加振力信号を時系列で取り込み、更には後述のように周波数応答関数を求めるものである。
【００１９】
ＦＦＴアナライザ１４においては、例えば、加速度センサ１６による加速度信号、およびインパクトハンマ１２の加振力信号について、それぞれ打撃してから、例えば、約１ミリ秒のサンプリング周期（４秒／４０９６間隔）で４秒間（サンプリング時間）取り込む。ＦＦＴアナライザ１４により、加速度センサ１６による加速度の時系列データ（時間波形）が取得されるとともに、インパクトハンマ１２による測定点ｐの打撃時の加振力の時系列データ（時間波形）が取得される。
【００２０】
また、ＦＦＴアナライザ１４は、加速度センサ１６による加速度の時系列データ（時間波形）、インパクトハンマ１２による加振力の時系列データ（時間波形）について、高速フーリエ変換（以下、ＦＦＴという）を行う機能を有する。
さらに、ＦＦＴアナライザ１４は、加速度センサ１６による加速度の時系列データ（時間波形）、およびインパクトハンマ１２による加振力の時系列データ（時間波形）に基づいて、周波数応答関数を求める機能を有する。その他、ＦＦＴアナライザ１４は、一般的なＦＦＴアナライザが備える機能を有する。
【００２１】
ＦＦＴアナライザ１４に取得された加速度センサ１６による加速度の時系列データ（時間波形）およびインパクトハンマ１２による加振力の時系列データ（時間波形）はＰＣ１９に出力される。
また、ＦＦＴアナライザ１４でなされたＦＦＴの結果および求められた周波数応答関数についてもＰＣ１９に出力される。
【００２２】
なお、ＦＦＴアナライザ１４による加速度信号、および加振力信号についてのサンプリング周期は、約１ミリ秒に限定されるものではなく、要求される振動の測定精度に応じて適宜変更可能である。さらには、ＦＦＴアナライザ１４による加速度信号、および加振力信号についての取り込み時間（サンプリング時間）も、４秒に限定されるものではなく、ゴルフクラブヘッド２０の特性、要求される測定精度に応じて適宜変更可能である。
【００２３】
ＰＣ１９は、処理装置であって、一般的なパーソナルコンピュータと同様の構成を有するものであり、ＣＰＵ（図示せず）、メモリ（図示せず）を備え、キーボード、マウスなどのコンピュータの入力に用いられる入力部（図示せず）と、入力部からの入力情報およびＣＰＵで情報処理された情報を表示するＬＣＤなどのモニタ１９ａとを有する。
【００２４】
ＰＣ１９は、ＦＦＴアナライザ１４からの加速度センサ１６による加速度信号の時系列データ（時間波形）、およびインパルスハンマ１２の加振力信号の時系列データ（時間波形）が入力されて、例えば、メモリ、ハードディスク（図示せず）等に記憶される。
また、ＰＣ１９には、更には後述するゴルフクラブヘッド２０のたわみ分布、たわみ分布の中心、フェース面２２上の重心位置等も記憶される。
さらには、ＰＣ１９は、ＦＦＴアナライザ１４からの加速度センサ１６による加速度信号の時系列データ、およびインパルスハンマ１２の加振力信号の時系列データ、周波数応答関数等、後述するゴルフクラブヘッド２０のたわみ分布、たわみ分布の中心、フェース面２２上の重心位置等をモニタ１９ａに、グラフ、数値、コンターマップ等の形態で表示させる機能を有する。
【００２５】
なお、ＰＣ１９は、ＦＦＴアナライザ１４を介することなく、加速度センサ１６による加速度信号をアンプ１８から直接、取り込むことができ、更にはインパクトハンマ１２の加振力信号を直接、取り込むことができる。この場合、加速度センサ１６による加速度信号、およびインパクトハンマ１２の加振力信号をＡＤ変換するＡＤ変換ボードをＰＣ１９に組み込み、アンプ１８と接続し、加速度信号、加振力信号をＰＣ１９に直接取り込む構成としてもよい。
【００２６】
本実施形態においては、測定点ｐは、図１（ｂ）に示すように格子の交点となる位置に予め設定されている。
ここで、ゴルフクラブヘッド２０のフェース面２２において図１（ｂ）に示すように、クラウンからヒールに向かう方向をＹ方向とする。このＹ方向と直交するトウからヒールに向かう方向をＸ方向とする。測定点ｐは、例えば、Ｙ方向に４行にわたって設定されている。各行においては、いずれも等間隔に測定点ｐが設定されており、２行目が１１点等間隔に、それ以外の１行目、３行目、４行目については１０点等間隔に設定されている。なお、本実施形態において、測定点ｐは、例えば、合計４１点設定されている。測定点ｐの設定数は４１点に限定されない。
【００２７】
また、図１（ｂ）に示すように、各測定点ｐに対して基準になる基準点Ｓが設けられている。この基準点Ｓは、第１行目のトウ側に設定されている。この基準点Ｓ（ｍ０、ｎ０）から間隔ｍで、Ｘ方向に測定点ｐが設定されており、Ｙ方向に基準点Ｓから間隔ｎで、測定点ｐが設定されている。例えば、基準点ＳのＹ方向のすぐ下の点の座標は、（ｍ０、ｎ１）である。
ゴルフクラブヘッド２０のフェース面２２に設定された各測定点ｐおよび基準点Ｓは、例えば、位置座標として、ＰＣ１９に記憶されている。
また、加速度センサ１６の取付位置Ｄは、Ｘ方向において最もヒール側に設定された測定点ｐからＸ方向側に１０ｍｍ離れた位置である。この加速度センサ１６の取付位置Ｄは、フェース面２２の各測定点ｐから外れる場所であれば、特に限定されるものではない。
【００２８】
次に、本実施形態のゴルフクラブヘッド２０のゴルフクラブヘッドのフェース面のたわみ分布測定方法について説明する。
まず、ゴルフクラブヘッド２０のフェース面２２の、例えば、第１行目のトウ側の測定点ｐを、インパクトハンマ１２で打撃する。この打撃時に、インパクトハンマ１２から加振力信号がＦＦＴアナライザ１４に出力される。測定点ｐの打撃によりゴルフクラブヘッド２０に生じた加速度が加速度センサ１６で測定されて、その加速度の加速度信号がＦＦＴアナライザ１４にアンプ１８を介して出力される。
【００２９】
ＦＦＴアナライザ１４において、加速度センサ１６による加速度信号、およびインパクトハンマ１２の加振力信号について、例えば、約１ミリ秒のサンプリング周期（４秒／４０９６間隔）で４秒間（サンプリング時間）取り込み、ゴルフクラブヘッド２０の加速度の時系列データおよびインパクトハンマ１２による測定点ｐの打撃時の加振力の時系列データが取得される。
【００３０】
次に、加速度の時系列データおよびインパクトハンマ１２による加振力の時系列データに基づいて、ＦＦＴアナライザ１４により測定点ｐにおける周波数応答関数が計算される。これらのゴルフクラブヘッド２０の加速度の時系列データおよびインパクトハンマ１２による加振力の時系列データならび周波数応答関数は、ＰＣ１９に出力され、ＰＣ１９のメモリ、ハードディスク等に記憶される。
【００３１】
周波数応答関数は、入力のフーリエスペクトルと出力のフーリエスペクトルとの比で表されるものである。入力のフーリエスペクトルは、インパクトハンマの加振力信号のフーリエスペクトルのことである。出力のフーリエスペクトルは、ゴルフクラブヘッドの加速度信号のフーリエスペクトルのことである。
本実施形態において、周波数応答関数をＨとするとき、この周波数応答関数Ｈは、下記数式１により表される。
【００３２】
Ｈ＝Ａ（ｆ）／Ｆ（ｆ）    （１）
【００３３】
なお、上記数式１において、Ｆ（ｆ）は、インパクトハンマ１２の加振力信号の複素フーリエスペクトルであり、Ａ（ｆ）は、ゴルフクラブヘッド２０の加速度の複素フーリエスペクトルである。
【００３４】
本実施形態において、周波数応答関数Ｈを求める場合、まず、ゴルフクラブヘッド２０の加速度の時系列データ、およびインパクトハンマ１２による加振力の時系列データについて、ＦＦＴアナライザ１４によりＦＦＴ（高速フーリエ変換）がなされる。これにより、ゴルフクラブヘッドの加速度についての複素フーリエスペクトル、すなわち、Ａ（ｆ）が得られる。更には、インパクトハンマ１２による加振力についての複素フーリエスペクトル、すなわち、Ｆ（ｆ）が得られる。
【００３５】
次に、ＦＦＴアナライザ１４において、加速度の複素フーリエスペクトルＡ（ｆ）および加振力の複素フーリエスペクトルＦ（ｆ）を用いて周波数応答関数Ｈが計算される。これにより、周波数応答関数Ｈが得られ、周波数応答関数ＨがＰＣ１９に出力されて記憶される。
【００３６】
次に、ＦＦＴアナライザ１４において、周波数応答関数Ｈから虚部が抽出されて、周波数応答関数Ｈの虚部と周波数応答関数Ｈの周波数との応答関数（以下、単に、応答関数という）が作成される。すなわち、周波数応答関数Ｈにおいて、虚部に関する関数が分離される。この応答関数もＰＣ１９に出力されて記憶される。例えば、応答関数は、図２（ａ）のグラフ３０に示す曲線３２のように表される。
【００３７】
次に、ＦＦＴアナライザ１４により、応答関数から１次共振周波数が求められる。これは、周波数応答関数の虚部がゴルフクラブヘッド２０のフェース面２２のモードシェイプを表すためである。
なお、高速フーリエ変換の処理自体は、公知の方法であり、市販のＦＦＴアナライザに内蔵されている機能を用いることができる。このため、その詳細な説明は省略する。
【００３８】
次に、図２（ａ）に示すグラフ３０の曲線３２（応答関数）において、周波数応答関数の虚部のレベルが最も高い周波数が１次共振周波数であることを利用して、ＦＦＴアナライザ１４により、図２（ａ）に示すグラフ３０に示す曲線３２において、周波数応答関数の虚部のレベルが最も高い周波数を抽出し、この抽出した周波数の値を１次共振周波数としてＰＣ１９に記憶させる。
【００３９】
なお、本実施形態においては、１次共振周波数を求める際、ＦＦＴアナライザ１４において、応答関数に対して、例えば、３０００Ｈｚ〜６０００Ｈｚの帯域の周波数分析がなされる。３０００Ｈｚ〜６０００Ｈｚの帯域において、周波数応答関数の虚部のピーク値が最も大きいものを１次共振周波数とする。
また、予めゴルフクラブヘッド２０について、ブランド別、商品別、またはゴルフクラブの番手、ドライバーもしくはスプーン等の種類別に１次共振周波数の発生帯域を、例えば、３２００〜５２００Ｈｚのように定めておいて、周波数応答関数の虚部の波形から、１次共振周波数を特定してもよい。
【００４０】
次に、ＦＦＴアナライザ１４により、図２（ｂ）に示すように、曲線３２（応答関数）を、１次共振周波数を中心周波数として所定の幅の周波数範囲に区画する。
本実施形態において、所定の幅の周波数範囲は、例えば、１／３オクターブバンドである。この場合、１次共振周波数をｆ１とするとき、１／３オクターブバンドで区画される周波数の最大値ｆmaxは下記数式２により、周波数の最小値ｆminは下記数式３により求めることができる。
【００４１】
ｆmax＝ｆ１・２（１／６）      （２）
【００４２】
ｆmin＝ｆ１・２（−１／６）      （３）
【００４３】
ここで、１次共振モードでは、フェース面２２のモードシェイプが太鼓状であるため、正の値をとる。そこで、本実施形態においては、図２（ｂ）に示す曲線３２が区画された周波数範囲において、負の値の領域３６を無視し、正の値の領域３４の合計値をたわみ量ｄとする。
そこで、次に、本実施形態においては、ＦＦＴアナライザ１４により、図２（ｂ）に示すように、例えば、１次共振周波数ｆ１を中心周波数として、１／３オクターブバンド幅に区画された周波数範囲において正の値の領域３４が合計された合計値が、たわみ量ｄとして求められる。
なお、本実施形態において、たわみ量を求めるための周波数範囲は、１／３オクターブバンド幅に限定されるものではない。
【００４４】
以上のようにして求めた測定点ｐのたわみ量ｄが、ＦＦＴアナライザ１４により、ＰＣ１９に出力されて記憶される。以降、順次、他の測定点ｐについてたわみ量ｄの測定を行い、合計４１個の測定点ｐについて、たわみ量ｄを求め、これらのたわみ量ｄがＰＣ１９に記憶される。
ＰＣ１９においては、各測定点ｐのたわみ量を、各測定点ｐに当てはめ、ゴルフクラブヘッド２０のフェース面２２のたわみ分布を得ることができる。
【００４５】
次に、ゴルフクラブヘッド２０のフェース面２２上のたわみ分布の中心Ｃを各測定点ｐのたわみ量ｄを用いて求める。本実施形態では、直交するＸ方向およびＹ方向におけるたわみ分布の中心Ｃを求める。本実施形態においては、例えば、基準点Ｓからの距離の加重平均を、たわみ分布の中心Ｃとする。
【００４６】
たわみ分布の中心Ｃにおいて、Ｘ方向の中心座標をＸＣとし、Ｙ方向の中心座標をＹＣとするとき、Ｘ方向の中心座標ＸＣは、下記数式４により表すことができる。また、Ｙ方向の中心座標をＹＣは、下記数式５により表すことができる。
なお、下記数式４、５において、ｍｉは測定点ｐのＸ方向における座標、ｎｉは測定点のＹ方向における座標であり、ｄｉは各測定点ｐでのたわみ量である。
【００４７】
【数１】
2011050440.tif 000003
【００４８】
【数２】
2011050440.tif 000004
【００４９】
以上のようにして、本実施形態においては、たわみ分布およびたわみ分布の中心Ｃを求めることができる。このたわみ分布の中心Ｃは、各測定点ｐのたわみ量ｄとともにＰＣ１９に記憶される。
【００５０】
本実施形態においては、周波数応答関数の虚部と周波数応答関数の周波数との応答関数から１次共振周波数ｆ１を求め、この１次共振周波数ｆ１を中心周波数として区画された周波数範囲における正の値を合計して、この正の値の合計値を測定点ｐにおけるたわみ量ｄとしている。これにより、たわみ量ｄが１次共振モードに応じたものとなるため、ゴルフクラブヘッド２０のフェース面２２のたわみ分布を高い精度で求めることができる。これにより、飛距離、スイートスポットの大きさ等のゴルフクラブヘッド２０の性能をより高い精度で評価することができる。
【００５１】
更には、本実施形態においては、各測定点ｐのたわみ量ｄを用いて、フェース面２２上におけるたわみ分布の中心Ｃを求めることができる。この場合、各測定点ｐのたわみ量ｄを高い精度で求めることができるため、たわみ分布の中心Ｃも高い精度で求めることができる。これによっても、スイートスポットの大きさ等のゴルフクラブヘッドの性能の評価をより高い精度で行うことができる。
また、インパクトハンマ１２で、フェース面２２の測定点ｐを打撃するだけなので、たわみ分布およびたわみ分布の中心Ｃを容易に、しかも、短時間で求めることができる。
【００５２】
なお、本実施形態においては、ＰＣ１９に、ゴルフクラブヘッド２０について得られたたわみ量ｄについて、所定の間隔で区画し、各区画について色を割り当てるか、または諧調を割り当てておく。これにより、ＰＣ１９により、各測定点ｐのたわみ量ｄに応じた色で区分けされるか、または諧調で区分けされたコンターマップを作成することができ、このコンターマップをモニタ１９ａに表示することができる。この場合、たわみ分布の中心Ｃも、コンターマップに表示することができる。
本実施形態においては、例えば、図３（ａ）〜（ｅ）に示すようなコンターマップが得られ、モニタ１９ａに表示される。これにより、視覚的にたわみ分布を知覚することができる。
【００５３】
なお、図３（ａ）〜（ｅ）に示すゴルフクラブヘッドのコンターマップにおいては、たわみ分布の中心Ｃおよび重心位置Ｇも示している。この重心位置Ｇは、ゴルフクラブヘッドのフェース面の重心位置を示すものである。重心位置Ｇは、例えば、特開２００８−２２０６６５号公報に記載されている重心高さＦＧＨの測定方法等の公知の方法により求めることができる。
【００５４】
なお、本実施形態においては、打撃されたときに加速度を測定しているが、応答信号として検出することができれば、加速度に限定されるものではない。例えば、加速度に代えて、歪ゲージを用いて測定された歪量を表す信号（歪信号）であっても、音圧計を用いて測定された音圧レベルを表す信号（音圧信号）であってもよい。なお、歪ゲージを用いる場合には、打撃による力が歪ゲージに直接的に作用しない取付位置とする。
【００５５】
本発明は、基本的に以上のようなものである。以上、本発明のゴルフクラブヘッドのフェース面のたわみ分布測定方法について詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良または変更をしてもよいのはもちろんである。
【符号の説明】
【００５６】
１０  測定装置
１２  インパルスハンマ
１４  ＦＦＴアナライザ
１６  加速度センサ
１８  アンプ
１９  コンピュータ（ＰＣ）
２０  ゴルフクラブヘッド
２２  フェース面
２４  ホーゼル部
３０  グラフ
【図１】
2011050440.tif 000005
【図２】
2011050440.tif 000006
【図３】
2011050440.tif 000007