最新版を正確に読む場合には pdf ファイル をどうぞ。これは web 検索のための簡易旧版です。
ここまでは,離散的に分布する質点がバネやダッシュポットで結合されている 系を対象としてきた。しかし,身の回りに,あるいは社会基盤構造には, 質量が離散的に分布しているそういった系はほとんど見当たらない。 実際の構造は,連続的に分布した質量を持つ連続体として捉えなければならないことは 明らかであろう。ただ,振動モードといった概念は,連続的な分布質量系でも 非常に重要な概念である。ここでは,拡がりをもった系の振動解析の基本として, まず曲げ剛性の無い系の振動を説明する。
吊橋や斜張橋に用いられているケーブルは弦としてモデル化できそうだ。
弦は,予め大きな張力で引張られて初めて安定する構造である。
ここでも重力の影響は無視し(つり合っている状態からの変動のみを
対象とし),支配方程式を誘導しておく。
まっすぐな弦が振動して変位したある瞬間の図-10.36で,
上下方向の運動方程式は
となる。は弦の密度で,は弦の断面積である。
弦の変位が十分小さい場合には,幾何学的な関係から
となる。これを上式に代入して整理すると
となる。あるいは
となる。ここには波の速度(位相速度)である。 あるいは式(10.130)は波動方程式 とも呼ばれる。
ここでは長さの両端を固定された弦の振動だけを対象とするので,
境界条件は両端で
と与えられることになる。初期条件については,
時間についても2階の微分方程式なので,初期形状と初速で
のように与えればいい。
ここでは変数分離で解を求めてみよう。つまり
と置き,運動方程式(10.130)に代入して整理すると
となる。プライムはに関する微分を表す。この式は
任意の場所と時刻で成立しなければならないが,左辺はのみの
関数であり,右辺はのみの関数になっており,この等号が成立するのは
この両辺ともに定数であるときのみである。つまり
と置けばいいので
が独立して成立しなければならない。境界条件も変数分離すると,
任意の時刻でそれが成立しなければならないので
となる。初期条件は分離できないので,しばらく放置しておこう。
運動方程式と境界条件だけが分離できたので,
式(10.133)のの方が
先に解けそうだ。まずの場合はと
置くと,一般解は
となるので,境界条件式(10.134)に代入すると
から積分定数を決めることになる。しかしで, そのときは常に正なので, 結局, のみが解となり,意味のあるは存在しないことになる。 では次にの場合を考えると,一般解は となるので, 境界条件が, を要求し,これも意味のある解にはならない。 したがっては負でなければならないことがわかる。
では
と置くと,一般解は
となり,境界条件から
を得る。ここでが零にならないためには
でなければならない。したがっては不定になり,解は
と求められる。最終的に式(10.133)のの
微分方程式に,式(10.135)のを代入して解くと
が解である。この形が固有振動数のsine, cosine関数に
なることから,弦の振動の固有振動数は
となる。つまり,長い弦や重く太い弦は振動数が小さくなるので,音は 低くなる。一方,その逆にするか,張力を大きくすると,音は高くなる。 グランドピアノの蓋を開けてみて欲しい。 低い音のピアノ線は太く長くなっているし, 高い音のピアノ線は細く短くなっている。 調律(チューニング)では弦の張力を調整する。
もう一度を求めるプロセスを書くと,解くべき問題は
と表すことができる。これは
が解の一つであり,
微分方程式で表された固有値問題になっている。は固有関数
と呼ばれ,あるいはが固有値である。
物理的な言葉を使うなら,が固有振動モード
であり,
が固有振動数である。連続体の場合は,
式(10.136)のように無限個の固有関数があり,
そのすべてが解の候補である。固有値もそれに対応して無限個ある。
ここで多自由度系の振動を解くときのモード解析法を思い出してみよう。
モード解析法では,まず自由振動解析をして
固有振動モード(固有ベクトル)と固有振動数(固有値)を
求め,そのモードを用いて解を級数解で仮定し,答を求めようとした。
弦の場合も同様に,上式(10.137)の固有値問題を解いて
固有関数と固有値を求めることができた。したがって,多自由度系の解の
表現式(10.107)を拡張し
と仮定して,問題を解くことを考えてみる。
まず多自由度系のモード解析法で重要な特性であった
固有振動モード(ベクトル)の直交性を思い出して,
弦の固有振動モード(関数)の直交性について確認しよう。
ある次の固有関数は,式(10.137)から
を満足している。
モードの直交性を検証するときには,別のモードを仮想変位と
考えて仮想仕事を算定してみることから始めていた。ここでも
同じアプローチをしてみる。つまり上式にを乗じて,
全領域積分すれば,弦全体の仮想仕事になる。
つまり,その仮想仕事式は
となる。第1項を部分積分すると
となるから,境界条件式(10.137)で両端の項が消え,結局上式は
という仮想仕事式が成立する。
同じプロセスを第次モードの方程式に適用するために,
固有関数をそれに乗じて仮想仕事式を求めると
になる。この2式を辺々差し引くと,第1項は同じなので消えてしまい
となる。式(10.137)から明らかなように,異なる固有値同士は
異なる値になるので,この問題の固有関数は
という性質を持っている。この左辺を二つの関数の内積 と呼んでいる。そして,二つの異なる関数の内積が零になることから, その関数同士は直交していると称するのである。 このような直交性 は,モード解析法にとって最も重要な特性である。
では,式(10.138)の級数解を運動方程式(10.130)に
代入すると
となるので,式(10.139)を左辺に代入すると
が成立する。この式と第次の固有関数と内積(仮想仕事)をとると,
結局直交性の式(10.140)によって,級数和の第次成分
以外はすべて零になり
が,係数を決定する方程式になる。この一般解は
と求められる。したがって,元の問題の一般解は
となる。
最後に,この一般解の式()を初期条件式(10.132)に
代入すると
となる。再度,第次の固有関数との内積をとると
となることから,積分定数が
と求められる。
文献[62]にある例題を転記しておく。
弦の中央をだけつまみ上げ,そぉっと手を離したときの
振動を求めてみる。初期条件は
のように書くことができる。これを式()に
代入すれば,積分定数を求めることができ,結局
となる。これを図-10.37に示した。 破線は級数の40項を用いた解で,実線は5項のみで描いた解である。 少ない項数でも比較的いい解が求められている。 厳密解は次の節の図-10.38に示されるように, 三角形や台形のような区分的な直線である。
もう一度,一般解の式()を眺めてみよう。 これを少し書き直すと
となることから,解はという関数とという関数で
のように表されていることがわかる。
ではという関数はどういった現象を表しているだろうか。
例えばある時刻にある場所で,この関数がという
値を持っていたとする。それから時間が
経った
のときの,
という場所の
この関数の値は
となる。 つまり,, における「物理状態」がそのままにで 観察されたことになる。したがって,という関数は, 「の正の方向にで伝播する」物理現象を表している。 同様に,は, 「の負の方向にで伝播する」物理現象を表す。 したがって弦の振動解は,正の方向にで伝播する波と, 負の方向にで伝播する波の重ね合わせと解釈してもいいことに なる。図-10.38には,前節で扱った例題の, 中央をつまみ上げて離した場合の,波動の解の式()の模式図を示した。 右に速度で伝播する実線と左に伝播する破線の重ね合わせが下にある図の 解であり,初期の三角形が台形になり,その高さが徐々に変化していくのが 正解である。 このように,波動と振動は同じ現象であり,その観察の仕方を変えただけで あることがわかる。 したがって,式(10.130)は波動方程式とも呼ばれているのだ。
最後に,モード解析法という観点から,多自由度系と弦を比較しておく。
まず固有値問題は
と表される。固有ベクトルと固有関数は,次のような内積に対して
直交性
を持っていることから,直交するベクトル,あるいは直交する
関数の集合である。そして解を
と仮定して問題を解こうとしている。
連続体は無数の自由度を持つと考えられることから,級数が無限級数になっている。
そして一般化された座標と呼んでもいいは,それぞれ
を解くことによって求められる。これは見かけ上,1自由度系の 運動方程式である。最終的に,積分定数は初期条件から 決定することができ,その際にも重要な性質が直交条件である。
膜というのは弦が2次元的に拡がったようなもので,これも比較的
大きな引張り張力を与えておかないと,安定な構造にはならない。
ここでは,あまり社会基盤構造では使わないが,ドームの天井を
設計するための頭の体操に
するためと,いわゆるFourier級数は必ずしも三角関数による
級数だけではないことを示すために,円形膜の振動を考えてみる。
その運動方程式は,式(10.130)を2次元に拡張した
になることは,すぐに推測できると思う。は
やはり張力と密度で定義された位相速度である。ただここでは,
円形領域の問題を対象としていることから,直角座標ではなく,
極座標で表現した方がいいことも容易に推測できる。
極座標でのラプラシアンは
となることは,数学の本を見ればわかるので,結局,
膜の運動方程式
は
となる。は膜の外径である。境界条件は,周囲が固定されているものとして
となるが,運動方程式(10.142)がで特異になる(分母が
零になる)ことを考慮して,円形膜の中央で「おかしなこと」が起きないことを
条件にするために
を付帯的な条件としておく。ただし,
もし内径がのドーナツ状の膜を対象とする場合には,
式(10.144)の付帯条件の代わりに
と
なることには注意して欲しい。さらに,膜は折れ曲がったり破れたりしないで
連続したままであるから,例えばにおける
という連続条件も課すことにする。初期条件は,時刻における
初期形状と初速を与えればいいから
と与えられるものとする。
弦の場合と同様に変数分離をすれば求めることができるが,
弦の固有値問題の拡張から,最も重要な固有値問題は
となりそうなことは,きっと容易に推測できると思う。
が
固有振動モード
である。疑問を持った人は変数分離をしてみて欲しい。が定数である。
ここでさらに変数分離をして
を上式に代入して整理すると
が成立しなければならないことがわかる。
プライムは,それぞれ引数についての微分である。
つまりはによる微分で,
はによる微分である。
これから,それぞれについての微分方程式は
と表現できる。一方,境界条件式(10.143) (10.144)
(10.145)は,については
であり,については
となる。
まず
は,微分方程式(10.149)を眺めた上で,
境界条件式(10.151)が周期解であることを要求するので,
最低限,この定数は正でなければならないことが明らかである。
また,その解は周期でない限り
境界条件(方向の連続条件)式(10.145)を満足しないので,
結局は非負の整数
でなければならず,対応する固有関数は
となる。を含むことに注意する。
このをの方程式(10.149)に代入すると
がを決める微分方程式になる。
まずの符号を検討するために,
まずそれが正()であると仮定してみる。
そこで式(10.154)をで
割ったものにを乗じて,からまで積分する。つまり
が成り立つ。右辺の被積分関数は非負なので,積分値も非負になる。
一方左辺については,部分積分をして境界条件を考慮すると
となる。結局,式()の左辺は負で,右辺は非負になることから,
仮定の()が成立しないことがわかる。
よって
と置こう。
こうすると,についての微分方程式は,式(10.154)から
となる。この式はBesselの微分方程式と呼ばれ,解はBessel関数
で
と表現することが多い。は次の第1種Bessel関数,は次の 第2種Bessel関数と呼ばれる。それぞれの例を図-10.39に 示した。
しかし,境界条件式(10.150)は,における有界性を
要求していることから,図からも明らかなように,
第2種のBessel関数は解の候補にはならない。
もちろん,ドーナツ状の膜の場合はを含まないので,第2種のBessel関数も
候補にはなる。したがって,ここの円形膜の場合には,は
が解になる。さらに,周上の境界条件式(10.150)から
が成立しなければならない。この条件が
固有値を決定する。図-10.39から
容易に推測できるように,Bessel関数は,一定周期ではないものの
周期的に零をまたぐ関数で,零点は無限個存在する。したがって,
上式()を満足するはそれぞれのに対して無限個
存在する。したがって,それを1個目から, 2, と記すことに
して,という固有値が
を満足しているものとして求められる。表-10.1に具体的な
零点を列挙したが,これを用いると例えば
のような値をは持つことになる。したがって,
固有関数の解は最終的に
となる。Bessel関数が一定周期ではないことから,も例えばの 整数倍のような表現にはならないことに注意して欲しい。
まず式(10.153)の固有関数
は,前節の弦の
場合の特性からも容易に予想できるように,に対して直交性
が成立する。すなわち直交関数列
である。
一方,Bessel関数の場合はどうなっているだろう。直交性の検討の場合には,
つり合い式(微分方程式)にある一つの固有関数を仮想変位と想定したときの
仮想仕事式を用いてきた。ここでも同様のアプローチを試みる。
あるBessel関数は,式(10.149)あるいは
式(10.155)から
を満足している。これにを仮想変位として乗じて
円形膜内で積分する。極座標の微分面積が
で
あることに注意すれば,その仮想仕事式は
としなければならないことがわかる。
被積分関数の第1項は,部分積分して境界条件を考慮すると
となるので,上の仮想仕事式は
となる。同じプロセスをが満足する微分方程式に適用するために,
それにを乗じて仮想仕事式を誘導すると
となる。この2式を差し引くと,結局第1, 2項は同じなので消えてしまい
となる。異なる固有値同士は値が違うから,最終的にのとき,
次式で定義した内積で
という関係が成立する。すなわち,Bessel関数も直交性 を有している。 弦の内積の定義式(10.140)とは若干異なり,`'という 重みつき内積になっていることに注意する。
最初の固有値問題に戻ると,式(10.148)の固有関数が
最終的に
のいずれかとして求められたことになる。
この二つの固有関数は,式(10.159)と
式(10.160)の直交性から,次の内積のもとで
のような直交性を有する。
そしてこの固有関数
は,式(10.147)から
を満足する。 簡単のため,しばらくはの上付き添え字SとCを省略する。
いよいよ元の問題に戻ろう。求められた固有関数を用いて,解を
と仮定する。これを膜の運動方程式(10.142)に代入すると
となる。左辺の被積分関数の中括弧に式(10.163)を代入すると
となる。この両辺にを乗じて内積(仮想仕事式)を求め,
式(10.162)の直交条件を考慮すると,
結局無限級数の中の, の項のみが非零になり
が,に対する微分方程式であり,解は
となる。したがって,元の問題の一般解は
となる。
最終的に,未定係数は初期条件で決定される。式(10.146)に 式(10.167)を代入すると
となる。そこで
あるいは
との
内積をとり,式(10.162)の直交条件を使えば
と求めることができる。図-10.40には各振動モードを 描いた。上半分の曲線が方向ののモードであり, 一点鎖線より外側にある部分ではたわみが正で,内側にある箇所では 負になっていると見て欲しい。 つまり,この図の並びの2行目は, までの部分が 下にたわむときに, の箇所は上に変位することを 示している。また,右軸方向に描いたのが方向のBessel関数で 表されたモードであり,この図の下向きが正のたわみと見て欲しい。 例えば2列目の図では,膜の約半分が下にたわむときに 残りの部分は上に変位することを示している。
になる。
になる。