% F=0:0.001:0.28;% F从0~0.28变化，步长为0.001
% w0=1;a=0.168;b=-1/2;c=1/2;T=2*pi/w0;
% axis([0 0.28 -0.4 1.6])% 坐标轴范围的选取
% hold on% 把图形画在一幅图上
% % 循环叠代
function []=fenchatu()
l=0:100:7000;% F从0~0.28变化，步长为0.001
Nb=3;Nr=63;w=pi*Nb*Nr/30;
T=2*pi/w;
% % w0=1;a=0.168;b=-1/2;c=1/2;T=2*pi/w0;
% a=190;
% w0=1.0;
% % F=400;
% T=2*pi/w0;
% b=-1000;
% c=1000;
% A=2000;
% rou=7850;E=210*10^9;r=0.127;Nb=3;l=7000;w=pi*Nb*Nr/30;D1=0.127;D0=0.1175;
% k=E*((D1)^2-(D0)^2)/(4*l);T=2*pi/w;
% alpha=0.1;beta=0.001;
% m1=1/4*rou*l*((D1)^2-(D0)^2);c=alpha*m1+beta*k;
% A=pi*r^2;
% B1=c/(rou*A);
% B2=k/(rou*A)+4*pi^2*E/(3*l^3);
% F=0.01/3*(w^2-k/(rou*A))*sin(w*t)-0.01/3*(w/(rou*A))*cos(w*t);
% axis([0 400 -0.4 1.6])% 坐标轴范围的选取
hold on% 把图形画在一幅图上
% 循环叠代
for j=1:length(l)% j从1到F的长度变化，步长为1
[t,y]=ode45(@zongxiangzhendong,[0:T/33:100*T],[0,0]);
% 龙格库塔法（子函数名，单调积分区间，初始条件，用odeset建立优化选项，参数传递）
plot(l(j),y(3000:3:end,1),'k.','marker','.','markersize',1);
% （画F-周期位移点对应图，蓝色，点，点的尺寸，1）
end
xlabel('T/N.m');% x轴坐标
ylabel('位移X/m');% y轴坐标
grid on
function f =zongxiangzhendong(t,y)
f = zeros(2,1);
% c=2.2*10^4;k=2.2*10^6;
c=20;k=900;
rou=7850;E=210*10^9;r=0.127;Nb=3;Nr=63;l=7000;w=pi*Nb*Nr/30;D1=0.127;D0=0.1175;
% k=E*((D1)^2-(D0)^2)/(4*l);
% alpha=0.1;beta=0.001;
% m1=1/4*rou*l*((D1)^2-(D0)^2);c=alpha*m1+beta*k;
A=pi*r^2;
B1=c/(rou*A);
B2=k/(rou*A)+4*pi^2*E/(3*l^3);
F=0.01/3*(w^2-k/(rou*A))*sin(w*t)-0.01/3*(w/(rou*A))*cos(w*t);
f(1) = y(2);
f(2) = -F-B1*y(2)-B2*y(1);
评论0