CN112898032A - 一种凝胶注模陶瓷生坯热脱脂速率控制方法 - Google Patents

一种凝胶注模陶瓷生坯热脱脂速率控制方法 Download PDF

Info

Publication number
CN112898032A
CN112898032A CN202110155639.7A CN202110155639A CN112898032A CN 112898032 A CN112898032 A CN 112898032A CN 202110155639 A CN202110155639 A CN 202110155639A CN 112898032 A CN112898032 A CN 112898032A
Authority
CN
China
Prior art keywords
thermal
alpha
green body
ceramic green
pyrolysis
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.)
Pending
Application number
CN202110155639.7A
Other languages
English (en)
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.)
Jiangxi University of Science and Technology
Original Assignee
Jiangxi University of Science and Technology
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 Jiangxi University of Science and Technology filed Critical Jiangxi University of Science and Technology
Priority to CN202110155639.7A priority Critical patent/CN112898032A/zh
Publication of CN112898032A publication Critical patent/CN112898032A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C04CEMENTS; CONCRETE; ARTIFICIAL STONE; CERAMICS; REFRACTORIES
    • C04BLIME, MAGNESIA; SLAG; CEMENTS; COMPOSITIONS THEREOF, e.g. MORTARS, CONCRETE OR LIKE BUILDING MATERIALS; ARTIFICIAL STONE; CERAMICS; REFRACTORIES; TREATMENT OF NATURAL STONE
    • C04B35/00Shaped ceramic products characterised by their composition; Ceramics compositions; Processing powders of inorganic compounds preparatory to the manufacturing of ceramic products
    • C04B35/622Forming processes; Processing powders of inorganic compounds preparatory to the manufacturing of ceramic products
    • CCHEMISTRY; METALLURGY
    • C04CEMENTS; CONCRETE; ARTIFICIAL STONE; CERAMICS; REFRACTORIES
    • C04BLIME, MAGNESIA; SLAG; CEMENTS; COMPOSITIONS THEREOF, e.g. MORTARS, CONCRETE OR LIKE BUILDING MATERIALS; ARTIFICIAL STONE; CERAMICS; REFRACTORIES; TREATMENT OF NATURAL STONE
    • C04B35/00Shaped ceramic products characterised by their composition; Ceramics compositions; Processing powders of inorganic compounds preparatory to the manufacturing of ceramic products
    • C04B35/622Forming processes; Processing powders of inorganic compounds preparatory to the manufacturing of ceramic products
    • C04B35/624Sol-gel processing
    • CCHEMISTRY; METALLURGY
    • C04CEMENTS; CONCRETE; ARTIFICIAL STONE; CERAMICS; REFRACTORIES
    • C04BLIME, MAGNESIA; SLAG; CEMENTS; COMPOSITIONS THEREOF, e.g. MORTARS, CONCRETE OR LIKE BUILDING MATERIALS; ARTIFICIAL STONE; CERAMICS; REFRACTORIES; TREATMENT OF NATURAL STONE
    • C04B2235/00Aspects relating to ceramic starting mixtures or sintered ceramic products
    • C04B2235/65Aspects relating to heat treatments of ceramic bodies such as green ceramics or pre-sintered ceramics, e.g. burning, sintering or melting processes
    • CCHEMISTRY; METALLURGY
    • C04CEMENTS; CONCRETE; ARTIFICIAL STONE; CERAMICS; REFRACTORIES
    • C04BLIME, MAGNESIA; SLAG; CEMENTS; COMPOSITIONS THEREOF, e.g. MORTARS, CONCRETE OR LIKE BUILDING MATERIALS; ARTIFICIAL STONE; CERAMICS; REFRACTORIES; TREATMENT OF NATURAL STONE
    • C04B2235/00Aspects relating to ceramic starting mixtures or sintered ceramic products
    • C04B2235/65Aspects relating to heat treatments of ceramic bodies such as green ceramics or pre-sintered ceramics, e.g. burning, sintering or melting processes
    • C04B2235/658Atmosphere during thermal treatment
    • CCHEMISTRY; METALLURGY
    • C04CEMENTS; CONCRETE; ARTIFICIAL STONE; CERAMICS; REFRACTORIES
    • C04BLIME, MAGNESIA; SLAG; CEMENTS; COMPOSITIONS THEREOF, e.g. MORTARS, CONCRETE OR LIKE BUILDING MATERIALS; ARTIFICIAL STONE; CERAMICS; REFRACTORIES; TREATMENT OF NATURAL STONE
    • C04B2235/00Aspects relating to ceramic starting mixtures or sintered ceramic products
    • C04B2235/65Aspects relating to heat treatments of ceramic bodies such as green ceramics or pre-sintered ceramics, e.g. burning, sintering or melting processes
    • C04B2235/658Atmosphere during thermal treatment
    • C04B2235/6583Oxygen containing atmosphere, e.g. with changing oxygen pressures

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Ceramic Engineering (AREA)
  • Inorganic Chemistry (AREA)
  • Materials Engineering (AREA)
  • Structural Engineering (AREA)
  • Organic Chemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

一种凝胶注模陶瓷生坯热脱脂速率控制方法,通过对凝胶注模陶瓷生坯进行非等温热失重分析,并采用热分析动力学方法获得聚合物的热解动力学参数,而后构建陶瓷生坯热脱脂过程分阶段升温控制算法,以实现对陶瓷生坯热脱脂速率的精准预测及可控。本发明适用于陶瓷生坯中的凝胶聚合物存在多个热解峰的复杂情况,对整个热脱脂过程进行分阶段保温并进一步优化,计算周期小于一分钟,适用范围广。

Description

一种凝胶注模陶瓷生坯热脱脂速率控制方法
技术领域
本发明涉及陶瓷材料胶态成型技术领域,尤其涉及一种凝胶注模陶瓷生坯热脱脂速率控制方法。
背景技术
高温结构陶瓷(如Si3N4、SiAlON等)因具有优异的综合性能,如高硬度、耐腐蚀和高化学稳定性等,使其在航空航天、军工等尖端领域备受关注。目前对于高温结构陶瓷材料而言,除改善其本质性脆的缺陷外,提高其可靠性和实现复杂形状陶瓷构件的低成本制备也是其未来发展的重要方向之一,同时这也是制约高温结构陶瓷进一步推广应用的关键技术问题。
近年来,随着科技的飞速发展,各应用领域对陶瓷构件的性能、形状复杂度、尺寸精度均提出了更高的要求。而传统的干粉压制成型技术,往往存在加工难度大、成本高,且在后续机械加工中易引入微裂纹等新缺陷的问题,从而限制了陶瓷构件的使用范围。目前常用的胶态成型技术,如注浆成型、注射成型,在制备过程中存在坯体尺寸收缩率高,且坯体在干燥、脱脂过程中极易出现开裂、变形甚至坍塌等缺陷问题,严重降低陶瓷材料的可靠性。解决该类问题的关键途径是在实现复杂形状陶瓷构件近净成型制备的同时,尽可能减少成型过程中的坯体缺陷,并降低生坯、烧结体的加工成本。
新兴的凝胶注模成型技术,因其工艺简单、坯体强度高、易于深加工、收缩形变小等一系列优势,在陶瓷材料制备领域引起了广泛关注。凝胶注模成型技术的原理是将粉体成型技术与有机化学聚合理论巧妙结合,基于有机单体和交联剂的交联反应形成具有较高强度和韧性的三维高分子聚合物网络结构,该网络结构将均匀分散在浆料中的陶瓷粉末颗粒紧密包裹并原位固定,最终获得包含陶瓷粉体颗粒和有机物且具有复合结构的特定形状陶瓷生坯,而后经干燥、脱脂、烧结获得优异综合性能的高温结构陶瓷。
但在使用凝胶注模成型技术制备陶瓷构件过程中,热脱脂作为该技术的重要步骤之一,其工艺恰当与否对成型坯体的质量至关重要。该过程对聚合物速率控制要求极为严格,若设计不合理,坯体内的热应力或残余应力以开裂、坍塌、变形或其他方式损坏陶瓷部件;除该类宏观缺陷外,热脱脂过程中引起的任何微观缺陷在随后的烧结过程中均会被进一步放大,最终影响烧结体的性能。为避免该类问题,通常使用较长的加热周期或反复尝试保温温度及保温时间以达到去除聚合物并减缓热解产物析出速率的目的,但这同时会降低生产效率,提高成本。
综上,在凝胶注模成型热脱脂过程中,完全去除聚合物而不引入诸如裂缝、翘曲之类的缺陷仍是一大挑战。因此,开发一种快速有效的从生坯中去除聚合物的热脱脂工艺,对于少乃至无缺陷结构陶瓷材料的凝胶注模制备至关重要。
发明内容
本发明要解决的技术问题是现有热脱脂工艺难以精准预测陶瓷生坯中聚合物的热解行为且不能有效控制热脱脂速率的难题,提出一种凝胶注模陶瓷生坯热脱脂速率控制方法,该方法能够在限定聚合物热解速率条件下,短时间内(<1min)完成多阶段热脱脂保温工艺的优化计算,从而快速获得低成本的优化脱脂方案,以解决上述背景技术中的问题。
本发明所解决的技术问题采用以下技术方案来实现:
一种凝胶注模陶瓷生坯热脱脂速率控制方法,具体步骤如下:
(1)采用同步热分析仪,在不同气氛下,对凝胶注模陶瓷生坯进行非等温热失重分析以获得热失重数据,并对获得的热失重数据进行多峰反卷积拟合处理,而后采用热分析动力学方法(变活化能模型、Málek法、
Figure BDA0002933319990000021
模型)获得聚合物的热解动力学参数(活化能、指前因子和机理函数),以实现对陶瓷生坯热脱脂过程中聚合物热解行为的精准预测;
(2)构建陶瓷生坯热脱脂过程分阶段升温控制算法,从而实现对陶瓷生坯热脱脂速率的精准控制。
在本发明中,所述步骤(1)中,陶瓷生坯为聚合物与陶瓷粉末所形成的复合材料,聚合物为涉及多种凝胶体系聚合形成的有机大分子,如丙烯酰胺及其衍生体系(如丙烯酰胺、基丙烯酸2-羟乙酯、甲基丙烯酰胺等单体,与之对应的交联剂为N,N′-亚甲基双丙烯酰胺、聚(乙二醇)甲基醚甲基丙烯酸酯等)和非水基凝胶体系(三羟甲基丙烷三丙烯酸酯/己二醇二丙烯酸酯)等。
在本发明中,所述步骤(1)中,陶瓷生坯需粉碎,过100~150目筛,并经充分干燥,以避免残余水分或溶剂影响聚合物的热失重数据,且每个测试条件至少重复测试3次以保证热失重数据的可靠性,整个实施过程所需陶瓷生坯样品仅需1~2g,成本较低。
在本发明中,所述步骤(1)中,热失重分析过程中,气氛包括惰性气氛(如高纯氩气、氮气)以及空气气氛。
在本发明中,所述步骤(1)的具体流程如下:
①首先连接测试仪器,确保测试仪器正常运行,而后点击同步热分析仪开始;
②由同步热分析仪对凝胶注模陶瓷生坯进行非等温热失重分析,以获得热失重数据;
③对步骤②)获得的热失重数据进行多峰反卷积拟合处理;
④对步骤③)反卷积拟合处理完毕的数据,采用热分析动力学方法以获得聚合物的热解动力学参数,热分析动力学方法包括采用多阶并行反应模型(M-PRM)对整个热脱脂过程进行热解阶段分析,并基于无模型方法(Friedman微分法)计算活化能E(α)、基于Málek法和可描述单步组合动力学反应机制的
Figure BDA0002933319990000031
模型确定机理函数f(α),由活化能E(α)结合变活化能模型计算变活化能模型系数,而后通过机理函数f(α)与变活化能模型系数计算指前因子k(α),实现对聚合物热解动力学行为预测;
⑤对步骤④)获得的聚合物热解动力学行为预测进行动力学参数验证,输出预测结果;
⑥结束。
上述步骤②)中,热失重分析实施过程至少采用三种不同的升温速率(如2.5℃/min、5℃/min、10℃/min为一组或5℃/min、10℃/min、15℃/min为一组等),以获得不同升温速率下凝胶注模陶瓷生坯的热重差分DTG曲线和反应速率dα/dT曲线。
上述步骤③)中,对获得的三组及以上的反应速率dα/dT曲线采用高斯函数进行多峰反卷积拟合,以获得多个高斯分布函数参数;多峰反卷积拟合采用Levenberg-Marquardt算法,以高斯分布函数公式(1)为拟合峰值函数,对dα/dT曲线中的多个重叠峰进行多峰拟合,最终获得匹配一个或多个热解峰的多组高斯分布函数参数:
Figure BDA0002933319990000041
式(1)中,A为峰面积,π为圆周率常数,w为半峰宽,y0和xc为实数常数。
上述步骤④)中,多阶并行反应模型(M-PRM)是一种有效的多峰拟合模型,其视整个热解过程为几个独立反应的加权和,并分别对每个子阶段进行动力学分析,确定各子阶段相应的动力学参数,假设每个子峰表示该过程中的单个假设物质的独立反应,加权因子是每个阶段的质量损失与总质量损失的比率,则热解子阶段i的动力学方程的微分形式的数学表达式为:
Figure BDA0002933319990000042
式(2)中,αi为生坯中聚合物i组分的转化率,T为温度,t为时间,k(αi)为反应子阶段i的指前因子,E(αi)为反应子阶段i的活化能,R为理想气体常数,β为升温速率,f(αi)为反应子阶段i的反应机理函数;
反应子阶段i的加权因子ri,采用以下公式计算:
Figure BDA0002933319990000043
则,整个阶段的转化率α和活化能E描述为:
Figure BDA0002933319990000044
Figure BDA0002933319990000045
式(3)和(4)中,N为反应子阶段数,mi0、mif分别为反应子阶段i中所对应反应物的初始和最终质量,m0、mf分别为聚合物的初始和最终质量,mT为聚合物在温度值为T时的质量。
上述步骤④)中,无模型方法(Friedman差分法)是一种基于速率方程的对数形式的等转换率方法,用于求解各热解子阶段的动力学参数E(αi),所涉及公式为:
Figure BDA0002933319990000051
具体过程为:基于步骤③)计算的各子阶段高斯分布函数,在等转换率条件下,绘制ln(βdαi/dT)对1/T的回归直线,通过直线的斜率从而获得对应αi的活化能E(αi)。
上述步骤④)中,变活化能模型用于描述活化能与转化率间的依赖性,并采用公式(7)所示的多项式,回归拟合获得模型参数(pi,5、pi,6、pi,7和pi,8):
E(αi)=pi,5+pi,6αi+pi,7αi 2+pi,8αi 3 (7)
式(7)中,pi,5、pi,6、pi,7和pi,8为多项式回归拟合获得模型参数。
上述步骤④)中,Málek法是确定反应机理函数f(αi)和G(αi)的一种有效方法,Málek方法认为f(αi)函数与y(αi)函数成一定比例,通过绘制实验曲线和理论曲线,若实验数据与理论曲线重合或者大部分值都落在该曲线上,则可认为其对应理论曲线的f(αi)和G(αi)为最概然机理函数,具体过程为:将αi,T,(dαi/dt)(i=1,2,…N)和αi=0.5,T0.5,(dαi/dt)0.5代入公式(8),通过(T/T0.5)2(dαi/dt)/(dαi/dt)0.5对αi作图,可绘制实验曲线:
Figure BDA0002933319990000052
式(8)中,y(αi)为定义函数,T0.5和(dαi/dt)0.5分别为转化率为0.5时的温度和反应速率。
上述步骤④)中,
Figure BDA0002933319990000053
模型是描述固态反应机理动力学的三参数转换函数,其数学表达式为:
Figure BDA0002933319990000054
式(9)中,ni、mi和pi为反应级数、幂律和扩散机制指数因子(通常为非整数);
具体过程为:将f(αi)、G(αi)、f(0.5)和G(0.5)代入公式(10),通过f(αi)·G(αi)/f(0.5)·G(0.5)对αi作图,可绘制理论曲线:
Figure BDA0002933319990000061
式(10)中,G(αi)为f(αi)机理函数的积分形式,随后基于广义简约梯度算法计算获得能满足实验数据曲线与理论曲线重合或者大部分值都落在该曲线上的
Figure BDA0002933319990000062
模型参数(mi、ni和pi)。
上述步骤④)中,指前因子在已经计算获得的机理函数f(αi)和活化能E(αi)后,以公式(2)为评判函数,基于广义简约梯度算法,以公式(11)为目标函数,计算获得指前因子k(αi):
k(αi)=exp(pi,1+pi,2αi+pi,3αi 2+pi,4αi 3) (11)
上述步骤⑤)中,对指前因子k(αi)进行动力学参数验证,具体实施流程为:依据除动力学计算用的升温制度外的加热条件获得的热失重实验数据,能否与动力学理论计算结果重合进行判断;若能基本重合,则所计算的动力学参数可行。
在本发明中,所述步骤(2)中,陶瓷生坯热脱脂过程分阶段升温控制算法通过步骤(1)获得的热解动力学参数,根据聚合物存在多热解峰的特性,采用程序设计语言构建限定聚合物反应速率的多段升温程序,以实现对聚合物热解行为的精准预测及可控;所述凝胶热脱脂过程分阶段升温控制算法的基本原理为:以最短脱脂时间为目标,通过构建限定反应速率(dα/dt)的分段升温程序实现,但在热脱脂线性升温过程中,当所处的温度对应的反应速率(dα/dt)值高于约定允许的dα/dt上限值(dα/dt)u,则线性升温过程进入保温步骤;在保温步骤中,温度保持不变,所使用的最小保温时间值为td,最大保温时间值为tu,即td≤tn≤tu,当保温时间tn所对应的反应速率(dα/dt)值低于约定允许的dα/dt下限值(dα/dt)d时,退出保温步转为线性升温步;若后续升温过程(dα/dt)值出现高于设定的反应速率(dα/dt)上限值(dα/dt)u的情况,则重复保温步骤;如此反复直到整个热脱脂过程始终维持反应速率(dα/dt)低于约定上限值(dα/dt)u;设置保温时间上限tu和保温时间下限值td,用于后续条件判断选择;同时不断调整约定允许的下限值(dα/dt)d重复上述升温-保温-升温计算步骤获得最短脱脂时间,获得不同保温段数的升温曲线控制策略。
所述步骤(2)中,多段升温程序通过给定起始dα/dt上限值(dα/dt)u和保温过程的终止的dα/dt下限值(dα/dt)d,并约束保温的最短时间td、最长时间tu以减小保温步数的方式来实现。
有益效果:本发明适用于采用多种凝胶体系制备的陶瓷生坯,可根据凝胶体系所形成聚合存在多热解峰的特性,采用上述热分析动力学方法,仅通过对非等温热失重数据计算,即可实现对其多热解阶段的有效描述和精准预测;同时,可通过限定聚合物的热解反应速率、约束保温时间的方式,在保证消除聚合物热解速率峰的前提下,实施多阶段保温策略,从而优化热脱脂工艺;更进一步,整个算法实施过程中,计算所消耗的时间小于1分钟;且由于对聚合物热解速率峰的消除,热脱脂过程中陶瓷生坯中热解气体产物的释放速率得到有效控制,能够降低坯体内应力的产生,从而减少缺陷的形成,有效提高凝胶注模陶瓷材料的可靠性;适用范围广,可推广至其他凝胶注模材料的热脱脂工艺方案优化设计;更进一步,适用于新型智能热脱脂炉的设计与研发,具有重要的理论意义和工业应用价值。
附图说明
图1为本发明中凝胶聚合物热解行为预测流程示意图。
图2为本发明中分段升温约束策略示意图。
图3为本发明中多段升温程序计算流程示意图。
图4为本发明中不同升温速率下氩气气氛热脱脂过程生坯的TG和DTG曲线示意图。
图5为本发明中处于氩气气氛下陶瓷生坯热脱脂过程聚合物热解行为预测结果示意图。
图6为本发明中不同升温速率下处于空气气氛热脱脂过程生坯的TG和DTG曲线示意图。
图7为本发明中处于空气气氛下陶瓷生坯热脱脂过程聚合物热解行为预测结果示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,以Si3N4、Al2O3和AlN粉末为原料,采用低毒水基N,N-二甲基丙烯酰胺/N,N′-亚甲基双丙烯酰胺(DMAA/MBAM)凝胶体系,凝胶注模制备SiAlON陶瓷生坯的热脱脂过程为例,进一步阐述本发明。
实施例1
采用同步热分析仪,在惰性气体气氛(高纯氩气)下,对陶瓷生坯进行非等温热失重分析,获得升温速率为2.5℃/min、5℃/min、15℃/min和20℃/min的热失重数据,如图4所示;
如图1所示的流程,根据同步热分析仪获取的热失重数据,将整个热脱脂过程分为三个热解阶段,求得热脱脂过程中SiAlON陶瓷生坯DMAA/MBAM聚合物热解动力学参数和机理函数,计算获得三个子阶段1、2和3的活化能函数,活化能随转化率的变化关系式分别为E(α)=139.862-110.481α+156.161α2-88.714α3kJ/mol、E(α)=160.791+152.496α-236.906α2+163.724α3kJ/mol和E(α)=72.132+452.830α-669.039α2+507.015α3kJ/mol;
三个子阶段1、2和3的动力学机理函数分别f(α)=(1-α)0.668α3.049(-ln(1-α))-3.874、f(α)=(1-α)0.700α3.177(-ln(1-α))-3.962和f(α)=(1-α)1.049α-0.161(-ln(1-α))0.518
依据上述获得的动力学参数,对SiAlON陶瓷生坯在惰性气氛热脱脂过程中DMAA/MBAM聚合物热解行为进行预测,如图5所示,预测结果与实验数据吻合度高,最大误差小于9%;
如图2、图3所示,计算热脱脂过程所需的保温阶段和保温时间,以控制DMAA/MBAM聚合物的最大热解速率为线性升温(1℃/min,起始温度25℃,终止温度600℃)下的30%为目标,即设置dα/dt上限0.3*(dα/dt)max和dα/dt下限为0.06*(dα/dt)max,设置保温时间上限为480min、保温时间下限值30min,通过陶瓷生坯热脱脂过程分阶段升温控制算法,计算获得惰性气氛优化脱脂工艺方案,计算结果如表1所示:
表1
Figure BDA0002933319990000091
实施例2
采用同步热分析仪,在空气气氛下,对陶瓷生坯进行非等温热失重分析,获得升温速率为5℃/min、8℃/min、10℃/min和15℃/min的热失重数据,如图6所示;
根据图1所示的流程,根据同步热分析仪获取的热失重数据,将整个热脱脂过程分为五个热解阶段,求得热脱脂过程中SiAlON陶瓷生坯中DMAA/MBAM聚合物热解动力学参数和机理函数,计算获得五个子阶段1、2、3、4和5的活化能函数,活化能随转换率的变化关系式分别为E(α)=115.027+11.950α-39.241α2+10.624α3kJ/mol、E(α)=139.595-66.162α+75.702α2-38.041α3kJ/mol、E(α)=190.854+135.755α-214.801α2+116.093α3kJ/mol、E(α)=64.068+280.086α-380.270α2+264.724α3kJ/mol和E(α)=188.257-77.086α+74.129α2-48.669α3kJ/mol,子阶段1、2、3和5的动力学机理函数均为f(α)=(1-α)0.711α2.970(-ln(1-α))-3.761,子阶段4的机理函数为f(α)=(1-α)1.050α-0.025(-ln(1-α))0.134
依据上述获得的动力学参数,对SiAlON陶瓷生坯在空气气氛热脱脂过程中DMAA/MBAM聚合物热解行为进行预测,如图7所示,预测结果与实验数据吻合度高,最大误差小于6%;
如图2、图3所示,计算热脱脂过程所需的保温阶段和保温时间,以控制聚合物的最大热解速率(dα/dt)max为线性升温(1℃/min,起始温度25℃,终止温度600℃)下的30%为目标,即设置dα/dt上限0.3*(dα/dt)max和dα/dt下限为0.06*(dα/dt)max,设置保温时间上限为480min和保温时间下限值30min,通过陶瓷生坯热脱脂过程分阶段升温控制算法,计算获得空气气氛优化脱脂工艺方案,计算结果如表2所示:
表2
Figure BDA0002933319990000101
在上述实施例1与实施例2中,对SiAlON陶瓷生坯在不同气氛热脱脂过程中DMAA/MBAM聚合物热解行为进行预测的具体流程如下:
①首先连接测试仪器,确保测试仪器正常运行,而后点击同步热分析仪开始(101);
②由同步热分析仪对凝胶注模陶瓷生坯进行非等温热失重分析(102),以获得热失重数据;
③对步骤102)获得的热失重数据进行多峰反卷积拟合处理(103);
④对步骤103)反卷积拟合处理完毕的数据,采用热分析动力学方法(104)以获得聚合物的热解动力学参数,热分析动力学方法(104)包括采用多阶并行反应模型(M-PRM)对整个热脱脂过程进行热解阶段分析,并基于无模型方法(Friedman微分法)计算活化能E(αi)(105)、基于Málek法和可描述单步组合动力学反应机制的
Figure BDA0002933319990000102
模型确定机理函数f(αi)(107),由活化能E(αi)结合变活化能模型计算变活化能模型系数(106),而后通过机理函数f(αi)(107)与变活化能模型系数(106)计算指前因子k(αi)(108),实现对聚合物热解动力学行为预测;
⑤对步骤④)获得的聚合物热解动力学行为预测进行动力学参数验证,输出预测结果(109);
⑥结束(110)。
上述步骤②)中,热失重分析实施过程至少采用三种不同的升温速率(如2.5℃/min、5℃/min、10℃/min为一组或5℃/min、10℃/min、15℃/min为一组等),以获得不同升温速率下凝胶注模陶瓷生坯的热重差分DTG曲线和反应速率dα/dT曲线。
上述步骤③)中,对获得的三组及以上的反应速率dα/dT曲线采用高斯函数进行多峰反卷积拟合,以获得多个高斯分布函数参数;多峰反卷积拟合采用Levenberg-Marquardt算法,以高斯分布函数公式(1)为拟合峰值函数,对dα/dT曲线中的多个重叠峰进行多峰拟合,最终获得匹配一个或多个热解峰的多组高斯分布函数参数:
Figure BDA0002933319990000111
式(1)中,A为峰面积,π为圆周率常数,w为半峰宽,y0和xc为实数常数。
上述步骤④)中,多阶并行反应模型(M-PRM)是一种有效的多峰拟合模型,其视整个热解过程为几个独立反应的加权和,并分别对每个子阶段进行动力学分析,确定各子阶段相应的动力学参数,假设每个子峰表示该过程中的单个假设物质的独立反应,加权因子是每个阶段的质量损失与总质量损失的比率,则热解子阶段i的动力学方程的微分形式的数学表达式为:
Figure BDA0002933319990000112
式(2)中,αi为生坯中聚合物i组分的转化率,T为温度,k(αi)为反应子阶段i的指前因子,E(αi)为反应子阶段i的活化能,R为理想气体常数,β为升温速率,f(αi)为反应子阶段i的反应机理函数;
反应子阶段i的加权因子ri,采用以下公式计算:
Figure BDA0002933319990000113
则,整个阶段的转化率α和活化能E描述为:
Figure BDA0002933319990000114
Figure BDA0002933319990000115
式(3)和(4)中,N为反应子阶段数,mi0、mif分别为反应子阶段i中所对应反应物的初始和最终质量,m0、mf分别为聚合物的初始和最终质量,mT为聚合物在温度值为T时的质量。
上述步骤④)中,无模型方法(Friedman差分法)是一种基于速率方程的对数形式的等转换率方法,用于求解各热解阶段的动力学参数E(αi),所涉及公式为:
Figure BDA0002933319990000121
具体过程为:基于步骤③)计算的各子阶段高斯分布函数,在等转换率条件下,绘制ln(βdαi/dT)对1/T的回归直线,通过直线的斜率从而获得对应αi的活化能E(αi)。
上述步骤④)中,变活化能模型用于描述活化能与转化率间的依赖性,并采用公式(7)所示的多项式,回归拟合获得模型参数(pi,5、pi,6、pi,7和pi,8):
E(αi)=pi,5+pi,6αi+pi,7αi 2+pi,8αi 3 (7)
式(7)中,pi,5、pi,6、pi,7和pi,8为多项式回归拟合获得模型参数。
上述步骤④)中,Málek法是确定反应机理函数f(αi)和G(αi)的一种有效方法,Málek方法认为f(αi)函数与y(αi)函数成一定比例,通过绘制实验曲线和理论曲线,若实验数据与理论曲线重合或者大部分值都落在该曲线上,则可认为其对应理论曲线的f(αi)和G(αi)为最概然机理函数,具体过程为:将αi,T,(dαi/dt)(i=1,2,…N)和αi=0.5,T0.5,(dαi/dt)0.5代入公式(8),通过(T/T0.5)2(dαi/dt)/(dαi/dt)0.5对αi作图,可绘制实验曲线:
Figure BDA0002933319990000122
式(8)中,y(αi)为定义函数,T0.5和(dαi/dt)0.5分别为转化率为0.5时的温度和反应速率。
上述步骤④)中,
Figure BDA0002933319990000123
模型是描述固态反应机理动力学的三参数转换函数,其数学表达式为:
Figure BDA0002933319990000124
式(9)中,ni、mi和pi为反应级数、幂律和扩散机制指数因子(通常为非整数);
具体过程为:将f(αi)、G(αi)、f(0.5)和G(0.5)代入公式(10),通过f(αi)·G(αi)/f(0.5)·G(0.5)对αi作图,可绘制理论曲线:
Figure BDA0002933319990000131
式(10)中,G(αi)为f(αi)机理函数的积分形式,随后基于广义简约梯度算法计算获得能满足实验数据曲线与理论曲线重合或者大部分值都落在该曲线上的
Figure BDA0002933319990000132
模型参数(mi、ni和pi)。
上述步骤④)中,指前因子在已经计算获得的机理函数f(αi)和活化能E(αi)后,以公式(2)为评判函数,基于广义简约梯度算法,以公式(11)为目标函数,计算获得指前因子k(αi):
k(αi)=exp(pi,1+pi,2αi+pi,3αi 2+pi,4αi 3) (11)
上述步骤⑤)中,对指前因子k(αi)进行动力学参数验证,具体实施流程为:依据除动力学计算用的升温制度外的加热条件获得的热失重实验数据,能否与动力学理论计算结果重合进行判断;若能基本重合,则所计算的动力学参数可行。
陶瓷生坯热脱脂过程分阶段升温控制算法通过获得的热解动力学参数,根据聚合物存在多热解峰的特性,采用程序设计语言开发构建限定聚合物反应速率的多段升温程序,以实现对聚合物热解行为的精准预测及可控;所述凝胶热脱脂过程分阶段升温控制算法的基本原理为:以最短脱脂时间为目标,通过构建限定反应速率(dα/dt)的分段升温程序实现,如图2所示,在热脱脂线性升温过程中,当所处的温度对应的反应速率(dα/dt)值高于约定允许的dα/dt上限值(dα/dt)u,则线性升温过程进入保温步骤;在保温步骤中,温度保持不变,所使用的最小保温时间值为td,最大保温时间值为tu,即td≤tn≤tu,当保温时间tn所对应的反应速率(dα/dt)值低于约定允许的dα/dt下限值(dα/dt)d时,退出保温步转为线性升温步;若后续升温过程(dα/dt)值出现高于设定的反应速率(dα/dt)上限值(dα/dt)u的情况,则重复保温步骤;如此反复直到整个热脱脂过程始终维持反应速率(dα/dt)低于约定上限值(dα/dt)u;设置保温时间上限和保温时间下限值,用于后续条件判断选择;同时不断调整约定允许的下限值(dα/dt)d重复上述升温-保温-升温计算步骤获得最短脱脂时间,获得不同保温段数的升温曲线控制策略。
所述多段升温程序通过给定起始dα/dt上限值(dα/dt)u和保温过程的终止的dα/dt下限值(dα/dt)d,并约束保温的最短时间td、最长时间tu以减小保温步数的方式来实现,具体计算流程如图3所示:
步骤301),初始状态时,时间t为0min,温度T为25℃;设置dα/dt上限值(dα/dt)u、dα/dt下限值(dα/dt)d、保温最短时间td、保温最长时间tu、时间步长dt和温度步长dT,当前升温制度保温步为0步;
步骤302),输入的动力学参数为多阶并行反应模型(M-PRM)中涉及各子阶段加权因子ri,各子阶段包含的机理函数f(αi)中的ni、mi、pi、活化能E(αi)的pi,5,pi,6,pi,7,pi,8和指前因子k(αi)中的pi,1,pi,2,pi,3,pi,4
步骤303),更新时间、温度是每循环迭代一次,时间和温度值增加一个时间步长dt和温度步长dT;
步骤304),动力学计算,采用的是多阶并行反应模型(M-PRM),通过公式(12)计算各子阶段的热解速率dαi/dt:
Figure BDA0002933319990000141
式(12)中,ln[k(αi)]=pi,1+pi,2αi+pi,3αi 2+pi,4αi 3,E(αi)=pi,5+pi,6αi+pi,7αi 2+pi,8αi 3
反应速率(dα/dt)和转化率α采用公式(13)和公式(14)计算,其中公式(14)中积分函数采用分段变步长梯形积分计算:
Figure BDA0002933319990000142
Figure BDA0002933319990000143
式(13)、(14)中,nc为拟物质总数量,nt为温度分段数,ci为物质加权因子,Δtj为第j温度段的时间;
步骤305),判断时间值,在判断当前计算过程中,所处的时间值是否达到预先设定的终止值,满足则进入步骤306)执行,否则进入步骤312)执行;
步骤306),用于通过文本输出、显示屏输出等方式打印出当前计算所用的保温制度,而后进入步骤307);
步骤307),用于更新保温段数、最短保温时间。通过对比程序中当前已存储的最短保温时间及最短保温段数数据,若比当前值更小,则替换当前值为本次计算所用的保温制度;
步骤308),用于判断当前计算步动态更新的反应速率下限值(dα/dt)d是否大于等于dα/dt上限值(dα/dt)u,若满足,则进入步骤309)执行,否则进入步骤311)执行;
步骤309),用于通过文本输出、显示屏输出等方式打印出所有计算步中不同保温段数条件的最短保温时间所用的保温制度,而后进入步骤310)结束,其输出数据可用于热脱脂炉的升温制度编制;
步骤310)结束;
步骤311),更新反应速率dα/dt下限值(dα/dt)d
步骤312),用于判断当前计算步的dα/dt值是否高于上限值(dα/dt)u,若满足,则进入步骤313),否则返回步骤303);
步骤313),为当前保温步数加1,而后进入步骤314);
步骤314),用于更新时间tn值,加上dt时间微元,温度值不变,而后进入步骤315);
步骤315),动力学计算,采用多阶并行反应模型(M-PRM)计算,随后进入步骤316);
步骤316),判断当前dα/dt值是否小于反应速率下限值(dα/dt)d,若满足则进入步骤317),否则进入步骤319);
步骤317),判断当前tn值是否大于等于保温时间下限td,若满足则进入步骤318),否则进入步骤314);
步骤318),记录当前计算步的保温温度、时间,用于后续结果输出;而后返回至步骤303);
步骤319),判断当前保温时间是否高于保温时间上限tu,若满足则进入步骤318),否则进入步骤314)。

Claims (9)

1.一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,具体步骤如下:
(1)采用同步热分析仪,在不同气氛下,对凝胶注模陶瓷生坯进行非等温热失重分析,以获得热失重数据,并对获得的热失重数据进行多峰反卷积拟合处理,而后采用热分析动力学方法获得聚合物的热解动力学参数,以实现对陶瓷生坯热脱脂过程中聚合物热解行为的精准预测;
(2)构建陶瓷生坯热脱脂过程分阶段升温控制算法,从而实现对陶瓷生坯热脱脂速率的精准控制。
2.根据权利要求1所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,所述陶瓷生坯为聚合物与陶瓷粉末所形成的复合材料。
3.根据权利要求2所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,所述聚合物含有的热解峰不少于2个。
4.根据权利要求1所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,所述非等温热失重分析中采用至少三种不同的升温速率,以获得不同升温速率下凝胶注模陶瓷生坯的热失重数据。
5.根据权利要求1所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,所述热分析动力学方法包括采用多阶并行反应模型对整个热脱脂过程进行热解阶段分析,并基于无模型方法计算活化能E(α)、基于Málek法和可描述单步组合动力学反应机制的
Figure FDA0002933319980000011
模型确定各阶段的机理函数f(α),由活化能E(α)结合变活化能模型计算变活化能模型系数,而后通过机理函数f(α)与变活化能模型系数计算指前因子k(α),以获得各阶段的反应动力学参数,实现对聚合物热解动力学行为预测。
6.根据权利要求1所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,所述陶瓷生坯热脱脂过程分阶段升温控制算法通过步骤(1)获得的热解动力学参数,根据聚合物存在多热解峰的特性,采用程序设计语言开发构建限定聚合物反应速率的多段升温程序,以实现对聚合物热解行为的精准预测及可控。
7.根据权利要求6所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,所述陶瓷生坯热脱脂过程分阶段升温控制算法以最短脱脂时间为目标,通过构建限定反应速率的分段升温程序实现。
8.根据权利要求7所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法,其特征在于,所述分段升温程序实现具体表述为:通过给定dα/dt的起始约束上限值和保温过程的终止约束下限值,并约束保温的最短时间与最长时间以减小保温步数的方式实现。
9.根据权利要求1所述的一种凝胶注模陶瓷生坯热脱脂速率控制方法的制备方法,其特征在于,所述气氛包括惰性气氛或空气气氛。
CN202110155639.7A 2021-02-04 2021-02-04 一种凝胶注模陶瓷生坯热脱脂速率控制方法 Pending CN112898032A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110155639.7A CN112898032A (zh) 2021-02-04 2021-02-04 一种凝胶注模陶瓷生坯热脱脂速率控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110155639.7A CN112898032A (zh) 2021-02-04 2021-02-04 一种凝胶注模陶瓷生坯热脱脂速率控制方法

Publications (1)

Publication Number Publication Date
CN112898032A true CN112898032A (zh) 2021-06-04

Family

ID=76122361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110155639.7A Pending CN112898032A (zh) 2021-02-04 2021-02-04 一种凝胶注模陶瓷生坯热脱脂速率控制方法

Country Status (1)

Country Link
CN (1) CN112898032A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116135818A (zh) * 2023-03-03 2023-05-19 上海科技大学 一种陶瓷坯体的排胶方法和陶瓷工件的制备方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002201075A (ja) * 2000-10-27 2002-07-16 Toshiba Corp 窒化けい素セラミックス基板およびそれを用いた窒化けい素セラミックス回路基板並びにその製造方法
WO2011009052A2 (en) * 2009-07-17 2011-01-20 University Of Florida Research Foundation, Inc. Electroactivated peptides and biosensors
CN108467272A (zh) * 2018-03-21 2018-08-31 浙江立泰复合材料股份有限公司 一种碳化硼/碳化硅层状复合陶瓷材料的制备方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002201075A (ja) * 2000-10-27 2002-07-16 Toshiba Corp 窒化けい素セラミックス基板およびそれを用いた窒化けい素セラミックス回路基板並びにその製造方法
WO2011009052A2 (en) * 2009-07-17 2011-01-20 University Of Florida Research Foundation, Inc. Electroactivated peptides and biosensors
CN108467272A (zh) * 2018-03-21 2018-08-31 浙江立泰复合材料股份有限公司 一种碳化硼/碳化硅层状复合陶瓷材料的制备方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
I.UL MOHSIN ET AL.: "Simulation and optimisation for thermal debinding of copper MIM parts using thermokinetic analysis", 《POWDER METALLURGY》 *
JING LI ET AL.: "Multistage kinetic analysis of DMAA/MBAM polymer removal from gelcast ceramic parts using a multistage parallel reaction model and model-free method", 《RSC ADVANCES》 *
JING LI ET AL.: "Thermal debinding behavior of a low-toxic DMAA polymer for gelcast ceramic parts based on TG-FTIR and kinetic modeling", 《RSC ADVANCES》 *
MEHDI SALEHI ET AL.: "Kinetic analysis of the polymer burnout in ceramic thermoplastic processing of the YSZ thin electrolyte structures using model free method", 《APPLIED ENERGY》 *
X.F.YANG ET AL.: "Dynamics of water debinding in ceramic injection moulding", 《ADVANCES IN APPLIED CERAMICS》 *
山香教师资格考试命题研究中心: "《信息技术学科知识与教学能力》", 31 March 2020, 首都师范大学出版社 *
高桥诚 著: "《创造技法实用手册》", 31 May 1992, 湖南科学技术出版社 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116135818A (zh) * 2023-03-03 2023-05-19 上海科技大学 一种陶瓷坯体的排胶方法和陶瓷工件的制备方法

Similar Documents

Publication Publication Date Title
Li et al. Modeling hot deformation behavior of low-cost Ti-2Al-9.2 Mo-2Fe beta titanium alloy using a deep neural network
Bae et al. Influence of residual monomer on cracking in ceramics fabricated by stereolithography
CN106119771B (zh) 一种料浆法制备铝硅共渗涂层的渗层厚度控制方法
CN112898032A (zh) 一种凝胶注模陶瓷生坯热脱脂速率控制方法
CN103886125B (zh) 一种钛合金热复合成形数值模拟方法
CN112632813B (zh) 一种大厚度树脂基复合材料固化制度的优化方法
CN106503389A (zh) 一种快速确定热冲压生产线工艺参数的数值方法
CN115015318B (zh) 一种大型构件热锻全流程宏微观分析方法及平台
CN116312902B (zh) 一种确定材料Johnson-Cook动态本构模型的方法
CN115659688B (zh) 提升增材铝合金及其复合材料可靠性的工艺确定方法
CN116959598A (zh) 受饱和晶粒尺寸控制的合金晶粒生长预测方法
CN114349490A (zh) 一种二氧化硅气凝胶隔热材料及其制备方法
CN113534661A (zh) 基于卡尔曼滤波和非最小状态空间的电阻炉温度控制方法
CN113343515A (zh) 一种小口径玻璃镜片的多工位精密热压成型模具轮廓偏移补偿方法
CN114927175B (zh) 一种用于材料拉伸本构参数拟合的高通量计算方法
Gao et al. Size effect in curing optimization for thick composite laminates
CN111506992A (zh) 一种获得玻璃材料结构松弛参数的方法
CN108694288B (zh) 快速获取步进梁式钢坯加热炉不同产率下设置温度的方法
CN106709129A (zh) 一种多孔材料隔热性能的优化控制方法
Wen et al. Check for updates Prediction and Control of Microstructure Evolution of a Novel P/M Nickel-Based Superalloy During Near-Isothermal Forging
CN110608660B (zh) 一种电涡流传感器位移、温度预测方法
Zhang et al. A rational square-root B-spline model approximation and control of output probability density functions for dynamic stochastic systems
CN117216891B (zh) 一种管式换热器设计方法
LIU et al. Numerical Simulation to Optimize Heating System in Vacuum Carburizing Furnace
CN118586186A (zh) 一种gh4706合金的高精度本构模型构建方法

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210604

RJ01 Rejection of invention patent application after publication