数学建模(3)-基于matlab的受轴弯剪的梁单元的结构求解-有限元-FEM

目录

1、此文是基于matlab的受轴弯剪的梁单元的结构求解，代码如下：

2、生成的结果 ：​

3、视频资料见链接：

---

1、此文是基于matlab的受轴弯剪的梁单元的结构求解，代码如下：

% 基于matlab的受轴弯剪的梁单元的结构求解；
% 输入量为：节点数目，坐标； 单元节点编号，自由度编号; 已知力，位移； 未知位移；
clc      
clear all 

num_nod=3;  %节点数目；
ele_nod=[1 2;2 3]; %单元节点编号；
nod_coor=[0 0;12 0; 24.73 12.73] ; %节点坐标，向右为x轴，向下为y轴；
ele_dof=[1 2 3 4 5 6;4 5 6 7 8 9]; %单元自由度编号，对应于子块搬家对应的家；
force=zeros(3*num_nod,1);
displacement=zeros(3*num_nod,1);
stiffness=zeros(3*num_nod);
ele_num=size(ele_nod,1);
force(5) = 1000000000.0; %已知节点力；
displacement(1)=0.0;  %约束条件；
displacement(2)=0.0;
displacement(3)=0.0;
displacement(7)=0.0;
displacement(8)=0.0;
displacement(9)=0.0;
unknown_f_a=[4;5;6]; %未知节点位移；

% 截面基本参数
b=0.5; 
h=0.5; 
E=2.1e11;


%初使化所要的矩阵，F K 位移矩阵；

A=b*h;
I=b*h^3/12; 

for e = 1 : ele_num
    L(e) = s