eval(demosdt('echoon'));comgui('closeall');

% this example illustrates the use of the FEMESH preprocessor to build a
% solid model of a U-beam, compute the associated modes, and display strain
% energy levels
%
% See also demos gartfe, beambar, d_truss, demo_fe, d_plate
%          doc   fem, feco, dfeplot

femesh('reset');
% Base nodes
FEnode=[1 0 0 0  -.5 -.5 0;2  0 0 0  -.5+1/6 -.5 0;3 0 0 0  -.5 .5-1/6 0
        4 0 0 0  -.5+1/6 .5-1/6 0;5 0 0 0  -.5 .5 0;6 0 0 0 -.5+1/6 .5 0
        7 0 0 0 .5-1/6 .5 0;8 0 0 0 .5 .5 0;9 0 0 0 .5-1/6 .5-1/6 0
        10 0 0 0 .5 .5-1/6 0;11 0 0 0 .5-1/6 -.5 0;12 0 0 0 .5 -.5 0];

% define U shaped base surface
FEelt=[Inf abs('quad4');4 6 5 3 1 1;9 10 8 7 1 1];
FEel0=[Inf abs('quad4');1 2 4 3 1 1;11 12 10 9 1 1];
femesh(';divide -new 1 5;addsel;');
FEel0=[Inf abs('quad4');4 6 7 9 1 1];
femesh(';divide -new 1 4;addsel;');

femesh('join group 1:3');

% extrude the base to form the beam model
femesh(';selgroup1;extrude 10 0 0 .25;orientel0');

model=femesh('model0');
model.pl = m_elastic('dbval 1 steel');
model.il = p_solid('dbval 1 d3 2');

% One can the compute the clamped free modes by 
% - defining boundary conditions,
% - solving for modes 
% - displaying them

model=fe_case(model,'FixDof','Clamped end','z==0')
model=stack_set(model,'info','EigOpt',[6 10 0 11]);
def = fe_eig(model,[6 10 0 11]);

% Finally, one can display the model and strain energy distribtion

cf=feplot; cf.model=model;cf.def=def; 
fecom(';view3;view s-90;colorfacew;showpatch');
fecom('colordata ener k');

%-----------------------------------------------------------------
eval(demosdt('echooff'))

%       Etienne Balmes
%       Copyright (c) 1990-2006 by SDTools, All Rights Reserved
%       $Revision: 1.6 $  $Date: 2006/03/21 17:13:40 $