基于MATLAB的船舶通风系统仿真
我们只能通过排风机的相似原理对风机进行计算，进而算出变容排风机的单位时间内所排出的气体体积。通过公式(2.1)公式(2.2)我们可以推到出变容风机容量与气压的关系，（由于进行对比的是同一个风机，环境温度相同，所以气体的密度，叶轮的半径都是相同的）推到公式如下：
PQ1?(1)2 (2.10) P2Q2通过公式(2.10)我们就可以得到变容排风机单位时间内的压强与气体量，进而通过公式就可以得到单位时间内的气体质量，从而得出变容排风机的数学模型。
2.3 各个密闭工作空间的数学建模
各个密闭空间的建模是根据图2.8所绘制的原理进行仿真建模的。
自然排风进入的气体气体质量的变化排出的气体
图2.7 密闭空间仿真原理图
由于从第一章我们可知，密闭空间排风分为自然排风和机械排风，所以影响空间气压的原因就由这两方面组成。可是，由于自然排风的不可控性，本课题没有对自然排风进行仿真研究，在建立模型时对自然排风只是给定了一个初始排风值。本课题主要是对机械排风进行仿真。有原理图可知，密闭空间的气压变化是由于空间内气体质量的变化所引起（在温度恒定的情况下）。所以，仿真的目的是计算出单位时间内密闭空间气体变化量，进而通过PV=m/M·RT这个气体方程变形式来计算出密闭空间的气压变化。
不同的密闭空间，数学模型不太一样。有排风机的密闭空间气压主要由送风机与排风机的气体质量差决定。无排风机的密闭空间气压由送风机送入的空气质量决定（自然排风影响小，本课题设置初始值影响，不过影响不大）。
12
基于MATLAB的船舶通风系统仿真
第3章 船舶通风系统的界面设计与仿真实现
3.1仿真软件介绍
MATLAB由一系列工具组成。这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。
MATLAB是一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序（M文件）后再一起运行。
MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。
本课题运用到的仿真软件是MATLAB中的GUI图形用户设计界面。在MATLAB强大的函数、矩阵运算支持下进行方针设计可以更加接近实际的情况。
此次仿真用到了GUI界面下的控件有Push Button、Static Text、Slider、Toggle Button、Edit Text等控件。Push Button控件是按键控件，是按下就弹起的按键；Static Text是静态的文本框，用来显示不需要更改的数据仿真结果。相对的Edit Text控件是用来显示可以根据实际情况进行更改的仿真数据。Toggle Button控件是按下后不弹起的按键，此次仿真用作为的开关按键。
3.2 仿真界面设计与实现
3.2.1通风管道网络界面设计
根据第一章介绍的某实船的通风网络管道布置，本课题运用MATLAB GUI中所提供的的控件进行模拟仿真。
（1）Push Button控件进行风机的仿真，按下按键表示风机开启，再按一次表示风机关闭，同时运用Static Text控件配合显示风机的状态。
（2）运用不同颜色的Static Text控件进行通风管道的布置，颜色区分管道是进风还是出风。
（3）各个密闭空间运用Push Button仿真，其目的是为了按下按键后可调出该空间的相关参数以及空间内气压的变化。
（4）运用Slider控件进行进风量，出风量大小的调节，即开关的仿真。
13
基于MATLAB的船舶通风系统仿真
（5）用Toggle Button控件设计系统的总开关，该仿真界面的所有运行程序都是在总开关的控制下进行的。
（6）Edit Text控件是设计温度，压强参数，目的是为了与实际的温度压强保持一致，进行更精确的仿真提供参数量。
（7）运用Static Text控件设置在甲板密闭空间以显示空间的压强变化。 通过以上各个控件的介绍以及每个控件所要模拟的单元，再通过第一章所给出的船舶通风系统的管道网络图，我们可以布置控件进行界面的设计，设计的界面如图3.1。各个控件初始化状态以运行时的界面状态为准。
图3.1 船舶通风管道网络仿真界面
3.2.2密闭空间界面设计
在MATLAB GUI 中运用参数传递的实现形式，把系统运行时间、每秒钟进入室内的气体体积、室内温度、室内压强、相对于外界的气压、进入室内的气体的总质量、每秒钟进入室内的气体质量、原机房气体质量、外界气体密度这些参数显示在密闭空间界面下，使得此界面能充分的体现出密闭空间的室内情况。同时，运用ActiveX Control 控件设计一个压力表，进行显示当前室内的相对气压，进而可以直观的展现室内气压状况，为控制气压提供依据。同时设计了退出按键，为方便的退出界面提供快捷方式。界面设计如图3.2。
14
基于MATLAB的船舶通风系统仿真
图3.2 船舶密闭空间仿真界面
3.2.3定容排风机系统模块的仿真实现
虽然每个风机系统模块是不同的，但是大体软件设计是相通的，所以，本文就针对两个相对复杂的系统进行软件设计实现思路的介绍。
首先先针对一个定容的排风机系统进行介绍。这里用了风机系统模块1，进行软件实现介绍。系统模块图如图2.2。仿真界面如图3.3。
图3.3 系统模块1仿真界面
15
基于MATLAB的船舶通风系统仿真
根据仿真图，实现定容排风机的仿真。仿真代码如下：
if(kC==0)
if (get(handles.TB7,'value')&anjian1C==1) if(anjian1C==1)
set(handles.e4,'string','3500');%设置定容排风机的容量
totalC1=str2double(r3C)-dC*str2double(r3C);%集控室剩余空气量 cC1 = dC*str2double(r3C);%集控室的空气量 totalC2=totalC1-totalC1*bC;%工作间剩余空气量 cC3=totalC1*bC;%工作间的空气量 totalC3=totalC2-totalC2*dCShang;
cC2=totalC2*dCShang;%主机舱上层的空气量 totalC4=totalC3-totalC3*dCDi;
cC4=totalC3*dCDi;%通舱底的空气量
totalC5=str2double(r3C)-dC*str2double(r3C)-totalC1*bC-totalC2*dCShang-totalC3*dCDi; %主机舱底层的空气量
set(handles.text74,'string',num2str(totalC5)); set(handles.text61,'string',num2str(cC2));
set(handles.text78,'string',num2str(cC4)); %分别显示出各个数值 Rv=str2double(r3C)*str2double(rC);
r11C=Rv/(str2double(r4)*str2double(r2));%每秒钟的进气量时间 r16=str2double(t5)*str2double(t6); r17=str2double(T0)+str2double(t5); r18=r16/r17;
R1=num2str(r18);%得到外界密度（可以不变） R2C=r11C*r18;%每秒钟的总质量
此代码是依据排风机的型号，我们先设置排风机的容量为‘3500’，并在e4控件中显示，同时设置r3C变量为风机的总容量，dC为集控室开关的开度，从而得到cC1，即空气的通气量。同理可以分别得到工作间、主机舱上层、通舱底以及主机舱底层的空气量并用text控件分别显示出来。然后进行固定数值的计算，如每秒钟的进气量以及外界的温度等。需要注意的是，这里涉及到了供风机的开关操作的功能设计，为了能使风机进行开关操作，本课题设计了一个全局变量global anjian1C，此变量与按键TB7共同管理此供风机的开关操作。
下面介绍此风机模块的两个密闭的工作空间，集控室与工作间的功能实现。集控室的代码实现如下：
if(k5C==1) %判断排风机开关是否开启
r20C=str2double(T1C)+str2double(t5);%室内温度34机房（变） r19C=R2C*dC;%计算进入室内的气体的质量
R2mC=num2str(r19C);%得到每秒钟进入的气体质量
MC=str2double('101000')*str2double(VC)*str2double(t8)/(str2double(t7)*r20C);
4机房的在一大气压下的空气质量（初始的室内空气质量）
m2C=str2double(m34C);%m34很重要！！！！！！ m2C=m2C+r19C-str2double(MzC);
16