错误81和92，希望哪位大神有类似情况，给个查错方向，谢谢！！

lj151019

model:
sets:
w/1..4/:ETi,LTi,q,t3;!第1个是配送中心，2-3是客户，有2个，4是充电站,1个、需求量;
truck/1..2/;!假设的三辆车、最大载重，最大电量;
link(w,w):d;!两点间的距离;
link1(w,w,truck):x;!0-1决策变量;
link2(w,truck):y;!0-1变量，k车是否服务于i点;
link3(w,truck):t1,t2,st,B1,B2,Q1,Q2;!车辆k到达i的时刻，车辆k离开点i的时刻,车辆k开始工作的时刻，车辆到达i处的剩余电能，离开i处的电能，Q表示k车在i处的剩余载重量;
link4(w,w,truck):tijk;!车辆k从点i到点j的行驶时间;
link5(w,truck):tf,B12;!第三个决策变量，车辆k在充电站的充电时长和电量，是个连续变量;
endsets
data:
a=1.6;
c1=200;
e1=0.6;
e21=0.8;
e22=36;
Qmax=1.5;
Bmax=80;
ETi=
0        0.2        0.3        0       
;
LTi=
15        1.4        1.9        15       
;
r=44;
p=100;
n=0.14;
m1=100;m2=300;
v=60;
d=
0        392.5        260        565       
392.5        0        1062.5        1878.5       
260        1062.5        0        221       
565        1878.5        221        0       
;!;
q=
0        0.2        0.37        0       
;
t3=
0        1        1.5        0       
;!卸货时间，已知已求,默认90;
k=1;!前面k不可以的原因是因为这里赋值了2，所以不可以作为变量出现在后面的程序当中;
enddata

!目标函数;
min=f1+f2+f3+f4+f5;
f1=c1*@sum(truck(b)@sum(w(j)|j#gt#1:x(1,j,b))));
f2=a*@sum(truck(b)@sum(link(i,j)|i#ne#j:e1*d(i,j)*x(i,j,b))));
f3=f31+f32;
f31=a*@sum(truck(b)@sum(link(i,j)|i#ne#j:e21*d(i,j)*x(i,j,b)*Q1(j,b))));
f32=a*@sum(truck(b)sum(w(i)|i#gt#1smax(ETi(i)-t1(i,b),0)+t3(i)));
f4=p*@sum(truck(b)sum(link(i,j)|i#ne#j1-@exp(-n*((d(i,j)/v)+t3(j))))*x(i,j,b)));
!f5=@sum(w(i)sum(link3(i,k):m1*@smax(t1(i,k)-Ei(i),0)))+@for(w(i)sum(link3(i,k):m2@smax(Li(i)-t1(i,k),0)));
f5=@sum(truck(b)sum(w(i):m1*@smax(t1(i,b)-ETi(i),0)+m2*@smax(t1(i,b)-LTi(i),0)));

!约束条件;
!!网络平衡参数;
@for(w(m)|m#eq#1 #or# m#eq#4:
        @sum(truck(b)sum(link1(i,m,b)|i#ne#m:x(i,m,b)))=
        @sum(truck(b)sum(link1(m,j,b)|m#ne#j:x(m,j,b))));!充电站/配送中心的流量平衡，但是车辆的访问次数不仅一次;

@for(w(m)|m#ne#1 #and# m#ne#4:
        @sum(truck(b)sum(link1(i,m,b)|i#ne#m:x(i,m,b)))=
        @sum(truck(b)sum(link1(m,j,b)|m#ne#j:x(m,j,b))));!客户的流量平衡，只能被访问一次;
@for(w(m)|m#ne#1 #and# m#ne#4:
        @sum(truck(b)sum(link1(i,m,b)|i#ne#m:x(i,m,b)))=1);

@for(w(i)|i#eq#1:
        @sum(w(j)|j#gt#1:
                @sum(truck(b):x(i,j,b)))<=3);!配送中心出去的车辆数，！！这里有个假设和隐形约束，车辆只有将自己的任务送完之后才可以返回配送中心;

!!车辆容量约束;
@for(truck(b)sum(w(i):q(i)*y(i,b))<=Qmax);!每辆车的容载约束;
@for(truck(b):
        @for(link1(i,j,b)|i#ne#j2(j,b)<=Q1(i,b)-q(i)*x(i,j,b)+Qmax*(1-x(i,j,b))));!车辆k在客户j处的剩余载重量与在上一点客户i处的关系;
                       
!!电能约束;
!从配送中心，充电站，客户点离开时的电量，;
@for(w(i)|i#eq#1for(truck(b):B2(i,b)=Bmax*y(i,b)));
@for(w(i)|i#eq#4for(truck(b):B2(i,b)=B1(i,b)+tf(i,b)*r));
@for(w(i)|i#ne#1 #and# i#ne#4:
        @for(truck(b):
                        B2(i,b)=B1(i,b)-e22*(@smax(ETi(i)-t1(i,b),0)+t3(i))*y(i,b)));
!从客户，充电站到达配送中心时的电量;
@for(truck(b):
        @for(link5(i,b)|i#ne#1:
                B1(1,b)<=B1(i,b)-e22*(@smax(ETi(i)-t1(i,b),0)+t3(i))*y(i,b)*x(i,1,b)-e1*d(i,1)*x(i,1,b)+Bmax*(1-x(i,1,b))));

@for(truck(b):
        @for(link1(i,j,b)|(i#eq#4) #and# j#eq#1:
                B1(j,b)<=B2(i,b)-e1*d(i,j)*x(i,j,b)+Bmax*(1-x(i,j,b))));

!从客户i到充电站/客户j的电量条件;
@for(truck(b):
        @for(link1(i,j,b)|(i#ne#1 #and# i#ne#4) #and# j#ne#1 #and# i#ne#j:
                B1(j,b)<=B1(i,b)-(e22*(@smax(ETi(i)-t1(i,b),0)+t3(i))+e1*d(i,j)+e21*d(i,j)*Q1(j,b))*x(i,j,b)+Bmax*(1-x(i,j,b))));

!从充电站i到客户j的电量;
@for(truck(b):
        @for(link1(i,j,b)|(i#eq#4 #and# j#ne#1 #and# j#ne#4):
                B1(j,b)<=B2(i,b)-e1*d(i,j)*x(i,j,b)-e21*d(i,j)*Q1(j,b)*x(i,j,b)+Bmax*(1-x(i,j,b))));

!从配送中心i到客户j的电量条件;
@for(truck(b):
        @for(link1(i,j,b)|i#eq#1 #and# j#ne#1 #and# j#ne#4:
                B1(j,b)<=Bmax-e1*d(i,j)*x(i,j,b)-e21*d(i,j)*Q1(j,b)*x(i,j,b)+Bmax*(1-x(i,j,b))));

!每个节点的剩余电量>下一配送路程上所消耗的电量;
@for(truck(b):
        @for(link1(i,j,b)|j#eq#1 #or# j#eq#4:
                B2(i,b)>=e1*d(i,j)*x(i,j,b)+e21*d(i,j)*Q1(j,b)*x(i,j,b)));

@for(truck(b):
        @for(link1(i,j,b)|j#ne#1 #and# j#ne#4:
                B2(i,b)>=e1*d(i,j)*x(i,j,b)+e21*d(i,j)*Q1(j,b)*x(i,j,b)+e22*(@smax(ETi(i)-t1(i,b),0)+t3(i))*x(i,j,b)));
!时间窗约束;
@for(truck(b):
        @for(link3(i,b)|i#eq#4:st(i,b)=t1(i,b)));

@for(truck(b):
        @for(link3(i,b)|i#ne#1 #and# i#ne#4:st(i,b)=t1(i,b)+@smax(ETi(i)-t1(i,b),0)));

@for(truck(b):
        @for(link1(i,j,b)|i#eq#4:t2(j,b)=t1(i,b)+(B2(i,b)-B1(i,b))/r+d(i,j)*x(i,j,b)/v));

@for(truck(b):
        @for(link1(i,j,b)|i#ne#1 #and# i#ne#4:t2(j,b)=st(i,b)+(d(i,j)*x(i,j,b))/v));

@for(link1bin(x));
@for(link2bin(y));

end

lj151019

上面的表情，都是@

lj151019

无解原因，与设置的数据有关。