Modelos mecânicos

Ver Todos os Exercícios

Teoria

Introdução

E ai galera! Tudo tranquilo?

Agora estamos naquele momento em que muitas vezes temos um sistema real físico, como por exemplo um carro numa estrada e precisamos modelar o sistema de amortecimento dele.

Esse sistema podem ser modelado das mais diversas formas, umas mais simples e outras mais complexas, isso vai depender do que precisamos avaliar e por consequência das hipóteses que vamos adotar. Como por exemplo, o sistema carro andando na estrada pode ser modelado com um grau de liberdade.

Com dois graus de liberdade.

Esse nível de detalhe vai depender da nossa necessidade, para nossa sorte, já foi observado que os graus de liberdade podem variar, mas existem elementos primários que aparecem representando características do sistema, que são: inércia, amortecimento e rigidez.

Além desses elementos primários, quando o sistema é de segunda ordem, o que ocorre normalmente, pois é o elemento primário de inércia que ocasiona isso, podemos relacionar esses para termos mais dois novos parâmetros: critério de amortecimento e frequência natural.

Agora vamos lá entender como eles estão relacionados.

Equação do movimento

Em sistemas mecânicos, conseguimos associar um sistema real com um sistema massa-mola-amortecedor, cada uma dessas forças são conhecidas como elementos primários, que podem ser de translação ou de rotação.

Para conseguir a correlação dessas forças usamos a segunda lei de Newton.

∑ F → t = d d t m d x → t d t

Como trabalhamos com sistemas invariantes no tempo, a massa ( m ) não varia. Se ela variasse, caso do foguete, não poderíamos utilizar a transformada de Laplace e por consequência não teríamos uma função de transferência e sem ela, nada de controle clássico.

∑ F → t = m d d t d x → t d t = m d 2 x → t d t 2

Sendo F sendo a força ( N ) , m a massa ( k g ) e x o deslocamento linear ( m ) .

Caso o sistema seja rotacional, não tem problema, só alterar os termos relacionados a translação pelos seus equivalentes relacionados a rotação, repaginando então a segunda lei de Newton:

∑ τ → t = J d d t d θ → t d t = J d 2 θ → t d t 2

Sendo τ sendo o torque ( N . m ) , J o momento de inércia ( k g . m 2 ) e θ o deslocamento angular r a d .

Elementos primários de translação

Os três tipos de força encontrados ao realizarmos nossa modelagem podem ser consideradas dentro desses três tipos

1 - I n é r c i a - F = m x ¨ = m d v d t

Sendo m a massa ( k g ) , x o deslocamento linear ( m ) , v a velocidade m s .

2 - R e s t i t u i ç ã o - F = - k x 1 - x 2 = - k ∫ v 1 - v 2 d t

Sendo k a constante elástica N m , x o deslocamento linear ( m ) , v a velocidade m s .

Em sistemas mecânicos de translação, a restituição é normalmente caracterizada por uma mola. Vale lembrar quando se trata de molas que elas podem ser arranjadas em paralelo (a) e em série (b).

E isso resulta num k e q u i v a l e n t e com a seguinte relação:

a k e q u i v a l e n t e = k 1 + k 2                                                                   b k e q u i v a l e n t e = k 1 k 2 k 1 + k 2

3 - A m o r t e c i m e n t o - F = - B x ˙ 1 - x ˙ 2 = - B v 1 - v 2

Sendo B o coeficiente de atrito viscoso N . s m , x o deslocamento linear ( m ) , v a velocidade m s .

Essa relação acima ocorre, pois admitimos que o amortecimento é linear. E o amortecimento representa as forças dissipativas do nosso sistema, ou seja, as que diminuem a energia armazenada no sistema.

Assim como as molas, é possível arranjar os amortecimentos em paralelo (a) e em série (b)

E eles tem um B e q u i v a l e n t e da seguinte forma:

a B e q u i v a l e n t e = B 1 + B 2                                                                   b B e q u i v a l e n t e = B 1 B 2 B 1 + B 2

Elementos primários de rotação

Apresentados os elementos primários de translação, os de rotação são idênticos com a alteração dos termos relacionados com a translação para termos relacionados com rotação, ou seja, o torque substitui a força e o deslocamento angular substitui o deslocamento translacional.

1 - I n é r c i a - τ = J θ ¨ = J d ω d t

Sendo J o momento de inércia ( k g . m 2 ) , θ o deslocamento angular ( r a d ) , ω a velocidade ( r a d s ) .

2 - R e s t i t u i ç ã o - τ = - k θ 1 - θ 2 = - k ∫ ω 1 - ω 2 d t

Sendo k a constante de mola torcional ( N . m r a d ) , θ o deslocamento angular ( r a d ) , ω a velocidade ( r a d s ) .

A relação para molas torcionais em série e paralelo são as mesmas do sistema de translação.

a k e q u i v a l e n t e = k 1 + k 2                                                                   b k e q u i v a l e n t e = k 1 k 2 k 1 + k 2

3 - A m o r t e c i m e n t o - τ = - D θ ˙ 1 - θ ˙ 2 = - D ω 1 - ω 2

Sendo D o coeficiente de atrito viscoso ( N . m r a d ) , θ o deslocamento angular ( r a d ) , ω a velocidade ( r a d s ) .

A relação para molas torcionais em série e paralelo são as mesmas do sistema de translação.

a D e q u i v a l e n t e = D 1 + D 2                                                                   b D e q u i v a l e n t e = D 1 D 2 D 1 + D 2

Equação do movimento

Elementos primários de translação

Elementos primários de rotação

Exercícios Resolvidos

Exercício Resolvido #1

Para o sistema mecânico da figura abaixo, sendo f ( t ) , o sinal de entrada e x ( t ) , o sinal de saída, obter a função de transferência do sistema.

Passo 1

Para esse exercício, precisaremos obter uma relação de equilíbrio entre as diferentes forças envolvidas.

F r = ∑ F

Temos três forças agindo sobre o nosso objeto: a força f ( t ) , a força da mola e a força de atrito.

A força da mola será:

F m o l a = K x ( t )

E a força de atrito:

F a t r i t o = B d x t d t

Passo 2

Vamos arbitrar que a força é positiva para baixo e negativa para cima. Sendo assim:

F r = f t - k x t - B d x t d t

A força resultante sobre o corpo, pode ser escrita usando-se a Segunda Lei de Newton:

F r = M a = M   d 2 x t d t 2

Sendo a a aceleração do corpor, e aceleração é a derivada dupla da posição. Substituindo na equação de equilíbrio:

M   d 2 x t d t 2 = f t - k x t - B d x t d t

Passo 3

Como o exercício pede a função de transferência do sistema, precisamos tirar a transformada de Laplace da equação acima e encontrar a relação entre a força f ( t ) e a posição x ( t ) .

Aplicando Laplace:

M s 2 X s = F s - k X s - B s X ( s )

Essa parte da transformada de Laplace é bem importante, se você estiver com dúvidas nesse passo seria legal dar uma revisada no tópico de transformadas de Laplace.

Rearrajando:

F s = M s 2 X s + B s X s + k X ( s )

F s = M s 2 + B s + k X ( s )

Logo:

X s F s = 1 M s 2 + B s + k

Resposta

X s F s = 1 M s 2 + B s + k

Exercício Resolvido #2

Escreva as equações de movimento para a velocidade do carro conforme mostrado na figura abaixo assumindo que o motor fornece uma força u . Obetnha a transformada de Laplace da equação diferencial resultante e encontre a função de transferência entre a entrada u e a saída v .

Passo 1

Precisamos encontrar a relação entre a força que o motor aplica no carro e a sua velocidade.

Nesse caso, precisamos considerar duas forças agindo sobre o carro, a força u fornececida pelo motor e a força de atrito.

F r = ∑ F

F r = u t - b d x t d t

Usando a segunda lei de newton:

F r = M a = M   d 2 x t d t 2

Sendo a a aceleração do corpo, e aceleração é a derivada dupla da posição. Substituindo na equação de equilíbrio:

M   d 2 x t d t 2 = u t - b d x t d t

Passo 2

Por enquanto, temos uma EDO que relaciona a força fornecida pelo motor com a posição do carro, entretanto, precisamos encontrar a relação com a velocidade e não com a posição. Como a velocidade é a derivada da posição, podemos usar essa relação:

v t = d x t d t

Logo:

M d d t d x t d t = u t - b d x t d t

M d v t d t = u t - b v ( t )

Rearranjando:

u t = M d v t d t + b v ( t )

Passo 3

Obtivemos então a EDO que relaciona a velocidade do carro com a força fornecida pelo motor:

u t = M d v t d t + b v ( t )

O próximo passo agora é aplicar a transformada de Laplace:

U s = M s V s + b V ( s )

Essa parte da transformada de Laplace é bem importante, se você estiver com dúvidas nesse passo seria legal dar uma revisada no tópico de transformadas de Laplace.

Rearrajando:

U s = M s + b V ( s )

Logo:

V s U s = 1 M s + b

Resposta

V s U s = 1 M s + b

Exercício Resolvido #3

Determine a amplitude e a frequência de oscilação da resposta a um impulso de torque no pêndulo de torção sem atrito representado abaixo, composto de um corpo rígido de inércia rotacional 𝐽 suspenso por um fio com rigidez rotacional 𝐾.

Passo 1

Nesse exercício temos um sistema mecânico rotacional.

Precisamos encontrar uma expressão para o ângulo de rotação θ ( t ) , para isso, vmos escrever a equação de equilíbrio de torque:

τ r = ∑ τ

O torque resultante é a soma de todos os torques, e temos que:

τ r = J d θ 2 d t 2

Logo:

J d θ 2 d t 2 = ∑ τ

Nesse caso, não temos atrito, mas temos um impulso de torque, T t , e também temos um fio com rigidez rotacional, K , que atrapalha a rotação.

J d θ 2 d t 2 = T t - K θ ( t )

Passo 2

Devemos agora aplicar transformada de Laplace na nossa EDO:

J s 2 θ s = 1 - K θ ( s )

Essa parte da transformada de Laplace é bem importante, se você estiver com dúvidas nesse passo seria legal dar uma revisada no tópico de transformadas de Laplace.

A transformada de Laplace da função de torque foi 1 devido ao fato de ser uma função impulso. Funções impulso sempre vão ter transformada de Laplace igual a 1.

Rearranjando:

J s 2 θ s + K θ s = 1

J s 2 + K θ s = 1

θ s = 1 J s 2 + K

Passo 3

Já obtivemos a expressão para o ângulo θ no domínio s :

θ s = 1 J s 2 + 1

Se manipularmos um pouco teremos:

θ s = 1 / J s 2 + K J = 1 K J   K J s 2 + K J

Fizemod essa manipulação para conseguir aplicar a inversa da transformada de Laplace. Se olharmos uma tabela com transformadas de Laplace, veremos que essa função ficou com o mesmo formato da transformada da função seno. Fazedno a transformada inversa teremos:

θ t = 1 K J sen ⁡ K J t

Passo 4

O exercício pede a amplitude e frequência de oscilação.

A amplitude é a constante que multiplica a função seno:

a m p l i t u d e =   1 K J

E a frequência é o argumento da função:

f r e q u ê n c i a = K J  

Resposta

a m p l i t u d e =   1 K J

f r e q u ê n c i a = K J

Exercício Resolvido #4

Encontre a função de transferência X s U s do sistema de massa molo abaixo:

Passo 1

Aqui nessa questão temos um corpor m e três forças agindo sobre ele: a força u ( t ) , a força da mola e a força proveniente do sistema de amortecimento. Sendo assim:

F r = ∑ F

F r = u t - k x t - b d x t d t

Usando a segunda lei de newton:

F r = M a = M   d 2 x t d t 2

Sendo M a massa do corpo e a a aceleração.

Aceleração é a derivada dupla da posição. Substituindo na equação de equilíbrio:

M   d 2 x t d t 2 = u t - k x t - b d x t d t

Passo 2

Vamos rearranjar para deixar tudo que tem a a posição x ( t ) do lado esquerdo da equação:

M   d 2 x t d t 2 + k x t + b d x t d t = u ( t )

Já temos a nossa EDO no jeitinho. O próximo passo é aplicar a transformada de Laplace nela para que possamos determinar a função de transferência:

Aplicando Laplace:

M s 2 X s + k X s + b s X s = U ( s )

Só lembrando: se você não está se lembrando como fazer esse passo, seria legal dar uma revisa em transformadas de Laplace. Esse passo é fundamental e precisa estar na ponta da língua hein.

Passo 3

Colocando X ( s ) em evidência do lado esquerdo da equação:

M s 2 + k + b s X s = U ( s )

Rearranjando:

M s U s = 1 M s 2 + b s + k

Resposta

M s U s = 1 M s 2 + b s + k

Exercício Resolvido #5

Encontre as funções de transferência Y 1 s U s e Y 2 s U s do sistema abaixo:

Passo 1

Vamos escrever a equação de equilíbrio de forças para cada um dos corpos.

Para o bloco 1:

m 1 d 2 y 1 d t 2 = - b d y 1 d t - k y 1 - y 2

O primeiro sofre influência da força de amortecimento e da força da mola.

Para o bloco 2:

m 2 d 2 y 2 d t 2 = u ( t ) - k y 2 - y 1

O segundo bloco sente o efeito apenas da força da mola e da força u ( t )

Vamos rearranjar essas expressões para ficar melhor para a gente trabalhar:

m 1 d 2 y 1 d t 2 + b d y 1 d t + k y 1 = k y 2

m 2 d 2 y 2 d t 2 + k y 2 = u t + k y 1

Passo 2

As EDOs já estão no jeito. Vamos aplicar Laplace nelas então:

m 1 s 2 Y 1 s + b s Y 1 s + k Y 1 s = k Y 2 s                       ( 1 )

m 2 s 2 Y 2 s + k Y 2 s = U s + k Y 1 S                             ( 2 )

Só lembrando: se você não está se lembrando como fazer esse passo, seria legal dar uma revisa em transformadas de Laplace. Esse passo é fundamental e precisa estar na ponta da língua hein.

Passo 3

Vamos isolar Y 1 ( s ) em ( 1 ) e substituir em ( 2 ) :

Y 1 s = k Y 2 s m 1 s 2 + b s + k         ( 3 )

Substituindo:

m 2 s 2 Y 2 s + k Y 2 s = U s + k 2 Y 2 s m 1 s 2 + b s + k  

Rearranjando:

m 2 s 2 + k - k 2 m 1 s 2 + b s + k Y 2 s = U ( s )

m 2 s 2 + k m 1 s 2 + b s + k - k 2 Y 2 s = m 1 s 2 + b s + k U ( s )

Por fim:

Y 2 s U s = m 1 s 2 + b s + k m 2 s 2 + k m 1 s 2 + b s + k - k 2             ( 4 )

Substituindo 4 em ( 3 ) :

Y 1 s = k U s m 2 s 2 + k m 1 s 2 + b s + k - k 2

Y 1 s U s = k m 2 s 2 + k m 1 s 2 + b s + k - k 2

Resposta

Y 1 s U s = k m 2 s 2 + k m 1 s 2 + b s + k - k 2

Y 2 s U s = m 1 s 2 + b s + k m 2 s 2 + k m 1 s 2 + b s + k - k 2

Exercício Resolvido #6

Encontre as função de transferência V 1 s R s para o sistema mecânico abaixo:

Passo 1

Vamos começar escrevendo a equação de equilíbrio para o bloco 1:

F r = ∑ F

m 1   d 2 x 1 d t 2 = r t - B 2 d x 2 d t - B 1 d x 1 t d t - d x 2 t d t

Como queremos uma expressão que relacione a velocidade com a força r ( t ) :

v t = d x t d t

m 1 d v 1 t d t = r t - B 2 v 1 t - B 1 v 1 t - v 2 t

Rearranjando:

m 1 d v 1 t d t + B 1 + B 2 v 1 t = r t + B 1 v 2 t                           ( 1 )

Analogamente para o segundo bloco teremos:

m 2   d 2 x 2 d t 2 = - k x 2 t - B 1 d x 2 t d t - d x 1 t d t

Fazendo:

v t = d x t d t

E, consequentemente:

x t = ∫ v t d t

m 2 d v 2 t d t = - k ∫ v 2 t d t - B 1 v 2 ( t ) - v 1 ( t )

Rearranjando:

m 2 d v 2 t d t + k ∫ v 2 t d t + B 1 v 2 t = B 1 v 1 t                     ( 2 )

Passo 2

Agora que já fizemos a modelagem do sistema mecânico, devemos aplicar a transformada de Laplace nas expressões (1) e (2).

Para a expressão (1):

m 1 s V 1 s + B 1 + B 2 V 1 s = R s + B 1 V 2 ( s )

Rearranjando:

m 1 s + B 1 + B 2 V 1 s = R s + B 1 V 2 s                           ( 3 )

Aplicando Laplace na expressão (2):

m 2 s V 2 s + k V 2 s s + B 1 V 2 s = B 1 V 1 s                         ( 4 )

Passo 3

Para terminar, vamos isolar V 2 ( s ) em (4) e substituir em (3):

m 2 s + k s + B 1 V 2 s = B 1 V 1 ( s )

V 2 s = B 1 V 1 s   m 2 s + k s + B 1

Substituindo em (3):

m 1 s + B 1 + B 2 V 1 s = R s + B 1 2 V 1 s   m 2 s + k s + B 1          

Logo:

m 1 s + B 1 + B 2 - B 1 2   m 2 s + k s + B 1 V 1 s = R ( s )

m 1 s + B 1 + B 2 m 2 s + k s + B 1 - B 1 2 V 1 s = m 2 s + k s + B 1 R ( s )  

Por fim:

V 1 s R s = m 2 s + k s + B 1 m 1 s + B 1 + B 2 m 2 s + k s + B 1 - B 1 2

Resposta

V 1 s R s = m 2 s + k s + B 1 m 1 s + B 1 + B 2 m 2 s + k s + B 1 - B 1 2

Exercício Resolvido #7

Encontre as funções de transferência X 1 s U s e X 2 s U s para o sistema mecânico abaixo:

Passo 1

Precisamos então encontrar a função de transferência que relaciona a posição de cada um dos blocos com a força u aplicada.

F r = ∑ F

Para o primeiro bloco teremos:

m 1   d 2 x 1 d t 2 = u t - k 1 x 1 t - k 2 x 1 t - x 2 t - B d x 1 t d t - d x 2 t d t

A força resultante ( m a s s a   x   a c e l e r a ç ã o ) é igual a soma de todas as forças atuantes no bloco. Para o primeiro bloco temos a força de entrada u ( t ) , a força proveniente das molas 1 e 2 e a força de atrito.

Analogamente para o segundo bloco teremos:

m 2   d 2 x 2 d t 2 = - k 3 x 2 t - k 2 x 2 t - x 1 t - B d x 2 t d t - d x 1 t d t

Vamos tentar reorganizar essas equações para conseguir trabalhar com elas um pouco mais facilmente. Para a primeira equação, vamos deixar tudo que tem x 1 do lado direito e todo o resto do lado esquerdo:

m 1   d 2 x 1 d t 2 + B d x 1 t d t + k 1 + k 2 x 1 t = B d x 2 t d t + k 2 x 2 t + u ( t )

Para a segunda equação, vamos deixar tudo que tem x 2 do lado esquerdo:

m 2   d 2 x 2 d t 2 + B d x 2 t d t + k 2 + k 3 x 2 t = B d x 1 d t + k 2 x 1 ( t )

Passo 2

Obtivemos então as EDOs que regem o sistema:

m 1   d 2 x 1 d t 2 + B d x 1 t d t + k 1 + k 2 x 1 t = B d x 2 t d t + k 2 x 2 t + u ( t )

m 2   d 2 x 2 d t 2 + B d x 2 t d t + k 2 + k 3 x 2 t = B d x 1 d t + k 2 x 1 ( t )

O próximo passo agora é aplicar a transformada de Laplace. Para a primeira equação:

m 1 s 2 X 1 s + B s X 1 s + k 1 + k 2 X 1 s = B s X 2 s + k 2 X 2 s + U ( s )

E para a segunda:

m 2 s 2 X 2 s + B s X 2 s + k 2 + k 3 X 2 s = B s X 1 s + k 2 X 1 s

Rearranjando as equações obtidas:

m 1 s 2 + B s + k 1 + k 2 X 1 s = B s + k 2 X 2 s + U s                               ( 1 )

m 2 s 2 + B s + k 2 + k 3 X 2 s = B s + k 2 X 1 s                                                         ( 2 )

Passo 3

Para achar as expressões de X 1 s U s e X 2 s U s precisaremos fazer um pouco de algebrismo. Vamos então isolar X 2 ( s ) na equação (2):

X 2 s = B s + k 2 X 1 s m 2 s 2 + B s + k 2 + k 3

Substituindo na equação (1), teremos:

m 1 s 2 + B s + k 1 + k 2 X 1 s = B s + k 2 2 X 1 s m 2 s 2 + B s + k 2 + k 3   + U s  

Rearranjando:

m 1 s 2 + B s + k 1 + k 2 - B s + k 2 2 m 2 s 2 + B s + k 2 + k 3   X 1 s = m 2 s 2 + B s + k 2 + k 3 U ( s )

m 1 s 2 + B s + k 1 + k 2 m 2 s 2 + B s + k 2 + k 3 - B s + k 2 2 X 1 s = ( m 2 s 2 + B s + k 2 + k 3 ) U ( s )

Finalmente:

X 1 s U s = ( m 2 s 2 + B s + k 2 + k 3 ) m 1 s 2 + B s + k 1 + k 2 m 2 s 2 + B s + k 2 + k 3 - B s + k 2 2

Já achamos X 1 s U s , agora só falta X 2 s U s .

Substituindo a equação para X 1 ( s ) U s que acabamos de encontrar, na equação (2):

m 2 s 2 + B s + k 2 + k 3 X 2 s = ( m 2 s 2 + B s + k 2 + k 3 ) B s + k 2 U ( s ) m 1 s 2 + B s + k 1 + k 2 m 2 s 2 + B s + k 2 + k 3 - B s + k 2 2

Podemos cancelar o termo m 2 s 2 + B s + k 2 + k 3 que aparece dos dois lados da equação:

X 2 s = B s + k 2 U ( s ) m 1 s 2 + B s + k 1 + k 2 m 2 s 2 + B s + k 2 + k 3 - B s + k 2 2

Finalmente:

X 2 s U s = B s + k 2 m 1 s 2 + B s + k 1 + k 2 m 2 s 2 + B s + k 2 + k 3 - B s + k 2 2

Resposta

X 1 s U s = ( m 2 s 2 + B s + k 2 + k 3 ) m 1 s 2 + B s + k 1 + k 2 m 2 s 2 + B s + k 2 + k 3 - B s + k 2 2

X 2 s U s = B s + k 2 m 1 s 2 + B s + k 1 + k 2 m 2 s 2 + B s + k 2 + k 3 - B s + k 2 2

Exercício Resolvido #8

Obtenha a função de transferência, θ s M s que representa o sistema mecânico abaixo:

Passo 1

Nesse exercício temos um sistema mecânico rotacional.

τ r = ∑ τ

O torque resultante é a soma de todos os torques, e temos que:

τ r = J d θ 2 d t 2

Logo:

J d θ 2 d t 2 = ∑ τ

Nesse caso temos um momento angular M ( t ) , mas também temos um momento de atrito e de resistência do fio que “atrapalham” a rotação:

J d θ 2 d t 2 = M t - B d θ d t - K θ ( t )

Passo 2

Para achar a função de transferência agora só precisamos aplicar a transformada de Laplace e depois manipular a equação obtida. Aplicando Laplace:

J s 2 θ s = M s - B s θ s - K θ ( s )

Rarranjando:

J s 2 + B s + K θ s = M s

θ s M s = 1 J s 2 + B s + K  

Resposta

θ s M s = 1 J s 2 + B s + K