Solution of spring mass equation with cubic stiffness

Question

Murat Yetkin 2018-10-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/423183-solution-of-spring-mass-equation-with-cubic-stiffness

重新打开： Walter Roberson 2018-12-22

I have below equation to solve with a given time series of f(t). Is there any idea how to do it in Matlab with ode solvers?

m=1; c=1; k=10; k2=20; k3=20;
my''+cy'+ky+k2*y^2+k3*y^3=f(t)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2018-10-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/423183-solution-of-spring-mass-equation-with-cubic-stiffness#answer_340673

在 MATLAB Online 中打开

function main
tf = ...;  % times for f-vector
f = ...;   % f vector
tspan=[0 10000];
y0=[0;1];
[t,y]=ode45(@(t,y) myode1(t,y,tf,f),tspan,y0);
end
function dy = myode1(t,y,tf,f)
m=1; c=1; k=10; k2=20; k3=20;
f_actual=interp1(tf,f,t);
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=(f_actual-c*y(2)-k*y(1)-k2*y(1)^2-k3*y(1)^3)/m;
end