首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The stress and strain-rate fields in a glacier or other type of rock with known rheological properties can be calculated by using a numerical solution technique. Calculations are based on force-balance equations and the constitutive relation for polar ice, and do not involve any mathematical approximations. The geometry is prescribed. Basal velocities also are prescribed, either by specifying their values or by relying on a sliding relation. All other quantities are determined numerically.  相似文献   

2.
冰川流速是冰川动力状况的重要标志,利用合成孔径雷达技术能快速获得大范围冰川的表面流速.利用日本高级陆地观测卫星(ALOS)相控阵型合成孔径雷达L波段(PALSAR)及欧洲太空局的ENVISAT/ASAR数据的特征匹配方法获得帕米尔高原公格尔山区冰川表面流速,并结合合成孔径雷达干涉相干与不同时期数字高程模型对公格尔山区典型冰川动力进行分析,获得研究区不同时间基线冰川表面相干性、表面流速以及基于不同时相DEM的典型冰川表面高程变化信息.结果表明:30 a来克拉牙依拉克冰川表面高程下降了(15±12.1)m,表碛区域近期运动速度变化不大;其木干冰川平均表面高程几乎无变化,但2007-2011年夏季表面流速明显减缓,靠近末端附近部分区域可能已经演化为非活动区;姜满加尔冰川位于西风带的迎风坡,积累区面积大,冰川流速较快,无表碛覆盖,但表面高程仍下降了(8.8±12.7)m.编号为5Y663D0009的冰川冰舌表碛覆盖区可能已经演化为非活动区,30 a来表面高程下降(8.6±12.0)m.综合分析表明,冰川规模特别是积累区面积的大小及所处位置、地形对冰川演化具有重要影响.  相似文献   

3.
The action of tunnel excavation reduces the in-situ stresses along the excavated circumference and can therefore be simulated by unloading of cavities from the in-situ stress state. Increasing evidence suggests that soil behavior in the plane perpendicular to the tunnel axis can be modelled reasonably by a contracting cylindrical cavity, while movements ahead of an advancing tunnel heading can be better predicted by spherical cavity contraction theory. In the past, solutions for unloading of cavities from in-situ stresses in cohesive-frictional soils have mainly concentrated on the small strain, cylindrical cavity model. Large strain spherical cavity contraction solutions with a non-associated Mohr–Coulomb model do not seem to be widely available for tunnel applications. Also, cavity unloading solutions in undrained clays have been developed only in terms of total stresses with a linear elastic-perfectly plastic soil model. The total stress analyses do not account for the effects of strain hardening/softening, variable soil stiffness, and soil stress history (OCR). The effect of these simplifying assumptions on the predicted soil behavior around tunnels is not known. In this paper, analytical and semi-analytical solutions are presented for unloading of both cylindrical and spherical cavities from in-situ state of stresses under both drained and undrained conditions. The non-associated Mohr-Coulomb model and various critical state theories are used respectively to describe the drained and undrained stress-strain behaviors of the soils. The analytical solutions presented in this paper are developed in terms of large strain formulations. These solutions can be used to serve two main purposes: (1) to provide models for predicting soil behavior around tunnels; (2) to provide valuable benchmark solutions for verifying various numerical methods involving both Mohr–Coulomb and critical state plasticity models. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

4.
Little is known about the detailed behavior of glaciers in the Karakoram Mountains. Advanced land observing satellite (ALOS) phased array type L-band synthetic aperture radar (PALSAR) data were used to obtain the surface velocity of the Yengisogat Glacier in the Karakoram Mountains. Four ALOS PALSAR data sets with 46?days temporal baseline acquired from 2007 to 2009 covered all four seasons and were used to extract the offset fields and estimate annual average surface velocity based on seasonal velocities. For the ALOS PALSAR data the synthetic aperture radar (SAR) feature-tracking procedures within the GAMMA software were utilized instead of SAR interferometry because of low coherence in case of fast-moving glaciers or large time intervals between the image acquisitions. The accuracy of the measurements is discussed, and the measurements were consistent with previous results from optical imagery feature tracking. It was revealed that the south tributaries contributed to the main flow of the glacier, with the glacier surface velocities of the south tributaries moving more rapidly than the north tributaries. This was mainly attributed to the effect of the glacier??s aspect in the glacier long-term development point of view. Seasonal and spatial variations of the glacier surface velocity imply that the tributary South Skamri Glacier is probably surging. This has previously been mentioned by some researchers such as Copland et al. The Equilibrium Line Altitude was found to be at about 5,000?m a.s.l for south tributaries, estimated from the surface velocity distribution along the glacier centerline.  相似文献   

5.
We propose an extension of the shifted Grünwald-Letnikov method to solve fractional partial differential equations in the Caputo sense with arbitrary fractional order derivative α and with an advective term. The method uses the relation between Caputo and Riemann-Liouville definitions, the shifted Grünwald-Letnikov, and the traditional backward and forward finite difference method. The stability of the method is investigated for the implicit and explicit scheme with homogeneous boundary conditions, and a stability criterion is found for the advective-dispersive equation. An application of the method is used to solve contaminant diffusion and advective-dispersive problems. The numerical solution for the fractional diffusion and fractional advection-dispersion is compared with their respective analytical solutions for different time and space grid refinements. The diffusion simulation exhibited a good fit between the analytical and numerical solutions, with the explicit scheme going from stable to unstable as the time and space refinement changes. The fractional advection-dispersion application produced small deviations from the analytical solution. These deviations, however, are analogous to the numerical dispersions encountered in conventional finite difference solutions of the advection-dispersion equation. The new method is also compared with the traditional L2 method. Notably, an example that involves asymmetrical fractional conditions, a fractional diffusivity that depends on time, and a source term show how the methods compare. Overall, this study assesses the quality and easiness of use of the numerical method.  相似文献   

6.
An efficient numerical technique has been used to compute the deformation of pores of arbitrary shape embedded in a homogeneous elastic solid under the influence of applied stresses. The scheme is based on the boundary-element method, where single linear elements are used to generate solutions that satisfy prescribed boundary conditions. These solutions can be employed to describe the behavior of elastic moduli and other petrophysical properties in porous rocks. The numerical algorithm allows computation of the stress field induced by the pores in the solid. In this way, the effect of the interactions between pores caused by stress concentrations can be studied from a quantitative point of view. To test the algorithm, some interesting results are compared with existing models, for special cases available in the literature. Also, a model for the compressibility and porosity of sedimentary rocks, as a function of applied hydrostatic stress, was generated by mixing some realistic pore geometries generated with the numerical algorithm. Results were in good agreement with data obtained from selected samples of sandstones.  相似文献   

7.
Arching is a well known phenomenon, which effects stress developments which were investigated and compared using analytical and numerical solutions. Marston’s (1930) solution was extended to a generalised 3-dimensional rectangular stope and later modified for square and circular stopes for comparison with FLAC results. Aubertin et al. (2003) & Li et al. (2003) models were improved significantly by placing the backfill within narrow stopes as lifts or layers in numerical modelling where the normal stress variation with depth were found to be more realistic. The FLAC results were compared with analytical solutions which were developed by previous researchers and modified by the authors to evaluate the arching effects in backfilled placed in narrow and circular stopes. It appeared from the investigation herein that δ = 0.67 ϕ and K = K o condition gives a very close match with the numerical model solutions obtained from FLAC. Many laboratory tests were conducted to find out friction angles for four Australian mines, which were between 30 and 49 degrees.  相似文献   

8.
Numerical simulation of debris flow with application on hazard area mapping   总被引：3，自引：0，他引：3  
A numerical program developed for field application is presented in this paper. We use the generalized Julien and Lan [8] rheological model to simulate debris flows. Due to the derivative discontinuous nature of the constitutive law, flow is separated into plug region and bottom region (with stress greater than yield stress). The program solves the plug flow layer solution first, and then corrects the solution with the bottom layer approximation. Numerical scheme with upwind method and central difference in space and Adam–Bashforth third-order scheme in time is used for both layers. The scheme is tested against analytical solutions and laboratory experiments with very good results. Application to a field case with more complicated geometry also achieves good agreement, with error less than 5% compared to field measurements. The final example demonstrates how this numerical program is used in a preliminary design.  相似文献   

9.
Biot theory predicts wave velocities in a saturated granular medium using the pore geometry, viscosity, densities, and elastic moduli of the solid skeleton and pore fluid, neglecting the interaction between constituent particles and local flow, which becomes essential as the wavelength decreases. Here, a hydro-micromechanical model, for direct numerical simulations of wave propagation in saturated granular media, is implemented by two-way coupling the lattice Boltzmann method (LBM) and the discrete element method (DEM), which resolve the pore-scale hydrodynamics and intergranular behavior, respectively. The coupling scheme is benchmarked with the terminal velocity of a single sphere settling in a fluid. In order to mimic a small amplitude pressure wave entering a saturated granular medium, an oscillating pressure boundary on the fluid is implemented and benchmarked with the one-dimensional wave equation. The effects of input waveforms and frequencies on the dispersion relations in 3D saturated poroelastic media are investigated with granular face-centered-cubic crystals. Finally, the pressure and shear wave velocities predicted by the numerical model at various effective confining pressures are found to be in excellent agreement with Biot analytical solutions, including his prediction for slow compressional waves.  相似文献   

10.
A new analytical proof is presented for steady‐state seepage in recharged heterogeneous unconfined aquifers. The paper also presents a detailed procedure and important rules for performing correctly numerical studies of unsaturated seepage. Once a numerical solution is calibrated with field data, using a set of spatially distributed values for hydraulic conductivity K and effective infiltration EI, any new numerical analysis with a set of αK and αEI values, where α is a constant, yields an equally good calibration. However, if the effective porosities of each layer are unchanged, the groundwater velocities are multiplied by α, whereas the travel times are divided by α, which may help to select α in order to match known travel time data. This is a clear example of multiple solutions to an inverse problem. The paper underlines the role and the need to finely mesh unsaturated zones and also contacts between layers to reach the asymptotic convergence range, as it was carried out to verify the proof and as it should be completed to study any seepage problem. A few consequences of the new analytical proof and the rigorous procedure are shown with examples. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

11.
Rock glaciers are an end unit of the coarse debris system that links frost weathering and rockfall to debris transport in mountainous environments. By examining topographic and climatic controls on creep, a better understanding of rock glacier formation and evolution could be obtained. A series of statistical tests were run comparing contributing area parameters (width, length, area, slope, headwall height, etc.) and rock glacier variables (width, length, area, thickness, slope, creep and temperature) in the Front Range of Colorado. Results showed that rock glacier width had the strongest correlation with contributing area width due to an abundance of lobate rock glaciers (r = 0.78). Rock glacier area and contributing area size were also related (r = 0.74), but are also a function of talus production factors. Mean surface velocity appears to be better correlated with thickness (r = 0.58) and length (r = 0.62), rather than slope (r = 0.33). Mean annual air temperature had a stronger exponential relationship with maximum velocities (r = 0.76) compared to mean annual velocity (r = 0.60). Front Range rock glacier velocities tend to decrease with warming, suggesting that the volume of ice and debris, rates of shear in plastic layers, or melt water may have a greater influence on deformation. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

12.
The thermal impact of several kilometre-thick magmatic underplating in the lower continental crust is studied with analytical and numerical methods. Simple analytical solutions are derived for the thermal transient in the case of an infinite depth below the underplate and also for the case of a finite depth (down to the asthenosphere). It is shown that these solutions lead to simple approximations for when the transient surface heat flow is at its maximum, what the maximum is, and for how long the transient lasts. Even though these solutions assume that the underplate is emplaced instantaneously, they are useful in the interpretation of underplating over finite time spans. A numerical scheme is suggested for the modelling of underplating that handles both short time intervals as well as long intervals. The scheme treats magmatic underplating in a mass and energy conservative manner, and it is compared against the analytical solutions. Finally, the analytical and numerical results for thermal transients are applied to a transect from the Vøring margin (NE Atlantic), with respect to various degrees of early Cenozoic magmatic intrusion. It appears that more than half of the lower crustal body (LCB) in the Vøring margin must be magmatic underplating for the vitrinite reflectance to be substantially higher than for the non-magmatic case, where the LCB is assumed to comprise Caledonian crust.  相似文献   

13.
基于SAR的表碛覆盖型冰川边界定位研究   总被引：2，自引：0，他引：2  
使用光学图像进行表碛覆盖型冰川边界判断相对比较困难。采用日本高级陆地观测卫星（ALOS）携带的L波段相控阵型合成孔径雷达（PALSAR）数据的干涉相干对表碛覆盖型冰川边界进行判断,并使用ALOS PALSAR数据的特征匹配方法获得表面流速进行验证分析,发现公格尔山区5Y663D0009冰川表碛覆盖区呈现高相干性且运动速度十分缓慢,表明该表碛区域可能已经演化成非活动区;而该冰川中碛覆盖区则表现出低相干性,运动速度比较高（5 m/a）,表明相干性是有效的判断依据,利用PALSAR数据相干性及获得的表面流速可以区分表碛覆盖型冰川活动与非活动区域,使气候波动情景下该类型冰川的动态变化监测成为可能并对该方法的可靠性与不确定性进行了探讨。  相似文献   

14.
Landsat-7机载扫描行校正器（Scan Line Corrector,SLC）失效后的ETM+影像（SLC-OFF影像）约有22%的数据缺失,严重限制了该影像在冰川研究中的应用,特别对长期缺乏高质量遥感影像的高亚洲地区冰川运动连续监测产生较大影响。以Karakoram 中部最大的锡亚琴冰川为例,初步评估了Landsat-7 ETM+SLC-OFF影像在山地冰川表面流速提取方面的适用性和可行性。选取2009年和2010年两景SLC-OFF影像,运用局部直方图匹配法（Local Linear Histogram Match,LLHM）和加权线性回归法（Weighted Liner Regression,WLR）修复缺失数据条带,并利用亚像元互相关方法对修复后两景影像进行冰川表面流速估算。结果表明,LLHM和WLR两种方法均能有效修复冰川区Landsat-7 SLC-OFF影像,其冰流估算结果与同期Landsat-5TM 影像的冰流结果较为一致,三者冰流速估算精度分别为±5.9m ·a-1、±6.3m ·a-1和±4.0m ·a-1,验证了Landsat-7 SLC-OFF 影像在山地冰川流速监测中的应用潜力。  相似文献   

15.
We consider a system of nonlinear partial differential equations that arises in the modeling of two-phase flows in a porous medium. The phase velocities are modeled using a Brinkman regularization of the classical Darcy’s law. We propose a notion of weak solution for these equations and prove existence of these solutions. An efficient finite difference scheme is proposed and is shown to converge to the weak solutions of this system. The Darcy limit of the Brinkman regularization is studied numerically using the convergent finite difference scheme in two space dimensions as well as using both analytical and numerical tools in one space dimension. The results suggest that the Brinkman regularization may not approximate the accepted entropy solutions of the Darcy model and raise fundamental questions about the use of Brinkman type models in two-phase flows.  相似文献   

16.
Coupled finite element analyses of the consolidation and deformation around stone columns have been performed to assess the accuracy of different analytical solutions. The numerical model reproduces the hypotheses of the closed-form solutions. In the model, a rigid load is applied to a unit cell formed by a fully penetrating column and its surrounding soil, and simple elastic or elasto-plastic soil models are used. The surface settlement, the dissipation of the pore pressure and the vertical stress concentration on the column are studied. These soil responses are accurately estimated with closed-form solutions that properly include the radial and plastic strains in the column. However, the surrounding soil does not yield for usual conditions, which reasonably justifies the elastic soil behavior assumed in the analytical solutions. The differences between drained and consolidation analyses are also evaluated. Comparing the numerical results with the closed-form solutions illustrates the implications of the assumptions of each approach.  相似文献   

17.
Summary The measurement of surface profiles is presented as a useful and simple approach to classifying statistically the essential features of rock joints. After introducing the reader to some existing analytical joint contact models for normal loading, a discrete numerical technique is developed. Using this technique the mechanical behaviour of a number of different slate joints is examined. The functional relationships between nominal stress, stiffness, true contact area and initial aperture are shown, for this class of joints, to be surprisingly simple. Experimental evidence is used to substantiate the numerical results. From the point of view of in-situ joint stiffness and hydraulic conductivity, numerical predictions seem feasible provided the degree of mating at some known stress level can be determined.  相似文献   

18.
In this article a numerical solution for a three‐dimensional isotropic, viscoelastic half‐space subjected to concentrated surface stress loadings is synthesized with the aid of the Radon and Fourier integral transforms. Dynamic displacement and stress fields are computed for points at the surface and inside the domain. The analysis is performed in the frequency domain. Viscoelastic effects are incorporated by means of the elastic–viscoelastic correspondence principle. The equations of motion are solved in the Radon–Fourier transformed domain. Inverse transformations to the physical domain are accomplished numerically. The scheme used to perform the numerical inverse transformations is addressed. The solution is validated by comparison with results available in the literature. A set of original dynamic displacement and stress solutions for points within the half‐space is presented. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

19.
2009年9月对祁连山冷龙岭宁缠河3号冰川外围建立控制网, 于冰川表面布设了13根标志杆, 随后分别于2010年7月、2010年9月再次对设立在冰川表面的花杆点进行测量, 获取2009/2010年度、2009年9月-2010年7月与2010年7-9月3个时段宁缠河3号冰川表面流速.结果显示: 2009/2010年度宁缠河3号冰川最大流速出现在海拔4 430 m附近, 为3.76 m·a-1;2009年9月-2010年7月表面流速最大值出现在海拔4 430 m附近, 为0.32 m·月-1;2010年7-9月最大流速出现在海拔4 380 m附近, 为0.47 m·月-1.总体来看, 2009/2010年度宁缠河3号冰川纵剖面上流速变化较为缓和, 显示出流速随海拔变化而变化的规律. 但不同季节表面流速在纵剖面上的分布情况不同, 横剖面上主流线附近流速最大, 向冰川两边逐渐递减, 各观测点均平行于主流线方向向冰川末端运动, 表现出冰川运动一般规律. 在冰川表面运动速度观测区域内东南边缘流速略大于西北边缘, 同时与规模相近的冰川运动速度相比, 宁缠河3号冰川运动速度较大.  相似文献   

20.
Muztag Ata is a mountain known as the “Father of glaciers” in the west of China. However, its glacier flow pattern has never been studied. In this paper, the velocity of mountain glaciers on Muztag Ata is mapped by two methods (InSAR and offset-tacking) using SAR data. This map provides a detailed view of the features of glacier motion on Muztag Ata. A special effort was also made to assess the accuracy of the glacier velocity estimates. To ensure the accuracy of the derived results, the validation of the results was conducted over two glaciers (Saliymek Glacier and Kuksay Glacier). The two methods yielded similar results that agree well with each other. The patterns of glacier velocity on Muztag Ata suggest that mountain topography, westerly winds and the monsoon have a strong influence. Our results show that motion of glaciers in the southern and western regions of Muztag Ata is faster than that in the northern and eastern regions. Eleven glaciers were identified, along which glacier velocities exhibited distinctive behavior in terms of the spatial variability of the glacier motion. The occurrence of the local flow maxima and minima at consistent locations over different parts of different glaciers suggests that the subglacial topography, glacier size and glacier orientation affect the overall flow patterns. The velocities are very low in the surface debris cover of some glaciers, which suggests debris cover has an impact on glacier motion.  相似文献   

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

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