huabanyu 发表于 2007-5-10 09:16

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

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

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

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=;
=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)];
为什么画出来的没有非线性阿？

咕噜噜 发表于 2007-5-10 14:09

原帖由 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=;
=ode45('zhihuif ...
什么叫做画出来的没有非线性？

xjzuo 发表于 2007-5-10 16:00

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

huabanyu 发表于 2007-5-10 22:31

回复 #4 xjzuo 的帖子

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

xjzuo 发表于 2007-5-11 08:47

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

handb 发表于 2007-5-11 15:48

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

huabanyu 发表于 2007-5-11 18:28

我基本明白了什么意思，谢谢各位阿！
页: [1]
查看完整版本: 请教高手关于非线性频响的问题