用matlab设计一平底直动从动件盘形凸轮机构凸轮的基圆半径𝑟0 =40mm从动件升程ℎ =10mm∅ =120°∅𝑠𝑠=30°∅𝑠𝑠′ =90°∅′=120°从动件在推程和回程均作简谐运动余弦加速度运动。请给出代码

由于需要进行余弦加速度运动，我们可以利用matlab中的sinc函数来实现。以下是设计这一平底直动从动件盘形凸轮机构的matlab代码：

clc clear all close all

r0 = 40; % 基圆半径 h = 10; % 从动件升程 phi = 120; % 凸轮旋转角度 phi_ss = 30; % 凸轮起始角度 phi_ss_prime = 90; % 凸轮起始角度的导数 phi_prime = 120; % 凸轮旋转角度的导数

% 计算从动件的升降速度 T = 2*pi; % 凸轮一周的弧长 v_max = h/(T/4); % 从动件最大速度 a_max = v_max/(T/4); % 从动件最大加速度

% 计算余弦加速度 t = linspace(0,1,1000); u = 2pit; s = h/2 * (1 - cos(pit)); v = h/2 * pisin(pit); a = h/2 * pi^2cos(pi*t);

% 计算从动件的位置 position = zeros(size(u)); for i=1:length(u) position(i) = r0 + s(i); end

% 绘制凸轮的轮廓 theta = linspace(0,2pi,100); camber = zeros(size(theta)); for i=1:length(theta) if (theta(i)>=deg2rad(phi_ss) && theta(i)<=deg2rad(phi_ss+phi)) camber(i) = r0 + h/2 + a_max/h * (r0 + h/2 - r0cos(theta(i)-deg2rad(phi_ss))); elseif (theta(i)>=deg2rad(phi_ss+phi) && theta(i)<=deg2rad(phi_prime+phi_ss)) camber(i) = r0 + h/2 + v_max/h * (theta(i)-deg2rad(phi_ss+phi)) + a_max/h * (r0 + h/2 - r0cos(theta(i)-deg2rad(phi_ss))); elseif (theta(i)>=deg2rad(phi_prime+phi_ss) && theta(i)<=2pi) camber(i) = r0 + h/2 + v_max/h * (2pi-phi-phi_ss-theta(i)) + a_max/h * (r0 + h/2 - r0cos(theta(i)-deg2rad(phi_ss))); end end

% 绘制从动件的运动轨迹 figure; hold on; plot(camber.*cos(theta),camber.*sin(theta),'b'); plot(position.*cos(u),position.*sin(u),'r'); axis equal; legend('凸轮轮廓','从动件轨迹'); xlabel('x'); ylabel('y');