inputdata.m

是个好东西

function inputdata(fid_in)
%-----读入数据
global node element material K bc bc_number nf_number nf  node_number element_number %（定义全局变量）
node_number = fscanf( fid_in, '%d', 1 ) ;  %read node number 
node = zeros( node_number, 3) ;
for i=1:1:node_number
   nn = fscanf( fid_in, '%d', 1 ) ;
   node( i, : ) = fscanf( fid_in, '%f', [1,3] ) ;  % read node definition
end
element_number = fscanf( fid_in, '%d', 1 ) ;  % read element number
element = zeros( element_number, 3 ) ;        
for i=1:1:element_number
   ne = fscanf( fid_in, '%d', 1 ) ;
   element( i,:) = fscanf( fid_in, '%d', [1,3] ) ; % read element definition
end
material_number = fscanf( fid_in, '%d', 1 );  % read material number
material = zeros( material_number, 2 ) ;
for i=1:1:material_number
   nm = fscanf( fid_in, '%d', 1 ) ;
   material( i, : ) = fscanf( fid_in, '%f', [1,2] )  ; % read materials definition
end
bc_number = fscanf( fid_in, '%d', 1 );   % read boundary conditions number
bc = zeros( bc_number, 3 ) ;
for i=1:1:bc_number                        % read boundary condition definition
   bc( i, 1 ) = fscanf( fid_in, '%d', 1 ); 
   bc( i, 2 ) = fscanf( fid_in, '%d', 1 ); 
   bc( i, 3 ) = fscanf( fid_in, '%f', 1 ); 
end
nf_number=fscanf( fid_in, '%d', 1 ); % read node force number
nf = zeros( nf_number, 3 );
for i=1:1:nf_number                      % read node force definition
   nf( i, 1 ) = fscanf( fid_in, '%d', 1 ); 
   nf( i, 2 ) = fscanf( fid_in, '%d', 1 ); 
   nf( i, 3 ) = fscanf( fid_in, '%f', 1 );
end
return