スプリングマスダンパーシステムの伝達関数を見つける


9

私は緒方近代制御工学の本を読み、基本的な制御原理の理解を深めるためにいくつかの演習を行ってきました。私は解決するのに苦労している次の例に出くわしました。

この振動治具をモデル化した伝達関数を考え出す必要があります。質問は次のとおりです。

この例では、振動試験装置を分析します(図1)。このシステムは、質量Mのテーブルと、質量がmのコイルで構成されています。地面にしっかりと取り付けられた永久磁石は、安定した磁場を提供します。磁場を通るコイルTheの動きは、式(1)のように、その速度proportionalに比例する電圧をコイルに誘導します。1.𝑒= 𝛼𝑦̇ [eq.1]

コイルに電流が流れると、コイルには電流に比例する磁力が発生します。2.𝐹= 𝛽𝑖 [eq.2]

ここに画像の説明を入力してください

質問:出力withから入力𝑉へのパラメトリック伝達関数を取得します。

私が答えるのは難しいが、TF全体に影響するいくつかの質問は次のとおりです。

  • K2とB2が距離Zだけ圧縮されている場合(
    コイルが磁場と相互作用するために上向きに移動する場合)、これはk1とb1が同じ距離Zだけ伸びていることを意味しますか?

  • 場合m(コイル)2センチメートルによって上方に移動しないM(表)は、2センチメートルによって上方に移動しますか?


私は何をする必要がありますか:

  • テーブルの質量Mとコイルの質量mの2つの独立した自由ボディダイアグラムを作成します。
  • 逆起電力を含む1つの回路図をスケッチします。
  • sドメインに変換します。
  • 同時に解決します。

これまでに行ったこと:

  • 自由体図を分離して描画し、方程式を抽出します。

  • 回路図を描き、方程式を抽出します。

  • sドメインに変換します。

MATLAB関数solveを使用して、2つの異なる5次伝達関数(以下で提案する方法ごとに1つ)を取得することができましたが、どちらが正しいのか、またその理由はわかりません。


全体システム:

これは、電気部品を除いて、振動試験治具をどのようにモデル化できるかを図式的に表したものです。

ここに画像の説明を入力してください


自由体図1-表-上向きの規則

スプリングk1k2ダンパーb1b2され、別々にモデル化。それらを一緒に追加して1つとして表示することはできないため、圧縮と拡張は別々です。

上向きの力から来ているk2b2コイルに取り付けられています。これらは上向きの動きを経験しています。

ここに画像の説明を入力してください

sドメインの方程式:

Ms^2X + b1sX + k1X = b2s(X-Y) + k2(X-Y)


フリーボディダイアグラム2-コイル-上向きの規則

コイルには上向きの力が加わっていますが、スプリングとダンパーがそれを抑制し、反対方向に作用しています。

ここに画像の説明を入力してください

sドメインの方程式:

Fem = Ms^2Y + b2s(X-Y) + k2(X-Y)


上記の表のFBDの2つの異なる方法は、s領域の異なる方程式と異なる伝達関数につながります。

テーブルとコイルの正しい自由体図は何ですか?


2
いい質問ですが、拡大するためにクリックせずに詳細が明確な写真を投稿してください。たとえば、これらのマイナス記号はほとんど識別できません。さらに、左下の方程式は部分的に切り取られています。シートを大きくするために使用するシートには、多くの空き領域があります。インターネットには多くの無料の画像編集プログラム(IrfanViewやFastSstone ImageViewerなど)があるので、シートの写真を何枚か撮り、必要な部分をカット/クロップして素敵な写真を投稿することもできます。
Lorenzo Donati-Codidact.org 2018年

@LorenzoDonati、提案に感謝し、すぐに編集します。左下の方程式については、自由体の図が問題なので、これは重要ではありません。それが正しければ、方程式は正しくなります。ただし、それに応じて編集を試みます。ご意見をいただきありがとうございます。
rrz0

自分が間違ったことを推測しないようにしてください。一連の考えに従って一連の適切に描画された方程式を投稿すると、あなたの努力が示され(したがって、質問が改善され、答えられる可能性が高まります)、考えられる間違いを指摘することもできます。目前の問題に関する関連情報は、回答者候補に役立つ可能性があります。
Lorenzo Donati-Codidact.org '27 / 10/27

ところで、LaTeX構文に慣れている場合、質問エディタはLaTeX数式の「ドル表記」を理解できます(オンラインヘルプを参照)。
Lorenzo Donati-Codidact.org '27 / 10/27

1
@LorenzoDonatiに感謝します。より構造化された読みやすい方法で質問を提示しようとしています。
rrz0 2018年

回答:


2

はじめに

Mとmの自由度は1つだけです。どちらも垂直方向にのみ移動できます。磁力は、質量Mではなく、磁石mに直接作用します。

絵を少しわかりやすくするために、磁石をテーブルの反対側に配置したと考えるとわかりやすいでしょう。画像はLTSPICEで描かれており、矢印はありません。したがって、矢印に最も近いのは出力ピンです。これらは水平方向にしか右を指すことができないため、画像全体が右に回転します。同じ理由で、矢印「-y」と「-F」は右を指しますが、矢印「y」と「F」を左に描画したいと思います。さらに、右側のはを読み取る必要があります。90ob1b2

これで、これは質量とそれらの間に動的要素を持つ直列接続であることは明らかです。したがって、V、y、Fを含むmの電気方程式から始めて、右から左に運動方程式を書き始めます。
その後、mとMの運動方程式を記述します
。Mは磁力の影響を受けないので、この最後の方程式はyをxの関数として与えます。これは、xを関連付ける最初の方程式で使用されます。 V.

電気

磁力と磁石の動きは、コイルの両端の電圧を介して結合されます。そして、 あり、Lがyから独立していると仮定すると、

e=αy˙,F=βi,Ve=Ri+Li˙
Ve=Vαy˙=Ri+Li˙=RβF+LβF˙

これで、(および)に関してられました。移動するオブジェクトにすべての力を加え、それらを強制的にゼロにすることで、運動方程式を(法則により)書くことができます。yFV

磁石

F+my¨+b2(y˙x˙)+k2(yx)=0
s-domainのFとyの間の上記の関係を解くことができます したがって、 磁石に加わる力の合計はゼロなので、(そして読みやすくするために、位置はまだs領域に変換されていません) sドメインに変換すると、この方程式は
Vαy˙=V(s)αsy=(R+Ls)i=(R+Ls)F/β
F=βR+Ls(V(s)αsy)
βV(s)R+LsαβR+Lssy+my¨+k2(yx)+b2(y˙x˙)=0
βV(s)R+LsαβR+Lssy+ms2y+k2(yx)+b2s(yx)=0
再グループ化した後、これは とを 分離すると、
ms2y+(b2αβR+Ls)sy+k2yb2sxk2x=βV(s)R+Ls
xy
(ms2+b2sαβsR+Ls+k2)y(b2s+k2)x=βV(s)R+Ls

移動テーブル

移動テーブルの場合、支配方程式は sに変換した後-domainこの方程式は ようになります。再グループ化すると、これは とを 分離する この方程式を書き直して、xに関してyを取得します。

Mx¨+k1x+b1x˙+k2(xy)+b2(x˙y˙)=0
Ms2x+k1x+b1sx+k2(xy)+b2s(xy)=0
b2syk2y+Ms2x+(b1+b2)sx+(k1+k2)x=0
xy
(b2s+k2)y+{Ms2+(b1+b2)s+k1+k2}x=0
y=Ms2+(b1+b2)s+k1+k2b2s+k2x

アンサンブル

置くとの間の関係に上記からの、及び磁石用: y=f(x)xyV

[(ms2+b2sαβsR+Ls+k2)Ms2+(b1+b2)s+k1+k2b2s+k2(b2s+k2)]x=βV(s)R+Ls

方程式の我々の乗算両側場合は、我々が得ますR+Ls

[{(R+Ls)(ms2+bs+k2)αβs}Ms2+(b1+b2)s+(k1+k2)b2s+k2(R+Ls)(b2s+k2)]x=βV(s)

我々乗算両側次とGETb2s+k2

[{(R+Ls)(ms2+bs+k2)αβs}{Ms2+(b1+b2)s+(k1+k2)}(R+Ls)(b2s+k2)2]x=(b2s+k2)βV(s)

目視検査から、分母の最大次数が1、分母の最大次数が伝達関数を期待できることがわかります。1つの零点が1つの極で相殺される可能性がありますが、これは推測に基づくものであり、見つけるにはさらに書き直す必要があります。x(s)/V(s)


コメントは詳細な議論のためのものではありません。この会話はチャットに移動しました
Dave Tweed
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.