fds手册
fds是常用的火灾动力学软件,在此对部分用法做基本介绍,望对各位学子有帮助。
一、语法特点
fds代码的的特点为 & 开头+ 函数名 + 变量 ,语句以 “/” 结尾,且字母均为大写。
&MASH XB = 0 6 0 4.5 0 2 COLOR = "RED" /
二、常用函数
(一)&MASH 函数
&MASH函数是建立一个基本区域,后续的建模全部基于此。XB代表其大小,COLOR代表其颜色,也可不指定颜色。
&MASH XB = 0 6 0 4.5 0 2 COLOR = "RED" /
(二)&VENT 函数
&VENT函数表示一个平面,可通过其的ID参数改变其特性。
1.打开空白平面
&VENT函数可在平面上打开一个空白区域,此情况常用在&MASH函数创建的区域上开口。
&VENT XB = 0 1 2 3 0 0 COLOR = "RED" ID = "OPEN" /
2.设置火源
&VENT函数可设置火源,只能是平面火源,但要设置火灾模型,HRRPUA是指每平方米十分的能量,FUEL表示燃料类型,SOOT_YIELD表示燃烧中产生碳烟的比例。
&SURF ID = 'FIRE' HRRPUA = 1000 /
&VENT XB = 0 1 2 3 0 0 SURF_ID = 'FIRE' /
&REAC FUEL = 'PROPANE' SOOT_YIELD = 0.05 /
(三)&OBST函数
&OBST函数表述设置障碍物,其可用来设置隔墙,也可用于设置火源,相较于&VENT函数,&OBST可设置立体火源。
1.设置隔墙
在用&OBST设置隔墙时其的大小不能超过&MASH的区域。
&OBST XB = 2.9 3.1 0 4.5 0 2.8 COLOR = 'GRAY' /
2.设置火源
在用&OBST设置火源时,基本操作于&VENT一致,注意Z坐标就行。
&SURF ID = 'FIRE' HRRPUA = 1000 /
&OBST XB = 0 1 2 3 1 2 SURF_ID = 'FIRE' COLOR = 'RED' /
&REAC FUEL = 'PROPANE' SOOT_YIELD = 0.05 /
(四)&TIME函数
&TIME函数用于项目中的时间类型,列如控制计算时间,其单位为秒。
&TIME T_END = 6 /
(五)&HOLE函数
如其字面意思,表示在物体上挖洞,其大小不能超过被挖洞物体的尺寸,大多数情况用于除&MASH区域的范围。
&OBST XB = 2.9 3.1 0 4.5 0 2.8 COLOR = 'GRAY' ID = '墙1'/
&HOLE XB = 2.9 3.1 1 1.6 0 2.8 ID = '墙1' /
(六)&HADE函数
此函数可设置输出的一系列文件的名称,设置之后开头都会以此。
&HADE CHID = 'first' /
(七)&TAIL函数
此函数没有任何参数,表示代码的结束。
&TAIL /
(八)&MULT函数
&MULT用于创造重复物体,是基于已存在的物体,DX,DY,DZ是指相邻物体同一点的三个坐标距离,而N_UPPER表示共创建的个数。
&MULT ID = '楼梯' DX = 0 DY = 0.3 DZ = 0.2 N_UPPER = 7 /
&OBST XB = 0 0.7 0 0.3 0 0.2 MULT_ID = '楼梯' COLOR = 'GRAY' /
(九)&DEVC函数
&DEVC函数用于测某一点随时间变化的数据,XYZ表示该点的坐标,QUANTITY表示测量该点的数据类型,ID则是为该数据列命名,运算结束后会输出一张csv表,但要注意点的坐标尽量不要出现在网格交线上。
&DEVC XYZ = 3.1 2.1 1.1 QUANTITY = 'TEMPERATURE' ID = 'T1' /
若是测该点的气体体积比,应指明气体类型。
&DEVC XYZ = 3.1 2.1 1.1 QUANTITY = 'VOLUME FRACTION' SPEC_ID = 'CARBEN DIOXIDE' ID = 'T1' /
(十)&SLCF函数
&SLCF此函数用于测量垂直于坐标轴平面的数据,PBY表示与Y轴垂直,距原点距离。
&SLCF PBY = 0.4 QUANTITY = 'TEMPERATURE' /
(十一)&ISOF函数
&ISOF是等值面命令,凸显数据相等的面,在下述代码中表示凸显温度为30,40,500度的平面,其中VALUE最多3个。
&ISOF QUANTITY = 'TEMPERATURE' VALUE = 30 40 500 /
(十二)&MATL函数
用于设置物体的性质,如密度,比热容,导热性等,但用此函数时要配合边界函数使用,并加上厚度,再加上物体函数.
&MATL ID = '钢' DENSITY = 7850 CONDUCTIVITY = 45 SPECIFIC_HEAT = 0.5 /
&SURF ID = '钢材' MATL_ID = '钢' THICKNESS = 0.1 /
&OBST XB = 0.1 1.1 0.1 1.1 0.2 1 SURF_ID = '钢材' COLOR = 'GRAY' /
(十三)&PROF命令
测物体内部的一条线的导热数据,其长度基于THICKNESS的值计算,且计算基于差分法,结果中的prof类的csv文件中,各列的意义,先是时间再到参考点的坐标,然后是其余差分法取的点到参考点的距离,然后是温度。
&PROF XYZ = 0.1 0.55 0.45 IOR = -1 QUANTITY = 'TEMPERATURE' /
(十四)&BNDF命令
此命令用于画物体相关数据的动画,在对实体物时,用QUANTITY = ‘WALL TEMPERATURE’。
&BNDF QUANTITY = 'WALL TEMPERATURE' /
(十五)&MISC函数
里面储存了函数的参数,可以对其进行修改。
&MISC BNDF_DEFAULT = .FALSE. /
(十六)&RMAP函数
&RMAP函数可用于定义两个变量的关系,且必须从起始值开始。
&MATL ID = '钢' DENSITY = 7850 CONDUCTIVITY = 45 SPECIFIC_HEAT_RMAP = 'RD' /
&RMAP ID = 'RD' T = 20 F = 0.4 /
&RMAP ID = 'RD' T = 100 F = 0.5 /
&RMAP ID = 'RD' T = 200 F = 0.6 /
三、常见操作
(一)t方火
火灾可根据燃烧速度分为:超快速火、快速火、中速火、低速火,其释放热的公式可为:Q = at^2。因为快速火最常见,在此介绍快速火:其系数为0.04689,其方程为Q = 0.04689t^2。
在这假设生成3MW的火,所以 3000 = 0.04689*t^2,得出 t = 253s。
&SURF ID = 'FIRE' HRRPUA = 3000 TAU_Q = -253 /
&VENT XB = 3 4 3 4 0 0 SURF_ID = 'FIRE' COLOR = 'RED' /
&REAC FUEL = 'PROPANE' SOOT_YIELD = 0.05 /
(二)导热系数
对一个物体来说,密度越小,比热容越小,升温越快,导热系数越大,物体的导热是按THICKNESS的值计算的,而对流的导热是按&OBST中的XB参数计算的。
利用&MATL函数设置导热参数
&MATL ID = '钢' DENSITY = 7850 CONDUCTIVITY = 45 SPECIFIC_HEAT = 0.5 /
&SURF ID = '钢材' MATL_ID = '钢' THICKNESS = 0.1 /
&OBST XB = 0.1 1.1 0.1 1.1 0.2 1 SURF_ID = '钢材' COLOR = 'GRAY' /
(三)固定边界温度
将某一个物体或者某一平面的表面温度设为固定值.
&SURF ID = 'HOT' TMP_FRONT = 800 /
&VENT XB = 0 0 0 1 0 1 SURF_ID = 'HOT' COLOR = 'RED' /
(四)物体表面温度
在此选用&DEVC函数输出物体表面温度,XYZ为坐标,必须在物体表面上,IOR代表相对于物体的位置,利用物体中心建立坐标轴,在X轴正半轴则为1反之为-1,Y轴正半轴为2,反之为-2,Z轴正半轴为3,反之为-3,由于被测的为固体,所以参数QUANTITY = ‘WALL TEMPERATURE’。
&DEVC XYZ = 0.1 0.55 0.45 QUANTITY = 'WALL TEMPERATURE' IOR = -1 ID = 'TT' /
(五)物体内部温度
和物体表面温度的测定方法差不多,XYZ一样在表面,就是QUANTITY = ‘INSIDE WALL TEMPERATURE’,然后DEPTH表示深度。
&DEVC XYZ = 0.1 0.55 0.45 QUANTITY = 'INSIDE WALL TEMPERATURE' DEPTH = 0.2 IOR = -1 ID = '内部温度' /
(六)物体表面温度动画
由于在使用&BNDF函数时会显示所有物体的动画,不方便观察,所以先关闭所有物体动画的显示,再打开目标物体的动画显示,其中&MISC为函数参数选项。
&OBST XB = 0.1 0.3 0.5 0.6 0.4 0.5 SURF_ID = '钢材' COLOR = 'GRAY' BNDF_OBST = .TRUE./
&BNDF QUANTITY = 'WALL TEMPERATURE' /
&MISC BNDF_DEFAULT = .FALSE. /
(七)质量损失比率
表示燃烧过程,可以用可燃物质量损失的比率,将热释放效率换为质量损失的比例就行。
&SURF ID = 'FIRE' MLRPUA = 0.05 /
(八)比热设置为与时间关联的函数
T代表温度,F代表比率。
&MATL ID = '钢' DENSITY = 7850 CONDUCTIVITY = 45 SPECIFIC_HEAT_RMAP = 'RD' /
&RMAP ID = 'RD' T = 20 F = 0.4 /
&RMAP ID = 'RD' T = 100 F = 0.5 /
&RMAP ID = 'RD' T = 200 F = 0.6 /
(九)热释放速率设置为时间的函数
这样操作时徐注意,HRRPUA需设置为目标的最大值,而&RMAP也应从最初值开始,而此时是关于时间的函数,所以都为0,而F是关于HRRPUA的比例。
&SURF ID = 'FIRE' HRRPUA = 1000 RMAP_Q = 'RT' /
&RMAP ID = 'RT' T = 0 F = 0.0 /
&RMAP ID = 'RT' T = 20 F = 0.5 /
&RMAP ID = 'RT' T = 40 F = 1.0 /