建築学生が学ぶ構造力学

建築学生が学ぶ構造力学
  1. HOME > 耐震設計の基礎 > 非減衰自由振動の解

非減衰自由振動の解とは|運動方程式の解法と固有振動数・周期の求め方を解説

この記事の要点

建物の地震応答を理解するために、まず最も単純な「外力なし・減衰なし」の振動から出発するのが構造力学の定石だ。

この解の形を理解しておくと、減衰付き・強制振動付きへの拡張が見通しやすくなる。

この記事では非減衰自由振動の運動方程式の解法・固有振動数・固有周期の導出を解説する。

固有円振動数ω=√(k/m)の定義とオイラーの公式を使った解法の流れを押さえておきましょう。

減衰自由振動の解と対比することで理解が深まります。

【管理人おすすめ!】セットで3割もお得!約1,100語の用語集+476点の図解集セット⇒ 建築構造がわかる基礎用語集&図解集セット


初めにもっとも単純な運動方程式である非減衰自由振動の解を導出してみましょう。


通常の建物には減衰という地震応答を低減させる復元力が発生しますが、これを無視して考えてみます。


自由振動であれば外力は0です。この状態での振動は、例えば、物理で習った振り子の問題で、


自分では手を加えずにある位置から玉を押さえておいて離したときに発生する振動と同様のものです。


振動している状態で、力のつり合いを考えると、図のようになりあます。まず、質量×重力加速度が物体の力として作用します。


この方向をプラスと定義しています。また、この力が作用したとき元の状態に戻ろうとするバネの力、復元力が逆方向に発生します。


この図は振動している一瞬を取り出したものと考えてください。つまり、この状態で力=復元力という等式が成り立つはずです。


よって、以下の式で表すことができます。


ma=-kx

ma+kx=0


ここで、mは質点の重量、aは加速度、kはバネ剛性xは変位を表しています。また、加速度は変位の二階微分と同様なので

非減衰自由振動のモデル

a=x''

mx''+kx=0

x''+(k/m)x=0

k/m=w2

x''+w2x=0


と変形できる。ここで、wは固有円振動数と言われています。なぜ、 (k/m)=w2としているかといえば計算しやすくするためですね。


さて、以上のように示した運動方程式は微分方程式であり非減衰自由振動解はこの方程式を解くことで得られます。数学的な手法からこのような問題を解く場合、


x=Aeλt


と考えます。この式はあくまでも物体が自由振動を起こしたときの応答を仮定したものです。


また、運動方程式のような動的な問題は時間に依存するために、tの時間関数であることに注意してください。


Aは適当に仮定した定数で、λは未知数です。


x=Aeλt

x'=Aλeλt

x''=2eλt

x''+w2x=0

2eλt+w2Aeλt=0

Aeλt(λ2+w2)=0


よって、以上の式が満たすようなλを求めればいいことになります。 Aeλtがゼロになるときはλの値が∞になる場合です。


このとき、変位xも無限となることから値は発散し解にはなりえません。よって、

(λ2+w2)=0

λ2=w2

λw

x=Aewit+Be-wit

となります。±の値が表れたのでこれを重ね合わせた上記に示す式となります。さらに、オイラーの公式を用いて三角関数を用いた形で表してみましょう。


※オイラーの公式とは、eがcosθとsinθの式で表すことができることをテイラ―展開で証明した等式です。

e±=cosθ±isinθ

x=Aeiwt+Be-wit

=A(coswt+isinwt)+B(coswt-isinwt)

=Acoswt+Bcoswt+Aisinwt-Bisinwt

=Acoswt+Bcoswt+Aisinwt-Bisinwt

=(A+B)coswt+i(A-B)sinwt

A1=(A+B),A2=(A-B)

=A1coswt+iA2sinwt


以上のように、非減衰自由振動の解を求めることが出来ました。理解できないところを集中的に勉強しましょうね。


また、掲示板にも質問を受け付けています。

気軽に書き込みしてください。

【管理人おすすめ!】セットで3割もお得!約1,100語の用語集+476点の図解集セット⇒ 建築構造がわかる基礎用語集&図解集セット

▶ 建築士試験ではどう問われる?
この論点は建築士(構造)で出題されます。実際の問題と解き方は 構造力学・荷重・耐震の過去問解説(分野別・全年度)(無料)で確認できます。

理解度チェック

Q.

非減衰自由振動の運動方程式は?

答えを見る

ma=-kxよりma+kx=0です。加速度a=x''なのでmx''+kx=0、変形してx''+(k/m)x=0となります(m:質量、k:バネ剛性)。

Q.

k/mは何と定義されますか?

答えを見る

k/m=ω2(ωは固有円振動数)と定義します。これによりx''+ω2x=0となります。

Q.

この微分方程式の解はどんな形になりますか?

答えを見る

x=Aeλtと仮定すると特性方程式λ2+ω2=0よりλ=±iωとなり、オイラーの公式(e±iθ=cosθ±isinθ)を用いて三角関数の形で表せます。

ハナダユキヒロ

この記事を書いた人

ハナダユキヒロミツメラボ

設計事務所に7年勤務。2010年より「建築学生が学ぶ構造力学」を運営(16年以上)。著書「わかる構造力学」「わかる構造力学(改訂版)」(工学社)。

▼用語の意味知らなくて大丈夫?▼

建築学生が学ぶ「構造力学」の用語集

▼同じカテゴリの記事一覧▼

▼カテゴリ一覧▼

▼他の勉強がしたい方はこちら▼

▼調べるたびに検索するのをやめたい人へ▼

▼過去問の答えは見たけど、解き方がわからない人へ▼

プロフィール

建築の本、紹介します。▼

すぐにわかる構造力学の本

同じカテゴリの記事一覧

Topへ >>

  1. HOME > 耐震設計の基礎 > 非減衰自由振動の解
  2. 1級の過去問(計算)解説
  3. わかる建築構造の用語集・図解集
  4. 1頁10円!PDF版の学習記事