一种硅碳基陶瓷涂层界面力学特性的跨尺度分析方法技术

本发明专利技术公开了一种硅碳基陶瓷涂层界面力学特性的跨尺度分析方法，采用基于混合法则的Tersoff三体势函数描述碳化硅基陶瓷材料内部各原子的相互作用；采用EAM势函数描述金属基体材料内部原子的相互作用；采用Buckingham函数和Morse势函数描述碳化硅基陶瓷与金属基体界面的相互作用；利用大规模MD对材料界面进行I型和II型断裂模拟，利用大规模MD断裂模拟中初始裂纹附近的应力和位移参数，提取CZM模型的TSL信息，获得硅碳基陶瓷和金属基体各自的界面TSL函数。本发明专利技术能更有效地对不同材料界面进行分析及设计，因而提高了涂层结构的设计效率并大大降低了设计成本。

【技术实现步骤摘要】

本专利技术涉及新型纳米结构的分析及设计领域，特别是一种硅碳基陶瓷涂层界面力学特性的跨尺度模拟方法。
技术介绍
随着近代超音速飞机、火箭、人造卫星以及原子能等尖端新技术的发展，相应的工作条件日益严格，对材料耐高温及超高温、耐腐蚀、抗震动、抗疲劳、抗温度急变以及耐火焰冲刷等性能要求越来越高。原来使用的合金钢、高温合金等金属因温度太高而产生过量的蠕变，或因高温气流的作用而产生过快的腐蚀与成片剥落等现象而寿命过短。为了消除或延缓这些现象的发生，在金属或其他高温材料表面上施加各种涂层的方法愈来愈被人们所重视。涂层技术在涂层材料、制备方法、性能表征、技术手段和防护机理等方面都取得了显著成就。纳米涂层以其优异的机械性能，良好的性价比优势在材料表面保护方面显示出广阔的应用前景。利用涂层技术改善材料表面的性能，对其起到防护、密封、抗磨、抗冲击、减振、隔热等作用，可提升材料的可靠性，延长使用寿命。其中，将耐高温陶瓷涂层应用于低热耗发动机是一个前景广阔的领域。使用陶瓷热障涂层制成的发动机涡轮叶片能在1600℃的高温下运行,这样可以提高发动机的热效率到60％以上。因此，耐高温陶瓷涂层的发展将对汽车和航天等工业领域起到巨大的推动作用。然而，由于陶瓷涂层与金属基体的热膨胀系数相差较大，在恶劣的使用工作环境条件下，冷热温度急变、震动、气流冲刷等因素可能造成涂层脱落，使很多耐高温陶瓷涂层无法完全获得实际应用。涂层中各种界面结构的力学性能对其可靠性的影响至关重要。在传统的涂层设计中，涂层的材料选择和结构设计主要采用试验分析和经验设计的方法。由于各种测量方法存在周期长、费用高等局限性，要想通过大量的测试研究来寻求最佳设计是非常困难的。分子动力学(MD)是纳米尺度应用最为普遍的建模和模拟方法。然而，对于MD模型而言，其计算时间随着原子数的增加而急剧增加。包含多种界面结构的涂层体系基本处于微米以上量级，采用MD模型描述整个系统完全超出了其建模范围。内聚力模型(CZM)是宏观力学中描述界面的主要方法，CZM模型通过traction-separation laws(TSL)描述断裂过程的主要特征，而精确获得TSL是通过大量实验的不断试错得到的，效率低且会消耗大量人力物力。
技术实现思路
本专利技术的目的在于，提供一种硅碳基陶瓷涂层界面力学特性的跨尺度分析方法。本专利技术能更有效地对不同材料界面进行分析及设计，因而提高了涂层结构的设计效率并大大降低了设计成本。本专利技术的技术方案：一种硅碳基陶瓷涂层界面力学特性的跨尺度分析方法，其特征在于：采用基于混合法则的Tersoff三体势函数描述碳化硅基陶瓷材料内部各原子的相互作用；采用EAM势函数描述金属基体材料内部原子的相互作用(常规方法)；采用Buckingham函数和Morse势函数描述碳化硅基陶瓷与金属基体界面的相互作用；利用大规模MD对材料界面进行I型和II型断裂模拟，利用大规模MD断裂模拟中初始裂纹附近的应力和位移参数，提取CZM模型的TSL信息，获得硅碳基陶瓷和金属基体各自的界面TSL函数。前述的硅碳基陶瓷涂层界面力学特性的跨尺度分析方法中，所述获得界面TSL函数的具体方法按下述步骤进行：1)在硅碳基陶瓷和金属基体界面处加入一条0.5nm×0.25nm的初始裂纹；2)以盒子底部和顶部1nm范围的原子区域作为加载层，拉伸和剪切加载通过在加载层原子施加作用力实现；3)I型和II型断裂时，每次施加作用力10-10N，并弛豫10ps使系统达到平衡状态；4)将沿初始裂纹方向±8nm范围作为TSL信息的采集区域，并将采集区域分为两层，第一层沿裂纹方向分为宽度0.5nm×0.5nm的子区域，用于收集裂尖附件应力应变信息；第二层沿裂纹方向分为宽度1nm×1nm的子区域，用于收集垂直裂纹方向的应力应变信息；5)每隔0.5ps，计算I型和II型断裂每个TSL信息采集子区域内的平均拉伸应力σyy、平均剪切应力τxy，以及子区域界面两边的原子在拉伸方向和剪切方向相对位移的平方和λ，直到界面完全断裂为止；6)根据模式I和模式II局部应力-张开曲线与载荷类型的相关性，以式(6)为基本形式构建CZM模型的TSL函数T(λ)：T(λ)＝A1σyyλ+A2τxyλexp(Bλ) (6)其中T为加载层的累积作用力，A1、A2、B是系数。前述的硅碳基陶瓷涂层界面力学特性的跨尺度分析方法中，所述采用Tersoff三体势函数描述碳化硅基陶瓷材料内部各原子的相互作用，其具体形式为： E = Σ i E i = 1 2 Σ i ≠ j V i j , V i j = f C ( r i j ) [ A i j exp ( - λ i j r i j ) - b i j B i j exp ( - μ i j r i j ) ] ]]> f C ( r i j ) = 1 , r i j 本文档来自技高网...

【技术保护点】
一种硅碳基陶瓷涂层界面力学特性的跨尺度分析方法，其特征在于：采用基于混合法则的Tersoff三体势函数描述碳化硅基陶瓷材料内部各原子的相互作用；采用EAM势函数描述金属基体材料内部原子的相互作用；采用Buckingham函数和Morse势函数描述碳化硅基陶瓷与金属基体界面的相互作用；利用大规模MD对材料界面进行I型和II型断裂模拟，利用大规模MD断裂模拟中初始裂纹附近的应力和位移参数，提取CZM模型的TSL信息，获得硅碳基陶瓷和金属基体各自的界面TSL函数。

【技术特征摘要】
1.一种硅碳基陶瓷涂层界面力学特性的跨尺度分析方法，其特征在于：采用基于混合法则的Tersoff三体势函数描述碳化硅基陶瓷材料内部各原子的相互作用；采用EAM势函数描述金属基体材料内部原子的相互作用；采用Buckingham函数和Morse势函数描述碳化硅基陶瓷与金属基体界面的相互作用；利用大规模MD对材料界面进行I型和II型断裂模拟，利用大规模MD断裂模拟中初始裂纹附近的应力和位移参数，提取CZM模型的TSL信息，获得硅碳基陶瓷和金属基体各自的界面TSL函数。2.根据权利要求1所述的硅碳基陶瓷涂层界面力学特性的跨尺度分析方法，其特征在于：所述获得界面TSL函数的具体方法按下述步骤进行：1)在硅碳基陶瓷和金属基体界面处加入一条0.5nm×0.25nm的初始裂纹；2)以盒子底部和顶部1nm范围的原子区域作为加载层，拉伸和剪切加载通过在加载层原子施加作用力实现；3)I型和II型断裂时，每次施加作用力10-10N，并弛豫10ps使系统达到平衡状态；4)将沿初始裂纹方向±8nm范围作为TSL信息的采集区域，并将采集区域分为两层，第一层沿裂纹方向分为宽度0.5nm×0.5nm的子区域，用于收集裂尖附件应力应变信息；第二层沿裂纹方向分为宽度1nm×1nm的子区域，用于收集垂直裂纹方向的应力应变信息；5)每隔0.5ps，计算I型和II型断裂每个TSL信息采集子区域内的平均拉伸应力σyy、平均剪切应力τxy，以及子区域界面两边的原子在拉伸方向和剪切方向相对位移的平方和λ，直到界面完全断裂为止；6)根据模式I和模式II局部应力-张开曲线与载荷类型的相关性，以式(6)为基本形式构建CZM模型的TSL函数T(λ)：T(λ)＝A1σyyλ+A2τxyλexp(Bλ) (6)其中T为加载层的累积作用力，A1、A2、B是系数。3.根据权利要求1所述的硅碳基陶瓷涂层界面力学特性的跨尺度分析方法，其特征在于：所述采用Tersoff三体势函数描述碳化硅基陶瓷材料内部各原子的相互作用，其具体形式为： E = Σ i E i = 1 2 Σ i ≠ j V i j , V i j = f C ( r i j ) [ A i j exp ( - λ i j r i j ) - b i j B i j exp ( - μ i j r i j ) ] ]]> f C ( r i j ) = 1 , r i j < R i j 1 2 + 1 2 cos [ π ( r i j - R i ...

【专利技术属性】
技术研发人员：廖宁波，周峰，薛伟，
申请(专利权)人：温州大学，
类型：发明
国别省市：浙江;33