how to represent a piecewise linear stiffness spring in simscape?
5 次查看(过去 30 天)
显示 更早的评论
I edited a .ssc file for simscape to represent a piecewise linear stiffness spring,but there is always something wrong, when the stiffness changed, the deformation of the spring changed suddenly too(diminishing), and here is the code, can someone help me? Thanks a lot!
component nlspring < foundation.mechanical.translational.branch
parameters
spr_rate1 = { 2500, 'N/m' }; % Spring rate
spr_rate2 = { 4000, 'N/m' }; % Spring rate
init_def = { 0, 'm' }; % Initial deformation
change_def = { 2.5e-3, 'm' }; % Stiffness changing deformation
end
variables
x = { 0, 'm '};
end
function setup
if spr_rate1 <= 0
pm_error('simscape:GreaterThanZero','spr_rate1' )
end
if spr_rate2 <= 0
pm_error('simscape:GreaterThanZero','spr_rate2' )
end
if change_def <= 0
pm_error('simscape:GreaterThanZero','change_def ' )
end
end
equations
if abs(x) < change_def
f == spr_rate1*x;
v == x.der;
else
f == spr_rate2*x;
v == x.der;
end
end
end
0 个评论
采纳的回答
Arnaud Miege
2011-3-24
I think you potentially have a discontinuity around +/-change_def. I would suggest changing the second equation to:
f == spr_rate2*(x-change_def)+spr_rate1*change_def;
to remove the discontinuity. I would also put the v == x.der; equation outside of the if...else statement. Finally, you need to add x = init_def; to the setup function.
HTH,
Arnaud
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Troubleshooting 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!