基于三角条块法的土体主、被动临界滑动场及土压力的计算

提出基于三角条块法的土体主、被动临界滑动场的数值模拟方法,并得到相应的土压力分布。首先,将土体划分为若干三角条块,引入带参数的条间力函数,建立土体极限平衡方程;然后,由最优控制理论得到一系列互不相交的临界滑动面,选定通过墙脚的临界滑动面,根据平衡方程求解格式,不断调整条间力函数中的参数,使土体满足力与力矩的平衡;最后,利用有限差分法,由墙面上各离散点的临界推力计算墙面上的土压力分布。数值计算结果表明,此方法适应性强,可计算任意倾斜角度的填土和墙体,且精度高、与现有解析解一致,可用于工程计算。     

Solution of random structural system subject to nonstationary excitation: transforming the equation with random coefficients to one with deterministic coefficients and random initial conditions

This paper deals with the lower order (first four) nonstationary statistical moments of the response of linear systems with random stiffness and random damping properties subject to random nonstationary excitation modeled as white noise multiplied by an envelope function. The method of analysis is based on a Markov approach using stochastic differential equations (SDE). The linear SDE with random coefficients subject to random excitation with deterministic initial conditions are transformed to an equivalent nonlinear SDE with deterministic coefficients and random initial conditions subject to random excitation. In this procedure, new SDE with random initial conditions, deterministic coefficients and zero forcing functions are introduced to represent the random variables. The joint statistical moments of the response are determined by considering an augmented dynamic system with state variables made up of the displacement and velocity vectors and the random variables of the structural system. The zero time-lag joint statistical moment equations for the augmented state vector are derived from the Itô differential formula. The statistical moment equations are ordinary nonlinear differential equations where hierarchy of moments appear. The hierarchy is closed by the cumulant neglect closure method applied at the fourth order statistical moment level. General formulation is given for multi-degree-of-freedom (MDOF) systems and the performance of the method in problems with nonstationary excitations and large variabilities is illustrated for a single-degree-of-freedom (SDOF) oscillator.     

高斯白噪声作用下刚性滑移结构稳态响应的概率密度函数

利用能量包线随机平均法，导出了库伦滑动摩擦刚性结构在高斯白噪声地面作用下响应的稳态联合概率密度函数的解析解，它与文献[2]的等效非线性系统法所得之解完全相同。与蒙特卡洛模拟结果对比表明，当库伦阻尼较小时，该解析解能够满意地预测研究系统响应的概率分布。     

工程结构地震破坏概率矩阵分析

本文提出了一种计算工程结构地震坡坏概率矩阵的方法，建立了地震地面运动模型和结构分析模型，对结构进行了随机地震反应分析，并获得了结构随机分应的统计量，进而采用双参数的结构破坏模型，给出了教育处结构地震破坏概率的表达式，利用此方法计算了一座按8度要求设计的钢筋混凝土框架型，给出了计算结构地震坡坏概率的表达式，利用此方法计算了一座按8度要求的钢筋混凝土框架结构的地震破坏概率矩阵，本文提出了方法可以在确定抗震设防标准和进行震害预测时采用。     

一种内蕴基线漂移校正的人工地震波 反应谱拟合方法

采用与传统频域法相反的思路,提出一种内蕴基线漂移校正且匹配目标反应谱的人工地震波直接拟合新方法。该方法首先构造非平稳位移时程表达式,进而求导得到速度、加速度时程表达式,由各表达式满足的归零条件确定相关的包络函数;其次,结合单自由度系统谐波响应的解析式,将匹配目标反应谱的拟合问题转化为关于待求幅值谱的非线性方程组;最后,结合高效的非线性方程迭代算法给出自然满足归零条件的加速度、速度及位移时程。分别以Rg1.60标准谱、某核岛楼层谱和《建筑抗震设计规范》中的设计谱为拟合算例验证所提方法的效率和精度。所提方法可成为人工地震波快速拟合的新途径。     

A new method for estimating BOD and the rate of biochemical oxidation based on modified Streeter-Phelps equations

An analytical solution of a modified system of Streeter-Phelps equations in the absence of aeration of the water body is given. Such solution is of interest in calculating the dynamics of biochemical oxidation of organic matter in a water body or a river channel covered by ice. The same solution can be used to describe processes taking place in a sealed flask in the determination of biochemical demand of the examined water sample over some period. The analytical solution of the modified system of Streeter-Phelps equations was used to derive formulas for calculating biochemical oxygen demand and the rate of biochemical oxidation.     

Moment analysis for compounds undergoing sequential chain reactions with first-order decay

Temporal and spatial moment analysis of one-dimensional equations governing fate and transport of parent compounds along with their transformation products is useful for parameter estimation of model parameters, and for understanding the average attributes of contaminant behavior. The objective of this paper is to present analytical expressions for the lower order moments of members in a sequential chain reaction, where members undergo first-order decay to produce the next member in the chain. Specifically, moments up to second order are discussed for the first two members. For the case of purely advective transport (Peclet number tending to infinity), temporal moment expressions are provided for more members of the chain. The sensitivity of temporal moments is examined with respect to Peclet number and transformation rates. Spatial moments are derived by two methods—one using Fourier transforms and another using moment generating differential equations. The behavior of lower order moments for the first few members of a chain can be elucidated from their mathematical expressions. However, expressions for higher order moments tend to be very complicated especially for members further down the chain.     

Generalized failure envelope for caisson foundations in cohesive soil: Static and dynamic loading

The response of massive caisson foundations to combined vertical (N), horizontal (Q) and moment (M) loading is investigated parametrically by a series of three-dimensional finite element analyses. The study considers foundations in cohesive soil, with due consideration to the caisson-soil contact interface conditions. The ultimate limit states are presented by failure envelopes in dimensionless and normalized forms and the effects of the embedment ratio, vertical load and interface friction on the bearing capacity are studied in detail. Particular emphasis is given on the physical and geometrical interpretation of the kinematic mechanisms that accompany failure, with respect to the loading ratio M/Q. Exploiting the numerical results, analytical expressions are derived for the capacities under pure horizontal, moment and vertical loading, for certain conditions. For the case of fully bonded interface conditions, comparison is given with upper bound limit equilibrium solutions based on Brinch Hansen theory for the ultimate lateral soil reaction. A generalized closed-form expression for the failure envelope in M–Q–N space is then proposed and validated for all cases examined. It is shown that the incremental displacement vector of the caisson at failure follows an associated flow rule, with respect to the envelope, irrespective of: (a) the caisson geometry, and (b) the interface conditions. A simplified geometrical explanation and physical interpretation of the associativity in M-Q load space is also provided. Finally, the derived failure envelope is validated against low (0.67 Hz) and high frequency (5 Hz) dynamic loading tests and the role of radiation damping on the response of the caisson at near failure conditions is unraveled.     

Characterizing the spatial variability of groundwater quality using the entropy theory: I. Synthetic data

This paper, the first in a series of two, applies the entropy (or information) theory to describe the spatial variability of synthetic data that can represent spatially correlated groundwater quality data. The application involves calculating information measures such as transinformation, the information transfer index and the correlation coefficient. These measures are calculated using discrete and analytical approaches. The discrete approach uses the contingency table and the analytical approach uses the normal probability density function. The discrete and analytical approaches are found to be in reasonable agreement. The analysis shows that transinformation is useful and comparable with correlation to characterize the spatial variability of the synthetic data set, which is correlated with distance. Copyright © 2004 John Wiley & Sons, Ltd.     

An improved response spectrum method for calculating seismic design response. Part 2: Non-classically damped structures

A response spectrum method which combines the analytical advantage of the mode acceleration formulation and the practical advantage of the mode displacement formulation is developed for seismic response calculation of non-classically damped structures. It reduces the error associated with the truncation of the high frequency modes without explicitly using them in the analysis. The method is especially effective for calculating the response of stiff structural systems and also for calculating the response quantities which are strongly affected by high frequency modes. Even with flexible structures, it is shown to provide more accurate response results than the results obtained with the mode displacement approach.     

Seismic analysis of structures with a fractional derivative model of viscoelastic dampers

Viscoelastic dampers are now among some of the preferred energy dissipation devices used for passive seismic response control. To evaluate the performance of structures installed with viscoelastic dampers, different analytical models have been used to characterize their dynamic force deformation characteristics. The fractional derivative models have received favorable attention as they can capture the frequency dependence of the material stiffness and damping properties observed in the tests very well. However, accurate analytical procedures are needed to calculate the response of structures with such damper models. This paper presents a modal analysis approach, similar to that used for the analysis of linear systems, for solving the equations of motion with fractional derivative terms for arbitrary forcing functions such as those caused by earthquake induced ground motions. The uncoupled modal equations still have fractional derivatives, but can be solved by numerical or analytical procedures. Both numerical and analytical procedures are formulated. These procedures are then used to calculate the dynamic response of a multi-degree of freedom shear beam structure excited by ground motions. Numerical results demonstrating the response reducing effect of viscoelastic dampers are also presented.     

AN ANALYTICAL SOLUTION FOR CALCULATING THE INITIATION OF SEDIMENT MOTION

This paper presents an analytical solution for calculating the initiation of sediment motion and the risk of river bed movement. It thus deals with a fundamental problem in sediment transport, for which no complete analytical solution has yet been found. The analytical solution presented here is based on forces acting on a single grain in state of initiation of sediment motion. The previous procedures for calculating the initiation of sediment motion are complemented by an innovative combination of optical surface measurement technology for determining geometrical parameters and their statistical derivation as well as a novel approach for determining the turbulence effects of velocity fluctuations. This two aspects and the comparison of the solution functions presented here with the well known data and functions of different authors mainly differ the presented solution model for calculating the initiation of sediment motion from previous approaches. The defined values of required geometrical parameters are based on hydraulically laboratory tests with spheres. With this limitations the derivated solution functions permit the calculation of the effective critical transport parameters of a single grain, the calculation of averaged critical parameters for describing the state of initiation of sediment motion on the river bed, the calculation of the probability density of the effective critical velocity as well as the calculation of the risk of river bed movement. The main advantage of the presented model is the closed analytical solution from the equilibrium of forces on a single grain to the solution functions describing the initiation of sediment motion.     

Rocking damage‐free steel column base with friction devices: design procedure and numerical evaluation

Earthquake‐resilient steel frames, such as self‐centering frames or frames with passive energy dissipation devices, have been extensively studied during the past decade, but little attention has been paid to their column bases. The paper presents a rocking damage‐free steel column base, which uses post‐tensioned high‐strength steel bars to control rocking behavior and friction devices to dissipate seismic energy. Contrary to conventional steel column bases, the rocking column base exhibits monotonic and cyclic moment–rotation behaviors that are easily described using simple analytical equations. Analytical equations are provided for different cases including structural limit states that involve yielding or loss of post‐tensioning in the post‐tensioned bars. A step‐by‐step design procedure is presented, which ensures damage‐free behavior, self‐centering capability, and adequate energy dissipation capacity for a predefined target rotation. A 3D nonlinear finite element (FE) model of the column base is developed in abaqus . The results of the FE simulations validate the accuracy of the moment–rotation analytical equations and demonstrate the efficiency of the design procedure. Moreover, a simplified model for the column base is developed in OpenSees . Comparisons among the OpenSees and abaqus models demonstrate the efficiency of the former and its adequacy to be used in nonlinear dynamic analysis. A prototype steel building is designed as a self‐centering moment‐resisting frame with conventional or rocking column bases. Nonlinear dynamic analyses show that the rocking column base fully protects the first story columns from yielding and eliminates the first story residual drift without any detrimental effect on peak interstory drifts. The study focuses on the 2D rocking motion and, thus, ignores 3D rocking effects such as biaxial bending deformations in the friction devices. The FE models, the analytical equations, and the design procedure will be updated and validated to cover 3D rocking motion effects after forthcoming experimental tests on the column base. Copyright © 2017 John Wiley & Sons, Ltd.     

An analytical solution for the probabilistic response of sdof non-linear random systems subjected to variable amplitude cyclic loading

Design for a specific ductile failure mode is assuming a rǒle of increasing importance for earthquake-resistant structures. This necessitates an accurate assessment of the distribution of overstrength in the structure, in order that the predefined failure mode can be realized. Consequently, the variability of the response for a given variability in the salient material properties, such as yield strength for steel structures, should be assessed and accounted for. In this paper an analytical method is proposed for the evaluation of the probability density function of the response of a single-degree-of-freedom hysteretic system with random parameters subject to a variable amplitude cyclic load history. A simple algorithm is derived which may be used to obtain the system response as a function of the system parameters. This response function may then be used to evaluate the displacement response probability density function when given the probability density function of the system parameters. Results derived from this procedure are verified against Monte Carlo simulation. It is shown that accurate response statistics are obtained at a fraction of the computing cost of simulation techniques.     

Divergence of solutions to perturbation-based advection-dispersion moment equations

Moment equation methods are popular and powerful tools for modeling transport processes in randomly heterogeneous porous media, but the application of these methods to advection-dispersion equations often leads to erroneous oscillations. Perturbative methods, required to close systems of moment equations, become inaccurate for large perturbations; however, little quantitative theory exists for determining when this occurs for advection-dispersion equations. We consider three different methods (asymptotic approximation, Eulerian truncation, and iterative solution) for closing and solving advection-dispersion moment equations describing transport in stratified porous media with random permeability. We obtain approximate analytical expressions for time above which the asymptotic approximation to the mean diverges, in particular quantifying the impact that dispersion has on delaying—but not eliminating—divergence. We demonstrate that Eulerian truncation and iterative solution methods do not eliminate divergent behavior either. Our divergence criteria provide a priori estimates that signal a warning to the practitioner of stochastic advection-dispersion equations to carefully consider whether to apply perturbative approaches.     

Vulnerability of ordinary moment resistant concrete frames

The effects of structural modeling (bar slip in lap splice), ground motion selection process (epsilon effect) and size of a structure (number of bays and stories) on the fragility of reinforced concrete ordinary moment resisting frames are investigated. An analytical model is developed to account for bar slip in lap splice, which exhibits good correlation with existing experimental data. Then, incremental dynamic analysis is used to derive the fragility curves for four model structures. The model structures simulate frames with a different number of bays and stories. Finally, the fragility curves are corrected for the epsilon effect. The results show that slip in the lap splice can significantly increase the failure probability, especially for the collapse prevention limit state. At the same time, it is found that spectral shape has a significant impact on the fragility curves. It is also found that accounting for or ignoring bar slip or epsilon effects increases the probability of failure for larger structures. This indicates an unconservative bias in the safety of larger structures.     

Order statistics of functionals of strong ground motion for a class of renewal processes

In this paper, probability distribution functions are derived for the order statistics of various functionals of strong ground motion at a site. These functionals can be: Modified Mercalli Intensity (MMI), peak ground acceleration (PGA), Fourier spectral amplitudes of acceleration, response spectrum amplitudes (spectral displacement, pseudo-spectral velocity and pseudo-spectral acceleration), and amplitudes of the peaks (local maxima and local minima) in the time historyof the response of SDOF and MDOF structures at the site. Three parameters of the response of a structure are considered: displacement, shear force and bending moment at each level (storey) of the structure. The earthquake sources contributing to the risk of ground motion at the site are a number of point, area or volume sources, each with defined frequency of occurence-magnitude relationship. The magnitudes of the possible events at these sources are discretized, and the occurrence of events of different magnitudes are assumed to be statistically independent. For each magnitude, it is assumed that the eartquakes occur in a Poissonian sequence or in a renewal process which is a generalization of the Poissonian. For these assumptions, the probability distribution functions are presented for the number of earthquakes, n, during which a given level of site or structural response is exceeded during the exposure time, and for the return period of the exceedances. For example, for single-degree- of-freedom: (SDOF) or multi-degree-of-freedom structures, (MDOF) n can be the number of earthquakes during which the response of a storey will exceed a given level at least m times(m = 1, 2, 3,…) during the exposure time. These probability distribution functions can be used to extend the concept of uniform probability functionals to more than one exceedance. A more important application is to generalize the uniform probability functionals method of site response (uniform probability Fourier or response spectra) to uniform probability envelopes of displacement, shears and bending moments of a given structure. The uniform probability envelopes can be for exceedance at least once during at least one earthquake, or, in general, for exceedance at least m times per earthquake (m = 1, 2,…) during at least n earthquakes. In other words, during at least n earthquakes at least m peaks in the response can be higher than the specified level. Such uniform probability envelopes can be used (1) to define new design guidelines for building codes based on cost-benefit analysis; (2) to construct more refined probability distribution functions for the damage and total economic losses caused by earthquakes; and (3) to develop planning and decision strategies on strengthening and retrofitting existing buildings.     

Horizontal dynamic response of a large-diameter pipe pile considering the second-order effect of axial force

The second-order effect of axial force on horizontal vibrating characteristics of a large-diameter pipe pile is theoretically investigated. Governing equations of the pile-soil system are established based on elastodynamics. Three-dimensional wave equations of soil are decoupled through differential transformation and variable separation. Consequently, expressions of soil displacements and horizontal resistances can be obtained. An analytical solution of the pile is derived based on continuity conditions between the pile and soil, subsequently from which expressions of the complex impedances are deduced. Analyses are carried out to examine the second-order effect of axial force on the horizontal vibrating behavior of the pipe pile. Some conclusions can be summarized as follows: stiffness and damping factors are decreased with the application of axial force on the pile head; distributions of the pile horizontal displacement and rotation angle are regenerated due to the second-order effect of the applied axial force; and redistributions of the bending moment and shearing force occur due to the second-order effect of the applied axial force.     

A review on stochastic differential equations for applications in hydrology

Fundamentals of the theory of stochastic calculus and stochastic differential equations (SDE's) which are finding increasing application in water resources engineering are reviewed. The basics of probability theory, mean square calculus and the Wiener, white Gaussian and compound Poisson processes are given in preparation for a discussion of the general Itô SDE with drift, diffusion and jump discontinuity terms driven by Gaussian white noise and compound Poissionian impulses. Also discussed are stochastic integration and the derivation of moment equations via the Itô differential rule. The lierature of SDE's is reviewed with an emphasis on the more accessible sources.     

Time-harmonic response of transversely isotropic multilayered half-space in a cylindrical coordinate system

With the aid of the analytical layer-element method, a comprehensive analytical derivation of the response of transversely isotropic multilayered half-space subjected to time-harmonic excitations is presented in a cylindrical coordinate system. Starting with the governing equations of motion and the constitutive equations of transversely isotropic elastic body, and based on the Fourier expansion, Hankel and Laplace integral transform, analytical layer-elements for a finite layer and a half-space are derived. Considering the continuity conditions on adjacent layers׳ interfaces and the boundary conditions, the global stiffness matrix equations for multilayered half-space are assembled and solved. Finally, some numerical examples are given to make a comparison with the existing solution and to demonstrate the influence of parameters on the dynamic response of the medium.