E=210e3;%弹性模量 t=2;%梁高度 b=100;%梁宽度 I=1/12*b*t^3;%惯性矩 L=1;%单元长度 TL=1000;%梁总长度 n=TL/L;%d单元个数 %创建单元刚度矩阵 for i=1:n ke(:,:,i)=danyuan(E,I,L); end; %创建全局刚度矩阵 K=zeros((n+1)*2); for i=1:n K=qvanjv(K,ke(:,:,i),i,i+1); end %加入边界条件悬臂 k=[K(3:(n+1)*2,3:(n+1)*2)]; %加入外力 f=zeros(2*n,1); for i=1:n f(2*(i-1)+1)=1*i*L;%线性分布力 end u=k\f; %节点位移 U=zeros((n+1)*2,1); U(3:(n+1)*2)=u; %求单元力 for i=1:n Fe(i,:)=ke(:,:,i)*U((i-1)*2+1:(i-1)*2+4); end F=K*U;%求节点力

这是一段MATLAB代码，用于求解悬臂梁在线性分布力作用下的节点位移和单元力。代码中包括以下步骤：

1. 定义悬臂梁的弹性模量、高度、宽度、惯性矩、长度等参数。

2. 创建单元刚度矩阵。这里采用了一个名为“danyuan”的函数，用于计算梁单元的刚度矩阵。

3. 创建全局刚度矩阵。将所有单元的刚度矩阵按照其位置组装成全局刚度矩阵。这里采用了一个名为“qvanjv”的函数，用于将两个单元的刚度矩阵组装到全局刚度矩阵中。

4. 加入边界条件，即悬臂端固定。

5. 加入外力，这里采用线性分布力的形式，即从悬臂端到梁另一端力逐渐增大。

6. 求解节点位移。将全局刚度矩阵和外力向量带入线性方程组求解函数“\”，得到节点位移。

7. 求解单元力。根据单元刚度矩阵和节点位移，计算每个单元的力。

8. 求解节点力。根据全局刚度矩阵和节点位移，计算每个节点的力。

这段代码可以用于悬臂梁的静力分析，对于其他类型的结构或者非线性问题，需要相应地修改代码。