Matlab在机械设计中的应用大作业
Matlab在机械设计中的应用大作业1、试用解析综合法设计一个曲柄摇杆机构。
已知机构行程速度变化系数 1.25
k=，摇杆CD的长度
3250mm
l=，摆角30
ψ=，要求机构
的最小传动角
min
γ≥40。
确定曲柄摇杆机构各构件杆长
调用函数的编制：
function f=funct(x)
k=1.25;
theta=pi*(k-1)/(k+1);
yg=250;
pis=pi/6
gamin=2*pi/9
f1=(x(2)+x(1))^2+(x(2)-x(1))^2-2*(x(2)+x(1))*(x(2)-x(1))*cos(theta)-(2*yg*sin(pis/2))^2;
f2=yg^2+x(3)^2-2*yg*x(3)*cos(x(4))-(x(2)-x(1))^2;
f3=yg^2 +x(3)^2-2*yg*x(3)*cos(x(4)+pis)-(x(2)+x(1))^2;
f4=yg^2+x(2)^2-2*yg*x(2)*cos(gamin)-(x(3)-x(1))^2;
f=[f1;f2;f3;f4];
主函数：
>> x0=[50 120 200 0.5];
>> k=1.25;
>> theta=pi*(k-1)/(k+1);
>> yg=250;
>> gamin=2*pi/9;
>> x=fsolve(@funct,x0)
求解结果：
Equation solved.
fsolve completed because the vector of funct values is near zero
as measured by the default value of the funct tolerance, and
the problem appears regular as measured by the gradient.
x =
62.9934 105.9045 245.0702 0.1724
2、四连杆机构如图1 所示，已知各构件的尺寸L1, L2, L3, L4及原动件1的角位移θ1和等角速度ω1，求构件2和3的角位移θ 2 ,θ3，角速度ω2,ω3，角加速度ε2,ε3。
解：设L1=304.8设定机架1长度，L2=101.6设定曲柄2长度，L3=254设定连杆3长度，L4=177.8设定摇杆4长度。
% 曲柄摇杆机构运动分析
% (1)-----为了计算连杆的输出角th3和摇杆的输出角th4
% 设定各杆的长度(单位:毫米)
Ntrpss函数
funct [th3,th4]=ntrpss(th,rs)
L=sqrt(rs(1)*rs(1)+rs(2)*rs(2)-2*rs(1)*rs(2)*cos(th(1)+0.4867));
phi=asin((rs(2)*sin(th(1))+90)/L);
beta=acos((rs(4)*rs(4)+L*L-rs(3)*rs(3))/(2*rs(4)*L));
if beta<0
beta=beta+pi;
end
th4=pi-phi-beta;
th3=asin((rs(4)*sin(th4)-rs(2)*sin(th(1))-90.00)/rs(3)); end
主函数：
clear
rs(1)=304.8;
rs(2)=101.6;
rs(3)=254;
rs(4)=177.8;
dr=pi/180.0;%
% 设定初始推测的输入
% 机构的初始位置
th(1)=0.0;
th(2)=45*dr;
th(3)=135*dr;
% 摇杆4的初始位置角可以用三角形的正弦定理确定
%th(3)=pi-asin(sin(th(2))*rs(3)/rs(4))
dth=5*dr;
for i=1:72
[th3,th4]=ntrps(th,rs);
th34(i,:)=[th(1)/dr th3/dr th4/dr]; %
th(1)=th(1)+dth;
th(2)=th3;
th(3)=th4;
end
subplot(2,2,1)
plot(th34(:,1),th34(:,2),th34(:,1),th34(:,3));
axis([0 360 0 170])
grid
xlabel('时间/秒')
ylabel('从动件角位移/度')
title('角位移线图')
text(110,110,'摇杆4角位移')
text(50,35,'连杆3角位移')
% (2)- - - - -计算连杆的角速度om3和摇杆的角速度om4
% Setting initial conditions
om2=250;
T=2*pi/om2;
for i=1:72
t=i*T/72;
ct(2)=i*dth;
A=[-rs(3)*sin(th34(i,2)*dr) rs(4)*sin(th34(i,3)*dr);rs(3)*cos(th34(i,2)*dr) -rs(4)*cos(th34(i,3)*dr)];
B=[om2*rs(2)*sin(ct(2));-om2*rs(2)*cos(ct(2))];
om=inv(A)*B;
om3=om(1);
om4=om(2);
om34(i,:)=[t om3 om4];
end
subplot(2,2,2) %
plot(om34(:,1),om34(:,2),om34(:,1),om34(:,3))
axis([0 0.026 -190 210])
grid
title('角速度线图')
xlabel('时间/秒')
ylabel('从动件角速度')
text(0.001,170,' 摇杆4角速度')
text(0.013,130,' 摇杆3角速度')
% (3)- - - - - 计算连杆的角加速度a3和摇杆的角加速度a4
a2=0;
for i=1:72
t=i*T/72;
c(2)=i*dth;
C=[-rs(3)*sin(th34(i,2)*dr) rs(4)*sin(th34(i,3)*dr);rs(3)*cos(th34(i,2)*dr) -rs(4)*cos(th34(i, 3)*dr)];
D(1)=a2*rs(2)*sin(c(2))+om2^2*rs(2)*cos(c(2))+om34(i,2)^2*rs(3)*cos(th34(i,2)*dr)-om34(i,3)^ 2*rs(4)*cos(th34(i,3)*dr);
D(2)=-a2*rs(2)*cos(c(2))+om2^2*rs(2)*sin(c(2))+om34(i,2)^2*rs(3)*sin(th34(i,2)*dr)-om34(i,3) ^2*rs(4)*sin(th34(i,3)*dr);
a=inv(C)*D';
a3=a(1);
a4=a(2);
a34(i,:)=[t a3 a4];
end
subplot(2,2,3)
plot(a34(:,1),a34(:,2),a34(:,1),a34(:,3))
axis([0 0.026 -60000 80000])
grid
title('角加速度线图')
xlabel('时间')
ylabel('从动件加速度')
text(0.003,62000,'摇杆4角加速度')
text(0.010,33000,'摇杆3角加速度')
disp '曲柄转角连杆转角-摇杆转角-连杆角速度-摇杆角速度-连杆加速度-摇杆加速度'
ydcs=[th34(:,1),th34(:,2),th34(:,3),om34(:,2),om34(:,3),a34(:,2),a34(:,3)];
disp(ydcs)
3、已知图中机构以下参数：
?1）、各构件尺寸（5个长度）；L AB=80mm, L BC=260mm, L CD=300mm, L CE=100mm, L EF=460mm。
?2）、各构件重量及重心位置；G1=36N，S1在A点，G2=60N，S2在BC中点，G3=72N，S3在C点，G4=85N，S4在EF中点，G5=85N，S5在F点。
?3）、各构件对重心的转动惯量（4个转动构件）；I1=0.03kgm2，I2=0.08kgm2，I3=0.1 kgm2，I4=0.12 kgm2。
?4）、阻力：注意推程为正、回程为负；作用在滑块5上的水平阻力（生产阻力）Pr=-4000N。当回程时，Pr=50N。
?5）、曲柄顺时针回转，n1=400r/min转动的转速（注意转速单位）
用matlab求解：
1)滑块5的位移线图（s-φ）、速度线图（v-φ）、加速度线图（a-φ）。这里φ指曲柄转过的
角度，曲柄从右边水平开始，逆时针计算角度0~2π.
2)作原动件1的平衡力矩线图。
3)编写界面程序，要求当输入曲柄某一角度φ值时，可计算出所有运动副反力R12，R23,
R34，R56, R61，R63。
要求最低实现以下基本功能：
其他格式自定。
clear
rs(1)=0.19235;
rs(2)=0.08;
rs(3)=0.26;
rs(4)=0.30;
rs(5)=0.40;
rs(6)=0.460;
e=0.090;
dr=pi/180.0000;
th(1)=0.0000;
th(2)=50.3800*dr;
th(3)=75.3600*dr;
dth=dr;
for i=1:360
[th3,th4]=ntrpss(th,rs);
th5=-asin((rs(5)*sin(th4)-e)/rs(6));
s=0.170+rs(5)*cos(th4)+rs(6)*cos(th5);
th34(i,:)=[th(1)/dr th3/dr th4/dr th5/dr s];
th(1)=th(1)+dth;
th(2)=th3;
th(3)=th4;
end
om2=400;
T=2*pi/om2;
for i=1:360
t=i*T/360;
ct(2)=i*dth;
A=[-rs(3)*sin(th34(i,2)*dr) rs(4)*sin(th34(i,3)*dr);rs(3)*cos(th34(i,2)*dr) -rs(4)*cos(th34(i,3)*dr)];
B=[om2*rs(2)*sin(ct(2));-om2*rs(2)*cos(ct(2))];
om=inv(A)*B;
om3=om(1);
om4=om(2);
C=[rs(6)*sin(th34(i,4)*dr) 1;-rs(6)*cos(th34(i,4)*dr) 0];
D=[-om(2)*rs(5)*sin(th34(i,3)*dr);om(2)*rs(5)*cos(th34(i,3)*dr)];
on=inv(C)*D;
om5=on(1);
om6=on(2);
om34(i,:)=[t om3 om4 om5 om6];
end
a2=0;
for i=1:360
t=i*T/360;
c(2)=i*dth;
E=[-rs(3)*sin(th34(i,2)*dr) rs(4)*sin(th34(i,3)*dr);rs(3)*cos(th34(i,2)*dr) -rs(4)*cos(th34(i,3)*dr)];
F(1)=a2*rs(2)*sin(c(2))+om2*om2*rs(2)*cos(c(2))+om34(i,2)*om34(i,2)*rs(3)*cos(th34(i,2)*dr) -om34(i,3)*om34(i,3)*rs(4)*cos(th34(i,3)*dr);
F(2)=-a2*rs(2)*cos(c(2))+om2*om2*rs(2)*sin(c(2))+om34(i,2)*om34(i,2)*rs(3)*sin(th34(i,2)*dr) -om34(i,3)*om34(i,3)*rs(4)*sin(th34(i,3)*dr);
a=inv(E)*F';
a3=a(1);
a4=a(2);
G=[rs(6)*sin(th34(i,4)*dr) 1;-rs(6)*cos(th34(i,4)*dr) 0];
H=[-om34(i,3).*om34(i,3)*rs(5)*cos(th34(i,3))-om34(i,4).*om34(i,4)*rs(6)*cos(th34(i,4));-om34 (i,3).*om34(i,3)*rs(5)*sin(th34(i,3))-om34(i,4).*om34(i,4)*rs(6)*sin(th34(i,4))];
b=inv(G)*H;
b1=b(1);
b2=b(2);
a34(i,:)=[t a3 a4 b1 b2];
end
syms p1 p2 p3 p4 jsd1x jsd1y jsd2x jsd2y jsd3x jsd3y jsd4x jsd4y jsd5x jsd5y F2x F2y M2 F3x F3y M3 F5x F5y M5 F6x F1 F2
for i=1:360
t=i*T/360;
p1=i*dth;
om2=400;
a2=0;
jsd2x=-0.08*om2*om2*cos(p1)-0.13*(om34(i,2)*om34(i,2)*cos(th34(i,2)*dr)+a34(i,2)*sin(th34( i,2)*dr));
jsd2y=-0.08*om2*om2*sin(p1)-0.13*(om34(i,2)*om34(i,2)*sin(th34(i,2)*dr)+a34(i,2)*cos(th34(i ,2)*dr));
jsd3x=-0.3*(om34(i,3)*om34(i,3)*cos(th34(i,3)*dr)+a34(i,3)*sin(th34(i,3)*dr));
jsd3y=-0.3*(om34(i,3)*om34(i,3)*sin(th34(i,3)*dr)+a34(i,3)*cos(th34(i,3)*dr));
jsd4x=-0.4*om34(i,3)*om34(i,3)*cos(th34(i,3)*dr)-0.23*(om34(i,4)*om34(i,4)*cos(th34(i,4)*dr) +a34(i,4)*sin(th34(i,4)*dr));
jsd4y=-0.4*om34(i,3)*om34(i,3)*sin(th34(i,3)*dr)-0.23*(om34(i,4)*om34(i,4)*sin(th34(i,4)*dr)+
a34(i,4)*cos(th34(i,4)*dr));
jsd5x=-0.4*om34(i,3)*om34(i,3)*cos(th34(i,3)*dr)-0.46*(om34(i,4)*om34(i,4)*cos(th34(i,4)*dr) +a34(i,4)*sin(th34(i,4)*dr));
F2x=-60/9.8*jsd2x;
F2y=-60/9.8*jsd2y;
M2=-0.08*a34(i,2);
F3x=-72/9.8*jsd3x;
F3y=-72/9.8*jsd3y;
M3=-0.1*a34(i,3);
F5x=-85/9.8*jsd4x;
F5y=-85/9.8*jsd4y;
M5=-0.12*a34(i,4);
F6x=-85/9.8*jsd5x;
F1=4000.0;
X=[0 1 0 1 0 0 0 0 0 0 0 0 0 0 0
1 0 1 0 0 0 0 0 0 0 0 0 0 0 0
0 0 -0.08*sin(p1) 0.08*cos(p1) 0 0 0 0 0 0 0 0 0 0 1
0 0 0 -1 0 1 0 0 0 0 0 0 0 0 0
0 0 -1 0 1 0 0 0 0 0 0 0 0 0 0
0 0 -0.13*sin(th34(i,2)*dr) 0.13*cos(th34(i,2)*dr) -0.13*sin(th34(i,2)*dr) 0.13*sin(th34(i,2)*dr) 0 0 0 0 0 0 0 0 0
0 0 0 0 0 -1 0 1 0 1 0 0 0 0 0
0 0 0 0 -1 0 1 0 1 0 0 0 0 0 0
0 0 0 0 0 0 0.3*sin(th34(i,3)*dr) -0.3*cos(th34(i,3)*dr) -0.1*sin(th34(i,3)*dr) -0.1*cos(th34(i,3)*dr) 0 0 0 0 0
0 0 0 0 0 0 0 0 -1 0 1 0 0 0 0
0 0 0 0 0 0 0 0 0 -1 0 1 0 0 0
0 0 0 0 0 0 0 0 -0.23*sin(th34(i,4)*dr) 0.23*cos(th34(i,4)*dr) 0.23*sin(th34(i,4)*dr) 0.23*cos(th34(i,4)*dr) 0 0 0
0 0 0 0 0 0 0 0 0 0 0 -1 0 1 0
0 0 0 0 0 0 0 0 0 0 0 0 -1 0 0];
Z=[0 0 0 -F2y -F2x -M2 -F3y -F3x -M3 -F5x -F5y -M5 0 F6x+F1 0];
t=X/Z;
t15=t(15)
power(i,:)=[t t15];
end
plot(power(:,1),power(:,2))
基于matlab的GUI设计——机械原理教学演示系统
机械原理教学演示系统——基于matlab的GUI设计 xxx 指导老师： 20年月日
目录 一、功能简介 (3) 二、总界面 (3) 三、凸轮模块 (4) 四、齿轮模块 (6) 五．参考书目 (6) 六．附录（部分程序源代码） (7)
一、功能简介 本系统能实现机械原理教学过程中凸轮模块与齿轮模块的设计与运动仿真，加深对机械原理课程学习的理解。 二、总界面 总界面标题设置：set(gcf,'name','机械原理教学演示系统 made by 翟鲁鑫'); 背景图片设置：ximg=imread('机械原理课本.jpg'); imshow(ximg); 背景声音播放：Fs=44100; [ywav,Fs]=wavread('夜的钢琴曲 - 六3.wav'); sound(ywav,Fs); 到各个模块：到凸轮模块 clc close(gcf); clear all
GUItulun 到齿轮模块 clc close(gcf); clear all GUIchilun 关闭系统：clc question='真的要退出吗？'; title='确认退出？'; button=questdlg(question,title,'是','否','是'); switch button case'是' clear all close case'否' return end 三、凸轮模块 设计要点： 1.背景声音设置方法同主界面
2.推程角、远休角、回程角之和不能大于360度的判别条件；基圆半径、滚子半径、行程不能为0的判别条件 sr0=get(handles.edit2,'string'); r0=str2num(sr0); if isequal(r0,0) errordlg('基圆半径不能为0,请重新输入','出错'); return end srr=get(handles.edit3,'string'); rr=str2num(srr); if isequal(rr,0) errordlg('滚子半径不能为0,请重新输入','出错'); return end sh=get(handles.edit4,'string'); h=str2num(sh); if isequal(h,0) errordlg('行程不能为0,请重新输入','出错'); return end n3=phi01+phi02+phi03; %推程角、远休止角与回程角的总和 if n3>360 errordlg('角度之和大于360，请重新输入','出错'); end 3.仿真程序。采用for 循环以及m(j)=getframe之前要先使用moviein函数Initialize movie frame memory，否则要提示错误 4.仿真之前要先清除绘图时留下的图像，命令如下cla(handles.axes1); 5.回主界面之前要先clear all，关闭音乐、清空global定义的全局变量，以防全局变量影响下一程序 6.图像保存。绘制出的图像可以保存供以后查看。主要命令有uiputfile()、imwrite() [sfilename,sfilepath]=uiputfile({'*.jpg';'*.bmp';'*.tif';'*.*'},... '保存图像文件','unititled.jpg'); if ~isequal([sfilename,sfilepath],[0,0]) sfilefullname=[sfilepath ,sfilename]; h_tulun = getframe(handles.axes1); imwrite(h_tulun.cdata,sfilefullname); else msgbox('您按了取消，保存失败','保存失败','error'); end
哈工大机械设计大作业V带传动设计完美版
哈工大机械设计大作业V带传动设计完美版
————————————————————————————————作者：————————————————————————————————日期： ?
Ｈarbin Insｔiｔutｅof Ｔｅchnoｌogｙ 机械设计大作业说明书 大作业名称：机械设计大作业 设计题目：V带传动设计 班级: 设计者: 学号: 指导教师: 设计时间: 2014.1０.25 哈尔滨工业大学
目录 一、大作业任务书 ........................................................................................................................... 1 二、电动机的选择 ........................................................................................................................... 1 三、确定设计功率d P ..................................................................................................................... 2 四、选择带的型号 ........................................................................................................................... 2 五、确定带轮的基准直径1d d 和2d d ............................................................................................. 2 六、验算带的速度 ........................................................................................................................... 2 七、确定中心距a 和V 带基准长度d L ......................................................................................... 2 八、计算小轮包角 ........................................................................................................................... 3 九、确定V 带根数Z ........................................................................................................................ 3 十、确定初拉力0F ......................................................................................................................... 3 十一、计算作用在轴上的压力 ....................................................................................................... 4 十二、小V 带轮设计 .. (4) 1、带轮材料选择 ............................................................................................................. 4 2、带轮结构形式 . (4) 十二、参考文献 ............................................................................................................................... 6 ?
春MATLAB仿真期末大作业
MATLAB仿真 期末大作业 姓名：班级：学号：指导教师：
2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。
仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。
哈尔滨工业大学机械设计大作业_带传动电算
H a r b i n I n s t i t u t e o f T e c h n o l o g y 上机电算说明书 课程名称：机械设计 电算题目：普通V带传动 院系：机电工程学院 班级： 设计者： 学号： 指导教师： 设计时间：2015.11.11-2015.12.1 哈尔滨工业大学
目录 一、普通V带传动的内容 (1) 二、变量标识符 (1) 三、程序框图 (2) 四、V带设计C程序 (3) 五、程序运行截图 (10) 参考文献 (11)
一、普通V带传动的内容 给定原始数据：传递的功率P，小带轮转速n1 传动比i及工作条件 设计内容：带型号，基准长度Ld，根数Z，传动中心距a，带轮基准直径dd1、dd2，带轮轮缘宽度B，初拉力F0和压轴力Q。 二、变量标识符 为了使程序具有较好的可读性易用性，应采用统一的变量标识符，如表1所示。表1变量标识符表。 表1 变量标识符表
三、程序框图
四、V带设计c程序 #include
MATLAB程序：已知三个位置设计平面四杆机构求解程序(位移矩阵法)
%MATLAB程序：已知三个位置设计平面四杆机构求解程序（位移矩阵法） clear;clc; %凡是变量名前带v的为数值变量，不带的是符号变量 vxp1=0; vyp1=0; vsita1=0*pi/180; vxp2=-2; vyp2=6; vsita2=40*pi/180; vxp3=-10; vyp3=8; vsita3=90*pi/180; %精确位置P1，P2，P3及各角度 vsita12=vsita2-vsita1; vsita13=vsita3-vsita1; vxa=-10; vya=-2; vxd=-5; vyd=-2; %选定A,D点 %所有数值均在此确定，更改此处即可解出不同数值的四杆机构位移矩阵方程 syms xp1 yp1 xp2 yp2 xp3 yp3 sita12 sita13; syms xa ya xb1 yb1 xb2 yb2 xb3 yb3; f1='(xb2-xa)^2+(yb2-ya)^2=(xb1-xa)^2+(yb1-ya)^2'; f2='(xb3-xa)^2+(yb3-ya)^2=(xb1-xa)^2+(yb1-ya)^2'; %前两个机构方程 f3='xb2=cos(sita12)*xb1-sin(sita12)*yb1+xp2-xp1*cos(sita12)+yp1*sin(sita12)'; f4='yb2=sin(sita12)*xb1+cos(sita12)*yb1+yp2-xp1*sin(sita12)-yp1*cos(sita12)'; %由第一个位移矩阵方程得出 f5='xb3=cos(sita13)*xb1-sin(sita13)*yb1+xp3-xp1*cos(sita13)+yp1*sin(sita13)'; f6='yb3=sin(sita13)*xb1+cos(sita13)*yb1+yp3-xp1*sin(sita13)-yp1*cos(sita13)'; %由第二个位移矩阵方程得出 f1=subs(f1,{xa,ya},{vxa,vya}); f2=subs(f2,{xa,ya},{vxa,vya}); f3=subs(f3,{xp1,xp2,yp1,sita12},{vxp1,vxp2,vyp1,vsita12}); f4=subs(f4,{xp1,yp1,yp2,sita12},{vxp1,vyp1,vyp2,vsita12}); f5=subs(f5,{xp1,xp3,yp1,sita13},{vxp1,vxp3,vyp1,vsita13}); f6=subs(f6,{xp1,yp1,yp3,sita13},{vxp1,vyp1,vyp3,vsita13}); %代入具体数值 [xb1,xb2,xb3,yb1,yb2,yb3]=solve(f1,f2,f3,f4,f5,f6); %解方程 vxb1=vpa(xb1); vyb1=vpa(yb1); vxb2=vpa(xb2); vyb2=vpa(yb2); vxb3=vpa(xb3); vyb3=vpa(yb3); (vxb1-vxa)^2+(vyb1-vya)^2; (vxb2-vxa)^2+(vyb2-vya)^2; (vxb3-vxa)^2+(vyb3-vya)^2; %去掉这三行分号可验证B点三个位置是否距离A点相等 syms xd yd xc1 yc1 xc2 yc2 xc3 yc3;
电机学matlab仿真大作业报告
. 基于MATLAB的电机学计算机辅助分析与仿真 实验报告
一、实验内容及目的 1.1 单相变压器的效率和外特性曲线 1.1.1 实验内容 一台单相变压器，N S =2000kVA, kV kV U U N N 11/127/21=，50Hz ，变压器的参数 和损耗为008.0* ) 75(=C k o R ，0725.0*=k X ，kW P 470=，kW P C KN o 160)75(=。 (1)求此变压器带上额定负载、)(8.0cos 2滞后=?时的额定电压调整率和额定效率。 (2)分别求出当0.1,8.0,6.0,4.0,2.0cos 2=?时变压器的效率曲线，并确定最大效率和达到负载效率时的负载电流。 (3)分析不同性质的负载（),(8.0cos 0.1cos ),(8.0cos 222超前，滞后===???）对变压器输出特性的影响。 1.1.2 实验目的 （1）计算此变压器在已知负载下的额定电压调整率和额定效率 （2）了解变压器效率曲线的变化规律 （3）了解负载功率因数对效率曲线的影响 （4）了解变压器电压变化率的变化规律 （5）了解负载性质对电压变化率特性的影响 1.1.3 实验用到的基本知识和理论 （1）标幺值、效率区间、空载损耗、短路损耗等概念 （2）效率和效率特性的知识 （3）电压调整率的相关知识 1.2串励直流电动机的运行特性 1.2.1实验内容 一台16kw 、220V 的串励直流电动机，串励绕组电阻为0.12Ω，电枢总电阻为0.2Ω。电动势常数为.电机的磁化曲线近似的为直线。其中为比例常数。假设电枢电流85A 时，磁路饱和(为比较不同饱和电流对应的效果，饱和电流可以自己改变）。
2013年机械设计大作业轴设计
大作业设计说明书 课程名称： 机 械 设 计 设计题目： 设计搅拌机用单级斜齿圆柱 齿轮减速器中的低速轴 院 系： 理 学 院 专业班级： 机械电子工程0211411班 设 计 者： 学 号： 设计时间： 2013年12月20日 湖 北 民 族 学 院 HUBEI MINZU UNIVERSITY
目录(宋体，三号，加粗，居中) 1、设计任务书 (1) 2、…………………………………………………………… 3、轴结构设计………………………………………………… 3.1轴向固定方式……………………………………………………… 3.2选择滚动轴承类型……………………………………………………… 3.3键连接设计………………………………………………… 3.4阶梯轴各部分直径确定…………………………………………………… 3.5阶梯轴各部段长度及跨距的确定……………………………………… 4、轴的受力分析…………………………………………………………… 4.1画轴的受力简图……………………………………………………… 4.2计算支反力……………………………………………………… 4.3画弯矩图……………………………………………………… 4.4画扭矩图……………………………………………………… 5、校核轴的弯扭合成强度…………………………………………………… 6、轴的安全系数校核计算……………………………………………… 7、参考文献…………………………………………… 注：其余小四，宋体。自己按照所需标题编号，排整齐。
设计任务书 1.已知条件 某搅拌机用单级斜齿圆柱减速器简图如上所示。已知：电动机额定功率P=4kW，转速n1=750r/min，低速轴转速n2=130r/min，大齿轮节圆直径d2=300mm，宽度B2=90mm，轮齿螺旋角β=12°，法向压力角αn=20°。 2．设计任务 设计搅拌机用单级斜齿圆柱减速器中的高速级/低速轴（包括选择两端的轴承及外伸端的联轴器）。 要求：（1）完成轴的全部结构设计； （2）根据弯扭合成理论验算轴的强度； （3）精确校核轴的危险截面是否安全。 - 1 -
matlab 大作业
上海电力学院 通信原理Matlab仿真 实验报告 实验名称： 8QAM误码率仿真 试验日期： 2014年 6月3日 专业：通信工程 姓名：罗侃鸣 班级： 2011112班 学号： 20112272
一、实验要求 写MATLAB程序，对图示的信号星座图完成M=8的QAM通信系统Monte Carlo仿真，在不同SNRindB=0：15时，对N=10000（3比特）个符号进行仿真。画出该QAM系统的符号误码率。 二、实验原理 1 QAM调制原理 QAM(Quadrature Amplitude Modulation)正交幅度调制技术，是用两路独立的基带信号对两个相互正交的同频载波进行抑制载波双边带调幅，利用这种已调信号的频谱在同一带宽内的正交性，实现两路并行的数字信息的传输。该调制方式通常有8QAM，16QAM，64QAM。 QAM调制实际上就是幅度调制和相位调制的组合，相位+ 幅度状态定义了一个数字或数字的组合。QAM的优点是具有更大的符号率，从而可获得更高的系统效率。通常由符号率确定占用带宽。因此每个符号的比特（基本信息单位）越多，频带效率就越高。 调制时，将输入信息分成两部分：一部分进行幅度调制；另一部分进行相位调制。对于星型8QAM信号，每个码元由3个比特组成，可将它分成第一个比特和后两个个比特两部分。前者用于改变信号矢量的振幅，后者用于差分相位调制，通过格雷编码来改变当前码元信号矢量相位与前一码元信号矢量相位之间的相位差。 QAM是一种高效的线性调制方式，常用的是8QAM，16QAM，64QAM等。当随着M 的增大，相应的误码率增高，抗干扰性能下降。 2 QAM星座图 QAM调制技术对应的空间信号矢量端点分布图称为星座图。QAM的星座图呈现星状分
简单的机械设计作品【机械设计大作业】
简单的机械设计作品【机械设计大作业】 《机械设计》齿轮设计程序 #include #include floatmin(floatx,floaty); floatmax(floatx,floaty); voidmain() {intB1,B2,b,z1,z2; floatP1,u,T1,Kt,FAId,CHlim1,CHlim2,S,n1,j,l,ZE,N1,N2,KHN 1,KHN2, d1t,v,mt,h,p,KA,KV,KH1,KH2,KF1,KF2,YFa1,YFa2,YSa1,YSa2,C FE1,CFE2, KFN1,KFN2,CH1,CH2,CH,CF1,CF2,a,d1,d2,m,w,K,bt,z1t,z2t,q; printf("请输入小齿轮的齿数z1和齿数比u:\n"); scanf("%f%f",&z1t,&u); printf("下面进行齿面强度计算\n"); printf("请输入输入功率P1,小齿轮的转速n1,j,齿轮工作时间l:\n"); scanf("%f%f%f%f",&P1,&n1,&j,&l); T1=(9550000*P1)/n1; N1=60*n1*j*l; N2=N1/u;
printf("T1=%10.4eN1=%10.4eN2=%10.4e\n",T1,N1,N2); printf("请输入接触疲劳寿命系数KHN1,KHN2:\n"); scanf("%f%f",&KHN1,&KHN2); CH1=KHN1*CHlim1/S; CH2=KHN2*CHlim2/S; CH=min(CH1,CH2); printf("CH1=%10.4f\nCH2=%10.4f\nCH=%10.4f\n",CH1,CH2,CH); printf("\n"); printf("请输入载荷系数Kt,齿宽系数FAId,弹性影响系数 ZE:\n"); scanf("%f%f%f",&Kt,&FAId,&ZE); q=pow(ZE/CH,2)*Kt*T1*(u+1)/(FAId*u); d1t=2.32*pow(q,1.0/3.0); v=3.1415926*d1t*n1/60000; bt=FAId*d1t; mt=d1t/z1t; h=2.25*mt; p=bt/h; printf("d1t=%10.4fv=%10.4fmt=%10.4fp=%10.4f\n",d1t,v,mt, p); printf("\n"); printf("根据v,P和精度等级查KV,KH1,KF1,KH2,KF2,KA:\n"); scanf("%f%f%f%f%f%f",&KV,&KH1,&KF1,&KH2,&KF2,&KA);
matlab机电系统仿真大作业
一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环
3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置
M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置
机械原理课程设计Matlab编程
/*Matlab程序*/ l1 = 59.1000; l2 = 263.9000; l3=120; l4=266.83; l5=180; l6=45; x2=170; y2=132.7289; w1=9.4248; N=42:10:402; ay=119:10:479 a=2*l1*l3*sin(N/180*pi); b=2*l3*(l1*cos(N/180*pi)-l4); c=l2^2-l1^2-l3^2-l4^2+2*l1*l4*cos(N/180*pi); jiao3=2*atan((a- sqrt(a.^2+b.^2-c.^2))./ (b-c))/pi*180+77 g=2*l1*l2*sin(N/180*pi); h=2*l2*(l1*cos(N/180*pi)-l4); m=l1^2+l2^2+l4^2-l3^2-2*l1*l4*cos(N/180*pi); jiao2=2*atan((g- sqrt(g.^2+h.^2-m.^2))./ (h-m))/pi*180+77 lof=-sqrt(l6^2-x2^2-l5^2+2*l5*x2*cos((180+jiao3)/180*pi)+l5^2*sin(j iao3/180*pi).^2)+y2-l5* sin((180+jiao3)/180*pi)
j12=N-(jiao2-77); j32=jiao3-jiao2; j13=(N-(jiao3-77)); j23=(jiao2-jiao3); w3=(w1*l1*sin(j12/180*pi))./ (l3*sin(j32/180*pi)) w2=(-1*w1*l1*sin(j13/180*pi))./(l2*sin(j23/180*pi)) a3=(w1^2*l1*cos(j12/180*pi)+w2.^2*l2-(w3.^2).*(l3*cos(j32/180*pi)) )./ (l3*sin(j32/180*pi)) a2=(-w1^2*l1*cos(j13/180*pi)-(w2.^2).*(l2*cos(j23/180*pi))+l3*w3.^ 2)./ (l3*sin(j23/180*pi)) jiao4=acos((x2-l5*cos((180+jiao3)/180*pi))/l6)/pi*180 w4=((-l5*sin((pi+jiao3)/180*pi)).*w3)./ (l6*sin(jiao4/180*pi)) vof=((l5*sin((180+jiao3-jiao4)/180*pi)).*w3).* sin(jiao4/180*pi) aof=(l6*w4.^2+(l5*w3.^2).*(cos((180+jiao3-jiao4)/180*pi))+l5*a3.*si n((180+jiao3-jiao4)/180*pi))./sin(jiao4/180*pi) 作图程序： /*F点的位移*/ plot(N+77,lof,'-xk') xlabel('AB杆的角度'),ylabel('F点的位移/(mm)') title('F点的位移曲线图') text(100, 171.3339,'初始值= 171.3339')
电机大作业(MATLAB仿真-电机特性曲线)
电机大作业 专业班级：电气XXXX 姓名：XXX 学号：XXX 指导老师：张威
一、研究课题（来源：教材习题 4-18 ） 1. 74 、R 2 0.416 、X 2 3.03 、R m 6. 2 X m 75 。电动机的机械损耗p 139W,额定负载时杂散损耗p 320W, 试求额定负载时的转差率、定子电流、定子功率因数、电磁转矩、输出转矩和效 率。 二、编程仿真 根据T 形等效电路: 3D - R Q 运用MATLAB 进行绘图。MATLAB 文本中，P N PN ，U N UN ，尺 R 1， X 1 X1 , R 2 R 2，X 2 X 2，R m Rm, X m Xm ，p pjixiesunh ao ， p pzasansunhao 。定子电流I11，定子功率因数 Cosangle1，电磁转矩Te ， 效率 Xiaolv 。 1.工作特性曲线绘制 MATLA 文本： R1=0.715;X 仁1.74;Rm=6.2;Xm=75;R2=0.416;X2=3.03;pjixiesu nhao=139; pzasa nsu nhao=320;p=2;m 仁 3; ns=1500;PN=17000;UN=380;fN=50; Z1=R1+j*X1; Zm=Rm+j*Xm; for i=1:2500 s=i/2500; nO=n s*(1-s); Z2=R2/s+j*X2; Z=Z1+Zm*Z2/(Zm+Z2); 有一台三相四极的笼形感应电动机， 参数为P N 17kW 、U N 380V (△联 Rm 结）、尺 0. 715 、X j lcr S
U1=UN; I1=U1/Z; l110=abs(l1); An gle 仁an gle(ll); Cosa ngle10=cos(A ngle1); P仁3*U1*l110*Cosa ngle10; l2=l1*Zm/(Zm+Z2); Pjixie=m1*(abs(I2))A2*(1-s)/s*R2; V=(1-s)*pi*fN; Te0=Pjixie/V; P20=Pjixie-pjixies un hao-pzasa nsun hao; Xiaolv0=P20/P1; P2(i)=P20; n (i)=n0; l11(i)=l110; Cosa ngle1(i)=Cosa ngle10; Te(i)=Te0; Xiaolv(i)=Xiaolv0; hold on; end figure(1) plot(P2, n); xlabel('P2[W]');ylabel(' n[rpm]'); figure(2) plot(P2,l11); xlabel('P2[W]');ylabel('l1[A]'); figure(3) plot(P2,Cosa nglel); xlabel('P2[W]');ylabel('go nglvyi nshu'); figure(4) plot(P2,Te); xlabel('P2[W]');ylabel('Te[Nm]'); figure(5) plot(P2,Xiaolv); xlabel('P2[W]');ylabel('xiaolv');
运动控制MATLAB仿真
大作业: 直流双闭环调速MATLAB仿真 运动控制技术课程名称： 名：姓电气学院院：学 自动化业：专 号：学 孟濬指导教师： 2012年6月2日
------------------------------------- -------------学浙大江 李超 一、Matlab仿真截图及模块功能描述 Matlab仿真截图如下，使用Matlab自带的直流电机模型： 模块功能描述： ⑴电机模块（Discrete DC_Machine）：模拟直流电机 ⑵负载转矩给定（Load Torque）：为直流电机添加负载转矩 ⑶Demux：将向量信号分离出输出信号 ⑷转速给定（Speed Reference）：给定转速 ⑸转速PI调节（Speed Controller）：转速PI调节器，对输入给定信号与实际信号
的差值进行比例和积分运算，得到的输出值作为电流给定信号。改变比例和积分运算系数可以得到不同的PI控制效果。 ⑹电流采样环节（1/z）：对电流进行采样，并保持一个采样周期 ⑺电流滞环调节（Current Controller）：规定一个滞环宽度，将电流采样值与给定值进行对比，若：采样值＞给定值+0.5*滞环宽度，则输出0； 若：采样值＜给定值—0.5*滞环宽度，则输出1； 若：给定值—0.5*滞环宽度＜采样值＜给定值+0.5*滞环宽度，则输出不变 输出值作为移相电压输入晶闸管斩波器控制晶闸管触发角 ：根据输入电压改变晶闸管触发角，从而改变电机端电压。GTO⑻晶闸管斩波．⑼续流二极管D1：在晶闸管关断时为电机续流。 ⑽电压传感器Vd：测量电机端电压 ⑾示波器scope：观察电压、电流、转速波形 系统功能概括如下：直流电源通过带GTO的斩波器对直流电机进行供电，输出量电枢电流ia和转速wm通过电流环和转速环对GTO的通断进行控制，从而达到对整个电机较为精确的控制。 下面对各个部分的功能加以详细说明： （1）直流电机 双击电动机模块，察看其参数：
哈工大机械设计大作业之千斤顶说明书
工业大学 机械设计作业设计计算说明书 题目螺旋起重器（千斤顶） 系别机电工程学院 班号 0908103 孟子航 学号 1090810314 日期 2011年9月13日
工业大学 机械设计作业任务书 题目 螺旋起重器（千斤顶） 设计原始数据： 表3.1 螺旋起重器的示意图及已知数据 题号 起重器/Q F kN 最大起重高度/H mm 3.1.1 30 180 3.1.2 40 200 3.1.3 50 150 设计要求： （1） 绘制装配图一，画出起重器的全比结构，按照比例装配图要求标注尺寸、序号及填 写明细栏、标题栏、编写技术要求。 （2） 撰写设计说明书一份，主要包括起重器各部分尺寸的计算，对螺杆和螺母螺纹牙强 度、螺纹副自锁性、螺杆的稳定性的校核等。
目录 一、设计题目---------------------------------------------------------------2 二、螺母、螺杆选材---------------------------------------------------------2 三、螺杆、螺母设计计算 3.1 耐磨性计算 ------------------------------------------------------------------------------------------------2 3.2 螺杆强度校核 ---------------------------------------------------------------------------------------------3
机械原理课程设计matlab程序及成果图
Wjr_main.m %1.输入已知数据 clear; l2=0.1605;%AB的长度单位m l4=0.6914;%CD的长度单位m l5=0.2074;%DE的长度单位m l1=0.370;%AC的长度单位m l1p=0.6572;%CF的长度单位m omg2=8.378; af2=0; hd=pi/180; du=180/pi; %2.调用子函数abc.m计算牛头刨机构位移，角速度，角加速度for n1=1:689; tt2(n1)=-0.4488+(n1-1)*hd; ll=[l2,l4,l5,l1,l1p]; [tt,omg,af]=abc(tt2(n1),omg2,af2,ll); s4(n1)=tt(1); tt4(n1)=tt(2); tt5(n1)=tt(3); sE(n1)=tt(4);
v34(n1)=omg(1); omg4(n1)=omg(2); omg5(n1)=omg(3); vE(n1)=omg(4); a3(n1)=af(1); af4(n1)=af(2); af5(n1)=af(3); aE(n1)=af(4); end %3.位移，角速度，角加速度 figure(1); n1=1:689; t=(n1-1)*pi/180; subplot(2,2,1); %绘角位移及位移线图plot(t,tt4*du,'r-.'); grid on; hold on; axis auto; [haxes,hline1,hine2]=plotyy(t,tt5*du,t,sE);
grid on; hold on; xlabel('时间/份'); axes(haxes(1)); ylabel('角位移/\circ'); axes(haxes(2)); ylabel('位移/m'); hold on; grid on; text(1.15,-0.65,'tt_4'); text(3.4,0.27,'tt_5'); text(2.25,-0.15,'s_E'); subplot(2,2,2); %绘角速度及速度线图plot(t,omg4,'r-.'); grid on; hold on; axis auto; [haxes,hline1,hline2]=plotyy(t,omg5,t,vE); grid on; hold on; xlabel('时间/份') axes(haxes(1));
华科matlab大作业
MATLAB语言、控制系统分析与设计 大作业 升降压斩波电路MATLAB仿真 专业：电气工程及其自动化 班级： 设计者： 学号： 评分： 华中科技大学电气与电子工程学院 2016 年1月
评分栏
目录 一、引言 (4) 二、电路设计与仿真 (4) 三、仿真结果及分析 (7) 四、深入讨论 (10) 五、总结 (10) 六、参考资料 (11)
升降压斩波电路MATLAB 仿真 一．引言 Buck/Boost 变换器是输出电压可低于或高于输入电压的一种单管直流变换器，其主电路与Buck 或Boost 变换器所用元器件相同，也有开关管、二极管、电感和电容构成。与Buck 和Boost 电路不同的是，电感L 在中间，不在输出端也不在输入端，且输出电压极性与输入电压相反。开关管也采用PWM 控制方式。Buck/Boost 变换器也有电感电流连续喝断续两种工作方式，在此只讨论电感电流在连续状态下的工作模式。 二．电路设计与仿真 1、电路原理 当可控开关V 处于通态时，电源E 经V 向电感L 供电使其储存能量，此时电流为i1。同时，电容C 维持输出电压基本恒定并向负载R 供电。此后，使V 关断，电感L 中储存的能量向负载L 释放，电流为i2。负载电压极性为上正下负，与电源电压极性相反。 稳态时，一个周期T 内电感L 两端电压UL 对时间积分为零，即 当V 处于通态期间，UL=E ；而当V 处于断态期间，UL=-Uo 。于是 所以，输出电压为 =?dt T L U off o on t U Et =E E t T t E t t U on on off on o α α -=-== 1