readdata.m

来自「教学程序」· M 代码 · 共 46 行

function readdata(fid_in)
%-----读入数据
global node element material K delta bc bc_number nfp nfe nfp_number nfe_number nf  node_number element_number %（定义全局变量）
node_number = fscanf( fid_in, '%d', 1 ) ;  %read node number 
node = zeros( node_number, 2) ;
for i=1:1:node_number
   nn = fscanf( fid_in, '%d', 1 ) ;
   node( i, : ) = fscanf( fid_in, '%f', [1,2] ) ;  % 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, 3 ) ;
for i=1:1:material_number
   nm = fscanf( fid_in, '%d', 1 ) ;
   material( i, : ) = fscanf( fid_in, '%f', [1,3] )  ; % 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
nfp_number=fscanf( fid_in, '%d', 1 );  % read node force number
nfp = zeros( nfp_number, 3 ) ;
for i=1:1:nfp_number                      % read node force definition
   nfp( i, 1 ) = fscanf( fid_in, '%d', 1 ) ;
   nfp( i, 2 ) = fscanf( fid_in, '%d', 1 ) ;
   nfp( i, 3 ) = fscanf( fid_in, '%f', 1 ) ;
end
nfe_number=fscanf( fid_in, '%d', 1 ); % read node force number
nfe = zeros( nfe_number, 6 ) ;
for i=1:1:nfe_number                      % read node force definition
   nfe( i, 1 ) = fscanf( fid_in, '%d', 1 ) ;
   nfe( i, 2 ) = fscanf( fid_in, '%d', 1 ) ;
   nfe( i, 3 ) = fscanf( fid_in, '%f', 1 ) ;
   nfe( i, 4 ) = fscanf( fid_in, '%f', 1 ) ;
   nfe( i, 5 ) = fscanf( fid_in, '%f', 1 ) ;
   nfe( i, 6 ) = fscanf( fid_in, '%f', 1 ) ;
end
return