frame.m

是个好东西

function frame  %  （定义新的函数truss）
% This function analyze the problem of plane stress use triangle elements 
% the variables lists:
%       node  -------  node definition
%       element -----  element definition
%       material ----  material definition
%       K -----------  global stiffness matrix                -|
%       f -----------  global node force vector                |  [K]{delta}={f}
%       delta -------  global node displacement vector -|
%       bc ----------  boundary condition 
%       nfp ----------  结点力信息
%       nfe-----------   非结点力信息
%       file_in -----  input file name
%       file_out ----  output file name
%       fid_in ------  input file ID
%       fid_out -----  output file ID
global node element material K delta bc bc_number nfp nfe nfp_number nfe_number nf  node_number element_number %（定义全局变量）
file_in='myinput.txt';           %input('input file name:','s' );
file_out='myoutput.txt';          %input('output file name:','s');
fid_in=fopen(file_in,'r');   % open input file 
readdata(fid_in);
fclose( fid_in ) ;    %close the input file which has been opened
K=sparse(node_number*3,node_number*3);%（创建稀疏矩阵K）
for ie=1:1:element_number
   k=StiffnessMatrix( ie ) ;  % (调用函数StiffnessMatrix)
   AssembleStiffnessMatrix( ie, k ) ;% (调用函数AssembleStiffnessMatrix)
end
f=loadff; %-------计算荷载列向量
bound;    %-------引入边界条件
delta = K\f ; %（解方程组,得位移delta）
fid_out = fopen( file_out, 'w' ) ; % open the output file
fprintf( fid_out, '节点位移\n');
for i=1:1:node_number
   fprintf( fid_out, 'node = %5d, x = %.8e, y = %.8e, theta = %.8e\n', ...
   i, delta((i-1)*3+1), delta((i-1)*3+2) , delta((i-1)*3+3)) ;  % output the node displacement
end
NQM(fid_out);
fclose(fid_out) ; % close the output file
return ;%（function truss结束）