CN115600461B - 一种计算四面受火下覆板复合保护方钢柱温度的方法 - Google Patents

一种计算四面受火下覆板复合保护方钢柱温度的方法 Download PDF

Info

Publication number
CN115600461B
CN115600461B CN202211259604.9A CN202211259604A CN115600461B CN 115600461 B CN115600461 B CN 115600461B CN 202211259604 A CN202211259604 A CN 202211259604A CN 115600461 B CN115600461 B CN 115600461B
Authority
CN
China
Prior art keywords
section
model
steel column
square steel
moment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211259604.9A
Other languages
English (en)
Other versions
CN115600461A (zh
Inventor
李�瑞
陈伟
尹亮
倪雅
叶继红
康壮苏
姜健
刘晴
叶薇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JIANGSU CUMT DAZHENG SURFACE ENGINEERING TECHNOLOGY CO LTD
Xuzhou Dehuai Technology Co ltd
China University of Mining and Technology CUMT
Tianjin Fire Research Institute of MEM
Original Assignee
JIANGSU CUMT DAZHENG SURFACE ENGINEERING TECHNOLOGY CO LTD
Xuzhou Dehuai Technology Co ltd
China University of Mining and Technology CUMT
Tianjin Fire Research Institute of MEM
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by JIANGSU CUMT DAZHENG SURFACE ENGINEERING TECHNOLOGY CO LTD, Xuzhou Dehuai Technology Co ltd, China University of Mining and Technology CUMT, Tianjin Fire Research Institute of MEM filed Critical JIANGSU CUMT DAZHENG SURFACE ENGINEERING TECHNOLOGY CO LTD
Priority to CN202211259604.9A priority Critical patent/CN115600461B/zh
Publication of CN115600461A publication Critical patent/CN115600461A/zh
Application granted granted Critical
Publication of CN115600461B publication Critical patent/CN115600461B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Building Environments (AREA)

Abstract

本发明公开了一种计算四面受火下覆板复合保护方钢柱温度的方法，本方法将覆板复合保护方钢柱断面内的各个保护层断面简化为阶梯状断面，对各阶梯状断面进行有限差分网格划分，建立有限差分显式或隐式分析模型；采用有限差分分析显式或隐式模型计算阶梯状断面各厚度处的温度；导出各时刻的各钢柱矩形断面与对应第n层保护层矩形断面交界处的最大温度，即得到四面受火下覆板复合保护方钢柱的方钢柱与复合保护交界处各时刻温度，本方法能够解决目前计算受火下方钢管柱温度需要繁琐建立虚拟仿真模型问题。

Description

一种计算四面受火下覆板复合保护方钢柱温度的方法
技术领域
本发明涉及受火下钢构件温度计算方法领域，具体涉及一种计算四面受火下覆板复合保护方钢柱温度的方法。
背景技术
钢柱是钢结构中重要构件，钢柱的破坏容易引起结构连续性倒塌，覆板复合保护是钢柱的常见防火保护。覆板复合保护为防火板(如石膏板、波特板等)与防火卷材(如岩棉、硅酸铝棉等)相互组合而成的防火保护，具有装饰性能好、无需养护、耐久性好、及施工过程环保无异味等优点。
当前计算覆板复合保护钢柱温度的方式一般为通过仿真模拟软件建立模型进行传热计算，此类方法缺点为需要专业人员进行建模且需要考虑复杂的热边界条件，同时计算耗时较长。
发明内容
针对上述存在的技术不足，本发明的目的是提供一种计算四面受火下覆板复合保护方钢柱温度的方法，其能够解决目前计算受火下方钢管柱温度需要繁琐建立虚拟仿真模型问题。
为解决上述技术问题，本发明采用如下技术方案：
一种计算四面受火下覆板复合保护方钢柱温度的方法，
所述覆板复合保护方钢柱断面由内向外包括方钢柱断面和复合保护断面，所述复合保护断面包括若干层保护层断面，该方法包括以下步骤：
步骤1：设覆板复合保护方钢柱断面的四个外角分别与对应的内角连线，将覆板复合保护方钢柱断面划分为四个梯形断面，每个梯形断面由外到内分别划分为第a层保护层梯形断面、第b层保护层梯形断面直至第n层保护层梯形断面以及方钢柱梯形断面；
步骤2：对所述梯形断面按断面面积相等的原则，沿厚度方向分别居中简化为第a层保护层矩形断面、第b层保护层矩形断面直至第n层保护层矩形断面以及方钢柱矩形断面，经过简化的各梯形断面均简化为阶梯状断面；
步骤3：对于阶梯状断面，其热量增量来源于对应的受火边界且除受火边界外各边界均简化为热绝缘边界，同时简化其热传导仅沿厚度方向传导；
步骤4：对各阶梯状断面进行有限差分网格划分，建立有限差分显式或隐式分析模型；
步骤5：采用有限差分分析显式或隐式模型计算阶梯状断面各厚度处的温度；
步骤6：导出各时刻的各方钢柱矩形断面与对应第n层保护层矩形断面交界处的最大温度，即得到四面受火下覆板复合保护方钢柱的方钢柱与复合保护交界处各时刻温度。
优选地，步骤4中所述的对阶梯状断面进行有限差分网格划分，建立有限差分显式或隐式分析模型；其具体包括以下步骤：
步骤4.1：阶梯状断面沿厚度方向均匀划分网格，形成若干单元，受火侧及柱腔侧边界单元的宽度为Δx/2，其余单元的宽度为Δx，每个单元的温度与相应节点温度一致，整个有限差分网格覆盖整个断面；
步骤4.2：当步骤4采用显式分析模型时，构建j时刻模型一维显式有限差分方程；
模型第a层保护层矩形断面受火侧边界单元节点：
Figure GDA0004143900150000031
模型第x层保护层矩形断面内部单元节点：
Figure GDA0004143900150000032
模型第x层保护层矩形断面与第x+1层保护层矩形断面边界单元节点：
Figure GDA0004143900150000033
模型第n层保护层矩形断面与方钢柱矩形断面边界单元节点：
Figure GDA0004143900150000034
模型方钢柱矩形断面内部单元节点：
Figure GDA0004143900150000035
模型方钢柱矩形断面的柱腔侧边界单元节点：
Figure GDA0004143900150000036
步骤4.3：当步骤4采用隐式分析模型时，构建j时刻模型一维隐式有限差分方程；
模型第a层保护层矩形断面受火侧边界单元节点：
Figure GDA0004143900150000041
模型第x层保护层矩形断面内部单元节点：
Figure GDA0004143900150000042
模型第x层保护层矩形断面与第x+1层保护层矩形断面边界单元节点：
Figure GDA0004143900150000043
模型第n层保护层矩形断面与方钢柱矩形断面边界单元节点：
Figure GDA0004143900150000044
模型方钢柱矩形断面内部单元节点：
Figure GDA0004143900150000045
模型方钢柱矩形断面的柱腔侧边界单元节点：
Figure GDA0004143900150000046
其中，Ti,j，Tg,j分别表示i节点及模型受火侧在j时刻的开尔文温度；kxi,j，cxi,j，ρxi,j分别表示第x层保护层矩形断面在i节点j时刻导热系数，比热容及密度；cxi,(x+1)i,j，ρxi,(x+1)i,j分别表示第x层保护层矩形断面与第x+1层保护层矩形断面的边界单元节点i在j时刻比热容及密度；cni,si,j，ρni,si,j分别表示第n层保护层矩形断面与方钢柱矩形断面的边界单元节点i在j时刻比热容及密度；ksi,j，csi,j，ρsi,j分别表示方钢柱矩形断面在i节点j时刻导热系数，比热容及密度；Δt为时间步长；εr为固体表面辐射率；σr为斯蒂芬-波尔茨曼常数；h表示j时刻模型受火侧表面对流换热系数；Dx，Ds分别表示第x层保护层矩形断面及方钢柱矩形断面边长。
优选地，步骤5中所述的采用有限差分分析隐式或显式模型计算阶梯状断面各厚度处的温度；其具体实现包括以下步骤：
步骤5.1：输入计算所需数据包括：受火侧各时刻温度、阶梯状断面各厚度位置的初始温度、各保护层的导热系数以及比热容和密度、钢柱的导热系数以及比热容和密度、总时长、时间步长、固体表面辐射率、斯蒂芬-波尔茨曼常数、受火侧表面对流换热系数、各保护层矩形断面及方钢柱矩形断面边长；
步骤5.2：当步骤4采用显式分析模型时，采用有限差分分析显式模型计算阶梯状断面各厚度处的温度：
步骤5.2.1：在第1时刻依次计算沿厚度方向从受火侧至柱腔侧各单元节点温度，得出第1时刻各单元节点温度；
步骤5.2.2：由上一时刻计算结果，再依次计算沿厚度方向从受火侧至柱腔侧各单元节点温度，得出下一时刻各单元节点温度，直至计算到总时长；
步骤5.3：当步骤4采用隐式分析模型时，采用有限差分分析隐式模型计算阶梯状断面各厚度处的温度：
步骤5.3.1：在第1时刻联立有限差分隐式模型中所有单元节点方程，计算得出第1时刻所有单元节点温度；
步骤5.3.2：由上一时刻计算结果，再联立下一时刻有限差分隐式模型中所有单元节点方程，计算得出下一时刻所有单元节点温度，直至计算到总时长。
优选地，步骤5.3中有限差分分析隐式模型计算单元节点温度中存在的高阶未知量，采用高斯-赛德尔迭代法对有限差分隐式模型单元节点方程求解。
本发明的有益效果在于：
1、本发明将计算四面受火下覆板复合保护钢柱温度的多维传热问题简化为一个一维传热问题，减少了计算难度。
2、本发明采用有限差分分析覆板复合保护钢柱传热模型，不同于以往需要通过建立仿真模型计算覆板复合保护钢柱温度，这节省了建模时间且易于操作，同时本发明给出的覆板复合保护钢柱一维隐式有限差分分析模型可以大幅提升计算效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的四面受火下覆板复合保护方钢柱断面示意图；
图2为本发明实施例提供的覆板复合保护方钢柱断面划分示意图；
图3为本发明实施例提供的各保护层矩形断面示意图；
图4为本发明实施例提供的阶梯状断面划分为若干单元的示意图。
附图标记说明：
1、覆板复合保护方钢柱断面，2、热源，3、复合保护断面，4、方钢柱断面，5、第a层保护层断面，6、第b层保护层断面，7、第n层保护层断面，8、梯形断面，9、第a层保护层梯形断面，10、第b层保护层梯形断面，11、第n层保护层梯形断面，12、方钢柱梯形断面，13、第a层保护层矩形断面，14、第b层保护层矩形断面，15、第n层保护层矩形断面，16、方钢柱矩形断面，17、阶梯状断面。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见图1-图4，一种计算四面受火下覆板复合保护方钢柱温度的简化方法，用于计算覆板复合保护方钢柱断面1的方钢柱与复合保护交界处温度，所述覆板复合保护方钢柱断面1外侧四面均设有热源2，所述覆板复合保护方钢柱断面1由外向内依次包括复合保护断面3及方钢柱断面4，所述复合保护断面3由若干层保护层断面组成，所述外层的保护层断面为防火板断面，所述内层的保护层断面为柔性毡断面，所述复合保护断面由外向内依次包括第a层保护层断面5、第b层保护层断面6直至第n层保护层断面7，所述第a层保护层断面、第b层保护层断面直至第n层保护层断面及方钢柱断面均为回形；
对所述覆板复合保护方钢柱断面的简化包括以下步骤1、步骤2及步骤3，对简化后的覆板复合保护方钢柱断面的钢柱温度计算包括以下步骤4、步骤5及步骤6：
步骤1：覆板复合保护方钢柱断面的四个外角分别与对应内角连线，将所述覆板复合保护方钢柱断面划分为四个梯形断面，同时将所述第a层保护层断面、第b层保护层断面直至第n层保护层断面及方钢柱断面分别划分为四个第a层保护层梯形断面9、四个第b层保护层梯形断面10直至四个第n层保护层梯形断面11及四个方钢柱梯形断面12；
步骤2：对于任意一个所述梯形断面，梯形断面的第a层保护层梯形断面、第b层保护层梯形断面直至第n层保护层梯形断面及方钢柱梯形断面按断面面积相等的原则，沿厚度方向分别居中简化为第a层保护层矩形断面13、第b层保护层矩形断面14直至第n层保护层矩形断面15及方钢柱矩形断面16，经过简化的各梯形断面均简化为阶梯状断面17；
步骤3：对于任意一个所述阶梯状断面，其热量增量来源于对应的受火边界且除受火边界外各边界均简化为热绝缘边界，同时简化其热传导仅沿厚度方向传导；
步骤4：对各阶梯状断面17进行有限差分网格划分，建立有限差分显式或隐式分析模型；
步骤4中所述的对阶梯状断面进行有限差分网格划分，建立有限差分显式或隐式分析模型；其具体实现包括以下步骤：
步骤4.1：阶梯状断面沿厚度方向均匀划分网格，受火侧及柱腔侧边界单元的宽度为Δx/2，其余单元的宽度为Δx，每个单元的温度与相应节点温度一致，整个有限差分网格覆盖整个断面；
步骤4.2：当步骤4采用显式分析模型时，构建j时刻模型一维显式有限差分方程；
模型第a层保护层矩形断面受火侧边界单元节点：
Figure GDA0004143900150000091
模型第x层保护层矩形断面内部单元节点：
Figure GDA0004143900150000092
模型第x层保护层矩形断面与第x+1层保护层矩形断面边界单元节点：
Figure GDA0004143900150000093
模型第n层保护层矩形断面与方钢柱矩形断面边界单元节点：
Figure GDA0004143900150000094
模型方钢柱矩形断面内部单元节点：
Figure GDA0004143900150000101
模型方钢柱矩形断面的柱腔侧边界单元节点：
Figure GDA0004143900150000102
步骤4.3：当步骤4采用隐式分析模型时，构建j时刻模型一维隐式有限差分方程；
模型第a层保护层矩形断面受火侧边界单元节点：
Figure GDA0004143900150000103
模型第x层保护层矩形断面内部单元节点：
Figure GDA0004143900150000104
模型第x层保护层矩形断面与第x+1层保护层矩形断面边界单元节点：
Figure GDA0004143900150000105
模型第n层保护层矩形断面与方钢柱矩形断面边界单元节点：
Figure GDA0004143900150000106
模型方钢柱矩形断面内部单元节点：
Figure GDA0004143900150000107
Figure GDA0004143900150000111
模型方钢柱矩形断面的柱腔侧边界单元节点：
Figure GDA0004143900150000112
其中，Ti,j，Tg,j分别表示i节点及模型受火侧在j时刻的开尔文温度；kxi,j，cxi,j，ρxi,j分别表示第x层保护层矩形断面在i节点j时刻导热系数，比热容及密度；cxi,(x+1)i,j，ρxi,(x+1)i,j分别表示第x层保护层矩形断面与第x+1层保护层矩形断面的边界单元节点i在j时刻比热容及密度；cni,si,j，ρni,si,j分别表示第n层保护层矩形断面与方钢柱矩形断面的边界单元节点i在j时刻比热容及密度；ksi,j，csi,j，ρsi,j分别表示方钢柱矩形断面在i节点j时刻导热系数，比热容及密度；Δt为时间步长；εr为固体表面辐射率；σr为斯蒂芬-波尔茨曼常数；h表示j时刻模型受火侧表面对流换热系数；Dx，Ds分别表示第x层保护层矩形断面及方钢柱矩形断面边长。
步骤5：采用有限差分分析显式或隐式模型计算阶梯状断面各厚度处的温度；
步骤5中所述的采用有限差分分析隐式或显式模型计算阶梯状断面各厚度处的温度；其具体实现包括以下步骤：
步骤5.1：输入计算所需数据包括：受火侧各时刻温度，阶梯状断面各厚度位置的初始温度，各保护层的导热系数、比热容及密度，钢柱的导热系数、比热容及密度，总时长，时间步长，固体表面辐射率，斯蒂芬-波尔茨曼常数，受火侧表面对流换热系数，各保护层矩形断面及方钢柱矩形断面边长；
步骤5.2：当步骤4采用显式分析模型时，采用有限差分分析显式模型计算阶梯状断面各厚度处的温度：
步骤5.2.1：在第一时刻依次计算沿厚度方向从受火侧至柱腔侧各单元节点温度，得出第一时刻各单元节点温度；
步骤5.2.2：由上一时刻计算结果，再依次计算沿厚度方向从受火侧至柱腔侧各单元节点温度，得出下一时刻各单元节点温度，直至计算到总时长；
步骤5.3：当步骤4采用隐式分析模型时，采用有限差分分析隐式模型计算阶梯状断面各厚度处的温度：
步骤5.3.1：在第一时刻联立有限差分隐式模型中所有单元节点方程，计算得出第一时刻所有单元节点温度；
步骤5.3.2：由上一时刻计算结果，再联立下一时刻有限差分隐式模型中所有单元节点方程，计算得出下一时刻所有单元节点温度，直至计算到总时长。
步骤5.3中有限差分分析隐式模型计算单元节点温度中存在的高阶未知量，采用高斯-赛德尔迭代法对有限差分隐式模型单元节点方程求解。
步骤6：导出各时刻的各方钢柱矩形断面与对应第n层保护层矩形断面交界处的最大温度，即得到四面受火下覆板复合保护方钢柱的方钢柱与复合保护交界处各时刻温度。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (3)

1.一种计算四面受火下覆板复合保护方钢柱温度的方法，
所述覆板复合保护方钢柱断面由内向外包括方钢柱断面和复合保护断面，所述复合保护断面包括若干层保护层断面，其特征在于，包括以下步骤：
步骤1：设覆板复合保护方钢柱断面的四个外角分别与对应的内角连线，将覆板复合保护方钢柱断面划分为四个梯形断面，每个梯形断面由外到内分别划分为第a层保护层梯形断面、第b层保护层梯形断面直至第n层保护层梯形断面以及方钢柱梯形断面；
步骤2：对所述梯形断面按断面面积相等的原则，沿厚度方向分别居中简化为第a层保护层矩形断面、第b层保护层矩形断面直至第n层保护层矩形断面以及方钢柱矩形断面，经过简化的各梯形断面均简化为阶梯状断面；
步骤3：对于阶梯状断面，其热量增量来源于对应的受火边界且除受火边界外各边界均简化为热绝缘边界，同时简化其热传导仅沿厚度方向传导；
步骤4：对各阶梯状断面进行有限差分网格划分，建立有限差分显式或隐式分析模型；其具体包括以下步骤：
步骤4.1：阶梯状断面沿厚度方向均匀划分网格，形成若干单元，受火侧及柱腔侧边界单元的宽度为Δx/2，其余单元的宽度为Δx，每个单元的温度与相应节点温度一致，整个有限差分网格覆盖整个断面；
步骤4.2：当步骤4采用显式分析模型时，构建j时刻模型一维显式有限差分方程；
模型第a层保护层矩形断面受火侧边界单元节点：
Figure FDA0004143900140000021
模型第x层保护层矩形断面内部单元节点：
Figure FDA0004143900140000022
模型第x层保护层矩形断面与第x+1层保护层矩形断面边界单元节点：
Figure FDA0004143900140000023
模型第n层保护层矩形断面与方钢柱矩形断面边界单元节点：
Figure FDA0004143900140000024
模型方钢柱矩形断面内部单元节点：
Figure FDA0004143900140000025
模型方钢柱矩形断面的柱腔侧边界单元节点：
Figure FDA0004143900140000026
步骤4.3：当步骤4采用隐式分析模型时，构建j时刻模型一维隐式有限差分方程；
模型第a层保护层矩形断面受火侧边界单元节点：
Figure FDA0004143900140000027
模型第x层保护层矩形断面内部单元节点：
Figure FDA0004143900140000031
模型第x层保护层矩形断面与第x+1层保护层矩形断面边界单元节点：
Figure FDA0004143900140000032
模型第n层保护层矩形断面与方钢柱矩形断面边界单元节点：
Figure FDA0004143900140000033
模型方钢柱矩形断面内部单元节点：
Figure FDA0004143900140000034
模型方钢柱矩形断面的柱腔侧边界单元节点：
Figure FDA0004143900140000035
其中，Ti,j，Tg,j分别表示i节点及模型受火侧在j时刻的开尔文温度；kxi,j，cxi,j，ρxi,j分别表示第x层保护层矩形断面在i节点j时刻导热系数，比热容及密度；cxi,(x+1)i,j，ρxi,(x+1)i,j分别表示第x层保护层矩形断面与第x+1层保护层矩形断面的边界单元节点i在j时刻比热容及密度；cni,si,j，ρni,si,j分别表示第n层保护层矩形断面与方钢柱矩形断面的边界单元节点i在j时刻比热容及密度；ksi,j，csi,j，ρsi,j分别表示方钢柱矩形断面在i节点j时刻导热系数，比热容及密度；Δt为时间步长；εr为固体表面辐射率；σr为斯蒂芬-波尔茨曼常数；h表示j时刻模型受火侧表面对流换热系数；Dx，Ds分别表示第x层保护层矩形断面及方钢柱矩形断面边长；
步骤5：采用有限差分分析显式或隐式模型计算阶梯状断面各厚度处的温度；
步骤6：导出各时刻的各方钢柱矩形断面与对应第n层保护层矩形断面交界处的最大温度，即得到四面受火下覆板复合保护方钢柱的方钢柱与复合保护交界处各时刻温度。
2.根据权利要求1所述的一种计算四面受火下覆板复合保护方钢柱温度的方法，其特征在于：步骤5中所述的采用有限差分分析隐式或显式模型计算阶梯状断面各厚度处的温度；其具体实现包括以下步骤：
步骤5.1：输入计算所需数据包括：受火侧各时刻温度、阶梯状断面各厚度位置的初始温度、各保护层的导热系数以及比热容和密度、钢柱的导热系数以及比热容和密度、总时长、时间步长、固体表面辐射率、斯蒂芬-波尔茨曼常数、受火侧表面对流换热系数、各保护层矩形断面及方钢柱矩形断面边长；
步骤5.2：当步骤4采用显式分析模型时，采用有限差分分析显式模型计算阶梯状断面各厚度处的温度：
步骤5.2.1：在第1时刻依次计算沿厚度方向从受火侧至柱腔侧各单元节点温度，得出第1时刻各单元节点温度；
步骤5.2.2：由上一时刻计算结果，再依次计算沿厚度方向从受火侧至柱腔侧各单元节点温度，得出下一时刻各单元节点温度，直至计算到总时长；
步骤5.3：当步骤4采用隐式分析模型时，采用有限差分分析隐式模型计算阶梯状断面各厚度处的温度：
步骤5.3.1：在第1时刻联立有限差分隐式模型中所有单元节点方程，计算得出第1时刻所有单元节点温度；
步骤5.3.2：由上一时刻计算结果，再联立下一时刻有限差分隐式模型中所有单元节点方程，计算得出下一时刻所有单元节点温度，直至计算到总时长。
3.根据权利要求2所述的一种计算四面受火下覆板复合保护方钢柱温度的方法，其特征在于：步骤5.3中有限差分分析隐式模型计算单元节点温度中存在的高阶未知量，采用高斯-赛德尔迭代法对有限差分隐式模型单元节点方程求解。
CN202211259604.9A 2022-10-14 2022-10-14 一种计算四面受火下覆板复合保护方钢柱温度的方法 Active CN115600461B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211259604.9A CN115600461B (zh) 2022-10-14 2022-10-14 一种计算四面受火下覆板复合保护方钢柱温度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211259604.9A CN115600461B (zh) 2022-10-14 2022-10-14 一种计算四面受火下覆板复合保护方钢柱温度的方法

Publications (2)

Publication Number Publication Date
CN115600461A CN115600461A (zh) 2023-01-13
CN115600461B true CN115600461B (zh) 2023-04-21

Family

ID=84845983

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211259604.9A Active CN115600461B (zh) 2022-10-14 2022-10-14 一种计算四面受火下覆板复合保护方钢柱温度的方法

Country Status (1)

Country Link
CN (1) CN115600461B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016020731A (ja) * 2014-07-15 2016-02-04 大成建設株式会社 すべり面温度推定方法、免震装置、および建物
CN106485034A (zh) * 2016-11-24 2017-03-08 南京工业大学 一种用于考虑胶缝导热效应的胶合木构件受火温度场的分析方法及其应用
CN106540967A (zh) * 2016-12-23 2017-03-29 东北大学 一种网格指数分布的宽厚板温度监控方法
CN114595606A (zh) * 2022-03-10 2022-06-07 福州大学 一种基于有限差分adi算法的喷雾冷却温度场求解方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7286972B2 (en) * 2001-04-17 2007-10-23 Livermore Software Technology Corporation Implicit-explicit switching for finite element analysis
CN101799843B (zh) * 2010-03-12 2012-09-26 哈尔滨工业大学深圳研究生院 注水空心钢管混凝土构件的防火性能评估方法及应用
CN104298884B (zh) * 2014-10-17 2017-06-16 武汉科技大学 一种快速计算轧件断面温度的有限元和有限差分耦合方法
JP6601762B2 (ja) * 2015-09-11 2019-11-06 株式会社日本製鋼所 鋼の熱処理シミュレーション方法および鋼の熱処理シミュレーションプログラム
CN108051118B (zh) * 2017-12-11 2019-07-05 国网新疆电力有限公司检修公司 一种测量复合空心绝缘子的玻璃钢管的温度的方法
CN111735307A (zh) * 2020-07-23 2020-10-02 应急管理部天津消防研究所 一种构件防火保护设计及防火监督用电炉装置及使用方法
CN114329967B (zh) * 2021-12-29 2024-04-16 中国人民警察大学 钢框架中柱温度应力的计算方法
CN115078442A (zh) * 2022-07-18 2022-09-20 南昌大学 一种火灾后混凝土构件表面受火温度检测方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016020731A (ja) * 2014-07-15 2016-02-04 大成建設株式会社 すべり面温度推定方法、免震装置、および建物
CN106485034A (zh) * 2016-11-24 2017-03-08 南京工业大学 一种用于考虑胶缝导热效应的胶合木构件受火温度场的分析方法及其应用
CN106540967A (zh) * 2016-12-23 2017-03-29 东北大学 一种网格指数分布的宽厚板温度监控方法
CN114595606A (zh) * 2022-03-10 2022-06-07 福州大学 一种基于有限差分adi算法的喷雾冷却温度场求解方法

Also Published As

Publication number Publication date
CN115600461A (zh) 2023-01-13

Similar Documents

Publication Publication Date Title
CN103886163B (zh) 一种涡轮叶片热障涂层的有限元模型的网格划分方法
CN114547790B (zh) 用于评估复杂多层热防护结构隔热性能的计算方法
CN115600461B (zh) 一种计算四面受火下覆板复合保护方钢柱温度的方法
Orabi et al. Integrated nonlinear structural simulation of composite buildings in fire
CN112983560B (zh) 一种涡轮气冷叶片温度场快速评估方法
Ren et al. Heterogeneous convective thermal and airborne pollutant removals from a partial building enclosure with a conducting baffle: Parametric investigations and steady transition flow solutions
CN108319775A (zh) 一种热传导问题中复合材料的近场动力学建模方法
CN116522457A (zh) 一种火灾后隧道衬砌爆裂损伤范围预测方法
CN109598059A (zh) 一种基于代理模型的热防护系统优化设计方法及设计系统
CN103886164A (zh) 含有多条冷却通道的涡轮叶片热障涂层有限元建模方法
Iqbal et al. An enriched finite element method for efficient solutions of transient heat diffusion problems with multiple heat sources
CN114218645B (zh) 物性时空变异的箱梁蒸养热力耦合场理论模型和计算方法
CN115186606A (zh) 一种多股流逆流板翅式换热器的校核方法、计算机程序产品
Yüncü et al. Laminar natural convective heat transfer in an air-filled parallelogramic cavity
Winwood et al. Advanced fabric energy storage II: Computational fluid dynamics modelling
CN116644565B (zh) 一种钢框架结构热-力耦合数值模拟方法
Singhal et al. Convective and radiative thermal analysis of composite wall with nonlinear temperature-dependent properties
Saitoh et al. Time‐space method for multidimensional melting and freezing problems
CN111832208B (zh) 一种混凝土复合自保温砌块热阻的确定方法
CN115422682A (zh) 一种风电机组混凝土基础水化热冷水管温控方法及系统
KR102413821B1 (ko) 원자로 용기의 외벽 냉각 성능을 평가하는 방법
Zeng et al. Study on the application of the hot-wall modified method in fluid-thermal-structure multi fields coupling strategy
JP4696345B2 (ja) ディスタンスブロックの配置方法
Urbikain et al. Determination of wall decay times by use of a polynomial equation
CN113961996B (zh) 一种基于参数化模型的整体结构耐火验算方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant