10个弹簧振子系统，受到光顺后库伦摩擦，结果有问题，

lex9317

10个弹簧振子系统，受到光顺后库伦摩擦，结果有问题。
function f=func4(t,y)
f=zeros(20,1)

f(19)=y(20)
f(17)=y(18)
f(15)=y(16)
f(13)=y(14)
f(11)=y(12)
f(9)=y(10)
f(7)=y(8)
f(5)=y(6)
f(3)=y(4)
f(1)=y(2)

f(20)=-20*(y(19)-y(17))-5*tanh(60*(y(20)+0.1))
f(18)=20*(y(19)-y(17))-20*(y(17)-y(15))-5*tanh(60*(y(18)+0.1))
f(16)=20*(y(17)-y(15))-20*(y(15)-y(13))-5*tanh(60*(y(16)+0.1))
f(14)=20*(y(15)-y(13))-20*(y(13)-y(11))-5*tanh(60*(y(14)+0.1))
f(12)=20*(y(13)-y(11))-20*(y(11)-y(9))-5*tanh(60*(y(12)+0.1))
f(10)=20*(y(11)-y(9))-20*(y(9)-y(7))-5*tanh(60*(y(10)+0.1))
f(8)=20*(y(9)-y(7))-20*(y(7)-y(5))-5*tanh(60*(y(8)+0.1))
f(6)=20*(y(7)-y(5))-20*(y(5)-y(3))-5*tanh(60*(y(6)+0.1))
f(4)=20*(y(5)-y(3))-20*(y(3)-y(1))-5*tanh(60*(y(4)+0.1))
%f(2)=20*(y(3)-y(1))-5*tanh(80*(y(2)+0.05))-200*sin(10*3.14*2*t)
f(2)=20*(y(3)-y(1))-5*tanh(60*(y(2)+0.1))+50*sin(10*3.14*t)




tspan=[0:0.1:40]
y0=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0]
[t,y]=ode15s('func4',tspan,y0)


问题和代码都在这里了啊，弹簧刚度都是200，物块质量10，摩阻系数0.5，右端施加一正弦力。不考虑图里面的速度，但是假的物块速度负方向的值超过0.1时候，摩擦力向右。求解出来，怎么是下面这样的呢？第一个物块的速度和第二个物块的速度差了那么多呢？后面，我查了，到第四个物块就没有振动了呢？到底咋回事啊？程序有问题么？求高手相助

lex9317

顺序有点问题哈，但是应该能看懂

lex9317

至于库伦摩擦光顺的问题，论坛里面有，rainboy曾经发过一篇这样的文章。

tammy

是不是F太小了？

lex9317

关键F的大小是可以让后面的物块震起来的，具体值我忘记了，但是应该大于前几个物块的摩擦力之和

happy

lex9317 发表于 2012-10-16 11:42

                               
登录/注册后可看大图

关键F的大小是可以让后面的物块震起来的，具体值我忘记了，但是应该大于前几个物块的摩擦力之和

模型中的V是什么含义？