JP6977817B1 - 自動車のパネル部品の振動騒音低減解析装置及び解析方法 - Google Patents
自動車のパネル部品の振動騒音低減解析装置及び解析方法 Download PDFInfo
- Publication number
- JP6977817B1 JP6977817B1 JP2020113763A JP2020113763A JP6977817B1 JP 6977817 B1 JP6977817 B1 JP 6977817B1 JP 2020113763 A JP2020113763 A JP 2020113763A JP 2020113763 A JP2020113763 A JP 2020113763A JP 6977817 B1 JP6977817 B1 JP 6977817B1
- Authority
- JP
- Japan
- Prior art keywords
- vibration
- model
- vibration noise
- weight
- noise reduction
- 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
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B62—LAND VEHICLES FOR TRAVELLING OTHERWISE THAN ON RAILS
- B62D—MOTOR VEHICLES; TRAILERS
- B62D65/00—Designing, manufacturing, e.g. assembling, facilitating disassembly, or structurally modifying motor vehicles or trailers, not otherwise provided for
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B62—LAND VEHICLES FOR TRAVELLING OTHERWISE THAN ON RAILS
- B62D—MOTOR VEHICLES; TRAILERS
- B62D25/00—Superstructure or monocoque structure sub-units; Parts or details thereof not otherwise provided for
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Manufacturing & Machinery (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Body Structure For Vehicles (AREA)
Abstract
【課題】自動車のパネル部品の振動騒音を低減させるために車体を軽量化する部位を特定する自動車のパネル部品の振動騒音低減解析装置及び解析方法を提供する。【解決手段】本発明は、車体モデル２００を取得し（１３）、車体モデル２００に複数の区分領域を設定し（１５）、振動騒音低減対象パネル部品モデルを設定し（１７）、車体モデル２００の振動解析により振動騒音低減対象パネル部品モデルの振動モードと等価放射パワーのピーク周波数を選定し（１９）、一つ又は複数の区分領域の重量を変更して振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得し（２１）、取得した等価放射パワーのピーク周波数に対する各区分領域の重量の寄与度を算出し（Ｓ２３）、各区分領域の重量の寄与度に基づいて、自動車のパネル部品からの振動騒音を低減するために車体を軽量化する部位を特定する（２５）、ものである。【選択図】 図１
Description
本発明は、自動車のパネル部品の振動騒音低減解析装置及び解析方法に関し、特に、自動車の起振源からの振動に起因するパネル部品の振動騒音を低減させ、かつ車体を軽量化できる部位を特定する自動車のパネル部品の振動騒音低減解析装置及び解析方法に関する。
フロアパネル、ダッシュロアパネル及びルーフパネル等の自動車のパネル部品は、鋼板、アルミニウム合金板、またはこれらにZn合金等をめっきしたものをプレス成形してなる。これらパネル部品の振動は、ロードノイズやこもり音の原因となり、車室内騒音を悪化させる。そのため、車室内騒音の低減は、自動車の静粛性や商品価値を向上させる上で課題となっている。
図２に例示するような、車体骨格を構成する骨格部品１０１とパネル部品１０３とを備えてなる自動車１００におけるパネル部品１０３の振動騒音は、（ａ）自動車１００のエンジン１０５ａや走行時に路面等からタイヤ１０５ｂを介して入力する振動が、（ｂ）骨格部品１０１を伝達し、（ｃ）パネル部品１０３を振動させることが原因である。
従来、（ｃ）パネル部品１０３の振動騒音を低減する技術としては、当該パネル部品にビードを付与することが有効であるとされている。そして、非特許文献１には、振動騒音の低減対象とするパネル部品に付与するビードの最適な位置を求める技術がトポグラフィー最適化として開示されている。
しかしながら、意匠性が重要視される自動車においては、ルーフパネルのような外板パネル部品にビードを直接付与することは困難であり、また、隣接する内板パネル部品に干渉するため問題であった。その結果、自動車の振動伝達経路にある（ｂ）骨格部品１０１に対策を施すことが求められていた。
しかしながら、意匠性が重要視される自動車においては、ルーフパネルのような外板パネル部品にビードを直接付与することは困難であり、また、隣接する内板パネル部品に干渉するため問題であった。その結果、自動車の振動伝達経路にある（ｂ）骨格部品１０１に対策を施すことが求められていた。
そこで、自動車の振動伝達経路の骨格部品に対策を施してパネル部品の振動を低減する技術に関しても、これまでにいくつか提案されている。
例えば、特許文献１には、車体前後方向に所定の間隔で複数設けたルーフボウの上にルーフパネルが接合される車体のルーフ構造において、ある一つのルーフボウと該ルーフボウに隣接する少なくとも１つのルーフボウの車体幅方向における中央部の幅を、他のルーフボウの幅よりも幅広に形成することで、エンジン振動等によって車体骨格に共振が生じた場合でもルーフパネルの振動を抑制する技術が開示されている。
例えば、特許文献１には、車体前後方向に所定の間隔で複数設けたルーフボウの上にルーフパネルが接合される車体のルーフ構造において、ある一つのルーフボウと該ルーフボウに隣接する少なくとも１つのルーフボウの車体幅方向における中央部の幅を、他のルーフボウの幅よりも幅広に形成することで、エンジン振動等によって車体骨格に共振が生じた場合でもルーフパネルの振動を抑制する技術が開示されている。
"構造最適化設計ソフトウェア Altair OptiStruct"、［online］、［令和１年７月１８日検索］、インターネット＜URL：https://www.terrabyte.co.jp/Hyper/OptiStruct-3.htm＞
しかしながら、特許文献１に開示されている技術は、ルーフボウの中央部を幅広に形成することで重量が大幅に増加してしまうという問題があり、また、骨格部品の形状を大きく変えると他の部品と干渉してしまって実現するのが難しい場合があった。
そのため、部品形状を大きく変えずに車体を軽量化しつつ、パネル部品の振動騒音を効率的に低減できる技術が望まれていた。
そのため、部品形状を大きく変えずに車体を軽量化しつつ、パネル部品の振動騒音を効率的に低減できる技術が望まれていた。
本発明は、上記のような課題を解決するためになされたものであり、自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させ、かつ前記自動車の車体の軽量化ができる部位を特定する自動車のパネル部品の振動騒音低減解析装置及び解析方法を提案することを目的とする。
（１）本発明に係る自動車のパネル部品の振動騒音低減解析装置は、自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させ、かつ前記自動車の車体の軽量化ができる部位を特定するものであって、
前記自動車の車体を構成する骨格部品とパネル部品のそれぞれをメッシュでモデル化した骨格部品モデルとパネル部品モデルを備えてなり、前記起振源からの振動を入力する加振部が設定された車体モデルを取得する車体モデル取得部と、
該車体モデルに、前記骨格部品モデル及び前記パネル部品モデルに基づいて区分した複数の区分領域を設定する区分領域設定部と、
操作者からの指示により、前記パネル部品モデルのうち、振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして前記車体モデルに設定する振動騒音低減対象パネル部品モデル設定部と、
前記車体モデルを用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定する振動モード・等価放射パワーピーク周波数選定部と、
前記車体モデルにおける一つ又は複数の区分領域の重量を変更し、前記車体モデルにおける前記区分領域の重量の組み合わせごとに振動解析を行い、前記振動モード・等価放射パワーピーク周波数選定部による処理を行って選定した振動モードにおける前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得する区分領域重量変更ピーク周波数取得部と、
前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、前記区分領域の重量を説明変数とする多変量解析を行い、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する前記区分領域の重量の寄与度を算出する区分領域重量寄与度算出部と、
該算出した前記各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために、前記自動車の車体における軽量化する部位を特定する振動騒音低減軽量化部位特定部と、を備えたことを特徴とするものである。
前記自動車の車体を構成する骨格部品とパネル部品のそれぞれをメッシュでモデル化した骨格部品モデルとパネル部品モデルを備えてなり、前記起振源からの振動を入力する加振部が設定された車体モデルを取得する車体モデル取得部と、
該車体モデルに、前記骨格部品モデル及び前記パネル部品モデルに基づいて区分した複数の区分領域を設定する区分領域設定部と、
操作者からの指示により、前記パネル部品モデルのうち、振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして前記車体モデルに設定する振動騒音低減対象パネル部品モデル設定部と、
前記車体モデルを用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定する振動モード・等価放射パワーピーク周波数選定部と、
前記車体モデルにおける一つ又は複数の区分領域の重量を変更し、前記車体モデルにおける前記区分領域の重量の組み合わせごとに振動解析を行い、前記振動モード・等価放射パワーピーク周波数選定部による処理を行って選定した振動モードにおける前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得する区分領域重量変更ピーク周波数取得部と、
前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、前記区分領域の重量を説明変数とする多変量解析を行い、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する前記区分領域の重量の寄与度を算出する区分領域重量寄与度算出部と、
該算出した前記各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために、前記自動車の車体における軽量化する部位を特定する振動騒音低減軽量化部位特定部と、を備えたことを特徴とするものである。
（２）本発明に係る自動車のパネル部品の振動騒音低減解析方法は、自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させ、かつ前記自動車の車体の軽量化ができる部位を特定するために、コンピュータが以下の各ステップを行うものであって、
前記自動車の車体を構成する骨格部品とパネル部品のそれぞれをメッシュでモデル化した骨格部品モデルとパネル部品モデルを備えてなり、前記起振源からの振動を入力する加振部が設定された車体モデルを取得する車体モデル取得ステップと、
該車体モデルに、前記骨格部品モデル及び前記パネル部品モデルに基づいて区分した複数の区分領域を設定する区分領域設定ステップと、
前記パネル部品モデルのうち、振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして前記車体モデルに設定する振動騒音低減対象パネル部品モデル設定ステップと、
前記車体モデルを用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定する振動モード・等価放射パワーピーク周波数選定ステップと、
前記車体モデルにおける一つ又は複数の区分領域の重量を変更し、前記車体モデルにおける前記区分領域の重量の組み合わせごとに振動解析を行い、前記振動モード・等価放射パワーピーク周波数選定ステップによる処理を行って選定した振動モードにおける前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得する区分領域重量変更ピーク周波数取得ステップと、
前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、前記区分領域の重量を説明変数とする多変量解析を行い、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する前記区分領域の重量の寄与度を算出する区分領域重量寄与度算出ステップと、
該算出した前記各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために、前記自動車の車体における軽量化する部位を特定する振動騒音低減軽量化部位特定ステップと、を備えたことを特徴とするものである。
前記自動車の車体を構成する骨格部品とパネル部品のそれぞれをメッシュでモデル化した骨格部品モデルとパネル部品モデルを備えてなり、前記起振源からの振動を入力する加振部が設定された車体モデルを取得する車体モデル取得ステップと、
該車体モデルに、前記骨格部品モデル及び前記パネル部品モデルに基づいて区分した複数の区分領域を設定する区分領域設定ステップと、
前記パネル部品モデルのうち、振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして前記車体モデルに設定する振動騒音低減対象パネル部品モデル設定ステップと、
前記車体モデルを用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定する振動モード・等価放射パワーピーク周波数選定ステップと、
前記車体モデルにおける一つ又は複数の区分領域の重量を変更し、前記車体モデルにおける前記区分領域の重量の組み合わせごとに振動解析を行い、前記振動モード・等価放射パワーピーク周波数選定ステップによる処理を行って選定した振動モードにおける前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得する区分領域重量変更ピーク周波数取得ステップと、
前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、前記区分領域の重量を説明変数とする多変量解析を行い、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する前記区分領域の重量の寄与度を算出する区分領域重量寄与度算出ステップと、
該算出した前記各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために、前記自動車の車体における軽量化する部位を特定する振動騒音低減軽量化部位特定ステップと、を備えたことを特徴とするものである。
（３）本発明に係る自動車のパネル部品の振動騒音低減解析方法は、上記（２）に記載の自動車のパネル部品の振動騒音低減解析方法の前記区分領域重量寄与度算出ステップにおいて、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出した関係式に少なくとも一つ以上の前記各区分領域の重量を変更して与えることにより、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を予測し、自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させるために前記自動車の車体における軽量化する部位を特定することを特徴とするものである。
（４）本発明に係る自動車のパネル部品の振動騒音低減解析方法は、上記（２）に記載の自動車のパネル部品の振動騒音低減解析方法の前記区分領域重量寄与度算出ステップにおいて、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出した関係式に板厚を算出したい一つの区分領域を除いた他の区分領域の重量と前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数とを与えて、前記選択した一つの区分領域の重量を算出し、該算出した前記一つの区分領域の重量に基づいて、該一つの区分領域を構成する前記骨格部品モデル及び／又はパネル部品モデルの板厚を算出することを特徴とするものである。
（５）上記（４）に記載のものにおいて、
前記関係式により前記選択した一つの区分領域の重量を算出するに際し、該一つの区分領域の重量の変化が前記等価放射パワーのピーク周波数に及ぼす影響を補正する補正値を前記関係式に与え、前記一つの区分領域の重量を算出することを特徴とするものである。
前記関係式により前記選択した一つの区分領域の重量を算出するに際し、該一つの区分領域の重量の変化が前記等価放射パワーのピーク周波数に及ぼす影響を補正する補正値を前記関係式に与え、前記一つの区分領域の重量を算出することを特徴とするものである。
本発明においては、振動騒音の低減対象とするパネル部品に対応するパネル部品モデルの等価放射パワーのピーク周波数に対する車体モデルに設定した区分領域の重量の寄与度を算出し、該算出した寄与度に基づいて、前記自動車の車体において軽量化することでパネル部品からの振動騒音の低減に寄与する部位を特定することができ、パネル部品からの振動騒音の低減と車体の軽量化とを両立させることができる。
本発明の実施の形態に係る自動車のパネル部品の振動騒音低減解析装置及び解析方法について説明するに先立ち、本発明で解析対象とする自動車について説明する。
＜対象とする自動車＞
本発明で解析対象とする自動車１００は、図２に示すように、骨格部品１０１及びパネル部品１０３と、自動車１００の車体に振動（周期的な荷重）を与える起振源１０５を有する。
本発明で解析対象とする自動車１００は、図２に示すように、骨格部品１０１及びパネル部品１０３と、自動車１００の車体に振動（周期的な荷重）を与える起振源１０５を有する。
骨格部品１０１は、自動車の車体骨格を構成する部品であり、ルーフレール、Ａピラ、Ｂピラ、Ｃピラ、サイドシルなどが例示される。
パネル部品１０３は、薄板構造の部品である外板パネルや内板パネルであり、ルーフパネルやフロアパネル（例えば、ミドルフロア、リアフロア）等が例示される。
起振源１０５は、エンジン１０５ａやタイヤ１０５ｂ等が例示され、自動車の走行時においては、路面からの振動がタイヤ１０５ｂを介してフロントサスペンション取付部１０６やリアサスペンション取付部１０７等の加振部に入力される。
パネル部品１０３は、薄板構造の部品である外板パネルや内板パネルであり、ルーフパネルやフロアパネル（例えば、ミドルフロア、リアフロア）等が例示される。
起振源１０５は、エンジン１０５ａやタイヤ１０５ｂ等が例示され、自動車の走行時においては、路面からの振動がタイヤ１０５ｂを介してフロントサスペンション取付部１０６やリアサスペンション取付部１０７等の加振部に入力される。
以下に説明する本発明の実施の形態１〜実施の形態３は、自動車１００のミドルフロア１０３ａを振動騒音の低減対象とするパネル部品とする例についてのものである。もっとも、本発明は、他のパネル部品（ルーフやリアフロア）を振動騒音の低減対象とするものであってもよい。
［実施の形態１］
＜自動車のパネル部品の振動騒音低減解析装置＞
本発明の実施の形態１に係る自動車のパネル部品の振動騒音低減解析装置（以下、単に「振動騒音低減解析装置」という。）の構成について、以下に説明する。
＜自動車のパネル部品の振動騒音低減解析装置＞
本発明の実施の形態１に係る自動車のパネル部品の振動騒音低減解析装置（以下、単に「振動騒音低減解析装置」という。）の構成について、以下に説明する。
本実施の形態に係る振動騒音低減解析装置１は、自動車１００（図２）における起振源１０５からの振動に起因するパネル部品１０３からの振動騒音を低減させ、かつ自動車１００の車体の軽量化ができる部位を特定するものであって、図１に例示するように、PC（パーソナルコンピュータ）等によって構成され、表示装置３、入力装置５、記憶装置７、作業用データメモリ９及び演算処理部１１を有している。
そして、表示装置３、入力装置５、記憶装置７及び作業用データメモリ９は、演算処理部１１に接続され、演算処理部１１からの指令によってそれぞれの機能が実行される。
以下、自動車１００のパネル部品１０３の一つであるミドルフロア１０３ａを振動騒音の低減対象とする場合を例として、振動騒音低減解析装置１の各構成について説明する。
以下、自動車１００のパネル部品１０３の一つであるミドルフロア１０３ａを振動騒音の低減対象とする場合を例として、振動騒音低減解析装置１の各構成について説明する。
≪表示装置≫
表示装置３は、解析結果の表示等に用いられ、液晶モニター等で構成される。
表示装置３は、解析結果の表示等に用いられ、液晶モニター等で構成される。
≪入力装置≫
入力装置５は、車体モデルファイル３０の表示指示や操作者の条件入力等に用いられ、キーボードやマウス等で構成される。
入力装置５は、車体モデルファイル３０の表示指示や操作者の条件入力等に用いられ、キーボードやマウス等で構成される。
≪記憶装置≫
記憶装置７は、車体モデルファイル３０等の各種ファイルの記憶等に用いられ、ハードディスク等で構成される。
記憶装置７は、車体モデルファイル３０等の各種ファイルの記憶等に用いられ、ハードディスク等で構成される。
車体モデル２００（図３参照）とは、車体骨格構造を構成する複数の骨格部品をメッシュ（平面要素及び／又は立体要素）でモデル化した複数の骨格部品モデルと、パネル部品をメッシュ（平面要素）でモデル化したパネル部品モデルと、を備えてなり、自動車の起振源からの振動が入力する加振部が設定されたものである。
そして、車体モデルファイル３０は、車体モデル２００の各種情報が格納されたものであり、各種情報として、骨格部品モデル及びパネル部品モデルの要素及び節点に関する情報や材料特性に関する情報等が格納されている。
そして、車体モデルファイル３０は、車体モデル２００の各種情報が格納されたものであり、各種情報として、骨格部品モデル及びパネル部品モデルの要素及び節点に関する情報や材料特性に関する情報等が格納されている。
≪作業用データメモリ≫
作業用データメモリ９は、演算処理部１１で使用するデータの一時保存（記憶）や演算に用いられ、RAM（Random Access Memory）等で構成される。
作業用データメモリ９は、演算処理部１１で使用するデータの一時保存（記憶）や演算に用いられ、RAM（Random Access Memory）等で構成される。
≪演算処理部≫
演算処理部１１は、図１に示すように、車体モデル取得部１３と、区分領域設定部１５と、振動騒音低減対象パネル部品モデル設定部１７と、振動モード・等価放射パワーピーク周波数選定部１９と、区分領域重量変更ピーク周波数取得部２１と、区分領域重量寄与度算出部２３と、振動騒音低減軽量化部位特定部２５と、を備え、PC等のCPU（中央演算処理装置）によって構成される。これらの各部は、CPUが所定のプログラムを実行することによって機能する。
演算処理部１１における上記の各部の機能を以下に説明する。
演算処理部１１は、図１に示すように、車体モデル取得部１３と、区分領域設定部１５と、振動騒音低減対象パネル部品モデル設定部１７と、振動モード・等価放射パワーピーク周波数選定部１９と、区分領域重量変更ピーク周波数取得部２１と、区分領域重量寄与度算出部２３と、振動騒音低減軽量化部位特定部２５と、を備え、PC等のCPU（中央演算処理装置）によって構成される。これらの各部は、CPUが所定のプログラムを実行することによって機能する。
演算処理部１１における上記の各部の機能を以下に説明する。
（車体モデル取得部）
車体モデル取得部１３は、自動車１００の骨格部品１０１とパネル部品１０３（図２）のそれぞれをメッシュ（平面要素及び／又は立体要素）でモデル化した骨格部品モデルとパネル部品モデルを備えるとともに、加振部が設定された車体モデル２００を取得するものである。
車体モデル取得部１３は、自動車１００の骨格部品１０１とパネル部品１０３（図２）のそれぞれをメッシュ（平面要素及び／又は立体要素）でモデル化した骨格部品モデルとパネル部品モデルを備えるとともに、加振部が設定された車体モデル２００を取得するものである。
図３〜図５に、骨格部品モデルとパネル部品モデルとを備えてなる車体モデル２００の一例を示す。
骨格部品モデルの例として、ルーフレールモデル２０１、Ａピラモデル２０３、Ｂピラモデル２０５、Ｃピラモデル２０７、ミドルフロアメンバモデル２０９、リアフロアサイドメンバモデル２１１、リアフロアクロスメンバモデル２１３を備え、これらは、平面要素及び／又は立体要素でモデル化されている。
パネル部品モデルの例として、本実施の形態において振動騒音の低減対象とするミドルフロアに対応するミドルフロアモデル２２１の他、ルーフモデル２２３、リアフロアモデル２２５を備え、これらは平面要素でモデル化されている。
加振部の例として、自動車のタイヤからの振動が入力するリアサスペンション取付部１０７に対応する部位２３１を備える。
加振部の例として、自動車のタイヤからの振動が入力するリアサスペンション取付部１０７に対応する部位２３１を備える。
車体モデル２００は、記憶装置７に記憶された車体モデルファイル３０から要素情報や材料特性情報を読み込むことにより取得することができる。
なお、車体モデル２００における加振部２３１は、予め車体モデル２００に設定された情報が車体モデルファイル３０に格納されたものであってもよいし、車体モデルファイル３０から車体モデル２００を読み込んだ後、操作者の指示により車体モデル２００に設定するものであってもよい。
なお、車体モデル２００における加振部２３１は、予め車体モデル２００に設定された情報が車体モデルファイル３０に格納されたものであってもよいし、車体モデルファイル３０から車体モデル２００を読み込んだ後、操作者の指示により車体モデル２００に設定するものであってもよい。
（区分領域設定部）
区分領域設定部１５は、骨格部品モデル及びパネル部品モデルに基づいて区分した複数の区分領域を車体モデル２００に設定するものである。
区分領域設定部１５は、骨格部品モデル及びパネル部品モデルに基づいて区分した複数の区分領域を車体モデル２００に設定するものである。
図６及び図７に、車体モデル２００に１２個の区分領域D1〜D12を設定した例を示す。
D1.ルーフは、ルーフモデル２２３により区分した区分領域である。
D2.ルーフレールは、ルーフレールモデル２０１により区分した区分領域である。
D3.リアクウォータサイドは、ホイールハウスモデル２４７と、クウォータインナモデル２４５により区分した区分領域である。
D4.リアクウォータアッパは、アッパバックパネルモデル２４１により区分した区分領域である。
D5.トランクルーム開口部は、ロアバックパネルモデル２４３により区分した区分領域である。
D6.Ｂピラは、Ｂピラモデル２０５により区分した区分領域である。
D1.ルーフは、ルーフモデル２２３により区分した区分領域である。
D2.ルーフレールは、ルーフレールモデル２０１により区分した区分領域である。
D3.リアクウォータサイドは、ホイールハウスモデル２４７と、クウォータインナモデル２４５により区分した区分領域である。
D4.リアクウォータアッパは、アッパバックパネルモデル２４１により区分した区分領域である。
D5.トランクルーム開口部は、ロアバックパネルモデル２４３により区分した区分領域である。
D6.Ｂピラは、Ｂピラモデル２０５により区分した区分領域である。
D7.リアフロアは、リアフロアモデル２２５により区分した区分領域である。
D8.ミドルフロアは、ミドルフロアモデル２２１により区分した区分領域である。
D9.Ａピラは、Ａピラモデル２０３により区分した区分領域である。
D10.ミドルフロアメンバは、ミドルフロアメンバモデル２０９により区分した区分領域である。
D11.リアフロアサイドメンバは、リアフロアサイドメンバモデル２１１により区分した区分領域である。
D12.リアフロアクロスメンバは、リアフロアクロスメンバモデル２１３により区分した区分領域である。
D8.ミドルフロアは、ミドルフロアモデル２２１により区分した区分領域である。
D9.Ａピラは、Ａピラモデル２０３により区分した区分領域である。
D10.ミドルフロアメンバは、ミドルフロアメンバモデル２０９により区分した区分領域である。
D11.リアフロアサイドメンバは、リアフロアサイドメンバモデル２１１により区分した区分領域である。
D12.リアフロアクロスメンバは、リアフロアクロスメンバモデル２１３により区分した区分領域である。
D1.ルーフ、D7.リアフロア、D8.ミドルフロアは、いずれも、パネル部品モデルのみ区分された区分領域である。
D2.ルーフレール、D6.Ｂピラ、D9.Ａピラ、D10.ミドルフロアメンバ、D11.リアフロアサイドメンバ、D12.リアフロアクロスメンバは、いずれも、骨格部品モデルのみで区分された区分領域である。
D3.リアクウォータサイド、D4.リアクウォータアッパ、D5.トランクルーム開口部、いずれも、パネル部品モデルと骨格部品モデルとを含む区分領域である。
D2.ルーフレール、D6.Ｂピラ、D9.Ａピラ、D10.ミドルフロアメンバ、D11.リアフロアサイドメンバ、D12.リアフロアクロスメンバは、いずれも、骨格部品モデルのみで区分された区分領域である。
D3.リアクウォータサイド、D4.リアクウォータアッパ、D5.トランクルーム開口部、いずれも、パネル部品モデルと骨格部品モデルとを含む区分領域である。
区分領域設定部１５により車体モデル２００に設定する区分領域は、上記のように車体を製造する工程を考慮したASSY単位や、骨格部品モデル単位やパネル部品モデル単位としてもよい。あるいは、感度解析等のCAE解析手法を用いて車体モデル２００に区分領域を設定するものであってもよい。
（振動騒音低減対象設定部）
振動騒音低減対象パネル部品モデル設定部１７は、操作者からの指示により、パネル部品モデルのうち振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして車体モデル２００に設定するものである。本実施の形態では、自動車１００のミドルフロア１０３ａ（図２）を振動騒音の低減対象とし、車体モデル２００におけるパネル部品モデルのうち、ミドルフロアモデル２２１（図５）を振動騒音低減対象パネル部品モデルとして設定した。
振動騒音低減対象パネル部品モデル設定部１７は、操作者からの指示により、パネル部品モデルのうち振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして車体モデル２００に設定するものである。本実施の形態では、自動車１００のミドルフロア１０３ａ（図２）を振動騒音の低減対象とし、車体モデル２００におけるパネル部品モデルのうち、ミドルフロアモデル２２１（図５）を振動騒音低減対象パネル部品モデルとして設定した。
（振動モード・等価放射パワーピーク周波数選定部）
振動モード・等価放射パワーピーク周波数選定部１９は、車体モデル２００を用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定するものである。
振動モード・等価放射パワーピーク周波数選定部１９は、車体モデル２００を用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定するものである。
振動解析としては、例えば、周波数応答解析、振動モード解析や振動エネルギー解析を適用すればよい。
振動解析により求められた振動騒音低減対象パネル部品モデルの振動挙動には、固有周波数ごとに振動モードが存在するが、振動モード・等価放射パワーピーク周波数選定部１９により選定する振動モードは、振動騒音への寄与の大きい振動モードとすればよい。
また、等価放射パワーとは、振動するパネル面から放射される音響特性を表す指標であり、音響を放射するパネル面の面積と、該パネル面の振動速度の2乗の積に比例する。
そこで、振動モード・等価放射パワーピーク周波数選定部１９は、以下の式（１）を用い、振動騒音低減対象パネル部品モデルの等価放射パワーを算出する。
そこで、振動モード・等価放射パワーピーク周波数選定部１９は、以下の式（１）を用い、振動騒音低減対象パネル部品モデルの等価放射パワーを算出する。
式（１）において、ηは放射損失係数、Cは音速、RHOは振動騒音低減対象パネル部品モデルの材料密度、Aiは振動騒音低減対象パネル部品モデルにおけるメッシュiの面積、viはメッシュiの振動速度、である。そして、振動速度viは、振動解析により各メッシュiについて求めた値を与える。
図８に車体モデル２００に生じた振動の変位を、図９にミドルフロアモデル２２１の等価放射パワーの周波数特性の結果の一例を示す。
本実施の形態では、図８に示すように、ミドルフロアモデル２２１の中央部が大きく振動する（変位が大きい）１次の振動モードを選定した。
そして、図９より、当該選定した振動モードに対応する等価放射パワーのピーク周波数（＝39.6Hz）を選定した。
本実施の形態では、図８に示すように、ミドルフロアモデル２２１の中央部が大きく振動する（変位が大きい）１次の振動モードを選定した。
そして、図９より、当該選定した振動モードに対応する等価放射パワーのピーク周波数（＝39.6Hz）を選定した。
（区分領域重量変更ピーク周波数取得部）
区分領域重量変更ピーク周波数取得部２１は、車体モデル２００における一つ又は複数の区分領域の重量を変更し、区分領域の重量の組み合わせごとに振動解析を行い、振動モード・等価放射パワーピーク周波数選定部１９による処理を行って選定した振動モードにおける振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得するものである。
表１に、車体モデル２００における一つ又は複数の区分領域の重量を変更した組み合わせ（水準）の一例を示す。
区分領域重量変更ピーク周波数取得部２１は、車体モデル２００における一つ又は複数の区分領域の重量を変更し、区分領域の重量の組み合わせごとに振動解析を行い、振動モード・等価放射パワーピーク周波数選定部１９による処理を行って選定した振動モードにおける振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得するものである。
表１に、車体モデル２００における一つ又は複数の区分領域の重量を変更した組み合わせ（水準）の一例を示す。
水準１は、区分領域D1〜区分領域D12のいずれについても重量を変更していない組み合わせである。
水準２〜水準１１は、いずれか一つの区分領域の重量を変化させた組み合わせである。
水準１２〜水準２２は、複数の区分領域の重量を同時に変化させた組み合わせである。
また、表１中の各区分領域の欄に記載の数値は、元の車体モデル２００における各区分領域の重量を基準とした重量比である。例えば、水準２は、区分領域D1.ルーフの重量について、元の車体モデル２００におけるD1.ルーフの重量の0.75倍としたものである。
水準２〜水準１１は、いずれか一つの区分領域の重量を変化させた組み合わせである。
水準１２〜水準２２は、複数の区分領域の重量を同時に変化させた組み合わせである。
また、表１中の各区分領域の欄に記載の数値は、元の車体モデル２００における各区分領域の重量を基準とした重量比である。例えば、水準２は、区分領域D1.ルーフの重量について、元の車体モデル２００におけるD1.ルーフの重量の0.75倍としたものである。
区分領域の重量の組み合わせの水準数は、車体モデル２００に設定された区分領域の総数以上であればよく、車体モデルに設定された区分領域について総組み合わせで設定してもよいし、実験計画法を用いて効率的に設定してもよい。
（区分領域重量寄与度算出部）
区分領域重量寄与度算出部２３は、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、区分領域の重量を説明変数とする多変量解析を行い、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する区分領域の重量の寄与度を算出するものである。
区分領域重量寄与度算出部２３は、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、区分領域の重量を説明変数とする多変量解析を行い、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する区分領域の重量の寄与度を算出するものである。
多変量解析に用いる回帰モデル（関係式）は、目的変数であるピーク周波数（Y）と説明変数である区分領域Diの重量（Mi）で表されるものであればよく、例えば、以下の式（２）を挙げることができる。なお、式（２）中のMiは、式（３）で示すように、重量を変更する前の区分領域Diの重量ρaを基準とした区分領域Diの重量の倍率である。
また、Miは、重量を変更する前の区分領域Diの密度ρaを基準とした区分領域Diの密度ρbの倍率としてもよい。
また、Miは、重量を変更する前の区分領域Diの密度ρaを基準とした区分領域Diの密度ρbの倍率としてもよい。
そして、多変量解析により求められた関係式（式（２））における偏回帰係数Aiの値を、ピーク周波数に対する各区分領域Diの寄与度として算出する。
なお、多変量解析での回帰モデルは、目的変数（Y）が説明変数（Mi）の線形和で表記される形式であればよく、例えば、以下の式（４）や式（５）を用いてもよい。
また、区分領域重量寄与度算出部２３により、回帰モデル（式（２）、式（４）、式（５）等）における偏回帰係数Aiを算出する際には、決定係数又は自由度調整済み決定係数を算出すると良い。決定係数又は自由度調整済み決定係数とは、回帰モデルに説明変数を代入して推定される目的変数の当てはまりの良さ（度合）を表す指標である。
本実施の形態においては、自由度調整済み決定係数が0.90以上であれば、区分領域の重量Miにより推定されるピーク周波数Yの推定において、回帰モデルの当てはまりは十分に良い。
本実施の形態においては、自由度調整済み決定係数が0.90以上であれば、区分領域の重量Miにより推定されるピーク周波数Yの推定において、回帰モデルの当てはまりは十分に良い。
（振動騒音低減軽量化部位特定部）
振動騒音低減軽量化部位特定部２５は、区分領域重量寄与度算出部２３により求められた各区分領域の寄与度に基づいて、振動騒音の低減対象とするパネル部品の振動騒音を低減するために自動車の車体における軽量化する部位を特定するものである。
すなわち、多変量解析により求めた関係式（例えば、式（２））の偏回帰係数Aiの値が大きい区分領域Diに対応する自動車における部位を、振動騒音の低減対象とするパネル部品の振動騒音を低減するために軽量化する部位として特定する。
振動騒音低減軽量化部位特定部２５は、区分領域重量寄与度算出部２３により求められた各区分領域の寄与度に基づいて、振動騒音の低減対象とするパネル部品の振動騒音を低減するために自動車の車体における軽量化する部位を特定するものである。
すなわち、多変量解析により求めた関係式（例えば、式（２））の偏回帰係数Aiの値が大きい区分領域Diに対応する自動車における部位を、振動騒音の低減対象とするパネル部品の振動騒音を低減するために軽量化する部位として特定する。
このように特定された軽量化する部位に対しては、例えば、当該部位に含まれる骨格部品及び／又はパネル部品の板厚を薄くすればよい。そして、このように軽量化することで、振動騒音を低減対象とするパネル部品からの振動騒音のピーク周波数が高くなり、起振源からの振動との共振を防ぐことができ、また、人が感じやすい音域を外すことができて、前記パネル部品からの振動騒音を低減することが可能となる。
＜自動車のパネル部品の振動騒音低減解析方法＞
本発明の実施の形態１に係る自動車のパネル部品の振動騒音低減解析方法（以下、単に「振動騒音低減解析方法」という。）の構成について、以下に説明する。
本発明の実施の形態１に係る自動車のパネル部品の振動騒音低減解析方法（以下、単に「振動騒音低減解析方法」という。）の構成について、以下に説明する。
本発明の実施の形態に係る振動騒音低減解析方法は、図２に示すように骨格部品とパネル部品とを備えてなる自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させ、かつ前記自動車の車体の軽量化できる部位を特定するものであり、図１０に示すように、車体モデル取得ステップＳ１と、区分領域設定ステップＳ３と、振動騒音低減対象パネル部品モデル設定ステップＳ５と、振動モード・等価放射パワーピーク周波数選定ステップＳ７と、区分領域重量変更ピーク周波数取得ステップＳ９と、区分領域重量寄与度算出ステップＳ１１と、振動騒音低減軽量化部位特定ステップＳ１３と、含むものである。
以下、図１０に基づいて、上記の各ステップについて説明する。なお、以下の説明では、上記の各ステップとも、コンピュータによって構成された本発明の実施の形態１に係る振動騒音低減解析装置１（図１）を用いて実行する。
以下、図１０に基づいて、上記の各ステップについて説明する。なお、以下の説明では、上記の各ステップとも、コンピュータによって構成された本発明の実施の形態１に係る振動騒音低減解析装置１（図１）を用いて実行する。
≪車体モデル取得ステップ≫
車体モデル取得ステップＳ１は、図３〜図５に一例として示すように、自動車１００（図２）の骨格部品１０１とパネル部品１０３（図２）のそれぞれをメッシュ（平面要素及び／又は立体要素）でモデル化した骨格部品モデルとパネル部品モデルを備えるとともに、起振源からの振動が入力する加振部が設定された車体モデル２００を取得するステップである。本実施の形態において、車体モデル取得ステップＳ１は、振動騒音低減解析装置１の車体モデル取得部１３が行う。
車体モデル取得ステップＳ１は、図３〜図５に一例として示すように、自動車１００（図２）の骨格部品１０１とパネル部品１０３（図２）のそれぞれをメッシュ（平面要素及び／又は立体要素）でモデル化した骨格部品モデルとパネル部品モデルを備えるとともに、起振源からの振動が入力する加振部が設定された車体モデル２００を取得するステップである。本実施の形態において、車体モデル取得ステップＳ１は、振動騒音低減解析装置１の車体モデル取得部１３が行う。
≪区分領域設定ステップ≫
区分領域設定ステップＳ３は、骨格部品モデル及びパネル部品モデルに基づいて区分した複数の区分領域を車体モデル２００に設定するステップである。本実施の形態において、区分領域設定ステップＳ３は、振動騒音低減解析装置１の区分領域設定部１５が行う。
区分領域設定ステップＳ３は、骨格部品モデル及びパネル部品モデルに基づいて区分した複数の区分領域を車体モデル２００に設定するステップである。本実施の形態において、区分領域設定ステップＳ３は、振動騒音低減解析装置１の区分領域設定部１５が行う。
区分領域設定ステップＳ３において設定される車体モデル２００の区分領域は、車体を製造する工程を考慮したASSY単位や、骨格部品モデルやパネル部品モデル単位とすればよい。あるいは、感度解析等のCAE解析手法を用いて車体モデル２００に区分領域を設定してもよい。
≪振動騒音低減対象パネル部品モデル設定ステップ≫
振動騒音低減対象パネル部品モデル設定ステップＳ５は、車体モデル２００におけるパネル部品モデルのうち振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして設定するステップである。本実施の形態において、振動騒音低減対象パネル部品モデル設定ステップＳ５は、振動騒音低減解析装置１の振動騒音低減対象パネル部品モデル設定部１７が行う。
振動騒音低減対象パネル部品モデル設定ステップＳ５は、車体モデル２００におけるパネル部品モデルのうち振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして設定するステップである。本実施の形態において、振動騒音低減対象パネル部品モデル設定ステップＳ５は、振動騒音低減解析装置１の振動騒音低減対象パネル部品モデル設定部１７が行う。
≪振動モード・等価放射パワーピーク周波数選定ステップ≫
振動モード・等価放射パワーピーク周波数選定ステップＳ７は、車体モデル２００を用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定するステップである。本実施の形態において、振動モード・等価放射パワーピーク周波数選定ステップＳ７は、振動騒音低減解析装置１の振動モード・等価放射パワーピーク周波数選定部１９が行う。
振動モード・等価放射パワーピーク周波数選定ステップＳ７は、車体モデル２００を用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定するステップである。本実施の形態において、振動モード・等価放射パワーピーク周波数選定ステップＳ７は、振動騒音低減解析装置１の振動モード・等価放射パワーピーク周波数選定部１９が行う。
振動解析としては、例えば、周波数応答解析、振動モード解析や振動エネルギー解析を適用することで、振動騒音低減対象パネル部品モデルの振動挙動を求めることができる。
振動解析により求められた振動騒音低減対象パネル部品モデルの振動挙動には、固有周波数ごとに振動モードが存在するが、振動モード・等価放射パワーピーク周波数選定ステップＳ７においては、振動騒音への寄与の大きい振動モードを選定するものとし、例えば、振動騒音低減対象パネル部品モデルの中央部が大きく振動する１次の振動モードを選定するとよい。
振動解析により求められた振動騒音低減対象パネル部品モデルの振動挙動には、固有周波数ごとに振動モードが存在するが、振動モード・等価放射パワーピーク周波数選定ステップＳ７においては、振動騒音への寄与の大きい振動モードを選定するものとし、例えば、振動騒音低減対象パネル部品モデルの中央部が大きく振動する１次の振動モードを選定するとよい。
また、等価放射パワーは、振動解析により求めた振動騒音低減対象パネル部品モデルの振動挙動の結果に基づいて、前述した式（１）により算出する。
≪区分領域重量変更ピーク周波数取得ステップ≫
区分領域重量変更ピーク周波数取得ステップＳ９は、車体モデル２００における一つ又は複数の区分領域の重量を変更し、区分領域の重量の組み合わせごとに振動解析を行い、振動モード・等価放射パワーピーク周波数選定ステップＳ７による処理を行って選定した振動モードにおける振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得するステップである。本実施の形態１において、区分領域重量変更ピーク周波数取得ステップＳ９は、振動騒音低減解析装置１の区分領域重量変更ピーク周波数取得部２１が行う。
区分領域重量変更ピーク周波数取得ステップＳ９は、車体モデル２００における一つ又は複数の区分領域の重量を変更し、区分領域の重量の組み合わせごとに振動解析を行い、振動モード・等価放射パワーピーク周波数選定ステップＳ７による処理を行って選定した振動モードにおける振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得するステップである。本実施の形態１において、区分領域重量変更ピーク周波数取得ステップＳ９は、振動騒音低減解析装置１の区分領域重量変更ピーク周波数取得部２１が行う。
車体モデル２００における一つ又は複数の区分領域の重量を変更するとは、１つの区分領域について重量を変更することや、複数の区分領域について同時に重量を変更することをいう。
そして、車体モデルにおける区分領域の重量の組み合わせの水準数は、車体モデル２００に設定された区分領域の総数以上であればよく、区分領域について総組み合わせで設定してもよいし、実験計画法を用いて効率的に設定してもよい。
≪区分領域重量寄与度算出ステップ≫
区分領域重量寄与度算出ステップＳ１１は、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、区分領域の重量を説明変数とする多変量解析を行い、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する区分領域の重量の寄与度を算出するステップである。本実施の形態１において、区分領域重量寄与度算出ステップＳ１１は、振動騒音低減解析装置１の区分領域重量寄与度算出部２３が行う。
区分領域重量寄与度算出ステップＳ１１は、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、区分領域の重量を説明変数とする多変量解析を行い、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する区分領域の重量の寄与度を算出するステップである。本実施の形態１において、区分領域重量寄与度算出ステップＳ１１は、振動騒音低減解析装置１の区分領域重量寄与度算出部２３が行う。
多変量解析により導出する関係式（回帰モデル）は、目的変数であるピーク周波数（Y）と説明変数である区分領域Diの重量（Mi）で表されるものであればよく、例えば、前述した式（２）や、式（４）又は式（５）を挙げることができる。
そして、多変量解析により求めた式（２）、式（４）又は式（５）における偏回帰係数Aiを、ピーク周波数に対する各区分領域Diの寄与度として取得する。
また、区分領域重量寄与度算出ステップＳ１１において回帰モデル（式（２）、式（４）、式（５）等）における偏回帰係数Aiを算出する際には、自由度調整済み決定係数が0.90以上であれば、区分領域の重量による等価放射パワーのピーク周波数の推定において、回帰モデルの当てはまりが十分に良い。
≪振動騒音低減軽量化部位特定ステップ≫
振動騒音低減軽量化部位特定ステップＳ１３は、区分領域重量寄与度算出ステップＳ１１により算出した各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために前記自動車の車体における軽量化する部位を特定するステップである。本実施の形態１において、振動騒音低減軽量化部位特定ステップＳ１３は、振動騒音低減解析装置１の振動騒音低減軽量化部位特定部２５が行う。
振動騒音低減軽量化部位特定ステップＳ１３は、区分領域重量寄与度算出ステップＳ１１により算出した各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために前記自動車の車体における軽量化する部位を特定するステップである。本実施の形態１において、振動騒音低減軽量化部位特定ステップＳ１３は、振動騒音低減解析装置１の振動騒音低減軽量化部位特定部２５が行う。
具体的には、区分領域重量寄与度算出ステップＳ１１における多変量解析により求めた関係式（式（２）等）の偏回帰係数Aiの値が大きい区分領域Diに対応する自動車における部位を、振動騒音を低減するために車体における軽量化する部位として特定する。
そして、特定された部位の軽量化としては、例えば、当該部位を構成する骨格部品及び／又はパネル部品の板厚を薄くして重量を軽くすればよい。そして、このような対策を施すことで、振動騒音を低減対象とするパネル部品からの振動騒音のピーク周波数を高くすることができ、起振源からの振動との共振を防ぎ、また、人が感じやすい音域を外すことで振動騒音を低減することが可能となる。
［実施の形態２］
本発明の実施の形態２に係る自動車のパネル部品の振動騒音低減解析方法は、前述した本発明の実施の形態１に係る自動車のパネル部品の振動騒音低減解析方法の区分領域重量寄与度算出ステップＳ１１において前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出した関係式に少なくとも一つ以上の前記各区分領域の重量を関係式導出時とは変更して与えて前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を予測し、自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させるために前記自動車の車体における軽量化する部位を特定するものである。
本発明の実施の形態２に係る自動車のパネル部品の振動騒音低減解析方法は、前述した本発明の実施の形態１に係る自動車のパネル部品の振動騒音低減解析方法の区分領域重量寄与度算出ステップＳ１１において前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出した関係式に少なくとも一つ以上の前記各区分領域の重量を関係式導出時とは変更して与えて前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を予測し、自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させるために前記自動車の車体における軽量化する部位を特定するものである。
これにより、予め車体モデルに設定した一つ又は複数の区分領域の重量を変更して振動騒音低減対象パネル部品モデルの等価放射パワーの周波数特性を求める振動解析を一回行って関係式を求めておくことで、自動車の骨格部品やパネル部品を含む部位の重量を減少することにより、振動騒音の低減対象とするパネル部品からの振動騒音の指標となる等価放射パワーのピーク周波数をどの程度高くすることができるかの指標を容易に得ることができる。
なお、振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と車体モデルに設定された区分領域の重量との関係式は、前述の実施の形態１で示した式（２）や式（４）、式（５）のいずれであってもよく、特に限定はない。
［実施の形態３］
本発明の実施の形態３に係る自動車のパネル部品の振動騒音低減解析方法は、前述した実施の形態１の区分領域重量寄与度算出ステップＳ１１において前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出した関係式を用いて板厚を算出したい一つの区分領域について、他の区分領域の重量と前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数とを与えて、前記一つの区分領域の重量を算出し、該算出した前記一つの区分領域の重量に基づいて、該一つの区分領域を構成する前記骨格部品モデル及び／又はパネル部品モデルの板厚を算出するものである。
本発明の実施の形態３に係る自動車のパネル部品の振動騒音低減解析方法は、前述した実施の形態１の区分領域重量寄与度算出ステップＳ１１において前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出した関係式を用いて板厚を算出したい一つの区分領域について、他の区分領域の重量と前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数とを与えて、前記一つの区分領域の重量を算出し、該算出した前記一つの区分領域の重量に基づいて、該一つの区分領域を構成する前記骨格部品モデル及び／又はパネル部品モデルの板厚を算出するものである。
このように、等価放射パワーのピーク周波数と区分領域の重量の関係式を用いることで、振動騒音の低減対策として軽量化する一つ又は複数の部位を選定しやすくなるだけでなく、振動騒音の低減対策として自動車の車体における一つ又は複数の部位の重量を減らすことによる等価放射パワーのピーク周波数の向上効果を推定することが可能となり、パネル部品からの振動騒音の低減と自動車の軽量化の双方を達成した車体設計を容易に行うことが可能となる。
なお、等価放射パワーのピーク周波数と車体モデルに設定された区分領域の重量との関係式は、前述した実施の形態２と同様に、振動騒音低減対象パネル部品モデルの前述の実施の形態１で示した式（２）や式（４）、式（５）のいずれであってもよく、特に限定はない。
もっとも、関係式に与える等価放射パワーのピーク周波数が、各区分領域の重量を変更する前の車体モデルを用いて求めた等価放射パワーのピーク周波数に比べて高すぎる場合は、区分領域の重量変更に伴う剛性の変化（例えば板厚を薄くしすぎて部品の剛性が低下するなど）が振動騒音低減対象パネル部品モデルの振動挙動に大きく影響しすぎて、式（２）等の関係式における等価放射パワーのピーク周波数（Y）と区分領域の重量（Mi）との相関（予測値の精度）が低下する場合がある。
そこで、本実施の形態３に係る自動車のパネル部品の振動騒音低減解析方法の他の態様として、等価放射パワーのピーク周波数と区分領域の重量との関係式である式（２）に、重量を変更する区分領域Diの軽量化による等価放射パワーのピーク周波数への影響を補正する補正値αiを与えた式（６）に、等価放射パワーのピーク周波数と、板厚を算出したい一つの区分領域を除いた他の区分領域の重量と、を与え、一つの区分領域の重量を算出することが好ましい。
この場合、区分領域の板厚を予め数条件変更して関係式から算出された予測値と、振動解析で得られた計算値との残差平方和が最小となる補正値αiを求めて設定するとよい。
この場合、区分領域の板厚を予め数条件変更して関係式から算出された予測値と、振動解析で得られた計算値との残差平方和が最小となる補正値αiを求めて設定するとよい。
補正値αiは、重量を求めるために選択する一つの区分領域ごとに、上述の方法等により求めることが望ましい。
式（６）において、他の区分領域の重量に関しては補正せずにαi＝1を与えてもよいし、大きく変更した重量を与える他の区分領域に関しては補正値αiを適宜与えてもよい。
式（６）において、他の区分領域の重量に関しては補正せずにαi＝1を与えてもよいし、大きく変更した重量を与える他の区分領域に関しては補正値αiを適宜与えてもよい。
このように、本実施の形態３に係る自動車のパネル部品の振動騒音低減解析方法の他の態様により、振動騒音の低減対策を施す部位の重量を減らしたことによる剛性の変化を考慮しつつ当該部位の板厚を精度良く推算し、パネル部品の振動騒音の低減と車体の軽量化の双方を達成することができる。
本発明に係る自動車のパネル部品の振動騒音低減解析装置及び解析方法の効果を検証するための数値実験を行ったので、以下、これについて説明する。
数値実験では、図３〜図５に示すように、自動車１００の骨格部品１０１とパネル部品１０３（図２）のそれぞれをメッシュ（平面要素及び／又は立体要素）でモデル化したパネル部品モデルと骨格部品モデルとともに、自動車１００のリアサスペンション取付部１０７（図２）に対応する部位を加振部２３１として設定した車体モデル２００を解析対象とした。
そして、自動車１００のミドルフロア１０３ａを振動騒音の低減対象とし、車体モデル２００におけるミドルフロアモデル２２１を振動騒音低減対象パネル部品モデルとした。
なお、車体モデル２００における骨格部品モデル及びパネル部品モデルには、いずれも、表２に示す物性値を設定した。
そして、自動車１００のミドルフロア１０３ａを振動騒音の低減対象とし、車体モデル２００におけるミドルフロアモデル２２１を振動騒音低減対象パネル部品モデルとした。
なお、車体モデル２００における骨格部品モデル及びパネル部品モデルには、いずれも、表２に示す物性値を設定した。
さらに、車体モデル２００に対し、図４及び図５に示すように、区分領域D1〜D12を設定した。
各区分領域D1〜D12を構成する骨格部品モデル及びパネル部品モデルについては、実施の形態１で述べたとおりである。
そして、自動車１００のミドルフロア１０３ａからの振動騒音の低減と車体の軽量化をするための解析として、Ａ．振動騒音低減軽量化部位の特定、Ｂ．車体モデルにおけるミドルフロアモデルの等価放射パワーのピーク周波数の予測、及び、Ｃ．振動騒音低減軽量化部位の板厚の算出、を行った。以下、Ａ〜Ｃについて、順に説明する。
各区分領域D1〜D12を構成する骨格部品モデル及びパネル部品モデルについては、実施の形態１で述べたとおりである。
そして、自動車１００のミドルフロア１０３ａからの振動騒音の低減と車体の軽量化をするための解析として、Ａ．振動騒音低減軽量化部位の特定、Ｂ．車体モデルにおけるミドルフロアモデルの等価放射パワーのピーク周波数の予測、及び、Ｃ．振動騒音低減軽量化部位の板厚の算出、を行った。以下、Ａ〜Ｃについて、順に説明する。
［Ａ．振動騒音低減軽量化部位の特定］
車体モデル２００のミドルフロアモデル２２１の振動騒音の低減に寄与する区分領域である自動車における車体を軽量化する部位を特定した。
まず、車体モデル２００を用いて周波数応答解析を行い、ミドルフロアモデルの振動挙動と等価放射パワーの周波数特性を求めた。
そして、振動解析の結果から、前述の実施の形態１で示したように、ミドルフロアモデルの中央部が大きく振動する振動モードを選定し、該振動モードに対応する等価放射パワーのピーク周波数を取得した。
車体モデル２００のミドルフロアモデル２２１の振動騒音の低減に寄与する区分領域である自動車における車体を軽量化する部位を特定した。
まず、車体モデル２００を用いて周波数応答解析を行い、ミドルフロアモデルの振動挙動と等価放射パワーの周波数特性を求めた。
そして、振動解析の結果から、前述の実施の形態１で示したように、ミドルフロアモデルの中央部が大きく振動する振動モードを選定し、該振動モードに対応する等価放射パワーのピーク周波数を取得した。
続いて、一つ又は複数の区分領域D1〜D12の重量を変更して周波数応答解析を行い、区分領域の重量の組み合わせごとに、等価放射パワーのピーク周波数を求めた。ここで、区分領域D1〜D12の重量の組み合わせの水準数は52とした。前述の表１に、区分領域の重量の組み合わせの一部を示す。
そして、振動解析の結果に基づいて求めた等価放射パワーのピーク周波数（Y）を目的変数とし、区分領域Diの重量（Mi）を説明変数として多変量解析を行い、等価放射パワーのピーク周波数Yと区分領域の重量Miとの関係式を求めた。
ここで、多変量解析として重回帰分析を行い、回帰モデルには式（２）を用いた。
多変量解析により求めた関係式を式（７）に示す。また、表３に、式（７）における偏回帰係数Aiと、多変量解析により求めた標準偏回帰係数をあわせて示す。
ここで、多変量解析として重回帰分析を行い、回帰モデルには式（２）を用いた。
多変量解析により求めた関係式を式（７）に示す。また、表３に、式（７）における偏回帰係数Aiと、多変量解析により求めた標準偏回帰係数をあわせて示す。
式（７）及び表３より、偏回帰係数Aiの値が大きい区分領域D8（ミドルフロア、A8＝1.81）と区分領域D10（ミドルフロアメンバ、A10＝3.57）の重量は、ミドルフロアモデル２２１の等価放射パワーのピーク周波数に対する寄与度が大きいことが分かる。
これに対し、偏回帰係数の値の小さい区分領域、例えば、D1（ルーフ、A1＝-0.07）やD11（リアフロアサイドメンバ、A11＝0.001）の重量は、ミドルフロアモデル２２１の等価放射パワーのピーク周波数に対する寄与度が小さいことが分かる。
以上の結果から、区分領域D8.ミドルフロア及び区分領域D10.ミドルフロアメンバのそれぞれに対応する車体における部位を、振動騒音の低減対策として軽量化する部位として特定することができた。
［Ｂ．等価放射パワーのピーク周波数の予測］
次に、等価放射パワーのピーク周波数の予測を行った。
等価放射パワーのピーク周波数と区分領域の重量との関係式として、前述した式（７）を用いた。
そして、車体モデル２００（図４及び図５参照）におけるミドルフロアモデル２２１の等価放射パワーのピーク周波数への寄与の高い区分領域であるD10.ミドルフロアメンバと、寄与の低い区分領域であるD11.リアフロアサイドメンバのそれぞれの重量を変更し、等価放射パワーのピーク周波数（予測値）を求めた。
一方、比較対象として、区分領域D10.ミドルフロアメンバ、又は、区分領域D11.リアフロアサイドメンバの重量を変更した車体モデル２００について周波数応答解析を行い、ミドルフロアモデル２２１の等価放射パワーの周波数応答を求め、ミドルフロアモデル２２１の中央部が大きく変形する振動モード（１次モード）となるピーク周波数（計算値）を求めた。
表４に、区分領域D10.ミドルフロアメンバ及び区分領域D11.リアフロアサイドメンバの重量と、式（７）によりピーク周波数を予測した結果（予測値）を示す。また、表４には、元の車体モデル２００を用いてミドルフロアモデル２２１の等価放射パワーのピーク周波数を求めた結果（計算値）もあわせて示す。
次に、等価放射パワーのピーク周波数の予測を行った。
等価放射パワーのピーク周波数と区分領域の重量との関係式として、前述した式（７）を用いた。
そして、車体モデル２００（図４及び図５参照）におけるミドルフロアモデル２２１の等価放射パワーのピーク周波数への寄与の高い区分領域であるD10.ミドルフロアメンバと、寄与の低い区分領域であるD11.リアフロアサイドメンバのそれぞれの重量を変更し、等価放射パワーのピーク周波数（予測値）を求めた。
一方、比較対象として、区分領域D10.ミドルフロアメンバ、又は、区分領域D11.リアフロアサイドメンバの重量を変更した車体モデル２００について周波数応答解析を行い、ミドルフロアモデル２２１の等価放射パワーの周波数応答を求め、ミドルフロアモデル２２１の中央部が大きく変形する振動モード（１次モード）となるピーク周波数（計算値）を求めた。
表４に、区分領域D10.ミドルフロアメンバ及び区分領域D11.リアフロアサイドメンバの重量と、式（７）によりピーク周波数を予測した結果（予測値）を示す。また、表４には、元の車体モデル２００を用いてミドルフロアモデル２２１の等価放射パワーのピーク周波数を求めた結果（計算値）もあわせて示す。
表４において、No.1は、元の車体モデル２００を用いて周波数応答解析行った場合、No.2は、区分領域D11.リアフロアメンバの重量を0.5倍とした場合、No.3は、区分領域D10.ミドルフロアメンバの重量を0.5倍とした場合である。
表４に示すように、式（７）により求めた等価放射パワーのピーク周波数の予測値は、車体モデル２００について振動解析により求めた等価放射パワーのピーク周波数の計算値とほぼ等しい結果であった。
寄与度の低い区分領域D1リアフロアサイドメンバの重量を変更（0.5倍）した場合（No.2）、等価放射パワーのピーク周波数の予測値及び計算値はともに、リアフロアサイドメンバの重量を変更する前の車体モデル２００の振動解析により求めた等価放射パワーのピーク周波数とほぼ同等であった。
一方、寄与度の高い区分領域D10.ミドルフロアメンバの重量を変更した場合、等価放射パワーのピーク周波数の予測値及び計算値はともに、元の車体モデルのピーク周波数よりも3.2Hz及び3.5Hz高い結果となった。
これよりも、ミドルフロアモデルの等価放射パワーのピーク周波数への寄与度の低いD11.リアフロアサイドメンバの重量を変更しても振動騒音のピーク周波数はほぼ変わらないのに対し、寄与度の高いD10.ミドルフロアメンバの重量を0.5倍とすることで、ピーク周波数を3Hz以上高くでき、車体の軽量化と振動騒音の低減の双方を達成できることが示された。
［振動騒音低減軽量化部位の板厚の算出］
続いて、振動騒音低減軽量化部位のピーク周波数に対する板厚の影響について検討した。
ここでは、等価放射パワーのピーク周波数と区分領域の重量との関係式として、前述した式（７）を用い、区分領域D10.ミドルフロアメンバモデルの板厚を変更したときのミドルフロアモデル２２１（図４及び図５参照）の等価放射パワーのピーク周波数を予測した。
さらに、区分領域D10.ミドルフロアメンバモデルを軽量化したことによる剛性変化を考慮する補正値αを与えた場合についても、ピーク周波数を予測した。
表５に、ピーク周波数の予測値と計算値を示す。
続いて、振動騒音低減軽量化部位のピーク周波数に対する板厚の影響について検討した。
ここでは、等価放射パワーのピーク周波数と区分領域の重量との関係式として、前述した式（７）を用い、区分領域D10.ミドルフロアメンバモデルの板厚を変更したときのミドルフロアモデル２２１（図４及び図５参照）の等価放射パワーのピーク周波数を予測した。
さらに、区分領域D10.ミドルフロアメンバモデルを軽量化したことによる剛性変化を考慮する補正値αを与えた場合についても、ピーク周波数を予測した。
表５に、ピーク周波数の予測値と計算値を示す。
表５において、No.1は、元の車体モデル２００におけるD10.ミドルフロアメンバの重量とした場合、No.4及びNo.5は、それぞれ、区分領域D10.ミドルフロアメンバを構成する部品の板厚を1.2mmから0.9mm及び0.8mmに薄くして軽量化した場合である。
表５よりも、区分領域D10.ミドルフロアメンバの板厚を変更したことによりピーク周波数の予測値と計算値は、いずれも、元の車体よりも高くなり、板厚を薄くするほど、すなわち、重量を減少するほどピーク周波数は高くなった。
また、予測値に関して、前述の式（６）の軽量化による剛性変化の補正値αiを与えることで、計算値に近づき良好な結果となった。
また、予測値に関して、前述の式（６）の軽量化による剛性変化の補正値αiを与えることで、計算値に近づき良好な結果となった。
上記の結果は、振動騒音低減対象パネル部品モデルの振動騒音への寄与度の大きい区分領域の板厚を変更してピーク周波数を求めるものであったが、ピーク周波数を与えて区分領域の板厚を求める場合であっても、所定のピーク周波数を満たすように板厚を求めることができ、振動騒音の低減と軽量化を達成することができる。さらに、区分領域の重量変更による剛性変化に係る補正値を与えることで、区分領域の板厚をより精度よく推算することが可能となる。
１ 振動騒音低減解析装置
３ 表示装置
５ 入力装置
７ 記憶装置
９ 作業用データメモリ
１１ 演算処理部
１３ 車体モデル取得部
１５ 区分領域設定部
１７ 振動騒音低減対象パネル部品モデル設定部
１９ 振動モード・等価放射パワーピーク周波数選定部
２１ 区分領域重量変更ピーク周波数取得部
２３ 区分領域重量寄与度算出部
２５ 振動騒音低減軽量化部位特定部
３０ 車体モデルファイル
１００ 自動車
１０１ 骨格部品
１０３ パネル部品
１０３ａ ミドルフロア
１０５ 起振源
１０５ａ エンジン
１０５ｂ タイヤ
１０６ フロントサスペンション取付部
１０７ リアサスペンション取付部
２００ 車体モデル
２０１ ルーフレールモデル
２０３ Ａピラモデル
２０５ Ｂピラモデル
２０７ Ｃピラモデル
２０９ ミドルフロアメンバモデル
２１１ リアフロアサイドメンバモデル
２１３ リアフロアクロスメンバモデル
２２１ ミドルフロアモデル
２２３ ルーフモデル
２２５ リアフロアモデル
２３１ 加振部
２４１ アッパバックパネルモデル
２４３ ロアバックパネルモデル
２４５ クウォータインナモデル
２４７ ホイールハウスモデル
３ 表示装置
５ 入力装置
７ 記憶装置
９ 作業用データメモリ
１１ 演算処理部
１３ 車体モデル取得部
１５ 区分領域設定部
１７ 振動騒音低減対象パネル部品モデル設定部
１９ 振動モード・等価放射パワーピーク周波数選定部
２１ 区分領域重量変更ピーク周波数取得部
２３ 区分領域重量寄与度算出部
２５ 振動騒音低減軽量化部位特定部
３０ 車体モデルファイル
１００ 自動車
１０１ 骨格部品
１０３ パネル部品
１０３ａ ミドルフロア
１０５ 起振源
１０５ａ エンジン
１０５ｂ タイヤ
１０６ フロントサスペンション取付部
１０７ リアサスペンション取付部
２００ 車体モデル
２０１ ルーフレールモデル
２０３ Ａピラモデル
２０５ Ｂピラモデル
２０７ Ｃピラモデル
２０９ ミドルフロアメンバモデル
２１１ リアフロアサイドメンバモデル
２１３ リアフロアクロスメンバモデル
２２１ ミドルフロアモデル
２２３ ルーフモデル
２２５ リアフロアモデル
２３１ 加振部
２４１ アッパバックパネルモデル
２４３ ロアバックパネルモデル
２４５ クウォータインナモデル
２４７ ホイールハウスモデル
Claims (5)
- 自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させ、かつ前記自動車の車体の軽量化ができる部位を特定する自動車のパネル部品の振動騒音低減解析装置であって、
前記自動車の車体を構成する骨格部品とパネル部品のそれぞれをメッシュでモデル化した骨格部品モデルとパネル部品モデルを備えてなり、前記起振源からの振動を入力する加振部が設定された車体モデルを取得する車体モデル取得部と、
該車体モデルに、前記骨格部品モデル及び前記パネル部品モデルに基づいて区分した複数の区分領域を設定する区分領域設定部と、
操作者からの指示により、前記パネル部品モデルのうち振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして前記車体モデルに設定する振動騒音低減対象パネル部品モデル設定部と、
前記車体モデルを用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定する振動モード・等価放射パワーピーク周波数選定部と、
前記車体モデルにおける一つ又は複数の区分領域の重量を変更し、前記車体モデルにおける前記区分領域の重量の組み合わせごとに振動解析を行い、前記振動モード・等価放射パワーピーク周波数選定部による処理を行って選定した振動モードにおける前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得する区分領域重量変更ピーク周波数取得部と、
前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、前記区分領域の重量を説明変数とする多変量解析を行い、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する前記区分領域の重量の寄与度を算出する区分領域重量寄与度算出部と、
該算出した前記各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために前記自動車の車体における軽量化する部位を特定する振動騒音低減軽量化部位特定部と、を備えたことを特徴とする自動車のパネル部品の振動騒音低減解析装置。 - 自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させ、かつ前記自動車の車体の軽量化ができる部位を特定するために、コンピュータが以下の各ステップを行う自動車のパネル部品の振動騒音低減解析方法であって、
前記自動車の車体を構成する骨格部品とパネル部品のそれぞれをメッシュでモデル化した骨格部品モデルとパネル部品モデルを備えてなり、前記起振源からの振動を入力する加振部が設定された車体モデルを取得する車体モデル取得ステップと、
該車体モデルに、前記骨格部品モデル及び前記パネル部品モデルに基づいて区分した複数の区分領域を設定する区分領域設定ステップと、
前記パネル部品モデルのうち、振動騒音の低減対象とするパネル部品モデルを振動騒音低減対象パネル部品モデルとして前記車体モデルに設定する振動騒音低減対象パネル部品モデル設定ステップと、
前記車体モデルを用いて振動解析を行い、振動騒音低減対象パネル部品モデルの振動挙動と振動騒音の指標である等価放射パワーの周波数特性とを求め、振動騒音への寄与の大きい振動モードと、該振動モードに対応する等価放射パワーのピーク周波数を選定する振動モード・等価放射パワーピーク周波数選定ステップと、
前記車体モデルにおける一つ又は複数の区分領域の重量を変更し、前記車体モデルにおける前記区分領域の重量の組み合わせごとに振動解析を行い、前記振動モード・等価放射パワーピーク周波数選定ステップによる処理を行って選定した振動モードにおける前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を取得する区分領域重量変更ピーク周波数取得ステップと、
前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を目的変数とし、前記区分領域の重量を説明変数とする多変量解析を行い、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数に対する前記区分領域の重量の寄与度を算出する区分領域重量寄与度算出ステップと、
該算出した前記各区分領域の寄与度に基づいて、振動騒音の低減対象とする前記パネル部品からの振動騒音を低減するために前記自動車の車体における軽量化する部位を特定する振動騒音低減軽量化部位特定ステップと、を備えたことを特徴とする自動車のパネル部品の振動騒音低減解析方法。 - 請求項２に記載の自動車のパネル部品の振動騒音低減解析方法の前記区分領域重量寄与度算出ステップにおいて、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出された関係式に少なくとも一つ以上の前記各区分領域の重量を変更して与えることにより前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数を予測し、自動車における起振源からの振動に起因するパネル部品の振動騒音を低減させるために前記自動車の車体における軽量化する部位を特定することを特徴とする自動車のパネル部品の振動騒音低減解析方法。
- 請求項２に記載の自動車のパネル部品の振動騒音低減解析方法の前記区分領域重量寄与度算出ステップにおいて、前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数と前記各区分領域の重量との関係式を導出し、該導出された関係式に板厚を算出したい一つの区分領域を除いた他の区分領域の重量と前記振動騒音低減対象パネル部品モデルの等価放射パワーのピーク周波数とを与えて、前記選択した一つの区分領域の重量を算出し、該算出した前記一つの区分領域の重量に基づいて、該一つの区分領域を構成する前記骨格部品モデル及び／又はパネル部品モデルの板厚を算出することを特徴とする自動車のパネル部品の振動騒音低減解析方法。
- 前記関係式により前記選択した一つの区分領域の重量を算出するに際し、該一つの区分領域の重量の変化が前記等価放射パワーのピーク周波数に及ぼす影響を補正する補正値を前記関係式に与え、前記一つの区分領域の重量を算出することを特徴とする請求項４記載の自動車のパネル部品の振動騒音低減解析方法。
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020113763A JP6977817B1 (ja) | 2020-07-01 | 2020-07-01 | 自動車のパネル部品の振動騒音低減解析装置及び解析方法 |
| PCT/JP2021/014571 WO2022004087A1 (ja) | 2020-07-01 | 2021-04-06 | 自動車のパネル部品の振動騒音低減解析装置及び解析方法 |
| MX2023000182A MX2023000182A (es) | 2020-07-01 | 2021-04-06 | Dispositivo de analisis de reduccion de vibracion y ruido y metodo de analisis para parte de panel de automovil. |
| KR1020227043521A KR20230005397A (ko) | 2020-07-01 | 2021-04-06 | 자동차의 패널 부품의 진동 소음 저감 해석 장치 및 해석 방법 |
| EP21833592.5A EP4148611A4 (en) | 2020-07-01 | 2021-04-06 | DEVICE FOR ANALYZING VIBRATIONS/NOISE REDUCTION OF AN AUTOMOBILE PANEL COMPONENT AND ANALYZING METHOD |
| US18/008,936 US20230237216A1 (en) | 2020-07-01 | 2021-04-06 | Vibration and noise reduction analysis device and analysis method for panel part of automobile |
| CN202180044836.XA CN115803249A (zh) | 2020-07-01 | 2021-04-06 | 汽车的板部件的振动噪音降低解析装置及解析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020113763A JP6977817B1 (ja) | 2020-07-01 | 2020-07-01 | 自動車のパネル部品の振動騒音低減解析装置及び解析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP6977817B1 true JP6977817B1 (ja) | 2021-12-08 |
| JP2022012149A JP2022012149A (ja) | 2022-01-17 |
Family
ID=78815513
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020113763A Active JP6977817B1 (ja) | 2020-07-01 | 2020-07-01 | 自動車のパネル部品の振動騒音低減解析装置及び解析方法 |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US20230237216A1 (ja) |
| EP (1) | EP4148611A4 (ja) |
| JP (1) | JP6977817B1 (ja) |
| KR (1) | KR20230005397A (ja) |
| CN (1) | CN115803249A (ja) |
| MX (1) | MX2023000182A (ja) |
| WO (1) | WO2022004087A1 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117113786B (zh) * | 2023-10-25 | 2024-02-23 | 泉州装备制造研究所 | 一种基于动力学振动响应分析的经编机的优化方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010230584A (ja) * | 2009-03-27 | 2010-10-14 | Mazda Motor Corp | 車両の音響解析装置 |
| JP2019114114A (ja) * | 2017-12-25 | 2019-07-11 | トヨタ自動車株式会社 | 車両のパネル部材の振動解析装置 |
| WO2020070922A1 (ja) * | 2018-10-05 | 2020-04-09 | Ｊｆｅスチール株式会社 | 車体部品の感度解析方法及び装置、車体部品の材料特性決定方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4752517B2 (ja) | 2006-01-13 | 2011-08-17 | 日産自動車株式会社 | 車体のルーフ構造 |
-
2020
- 2020-07-01 JP JP2020113763A patent/JP6977817B1/ja active Active
-
2021
- 2021-04-06 CN CN202180044836.XA patent/CN115803249A/zh active Pending
- 2021-04-06 EP EP21833592.5A patent/EP4148611A4/en not_active Withdrawn
- 2021-04-06 WO PCT/JP2021/014571 patent/WO2022004087A1/ja unknown
- 2021-04-06 US US18/008,936 patent/US20230237216A1/en active Pending
- 2021-04-06 KR KR1020227043521A patent/KR20230005397A/ko unknown
- 2021-04-06 MX MX2023000182A patent/MX2023000182A/es unknown
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010230584A (ja) * | 2009-03-27 | 2010-10-14 | Mazda Motor Corp | 車両の音響解析装置 |
| JP2019114114A (ja) * | 2017-12-25 | 2019-07-11 | トヨタ自動車株式会社 | 車両のパネル部材の振動解析装置 |
| WO2020070922A1 (ja) * | 2018-10-05 | 2020-04-09 | Ｊｆｅスチール株式会社 | 車体部品の感度解析方法及び装置、車体部品の材料特性決定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20230237216A1 (en) | 2023-07-27 |
| KR20230005397A (ko) | 2023-01-09 |
| JP2022012149A (ja) | 2022-01-17 |
| WO2022004087A1 (ja) | 2022-01-06 |
| EP4148611A4 (en) | 2023-11-08 |
| EP4148611A1 (en) | 2023-03-15 |
| CN115803249A (zh) | 2023-03-14 |
| MX2023000182A (es) | 2023-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2021065064A1 (ja) | 自動車のパネル部品の振動騒音低減解析方法及び解析装置 | |
| US20100299107A1 (en) | Acoustic analysis apparatus for vehicle | |
| Lee et al. | Structural optimization of an automotive door using the kriging interpolation method | |
| JP6935743B2 (ja) | 車両のパネル部材の振動解析装置 | |
| JP6977817B1 (ja) | 自動車のパネル部品の振動騒音低減解析装置及び解析方法 | |
| Abdullah et al. | Computational modal analysis on finite element model of body-in-white structure and its correlation with experimental data | |
| Zargar et al. | Investigation and optimization for the dynamical behaviour of the vehicle structure | |
| WO2021053872A1 (ja) | 自動車のパネル部品の振動騒音低減解析方法及び解析装置 | |
| Schäffer et al. | Methodological approach for reducing computational costs of vehicle frontal crashworthiness analysis by using simplified structural modelling | |
| Dhandole et al. | On improving weakly coupled cavity models for vibro-acoustic predictions and design | |
| Tebby et al. | Optimization of vehicle structure considering torsion stiffness using simple structural beam frame-approach | |
| JP7327424B2 (ja) | 自動車の車体部品の板厚感度解析方法、板厚感度解析装置及び板厚感度解析プログラム | |
| Butkunas | Analysis of vehicle structures | |
| Donders et al. | Cae technologies for efficient vibro-acoustic vehicle design modification and optimization | |
| Christian | Using Frequency Dependent Causality Analysis and Automated Tuning with Broadband ANC Systems to Optimize the Performance of the 3D Sound Field in a Passenger Vehicle | |
| Huang et al. | Development of a luxury vehicle acoustic package using SEA full vehicle model | |
| Galante et al. | Acoustic Package Development and Target Setting of an Autonomous Vehicle | |
| Pedersen | Conceptual Dynamic Analysis of a Vehicle Body | |
| Vlahopoulos et al. | Design of Rotorcraft Gearbox Foundation for Reduced Vibration and Increased Crashworthiness Characteristics | |
| Haider | Efficient reanalysis methods for size, shape and topology optimization of linear dynamic systems | |
| Echempati et al. | Integrated metal forming and vibration analysis of sheet metal parts | |
| Tebby et al. | Analysis and Optimization of Automotive Structure’s Bending Stiffness Using Beam Elements | |
| CN117313342A (zh) | 车辆噪声计算方法、车载设备和车辆 | |
| Moron et al. | Use of SEA for vehicle target setting and efficient realization of vehicle acoustic goals | |
| Liang et al. | Numerical Study of the Noise in Compartment after Auto-Body Lightweighting |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210714 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210714 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20210714 |
|
| 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: 20211012 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20211025 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6977817 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |