exs5.m

用于在matlab平台上进行有限元分析

% example exs5
%----------------------------------------------------------------
% PURPOSE 
%    Analysis of a simply supported beam.
%----------------------------------------------------------------

% REFERENCES
%     G"oran Sandberg 94-03-08 
%     Karl-Gunnar Olsson 95-09-28
%     Ola Dahlblom 2004-09-21
%----------------------------------------------------------------
 echo on

%----- Topology -------------------------------------------------

 Edof=[1  1  2  3  4  5  6;
       2  4  5  6  7  8  9;  
       3  7  8  9 10 11 12];

%----- Stiffness matrix K and load vector f ---------------------

 K=zeros(12);   f=zeros(12,1); f(5)=-10000;

%----- Element stiffness matrices  ------------------------------

 E=2.1e11;     A=45.3e-4;     I=2510e-8;      ep=[E A I];
 ex=[0 3];     ey=[0 0];

 Ke=beam2e(ex,ey,ep)

%----- Assemble Ke into K ---------------------------------------

 K=assem(Edof,K,Ke);

%----- Solve the system of equations and compute support forces -

 bc=[1 0; 2 0; 11 0];
 [a,r]=solveq(K,f,bc);

%----- Section forces -------------------------------------------

 Ed=extract(Edof,a);

 es1=beam2s(ex,ey,ep,Ed(1,:));
 es2=beam2s(ex,ey,ep,Ed(2,:));
 es3=beam2s(ex,ey,ep,Ed(3,:));

%----- Results --------------------------------------------------

 a
 r
 es1
 es2
 es3

%------------------------ end -----------------------------------
 echo off