首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 41 毫秒
1.
Large deformations and discontinuous problems can be calculated using the discontinuous deformation analysis (DDA) method by solving time steps, and this method is suitable for simulating the seismic dynamic response of engineering rock mass structures. However, the boundary setting must be carefully analyzed. In this paper, four boundary settings for the DDA method are investigated. First, the contributions to the DDA equations for nonreflecting boundaries (including the viscous boundary and the viscoelastic boundary) are deduced based on the Newmark method. Second, a free‐field boundary is introduced in the DDA method with boundary grid generation and coupling calculation algorithms to accurately simulate external source wave motion, such as earthquakes. Third, seismic input boundary treatments are intensively examined, and the force input method is introduced based on nonreflecting boundaries. Finally, the static‐dynamic unified boundary is implemented to ensure consistent boundary transformation. The boundary setting method in the DDA method is discussed, and the suggested treatments are used to analyze the seismic dynamic response of underground caverns. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

2.
边界条件对土层粘弹性地震反应的影响   总被引:2,自引:0,他引:2  
刘汉龙  高玉峰  朱伟  费康 《岩土力学》2001,22(4):408-412
边界条件是土层粘弹性地震反应必须考虑的问题。对于地表面,一般作为自由边界,而对于地下深处基岩面,边界条件有两种处理方法:一是静止边界,二是非静止边界的。首先给出了在时间域内单层地基一维土层粘弹性地震反应封闭解析解,然后从理论上分析边界条件对土层粘弹性地震反应的影响,认为采用静止边界进行土层地震反应计算得到的地震动加速度时程符合天然地震记录形式,而采用非静止边界进行土层地震反应计算得到的地震动加速度时程初始阶段出现较大幅度的振荡现象,随着土层厚度增加,这种振荡现象逐渐减弱直至消失。当土层阻尼较小或剪切波速增长时,振荡现象加剧;当土层阻尼比较大或剪切波速减小时,振荡现象减弱。此外,采用非静止边界进行土层地震反应计算,表现出随着土层厚度的增加,土层对基岩输入地震动的放大作用逐步转化为吸收作用。  相似文献   

3.
In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

4.
夏成志  石振明  李博  郑鸿超  刘毛毛 《地球科学》2022,47(12):4469-4483
开发一种改进损伤框架的粒子流算法,被称为核断裂的光滑粒子流法(kernel?broken smoothed particle hydrodynamics,KBSPH),用于模拟地震条件下岩质边坡的裂纹扩展和变形破坏过程.在KBSPH中,提出一种改进的损伤框架,通过引入断裂标志来改进损伤粒子的核函数,使损伤粒子的虚拟应力键直接断裂,裂纹在断裂的应力键间生成,从而模拟岩石的裂纹扩展过程.在地震边界上采用了双层边界,将动力输入边界与黏滞边界分离.首先通过薄板振动实验验证KBSPH的动力特性.其次以单裂隙岩体单轴压缩试验验证KBSPH的断裂力学特性.最后模拟地震条件下多节理岩质边坡中裂纹扩展过程和动力响应.薄板振动实验验证了KBSPH的动力特性的准确性.单裂隙岩体单轴压缩试验,证明了KBSPH可以正确模拟预制裂隙尖端的翼型裂纹.通过对比以往数值模拟方法和现场案例,表明KBSPH正确揭示了加速度放大效应以及地震条件下岩质边坡的裂纹扩展过程.KBSPH避免了传统算法的网格畸变,损伤粒子应力分量重新分配的问题,降低了编程难度,提高了运行速率,可为SPH在地震条件下岩石力学中的应用和理解岩石断裂机理提供一定的参考.   相似文献   

5.
首先,介绍了基于OpenSees独立开发的一套用于挡土墙-土地震反应相互作用有限元分析计算软件RW_2DPS.据此建立了俯斜式混凝土重力挡土墙-土强震相互作用有限元模型.模型中,引入非线性有限元计算方法,选用多屈服面弹塑性本构模型模拟砂土的动力属性,应用零长度接触单元模拟墙与土体之间的接触特性,且采用一致耗能阻尼边界与速度边界条件.最后,输入随机地震动,进行挡土墙-土强震反应分析,并重点探讨墙背地震土压力和水平地震惯性力沿挡土墙高度分布规律.结果表明,墙背动土压力峰值出现在距挡土墙底约1/3墙高处;挡土墙背加速度具有放大效应,加速度峰值出现在挡土墙顶部;不同地震动作用下,加速度放大系数沿墙高分布规律不同,动土压力沿墙高变化规律基本一致.  相似文献   

6.
Rock slope failure is a complex process that usually involves both opening/sliding along pre‐existing discontinuities as well as fracturing of intact rock bridges. Discontinuity persistence is an important factor governing rock slope instabilities. However, traditional slope failure analysis assumes persistent discontinuities, and rock slope fails along a predefined persistent continuous potential failure surface because of the limitations of the analysis tools. This paper proposes the numerical manifold method (NMM) incorporated with a Mohr–Coulomb criterion‐based fracturing algorithm to simulate the progressive failure of rock slopes with non‐persistent joints. Detailed fracturing algorithm is first presented. Then, the NMM enabling fracturing is calibrated through simulating an edge‐cracked plate and the Brazilian test. Lastly, the developed code is applied to investigate the failure process of rock slopes involving non‐persistent joints. Numerical results indicate that the proposed method can capture the opening/sliding along existing discontinuities, the fracturing in intact rock bridges and the final kinematic release. Progressive slope failure is well exhibited. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

7.
陈国兴  李磊  丁杰发  赵凯 《岩土力学》2020,41(9):3056-3065
地震基岩深度和土体动力本构模型的选取对核岛场地地震效应评价结果的合理性具有重要影响。以拟建某沿海核电厂深度470 m沉积土夹火山岩层场地的3个钻孔剖面为研究对象,采用一维等效线性波传分析(ELA)法、基于Matasovic本构模型和Davidenkov-Chen-Zhao(DCZ)本构模型的一维非线性分析(NLA)法,选取不同剪切波速的5个岩土层作为地震基岩,研究了输入地震动特性、地震基岩深度和土体动力本构模型的选取对巨厚沉积土夹火山岩层场地非线性地震反应特性的影响。结果表明:(1)以浅层硬岩夹层或深部土层作为地震基岩,NLA法计算的5%阻尼比的地表谱加速度SA的短周期部分较之ELA法的计算值大,但两者计算的地表SA谱的长周期部分几乎一致;(2)基于Matasovic模型和DCZ模型的NLA法计算的地表SA谱谱形和峰值加速度随深度的变化趋势基本一致;(3)从NLA法计算的地表峰值加速度和累积绝对速度而言,以剪切波速约2 500 m/s的浅层硬岩夹层作为地震基岩是适宜的。  相似文献   

8.
三维黏弹性人工边界地震波动输入方法研究   总被引:3,自引:0,他引:3  
对三维黏弹性人工边界给出了一种实用的地震波动输入方法。用等效边界力法进行地震波动输入时,把入射地震波假设为竖直向上入射的平面波,使施加在人工边界上的总荷载可以由输入波位移和速度的关系式完全表达,使地震动输入更容易在有限元计算程序中实现。解决了三维黏弹性人工边界采用球面波进行地震波动输入会出现尚待确定的未知项,不易实现地震波动的输入的困难。数值算例表明,该地震波动输入方法是十分有效的。  相似文献   

9.
柳国环  练继建  于通顺 《岩土力学》2014,35(9):2651-2658
近海风力发电结构是由无限域地基、基础及上部结构组成的体系,有限域透射人工边界和地震动输入模式直接影响结构体系的动力响应。首先,建立透射人工边界-地基-基础-塔筒结构和固定边界-地基-基础-塔筒结构的有限元模型;然后,简要澄清了透射边界地震动输入之所以采用外力而不采取地震动物理量(加速度、速度和位移)的原因所在;最后,对给出的3种作用模式:(1)固定边界地震动输入,(2)只考虑外源输入波作用,(3)同时考虑外源输入波和内源振动的散射波作用,分别进行计算分析和比较。通过数值计算并结合理论分析,(1)澄清了对透射边界采用外力进行地震激励输入的方法的原因,这种输入方法合理可靠,符合实际情况;(2)在地震波等效荷载的生成中,针对圆形横截面地基提出了一种阻尼力、刚度力的便捷性生成方法,大大减小了工作量;(3)从自身振动特点及阻尼设置角度,解释了地震作用下采用透射边界比采用固定边界时风电结构响应减小的原因;(4)指出了地震作用下透射边界地基模型在结构动力响应与倒塌计算中内源振动反射的不容忽视性。  相似文献   

10.
The primary goals of seismic interpretation and quantification are to understand and define reservoir architecture and the distribution of petrophysical properties. Since seismic interpretation is associated with major uncertainties, outcrop analogues are used to support and improve the resulting conceptual models. In this study, the Miocene carbonates of Cerro de la Molata (Las Negras, south‐east Spain) have been selected as an outcrop analogue. The heterogeneous carbonate rocks of the Cerro de la Molata Platform were formed by a variety of carbonate‐producing factories, resulting in various platform morphologies and a wide range of physical properties. Based on textural (thin sections) and petrophysical (porosity, density, carbonate content and acoustic properties) analyses of the sediments, eleven individual facies types were determined. The data were used to produce synthetic seismic profiles of the outcrop. The profiles demonstrate that the spatial distribution of the facies and the linked petrophysical properties are of key importance in the appearance of the synthetic seismic sections. They reveal that carbonate factory and facies‐specific reflection patterns are determined by porosity contrasts, diagenetic modifications and the input of non‐carbonate sediment. The reflectors of the seismograms created with high‐frequency wavelets are coherent with the spatial distribution of the predefined facies within the depositional sequences. The synthetic seismograms resulting from convolution with lower frequency wavelets do not show these details – the major reflectors coincide with: (i) the boundary between the volcanic basement and the overlying carbonates; (ii) the platform geometries related to changes in carbonate factories, thus sequence boundaries; and (iii) diagenetic zones. Changes in seismic response related to diagenesis, switching carbonate producers and linked platform geometries are important findings that need to be considered when interpreting seismic data sets.  相似文献   

11.
综合管廊由于埋深较浅,在抗震分析中应考虑Rayleigh波的作用,为研究Rayleigh波与底部地震加速度共同作用下综合管廊的动力响应特征,建立双仓的综合管廊三维动力有限元数值模型,土体采用考虑滞回环特性的高级本构模型(HSS模型),通过边界上多次脉冲荷载生成Rayleigh波,模型底部横向分别作用Upland波、Kobe波、Taft波,并与仅考虑底部横向作用的常规时程分析进行对比。研究表明:综合管廊结构的横向动力响应主要受横向地震波影响,结构纵向动力响应受沿其轴向入射的Rayleigh波影响相对较大;采用Rayleigh波+底部地震波的输入方法比单独底部地震波输入得到的结构动力响应整体上要更显著一些;输入不同的地震加速度时程,管廊动力反应规律相似,但综合管廊结构影响大小有差异,可见底部地震波与地表Rayleigh波作用的匹配程度对结构动力响应结果有一定影响。研究成果可供地下综合管廊结构地震动力响应精细化数值分析及抗震设计参考。  相似文献   

12.
根据三维地震地质模型对地震数据进行模拟是从勘探到生产的周期内决策过程中的一个不可或缺的组成部分。虽然对于在储层内的动力过程和地震地质的模型表述已经取得很大进展,但如何从这些模型得到地震数据的精确模拟仍面临很多挑战。通常是在地球模型范围内根据物性用一维褶积方法来模拟地震数据。然而这个过程一般不考虑地震勘探布局和盖层对地震信号的影响。我们审视了为什么这些因素会制约三维地球模型的有效性,并考虑了为什么需要把盖层和地震勘探布局对三维覆盖和分辨率的影响加进模拟过程之中。我们提出了一种新方法,把建立物性模型和一种新的地震模拟技术结合起来,给出一个工作流程;利用这个流程,勘探工作者可以很快模拟出三维的PSDM数据,这些数据加进了盖层和地震勘探布局对覆盖及分辨率的影响。我们利用从远离挪威海岸的一个油田得到的数据,在考虑覆盖和分辨率效应的地震数据模拟之前,对岩石物性做了一些扰动,然后进行地震数据模拟,以此来说明如何可以用这种方法提高三维地球模型的精确性和增进我们对储层的了解。  相似文献   

13.
考虑地震动的随机性,利用复反应分析技术,采用随机地震反应计算方法对某一特大型桥梁群桩基础与土动力相互作用效应进行了数值试验研究。将土与群桩体系视为一个整体进行有限元离散,采用等效线性化方法考虑土体的动力非线性性能。将桥墩-群桩-土相互作用体系与自由场随机地震反应进行比较,结果表明:相互作用效应的影响与桩土模量差异以及土体与群桩基础距离有关;软土层剪应变水平及分布发生了显著变化,群桩基础两侧附近土域剪应变呈现明显的弧形分布;地表及浅层土体最大地震加速度反应有所减小,但覆盖层中下部土体加速度反应峰值明显增大,增幅达5 %~30 %左右;此外,地震地面运动的频谱成分存在显著差别。桥梁桩基础抗震设计中应充分考虑桥梁结构-群桩-土相互作用效应。  相似文献   

14.
A major challenge in seepage analysis is to locate the phreatic surface in an unconfined aquifer. The phreatic surface is unknown and assumed as a discontinuity separating the seepage domain into dry and wet parts, thus should be determined iteratively with special schemes. In this study, we systematically developed a new numerical manifold method (NMM) model for unconfined seepage analysis. The NMM is a general numerical method for modeling continuous and discontinuous deformation in a unified mathematical form. The novelty of our NMM model is rooted in the NMM two‐cover‐mesh system: the mathematical covers are fixed and the physical covers are adjusted with iterations to account for the discontinuity feature of the phreatic surface. We developed an energy‐work seepage model, which accommodates flexible approaches for boundary conditions and provides a form consistent with that in mechanical analysis with clarified physical meaning of the potential energy. In the framework of this energy‐work seepage model, we proposed a physical concept model (a pipe model) for constructing the penalty function used in the penalty method to uniformly deal with Dirichlet, Neumann, and material boundaries. The new NMM model was applied to study four example problems of unconfined seepage with varying geometric shape, boundary conditions, and material domains. The comparison of our simulation results to those of existing numerical models for these examples indicates that our NMM model can achieve a high accuracy and faster convergence speed with relatively coarse meshes. This NMM seepage model will be a key component of our future coupled hydro‐mechanical NMM model. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

15.
This paper proposes a numerical model for jointed rock masses within the 3‐D numerical manifold method (NMM) framework equipped with a customized contact algorithm. The strength of rock sample containing a few sets of discontinuities is first investigated. The results of models with simple geometries are compared with the available analytical solutions to verify the developed computer code, whereas models with complex geometries are simulated to better understand the fundamental behavior and failure mechanism of jointed rock mass. Furthermore, the stability of jointed rock mass in an underground excavation is studied, where rock failure process is determined by the 3‐D NMM simulation. The simulation results provide valuable guidance on excavation process design and stabilization design in rock engineering practice. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

16.
Seismic stability analysis of a rock block is an important issue in the field of rock mechanics and rock engineering. To solve this problem, the block theory and the Newmark method are combined, and a general method for seismic response of a rock block is introduced. First, the formation method of a three-dimensional rock block, which includes establishing a topological relationship among the block-polygon face-edge-vertex and dividing a complex block into convex subblocks, is presented, and the simplex integration method is employed to calculate the volume of a rock block and the areas of the polygonal faces. Second, the assumptions and algorithms of the general method for seismic response of a rock block are detailed. The dynamic analysis is carried out in time step. In each time step, the key technologies including analysis of the seismic force and motion mode, trial of the incremental displacement, check for the block entrance, and update of the motion parameters are performed in order. Last, two verification examples and three application examples including a wedge, a rock block of slope, and a combined rock block are used to analyze the correctness and practicality of the general method. The results show that, for a given ground motion record, the Newmark program can effectively simulate the dynamic response process of the motion mode, velocity, and safety factor of the rock block, and the permanent displacement under the earthquake action is obtained, which provides a quantitative parameter to evaluate the dynamic stability of a rock block.  相似文献   

17.
为研究注浆措施在地震区隧道洞口段的减震效果,基于围岩-结构三维模型,改变注浆的范围及参数,采用有限差分法软件FLAC-3D,分析了某隧道洞口段的动力响应。采用人工拟合地震波、黏性边界条件,输入三个方向的地震波,模拟了原设计、全环间隔注浆措施、局部加强注浆措施三种工况。研究分析了各工况的内力、位移、周边围岩塑性区等,并与全环间隔注浆、局部加强注浆两种减震措施的效果作了对比,得出了地震区隧道洞口段的动力响应规律、减震措施的有效范围及全环间隔注浆措施能起到阻尼减震的效果、局部加强注浆措施能改善衬砌结构内力分布的作用等结论。对高烈度地震区隧道洞口段抗减震措施研究有一定的借鉴意义。  相似文献   

18.
戴永浩  陈卫忠  谭贤君  黄胜 《岩土力学》2010,31(Z2):236-241
将盾构隧道管片视为分段曲梁、管片之间为铰接关系、管片与地层之间为非线性接触关系,边界采用无限元吸收地震波的反向与折射,地层服从Mohr-Coulomb弹塑性本构模型,地震波采用EI-Centro波,采用地层-结构法,应用隐式动力有限元分析地震荷载作用下管片轴力与弯矩分布情况。根据有限元计算结果,结合结构力学方法,校核管片配筋是否合理。与荷载结构法相比,其方法可模拟围岩与管片、管片与管片之间的相互关系,并可考虑锚杆等支护结构对管片受力的影响,因而,计算结果更为精确。研究为盾构隧道管片配筋优化设计,降低隧道支护成本,提供了一种管片配筋校核计算方法与设计依据。  相似文献   

19.
Practical civil engineering problems are usually formulated in an infinite half-space domain, and a selected finite domain is required to analyze the dynamic responses of a fluid-saturated porous medium by the finite element method (FEM). Devising a method to deal with the boundaries of the finite domain is the key issue for this open system. In this paper, a two-dimensional spring–dashpot artificial boundary (SDAB) for transient analysis in a fluid-saturated porous media is developed. Based on Biot’s dynamic theory of fluid-saturated porous media, the normal and tangential boundary stress formulae are deduced for out-going cylindrical body waves. The boundary stress is proportional to displacement and velocity, thus continuously distributed dashpots and springs can be placed on the artificial boundaries in the normal and tangential directions to simulate the energy absorption of the infinite media outside of the finite domain for the interior distributed source problems. In this paper, the input seismic motion can be realized by applying an equivalent load on the SDAB for the seismic scattering problems of exterior distributed sources. Numerical examples are given and the analyzed results show that the SDAB and the method of wave motion input have good stability and acceptable accuracy.  相似文献   

20.
姚虞  王睿  刘天云  张建民 《岩土力学》2018,39(6):2259-2266
在多点输入方法和等价黏弹性模型的基础上,采用半解析的波函数组合法实现了面板堆石坝在非一致输入下的动力计算。从频域散射角度对计算方法的合理性进行了验证。对比了非一致输入与一致输入下高面板堆石坝的动力响应,发现在采用基岩自由表面点的振动过程相同的对比标准下,非一致输入的整体动力响应较小;在面板坝防渗系统相对薄弱的止水结构附近,非一致输入的动拉应力最大值比一致输入的结果更大;非一致输入下大坝动力响应值的分布相比一致性输入呈现出中间小、周边大的特性。并基于波动理论进一步分析了不同种类地震波入射角度对高面板堆石坝动力响应的影响,揭示了P波、SV波和SH波入射下高面板堆石坝动力响应规律:随着入射角的增大,SH波入射时,高面板坝动力响应强度基本不变; SV波入射时,存在一个临界角,当入射角在临界角左右时高面板坝动力响应强度急剧增大和减小,之前基本不变,之后一直减小;P波入射时,存在一个特征角度,在入射角小于特征角度时高面板坝动力响应强度基本不变,大于特征角度时高面板坝动力响应强度减小。  相似文献   

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

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