Pub Date : 2026-01-01Epub Date: 2025-12-23DOI: 10.1016/j.probengmech.2025.103881
Xinzhou Qiao , Yan Li , Isaac Elishakoff , Naigang Hu
The currently available multidimensional parallelepiped models in non-probabilistic convex treatment of uncertainty are applicable for the multi-source uncertainty problem with the coexistence of dependent and independent uncertain parameters. However, these may encounter obstacles either in possessing the minimum volume or in enclosing all experimental data, or in both. In this paper, a novel multidimensional parallelepiped model, which is defined as the intersection of a finite number of halfspaces, is proposed to bound the uncertainty domain to address the above issue. Based on the proposed model, two uncertainty propagation analysis methods, namely the Monte Carlo simulation method and the sub-multidimensional parallelepiped analysis method, are developed to predict the structural response interval. Three numerical examples are provided to demonstrate the superiority of the proposed model over the existing ones and to illustrate the effectiveness and validity of the proposed methods.
{"title":"A novel multidimensional parallelepiped model for structural uncertainty quantification and propagation analysis","authors":"Xinzhou Qiao , Yan Li , Isaac Elishakoff , Naigang Hu","doi":"10.1016/j.probengmech.2025.103881","DOIUrl":"10.1016/j.probengmech.2025.103881","url":null,"abstract":"<div><div>The currently available multidimensional parallelepiped models in non-probabilistic convex treatment of uncertainty are applicable for the multi-source uncertainty problem with the coexistence of dependent and independent uncertain parameters. However, these may encounter obstacles either in possessing the minimum volume or in enclosing all experimental data, or in both. In this paper, a novel multidimensional parallelepiped model, which is defined as the intersection of a finite number of halfspaces, is proposed to bound the uncertainty domain to address the above issue. Based on the proposed model, two uncertainty propagation analysis methods, namely the Monte Carlo simulation method and the sub-multidimensional parallelepiped analysis method, are developed to predict the structural response interval. Three numerical examples are provided to demonstrate the superiority of the proposed model over the existing ones and to illustrate the effectiveness and validity of the proposed methods.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103881"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145840780","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-12-23DOI: 10.1016/j.probengmech.2025.103884
Zonghui Wu , Jian He , Chenyang Wang , Xiaodan Sun , Di Yao
The fundamental purpose of structural reliability analysis is defined as the quantitative measurement of structural failure possibilities. The surrogate model method is currently regarded as the most widely used reliability evaluation method, but it has problems such as low fitting accuracy, high computational cost, low convergence efficiency and high parameter sensitivity when dealing with small probability events. Although there are some methods to accelerate the analysis, such as Adaptive Kriging combined with Monte Carlo Simulation (AK-MCS), the construction of the model still requires a large number of samples, resulting in a very large amount of calculation of the surrogate model. Therefore, this study combines the adaptive Kriging model with advanced subset simulation (AK-ASS) to solve these problems. In this paper, through the verification of mathematical examples and engineering examples, it is proved that this method reduces the analysis time required to deal with the problem of small probability failure, and overcomes some limitations of subset simulation. Furthermore, it has the potential to be used in combination with new efficient learning functions in the future.
结构可靠度分析的根本目的是对结构破坏可能性进行定量测量。代理模型法是目前应用最广泛的可靠性评估方法,但在处理小概率事件时存在拟合精度低、计算成本高、收敛效率低和参数灵敏度高等问题。虽然有一些加速分析的方法,如Adaptive Kriging结合Monte Carlo Simulation (AK-MCS),但模型的构建仍然需要大量的样本,导致代理模型的计算量非常大。因此,本研究将自适应Kriging模型与先进子集仿真(AK-ASS)相结合来解决这些问题。本文通过数学实例和工程实例的验证,证明了该方法减少了处理小概率故障问题所需的分析时间,克服了子集仿真的一些局限性。此外,它在未来有可能与新的高效学习函数结合使用。
{"title":"AK-ASS: An improvement of the Kriging model for dealing with small failure probability problems","authors":"Zonghui Wu , Jian He , Chenyang Wang , Xiaodan Sun , Di Yao","doi":"10.1016/j.probengmech.2025.103884","DOIUrl":"10.1016/j.probengmech.2025.103884","url":null,"abstract":"<div><div>The fundamental purpose of structural reliability analysis is defined as the quantitative measurement of structural failure possibilities. The surrogate model method is currently regarded as the most widely used reliability evaluation method, but it has problems such as low fitting accuracy, high computational cost, low convergence efficiency and high parameter sensitivity when dealing with small probability events. Although there are some methods to accelerate the analysis, such as Adaptive Kriging combined with Monte Carlo Simulation (AK-MCS), the construction of the model still requires a large number of samples, resulting in a very large amount of calculation of the surrogate model. Therefore, this study combines the adaptive Kriging model with advanced subset simulation (AK-ASS) to solve these problems. In this paper, through the verification of mathematical examples and engineering examples, it is proved that this method reduces the analysis time required to deal with the problem of small probability failure, and overcomes some limitations of subset simulation. Furthermore, it has the potential to be used in combination with new efficient learning functions in the future.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103884"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884335","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2026-01-30DOI: 10.1016/j.probengmech.2026.103895
Hengchao Li , Zhenzhou Lu , Yuhua Yan , Yixin Lu
When distribution parameters are uncertain, it is necessary to employ a continuous estimation method for posterior failure probabilities (PFPs) to efficiently track changes in structural reliability as new observations become available. This study proposes an efficient method for continuously estimating PFPs, thereby addressing the lack of efficient and robust methods. In the proposed method, the integrand of PFP can be equivalently expressed as an explicit updating factor operation related to the observations and the estimations of the conditional failure probabilities based on a single reliability analysis in an augmented space. Thus, repeated reliability analyses are prevented when new observations are gradually incorporated. An interpolation-based method is designed to estimate all necessary conditional failure probabilities. This method first leverages the sifting property of the Dirac delta function to derive the transformed expression of the conditional probability density for the random input vector. A normal density approximating the Dirac delta function is then adopted as the interpolation weight function to solve the integral. This enables sifting the same set of sample information in the augmented space to calculate all conditional failure probabilities. An adaptive Kriging model of the performance function is introduced to further enhance the efficiency of the reliability analysis. Examples demonstrate that compared with existing advanced methods, the proposed method notably improves the efficiency of continuous PFP estimation while maintaining high accuracy.
{"title":"An efficient method for continuously estimating posterior failure probabilities based on a single reliability analysis","authors":"Hengchao Li , Zhenzhou Lu , Yuhua Yan , Yixin Lu","doi":"10.1016/j.probengmech.2026.103895","DOIUrl":"10.1016/j.probengmech.2026.103895","url":null,"abstract":"<div><div>When distribution parameters are uncertain, it is necessary to employ a continuous estimation method for posterior failure probabilities (PFPs) to efficiently track changes in structural reliability as new observations become available. This study proposes an efficient method for continuously estimating PFPs, thereby addressing the lack of efficient and robust methods. In the proposed method, the integrand of PFP can be equivalently expressed as an explicit updating factor operation related to the observations and the estimations of the conditional failure probabilities based on a single reliability analysis in an augmented space. Thus, repeated reliability analyses are prevented when new observations are gradually incorporated. An interpolation-based method is designed to estimate all necessary conditional failure probabilities. This method first leverages the sifting property of the Dirac delta function to derive the transformed expression of the conditional probability density for the random input vector. A normal density approximating the Dirac delta function is then adopted as the interpolation weight function to solve the integral. This enables sifting the same set of sample information in the augmented space to calculate all conditional failure probabilities. An adaptive Kriging model of the performance function is introduced to further enhance the efficiency of the reliability analysis. Examples demonstrate that compared with existing advanced methods, the proposed method notably improves the efficiency of continuous PFP estimation while maintaining high accuracy.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103895"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395575","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The assessment of the structural behavior of masonry towers often involves developing and identifying computational models to be used to perform static and/or time-history nonlinear analyses. Such models are frequently developed assuming isotropic masonry behavior, with model identification carried out - based on available experimental data - through deterministic tuning of a limited set of representative parameters. However, masonry exhibits complex structural textures that often deviate significantly from the commonly made assumption of isotropic behaviour. This paper investigates the effects of this assumption by focusing on the masonry shear modulus. A two-dimensional rigid body–spring model is adopted as computational approach, as its fully discrete formulation allows overcoming the limitations of the Cauchy continuum while maintaining a reduced computational cost. A Bayesian updating framework based on dynamic experimental data is employed to account for multiple sources of uncertainty, including parameter uncertainty, model inadequacy, and observation errors. Three masonry towers with different slenderness ratios are considered as representative case studies. The adopted Bayesian model updating approach allows for the estimation of their shear modulus while accounting for uncertainties, and the results show that the common assumption of isotropic behaviour in masonry numerical modelling does not hold. Using a fully discrete computational approach in combination with experimental frequency data, this behaviour has been observed for squat masonry towers.
{"title":"Refining the masonry shear modulus in masonry towers via Bayesian model updating","authors":"Silvia Monchetti , Gianni Bartoli , Michele Betti , Siro Casolo , Francesco Clementi","doi":"10.1016/j.probengmech.2026.103903","DOIUrl":"10.1016/j.probengmech.2026.103903","url":null,"abstract":"<div><div>The assessment of the structural behavior of masonry towers often involves developing and identifying computational models to be used to perform static and/or time-history nonlinear analyses. Such models are frequently developed assuming isotropic masonry behavior, with model identification carried out - based on available experimental data - through deterministic tuning of a limited set of representative parameters. However, masonry exhibits complex structural textures that often deviate significantly from the commonly made assumption of isotropic behaviour. This paper investigates the effects of this assumption by focusing on the masonry shear modulus. A two-dimensional rigid body–spring model is adopted as computational approach, as its fully discrete formulation allows overcoming the limitations of the Cauchy continuum while maintaining a reduced computational cost. A Bayesian updating framework based on dynamic experimental data is employed to account for multiple sources of uncertainty, including parameter uncertainty, model inadequacy, and observation errors. Three masonry towers with different slenderness ratios are considered as representative case studies. The adopted Bayesian model updating approach allows for the estimation of their shear modulus while accounting for uncertainties, and the results show that the common assumption of isotropic behaviour in masonry numerical modelling does not hold. Using a fully discrete computational approach in combination with experimental frequency data, this behaviour has been observed for squat masonry towers.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103903"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395577","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-11-19DOI: 10.1016/j.probengmech.2025.103865
M. Grigoriu
Translation models of stationary non-Gaussian processes, referred to as stationary translation models, are extended to models for nonstationary non-Gaussian processes, referred to as nonstationary translation models. The stationary/nonstationary translation models are time-invariant/time-variant memoryless transformations of stationary/nonstationary Gaussian processes. We give and prove properties of nonstationary translation models and specialize them to stationary translation models. Numerous examples are presented to illustrate the construction of nonstationary translation models and highlight some of their properties.
{"title":"Nonstationary translation processes","authors":"M. Grigoriu","doi":"10.1016/j.probengmech.2025.103865","DOIUrl":"10.1016/j.probengmech.2025.103865","url":null,"abstract":"<div><div>Translation models of stationary non-Gaussian processes, referred to as stationary translation models, are extended to models for nonstationary non-Gaussian processes, referred to as nonstationary translation models. The stationary/nonstationary translation models are time-invariant/time-variant memoryless transformations of stationary/nonstationary Gaussian processes. We give and prove properties of nonstationary translation models and specialize them to stationary translation models. Numerous examples are presented to illustrate the construction of nonstationary translation models and highlight some of their properties.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103865"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395574","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2026-02-12DOI: 10.1016/j.probengmech.2026.103898
Tengfei Xing , Xiaodan Ren , Jie Li
Biological tissue materials exhibit high tensile properties and a typical “J-shaped” stress-strain response under external loading due to their unique microstructural architecture. Inspired by this, soft network metamaterials have been widely applied due to their outstanding tensile performance and nonlinear mechanical response. However, achieving both high fracture toughness and defect insensitivity in two-dimensional random soft network metamaterials (2D-RSMs) remains a significant challenge in biomimetic design. In this study, we propose a 2D-RSM with tunable fracture toughness and defect insensitivity. To comprehensively investigate the probabilistic mechanical response of the material post-failure, we employ the Probability Density Evolution Method (PDEM) to analyze the stochastic toughening mechanism of 2D-RSMs under various characteristic conditions, including random control parameters, microstructural geometric configurations, and scale control vectors. We found that, compared to two-dimensional regular soft network metamaterials (2D-SMs), 2D-RSMs leverage randomness to enhance robustness, resulting in remarkable defect insensitivity. This characteristic significantly mitigates the impact of defects on the mechanical response of the material, providing critical reliability for material fabrication and practical applications under complex conditions.
{"title":"Stochastic toughening mechanism of two-dimensional random soft network metamaterials","authors":"Tengfei Xing , Xiaodan Ren , Jie Li","doi":"10.1016/j.probengmech.2026.103898","DOIUrl":"10.1016/j.probengmech.2026.103898","url":null,"abstract":"<div><div>Biological tissue materials exhibit high tensile properties and a typical “J-shaped” stress-strain response under external loading due to their unique microstructural architecture. Inspired by this, soft network metamaterials have been widely applied due to their outstanding tensile performance and nonlinear mechanical response. However, achieving both high fracture toughness and defect insensitivity in two-dimensional random soft network metamaterials (2D-RSMs) remains a significant challenge in biomimetic design. In this study, we propose a 2D-RSM with tunable fracture toughness and defect insensitivity. To comprehensively investigate the probabilistic mechanical response of the material post-failure, we employ the Probability Density Evolution Method (PDEM) to analyze the stochastic toughening mechanism of 2D-RSMs under various characteristic conditions, including random control parameters, microstructural geometric configurations, and scale control vectors. We found that, compared to two-dimensional regular soft network metamaterials (2D-SMs), 2D-RSMs leverage randomness to enhance robustness, resulting in remarkable defect insensitivity. This characteristic significantly mitigates the impact of defects on the mechanical response of the material, providing critical reliability for material fabrication and practical applications under complex conditions.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103898"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395578","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2026-02-12DOI: 10.1016/j.probengmech.2026.103894
Haolin Zhang , Linfang Qian , Longmiao Chen , Taisu Liu , Weiwei Chen , Liu Yang
Time-dependent reliability analysis of mechanisms involving stochastic processes and random variables has always been a challenge. Enhancing the efficiency of reliability analysis for complex nonlinear systems while ensuring accuracy is a critical focus recently. In this study, a double-layer surrogate model-based reliability analysis method is proposed: the first layer constructs a novel deep learning operator net surrogate model to implement uncertainty propagation; the second layer estimates failure probability based on an active learning kriging model. The extended optimal linear estimation discretization method is employed to transform stochastic processes into time-independent random variables. To improve the training efficiency, a sample selection iterative training method is applied to the first-layer model training. Two numerical examples and an engineering application demonstrate the efficiency and accuracy of the proposed method. The methodology proposed in this paper can also be further extended to other reliability analysis problems with high nonlinearity and without an explicit performance function.
{"title":"A novel ESN-DeepONet and kriging double layer model for time-dependent reliability analysis with random process","authors":"Haolin Zhang , Linfang Qian , Longmiao Chen , Taisu Liu , Weiwei Chen , Liu Yang","doi":"10.1016/j.probengmech.2026.103894","DOIUrl":"10.1016/j.probengmech.2026.103894","url":null,"abstract":"<div><div>Time-dependent reliability analysis of mechanisms involving stochastic processes and random variables has always been a challenge. Enhancing the efficiency of reliability analysis for complex nonlinear systems while ensuring accuracy is a critical focus recently. In this study, a double-layer surrogate model-based reliability analysis method is proposed: the first layer constructs a novel deep learning operator net surrogate model to implement uncertainty propagation; the second layer estimates failure probability based on an active learning kriging model. The extended optimal linear estimation discretization method is employed to transform stochastic processes into time-independent random variables. To improve the training efficiency, a sample selection iterative training method is applied to the first-layer model training. Two numerical examples and an engineering application demonstrate the efficiency and accuracy of the proposed method. The methodology proposed in this paper can also be further extended to other reliability analysis problems with high nonlinearity and without an explicit performance function.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103894"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395673","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2025-12-23DOI: 10.1016/j.probengmech.2025.103885
Wenchao Ding, Feng Li, Zhaojie Yu, Fei Cheng
An adaptive subinterval univariate dimension reduction method is proposed to address uncertain problems involving large interval parameters, where the traditional subinterval method often suffers from prohibitive computing costs. The key innovation lies in the development of an adaptive iterative partitioning strategy guided by sensitivity analysis, which dynamically decomposes the original large interval into smaller subintervals. Within each subinterval, the univariate dimension reduction method is used to estimate the response bounds, leveraging the high-order terms in the Taylor expansion series to improve both precision and convergence speed. Three numerical examples demonstrate that the proposed method significantly reduces computational cost while achieving higher accuracy.
{"title":"An adaptive subinterval univariate dimension reduction method for uncertain problems with large interval parameters","authors":"Wenchao Ding, Feng Li, Zhaojie Yu, Fei Cheng","doi":"10.1016/j.probengmech.2025.103885","DOIUrl":"10.1016/j.probengmech.2025.103885","url":null,"abstract":"<div><div>An adaptive subinterval univariate dimension reduction method is proposed to address uncertain problems involving large interval parameters, where the traditional subinterval method often suffers from prohibitive computing costs. The key innovation lies in the development of an adaptive iterative partitioning strategy guided by sensitivity analysis, which dynamically decomposes the original large interval into smaller subintervals. Within each subinterval, the univariate dimension reduction method is used to estimate the response bounds, leveraging the high-order terms in the Taylor expansion series to improve both precision and convergence speed. Three numerical examples demonstrate that the proposed method significantly reduces computational cost while achieving higher accuracy.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103885"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145840782","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-01Epub Date: 2026-01-13DOI: 10.1016/j.probengmech.2026.103893
Xiaobo Zhang
Local reliability sensitivity (LRS), defined as the partial derivatives of structural failure probability with respect to limit state parameters and distribution parameter, provides critical gradient information for reliability-based design optimization and safety-informed decision-making. LRS analysis typically relies on a post-processing step of existing reliability analysis results. In this study, LRS analysis in complex failure domains is investigated using the Directional Sampling (DS) framework. DS achieves computational efficiency through unit-hypersphere sampling and conditional one-dimensional reliability analysis. However, complex failure domains, characterized by disconnected regions, high nonlinearity, or multiple failure modes, may have multiple intersections between a sampled direction and the limit-state surface, posing significant challenges to both reliability and LRS analysis. Existing DS-based sensitivity approaches lack consideration of complex failure domains. Therefore, for LRS with respect to limit-state parameters, an improved domain integral method within the DS framework is developed, transforming the surface integral over complex boundaries into a tractable domain integral coupled with sequential conditional sensitivity analysis. For LRS with respect to distribution parameters, an extended score function method within the DS framework is derived by introducing a directional score function, enabling sensitivity estimation without requiring additional limit-state function calls. The accuracy and efficiency of the proposed methods are validated through challenging examples with complex failure domains, including the modified Rastrigin function and a nonlinear oscillator under white noise base excitation.
{"title":"Local reliability sensitivity analysis with directional sampling framework in complex failure domain","authors":"Xiaobo Zhang","doi":"10.1016/j.probengmech.2026.103893","DOIUrl":"10.1016/j.probengmech.2026.103893","url":null,"abstract":"<div><div>Local reliability sensitivity (LRS), defined as the partial derivatives of structural failure probability with respect to limit state parameters and distribution parameter, provides critical gradient information for reliability-based design optimization and safety-informed decision-making. LRS analysis typically relies on a post-processing step of existing reliability analysis results. In this study, LRS analysis in complex failure domains is investigated using the Directional Sampling (DS) framework. DS achieves computational efficiency through unit-hypersphere sampling and conditional one-dimensional reliability analysis. However, complex failure domains, characterized by disconnected regions, high nonlinearity, or multiple failure modes, may have multiple intersections between a sampled direction and the limit-state surface, posing significant challenges to both reliability and LRS analysis. Existing DS-based sensitivity approaches lack consideration of complex failure domains. Therefore, for LRS with respect to limit-state parameters, an improved domain integral method within the DS framework is developed, transforming the surface integral over complex boundaries into a tractable domain integral coupled with sequential conditional sensitivity analysis. For LRS with respect to distribution parameters, an extended score function method within the DS framework is derived by introducing a directional score function, enabling sensitivity estimation without requiring additional limit-state function calls. The accuracy and efficiency of the proposed methods are validated through challenging examples with complex failure domains, including the modified Rastrigin function and a nonlinear oscillator under white noise base excitation.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"83 ","pages":"Article 103893"},"PeriodicalIF":3.5,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022989","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-01Epub Date: 2025-11-21DOI: 10.1016/j.probengmech.2025.103867
Kaixuan Feng , Zhenzhou Lu
The Bayesian updating method is an effective approach for calibrating model characteristics when new observational data become available. In most existing Bayesian updating methods, the distribution parameters of random model inputs are considered constants. However, these parameters may themselves be uncertain due to limited knowledge of the model inputs. Therefore, a new Bayesian updating method is developed herein considering the uncertainty in the distribution parameters of random model inputs. In the proposed method, new observations may include input data, output data, or a combination of both. The principal contribution of this work lies in the adaptive construction of the likelihood function for the distribution parameters based on different sources of observations. Using the likelihood function and the prior probability density function (PDF) of the distribution parameters, the posterior PDF of these parameters is first obtained. Subsequently, the posterior PDF of the model output can be derived via either a direct or an indirect approach. The theoretical equivalence of these two perspectives is demonstrated. Finally, an example is provided to illustrate the feasibility and validity of the proposed method.
{"title":"Bayesian updating method considering the uncertainty of distribution parameters of random model inputs","authors":"Kaixuan Feng , Zhenzhou Lu","doi":"10.1016/j.probengmech.2025.103867","DOIUrl":"10.1016/j.probengmech.2025.103867","url":null,"abstract":"<div><div>The Bayesian updating method is an effective approach for calibrating model characteristics when new observational data become available. In most existing Bayesian updating methods, the distribution parameters of random model inputs are considered constants. However, these parameters may themselves be uncertain due to limited knowledge of the model inputs. Therefore, a new Bayesian updating method is developed herein considering the uncertainty in the distribution parameters of random model inputs. In the proposed method, new observations may include input data, output data, or a combination of both. The principal contribution of this work lies in the adaptive construction of the likelihood function for the distribution parameters based on different sources of observations. Using the likelihood function and the prior probability density function (PDF) of the distribution parameters, the posterior PDF of these parameters is first obtained. Subsequently, the posterior PDF of the model output can be derived via either a direct or an indirect approach. The theoretical equivalence of these two perspectives is demonstrated. Finally, an example is provided to illustrate the feasibility and validity of the proposed method.</div></div>","PeriodicalId":54583,"journal":{"name":"Probabilistic Engineering Mechanics","volume":"82 ","pages":"Article 103867"},"PeriodicalIF":3.5,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145617948","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}