CN110411863A - 一种基于蠕变延性的高温蠕变寿命预测方法 - Google Patents

一种基于蠕变延性的高温蠕变寿命预测方法 Download PDF

Info

Publication number
CN110411863A
CN110411863A CN201810388364.XA CN201810388364A CN110411863A CN 110411863 A CN110411863 A CN 110411863A CN 201810388364 A CN201810388364 A CN 201810388364A CN 110411863 A CN110411863 A CN 110411863A
Authority
CN
China
Prior art keywords
creep
life
temperature
ductility
stress
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
CN201810388364.XA
Other languages
English (en)
Other versions
CN110411863B (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201810388364.XA priority Critical patent/CN110411863B/zh
Publication of CN110411863A publication Critical patent/CN110411863A/zh
Application granted granted Critical
Publication of CN110411863B publication Critical patent/CN110411863B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/28—Investigating ductility, e.g. suitability of sheet metal for deep-drawing or spinning
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001—Type of application of the stress
    • G01N2203/0003—Steady
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/006—Crack, flaws, fracture or rupture
    • G01N2203/0067—Fracture or rupture
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0071—Creep
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/022—Environment of the test
    • G01N2203/0222—Temperature
    • G01N2203/0226—High temperature; Heating means

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供了一种基于蠕变延性的高温蠕变寿命预测模型，基于能量守恒和熵守恒建立的蠕变寿命与蠕变延性的相关性，并考虑蠕变延性随蠕变应变速率的变化，提出以下预测模型：本发明可实现不同蠕变机制下高温蠕变寿命的预测，预测可靠性高。

Description

一种基于蠕变延性的高温蠕变寿命预测方法
技术领域
本发明涉及高温构件蠕变寿命技术领域，特别是涉及一种基于蠕变延性的高温蠕变寿命预测方法，考虑蠕变延性随蠕变应变速率的变化，实现金属材料长时蠕变寿命的预测。
背景技术
应用在火力发电站等大型设备中的高温构件，在高压高温等环境中，长时间承受各种应力的作用下，易于发生蠕变断裂，造成构件的失效，从而危及大型设备的正常安全使用，因此，精准预测蠕变寿命对高温构件的使用和设计具有非常重要的意义。寿命预测主要是以传统持久强度试验外推为主的方法，通过提高服役温度或加载应力的方法达到缩短实验时间的目的。蠕变寿命评估研究方法经过长期的发展已逐渐成熟大体经历了不同的阶段。最早是等温线外推法，但大量的实验结果表明在恒定的温度，双对数坐标下，时间和应力之间呈直线关系。在部分耐热钢的数据中是符合这一关系，但实际中有相当部分的耐热钢的持久性能试验点并不符合线性关系，低应力、长时间数据经常偏离原有的直线关系，相关研究也表明等温线外推法在评估新型耐热钢或不锈钢的持久性能时有高估的情况。为综合考虑温度和应力与蠕变持久断裂时间三者之间的关系，工程上将断裂时间和试验温度表示为时间-温度互补的参数，并和应力相关联。此方法简称TTP参数法，国内研究运用最多的是L-M参数法，虽然被广泛应用，但对其模型基础及预测精度一直存在争议，由于该方法比较难以描述蠕变机制变化引起的寿命变化。蠕变机制的变化与蠕变激活能有关，即蠕变激活能随持久时间而变化，提出分区处理的拟合方法，可以提高长时低应力下预测机制变化引起的寿命变化。随着寿命预测技术的发展，国内外开始以基于蠕变曲线的蠕变数据来进行寿命外推，利用短时的蠕变数据通过精确的模型预测长时的蠕变曲线从而进行精确的寿命预测，运用最多的是θ投影法，在不同断裂模式下，均可运用θ投影法处理蠕变性能数据，其优点是可以不考虑断裂模式的变化对持久寿命外推的影响，但更适用于短时蠕变寿命的预测。
近几年，从自然界事物普遍适用的热力学基本定律推导而出的蠕变断裂寿命计算方法越来越受到研究学者们的重视，此类模型具有材料参数拟合方法简单，而且一般具有较高的预测精度，是高温构件损伤评估和寿命预测的一个新的方向。
发明内容
本发明的目的是为应用在高温高压环境中的大型重要设备上的高温承压构件的蠕变断裂寿命预测提供一种新的计算方法；基于蠕变延性并考虑蠕变延性与蠕变应变速率相关性，提出的一种新的金属材料高温蠕变寿命的计算方法。
为实现本发明的目的所采用的技术方案是：
一种基于蠕变延性的高温蠕变寿命计算方法，包括以下步骤：
步骤1，获取高温条件下材料不同应力水平的蠕变性能的数据，每个试验点包括材料的应力σ(单位为MPa)、断裂时间tf(单位为h)、最小蠕变应变速率(单位为h-1)、蠕变延性εf(单位为％)、最大蠕变断裂应变εmax(单位为％)和最小蠕变断裂应变εmin(单位为％)；
步骤2，将步骤1中试验数据按照式利用数学分析软件(可采用Matlab， Origin等软件)，按最小二乘法回归，求出待定系数B和n；
步骤3，将步骤1中的试验数据按照式利用数学分析软件(可采用Matlab，Origin等软件)，按最小二乘法回归，将不同应力水平下的最小蠕变应变速率蠕变延性εf、最大蠕变断裂应变εmax、最小蠕变断裂应变εmin输入计算机，运行数学分析软件(可采用Matlab，Origin等软件)，求得待定系数α和
步骤4，根据能量守恒定律和熵守恒定律的蠕变断裂时间与应力和蠕变延性的关系式利用步骤1中的tf和σ、εf的值，采用非线性数值拟合方法(可采用Matlab， Origin等软件)求得系数K和p；
步骤5，利用寿命预测模型预测材料的高温蠕变寿命，所述寿命预测模型中的B和n由步骤2得到，所述寿命预测模型中的α和由步骤3得到，所述寿命预测模型中的K和p由步骤4得到。
在上述技术方案中，所述高温蠕变寿命计算方法适用应力水平为0.2σys-σys，其中σys为材料屈服强度。
在上述技术方案中，所述高温蠕变寿命计算方法适用作业温度为400-1200℃。
与现有技术相比，本发明的有益效果是：
1.本发明提出的寿命预测方法考虑蠕变延性与最小蠕变速率的相关性，克服了由于蠕变应变变化引起的最小蠕变速率快速衰减导致的无法实现金属材料长期寿命的准确预测的难题。
2.提高了材料蠕变寿命预测的精度，适用于各种高温金属材料。
3.预测方法简单，所需要的数据采用常规的材料蠕变性能的测试获得。
附图说明
图1是参数B和n值的拟合曲线图。
图2是蠕变延性与最小蠕变应变速率的拟合曲线，用以确定公式中和α的值。
图3是参数K和p值的拟合曲线图。
图4是基于蠕变延性的高温蠕变寿命预测方法的寿命预测曲线与试验值的比较。
具体实施方式
采用本发明提供的一种更加精确地基于蠕变延性的高温蠕变寿命预测方法，对Sanicro25钢在700℃和750℃不同应力水平下进行蠕变寿命预测分析。
第一步：获取材料不同温度、不同应力水平下蠕变性能的数据。
试验按照GB/T2039-2012《金属拉伸蠕变持久试验方法》进行。试样尺寸：直径为5mm 的标准圆棒试样，标距为50mm。试验设备为高温蠕变持久强度试验机。该试验机的组成如下：主机；加热炉；温度测控系统；变形测量系统。其负荷范围为0.3-30KN，载荷误差小于等于±1％。蠕变自动记录仪的量程为：0-10mm，测量误差不超过±0.1％。首先将试样安装在试验机上，安装好引伸计，检查试样的同轴度在规定的范围内，若超出应当按要求调整。试样安装好之后，先施加预载荷200N，并开始升温加热到预定温度然后保持60min。最后施加总载荷，记录断裂时间。在此试验中，试验温度为700℃和750℃，应力水平为:87-240Mpa。试验结束，整理得到不同温度和应力水平下材料的蠕变性能数据，包括：材料的应力σ(单位为MPa)、断裂时间tf(单位为h)、最小蠕变应变速率(单位为h-1)，蠕变延性εf(单位为％)、最大蠕变断裂应变εmax(单位为％)和最小蠕变断裂应变εmin(单位为％)。
第二步：采用诺顿函数来表示最小蠕变应变速率与应力σ的关系。由表 1中蠕变试验数据中的与σ的值，按照式利用数学分析软件(可采用Matlab， Origin等软件)，按最小二乘法回归，求出待定系数B和n。B和n参数的拟合曲线如图2 所示，B和n的值列于表2。
表1不同温度和应力水平下材料的σ、和εf数据
表2 700℃和750℃条件下分析出的B和n参数的值
第三步：按照式利用数学分析软件(可采用Matlab，Origin 等软件)，按最小二乘法回归可得到α的值。具体可将实验数据(如表1所示)中得到的最小蠕变应变速率与蠕变延性值的对应关系做成散点图(其中εmax对应εf的最大值，εmin对应εf的最小值)，对数据点在数学分析软件中进行非线性拟合，得到Sanicro25钢的蠕变延性与最小蠕变应变速率的拟合曲线，如图3所示，并可获得参数α和的值，结果如表3所示。
表3 700℃和750℃条件下分析出的参数α和的值
第四步：基于蠕变断裂时间与应力和破断蠕变量的关系式利用试验数据，采用非线性数值拟合方法(可采用Matlab，Origin等软件)求系数K和p。根据蠕变试验数据中的tf和σ、εf的值(见表4)，拟合得到如图3所示的曲线，由此可以获得参数K 与p的值，其结果列于表5。至此可以得到金属材料高温寿命预测所需的所有参数。
表4不同温度和应力水平下材料的σ、εf、tf和数据
表5 700℃和750℃条件下分析出的参数α和的值
第五步：本发明提出的考虑蠕变延性的金属材料高温蠕变寿命预测方法如式1所示：
式1中所需参数可由步骤1和步骤4获得。
将步骤1和步骤4获得参数代入公式(1)，则得到Sanicro25钢不同温度下基于蠕变延性的高温蠕变寿命预测模型。
即在700℃下：
在750℃下：
图4是利用本发明提出的蠕变寿命预测方法对Sanicro25钢的高温蠕变寿命预测结果，图中散点是文献(Chai G,Hernblom J,Peltola T,et al.Creep behavior in anewly developed heat resistant austenitic stainless steel[J].BHM Berg-undMonatshefte,2015,160(9):400-405.)报道的蠕变寿命数值，曲线是利用本发明的预测方法模拟出的曲线。和实际的Sanicro25钢高温蠕变寿命对比，实现了低应力下蠕变寿命快速衰减趋势准确模拟。可以发现本发明的基于蠕变延性的高温蠕变寿命方法，能够简便并且较高精度地计算得到蠕变寿命。
以上所述仅是本发明的优选实施方式，应当指出的是，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种基于蠕变延性的高温蠕变寿命预测方法，其特征在于，包括以下步骤：
步骤1，获取高温条件下材料不同应力水平的蠕变性能的数据，每个试验点包括材料的应力σ、断裂时间tf、最小蠕变应变速率蠕变延性εf、最大蠕变断裂应变εmax和最小蠕变断裂应变εmin；
步骤2，将步骤1中试验数据按照式按最小二乘法回归，求出待定系数B和n；
步骤3，按照式将步骤1中的不同应力水平下的最小蠕变应变速率蠕变延性εf、最大蠕变断裂应变εmax和最小蠕变断裂应变εmin，按最小二乘法回归，求得待定系数α和
步骤4，根据能量守恒定律和熵守恒定律的蠕变断裂时间与应力和蠕变延性的关系式利用步骤1中的tf和σ、εf的值，采用非线性数值拟合方法求得系数K和p；
步骤5，利用寿命预测模型预测材料的高温蠕变寿命，所述寿命预测模型中的B和n由步骤2得到，所述寿命预测模型中的α和由步骤3得到，所述寿命预测模型中的K和p由步骤4得到。
2.如权利要求1所述的一种基于蠕变延性的高温蠕变寿命预测方法，其特征在于，包括以下步骤：所述高温蠕变寿命计算方法适用应力水平为0.2σys-σys，其中σys为材料屈服强度。
3.如权利要求1所述的一种基于蠕变延性的高温蠕变寿命预测方法，其特征在于，包括以下步骤：所述高温蠕变寿命计算方法适用作业温度为400-1200℃。
CN201810388364.XA 2018-04-26 2018-04-26 一种基于蠕变延性的高温蠕变寿命预测方法 Active CN110411863B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810388364.XA CN110411863B (zh) 2018-04-26 2018-04-26 一种基于蠕变延性的高温蠕变寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810388364.XA CN110411863B (zh) 2018-04-26 2018-04-26 一种基于蠕变延性的高温蠕变寿命预测方法

Publications (2)

Publication Number Publication Date
CN110411863A true CN110411863A (zh) 2019-11-05
CN110411863B CN110411863B (zh) 2022-02-11

Family

ID=68346017

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810388364.XA Active CN110411863B (zh) 2018-04-26 2018-04-26 一种基于蠕变延性的高温蠕变寿命预测方法

Country Status (1)

Country Link
CN (1) CN110411863B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112730092A (zh) * 2020-12-24 2021-04-30 南京航空航天大学 一种基于失效模式一致的蠕变载荷等效加速方法
CN113008677A (zh) * 2021-03-02 2021-06-22 中南大学 一种镍基高温合金的蠕变持久性能预测方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60260839A (ja) * 1984-06-07 1985-12-24 Japan Steel Works Ltd:The 電気化学的劣化度測定法
CN101710053A (zh) * 2009-11-06 2010-05-19 上海师范大学 高温材料的蠕变寿命预测方法
CN102976700A (zh) * 2012-11-13 2013-03-20 云南磷化集团有限公司 用于模拟磷矿床软夹层的模型试验材料及其制备方法
CN103091164A (zh) * 2013-01-15 2013-05-08 天津大学 一种适用于微纳米薄膜材料的双系统拉伸装置
TW201348452A (zh) * 2012-05-16 2013-12-01 Nippon Steel & Sumitomo Metal Corp 金屬材料之塑性加工方法及塑性加工裝置
CN104156577A (zh) * 2014-07-31 2014-11-19 广东电网公司电力科学研究院 一种超超临界锅炉异种钢管焊接接头寿命评估方法
CN105158084A (zh) * 2015-09-15 2015-12-16 华东理工大学 一种材料的蠕变-疲劳寿命预测方法
CN106934168A (zh) * 2017-03-21 2017-07-07 中国石油大学(华东) 一种材料多轴蠕变失效应变预测方法
CN107391903A (zh) * 2017-06-07 2017-11-24 上海电力学院 用于马氏体钢的蠕变‑疲劳寿命预测及其可靠性分析方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60260839A (ja) * 1984-06-07 1985-12-24 Japan Steel Works Ltd:The 電気化学的劣化度測定法
CN101710053A (zh) * 2009-11-06 2010-05-19 上海师范大学 高温材料的蠕变寿命预测方法
TW201348452A (zh) * 2012-05-16 2013-12-01 Nippon Steel & Sumitomo Metal Corp 金屬材料之塑性加工方法及塑性加工裝置
CN102976700A (zh) * 2012-11-13 2013-03-20 云南磷化集团有限公司 用于模拟磷矿床软夹层的模型试验材料及其制备方法
CN103091164A (zh) * 2013-01-15 2013-05-08 天津大学 一种适用于微纳米薄膜材料的双系统拉伸装置
CN104156577A (zh) * 2014-07-31 2014-11-19 广东电网公司电力科学研究院 一种超超临界锅炉异种钢管焊接接头寿命评估方法
CN105158084A (zh) * 2015-09-15 2015-12-16 华东理工大学 一种材料的蠕变-疲劳寿命预测方法
CN106934168A (zh) * 2017-03-21 2017-07-07 中国石油大学(华东) 一种材料多轴蠕变失效应变预测方法
CN107391903A (zh) * 2017-06-07 2017-11-24 上海电力学院 用于马氏体钢的蠕变‑疲劳寿命预测及其可靠性分析方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
A•MEHMANPARAST: "Creep crack growth rate predictions in 316H steel using stress dependent creep ductility", 《MATERIALS AT HIGH TEMPERATURES》 *
DONGQUAN WU: "Theoretical and numerical analysis of creep crack initiation combined with primary and secondary stresses", 《THEORETICAL AND APPLIED FRACTURE MECHANICS》 *
M.T. WHITTAKER 等: "Creep and creep fracture of 2.25Cr-1.6W steels", 《MATERIALS SCIENCE AND ENGINEERING A》 *
WARWICK M.PAYTEN 等: "A strain energy density method for the prediction of creep–fatigue damage in high temperature components", 《MATERIALS SCIENCE AND ENGINEERING: A》 *
张俊善: "《材料强度学》", 31 December 2014, 哈尔滨工业大学出版社 *
李聪成 等: "蠕变疲劳交互作用下裂纹萌生的有限元模拟", 《焊接学报》 *
温建锋: "基于应变的损伤力学模型及其在蠕变裂纹扩展数值模拟中的应用", 《中国博士学位论文全文数据库基础科学辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112730092A (zh) * 2020-12-24 2021-04-30 南京航空航天大学 一种基于失效模式一致的蠕变载荷等效加速方法
US11965861B2 (en) 2020-12-24 2024-04-23 Nanjing University Of Aeronautics And Astronautics Equivalent acceleration method of creep loads based on consistent failure mode
CN113008677A (zh) * 2021-03-02 2021-06-22 中南大学 一种镍基高温合金的蠕变持久性能预测方法
CN113008677B (zh) * 2021-03-02 2022-06-10 中南大学 一种镍基高温合金的蠕变持久性能预测方法

Also Published As

Publication number Publication date
CN110411863B (zh) 2022-02-11

Similar Documents

Publication Publication Date Title
Busby et al. The relationship between hardness and yield stress in irradiated austenitic and ferritic steels
Chen et al. Experimental and modeling results of creep–fatigue life of Inconel 617 and Haynes 230 at 850° C
CN107843510B (zh) 基于室温布氏硬度预测超临界机组t/p91耐热钢剩余持久寿命评估方法
Kamaya Fatigue crack tolerance design for stainless steel by crack growth analysis
Heldt et al. Stress corrosion cracking of low-alloy, reactor-pressure-vessel steels in oxygenated, high-temperature water
CN110411863A (zh) 一种基于蠕变延性的高温蠕变寿命预测方法
CN110411864A (zh) 一种基于蠕变激活能的高温蠕变寿命预测解析计算方法
Fargione et al. The influence of the load frequency on the high cycle fatigue behaviour
Wang et al. A method for directly measuring fracture toughness and determining reference temperature for RPV steels by Charpy impact test
Zhao et al. Analysis on stress‐strain behavior and life prediction of P92 steel under creep‐fatigue interaction conditions
Solin et al. Fatigue of stainless steel in simulated operational conditions: effects of PWR water, temperature and holds
Sawada et al. Analysis of long-term creep curves by constitutive equations
Majumdar Designing against low-cycle fatigue at elevated temperature
Wilshire et al. Acquisition and analysis of creep data
Kharchenko et al. Analysis of the methods for determination of strength characteristics of NPP main equipment metal from the results of hardness and indentation measurements
Kerezsi et al. A two-stage model for predicting crack growth due to repeated thermal shock
Evans Assessing the predictive performance of creep models using absolute rather than squared prediction errors: an application to 2.25 Cr-1Mo steel and 316H stainless steel
Tao et al. Multiaxial notch fatigue life prediction based on the dominated loading control mode under variable amplitude loading
CN110008527A (zh) 一种耐热钢材料剩余寿命评估方法
JP6582753B2 (ja) 耐熱鋼材の寿命予測方法
Jaske et al. Long-life fatigue of type 316 stainless steel at temperatures up to 593 C
Le Duff et al. Effects of surface finish and loading conditions on the low cycle fatigue behavior of austenitic stainless steel in PWR environment: comparison of LCF test results with NUREG/CR-6909 life estimations
Spindler et al. Creep fatigue behaviour of type 321 stainless steel at 650 C
Hou et al. Notch Effect on Creep-Fatigue Behavior of Alloy 617 at Elevated Temperature
Evans Incorporating the Wilshire equations for time to failure and the minimum creep rate into a continuum damage mechanics for the creep strain of Waspaloy

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
CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: 300452 Binhai Industrial Research Institute Campus of Tianjin University, No. 48 Jialingjiang Road, Binhai New Area, Tianjin

Patentee after: Tianjin University

Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92

Patentee before: Tianjin University