共查询到20条相似文献，搜索用时 109 毫秒
1.
基于DDA方法一种流-固耦合模型的建立及裂隙体渗流场分析和应用 总被引：1，自引：0，他引：1
以非连续变形分析方法（DDA）为基础并采用稳态流体计算方法将二者结合进行裂隙岩体流-固耦合分析。利用DDA方法生成裂隙岩体模型,在此基础上采用矩阵搜索等方法形成新的裂隙水通网络模型。采用稳态迭代算法和立方定律求得裂隙水压力,并把裂隙水压力作为线载荷施加到块体边界,在DDA算法中每个迭代步完成后更新裂隙开度和水压值,与DDA算法结合研究裂隙水与块体之间相互作用关系。利用以上裂隙岩体流-固耦合计算方法研究了某水封油库开挖和运行过程洞室围岩流量和密封性,为该工程预测水封效果提供了有益的主要依据,也是国内首次采用DDA方法做大型工程的流-固耦合模型分析。 相似文献
2.
To investigate the movement mechanism of debris flow, a two‐dimensional, two‐phase, depth‐integrated model is introduced. The model uses Mohr‐Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second‐order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement. 相似文献
3.
One dimensional solutions for the classic critical upward seepage gradient/quick condition and the time rate of consolidation problems are obtained using coupled routines for the finite volume method (FVM) and discrete element method (DEM), and the results compared with the analytical solutions. The two phase flow in a system composed of fluid and solid is simulated with the fluid phase modeled by solving the averaged Navier–Stokes equation using the FVM and the solid phase is modeled using the DEM. A framework is described for the coupling of two open source computer codes: YADE-OpenDEM for the discrete element method and OpenFOAM for the computational fluid dynamics. The particle–fluid interaction is quantified using a semi-empirical relationship proposed by Ergun [12]. The two classical verification problems are used to explore issues encountered when using coupled flow DEM codes, namely, the appropriate time step size for both the fluid and mechanical solution processes, the choice of the viscous damping coefficient, and the number of solid particles per finite fluid volume. 相似文献
4.
In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
5.
非连续变形分析（DDA）是一种隐式求解的动力学计算方法,且采用在块体界面加减刚硬弹簧的方式来满足块体界面无张拉和无嵌入的接触准则,其中时间步长和弹簧刚度两个物理量的取值直接影响DDA的计算结果。基于对DDA时间步和弹簧刚度在程序运行过程中的调整策略和块体接触的简化力学概念模型,研究了惯性力在DDA收敛求解中的作用过程。采用数值模拟试验对自由落体和斜面单滑块模型在3种力学状态下的相关力学问题进行了数值模拟研究,通过对自由落体运动的模拟,研究了时间步长单一因素对计算结果的影响规律,并初步确定了时间步长的合理取值区间。在此基础上,采用斜面单滑块模型,研究了时间步长和弹簧刚度对计算结果的共同影响,确定了不同时间步长条件下弹簧刚度的合理取值区间。研究成果表明,合适的时间步长和弹簧刚度的取值组合构成一个单连通参数取值域,当时间步和弹簧刚度的取值组合位于此“域”范围内时,DDA的计算结果是合理的。 相似文献
6.
A coupled 3‐dimensional bonded discrete element and lattice Boltzmann method for fluid‐solid coupling in cohesive geomaterials 
下载免费PDF全文
下载免费PDF全文 This paper presents a 3D bonded discrete element and lattice Boltzmann method for resolving the fluid‐solid interaction involving complicated fluid‐particle coupling in geomaterials. In the coupled technique, the solid material is treated as an assembly of bonded and/or granular particles. A bond model accounting for strain softening in normal contact is incorporated into the discrete element method to simulate the mechanical behaviour of geomaterials, whilst the fluid flow is solved by the lattice Boltzmann method based on kinetic theory and statistical mechanics. To provide a bridge between theory and application, a 3D algorithm of immersed moving boundary scheme was proposed for resolving fluid‐particle interaction. To demonstrate the applicability and accuracy of this coupled method, a benchmark called quicksand, in which particles become fluidised under the driving of upward fluid flow, is first carried out. The critical hydraulic gradient obtained from the numerical results matches the theoretical value. Then, numerical investigation of the performance of granular filters generated according to the well‐acknowledged design criteria is given. It is found that the proposed 3D technique is promising, and the instantaneous migration of the protected soils can be readily observed. Numerical results prove that the filters which comply with the design criteria can effectively alleviate or eliminate the appearance of particle erosion in dams. 相似文献
7.
《Geomechanics and Geoengineering》2013,8(4):293-302
This paper presents a comparative study of two methods, Sarma's method and the discontinuous deformation analysis (DDA), for rock slope stability analysis. The comparison concerns the stability analysis of two classic rock slopes. The study shows that the DDA, which accounts for the block kinematics, provides a very different factor of safety as compared with Sarma's method. More realistic reaction forces around each rock block can be obtained by the DDA, including the thrust forces between rock blocks and the forces between the base and the blocks. The DDA's result shows two possible directions for the relative movement between two contiguous blocks at the initiation of slope failure. It also indicates that the limit equilibrium condition may not occur along the interfaces of rock blocks at the initiation of slope failure. The determination of realistic interaction forces around each block will be very important in rock slope stability analysis if nonlinear failure criteria are considered. 相似文献
8.
A Boundary Element based Discontinuous Deformation Analysis (BE‐DDA) method is developed by implementing the improved dual reciprocity boundary element method into the open close iterations based DDA. This newly developed BE‐DDA is capable of simulating both the deformation and movement of blocks in a blocky system. Based on geometry updating, it adopts an incremental dynamic formulation taking into consideration initial stresses and dealing with external concentrated and contact forces conveniently. The boundaries of each block in the discrete blocky system are discretized with boundary elements while the domain of each block is divided into internal cells only for the integration of the domain integral of the initial stress term. The contact forces among blocks are treated as concentrated forces and the open–close iterations are applied to ensure the computational accuracy of block interactions. In the current method, an implicit time integration scheme is adopted for numerical stability. Three examples are used to show the effectiveness of the algorithm in simulating block movement, sliding, deformation and interaction of blocks. At last, block toppling and tunnel stability examples are conducted to demonstrate that the BE‐DDA is applicable for simulation of blocky systems. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
9.
岩体饱和–非饱和渗流、应力耦合作用对工程岩体的强度和稳定性有十分重要的影响。目前对于裂隙岩体饱和渗流应力耦合的研究取得了一些进展,但在很多工程领域不能简单地采用饱和渗流分析。根据DDA力学计算和非饱和渗流计算原理,提出了新的基于DDA方法的非饱和渗流应力耦合模型。并给出了在库水位骤降工况下的边坡水力耦合算例,其计算结果显示在库水位骤降情况下：考虑水力耦合且库水位下降较快时的安全系数要小于库水位下降慢时的安全系数;考虑水力耦合的边坡稳定安全系数要小于不考虑耦合时的安全系数。仿真实验和工程应用表明其计算成果是符合实践规律的。 相似文献
10.
数值流形方法是包含流形元、有限元及DDA在内的数值方法体系,建立流形元与DDA块体的接触方程,则可实现流形方法框架下的连续介质和散体系统共同作用模拟。针对填石路堤工程,编制了大型数值计算程序,采用块体随机生成、块体粒径控制及块体自然堆积的方法建立散体系统的DDA模型,对路堤的分层铺设、碾压及工后沉降变形等进行模拟分析。通过算例表明,在数值流形方法框架下,采用流形元与DDA共同作用的方法,可以很好地对同时存在连续变形和散体大变形的体系进行计算分析,其对该类问题的模拟更接近分析对象的实际情况,有助于从根本上揭示分析对象变形的细观机制和规律,并能考察更多因素对工程问题的影响。 相似文献
11.
Three porous media flow problems, in which the fluid mechanical interactions are critical, are studied in a mesoscopic–microscopic coupling system. In this system, fluid flow in the pore space is explicitly modeled at mesoscopic level by the lattice Boltzmann method, the geometrical representation and the mechanical behavior of the solid skeleton are modeled at microscopic level by the particulate distinct element method (DEM), and the interfacial interaction between the fluid and the solids is resolved by an immersed boundary scheme. In the first benchmark problem, the well‐known and frequently utilized Ergun equation is validated in periodic particle and periodic pore models. In the second problem, the upward seepage problem is simulated over three stages: The settlement of the column of sphere under gravity loading is measured to illustrate the accuracy of the DEM scheme; the system is solved to hydrostatic state with pore space filled with fluid, showing that the buoyancy effect is captured correctly in the mesoscopic–microscopic coupling system; then, the flow with constant rate is supplied at the bottom of the column; the swelling of the ground surface and pore pressure development from the numerical simulation are compared with the predictions of the macroscopic consolidation theory. In the third problem, the fluid‐flow‐induced collapse of a sand arch inside a perforation cavity is tested to illustrate a more practical application of the developed system. Through comparing simulation results with analytical solutions, empirical law and physical laboratory observations, it is demonstrated that the developed lattice Boltzmann–distinct element coupling system is a powerful fundamental research tool for investigating hydromechanical physics in porous media flow. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
12.
水库淤积形态是影响库容分布、水库排沙的一项重要因素。小浪底水库近坝段淤积泥沙粒径极细,具有流动性,针对其细颗粒泥沙淤积特点,揭示了水库细颗粒淤积物的流变特性与流型特征;通过引入水、淤积物、床面之间的界面受力分析,构建了细颗粒淤积物失稳流动描述模式,并与水沙输移模型相耦合,建立了考虑细颗粒淤积物流动特性的水库淤积形态模拟方法,在此基础上对小浪底水库淤积形态进行了验证分析。研究结果表明:低密度细颗粒淤积物为宾汉型流体,淤积平衡坡降较小,当其密度大于1.25 g/cm3后,流动性快速减弱;考虑细颗粒淤积物流动特征的水库淤积形态模拟结果与实测结果吻合较好。研究成果可为水库淤积形态形成机理及其对水沙调控的响应研究提供技术支撑。 相似文献
13.
14.
在明渠中开展不同水流条件下低强度均匀沙平衡输沙试验,基于灰度相减方法分析整个床面推移质运动特征,以探寻紊流结构与推移质运动之间的作用机理。结果表明:①受水槽横向方向水流强度分布特征影响,推移质运动概率从水槽中线到两侧壁逐渐变小,且基本呈对称分布;②在紊流相干结构作用下,推移质床面在紊流低速条带区形成凸槽,高速条带区形成凹槽,推移质运动概率沿水槽横向方向存在高低相间的带状分区;③随摩阻雷诺数增大,相邻两推移质运动高概率区域的间距值变化范围为0.13~0.24倍水槽宽度,其值随摩阻雷诺数增大而增大;④不同水流条件下,推移质运动高概率区域间距值约为水深的2倍,这与流向涡模型吻合,表明流向涡是诱导床面出现凹凸相间形态的重要因素。 相似文献
15.
水库淤积形态是影响库容分布、水库排沙的一项重要因素。小浪底水库近坝段淤积泥沙粒径极细,具有流动性,针对其细颗粒泥沙淤积特点,揭示了水库细颗粒淤积物的流变特性与流型特征;通过引入水、淤积物、床面之间的界面受力分析,构建了细颗粒淤积物失稳流动描述模式,并与水沙输移模型相耦合,建立了考虑细颗粒淤积物流动特性的水库淤积形态模拟方法,在此基础上对小浪底水库淤积形态进行了验证分析。研究结果表明:低密度细颗粒淤积物为宾汉型流体,淤积平衡坡降较小,当其密度大于1.25 g/cm 3后,流动性快速减弱;考虑细颗粒淤积物流动特征的水库淤积形态模拟结果与实测结果吻合较好。研究成果可为水库淤积形态形成机理及其对水沙调控的响应研究提供技术支撑。 相似文献
16.
In engineering practices, different numerical methods for fluid flow simulation and solid deformation/stress simulation are adopted to model fluid–structure interaction problems in porous media. Cell‐centered finite volume method is widely used in fluid flow simulation, while the solid deformation/stress simulation is usually accomplished by using the Galerkin vertex‐centered finite element method, which leads to the incompatibility between cell variables with nodal variables. Therefore, the data transfer between cell variables and nodal variables is inevitable. Consequently, this kind of transfer will lead to extra artificial error. Hence, the major concern is how to minimize the error due to cell to node projections. In this paper, a problem of pore pressure diffusion within a one‐dimensional heterogeneous porous medium is investigated. We present a new projection scheme and corresponding error formula, where the error control factor is introduced. The new projection scheme is based on piecewise linear interpolations. Results demonstrate that if the error control factor is chosen properly, the error due to the projection from cell to node can be controlled effectively, and the most desired zero error can be achieved. Finally, we analyze some practical cases in consideration of permeability contrast and mesh uniformity. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
17.
《International Journal of Mineral Processing》2005,76(3):163-172
The electrokinetics and flotation behavior of kaolinite have been investigated through flotation experiments, zeta potential measurements, and FTIR spectral analysis. These results have been explained based on crystal structure considerations and quantum mechanical calculations. It has been shown that the PZC of kaolinite is at pH 4.2. Kaolinite exhibits good floatability in weakly acidic solution with dodecyl amine (DDA) collector. The negative zeta potential of kaolinite increases with an increase in pH, and hence the adsorption of DDA on kaolinite increases. However, the flotation of kaolinite decreases with an increase in pH and an increase in DDA adsorption. This anomalous flotation behavior is rationalized based on crystal structure considerations and particle aggregation phenomena. It is considered that the silica (001) and the alumina (001¯) basal planes of kaolinite are negatively charged. From quantum chemistry calculations, the interaction between DDA and the (001) plane was found to be stronger than at the (001¯) alumina plane due to the difference in the structure at the (001) and (001¯) planes. The self-aggregation between (001¯) faces and the edge planes and the adsorption of DDA at the silica (001) plane make the kaolinite aggregates hydrophobic, and good floatability is achieved in acidic solution. In alkaline solution, the kaolinite particles are dispersed. In the presence of DDA, hydrophobic aggregation appears to occur in alkaline solution between the (001) planes due to adsorbed DDA, and thus the hydrophilic (001¯) faces are exposed and flotation is not achieved. 相似文献
18.
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. 相似文献
19.
Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method 下载免费PDF全文
下载免费PDF全文 A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
20.
M. Rouainia H. Lewis C. Pearce N. Bicanic G.D. Couples M.A. Reynolds 《Engineering Geology》2006,82(4):222-233
A coupled hydro-geomechanical modelling environment, developed to evaluate the coupled responses of fluid flow in deforming discontinuous media, is described. A staggered computational framework is presented, where the two simulations tools, HYDRO and DDA, communicate via the mapping of an equivalent porosity (and related permeabilities) from the rock system to the fluid phase and an inverse mapping of the pressure field. Several algorithmic and modelling issues are discussed, in particular the computational procedure to map the current geometry of the discontinuous rock blocks assembly into an equivalent porosity (and permeability) field. A generic, geometrically simple, overpressured reservoir/seal system is analysed for illustration. Further examples investigate discontinuous, fractured configurations in flexure causing a degree of spatial variability in the induced stresses. Model predictions show that the combination of hydraulic and mechanical loads causes a dilational opening of some pre-existing fractures and closure of others, with strong localisation of the modified flow pattern along wider fracture openings. 相似文献