beam_2d.m

结构力学中的有限元例子,包含了7个分类文件夹

function in_data = beam_2d (in_data)

% beam_2d.m  2-D beams

% nodes: [ node# x y ] 
in_data.ND = [   
         1   0    0;
         2   16   0;
         3   32   0;
         4   16   12;
         5   32   12];

in_data.mater.E = 3e7;
in_data.mater.A = 10;
in_data.mater.I = 3e3;

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   4   E  A  I;
        2  0   1   2   E  A  I;
        3  0   2   4   E  A  I;
        4  0   4   5   E  A  I;
        5  0   3   4   E  A  I;
        6  0   2   5   E  A  I;
        7  0   2   3   E  A  I;
        8  0   3   5   E  A  I];
% constrains: [ node# x y ] 
in_data.CON = [  
         1  0  0  0;
         3  0  0  0];
% loads: [node# x y ] 
in_data.LOAD_ = [  
           2      0  -100e2   0;
           5   50e2       0   0];

in_data.MASS = [  
          2  10e4  30e4  10e6;
          4  10e4  30e4  10e6;   
          5  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     = [2];                    % plot mode of vibration