请教高手关于非线性频响的问题

huabanyu

在matlab中，如何用数值方法做出非线性的频响曲线？

[ 本帖最后由 eight 于 2007-5-10 10:01 编辑 ]

huabanyu

主程序：
global f E d33 u33 w w0
f=10;
E=40*10^9;
d33=2.97*10^-8;
u33=5*4*pi*10^-7;
figure
for w=5000:20.0:20000
     t0=0:0.001:0.02
     y0=[0,0.01,0.0000001];
[t,y]=ode45('zhihuifun',t0,y0,[])
aa=max((y(:,1)));
plot(w/w0,aa,'r.')
hold on
end

调用程序：
function dy=zhihuifun(t,y,flag)
global f E d33 u33 w w0
m=1.38;k=2*10^7;l=0.4;A=1.21*pi*10^-4;kesai=0.1;
w0=sqrt(k/m-(E*A/(l*m))*(1-E*d33^2/u33));
dy=[y(2);
    -w0^2*y(3)-2*kesai^2*w0*y(2)+f*sin(w*t);
    y(2)-0.15*sqrt(y(2)^2)*y(3)-0.85*y(2)*sqrt(y(3)^2)];
为什么画出来的没有非线性阿？

咕噜噜

原帖由 huabanyu 于 2007-5-10 09:54 发表
主程序：
global f E d33 u33 w w0
f=10;
E=40*10^9;
d33=2.97*10^-8;
u33=5*4*pi*10^-7;
figure
for w=5000:20.0:20000
     t0=0:0.001:0.02
     y0=[0,0.01,0.0000001];
[t,y]=ode45('zhihuif ...

什么叫做画出来的没有非线性？

xjzuo

我画了一下图形(稍微修改了一下程序),应该没有什么问题.

huabanyu

能运行出来，就是做出来的图应该是一对多的图形，是不是频响程序不对啊？谢谢

xjzuo

这恐怕要将原问题贴一下.  
"一对多的图形" 具体是指...?

handb

是指频响曲线出现了跳跃（频响曲线发生弯曲，向左弯表示振动系统具有软特性、向右弯表示具有硬特性）。这个过程刚度要发生变化，而程序中的刚度是个不变量，不可能出现多值的情形。

huabanyu

我基本明白了什么意思，谢谢各位阿！