目录
1、此文是基于matlab的受轴弯剪的梁单元的结构求解,代码如下:
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