非线性油膜力的程序如何写

oiloil

用matlab编非线性油膜力的表达式，可是其中的导数部分咋编阿？我觉得应该用循环语句变，可不知道咋编？烦请大侠指点一二，最简单的长轴承模型就可以。

[ 本帖最后由 eight 于 2008-3-18 20:59 编辑 ]

kangtj

你应该把公式贴出来，然后别人才知道具体什么问题。有什么解决办法啊！

eight

没办法，谁让楼主（包括为数不少的其他会员）认为我们（包括管理员和高手）都是神仙啊，以为这样一摆就可以完事

[ 本帖最后由 eight 于 2008-3-19 13:02 编辑 ]

sogooda

原帖由 eight 于 2008-3-19 13:01 发表

                               
登录/注册后可看大图

没办法，谁让楼主（包括为数不少的其他会员）认为我们（包括管理员和高手）都是神仙啊，以为这样一摆就可以完事

呵呵。。。老八幽默

oiloil

不好意思，那天有些匆忙，急匆匆的就把问题发了上去了，没考虑大侠的感受阿，哎，惭愧惭愧:@L ，现请听我慢慢道来：

程序目的是求解方程[M]{U}"+[C]{U}'+[K]{U}={Q}，广义力考虑不平衡力和油膜力，一开始油膜力我用的是八系数线性模型，都是常数，没啥问题，可是当油膜力改成非线性模型时，就不知道该如何用matlab实现了，因为非线性油膜力模型是变量，含有导数项。比如无限长轴承模型(直角坐标下无量纲表达式)，如附件。

初来乍到，烦请大侠指点一二。。。



[ 本帖最后由 sigma665 于 2008-3-20 21:09 编辑 ]