三峡库区龙门寨危岩体崩塌产生涌浪研究   总被引：1，自引：1，他引：0  

郑嘉豪  黄波林  张全  赵海林  冯万里  王健  陈小婷 《地质力学学报》2020,26(4):533-543

长江两岸高耸的危岩体对航道、沿岸居民带来巨大安全隐患。大宁河属于长江一级支流，龙门寨危岩体位于大宁河上，距离巫山县城仅1 km。利用FLOW-3D软件，模拟了145 m、175 m两种水位工况下龙门寨危岩体崩塌产生涌浪过程和涌浪传播过程。模拟结果表明，涌浪在145 m水位工况下最大浪高约为17.9 m，175 m水位工况下最大浪高约为11.6 m；在巫山县的五个码头处，两种水位工况最大涌浪爬高分别约为10.9 m、3.8 m；根据涌浪高度，对大宁河进行危险分区，145 m水位工况下极高危险区长度约4.4 km，很高危险区长度约1.9 km；175 m水位工况下极高危险区长度约3.0 km，很高危险区长度约1.0 km。研究结果有助于防控龙门寨危岩体潜在涌浪灾害危害，保障大宁河航道和巫山县码头安全，同时也为三峡库区滑坡涌浪灾害提供了预警依据。   相似文献   

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses          下载免费PDF全文

J.C. Li  N.N. Li  S.B. Chai  H.B. Li 《国际地质力学数值与分析法杂志》2018,42(1):95-109

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.  相似文献   

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials          下载免费PDF全文

Chao Jiang  Gao‐Feng Zhao 《国际地质力学数值与分析法杂志》2018,42(4):674-693

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.  相似文献   

南亚高压区Rossby波能量向上穿越对流层顶传播的特征与机制          下载免费PDF全文

俞子闲  郭栋  李丽平  施春华  金鑫  李震坤  陈丹  林祥 《地球物理学报》2020,63(1):73-85

夏季平流层盛行强东风,Rossby波能量难以从对流层向上传播至平流层,而冬季平流层盛行西风,Rossby波能量容易上传,因此以往对Rossby波能量向平流层传播的研究多考虑冬季的情况.而事实上,因为夏季高原上空南亚高压反气旋环流,并非只有强东风存在,所以Rossby波能量也可能在南亚高压区向上传播,从而影响平流层的温度、风场及大气成分等.因此,本文利用ERA-interim逐日再分析资料,分析了1979—2015年夏季南亚高压区Rossby波能量穿越对流层顶传播的特征与机制.结果表明:Rossby波能量可以从南亚高压西北部的窗口区上传至平流层,最高可到达平流层顶,而在南亚高压的其他部分,Rossby波能量均不能穿越对流层顶上传或穿越对流层顶后无法继续上传.南亚高压西北区Rossby波能量可以穿越对流层顶传播的原因是盛行西风,且西风急流出现的频率很小,同时涡动热量通量异常引起的垂直分量的第一项对其上传有很大贡献.南亚高压东北区也盛行西风,然而Rossby波能量不能向上穿越对流层顶的原因是强西风出现频率较高,且温度脊与高度脊位相相近,不利于上传.南亚高压南部均盛行东风,在平流层中下层均为稳定层结,因此Rossby波能量很难上传.南亚高压西南区在对流层位于青藏高原环流的伊朗高原下沉区附近,层结稳定,并且温度脊超前于高度脊,所以Rossby波能量很难上传.而南亚高压东南区在对流层位于南海-西太平洋热带幅合带,层结不稳定,存在Rossby波能量较弱的上传,达到对流层顶后无法继续上传,该区域温度脊落后于高度脊的温压场配置也为Rossby波能量在对流层内的传播提供了条件.  相似文献   

SH波地表散射的局域边界元模拟          下载免费PDF全文

盖增喜  陈晨旭 《地球物理学报》2020,63(6):2274-2280

本文提出了一种计算不规则起伏地形中SH波散射的有效方法——局域边界元法.本方法基于传统边界元法,为计算复杂地表散射问题提供了一种更加高效的解决方案.根据地震波满足的边界积分方程中牵引力格林函数的特性,我们将自由边界分解成水平部分和起伏部分.通过公式推导,可将水平部分的位移由起伏部分的位移通过格林函数线性叠加表示,因此只需对起伏部分的位移进行直接求解,从而极大地减少了待求解的未知数个数,显著提高了计算效率.通过与半圆形山谷SH波平面波入射的解析解比较,验证了方法的正确性.数值模型比较显示,局域边界元模拟结果与传统边界元数值解完全吻合,但是大幅提高了计算效率.因此,局域边界元法可以作为模拟不规则地形中地震波散射的有效工具.  相似文献   

Modelling of viscoacoustic wave propagation in transversely isotropic media using decoupled fractional Laplacians     

Zhihao Qiao  Chengyu Sun  Jie Tang 《Geophysical Prospecting》2020,68(8):2400-2418

A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media.  相似文献   

木本植被覆盖岸坡上波浪爬升过程的数值模拟研究   总被引：1，自引：1，他引：0  

徐海珏  胡萍  白玉川  杨波 《海洋学报》2020,42(3):10-24

近岸木本植物构成的生态缓冲带作为新型的海岸软防护结构，兼具功能性和生态友好性，在沿海工程建设中愈发受到关注，如何深入开展其防护效果的机理研究是目前亟待解决的问题。本文采用数值模拟方法，在N-S方程中分别考虑树枝和树干的拖曳力影响，提出了木本植被作用下波浪沿斜坡爬升的表面波衰减的连续介质等效模型，并采用MAC法来跟踪自由曲面上的水颗粒轨迹。本文以波浪沿1/30的斜坡爬升为算例，对比讨论了有无植被作用下波浪的传播过程，并将算例结果与以往试验结果规律进行对照，验证了数值模型的有效性。最后，分别讨论了植物枝干的高度、密度、树枝倾斜角度等植被特性和波浪因素对植被消浪效果的影响，得到植被消浪的基本规律。文中的计算结果也可为实际的护岸工程和生态景观设计提供参考。  相似文献   

基于GPU加速的Boussinesq类波浪传播变形数值模型     

孙家文  朱桐  房克照  刘忠波 《海洋工程》2020,38(2):111-119

Boussinesq波浪模型是一类相位解析模型,在时域内求解需要较高的空间和时间分辨率以保证计算精度。为提高计算效率,有必要针对该类模型开展并行算法的研究。与传统的中央处理器(CPU)相比,图形处理器(GPU)有大量的运算器,可显著提高计算效率。基于统一计算设备架构CUDA C语言和图形处理器,实现了Boussinesq模型的并行运算。将本模型的计算结果同CPU数值模拟结果和解析解相比较,发现得到的结果基本一致。同时也比较了CPU端与GPU端的计算效率,结果表明,GPU数值模型的计算效率有明显提升,并且伴随数值网格的增多,提升效果更为明显。  相似文献   

基于离散元方法对走滑拉分盆地演化及次级断裂扩展过程研究     

刘源  KONIETZKY Heinz 《地质力学学报》2019,25(5):840-852

拉分盆地是一种与走滑断裂带密切相关的特殊拉张构造，因其重要的构造意义，及其与火山活动、中小地震群集、特殊的成矿作用间的伴生关系而受到研究者的高度重视。关于拉分盆地的形成演化过程，已有较多的研究成果，但是由于研究手段的限制，缺少对盆地演化中次级断裂扩展过程的研究。基于离散元的数值计算方法是研究断裂扩展方式的理想方法。本文采用基于离散元的颗粒流方法，揭示纯走滑拉分盆地发育过程中的断裂扩展和连接过程，为拉分盆地演化机理和断裂扩展提供新的研究方法。同时，根据主走滑断层与块体运动方向的夹角不同，建立不同的张扭性拉分盆地模型，系统研究张扭性盆地的断裂扩展和演化机理。将上述理论研究结果与死海盆地等经典拉分盆地实例相结合，探讨了死海盆地、土耳其Cinarcik盆地、哥伦比亚El Paraiso盆地等的形成演化机理和断裂扩展方式。   相似文献