你找到解决的办法了吗?我面临着同样的问题:
有限元法:策划应力场在每个T3元素
22 Ansichten(30天)的
Altere Kommentare anzeigen
下午好,我写了在matlab有限元法。我的代码找到一个孔板的位移事业中提交一个内部压力。我的网是由53个节点。所有节点的坐标(xi, yi) 53 x2矩阵中包含“
coordonnees
”。那些53 80 T3三角形元素节点形式,连接80 x3矩阵中包含被称为“
connectivite
”。每80个元素我都发现,有限元模型,强调80年提交那些元素。强调保存在两个80 x1向量称为“
sigma_x
”和“
sigma_y
”。现在,我想画网格,和颜色每个三角形元素根据他们的σ值,但我不知道如何做到这一点呢?之一,你有一个想法如何做到这一点呢?(我发现其他有限元用户共享的代码,做类似的事情,但在他们的代码,他们阴谋位移场为例,给出了每一个网格的节点,因此,卡拉季奇的节点的颜色。我的问题是向量的不同的方式
sigma_x
给出了应力为每一个元素,而不是每一个节点)。
非常感谢你的时间和帮助!
达米安•萨默
0 Kommentare
Antworten (2)
darova
我2021年7月26日
这是一个例子:
(x, y) = meshgrid (2:0.2:2);% x, y的范围
z = x ^ 2 + y ^ 2 - 4;% z值
h0 =冲浪(x, y, z,“可见”,“关闭”);%创建冲浪(不显示)
h1 = surf2patch (h0,“三角形”);%转换为表面修补(三角形)
h1。facevertexCData = h1.vertices (: 3);%每个节点分配的颜色作为Z值
补丁(h1,“facecolor”,“平”)%可视化
1 Kommentar
darova
我2021年7月26日
更好的模型
%建立模型
(x1, y1) = pol2cart((0:90:270) * 180(π/ 3);%广场/板(未关的轮廓)
(x2, y2) = pol2cart (linspace(0, 2 *π- 0.1 (20),1);%圆(未关的轮廓)
阻止gd2 =[2;长度(x2); x2 (:), y2 (:));%圆几何描述
gd1 =阻止gd2 * 0;
gd1(1:2 + 2 *长度(x1)) =(2;长度(x1); x1(:);日元(:));%板几何描述
科幻小说=“c1”;%固体fomrula(减圆板)
ns = char (C2的,“C1”)”;%名称空间(圆板)
dl = decsg([阻止gd2 gd1],科幻,ns);%分解
[p, e t] = initmesh (dl);%创建网格
%创建补丁数据
z = p (1:) ^ 2 + p (2:)。^ 2 - 4;% calcualte z坐标
h1。面临= t (1:3:)”;%从三角形创建的脸
h1。顶点= [p;z]”;%顶点
h1。faceVertexCData = z (:);%根据z coord颜色
补丁(h1,“facecolor”,“平”)
视图(60)
翻译的
您可以在欧什一张网站来自der folgenden Liste auswahlen: