Recherche de la fonction de transfert du système d'amortisseur de masse à ressort

9

J'ai parcouru le livre Ogata Modern Control Engineering et travaillé à travers plusieurs exercices pour améliorer ma compréhension des principes de base du contrôle. Je suis tombé sur l'exemple suivant que j'ai du mal à résoudre.

Je dois trouver la fonction de transfert qui modélise ce gabarit de vibration. Les questions sont les suivantes:

Dans cet exemple, vous allez analyser un banc d'essai de vibration (Fig. 1). Ce système est constitué d'une table de masse M et d'une bobine dont la masse est m. Un aimant permanent fixé rigidement au sol fournit un champ magnétique constant. Le mouvement de la bobine, 𝑦, à travers le champ magnétique induit une tension dans la bobine qui est proportionnelle à sa vitesse, 𝑦̇, comme dans Eq. 1. 𝑒 = 𝛼𝑦̇ [éq.1]

Le passage du courant à travers la bobine lui fait subir une force magnétique proportionnelle au courant comme dans l'équation. 2. 𝐹 = 𝛽𝑖 [éq.2]

entrez la description de l'image ici

Question: Obtenez une fonction de transfert paramétrique avec la sortie 𝑥 vers l'entrée 𝑉.

Certaines questions auxquelles j'ai du mal à répondre mais qui affectent l'ensemble du TF sont:

  • Si K2 et B2 sont compressés d'une distance Z, (en se déplaçant vers le haut en
    raison de l'interaction de la bobine avec le champ magnétique), cela signifie-t-il que k1 et b1 sont étendus de la même distance Z?

  • Si m(la bobine) se déplace de 2 cm vers le haut, M(la table) se déplace-t-elle également de 2 cm vers le haut?


Ce que je dois faire:

  • Trouvez deux diagrammes de corps libres séparés, un pour la masse M de la table et un pour la masse m de la bobine.
  • Esquissez un schéma de circuit, y compris la FEM arrière.
  • Transformez-vous en domaine s.
  • Résolvez simultanément.

Ce que j'ai fait jusqu'à présent:

  • Dessinez pour séparer les diagrammes de corps libres et extraire les équations.

  • Dessinez le schéma de circuit et extrayez l'équation.

  • Convertissez en s-domaine.

En utilisant la fonction MATLAB, solvej'ai réussi à obtenir 2 fonctions de transfert de 5e ordre différentes (une pour chaque méthode que je propose ci-dessous), mais je ne sais pas laquelle est correcte et pourquoi.


Système global :

Ceci est une représentation schématique de la façon dont je pense que le gabarit d'essai de vibration peut être modélisé, à l'exclusion de la partie électrique.

entrez la description de l'image ici


Diagramme de corps libre 1 - Tableau - Convention ascendante

Ressorts k1et k2et des amortisseurs b1et b2sont modélisés séparément . Puisqu'ils ne peuvent pas être additionnés et considérés comme un, leur compression et leur extension sont distinctes.

La force ascendante vient de k2et b2qui est attachée à la bobine. Ceux-ci connaissent un mouvement ascendant.

entrez la description de l'image ici

Équation dans le domaine s:

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


Diagramme de corps libre 2 - Bobine - Convention ascendante

La bobine subit une force vers le haut, mais le ressort et l'amortisseur la retiennent, agissant ainsi dans la direction opposée.

entrez la description de l'image ici

Équation dans le domaine s:

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


Les deux méthodes différentes sont indiquées ci-dessus pour le FBD du tableau conduisent à différentes équations dans le domaine s et différentes fonctions de transfert.

Quel est le bon diagramme de corps libre pour la table et la bobine?

rrz0
la source
2
Jolie question, mais merci de poster une photo où les détails sont clairs sans nous forcer à cliquer dessus pour l'agrandir. Par exemple, ces signes négatifs sont à peine perceptibles. De plus, l'équation en bas à gauche a été partiellement recadrée. Il y a beaucoup d'espace libre à utiliser sur votre feuille pour être utilisé pour agrandir les choses. Il existe de nombreux programmes d'édition d'images gratuits sur Internet (par exemple, IrfanView ou FastSstone ImageViewer), vous pouvez donc prendre également plusieurs photos de vos feuilles et couper / recadrer les pièces dont vous avez besoin pour publier de belles photos.
Lorenzo Donati - Codidact.org
@LorenzoDonati, merci pour la suggestion, modifiera immédiatement. En ce qui concerne l'équation en bas à gauche, cela n'a pas d'intérêt puisque ma préoccupation est le diagramme du corps libre. Si c'est correct, alors l'équation sera correcte. Cependant, je vais essayer de modifier en conséquence. Merci pour vos commentaires.
rrz0
Essayez de ne pas faire d'hypothèses sur ce que vous avez fait de mal. La publication d'un ensemble d'équations bien dessinées suivant votre ligne de pensée montrera vos efforts (et améliorera ainsi votre question - en lui donnant plus de chances d'y répondre) et peut également signaler d'éventuelles erreurs. Toute information pertinente concernant votre problème pourrait être utile pour le répondeur potentiel.
Lorenzo Donati - Codidact.org
BTW, si vous êtes à l'aise avec la syntaxe LaTeX, l'éditeur de questions peut comprendre la "notation en dollars" des formules LaTeX (voir l'aide en ligne).
Lorenzo Donati - Codidact.org
1
Merci @LorenzoDonati, j'essaie de présenter la question de manière plus structurée et lisible.
rrz0

Réponses:

2

Intro

M et m n'ont qu'un seul degré de liberté; les deux ne peuvent se déplacer que verticalement. La force magnétique agit directement sur l'aimant m, pas sur la masse M.

90ob1b2

Maintenant, il est clair qu'il s'agit d'une connexion en série de masses avec des éléments dynamiques entre elles, nous commençons donc à écrire les équations de mouvement de droite à gauche, en commençant par l'équation électrique de m en premier, qui contiendra V, y et F.
Après cela, nous écrirons l'équation de mouvement pour m et pour M.
Comme M n'est pas affecté par une force magnétique, cette dernière équation nous donnera y en fonction de x, qui sera utilisé dans la première équation pour relier x à V.

Électrique

e=αy˙,F=βje,V-e=Rje+Lje˙
V-e=V-αy˙=Rje+Lje˙=RβF+LβF˙

yFV

L'aimant

F+my¨+b2(y˙-X˙)+k2(y-X)=0
V-αy˙=V(s)-αsy=(R+Ls)je=(R+Ls)F/β
F=βR+Ls(V(s)-αsy)
βV(s)R+Ls-αβR+Lssy+my¨+k2(y-X)+b2(y˙-X˙)=0
βV(s)R+Ls-αβR+Lssy+ms2y+k2(y-X)+b2s(y-X)=0
ms2y+(b2-αβR+Ls)sy+k2y-b2sX-k2X=-βV(s)R+Ls
Xy
(ms2+b2s-αβsR+Ls+k2)y-(b2s+k2)X=-βV(s)R+Ls

La table mobile

MX¨+k1X+b1X˙+k2(X-y)+b2(X˙-y˙)=0
Ms2X+k1X+b1sX+k2(X-y)+b2s(X-y)=0
-b2sy-k2y+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

Ensemble

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)

b2s+k2

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

X(s)/V(s)

joe electro
la source
Les commentaires ne sont pas pour une discussion approfondie; cette conversation a été déplacée vers le chat .
Dave Tweed