JP7478710B2 - Noise source contribution analysis device, noise source contribution analysis method, and noise source contribution analysis program - Google Patents

Noise source contribution analysis device, noise source contribution analysis method, and noise source contribution analysis program Download PDF

Info

Publication number
JP7478710B2
JP7478710B2 JP2021130060A JP2021130060A JP7478710B2 JP 7478710 B2 JP7478710 B2 JP 7478710B2 JP 2021130060 A JP2021130060 A JP 2021130060A JP 2021130060 A JP2021130060 A JP 2021130060A JP 7478710 B2 JP7478710 B2 JP 7478710B2
Authority
JP
Japan
Prior art keywords
sound pressure
noise source
frequency band
noise
observation point
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
JP2021130060A
Other languages
Japanese (ja)
Other versions
JP2023024023A (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.)
Toyota Technical Development Corp
Original Assignee
Toyota Technical Development Corp
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 Toyota Technical Development Corp filed Critical Toyota Technical Development Corp
Priority to JP2021130060A priority Critical patent/JP7478710B2/en
Publication of JP2023024023A publication Critical patent/JP2023024023A/en
Application granted granted Critical
Publication of JP7478710B2 publication Critical patent/JP7478710B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

本発明は、騒音源寄与解析装置、騒音源寄与解析方法、及び騒音源寄与解析プログラムに関し、特に騒音対策を講ずるべき場所での騒音の音圧に、複数の騒音源の音圧がそれぞれどのように寄与しているかを解析する騒音源寄与解析装置、騒音源寄与解析方法、及び騒音源寄与解析プログラムに関するものである。 The present invention relates to a noise source contribution analysis device, a noise source contribution analysis method, and a noise source contribution analysis program, and in particular to a noise source contribution analysis device, a noise source contribution analysis method, and a noise source contribution analysis program that analyze how the sound pressure of multiple noise sources contributes to the sound pressure of noise in a location where noise control measures should be taken.

騒音対策を行う際に、騒音源が複数あるなかで観測点での騒音の音圧に最も寄与している騒音源に対して対策を講じることは有効である。観測点の騒音の音圧に対する、騒音源の音圧の寄与の度合を解析する装置としては、従来、騒音源寄与解析装置が知られている(例えば、特許文献1参照)。 When implementing noise control measures, it is effective to take measures against the noise source that contributes most to the sound pressure of the noise at the observation point among multiple noise sources. Noise source contribution analysis devices have been known as devices that analyze the degree of contribution of the sound pressure of a noise source to the sound pressure of the noise at the observation point (see, for example, Patent Document 1).

従来の騒音源寄与解析装置において、観測点での騒音の音圧は、複数の騒音源の観測点における音圧エネルギーの加算により算出され、この算出された観測点での騒音の音圧に基づいて観測点における複数の騒音源の寄与の度合を解析していた。 In conventional noise source contribution analysis devices, the sound pressure of the noise at an observation point was calculated by adding up the sound pressure energy of multiple noise sources at the observation point, and the degree of contribution of multiple noise sources at the observation point was analyzed based on this calculated sound pressure of the noise at the observation point.

上記した方法で算出された観測点での騒音の音圧は、観測点で実際に測定した音圧に対して、高周波領域では概ね一致しているが低周波領域では乖離が大きいことが分かっている。これにより従来の騒音源寄与解析装置では、観測点における複数の騒音源の寄与の度合を正確に解析することが困難になる恐れがあった。 It is known that the sound pressure of noise at an observation point calculated using the above method generally matches the sound pressure actually measured at the observation point in the high frequency range, but deviates significantly in the low frequency range. This means that with conventional noise source contribution analysis devices, it may be difficult to accurately analyze the degree of contribution of multiple noise sources at an observation point.

特開平5-26722号公報Japanese Patent Application Laid-Open No. 5-26722

そこで、本開示は、観測点における複数の騒音源の寄与の度合をより正確に解析するために必要となる正確な観測点の総音圧を推定できる騒音源寄与解析装置、騒音源寄与解析方法、及び騒音源寄与解析プログラムを提供することを目的とする。 The present disclosure therefore aims to provide a noise source contribution analysis device, a noise source contribution analysis method, and a noise source contribution analysis program that can accurately estimate the total sound pressure at an observation point, which is necessary to more accurately analyze the degree of contribution of multiple noise sources at the observation point.

すなわち、第1の態様に係る騒音源寄与解析装置は、複数の騒音源が寄与する観測点における総音圧に対し、当該騒音源ごとの騒音源音圧の寄与の度合を解析する騒音源寄与解析装置において、複数の騒音源の各々の近傍で検出された騒音源音圧を取得する騒音源音圧取得部と、複数の騒音源の各々から観測点までの騒音の伝達特性を示す伝達関数を複数の騒音源ごとに取得する伝達関数取得部と、伝達関数取得部により取得された複数の騒音源ごとの伝達関数に、当該伝達関数に対応する騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する観測点音圧算出部と、観測点音圧算出部により算出された観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する高周波帯域抽出部と、観測点音圧算出部により算出された観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する低周波帯域抽出部と、複数の騒音源の各々の高周波帯域音圧の2乗和と複数の騒音源の各々の低周波帯域音圧の和とに基づき、観測点における総音圧を算出する総音圧算出部と、を備える。 That is, the noise source contribution analysis device according to the first aspect is a noise source contribution analysis device that analyzes the degree of contribution of the noise source sound pressure of each of the multiple noise sources to the total sound pressure at an observation point to which the multiple noise sources contribute, and includes a noise source sound pressure acquisition unit that acquires the noise source sound pressure detected in the vicinity of each of the multiple noise sources, a transfer function acquisition unit that acquires, for each of the multiple noise sources, a transfer function that indicates the transfer characteristics of the noise from each of the multiple noise sources to the observation point, and a noise source sound pressure corresponding to the transfer function acquired for each of the multiple noise sources by the transfer function acquisition unit. The system includes an observation point sound pressure calculation unit that receives each of the noise sources and outputs the respective outputs as observation point sound pressures, a high frequency band extraction unit that extracts the high frequency band of the observation point sound pressures calculated by the observation point sound pressure calculation unit as high frequency band sound pressures, a low frequency band extraction unit that extracts the low frequency band of the observation point sound pressures calculated by the observation point sound pressure calculation unit as low frequency band sound pressures, and a total sound pressure calculation unit that calculates the total sound pressure at the observation point based on the sum of the squares of the high frequency band sound pressures of each of the multiple noise sources and the sum of the low frequency band sound pressures of each of the multiple noise sources.

第2の態様は、第1の態様に係る振動源寄与解析装置において、算出された総音圧に対して、複数の騒音源それぞれの寄与の度合を算出する寄与度算出部と、をさらに備えることとしてもよい。 In a second aspect, the vibration source contribution analysis device according to the first aspect may further include a contribution calculation unit that calculates the degree of contribution of each of a plurality of noise sources to the calculated total sound pressure.

第3の態様は、第1または第2の態様に係る振動源寄与解析装置において、高周波帯域と低周波帯域との境界となるカットオフ周波数は、騒音源と観測点との間の距離に基づいて設定されることとしてもよい。 In a third aspect, in the vibration source contribution analysis device according to the first or second aspect, the cutoff frequency that is the boundary between the high frequency band and the low frequency band may be set based on the distance between the noise source and the observation point.

第4の態様は、第3の態様に係る振動源寄与解析装置において、カットオフ周波数は、騒音源と観測点との間の距離を1波長とする周波数であることとしてもよい。 In a fourth aspect, in the vibration source contribution analysis device according to the third aspect, the cutoff frequency may be a frequency at which the distance between the noise source and the observation point is one wavelength.

第5の態様は、第1または第2の態様に係る振動源寄与解析装置において、カットオフ周波数は、騒音源が存在する面の寸法に基づいて設定されることとしてもよい。 In a fifth aspect, in the vibration source contribution analysis device according to the first or second aspect, the cutoff frequency may be set based on the dimensions of the surface on which the noise source is present.

第6の態様に係る騒音源寄与解析方法は、複数の騒音源が寄与する観測点における総音圧に対し、当該騒音源ごとの騒音源音圧の寄与の度合を解析する騒音源寄与解析方法において、騒音源寄与解析装置に用いられるコンピュータに、複数の騒音源の各々の近傍で検出された騒音源音圧を取得する騒音源音圧取得ステップと、複数の騒音源の各々から観測点までの騒音の伝達特性を示す伝達関数を複数の騒音源ごとに取得する伝達関数取得ステップと、伝達関数取得ステップにおいて取得された複数の騒音源ごとの伝達関数に、当該伝達関数に対応する騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する観測点音圧算出ステップと、観測点音圧算出ステップにおいて算出された観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する高周波帯域抽出ステップと、観測点音圧算出ステップにおいて算出された観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する低周波帯域抽出ステップと、複数の騒音源の各々の高周波帯域音圧の2乗和と複数の騒音源の各々の低周波帯域音圧の和とに基づき、観測点における総音圧を算出する総音圧算出ステップと、を実行させる。 In a noise source contribution analysis method according to a sixth aspect, in a noise source contribution analysis method for analyzing the degree of contribution of a noise source sound pressure for each of a plurality of noise sources to a total sound pressure at an observation point to which the plurality of noise sources contribute, a computer used in a noise source contribution analysis device includes a noise source sound pressure acquisition step for acquiring a noise source sound pressure detected in the vicinity of each of the plurality of noise sources, a transfer function acquisition step for acquiring, for each of the plurality of noise sources, a transfer function indicating the transfer characteristics of the noise from each of the plurality of noise sources to the observation point, and a noise source sound pressure corresponding to the transfer function acquired for each of the plurality of noise sources in the transfer function acquisition step. The system executes an observation point sound pressure calculation step in which each of the sound pressures is input and the output is output as the observation point sound pressure, a high frequency band extraction step in which high frequency band sound pressures from the observation point sound pressures calculated in the observation point sound pressure calculation step are extracted as high frequency band sound pressures, a low frequency band extraction step in which low frequency band sound pressures from the observation point sound pressures calculated in the observation point sound pressure calculation step are extracted as low frequency band sound pressures, and a total sound pressure calculation step in which a total sound pressure at the observation point is calculated based on the sum of the squares of the high frequency band sound pressures of each of the multiple noise sources and the sum of the low frequency band sound pressures of each of the multiple noise sources.

第7の態様に係る騒音源寄与解析プログラムにおいて、複数の騒音源が寄与する観測点における総音圧に対し、当該騒音源ごとの騒音源音圧の寄与の度合を解析する騒音源寄与解析プログラムにおいて、騒音源寄与解析装置に用いられるコンピュータに、複数の騒音源の各々の近傍で検出された騒音源音圧を取得する騒音源音圧取得機能と、複数の騒音源の各々から観測点までの騒音の伝達特性を示す伝達関数を複数の騒音源ごとに取得する伝達関数取得機能と、伝達関数取得機能により取得された複数の騒音源ごとの伝達関数に、当該伝達関数に対応する騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する観測点音圧算出機能と、観測点音圧算出機能により算出された観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する高周波帯域抽出機能と、観測点音圧算出機能により算出された観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する低周波帯域抽出機能と、複数の騒音源の各々の高周波帯域音圧の2乗和と複数の騒音源の各々の低周波帯域音圧の和とに基づき、観測点における総音圧を算出する総音圧算出機能と、を発揮させる。 In a noise source contribution analysis program according to a seventh aspect, which analyzes the degree of contribution of the noise source sound pressure of each of a plurality of noise sources to the total sound pressure at an observation point to which the plurality of noise sources contribute, a computer used in the noise source contribution analysis device includes a noise source sound pressure acquisition function for acquiring the noise source sound pressure detected in the vicinity of each of the plurality of noise sources, a transfer function acquisition function for acquiring, for each of the plurality of noise sources, a transfer function indicating the transfer characteristics of the noise from each of the plurality of noise sources to the observation point, and a transfer function for each of the plurality of noise sources acquired by the transfer function acquisition function, It performs an observation point sound pressure calculation function that inputs the corresponding noise source sound pressures and outputs the output as the observation point sound pressure, a high frequency band extraction function that extracts the high frequency band of the observation point sound pressures calculated by the observation point sound pressure calculation function as high frequency band sound pressure, a low frequency band extraction function that extracts the low frequency band of the observation point sound pressures calculated by the observation point sound pressure calculation function as low frequency band sound pressure, and a total sound pressure calculation function that calculates the total sound pressure at the observation point based on the sum of the squares of the high frequency band sound pressures of each of the multiple noise sources and the sum of the low frequency band sound pressures of each of the multiple noise sources.

本開示によれば、観測点における複数の騒音源の寄与の度合をより正確に解析するために必要となる正確な観測点の総音圧を推定できる騒音源寄与解析装置、騒音源寄与解析方法、及び騒音源寄与解析プログラムを提供することができる。 The present disclosure provides a noise source contribution analysis device, a noise source contribution analysis method, and a noise source contribution analysis program that can accurately estimate the total sound pressure at an observation point, which is necessary to more accurately analyze the degree of contribution of multiple noise sources at the observation point.

本実施形態に係る自動車の内部における騒音を説明するための図である。1 is a diagram for explaining noise inside the automobile according to the present embodiment. FIG. 本実施形態に係る耳位置(観測点)での騒音について説明する図である。1 is a diagram illustrating noise at an ear position (observation point) according to an embodiment of the present invention. FIG. 本実施形態に係る耳位置(観測点)での騒音について説明する図である。1 is a diagram illustrating noise at an ear position (observation point) according to an embodiment of the present invention. FIG. 本実施形態に係る騒音の周波数と位相差との関係を説明する図である。5 is a diagram illustrating the relationship between the frequency and phase difference of noise according to the embodiment. FIG. 本実施形態に係る位相差による音波干渉の特徴を説明するための図である。1A to 1C are diagrams for explaining characteristics of acoustic wave interference due to a phase difference according to the present embodiment. 本実施形態に係る耳位置(観測点)における、騒音の音圧の実測値と従来手法での計算値との乖離を説明するための図である。1 is a diagram for explaining the discrepancy between the actual measured sound pressure value of noise at the ear position (observation point) according to the present embodiment and the calculated value according to the conventional method. FIG. 本実施形態に係る騒音源寄与解析装置の計算値の一例を説明するための図である。5 is a diagram for explaining an example of a calculated value of the noise source contribution analysis device according to the embodiment. FIG. 本実施形態に係る騒音源寄与解析装置の機械的構成の一例を説明するためのブロック図である。FIG. 2 is a block diagram for explaining an example of a mechanical configuration of the noise source contribution analysis device according to the present embodiment. 本実施形態に係る騒音源寄与解析装置の機能的構成の一例を説明するためのブロック図である。1 is a block diagram for explaining an example of a functional configuration of a noise source contribution analysis device according to an embodiment of the present invention. 本実施形態に係る騒音源寄与解析装置の内部の処理の流れについて説明するための図である。FIG. 2 is a diagram for explaining the flow of processing inside the noise source contribution analysis device according to the present embodiment. 本実施形態に係る耳位置(観測点)における正確な総音圧が必要な理由を説明するための図である。1 is a diagram for explaining why an accurate total sound pressure at the ear position (observation point) is required according to the present embodiment. FIG. 本実施形態に係る耳位置(観測点)における正確な総音圧が必要な理由を説明するための図である。1 is a diagram for explaining why an accurate total sound pressure at the ear position (observation point) is required according to the present embodiment. FIG. 本実施形態に係る騒音源寄与解析プログラムのフローチャートである。4 is a flowchart of a noise source contribution analysis program according to the present embodiment.

図1乃至図12を参照して本実施形態に係る騒音源寄与解析装置50について説明する。
騒音源寄与解析装置50は、複数の騒音源が寄与する観測点における総音圧に対し、当該騒音源ごとの騒音源音圧の寄与の度合を解析するものである。以下、自動車10の居住空間の内部に向かって放射される騒音を例にして説明する。なお、自動車10に限定されるものではなく、鉄道、船舶、飛行機などの自動車以外の乗り物、若しくはビル、家などの建物などの内部に侵入する騒音に対しても本実施形態に係る騒音源寄与解析装置50は適用可能である。
A noise source contribution analysis device 50 according to this embodiment will be described with reference to FIGS.
The noise source contribution analysis device 50 analyzes the degree of contribution of the noise source sound pressure of each noise source to the total sound pressure at an observation point to which the multiple noise sources contribute. Hereinafter, an example of noise radiated toward the interior of the living space of an automobile 10 will be described. Note that the noise source contribution analysis device 50 according to this embodiment is not limited to automobiles 10, and can also be applied to noise that penetrates into the interior of vehicles other than automobiles, such as trains, ships, and airplanes, or buildings such as buildings and houses.

図1を参照して、自動車10の居住空間での騒音について説明する。図1は本実施形態に係る自動車10の内部における騒音を説明するための図である。
自動車10の居住空間の内部で感じる騒音には、エンジン音、風きり音、ロードノイズ、及びこもり音などがある。これらの騒音は、自動車10の居住空間を構成する壁面(パネル12)が振動することで騒音源になり、自動車10の居住空間内に向けて放射されて搭乗者の耳に騒音として届く。
Noise in the living space of an automobile 10 will be described with reference to Fig. 1. Fig. 1 is a diagram for explaining noise inside an automobile 10 according to this embodiment.
Noises felt inside the living space of the automobile 10 include engine noise, wind noise, road noise, muffled noise, etc. These noises become noise sources when the wall surfaces (panels 12) constituting the living space of the automobile 10 vibrate, and are radiated into the living space of the automobile 10 and reach the ears of the passengers as noise.

エンジン音は、自動車10のエンジンルームに設置されたエンジン14を起源とする振動・騒音が固体伝搬21及び空気伝搬22によりパネル12に伝わる。風きり音は、自動車10のボディー11の外表面に沿って流れる空気20がパネル12を振動させることによって生じる。ロードノイズは、自動車10の走行中の路面とタイヤ13との間に生じる振動・騒音が固体伝搬21及び空気伝搬22によりパネル12に伝わる。こもり音は、プロペラシャフト16に起因する振動・騒音が固体伝搬21及び空気伝搬22によりパネル12に伝わる。 Engine sound is vibration and noise originating from the engine 14 installed in the engine room of the automobile 10, which is transmitted to the panel 12 via solid propagation 21 and air propagation 22. Wind noise is caused by air 20 flowing along the outer surface of the body 11 of the automobile 10 vibrating the panel 12. Road noise is vibration and noise generated between the road surface and the tires 13 while the automobile 10 is traveling, which is transmitted to the panel 12 via solid propagation 21 and air propagation 22. Booming sound is vibration and noise caused by the propeller shaft 16 which is transmitted to the panel 12 via solid propagation 21 and air propagation 22.

パネル12に複数の騒音源が発生した場合、搭乗者の耳位置(観測点)で聞こえる騒音は複数の騒音源から放射される騒音が重ね合わさったものである。
自動車10の居住空間内の騒音対策の手段として、搭乗者の耳位置(観測点)で聞こえる騒音に大きく影響するパネル12の振動を特定し、当該振動を抑えることが有効である。そこで効果的に耳位置(観測点)で聞こえる騒音に大きく影響するパネル12を特定・対策するために、各パネル12の振動が耳位置(観測点)の騒音にどの程度の割合で寄与しているかを算出する騒音源寄与解析が実施される。
When multiple noise sources occur on the panel 12, the noise heard at the passenger's ear position (observation point) is a mixture of noises radiated from the multiple noise sources.
As a means for noise control within the living space of the automobile 10, it is effective to identify the vibrations of the panels 12 that have a large effect on the noise heard at the ear position (observation point) of the passengers and suppress the vibrations. Therefore, in order to effectively identify and take measures against the panels 12 that have a large effect on the noise heard at the ear position (observation point), a noise source contribution analysis is performed to calculate the percentage to which the vibration of each panel 12 contributes to the noise at the ear position (observation point).

次に、図2乃至図7を参照して、異なる騒音源から放射される2つの騒音の耳位置(観測点)24での重ね合わせについて説明する。高周波帯域の騒音の重ね合わせについて図2を参照して説明する。図2は、本実施形態に係る耳位置(観測点)での騒音について説明する図である。 Next, the superposition of two noises radiated from different noise sources at ear position (observation point) 24 will be described with reference to Figures 2 to 7. The superposition of high-frequency band noises will be described with reference to Figure 2. Figure 2 is a diagram explaining noises at ear position (observation point) according to this embodiment.

図2において、騒音源は第1騒音源28と第2騒音源29とがある。第1騒音源28と第2騒音源29との音圧は、それぞれの騒音源の近傍に設置された第1マイク32aと第2マイク32bによって検出される。第1マイク32aと第2マイク32bとに音響粒子速度センサのPUプローブを用いた。 In FIG. 2, there are a first noise source 28 and a second noise source 29. The sound pressures of the first noise source 28 and the second noise source 29 are detected by a first microphone 32a and a second microphone 32b installed near the respective noise sources. The first microphone 32a and the second microphone 32b are PU probes of an acoustic particle velocity sensor.

第1騒音源28、及び第2騒音源29の近傍で計測された音圧をそれぞれ音圧a1、a2とする。第1騒音源28、及び第2騒音源29から耳位置(観測点)までの伝達関数をそれぞれ伝達関数(f1)、及び(f2)とする。第1騒音源28、及び第2騒音源29から放射された騒音の耳位置(観測点)の音圧の計算値をそれぞれ音圧b1、b2と定義すると、以下の式(1)、(2)が成立する。 Let the sound pressures measured near the first noise source 28 and the second noise source 29 be sound pressures a1 and a2, respectively. Let the transfer functions from the first noise source 28 and the second noise source 29 to the ear position (observation point) be transfer functions (f1) and (f2), respectively. If the calculated sound pressures at the ear position (observation point) of the noise radiated from the first noise source 28 and the second noise source 29 are defined as sound pressures b1 and b2, respectively, then the following equations (1) and (2) hold.

Figure 0007478710000001

耳位置(観測点)の音圧を(c)と定義すると以下の式(3)が成立する。
Figure 0007478710000001

If the sound pressure at the ear position (observation point) is defined as (c), the following equation (3) holds.

Figure 0007478710000002
Figure 0007478710000002

式(3)は、高周波帯域において耳位置(観測点)での音圧エネルギーは、寄与する音のエネルギーの和にほぼ等しいことに由来する。式(3)はさらに以下の式(4)に拡張することができる。 Equation (3) comes from the fact that in the high frequency band, the sound pressure energy at the ear position (observation point) is approximately equal to the sum of the energy of the contributing sounds. Equation (3) can be further expanded to the following equation (4).

Figure 0007478710000003
Figure 0007478710000003

式(4)は、高周波帯域における耳位置(観測点)での音圧エネルギーは寄与する音圧エネルギーの総和にほぼ等しいことを意味する。これは、高周波帯域における音同士は干渉の度合が小さく独立していると見なせることによる。 Equation (4) means that the sound pressure energy at the ear position (observation point) in the high frequency band is approximately equal to the sum of the contributing sound pressure energies. This is because sounds in the high frequency band have little interference with each other and can be considered independent.

次に、低周波帯域の騒音の重ね合わせについて図3乃至図7を参照して説明する。図3は本実施形態に係る耳位置(観測点)での騒音について説明する図であり、図4は本実施形態に係る騒音の周波数と位相差との関係を説明する図であり、図5は本実施形態に係る位相差による音波干渉の特徴を説明するための図であり、図6は本実施形態に係る耳位置(観測点)における、騒音の音圧の実測値と従来手法での計算値との乖離を説明するための図であり、図7は本実施形態に係る騒音源寄与解析装置の計算値の一例を説明するための図である。 Next, the superposition of low-frequency band noise will be described with reference to Figs. 3 to 7. Fig. 3 is a diagram for explaining noise at ear position (observation point) according to this embodiment, Fig. 4 is a diagram for explaining the relationship between noise frequency and phase difference according to this embodiment, Fig. 5 is a diagram for explaining the characteristics of sound wave interference due to phase difference according to this embodiment, Fig. 6 is a diagram for explaining the deviation between the actual measured value of the sound pressure of the noise at ear position (observation point) according to this embodiment and the calculated value using a conventional method, and Fig. 7 is a diagram for explaining an example of a calculated value of the noise source contribution analysis device according to this embodiment.

第1騒音源28と第2騒音源29との間の距離が近い場合は、この2つの騒音源は1つの騒音源33とみなすことができる(図3参照)。この場合、耳位置(観測点)での騒音の音圧(c)は、以下の式(5)により算出することができる。 When the distance between the first noise source 28 and the second noise source 29 is short, these two noise sources can be regarded as one noise source 33 (see FIG. 3). In this case, the sound pressure (c) of the noise at the ear position (observation point) can be calculated by the following formula (5).

Figure 0007478710000004
Figure 0007478710000004

式(5)のPcont(yi)は、騒音源yiから放射された騒音の耳位置(観測点)における音圧の計算値であり、前述の伝達関数に騒音源yiの近傍の音圧を乗じて得た値となる。Pref(x)は、耳位置(観測点x)における音圧の実測値である。 Pcont(yi) in equation (5) is the calculated sound pressure at the ear position (observation point) of the noise radiated from noise source yi, and is the value obtained by multiplying the transfer function described above by the sound pressure near noise source yi. Pref(x) is the actual measured sound pressure at the ear position (observation point x).

低周波帯域では、耳位置(観測点)での音圧は、当該音に寄与する音同士の干渉を考慮する必要があり、各騒音の位相が同相であれば振幅の和、逆相であれば振幅の差として計算する。 In the low-frequency band, the sound pressure at ear position (observation point) must take into account the interference between the sounds that contribute to the sound in question, and is calculated as the sum of the amplitudes if the noises are in phase, and as the difference in amplitudes if they are out of phase.

さらに、第1騒音源28及び第2騒音源29と耳位置(観測点)との間の距離に対して騒音の波長が長い場合(即ち、低周波帯域の場合)についても、耳位置(観測点)での騒音の音圧(c)は、式(5)により算出することができる。 Furthermore, even when the wavelength of the noise is long relative to the distance between the first noise source 28 and the second noise source 29 and the ear position (observation point) (i.e., in the case of a low frequency band), the sound pressure (c) of the noise at the ear position (observation point) can be calculated using formula (5).

経路長の差による位相のずれの影響は、低周波帯域より高周波帯域のほうが大きいことが知られている(図4参照)。さらに、2つの騒音における相互の干渉度については、低周波帯域は高周波帯域と比較して、その影響は大きいことが知られている(図5参照)。 It is known that the effect of phase shift due to differences in path length is greater in the high frequency band than in the low frequency band (see Figure 4). Furthermore, it is known that the degree of mutual interference between two noises is greater in the low frequency band than in the high frequency band (see Figure 5).

図6を参照して従来の計算手法の問題点について説明する。図6に示す様に、従来の計算手法の計算結果(c)39と実測値(d)38とは、1400Hz以下で乖離が拡大している。従来の計算手法では、高周波帯域及び低周波帯域を含む全ての帯域で式(4)を用いて(c)を算出する。式(4)では低周波帯域における騒音相互間の干渉について考慮されておらず、これが計算結果(c)39と実測値(d)38との乖離の原因であることを本開示の発明者らは見出した。 The problems with the conventional calculation method will be explained with reference to FIG. 6. As shown in FIG. 6, the divergence between the calculation result (c) 39 and the actual measurement value (d) 38 of the conventional calculation method increases below 1400 Hz. In the conventional calculation method, (c) is calculated using equation (4) in all bands including the high-frequency band and the low-frequency band. Equation (4) does not take into account interference between noises in the low-frequency band, and the inventors of the present disclosure have found that this is the cause of the divergence between the calculation result (c) 39 and the actual measurement value (d) 38.

本開示の発明者らは、従来の計算手法を調査し、当該計算手法は各騒音源の振動による騒音が耳位置(観測点)で無相関に干渉することを前提としていることをつきとめた。耳位置(観測点)での各騒音の相互間では、各騒音源の振動による騒音は耳位置(観測点)に届くまでの経路長に差がある等の原因のため、耳位置(観測点)の各騒音の相互間に位相差ができる(図4参照)。さらに低周波帯域では各騒音の相互間の位相差が小さいため、波の重ね合わせの原理により耳位置(観測点)における各騒音の相互間の干渉度が強くなり、高周波帯域では逆に干渉度が弱くなる(図5参照)。 The inventors of the present disclosure investigated conventional calculation methods and found that the calculation methods are based on the premise that noises caused by vibrations of each noise source interfere with each other in an uncorrelated manner at the ear position (observation point). Because noises caused by vibrations of each noise source have different path lengths to reach the ear position (observation point), there is a phase difference between each noise at the ear position (observation point) (see Figure 4). Furthermore, because the phase difference between each noise is small in the low frequency band, the degree of interference between each noise at the ear position (observation point) becomes strong due to the principle of wave superposition, and conversely, the degree of interference becomes weaker in the high frequency band (see Figure 5).

そこで各騒音源での音響粒子速度および音圧の計測値を複素数として求め、騒音源の間で発生する位相差を保持して寄与度を計算する計算手法(式(5))に改修した結果、計算値71と実測値70との最大乖離幅を12dBから6dBに低減することができ、精度の向上が見られた(図7参照)。以上より、従来の解析手法で発生していた計測値72と実測値70の乖離を約6dB低減し精度向上を実現することができた。 Therefore, we modified the calculation method (equation (5)) to obtain the measured values of acoustic particle velocity and sound pressure at each noise source as complex numbers and calculate the contribution by retaining the phase difference generated between the noise sources. As a result, the maximum deviation between the calculated value 71 and the actual measured value 70 could be reduced from 12 dB to 6 dB, and accuracy was improved (see Figure 7). As a result, the deviation between the measured value 72 and the actual measured value 70, which occurred with the conventional analysis method, could be reduced by approximately 6 dB, and accuracy was improved.

次に図8を参照して、騒音源寄与解析装置50の機械的構成について説明する。本実施形態に係る騒音源寄与解析装置50の機械的構成の一例を説明するためのブロック図である。騒音源寄与解析装置50は通信ネットワーク51に接続されてもよい。
騒音源寄与解析装置50は、Read Only Memory(ROM)50a、Random Access Memory(RAM)50b、記憶部50c、Central Processing Unit(CPU)50d、入出力インターフェース50e、通信インターフェース50fなどを備えている。また、騒音源寄与解析装置50は、その外部装置として入力装置50g及び出力装置50hを備えている。
Next, a mechanical configuration of the noise source contribution analysis device 50 will be described with reference to Fig. 8. It is a block diagram for explaining an example of the mechanical configuration of the noise source contribution analysis device 50 according to the present embodiment. The noise source contribution analysis device 50 may be connected to a communication network 51.
The noise source contribution analysis device 50 includes a Read Only Memory (ROM) 50a, a Random Access Memory (RAM) 50b, a storage unit 50c, a Central Processing Unit (CPU) 50d, an input/output interface 50e, a communication interface 50f, etc. The noise source contribution analysis device 50 also includes an input device 50g and an output device 50h as its external devices.

記憶部50cは、記憶装置として利用でき、騒音源寄与解析装置50が動作する上で必要となる後述の騒音源寄与解析プログラム、各種アプリケーション及び当該アプリケーションによって利用される各種データなどが記録される。 The memory unit 50c can be used as a storage device, and stores the noise source contribution analysis program (described below) required for the operation of the noise source contribution analysis device 50, various applications, and various data used by the applications.

入出力インターフェース50eは、入力装置50g、及び出力装置50hに対してデータなどの送受信を行う。入力装置50gは、キーボード52、マウス53、Webカメラ54などのことであり、出力装置50hはモニター55、スピーカ56、及びプリンタ57などのことであり、これらはいわゆる情報処理装置の周辺機器である。 The input/output interface 50e transmits and receives data to and from the input device 50g and the output device 50h. The input device 50g includes a keyboard 52, a mouse 53, a web camera 54, etc., and the output device 50h includes a monitor 55, a speaker 56, a printer 57, etc., which are so-called peripheral devices of the information processing device.

図9を参照して、騒音源寄与解析装置50の機能的構成について説明する。図9は本実施形態に係る騒音源寄与解析装置50の機能的構成の一例を説明するためのブロック図である。
騒音源寄与解析装置50は、動作する上で必要となる騒音源寄与解析プログラムをROM50a若しくは記憶部50cに保存し、RAM50bなどで構成されるメインメモリに騒音源寄与解析プログラムを取り込む。CPU50dは、騒音源寄与解析プログラムを取り込んだメインメモリにアクセスして、騒音源寄与解析プログラムを実行する。
The functional configuration of the noise source contribution analysis device 50 will be described with reference to Fig. 9. Fig. 9 is a block diagram for explaining an example of the functional configuration of the noise source contribution analysis device 50 according to this embodiment.
The noise source contribution analysis device 50 stores a noise source contribution analysis program required for operation in a ROM 50a or a storage unit 50c, and loads the noise source contribution analysis program into a main memory constituted by a RAM 50b etc. The CPU 50d accesses the main memory into which the noise source contribution analysis program has been loaded, and executes the noise source contribution analysis program.

騒音源寄与解析装置50は、騒音源寄与解析プログラムを実行することで、CPU50dに騒音源音圧取得部60、伝達関数取得部61、観測点音圧算出部62、高周波帯域抽出部63、低周波帯域抽出部64、総音圧算出部65、及び寄与度算出部66などを機能部として備える。 By executing the noise source contribution analysis program, the noise source contribution analysis device 50 provides the CPU 50d with a noise source sound pressure acquisition unit 60, a transfer function acquisition unit 61, an observation point sound pressure calculation unit 62, a high frequency band extraction unit 63, a low frequency band extraction unit 64, a total sound pressure calculation unit 65, and a contribution degree calculation unit 66 as functional units.

騒音源音圧取得部60は、複数の騒音源の各々の近傍で検出された騒音源音圧を取得する。
騒音源音圧は、複数の騒音源の各々の近傍に設置されたマイクにより検出される。騒音源音圧取得部60は、当該マイクから騒音源音圧を取得する。
The noise source sound pressure acquisition unit 60 acquires the noise source sound pressure detected in the vicinity of each of the multiple noise sources.
The noise source sound pressure is detected by microphones installed near each of the multiple noise sources. The noise source sound pressure acquisition unit 60 acquires the noise source sound pressure from the microphones.

伝達関数取得部61は、複数の騒音源の各々から観測点までの騒音の伝達特性を示す伝達関数を複数の騒音源ごとに取得する。
伝達関数は騒音源から発せられた騒音が、観測点に伝わるまでにどのように変化しているかを示す関数である。
The transfer function acquisition unit 61 acquires, for each of the plurality of noise sources, a transfer function that indicates the transfer characteristics of noise from each of the plurality of noise sources to the observation point.
The transfer function is a function that indicates how noise emitted from a noise source changes before it is transmitted to an observation point.

観測点音圧算出部62は、伝達関数取得部61により取得された複数の騒音源ごとの伝達関数に、当該伝達関数に対応する騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する。 The observation point sound pressure calculation unit 62 inputs the noise source sound pressure corresponding to each transfer function acquired by the transfer function acquisition unit 61 for each of the multiple noise sources, and outputs the output as the observation point sound pressure.

高周波帯域抽出部63は、観測点音圧算出部62により算出された観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する。
高周波帯域音圧を抽出するカットオフ周波数は、音場の環境によって異なる。従って、総音圧算出部65の算出結果と観測点音圧の実測値との間の乖離が最小となるカットオフ周波数を予め調査し、カットオフ周波数を設定する必要がある。カットオフ周波数は、騒音源と耳位置(観測点)との間の距離に基づいて設定してもよく、当該距離を1波長とする周波数に設定してもよい。
The high frequency band extraction unit 63 extracts, from the observation point sound pressures calculated by the observation point sound pressure calculation unit 62, those in the high frequency band as high frequency band sound pressures.
The cutoff frequency for extracting the high-frequency band sound pressure varies depending on the sound field environment. Therefore, it is necessary to investigate in advance the cutoff frequency at which the deviation between the calculation result of the total sound pressure calculation unit 65 and the actual measurement value of the observation point sound pressure is the smallest, and set the cutoff frequency. The cutoff frequency may be set based on the distance between the noise source and the ear position (observation point), or may be set to a frequency that takes the distance as one wavelength.

さらにカットオフ周波数は、騒音源が存在する面の寸法に基づいて設定してもよい。 The cutoff frequency may further be set based on the dimensions of the surface on which the noise source resides.

低周波帯域抽出部64は、観測点音圧算出部62により算出された観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する。
低周波帯域音圧を抽出するカットオフ周波数は、高周波帯域音圧を抽出するカットオフ周波数と同じものを用いる。
The low-frequency band extraction unit 64 extracts low-frequency band sound pressures from the observation point sound pressures calculated by the observation point sound pressure calculation unit 62 as low-frequency band sound pressures.
The cutoff frequency for extracting the low-frequency band sound pressure is the same as the cutoff frequency for extracting the high-frequency band sound pressure.

総音圧算出部65は、複数の騒音源の各々の高周波帯域音圧の2乗和と複数の騒音源の各々の低周波帯域音圧の和とに基づき、観測点における総音圧を算出する。 The total sound pressure calculation unit 65 calculates the total sound pressure at the observation point based on the sum of the squares of the high-frequency band sound pressures of each of the multiple noise sources and the sum of the low-frequency band sound pressures of each of the multiple noise sources.

総音圧算出部65は、式(4)を用いて高周波帯域における総音圧を算出し、式(5)を用いて低周波帯域における総音圧を算出する。 The total sound pressure calculation unit 65 calculates the total sound pressure in the high frequency band using equation (4) and calculates the total sound pressure in the low frequency band using equation (5).

寄与度算出部66は、算出された総音圧に対して、複数の騒音源それぞれの寄与の度合を算出する。
寄与度算出部66は、耳位置(観測点)での総音圧に対する当該各々の音圧の割合により、各々の騒音源近傍で計測された音圧の寄与の度合を算出する。
The contribution degree calculation unit 66 calculates the degree of contribution of each of the multiple noise sources to the calculated total sound pressure.
The contribution calculation unit 66 calculates the degree of contribution of the sound pressure measured near each noise source based on the ratio of each sound pressure to the total sound pressure at the ear position (observation point).

次に図10を参照して、騒音源寄与解析装置50の処理の流れについて説明する。
図10において説明の便宜上、騒音源は4箇所とする。4箇所の騒音源近傍には第1マイク32a、第2マイク32b、第3マイク32c、及び第4マイク32dがそれぞれ設置され、4箇所の騒音源の音圧が検出される。実施形態において、高周波帯域と低周波帯域とは波長の長さにより区分される。概ね1m以上の波長の音については低周波帯域としている。これは、車内空間(図1参照)において、騒音源から運転手の耳までの距離が1mであることに基づいている。
Next, the flow of processing by the noise source contribution analysis device 50 will be described with reference to FIG.
For ease of explanation, four noise sources are shown in Fig. 10. A first microphone 32a, a second microphone 32b, a third microphone 32c, and a fourth microphone 32d are installed near the four noise sources, respectively, to detect the sound pressure of the four noise sources. In the embodiment, the high frequency band and the low frequency band are distinguished by the wavelength. Sounds with wavelengths of approximately 1 m or more are considered to be in the low frequency band. This is based on the fact that the distance from the noise source to the driver's ear in the vehicle interior space (see Fig. 1) is 1 m.

なお、高周波帯域と低周波帯域とを区分する波長の長さは、騒音源が存在する面(パネル12)の寸法に基づいて設定してもよい。一例として、騒音源が存在する面が長方形(正方形を含む)の場合、高周波帯域と低周波帯域とを区分する波長の長さは、騒音源が存在する面の長辺の長さとしてもよい。第二の例として、高周波帯域と低周波帯域とを区分する波長の長さは、騒音源が存在する面の径に基づいて、又は径と同じ長さに設定してもよい。ここで、図形の「径」とは、当該図形の両側から接する二本の平行線の間の最長距離をいい、円や楕円のみならず凸図形一般に対して定義される概念である。なお、凸図形とは、図形上の任意の2点を結ぶ線分上の点はまた図形上にある図形のことをいう。 The wavelength length that separates the high frequency band from the low frequency band may be set based on the dimensions of the surface (panel 12) on which the noise source is present. As an example, if the surface on which the noise source is present is rectangular (including square), the wavelength length that separates the high frequency band from the low frequency band may be the length of the long side of the surface on which the noise source is present. As a second example, the wavelength length that separates the high frequency band from the low frequency band may be set based on the diameter of the surface on which the noise source is present, or may be set to the same length as the diameter. Here, the "diameter" of a figure refers to the longest distance between two parallel lines that contact the figure from both sides, and is a concept that is defined not only for circles and ellipses but also for convex figures in general. Note that a convex figure refers to a figure in which a point on a line segment connecting any two points on the figure is also on the figure.

4箇所の騒音源近傍の検出された音圧は、騒音源音圧取得部60によりCPU50dに取り込まれる。次に、観測点音圧算出部62は、4箇所の騒音源近傍の検出された音圧を伝達関数f1~f4(67a~67d)に乗じることで、4箇所の騒音源から放射された騒音の観測点における音圧をそれぞれ算出する。 The detected sound pressures near the four noise sources are input to the CPU 50d by the noise source sound pressure acquisition unit 60. Next, the observation point sound pressure calculation unit 62 multiplies the detected sound pressures near the four noise sources by transfer functions f1 to f4 (67a to 67d) to calculate the sound pressures at the observation points of the noise radiated from the four noise sources.

観測点音圧算出部62により算出された4個の音圧の算出値は、高周波帯域抽出部63であるHPF(High Pass Filter)1~HPF4(63a~63d)、及び低周波帯域抽出部64であるLPF(Low Pass Filter)1~LPF4(64a~64d)へと送られる。 The four sound pressure values calculated by the observation point sound pressure calculation unit 62 are sent to the high frequency band extraction unit 63, which consists of HPF (High Pass Filter) 1 to HPF 4 (63a to 63d), and the low frequency band extraction unit 64, which consists of LPF (Low Pass Filter) 1 to LPF 4 (64a to 64d).

総音圧算出部65は、高周波帯域加算部65a及び低周波帯域加算部65bを備える。
高周波帯域加算部65aは、HPF1~HPF4(63a~63d)より抽出された高周波帯域音圧をそれぞれ(b1~b4)として式(4)に代入し高周波帯域における観測点の音圧(c)を算出する。
The total sound pressure calculation unit 65 includes a high frequency band addition unit 65a and a low frequency band addition unit 65b.
The high frequency band adder 65a substitutes the high frequency band sound pressures extracted from HPF1 to HPF4 (63a to 63d) as (b1 to b4) into equation (4) to calculate the sound pressure (c) at the observation point in the high frequency band.

低周波帯域加算部65bは、LPF1~LPF4(64a~64d)より抽出された低周波帯域音圧をそれぞれ(Pcont(y1)~Pcont(y4))として式(5)に代入し低周波帯域における観測点の音圧(c)を算出する。 The low-frequency band addition unit 65b substitutes the low-frequency band sound pressures extracted from LPF1 to LPF4 (64a to 64d) as (Pcont(y1) to Pcont(y4)) into equation (5) to calculate the sound pressure (c) at the observation point in the low-frequency band.

総音圧算出部65の算出結果である音圧P(t)は、高周波帯域における観測点の音圧(c)と低周波帯域における観測点の音圧(c)とを足し合わせることで得られる。 The sound pressure P(t), which is the calculation result of the total sound pressure calculation unit 65, is obtained by adding together the sound pressure (c) at the observation point in the high frequency band and the sound pressure (c) at the observation point in the low frequency band.

次に、図11及び図12を参照して、寄与解析の考え方について説明する。
図11及び図12は、本実施形態に係る耳位置(観測点)における正確な総音圧が必要な理由を説明するための図である。
Next, the concept of contribution analysis will be described with reference to FIGS.
11 and 12 are diagrams for explaining the reason why an accurate total sound pressure at the ear position (observation point) according to this embodiment is necessary.

騒音源寄与解析は、耳位置(観測点)での総音圧に対する各騒音源の音圧の割合を求めるだけではなく、次の観点が必要となる。すなわち、各騒音源から耳位置(観測点)に到達した騒音の音圧b1、b2、総音圧cとすると、音圧b1若しくはb2を無くした場合の総音圧cの変化について知ることである。上述した従来の計算手法と本開示の計算手法とでは、総音圧cへの影響が異なって算出される。図11及び図12において、音場1に従来と本開示との間に大差がない例を示し、一方で音場2に従来と本開示の間に大きな差が生じる例を示す。 Noise source contribution analysis does not only require finding the ratio of the sound pressure of each noise source to the total sound pressure at the ear position (observation point), but also requires the following viewpoints. That is, if the sound pressures of the noise reaching the ear position (observation point) from each noise source are b1, b2, and the total sound pressure is c, it is necessary to know the change in the total sound pressure c when sound pressure b1 or b2 is eliminated. The above-mentioned conventional calculation method and the calculation method of the present disclosure calculate the effect on the total sound pressure c differently. Figures 11 and 12 show an example in which there is not much difference between the conventional method and the present disclosure in sound field 1, while an example in which there is a large difference between the conventional method and the present disclosure in sound field 2.

音場1(図11参照)は、t=0での振幅が1.0の余弦波が第1騒音源28から放射され、t=0での振幅が0.5の余弦波が第2騒音源29から放射されているものとする。 In sound field 1 (see FIG. 11), a cosine wave with an amplitude of 1.0 at t=0 is radiated from first noise source 28, and a cosine wave with an amplitude of 0.5 at t=0 is radiated from second noise source 29.

音場2(図12参照)は、t=0での振幅が1.0の余弦波が第3騒音源30から放射され、t=0での振幅が-0.5の余弦波が第4騒音源31から放射されているものとする。 In sound field 2 (see Figure 12), a cosine wave with an amplitude of 1.0 at t = 0 is radiated from the third noise source 30, and a cosine wave with an amplitude of -0.5 at t = 0 is radiated from the fourth noise source 31.

音場1において、上記した従来の計算手法を用いた解析結果を以下に述べる。
総音圧cの最大振幅は式(4)より1.12である。第1騒音源28を無くした場合の総音圧cの最大振幅は0.5(0.62減)となる。従って、第1騒音源28を無くしたことにより、総音圧cの最大振幅は0.45倍になり、-7dBの減少となった。第2騒音源29を無くした場合の総音圧cの最大振幅は1.0(0.12減)となる。従って、第2騒音源29を無くしたことにより、総音圧cの最大振幅は0.89倍になり、-1dBの減少となった。
The analysis results of the sound field 1 using the above-mentioned conventional calculation method will be described below.
The maximum amplitude of the total sound pressure c is 1.12 according to formula (4). When the first noise source 28 is removed, the maximum amplitude of the total sound pressure c becomes 0.5 (a decrease of 0.62). Therefore, by removing the first noise source 28, the maximum amplitude of the total sound pressure c becomes 0.45 times, a decrease of -7 dB. When the second noise source 29 is removed, the maximum amplitude of the total sound pressure c becomes 1.0 (a decrease of 0.12). Therefore, by removing the second noise source 29, the maximum amplitude of the total sound pressure c becomes 0.89 times, a decrease of -1 dB.

第1騒音源28を無くした場合の総音圧cの最大振幅の変化した値は、0.62(1.12-0.5)の減少である。第2騒音源29を無くした場合の総音圧cの最大振幅の変化した値は、0.12(1.12-1.0)の減少である。従って、耳位置(観測点)における第1騒音源28及び第2騒音源29の寄与度それぞれの算出を以下の様に行ってもよい。
第1騒音源28の寄与度=100×0.62/(0.62+0.12)=83.8%
第2騒音源29の寄与度=100×0.12/(0.62+0.12)=16.2%
The change in the maximum amplitude of the total sound pressure c when the first noise source 28 is removed is a decrease of 0.62 (1.12-0.5). The change in the maximum amplitude of the total sound pressure c when the second noise source 29 is removed is a decrease of 0.12 (1.12-1.0). Therefore, the contributions of the first noise source 28 and the second noise source 29 at the ear position (observation point) may be calculated as follows.
Contribution of first noise source 28=100×0.62/(0.62+0.12)=83.8%
Contribution of second noise source 29=100×0.12/(0.62+0.12)=16.2%

次に、音場1において、本開示の計算手法を用いた解析結果を以下に述べる。
総音圧cの最大振幅は式(5)より1.5である。第1騒音源28を無くした場合の総音圧cの最大振幅は0.5(1.0減)となる。従って、第1騒音源28を無くしたことにより、総音圧cの最大振幅は0.33倍になり、-9.6dBの減少となった。第2騒音源29を無くした場合の総音圧cの最大振幅は1.0(0.5減)となる。従って、第2騒音源29を無くしたことにより、総音圧cの最大振幅は0.7倍になり、-3.1dBの減少となった。
Next, the analysis results of the sound field 1 using the calculation method of the present disclosure will be described below.
The maximum amplitude of the total sound pressure c is 1.5 according to formula (5). When the first noise source 28 is removed, the maximum amplitude of the total sound pressure c is 0.5 (a decrease of 1.0). Therefore, by removing the first noise source 28, the maximum amplitude of the total sound pressure c becomes 0.33 times, a decrease of -9.6 dB. When the second noise source 29 is removed, the maximum amplitude of the total sound pressure c is 1.0 (a decrease of 0.5). Therefore, by removing the second noise source 29, the maximum amplitude of the total sound pressure c becomes 0.7 times, a decrease of -3.1 dB.

第1騒音源28を無くした場合の総音圧cの最大振幅の変化した値は、1.0(1.5--0.5)の減少である。第2騒音源29を無くした場合の総音圧cの最大振幅の変化した値は、0.5(1.5-1.0)の減少である。従って、耳位置(観測点)における第1騒音源28及び第2騒音源29の寄与度それぞれの算出を以下の様に行ってもよい。
第1騒音源28の寄与度=100×1.0/(1.0+0.5)=66.7%
第2騒音源29の寄与度=100×0.5/(1.0+0.5)=33.3%
The change in the maximum amplitude of the total sound pressure c when the first noise source 28 is removed is a decrease of 1.0 (1.5--0.5). The change in the maximum amplitude of the total sound pressure c when the second noise source 29 is removed is a decrease of 0.5 (1.5-1.0). Therefore, the contributions of the first noise source 28 and the second noise source 29 at the ear position (observation point) may be calculated as follows.
Contribution of first noise source 28=100×1.0/(1.0+0.5)=66.7%
Contribution of second noise source 29=100×0.5/(1.0+0.5)=33.3%

音場1において、従来の計算手法の結果と本開示の計算手法の結果ともに、第1騒音源28、第2騒音源29の両方に騒音対策の効果があるという結果となる。さらに、従来の計算手法の結果と本開示の計算手法の結果ともに、第1騒音源28の対策に、より大きな効果があるという示唆が得られる。しかしながらその一方で、以下に示す音場2の場合では、従来の計算手法よりも本開示の計算手法の結果の方が妥当である。 In sound field 1, both the results of the conventional calculation method and the calculation method disclosed herein show that noise control measures are effective for both the first noise source 28 and the second noise source 29. Furthermore, both the results of the conventional calculation method and the calculation method disclosed herein suggest that measures for the first noise source 28 have a greater effect. However, on the other hand, in the case of sound field 2 shown below, the results of the calculation method disclosed herein are more appropriate than those of the conventional calculation method.

次に、音場2において、上記した従来の計算手法を用いた解析結果を以下に述べる。
総音圧cの最大振幅は式(4)より1.12である。第3騒音源30を無くした場合の総音圧cの最大振幅は0.5(0.62減)となる。従って、第3騒音源30を無くしたことにより、総音圧cの最大振幅は0.45倍になり、-7dBの減少となった。第4騒音源31を無くした場合の総音圧cの最大振幅は1.0(0.12減)となる。従って、第4騒音源31を無くしたことにより、総音圧cの最大振幅は0.89倍になり、-1dBの減少となった。
Next, the analysis results for the sound field 2 using the above-mentioned conventional calculation method will be described below.
The maximum amplitude of the total sound pressure c is 1.12 according to formula (4). When the third noise source 30 is removed, the maximum amplitude of the total sound pressure c becomes 0.5 (a decrease of 0.62). Therefore, by removing the third noise source 30, the maximum amplitude of the total sound pressure c becomes 0.45 times, a decrease of -7 dB. When the fourth noise source 31 is removed, the maximum amplitude of the total sound pressure c becomes 1.0 (a decrease of 0.12). Therefore, by removing the fourth noise source 31, the maximum amplitude of the total sound pressure c becomes 0.89 times, a decrease of -1 dB.

第3騒音源30を無くした場合の総音圧cの最大振幅の変化した値は、0.62(1.12-0.5)の減少である。第4騒音源31を無くした場合の総音圧cの最大振幅の変化した値は、0.12(1.12-1.0)の減少である。従って、耳位置(観測点)における第3騒音源30及び第4騒音源31の寄与度それぞれの算出を以下の様に行ってもよい。
第3騒音源30の寄与度=100×0.62/(0.62+0.12)=83.8%
第4騒音源31の寄与度=100×0.12/(0.62+0.12)=16.2%
The change in the maximum amplitude of the total sound pressure c when the third noise source 30 is removed is a decrease of 0.62 (1.12-0.5). The change in the maximum amplitude of the total sound pressure c when the fourth noise source 31 is removed is a decrease of 0.12 (1.12-1.0). Therefore, the contributions of the third noise source 30 and the fourth noise source 31 at the ear position (observation point) may be calculated as follows.
Contribution of third noise source 30=100×0.62/(0.62+0.12)=83.8%
Contribution of fourth noise source 31=100×0.12/(0.62+0.12)=16.2%

次に、音場2において、本開示の計算手法を用いた解析結果を以下に述べる。
総音圧cの最大振幅は式(5)より0.5である。第3騒音源30を無くした場合の総音圧cの最大振幅は0.5(0.0減)となる。従って、第3騒音源30を無くしたことにより、総音圧cの最大振幅は1.0倍になり、効果なし。第4騒音源31を無くした場合の総音圧cの最大振幅は1.0(0.5増)となる。従って、第4騒音源31を無くしたことにより、総音圧cの最大振幅は2倍になり、6.0dBの増加となった。
Next, the analysis results for sound field 2 using the calculation method of the present disclosure will be described below.
The maximum amplitude of the total sound pressure c is 0.5 according to formula (5). When the third noise source 30 is removed, the maximum amplitude of the total sound pressure c is 0.5 (a decrease of 0.0). Therefore, by removing the third noise source 30, the maximum amplitude of the total sound pressure c becomes 1.0 times, and there is no effect. When the fourth noise source 31 is removed, the maximum amplitude of the total sound pressure c becomes 1.0 (an increase of 0.5). Therefore, by removing the fourth noise source 31, the maximum amplitude of the total sound pressure c becomes double, an increase of 6.0 dB.

第3騒音源30を無くした場合の総音圧cの最大振幅の変化した値は、0.0(0.5-0.5)の減少である。第4騒音源31を無くした場合の総音圧cの最大振幅の変化した値は、-0.5(0.5-1.0)の減少、即ち0.5の増加である。従って、耳位置(観測点)における第3騒音源30及び第4騒音源31の寄与度それぞれの算出を以下の様に行ってもよい。
第3騒音源30の寄与度=100×0.0/(0.0+0.5)=0.0%
第4騒音源31の寄与度=100×(-0.5)/(0.0+0.5)=-100.0%
The change in the maximum amplitude of the total sound pressure c when the third noise source 30 is removed is a decrease of 0.0 (0.5-0.5). The change in the maximum amplitude of the total sound pressure c when the fourth noise source 31 is removed is a decrease of -0.5 (0.5-1.0), that is, an increase of 0.5. Therefore, the contributions of the third noise source 30 and the fourth noise source 31 at the ear position (observation point) may be calculated as follows.
Contribution of third noise source 30=100×0.0/(0.0+0.5)=0.0%
Contribution of fourth noise source 31=100×(−0.5)/(0.0+0.5)=−100.0%

実際に音場2では第3騒音源30を対策しても効果が無く、第4騒音源31を対策すると逆に騒音が増加する。本開示の計算手法の結果は当該実際の結果を表しており妥当である。一方で従来の計算手法では第3騒音源30、第4騒音源31ともに対策の効果があり、さらに第3騒音源30の対策に、より大きな効果があるという結果となるため、間違った示唆を与える。 In reality, in sound field 2, taking measures against the third noise source 30 has no effect, and taking measures against the fourth noise source 31 actually increases the noise. The results of the calculation method disclosed herein represent the actual results and are appropriate. On the other hand, the conventional calculation method gives a result in which measures against both the third noise source 30 and the fourth noise source 31 are effective, and furthermore, measures against the third noise source 30 have a greater effect, which gives an erroneous suggestion.

なお、上記した各騒音源の寄与度の算出は、次のように拡張することができる。
すなわち、総音圧cに寄与している騒音源が第1から第Nまで有り、第n(nは1からNまでの任意の整数)騒音源を無くした場合の総音圧cの最大振幅の変化した値をanとする。この場合、第n騒音源の寄与度は式(6)により算出することができる。anは、総音圧cの最大振幅が減少した場合を正とし、総音圧cの最大振幅が増加した場合を負とする。
The calculation of the contribution of each noise source described above can be expanded as follows.
That is, there are first to N noise sources that contribute to the total sound pressure c, and the change in the maximum amplitude of the total sound pressure c when the nth noise source (n is any integer from 1 to N) is removed is denoted as an. In this case, the contribution of the nth noise source can be calculated by equation (6). an is considered positive when the maximum amplitude of the total sound pressure c decreases, and is considered negative when the maximum amplitude of the total sound pressure c increases.

Figure 0007478710000005
Figure 0007478710000005

次に図13を参照して、本実施形態に係る騒音源寄与解析方法について騒音源寄与解析プログラムとともに説明する。図13は、本実施形態に係る騒音源寄与解析プログラムのフローチャートである。 Next, the noise source contribution analysis method according to this embodiment will be described together with the noise source contribution analysis program with reference to FIG. 13. FIG. 13 is a flowchart of the noise source contribution analysis program according to this embodiment.

図13に示す様に、騒音源寄与解析プログラムは、騒音源音圧取得ステップS60、伝達関数取得ステップS61、観測点音圧算出ステップS62、高周波帯域抽出ステップS63、低周波帯域抽出ステップS64、総音圧算出ステップS65、及び寄与度算出ステップS66などを含む。 As shown in FIG. 13, the noise source contribution analysis program includes a noise source sound pressure acquisition step S60, a transfer function acquisition step S61, an observation point sound pressure calculation step S62, a high frequency band extraction step S63, a low frequency band extraction step S64, a total sound pressure calculation step S65, and a contribution degree calculation step S66.

騒音源寄与解析装置50は、ROM50a若しくは記憶部50cに保存された騒音源寄与解析プログラムをメインメモリに取り込み、CPU50dにより騒音源寄与解析プログラムを実行する。 The noise source contribution analysis device 50 loads the noise source contribution analysis program stored in the ROM 50a or the memory unit 50c into the main memory, and executes the noise source contribution analysis program using the CPU 50d.

騒音源寄与解析プログラムは、騒音源寄与解析装置50のCPU50dに対して、騒音源音圧取得機能、伝達関数取得機能、観測点音圧算出機能、高周波帯域抽出機能、低周波帯域抽出機能、総音圧算出機能、及び寄与度算出機能などの機能を実現させる。 The noise source contribution analysis program enables the CPU 50d of the noise source contribution analysis device 50 to realize functions such as a noise source sound pressure acquisition function, a transfer function acquisition function, an observation point sound pressure calculation function, a high frequency band extraction function, a low frequency band extraction function, a total sound pressure calculation function, and a contribution degree calculation function.

これらの機能は図13のフローチャートに示す順序で処理を行う場合を例示したが、これに限らず、これらの順番を適宜入れ替えて騒音源寄与解析プログラムを実行してもよい。 These functions are illustrated as being processed in the order shown in the flowchart of Figure 13, but this is not limiting, and the noise source contribution analysis program may be executed by changing the order as appropriate.

なお、上記した各機能は、前述の騒音源寄与解析装置50の騒音源音圧取得部60、伝達関数取得部61、観測点音圧算出部62、高周波帯域抽出部63、低周波帯域抽出部64、総音圧算出部65、及び寄与度算出部66の説明と重複するため、その詳細な説明は省略する。 The above-mentioned functions overlap with the explanations of the noise source sound pressure acquisition unit 60, transfer function acquisition unit 61, observation point sound pressure calculation unit 62, high frequency band extraction unit 63, low frequency band extraction unit 64, total sound pressure calculation unit 65, and contribution calculation unit 66 of the noise source contribution analysis device 50 described above, so detailed explanations thereof will be omitted.

騒音源音圧取得機能は、複数の騒音源の各々の近傍で検出された騒音源音圧を取得する(S60:騒音源音圧取得ステップ)。 The noise source sound pressure acquisition function acquires the noise source sound pressure detected in the vicinity of each of the multiple noise sources (S60: noise source sound pressure acquisition step).

伝達関数取得機能は、複数の騒音源の各々から観測点までの騒音の伝達特性を示す伝達関数を複数の騒音源ごとに取得する(S61:伝達関数取得ステップ)。 The transfer function acquisition function acquires a transfer function for each of the multiple noise sources that indicates the transfer characteristics of the noise from each of the multiple noise sources to the observation point (S61: transfer function acquisition step).

観測点音圧算出機能は、伝達関数取得部61により取得された複数の騒音源ごとの伝達関数に、当該伝達関数に対応する騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する(S62:観測点音圧算出ステップ)。 The observation point sound pressure calculation function inputs the noise source sound pressure corresponding to each transfer function acquired by the transfer function acquisition unit 61 for each of the multiple noise sources, and outputs the output as the observation point sound pressure (S62: Observation point sound pressure calculation step).

高周波帯域抽出機能は、観測点音圧算出部62により算出された観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する(S63:高周波帯域抽出ステップ)。 The high frequency band extraction function extracts the high frequency band of the observation point sound pressures calculated by the observation point sound pressure calculation unit 62 as high frequency band sound pressures (S63: high frequency band extraction step).

低周波帯域抽出機能は、観測点音圧算出部62により算出された観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する(S64:低周波帯域抽出ステップ)。 The low-frequency band extraction function extracts the low-frequency band from the observation point sound pressures calculated by the observation point sound pressure calculation unit 62 as low-frequency band sound pressures (S64: low-frequency band extraction step).

総音圧算出機能は、複数の騒音源の各々の高周波帯域音圧の2乗和と複数の騒音源の各々の低周波帯域音圧の和とに基づき、観測点における総音圧を算出する(S65:総音圧算出ステップ)。 The total sound pressure calculation function calculates the total sound pressure at the observation point based on the sum of the squares of the high-frequency band sound pressures of each of the multiple noise sources and the sum of the low-frequency band sound pressures of each of the multiple noise sources (S65: total sound pressure calculation step).

寄与度算出機能は、算出された総音圧に対して、複数の騒音源それぞれの寄与の度合を算出する(S66:寄与度算出ステップ)。 The contribution calculation function calculates the degree of contribution of each of the multiple noise sources to the calculated total sound pressure (S66: contribution calculation step).

上述した本実施形態によれば、耳位置(観測点)での総音圧を算出する計算式を高周波帯域と低周波帯域とでそれぞれ異なる計算式を用いた。低周波帯域での総音圧の算出では、騒音の相互間の干渉の影響を考慮した。これにより、耳位置(観測点)の総音圧と計算値との乖離が減少し、約6dB低減し精度を向上させることができる。 According to the present embodiment described above, different formulas are used to calculate the total sound pressure at ear position (observation point) for the high frequency band and the low frequency band. When calculating the total sound pressure in the low frequency band, the effect of interference between noises is taken into account. This reduces the discrepancy between the total sound pressure at ear position (observation point) and the calculated value by approximately 6 dB, improving accuracy.

本開示は上記した実施形態に係る騒音源寄与解析装置50、騒音源寄与解析方法、及び騒音源寄与解析プログラムに限定されるものではなく、特許請求の範囲に記載した本開示の要旨を逸脱しない限りにおいて、その他種々の変形例、若しくは応用例により実施可能である。 The present disclosure is not limited to the noise source contribution analysis device 50, the noise source contribution analysis method, and the noise source contribution analysis program according to the above-described embodiments, and can be implemented by various other modified examples or application examples without departing from the gist of the present disclosure described in the claims.

10 自動車
11 ボディー
12 パネル
13 タイヤ
14 エンジン
15 シート
16 プロペラシャフト
17 フロントガラス
18 リアガラス
20 風の流れ
21 固体伝搬
22 空気伝搬
23 音の放射
24 耳位置(観測点)
28 第1騒音源
29 第2騒音源
30 第3騒音源
31 第4騒音源
32a 第1マイク
32b 第2マイク
32c 第3マイク
32d 第4マイク
33 騒音源
34 第1パネル
35 第2パネル
36 低周波
37 高周波
38 実測値(d)
39 計算値(c)
50 騒音源寄与解析装置
50a Read Only Memory(ROM)
50b Random Access Memory(RAM)
50c 記憶部
50d Central Processing Unit(CPU)
50e 入出力インターフェース
50f 通信インターフェース
50g 入力装置
50h 出力装置
51 通信ネットワーク
52 キーボード
53 マウス
54 Webカメラ
55 モニター
56 スピーカ
57 プリンタ
60 騒音源音圧取得部
61 伝達関数取得部
62 観測点音圧算出部
63 高周波帯域抽出部
63a 高周波帯域抽出部
63b 高周波帯域抽出部
63c 高周波帯域抽出部
63d 高周波帯域抽出部
64 低周波帯域抽出部
64a 低周波帯域抽出部
64b 低周波帯域抽出部
64c 低周波帯域抽出部
64d 低周波帯域抽出部
65 総音圧算出部
65a 高周波帯域加算部
65b 低周波帯域加算部
66 寄与度算出部
67a 伝達関数f1
67b 伝達関数f2
67c 伝達関数f3
67d 伝達関数f4
70 実測値
71 今回計算値
72 従来計算値
75 騒音源
10 Automobile 11 Body 12 Panel 13 Tire 14 Engine 15 Seat 16 Propeller shaft 17 Windshield 18 Rear window 20 Wind flow 21 Solid propagation 22 Air propagation 23 Sound radiation 24 Ear position (observation point)
28 First noise source 29 Second noise source 30 Third noise source 31 Fourth noise source 32a First microphone 32b Second microphone 32c Third microphone 32d Fourth microphone 33 Noise source 34 First panel 35 Second panel 36 Low frequency 37 High frequency 38 Actual measurement value (d)
39 Calculated value (c)
50 Noise source contribution analysis device 50a Read Only Memory (ROM)
50b Random Access Memory (RAM)
50c Memory unit 50d Central Processing Unit (CPU)
50e Input/Output interface 50f Communication interface 50g Input device 50h Output device 51 Communication network 52 Keyboard 53 Mouse 54 Web camera 55 Monitor 56 Speaker 57 Printer 60 Noise source sound pressure acquisition unit 61 Transfer function acquisition unit 62 Observation point sound pressure calculation unit 63 High frequency band extraction unit 63a High frequency band extraction unit 63b High frequency band extraction unit 63c High frequency band extraction unit 63d High frequency band extraction unit 64 Low frequency band extraction unit 64a Low frequency band extraction unit 64b Low frequency band extraction unit 64c Low frequency band extraction unit 64d Low frequency band extraction unit 65 Total sound pressure calculation unit 65a High frequency band addition unit 65b Low frequency band addition unit 66 Contribution degree calculation unit 67a Transfer function f1
67b Transfer function f2
67c Transfer function f3
67d Transfer function f4
70 Actual measurement value 71 Current calculation value 72 Previous calculation value 75 Noise source

Claims (7)

複数の騒音源が寄与する観測点における総音圧に対し、当該騒音源ごとの騒音源音圧の寄与の度合を解析する騒音源寄与解析装置において、
前記複数の騒音源の各々の近傍で検出された前記騒音源音圧を取得する騒音源音圧取得部と、
前記複数の騒音源の各々から前記観測点までの騒音の伝達特性を示す伝達関数を前記複数の騒音源ごとに取得する伝達関数取得部と、
前記伝達関数取得部により取得された前記複数の騒音源ごとの前記伝達関数に、当該伝達関数に対応する前記騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する観測点音圧算出部と、
前記観測点音圧算出部により算出された前記観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する高周波帯域抽出部と、
前記観測点音圧算出部により算出された前記観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する低周波帯域抽出部と、
前記複数の騒音源の各々の高周波帯域音圧の2乗和と前記複数の騒音源の各々の低周波帯域音圧の和とに基づき、前記観測点における前記総音圧を算出する総音圧算出部と、
を備える騒音源寄与解析装置。
A noise source contribution analysis device that analyzes the degree of contribution of a noise source sound pressure of each of a plurality of noise sources to a total sound pressure at an observation point to which the noise sources contribute,
a noise source sound pressure acquisition unit that acquires the noise source sound pressure detected in the vicinity of each of the plurality of noise sources;
a transfer function acquisition unit that acquires, for each of the plurality of noise sources, a transfer function that indicates a transfer characteristic of noise from the each of the plurality of noise sources to the observation point;
an observation point sound pressure calculation unit that inputs the noise source sound pressure corresponding to each of the transfer functions acquired for each of the plurality of noise sources by the transfer function acquisition unit, and outputs the outputs as observation point sound pressures;
a high frequency band extraction unit that extracts, as a high frequency band sound pressure, a high frequency band sound pressure from the observation point sound pressures calculated by the observation point sound pressure calculation unit;
a low-frequency band extraction unit that extracts, as a low-frequency band sound pressure, a low-frequency band sound pressure from the observation point sound pressures calculated by the observation point sound pressure calculation unit;
a total sound pressure calculation unit that calculates the total sound pressure at the observation point based on a sum of squares of high-frequency band sound pressures of each of the plurality of noise sources and a sum of low-frequency band sound pressures of each of the plurality of noise sources;
A noise source contribution analysis device comprising:
算出された前記総音圧に対して、前記複数の騒音源それぞれの寄与の度合を算出する寄与度算出部と、
をさらに備えることを特徴とする請求項1に記載の騒音源寄与解析装置。
a contribution calculation unit that calculates a degree of contribution of each of the plurality of noise sources to the calculated total sound pressure;
The noise source contribution analysis device according to claim 1, further comprising:
前記高周波帯域と前記低周波帯域との境界となるカットオフ周波数は、前記騒音源と前記観測点との間の距離に基づいて設定されることを特徴とする請求項1または2に記載の騒音源寄与解析装置。 The noise source contribution analysis device according to claim 1 or 2, characterized in that the cutoff frequency that is the boundary between the high frequency band and the low frequency band is set based on the distance between the noise source and the observation point. 前記カットオフ周波数は、前記騒音源と前記観測点との間の距離を1波長とする周波数であることを特徴とする請求項3に記載の騒音源寄与解析装置。 The noise source contribution analysis device according to claim 3, characterized in that the cutoff frequency is a frequency in which the distance between the noise source and the observation point is one wavelength. 前記高周波帯域と前記低周波帯域との境界となるカットオフ周波数は、前記騒音源が存在する面の寸法に基づいて設定されることを特徴とする請求項1または2に記載の騒音源寄与解析装置。 The noise source contribution analysis device according to claim 1 or 2, characterized in that the cutoff frequency that is the boundary between the high frequency band and the low frequency band is set based on the dimensions of the surface on which the noise source exists. 複数の騒音源が寄与する観測点における総音圧に対し、当該騒音源ごとの騒音源音圧の寄与の度合を解析する騒音源寄与解析方法において、
騒音源寄与解析装置に用いられるコンピュータに、
前記複数の騒音源の各々の近傍で検出された前記騒音源音圧を取得する騒音源音圧取得ステップと、
前記複数の騒音源の各々から前記観測点までの騒音の伝達特性を示す伝達関数を前記複数の騒音源ごとに取得する伝達関数取得ステップと、
前記伝達関数取得ステップにおいて取得された前記複数の騒音源ごとの前記伝達関数に、当該伝達関数に対応する前記騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する観測点音圧算出ステップと、
前記観測点音圧算出ステップにおいて算出された前記観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する高周波帯域抽出ステップと、
前記観測点音圧算出ステップにおいて算出された前記観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する低周波帯域抽出ステップと、
前記複数の騒音源の各々の高周波帯域音圧の2乗和と前記複数の騒音源の各々の低周波帯域音圧の和とに基づき、前記観測点における前記総音圧を算出する総音圧算出ステップと、
を実行させる騒音源寄与解析方法。
A noise source contribution analysis method for analyzing the degree of contribution of each noise source to a total sound pressure at an observation point to which a plurality of noise sources contribute, comprising:
The computer used in the noise source contribution analysis device is
a noise source sound pressure acquisition step of acquiring the noise source sound pressure detected in the vicinity of each of the plurality of noise sources;
a transfer function acquisition step of acquiring, for each of the plurality of noise sources, a transfer function indicating a transfer characteristic of noise from each of the plurality of noise sources to the observation point;
an observation point sound pressure calculation step of inputting the noise source sound pressure corresponding to each of the transfer functions acquired for the plurality of noise sources in the transfer function acquisition step, and outputting the outputs as observation point sound pressures;
a high frequency band extraction step of extracting, as a high frequency band sound pressure, a high frequency band sound pressure from the observation point sound pressures calculated in the observation point sound pressure calculation step;
a low-frequency band extraction step of extracting, as a low-frequency band sound pressure, a low-frequency band sound pressure from the observation point sound pressures calculated in the observation point sound pressure calculation step;
a total sound pressure calculation step of calculating the total sound pressure at the observation point based on a sum of squares of high-frequency band sound pressures of each of the plurality of noise sources and a sum of low-frequency band sound pressures of each of the plurality of noise sources;
A noise source contribution analysis method that performs the above.
複数の騒音源が寄与する観測点における総音圧に対し、当該騒音源ごとの騒音源音圧の寄与の度合を解析する騒音源寄与解析プログラムにおいて、
騒音源寄与解析装置に用いられるコンピュータに、
前記複数の騒音源の各々の近傍で検出された前記騒音源音圧を取得する騒音源音圧取得機能と、
前記複数の騒音源の各々から前記観測点までの騒音の伝達特性を示す伝達関数を前記複数の騒音源ごとに取得する伝達関数取得機能と、
前記伝達関数取得機能により取得された前記複数の騒音源ごとの前記伝達関数に、当該伝達関数に対応する前記騒音源音圧をそれぞれ入力し、その出力をそれぞれ観測点音圧として出力する観測点音圧算出機能と、
前記観測点音圧算出機能により算出された前記観測点音圧のうち高周波帯域のものを高周波帯域音圧として抽出する高周波帯域抽出機能と、
前記観測点音圧算出機能により算出された前記観測点音圧のうち低周波帯域のものを低周波帯域音圧として抽出する低周波帯域抽出機能と、
前記複数の騒音源の各々の高周波帯域音圧の2乗和と前記複数の騒音源の各々の低周波帯域音圧の和とに基づき、前記観測点における前記総音圧を算出する総音圧算出機能と、
を発揮させる騒音源寄与解析プログラム。
In a noise source contribution analysis program that analyzes the degree of contribution of each noise source to the total sound pressure at an observation point to which multiple noise sources contribute,
The computer used in the noise source contribution analysis device is
a noise source sound pressure acquisition function for acquiring the noise source sound pressure detected in the vicinity of each of the plurality of noise sources;
a transfer function acquisition function that acquires, for each of the plurality of noise sources, a transfer function that indicates a transfer characteristic of noise from each of the plurality of noise sources to the observation point;
an observation point sound pressure calculation function that inputs the noise source sound pressure corresponding to each of the transfer functions acquired for each of the plurality of noise sources by the transfer function acquisition function, and outputs the outputs as observation point sound pressures;
a high frequency band extraction function for extracting, as a high frequency band sound pressure, a high frequency band sound pressure from the observation point sound pressures calculated by the observation point sound pressure calculation function;
a low-frequency band extraction function for extracting, as a low-frequency band sound pressure, a low-frequency band sound pressure from the observation point sound pressures calculated by the observation point sound pressure calculation function;
a total sound pressure calculation function that calculates the total sound pressure at the observation point based on the sum of squares of the high-frequency band sound pressures of each of the plurality of noise sources and the sum of the low-frequency band sound pressures of each of the plurality of noise sources;
A noise source contribution analysis program that demonstrates
JP2021130060A 2021-08-06 2021-08-06 Noise source contribution analysis device, noise source contribution analysis method, and noise source contribution analysis program Active JP7478710B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021130060A JP7478710B2 (en) 2021-08-06 2021-08-06 Noise source contribution analysis device, noise source contribution analysis method, and noise source contribution analysis program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021130060A JP7478710B2 (en) 2021-08-06 2021-08-06 Noise source contribution analysis device, noise source contribution analysis method, and noise source contribution analysis program

Publications (2)

Publication Number Publication Date
JP2023024023A JP2023024023A (en) 2023-02-16
JP7478710B2 true JP7478710B2 (en) 2024-05-07

Family

ID=85204288

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021130060A Active JP7478710B2 (en) 2021-08-06 2021-08-06 Noise source contribution analysis device, noise source contribution analysis method, and noise source contribution analysis program

Country Status (1)

Country Link
JP (1) JP7478710B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116025544B (en) * 2023-03-27 2023-06-30 江西五十铃汽车有限公司 Method, system and equipment for reducing working noise of in-vehicle compressor

Also Published As

Publication number Publication date
JP2023024023A (en) 2023-02-16

Similar Documents

Publication Publication Date Title
JP6325663B2 (en) Method for determining noise sound contribution of noise source of motor driven moving body
US11361782B2 (en) Abnormal noise determination apparatus and method
JP7478710B2 (en) Noise source contribution analysis device, noise source contribution analysis method, and noise source contribution analysis program
JP2006523828A (en) Method and device for determining acoustic transfer impedance
JP6754738B2 (en) Spatial acoustic analysis method and its system
US6324290B1 (en) Method and apparatus for diagnosing sound source and sound vibration source
Moron et al. A CFD/SEA approach for prediction of vehicle interior noise due to wind noise
JP2013079953A (en) Method for analyzing sound transmission paths in system
CN106258000A (en) Acceleration detector and active noise controller
Wang et al. CSA-based acoustic beamforming for the contribution analysis of air-borne tyre noise
JP6433630B2 (en) Noise removing device, echo canceling device, abnormal sound detecting device, and noise removing method
JP2009073417A (en) Apparatus and method for controlling noise
JP2024502827A (en) Apparatus, system, and method for testing acoustic equipment
Bremner et al. A transient sea model for transmission of non-stationary wind noise
JP7196504B2 (en) Acoustic property measurement system, acoustic property measurement method, and program
JP2021152499A (en) Abnormal noise generating device and abnormal noise generation method
JP2017041752A (en) Acoustic device for vehicle
CN112334971A (en) Concurrent FXLMS system using common reference signal and error signal
Fonseca et al. Pass-by noise measurements using beamforming technique
Sebald Quantifying Door Closing Sound Quality, Measurement Techniques and Influence of Selected Parameters
KR20230092421A (en) Method And Apparatus for Noise Canceling Based on Partial Correlation Function
Pierro et al. Using pu probes for the experimental vibro-acoustical modal analysis of a helicopter
JPH0758223B2 (en) Method of measuring sound deadening / damping effect, measuring device, and signal source search device
Padois et al. Acoustic imaging using different weighting functions with the generalized cross correlation based on the generalized mean
Oettle et al. Beamforming Quantification of Acoustic Transmission Paths for Passenger Vehicles Using a Reciprocal Approach

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230810

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240315

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240422

R150 Certificate of patent or registration of utility model

Ref document number: 7478710

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150