(57)【要約】【目的】現実的に実現可能な時間内で実用的かつ大規模な振動の解析が行なえるような方法を提供することを目的とする。【構成】従来の方法によれば，実際の問題に応用可能な大規模な構造・音場連成系問題を解くためには極めて長時間にわたる計算時間と多くのコンピュータ資源（メモリやハードディスク）が必要となる。本発明では構造と音場連成系を分解し，構造と音場それぞれの固有ペア（固有値と固有モード）の計算結果を入力として，求める精度に応じた最小限の固有ペアの数で連成系の近似固有値を求めることが可能となり，従来より短時間で計算を行なうことが可能となる。【選択図】　　図３
2004094331.tif 000002
【特許請求の範囲】
【請求項１】
振動を解析すべき領域の音場と構造のそれぞれの非連成振動方程式から正規化条件としてＫ直交条件を用いて算出された固有値と固有モードをコンピュータに入力してメモリに記憶し，
メモリに記憶した固有値と固有モードの固有ペアについて，固有ペアにおける固有値と所与の基準固有値との差の絶対値が小さい順番にメモリ上で並び替えを行い，
振動を解析すべき上記領域の音場と構造の連成系における固有値と固有モードを連成系の結合パラメータεの級数で表現した摂動法による展開式の係数をコンピュータにより算出するにあたって，上記係数を，離散化された連成系振動方程式に基づいて，連成系振動方程式内の固有値を周波数の逆数で置きかえた上で，連成系振動方程式を対称化することなく，非対称のまま，上記の並び替えされた固有ペアを用いて算出して，最終的に上記連成系の固有ペアを導出する，
構造物の設計支援方法。
【請求項２】
上記コンピュータのメモリに記憶された固有ペア内の固有値において，上記基準固有値を指定する，請求項１に記載の方法。
【請求項３】
許容誤差ｅをコンピュータに入力し，上記コンピュータによる上記係数の近似計算において，前回の計算結果と今回の計算結果との差が上記許容誤差ｅよりも小さくなったときに近似計算を停止する，請求項１または２に記載の方法。
【請求項４】
打切り回数ｍをコンピュータに入力し，上記コンピュータによる上記係数の近似計算において，前回の計算結果と今回の計算結果との差が上記許容誤差ｅよりも大きい場合であっても，繰返し計算回数が上記打切り回数ｍに達した場合に，近似計算を停止する，請求項３に記載の方法。
【請求項５】
振動を解析すべき領域の音場と構造のそれぞれの非連成振動方程式から正規化条件としてＫ直交条件を用いて算出された固有値と固有モードの入力を受入れてメモリに記憶する入力手段，
上記メモリに記憶される固有値と固有モードの固有ペアについて，固有ペアにおける固有値と所与の基準固有値との差の絶対値が小さい順番にメモリ上で並び替えを行う並び替え手段，
振動を解析すべき上記領域の音場と構造の連成系における固有値と固有モードを連成系の結合パラメータεの級数で表現した摂動法による展開式の係数をコンピュータにより算出するにあたって，上記係数を，離散化された連成系振動方程式に基づいて，連成系振動方程式内の固有値を周波数の逆数で置きかえた上で，連成系振動方程式を対称化することなく，非対称のまま，上記の並び替えされた固有ペアを用いて算出して，最終的に上記連成系の固有ペアを導出する処理手段，ならびに
導出された上記固有ペアを出力する出力手段，
を備えた構造物の設計支援システム。
【請求項６】
振動を解析すべき領域の音場と構造のそれぞれの非連成振動方程式から正規化条件としてＫ直交条件を用いて算出された固有値と固有モードの入力を受付けてメモリに記憶し，
メモリに記憶した固有値と固有モードの固有ペアについて，固有ペアにおける固有値と所与の基準固有値との差の絶対値が小さい順番にメモリ上で並び替えを行い，
振動を解析すべき上記領域の音場と構造の連成系における固有値と固有モードを連成系の結合パラメータεの級数で表現した摂動法による展開式の係数を算出するにあたって，上記係数を，離散化された連成系振動方程式に基づいて，連成系振動方程式内の固有値を周波数の逆数で置きかえた上で，連成系振動方程式を対称化することなく，非対称のまま，上記の並び替えされた固有ペアを用いて算出して，最終的に上記連成系の固有ペアを導出するようにコンピュータを制御する連成系振動解析のためのプログラム。
【発明の詳細な説明】
【０００１】
【技術分野】
この発明は振動のふるまいを考慮した，構造物の設計を支援するシステムおよび方法，または音場と構造の連成系振動解析システムおよび方法，ならびにそのシステムを制御するプログラムに関する。
【０００２】
【技術的背景】
さまざまな構造物の設計において振動を考慮することはきわめて重要である。たとえば自動車の設計においては，乗車時の快適性を向上させるため振動を極力抑えることが要求される。またスピーカの設計においてはスピーカに対する要求仕様に従い，必要な周波数帯域において最大の電気／音響変換効率が得られるような特性が求められる。
【０００３】
ここで自動車の設計を例にとると，設計仕様として要求される項目としては自動車の運動性能の向上，製造コストの削減等があげられるが，これに加えて快適性確保のため，車室内における振動（騒音）を低減するための対策は欠かすことのできない重要な要求項目のひとつである。自動車の車室内でシートに座った人間が感じる振動（騒音）は，いくつかの条件によって決まる。音場の条件すなわち車室内空間に関する条件，構造の条件すなわち車体やシートに関する条件，等である。音場の条件としては，具体的には車室内の広さ（容積），空間の形状等があげられる。構造の条件としては，車体を構成する部材の材質や剛性，部材の取付け方法等がある。車室内の空間は車体によって規定されるなど，これらの条件は相互に関連している。車室内の振動（騒音）はこのような音場の条件と構造の条件に依存する。このように相互に関連する複数の条件に依存して生じる振動を連成振動現象という。このような騒音ないしは振動現象を解析するための方法は，自動車等の設計にとって極めて重要な手法であり，その設計のための支援と位置付けられる。
【０００４】
数値計算において連成振動現象を表現する行列は非対称形である。非対称形の行列をコンピュータを用いて解くのには膨大な時間がかかる。計算を容易にする方法として非対称形行列を対称形行列に誘導する方法があるが，対称形行列を誘導できたとしても対称形行列にはバンド性がないためにその固有値解析に要する計算時間は長く，結局，行列を解くには時間がかかる。ここでバンド性とは，行列において０（ゼロ）でない項が行列の対角線に集中する性質をいう。
【０００５】
振動には固有値と固有モード（固有ベクトルともいう）がある。固有値と固有モードを合わせて固有ペアと呼ぶ。振動の解析時間短縮のために，非連成の固有ペアを用いて連成系の固有ペアを表現する試みがなされているが，途中で発散して解けないか，解けたとして時間がかかり実用的ではなかった。
【０００６】
【発明の開示】
この発明は，現実的に実行が可能な時間内で振動（騒音）解析を行うことができる構造物の設計支援システムないしは振動の解析システムおよび方法，ならびに同システムのためのプログラムを提供することを目的とする。
【０００７】
この発明による構造物の設計支援方法（または連成系振動解析方法）は，振動を解析すべき領域の音場と構造のそれぞれの非連成振動方程式から正規化条件としてＫ直交条件を用いて算出された固有値と固有モードをコンピュータに入力してメモリに記憶し，メモリに記憶した固有値と固有モードの固有ペアについて，固有ペアにおける固有値と所与の基準固有値との差の絶対値が小さい順番にメモリ上で並び替えを行い，振動を解析すべき上記領域の音場と構造の連成系における固有値と固有モードを連成系の結合パラメータεの級数で表現した摂動法による展開式の係数をコンピュータにより算出するにあたって，上記係数を，離散化された連成系振動方程式に基づいて，連成系振動方程式内の固有値を周波数の逆数で置きかえた上で，連成系振動方程式を対称化することなく，非対称のまま，上記の並び替えされた固有ペアを用いて算出して，最終的に上記連成系の固有ペアを導出するものである。
【０００８】
この発明による構造物の設計支援システム（または連成系振動解析システム）は，振動を解析すべき領域の音場と構造のそれぞれの非連成振動方程式から正規化条件としてＫ直交条件を用いて算出された固有値と固有モードの入力を受入れてメモリに記憶する入力手段，上記メモリに記憶される固有値と固有モードの固有ペアについて，固有ペアにおける固有値と所与の基準固有値との差の絶対値が小さい順番にメモリ上で並び替えを行う並び替え手段，振動を解析すべき上記領域の音場と構造の連成系における固有値と固有モードを連成系の結合パラメータεの級数で表現した摂動法による展開式の係数をコンピュータにより算出するにあたって，上記係数を，離散化された連成系振動方程式に基づいて，連成系振動方程式内の固有値を周波数の逆数で置きかえた上で，連成系振動方程式を対称化することなく，非対称のまま，上記の並び替えされた固有ペアを用いて算出して，最終的に上記連成系の固有ペアを導出する処理手段，ならびに導出された上記固有ペアを出力する出力手段を備えているものである。
【０００９】
この発明による構造物の設計支援プログラム（または連成系振動解析プログラム）は，振動を解析すべき領域の音場と構造のそれぞれの非連成振動方程式から正規化条件としてＫ直交条件を用いて算出された固有値と固有モードをコンピュータに入力してメモリに記憶し，メモリに記憶した固有値と固有モードの固有ペアについて，固有ペアにおける固有値と所与の基準固有値との差の絶対値が小さい順番にメモリ上で並び替えを行い，振動を解析すべき上記領域の音場と構造の連成系における固有値と固有モードを連成系の結合パラメータεの級数で表現した摂動法による展開式の係数をコンピュータにより算出するにあたって，上記係数を，離散化された連成系振動方程式に基づいて，連成系振動方程式内の固有値を周波数の逆数で置きかえた上で，連成系振動方程式を対称化することなく，非対称のまま，上記の並び替えされた固有ペアを用いて算出して，最終的に上記連成系の固有ペアを導出するものである。
【００１０】
この発明においては，振動を解析すべき領域の音場と構造のそれぞれの非連成系の固有ペアを用いて，連成系の固有ペアを表現することを前提としている。音場と構造のそれぞれの非連成系の固有ペアを得るために，この発明では非連成振動方程式を正規化条件としてＫ直交条件を用いて算出している。このことは，摂動法による展開式の係数の算出における収束を保証するものとなる。
【００１１】
加えてこの発明では，上記の通り算出された固有ペアを所与の基準固有値（着目する固有値）との差の絶対値が小さい順番に並び替えを行っており，この並び替えの順序で固有ペアを上記係数の近似計算に用いている。上述のように展開式は収束することが保証されているので，固有ペアが多数存在したとしても上記近似計算を途中で打切ることができ，しかも上記の順序で固有ペアを用いているから，計算を途中で打切ったとしても得られる値は真の値に非常に近いものとなっている。
【００１２】
このようにして，この発明によると，現実的に実行が可能な時間内で振動解析の結果，すなわち連成系の固有ペアを得ることができる。
【００１３】
振動方程式内の固有値を周波数の逆数で置きかえているので，有限要素法により離散化された連成系振動方程式を非対称のまま取扱うことが可能である。
【００１４】
上記近似計算の途中打切りのために，許容誤差ｅと打切り回数ｍとをコンピュータに入力するようにするとよい。上記コンピュータによる上記係数の近似計算において，前回の計算結果と今回の計算結果との差が上記許容誤差ｅよりも小さくなったときに近似計算が停止される，または上記差が上記許容誤差ｅよりも大きくても，繰返し計算回数が上記打切り回数ｍに達した場合には，近似計算が停止される。
【００１５】
上記基準固有値は，上記コンピュータのメモリに記憶された固有ペア内の固有値において指定される。
【００１６】
【実施例】
図１は自動車の断面図であり，ボディ１の内側にエンジン２とシート３が設置されている状態を表わしている。
【００１７】
自動車の車室内の振動（音を含む）には，エンジンを音源とする音，自動車が停止しているときのこもり音等の人間が不快と感じる騒音（６０〜　３００Ｈｚ程度）がある。振動現象の解析の目的はこれらの騒音の固有周波数，固有周波数における音圧および変位を求めることである。これらを求めることにより，固有周波数が不快に感じる帯域内にあれば，固有周波数をシフトさせる，音圧レベル，変位レベルを下げる等の対策が可能となる。
【００１８】
上記連成振動現象を解くために，音場と構造を考える。音場は，解析対象の領域としてユーザーが定める空間である。また，構造は，弾性体，固体，構造体であって，この実施例ではシートなどが該当する。
【００１９】
一例として，運転者のいる空間Ａについて解析する。図２（ａ）　は図１の自動車を３次元でモデル化したものである。車体のシート３より前方の部分を実線で，後方部分を破線でそれぞれ示している。車体の一部は１辺がπである立方体で表わされており，これが構造である（シートに対応するＳ０を含む）。Ω０は立方体の内部を示しており，図１でいうと車室内の空間Ａを表わす音場の領域である。
【００２０】
説明を容易にし，且つ検証を可能にするために，図２（ｂ）　に示すような２次元のモデルを考える。
【００２１】
図２（ｂ）　ではΓおよびＳで囲まれたπ×πの正方形を音場の境界とし，正方形の内部Ωを音場の領域とし，音場の境界（ΓとＳ）を構造とする。Γがボディで，Ｓがシートに対応する。
【００２２】
シートＳは支持点Ｐ１および支持点Ｐ２で固定され，シートＳ自体は微小にたわむことは許されるものとする。Ｐ１，Ｐ２の両支持点ではｘ，ｙ方向の変位は許されないが，支持点はシートＳのたわみによって支持点を中心に回転することは許される。（この条件を単純支持という）
【００２３】
図２（ｃ）　は，シートＳが音場Ω内の音源から発せられる音圧Ｐによって微小変化した状態（たわみが生じた状態）を示すものである。ｕは，音圧Ｐにより実線で表わすシートＳが一点鎖線で表わすＳ’の状態になり，シートＳ上の点Ｐｘ１　がＳ’上の点Ｐｘ２に変位したときに生ずるｘ方向の変位量を表わすものである。
【００２４】
【数１】
2004094331.tif 000003
【００２５】
式（１）　は図２の（ｂ）　に示す音場と構造についての関係を表わしたものである。式（１）　に関連して，Ωは図２の（ｂ）　でいうところの音場の領域，Γは音場の境界であり，構造である（具体的にはボディ）。Ｓは構造であり，音場の境界でもある（具体的にはシート）。ＰはΩ内の音圧，ｕは構造の変位，ｃは音速，ρ０は気体の質量密度，Ｄ　は構造の剛性率，ρ１は構造の質量密度，ｍはＺ方向のフーリエ（Ｆｏｕｒｉｅｒ）モード数，ω　は固有振動数（固有値は１／ω２　で表わされる）である。音場の境界条件は，Γ上で音圧が０，構造の境界条件は単純支持とする。
【００２６】
式（１−１）　は空気中における固有振動数ωと音圧Ｐの関係を表わす音場振動方程式である。
【００２７】
式（１−２）　は音場の境界条件であり，境界Γ上で音圧が０（検証を容易にするために０としたが，実際の計算では，式（１−２）　に代えて剛壁条件式　（１−２）’を用いてもよい）であることを表わす。
【００２８】
式（１−３）　も音場の境界条件である。音圧Ｐによる境界Ｓの法線方向の変化率が，空気質量密度ρ０と，構造の変位ｕと，固有振動数の２乗ω２との積に等しいことを表わす。構造の変位ｕは，振動によって生ずる構造Ｓのｘ方向への変位を表わすものである。
【００２９】
式（１−４）　は図２（ｂ）　に示すような境界Ｓが，音圧Ｐを受けた時の構造の固有振動方程式である。
【００３０】
式（１−５）　は構造境界条件で，境界Ｓの両端にある支持点は単純支持であるということを表わしている。式（１−１）　から式（１−５）　をまとめて式（１）　とする。
【００３１】
式（１−３）　は構造が音場に与える影響を表わす式（音場系の式）であり，式（１−４）　は音場が構造に与える影響を表わす式（構造系の式）である。これら２つの式によって音場と構造が結びつけられている。
【００３２】
音場と構造の結びつきの強さを表わすパラメータをεとし，式（１−３）　と式（１−４）　の右辺にパラメータεを乗じたものをそれぞれ式（２−３）　，（２−４）　とする。式（１−１）　，（１−２）　，（１−５）　にそれぞれ対応する式（２−１）　，（２−２）　，（２−５）　と上記の式（２−３）　，（２−４）　の全体を式（２）　とする。
【００３３】
ε＝１の場合は完全連成を表わす。ε＝０の場合，音場と構造を関係づける因子はなくなり，完全非連成となり，この場合には，式（２）　の音場系の振動方程式，式（２−１）　〜式（２−３）　は，次の式（３）　で表わされ，構造系の振動方程式，式（２−４）　〜式（２−５）　は，次の式（４）　で表わされる。
【００３４】
【数２】
2004094331.tif 000004
【００３５】
【数３】
2004094331.tif 000005
【００３６】
式（２）　を摂動法の原理を用いて解く。まず，一般的な摂動法を説明する。
【００３７】
摂動法とは，ある複雑な系（系とは解析しようとする対象を指す。例えば，自動車の車室の騒音を解析するときの構造であり，音楽ホールの共鳴を解析するときのホールの構造である。）の状態を求める場合に，厳密に求めるのではなく，この系を構成する既知で簡単な系の状態から近似的に求める方法である。このとき，複雑な系は既知の系に小さな変化を与えることによって実現される状態であることを前提としている。摂動法によって求められるある系に微小変動を加えたときの固有値は，既知の系の固有値の級数で表わされる。
【００３８】
式（２）は理論解を求めるものである。実際の解析に応用するために，一般的には，連続系の式は離散化され，有限要素法を用いて解かれる。式（２）を離散化すると，式（５）が得られる。
【００３９】
【数４】
2004094331.tif 000006
【００４０】
ここで，式（５）について，次の式（６）を用いた置き換えを行うと，
【数５】
2004094331.tif 000007
式（７）のように書くことができる。式（７）は非対称のままである。上述したように，非対称のマトリックスを解くのは対称マトリックスを解くのに比べ，所要時間もコンピュータ上の一次記憶域も余分に必要となる。
【００４１】
【数６】
2004094331.tif 000008
【００４２】
最終的な目的は連成系における固有値と固有モードを求めることにあり，連成系における固有値と固有モードは非連成系の（音場と構造の）固有値と固有モードで表現することができる。この表現において摂動法を用い，式（８−１），（８−２）に示すように，固有値λｉ　と固有モードφｉ　をそれぞれ結合パラメータεの級数で表わす（摂動法による展開式）。
【００４３】
そこで，式（５）に示す非対称のマトリックスについて，以下に，摂動法を用いて近似解を求める。式（７）のλｉ（ε）とφｉ（ε）について摂動法を用いて近似解を求める式は，上記の通り式（８−１），（８−２）のようになる。これらの式における高次項の係数は，下記式（９）（式（９−１）〜式（９−４））で求められる。
【００４４】
【数７】
2004094331.tif 000009
【００４５】
【数８】
2004094331.tif 000010
【００４６】
式（８）（式（８−１）〜式（８−２））　において左辺は連成系における固有値λｉと固有モードφｉである。右辺のεの級数の各係数は非連成系における（音場と構造の）固有値と固有モードを用いて算出される。
【００４７】
非連成系における固有モードを求めるための正規化条件としてＫ直交条件を用いることにより（従来はＭ直交条件を用いていた），式（８）　が収束すること　（項数が増大すればするほど，各項の係数の値が小さくなること）が保証されることが分った。
【００４８】
式（８）が収束することを説明する。下記式（ａ）は，式（７）のεを「０」とした場合，即ち連成前の状態を示すものである。
λｉＫφｉ　＝　Ｍφｉ　…（ａ）
式（ａ）の左からφｉＴを乗ずると下記式（ｂ）のようになる。
λｉφｉＴＫφｉ　＝　φｉＴＭφｉ　…（ｂ）
【００４９】
ここでＫ直交とは次の式（１０）を満たすことである。ただしＫは剛性マトリクスである。
【００５０】
【数９】
2004094331.tif 000011
【００５１】
Ｋ直交条件を満たすということは，式（ｂ）の「φｉＴＫφｉ　」が「１」に等しいことを意味する。また，高次の固有値「λｉ」は「０」に近づくことが分っている。したがって式（ｂ）は「０」に収束することとなり，構造解析のための計算が収束に向かうことが分かる。
【００５２】
連成系の固有値，固有モードは，下記のとおり既知の非連成系の固有値λｉ（０），固有モードφｉ（０）を第１項とする級数で表わされる（ε＝１の場合）。
【数１０】
2004094331.tif 000012
（λｉ（０）　は連成前の構造または音場の固有値であり，奇数項は全て０になることが分かっている。φｉ（０）は連成前の構造または音場の固有モードである。）
【００５３】
上記式（ｃ）　の第２項以降（補正後）はλｉ（２ｎ）　（但しｎ＝１，２，３…）と表わすことができ，図３に示す処理フローの第ｎステップの近似計算により順次求られめる。
【００５４】
固有値と同様に，連成系の固有モードも既知の非連成系の固有モードを第１項とする級数で表わされる。φｉ（ｎ）　も同様に近似計算によって求められるものである。
【００５５】
固有値の近似計算をする際には，式（ｃ）と式（ｄ）の初項を決定しておく必要がある。初項とは，λｉ（０）　やφｉ（０）　である。自動車の車室内の騒音について解析するような場合は，人間の耳に不快に感じる周波数帯に着目する。したがって，連成系の固有値をコンピュータにより求める（計算する）に先立ち，音場の固有値と構造の固有値をコンピュータに読み込むが，音場の固有値と構造の固有値を比較して，より不快と感ずる周波数帯をもつ方を上記式の初項とする（これを基準固有値という）。
【００５６】
式（ｃ）　はλｉ（０）　を除く右辺の項を無限個求めれば等号が成り立つものである。しかし，現実には有限回で計算を打ち切る必要があるので，図３において計算される第ｎステップ終了後の固有値は式（９−２）　のように表わすことができる。計算によって求められる固有値は式（９−２）　のｎが大きいほど（式（ｃ）　の項をより多く求めるほど）真の値に近づく。固有モードについても同様である（式（９−４））。
【００５７】
図５（ａ）　は，図２（ｂ）　に示す音場の領域Ωを１辺がπ／４である３２個の正方形に分割し，各正方形の頂点の座標を示すものである。（横軸をｘ，縦軸をｙとし，座標は（ｘ，ｙ）で表わされる。また，π／４を１とする。）
【００５８】
一般的には，図５に小さな黒丸で示す２５個の点における音圧の集合が音場の固有モード（固有ベクトル）を表わすものとする。一方，構造については計算結果の精度を確保するために，音場よりも細かい分割をして固有モードを計算することが多い。境界Ｓ上において音場のデータを取得する点は○で囲まれた５点であるのに対して，構造のデータを求める点については構造Ｓ上で三角で表示される９点が設けられている。
【００５９】
図５（ｂ）　は，連成系のデータとして採用するデータＤ３と，境界Ｓ上における構造系のデータＤ１および音場系のデータＤ２の関係を図示したものである。構造系のデータはｕ（１）　からｕ（９）　まで９個のデータがあるが，連成系のデータとして採用するのはｕ（１）　，ｕ（３）　，ｕ（５）　，ｕ（７）　，ｕ（９）　の５個であり，音場系のデータは図５（ａ）　で黒丸で示す２５個の点のうち，連成系のデータとして採用するのはｐ（１，５）　，ｐ（２，５）　，ｐ（３，５）　，ｐ（４，５）　，ｐ（５，５）　の５個である。
【００６０】
上述の説明でわかるように，連成振動現象を解く場合に必要なデータは，音場と構造の接する部分（構造表面）だけでよく，さらに，構造のデータは音場のデータと重なる点のデータだけでよい。
【００６１】
図６（ａ）　は，上記の記述に基づいて，図３のフローチャートの処理１０で読み込むデータの並びを示したものである（このデータはコンピュータのメモリに記憶される）。音場の固有値をλｓ＿ｎ　で表し，データＤ３のように選択された固有モードを音圧ｐ（ｉ，ｊ）＿ｎ　で表す。構造の固有値はλｐ＿ｎ　で表し，固有モードを変位ｕ（ｉ）＿ｎ　で表す。ｎは入力するｎ個の固有値の順番を表わす。図６（ｂ）　は具体的数値例を示す。
【００６２】
音場の固有値λｓ＿ｎ　と音場の固有モードｐ（ｉ，ｊ）＿ｎ　は式（３）　によってあらかじめ求められ，構造の固有値λｐ＿ｎ　と構造の固有モードｕ（ｉ）＿ｎ　は式（４）　によってあらかじめ求められる。このとき，上述したように正規化条件としてＫ直交条件が用いられる。
【００６３】
図６（ａ）により具体的に説明すると，音場の固有値はきわめて多数存在するが，簡単のために５つ（λｓ＿１，λｓ＿２，λｓ＿３，λｓ＿４，λｓ＿５）が示されている。固有値λｓ＿１とペアをなす固有モードがｐ（１，５）＿１，ｐ（２，５）＿１，　・・・・・，ｐ（５，５）＿１であり，これらは図５（ｂ）　の縦に配列された５つの点の音圧によって表わされる固有モードである。
【００６４】
図６（ｂ）は図６（ａ）のそれぞれについて具体的数値例を示すものである。
【００６５】
図３は実施例の処理フロー図である。この処理はコンピュータによって実行される。処理１０では予め計算しておいた構造系と音場系それぞれの固有値と図５のＤ３で示すように選択された固有モードをコンピュータに，またはメモリのワークエリアにより読み込む（図６（ａ），（ｂ）参照）。図６（ａ），（ｂ）において，上段が音場系の固有ペアであり，下段が構造系の固有ペアである。固有ペアの値がＦＤ等に記憶されている場合には，ＦＤからワークエリアに読込めばよいし，ハードディスクに記憶されている場合には，それをメモリのワークエリアに読込めばよい。固有ペアの値を他のコンピュータから伝送してもよい。
【００６６】
処理２０では計算の基準とした（構造か音場のいずれか）固有値λ（基準固有値）と，許容誤差ｅおよび繰返し計算打ち切り回数ｍとをコンピュータに入力する。この入力は一般には手操作により行われるが，ＦＤから入力，またはオンラインで入力してもよい。固有値λは周波数ωとの間にλ＝１／ω２　の関係がある。許容誤差とは近似計算の結果，ｎ−１回目の計算結果の固有値とｎ回目の計算結果の固有値との差が指定した許容誤差ｅよりも小さければ計算を打ち切らせるための値である。繰返し計算打ち切り回数ｍは，収束計算が期待する許容誤差範囲に無い場合であっても（すなわち，上記のように，ｎ−１回目の値とｎ回目の値との差が許容誤差ｅより小さくなくても）計算を中止させるためのものである。なお，処理における変数ｎは，上記のｎ個の固有値の順番を表わすｎとは意味が異なるものである。
【００６７】
処理３０では処理１０で入力した複数の固有ペアを固有値をキーにしてメモリ（ワークエリア）上で並べ替えを行う。並べ替えは各固有ペアの固有値と処理２０で読み込んだ基準固有値λとの差の絶対値が小さい順番に行う。並べ替えは構造系と音場系それぞれについて行う。その後，音場と構造の共通節点（図５で□を付けた部分）の固有モードを選び出し，メモリに保存する。（これにより，保存量を減らすことができる。図５の場合では，音場の全体節点数２５が５になり，構造の全体節点数９が５になる。）
【００６８】
一般に自動車の車室内の騒音について解析するような場合は，人間にとって不快な騒音と感ずる周波数の音を中心に分析を行う。図６の例で，構造固有値λｓ＿４　が不快な音になる周波数の固有値（基準固有値）であるとした場合，処理３０ではまず第１番目にλｓ＿４　（これが基準固有値である）を置き，残りの固有値についてλｓ＿４　に近い順番（各固有値とλｓ＿４　との差の絶対値が小さい順）メモリ上でに並べ替える。
【００６９】
処理４０ではカウンタを＝０とする。処理５０ではカウンタを１プラス（インクレメント）する。処理６０では式（９）での近似計算λｉ（２ｎ）　を上記において並べ替えした順序で固有ペアを用いて，行う。式（９）　は式（８）　における各係数を求めるための計算式である。式（９）において，λｉ（２ｎ）を計算するためには（式（９−２）），λｉとφｉを用いる。φｉは式（９−３）または式（９−４）を用いて計算する。このときλｉとλｊ（ｉ≠ｊ）を用いる。λｊのｊは，ｉを固定したときにカウンタによりカウントアップされていく（下記のｎに相当）。このカウンタにより指定される固有値の順番が上記の並べ替えの順番である。
【００７０】
処理７０では処理６０で求めたλｉ（２ｎ）の絶対値を許容誤差ｅと比較する（λｉ（２ｎ）はｎ回目の計算結果の固有値とｎ−１回目の計算結果の固有値との差である。λｉ（２ｎ）　が許容誤差ｅより大きければ（処理７０でＮＯ），処理７２に進み，ｎ（カウンタの値）がｍに達したかどうかをチェックする。ｎがｍに達していなければ（処理７２でＮＯ），処理５０に戻り，カウンタの値をインクレメントして再び近似計算を行う（処理６０）。
λｉ（２ｎ）　が許容誤差ｅより小さければ処理８０に進む。λｉ（２ｎ）　が許容誤差ｅより小さくならなくても，ｍがｎに一致した場合（ｎ＝ｍになったとき）には（処理７２でＹＥＳ），処理７４　に進み，回数制限で計算打切られたことを示すメッセージを出力し，処理８０に進む。近似計算は途中で打切られるので計算時間は短くてすむが，下に説明するように収束が確保される。
【００７１】
処理８０では連成系の近似固有値λ＝λｉ＋λｉ（２）＋λｉ（４）＋…λｉ（２ｎ）を計算して出力する。近似固有モードも同様に求められるのはいうまでもない。
【００７２】
上記の処理１０〜８０はパラメータεを適当な値に固定して実行される。もし，所望の，または適切な結果が得られない場合には，パラメータεを変更して再度上記の処理を繰返すことが望ましい。これは，たとえば構造の材質を変更するようなことに相当するので，最も適切な材質の選定にも役立つ。
【００７３】
図４（ａ），（ｂ）は式（８−１）　の第２項以降の項のうちの１つの項について，縦軸にλ，横軸に計算回数を取ったときの収束状況を示す図である。固有値をλに近い順に並べ替えずに計算した場合は，図４（ｂ）の曲線（ｂ）で示すように誤差が大小に変動し，何回で収束したと見なすかの見極めができない。一方，並べ替えを行った場合は，図４（ａ）　の曲線（ａ１）または（ａ２），（ａ３）のように確実に回を追うごとに収束に向かう形となる。
【００７４】
以下に具体例を示す。
【００７５】
第１の例　（構造からの連成系固有値）
許容誤差ｅ＝１０−６とする。
ｉ）現実的なパラメータを使うケース
Ｄ＝１９．２３０７６９，ρ０＝１．２９３００００，ρ１＝７８５０，ｃ＝３４０，ｍ＝１，ｉ＝１
理論解λｉ（１）＝１０２．０６１８８，λｉ（０）＝１０２．０５０００
ここで，理論解λｉ（１）は下記式（１１）のタイプ１で計算されたものである。
理論解λｉ（０）は式（３）で計算されたものである。
１回近似結果＝λｉ（０）＋λｉ（２）＝１０２．０６１８８，
｜λｉ（２）｜＝｜０．０１１８８｜＞ｅ
２回近似結果＝λｉ（０）＋λｉ（２）＋λｉ（４）＝１０２．０６１８８，
｜λｉ（４）｜＝｜０．０００００００｜＜ｅ
３回近似結果＝λｉ（０）＋λｉ（２）＋λｉ（４）＋λｉ（６）＝１０２．０６１８８，
｜λｉ（６）｜＝｜０．０００００００｜＜ｅ
４回近似結果＝λｉ（０）＋λｉ（２）＋λｉ（４）＋λｉ（６）＋λｉ（８）＝１０２．０６１８８，
｜λｉ（８）｜＝｜０．０００００００｜＜ｅ
【００７６】
ｉｉ）非現実的なパラメータを使うケース
理論解λｉ（１）＝６．６９４４７６５７５６，λｉ（０）＝６．２５
Ｄ＝２，ρ０＝５，ρ１＝５０，ｃ＝２．５，ｍ＝１，ｉ＝１
ここで，理論解λｉ（１）は式（１１）のタイプ１で計算されたものである。
理論解λｉ（０）は式（３）で計算されたものである。
１回近似結果＝λｉ（０）＋λｉ（２）＝６．６９４６６７３３４３６１３７７，
｜λｉ（２）｜＝｜０．４４４６６７３３｜＞ｅ
２回近似結果＝λｉ（０）＋λｉ（２）＋λｉ（４）＝６．６９４４６２７６４９８２１８７，
｜λｉ（４）｜＝｜−０．０００２０４５６９｜＞ｅ
３回近似結果＝λｉ（０）＋λｉ（２）＋λｉ（４）＋λｉ（６）＝６．６９４４７７５５４０００３７１，
｜λｉ（６）｜＝｜０．００００１４７８９｜＞ｅ
４回近似結果
＝λｉ（０）＋λｉ（２）＋λｉ（４）＋λｉ（６）＋λｉ（８）＝６．６９４４７６４７８０８１４９８，
｜λｉ（８）｜＝｜０．０００００１０７５９｜≒ｅ
【００７７】
第２の例　（音場からの連成系固有値）
許容誤差ｅ＝１０−４とする。
理論解λｉ（１）＝０．０６９２５１５９９，λｉ（０）＝０．０７１１１１１１
ここで，理論解λｉ（１）は式（１１）のタイプ２で計算されたものである。
理論解λｉ（０）は式（２）で計算されたものである。
１回近似結果＝λｉ（０）＋λｉ（２）＝６．９０７５９２２３２４４８７５０７×１０−２，
｜λｉ（２）｜＝｜−２．０３５１８８×１０−３｜＞ｅ，
２回近似結果＝λｉ（０）＋λｉ（２）＋λｉ（４）＝６．９２６５８９４６１４１７９７４２×１０−２，
｜λｉ（４）｜＝｜−４．９３３２９×１０−４｜≒ｅ
【００７８】
【数１１】
2004094331.tif 000013
【００７９】
図７は構造物の設計支援システム（または連成系振動解析装置）の構成側を示すものであり，コンピュータシステムにより実現される。
【００８０】
処理装置（並び替え手段，処理手段）１はコンピュータ本体であり，図３ステップ１０〜ステップ８０の処理を実行するプログラム１ａを格納するとともに，図６に示すようなデータを格納するメモリ１ｂを備えている。
【００８１】
入力装置（入力手段）２はキーボード，マウス，ＦＤドライブ，通信装置等により実現され，ステップ１０，２０におけるデータ等の入力のために用いられる。
【００８２】
出力装置（出力手段）３は表示装置，プリンタ，ＦＤドライブ等により実現され，ステップ８０の固有値出力に用いられる。
【図面の簡単な説明】
【図１】自動車の断面図であり，この発明に適する現実問題の一例を示すものである。
【図２】図１をモデル化した図であり，（ａ）は３次元モデルを，（ｂ）は２次元モデルをそれぞれ示す。
【図３】コンピュータにおける処理手順を示す処理フロー図である。
【図４】演算の収束状況を表わすグラフであり，（ａ）は収束する場合を，（ｂ）は収束しない場合をそれぞれ示す。
【図５】（ａ），（ｂ）は音場系と構造系の共通節点を表わすイメージ図である。
【図６】（ａ），（ｂ）は音場系と構造系の固有ペアのメモリ上の保存状態を表わす図である。
【図７】構造物の設計支援システムまたは連成系振動解析システムの構成例を示すブロック図である。
【符号の説明】
１　　処理装置
１ａ　　　プログラム
１ｂ　　　メモリ
２　　入力装置
３　　出力装置
【図１】
2004094331.tif 000014
【図２】
2004094331.tif 000015
【図３】
2004094331.tif 000016
【図４】
2004094331.tif 000017
【図５】
2004094331.tif 000018
【図６】
2004094331.tif 000019
【図７】
2004094331.tif 000020