Hello, I have similar problem, gams shows me **** MODEL STATUS 10 Integer Infeasible.
SET
I clientes/0*5/
N0(i) deposito /0/
Nd(i) N demanda(TT)conjunto de nodos de demanda no atendidos en el momento TT /1*5/
Ns(i) N serviced(TT) conjunto de nodos que se sirven en el TT
Nar(i) N artificial(TT)conjunto de nodos intermedios artificiales para vehiculos que se dirigen en el momento TT
N0d(i) N0+Nd
Nds(i) Nd+Ns
N0sar(i) N0+Ns+Nar
Ndsar(i) Nd+Ns+Nar
N0dsar(i) N0+Nd+Ns+Nar
k vehiculos /1*3/
t tiempos /0*660/
Alias(i,j)
;
Ns(i)=no;
Nar(i)=no;
N0d(i)=N0(i)+Nd(i);
Nds(i)=Nd(i)+Ns(i);
N0sar(i)=N0(i)+Ns(i)+Nar(i);
Ndsar(i)=Nd(i)+Ns(i)+Nar(i);
N0dsar(i)=N0(i)+Nd(i)+Ns(i)+Nar(i);
PARAMETERS
q(i) demanda del nodo i(articulos)
/1 115
2 165
3 150
4 63
5 30/
QQ(k) capacidad del vehiculo (articulos)
/1 300
2 300
3 300/
QQTT(k) carga del vehiculo K en el momento TT(articulos)
/1 300
2 300
3 300/
e(i) ventana de tiempo de inicio en el nodo i
/0 0
1 60
2 60
3 120
4 0
5 180/
l(i) finalizacion de la ventana de tiempo en el nodo i
/0 660
1 600
2 540
3 600
4 300
5 600/
s(i) tiempo de servicio en el nodo i
/0 45
1 20
2 35
3 90
4 20
5 40/
;
TABLE
tTT(i,j) estimacion del tiempo de viaje del enlace entre los nodos i y j en el momento TT (minutos)
0 1 2 3 4 5
0 0 34 40 28 32 24
1 34 0 36 28 12 25
2 40 32 0 51 30 24
3 25 30 48 0 32 33
4 32 9 32 32 0 17
5 23 20 27 32 17 0
;
SCALAR
TT hora de la ultima determinación o actualización de rutas /0/
*k(i) vehiculo que en el momento TT se encuentra en el nodo i+Ns+Nar/1/
M un numero positivo muy grande/660/
;
POSITIVE VARIABLE
a(i,k) inicio del servicio en el nodo i por el vehiculo k ;
free variable
z variable de la funcion objetivo;
binary variable
x(i,j,k,t) es igual a 1 si el vehiculo k se aleja del nodo i al nodo j en el momento t de lo contrario es igual a cero;
Equations
FO,eq2,eq3,eq4,eq5,eq8,eq9,eq10,eq11a,eq11b,eq12;
FO.. z=e= SUM((i,j,k,t)${N0dsar(i) and N0d(j)},tTT(i,j)*x(i,j,k,t));
eq2(i)$Ndsar(i).. sum((j,k,t)$N0d(j),x(i,j,k,t)) =e= 1;
eq3(j,k)$Nd(j).. sum((i,t)$N0dsar(i),x(i,j,k,t))-sum((i,t)$N0d(i),x(j,i,k,t)) =e= 0;
eq4(k).. sum((j,t)$Nd(j),x('0',j,k,t)) =l= 1;
eq5(k).. sum((i,j,t)${N0sar(i) and N0d(j)},x(i,j,k,t))-sum((i,t)$Ndsar(i),x(i,'0',k,t)) =e= 0;
*eq6(i)$Nar(i).. sum((j)$N0d(j),x(i,j,k,'0'))=e=1;
*eq7(i)$Ns(i).. sum((j)$N0d(j),x(i,j,ks(i)))=e=1;
eq8(i,j,k,t)${N0dsar(i) and N0d(j)}.. a(i,k)+s(i)+tTT(i,j)-a(j,k)=l=(1-x(i,j,k,t))*M;
eq9(i,k,t)$Nds(i).. e(i)* sum((j)$N0d(j),x(i,j,k,t))=l=a(i,k);
eq10(i,k,t)$Nds(i).. a(i,k)=l=(l(i)-s(i))*sum((j)$N0d(j),(x(i,j,k,t)));
eq11a(k).. e('0')=l=a('0',k);
eq11b(k).. a('0',k)=l=l('0');
eq12(k).. QQTT(k)+sum((i,j,t)${Nds(i) and N0d(j)},(q(i)*x(i,j,k,t)))=l=QQ(k);
MODEL MODELO MATEMATICO VRP DINAMICO /all/;
SOLVE MODELO using MIP minimizing z;
the error shows:
EMAINING 71385 ENTRIES SKIPPED
GAMS Rev 237 WEX-VS8 23.7.3 x86/MS Windows 03/05/21 18:41:47 Page 4
G e n e r a l A l g e b r a i c M o d e l i n g S y s t e m
Model Statistics SOLVE MODELO Using MIP From line 127
MODEL STATISTICS
BLOCKS OF EQUATIONS 11 SINGLE EQUATIONS 91,254
BLOCKS OF VARIABLES 3 SINGLE VARIABLES 71,407
NON ZERO ELEMENTS 626,635 DISCRETE VARIABLES 71,388
GENERATION TIME = 2.980 SECONDS 37 Mb WIN237-237 Aug 23, 2011
EXECUTION TIME = 2.980 SECONDS 37 Mb WIN237-237 Aug 23, 2011
GAMS Rev 237 WEX-VS8 23.7.3 x86/MS Windows 03/05/21 18:41:47 Page 5
G e n e r a l A l g e b r a i c M o d e l i n g S y s t e m
Solution Report SOLVE MODELO Using MIP From line 127
S O L V E S U M M A R Y
MODEL MODELO OBJECTIVE z
TYPE MIP DIRECTION MINIMIZE
SOLVER CPLEX FROM LINE 127
**** SOLVER STATUS 1 Normal Completion
**** MODEL STATUS 10 Integer Infeasible
**** OBJECTIVE VALUE 0.0000
RESOURCE USAGE, LIMIT 0.562 1000.000
ITERATION COUNT, LIMIT 0 2000000000
IBM ILOG CPLEX Jul 14, 2011 23.7.3 WIN 27723.27726 VS8 x86/MS Windows
--- GAMS/Cplex licensed for continuous and discrete problems.
Cplex 12.3.0.0
MIP status(119): integer infeasible or unbounded
Problem is integer infeasible.
Unexpected solution discarded