共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
Computational fluid dynamics and discrete element method (CFD–DEM) is extended with the volume of fluid (VOF) method to model free‐surface flows. The fluid is described on coarse CFD grids by solving locally averaged Navier–Stokes equations, and particles are modelled individually in DEM. Fluid–particle interactions are achieved by exchanging information between DEM and CFD. An advection equation is applied to solve the phase fraction of liquid, in the spirit of VOF, to capture the dynamics of free fluid surface. It also allows inter‐phase volume replacements between the fluid and solid particles. Further, as the size ratio (SR) of fluid cell to particle diameter is limited (i.e. no less than 4) in coarse‐grid CFD–DEM, a porous sphere method is adopted to permit a wider range of particle size without sacrificing the resolution of fluid grids. It makes use of more fluid cells to calculate local porosities. The developed solver (cfdemSolverVOF) is validated in different cases. A dam break case validates the CFD‐component and VOF‐component. Particle sedimentation tests validate the CFD–DEM interaction at various Reynolds numbers. Water‐level rising tests validate the volume exchange among phases. The porous sphere model is validated in both static and dynamic situations. Sensitivity analyses show that the SR can be reduced to 1 using the porous sphere approach, with the accuracy of analyses maintained. This allows more details of the fluid phase to be revealed in the analyses and enhances the applicability of the proposed model to geotechnical problems, where a highly dynamic fluid velocity and a wide range of particle sizes are encountered. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
3.
The flow of water–kaolinite mixtures exhibits a non‐Newtonian nature that differs from the flow of Newtonian fluid. The varying viscosities and shear history of non‐Newtonian fluid flows necessitate the use of a rheology model in moving particle semi‐implicit (MPS) for the numerical studies. On the other hand, the Lagrangian method has the advantage of handling free surface flows with large deformation and fragmentation. This study proposes a mesh‐free Lagrangian method, namely, the MPS method, together with a simple rheology model to investigate the non‐Newtonian free surface flows. The rheological parameters required in the rheology model are determined based upon experiments. The proposed method is applied to a water–kaolinite mixture collapse problem and is proved to be capable of reproducing the significant flow features observed in laboratory experiments. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
4.
A fully coupled transient two‐dimensional model was employed to study fundamentals of flood‐induced surface erosion in a particle bed. The interaction of the liquid and solid phases is the key mechanism related to surface erosion. The solid phase was idealized at a particle scale by using the discrete element method. The fluid phase was modeled at a mesoscale level and solved using the lattice Boltzmann method. The fluid forces applied on the particles were calculated on the basis of the momentum the fluid exchanges with the particle. The proposed approach was used to model both single particles and particle beds subjected to Couette flow conditions. The behavior of both the single particle and the particle bed depended on particle diameter and surface shear fluid velocity. The conducted simulations show that the fluid flow profile penetrates the bed for a small distance. This penetration initiates sheet‐flow and surface erosion as the fluid interacts with particles. The effect of suppressing particle rotation on the fluid‐induced forces on the particle was also examined. Suppressing particle spinning may lead to underestimated erosion rate. Results of fluid and particle velocities were compared against experimental results and appeared to agree with the observed trends.Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
5.
A three‐phase soil model is proposed to simulate stress wave propagation in soil mass to blast loading. The soil is modelled as a three‐phase mass that includes the solid particles, water and air. It is considered as a structure that the solid particles form a skeleton and their voids are filled with water and air. The equation of state (EOS) of the soil is derived. The elastic–plastic theory is adopted to model the constitutive relation of the soil skeleton. The damage of the soil skeleton is also modelled. The Drucker–Prager strength model including the strain rate effect is used to describe the strength of the soil skeleton. The model is implemented into a hydrocode Autodyn. The recorded results obtained by explosion tests in soil are used to validate the proposed model. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
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.
A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
8.
Time‐domain analysis of dynamic soil–structure interaction based on the substructure method plays an increasing role in practical applications as compared with the frequency‐domain analysis. Efficient and accurate modelling of the unbounded soil or rock medium has been a key issue in such an analysis. This paper presents a subregional stepwise damping‐solvent extraction formulation for solving large‐scale dynamic soil–structure problems in the time domain. Accuracy and efficiency of the formulation are evaluated in detail for a classical problem involving a rigid strip foundation embedded in a half‐space. A practical large‐scale soil–structure interaction problem, which represents a high concrete gravity dam subjected to seismic load, is then analysed using the proposed method. Various responses of the dam, including time histories of the crest displacement and acceleration and contours of the peak principal stresses within the dam body, are presented. Comparisons are also made between these results with those obtained using other models for the unbounded medium. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
9.
A modulus‐multiplier approach, which applies a reduction factor to the modulus of single pile p–y curves to account for the group effect, is presented for analysing the response of each individual pile in a laterally loaded pile group with any geometric arrangement based on non‐linear pile–soil–pile interaction. The pile–soil–pile interaction is conducted using a 3D non‐linear finite element approach. The interaction effect between piles under various loading directions is investigated in this paper. Group effects can be neglected at a pile spacing of 9 times the pile diameter for piles along the direction of the lateral load and at a pile spacing of 6 times the pile diameter for piles normal to the direction of loading. The modulus multipliers for a pair of piles are developed as a function of pile spacing for departure angle of 0, 90, and 180sup>/sup> with respect to the loading direction. The procedure proposed for computing the response of any individual pile within a pile group is verified using two well‐documented full‐scale pile load tests. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
10.
This paper presents a non‐linear soil–structure interaction (SSI) macro‐element for shallow foundation on cohesive soil. The element describes the behaviour in the near field of the foundation under cyclic loading, reproducing the material non‐linearities of the soil under the foundation (yielding) as well as the geometrical non‐linearities (uplift) at the soil–structure interface. The overall behaviour in the soil and at the interface is reduced to its action on the foundation. The macro‐element consists of a non‐linear joint element, expressed in generalised variables, i.e. in forces applied to the foundation and in the corresponding displacements. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. Mechanisms of yielding and uplift are modelled through a global, coupled plasticity–uplift model. The cyclic model is dedicated to modelling the dynamic response of structures subjected to seismic action. Thus, it is especially suited to combined loading developed during this kind of motion. Comparisons of cyclic results obtained from the macro‐element and from a FE modelization are shown in order to demonstrate the relevance of the proposed model and its predictive ability. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
11.
This paper presents a unified modeling framework to investigate the impacts of debris flow on flexible barriers, based on coupled computational fluid dynamics and discrete element method (CFD‐DEM). We consider a debris flow as a mixture of fluid and particles where the fluid and particle phases are modeled by the CFD and the DEM, respectively. The fluid‐particle coupling is considered by the exchange of interaction forces between CFD and DEM calculations. The flexible barrier is simulated by the DEM as a network of bonded particles with remote interactions. The proposed coupled CFD‐DEM approach enables us to conveniently handle the complicated three‐way interactions among the fluid, the particles, and the flexible barrier structure for debris flow impact simulations. The proposed approach is first used to investigate the influences of channel inclination and the volumetric solid fraction in a debris mixture on the impact force, the resultant deformation, and the retained mass in a flexible barrier. The predictions agree well with existing experimental and numerical studies. We further examine the possible failure modes of a flexible barrier under debris flow impact and their underlying mechanisms. The performance of different components in a flexible barrier system, including single wires, double twists and cables, and their load sharing mechanisms, are carefully evaluated. The proposed unified framework offers a novel, promising pathway towards physically based, quantitative analysis and design of flexible barriers for debris flow mitigation. 相似文献
12.
A novel three‐dimensional particle‐based technique utilizing the discrete element method is proposed to analyze the seismic response of soil‐foundation‐structure systems. The proposed approach is employed to investigate the response of a single‐degree‐of‐freedom structure on a square spread footing founded on a dry granular deposit. The soil is idealized as a collection of spherical particles using discrete element method. The spread footing is modeled as a rigid block composed of clumped particles, and its motion is described by the resultant forces and moments acting upon it. The structure is modeled as a column made of particles that are either clumped to idealize a rigid structure or bonded to simulate a flexible structure of prescribed stiffness. Analysis is done in a fully coupled scheme in time domain while taking into account the effects of soil nonlinear behavior, the possible separation between foundation base and soil caused by rocking, the possible sliding of the footing, and the dynamic soil‐foundation interaction as well as the dynamic characteristics of the superstructure. High fidelity computational simulations comprising about half a million particles were conducted to examine the ability of the proposed technique to model the response of soil‐foundation‐structure systems. The computational approach is able to capture essential dynamic response patterns. The cyclic moment–rotation relationships at the base center point of the footing showed degradation of rotational stiffness by increasing the level of strain. Permanent deformations under the foundation continued to accumulate with the increase in number of loading cycles. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
13.
Coupled bonded particle and lattice Boltzmann method for modelling fluid–solid interaction 下载免费PDF全文
下载免费PDF全文 This paper presents a two‐dimensional coupled bonded particle and lattice Boltzmann method (BPLBM) developed to simulate the fluid–solid interactions in geomechanics. In this new technique, the bonded particle model is employed to describe the inter‐particle movement and forces, and the bond between a pair of contacting particles is assumed to be broken when the tensile force or tangential force reaches a certain critical value. As a result the fracture process can be delineated based on the present model for the solid phase comprising particles, such as rocks and cohesive soils. In the meantime, the fluid phase is modelled by using the LBM, and the immersed moving boundary scheme is utilized to characterize the fluid–solid interactions. Based on the novel technique case studies have been conducted, which show that the coupled BPLBM enjoys substantially improved accuracy and enlarged range of applicability in characterizing the mechanics responses of the fluid–solid systems. Indeed such a new technique is promising for a wide range of application in soil erosion in Geotechnical Engineering, sand production phenomenon in Petroleum Engineering, fracture flow in Mining Engineering and fracture process in a variety of engineering disciplines. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
14.
Mohammad Reza Maleki Javan Asadollah Noorzad Manouchehr Latifi Namin 《国际地质力学数值与分析法杂志》2008,32(6):681-699
The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
15.
Artificial ground freezing (AGF) is a commonly used technique in geotechnical engineering for ground improvement such as ground water control and temporary excavation support during tunnel construction in soft soils. The main potential problem connected with this technique is that it may produce heave and settlement at the ground surface, which may cause damage to the surface infrastructure. Additionally, the freezing process and the energy needed to obtain a stable frozen ground may be significantly influenced by seepage flow. Evidently, safe design and execution of AGF require a reliable prediction of the coupled thermo‐hydro‐mechanical behavior of freezing soils. With the theory of poromechanics, a three‐phase finite element soil model is proposed, considering solid particles, liquid water, and crystal ice as separate phases and mixture temperature, liquid pressure, and solid displacement as the primary field variables. In addition to the volume expansion of water transforming into ice, the contribution of the micro‐cryo‐suction mechanism to the frost heave phenomenon is described in the model using the theory of premelting dynamics. Through fundamental physical laws and corresponding state relations, the model captures various couplings among the phase transition, the liquid transport within the pore space, and the accompanying mechanical deformation. The verification and validation of the model are accomplished by means of selected analyses. An application example is related to AGF during tunnel excavation, investigating the influence of seepage flow on the freezing process and the time required to establish a closed supporting frozen arch. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
16.
In this article, we present a numerical method to deal with fluid–solid interactions and simulate particle–fluid systems as encountered in soils. This method is based on a coupling between two methods, now widely used in mechanics of granular media and fluid dynamics respectively: the discrete element (DE) method and the lattice Boltzmann (LB) method. The DE method is employed to model interactions between particles, whereas the LB method is used to describe an interstitial Newtonian fluid flow. The coupling presented here is a full one in the sense that particle motions act on fluid flow and reciprocally. This article presents in details each of the two methods and the principle of the coupling scheme. Determination of hydrodynamic forces and torques is also detailed, and the treatment of boundaries is explained. The coupled method is finally illustrated on a simple example of piping erosion, which puts in evidence that the combined LB–DE scheme constitutes a promising tool to study coupled problems in geomechanics. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
17.
Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
18.
Coupled hydro‐mechanical processes in granular media represent the interaction between solid particles movement and fluid flow during external and/or internal loading. The processes attract attention in geotechnical engineering since they cause many disasters such as landslide, slope collapse, boiling/quicksand and soil liquefaction, etc. This study presents a numerical method combining discontinuous deformation analysis (DDA) for mechanical calculation and finite element method for fluid flow simulation to model the interaction between solid particles' movement and fluid flow from microscopic point of view. The term of hydraulic pressure was newly formulated and introduced to the original DDA. The simplified examples are given to verify the new method, and the computational results correlate well with the theoretical calculations, further development is also considered. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
19.
20.
Analysis of large deformation of geomaterials subjected to time‐varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator‐split arbitrary Lagrangian–Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid–fluid coupling and strong material non‐linearity. Each time step of the operator‐split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one‐dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献