基于损伤的P92钢K-R蠕变本构模型拟合方法

Fitting Method of K-R Creep Constitutive Model for P92 Steel Based on Damage

  • 摘要: 为建立适用于反应堆的P92钢高温蠕变本构模型,解决现有方法拟合精度不足的问题,本文通过高温蠕变实验获取不同应力水平下的蠕变数据,采用解析法与多目标遗传算法(MOGA)识别Kachanov-Rabotnov蠕变本构模型(简称K-R模型)材料常数,分析两种方法的局限性。针对模型存在的系统性偏差,进一步构建了一种基于断裂寿命的修正函数以优化预测精度。结果表明:MOGA具有良好的鲁棒性,且有效避免了解析法的误差积累,提高了K-R模型在180 MPa和220 MPa下的拟合精度( R^2 ≥0.93),但两者均无法覆盖200 MPa下的实验数据;修正后的模型在全应力范围拟合精度显著提升,解析法修正后 R^2 ≥0.91,MOGA修正后 R^2 ≥0.97,残差平方和(RSS)降低1个数量级,在超过6倍时间尺度的160 MPa长时外推验证中,修正模型表现出极强的稳定性,断裂寿命预测相对误差为2.45%,蠕变变形 R^2 为0.806。MOGA适用于小样本多目标材料常数识别,修正方法可有效弥补K-R模型固有偏差,建立的本构模型可为反应堆结构完整性评估及其他高温材料相关研究提供参考。

     

    Abstract: P92 steel is a critical structural material for high-temperature components in advanced nuclear reactors and power plants, where its long-term performance under creep conditions is a primary concern for operational safety and structural integrity. The Kachanov-Rabotnov creep constitutive model (K-R model) is a widely accepted framework for describing the creep damage evolution, particularly the tertiary creep stage leading to fracture. However, the predictive accuracy of K-R model is critically dependent on the precise identification of its five material constants, a task that is often challenged by limited experimental data and the inherent limitations of conventional fitting methods. This study aims to develop a comprehensive and high-fidelity framework for establishing the K-R model for P92 steel at 600 ℃, focusing on overcoming the shortcomings of both parameter identification techniques and the model’s intrinsic structural deficiencies. To achieve this objective, a systematic investigation was conducted. First, a series of uniaxial constant-load creep tests were performed on P92 steel specimens at 600 ℃ under three distinct stress levels (180, 200, and 220 MPa) to obtain the complete creep strain-time curves until rupture. The conventional analytical method, which involves a decoupled, step-by-step fitting procedure, was initially employed to identify the K-R model constants. Subsequently, to address the error propagation issues inherent in the analytical approach, a multi-objective genetic algorithm (MOGA) was implemented. This global optimization technique simultaneously identified all five material constants by minimizing the discrepancies between model predictions and experimental data across three objectives: the overall creep curve, the creep rupture time, and the creep rupture strain. Upon discovering a consistent systematic deviation between the model predictions and experimental results even with optimized parameters, a data-driven time-scaling correction methodology was proposed. This method introduced a stress-dependent correction function, derived from the ratio of experimental time to predicted rupture time, to adjust the model’s temporal axis without altering its underlying physical formulation. The results demonstrate a clear hierarchy in the performance of the methods. The analytical method provides a poor fit to the experimental data, yielding R2 values as low as 0.29 and failing to capture the creep behavior accurately due to accumulated errors. In contrast, the MOGA method significantly enhances the parameter identification accuracy and robustness, with R2 for the 180 MPa and 220 MPa stress levels improving to above 0.93. Crucially, both methods reveal a shared systematic deviation at the 200 MPa stress level, suggesting an inherent limitation of the classic K-R model in capturing the specific creep mechanisms of P92 steel under these conditions. The application of the time-scaling correction function proves highly effective. The final, MOGA-based corrected K-R model achieves exceptional predictive accuracy across all tested stress levels. It successfully covers the entire dataset, with R2 consistently exceeding 0.97 and a substantial reduction in the residual sum of squares (RSS), indicating an excellent agreement with the experimental observations. And the model’s extrapolation capability was validated using an independent 160 MPa dataset. For a rupture life of 2 974 h, which exceeds the training set duration by more than six times, the corrected model reduces the life prediction relative error from 24.48% to 2.45%, and improved the R2 from 0.039 to 0.806. In conclusion, this research successfully establishes a robust hybrid framework that integrates an advanced optimization algorithm (MOGA) for reliable parameter identification and a pragmatic correction function to compensate for the constitutive model’s intrinsic limitations. This framework not only delivers a highly accurate and reliable K-R model for P92 steel but also presents a valuable method applicable to other materials and constitutive models facing similar challenges. The developed tool offers a precise basis for the creep life assessment and structural integrity analysis of P92 steel components, thereby contributing to the enhanced design reliability and operational safety of high-temperature engineering systems.

     

/

返回文章
返回