matlab编写的有限元杆的计算

这是我们上有限元程序这门课的时候，老师给我们的一个有限元程序，计算杆或梁的受力分析，包括完整的程序，感觉很详细，还有例子用word说明。

% ------------------------------------------------------------------------|
%                                                                         |
% MAE4700-5700 Finite Element Analysis for Mechanical & Aerospace Design |
%                                                                         |
% Copyright: Cornell University （this software should not be used without |
% written permission）                                                     |
%                                                                         |
% Authors: N. Zabaras （zabaras@cornell.edu） & Xiang Ma （xm25@cornell.edu） |
%                                                                         |
% ------------------------------------------------------------------------|
%

function Assemble;
include_variables;

% We go through each truss element compute its stiffness and then perform
% the assebmly process 

ldof = nen*ndof;       % degrees of freedom in each element 
                       % （number of nodes x number of degrees of freedom
                       % in each node）

for elmID = 1 : nel
   
   Ke = zeros（ldof）;                  % initialize element stiffness matrix
   fe = zeros（ldof1）;                % initialize element load vector
   
   [Kefe] = BeamElement（elmID Ke fe）;
      
   AssembleGlobalMatrix（elmIDKefe）; % Place the element stiffness in the 
                                      % global matrix
   
end

clear Ke fe;

% Once the global stiffness is computed we generate its sparse format
% Constructs the global stiffness matrix in sparse format using IJK

K = sparse（I（1:ntriplets） J（1:ntriplets） X（1:ntriplets） neq neq）;

clear global I J X

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1724  2008-09-06 13:42  successful\Assemble.m

     文件       2996  2008-09-06 13:29  successful\AssembleGlobalMatrix.m

     文件       1846  2010-06-29 16:15  successful\BeamElement.m

     文件       2219  2008-09-23 14:36  successful\exact.m

     文件       2654  2010-06-29 16:15  successful\include_variables.m

     文件       3981  2010-06-29 17:28  successful\InputData.m

     文件       2053  2008-09-23 14:35  successful\InputGrid.m

     文件       1335  2008-09-06 13:28  successful\Main.m

     文件       2740  2008-08-21 09:24  successful\NodalSoln.m

     文件       1366  2008-09-06 13:06  successful\PlotBeam.m

     文件       5024  2010-06-29 16:47  successful\postprocessor.m

     文件       1934  2008-09-09 22:06  successful\preprocessor.m

    ..A..H.       162  2010-06-30 08:58  successful\~$ampleProblem-Beam.doc

     目录          0  2010-06-30 08:59  successful

     文件     101888  2008-09-10 15:37  successful\ExampleProblem-Beam.doc

----------- ---------  ---------- -----  ----

               131922                    15