我正在研究fem弹性平面应变问题。我现在需要使用3种不同的材料,我很难进行刚度矩阵的组装。
这是我的装配的一个小鳕鱼,有一种材料:
edof = length(index); % dof per element (x,y)
for i = 1:edof
ii = index(i);
for j = 1:edof
jj = index(j);
K(ii,jj) = K(ii,jj) + k(i,j); % k elem stiffness matrix
end
end
欢迎提出建议。
答案 0 :(得分:0)
如果你想使用有限元方法,你首先需要找到Hilbert-1空间的有限子空间V(具有紧致支持)。
让a成为问题的微分算子。然后你的刚度矩阵就足够了:
A(i,j) = a(phi(i),phi(j))
其中phi(i)
和phi(j)
是V(您的子空间)的元素。
通过这个lecture notes我学习了有关FEM的方法(第23页将更详细地解释上述内容),在这里您可以找到关于numerical PDEs的更多注释。
希望这有帮助。