| 查看: 809 | 回复: 4 | |||
[求助]
刚度耦合振动系统结果不收敛
|
|
小弟做轮胎动力学的，得到振动方程如图所示，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  Matlab程序和自由度c、z的计算结果.jpg [ Last edited by feixiaolin on 2015-7-7 at 09:18 ] |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
ansys接触分析不收敛
已经有8人回复- 流固耦合问题不收敛，流场计算不对称
已经有7人回复
- 接触分析不收敛
已经有3人回复
- ansys接触分析时，面面接触的接触刚度有什么物理意义吗？
已经有15人回复
- ansys接触分析中接触刚度及容忍穿透调整
已经有4人回复
- PCG求解器迭代不收敛
已经有4人回复
- ansys接触计算不收敛问题
已经有16人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 2930
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2015-07-08 17:21:36
3楼2015-07-08 17:56:05
4楼2015-07-08 17:58:16
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 2930
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
5楼2015-07-08 18:04:58