请高手帮忙解决lingo问题

buct

min=2*3.1415926*(t0-ta)/(2/(d1*hin)+1/lamuda0*@log(d2/d1)+1/lamuda1*@log(d3/d2)+1/lamuda2*@log(d4/d3)+1/lamuda3*@log(d5/d4)+1/lamuda4*@log(d6/d5)+2/(d6*hout))*T*Pe*3600*L/10^9+dan*(Qb+Qf+Qs);
T=8000;
Pe=26;
lamuda0=45;
lamuda1=0.03;
lamuda2=0.037;
lamuda3=0.072;
lamuda4=0.1;
Nu=0.0023*re^0.8*pr^0.3*(t0-t1);
pr=0.70681-2.08938*10^(-4)*t0+3.43447*10^(-7)*t0^2;
t0=350+273.15;
ta=25+273.15;
re=d1*um/v;
um=25;
d1=0.108;
d2=0.114;
L=1000;
B1=15000;
B2=950;
B3=1400;
B4=20;
q1=6000;
W=10;
v=-2.40084*10^(-5)+1.2834*10^(-7)*t0;
lamuda=0.00648+6.79112*10^(-5)*t0;
hin=Nu*lamuda/d1;
hout=1.163*(10+@sqrt(W));
!(t0-t1)*(3.1415926*d1*hin)=2*3.1415926*(t0-ta)/(2/(d1*hin)+1/lamuda0*@log(d2/d1)+1/lamuda1*@log(d3/d2)+1/lamuda2*@log(d4/d3)+1/lamuda3*@log(d5/d4)+1/lamuda4*@log(d6/d5)+2/(d6*hout));
(t0-t1)*(3.1415926*d1*hin)=2*3.1415926*lamuda0*(t1-t2)/@log(d2/d1);
2*3.1415926*lamuda0*(t1-t2)/@log(d2/d1)=2*3.1415926*lamuda1*(t2-t3)/@log(d3/d2);
2*3.1415926*lamuda1*(t2-t3)/@log(d3/d2)=2*3.1415926*lamuda2*(t3-t4)/@log(d4/d3);
2*3.1415926*lamuda2*(t3-t4)/@log(d4/d3)=2*3.1415926*lamuda3*(t4-t5)/@log(d5/d4);
2*3.1415926*lamuda3*(t4-t5)/@log(d5/d4)=2*3.1415926*lamuda4*(t5-t6)/@log(d6/d5);
2*3.1415926*lamuda4*(t5-t6)/@log(d6/d5)=(t6-ta)*(3.1415926*d6*hout);
Qb=(3.1415926*d3^2/4-3.1415926*d2^2/4)*L*B1+(3.1415926*d4^2/4-3.1415926*d3^2/4)*B2*L+(3.1415926*d5^2/4-3.1415926*d4^2/4)*B3*L;
Qf=q1+3.1415926*((d5+d6)/2)*B4*L;
Qs=beita*(Qb+Qf);
beita=0.1;
2*3.1415926*(t0-ta)/(2/(d1*hin)+1/lamuda0*@log(d2/d1)+1/lamuda1*@log(d3/d2)+1/lamuda2*@log(d4/d3)+1/lamuda3*@log(d5/d4)+1/lamuda4*@log(d6/d5)+2/(d6*hout))<=186;
n=10;
i=0.11;
dan=i*(1+i)/((1+i)^n-1);
(d3-d2)/2<=0.1;
d3-d2>0.001;
d4-d3>0.0010;
d5-d4>0.0010;
d6>d5+0.001;
d2>d1+0.001;
(d4-d3)/2<=0.1;
(d5-d4)/2<=0.1;
t6<=333;
d1>0;
d2>0;
d3>0;
d4>0;
d5>0;
d6>0;
t1>0;
t2>0;
t3>0;
t4>0;
t5>0;
t6>0;
t0>t1;
t1>t2+1;
t2>t3+1;
t3>t4+1;
t4>t5+1;
t5>t6+1;
t6>ta+1;

buct

其中参数可以调整，使整个程序可以运行

buct

已解决