材料的装配刚度矩阵

时间:2016-07-21 14:29:16

标签: matlab

我正在研究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

欢迎提出建议。

1 个答案:

答案 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的更多注释。

希望这有帮助。