brick8assemble.m

来自「matlab编程实现有限元计算,对有限元学习很有帮助」· M 代码 · 共 14 行

function y=Brick8Assemble(K,k,i,j,l,m,n,o,p,q);
node=[i,j,l,m,n,o,p,q];
for a=1:8
    for b=1:8
        for c=1:3
            for d=1:3
              s=node(a);t=node(b);
              K(3*s+c-3,3*t+d-3)=K(3*s+c-3,3*t+d-3)+k(3*a+c-3,3*b+d-3);
            end
        end
    end
end
y=K;