step.m

基于力密度法的索杆张力结构找形程序

function stepp=Step(xx0,yy0,zz0,xx1,yy1,zz1,xx2,yy2,zz2,xx3,yy3,zz3,stepp)
% if stepp<0.2
    aa=stepp;
% else
%     aa=0.12;
% end

if max(abs(zz0-zz2))<0.9*aa&max(abs(zz1-zz3))<0.9*aa
    if max(abs(yy0-yy2))<0.9*aa&max(abs(xx0-xx2))<0.9*aa&max(abs(xx1-xx3))<0.9*aa&max(abs(yy1-yy3))<0.9*aa
        stepp=stepp*0.5;
    end
end


% if time>40&time<100
%     stepp=0.1;
% end
% if time>=100&time<300
%     stepp=0.01;
% end
% if time>=300&time<800
%     stepp=0.001;
% end
% if time>=800&time<1600
%     stepp=0.0002;
% end
% if time>16000
%     stepp=0.000001;
% end
% if max(abs(dertX))<0.1*stepp&max(abs(dertY))<0.1*stepp&max(abs(dertZ))<0.1*stepp
%     stepp=0.1*stepp;
% end