Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method          下载免费PDF全文

Mohammad Hossein Bazyar  Abbas Talebi 《国际地质力学数值与分析法杂志》2015,39(1):1-22

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Liquefaction potential evaluations by energy-based method and stress-based method for various ground motions     

《Soil Dynamics and Earthquake Engineering》2015

An energy-based liquefaction potential evaluation method (EBM) previously developed was applied to a uniform sand model shaken by seismic motions recorded at different sites during different magnitude earthquakes. It was also applied to actual liquefaction case histories in Urayasu city during the 2011 M9.0 Tohoku earthquake and in Tanno-cho during the 2003 M8.0 Tokachi-oki earthquake. In all these evaluations, the results were compared with those by the currently used stress-based method (SBM) under exactly the same seismic and geotechnical conditions. It was found that EBM yields similar results with SBM for several ground motions of recent earthquakes but has easier applicability without considering associated parameters. In Urayasu city, the two methods yielded nearly consistent results by using an appropriate coefficient in SBM for the M9.0 earthquake, though both overestimated the actual liquefaction performance, probably because effects of plasticity and aging on in situ liquefaction strength were not taken into account. In Tanno-cho, EBM could evaluate actual liquefaction performance due to a small-acceleration motion during a far-field large magnitude earthquake while SBM could not.  相似文献   

用随机有限断层方法模拟2016-11-13新西兰Mw7.8地震     

李启成  何书耕  闵也  郑新娟 《大地测量与地球动力学》2019,39(12):1237-1242

在未来大地震发震构造已知的前提下，提出用经验公式确定大地震断层长度、宽度、地震矩和滑动分布等参数，并用随机有限断层方法预测未来大地震。用随机有限断层方法模拟2016-11-13新西兰M_W7.8地震的12个基岩台站地震动时程和反应谱，并用模拟的频谱幅值与记录频谱幅值平均比值确定模拟误差。结果表明，在周期为0~10 s的范围内，模拟误差在0.92~1.08之间，不同频率模拟误差的标准差变化范围不超过1。95%的置信区间宽度随频率无明显变化，模拟结果反映了地震动记录的平均效果。用准随机方法重新获取新西兰大地震断层初始破裂点和位错滑动分布，模拟地震动并计算模拟误差随频率的变化，得到与最初模型相似的结论，从而进一步证实了本文提出的获取震源参数并用随机有限断层方法预测未来大地震的可靠性，且其特别适用于大地震远场模拟。  相似文献   

Hydrodynamic Tidal Model of the Persian Gulf Based on Spatially Variable Bed Friction Coefficient     

Mostafa Sohrabi Athar  Alireza A. Ardalan  Roohollah Karimi 《Marine Geodesy》2019,42(1):25-45

The aim of this study is to develop a two-dimensional hydrodynamic tidal model for the Persian Gulf (PG2017) using 2D-MIKE21 software. The advantages of present study is accounting for the spatial variation of bed friction coefficient besides a precise bathymetry together with a 23-year of combined records of satellite altimetry data. We found that the bed friction coefficient has a significant effect on sea level changes in the region under our modeling consideration. Since the tidal behavior in the northern part of the Qeshm Island is significantly different from the other parts of the Persian Gulf, to present a more accurate hydrodynamic tidal model, the Gulf is divided into two regions where the bed friction coefficient is modeled separately for each region. The root mean square value of the differences between the amplitude of dominant constituents; M2, S2, K1, and O1 derived from the PG2017 model and that of 98 altimetry and coastal tide gauge stations are respectively equal to 1.6, 1.9, 2.8, and 1.3?cm. Moreover, comparing the PG2017 model efficiency with the FES2014, OSU12, EOT11a, DTU10, and Admiralty models shows that the PG2017 model has an improvement of 22.1%, 47.2%, 43.2%, 44.2%, and 57.6% in terms of relative error, respectively.  相似文献   

Designing geotechnical structures with a proper stability criterion as a safety factor     

《Computers and Geotechnics》2016

  相似文献   

Discrete element model for crack propagation in brittle materials          下载免费PDF全文

Ba Danh Le  Georg Koval  Cyrille Chazallon 《国际地质力学数值与分析法杂志》2016,40(4):583-595

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

青藏高原东南缘构造应力场重构     

侯强  丁小军  赵宏  欧明霖  张博康 《大地测量与地球动力学》2019,39(5):458-463

采用有限元方法，针对青藏高原东南缘建立更细致、更精确的三维有限元弹性模型。选取9种不同的应力边界条件，分别进行优化分析后处理，将对应台站形变模拟值与GPS实测值进行误差分析，最终选取最佳方案作为古构造应力场。结果表明，青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用，特别是印度板块NNE向强烈碰撞作用，成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源，控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括几个力源：西北部青藏高原侧向挤压造成的WE向应力约105 MPa；西南部直接来自于印度板块的NE-WS向应力约70 MPa；南部NS向作用力33 MPa；东南部扬子块体侧向NW-SE阻挡力56 MPa；北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘，形成现今复杂应力场。
  相似文献   

An efficient hybrid integral-equation method for point-absorber wave energy converters with a vertical axis of symmetry     

《Applied Ocean Research》2019

To assist in the prototyping and controller design of point-absorber wave energy converters (WECs), an easy-to-implement hybrid integral-equation method is presented for computing the frequency-domain hydrodynamic properties of bodies with a vertical axis of symmetry in waves. The current hybrid method decomposes the flow domain into two parts: an inner domain containing the body and an outer domain extending to infinity. The solution in the inner domain is computed using the boundary-element method, and the outer-domain solution is expressed using eigenfunctions. Proper matching at the domain boundary is achieved by enforcing continuity of velocity potential and its normal derivative. Body symmetry allows efficient computation using ring sources in the inner domain. The current method is successfully applied to three different body geometries including a vertical truncated floating cylinder, the McIver toroid, and the coaxial-cylinder WEC being developed in the authors’ laboratory. In particular, the current results indicate that, by replacing the flat bottom of the coaxial-cylinder WEC with the Berkeley-Wedge (BW) shape, viscous effect can be significantly reduced with only minor negative impact on wave-exciting force, thus increasing WEC efficiency. Finally, by comparing to experimental measurements, the current method is demonstrated to accurately predict the heave added mass and wave-exciting force on the coaxial-cylinder WEC with BW geometry. If a viscous damping correction factor is used, the heave motion amplitude can also be accurately computed.  相似文献   

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution          下载免费PDF全文

Wencheng Jin  Hao Xu  Chloé Arson  Seth Busetti 《国际地质力学数值与分析法杂志》2017,41(2):223-250

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Macropores as preferential flow paths in meander bends          下载免费PDF全文

Garrett T. Menichino  Adam S. Ward  Erich T. Hester 《水文研究》2014,28(3):482-495

Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献