材料のアセンブリ剛性マトリックス

Question

私はfemの弾性面歪み問題に取り組んでいます。 3つの異なる素材を使用する必要があり、剛性マトリックスのアセンブリを行うのが難しいです。ここで

は、材料の1種類と私のアセンブリの小さなタラです：

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 

提案は歓迎されています。

Answer 1

有限要素法を使用する場合は、最初にHilbert-1空間の有限部分空間V（コンパクトサポートあり）を見つける必要があります。

aを問題の微分演算子とします。次に、あなたの剛性マトリックスで十分：

A(i,j) = a(phi(i),phi(j))

どこphi(i)とphi(j)はV（あなたの部分空間）の要素です。

このlecture notesでは、FEMの方法について学びました（23ページで詳細を説明します）。ここではnumerical PDEsに関する詳細をご覧いただけます。

これが役に立ちます。