共查询到20条相似文献，搜索用时 140 毫秒
1.
The Finite Element Method (FEM) is widely used in the simulation of geotechnical applications. Owing to the limitations of FEM to model problems involving large deformations, many efforts have been made to develop methods free of mesh entanglement. One of these methods is the Material Point Method (MPM) which models the material as Lagrangian particles capable of moving through a background computational mesh in Eulerian manner. Although MPM represents the continuum by material points, solution is performed on the computational mesh. Thus, imposing boundary conditions is not aligned with the material representation. In this paper, a non-zero kinematic condition is introduced where an additional set of particles is incorporated to track the moving boundary. This approach is then applied to simulate the seismic motion resulting in failure of slopes. To validate this simulation procedure, two geotechnical applications are modelled using MPM. The first is to reproduce a shaking table experiment where the results of another numerical method are available. After validating the present numerical scheme for relatively large deformation problem, it is applied to simulate progression of a large-scale landslide during the Chi-Chi earthquake of Taiwan in which excessive material deformation and transportation is taking place. 相似文献
2.
Jian-Hong Wu 《Computers and Geotechnics》2010,37(1-2):153-163
This work calculates the integrations of convex and concave polyhedrons in three-dimensional discontinuous deformation analysis (3D DDA) using a novel method. With this novel method, mathematical topology is applied to the integration, and the numerical mesh used by the Finite Element Method (FEM) is not needed to subdivide polyhedrons into tetrahedrons or hexahedrons. Examples demonstrate the accuracy of the proposed method for integrating convex and concave polyhedrons with planar boundaries. 相似文献
3.
G. Beer 《国际地质力学数值与分析法杂志》1983,7(1):1-7
The Finite Element Method is frequently used to analyse problems involving an ‘infinite domain’. A typical problem is an underground excavation in prestressed infinite medium in either tunnelling or mining operations or a foundation in an infinite half space. This paper examines the implications of mesh truncation on the accuracy of the accuracy of the solution. At the same time, a more economical and accurate method of analysing these problems using special ‘infinite domain’ Finite Elements is presented. 相似文献
4.
Coupled formulation and algorithms for the simulation of non‐planar three‐dimensional hydraulic fractures using the generalized finite element method 
下载免费PDF全文
下载免费PDF全文 This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
5.
A. M. Oskoorouchi B. Novrouzian
G. De Roeck
J. Van Den Broeck 《Computers and Geotechnics》1991,11(4):265-294The technique of Finite Strip method has been extended for zoned strips. One strip can have several different materials in horizontal direction including dummy materials. The new developed zoned Finite Strip (Z.F.S.) method has many engineering applications, particularly in geomechanical problems. The formulation is presented and three different geotechnical problems have been solved using the new Z.F.S. method, which could not be solved by the original Finite Strip method. The results are compared with the solution for the same problems obtained by analytical and/or numerical (Finite Element) methods. The reduced input and output for the Z.F.S. method compared to the Finite Element method makes this method attractive particularly for practising engineers. However the results can be as good as the well known Finite Element method. 相似文献
6.
This paper is devoted to the numerical reliability and time requirements of the Mixed Finite Element (MFE) and Mixed-Hybrid Finite Element (MHFE) methods. The behavior of these methods is investigated under the influence of two factors: the mesh discretization and the medium heterogeneity. We show that, unlike the MFE, the MHFE suffers with the presence of badly shaped discretized elements. Thereat, a numerical reliability analyzing software (Aquarels) is used to detect the instability of a matrix-inversion code generated automatically by a symbolic manipulator. We also show that the spectral condition number of the algebraic systems furnished by both methods in heterogeneous media grows up linearly according to the smoothness of the hydraulic conductivity. Furthermore, it is found that the MHFE could accumulate numerical errors if large jumps in the tensor of conductivity take place. Finally, we compare running-times for both algorithms by giving various numerical experiments. 相似文献
7.
Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. 相似文献
8.
The numerical analysis of shallow tunnels using a combined Finite Element-Boundary Element method is discussed. The theory of coupling the Boundary and Finite Element discretisation is presented and solutions obtained by the method compared with analytical results for near surface circular excavations. The coupled method is then applied to a practical problem in tunneling. 相似文献
9.
Seepage problems with complex drainage systems are commonly encountered in civil engineering, with strong non-linearity. A numerical solution based on the Finite Element Method combining the substructure technique with a variational inequality formulation of Signorini’s type is proposed to solve these problems. The aims of this work are to accurately characterize the boundary conditions of the drainage systems, to reduce the difficulty in mesh generation resulting from the drainage holes with small radius and dense spacing, and to eliminate the singularity at the seepage points and the resultant mesh dependency. Numerical stability and robustness of the proposed method are guaranteed by an adaptive procedure for progressively relaxing the penalized Heaviside function associated with the formulation of the discrete variational inequality. Two challenging numerical examples are presented to validate the effectiveness and robustness of the proposed method. 相似文献
10.
有限元的计算结果常需要用等值线表示出来。Surfer是一应用得十分广泛地作等值线的软件,但它仍有不足,如有的等值线所代表的值无法判断等。基于有限元网格和线性插值的方法及其他一些处理技巧,可有效地克服Surfer等软件的缺点,通过VBA语言编程可直接将等值线画在常用的制图软件AutoCAD界面里。所编的程序在平面有限元计算结果的等值线绘制中运用效果良好。 相似文献
11.
A Lagrangian numerical approach for the simulation of rapid landslide runouts is presented and discussed. The simulation approach is based on the so‐called Particle Finite Element Method. The moving soil mass is assumed to obey a rigid‐viscoplastic, non‐dilatant Drucker–Prager constitutive law, which is cast in the form of a regularized, pressure‐sensitive Bingham model. Unlike in classical formulations of computational fluid mechanics, where no‐slip boundary conditions are assumed, basal slip boundary conditions are introduced to account for the specific nature of the landslide‐basal surface interface. The basal slip conditions are formulated in the form of modified Navier boundary conditions, with a pressure‐sensitive threshold. A special mixed Eulerian–Lagrangian formulation is used for the elements on the basal interface to accommodate the new slip conditions into the Particle Finite Element Method framework. To avoid inconsistencies in the presence of complex shapes of the basal surface, the no‐flux condition through the basal surface is relaxed using a penalty approach. The proposed model is validated by simulating both laboratory tests and a real large‐scale problem, and the critical role of the basal slip is elucidated. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
12.
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. 相似文献
13.
采用迦辽金有限单元法对龙口市区地下水资源进行评价。按当地设计的开采量预报 ,本区未来地下水位降落漏斗范围有所缩小 ,地下水位负值区面积基本没有变化。只有采取如下措施 :增开新的水源地 ,充分利用地表水及加强水库工程 ,合理布局开采井 ,再生水利用及引进容水 ,方能达到最优的环境和社会效益。 相似文献
14.
给出分析各向异性非均质稳定随机渗流场问题的三维有限元模型;结合实际工程问题,统计分析长江荆南干堤士性参数的分布特征,通过Kolomogorov-Smirnov统计检验表明,渗透系数呈高斯分布假设可以接受;通过对长江荆南干堤随机渗流场的三维有限元统计模拟的数值分析,研究长江荆南干堤渗流场的各种随机特性,并进一步对随机模拟结果进行统计检验,验证模拟结果的合理性;在实际的分析研究中把上下游水位的随机波动引入三维有限元的随机分析模型,分析上下游水位的变异性对渗流场矢量的随机干扰和边界条件的随机性对随机渗流场分析结果变异性的影响。在此基础上进一步考虑施加诸如垂直截渗墙、下游导渗沟等抗渗措施后,它们作为复杂边界条件的扰动,在与场内土性参数的变异性共同影响下,对渗流场水头势分布的随机干扰特性,并与相应的确定性稳定渗流场问题的结果对比,证实随机渗流场研究的必要性、可行性及实用性。实现了对长江荆南干堤的三维渗流场的较为全面的随机场模拟及特性分析,分析得到的结论通过统计检验并结合实测工程数据对照证明是可靠的,所研制的程序是适用的。 相似文献
15.
Computational Geosciences - In this work, we consider an online enrichment procedure in the context of the Generalized Multiscale Finite Element Method (GMsFEM) for the two-phase flow model in... 相似文献
16.
A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
17.
G. Bernard-Michel C. Le Potier A. Beccantini S. Gounand M. Chraibi 《Computational Geosciences》2004,8(2):187-201
We present here results for the Andra Couplex 1 test case, obtained with the code Cast3m. This code is developped at the CEA (Commissariat l'nergie atomique) and is used mainly to solve problems of solid mechanics, fluid mechanics and heat transfers. Different types of discretization are available, among them finite element, finite volume and mixed hybrid finite element method. Cast3m is also a componant of the platteform Alliances (co-developped by Andra, CEA), which will be used by Andra for the safety calculation of an underground waste disposal in year 2004. We solve the Darcy equation for the water flow and a convection–diffusion transport equation for the Iodine 129 which escapes from a repository cave into the water. The water flow is calculated with a MHFE discretization. It is shown that this method provides sharp results even on relatively coarse grids. The convection–diffusion transport equation is discretized with FE (Finite Element), MHFE (Mixed Hybrid Finite Element) and FV (Finite Volume) methods. In our comparison, we point out the differences of these methods in term of accuracy, respect of the maximum principle and calculations cost. Neither the finite element nor the mixed hybrid finite element approach respects the maximum principle. This results in the presence of negative concentrations near the repository cave, whereas FV calculations respect the monotonicity. We show that mass lumping techniques suppress this problem but with strong restrictions on the grid. FE and MHFE approaches are more accurate than FV for the diffusion equation, but the overall results are equivalent since the advective terms are dominant in the far field and are discretized with centered schemes. We conclude by studying the influence of the grid: a very fine grid near the repository solves almost all the problems of monotonicity, without employing mass lumping techniques. We also observed a very important increase of the accuracy on a structured grid made up of rectangles. 相似文献
18.
基于渗流和管流耦合的管涌数值模拟 总被引：6，自引：0，他引：6
堤坝地基的渗透变形过程实际上是“土中水”转变为“水中土”的过程。在渗透变形发生的集中管涌通道区域,采用常规渗流分析理论,单纯增大管涌通道渗透系数的方法是不太合适的。在未发生渗透变形的区域,用常规渗流理论计算;在管涌通道区域,用管流理论,公共边界上两者之间水头相等、流量大小相等且方向相反,能够较好符合渗透变形的发展规律。为了适应计算过程中内部边界条件不断变化的特点,采用无网格法伽辽金法（element free Galerkin method,EFG）对渗流场进行计算。算例计算表明,这种渗流-管流耦合的方法能够模拟管涌通道绕过防渗墙等复杂的发展过程。 相似文献
19.
A new numerical method to solve the system of equations describing two phase flow in a Hele-Shaw cell is presented. It combines
a mixed finite element method, the method of subtraction of the singularity and a front tracking grid in a single computational
strategy. This choice of discretization techniques is well motivated by the difficulties present in the system of equations
and the physics of the problem. The new method was tested against analytical solutions and also by solving the Saffman–Taylor
viscous fingering problem for finite and infinite mobility ratios. In both cases convergence under mesh refinement is achieved
for the fingers developed from an initial sinusoidal interface. Finger splitting is observed for low values of the surface
tension and high mobility ratio. Different explanations, based in our results, are provided for this phenomenon.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
20.
In this paper, numerical simulation of 3-dimensional assemblies of 1000 polydisperse sphere particles using Discrete Element
Method (DEM) is used to study the liquefaction behaviour of granular materials. Numerical simulations of cyclic triaxial shear
tests under undrained conditions are performed at different confining pressures under constant strain amplitude. Results obtained
in these numerical simulations indicate that with increase in confining pressure there is an increase in liquefaction resistance. 相似文献