共查询到20条相似文献，搜索用时 31 毫秒
1.
Soil–structure interaction problems are commonly encountered in engineering practice, and the resulting linear systems of equations are difficult to solve due to the significant material stiffness contrast. In this study, a novel partitioned block preconditioner in conjunction with the Krylov subspace iterative method symmetric quasiminimal residual is proposed to solve such linear equations. The performance of these investigated preconditioners is evaluated and compared on both the CPU architecture and the hybrid CPU–graphics processing units (GPU) computing environment. On the hybrid CPU–GPU computing platform, the capability of GPU in parallel implementation and high-intensity floating point operations is exploited to accelerate the iterative solutions, and particular attention is paid to the matrix–vector multiplications involved in the iterative process. Based on a pile-group foundation example and a tunneling example, numerical results show that the partitioned block preconditioners investigated are very efficient for the soil–structure interaction problems. However, their comparative performances may apparently depend on the computer architecture. When the CPU computer architecture is used, the novel partitioned block symmetric successive over-relaxation preconditioner appears to be the most efficient, but when the hybrid CPU–GPU computer architecture is adopted, it is shown that the inexact block diagonal preconditioners embedded with simple diagonal approximation to the soil block outperform the others. 相似文献
2.
The finite element (FE) simulation of large‐scale soil–structure interaction problems (e.g. piled‐raft, tunnelling, and excavation) typically involves structural and geomaterials with significant differences in stiffness and permeability. The symmetric quasi‐minimal residual solver coupled with recently developed generalized Jacobi, modified symmetric successive over‐relaxation (MSSOR), or standard incomplete LU factorization (ILU) preconditioners can be ineffective for this class of problems. Inexact block diagonal preconditioners that are inexpensive approximations of the theoretical form are systematically evaluated for mitigating the coupled adverse effects because of such heterogeneous material properties (stiffness and permeability) and because of the percentage of the structural component in the system in piled‐raft foundations. Such mitigation led the proposed preconditioners to offer a significant saving in runtime (up to more than 10 times faster) in comparison with generalized Jacobi, modified symmetric successive over‐relaxation, and ILU preconditioners in simulating piled‐raft foundations. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
3.
An efficient finite–discrete element method applicable for the analysis of quasi‐static nonlinear soil–structure interaction problems involving large deformations in three‐dimensional space was presented in this paper. The present method differs from previous approaches in that the use of very fine mesh and small time steps was not needed to stabilize the calculation. The domain involving the large displacement was modeled using discrete elements, whereas the rest of the domain was modeled using finite elements. Forces acting on the discrete and finite elements were related by introducing interface elements at the boundary of the two domains. To improve the stability of the developed method, we used explicit time integration with different damping schemes applied to each domain to relax the system and to reach stability condition. With appropriate damping schemes, a relatively coarse finite element mesh can be used, resulting in significant savings in the computation time. The proposed algorithm was validated using three different benchmark problems, and the numerical results were compared with existing analytical and numerical solutions. The algorithm performance in solving practical soil–structure interaction problems was also investigated by simulating a large‐scale soft ground tunneling problem involving soil loss near an existing lining. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
4.
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. 相似文献
5.
This paper describes a soil‐structure coupling method to simulate blast loading in soil and structure response. For the last decade, simulation of soil behavior under blast loading and its interaction with semi buried structure in soil becomes the focus of computational engineering in civil and mechanical engineering communities. In current design practice, soil‐structure interaction analysis often assumes linear elastic properties of the soil and uses small displacement theory. However, there are numerous problems, which require a more advanced approach that account for soil‐structure interaction and appropriate constitutive models for soil. In simplified approaches, the effect of soil on structure is considered using spring‐dashpot‐mass system, and the blast loading is modeled using linearly decaying pressure–time history based on equivalent trinitrotoluene and standoff distance, using ConWep, a computer program based on semi‐empirical equations. This strategy is very efficient from a CPU time computing point of view but may not provide accurate results for the dynamic response of the structure, because of its significant limitations, mainly when soil behavior is strongly nonlinear and when the buried charge is close to the structure. In this paper, both soil and explosive are modeled using solid elements with a constitutive material law for soil, and a Jones–Wilkins–Lee equation of state for explosive. One of the problems we have encountered when solving fluid structure interaction problems is the high mesh distortion at the contact interface because of high fluid nodal displacements and velocities. Similar problems have been encountered in soil structure interaction problems. To prevent high mesh distortion for soil, a new coupling algorithm is performed at the soil structure interface for structure loading. The coupling method is commonly used for fluid structure interaction problems in automotive and aerospace industry for fuel sloshing tank, and bird impact problems, but rarely used for soil structure interaction problems, where Lagrangian contact type algorithms are still dominant. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
6.
Preconditioners in computational geomechanics: A survey 总被引：1，自引：0，他引：1
The finite element (FE) solution of geomechanical problems in realistic settings raises a few numerical issues depending on the actual process addressed by the analysis. There are two basic problems where the linear solver efficiency may play a crucial role: 1. fully coupled consolidation and 2. faulted uncoupled consolidation. A class of general solvers becoming increasingly popular relies on the Krylov subspace (or Conjugate Gradient‐like) methods, provided that an efficient preconditioner is available. For both problems mentioned above, the possible preconditioners include the diagonal scaling (DS), the Incomplete LU decomposition (ILU), the mixed constraint preconditioning (MCP) and the multilevel incomplete factorization (MIF). The development and the performance of these algorithms have been the topic of several recent works. The present paper aims at providing a survey of the preconditioners available to date in computational geomechanics. In particular, a review and a critical discussion of DS, ILU, MCP and MIF are given along with some comparative numerical results. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
7.
Constraint preconditioners have proved very efficient for the solution of ill-conditioned finite element (FE) coupled consolidation
problems in a sequential computing environment. Their implementation on parallel computers, however, is not straightforward
because of their inherent sequentiality. The present paper describes a novel parallel inexact constraint preconditioner (ParICP)
for the efficient solution of linear algebraic systems arising from the FE discretization of the coupled poro-elasticity equations.
The ParICP implementation is based on the use of the block factorized sparse approximate inverse incomplete Cholesky preconditioner,
which is a very recent and effective development for the parallel preconditioning of symmetric positive definite matrices.
The ParICP performance is experimented with in real 3D coupled consolidation problems, proving a scalable and efficient implementation
of the constraint preconditioning for high-performance computing. ParICP appears to be a very robust algorithm for solving
ill-conditioned large-size coupled models in a parallel computing environment. 相似文献
8.
This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
9.
In this paper we present a study of the performance of sparse iterative solvers regarding the resolution of three-dimensional and non-linear problems encountered in soil/structure interaction. It is composed of two parts. In the first one, we present briefly iterative methods and preconditioners used in this study, then we analyse their performance on three soil/structure interaction problems: a shallow foundation under a vertical loading, a single pile subjected to a lateral loading and the construction of a lined tunnel in a soft soil. Tests are performed assuming an elastic–perfectly plastic constitutive law for the soil material with a non-associated Mohr–Coulomb flow rule. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
10.
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. 相似文献
11.
12.
Contact between stiff structural elements and soil is encountered in many applications in geotechnical engineering. Modelling of such contact is challenging as it often involves impact that would lead to large deformation and failure of the soil. The Material Point Method (MPM) is a mesh‐free method that has been applied to simulate such phenomena. However, the frictional contact algorithm commonly used in MPM only supports Coulomb friction and cannot model fully or partially rough contact conditions in terms of geotechnical engineering. Moreover, because of very different stiffness of contacting materials, the contact force predicted by the previous frictional contact algorithms usually suffers from severe oscillation when applied in structure–soil interaction. This paper presents a new contact algorithm, termed Geo‐contact, designed for geotechnical engineering. In Geo‐contact, a penalty function is incorporated to reduce the oscillation in contact computation, and a limited shear stress is specified along the contact interface. The proposed Geo‐contact algorithm has been implemented to simulate smooth, partially rough and rough contact in typical large deformation penetration problems. The resistance–displacement curves obtained using the Geo‐contact are compared with analytical solutions of limit analysis and large deformation finite element results to verify the accuracy and robustness of the proposed contact algorithm. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
13.
Foundation settlements and soil–structure interaction are important problems to structural and geotechnical engineers. This study introduces a novel elastoplastic three‐degree‐of‐freedom medium which models foundations settlements under combined loadings. A soil–structure interaction problem can then be solved by replacing the soil mass with this three‐degree‐of‐freedom elastoplastic medium, thus reducing significantly the size of the problem. The model was developed by extending the classical plasticity concepts to the force‐deformation level. Its ability to predict foundation deformations was evaluated using finite element solutions of a typical shallow foundation problem and was found reasonably accurate while producing significant time savings. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
14.
15.
This paper presents a superposition method expanded for computing impedance functions (IFs) of inclined‐pile groups. Closed‐form solutions for obtaining horizontal, vertical, and rocking IFs, estimated by using pile‐to‐pile interaction factors, are proposed. IFs of solitary inclined piles, crossed IFs, and explicit incorporation of compatibility conditions for pile‐head movements are also appropriately taken into consideration. All of these factors should be known in advance and will be computed and shown for the most relevant cases. The accuracy of the proposed closed‐form solutions is verified for 2 × 2 and 3 × 3 square inclined‐pile groups embedded in an isotropic viscoelastic homogeneous half‐space soil medium, with hysteretic damping. The pile‐to‐pile interaction factors are computed by means of a three‐dimensional time‐harmonic boundary elements–finite elements coupling formulation. The results indicate that the IFs obtained from the proposed method are in good agreement with those obtained from the coupling formulation. Furthermore, crossed vertical‐rocking IFs of solitary piles need to be appropriately considered for obtaining rocking IFs when the number of piles is small. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
16.
The focus of this work is efficient solution methods for mixed finite element models of variably saturated fluid flow through
deformable porous media. In particular, we examine preconditioning techniques to accelerate the convergence of implicit Newton–Krylov
solvers. We highlight an approach in which preconditioners are built from block-factorizations of the coupled system. The
key result of the work is the identification of effective preconditioners for the various sub-problems that appear within
the block decomposition. We use numerical examples drawn from both linear and nonlinear hydromechanical models to test the
robustness and scalability of the proposed methods. Results demonstrate that an algebraic multigrid variant of the block preconditioner
leads to mesh-independent convergence, good parallel efficiency, and insensitivity to the material parameters of the medium. 相似文献
17.
Large‐scale simulations of flow in deformable porous media require efficient iterative methods for solving the involved systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in geological applications, such as basin evolution at regional scales. The success of iterative methods for this type of problems depends strongly on finding effective preconditioners. This paper investigates how the block‐structured matrix system arising from single‐phase flow in elastic porous media should be preconditioned, in particular for highly discontinuous permeability and significant jumps in elastic properties. The most promising preconditioner combines algebraic multigrid with a Schur complement‐based exact block decomposition. The paper compares numerous block preconditioners with the aim of providing guidelines on how to formulate efficient preconditioners. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
18.
Guang Y. Sheu 《国际地质力学数值与分析法杂志》2006,30(8):823-842
This study presents a preliminary development of a direct back analysis procedure by the meshless local Petrov–Galerkin (MLPG) method and Bayesian statistics and the application of resulting procedure to characterize soil properties using laboratory results. As compared to direct back analyses that are based on the finite element method (FEM), it is intended to show that MLPG‐based direct back analyses may be more suitable for some kinds of characterization problems; for example, involving a complex subsurface stratification or the characterization of soil properties of just an inclusion of a soil profile. The existing MLPG method is first slightly modified to analyse time‐dependent problems. Using the resulting method, quantities to be characterized are evaluated so that they give numerical results as close to measured data as possible. The Akaike information criterion is introduced for simplifying the evaluation. A one‐dimensional finite strain consolidation problem is introduced to do an error analysis for prediction by the proposed MLPG method. Another example illustrates experiences of performing an MLPG‐based direct back analysis. Comparison of MLPG‐based and FEM‐based direct back analyses is taken. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
19.
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. 相似文献
20.
An iterative process based upon a hybrid ‘residual force’ method is presented for solving elasto–plastic soil–structure interaction problems. In this approach the soil and the structure are treated as separate bodies and related only by compatibility of displacements and equilibrium of forces at the soil–structure interface. This scheme enables a significant improvement in numerical stability and rate of convergence over the conventional initial stress method. It is also shown that various interface conditions such as shear failure, slip and breakaway, and frictional and dilatant behaviour can be readily accounted for. Some practical aspects associated with the proposed scheme are emphasized for a number of numerical examples. 相似文献