- Ràng buộc biến
PHỤ LỤC Chương trình tối ưu trên Matlab
Chương trình tối ưu trên Matlab
Với thép 9XC
function M=timnghiem1(x) % Ham do nham be mat
R=1.1252*x(1).^0.1833*x(2).^0.2572*x(3).^0.4484*x(4).^-0.186; % Cac bien
% x(1) luong chay dao doc Sd
% x(2) toc do quay cua chi tiet nw % x(3) chieu sau cat t
% x(4) do cung vat leu gia cong HRC % x(1)>=0.3 % x(1)<=0.5 % x(2)>=100 % x(2)<=200 % x(3)>=0.0025 % x(3)<=0.025 %x(4)>=40 %x(4)<60
% Gioi han ve do nham be mat R1=0.74;
R<=R1;
% Ham luc cat Py
Py=22.31*x(1).^0.8608*x(3).^0.3431*x(4).^0.4712; % Gioi han ve luc cat Py
P1=100; Py<=P1; Py<=P1;
% Ham luc cat Pz
Pz=6.887*x(1).^0.8623*x(3).^0.3635*x(4).^0.5315; % Gioi han ve luc cat Pz
P2=50; Pz<=P2; Pz<=P2;
% Ham rung dong
A=0.9564*x(1).^0.298*x(2).^0.2208*x(3).^0.1717; % Gioi han ve bien do rung dong
A1=2; A<=A1; A<=A1;
% Ham chi tieu boc tach vat lieu Q=3.14*30*x(2)*x(3)/60;
% Gioi han ve chi tieu boc tach vat lieu Q1=0.26;
Q>=Q1;
% Ham da muc tieu M=0.5*R/R1-0.5*Q/Q1; % M -> Min
end
Với thép 40X
function M=timnghiem2(x) % Ham do nham be mat
% Cac bien
% x(1) luong chay dao doc Sd
% x(2) toc do quay cua chi tiet nw % x(3) chieu sau cat t
% x(4) do cung vat leu gia cong HRC % x(1)>=0.3 % x(1)<=0.5 % x(2)>=100 % x(2)<=200 % x(3)>=0.0025 % x(3)<=0.025 %x(4)>=40 %x(4)<60
% Gioi han ve do nham be mat R1=0.77;
R<=R1; (adsbygoogle = window.adsbygoogle || []).push({});
% Ham luc cat Py
Py=28.663*x(1).^0.8682*x(3).^0.3556*x(4).^0.3739; % Gioi han ve luc cat Py
P1=100; Py<=P1; Py<=P1;
% Ham luc cat Pz
Pz=10.837*x(1).^0.8133*x(3).^0.3585*x(4).^0.3918; % Gioi han ve luc cat Pz
P2=50; Pz<=P2; Pz<=P2;
% Ham rung dong
A=0.9555*x(1).^0.3451*x(2).^0.2309*x(3).^0.1696; % Gioi han ve bien do rung dong
A1=2; A<=A1; A<=A1;
% Ham chi tieu boc tach vat lieu Q=3.14*30*x(2)*x(3)/60;
% Gioi han ve chi tieu boc tach vat lieu Q1=0.26;
Q>=Q1;
% Ham da muc tieu M=0.5*R/R1-0.5*Q/Q1; % M -> Min
end
Với thép 65
function M=timnghiem3(x) % Ham do nham be mat
R=1.8515*x(1).^0.1453*x(2).^0.2545*x(3).^0.4337*x(4).^-0.1518; % Cac bien
% x(1) luong chay dao doc Sd
% x(2) toc do quay cua chi tiet nw % x(3) chieu sau cat t
% x(4) do cung vat leu gia cong HRC % x(1)>=0.3
% x(1)<=0.5 % x(2)>=100 % x(2)<=200 % x(3)>=0.0025
% x(3)<=0.025 %x(4)>=40 %x(4)<60
% Gioi han ve do nham be mat R1=0.75;
R<=R1;
% Ham luc cat Py
Py=28.777*x(1).^0.8913*x(3).^0.3763*x(4).^0.3989; % Gioi han ve luc cat Py
P1=100; Py<=P1; Py<=P1;
% Ham luc cat Pz
Pz=11.498*x(1).^0.8412*x(3).^0.3779*x(4).^0.4096; % Gioi han ve luc cat Pz
P2=50; Pz<=P2; Pz<=P2; (adsbygoogle = window.adsbygoogle || []).push({});
% Ham rung dong
A=0.8746*x(1).^0.3309*x(2).^0.2226*x(3).^0.1427; % Gioi han ve bien do rung dong
A1=2; A<=A1; A<=A1;
% Ham chi tieu boc tach vat lieu Q=3.14*30*x(2)*x(3)/60;
% Gioi han ve chi tieu boc tach vat lieu Q1=0.26;
Q>=Q1;
% Ham da muc tieu M=0.5*R/R1-0.5*Q/Q1; % M -> Min
end
Với thép P18
function M=timnghiem4(x) % Ham do nham be mat
R=2.2463*x(1).^0.1559*x(2).^0.2204*x(3).^0.4277*x(4).^-0.1605; % Cac bien
% x(1) luong chay dao doc Sd
% x(2) toc do quay cua chi tiet nw % x(3) chieu sau cat t
% x(4) do cung vat leu gia cong HRC % x(1)>=0.3 % x(1)<=0.5 % x(2)>=100 % x(2)<=200 % x(3)>=0.0025 % x(3)<=0.025 %x(4)>=40 %x(4)<60
% Gioi han ve do nham be mat R1=0.74;
R<=R1;
% Ham luc cat Py
Py=36.445*x(1).^0.9324*x(3).^0.3507*x(4).^0.3373; % Gioi han ve luc cat Py
P1=100; Py<=P1; Py<=P1;
% Ham luc cat Pz
Pz=13.397*x(1).^0.8176*x(3).^0.3655*x(4).^0.3662; % Gioi han ve luc cat Pz
P2=50; Pz<=P2; Pz<=P2;
% Ham rung dong
A=1.2409*x(1).^0.3397*x(2).^0.1908*x(3).^0.1809; % Gioi han ve bien do rung dong
A1=2; A<=A1; A<=A1;
% Ham chi tieu boc tach vat lieu Q=3.14*30*x(2)*x(3)/60;
% Gioi han ve chi tieu boc tach vat lieu Q1=0.26;
Q>=Q1;
% Ham da muc tieu M=0.5*R/R1-0.5*Q/Q1; % M -> Min
(adsbygoogle = window.adsbygoogle || []).push({});