Engee 文档
Notebook

具有刚性平移运动限制器的机械系统

这个例子显示了两个质量通过一个刚性限位器连接. 质量1由理想的速度源驱动。 当输入速度的方向改变时，质量块2不会移动，直到质量块1到达由平移刚性限制器建模的间隙的另一端。

模型图:

mechanical_system_translational_hardstop--1733299118648.png

定义加载和运行模型的函数:

In [ ]: 
function start_model_engee()
    try
        engee.close("mechanical_system_translational_hardstop", force=true) # закрытие модели 
        catch err # в случае, если нет модели, которую нужно закрыть и engee.close() не выполняется, то будет выполнена её загрузка после catch
            m = engee.load("$(@__DIR__)/mechanical_system_translational_hardstop.engee") # загрузка модели
        end;

    try
        engee.run(m) # запуск модели
        catch err # в случае, если модель не загружена и engee.run() не выполняется, то будут выполнены две нижние строки после catch
            m = engee.load("$(@__DIR__)/mechanical_system_translational_hardstop.engee") # загрузка модели
            engee.run(m) # запуск модели
        end
end
Out[0]:
start_model_engee (generic function with 1 method)

启动模型:

In [ ]: 
try
    start_model_engee() # запуск симуляции с помощью специальной функции, реализованной выше
    catch err
    end;

从simout变量输出仿真结果:

In [ ]: 
res = collect(simout)
Out[0]:
8-element Vector{WorkspaceArray}:
 WorkspaceArray{Float64}("Mass 2.flange.v")
 WorkspaceArray{Float64}("Mass 2.F")
 WorkspaceArray{Float64}("Mass 2.v")
 WorkspaceArray{Float64}("Mass 1.flange.v")
 WorkspaceArray{Float64}("mechanical_system_translational_hardstop/SysOutput")
 WorkspaceArray{Float64}("Mass 1.F")
 WorkspaceArray{Float64}("mechanical_system_translational_hardstop/Ideal Translational Motion Sensor.2")
 WorkspaceArray{Float64}("Mass 1.v")

将结果写入变量:

In [ ]: 
P1 = collect(res[7]); # Перемещение массы 1 "mechanical_system_translational_hardstop/Ideal Translational Motion Sensor.2")
P2 = collect(res[5]); # Перемещение массы 2 ("mechanical_system_translational_hardstop/SysOutput")
v1 = collect(res[8]); # Скорость массы 1
v2 = collect(res[3]); # Скорость массы 2

模拟结果:

In [ ]: 
using Plots
plot(P1[:,1], P1[:,2], linewidth=3, xlabel= "Время, с", ylabel= "Перемещение, м", legend=:bottomright, label="Масса 1")
plot!(P2[:,1], P2[:,2], linewidth=3, label="Масса 2")
Out[0]:
In [ ]: 
using Plots
plot(v1[:,1], v1[:,2], linewidth=3, xlabel= "Время, с", ylabel= "Скорость, м/с", legend=:bottomright, label="Масса 1")
plot!(v2[:,1], v2[:,2], linewidth=3, label="Масса 2")
Out[0]:

示例中使用的块