Encontrar la función de transferencia del sistema Spring Mass Damper


9

He estado revisando el libro Ogata Modern Control Engineering y trabajando en varios ejercicios para mejorar mi comprensión de los principios básicos de control. Encontré el siguiente ejemplo que estoy luchando por resolver.

Necesito idear la función de transferencia que modela esta plantilla de vibración. Las preguntas son las siguientes:

En este ejemplo, analizará una plataforma de prueba de vibración (Fig. 1). Este sistema consiste en una tabla de masa M y una bobina cuya masa es m. Un imán permanente rígidamente unido al suelo proporciona un campo magnético estable. El movimiento de la bobina, 𝑦, a través del campo magnético induce un voltaje en la bobina que es proporcional a su velocidad, 𝑦̇, como en la ecuación. 1. 𝑒 = 𝛼𝑦̇ [eq.1]

El paso de corriente a través de la bobina hace que experimente una fuerza magnética proporcional a la corriente como en la ecuación. 2. 𝐹 = 𝛽𝑖 [eq.2]

ingrese la descripción de la imagen aquí

Pregunta: Obtenga una función de transferencia paramétrica con salida 𝑥 a entrada 𝑉.

Algunas preguntas que encuentro difíciles de responder pero que afectan a todo el TF son:

  • Si K2 y B2 están comprimidos por una distancia Z (cuando se mueve hacia arriba
    debido a que la bobina interactúa con el campo magnético), ¿significa esto que k1 y b1 se extienden en la misma distancia Z?

  • Si m(bobina) se mueve hacia arriba 2 cm, ¿se mueve M(mesa) también hacia arriba 2 cm?


Lo que necesito hacer:

  • Cree dos diagramas de cuerpo libre separados, uno para la masa M de la tabla y otro para la masa m de la bobina.
  • Dibuje un diagrama de circuito que incluya la fem posterior.
  • Transformar a dominio s.
  • Resuelve simultáneamente.

Lo que he hecho hasta ahora:

  • Dibuje para separar diagramas de cuerpo libre y extraer ecuaciones.

  • Dibuja el diagrama del circuito y extrae la ecuación.

  • Convierte a dominio s.

Usando la función MATLAB solve, logré obtener 2 funciones de transferencia de quinto orden diferentes (una para cada método que propongo a continuación), sin embargo, no estoy seguro de cuál es la correcta y por qué.


Sistema en general :

Esta es una representación esquemática de cómo creo que se puede modelar la plantilla de prueba de vibración, excluyendo la parte eléctrica.

ingrese la descripción de la imagen aquí


Diagrama de cuerpo libre 1 - Tabla - Convención ascendente

Muelles k1y k2amortiguadores b1y b2se modelan por separado . Como no pueden agregarse y verse como uno solo, su compresión y extensión son independientes.

La fuerza ascendente proviene k2y b2se unen a la bobina. Estos están experimentando un movimiento ascendente.

ingrese la descripción de la imagen aquí

Ecuación en el dominio s:

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


Diagrama de cuerpo libre 2 - Bobina - Convención ascendente

La bobina está experimentando una fuerza hacia arriba, sin embargo, el resorte y el amortiguador lo retienen, actuando en la dirección opuesta.

ingrese la descripción de la imagen aquí

Ecuación en el dominio s:

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


Los dos métodos diferentes que se muestran arriba para el FBD de la tabla conducen a diferentes ecuaciones en el dominio s y diferentes funciones de transferencia.

¿Cuál es el diagrama de cuerpo libre correcto para la mesa y la bobina?


2
Buena pregunta, pero, por favor, publique una foto donde los detalles sean claros sin obligarnos a hacer clic en ella para ampliarla. Por ejemplo, esos signos menos son apenas perceptibles. Además, la ecuación en la parte inferior izquierda se ha recortado parcialmente. Hay mucho espacio libre para usar en tu hoja para agrandar las cosas. Hay muchos programas gratuitos de edición de imágenes en Internet (por ejemplo, IrfanView o FastSstone ImageViewer), por lo que también puede tomar varias fotos de sus hojas y cortar / recortar las partes que necesita para publicar buenas fotos.
Lorenzo Donati - Codidact.org

@LorenzoDonati, gracias por la sugerencia, se editará de inmediato. Con respecto a la ecuación en la parte inferior izquierda, eso no es interesante ya que mi preocupación es el diagrama de cuerpo libre. Si eso es correcto, entonces la ecuación será correcta. Sin embargo, intentaré editar en consecuencia. Gracias por sus comentarios.
rrz0

Trate de no hacer suposiciones sobre lo que hizo mal. Publicar un conjunto de ecuaciones bien dibujadas siguiendo su línea de pensamiento mostrará sus esfuerzos (y por lo tanto mejorará su pregunta, dándole más oportunidades de ser respondido) y también puede señalar posibles errores. Cualquier información relevante sobre su problema en cuestión podría ser útil para el posible respondedor.
Lorenzo Donati - Codidact.org

Por cierto, si se siente cómodo con la sintaxis de LaTeX, el editor de preguntas puede comprender la "notación en dólares" de las fórmulas de LaTeX (consulte la ayuda en línea).
Lorenzo Donati - Codidact.org

1
Gracias @LorenzoDonati, estoy tratando de presentar la pregunta de una manera más estructurada y legible.
rrz0

Respuestas:


2

Introducción

M ym tienen solo un grado de libertad; ambos pueden moverse verticalmente solamente. La fuerza magnética actúa directamente sobre el imán m, no sobre la masa M.

90ob1b2

Ahora está claro que esta es una conexión en serie de masas con elementos dinámicos entre ellas, por lo que comenzamos a escribir las ecuaciones de movimiento de derecha a izquierda, comenzando con la ecuación eléctrica para m primero, que contendrá V, y y F.
Después de eso, escribiremos la ecuación de movimiento para m y para M.
Como M no se ve afectada por una fuerza magnética, esta última ecuación nos dará y como una función de x, que se usará en la primera ecuación para relacionar x con V.

Eléctrico

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

yFV

El imán

F+my¨+b2(y˙x˙)+k2(yx)=0
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

La mesa móvil

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

Conjunto

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)


Los comentarios no son para discusión extendida; Esta conversación se ha movido al chat .
Dave Tweed
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.