基于matlab的多类结构动力学，凸轮机构、双凸轮、弦振动模拟、阻尼振动 、四连杆机构 、套杆运动 、三根弹簧作用的振子。 程序已调通，可直接运行。
最近在Matlab里玩了不少有趣的结构动力学模拟，感觉像是打开了一个神奇的物理世界大门😃。今天就来跟大家分享一下基于Matlab实现的多类结构动力学模拟，这里面有凸轮机构、双凸轮、弦振动模拟、阻尼振动 、四连杆机构 、套杆运动 、三根弹簧作用的振子，程序都已经调通，可以直接运行哦👏。
凸轮机构模拟
先来说说凸轮机构吧。凸轮机构在很多机械装置中都有重要应用，通过Matlab可以很直观地模拟它的运动。
% 定义凸轮的基本参数 r = 0.1; % 凸轮半径 theta = 0:0.01:2*pi; % 角度范围 x = r*cos(theta); % 凸轮的x坐标 y = r*sin(theta); % 凸轮的y坐标 % 绘制凸轮 plot(x,y,'b','LineWidth',2); xlabel('X (m)'); ylabel('Y (m)'); title('Cam Mechanism'); grid on;代码分析：这里首先定义了凸轮的半径r，然后通过theta的取值范围生成了凸轮在不同角度下的坐标x和y。最后用plot函数绘制出凸轮的形状，设置了线条颜色为蓝色，线宽为2，还添加了坐标轴标签和标题，打开了网格，让图形看起来更直观。
双凸轮模拟
双凸轮的模拟稍微复杂一些，但也很有意思。
% 双凸轮参数 r1 = 0.1; r2 = 0.08; theta1 = 0:0.01:2*pi; theta2 = 0:0.01:2*pi; x1 = r1*cos(theta1); y1 = r1*sin(theta1); x2 = r2*cos(theta2+pi/2); y2 = r2*sin(theta2+pi/2); % 绘制双凸轮 figure; plot(x1,y1,'b','LineWidth',2); hold on; plot(x2,y2,'r','LineWidth',2); xlabel('X (m)'); ylabel('Y (m)'); title('Double Cam Mechanism'); grid on; legend('Cam 1', 'Cam 2'); hold off;代码分析：这里分别定义了两个凸轮的半径r1和r2，以及角度范围theta1和theta2。通过三角函数计算出两个凸轮在不同角度下的坐标x1、y1、x2和y2。然后使用plot函数分别绘制两个凸轮，并且用hold on保持图形，添加了第二个凸轮后，设置了坐标轴标签、标题、网格和图例，最后用hold off关闭保持状态。
弦振动模拟
弦振动是一个经典的物理现象，Matlab也能很好地模拟出来。
% 弦振动模拟参数 L = 1; % 弦长 c = 1; % 波速 T = 2; % 总时间 dt = 0.01; % 时间步长 dx = 0.01; % 空间步长 Nx = round(L/dx); % x方向节点数 Nt = round(T/dt); % 时间节点数 u = zeros(Nx,Nt); % 初始化位移矩阵 % 初始条件 u(:,1) = 0.1*sin(pi*[1:Nx]'/Nx); % 边界条件 u(1,:) = 0; u(Nx,:) = 0; % 迭代求解 for n = 1:Nt-1 for i = 2:Nx-1 u(i,n+1) = 2*u(i,n) - u(i,n-1) + c^2*dt^2/dx^2*(u(i+1,n)-2*u(i,n)+u(i-1,n)); end end % 绘制弦振动 figure; for n = 1:1:Nt plot([0:dx:L],u(:,n)); xlabel('Position (m)'); ylabel('Displacement (m)'); title(['String Vibration at t = ',num2str((n-1)*dt)]); drawnow; pause(0.05); end代码分析：这段代码首先定义了弦振动模拟的各种参数，比如弦长L、波速c、总时间T、时间步长dt、空间步长dx等。然后初始化了位移矩阵u，设置了初始条件和边界条件。通过两层循环，使用有限差分法迭代求解弦振动的位移。最后通过循环绘制出不同时刻弦的振动形状，每次绘制后用drawnow更新图形，并暂停0.05秒，让我们能清楚看到弦振动的动态过程。
阻尼振动模拟
阻尼振动也是很常见的物理现象。
% 阻尼振动模拟参数 m = 1; % 质量 k = 1; % 弹簧常数 c = 0.1; % 阻尼系数 t = 0:0.01:10; % 时间范围 % 初始条件 x0 = 1; v0 = 0; % 求解阻尼振动方程 [~,x] = ode45(@(t,x)[x(2);-(k/m)*x(1)-(c/m)*x(2)],t,[x0;v0]); % 绘制阻尼振动 figure; plot(t,x(:,1)); xlabel('Time (s)'); ylabel('Displacement (m)'); title('Damped Vibration'); grid on;代码分析：这里定义了质量m、弹簧常数k、阻尼系数c和时间范围t。通过ode45函数求解阻尼振动方程，初始条件为初始位移x0和初始速度v0。最后用plot函数绘制出位移随时间的变化曲线，添加了坐标轴标签、标题和网格。
四连杆机构模拟
四连杆机构在机械设计中应用广泛。
% 四连杆机构参数 L1 = 1; L2 = 1; L3 = 1; L4 = 1; theta1 = 0:0.01:2*pi; theta2 = acos((L1^2+L2^2-L3^2-L4^2)/(2*(L1*L2+L3*L4)))+atan2(2*L1*L2*sin(theta1),(L1^2+L2^2-L3^2-L4^2)); x1 = L1*cos(theta1); y1 = L1*sin(theta1); x2 = L2*cos(theta2); y2 = L2*sin(theta2); x3 = L3*cos(theta2); y3 = L3*sin(theta2); x4 = L4*cos(theta1); y4 = L4*sin(theta1); % 绘制四连杆机构 figure; plot([0 x1],[0 y1],'b','LineWidth',2); hold on; plot([x1 x2],[y1 y2],'r','LineWidth',2); plot([x2 x3],[y2 y3],'g','LineWidth',2); plot([x3 x4],[y3 y4],'m','LineWidth',2); plot([x4 0],[y4 0],'k','LineWidth',2); xlabel('X (m)'); ylabel('Y (m)'); title('Four - Bar Linkage Mechanism'); grid on; legend('Link 1', 'Link 2', 'Link 3', 'Link 4'); hold off;代码分析：这段代码定义了四连杆机构各杆的长度L1、L2、L3、L4，通过一系列三角函数计算出各杆端点在不同角度下的坐标x1、y1、x2、y2、x3、y3、x4、y4。然后用plot函数依次绘制出四连杆机构的四根杆，设置了线条颜色和线宽，添加了坐标轴标签、标题、网格和图例。
套杆运动模拟
套杆运动模拟也很有趣。
% 套杆运动模拟参数 L = 1; % 杆长 omega = 1; % 角速度 t = 0:0.01:10; % 时间范围 theta = omega*t; x = L*cos(theta); y = L*sin(theta); % 绘制套杆运动 figure; plot(x,y,'b','LineWidth',2); xlabel('X (m)'); ylabel('Y (m)'); title('Slider Movement'); grid on;代码分析：这里定义了杆长L、角速度omega和时间范围t。通过计算得到角度theta，进而算出杆端点的坐标x和y。最后用plot函数绘制出套杆的运动轨迹，设置了线条颜色和线宽，添加了坐标轴标签、标题和网格。
三根弹簧作用的振子模拟
最后看看三根弹簧作用的振子模拟。
% 三根弹簧作用的振子模拟参数 m = 1; % 振子质量 k1 = 1; k2 = 1; k3 = 1; t = 0:0.01:10; % 时间范围 % 初始条件 x0 = 1; v0 = 0; % 求解运动方程 [~,x] = ode45(@(t,x)[x(2);-(k1 + k2 + k3)/m*x(1)],t,[x0;v0]); % 绘制位移随时间变化 figure; plot(t,x(:,1)); xlabel('Time (s)'); ylabel('Displacement (m)'); title('Oscillator with Three Springs'); grid on;代码分析：此代码定义了振子质量m和三根弹簧的弹簧常数k1、k2、k3以及时间范围t。通过ode45函数求解振子在三根弹簧作用下的运动方程，初始条件为初始位移x0和初始速度v0。最后用plot函数绘制出振子位移随时间的变化曲线，添加了坐标轴标签、标题和网格。
通过这些基于Matlab的模拟，能很直观地看到各种结构动力学现象，真的超级好玩😜！大家也可以自己试试，看看能不能发现更多有趣的东西~
希望这篇博文对大家了解Matlab在结构动力学模拟中的应用有所帮助🤗。
以上就是我这次分享的全部内容啦，下次有好玩的再跟大家继续唠😃。
