08. SIMULANDO MODELOS QUANTITATIVOS
OBJETIVOS:
1. Escrever uma equação para a variação da quantidade, em um
tanque de depósito onde existe uma drenagem;
2. Usando uma tabela, calcular a quantidade e fluxo por hora,
plotando os resultados em um gráfico de quantidade versus tempo,
explicando como estes cálculos são capazes de simular processos
representados por modelos;
3. Representar a simulação com um fluxograma;
4. Escrever um programa em linguagem BASIC que faça a mesma
simulação em um computador;
5. Preparar um modelo de produção-consumo, fazer uma
simulação e comparar com a simulação do computador.
08. SIMULANDO MODELOS QUANTITATIVOS
Neste capítulo se introduzirá nas técnicas para simulação
de modelos quantitativos de sistemas.
Primeiro, usando um modelo de tanque com água, serão
feitos os cálculos a mão e será preparada uma tabela e
um gráfico com as mudanças em quantidade de água
armazenada em um tanque versus o tempo.
Aqueles que dispõem de um computador poderão fazer
os cálculos nele ("rodar o programa“).
8.1 INTRODUÇÃO.
Os diagramas de energia são uma maneira de visualizar o
comportamento os sistemas. Seis modelos de sistemas com
diferentes tipos de fontes de energia e depósito foram
introduzidos nos capítulos 6 e 7.
Logo após a leitura do capítulo e seus exercícios, estará
apto para fazer gráficos de como cada ecossistema
responde em função do tempo.
A linguagem simbólica de energia se tornará muito mais
compreensível, sobretudo quando se usar números reais
para mostrar o comportamento do sistema.
8.1 INTRODUÇÃO.
A linguagem de diagramas de energia, que temos estado
utilizando, são expressões matemáticas. Com ele temos
representado os processos e suas relações.
Eles se converterão em expressões matemáticas quando
colocaremos números nos símbolos de processo, depósito
e nas razões de fluxo, em cada trajetória de energia.
Se fizermos isto, teremos uma linguagem próxima a
aquela que os computadores conseguem compreender.
8.1 INTRODUÇÃO.
Se tiver acesso a um computador, poderá fazer exercícios
de simulação desta seção usando a máquina.
Se você tem um acesso restringido, pode obter uma
demonstração de simulação de computador usando os
programas do Capítulo 8.
Se não tem acesso a um computador, não se preocupe,
poderá fazer os cálculos manualmente.
8.2 Coeficientes para parâmetros simples.
Para representar quantitativamente o que está sucedendo em um
modelo, a qualquer hora, se escreve números sobre diagramas. As
razões de fluxo se escrevem sobre parâmetros de linhas e
quantidades em depósito se escrevem nos símbolos de depósito.
Imagine um tanque contendo 20 litros (20 l) de água. Uma
mangueira drena 10% do fluxo de água restante por hora. Durante as
primeiras horas, o tanque drenará 2 l de água. Os diagramas para
quantidade de energia destes sistemas são assim:
Figura 8.1 Modelo de tanque de drenagem.
Vazão de saída proporcional ao estoque
8.2 Coeficientes para parâmetros simples.
Este diagrama descreve quantitativamente o sistema
durante a primeira hora.
Não obstante, ao início da segunda hora, as coisas
mudam; a quantidade de água restante no tanque é
agora 18 l e a razão de fluxo de saída é 10% deste, ou
seja 1.8 l/h.
Como a cada hora que passa, mais água corre, e os
valores diminuem.
8.2 Coeficientes para parâmetros simples.
Devemos encontrar uma equação para estes cálculos:
•
Primeiro, representamos a quantidade de água em depósito com
Q (e admitimos que Q mudará com o tempo).
•
Segundo, descrevemos o fluxo como um coeficiente de
parâmetro, chamado k, que indica a fração de água restante
drenada por hora. Quanto maior o diâmetro da mangueira,
maior é o coeficiente. Neste exemplo o coeficiente de
parâmetro é 0.1 (ou 10%).
8.2 Coeficientes para parâmetros simples.
Note que:
Fluxo = coeficiente de fluxo x quantidade no depósito
Fluxo= k x Q,
onde
Coeficiente de fluxo = fração decimal do fluxo de depósito
por unidade de tempo.
O modelo para o tanque com drenagem seria assim:
Figura 8.2 Modelo de tanque mostrando Q e kQ,
sendo k = 0.1.
8.3 SIMULAÇÃO MANUAL.
Tabela 8.1 Cálculos do fluxo do depósito por hora.
Tempo
(horas)
Quantidade Q
(litros)
Taxa de saída
(litros/hora)
0
20
2
1
18
1.8
2
16.20
1.62
3
14.58
1.46
4
13.12
1.31
5
11.81
1.18
8.3 SIMULAÇÃO MANUAL.
Vamos examinar os cálculos que foram feitos para
formar a Tabela 8.1. O processo se calculou repetindo
uma série de subtrações, passo a passo, do depósito.
O processo de fazer cálculos repetidos como este se
chama interação. O processo nas Figuras 8.1 e 8.2
poderia ser exposto como segue:
A quantidade, a um intervalo de tempo próximo, é a
quantidade no momento presente, menos o fluxo:
(novo Q) = (anterior Q) - (k x Q)
8.3 SIMULAÇÃO MANUAL.
Em outras palavras, a equação diz:
“O novo Q igualará ao anterior Q menos k vezes o
anterior Q.”
Vamos escrever esta equação da forma como
aparece na tela de um computador, onde:
* significa multiplicação,
= significa "será igual".
Q=Q-k*Q
8.3 SIMULAÇÃO MANUAL.
No modelo do tanque com drenagem (Figura 8.1) o valor inicial
para Q é 20 lt e k = 0.1 por hora. Na tabela 8.1 o início é indicado
pela primeira linha, quando o tempo é 0, e Q é 20 l.
1.Para a primeira hora:
Fluxo de saída = K*Q
= 0.1/h * 20 lt
= 2 l/h (2 litros por hora)
Ao final da primeira hora a água no depósito se calcula pela
subtração do fluxo por hora:
Novo depósito = depósito anterior - fluxo
= 20 l - 2 l = 18 l
8.3 SIMULAÇÃO MANUAL.
2.Para a segunda hora:
Fluxo de saída = k*Q
= 0.1 l/h * 18 l
= 1.8 l/h
Ao final da segunda hora:
Novo depósito = depósito anterior- fluxo
= 18 l - 1.8 l
= 16.2 lt
A Tabela 8.1 mostra os cálculos para as primeiras cinco horas. Se
estendemos esta tabela para 20 horas e se plotamos os pontos
sobre o gráfico se obteria a seguinte figura:
Figura 8.3
Gráfico da quantidade em depósito (Q) versus tempo, como foi
calculado na Tabela 8.1, simulando o modelo da Figura 8.2.
8.4 DIAGRAMA DE FLUXO.
Quando uma série de cálculos são feitos uma e outra vez, podese escrever os passos do procedimento de cálculo, como um
diagrama de fluxo.
A Figura 8.4 é o diagrama de fluxo para os cálculos que são
feitos na Tabela 8.1.
Lendo desde o início até o fim, é necessário fornecer os números
iniciais; transladar os valores sobre um gráfico, calcular os
valores depois do intervalo de tempo, retornar, graficar e
calcular novamente, assim sucessivamente até chegar a 20
horas.
8.4 DIAGRAMA DE FLUXO.
A lista de passos, no procedimento, é chamado
programa.
O gráfico de fluxo (Figura 8.4) é uma maneira de
escrever um programa.
Escreve-se uma lista de instruções para um
computador, conhecido como programa de
computador.
Figura 8.4
Fluxograma para
simular o modelo da
Figura 8.2 e
calcular os valores
da tabela 8.1
8.5 SIMULAÇÃO POR COMPUTADOR.
Para fazer que um computador realize cálculos,
devemos dar uma lista de instruções na
linguagem na qual ele esteja programado para
responder.
As palavras e símbolos que necessitamos usar
para instruir um computador, estão na Tabela
8.2. Eles são parte da linguagem BASIC.
8.5 SIMULAÇÃO POR COMPUTADOR.
As instruções necessárias se dão na Tabela 8.3, estas se
numeram: 10, 20, 30, etc.
Depois de que o programa tenha sido digitado e esteja
armazenado na memória
de trabalho e na tela, digite RUN e o computador
seguirá a lista de instruções até que os cálculos estejam
completos.
Os números calculados a mão na Tabela 8.1 serão
listados na tela. Para que o programa apareça na tela,
digite LIST. Para salvar o programa em um disquete,
digite SAVE e o nome do programa.
8.5 SIMULAÇÃO POR COMPUTADOR.
Tabela 8.2. Algumas instruções em linguagem
Comando
O que faz
run
Roda o programa, listando as instruções em ordem
numérica
go to
Vai à instrução designada pelo número e o executa
no texto
if
Dá uma instrução para realizar alguma coisa, como
ir a outra linha ( ex: IF T é menor que 20, GO TO )
print
Mostra na tela o valor numérico das quantidades
que se listaram depois do comando PRINT
pset
Mostra na tela um ponto relativo aos novos valores
das variáveis
8.5 SIMULAÇÃO POR COMPUTADOR.
Tabela 8.2. Algumas instruções em linguagem
Comando
end
=
+
*
/
>
<
O que faz
Detém o programa
Dispõe uma parte igual para o que é especificada
Adiciona a próxima quantidade
Multiplica a próxima quantidade
Divide entre a próxima quantidade
Menor que
Maior que
8.5 SIMULAÇÃO POR COMPUTADOR.
Se o computador é compatível com IBM-PC,
imprima o programa teclando LLIST.
Para imprimir os cálculos,
tecle CTRL PRTSC e então RUN.
8.5 SIMULAÇÃO POR COMPUTADOR.
O que segue, é uma explicação das instruções no
programa (Tabela 8.3 e Figura 8.2).
Primeiro, dissemos ao computador o tamanho das
quantidades com as que se trabalharam ao princípio.
Assim, teremos (na Tabela 8.3):
10 Q = 20 (quantidade em depósito = 20).
20 k = 0.1 (coeficiente de parâmetro = 0.1).
30 T = 0 (tempo = 0).
8.5 SIMULAÇÃO POR COMPUTADOR.
Logo diremos ao computador que imprima estes
números:
40 PRINT T, Q, k*Q
A seguir, diremos ao computador o que fazer com estes
números:
50 Q = Q - k*Q
o qual significa, "novo Q é igual ao anterior Q menos k
multiplicado pelo anterior Q".
(Note: * significa multiplicar, para evitar confusão
sobre o significado de x).
8.5 SIMULAÇÃO POR COMPUTADOR.
Havendo feito isto, diremos ao computador para
avançar no tempo uma unidade:
60 T = T + 1
e então, se T é menor de que 20, repetirá as instruções
40, 50 e 60:
70 IF T < 20 GO TO 40
O computador repete os cálculos para cada novo
intervalo de tempo, imprime os resultados e avança o
tempo até chegar a T = 20. Neste ponto, quando chega
à instrução 70, não volta a 40, em lugar disto, vai à
linha 80 a qual diz:
80 END
8.5 SIMULAÇÃO POR COMPUTADOR.
A seqüência inteira dos cálculos leva poucos segundos;
os resultados são listados na tela em forma de tabela.
Agora pode simular a mão.
Tabela 8.3. Programa em BASIC para uma simulação de modelo
na Figura 8.1.
10 Q = 20
20 K = .1
30 T = 20
40 PRINT T, Q, K*Q
50 Q = Q – K*Q
60 T = T + 1
70 IF T < 20 GO TO 40
80 END
8.5 SIMULAÇÃO POR COMPUTADOR.
Para pôr o computador em modo gráfico e mudar a cor
do fundo para branco, é necessário outra instrução:
Em IBM PC as instruções são:
SCREEN 1,0: COLOR 0,0
e
PSET (T/0.07, 180-Q/200) , 3
A curva na Figura 8.3 mostra que a taxa de fluxo
diminui proporcionalmente à diminuição da pressão de
água no tanque. Um programa similar está no Apêndice
Tabela A.8.
8.6 UMA SIMULAÇÃO MAIS COMPLEXA.
Se puder trabalhar na simulação do tanque com drenagem sem
dificuldade, está pronto para um modelo mais complexo.
Figura 8.5 – Diagrama de ecossistema
8.6 UMA SIMULAÇÃO MAIS COMPLEXA.
Este modelo representa qualquer ecossistema.
Mostra que, a luz solar que chega aos produtores, é
capturada durante a fotossíntese e se armazena
como biomassa até ser consumida por tecidos
animais ou tecidos vegetais durante a noite.
Agora vamos adicionar alguns dados.
8.6 UMA SIMULAÇÃO MAIS COMPLEXA.
A incidência da luz do sol varia durante o ano, pode prover
as seguintes quantidades de energia:
Ano
Estação
Luz Solar
(E3 J/m²/estação)
1
Inverno
5.000
1
Primavera
10.000
1
Verão
15.000
1
Outono
10.000
8.6 UMA SIMULAÇÃO MAIS COMPLEXA.
Se as plantas capturam e armazenam 1% da energia solar disponível
(k1 = 0.001), e se os animais consomem 20% da energia total
armazenada nos tecidos das plantas (k2 = 0.2), então o modelo
quantitativo se parece à Figura 8.6.
Diagrama de ecossistema com coeficientes de parâmetros.
8.6 UMA SIMULAÇÃO MAIS COMPLEXA
A simulação manual do sistema poderia
começar com Q=0.1.
Faça os cálculos em cada linha de esquerda a
direita, para ver se pode reproduzir os números
da Tabela 8.4.
Trabalhe com duas casas decimais.
8.6 UMA SIMULAÇÃO MAIS COMPLEXA
Tabela 8.4. Cálculos para simulação manual do diagrama na
Figura 8.6. As unidades são E3 joules/m2/estação.
Vegetal
P=
0.001*S
Quantidade de
Consumo Animal matéria vegetal
C=0.2*(anterior
Novo Q= anterior
Q)
Q+P-C
Ano
Estação
Luz
Solar
S
0
Começo
-
-
-
0.1
1
Inverno
Primavera
Verão
Outono
5.000
10.000
15.000
.
5
10
15
.
0.2*0.1=0.02
0.2*5.1=1
0.214=2.8
.
0.1+5-0.02=5.1
5+10-1=14
14+15-2.8=26
.
Inverno
Primavera
Etc...
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
2
8.6 UMA SIMULAÇÃO MAIS COMPLEXA.
Continue os cálculos em outra folha de papel até
obter dados de cinco anos.
Grafique os valores para luz solar e quantidade de
matéria vegetal sobre a Figura 8.7, continuando o
gráfico iniciado acima.
8.6 UMA SIMULAÇÃO MAIS COMPLEXA.
Figura 8.7 Gráfico da simulação do modelo do sistema P-R versus tempo.
8.7 COMPUTADOR PARA SIMULAÇÃO DE
PRODUÇÃO E MODELO DE CONSUMO
O modelo na Figura 8.5, que foi "simulado
manualmente" na Figura 8.7 e que pode ser
escrito para simulações em computador como
se mostra na Tabela 8.5.
Os programas para este e outros modelos estão
listados no Apêndice A.
8.7 SIMULAÇÃO DE PRODUÇÃO E CONSUMO
Tabela 8.5. Programa de computador em BASIC para o modelo P-R na Figura
8.6. (resultados multiplicados por 1000.)
10 Q = 0.1
90 N = N + 1
20 K1 = 0.001
100 IF N=5 THEN N=1
30 K2 = 0.02
110 PRINT T, S, P, C, Q
40 N = 1
120 P = K1*S
50 IF N=1THEN S=5.000
130 C = K2*Q
60 IF N=2 THEN S=10.000
140 Q = Q + P – C
70 IF N=3 THEN S=15.000
150 T = T + 1
80 IF N=4 THEN S=10.000
200 IF T>20 GO TO 50
8.7 COMPUTADOR PARA SIMULAÇÃO DE
PRODUÇÃO E MODELO DE CONSUMO
Para o IBM PC, as mudanças seriam:
5 SCREEN 1,0: COLOR 0,0
6 LINE (0,0)-(320,180),1,B
110 PSET (T/.07, 50-S/350)
115 PSET (T/.07, 180-Q/.5)
200 IF T/.07<320 GO TO 50
8.7 COMPUTADOR PARA SIMULAÇÃO DE
PRODUÇÃO E MODELO DE CONSUMO
Para imprimir o gráfico na folha de papel, no
IBM PC pressione a tecla SHIFT com a tecla
PRINT SCREEN (PrtSc).
8.8 DISCUSSÃO.
O anterior gráfico de quantidade, mostra o
crescimento e o estado estacionário como no Modelo 3
(Figura 6.3). O sol é uma fonte renovável fixa com
fluxo constante.
Em conseqüência, a produção de material vegetal
aumenta rapidamente ao princípio, mas desde que o
consumo animal é uma porcentagem fixa de material
vegetal disponível, os consumidores começam a
aumentar rapidamente até que a produção e consumo
são iguais.
8.8 DISCUSSÃO.
O crescimento não é uniforme por causa das variações
da luz solar e o pico de crescimento vegetal está depois
do pico de luz solar porque há um atraso na formação do
depósito da energia (armazenamento).
Um exemplo deste tipo de crescimento é a sucessão
ecológica. O crescimento rápido de plantas em um
campo aberto muda a um crescimento líquido mais lento
de arbustos e logo árvores, e culmina em um estado
estacionário onde árvores e outros produtores estão em
balanço com os consumidores.
8.9 Perguntas e atividades para o capítulo 8
1. Defina os seguintes termos:
a. simulação
b. quantitativo
c. coeficiente
d. programa
e. interação
f. equação diferencial
g. BASIC
2. Faça o diagrama de uma descrição quantitativa de um
tanque de depósito mostrando depósito e fluxos.
8.9 Perguntas e atividades para o capítulo 8
3. Calcule os coeficientes de parâmetro de seu tanque
de depósito, se a quantidade armazenada é 100
litros e o primeiro fluxo é de 5 litros/hora.
4. Usando Tabela 8.1 como guia, faça uma lista dos
dados da pergunta #3 em forma de tabela, esta
será estendida para um mínimo de 15 horas.
5. Use seus dados do #4 para graficar quantidade (Q)
versus tempo (T) sobre eixos cartesianos.
8.9 Perguntas e atividades para o capítulo 8
6. Grafique os valores para luz solar e quantidade de
matéria vegetal na Figura 8.7, continuando o
gráfico já começado.
7. Explique porque o gráfico que completou na Figura
8.7 é semelhante a um gráfico de sucessão
ecológica.
8. Como provavelmente você sabe, há muito mais
sobre programação de computadores em
linguagem BASIC, mas você conhece o suficiente
para simular os Modelos nos Capítulos 6 e 7.
Tente então.
Download

Cap. 08 - Unicamp