GAMSExampleCode

XXI Seminario de Ingeniería Química

2003

Instituto Tecnológico de Celaya

Código GAMS
El código de GAMS se puede escribir con cualquier procesador de texto o a
través de la interfase de GAMS. Si se utilizan procesadores especializados como
Word, FrameMaker, PageMaker, etc., asegúrese de guardar el archivo sin
formato (como texto, código ASCII).
Los archivos de GAMS deberán tener la extensión*.gms
Luego de la solución de algún modelo, GAMS crea un archivo de resultados
también en formato de texto y con el mismo nombre que el archivo del código,
pero con extensión *.lst

Como regla general, un modelo de GAMS debe contener las siguientes partes
(se muestra un caso ilustrativo):
1) Título
$TITLE MIMODELO
2) Declaración de Conjuntos
SETS
J COMPONENTES /1*3/
3) Declaración de ParámetrosPARAMETERS ALFA(J), Z(J), Q;
4) Declaración de Variables (positivas y generales)
POSITIVE VARIABLES TETA(I), XD(J), RMIN;
VARIABLES C;
5) Declaración de Ecuaciones
EQUATIONS EQ1(I), EQ2(I), EQ3, OBJ;
6) Ecuaciones del Sistema
EQ1(I).. SUM(J,((ALFA(J)*Z(J))/(ALFA(J)-TETA(I))))=E= 1.0 - Q;

Optimización de Procesos

36

XXI Seminario de Ingeniería Química

2003

Instituto Tecnológico de Celaya

Noteque el identificador de la ecuación va precedido de dos puntos. En las
ecuaciones el símbolo =E= significa igual, =G= significa mayor que y =L=
significa menor que.
7) Definición de una función objetivo (“Dummy” o verdadera)
OBJ..

C =E= 1.0;

8) Establecimiento de que ecuaciones componen un modelo en particular
MODEL UNDEQN /ALL/;
9) Valores de parámetros, estimados iniciales, límites de lasvariables
Q = 1.0;
TETA.L('1')= 1.05;
10) Llamado a la técnica de solución de acuerdo al tipo de problema
SOLVE UNDEQN USING NLP MINIMIZING C;

Optimización de Procesos

37

XXI Seminario de Ingeniería Química

2003

Instituto Tecnológico de Celaya

Ejemplos Ilustrativos del Uso de GAMS
1. Para el sistema de extracción mostrado en la Figura, utilice el sistema de
modelación GAMS para determinar losvalores de las variables W1 y x1 que
maximizan la función:

Q( xF − x1 ) − λW1
donde λ = 0.05. Considere que la relación de equilibrio entre y1 y x1 está
dada por la expresión:

y1 =

Hx1
(H − 1)x1 + 1

Use un valor de H = 1.2. Note también que el balance de masa en el
sistema resulta en la ecuación:

Qx F = Qx1 + Wy1

W1
y0 = 0

Q=1000 lb/hr
xF = 0.2

Etapa de
Extracción

Q=1000 lb/hr
x1

W1
y1Figura

Optimización de Procesos

38

XXI Seminario de Ingeniería Química

2003

Instituto Tecnológico de Celaya

Código GAMS del Ejemplo 1
$TITLE EXTRACCION
$OFFSYMXREF
$OFFSYMLIST
*
*DEFINICION DE VARIABLES, PARAMETROS Y ECUACIONES
*
VARIABLES F;
POSITIVE VARIABLES X1, Y1, W1;
PARAMETERS Q, XF, LAMBDA, H;
EQUATIONS MASBAL, EQUILIBRIO, OBJ;
*
*ECUACIONES
*
MASBAL.. Q * XF =E= Q * X1 + W1 * Y1;EQUILIBRIO.. Y1 =E= (H * X1)/(((H - 1.0) * X1) + 1.0);
OBJ.. F =E= Q * ( XF -X1) - LAMBDA * W1;
*
* DEFINICION DE LAS ECUACIONES QUE FORMAN PARTE DEL MODELO
*
MODEL EXTRACTOR /ALL/;
*
* ASIGNACION DE VALORES A LOS PARAMETROS
*
Q = 1000;
XF = 0.2;
LAMBDA = 0.05;
H = 1.2;
*
* LIMITES Y VALORES INICIALES
*
Y1.L = 0.1;
Y1.UP = 1.0;
X1.L = 0.1;
X1.UP = 0.2;
W1.L = 500;
OPTION LIMROW=0;
OPTION LIMCOL=0;
*
*LLAMADO A LA TECNICA DE SOLUCION
*
SOLVE EXTRACTOR USING NLP MAXIMIZING F;

Optimización de Procesos

39

XXI Seminario de Ingeniería Química

2003

Instituto Tecnológico de Celaya

Resultados GAMS del Ejemplo 1
COMPILATION TIME
116
Model Statistics

=

0.000 SECONDS

0.7 Mb

WIN194-

SOLVE EXTRACTOR USING NLP FROM LINE 55

MODEL STATISTICS
BLOCKS OF EQUATIONS
BLOCKS OF VARIABLES
NON ZEROELEMENTS
DERIVATIVE POOL
CODE LENGTH

3
4
8
5
40

SINGLE EQUATIONS
SINGLE VARIABLES
NON LINEAR N-Z
CONSTANT POOL

3
4
3
10

GENERATION TIME
116

=

0.110 SECONDS

1.9 Mb

WIN194-

EXECUTION TIME
116

=

0.110 SECONDS

1.9 Mb

WIN194-

S O L V E
MODEL
TYPE
SOLVER

S U M M A R Y

EXTRACTOR
NLP
CONOPT

**** SOLVER STATUS
**** MODEL STATUS
**** OBJECTIVE VALUE

OBJECTIVE
DIRECTION
FROM LINE
1 NORMAL...