首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
Modeling of wave propagation induced by underground explosion   总被引:5,自引:0,他引:5  
A piecewise linear Drucker–Prager strength criterion and an isotropic continuum damage model with the damage scalar depending on an equivalent tensile strain are suggested to model rock mass behavior under blast loading. A rate-dependent constitutive relation is employed to model the energy dissipation caused by two sources, namely irreversible degradation of damage and permanent deformation caused by plasticity. The suggested model is incorporated with a commercially available software AUTODYN through its user’s subroutine function. Coupling of Euler and Lagrange processors are used to include all the materials under consideration such as explosive, air and rock mass, in the calculation. Using AUTODYN and the suggested model, shock wave propagation in rock mass induced by an underground explosion is simulated. Numerical results obtained agree favorably well with those obtained from an independently conducted field test. It demonstrates that the suggested model can be used to predict the damage area, plastic zone and ground motions generated by underground explosions.  相似文献   

2.
The year 2015 marks the bicentenary of the largest eruption in recent historic times: the 10–11 April 1815 eruption of Mount Tambora, Indonesia. Two hundred years after the eruption, an incomplete or inaccurate record of large eruptions over the past millennia, and uncertainties in determining the true sizes of eruptions, hamper our ability to predict when the next eruption of this scale may occur. Such events would have catastrophic effects locally and, possibly, world‐wide. The problem is compounded by a lack of detailed knowledge of how and over what timescales large magma reservoirs that feed such eruptions grow and assemble, and of the surface manifestations of these processes recorded through geophysical or geochemical monitoring techniques.  相似文献   

3.
Engineered barriers are basic elements in the design of repositories for the isolation of high‐level radioactive waste. This paper presents the thermo‐hydro‐mechanical (THM) analysis of a clay barrier subjected to heating and hydration. The study focuses on an ongoing large‐scale heating test, at almost full scale, which is being carried out at the CIEMAT laboratory under well‐controlled boundary conditions. The test is intensely instrumented and it has provided the opportunity to study in detail the evolution of the main THM variables over a long period of time. Comprehensive laboratory tests carried out in the context of the FEBEX and NF‐PRO projects have allowed the identification of the model parameters to describe the THM behaviour of the compacted expansive clay. A conventional THM approach that assumes the swelling clay as a single porosity medium has been initially adopted to analyse the evolution of the test. The model was able to predict correctly the global THM behaviour of the clay barrier in the short term (i.e. for times shorter than three years), but some model limitations were detected concerning the prediction of the long‐term hydration rate. An additional analysis of the test has been carried out using a double structure model to describe the actual behaviour of expansive clays. The double structure model explicitly considers the two dominant pore levels that actually exist in the FEBEX bentonite and it is able to account for the evolution of the material fabric. The simulation of the experiment using this enhanced model provides a more satisfactory reproduction of the long‐term experimental results. It also contributes to a better understanding of the observed test behaviour and it provides a physically based explanation for the very slow hydration of the barrier. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
Large‐scale soft‐sediment deformation structures occur within fluvial sandstone bodies of the Upper Cretaceous Wahweap Formation in the Kaiparowits basin, southern Utah, USA. These structures represent an exceptional example of metre‐scale fault‐proximal, seismogenic load structures in nearly homogenous sandstones. The load structures consist of two types: large‐scale load casts and wedge‐shaped load structures. Large‐scale load casts penetrate up to 4·5 m into the underlying sandstone bed. Wedge‐shaped load structures include metre‐scale, parallel, sub‐vertical features and decimetre‐scale features along the periphery of the large‐scale load casts or other wedge‐shaped load structures. Wedge‐shaped load structures contain well‐developed, medial cataclastic shear deformation bands. All load structures contain pervasive well‐defined millimetre‐thick to centimetre‐thick internal laminae, oriented parallel to the outside form of the load structures and asymptotic to deformation bands. Both types of load structures formed because of an inverted density profile, earthquake‐triggered liquefaction and growth of irregularities (a Rayleigh–Taylor instability) on the sandstone–sandstone erosional contact. The internal laminae and deformation bands formed during deformation and clearly demonstrate polyphase deformation, recording a transition from liquefied to hydroplastic to brittle modes of deformation. Decimetre‐scale wedge‐shaped load structures on the edge of the large‐scale load casts probably formed towards the end of a seismic event after the sediment dewatered and increased the frictional contact of grains enough to impart strength to the sands. Metre‐scale wedge‐shaped load structures were created as the tips of downward foundering sediments were driven into fractures, which widened incrementally with seismic pulsation. With each widening of the fracture, gravity and a suction effect would draw additional sediment into the fracture. Superimposed laminae indicate a secondary syndeformational origin for internal laminae, probably by flow‐generated shearing and vibrofluidization mechanisms. Large‐scale and wedge‐shaped load structures, polyphase deformation and secondary laminae may characterize soft‐sediment deformation in certain fault‐proximal settings.  相似文献   

5.
In this paper, a novel combination of well‐established numerical procedures is explored in order to accelerate the simulation of sequential excavation. Usually, large‐scale models are used to represent these problems. Due to the high number of equations involved, the solver algorithm represents the critical aspect which makes the simulation very time consuming. The mutable nature of the excavation models makes this problem even more pronounced. To accomplish the representation of geometrical and mechanical aspects in an efficient and simple manner, the proposed solution employs the boundary element method with a multiple‐region strategy. Together with this representational system, a segmented storage scheme and a time‐ordered tracking of the changes form an adequate basis for the usage of fast updating methods instead of frontal solvers. The present development employs the Sherman–Morrison–Woodbury method to speed up the calculation due to sequential changes. The efficiency of the proposed framework is illustrated through the simulation of test examples of 2D and 3D models. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

6.
The problem of finite element simulation of incompressible fluid flow in porous medium is considered. The porous medium is characterized by the X‐ray microtomography technique in three dimensions. The finite calculus‐based stabilization technique is reviewed to implement the equal order finite element interpolation functions for both velocity and pressure. A noble preconditioner, the nodal block diagonal preconditioner, is considered whose performance is thoroughly investigated. Combining this preconditioner with a standard iterative solver during the computational homogenization procedure, it is possible to carry out the large‐scale fluid flow simulation for estimating permeability of the porous medium with reasonable accuracy and reliability. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

7.
Soil–structure interaction problems are commonly encountered in geotechnical practice and remarkably characterized with significant material stiffness contrast. When solving the soil–structure interaction problems, the employed Krylov subspace iterative method may converge slowly or even fail, indicating that the adopted preconditioning method may not suit for such problems. The inexact block diagonal preconditioners proposed recently have been shown effective for the soil–structure interaction problems; however, they haven't been exploited to full capabilities. By using the same partition strategy according to the structure elements and soil elements, the partitioned block symmetric successive over‐relaxation preconditioners or partitioned block constraint preconditioners are proposed. Based on two pile‐group foundation problems and a tunnel problem, the proposed preconditioners are evaluated and compared with the available preconditioners for the consolidation analysis and the drained analysis, respectively. In spite of one additional solve associated with the structure block and multiplications with off‐diagonal blocks in the preconditioning step, numerical results reveal that the proposed preconditioners obviously possess better performance than the recently developed inexact block preconditioners. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

8.
Time‐domain analysis of dynamic soil–structure interaction based on the substructure method plays an increasing role in practical applications as compared with the frequency‐domain analysis. Efficient and accurate modelling of the unbounded soil or rock medium has been a key issue in such an analysis. This paper presents a subregional stepwise damping‐solvent extraction formulation for solving large‐scale dynamic soil–structure problems in the time domain. Accuracy and efficiency of the formulation are evaluated in detail for a classical problem involving a rigid strip foundation embedded in a half‐space. A practical large‐scale soil–structure interaction problem, which represents a high concrete gravity dam subjected to seismic load, is then analysed using the proposed method. Various responses of the dam, including time histories of the crest displacement and acceleration and contours of the peak principal stresses within the dam body, are presented. Comparisons are also made between these results with those obtained using other models for the unbounded medium. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

9.
10.
地应力场下岩体爆体的数值模拟   总被引:2,自引:0,他引:2  
刘艳  许金余 《岩土力学》2007,28(11):2485-2488
建立了地应力作用下岩体在爆炸荷载作用下的运动微分方程,并对岩体中的地应力采用隐式分析方法,对爆炸的高速动载作用采用显式分析方法,通过隐式-显式连续求解完成地应力作用下岩体爆炸的数值模拟。数值算例结果表明,近爆腔处,较低的静应力场与较高的动应力场作用时爆炸动态应力波占主导地位,初始应力对爆破应力波影响较小。在远离爆腔处,高地应力存在改变了爆破应力波的传播规律,对爆炸岩体有很大影响。  相似文献   

11.
The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

12.
This study presents a synthesis of the geomorphology, facies variability and depositional architecture of ice‐marginal deltas affected by rapid lake‐level change. The integration of digital elevation models, outcrop, borehole, ground‐penetrating radar and high‐resolution shear‐wave seismic data allows for a comprehensive analysis of these delta systems and provides information about the distinct types of deltaic facies and geometries generated under different lake‐level trends. The exposed delta sediments record mainly the phase of maximum lake level and subsequent lake drainage. The stair‐stepped profiles of the delta systems reflect the progressive basinward lobe deposition during forced regression when the lakes successively drained. Depending on the rate and magnitude of lake‐level fall, fan‐shaped, lobate or more digitate tongue‐like delta morphologies developed. Deposits of the stair‐stepped transgressive delta bodies are buried, downlapped and onlapped by the younger forced regressive deposits. The delta styles comprise both Gilbert‐type deltas and shoal‐water deltas. The sedimentary facies of the steep Gilbert‐type delta foresets include a wide range of gravity‐flow deposits. Delta deposits of the forced‐regressive phase are commonly dominated by coarse‐grained debrisflow deposits, indicating strong upslope erosion and cannibalization of older delta deposits. Deposits of supercritical turbidity currents are particularly common in sand‐rich Gilbert‐type deltas that formed during slow rises in lake level and during highstands. Foreset beds consist typically of laterally and vertically stacked deposits of antidunes and cyclic steps. The trigger mechanisms for these supercritical turbidity currents were both hyperpycnal meltwater flows and slope‐failure events. Shoal‐water deltas formed at low water depths during both low rates of lake‐level rise and forced regression. Deposition occurred from tractional flows. Transgressive mouthbars form laterally extensive sand‐rich delta bodies with a digitate, multi‐tongue morphology. In contrast, forced regressive gravelly shoal‐water deltas show a high dispersion of flow directions and form laterally overlapping delta lobes. Deformation structures in the forced‐regressive ice‐marginal deltas are mainly extensional features, including normal faults, small graben or half‐graben structures and shear‐deformation bands, which are related to gravitational delta tectonics, postglacial faulting during glacial‐isostatic adjustment, and crestal collapse above salt domes. A neotectonic component cannot be ruled out in some cases.  相似文献   

13.
A standing wave in front of a seawall may reach a height more than twice of its incident component. When excess pore pressure occurs, it may even induce seabed instability, hence endangering the structure. This issue was studied previously using only linear wave theory. In this paper, standing‐wave theory to a second‐order approximation is applied, in order to demonstrate the differences between these two solutions. The spatial and temporal variations in the instantaneous pore pressure are first calculated, in addition to their vertical distributions. The effects of wave height, water depth and the degree of soil saturation on pore pressure distributions are then discussed, followed by the net pore pressure averaged over one wave cycle. The results suggest the existence of a residual pore pressure in the seabed and its net pore pressure can be used to estimate the wave‐induced liquefaction potential in a soil column. It also indicates that, in deep water, the second‐order solution predicts that a negative pore pressure at an antinode which may be greater than a positive pressure. Overall, the second‐order solution is found to agree better with the experimental results of the pore pressures available, compared to the linear solution. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

14.
Evaluating the induced subsidence is a critical step in multi‐seam longwall mining. Numerical modelling can be a cost‐effective approach to this problem. Numerical evaluation of longwall mining‐induced subsidence is much more complicated when more than one seam is to be extracted. Only a few research works have dealt with this problem. This paper discusses the essential requirements of a robust numerical modelling approach to simulation of multi‐seam longwall mining‐induced subsidence. In light of these requirements, the previous works on this topic are critically reviewed. A simple yet robust FEM‐based modelling approach is also proposed that is capable of simulating caving process, rock mass deterioration and subsidence around multi‐seam excavations. The effectiveness of this approach in comparison with two other conventional FEM approaches is demonstrated through numerical examples of two different multi‐seam mining configurations. Results show that the proposed numerical modelling approach is the only robust method, which is capable of simulating multi‐seam subsidence in both demonstrated cases. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

15.
条形药包在隧道爆破中产生的应力场的实测分析   总被引:3,自引:0,他引:3  
傅洪贤 《岩土力学》2009,30(2):483-486
基于条形药包的爆破特性,针对隧道爆破的装药形式,把隧道爆破简化成几个条形药包的微差爆破,在隧道爆破区域的径向应力波为柱形,隧道掘进的轴向应力波为几个并列半圆形。应力波强度分布很不均匀,爆破区域的径向应力强度最大,向爆区端部延伸径向应力强度逐渐减小。进入爆区端部后,应力强度减小更明显。随着与爆区轴线夹角的减小,端部应力强度再次减小,爆区端部轴线方向应力强度最小。  相似文献   

16.
Extreme waves caused by tsunamis and storm surges can lead to soil failures in the near‐shore region, which may have severe impact on coastal environments and communities. Multiphase flows in deformable porous media involve several coupled processes and multiple time scales, which are challenging for numerical simulations. The objective of this study is to investigate the roles of the various processes and their interactions in multiphase flows in unsaturated soils under external wave loading, via theoretical time‐scale analysis and numerical simulations. A coupled geomechanics–multiphase flow model based on conservation laws is used. Theoretical analysis based on coupled and decoupled models demonstrates that transient and steady‐state responses are governed by pore pressure diffusion and saturation front propagation, respectively, and that the two processes are essentially decoupled. Numerical simulations suggest that the compressibility of the pore fluids and the deformation of the soil skeleton are important when the transient responses of the media are of concern, while the steady‐state responses are not sensitive to these factors. The responses obtained from the fully coupled numerical simulations are explained by a simplified time‐scale analysis based on coupled and decoupled models. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

17.
The beach‐bar reservoir play has become an important exploration target within the Bohai Bay Basin, especially in the Boxing Sag within the Dongying Depression, where a large‐scale lacustrine beach‐bar oil pool has been discovered recently. The sedimentary characteristics, distribution patterns and formation mechanisms of beach‐bar sand bodies in the upper fourth member of the Eocene Shahejie Formation (Es4s) in the Boxing Sag were studied in detail based on seismic, well log data and core data. The Es4s in the Boxing Sag is composed of a third‐order sequence consisting of three systems tracts, i.e. a lowstand systems tract, a transgressive systems tract and a highstand systems tract. Beach‐bar sand bodies were deposited widely in the basin during the lowstand systems tract period. The sandy beach‐bars are characterized by siltstones, fine‐grained silty sandstones interbedded with thin mudstone units. The presence of well‐developed sedimentary structures, such as swash bedding, parting lineation, parallel bedding, ellipsoidal mud clasts, ripples, terrestrial plant debris and vertical burrows, suggests that beach‐bars were deposited in a relatively shallow water environment under the influence of strong hydrodynamics. Laterally, the sandy beach facies occurred as a more continuous sheet‐like body around the sandy bar in most parts of the sag. Stratigraphically, beach‐bars were distributed mainly in the lowstand systems tract and they were less well‐developed in the transgressive systems tract and highstand systems tract. Several factors were probably responsible for the occurrence of the large‐scale beach‐bars during the lowstand systems tract period, including: (i) a gentle palaeoslope and relatively weak structural activities; (ii) a shallow‐water condition with a strong hydrodynamic environment; (iii) high‐frequency oscillations of the lake level; and (iv) an abundant terrigenous clastic feeding system with multiple‐point and linear sediment sources.  相似文献   

18.
Natural evaporite dissolution in the subsurface can lead to cavities having critical dimensions in the sense of mechanical stability. Geomechanical effects may be significant for people and infrastructures because the underground dissolution may lead to subsidence or collapse (sinkholes). The knowledge of the cavity evolution in space and time is thus crucial in many cases. In this paper, we describe the use of a local nonequilibrium diffuse interface model for solving dissolution problems involving multimoving interfaces within three phases, that is, solid–liquid–gas as found in superficial aquifers and karsts. This paper generalizes developments achieved in the fluid–solid case, that is, the saturated case [1]. On one hand, a local nonequilibrium dissolution porous medium theory allows to describe the solid–liquid interface as a diffuse layer characterized by the evolution of a phase indicator (e.g., porosity). On the other hand, the liquid–gas interface evolution is computed using a classical porous medium two‐phase flow model involving a phase saturation, that is, generalized Darcy's laws. Such a diffuse interface model formulation is suitable for the implementation of a finite element or finite volume numerical model on a fixed grid without an explicit treatment of the interface movement. A numerical model has been implemented using a finite volume formulation with adaptive meshing (e.g., adaptive mesh refinement), which improves significantly the computational efficiency and accuracy because fine gridding may be attached to the dissolution front. Finally, some examples of three‐phase dissolution problems including density effects are also provided to illustrate the interest of the proposed theoretical and numerical framework. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

19.
3D finite strain analyses and kinematic vorticity measurements were carried out on the Loghon Anticline within the HP‐LT Sanandaj–Sirjan metamorphic belt (Neyriz area, SW Iran). Rƒ/φ and Fry methods were used on the strain markers (e.g. deformed fossils) to interpret geometric relationships between the fold axis, strain ellipsoid axes and shear zone boundaries. The results indicate the predominance of prolate strain in the anticline. Quantitative kinematic analyses show that the Wk parameter is 0. 67 ± 0. 06 (i.e. pure‐shear dominated non‐coaxial flow). This study quantitatively supports the establishment of a dextral transpressive system, which is responsible for the development of the large‐scale right‐lateral shear zones that strike sub‐parallel to the major folds. Flexural shear combined with regional dextral‐shear is suggested to be the most common mechanism of folding in this area. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

20.
地下洞室开挖爆破围岩松动圈的数值分析计算   总被引:2,自引:2,他引:2  
肖明  张雨霆  陈俊涛  田华 《岩土力学》2010,31(8):2613-2618
根据数值计算提出了在地下洞室开挖爆破中,确定围岩松动圈的方法。给出了地下洞室开挖爆破的三维弹塑性损伤有限元计算方法。该方法论述了爆破荷载、重力荷载和开挖荷载的计算、施加以及迭代的方法,可以合理反应爆破对围岩稳定的影响和锚固支护施加的效果;结合岩石松动圈的实测方法,推导了判定围岩开始松动的损伤系数阈值公式,可以结合数值计算得到的洞周围岩损伤系数分布,确定在开挖爆破作用下围岩松动圈的范围。对实际工程的分析计算表明,根据数值计算判定的围岩松动范围与实测松动范围的规律基本一致,量值吻合较好,其结果为地下洞室开挖爆破的围岩松动圈的确定提供了有效的方法。  相似文献   

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

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