fem_inp_csq1.m

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

function in_data = fem_inp_CSQ1 (in_data)

% fem_inp_CSQ1.m - input file for FEM analysis using CSQ element 

% nodes: [ node# x y ] 
%in_data = struct('ND',[],'mater',[],'EL',[],'CON',[],'LOAD_',[],'MASS',[],'dynam',[]);
in_data.ND = [ 1 0 0;
   2 7  0;
   3 10 0;
   4 0  4;
   5 7  4;
   6 10 4;
   7 0  10;
   8 7  10;
   9 10 10];

in_data.mater.E     = 3e7;
in_data.mater.h     = 0.6;
in_data.mater.miu   = 0.3;
in_data.mater.rhoX  = 300; % density of material - X
in_data.mater.rhoY  = 300; % density of material - Y

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, 5 - CSQ) ] 
in_data.EL = [
   1  5  1  2  5  4   E  ;
   2  5  2  3  6  5   E  ;
   3  5  4  5  8  7   E  ;
   4  5  5  6  9  8   E  ];

% constrains: [ node# x y ] 

in_data.CON = [ 1 0 0;
   2 0 0;
   3 0 0];
% loads: [node# x y ] 
in_data.LOAD_ = [4  00e1  00e1;
         7  00e1 00e1;
         8  00e1 00e1;
         9  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 1];                % directivity, dof
in_data.dynam.TIMEHM   = [4:9];                % apply to masses
in_data.dynam.TIMEHPL  = [17];                % 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     = [1];                  % plot mode of vibration