CN106934109B - 一种Ni3Al基合金代表体积单元模型的构建方法 - Google Patents

一种Ni3Al基合金代表体积单元模型的构建方法 Download PDF

Info

Publication number
CN106934109B
CN106934109B CN201710089897.3A CN201710089897A CN106934109B CN 106934109 B CN106934109 B CN 106934109B CN 201710089897 A CN201710089897 A CN 201710089897A CN 106934109 B CN106934109 B CN 106934109B
Authority
CN
China
Prior art keywords
phase
rve
strengthening
following
steps
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
CN201710089897.3A
Other languages
English (en)
Other versions
CN106934109A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710089897.3A priority Critical patent/CN106934109B/zh
Publication of CN106934109A publication Critical patent/CN106934109A/zh
Application granted granted Critical
Publication of CN106934109B publication Critical patent/CN106934109B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种Ni3Al基合金代表体积单元模型的构建方法,包括如下步骤:1)通过扫描电子显微镜观测合金微细观结构:2)分析统计合金微细观结构微观几何特征参量的分布规律:3)设置RVE中强化相的约束条件:4)生成RVE中强化相关键点的坐标:5)生成RVE。本发明的有益效果:本发明所生成的RVE包含了足够的微观信息,涵盖了材料的微观不确定性,可以有效的描述材料内部的力学性能,运用本专利生成的RVE开展力学分析可以更有效的预测材料受载后的力学响应,具有一定的科学理论意义和工程应用价值。

Description

一种Ni3Al基合金代表体积单元模型的构建方法
技术领域
本发明属于航空材料技术领域,具体涉及航空发动机用定向凝固高温合金材料微细观代表体积单元(RVE)模型的构建,尤其是一种Ni3Al基合金代表体积单元(RVE) 模型的构建方法。
背景技术
在研究Ni3Al基合金内部γ′/γ结构时,研究者通常都是以单胞的形式进行研究,但实际上材料内部强化相γ′相的分布并不是完全有序的,而采用单胞进行分析就忽略了强化相分布的不均匀性对材料性能带来的影响。代表体积单元由大量的强化相和基体相组成,其尺度比单胞的尺度更高,包含了更多的微观信息,为考虑材料微观信息的不确定性,可以基于强化相的大小和形状特征参量的统计分布规律来确定各强化相的大小和形状分布,对代表性体积单元(RVE)进行分析,可以更有效的描述材料的力学性能。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种Ni3Al基合金代表体积单元(RVE)模型的构建方法,能够更加有效的预测材料在受载后的力学响应。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种Ni3Al基合金代表体积单元模型的构建方法,包括如下步骤:
1)运用扫描电子显微镜,获取材料的微观结构观测图;
2)分析统计合金微细观结构微观几何特征参量的分布规律;
3)设置RVE中强化相的约束条件;
4)生成RVE中强化相关键点的坐标;
5)生成RVE。
进一步的,步骤2)的具体方法为:运用强化相体积分数wp、强化相大小lp以及强化相形状参数np三个微观几何特征参量来描述材料的微观结构观测图。
进一步的,所述强化相的体积分数wp的具体方法为:
采用颜色识别的方法对Ni3Al基合金的微观结构图进行识别,通过强化相γ′相与基体γ相之间颜色的差异对强化相γ′进行识别,截面上强化相的面积分数Af为:
其中,a、b分别为图片的长度与宽度,Ai为各个强化相的面积。
强化相的体积分数wp等于各个截面强化相面积分数的平均值:
其中,为第i个截面的强化相面积分数。
进一步的,所述强化相的大小lp的具体方法为:
以椭圆的面积计算公式以及la和lb定义强化相的等效长度为:
其中,la为图形的最大弗里特长度,lb为图形的最小弗里特长度。
进一步的,所述强化相形状参数np为:
其中,B为某一方向强化相的最大宽度,A为边界处这一方向近似直线段的宽度,下标1和2指的是两个方向。
进一步的,所述强化相大小lp和强化相形状参数np的微观概率分布形式根据统计结果以及MATLAB软件确定。
进一步的,步骤3)中RVE中强化相的约束条件为:以RVE的左下角为原点,三条边为坐标轴建立笛卡尔坐标系,设RVE的总尺寸为a×a×a,第i个强化相中心Oi的坐标为(xi,yi,zi),最大宽度为li;第j个强化相的中心Oj坐标为(xj,yj,zj),最大宽度为lj;那么,生成的强化相需要满足以下几个条件:
[1]强化相的大小li、位置(xi,yi,zi)以及形状参数ni满足微观概率模型统计出的分布规律;
[2]强化相的位置不能超出RVE的范围;
[3]各个强化相之间不能出现相互干涉的问题;
[4]对于本发明的镍基合金,根据试验测定,体积分数处于40%~52%之间强化相体积含量wp满足:
wpmin≤wp≤wpmax
其中,wpmin=40%、wpmax=52%。
进一步的,令相邻的强化相以及RVE边界之间的间距值为d,分别满足以下条件:
不超出RVE的边界:
各个强化相相互之间不干涉:
进一步的,步骤4)的具体方法为:根据所得的材料微观几何特征参数以及RVE中强化相的约束,通过MATLAB软件生成各强化相的关键点坐标,并将关键点的坐标写入指定的文件。
有益效果:本发明提供的Ni3Al基合金代表体积单元(RVE)模型的构建方法,与现有技术比较,具有以下优势:本发明所生成的RVE包含了足够的微观信息,涵盖了材料的微观不确定性,可以有效的描述材料内部的力学性能,运用本专利生成的RVE 开展力学分析可以更有效的预测材料受载后的力学响应,具有一定的科学理论意义和工程应用价值。
附图说明
图1为代表体积单元的构建流程图;
图2为IC10合金的微观结构图;
图3为IPP软件对微观结构图进行颜色识别的结果示意图;
图4为不同数目的截面时强化相面积分数平均值随截面数目的变化情况示意图;
图5为强化相大小定义方式示意图;
图6为形状参数的定义方法示意图;
图7为强化相形状参数的测量方法示意图;
图8为完全随机的RVE示意图;
图9为组合单胞式RVE的效果图;
图10为三维及二维RVE几何模型示意图;
图11为两张试验结果图。
具体实施方式
本发明为一种Ni3Al基合金代表体积单元(RVE)模型的构建方法,包括如下步骤:1)通过扫描电子显微镜(SEM)观测合金微细观结构:2)分析统计合金微细观结构微观几何特征参量的分布规律:3)设置RVE中强化相的约束条件:4)生成RVE中强化相关键点的坐标:5)生成RVE。本发明的有益效果:本发明所生成的RVE包含了足够的微观信息,涵盖了材料的微观不确定性,可以有效的描述材料内部的力学性能,运用本专利生成的RVE开展力学分析可以更有效的预测材料受载后的力学响应,具有一定的科学理论意义和工程应用价值。
下面结合附图和实施例对本发明作更进一步的说明。
实施例
1.一种Ni3Al基合金代表体积单元(RVE)模型的构建方法,如图1所示,包括如下步骤:
1),运用扫描电子显微镜,获取材料的微观结构观测图,如图2所示为IC10合金的微观结构图;
2),通过强化相体积分数wp、强化相大小lp以及强化相形状参数np三个微观几何特征参量来描述材料的微观结构观测图。
2-1),强化相的体积分数(wp)
在ImagePro plus(IPP)软件中,图片中需要进行的测量的区域被称为AOI(areaof interest)即感兴趣区域,颜色识别是其最有效的方法之一。本文就是采用的颜色识别的方法对Ni3Al基合金的微观结构图进行识别,通过强化相γ′相与基体γ相之间颜色的差异对强化相γ′进行识别。
图3即为IPP软件对IC10合金γ′相进行颜色识别的结果,在其测量模块选择测量每个块的面积。同时运用手动测量工具测量整个图片的长度a与宽度b,则此截面上强化相的面积分数可以用下式求出:
实际上材料都是三维结构,而通过SEM实验只能计算到二维表面上的面积分数,本文通过微元分析的方法找到了三维体积分数与二维面积分数之间的关系,发现强化相的体积分数可以用多个截面的面积分数平均值来代替。
首先在材料内部取一个立方体代表体积单元(RVE),V为RVE体积,A为RVE截面面积,l为RVE边长,将RVE沿其中一个方向切为厚度相同的n份薄片。
假设当n足够大时,薄片在厚度方向上强化相的形状分布不变,则第i个薄片强化相的体积Vi等于第i个薄片前表面强化相面积Api乘以微元段的厚度Δl,即:
Vi=ΔlApi (2)
强化相体积分数可以用下式表示:
其中Vp为强化相体积。同时可知:
第i薄片前表面上强化相面积分数为:
将式(4)和5)代入式(6)可得:
由(6)可以看出,当n足够大时,强化相体积分数等于各个截面强化相面积分数的平均值。
下一步将找到一个合适的n,使得切片数大于n时,式(6)成立,拟采用蒙特卡洛的概率方法。
首先,建立合适的代表体积单元。建立的RVE大小为1,内部含有18个随机生成的直径为0.3的球状强化相,强化相的体积分数wp为25.43%。
接着,选择指定的截面,本文的切片方向选择为z方向(即z坐标固定),x,y坐标按照均匀分布产生,共生成10000个点,将生成的点i(x,y,z)进行判断,统计落在强化相内部点的数目m,则这个面上强化相面积分数就为m/10000。
第三,选取一定的截面数目,计算面积分数的平均值。
最后,对不同数目的截面进行10次操作,得到10个强化相面积分数平均值。
为了进一步验证式(5)的实用性,本文在选取截面坐标时也采取生成随机数的方式,图4给出了选取不同数目的截面时强化相面积分数平均值随截面数目的变化情况。其中中间黑色实线为强化相的体积分数,上下两条黑色实线线为强化相体积分数10%误差带(22.887%~27.973%),散点为各数目下强化相面积分数的平均值,其中每个数目下共进行了10组随机取样。
从图4可以看出,当选取的有效截面超过18时,强化相面积分数的平均值与体积分数之间的差值就基本在±2.54%上下浮动,并且图4的算例中强化相体积分数仅为25%,随着强化相体积分数的增加,由于强化相体积的增加,这之间的差值还会进一步减少。这也说明强化相的体积分数可以用多个截面(视场)的强化相面积分数的平均值表示。
2-2),强化相的大小(lp)
由图2可以看出,强化相的形状都是不规则的,并且由于IPP软件在进行颜色识别时存在一定误差,如果用强化相面积直接换算大小会产生较大的误差,为了更好的描述强化相的大小,本文综合考虑IPP软件的测量功能以及IC10合金强化相的实际形状,为强化相的大小lp给出了如下的定义方式。
假设强化相的形状为图5显示的不规则形状,由于其形状复杂,不能用一个参数来确定其大小,但在后期的建模过程中所需生成的强化相数目较多,所以需要简化参数数目,用尽可能少的参数去形容强化相的大小。
在IPP软件中,其可以自己测量观测物最大的以及最小的弗里特长度(Feretlength),如图5所示,la即为图形的最大弗里特长度,lb为图形的最小弗里特长度。
从图2可以看出,实际结构中IC10合金的强化相基本上都是近似椭圆的形状,本文以椭圆的面积计算公式以及la和lb定义强化相的等效长度为:
2-3),强化相形状参数(np)
图6给出了强化相形状参数np的定义方法,其中B为强化相的最大宽度,A为此方向强化相边缘处直线段的长度,强化相的形状参数就为A与B的比值;图7给出了np的测量方法,首先按照图6所示的强化相形状参数的定义方法将强化相按照图7(a)的方法进行分割,测得强化相最大宽度以及此方向边界处近似直线段长度。图7(b)给出了实际测量过程中的效果图。在实际测量过程中,本文选择了强化相的两个方向进行测量参数 A和B,最后强化相形状参数np的值为:
其中下标1和2指的是两个方向。
2-4),强化相大小lp和强化相形状参数np的微观概率分布形式根据统计结果以及MATLAB软件确定。
3),RVE中强化相的约束条件。图8给出了一个完全随机RVE的示意图,为方便叙述,以RVE的左下角为原点,三条边为坐标轴建立笛卡尔坐标系。
设RVE的总尺寸为a×a×a,第i个强化相中心Oi的坐标为(xi,yi,zi),最大宽度为li;第j个强化相的中心Oj坐标为(xj,yj,zj),最大宽度为lj
根据微观概率模型和几何原则,生成的强化相需要满足以下几个条件:
[1]强化相的大小li、位置(xi,yi,zi)以及形状参数ni满足微观概率模型统计出的分布规律。
[2]强化相的位置不能超出RVE的范围。
[3]各个强化相之间不能出现相互干涉的问题。
[4]需要满足一定的强化相体积含量wp
同时为了防止相邻的强化相以及RVE边界之间的间距太小,造成后续有限元网格畸形的问题,所以设置了一个合适的间距值d。将上述的几个条件用数学公式进行如下的表述。
不超出RVE的边界:
各个强化相相互之间不干涉:
满足一定的强化相体积含量:
wpmin≤wp≤wpmax (15)
其中wpmin=40%、wpmax=52%。
当强化相的体积分数和强化相的数目达到一定程度,新的强化相生成的约束情况会急剧生成,导致RVE的生成难度大大上升。图9给出了一种组合单胞式的RVE生成效果图,从图9可以看出组合单胞式的RVE忽略了强化相位置的随机性,每个强化相的位置都被限制在了每个立方体单胞内。虽然在组合单胞式的RVE中,强化相位置的随机性未被考虑,但由于其位置被限制在各自的立方体单胞内,也使得各个强化相之间不会出现相互干涉的情况,所以这种生成RVE方法在强化相体积分数高以及强化相数量多的情况下适用。
组合单胞式RVE的生成方法与完全随机法类似,但是由于其强化相位置被限定,所以相比于完全随机的生成方法,其少了位置约束条件,这极大的减小了有效随机数的生成难度,使得其在强化相体积分数高以及强化相数量多的情况下表现良好。
4),根据上文给出的材料微观几何特征参数以及RVE中强化相的约束,通过MATLAB软件生成各强化相的关键点坐标。将关键点的坐标写入指定的文件。
5),通过有限元软件ABAQUS读取关键点的信息,生成RVE。
图10为通过本文所描述方法建立的二维(图10(b))及三维(图10(a))RVE的几何模型。
本发明以镍基合金晶粒内部多相组织结构为侧重点,运用三个微观几何特征参数表征了材料的微细观结构,并通过统计分析得到了微观几何特征参数的大小及分布规律。随后根据得到的微观几何特征参数通过完全随机以及组合单胞两种方法建立材料的微观代表体积单元。通过三个微观几何特征参数架构了材料实际的微观结构与用于分析的代表体积单元之间的关系,最大程度的还原了材料的微观结构。
为了进一步说明本发明的的有效性,如图11所示,取两张试验所得微观结构图与二维RVE进行对比。表1给出了对比结果:
表1试验结果与RVE的对比
从表1可以看出,最终建立的RVE三个几何特征参数与原试验结果的十分想近,并且由于强化相分布不均匀,误差也是不可避免的。表2给出了最终建立的RVE与试验结果的误差分析。从表2可以看出最终建立的2为RVE与两个试验结果的误差均不大,这也说明本发明所提出的RVE构建方法的有效性。
表2误差分析
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (2)

1.一种Ni3Al基合金代表体积单元模型的构建方法,其特征在于:包括如下步骤:
1)运用扫描电子显微镜,获取材料的微观结构观测图;
2)分析统计合金微观结构微观几何特征参量的分布规律;
3)设置RVE中强化相的约束条件;
4)生成RVE中强化相关键点的坐标;
5)生成RVE;
步骤2)的具体方法为:运用强化相体积分数wp、强化相大小lp以及强化相形状参数np三个微观几何特征参量来描述材料的微观结构观测图,具体方法为:
2-1),所述强化相的体积分数wp的具体方法为:
采用颜色识别的方法对Ni3Al基合金的微观结构图进行识别,通过强化相γ′相与基体γ相之间颜色的差异对强化相γ′进行识别,截面上强化相的面积分数为:
其中,a、b分别为图片的长度与宽度,Ai为各个强化相的面积;
强化相的体积分数wp等于各个截面强化相面积分数的平均值:
其中,为第i个截面的强化相面积分数;
2-2),所述强化相的大小lp的具体方法为:
以椭圆的面积计算公式以及la和lb定义强化相的等效长度为:
其中,la为图形的最大弗里特长度,lb为图形的最小弗里特长度;
2-3),所述强化相形状参数np为:
其中,B1、B2为某一方向强化相的最大宽度,A1、A2为边界处这一方向近似直线段的宽度;
2-4),所述强化相大小lp和强化相形状参数np的微观概率分布形式根据统计结果以及MATLAB软件确定,步骤3)中RVE中强化相的约束条件为:以RVE的左下角为原点,三条边为坐标轴建立笛卡尔坐标系,设RVE的总尺寸为a×a×a,第i个强化相中心Oi的坐标为(xi,yi,zi),最大宽度为li;第j个强化相的中心Oj坐标为(xj,yj,zj),最大宽度为lj;那么,生成的强化相需要满足以下几个条件:
[1]强化相的大小li、位置(xi,yi,zi)以及形状参数ni满足微观概率模型统计出的分布规律;
[2]强化相的位置不能超出RVE的范围;
[3]各个强化相之间不能出现相互干涉的问题;
[4]对于本发明的镍基合金,根据试验测定,体积分数处于40%~52%之间强化相体积分数wp满足:
wpmin≤wp≤wpmax
其中,wpmin=40%、wpmax=52%;
令相邻的强化相以及RVE边界之间的间距值为d,生成的强化相满足以下条件:
强化相的位置不超出RVE的边界:
各个强化相相互之间不干涉:
2.根据权利要求1所述的Ni3Al基合金代表体积单元模型的构建方法,其特征在于:步骤4)的具体方法为:根据所得的材料微观几何特征参数以及RVE中强化相的约束,通过MATLAB软件生成各强化相的关键点坐标,并将关键点的坐标写入指定的文件。
CN201710089897.3A 2017-02-20 2017-02-20 一种Ni3Al基合金代表体积单元模型的构建方法 Active CN106934109B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710089897.3A CN106934109B (zh) 2017-02-20 2017-02-20 一种Ni3Al基合金代表体积单元模型的构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710089897.3A CN106934109B (zh) 2017-02-20 2017-02-20 一种Ni3Al基合金代表体积单元模型的构建方法

Publications (2)

Publication Number Publication Date
CN106934109A CN106934109A (zh) 2017-07-07
CN106934109B true CN106934109B (zh) 2019-10-11

Family

ID=59423356

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710089897.3A Active CN106934109B (zh) 2017-02-20 2017-02-20 一种Ni3Al基合金代表体积单元模型的构建方法

Country Status (1)

Country Link
CN (1) CN106934109B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108256258B (zh) * 2018-01-31 2019-01-25 西安科技大学 基于sem图像的胶结充填体力学响应特性预测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103714216A (zh) * 2013-12-31 2014-04-09 北京理工大学 基于三维微观晶须构成的复合材料力学性能评估方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103714216A (zh) * 2013-12-31 2014-04-09 北京理工大学 基于三维微观晶须构成的复合材料力学性能评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Study on a statistical unit cell model for Ni3Al-base superalloy;Hongjian Zhang, Jianfeng Xiao, Weidong Wen, etc.;《Mechanics of Materials》;20160407;全文 *
多晶镍基合金循环塑性细观本构关系;蔚夺魁,杨晓光,张克实;《航空动力学报》;20131031;全文 *

Also Published As

Publication number Publication date
CN106934109A (zh) 2017-07-07

Similar Documents

Publication Publication Date Title
Quey et al. Large-scale 3D random polycrystals for the finite element method: Generation, meshing and remeshing
Hsu et al. Direct immersogeometric fluid flow analysis using B-rep CAD models
Falco et al. Generation of 3D polycrystalline microstructures with a conditioned Laguerre-Voronoi tessellation technique
Feather et al. A numerical study into element type and mesh resolution for crystal plasticity finite element modeling of explicit grain structures
Saylor et al. Statistically representative three-dimensional microstructures based on orthogonal observation sections
Rollett et al. Stress hot spots in viscoplastic deformation of polycrystals
Latypov et al. Application of chord length distributions and principal component analysis for quantification and representation of diverse polycrystalline microstructures
Rollett et al. Three-dimensional characterization of microstructure by electron back-scatter diffraction
Cairney et al. Mining information from atom probe data
Hallberg et al. Modeling of grain growth under fully anisotropic grain boundary energy
Qian et al. Quality improvement of non‐manifold hexahedral meshes for critical feature determination of microstructure materials
Li et al. Topological and atomic scale characterization of grain boundary networks in polycrystalline and nanocrystalline materials
CN103714216B (zh) 基于三维微观晶须构成的复合材料力学性能评估方法
JP2008197852A (ja) 塑性加工における被加工材の組織発展の解析装置、解析システムおよび記録媒体
US20020042697A1 (en) Mesh generation system, design support system, analysis system, analysis method, mesh generation method, and storage medium and program transmission apparatus therefor
CN111539139A (zh) 一种颗粒随机分布的复合材料2d细观结构建模方法
Falco et al. A new method for the generation of arbitrarily shaped 3D random polycrystalline domains
CN106934109B (zh) 一种Ni3Al基合金代表体积单元模型的构建方法
Ghosh et al. Statistically equivalent representative volume elements (SERVE) for material behaviour analysis and multiscale modelling
Kang et al. PolyProc: A modular processing pipeline for X-ray diffraction tomography
Critchfield et al. Representative and statistical volume elements for grain boundary networks: A stereological approach
Li et al. A texture evolution model in cubic-orthotropic polycrystalline system
CN115356703B (zh) 基于面元分布的粗糙目标rcs缩比测量方法及装置
US10371650B2 (en) Macrotexture map visualizing texture heterogeneity in polycrystalline parts
Ghosh et al. A domain partitioning based pre-processor for multi-scale modelling of cast aluminium alloys

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