animate_response_tri.m

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

% animate_response_tri.m     -  animate response of CST FEM model 
function animate_response_tri (resp,in_data);

%cd D:\matlabR12\work\fem\anim\3;

dofN = 2;   % dof's per node
dof_ = size(in_data.ND)*dofN;
warning off;
    XXX = max(max(abs( resp.dynamic.DYN_D(:,1:2:end))));
    YYY = max(max(abs( resp.dynamic.DYN_D(:,2:2:end))));

for k=1:size(resp.dynamic.DYN_D,1)
%    if k<2 SIGsys_dyn = CST_CSQ_stress (EL, ND, resp.DYN_D(k,:).*5e3,E, h, miu); end;% stress at this time
    plot_deform_tri( in_data, resp, k, dof_, XXX, YYY, 2);
    figure (2); hold on;
    jo = num2str(in_data.dynam.delta_tm*k); jo = strcat(jo, '  s');  % time, s 
    title(jo); hold off;
    pause(0.00001);
% write movie 'M' (<300) frames into jpg files
    if k<100 A='a'; end; if k<200 & k>=100 A='b'; end;
    if k<300 & k>=200 A='c'; end;
    eXt='.jpg'; B=strcat(A,int2str(k),eXt);
%    print('-djpeg','-r110',B);
%    IMWRITE(2,B,'jpg');
end


warning on;