exs3.m

来自「用于有限元分析的matlab工具箱」· M 代码 · 共 56 行

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

% REFERENCES
%     G"oran Sandberg 94-03-08 
%     Karl-Gunnar Olsson 95-09-28
%----------------------------------------------------------------
 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 reactions ------

 bc=[1 0; 2 0; 11 0];	[a,Q]=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
 Q
 es1
 es2
 es3

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