第五章Matlab应用之机械零件设计.pptVIP

（一）机械零件设计： 向量、受力和刚体平衡 梁的挠度计算 四连杆机构运动设计与分析 凸轮轮廓 （二）动力学与振动： 轨迹 单自由度系统 两自由度系统 （三）优化设计（重点）： 线性规划 无约束优化 单目标约束优化 约束优化问题 （四）工程统计 教学目标 介绍运用Matlab进行各种类型机械零件设计的分析方法 目录 5.1 向量、受力和刚体平衡 5.2 梁的挠度计算 5.3 连杆机构运动设计与分析 5.4 凸轮轮廓 习题 5.1 向量、受力和刚体平衡 5.1.1 向量 5.1.2 力的分解与合成 5.1.3 刚体的平衡 5.1 向量、受力和刚体平衡 5.1.1 向量 5.1 向量、受力和刚体平衡 5.1 向量、受力和刚体平衡 5.1.2 力的合成及分解 5.1 向量、受力和刚体平衡 5.1 向量、受力和刚体平衡 5.1 向量、受力和刚体平衡 5.1.3 刚体的平衡 5.1 向量、受力和刚体平衡 如下图所示杆系，设已知：G1=200; G2=100; L1= 2; L2=sqrt(2);theta1 =30*pi/180; theta2 = 45*pi/180;求其支撑反力Na,Nb,Nc。 5.1 向量、受力和刚体平衡 ΣX=0 Nax + Ncx = 0 ΣY=0 Nay + Ncy - G1 = 0; ΣM=0 Ncy*L1*cos(theta1)-Ncx*L1*sin(theta1)- G1*L1/2*cos(theta1)=0; ΣX=0 Nbx - Ncx = 0; ΣY=0 Nby - Ncy - G2 = 0 ΣM=0 Ncy*L2*cos(theta2)+ Ncx*L2*sin(theta2)+ G2*L2/2*cos(theta2)=0; 这是一组包含六个未知数Nax, Nay, Nbx, Nby, Ncx, Ncy的六个线性代数方程,通常是要寻找简化的方法,但利用MATLAB工具，就可以列出矩阵方程AX=B,(其中X=[Nax,Nay,Nbx,Nby,Ncx,Ncy]T，可用矩阵除法直接来解。 5.1 向量、受力和刚体平衡 程序如下： %给原始参数赋值 G1=200; G2=100; L1= 2; L2 = sqrt (2) ; % 将度化为弧度 theta1 =30*pi/180; theta2 = 45*pi/180; % 则按此次序，系数矩阵A,B可写成下式 A=[1,0,0,0,1,0; 0,1,0,0,0,1; 0,0,0,0,-L1*sin(theta1), L1* cos(theta1); 0,0,1,0,-1,0; 0,0,0,1,0,-1; 0,0,0,0, L2* sin(theta2), L2* cos(theta2)] B=[0;G1; G1*L1/2*cos(theta1);0;G2;-G2*L2/2*cos(theta2)] X = A\B; %用左除求解线性方程组 程序（compute_naodu.m）如下 L=2; P=2000; L1=1.5; %给出常数 E = 200e9; I=2e-5; x=linspace(0,L,101); dx=L/100; %将L分100段 n1=L1/dx+1; % 确定x=L1处对应的下标 M1=-P*( L1-x(1:n1)); % 第一段弯矩赋值 M2=zeros(1,101-n1); % 第二段弯矩赋值(全为零) M=[M1,M2]; % 全梁的弯矩 A=cumsum(M)*dx/(E*I); % 对弯矩积分求转角 Y=cumsum(A)*dx; % 对转角积分求挠度 subplot(3,1,1),plot(x,M),grid,title(弯矩图) % 绘弯矩图 subplot(3,1,2),plot(x,A), grid,title(转角图) % 绘转角图 subplot(3,1,3),plot(x,Y),grid,title(挠度图) % 绘挠度图 所得的结果见右，注意几根曲线之