k_con.m

来自「结构力学中的有限元例子,包含了7个分类文件夹」· M 代码 · 共 21 行

function [Kred,L,Rch] = K_CON (dofN, dof_, obj, in_data)

% eliminate constrained dof's from the global stiffness matrix

mark = 121212125.2222111; % marker for constrains 
for i=1:size(in_data.CON,1)
    for k=1:dofN
        if in_data.CON(i,1+k)==0 
            obj.Ksys.Kgl(in_data.CON(i,1)*dofN-(dofN-k),:)=mark; 
            obj.Ksys.Kgl(:,in_data.CON(i,1)*dofN-(dofN-k))=mark; end;
    end;
end;
for k=1:length(obj.Ksys.Kgl)/dofN % # of nodes
    % along "Rch" row & column "mark" will be searched
    if k~=in_data.CON(:,1)   Rch=k*dofN-(dofN-1); break;  end;
end;
L  = (obj.Ksys.Kgl(Rch,:)~=mark);
Y  = obj.Ksys.Kgl(L,:);
Kred = Y(:,L);           % reduced global stiff matrix