用赫兹理论计算轴承刚度的MATLAB（牛顿迭代法解）

liumeironglmr

您好：
我想做的是根据赫兹理论、弹性接触理论推导径向力、转速对径向刚度的影响规律。次方法已经进过n多人使用。
实质上就是解一个非线性方程组，未知量为
X1,X2,bi,bo,br,calfai
方程为
F(1)=(A1-X1)^2+(A2-X2)^2-((fi-0.5)*D1+bi)^2;
F(2)=X1^2+X2^2-((fo-0.5)*D1+bo)^2;
F(3)=(Mg*X2/D1-Ko*bo^1.5*X1)/((fo-0.5)*D1+bo)+(Ki*bi^1.5*(A1-X1)-Mg*(A2-X2)/D1)/((fi-0.5)*D1+bi);
F(4)=(Mg*X1/D1+Ko*bo^1.5*X2)/((fo-0.5)*D1+bo)-(Fc+(Ki*bi^1.5*(A2-X2)+Mg*(A1-X1)/D1)/((fi-0.5)*D1+bi));
F(5)=Fri-(Ki*bi^1.5*(A2-X2)+Mg*(A1-X1)/D1)/((fi-0.5)*D1+bi);
F(6)=calfai-(A2-X2)/((fi-0.5)*D1+bi);
由于是非线性方程组，所以大家惯用牛顿迭代法，但是由于由于都是复合函数，中间参数过多，使得我的结果不准确。
我的第四个.m 文件 是利用求得得X1,X2,bi,bo,br,calfai来计算内外圈接触角、接触力和径向刚度，结果都是不正确的。因为X1,X2,bi,bo,br,calfai算出来就不靠谱了，前5个应该在1微米 -100微米之间 ，第6个是在0.9左右，可是我得到的结果根本不是一个数量级。导致后面算角度、刚度也都错的。
请高手帮忙看看我的程序出什么问题了

zhouzhouxiwei

我也在做滚动轴承刚度计算。有机会多交流。

zhouzhouxiwei

八六七一三七七一〇。

等待苍老

你好，你的程序可以给我参考一下吗