JP5229073B2 - Method for evaluating rod-shaped body and system for evaluating rod-shaped body - Google Patents

Method for evaluating rod-shaped body and system for evaluating rod-shaped body Download PDF

Info

Publication number
JP5229073B2
JP5229073B2 JP2009090243A JP2009090243A JP5229073B2 JP 5229073 B2 JP5229073 B2 JP 5229073B2 JP 2009090243 A JP2009090243 A JP 2009090243A JP 2009090243 A JP2009090243 A JP 2009090243A JP 5229073 B2 JP5229073 B2 JP 5229073B2
Authority
JP
Japan
Prior art keywords
vibration
peak
value
absolute value
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2009090243A
Other languages
Japanese (ja)
Other versions
JP2010243251A (en
Inventor
一則 小野
真生 山本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yokohama Rubber Co Ltd
Original Assignee
Yokohama Rubber Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yokohama Rubber Co Ltd filed Critical Yokohama Rubber Co Ltd
Priority to JP2009090243A priority Critical patent/JP5229073B2/en
Publication of JP2010243251A publication Critical patent/JP2010243251A/en
Application granted granted Critical
Publication of JP5229073B2 publication Critical patent/JP5229073B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

本発明は、棒状体の振動減衰性能等の振動特性を評価する評価方法およびその棒状体の評価システムに関し、特に、ゴルフクラブシャフト、釣り竿、アイスホッケーのスティック、グランドホッケーのスティック等について振動減衰性能等の振動特性を評価する評価方法および棒状体の評価システムに関する。   The present invention relates to an evaluation method for evaluating vibration characteristics such as vibration damping performance of a rod-shaped body, and an evaluation system for the rod-shaped body, and in particular, vibration damping performance for golf club shafts, fishing rods, ice hockey sticks, ground hockey sticks, and the like. The present invention relates to an evaluation method for evaluating vibration characteristics such as a rod-like body.

現在、ゴルフクラブシャフトとして、金属により形成された金属製のゴルフクラブシャフトと、繊維強化複合材（ＦＲＰ）で形成されたＦＲＰ製のゴルフクラブシャフトがある。金属製のゴルフクラブシャフトのうち、特に、スチール製のものはスチールシャフトと呼ばれ、ＦＲＰ製のゴルフクラブシャフトのうち、カーボン繊維強化複合材（ＣＦＲＰ）で形成されたゴルフクラブシャフトはカーボンシャフトと呼ばれている。
このＦＲＰ製のゴルフクラブシャフトは、例えば、シートワインディング法、フィラメントワインディング法、３軸ブレーディング法等の製造方法によって製造されている。
ＦＲＰ製のゴルフクラブシャフトにおいては、シャフトの周方向における変形特性が異なることが知られている。この変形特性の異方性は、上記製造方法、またはもともとＦＲＰ自体に異方性があることによるものである。
また、ＦＲＰ製のゴルフクラブシャフトは、強化繊維の配向方向、積層数などを変えることにより、その長手方向の剛性分布、振動減衰等の振動特性を調整することができる。 Currently, as golf club shafts, there are metal golf club shafts made of metal and FRP golf club shafts made of fiber reinforced composite material (FRP). Of the golf club shafts made of metal, in particular, those made of steel are called steel shafts. Among golf club shafts made of FRP, golf club shafts made of carbon fiber reinforced composite (CFRP) are called carbon shafts. It is.
This FRP golf club shaft is manufactured by a manufacturing method such as a sheet winding method, a filament winding method, or a triaxial braiding method.
It is known that the FRP golf club shaft has different deformation characteristics in the circumferential direction of the shaft. The anisotropy of the deformation characteristics is due to the above-described manufacturing method or the fact that the FRP itself has anisotropy.
In addition, by changing the orientation direction of the reinforcing fibers, the number of laminated layers, etc., the FRP golf club shaft can adjust the vibration characteristics such as rigidity distribution and vibration attenuation in the longitudinal direction.

ＦＲＰ製のゴルフクラブシャフトについて、振動特性が設計された通りになっているかを確認するために振動特性を評価することがなされる。この振動特性の評価のために、ゴルフクラブシャフトの一端を固定して、他端を上下に振動させると、上述のシャフトの周方向における変形特性の異方性により、ゴルフクラブシャフトが左右にも振動してしまい、振動特性を正確に評価できないという問題点がある。更には、製造されたＦＲＰ製のゴルフクラブシャフトの振動減衰等の振動特性についての評価方法に特に決まったものがないのが現状である。   For the FRP golf club shaft, the vibration characteristics are evaluated in order to confirm whether the vibration characteristics are as designed. In order to evaluate this vibration characteristic, when one end of the golf club shaft is fixed and the other end is vibrated up and down, the golf club shaft is also moved left and right due to the anisotropy of the deformation characteristic in the circumferential direction of the shaft described above. There is a problem that the vibration characteristics are not accurately evaluated. Furthermore, there is currently no specific evaluation method for vibration characteristics such as vibration attenuation of the manufactured FRP golf club shaft.

本発明の目的は、前記従来技術に基づく問題点を解消し、振動減衰を容易にかつ高い精度で評価できる棒状体の評価方法および棒状体の評価システムを提供することにある。   An object of the present invention is to provide a rod-like body evaluation method and a rod-like body evaluation system that can solve the problems based on the prior art and can easily evaluate vibration damping with high accuracy.

上記目的を達成するために、本発明は、棒状体の一端部を固定し、他端側を振動させ、前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測し、上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得する工程と、前記計測された上下方向の振動の時系列データおよび前記左右方向の振動の時系列データについて、それぞれ、高速フーリエ変換し、振動の周波数に対する振動の振幅レベルを示す振幅特性データと、前記振動の周波数に対する振動の位相を示す位相特性データをそれぞれ作成する工程と、前記上下方向の振動と前記左右方向の振動に対して、前記各振幅特性データから、それぞれ少なくとも１次共振周波数を求める工程と、前記求められた各１次共振周波数を中心周波数としたバンドパスフィルタを設定する工程と、前記各バンドパスフィルタを、それぞれ前記各振幅特性データにかけた後、それぞれ逆高速フーリエ変換し、１次共振周波数における時間波形データを、それぞれ求める工程と、前記１次共振周波数の時間波形データについて、山側と谷側のピークのピーク時間およびピークレベルを求める工程と、前記上下方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＵｐとし、前記左右方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出する工程と、前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求める工程と、前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する工程とを有することを特徴とする棒状体の評価方法を提供するものである。
本発明において、上下方向および左右方向の振動の時系列データは、例えば、上下方向および左右方向の加速度の時系列データである。 In order to achieve the above object, the present invention fixes one end of a rod-shaped body, vibrates the other end, and simultaneously measures the vertical vibration and the horizontal vibration of the other end for a predetermined time. The step of acquiring time series data of vibration in the vertical direction and time series data of vibration in the left and right direction, and the time series data of the measured vertical vibration and the time series data of the left and right vibration, respectively, A step of performing fast Fourier transform to generate amplitude characteristic data indicating a vibration amplitude level with respect to a vibration frequency, and phase characteristic data indicating a vibration phase with respect to the vibration frequency; and With respect to vibration, a step of obtaining at least a primary resonance frequency from each amplitude characteristic data, and a bandpass having the obtained primary resonance frequency as a center frequency. A step of setting a filter, a step of applying each of the bandpass filters to each of the amplitude characteristic data, and then performing inverse fast Fourier transform, respectively, to obtain time waveform data at a primary resonance frequency, and the primary resonance For the time waveform data of the frequency, the step of obtaining the peak time and peak level of the peak on the peak side and the valley side, the absolute value of the peak value in the time waveform data on the primary resonance frequency in the vertical direction is Up, and the horizontal direction The peak value in the left-right direction of the peak time closest to the absolute value Up of each peak value with Lp as the absolute value of the peak value in the time waveform data of the primary resonance frequency Extracting the absolute value Lp of the peak value, the absolute value Up of the peak value, and the absolute value Lp of the extracted peak value Evaluation of the rod-shaped body, comprising: a step of obtaining a square root value Vp of the first and second steps; and a step of calculating an attenuation factor using the obtained value Vp and calculating a loss factor based on the attenuation factor. A method is provided.
In the present invention, the time series data of vertical and horizontal vibrations is, for example, time series data of vertical and horizontal accelerations.

また、本発明は、棒状体の一端部を固定し、他端側を振動させ、前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測し、上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得する工程と、前記上下方向の振動の時系列データおよび前記左右方向の振動の時系列データについて、それぞれ山側と谷側のピークのピーク時間およびピークレベルを求める工程と、前記上下方向の振動の時系列データにおけるピーク値の絶対値をＵｐとし、前記左右方向の振動の時系列データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出する工程と、前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求める工程と、前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する工程とを有することを特徴とする棒状体の評価方法を提供するものである。
本発明において、上下方向および左右方向の振動の時系列データは、例えば、上下方向および左右方向の変位の時系列データである。 The present invention also fixes one end of the rod-like body, vibrates the other end, and simultaneously measures the vertical vibration and the horizontal vibration of the other end for a predetermined time, A step of obtaining time series data and time series data of left and right vibrations, and peak time of peak and valley sides respectively for the time series data of vertical vibrations and time series data of left and right vibrations, and A step of obtaining a peak level; an absolute value of a peak value in the time-series data of the vertical vibrations is Up; an absolute value of the peak value in the time-series data of the horizontal vibrations is Lp; Using the absolute value Up as a reference, a step of extracting the absolute value Lp of the peak value in the left-right direction at the peak time closest to the absolute value Up of each peak value, and the absolute value Up of the peak value and the previous A step of obtaining a square root value Vp of the absolute value Lp of the extracted peak value, a step of calculating an attenuation factor using the obtained value Vp, and calculating a loss factor based on the attenuation factor. It provides the evaluation method of the rod-shaped body characterized by having.
In the present invention, the time-series data of vertical and horizontal vibrations is, for example, time-series data of vertical and horizontal displacements.

本発明において、前記各振幅特性データから求める共振周波数は、１次共振周波数、２次共振周波数および３次共振周波数であることが好ましい。
また、本発明において、前記バンドパスフィルタは、１／３オクターブバンドパスフィルタであることが好ましい。
さらに、本発明において、前記減衰率は、例えば、対数減衰率である。
さらにまた、本発明において、前記棒状体は、例えば、ゴルフクラブシャフトである。 In the present invention, the resonance frequency obtained from each amplitude characteristic data is preferably a primary resonance frequency, a secondary resonance frequency, and a tertiary resonance frequency.
In the present invention, the bandpass filter is preferably a 1/3 octave bandpass filter.
Furthermore, in the present invention, the attenuation rate is, for example, a logarithmic attenuation rate.
Furthermore, in the present invention, the rod-shaped body is, for example, a golf club shaft.

本発明は、棒状体の一端部を固定し、他端側を振動させたときの前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測する計測装置と、前記計測装置で計測された上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得するデータ取得装置と、前記計測された上下方向の振動の時系列データおよび前記左右方向の振動の時系列データについて、それぞれ、高速フーリエ変換し、振動の周波数に対する振動の振幅レベルを示す振幅特性データと、前記振動の周波数に対する振動の位相を示す位相特性データをそれぞれ作成し、前記上下方向の振動と前記左右方向の振動に対して、前記各振幅特性データから、それぞれ少なくとも１次共振周波数を求め、前記求められた各１次共振周波数を中心周波数としたバンドパスフィルタを設定し、前記各バンドパスフィルタを、それぞれ前記各振幅特性データにかけた後、それぞれ逆高速フーリエ変換し、１次共振周波数における時間波形データを、それぞれ求め、前記１次共振周波数の時間波形データについて、山側と谷側のピークのピーク時間およびピークレベルを求め、前記上下方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＵｐとし、前記左右方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出し、前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求め、前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する処理装置とを有することを特徴とする棒状体の評価システムを提供するものである。   The present invention provides a measuring device for measuring a vertical vibration and a horizontal vibration of the other end at the same time for a predetermined time when one end of a rod-like body is fixed and the other end is vibrated, and the measurement A data acquisition device for acquiring time series data of vertical vibration and time series data of horizontal vibration measured by the apparatus, time series data of the vertical vibration measured and time of the horizontal vibration The series data is subjected to fast Fourier transform, respectively, and amplitude characteristic data indicating the vibration amplitude level with respect to the vibration frequency, and phase characteristic data indicating the vibration phase with respect to the vibration frequency are respectively created, and the vertical vibration and For each vibration in the left-right direction, at least a primary resonance frequency is obtained from each amplitude characteristic data, and each of the obtained primary resonance frequencies is set as a center frequency. A band pass filter is set, and each band pass filter is applied to each amplitude characteristic data, and then inverse fast Fourier transform is performed to obtain time waveform data at the primary resonance frequency. For the time waveform data, the peak time and peak level of peaks on the peak side and the valley side are obtained, the absolute value of the peak value in the time waveform data of the primary resonance frequency in the vertical direction is set as Up, and the primary value in the horizontal direction is set. The absolute value Lp of the peak value in the left and right direction of the peak time closest to the absolute value Up of each peak value, with the absolute value Up of each peak value as a reference, in the absolute value Up of the peak value in the time waveform data of the resonance frequency , And the square root value Vp of the absolute value Up of the peak value and the absolute value Lp of the extracted peak value. Calculated, to calculate the attenuation rate by using the obtained value Vp, there is provided an evaluation system for rod-like body; and a processor for calculating a loss factor, based on the attenuation factor.

本発明においては、前記計測装置は、前記他端部の上下方向の加速度および左右方向の加速度を計測する加速度センサを有することが好ましい。   In the present invention, it is preferable that the measurement apparatus includes an acceleration sensor that measures the vertical acceleration and the horizontal acceleration of the other end.

本発明は、棒状体の一端部を固定し、他端側を振動させたときの前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測する計測装置と、前記計測装置で計測された上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得するデータ取得装置と、前記上下方向の振動の時系列データおよび前記左右方向の振動の時系列データについて、それぞれ山側と谷側のピークのピーク時間およびピークレベルを求め、前記上下方向の振動の時系列データにおけるピーク値の絶対値をＵｐとし、前記左右方向の振動の時系列データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出し、前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求め、前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する処理装置とを有することを特徴とする棒状体の評価システムを提供するものである。
本発明においては、前記計測装置は、前記他端部の上下方向の変位および左右方向の変位を計測するレーザ変位計を有することが好ましい。 The present invention provides a measuring device for measuring a vertical vibration and a horizontal vibration of the other end at the same time for a predetermined time when one end of a rod-like body is fixed and the other end is vibrated, and the measurement Data acquisition device for acquiring time series data of vertical vibration and time series data of horizontal vibration measured by the apparatus, time series data of vertical vibration and time series data of left and right vibration The peak time and peak level of the peak on the mountain side and the valley side are obtained, the absolute value of the peak value in the time series data of the vertical vibration is taken as Up, and the absolute value of the peak value in the time series data of the horizontal vibration Using the absolute value Up of each peak value as a reference, the absolute value Lp of the peak value in the left-right direction of the peak time closest to the absolute value Up of each peak value is extracted. A square root value Vp of the absolute value Up of the peak value and the absolute value Lp of the extracted peak value is obtained, an attenuation rate is calculated using the obtained value Vp, and a loss is calculated based on the attenuation rate. The present invention provides a rod-like body evaluation system comprising a processing device for calculating a coefficient.
In this invention, it is preferable that the said measuring apparatus has a laser displacement meter which measures the displacement of the up-down direction of the said other end part, and the displacement of the left-right direction.

本発明によれば、振動減衰を容易に評価することができる。しかも、棒状体を上下方向に振動させた場合、左右方向にも振動してしまうものであっても、高い精度で振動減衰を評価することができる。
このため、本発明によれば、例えば、ゴルフクラブシャフトの振動減衰の評価を、容易かつ高い精度でできる。 According to the present invention, vibration damping can be easily evaluated. Moreover, when the rod-like body is vibrated in the vertical direction, vibration attenuation can be evaluated with high accuracy even if it vibrates in the horizontal direction.
For this reason, according to the present invention, for example, evaluation of vibration damping of a golf club shaft can be performed easily and with high accuracy.

（ａ）は、本発明の第１の実施形態の棒状体の評価方法に用いられる評価システムを示す模式図であり、（ｂ）は、図１（ａ）に示す評価システムにおける第１の加速度センサ、および第２の加速度センサの取り付け位置を拡大して示す模式図である。(A) is a schematic diagram which shows the evaluation system used for the evaluation method of the rod-shaped body of the 1st Embodiment of this invention, (b) is the 1st acceleration in the evaluation system shown to Fig.1 (a). It is a schematic diagram which expands and shows the attachment position of a sensor and a 2nd acceleration sensor. 本発明の第１の実施形態の棒状体の評価方法によるゴルフクラブシャフトの評価方法を工程順に示すフローチャートである。It is a flowchart which shows the evaluation method of the golf club shaft by the evaluation method of the rod-shaped body of the 1st Embodiment of this invention in order of a process. （ａ）は、縦軸に加速度をとり、横軸に時間をとって、ゴルフクラブシャフトの上下方向の加速度の時系列データを示すグラフであり、（ｂ）は、縦軸に加速度をとり、横軸に時間をとって、ゴルフクラブシャフトの左右方向の加速度の時系列データを示すグラフである。(A) is a graph showing time-series data of acceleration in the vertical direction of the golf club shaft, with acceleration on the vertical axis and time on the horizontal axis, and (b) takes acceleration on the vertical axis, It is a graph which shows the time series data of the acceleration of the horizontal direction of a golf club shaft, taking time on a horizontal axis. 縦軸にパワースペクトルをとり、横軸に周波数をとって、振幅特性データを示すグラフである。It is a graph which shows an amplitude characteristic data by taking a power spectrum on a vertical axis | shaft and taking a frequency on a horizontal axis. 縦軸にフィルタ値をとり、横軸に周波数をとって、本実施形態で用いたバンドパスフィルタを説明するためのグラフである。It is a graph for demonstrating the band pass filter used by this embodiment, taking a filter value on a vertical axis | shaft and taking a frequency on a horizontal axis. 縦軸に振動ピークレベルをとり、横軸に時間をとって、１次共振周波数における減衰率を説明するためのグラフである。It is a graph for demonstrating the damping factor in a primary resonant frequency, taking a vibration peak level on a vertical axis | shaft and taking time on a horizontal axis. 縦軸に振動ピークレベルをとり、横軸に時間をとって、２次共振周波数における減衰率を説明するためのグラフである。It is a graph for explaining the attenuation rate at the secondary resonance frequency with the vertical axis representing the vibration peak level and the horizontal axis representing time. 縦軸に振動ピークレベルをとり、横軸に時間をとって、３次共振周波数における減衰率を説明するためのグラフである。It is a graph for explaining the attenuation factor at the third resonance frequency with the vertical axis representing the vibration peak level and the horizontal axis representing time. 縦軸に損失係数をとり、横軸に周波数をとって、各共振周波数における損失係数の結果を示すグラフである。It is a graph which shows the result of the loss coefficient in each resonance frequency, taking a loss coefficient on a vertical axis | shaft and taking a frequency on a horizontal axis. 本発明の第２の実施形態の棒状体の評価方法に用いられる評価システムを示す模式的斜視図である。It is a typical perspective view which shows the evaluation system used for the evaluation method of the rod-shaped body of the 2nd Embodiment of this invention. 本発明の第２の実施形態の棒状体の評価方法によるゴルフクラブシャフトの評価方法を工程順に示すフローチャートである。It is a flowchart which shows the evaluation method of the golf club shaft by the evaluation method of the rod-shaped body of the 2nd Embodiment of this invention in order of a process.

以下に、添付の図面に示す好適実施形態に基づいて、本発明の棒状体の評価方法および棒状体の評価システムを詳細に説明する。
図１は、（ａ）は、本発明の第１の実施形態の棒状体の評価方法に用いられる評価システムを示す模式図であり、（ｂ）は、図１（ａ）に示す評価システムにおける第１の加速度センサ、および第２の加速度センサの取り付け位置を拡大して示す模式図である。 The rod-like body evaluation method and rod-like body evaluation system of the present invention will be described below in detail based on preferred embodiments shown in the accompanying drawings.
1A is a schematic diagram showing an evaluation system used in the rod-like body evaluation method according to the first embodiment of the present invention, and FIG. 1B is a schematic diagram of the evaluation system shown in FIG. It is a schematic diagram which expands and shows the attachment position of a 1st acceleration sensor and a 2nd acceleration sensor.

図１に示す評価システム１０は、棒状体として、例えば、ゴルフクラブシャフト３０の減衰性能などの振動特性を評価するものである。
図１に示すように、評価システム１０は、ゴルフクラブシャフト３０を固定する固定部１２と、ゴルフクラブシャフト２２の振動を検出する第１の加速度センサ１４および第２の加速度センサ１６と、チャージアンプ１８と、ＦＦＴアナライザ２０と、コンピュータ（以下、ＰＣという）２２とを有する。
また、第１の加速度センサ１４および第２の加速度センサ１６は、チャージアンプ１８を介してＦＦＴアナライザ２０に接続されている。ＦＦＴアナライザ２０は、ＰＣ２２により制御される。 An evaluation system 10 shown in FIG. 1 evaluates vibration characteristics such as a damping performance of a golf club shaft 30 as a rod-like body.
As shown in FIG. 1, the evaluation system 10 includes a fixing portion 12 that fixes a golf club shaft 30, a first acceleration sensor 14 and a second acceleration sensor 16 that detect vibrations of the golf club shaft 22, and a charge amplifier. 18, an FFT analyzer 20, and a computer (hereinafter referred to as PC) 22.
The first acceleration sensor 14 and the second acceleration sensor 16 are connected to the FFT analyzer 20 via the charge amplifier 18. The FFT analyzer 20 is controlled by the PC 22.

固定部１２は、ゴルフクラブシャフト３０の後端部３２ａを固定するものである。この固定部１２は、例えば、水平面Ｂの上に設けられている。固定部１２にゴルフクラブシャフト３０が固定された場合、ゴルフクラブシャフト３０の中心軸Ｃが水平面Ｂと平行になる。この固定部１２におけるゴルフクラブシャフト３０の後端部３２ａの固定長さＬ２は、例えば、１７０ｍｍである。
なお、固定部１２のゴルフクラブシャフト３０の固定方法は、特に限定されるものではない。 The fixing portion 12 fixes the rear end portion 32 a of the golf club shaft 30. The fixing portion 12 is provided on the horizontal plane B, for example. When the golf club shaft 30 is fixed to the fixing portion 12, the center axis C of the golf club shaft 30 is parallel to the horizontal plane B. Fixed length L 2 of the rear end portion 32a of the golf club shaft 30 in the fixing unit 12 is, for example, 170 mm.
In addition, the fixing method of the golf club shaft 30 of the fixing | fixed part 12 is not specifically limited.

ここで、ゴルフクラブシャフト３０の後端部３２ａとは、ゴルフクラブシャフト３０をゴルフクラブヘッド（図示せず）に取り付け、ゴルフクラブ（図示せず）とする場合、グリップ（図示せず）が取り付けられる側の端部のことである。
また、ゴルフクラブシャフト３０の先端部３２ｂとは、ゴルフクラブヘッド（図示せず）を取り付ける側の端部のことである。 Here, when the golf club shaft 30 is attached to a golf club head (not shown) and used as a golf club (not shown), a grip (not shown) is attached to the rear end portion 32a of the golf club shaft 30. It is the end of the side to be used.
Further, the tip portion 32b of the golf club shaft 30 is an end portion on the side where a golf club head (not shown) is attached.

本実施形態においては、ゴルフクラブシャフト３０は、単体であっても、ゴルフクラブヘッドが取り付けられた状態であってもよい。ゴルフクラブヘッドが取り付けられている場合、グリップ側の端部を、例えば、固定長さ１７０ｍｍで固定する。
また、本実施形態においては、後述するように、ゴルフクラブシャフト３０の評価には、例えば、評価するゴルフクラブシャフト３０の後端部３２ａを固定し、先端部３２ｂを振動させて、ゴルフクラブシャフト３０の振動特性を測定して評価する。 In the present embodiment, the golf club shaft 30 may be a single body or a state in which a golf club head is attached. When the golf club head is attached, the grip side end is fixed at a fixed length of 170 mm, for example.
In the present embodiment, as will be described later, for evaluation of the golf club shaft 30, for example, the golf club shaft 30 is evaluated by fixing the rear end portion 32a of the golf club shaft 30 to be evaluated and vibrating the front end portion 32b. 30 vibration characteristics are measured and evaluated.

第１の加速度センサ１４は、ゴルフクラブシャフト２２の長手方向Ｌと直交する上下方向Ｖにおけるゴルフクラブシャフト２２の先端部２４ｂの振動を測定するものである。なお、上下方向Ｖは、水平面Ｂに対して垂直な方向である。
この第１の加速度センサ１４は、例えば、ゴルフクラブシャフト２２の先端２２ａから後端部２４ａ側に距離Ｌ１、離れた位置で、かつゴルフクラブシャフト２２の周面２２ｂにおいて、水平面Ｂから上記上下方向Ｖにおける距離が最も離れた位置に設けられている。換言すれば、第１の加速度センサ１４は、ゴルフクラブシャフト２２の先端２２ａから後端部２４ａ側にＬ１離れた位置で、かつゴルフクラブシャフト２２の周面２２ｂにおいて、ゴルフクラブシャフト２２の中心軸Ｃを通り上記上下方向Ｖと平行な線Ｙ上に設けられている。この第１の加速度センサ１４が設けられる位置を第１の取付位置αという。この距離Ｌ１は、例えば、１００ｍｍである。 The first acceleration sensor 14 measures the vibration of the tip portion 24 b of the golf club shaft 22 in the vertical direction V perpendicular to the longitudinal direction L of the golf club shaft 22. The vertical direction V is a direction perpendicular to the horizontal plane B.
For example, the first acceleration sensor 14 is located above the horizontal plane B at a distance L 1 away from the front end 22a of the golf club shaft 22 toward the rear end portion 24a and on the peripheral surface 22b of the golf club shaft 22. The distance in the direction V is provided at the farthest position. In other words, the first acceleration sensor 14 is located at the center of the golf club shaft 22 at a position L 1 away from the front end 22 a of the golf club shaft 22 toward the rear end portion 24 a and on the peripheral surface 22 b of the golf club shaft 22. It is provided on a line Y passing through the axis C and parallel to the vertical direction V. The position where the first acceleration sensor 14 is provided is referred to as a first attachment position α. The distance L 1 is, for example, 100 mm.

第２の加速度センサ１６は、ゴルフクラブシャフト２２の長手方向Ｌおよび上下方向Ｖと直交する左右方向Ｈにおけるゴルフクラブシャフト２２の先端部２４ｂの振動を測定するものである。
この第２の加速度センサ１６は、第１の加速度センサ１４を、ゴルフクラブシャフト２２の中心軸Ｃを回転中心にして９０°回転させた位置に設けられている。換言すれば、第２の加速度センサ１６は、ゴルフクラブシャフト２２の先端２２ａから後端部２４ａ側にＬ１離れた位置で、かつゴルフクラブシャフト２２の周面２２ｂにおいて、ゴルフクラブシャフト２２の中心軸Ｃを通り長手方向Ｌおよび上下方向Ｖと直交する左右方向Ｈと平行な線Ｘ上に設けられている。この第２の加速度センサ１６が設けられる位置を第２の取付位置βという。
第１の加速度センサ１４および第２の加速度センサ１６には、例えば、軽量加速度ピックアップが用いられる。この軽量加速度ピックアップは、加速度信号を電荷の形態で、出力する。なお、加速度ピックアップは超軽量のものが好ましい。例えば、質量が０．６ｇのものを用いるとよい。これは、質量が大きいと共振周波数への影響が大きくなるためである。 The second acceleration sensor 16 measures the vibration of the front end portion 24 b of the golf club shaft 22 in the left-right direction H perpendicular to the longitudinal direction L and the up-down direction V of the golf club shaft 22.
The second acceleration sensor 16 is provided at a position obtained by rotating the first acceleration sensor 14 by 90 ° about the central axis C of the golf club shaft 22. In other words, the second acceleration sensor 16 is located at the center of the golf club shaft 22 at a position L 1 away from the front end 22 a of the golf club shaft 22 toward the rear end portion 24 a and on the peripheral surface 22 b of the golf club shaft 22. It is provided on a line X parallel to the horizontal direction H passing through the axis C and orthogonal to the longitudinal direction L and the vertical direction V. The position where the second acceleration sensor 16 is provided is referred to as a second attachment position β.
As the first acceleration sensor 14 and the second acceleration sensor 16, for example, a lightweight acceleration pickup is used. This lightweight acceleration pickup outputs an acceleration signal in the form of electric charges. The acceleration pickup is preferably an ultralight one. For example, a material having a mass of 0.6 g may be used. This is because a large mass increases the influence on the resonance frequency.

本実施形態においては、ゴルフクラブシャフト２２の先端２２ａに、水平面Ｂに向う下方向ＶＬに、例えば、２０ｍｍの変位を与えてゴルフクラブシャフト２２を自由減衰振動させる。このとき、第１の加速度センサ１４により、上下方向Ｖのゴルフクラブシャフト２２の自由減衰振動が加速度として測定され、第２の加速度センサ１６により、左右方向Ｈのゴルフクラブシャフト２２の自由減衰振動が加速度として測定される。 In the present embodiment, the golf club shaft 22 is subjected to free-damping vibration by applying a displacement of, for example, 20 mm to the front end 22a of the golf club shaft 22 in the downward direction VL toward the horizontal plane B. At this time, the first acceleration sensor 14 measures free damping vibration of the golf club shaft 22 in the vertical direction V as acceleration, and the second acceleration sensor 16 causes free damping vibration of the golf club shaft 22 in the left-right direction H. Measured as acceleration.

チャージアンプ１８は、第１の加速度センサ１４および第２の加速度センサ１６に接続されている。このチャージアンプ１８には、第１の加速度センサ１４および第２の加速度センサ１６のそれぞれから得られた加速度信号が、例えば、電荷の形態で入力される。チャージアンプ１８は、その電荷の形態で入力された各加速度信号を、それぞれ電圧に変換してＦＦＴアナライザ２０に加速度信号として出力する。   The charge amplifier 18 is connected to the first acceleration sensor 14 and the second acceleration sensor 16. An acceleration signal obtained from each of the first acceleration sensor 14 and the second acceleration sensor 16 is input to the charge amplifier 18 in the form of electric charge, for example. The charge amplifier 18 converts each acceleration signal input in the form of the charges into a voltage and outputs the voltage to the FFT analyzer 20 as an acceleration signal.

ＦＦＴアナライザ２０は、チャージアンプ１８に接続されており、このチャージアンプ１８から電圧の形態で出力される第１の加速度センサ１４の加速度信号および第２の加速度センサ１６の加速度信号を時系列で取り込み、ＰＣ２２に出力するものである。
このＦＦＴアナライザ２０においては、例えば、第１の加速度センサ１４の加速度信号、および第２の加速度センサ１６の加速度信号について、それぞれ振動を開始してから例えば、約１ミリ秒のサンプリング周期（４秒／４０９６間隔）で４秒間（サンプリング時間）取り込む。ＦＦＴアナライザ２０により、第１の加速度センサ１４による上下方向Ｖの加速度の時系列データ（時間波形）が得られ、第２の加速度センサ１６による左右方向Ｈの加速度の時系列データ（時間波形）が得られる。 The FFT analyzer 20 is connected to the charge amplifier 18 and takes in time series the acceleration signal of the first acceleration sensor 14 and the acceleration signal of the second acceleration sensor 16 output from the charge amplifier 18 in the form of voltage. Are output to the PC 22.
In the FFT analyzer 20, for example, the acceleration signal of the first acceleration sensor 14 and the acceleration signal of the second acceleration sensor 16 start to vibrate, for example, about 1 millisecond sampling period (4 seconds). / 4096 interval) for 4 seconds (sampling time). The FFT analyzer 20 obtains time series data (time waveform) of acceleration in the vertical direction V by the first acceleration sensor 14, and time series data (time waveform) of acceleration in the left and right direction H by the second acceleration sensor 16. can get.

なお、ＦＦＴアナライザ２０による加速度信号についてのサンプリング周期は、約１ミリ秒に限定されるものではなく、要求される振動の測定精度に応じて適宜変更可能である。さらには、ＦＦＴアナライザ２０による変位についての取り込み時間（サンプリング時間）も、４秒に限定されるものではなく、ゴルフクラブシャフト３０の特性、要求される振動の測定精度に応じて適宜変更可能である。   The sampling period for the acceleration signal by the FFT analyzer 20 is not limited to about 1 millisecond, and can be changed as appropriate according to the required measurement accuracy of vibration. Furthermore, the capture time (sampling time) for displacement by the FFT analyzer 20 is not limited to 4 seconds, and can be changed as appropriate according to the characteristics of the golf club shaft 30 and the required measurement accuracy of vibration. .

ＰＣ２２は、処理装置であって、一般的なパーソナルコンピュータと同様の構成を有するものであり、ＣＰＵ（図示せず）、メモリ（図示せず）を備え、キーボード、マウスなどのコンピュータの入力に用いられる入力部（図示せず）と、入力部からの入力情報およびＣＰＵで情報処理された情報を表示するＬＣＤなどのモニタ２４とを有する。   The PC 22 is a processing device and has a configuration similar to that of a general personal computer, and includes a CPU (not shown) and a memory (not shown), and is used for input of a computer such as a keyboard and a mouse. Input unit (not shown), and a monitor 24 such as an LCD for displaying input information from the input unit and information processed by the CPU.

ＰＣ２２は、ＦＦＴアナライザ２０から出力された第１の加速度センサ１４による上下方向Ｖの加速度の時系列データ（時間波形）、第２の加速度センサ１６による左右方向Ｈの加速度の時系列データ（時間波形）がメモリ（図示せず）に記憶される。
ＰＣ２２は、上下方向Ｖの加速度波形のデータ、左右方向Ｈの加速度波形のデータについて、ＦＦＴ（高速フーリエ変換）をする機能を有する。 The PC 22 outputs time series data (time waveform) of acceleration in the vertical direction V by the first acceleration sensor 14 output from the FFT analyzer 20 and time series data (time waveform) of acceleration in the left and right direction H by the second acceleration sensor 16. ) Is stored in a memory (not shown).
The PC 22 has a function of performing FFT (Fast Fourier Transform) on the acceleration waveform data in the vertical direction V and the acceleration waveform data in the horizontal direction H.

また、ＰＣ２２は、１次共振周波数、２次共振周波数および３次共振周波数等の共振周波数を算出する機能、例えば、各共振周波数を中心にした１／３オクターブバンドフィルタ等のバンドパスフィルタの設定する機能を有する。
また、ＰＣ２２は、１／３オクターブバンドパスフィルタを用いた共振周波数別の振幅特性データを算出する機能、各共振周波数の振幅特性データに対して逆ＦＦＴを行い、加速度波形のデータを作成する機能を有する。 The PC 22 has a function of calculating a resonance frequency such as a primary resonance frequency, a secondary resonance frequency, and a tertiary resonance frequency, for example, setting of a band pass filter such as a 1/3 octave band filter centering on each resonance frequency. It has the function to do.
Further, the PC 22 has a function of calculating amplitude characteristic data for each resonance frequency using a 1/3 octave bandpass filter, and a function of performing inverse FFT on the amplitude characteristic data of each resonance frequency to create acceleration waveform data. Have

また、ＰＣ２２は、波形のデータについて、山側および谷側のピークを抽出し、かつそのピークの時間を求める機能を有する。
さらには、ＰＣ２２は、各共振周波数における減衰率、および損失係数を求める機能を有する。 The PC 22 has a function of extracting peaks and valleys from the waveform data and obtaining the peak time.
Furthermore, the PC 22 has a function of obtaining an attenuation factor and a loss coefficient at each resonance frequency.

なお、ＰＣ２２においては、得られたデータの数値についての四則演算の機能を有することは言うまでもなく、更には上記機能を発揮する段階で作成された各種のデータについては、適宜メモリに記憶させるとともに、モニタ２４に、グラフ、数値等の形態で表示させる機能を有する。
また、ＰＣ２２における上記各種の機能は、ＦＦＴアナライザを利用した周波数解析をするときになされる各種の信号処理をする際に用いられる市販のソフトウェアによる実現できるものである。 In addition, it goes without saying that the PC 22 has a function of four arithmetic operations on the numerical value of the obtained data, and further, various data created at the stage of exhibiting the above functions are appropriately stored in a memory, The monitor 24 has a function of displaying in the form of a graph, a numerical value, or the like.
The various functions in the PC 22 can be realized by commercially available software used when performing various signal processing performed when performing frequency analysis using an FFT analyzer.

なお、ＰＣ２２は、ＦＦＴアナライザ２０を介することなく、第１の加速度センサ１４の加速度信号、および第２の加速度センサ１６の加速度信号をチャージアンプ１８から直接、取り込むこともできる。この場合、第１の加速度センサ１４の加速度信号、および第２の加速度センサ１６の加速度信号をＡＤ変換するＡＤ変換ボードをＰＣ２２に組み込み、チャージアンプ１８と接続し、各加速度信号をコンピュータ２０に直接取り込む構成としてもよい。   Note that the PC 22 can directly take in the acceleration signal of the first acceleration sensor 14 and the acceleration signal of the second acceleration sensor 16 from the charge amplifier 18 without going through the FFT analyzer 20. In this case, an AD conversion board for AD converting the acceleration signal of the first acceleration sensor 14 and the acceleration signal of the second acceleration sensor 16 is incorporated in the PC 22 and connected to the charge amplifier 18, and each acceleration signal is directly connected to the computer 20. It is good also as a structure to take in.

本発明の評価方法および評価システムにおいては、評価対象は、棒状体のものであれば、断面形状が円、矩形であってもよく、特に限定されるものではない。また、ゴルフクラブシャフト３０に限定されるものではなく、釣り竿、アイスホッケーのスティック、グランドホッケーのスティック等についても、振動減衰等の振動特性を評価することができる。   In the evaluation method and the evaluation system of the present invention, as long as the evaluation target is a rod-like body, the cross-sectional shape may be a circle or a rectangle, and is not particularly limited. Further, the present invention is not limited to the golf club shaft 30, and vibration characteristics such as vibration damping can be evaluated for fishing rods, ice hockey sticks, ground hockey sticks, and the like.

次に、本実施形態のゴルフクラブシャフトの評価方法について、ゴルフクラブシャフト３０として、長さが１１７０ｍｍのカーボンシャフトを例にして説明する。
図２は、本発明の第１の実施形態の棒状体の評価方法によるゴルフクラブシャフトの評価方法を工程順に示すフローチャートである。 Next, the golf club shaft evaluation method of the present embodiment will be described by taking a carbon shaft having a length of 1170 mm as an example as the golf club shaft 30.
FIG. 2 is a flowchart showing the golf club shaft evaluation method according to the rod-like body evaluation method of the first embodiment of the present invention in the order of steps.

先ず、図１に示す評価システム１０の固定部１２にゴルフクラブシャフト３０の後端部３２ａを、例えば、１７０ｍｍの固定長さで固定する。この固定部１２によるゴルフクラブシャフト３０の固定圧力は、例えば、１９．６Ｎ（２ｋｇｆ）である。   First, the rear end portion 32a of the golf club shaft 30 is fixed to the fixing portion 12 of the evaluation system 10 shown in FIG. 1 with a fixing length of 170 mm, for example. The fixing pressure of the golf club shaft 30 by the fixing portion 12 is, for example, 19.6 N (2 kgf).

次に、第１の加速度センサ１４を、距離Ｌ１が１００ｍｍの第１の取付位置αに取り付ける。第２の加速度センサ１６を、第２の取付位置βに取り付ける。
第１の加速度センサ１４および第２の加速度センサ１６は、チャージアンプ１８に接続される。 Next, the first acceleration sensor 14, the distance L 1 is attached to the α first mounting position of 100 mm. The second acceleration sensor 16 is attached to the second attachment position β.
The first acceleration sensor 14 and the second acceleration sensor 16 are connected to a charge amplifier 18.

次に、ゴルフクラブシャフト３０の先端３０ａを下方向ＶＬに２０ｍｍ変位させた後、先端３０ａを解放する。これにより、ゴルフクラブシャフト３０の先端部３２ｂが振動する（ステップＳ１０）。
このとき、第１の加速度センサ１４による左右方向Ｈの加速度の計測と、第２の加速度センサ１６による上下方向Ｖの加速度の計測とが同時にされて、第１の加速度センサ１４による左右方向Ｈの加速度信号がチャージアンプ１８に入力されて変換された後、加速度信号がＦＦＴアナライザ２０に取り込まれる。また、第２の加速度センサ１６による上下方向Ｖの加速度信号もチャージアンプ１８に入力された後、チャージアンプ１８の出力信号がＦＦＴアナライザ２０に取り込まれる。 Next, after the tip 30a of the golf club shaft 30 is displaced 20 mm downward VL , the tip 30a is released. Thereby, the front-end | tip part 32b of the golf club shaft 30 vibrates (step S10).
At this time, the measurement of the acceleration in the horizontal direction H by the first acceleration sensor 14 and the measurement of the acceleration in the vertical direction V by the second acceleration sensor 16 are performed at the same time, and the measurement in the horizontal direction H by the first acceleration sensor 14 is performed. After the acceleration signal is input to the charge amplifier 18 and converted, the acceleration signal is taken into the FFT analyzer 20. Further, the acceleration signal in the vertical direction V from the second acceleration sensor 16 is also input to the charge amplifier 18, and then the output signal of the charge amplifier 18 is taken into the FFT analyzer 20.

本実施形態においては、ＦＦＴアナライザ２０には、先端３０ａを解放した後、例えば、約１ミリ秒のサンプリング周期（４秒／４０９６間隔）で４秒間（サンプリング時間）の加速度信号が取得される。これにより、図３（ａ）に示すように、上下方向Ｖの加速度の時系列データ４０ａ、および図３（ｂ）に示すように、左右方向Ｈの加速度の時系列データ４０ｂが得られる（ステップＳ１２）。
図３（ａ）に示すように上下方向Ｖの加速度は、１．４秒位まで減少していき、その後、変位がまた大きくなっている。図３（ｂ）に示す左右方向Ｈの加速度は、０．５秒位から１．５秒位まで減少している。
上下方向Ｖの加速度の時系列データ４０ａ、および左右方向Ｈの加速度の時系列データ４０ｂがＦＦＴアナライザ２０からＰＣ２２に出力されて、ＰＣ２２のメモリに記憶される。 In the present embodiment, the FFT analyzer 20 acquires an acceleration signal of 4 seconds (sampling time) at a sampling period of about 1 millisecond (4 seconds / 4096 interval) after releasing the tip 30a. Thereby, as shown in FIG. 3A, time series data 40a of acceleration in the vertical direction V and time series data 40b of acceleration in the horizontal direction H are obtained as shown in FIG. S12).
As shown in FIG. 3A, the acceleration in the vertical direction V decreases to about 1.4 seconds, and then the displacement increases again. The acceleration in the left-right direction H shown in FIG. 3B decreases from about 0.5 seconds to about 1.5 seconds.
The time series data 40a of the acceleration in the vertical direction V and the time series data 40b of the acceleration in the horizontal direction H are output from the FFT analyzer 20 to the PC 22 and stored in the memory of the PC 22.

次に、ＰＣ２２において、上下方向Ｖの加速度の時系列データ４０ａ、および左右方向Ｈの加速度の時系列データ４０ｂについて、ＦＦＴ（高速フーリエ変換）を行う（ステップＳ１４）。これにより、上下方向Ｖの加速度の時系列データ４０ａから、例えば、図４に示すような振幅特性データ４２が得られる。図４の縦軸のパワースペクトル（ｄＢ）は、振動の振幅の大きさを示すものである。
なお、左右方向Ｈの加速度の時系列データ４０ｂについてもＦＦＴを行うため、図示はしないが、時系列データ４０ａと同様に図４に示すような振幅特性データ４２が得られる。 Next, in the PC 22, FFT (Fast Fourier Transform) is performed on the time series data 40a of acceleration in the vertical direction V and the time series data 40b of acceleration in the horizontal direction H (step S14). Thereby, for example, amplitude characteristic data 42 as shown in FIG. 4 is obtained from the time series data 40a of the acceleration in the vertical direction V. The power spectrum (dB) on the vertical axis in FIG. 4 indicates the magnitude of the vibration amplitude.
Since the FFT is also performed on the time series data 40b of acceleration in the horizontal direction H, although not shown, amplitude characteristic data 42 as shown in FIG. 4 is obtained similarly to the time series data 40a.

ここで、ＦＦＴをして得られるフーリエ変換データは、振動の周波数に対する振動の振幅レベルを示す振幅特性データと、振動の周波数に対する振動の位相を示す位相特性データとを含むものである。
高速フーリエ変換の処理自体は、公知の方法であるため、その詳細な説明を省略する。この高速フーリエ変換は、例えば、振動の時間波形データである変位の時系列データからフーリエ係数を求め、このフーリエ係数を演算することで周波数毎に振幅と位相を求めるものである。周波数と振幅との関係を表すものが、振動の周波数に対する振動の振幅レベルを示す振幅特性データとなる。 Here, the Fourier transform data obtained by performing FFT includes amplitude characteristic data indicating the vibration amplitude level with respect to the vibration frequency and phase characteristic data indicating the vibration phase with respect to the vibration frequency.
Since the fast Fourier transform process itself is a known method, its detailed description is omitted. In this fast Fourier transform, for example, a Fourier coefficient is obtained from time series data of displacement, which is vibration time waveform data, and the amplitude and phase are obtained for each frequency by calculating the Fourier coefficient. What represents the relationship between the frequency and the amplitude is amplitude characteristic data indicating the amplitude level of the vibration with respect to the frequency of the vibration.

次に、上述の振幅特性データに基づいて、ＰＣ２２は、ゴルフクラブシャフト３０の上下方向Ｖおよび左右方向Ｈにおける１次共振周波数、２次共振周波数、および３次共振周波数を含む複数の共振周波数を求める（ステップＳ１６）。
例えば、ＰＣ２２において、図４に示す振幅特性データ４２について、周波数解析を行い、振幅特性データ４２に含まれる複数のピーク４２ａ、４２ｂ、４２ｃの値に対応する周波数の値を、それぞれ求めることにより、１次共振周波数ｆ１、２次共振周波数ｆ２、３次共振数波数ｆ３を得ることができる。
また、図４では、振幅特性データ４２を１つだけ示しているが、本実施形態においては、左右方向Ｈの振幅特性データについても、上下方向Ｖと同様に、１次共振周波数〜３次共振周波数を求め、合計６個の共振周波数を求める。 Next, based on the amplitude characteristic data described above, the PC 22 has a plurality of resonance frequencies including a primary resonance frequency, a secondary resonance frequency, and a tertiary resonance frequency in the vertical direction V and the horizontal direction H of the golf club shaft 30. Obtained (step S16).
For example, the PC 22 performs frequency analysis on the amplitude characteristic data 42 shown in FIG. 4 and obtains the frequency values corresponding to the values of the plurality of peaks 42a, 42b, and 42c included in the amplitude characteristic data 42, respectively. The primary resonance frequency f 1 , the secondary resonance frequency f 2 , and the tertiary resonance frequency wave number f 3 can be obtained.
In FIG. 4, only one amplitude characteristic data 42 is shown, but in this embodiment, the amplitude characteristic data in the horizontal direction H is also similar to the vertical direction V in the primary resonance frequency to the tertiary resonance. The frequency is obtained and a total of 6 resonance frequencies are obtained.

また、本実施形態の評価方法においては、１次共振周波数が少なくとも求めることができればよいため、振幅特性データ４２から、２次共振周波数および３次共振周波数を、必ずしも求める必要はない。
本実施形態においては、振動特性に与える影響が大きい１次共振周波数〜３次共振周波数ｆ１、ｆ２、ｆ３について求めることが好ましい。なお、何次の共振周波数まで求めるかについては、特に限定されるものではない。 Further, in the evaluation method of the present embodiment, it is only necessary to obtain at least the primary resonance frequency, and therefore it is not always necessary to obtain the secondary resonance frequency and the tertiary resonance frequency from the amplitude characteristic data 42.
In the present embodiment, it is preferable to obtain the primary resonance frequency to the tertiary resonance frequencies f 1 , f 2 , and f 3 that have a great influence on the vibration characteristics. In addition, it is not specifically limited about what resonance frequency is calculated | required.

次に、ＰＣ２２は、ステップＳ１８で求めた合計６個の１次共振周波数〜３次共振周波数ｆ１、ｆ２、ｆ３のそれぞれに対応して、各共振周波数を中心周波数とするバンドパスフィルタを設定する（ステップＳ１８）。
ここで、図５は、縦軸にフィルタ値をとり、横軸に周波数をとって、本実施形態で用いたバンドパスフィルタを説明するためのグラフである。
図５の縦軸のフィルタ値とは、ゲインに相当するものであり、例えば、フィルタ値が０ならば入力信号が１００％阻止され、フィルタ値が０．５ならば入力信号が５０％通過され、フィルタ値が１ならば入力信号が１００％通過される。 Next, the PC 22 corresponds to each of the total of six primary resonance frequencies to tertiary resonance frequencies f 1 , f 2 , and f 3 obtained in step S 18, and a band pass filter having each resonance frequency as a center frequency. Is set (step S18).
Here, FIG. 5 is a graph for explaining the bandpass filter used in this embodiment, with the filter value on the vertical axis and the frequency on the horizontal axis.
The filter value on the vertical axis in FIG. 5 corresponds to the gain. For example, if the filter value is 0, the input signal is blocked by 100%, and if the filter value is 0.5, the input signal is passed by 50%. If the filter value is 1, 100% of the input signal is passed.

図５に示すように、バンドパスフィルタの特性は、中心周波数ｆｃと、下限遮断周波数ｆＬと、上限遮断周波数ｆＨとによって規定される。
したがって、ステップＳ１８で求めた各共振周波数ｆ１、ｆ２、ｆ３に対応して、図５に示すように、中心周波数ｆｃ＝ｆ１のバンドパスフィルタ４４ａと、中心周波数ｆｃ＝ｆ２のバンドパスフィルタ４４ｂと、中心周波数ｆｃ＝ｆ３のバンドパスフィルタとを設定する。なお、図５において、３次共振周波数ｆ３のバンドパスフィルタの中心周波数ｆｃと、下限遮断周波数ｆＬと、上限遮断周波数ｆＨとが周波数の表示範囲を超えているため、その図示を省略している。 As shown in FIG. 5, the characteristics of the bandpass filter is defined and the center frequency fc, and lower cut-off frequency f L, by the upper cut-off frequency f H.
Accordingly, as shown in FIG. 5, the band-pass filter 44a having the center frequency fc = f 1 and the center frequency fc = f 2 corresponding to the resonance frequencies f 1 , f 2 , and f 3 obtained in step S18. to set a bandpass filter 44b, a bandpass filter having a center frequency fc = f 3. Incidentally, omitted in FIG. 5, the center frequency fc of the band-pass filter of the third-order resonance frequency f 3, and the lower cutoff frequency f L, because the the upper cut-off frequency f H is beyond the display range of frequencies, the illustrated doing.

本実施形態において、バンドパスフィルタとして、例えば、１／３オクターブバンドパスフィルタを用いる。この１／３オクターブバンドパスフィルタにおける下限遮断周波数ｆＬ、上限遮断周波数ｆＨについては、下記（１）式と（２）式により求めることができる。 In the present embodiment, for example, a 1/3 octave bandpass filter is used as the bandpass filter. The lower limit cutoff frequency f L and the upper limit cutoff frequency f H in the 1/3 octave bandpass filter can be obtained by the following formulas (1) and (2).

ｆＬ＝ｆｃ・２（−１／６） （１）
ｆＨ＝ｆｃ・２（１／６） （２） f L = fc · 2 (−1/6) (1)
f H = fc · 2 (1/6) (2)

本実施形態における１／３オクターブバンドパスフィルタは、ＪＩＳ Ｃ１５１３に規定されている１／３オクターブの規格の考え方と同様である。ただし、上記ＪＩＳ規格では、中心周波数ｆｃが予め決定されているのに対して、本実施形態では、中心周波数ｆｃがステップＳ１８で求められた各共振周波数である点が異なる。   The 1/3 octave bandpass filter in this embodiment is the same as the concept of the 1/3 octave standard defined in JIS C1513. However, the center frequency fc is determined in advance in the JIS standard, whereas the present embodiment is different in that the center frequency fc is each resonance frequency obtained in step S18.

なお、本実施形態においては、バンドパスフィルタに１／３オクターブバンドパスフィルタを用いたが、これに限定されるものではない。バンドパスフィルタとしては、例えば、図４に示す振幅特性データ４２から、目的とする共振周波数成分の振動波形を正確に分離することができればよく、バンドパスフィルタの下限遮断周波数ｆＬ、上限遮断周波数ｆＨをどのように設定するかは、特に限定されるものではない。 In the present embodiment, a 1/3 octave bandpass filter is used as the bandpass filter, but the present invention is not limited to this. As the band-pass filter, for example, it is only necessary that the vibration waveform of the target resonance frequency component can be accurately separated from the amplitude characteristic data 42 shown in FIG. 4, and the lower-limit cutoff frequency f L and the upper-limit cutoff frequency of the band-pass filter are used. How f H is set is not particularly limited.

次に、ＰＣ２２は、図４に示す振幅特性データ４２を、ステップＳ２０で設定した３つのバンドパスフィルタにそれぞれかける（ステップＳ２０）。これにより、各バンドパスフィルタから３つの共振周波数（１次共振周波数ｆ１、２次共振周波数ｆ２、３次共振周波数ｆ３）のそれぞれに対応した３つの共振周波数別振幅特性データが得られる。
なお、左右方向Ｈの振幅特性データ（図示せず）についても、３つの共振周波数（１次共振周波数ｆ１、２次共振周波数ｆ２、３次共振周波数ｆ３）のそれぞれに対応した３つの共振周波数別振幅特性データが得られる。 Next, the PC 22 applies the amplitude characteristic data 42 shown in FIG. 4 to each of the three bandpass filters set in step S20 (step S20). Thus, three resonance frequency-specific amplitude characteristic data corresponding to the three resonance frequencies (primary resonance frequency f 1 , secondary resonance frequency f 2 , and tertiary resonance frequency f 3 ) are obtained from each bandpass filter. .
Note that the amplitude characteristic data (not shown) in the left-right direction H also includes three resonance frequencies (primary resonance frequency f 1 , secondary resonance frequency f 2 , tertiary resonance frequency f 3 ) corresponding to each of the three resonance frequencies. Amplitude characteristic data by resonance frequency is obtained.

次に、ＰＣ２２は、ステップＳ２０で得られた３つの共振周波数別振幅特性データのそれぞれと、ステップＳ１４で求めた位相特性データとに基づいて逆ＦＦＴ（高速フーリエ変換）を行う（ステップＳ２２）。
これにより、３つの共振周波数（１次共振周波数ｆ１、２次共振周波数ｆ２、３次共振周波数ｆ３）のそれぞれに対応した共振周波数別の時間波形データを得る（ステップＳ２４）。この場合、上下方向および左右方向について、それぞれ、図３（ａ）に示すような時系列データおよび図３（ｂ）に示すような時系列データが得られる。このようにして、ステップＳ２４においては、各共振周波数別に分離された時間波形データが得られる。 Next, the PC 22 performs inverse FFT (fast Fourier transform) based on each of the three resonance frequency-specific amplitude characteristic data obtained in step S20 and the phase characteristic data obtained in step S14 (step S22).
Thereby, time waveform data for each resonance frequency corresponding to each of the three resonance frequencies (primary resonance frequency f 1 , secondary resonance frequency f 2 , tertiary resonance frequency f 3 ) is obtained (step S 24). In this case, time-series data as shown in FIG. 3A and time-series data as shown in FIG. 3B are obtained in the vertical direction and the horizontal direction, respectively. In this way, time waveform data separated for each resonance frequency is obtained in step S24.

次に、ステップＳ２４で得られた各共振周波数毎の時間波形データについて、０クロス点を利用して、０点を境にプラス側のピークについてその時間とピークレベルを求め、更にはマイナス側のピークについてその時間とピークレベルを求める（ステップＳ２６）。   Next, with respect to the time waveform data for each resonance frequency obtained in step S24, the time and peak level are obtained for the plus side peak at the 0 point using the 0 cross point, and further, the minus side is obtained. The peak time and peak level are obtained (step S26).

各共振周波数において、上下方向Ｖにおける時間波形データのピーク値の絶対値をＵｐとし、左右方向Ｈにおける時間波形データのピーク値の絶対値をＬｐとするとき、時間波形データのピーク値の絶対値Ｕｐを基準として、このピーク値の絶対値Ｕｐに、最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを抽出する。
次に、基準となるピーク値の絶対値Ｕｐと、抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求める（ステップＳ２８）。なお、自乗平方根の値Ｖｐとは、以下の式により表されるものである。Ｖｐ＝（Ｕｐ２＋Ｌｐ２）１／２ At each resonance frequency, when the absolute value of the peak value of the time waveform data in the vertical direction V is Up and the absolute value of the peak value of the time waveform data in the horizontal direction H is Lp, the absolute value of the peak value of the time waveform data With reference to Up, the absolute value Lp of the peak value in the left-right direction at the nearest peak time is extracted from the absolute value Up of the peak value.
Next, a square root value Vp between the absolute value Up of the peak value serving as a reference and the absolute value Lp of the extracted peak value is obtained (step S28). The square root value Vp is represented by the following equation. Vp = (Up 2 + Lp 2 ) 1/2

次に、１次共振周波数、２次共振周波数、３次共振周波数毎に、減衰率を求める（ステップＳ３０）。本実施形態において、減衰率に、例えば、対数減衰率を用いる。
この場合、図６に示すように、例えば、縦軸に振動ピークレベルをとり、横軸に時間をとって、１次共振周波数について得られたＶｐの値を二乗して常用対数に変換した後、１０倍した数値をプロットする。これをプロット群５０ａとして図６に示す。このプロット群５０ａの各数値は、１０ｌｏｇＶｐ２により求められる。
そして、このプロット群５０ａについて、例えば、最小自乗法を用いて直線近似し、近似直線５２ａを求める。この近似直線５２ａの傾きが減衰率Ｄ（ｄＢ／秒）である。 Next, an attenuation factor is obtained for each primary resonance frequency, secondary resonance frequency, and tertiary resonance frequency (step S30). In the present embodiment, for example, a logarithmic attenuation rate is used as the attenuation rate.
In this case, as shown in FIG. 6, for example, after taking the vibration peak level on the vertical axis and taking the time on the horizontal axis, the value of Vp obtained for the primary resonance frequency is squared and converted to the common logarithm. Plot the number multiplied by 10. This is shown in FIG. 6 as a plot group 50a. Each number of the group of plots 50a is obtained by 10logVp 2.
The plot group 50a is approximated by a straight line using, for example, the least square method to obtain an approximate straight line 52a. The inclination of the approximate straight line 52a is the attenuation rate D (dB / second).

図６には、上下方向のピーク値Ｕｐを二乗して常用対数に変換した後、１０倍した数値をプロットしたプロット群５４ａを示している。
このプロット群５４ａの数値は、１０ｌｏｇＵｐ２により求められる。
さらに、図６には、左右方向のピーク値Ｌｐも二乗して常用対数に変換した後、１０倍した数値をプロットしたプロット群５６ａを示している。
このプロット群５６ａの数値は、１０ｌｏｇＬｐ２により求められる。
図６に示すように、自乗平方根の値Ｖｐに比して、上下方向のピーク値Ｕｐおよび左右方向のピーク値Ｌｐは直線性が劣っている。 FIG. 6 shows a plot group 54 a in which the peak value Up in the vertical direction is squared and converted to a common logarithm, and then a value multiplied by 10 is plotted.
The numbers in this plot group 54a is obtained by 10logUp 2.
Furthermore, FIG. 6 shows a plot group 56a in which the peak value Lp in the left-right direction is also squared and converted to a common logarithm, and then a value multiplied by 10 is plotted.
The numbers in this plot group 56a is obtained by 10logLp 2.
As shown in FIG. 6, the peak value Up in the vertical direction and the peak value Lp in the horizontal direction are inferior in linearity compared to the square root value Vp.

また、図７に示すように、例えば、縦軸に振動ピークレベルをとり、横軸に時間をとって、２次共振周波数について得られたＶｐの値を二乗して常用対数に変換した後、１０倍した数値をプロットする。これをプロット群５０ｂとして図７に示す。
そして、このプロット群５０ｂについて、例えば、最小自乗法を用いて直線近似し、近似直線５２ｂを求める。この近似直線５２ｂの傾きが減衰率Ｄ（ｄＢ／秒）である。
図７には、上下方向のピーク値Ｕｐを二乗して常用対数に変換した後、１０倍した数値をプロットしたプロット群５４ｂを示すとともに、左右方向のピーク値Ｌｐも二乗して常用対数に変換した後、１０倍した数値をプロットしたプロット群５６ｂを示している。
２次共振周波数においても、図７に示すように、自乗平方根の値Ｖｐに比して、上下方向のピーク値Ｕｐおよび左右方向のピーク値Ｌｐは直線性が劣っている。 Further, as shown in FIG. 7, for example, after taking the vibration peak level on the vertical axis and taking the time on the horizontal axis, the value of Vp obtained for the secondary resonance frequency is squared and converted to the common logarithm, Plot the number multiplied by 10. This is shown in FIG. 7 as a plot group 50b.
The plot group 50b is approximated by a straight line using, for example, the least square method to obtain an approximate straight line 52b. The slope of the approximate straight line 52b is the attenuation rate D (dB / second).
FIG. 7 shows a plot group 54b in which the peak value Up in the vertical direction is squared and converted to a common logarithm, and then a value multiplied by 10 is plotted. After that, the plot group 56b in which the numerical values multiplied by 10 are plotted is shown.
Also at the secondary resonance frequency, as shown in FIG. 7, the peak value Up in the vertical direction and the peak value Lp in the horizontal direction are inferior in linearity compared to the square root value Vp.

また、図８に示すように、例えば、縦軸に振動ピークレベルをとり、横軸に時間をとって、３次共振周波数について得られたＶｐの値を二乗して常用対数に変換した後、１０倍した数値をプロットする。これをプロット群５０ｃとして図８に示す。
そして、このプロット群５０ｃについて、例えば、最小自乗法を用いて直線近似し、近似直線５２ｃを求める。この近似直線の傾き５２ｃが減衰率Ｄ（ｄＢ／秒）である。
図８には、上下方向のピーク値Ｕｐを二乗して常用対数に変換した後、１０倍した数値をプロットしたプロット群５４ｃを示すとともに、左右方向のピーク値Ｌｐも二乗して常用対数に変換した後、１０倍した数値をプロットしたプロット群５６ｃを示している。
３次共振周波数においても、図８に示すように、自乗平方根の値Ｖｐに比して、上下方向のピーク値Ｕｐおよび左右方向のピーク値Ｌｐは直線性が劣っている。 Also, as shown in FIG. 8, for example, after taking the vibration peak level on the vertical axis and taking the time on the horizontal axis, the value of Vp obtained for the third-order resonance frequency is squared and converted to a common logarithm, Plot the number multiplied by 10. This is shown in FIG. 8 as a plot group 50c.
The plot group 50c is approximated by a straight line using, for example, the least square method to obtain an approximate straight line 52c. The slope 52c of this approximate line is the attenuation rate D (dB / second).
FIG. 8 shows a plot group 54c in which the peak value Up in the vertical direction is squared and converted to the common logarithm, and then the value multiplied by 10 is plotted, and the peak value Lp in the horizontal direction is also squared and converted to the common logarithm. After that, a plot group 56c is shown in which the values multiplied by 10 are plotted.
Also at the third-order resonance frequency, as shown in FIG. 8, the peak value Up in the vertical direction and the peak value Lp in the horizontal direction are inferior in linearity as compared to the square root value Vp.

ＰＣ２２においては、上記図６に示す近似直線５２ａ、図７に示す近似直線５２ｂ、図８に示す近似直線５２ｃを求める。各近似直線５２ａ、近似直線５２ｂおよび近似直線５２ｃの傾きを求め、減衰率Ｄを得る（ステップＳ３０）。
そして、各共振周波数における減衰率Ｄから、各共振周波数について、損失係数ηを求める（ステップＳ３２）。
なお、減衰率Ｄと損失係数ηとは式（３）で示される関係にある。このため、減衰率Ｄを求めることにより、損失係数ηを求めることができる。
η＝Ｄ／（２７．３・ｆｘ） （３） （ｆｘは各共振周波数） In the PC 22, the approximate straight line 52a shown in FIG. 6, the approximate straight line 52b shown in FIG. 7, and the approximate straight line 52c shown in FIG. 8 are obtained. The slopes of the approximate straight lines 52a, approximate straight lines 52b, and approximate straight lines 52c are obtained, and the attenuation rate D is obtained (step S30).
Then, the loss coefficient η is obtained for each resonance frequency from the attenuation rate D at each resonance frequency (step S32).
The attenuation rate D and the loss coefficient η are in a relationship represented by the equation (3). Therefore, the loss factor η can be obtained by obtaining the attenuation rate D.
η = D / (27.3 · fx) (3) (fx is each resonance frequency)

共振周波数ｆｘには、例えば、ステップＳ１６で求めた上下方向Ｖ（加振方向）の各共振周波数が用いられる。
なお、共振周波数別振動波形の波数をＮとし、この波数Ｎの発生時間をＴとしたときに、共振周波数ｆｎは、Ｎ／Ｔで求めることができる。このようにして得られた共振周波数の値を共振周波数ｆｘとして用いてもよい。 As the resonance frequency fx, for example, each resonance frequency in the vertical direction V (excitation direction) obtained in step S16 is used.
The resonance frequency fn can be obtained by N / T, where N is the wave number of the vibration waveform for each resonance frequency and T is the generation time of the wave number N. The value of the resonance frequency obtained in this way may be used as the resonance frequency fx.

本実施形態において、ブレーディング法により形成されたゴルフクラブシャフトについて取付角度を０°、４５°、９０°、１３５°と変えて、それぞれ１次共振周波数〜３次共振周波数における損失係数を求めた。さらには、取付角度が０°、４５°、９０°、１３５°のうち、左右方向の振動が最も小さいものにおける上下方向の振動波形データだけを用いて損失係数を求めた。この結果を図９に示す。
なお、取付角度０°とは、計測の最初に取り付けた状態のことである。また、取付角度４５°、９０°、１３５°とは、取付角度０°を基準として、時計周りにゴルフクラブシャフトを回転させた角度のことである。 In this embodiment, the golf club shaft formed by the braiding method was changed in the mounting angle to 0 °, 45 °, 90 °, and 135 °, and the loss coefficients at the primary resonance frequency to the tertiary resonance frequency were obtained. . Further, the loss coefficient was obtained using only the vibration waveform data in the vertical direction when the vibration in the horizontal direction is smallest among the mounting angles of 0 °, 45 °, 90 °, and 135 °. The result is shown in FIG.
Note that the attachment angle of 0 ° is a state of attachment at the beginning of measurement. Further, the attachment angles of 45 °, 90 °, and 135 ° are angles obtained by rotating the golf club shaft clockwise with the attachment angle of 0 ° as a reference.

図９には、１次共振周波数における損失係数をプロットしたプロット群６０ａと、２次共振周波数における損失係数をプロットしたプロット群６０ｂと、３次共振周波数における損失係数をプロットしたプロット群６０ｃとを示している。
本実施形態の評価方法によれば、図９に示すように、取付角度によらず、損失係数の値に差がない。また、左右方向の振動が最も小さいものにおける上下方向の振動波形データだけを用いて得られた損失係数と、本実施形態の評価方法による損失係数との差もない。 FIG. 9 shows a plot group 60a in which the loss coefficient at the primary resonance frequency is plotted, a plot group 60b in which the loss coefficient at the secondary resonance frequency is plotted, and a plot group 60c in which the loss coefficient at the tertiary resonance frequency is plotted. Show.
According to the evaluation method of this embodiment, as shown in FIG. 9, there is no difference in the value of the loss coefficient regardless of the mounting angle. Further, there is no difference between the loss coefficient obtained by using only the vertical vibration waveform data with the smallest vibration in the left-right direction and the loss coefficient according to the evaluation method of the present embodiment.

本実施形態の評価システムおよび評価方法によれば、ゴルフクラブシャフト３０について、取付角度に依存することなく、損失係数を求めることができる。しかも、上下方向Ｖと左右方向Ｈについて、振動を同時に計測するだけなので容易である。   According to the evaluation system and the evaluation method of the present embodiment, the loss coefficient can be obtained for the golf club shaft 30 without depending on the mounting angle. Moreover, it is easy because vibrations are only measured simultaneously in the vertical direction V and the horizontal direction H.

また、本実施形態の評価方法は、例えば、ゴルフクラブシャフトが設計した振動減衰能を有するかの等の品質管理に利用することができる。このとき、製造されたゴルフクラブ全品、または抜き取り検査する場合には、製造されたゴルフクラブシャフトのうち、適用に引き抜かれたものについて、損失係数を計測する。そして、ゴルフクラブシャフトが不良品で、剥離等が生じていると、全体の剛性が低下し、これにより共振周波数が低下する。また、剥離した箇所により内部摩擦が増え、損失係数が高くなる。このことから、図９に示すプロット群６０ａを良品とした場合、不良品のゴルフクラブシャフトは、符号６２ａで示すような位置にプロットされることになる。このようにして、不良品と良品とを区別することができる。
なお、２次共振周波数についても、図９に示すプロット群６０ｂを良品とした場合、不良品のゴルフクラブシャフトは、符号６２ｂで示すような位置にプロットされることになり、３次共振周波数についても、図９に示すプロット群６０ｃを良品とした場合、不良品のゴルフクラブシャフトは、符号６２ｃで示すような位置にプロットされることになる。 Moreover, the evaluation method of this embodiment can be used for quality control, such as whether the golf club shaft has the vibration damping ability designed, for example. At this time, when all manufactured golf clubs or a sampling inspection is performed, a loss factor is measured for one of the manufactured golf club shafts that has been pulled out for application. When the golf club shaft is a defective product and peeling or the like occurs, the overall rigidity is lowered, thereby reducing the resonance frequency. Further, the internal friction increases due to the peeled portion, and the loss factor increases. From this, when the plot group 60a shown in FIG. 9 is a non-defective product, the defective golf club shaft is plotted at a position indicated by the reference numeral 62a. In this way, it is possible to distinguish defective products from non-defective products.
In addition, regarding the secondary resonance frequency, when the plot group 60b shown in FIG. 9 is a non-defective product, the defective golf club shaft is plotted at a position indicated by reference numeral 62b. However, if the plot group 60c shown in FIG. 9 is a non-defective product, the defective golf club shaft is plotted at a position indicated by reference numeral 62c.

次に、本発明の第２の実施形態について説明する。
図１０は、本発明の第２の実施形態の棒状体の評価方法に用いられる評価システムを示す模式的斜視図である。
図１１は、本発明の第２の実施形態の棒状体の評価方法によるゴルフクラブシャフトの評価方法を工程順に示すフローチャートである。
本実施形態においては、図１に示す第１の実施形態の評価システム１０と同一構成物には、同一符号を付して、その詳細な説明は省略する。 Next, a second embodiment of the present invention will be described.
FIG. 10 is a schematic perspective view showing an evaluation system used in the rod-like body evaluation method of the second embodiment of the present invention.
FIG. 11 is a flowchart showing the golf club shaft evaluation method according to the rod-like body evaluation method of the second embodiment of the present invention in the order of steps.
In the present embodiment, the same components as those in the evaluation system 10 of the first embodiment shown in FIG. 1 are denoted by the same reference numerals, and detailed description thereof is omitted.

本実施形態の評価システム１１は、第１の実施形態の評価システム１０（図１参照）に比して、第１の加速度センサ１４および第２の加速度センサ１６にかえて第１のレーザ変位計１４および第２のレーザ変位計１６を用いている点、およびチャージアンプ１８が設けられていない点が異なり、それ以外の構成は、第１の実施形態と同様の構成であるため、その詳細な説明は省略する。   As compared with the evaluation system 10 (see FIG. 1) of the first embodiment, the evaluation system 11 of the present embodiment replaces the first acceleration sensor 14 and the second acceleration sensor 16 with a first laser displacement meter. 14 and the second laser displacement meter 16 and the point that the charge amplifier 18 is not provided, and the other configurations are the same as those in the first embodiment, and the details thereof are the same. Description is omitted.

本実施形態の評価システム１１において、第１のレーザ変位計７０および第２のレーザ変位計７４はＦＦＴアナライザ２０に接続されている。第１のレーザ変位計７０、第２のレーザ変位計７４およびＦＦＴアナライザ２０は、いずれもＰＣ２２により制御される。   In the evaluation system 11 of this embodiment, the first laser displacement meter 70 and the second laser displacement meter 74 are connected to the FFT analyzer 20. The first laser displacement meter 70, the second laser displacement meter 74, and the FFT analyzer 20 are all controlled by the PC 22.

第１のレーザ変位計７０は、左右方向Ｈにおけるゴルフクラブシャフト３０の先端部３２ｂの振動の変位を測定するものである。
第１のレーザ変位計７０は、例えば、透過式のレーザ変位計であり、投光部７０ａと、受光部７０ｂとを有し、投光部７０ａと、受光部７０ｂとが上下方向Ｖに所定の間隔をあけて対向して配置されている。本実施形態においては、投光部７０ａと受光部７０ｂとの間に評価対象物であるゴルフクラブシャフト３０が配置される。 The first laser displacement meter 70 measures the vibration displacement of the front end portion 32 b of the golf club shaft 30 in the left-right direction H.
The first laser displacement meter 70 is, for example, a transmission type laser displacement meter, and includes a light projecting unit 70a and a light receiving unit 70b. The light projecting unit 70a and the light receiving unit 70b are predetermined in the vertical direction V. Are arranged to face each other with an interval of. In the present embodiment, the golf club shaft 30 that is an evaluation object is disposed between the light projecting unit 70a and the light receiving unit 70b.

第１のレーザ変位計７０においては、投光部７０ａから帯状のレーザ光７２が受光部７０ｂに向けて照射され、レーザ光７２の一部がゴルフクラブシャフト３０で遮られて受光部７０ｂで受光する。受光部７０ｂは、受光したレーザ光７２に基づいて、左右方向Ｈのゴルフクラブシャフト３０の周面３０ｂのエッジを検出する。この左右方向Ｈの周面３０ｂのエッジの位置の時間変化を検出することにより、振動の際の左右方向Ｈの変位が得られる。第１のレーザ変位計７０は、この検出された変位を、例えば、アナログ信号の形態で変位信号としてＦＦＴアナライザ２０に出力する。   In the first laser displacement meter 70, a band-shaped laser beam 72 is emitted from the light projecting unit 70a toward the light receiving unit 70b, and a part of the laser beam 72 is blocked by the golf club shaft 30 and received by the light receiving unit 70b. To do. The light receiving unit 70 b detects the edge of the peripheral surface 30 b of the golf club shaft 30 in the left-right direction H based on the received laser light 72. By detecting the temporal change in the position of the edge of the peripheral surface 30b in the left-right direction H, a displacement in the left-right direction H during vibration is obtained. The first laser displacement meter 70 outputs the detected displacement to the FFT analyzer 20 as a displacement signal, for example, in the form of an analog signal.

第２のレーザ変位計７４は、第１のレーザ変位計７０の投光部７０ａおよび受光部７０ｂに対して９０°回転させて配置されており、投光部７４ａと受光部７４ｂとが左右方向Ｈに所定の間隔をあけて対向して配置されている。本実施形態においては、投光部７４ａと受光部７４ｂとの間に評価対象物であるゴルフクラブシャフト３０が配置される。   The second laser displacement meter 74 is disposed by being rotated by 90 ° with respect to the light projecting unit 70a and the light receiving unit 70b of the first laser displacement meter 70, and the light projecting unit 74a and the light receiving unit 74b are arranged in the left-right direction. It is arranged to face H with a predetermined interval. In the present embodiment, the golf club shaft 30 that is an evaluation object is disposed between the light projecting unit 74a and the light receiving unit 74b.

なお、第２のレーザ変位計７４の投光部７４ａと受光部７４ｂとは、第１のレーザ変位計７０の投光部７０ａおよび受光部７０ｂと同様の構成であるため、その詳細な説明は省略する。第２のレーザ変位計７４においては、周面３０ｂの上下方向Ｖのエッジの位置の時間変化を検出することにより、振動の際の上下方向Ｖの変位が得られる。第２のレーザ変位計７４は、この検出された変位を、例えば、アナログ信号の形態で変位信号としてＦＦＴアナライザ２０に出力する。   The light projecting unit 74a and the light receiving unit 74b of the second laser displacement meter 74 have the same configuration as the light projecting unit 70a and the light receiving unit 70b of the first laser displacement meter 70, and the detailed description thereof will be described. Omitted. In the second laser displacement meter 74, the displacement in the vertical direction V at the time of vibration is obtained by detecting the temporal change in the position of the edge in the vertical direction V of the peripheral surface 30b. The second laser displacement meter 74 outputs the detected displacement to the FFT analyzer 20 as a displacement signal in the form of an analog signal, for example.

この第１のレーザ変位計７０は、第１の加速度センサ１４が取り付けられる第１の位置αにレーザ光７２を照射する。また、第２のレーザ変位計７４は、第２の加速度センサ１６が取り付けられる第２の位置βにレーザ光７６を照射する。
第１のレーザ変位計７０および第２のレーザ変位計７４は、ゴルフクラブシャフト３０の左右方向Ｈと上下方向Ｖの変位を同時に計測するものであり、計測タイミングなどはＰＣ２２により制御される。 The first laser displacement meter 70 irradiates a laser beam 72 to a first position α where the first acceleration sensor 14 is attached. The second laser displacement meter 74 irradiates the laser beam 76 to the second position β where the second acceleration sensor 16 is attached.
The first laser displacement meter 70 and the second laser displacement meter 74 measure the displacement in the horizontal direction H and the vertical direction V of the golf club shaft 30 at the same time, and the measurement timing and the like are controlled by the PC 22.

第１のレーザ変位計７０および第２のレーザ変位計７４は、いずれも、投光部７０ａ、７４ａと受光部７０ｂ、７４ｂとの間の距離を変えることができるとともに、投光部７０ａ、７４ａおよび受光部７０ｂ、７４ｂは、ゴルフクラブシャフト３０の長手方向Ｌにおける位置も適宜調整することができる。   Both the first laser displacement meter 70 and the second laser displacement meter 74 can change the distance between the light projecting units 70a and 74a and the light receiving units 70b and 74b, and the light projecting units 70a and 74a. Further, the positions of the light receiving portions 70b and 74b in the longitudinal direction L of the golf club shaft 30 can be adjusted as appropriate.

本実施形態において、ＦＦＴアナライザ２０は、例えば、約１ミリ秒のサンプリング周期で４秒間、第１のレーザ変位計７０からアナログ信号の形態で出力された左右方向Ｈの変位信号と、第２のレーザ変位計７４からアナログ信号の形態で出力された上下方向Ｖの変位信号とを取り込む。   In this embodiment, the FFT analyzer 20 includes, for example, a displacement signal in the left-right direction H output in the form of an analog signal from the first laser displacement meter 70 for 4 seconds at a sampling period of about 1 millisecond, The displacement signal in the vertical direction V outputted from the laser displacement meter 74 in the form of an analog signal is taken in.

ＦＦＴアナライザ２０により、第１のレーザ変位計７０によるゴルフクラブシャフト３０の左右方向Ｈの変位の時系列データ（時間波形）、すなわち、左右方向Ｈの振動の時系列データが得られる。さらには、ＦＦＴアナライザ２０により、第２のレーザ変位計７４によるゴルフクラブシャフト３０の上下方向Ｖの変位の時系列データ（時間波形）、すなわち、上下方向Ｖの振動の時系列データが得られる。
このようにして、ゴルフクラブシャフト３０を振動させたときの上下方向Ｖの変位の時系列データ（以下、振動波形データともいう）および左右方向Ｈの変位の時系列データ（振動波形データ）が得られる。
これらの上下方向Ｖの振動の時系列データおよび左右方向Ｈの振動の時系列データがＰＣ２２に出力される。 The FFT analyzer 20 obtains time series data (time waveform) of the horizontal displacement H of the golf club shaft 30 by the first laser displacement meter 70, that is, time series data of vibration in the horizontal direction H. Furthermore, the FFT analyzer 20 obtains time-series data (time waveform) of the vertical displacement V of the golf club shaft 30 by the second laser displacement meter 74, that is, time-series data of vibration in the vertical direction V.
Thus, time series data of displacement in the vertical direction V (hereinafter also referred to as vibration waveform data) and time series data of vibration in the horizontal direction H (vibration waveform data) when the golf club shaft 30 is vibrated are obtained. It is done.
The time series data of the vibration in the vertical direction V and the time series data of the vibration in the horizontal direction H are output to the PC 22.

なお、ＦＦＴアナライザ２０による変位についてのサンプリング周期は、約１ミリ秒に限定されるものではなく、要求される振動の測定精度に応じて適宜変更可能である。
さらには、ＦＦＴアナライザ２０による変位についての取り込み時間（サンプリング時間）も、４秒に限定されるものではなく、ゴルフクラブシャフト３０の特性、要求される振動の測定精度に応じて適宜変更可能である。 Note that the sampling period for displacement by the FFT analyzer 20 is not limited to about 1 millisecond, and can be appropriately changed according to the required measurement accuracy of vibration.
Furthermore, the capture time (sampling time) for displacement by the FFT analyzer 20 is not limited to 4 seconds, and can be changed as appropriate according to the characteristics of the golf club shaft 30 and the required measurement accuracy of vibration. .

本実施形態において、ＰＣ２２は、ＦＦＴアナライザ２０から出力された第１のレーザ変位計７０で得られた左右方向Ｈの変位の時系列データ（振動波形データ）、第２のレーザ変位計７４で得られた上下方向Ｖの変位の時系列データ（振動波形データ）がメモリ（図示せず）に記憶される。ＰＣ２２は、振動波形のデータについて、山側および谷側のピークを抽出し、かつそのピークの時間を求める機能を有する。   In this embodiment, the PC 22 obtains the time series data (vibration waveform data) of the displacement in the left-right direction H obtained by the first laser displacement meter 70 output from the FFT analyzer 20 and obtained by the second laser displacement meter 74. The time series data (vibration waveform data) of the displacement in the vertical direction V is stored in a memory (not shown). The PC 22 has a function of extracting peak and valley peaks from the vibration waveform data and obtaining the peak time.

なお、ＰＣ２２は、ＦＦＴアナライザ２０による第１のレーザ変位計７０の変位信号、および第２のレーザ変位計７４の変位信号を直接取り込むこともできる。この場合、第１のレーザ変位計７０の変位信号、および第２のレーザ変位計７４の変位信号をＡ／Ｄ変換するＡ／Ｄ変換ボードをＰＣ２２に組み込み、各変位信号をコンピュータ２０に直接取り込む。   The PC 22 can also directly capture the displacement signal of the first laser displacement meter 70 and the displacement signal of the second laser displacement meter 74 by the FFT analyzer 20. In this case, an A / D conversion board for A / D converting the displacement signal of the first laser displacement meter 70 and the displacement signal of the second laser displacement meter 74 is incorporated in the PC 22 and each displacement signal is directly captured in the computer 20. .

本実施形態において、第１のレーザ変位計７０および第２のレーザ変位計７４は、透過式のレーザ変位計としたが、これに限定されるものではない。例えば、第１のレーザ変位計７０および第２のレーザ変位計７４は、反射型のレーザ変位計であっても、レーザーフォーカス型のレーザ変位計でもよい。   In the present embodiment, the first laser displacement meter 70 and the second laser displacement meter 74 are transmissive laser displacement meters, but are not limited thereto. For example, the first laser displacement meter 70 and the second laser displacement meter 74 may be a reflection type laser displacement meter or a laser focus type laser displacement meter.

次に、本実施形態のゴルフクラブシャフトの評価方法について、図１１に示すフローチャートを参照して説明する。
本実施形態のように、第１のレーザ変位計７０および第２のレーザ変位計７４を用いた場合、得られる振動波形のデータは、１次共振周波数が支配的であり、高次の振動が小さい。このため、以下のようにして、上下方向の振動波形データと、左右方向の振動波形データとについて、ＦＦＴなどの信号処理を加えることなく、上下方向の振動波形データおよび左右方向の振動波形データについて、ピーク値の絶対値Ｕｐ、Ｌｐを求めて、上記Ｖｐを求めて、損失係数を求めることができる。 Next, the golf club shaft evaluation method of this embodiment will be described with reference to the flowchart shown in FIG.
When the first laser displacement meter 70 and the second laser displacement meter 74 are used as in the present embodiment, the data of the vibration waveform obtained is dominated by the primary resonance frequency, and higher-order vibrations are present. small. Therefore, with respect to the vertical vibration waveform data and the horizontal vibration waveform data, the vertical vibration waveform data and the horizontal vibration waveform data are not added to the vertical vibration waveform data and the horizontal vibration waveform data as follows. The absolute values Up and Lp of the peak value can be obtained, and the loss factor can be obtained by obtaining the Vp.

本実施形態のゴルフクラブシャフトの評価方法は、図２に示す第１の実施形態の評価方法と加振する工程（Ｓ４０）が第１の実施形態のステップＳ１０と同じ工程である。   The golf club shaft evaluation method of the present embodiment is the same as step S10 of the first embodiment in the step of vibrating (S40) with the evaluation method of the first embodiment shown in FIG.

ステップＳ４０において、先端３０ａを解放した後、例えば、約１ミリ秒のサンプリング周期（４秒／４０９６間隔）で４秒間（サンプリング時間）の変位信号が取得される。これにより、上下方向Ｖの変位の時系列データ、および左右方向Ｈの変位の時系列データが得られる（ステップＳ４２）。
上下方向Ｖの変位の時系列データ（図示せず）および左右方向Ｈの変位の時系列データ（図示せず）がＦＦＴアナライザ２０からＰＣ２２に出力されて、ＰＣ２２のメモリに記憶される。 In step S40, after the tip 30a is released, for example, a displacement signal of 4 seconds (sampling time) is acquired at a sampling period (4 seconds / 4096 interval) of about 1 millisecond. Thereby, time series data of displacement in the vertical direction V and time series data of displacement in the horizontal direction H are obtained (step S42).
Time series data (not shown) of displacement in the vertical direction V and time series data (not shown) of displacement in the horizontal direction H are output from the FFT analyzer 20 to the PC 22 and stored in the memory of the PC 22.

次に、ＰＣ２２において、上下方向Ｖの変位の時系列データおよび左右方向Ｈの変位の時系列データについて、０クロス点を利用して、０点を境にプラス側のピークについてその時間とピークレベルを求め、更にはマイナス側のピークについてその時間とピークレベルを求める（ステップＳ４４）。   Next, in the PC 22, with respect to time series data of displacement in the vertical direction V and time series data of displacement in the horizontal direction H, using the 0 cross point, the time and peak level of the plus side peak at the 0 point Further, the time and peak level of the minus peak are obtained (step S44).

次に、上下方向Ｖにおける時間波形データのピーク値の絶対値をＵｐとし、左右方向Ｈにおける時間波形データのピーク値の絶対値をＬｐとするとき、時間波形データのピーク値の絶対値Ｕｐを基準として、このピーク値の絶対値Ｕｐに、最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを抽出する。
次に、基準となるピーク値の絶対値Ｕｐと、抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求める（ステップＳ４６）。なお、自乗平方根の値Ｖｐは、第１の実施形態と同じ式で求める。 Next, when the absolute value of the peak value of the time waveform data in the vertical direction V is Up and the absolute value of the peak value of the time waveform data in the horizontal direction H is Lp, the absolute value Up of the peak value of the time waveform data is As a reference, the absolute value Lp of the peak value in the left-right direction at the nearest peak time is extracted as the absolute value Up of the peak value.
Next, a square root value Vp of the absolute value Up of the peak value serving as a reference and the absolute value Lp of the extracted peak value is obtained (step S46). Note that the square root value Vp is obtained by the same equation as in the first embodiment.

次に、減衰率を求める（ステップＳ４８）。本実施形態においては、減衰率に、例えば、対数減衰率を用いる。
この場合、図示はしないが、例えば、縦軸に振動ピークレベルをとり、横軸に時間をとって、得られたＶｐの値を二乗して常用対数に変換した後、１０倍した数値をプロットする。これらの各数値は、１０ｌｏｇＶｐ２により求められる。
そして、ＰＣ２２においては、例えば、最小自乗法を用いて直線近似し、近似直線を求める。この近似直線の傾きを求める。この傾きが減衰率Ｄ（ｄＢ／秒）である。
そして、共振周波数における減衰率Ｄから、損失係数ηを求める（ステップＳ５０）。この損失係数ηは、上記式（３）により求められる。
なお、共振周波数ｆｘは、上下方向Ｖの振動波形の波数をＮとし、この波数Ｎの発生時間をＴとしたときに、Ｎ／Ｔで求められるものである。 Next, an attenuation rate is obtained (step S48). In the present embodiment, for example, a logarithmic attenuation rate is used as the attenuation rate.
In this case, although not shown, for example, the vertical axis is the vibration peak level, the horizontal axis is time, the obtained Vp value is squared and converted to the common logarithm, and then the value multiplied by 10 is plotted. To do. Each of these figures are obtained by 10logVp 2.
In the PC 22, for example, a straight line approximation is performed using the least square method to obtain an approximate straight line. The inclination of this approximate straight line is obtained. This inclination is the attenuation rate D (dB / second).
Then, the loss coefficient η is obtained from the attenuation rate D at the resonance frequency (step S50). This loss coefficient η is obtained by the above equation (3).
The resonance frequency fx is obtained by N / T, where N is the wave number of the vibration waveform in the vertical direction V and T is the generation time of the wave number N.

このように、本実施形態の評価方法においては、第１の実施形態と同様に、振動減衰等の振動特性について、容易かつ高い精度で評価することができる。
しかも、本実施形態においては、レーザ変位計を用いた場合、非接触としているため、センサを取り付ける手間がなく、更に容易に評価することができる。
なお、上述のいずれの実施形態においても、減衰率に、対数減衰率を用いているが、特にこれに限定されるものではなく、また、減衰率を求める方法も、特に限定されるものではない。 As described above, in the evaluation method of the present embodiment, as in the first embodiment, vibration characteristics such as vibration attenuation can be evaluated easily and with high accuracy.
In addition, in the present embodiment, when a laser displacement meter is used, since it is non-contact, there is no need to attach a sensor, and evaluation can be performed more easily.
In any of the embodiments described above, the logarithmic attenuation factor is used as the attenuation factor. However, the present invention is not particularly limited to this, and the method for obtaining the attenuation factor is not particularly limited. .

本発明は、基本的に以上のようなものである。以上、本発明の棒状体の評価方法および棒状体の評価システムについて詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良または変更をしてもよいのはもちろんである。   The present invention is basically as described above. The rod-like body evaluation method and the rod-like body evaluation system of the present invention have been described in detail above. However, the present invention is not limited to the above-described embodiment, and various improvements or modifications can be made without departing from the spirit of the present invention. Of course.

１０ 評価システム
１２ 固定部
１４ 第１の加速度センサ
１６ 第２の加速度センサ
１８ ＦＦＴアナライザ
２０ コンピュータ（ＰＣ）
２２ モニタ
３０ ゴルフクラブシャフト
３０ａ 先端
３０ｂ 周面
３２ａ 後端部
３２ｂ先端部
４２ 振幅特性データ
４２ａ、４２ｂ、４２ｃ ピーク
７０ 第１のレーザ変位計
７４ 第２のレーザ変位計 DESCRIPTION OF SYMBOLS 10 Evaluation system 12 Fixed part 14 1st acceleration sensor 16 2nd acceleration sensor 18 FFT analyzer 20 Computer (PC)
22 monitor 30 golf club shaft 30a tip 30b peripheral surface 32a rear end 32b tip 42 amplitude characteristic data 42a, 42b, 42c peak 70 first laser displacement meter 74 second laser displacement meter

Claims (10)

棒状体の一端部を固定し、他端側を振動させ、前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測し、上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得する工程と、
前記計測された上下方向の振動の時系列データおよび前記左右方向の振動の時系列データについて、それぞれ、高速フーリエ変換し、振動の周波数に対する振動の振幅レベルを示す振幅特性データと、前記振動の周波数に対する振動の位相を示す位相特性データをそれぞれ作成する工程と、
前記上下方向の振動と前記左右方向の振動に対して、前記各振幅特性データから、それぞれ少なくとも１次共振周波数を求める工程と、
前記求められた各１次共振周波数を中心周波数としたバンドパスフィルタを設定する工程と、
前記各バンドパスフィルタを、それぞれ前記各振幅特性データにかけた後、それぞれ逆高速フーリエ変換し、１次共振周波数における時間波形データを、それぞれ求める工程と、
前記１次共振周波数の時間波形データについて、山側と谷側のピークのピーク時間およびピークレベルを求める工程と、
前記上下方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＵｐとし、前記左右方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出する工程と、
前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求める工程と、
前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する工程とを有することを特徴とする棒状体の評価方法。 One end of the rod-shaped body is fixed, the other end is vibrated, and the vertical vibration and the horizontal vibration of the other end are simultaneously measured for a predetermined time, and time series data of the vertical vibration and the horizontal direction are measured. Obtaining time series data of vibrations of
About the measured time series data of the vertical vibration and the time series data of the left and right vibration, respectively, fast Fourier transform, amplitude characteristic data indicating the amplitude level of the vibration with respect to the vibration frequency, and the frequency of the vibration Creating phase characteristic data indicating the phase of vibration with respect to
Obtaining at least a primary resonance frequency from each amplitude characteristic data for the vertical vibration and the horizontal vibration;
Setting a band-pass filter having each of the obtained primary resonance frequencies as a center frequency;
Applying each of the bandpass filters to each of the amplitude characteristic data, and then inverse fast Fourier transform to obtain time waveform data at the primary resonance frequency, respectively,
For the time waveform data of the primary resonance frequency, obtaining peak time and peak level of peaks on the peak side and the valley side;
The absolute value of the peak value in the time waveform data of the primary resonance frequency in the vertical direction is Up, the absolute value of the peak value in the time waveform data of the primary resonance frequency in the horizontal direction is Lp, and each peak value Extracting the absolute value Lp of the peak value in the left-right direction of the peak time closest to the absolute value Up of each peak value with the absolute value Up of
Obtaining a square root value Vp of the absolute value Up of the peak value and the absolute value Lp of the extracted peak value;
And a step of calculating an attenuation rate using the obtained value Vp and calculating a loss coefficient based on the attenuation rate. 棒状体の一端部を固定し、他端側を振動させ、前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測し、上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得する工程と、
前記上下方向の振動の時系列データおよび前記左右方向の振動の時系列データについて、それぞれ山側と谷側のピークのピーク時間およびピークレベルを求める工程と、
前記上下方向の振動の時系列データにおけるピーク値の絶対値をＵｐとし、前記左右方向の振動の時系列データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出する工程と、
前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求める工程と、
前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する工程とを有することを特徴とする棒状体の評価方法。 One end of the rod-shaped body is fixed, the other end is vibrated, and the vertical vibration and the horizontal vibration of the other end are simultaneously measured for a predetermined time, and time series data of the vertical vibration and the horizontal direction are measured. Obtaining time series data of vibrations of
For the time series data of the vibration in the vertical direction and the time series data of the vibration in the left and right direction, obtaining peak time and peak level of the peak on the mountain side and the valley side, respectively,
The absolute value of the peak value in the time series data of the vertical vibration is Up, the absolute value of the peak value in the time series data of the horizontal vibration is Lp, and the absolute value Up of each peak value is a reference. Extracting the absolute value Lp of the peak value in the left-right direction of the peak time closest to the absolute value Up of each peak value;
Obtaining a square root value Vp of the absolute value Up of the peak value and the absolute value Lp of the extracted peak value;
And a step of calculating an attenuation rate using the obtained value Vp and calculating a loss coefficient based on the attenuation rate. 前記各振幅特性データから求める共振周波数は、１次共振周波数、２次共振周波数および３次共振周波数である請求項１に記載の棒状体の評価方法。   The method for evaluating a rod-shaped body according to claim 1, wherein the resonance frequencies obtained from the amplitude characteristic data are a primary resonance frequency, a secondary resonance frequency, and a tertiary resonance frequency. 前記バンドパスフィルタは、１／３オクターブバンドパスフィルタである請求項１または３に記載の棒状体の評価方法。   The rod-shaped body evaluation method according to claim 1, wherein the band-pass filter is a 1/3 octave band-pass filter. 前記減衰率は、対数減衰率である請求項１〜４のいずれか１項に棒状体の評価方法。   The rod-like body evaluation method according to claim 1, wherein the attenuation rate is a logarithmic attenuation rate. 前記棒状体は、ゴルフクラブシャフトである請求項１〜５のいずれか１項に記載の棒状体の評価方法。   The rod-like body evaluation method according to claim 1, wherein the rod-like body is a golf club shaft. 棒状体の一端部を固定し、他端側を振動させたときの前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測する計測装置と、
前記計測装置で計測された上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得するデータ取得装置と、
前記計測された上下方向の振動の時系列データおよび前記左右方向の振動の時系列データについて、それぞれ、高速フーリエ変換し、振動の周波数に対する振動の振幅レベルを示す振幅特性データと、前記振動の周波数に対する振動の位相を示す位相特性データをそれぞれ作成し、
前記上下方向の振動と前記左右方向の振動に対して、前記各振幅特性データから、それぞれ少なくとも１次共振周波数を求め、
前記求められた各１次共振周波数を中心周波数としたバンドパスフィルタを設定し、
前記各バンドパスフィルタを、それぞれ前記各振幅特性データにかけた後、それぞれ逆高速フーリエ変換し、１次共振周波数における時間波形データを、それぞれ求め、
前記１次共振周波数の時間波形データについて、山側と谷側のピークのピーク時間およびピークレベルを求め、
前記上下方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＵｐとし、前記左右方向の前記１次共振周波数の時間波形データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出し、
前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求め、
前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する処理装置とを有することを特徴とする棒状体の評価システム。 A measuring device that fixes one end of the rod-shaped body and measures the vertical vibration and the horizontal vibration at the other end when the other end is vibrated simultaneously for a predetermined time;
A data acquisition device that acquires time series data of vertical vibrations and time series data of horizontal vibrations measured by the measurement device;
About the measured time series data of the vertical vibration and the time series data of the left and right vibration, respectively, fast Fourier transform, amplitude characteristic data indicating the amplitude level of the vibration with respect to the vibration frequency, and the frequency of the vibration Create phase characteristic data indicating the phase of vibration with respect to
For each of the vibration in the vertical direction and the vibration in the horizontal direction, obtain at least a primary resonance frequency from each amplitude characteristic data,
A bandpass filter having each of the obtained primary resonance frequencies as a center frequency is set;
After applying each bandpass filter to each amplitude characteristic data, respectively, inverse fast Fourier transform, respectively, to obtain time waveform data at the primary resonance frequency, respectively,
About the time waveform data of the primary resonance frequency, the peak time and peak level of the peak on the peak side and the valley side are obtained,
The absolute value of the peak value in the time waveform data of the primary resonance frequency in the vertical direction is Up, the absolute value of the peak value in the time waveform data of the primary resonance frequency in the horizontal direction is Lp, and each peak value The absolute value Lp of the peak value in the left-right direction at the peak time closest to the absolute value Up of each peak value is extracted with respect to the absolute value Up of
A square root value Vp of the absolute value Up of the peak value and the absolute value Lp of the extracted peak value is obtained;
An evaluation system for a rod-shaped body, comprising: a processing device that calculates an attenuation rate using the obtained value Vp and calculates a loss coefficient based on the attenuation rate. 棒状体の一端部を固定し、他端側を振動させたときの前記他端部の上下方向の振動と左右方向の振動とを同時に所定の時間計測する計測装置と、
前記計測装置で計測された上下方向の振動の時系列データおよび左右方向の振動の時系列データを取得するデータ取得装置と、
前記上下方向の振動の時系列データ、および前記左右方向の振動の時系列データについて、それぞれ山側と谷側のピークのピーク時間およびピークレベルを求め、
前記上下方向の振動の時系列データにおけるピーク値の絶対値をＵｐとし、前記左右方向の振動の時系列データにおけるピーク値の絶対値をＬｐとし、前記各ピーク値の絶対値Ｕｐを基準として、各ピーク値の絶対値Ｕｐに最も近いピーク時間の左右方向のピーク値の絶対値Ｌｐを、それぞれ抽出し、
前記ピーク値の絶対値Ｕｐと前記抽出されたピーク値の絶対値Ｌｐとの自乗平方根の値Ｖｐを求め、
前記求められた値Ｖｐを用いて減衰率を算出し、前記減衰率に基づいて損失係数を算出する処理装置とを有することを特徴とする棒状体の評価システム。 A measuring device that fixes one end of the rod-shaped body and measures the vertical vibration and the horizontal vibration at the other end when the other end is vibrated simultaneously for a predetermined time;
A data acquisition device that acquires time series data of vertical vibrations and time series data of horizontal vibrations measured by the measurement device;
For the time series data of the vibration in the vertical direction and the time series data of the vibration in the left and right direction, the peak time and peak level of the peak on the mountain side and the valley side are obtained,
The absolute value of the peak value in the time series data of the vertical vibration is Up, the absolute value of the peak value in the time series data of the horizontal vibration is Lp, and the absolute value Up of each peak value is a reference. Extracting the absolute value Lp of the peak value in the horizontal direction of the peak time closest to the absolute value Up of each peak value, respectively,
A square root value Vp of the absolute value Up of the peak value and the absolute value Lp of the extracted peak value is obtained;
An evaluation system for a rod-shaped body, comprising: a processing device that calculates an attenuation rate using the obtained value Vp and calculates a loss coefficient based on the attenuation rate. 前記計測装置は、前記他端部の上下方向の加速度および左右方向の加速度を計測する加速度センサを有する請求項７に記載の棒状体の評価システム。   The rod-like body evaluation system according to claim 7, wherein the measuring device includes an acceleration sensor that measures the vertical acceleration and the horizontal acceleration of the other end portion. 前記計測装置は、前記他端部の上下方向の変位および左右方向の変位を計測するレーザ変位計を有する請求項８に記載の棒状体の評価システム。   The rod-like body evaluation system according to claim 8, wherein the measuring device includes a laser displacement meter that measures a vertical displacement and a horizontal displacement of the other end portion.
JP2009090243A 2009-04-02 2009-04-02 Method for evaluating rod-shaped body and system for evaluating rod-shaped body Active JP5229073B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009090243A JP5229073B2 (en) 2009-04-02 2009-04-02 Method for evaluating rod-shaped body and system for evaluating rod-shaped body

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009090243A JP5229073B2 (en) 2009-04-02 2009-04-02 Method for evaluating rod-shaped body and system for evaluating rod-shaped body

Publications (2)

Publication Number Publication Date
JP2010243251A JP2010243251A (en) 2010-10-28
JP5229073B2 true JP5229073B2 (en) 2013-07-03

Family

ID=43096424

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009090243A Active JP5229073B2 (en) 2009-04-02 2009-04-02 Method for evaluating rod-shaped body and system for evaluating rod-shaped body

Country Status (1)

Country Link
JP (1) JP5229073B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6438746B2 (en) * 2014-11-14 2018-12-19 株式会社熊谷組 Estimation method for inter-building displacement
JP6438745B2 (en) * 2014-11-14 2018-12-19 株式会社熊谷組 Estimation method for inter-building displacement
JP6721306B2 (en) * 2015-09-10 2020-07-15 株式会社ミツトヨ Hardness tester
CN106644334B (en) * 2016-09-14 2019-02-15 中国汽车工业工程有限公司 Hydraulic vibration gen exciting force tri linear logit analysis method
JPWO2022172514A1 (en) * 2021-02-10 2022-08-18

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62232559A (en) * 1986-04-01 1987-10-13 Osaka Pref Gov Method and apparatus and measuring loss factor of string
JPH0785047B2 (en) * 1987-04-06 1995-09-13 大阪府 Measuring method and measuring device for loss coefficient, dynamic elastic coefficient, dynamic shear elastic coefficient and dynamic Poisson's ratio of rod or plate
JPH11287702A (en) * 1998-03-31 1999-10-19 Fujikura Rubber Ltd Apparatus for measuring natural frequency of golf club
JP4557202B2 (en) * 2003-04-01 2010-10-06 花王株式会社 Polymer concrete composition for audio equipment
JP4121426B2 (en) * 2003-07-18 2008-07-23 協同油脂株式会社 Method and apparatus for measuring coefficient for vibration energy loss
JP3802534B2 (en) * 2004-01-30 2006-07-26 東エン株式会社 Method and apparatus for measuring sound absorption coefficient of road pavement
JP2006275557A (en) * 2005-03-28 2006-10-12 Jtekt Corp Vibration measuring structure and method for measuring vibration using this
JP4715354B2 (en) * 2005-07-20 2011-07-06 株式会社ジェイテクト Bearing device
JP4399438B2 (en) * 2006-06-16 2010-01-13 株式会社日立製作所 Elevator equipment
JP5046613B2 (en) * 2006-11-01 2012-10-10 株式会社日立製作所 Elevator equipment
JP5176464B2 (en) * 2007-10-04 2013-04-03 横浜ゴム株式会社 Method for evaluating vibration damping performance of objects

Also Published As

Publication number Publication date
JP2010243251A (en) 2010-10-28

Similar Documents

Publication Publication Date Title
JP5229073B2 (en) Method for evaluating rod-shaped body and system for evaluating rod-shaped body
CN110307994B (en) Abnormal sound detection device and abnormal sound detection method
US10180410B2 (en) Ultrasonic test system, ultrasonic test method and aircraft structural object
JP5924929B2 (en) Finger joint damage judgment method
US5571966A (en) Method and apparatus for predicting lifetime of measured object
CN111649816A (en) Structural vibration mode testing system and method based on digital image recognition
JP5176464B2 (en) Method for evaluating vibration damping performance of objects
JP5404218B2 (en) Method and apparatus for measuring crack condition of concrete structure
JP5280407B2 (en) Damage detection method, damage detection apparatus and program for columnar structure
JP6973282B2 (en) Vibration test device and vibration test method
JP5979755B2 (en) Finger joint damage judgment method
JP4121426B2 (en) Method and apparatus for measuring coefficient for vibration energy loss
JP3895573B2 (en) Elastic wave propagation velocity measurement calculation method and nondestructive compressive strength test apparatus using the method
JP5298542B2 (en) Golf club shaft evaluation method
JP4606236B2 (en) Weight measurement method
JP4033119B2 (en) Material testing method, material testing machine
JP6933613B2 (en) Temperature distribution measuring device and temperature distribution measuring method
CN113008174A (en) Electromagnetic ultrasonic sound time measuring method and device
JP3055788B2 (en) Vibration characteristic analysis method and device
Asada Method to visualize the effect of violin string pressure on the body structure by computed tomography measurements
Moreno-Oliva et al. Analysis of the frame rate limit for the estimation of the vibration measurements in a mechanical system using optical techniques
CN114791319B (en) Detection method for structure with natural frequency
JP2015158441A (en) Image creation device for structure deterioration diagnosis
JP7297696B2 (en) Detection device and focusing member
JP2000055776A (en) Measuring method for bending rigidity of structure

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120402

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130208

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130219

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130304

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160329

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5229073

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250