求解方法，变刚度

dw04116

静刚度k是x的函数，如何求解激励响应。我用线性加速度法求解，貌似有点问题，请各位大侠指点。

无水1324

用ode45可以不？
另外，有什么问题最好说明清楚，否则怎么帮你看呢

dw04116

回复 2 # 无水1324 的帖子

ode45,我感觉不容易控制。我的问题是，系统的刚度有两个值，当，位移<u0，时，k=k1，当位移>u0时，k=k2。我用newmark，积分编写的程序，明显有问题，一旦位移大于u0，振动就会出现非常诡异的变化。。所以，我猜想，应该是这种线性加速度法不适合这个问题。

无水1324

ode45,很好控制这个变化的参数，你可以把你的问题描述清楚！！

dw04116

回复 4 # 无水1324 的帖子

谢谢，我尝试一下。