Scribd
Carregar
5 visualizações4 páginas

PHD3002 Aula 05

Enviado por

vinicius.lino.poli
Direitos autorais
© © All Rights Reserved
Levamos muito a sério os direitos de conteúdo. Se você suspeita que este conteúdo é seu, reivindique-o aqui.
Formatos disponíveis
Baixe no formato PDF, TXT ou leia on-line no Scribd
5 visualizações4 páginas

PHD3002 Aula 05

Enviado por

vinicius.lino.poli
Direitos autorais
© © All Rights Reserved
Levamos muito a sério os direitos de conteúdo. Se você suspeita que este conteúdo é seu, reivindique-o aqui.
Formatos disponíveis
Baixe no formato PDF, TXT ou leia on-line no Scribd
Você está na página 1/ 4

02/04/2018

Modelagem em
Engenharia C & A
Aula 5- Sistemas Dinâmicos
Tanque oscilador

Controle de fluxo
• Contador
For i=1 to x step m
....
Next i

• Loop
While
.....
Wend

1
02/04/2018

Controle de Fluxo
• IF ... Then......

• Goto label

Problema do Tanque
Oscilante

Pontos da trajetória do jato 30


x dos pontos 0.1 m
Ho 0.5 m
h 0.5 m
T (frequência de oscilação 20 s
t 0.02 s
Cv (coeficiente de velocidade do jato 0,98

2
02/04/2018

Modelo matemático
Equação da aceleração total em x:

DV dVx dVx
  Vx 0
Dt dt Dx
Aproximação numérica da velocidade pontual em x:

(Vi t  Vi t1 ) t t
Vi t 1  Vi t  (Vi  Vi t1 )
2 x
Equação da velocidade do jato:

Vx  Cv 2 gh
Equação oscilante sobre o orifício:
2
h  ho  h cos( t )
T
Velocidade pontual em y:

Vy  Vyo  gt
Posição temporal em y:

1
y  yo  Vy  t  g t2
2

Gráfico

1 2 3 4 5 2
-11 Pontos 30
-10 Dx 0.1 m 1.5
-9 ho 0.5 m
Avança 1
-8 ah 0.5 m
-7 T 20 s roda 0.5
-6 dt 0.02 s
para
-5 Ang 45 Graus -0.1 00.027468
-4 0-0.50.027468
y

-3 k 2.967729 m
-2 h 0.027468 m -1
-1 t 8.94
-1.5
0 x V y v t+1
1 0 0.347798 0 0.347798 0.027468 -2
0.1 0.420291 -0.17768 0.414723 0.027468
0.2 0.477178 -0.66166 0.472073 -2.5
0.3 0.525262 -1.30004 0.520442
-3
0.4 0.567551 -2.0364 0.56293
-0.1 0.4 0.9 1.4 1.9
0.5 0.60567 -2.84277 0.601198
x
0.6 0.640607 -3.70287 0.636253
0.7 0.673013 -4.60625 0.668756
0.8 0.703346 -5.54572 0.699171
0.9 0.731938 -6.5161 0.727834
1 0.759042 -7.51348 0.755001
1.1 0.784855 -8.53486 0.78087
1.2 0.809533 -9.57785 0.805599
1.3 0.833203 -10.6405 0.829315
1.4 0.85597 -11.7213 0.852124
1.5 0.87792 -12.819 0.874114
1.6 0.899129 -13.9323 0.89536

3
02/04/2018

'Interpola/Extrapola a partir de valores em um par de colunas


'
Function Pint(xf As Range, yf As Range, xint As Variant) As Variant
Dim j As Integer
For j = 2 To xf.Rows.Count - 1
If xint <= xf.Cells(j, 1) Then Exit For
Next j
Pint = yf.Cells(j - 1, 1) + (yf.Cells(j, 1) - yf.Cells(j - 1, 1)) / (xf.Cells(j, 1) -
xf.Cells(j - 1, 1)) * (xint - xf.Cells(j - 1, 1))
End Function
'------------------------------
'------------------------------
'Interpola valores em um par de linhas
'
Function pintR(xf As Range, yf As Range, xint As Variant) As Variant
Dim j As Integer
For j = 1 To xf.Columns.Count - 1
If xint <= xf.Cells(1, j) Then Exit For
Next j
pintR = yf.Cells(1, j - 1) + (yf.Cells(1, j) - yf.Cells(1, j - 1)) / (xf.Cells(1, j) -
xf.Cells(1, j - 1)) * (xint - xf.Cells(1, j - 1))
End Function
'------------------------------

'Esta rotina atualiza o gráfico


Sub fncAtualizaGraficos()
Dim myChart As ChartObject
Dim myCharts As ChartObjects
Dim myChartname As String

Set myCharts = ActiveSheet.ChartObjects

For Each myChart In myCharts


myChartname = myChart.Name
ActiveSheet.ChartObjects(myChartname).Chart.Refresh
Next

End Sub

Você também pode gostar