共查询到19条相似文献，搜索用时 625 毫秒
1.
针对光滑粒子流体动力学方法（SPH）在滑坡模拟中建立粒子模型的难题,提出了基于地理信息系统（GIS）栅格数据的粒子排列与插入方法。根据该方法,建立了滑坡SPH粒子模型及相关粒子生成程序,进一步以结合摩尔-库仑破坏准则的SPH宾汉流体模型为核心,实现了运用SPH方法模拟滑坡破坏后三维运动的过程。该SPH模型在对唐家山滑坡的模拟中得到了验证,并预测了金坪子滑坡破坏后的影响范围。结果表明：基于GIS空间数据的滑坡SPH粒子模型具有可行性与良好的适用性。以GIS数据库为基础,开展滑坡灾害的模拟研究,将大大提高对滑坡等地质灾害的仿真分析,为滑坡灾害的预测与防治提供参考。 相似文献
2.
土体的大变形流滑导致了许多地质灾害的发生,对人们的生命财产安全构成了极大的威胁,因此越来越多的研究开始关注土体的大变形流滑特性。其中,光滑粒子流体动力学（SPH）方法是常用的模拟方法之一,但SPH方法的粒子特性导致其计算时间过长,影响了在工程地质领域的进一步应用。对此,本研究基于SPH方法的基本原理、非牛顿流体理论和等效黏度概念,提出了适用于土体大变形流滑分析的三维SPH仿真模型。结合OpenMP并行计算原理,实现了SPH算法的并行优化。在此基础上,对土体流滑模型试验进行了二维和三维分析,得到了滑动距离、滑动冲击力和冲击力峰值等动力学参数,分析了计算维数和边界条件对流滑特性的影响机制。通过不同线程数下计算时间的对比,获得了计算效率随线程数的变化规律。结果证明了本文的OpenMP并行优化具有较高的计算效率,显著降低了三维SPH模拟的计算耗时,对工程地质数值方法的效率提升具有重要的借鉴意义。 相似文献
3.
溃坝洪水演进模拟的准确性是制约水库洪水预演有效性的关键。基于光滑粒子流体动力学(Smoothed Particle Hydrodynamics, SPH)方法提出了适用于溃坝洪水演进分析的数值模拟方法。通过设置溃口粒子与粒子库,基于黎曼不变量对SPH粒子状态进行修正,构建施加边界条件的改进SPH溃坝洪水演进模型,将SPH瞬时全溃整体模型转换为考虑溃口水流变化的入流边界模型,实现SPH方法与溃口计算模型的耦合。以Malpasset溃坝事件为例,检验了该模型计算溃坝洪水的精度,结果表明该模型精度相对较高,与实测值吻合较好;应用该模型模拟了某水库溃坝洪水演进预演过程,评估其对下游输水干渠及交叉建筑物排水倒虹吸的洪水冲击风险,结果表明在上游水库遭遇超标准洪水漫顶溃坝工况下,洪水演进至排水倒虹吸处的最大洪水位未超过校核洪水位。改进SPH模型精度高,可靠性强,与溃口计算模型耦合性好,可作为溃坝洪水演进模拟的通用手段之一。 相似文献
4.
考虑海冰热力因素对其厚度、密集度的影响,在光滑质点流体动力学(SPH)基础上发展了一个海冰热力-动力数值模式。该模式既解决了传统欧拉有限差分法和质点网格法存在的数值扩散问题,同时弥补了光滑质点动力学海冰动力模式未考虑热力因素的不足,具有精确模拟冰缘线运动、计算精度高等优点。首先介绍了光滑质点流体动力学的基本原理,并对海冰生消的热力因素进行了分析,将影响冰厚和密集度的热力因素引入到光滑质点流体动力学的海冰动力模式中,得到该热力-动力模式的控制方程。应用该数值模式对渤海海冰进行了48 h数值模拟,得到了海冰厚度和速度矢量的分布规律;对JZ20-2海域的海冰厚度、冰内温度场分布以及热力因素的变化特性进行了讨论。数值模拟结果表明,该数值模式能够很好地适用于渤海海冰数值模拟,是一种有效的海冰数值模拟方法。 相似文献
5.
海冰动力学数值方法研究进展 总被引：1，自引：0，他引：1
在海冰动力学数值模拟和预测研究中,人们将海冰视为连续介质分别建立了欧拉坐标下的有限差分（FD）方法、拉格朗日坐标下的光滑质点流体动力学（SPH）方法、欧拉和拉格朗日坐标相结合的质点网格法（PIC）,近年来又发展了基于非连续介质的颗粒流（GF）方法。对以上几种海冰动力学数值方法的特点和适用性进行了讨论,结果表明：FD、PIC和SPH方法可适用于中长期海冰动力学数值模拟,但SPH方法的计算效率需进一步提高;GF方法在不同尺度下的海冰动力学数值模拟中的计算精度均有很强的适用性,但目前较适用于小尺度下海冰动力学基本特性的数值试验研究,计算时效还不能满足实际海冰数值模拟和预测的要求。为进一步提高海冰动力学模拟的精度和适用性,在不同时空尺度下分别发展与其相适应的数值方法是必要的。 相似文献
6.
光滑粒子流体动力学(SPH)由于无需网格生成和拉格朗日特性,对求解带有自由表面和大变形的力学问题有优势。但是该方法存在计算精度不高,计算效率较低等缺点。为此重点对SPH方法的精度提高进行研究。介绍了传统算法的基本公式,根据误差分析指出该算法精度不高的原因,提出了SPH二阶精度算法。通过精度验证分析,证明了该方法的精度的确能够达到二阶。通过二维计算实例,给出传统方法和二阶方法在粒子均匀分布和非均匀分布时函数值以及函数的一、二阶导数的误差分布,证明二阶算法能够克服传统算法的一些缺点,且计算精度有较大提高。 相似文献
7.
《岩土力学》2017,(4):1226-1232
采用光滑粒子流体动力学(SPH)方法和离散单元法(DEM),并基于达西渗透试验原理,提出一种SPH-DEM耦合算法用于处理宏观尺度下离散体与流体之间的相互作用问题,并基于该流-固耦合方法采用FORTRAN语言建立滑坡涌浪数学模型。模拟了块体滑坡问题,并与试验数据进行对比以验证该耦合方法的有效性。模拟分析了离散体滑坡产生涌浪以及涌浪传播的过程,分析比较了数值模拟得到的最大涌浪高度与经验公式结果间的关系。研究表明:采用该耦合方法模拟水下块体滑坡问题得到的结果与试验数据吻合良好;模拟离散体滑坡涌浪问题时,较清晰地反映了滑坡体对水的排挤到涌浪的产生过程、水在滑坡体中的渗透过程和滑坡体与水耦合作用发生变形的过程;不同公式分析方法得到的最大涌浪高度之间均有一定差别,算例中滑坡体受重力作用沿倾斜坡面滑动入水过程更接近于潘家铮法的垂直运动模式,因此,模拟结果与之最为接近。 相似文献
8.
本构模型是描述泥石流流变特性的关键,也是决定其动力过程数值模拟准确性的核心问题之一。泥石流流体属多相混合物,现有的研究已证实其存在剪切增稠或剪切变稀的现象,传统基于Bingham及Cross线性本构关系的数值模型难以准确描述泥石流流变特性。文中探讨了Bingham模型在低剪应变率下的数值发散问题,在光滑粒子流体动力学（SPH）方法框架上建立了整合Herschel-Bulkley-Papanastasiou(HBP）本构关系的稀性泥石流动力过程三维数值模型。相比传统基于浅水波假设的二维数值模型,所述方法从三维尺度建立SPH形式下的泥石流浆体纳维?斯托克斯方程并进行数值求解,可获取泥石流速度场时空分布及堆积形态,同时采用HBP本构关系描述泥石流流变特性,能在确保数值收敛的前提下反映泥石流流体在塑性屈服过渡段及大变形状态下应力?应变的非线性变化。为验证提出方法的合理性,结合小型模型槽实验观测进行了对比,结果表明数值模拟与实测结果基本吻合。 相似文献
9.
为了研究热源温度和外界水压对缓冲层中水-热迁移规律的影响,以GMZ膨润土为例,从基于势能的非饱和土的水-热迁移控制方程出发,考虑了蒸发效应的影响,得到了水-热耦合的方程组,采用改进的光滑粒子流体动力学(SPH)算法,能够对每一处土体根据不同时刻的不同状态实时更新计算参数,得到参数变化的水-热耦合解。计算结果表明：土的物理性质参数与土体的温度和饱和度密切相关,是否考虑这些参数的变化会对计算结果产生较大影响;核废料释放的热量能够在较短的时间内扩散到外边界,水分迁移的速度则相对慢很多;缓冲层温度的升高会加快水分的迁移速度,外界水压对温度的分布则影响较小。 相似文献
10.
在1:60溢流坝表、中泄水孔断面模型上对密松水电站泄水建筑物的布置方案进行了泄流能力、流速流态、压力分布等相关水力学问题的试验,并通过十多个方案的比较优化,提出了基本可行的泄水建筑物体型和消力池型式。研究成果表明:表孔采用合适的宽尾墩形式,对增加消能率、缩短消力池池长有利,但无法消除宽尾墩形成的水翅击打两侧导墙现象;表孔取消宽尾墩,采用跌坎式消力池方案,可有效降低底流消能时的池底板流速,同时消除宽尾墩形成的水翅击打两侧导墙现象,达到较满意的消能效果。本文中十数个消能方案的研究成果,对类似工程的设计和科研均有参考价值。 相似文献
11.
A Lagrangian particle‐based method, smooth particle hydrodynamics (SPH), is used in this paper to model the flow of self‐compacting concretes (SCC) with or without short steel fibres. An incompressible SPH method is presented to simulate the flow of such non‐Newtonian fluids whose behaviour is described by a Bingham‐type model, in which the kink in the shear stress vs shear strain rate diagram is first appropriately smoothed out. The viscosity of the SCC is predicted from the measured viscosity of the paste using micromechanical models in which the second phase aggregates are treated as rigid spheres and the short steel fibres as slender rigid bodies. The basic equations solved in the SPH are the incompressible mass conservation and Navier–Stokes equations. The solution procedure uses prediction–correction fractional steps with the temporal velocity field integrated forward in time without enforcing incompressibility in the prediction step. The resulting temporal velocity field is then implicitly projected on to a divergence‐free space to satisfy incompressibility through a pressure Poisson equation derived from an approximate pressure projection. The results of the numerical simulation are benchmarked against actual slump tests carried out in the laboratory. The numerical results are in excellent agreement with test results, thus demonstrating the capability of SPH and a proper rheological model to predict SCC flow and mould‐filling behaviour. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
12.
Weijie Zhang Kenichi Maeda Hiroshi Saito Zhaoqing Li Yu Huang 《Acta Geotechnica》2016,11(6):1401-1418
For seepage failures of dike due to water level-up and rainfall, surface infiltration and strength change induced by suction reduction are important factors; thus, numerical analysis should consider the coupling of water and soil, as well as the effect of saturation to obtain more precise failure mechanism. Based on the advanced smoothed particle hydrodynamics (SPH) method, this work proposed a two-phase-coupled SPH model in coordination with a novel constitutive model for unsaturated soils. Then, a triaxial compression test is simulated to check the applicability of the SPH method on the soil phase. After that, the failure test of a dike due to water level-up is discretized and simulated, from which the seepage process, the distribution of maximum shear strain, the slip surface, and pore water pressure are obtained. The two-phase-coupled SPH model is also applied to a slope failure test of heavy rainfall, and the results are compared to the model test. Finally, a dike failure test due to rainfall is analyzed using the proposed SPH model to reproduce the surface infiltration and suction reduction. The proposed SPH model provides several insights of seepage failures and can be a helpful tool for the analysis of dike failures induced by water level-up and rainfall. 相似文献
13.
Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
14.
Smoothed particle hydrodynamics (SPH) is a versatile technique which can be applied to single and multiphase flow through porous media. The versatility of SPH is offset by its computational expense which limits the practicability of SPH for large problems involving low Reynolds number flow. A parallel pore-scale numerical model based on SPH is described for modeling flow phenomena in porous media. Aspects of SPH which complicate parallelization are emphasized. The speed of the method is demonstrated to be proportional to the number of processors for test cases where load balance was achieved. The parallel algorithm permits the application of SPH to more complicated porous media problems than previously considered. For such problems, best performance is achieved when several soil grains are simulated by each processor. Finally, future applications of the method and possible extensions are discussed. 相似文献
15.
This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow. 相似文献
16.
Natural Hazards - The paper presented a smoothed particle hydrodynamics (SPH) method to study the three-dimensional (3D) tidal bore scenarios. The SPH method is a mesh-free particle modeling... 相似文献
17.
The present study focuses on the emergency response measures for landslide dams. This work presents a series of centrifuge model tests conducted on the draining processes of barrier dams that are based on the grain composition of the Tangjiashan landslide dam. The effects of diversion channels with trapezoid, triangular, and compound sections on the discharge, process, and size of the residual dam are discussed. The characteristics of the flow erosion during the process of discharge in the different channels are analyzed based on hydrodynamics. The results suggest that a diversion channel with a compound section has a higher initial discharge efficiency and lower peak flow, and the flow process curve corresponds to a “chunky-type.” The draining from this type of a diversion channel could clearly reduce the flood pressure of the river downstream and make the entire process smoother. Thus, the excavation of a diversion channel with a compound section is an efficient and safe method for landslide dam emergency mitigation. 相似文献
18.
To overcome the disadvantages of traditional flow analysis methods for liquefied soils that exhibit fluidization and large
deformation characteristics, Smoothed particle hydrodynamics (SPH) is adopted in this study to analyze the flow processes
of liquefied soils. Bingham model with the use of the Mohr–Coulomb yield criterion, the concepts of equivalent Newtonian viscosity,
and the Verlet neighbor list method are introduced into the framework of SPH to build an algorithm for the analysis of flowing
liquefied soils. This modeling involves a simulation of physical model test of flowing liquefied soils that can be compared
with numerical results. In addition, a shaking table test is selected from the literature for SPH analysis to verify the validation
of the SPH method and extend its applications. The SPH simulation can reproduce the flow processes of liquefied soils and
constrain estimates of the horizontal displacement, vertical displacement, and velocity of soils after liquefaction. According
to the dynamic behaviors of the materials involved, designs can be implemented to improve the seismic safety of structures. 相似文献
19.
The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献