应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 刚度耦合振动系统结果不收敛
31/1返回列表
查看: 814  |  回复: 4
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

傻傻先生

新虫 (小有名气)

[求助] 刚度耦合振动系统结果不收敛

小弟做轮胎动力学的，得到振动方程如图所示，matlab程序和计算结果如图所示。

计算结果c、z随时间呈线性增加并不收敛，见第2个图，不是我想要的。
查找原因后发现问题在于刚度矩阵力k1=k2，当k1>k2时结果就像阻尼系统里收敛。
重新计算模型后发现k1的确与k2相等。

请问：c、z随时间呈线性增加是唯一解吗？我希望结果是像像阻尼系统里收敛。
如果是唯一解，那么这个解得物理含义是怎么样的？
求大神指导~会的虫友请看我发的帖子，那里有50个金币悬赏，还可以再加，希望虫友们
为小弟我解惑！！！万分感谢~

为方便各位虫友，特附上matlab程序
function z=odefun3(t,x)
m1=44.8;
k1=7.5091e6;
ma=70.7945;
k2=7.5091e6;
lamida=0.5;%阻尼比
c1=2*lamida*sqrt(m1*k1);
ca=2*lamida*sqrt(ma*k2);
z1=x(3);
z2=x(4);
z3=-c1/m1*x(3)-k1/m1*x(1)+k2/m1*x(2);
z4=-ca/ma*x(4)+k2/ma*x(1)-k2/ma*x(2)+943.9270/ma;
z=[z1;z2;z3;z4];
end

[t f]=ode45(@odefun3,0:0.01:1,[0 0 0 0])
figure(1)
plot(t,f(:,1))
figure(2)
plot(t,f(:,2))

振动方程
刚度耦合振动系统结果不收敛
666.PNG


刚度耦合振动系统结果不收敛-1
Matlab程序和自由度c、z的计算结果.jpg

[ Last edited by feixiaolin on 2015-7-7 at 09:18 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2015-07-07 00:16:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
塑性形变？特殊过程【自举；自激】？
赞一下 回复此楼
2楼2015-07-08 17:21:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
3楼: Originally posted by 傻傻先生 at 2015-07-08 17:56:05
虫友，能详细说一下吗？我把时间设为0-100秒，响应也是线性增大的

那可能是方程本身不定或多解。

[ 发自小木虫客户端 ]
赞一下(1人) 回复此楼
5楼2015-07-08 18:04:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 傻傻先生 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭