強制振動の解
次に外力が加わったときの運動方程式の解を導出してみましょう。地震動はデジタルデータと呼ばれ、関数としては意味をもたないランダムなものです。
よって、地震動が外力の場合はニューマークβ法と呼ばれる近似解を得る数値解析で応答解を求めます。
このニューマークβ法は後に勉強するとして、ここでは、調和外力という周期性を持つ外力が作用したときの応答解を求めてみましょう。
通常の建物には減衰という地震応答を低減させる復元力が発生します。調和外力は周期性を持つ外力で、これはsin波やcos波にあたります。
振動している状態で、力のつり合いを考えると、図のようになります。まず、質量×重力加速度が物体の力として作用します。この方向をプラスと定義しています。また、この力が作用したとき元の状態に戻ろうとするバネの力、復元力が逆方向に発生します。この図は振動している一瞬を取り出したものと考えてください。つまり、この状態で力=復元力+減衰力という等式が成り立つはずです。よって、以下の式で表すことができます。
mx''+cx'+kx=f(t)
a=x''
f(t)=Fsinpt
mx''+cx'+kx=Fsinpt
x''+(c/m)x'+(k/m)x=(F/m)sinpt
k/m=w2
c/m=2hwとすれば
x''+2hwx'+w2x=(F/m)sinpt
と変形できます。ここで、wは固有円振動数と言われています。さて、(k/m)=w2、(c/m)=2hwという等式も、ここでは定義としているので深く考えないように進みましょう。
さて、以上のように示した微分方程式の一般解は「斉次方程式の一般解+特殊解」で示されます。
スポンサーリンク
 
※斉次方程式とは、右辺の値が0である微分方程式です。例えば、自由振動を表す振動方程式はそれに当たります。
さて、周期的な外力sinptが作用した場合の応答解を以下のように仮定します。
x=
Asin(
pt-
θ)
また、運動方程式のような動的な問題は時間に依存するために、tの時間関数であることに注意してください。Aは適当に仮定した定数です。
先ほどの式は加法定理を用いれば次のように表すこともできますね。
※加法定理:sin(α±β)=sinαcosβ±cosαsinβ
Asin(pt-θ)= Asinptcosθ-Acosptsinθ
ここで、定数はそれぞれ次のように示します。
Acosθ=a
Asinθ=b
よって、
Asin(pt-θ)=asinpt-bcospt
です。さて、何がしたいかわかってきましたか?要は応答解を求めたいわけですよね。応答解は以下の式で仮定しました。また、以上の関係から、
tan
θ=
b/
a
です。
以下の式で未知数なのはa,bです。よって、a,bを求めれば特殊解が求められるわけですね。まず、振動方程式に代入するために、1,2階の微分を行います。
x=
Asin(
pt-
θ)=
asin
pt-
bcos
pt
x'={
Asin(
pt-
θ)}' =
apcos
pt+
bpsin
pt
x''={
Asin(
pt-
θ)}'' = -
ap2sinpt+
bp2cos
pt
ですから、以上の式を代入すると、
x''+(
c/
m)
x'+(
k/m)x=(F/m)sin
pt
-
ap2sinpt+
bp2cos
pt+2
hw(
apcos
pt+
bpsin
pt)+w
2(
asin
pt-
bcos
pt) = (F/m)sinpt
ですね。さて、ここでsinptとcosptに分けてみましょう。上記の式は、以下の式を足し合わせたものなので、分けても問題ないですよね。
bp2cos
pt+2
hwapcos
pt-
w2bcos
pt = 0
-
ap2sin
pt+2
hwbpsin
pt+
w2asin
pt= (
F/
m)sin
pt
cos
pt及びsin
ptを消します。
bp2+2
hwap-
w2b = 0
-
ap2+2
hwbp+
w2a= (
F/
m)
式を変形すれば、
b=-2
hwap/(
p2-
w2)
です。
bに代入します。
-
ap2+2
hwp(-2
hwap/(
p2-
w2))+
w2a= (
F/
m)
次にaの結果をbに代入すると、
b=-2
hwap/(
p2-
w2)
となります。
なので、aとbを代入します。
となり、この式が強制外力に対する特殊解です。斉次方程式の一般解は減衰自由振動の解ですから、それを足し合わせたものが、強制外力に対する一般解となります。
よって、強制外力の一般解は

です。
また、もうひとつの未知数θを求めます。
tanθ=b/a
以上が強制外力に対する一般解の求め方です。いかがでしたか?自由振動の
導出よりは少しだけ難解のように感じるかも知れませんが、要は慣れです。
何回も読んでいるうちに自然と理解できるようになると思います。
また、用語や導出過程は、かなりくどく説明しているつもりです。
急に、違う数式に飛ぶのではなく、順序を追って数式を変化させているので、
じっくりと読めば必ず理解できるので、頑張ってみましょう!
【管理人おすすめ!】セットで3割もお得!大好評の用語集と図解集のセット⇒ 建築構造がわかる基礎用語集&図解集セット(※既に26人にお申込みいただきました!)
▼こちらも人気の記事です▼
▼人気の記事ベスト3▼
▼用語の意味知らなくて大丈夫?▼
▼同じカテゴリの記事一覧▼
▼カテゴリ一覧▼
▼他の勉強がしたい方はこちら▼
公式LINEで気軽に学ぶ構造力学!
一級建築士の構造・構造力学の学習に役立つ情報を発信中。
【フォロー求む!】Pinterestで図解をまとめました
図解で構造を勉強しませんか?⇒ 当サイトのPinterestアカウントはこちら