JP6765754B2 - Deterioration judgment method for structures - Google Patents

Deterioration judgment method for structures Download PDF

Info

Publication number
JP6765754B2
JP6765754B2 JP2016167142A JP2016167142A JP6765754B2 JP 6765754 B2 JP6765754 B2 JP 6765754B2 JP 2016167142 A JP2016167142 A JP 2016167142A JP 2016167142 A JP2016167142 A JP 2016167142A JP 6765754 B2 JP6765754 B2 JP 6765754B2
Authority
JP
Japan
Prior art keywords
sound
frequency spectrum
data
deterioration
item
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
JP2016167142A
Other languages
Japanese (ja)
Other versions
JP2018036073A (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.)
Kumagai Gumi Co Ltd
Original Assignee
Kumagai Gumi 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 Kumagai Gumi Co Ltd filed Critical Kumagai Gumi Co Ltd
Priority to JP2016167142A priority Critical patent/JP6765754B2/en
Publication of JP2018036073A publication Critical patent/JP2018036073A/en
Application granted granted Critical
Publication of JP6765754B2 publication Critical patent/JP6765754B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Description

本発明は、例えば、橋梁やトンネルの内壁などの構造物が劣化しているか否かを判定する方法に関するものである。 The present invention relates to, for example, a method for determining whether or not a structure such as a bridge or an inner wall of a tunnel has deteriorated.

従来、橋梁やトンネルの内壁などのコンクリート構造物にひび割れや空隙などの非健全部があるか否かを検査する方法として、構造物表面をハンマーなどで打撃して加振したときに発生する音を作業者が聞いて、非健全部があるか否かを判断していた。しかし、このような方法では、加振力が一定でないだけでなく、非健全部があるか否かの判断が作業者により異なる場合があるといった問題点があった。
そこで、ハンマー加振したときに発生する音をマイクロフォンで採取し、採取された音の周波数スペクトルと予め記憶しておいた健全部の周波数スペクトルである基準周波数スペクトルを比較することで、コンクリート構造物の剥離の有無を判定する方法が提案されている(例えば、特許文献1参照)。
この特許文献1では、ハンマーとマイクロフォンとを一体に構成した打撃・集音装置と、この打撃・集音装置をトンネルの内壁近傍に移動させる台車とを備えたコンクリートの剥離診断装置をトンネル内に設置するとともに、打撃・集音装置をトンネルの内壁に沿って移動させることで、トンネルの内壁のコンクリートの剥離を検出するようにしている。
Conventionally, as a method of inspecting concrete structures such as bridges and tunnel inner walls for unhealthy parts such as cracks and voids, the sound generated when the surface of the structure is hit with a hammer and vibrated. The worker listened to and judged whether there was an unhealthy part. However, such a method has a problem that not only the exciting force is not constant, but also the judgment as to whether or not there is an unhealthy part may differ depending on the worker.
Therefore, by collecting the sound generated when the hammer is vibrated with a microphone and comparing the frequency spectrum of the collected sound with the reference frequency spectrum, which is the frequency spectrum of the sound part stored in advance, the concrete structure A method for determining the presence or absence of peeling has been proposed (see, for example, Patent Document 1).
In Patent Document 1, a concrete peeling diagnostic device including a striking / sound collecting device in which a hammer and a microphone are integrally formed and a trolley for moving the striking / sound collecting device near the inner wall of the tunnel is installed in the tunnel. Along with the installation, the striking / sound collecting device is moved along the inner wall of the tunnel to detect the peeling of concrete on the inner wall of the tunnel.

特開2001−249117号公報Japanese Unexamined Patent Publication No. 2001-249117

しかしながら、前記従来のコンクリートの剥離診断装置は、設置や移動が困難な場所では使用できないため、橋梁や河川に渡された道路などのように足場が安定しない場所では、作業者が構造物表面をハンマーなどで打撃しなければならなかった。
作業者が加振する場合には、加振力を一定にすることが困難であるため、前記従来の方法のように、採取した音の周波数スペクトルを用いて剥離診断したとしても、診断精度が低下してしまうといった問題点があった。
However, since the conventional concrete peeling diagnostic device cannot be used in a place where it is difficult to install or move, an operator can touch the surface of the structure in a place where the scaffolding is not stable, such as a bridge or a road passed to a river. I had to hit it with a hammer.
When the operator vibrates, it is difficult to keep the vibrating force constant. Therefore, even if the peeling diagnosis is performed using the frequency spectrum of the collected sound as in the conventional method, the diagnosis accuracy is high. There was a problem that it would decrease.

本発明は、従来の問題点に鑑みてなされたもので、加振力がばらついた場合でも、構造物が劣化しているか否かを確実に判定する方法を提供することを目的とする。 The present invention has been made in view of the conventional problems, and an object of the present invention is to provide a method for surely determining whether or not a structure has deteriorated even when the exciting force varies.

本発明は、構造物の表面を加振したときに構造物が発生する発生音の音圧信号の周波数スペクトルから、前記構造物が劣化しているか否かを判定する方法であって、前記構造物の基準健全領域を設定するステップと、前記基準健全領域内の複数の加振点の発生音の音圧信号を採取するステップと、前記基準健全領域内の加振点の音圧信号の周波数スペクトルを算出するステップと、前記算出された周波数スペクトルを、前記周波数スペクトルのパワー平均の値が所定の値になるように基準化した基準化健全周波数スペクトルを求めるステップと、前記基準化健全部周波数スペクトルから、加振点毎に判定項目のデータを算出して基準データ群を構築するステップと、前記判定項目のデータから前記判定項目のデータ間の相関関数を算出するステップと、前記構造物の劣化判定箇所である計測点における発生音の音圧信号を採取するステップと、前記計測点の音圧信号の周波数スペクトルを算出するステップと、前記算出された音圧信号の周波数スペクトルを、前記周波数スペクトルのパワー平均の値が所定の値になるように前記周波数スペクトルを基準化した基準化計測部周波数スペクトルを求めるステップと、前記基準化計測部周波数スペクトルから、前記判定項目のデータである計測データを算出するステップと、前記計測データと前記相関関数とを用いて、前記計測データと前記基準データ群との間のマハラノビスの距離を算出するステップと、前記算出されたマハラノビスの距離を用いて前記計測点の劣化判定を行うステップと、を備え、前記基準健全領域が、前記構造物の表面を加振したときの発生音の高低、もしくは、音色が変化しない連続した領域であり、前記判定項目が、前記基準化された周波数スペクトルの音圧レベル、1/3オクターブバンドの各レベル、周波数スペクトルの変化量と存在量のいずれか、または、複数であることを特徴とする。
このように、計測する構造物を加振したときに音が変わらない領域を基準健全領域とし、この基準健全領域のデータを基準データ群とするとともに、計測データと基準データ群との間のマハラノビスの距離により構造物が劣化しているか否かを判定するようにしたので、構造物の劣化判定を容易にかつ確実に行うことができる。
また、周波数スペクトルのパワー平均の値が所定の値になるように基準化した周波数スペクトルを用いて構造物が劣化しているか否かを判定したので、加振力がばらついた場合でも、構造物が劣化しているか否かを精度よく判定することができる。
また、本発明は、前記所定の値を、予め設定された前記構造物の基準健全領域内の加振点の発生音の周波数スペクトルのうちの予め設定された周波数範囲のパワー平均の値としたことを特徴とする。
これにより、前記所定の値が、劣化判定を行う現場の構造物の形状や配筋の情報を含んだ値となるので、構造物の劣化判定精度を更に向上させることができる。
また、前記周波数スペクトル(基準化された周波数スペクトルも含む)の周波数範囲、及び、前記予め設定された周波数範囲を、800Hz〜10000Hzとしたので、道路交通量の影響により生じる不要な低周波成分を排除できるので、劣化判定精度を向上させることができるとともに、劣化判定を効率よく行うことができる。
また、前記周波数スペクトルの周波数範囲、及び、前記予め設定された周波数範囲を、3000Hz〜6500Hzとすれば、劣化判定精度を更に向上させることができる。
The present invention is a method for determining whether or not the structure is deteriorated from the frequency spectrum of the sound pressure signal of the generated sound generated by the structure when the surface of the structure is vibrated. A step of setting a reference sound region of an object, a step of collecting sound pressure signals of sounds generated by a plurality of vibration points in the reference sound region, and a frequency of a sound pressure signal of vibration points in the reference sound region. A step of calculating a spectrum, a step of obtaining a standardized sound part frequency spectrum obtained by standardizing the calculated frequency spectrum so that the power average value of the frequency spectrum becomes a predetermined value, and the standardized sound part. A step of calculating the data of the determination item for each excitation point from the frequency spectrum to construct a reference data group, a step of calculating the correlation function between the data of the determination item from the data of the determination item, and the structure. The step of collecting the sound pressure signal of the generated sound at the measurement point, which is the deterioration determination point of the above, the step of calculating the frequency spectrum of the sound pressure signal of the measurement point, and the frequency spectrum of the calculated sound pressure signal are described above. A step of obtaining a standardized measurement unit frequency spectrum in which the frequency spectrum is standardized so that the power average value of the frequency spectrum becomes a predetermined value, and a measurement which is data of the determination item from the standardized measurement unit frequency spectrum. Using the step of calculating the data, the step of calculating the distance of the maharanobis between the measurement data and the reference data group using the measurement data and the correlation function, and the step of calculating the distance of the maharanobis calculated. A step of determining deterioration of the measurement point is provided, and the reference sound region is a continuous region in which the pitch of the generated sound when the surface of the structure is vibrated or the tone color does not change, and the determination is made. The item is characterized in that the sound pressure level of the standardized frequency spectrum, each level of the 1/3 octave band, the amount of change and the amount of abundance of the frequency spectrum, or a plurality of items.
In this way, the region where the sound does not change when the structure to be measured is vibrated is set as the reference sound region, the data in this reference sound region is used as the reference data group, and the Mahalanobis between the measurement data and the reference data group is used. Since it is determined whether or not the structure is deteriorated based on the distance between the two, the deterioration of the structure can be easily and surely determined.
In addition, since it was determined whether or not the structure was deteriorated by using the frequency spectrum standardized so that the power average value of the frequency spectrum became a predetermined value, the structure was determined even if the exciting force varied. Can be accurately determined whether or not the frequency is deteriorated.
Further, in the present invention, the predetermined value is set as the value of the power average of the preset frequency range in the frequency spectrum of the sound generated at the vibration point in the reference sound region of the structure. It is characterized by that.
As a result, the predetermined value becomes a value including information on the shape and bar arrangement of the structure at the site where the deterioration determination is performed, so that the deterioration determination accuracy of the structure can be further improved.
Further, since the frequency range of the frequency spectrum (including the standardized frequency spectrum) and the preset frequency range are set to 800 Hz to 10000 Hz, unnecessary low frequency components generated by the influence of road traffic volume can be obtained. Since it can be eliminated, the deterioration determination accuracy can be improved and the deterioration determination can be efficiently performed.
Further, if the frequency range of the frequency spectrum and the preset frequency range are set to 3000 Hz to 6500 Hz, the deterioration determination accuracy can be further improved.

また、前記劣化判定を行うステップで、予め求めておいた劣化していない試験体の複数個所を加振したときの発生音のデータを用いて算出したマハラノビスの距離に基づいて設定した閾値を用いて前記計測点の劣化判定を行うようにしたので、劣化判定の精度をより高めることができる。
また、前記劣化判定を行うステップで、予め求めておいた内部に空洞のある試験体を加振したときの発生音のデータを用いて算出したマハラノビスの距離に基づいて設定した閾値と、ひび割れによる浮きが存在する試験体を加振したときの発生音のデータを用いて算出したマハラノビスの距離に基づいて設定した閾値とのいずれか一方、または、両方を用いて前記計測点の劣化判定を行うようにしたので、計測点の劣化状態についても判定できる。
In addition, in the step of determining the deterioration, a threshold value set based on the Mahalanobis distance calculated using the data of the sound generated when a plurality of non-deteriorated test specimens are vibrated is used. Since the deterioration determination of the measurement point is performed, the accuracy of the deterioration determination can be further improved.
Further, in the step of determining the deterioration, the threshold value set based on the Mahalanobis distance calculated using the data of the sound generated when the test body having a cavity inside is vibrated, which is obtained in advance, and the crack are used. Deterioration of the measurement point is determined using either or both of the threshold values set based on the Mahalanobis distance calculated using the data of the sound generated when the test body having a float is vibrated. Therefore, it is possible to determine the deterioration state of the measurement point.

なお、前記発明の概要は、本発明の必要な全ての特徴を列挙したものではなく、これらの特徴群のサブコンビネーションもまた、発明となり得る。 The outline of the present invention does not enumerate all the necessary features of the present invention, and subcombinations of these feature groups can also be inventions.

本実施の形態に係るコンクリート構造物の劣化判定システムの構成を示す図である。It is a figure which shows the structure of the deterioration determination system of the concrete structure which concerns on this embodiment. 基準健全領域内の複数点を加振したときの周波数スペクトルのバラつきを示す図である。It is a figure which shows the variation of the frequency spectrum when a plurality of points in a reference sound region are vibrated. 周波数スペクトルを基準化する方法を説明するための図である。It is a figure for demonstrating the method of standardizing a frequency spectrum. 非健全部の周波数スペクトルの一例を示す図である。It is a figure which shows an example of the frequency spectrum of an unhealthy part. コンクリート構造物の劣化判定方法を示すフローチャートである。It is a flowchart which shows the deterioration determination method of a concrete structure. 1/3オクターブバンドの各レベルの算出方法を示す図である。It is a figure which shows the calculation method of each level of a 1/3 octave band. 基準データ群と正規化された基準データ群を示す図である。It is a figure which shows the reference data group and the normalized reference data group. 空隙部を有する試験体のマハラノビスの距離の分布例を示す図である。It is a figure which shows the distribution example of the Mahalanobis distance of the test body which has a void part. 周波数スペクトルにおける存在量と変化量の抽出方法を示す図である。It is a figure which shows the extraction method of the abundance amount and the change amount in a frequency spectrum.

図1は、本実施の形態に係るコンクリート構造物の劣化判定システム(以下、劣化判定システムという)1の構成を示す図で、劣化判定システム1は、加振手段としてのハンマー2と、音採取手段としてのマイクロフォン3と、劣化判定装置4と、表示装置5とを備える。
同図に示すように、本例では、作業者W(腕のみを記載)が、劣化診断の対象物であるコンクリート構造物としての建造物の壁部(以下、壁部6という)の表面を、ハンマー2で打撃して加振した際に発生する音(以下、発生音という)をマイクロフォン3で採取することで、壁部6の発生音の音圧信号を検出する。
劣化判定装置4は、音データ取込手段41と、記憶手段42と、周波数スペクトル算出手段43と、周波数スペクトル基準化手段44と、統計データ作成手段45と、相関行列生成手段46と、劣化判定手段47とを備える。
劣化判定装置4の記憶手段42〜劣化判定手段47までの各手段は、例えば、パーソナルコンピュータのソフトウェアとメモリーとにより構成される。
音データ取込手段41は、増幅器41aとA/D変換器41bとを備える。
増幅器41aは、マイクロフォン3で測定された音圧信号から高周波ノイズ成分を除去するとともに、音圧信号を増幅する。A/D変換器41bは、増幅された音圧信号をA/D変換し、A/D変換された音圧信号である音圧波形データを記憶手段42に送る。
記憶手段42は、基準発生音記憶部42aと計測発生音記憶部42bとを備える。
基準発生音記憶部42aは、「基準健全領域」における加振点の音圧波形データを記憶し、計測発生音記憶部42bは、劣化判定する加振点の音圧波形データ(以下、計測音データという)を記憶する。ここで、「基準健全領域」とは、劣化判定する壁部6の表面を打撃して加振した際に発生する発生音の高低(もしくは、音色)が殆ど変化しない連続した領域を指す。以下、基準健全領域における加振点の打音データを健全部データという。
なお、基準健全領域内の点を劣化判定する場合には、再度加振した基準健全領域内の加振点の音圧波形データは、計測音データとして計測発生音記憶部42bに記憶される。
FIG. 1 is a diagram showing a configuration of a deterioration determination system (hereinafter referred to as a deterioration determination system) 1 for a concrete structure according to the present embodiment. The deterioration determination system 1 includes a hammer 2 as a vibration exciting means and sound sampling. A microphone 3 as means, a deterioration determination device 4, and a display device 5 are provided.
As shown in the figure, in this example, the worker W (only the arm is described) covers the surface of the wall portion (hereinafter referred to as wall portion 6) of the building as a concrete structure which is the object of deterioration diagnosis. , The sound pressure signal of the sound generated by the wall portion 6 is detected by collecting the sound generated when the hammer 2 hits and vibrates (hereinafter referred to as the generated sound) with the microphone 3.
The deterioration determination device 4 includes sound data acquisition means 41, storage means 42, frequency spectrum calculation means 43, frequency spectrum standardization means 44, statistical data creation means 45, correlation matrix generation means 46, and deterioration determination. The means 47 is provided.
Each means from the storage means 42 to the deterioration judgment means 47 of the deterioration determination device 4 is composed of, for example, software and a memory of a personal computer.
The sound data acquisition means 41 includes an amplifier 41a and an A / D converter 41b.
The amplifier 41a removes a high-frequency noise component from the sound pressure signal measured by the microphone 3 and amplifies the sound pressure signal. The A / D converter 41b A / D-converts the amplified sound pressure signal, and sends the sound pressure waveform data, which is the A / D-converted sound pressure signal, to the storage means 42.
The storage means 42 includes a reference generated sound storage unit 42a and a measurement generated sound storage unit 42b.
The reference generated sound storage unit 42a stores the sound pressure waveform data of the vibration point in the “reference sound region”, and the measurement generated sound storage unit 42b stores the sound pressure waveform data of the vibration point for determining deterioration (hereinafter, measured sound). Data) is memorized. Here, the "reference sound region" refers to a continuous region in which the pitch (or timbre) of the generated sound generated when the surface of the wall portion 6 for which deterioration is determined is hit and vibrated is hardly changed. Hereinafter, the tapping sound data of the vibration point in the reference sound region is referred to as sound part data.
When determining the deterioration of a point in the reference sound region, the sound pressure waveform data of the vibration point in the reference sound region that has been vibrated again is stored in the measurement generation sound storage unit 42b as measurement sound data.

周波数スペクトル算出手段43は、音圧波形抽出部43aと、周波数スペクトル算出部43bとを備える。
音圧波形抽出部43aは、基準発生音記憶部42aと計測発生音記憶部42bとに記憶された加振点毎の音圧波形データである健全部データと計測音データとを抽出し、周波数スペクトル算出部43bは、抽出した健全部データと計測音データとをそれぞれFFT処理して、周波数スペクトルを算出する。
周波数スペクトル基準化手段44は、基準パワー値算出部44aと周波数スペクトル基準化部44bとを備える。
基準パワー値算出部44aは、健全部データの周波数スペクトルのそれぞれについて、加振点毎に、所定の周波数範囲(ここでは、特定の卓越周波数を持たない、比較的高い周波数領域である3000Hz〜6500Hz)の音圧信号のパワー平均(エネルギー平均)を求めるとともに、これらパワー平均の平均値である基準パワー値を算出する。
周波数スペクトル基準化部44bは、健全部データの周波数スペクトルと計測音データの周波数スペクトルを、基準パワー値算出部44aで算出された基準パワー値を用いて基準化する。
ところで、800Hz未満の低周波数領域の音圧信号は、道路交通量の影響などのため、不要な成分が入る恐れがあり、10000Hz以上の高周波領域では、健全部と非健全部との差がないので、フーリエ変換において対象とする周波数範囲は800Hz〜10000Hzの範囲とすることが好ましい。なお、本例では、基準化する周波数スペクトルの範囲を前記の所定の周波数範囲よりも広い890Hz〜7070Hzとした。
The frequency spectrum calculation means 43 includes a sound pressure waveform extraction unit 43a and a frequency spectrum calculation unit 43b.
The sound pressure waveform extraction unit 43a extracts sound pressure waveform data for each excitation point stored in the reference generated sound storage unit 42a and the measurement generated sound storage unit 42b, and the measured sound data, and extracts the frequency. The spectrum calculation unit 43b calculates the frequency spectrum by performing FFT processing on the extracted sound part data and the measurement sound data, respectively.
The frequency spectrum standardizing means 44 includes a reference power value calculating unit 44a and a frequency spectrum standardizing unit 44b.
The reference power value calculation unit 44a has a predetermined frequency range (here, 3000 Hz to 6500 Hz, which is a relatively high frequency region having no specific dominant frequency) for each vibration point for each frequency spectrum of the sound part data. ), The power average (energy average) of the sound pressure signal is obtained, and the reference power value, which is the average value of these power averages, is calculated.
The frequency spectrum standardization unit 44b standardizes the frequency spectrum of the sound unit data and the frequency spectrum of the measurement sound data using the reference power value calculated by the reference power value calculation unit 44a.
By the way, the sound pressure signal in the low frequency region of less than 800 Hz may contain unnecessary components due to the influence of road traffic, etc., and there is no difference between the healthy part and the unhealthy part in the high frequency region of 10,000 Hz or more. Therefore, the target frequency range in the Fourier transform is preferably in the range of 800 Hz to 10000 Hz. In this example, the frequency spectrum to be standardized is set to 890 Hz to 7070 Hz, which is wider than the predetermined frequency range.

ここで、周波数スペクトルの基準化について説明する。
図2は、基準健全領域内の表面136点をハンマー2で加振したときの周波数スペクトルを示す図である。同図の、上側の細い実線が最大値(max)、下側の細い実線が最小値(min)、中央の太い実線が平均値(ave)である。また、上側の破線が平均値+2×標準偏差(ave+2σ)、下側の破線が平均値−2×標準偏差(ave−2σ)である。
このグラフから、平均値±2×標準偏差(95.5%の範囲に相当)における周波数スペクトルの「形」は概ね保たれていることが分かる。
したがって、加振音を基準化しても、周波数スペクトルの形、すなわち、加振音の周波数特性は同一であると見做すことができる。
周波数スペクトルの基準化は、周波数スペクトルの大きさを、当該周波数スペクトルのパワー平均値(Power Average)を基準とした相対値とすることである。すなわち、図3に示すように、算出したパワー平均値の値を、同図の四角い枠で囲った基準パワー値(例えば、50dB)とすれば周波数スペクトルを基準化できる。なお、基準パワー値は、上記のように、基準健全領域の加振音のパワー平均の平均値である。
図3に示すように、健全部では、対象としている周波数範囲(3000Hz〜6500Hz)では、大きく卓越した周波数成分は見られず、周波数特性は概ね平坦となるが、図4(a)に示す空隙部や、図4(b)に示す浮き部などの非健全部では、卓越した周波数成分が見られる。
したがって、周波数スペクトルの基準化すれば、健全部と非健全部の周波数スペクトルの「形」の差をより明確にできるとともに、加振力によるレベル変動に関係なく加振点が健全部か、非健全部(空隙部もしくは浮き部)かを判定することができる。このことは、周波数範囲を890Hz〜7070Hzとしても同様である。
Here, the standardization of the frequency spectrum will be described.
FIG. 2 is a diagram showing a frequency spectrum when 136 points on the surface in the reference sound region are vibrated by the hammer 2. In the figure, the thin solid line on the upper side is the maximum value (max), the thin solid line on the lower side is the minimum value (min), and the thick solid line in the center is the average value (ave). The upper dashed line is the mean value + 2 × standard deviation (ave + 2σ), and the lower dashed line is the mean value-2 × standard deviation (ave-2σ).
From this graph, it can be seen that the "shape" of the frequency spectrum at the mean value ± 2 × standard deviation (corresponding to the range of 95.5%) is generally maintained.
Therefore, even if the vibrato is standardized, it can be considered that the shape of the frequency spectrum, that is, the frequency characteristics of the vibrato are the same.
The standardization of the frequency spectrum is to make the magnitude of the frequency spectrum a relative value based on the power average value (Power Average) of the frequency spectrum. That is, as shown in FIG. 3, if the calculated power average value is set to the reference power value (for example, 50 dB) surrounded by the square frame in the figure, the frequency spectrum can be standardized. As described above, the reference power value is the average value of the power averages of the vibrato sounds in the reference sound region.
As shown in FIG. 3, in the sound portion, in the target frequency range (3000 Hz to 6500 Hz), no large outstanding frequency component is observed, and the frequency characteristics are generally flat, but the voids shown in FIG. 4 (a) are observed. Excellent frequency components can be seen in the unhealthy parts such as the parts and the floating parts shown in FIG. 4 (b).
Therefore, if the frequency spectrum is standardized, the difference in the "shape" of the frequency spectrum between the healthy part and the unhealthy part can be clarified, and the excitation point is the healthy part or the unhealthy part regardless of the level fluctuation due to the exciting force. It can be determined whether it is a sound part (a gap part or a floating part). This also applies when the frequency range is set to 890 Hz to 7070 Hz.

統計データ作成手段43は、周波数スペクトルの特徴を比較するための統計データを作成する。本例では、統計データとして、周波数スペクトル毎に矩形フィルタを用いて1/3オクターブバンドの各レベルを加振点P毎に求める。なお、周波数スペクトルは、周波数スペクトル基準化手段44で基準化された、健全部の加振点の周波数スペクトルである基準化健全部周波数スペクトルと、計測点の周波数スペクトルである基準化計測部周波数スペクトルである。
具体的には、基準健全領域内の加振点をPi(i=1〜n)とし、1/3オクターブバンドの中心周波数をf(j=1〜k)とすると、基準健全領域内の加振点Piにおける1/3オクターブバンドの各レベルが基準データ群を構築するデータ(基準データyij)で、劣化判定箇所における加振点Prにおける1/3オクターブバンドの各レベルが計測データyrjである。以下、劣化判定箇所における加振点Prを計測点Prという。
相関行列生成手段46は、基準データ群のデータである基準データyijを用いて、1/3オクターブバンドの各レベル間の相関行列Rを生成する。
劣化判定手段47は、マハラノビスの距離算出部47aと、判定部47bとを備える。
マハラノビスの距離算出部47aは、計測点Prのデータである1/3オクターブバンドの各レベルyrjと、相関行列Rとを用いて、計測点Prにおける発生音のデータと基準データ群とのマハラノビスの距離D2を求める。
相関行列の生成方法及びマハラノビスの距離D2の算出方法については後述する。
判定部47bは、マハラノビスの距離D2と予め設定された閾値とを比較することで、計測点Prが健全部であるか否かを判定する。
表示装置5は、劣化判定手段45の判定結果を、ディスプレイなどの表示画面5a上に表示する。
The statistical data creating means 43 creates statistical data for comparing the characteristics of the frequency spectrum. In this example, as statistical data, each level of the 1/3 octave band is obtained for each excitation point P by using a rectangular filter for each frequency spectrum. The frequency spectra are the standardized sound part frequency spectrum, which is the frequency spectrum of the excitation point of the sound part, and the standardized measurement part frequency spectrum, which is the frequency spectrum of the measurement point, standardized by the frequency spectrum standardizing means 44. Is.
Specifically, the vibration point of reference healthy region and P i (i = 1~n), when the center frequency of 1/3-octave-band and f j (j = 1~k), the reference sound region Each level of the 1/3 octave band at the vibration point P i is the data (reference data y ij ) that constructs the reference data group, and each level of the 1/3 octave band at the vibration point P r at the deterioration judgment point is The measurement data is y rj . Hereinafter, the vibration point P r in the deterioration determination portion that measurement point P r.
The correlation matrix generation means 46 uses the reference data y ij , which is the data of the reference data group, to generate the correlation matrix R between each level of the 1/3 octave band.
The deterioration determination means 47 includes a Mahalanobis distance calculation unit 47a and a determination unit 47b.
Mahalanobis distance calculating unit 47a includes a respective level y rj 1/3 octave band is the data of the measuring points P r, by using the correlation matrix R, the data and the reference data group generated sound at the measurement point P r Find the Mahalanobis distance D 2 of.
The method of generating the correlation matrix and the method of calculating the Mahalanobis distance D 2 will be described later.
The determination unit 47b determines whether or not the measurement point Pr is a sound unit by comparing the Mahalanobis distance D 2 with a preset threshold value.
The display device 5 displays the determination result of the deterioration determination means 45 on a display screen 5a such as a display.

次に、劣化判定方法について、図5のフローチャートを参照して説明する。
まず、基準健全領域を設定する(ステップS11)。
基準健全領域は、劣化判定する壁部6の表面を打撃して加振した際に発生する発生音の高低(もしくは、音色)が変化しない連続した領域を指す。以下、基準健全領域における加振点の打音データを健全部データという。
すなわち、健全部と非健全部とでは、打音の周波数特性が異なるので、予め設定された1つもしくは複数の特定周波数領域の音圧レベルを調べることで、加振点が健全部か否かを判定できる。
また、健全部における打音と非健全部における打音とでは、発生音の高低(もしくは、音色)が明らかに異なるので、加振点が健全部か否かの判定を作業者Wが行っても良い。
基準健全領域の設定が終了した後には、健全部データを取得(ステップS12)した後、取得された健全部データ(複数の加振点の音圧波形データ)の周波数スペクトルを算出する(ステップS13)。
ステップ14では、周波数スペクトル毎に、所定の周波数範囲の音圧信号のパワー平均を求めるとともに、これらパワー平均の平均値である基準パワー値を算出し、ステップ15では、この基準パワー値を用いて健全部データの周波数スペクトル基準化して基準化健全部周波数スペクトルを求める。
Next, the deterioration determination method will be described with reference to the flowchart of FIG.
First, the reference sound region is set (step S11).
The reference sound region refers to a continuous region in which the pitch (or timbre) of the generated sound generated when the surface of the wall portion 6 for which deterioration is determined is hit and vibrated does not change. Hereinafter, the tapping sound data of the vibration point in the reference sound region is referred to as sound part data.
That is, since the frequency characteristics of the tapping sound are different between the sound part and the unhealthy part, whether or not the excitation point is the sound part by examining the sound pressure level of one or a plurality of preset specific frequency regions. Can be determined.
Further, since the pitch (or timbre) of the generated sound is clearly different between the tapping sound in the sound portion and the tapping sound in the unhealthy portion, the operator W makes a determination as to whether or not the excitation point is the healthy portion. Is also good.
After the setting of the reference sound region is completed, the sound region data is acquired (step S12), and then the frequency spectrum of the acquired sound region data (sound pressure waveform data of a plurality of excitation points) is calculated (step S13). ).
In step 14, the power average of the sound pressure signal in a predetermined frequency range is obtained for each frequency spectrum, and the reference power value which is the average value of these power averages is calculated. In step 15, this reference power value is used. The frequency spectrum of the sound part data is standardized to obtain the standardized sound part frequency spectrum.

次に、相関行列を算出するための基準データ群を構築する(ステップS16)。
基準データ群は、上記したように、加振点Pi(i=1〜n)毎に求められる音圧波形データの変量であるk個の判定項目の値yij(j=1〜k)から構成される。
判定項目としては、音圧レベル、1/3オクターブバンドの各レベル、時間軸波形の変化量と存在量、周波数特性の変化量と存在量、加振点の位置と音源方向との差などが挙げられる。
本例では、1/3オクターブバンドの各レベルを判定項目とするとともに、取得された健全部データの1/3オクターブバンドの各レベルを求めることで、相関行列を算出するための基準データ群を構築する。
具体的には、図6(a)に示すように、健全部データである基準健全領域における打音の時系列波形を、125msec(サンプリング周波数16384Hz、データ長2048個)間隔で切り出し、それぞれの時系列波形をフーリエ変換して、図6(b)に示すような、周波数スペクトルを求めた後、この周波数スペクトルを基準化した基準化健全部周波数スペクトルを求め、基準化健全部周波数スペクトル毎に矩形フィルタを用いて1/3オクターブバンドの各レベルを求める。
本例では、フーリエ変換において対象とする周波数範囲を890Hz〜7070Hzとしたので、図6(c)に示すように、音圧レベルを求める1/3オクターブバンドの中心周波数fj(j=1〜k)は、1000Hz,1250Hz,1600Hz,2000Hz,2500Hz,3150Hz,4000Hz,5000Hz,6300Hzとなる。すなわち、基準データ群を構成する変量の数である判定項目の数(以下、項目数kという)は、k=9である。
また、本例では、基準データ群を構築するためのデータ数である健全部のデータ数nをn=20とした。多変量解析ではデータ数は多い方が望ましいが、健全部では変数間の相関が高いので、データ数としてはn>kであればよい。
このように、打音の取得位置であるデータ番号をi(i=1〜n)、判定項目である中心周波数fjをj(j=1〜k)とすると、基準データ群は、n×k個の基準データyijから構築される。
Next, a reference data group for calculating the correlation matrix is constructed (step S16).
Reference data groups, as described above, a vibration point P i (i = 1~n) is a variable of the sound pressure waveform data obtained for each of k judgment item values y ij (j = 1~k) Consists of.
Judgment items include sound pressure level, 1/3 octave band level, time axis waveform change and abundance, frequency characteristic change and abundance, difference between vibration point position and sound source direction, etc. Can be mentioned.
In this example, each level of the 1/3 octave band is used as a judgment item, and each level of the 1/3 octave band of the acquired sound part data is obtained to obtain a reference data group for calculating the correlation matrix. To construct.
Specifically, as shown in FIG. 6A, the time-series waveform of the tapping sound in the reference sound region, which is the sound part data, is cut out at intervals of 125 msec (sampling frequency 16384 Hz, data length 2048 pieces), and each time. After Fourier transforming the series waveform to obtain the frequency spectrum as shown in FIG. 6 (b), the standardized sound part frequency spectrum based on this frequency spectrum is obtained, and each standardized sound part frequency spectrum is rectangular. Use a filter to determine each level of the 1/3 octave band.
In this example, the target frequency range in the Fourier transform is set to 890 Hz to 7070 Hz, so as shown in FIG. 6 (c), the center frequency f j (j = 1 to 1) of the 1/3 octave band for obtaining the sound pressure level. k) is 1000 Hz, 1250 Hz, 1600 Hz, 2000 Hz, 2500 Hz, 3150 Hz, 4000 Hz, 5000 Hz, 6300 Hz. That is, the number of determination items (hereinafter referred to as the number of items k), which is the number of variables constituting the reference data group, is k = 9.
Further, in this example, n = 20 is set as the number of data n in the sound part, which is the number of data for constructing the reference data group. In multivariate analysis, it is desirable that the number of data is large, but since the correlation between variables is high in the sound part, the number of data may be n> k.
As described above, assuming that the data number that is the acquisition position of the tapping sound is i (i = 1 to n) and the center frequency f j that is the determination item is j (j = 1 to k), the reference data group is n ×. It is constructed from k reference data y ij .

次に、基準データyijを正規化して相関行列を算出する(ステップS17)。
詳細には、図7に示すように、基準データyijの判定項目毎の平均mjと標準偏差σjとを求めた後、平均mjと標準偏差σjとを用いて基準データyijを正規化する。正規化されたデータをYijとすると、Yij=(yij−mj)/σjである。正規化されたデータYijの平均は0で、標準偏差は1である。
次に、この正規化されたデータYijから、以下の[数1]に示す式を用いてYipとYiqとの相関係数rpqを算出することで、以下の[数2]に示すような、対角要素が1、p行q列の要素がrpqである、k×kの相関行列Rを算出する。

Figure 0006765754
Figure 0006765754
Next, the reference data y ij is normalized to calculate the correlation matrix (step S17).
Specifically, as shown in FIG. 7, after obtaining the average m j and the standard deviation sigma j for each determination item reference data y ij, reference data y ij by using the mean m j and the standard deviation sigma j Normalize. If the normalized data is Y ij , then Y ij = (y ij −m j ) / σ j . The mean of the normalized data Y ij is 0 and the standard deviation is 1.
Next, from this normalized data Y ij , the correlation coefficient r pq between Y ip and Y iq is calculated using the formula shown in the following [Equation 1] to obtain the following [Equation 2]. As shown, the correlation matrix R of k × k is calculated, in which the diagonal elements are 1 and the elements in p rows and q columns are r pq .
Figure 0006765754
Figure 0006765754

相関行列の算出が終了した後には、劣化判定箇所における加振点である計測点Prにおける発生音を取得して基準化計測部周波数スペクトルを求め(ステップS18)た後、ステップS19に進んで、発生音の1/3オクターブバンドの各レベルである計測データyrj(j=1〜k)を算出する。
次に、算出された計測データyrj(j=1〜k)を正規化し、正規化された計測データである1/3オクターブバンドの各レベルYrj(j=1〜k)とステップS17で算出した相関行列Rの逆行列R-1とを用いて、計測点Prにおける計測データと基準データ群とのマハラノビスの距離D2を求める(ステップS20)。
マハラノビスの距離D2は、以下の[数3]に示す式を用いて求められる。

Figure 0006765754
ここで、kは項目数である。 After calculation of the correlation matrix has been completed, after was to get the generated sound at the measuring point P r is a vibration point of degradation determination portion obtains a normalized measurement unit frequency spectrum (step S18), and proceeds to step S19 , The measurement data y rj (j = 1 to k), which is each level of the 1/3 octave band of the generated sound, is calculated.
Next, the calculated measurement data y rj (j = 1 to k) is normalized, and each level Y rj (j = 1 to k) of the 1/3 octave band, which is the normalized measurement data, and step S17. calculated by using the inverse matrix R -1 of the correlation matrix R, Mahalanobis distance D 2 between the measurement data and the reference data group at the measuring point P r (step S20).
The Mahalanobis distance D 2 is calculated using the formula shown in [Equation 3] below.
Figure 0006765754
Here, k is the number of items.

次に、ステップS20で求めたマハラノビスの距離D2と予め設定された閾値Kとを用いて、計測点Prの劣化判定を行う(ステップS21)。
2>Kである場合には計測点Prが非健全部であると判定し、D2≦Kである場合には計測点Prが健全部であると判定する。本例では、K=4とした。
判定後は、判定結果をディスプレイなどの表示画面5a上に表示する(ステップS22)。
ステップS23では、計測が終了したか否かを判定する。
計測が終了していない場合には、ステップS18に戻って、次の計測点Pr+1の発生音を計測し、計測が終了した場合には、本処理を終了する。
なお、計測点Prが非健全部であると判定された場合には、作業者Wが、スプレーなどで、計測点Prに印をつけるようにすればよい。これにより、非健全部の位置だけでなく、壁6における非健全部の分布状態についても知ることができる。
Next, the deterioration of the measurement point Pr is determined by using the Mahalanobis distance D 2 obtained in step S20 and the preset threshold value K (step S21).
When D 2 > K, it is determined that the measurement point Pr is an unhealthy part, and when D 2 ≤ K, it is determined that the measurement point Pr is a healthy part. In this example, K = 4.
After the determination, the determination result is displayed on a display screen 5a such as a display (step S22).
In step S23, it is determined whether or not the measurement is completed.
If the measurement is not completed, the process returns to step S18, the sound generated at the next measurement point Pr + 1 is measured, and if the measurement is completed, this process is terminated.
Incidentally, when the measurement point P r is determined to be unhealthy part, worker W is spray etc., it is sufficient to mark the measurement point P r. As a result, not only the position of the unhealthy part but also the distribution state of the unhealthy part on the wall 6 can be known.

[実験例]
図7は、空隙部と浮き部とを有する直方体状のコンクリート試験体を加振したときに発生する発生音のマハラノビスの距離D2の分布例を示す図で、基準データ群としては、空隙部及び浮き部から所定距離離れた健全部を加振したときの発生音のデータを用いた。
同図において、横軸は加振位置、縦軸はマハラノビスの距離D2で、D2>100のものは、D2=100とした。
図7から明らかなように、空隙部及び浮き部では、マハラノビスの距離D2が非常に大きく、計測データが基準データ群から大きく離れていることがわかる。
これにより、マハラノビスの距離D2を用いてコンクリート構造物の劣化判定を行うことができることが確認された。
[Experimental example]
FIG. 7 is a diagram showing a distribution example of the Mahalanobis distance D 2 of the sound generated when a rectangular parallelepiped concrete test piece having a gap portion and a floating portion is vibrated, and the gap portion is used as a reference data group. And the data of the sound generated when the sound part separated from the floating part by a predetermined distance was vibrated was used.
In the figure, the horizontal axis is the vibration position, the vertical axis is the Mahalanobis distance D 2 , and when D 2 > 100, D 2 = 100.
As is clear from FIG. 7, it can be seen that the Mahalanobis distance D 2 is very large in the gap and the floating portion, and the measurement data is far away from the reference data group.
As a result, it was confirmed that the deterioration of the concrete structure can be determined using the Mahalanobis distance D 2 .

以上、本発明を実施の形態及び実験例を用いて説明したが、本発明の技術的範囲は前記実施の形態に記載の範囲には限定されない。前記実施の形態に、多様な変更または改良を加えることが可能であることが当業者にも明らかである。そのような変更または改良を加えた形態も本発明の技術的範囲に含まれ得ることが、特許請求の範囲から明らかである。 Although the present invention has been described above with reference to embodiments and experimental examples, the technical scope of the present invention is not limited to the scope described in the embodiments. It will be apparent to those skilled in the art that various changes or improvements can be made to the above embodiments. It is clear from the claims that forms with such modifications or improvements may also be included in the technical scope of the invention.

例えば、前記実施の形態では、基準パワー値として、健全部の周波数スペクトルのパワー平均の平均値としたが、例えば、55dBなど、単に所定の値であってもよい。なお、この所定の値は、実験室等で空隙部や浮き部のないコンクリートの表面を複数個所加振して求めた値とすることが好ましい。
また、前記実施の形態では、建造物の壁部6の劣化状態を判定したが、本発明は橋梁やトンネルの内壁、更には建造物の躯体などの他のコンクリート構造物にも適用可能である。また、本発明は、コンクリート構造物だけでなく、モルタルやタイルなどから成る構造物の劣化診断にも適用できる。
また、前記実施の形態では、フーリエ変換において対象とする周波数範囲を890Hz〜7070Hzとしたが、健全部と非健全部の周波数スペクトルの「形」の差がより明確な3000Hz〜6500Hzとすれば、劣化判定精度を更に向上させることができる。
For example, in the above-described embodiment, the reference power value is the average value of the power average of the frequency spectrum of the sound portion, but it may be simply a predetermined value such as 55 dB. It is preferable that this predetermined value is a value obtained by vibrating a plurality of concrete surfaces having no voids or floating portions in a laboratory or the like.
Further, in the above-described embodiment, the deteriorated state of the wall portion 6 of the building is determined, but the present invention can be applied to the inner wall of a bridge or a tunnel, and further to other concrete structures such as a building frame. .. Further, the present invention can be applied not only to a concrete structure but also to a deterioration diagnosis of a structure made of mortar, tile, or the like.
Further, in the above-described embodiment, the frequency range targeted in the Fourier transform is set to 890 Hz to 7070 Hz, but if the difference in the "shape" of the frequency spectra of the sound portion and the unhealthy portion is set to 3000 Hz to 6500 Hz, which is clearer. The deterioration determination accuracy can be further improved.

また、前記実施の形態では、1/3オクターブバンドの各レベルを判定項目としたが、音圧レベル、時間軸波形の変化量と存在量、周波数特性の変化量と存在量などを判定項目としてもよい。
例えば、周波数特性の変化量と存在量とを判定項目とする場合には、健全部データの音圧信号を、サンプリング周波数16384Hz、データ長2048個(周波数分解能8Hz)でフーリエ変換して得られた周波数スペクトルを用いればよい。この場合も、道路交通音などの影響を考慮して、変化量と存在量とを求める打音の周波数スペクトルの周波数範囲は890Hz〜7070Hzとすることが好ましい。
変化量と存在量とは周波数スペクトルの特徴を抽出するために用いられるもので、図8に示すように、周波数軸(横軸)に平行な複数の標本線を設定し、周波数特性値が標本線の値よりも大きい周波数のカウント数が存在量で、周波数特性値が標本線を横切った周波数のカウント数が変化量である。
例えば、標本線の数を7本とすれば、測定項目数はk=7×2=14となる。したがって、健全部のデータ数としては、n>14であればよい。
なお、判定項目を求めるための周波数スペクトルとしては、基準化健全部周波数スペクトル及び基準化計測部周波数スペクトルを用いることはいうまでもない。
Further, in the above-described embodiment, each level of the 1/3 octave band is used as a determination item, but the sound pressure level, the amount of change and the abundance of the time axis waveform, the amount of change and the abundance of the frequency characteristic, and the like are used as the determination items. May be good.
For example, when the amount of change and the amount of abundance of frequency characteristics are used as determination items, the sound pressure signal of the sound part data is obtained by Fourier transform with a sampling frequency of 16384 Hz and a data length of 2048 (frequency resolution of 8 Hz). The frequency spectrum may be used. Also in this case, the frequency range of the frequency spectrum of the tapping sound for obtaining the amount of change and the amount of abundance is preferably set to 890 Hz to 7070 Hz in consideration of the influence of road traffic sound and the like.
The amount of change and the amount of abundance are used to extract the characteristics of the frequency spectrum. As shown in FIG. 8, a plurality of sample lines parallel to the frequency axis (horizontal axis) are set, and the frequency characteristic value is sampled. The number of counts of frequencies larger than the value of the line is the abundance amount, and the number of counts of frequencies whose frequency characteristic value crosses the sample line is the amount of change.
For example, if the number of sample lines is 7, the number of measurement items is k = 7 × 2 = 14. Therefore, the number of data in the sound part may be n> 14.
Needless to say, the standardized sound part frequency spectrum and the standardized measurement part frequency spectrum are used as the frequency spectrum for obtaining the determination item.

また、前記実施の形態では、劣化判定の閾値をK=4としたが、予め求めておいた劣化していない試験体の複数個所を加振したのデータysjを用いて算出したマハラノビスの距離Ds 2に基づいて設定した閾値Ksを用いて計測点の劣化判定を行ってもよい。閾値Ksとしては、例えば、Ds 2の平均値Dsm 2、もしくは、Ds 2の最大値Dsmax 2に、a>0である係数をかけたものを用いるなどすればよい。
また、予め求めておいた内部に空洞のある試験体を加振したときの発生音のデータyvjを用いて算出したマハラノビスの距離Dvに基づいて設定した閾値Kvと用いて、計測点の劣化状態が空洞であるか否かを判定するようにしてもよい。
また、ひび割れによる浮きが存在する試験体を加振したときの発生音のデータyfjを用いて算出したマハラノビスの距離Dに基づいて設定した閾値Kfと用いて、計測点の劣化状態が浮きであるか否かを判定することも可能である。
Further, in the above-described embodiment, the threshold value for determining deterioration is set to K = 4, but the distance of Mahalanobis calculated using the data y sj obtained in advance by vibrating a plurality of non-deteriorated test specimens. Deterioration of the measurement point may be determined using the threshold value K s set based on D s 2 . The threshold K s, for example, the average value D sm 2 of D s 2, or, to a maximum value D smax 2 of D s 2, may be such as used multiplied by the coefficient is a> 0.
In addition, the measurement point is used as the threshold value K v set based on the Mahalanobis distance D v calculated using the data y vj of the sound generated when the test body having a cavity inside is vibrated, which is obtained in advance. It may be determined whether or not the deteriorated state of the is hollow.
Further, the deterioration state of the measurement point is determined by using the threshold value K f set based on the Mahalanobis distance D f calculated by using the data y fj of the sound generated when the test body having the floating due to the crack is vibrated. It is also possible to determine whether or not it is floating.

1 コンクリート構造物の劣化判定システム、2 ハンマー、3 マイクロフォン、
4 劣化判定装置、5 表示装置、6 壁部、
41 音データ取込手段、42 記憶手段、43 周波数スペクトル算出手段、
44 周波数スペクトル基準化手段、45 統計データ作成手段、
46 相関行列生成手段、47 劣化判定手段。
1 Deterioration judgment system for concrete structures, 2 hammers, 3 microphones,
4 Deterioration judgment device, 5 Display device, 6 Wall part,
41 Sound data acquisition means, 42 Storage means, 43 Frequency spectrum calculation means,
44 Frequency spectrum standardization means, 45 Statistical data creation means,
46 Correlation matrix generation means, 47 Deterioration determination means.

Claims (6)

構造物の表面を加振したときに構造物が発生する発生音の音圧信号の周波数スペクトルから、前記構造物が劣化しているか否かを判定する方法であって、
前記構造物の基準健全領域を設定するステップと、
前記基準健全領域内の複数の加振点の発生音の音圧信号を採取するステップと、
前記基準健全領域内の加振点の音圧信号の周波数スペクトルを算出するステップと、
前記算出された周波数スペクトルを、前記周波数スペクトルのパワー平均の値が所定の値になるように基準化した基準化健全周波数スペクトルを求めるステップと、
前記基準化健全部周波数スペクトルから、加振点毎に判定項目のデータを算出して基準データ群を構築するステップと、
前記判定項目のデータから前記判定項目のデータ間の相関関数を算出するステップと、
前記構造物の劣化判定箇所である計測点における発生音の音圧信号を採取するステップと、
前記計測点の音圧信号の周波数スペクトルを算出するステップと、
前記算出された音圧信号の周波数スペクトルを、前記周波数スペクトルのパワー平均の値が所定の値になるように前記周波数スペクトルを基準化した基準化計測部周波数スペクトルを求めるステップと、
前記基準化計測部周波数スペクトルから、前記判定項目のデータである計測データを算出するステップと、
前記計測データと前記相関関数とを用いて、前記計測データと前記基準データ群との間のマハラノビスの距離を算出するステップと、
前記算出されたマハラノビスの距離を用いて前記計測点の劣化判定を行うステップと、
を備え、
前記基準健全領域が、前記構造物の表面を加振したときの発生音の高低、もしくは、音色が変化しない連続した領域であり、
前記判定項目が、前記基準化された周波数スペクトルの音圧レベル、1/3オクターブバンドの各レベル、周波数スペクトルの変化量と存在量のいずれか、または、複数であることを特徴とする構造物の劣化判定方法。
It is a method of determining whether or not the structure is deteriorated from the frequency spectrum of the sound pressure signal of the sound generated by the structure when the surface of the structure is vibrated.
The step of setting the reference healthy area of the structure and
A step of collecting sound pressure signals of sounds generated at a plurality of excitation points in the reference sound region, and
The step of calculating the frequency spectrum of the sound pressure signal at the excitation point in the reference sound region, and
The calculated frequency spectrum, the steps of the value of the mean power of the frequency spectrum to seek normalized sound unit frequency spectrum normalized to a predetermined value,
From the standardized sound part frequency spectrum, the step of calculating the data of the judgment item for each excitation point and constructing the reference data group, and
The step of calculating the correlation function between the data of the judgment item from the data of the judgment item, and
The step of collecting the sound pressure signal of the generated sound at the measurement point which is the deterioration determination point of the structure, and
The step of calculating the frequency spectrum of the sound pressure signal at the measurement point and
A step of obtaining the frequency spectrum of the standardized measurement unit based on the frequency spectrum of the calculated sound pressure signal so that the power average value of the frequency spectrum becomes a predetermined value.
A step of calculating measurement data, which is data of the determination item, from the frequency spectrum of the standardized measurement unit, and
A step of calculating the Mahalanobis distance between the measurement data and the reference data group using the measurement data and the correlation function, and
The step of determining the deterioration of the measurement point using the calculated Mahalanobis distance, and
With
The reference sound region is a continuous region in which the pitch of the sound generated when the surface of the structure is vibrated or the timbre does not change.
The structure is characterized in that the determination item is one or more of the sound pressure level of the standardized frequency spectrum, each level of the 1/3 octave band, the amount of change and the abundance of the frequency spectrum. Deterioration judgment method.
前記所定の値が、予め設定された前記構造物の基準健全領域内の加振点の発生音の周波数スペクトルのうちの予め設定された周波数範囲のパワー平均の値であることを特徴とする請求項1に記載の構造物の劣化判定方法。 A claim characterized in that the predetermined value is a power average value of a preset frequency range in a frequency spectrum of a sound generated at a vibration point within a reference sound region of the structure. Item 1. The method for determining deterioration of a structure according to Item 1. 前記周波数スペクトルの周波数範囲、及び、前記予め設定された周波数範囲を、800Hz〜10000Hzとしたことを特徴とする請求項2に記載の構造物の劣化判定方法。 Frequency range of the frequency spectrum, and wherein the predetermined frequency range, the deterioration determination method of the structure according toMotomeko 2 you, characterized in that the 800Hz~10000Hz. 前記周波数スペクトルの周波数範囲、及び、前記予め設定された周波数範囲を、3000Hz〜6500Hzとしたことを特徴とする請求項3に記載の構造物の劣化判定方法。 The method for determining deterioration of a structure according to claim 3, wherein the frequency range of the frequency spectrum and the preset frequency range are set to 3000 Hz to 6500 Hz. 前記劣化判定を行うステップでは、
予め求めておいた劣化していない試験体の複数個所を加振したときの発生音のデータを用いて算出したマハラノビスの距離に基づいて設定した閾値を用いて前記計測点の劣化判定を行うことを特徴とする請求項1〜請求項4のいずれかに記載の構造物の劣化判定方法。
In the step of performing the deterioration determination,
Deterioration judgment of the measurement point is performed using a threshold value set based on the Mahalanobis distance calculated using the data of the sound generated when a plurality of non-deteriorated test specimens have been obtained in advance. The method for determining deterioration of a structure according to any one of claims 1 to 4.
前記劣化判定を行うステップでは、
予め求めておいた内部に空洞のある試験体を加振したときの発生音のデータを用いて算出したマハラノビスの距離に基づいて設定した閾値と、ひび割れによる浮きが存在する試験体を加振したときの発生音のデータを用いて算出したマハラノビスの距離に基づいて設定した閾値とのいずれか一方、または、両方を用いて前記計測点の劣化判定を行うことを特徴とする請求項1〜請求項5のいずれかに記載の構造物の劣化判定方法。
In the step of performing the deterioration determination,
The threshold value set based on the Mahalanobis distance calculated using the data of the sound generated when the test piece with a cavity inside was obtained in advance and the test piece with the floating due to cracks were shaken. Claims 1 to claim that the deterioration determination of the measurement point is performed using either one or both of the threshold values set based on the Mahalanobis distance calculated using the data of the generated sound at the time. Item 5. The method for determining deterioration of a structure according to any one of Item 5.
JP2016167142A 2016-08-29 2016-08-29 Deterioration judgment method for structures Active JP6765754B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016167142A JP6765754B2 (en) 2016-08-29 2016-08-29 Deterioration judgment method for structures

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016167142A JP6765754B2 (en) 2016-08-29 2016-08-29 Deterioration judgment method for structures

Publications (2)

Publication Number Publication Date
JP2018036073A JP2018036073A (en) 2018-03-08
JP6765754B2 true JP6765754B2 (en) 2020-10-07

Family

ID=61565607

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016167142A Active JP6765754B2 (en) 2016-08-29 2016-08-29 Deterioration judgment method for structures

Country Status (1)

Country Link
JP (1) JP6765754B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7166230B2 (en) * 2019-08-30 2022-11-07 株式会社奥村組 Apparatus and method for acquiring tile soundness standard
JP7379974B2 (en) * 2019-09-12 2023-11-15 東ソー株式会社 Method for manufacturing metal component-polyarylene sulfide resin component composite

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1994028388A1 (en) * 1993-05-25 1994-12-08 Commonwealth Scientific And Industrial Research Organisation Method and system for detecting a discontinuity in a structure
JPH09257765A (en) * 1996-03-19 1997-10-03 Fujita Corp Apparatus for diagnosing peeling of wall surface
JPH11304773A (en) * 1998-04-21 1999-11-05 Makoto Kajitani Method for discriminating peeling of external wall tile
JP2002082101A (en) * 2000-09-08 2002-03-22 Showa Denki Kenkyusho:Kk Method for inspecting void in boundary of two-layer structure
JP2003014708A (en) * 2001-07-05 2003-01-15 Mitsubishi Electric Corp Nondestructive inspection apparatus for plate
JP4869091B2 (en) * 2007-01-30 2012-02-01 株式会社東芝 Insulation abnormality diagnosis system for electrical equipment
JP5208625B2 (en) * 2008-09-01 2013-06-12 佐藤工業株式会社 How to determine the quality type of a structure
JP6253903B2 (en) * 2013-06-25 2017-12-27 公益財団法人鉄道総合技術研究所 Structural hammering evaluation inspection system
JP6567268B2 (en) * 2014-11-18 2019-08-28 株式会社東芝 Signal processing device, server device, detection system, and signal processing method
JP6501673B2 (en) * 2015-08-20 2019-04-17 株式会社熊谷組 Method of judging deterioration of structure

Also Published As

Publication number Publication date
JP2018036073A (en) 2018-03-08

Similar Documents

Publication Publication Date Title
JP6501673B2 (en) Method of judging deterioration of structure
JP6765754B2 (en) Deterioration judgment method for structures
CN111487315A (en) Audio frequency nondestructive testing method for tunnel lining thickness and void
CN113326550B (en) Real-time verification method for vibration quality of concrete precast bridge
JP2013134090A (en) Breakage determination method for finger joint
JP5979755B2 (en) Finger joint damage judgment method
KR100707389B1 (en) Apparatus and method for evaluating integrity of concrete structure using sound wave analysis
JP2018096858A (en) Method for non-contact acoustic probing and non-contact acoustic probing system
Elie et al. Estimation of mechanical properties of panels based on modal density and mean mobility measurements
JP4598809B2 (en) Soundness diagnosis method by sound analysis
JP2549482B2 (en) Structure soundness determination device
JP2019032242A (en) Attenuation time analytic method, device, and program
US11573135B2 (en) Tension in post-tensioned rods
JP3922459B2 (en) Separation and cavity detection method and apparatus by percussion method
JP2020190460A (en) Noncontact acoustic analysis system
JPH09257765A (en) Apparatus for diagnosing peeling of wall surface
JP7166230B2 (en) Apparatus and method for acquiring tile soundness standard
US6757620B1 (en) Method for examining structures having high natural vibration frequency using alternating manual vibration-exciting method
JP2001208733A (en) Device for measuring degradation of concrete structure
JP2941150B2 (en) Rock property measurement device
JP2007051991A (en) Method and system for measuring wall thickness of metal pipe
JP3741702B2 (en) Floor impact sound level analysis method
Radhika et al. Time-Frequency analysis of acoustic signals from tyre-pavement interaction
JP7175861B2 (en) TILE HEALTH DETERMINATION DEVICE AND METHOD
CN116773661A (en) Method, system and equipment for detecting internal cavity of structure based on sounding method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190510

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200318

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200707

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200903

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: 20200915

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200915

R150 Certificate of patent or registration of utility model

Ref document number: 6765754

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350