Programacion no resultante
Publicado por Ana (1 intervención) el 11/05/2017 17:02:56
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
Dim CA, CB, deltaP, deltaCA, deltaCB As Double
Dim Fila As Integer
Sub Rendimientodereacciónsinrecirculación()
'Datos de entrada
'Parámetros del problema
Fe = Worksheets("Hoja1").Cells(5, 3)
Fs = Worksheets("Hoja1").Cells(6, 3)
V = Worksheets("Hoja1").Cells(7, 3)
R = Worksheets("Hoja1").Cells(8, 3)
'Condiciones iniciales
Cao = Worksheets("Hoja1").Cells(11, 3)
Cbo = Worksheets("Hoja1").Cells(12, 3)
'Parámetros de integración
tmax = Worksheets("Hoja1").Cells(15, 3)
deltat = Worksheets("Hoja1").Cells(16, 3)
'Inicialización de variables
CA = CaoCB = CboFila = 4deltaP = tmax / 1000
'Integración de las ecuaciones de balance de masa
For t = 0 To tmax Step deltatMeA = (Fe * Cao) / V
MsA = (Fs * CA) / V
McA = R * CA
deltaCA = deltat(MeA - MsA - McA)
CA = CA + deltaCA
MsB = (Fs * CB) / V
deltaCB = deltat(-MsB + McA)
CB = CB + deltaCB
'Lazo para imprimir Resultados
If t = 0 ThenWorksheets("Hoja2").Cells(Fila, 2) = 0
Worksheets("Hoja2").Cells(Fila, 3) = Cao
Worksheets("Hoja2").Cells(Fila, 4) = Cbo
deltaP = tmax / 1000 + deltaP
Fila = Fila + 1
ElseIf t >= deltaP Then
Worksheets("Hoja2").Cells(Fila, 2) = t
Worksheets("Hoja2").Cells(Fila, 3) = CA
Worksheets("Hoja2").Cells(Fila, 4) = CB
deltaP = tmax / 1000 + deltaP
Fila = Fila + 1
End If
Next t
End Sub
Valora esta pregunta
0