铣削稳定性分析

tic
clc;close all;clear
mux1 = [600,0.02,12500,800e6,2];% wn(Hz) fai k(N/mm) kt(N/mm) ae(mm)
mux2 = [600,0.05,12500,800e6,2];
mux3 = [600,0.02,25000,800e6,2];
mux4 = [600,0.05,25000,800e6,2];
sigmax1 = 0.005*mux1;
sigmax2 = 0.005*mux2;
sigmax3 = 0.005*mux3;
sigmax4 = 0.005*mux4;
%% 第一种情况
[zhuansu_n,aplim_n] = ae_in1(mux1(1),mux1(2),mux1(3),mux1(4),mux1(5));% 求解基本的flimap
disp('请输入抽样次数')
number=input('');%1e7
w1 = normrnd(mux1(1),sigmax1(1),number,1);
fai1 = normrnd(mux1(2),sigmax1(2),number,1);
k1 = normrnd(mux1(3),sigmax1(3),number,1);
kt1 = normrnd(mux1(4),sigmax1(4),number,1);
ae1 = normrnd(mux1(5),sigmax1(5),number,1);
num = find(zhuansu_n >= 5e3 & zhuansu_n <= 7e3);
zhuansu_base1 = zhuansu_n(num(1):num(end));
aplim_base = aplim_n(num(1):num(end))*1e3;
figure('Name','基本转速1','NumberTitle','off')
plot(zhuansu_base1,aplim_base)
aplim_chazhi = zeros(length(num),number);
for i = 1 : number
    aplim_chazhi(:,i) = ae_in2(w1(i),fai1(i),k1(i),kt1(i),ae1(i),zhuansu_base1,aplim_base);
end
reliability1 = zeros(length(num),1);
for i = 1:length(num)
    reliability1(i) = length(find(aplim_chazhi(i,:)>=aplim_base(i)))/number;
end
figure('Name','可靠性指数1','NumberTitle','off')
plot(zhuansu_base1,reliability1)
toc

%%------绘制刚度对切削深度的影响-------------
clc;close all;
ae=0.1;  %径向切深（mm）
disp('请选择逆铣还是顺铣，0代表逆铣，1代表顺铣')
flag=input('');
[zhuansu_n,aplim_n]=ae_in(ae,flag);
%% 
figure(2)
plot(zhuansu_n1,aplim_n1,'--',zhuansu_n2,aplim_n2,'*',zhuansu_n3,aplim_n3,'+') 
title({'薄壁件高速铣削三维稳定性'});
xlabel({'主轴转速（rpm）'});
ylabel({'轴向切深（mm）'});
grid on
ylim([0 0.3])
xlim([2000 5000])
set(gcf,'color','w');
legend('k=100e3','k=200e3','k=300e3')