Encontrar a função de transferência do sistema de amortecedor de massa de mola

9

Venho estudando o livro Ogata Modern Control Engineering e trabalhando em vários exercícios para melhorar minha compreensão dos princípios básicos de controle. Me deparei com o exemplo a seguir, que estou tentando resolver.

Preciso criar a função de transferência que modela esse gabarito de vibração. As perguntas são as seguintes:

Neste exemplo, você analisará uma plataforma de teste de vibração (Fig. 1). Este sistema consiste em uma tabela de massa M e uma bobina cuja massa é m. Um ímã permanente rigidamente fixado ao solo fornece um campo magnético constante. O movimento da bobina, 𝑦, através do campo magnético induz uma tensão na bobina que é proporcional à sua velocidade, 𝑦̇, como na Eq. 1. 𝑒 = 𝛼𝑦̇ [eq.1]

A passagem de corrente através da bobina faz com que experimente uma força magnética proporcional à corrente, como na Eq. 2. 𝐹 = 𝛽𝑖 [eq.2]

insira a descrição da imagem aqui

Pergunta: Obtenha uma função de transferência paramétrica com a saída 𝑥 para entrada 𝑉.

Algumas perguntas que estou achando difícil de responder, mas afetam todo o TF são:

  • Se K2 e B2 são comprimidos por uma distância Z (quando se move para cima
    devido à interação da bobina com o campo magnético), isso significa que k1 e b1 são estendidos pela mesma distância Z?

  • Se m(bobina) se move para cima 2 cm, a M(mesa) também se move para cima 2 cm?


O que eu preciso fazer:

  • Crie dois diagramas de corpo livre separados, um para a massa M da mesa e outro para a massa m da bobina.
  • Esboce um diagrama de circuito, incluindo a fem reversa.
  • Transforme em domínio s.
  • Resolva simultaneamente.

O que eu fiz até agora:

  • Desenhe para separar diagramas de corpo livre e extrair equações.

  • Desenhe o diagrama do circuito e extraia a equação.

  • Converta para domínio s.

Usando a função MATLAB solve, consegui obter 2 funções diferentes de transferência de 5ª ordem (uma para cada método que proponho abaixo), no entanto, não tenho certeza de qual está correta e por quê.


Sistema Geral :

Esta é uma representação esquemática de como acho que o gabarito de teste de vibração pode ser modelado, excluindo a parte elétrica.

insira a descrição da imagem aqui


Diagrama do corpo livre 1 - Tabela - Convenção ascendente

Molas k1e k2amortecedores b1e b2são modelados separadamente . Como eles não podem ser adicionados e vistos como um só, sua compactação e extensão são separadas.

A força ascendente é proveniente k2e b2está ligada à bobina. Estes estão experimentando um movimento ascendente.

insira a descrição da imagem aqui

Equação no domínio s:

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


Livre Corpo Diagrama 2 - Coil - Convenção Upward

A bobina está sofrendo uma força para cima, no entanto, a mola e o amortecedor estão retendo-a, agindo assim na direção oposta.

insira a descrição da imagem aqui

Equação no domínio s:

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


Os dois métodos diferentes são mostrados acima para o FBD da tabela levam a diferentes equações no domínio s e a diferentes funções de transferência.

Qual é o diagrama de corpo livre correto para a mesa e a bobina?

rrz0
fonte
2
Boa pergunta, mas por favor, poste uma foto onde os detalhes sejam claros sem nos forçar a clicar nela para ampliá-la. Por exemplo, esses sinais negativos são quase imperceptíveis. Além disso, a equação no canto inferior esquerdo foi parcialmente cortada. Há muito espaço livre a ser usado em sua planilha para aumentar as coisas. Existem muitos programas gratuitos de edição de imagens na Internet (por exemplo, IrfanView ou FastSstone ImageViewer), para que você também possa tirar várias fotos de suas folhas e cortar / cortar as peças necessárias para publicar boas fotos.
Lorenzo Donati - Codidact.org 27/10
@LorenzoDonati, obrigado pela sugestão, será editado imediatamente. Em relação à equação no canto inferior esquerdo, isso não interessa, pois minha preocupação é o diagrama de corpo livre. Se isso estiver correto, a equação estará correta. No entanto, tentarei editar adequadamente. Obrigado pelo seu feedback!
rrz0
Tente não fazer suposições sobre o que você fez de errado. A publicação de um conjunto de equações bem elaboradas seguindo sua linha de pensamento mostrará seus esforços (e, portanto, melhorará sua pergunta - dando mais chances de resposta) e também poderá apontar possíveis erros. Qualquer informação relevante sobre o seu problema em questão pode ser útil para o candidato a responder.
Lorenzo Donati - Codidact.org 27/10
Por outro lado, se você se sentir confortável com a sintaxe do LaTeX, o editor de perguntas poderá entender a "notação em dólar" das fórmulas do LaTeX (consulte a ajuda on-line).
Lorenzo Donati - Codidact.org 27/10
1
Obrigado @LorenzoDonati, estou tentando apresentar a pergunta de uma maneira mais estruturada e legível.
rrz0

Respostas:

2

Introdução

M e m têm apenas um grau de liberdade; ambos podem se mover apenas na vertical. A força magnética atua diretamente sobre o ímã m, não sobre a massa M.

90ob1b2

Agora está claro que esta é uma conexão em série de massas com elementos dinâmicos entre elas; portanto, começamos a escrever as equações de movimento da direita para a esquerda, começando com a equação elétrica de m primeiro, que conterá V, ye F.
Depois disso, escreveremos a equação motional para me para M.
Como M não é afetado por uma força magnética, essa última equação nos dará y em função de x, que será usado na primeira equação para relacionar x a V.

Elétrico

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

yFV

O ímã

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

A mesa móvel

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)

joe electro
fonte
Comentários não são para discussão prolongada; esta conversa foi movida para o bate-papo .
Dave Tweed