function in_data = beam_2d (in_data)

% beam_2d.m  2-D beams

% nodes: [ node# x y ] 
in_data.ND = [   
         1   6    0;
         2   12   0;
         3   18   0;
         4   24   0;
         5   30   0;
         6   0   6;
         7   6   6;
         8   12  6;
         9   18  6;
         10  24  6;
         11  30  6;
         12  36  6];

in_data.mater.E = 9e9;
in_data.mater.A = .4;
in_data.mater.I = 1e-4;

E=in_data.mater.E; A=in_data.mater.A; I=in_data.mater.I;
% elements: [ element# el_type node1# node2# E A I  (el_type: 0 - FF, 1 - FP, 2- PF) ] 
in_data.EL = [  
        1  0   1   2   E  A  I;
        2  0   2   3   E  A  I;
        3  0   3   4   E  A  I;
        4  0   4   5   E  A  I;
        5  0   1   6   E  A  I;
        6  0   1   7   E  A  I;
        7  0   1   8   E  A  I;
        8  0   2   8   E  A  I;
        9  0   3   8   E  A  I;
       10  0   3   9   E  A  I;
       11  0   3   10  E  A  I;
       12  0   4   10  E  A  I;
       13  0   5   10  E  A  I;
       14  0   5   11  E  A  I;
       15  0   5   12  E  A  I;
       16  0   6    7  E  A  I;
       17  0   7    8  E  A  I;
       18  0   8    9  E  A  I;
       19  0   9   10  E  A  I;
       20  0   10  11  E  A  I;
       21  0   11  12  E  A  I];
% constrains: [ node# x y ] 
in_data.CON = [  
         6  0  0  0;
        12  0  0  0];
% loads: [node# x y ] 
in_data.LOAD_ = [  
           1      0  -0e2   0;
           2      0  -0e2   0;
           3      0  -200e2 0;
           4      0  -0e2   0;
           5      0  -0e2   0];

in_data.MASS = [  
          1  10e4  30e4  10e6;
          2  10e4  30e4  10e6;
          3  10e4  30e4  10e6;   
          4  10e4  30e4  10e6;
          5  10e4  30e4  10e6;
          6  10e4  30e4  10e6;
          7  10e4  30e4  10e6;
          8  10e4  30e4  10e6;
          9  10e4  30e4  10e6;
         10  10e4  30e4  10e6;
         11  10e4  30e4  10e6;
         12  10e4  30e4  10e6];
% dynamics:       
in_data.dynam.TIMEH    = [ 'bedr.txt' ];         % file name
in_data.dynam.delta_tm = [0.0079];               % times step, s
in_data.dynam.TIMEHDIR = [1 1 0];                % directivity, dof
in_data.dynam.TIMEHM   = [1:length(in_data.EL)*3];     % apply to masses
in_data.dynam.TIMEHPL  = [2*3-1];                % plot reponse displ for these dof's
in_data.dynam.DAMP_C   = [0.06 0.06];            % damping for the first DAMP_F modes
in_data.dynam.DAMP_F   = [3];                    % 1st and this mode will be taken for C calculation
% modal analysis
in_data.dynam.MODA     = [1];                    % plot mode of vibration