一种基于参数化模型的整体结构耐火验算方法与流程

1.本发明涉及钢结构抗火设计技术领域，尤其是一种基于区域升温分析的整体结构耐火极限计算方法。


背景技术：

2.建筑火灾是人类生产生活中危害最严重的灾害之一，对人类生命财产安全构成了直接的危害，而发生的频率则居各灾害之首。由于钢材的耐火性较差，高温下会丧失大部分强度，因此钢结构在火灾下的安全性常常会受到威胁，甚至会造成结构倒塌并形成二次灾害。为了达到减少火灾直接经济损失、避免或减少人员伤亡的建筑防火设计目标，结构防火设计需要避免火灾中结构局部倒塌造成灭火及人员疏散困难，以及避免结构在火灾中整体倒塌造成人员伤亡。
3.《建筑钢结构防火技术规范》gb51249-2017给出了结构防火设计的具体实施办法，根据验算对象和层次的不同，分为基于整体结构耐火验算和基于构件耐火验算的防火设计方法。其中，预应力钢结构和跨度不小于120m的大跨度建筑中的钢结构，应采用基于整体结构耐火验算的防火设计方法。对于基于构件耐火验算的防火设计方法，gb51249中给出了明确的构件升温计算方法以及耐火验算方法，而针对整体结构耐火验算的方法则仅限于总流程和规定，详细具体的指导则较少。
4.现有技术的整体结构耐火验算存在着以下问题：
5.1)建筑方案的一般改动较频繁，尤其是形态复杂的公共建筑，而火灾场景的网格设置需要及时适应建筑形态的变化，y、z轴网格划分的数量需采用2k
×
3m
×
5n来提升fds运行速度，同时火源位置还需要采用节点与整体网格一致的局部精细网格来提升模拟精度。
6.2)火灾场景得到的空气升温情况，需要为构件的升温提供数据；而构件的升温，除受所在位置空气升温情况的影响外，还受到防火层厚度、截面尺寸、材料参数的影响，若采用手动计算和设置，简化模型会存在偏差较大的问题，而精细程度较高的模型则会耗费大量时间和人力。
7.3)由于火灾发生的不确定性，一般建筑会有多个较不利的火灾场景，火灾场景的变化也会使上述模型调整的工作量成倍增加，严重影响设计效率。


技术实现要素：

8.本发明的目的是针对现有技术的不足而设计的一种基于参数化模型的整体结构耐火验算方法，采用建筑表皮模型和结构分析模型的方法，形成参数控制的火灾场景分析模型和整体结构耐火验算模型，自动进行火灾场景分析模型的建立、空气升温曲线的提取、各构件升温曲线的求解以及整体结构耐火验算模型的建立，为判断整体结构在设定火灾场景下的安全性评估提供依据，方法简便，可快速求得整体结构的耐火极限，不但解决了手动计算和设置偏差大的问题，而且大大降低了人力和大量时间的耗费，进一步提高设计效率，尤其适用于高层钢结构的耐火验算和安全性能分析。
9.本发明的目的是这样实现的：一种基于参数化模型的整体结构耐火验算方法，其特征在于整体结构耐火验算具体包括以下步骤：
10.a、形成建筑结构一体化信息模型
11.将结构分析模型读入到参数化平台rhino grasshopper中，提取其中的节点信息、单元信息、截面信息和材料信息。同时，在rhino软件中导入火灾场景分析需要用到的建筑模型，包括建筑表皮、楼屋面、门窗的曲面，并使建筑与结构模型的空间相对位置关系与实际情况保持一致。
12.b、生成参数化的火灾场景分析模型
13.根据建筑信息生成用于火灾场景分析的fds文本文件，具体为根据建筑总体尺寸信息形成整体网格的mesh命令，根据建筑表皮、楼屋面曲面信息等形成obst命令。然后根据防火分区内可能出现的最不利火灾场景，确定火源功率、火灾增长类型和火床面积信息，形成火源处局部精细网格的mesh命令、火源位置的obst命令、热边界条件的surf命令、以及火灾时间的time命令。
14.根据结构分析模型中提取的节点信息、单元信息，形成获取构件所在位置空气升温的devc命令，并在命令中设定带结构单元编号的id信息。
15.c、火灾场景模拟
16.运行生成的fds文件，进行火灾场景分析，并完成求解。
17.d、计算各构件升温曲线
18.设定不同截面构件的防火层厚度，再根据结构分析模型中提取的单元信息和截面信息、fds求解得到的每根构件位置处的空气升温曲线，按《建筑钢结构防火技术规范》gb51249中构件的升温计算公式，求解达到整体结构耐火极限时各构件的温度。
19.e、整体结构耐火验算模型的建立及求解
20.根据结构分析模型建立整体结构耐火验算模型，并将求解得到的各构件温度作为荷载输入到验算模型中，同时采用与温度相关的材料进行求解。
21.f、整体结构耐火验算及防火保护措施调整
22.根据结构响应情况判断是否达到耐火承载力极限状态，判断依据为：结构在高温下产生足够的塑性铰形成可变机构，或结构整体丧失稳定。若结构达到耐火承载力极限状态，则表明设置的防火层厚度不足，调整防火层厚度，重复第4～6步，直至整体结构不会达到耐火承载力极限状态，由此得到设定火灾场景下所需的防火层厚度，形成防火保护设计方案。
23.本发明与现有技术相比具有以下显著的技术进步和优点：
24.1)采用建筑结构一体化信息模型，并结合参数化控制，能够准确、自动地生成火灾场景命令文本，使火灾场景分析时fds能自动输出每根构件处空气的升温情况，用于进行每根结构构件的升温计算。
25.2)通过参数控制，能够根据fds输出的空气温度数据文件、结构的截面材料信息以及防火层设定，自动计算每根构件的升温情况，并形成分析所需的命令文本，开展精细的整体结构耐火验算。
26.3)采用参数化平台进行插件开发，调整体网格尺寸、建筑信息(建筑表皮、火床位置及大小、门窗情况)、火源设定(火源功率、火灾增长类型、求解时间)、整体结构耐火极限、
防火保护材料类型及设计指标等设定，火灾场景模型、整体结构耐火验算模型会在参数驱动下自动进行调整，可减少人力，提高工作效率。
附图说明
27.图1为本发明流程示意图；
28.图2为建筑结构一体化信息模型示意图；
29.图3为实施例1的建筑壁面障碍物模型示意图；
30.图4为实施例1的火灾场景模拟示意图；
31.图5为实施例1各构件达到耐火极限时的升温示意图；
32.图6为实施例1整体结构耐火验算模型的温度分布示意图；
33.图7实施例1整体结构在设定火灾场景下的构件应变示意图。
具体实施方式
34.下面以具体实施例对本发明做进一步的阐述。
35.实施例1
36.参阅附图1，整体结构耐火验算按下述步骤进行：
37.1、形成建筑结构一体化信息模型
38.在midas/gen中建立结构模型，导出mgt文本文件。在参数化平台rhino grasshopper中读入mgt文件，提取其中的节点信息*node、单元信息*element、截面信息*section和材料信息*matreial。
39.参阅附图2，在rhino软件中导入火灾场景分析需要用到的建筑模型，包括关注的防火分区内的建筑表皮、楼屋面、门窗、火床的曲面，并使建筑与结构模型的空间相对位置关系与实际情况保持一致。
40.2、生成参数化的火灾场景分析模型
41.根据建筑信息生成用于火灾场景分析的fds文本文件，具体为：
42.1)提取所有建筑表皮边缘及节点信息
43.确定建筑总体尺寸为：x向27.15m～69.35m，y向34.03m～59.42m，z向为-0.01m～16.9m。设定整体网格尺寸为0.4m，保留小数点后1位，采用fds推荐的单轴2k
×
3m
×
5n的网格数量划分网格，形成mesh命令如下：&meshxb＝27.1,70.3,34.0,59.6,-0.1,17.9,ijk＝108,64,45/building
44.2)根据建筑表皮、楼屋面曲面信息等形成obst命令
45.参阅附图3，根据建筑表皮、楼屋面曲面信息等形成obst命令如下：
46.&obstid＝'gh-xy-0',xb＝45.1,45.5,40,40,-0.1,-0.1,rgb＝61,129,161/；
47.&obstid＝'gh-xz-0',xb＝36.3,36.3,40,40.4,11.1,11.5,rgb＝161,157,61/；
48.&obstid＝'gh-xz-0'，xb＝27.1，27.1，39.2，39.2，12.7，13.5，rgb＝161,61,61/；
49.3)设定最不利火灾场景，火源功率取qs＝8mw，火灾增长类型选择“快速”，火床面积设为2.4m
×
2.4m，并给定火床位置。根据整体网格划分情况，形成局部精细网格，其命令如下：
50.&mesh xb＝55.9,58.3，45.2,47.6,-0.1,2.3,ijk＝12,12,12/fire；
51.火源位置的obst命令、热边界条件的surf命令、火灾时间的time命令如下：
52.&obst xb＝56.7，59.1,42.0,44.4,-0.1,-0.1,surf_id＝'fire'/；
53.&surfid＝'fire',color＝'red',hrrpua＝1388.9,tau_q＝-413/；
54.&time t_end＝1800.0/；
55.根据结构分析模型中提取的节点信息、单元信息，形成获取构件所在位置空气升温的devc命令，并在命令中设定带结构单元编号的id信息，形成的部分命令如下：
56.&devcid＝'device1a',quantity＝'temperature',xyz＝3.2,51.2,5.7/；
57.&devcid＝'device2a',quantity＝'temperature',xyz＝9.2,51.2,5.7/；
58.&devcid＝'device3a',quantity＝'temperature',xyz＝18.2,51.2,5.7/。
59.3、火灾场景模拟
60.参阅附图4，运行生成的fds文件，进行火灾场景分析，并完成求解，得到的空间温度分布情况。
61.4、计算各构件升温曲线
62.统一采用非膨胀型防火层，等效热传导系数为0.1w/(m
·
℃)，厚度为30mm。根据mgt文件中提取的单元信息，确定每根构件的截面编号；根据mgt文件中提取的截面信息，确定各编号截面的截面形状系数fi/v，其中，fi为有防火保护钢构件单位长度的受火表面积(m2)，v为单位长度钢构件的体积(m3)。从fds求解得到的_devc.csv文件中，提取各测点温度时程数据，并根据id号中的结构单元编号，将温度数据定位到每一根构件，按《建筑钢结构防火技术规范》gb51249-2017第6.2.2条计算火灾下有防火保护构件的温度情况，即下述a式计算各构件的升温曲线：
[0063][0064]
当防火保护层为非轻质防火保护层，即2ρicidifi》ρiciv时，α则由下述b式计算：
[0065][0066]
当防火保护层为轻质防火保护层，即2ρicidifi≤ρiciv时，对非膨胀型防火涂料、防火板等防火保护层，α则由下述c式计算：
[0067][0068]
式中，δts为钢构件在时间(t,t+δt)内的温升(℃)；α为综合热传递系数[w/(m2·
℃)]；ρs、cs分别为钢材的密度(kg/m3)和比热[j/(kg
·
℃)]；ts、tg为t时刻钢构件的内部温度和热烟气的平均温度(℃)；δt为时间步长(s)；ci为防火保护材料的比热容[j/(kg
·
℃)]；ρi为防火保护材料的密度(kg/m3)；λi为防火保护材料的等效热传导系数，取0.1w/(m
·
℃)；di为防火保护层的厚度，取0.03m。
[0069]
参阅附图5，设定受火时间为1.5小时，空气温度数据tg取包络值，求解得到每根构件在达到整体结构耐火极限时的温度情况。
[0070]
5、整体结构耐火验算模型的建立及求解
[0071]
根据结构mgt文本文件中的结构布置及荷载信息，采用ansys ls-dyna软件建立整体结构耐火验算模型，并将自动计算得到的各构件温度，形成单元温度荷载命令，添加到模
型中，进行整体结构在火灾荷载工况下的结构响应分析。
[0072]
6、整体结构耐火验算及防火保护措施调整
[0073]
参阅附图6，整体结构耐火验算模型中的温度分布情况与附图5中各构件升温求解情况一致。
[0074]
参阅附图7，在该工况下构件的应变最大值为0.0006，所有构件均保持弹性状态而未形成塑性铰，结构保持稳定，没有达到耐火承载力极限状态。由分析结果可见，采用的防火层厚度能够确保该火灾场景下整体结构的安全性。
[0075]
以上实施例只是对本发明做进一步说明，并非用以限制本发明专利，凡为本发明等效实施，均应包含于本发明专利的权利要求范围之内。