MATLAB的答案

计算四面体有限元的刚度矩阵

36视图(30天)
内森•韦尔奇
内森•韦尔奇 2020年12月4日
评论道: Biswabhanu Puhan2023年3月30日
我写一些有限元模拟和网格排序。
我现在把刚度矩阵在一起,肯定觉得这是一个通用的任务已经完成后,但谷歌和查看文件交换似乎不是!有几件事,但没有提到桁架和东西需要一个明显的输入和给一个简单的刚度矩阵。
有人知道任何好的代码,一个网状tetrahedralisation矩阵T (Ne, 4)和顶点V (Nv), 3)并产生线性刚度矩阵K (Nv, Nv) ?显然更高的订单也会大但线性足以开始。
2的评论
显示 1年长的评论
内森•韦尔奇
内森•韦尔奇 2020年12月29日
谢谢,是的,这似乎是标准的方法。也相当不错的介绍有限元法解决如何这么做我想的!

登录置评。

在回答这个问题。

答案(1)

内森•韦尔奇
内森•韦尔奇 2020年12月29日
似乎没有任何明显的来源,我创造了我自己。它计算刚度矩阵的线性元素四面体网格。与所有在线代码,使用在你的自己的风险。如果任何人都可以描述如何更好的计算,特别是如何扩展到二阶或(梦想)任意的顺序,我仍然非常感兴趣。
函数K = LinearStiffnessMatrix3(节点,元素)
%得到大小
Nn =大小(节点,1);
不=大小(元素,1);
%创建稀疏矩阵的存储(i, j)指数和价值观
第1 = 0;
如果= 0 (16 * Ne, 1);% i刚度矩阵的指数
Sj = 0 (16 * Ne, 1);% j刚度矩阵的指数
Sk = 0 (16 * Ne, 1);% (i, j)的价值元素的刚度矩阵
% 4指数与指数
2 = 1:4;霁= 1:4;[ii,霁]= ndgrid (ii, ji);
2 = 2 (:);霁=霁(:);
%网格遍历每个元素——这样做就好了
%循环——或者这可能最好使用汇编语言的原因
因为我= 1:不
eleNodes =元素(我,1:4);
%使用简单的线性项方法的集成
B =[节点(eleNodes:) ';(1,1,1,1)];
C =发票(B);
ve = abs(侦破(B));
Kl2 = C (: 1:3) * C (: 1:3)。* abs (ve) / 6;
%存储所有16项
第1 =第1(结束)+ (1:16);
Si(第1)= eleNodes (ii);
Sj(第1)= eleNodes (ji);
Sk(第1)= Kl2 (:);
结束
%建立刚度矩阵——matlab稀疏运营商一起自然地补充道
% i和j相同的所有条款,即所有贡献第i个节点
%从节点j
K =稀疏(Si, Sj, Sk,神经网络,神经网络);
结束
1评论
Biswabhanu Puhan
Biswabhanu Puhan 2023年3月30日
谢谢韦尔奇先生,你觉得二阶元素的配方(C3D10有限元分析)?请让我知道。

登录置评。

类别

找到更多的在结构分析在帮助中心和文件交换

标签

下载188bet金宝搏

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!