车辆制动建模
本例将演示车辆制动建模。
汽车运动及其惯性由 Car 块描述。该程序块的左端口接收来自与制动机构相连的两对车轮的合力。右端口接收由同名子系统描述的空气阻力。
模型示意图：
定义加载和运行模型的函数：
In [ ]:
function start_model_engee()
try
engee.close("braking_distance", force=true) # закрытие модели
catch err # в случае, если нет модели, которую нужно закрыть и engee.close() не выполняется, то будет выполнена её загрузка после catch
m = engee.load("$(@__DIR__)/braking_distance.engee") # загрузка модели
end;
try
engee.run(m) # запуск модели
catch err # в случае, если модель не загружена и engee.run() не выполняется, то будут выполнены две нижние строки после catch
m = engee.load("$(@__DIR__)/braking_distance.engee") # загрузка модели
engee.run(m) # запуск модели
end
end
Out[0]:
运行模拟
In [ ]:
try
start_model_engee() # запуск симуляции с помощью специальной функции, реализованной выше
catch err
end;
从 simout 向变量写入速度、位移和力信号：
In [ ]:
t = simout["braking_distance/Скорость в км/ч"].time[:]
distance = simout["braking_distance/Перемещение в км"].value[:]
speed = simout["braking_distance/Скорость в км/ч"].value[:]
forces_left = simout["braking_distance/Сила со стороны трансмиссии"].value[:]
forces_right = simout["braking_distance/Сила сопротивления воздуха"].value[:];
模拟结果可视化
In [ ]:
using Plots
In [ ]:
plot(t, forces_left, label="Усилие от трансмиссии, Н", linewidth=2)
plot!(t, forces_right, label="Сопротивление воздуха, Н", linewidth=2)
plot!(t, (forces_left+forces_right), label="Сумма сил, Н", linewidth=2, title="Силы, действующие на автомобиль")
Out[0]:
In [ ]:
plot(t, speed, label="Скорость, км/ч", linewidth=2)
plot!(t, distance*1000, label="Перемещение, м", linewidth=2)
Out[0]:
定义一个计算汽车制动距离的函数，并将其应用于位移数据：
In [ ]:
function braking_distance(distance, t)
index = findfirst(x -> x < 0.00001, distance)
if index != nothing
return println("Тормозной путь: ", round(maximum(distance)*1000, digits=2), " м")
end
end
braking_distance(distance, t)
结论
本例演示了车辆制动建模。使用函数braking_distance 计算了一辆初始速度为 100 km/h 的汽车的制动距离。