FEM toolbox for solid mechanics: main_form_data_plot (handles,in_data,pax,plot_dyn) - File Exchange

Code covered by the BSD License

Highlights from
FEM toolbox for solid mechanics

fem(varargin) FEM Application M-file for fem.fig
BRICK_stress (in_data, resp)
BRICK_stress_main (in_data, S...
CST_CSQ_stress (in_data,resp)
Cantilever_Beam(in_data)
D2_CSQ (x1, y1, x2, y2, x3, y... returns a 2D linear constant stress quadrilateral (CSQ) element stiffness
D2_CST (x1, y1, x2, y2, x3, y... returns a 2D linear constant strain triangle (CST) element K & M matrix
D3_BEAM (x1, y1, z1, x2, y2, ... returns a 3D fixed-fixed beam element stiffness matrix in global coordinates
D3_LB (x1, y1, z1, x2, y2, z2... returns a 3D linear constant stress brick (LB) element stiffness matrix
D3_TETRAH (x1, y1, z1, x2, y2... returns a 3D tetrahedron (CST) element K & M matrix
D3_TRUSS (x1, y1, z1, x2, y2,...
DBC_solve (K,L,force,in_data)
FFTsig (T,T2,dt,FNO)
FF_BEAM (x1, y1, x2, y2, E, A... returns a 2D fixed-fixed beam element stiffness & mass matrix in global coordinates
FP_BEAM (x1, y1, x2, y2, E, A... returns a 2D fixed-pinned beam element stiffness & mass matrix in global coordinates
FT_freq (T,dt)
KM_info (K, M)
K_CON (dofN, dof_, obj, in_da...
K_CST (ND, EL, rhoX, rhoY)
K_assembly (in_data, type_ana...
K_pre_allocate (in_data, Nrow...
L=cholesky(K)
LOAD_V (dofN, dof_, in_data, ...
MFC_master_slave (K,L,force,s...
M_2D_beam_assembly (in_data)
M_3D_beam_assembly (in_data)
PF_BEAM (x1, y1, x2, y2, E, A... returns a 2D pinned-fixed beam element stiffness & mass matrix in global coordinates
animate_2D_frame (resp,in_dat...
animate_3D_brick (resp,in_dat...
animate_3D_frame (resp,in_dat...
animate_3D_tetrah (resp,in_da...
animate_3D_truss (resp,in_dat...
animate_dyn_response (resp,in...
animate_mode_shape_2( in_data...
animate_response_tri (resp,in... animate_response_tri.m
beam_2d (in_data)
beam_2d_2 (in_data)
build_object_KM (in_data,type...
dxdtU(t,x,u,params)
dxxdt(t,x)
dxxdt_p(t,x)
dyn_FEM ( in_data, obj, L, LK...
eigFEM (K, M, damp)
eigFEM_s (K, M, N)
fea_dynamic_analysis (obj,in_...
fea_heat_analysis (obj,in_dat...
fea_modal_analysis (obj,in_da...
fea_static_analysis (obj,in_d...
fem3D_brick
fem3D_inp_frame (in_data)
fem3D_tetrah
fem3D_tetrah_cyl (in_data)
fem_inp_CSQ1
fem_inp_CSQ2 (in_data) modal analysis
fem_inp_CST modal analysis
fem_run (in_data,type_analysi...
frame2d_uniform_loads (in_dat...
freq_damp_FEM(K,M,a,b)
in_data_tetr_modal_(in_data) beam_tetr_3d.m 2-D beams
main_form_data_plot (handles,... plot location of DOF for output from dynamic analysis on the main form
matrx_info (A)
modal_analysis ( obj , L, in_...
modal_analysis_large ( obj , ...
odeRK6( dydt, time, y0, u, ep...
plot3D_results(in_data,resp)
plot3D_results_TRUSS(in_data,...
plot_2dbeam_stat_diagram (in_... PLOT STATIC DISPLACEMENTS OF 2D BEAMS
plot_2dbeam_static (in_data,r... PLOT STATIC DISPLACEMENTS OF 2D BEAMS
plot_deform_tri ( in_data, D,...
plot_dyn_2Dbeam (in_data,D,in... plot 2D-FRAME dynamic displacements
plot_dyn_3Dframe(in_data,D,in...
plot_dyn_3Dtruss(in_data,D,in...
plot_dyn_brick ( in_data, D, ...
plot_dyn_tetrah ( in_data, D,...
plot_dyn_tri ( in_data, D, in...
plot_fem_data(in_data,PND) plot input object on GUI main form
plot_fem_static(in_data, obj,... plot static deflections
plot_modal_shape ( in_data, D...
plot_results_brick ( in_data,...
plot_results_heat ( in_data, ... PLOT temperatures (2D)
plot_results_tetrah ( in_data...
plot_results_tri ( in_data, S...
rayleigh_coef (K,M,damp,N)
saveOBJECTeps
surfbrick(p,t) SURFbrick: find surface faces from quad 3d mesh
surftri(p,t)
truss2d (x1, y1, x2, y2, E, A... 2D truss element stiffness and mass matrix in global coordinates
truss_3d (in_data)
type_of_elem (in_data)
View all files

from FEM toolbox for solid mechanics by Anton Zaicenco
The finite element toolbox for solid mechanics with GUI.

main_form_data_plot (handles,in_data,pax,plot_dyn)

function main_form_data_plot (handles,in_data,pax,plot_dyn)

% plot location of DOF for output from dynamic analysis on the main form

axes(pax);

dofN = handles.in_data.dofN;
if isfield(handles.in_data.dynam,'TIMEHPL')
    if length(in_data.dynam.TIMEHPL)~=0
        if handles.in_data.type == '2dbeam'
            vr = ceil(in_data.dynam.TIMEHPL(1)/dofN);
            node = find(in_data.ND(:,1)==vr); hold on;
            if  plot_dyn
                plot(in_data.ND(node,2),in_data.ND(node,3),'b.','MarkerSize',14);
                plot(in_data.ND(node,2),in_data.ND(node,3),'bo','MarkerSize',18);
            end
            labx = (max(in_data.ND(:,2)) / 13); laby = (max(in_data.ND(:,3)) / 13);
            labx = min([labx laby]); laby = labx;
            if in_data.dynam.TIMEHPL(1)-dofN*vr==0
                if  plot_dyn
                    plot(in_data.ND(node,2),in_data.ND(node,3),'bs','MarkerSize',15);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-1
                if  plot_dyn
                    plot([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)-laby*2  in_data.ND(node,3)+laby*2],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-2
                if  plot_dyn
                    plot([in_data.ND(node,2)-labx*2 in_data.ND(node,2)+labx*2], ...
                        [in_data.ND(node,3)  in_data.ND(node,3)],'b-','LineWidth',2);
                end
            end
            rotate3d(gca); set(gcf,'Pointer','arrow');
            pause(0.001); hold off;
        end;
        if handles.in_data.type == '3dbeam'
            vr = ceil(in_data.dynam.TIMEHPL(1)/dofN);
            node = find(in_data.ND(:,1)==vr); hold on;
            if  plot_dyn
                plot3(in_data.ND(node,2),in_data.ND(node,3),in_data.ND(node,4),'b.','MarkerSize',14);
                plot3(in_data.ND(node,2),in_data.ND(node,3),in_data.ND(node,4),'bo','MarkerSize',18);
            end
            labx = (max(in_data.ND(:,2))/9); laby=(max(in_data.ND(:,3))/9); labz=(max(in_data.ND(:,4))/9);
            labx = min([labx laby labz]); laby = labx;  labz=labx;
            if in_data.dynam.TIMEHPL(1)-dofN*vr==0 | ...
                    in_data.dynam.TIMEHPL(1)-dofN*vr==-1 | in_data.dynam.TIMEHPL(1)-dofN*vr==-2
                if  plot_dyn
                    plot3(in_data.ND(node,2),in_data.ND(node,3),in_data.ND(node,4),'bs','MarkerSize',15);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-3
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)  in_data.ND(node,3)],...
                        [in_data.ND(node,4)-labz in_data.ND(node,4)+labz],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-4
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)-laby  in_data.ND(node,3)+laby],...
                        [in_data.ND(node,4) in_data.ND(node,4)],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-5
                if  plot_dyn
                    plot3([in_data.ND(node,2)-labx in_data.ND(node,2)+labx], ...
                        [in_data.ND(node,3) in_data.ND(node,3)],...
                        [in_data.ND(node,4) in_data.ND(node,4)],'b-','LineWidth',2);
                end
            end
            rotate3d(gca); set(gcf,'Pointer','arrow');
            pause(0.01); hold off;
        end;
        if handles.in_data.type == 'cstcsq'
            vr = ceil(in_data.dynam.TIMEHPL(1)/dofN);
            node = find(in_data.ND(:,1)==vr); hold on;
            if  plot_dyn
                plot(in_data.ND(node,2),in_data.ND(node,3),'bo','MarkerSize',14);
                plot(in_data.ND(node,2),in_data.ND(node,3),'b.','MarkerSize',18);
            end
            labx = (max(in_data.ND(:,2)) / 22); laby = (max(in_data.ND(:,3)) / 22);
            labx = min([labx laby]); laby = labx;
            if in_data.dynam.TIMEHPL(1)-dofN*vr==0
                if  plot_dyn
                    plot([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)-laby  in_data.ND(node,3)+laby],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-1
                if  plot_dyn
                    plot([in_data.ND(node,2)-labx in_data.ND(node,2)+labx], ...
                        [in_data.ND(node,3)  in_data.ND(node,3)],'b-','LineWidth',2);
                end
            end
            rotate3d(gca); set(gcf,'Pointer','arrow');
            pause(0.01); hold off;
        end;
        if handles.in_data.type == '3dbrik'
            vr = ceil(in_data.dynam.TIMEHPL(1)/dofN);
            node = find(in_data.ND(:,1)==vr); hold on;
            if  plot_dyn
                plot3(in_data.ND(node,2),in_data.ND(node,3),in_data.ND(node,4),'b.','MarkerSize',14);
                plot3(in_data.ND(node,2),in_data.ND(node,3),in_data.ND(node,4),'bo','MarkerSize',18);
            end
            labx = (max(in_data.ND(:,2))/9); laby=(max(in_data.ND(:,3))/9); labz=(max(in_data.ND(:,4))/9);
            labx = min([labx laby labz]); laby = labx;  labz=labx;
            if in_data.dynam.TIMEHPL(1)-dofN*vr==0
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)  in_data.ND(node,3)],...
                        [in_data.ND(node,4)-labz in_data.ND(node,4)+labz],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-1
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)-laby  in_data.ND(node,3)+laby],...
                        [in_data.ND(node,4) in_data.ND(node,4)],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-2
                if  plot_dyn
                    plot3([in_data.ND(node,2)-labx in_data.ND(node,2)+labx], ...
                        [in_data.ND(node,3) in_data.ND(node,3)],...
                        [in_data.ND(node,4) in_data.ND(node,4)],'b-','LineWidth',2);
                end
            end
            rotate3d(gca); set(gcf,'Pointer','arrow');
            pause(0.01); hold off;
        end;
        if handles.in_data.type == 'bciz  '
            vr = ceil(in_data.dynam.TIMEHPL(1)/dofN);
            node = find(in_data.ND(:,1)==vr); hold on;
            if  plot_dyn
                plot(in_data.ND(node,2),in_data.ND(node,3),'b.','MarkerSize',14);
                plot(in_data.ND(node,2),in_data.ND(node,3),'bo','MarkerSize',18);
            end
            labx = (max(in_data.ND(:,2)) / 9); laby = (max(in_data.ND(:,3)) / 9);
            labx = min([labx laby]); laby = labx; SCz = .5; labz = SCz/2;
            if in_data.dynam.TIMEHPL(1)-dofN*vr==0
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)-laby  in_data.ND(node,3)+laby],...
                        [0 0],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-1
                if  plot_dyn
                    plot3([in_data.ND(node,2)-labx in_data.ND(node,2)+labx], ...
                        [in_data.ND(node,3) in_data.ND(node,3)],...
                        [0 0],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-2
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)  in_data.ND(node,3)],...
                        [0-labz 0+labz],'b-','LineWidth',2);
                end
            end
            rotate3d(gca); set(gcf,'Pointer','arrow');
            pause(0.01); hold off;
        end;
        %----
        
        
        if handles.in_data.type == '3dtrus'
            vr = ceil(in_data.dynam.TIMEHPL(1)/dofN);
            node = find(in_data.ND(:,1)==vr); hold on;
            if  plot_dyn
                plot3(in_data.ND(node,2),in_data.ND(node,3),in_data.ND(node,4),'b.','MarkerSize',14);
                plot3(in_data.ND(node,2),in_data.ND(node,3),in_data.ND(node,4),'bo','MarkerSize',18);
            end
            labx = (max(in_data.ND(:,2))/3); laby=(max(in_data.ND(:,3))/3); labz=(max(in_data.ND(:,4))/3);
            labx = min([labx laby labz]); laby = labx;  labz=labx;
            if in_data.dynam.TIMEHPL(1)-dofN*vr==0
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)  in_data.ND(node,3)],...
                        [in_data.ND(node,4)-labz in_data.ND(node,4)+labz],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-1
                if  plot_dyn
                    plot3([in_data.ND(node,2) in_data.ND(node,2)], ...
                        [in_data.ND(node,3)-laby  in_data.ND(node,3)+laby],...
                        [in_data.ND(node,4) in_data.ND(node,4)],'b-','LineWidth',2);
                end
            end
            if in_data.dynam.TIMEHPL(1)-dofN*vr==-2
                if  plot_dyn
                    plot3([in_data.ND(node,2)-labx in_data.ND(node,2)+labx], ...
                        [in_data.ND(node,3) in_data.ND(node,3)],...
                        [in_data.ND(node,4) in_data.ND(node,4)],'b-','LineWidth',2);
                end
            end
            rotate3d(gca); set(gcf,'Pointer','arrow');
            pause(0.01); hold off;
        end;
            
        %----

        end;
end;

Highlights from FEM toolbox for solid mechanics

Highlights from
FEM toolbox for solid mechanics