回复 #15 无水1324 的帖子
初始值不能为零要不e就等于零了
%判断是否碰摩
if e<dd
rud=0;
else
rud=1;
end
if e==0
Px=0;
Py=0;
else
%碰摩力
Px=-c*(e-dd)*kc/e*(x(5)-f*x(7));
Py=-c*(e-dd)*kc/e*(f*x(5)+x(7));
end
这样就不会有哪个警告了
谢谢啊
回复 #14 octopussheng 的帖子
ode15不能运行啊? 应该是下面几个 ode23, ode113, ode15s, ode23s, ode23t, ode23tb没有ode15的,呵呵!
回复 #18 octopussheng 的帖子
都不行啊!怎么办呢??帮忙想想办法?:'( 那没有办法了,你只能把你的参数一个个试着调整调整了! 有谁做过着方面的,帮忙看看!谢谢了!
回复 #21 秋月 的帖子
这样的问题经常遇到,以前也总结过一些,但是具体是什么因素引起的就不清楚了建议:
1、检查里面的参数;
2、看一下表达式是否写错;
3、ode算法换一下。
仅仅是建议,不一定能够解决问题。我以前也是这样,但是都是出在上面三个问题上
回复 #22 无水1324 的帖子
谢谢你的提议!newmark方法怎么计算啊?
回复 #23 秋月 的帖子
具体的没有做过,你可以看一下论坛相关的帖子回复 #24 无水1324 的帖子
请问你有转子碰磨这方面的例子吗?能否给我个例子,
我怎么做不出来,
谢谢了!
回复 #5 sssssxxxxx921 的帖子
你好!能帮我写个转子碰摩程序吗?
我怎么作不出来啊!都快郁闷死了!
谢谢了! 上边的程序无水、oct这些牛人已经给你提出了很好的建议
不过我昨天运行程序时偶然发现转子质量也会极大的影响分岔和混沌是否出现的,
你试着调整一下转子质量
我现在两台机子在运行程序(占着别人一台),写出来也没地方试运行调整
另外能力有限对碰摩也不太熟悉啊呵呵
[ 本帖最后由 无水1324 于 2007-10-7 13:33 编辑 ]
回复 #25 秋月 的帖子
程序没有了,只是你可以照书上写一个,搞清楚参数之间的关系应该比较容易。遇到问题可以再一起讨论。 原帖由 无水1324 于 2007-10-7 13:33 发表 http://www.chinavib.com/forum/images/common/back.gif程序没有了,只是你可以照书上写一个,搞清楚参数之间的关系应该比较容易。遇到问题可以再一起讨论。
帮忙看看这个刚性支撑的分岔图有什么问题?
没有考虑轴承.
也没有考虑碰摩条件,加上碰摩条件就得不到这样得结果了!
回复 #27 sssssxxxxx921 的帖子
你可以给我一个简单得碰摩例子吗?我好不容易得到一个上面一个结果
但是发现程序没有加碰摩条件!
谢谢了,帮帮忙!!