%这Matlab代码帮助你计算的位移和反应
%梁。您可能需要改变边界条件不同的情况下
%的梁。如果你正在改变元素的数量,你可能需要
力向量变化百分比(F_udl & F_pl) 16日和17日。
clc;清晰;关闭所有;

L = 1;%在米长度
E = 2.1 * 10 ^ 8;%弹性模量KN / m2
我= 2120/100 ^ 4;%惯性矩m4
w = 1000;%均匀分布载荷kN / m
wp = 1000;%集中负荷kN

公元前= (1 2 5);%边界条件

F_udl = w * [- L / 2; (L ^ 2) / 48 - (L ^ 2/48); (L ^ 2/48)];
F_pl = wp * [1, -0.05, 0];

Ne =输入('类型的元素或按“输入”使用默认2元素”);
如果(isempty (Ne))
不= 2;
其他的
结束

Ndof = 2 *(不+ 1);
景深= 1:1:Ndof;

勒= L /东北;

%计算结构的刚度矩阵
Kj = 0 (4, 4);
Kj (1,1) = 6;
Kj(1、2) = 3 *勒;
Kj (1、3) = 6;
Kj(1、4) = 3 *勒;
Kj(2, 2) = 2 *勒^ 2;
Kj(2、3) = 3 *勒;
Kj(2、4) =勒^ 2;
Kj (3) = 6;
Kj(3、4) = 3 *勒;
Kj(4, 4) = 2 *勒^ 2;

2 = 2:4,
jj = 1:第二,
Kj (ii, jj) = Kj (jj, ii);
结束
结束
克= (2 * E *我/ le ^ 3)。* Kj;

公斤= 0 (Ndof Ndof);
2 = 1:不
aa = 2 * ii-1;
公斤(aa aa: aa + 3: aa + 3) =公斤(aa aa: aa + 3: aa + 3) +客;
结束

K_temp =公斤;
jj = 1:长度(BC)
公元前K_temp (:, (jj)) = 0;
K_temp (BC (jj):) = 0;
结束
K_temp (~ (K_temp 2):) = [];
K_temp (:, ~ (K_temp, 1)) = [];
基米-雷克南= K_temp;

FR = F_udl + F_pl;

你= KR \ FR;
UG = [0 0 UR(1)你(2)0 UR (3)] ';
Rc =公斤* (UG);

clc
disp (“mm & rad位移”)
你= (UR(1) * 1000;你的(2),你(3)];
disp (UR)
disp (kN的反应)
反应= [Rc (1) + w * L / 4;Rc (2) + w * (L / 2) ^ 2/12;Rc (5) + w * L / 4);
disp(反应)