首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
传统基于瑞典法的边坡加固抗震性能参数分析,通过建立强度参数与临界滑动面对应关系完成推导,由于因荷载过大,边坡建筑物抗震性能降低,准确率不高等问题,提出新的边坡加固抗震性能参数分析方法,对边坡加固的抗震性能参数进行分析,更好地通过加固提高其抗震性能。运用新型边坡加固抗震性能参数分析方法,采用基于极向条分的极限分析上限法和非线性Mohr-Coulomb破坏准则,分析边坡加固措施对边坡抗震性能的影响,通过对比有、无加固措施状态下不同能量功率的运算确定边坡抗震性能,利用MATLAB软件获取边坡加固抗震性能参数的最优解。经实验证明,边坡潜在破坏范围受边坡坡顶受荷和边坡加固情况的作用较大,地震荷载系数及非线性系数越大,土体强度非线性越小,边坡抗震性能受地震荷载的作用也越大。边坡抗滑桩最优抗震支护方位位于X_F/L_x=0.7处,在地震荷载从0上升至0.2的情况下,非线性系数为1.2、1.4、1.6以及2.0时的边坡抗震性能分别下降了40.1%、46.8%、57.5%以及61.5%。新型的抗震性能参数分析方法有效地提高了结果准确率,对边坡加固的抗震性能能达到准确分析。  相似文献   

2.
A method is proposed for evaluating rock hydraulic conductivity (k x;y , LT?1) in grid model blocks to be used in finite-difference simulation of mass exchange in gypcified fractured-porous rocks containing groundwater. The values of k x and k y were determined taking into account the dominating form of transport (advection or transverse dispersion), contributing most to gypsum dissolution. The method was used to simulate groundwater flow in a vertical section of gypsified rocks in the base of the dam at the Lower-Kafirniganskii Hydropower System. The changes in rock permeability in the base of the dam were used to evaluate whether there is a need to take seepage-control measures (cement-grout curtain, blanket, and a gypsum fill over the blanket). A leaching zone was found to form near the curtain and blanket during the simulation period (100 years) under background seepage conditions. In this period, the total water flow in the lower pool can increase by no more than 10–15%. Under such conditions, the “seepage-control” measures can enhance the reliability of the dam at the Lower-Kafirniganskii Hydropower System.  相似文献   

3.
 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.  相似文献   

4.
Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold‐and‐thrust belts) and in subsalt exploration. Here, we introduce a methodology for P‐wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters ε and δ and linearly varying symmetry‐direction velocity VP0 controlled by the vertical (kz) and lateral (kx) gradients. Since determination of tilt from P‐wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity VP0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients kz and kx and the anisotropy parameters ε and δ in the layer‐stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry‐axis direction is fixed and VP0 is known, the parameters kz, kx, ε and δ can be resolved from reflection data. It should be emphasized that estimation of ε in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas.  相似文献   

5.
Starting from the basic erosion principles, an upland soil erosion model to predict soil loss by overland flow from individual storms on forested hillslopes can be derived in the form where Qs is total soil loss for a storm event, n is roughness coefficient, x is down slope distance, Kf is soil erodibility factor, S is slope, α is slope exponent and Q is runoff. Values of n and α are to be determined for different environments and are 0·58 and 2·1 for a mixed pine forest ecosystem. A significant correlation (r = 0·933, n = 96) fits between the observed and predicted values using this expression, and the model fitting is good.  相似文献   

6.
A one-dimensional non-linear diffusion wave equation is derived from the Saint Venant equations with neglect of the inertia terms. This non-linear equation has no general analytical solution. Numerical schemes are therefore employed to discretize the space and time axes and convert the differential equation to difference form. In this study, the mixing cell method is used to convert the diffusion wave equation to difference form, in which the difference term can be eliminated by selecting an optimal space step size Δx when time step size Δt is given. When the time step size Δt→0, the space step size Δx=Q/(2S0BC]k) where Q is discharge, S0 is bed slope, B is channel width and Ck is kinematic wave celerity, which is the same as the characteristic length proposed by Kalinin and Milyukov. The results of application to two cases show that the mixing cell and linear channel flow routing methods produce hydrographs that are in agreement with the observed flood hydrographs. © 1997 John Wiley & Sons, Ltd.  相似文献   

7.
Slope frequency distributions are computed from Digital Elevation Models (DEMs) of 18 areas from five different geographic and physiographic regions. Over 160,000 slope estimates are used in the definition of the slope frequency distribution for each area. Analysis of several transformations of the slope data indicates that no single transformation is capable of normalizing all slope distributions. This is, in part, due to the arbitrary placement of the boundaries of the map sheets used in the data collection. The square root of sine, however, provides better results than sine, log-tangent, or several other transformations. The degree of (positive) skewness of gradients in degrees correlates with the proportion of slopes below 2° or 5° (rs = + 0·93), as does the kurtosis (rs = + 0·91).  相似文献   

8.
Three active earth-slide slopes of Tertiary mudstone were investigated: Slope 1 has an angle of 17 ?4°, Slope 2 of 12.9°, and Slope 3 of 11.6°. Infinite slope analysis indicated that the instability of these three slopes can be well explained by using the residual strength parameters of earth-slide soils near the sliding surface in conjunction with the highest ground water table. The residual angle of shearing resistance, ø'r, plays an important role in the determination of slope angle because it differs greatly among slopes. Mineralogical studies and X-ray diffraction analyses were performed for the clay minerals included in the slope material. The results showed that illite and chlorite were found in Slope 1, and also in Slopes 2 and 3 together with interstratified illite/montmorillonite and montmorillonite. The degree of weathering is progressive in order of Slopes 1, 2 and 3. The alteration of clay minerals by weathering causes the reduction in ø'r -values, i.e., 19.4° in Slope 1 (steep, less weathered) and 12.1-9.2° in Slopes 2 and 3 (gentle, much more weathered). This result indicates that the degree of weathering has a great influence on the value of ø'r, which in turn determines the slope angle.  相似文献   

9.
M. Barrios  F. Francés 《水文研究》2012,26(7):1022-1033
Nonlinear dynamics and spatial variability in hydrological systems make the formulation of scaling theories difficult. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter Hu and the saturated hydraulic conductivity parameter ks from microscale (sub‐grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of Hu and ks was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre‐established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the nonlinearity of the system from microscale to the mesoscale via non‐stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

10.
Abstract

The process of wave steepening in Long's model of steady, two-dimensional stably stratified flow over orography is examined. Under conditions of the long-wave approximation, and constant values of the background static stability and basic flow, Long's equation is cast into the form of a nonlinear advection equation. Spectral properties of this latter equation, which could be useful for the interpretation of data analyses under mountain wave conditions, are presented. The principal features, that apply at the onset of convective instability (density constant with height), are:

i) a power spectrum for available potential energy that exhibits a minus eight-thirds decay, in terms of the vertical wavenumber k z -;

ii) a rate of energy transfer across the spectrum that is inversely proportional to the wavenumber for large k z -;

iii) an equipartition between the kinetic energy of the horizontal motion and the available potential energy, under the longwave approximation, although all the disturbance energy is kinetic at the point where convective instability is initiated. It is also shown that features i) and ii) apply to more general conditions that are appropriate to Long's model, not just the long-wave approximation. Application to fully turbulent flow or to conditions at the onset of shearing instability are not considered to be warranted, since the development only applies to conditions at the onset of convective instability.  相似文献   

11.
The temporal evolution of simple landforms such as cinder cones by nonlinear diffusive processes is studied through the use of a new 2D numerical model using well‐established and accurate numerical mathematics and high‐resolution digital elevation models (DEMs). Extending 1D (profile) nonlinear diffusion analyses used in cinder cone, hillslope and fault scarp evolution studies, we have implemented a 2D numerical model with a spatially and temporally varying sediment transport rate coefficient scaled nonlinearly by the ratio of local slope to critical slope. The high accuracy and efficient numerical implementation are documented in the paper and the MATLAB toolkit developed is used to solve for the developmentof an initial 2D cone form. First, we examine the nonlinear transport rule and suggest a refinement that accounts explicitly for flux at threshold slopes. We find that the maximum diffusion (necessarily introduced in the numerical model to avoid infinite rates) at the critical slope controls the final morphology, especially approaching steady state. Secondly, solving the landscape evolution problem in 2D enables a natural accounting for sediment flux convergence or divergence in the profile. Thirdly, the boundary behavior of a given landscape element controls much of what happens in that domain and so we allow for arbitrary flux magnitude or elevation boundary conditions. Fourthly, landscapes are heterogeneous in their surface cover and so we allow for spatially and temporally varying transport rate k and we permit an arbitrary vertical displacement field within the model domain. To test the new formulation for the nonlinear term, the effect of variable diffusivity k and the numerical schemes implemented, we apply the model to cinder cones built on the flanks of Mount Etna in 2001 and 2002–2003. We explore the effects of DEM resolution with data from the 2001 cone and the utility of spatially variable diffusivity to explain the variation in erosion measured by differencing repeat light detection and ranging (LIDAR) surveys gathered in 2004 and 2007 over the 2002–2003 cone complex. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

12.
The catalytic destruction of stratospheric ozone by the oxides of nitrogen is believed to be an important part of the global ozone balance. The lack of sufficient measurements of NO x concentrations has impeded efforts to quantify this process. Recent measurements of stratospheric nitrogen dioxide from ground-based stations as well as aircraft and balloons have provided a first approximation to a global distribution of NO2 vertical columns at sunset. These observed vertical columns have been translated into time-dependent vertical NO2 profiles by means of a one-dimensional atmospheric photochemical model. Using recent observations of air temperature and ozone along with this information, the independent instantaneous (one second) rates of ozone production from oxygen photolysis P(O3), of ozone destruction from pure oxygen species (Chapman reactions) L(O x ), and of ozone destruction by nitrogen oxides L(NO x ) were estimated over the three-dimensional atmosphere. These quantities are displayed as zonal average contour maps, summed over various latitude zones, summed over various altitude bands, and integrated globally between 15 and 45 km. Although the global summation between 15 and 45 km by no means tells the complete story, these numbers are of some interest, and the relative values are: P(O3), 100; L(O x ), 15; L(NO x ), 45±15. It is to be emphasized that this relative NO x contribution to the integrated ozone balance is not a measure of the sensitivity of ozone to possible perturbations of stratospheric NO x ; recent model results must be examined for current estimates of this sensitivity.  相似文献   

13.
We look at the large-scale dynamo properties of spatially periodic, time dependent, helical 2D flows of the form u(x, t)?=?(? y ?ψ?(x, y, t), ?? x ?ψ?(x, y, t), ?ψ (x, y, t). These flows act as kinematic fast dynamos and are able to generate a mean magnetic field uniform and constant in the xy-plane but whose direction varies periodically along z with wavenumber k. Using Mean Field Electrodynamics, the generation mechanism can be understood in terms of a k-dependent α-effect, which depends on the magnetic Reynolds number, R m . We calculate this effect for different motions and investigate how its limit as k?→?0 depends on R m and on the properties of the flows such as their spatial structure or correlation time. This work generalises earlier studies based on 2D steady flows to motions with time dependence.  相似文献   

14.
Abstract

The influence of mesoscale topography on the baroclinic instability of a two-layer model of the open ocean is considered. For westward velocities in the top layer (U), and for a sinusoidal topography independent of x or longitude (a cross-stream topography), the critical value of U (Uc ) leading to instability is the same as when there is no topography. The wavelength of the unstable perturbation corresponding to U c is shortened. For a given wavevector (k) of the perturbation the system becomes stable (as also in the absence of topography) for large values of |U|. The minimum value of the shear leading to stability is, however, significantly reduced by the topography.

For sufficiently large values of the height of the topographic features, instabilities appear which are localized within a narrow range of the shear. These instabilities are studied for a topography that depends both on x and y.

For a cross-stream topography the growth rates are somewhat smaller than those without topography and they depend only weakly on ky . For the topographies considered here which depend both on x and y, perturbations with different values of ky can again have roughly the same growth rate.

In the case of stable oscillations, variations in the eddy energy with very long periods are made possible by the coexistence of topographic modes with closely lying periods.  相似文献   

15.
The time required at a field site to obtain a few measurements of saturated hydraulic conductivity (Ks) will allow for many measurements of soil air permeability (ka). This study investigates if ka measured in situ (ka, in situ) can be a substitute for measurement of Ks in relation to infiltration and surface runoff modelling. Measurements of ka, in situ were carried out in two small agricultural catchments. A spatial correlation of the log‐transformed values existed having a range of approximately 100 m. A predictive relationship between Ks and ka measured on 100‐cm3 soil samples in the laboratory was derived for one of the field slopes and showed good agreement with an earlier suggested predictive Ks–ka relationship. In situ measurements of Ks and ka suggested that the predictive relationships also could be used at larger scale. The Ks–ka relationships together with the ka, in situ data were applied in a distributed surface runoff (DSR) model, simulating a high‐intensity rainfall event. The DSR simulation results were highly dependent on whether the geometric average of ka, in situ or kriged values of ka, in situ was used as model input. When increasing the resolution of Ks in the DSR model, a limit of 30–40 m was found for both field slopes. Below this limit, the simulated runoff and hydrograph peaks were independent of resolution scale. If only a few randomly chosen values of Ks were used to represent the spatial variation within the field slope, very large deviations in repeated DSR simulation results were obtained, both with respect to peak height and hydrograph shape. In contrast, when using many predicted Ks values based on a Ks–ka relationship and measured ka, in situ data, the DSR model generally captured the correct hydrograph shape although simulations were sensitive to the chosen Ks–ka relationship. As massive measurement efforts normally will be required to obtain a satisfactory representation of the spatial variability in Ks, the use of ka, in situ to assess spatial variability in Ks appears a promising alternative. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

16.
In order to reconcile the larger scatter and avoid the biased estimate from deterministic predictions for the shear strength of reinforced concrete (RC) squat structural walls, a probabilistic shear strength model is developed in this paper based on the strut‐and‐tie model and the generalized likelihood uncertainty estimation (GLUE) method. The strut‐and‐tie model is used to derive an appropriate function form for the probabilistic shear strength model, where four unknown model parameters (e.g. k1, k2, k3 and k4) are defined carefully to guarantee them having a clear physical‐based meaning so that the corresponding prior distribution ranges can be specified reasonably. Then, the GLUE method is adopted to estimate the posterior cumulative distribution of k1, k2, k3 and k4 with an available experimental database. Furthermore, to demonstrate the stability of the estimated posterior cumulative distribution, the sensitivity of three major aspects in GLUE method is investigated. Finally, based on the estimated cumulative distribution of k1, k2, k3 and k4, the developed probabilistic shear strength model is simplified as a mean prediction model and a standard deviation prediction model for facilitate using in engineering practice. Therefore, with the developed probabilistic shear strength model, not only can the squat structural walls be designed in confidence, but the accuracy of those deterministic predictions can be evaluated in a probabilistic manner. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

17.
Song‐Bae Kim 《水文研究》2005,19(20):4069-4079
  相似文献   

18.
  相似文献   

19.
Fermat's variational principle states that the signal propagates from point S to R along a curve which renders Fermat's functional (l) stationary. Fermat's functional (l) depends on curves l which connect points S and R, and represents the travel times from S to R along l. In seismology, it is mostly expressed by the integral (l) = (x k,x k ')du, taken along curve l, where (x k,x k ') is the relevant Lagrangian, x k are coordinates, u is a parameter used to specify the position of points along l, and x k ' = dx k÷du. If Lagrangian (x k,x k ') is a homogeneous function of the first degree in x k ', Fermat's principle is valid for arbitrary monotonic parameter u. We than speak of the first-degree Lagrangian (1)(x k,x k '). It is shown that the conventional Legendre transform cannot be applied to the first-degree Lagrangian (1)(x k,x k ') to derive the relevant Hamiltonian (1)(x k,p k), and Hamiltonian ray equations. The reason is that the Hessian determinant of the transform vanishes identically for first-degree Lagrangians (1)(x k,x k '). The Lagrangians must be modified so that the Hessian determinant is different from zero. A modification to overcome this difficulty is proposed in this article, and is based on second-degree Lagrangians (2). Parameter u along the curves is taken to correspond to travel time , and the second-degree Lagrangian (2)(x k, k ) is then introduced by the relation (2)(x k, k ) = [(1)(x k, k )]2, with k = dx k÷d. The second-degree Lagrangian (2)(x k, k ) yields the same Euler/Lagrange equations for rays as the first-degree Lagrangian (1)(x k, k ). The relevant Hessian determinant, however, does not vanish identically. Consequently, the Legendre transform can then be used to compute Hamiltonian (2)(x k,p k) from Lagrangian (2)(x k, k ), and vice versa, and the Hamiltonian canonical equations can be derived from the Euler-Lagrange equations. Both (2)(x k, k ) and (2)(x k,p k) can be expressed in terms of the wave propagation metric tensor g ij(x k, k ), which depends not only on position x k, but also on the direction of vector k . It is defined in a Finsler space, in which the distance is measured by the travel time. It is shown that the standard form of the Hamiltonian, derived from the elastodynamic equation and representing the eikonal equation, which has been broadly used in the seismic ray method, corresponds to the second-degree Lagrangian (2)(x k, k ), not to the first-degree Lagrangian (1)(x k, k ). It is also shown that relations (2)(x k, k ) = ; and (2)(x k,p k) = are valid at any point of the ray and that they represent the group velocity surface and the slowness surface, respectively. All procedures and derived equations are valid for general anisotropic inhomogeneous media, and for general curvilinear coordinates x i. To make certain procedures and equations more transparent and objective, the simpler cases of isotropic and ellipsoidally anisotropic media are briefly discussed as special cases.  相似文献   

20.
This study develops a straightforward approximate method to estimate inelastic displacement ratio, C1 for base‐isolated structures subjected to near‐fault and far‐fault ground motions. Taking into account the inelastic behavior of isolator and superstructure, a 2 degrees of freedom model is employed. A total of 90 earthquake ground motions are selected and classified into different clusters according to the frequency content features of records represented by the peak ground acceleration to peak ground velocity ratio, Ap/Vp. A parametric study is conducted, and effective factors in C1 (i.e., fundamental vibration period of the superstructure, Ts; postyield stiffness ratio of the superstructure, αs; strength reduction ratio, R; vibration period of the isolator, Tb; strength of the isolator, Q; ratio of superstructure mass to total mass of the system, γm) are recognized. The results indicate that the practical range of C1 values could be expected for base‐isolated structures. Subsequently, effective parameters are included in simple predictive equations. Finally, the accuracy of the proposed approximate equations is evaluated and verified through error measurement, and comparisons are made in the analyses.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号