Recent variation in reach‐scale bankfull discharge in the Lower Yellow River

Bankfull discharge is a key parameter in the context of river engineering and geomorphology, as an indicator of flood discharge capacity in alluvial rivers, and varying in response to the incoming flow and sediment regimes. Bankfull channel dimensions have significantly adjusted along the Lower Yellow River (LYR) due to recent channel degradation, caused by the operation of the Xiaolangdi Reservoir, which has led to longitudinal variability in cross‐sectional bankfull discharges. Therefore, it is more representative to describe the flood discharge capacity of the LYR, using the concept of reach‐averaged bankfull discharge. Previous simple mean methods to estimate reach‐scale bankfull discharge cannot meet the condition of flow continuity or account for the effect of different spacing between two sections. In this study, a general method to calculate cross‐sectional bankfull discharge using the simulated stage‐discharge relation is outlined briefly, and an integrated method is then proposed for estimating reach‐scale bankfull discharge. The proposed method integrates a geometric mean based on the log‐transformation with a weighted average based on the spacing between two consecutive sections, which avoids the shortcomings of previous methods. The post‐flood reach‐scale bankfull discharges in three different channel‐pattern reaches of the LYR were estimated annually during the period from 1999 to 2010 using the proposed method, based on surveyed post‐flood profiles at 91 sedimentation sections and the measured hydrological data at seven hydrometric sections. The calculated results indicate that: (i) the estimated reach‐scale bankfull discharges can effectively represent the flood discharge capacity of different reaches, with their ranges of variation being less than those of typical cross‐sectional bankfull discharges; and (ii) the magnitude of the reach‐scale bankfull discharge in each reach can respond well to the accumulative effect of incoming flow and sediment conditions. Finally, empirical relationships for different reaches in the LYR were developed between the reach‐scale bankfull discharge and the previous four‐year average discharge and incoming sediment coefficient during flood seasons, with relatively high correlation coefficients between them being obtained, and the reach‐scale bankfull discharges in different reaches predicted by the delayed response model were also presented for a comparison. These relations for the prediction of reach‐scale bankfull discharges were validated using the cross‐sectional profiles and hydrological data measured in 2011. Copyright © 2013 John Wiley & Sons, Ltd.     

Towards error‐free hybrid simulation using mixed variables

Two procedures are developed and implemented in a hybrid simulation system (HSS) with the aim of enhancing the accuracy and reliability of the online, i.e. pseudo‐dynamic, test results. The first procedure aims at correcting the experimental systematic error in executing the displacement command signal. The error is calculated as the difference between command and feedback signals and correlated to the actuator velocity using the least‐squares method. A feed‐forward error compensation scheme is devised leading to a more accurate execution of the test. The second procedure employs mixed variables with mode switching between displacement and force controls. The newly derived force control algorithm is evaluated using a parametric study to assess its stability and accuracy. The implementation of the mixed variables procedure is designed to adopt force control for high stiffness states of the structural response and displacement control otherwise, where the resolution of the involved instruments may favour this type of mixed control. A simple pseudo‐dynamic experiment of steel cantilever members is used to validate the HSS. Moreover, two experiments as application examples for the two developed procedures are presented. The two experiments focus on the seismic response of (a) timber shear walls and (b) reinforced concrete frames with and without unreinforced masonry infill wall. Copyright © 2007 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.     

Validation of (not‐historical) large‐event near‐fault ground‐motion simulations for use in civil engineering applications

Ground‐motion simulations generated from physics‐based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance‐based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground‐motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four‐step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground‐Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near‐field ground motions (<10km) from large earthquakes (M_w7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near‐field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse‐like near‐field real records from multiple events and, when applicable, with NGA‐W2 GMPEs. The proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground‐motion simulations, and for building confidence in their use in engineering domains.     

Seismic response of a cut‐and‐cover tunnel isolated by polymer material

An effective method for the seismic retrofitting of the sidewalls of cut‐and‐cover tunnels has not yet been established. Thus, a new method of seismic retrofitting for cut‐and‐cover tunnels, called the ‘polymer isolation method’, is developed here. In this method, thin walls made of polymer materials, called ‘isolation walls’, are inserted between the ground and the sidewalls of a tunnel. We demonstrate the effectiveness of the proposed method in reducing the seismic response of tunnels by using some numerical simulations. It is found that the proposed method depends on the thickness of the soil cover over the tunnel and the ratio of the stiffness between the soil and the structure. Furthermore, a simple chart is proposed for convenience to represent the applicability of the polymer isolation method to the design of seismic retrofits for cut‐and‐cover tunnels. Although this particular chart has been obtained through a limited case, some other typical cases in which the chart can also be implemented for special conditions are addressed, and the applicability and the limitations of this chart are suggested for possible conditions of the ground and of tunnel structures. Copyright © 2012 John Wiley & Sons, Ltd.     

Compensation of actuator delay and dynamics for real‐time hybrid structural simulation

Compensation of delay and dynamic response of servo‐hydraulic actuators is critical for stability and accuracy of hybrid experimental and numerical simulations of seismic response of structures. In this study, current procedures for compensation of actuator delay are examined and improved procedures are proposed to minimize experimental errors. The new procedures require little or no a priori information about the behavior of the test specimen or the input excitation. First, a simple approach is introduced for rapid online estimation of system delay and actuator command gain, thus capturing the variability of system response through a simulation. Second, an extrapolation procedure for delay compensation, based on the same kinematics equations used in numerical integration procedures is examined. Simulations using the proposed procedures indicate a reduction in high‐frequency noise in force measurements that can minimize the excitation of high‐frequency modes. To further verify the effectiveness of the compensation procedures, the artificial energy added to a hybrid simulation as a result of actuator tracking errors is measured and used for demonstrating the improved accuracy in the simulations. Copyright © 2007 John Wiley & Sons, Ltd.     

System identification applied to long‐span cable‐supported bridges using seismic records

This paper presents the application of system identification (SI) to long‐span cable‐supported bridges using seismic records. The SI method is based on the System Realization using Information Matrix (SRIM) that utilizes correlations between base motions and bridge accelerations to identify coefficient matrices of a state‐space model. Numerical simulations using a benchmark cable‐stayed bridge demonstrate the advantages of this method in dealing with multiple‐input multiple‐output (MIMO) data from relatively short seismic records. Important issues related to the effects of sensor arrangement, measurement noise, input inclusion, and the types of input with respect to identification results are also investigated. The method is applied to identify modal parameters of the Yokohama Bay Bridge, Rainbow Bridge, and Tsurumi Fairway Bridge using the records from the 2004 Chuetsu‐Niigata earthquake. Comparison of modal parameters with the results of ambient vibration tests, forced vibration tests, and analytical models are presented together with discussions regarding the effects of earthquake excitation amplitude on global and local structural modes. Copyright © 2007 John Wiley & Sons, Ltd.     

Most frequent value‐based factor analysis of direct‐push logging data

An improved, iteratively re‐weighted factor analysis procedure is presented to interpret engineering geophysical sounding logs in shallow unsaturated sediments. We simultaneously process cone resistance, electric resistivity, and nuclear data acquired by direct‐push tools to give robust estimates of factor variables and water content in unconsolidated heterogeneous formations. The statistical procedure is based on the iterative re‐weighting of the deviations between the measured and calculated data using the most frequent value method famous for its robustness and high statistical efficiency. The iterative approach improves the result of factor analysis for not normally distributed data and extremely noisy measurements. By detecting a strong regression relation between one of the extracted factors and the fractional volume of water, we establish an independent method for water content estimation along the penetration hole. We verify the estimated values of water volume by using a highly over‐determined, quality‐checked interval inversion procedure. The multidimensional extension of the statistical method allows the estimation of water content distribution along both the vertical and the horizontal coordinates. Numerical tests using engineering geophysical sounding data measured in a Hungarian loessy–sandy formation demonstrate the feasibility of the most frequent value‐based factor analysis, which can be efficiently used for a more reliable hydrogeophysical characterisation of the unsaturated zone.     

A model‐based data‐driven dictionary learning for seismic data representation

Planar waves events recorded in a seismic array can be represented as lines in the Fourier domain. However, in the real world, seismic events usually have curvature or amplitude variability, which means that their Fourier transforms are no longer strictly linear but rather occupy conic regions of the Fourier domain that are narrow at low frequencies but broaden at high frequencies where the effect of curvature becomes more pronounced. One can consider these regions as localised “signal cones”. In this work, we consider a space–time variable signal cone to model the seismic data. The variability of the signal cone is obtained through scaling, slanting, and translation of the kernel for cone‐limited (C‐limited) functions (functions whose Fourier transform lives within a cone) or C‐Gaussian function (a multivariate function whose Fourier transform decays exponentially with respect to slowness and frequency), which constitutes our dictionary. We find a discrete number of scaling, slanting, and translation parameters from a continuum by optimally matching the data. This is a non‐linear optimisation problem, which we address by a fixed‐point method that utilises a variable projection method with ?₁ constraints on the linear parameters and bound constraints on the non‐linear parameters. We observe that slow decay and oscillatory behaviour of the kernel for C‐limited functions constitute bottlenecks for the optimisation problem, which we partially overcome by the C‐Gaussian function. We demonstrate our method through an interpolation example. We present the interpolation result using the estimated parameters obtained from the proposed method and compare it with those obtained using sparsity‐promoting curvelet decomposition, matching pursuit Fourier interpolation, and sparsity‐promoting plane‐wave decomposition methods.     

Kinematic characteristics and the influence of reference velocities of phase‐shift‐plus‐interpolation and extended‐split‐step‐Fourier migration methods

The phase‐shift‐plus‐interpolation and extended‐split‐step‐Fourier methods are wavefield‐continuation algorithms for seismic migration imaging. These two methods can be applied to regions with complex geological structures. Based on their unified separable formulas, we show that these two methods have the same kinematic characteristics by using the theory of pseudodifferential operators. Numerical tests on a Marmousi model demonstrate this conclusion. Another important aspect of these two methods is the selection of reference velocities and we explore the influence of the selection of reference velocities by comparing the geometric progression method and the statistical method. We show that the geometric progression method is simple but does not take into account the velocity distribution while the statistical approach is relatively complex but reflects the velocity distribution.     

Earthquake soil‐structure interaction of nuclear power plants,differences in response to 3‐D, 3 × 1‐D,and 1‐D excitations

In soil‐structure interaction modeling of systems subjected to earthquake motions, it is classically assumed that the incoming wave field, produced by an earthquake, is unidimensional and vertically propagating. This work explores the validity of this assumption by performing earthquake soil‐structure interaction modeling, including explicit modeling of sources, seismic wave propagation, site, and structure. The domain reduction method is used to couple seismic (near‐field) simulations with local soil‐structure interaction response. The response of a generic nuclear power plant model computed using full earthquake soil‐structure interaction simulations is compared with the current state‐of‐the‐art method of deconvolving in depth the (simulated) free‐field motions, recorded at the site of interest, and assuming that the earthquake wave field is spatially unidimensional. Results show that the 1‐D wave‐field assumption does not hold in general. It is shown that the way in which full 3‐D analysis results differ from those which assume a 1‐D wave field is dependent on fault‐to‐site geometry and motion frequency content. It is argued that this is especially important for certain classes of soil‐structure systems of which nuclear power plants subjected to near‐field earthquakes are an example.     

Out‐of‐plane seismic behaviour of rocking masonry walls

The evaluation of the out‐of‐plane behaviour of unreinforced walls is one of the most debated topics in the seismic assessment of existing masonry buildings. The discontinuous nature of masonry and its interaction with the remainder of the building make the dynamic modelling of out‐of‐plane response troublesome. In this paper, the results of a shaking table laboratory campaign on a tuff masonry, natural scale, U‐shaped assemblage (façade adjacent to transverse walls) are presented. The tests, excited by scaled natural accelerograms, replicate the behaviour of external walls in existing masonry buildings, from the beginning of rocking motion to overturning. Two approaches have been developed for modelling the out‐of‐plane seismic behaviour: the discrete element method and an SDOF analytic model. Both approaches are shown to be capable of reproducing the experimental behaviour in terms of maximum rotation and time history dynamic response. Finally, test results and numerical time history simulations have been compared with the Italian seismic code assessment procedures. Copyright © 2011 John Wiley & Sons, Ltd.     

A noniterative design procedure for supplemental brace–damper systems in single‐degree‐of‐freedom systems

In this paper, a method for designing supplemental brace–damper systems in single‐degree‐of‐freedom (SDOF) structures is presented. We include the effects of the supporting brace stiffness in the dynamic response by using a viscoelastic Maxwell model. On the basis of the study of an SDOF under ground excitation, we propose a noniterative design procedure for simultaneously specifying both the damper and the brace while assuring a desired structural performance. It is shown that to increase the damper size beyond the value delivered by the proposed criteria will not provide any improvement but actually worsen the structural response. The design method presented here shows excellent agreement with the FEMA 273 design approach but offers solutions closer to optimality. Copyright © 2013 John Wiley & Sons, Ltd.     

Free‐vibration analysis of a three‐dimensional soil–structure system

A procedure which involves a non‐linear eigenvalue problem and is based on the substructure method is proposed for the free‐vibration analysis of a soil–structure system. In this procedure, the structure is modelled by the standard finite element method, while the unbounded soil is modelled by the scaled boundary finite element method. The fundamental frequency, and the corresponding radiation damping ratio as well as the modal shape are obtained by using inverse iteration. The free vibration of a dam–foundation system, a hemispherical cavity and a hemispherical deposit are analysed in detail. The numerical results are compared with available results and are also verified by the Fourier transform of the impulsive response calculated in the time domain by the three‐dimensional soil–structure–wave interaction analysis procedure proposed in our previous paper. The fundamental frequency obtained by the present procedure is very close to that obtained by Touhei and Ohmachi, but the damping ratio and the imaginary part of modal shape are significantly different due to the different definition of damping ratio. This study shows that although the classical mode‐superposition method is not applicable to a soil–structure system due to the frequency dependence of the radiation damping, it is still of interest in earthquake engineering to evaluate the fundamental frequency and the corresponding radiation damping ratio of the soil–structure system. Copyright © 2001 John Wiley & Sons, Ltd.     

Compensation of time‐delay for control of civil engineering structures

Time‐delay is an important issue in structural control. Applications of unsynchronized control forces due to time‐delay may result in a degradation of the control performance and it may even render the controlled structures to be unstable. In this paper, a state‐of‐the‐art review for available methods of time‐delay compensation is presented. Then, five methods for the compensation of fixed time‐delay are presented and investigated for active control of civil engineering structures. These include the recursive response method, state‐augmented compensation method, controllability based stabilization method, the Smith predictor method and the Pade approximation method, all are applicable to any control algorithm to be used for controlled design. Numerical simulations have been conducted for MDOF building models equipped with an active control system to demonstrate the stability and control performance of these time‐delay compensation methods. Finally, the stability and performance of the phase shift method, that is well‐known in civil engineering applications, have also been critically evaluated through numerical simulations. Copyright © 2000 John Wiley & Sons, Ltd.     

Large‐scale real‐time hybrid simulation of a three‐story steel frame building with magneto‐rheological dampers

A series of large‐scale real‐time hybrid simulations (RTHSs) are conducted on a 0.6‐scale 3‐story steel frame building with magneto‐rheological (MR) dampers. The lateral force resisting system of the prototype building for the study consists of moment resisting frames and damped brace frames (DBFs). The experimental substructure for the RTHS is the DBF with the MR dampers, whereas the remaining structural components of the building including the moment resisting frame and gravity frames are modeled via a nonlinear analytical substructure. Performing RTHS with an experimental substructure that consists of the complete DBF enables the effects of member and connection component deformations on system and damper performance to be accurately accounted for. Data from these tests enable numerical simulation models to be calibrated, provide an understanding and validation of the in‐situ performance of MR dampers, and a means of experimentally validating performance‐based seismic design procedures for real structures. The details of the RTHS procedure are given, including the test setup, the integration algorithm, and actuator control. The results from a series of RTHS are presented that includes actuator control, damper behavior, and the structural response for different MR control laws. The use of the MR dampers is experimentally demonstrated to reduce the response of the structure to strong ground motions. Comparisons of the RTHS results are made with numerical simulations. Based on the results of the study, it is concluded that RTHS can be conducted on realistic structural systems with dampers to enable advancements in resilient earthquake resistant design to be achieved. Copyright © 2014 John Wiley & Sons, Ltd.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 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.     

Tangential‐projection algorithm for manifold representation in unidentifiable model updating problems

The problem of updating a structural model and its associated uncertainties by utilizing structural response data is addressed. In an identifiable case, the posterior probability density function (PDF) of the uncertain model parameters for given measured data can be approximated by a weighted sum of Gaussian distributions centered at a number of discrete optimal values of the parameters at which some positive measure‐of‐fit function is minimized. The present paper focuses on the problem of model updating in the general unidentifiable case for which certain simplifying assumptions available for identifiable cases are not valid. In this case, the PDF is distributed in the neighbourhood of an extended and usually highly complex manifold of the parameter space that cannot be calculated explicitly. The computational difficulties associated with calculating the highly complex posterior PDF are discussed and a new adaptive algorithm, referred to as the tangential‐projection (TP) algorithm, allowing for an efficient approximate representation of the above manifold and the posterior PDF is presented. Using this approximation, expressions for calculating the uncertain predictive response are established. A numerical example involving noisy data is presented to demonstrate the proposed method. Copyright © 2002 John Wiley & Sons, Ltd.     

High‐accuracy practical spline‐based 3D and 2D integral transformations in potential‐field geophysics

Potential, potential field and potential‐field gradient data are supplemental to each other for resolving sources of interest in both exploration and solid Earth studies. We propose flexible high‐accuracy practical techniques to perform 3D and 2D integral transformations from potential field components to potential and from potential‐field gradient components to potential field components in the space domain using cubic B‐splines. The spline techniques are applicable to either uniform or non‐uniform rectangular grids for the 3D case, and applicable to either regular or irregular grids for the 2D case. The spline‐based indefinite integrations can be computed at any point in the computational domain. In our synthetic 3D gravity and magnetic transformation examples, we show that the spline techniques are substantially more accurate than the Fourier transform techniques, and demonstrate that harmonicity is confirmed substantially better for the spline method than the Fourier transform method and that spline‐based integration and differentiation are invertible. The cost of the increase in accuracy is an increase in computing time. Our real data examples of 3D transformations show that the spline‐based results agree substantially better or better with the observed data than do the Fourier‐based results. The spline techniques would therefore be very useful for data quality control through comparisons of the computed and observed components. If certain desired components of the potential field or gradient data are not measured, they can be obtained using the spline‐based transformations as alternatives to the Fourier transform techniques.