fem3d_inp_frame.m

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

function in_data = fem3D_inp_frame (in_data)

% fem3D_inp.m - input file for FEM analysis using BEAM element 
% nodes: [ node# x y z] 
in_data.ND = [   1   0   0   0;
   2   3   0   0;
   3   6   0   0;
   4   0   2   0;
   5   3   2   0;
   6   6   2   0;
   7   0   4   0;
   8   3   4   0;
   9   6   4   0;
   10  0   0   2;
   11  3   0   2;
   12  6   0   2;
   13  0   2   2;
   14  3   2   2;
   15  6   2   2;
   16  0   4   2;
   17  3   4   2;
   18  6   4   2;
   19  0   0   4;
   20  3   0   4;
   21  6   0   4;
   22  0   2   4;
   23  3   2   4;
   24  6   2   4;
   25  0   4   4;
   26  3   4   4;
   27  6   4   4];


in_data.mater.E = 80e8;
in_data.mater.G = 3e8;
in_data.mater.b = 0.3;
in_data.mater.ho = 0.3;

E=in_data.mater.E; G=in_data.mater.G; b=in_data.mater.b; ho=in_data.mater.ho;
% elements: [ element# node1# node2# E A I  (1 - FF, 2 - FP, 3- PF) ] 
in_data.EL = [  
   1 3   1   10   E  G  b  ho;
   2 3   2   11   E  G  b  ho;
   3 3   3   12   E  G  b  ho;
   4 3   4   13   E  G  b  ho;
   5 3   5   14   E  G  b  ho;
   6 3   6   15  E  G  b   ho;
   7 3   7   16  E  G  b   ho;
   8 3   8   17  E  G  b   ho;
   9 3   9   18  E  G  b   ho;
   10 3  10  19  E  G  b   ho;
   11 3  11  20  E  G  b   ho;
   12 3  12  21  E  G  b   ho;
   13 3  13  22  E  G  b   ho;
   14 3  14  23  E  G  b   ho;
   15 3  15  24  E  G  b   ho;
   16 3  16  25  E  G  b   ho;
   17 3  17  26  E  G  b   ho;
   18 3  18  27  E  G  b   ho;
   19 3  10  11  E  G  b   ho;
   20 3  11  12  E  G  b   ho;
   21 3  13  14  E  G  b   ho;
   22 3  14  15  E  G  b   ho;
   23 3  16  17  E  G  b   ho;
   24 3  17  18  E  G  b   ho;
   25 3  19  20  E  G  b   ho;
   26 3  20  21  E  G  b   ho;
   27 3  22  23  E  G  b   ho;
   28 3  23  24  E  G  b   ho;
   29 3  25  26  E  G  b   ho;
   30 3  26  27  E  G  b   ho;
   31 3  10  13  E  G  b   ho;
   32 3  13  16  E  G  b   ho;
   33 3  11  14  E  G  b   ho;
   34 3  14  17  E  G  b   ho;
   35 3  12  15  E  G  b   ho;
   36 3  15  18  E  G  b   ho;
   37 3  19  22  E  G  b   ho;
   38 3  22  25  E  G  b   ho;
   39 3  20  23  E  G  b   ho;
   40 3  23  26  E  G  b   ho;
   41 3  21  24  E  G  b   ho;
   42 3  24  27  E  G  b   ho];
   
% constrains: [ node# x y z xr yr zr] 
in_data.CON = [  
         1  0  0  0  0  0  0;
         2  0  0  0  0  0  0;
         3  0  0  0  0  0  0;
         4  0  0  0  0  0  0;
         5  0  0  0  0  0  0;
         6  0  0  0  0  0  0;
         7  0  0  0  0  0  0;
         8  0  0  0  0  0  0;
         9  0  0  0  0  0  0];
% loads: [node# x y z xr yr zr] 
in_data.LOAD_ = [  
           22   0e1      10e1    0    0         0    0;
           23   0e1      0e1     0    0         0    0;
           24   0e1      10e1    0    0         0    0];
% masses: [node# x y z xr yr zr] 
in_data.MASS = [  
          10  10e3  30e3  10e3  10e6  10e6  10e6;
          11  10e3  30e3  10e3  10e6  10e6  10e6;
          12  10e3  30e3  10e3  10e6  10e6  10e6;
          13  10e3  30e3  10e3  10e6  10e6  10e6;
          14  10e3  30e3  10e3  10e6  10e6  10e6;
          15  10e3  30e3  10e3  10e6  10e6  10e6;
          16  10e3  30e3  10e3  10e6  10e6  10e6;
          17  10e3  30e3  10e3  10e6  10e6  10e6;
          18  10e3  30e3  10e3  10e6  10e6  10e6;
          19  10e3  30e3  10e3  10e6  10e6  10e6;
          20  10e3  30e3  10e3  10e6  10e6  10e6;
          21  10e3  30e3  10e3  10e6  10e6  10e6;
          22  10e3  30e3  10e3  10e6  10e6  10e6;
          23  10e3  30e3  10e3  10e6  10e6  10e6;
          24  10e3  30e3  10e3  10e6  10e6  10e6;
          25  10e3  30e3  10e3  10e6  10e6  10e6;
          26  10e3  30e3  10e3  10e6  10e6  10e6;
          27  10e3  30e3  10e3  10e6  10e6  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 0 0 0];       % directivity, dof
in_data.dynam.TIMEHM   = [10:1:27];           % apply to masses
in_data.dynam.TIMEHPL  = [19*6-3];            % 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