[有限元]桁架结构中杆单元的节点序号和单元序号的确定依据

本文链接：https://blog.csdn.net/PriceCheap/article/details/121068743

在使用Matlab进行桁架结构的有限元分析时，确定节点和单元编号的关键是将结构分解为相同单元块。通过避免不同类型的杆件混合，可以简化编号规律。例如，将直杆和斜杆分开处理，确保单元块内部编号的一致性。通过这样的方法，可以有效地编写简洁的赋值语句，如图1和图2结构的解决方案所示。

摘要由CSDN通过智能技术生成

用 Matlab 做简单的桁架结构的有限元分析的时候，会需要一个问题就是，怎么给每个节点和单元编号，以找出最简单的赋值规律，写出最简单的赋值语句
要遵循的原则是，你必须把你的桁架结构分解成若干个最小的完全相同的单元块，并且这些单元块内的各个单元的赋值规律应该是类似的
在这里插入图片描述图1 例子

比如图1结构中，最小的单元块就是由②③④⑤⑥号单元组成的
图2 例子

但是对于图2的结构，你会发现，如果你把这个桁架分成四个块，然后把斜杆和竖直杆水平杆放在一起，你会发现你找出来的规律特别繁琐

在这里插入图片描述

图3 斜杆和直杆放在一个单元块内在这里插入图片描述图4 找规律

%节点坐标矩阵
node=zeros(4*Nspan+4,3);
%x方向坐标
node(1:4:end,1)=(0:Nspan)'*a;
node(2:4:end,1)=(0:Nspan)'*a;
node(3:4:end,1)=(0:Nspan)'*a;
node(4:4:end,1)=(0:Nspan)'*a;
%y方向坐标
node(3:4:end,1)=a;
node(4:4:end,1)=a;
%z方向坐标
node(2:2:end,1)=a;
%单元的端点节点矩阵
element=zeros(12*Nspan+7,2);
element(1)=[1,4]
element(55)=[19.18]

%xz面上的单元

%特殊部分(懒得找规律
element(4)=[1,2]
element(5)=[6,2]
element(6)=[5,1]
element(7)=[5,2]
element(8)=[5,6]
element(9)=[3,4]
element(10)=[8,4]
element(11)=[7,3]
element(12)=[