西松模型的MATLAB程序编制

 

 

 

程序：
%西松模型
function F=xscuttermodel(~)
a=input('请输入剪切强度值');%通过实验测量
b=input('请输入切削深度值');%可调整的数值
c=input('请输入内摩擦角');%通过实验测量
d=input('请输入刀具倾角');%可调整的数值
e=input('也称岩石切削摩擦角');%切削合力方向与刀具法线之间的夹角
F(1)=(2*a*b*cos(c))./((10.3-0.8*d)*(1-sin(c-d+e)));%切削力的值
disp('切削力F为:');
disp(F);
P=F*cos(e-d);%切向力的水平分量P
Q=F*sin(e-d);%切削力的法相分量Q
disp('切削力F水平分量，水平切削力P为：');
disp(P);
disp('切削力F法相分量，推力Q为：');
disp(Q);

作者：GaiJM