fea_modal_analysis.m

来自「FEM tools for caculation of nonlinear pr」· M 代码 · 共 23 行

function  fea_modal_analysis (obj,in_data,L,dof_)

% =========================================================================
% FEM analysis toolbox for solid mechanics. Project started: Aug 2004
% Anton Zaicenco <[e-mail]>
% =========================================================================


if in_data.EL(1,2)==4 | in_data.EL(1,2)==5 | in_data.EL(1,2)==3 | ...
        in_data.EL(1,2)==6 | in_data.EL(1,2)==0 | in_data.EL(1,2)==1 | ...
        in_data.EL(1,2)==2 | in_data.EL(1,2)==9 | in_data.EL(1,2)==51 | ...
        in_data.EL(1,2)==31
    [E_Vec_ok, wn] = modal_analysis ( obj , L, in_data);
    resp.modal.wn = 2*pi./wn; 
    i = size(E_Vec_ok);
    plot_modal_shape(in_data,E_Vec_ok(:,i(2)+1-in_data.dynam.MODA)',dof_,3413);
    mode_T_sec = 2*pi/wn(in_data.dynam.MODA); 

    sp_=sprintf('\n');
    title([sp_ 'Mode shape:  '  num2str(in_data.dynam.MODA) '.  T =  ' num2str(mode_T_sec) '  s']);
    disp([sp_ '--------- MODAL ANALYSIS DONE ---------']);
end;