车削端面的表面微观形貌仿真（matlab）

本文主要是介绍车削端面的表面微观形貌仿真（matlab），希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

代码段禁止出现在毕业设计、商业活动中！仅供交流。

代码段

相关参数已在代码中标注，由于是自己参考论文编写的并且才接触matlab不到一周，因此采用的比较多的循环结构，程序某些地方可能多余（懒得删除），而且程序运行比较慢，大家可以参考自己进行优化。

%% initialization
clear;
close all;
clc;%% known conditions
n = 500; %Spindle speed of lathe，r/min		主轴转速
Rdt = 1.554; %Radius of arc knife，mm			刀尖圆弧半径
s = 15; %Spindle feed speed，mm/min			主轴进给速度
h = 2/1000; %Depth of cut，mm					加工深度%% Assumptions
Rw = 0.5;%Workpiece radius，mm				工件半径
Lx = 0.25;%The length of the simulation area，mm		加工长度
Ly = 0.25;%Width of simulation area，mm						加工宽度
Sf = s/n; %Feed per revolution,mm/r			进给速度（注意单位）
mx = Sf/10;
my = Sf/10;
N = ceil(Lx/mx);
M = ceil(Rw/Sf);% Overall surface roughness
Rt = Rdt - sqrt(Rdt^2-(Sf^2)/4);			%总的表面粗糙度%% Set matrix
AA = [];
k = 0:1:M;
rk = [];%% Calculate profile height
for i = 0:Nfor j = 0:N%Cartesian coordinates of pointsXij = i*mx-Lx/2;Yij = j*my-Ly/2;%Point polar coordinatesr = sqrt(Xij^2+Yij^2);theta = atan(Yij/Xij);% jidge thetaif theta < 0theta = theta +2*pi;elseif theta >= 2*pitheta = theta -2*pi;elsetheta = theta;end%compute the height of the pointhk = [];for k = 1:Mif i < N/2rk = Rw - Sf*(k + theta/(2*pi));elserk = Rw - Sf*((k-0.5) + theta/(2*pi));endhk1 = Rdt*(1-(sqrt(1-((r-rk)/Rdt)^2)));hk = [hk,hk1];endhklow = min(min(hk));% Judge height and processing depthif hklow > hhklow = h;end% save the min heightAA = [AA,hklow];end
end
%List the coordinate range
Ax = linspace(0,0.25,N+1);
Ay = Ax;% Into a square
AA = reshape(AA,[N+1,N+1]);% Drawing
surf(Ay,Ax,AA);
title('Surface microtopography');
xlabel('x-axis');
ylabel('y-axis');
zlabel('z-axis');