plot_fem_static.m

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

function plot_fem_static(in_data, resp)

% plot static deflections

if in_data.EL(1,2)==0 | in_data.EL(1,2)==1 | in_data.EL(1,2)==2 dofN = 3; end; % 2D-beam
if in_data.EL(1,2)==3 dofN = 6; end; % 3D-beam  elements
if in_data.EL(1,2)==4 dofN = 2; end; % CST  elements
if in_data.EL(1,2)==5 dofN = 2; end; % CSQ  elements
if in_data.EL(1,2)==6 dofN = 3; end; % "6" - 3D-BRICK element (8-nodes)
if in_data.EL(1,2)==9 dofN = 3; end; % BCIZ elements

dof_ = size(in_data.ND,1)*dofN; % total dof

%    if plot_bend_tri (..., A )<100 - don't plot loads & constrains
if in_data.EL(1,2)==9 plot_bend_tri ( in_data, resp, dof_, 120); end; % BCIZ
if in_data.EL(1,2)==4  |  in_data.EL(1,2)==5                  % 4 - CST, 5 - CSQ
    SIGsys = CST_CSQ_stress (in_data, resp);
    SF_ = zeros(size(in_data.EL,1), 4);
    SF_(1:size(in_data.EL,1),1)=[1:size(in_data.EL,1)]';
    for i=1:size(in_data.EL,1)
        % disp([ sprintf('\t%d\t%1.7f\t',i,SIGsys(i*3-(3-1):i*3)) ]); % plot static stresses
        SF_(i,2:end) = [SIGsys(i*3-(3-1):i*3)]';
    end;
    save FEM_STRESS.txt SF_ -ascii; format;
    disp('---  Elements Stresses saved to "FEM_STRESS.txt": [elem#  x  y  angle(deg)]');
    plot_results_tri(  in_data, SIGsys, resp, dof_, 110 );
end;  
if in_data.EL(1,2)==3  plot3D_results(in_data,resp);  end;    % 3 - 3D-beam
if in_data.EL(1,2)==6  % 3D-brick
    [SIGsys] = BRICK_stress (in_data, resp);
    [SIG_main] = BRICK_stress_main (in_data, SIGsys);
    SF_ = zeros(size(in_data.EL,1), 4);
    SF_(1:size(in_data.EL,1),1)=[1:size(in_data.EL,1)]';
    for i=1:size(in_data.EL,1)
        %disp([ sprintf('\t%d\t%1.7f\t',i,SIG_main(i*3-(3-1):i*3)) ]);
        SF_(i,2:end) = [SIGsys(i*3-(3-1):i*3)]';
    end;
    save FEM_STRESS.txt SF_ -ascii; format;
    disp('--- "BRICK" Principal Stresses saved to "FEM_STRESS.txt"');
    plot_results_brick (in_data, resp, dof_, SIG_main);
end;
if in_data.EL(1,2)==0 | in_data.EL(1,2)==1 | in_data.EL(1,2)==2
    plot_2dbeam_static (in_data,resp);  % "0/1/2" "FF/FP/PF" 2D-beams
end;