• ベストアンサー

振動方程式の解き方

現在0から独学で内部摩擦について勉強中なのですが、KGLモデルによると内部摩擦による振動方程式は下式のようになるそうですが、なぜこのような式になるのか分かりません。簡単で良いのでこの式の説明と解の求め方をご教授ください。 参考書やホームページの紹介やアドバイスなど何でも良いのでよろしく御願いします。 ちなみに波動方程式の解き方などはだいぶ前に習ったためすっかり忘れており、現在勉強中です。 A*∂^2u/∂t^2 + B*∂u/∂t - C*∂^2u/∂x^2 = σb u:displacemet A:effective mass per unit length B:viscous resistance C:line tension b:バーガーズベクトル 解は下式のようになるそうです。 u = 4bσ_0/πA*sin(πx/L)*exp[i(ωt-φ_0)]/√{(ω_0^2-ω^2)^2+(ωd)^2} (ω_0 = (π/L)√(C/A) , tanφ_0 = ωd/(ω_0^2-ω^2) , d=B/A)

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

  • ベストアンサー
回答No.1

速度に比例する抵抗と外力が働く場合の波動方程式という程度の 質問ではないんでしょうね。 検索しても「KGLモデル」では一件もかかってこないので、 その方面の専門書を見るしかないんじゃないでしょうか？ 解法は偏微分方程式関係の本を見ればよいのではないでしょうか。 多分、σは応力なんだろうと思いますが、解の形から見ると、 e^[-i ω_0 t]の強制振動っぽいですね。 u=(定数)e^{i(kx-ωt)}とおいて代入し、初期条件と境界条件を 適切に設定すれば、解けるような気はしますが。

その他の回答 (1)

回答No.2

KGLモデルのことはよく知りません．が， google scholar で"KGL Model"と検索するとたくさん論文が引っかかります．ぱっと見たところ日本語で解説されているものはないようなので，一生懸命よんでみるしかないと思いました． こういう微分方程式は大抵フーリエ変換というのを使って解くと楽なことが多いですよ．実際には計算してませんが，解の形からもきっとそうです．っというわけで，物理数学か解析学かなんかの教科書でフーリエ変換を使って微分方程式を解いているところをみればOKだとおもいますよ．

参考URL：
http://scholar.google.co.jp/schhp?hl=ja
  1. カテゴリ
  2. 学問・教育
  3. 自然科学
  4. 物理学

関連するQ&A

  • ３次方程式の解き方

    x^3+11×10^-2x^2-10^-14ｘ-10^-16=0 この方程式を解きたいのですが、解き方を調べてみたところ、p=1/a(c-b^2/3a) q=1/a(d-bc/3a+2b^3/27a^2) これからｐ・ｑをだして t^2+qt-p^3/27=0 このｔの３乗根をAとして、B=-p/3AとしてA・Bを求める。そして、A+B-b/3a Aω+Bω^2-b/3a Aω^2+Bω-b/3a ただし、ω=-1/2+√3i/2 この３つの式を解けばそれが解となる。と書いてあったんですが、この解き方で解けるんでしょうか？

  • 減衰振動の微分方程式の解

    先生から配られたプリントには減衰振動の微分方程式が「m(dx/dt)^2+2γ・dx/dt+ω^2x=0」の時、解が「x=A・EXP（-γt）cos（ω´t＋φ）」って書かれてます。 摩擦：гならг/m=2γ、バネ定数：kならk/m=ω^2、A=√C1^2+C2^2、φ=C1/C2、ω´^2=ω^2-γ^2です。 解の式で、cosじゃなくてsinではないのですか？単振動・強制振動の場合も同様にcosでした。 誰かよろしくお願いいたします。

  • 3次方程式について

    a,b,c,に関して a＋b＋c=(a＾2)＋(b＾2)＋(c＾2)=(a＾3)＋(b＾3)＋(c＾3)=n が成り立つとき。(a＾4)＋(b＾4)＋(c＾4)=nが成り立つとき、nの値を求める。 問題で、ab + bc + ca = n(n-1)/2　　 から３次方程式の解と係数の関係を用いて。 a + b + c = n。更にab + bc + ca = u, abc = vと置いておきます。 解と係数の関係より、a, b, cはtの３次方程式t3-nt2+ut-v=0 t3-nt2+ut-v=0 の式になる方法がわかりません。 お願いします

  • 平面の方程式

    3つの点、S(1)(3)(1)、T(1)(2)(1)、U(2)(-1)(1)を通る平面の方程式を求めよ。 という問題なのですが、下の解答で合っているでしょうか？ どこか間違えていたりしていたら、教えて頂きたいのですが…。 解）求める平面の方程式をax+by+cz+d=0とおく。…(1) (a≠0またはb≠0またはc≠0とする。) 点Sを通るから、a+3b+c+d=0…(2) 点Tを通るから、a+2b+c+d=0…(3) 点Uを通るから、2a-b+c+d=0…(4) (2),(3),(4)の連立方程式を解く。 同次ゆえ、係数行列に行の基本変形をしても、解の集合は変わらない。 (1 3 1 1)(1 2 1 1)(2 -1 1 1) →(0 1 0 0)(1 2 1 1)(1 -3 0 0) →(0 1 0 0)(0 5 1 1)(1 -3 0 0) →(0 1 0 0)(0 0 1 1)(1 0 0 0) 方程式に直すと b=0,c+d=0,a=0 よって b=0,c=-d,a=0 ここで、d=k（k:自由定数）とおくと c=-k K=0とすると、a=0かつb=0かつc=0となり、不適当。 よって、k≠0である。 (1)に代入して -kz+k=0 両辺をkで割って -z+1=0 z=1 これが求める平面の方程式である。

  • 単振動の運動方程式の解（複素数表示？）

    中学の者ですが、なんとか独学でここまで理解しています。 答えの載っていない参考書を持っていて、 それを読みながら勉強しているのですが、 ある問で、 微分方程式m(d^2x/dt^2)＝-kt の解は x=A*exp(iωt)+B*exp(-iωt)　　　（ただしω＝√(k/m) ） の形で表わされることを示せ というのがありました。 微分方程式の解き方は分かっていたので、素直に x＝C*sin(ωt+D)　　　（CとDは定数） としました。 ここからどうやって示すべき式に持っていくのでしょうか。 見当がつきません。 それから、速度をあらわす式vを時間tで表わし、 t=0のときx=a、t=0のときv=0という条件で、A,Bを aとωを用いて表せというのもありました。 これは与えられた x=A*exp(iωt)+B*exp(-iωt)　 をそのままtで微分していいということなのでしょうか。 機械的にやってみたのですが、AもBもa/2という結果になり、 ωは出てきませんでした。 いま一つ問題が何を読者に気付いてほしいのかわからないことと、 最初に書きました表示の部分が分かりません。 どなたか詳しく教えていただけませんか。お願いします。

  • 二次方程式

    方程式x^2-2x-5=0の解をa,b（a<b)、方程式x^2-2x-7=0の解をc,d(c<d)とするとき (a-c)(a-d)(c-a)(c-b)の値を求めなさいとあります。 ここで答えを見ると {a^2-(c+d)a+cd}{c^2-(a+b)+ab}・・・・(1) ここで解と係数の関係より a+b=2, ab=-5;c+d=2,cd=-7 であるから、(1)=(a^2-2a-7)(c^2-2c-5)・・・・・・・・(2) さらにa^2-2a-5=0, c^2-2c-7=0 が成り立つことから a^2-2a=5, c^2-2c=7 ∴(2)=(5-7)X(7-5)=-4(答） と書いてありますがいくつかわからない部分があります。 まず解と係数の関係のはなしですが、a+b=2　とありますが　a+b=-2/xで-2だと思うのですがなぜ2なのですか？ それと a^2-2a-5=0, c^2-2c-7=0 の部分がわかりません。　 助けてください

  • 線形代数　連立方程式

    連立方程式の問題なのですが、方向性はわかっていると思うのですが、うまい形にもっていけません。教えてください。 次の連立方程式が解を持つようにtを定め、そのときの解を答えよ。 a-2b+2c-d+2e=2 4a-9b+8c-2d+7e=8 3a-7b+4c-d+7e=t a-2b+3c-d-e=2

  • 五次の不定方程式の問題

    A,B,C,D,Eを整数（負でもよい）とします．この時，五次の不定方程式 A^5 + B^5 + C^5 + D^5 = E^5 ・・・・・・・式(1) の整数解として， 27^5 + 84^5 + 110^5 + 133^5 = 144^5 ・・・・・・・式(2) が知られています．では，式(1)の解は一般的に解かれているのでしょうか？ 私は，式(2)の数値解しか知らないのですが，どなたか教えて下さい． 余談ですが，因みに，以下のような不定方程式もあります． 1020^4 = 400000^2 + 980^4 + 20^6, 18^4 = 14^4 + 16^4 + 4^5, 51^4 = 10^2 + 49^4 + 10^6, 4020^4 = 16000^2 + 3980^4 + 400^5, 9^4 = 2^6 + 7^4 + 8^4, 15^4 = 200^2 + 5^4 + 10^4, 35^4 = 800^2 + 15^4 + 30^4, 12^5 = 4^3 + 60^3 + 8^5, 85^4 = 3200^2 + 5^4 + 80^4 + 10^6, 4^3 = 3^3 + 6^2 + 1, 5^4 = 19^2 + 2^3 + 4^4, 7^6 = 10^3 + 4^5 + 10^5 + 5^6.

  • 強制振動

    m×d^2x/dt^2＋mω。^2x＋2mν×dx/dt=Fcosωt...(1) 1.強制振動の場合の一般解が、F=0とした斉次方程式の解と、Fの入った方程式の特解の和で与えられることを示せ。 2.特解を求めよ。 3.特解の振幅について、外力の振動数を変えたときどうなるかしらべよ。 4.特解の位相と外力の位相の関係を、外力の振動数が小さいときから大きくしていく場合についてどうなるかを議論せよ。 自分の解等 1.わかりません 2.(1)×1/mより 　d^2x/dt^2＋ω。^2x＋2ν×dx/dt=F/m×cosωt 特解をx=Acosωtとおき上の式に代入する 　-ω^2Acosωt＋ω。^2Acosωt-2νωAcosωt=F/m×cosωt F/m=fとすると 　A=f/-ω^2＋ω。^2-2νω よって特解は 　　　x=f/-ω^2＋ω。^2-2νω×cosωt 3. 振幅はω<ω。で正、ω>ω。で負の値をとるがωに近づくにつれ、そ 　の絶対値は無限大に発散する。 4.わかりません。

  • 物体に摩擦力が働く調和振動

    <<問題>> 水平な床の上を、ばね定数kのばねで繋がれた質量mの物体が運動する場合を考える。ばねの自然の長さからの変位をxとし、x軸を右方向が正となるように選ぶことにする。時刻t＝0においてx=X（X≧０）で静かに手を離す場合の物体の運動について、次の問いに答えなさい。動摩擦係数をμ 重力加速度をgとする。 摩擦力が働いている場合、運動方程式は非同次方程式となる。その特別解をx=Aの形に仮定し,xが解となるようにAをμ,g,m,kで表しなさい。 <<解法>> 摩擦力が０の場合の一般解を求めて、それを摩擦力≠０の時の運動方程式に代入していくと 　2 mω (Acosωt＋Bsinωt)＝-k(Acosωt＋Bsinωt)＋F ここまでは分かるのですが、ここからどのようにして Aを表していけばいいかがわかりません。 お願いします。

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する