This application claims the priority benefit of China application serial no. 202311544897.X, filed on Nov. 20, 2023. The entirety of the above-mentioned patent application is hereby incorporated by reference herein and made a part of this specification.
The invention relates to the technical field of seismic performance evaluation of buildings and, in particular, to a ground motion intensity measure optimization method for seismic response prediction.
The safety of buildings is a key issue in the advancement of urbanization in China, due to frequent earthquake disasters. Once a building structure is damaged, it can cause significant economic losses and even casualties. Therefore, there is an urgent need to develop a reliable and rapid seismic damage assessment method for buildings to provide technical support for post-earthquake emergency management decisions. Machine learning-based seismic damage prediction of buildings, with high accuracy and greatly shortens the assessment time, is the future development direction of design and risk assessment management for urban buildings. Although the machine learning model is not mathematically limited by the dimensions of ground motion intensity measures, it is still very important to optimize intensity measures. On the one hand, considering the complexity of ground motions, numerous intensity measures are needed to characterize the features of ground shaking, so as to train practical machine learning models to predict the complex response of buildings during strong earthquake events. On the other hand, considering numerous intensity measures as input variables will lead to a high-dimensional probability space mapping problem, it requires a large number of samples to train machine learning models, and even the models cannot be accurately constructed. Existing ground motion intensity measure optimization methods often focus on a single intensity measure, and rely on an assumption that intensity measures and structural demand parameters exhibit a logarithmic linear correlation. However, structural damage often exhibits inelastic behaviors, and especially for complex high-rise buildings, it greatly limits the applicability of traditional ground motion intensity optimization methods. Therefore, by making full use of the advantages of data drive, a combination of ground motion intensity measures that are closely related to the responses of buildings are selected to improve the construction efficiency and prediction accuracy of machine learning models.
Objectives of the invention: An objective of the invention is to provide a ground motion intensity measure optimization method for seismic response prediction, which avoids the assumption of logarithmic linearity in traditional methods and the limitation of intensity measure dimensions. The method can be performed by a computer including a processor and a memory. The memory stores instructions which are executed by the processor to perform the method.
Technical solution: The ground motion intensity measure optimization method for seismic response prediction according to the invention comprises the following steps:
Further, in step (1), the candidate ground motion intensity measures IMs are divided into structure-independent or structure-dependent related to acceleration, related to velocity, related to displacement, and related to time based on their physical definition; according to the fortification level, site condition and design characteristic period, actual ground motion records that match the site-specific target response spectrum are obtained.
Further, in an embodiment, in step (2), there are 12 structural parameters SPs of the building, including: bulk density of concrete, compressive strength of core concrete, peak strain of core concrete, ultimate strain of core concrete, compressive strength of cover concrete, ultimate strain of cover concrete, diameter of steel bars in columns, diameter of Φ20 steel bars in beams, diameter of Φ18 steel bars in beams, yield strength of steel bars, elastic modulus of steel bars, and hardening rate of steel bars; the probability distribution parameters include mean, variance and distribution type, and then 12 structural parameter variables are sampled respectively by LHS to generate N random samples, which are finally combined into N pieces of building sample data.
Further, in step (3), the N numerical models of the building are randomly matched with N ground motion records; the damage measures DM of the building refer to a maximum inter-story drift ratio.
Further, the post-earthquake building damage database in step (3) is built by: using the ground motion intensity measures and the structural parameters of the building as input samples; and the damage measures of the building as output samples.
Further, in an embodiment, step 4 is implemented as follows:
Further, in step (5), the efficiency of an intensity measure is represented by the goodness of fit (R2) of the elastic-net model, which is calculated by:
where y and ŷ are a true value and a predicted value, respectively.
The practicality of an intensity measure is represented by a variable regression coefficient β of the elastic-net model; a linear regression model over each intensity measure and a corresponding damage measure of the building is constructed to obtain a relative error (i.e., residual) between the predicted value and the true value; the adequacy of the intensity measure is characterized according to a significance test between the corresponding residual of the intensity measure and seismic characteristics (i.e., magnitude and rupture distance). It is generally believed that when the p-value of the significance test level is greater than 0.05, the intensity measure meets the adequacy requirement.
Further, the machine learning method in step (6) includes K nearest neighbor (KNN), Gaussian process regression (GPR), support vector machine (SVM), artificial neural network (ANN), random forest tree (RF) and extreme gradient boosting tree (XGBoost).
Further, in step (6), the input samples of the database are x=[xIMs*, xSPs], including the optimal ground motion intensity measures and the structural parameters, and the output samples are damage measures yDM of the building.
the database is standardized and then randomly divided into a training set and a test set according to a P:Q ratio for training.
Further, in step (7), mean square error (MSE) is used as an assessment measure to determine the optimal hyperparameters of the machine learning models as follows:
where y and ŷ are a true value and a predicted value, respectively.
The goodness of fit (R2) and the root mean square error (RMSE) of the test set are used to demonstrate the prediction performance of the machine learning models, wherein the RMSE is calculated by:
where y and ŷ are a true value and a predicted value, respectively.
Beneficial effects: Compared with the prior art, the invention has the following significant advantages: based on a damage database of the buildings under strong earthquake excitations, a combination of optimal ground motion intensity measures is selected by an elastic-net technology, avoiding the assumption of logarithmic linearity in traditional methods and the limitation of intensity measure dimensions. Using the optimal ground motion intensity measures as input can significantly improve the prediction accuracy of the machine learning models as compared with using a single intensity measure as input, and can greatly reduce the number of samples required for training as compared with using all candidate intensity measures as input. The application of the invention can be easily extended to the rapid prediction of post-earthquake damage states of bridges, dams, tunnels and other structures.
The technical solution of the invention is further described below in conjunction with the accompanying drawings.
Referring to
(1) Determining candidate intensity measures IMs by the collection and sorting of literature, obtaining N actual ground motion records that match the site-specific target response spectrum, and calculating candidate intensity measures corresponding to each ground motion record. In an embodiment, according to the site design requirements (i.e., fortification level, site condition, and design characteristic period) of a study area, 320 ground motion records are selected from an earthquake database of the Pacific Earthquake Engineering Research Center, and the acceleration response spectrum of the ground motion records are shown in
(2) Determining structural parameters SPs of a building and probability distribution, and carrying out Latin Hypercube Sampling (LHS) to obtain N building samples that meet parameter probability distribution. In an embodiment, taking a concrete frame shown in
(3) Using OpenSees software to establish 320 numerical models of the building corresponding to the building samples and obtaining N damage measures DMs of the building by nonlinear time-history analysis to build a post-earthquake building damage database. In an embodiment, a finite element model is established by means of OpenSees software, as shown in
(4) Constructing a variable selection model based on the elastic-net algorithm according to the post-earthquake building damage database. In an embodiment:
(5) Calculating the efficiency and practicality of input variables by means of an elastic-net model, then selecting a combination of optimal intensity measures IMs*, and testing the adequacy of the intensity measures. In an embodiment, the hyperparameter λ for the elastic-net algorithm is determined according to the goodness of fit (R2) of the test set to ensure the efficiency of optimal intensity measure selection. R2 is calculated by:
where y and ŷ are a true value and a predicted value, respectively.
The practicality of an intensity measure is assessed according to the regression coefficient corresponding to the intensity measure. A linear regression model over each intensity measure and a corresponding damage measure of the building is constructed to obtain a relative error (i.e., residual) between the predicted value and the true value of the damage measure; the p-value of the significance test level between the corresponding residual of the intensity measure and seismic characteristics (i.e., seismic magnitude and seismic distance) is calculated. intensity measures with p-values less than or equal to 0.05 are screened out to meet the adequacy requirement. Table 3 lists the combination of optimal ground motion intensity measures and the practical ranking of the ground motion intensity measures.
(6) Based on the post-earthquake building damage database, performing model training by means of a machine learning method to obtain prediction accuracy measures. In an embodiment,
Taking the combination of 320 optimal ground motion intensity measures xIMs* and corresponding structural parameters xSPs of the building as input samples, and the corresponding damage measures yDM of the building as output samples, the database is randomly divided into a training set and a test set according to a ratio of 7:3. Based on the training set, KNN, GPR, SVM, ANN, RF and XGBoost are used for model training respectively. The optimal values of the key hyperparameters listed in Table 4 are determined by grid search technique to obtain a most accurate machine learning prediction model.
In order to demonstrate the advantages of using the combination of optimal intensity measures as input for the machine learning model to predict the seismic damage response of a building, a common intensity measure PGA, a single optimal intensity measure Sa−1.0 s, a combination of optimal intensity measures and all candidate intensity measures are respectively used as input features to train a machine learning model. Taking an SVM model as an example, as shown in
(7) Based on assessment measures for prediction model performance, determining optimal hyperparameters and corresponding machine learning models by grid search technique. In an embodiment, mean square error MSE is used as an assessment measure to determine the optimal hyperparameters of the machine learning model as follows:
where y and ŷ are a true value and a predicted value, respectively.
the goodness of fit R2 and the root mean square error RMSE of the test set are used to demonstrate the prediction performance of the machine learning models, wherein the RMSE is calculated by:
where y and ŷ are a true value and a predicted value, respectively.
The goodness of fit (R2) and root mean square error RMSE of six machine learning models on the test set are shown in
(8) Rapidly predicting a structural damage state of the building under earthquake damage by means of the trained machine learning models. In an embodiment, by inputting the optimal ground motion intensity measures and basic information data of the building into the trained machine learning models, the maximum inter-story drift ratio of the target building can be quickly obtained. According to damage states thresholds recommended by the specification: 0-0.002 means no damage; 0.002-0.01 means slight damage; 0.01-0.02 means moderate damage; 0.02-0.04 means severe damage; and greater than 0.04 means complete collapse, the damage states of the building under given ground motion intensity measures can be assessed.
Number | Date | Country | Kind |
---|---|---|---|
202311544897.X | Nov 2023 | CN | national |
Number | Name | Date | Kind |
---|---|---|---|
20200025957 | Guo | Jan 2020 | A1 |
Number | Date | Country |
---|---|---|
109359311 | Feb 2019 | CN |
Entry |
---|
Chen, Libo et al., “Comparison and selection of ground motion intensity measures in probabilistic seismic demand analysis”, Journal of Earthquake Engineering and Engineering Vibration, vol. 32, No. 6, Dec. 15, 2012, with English abstract thereof, pp. 23-31. |
Wentao Wang et al., “Machine learning-based collapse prediction for post-earthquake damaged RC columns under subsequent earthquakes”, Soil Dynamics and Earthquake Engineering, vol. 172, May 26, 2023, pp. 1-16. |