共查询到20条相似文献，搜索用时 919 毫秒
1.
Direct, partitioned, and projected (conjugate gradient‐like) solution approaches are compared on unsymmetric indefinite systems arising from the finite element integration of coupled consolidation equations. The direct method is used in its most recent and computationally efficient implementations of the Harwell Software Library. The partitioned approach designed for coupled problems is especially attractive as it addresses two separate positive definite problems of a smaller size that can be solved by symmetric conjugate gradients. However, it may stagnate and when converging it does not prove competitive with a global projection method such as Bi‐CGSTAB, which may take full advantage of its flexibility in working on scaled and reordered equations, and thus may greatly improve its computational performance in terms of both robustness and convergence rate. The Bi‐CGSTAB superiority to the other approaches is discussed and demonstrated with a few representative examples in two‐dimensional (2‐D) and three‐dimensional (3‐D) coupled consolidation problems. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
2.
Preconditioned projection (or conjugate gradient like) methods are increasingly used for the accurate and efficient solution to finite element (FE) coupled consolidation equations. Theory indicates that preliminary row/column scaling does not affect the eigenspectrum of the iteration matrix controlling convergence as long as the preconditioner relies on the incomplete factorization of the FE coefficient matrix. However, computational experience with mid‐large size problems shows that the above inexpensive operation can significantly accelerate the solver convergence, and to a minor extent also improve the final accuracy, as a result of a better solver stability to the accumulation and propagation of floating point round‐off errors. This is demonstrated with the aid of the least square logarithm (LSL) scaling algorithm on FE consolidation problems of increasing size up to more than 100 000. It is shown that a major source of numerical instability rests with the sub‐matrix which couples the structural to the fluid part of the underlying mathematical model. It is concluded that for mid‐large size, possibly difficult, FE consolidation problems left/right LSL scaling is to be always recommended when the incomplete factorization is used as a preconditioning technique. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
3.
This paper identifies imbalanced columns (or rows) as a significant source of ill‐conditioning in the preconditioned coefficient matrix using the standard Jacobi preconditioner, for finite element solution of Biot's consolidation equations. A simple and heuristic preconditioner is proposed to reduce this source of ill‐conditioning. The proposed preconditioner modifies the standard Jacobi preconditioner by scaling the excess pore pressure degree‐of‐freedoms in the standard Jacobi preconditioner with appropriate factors. The performance of such preconditioner is examined using the symmetric quasi‐minimal residual method. To alleviate storage requirements, element‐by‐element iterative strategies are implemented. Numerical experiment results show that the proposed preconditioner reduces both the number of iteration and CPU execution time significantly as compared with the standard Jacobi preconditioner. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
4.
Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix Kep (ILU(0)‐Kep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐Kep. To further reduce computational cost, the global stiffness matrix Kep is divided into two parts. The first part is the linear elastic stiffness matrix Ke, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐Ke preconditioner is more time effective than IDR(s = 6) with ILU(0)‐Kep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
5.
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. 相似文献
6.
This paper presents a complete finite‐element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress–strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time‐stepping scheme which automatically adjusts the time‐step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non‐linearity caused by suction‐dependent plastic yielding, suction‐dependent degree of saturation, and saturation‐dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress–strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
7.
Giuseppe Gambolati Massimiliano Ferronato Pietro Teatini Roberto Deidda Giuditta Lecca 《国际地质力学数值与分析法杂志》2001,25(4):307-327
The solution of the poroelastic equations for predicting land subsidence above productive gas/oil fields may be addressed by the principle of virtual works using either the effective intergranular stress, with the pore pressure gradient regarded as a distributed body force, or the total stress incorporating the pore pressure. In the finite element (FE) method both approaches prove equivalent at the global assembled level. However, at the element level apparently the equivalence does not hold, and the strength source related to the pore pressure seems to generate different local forces on the element nodes. The two formulations are briefly reviewed and discussed for triangular and tetrahedral finite elements. They are shown to yield different results at the global level as well in a three‐dimensional axisymmetric porous medium if the FE integration is performed using the average element‐wise radius. A modification to both formulations is suggested which allows to correctly solve the problem of a finite reservoir with an infinite pressure gradient, i.e. with a pore pressure discontinuity on its boundary. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
8.
This paper examines the performance of the Jacobi preconditioner when used with two Krylov subspace iterative methods. The number of iterations needed for convergence was shown to be different for drained, undrained and consolidation problems, even for similar condition number. The differences were due to differences in the eigenvalue distribution, which cannot be completely described by the condition number alone. For drained problems involving large stiffness ratios between different material zones, ill‐conditioning is caused by these large stiffness ratios. Since Jacobi preconditioning operates on degrees‐of‐freedom, it effectively homogenizes the different spatial sub‐domains. The undrained problem, modelled as a nearly incompressible problem, is much more resistant to Jacobi preconditioning, because its ill‐conditioning arises from the large stiffness ratios between volumetric and distortional deformational modes, many of which involve the similar spatial domains or sub‐domains. The consolidation problem has two sets of degrees‐of‐freedom, namely displacement and pore pressure. Some of the eigenvalues are displacement dominated whereas others are excess pore pressure dominated. Jacobi preconditioning compresses the displacement‐dominated eigenvalues in a similar manner as the drained problem, but pore‐pressure‐dominated eigenvalues are often over‐scaled. Convergence can be accelerated if this over‐scaling is recognized and corrected for. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
9.
Finite element discretization of Biot's consolidation equations can produce a symmetric indefinite system (commonly used in geomechanics) or a non‐symmetric system. While this difference appears to be minor, however, it will require the selection of entirely different Krylov subspace solvers with potentially significant impact on solution efficiency. The former is solved using the symmetric quasi‐minimal residual whereas the latter is solved using the popular bi‐conjugate gradient stabilized. This paper presents an extensive comparison of the symmetric and non‐symmetric forms by varying the time step, size of the spatial domain, choice of physical units, and left versus left–right preconditioning. The generalized Jacobi (GJ) preconditioner is able to handle the non‐symmetric version of Biot's finite element method equation, although there are no practical incentives to do so. The convergence behaviour of GJ‐preconditioned systems and its relation to the spectral condition number or the complete spectrum are studied to clarify the concept of ill‐conditioning within the context of iteration solvers. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
10.
Daniel T. C. Yao Waldyr Lopes de Oliveira‐Filho Xiao‐Chuan Cai Dobroslav Znidarcic 《国际地质力学数值与分析法杂志》2002,26(2):139-161
The consolidation and desiccation behaviour of soft soils can be described by two time‐dependent non‐linear partial differential equations using the finite strain theory. Analytical solutions do not exist for these governing equations. In this paper, we develop efficient numerical methods and software for finding the numerical solutions. We introduce a semi‐implicit time integration scheme, and show numerically that our method converges. In addition, the numerical solution matches well with the experimental result. A boundary refinement method is also developed to improve the convergence and stability for the case of Neumann type boundary conditions. Interface governing equations are derived to maintain the continuity of consolidation and desiccation processes. This is useful because the soil column can undergo desiccation on top and consolidation on the bottom simultaneously. The numerical algorithms has been implemented into a computer program and the results have been verified with centrifuge test results conducted in our laboratory. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
11.
12.
The repeated solution in time of the linear system arising from the finite element integration of coupled consolidation equations is a major computational effort. This system can be written in either a symmetric or an unsymmetric form, thus calling for the implementation of different preconditioners and Krylov subspace solvers. The present paper aims at investigating when either a symmetric or an unsymmetric approach should be better used. The results from a number of representative numerical experiments indicate that a major role in selecting either form is played by the preconditioner rather than by the Krylov subspace method itself. Two other important issues addressed are the size of the time integration step and the possible lumping of the flow capacity matrix. It appears that ad hoc block constrained preconditioners provide the most robust algorithm independently of the time step size, lumping, and symmetry. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
13.
由于非饱和土的渗透系数是基质吸力的函数,使得控制方程带有强非线性的特征,进而使得控制方程的解析求解变得十分困难。同伦分析法对级数基函数和辅助线性算子的选择具有更大的自由性、灵活性,且收敛性的控制和调节更加容易实现,求解强非线性微分方程时在选择线性算子以及辅助参数上具有明显的优势。因此,针对非饱和土固结方程的非线性特征,对于处于地表浅层的非饱和土层,假设孔隙气压力为大气压力,在Richard经验公式与非饱和土一维固结理论的基础上,推导了非饱和一维固结无量纲控制方程;应用同伦分析法,通过选取适当的初始猜测解与辅助参数,将该非线性方程转换为线性的微分方程组并求解得到固结问题的级数解。此外,以压实高岭土为研究对象,在收集相关试验参数基础之上,将由同伦分析法求得的固结问题的近似解析解与有限差分法数值结果相对比,分析结果验证了解析解的正确性。 相似文献
14.
随着地球物理设备和探测技术的不断发展,快速处理大规模地球物理数据的需求也随之增长。为了解决三维重力数据密度反演的耗时问题,提出一种并行的预处理共轭梯度算法来提高计算效率。本文分别采用两种不同的预处理算子通过组合模型数据反演进行测试比较,并利用迭代残差和计算用时共同评价其加速效果。结果表明：对称逐次超松弛预处理方法比对角预处理方法反演计算速度快,密度结果更贴近实际模型;与传统串行的共轭梯度算法相比,本文并行预处理快速算法可以获得近19倍的加速比。将该算法应用于美国Vinton盐丘的实测重力数据中,反演结果能够很好地圈定出岩体的位置,验证了本文并行预处理共轭梯度法在三维重力数据快速反演中的高效性和可行性。 相似文献
15.
Parallel computers are potentially very attractive for the implementation of large size geomechanical models. One of the main difficulties of parallelization, however, relies on the efficient solution of the frequently ill‐conditioned algebraic system arising from the linearization of the discretized equilibrium equations. While very efficient preconditioners have been developed for sequential computers, not much work has been devoted to parallel solution algorithms in geomechanics. The present study investigates the state‐of‐the‐art performance of the factorized sparse approximate inverse (FSAI) as a preconditioner for the iterative solution of ill‐conditioned geomechanical problems. Pre‐and post‐filtration strategies are experimented with to increase the FSAI efficiency. Numerical results show that FSAI exhibits a promising potential for parallel geomechanical models mainly because of its almost ideal scalability. With the present formulation, however, at least 4 or 8 processors are required in the selected test cases to outperform one of the most efficient sequential algorithms available for FE geomechanics, i.e. the multilevel incomplete factorization (MIF). Further research is needed to improve the FSAI efficiency with a more effective selection of the preconditioner non‐zero pattern. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
16.
The mixed finite-element approximation to a second-order elliptic PDE results in a saddle-point problem and leads to an indefinite
linear system of equations. The mixed system of equations can be transformed into coupled symmetric positive-definite matrix
equations, or a Schur complement problem, using block Gauss elimination. A preconditioned conjugate-gradient algorithm is
used for solving the Schur complement problem. The mixed finite-element method is closely related to the cell-centered finite
difference scheme for solving second-order elliptic problems with variable coefficients. For the cell-centered finite difference
scheme, a simple multigrid algorithm can be defined and used as a preconditioner. For distorted grids, an additional iteration
is needed. Nested iteration with a multigrid preconditioned conjugate gradient inner iteration results in an effective numerical
solution technique for the mixed system of linear equations arising from a discretization on distorted grids. Numerical results
show that the preconditioned conjugate-gradient inner iteration is robust with respect to grid size and variability in the
hydraulic conductivity tensor. 相似文献
17.
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. 相似文献
18.
In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
19.
The problem of finite element simulation of incompressible fluid flow in porous medium is considered. The porous medium is characterized by the X‐ray microtomography technique in three dimensions. The finite calculus‐based stabilization technique is reviewed to implement the equal order finite element interpolation functions for both velocity and pressure. A noble preconditioner, the nodal block diagonal preconditioner, is considered whose performance is thoroughly investigated. Combining this preconditioner with a standard iterative solver during the computational homogenization procedure, it is possible to carry out the large‐scale fluid flow simulation for estimating permeability of the porous medium with reasonable accuracy and reliability. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
20.
In this study, dynamic behavior and earthquake resistance of Alibey earth dam was investigated. The dam was modeled with four
node plane-strain finite elements (FE) and displacement-pore pressure coupled FE analyses were performed. Nonlinear material
models such as pressure dependent and independent multi yield materials were implemented during the analyses. Transient dynamic
FE analyses were performed with Newmark method. The Newton-Raphson solution scheme was adopted to solve the equations. Liquefaction
and/or cyclic mobility effects were considered during the analysis. For the FE analyses, OpenSees (Open System for Earthquake
Engineering Simulation) framework was adopted. 相似文献