dynamicstability.m

薄壁结构有限元计算

R=Nc(1)^2+Nc(2)^2+(Ix+Iy)/A;
px=(pi/n/l)^2*E*Ix;
py=(pi/n/l)^2*E*Iy;
mju=1-J/Ip;
pw=G*Ip*(pi^2*E*If+mju*G*J*n^2*l^2)/(pi^2*E*If+mju*G*Ip*n^2*l^2)/R;
p=zeros(1,1000);
pcr=p;
pm=zeros(1000,4);
for i=1:1000
    pm(i,:)=[(Nc(2)-x(129,i))^2+(Nc(1)-x(128,i))^2-R (px+py+pw)*R-px*(Nc(2)-x(129,i))^2-py*(Nc(1)-x(128,i))^2 -R*(px*py+py*pw+pw*px) px*py*pw*R];
    pcr(i)=min(roots(pm(i,:)));
    if (x(128,i)^2+x(129,i)^2)>0
        for j=1:n
             p(i)=p(i)+mem*9.8*(x(7*j-4,i)^2+x(7*j-5,i)^2)^0.5/(x(128,i)^2+x(129,i)^2)^0.5;
        end
    else
        p(i)=mem*9.8*n/2;
    end
end