共查询到20条相似文献，搜索用时 218 毫秒
1.
S. Amir Reza Beyabanaki Ahmad Jafari S. Omid Reza Biabanaki M. Ronald Yeung 《Computers and Geotechnics》2009
This paper presents a new numerical model that can add a finite element mesh into each block of the three-dimensional discontinuous deformation analysis (3-D DDA), originally developed by Gen-hua Shi. The main objectives of this research are to enhance DDA block’s deformability. Formulations of stiffness and force matrices in 3-D DDA with conventional Trilinear (8-node) and Serendipity (20-node) hexahedral isoparametric finite elements meshed block system due to elastic stress, initial stress, point load, body force, displacement constraints, inertia force, normal and shear contact forces are derived in detail for program coding. The program code for the Trilinear and Serendipity hexahedron elements have been developed, and it has been applied to some examples to show the advantages achieved when finite element is associated with 3-D DDA to handle problems under large displacements and deformations. Results calculated for the same models by use of the original 3-D DDA are far from the theoretical solutions while the results of new numerical model are quite good in agreement with theoretical solutions; however, for the Trilinear elements, more number of elements are needed. 相似文献
2.
The discontinuous deformation analysis (DDA) with second‐order displacement functions was derived based on six‐node triangular mesh in order to satisfy the requirement for the accurate calculations in practical applications. The matrices of equilibrium equations for the second‐order DDA were given in detail for program coding. By close comparison with widely used finite element method and closed form solutions, the advantages of the modified DDA were illustrated. The program coding was carried out in C++ environment and the new code applied to three examples with known analytical solutions. A very good agreement was achieved between the analytical and numerical results produced by the modified DDA code. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
3.
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. 相似文献
4.
Shunde Yin 《国际地质力学数值与分析法杂志》2013,37(15):2523-2538
Integration of poromechanics and fracture mechanics plays an important role in understanding a series of thermal fracturing phenomena in subsurface porous media such as cold water flooding for enhanced oil recovery, produced‐water reinjection for waste disposal, cold water injection for geothermal energy extraction, and CO2 injection for geosequestration. Thermal fracturing modeling is important to prevent the potential risks when fractures propagate into undesired zones, and it involves the coupling of heat transfer, mass transport, and stress change as well as the fracture propagation. Analytical method, finite element method, and finite difference method as well as boundary element method have been used to perform the thermal fracturing modeling considering different degrees and combinations of coupling. In this paper, extended finite element method is employed for the thermal fracturing modeling in a fully coupled fashion with remeshing avoided, and the stabilized finite element method is employed to account for the convection‐dominated heat transfer in the fracturing process with numerical oscillation circumvented. With the thermal fracturing model, a hypothetical numerical experiment on cold water injection into a deep warm aquifer is conducted. Results show that parameters such as injection rate, injection temperature, aquifer stiffness, and permeability can affect the fracture development in different ways and extended finite element method and stabilized finite element method provide effective tools for thermal fracturing simulation. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
5.
We present a stabilized extended finite element formulation to simulate the hydraulic fracturing process in an elasto‐plastic medium. The fracture propagation process is governed by a cohesive fracture model, where a trilinear traction‐separation law is used to describe normal contact, cohesion and strength softening on the fracture face. Fluid flow inside the fracture channel is governed by the lubrication equation, and the flow rate is related to the fluid pressure gradient by the ‘cubic’ law. Fluid leak off happens only in the normal direction and is assumed to be governed by the Carter's leak‐off model. We propose a ‘local’ U‐P (displacement‐pressure) formulation to discretize the fluid‐solid coupled system, where volume shape functions are used to interpolate the fluid pressure field on the fracture face. The ‘local’ U‐P approach is compatible with the extended finite element framework, and a separate mesh is not required to describe the fluid flow. The coupled system of equations is solved iteratively by the standard Newton‐Raphson method. We identify instability issues associated with the fluid flow inside the fracture channel, and use the polynomial pressure projection method to reduce the pressure oscillations resulting from the instability. Numerical examples demonstrate that the proposed framework is effective in modeling 3D hydraulic fracture propagation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
6.
A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows
The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
7.
非连续变形分析方法（DDA）是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。 相似文献
8.
传统的非连续变形分析法（DDA）法采用简单的线性位移模式计算效率高,描述大块体的高阶多项式位移模式在一定程度保留了该特点,并提高了计算精度。近年来流行的耦合有限元、自然单元的DDA法实质上是引入相应的插值形函数构成块体位移函数,计算相对低效,但具有计算更精细、更容易施加边界条件等优点。为结合传统DDA法与DDA耦合法各自的优点,建立了一种同时利用传统DDA法线性位移模式与耦合型DDA法非线性位移模式的混合法。该方法非线性模式主要针对大块体,采用了自然单元插值,缘于其具有一定无网格特征,且效率比有限元高。建立了混合模式下的整体矩阵并推导出接触等因素刚度子矩阵和荷载子向量的具体表达式。该方法建模更加方便合理,计算精度、效率介于线性模式的传统DDA法和非线性位移模式的耦合法之间。通过基本算例验证了混合法的有效性,并给出了节理围岩-隧道衬砌整体分析模型的计算结果,体现了新方法的优越性。 相似文献
9.
Discontinuous deformation analysis (DDA), a discrete numerical analysis method, is used to simulate the behaviour of falling rock by applying a linear displacement function in the computations. However, when a block rotates, this linear function causes a change in block size called the free expansion phenomenon. In addition, this free expansion results in contact identification problems when the rotating blocks are close to each other. To solve this problem of misjudgment and to obtain a more precise simulation of the falling rock, a new method called Post‐Contact Adjustment Method has been developed and applied to the program. The basic procedure of this new method can be divided into three stages: using the linear displacement function to generate the global matrix, introducing the non‐linear displacement function to the contact identification, and applying it to update the co‐ordinates of block vertices. This new method can be easily applied to the original DDA program, demonstrating better contact identification and size conservation results for falling rock problems than the original program. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
10.
Recent attempts to couple discontinuous deformation analysis (DDA) blocks with finite element meshes have been regarded as fruitful. The hybrid model has been proven to be reliable in modeling continuum–discontinuum problems. In this paper, a hydraulic crack initiation–propagation extension complete with a coupled hydro-mechanical analysis algorithm is introduced to the hybrid model for the simulation of fracturing problems initiated by hydraulic pressure. Several simulation cases, of which the results appear to be consistent with existing theories of hydraulic fracturing mechanisms, are presented to demonstrate the capability of the proposed algorithm in modeling hydraulic fracturing processes. 相似文献
11.
Softening solids are analysed under impact loading using a new numerical method which allows displacement discontinuities to propagate arbitrarily through a finite element mesh. The Dirac‐delta distributions that arise in the strain field of classical continuum theory in the presence of strain softening are interpreted as discontinuities in the displacement field. A new finite element procedure with Heaviside jumps added to the underlying displacement interpolation basis is able to capture displacement jumps independent of the spatial discretisation. The amplitudes of displacement jumps are represented by extra degrees of freedom at existing nodes. Numerical results for mode‐I and mode‐II failure due to impact loading are presented. The numerical results highlight the objectivity of the approach with respect to spatial discretisation under dynamic loading conditions. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
12.
In the original formulation of DDA by Shi, a linear displacement function term is used. This has the limitations of uncontrolled block and stress distortion due to rigid body rotation. In the present paper, the authors propose a new iterative method which can avoid the distortion due to the rotation even when the rotation or number of time‐step is large. Furthermore, the authors propose a simple internal discretization scheme which is applicable for both concave and convex polygon which is particularly important for a large block. The stress and strain distribution with a large block can be obtained with ease under this scheme. The numerical examples as shown have demonstrated the advantages of the present proposal in DDA analysis. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
13.
The driving response of thin‐walled open‐ended piles is studied using numerical simulation of the wave propagation inside the soil plug and the pile. An elastic finite element analysis is carried out to identify the stress wave propagation in the vicinity of the pile toe. It is found that the shear stress wave has the highest magnitude above the bottom of the soil plug. Below the bottom of the soil plug, the vertical stress wave has the highest magnitude. Although the shear stress wave propagating in the radial direction is similar in magnitude to the vertical stress wave at the bottom of the soil plug, it decays rapidly while travelling downwards. The highest vertical stress at the bottom of the soil plug appears after the vertical stress wave interacts with the shear stress wave travelling in the radial direction. Initially, the vertical stress wave propagates with the dilation wave velocity in both the radial and vertical directions. After it interacts with the shear stress wave, the vertical stress wave starts to propagate with the shear wave velocity in the radial direction and with the axial wave velocity downwards. It is concluded that at the bottom of the soil plug, the interaction between the waves travelling in radial and vertical directions is important. The capabilities of several one‐dimensional pile‐in‐pile models to reproduce the driving response given by a two‐dimensional axisymmetric finite element model is studied. It is seen that when the base of the soil plug fails, a one‐dimensional pile‐in‐pile model can be used to achieve results in agreement with the finite element model. However, when the pile is unplugged, where the base of the soil plug does not fail, a reduced finite element mesh that permits the radial wave propagation inside the soil plug must be used. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
14.
15.
A finite element, variable mesh analysis of unconfined steady-state seepage problems is presented based on a nonlinear programming algorithm. It is shown that the minimization of an objective function which merely represents a measure of the total flux leaving or entering the mesh at the free surface nodes (except those that belong also to pervious boundaries) does not permit a unique definition of the free surface geometry. This problem, which is apparently related to the numerical instabilities often met when using variable mesh approaches, can be eliminated by adding to the objective function a term representing a sort of overall ‘regularity’ condition for the shape of the free surface. The modified solution procedure turns out to be stable and able to provide meaningful results for practical problems even when rather coarse meshes are adopted. 相似文献
16.
This paper presents a ‘Eulerian‐like’ finite element technique to simulate the large accumulated displacements of piles subjected to multiple hammer blows. For each hammer blow, results are obtained using a standard small strain finite element model and, at the end of each hammer blow, material flow is taken into account with reference to a fixed finite element mesh. Residual stresses calculated at the Gauss integration points of the deformed finite element mesh are mapped on to the fixed finite element mesh, and these stresses are used as initial stresses for the next hammer blow. At the end of each hammer blow, stiffness and mass matrices are recalculated for the volume of material remaining inside the fixed finite element mesh. Results obtained with and without allowing material to flow through the fixed mesh are compared for several hammer blows. Build up of residual stresses, soil flow and yielded points around the pile are presented for plugged, partially‐plugged and unplugged piles. Using the new finite element technique, the driving of a pile from the soil surface is studied. The ability to analyse this and other large deformation problems is the main advantage of the new finite element technique. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
17.
Coupled hydromechanical‐fracture simulations of nonplanar three‐dimensional hydraulic fracture propagation 
下载免费PDF全文
下载免费PDF全文 This paper presents an algorithm and a fully coupled hydromechanical‐fracture formulation for the simulation of three‐dimensional nonplanar hydraulic fracture propagation. The propagation algorithm automatically estimates the magnitude of time steps such that a regularized form of Irwin's criterion is satisfied along the predicted 3‐D fracture front at every fracture propagation step. A generalized finite element method is used for the discretization of elasticity equations governing the deformation of the rock, and a finite element method is adopted for the solution of the fluid flow equation on the basis of Poiseuille's cubic law. Adaptive mesh refinement is used for discretization error control, leading to significantly fewer degrees of freedom than available nonadaptive methods. An efficient computational scheme to handle nonlinear time‐dependent problems with adaptive mesh refinement is presented. Explicit fracture surface representations are used to avoid mapping of 3‐D solutions between generalized finite element method meshes. Examples demonstrating the accuracy, robustness, and computational efficiency of the proposed formulation, regularized Irwin's criterion, and propagation algorithm are presented. 相似文献
18.
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. 相似文献
19.
20.
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. 相似文献