具有平移摩擦的机械系统
这个例子显示了一个质量连接到一个弹簧和一个粘性阻尼器. 质量由通过摩擦元件的理想速度源驱动。 速度源的运动轮廓以这样的方式选择，即当绘制质量位移对源提供的位移的依赖性时，获得典型的滞后曲线。
模型图:
定义加载和运行模型的函数:
In [ ]:
function start_model_engee()
try
engee.close("mechanical_system_translational_friction", force=true) # закрытие модели
catch err # в случае, если нет модели, которую нужно закрыть и engee.close() не выполняется, то будет выполнена её загрузка после catch
m = engee.load("$(@__DIR__)/mechanical_system_translational_friction.engee") # загрузка модели
end;
try
engee.run(m) # запуск модели
catch err # в случае, если модель не загружена и engee.run() не выполняется, то будут выполнены две нижние строки после catch
m = engee.load("$(@__DIR__)/mechanical_system_translational_friction.engee") # загрузка модели
engee.run(m) # запуск модели
end
end
Out[0]:
In [ ]:
Pkg.add("Plots");
启动模型:
In [ ]:
try
start_model_engee() # запуск симуляции с помощью специальной функции, реализованной выше
catch err
end;
从simout变量输出仿真结果:
In [ ]:
res = collect(simout)
Out[0]:
将结果写入变量:
In [ ]:
P1 = collect(res[1]); # Перемещение массы
P2 = collect(res[4]); # Перемещение источника скорости
v1 = collect(res[2]); # Скорость источника
模拟结果:
In [ ]:
using Plots
plot(P1[:,1], P1[:,2], linewidth=3, xlabel= "Время, с", ylabel= "Перемещение, м", legend=:bottomright, label="Масса")
plot!(P2[:,1], P2[:,2], linewidth=3, label="Источник скорости")
Out[0]:
In [ ]:
plot(v1[:,1], v1[:,2], linewidth=3, xlabel= "Время, с", ylabel= "Скорость, м/с", legend=:bottomright, label="Источник скорости")
Out[0]:
In [ ]:
plot(P2[:,2], P1[:,2], linewidth=3, xlabel="Перемещение источника скорости, м", ylabel="Перемещение массы, м", title="Кривая гистерезиса перемещения", legend=false)
Out[0]: