【技术实现步骤摘要】
本专利技术涉及多晶陶瓷先微结构分析，尤其是涉及一种考虑晶粒微结构的结构陶瓷本征强度模拟方法。
技术介绍
1、众所周知，脆性是陶瓷材料典型特点，即陶瓷材料对缺陷非常敏感，阐明微观结构对陶瓷材料强度的影响是材料研发的关键。
2、由于在实际生产中，缺陷种类众多，很难用单一的模型分析所有种类的缺陷，压痕在以前的研究中被广泛用于模拟初始缺陷裂纹。研究结果表明，当缺陷尺寸较大时，裂纹尺寸与抗弯强度的对数分布关系是一条直线。随着压痕裂纹的减少，强度和压痕裂纹尺寸之间的关系变为非线性关系，直到达到强度对初始裂纹不再敏感的平台。基于这一实验可以发现，材料的强度和缺陷尺寸息息相关，那么多大的缺陷对应的强度为材料的本征强度。由于材料检测放大的局限，实际材料中的所有缺陷很难全部检测得到。因此，材料强度的定义似乎并不严格，材料的本征强度是多少以及初始裂纹如何影响强度测试值均不可知。根据断裂力学中强度的定义(断裂韧性)，材料断裂的临界状态是裂纹不稳定扩展，因此，与裂纹的临界不稳定扩展相对应的强度才是材料的本征强度。
技术实现思路
1、本专利技术的目的是提供一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，通过计算裂纹扩展阻力研究材料微观结构与本征断裂强度之间的关系，为高强度、高韧性结构陶瓷的本征强度调控提供依据。
2、为实现上述目的，本专利技术提供了一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，包括以下步骤：
3、s1、建立双尺度模型，双尺度模型包括微结构模型和宏观模型；>
4、s2、计算微结构模型的损伤，获取材料的损伤力学参数；
5、s3、将微结构模型的损伤力学参数传递给宏观模型，通过宏观模型的断裂，计算材料中裂纹扩展阻力曲线；
6、s4、通过裂纹扩展阻力曲线计算材料的本征强度。
7、优选的，所述s1中双尺度模型中的微结构模型的建立包括以下步骤：
8、步骤一、利用python生成voronoi多边形，voronoi多边形用于晶粒模拟，voronoi多边形的边用于晶界模拟；
9、步骤二、通过软件abaqus对微结构模型划分网格；为实现任意裂纹扩展，基于abaqus二次开发接口，给已划分好网格的模型生成覆盖全模型的零厚度内聚力单元，在内聚力单元损伤起始和演化前起遵循线弹性响应关系，输入材料力学参数并施边界条件，通过静态分析进行模拟拉伸、剪切测试，得到微观模型的损伤力学参数。
10、优选的，所述s2中的损伤力学参数包括在拉伸和剪切应力作用下的拉伸强度、剪切强度、拉伸断裂能、剪切断裂能和材料参数。
11、优选的，所述s1中双尺度模型中的宏观模型为具有预制裂纹的宏观有限元模型，其建立包括以下步骤：建立长宽分别为40mm*4mm的三点抗弯有限元模型，为模型布置全局网格种子点，生成合适大小的网格。
12、优选的，所述宏观模型为了实现任意裂纹扩展，扩展有限元技术用于模拟裂纹的扩展，损伤起始采用最大正应力准则，损伤演化采用线性退化准则，具体参考微观模型的最大正应力准则和线性退化准则；其输入参数来源于微观代表性体积元的损伤力学参数，对宏观模型输入材料力学参数并施加载荷，赋予正确的网格类型，并通过静态分析进行三点抗弯工况，得到裂纹开始扩展时的载荷。
13、优选的，所述s3中的裂纹扩展的阻力曲线基于具有预制裂纹的宏观有限元模型裂纹开始扩展时的载荷，根据astm标准中的三点抗弯公式计算应力强度因子，通过不同初始裂纹长度的应力强度因子，从而得到裂纹扩展的阻力曲线，裂纹扩展的阻力曲线计算公式如下：
14、
15、其中kmax是应力强度因子饱和值，δkmax是由于裂纹稳定扩展导致的断裂韧性增长量，a0是初始裂纹长度，λ是与增加的裂纹扩展阻力相关的因子。
16、优选的，所述宏观模型的裂纹扩展过程中，当施加在裂纹尖端的应力强度因子与裂纹扩展阻力相等，同时裂纹尖端的应力强度因子与裂纹扩展阻力随着裂纹尺寸增加的速率相等时，被视为裂纹失稳扩展的临界条件，不稳定裂纹扩展条件通过如下关系进行评估：
17、
18、其中a是裂纹长度，kr为裂纹扩展阻力曲线，kapp为施加在裂纹尖端的应力强度因子，由下式给出：
19、
20、其中σ是施加的应力，y是几何因子，裂纹的几何因子参考：
21、
22、q＝1+1.464(2b/a)1.65 (5)
23、
24、h＝1-(1.22+0.24b/a)b/h+[0.55-1.05(2b/a)0.75+0.47(2b/a)1.5](b/h)2 (7)
25、其中b是裂纹宽度，h为样品高度，q为裂纹形状因子，m为表面修正因子，h为用于修正应力强度因子的多项式。
26、优选的，所述s4中计算材料的本征强度时联合公式(1)-(3)得：
27、
28、
29、所以：
30、
31、因为远小于和所以：
32、
33、当初始裂纹小于λ/2时，强度随着裂纹的扩展而增加，当裂纹扩展达到λ/2时，强度达到其最大值，如果裂纹需要进一步稳定扩展，则需要减少载荷才能使裂纹再次稳定，如果载荷继续增加，裂纹将发生不稳定扩展，如果初始裂纹大于λ/2，随着载荷的增加，裂纹不会稳定扩展，直接进入不稳定扩展阶段，最大强度即本征强度为λ和y的函数：
34、
35、因此，本专利技术采用上述一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，具有以下有益效果：
36、(1)通过建立双尺度模型，定量分析了显微结构对断裂韧性的影响，解决了传统微尺度模型由于预制裂纹过短、裂纹扩展长度有限、参与计算晶粒有限等缺点导致的断裂韧性计算不准确的问题，同时也解决了宏观模型与显微结构脱离的问题。
37、(2)通过计算可得到材料的本征强度，解决了实验中由于初始缺陷大小难以控制导致的强度测试不准确的问题，为多晶陶瓷材料的研发、制备、应用提供了一种快速评估的方法，起到省时省力、快速高效且准确的作用效果。
38、下面通过附图和实施例，对本专利技术的技术方案做进一步的详细描述。
本文档来自技高网...【技术保护点】
1.一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：包括以下步骤：
2.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述S1中双尺度模型中的微结构模型的建立包括以下步骤：
3.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述S2中的损伤力学参数包括在拉伸和剪切应力作用下的拉伸强度、剪切强度、拉伸断裂能、剪切断裂能和材料参数。
4.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述S1中双尺度模型中的宏观模型为具有预制裂纹的宏观有限元模型，其建立包括以下步骤：建立长宽分别为40mm*4mm的三点抗弯有限元模型，为模型布置全局网格种子点，生成合适大小的网格。
5.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述宏观模型为了实现任意裂纹扩展，扩展有限元技术用于模拟裂纹的扩展，损伤起始采用最大正应力准则，损伤演化采用线性退化准则，具体参考微观模型的最大正应力准则和线性退化准则；其输入参数来源于微观代
6.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述S3中的裂纹扩展的阻力曲线基于具有预制裂纹的宏观有限元模型裂纹开始扩展时的载荷，根据ASTM标准中的三点抗弯公式计算应力强度因子，通过不同初始裂纹长度的应力强度因子，从而得到裂纹扩展的阻力曲线，裂纹扩展的阻力曲线计算公式如下：
7.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述宏观模型的裂纹扩展过程中，当施加在裂纹尖端的应力强度因子与裂纹扩展阻力相等，同时裂纹尖端的应力强度因子与裂纹扩展阻力随着裂纹尺寸增加的速率相等时，被视为裂纹失稳扩展的临界条件，不稳定裂纹扩展条件通过如下关系进行评估：
8.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述S4中计算材料的本征强度时联合公式(1)-(3)得：
...【技术特征摘要】
1.一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：包括以下步骤：
2.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述s1中双尺度模型中的微结构模型的建立包括以下步骤：
3.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述s2中的损伤力学参数包括在拉伸和剪切应力作用下的拉伸强度、剪切强度、拉伸断裂能、剪切断裂能和材料参数。
4.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述s1中双尺度模型中的宏观模型为具有预制裂纹的宏观有限元模型，其建立包括以下步骤：建立长宽分别为40mm*4mm的三点抗弯有限元模型，为模型布置全局网格种子点，生成合适大小的网格。
5.根据权利要求1所述的一种考虑晶粒微结构的结构陶瓷本征强度模拟方法，其特征在于：所述宏观模型为了实现任意裂纹扩展，扩展有限元技术用于模拟裂纹的扩展，损伤起始采用最大正应力准则，损伤演化采用线性退化准则，具体参考微观模型的最大正应力准则和线性退化...
【专利技术属性】
技术研发人员：关康，崔金平，饶平根，曾庆丰，刘建涛，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。