📄 11-主轴支承静不定结构计算.m__.htm
字号:
<html><head><title>MATLAB.rar 11-主轴支承静不定结构计算.m</title>
<LINK href="/inc/read_style.css" type=text/css rel=stylesheet></head>
<body>
<p><a href=http://www.pudn.com>www.pudn.com</a> > <a href="http://www.pudn.com/downloads111/sourcecode/graph/detail461284.html">MATLAB.rar</a> > 11-主轴支承静不定结构计算.m</p><!-- saved from http://www.pudn.com -->
<script src="/inc/gg_read1.js"></script><BR>
<pre name="code" class="m">
% 主轴支承静不定结构的计算(调用zzjbd.m)
disp ' '
disp ' ******** 主轴支承静不定结构的计算 ********'
disp ' '
ZCLX = input(' 滚动轴承类型：角接触球轴承-"Q";圆锥滚子轴承-"Z" == ');
disp ' '
if ZCLX=='Q'
disp ' ======== 角接触球轴承 ========'
lmd=0.3636;
Qe = input(' 载荷转换参数 e = ');
ctgalf=1.25/Qe;
m=0.3608;
n=0.6622;
h=0.32;
elseif ZCLX=='Z'
disp ' ======== 圆锥滚子轴承 ========'
lmd=0.3148;
Zy = input(' 轴向载荷系数 Y = ');
ctgalf=2.5*Zy;
m=0.3234;
n=0.6872;
h=0.34;
end
fprintf (1,' 相对位置系数 lmd = %3.4f \n',lmd)
fprintf (1,' 接触角余切 ctgalf = %3.4f \n',ctgalf)
fprintf (1,' 曲线拟合常数 m = %3.4f \n',m)
fprintf (1,' 曲线拟合常数 n = %3.4f \n',n)
disp ' '
% 主轴结构尺寸
l=525;a=203;b=130;c=194;d=271.31;D=180;
disp ' ======== 主轴的结构尺寸 ========'
fprintf (1,' 支承跨度 l = %3.2f mm \n',l)
fprintf (1,' 悬臂长度 a = %3.2f mm \n',a)
fprintf (1,' 前支承宽度 b = %3.2f mm \n',b)
fprintf (1,' 齿轮位置 c = %3.2f mm \n',c)
fprintf (1,' 齿轮分度圆直径 d = %3.2f mm \n',d)
fprintf (1,' 机床切削直径 D = %3.2f mm \n',D)
disp ' '
% 齿轮传动力
Qt=4839;Qr=1820;Qa=1294;fai=128;
disp ' ======== 齿轮传动力 ========'
fprintf (1,' 圆周力 Qt = %3.2f N \n',Qt)
fprintf (1,' 径向力 Qr = %3.2f N \n',Qr)
fprintf (1,' 轴向力 Qa = %3.2f N \n',Qa)
fprintf (1,' 载荷方位角 fai = %3.2f 度 \n',fai)
disp ' '
% 机床切削力
Pz=6978;Py=2791;Px=4400;
disp ' ======== 机床切削力 ========'
fprintf (1,' 主切削力 Pz = %3.2f N \n',Pz)
fprintf (1,' 径向力 Py = %3.2f N \n',Py)
fprintf (1,' 轴向力 Px = %3.2f N \n',Px)
disp ' '
% 计算力矩参数
hd=pi/180;
u=-Py*(l+a)+0.5*Px*D-(l-c)*(Qr*cos(fai*hd)-Qt*sin(fai*hd))-0.5*Qa*d*cos(fai*hd);
v=-Pz*(l+a)+(l-c)*(Qr*sin(fai*hd)+Qt*cos(fai*hd))+0.5*Qa*d*sin(fai*hd);
w=sqrt(u^2+v^2);
disp ' ======== 力矩参数 ========'
fprintf (1,' XOY面参数 u = %3.0f Nmm \n',u)
fprintf (1,' ZOY面参数 v = %3.0f Nmm \n',v)
fprintf (1,' 综合力矩参数 w = %3.0f Nmm \n',w)
Fa=Px-Qa;
b1b=lmd*Fa*l*ctgalf/(w-lmd*Fa*ctgalf);
if b1b<=h
Fr=lmd*Fa*ctgalf/b1b;
elseif b1b>h
mu=w/(n*b);
lo=l/(n*b)+m/n+log(Fa*ctgalf);
x0=w/(l+n*b); % 前支承径向载荷Fr初值
Fr=fsolve(@zzjbd,x0); % 利用函数fsolve解非线性方程
end
b1=b1b*b; % 前支承径向载荷作用位置
% 计算后支承径向载荷
Fry=u/(l+b1); % 前支承XOY面-径向载荷
Frz=v/(l+b1); % XOZ面-径向载荷
Fray=-Qr*cos(fai*hd)+Qt*sin(fai*hd)-Fry-Py; % 后支承XOY面-径向载荷
Fraz=Qr*sin(fai*hd)+Qt*cos(fai*hd)-Frz-Pz; % XOZ面-径向载荷
Fra=sqrt(Fray^2+Fraz^2); % 后支承径向载荷
disp ' '
disp ' ======== 主轴载荷和作用位置 ========'
fprintf (1,' 主轴轴向载荷 Fa = %3.0f N \n',Fa)
fprintf (1,' 相对位置参数 b1b = %3.4f \n',b1b)
fprintf (1,' 前支承B-径向载荷 FrB = %3.0f N \n',Fr)
fprintf (1,' 前支承径向载荷作用位置 b1 = %3.2f mm \n',b1)
fprintf (1,' 后支承A-径向载荷 FrA = %3.0f N \n',Fra)
% 主轴支承静不定结构超越非线性方程
function f=zzjbd(x)
global mu lo % 定义全局变量
f=log(x)-mu/x-lo;
</pre>
<script src="/inc/gg_read2.js"></script><BR>
<script src="http://s117.cnzz.com/stat.php?id=1236358&web_id=1236358&show=pic" language="JavaScript" charset="gb2312"></script>
</body></html>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -