共查询到20条相似文献，搜索用时 672 毫秒
1.
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. 相似文献
2.
Qingli Dai 《国际地质力学数值与分析法杂志》2010,34(13):1321-1345
This paper presents parallel and serial viscoelasto‐plastic models to simulate the rate‐independent and the rate‐dependent permanent deformation of stone‐based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto‐plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement‐based incremental formulation for the viscoelastic part and an elastic predictor—plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto‐plastic model employed a viscoelastic predictor—plastic corrector algorithm. The stone‐based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto‐plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto‐plastic matrix models. The simulation results showed that the serial viscoelasto‐plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto‐plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate‐independent and rate‐dependent viscoelasto‐plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto‐plastic behaviors of the stone‐based materials. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
3.
研究了不同测量方式(井-地,地-井,井-井)下点源场井中电法的三维有限元数值模拟。考虑到深度方向上大范围的网格剖分和井眼的影响及井-井测量等因素,采用放射状三棱柱单元的网格剖分方式,以提高网格质量,减少剖分单元数;给出了三棱柱单元的坐标变换公式,进行精确的单元积分,减少了单元积分时间;结合非结构化网格技术,实现了复杂模型的模拟;开发出相应的程序实现了复杂条件下(如考虑井眼影响、井井测量、倾斜井情形、地形起伏等)电法测井的三维有限元模拟,数值算例验证了方法的可靠性及计算效率,并对不同情形下的异常响应进行了分析,为进一步的反演工作奠定了基础。 相似文献
4.
Konstantinos I. Andrianopoulos Achilleas G. Papadimitriou George D. Bouckovalas 《国际地质力学数值与分析法杂志》2010,34(15):1586-1614
This paper presents a new plasticity model developed for the simulation of monotonic and cyclic loading of non‐cohesive soils and its implementation to the commercial finite‐difference code FLAC, using its User‐Defined‐Model (UDM) capability. The new model incorporates the framework of Critical State Soil Mechanics, while it relies upon bounding surface plasticity with a vanished elastic region to simulate the non‐linear soil response. Stress integration of constitutive relations is performed using a recently proposed explicit scheme with automatic error control and substepping, which so far has been employed in the literature only for constitutive models aiming at monotonic loading. The overall accuracy of this scheme is evaluated at element level by simulating cyclic loading along complex stress paths and by using iso‐error maps for paths involving change of the Lode angle. The performance of the new constitutive model and its stress integration scheme in complex boundary value problems involving earthquake‐induced liquefaction is evaluated, in terms of accuracy and computational cost, via a number of parametric analyses inspired by the successful simulation of the VELACS centrifuge Model Test No. 2 studying the lateral spreading response of a liquefied sand layer. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
5.
This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
6.
Filip Nilenius Fredrik Larsson Karin Lundgren Kenneth Runesson 《国际地质力学数值与分析法杂志》2013,37(11):1535-1551
Effective moisture and chloride ion diffusivity coefficients for concrete are determined by computational homogenization, where concrete is modeled on the mesoscale as a heterogenous three‐phase composite material. By imposing moisture and chloride ion gradients on a representative volume element, effective macroscale properties are obtained through finite element analysis. A parametric study of the effects of the ballast content was carried out. The numerical results correspond well with an estimate of the Hashin–Shtrikman type available in the literature. The computational homogenization strategy proposed here also includes the interfacial transition zone, and its influence on the effective diffusivity coefficients is assessed. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
7.
8.
A simplified computational technique based on a refined global–local method is applied to the failure analysis of concrete structures. The technique distinguishes the scale of the structure, modelled with large size finite elements, from the scale at which material non‐linearity occurs due to progressive cracking and macro‐crack propagation. The finite element solution is split into two parts: a linear elastic analysis on a coarse mesh over the entire structure and a non‐linear analysis over a small part of the structure where a dense finite element grid is employed. In the non‐linear calculation, a non‐local damage model is implemented. These two computations are coupled with the help of an iterative scheme. The size and location of the region where a non‐linear analysis is performed, is adapted to follow the development of the damage zone. Numerical examples of mode I fracture of concrete specimens with straight and curved cracks are presented. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
9.
This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
10.
可液化土中地铁结构的地震响应 总被引：14，自引：3，他引：11
在饱和土耦合作用与土和结构相互作用理论基础上,以地铁车站为例,用有限元法研究地下结构在地震液化作用下的响应。所采用的软件为动力两相体非线性有限元软件Dyna-Swandyne-II,该软件可以应用先进的Pastor-Zienkiewicz III广义塑性模型模拟可液化土的动力特性,应用u-p形式的Biot方程,在有限元分析中充分考虑孔隙水与土之间的耦合,同时考虑地下结构与饱和土在动力作用下的非线性相互作用。分析了地铁车站的动力响应,包括地铁内力、加速度以及地铁位移。研究结果表明,地铁结构在地震液化作用下会产生较大的上浮,从而对结构造成比较严重的破坏;地铁结构在地震作用下的最大内力位于结构的交接处。因此,结构交接处的配筋应该格外小心。 相似文献
11.
The material point method (MPM), which is a combination of the finite element and meshfree methods, suffers from significant computational workload due to the fine mesh that is required in spite of its advantages in simulating large deformations. This paper presents a parallel computing strategy for the MPM on the graphics processing unit (GPU) to boost the method’s computational efficiency. The interaction between a structural element and soil is investigated to validate the applicability of the parallelisation strategy. Two techniques are developed to parallelise the interpolation from soil particles to nodes to avoid a data race; the technique that is based on workload parallelisation across threads over the nodes has a higher computational efficiency. Benchmark problems of surface footing penetration and a submarine landslide are analysed to quantify the speedup of GPU parallel computing over sequential simulations on the central processing unit. The maximum speedup with the GPU used is ∼30 for single-precision calculations and decreases to ∼20 for double-precision calculations. 相似文献
12.
Joachim Berdal Haga Harald Osnes Hans Petter Langtangen 《Computational Geosciences》2012,16(3):723-734
Large-scale simulations of coupled flow in deformable porous media require iterative methods for solving the 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 realistic geological applications, such as basin evolution at regional
scales. The success of iterative methods for such problems depends strongly on finding effective preconditioners with good
parallel scaling properties, which is the topic of the present paper. We present a parallel preconditioner for Biot’s equations
of coupled elasticity and fluid flow in porous media. The preconditioner is based on an approximation of the exact inverse
of the two-by-two block system arising from a finite element discretisation. The approximation relies on a highly scalable
approximation of the global Schur complement of the coefficient matrix, combined with generally available state-of-the-art
multilevel preconditioners for the individual blocks. This preconditioner is shown to be robust on problems with highly heterogeneous
material parameters. We investigate the weak and strong parallel scaling of this preconditioner on up to 512 processors and
demonstrate its ability on a realistic basin-scale problem in poroelasticity with over eight million tetrahedral elements. 相似文献
13.
基于有限元理论,进行了缓坡土层地震液化引起大变形的数值方法研究,即采用二维有效应力动力有限元方法进行分析,在液化分析过程的每一时段考虑地震液化和振动软化得到土单元的模量,通过非线性静力方法计算每时段地震液化引起的大变形,得到土层各深度处的水平和竖向位移。由算例分析了地震动和土层坡度等因素的影响,通过对比分析表明了该方法的有效性,可为工程场地地震地质灾害评价提供参考数据。 相似文献
14.
弹簧元法是一种将单元离散为一系列弹簧的数值计算方法。不同的单元具有不同的离散方式,确定相应的离散弹簧的刚度系数表达式是弹簧元法的关键。将四节点矩形单元离散为6个基本弹簧,每个基本弹簧包括法向弹簧和切向弹簧两个派生弹簧,并用泊松弹簧和纯剪弹簧描述单元的泊松效应和剪切效应,用有限元的单元刚度矩阵标定各弹簧的刚度系数,实现了一种四节点矩形弹簧元的构造形式。该单元的同类弹簧具有相同的表达形式。法向与切向弹簧的刚度表达式中分别含有法向和切向弹簧刚度待定系数。通过改变待定系数的值可使该单元分别对应于有限元的常应变、双线性及Wilson非协调单元。将上述弹簧元方法进行理论推导,并应用于基于连续介质的离散单元法（CDEM）的核心计算进行简单算例验证,证明了提出方法的正确性。通过以上研究发现,四节点矩形弹簧元有以下特点：对于相同问题,不同单元有不同的计算精度;对于梁弯曲问题,应用该单元可显著提高离散单元法的求解精度;改变待定系数的值,可得到更高或者更低精度的单元。 相似文献
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.
A four‐node plane parametric element based on quadrilateral area coordinate and its application to coupled solid‐deformation/fluid‐flow simulation for porous geomaterials 
下载免费PDF全文
下载免费PDF全文 Gen Li 《国际地质力学数值与分析法杂志》2015,39(3):251-276
A four‐node plane parametric element AQGβ6‐I is constructed on the basis of the quadrilateral area coordinate, the generalized conforming principle and the projection technique with a penalty factor β within an interval of 0–1. When β = 0, the element has excellent bending performance. When β = 1, the element can pass patch test strictly; its performance is as good as many famous elements. When β value is between 0 and 1, such as β = 0.5, the element can arrive at a compromise between (relatively) low sensitivity to mesh distortion and perfect convergence. The work provides an illuminating method to alleviate a difficult problem in finite element modelling using the four‐node quadrilateral element, which can pass the strict patch test, but has poor performance in bending dominated problem; on the contrary, it has excellent performance in bending dominated problem but cannot pass the strong patch test. The AQGβ6‐I with the convergence formulation (β = 1) is then applied to coupled solid‐deformation/fluid‐flow simulation for porous geomaterials. The computational examples are carried out to demonstrate that the AQGβ6‐I (β = 1) element is not only stable, reliable and efficient but also of high accuracy. The present study provides a good applicable element for finite element simulations of solid‐deformation/fluid‐flow for porous geomaterials. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
17.
In modeling of many geomechanics problems such as underground openings, soil-foundation structure interaction problems, and in wave propagation problems through semi-infinite soil medium the soil is represented as a region of either infinite or semi-infinite extent. Numerical modeling of such problems using conventional finite elements involves a truncation of the far field in which the infinite boundary is terminated at a finite distance. In these problems, appropriate boundary conditions are introduced to approximate the solution of the infinite or semi-infinite boundaries as closely as possible. However, the task of positioning the finite boundary in conventional finite element discretization and the definition of the boundary and its conditions is very delicate and depends on the modeller's skill and intuition. Moreover, such a choice is influenced by the size of the domain to be discretized. Consequently, the dimensions of the global matrices and the time required for solution of the problem will increase considerably and also selection of the arbitrary location of truncated boundary may lead to erroneous result. In order to over come these problems, mapped infinite elements have been developed by earlier researchers (Simoni and Schrefier, 1987). In the present work the applicability of infinite element technique is examined for different geomechanics problems. A computer program INFEMEP is developed based on the conventional finite element and mapped infinite element technique. It is then validated using selected problems such as strip footing and circular footing. CPU time taken to obtain solutions using finite element approach and infinite element approach was estimated and presented to show the capability of coupled modeling in improving the computational efficiency. Mesh configurations of different sizes were used to explore the enhancement of both computational economy and solution accuracy achieved by incorporation of infinite elements to solve elastic and elasto-plastic problems in semi-infinite/finite domain as applied to geotechnical engineering. © Rapid Science Ltd. 1998 相似文献
18.
In an urban regional groundwater flow field, the presence of thousands of channels between impervious structures makes the field difficult to simulate using the finite-element method (FEM), because the scale of the field is usually several orders of magnitude larger than that of the channels. To overcome this problem, a simple element for the simulation of potential channel flows has been developed. This element works with linear triangular elements and can be easily implemented in a finite-element code to simulate the channel flows with a sparse mesh without a loss of global accuracy. The transmissivity matrix of this element is deduced from the analytical solution for channel flow. The application of the element is discussed, and the accuracy of the element is assessed. The element makes it easy to merge small structures that are close to each other into a larger one that can be modeled with a sparse mesh. 相似文献
19.
General coupling extended multiscale FEM for elasto‐plastic consolidation analysis of heterogeneous saturated porous media 下载免费PDF全文
下载免费PDF全文 This paper presents a general coupling extended multiscale FEM (GCEMs) for solving the coupling problem of elasto‐plastic consolidation of heterogeneous saturated porous media. In the GCEMs, the numerical multiscale base functions for the solid skeleton and fluid phase of the coupling system are all constructed on the basis of the equivalent stiffness matrix of the unit cell, which not only contain the interaction between the solid and fluid phases but also consider the time effect. Furthermore, in order to improve the computational accuracy for two‐dimensional problems, a multi‐node coarse element strategy for the GCEMs is proposed, and a two‐scale iteration algorithm for the elasto‐plastic consolidation analysis is developed. Some one‐dimensional and two‐dimensional homogeneous and heterogeneous numerical examples are carried out to validate the proposed method through the comparison with the coupling multiscale FEM and standard FEM. Numerical results show that the newly developed GCEMs can almost preserve the same convergent property as the standard FEM and also possesses the advantages of high computational efficiency. In addition, the GCEMs can be easily applied to other coupling multifield and multiphase transient problems. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
20.
A computational method, based on an advanced elasto-plastic large strain formulation, well suited for the analysis of the cone penetration problem, is presented. A new approach of finite strain elasto-plastic analysis is employed.13 The basic (non-rate) constitutive relations are developed in a spatial reference space to preserve their physical significance. They are subsequently transformed in Lagrangian co-ordinates, and through simple time differentiation, their rate equations are introduced. The method is computationally implemented with the finite element method and special provisions are taken to allow for the moving boundary conditions of the problem. 相似文献