dbc.m

来自「有限元法计算平板受集中载荷的程序」· M 代码 · 共 27 行

function [kk,ff]=dbc2(kk,ff,bcdof,bcval)

%----------------------------------------------------------
%  Purpose:
%     为方程[kk]{x}={ff}确定边界条件
%
%  变量:
%     kk - 系统刚度矩阵
%     ff - 载荷
%     bcdof - 边界节点对应的自由度
%     bcval - 边界节点对应的自由度的值
%-----------------------------------------------------------
 
 n=length(bcdof);
 sdof=size(kk);

 for i=1:n
    c=bcdof(i);
    for j=1:sdof
       kk(c,j)=0;
    end

    kk(c,c)=1;
    ff(c)=bcval(i);
 end