?? 07-攪拌機工作頭旋輪軌跡分析與仿真-1.m__.htm
字號:
<html><head><title>MATLAB.rar 07-攪拌機工作頭旋輪軌跡分析與仿真-1.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> > 07-攪拌機工作頭旋輪軌跡分析與仿真-1.m</p><!-- saved from http://www.pudn.com -->
<script src="/inc/gg_read1.js"></script><BR>
<pre name="code" class="m">
% 旋輪線軌跡模擬
% 圓錐齒輪傳動參數
m=2;z7=11;z8=36;
dt7=10.5;dt8=35;
bt=65.25; % 行星輪軸線與XOY平面夾角
% 攪拌桿外點旋轉半徑
l=65;
hd=pi/180;
% 圓錐齒輪幾何尺寸
r7=0.5*m*z7;
% r7v=r7/cos(dt7*hd);
r8=0.5*m*z8;
% r8v=r8/cos(dt8*hd);
k1=r7+r8;
k2=(r7+r8)/r7;
% 1----二維旋輪線參數計算
for i=1:1:360
x=k1*cos(i*hd)+l*cos(k2*i*hd);
y=k1*sin(i*hd)+l*sin(k2*i*hd);
xl1(i,:)=[i x y];
end
for i=1:1:720
x=k1*cos(i*hd)+l*cos(k2*i*hd);
y=k1*sin(i*hd)+l*sin(k2*i*hd);
xl2(i,:)=[i x y];
end
for i=1:1:1080
x=k1*cos(i*hd)+l*cos(k2*i*hd);
y=k1*sin(i*hd)+l*sin(k2*i*hd);
xl3(i,:)=[i x y];
end
for i=1:1:1440
x=k1*cos(i*hd)+l*cos(k2*i*hd);
y=k1*sin(i*hd)+l*sin(k2*i*hd);
xl4(i,:)=[i x y];
end
% 繪制二維旋輪線
figure(1); % 生成第1個圖形窗口
subplot(2,2,1); % 選擇第1個子窗口
plot(xl1(:,2),xl1(:,3))
grid % 繪制網格線
title('二維旋輪線(\theta =360度)')
subplot(2,2,2); % 選擇第2個子窗口
plot(xl2(:,2),xl2(:,3))
grid % 繪制網格線
title('二維旋輪線(\theta =720度)')
subplot(2,2,3); % 選擇第3個子窗口
plot(xl3(:,2),xl3(:,3))
grid % 繪制網格線
title('二維旋輪線(\theta =1080度)')
subplot(2,2,4); % 選擇第4個子窗口
plot(xl4(:,2),xl4(:,3))
grid % 繪制網格線
title('二維旋輪線(\theta =1440度)')
% 2----三維旋輪線參數計算
for i=1:1:360
x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
z=sqrt(x^2+y^2)*tan(bt*hd);
xlx1(i,:)=[i x y z];
end
for i=1:1:720
x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
z=sqrt(x^2+y^2)*tan(bt*hd);
xlx2(i,:)=[i x y z];
end
for i=1:1:1080
x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
z=sqrt(x^2+y^2)*tan(bt*hd);
xlx3(i,:)=[i x y z];
end
for i=1:1:1440
x=(k1*cos(i*hd)+l*cos(k2*i*hd))*cos(bt*hd);
y=(k1*sin(i*hd)+l*sin(k2*i*hd))*cos(bt*hd);
z=sqrt(x^2+y^2)*tan(bt*hd);
xlx4(i,:)=[i x y z];
end
% 矩陣cs(:,j)表示第j列的各行元素
% 繪制三維旋輪線
figure(2); % 生成第2個圖形窗口
subplot(2,2,1); % 選擇第1個子窗口
plot3(xlx1(:,2),xlx1(:,3),xlx1(:,4))
grid
title('三維旋輪線(\theta =360度)')
subplot(2,2,2); % 選擇第2個子窗口
plot3(xlx2(:,2),xlx2(:,3),xlx2(:,4))
grid
title('三維旋輪線(\theta =720度)')
subplot(2,2,3); % 選擇第3個子窗口
plot3(xlx3(:,2),xlx3(:,3),xlx3(:,4))
grid
title('三維旋輪線(\theta =1080度)')
subplot(2,2,4); % 選擇第4個子窗口
plot3(xlx4(:,2),xlx4(:,3),xlx4(:,4))
grid
title('三維旋輪線(\theta =1440度)')
% 輸出數據
disp ' '
disp ' ========== 圓錐齒輪傳動參數 ==========';
fprintf(1,' 行星輪齒數 z7 = %3.0f \n',z7);
fprintf(1,' 分度圓半徑 r7 = %3.3f 度 \n',r7);
fprintf(1,' 分度圓錐角 dt7 = %3.3f 度 \n',dt7);
fprintf(1,' 中心輪齒數 z8 = %3.0f \n',z8);
fprintf(1,' 分度圓半徑 r8 = %3.3f 度 \n',r8);
fprintf(1,' 分度圓錐角 dt8 = %3.3f 度 \n',dt8);
fprintf(1,' 中心輪與行星輪當量半徑之比 rb = %3.0f \n',r8/r7);
disp ' '
[' 行星輪公轉',' 動點x坐標',' 動點y坐標',' 動點z坐標']
[xlx1(:,1),xlx1(:,2),xlx1(:,3),xlx1(:,4)]
</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 + -