CN106503390A - 一种板翅式换热器的蠕变疲劳强度设计方法 - Google Patents

一种板翅式换热器的蠕变疲劳强度设计方法 Download PDF

Info

Publication number
CN106503390A
CN106503390A CN201610984992.5A CN201610984992A CN106503390A CN 106503390 A CN106503390 A CN 106503390A CN 201610984992 A CN201610984992 A CN 201610984992A CN 106503390 A CN106503390 A CN 106503390A
Authority
CN
China
Prior art keywords
heat exchanger
fin heat
plate
equivalent
fin
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.)
Granted
Application number
CN201610984992.5A
Other languages
English (en)
Other versions
CN106503390B (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 Gelan Environmental Technology Co ltd
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201610984992.5A priority Critical patent/CN106503390B/zh
Publication of CN106503390A publication Critical patent/CN106503390A/zh
Application granted granted Critical
Publication of CN106503390B publication Critical patent/CN106503390B/zh
Priority to PCT/CN2017/101922 priority patent/WO2018086411A1/zh
Priority to US15/737,008 priority patent/US10289772B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21D—WORKING OR PROCESSING OF SHEET METAL OR METAL TUBES, RODS OR PROFILES WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
    • B21D53/00—Making other particular articles
    • B21D53/02—Making other particular articles heat exchangers or parts thereof, e.g. radiators, condensers fins, headers
    • 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]
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23P—METAL-WORKING NOT OTHERWISE PROVIDED FOR; COMBINED OPERATIONS; UNIVERSAL MACHINE TOOLS
    • B23P15/00—Making specific metal objects by operations not covered by a single other subclass or a group in this subclass
    • B23P15/26—Making specific metal objects by operations not covered by a single other subclass or a group in this subclass heat exchangers or the like
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F28—HEAT EXCHANGE IN GENERAL
    • F28F—DETAILS OF HEAT-EXCHANGE AND HEAT-TRANSFER APPARATUS, OF GENERAL APPLICATION
    • F28F3/00—Plate-like or laminated elements; Assemblies of plate-like or laminated elements
    • F28F3/02—Elements or assemblies thereof with means for increasing heat-transfer area, e.g. with fins, with recesses, with corrugations
    • F28F3/025—Elements or assemblies thereof with means for increasing heat-transfer area, e.g. with fins, with recesses, with corrugations the means being corrugated, plate-like elements
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • 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
    • 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/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • Evolutionary Computation (AREA)
  • Thermal Sciences (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
  • Pressure Welding/Diffusion-Bonding (AREA)

Abstract

一种板翅式换热器的蠕变疲劳强度设计方法，属于换热器技术领域。根据服役要求对板翅式换热器进行初步设计，再对板翅式换热器进行一次应力评定，计算满足许用应力要求的板翅式换热器芯体的等效力学参数和等效热物性参数，以此为基础对板翅式换热器进行热疲劳分析，再计算板翅式换热器的疲劳寿命和蠕变寿命，实现高温服役板翅式换热器的综合设计。本板翅式换热器的蠕变疲劳强度设计方法解决了板翅式换热器由于复杂的周期性结构而无法直接运用有限元软件进行高温强度设计的问题，为板翅式换热器的高温强度设计提供了理论基础，从而对在高温、交变载荷下服役的板翅式换热器进行有效的寿命预测，为板翅式换热器的高温设计提供了有效的方法。

Description

一种板翅式换热器的蠕变疲劳强度设计方法
技术领域
一种板翅式换热器的蠕变疲劳强度设计方法，属于换热器技术领域。
背景技术
随着科技的发展，能源消耗也越来越多，节能减排，提高能源的利用率就成为人们关注的焦点。换热设备，作为高温系统的核心一员，不仅仅要求具有高效的换热性能，而且要求有紧凑的结构。但是，现有的换热设备，多为管壳式换热器，占用体积大，换热效率低，难以满足在航空航天、高温气冷堆、燃气轮机等领域的要求。
板翅式换热器具有结构紧凑，换热效率高的特点，研究板翅式换热器很有前景。然而，板翅式换热器的服役环境越来越恶劣，高温高压的环境使得板翅式换热器的设计要求越来越高。在高温、交变载荷下服役，破坏模式与时间相关，现有的压力容器设计规范局限于管壳式换热器，以弹塑性破坏理论为基础，未考虑蠕变与疲劳破坏的特点，没有涉及钎焊过程、服役环境等因素的影响，不能直接用于板翅式换热器的设计。
发明内容
本发明要解决的技术问题是：克服现有技术的不足，提供一种能够结合钎焊过程、服役环境、失效位置等进行板翅式换热器高温强度设计的板翅式换热器的蠕变疲劳强度设计方法。
本发明解决其技术问题所采用的技术方案是：该板翅式换热器的蠕变疲劳强度设计方法，其特征在于：包括如下步骤：
步骤1，根据设计温度、设计压力要求对板翅式换热器结构进行初步设计，并明确板翅式换热器的工作温度、运行循环次数以及服役寿命；
步骤2，通过有限元软件进行板翅结构一次应力分析，确定应力集中部位，并确定许用应力；
步骤3，判断应力集中部位的应力水平是否满足如下条件：
；；
其中，为一次薄膜应力，为局部薄膜应力，为一次弯曲应力，为与时间有关的许用应力，的取值范围为1.05~1.16；
若满足条件，则执行步骤4；若一次应力评定不满足条件，则改变板翅式换热器芯体的结构、板材厚度，返回步骤2；
步骤4，在服役环境下，对板翅结构进行蠕变断裂实验和疲劳实验，对老化母材进行蠕变断裂实验和疲劳实验，计算应力放大系数和应变放大系数，并根据实验结果对母材的疲劳设计曲线和蠕变断裂设计曲线进行修正；
，，
其中，、分别为相同蠕变断裂时间下母材和板翅结构的蠕变断裂强度，
、分别为相同疲劳寿命下母材和板翅结构的宏观应变范围；
步骤5，获取板翅结构的等效力学参数和等效热物性参数，从而对板翅式换热器进行热疲劳有限元分析，得出板翅式换热器芯体沿高度方向的宏观应力的时间历程，计算钎角处总应变，
，
其中，为由一次应力分析所得应力范围得出的应变范围；
为热疲劳分析所得宏观应力的最大值与最小值的差值与板翅式换热器芯体高度方向弹性模量的比值；
步骤6，计算板翅式换热器的疲劳损伤和蠕变损伤，
，
其中，为疲劳周期数，
为修正后的疲劳设计曲线上，应变范围为ε时对应的疲劳寿命；
，
其中，为疲劳周期数，
为应变保持时间，
为t时刻的宏观应力，
为修正后的蠕变断裂设计曲线上，应力为σ时对应的蠕变断裂寿命；
步骤7，如果小于1，则执行步骤8；如果大于等于1，则执行步骤1；
步骤8，板翅式换热器设计完成。
优选的，步骤2和步骤3中所述的许用应力包括翅片区许用应力以及封条区许用应力。
优选的，步骤5中所述的获取板翅式换热器芯体的等效力学参数和等效热物性参数包括如下步骤：
a，将板翅式换热器芯体划分为多个形状相同的板翅胞元；
b，将板翅胞元等效成均匀固态板；
c，获取一个板翅胞元的等效力学参数和等效热物性参数，从而获得整个板翅式换热器芯体的等效力学参数和等效热物性参数。
优选的，所述的等效力学参数包括各向异性等效弹性模量、等效剪切模量、泊松比；所述的等效热物性参数包括等效导热系数、等效热膨胀系数、等效密度和等效比热。
优选的，所述的各向异性等效弹性模量的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道轴线平行的方向为x轴，垂直于流道轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，、、分别为x轴、y轴、z轴方向的等效弹性模量，
为母材的弹性模量，
d为板翅式换热器芯体的流道上侧的宽度，
为流道侧面与竖直面的夹角，
为流道的侧面的长度，
t为板翅式换热器芯体平板的厚度，
δ为板翅式换热器芯体的翅片的厚度。
优选的，泊松比的计算公式如下：
，
，
，
其中，为在y轴方向载荷作用下x轴方向与y轴方向应变的比值，
为在z轴方向载荷作用下x轴方向与z轴方向应变的比值，
为在z轴方向载荷作用下y轴方向与z轴方向应变的比值，
为母材的泊松比。
优选的，所述的等效剪切模量的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道轴线平行的方向为x轴，垂直于流道轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，、分别为母材的泊松比和弹性模量，
d为板翅式换热器芯体的流道上侧的宽度，
为流道侧面与竖直面的夹角，
为流道的侧面的长度，
t为板翅式换热器芯体平板的厚度，
δ为板翅式换热器芯体的翅片的厚度。
优选的，所述的等效导热系数的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道轴线平行的方向为x轴，垂直于流道轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，、、分别为x轴方向、y轴方向、z轴方向的等效导热系数，
、分别为母材和空气的导热系数，
d为板翅式换热器芯体的流道上侧的宽度，
为流道侧面与竖直面的夹角，
为流道的侧面的长度，
t为板翅式换热器芯体平板的厚度，
δ为板翅式换热器芯体的翅片的厚度。
优选的，所述的等效热膨胀系数的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道轴线平行的方向为x轴，垂直于流道轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，为等效母材热膨胀系数，
d为板翅式换热器芯体的流道上侧的宽度，
为流道侧面与竖直面的夹角，
为流道的侧面的长度，
t为板翅式换热器芯体平板的厚度，
δ为板翅式换热器芯体的翅片的厚度。
优选的，所述的等效密度和等效比热的计算公式如下：
，
，
，
，
其中，、分别为母材体积分数和空气体积分数，
、为母材密度和空气密度，
、分别为母材和空气的比热，
、分别为等效比热和等效密度，
d为板翅式换热器芯体的流道上侧的宽度，
为流道侧面与竖直面的夹角，
为流道的侧面的长度，
t为板翅式换热器芯体平板的厚度，
δ为板翅式换热器芯体的翅片的厚度。
与现有技术相比，本发明所具有的有益效果是：
1、本板翅式换热器的蠕变疲劳强度设计方法综合考虑了板翅式换热器的钎焊过程、服役环境、失效位置等影响，并对板翅式换热器芯体进行了等效均匀化，计算了板翅式换热器芯体等效力学参数和等效热物性参数，解决了板翅式换热器由于复杂的周期性结构而无法直接有效运用有限元软件进行高温强度设计的问题，为板翅式换热器的高温强度设计提供了理论基础，从而能够对在高温、交变载荷下服役的板翅式换热器进行有效的寿命预测，为高温服役的板翅式换热器的设计提供了有效的方法。
2、等效均匀化，将板翅结构分为相同的板翅胞元，并将板翅胞元等效成均匀固态板，进而计算板翅结构的等效力学参数计算和等效热物性参数，从而解决了由于板翅式换热器周期性复杂结构而导致的难以对板翅式换热器进行有限元模拟的问题，进而方便了后续通过有限元分析软件对板翅式换热器进行热疲劳分析。
3、板翅式换热器芯体的等效力学参数和等效热物性参数均以解析式的方式给出，解决了之前只能通过复杂的计算机模拟或实验获取参数的问题，方便了等效参数的计算，大大提高了板翅式换热器设计的效率。
附图说明
图1为板翅式换热器芯体的主视示意图。
图2为板翅胞元的主视示意图。
图3为蠕变疲劳寿命评价折线图。
图中：1、平板 2、翅片 3、流道。
具体实施方式
图1~3是本发明的最佳实施例，下面结合附图1~3对本发明做进一步说明。
该板翅式换热器的蠕变疲劳强度设计方法包括如下步骤：
步骤1，根据设计温度、设计压力要求对板翅式换热器结构进行初步设计，并明确板翅式换热器芯体的工作温度、运行循环次数以及服役寿命；
运行循环次数为设计使用年限与每年停机次数的乘积；服役寿命即设计寿命。
如图1~2所示：板翅式换热器芯体包括平板1和翅片2，每相邻的两块平板1之间设有翅片2，由多块平板1和翅片2间隔叠加钎焊而成，，从而在每相邻的两块平板1之间形成多个流道3，流道3的截面为等腰梯形。
步骤2，通过有限元软件进行板翅结构一次应力分析，确定应力集中部位，并确定许用应力；
在分析过程中不考虑钎焊焊缝对结构应力的影响，并且假设钎焊粘结率为100%。然后考虑材料在钎焊过程中的热力老化、服役环境（如高温气冷堆中氦气环境对材料强度的影响）、封条结构，进行实验分析来确定许用应力。
许用应力为与时间有关的许用应力，许用应力包括翅片区许用应力以及封条区许用应力。在获取翅片区许用应力时，针对钎焊高温和服役环境影响，对老化母材在服役环境中（如氦气环境）进行单轴拉伸和蠕变断裂实验；在获取封条区许用应力时，需要进行钎料的拉伸强度实验。最后根据ASME设计准则及修正结果，分别确定翅片区许用应力以及封条区许用应力。
翅片区许用应力取以下四个因子的最小值：
①设计温度下屈服应力*（1/1.1）*p，
②引起蠕变断裂的最小应力的67%*q，
③第三阶段蠕变开始的最小应力的80%*q，
④达到总应变（弹性、塑性、蠕变）1%的最小应力，
在封条结构应变进行限定，许用应变取钎料断裂延伸率的三分之一，封条区许用应力取以下三个因子的最小值：
①设计温度下母材屈服应力*（1/1.1）*p，
②许用应变*母材弹性模量，
③引起蠕变断裂的最小应力的67%*q，
其中，p为老化母材与未老化母材屈服应力的比值；
q为老化母材蠕变断裂强度和未老化母材的蠕变断裂强度的比值。
步骤3，判断应力集中部位的应力水平是否满足如下条件：
；；
其中，为一次薄膜应力，为局部薄膜应力，为一次弯曲应力，为与时间有关的许用应力，的取值范围为1.05~1.16；
若满足条件，则执行步骤4；若一次应力评定不满足条件，则改变板翅式换热器芯体的结构、板材厚度，返回步骤2。
上述判断条件成立的条件是：当翅片区应力集中部位应力达到翅片区的许用应力时，板翅结构失效；当封条结构应力集中部位达到封条区许用应力时，封条结构失效。板翅结构和封条结构任一个失效均认为板翅式换热器芯体失效，需重新设计降低结构应力水平。
在进行应力评价时，分别利用翅片区许用应力以及封条区的许用应力的进行评价，当翅片区应力满足许用应力以及封条区应力满足许用应力时，即作出一次应力评价合格的判断，当翅片区应力不满足许用应力或封条区的应力不满足许用应力时，则作出不满足一次应力评价的判断。
步骤4，在服役环境下，对板翅结构进行蠕变断裂实验和疲劳实验，对老化母材也进行蠕变断裂实验和疲劳实验，计算应力放大系数和应变放大系数，并根据实验结果对母材的疲劳设计曲线和蠕变断裂设计曲线进行修正；
，，
其中，、分别为相同蠕变断裂时间下母材和板翅结构的蠕变断裂强度，
、分别为相同疲劳寿命下母材和板翅结构的宏观应变范围；
步骤5，获取板翅结构的等效力学参数和等效热物性参数，从而对板翅换热器进行热疲劳有限元分析。
板翅式换热器芯体的平板1和翅片2通过钎焊而成，结构紧凑，不同于常规均质材料，板翅式换热器芯体具有周期性多孔结构的特征，板翅式换热器芯体结构的周期性和复杂性，很难直接进行有限元分析，需利用等效均匀化的方法进行有限元分析。为了能过对板翅式换热器芯体进行有限元分析，引入了均匀化方法。
均匀化方法是指复合材料具有规则或近似规则的结构，这种相当规则的异质性材料可以假设为具有周期性的结构，应该强调的是，相比于复合体的尺寸规模，这些非均匀材质是很小的。鉴于此，这些类型的材料有时候被称为具有周期性微观结构的复合材料。但是即使运用现代高速计算机分析这些包含大量异质性材料的边界值问题也是相当困难的。克服这种困难就需要找到一种方法，利用一种等效的材料模型代替复合材料，这个过程被称为均匀化。均匀化的实质是将等效材料代替周期性结构的复合材料，获得等效材料的性能参数，这是进行均匀化的关键步骤。
获取板翅式换热器芯体的等效力学参数和等效热物性参数包括如下步骤：
步骤a，将板翅式换热器芯体划分为多个形状相同的板翅胞元；
在本实施例中，如图2所示的结构为一个板翅胞元，从而可以将板翅式换热器芯体看做多个板翅胞元组合而成。
步骤b，将板翅胞元等效成均匀固态板；
由于板翅胞元的结构不是均匀的，将板翅胞元看成是均质材料，即将板翅胞元等效成均匀固态板，以等效固态板来代替不均匀的板翅结构。
步骤c，获取一个板翅胞元的等效力学参数和等效热物性参数，从而获得整个板翅式换热器芯体的等效力学参数和等效热物性参数。这里的等效力学参数和等效热物性参数也可以通过有限元分析软件或实验的方法获得。
等效力学参数包括各向异性等效弹性模量、等效剪切模量、泊松比；等效热物性参数包括等效导热系数、等效热膨胀系数、等效密度和等效比热。以板翅式换热器前侧底部的中点为原点，以水平面内与流道3轴线平行的方向为x轴，垂直于流道3轴线的方向为y轴，竖直方向为z轴建立坐标系，从而对板翅式换热器芯体的等效力学参数和等效热物性参数进行计算。
板翅式换热器芯体的等效弹性模量的计算方法如下：
计算z轴方向的等效弹性模量，利用平板1上所受的力与翅片2垂直部分所受的力平衡，
，
计算x轴方向的等效弹性模量，利用等效应力和实际应变的概念，
，
计算y轴方向的等效弹性模量，利用等效应力和实际应变的概念，
，
其中，、、分别为x轴、y轴、z轴方向的等效弹性模量，
为母材的弹性模量，
d为板翅式换热器芯体的流道3上侧的宽度，
为流道3侧面与竖直面的夹角，
为流道3的侧面的长度，
t为板翅式换热器芯体平板1的厚度，
δ为板翅式换热器芯体的翅片2的厚度。
板翅式换热器芯体的泊松比的计算方法如下：
由于平板加强作用，计算，先计算，即先计算在x轴方向载荷作用下y轴方向与x轴方向应变的比值，然后根据弹性模量和泊松比的关系即可得出，
即：由，
得出，
同理可得：
，
，
其中，为在y轴方向载荷作用下x轴方向与y轴方向应变的比值，
为在z轴方向载荷作用下x轴方向与z轴方向应变的比值，
为在z轴方向载荷作用下y轴方向与z轴方向应变的比值，
为母材的泊松比。
板翅式换热器芯体的等效剪切模量的计算方法如下：
计算，对各向同性均质材料，有，
从而得出：，
计算和，、均为等效切应力与实际切应变的比值，
从而得出：
，
，
其中，、分别为母材的泊松比和弹性模量。
板翅式换热器芯体的等效导热系数的计算方法如下：
最小热阻力法则，也称为并联法则，热量在物体内传递时，热流会沿阻力最小的通道传递，或通道在流过定向热流量时呈最小热阻力状态，相应通道的总热阻即为最小热阻，也称等效热阻。等效导热系数法则，只要复合材料的单元体在与总体有相等的比等效热阻，不论单元体尺寸大小，在只考虑热传导时，这种单元体与总体的等效导热系数相等。
由以上可知，要求板翅式换热器芯体整体的导热系数，只需要求一个板翅胞元的等效导热系数即可。
，
，
，
其中，、、分别为x轴方向、y轴方向、z轴方向的等效导热系数，
、分别为母材和空气的导热系数。
板翅式换热器芯体的等效热膨胀系数的计算方法如下：
由于一个板翅胞元的上侧和下侧的平板1可自由膨胀，翅片2部分水平部分和倾斜部分由于在z轴方向的膨胀量不同，有相互作用，从而得出：
，
由于平板1和翅片2在y轴方向和x轴方向的热膨胀量相同，所以有
，
，
其中，为等效母材热膨胀系数。
板翅式换热器芯体的等效密度和等效比热的计算方法如下：
，
，
，
，
其中，、分别为母材体积分数和空气体积分数，
、为母材密度和空气密度，
、分别为母材和空气的比热，
、分别为等效比热和等效密度。
利用计算的等效力学参数和等效热物性参数，通过有限元分析软件，运用均匀化方法进行热疲劳的各向异性弹性分析，从热应力分析的结果，得出板翅结构叠加方向，即z轴方向的宏观应力的时间历程，是宏观应力最大值与最小值的差值，其与z轴方向弹性模量的比值，即为应变范围；由一次应力范围得出的应变范围，从而计算钎角处总应变，
。
步骤6，计算板翅式换热器的疲劳损伤和蠕变损伤，
，
其中， 为疲劳周期数，
为修正后的疲劳设计曲线上，应变范围为ε时对应的疲劳寿命；
，
其中，为疲劳周期数，
为应变保持时间，
为t时刻的宏观应力，
为修正后的蠕变断裂设计曲线上应力为σ时对应的蠕变断裂寿命。
疲劳损伤为老化母材在温度最大时的疲劳设计曲线中板翅式换热器芯体各点运行循环数与允许循环数的比值。计算蠕变损伤时，由于保持时间内，会发生应力松弛，最好可得出结构的应力松弛曲线，即的变化曲线。
步骤7，如果小于1，则执行步骤8；如果大于等于1，则执行步骤1；
根据ASME中蠕变-疲劳损伤评定准则，如图3所示，根据算出的总蠕变损伤和总疲劳损伤，以为横坐标，为纵坐标，为裂纹萌生的包络线，如果小于1，即位于包络线之下，则说明板翅式换热器在设计温度和压力下，整个设计寿命内不会失效，符合设计要求，执行步骤8。如果大于等于1，说明板翅式换热器不满足设计要求，此时需改进结构，更换材料，在允许条件下降低服役压力、温度等，重新执行步骤1，直到位于包络线下，完成板翅式换热器的高温强度设计。
步骤8，板翅式换热器设计完成。
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (10)

1.一种板翅式换热器的蠕变疲劳强度设计方法，其特征在于：包括如下步骤：
步骤1，根据设计温度、设计压力要求对板翅式换热器结构进行初步设计，并明确板翅式换热器的工作温度、运行循环次数以及服役寿命；
步骤2，通过有限元软件进行板翅结构一次应力分析，确定应力集中部位，并确定许用应力；
步骤3，判断应力集中部位的应力水平是否满足如下条件：
；；
其中，为一次薄膜应力，为局部薄膜应力，为一次弯曲应力，为与时间有关的许用应力，的取值范围为1.05~1.16；
若满足条件，则执行步骤4；若一次应力评定不满足条件，则改变板翅式换热器芯体的结构、板材厚度，返回步骤2；
步骤4，在服役环境下，对板翅结构进行蠕变断裂实验和疲劳实验，对老化母材进行蠕变断裂实验和疲劳实验，计算应力放大系数和应变放大系数，并根据实验结果对母材的疲劳设计曲线和蠕变断裂设计曲线进行修正；
，，
其中，、分别为相同蠕变断裂时间下母材和板翅结构的蠕变断裂强度，
、分别为相同疲劳寿命下母材和板翅结构的宏观应变范围；
步骤5，获取板翅结构的等效力学参数和等效热物性参数，从而对板翅式换热器进行热疲劳有限元分析，得出板翅式换热器芯体沿高度方向的宏观应力的时间历程，计算钎角处总应变，
，
其中，为由一次应力分析所得应力范围得出的应变范围；
为热疲劳分析所得宏观应力的最大值与最小值的差值与板翅式换热器芯体高度方向弹性模量的比值；
步骤6，计算板翅式换热器芯体的疲劳损伤和蠕变损伤，
，
其中， 为疲劳周期数，
为修正后的疲劳设计曲线上，应变范围为ε时对应的疲劳寿命；
，
其中，为疲劳周期数，
为应变保持时间，
为t时刻的宏观应力，
为修正后的蠕变断裂设计曲线上，应力为σ时对应的蠕变断裂寿命；
步骤7，如果小于1，则执行步骤8；如果大于等于1，则执行步骤1；
步骤8，板翅式换热器设计完成。
2.根据权利要求1所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：步骤2和步骤3中所述的许用应力包括翅片区许用应力以及封条区许用应力。
3.根据权利要求1所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：步骤5中所述的获取板翅式换热器芯体的等效力学参数和等效热物性参数包括如下步骤：
a，将板翅式换热器芯体划分为多个形状相同的板翅胞元；
b，将板翅胞元等效成均匀固态板；
c，获取一个板翅胞元的等效力学参数和等效热物性参数，从而获得整个板翅式换热器芯体的等效力学参数和等效热物性参数。
4.根据权利要求3所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：所述的等效力学参数包括各向异性等效弹性模量、等效剪切模量、泊松比；所述的等效热物性参数包括等效导热系数、等效热膨胀系数、等效密度和等效比热。
5.根据权利要求4所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：所述的各向异性等效弹性模量的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道（3）轴线平行的方向为x轴，垂直于流道（3）轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，、、分别为x轴、y轴、z轴方向的等效弹性模量，
为母材的弹性模量，
d为板翅式换热器芯体的流道（3）上侧的宽度，
为流道（3）侧面与竖直面的夹角，
为流道（3）的侧面的长度，
t为板翅式换热器芯体平板（1）的厚度，
δ为板翅式换热器芯体的翅片（2）的厚度。
6.根据权利要求5所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：泊松比的计算公式如下：
，
，
，
其中，为在y轴方向载荷作用下x轴方向与y轴方向应变的比值，
为在z轴方向载荷作用下x轴方向与z轴方向应变的比值，
为在z轴方向载荷作用下y轴方向与z轴方向应变的比值，
为母材的泊松比。
7.根据权利要求4所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：所述的等效剪切模量的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道（3）轴线平行的方向为x轴，垂直于流道（3）轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，、分别为母材的泊松比和弹性模量，
d为板翅式换热器芯体的流道（3）上侧的宽度，
为流道（3）侧面与竖直面的夹角，
为流道（3）的侧面的长度，
t为板翅式换热器芯体平板（1）的厚度，
δ为板翅式换热器芯体的翅片（2）的厚度。
8.根据权利要求4所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：所述的等效导热系数的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道（3）轴线平行的方向为x轴，垂直于流道（3）轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，、、分别为x轴方向、y轴方向、z轴方向的等效导热系数，
、分别为母材和空气的导热系数，
d为板翅式换热器芯体的流道（3）上侧的宽度，
为流道（3）侧面与竖直面的夹角，
为流道（3）的侧面的长度，
t为板翅式换热器芯体平板（1）的厚度，
δ为板翅式换热器芯体的翅片（2）的厚度。
9.根据权利要求4所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：所述的等效热膨胀系数的计算公式如下：
以板翅式换热器前侧底部的中点为原点，以水平面内与流道（3）轴线平行的方向为x轴，垂直于流道（3）轴线的方向为y轴，竖直方向为z轴建立坐标系，
，
，
，
其中，为等效母材热膨胀系数，
d为板翅式换热器芯体的流道（3）上侧的宽度，
为流道（3）侧面与竖直面的夹角，
为流道（3）的侧面的长度，
t为板翅式换热器芯体平板（1）的厚度，
δ为板翅式换热器芯体的翅片（2）的厚度。
10.根据权利要求4所述的板翅式换热器的蠕变疲劳强度设计方法，其特征在于：所述的等效密度和等效比热的计算公式如下：
，
，
，
，
其中，、分别为母材体积分数和空气体积分数，
、为母材密度和空气密度，
、分别为母材和空气的比热，
、分别为等效比热和等效密度，
d为板翅式换热器芯体的流道（3）上侧的宽度，
为流道（3）侧面与竖直面的夹角，
为流道（3）的侧面的长度，
t为板翅式换热器芯体平板（1）的厚度，
δ为板翅式换热器芯体的翅片（2）的厚度。
CN201610984992.5A 2016-11-09 2016-11-09 一种板翅式换热器的蠕变疲劳强度设计方法 Active CN106503390B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201610984992.5A CN106503390B (zh) 2016-11-09 2016-11-09 一种板翅式换热器的蠕变疲劳强度设计方法
PCT/CN2017/101922 WO2018086411A1 (zh) 2016-11-09 2017-09-15 一种板翅式换热器的蠕变疲劳强度设计方法
US15/737,008 US10289772B2 (en) 2016-11-09 2017-09-15 Method for creep-fatigue strength of plate-fin heat exchanger

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610984992.5A CN106503390B (zh) 2016-11-09 2016-11-09 一种板翅式换热器的蠕变疲劳强度设计方法

Publications (2)

Publication Number Publication Date
CN106503390A true CN106503390A (zh) 2017-03-15
CN106503390B CN106503390B (zh) 2017-08-25

Family

ID=58323957

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610984992.5A Active CN106503390B (zh) 2016-11-09 2016-11-09 一种板翅式换热器的蠕变疲劳强度设计方法

Country Status (3)

Country Link
US (1) US10289772B2 (zh)
CN (1) CN106503390B (zh)
WO (1) WO2018086411A1 (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330208A (zh) * 2017-07-06 2017-11-07 沈阳透平机械股份有限公司 一种增压机定子件分析方法
CN107842979A (zh) * 2017-11-30 2018-03-27 苏州惠林节能材料有限公司 一种节能板翅式热回收芯体及其制作方法
WO2018086411A1 (zh) * 2016-11-09 2018-05-17 中国石油大学(华东) 一种板翅式换热器的蠕变疲劳强度设计方法
CN109670241A (zh) * 2018-12-19 2019-04-23 中国石油大学(华东) 一种基于蠕变屈曲失效判定准则的有机玻璃承压结构的长期服役寿命预测方法
CN109766624A (zh) * 2019-01-04 2019-05-17 北京航空航天大学 一种胶接结构粘接层在高低温冷热循环条件下疲劳寿命的预测方法
CN109799104A (zh) * 2018-12-18 2019-05-24 上海市特种设备监督检验技术研究院 高压加热器弹塑性检测缩比样机
CN110581292A (zh) * 2019-08-13 2019-12-17 中国矿业大学 一种用于高温燃料电池堆的冷却器及热管理方法
CN113221238A (zh) * 2021-06-02 2021-08-06 沪东中华造船(集团)有限公司 一种船舶换热器可靠性评估方法
CN116822167A (zh) * 2023-06-08 2023-09-29 西安交通大学 换热器热力耦合性能多尺度分析方法、系统、介质及设备
CN109799104B (zh) * 2018-12-18 2024-05-31 上海市特种设备监督检验技术研究院 高压加热器弹塑性检测缩比样机

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110222363B (zh) * 2019-04-26 2023-04-07 华北水利水电大学 正交各向异性材料三维蠕变特性的表征及应用
CN110472270B (zh) * 2019-06-28 2022-09-27 合肥通用机械研究院有限公司 一种钠冷堆非能动余热排出系统的建模优化方法
CN111428184B (zh) * 2020-04-30 2023-06-23 庆安集团有限公司 一种板翅式换热器芯体尺寸计算方法
CN112100765B (zh) * 2020-08-28 2022-08-26 北京航空航天大学 一种涡轮盘冷挤压强化孔结构的高低周复合疲劳试验件及设计方法
CN112231953A (zh) * 2020-10-14 2021-01-15 西安热工研究院有限公司 一种提高核电站应急柴油机顶头座疲劳强度的方法
CN112347671A (zh) * 2020-10-28 2021-02-09 中国核动力研究设计院 一种反应堆压力容器密封结构三维数值模拟方法
CN116680838B (zh) * 2023-07-27 2024-04-26 东莞市鹏锦机械科技有限公司 一种板翅式换热器的传热计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5143563A (en) * 1989-10-04 1992-09-01 General Electric Company Creep, stress rupture and hold-time fatigue crack resistant alloys
JP2003114177A (ja) * 2001-10-02 2003-04-18 Mitsubishi Heavy Ind Ltd 熱サイクルクリープ疲労の試験片、その試験装置、及び、その試験方法
CN101625351A (zh) * 2008-07-08 2010-01-13 华东理工大学 一种蠕变数据转换为材料高温应力松弛数据的方法
CN103761365A (zh) * 2013-12-28 2014-04-30 合肥通用机械研究院 一种基于寿命的高温压力容器蠕变疲劳强度设计方法
CN105158084A (zh) * 2015-09-15 2015-12-16 华东理工大学 一种材料的蠕变-疲劳寿命预测方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6957695B2 (en) * 2003-05-13 2005-10-25 H2Gen Innovations, Inc. Heat exchanger housing and seals
JP2005227065A (ja) 2004-02-12 2005-08-25 Hitachi Ltd 構造物の推定寿命表示方法、寿命分布作成方法、画像表示装置
DE102006013503A1 (de) * 2006-03-23 2008-01-24 Esk Ceramics Gmbh & Co. Kg Plattenwärmetauscher, Verfahren zu dessen Herstellung und dessen Verwendung
US9764435B2 (en) * 2013-10-28 2017-09-19 Honeywell International Inc. Counter-flow heat exchange systems
CN106503390B (zh) 2016-11-09 2017-08-25 中国石油大学(华东) 一种板翅式换热器的蠕变疲劳强度设计方法
CN106355042B (zh) 2016-11-09 2017-12-01 中国石油大学(华东) 板翅式换热器的均匀化设计方法
CN106338431B (zh) 2016-11-09 2017-09-29 中国石油大学(华东) 用于确定板翅式换热器力学参数的方法及其应用设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5143563A (en) * 1989-10-04 1992-09-01 General Electric Company Creep, stress rupture and hold-time fatigue crack resistant alloys
JP2003114177A (ja) * 2001-10-02 2003-04-18 Mitsubishi Heavy Ind Ltd 熱サイクルクリープ疲労の試験片、その試験装置、及び、その試験方法
CN101625351A (zh) * 2008-07-08 2010-01-13 华东理工大学 一种蠕变数据转换为材料高温应力松弛数据的方法
CN103761365A (zh) * 2013-12-28 2014-04-30 合肥通用机械研究院 一种基于寿命的高温压力容器蠕变疲劳强度设计方法
CN105158084A (zh) * 2015-09-15 2015-12-16 华东理工大学 一种材料的蠕变-疲劳寿命预测方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018086411A1 (zh) * 2016-11-09 2018-05-17 中国石油大学(华东) 一种板翅式换热器的蠕变疲劳强度设计方法
US10289772B2 (en) 2016-11-09 2019-05-14 China University Of Petroleum (East China) Method for creep-fatigue strength of plate-fin heat exchanger
CN107330208A (zh) * 2017-07-06 2017-11-07 沈阳透平机械股份有限公司 一种增压机定子件分析方法
CN107330208B (zh) * 2017-07-06 2021-09-24 沈阳透平机械股份有限公司 一种增压机定子件分析方法
CN107842979A (zh) * 2017-11-30 2018-03-27 苏州惠林节能材料有限公司 一种节能板翅式热回收芯体及其制作方法
CN109799104A (zh) * 2018-12-18 2019-05-24 上海市特种设备监督检验技术研究院 高压加热器弹塑性检测缩比样机
CN109799104B (zh) * 2018-12-18 2024-05-31 上海市特种设备监督检验技术研究院 高压加热器弹塑性检测缩比样机
CN109670241A (zh) * 2018-12-19 2019-04-23 中国石油大学(华东) 一种基于蠕变屈曲失效判定准则的有机玻璃承压结构的长期服役寿命预测方法
CN109766624A (zh) * 2019-01-04 2019-05-17 北京航空航天大学 一种胶接结构粘接层在高低温冷热循环条件下疲劳寿命的预测方法
CN110581292A (zh) * 2019-08-13 2019-12-17 中国矿业大学 一种用于高温燃料电池堆的冷却器及热管理方法
CN113221238A (zh) * 2021-06-02 2021-08-06 沪东中华造船(集团)有限公司 一种船舶换热器可靠性评估方法
CN116822167A (zh) * 2023-06-08 2023-09-29 西安交通大学 换热器热力耦合性能多尺度分析方法、系统、介质及设备
CN116822167B (zh) * 2023-06-08 2024-04-16 西安交通大学 换热器热力耦合性能多尺度分析方法、系统、介质及设备

Also Published As

Publication number Publication date
CN106503390B (zh) 2017-08-25
US20180365356A1 (en) 2018-12-20
US10289772B2 (en) 2019-05-14
WO2018086411A1 (zh) 2018-05-17

Similar Documents

Publication Publication Date Title
CN106503390B (zh) 一种板翅式换热器的蠕变疲劳强度设计方法
CN106338431B (zh) 用于确定板翅式换热器力学参数的方法及其应用设计方法
CN106355042B (zh) 板翅式换热器的均匀化设计方法
Wang et al. Turbulence modeling of deep dynamic stall at relatively low Reynolds number
CN110991104B (zh) 焊接空间网架节点焊缝风致多轴高周疲劳损伤评定方法
CN104316388A (zh) 一种对各向异性材料结构件进行疲劳寿命测定的方法
CN107025315A (zh) 一种核电站蒸汽发生器u型传热管流致振动及微动磨损耦合分析计算方法
CN109543258A (zh) 一种基于模态叠加法的风电机组塔架应力状态计算方法
Li et al. Life prediction model of creep-rupture and creep-buckling of a pyramidal lattice truss panel structure by analytical and finite element study
Liu et al. Numerical simulation of flow past stationary and oscillating deformable circles with fluid-structure interaction
CN106354929B (zh) 基于刚度变化原理的承载结构载荷传递路径可视化方法
CN108009311A (zh) 一种针对蠕变试验的蠕变本构模型参数识别方法
Park et al. On the unsteady Reynolds-averaged Navier–Stokes capability of simulating turbulent boundary layers under unsteady adverse pressure gradients
Lotfi et al. 3D fluid-structure interaction (FSI) simulation of new type vortex generators in smooth wavy fin-and-elliptical tube heat exchanger
Zhang et al. Experimental and numerical investigations of aerodynamic loads and 3D flow over non‐rotating MEXICO blades
Gaheen et al. Fluid and structure analysis of wind turbine blade with winglet
Liu et al. Impact of saddle membrane structure by hail with combined particle sizes: Numerical simulation and experimental investigation
CN103177162A (zh) 一种基于交错迭代耦合技术的薄壁结构动力学热性预测方法
Tu et al. Creep of brazed plate-fin structures in high temperature compact heat exchangers
Ridluan et al. Benchmark simulation of turbulent flow through a staggered tube bundle to support CFD as a reactor design tool. Part I: SRANS CFD simulation
Yeh et al. Stress Analysis of Composite Wind Turbine Blade with Different Stacking Angle and Different Skin Thickness
Wang et al. Study on heat transfer and resistance characteristics of H-type finned tube
CN110032789A (zh) 气膜孔排布方式对镍基单晶合金性能影响的预测方法
Zhang et al. STRESS AND FATIGUE ANALYSIS OF AIR-TO-REFRIGERANT HEAT EXCHANGERS WITH NON-ROUND TUBE SHAPES
Rahman et al. Finite Element Structural Analysis of Commonly Used Horizontal Axis Wind Turbine Airfoils of Various Geometries

Legal Events

Date Code Title Description
C06 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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220921

Address after: 214500 1 Jianxin Road, Xinqiao Park, Jingjiang Economic Development Zone, Taizhou, Jiangsu

Patentee after: JIANGSU GELAN ENVIRONMENTAL TECHNOLOGY CO.,LTD.

Address before: No. 66, Changjiang West Road, Qingdao Economic Development Zone, Qingdao City, Shandong Province, 266580

Patentee before: CHINA University OF PETROLEUM (EAST CHINA)