1.=jωρ{Aej(ωtkx)ωt+kx.pdfVIP

7-185 土木学会第60回年次学術講演会（平成17年9月） セルオートマトン法を用いた 1 次元音響管の解析 長岡技術科学大学 学生会員 ○小野 泰和 長岡技術科学大学 学生会員 富樫 孝介 長岡技術科学大学 環境?建設系 正会員 宮木 康幸 1.はじめに 平面波の場合,一般的に次式で求まる． 騒音伝播解析には，有限要素法や差分法など，現 p(x,t) = jωρ{Ae j(ωt?kx) Be j(ωt+kx) } (2) 象を支配方程式によりモデル化し解析する手法が主 ここで，A,B は境界条件より定まる係数であり，ω に用いられてきた．しかし，音の伝播は複雑であり， は音源の角周波数，k は波数(k =ω/c)である．一端に 結果を得ることは難しい．そこで本研究では，セル 加振端を持つ音響管の音圧に関する解は，境界条件 として粒子速度 jωt を与え，さらに他端が オートマトン法(Cellular Automata：以下CA) を用い u(0,t) = u0e て解析を行うことにした．本モデルの妥当性の検討 開放端の場合， p(l,t) = 0 の条件より， に際しては，1 次元に関して，1次元音響管を用いて sin k(l ? x) jωt (3) p(x,t) = jρcu0 e 開管と閉管の2つの場合について，波動方程式の解と cos kl CA による解析解を比較する． と求まる．ここでｌは管の長さを表す．また，閉端 2.セルオートマトン法 の場合も同様に，境界条件u(l,t) = 0 を考慮して， CAは解析対象をセルと称する均一な大きさの区 cos k(l ? x) jωt (4) p(x,t) = ? jρcu0 e 分領域に分割し，各セル上に定義された離散的な状 sin kl 態量をごく近傍のセル間に設けた局所近傍則によっ