matlab计算惯性矩,梁单元有限元计算程序(matlab)

%此程序计算杆的总刚度矩阵及节点位移，分五步完成：输入各单元数据、计算单元刚度矩阵

%组集总刚度矩阵、计算输出总刚度矩阵 、计算输出节点位移

%2011.4

%输入个单元数据

%输入单元节点编号矩阵，每一行代表该单元的节点编号，即用分号将各单元分开

%用逗号将每个单元内的节点分开，按单元编号顺序排列

%如[1,2;2,3;1,3]表示三个杆中的节点编号分别为(1,2)、(2,3)、(1,3)

cod=input('please input the node of each element in order:');

%计算单元个数，nm为单元个数

[nm,nmn]=size(cod);

%输入各单元的角度

alpha=input('please input the angle (degree) of each element in order:');

%输入节点坐标，每一行代表该节点的坐标，按节点编号顺序排列，即用分号将节点分开

%用逗号将每个节点的坐标分开，按单元编号顺序排列

%如[1,2;2,3;1,3]表示三个节点的横纵坐标分别为(1,2)、(2,3)、(1,3)

con=input('please input the coordinates (m) of each node in order:');

%计算结点个数，nn为结点个数

[nn,nnn]=size(con);

%输入单元的弹性模量

E=1e9*input('please input E array (GPa) of each element