A one‐dimensional model for simulating armouring and erosion on hillslopes: 2. Long term erosion and armouring predictions for two contrasting mine spoils

This paper investigates the dynamics of soil armouring as a result of fluvial erosion for a non‐cohesive sandy gravel spoil from the Ranger Mine, Australia, and a cohesive silt loam spoil from the Northparkes Mine, Australia, using a model for hillslope soil armouring. These long term predictions concentrate on the temporal and spatial changes of the spoil grading and erosion over 100–200 years for the flat cap regions (1–2%) and steep batter edges (10–30%) typically encountered on waste rock dumps. The existence of a significant rock fragment fraction in the Ranger spoil means that it armours readily, while Northparkes does not. For Ranger the waste rock showed reductions in (1) cumulative erosion of up to 81% from that obtained by extrapolating the initial erosion rate out 100 years and (2) the erosion/year by more than 10‐fold. For Northparkes reductions were less marked, with the maximum reduction in erosion/year being 37% after 200 years. For Ranger the reductions were greatest and fastest for intermediate gradient hillslopes. For the steepest hillslopes the armouring decreased because the flow shear stresses were large enough to mobilize all material in the armour layer. Model uncertainty was assessed with probabilistic confidence limits demonstrating that these erodibility reductions were statistically significant. A commonly used hillslope erosion model (sediment flux = β₁ discharge m₁ slope n₁) was fitted to these predictions. The erodibility, β₁, and m₁ decreased with time, which was consistent with our physical intuition about armouring. At Ranger the parameter m₁ asymptoted to 1·5–1·6 while at Northparkes it asymptoted to 1·2–1·3. At Ranger transient spatial trends in armouring led to a short term (50–200 years in the future) reduction in n₁, to below zero under certain circumstances, recovering to an asymptote of about 0·5–1. At Northparkes n₁ asymptoted to about 0·6, with no negative transients predicted. The m₁ and n₁ parameters predicted for Ranger were shown to be consistent with field data from a 10‐year‐old armoured hillslope and consistent with published relationships between erodibility and rock content for natural hillslopes. Copyright © 2007 John Wiley & Sons, Ltd.     

Strength,stiffness and cyclic deformation capacity of RC frames converted into walls by infilling with RC

In seismic retrofitting of concrete buildings, frame bays are converted into reinforced concrete (RC) walls by infilling the space between the frame members with RC of a thickness of not more than their width. The cyclic behavior of the resulting wall depends on the connection between the RC infill and the surrounding RC members. The paper uses the results from 56 cyclic tests on such composite walls to express their properties in terms of the geometry, the reinforcement and the connection. Properties addressed are: (a) the yield moment at the story base; (b) the secant-to-yield-point stiffness over the shear span of the wall in a story; (c) the deflection at flexural failure in cyclic loading; (d) the cyclic shear resistance, including a sliding shear failure mode. Separate models are given for squat walls failing in shear and for those where the top of the column shears-off. The proposals are modifications of models developed in the past for monolithic RC walls from several hundred cyclic tests; blind application of these latter models as though the walls were monolithic gives, in general, unsafe predictions. By contrast, the diagonal compression strut approach in ASCE41-06 is safe-sided, but gives unacceptably large prediction scatter.     

A probabilistic model of seismicity: Kamchatka earthquakes

The catalog of Kamchatka earthquakes is represented as a probability space of three objects {Ω, $ \tilde F $ \tilde F P}. Each earthquake is treated as an outcome ω _i in the space of elementary events Ω whose cardinality for the period under consideration is given by the number of events. In turn, ω _i is characterized by a system of random variables, viz., energy class k_i, latitude φ _i , longitude λ _i , and depth h _i . The time of an outcome has been eliminated from this system in this study. The random variables make up subsets in the set $ \tilde F $ \tilde F and are defined by multivariate distributions, either by the distribution function $ \tilde F $ \tilde F (φ, λ, h, k) or by the probability density f(φ, λ, h, k) based on the earthquake catalog in hand. The probabilities P are treated in the frequency interpretation. Taking the example of a recurrence relation (RR) written down in the form of a power law for probability density f(k), where the initial value of the distribution function f(k ₀) is the basic data [Bogdanov, 2006] rather than the seismic activity A ₀, we proceed to show that for different intervals of coordinates and time the distribution f _elim(k) of an earthquake catalog with the aftershocks eliminated is identical to the distribution f _full(k), which corresponds to the full catalog. It follows from our calculations that f ₀(k) takes on nearly identical numeral values for different initial values of energy class k ₀ (8 ≤ k ₀ ≤ 12) f(k ₀). The difference decreases with an increasing number of events. We put forward the hypothesis that the values of f(k ₀) tend to cluster around the value 2/3 as the number of events increases. The Kolmogorov test is used to test the hypothesis that statistical recurrence laws are consistent with the analytical form of the probabilistic RR based on a distribution function with the initial value f(k ₀) = 2/3. We discuss statistical distributions of earthquake hypocenters over depth and the epicenters over various areas for several periods     

基于简化拉压杆模型的钢筋混凝土框架异型节点抗剪性能研究

在以往试验研究的基础上,建立了框架异型节点等效核心区在剪、压复合作用下的计算模型,应用简化拉压杆模型方法对等效核心区进行了受力分析.研究表明:基于简化拉压杆模型方法计算得到异型节点的抗剪承载力与现行规范建议方法计算结果相当,但简化拉压杆模型方法有明确的力学计算模型,基于该模型计算得到的抗剪承载力与试验结果吻合较好.在此...     

Optimum strength ratio of hysteretic damper

This paper deals with steel moment frames combined with hysteretic dampers. Specifically, it examines the strength of dampers required to minimize damage to moment frames under earthquake loading. The important structural variables are the ratio of the dampers' shear strength to the maximum resistance, termed β in this study, and the relative stiffness between the damper and the main frame, termed k. The equivalent viscous damping ratio of the entire structural system is affected by the quantities β and k and the ductility ratio of the main frame. The optimum dampers' strength ratio (β_opt), which maximizes the equivalent viscous damping ratio, was formulated as a function of k, and it was found that β_opt is independent of the ductility ratio of the main frame. Earthquake response analysis confirmed that damage to the main frame can be significantly reduced by hysteretic dampers and minimized at the value of β_opt. © 1998 John Wiley & Sons, Ltd.     

In‐plane and out‐of‐plane behavior of confined masonry walls for various toothing and openings details and prediction of their strength and stiffness

Eight half‐scale brick masonry walls were tested to study two important aspects of confined masonry (CM) walls related to its seismic behavior under in‐plane and out‐of‐plane loads. Four solid wall specimens tested to investigate the role of type of interface between the masonry and tie‐columns, such as toothing varying from none to every course. The other four specimens with openings were tested to study the effectiveness of various strengthening options around opening to mitigate their negative influence. In the set of four walls, one wall was infilled frame while the other three were CM walls of different configurations. The experimental results were further used to determine the accuracy of various existing models in predicting the in‐plane response quantities of CM walls. Confined masonry walls maintained structural integrity even when severely damaged and performed much better than infill frames. No significant effect of toothing details was noticed although toothing at every brick course was preferred for better post‐peak response. For perforated walls, provision of vertical elements along with continuous horizontal bands around openings was more effective in improving the overall response. Several empirical and semi‐empirical equations are available to estimate the lateral strength and stiffness of CM walls, but those including the contribution of longitudinal reinforcement in tie‐columns provided better predictions. The available equations along with reduction factors proposed for infills could not provide good estimates of strength and stiffness for perforated CM walls. However, recently proposed relations correlating strength/stiffness with the degree of confinement provided reasonable predictions for all wall specimens. Copyright © 2016 John Wiley & Sons, Ltd.     

Direct use of PGV for estimating peak nonlinear oscillator displacements

A predictive model is presented for estimating the peak inelastic oscillator displacements (S_d,ie) from peak ground velocity (PGV). The proposed model accounts for the variation of S_d,ie for bilinear hysteretic behavior under constant ductility (µ) and normalized lateral strength ratio (R) associated with postyield stiffness ratios of α=0 and 5%. The regression coefficients are based on a ground‐motion database that contains dense‐to‐stiff soil site recordings at distances of up to 30 km from the causative fault. The moment magnitude ( M ) range of the database is 5.2? M ?7.6 and the ground motions do not exhibit pulse‐dominant signals. Confined to the limitations imposed by the ground‐motion database, the model can estimate S_d,ie by employing the PGV predictions obtained from the attenuation relationships (ground‐motion prediction equations). In this way, the influence of important seismological parameters can be incorporated to the variation of S_d,ie in a fairly rationale manner. This feature of the predictive model advocates its implementation in the probabilistic seismic hazard analysis that employs scalar ground‐motion intensity indices. Various case studies are presented to show the consistent estimations of S_d,ie by the proposed model. The error propagation in the S_d,ie estimations is also discussed when the proposed model is associated with attenuation relationships. Copyright © 2008 John Wiley & Sons, Ltd.     

CPT-based seismic stability assessment of a hazardous waste site

In areas of high seismicity, seismic stability often controls hazardous and solid waste landfill closure design. The undrained shear strength (S_u) of the waste mass is fundamental to seismic slope stability analyses. The value of S_u for hazardous waste fill is often difficult to characterize. The physical and chemical natures of the waste fill typically preclude laboratory testing of the materials. In certain cases, Cone Penetration Test (CPT) soundings can provide a viable technique for evaluation of S_u provided that the cone shear strength factor N_k can be established. If hazardous waste materials laboratory testing is not an option, N_k may be evaluated based upon results of non-intrusive in situ testing. This paper presents a case history of the seismic stability assessment of a hazardous waste site in which N_k was established from the results of non-intrusive Spectral Analysis of Surface Waves (SASW) soundings and empirical correlations to shear strength of soils. Generalization of the proposed methodology to other sites should be done with caution owing to variability among the parameters used in the analyses.     

Experimental verification of the fracture density and shear‐wave splitting relationship using synthetic silica cemented sandstones with a controlled fracture geometry

We present laboratory ultrasonic measurements of shear‐wave splitting from two synthetic silica cemented sandstones. The manufacturing process, which enabled silica cementation of quartz sand grains, was found to produce realistic sandstones of average porosity 29.7 ± 0.5% and average permeability 29.4 ± 11.3 mD. One sample was made with a regular distribution of aligned, penny‐shaped voids to simulate meso‐scale fractures in reservoir rocks, while the other was left blank. Ultrasonic shear waves were measured with a propagation direction of 90° to the coincident bedding plane and fracture normal. In the water saturated blank sample, shear‐wave splitting, the percentage velocity difference between the fast and slow shear waves, of <0.5% was measured due to the bedding planes (or layering) introduced during sample preparation. In the fractured sample, shear‐wave splitting (corrected for layering anisotropy) of 2.72 ± 0.58% for water, 2.80 ± 0.58% for air and 3.21 ± 0.58% for glycerin saturation at a net pressure of 40 MPa was measured. Analysis of X‐ray CT scan images was used to determine a fracture density of 0.0298 ± 0.077 in the fractured sample. This supports theoretical predictions that shear‐wave splitting (SWS) can be used as a good estimate for fracture density in porous rocks (i.e., SWS = 100ε_f, where ε_f is fracture density) regardless of pore fluid type, for wave propagation at 90° to the fracture normal.     

Experimental seismic behavior of squat shear walls with precast concrete hollow moulds

This study proposes an innovative precast shear wall system, called an EVE precast hollow shear wall structure (EVE-PHSW). Precast panels in EVE-PHSW are simultaneously precast with vertical and horizontal holes. Noncontact lap splices of rebars are used in vertical joints connecting adjacent precast panels for automated prefabrication and easy in situ erection. The seismic behavior of EVE walls was examined through a series of tests on six wall specimens with aspect ratios of 1.0∼1.3. Test results showed that EVE wall specimens with inside cast-in situ concrete achieved the desired “strong bending and weak shear” and failed in shear mode. Common main diagonal cracks and brittle shear failure in squat cast-in situ walls were prevented. Inside cast-in situ concrete could significantly improve the shear strength and stiffness of EVE walls. The details of boundary elements (cast-in situ or prefabricated) and vertical joints (contiguous or spaced) had little effect on the global behavior of EVE walls. Noncontact lap splices in vertical joints could enable EVE walls to exhibit stable load-carrying capacity through extensive deformations. Evaluation on design codes revealed that both JGJ 3-2010 and ACI 318-14 provide conservative estimation of shear strength of EVE walls, and EVE walls achieved shear strength reserves comparative to cast-in situ walls. The recommended effective stiffness for cast-in situ walls in ASCE 41–17 appeared to be appropriate for EVE walls.

     

Probabilistic SWE reanalysis as a generalization of deterministic SWE reconstruction techniques

Snow accumulation and melt is highly variable in space and time in complex mountainous environments. Therefore, it is necessary to provide high‐resolution spatially and temporally distributed estimates of sub‐basin snow water equivalent (SWE) to accurately predict the timing and magnitude of snowmelt runoff. In this study, we compare two reconstruction techniques (a commonly used deterministic reconstruction vs a probabilistic data assimilation framework). The methods retrospectively estimate SWE from a time series of remotely sensed maps of fractional snow‐covered area (FSCA). In testing both methods over the Tokopah watershed in the Sierra Nevada (California), the probabilistic reconstruction approach is shown to be a more robust generalization of the deterministic reconstruction. Under idealized conditions, both probabilistic and deterministic approaches perform reasonably well and yield similar results when compared with in situ verification data, whereas the probabilistic reconstruction was found to be in slightly better agreement with snow‐pit observations. More importantly, the probabilistic approach was found to be more robust: unaccounted for biases in solar radiation impacted the probabilistic SWE estimates less than the deterministic case (4% vs 7% errors for water year (WY)1997 and 0% vs 3% errors for WY1999); the probabilistic reconstruction was found to be less sensitive to the number of available observations (6% vs 10% errors in WY1997 and 13% vs 44% errors in WY1999 from the nominal cases when four fewer FSCA images were available). Finally, results from the probabilistic reconstruction approach, which requires precipitation inputs (unlike the deterministic approach), were found to be relatively robust to bias in prior precipitation estimates, where the nominal case mean estimates were recovered even when an underestimated prior precipitation was used. The additional robustness of the probabilistic SWE reconstruction technique should prove useful in future applications over larger basins and longer periods in mountainous terrain. Copyright © 2013 John Wiley & Sons, Ltd.     

Control strategy for variable damping element considering near‐future excitation influence

A control strategy is proposed for variable damping elements (VDEs) used together with auxiliary stiffness elements (ASEs) that compose a time‐varying non‐linear Maxwell (NMW) element, considering near‐future excitation influence. The strategy first composes a state equation for the structural dynamics and the mechanical balance in the NMW elements. Next, it establishes a cost function for estimating future responses by the weighted quadratic norms of the state vector, the controlled force and the VDEs' damping coefficients. Then, the Euler equations for the optimum values are introduced, and also approximated by the first‐order terms under the autoregressive (AR) model of excitation information. Thus, at each moment t_k, the strategy conducts the following steps: (1) identify the obtained seismic excitation information to an AR model, and convert it to a state equation; and (2) determine VDEs' damping coefficients under the initial conditions at t_k and the final state at t_k+L, using the first‐order approximation of the Euler equations. The control effects are examined by numerical experiments. Copyright © 2000 John Wiley & Sons, Ltd.     

Dispersions for the pushover‐based risk assessment of reinforced concrete frames and cantilever walls

The paper presents the results of an investigation into the dispersion values, expressed in terms of limit‐state spectral accelerations, which could be used for the pushover‐based risk assessment of low‐height to mid‐height reinforced concrete frames and cantilever walls. The results of an extensive parametric study of a portfolio of test structures indicated that the dispersion values due to record‐to‐record variability and modelling uncertainty (β_LS,RU) are within the range from 0.3 to 0.55 for the near collapse limit state, and between 0.35 and 0.60 for the collapse limit state. The dispersions β_LS,RU proposed for the code‐conforming and the majority of old (non code‐conforming) frames are in between these values. On the other hand, the dispersions proposed for the old frames with a soft storey and an invariant plastic mechanism, and for the code‐conforming cantilever walls, are at the lower and upper bounds of the presented values, respectively. The structural parameters that influence these dispersions were identified, and the influence of different ground motion sets, and of the models used for the calculation of the rotation capacities of the columns, on the calculated fragility parameters was examined and quantified. The proposed dispersion values were employed in a practice‐oriented pushover‐based method for the estimation of failure probability for eight selected examples. The pushover‐based risk assessment method, although extremely simple and economical when compared with more rigorous probabilistic methods, was able to predict seismic risk with reasonable accuracy, thus showing it to be a practical tool for engineers. Copyright © 2016 John Wiley & Sons, Ltd.     

Parameter identification and uncertainty analysis for simulating streamflow in a river basin of Eastern India

Uncertainty is inherent in modelling studies. However, the quantification of uncertainties associated with a model is a challenging task, and hence, such studies are somewhat limited. As distributed or semi‐distributed hydrological models are being increasingly used these days to simulate hydrological processes, it is vital that these models should be equipped with robust calibration and uncertainty analysis techniques. The goal of the present study was to calibrate and validate the Soil and Water Assessment Tool (SWAT) model for simulating streamflow in a river basin of Eastern India, and to evaluate the performance of salient optimization techniques in quantifying uncertainties. The SWAT model for the study basin was developed and calibrated using Parameter Solution (ParaSol), Sequential Uncertainty Fitting Algorithm (SUFI‐2) and Generalized Likelihood Uncertainty Estimation (GLUE) optimization techniques. The daily observed streamflow data from 1998 to 2003 were used for model calibration, and those for 2004–2005 were used for model validation. Modelling results indicated that all the three techniques invariably yield better results for the monthly time step than for the daily time step during both calibration and validation. The model performances for the daily streamflow simulation using ParaSol and SUFI‐2 during calibration are reasonably good with a Nash–Sutcliffe efficiency and mean absolute error (MAE) of 0.88 and 9.70 m³/s for ParaSol, and 0.86 and 10.07 m³/s for SUFI‐2, respectively. The simulation results of GLUE revealed that the model simulates daily streamflow during calibration with the highest accuracy in the case of GLUE (R² = 0.88, MAE = 9.56 m³/s and root mean square error = 19.70 m³/s). The results of uncertainty analyses by SUFI‐2 and GLUE were compared in terms of parameter uncertainty. It was found that SUFI‐2 is capable of estimating uncertainties in complex hydrological models like SWAT, but it warrants sound knowledge of the parameters and their effects on the model output. On the other hand, GLUE predicts more reliable uncertainty ranges (R‐factor = 0.52 for daily calibration and 0.48 for validation) compared to SUFI‐2 (R‐factor = 0.59 for daily calibration and 0.55 for validation), though it is computationally demanding. Although both SUFI‐2 and GLUE appear to be promising techniques for the uncertainty analysis of modelling results, more and more studies in this direction are required under varying agro‐climatic conditions for assessing their generic capability. Copyright © 2015 John Wiley & Sons, Ltd.     

Erodibility of cohesive streambeds in the loess area of the midwestern USA

Excess stress parameters, critical shear stress (τ_c) and erodibility coefficient (k_d), for degrading channels in the loess areas of the midwestern USA are presented based on in situ jet‐testing measurements. Critical shear stress and k_d are used to define the erosion resistance of the streambed. The jet‐testing apparatus applies hydraulic stresses to the bed and the resulting scour due to the impinging jet is related to the excess stress parameters. Streams tested were primarily silt‐bedded in texture with low densities, which is typical of loess soils. Results indicate that there is a wide variation in the erosion resistance of streambeds, spanning six orders of magnitude for τ_c and four orders of magnitude for k_d. Erosion resistance was observed to vary within a streambed, from streambed to streambed, and from region to region. An example of the diversity of materials within a river system is the Yalobusha River Basin in Mississippi. The median value of τ_c for the two primary bed materials, Naheola and Porters Creek Clay Formations, was 1·31 and 256 Pa, respectively. Streambeds composed of the Naheola Formation are readily eroded over the entire range of shear stresses, whereas only the deepest flows generate boundary stresses great enough to erode streambeds composed of the Porters Creek Clay Formation. Therefore, assessing material resistance and location is essential in classifying and modelling streambed erosion processes of these streams.     

The null distribution of sample serial correlation coefficient

 The null distribution of the lag-k sample serial correlation coefficient (r _k , k=1,2,3) was investigated by Monte Carlo simulation. For a time series with normal, exponential, Pearson 3, EV1 (Gumbel), or generalized Pareto (GP) distribution type, the null distribution of its r _k can be approximated by the normal distribution with mean −1/(n−k) and variance 1/(n−1). But for a time series with the lognormal, EV2 or EV3 (Weibull) distribution type, the null distribution of r _k is skewed distributed. In such cases, a simulation technique is suggested to construct percentile confidence intervals at a given significance level.     

Experimental and analytical studies on earthquake resisting behaviour of confined concrete block masonry structures

To improve the seismic performance of masonry structures, confined masonry that improves the seismic resistance of masonry structures by the confining effect of surrounding bond beams and tie columns is constructed. This study investigated the earthquake resisting behaviour of confined masonry structures that are being studied and constructed in China. The structural system consists of unreinforced block masonry walls with surrounding reinforced concrete bond beams and tie columns. The characteristics of the structure include: (1) damage to blocks is reduced and brittle failure is avoided by the comparatively lower strength of the joint mortar than that of the blocks, (2) the masonry walls and surrounding reinforced concrete bond beams and tie columns are securely jointed by the shear keys of the tie columns. In this study, wall specimens made of concrete blocks were tested under a cyclic lateral load and simulated by a rigid body spring model that models non‐linear behaviour by rigid bodies and boundary springs. The results of studies outline the resisting mechanism, indicating that a rigid body spring model is considered appropriate for analysing this type of structure. Copyright © 2006 John Wiley & Sons, Ltd.     

A numerical model of the wave boundary layer

A numerical model of deep, uniform, oscillatory, rough-turbulent boundary-layer flow is described. The model is based upon the governing horizontal momentum equation and a closure scheme involving the turbulent-energy equation and various turbulence-scaling laws. Finite difference solutions of these equations are obtained for a range of values of the ‘relative roughness’ (A₀/k_s), whereA₀ is the excursion amplitude of the water particles in the free-stream flow andk_s is the ‘equivalent bed roughness’. Typical vertical profiles of horizontal velocity, turbulence energy and eddy viscosity, and time-series of the bed shear stress are presented. The model results are then used to determine the wave drag coefficient, boundary-layer thickness and phase lead of the bed shear stress over the free-stream velocity, each as a function ofA₀/k_s. These results are shown to be in generally good agreement with previous experimental and theoretical results. Finally, the model is used to test for the existence of a universal velocity distribution for uniform oscillatory (sinusoidal) rough-turbulent flow. The ‘law of the wall’ and the ‘defect law’ proposed by Jonsson (1980, Ocean Engineering, 7, 109–152) are well supported by the model, and the existence of a logarithmic ‘velocity overlap layer’ in which both of these laws are valid is demonstrated forA_o/k_s30.     

Scaling of ductility and damage‐based strength reduction factors for horizontal motions

The conventional approach of obtaining the inelastic response spectra for the aseismic design of structures involves the reduction of elastic spectra via response modification factors. A response modification factor is usually taken as a product of (i) strength factor, R_S, (ii) ductility factor, R_μ, and (iii) redundancy factor, R_R. Ductility factor, also known as strength reduction factor (SRF), is considered to primarily depend on the initial time period of the single‐degree‐of‐freedom (SDOF) oscillator and the displacement ductility demand ratio for the ground motion. This study proposes a preliminary scaling model for estimating the SRFs of horizontal ground motions in terms of earthquake magnitude, strong motion duration and predominant period of the ground motion, geological site conditions, and ductility demand ratio, with a given level of confidence. The earlier models have not considered the simultaneous dependence of the SRFs on various governing parameters. Since the ductility demand ratio is not a complete measure of the cumulative damage in the structure during the earthquake‐induced vibrations, the existing definition of the SRF is sought to be modified with the introduction of damage‐based SRF (in place of ductility‐based SRF). A parallel scaling model has been proposed for estimating the damage‐based SRFs. This model considers damage and ductility supply ratio as parameters instead of ductility demand ratio. Through a parametric study on ductility‐based SRFs, it has been shown that the hitherto assumed insensitivity of earthquake magnitude and strong motion duration may not be always justified and that the initial time period of the oscillator plays an important role in the dependence of SRF on these parameters. Further, the damage‐based SRFs are found to show similar parametric dependence as observed in the case of the ductility‐based SRFs. Copyright © 2000 John Wiley & Sons, Ltd.