声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1887|回复: 4

[非线性振动] 关于BOUC-WEN滞回曲线不闭合的问题

[复制链接]
发表于 2011-11-18 17:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
请各位大侠帮助看下我的程序,为什么出来的图形不闭合??
动力学模型如下:见附件

***绘图程序******
tspan=[0:0.01:2];
y0=zeros(3,1);
[t,y]=ode45(@bouc_wenfun,tspan,y0);
plot(y(:,1),y(:,3));
hold on;
box off;
******函数******
function dy=bouc_wenfun(t,y)
%系统参数
A=1;k=1;c=0.1;
f=5*sin(3*t);
n=1;
belta=0.25;garma=0.65;
f=5*cos(3*t);
%函数定义
dy=zeros(3,1);
dy(1)=y(2);
dy(2)=(f-(1-a)*k*y(3)-a*k*y(1)-c*y(2))/m;
dy(3)=A*y(2)-belta*abs(y(2))*(abs(y(3))^(n-1)*y(3)-garma*y(2)*abs(y(3))^(n));
end


bouc-wen动力学模型图.jpeg

bouc-wen动力学模型.doc

34.5 KB, 下载次数: 9

回复
分享到:

使用道具 举报

 楼主| 发表于 2011-11-18 17:27 | 显示全部楼层
回复 1 # zhengshuhe 的帖子

我自己细想一下,提出两种途径是否可以解决?
1、把方程分阶段展开解析解在分段绘图?
2、不调用ODE45函数,自己编写4阶RK算法?
 楼主| 发表于 2011-11-21 17:22 | 显示全部楼层
有没高手帮助看下啊,急~~~~~~~~~~~
 楼主| 发表于 2011-11-23 11:30 | 显示全部楼层
回复 3 # zhengshuhe 的帖子

程序出了点错,搞定了,可以闭合
发表于 2012-3-28 16:51 | 显示全部楼层
回复 4 # zhengshuhe 的帖子

我的程序也是这样的,是哪里有问题呢?急死了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 13:43 , Processed in 0.067946 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表