• ベストアンサー

固有振動数を求めたい

あるシステムの運動方程式を導くと、 mx"1+k(2x1-x2)=f1 mx"2-k(x1-2x2)=f2 となりました。("はツードットを意味します。xの右の数字は添え字です。) このシステムの固有振動数を求めたいのですが方法がわかりません。 参考サイトでもよいので教えていただけませんか? また、この固有振動数の時のシステムの振る舞いについて述べよという問題もあるのですが、こういう場合、文章で書くのですか?

質問者が選んだベストアンサー

  • ベストアンサー
  • reiman
  • ベストアンサー率62% (102/163)
回答No.2

#1でもいいのですが 一般的には複雑な線形結合にしなければならないので 本筋の解き方は x= ┌ ┐ │x1│ │x2│ └ ┘ A= ┌-- --┐ │ 2 -1│ │-1 2│ └-- --┘ f= ┌ ┐ │f1│ │f2│ └ ┘ とおくと mx"+kAx=f あとは 正則行列Pを使って座標変換 y=Px を行い g=P^-1f とおくと my"+kP^-1APy=g P^-1APが対角行列ΛになるようにPを決定すれば良い。 すると my"+kΛy=g すなわち my1"+kλ1y1=g1 my2"+kλ2y2=g2 となり簡単な非ベクトルの定係数2階線形非斉次線形方程式になる。 ただし y1:yの第1成分 y2:yの第2成分 g1:gの第1成分 g2:gの第2成分 λ1:Λの1行1列 λ2:Λの2行2列 要するにこの問題は行列の対角化の問題です。 固有振動数は斉次の解を求めればすぐに求まります。

charmin
質問者

お礼

回答ありがとうございます。 なんとかできました。

その他の回答 (1)

  • yokkun831
  • ベストアンサー率74% (674/908)
回答No.1

辺々引けば,x1-x2について, 辺々足せば,x1+x2について強制振動の方程式が 得られますよね。それぞれ,相対運動および 重心運動を示す方程式になっていると思います。 それぞれの固有振動数は明らかです。後半は, 運動方程式をもとに,それぞれの振動数のときの x1,x2のふるまいをのべればいいのでは? 外力があるので,共振として議論することに なるのかな?

charmin
質問者

お礼

回答していただきありがとうございます。 振動工学なるものを受講したことがなく詳しいことが分からないのですが、書いていただいたことを参考にして解いてみます。

関連するQ&A

  • 固有円振動数と固有振動数

    固有円振動数と固有振動数、違いが理解できません。 固有円振動数とは、簡単に言うと何ですか? そして、固有振動数とは、簡単に言うと何ですか? 共振は、「振動体に固有振動数と等しい振動が外部から加わると、振動の幅が大きくなる現象」とあります。 共振曲線を見ると、縦軸が応答倍率で、横軸が振動数ω/ω0とあります。 なので、だいたい ω/ω0=1 のところで応答倍率は最大となるので、振動体・外部からの振動の固有円振動数が一致したときを「共振」というのだと思っていました。しかも何故横軸が「振動数f/f0」ではなく「振動数ω/ω0」なのでしょうか。 よく理解できていないので、変な質問になってしまっているかもしれませんが、よろしくお願い致します。

  • 振動です

    図のように、2つのばねk1,k2と2つの質量m1,m2がとり付けられています。xo(t)=Xosinωtで支持部が変位するとき (1)この振動系の運動方程式を導出しなさい 私の回答 m1(d^2 x1 /dt^2)=-k1(x1-xo)+k2(x2-x1) と m2(d^2 x2 /dt^2)=-k2(x2-x1) (2)k1=2k , k2=k ,m1=m ,m2=m/2 の時、設問(1)の運動方程式はどう書き直せるか。ωo=√(k/m)を用いて記述せよ 私の回答 (d^2 x1 /dt^2)=-2(ωo^2)(x1-xo)+(ωo^2)(x2-x1) (d^2 x2 /dt^2)=-2(ωo^2)(x2-x1) (3)ωo=1rad/sのとき、設問(2)の運動方程式を用い、固有角振動数ωnを求めよ。 この問題を見た瞬間、あれωoが(不減衰)固有角振動数でないの?って思いました。ωoとωn何が違うのですか? それと設問(1)(2)は正しいですか?教えてください

  • 運動方程式 m(d^2x/dt^2)+kx=0 の固有角振動数、固有周

    運動方程式 m(d^2x/dt^2)+kx=0 の固有角振動数、固有周期、固有振動数ってどうやって求めるんですか?

  • いろいろな物の固有振動数が知りたいです。

    いろいろな物体の固有振動数が知りたいです。 身の回りにあるいろいろな物の固有振動数が知りたいのですが、調べてみても「人間の固有振動数は8Hz」というものしか出てきませんでした。 皆さんが知っているいろいろな物の固有振動数を教えてくださりませんか?なんでもいいし、正確な数字じゃなくてもいいです。

  • バネの固有振動数を求める問題が分かりません。

    バネの固有振動数を求める問題が分かりません。 下の図の(a)~(d)に示す系の固有振動数を求める問題が分かりません。ばね x=Ae^jωtと置いて計算していき (a)は多分 ω=√(k1/m) となり f=1/2π・√(k1/m) になりました。 (b)も同様に ω=√{k1k2/m(k1+k2)} となり f=1/2π・√{k1k2/m(k1+k2)} (c)(d)が分かりません。 よろしくお願いします。

  • 大学の数学の勉強をしていてどうしても分からない問題があるので教えてください。

    運動方程式まででも良いのでお願いします。 同じ質量mを持つ3つの質点A、B、Cが直線上に左から右へと並んでおり AとB、BとCはバネ定数kのバネで繋がれている。 さらにAは左側の壁と、Cは右側の壁との間にもバネ定数kのバネで繋がれている。 (重力は考えなくてよい) (1) A、B、Cそれぞれの釣り合いの位置からの変位をx1、x2、x3とする。 A、B、Cそれぞれの運動方程式(微分方程式)を書き下ろし ω=k/mX=(x1 x2 x3) として3つの運動方程式を((d^2)X)/(dt^2) = -(ω^2)AX の形にまとめると、 A=( 2 -1 0 -1 2 -1 0 -1 2 ) となることを示せ。 (2) この系の固有振動の振動数と各固有振動における各質点の相対的変位 (振動モード) を求めよ。

  • 調和振動子の固有値

    壁|---○---○---|壁 上の図のように、ばね定数Kの3つのバネ3本で連結された二つの粒子(質量はともにm)の運動のうちで、バネに平行な方向の成分だけに着目して、そのハミルトニアンを     重心座標 X=(x1+x2)/2 相対座標 x=x1-x2 を用いて書き直すと H=[{-(hbar)^2/4m}(∂/∂X)^2 +KX^2]-[{(hbar)^2/m}(∂/∂X)^2 + 3Kx^2/4}] となる。 ここから固有値を求めるわけですが、解説には 重心運動は 質量 2m ばね定数 2K 相対運動は 質量 m/2 ばね定数 3K/2 となってωo=√(K/m) ω=√(3K/m) 固有関数は、それぞれの固有関数の積、固有値は和であらわされるから、 εNn=(N+1/2)(hbar)ωo + (n+1/2)(hbar)ω と書いてあります。 ここで二つ質問があって、 (1)重心運動のばね定数が 2K 相対運動の 質量 m/2 ばね定数 3K/2   という風にどうして表されるのか、 (2)なぜ調和振動子の固有値はそれぞれの和で表されるのか というのが分かりませんでした。 どなたかご教授ください。お願いします。

  • 入力の振動数と、固有振動数と、出力の振動数の関係

    入力の振動数をf1 振動体の固有振動数をf0 出力の振動数をf2とすると 振動体が自由振動をするときの この3つの振動数の関係を教えてもらえないでしょうか?お願いします

  • 固有振動数 二重振り子

    ラグランジュを使って二重振り子の固有振動数を求めたいです。 運動方程式をつくることは出来ました。 しかし、ここからどうやって固有振動数を求めていいのかがわかりません。 1つ目の振り子が長さl1重さm1角度シータ1 2つめの振り子が長さl2重さm2角度シータ2で 微小なゆれなので (m1+m2)l1シータ1(ドット2つ)+m2l2シータ2ドット=ー(m1+m2)gシータ1 l2シータ2(ドット2つ)+l1シータ1(ドット2つ)=ーgシータ2 ググルと色々出てくるのですが、 運動方程式→ω= みたいにすぐ出ているので、ωを出す過程がわかりません。 1つだけ次の振動数方程式を解く、と (m1+m2)(ーω^2l1+g)(ーω^2l2+g)ーω^4m2l1l2=0 という式が乗っていたのですが、 何でこれをとくのかがさっぱりです。 これは他のでは使えませんよね?? 万能な振動数方程式というのは存在するのですか? あと超初歩ですが シータ(ドット)=ω シータ(ドット2つ)=ωドット=角加速度 だと思っているのですが、振動数方程式には角加速度は出てきません。どうしたらいいのでしょうか? 色々質問がありますが、解く過程が知りたいです。 詳しくお願いします!!!!

  • 振動数方程式

    3自由度系の振動の問題で 振動数方程式を立てたのですが、解けなくて困っています。 (行列式=0) |Iλ^2+2k       -k       0 |        | -k      I'λ^2+2k        -k  | =0      | 0        -k         Iλ^2+2k  |            λ=...... (I=ma^2/2, I'=8ma^2/3π) 慣性モーメント これから固有振動数(低い順にw1,w2,w3として)を求めたいのですが、 方程式が解けないので分かりません。 できれば固有振動数も教えてください。お願いします 結果は(m:質量,k:ねじり剛性,a:半径)剛体の寸法でお願いします。