fem_inp_bend.m

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

function in_data = fem_inp_bend (in_data)

% fem_inp_bend.m - input file for FEM analysis using CST element 
% nodes: [ node# x y ] 

in_data = struct('ND',[],'mater',[],'EL',[],'CON',[],'LOAD_',[],'MASS',[],'dynam',[]);
in_data.ND = [   1   0   0;
   2   1  0;
   3   2  0;
   4   3  0;
   5   4  0;
   6   5  0;
   7   6  0;
   8   0  1;
   9   1  1;
   10  2  1;
   11  3  1;
   12  4  1;
   13  5  1;
   14  6  1;
   15  0  2;
   16  1  2;
   17  2  2;
   18  3  2;
   19  4  2;
   20  5  2;
   21  6  2;
   22  0  3;
   23  1  3;
   24  2  3;
   25  3  3;
   26  4  3;
   27  5  3;
   28  6  3];

in_data.mater.E   = 3e9;
in_data.mater.h   = 0.03;
in_data.mater.miu = 0.33;
in_data.mater.rhoX  = 2000; % density of material - X
in_data.mater.rhoY  = 2000; % density of material - Y

E = in_data.mater.E; h = in_data.mater.h; miu = in_data.mater.miu;
% elements: [ element# node1# node2# node3# E h miu 'type: 9 - BCIZ' ] 
in_data.EL = [   1   9   1   9   8   E  h  miu;
   2   9   1   2   9   E  h  miu;
   3   9   2   10  9   E  h  miu;
   4   9   2   3   10  E  h  miu;
   5   9   3   11  10  E  h  miu;
   6   9   3   4   11  E  h  miu;
   7   9   4   12  11  E  h  miu;
   8   9   4   5   12  E  h  miu;
   9   9   5   13  12  E  h  miu;
   10  9   5   6   13  E  h  miu;
   11  9   6   14  13  E  h  miu;
   12  9   6   7   14  E  h  miu;
   13  9   8   16  15  E  h  miu;
   14  9   8   9   16  E  h  miu;
   15  9   9   17  16  E  h  miu;
   16  9   9   10  17  E  h  miu;
   17  9   10  18  17  E  h  miu;
   18  9   10  11  18  E  h  miu;
   19  9   11  19  18  E  h  miu;
   20  9   11  12  19  E  h  miu;
   21  9   12  20  19  E  h  miu;
   22  9   12  13  20  E  h  miu;
   23  9   13  21  20  E  h  miu;
   24  9   13  14  21  E  h  miu;
   25  9   15  23  22  E  h  miu;
   26  9   15  16  23  E  h  miu;
   27  9   16  24  23  E  h  miu;
   28  9   16  17  24  E  h  miu;
   29  9   17  25  24  E  h  miu;
   30  9   17  18  25  E  h  miu;
   31  9   18  26  25  E  h  miu;
   32  9   18  19  26  E  h  miu;
   33  9   19  27  26  E  h  miu;
   34  9   19  20  27  E  h  miu;
   35  9   20  28  27  E  h  miu;
   36  9   20  21  28  E  h  miu];



% constrains: [ node# x y ] 
in_data.CON = [  1  0  0  0;
         8  0  0  0;
         15 0  0  0;
         22 0  0  0;
         2  0  0  0;
         3 0  0  0;
         4 0  0  0;
         5 0  0  0;
         6 0  0  0;
         7 0  0  0];
         
% loads: [node# x y ] 
in_data.LOAD_ = [ 
          28   7e1         00e1   00e1;
          26   00e1        -1e1   00e1];


% masses: [node# x y z] 
in_data.MASS = [  
          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;
          13 10e4  30e4  10e6;
          14 10e4  30e4  10e6;
          15 10e4  30e4  10e6];
% accelerations at nodes       
in_data.dynam.TIMEH    = [ 'bedr.txt' ];               % file name
in_data.dynam.delta_tm = [0.0079];                     % times step, s
in_data.dynam.TIMEHDIR = [1 0 0];                      % directivity, dof
in_data.dynam.TIMEHM   = [1:1:length(in_data.EL)*3];   % apply to masses
in_data.dynam.TIMEHPL  = [14*3-2];                        % 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     = [3];                           % plot mode of vibration 

% displacement vector: [ w1 taux1 tauy1 w2 taux2 tauy2 w3 taux3 tauy3 ]