悬臂梁挠度与转角Matlab程序

以下是悬臂梁挠度和转角的 Matlab 程序：

clear all; % 清空所有变量
clc; % 清空命令窗口

L = 10; % 梁的长度
E = 2.1e11; % 梁的杨氏模量
I = 1.333e-5; % 梁的截面惯性矩
w = 1000; % 分布载荷，单位：N/m

x = 0:0.01:L; % 离散化梁的长度
M = w * x.^2 / 2; % 弯矩分布
theta = M ./ (E * I); % 转角分布
y = (w * x.^4) / (24 * E * I) .* (L - x); % 挠度分布

% 绘制梁的弯矩分布图
subplot(2,1,1);
plot(x, M);
title('弯矩分布图');
xlabel('梁的长度（m）');
ylabel('弯矩（N·m）');

% 绘制梁的转角分布图
subplot(2,1,2);
plot(x, theta);
title('转角分布图');
xlabel('梁的长度（m）');
ylabel('转角（rad）');

% 绘制梁的挠度分布图
figure;
plot(x, y);
title('挠度分布图');
xlabel('梁的长度（m）');
ylabel('挠度（m）');

在这个程序中，我们首先定义了梁的长度、杨氏模量、截面惯性矩和分布载荷等参数。然后，我们使用 `x = 0:0.01:L` 将梁的长度离散化，并计算出梁的弯矩、转角和挠度分布。最后，我们使用 `plot` 函数绘制了梁的弯矩、转角和挠度分布图。

注意：这个程序中的计算结果仅适用于悬臂梁，对于其他类型的梁，需要根据实际情况进行修改。