CN109255175B - 一种基于临界再结晶残余应力的再结晶控制方法 - Google Patents
一种基于临界再结晶残余应力的再结晶控制方法 Download PDFInfo
- Publication number
- CN109255175B CN109255175B CN201811009741.0A CN201811009741A CN109255175B CN 109255175 B CN109255175 B CN 109255175B CN 201811009741 A CN201811009741 A CN 201811009741A CN 109255175 B CN109255175 B CN 109255175B
- Authority
- CN
- China
- Prior art keywords
- stress
- recrystallization
- temperature
- critical
- strain
- 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
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
-
- C—CHEMISTRY; METALLURGY
- C30—CRYSTAL GROWTH
- C30B—SINGLE-CRYSTAL GROWTH; UNIDIRECTIONAL SOLIDIFICATION OF EUTECTIC MATERIAL OR UNIDIRECTIONAL DEMIXING OF EUTECTOID MATERIAL; REFINING BY ZONE-MELTING OF MATERIAL; PRODUCTION OF A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; SINGLE CRYSTALS OR HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; AFTER-TREATMENT OF SINGLE CRYSTALS OR A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; APPARATUS THEREFOR
- C30B11/00—Single-crystal growth by normal freezing or freezing under temperature gradient, e.g. Bridgman-Stockbarger method
-
- C—CHEMISTRY; METALLURGY
- C30—CRYSTAL GROWTH
- C30B—SINGLE-CRYSTAL GROWTH; UNIDIRECTIONAL SOLIDIFICATION OF EUTECTIC MATERIAL OR UNIDIRECTIONAL DEMIXING OF EUTECTOID MATERIAL; REFINING BY ZONE-MELTING OF MATERIAL; PRODUCTION OF A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; SINGLE CRYSTALS OR HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; AFTER-TREATMENT OF SINGLE CRYSTALS OR A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; APPARATUS THEREFOR
- C30B29/00—Single crystals or homogeneous polycrystalline material with defined structure characterised by the material or by their shape
- C30B29/10—Inorganic compounds or compositions
- C30B29/52—Alloys
-
- 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/06—Power analysis or power optimisation
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Organic Chemistry (AREA)
- Physics & Mathematics (AREA)
- Metallurgy (AREA)
- Materials Engineering (AREA)
- Crystallography & Structural Chemistry (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Inorganic Chemistry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种基于临界再结晶残余应力的再结晶控制方法,通过试验的方式建立一种基于临界再结晶应力控制再结晶的模型,根据该模型结合定向凝固的数值模拟工艺,实现对再结晶缺陷的预测,从而达抑制与消除再结晶现象的目的。本发明能够实现在初始设计过程中,避免由于结构的不合理设计引起的局部应力集中,导致单晶铸件在后期的热处理中出现再结晶的现象,该技术从由结构的不合理设计引起再结晶缺陷的源头出发,从根本上解决了由于细节结构的不完善引起的再结晶,为叶片设计工作者提供一种系统解。
Description
技术领域
本发明涉及熔模铸件再结晶领域,尤其是一种残余应力的再结晶控制方法。
背景技术
镍基单晶高温合金空心涡轮叶片作为先进航空发动机中的关键热端部件,被誉为“工业之花”。由于其为净成形零件,主要通过熔模铸造的方式进行生产。然而在生产过程中存在着严重再结晶现象。
为了控制与消除再结晶缺陷,文献“申请公布号是CN103451736A中国发明专利”提出通过分级回复退火工艺对单晶高温合金铸件进行预出来,以降低在后续热处理过程中发生再结晶缺陷。文献“申请公布号是CN104550731A中国发明专利”提出通过在定位铂丝露出部分(尾部)加球形蜡件,避免单晶凝固过程中铂丝露出部分形成局部残余应力,从而降低再结晶的倾向。这些方法均是针对再结晶已存在后进行的消除或抑制措施,缺乏为叶片设计工作者提供一种系统解决在设计前期的如何选取合适的几何结构,达到预防与降低再结晶倾向的方法。因此,需要研发一种针对叶片设计中适用,将再结晶缺陷通过数值模拟的方式在设计阶段实现再结晶缺陷的抑制与消除。
发明内容
为了克服现有技术的不足,根据单晶叶片铸造过程的真实情况,通过试验的方式建立一种基于临界再结晶应力控制再结晶的模型,根据该模型结合定向凝固的数值模拟工艺,实现对再结晶缺陷的预测,从而达抑制与消除再结晶现象的目的。
为了解决缺乏叶片设计工作者提供一种系统解决在设计前期的如何选取合适的几何结构,达到预防与降低再结晶倾向的方法的问题。本发明提出了一种基于临界再结晶应力控制再结晶缺陷的模型,结合该模型通过数值模拟得到的应力值,可以预测该区域是否出现再结晶。
本发明解决其技术问题所采用的技术方案具体包括以下步骤:
步骤一:利用X衍射法对铸造试棒进行晶向检查,沿着晶粒取向应用线切割,将试棒按照Gleeble压缩试样标准,切割成高度和直径之比为1.5的圆柱,通过热拉伸试验构建不同温度及变形速率下的流变应力-应变曲线,试验中根据高温合金工作温度特征以及变形的特征,在850℃~1200℃温度区间及0.0003S-1~1S-1应变速率区间,建立2因素n水平的正交试验,即在850℃~1200℃温度区间任意取n个温度值,在0.0003S-1~1S-1应变速率区间任意取n个应变速率,根据全尺度的正交试验设计表,通过热压缩试验机建立流变应力-应变曲线;
步骤二:根据步骤1中的实验得到的流变应力-应变曲线计算峰值应力与变形温度及应变速率之间的关系如下:
步骤三:根据步骤1得到的流变应力-应变曲线,结合方程(2)计算硬化率:
式中θi——材料发生变形的硬化率;σi——压缩过程中对应的应力;εi——压缩过程中的变形量;
σp=0.9479σc+83.8544 (3)
式中σc——再结晶临界应力;
步骤七:基于临界应力与峰值应力关系结合峰值应力与变形速率及温度之间的关系模型,将公式(3)代入公式(1)中,建立再结晶临界应力与温度关系模型,选择模型适用的温度区间为(1000-1200℃):
步骤八:应用临界应力与温度关系模型结合数值模拟,对定向凝固过程中的由于几何细节结构设计不合理,引起的局部应力集中带来的再结晶问题进行预测,具体预测方法是通过数值模拟的方式,依次选择m个不同的倒圆角半径,根据倒圆角半径尺寸建立叶片的CAE分析模型,根据CAE分析模型得到的应力云图,分析倒圆角半径改变后该区域的应力场分布,取应力场中的最大值作为该结构引起的应力水平值,将不同尺寸变化后的应力水平值进行线性回归分析,建立尺寸结构与应力水平之间的回归方程,根据建立的回归方程,配合公式(4)再结晶临界应力约束条件,通过修改倒圆角的半径尺寸实现对应力的调控,完成对倒圆角结构的优化过程。
本发明的有益效果是通过试验的方式建立一种基于临界再结晶应力控制再结晶的模型,根据该模型结合定向凝固的数值模拟工艺,实现对再结晶缺陷的预测。结合该模型,能够实现在初始设计过程中,避免由于结构的不合理设计引起的局部应力集中,导致单晶铸件在后期的热处理中出现再结晶的现象,该技术从由结构的不合理设计引起再结晶缺陷的源头出发,从根本上解决了由于细节结构的不完善引起的再结晶,为叶片设计工作者提供一种系统解决在设计前期的如何选取合适的几何结构,达到预防与降低再结晶倾向的方法。
附图说明
图1是本发明所述方法的步骤流程图方法的流程图。
图2是流变应力和应变速率之间的关系图,其中图2(a)是在低应力水平下流变应力和应变速率之间的关系图;图2(b)是在高应力水平下流变应力和应变速率之间的关系图。
图4是ln[sinh(ασ)]-T-1关系图
图5是ln[sinh(ασ)]-lnZ关系图
图6是硬化率θ与应力σ关系图,其中图6(a)是850℃硬化率θ与应力σ关系图;图6(b)是980℃硬化率θ与应力σ关系图;图6(c)是1080℃硬化率θ与应力σ关系图;图6(d)是1180℃硬化率θ与应力σ关系图。
图8是不同倒圆角半径下残余应力分布对比图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
参照图1-图8。本发明一种基于临界再结晶残余应力的再结晶控制方法具体步骤如下:
步骤一:利用X衍射法对铸造试棒进行晶向检查,沿着晶粒取向应用线切割,将试棒按照Gleeble压缩试样标准,切割成高度和直径之比为1.5的圆柱,本发明切割一个Φ6mm×9mm的圆柱,通过热拉伸试验构建不同温度及变形速率下的流变应力-应变曲线,试验中根据高温合金工作温度特征以及变形的特征,在850℃~1200℃温度区间及0.0003S-1~1S-1应变速率区间,建立2因素n水平的正交试验,即在850℃~1200℃温度区间任意取n个温度值,在0.0003S-1~1S-1应变速率区间任意取n个应变速率,本发明取温度的4个水平依次为850℃,980℃,1080℃以及1180℃,应变速率的4水平依次为0.0003S-1,0.001S-1,0.01S-1以及1S-1,根据全尺度的正交试验设计表,通过热压缩试验机建立流变应力-应变曲线;
步骤二:由于应变速率与峰值应力σp在不同的应变速率下可以表达为与将不同温度条件压缩状态下的DD6单晶高温合金的峰值应力和应变速率的数据代入到两个表达式中,分别绘制出关系图,进行线性回归处理。根据拟合的系数,求解各条直线斜率倒数的平均值(n1=14.4422)。根据拟合结果,可以求出材料系数β的值,其大小为图中直线斜率倒数的平均值(β=0.02527)。由此可求出公式α=β/n1所需的α值的大小,即a=β/n1=0.0017497MPa-1。然后绘出的关系图。根据绘制关系曲线,由于其关系成线性,由此可以得到其斜率的倒数即应力指数n的取值,应力指数取各温度下应力指数的平均值,即n=10.3025。
根据余下试验数据可拟合出ln[sinh(ασ)]-1000/T关系曲线。当应变速率不变时,峰值应力的ln[sinh(ασ)]与绝对温度的倒数1000/T之间成线性关系,对其进行一元线性回归处理,其中n值为各应变速率所对应直线斜率的倒数。由此可以表明DD6单晶合金高温塑性变形时的流变应力和变形温度之间满足关系。
根据绘制ln[sinh(ασ)]-T-1曲线,通过线性回归拟合,其斜率值即为Q/1000nR,R=8.314。由于在不同的变形温度和应变速率下各直线的斜率略有所不同,因此激活能的计算值也因变形条件的不同而有所区别。根据实验数据求得不同温度及应变速率下的Q值,其平均值为862.787KJ/mol。
采用各个变形条件下的Z参数,绘制ln[sinh(ασ)]-lnZ曲线。根据ln[sinh(ασ)]-lnZ关系曲线,DD6合金峰值应力双曲正弦的自然对数和Z参数的自然对数之间满足线性关系。因此说明DD6单晶合金的流变应力与变形条件满足双曲正弦函数关系,即ln[sinh(ασ)]=0.07754lnZ-7.089。
因此可以求得:A=5.0681×1039S-1与n=12.899;结合该n值对α进行多次迭代,计算得到A,n,α:将求得的A=1.9947×1038S-1、Q=862.787KJ/mol、α=0.0063Mpa-1以及n=12.0599代入中,可以得到峰值应力与变形温度及应变速率之间的关系。
步骤三:根据流变应力应变曲线,结合方程(2)计算硬化率与应力之间的关系。依据DD6单晶合金在应变速率0.001s-1,变形温度850℃、980℃、1080℃和1180℃的应力—应变曲线,根据计算得到硬化率与应力的关系曲线。
式中θi——材料发生变形的硬化率;σi——压缩过程中对应的应力;εi——压缩过程中的变形量;
步骤五:由于铸件在设计过程中产生变形的变形速率可以忽略不计,再结合不同应变速率下(具体要包含高应变速率>1S-1,低应变速率<0.01S-1)的微观组织图片,可以得到在相同温度下,当应变速率低于0.001S-1相同热处理制度后再结晶区域的面积近似,即应变速率对再结晶的影响可以忽略,因此确定该模型应用的应变速率为一恒定值,即
步骤六:根据步骤五中确定的应变速率,结合峰值应力与发生再结晶的临界应力数据,建立再结晶临界应力与峰值应力的关系。
σp=0.9479σc+83.8544 (3)
式中σc——再结晶临界应力。
步骤七:基于临界应力与峰值应力关系结合峰值应力与变形速率及温度之间的关系模型,将公式(3)代入公式(1)中,建立再结晶临界应力与温度关系模型,由于再结晶的发生一般在1000℃以上,因此选择模型适用的温度区间为(1000-1200℃),建立再结晶临界应力与温度关系模型。
步骤八:应用临界应力与温度关系模型结合数值模拟,对定向凝固过程中的由于几何细节结构设计不合理,引起的局部应力集中带来的再结晶问题进行预测,具体预测方法是通过数值模拟的方式,依次选择m个不同的倒圆角半径,本发明取倒圆角半径分别为0.5mm,0.8mm,1.0mm,1.2mm和1.4mm,根据倒圆角半径尺寸建立叶片的CAE分析模型,根据CAE分析模型得到的应力云图,分析倒圆角半径改变后该区域的应力场分布,取应力场中的最大值作为该结构引起的应力水平值,将不同尺寸变化后的应力水平值进行线性回归分析,建立尺寸结构与应力水平之间的回归方程,根据建立的回归方程,配合公式(4)再结晶临界应力约束条件,通过修改倒圆角的半径尺寸实现对应力的调控,完成对倒圆角结构的优化过程。
应用临界应力与温度关系模型结合数值模拟,对定向凝固过程中由于榫头进气窗口由于圆角设计的不合理引起的局部应力集中,造成的再结晶问题,结合数值模拟的方式进行了优化分析,结果显示,当倒圆角为1.23mm与倒圆角为1.2mm时均在再结晶临界应力之下,由于倒圆角为1.2mm符合铸造圆角设计要求,因此确定倒圆角为1.2mm,结合实验结果,发现该区域再结晶现象得到消除,这表明基于该模型进行圆角的优化实现了对该区域的再结晶的控制。
Claims (1)
1.一种基于临界再结晶残余应力的再结晶控制方法,其特征在于包括下述步骤:
步骤一:利用X衍射法对铸造试棒进行晶向检查,沿着晶粒取向应用线切割,将试棒按照Gleeble压缩试样标准,切割成高度和直径之比为1.5的圆柱,通过热拉伸试验构建不同温度及变形速率下的流变应力-应变曲线,试验中根据高温合金工作温度特征以及变形的特征,在850℃~1200℃温度区间及0.0003S-1~1S-1应变速率区间,建立2因素n水平的正交试验,即在850℃~1200℃温度区间任意取n个温度值,在0.0003S-1~1S-1应变速率区间任意取n个应变速率,根据全尺度的正交试验设计表,通过热压缩试验机建立流变应力-应变曲线;
步骤二:根据步骤一中的实验得到的流变应力-应变曲线计算峰值应力与变形温度及应变速率之间的关系如下:
步骤三:根据步骤一得到的流变应力-应变曲线,结合方程(2)计算硬化率:
式中θi——材料发生变形的硬化率;σi——压缩过程中对应的应力;εi——压缩过程中的变形量;
σp=0.9479σc+83.8544 (3)
式中σc——再结晶临界应力;
步骤七:基于临界应力与峰值应力关系结合峰值应力与变形速率及温度之间的关系模型,将公式(3)代入公式(1)中,建立再结晶临界应力与温度关系模型,选择模型适用的温度区间为1000-1200℃:
步骤八:应用临界应力与温度关系模型结合数值模拟,对定向凝固过程中的由于几何细节结构设计不合理,引起的局部应力集中带来的再结晶问题进行预测,具体预测方法是通过数值模拟的方式,依次选择m个不同的倒圆角半径,根据倒圆角半径尺寸建立叶片的CAE分析模型,根据CAE分析模型得到的应力云图,分析倒圆角半径改变后区域的应力场分布,取应力场中的最大值作为该结构引起的应力水平值,将不同尺寸变化后的应力水平值进行线性回归分析,建立尺寸结构与应力水平之间的回归方程,根据建立的回归方程,配合公式(4)再结晶临界应力约束条件,通过修改倒圆角的半径尺寸实现对应力的调控,完成对倒圆角结构的优化过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811009741.0A CN109255175B (zh) | 2018-08-31 | 2018-08-31 | 一种基于临界再结晶残余应力的再结晶控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811009741.0A CN109255175B (zh) | 2018-08-31 | 2018-08-31 | 一种基于临界再结晶残余应力的再结晶控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109255175A CN109255175A (zh) | 2019-01-22 |
CN109255175B true CN109255175B (zh) | 2021-04-02 |
Family
ID=65049833
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811009741.0A Active CN109255175B (zh) | 2018-08-31 | 2018-08-31 | 一种基于临界再结晶残余应力的再结晶控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109255175B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110542503A (zh) * | 2019-08-29 | 2019-12-06 | 西安理工大学 | 一种高温合金表面残余应力的表征方法 |
CN111445960B (zh) * | 2020-05-18 | 2023-02-21 | 上海电机学院 | 14Cr17Ni2钢的锻造工艺参数的优化方法 |
CN113702613B (zh) * | 2021-08-19 | 2024-02-13 | 鞍钢股份有限公司 | 一种确定材料发生动态再结晶临界条件的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102519801B (zh) * | 2011-12-16 | 2015-07-29 | 南京航空航天大学 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
CN104928605A (zh) * | 2015-07-20 | 2015-09-23 | 中南大学 | 一种预测镍基合金高温流变应力和动态再结晶行为的方法 |
-
2018
- 2018-08-31 CN CN201811009741.0A patent/CN109255175B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102519801B (zh) * | 2011-12-16 | 2015-07-29 | 南京航空航天大学 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
CN104928605A (zh) * | 2015-07-20 | 2015-09-23 | 中南大学 | 一种预测镍基合金高温流变应力和动态再结晶行为的方法 |
Non-Patent Citations (4)
Title |
---|
Recrystallization and fatigue fracture of single crystal turbine blades;Y.H.He等;《Engineering Failure Analysis》;20101204;第18卷;第944-949页 * |
单晶再结晶临界应力与涡轮叶片结构设计参数间的映射模型;李世峰等;《南京航空航天大学学报》;20141231;第46卷(第6期);第909-914页 * |
基于残余应力的单晶叶片榫头进气窗结构优化;邱飞等;《哈尔滨工程大学学报》;20170228;第38卷(第2期);第235-240页 * |
镍基单晶叶片制造技术及再结晶研究进展;卜昆等;《航空制造技术》;20161125;第516卷(第21期);第34-40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109255175A (zh) | 2019-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109255175B (zh) | 一种基于临界再结晶残余应力的再结晶控制方法 | |
CN109740205B (zh) | 基于耦合晶界扩散相场法的多元合金固态相变的模拟方法 | |
US20160279878A1 (en) | Rapid synthetic material prototyping process | |
Fatemi-Varzaneh et al. | Deformation homogeneity in accumulative back extrusion processing of AZ31 magnesium alloy | |
CN109855958B (zh) | 一种金属材料拉伸性能的预测方法 | |
DE102010048969B4 (de) | Systeme und Verfahren zum rechnergestützten Vorhersagen von Wärmeübergangskoeffizienten während des Abschreckens und zum Abschrecken eines Gussteils | |
KR102028502B1 (ko) | 세로방향에서 상이한 두께를 가진 판재의 압연방법 | |
CN104504185A (zh) | 一种基于开裂准则和有限元优化钛合金开坯锻造工艺方法 | |
Tang et al. | Behavior and modeling of microstructure evolution during metadynamic recrystallization of a Ni-based superalloy | |
CN107220485B (zh) | 一种适用于多道次压缩的本构模型的建立方法 | |
CN107937842B (zh) | 一种Al-Zn-Mg-Cu铝合金的热处理方法 | |
CN110929438B (zh) | 一种考虑应变率效应的本构曲线耦合外延方法 | |
JPH0647463A (ja) | 時効成形用治具輪郭の決定方法 | |
CN115323298A (zh) | 一种镍基高温合金锻造过程的物理模拟方法 | |
CN115015318A (zh) | 一种大型构件热锻全流程宏微观分析方法及平台 | |
CN114912307A (zh) | 一种镍基单晶高温合金的高周疲劳寿命预测方法 | |
CN106202937B (zh) | M50钢锻造组织中碳化物尺寸预测方法 | |
Hudak et al. | On comparing the shape parameters of two Weibull distributions | |
CN104174803B (zh) | 高压共轨预锻模具 | |
CN111879910A (zh) | 一种优化锻造工艺参数与组织性能的试验方法 | |
CN105527309A (zh) | 一种测试热处理工艺对2024铝合金力学性能影响的方法 | |
CN110253221A (zh) | 一种用于连杆毛坯去应力的方法 | |
Vorel et al. | Material-technological modelling of controlled cooling of closed die forgings from finish forging temperature | |
CN110018046B (zh) | 一种节约型双相不锈钢中trip效应致塑性增量的表征方法 | |
CN113673055B (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 |