共查询到20条相似文献，搜索用时 15 毫秒
1.
Four classical geomechanics problems involving semi-infinite linear elastic media have been solved numerically using recently developed mapped infinite elements coupled to finite elements.The effect of the remoteness of the truncated boundary and the location of infinite element coupling on solution accuracy has been studied. The results of conventional analyses using finite elements over a relatively large but restricted region are compared to the coupled analyses. Comparison of the results shows that for the same number of degrees of freedom the performance of the coupled solutions is superior to the conventional approach with respect to accuracy of solution and computational efficiency. Finally, some general guidelines are proposed for the efficient numerical solution of these types of problems using the coupled finite/infinite element approach. 相似文献
2.
Large deformation soil behavior underpins the operation and performance for a wide range of key geotechnical structures and needs to be properly considered in their modeling, analysis, and design. The material point method (MPM) has gained increasing popularity recently over conventional numerical methods such as finite element method (FEM) in tackling large deformation problems. In this study, we present a novel hierarchical coupling scheme to integrate MPM with discrete element method (DEM) for multiscale modeling of large deformation in geomechanics. The MPM is employed to treat a typical boundary value problem that may experience large deformation, and the DEM is used to derive the nonlinear material response from small strain to finite strain required by MPM for each of its material points. The proposed coupling framework not only inherits the advantages of MPM in tackling large deformation engineering problems over the use of FEM (eg, no need for remeshing to avoid mesh distortion in FEM), but also helps avoid the need for complicated, phenomenological assumptions on constitutive material models for soil exhibiting high nonlinearity at finite strain. The proposed framework lends great convenience for us to relate rich grain-scale information and key micromechanical mechanisms to macroscopic observations of granular soils over all deformation levels, from initial small-strain stage en route to large deformation regime before failure. Several classic geomechanics examples are used to demonstrate the key features the new MPM/DEM framework can offer on large deformation simulations, including biaxial compression test, rigid footing, soil-pipe interaction, and soil column collapse. 相似文献
3.
We pay a revisit to some classical geomechanics problems using a novel computational multiscale modelling approach. The multiscale approach employs a hierarchical coupling of the finite element method (FEM) and the discrete element method. It solves a boundary value problem at the continuum scale by FEM and derives the material point response from the discrete element method simulation attached to each Gauss point of the FEM mesh. The multiscale modelling framework not only helps successfully bypass phenomenological constitutive assumptions as required in conventional modelling approaches but also facilitates effective cross‐scale interpretation and understanding of soil behaviour. We examine the classical retaining wall and footing problems by this method and demonstrate that the simulating results can be well validated and verified by their analytical solutions. Furthermore, the study sheds novel multiscale insights into these classical problems and offers a new tool for geotechnical engineers to design and analyse geotechnical applications based directly upon particle‐level information of soils. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
4.
B. Nath 《Computers and Geotechnics》1990,10(4):265-285
What may be called a ‘Continuum’ method of finite element analysis is used in this paper to predict the behaviour of a pile during driving. In this both the pile and the soil are treated as two distinct parts of the same solid continuum, but with different properties. The behaviour of the soil medium, assumed to be semi-infinite and nonlinear, is represented by the hyperbolic stress-strain relationship. Discretizing the pile-soil system in turn by conventional axi-symmetric and mapped finite elements, the problem is solved in the time-domain using the central difference scheme.
The example considered is that of a fully embedded steel pipe pile for which both field test data and Wave Equation solution are available. Results show that: the Continuum Method is capable of a greater degree of accuracy than the conventional Wave Equation Method, but it is far more expensive than the latter in terms of computational effort needed; the effects of radiation damping and wave dispersion in the soil are found to be small; and the mapped finite elements give significantly better results than conventional elements. 相似文献
5.
A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
6.
Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method 
下载免费PDF全文
下载免费PDF全文 The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
7.
Twana Kamal Haji Asaad Faramarzi Nicole Metje David Chapman Farough Rahimzadeh 《国际地质力学数值与分析法杂志》2020,44(3):418-431
The accurate modelling of gravity is of crucial importance for a variety of issues including, but not restricted to, the identification of buried objects. Gravity is an unbounded problem, which causes challenges when applying numerical models, i..e.., it results in computational difficulties when specifying the relevant boundary conditions. In order to address this, previous research has tended to generate artificial boundary conditions, e.g., truncating the simulated domain and adding many unrealistic zero-density layers, which introduces more unknown parameters and unnecessarily excessive computational time. In order to overcome such inaccuracies, this paper proposes an innovative development of the finite element modelling technique, which represents a step change in the field of gravity forward modelling. A comprehensive formulation of an infinite element to reproduce the far-field boundary effect using only one layer of infinite elements is presented. The developed model considerably reduces the computational time while obtaining high degrees of accuracy. The model is validated against the exact solution of the problem, and its results show an excellent performance. The proposed method can significantly improve the postprocessing and interpretation stages of data analysis relevant to micro-gravity sensors. The new method is applied to subsurface civil engineering although its applicability is manifold. 相似文献
8.
Unbounded plane stress and plane strain domains subjected to static loading undergo infinite displacements, even when the zero displacement boundary condition at infinity is enforced. However, the stress and strain fields are well behaved, and are of practical interest. This causes significant difficulty when analysis is attempted using displacement‐based numerical methods, such as the finite‐element method. To circumvent this difficulty problems of this nature are often changed subtly before analysis to limit the displacements to finite values. Such a process is unsatisfactory, as it distorts the solution in some way, and may lead to a stiffness matrix that is nearly singular. In this paper, the semi‐analytical scaled boundary finite‐element method is extended to permit the analysis of such problems without requiring any modification of the problem itself. This is possible because the governing differential equations are solved analytically in the radial direction. The displacement solutions so obtained include an infinite component, but relative motion between any two points in the unbounded domain can be computed accurately. No small arbitrary constants are introduced, no arbitrary truncation of the domain is performed, and no ill‐conditioned matrices are inverted. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
9.
A set of mapping functions in the form of convergent series for an infinite element, which is capable to include the infinitely distanced constant head boundary condition from the area of disturbance (e.g. pumping), is proposed based on the asymptotic far-field behaviour of typical seepage flow problems. The derived mapping functions have been successfully used in three-dimensional point symmetric, two-dimensional axi-symmetric and one-dimensional unidirectional flow for the fixed head boundary at infinite distance. The result shows excellent agreement with analytical solution. For the first time, the mapping function of an infinite element is presented in the form of a convergent series. The infinite elements are really capable of reducing the cost and efficiency of conventional finite element analysis. Finally, a figure is also proposed to indicate the required size of the near field to obtain accurate drawdown at specified locations based on some calculations for two-dimensional radial flow case. 相似文献
10.
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. 相似文献
11.
S. K. Sharan 《国际地质力学数值与分析法杂志》1989,13(5):565-570
A highly efficient technique is presented for the finite element analysis of stresses around underground openings surrounded by an infinite extent of linearly elastic rock mass. The effect of unbounded rock is incorporated in the analysis by using elastic supports around the finite element model. Expressions for stiffnesses of the elastic supports are derived and these are found to depend on the location of elastic supports, the shear modulus and Poisson's ratio of the rock, and the ratio of horizontal to vertical initial stresses. With the use of the proposed technique, the extent of the finite domain to be considered in the analysis is highly reduced. This results in a great computational advantage. The other principal merit of the proposed technique is that a standard finite element code for stress analysis may be used without making any modification. Results of some numerical tests are reported to demonstrate the effectiveness and efficiency of the proposed method. The technique has the potential of being applied to more complex problems of unbounded domains in geomechanics. 相似文献
12.
B. Nath 《国际地质力学数值与分析法杂志》1981,5(2):139-163
A novel finite element method has been proposed in this paper for the solution of seepage problems economically and accurately. In this method the governing equation and the prescribed boundary conditions are transformed so that they refer to a suitable logarithmically condensed ‘image’ space; the physical problem domain is also mapped into the image space. The transformed equation is then solved in the image space using standard finite elements, subject to the transformed boundary conditions. Because physical space is logarithmically condensed in the image space, the proposed method is capable of dealing with large or very large aspect ratio seepage problems economically and accurately. The validity of the method has been demonstrated by means of a number of examples including anisotropy and non-linearity. In all cases an excellent degree of accuracy was achieved, efficiently and economically. 相似文献
13.
In this paper, a large deformation finite element (LDFE) approach termed ‘remeshing and interpolation technique with small strain (RITSS)’ is extended from static to dynamic soil-structure interaction applications. In addition, a technique termed ‘element addition’ is developed to improve the computational efficiency of both static and dynamic LDFE analyses that involve moving boundaries. The RITSS approach is based on frequent mesh generation to avoid element distortion. In dynamic RITSS, the field variables mapped from the old to the new mesh involve not only the stresses and material properties, but also the nodal velocities and accelerations. Using the element addition technique, new soil elements are attached to the domain boundaries periodically when the soil near the boundaries becomes affected by large displacements of the structure. The procedures of this Abaqus-based dynamic LDFE analysis and element addition technique are detailed, and the robustness of the techniques is validated and assessed through three example analyses: penetration of a flat footing into a half-space and movement of rigid and deformable landslides down slopes. 相似文献
14.
A new finite element scheme is proposed, in this paper, for solving two-dimensional wave propagation problems in multilayered soils resting on a rigid base. The multilayered soils are treated as multiple horizontal layers of lateral infinite extension in geometry. Since these horizontal layers can be truncated by two artificially truncated vertical boundaries, two high-order artificial boundary conditions are applied for propagating the incoming waves from the interior domain into the far field of the system. Both the semi-analytical method and the truncated boundary migration procedure are used to derive the high-order artificial boundary conditions, which are comprised of a physically meaningful dashpot and a generalized energy absorber. The main advantage of using the proposed finite element scheme is that the derived artificial boundary condition can be straightforwardly implemented in the finite element analysis, without violating the band/sparse structure of the conventional finite element equation. The related numerical examples have demonstrated that the proposed finite element scheme is of high accuracy in dealing with wave propagation problems in multiple horizontal layers. 相似文献
15.
The coupled discrete element method and lattice Boltzmann method (DEMLBM) has increasingly drawn attention of researchers in geomechanics due to its mesoscopic nature since 2000. Immersed boundary method (IBM) and immersed moving boundary (IMB) are two popular schemes for coupling fluid particle in DEMLBM. This work aims at coupling DEM and LBM using the latest IBM algorithm and investigating its accuracy, computational efficiency, and applicability. Two benchmark tests, interstitial fluid flow in an ideal packing and single particle sedimentation in viscous fluid, are carried out to demonstrate the accuracy of IBM through semi-empirical Ergun equation, finite element method (FEM), and IMB. Then, simulations of particle migration with relatively large velocity in Poiseuille flow are utilized to address limitations of IBM in DEMLBM modeling. In addition, advantages and deficiencies of IBM are discussed and compared with IMB. It is found that the accuracy of IBM can be only guaranteed when sufficient boundary points are used and it is not suitable for geomechanical problems involving large fluid or particle velocity. 相似文献
16.
A local artificial boundary for transient seepage problems with unsteady boundary conditions in unbounded domains 下载免费PDF全文
下载免费PDF全文 Numerical analysis of transient seepage in unbounded domains with unsteady boundary conditions requires a more sophisticated artificial boundary approach to deal with the infinite character of the domain. To that end, a local artificial boundary is established by simplifying a global artificial boundary. The global artificial boundary conditions (ABCs) at the truncated boundary are derived from analytical solutions for one‐dimensional axisymmetric diffusion problems. By applying Laplace transforms and introducing some specially defined auxiliary variables, the global ABCs are simplified to local ABCs to significantly enhance the computational efficiency. The proposed local ABCs are implemented in a finite element computer program so that the solutions to various seepage problems can be calculated. The proposed approach is first verified by the computation of a one‐dimensional radial flow problem and then tentatively applied to more general two‐dimensional cylindrical problems and planar problems. The solutions obtained using the local ABCs are compared with those obtained using a large element mesh and using a previously proposed local boundary. This comparison demonstrates the satisfactory performance and obvious superiority of the newly established boundary to the other local boundary. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
17.
An infinite element is presented to treat wave propagation problems in unbounded saturated porous media. The porous media is modeled by Biot's theory. Conventional finite elements are used to model the near field, whereas infinite elements are used to represent the behavior of the far field. They are constructed in such a way that the Sommerfeld radiation condition is fulfilled, i.e. the waves decay with distance and are not reflected at infinity. To provide the wave information the infinite elements are formulated in Laplace domain. The time domain solution is obtained by using the convolution quadrature method as the inverse Laplace transformation. The temporal behavior of the near field is calculated using standard time integration schemes, e.g. the Newmark method. Finally, the near and far field are combined using a substructure technique for any time step. The accuracy as well as the necessity of the proposed infinite elements, when unbounded domains are considered, is demonstrated by different examples. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
18.
The paper presents total-stress numerical analyses of large-displacement soil-structure interaction problems in geomechanics using the Particle Finite Element Method (PFEM). This method is characterized by frequent remeshing and the use of low order finite elements to evaluate the solution. Several important features of the method are: (i) a mixed formulation (displacement-mean pressure) stabilized numerically to alleviate the volumetric locking effects that are characteristic of low order elements when the medium is incompressible, (ii) a penalty method to prescribe the contact constraints between a rigid body and a deformable media combined with an implicit scheme to solve the tangential contact constraint, (iii) an explicit algorithm with adaptive substepping and correction of the yield surface drift to integrate the finite-strain multiplicative elasto-plastic constitutive relationship, and (iv) the mapping schemes to transfer information between successive discretizations. The performance of the method is demonstrated by several numerical examples, of increasing complexity, ranging from the insertion of a rigid strip footing to a rough cone penetration test. It is shown that the proposed method requires fewer computational resources than other numerical approaches addressing the same type of problems. 相似文献
19.
20.
A new development of the scaled boundary finite element method for wave motion in layered half-space
The scaled boundary finite element method (SBFEM) developed by Wolf and Song has shown certain parallels to the finite element method (FEM) and boundary element method (BEM). Because of its semi-analytical nature, SBFEM is particularly suitable for the analysis of wave propagation in unbounded domains. This paper makes a certain modification of the standard SBFEM. A new idea of scaling surface instead of a scaling center is introduced to formulate the governing SBFE equations for the analysis of wave propagation in multilayered half-space, which leads to simplifying the modeling and saving considerably the computational effort. In addition, by employing the proposed approach, some problems encountered in engineering practice, which are difficult to deal with by the conventional SBFEM, for example, 3D foundation impedance on half-space with irregular geographical features, can be effectively solved. The proposed approach also helps to simplify the solution of shell structures. Numerical examples are provided to validate the accuracy and efficiency of the proposed approach. 相似文献