Complex Creek -Wikipedia、無料百科事典

before-content-x4

複雑な高調波運動 これは、単純な高調波運動の線形オーバーラップの動きです。単純な高調波運動は常に新聞ですが、複雑な高調波運動は必ずしも新聞ではありませんが、フーリエの調和のとれた分析によって分析できます。複雑な高調波運動は、周波数がすべて基本周波数の合理的な倍数である単純な高調波運動の組み合わせである場合にのみ新聞です。

複雑な高調波運動の運動学 [ 編集します ]

調和振動を提示するシステム n 一般的な自由度には伸びがあります バツ oフォームの独立した方向全体の動き:

1a ))

バツ t )) = j=1nc jAjcos おお jt + ϕ j)) {displaystyle mathbf {x}(t)= sum _ {j = 1}^{n} c_} mathbf {a} _ _ {j} cos(omega _ {j} t+phi _ {j})}}

または詳細:

1b ))

(x1(t)xn(t))= (A11A1nAn1Ann)(C1cos(ω1t+ϕ1)Cncos(ωnt+ϕn))= j(A1jAnj)c jcos おお jt + ϕ j)) {displaystyle {begin {pmatrix} x_ {1}(t)\ vdots \ x_ {n}(t)end {pmatrix}} = {begin {pmatrix} a_ {11}&cdots&a_ {1n} \ vdots&vdots&vdots&vdots&a_ {n1} {pmatrix}} {begin {pmatrix} c_ {1} cos(omega _ {1} t+phi _ {1})\ vdots \ c_ {n} cos(omega _ {n} t+phi _ {n} _ {1j} \ vdots \ a_ {nj} end {pmatrix}} c_ {j} cos(omega _ {j} t+phi _ {j})}

どこ、

{ ωi} = 初め n {displaystyle leftlbrace omega _ {i} rightrbrace _ {i = 1、…、n}}

です 独自の周波数 システムの、

{ ϕi} = 初め n {displaystyle leftlbrace phi _ {i} rightrbrace _ {i = 1、…、n}}}

初期フェーズ。各マトリックス列ベクトル a 呼ばれています 独自のやり方 振動の、そして c それらは、各片方の相対的な振幅です。そのために見ることができます n = 1複雑な高調波運動は、単に単純な高調波運動の合計です。

一般的な複雑な高調波運動の速度と加速は、時間に関して導出することによって得られ、また、同じ周波数の動きの複雑な高調波の動きであることが判明します。ただし、単純な高調波の動きのように、ゼロ速度ポイントは現在ありません。

after-content-x4

周期性 [ 編集します ]

運動は、定期的に時間の間隔で繰り返されると新聞が言われています。つまり、一定の時間間隔の後、同じ位置を通り、同じ速度で戻る場合です。周期性には、位置ベクトルが必要です バツ t )= バツ t+t ) すべてのために t の価値のために t 。次のような複雑な高調波運動の場合 1a )それはすべてのために必要です

cos おお it + t )) )) = cos おお it )) { k1,...kn} t = 2πkiωi= = 2πknωn{displaystyle cos(omega _ {i}(t+t))= cos(omega _ {i} t)rightarrow leftlbraceはexists k_ {1}、… k_ {n} rightrbrace Subset Mathbb {z}:t = {frac {{i} {i}} {} {{} {i} {{i}} {{i}} {frac {2pi k_ {n}} {omega _ {n}}}}}

周期性は、どんな周波数でも、その商が合理的な数である場合にのみ可能です。合理的な数値がゼロまたはヌルセットセットであるため、すべての周波数の商が合理的な数値である可能性はゼロであるため、実際の複雑な高調波の動きはQuasiperiódicosですが、新聞ではありません。

動き方程式 [ 編集します ]

によって与えられた複雑な高調波運動( 1a )o( 1b )これは、タイプの小さな振動の問題の方程式に対する解決策です。

m x¨t )) + k バツ t )) = 0 {displaystyle mathbf {m} {ddot {mathbf {x}}}(t)+mathbf {k} mathbf {x}(t)= mathbf {0}}

どこ:

線形減衰と内部励起力を持つシステムの最も一般的なケースでは、動きの方程式はより一般的です。

m x¨t )) + c x˙t )) + k バツ t )) = f t )) {displaystyle mathbf {m} {ddot {mathbf {x}}}(t)+mathbf {c}} {dot {mathbf {x}}}(t)+mathbf {k} mathbf {x}(t)= mathbf {f}(t)}

マトリックスが追加された場所

c {displaystyle mathbf {c}}

それは減衰を占めています。

結合振動 [ 編集します ]

複雑な高調波運動の一般的なケースは、 結合振動 。結合された振動のこの問題は、たとえば、結晶の熱振動、地震での建物の水平方向の動き、ドックまたはスプリングによる大規模な質量システムの動きに現れます。これらの問題は、次のタイプの方程式のシステムにつながります。

2 ))

それはマトリックスの方法で、次のように書くことができます。

2 ‘ ))

(m1000m2000mN)(x¨1x¨2x¨N)+ (k11k12k1Nk21k22k2NkN1kN2kNN)(x1x2xN)= 0 {displaystle {begin {pmatrix} m_ {1}&0&0&cdots&0 \ 0&m_ {2}&cdots&0 \ vdots&vdots&dots&vdots&vdots \ 0&0&0&cdots&m_ {n} end {pmatrix}} {{Pmatrix} {x} {x} {x} {x} {x}} _ {2} \ vdots \ {dot {x}} _ {n} end {pmatrix}}+{begin {pmatrix} k_ {11}&k_ {12}&k_ {1n} \ k_ {21} {22} {22} {2n} {2n} {2n} {2n} {22} {2n} {2n} {2n} { &vdots \ k_ {n1}&k_ {n2}&cdots&k_ {nn} end {pmatrix}}} {begin {pmatrix} x_ {1} \ x_ {2} \ vdots} end {pmatrix}}} = 0}}

問題は、通常の座標または振幅につながる変数の特定の変更によって解決できます。 独自のモード 実際、元の機械的問題の一般化された座標の特定の形態である振動の。

独自の周波数とモード [ 編集します ]

独自のモードは問題解決策を提供します( 2 ‘ ) 形 ( 初め )。このためには、次のように計算できるシステムの一連の自然周波数を決定する必要があります。

|k11ωj2m1k12k1Nk21k22ωj2m2k2NkN1kN2kNNωj2mN|= 0 {displaystyle {begin {vmatrix} k_ {11} -omega _ {j}^{2} m_ {1}&k_ {12}&cdots&k_ {1n} \ vdots&vdots&ddots&vdots \ k_ {n1}&k_ {n2}&cdots&k_ {nn} -omega _ {j}^{2} m_ {n} end {vmatrix}} = 0} = 0}

これは提供します n 固有周波数の正方形のソリューション。これらのソリューションのそれぞれについて、ユニットベクトルが求められ、呼ばれます 独自のやり方 、それは不確定な互換性のある方程式を満たします:

(k11ωj2m1k12k1Nk21k22ωj2m2k2NkN1kN2kNNωj2mN)(a1ja2jaNj)= 0 {displaystyle {begin {pmatrix} k_ {11} -omega _ {j}^{2} m_ {1}&k_ {12}&cdots&k_ {1n} vdots&vdots&ddots&vdots \ k_ {n1}&k_ {n2}&cdots&k_ {nn} -omega _ {j}^{2} m_ {n} end {pmatrix}} {begin {pmatrix} a_ {a_ a_ {a_ \ a_ \ a_ {a_ \ a_ \ a_ {pmatrix} {pmatrix} } end {pmatrix}} = 0}

システムのさまざまなモードを表すこれらのベクトルが互いに直交しているため、それらすべてによって形成されるマトリックスは直交行列であることを確認できます。

a = {aij}i=1…nj=1…n= (a11a1nan1ann){displaystyle mathbf {a} = {leftlbrace a_ {ij} rightrbrace} _ {i = 1 … n}^{j = 1 … n} = {begin {pmatrix} a_ {11}&cdots&a_ {1n} \ vdots&ddots&ddots&vdots} End {pmatrix}}}

通常の座標 、独自のモードに関連付けられており、従来の座標からの線形変化によって得られます。

(q1(t)qN(t))= (b11b1nbn1bnn)(x1(t)xN(t)){displaystyle {begin {pmatrix} q_ {1}(t)\ vdots \ q_ {n}(t)end {pmatrix}} = {begin {pmatrix} {11}&cdots&b_ {1n} \ vdots&vdots&vdots&vdots&vdots \ b_ {n1} {pmatrix}} {begin {pmatrix} x_ {1}(t)\ vdots \ x_ {n}(t)end {pmatrix}}}}

どこ

b = { b ij} = a t {displaystyle mathbf {b} = leftlbrace b_ {ij} rightrbrace = mathbf {a} ^{t}}

、それを満たします b ‘ の逆マトリックスです a a・b = b・a = )。

自由振動の問題の解決 [ 編集します ]

一般的なソリューションは、通常の座標で問題を解決することで簡単に取得できます。これらの座標を使用して、それを考慮すると簡単に取得できます

バツ = a Q {displaystyle mathbf {x} = mathbf {a} mathbf {q}}

m x¨+ k バツ = 0 q¨+ [ ATM1KA] Q = 0 {displaystyle mathbf {m} {ddt {xbf {x}}}+mathbf {k} mathbf {x} = 0qqquad rightarrow qud {ddot {q}}}+left [mathbf {a} ^{mathbf {a} mathbf {mathbf {-1bf {{-1bf {-1bf {-1bf {mathbf {{-1bf {{-1bf {-1bf} {a} right] mathbf {q} = 0}

しかし、マトリックスのプロパティによる

a {displaystyle mathbf {a}}

括弧内のマトリックスは斜めのマトリックスであることが判明しているため、その最後のシステムの解は解決することによって得られます n 型方程式の各セットの方程式:

q¨j+ おお j2Q j= 0 Q jt )) = c jcos おお jt + ϕ j)) {displaystyle {ddot {q}} _ {j}+omega _ {j}^{2} q_ {j} = 0qquad right arrow qquad q_ {j}(t)= c_ {j} cos _ {j} t+phi_ {j}}

通常の座標と独自のモードのマトリックスに関しては、システムの一般的なソリューションが書かれています。

(x1(t)xN(t))= (a11a1nan1ann)(q1(t)qN(t))= (a11a1nan1ann)(C1cos(ω1t+ϕ1)Cncos(ωnt+ϕn)){displaystyle {begin {pmatrix} x_ {1}(t)\ vdots \ x_ {n}(t)end {pmatrix}} = {begin {pmatrix} a_ {11}&cdots&a_ {1n} \ vdots&vdots&vdots&vdots&a_ {n1} {pmatrix}} {begin {pmatrix} q_ {1}(t)\ vdots \ q_ {n}(t)end {pmatrix}} = {begin {pmatrix} a_ {11}&cdots&a_ {1n} \ vdots&ddots&a_ a_ {n1} } end {pmatrix}}} {begin {pmatrix} c_ {1} cos(omega _ {1} t+phi _ {1})\ vdots \ c_ {n}

強制振動の問題の解決策が緩衝されました [ 編集します ]

一般的なソリューションは、通常の座標を使用して前と同様に取得されます

バツ = a Q {displaystyle mathbf {x} = mathbf {a} mathbf {q}}

m x¨+ c x˙+ k バツ = f q¨+ [ ATM1CA] q˙+ [ ATM1KA] Q = [ ATM1F] {Displaystyle Mathbf {m} {ddot {ddot {xbf {xbf {x} }+Mathbf {c} {dot {dot {dot {dot {xbf {x {x} }+mathbf {k} mathbf {f} Quad Rightarrow quad {ddf {Q} }+Left[ Mathbf {a} ^{t}mathbf {m} ^{1}mathbf {c} mathbf {a} right]{dot {a} right]{mathbf {Q} }+Left[mathbf {a} ^{t}mathbf {m} {-1}mathbf {k} mathbf {a} mathbf {a} mathbf {a} mathbf {a} mathbf =Left[mathbf {a} ^{t}mathbf {m} ^{-1}mathbf {f} right]}

乗算する直交座標マトリックスの構築により

q¨{displaystyle scriptStyle {ddot {mathbf {q}}}}}

q˙{displaystyle scriptStyle {dot {mathbf {q}}}}}

それらは対角線です(これには追加の条件として必要です

k c = c k {displaystyle scriptStyle mathbf {kc} = mathbf {ck}}

)、したがって、最後のシステムはに縮小されます n タイプの独立した方程式:

q¨j+ 2 n おお q˙j+ おお j2Q j= 0 Q jt )) = c j0tQjωj1νjそうです νjωj(tτ)おお j1νjt t )) t + ϕ j)) d t {displaystyle {ddot {q}} _ {j}+2nu omega {dot {q}} _ {j}+omega _ {j}^{2} q_ {j} = 0qquad rightarrow qqquad q_ {j}(j}(t)= c_ {j} {j} {j} {-q_ {j}} {omega _ {j} {sqrt {1-nu _ {j}}}}} e^{ – nu _ {j}} omega _ {j}(t-tau)} sin(omega _ {j} {sqr} {x} {x} {x} {x} {j} {j} {j} {j}) t+phi _ {j})dtau}

したがって、解決策は次のとおりです。

(x1(t)xN(t))= (a11a1nan1ann)(q1(t)qN(t)){displaystyle {begin {pmatrix} x_ {1}(t)\ vdots \ x_ {n}(t)end {pmatrix}} = {begin {pmatrix} a_ {11}&cdots&a_ {1n} \ vdots&vdots&vdots&vdots&a_ {n1} {pmatrix}} {begin {pmatrix} q_ {1}(t)\ vdots \ q_ {n}(t)end {pmatrix}}}}

つまり、の線形の組み合わせです n 強制高調波運動は緩和されました。

after-content-x4