手把手教你学simulink--基于Simulink的建筑结构自适应振动控制仿真实例

小蘑菇二号

已于 2025-02-05 10:45:13 修改

阅读量696

点赞数 4

文章标签： matlab simulink

于 2025-02-05 01:03:09 首次发布

本文链接：https://blog.csdn.net/MHD0815/article/details/145367364

版权

步骤2：在Simulink中实现自适应控制器

基于Simulink的建筑结构自适应振动控制仿真实例，我们将从系统建模、控制器设计、Simulink实现等方面进行详细说明。假设我们要设计一个用于建筑结构（如高层建筑或桥梁）的自适应振动控制系统，并使用Simulink进行仿真验证。

项目背景

系统描述

建筑物在受到外部激励（如地震、风荷载等）时会发生振动，这些振动可能会导致结构损坏或影响建筑物内人员的安全和舒适度。因此，设计有效的振动控制系统对于提高建筑物的安全性和耐久性至关重要。

建筑结构的振动控制可以通过安装主动或半主动控制系统来实现，例如主动质量阻尼器（Active Mass Damper, AMD）或磁流变阻尼器（Magnetorheological Damper, MR）。这些控制系统通过实时调节阻尼力来抑制结构的振动。

建筑结构的动力学模型可以表示为：

单自由度（SDOF）模型：

mx¨(t)+cx˙(t)+kx(t)=Fext(t)−Fctrl(t)mx¨(t)+cx˙(t)+kx(t)=Fext(t)−Fctrl(t)

其中：
- x(t)x(t) 是结构的位移。
- mm 是结构的质量。
- cc 是结构的阻尼系数。
- kk 是结构的刚度系数。
- Fext(t)Fext(t) 是外部激励力（如地震或风荷载）。
- Fctrl(t)Fctrl(t) 是控制输入力（如AMD或MR阻尼器提供的力）。
振动控制目标：控制目标是通过适当的控制策略使建筑结构能够在外部激励下保持较小的振动幅度，从而提高建筑物的安全性和舒适性。

我们的目标是设计一个自适应控制器，使建筑结构能够在不确定环境下实现精确的振动控制，并有效抑制振动。

步骤1：设计自适应控制器

1.1 定义系统参数

首先，在MATLAB命令窗口中定义系统的参数和状态空间模型。

matlab

深色版本

% 系统参数
m = 1000; % 结构质量 (kg)
c = 50; % 阻尼系数 (Ns/m)
k = 1e6; % 刚度系数 (N/m)

% 外部激励参数
F_ext_amplitude = 1000; % 激励力幅值 (N)
F_ext_frequency = 1; % 激励力频率 (Hz)

% 初始条件
x0 = 0; % 初始位移 (m)
v0 = 0; % 初始速度 (m/s)

1.2 设计自适应控制器

自适应控制器的设计包括选择合适的控制律和自适应机制。

假设我们采用基于滑模控制（Sliding Mode Control, SMC）的自适应控制策略：

滑模控制（SMC）：通过引入滑模面来设计控制器，以实现对期望位移的精确跟踪并抑制振动。

控制律可以表示为：

Fctrl(t)=Kpex(t)+Ki∫ex(t)dt+Kddex(t)dt+λsign(s)Fctrl(t)=Kpex(t)+Ki∫ex(t)dt+Kddtdex(t)+λsign(s)

其中：

ex(t)=xd(t)−x(t)ex(t)=xd(t)−x(t) 是位移误差。
s=e˙x(t)+λ1ex(t)s=e˙x(t)+λ1ex(t) 是滑模面。
λ1λ1 和 λλ 是控制器增益。
sign(s)sign(s) 是符号函数。

步骤2：在Simulink中实现自适应控制器

2.1 创建新的Simulink模型

在MATLAB命令窗口中创建一个新的Simulink模型。

matlab

深色版本

modelName = 'Building_Structure_Adaptive_Vibration_Control';
new_system(modelName);
open_system(modelName);

2.2 添加被控对象模型

在Simulink中添加被控对象的状态方程模型。

建筑结构动力学模型

我们可以使用State-Space块来表示建筑结构的动力学模型。

matlab

深色版本

% 添加State-Space块表示建筑结构动力学
add_block('simulink/Continuous/State-Space', [modelName '/Structure Dynamics']);
set_param([modelName '/Structure Dynamics'], 'A', '[0 1; -k/m -c/m]');
set_param([modelName '/Structure Dynamics'], 'B', '[0; 1/m]');
set_param([modelName '/Structure Dynamics'], 'C', '[1 0]');
set_param([modelName '/Structure Dynamics'], 'D', '[0]');

外部激励模型

我们可以使用Sine Wave块来模拟外部激励力（如地震或风荷载）。

matlab

深色版本

% 使用Sine Wave块模拟外部激励力
add_block('simulink/Sources/Sine Wave', [modelName '/External Excitation']);
set_param([modelName '/External Excitation'], 'Amplitude', 'F_ext_amplitude');
set_param([modelName '/External Excitation'], 'Frequency', 'F_ext_frequency');
set_param([modelName '/External Excitation'], 'Phase', '0');

2.3 实现自适应控制器

使用MATLAB Function块实现自适应控制器。

matlab

深色版本

% 自适应控制器的MATLAB函数
function F_ctrl = adaptive_vibration_control(e_x, de_x, integral_e_x, lambda1, lambda, Kp, Ki, Kd)
    % 计算滑模面
    s = de_x + lambda1 * e_x;
    
    % 控制律
    F_ctrl = Kp * e_x + Ki * integral_e_x + Kd * de_x + lambda * sign(s);
end

% 添加MATLAB Function块实现自适应控制器
add_block('simulink/User-Defined Functions/MATLAB Function', [modelName '/Adaptive Controller']);
set_param([modelName '/Adaptive Controller'], 'Function', 'adaptive_vibration_control');
set_param([modelName '/Adaptive Controller'], 'Inputs', 'e_x, de_x, integral_e_x, lambda1, lambda, Kp, Ki, Kd');
set_param([modelName '/Adaptive Controller'], 'Outputs', 'F_ctrl');

2.4 添加参考轨迹生成器

使用Step块生成参考位移信号。

matlab

深色版本

% 使用Step块生成参考位移信号
add_block('simulink/Sources/Step', [modelName '/Reference Displacement']);
set_param([modelName '/Reference Displacement'], 'Time', '5'); % 在5秒时发生位移变化
set_param([modelName '/Reference Displacement'], 'Before', '0'); % 变化前的位移 (m)
set_param([modelName '/Reference Displacement'], 'After', '0.01'); % 变化后的位移 (m)

2.5 添加积分器和微分器

使用Integrator块计算位移偏差的积分，使用Derivative块计算位移偏差的微分。

matlab

深色版本

% 添加Integrator块计算位移偏差的积分
add_block('simulink/Continuous/Integrator', [modelName '/Integral of Displacement Error']);

% 添加Derivative块计算位移偏差的微分
add_block('simulink/Continuous/Derivative', [modelName '/Derivative of Displacement Error']);

2.6 连接反馈回路

连接各个模块形成闭环控制系统。

matlab

深色版本

% 连接Structure Dynamics的输出到Adaptive Controller的输入端
connect_blocks(modelName, {'/Structure Dynamics', '/Adaptive Controller'});

% 连接Reference Displacement到Adaptive Controller的输入端
connect_blocks(modelName, {'/Reference Displacement', '/Adaptive Controller'});

% 连接Integral of Displacement Error到Adaptive Controller的输入端
connect_blocks(modelName, {'/Integral of Displacement Error', '/Adaptive Controller'});

% 连接Derivative of Displacement Error到Adaptive Controller的输入端
connect_blocks(modelName, {'/Derivative of Displacement Error', '/Adaptive Controller'});

% 连接Adaptive Controller的输出到Structure Dynamics的输入端
add_line(modelName, '/Adaptive Controller/1', '/Structure Dynamics/2', 'autorouting', 'on');

% 如果有外部激励模型，连接External Excitation到Structure Dynamics的输入端
if exist([modelName '/External Excitation'], 'file')
    add_line(modelName, '/External Excitation/1', '/Structure Dynamics/2', 'autorouting', 'on');
end

2.7 添加控制增益

使用Constant块设置比例、积分和微分增益。

matlab

深色版本

% 添加Constant块设置比例增益
add_block('simulink/Sources/Constant', [modelName '/Kp']);
set_param([modelName '/Kp'], 'Value', '1');

% 添加Constant块设置积分增益
add_block('simulink/Sources/Constant', [modelName '/Ki']);
set_param([modelName '/Ki'], 'Value', '0.1');

% 添加Constant块设置微分增益
add_block('simulink/Sources/Constant', [modelName '/Kd']);
set_param([modelName '/Kd'], 'Value', '0.05');

% 添加Constant块设置滑模面增益
add_block('simulink/Sources/Constant', [modelName '/Lambda1']);
set_param([modelName '/Lambda1'], 'Value', '1');

% 添加Constant块设置符号函数增益
add_block('simulink/Sources/Constant', [modelName '/Lambda']);
set_param([modelName '/Lambda'], 'Value', '0.5');

2.8 添加输出显示

使用Scope块显示系统输出。

matlab

深色版本

% 添加Scope块显示系统输出
add_block('simulink/Sinks/Scope', [modelName '/Scope']);
connect_blocks(modelName, {'/Structure Dynamics', '/Scope'});
connect_blocks(modelName, {'/Reference Displacement', '/Scope'});
connect_blocks(modelName, {'/External Excitation', '/Scope'});

步骤3：运行仿真并分析结果

3.1 设置仿真参数

在上述代码中，已经设置了仿真时间为20秒，并保存了模型。

matlab

深色版本

% 设置仿真时间
set_param(modelName, 'StopTime', '20');

% 保存模型
save_system(modelName);

3.2 运行仿真

打开Simulink模型，点击工具栏中的“Run”按钮运行仿真。

matlab

深色版本

% 打开Simulink模型并运行仿真
open_system(modelName);
sim(modelName);

3.3 分析仿真结果

在Simulink模型中，打开Scope块查看系统响应曲线。

matlab

深色版本

% 打开Scope窗口
open_system([modelName '/Scope']);

分析仿真结果

在Scope窗口中，你可以看到以下主要曲线：

位移（Displacement）：展示建筑结构的位移随时间的变化情况。
参考位移（Reference Displacement）：展示预设的目标位移随时间的变化情况。
外部激励力（External Excitation）：展示外部激励力的变化情况。

通过这些曲线，你可以评估控制器的性能：

位移跟踪性能：检查建筑结构的位移是否能够准确跟随预设的目标位移。
振动抑制性能：观察位移响应曲线是否存在明显的振动，并评估振动的幅度和频率。

示例分析

假设你希望建筑结构在5秒时将位移从0米调整为0.01米，并观察位移响应。通过观察Scope窗口中的曲线，可以得出以下结论：

位移跟踪性能：如果位移曲线能够在指定的时间点快速且平滑地达到目标值，并保持稳定，则表明控制器具有良好的跟踪性能。
振动抑制性能：如果位移曲线在调整后没有明显的振荡，并且最终趋于目标位移，则表明控制器具有良好的振动抑制性能。