fem_inp_cst.m

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

function in_data = fem_inp_CST (in_data)

% fem_inp_CST.m - input file for FEM analysis using CST element 
% wing-shape structure
% nodes: [ node# x y ] 
in_data.ND = [   
   1   1   1;
   2   2.6 1;
   3   4.8 1;
   4   7   1;
   5   8   1;
   6   1.6 2.5;
   7   3   2.5;
   8   4.9 2.5;
   9   7   2.5;
   10  8   2.5;
   11  2.3 4;
   12  3.4 4;
   13  5   4;
   14  7   4;
   15  8   4;
   16  2.8 5;
   17  3.8 5;
   18  5.1 5;
   19  7   5;
   20  8   5;
   21  3.2 6;
   22  4.2 6;
   23  5.4 6;
   24  7.2 6;
   25  8.2 6;
   26  3.7 7;
   27  4.7 7;
   28  5.7 7;
   29  7.3 7;
   30  8.3 7;
   31  4.2 8;
   32  5.2 8;
   33  6.2 8;
   34  7.6 8;
   35  8.6 8;
   36  4.6 9;
   37  5.5 9;
   38  6.5 9;
   39  7.8 9;
   40  8.8 9;
   41  5   10;
   42  6   10;
   43  7   10;
   44  8   10;
   45  9   10];
   

in_data.mater.E = 9e6;
in_data.mater.h = 0.2;
in_data.mater.miu = 0.3;
in_data.mater.rhoX = 2000; in_data.mater.rhoY = 2000;

E = in_data.mater.E; h = in_data.mater.h; miu = in_data.mater.miu;
% elements: [ element# node1# node2# E h miu type: 
%            (1 - FF, 2 - FP, 3- PF, 4 - CST) ] 
in_data.EL = [  
   1   4   1   2   6   E  h  miu;
   2   4   6   2   7   E  h  miu;
   3   4   2   3   7   E  h  miu;
   4   4   7   3   8   E  h  miu;
   5   4   3   4   8   E  h  miu;
   6   4   8   4   9   E  h  miu;
   7   4   9   4   5   E  h  miu;
   8   4   9   5  10   E  h  miu;
   9   4   6   7  11   E  h  miu;
  10   4   11  7  12   E  h  miu;
  11   4   7   8  12   E  h  miu;
  12   4   12  8  13   E  h  miu;
  13   4   8   9  13   E  h  miu;
  14   4   13   9  14   E  h  miu;
  15   4   14   9  10   E  h  miu;
  16   4   14  10  15   E  h  miu;
  17   4   11  12  16   E  h  miu;
  18   4   16  12  17   E  h  miu;
  19   4   12  13  17   E  h  miu;
  20   4   17  13  18   E  h  miu;
  21   4   13  14  18   E  h  miu;
  22   4   18  14  19   E  h  miu;
  23   4   19  14  15   E  h  miu;
  24   4   19  15  20   E  h  miu;
  25   4   16  17  21   E  h  miu;
  26   4   21  17  22   E  h  miu;
  27   4   17  18  22   E  h  miu;
  28   4   22  18  23   E  h  miu;
  29   4   18  19  23   E  h  miu;
  30   4   23  19  24   E  h  miu;
  31   4   19  20  24   E  h  miu;
  32   4   24  20  25   E  h  miu;
  33   4   21  22  26   E  h  miu;
  34   4   26  22  27   E  h  miu;
  35   4   22  23  27   E  h  miu;
  36   4   27  23  28   E  h  miu;
  37   4   23  24  28   E  h  miu;
  38   4   28  24  29   E  h  miu;
  39   4   24  25  29   E  h  miu;
  40   4   29  25  30   E  h  miu;
  41   4   26  27  31   E  h  miu;
  42   4   31  27  32   E  h  miu;
  43   4   27  28  32   E  h  miu;
  44   4   32  28  33   E  h  miu;
  45   4   28  29  33   E  h  miu;
  46   4   33  29  34   E  h  miu;
  47   4   34  29  30   E  h  miu;
  48   4   34  30  35   E  h  miu;
  49   4   36  31  32   E  h  miu;
  50   4   36  32  37   E  h  miu;
  51   4   32  33  37   E  h  miu;
  52   4   37  33  38   E  h  miu;
  53   4   33  34  38   E  h  miu;
  54   4   38  34  39   E  h  miu;
  55   4   39  34  35   E  h  miu;
  56   4   39  35  40   E  h  miu;
  57   4   36  37  41   E  h  miu;
  58   4   41  37  42   E  h  miu;
  59   4   37  38  42   E  h  miu;
  60   4   42  38  43   E  h  miu;
  61   4   38  39  43   E  h  miu;
  62   4   43  39  44   E  h  miu;
  63   4   44  39  40   E  h  miu;
  64   4   44  40  45   E  h  miu];
  
in_data.mater.E=0; in_data.mater.rho=0; in_data.mater.miu=0;
  
% constrains: [ node# x y ] 
in_data.CON = [  
   1  0  0;
   2  0  0;
   3  0  0;
   4  0  0;
   5  0  0];
% loads: [node# x y ] 
in_data.LOAD_ = [ 
    6   100e1       00e1;
   11   100e1       00e1;
   16   100e1       00e1;
   21   100e1       00e1;
   26   100e1        00e1;
   31   100e1        00e1;
   36   100e1        00e1;
   25   00e1         00e1;
   42   00e1        00e1;
   43   00e1        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];
% dynamics:       
in_data.dynam.TIMEH    = [ 'bedr.txt' ];       % file name
in_data.dynam.delta_tm = [0.0079];             % times step, s
in_data.dynam.TIMEHDIR = [1 0];                % directivity, dof
in_data.dynam.TIMEHM   = [1:length(in_data.EL)*2];     % apply to masses
in_data.dynam.TIMEHPL  = [41];                 % 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