共查询到20条相似文献，搜索用时 15 毫秒
1.
为了充分发挥有限元与颗粒离散元各自的优势,提出了一种由有限元转化为颗粒流的方法。数值模型首先用较粗的有限元网格进行离散,并在单元上引入连续介质本构模型。力学计算开始后,实时跟踪各单元的应力状态。一旦某单元的应力满足Mohr-Coulomb准则或最大拉应力准则,删除该单元,同时创建具有一定数目、随机分布且微嵌入的颗粒簇。其后,该单元所在区域的非连续变形及失稳断裂由颗粒簇演化获得。各颗粒的质量、材料参数、速度、位移、接触力等信息根据插值从有限元单元中继承。为了实现有限元与颗粒流接触面的耦合计算,引入了点-棱（二维）及点-面（三维）接触模型,通过法向及切向弹簧实现接触力的计算。颗粒球与有限元板的碰撞分析、单轴压缩、岩石切割等案例展示了上述方法的精确性及合理性。 相似文献
2.
圆柱体点源场有限单元模拟 总被引：3，自引：0，他引：3
我们对有限长度园柱形电性不均匀体的点源场的有限单元法数值模拟进行了研究和讨论。利用“镜像法”原理将三维数值模拟问题转化成二维问题,从理论上严格地推导出了求“旁测剖面”电位的公式。依据此算法编制出的程序对各种模型进行了计算,结果令人满意。 相似文献
3.
Numerical modelling of multistage caving processes: insights from multi‐seam longwall mining‐induced subsidence 
下载免费PDF全文
下载免费PDF全文 Evaluating the induced subsidence is a critical step in multi‐seam longwall mining. Numerical modelling can be a cost‐effective approach to this problem. Numerical evaluation of longwall mining‐induced subsidence is much more complicated when more than one seam is to be extracted. Only a few research works have dealt with this problem. This paper discusses the essential requirements of a robust numerical modelling approach to simulation of multi‐seam longwall mining‐induced subsidence. In light of these requirements, the previous works on this topic are critically reviewed. A simple yet robust FEM‐based modelling approach is also proposed that is capable of simulating caving process, rock mass deterioration and subsidence around multi‐seam excavations. The effectiveness of this approach in comparison with two other conventional FEM approaches is demonstrated through numerical examples of two different multi‐seam mining configurations. Results show that the proposed numerical modelling approach is the only robust method, which is capable of simulating multi‐seam subsidence in both demonstrated cases. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
4.
In this paper a typical soil–structure interaction problem is considered, the case of a vertical pile installed in sand and submitted to an axial compression loading. Results from two full scale pile tests are analysed and the tests are reproduced by numerical simulations via finite elements method (FEM). The choice of the mechanical parameters for the soil and the sand–pile interface and the modelling approach are first described. A new numerical strategy is outlined to account for pile installation effects due to jacking and driving via FEM. The proposed approach is based on the application of existing empirical correlations available for the quantification of residual radial and shear stresses along the pile shaft as well as residual pressures around the pile base after the installation. This approach is proposed as an alternative to more complex methods based on the numerical modelling of the pile penetration problem. The role of the constitutive modelling of the interface is also discussed. Finally, comparative analyses of pile loading tests using FEM are provided and the comparisons between numerical and experimental results are presented and discussed. 相似文献
5.
为了模拟岩体中裂纹的萌生、扩展,Munjiza提出了有限元法/离散元法(FEM/DEM)耦合分析方法。因为裂纹是沿单元边界进行扩展的,亦即裂纹扩展具有网格依赖性,为获得较好的裂纹扩展形态,需要划分密集的初始网格。为解决上述难题,基于FEM/DEM耦合分析方法,提出了基于局部单元动态劈裂的FEM/DEM自适应分析方法,以克服裂纹扩展形态对网格的依赖性。该方法在最初建模时无需划分很密的初始网格,随着荷载的施加,对裂纹尖端附近的局部单元进行动态劈裂,为裂纹的后续扩展提供了更多可能的扩展方向,使得裂纹扩展不必沿着初始网格的单元边界扩展,即可以沿着单元内部进行扩展,裂纹扩展形态更为平滑,与实际情况更为接近。同时相对原FEM/DEM耦合分析方法一开始就划分很密的网格而言,新方法可以划分较为稀疏的初始网格,计算成本降低。最后,通过巴西劈裂算例与原FEM/DEM耦合分析方法对比,分析表明,新方法在一定程度上克服了裂纹扩展形态对初始网格的依赖性。 相似文献
6.
Strain gradient implies an important characteristic in localized damage deformation, which can be observed in the softening state of brittle materials, and strain gradients constitute the basic behaviours of localization failure area of the materials. The most important point in strain gradient is its damaging function including an internal length scale, which can be used to express the scale effects of mechanical responses of brittle rock mass. By extending the strain gradient theory and introducing an intrinsic material length scale into the constitutive law, the authors develop an isotropic damage model as well as a micro‐crack‐based anisotropic damage model for rock‐like materials in this paper. The proposed models were used to simulate the damage localization under uniaxial tension and plain strain compression, respectively. The simulated results well illustrated the potential of these models in dealing with the well‐known mesh‐sensitivity problem in FEM. In the computation, elements with C1 continuity have been implemented to incorporate the proposed models for failure localization. When regular rectangle elements are encountered, the coupling between finite difference method (FDM) and conventional finite element method (FEM) is used to avoid large modification to the existing FEM code, and to obtain relatively higher efficiency and reasonably good accuracy. Application of the anisotropic model to the 3D‐non‐linear FEM analysis of Ertan arch dam has been conducted and the results of its numerical simulation coincide well with those from the failure behaviours obtained by Ertan geophysical model test. In this paper, new applications of gradient theories and models for a feasible approach to simulate localized damage in brittle materials are presented. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
7.
The lattice Boltzmann (LB) method is an efficient technique for simulating fluid flow through individual pores of complex porous media. The ease with which the LB method handles complex boundary conditions, combined with the algorithm’s inherent parallelism, makes it an elegant approach to solving flow problems at the sub-continuum scale. However, the realities of current computational resources can limit the size and resolution of these simulations. A major research focus is developing methodologies for upscaling microscale techniques for use in macroscale problems of engineering interest. In this paper, we propose a hybrid, multiscale framework for simulating diffusion through porous media. We use the finite element (FE) method to solve the continuum boundary-value problem at the macroscale. Each finite element is treated as a sub-cell and assigned permeabilities calculated from subcontinuum simulations using the LB method. This framework allows us to efficiently find a macroscale solution while still maintaining information about microscale heterogeneities. As input to these simulations, we use synchrotron-computed 3D microtomographic images of a sandstone, with sample resolution of 3.34 μm. We discuss the predictive ability of these simulations, as well as implementation issues. We also quantify the lower limit of the continuum (Darcy) scale, as well as identify the optimal representative elementary volume for the hybrid LB–FE simulations. 相似文献
8.
Many problems in mining and civil engineering require using numerical stress analysis methods to repeatedly solve large models. Widespread acceptance of tunneling methods, such as New Austrian Tunneling Method, which depend heavily on numerical stress analysis tools and the fact that the effects of excavation at the face of a tunnel are distinctively three–dimensional (3D), necessitates the use of 3D numerical analysis for these problems. Stress analysis of a practical mining problem can be very lengthy, and the processing time can be measured in days or weeks at times. A framework is developed to facilitate efficient modeling of underground excavations and to create an optimal 3D mesh by reducing the number of surface and volume elements while keeping the result of stress analysis accurate enough at the region of interest, where a solution is sought. Fewer surface and volume elements mean fewer degrees of freedom in the numerical model, which directly translates into savings in computational time and resources. The mesh refinement algorithm is driven by a set of criteria that are functions of distance and visibility of points from the region of interest, and the framework can be easily extended by adding new types of criteria. This paper defines the framework, whereas a second companion paper will investigate its efficiency, accuracy and application to a number of practical mining problems. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
9.
This paper presents a computational method able to effectively model both the simultaneous processes typically observed in backward erosion piping, ie, the pipe tip propagation and the conduit cross-section enlargement. The numerical method is based on the novel formulation of a problem of localized erosion along a line propagating in a multidimensional porous medium. In this line, a conduit with evolving transverse size is embedded, which conveys a multiphase flow. The two systems, porous medium and pipe, are bridged by exchange terms of multiphase fluid mass and by a shared fluid pressure field. On the contrary, different fields are considered to describe flows, which are assumed as Darcian in the porous medium and turbulent in the conduit. These two flows drive pipe propagation and enlargement, respectively, as modeled by means of proper erosion kinetic laws. The corresponding numerical formulation is based on the combination between one- and multidimensional finite elements, to model the erosion conduit and the porous medium, respectively. Several simulations are proposed to demonstrate the ability of the proposed approach in reproducing available experimental data of real-scale tests on levees. Our results point out the crucial role played by the combined influence of pipe propagation and enlargement, as well as of three-dimensional (3D) effects. We also assess the mesh independence of the proposed numerical solution, particularly as concerns the calculated pipe propagation history. 相似文献
10.
井孔与井列的渗流数值计算 总被引：3，自引：2，他引：3
对有限元和有限差分的数值计算提出了井孔奇眯修正和以沟代井列的计算方法，从而提高了离散化网络的计算精度，并可简化取代三维问题为二维。为推广电阻网模拟井孔等问题的优点，已综合编制了网络程序NETW，拓宽数值法的领域。 相似文献
11.
For non-linear dynamic problems, it has been recognized that an explicit time-integration method of approach is a very efficient way of solving the dynamic equations of motion. The numerical formulation and computation for such problems fall into the two general categories of finite elements and finite differences. Over the years, there have been many arguments between schools which adopt the finite element approach and those which adopt the finite difference approach. At one extreme, arguments areconcerned with the superiority of each approach and at the other end of the spectrum the arguments are about which approach is a subset of the other. The most common of these arguments are concerned with efficiency and accuracy. This publication addresses the accuracy issue with specific reference to explicit calculations in which the analysis domain is discretized into triangular or quadrilateral plane-strain elements. It concludes that if the same basic assumptions are made in the two approaches, they, will give identical answers for problems in this category. 相似文献
12.
The numerical simulation of soil-pile interaction problems, by means of full 3D finite element models, involves a large number of degrees of freedom (DOF) and difficulties during the mesh generation process. In order to reduce the unknowns and simplify and properly analyze such class of geotechnical problems, the so-called embedded beam elements (EBE) have recently been developed. In a preceding contribution of the authors, an improved EBE formulation, which brings into play the soil-pile interaction surface, was proposed with the aim to localize material plasticity in the soil surrounding the pile. This embedded beam model couples two different finite elements, each described by distinct kinematics (ie, solid and beam). The coupling is incorporated in the formulation by means of kinematical constrains established over the solid and beam displacement fields on the interaction surface. One of the main advantages of the embedded elements is that the addition of beams structural members immersed within the 3D soil model does not represent a constraint for the solid mesh, which can be adopted independently from the beam mesh. In this paper, the lateral loading of pile groups is studied by means of the proposed EBE approach with elasto-plastic interfaces. In order to represent a rigid cap, a master node and a special set of kinematical restrictions are incorporated into the formulation. The paper presents results obtained by means of the present formulation compared against other well-established analysis methods and test results published in the literature, for both elastic and elasto-plastic cases. 相似文献
13.
This paper presents a second-order work analysis in application to geotechnical problems by using a novel effective multiscale approach. To abandon complicated equations involved in conventional phenomenological models, this multiscale approach employs a micromechanically-based formulation, in which only four parameters are involved. The multiscale approach makes it possible a coupling of the finite element method (FEM) and the micromechanically-based model. The FEM is used to solve the boundary value problem (BVP) while the micromechanically-based model is utilized at the Gauss point of the FEM. Then, the multiscale approach is used to simulate a three-dimensional triaxial test and a plain-strain footing. On the basis of the simulations, material instabilities are analyzed at both mesoscale and global scale. The second-order work criterion is then used to analyze the numerical results. It opens a road to interpret and understand the micromechanisms hiding behind the occurrence of failure in geotechnical issues. 相似文献
14.
A simple three-dimensional heat transfer model is developed to consider the hindering effect of cracks on heat transfer. The 3D heat transfer model can also be applied to numerical methods such as the combined finite-discrete element method (FDEM), discrete element method (DEM), discontinuous deformation analysis (DDA), the numerical manifold method (NMM), and the finite element method (FEM) to construct thermo-mechanical coupling models that allow these methods to solve thermal cracking problems and dynamically consider the hindering effect of cracks on heat transfer. In the 3D heat transfer model, the continuous-discontinuous medium is discretized into independent tetrahedral elements, and joint elements are inserted between adjacent tetrahedral elements. Heat transfer calculations for continuous-discontinuous media are converted to heat conduction in tetrahedral elements and the heat exchange between the adjacent tetrahedral elements through the joint element. If the joint element between adjacent tetrahedral elements breaks (ie, a crack generates), the heat exchange coefficient of the joint element is reduced to account for the hindering effect of cracks on heat conduction. Then the model and the FDEM are combined to build a thermo-mechanical coupling model to simulate thermal cracking. The thermally induced deformation, stress, and cracking are investigated by the thermo-mechanical coupling model, and the numerical results are compared with analytical solutions or experimental results. The 3D heat transfer model and thermo-mechanical model can provide a powerful tool for simulating heat transfer and thermal cracking in a continuous-discontinuous medium. 相似文献
15.
用三维梁-颗粒模型BPM3D(beam-particlemodelinthreedimensions)对岩石类非均质脆性材料的力学性质和破坏过程进行了数值模拟。梁-颗粒模型是在离散单元法基础上,结合有限单元法中的网格模型提出的用于模拟岩石类材料损伤破坏过程的数值模型。在模型中,材料在细观层次上被离散为颗粒单元集合体,相邻颗粒单元由有限单元法中的弹脆性梁单元联结。梁单元的力学性质均按韦伯(Weibull)分布随机赋值,以模拟岩石类材料力学参数的空间变异性。材料内部裂纹通过断开梁单元来模拟。通过自动生成的非均质材料模型对岩石类材料的破坏机理进行研究。岩石类非均质脆性材料在单轴压缩状态下破坏过程细观数值模拟结果显示,岩石材料宏观破坏是由于其内部细观裂纹产生、扩展、贯通的结果。通过数值模拟结果之间的对比分析,揭示出岩石试样宏观破坏模式随细观层次上韦伯分布参数的变化而不同。与实际矿柱破坏形态的对比分析表明了模型的适用性。根据数值模拟结果对岩石类非均质材料的破坏机理进行了探讨。 相似文献
16.
Numerical limit analysis of three-dimensional slope stability problems in catchment areas 总被引：1，自引：1，他引：0
Risk analysis of existing slopes in catchment areas requires quantification of their stability. This quantification becomes particularly difficult when dealing with larger areas under 3D conditions and including saturated and unsaturated water flow. This paper proposes the use of an effective numerical procedure to solve three-dimensional slope stability problems in large areas subjected to pore pressure effects. This numerical approach, numerical limit analysis, utilizes the finite element method and mathematical programming techniques. Mathematical programming is needed because the basic plasticity theorems for limit analysis can be cast as optimization problems. The generated optimization problem is formulated under a second-order cone programming framework, which is known to solve large-scale problems with great computational efficiency. The main objective of this work was to determine the slope safety factor and the collapse mechanism of soils governed by the Drucker–Prager yield criterion for large-scale 3D problems including pore pressure effects. This approach is applied to an experimental catchment in the Oregon Coast Range that failed after an intense rainfall. The results were compared with a previous stability analysis of the area available in the literature that used a novel 3D limit equilibrium method. 相似文献
17.
This paper concerns the buckling analysis of three-dimensional (3D) layered continua using the Cosserat approach. The finite element method (FEM) is used to implement the adopted Cosserat formulation. As a result, the interfaces between the layers need not be explicitly modelled. Instead, the internal characteristic length, i.e., the layer thickness and also the interaction conditions between the layers are incorporated into the governing equations of the solution. This paper introduces a new 3D geometric stiffness matrix based on the principle of virtual work. The proposed geometric stiffness matrix is applied to the linear buckling analysis of a number of benchmark problems with various geometries, boundary conditions, and interaction conditions between the layers. In all cases, the FEM Cosserat solution exhibits a high level of consistency with the analytical solution. 相似文献
18.
This paper presents a novel approach to the simulation of NATM tunnel construction using the Boundary Element Method (BEM) as principal numerical method. This new approach has the advantage that only the excavation surface, the possible plastic zones and the tunnel lining have to be discretised. The whole rock mass is represented by the BEM whereas the Finite Element Method (FEM) is used to represent the tunnel lining only. Thus, a general coupling strategy for coupling three-dimensional boundary elements with shell finite elements (shotcrete) and beam finite elements (steel arches) is presented. To achieve realistic results the effect of hydration of the shotcrete and yielding of the steel arches is considered in the excavation process. Furthermore, the nonlinear rock behaviour is modelled more realistically by using a powerful hierarchical constitutive model which considers a large range of rock materials. The combination of these ideas results in higher user-friendliness and efficiency. Some verification tests and practical applications in tunnelling are presented. 相似文献
19.
Explicit dynamic relaxation is an efficient tool that has been used to solve problems involving highly non-linear differential equations. The key feature of this method is the ability to use explicit dynamic algorithms in solving static problems. Few attempts have been made to date to apply this technique in conventional geotechnical engineering. In this study, an algorithm that incorporates the application of a stiffness dependent time step scheme is proposed. The algorithm has been successfully used to solve 2D and 3D non-linear geotechnical engineering problems. To calibrate the developed algorithm, numerical simulations have been conducted for a strip and square footings supported by Mohr–Coulomb material. Performance of four different types of brick elements used in collapse load calculation is examined in terms of convergence speed and accuracy. In addition, the role of employing adaptive time steps in reducing the number of iterations needed for convergence is also evaluated. 相似文献
20.
Nathan Benkemoun Hayder Al Khazraji Philippe Poullain Marta Choinska Abdelhafid Khelidj 《国际地质力学数值与分析法杂志》2018,42(3):449-468
In this paper, mesoscale hydromechanical simulations are performed to study (1) fracture features and (2) crack‐gas permeability coupling evolution in the context of the tensile splitting test. The mesostructure is based on a 2‐phase 3‐D representation of heterogeneous materials, such as concrete, where stiff aggregates are embedded into a mortar matrix. To take into account these heterogeneities without any mesh adaptation, a weak discontinuity is introduced into the strain field. In addition, a strong discontinuity is also added to take into account microcracking. This mechanical model is cast into the framework of the enhanced finite element method. Concerning the coupling with gas permeability, a double‐porosity method is used to simulate the flow through the cracks and the porosity. The apparent gas permeability is afterwards evaluated by a homogenization method. On the basis of finite element simulations, influence of aggregate size on ultimate crack opening, macroscopic ultimate tensile stress, total dissipated energy, and gas permeability evolution is numerically investigated. Furthermore, gas permeability evolution is also compared with experimental results from the literature. In addition, in the spirit of a sequential multiscale approach, macroscale gas permeability equations are identified from the hydromechanical results coming from the mesoscale computations. These equations lead to a relation between macroscale gas permeability evolution and crack opening. Besides, we show how the aggregate size influences the percolation threshold and that after this threshold, a cubic relation between macroscale gas permeability and crack opening is obtained. 相似文献