Fracture mechanics analysis of multiple cracks in anisotropic media     

Chia‐Huei Tu  C.‐S. Chen  T.‐T. Yu 《国际地质力学数值与分析法杂志》2011,35(11):1226-1242

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading          下载免费PDF全文

K. Duan  C.Y. Kwok  M. Pierce 《国际地质力学数值与分析法杂志》2016,40(8):1150-1183

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures          下载免费PDF全文

Wenjun Gong  Yunsheng Wang  Youjun Ning  Yonghong Luo  Yujie Wang 《国际地质力学数值与分析法杂志》2017,41(10):1267-1282

The effects of fractures on wave propagation problems are increasingly abstracting the attention of scholars and engineers in rock engineering field. This study aims to fully validate the ability of discontinuous deformation analysis (DDA) to model normal P‐wave propagation across rock fractures. The effects of a single fracture and multiple parallel fractures are all tested. The results indicate that DDA can accurately reflect the fracture effects, including the fractures stiffness, the fracture spacing and the fracture number, and the effects of incident wave frequency on one‐dimensional P‐wave propagation problems. Thus, DDA is able to deal well with normal incident P‐wave propagation problems. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

Investigation into the effect of fault properties on wave transmission          下载免费PDF全文

H.R. Mohammadi  H. Mansouri  M. Bahaaddini  H. Jalalifar 《国际地质力学数值与分析法杂志》2017,41(17):1741-1757

The estimation of wave transmission across the fractured rock masses is of great importance for rock engineers to assess the stability of rock slopes in open pit mines. Presence of fault, as a major discontinuity, in the jointed rock mass can significantly impact on the peak particle velocity and transmission of blast waves, particularly where a fault contains a thick infilling with weak mechanical properties. This paper aims to study the effect of fault properties on transmission of blasting waves using the distinct element method. First, a validation study was carried out on the wave transmission across a single joint and different rock mediums through undertaking a comparative study against analytical models. Then, the transmission of blast wave across a fault with thick infilling in the Golgohar iron mine, Iran, was numerically studied, and the results were compared with the field measurements. The blast wave was numerically simulated using a hybrid finite element and finite difference code which then the outcome was used as the input for the distinct element method analysis. The measured uplift of hanging wall, as a result of wave transmission across the fault, in the numerical model agrees well with the recorded field measurement. Finally, the validated numerical model was used to study the effect of fault properties on wave transmission. It was found that the fault inclination angle is the most effective parameter on the peak particle velocity and uplift. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data          下载免费PDF全文

Antonio‐Juan Collados‐Lara  Eulogio Pardo‐Igúzquiza  David Pulido‐Velazquez 《水文研究》2017,31(10):1966-1982

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.  相似文献   

Using raw regional climate model outputs for quantifying climate change impacts on hydrology          下载免费PDF全文

Jie Chen  François P. Brissette  Pan Liu  Jun Xia 《水文研究》2017,31(24):4398-4413

General circulation model outputs are rarely used directly for quantifying climate change impacts on hydrology, due to their coarse resolution and inherent bias. Bias correction methods are usually applied to correct the statistical deviations of climate model outputs from the observed data. However, the use of bias correction methods for impact studies is often disputable, due to the lack of physical basis and the bias nonstationarity of climate model outputs. With the improvement in model resolution and reliability, it is now possible to investigate the direct use of regional climate model (RCM) outputs for impact studies. This study proposes an approach to use RCM simulations directly for quantifying the hydrological impacts of climate change over North America. With this method, a hydrological model (HSAMI) is specifically calibrated using the RCM simulations at the recent past period. The change in hydrological regimes for a future period (2041–2065) over the reference (1971–1995), simulated using bias‐corrected and nonbias‐corrected simulations, is compared using mean flow, spring high flow, and summer–autumn low flow as indicators. Three RCMs driven by three different general circulation models are used to investigate the uncertainty of hydrological simulations associated with the choice of a bias‐corrected or nonbias‐corrected RCM simulation. The results indicate that the uncertainty envelope is generally watershed and indicator dependent. It is difficult to draw a firm conclusion about whether one method is better than the other. In other words, the bias correction method could bring further uncertainty to future hydrological simulations, in addition to uncertainty related to the choice of a bias correction method. This implies that the nonbias‐corrected results should be provided to end users along with the bias‐corrected ones, along with a detailed explanation of the bias correction procedure. This information would be especially helpful to assist end users in making the most informed decisions.  相似文献   

碎冰条件下单点系泊船舶动态响应特性研究     

杨佳东  吕海宁  王晋 《海洋工程》2019,37(2):29-39

为了研究不同单点系泊位置对极地系泊船舶在来冰角度改变时动态响应的影响,运用离散元方法(DEM),计算了船舶在60%密集度下,冰速为1.0 m/s、冰厚为1 m时不同角度的船冰阻力,结合系泊分析软件二次开发功能,计算分析系泊船舶在不同来冰角度和不同系泊位置下受上述冰情影响下的响应。通过对模拟结果与试验结果的对比,认为离散元方法(DEM)的数值模拟结果与试验结果吻合较好;单点系泊位置为距船头L/4时,系泊船舶能以最快速、最稳定的方式对来冰方向的变化做出响应。  相似文献   

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%的置信区间宽度随频率无明显变化，模拟结果反映了地震动记录的平均效果。用准随机方法重新获取新西兰大地震断层初始破裂点和位错滑动分布，模拟地震动并计算模拟误差随频率的变化，得到与最初模型相似的结论，从而进一步证实了本文提出的获取震源参数并用随机有限断层方法预测未来大地震的可靠性，且其特别适用于大地震远场模拟。  相似文献