fem_inp_csq2.m

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

function in_data = fem_inp_CSQ2 (in_data)

% fem_inp_CSQ2.m - input file for FEM analysis using CST element 
% nodes: [ node# x y ] 
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;
  29 0 4;
  30 1 4;
  31 2 4;
  32 3 4;
  33 4 4;
  34 5 4;
  35 6 4;
  36 0 5;
  37 1 5;
  38 2 5;
  39 3 5;
  40 4 5;
  41 5 5;
  42 6 5;
  43 0 6;
  44 1 6;
  45 2 6;
  46 3 6;
  47 4 6;
  48 5 6;
  49 6 6;
  50 0 7;
  51 1 7;
  52 2 7;
  53 3 7;
  54 4 7;
  55 5 7;
  56 6 7];

in_data.mater.E = 3e7;
in_data.mater.h = 0.6;
in_data.mater.miu = 0.3;
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;
rhoX = in_data.mater.rhoX;  rhoY = in_data.mater.rhoY;
% elements: [ element# node1# node2# E h miu type: 
%             (1 - FF, 2 - FP, 3- PF, 4 - CST, 5 - CSQ) ] 
in_data.EL = [
   1  5  1  2  9  8   E  [1 1 1 1]  miu rhoX rhoY;
   2  5  2  3 10  9   E  [1 1 1 1]  miu rhoX rhoY;
   3  5  3  4 11 10   E  [1 1 1 1]  miu rhoX rhoY;
   4  5  4  5 12 11   E  [1 1 1 1]  miu rhoX rhoY;
   5  5  5  6 13 12   E  [1 1 1 1]  miu rhoX rhoY;
   6  5  6  7 14 13   E  [1 1 1 1]  miu rhoX rhoY;
   7  5  8  9 16 15   E  [1 1 1 1]  miu rhoX rhoY;
   8  5  9 10 17 16   E  [1 1 1 1]  miu rhoX rhoY;
   9  5 10 11 18 17   E  [1 1 1 1]  miu rhoX rhoY;
  10  5 11 12 19 18   E  [1 1 1 1]  miu rhoX rhoY;
  11  5 12 13 20 19   E  [1 1 1 1]  miu rhoX rhoY;
  12  5 13 14 21 20   E  [1 1 1 1]  miu rhoX rhoY;
  13  5 15 16 23 22   E  [1 1 1 1]  miu rhoX rhoY;
  14  5 16 17 24 23   E  [1 1 1 1]  miu rhoX rhoY;
  15  5 17 18 25 24   E  [1 1 1 1]  miu rhoX rhoY;
  16  5 18 19 26 25   E  [1 1 1 1]  miu rhoX rhoY;
  17  5 19 20 27 26   E  [1 1 1 1]  miu rhoX rhoY;
  18  5 20 21 28 27   E  [1 1 1 1]  miu rhoX rhoY;
  19  5 22 23 30 29   E  [1 1 1 1]  miu rhoX rhoY;
  20  5 23 24 31 30   E  [1 1 1 1]  miu rhoX rhoY;
  21  5 24 25 32 31   E  [1 1 1 1]  miu rhoX rhoY;
  22  5 25 26 33 32   E  [1 1 1 1]  miu rhoX rhoY;
  23  5 26 27 34 33   E  [1 1 1 1]  miu rhoX rhoY;
  24  5 27 28 35 34   E  [1 1 1 1]  miu rhoX rhoY;
  25  5 29 30 37 36   E  [1 1 1 1]  miu rhoX rhoY;
  26  5 30 31 38 37   E  [1 1 1 1]  miu rhoX rhoY;
  27  5 31 32 39 38   E  [1 1 1 1]  miu rhoX rhoY;
  28  5 32 33 40 39   E  [1 1 1 1]  miu rhoX rhoY;
  29  5 33 34 41 40   E  [1 1 1 1]  miu rhoX rhoY;
  30  5 34 35 42 41   E  [1 1 1 1]  miu rhoX rhoY;
  31  5 36 37 44 43   E  [1 1 1 1]  miu rhoX rhoY;
  32  5 37 38 45 44   E  [1 1 1 1]  miu rhoX rhoY;
  33  5 38 39 46 45   E  [1 1 1 1]  miu rhoX rhoY;
  34  5 39 40 47 46   E  [1 1 1 1]  miu rhoX rhoY;
  35  5 40 41 48 47   E  [1 1 1 1]  miu rhoX rhoY;
  36  5 41 42 49 48   E  [1 1 1 1]  miu rhoX rhoY;
  37  5 43 44 51 50   E  [1 1 1 1]  miu rhoX rhoY;
  38  5 44 45 52 51   E  [1 1 1 1]  miu rhoX rhoY;
  39  5 45 46 53 52   E  [1 1 1 1]  miu rhoX rhoY;
  40  5 46 47 54 53   E  [1 1 1 1]  miu rhoX rhoY;
  41  5 47 48 55 54   E  [1 1 1 1]  miu rhoX rhoY;
  42  5 48 49 56 55   E  [1 1 1 1]  miu rhoX rhoY];

% constrains: [ node# x y ] 

in_data.CON = [ 1 0 0;
   2 0 0;
   3 0 0;
   4 0 0;
   5 0 0;
   6 0 0;
   7 0 0];
% loads: [node# x y ] 
in_data.LOAD_ = [22  00e1  00e1;
         28  00e1 00e1;
         53  00e1 00e1;
         50  100e1 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];                % directivity, dof
in_data.dynam.TIMEHM   = [8:length(in_data.EL)*2];             % apply to masses
in_data.dynam.TIMEHPL  = [21*2-1];                 % plot reponse displ for these dof's
in_data.dynam.DAMP_C   = [0.04 0.04];               % 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