1. Realizar una representación esquemática del problema 2




descargar 1.52 Mb.
título1. Realizar una representación esquemática del problema 2
página7/8
fecha de publicación23.11.2015
tamaño1.52 Mb.
tipoDocumentos
med.se-todo.com > Química > Documentos
1   2   3   4   5   6   7   8



Coeficiente global de transferencia de calor

Userviciocalentamiento = 0.3505 kW/m2 K Uservicioenfriamiento = 0.2629 kW/m2 K Uproceso-proceso = 0.2629 kW/m2 K




Costo de los cambiadores de calor:

Cp = 1456.3 A0.6

($, m 2 )

Tiempo de operación: 8500 hr/año Rapidez de retorno = 0.1

IIII. [Secuencias de separación] Calcule la secuencia de separación con el costo mínimo, la segunda mejor, la tercera mejor y la peor



Separador único

Costo

[$/año]

A/B

26,100

B/C

94,900

C/D

59,300

A/BC

39,500

AB/C

119,800

B/CD

112,600

BC/D

76,800

A/BCD

47,100

AB/CD

140,500

ABC/D

94,500



la alimentación al sistema de separacion es la siguiente:

propano:100, iso-butano: 300, n-butano:500, iso-pentano: 400, todos estos valores están expresados en lbmol/h

6. Códigos (GAMS) para la resolución de los ejemplos.



Solución del ejemplo1
Modelo

OPTION LIMROW=0; OPTION LIMCOL=0;
* declaracion de las variables del problema


A0

max G = 50CBV

[10C
V + 50CA
V + 30(e0.5CC

-1)]
positive variables Ca, Cb, Cc, t, ra, rb, rc;

* declaracion de la varibale de la funcion objetivo free variable G;

* declaracion del conjunto de ecuaciones del modelo

sujeta a las restricciones

equations obj, EcCa, EcCb, EcCc;


C A =
C =

C A0e k1t

k1C A0 [e k1t
e k2t ]

* declaracion de las expresiones de las ecuaciones
* expresion de la funcion objetivo G [=] $ / lote

obj.. G =e= 500*50*Cb - (500*10*1.6 + 500*50*Ca +

30*(exp(0.5*Cc)-1));

B k2

k1




k e k1t k


e k2t
* expresiones para las concentraciones

EcCa.. Ca =e= 1.6*exp(-0.4*t);

EcCb.. Cb =e= 0.4*1.6*(exp(-0.4*t)-exp(-0.01*t))/(0.01-0.4); EcCc.. Cc =e= 1.6*(1-(0.01*exp(-0.4*t)-0.4*exp(-0.01*t))/(0.01-

CC =

C A0 1 2 1

0.4));

k2 k1

el código en GAMS para resolver esto es
* inicio del archivo abcprofit.gms

* modelo para la maximizacion de la ganancia

* para un sistema de reaccion consecutiva

* A -> B -> C en un reactor por lotes

$ TITLE abcprofit

$ OFFSYMXREF

$ OFFSYMLIST
* asignacion del nombre al problema model abcprofit /all/;

* declaracion de solucion del problema solve abcprofit using NLP maximizing G;

*fin de la codificacion del archivo abcprofit gms




Con solución



















g1.. -2*x1 - x2 =g= -5;

g2 .. -x1 - x3 =g= -2;




LOWER

LEVEL

UPPER










g3.. -x1 - 2*x2 - x3 =g= -10;

VAR Ca

.

0.018

+INF










h1.. 2*x1 - 2*x2 + x3 =e= -2;

VAR Cb

.

1.449

+INF










h2.. 10*x1 + 8*x2 - 14*x3 =e= 26;

VAR Cc

.

0.133

+INF










h3.. -4*x1 + 5*x2 - 6*x3 =e= 6;

VAR t

horas

.

11.202

+INF

[

=

]


* declaracion de los limites de las variables


VAR G -INF 27758.374 +INF [ = ] $/carga
a. Ejemplo de codificación para la solucion del problema 8.26 del libro de Himmelblau (2001)
***********************************************************

******************

$ TITLE problema 8 26

$ OFFSYMXREF

$ OFFSYMLIST
OPTION LIMROW=0; OPTION LIMCOL=0;
* ejemplo para resolver el problema 8 26 del libro

* declaracion de las variables del problema positive variables x1, x2, x3 ;
* declaracion de la varibale de la funcion objetivo free variable f;
* delcaracion del conjunto de ecuaciones del modelo equations obj, g1, g2, g3, h1, h2, h3;
* declaracion de las expresiones de las ecuaciones obj.. f =e= x1**2 + x2 **2 + x3 **2;
x1.l = 1.; x2.l = 2.; x3.l = 0.;
* asignacion del nombre al problema model p826 /all/;

* declaracion de solucion del problema
solve p826 using NLP minimizing f;
* fin de la codificacion del archivo p826 gms


LOWER LEVEL UPPER MARGINAL
---- VAR x1 . 1.000 +INF .


---- VAR x3

.

.

+INF

.




---- VAR f

-INF

5.000

+INF




.



---- VAR x2 . 2.000 +INF .

Solución ejemplo 2 por medio de Gams
Actividades
1. Representación esquemática

2. Definición de variables

3. Formulación de la función objetivo

4. Formulación de la serie de restricciones

5. Desarrollo o selección y aplicación del algoritmo de cálculo

6. Evaluación de resultados

¿ Qué proceso se construye (II o III)? El II

¿ Cómo se obtiene el químico B? Se produce a través del proceso

I

¿ Qué cantidad deberá producirse del químico C? 10 ton/h


% inicio del archivo ehlp8_26.m

% codigo para calular el minimo de una funcion

% con restricciones correspondiente al problema

% 8.26 del libro de Edgar, Himmelblau & Lasdon

% min f(x)=x1^2+x2^2+x3^2

% sujeta a las restricciones siguientes

% g1(x)=2x1+x2<=5

% g2(x)=X1+X3<=2

% g3(X)=x1+2x2+x3<=10

% h1(x)=2x1-2x2+x3=-2

% h2(x)=10x1+8x2-14x3=26

% h3(x)=-4x1+5x2-6x3=6

% x1>=1 x2>=2 x3>=0

% con un estimado inicial de

% x0=[1 1 1]

clc; clear all; format compact;
x0=[1 1 1]; % estimado inicial

A=[2 1 0; 1 0 1; 1 2 1]; % coeficientes de las eciuaciones de

% restricciones de desigualdad

B=[5;2;10]; % termino independiente

Aeq=[2 -2 1; 10 8 -14; -4 5 -6]; % coeficientes de las restricciones de igualdad

Beq=[-2; 26; 6]; % termino independiente

Lb=[1 2 0];

x0=[-1 -1 -1 ]; % estimado inicial options=optimset('LargeScale','off','Display','iter'); [x,fmin,flag,output,lambda,dfdxi,d2fdxi2]=fmincon('p8_26',x0,A, B,Aeq,Beq,Lb,[],[],options)

x1=x(1); x2=x(2); x3=x(3);

lambda_mu=eig(d2fdxi2); fprintf('x1 = %10.4f\n',x1) fprintf('x2 = %10.4f\n',x2) fprintf('x3 = %10.4f\n',x3) fprintf('fmin = %10.4f\n',fmin)

fprintf('Valores propios de la Hessiana = %f\n',lambda_mu)

% fin del archivo ehlp8_26.m

% inicio del archivo p8_26.m function f=p8_26(x)

x1=x(1); x2=x(2); x3=x(3);

f = x1^2 + x2^2 + x3^2;

% fin del archivo p8_26.m

max Directional

First-order

Iter F-count f(x) constraint Step-size derivative optimality Procedure


1

9

5

3.553e-015

1

-0.0006

Valores propios de la Hessiana = 1.999586

23.2



















2

14

5

3.553e-015

1

4.21e-015

Template GAMS [ regresar ]

23.2

Hessian




















Active Constraints: 1 2 3 x = 1.0000 2.0000 0.0000

fmin = 5.0000

flag = 1

* inicio del archivo template.gms

m * archivo base para resolver problemas de optimizacion o

di * definicion del titulo del problema

fi $ TITLE nombre del problema e $ OFFSYMXREF

d $ OFFSYMLIST
* se recomienda no quitar las lineas siguientes

OPTION LIMROW=0; OPTION LIMCOL=0;

* solver

output = iterations: 2 funcCount: 14 stepsize: 1 algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'

firstorderopt: 1.7764e-015 cgiterations: [ ]

lambda = lower: [3x1 double] upper: [3x1 double]

eqlin: [3x1 double] eqnonlin: [0x1 double]

ineqlin: [3x1 double] ineqnonlin: [0x1 double]

dfdxi =

2.0000

4.0000

0 d2fdxi2 =

1.3324 0.3330 0.3330

0.3330 1.3336 0.3336

0.3330 0.3336 1.3336
x1 = 1.0000 x2 = 2.0000 x3 = 0.0000
fmin = 5.0000
Valores propios de la Hessiana = 1.000000

option nlp = conopt; option mip = cplex; option rminlp = conopt2; option minlp = dicopt;

*------------------------------------------------------

* declaracion del conjunto de variables del problema

* en estas lineas se declaran las variables positivas

*------------------------------------------------------

* free variable

* positive variable

* negative variable

* binary variable
*------------------------------------------------------

* declaracion del nombre de las funciones

*------------------------------------------------------

* equations
*------------------------------------------------------

* declaracion de las expresiones de las ecuaciones

*------------------------------------------------------

* ejemplo se tiene f = 2 x2 + x1



*------------------------------------------------------

* asignacion del nombre al problema

*------------------------------------------------------

*model nombre_problema /all/;

*------------------------------------------------------

* declaracion de solucion del problema

*------------------------------------------------------

* solve nombre_problema using mip o lp o nlp minimizing o maximaxing variable;
*fin del archivo template.gms

Código en Gams para la solución del problema de optimización de un reactor por lotes


* inicio del archivo optimbatch.gms

* archivo base para resolver problemas de optimizacion
* definicion del titulo del problema

*=============================================

$ TITLE optimizacion de un reactor batch

*=============================================

$ OFFSYMXREF

$ OFFSYMLIST
$ontext Optimizacion de un reactor por lotes

Se lleva a cabo una reaccion en serie A -> B -> C en fase liquida en un

reactor por lotes de 500 dm3. La concentracion inicial de A

es de 1.6 mol/dm3.

El producto deseado es B y la separacion del producto no deseado C es muy

dificil y costosa. Dado que la reaccion se efectua a una

relativamente alta, la reaccion se puede extinguir facilmente.
Informacion adicional:

Costo del reactivo A puro: $10 / mol A Precio de venta de B puro: $50 / mol B

Costo de separar A de B: $50 / mol de A Costo de separar C de B: $30 [exp(0.5Cc)-1] k1 = 0.4 h^-1

k2 = 0.01 h^-1 a 100 C

$offtext


* se recomienda no quitar las lineas siguientes

OPTION LIMROW=0; OPTION LIMCOL=0;

* solver

option lp = minos; option nlp = conopt; option mip = cplex; option rminlp = conopt2; option minlp = dicopt;


$set Np 700

*============================================= Set Np Numero de subintervalos / 0*%Np% /;

Alias (Np,i);
*=============================================

* parametros del sistema

Scalar Cmp costo quimico A /10/

Pvb precio de venta quimico B /50/

Csab costo de separacion del quimico A de B /50/

k1 coeficiente cinetico 1 h-1 /0.4/

k2 coeficiente cinetico 2 h-1 /0.01/ V volumen del reactor dm3/500/ Ca0 mol por dm3 /1.6/

Cc0 mol por dm3 /0/;

*------------------------------------------------------

* declaracion del conjunto de variables del problema

* en estas lineas se declaran las variables positivas

*------------------------------------------------------ free variable G ganancia usd por lote;

*=============================================

positive variable Ca(Np) mol A por dm3

Cb(Np) mol B por dm3

Cc(Np) mol C por dm3 t(Np) horas

dt horas tf horas;

*=============================================

* negative variable

* binary variable

*------------------------------------------------------

* declaracion del nombre de las funciones

*------------------------------------------------------

equations obj, bma(Np), bmb(Np), bmc(Np), dt_time, timevec;

* ejemplo se tiene f = 2 x2 + x1

* nombre.. f =e= 2*x2 + x1;

*=============================================

* declaracion de las expresiones de las ecuaciones

*=============================================

obj.. G =e= Pvb*Cb['%Np%']*V - ( Cmp*Ca0*V + Csab*Ca['%Np%']*V +

Cb.fx['0']=Cb0; Cc.fx['0']=Cc0;

*------------------------------------------------------

* asignacion del nombre al problema

*------------------------------------------------------ model reactorabc /all/;

*------------------------------------------------------

* declaracion de solucion del problema

*------------------------------------------------------ solve reactorabc using nlp maximaxing G;

*=============================================

* escritura de la serie de resultados

*=============================================

* definicion del nombre del archivo file res /reactor.dat/;

put res;

* escritura de la serie de valores por medio del uso de un ciclo

loop (i, put @1 t.l(i) :< 7:4

@10 Ca.l(i) :<12:8

@22 Cb.l(i) :<12:8

@34 Cc.l(i) :<12:8/

);

putclose res;

* fin del archivo optimbatch.gms
Código para determinar las secuencias de separacion

30*( exp(Cc['%Np%']) - 1) );

bma(Np(i+1)).. Ca[i+1] =e= Ca[i]- k1*Ca[i]*dt; bmb(Np(i+1)).. Cb[i+1] =e= Cb[i]+ (k1*Ca[i]-k2*Cb[i])*dt; bmc(Np(i+1)).. Cc[i+1] =e= Cc[i]+ k2*Cb[i]*dt;

dt_time.. (%Np%-1)*dt =e= tf;

timevec(Np(i+1)).. t[i+1] =e= t[i] + dt;

*=============================================

* especificacion de las condiciones iniciales

*=============================================

* inicio del archivo secuencia_separacion.gms

* archivo base para resolver problemas de optimizacion

* definicion del titulo del problema

$ TITLE secuencia de separacion

$ OFFSYMXREF

$ OFFSYMLIST

* se recomienda no quitar las lineas siguientes

OPTION LIMROW=0; OPTION LIMCOL=0;

* solver

option mip = cplex; option rminlp = conopt2; option minlp = dicopt;

*------------------------------------------------------

* declaracion del conjunto de variables del problema

* en estas lineas se declaran las variables positivas

*------------------------------------------------------ scalar U kgmol por hora/1000/

Ft kgmol por hora/1000/;

*

free variable Ct usd por anio;

positive variable F1,F2,F3,F4,F5,F6,F7,F8,F9,F10 kgmol por hora

Q1,Q2,Q3,Q4,Q5,Q6,Q7,Q8,Q9,Q10, Qt kJ por

kgmol por hora;

* negative variable

binary variable Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10;

*------------------------------------------------------

* declaracion del nombre de las funciones

*------------------------------------------------------ equations OBJ,

BM1,BM2,BM3,BM4,BM5,BM6,
BE1,BE2,BE3,BE4,BE5,BE6,BE7,BE8,BE9,BE10,BE11, DJ1,DJ2,DJ3,DJ4,DJ5,DJ6,DJ7,DJ8,DJ9,DJ10, RL1,RL2,RL3,

*RL4,RL5, RL6,RL7,RL8,RL9,RL10,RL11,RL12,

RL13 nueva restriccion para eliminar la mejor

secuencia

RL14 nueva restriccion para eliminar la segunda mejor secuencia

RL15 nueva restriccion para eliminar la tercera mejor

secuencia

RL11A;
*------------------------------------------------------

*------------------------------------------------------

* ejemplo se tiene f = 2 x2 + x1

* nombre.. f =e= 2*x2 + x1;

OBJ.. Ct =E=

145*Y1+52*Y2+76*Y3+38*Y4+66*Y5+125*Y6+44*Y7+58*Y

8+37*Y9+112*Y10 +
0.42*F1+0.12*F2+0.25*F3+0.14*F4+0.21*F5+0.78*F6+0.11* F7+0.19*F8+0.08*F9+0.39*F10 +

35.3*Qt;

BM1.. F4 + F5 =E= 0.85*F1; BM2.. F6 + F7 =E= 0.8*F3;

BM3.. F10 =E= 0.45*F2 + 0.563*F7;

BM4.. F9 =E= 0.765*F5 + 0.812*F6; BM5.. F8 =E= 0.55*F2 + 0.647*F4;

BM6.. F1 +F2 +F3 =E= Ft;

* balances de energia BE1.. Q1 =E= 0.028*F1; BE2.. Q2 =E= 0.042*F2; BE3.. Q3 =E= 0.054*F3; BE4.. Q4 =E= 0.040*F4; BE5.. Q5 =E= 0.047*F5; BE6.. Q6 =E= 0.024*F6; BE7.. Q7 =E= 0.039*F7; BE8.. Q8 =E= 0.044*F8; BE9.. Q9 =E= 0.036*F9; BE10.. Q10 =E= 0.022*F10;

BE11.. Qt =E= Q1+Q2+Q3+Q4+Q5+Q6+Q7+Q8+Q9+Q10;

* disjuncion M grande


DJ5..

F5

=L= U*Y5;

DJ6..

DJ7..

F6

F7

=L= U*Y6;

=L= U*Y7;

DJ8..

F8

=L= U*Y8;



DJ1.. F1 =L= U*Y1; DJ2.. F2 =L= U*Y2; DJ3.. F3 =L= U*Y3; DJ4.. F4 =L= U*Y4;

DJ10.. F10 =L= U*Y10;

* restricciones logicas

RL1.. Y1 + Y2 + Y3 =E=1; RL2.. Y4 + Y5 =E= Y1;

RL3.. Y6 + Y7 =E= Y3;

*RL4.. Y2 =G= Y3;

*RL5.. Y2 =G= Y10; RL6.. Y2 + Y4 =E= Y8; RL7.. Y5 + Y6 =E= Y9; RL8.. Y2 + Y7 =E= Y10;

RL9.. Y4 + Y5 + Y6 + Y7 =L= 1; RL10.. Y4 + Y9 =L= 1;

RL11.. Y9 + Y10 =L= 1;

RL11A.. Y8 + Y9 + Y10 =L= 2;

RL12.. Y1+Y2+Y3+Y4+Y5+Y6+Y7+Y8+Y9+Y10 =E= 3;

* RESTRICCION PARA ELIMINAR LA MEJOR SOLUCION

RL13.. Y2 + Y8 +Y10 =L= 2;

* RESTRICCION PARA ELIMINAR LA segunda MEJOR SOLUCION

RL14.. Y1 + Y4 +Y8 =L= 2;

* RESTRICCION PARA ELIMINAR LA tercera MEJOR SOLUCION

RL15.. Y1 + Y5 +Y9 =L= 2;


*------------------------------------------------------

* asignacion del nombre al problema

*------------------------------------------------------ model mejor /OBJ,

BM1,BM2,BM3,BM4,BM5,BM6,
BE1,BE2,BE3,BE4,BE5,BE6,BE7,BE8,BE9,BE10,BE11, DJ1,DJ2,DJ3,DJ4,DJ5,DJ6,DJ7,DJ8,DJ9,DJ10, RL1,RL2,RL3,

*RL4,RL5, RL6,RL7,RL8,RL9,RL10,RL11,RL12,RL11A/;

* Ct: 3308.330 Y2= Y8= Y10=1

model peor /OBJ, BM1,BM2,BM3,BM4,BM5,BM6,
BE1,BE2,BE3,BE4,BE5,BE6,BE7,BE8,BE9,BE10,BE11, DJ1,DJ2,DJ3,DJ4,DJ5,DJ6,DJ7,DJ8,DJ9,DJ10, RL1,RL2,RL3,

*RL4,RL5, RL6,RL7,RL8,RL9,RL10,RL11,RL12,RL11A/;

* Ct: 4123.493 Y1= Y5= Y9=1

model segunda_mejor /OBJ, BM1,BM2,BM3,BM4,BM5,BM6,
BE1,BE2,BE3,BE4,BE5,BE6,BE7,BE8,BE9,BE10,BE11, DJ1,DJ2,DJ3,DJ4,DJ5,DJ6,DJ7,DJ8,DJ9,DJ10, RL1,RL2,RL3,

*RL4,RL5, RL6,RL7,RL8,RL9,RL10,RL11,RL12,RL13,RL11A/;

* Ct: 3927.273 Y1= Y4= Y8=1

model tercera_mejor /OBJ, BM1,BM2,BM3,BM4,BM5,BM6,
BE1,BE2,BE3,BE4,BE5,BE6,BE7,BE8,BE9,BE10,BE11, DJ1,DJ2,DJ3,DJ4,DJ5,DJ6,DJ7,DJ8,DJ9,DJ10, RL1,RL2,RL3,

*RL4,RL5,
RL6,RL7,RL8,RL9,RL10,RL11,RL12,RL13,RL14,RL11A/;

* Ct: 4123.493 Y1= Y5= Y9=1 model cuarta_mejor /all/;

* Ct: 4123.493 Y1= Y5= Y9=1
*------------------------------------------------------

* declaracion de solucion del problema

*------------------------------------------------------ solve mejor using mip minimizing Ct;

solve peor using mip maximizing Ct;

solve segunda_mejor using mip minimizing Ct;

solve tercera_mejor using mip minimizing Ct;

solve cuarta_mejor using mip minimizing Ct;

*fin del archivo secuencia_separacion.gms


Código en GAMS para las redes de intercambo de calor
* inicio del archivo serv_min.gms

* calculo del consumo minimo de servicios aplicando

* el metodo de transbordo expandido

$ TITLE calculo de la red de intercambio de calor para costo minimo de servicios

$ OFFSYMXREF

$ OFFSYMLIST OPTION LIMROW=0;

OPTION LIMCOL=0;

*

*option lp = snopt; option lp = xpress; option nlp = conopt; option mip = cplex;

option rminlp = conopt2;

option minlp = dicopt;

* declaracion del conjunto de variables del problema positive variables Qs, Qw MW,

Qs11, Qs12, Rs1, Rs2 MW,

Q112, Q113, Q114, Q124, Q125, R12, R13, R14 MW, Q213, Q214, Q224, Q225, R23, R24 MW,

Q1w5, Q2w5 MW;

binary variables

Ys1a,Y11a,Y11b,Y12b,Y1wb,Y21b,Y22b,Y2wb;

* declaracion de la variables de la funcion objetivo z[=] $/año free variable z usd por anio;

* declaracion de las ecuaciones del modelo

equations obj1, obj2,

* para el intervalo 1 s1, c11,

* para el intervalo 2 s2, h12, c12,

* para el intervalo 3 h13, h23, c13,

* para el intervalo 4 h14, h24, c14, c24,

* para el intervalo 5

h15, h25, c25,w5,

* para las restricciones logicas

RL1,RL2,RL3,RL4,RL5,RL6,RL7,RL8;

* expresiones del modelo

obj1.. z =e= 80000 * Qs + 20000 * Qw;

obj2.. z =e= Ys1a + Y11a + Y11b + Y12b + Y1wb + Y21b + Y22b + Y2wb;

* balance de energia para el intervalo 1

s1.. Qs =e= Qs11 + Rs1;

c11.. Qs11 =e= 30;

* balance de energia para el intervalo 2 s2.. Rs1 =e= Qs12 + Rs2;

h12.. 60 =e= Q112 + R12;

c12.. Qs12 + Q112 =e= 90;

* punto de pliegue

*R12.l = 0;

*Rs2.l = 0;

* balance de energia para el intervalo 3 h13.. 70 + R12 =e= Q113 + R13;

h23.. 140 =e= Q213 + R23;

c13.. Q113 + Q213 =e= 105;

* balance de energia para el intervalo 4 h14.. R13 + 90 =e= Q114 + Q124 + R14; h24.. R23 + 180 =e= Q214 + Q224 + R24; c14.. Q114 + Q214 =e= 135;

c24.. Q124 + Q224 =e= 117;

* balance de energia para el intervalo 5 h15.. 60 + R14 =e= Q125 + Q1w5;

h25.. 120 + R24 =e= Q225 + Q2w5 ;

c25.. Q125 + Q225 =e= 78; w5.. Q1w5 + Q2w5 =e= Qw; RL1.. Qs11 + Qs12 =L= 60*Ys1a; RL2.. Q112 =L= 60*Y11a;

RL3.. Q113 + Q114 =L= 220*Y11b; RL4.. Q124 + Q125 =L= 195*Y12b; RL5.. Q1w5 =L= 220*Y1wb;

RL6.. Q213 + Q214 =L= 240*Y21b; RL7.. Q224 + Q225 =L= 195*Y22b; RL8.. Q2w5 =L= 225*Y2wb;

* para resolver la red con el consumo minimo de servicios model servmin /obj1,

s1, c11,

s2, h12, c12, h13, h23, c13,

h14, h24, c14, c24, h15, h25, c25, w5/;

* para resolver la red con el numero minimo de unidades

model unitmin/obj2, s1, c11,

s2, h12, c12, h13, h23, c13,

h14, h24, c14, c24,

h15, h25, c25, w5, RL1,RL2,RL3,RL4,RL5,RL6,RL7,RL8/;

* declaracion de solucion del problema solve servmin using lp minimizing z; display z.l,Qs.l,Qw.l,

Qs11.l,Rs1.l, Qs12.l,Q112.l,Rs2.l,R12.l,

Q113.l,Q213.l,R13.l,R23.l, Q114.l,Q124.l,Q214.l,Q224.l,R14.l,R24.l, Q125.l,Q225.l,

Q1w5.l,Q2w5.l;

* impresion de resultados al archivo min_utility.out file soln1 /min_utility.out/;

put soln1;

put "z: "; put z.l; put /;

put "Qs11: "; put Qs11.l; put /; put "Rs1: "; put Rs1.l; put /; put "Qs12: "; put Qs12.l; put /;

put "Rs2: "; put Rs2.l; put /; put "R12: "; put R12.l; put /; put "Q113: "; put Q113.l; put /; put "Q213: "; put Q213.l; put /; put "R13: "; put R13.l; put /; put "R23: "; put R23.l; put /; put "Q114: "; put Q114.l; put /; put "Q124: "; put Q124.l; put /; put "Q224: "; put Q224.l; put /; put "R14: "; put R14.l; put /; put "R24: "; put R24.l; put /; put "Q125: "; put Q125.l; put /; put "Q225: "; put Q225.l; put /;

put "Q1w5: "; put Q1w5.l; put /;

put "Q2w5: "; put Q2w5.l; put /;

*

solve unitmin using mip minimizing z;

display z.l,Qs.l,Qw.l, Qs11.l,Rs1.l,

Qs12.l,Q112.l,Rs2.l,R12.l, Q113.l,Q213.l,R13.l,R23.l,

Q114.l,Q124.l,Q214.l,Q224.l,R14.l,R24.l, Q125.l,Q225.l,

Q1w5.l,Q2w5.l,
Ys1a.l,Y11a.l,Y11b.l,Y12b.l,Y1wb.l,Y21b.l,Y22b.l,Y2wb.l;

* impresion de resultados al archivo units_min.out file soln2 /units_min.out/;

put soln2;

put "z: "; put z.l; put /;

put "Qs11: "; put Qs11.l; put /;

put "Rs1: "; put Rs1.l; put /; put "Qs12: "; put Qs12.l; put /; put "Q112: "; put Q112.l; put /; put "Rs2: "; put Rs2.l; put /;

put "R12: "; put R12.l; put /;

put "Q113: "; put Q113.l; put /;

put "R13: "; put R13.l; put /; put "R23: "; put R23.l; put /; put "Q114: "; put Q114.l; put /; put "Q124: "; put Q124.l; put /; put "Q224: "; put Q224.l; put /; put "R14: "; put R14.l; put /; put "R24: "; put R24.l; put /; put "Q125: "; put Q125.l; put /; put "Q225: "; put Q225.l; put /;

put "Q1w5: "; put Q1w5.l; put /; put "Q2w5: "; put Q2w5.l; put /; put "Ys1a: "; put Ys1a.l; put /; put "Y11a: "; put Y11a.l; put /; put "Y11b: "; put Y11b.l; put /; put "Y12b: "; put Y12b.l; put /; put "Y1wb: "; put Y1wb.l; put /; put "Y21b: "; put Y21b.l; put /; put "Y22b: "; put Y22b.l; put /; put "Y2wb: "; put Y2wb.l; put /;

*fin del archivo serv_min.gms
1   2   3   4   5   6   7   8

similar:

1. Realizar una representación esquemática del problema 2 iconResumen -se presenta una solución sencilla a un problema de investigación...

1. Realizar una representación esquemática del problema 2 iconII. hacia una representación gráfica del amor

1. Realizar una representación esquemática del problema 2 iconLos fragmentos reunidos en este libro constituyen una he-rramienta...
«lenguaje», está contenida una deter­minada concep­ción del mun-do*23, se pasa al segundo momento, el de la crítica y el conoci­miento,...

1. Realizar una representación esquemática del problema 2 iconEl problema del conocimiento esta vinculado al problema de la verdad

1. Realizar una representación esquemática del problema 2 iconEl problema del conocimiento esta vinculado al problema de la verdad

1. Realizar una representación esquemática del problema 2 iconResumen : se pretende aportar una visión del ciclo completo del problema...

1. Realizar una representación esquemática del problema 2 icon1 El teatro latino. Los géneros: tragedia y comedia. Principales...

1. Realizar una representación esquemática del problema 2 iconEl profesional del área se encuentra con el problema de lo que influye...

1. Realizar una representación esquemática del problema 2 iconResumen Existen diversas propuestas para modelar moléculas orgánicas,...

1. Realizar una representación esquemática del problema 2 iconNota de prensa /10/2013 Una nueva solución para el problema del lindano


Medicina



Todos los derechos reservados. Copyright © 2015
contactos
med.se-todo.com