首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
提出了一种新的波函数组合法,将体波与面波特征函数组合求解不规则河谷地形在不同波型以不同角度入射的非一致位移输入问题.新方法严格满足不规则河谷散射表面自由条件,克服了目前波函数函数展开法不能完备表达二维不规则地形散射问题的缺点.数值算例表明波函数组合方法具有精确、稳定、计算量小的特点,可以应用于结构-地基动力相互作用的非一致输入问题.  相似文献   

2.
This study focuses on investigating spatial variation of ground motion that has great influence on the dynamic behavior of the large structures located on the surface topography. One of the most effective parameters on the spatial variation of ground motion is the difference between the arrival time of seismic waves in different points located on the abutments. In this research, a three-dimensional model of the Pacoima Dam site is prepared. The time domain 3 D boundary element method is used to apply non-uniform excitation at the dam supports. This model is subjected to vertically propagating incident SH and P waves. The time delay can be characterized by calculating the value of the time delay for which the cross-correlation between two records is maximized. Finally, to obtain the time delay in a topographic site, a function considering effective parameters such as the height from the canyon base, wave velocity and predominant frequency, is presented. Furthermore, a code was developed for generating the spatially variation of seismic ground motions. The results show that the proposed functions have an acceptable accuracy in estimating the time delay to generate non-uniform ground motion.  相似文献   

3.
Bed topography and grain size are predicted for steady, uniform flow in circular bends by consideration of the balance of fluid, gravity and frictional forces acting on bed load particles. Uniform flow pattern is adequately described by conventional hydraulic equations, with bed shear defined as that effectively acting on bed load grains. This analysis is used as a basis to predict bed topography and grain size for steady, non-uniform flow in non-circular bends (represented by a ‘sine-generated’ curve). The non-uniform flow pattern is calculated using the method of Engelund (1974a). Equilibrium bed form, hence sedimentary structure, is found by comparison of existing flow conditions with one of the schemes describing the hydraulic stability limits of the various bed forms. The model was compared with bankfull flow observations from a channel bend on the River South Esk, Scotland. Theoretical bed topography and velocity distribution were very close to the observed data. However, bed shear stress showed only a broad agreement, probably because of the use a constant friction coefficient value. Mean grain size distribution showed good agreement, but theory did not account adequately for gravel sizes in the talweg region and on the upstream, inner part of the bar, possibly due to theoretical underestimation of effective bed shear. Bed form and sedimentary structure are predicted well using the familiar stream power-grain size scheme. The behaviour of the model under unsteady uniform flow conditions in circular bends was analyzed, and suggests that any variation of grain size and bed topography with stage is likely to be limited to deeper parts of the channel.  相似文献   

4.
The last two decades have witnessed the development and application of well-balanced numerical models for shallow flows in natural rivers.However,until now there have been no such models for flows with non-uniform sediment transport.This paper presents a 1D well-balanced model to simulate flows and non-capacity transport of non-uniform sediment in alluvial rivers.The active layer formulation is adopted to resolve the change of bed sediment composition.In the framework of the finite volume Slope Llmiter Centred(SLIC) scheme,a surface gradient method is incorporated to attain well-balanced solutions to the governing equations.The proposed model is tested against typical cases with irregular topography,including the refilling of dredged trenches,aggradation due to sediment overloading and flood flow due to landslide dam failure.The agreement between the computed results and measured data is encouraging.Compared to a non-well-balanced model,the well-balanced model features improved performance in reproducing stage,velocity and bed deformation.It should find general applications for non-uniform sediment transport modelling in alluvial rivers,especially in mountain areas where the bed topography is mostly irregular.  相似文献   

5.
The Greek coastal town of Egion on 15 June 1995 was shaken by a strong, small epicentral distance, earthquake that caused heavy damages to buildings and loss of life. The damages were concentrated in the central elevated part of the town whereas the flat coastal region remained almost intact. This non-uniform distribution of damage is studied in this article in terms of surface topography effects by conducting seismic response analyses of a simplified 2-D profile of the town. A dynamic finite element code implementing the equivalent-linear soil behavior (FLUSHPLUS) was used for the analyses and it was found that the step-like topography amplified greatly the intensity of motion without affecting its frequency content. The analyses showed that the motion recorded by an accelerograph installed at the center of the town is in agreement with the computed values; they also indicated a particularly intense amplification close to the crest of the steep slope, where a multi-story RC residential building partially collapsed. In contrast, the level of motion was found to be low at the flat coastal zone of the town where the earthquake damages were insignificant. It is concluded that the characteristic surface topography of the town played an important role in modifying the intensity of base motion.  相似文献   

6.
Summary Theory for the transmission of electromagnetic waves in a non-uniform troposphere is presented. The atmospheric refractive index is allowed to vary both with height and range. The formulation is facilitated by the fact that the field within a piecewise (laterally) uniform section can be represented by a discrete sum of modes that are orthogonal in the transverse plane. Using the appropriate radial wave functions, mode conversion at a junction between two laterally uniform sections is considered. The analysis is then extended to multi-section structures. A few typical profiles of refractive index are examined numerically for both uniform and non-uniform models. One of the important findings is that the lateral non-uniformities in an elevated duct will enhance the overall transmission for ground-based transmitting and receiving antennas.  相似文献   

7.
The restriction of the Babuska–Brezzi stability criteria for the mixed formulation requires non-uniform interpolation of displacement and pore-pressure variables, which complicates the coding in the finite element analysis. This will certainly preclude the use of many otherwise useful elements. A possible solution to overcome this problem introduced earlier by the authors requires a semi-explicit form. However, the effectiveness of the semi-explicit procedure is dependent on the time step length used, which presents difficulties for transient analysis. In this paper, an unconditionally stable staggered implicit-implicit algorithm is developed, which is generalised from the semi-explicit form. For non-uniform mesh configurations, the results indicate a significant improvement compared with those obtained by using the semi-explicit procedure.  相似文献   

8.
An unconventional scheme is used to estimate the flexural rigidity, or equivalently the elastic thickness of the lithosphere, given the topography and gravity data. The flexural rigidity is the parameter that governs the flexural response of the lithosphere in the frame of the thin plate flexure model. The scheme is an alternative to the widely used calculation of admittance of topography (sea-floor or continental topography) and gravity, bearing some advantages which are explained in the paper. The scheme involves the inversion of the gravity data in order to formulate a model of the crust-mantle interface (CMI) undulations. In a second step the flexure parameter is then evaluated from the relation between topography and CMI variations. Instead of calculating the admittance function using a spectral analysis, a set of point-load response functions are used in order to retrieve the optimal flexure parameter. This has two main advantages: instabilities of the numerical admittance evaluation at wavenumbers with low spectral energy in the topography are overcome and the analysis can be made over an area which is not necessarily rectangular, as required for the spectral analysis. The proposed method allows a higher space resolution of elastic thickness than any spectral method. For validation, the numerical strategy is applied to the situation of a realistic synthetic model, where all inputs and outputs are known a priori. Finally the spatial variations of the elastic thickness are studied in an area across the Eastern Alps.  相似文献   

9.
Experiments on initiation of motion of different fractions of non-uniform sediment mixtures are reported. Four sediment mixes and reference transport method to define the threshold of bed particle movement are used in this analysis. Factors controlling the initiation of motion are analyzed. It is found that relative size of different grains with respect to a median size greatly affects the initiation of movement. The incipient motion of non-uniform sediment is very much dependent on the absolute size for coarser fractions. Considering the factors controlling the initiation of motion the experimental data sets along with the other available data yield a new relationship for the calculation of critical shear stress of non-uniform sediments. A set of equations have found that all the data points fall into a single line for the finer and for coarser fractions.  相似文献   

10.
Iceland is the type example of a ridge-centered hotspot. It is controversial whether the seismic anomaly beneath it originates in the lower mantle or the upper mantle. Some recent studies reported that the 660-km discontinuity beneath central Iceland is shallow relative to peripheral regions and this was interpreted as an effect of elevated temperature at that depth. We investigate topography of the major upper mantle discontinuities by separating the effects of the topography and volumetric velocity heterogeneity in P receiver functions from 55 seismograph stations. Our analysis demonstrates that a significant (at least 10-km) shallowing of the 660-km discontinuity is only possible in the case of improbably low seismic velocities in the mantle transition zone beneath central Iceland. If, as in previous studies, lateral velocity variations in the mantle transition zone are neglected, the data require a depressed rather than an uplifted 660-km discontinuity. For a reasonable S-wave velocity anomaly in the mantle transition zone (around − 3%) no topography on the 660-km discontinuity is required. This can be explained by the lack of temperature anomaly or an effect of two phase transitions with opposite Clapeyron slopes.  相似文献   

11.
Non-uniform ground motions are generated based on a single record available at a site and seismic wave scattering analysis. The Chino Hills 2008 earthquake records at the Pacoima Dam site are used to indicate the accuracy of the method. Dynamic analysis of the Pacoima dam-reservoir-foundation under uniform and non-uniform ground motions is carried out using the EACD-3D2008 software, and the results are compared to recorded responses at different locations on the dam. There is good agreement between computed and recorded displacements of the dam for non-uniform excitation. For uniform excitation, the displacements are underestimated in comparison with those obtained from recorded excitation. Significant intensification of stresses, especially near the foundation, and different patterns of stress distribution are observed for non-uniform excitation in comparison with uniform excitation. For uniform excitation maximum stresses occur in the crown cantilever near the crest, but for non-uniform excitation the maximum stresses occur along the sides and near the foundation.  相似文献   

12.
三维起伏地形条件下航空瞬变电磁响应特征研究   总被引:1,自引:0,他引:1       下载免费PDF全文
航空瞬变电磁法以其速度快、成本低、通行性好等的优势能够有效的应用于地质地形条件复杂的地区.目前对于航空瞬变电磁法的研究主要基于平坦地形的理想情况,对于地形效应的研究相对较少,然而实际应用中地形不可避免,若忽略地形影响将对资料解释造成较大的误差,从而制约航空电磁方法的进一步发展.本文基于交错网格的时域有限差分方法对三维起伏地形条件下航空瞬变电磁进行正演模拟,在保证算法准确性的前提下给出大量模型算例.以经典地形模型为例,利用所给方法计算三维正演响应,结果显示起伏对于航空瞬变电磁数据有着显著的影响且影响主要集中在早期.而后,以实际地质资料为基础,构建起伏地形条件下包含多个异常体的三维复杂模型,计算了复杂模型的航空瞬变电磁响应,并给出三维全域视电阻率曲线,从而对地形效应的影响有了更加直观的认知.最后,通过大量模型讨论了地形的尺寸参数、电性参数、飞行轨迹与飞行高度等因素变化对于航空瞬变电磁数据的影响情况,并得出有价值的结论.  相似文献   

13.
This paper presents a lattice Boltzmann model (LBM) for 2-D advection and anisotropic dispersion equation (AADE) based on the Bhatnagar, Gross and Krook (BGK) model. In the proposed model, the particle speed space is discretized using a rectangular lattice that has four speeds in nine directions, and the single relaxation time is assumed to be directionally dependent. To ensure that the collision is mass-invariant when the relaxation time is directionally dependent, the concentration is calculated from a weighted summation of the particle distribution functions. The proposed model was verified against benchmark problems and the finite difference solution of solute transport with spatially variable dispersion coefficients and non-uniform velocity field. The significant results are that it conserves mass perfectly and offers accurate and efficient solutions for both dispersion-dominated and advection-dominated problems.  相似文献   

14.
In this work, a hybrid boundary integral equation method (BIEM) is developed, based on both displacement and hypersingular traction formulations, for the analysis of time-harmonic seismic waves propagating through cracked, multi-layered geological regions with surface topography and under plane strain conditions. Specifically, the displacement-based BIEM is used for a multi-layered deposit with interface cracks, while the regularized, traction-based BIEM is used when internal cracks are present within the layers. The standard uni-dimensional boundary element with parabolic shape functions is employed for discretizing the free surface and the layer interfaces, while special discontinuous boundary elements are placed near the crack tips to model the asymptotic behaviour of both displacements and tractions. This formulation yields displacement amplitudes and phase angles on the free surface of a geological deposit, as well as stress intensity factors near the tips of the cracks. Finally, in the companion paper, numerical results are presented which show that both scattered wave and stress concentration fields are sensitive to the incidence seismic wave parameters and to specific site conditions such as surface topography, layering, the presence of cracks and crack interaction.  相似文献   

15.
Abstract

A class of long planetary waves in a zonal channel analogous to the solitary and cnoidal waves of surface and internal gravity wave theory is discussed. On a mid-latitude β-plane, such waves exist as the result of divergence, non-uniform zonal velocity fields or bottom topography. In all cases studied the wave profile along the channel was found to satisfy the Korteweg-de Vries equation.  相似文献   

16.
Due to the increasing popularity of analyzing empirical Green’s functions obtained from ambient seismic noise, more and more regional tomographical studies based on short-period surface waves are published. Results could potentially be biased in mountainous regions where topography is not small compared to the wavelength and penetration depth of the considered waves. We investigate the effect of topography on the propagation of short-period Rayleigh waves empirically by means of synthetic data using a spectral element code and a 3-D model with real topography. We show that topography along a profile through the studied area can result in an underestimation of phase velocities of up to about 0.7% at the shortest investigated period (3 s). Contrary to the expectation that this bias results from the increased surface distance along topography, we find that this error can be estimated by local topographic contrasts in the vicinity of the receiver alone. We discuss and generalize our results by considering topographic profiles through other mountain ranges and find that southern Norway is a good proxy to assess the topography effect. Nevertheless, topographic bias on phase velocity measurements is in general not large enough to significantly affect recovered velocity variations in the ambient noise frequency range.  相似文献   

17.
ABSTRACT

An infiltration model based on the Green-Ampt assumptions is developed for a class of non-uniform soils in which saturated hydraulic conductivity decreases as an exponential function of depth, and the storage-suction factor C = ΔθΔΨ is a constant. An analysis of measured hydraulic conductivity and porosity data suggests that the model may be useful in some soils where the changes of hydraulic characteristics with depth are of this form. A method for using the model with time variable rainfall rates is given. The model is applied to simulate the experimental results of Childs &; Bybordi (1969) who measured infiltration into layered sand profiles. Agreement was satisfactory even for this case.  相似文献   

18.
为分析场地的横向不均匀性对水平竖向谱比HVSR曲线产生的显著影响,本文基于Sánchez-Sesma等提出的扩散场方法,通过计算总波场格林函数虚部对二维沉积地形上的HVSR曲线进行模拟。格林函数虚部则通过刚度矩阵和平面内斜线格林函数采用间接边界元方法进行求解。对二维沉积地形和相应的一维层状半空间的HVSR曲线进行了参数分析,着重讨论了沉积地形的形状、相对计算点位置等因素对HVSR曲线的影响规律。结果表明:沉积地形内外材料阻抗比对HVSR曲线的影响最为显著;随着沉积地形内外材料阻抗差异和沉积侧界面坡度的增大,沉积地形上HVSR曲线的第一峰值点的频率显著增大至相应层状半空间结果的3.3倍,同时HVSR曲线的形态呈现出平台现象;随着计算点到沉积边界距离的减小,HVSR曲线高频段幅值相对较大。根据本文得到的局部地形对HVSR曲线的影响规律,在进行场地勘探时可采用HVSR方法初步确定局部地形的分布位置以降低勘探成本。   相似文献   

19.
A generalization of the Timoshenko beam model is presented which accounts for the influence of the shear strains, due to non-uniform bending and torsion, on the flexural-torsional vibrations of thin-walled cores with open or closed cross-section. The axial displacement field incorporates the torsion secondary warping as well as the warping terms depending on the shear resultants. It is shown that exact solutions for the interior domain problem can be obtained under proper load conditions. Moreover, a discrete model for the free-vibration analysis is derived by adopting a linear interpolation of the unknown functions and a reduced integration in order to avoid locking phenomena. Various applications are developed, including the case of the coupled vibrations of a shear-core-steel-frame building.  相似文献   

20.
Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.  相似文献   

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

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