etsd.m

力学计算 力学计算 力学计算 力学计算 力学计算

function [index]=etsd(nd,nnel,ndof)
%----------------------------------------------------------
% 目的:
%     单元对应的系统自由度号
%
%  变量：
%     index - 和单元"iel"对应的系统自由度数
%     iel - 要确定系统自由度的单元
%     nnel - 每个单元的节点数
%     ndof - 每个节点的自由度
%-----------------------------------------------------------
 
 edof = nnel*ndof;
   k=0;
   for i=1:nnel
     start = (nd(i)-1)*ndof;
       for j=1:ndof
         k=k+1;
         index(k)=start+j;
       end
   end