最新版を正確に読む場合には pdf ファイル をどうぞ。これは web 検索のための簡易旧版です。
まず,3次元の連続体のつり合い式は式(3.22)であった。
したがって,連続体の運動方程式は,このつり合った
力の総和がNewtonの法則を満足する条件として得ることができる。
ここ以下では,p.で
規定した総和規約を用いて運動方程式を表現すると
となり,これは式(3.96)と同じである。
添え字にあるコンマは,その次の添え字に相当する独立変数で微分することを
意味する。
ここには体積力であるが,つり合っている弾性体への
付加的な作用によって波が生じる場合を対象とするので,以下では無視する。
またHookeの法則は式(3.42)に示したように
と書くことができる。
ここに
はひずみテンソルであり,とはLaméの
定数である。またはKroneckerのデルタである。
さらにひずみは式(3.6)で定義したように,変位と
のように関係付けられている。式(11.14)を式(11.13)に
代入したあと,それを式(11.12)に代入して整理すると,
変位で表した運動方程式は
となる。あるいは
のような表現も可能である。 各自がわかり易い方を使うといい。節-3.5.1では2種類の ポテンシャルを用いて波動方程式を誘導したが,ここでは, もう少し物理的なアプローチを用いて最も基本的な波である 平面的に伝わる波の波動方程式を誘導し,それを解くことにしよう。
ある観測地から非常に遠い活断層で地震が発生したとすると,
その震源は観測地から見たとき,ほとんど点源にしか見えないだろう。
その点源からは球面状の波動が発生するのではあるが,
観測地が震源からとても遠いため,そこに到達する波は点源からの
曲率半径がとても大きな球面波であり,それは
ほぼ平面状の現象が伝播する波として近似しても,
工学的には問題は生じないであろう。さて,
式(11.5)で示したように,1次元の進行波はという
関数形で表現できる。この関数は,
の面の物理現象が,
その面の法線方向つまり方向にの速さで伝播していることを
示している。
一方3次元空間の任意の平面は,図-11.7に
示したように
で定義できる。
つまり,をある平面の単位法線ベクトルとすると,
その平面は
と表すことができる。
したがって,1次元の進行波からの類推で,3次元空間の平面波の場合は
という関数で,変位を表現できることは容易に推測できるだろう。 ここに,は変位ベクトルの向きを表す単位ベクトルであり,は 考えている波の位相速度である。 を 位相と呼ぶこともある。
式(11.16)のような平面状の波が存在するかどうかを
調べるために,それを式(11.15)の運動方程式に代入してみよう。
そのために式(11.16)を微分すると
等という関係が成り立つ。第2式ではが単位ベクトルで
あることを用いた。これを運動方程式に代入して整理すると
となるので,この関数形を持った解が存在するためには,係数が零で
なければならないから
が成立していなければならない。
この式(11.17)が成り立つ可能性は二つある。
となる。つまり,位相速度が式(3.100a)で定義したものと同じく
を満たす運動,つまり非回転波 になっている。 これはせん断変形が伝播するのではなく,1次元の波動のような疎密波に相当する。
となるから,位相速度は式(3.100b)で定義された横波(S波)
の速度になる。この場合は
以下では,-平面が水平面と平行で, 波の進行方向が-平面内にある場合を対象とする。 その条件下で,図-11.8に示したような, 次の3種類の波の挙動を検討する。
これまではの関数形は任意として扱い,
初期条件や境界条件で決まるものだったが,
ここからは定常運動を対象として,定周期運動とする。
つまり振幅の調和振動運動として
とする。ここに は 波数であり,は振動数である。 また,波が伝播する領域は図-11.9のように, 水平面で上下に別々の領域に分割されているものを対象とする。
SH波は,図-11.9に示した
の波である。
後述のような境界条件の検討から明らかになることではあるが,
このSH波の反射ではSH波しか生じない。したがって,
入射する波も反射波も,生じる変位は方向成分のみである。
そこで図のように角度と波の進むベクトルを定義すると
となるので,入射波は
であり,反射波は
となる。, は複素数の振幅である。
例えば比較的硬いものの下で生じる反射を
モデル化したものとして,より上が剛な場合をまず考えよう。
この場合,の境界条件はであるから,入射波と反射波の
重ね合わせから
が任意の, で成立しなければならない。
つまり,二つの指数関数()の形は共に同じでなければならず,の
係数同士との係数同士から
が成り立つ必要がある。よって
となるから,境界条件式(11.24)に代入し直せば
と求められる。したがって,定常状態にある波の方向変位は,
この結果を反射波の式(11.23)に代入したものと
入射波の式(11.22)とを重ね合わせて
となる。解の最初の関数 の 部分は,深さ方向の方向にはsine関数状の振幅分布をすることを 示しており,このような波を重複波 と呼んでいる。また二つ目の指数関数は, その引数が であることから,の正の 方向に速さ で進む進行波 であることを示している。
次に,例えば水平面より上が大気(真空)であるような
場合の反射をモデル化したものとして,が自由表面である場合を考えよう。
このときのの
境界条件は (, 2, 3)である。SH波は,
,
なので,各ひずみ成分は
,
,
,
となるので,Hookeの
法則からとは常に成立する。
したがって残った条件は
となる。このことから自由表面の境界条件は
で与えられる。剛体の場合と同様に,入射波と反射波の重ね合わせで
境界条件式(11.28)を表すと
が,任意の, で成り立たなければならない。したがって
が解になる。よって,定常状態の変位は,剛体の場合とほぼ同様に
と求められる。境界条件の違いは,重複波の部分のみに現れている。
SH波は,境界条件によらずSH波の反射しか生じなかったが,P波の場合には, 例えば
のような反射が生じることがわかっている。
ここでは代表的な例として,
自由表面が存在する場合の半無限体を対象として解くので,
その他の場合については各自試してみて欲しい。図-11.10に
示したように角度とベクトルを定義すると,
三つの波の変位を, 1, 2で区別して
と表すことができる。ここに
である。ひずみ成分はそれぞれ
と書くことができるので,Hookeの法則に代入して,各応力成分も
となる。
自由表面の境界条件はで, なので,
式(11.35)の重ね合わせから
となる。ここに
と置いた。
この条件式に式(11.33)を代入して整理すると
となる。この式が任意の, で成立するためには,まず指数関数同士が
等価でなければならないから
の条件が成り立つ必要がある。後者はSnellの法則
と呼ばれる。したがって
となり,角度は
と求められる。が1より大きいことから,が必ず
実数として存在し,最後の式で示したように,その値はより
小さいことには注意する。この結果を式(11.37)に代入して整理すると
となる。第1式をで除して
を使うと
と書き換えられ,さらに
を
使えば,
となるので,
この式は
となり,結局上の二つの式()は
と書き換えられる。したがって
となることから
と求められる。これから
といった特徴が明らかになる。
次に,上半分の領域の材料が下半分のそれとは異なる場合の無限領域における,
その二つの領域境界を通した反射と透過について,
その概要を説明しておく。図-11.11に
示したように角度とベクトルを定義すると,
上半分の領域の透過波についても, 4で区別すれば,その変位を
同じ式(11.32)で表すことができる。
また透過波についての諸量は
である。ここに,,
が
上半分の領域のLaméの定数であり,が密度である。
これによる応力成分も前節同様に表現でき
となる。
このにおける条件は,二つの変位成分と,応力の
二つの成分, が,
それぞれ連続していることになるので
式(11.35) (11.43)の重ね合わせから
となる。ここでも前節同様,
と置いてある。
この条件式に式(11.33) (11.42)を代入して整理すれば
角度や振幅の間の関係式を得る。
まず指数関数同士が等価でなければならないから
となる。ここに
であり,は式(11.39)で定義した。
また角度は
と求められる。これを上式(11.44)に代入し直せば, 各振幅比を求めることができる。 以下面倒なので結果は省略する。
図-11.10に示した左からの入射波がSV波である場合を
考えよう。この場合には
となる以外は,式(11.33)と同じになっている。
これを用いて自由表面の境界条件に代入すると
となる。したがって,指数関数の形から
でなければならないため
(11.51)
と求められる。この結果を式(11.49)に代入して整理すると
となる。これを解けば
と求められる。この結果からは
といった特徴が明らかになる。
ところが,式(11.39)で
示したようにであることから,式(11.51d)からは
の条件を満足するように入射しない限り,が
実数にはならないことがわかる。では,
のときには
どのような現象になっているのか考察してみよう。この場合には
と置こう。これより
となるので,反射P波の変位成分は
となる。この波は,の正方向には 速さ で進む波になってはいるが, 深さ方向のの領域では, のように 指数関数的に振幅が小さくなる波になっている。 すなわち,浅い部分においてのみ顕著な振幅を持つ波になっている。 表面のみで顕著な波なので,表面波と呼ばれている。これについては 次節で検討しよう。
特に,透過する波の振幅の 増幅率 と角度を用いて, 地震波の増幅の仕方と角度の変化について議論したい。 文献[85]には右表に示したような地盤の基盤の加速度記録と 表層の加速度記録が示されており,その加速度振幅の増幅率は2.8倍になっている。 これを,最初に表面に到達するSH波の増幅率を用いて計算してみよ。 なお,簡単のために最初以外の反射波は無視していい。 入射角度はと60度の場合について求めてみよ。
層 | (kg/m) | (MN/m) |
表層 | 1500 | 12 |
第2層 | 1670 | 28 |
第1層 | 1850 | 56 |
基盤 | 1950 | 490 |