Large‐scale real‐time hybrid simulation involving multiple experimental substructures and adaptive actuator delay compensation

Real‐time hybrid simulation provides a viable method to experimentally evaluate the performance of structural systems subjected to earthquakes. The structural system is divided into substructures, where part of the system is modeled by experimental substructures, whereas the remaining part is modeled analytically. The displacements in a real‐time hybrid simulation are imposed by servo‐hydraulic actuators to the experimental substructures. Actuator delay compensation has been shown by numerous researchers to vitally achieve reliable real‐time hybrid simulation results. Several studies have been performed on servo‐hydraulic actuator delay compensation involving single experimental substructure with single actuator. Research on real‐time hybrid simulation involving multiple experimental substructures, however, is limited. The effect of actuator delay during a real‐time hybrid simulation with multiple experimental substructures presents challenges. The restoring forces from experimental substructures may be coupled to two or more degrees of freedom (DOF) of the structural system, and the delay in each actuator must be adequately compensated. This paper first presents a stability analysis of actuator delay for real‐time hybrid simulation of a multiple‐DOF linear elastic structure to illustrate the effect of coupled DOFs on the stability of the simulation. An adaptive compensation method then proposed for the stable and accurate control of multiple actuators for a real‐time hybrid simulation. Real‐time hybrid simulation of a two‐story four‐bay steel moment‐resisting frame with large‐scale magneto‐rheological dampers in passive‐on mode subjected to the design basis earthquake is used to experimentally demonstrate the effectiveness of the compensation method in minimizing actuator delay in multiple experimental substructures. Copyright © 2011 John Wiley & Sons, Ltd.     

应用人工神经网络技术的大型斜拉桥子结构损伤识别研究

本文应用人工神经网络技术对大型斜拉桥结构进行了子结构损伤识别研究。文中首先介绍了子结构损伤识别的基本方法，然后应用自组织竞争神经网络建立了对于大型桥梁结构识别子结构损伤情况的子结构损伤识别方法，并且应用BP网络进一步建立了大型桥梁结构各子结构内部的损伤位置和损伤程度的识别方法，数值模拟了一大跨度斜拉桥子结构损伤以及子结构内部损伤的识别过程，最后得出结论：(1)基于自组织竞争网络的子结构损伤识别方法能迅速准确地识别大型结构的损伤情况；(2)基于BP网络所建立的结构损伤识别方法，能对子结构中结构损伤的位置和程度进行进一步的识别；(3)基于人工神经网络技术的结构损伤识别方法是大型土木工程结构损伤识别的有效方法，可在工程结构损伤识别中广泛应用。     

Vibration‐based damage assessment of steel structure using global and local response measurements

Damage assessment of a structure involves acquiring and identifying dynamic characteristics of the structure and using these characteristics to evaluate behavior and performance. In this study, an unsymmetrical three‐story steel structure (fabricated with one weak column in the first floor) was tested on shaking table and subjected to a series of earthquake excitations with increasing level of excitation back to back. Besides, white noise excitation was also applied in between the earthquake excitation to serve as the reference state. Both the traditional sensing system (accelerometer and linear variable differential transformer) and the local optical tracker system were implemented in the structure to collect the vibration‐based responses. For operational modal analysis, structural response from white noise excitation will be used in this study. First, the traditional system identification using global response data is used (multivariate autoregressive (AR)‐model) to extract system natural frequencies and mode shapes from all different set of white noise responses after earthquake excitation. The migration of AR‐coefficient ellipse error from each sensor response was used to identify the damage location. Second, blind source separation technique was used to identify the modal contribution of the structure from each test, which provide information to detect the damage severity. Finally, from the local optical tracker array data, the principal component analysis was applied to quantify the earthquake‐induce local stress of the structural member. Combine the result from damage detection using global measurement and the identified local element stress, one can locate and quantify the damage. Copyright © 2015 John Wiley & Sons, Ltd.     

Post‐earthquake bridge repair cost and repair time estimation methodology

While structural engineers have traditionally focused on individual components (bridges, for example) of transportation networks for design, retrofit, and analysis, it has become increasingly apparent that the economic costs to society after extreme earthquake events are caused at least as much from indirect costs as direct costs due to individual structures. This paper describes an improved methodology for developing probabilistic estimates of repair costs and repair times that can be used for evaluating the performance of new bridge design options and existing bridges in preparation for the next major earthquake. The proposed approach in this paper is an improvement on previous bridge loss modeling studies—it is based on the local linearization of the dependence between repair quantities and damage states so that the resulting model follows a linear relationship between damage states and repair points. The methodology uses the concept of performance groups (PGs) that account for damage and repair of individual bridge components and subassemblies. The method is validated using two simple examples that compare the proposed method to simulation and previous methods based on loss models using a power–law relationship between repair quantities and damage. In addition, an illustration of the method is provided for a complete study on the performance of a common five‐span overpass bridge structure in California. Intensity‐dependent repair cost ratios (RCRs) and repair times are calculated using the proposed approach, as well as plots that show the disaggregation of repair cost by repair quantity and by PG. This provides the decision maker with a higher fidelity of data when evaluating the contribution of different bridge components to the performance of the bridge system, where performance is evaluated in terms of repair costs and repair times rather than traditional engineering quantities such as displacements and stresses. Copyright © 2009 John Wiley & Sons, Ltd.     

Residual structural capacity evaluation of steel moment‐resisting frames with dynamic‐strain‐based model updating method

This paper presents a method for evaluating the residual structural capacity of earthquake‐affected steel structures. The method first quantifies the damage severity of a beam by computing the dynamic‐strain‐based damage index. Next, the model used to analyze the structure is updated based on the damage index, to reflect the observed damage conditions. The residual structural capacity is then estimated in terms of changes in stiffness and strength, which can be applied by structural engineers, via a nonlinear static analysis of the updated model. The main contributions of this paper are in performance evaluation of the dynamic‐strain‐based damage index for seismically induced damage using a newly developed substructure testing environment, consideration of various damage patterns in composite beams, and extension of a local damage evaluation technique to a residual capacity estimation procedure by incorporating the model‐updating technique. In laboratory testing, the specimens were damaged quasi‐statically, and vibration tests were conducted as the damage proceeded. First, a bare steel beam–column connection was tested, and then a similar one with a floor slab was used for a more realistic case. The estimated residual structural capacities for these specimens were compared with the static test results. The results verified that the proposed method can provide fine estimates of the stiffness and strength deteriorations within 10% for the specimen without the floor slab and within 30% for that with the floor slab. Copyright © 2017 John Wiley & Sons, Ltd.     

LMS‐based structural health monitoring of a non‐linear rocking structure

A structure's health or level of damage can be monitored by identifying changes in structural or modal parameters. This research directly identifies changes in structural stiffness due to modelling error or damage for a post‐tensioned pre‐cast reinforced concrete frame building with rocking beam column connections and added damping and stiffness (ADAS) elements. A structural health monitoring (SHM) method based on adaptive least mean squares (LMS) filtering theory is presented that identifies changes from a simple baseline model of the structure. This method is able to track changes in the stiffness matrix, identifying when the building is (1) rocking, (2) moving in a hybrid rocking–elastic regime, or (3) responding linearly. Results are compared for two different LMS‐based SHM methods using an L ₂ error norm metric. In addition, two baseline models of the structure, one using tangential stiffness and the second a more accurate bi‐linear stiffness model, are employed. The impact of baseline model complexity is then delineated. The LMS‐based methods are able to track the non‐linearity of the system to within 15% using this metric, with the error due primarily to filter convergence rates as the structural response changes regimes while undergoing the El Centro ground motion record. The use of a bi‐linear baseline model for the SHM problem is shown to result in error metrics that are at least 50% lower than those for the tangential baseline model. Errors of 5–15% with this L ₂ error norm are fairly stringent compared to the greater than 2 × changes in stiffness undergone by the structure, however, in practice the usefulness of the results is dependent on the resolution required by the user. The impact of sampling rate is shown to be negligible over the range of 200–1000Hz, along with the choice of LMS‐based SHM method. The choice of baseline model and its level of knowledge about the actual structure is seen to be the dominant factor in achieving good results. The methods presented require 2.8–14.0 Mcycles of computation and therefore could easily be implemented in real time. Copyright © 2005 John Wiley & Sons, Ltd.     

Damage identification of a 3D full scale steel–concrete composite structure with partial‐strength joints at different pseudo‐dynamic load levels

Partial‐strength composite steel–concrete moment‐resisting (MR) frame structures represent an open research field in seismic design from both a theoretical and an experimental standpoint. Among experimental techniques, vibration testing is a well‐known and powerful technique for damage detection, localization and quantification, where actual modal parameters of a structure at different states can be determined from test data by using system identification methods. However, the identification of semi‐rigid connections in framed structures is limited, and hence this paper focuses on a series of vibration experiments that were carried out on a realistic MR frame structure, following the application of pseudo‐dynamic and quasi‐static cyclic loadings at the European laboratory for structural assessment of the Joint Research Centre at Ispra, Italy, with the scope of understanding the structural behaviour and identifying changes in the dynamic response. From the forced vibration response, natural frequencies, damping ratios, modal displacements and rotations were extracted using the circle fitting technique. These modal parameters were used for local and global damage identification by updating a 3D finite element model of the intact structure. The identified results were then correlated with observations performed on the structure to understand further the underlying damage mechanisms. Finally, the latin hypercube sampling technique, a variant of the Monte Carlo method, was employed in order to study the sensitivity of the updated parameters of the 3D model to noise on the modal inputs. Copyright © 2009 John Wiley & Sons, Ltd.     

Two‐stage damage detection algorithms of structure using modal parameters identified from recursive subspace identification

Structural damage assessment under external loading, such as earthquake excitation, is an important issue in structural safety evaluation. In this regard, appropriate data analysis and feature extraction techniques are required to interpret the measured data and to identify the state of the structure and, if possible, to detect the damage. In this study, the recursive subspace identification with Bona‐fide LQ renewing algorithm (RSI‐BonaFide‐Oblique) incorporated with moving window technique is utilized to identify modal parameters such as natural frequencies, damping ratios, and mode shapes at each instant of time during the strong earthquake excitation. From which the least square stiffness method (LSSM) combined with the model updating technique, called efficient model correction method (EMCM), is used to estimate the first‐stage system stiffness matrix using the simplified model from the previously identified modal parameters (nominal model). In the second stage, 2 different damage assessment algorithms related to the nominal system stiffness matrix were derived. First, the model updating technique, called EMCM, is applied to correct the nominal model by the newly identified modal parameters during the strong motion. Second, the element damage index can be calculated using element damage index method (EDIM) to quantify the damage extent in each element. Verification of the proposed methods through the shaking table test data of 2 different types of structures and a building earthquake response data is demonstrated to specify its corresponding damage location, the time of occurrence during the excitation, and the percentage of stiffness reduction.     

Analytical model for the in‐plane seismic performance of cold‐formed steel‐framed gypsum partition walls

This paper proposes an experimentally verified procedure to analytically model cold‐formed steel‐framed gypsum nonstructural partition walls considering all the critical components. In this model, the nonlinear behaviors of the connections are represented by hysteretic load‐deformation springs, which have been calibrated using the component‐level experimental data. The studs and tracks are modeled adopting beam elements with their section properties accounting for nonlinear behavior. The gypsum boards are simulated by linear four‐node shell elements. The proposed procedure is implemented to generate the analytical models of three full‐scale partition wall specimens in the OpenSees platform. The specimens were tested as a part of the NEESR‐GC Project on Simulation of the Seismic Performance of Nonstructural Systems. Force‐displacement responses, cumulative dissipated energy, and damage mechanisms from the analytical simulation are compared to the experimental results. The comparison shows that the analytical model accurately predicts the trend of the response as well as the possible damage mechanisms. The procedure proposed here can be adopted in future studies by researchers and also engineers to assess the seismic performance of partition walls with various dimensions and construction details, especially where test data are not available. Copyright © 2015 John Wiley & Sons, Ltd.     

Reliability‐based control algorithms for nonlinear hysteretic systems based on enhanced stochastic averaging of energy envelope

Current reliability‐based control techniques have been successfully applied to linear systems; however, incorporation of stochastic nonlinear behavior of systems in such control designs remains a challenge. This paper presents two reliability‐based control algorithms that minimize failure probabilities of nonlinear hysteretic systems subjected to stochastic excitations. The proposed methods include constrained reliability‐based control (CRC) and unconstrained reliability‐based control (URC) algorithms. Accurate probabilistic estimates of nonlinear system responses to stochastic excitations are derived analytically using enhanced stochastic averaging of energy envelope proposed previously by the authors. Convolving these demand estimates with capacity models yields the reliability of nonlinear systems in the control design process. The CRC design employs the first‐level and second‐level optimizations sequentially where the first‐level optimization solves the Hamilton–Jacobi–Bellman equation and the second‐level optimization searches for optimal objective function parameters to minimize the probability of failure. In the URC design, a single optimization minimizes the probability of failure by directly searching for the optimal control gain. Application of the proposed control algorithms to a building on nonlinear foundation has shown noticeable improvements in system performance under various stochastic excitations. The URC design appears to be the most optimal method as it reduced the probability of slight damage to 8.7%, compared with 11.6% and 19.2% for the case of CRC and a stochastic linear quadratic regulator, respectively. Under the Kobe ground motion, the normalized peak drift displacement with respect to stochastic linear quadratic regulator is reduced to 0.78 and 0.81 for the URC and CRC cases, respectively, at comparable control force levels. Copyright © 2017 John Wiley & Sons, Ltd.     

Identification of modal parameters of non‐classically damped linear structures under multi‐component earthquake loading

A method for parametric system identification of classically damped linear system in frequency domain is adopted and extended for non‐classically damped linear systems subjected up to six components of earthquake ground motions. This method is able to work in multi‐input/multi‐output (MIMO) case. The response of a two‐degree‐of‐freedom model with non‐classical damping, excited by one‐component earthquake ground motion, is simulated and used to verify the proposed system identification method in the single‐input/multi‐output case. Also, the records of a 10 storey real building during the Northridge earthquake is used to verify the proposed system identification method in the MIMO case. In this case, at first, a single‐input/multi‐output assumption is considered for the system and modal parameters are identified, then other components of earthquake ground motions are added, respectively, and the modal parameters are identified again. This procedure is repeated until all four components of earthquake ground motions which are measured at the base level of the building are included in the identification process. The results of identification of real building show that consideration of non‐classical damping and inclusion of the multi‐components effect of earthquake ground motions can improve the least‐squares match between the finite Fourier transforms of recorded and calculated acceleration responses. Copyright © 2006 John Wiley & Sons, Ltd.     

Performance of base‐isolated reinforced concrete buildings under bidirectional seismic excitation considering pounding with retaining walls including friction effects

Seismic pounding of base‐isolated buildings has been mostly studied in the past assuming unidirectional excitation. Therefore, in this study, the effects of seismic pounding on the response of base‐isolated reinforced concrete buildings under bidirectional excitation are investigated. For this purpose, a three‐dimensional finite element model of a code‐compliant four‐story building is considered, where a newly developed contact element that accounts for friction and is capable of simulating pounding with retaining walls at the base, is used. Nonlinear behavior of the superstructure as well as the isolation system is considered. The performance of the building is evaluated separately for far‐fault non‐pulse‐like ground motions and near‐fault pulse‐like ground motions, which are weighted scaled to represent two levels of shaking viz. the design earthquake (DE) level and the risk‐targeted maximum considered earthquake (MCE_R) level. Nonlinear time‐history analyses are carried out considering lower bound as well as upper bound properties of isolators. The influence of separation distance between the building and the retaining walls at the base is also investigated. It is found that if pounding is avoided, the performance of the building is satisfactory in terms of limiting structural and nonstructural damage, under DE‐level motions and MCE_R‐level far‐fault motions, whereas unacceptably large demands are imposed by MCE_R‐level near‐fault motions. In the case of seismic pounding, MCE_R‐level near‐fault motions are found to be detrimental, where the effect of pounding is mostly concentrated at the first story. In addition, it is determined that considering unidirectional excitation instead of bidirectional excitation for MCE_R‐level near‐fault motions provides highly unconservative estimates of superstructure demands. Copyright © 2014 John Wiley & Sons, Ltd.     

Stability analysis of MDOF real‐time dynamic hybrid testing systems using the discrete‐time root locus technique

The time delay resulting from the servo hydraulic systems can potentially destabilize the real‐time dynamic hybrid testing (RTDHT) systems. In this paper, the discrete‐time root locus technique is adopted to investigate the delay‐dependent stability performance of MDOF RTDHT systems. Stability analysis of an idealized two‐story shear frame with two DOFs is first performed to illustrate the proposed method. The delay‐dependent stability condition is presented for various structural properties, time delay, and integration time steps. Effects of delay compensation methods on stability are also investigated. Then, the proposed method is applied to analyze the delay‐dependent stability of a single shaking table RTDHT system with an 18‐DOF finite element numerical substructure, and corresponding RTDHTs are carried out to verify the theoretical results. Furthermore, the stability behavior of a finite element RTDHT system with two physical substructures, loaded by twin shaking tables, is theoretically and experimentally investigated. All experimental results convincingly demonstrate that the delay‐dependent stability analysis on the basis of the discrete‐time root locus technique is feasible. Copyright © 2014 John Wiley & Sons, Ltd.     

Three‐dimensional vibration control of high‐tech facilities against earthquakes and microvibration using hybrid platform

High‐tech equipments engaged in the production of ultra‐precision products have very stringent vibration criteria for their functionality in normal operation conditions and their safety during an earthquake. Most previous investigations were based on simplified planar models of building structures, despite the fact that real ground motions and structures are always three‐dimensional. This paper hence presents a three‐dimensional analytical study of a hybrid platform on which high‐tech equipments are mounted for their vibration mitigation. The design methodology of the hybrid platform proposed in this study is based on dual‐level performance objectives for high‐tech equipments: safety against seismic hazard and functionality against traffic‐induced microvibration. The passive devices (represented by springs and viscous dampers) and the active actuators are designed, respectively, to meet vibration criteria corresponding to safety level and functionality level. A prototype three‐story building with high‐tech equipments installed on the second floor is selected in the case study to evaluate the effectiveness of the hybrid platform. The optimal location of the platform on the second building floor is determined during the design procedure in terms of the minimal H ₂ cost function of absolute velocity response. The simulation of the coupled actuator‐platform‐building system subjected to three‐dimensional ground motions indicates that the optimally designed hybrid platform can well achieve the dual target performance and effectively mitigate vibration at both ground motion levels. Copyright © 2009 John Wiley & Sons, Ltd.     

Screening‐Level Sensitivity Analysis for the Design of Pump‐and‐Treat Systems

Pump‐and‐treat systems can prevent the migration of groundwater contaminants and candidate systems are typically evaluated with groundwater models. Such models should be rigorously assessed to determine predictive capabilities and numerous tools and techniques for model assessment are available. While various assessment methodologies (e.g., model calibration, uncertainty analysis, and Bayesian inference) are well‐established for groundwater modeling, this paper calls attention to an alternative assessment technique known as screening‐level sensitivity analysis (SLSA). SLSA can quickly quantify first‐order (i.e., main effects) measures of parameter influence in connection with various model outputs. Subsequent comparisons of parameter influence with respect to calibration vs. prediction outputs can suggest gaps in model structure and/or data. Thus, while SLSA has received little attention in the context of groundwater modeling and remedial system design, it can nonetheless serve as a useful and computationally efficient tool for preliminary model assessment. To illustrate the use of SLSA in the context of designing groundwater remediation systems, four SLSA techniques were applied to a hypothetical, yet realistic, pump‐and‐treat case study to determine the relative influence of six hydraulic conductivity parameters. Considered methods were: Taguchi design‐of‐experiments (TDOE); Monte Carlo statistical independence (MCSI) tests; average composite scaled sensitivities (ACSS); and elementary effects sensitivity analysis (EESA). In terms of performance, the various methods identified the same parameters as being the most influential for a given simulation output. Furthermore, results indicate that the background hydraulic conductivity is important for predicting system performance, but calibration outputs are insensitive to this parameter (K_BK). The observed insensitivity is attributed to a nonphysical specified‐head boundary condition used in the model formulation which effectively “staples” head values located within the conductivity zone. Thus, potential strategies for improving model predictive capabilities include additional data collection targeting the K_BK parameter and/or revision of model structure to reduce the influence of the specified head boundary.     

Stochastic procedure for the simulation of synthetic main shock‐aftershock ground motion sequences

According to the current seismic codes, structures are designed to resist the first damaging earthquake during their service life. However, after a strong main shock, a structure may still face damaging aftershocks. The main shock‐aftershock sequence may result in major damage and eventually the collapse of a structure. Current studies on seismic hazard mainly focus on the modeling and simulation of main shocks. This paper proposes a 3‐step procedure to generate main shock‐aftershock sequences of pairs of horizontal components of a ground motion at a site of interest. The first step generates ground motions for the main shock using either a source‐based or site‐based model. The second step generates sequences of aftershocks' magnitudes, locations, and times of occurrence using either a fault‐based or seismicity‐based model. The third step simulates pairs of ground motion components using a new empirical model proposed in this paper. We develop prediction equations for the controlling parameters of a ground motion model, where the predictors are the site condition and the aftershock characteristics from the second step. The coefficients in the prediction equations and the correlation between the model parameters (of the 2 horizontal components of 1 record and of several records in 1 sequence) are estimated using a database of aftershock accelerograms. A backward stepwise deletion method is used to simplify the initial candidate prediction equations and avoid overfitting the data. The procedure, based on easily identifiable engineering parameters, is a useful tool to incorporate effects of aftershocks into seismic analysis and design.     

Improved risk‐targeted performance‐based seismic design of reinforced concrete frame structures

This paper presents a procedure for seismic design of reinforced concrete structures, in which performance objectives are formulated in terms of maximum accepted mean annual frequency (MAF) of exceedance, for multiple limit states. The procedure is explicitly probabilistic and uses Cornell's like closed‐form equations for the MAFs. A gradient‐based constrained optimization technique is used for obtaining values of structural design variables (members' section size and reinforcement) satisfying multiple objectives in terms of risk levels. The method is practically feasible even for real‐sized structures thanks to the adoption of adaptive equivalent linear models where element‐by‐element stiffness reduction is performed (2 linear analyses per intensity level). General geometric and capacity design constraints are duly accounted for. The procedure is applied to a 15‐storey plane frame building, and validation is conducted against results in terms of drift profiles and MAF of exceedance, obtained by multiple‐stripe analysis with records selected to match conditional spectra. Results show that the method is suitable for performance‐based seismic design of RC structures with explicit targets in terms of desired risk levels.     

Shake‐table tests of a reinforced concrete frame designed following modern codes: seismic performance and damage evaluation

This paper presents shake‐table tests conducted on a two‐fifths‐scale reinforced concrete frame representing a conventional construction design under current building code provisions in the Mediterranean area. The structure was subjected to a sequence of dynamic tests including free vibrations and four seismic simulations in which a historical ground motion record was scaled to levels of increasing intensity until collapse. Each seismic simulation was associated with a different level of seismic hazard, representing very frequent, frequent, rare and very rare earthquakes. The structure remained basically undamaged and within the inter‐story drift limits of the ‘immediate occupancy’ performance level for the very frequent and frequent earthquakes. For the rare earthquake, the specimen sustained significant damage with chord rotations of up to 28% of its ultimate capacity and approached the upper bound limit of inter‐story drift associated with ‘life safety’. The specimen collapsed at the beginning of the ‘very rare’ seismic simulation. Besides summarizing the experimental program, this paper evaluates the damage quantitatively at the global and local levels in terms of chord rotation and other damage indexes, together with the energy dissipation demands for each level of seismic hazard. Further, the ratios of column‐to‐beam moment capacity recommended by Eurocode 8 and ACI‐318 to guarantee the formation of a strong column‐weak beam mechanism are examined. Copyright © 2013 John Wiley & Sons, Ltd.     

Parameter identification for on‐line model updating in hybrid simulations using a gradient‐based method

To improve the efficiency of model fitting, parameter identification techniques have been actively investigated. Recently, the applications of parameter identification migrated from off‐line model fitting to on‐line model updating. The objective of this study is to develop a gradient‐based method for model updating to advance hybrid simulation also called hybrid test. A novel modification of the proposed method, which can reduce the number of design variables to improve the identification efficiency, is illustrated in detail. To investigate the model updating, simulated hybrid tests were conducted with a 5‐story steel frame equipped with buckling‐restrained braces (BRBs) utilized in the shaking table tests conducted in E‐Defense in Japan in 2009. The calibrated analytical model that was verified with the test results can serve as the reference model. In the simulated hybrid tests, the physical BRB substructure is numerically simulated by utilizing a truss element with the 2‐surface model identical to the part of the reference model. Such numerical verification allows simulation of measurement errors for investigation on the performance of the proposed method. Moreover, the feasibility of sharing the identified parameter values, which were obtained from the physical substructure responses, with the relevant numerical models is also verified with the artificial component responses derived from the physical experiments.