eval(demosdt('echoon'))

%
% This example illustrates the use of mixed element types with FEMESH
%
% See also demos gartfe, d_truss, demo_fe, d_plate, d_ubeam
%          doc   fem, feco, dfeplot


femesh('reset'); % initialize FEMESH


FEnode = [1  0 0 0  0 0 0;2  0 0 0  1 0 0;
          3  0 0 0  0 1 0;4  0 0 0  1 1 0];
FEelt=[];

% make an outer frame with beam elements

femesh(';object beamline 1 2 0 4 3 1;addsel;info FEelt');

% make diagonals with bar elements

femesh(';object beamline 1 4 0 2 3;addsel;set group2 name bar1');

% repeat the cell 10 times and keep the result

femesh(';selgroup1:2;repeatsel 10 1 0 0;');
FEelt=FEel0;

% add a beam at the end of the last bay

femesh('objectbeamline',femesh('findnode x==10'));
femesh addsel;


% Define the properties of both groups - - - - - - - - - - - - - - -

% Set the MatID and SecID values for elements
femesh(';setgroup1 3 mat1sec1;setgroup2 mat2 sec2');

model=femesh;


model.pl=[m_elastic('dbval 1 aluminum');
          m_elastic('dbval 2 steel')];
          
model.il=[p_beam('dbval 1 rectangle .1 .1')
          p_beam('dbval 2 circle .03')];

% compute modes and display the result

def= fe_eig(model,[6 20 1e3]);

cf=feplot; cf.model=model; cf.def=def;
fecom(';scaledef1;view2;ch4');

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

%	Etienne Balmes
%       Copyright (c) 1990-2003 by SDTools, All Rights Reserved
%       $Revision: 1.4 $  $Date: 2005/10/26 09:51:08 $