单向陶瓷基复合材料任意应变加卸载本构关系预测方法技术

本发明专利技术属于复合材料应力应变行为预测方法技术领域，涉及一种单向陶瓷基复合材料任意应变加卸载本构关系预测方法，目的在于提供一种能快速预测单向陶瓷基复合材料任意应变加卸载应力应变行为的预测方法。本发明专利技术通过将应变表示为正向滑移区长度或反向滑移区长度的二次多项式，推导出了正向滑移区长度及反向滑移区长度的解析式，进而计算得对应的应力。本发明专利技术无需多次反复迭代计算，即可实现单向陶瓷基复合材料任意应变加卸载应力应变行为的快速预测。

Prediction Method of Constitutive Relation of Unidirectional Ceramic Matrix Composites under Arbitrary Strain Loading and Unloading

The invention belongs to the technical field of stress-strain behavior prediction methods for composite materials, and relates to a method for predicting the constitutive relationship of unidirectional ceramic matrix composites under arbitrary strain loading and unloading. The purpose is to provide a method for predicting the stress-strain behavior of unidirectional ceramic matrix composites under arbitrary strain loading and unloading. By expressing the strain as a quadratic polynomial of the length of the forward slip zone or the reverse slip zone, the analytical expressions of the length of the forward slip zone and the reverse slip zone are derived, and the corresponding stresses are calculated. The invention can realize fast prediction of stress-strain behavior of unidirectional ceramic matrix composites under arbitrary strain loading and unloading without repeated iteration calculation.

【技术实现步骤摘要】
单向陶瓷基复合材料任意应变加卸载本构关系预测方法
本专利技术属于复合材料应力应变行为预测方法
，具体涉及单向陶瓷基复合材料任意应变加卸载本构关系预测方法。
技术介绍
陶瓷基复合材料具有比刚度大、比强度高、耐高温、耐腐蚀、密度低等诸多优点，在航空、航天、核电等行业的热端部件上具有广泛的应用前景。为了更加高效、安全地应用陶瓷基复合材料，需要对其应力应变行为进行预测。在实际应用过程中，陶瓷基复合材料不得不承受变幅循环载荷的作用，且载荷过程无明显规律。此外，在位移型有限元计算的过程中，由应变计算应力，是形成刚度矩阵所必需的。因而，需要一种陶瓷基复合材料任意应变加卸载本构关系的预测方法，本构关系即应力应变行为。目前国内外尚无陶瓷基复合材料任意应变加卸载应力应变行为预测方法的报道。宋迎东等(宋迎东，张盛，高希光，方光武，单向陶瓷基复合材料任意加卸载应力应变行为预测方法，中国专利技术专利，授权号：CN104866690B)提供了一种单向陶瓷基复合材料任意加卸载应力应变行为预测方法，但是该方法只提供了依据应力载荷过程计算应变载荷过程的解析表达式，而未提供由应变载荷过程计算应力载荷过程的解析表达式，因而无法实现依据应变载荷过程快速计算出应力载荷过程。由此可见，单向陶瓷基复合材料任意应变加卸载应力应变行为预测方法是本
一项重要而难以解决的关键技术。
技术实现思路
本专利技术针对现有技术中的不足，提供一种单向陶瓷基复合材料任意应变加卸载本构关系预测方法，适用于任意的应变加卸载过程，无需反复迭代求解，计算效率高。为实现上述目的，本专利技术采用以下技术方案：单向陶瓷基复合材料任意应变加卸载本构关系预测方法，其特征在于，包括以下步骤：步骤1：判断复合材料损伤情况；步骤2：计算基体裂纹平均间距；步骤3：假设滑移区分布情况；步骤4：将复合材料的平均应变表示为滑移区长度的函数；步骤5：整理出关于滑移区长度的一元二次方程；步骤6：通过求根公式求出步骤5中的一元二次方程的解，即滑移区长度；步骤7：通过滑移区长度计算出当前应力；步骤8：验证滑移区分布情况。为优化上述技术方案，采取的具体措施还包括：进一步地，所述步骤1中，若无损伤，则采用混合率公式计算复合材料应力σ，混合率公式为：σ＝(EfVf+EmVm)ε其中，ε为应变，E为弹性模量，V为体积分数；下标f表示纤维，下标m表示基体；若采用混合率公式计算得的应力小于初始基体开裂应力，则复合材料无损伤，其应力等于采用混合率公式计算结果；否则，复合材料发生损伤，执行步骤2。进一步地，所述步骤2中，基体裂纹平均间距L由以下公式计算得到：其中，Lsat为基体裂纹饱和间距，σ0及m为概率统计参数。进一步地，所述步骤3中，假设滑移区分布情况如下：a)当已存在n个正向滑移区，n个反向滑移区时：若当前应变大于第n-1个正向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-2个反向滑移区；若当前应变小于第n-1个正向滑移区对应的应变，而大于第n个正向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-1个反向滑移区；若当前应变小于第n个正向滑移区对应的应变，而大于第n个反向滑移区对应的应变时，滑移区分布假设为n+1个正向滑移区，n个反向滑移区；若当前应变小于第n个反向滑移区对应的应变，而大于第n-1个反向滑移区对应的应变，滑移区分布仍为n个正向滑移区，n个反向滑移区；若当前应变小于第n-1个反向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-1个反向滑移区；b)当已存在n个正向滑移区，n-1个反向滑移区时：若当前应变大于第n-1个正向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-2个反向滑移区；若当前应变小于第n-1个正向滑移区对应的应变，而大于第n个正向滑移区对应的应变，则滑移区分布仍为n个正向滑移区，n-1个反向滑移区；若当前应变小于第n个正向滑移区对应的应变，而大于第n-1个反向滑移区对应的应变时，滑移区分布假设为n个正向滑移区，n个反向滑移区；若当前应变小于第n-1个反向滑移区对应的应变时，滑移区分布假设为n-1个正向滑移区，n-1个反向滑移区。进一步地，所述步骤4中，在不考虑纤维断裂的条件下，复合材料的平均应变即等于纤维的平均应变其中，σf(x)表示纤维轴向应力分布，x表示纤维轴向坐标，下标c表示复合材料，α为热膨胀系数，ΔT为环境温度与复合材料制备温度之差；关于滑移区的分布情况，有如下两种情况：a)当假设存在n个正向滑移区，n-1个反向滑移区时，将关于未知量及与呈线性关系的σ项单独积分，复合材料应变表示为：其中，表示第n个正向滑移区长度，rf表示纤维半径，τi表示界面剪应力，表示第1个正向滑移区长度，表示第n-1个反向滑移区长度；b)当假设存在n个正向滑移区，n个反向滑移区时，将关于未知量及与呈线性关系的σ项单独积分，复合材料应变表示为：其中，表示第n个反向滑移区长度。进一步地，所述步骤5中，关于滑移区长度的一元二次方程如下：a)当假设存在n个正向滑移区，n-1个反向滑移区时，复合材料应变表达式整理为关于的一元二次方程：设则式中a、b、c分别为：b)当假设存在n个正向滑移区，n个反向滑移区时，复合材料应变表达式整理为关于的一元二次方程：设则式中a、b、c分别为：进一步地，所述步骤6中，求解如下：a)当假设存在n个正向滑移区，n-1个反向滑移区时，由一元二次方程求根公式，可得：b)当假设存在n个正向滑移区，n个反向滑移区时，由一元二次方程求根公式，可得：进一步地，所述步骤7中，通过滑移区长度计算出当前应力如下：a)当假设存在n个正向滑移区，n-1个反向滑移区时，可得：b)当假设存在n个正向滑移区，n个反向滑移区时，可得：进一步地，所述步骤8中，验证滑移区分布情况如下：a)当已存在n个正向滑移区，n-1个反向滑移区时：若当前计算得的应力大于第n-1个正向滑移区对应的应力，则滑移区分布变为n-1个正向滑移区，n-2个反向滑移区；若当前计算得的应力小于第n-1个正向滑移区对应的应力，而大于第n个正向滑移区对应的应力，则滑移区分布仍为n个正向滑移区，n-1个反向滑移区；若当前计算得的应力小于第n个正向滑移区对应的应力，而大于第n-1个反向滑移区对应的应力时，滑移区分布变为n个正向滑移区，n个反向滑移区；若当前计算得的应力小于第n-1个反向滑移区对应的应力时，滑移区分布变为n-1个正向滑移区，n-1个反向滑移区；b)当已存在n个正向滑移区，n个反向滑移区时：若当前计算得的应力大于第n-1个正向滑移区对应的应力，则滑移区分布变为n-1个正向滑移区，n-2个反向滑移区；若当前计算得的应力小于第n-1个正向滑移区对应的应力，而大于第n个正向滑移区对应的应力，则滑移区分布变为n-1个正向滑移区，n-1个反向滑移区；若当前计算得的应力小于第n个正向滑移区对应的应力，而大于第n个反向滑移区对应的应力时，滑移区分布变为n+1个正向滑移区，n个反向滑移区；若当前计算得的应力小于第n个反向滑移区对应的应力，而大于第n-1个反向滑移区对应的应力，滑移区分布仍为n个正向滑移区，n个反向滑移区；若当前计算得的应力小于第n-1个反向滑移区对应的应力，则滑移区分布变为n-1个正向滑移区，n本文档来自技高网...

【技术保护点】
1.单向陶瓷基复合材料任意应变加卸载本构关系预测方法，其特征在于，包括以下步骤：步骤1：判断复合材料损伤情况；步骤2：计算基体裂纹平均间距；步骤3：假设滑移区分布情况；步骤4：将复合材料的平均应变表示为滑移区长度的函数；步骤5：整理出关于滑移区长度的一元二次方程；步骤6：通过求根公式求出步骤5中的一元二次方程的解，即滑移区长度；步骤7：通过滑移区长度计算出当前应力；步骤8：验证滑移区分布情况。

【技术特征摘要】
1.单向陶瓷基复合材料任意应变加卸载本构关系预测方法，其特征在于，包括以下步骤：步骤1：判断复合材料损伤情况；步骤2：计算基体裂纹平均间距；步骤3：假设滑移区分布情况；步骤4：将复合材料的平均应变表示为滑移区长度的函数；步骤5：整理出关于滑移区长度的一元二次方程；步骤6：通过求根公式求出步骤5中的一元二次方程的解，即滑移区长度；步骤7：通过滑移区长度计算出当前应力；步骤8：验证滑移区分布情况。2.如权利要求1所述的单向陶瓷基复合材料任意应变加卸载本构关系预测方法，其特征在于：所述步骤1中，若无损伤，则采用混合率公式计算复合材料应力σ，混合率公式为：σ＝(EfVf+EmVm)ε其中，ε为应变，E为弹性模量，V为体积分数；下标f表示纤维，下标m表示基体；若采用混合率公式计算得的应力小于初始基体开裂应力，则复合材料无损伤，其应力等于采用混合率公式计算结果；否则，复合材料发生损伤，执行步骤2。3.如权利要求2所述的单向陶瓷基复合材料任意应变加卸载本构关系预测方法，其特征在于：所述步骤2中，基体裂纹平均间距L由以下公式计算得到：其中，Lsat为基体裂纹饱和间距，σ0及m为概率统计参数。4.如权利要求3所述的单向陶瓷基复合材料任意应变加卸载本构关系预测方法，其特征在于：所述步骤3中，假设滑移区分布情况如下：a)当已存在n个正向滑移区，n个反向滑移区时：若当前应变大于第n-1个正向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-2个反向滑移区；若当前应变小于第n-1个正向滑移区对应的应变，而大于第n个正向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-1个反向滑移区；若当前应变小于第n个正向滑移区对应的应变，而大于第n个反向滑移区对应的应变时，滑移区分布假设为n+1个正向滑移区，n个反向滑移区；若当前应变小于第n个反向滑移区对应的应变，而大于第n-1个反向滑移区对应的应变，滑移区分布仍为n个正向滑移区，n个反向滑移区；若当前应变小于第n-1个反向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-1个反向滑移区；b)当已存在n个正向滑移区，n-1个反向滑移区时：若当前应变大于第n-1个正向滑移区对应的应变，则滑移区分布假设为n-1个正向滑移区，n-2个反向滑移区；若当前应变小于第n-1个正向滑移区对应的应变，而大于第n个正向滑移区对应的应变，则滑移区分布仍为n个正向滑移区，n-1个反向滑移区；若当前应变小于第n个正向滑移区对应的应变，而大于第n-1个反向滑移区对应的应变时，滑移区分布假设为n个正向滑移区，n个反向滑移区；若当前应变小于第n-1个反向滑移区对应的应变时，滑移区分布假设为n-1个正向滑移区，n-1个反向滑移区。5.如权利要求4所述的单向陶瓷基复合材料任意应变加卸载本构关系预测方法，其特征在于：所述步骤4中，在不考虑纤维断裂的条件下，复合材料的平均应变即等于纤维的平均应变其中，σf(x)表示纤维轴向应力分布，x表示纤维轴向坐标，下标c表示复合材料，α为热膨胀系数，ΔT为环境温度与复合材料制备温度之差；关于滑移区的分布情况，有如下两种情况：a)当假设存在n个正向滑移区，n-1个反向滑移区时，将关于未知量及...

【专利技术属性】
技术研发人员：高希光，董洪年，宋迎东，张盛，于国强，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32