eval(demosdt('echoon'));

% Demonstration of direct finite element model declaration and basic uses
% of FEMESH on truss examples.
%
% Type doc('sdt/fem') at the MATLAB prompt to access the documentation 
% associated to this demo.


femesh('reset');; 

% This is a direct declaration of a 2 bay truss model

%      NodeID   unused   x y z
FEnode=[ 1      0 0 0    0 1 0;
         2      0 0 0    0 0 0;
         3      0 0 0    1 1 0;
         4      0 0 0    1 0 0;
         5      0 0 0    2 0 0;
         6      0 0 0    2 1 0;
         7      0 0 0    1 1 1]; % reference node

FEelt=[
 % declaration of element group for longerons
    Inf     abs('beam1') 
 %node1 node2   MatID ProID nodeR ,  zeros to fill the matrix 
    1       3      1    1     7       0
    3       6      1    1     7       0
    2       4      1    1     7       0
    4       5      1    1     7       0
 % declaration of element group for diagonals
    Inf     abs('beam1')
    2       3      1    2     7       0
    4       6      1    2     7       0
 % declaration of element group for battens
    Inf     abs('beam1')
    3       4      1    3     7       0
    5       6      1    3     7       0];

femesh('plotelt'); fecom(';textnode;view2');

demosdt('pause'); %----------------------------------------------------------

% In practice you will generally prefer to generate the model using FEMESH
% in this case this can be done as follows

 femesh('reset'); % just to make sure that variables are declared as global
 FEelt =[]; % we want to erase the prevous model

 % build the first bay

 FEnode=[1 0 0 0  0 0 0;2 0 0 0    0 1 0;
         3 0 0 0  1 0 0;4 0 0 0    1 1 0];
 femesh('objectbeamline 1 3 0 2 4 0 3 4 0 1 4')

 % put it in the main model, translate, add the second bay, see result

 femesh(';addsel;transsel 1 0 0;addsel;info;plotelt');

demosdt('pause'); %----------------------------------------------------------

% You can do fairly complicated things with FEMESH. The following creates
% a 3-D truss based on the same basic cell.

 femesh('removeelt group2');
 femesh('divide group 1 InNode 1 4');
 femesh('set group1 name bar1');
 femesh(';selgroup2 1;repeatsel 10 1 0 0;addsel');
 femesh(';rotatesel 1 60 1 0 0;addsel;');
 femesh(';selgroup3:4;rotatesel 2 -60 1 0 0;addsel;');
 femesh(';selgroup3:8;plotel0');
 fecom(';triaxon;view3;view y+180');

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

%       Etienne Balmes  02/01/92, 07/20/98
%       Copyright (c) 1990-2004 by SDTools,
%       All Rights Reserved.