<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">CMC</journal-id>
<journal-id journal-id-type="nlm-ta">CMC</journal-id>
<journal-id journal-id-type="publisher-id">CMC</journal-id>
<journal-title-group>
<journal-title>Computers, Materials &#x0026; Continua</journal-title>
</journal-title-group>
<issn pub-type="epub">1546-2226</issn>
<issn pub-type="ppub">1546-2218</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">69194</article-id>
<article-id pub-id-type="doi">10.32604/cmc.2025.069194</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>An Integrated Approach to Condition-Based Maintenance Decision-Making of Planetary Gearboxes: Combining Temporal Convolutional Network Auto Encoders with Wiener Process</article-title>
<alt-title alt-title-type="left-running-head">An Integrated Approach to Condition-Based Maintenance Decision-Making of Planetary Gearboxes: Combining Temporal Convolutional Network Auto Encoders with Wiener Process</alt-title>
<alt-title alt-title-type="right-running-head">An Integrated Approach to Condition-Based Maintenance Decision-Making of Planetary Gearboxes: Combining Temporal Convolutional Network Auto Encoders with Wiener Process</alt-title>
</title-group>
<contrib-group>
<contrib id="author-1" contrib-type="author">
<name name-style="western"><surname>Zhu</surname><given-names>Bo</given-names></name><xref ref-type="aff" rid="aff-1">1</xref><xref ref-type="author-notes" rid="afn1">#</xref></contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western"><surname>Dong</surname><given-names>Enzhi</given-names></name><xref ref-type="aff" rid="aff-1">1</xref><xref ref-type="author-notes" rid="afn1">#</xref></contrib>
<contrib id="author-3" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Cheng</surname><given-names>Zhonghua</given-names></name><xref ref-type="aff" rid="aff-1">1</xref><email>a15032073178@sina.com</email></contrib>
<contrib id="author-4" contrib-type="author">
<name name-style="western"><surname>Zhan</surname><given-names>Xianbiao</given-names></name><xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-5" contrib-type="author">
<name name-style="western"><surname>Jiang</surname><given-names>Kexin</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-6" contrib-type="author">
<name name-style="western"><surname>Wang </surname><given-names>Rongcai</given-names></name><xref ref-type="aff" rid="aff-3">3</xref></contrib>
<aff id="aff-1"><label>1</label><institution>Shijiazhuang Campus of Army Engineering University of PLA</institution>, <addr-line>Shijiazhuang, 050003</addr-line>, <country>China</country></aff>
<aff id="aff-2"><label>2</label><institution>School of Electronic and Control Engineering, North China Institute of Aerospace Engineering</institution>, <addr-line>Langfang, 065000</addr-line>, <country>China</country></aff>
<aff id="aff-3"><label>3</label><institution>No. 32181 Unit of PLA</institution>, <addr-line>Xi&#x2019;an, 710061</addr-line>, <country>China</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Zhonghua Cheng. Email: <email>a15032073178@sina.com</email></corresp>
<fn id="afn1">
<p><sup>#</sup>These authors contributed equally to this work</p>
</fn>
</author-notes>
<pub-date date-type="collection" publication-format="electronic">
<year>2025</year></pub-date>
<pub-date date-type="pub" publication-format="electronic">
<day>10</day>
<month>11</month>
<year>2025</year>
</pub-date>
<volume>86</volume>
<issue>1</issue>
<fpage>1</fpage>
<lpage>26</lpage>
<history>
<date date-type="received">
<day>17</day>
<month>6</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>30</day>
<month>7</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2025 The Authors.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Published by Tech Science Press.</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_CMC_69194.pdf"></self-uri>
<abstract>
<p>With the increasing complexity of industrial automation, planetary gearboxes play a vital role in large-scale equipment transmission systems, directly impacting operational efficiency and safety. Traditional maintenance strategies often struggle to accurately predict the degradation process of equipment, leading to excessive maintenance costs or potential failure risks. However, existing prediction methods based on statistical models are difficult to adapt to nonlinear degradation processes. To address these challenges, this study proposes a novel condition-based maintenance framework for planetary gearboxes. A comprehensive full-lifecycle degradation experiment was conducted to collect raw vibration signals, which were then processed using a temporal convolutional network autoencoder with multi-scale perception capability to extract deep temporal degradation features, enabling the collaborative extraction of long-period meshing frequencies and short-term impact features from the vibration signals. Kernel principal component analysis was employed to fuse and normalize these features, enhancing the characterization of degradation progression. A nonlinear Wiener process was used to model the degradation trajectory, with a threshold decay function introduced to dynamically adjust maintenance strategies, and model parameters optimized through maximum likelihood estimation. Meanwhile, the maintenance strategy was optimized to minimize costs per unit time, determining the optimal maintenance timing and preventive maintenance threshold. The comprehensive indicator of degradation trends extracted by this method reaches 0.756, which is 41.2% higher than that of traditional time-domain features; the dynamic threshold strategy reduces the maintenance cost per unit time to 55.56, which is 8.9% better than that of the static threshold optimization. Experimental results demonstrate significant reductions in maintenance costs while enhancing system reliability and safety. This study realizes the organic integration of deep learning and reliability theory in the maintenance of planetary gearboxes, provides an interpretable solution for the predictive maintenance of complex mechanical systems, and promotes the development of condition-based maintenance strategies for planetary gearboxes.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>Temporal convolutional network autoencoder</kwd>
<kwd>full lifecycle degradation experiment</kwd>
<kwd>nonlinear Wiener process</kwd>
<kwd>condition-based maintenance decision-making</kwd>
<kwd>fault monitoring</kwd>
</kwd-group>
<funding-group>
<award-group id="awg1">
<funding-source>scientific research projects</funding-source>
<award-id>JY2024B011</award-id>
</award-group>
</funding-group>
</article-meta>
</front>
<body>
<sec id="s1">
<label>1</label>
<title>Introduction</title>
<sec id="s1_1">
<label>1.1</label>
<title>Background</title>
<p>In aerospace, energy development, and industrial manufacturing sectors, planetary gearboxes have emerged as critical drive components in wind turbines, aircraft control systems, machine tools, and industrial robots due to their unique structural advantages and performance characteristics. They are indispensable in applications such as wind turbines [<xref ref-type="bibr" rid="ref-1">1</xref>] and helicopter main rotors [<xref ref-type="bibr" rid="ref-2">2</xref>]. While Industry 4.0 has significantly improved productivity, it has also heightened equipment failure risks [<xref ref-type="bibr" rid="ref-3">3</xref>]. As core components of mechanical transmission systems, planetary gearboxes often operate under heavy loads and high-speed conditions, rendering their key parts susceptible to gear wear, tooth surface fatigue pitting, tooth breakage, and bearing wear [<xref ref-type="bibr" rid="ref-4">4</xref>]. Such failures not only compromise transmission performance but may also lead to severe economic and safety consequences. Thus, designing scientifically sound and efficient maintenance strategies is imperative to reduce operational costs while enhancing system reliability and safety. The current maintenance of planetary gearboxes primarily relies on manufacturer-prescribed fixed intervals. This conventional maintenance approach exhibits significant drawbacks: on one hand, it may lead to &#x201C;over-maintenance&#x201D;, while on the other hand, it could result in &#x201C;under-maintenance&#x201D;, representing one of the key industry pain points. Traditional maintenance strategies based on age or fixed thresholds have been progressively supplanted by CBM decision-making. The widespread adoption of monitoring devices further highlights its broad application prospects in the operation and maintenance of offshore wind turbines [<xref ref-type="bibr" rid="ref-5">5</xref>,<xref ref-type="bibr" rid="ref-6">6</xref>] and gearboxes [<xref ref-type="bibr" rid="ref-7">7</xref>,<xref ref-type="bibr" rid="ref-8">8</xref>]. The CBM employs sensor data acquired through comprehensive condition monitoring systems to establish data-driven reliability models. This proactive strategy enables maintenance personnel to perform continuous health assessments and implement timely, evidence-based maintenance decisions when operational conditions dictate intervention.</p>
<p>With the rapid advancement of information technology and big data analytics, condition monitoring data in modern systems demonstrates increasingly prominent characteristics including high diversity, rapid generation velocity, and massive volume [<xref ref-type="bibr" rid="ref-9">9</xref>]. Traditional machine learning approaches face significant limitations when processing such large-scale datasets, as they excessively rely on signal processing techniques and domain expertise for manual feature extraction. In contrast, deep learning methodologies eliminate the need for manual feature engineering by supporting automated natural feature processing, enabling direct learning of complex feature representations from raw data [<xref ref-type="bibr" rid="ref-10">10</xref>,<xref ref-type="bibr" rid="ref-11">11</xref>]. Following an end-to-end training paradigm from data input to health indicator output, the deep learning effectively circumvents human-induced biases [<xref ref-type="bibr" rid="ref-12">12</xref>]. Current CBM methods typically use statistical models or fault pattern recognition to assess equipment health and predict failures. Unlike traditional approaches relying on manual feature extraction, our TCNAE method automatically learns deep temporal degradation features directly from vibration signals.</p>
<p>The advent of deep learning technology has provided robust technical support for advancing CBM decision-making. Current studies leverage deep neural networks to analyze operational data [<xref ref-type="bibr" rid="ref-13">13</xref>&#x2013;<xref ref-type="bibr" rid="ref-16">16</xref>], extracting features from large-scale sensor measurements to build predictive models. This enables precise equipment degradation forecasting and optimized maintenance strategies [<xref ref-type="bibr" rid="ref-17">17</xref>]. To achieve the practical implementation of deep learning-based CBM technology for planetary gearboxes, this study innovatively establishes the complete lifecycle degradation dataset for planetary gearboxes and develops an end-to-end intelligent maintenance framework that integrates temporal convolutional network autoencoders (TCNAE) with nonlinear Wiener processes, achieving closed-loop optimization from condition monitoring to maintenance decision-making.</p>
</sec>
<sec id="s1_2">
<label>1.2</label>
<title>Related Work</title>
<p>Benefiting from rapid advancements in sensing and condition monitoring technologies, the CBM has gradually attracted significant research attention. The CBM refers to an approach that evaluates current equipment health status by detecting and analyzing condition indicators (such as temperature, pressure, metal content in lubricants, etc.) closely related to equipment health, thereby enabling optimal maintenance decision-making. This approach fundamentally overcomes the limitations of relying solely on expert experience to determine maintenance timing. By effectively monitoring relevant parameters during the degradation process, it becomes possible to infer system reliability information without actual failure occurrence [<xref ref-type="bibr" rid="ref-18">18</xref>&#x2013;<xref ref-type="bibr" rid="ref-21">21</xref>]. Given the stochastic nature of failures, Markov Chains (MC) are frequently employed for effective fault modeling. Naga Srinivasa Rao and Achutha Naikan [<xref ref-type="bibr" rid="ref-22">22</xref>] proposed a Semi-Markov Decision Process (SMDP) to optimize dynamic threshold strategies. Cheng [<xref ref-type="bibr" rid="ref-23">23</xref>] leveraged reinforcement learning to develop a framework for CBM strategy optimization based on ship fatigue damage mechanisms, which has been successfully applied in practice. Kang et al. [<xref ref-type="bibr" rid="ref-24">24</xref>] addressed maintenance uncertainties and weather impacts by introducing a condition-based maintenance policy for offshore wind turbines, optimizing schedules through maintenance bundling and timing to reduce operational expenditures. Zhao et al. [<xref ref-type="bibr" rid="ref-25">25</xref>] demonstrated a conditional opportunistic maintenance strategy for multi-component systems in wind turbines, utilizing Weibull proportional hazard models to characterize individual degradation processes. However, most current CBM research still relies on strong assumptions (e.g., fixed Weibull distribution, linear degradation, etc.). This assumption-driven approach often demonstrates poor adaptability to nonlinear/abrupt degradation patterns (such as stepwise wear in gearboxes), posing significant challenges for practical engineering implementation. Beyond this, research on CBM for gearboxes remains in its nascent stages, with limited studies available in this area.</p>
<p>Accurate health state evaluation constitutes the fundamental premise for implementing effective CBM strategies. Health Indicator (HI) curves, typically exhibiting monotonic trends, enable equipment health assessment through threshold-based evaluation methods [<xref ref-type="bibr" rid="ref-26">26</xref>]. Shen et al. [<xref ref-type="bibr" rid="ref-27">27</xref>] proposed an approach utilizing Adaptive Order Bispectral Slices (AOBS) and Fault Characteristic Energy Ratio (FCER) to investigate vibration signal characteristics. Huang and Chang [<xref ref-type="bibr" rid="ref-28">28</xref>] developed a wireless measurement system to identify vibration signal models featuring Amplitude Modulation (AM) and Frequency Modulation (FM) induced by gear damage. Al-Bedhany et al. [<xref ref-type="bibr" rid="ref-29">29</xref>] and Huang et al. [<xref ref-type="bibr" rid="ref-30">30</xref>] respectively investigated premature failure of wind turbine gearbox bearings, along with the impact of lubricant selection on gearbox efficiency and energy production. Within the context of CBM for gearboxes, current research exhibits a notable gap in literature integrating fault diagnosis with maintenance decision-making, failing to establish a closed-loop system from data-driven analysis to maintenance strategies.</p>
<p>Accurate degradation modeling represents a pivotal element in CBM systems. The Wiener process has demonstrated superior capability in characterizing degradation patterns of various mechanical components (e.g., rolling bearings, lasers, and gyroscopes) during practical maintenance operations, establishing itself as a predominant model for degradation trend analysis [<xref ref-type="bibr" rid="ref-31">31</xref>]. Li et al. [<xref ref-type="bibr" rid="ref-32">32</xref>] developed a two-phase degradation model accounting for bearing performance deterioration characteristics, employing an autoregressive model coupled with nonlinear Wiener processes to delineate performance degradation in each phase, while utilizing extended Kalman filtering to disregard state increment dynamics during state updates. Sun et al. [<xref ref-type="bibr" rid="ref-33">33</xref>] proposed a comprehensive framework for multi-performance, multi-phase Wiener process modeling and reliability analysis, incorporating a biphasic nonlinear Wiener degradation model with Schwarz Information Criterion (SIC) for change-point identification. The Wiener process is defined by its memoryless property, Markovian nature, and Brownian motion, where stochasticity is key. These traits enable broad applications in modeling stochastic degradation. Due to its intuitive physical interpretation and computational efficiency, the Wiener process effectively describes degradation in critical mechanical components (e.g., rolling bearings, lasers, and gyroscopes).</p>
<p>In recent years, deep learning techniques have achieved remarkable progress in the field of the CBM, emerging as a core methodology for prognostic and health management [<xref ref-type="bibr" rid="ref-34">34</xref>]. The inherent uncertainty and stochastic nature of equipment failures have led to widespread adoption of deep learning models for condition assessment and decision-making. These models demonstrate superior capability in processing high-dimensional, nonlinear, and non-stationary sensor data, thereby delivering precise predictive outcomes. Huang et al. [<xref ref-type="bibr" rid="ref-35">35</xref>] proposed a fault diagnosis method based on cloud computing, wavelet transform, and CADRN deep residual networks. By utilizing the channel attention mechanism (CAM), this approach addresses the issues of information redundancy in training samples and the difficulty of conventional networks in extracting subtle fault features. Comparative experiments demonstrated that incorporating the channel attention module into deep residual networks can effectively improve fault diagnosis accuracy. Under different load or rotational speed conditions, the Collaborative Adversarial Domain Adaptation network (CADA) can reduce domain shift through adversarial training, thereby enhancing the generalization ability of degradation feature extraction. To address the limitation that unsupervised domain adaptation (UDA) models struggle to dynamically adjust according to changes in target tasks, Wang et al. [<xref ref-type="bibr" rid="ref-36">36</xref>] proposed a Dynamic Collaborative Adversarial Domain Adaptation Network (DCADAN) for unsupervised fault diagnosis of rotating machinery. The advantage of the Multiscale Lifting Wavelet Contrastive Network (MLWCN) in multi-resolution analysis of vibration signals lies in its ability to capture both high-frequency fault components and low-frequency degradation trends simultaneously. Addressing the issues of insufficient interpretability of deep models and inadequate capability in mining features from small samples, Dong et al. [<xref ref-type="bibr" rid="ref-37">37</xref>] proposed an interpretable multi-scale lifting wavelet contrastive network for fault diagnosis of planetary gearboxes. The Meta-Learning Network with Sensitivity Penalty (MLSP) discusses its robustness with a small amount of labeled data, especially the improved sensitivity to early weak faults in planetary gearboxes. Mu et al. [<xref ref-type="bibr" rid="ref-38">38</xref>] developed an Enhanced Meta-Learning Network with Sensitivity Penalty (EMLN-SP) for the diagnosis of rare faults under severe domain bias. Nevertheless, the opaque &#x201C;black-box&#x201D; nature of deep learning models frequently impedes interpretability in maintenance decision-making, revealing significant gaps in integrating these models with condition information for maintenance strategy formulation.</p>
<p>Current research still exhibits several limitations, which can be categorized into three key issues: First, existing studies generally lack systematic full-lifecycle experimental data support, predominantly relying on fragmented operational data or simplified accelerated aging tests under idealized conditions, which fail to accurately reflect the true degradation patterns in actual complex operating scenarios. Second, there exists a severe disconnection in the workflow from data acquisition to maintenance decision-making, where independent development of each stage leads to substantial information loss, resulting in a significant performance gap between laboratory models and real-world engineering applications. Third, purely data-driven methods suffer from insufficient interpretability, while traditional reliability theories oversimplify actual degradation processes, making current models unable to simultaneously meet the dual requirements of prediction accuracy and engineering applicability.</p>
</sec>
<sec id="s1_3">
<label>1.3</label>
<title>Overview</title>
<p>To address these challenges, this study proposes a novel deep learning-based CBM framework for planetary gearboxes. A comprehensive lifecycle degradation experiment was conducted to collect raw vibration signals, which were processed using the TCNAE to extract deep temporal degradation features. Kernel principal component analysis (KPCA) was employed to fuse and normalize these features, enhancing the degradation process characterization. A nonlinear Wiener process model was implemented to characterize the degradation trajectory, with parameters optimized through maximum likelihood estimation. The maintenance strategy was optimized to minimize costs per unit time, determining both the optimal maintenance timing and preventive maintenance thresholds. The principal innovations of this research can be summarized as follows:
<list list-type="simple">
<list-item><label>(1)</label>
<p><italic>Lifecycle experimentation and high-quality dataset construction</italic>. This study pioneers systematic full lifecycle degradation experiments for planetary gearboxes, acquiring comprehensive time-series vibration data that fills a critical gap in high-quality degradation datasets for this field. The dataset is conducive to discovering the degradation law of the gearbox, providing an essential foundation for data-driven degradation modeling research.</p></list-item>
<list-item><label>(2)</label>
<p><italic>End-to-end intelligent maintenance technology chain</italic>. We innovatively establish a complete technical workflow encompassing &#x201C;condition monitoring-feature extraction-degradation modeling-decision optimization&#x201D;. The framework employs TCNAE for automated feature extraction, integrates KPCA for feature fusion, and utilizes nonlinear Wiener processes for degradation modeling, ultimately forming a closed-loop intelligent maintenance decision system.</p></list-item>
<list-item><label>(3)</label>
<p><italic>Innovative integration of deep learning and traditional reliability theory</italic>. This work represents the first successful combination of deep learning&#x2019;s automated feature extraction capabilities with Wiener process-based reliability theory. This hybrid approach overcomes the limitations of manual feature dependency in conventional methods while preserving model interpretability, achieving both enhanced prediction accuracy and engineering practicality.</p></list-item>
</list></p>
<p>The remainder of this paper is organized as follows: <xref ref-type="sec" rid="s2">Section 2</xref> introduces the TCNAE-based method for extracting degradation features from planetary gearboxes. <xref ref-type="sec" rid="s2_1">Section 2.1</xref> details the parameter estimation procedures. <xref ref-type="sec" rid="s3">Section 3</xref> develops the nonlinear Wiener process-based CBM decision model and details the parameter estimation procedures. <xref ref-type="sec" rid="s4">Section 4</xref> describes the full life-cycle degradation experiment and data processing workflow. <xref ref-type="sec" rid="s5">Section 5</xref> presents the application results integrating the methodologies from <xref ref-type="sec" rid="s3">Sections 3</xref>, <xref ref-type="sec" rid="s4">4</xref>. <xref ref-type="sec" rid="s6">Section 6</xref> concludes the study.</p>
</sec>
</sec>
<sec id="s2">
<label>2</label>
<title>Equipment Degradation Feature Extraction</title>
<sec id="s2_1">
<label>2.1</label>
<title>Procedure</title>
<p>Current research on equipment degradation trend extraction reveals that mechanical systems typically exhibit a biphasic life cycle, as illustrated in <xref ref-type="fig" rid="fig-1">Fig. 1</xref>. The initial phase represents normal operation, followed by a degradation phase. The transition point, termed First Predicting Time (FPT), is identified when the signal trend crosses a predetermined threshold, marking the onset of accelerated degradation. This critical juncture serves as the optimal timing for initiating CBM decision-making. <xref ref-type="fig" rid="fig-2">Fig. 2</xref> presents the systematic implementation workflow for equipment degradation trend extraction methodology.</p>
<fig id="fig-1">
<label>Figure 1</label>
<caption>
<title>Equipment degradation trend extraction process</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-1.tif"/>
</fig><fig id="fig-2">
<label>Figure 2</label>
<caption>
<title>Implementation flowchart of equipment degradation trend extraction method</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-2.tif"/>
</fig>
<p>The proposed methodology enables precise monitoring of equipment health status, facilitating the detection of subtle operational variations by maintenance personnel. This enhanced monitoring capability supports more accurate maintenance intervention decisions. Furthermore, the extracted degradation trends allow for predictive assessment of equipment health progression, thereby providing critical data support for CBM decision-making processes.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Method</title>
<p>This study employs the computationally efficient Fast Fourier Transform (FFT) algorithm to perform spectral transformation of raw vibration signals. The transformed data is input into TCNAE network. The vibration signals of mechanical systems are typical time-series data, containing rich temporal features and local temporal dependencies. The TCNAE structure can automatically capture these local temporal features through convolutional layers. while its autoencoder structure enables effective feature dimensionality reduction and reconstruction, extracting more representative and robust degradation features. Compared with traditional feature extraction methods or some recurrent neural networks (such as LSTM and GRU), TCNAE does not suffer from gradient vanishing or explosion problems when dealing with long sequence data, and achieves higher computational efficiency. It can better meet the needs of degradation feature extraction for complex mechanical systems like planetary gearboxes, providing a high-quality feature basis for subsequent degradation modeling and maintenance decision making.</p>
<p>The TCN was first proposed by Lea et al. in 2016 [<xref ref-type="bibr" rid="ref-39">39</xref>] and has since gained widespread application in time series analysis. The causal convolution, a unidirectional architecture, determines the current output based solely on present and past inputs. Given an input time series <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-2"><mml:math id="mml-ieqn-2"><mml:mi>X</mml:mi><mml:mo>&#x2208;</mml:mo><mml:msup><mml:mi>R</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> and filter <inline-formula id="ieqn-3"><mml:math id="mml-ieqn-3"><mml:mi>F</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>K</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, the causal convolution operation at position <italic>x</italic><sub><italic>T</italic></sub> can be expressed as:
<disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:mi>Y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>T</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:mo>&#x2217;</mml:mo><mml:mi>f</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>T</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:munderover><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>i</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x22C5;</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p>In <xref ref-type="disp-formula" rid="eqn-1">Eq. (1)</xref>, <italic>k</italic> denotes the filter size, representing the number of time steps involved in each convolution operation. The term <italic>T</italic> &#x2212; <italic>i</italic> encodes historical information by considering past temporal directions.</p>
<p><xref ref-type="fig" rid="fig-3">Fig. 3</xref> illustrates the sampling patterns of dilated convolution vs. standard convolution. As shown in <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, the architecture of dilated convolution employs a base dilation rate <italic>d</italic> &#x003D; 1 at the bottom layer, indicating direct sampling of each input element. When <italic>d</italic> &#x003E; 1, the hierarchical structure ensures that even shallow networks can cover extensive input ranges. From <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, given an input time series <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-5"><mml:math id="mml-ieqn-5"><mml:mi>X</mml:mi><mml:mo>&#x2208;</mml:mo><mml:msup><mml:mi>R</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> and filter <inline-formula id="ieqn-6"><mml:math id="mml-ieqn-6"><mml:mi>F</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>K</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, the dilated convolution operation at position <italic>x</italic><sub><italic>T</italic></sub> is formulated as:
<disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:mi>Y</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>T</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>X</mml:mi><mml:msub><mml:mo>&#x2217;</mml:mo><mml:mrow><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mi>f</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>T</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:munderover><mml:mi>F</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>i</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x22C5;</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>d</mml:mi><mml:mi>g</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula>where <italic>d</italic> represents the dilation rate, <italic>k</italic> is the filter size, and <italic>T</italic> &#x2212; <italic>dgi</italic> maintains temporal causality by exclusively encoding past-oriented information.</p>
<fig id="fig-3">
<label>Figure 3</label>
<caption>
<title>Convolution diagrams: (<bold>a</bold>) standard convolution (3 &#x00D7; 3 kernel) (<bold>b</bold>) dilated convolution (3 &#x00D7; 3 kernel, dilation rate of 2)</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-3.tif"/>
</fig><fig id="fig-4">
<label>Figure 4</label>
<caption>
<title>Architecture of dilated convolution</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-4.tif"/>
</fig>
<p>TCN incorporates universal residual modules, whose architecture is depicted in <xref ref-type="fig" rid="fig-5">Fig. 5</xref>. The residual module replaces simple inter layer connections in conventional networks, enhancing model stability during deep layer propagation while improving generalization capability. The residual module&#x2019;s operation can be formalized as:
<disp-formula id="eqn-3"><label>(3)</label><mml:math id="mml-eqn-3" display="block"><mml:mi>O</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi><mml:mi>p</mml:mi><mml:mi>u</mml:mi><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mi>A</mml:mi><mml:mi>c</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>v</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>i</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>+</mml:mo><mml:mi>f</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>Activation</italic> denotes the activation function and <italic>f</italic> represents the transformation performed within the residual module. The residual module allows the network to learn identity mappings more effectively, mitigating issues such as gradient vanishing and degradation in deep networks.</p>
<fig id="fig-5">
<label>Figure 5</label>
<caption>
<title>Residual module</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-5.tif"/>
</fig>
<p>The AE architecture comprises two symmetric components: an encoder that compresses high-dimensional inputs into low-dimensional latent features, and a decoder that reconstructs the original input from these features. In this study, our designed AE processes vibration signals transformed via FFT and progressively compresses the features through multiple layers of dilated causal convolution. This yields compressed bottleneck-layer features containing essential degradation-related information. The decoder then expands these features via a fully connected layer, restores temporal structure using a RepeatVector layer, and fine-tunes temporal correlations with a lightweight single-layer TCN. Finally, a multilayer perceptron (MLP) efficiently reconstructs the input signal, balancing computational cost and processing efficiency.</p>
<p>The TCNAE network&#x2019;s final residual block employs multiple neurons to preserve feature diversity and discriminative power, enabling multidimensional temporal feature extraction. Experimental validation confirms that a single-neuron configuration performs poorly, as unidimensional features cannot effectively capture complex degradation patterns. Since predictive maintenance models require a single degradation trajectory, this study utilizes KPCA [<xref ref-type="bibr" rid="ref-40">40</xref>,<xref ref-type="bibr" rid="ref-41">41</xref>] to fuse multidimensional outputs into a one-dimensional representation, and the polynomial kernel function is adopted to address nonlinearity:
<disp-formula id="eqn-4"><label>(4)</label><mml:math id="mml-eqn-4" display="block"><mml:mi mathvariant="bold">&#x03A6;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext mathvariant="bold">x</mml:mtext></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mtext mathvariant="bold">y</mml:mtext></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B5;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mrow><mml:mtext mathvariant="bold">x</mml:mtext></mml:mrow><mml:mo>&#x22C5;</mml:mo><mml:mrow><mml:mtext mathvariant="bold">y</mml:mtext></mml:mrow><mml:mo>+</mml:mo><mml:mi>&#x03B7;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>&#x03C5;</mml:mi></mml:mrow></mml:msup></mml:math></disp-formula>where <inline-formula id="ieqn-7"><mml:math id="mml-ieqn-7"><mml:mi>&#x03B5;</mml:mi></mml:math></inline-formula> &#x003D; 1 and <inline-formula id="ieqn-8"><mml:math id="mml-ieqn-8"><mml:mi>&#x03B7;</mml:mi></mml:math></inline-formula> &#x003D; 1. The parameter <inline-formula id="ieqn-9"><mml:math id="mml-ieqn-9"><mml:mi>&#x03BD;</mml:mi></mml:math></inline-formula> is determined through distance matrix computation: (1) replacing zero values with infinity, (2) identifying the global minimum distance, (3) computing the average minimum distance, and (4) scaling the average by a factor of 5. The distance matrix itself is obtained by calculating pairwise Euclidean distances between sample points.</p>
<p>The unique temporal features after KPCA processing undergo normalization per <xref ref-type="disp-formula" rid="eqn-5">Eq. (5)</xref>:
<disp-formula id="eqn-5"><label>(5)</label><mml:math id="mml-eqn-5" display="block"><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>x</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:math></disp-formula>where <inline-formula id="ieqn-10"><mml:math id="mml-ieqn-10"><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula id="ieqn-11"><mml:math id="mml-ieqn-11"><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:math></inline-formula> represent the maximum and minimum values of the temporal feature, respectively.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Evaluation Criteria</title>
<p>This study employs monotonicity, correlation, and robustness as evaluation criteria to quantitatively assess the performance of the extracted degradation trend [<xref ref-type="bibr" rid="ref-42">42</xref>]. First, the degradation trend is decomposed into a smooth trend and random error using polynomial fitting:
<disp-formula id="eqn-6"><label>(6)</label><mml:math id="mml-eqn-6" display="block"><mml:mi>H</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>R</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <inline-formula id="ieqn-12"><mml:math id="mml-ieqn-12"><mml:mi>H</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> represents the degradation trend value at time <italic>t</italic><sub><italic>n</italic></sub>, <inline-formula id="ieqn-13"><mml:math id="mml-ieqn-13"><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> denotes the smoothed trend, and <inline-formula id="ieqn-14"><mml:math id="mml-ieqn-14"><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>R</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is the random error.</p>
<p>Monotonicity evaluates whether the degradation trend exhibits a consistent increasing or decreasing pattern. It is calculated as:
<disp-formula id="eqn-7"><label>(7)</label><mml:math id="mml-eqn-7" display="block"><mml:mi>M</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>H</mml:mi><mml:mi>I</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x003E;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mfrac><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x003C;</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mfrac><mml:mo>|</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>Mon(HI)</italic> is the monotonicity index and <italic>m</italic> is the length of the degradation trend vector. A higher <italic>Mon</italic> value indicates a stronger monotonic trend.</p>
<p>The linear correlation between the degradation trend and time is measured using:
<disp-formula id="eqn-8"><label>(8)</label><mml:math id="mml-eqn-8" display="block"><mml:mi>C</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>H</mml:mi><mml:mi>I</mml:mi><mml:mo>,</mml:mo><mml:mi>T</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mo>|</mml:mo><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mrow><mml:mo>(</mml:mo><mml:mi>H</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>H</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>T</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:msup><mml:msqrt><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>H</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>H</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>T</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow></mml:msqrt><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mfrac></mml:math></disp-formula></p>
<p>In <xref ref-type="disp-formula" rid="eqn-8">Eq. (8)</xref>, <italic>Cor(HI)</italic> represents the correlation between the degradation trend and time, where <inline-formula id="ieqn-15"><mml:math id="mml-ieqn-15"><mml:mover><mml:mi>H</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mi>H</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-16"><mml:math id="mml-ieqn-16"><mml:mover><mml:mi>T</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>. A higher <italic>Cor</italic> value indicates a stronger linear relationship.</p>
<p>Robustness evaluates the tolerance of the Health Indicator (HI) to outliers, computed as:
<disp-formula id="eqn-9"><label>(9)</label><mml:math id="mml-eqn-9" display="block"><mml:mi>R</mml:mi><mml:mi>o</mml:mi><mml:mi>b</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>H</mml:mi><mml:mi>I</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mi>R</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>H</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>|</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>Rob(HI)</italic> is the robustness index. A higher <italic>Rob</italic> value indicates better resistance to noise and outliers.</p>
<p>To holistically evaluate degradation trend performance, a Comprehensive Indicator (CI) integrating all three metrics is defined:
<disp-formula id="eqn-10"><label>(10)</label><mml:math id="mml-eqn-10" display="block"><mml:mi>C</mml:mi><mml:mi>I</mml:mi><mml:mo>=</mml:mo><mml:mn>0.4</mml:mn><mml:mo>&#x2217;</mml:mo><mml:mi>M</mml:mi><mml:mi>o</mml:mi><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn>0.3</mml:mn><mml:mo>&#x2217;</mml:mo><mml:mi>C</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mo>+</mml:mo><mml:mn>0.3</mml:mn><mml:mo>&#x2217;</mml:mo><mml:mi>R</mml:mi><mml:mi>o</mml:mi><mml:mi>b</mml:mi></mml:math></disp-formula></p>
<p>By combining these three metrics, the degradation trend can be comprehensively assessed. A well-performing trend accurately reflects the equipment&#x2019;s health state, facilitating precise degradation analysis and CBM decision-making.</p>
</sec>
</sec>
<sec id="s3">
<label>3</label>
<title>CBM Decision-Making Model Based on Nonlinear Wiener</title>
<p>While the proposed TCNAE feature extraction method provides high-precision characterization of equipment degradation, effective CBM decisions require translating these deep temporal features into quantifiable degradation models. Therefore, <xref ref-type="sec" rid="s3">Section 3</xref> establishes a degradation trajectory model based on the nonlinear Wiener process.</p>
<p>The Wiener process, characterized by linear drift and stochastic fluctuations, is inherently suitable for modeling degradation processes that combine deterministic trends with random noise. Its non-monotonic nature accommodates temporary &#x201C;recovery&#x201D; in degradation measurements, making it particularly applicable to multi-phase degradation scenarios. Mechanical system degradation typically exhibits complex nonlinear behavior due to mechanisms like wear, lubrication failure, and fatigue crack propagation. The nonlinear Wiener process captures these characteristics through nonlinear functions describing degradation rate changes, accurately reflecting actual degradation trajectories. Additionally, this process effectively models stochastic influences from environmental fluctuations, load variations, and material defects by combining nonlinear functions with random disturbances, providing a comprehensive simulation of degradation uncertainties.</p>
<p>As industrial equipment grows increasingly complex, the linear Wiener process proves inadequate for characterizing performance evolution. This study adopts a nonlinear Wiener process with adaptive time-scale transformation to better capture dynamic degradation characteristics and failure mechanisms. Through maximum likelihood estimation, we develop an analytical nonlinear degradation model that accurately reflects real-world trends, enabling precise CBM decision-making.</p>
<p>Let <italic>G(t)</italic> denote the degradation amount at time <italic>t</italic>. The degradation process based on Wiener process can be mathematically expressed as:
<disp-formula id="eqn-11"><label>(11)</label><mml:math id="mml-eqn-11" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mi>&#x03C3;</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>&#x03BB;</italic> represents the drift coefficient, <italic>&#x03C3;</italic> denotes the diffusion coefficient, <italic>B(t)</italic> is a standard Brownian motion characterizing the stochastic nature of the degradation process, satisfying <italic>B</italic>(<italic>t</italic>)&#x007E;<italic>N</italic>(0, <italic>t</italic>).</p>
<p>Given the current degradation level <italic>X</italic><sub><italic>k</italic></sub> at time <italic>t</italic><sub><italic>k</italic></sub>, the degradation state at the next monitoring instant <italic>t</italic><sub><italic>k</italic>&#x002B;1</sub> can be expressed as:
<disp-formula id="eqn-12"><label>(12)</label><mml:math id="mml-eqn-12" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mi>&#x03C3;</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>To account for nonlinear characteristics, we incorporate a time-scale transformation function into <xref ref-type="disp-formula" rid="eqn-12">Eq. (12)</xref>, yielding the modified model:
<disp-formula id="eqn-13"><label>(13)</label><mml:math id="mml-eqn-13" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mi mathvariant="normal">&#x03A5;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>;</mml:mo><mml:mi>&#x03C2;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <inline-formula id="ieqn-17"><mml:math id="mml-ieqn-17"><mml:mi mathvariant="normal">&#x03A5;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>;</mml:mo><mml:mi>&#x03C2;</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> is a nonlinear function of time <italic>t</italic>, <inline-formula id="ieqn-18"><mml:math id="mml-ieqn-18"><mml:mi>&#x03C2;</mml:mi></mml:math></inline-formula> represents unknown time-dependent parameters in the nonlinear function.</p>
<p>At the (<italic>k</italic> &#x002B; 1)-th monitoring point, the evolution of the degradation process follows:
<disp-formula id="eqn-14"><label>(14)</label><mml:math id="mml-eqn-14" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>g</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi mathvariant="normal">&#x03A5;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>;</mml:mo><mml:mi>&#x03C2;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi mathvariant="normal">&#x03A5;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>;</mml:mo><mml:mi>&#x03C2;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi mathvariant="normal">&#x03A5;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>;</mml:mo><mml:mi>&#x03C2;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi mathvariant="normal">&#x03A5;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>;</mml:mo><mml:mi>&#x03C2;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>When the drift and diffusion coefficients are assumed to follow non-normal distributions, we adopt the parameter estimation method from reference [<xref ref-type="bibr" rid="ref-43">43</xref>]. Suppose there are <italic>N</italic> sample degradation trajectories, and the inspection time points for each sample are denoted as <italic>t</italic><sub><italic>n</italic>,1</sub>, <italic>t</italic><sub><italic>n</italic>,2</sub>, &#x2026;, <inline-formula id="ieqn-19"><mml:math id="mml-ieqn-19"><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:math></inline-formula>. Here, <italic>m</italic><sub><italic>n</italic></sub> represents the number of degradation values available for the <italic>n</italic>-th sample trajectory, when <italic>n</italic> &#x003D; 1, 2, &#x2026;, <italic>N</italic>. Therefore, the degradation amount of the <italic>n</italic>-th sample trajectory at time <italic>t</italic><sub><italic>n,j</italic></sub> is given by:
<disp-formula id="eqn-15"><label>(15)</label><mml:math id="mml-eqn-15" display="block"><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>&#x22C5;</mml:mo><mml:mi>&#x03B3;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mi>&#x03C3;</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>In <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref>, <inline-formula id="ieqn-20"><mml:math id="mml-ieqn-20"><mml:mi>&#x03B3;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> represents the nonlinear function of time <italic>t</italic> and <inline-formula id="ieqn-21"><mml:math id="mml-ieqn-21"><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>. Let <inline-formula id="ieqn-22"><mml:math id="mml-ieqn-22"><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>, then <inline-formula id="ieqn-23"><mml:math id="mml-ieqn-23"><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>&#x03B3;</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>, while also <inline-formula id="ieqn-24"><mml:math id="mml-ieqn-24"><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>. Here, <bold>X</bold><sub><bold><italic>n</italic></bold></sub> represents the vector of degradation values for the <italic>n</italic>-th sample trajectory, and <inline-formula id="ieqn-25"><mml:math id="mml-ieqn-25"><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mo>&#x22C5;</mml:mo><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> denotes the vector transpose.</p>
<p>To derive the maximum likelihood estimates for parameters <inline-formula id="ieqn-26"><mml:math id="mml-ieqn-26"><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-27"><mml:math id="mml-ieqn-27"><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>, it is essential to formulate the joint probability density function of <inline-formula id="ieqn-28"><mml:math id="mml-ieqn-28"><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-29"><mml:math id="mml-ieqn-29"><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>. According to the theory of multivariate normal distribution, <bold>X</bold><sub><italic>n</italic></sub> obeys a multivariate normal distribution with mean <inline-formula id="ieqn-30"><mml:math id="mml-ieqn-30"><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03BC;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> and variance <inline-formula id="ieqn-31"><mml:math id="mml-ieqn-31"><mml:msub><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, where:
<disp-formula id="eqn-16"><label>(16)</label><mml:math id="mml-eqn-16" display="block"><mml:msub><mml:mrow><mml:mover><mml:mi>&#x03BC;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mtext mathvariant="bold">n</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mi mathvariant="bold">&#x03A9;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula>
<disp-formula id="eqn-17"><label>(17)</label><mml:math id="mml-eqn-17" display="block"><mml:msub><mml:mi mathvariant="bold">&#x03A9;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mtable columnalign="left left left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mo>&#x22EF;</mml:mo></mml:mtd><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mo>&#x22EF;</mml:mo></mml:mtd><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo>&#x22EE;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x22EE;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x22F1;</mml:mo></mml:mtd><mml:mtd><mml:mo>&#x22EE;</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mo>&#x22EF;</mml:mo></mml:mtd><mml:mtd><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>Assuming the measurement values from different degradation trajectories are mutually independent, the joint probability density function for the parameter set <inline-formula id="ieqn-32"><mml:math id="mml-ieqn-32"><mml:mi mathvariant="bold">&#x0398;</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> can be expressed as:
<disp-formula id="eqn-18"><label>(18)</label><mml:math id="mml-eqn-18" display="block"><mml:mrow><mml:mrow><mml:mi>&#x1D4AB;</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi mathvariant="bold">&#x0398;</mml:mi><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mroot><mml:mrow><mml:mn>2</mml:mn><mml:mi>&#x03C0;</mml:mi></mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mroot></mml:mfrac><mml:msup><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:mrow></mml:msup><mml:mi>exp</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msup><mml:msubsup><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>Taking the natural logarithm of both sides of <xref ref-type="disp-formula" rid="eqn-18">Eq. (18)</xref>, we obtain the log-likelihood function:
<disp-formula id="eqn-19"><label>(19)</label><mml:math id="mml-eqn-19" display="block"><mml:mrow><mml:mtext>In</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mi>&#x1D4AB;</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi mathvariant="normal">&#x0398;</mml:mi><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mtext>In</mml:mtext></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mn>2</mml:mn><mml:mi>&#x03C0;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mtext>In</mml:mtext></mml:mrow><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi mathvariant="normal">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msup><mml:msubsup><mml:mi mathvariant="normal">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>By taking partial derivatives with respect to <inline-formula id="ieqn-33"><mml:math id="mml-ieqn-33"><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-34"><mml:math id="mml-ieqn-34"><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>, respectively, we derive:
<disp-formula id="eqn-20"><label>(20)</label><mml:math id="mml-eqn-20" display="block"><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mtext>In</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mi>&#x1D4AB;</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi mathvariant="bold">&#x0398;</mml:mi><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:msubsup><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msubsup><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msubsup><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msubsup><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula>
<disp-formula id="eqn-21"><label>(21)</label><mml:math id="mml-eqn-21" display="block"><mml:mfrac><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:mrow><mml:mtext>In</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mi>&#x1D4AB;</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi mathvariant="bold">&#x0398;</mml:mi><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mi>N</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>2</mml:mn><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msup><mml:msup><mml:mi mathvariant="normal">&#x03A9;</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>Setting <xref ref-type="disp-formula" rid="eqn-20">Eq. (20)</xref> to zero, the maximum likelihood estimate of the drift coefficient <inline-formula id="ieqn-35"><mml:math id="mml-ieqn-35"><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> is obtained as:
<disp-formula id="eqn-22"><label>(22)</label><mml:math id="mml-eqn-22" display="block"><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msubsup><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msubsup><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msubsup><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msup><mml:mi mathvariant="normal">&#x03A9;</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msubsup><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup><mml:msup><mml:mi mathvariant="normal">&#x03A9;</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>Similarly, setting <xref ref-type="disp-formula" rid="eqn-21">Eq. (21)</xref> to zero, the maximum likelihood estimate of the diffusion coefficient <inline-formula id="ieqn-36"><mml:math id="mml-ieqn-36"><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> is obtained as:
<disp-formula id="eqn-23"><label>(23)</label><mml:math id="mml-eqn-23" display="block"><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msup><mml:msup><mml:mi mathvariant="normal">&#x03A9;</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mrow><mml:mtext mathvariant="bold">X</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>In <xref ref-type="disp-formula" rid="eqn-20">Eqs. (20)</xref> and <xref ref-type="disp-formula" rid="eqn-21">(21)</xref>, <inline-formula id="ieqn-37"><mml:math id="mml-ieqn-37"><mml:msubsup><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2032;</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> represents the transpose of the matrix <inline-formula id="ieqn-38"><mml:math id="mml-ieqn-38"><mml:msub><mml:mrow><mml:mtext mathvariant="bold">V</mml:mtext></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>, and <inline-formula id="ieqn-39"><mml:math id="mml-ieqn-39"><mml:msubsup><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula> denotes the inverse of the vector matrix <inline-formula id="ieqn-40"><mml:math id="mml-ieqn-40"><mml:msub><mml:mi mathvariant="bold">&#x03A3;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>.</p>
<p>Based on the equipment degradation trend curves presented in <xref ref-type="sec" rid="s2">Section 2</xref>, we establish four distinct degradation modeling phases as follows:</p>
<p>When <italic>t</italic><sub><italic>k</italic></sub> &#x2264; <italic>t</italic><sub>1</sub>, the degradation process is modeled by:
<disp-formula id="eqn-24"><label>(24)</label><mml:math id="mml-eqn-24" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo>[</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>s</italic> denotes the dummy variable for integration, representing the integration variable when differentiating Brownian motion over the time domain, <italic>a</italic><sub>1</sub>, <italic>a</italic><sub>2</sub>, <italic>a</italic><sub>3</sub>, <italic>a</italic><sub>4</sub> &#x003E; 0. The drift and diffusion terms are respectively:
<disp-formula id="eqn-25"><label>(25)</label><mml:math id="mml-eqn-25" display="block"><mml:mtable columnalign="right left right left right left right left right left right left" rowspacing="3pt" columnspacing="0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em" displaystyle="true"><mml:mtr><mml:mtd /><mml:mtd><mml:msub><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BC;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msup><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="eqn-26"><label>(26)</label><mml:math id="mml-eqn-26" display="block"><mml:mtable columnalign="right left right left right left right left right left right left" rowspacing="3pt" columnspacing="0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em" displaystyle="true"><mml:mtr><mml:mtd /><mml:mtd><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p>
<p>When <italic>t</italic><sub>1</sub> &#x003C; <italic>t</italic><sub><italic>k</italic></sub> &#x2264; <italic>t</italic><sub>2</sub>, the degradation process is modeled by:
<disp-formula id="eqn-27"><label>(27)</label><mml:math id="mml-eqn-27" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mi>arctan</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>5</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where coefficients <italic>p</italic><sub>1</sub>, <italic>p</italic><sub>2</sub>, <italic>p</italic><sub>3</sub>, <italic>p</italic><sub>4</sub>, <italic>p</italic><sub>5</sub> &#x003E; 0. The corresponding drift and diffusion terms are:
<disp-formula id="eqn-28"><label>(28)</label><mml:math id="mml-eqn-28" display="block"><mml:msub><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BC;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:mfrac><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>[</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>]</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>+</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>5</mml:mn></mml:mrow></mml:msub><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>p</mml:mi><mml:mrow><mml:mn>5</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>
<disp-formula id="eqn-29"><label>(29)</label><mml:math id="mml-eqn-29" display="block"><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>When <italic>t</italic><sub>2</sub> &#x003C; <italic>t</italic><sub><italic>k</italic></sub> <italic>&#x2264; t</italic><sub>3</sub>, the degradation process is modeled by:
<disp-formula id="eqn-30"><label>(30)</label><mml:math id="mml-eqn-30" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where coefficients <italic>q</italic><sub>1</sub>, <italic>q</italic><sub>2</sub>, <italic>q</italic><sub>3</sub>, <italic>q</italic><sub>4</sub> &#x003E; 0. The drift and diffusion terms are:
<disp-formula id="eqn-31"><label>(31)</label><mml:math id="mml-eqn-31" display="block"><mml:mtable columnalign="right left right left right left right left right left right left" rowspacing="3pt" columnspacing="0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em" displaystyle="true"><mml:mtr><mml:mtd /><mml:mtd><mml:msub><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BC;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>q</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="eqn-32"><label>(32)</label><mml:math id="mml-eqn-32" display="block"><mml:mtable columnalign="right left right left right left right left right left right left" rowspacing="3pt" columnspacing="0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em" displaystyle="true"><mml:mtr><mml:mtd /><mml:mtd><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p>
<p>When <italic>t</italic><sub>3</sub> &#x003C; <italic>t</italic><sub><italic>k</italic></sub> <italic>&#x2264; t</italic><sub>4</sub>, the degradation process is modeled by:
<disp-formula id="eqn-33"><label>(33)</label><mml:math id="mml-eqn-33" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:msup><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msubsup><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where coefficients <italic>r</italic><sub>1</sub>, <italic>r</italic><sub>2</sub>, <italic>r</italic><sub>3</sub>, <italic>r</italic><sub>4</sub> &#x003E; 0. The drift and diffusion terms are:
<disp-formula id="eqn-34"><label>(34)</label><mml:math id="mml-eqn-34" display="block"><mml:mtable columnalign="right left right left right left right left right left right left" rowspacing="3pt" columnspacing="0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em" displaystyle="true"><mml:mtr><mml:mtd /><mml:mtd><mml:msub><mml:mi>&#x03BC;</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03BC;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>r</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<disp-formula id="eqn-35"><label>(35)</label><mml:math id="mml-eqn-35" display="block"><mml:mtable columnalign="right left right left right left right left right left right left" rowspacing="3pt" columnspacing="0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em" displaystyle="true"><mml:mtr><mml:mtd /><mml:mtd><mml:msub><mml:mi>&#x03C3;</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:msup><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula></p>
<p>The developed Wiener process with nonlinear structure effectively captures actual degradation trends, yielding precise degradation trajectories.</p>
<p>This study investigates the degradation performance of planetary gearbox systems, employing variable-period inspection as an indicator for equipment performance degradation [<xref ref-type="bibr" rid="ref-44">44</xref>&#x2013;<xref ref-type="bibr" rid="ref-46">46</xref>]. Let <italic>G(t)</italic> denote the system&#x2019;s degradation amount at time <italic>t</italic>, with the following assumptions:
<list list-type="simple">
<list-item><label>(1)</label>
<p>The initial degradation amount is 0, with a continuous and monotonic degradation trend;</p></list-item>
<list-item><label>(2)</label>
<p>Inspections do not alter the system&#x2019;s degradation state, and inspection duration is negligible;</p></list-item>
<list-item><label>(3)</label>
<p>The degradation rate accelerates with both operational duration and preventive maintenance frequency;</p></list-item>
<list-item><label>(4)</label>
<p>After preventive maintenance, the system&#x2019;s health state returns to its initial condition.</p></list-item>
</list></p>
<p>Under these assumptions, we classify system states based on inspection interval <inline-formula id="ieqn-41"><mml:math id="mml-ieqn-41"><mml:mi>&#x03D1;</mml:mi></mml:math></inline-formula> and preventive maintenance threshold <italic>Th</italic><sub><italic>pm</italic></sub>.
<list list-type="simple">
<list-item><label>(1)</label>
<p>If <inline-formula id="ieqn-42"><mml:math id="mml-ieqn-42"><mml:mi>G</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x003C;</mml:mo><mml:mi>T</mml:mi><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>: The system operates normally, incurring only inspection cost <italic>C</italic><sub><italic>in</italic></sub>;</p></list-item>
<list-item><label>(2)</label>
<p>If <inline-formula id="ieqn-43"><mml:math id="mml-ieqn-43"><mml:mi>T</mml:mi><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> <italic>&#x2264; X</italic>(<italic>t</italic><sub><italic>i</italic></sub>) &#x003C; <italic>Th</italic><sub><italic>cm</italic></sub>: The degradation exceeds preventive maintenance threshold but remains below failure threshold, requiring preventive maintenance with combined costs (<italic>C</italic><sub><italic>in</italic></sub> &#x002B; <italic>C</italic><sub><italic>pm</italic></sub>);</p></list-item>
<list-item><label>(3)</label>
<p>If <inline-formula id="ieqn-44"><mml:math id="mml-ieqn-44"><mml:mi>G</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x003E;</mml:mo><mml:mi>T</mml:mi><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>: The degradation surpasses failure threshold, necessitating corrective maintenance with total costs (<italic>C</italic><sub><italic>in</italic></sub> &#x002B; <italic>C</italic><sub><italic>cm</italic></sub> &#x002B; <italic>C</italic><sub><italic>puni</italic></sub>), where <italic>C</italic><sub><italic>puni</italic></sub> represents the time-dependent penalty cost.</p></list-item>
</list></p>
<p>From an economic perspective, we establish the maintenance cost rate <italic>CR</italic> over the system&#x2019;s total life cycle <italic>&#x03C4;</italic> by minimizing the objective function:
<disp-formula id="eqn-36"><label>(36)</label><mml:math id="mml-eqn-36" display="block"><mml:mo movablelimits="true" form="prefix">min</mml:mo><mml:mi>E</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>C</mml:mi><mml:mi>R</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mo movablelimits="true" form="prefix">min</mml:mo><mml:mfrac><mml:mrow><mml:mi>E</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>C</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>E</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C4;</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>The total maintenance cost <italic>C(t)</italic> of the planetary gearbox system can be mathematically expressed as:
<disp-formula id="eqn-37"><label>(37)</label><mml:math id="mml-eqn-37" display="block"><mml:mi>C</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:munderover><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:munderover><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:munderover><mml:mo>&#x2211;</mml:mo><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:munderover><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>u</mml:mi><mml:mi>n</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>C</italic><sub><italic>in</italic></sub>, <italic>C</italic><sub><italic>pm</italic></sub>, <italic>C</italic><sub><italic>cm</italic></sub> denote inspection, preventive maintenance, and corrective maintenance cumulative counts respectively, &#x0394;<italic>t</italic><sub><italic>k</italic></sub> represents the downtime duration for the <italic>k</italic>-th corrective maintenance.</p>
<p>Considering the actual situation that the preventive maintenance threshold decreases with the increase of maintenance times in maintenance activities, a threshold decay function is established as follows:
<disp-formula id="eqn-38"><label>(38)</label><mml:math id="mml-eqn-38" display="block"><mml:msub><mml:mi>&#x03F1;</mml:mi><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x03F1;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&#x22C5;</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03C1;</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msup><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>&#x03B6;</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>In <xref ref-type="disp-formula" rid="eqn-38">Eq. (38)</xref>, <italic>n</italic> represents the number of maintenance times, <inline-formula id="ieqn-45"><mml:math id="mml-ieqn-45"><mml:msub><mml:mi>&#x03F1;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> represents the initial preventive maintenance threshold, <italic>&#x03C1;</italic> represents the decay coefficient, and <inline-formula id="ieqn-46"><mml:math id="mml-ieqn-46"><mml:mi>&#x03B6;</mml:mi></mml:math></inline-formula> denotes the nonlinear factor.</p>
<p>The objective function derived from <xref ref-type="disp-formula" rid="eqn-36">Eqs. (36)</xref> and <xref ref-type="disp-formula" rid="eqn-37">(37)</xref> is formulated as:
<disp-formula id="eqn-39"><label>(39)</label><mml:math id="mml-eqn-39" display="block"><mml:mrow><mml:mo>(</mml:mo><mml:msup><mml:mi>&#x03D1;</mml:mi><mml:mrow><mml:mo>&#x2217;</mml:mo></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mi>T</mml:mi><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x2217;</mml:mo></mml:mrow></mml:msubsup><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>arg</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mo movablelimits="true" form="prefix">min</mml:mo><mml:mi>E</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>C</mml:mi><mml:mi>R</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03D1;</mml:mi><mml:mo>,</mml:mo><mml:mi>T</mml:mi><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mi>s</mml:mi><mml:mo>.</mml:mo><mml:mi>t</mml:mi><mml:mo>.</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mi>&#x03D1;</mml:mi><mml:mo>&#x003E;</mml:mo><mml:mn>0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>T</mml:mi><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>&#x003E;</mml:mo><mml:mi>T</mml:mi><mml:msub><mml:mi>h</mml:mi><mml:mrow><mml:mi>p</mml:mi><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>&#x003E;</mml:mo><mml:mn>0</mml:mn></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>As indicated in <xref ref-type="disp-formula" rid="eqn-39">Eq. (39)</xref>, the Condition-Based Maintenance decision model employs simulation methods to analyze and solve for the optimal decision variables.</p>
</sec>
<sec id="s4">
<label>4</label>
<title>Case Study</title>
<p>To validate the practical utility of the theoretical model, <xref ref-type="sec" rid="s4">Section 4</xref> presents a full life-cycle experiment on planetary gearboxes, where vibration signals are collected to implement and verify the TCNAE feature extraction and Wiener modeling methods proposed earlier, thereby providing data support for subsequent maintenance decision optimization.</p>
<sec id="s4_1">
<label>4.1</label>
<title>Experimental Setup</title>
<p>A 900-h full-life degradation experiment was conducted, and the main components of the experimental platform are as follows: the main motor is a three-speed electromagnetic speed-regulating motor (model YCT180-4A) with a rated power of 4 kW; the speed-torque sensor (model JN338) is used to collect the speed and torque signals of the input shaft of the planetary gearbox; the magnetic powder brake (model FZJ-5) provides load to the planetary gearbox, and the SC-1D (3A) load controller adjusts the load value of the magnetic powder brake by controlling the current [<xref ref-type="bibr" rid="ref-47">47</xref>], as shown in <xref ref-type="fig" rid="fig-6">Fig. 6</xref>.</p>
<fig id="fig-6">
<label>Figure 6</label>
<caption>
<title>Experimental platform of planetary gearbox (1&#x2014;main motor, 2&#x2014;speed and torque sensor, 3&#x2014;planetary gearbox, 4&#x2014;magnetic powder brake, 5&#x2014;base, 6&#x2014;coupling)</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-6.tif"/>
</fig>
<p>The test employed an NGW11 single-stage planetary gearbox, with its internal structure schematic shown in <xref ref-type="fig" rid="fig-7">Fig. 7</xref> and detailed component parameters provided in <xref ref-type="table" rid="table-1">Table 1</xref>. During operation, the gearbox achieved a transmission ratio of 12.5. Four vibration acceleration sensors were strategically mounted on the gearbox housing at locations depicted in <xref ref-type="fig" rid="fig-8">Fig. 8</xref>, enabling comprehensive vibration signal acquisition from multiple orientations to effectively monitor degradation progression.</p>
<fig id="fig-7">
<label>Figure 7</label>
<caption>
<title>Internal structure diagram of the planetary gearbox</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-7.tif"/>
</fig><table-wrap id="table-1">
<label>Table 1</label>
<caption>
<title>Key component parameters of planetary gearbox</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th>Main components</th>
<th>Sun gear</th>
<th>Planetary gear</th>
<th>Ring gear</th>
<th>Planet carrier</th>
<th>Speed ratio</th>
</tr>
</thead>
<tbody>
<tr>
<td>Number of teeth</td>
<td>13</td>
<td>64</td>
<td>146</td>
<td>/</td>
<td rowspan="2">12.5</td>
</tr>
<tr>
<td>Number of plates</td>
<td>1</td>
<td>3</td>
<td>1</td>
<td>1</td>
</tr>
</tbody>
</table>
</table-wrap><fig id="fig-8">
<label>Figure 8</label>
<caption>
<title>Sensor distribution locations</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-8.tif"/>
</fig>
<p>The main operating parameters of the experiment are as follows: (1) Input shaft speed: 1000 revolutions per minute (rpm). (2) Load current applied to the magnetic powder brake: 1 ampere (A), equivalent to a load of approximately 340 Newton-meters (N&#x00B7;m) on the planetary gearbox. (3) Vibration signal sampling frequency: 20 kHz, with a single sampling duration of 1 s and a sampling interval of 5 min. (4) The ambient temperature in the laboratory is controlled at 20 &#x00B1; 1&#x00B0;C, and the relative humidity is maintained at 55 &#x00B1; 5%.</p>
<p>This study selects the sample data collected by Sensor 2 for feature analysis. The planetary gearbox dataset obtained from the PHM laboratory consists of 900 CSV files with a total storage capacity of approximately 2.09 GB. Each file is sequentially named from data-T155412Z to data-T161513Z in chronological order, precisely recording continuous monitoring data throughout the experimental process to ensure the integrity of signal details. Every CSV file contains vibration acceleration data organized in a single column with 240,000 rows, capturing the vertical vibration signals of the planetary gearbox with a precision of five decimal places.</p>
<p>An electromagnetic brake provided a stable load during the experiment, with a load controller precisely adjusting the current to regulate braking torque. The magnetic powder brake operates via electromagnetic effects: when the controller supplies a 1A DC current, the magnetic powder forms chains, coupling the inner and outer rotors to produce a braking torque proportional to the current. To minimize instability effects, vibration signals were collected only during steady-state operation, excluding transient data from motor start-up and shutdown phases.</p>
<p>Following the complete life-cycle testing, the planetary gearbox exhibited multiple failure modes across its primary components, including inner and outer race wear, cracking damage, and varying degrees of tooth surface wear on major gears, as shown in <xref ref-type="fig" rid="fig-9">Fig. 9</xref>. From <xref ref-type="fig" rid="fig-9">Fig. 9</xref>, it can be observed that there are obvious signs of fatigue pitting and wear on the tooth surface of the sun gear. A small number of tooth surfaces on the planet gears exhibit wear and local fractures. The tooth surface wear of the ring gear is relatively uniform, but fracture traces were found in some areas. These damage patterns match typical planetary gearbox degradation under prolonged high-load operation, confirming the experimental data&#x2019;s reliability.</p>
<fig id="fig-9">
<label>Figure 9</label>
<caption>
<title>Schematic diagram of wear components in planetary gearbox</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-9.tif"/>
</fig>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Data Processing</title>
<sec id="s4_2_1">
<label>4.2.1</label>
<title>TCNAE Model Training</title>
<p>After multiple pre-training sessions with the TCNAE model and referencing existing literature methods, the optimal hyperparameters are determined as follows: the encoding part consists of a stacked architecture with 3 residual blocks and 2 fully connected layers, where the number of neurons in the 3 residual blocks is 256, 128, and 64, respectively. Based on the time resolution of vibration signals and the requirement for receptive fields, the kernel sizes are set to 7, 5, and 3, which can capture long-period vibrations (such as meshing frequencies), obtain medium-frequency features, and extract transient features (such as impact signals), respectively, with reference to the practical applications of TCN in mechanical signal processing as described in [<xref ref-type="bibr" rid="ref-39">39</xref>].</p>
<p>To cover multi-scale temporal features, the dilation rates are set to 1, 2, and 4, respectively. Since the vibration signals of gearboxes contain high-frequency bearing fault components and low-frequency gear wear components, the design of exponentially increasing dilation rates can cover different time scales. To address the issue that vibration signals involve positive and negative oscillations, and the negative information is prone to loss with the traditional ReLU, we have selected the Leaky ReLU activation function, which adapts to the bidirectional oscillation characteristics of vibration signals by retaining the information on the negative semi-axis and alleviates the problem of gradient vanishing.</p>
<p>The two convolutional layers within each residual block undergo kernel initialization, and layer normalization is used for normalization processing. The decoder part is composed of a residual block with 64 neurons (dilation rate &#x003D; 1) and an MLP network. Feature integrity is maintained through skip connections. An Adam optimizer with an initial learning rate of 0.001 and an MSE loss function are adopted to form the baseline for the regression task. The batch size is set to 128 to balance GPU memory usage and training efficiency, while the number of training epochs is set to 100 to ensure loss convergence. To avoid overfitting, Dropout is set to 0.25, and shortcut connections are included. The training status of the TCNAE model is shown in <xref ref-type="fig" rid="fig-10">Fig. 10</xref>, which indicates that the MSE curve gradually declines and tends to stabilize as the epoch increases, with the final loss value remaining at 0.0163.</p>
<fig id="fig-10">
<label>Figure 10</label>
<caption>
<title>Training status of the TCNAE network model</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-10.tif"/>
</fig>
<p>For high-dimensional time-series data like the 16,383-dimensional vibration signals in our gearbox dataset, TCNAE demonstrates superior suitability for industrial feature extraction tasks owing to its low parameter requirements, high parallelism, and efficient memory utilization. However, LSTM/GRU struggle to handle high dimensional data due to memory limitations and sequence calculation characteristics.</p>
</sec>
<sec id="s4_2_2">
<label>4.2.2</label>
<title>Degradation Trend Extraction</title>
<p>The TCNAE model autonomously extracted 64 deep temporal features characterizing planetary gearbox degradation. Among these, 53 features exhibited values below zero, resulting in 11 valid deep temporal features. These valid features were smoothed using Savitzky-Golay filtering (window length &#x003D; 301, polynomial order &#x003D; 30) to reduce noise, as shown in <xref ref-type="fig" rid="fig-11">Fig. 11</xref>.</p>
<fig id="fig-11">
<label>Figure 11</label>
<caption>
<title>Effective deep temporal features extracted by the TCNAE network model</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-11.tif"/>
</fig>
<p>The KPCA method described in 2.2 fused effective deep time-series features, with a polynomial function. After normalization, the fused deep time-series features are used as the final degradation trend of the planetary gearbox throughout its entire life cycle, as presented in <xref ref-type="fig" rid="fig-12">Fig. 12</xref>.</p>
<fig id="fig-12">
<label>Figure 12</label>
<caption>
<title>Fused deep temporal features</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-12.tif"/>
</fig>
</sec>
<sec id="s4_2_3">
<label>4.2.3</label>
<title>Comparison of Degradation Trends</title>
<p>As a traditional approach for degradation trend construction, statistical parameter-based methods analyze time-domain features from vibration signals [<xref ref-type="bibr" rid="ref-48">48</xref>]. In planetary gearbox degradation, signal amplitude and probability distribution variations directly affect these parameters, which effectively characterize vibration features. These statistical parameters can intuitively characterize the features of vibration signals. <xref ref-type="table" rid="table-2">Table 2</xref> lists six key time-domain features among statistical parameters, with their extracted trends displayed in <xref ref-type="fig" rid="fig-13">Fig. 13</xref>. Here, <italic>x</italic><sub><italic>i</italic></sub> is the sampled time signal, <italic>i</italic> is the sample index, and <italic>N</italic> is the number of samples.</p>
<table-wrap id="table-2">
<label>Table 2</label>
<caption>
<title>Formulas of 6 time-domain feature parameters</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th>Characteristic parameters</th>
<th>Formula</th>
<th>Characteristic parameters</th>
<th>Formula</th>
</tr>
</thead>
<tbody>
<tr>
<td>Mean value</td>
<td><inline-formula id="ieqn-47"><mml:math id="mml-ieqn-47"><mml:mover><mml:mi>x</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula></td>
<td>Root mean square</td>
<td><inline-formula id="ieqn-48"><mml:math id="mml-ieqn-48"><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>M</mml:mi><mml:mi>S</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>/</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:math></inline-formula></td>
</tr>
<tr>
<td>Margin factor</td>
<td><inline-formula id="ieqn-49"><mml:math id="mml-ieqn-49"><mml:mfrac><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:mrow><mml:msup><mml:mrow><mml:mo>|</mml:mo><mml:mover><mml:mi>x</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>|</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mfrac></mml:math></inline-formula></td>
<td>Peak to Peak value</td>
<td><inline-formula id="ieqn-50"><mml:math id="mml-ieqn-50"><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:math></inline-formula></td>
</tr>
<tr>
<td>Kurtosis</td>
<td><inline-formula id="ieqn-51"><mml:math id="mml-ieqn-51"><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>x</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msup><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>M</mml:mi><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msubsup></mml:mfrac></mml:math></inline-formula></td>
<td>Skewness</td>
<td><inline-formula id="ieqn-52"><mml:math id="mml-ieqn-52"><mml:mfrac><mml:mn>1</mml:mn><mml:mi>N</mml:mi></mml:mfrac><mml:munderover><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:munderover><mml:mfrac><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mover><mml:mi>x</mml:mi><mml:mo accent="false">&#x00AF;</mml:mo></mml:mover><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>M</mml:mi><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msubsup></mml:mfrac></mml:math></inline-formula></td>
</tr>
</tbody>
</table>
</table-wrap><fig id="fig-13">
<label>Figure 13</label>
<caption>
<title>Time-domain characteristics of vibration signals over the full life cycle of the planetary gearbox</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-13.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Results</title>
<p>Based on processed experimental metrics, this section systematically evaluates the feature extraction performance of TCNAE, parameter estimation accuracy, and the superiority of dynamic threshold strategies, offering comprehensive assessment of the framework&#x2019;s engineering value.</p>
<sec id="s5_1">
<label>5.1</label>
<title>Data Analysis</title>
<p>According to the experimental protocol outlined in <xref ref-type="sec" rid="s4">Section 4</xref>, the acceleration values of the vibration signals in the vertical direction of the power input shaft are presented in <xref ref-type="fig" rid="fig-14">Fig. 14</xref>. As shown in <xref ref-type="fig" rid="fig-14">Fig. 14</xref>, the planetary gearbox exhibited stable operation with minor fluctuations, followed by progressive degradation. A sharp signal change at 750 h indicated component failure, culminating in severe vibration and audible noise at 900 h that marked complete gearbox failure.</p>
<fig id="fig-14">
<label>Figure 14</label>
<caption>
<title>Full life cycle vibration signals of the planetary gearbox-sensor 2</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-14.tif"/>
</fig>
<p>Comparative analysis of <xref ref-type="fig" rid="fig-15">Fig. 15</xref> and <xref ref-type="table" rid="table-3">Table 3</xref> clearly demonstrates TCNAE&#x2019;s superior performance over LSTM and GRU in three key aspects. First, it exhibits stronger feature representation capability, with normalized feature values significantly higher than LSTM/GRU across most sample intervals, indicating more effective capture of vibration signal patterns and higher compression efficiency for high-dimensional FFT signals in its hidden states. Second, it shows better temporal stability, with notably smaller curve fluctuations compared to the pronounced oscillations observed in LSTM/GRU. The fixed receptive field in TCN stabilizes long range feature extraction, making it more suitable for stable industrial monitoring. Third, it demonstrates higher fault sensitivity, displaying significant feature value increases in tail sample intervals through multi-scale convolutional kernels that amplify fault-band energy, whereas LSTM/GRU shows weak responses.</p>
<fig id="fig-15">
<label>Figure 15</label>
<caption>
<title>Comparison of time-series features trained by three models</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-15.tif"/>
</fig><table-wrap id="table-3">
<label>Table 3</label>
<caption>
<title>Comparison of indicators among three models</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th>Indicators</th>
<th>TCNAE</th>
<th>LSTM</th>
<th>GRU</th>
</tr>
</thead>
<tbody>
<tr>
<td>Average feature value (in the 600&#x2013;900 interval)</td>
<td>0.72</td>
<td>0.51</td>
<td>0.48</td>
</tr>
<tr>
<td>Fluctuation standard deviation (&#x00D7;10<sup>&#x2212;2</sup>)</td>
<td>1.81</td>
<td>3.25</td>
<td>2.95</td>
</tr>
<tr>
<td>Amplitude of increase in the fault interval</td>
<td>0.25</td>
<td>0.08</td>
<td>0.06</td>
</tr>
<tr>
<td>Performance comparison</td>
<td>0.039</td>
<td>0.027</td>
<td>0.012</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Using the comprehensive evaluation criteria for degradation trend performance proposed in <xref ref-type="sec" rid="s2_3">Section 2.3</xref>, the CI values (also known as health factors) of the degradation trends derived from six time-domain features and the TCNAE method were calculated, as shown in <xref ref-type="table" rid="table-4">Table 4</xref>.</p>
<table-wrap id="table-4">
<label>Table 4</label>
<caption>
<title>Various indicators reflecting degradation trends</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th>Indicators</th>
<th>Mon</th>
<th>Cor</th>
<th>Rob</th>
<th>CI</th>
</tr>
</thead>
<tbody>
<tr>
<td>CI</td>
<td>0.551</td>
<td>0.781</td>
<td>0.998</td>
<td>0.756</td>
</tr>
<tr>
<td>Mean value</td>
<td>0.258</td>
<td>0.201</td>
<td>0.999</td>
<td>0.463</td>
</tr>
<tr>
<td>Margin factor</td>
<td>0.010</td>
<td>0.235</td>
<td>0.744</td>
<td>0.294</td>
</tr>
<tr>
<td>Root mean square</td>
<td>0.227</td>
<td>0.865</td>
<td>0.976</td>
<td>0.643</td>
</tr>
<tr>
<td>Kurtosis</td>
<td>0.044</td>
<td>&#x2212;0.307</td>
<td>0.967</td>
<td>0.216</td>
</tr>
<tr>
<td>Peak to Peak value</td>
<td>0.025</td>
<td>0.744</td>
<td>0.334</td>
<td>0.333</td>
</tr>
<tr>
<td>Skewness</td>
<td>0.038</td>
<td>0.369</td>
<td>0.993</td>
<td>0.424</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>As indicated in <xref ref-type="table" rid="table-4">Table 4</xref>, most time-domain features exhibit severe fluctuations, while the degradation trend extracted by the TCNAE method yields the highest CI value, demonstrating the best characterization effect for the degradation trend of the planetary gearbox.</p>

<p>Based on the parameter estimation method described in <xref ref-type="sec" rid="s3">Section 3</xref>, data fitting was performed on the deep time-series features extracted by TCNAE. Through operations such as nonlinear least squares, matrix computation, and vector manipulation, the estimated value of the drift coefficient <inline-formula id="ieqn-53"><mml:math id="mml-ieqn-53"><mml:mrow><mml:mover><mml:mi>&#x03BB;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> was determined to be 4.5 &#x00D7; 10<sup>&#x2212;3</sup>, and the estimated value of the diffusion coefficient squared <inline-formula id="ieqn-54"><mml:math id="mml-ieqn-54"><mml:msup><mml:mrow><mml:mover><mml:mi>&#x03C3;</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> was 1.4256 &#x00D7; 10<sup>&#x2212;4</sup>. The estimation process is illustrated in <xref ref-type="fig" rid="fig-16">Fig. 16</xref>.</p>
<fig id="fig-16">
<label>Figure 16</label>
<caption>
<title>Parameter estimation process</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-16.tif"/>
</fig>
</sec>
<sec id="s5_2">
<label>5.2</label>
<title>Decision-Making of CBM</title>
<p>The planetary gearboxes degradation progresses through two distinct phases: an initial stable operation phase where degradation remains below the preventive maintenance threshold, followed by an accelerated degradation phase where deterioration exceeds this threshold and may reach failure levels. Based on fused degradation features, we identified four characteristic stages, with degradation paths shown in <xref ref-type="fig" rid="fig-17">Fig. 17</xref> and parameters listed in <xref ref-type="table" rid="table-5">Table 5</xref>.</p>
<fig id="fig-17">
<label>Figure 17</label>
<caption>
<title>Degradation path fitted by nonlinear Wiener process</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-17.tif"/>
</fig><table-wrap id="table-5">
<label>Table 5</label>
<caption>
<title>Parameter values in each stage</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th></th>
<th><italic>a</italic> <sub><italic>i</italic></sub></th>
<th><italic>p</italic> <sub><italic>i</italic></sub></th>
<th><italic>q</italic> <sub><italic>i</italic></sub></th>
<th><italic>r</italic> <sub><italic>i</italic></sub></th>
</tr>
</thead>
<tbody>
<tr>
<td>1</td>
<td>0.1008</td>
<td>0.1011</td>
<td>0.9851</td>
<td>0.1500</td>
</tr>
<tr>
<td>2</td>
<td>0.0022</td>
<td>0.0917</td>
<td>0.0019</td>
<td>0.0100</td>
</tr>
<tr>
<td>3</td>
<td>0.1757</td>
<td>47.0933</td>
<td>0.8304</td>
<td>0.0450</td>
</tr>
<tr>
<td>4</td>
<td>0.1439</td>
<td>0.0266</td>
<td>0.0173</td>
<td>0.1400</td>
</tr>
<tr>
<td>5</td>
<td>/</td>
<td>0.0894</td>
<td>/</td>
<td>/</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Using cost parameters <italic>C</italic><sub><italic>in</italic></sub> &#x003D; 100, <italic>C</italic><sub><italic>pm</italic></sub> &#x003D; 3500, <italic>C</italic><sub><italic>cm</italic></sub> &#x003D; 4800, and <italic>C</italic><sub><italic>puni</italic></sub> &#x003D; 50, Monte Carlo simulation yielded the following results: Under static thresholds, the minimum objective function CR value was 60.9736 with optimal decision variables of 155 and 0.907. With a threshold decay rate of 0.2, the minimum CR value reduced to 55.5556, with optimal decision variables of 151 and 0.752. The corresponding 3D surface plots and contour maps at estimated values are presented in <xref ref-type="fig" rid="fig-18">Figs. 18</xref> and <xref ref-type="fig" rid="fig-19">19</xref>, respectively.</p>
<fig id="fig-18">
<label>Figure 18</label>
<caption>
<title>3D Surface Plot (<bold>a</bold>) and Contour Plot (<bold>b</bold>) of CR under static threshold</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-18.tif"/>
</fig><fig id="fig-19">
<label>Figure 19</label>
<caption>
<title>3D Surface Plot (<bold>a</bold>) and Contour Plot (<bold>b</bold>) of CR when the threshold decay rate is 0.2</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="CMC_69194-fig-19.tif"/>
</fig>
<p>A comparison of <xref ref-type="fig" rid="fig-18">Figs. 18</xref> and <xref ref-type="fig" rid="fig-19">19</xref> reveals the impact of inspection intervals on cost rates. Short intervals yield high costs due to excessive expenses from frequent inspections, leading to &#x201C;over-maintenance&#x201D;. While extending intervals initially reduces costs, excessively long intervals risk undetected degradation, causing &#x201C;under-maintenance&#x201D; and increasing failure risk and repair costs. <xref ref-type="fig" rid="fig-18">Fig. 18a</xref> exhibits a smooth surface, suitable for idealized scenario. In contrast, after introducing threshold decay, <xref ref-type="fig" rid="fig-19">Fig. 19a</xref> displays a rugged surface, accurately reflecting real-world performance degradation from aging and environmental factors. This model simulates a closed loop interaction between dynamic maintenance decisions and cost fluctuations, showing lower final cost rates with dynamic thresholds compared to fixed ones notably. In practical operations, equipment performance and degradation thresholds are not static but naturally deteriorate over time. <xref ref-type="fig" rid="fig-19">Fig. 19</xref> more accurately depicts the real-world relationship among &#x201C;maintenance strategy, dynamic thresholds, and cost control&#x201D;.</p>
<p>Based on the above research findings, to ensure the effectiveness and reliability of the degradation model for planetary gearboxes of the same type, the following conditions for model applicability can be summarized: First, the model is applicable to gearboxes with the same or similar design characteristics (such as gear size, module, and number of teeth) and material properties (such as material type, hard-ness, and toughness). These characteristics and properties directly affect the degradation mechanisms and rates. Second, the model is applicable to gearboxes operating under similar working conditions, including load levels, speed ranges, operating temperatures, and lubrication conditions. Significantly different working conditions may lead to substantial differences in degradation behavior. Third, the model is applicable to gearboxes operating under similar environmental conditions, such as humidity, temperature, and vibration environments. These factors may influence gearboxes&#x2019; degradation rates and patterns.</p>
</sec>
</sec>
<sec id="s6">
<label>6</label>
<title>Conclusion</title>
<p>This study proposes a deep learning-based CBM decision framework, establishing a complete closed loop system integrating data-driven analysis with maintenance decision support. Compared with conventional time-based maintenance approaches, the proposed strategy demonstrates superior performance in reducing maintenance costs and preventing unexpected failures. Experimental results validate the effectiveness and practicality of the proposed method, demonstrating its significance in improving maintenance timing and strategies for planetary gearboxes, enhancing system reliability and safety. Future research will explore the potential application of this method to other complex equipment systems, incorporating additional maintenance resources and environmental factors for enhanced CBM decisions. Furthermore, exploring degradation trend extraction under transfer learning frameworks and jointly optimizing maintenance scheduling with spare parts procurement presents a promising direction.</p>
</sec>
</body>
<back>
<ack>
<p>The authors sincerely appreciate the Planetary Gearbox experimental equipment and materials provided by the PHM Laboratory at the Shijiazhuang Campus of Army Engineering University of PLA. The laboratory&#x2019;s advanced testing platform and professional technical support were instrumental in data acquisition and analysis, for which we are deeply grateful.</p>
</ack>
<sec>
<title>Funding Statement</title>
<p>This research was funded by scientific research projects under Grant JY2024B011.</p>
</sec>
<sec>
<title>Author Contributions</title>
<p>The authors confirm contribution to the paper as follows: Conceptualization, Bo Zhu, Zhonghua Cheng; methodology, Xianbiao Zhan; software, Enzhi Dong; validation, Kexin Jiang; formal analysis, Bo Zhu; investigation, Enzhi Dong; resources, Xianbiao Zhan, Rongcai Wang; data curation, Xianbiao Zhan; writing&#x2014;original draft preparation, Bo Zhu; writing&#x2014;review and editing, Enzhi Dong, Kexin Jiang, Rongcai Wang, Zhonghua Cheng; visualization, Kexin Jiang; supervision, Zhonghua Cheng; project administration, Zhonghua Cheng. All authors reviewed the results and approved the final version of the manuscript.</p>
</sec>
<sec sec-type="data-availability">
<title>Availability of Data and Materials</title>
<p>Data available on request from the authors. The data that support the findings of this study are available from the corresponding author, Zhonghua Cheng, upon reasonable request.</p>
</sec>
<sec>
<title>Ethics Approval</title>
<p>Not applicable.</p>
</sec>
<sec sec-type="COI-statement">
<title>Conflicts of Interest</title>
<p>The authors declare no conflicts of interest to report regarding the present study.</p>
</sec>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1"><label>[1]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname> <given-names>D</given-names></string-name>, <string-name><surname>Zhao</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Zhao</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>A dynamic-model-based fault diagnosis method for a wind turbine planetary gearbox using a deep learning network</article-title>. <source>Prot Control Mod Power Syst</source>. <year>2022</year>;<volume>7</volume>:<fpage>22</fpage>. doi:<pub-id pub-id-type="doi">10.1186/s41601-022-00244-z</pub-id>.</mixed-citation></ref>
<ref id="ref-2"><label>[2]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhang</surname> <given-names>L</given-names></string-name>, <string-name><surname>Hu</surname> <given-names>N</given-names></string-name></person-group>. <article-title>Fault diagnosis of sun gear based on continuous vibration separation and minimum entropy deconvolution</article-title>. <source>Measurement</source>. <year>2019</year>;<volume>141</volume>:<fpage>332</fpage>&#x2013;<lpage>44</lpage>. doi:<pub-id pub-id-type="doi">10.1016/j.measurement.2019.04.049</pub-id>.</mixed-citation></ref>
<ref id="ref-3"><label>[3]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Mourtzis</surname> <given-names>D</given-names></string-name></person-group>. <article-title>Special issue &#x201C;towards digitalized manufacturing 4.0</article-title>. <source>Int J Comput Integr Manuf</source>. <year>2023</year>;<volume>36</volume>(<issue>12</issue>):<fpage>1735</fpage>&#x2013;<lpage>8</lpage>. doi:<pub-id pub-id-type="doi">10.1080/0951192x.2023.2290394</pub-id>.</mixed-citation></ref>
<ref id="ref-4"><label>[4]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tong</surname> <given-names>S</given-names></string-name>, <string-name><surname>Li</surname> <given-names>J</given-names></string-name>, <string-name><surname>Cong</surname> <given-names>F</given-names></string-name>, <string-name><surname>Fu</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Tong</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>Vibration separation methodology compensated by time-varying transfer function for fault diagnosis of non-hunting tooth planetary gearbox</article-title>. <source>Sensors</source>. <year>2022</year>;<volume>22</volume>(<issue>2</issue>):<fpage>557</fpage>. doi:<pub-id pub-id-type="doi">10.3390/s22020557</pub-id>; <pub-id pub-id-type="pmid">35062518</pub-id></mixed-citation></ref>
<ref id="ref-5"><label>[5]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Dao</surname> <given-names>CD</given-names></string-name>, <string-name><surname>Kazemtabrizi</surname> <given-names>B</given-names></string-name>, <string-name><surname>Crabtree</surname> <given-names>CJ</given-names></string-name>, <string-name><surname>Tavner</surname> <given-names>PJ</given-names></string-name></person-group>. <article-title>Integrated condition-based maintenance modelling and optimisation for offshore wind turbines</article-title>. <source>Wind Energy</source>. <year>2021</year>;<volume>24</volume>(<issue>11</issue>):<fpage>1180</fpage>&#x2013;<lpage>98</lpage>. doi:<pub-id pub-id-type="doi">10.1002/we.2625</pub-id>.</mixed-citation></ref>
<ref id="ref-6"><label>[6]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Luo</surname> <given-names>J</given-names></string-name>, <string-name><surname>Luo</surname> <given-names>X</given-names></string-name>, <string-name><surname>Ma</surname> <given-names>X</given-names></string-name>, <string-name><surname>Zan</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Bai</surname> <given-names>X</given-names></string-name></person-group>. <article-title>An integrated condition-based opportunistic maintenance framework for offshore wind farms</article-title>. <source>Reliab Eng Syst Saf</source>. <year>2025</year>;<volume>256</volume>:<fpage>110701</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ress.2024.110701</pub-id>.</mixed-citation></ref>
<ref id="ref-7"><label>[7]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname> <given-names>S</given-names></string-name>, <string-name><surname>Zhao</surname> <given-names>M</given-names></string-name>, <string-name><surname>Ou</surname> <given-names>S</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>D</given-names></string-name>, <string-name><surname>Wei</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>A periodic anomaly detection framework based on matrix profile for condition monitoring of planetary gearboxes</article-title>. <source>Measurement</source>. <year>2023</year>;<volume>218</volume>:<fpage>113243</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.measurement.2023.113243</pub-id>.</mixed-citation></ref>
<ref id="ref-8"><label>[8]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhang</surname> <given-names>M</given-names></string-name>, <string-name><surname>Xie</surname> <given-names>R</given-names></string-name>, <string-name><surname>Kang</surname> <given-names>T</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>J</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Feng</surname> <given-names>X</given-names></string-name>, <etal>et al</etal></person-group>. <article-title>A novel co-modulation and hybrid resolution strategy (CHRS) for fault diagnosis of planetary gearboxes</article-title>. <source>Mech Syst Signal Process</source>. <year>2025</year>;<volume>230</volume>:<fpage>112573</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ymssp.2025.112573</pub-id>.</mixed-citation></ref>
<ref id="ref-9"><label>[9]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wei</surname> <given-names>J</given-names></string-name>, <string-name><surname>Dong</surname> <given-names>G</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>Remaining useful life prediction and state of health diagnosis for lithium-ion batteries using particle filter and support vector regression</article-title>. <source>IEEE Trans Ind Electron</source>. <year>2018</year>;<volume>65</volume>(<issue>7</issue>):<fpage>5634</fpage>&#x2013;<lpage>43</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TIE.2017.2782224</pub-id>.</mixed-citation></ref>
<ref id="ref-10"><label>[10]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Chen</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Cao</surname> <given-names>S</given-names></string-name>, <string-name><surname>Mao</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>Remaining useful life estimation of aircraft engines using a modified similarity and supporting vector machine (SVM) approach</article-title>. <source>Energies</source>. <year>2018</year>;<volume>11</volume>(<issue>1</issue>):<fpage>28</fpage>. doi:<pub-id pub-id-type="doi">10.3390/en11010028</pub-id>.</mixed-citation></ref>
<ref id="ref-11"><label>[11]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Pei</surname> <given-names>H</given-names></string-name>, <string-name><surname>Hu</surname> <given-names>C</given-names></string-name>, <string-name><surname>Si</surname> <given-names>X</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>J</given-names></string-name></person-group>. <article-title>Remaining useful life prediction methods for equipment based on machine learning: a review</article-title>. <source>J Mech Eng</source>. <year>2019</year>;<volume>55</volume>(<issue>8</issue>):<fpage>1</fpage>&#x2013;<lpage>13</lpage>. doi:<pub-id pub-id-type="doi">10.3901/JME.2019.08.001</pub-id>.</mixed-citation></ref>
<ref id="ref-12"><label>[12]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhang</surname> <given-names>W</given-names></string-name>, <string-name><surname>Yang</surname> <given-names>D</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>H</given-names></string-name></person-group>. <article-title>Data-driven methods for predictive maintenance of industrial equipment: a survey</article-title>. <source>IEEE Syst J</source>. <year>2019</year>;<volume>13</volume>(<issue>3</issue>):<fpage>2213</fpage>&#x2013;<lpage>27</lpage>. doi:<pub-id pub-id-type="doi">10.1109/JSYST.2019.2905565</pub-id>.</mixed-citation></ref>
<ref id="ref-13"><label>[13]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Avci</surname> <given-names>O</given-names></string-name>, <string-name><surname>Abdeljaber</surname> <given-names>O</given-names></string-name>, <string-name><surname>Kiranyaz</surname> <given-names>S</given-names></string-name>, <string-name><surname>Hussein</surname> <given-names>M</given-names></string-name>, <string-name><surname>Gabbouj</surname> <given-names>M</given-names></string-name>, <string-name><surname>Inman</surname> <given-names>DJ</given-names></string-name></person-group>. <article-title>A review of vibration-based damage detection in civil structures: from traditional methods to machine learning and deep learning applications</article-title>. <source>Mech Syst Signal Process</source>. <year>2021</year>;<volume>147</volume>:<fpage>107077</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ymssp.2020.107077</pub-id>.</mixed-citation></ref>
<ref id="ref-14"><label>[14]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhao</surname> <given-names>D</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>T</given-names></string-name>, <string-name><surname>Chu</surname> <given-names>F</given-names></string-name></person-group>. <article-title>Deep convolutional neural network based planet bearing fault classification</article-title>. <source>Comput Ind</source>. <year>2019</year>;<volume>107</volume>:<fpage>59</fpage>&#x2013;<lpage>66</lpage>. doi:<pub-id pub-id-type="doi">10.1016/j.compind.2019.02.001</pub-id>.</mixed-citation></ref>
<ref id="ref-15"><label>[15]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhao</surname> <given-names>M</given-names></string-name>, <string-name><surname>Kang</surname> <given-names>M</given-names></string-name>, <string-name><surname>Tang</surname> <given-names>B</given-names></string-name>, <string-name><surname>Pecht</surname> <given-names>M</given-names></string-name></person-group>. <article-title>Deep residual networks with dynamically weighted wavelet coefficients for fault diagnosis of planetary gearboxes</article-title>. <source>IEEE Trans Ind Electron</source>. <year>2018</year>;<volume>65</volume>(<issue>5</issue>):<fpage>4290</fpage>&#x2013;<lpage>300</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TIE.2017.2762639</pub-id>.</mixed-citation></ref>
<ref id="ref-16"><label>[16]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>J</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>An intelligent diagnosis scheme based on generative adversarial learning deep neural networks and its application to planetary gearbox fault pattern recognition</article-title>. <source>Neurocomputing</source>. <year>2018</year>;<volume>310</volume>:<fpage>213</fpage>&#x2013;<lpage>22</lpage>. doi:<pub-id pub-id-type="doi">10.1016/j.neucom.2018.05.024</pub-id>.</mixed-citation></ref>
<ref id="ref-17"><label>[17]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ali</surname> <given-names>A</given-names></string-name>, <string-name><surname>Abdelhadi</surname> <given-names>A</given-names></string-name></person-group>. <article-title>Condition-based monitoring and maintenance: state of the art review</article-title>. <source>Appl Sci</source>. <year>2022</year>;<volume>12</volume>(<issue>2</issue>):<fpage>688</fpage>. doi:<pub-id pub-id-type="doi">10.3390/app12020688</pub-id>.</mixed-citation></ref>
<ref id="ref-18"><label>[18]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Alaswad</surname> <given-names>S</given-names></string-name>, <string-name><surname>Xiang</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>A review on condition-based maintenance optimization models for stochastically deteriorating system</article-title>. <source>Reliab Eng Syst Saf</source>. <year>2017</year>;<volume>157</volume>:<fpage>54</fpage>&#x2013;<lpage>63</lpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ress.2016.08.009</pub-id>.</mixed-citation></ref>
<ref id="ref-19"><label>[19]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Michalak</surname> <given-names>A</given-names></string-name>, <string-name><surname>Wodecki</surname> <given-names>J</given-names></string-name>, <string-name><surname>Drozda</surname> <given-names>M</given-names></string-name>, <string-name><surname>Wy&#x0142;oma&#x0144;ska</surname> <given-names>A</given-names></string-name>, <string-name><surname>Zimroz</surname> <given-names>R</given-names></string-name></person-group>. <article-title>Model of the vibration signal of the vibrating sieving screen suspension for condition monitoring purposes</article-title>. <source>Sensors</source>. <year>2020</year>;<volume>21</volume>(<issue>1</issue>):<fpage>213</fpage>. doi:<pub-id pub-id-type="doi">10.3390/s21010213</pub-id>; <pub-id pub-id-type="pmid">33396259</pub-id></mixed-citation></ref>
<ref id="ref-20"><label>[20]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Chen</surname> <given-names>M</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>H</given-names></string-name>, <string-name><surname>Cheng</surname> <given-names>J</given-names></string-name>, <string-name><surname>Zhao</surname> <given-names>L</given-names></string-name>, <string-name><surname>Yang</surname> <given-names>H</given-names></string-name>, <string-name><surname>Liu</surname> <given-names>Q</given-names></string-name>, <etal>et al</etal></person-group>. <article-title>Progress in detection and suppression techniques for processing-induced subsurface defects of fused silica optical elements</article-title>. <source>J Mech Eng</source>. <year>2021</year>;<volume>57</volume>(<issue>20</issue>):<fpage>1</fpage>&#x2013;<lpage>19</lpage>. doi:<pub-id pub-id-type="doi">10.3901/JME.2021.20.001</pub-id>.</mixed-citation></ref>
<ref id="ref-21"><label>[21]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liu</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>Smart sensors for structural health monitoring and nondestructive evaluation</article-title>. <source>Sensors</source>. <year>2024</year>;<volume>24</volume>(<issue>2</issue>):<fpage>603</fpage>. doi:<pub-id pub-id-type="doi">10.3390/s24020603</pub-id>; <pub-id pub-id-type="pmid">38257695</pub-id></mixed-citation></ref>
<ref id="ref-22"><label>[22]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Naga Srinivasa Rao</surname> <given-names>P</given-names></string-name>, <string-name><surname>Achutha Naikan</surname> <given-names>VN</given-names></string-name></person-group>. <article-title>An algorithm for simultaneous optimization of parameters of condition-based preventive maintenance</article-title>. <source>Struct Health Monit</source>. <year>2009</year>;<volume>8</volume>(<issue>1</issue>):<fpage>83</fpage>&#x2013;<lpage>94</lpage>. doi:<pub-id pub-id-type="doi">10.1177/1475921708094793</pub-id>.</mixed-citation></ref>
<ref id="ref-23"><label>[23]</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name><surname>Cheng</surname> <given-names>J</given-names></string-name></person-group>. <article-title>Research on full life cycle maintenance strategy of fatigue damage structure based on reinforcement learning [master&#x2019;s thesis]. Wuhan, China: Huazhong University of Science and Technology</article-title>; <year>2021</year>. doi:<pub-id pub-id-type="doi">10.27157/d.cnki.ghzku.2021.003061</pub-id>.</mixed-citation></ref>
<ref id="ref-24"><label>[24]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Kang</surname> <given-names>J</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Guedes Soares</surname> <given-names>C</given-names></string-name></person-group>. <article-title>Condition-based maintenance for offshore wind turbines based on support vector machine</article-title>. <source>Energies</source>. <year>2020</year>;<volume>13</volume>(<issue>14</issue>):<fpage>3518</fpage>. doi:<pub-id pub-id-type="doi">10.3390/en13143518</pub-id>.</mixed-citation></ref>
<ref id="ref-25"><label>[25]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhao</surname> <given-names>H</given-names></string-name>, <string-name><surname>Xu</surname> <given-names>F</given-names></string-name>, <string-name><surname>Liang</surname> <given-names>B</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>J</given-names></string-name>, <string-name><surname>Song</surname> <given-names>P</given-names></string-name></person-group>. <article-title>A condition-based opportunistic maintenance strategy for multi-component system</article-title>. <source>Struct Health Monit</source>. <year>2019</year>;<volume>18</volume>(<issue>1</issue>):<fpage>270</fpage>&#x2013;<lpage>83</lpage>. doi:<pub-id pub-id-type="doi">10.1177/1475921717751871</pub-id>.</mixed-citation></ref>
<ref id="ref-26"><label>[26]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhan</surname> <given-names>X</given-names></string-name>, <string-name><surname>Liu</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Yan</surname> <given-names>H</given-names></string-name>, <string-name><surname>Wu</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Guo</surname> <given-names>C</given-names></string-name>, <string-name><surname>Jia</surname> <given-names>X</given-names></string-name></person-group>. <article-title>A novel method of health indicator construction and remaining useful life prediction based on deep learning</article-title>. <source>Eksploat I Niezawodn</source>. <year>2023</year>;<volume>25</volume>(<issue>4</issue>):<fpage>171374</fpage>. doi:<pub-id pub-id-type="doi">10.17531/ein/171374</pub-id>.</mixed-citation></ref>
<ref id="ref-27"><label>[27]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Shen</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Shi</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Zhen</surname> <given-names>D</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>H</given-names></string-name>, <string-name><surname>Gu</surname> <given-names>F</given-names></string-name></person-group>. <article-title>Fault diagnosis of planetary gearbox based on adaptive order bispectrum slice and fault characteristics energy ratio analysis</article-title>. <source>Sensors</source>. <year>2020</year>;<volume>20</volume>(<issue>8</issue>):<fpage>2433</fpage>. doi:<pub-id pub-id-type="doi">10.3390/s20082433</pub-id>; <pub-id pub-id-type="pmid">32344737</pub-id></mixed-citation></ref>
<ref id="ref-28"><label>[28]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Huang</surname> <given-names>LT</given-names></string-name>, <string-name><surname>Chang</surname> <given-names>JY</given-names></string-name></person-group>. <article-title>Vibration characterization and fault diagnosis of a planetary gearbox with a wireless embedded sensor</article-title>. <source>Appl Sci</source>. <year>2023</year>;<volume>13</volume>(<issue>2</issue>):<fpage>729</fpage>. doi:<pub-id pub-id-type="doi">10.3390/app13020729</pub-id>.</mixed-citation></ref>
<ref id="ref-29"><label>[29]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Al-Bedhany</surname> <given-names>JH</given-names></string-name>, <string-name><surname>Mankhi</surname> <given-names>TA</given-names></string-name>, <string-name><surname>Legutko</surname> <given-names>S</given-names></string-name></person-group>. <article-title>A surface study of failed planetary wind turbine gearbox Bearings to investigate the causes of the Bearing premature failure issue</article-title>. <source>Heliyon</source>. <year>2024</year>;<volume>10</volume>(<issue>4</issue>):<fpage>e25860</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.heliyon.2024.e25860</pub-id>; <pub-id pub-id-type="pmid">38390070</pub-id></mixed-citation></ref>
<ref id="ref-30"><label>[30]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Huang</surname> <given-names>X</given-names></string-name>, <string-name><surname>Li</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Chai</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>Intelligent fault diagnosis method of wind turbines planetary gearboxes based on a multi-scale dense fusion network</article-title>. <source>Front Energy Res</source>. <year>2021</year>;<volume>9</volume>:<fpage>747622</fpage>. doi:<pub-id pub-id-type="doi">10.3389/fenrg.2021.747622</pub-id>.</mixed-citation></ref>
<ref id="ref-31"><label>[31]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wen</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Wu</surname> <given-names>J</given-names></string-name>, <string-name><surname>Das</surname> <given-names>D</given-names></string-name>, <string-name><surname>Tseng</surname> <given-names>TL</given-names></string-name></person-group>. <article-title>Degradation modeling and RUL prediction using Wiener process subject to multiple change points and unit heterogeneity</article-title>. <source>Reliab Eng Syst Saf</source>. <year>2018</year>;<volume>176</volume>:<fpage>113</fpage>&#x2013;<lpage>24</lpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ress.2018.04.005</pub-id>.</mixed-citation></ref>
<ref id="ref-32"><label>[32]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname> <given-names>J</given-names></string-name>, <string-name><surname>Fan</surname> <given-names>J</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Pang</surname> <given-names>X</given-names></string-name>, <string-name><surname>Qiu</surname> <given-names>M</given-names></string-name></person-group>. <article-title>Enhanced RUL predictions of rolling bearings using a nonlinear Wiener model with an extended incremental Kalman filter</article-title>. <source>Meas Sci Technol</source>. <year>2025</year>;<volume>36</volume>(<issue>1</issue>):<fpage>016196</fpage>. doi:<pub-id pub-id-type="doi">10.1088/1361-6501/ad9bd3</pub-id>.</mixed-citation></ref>
<ref id="ref-33"><label>[33]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Sun</surname> <given-names>L</given-names></string-name>, <string-name><surname>Liang</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Yan</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>Remaining useful life (RUL) prediction based on the bivariant two-phase nonlinear Wiener degradation process</article-title>. <source>Entropy</source>. <year>2025</year>;<volume>27</volume>(<issue>4</issue>):<fpage>349</fpage>. doi:<pub-id pub-id-type="doi">10.3390/e27040349</pub-id>; <pub-id pub-id-type="pmid">40282584</pub-id></mixed-citation></ref>
<ref id="ref-34"><label>[34]</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name><surname>Huang</surname> <given-names>Q</given-names></string-name></person-group>. <article-title>Research on health state assessment and remaining useful life prediction method of rolling bearing [master&#x2019;s thesis]. Suzhou, China: Soochow University</article-title>; <year>2022</year>. doi:<pub-id pub-id-type="doi">10.27351/d.cnki.gszhu.2022.002373</pub-id>.</mixed-citation></ref>
<ref id="ref-35"><label>[35]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Huang</surname> <given-names>X</given-names></string-name>, <string-name><surname>Qi</surname> <given-names>G</given-names></string-name>, <string-name><surname>Mazur</surname> <given-names>N</given-names></string-name>, <string-name><surname>Chai</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>Deep residual networks-based intelligent fault diagnosis method of planetary gearboxes in cloud environments</article-title>. <source>Simul Model Pract Theory</source>. <year>2022</year>;<volume>116</volume>:<fpage>102469</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.simpat.2021.102469</pub-id>.</mixed-citation></ref>
<ref id="ref-36"><label>[36]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname> <given-names>X</given-names></string-name>, <string-name><surname>Jiang</surname> <given-names>H</given-names></string-name>, <string-name><surname>Mu</surname> <given-names>M</given-names></string-name>, <string-name><surname>Dong</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>A dynamic collaborative adversarial domain adaptation network for unsupervised rotating machinery fault diagnosis</article-title>. <source>Reliab Eng Syst Saf</source>. <year>2025</year>;<volume>255</volume>:<fpage>110662</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ress.2024.110662</pub-id>.</mixed-citation></ref>
<ref id="ref-37"><label>[37]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Dong</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Jiang</surname> <given-names>H</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>X</given-names></string-name>, <string-name><surname>Mu</surname> <given-names>M</given-names></string-name>, <string-name><surname>Jiang</surname> <given-names>W</given-names></string-name></person-group>. <article-title>An interpretable multiscale lifting wavelet contrast network for planetary gearbox fault diagnosis with small samples</article-title>. <source>Reliab Eng Syst Saf</source>. <year>2024</year>;<volume>251</volume>:<fpage>110404</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ress.2024.110404</pub-id>.</mixed-citation></ref>
<ref id="ref-38"><label>[38]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Mu</surname> <given-names>M</given-names></string-name>, <string-name><surname>Jiang</surname> <given-names>H</given-names></string-name>, <string-name><surname>Jiang</surname> <given-names>W</given-names></string-name>, <string-name><surname>Dong</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Wu</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>An enhanced meta-learning network with sensitivity penalty for cross-domain few-shot fault diagnosis</article-title>. <source>Meas Sci Technol</source>. <year>2024</year>;<volume>35</volume>(<issue>9</issue>):<fpage>095106</fpage>. doi:<pub-id pub-id-type="doi">10.1088/1361-6501/ad5039</pub-id>.</mixed-citation></ref>
<ref id="ref-39"><label>[39]</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name><surname>Lea</surname> <given-names>C</given-names></string-name>, <string-name><surname>Vidal</surname> <given-names>R</given-names></string-name>, <string-name><surname>Reiter</surname> <given-names>A</given-names></string-name>, <string-name><surname>Hager</surname> <given-names>GD</given-names></string-name></person-group>. <article-title>Temporal convolutional networks: a unified approach to action segmentation</article-title>. <comment>arXiv:1608.08242. 2016</comment>. doi:<pub-id pub-id-type="doi">10.48550/arXiv.1608.08242</pub-id>.</mixed-citation></ref>
<ref id="ref-40"><label>[40]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Filho</surname> <given-names>JBOS</given-names></string-name>, <string-name><surname>Diniz</surname> <given-names>PSR</given-names></string-name></person-group>. <article-title>A sparse fixed-point online KPCA extraction algorithm</article-title>. <source>IEEE Trans Signal Process</source>. <year>2024</year>;<volume>72</volume>:<fpage>4604</fpage>&#x2013;<lpage>17</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TSP.2024.3446512</pub-id>.</mixed-citation></ref>
<ref id="ref-41"><label>[41]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Hofmann</surname> <given-names>T</given-names></string-name>, <string-name><surname>Sch&#x00F6;lkopf</surname> <given-names>B</given-names></string-name>, <string-name><surname>Smola</surname> <given-names>AJ</given-names></string-name></person-group>. <article-title>Kernel methods in machine learning</article-title>. <source>Ann Statist</source>. <year>2008</year>;<volume>36</volume>(<issue>3</issue>):<fpage>1171</fpage>&#x2013;<lpage>220</lpage>. doi:<pub-id pub-id-type="doi">10.1214/009053607000000677</pub-id>.</mixed-citation></ref>
<ref id="ref-42"><label>[42]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liu</surname> <given-names>S</given-names></string-name>, <string-name><surname>Gao</surname> <given-names>L</given-names></string-name>, <string-name><surname>Du</surname> <given-names>J</given-names></string-name>, <string-name><surname>Liu</surname> <given-names>C</given-names></string-name></person-group>. <article-title>Optimal feature selection and life prediction of rolling bearing based on adaptive order</article-title>. <source>Ship Sci Technol</source>. <year>2019</year>;<volume>41</volume>(<issue>11</issue>):<fpage>71</fpage>&#x2013;<lpage>6</lpage>. doi:<pub-id pub-id-type="doi">10.3404/j.issn.1672-7649.2019.11.014</pub-id>.</mixed-citation></ref>
<ref id="ref-43"><label>[43]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Si</surname> <given-names>XS</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>W</given-names></string-name>, <string-name><surname>Hu</surname> <given-names>CH</given-names></string-name>, <string-name><surname>Zhou</surname> <given-names>DH</given-names></string-name>, <string-name><surname>Pecht</surname> <given-names>MG</given-names></string-name></person-group>. <article-title>Remaining useful life estimation based on a nonlinear diffusion degradation process</article-title>. <source>IEEE Trans Reliab</source>. <year>2012</year>;<volume>61</volume>(<issue>1</issue>):<fpage>50</fpage>&#x2013;<lpage>67</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TR.2011.2182221</pub-id>.</mixed-citation></ref>
<ref id="ref-44"><label>[44]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>He</surname> <given-names>R</given-names></string-name>, <string-name><surname>Tian</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Zuo</surname> <given-names>M</given-names></string-name>, <string-name><surname>Guo</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>Condition-based maintenance optimization for multi-component systems considering prognostic information and degraded working efficiency</article-title>. <source>Reliab Eng Syst Saf</source>. <year>2023</year>;<volume>234</volume>:<fpage>109167</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.ress.2023.109167</pub-id>.</mixed-citation></ref>
<ref id="ref-45"><label>[45]</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name><surname>Yao</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>Research on reliability modeling and condition-based maintenance strategy for two-stage degradation system [master&#x2019;s thesis]. Tianjin, China: Tianjin University</article-title>; <year>2020</year>. doi:<pub-id pub-id-type="doi">10.27356/d.cnki.gtjdu.2020.001899</pub-id>.</mixed-citation></ref>
<ref id="ref-46"><label>[46]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Oh</surname> <given-names>SY</given-names></string-name>, <string-name><surname>Joung</surname> <given-names>C</given-names></string-name>, <string-name><surname>Lee</surname> <given-names>S</given-names></string-name>, <string-name><surname>Shim</surname> <given-names>YB</given-names></string-name>, <string-name><surname>Lee</surname> <given-names>D</given-names></string-name>, <string-name><surname>Cho</surname> <given-names>GE</given-names></string-name>, <etal>et al</etal></person-group>. <article-title>Condition-based maintenance of wind turbine structures: a state-of-the-art review</article-title>. <source>Renew Sustain Energy Rev</source>. <year>2024</year>;<volume>204</volume>:<fpage>114799</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.rser.2024.114799</pub-id>.</mixed-citation></ref>
<ref id="ref-47"><label>[47]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Dong</surname> <given-names>E</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Zhan</surname> <given-names>X</given-names></string-name>, <string-name><surname>Bai</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Cheng</surname> <given-names>Z</given-names></string-name></person-group>. <article-title>A novel dynamic predictive maintenance framework for gearboxes utilizing nonlinear Wiener process</article-title>. <source>Meas Sci Technol</source>. <year>2024</year>;<volume>35</volume>(<issue>12</issue>):<fpage>126210</fpage>. doi:<pub-id pub-id-type="doi">10.1088/1361-6501/ad762e</pub-id>.</mixed-citation></ref>
<ref id="ref-48"><label>[48]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tang</surname> <given-names>H</given-names></string-name>, <string-name><surname>Tang</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Su</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Feng</surname> <given-names>W</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>B</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>P</given-names></string-name>, <etal>et al</etal></person-group>. <article-title>Feature extraction of multi-sensors for early bearing fault diagnosis using deep learning based on minimum unscented Kalman filter</article-title>. <source>Eng Appl Artif Intell</source>. <year>2024</year>;<volume>127</volume>:<fpage>107138</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.engappai.2023.107138</pub-id>.</mixed-citation></ref>
</ref-list>
</back></article>