In engineering practice, Bayesian model updating using field data is often conducted to reduce the substantial inherent epistemic uncertainties in geomaterial properties resulting from complex geological processes. The Bayesian Updating with Subset simulation (BUS) method is commonly employed for this purpose. However, the wealth of field data available for engineers to interpret can lead to challenges associated with the “curse of dimensionality”. Specifically, the value of the likelihood function in the BUS method can become extremely small as the volume of field data increases, potentially falling below the accuracy threshold of computer floating-point operations. This undermines both the computational efficiency and accuracy of Bayesian model updating. To effectively address this technical challenge, this paper proposes an improved BUS method developed based on parallel system reliability analysis. Leveraging the Cholesky decomposition-based midpoint method, the total failure domain in the original BUS method, which involves a low acceptance rate, is subdivided into several sub-failure domains with a high acceptance rate. Facilitated with an improved Metropolis-Hastings algorithm, the improved BUS method enables the consideration of a large volume of field data and spatial variability of geomaterial properties in the probabilistic back analysis. The results of an illustrative soil slope, involving spatially variable undrained shear strength, demonstrate that the improved BUS method is effective in simultaneously incorporating a substantial volume of field measurements and observations in the model updating process. Through a comparison with the original BUS method, the improved BUS method is shown to be useful for Bayesian model updating of high-dimensional spatially variable geomaterial properties and slope reliability assessment.
The seismic engineering demand parameters (EDPs) of building clusters exhibit significant spatial correlations and need full consideration in regional risk and reliability assessments. This study presents an efficient scheme to determine the joint distribution of multi-structure EDPs, which captures all EDP correlations and enables direct calculation of system reliability for building clusters. This scheme generates spatially correlated random ground motion fields through ground motion cross power spectrum density (PSD) models with stochastic harmonic function simulations. Subsequently, the decoupled multi-probability density evolution method (M−PDEM) is integrated to conduct seismic analysis of building clusters under random ground motion fields to determine their EDP joint distribution. An example of three linear single-degree-of-freedom (SDOF) models shows that the proposed scheme requires only hundreds of analyses to achieve the same accuracy as 105 Monte Carlo Simulation (MCS) analyses, while also capturing the nonlinear correlations among EDPs. Finally, an engineering application of three reinforced concrete (RC) frame shear-wall buildings under a rare earthquake scenario is investigated, and the joint collapse probability by the scheme is compared with that by commonly-adopted assumptions of mutual independence and linear correlation. The results reveal that relative errors by the two assumptions can reach up to 39 % and 22 %, respectively.
The life-cycle seismic resilience assessment of sea-crossing highway bridges plays a crucial role in guiding decisions for their long-term operation, maintenance, and rehabilitation. Due to the inherently stochastic nature of marine environments, evaluating the resilience of bridges while considering all possible environmental scenarios throughout their service life necessitates substantial computational efforts and presents practical challenges. Thus, this study develops a three-stage framework for predicting the life-cycle seismic resilience of sea-crossing highway bridges. Stochastic models for marine environmental conditions and bridge durability are developed and validated using experimental measurement data. A modified Good Lattice Point-Partially Stratified Sampling (GLP-PSS) method is employed to generate a uniform and limited number of samples. A typical prestressed concrete sea-crossing highway bridge is selected as the benchmark bridge, and parameterized numerical models are established using 460 representative environmental parameter samples on the OpenSees platform. Leveraging the environmental model and material properties, the durability of the bridge is predicted over its service life. Nonlinear time history analyses are carried out for each bridge model using 120 real ground motion records, which allow the identification of variations in seismic demands, capacities, and system fragilities at different time intervals. Subsequently, the life-cycle seismic resilience of the bridge is predicted utilizing surrogate models based on the response surface method (RSM) and artificial neural networks (ANN), respectively. Finally, the time-dependent probabilistic characteristics of seismic resilience are thoroughly discussed. Results indicate that ANN demonstrates a higher degree of generalization capability in predicting the life-cycle seismic resilience. Focusing solely on changes in mean resilience over the service time may lead to an underestimation of bridge resilience, as it may ignore the tails of its distribution, potentially resulting in an overestimation of bridge resilience. Furthermore, global warming may expedite the decline in resilience.
Acquiring engineering data is frequently expensive, resulting in sparse data that may lead to a lack of knowledge for design and analysis. Thus, it is not always feasible to precisely determine the probability density functions (PDFs) of uncertain model parameters. Under such circumstances that involve simultaneous aleatory and epistemic uncertainties, repeated uncertainty propagation (UP) analysis is generally required. In this paper, a novel approach for hybrid UP is proposed by integrating B-spline chaos and augmented change of probability measure (aCOM) for meeting different goals. The B-spline chaos is adopted to represent the complicated computational model as a function of an arbitrary input random variable, while the aCOM is employed to reconstruct the PDF of the model output when the input PDF is changed due to epistemic uncertainty. In the case of small epistemic uncertainty, hybrid UP can be achieved directly by changing the assigned probabilities of existing sample results. While in the case of large epistemic uncertainty, additional samples from an augmenting PDF are generated. The proposed method is compatible with both cases. The numerical algorithm of the proposed method is presented and illustrated by four benchmark problems. Further, the accuracy and efficiency of the proposed method are substantiated by four numerical examples compared with analytical solutions or Monte Carlo simulations. An attempt to enhance the proposed method with the aid of active subspace methods to handle high-dimensional problems is also discussed in this work. The limitations and potential improvements of the proposed approach are outlined as well.
A novel spectral incremental dynamic analysis methodology for analysing structural response in nonlinear systems with fractional derivative elements is presented, aligning with modern seismic design codes, like Eurocode 8. Drawing inspiration from the concept of fully non-stationary stochastic processes, the vector of the imposed seismic excitations is characterised by time and frequency evolving power spectra stochastically compatible with elastic response spectra of specified damping ratio and ground acceleration. The proposed method efficiently determines the nonlinear system time-dependent probability density functions for the non-stationary system response amplitude by employing potent nonlinear stochastic dynamics concepts, such as stochastic averaging and statistical linearisation. Unlike traditional incremental dynamic analysis curves found in the literature, the herein proposed method introduces a three-dimensional alternative counterpart, that of stochastic engineering demand parameter surfaces, providing with higher-order statistics of the system response. An additional noteworthy aspect involves the derivation of response evolutionary power spectra as function of spectral acceleration, offering a deeper insight into the underlying system dynamics. Besides its capabilities, the method maintains the coveted element of a particularly low associated computational cost, increasing its attractiveness and practicality among diverse applications of engineering interest. Numerical examples comprising the bilinear hysteretic model endowed with fractional derivative elements subject to an Eurocode 8 elastic design spectrum demonstrate the capabilities and reliability of the proposed methodology. Its accuracy is assessed by juxtaposing the derived results with germane Monte Carlo Simulation data.
Due to the characteristics of high stiffness-weight ratio, timber floors are prone to annoying vibrations under human excitation. Given the natural origin of timber, its mechanical properties exhibit significant variability. The randomness inherent in human excitation cannot be overlooked during structural dynamic analysis. Consequently, the adoption of a stochastic approach is imperative for attaining reliable serviceability evaluation results. However, current research on human-induced vibrations in the timber floor, accounting for this randomness, remains inadequate. In this paper, an experimental investigation is conducted on the dynamic properties and human-induced responses of a timber floor composed of glued laminated timber and oriented strand board. A finite element model is developed and subsequently validated for accuracy in terms of modal properties and dynamic responses. The probability density evolution method is introduced for stochastic analysis, which demonstrates that both the natural frequency and dynamic responses of the floor exhibit considerable variability when uncertainty factors are considered. The Kullback–Leibler divergence indices are used to assess the impact of each uncertain variable quantitatively. The results indicate that the longitudinal elastic modulus has the greatest influence on the natural frequency, while the first dynamic load factor, αz1, exerts the most significant impact on dynamic responses. Notably, both material mechanical properties and load model parameters contribute to the uncertainty of dynamic responses, with the influence of the load model parameters being significantly greater than that of material mechanical properties.
The reliability assessment of structural systems presents a significant challenge in structural engineering. A commonly employed approximation is the First-Order System Reliability Method (FOSRM), which estimates system reliability using the FORM component reliabilities and sensitivity factors. An essential step in FORM involves transforming the random vector into the standard vector , often using the Rosenblatt transformation (RT). Several studies demonstrated that different conditioning orders in the RT yield different FORM component results. This study investigates how these differences on component level propagate into the FOSRM system level. We conducted several typical engineering case studies with various failure probabilities, system sizes, and dependency structures (Gaussian and Frank Copula). For the Frank Copula, different Rosenblatt conditioning orders systematically yielded different FOSRM results, with most cases showing differences between 10% and 30% in estimated failure probability. For some systems, these differences increased with system size, suggesting that greater variations might be observed for larger systems. Notably, systems with Gaussian Copula functions also proved vulnerable to the Rosenblatt conditioning order when different components were assessed with different conditioning orders. The observed differences were larger than previously reported and should be carefully considered in uniform safety assessments.