A uniform multiscale method for 2D static and dynamic analyses of heterogeneous materials

A uniform extended multiscale finite element method is developed for solving the static and dynamic problems of heterogeneous materials in elasticity. To describe the complex deformation, a multinode two‐dimensional coarse element is proposed, and a new approach is elaborated to construct the displacement base functions of the coarse element. In addition, to improve the computational accuracy, the mode base functions are introduced to consider the effect of the inertial forces of the structure for dynamic problems. Furthermore, the orthogonality between the displacement and mode base functions is proved theoretically, which indicates that the proposed multiscale method can be used for the static and dynamic analyses uniformly. Numerical experiments show that the mode base functions almost do not work for the static problems, while they can improve the computational accuracy of the dynamic problems significantly. On the other hand, it is also found that the number of the macro nodes of the multinode coarse element has a great influence on the accuracy of the numerical results for both the static and dynamic analyses. Numerical examples also indicate that the uniform extended multiscale finite element method can obtain sufficiently accurate results with less computational cost compared with the standard FEM. Copyright © 2012 John Wiley & Sons, Ltd.     

A fully coupled finite element analysis of water‐table fluctuation and land deformation in partially saturated soils due to surface loading

A fully coupled numerical model is presented for the water‐table fluctuation and land deformation in partially saturated soils due to surface loading. This numerical model is developed based on the poroelastic governing equations for groundwater flow in deforming variably saturated porous media and the Galerkin finite element method. The numerical model is verified and validated against a one‐dimensional consolidation problem concerning surface loading on a soil column which has six different initial water‐table elevations. The numerical model is then applied to a two‐dimensional consolidation problem of surface loading on a partially saturated soil at a construction site. Results from the numerical simulations of both problems show that the water table fluctuates in the partially saturated soils, and the unsaturated zone above the water table has significant effects on the consolidation behaviour of the partially saturated soils under surface loading. Such effects are caused by the permanent absorption of a portion of the mechanical loading stress and the weak hydromechanical coupling between the solid skeleton deformation field and the groundwater flow field in the unsaturated zone due to its partial saturation. Copyright © 2000 John Wiley & Sons, Ltd.     

基于耦合扩展多尺度有限元方法的功能梯度材料热应力分析

以高效模拟功能梯度材料(FGM)微观非均质性对整体热力学性能的影响为研究目的,通过随机形态描述函数(RMDF)法和体积分数的指数分布建立FGM二维微结构,在此基础上,发展了FGM热应力分析的耦合扩展多尺度有限元方法(CEMsFEM)。该方法基于扩展多尺度有限元方法(EMsFEM)的基本思想,对温度场和位移场构造数值基函数,以把微观非均质材料性质带到宏观响应中。同时为了考虑泊松效应导致的不同方向间的耦合作用,在位移场数值基函数中增加了耦合附加项。通过数值基函数建立宏微观单元信息的映射关系,在宏观尺度求解有效方程,节约计算量。为了更好地考虑微观载荷的影响,把结构的真实响应分解为宏观响应和微观扰动,进一步推导出修正的宏观载荷向量。通过不同体积分数分布的FGM在不同载荷工况下的热应力分析算例验证了本文中方法的正确性和有效性,最后讨论了微结构的尺寸效应对结构热力学响应的影响。     

Fully coupled hydromechanical multiscale model with microdynamic effects

In this paper, a multiscale finite element framework is developed based on the first‐order homogenization method for fully coupled saturated porous media using an extension of the Hill‐Mandel theory in the presence of microdynamic effects. The multiscale method is employed for the consolidation problem of a 2‐dimensional saturated soil medium generated from the periodic arrangement of circular particles embedded in a square matrix, which is compared with the direct numerical simulation method. The effects of various issues, including the boundary conditions, size effects, particle arrangements, and the integral domain constraints for the microscale boundary value problem, are numerically investigated to illustrate the performance of a representative volume element in the proposed computational homogenization method of fully coupled saturated porous media. This study is aimed to clarify the effect of scale separation and size dependence, and to introduce characteristics of a proper representative volume element in multiscale modeling of saturated porous media.     

An efficient finite element procedure for analyzing three‐phase porous media based on the relaxed Picard method

Effective simulation of the solid‐liquid‐gas coupling effect in unsaturated porous media is of great significance in many diverse areas. Because of the strongly nonlinear characteristics of the fully coupled formulations for the three‐phase porous media, an effective numerical solution scheme, such as the finite element method with an efficient iterative algorithm, has to be employed. In this paper, an efficient finite element procedure based on the adaptive relaxed Picard method is developed for analyzing the coupled solid‐liquid‐gas interactions in porous media. The coupled model and the finite element analysis procedure are implemented into a computer code PorousH2M, and the proposed procedure is validated through comparing the numerical simulations with the experimental benchmarks. It is shown that the adaptive relaxed Picard method has salient advantage over the traditional one with respect to both the efficiency and the robustness, especially for the case of relatively large time step sizes. Compared with the Newton‐Raphson scheme, the Picard method successfully avoids the unphysical ‘spurious unloading’ phenomenon under the plastic deformation condition, although the latter shows a better convergence rate. The proposed procedure provides an important reference for analyzing the fully coupled problems related to the multi‐phase, multi‐field coupling in porous media. Copyright © 2014 John Wiley & Sons, Ltd.     

Finite–infinite element for dynamic analysis of axisymmetrically saturated composite foundations

The formulations of axisymmetrically infinite elements for dynamic analysis of vertical vibration problems in unbounded saturated composite foundations are presented. The theoretical basis as well as the implementation of the elements is discussed, and the element decay functions are derived using the analytical solutions of axially symmetric configurations. Using the proposed finite–infinite element method, the surface vertical displacements of air‐saturated soil (‘dry’ soil) and of water‐saturated soil with extremely low permeability subjected to a surface point excitation (called as the Lamb's problem) are calculated and the results agree very well with the existing theoretical solutions of single‐phase elastic media. As an application, the velocity admittances of a concrete block resting on cement mixing‐pile or gravel‐pile saturated composite foundations are calculated. The influence of soil permeability and pile rigidity on the dynamic response of the composite foundations is investigated. The method proposed by this paper is a simple and reliable numerical one that could be used to study axisymmetrically dynamic problems of layered saturated media and to get the mechanism of dynamic testing on single‐pile saturated composite foundations. Copyright © 2006 John Wiley & Sons, Ltd.     

Computational simulation of chemical dissolution‐front instability in fluid‐saturated porous media under non‐isothermal conditions

This paper primarily deals with the computational aspects of chemical dissolution‐front instability problems in two‐dimensional fluid‐saturated porous media under non‐isothermal conditions. After the dimensionless governing partial differential equations of the non‐isothermal chemical dissolution‐front instability problem are briefly described, the formulation of a computational procedure, which contains a combination of using the finite difference and finite element method, is derived for simulating the morphological evolution of chemical dissolution fronts in the non‐isothermal chemical dissolution system within two‐dimensional fluid‐saturated porous media. To ensure the correctness and accuracy of the numerical solutions, the proposed computational procedure is verified through comparing the numerical solutions with the analytical solutions for a benchmark problem. As an application example, the verified computational procedure is then used to simulate the morphological evolution of chemical dissolution fronts in the supercritical non‐isothermal chemical dissolution system. The related numerical results have demonstrated the following: (1) the proposed computational procedure can produce accurate numerical solutions for the planar chemical dissolution‐front propagation problem in the non‐isothermal chemical dissolution system consisting of a fluid‐saturated porous medium; (2) the Zhao number has a significant effect not only on the dimensionless propagation speed of the chemical dissolution front but also on the distribution patterns of the dimensionless temperature, dimensionless pore‐fluid pressure, and dimensionless chemical‐species concentration in a non‐isothermal chemical dissolution system; (3) once the finger penetrates the whole computational domain, the dimensionless pore‐fluid pressure decreases drastically in the non‐isothermal chemical dissolution system.     

Stabilized low‐order finite elements for strongly coupled poromechanical problems

Finite element solutions of poromechanical problems often exhibit oscillating pore pressures in the limits of low permeability, fast loading rates, coarse meshes, and/or small time step sizes. To suppress completely the pore pressure oscillations, a stabilized finite element scheme with a better performance on monotonicity is proposed for modeling compressible fluid‐saturated porous media. This method, based on the polynomial pressure projection technique, allows the use of linear equal‐order interpolation for both displacement and pore pressure fields, which is more straightforward for both code development and maintenance compared to others. By employing the discrete maximum principle, a proper stabilization parameter is deduced, which is efficient to guarantee the monotonicity and optimal in theory in the 1‐dimensional case. An appealing feature of the method is that the stabilization parameter is evaluated in terms of the properties of porous material only, while no mesh or time step size is involved. Through comparing the numerical simulations with the analytical benchmarks, the efficiency of the proposed stabilization scheme is confirmed.     

考虑土-结构相互作用大型动力机器基础三维有限元分析

本文提出了考虑土-结构相互作用计算大型机器基础动力反应的三维有限元方法.该方法将动力机器基础与周围地基作为一个完整的研究对象,将基础作为粘弹性体反映基础变形,建立粘-弹性人工边界模拟半无限地基的辐射阻尼和远场地球介质的弹性恢复性能,并通过构造等效弹簧实体元使这一人工边界在商业有限元程序中实现。文中进一步将上述方法用于分析燕化高压聚乙稀装置压缩机基础的动力反应，并通过数值试验定量讨论了影响大型机器基础动力反应的主要因素，得出了一些有用的结论。     

A dual‐primal FETI method for solving a class of fluid–structure interaction problems in the frequency domain

The dual‐primal finite element tearing and interconnecting method (FETI‐DP) is extended to systems of linear equations arising from a finite element discretization for a class of fluid–structure interaction problems in the frequency domain. A preconditioned generalized minimal residual method is used to solve the linear equations for the Lagrange multipliers introduced on the subdomain boundaries to enforce continuity of the solution. The coupling between the fluid and the structure on the fluid–structure interface requires an appropriate choice of coarse level degrees of freedom in the FETI‐DP algorithm to achieve fast convergence. Several choices are proposed and tested by numerical experiments on three‐dimensional fluid–structure interaction problems in the mid‐frequency regime that demonstrate the greatly improved performance of the proposed algorithm over the standard FETI‐DP method. Copyright © 2011 John Wiley & Sons, Ltd.     

Mixed finite element method for saturated poroelastoplastic media at large strains

A mixed finite element for hydro‐dynamic analysis in saturated porous media in the frame of the Biot theory is proposed. Displacements, effective stresses, strains for the solid phase and pressure, pressure gradients, and Darcy velocities for the fluid phase are interpolated as independent variables. The weak form of the governing equations of coupled hydro‐dynamic problems in saturated porous media within the element are given on the basis of the Hu–Washizu three‐field variational principle. In light of the stabilized one point quadrature super‐convergent element developed in solid continuum, the interpolation approximation modes for the primary unknowns and their spatial derivatives of the solid and the fluid phases within the element are assumed independently. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of pressure‐dependent non‐associated plasticity. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elastoplastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is used. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization due to strain softening in poroelastoplastic media subjected to dynamic loading at large strain. Copyright © 2003 John Wiley & Sons, Ltd.     

Multiscale analysis method for thermo‐mechanical performance of periodic porous materials with interior surface radiation

This study develops a novel multiscale analysis method to predict thermo‐mechanical performance of periodic porous materials with interior surface radiation. In these materials, thermal radiation effect at microscale has an important impact on the macroscopic temperature and stress field, which is our particular interest in this paper. Firstly, the multiscale asymptotic expansions for computing the dynamic thermo‐mechanical coupling problem, which considers the mutual interaction between temperature and displacement field, are given successively. Then, the corresponding numerical algorithm based on the finite element‐difference method is brought forward in details. Finally, some numerical results are presented to verify the validity and relevancy of the proposed method by comparing it with a direct finite element analysis with detailed numerical models. The comparison shows that the new method is effective and valid for predicting the thermo‐mechanical performance and can capture the microstructure behavior of periodic porous materials exactly.s Copyright © 2015 John Wiley & Sons, Ltd.     

A two‐scale approach for fluid flow in fractured porous media

A two‐scale numerical model is developed for fluid flow in fractured, deforming porous media. At the microscale the flow in the cavity of a fracture is modelled as a viscous fluid. From the micromechanics of the flow in the cavity, coupling equations are derived for the momentum and the mass couplings to the equations for a fluid‐saturated porous medium, which are assumed to hold on the macroscopic scale. The finite element equations are derived for this two‐scale approach and integrated over time. By exploiting the partition‐of‐unity property of the finite element shape functions, the position and direction of the fractures is independent from the underlying discretization. The resulting discrete equations are non‐linear due to the non‐linearity of the coupling terms. A consistent linearization is given for use within a Newton–Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach, and show that faults in a deforming porous medium can have a significant effect on the local as well as on the overall flow and deformation patterns. Copyright © 2006 John Wiley & Sons, Ltd.     

A coupled molecular dynamics and extended finite element method for dynamic crack propagation

A multiscale method is presented which couples a molecular dynamics approach for describing fracture at the crack tip with an extended finite element method for discretizing the remainder of the domain. After recalling the basic equations of molecular dynamics and continuum mechanics, the discretization is discussed for the continuum subdomain where the partition‐of‐unity property of finite element shape functions is used, since in this fashion the crack in the wake of its tip is naturally modelled as a traction‐free discontinuity. Next, the zonal coupling method between the atomistic and continuum models is recapitulated. Finally, it is discussed how the stress has been computed in the atomic subdomain, and a two‐dimensional computation is presented of dynamic fracture using the coupled model. The result shows multiple branching, which is reminiscent of recent results from simulations on dynamic fracture using cohesive‐zone models. Copyright © 2009 John Wiley & Sons, Ltd.     

饱和多孔介质动力分析的数值流形单元

基于数值流形方法中覆盖函数的基本思想,构造了适用于饱和多孔介质动力耦合分析的三节点平面流形单元,该单元满足Babuska-Brezzi稳定性准则与Zienkiewicz-Taylor分片试验条件,对于位移和孔隙压力具有不等阶的插值函数,且所有节点上具有相同自由度。用标准Galerkin法和Newmark法将饱和多孔介质动力基本方程在空间和时间上离散,得到饱和多孔介质动力分析的流形元离散的算法公式。数值结果表明,与传统有限元相比在孔隙流体不可压缩且非渗流的条件下,数值流形单元对于压力场的计算具有良好的数值稳定性。     

Lower‐dimensional interface elements with local enrichment: application to coupled hydro‐mechanical problems in discretely fractured porous media

In this study, we develop lower‐dimensional interface elements to represent preexisting fractures in rock material, focusing on finite element analysis of coupled hydro‐mechanical problems in discrete fractures–porous media systems. The method adopts local enrichment approximations for a discontinuous displacement and a fracture relative displacement function. Multiple and intersected fractures can be treated with the new scheme. Moreover, the method requires less mesh dependencies for accurate finiteelement approximations compared with the conventional interface element method. In particular, for coupled problems, the method allows for the use of a single mesh for both mechanical and other related processes such as flow and transport. For verification purposes, several numerical examples are examined in detail. Application to a coupled hydro‐mechanical problem is demonstrated with fluid injection into a single fracture. The numerical examples prove that the proposed method produces results in strong agreement with reference solutions. Copyright © 2012 John Wiley & Sons, Ltd.     

On a two‐level finite element method for the incompressible Navier–Stokes equations

We consider the Galerkin finite element method for the incompressible Navier–Stokes equations in two dimensions, where the finite‐dimensional space(s) employed consist of piecewise polynomials enriched with residual‐free bubble functions. To find the bubble part of the solution, a two‐level finite element method (TLFEM) is described and its application to the Navier–Stokes equation is displayed. Numerical solutions employing the TLFEM are presented for three benchmark problems. We compare the numerical solutions using the TLFEM with the numerical solutions using a stabilized method. Copyright © 2001 John Wiley & Sons, Ltd.     

斜入射下水下地基场地地震输入的一维化时域方法

地震波斜入射下水下地基场地地震动输入问题,是涉及水层、饱和多孔介质和基岩三种不同性质介质耦合的复杂问题,目前理论成果还很匮乏。该文在现有的成层饱和多孔介质平面波斜入射的一维化时域方法的基础上,建立了水下地基场地的一维化时域计算方法。该计算方法依据Snell定理,将地震波斜入射下水层波动的空间二维问题转化为简单的一维问题,通过考虑水层与饱和多孔介质层交界面、不同饱和多孔介质交界面以及基岩与饱和多孔介质层交界面的边界条件,将已有的成层饱和土的一维化有限元方程与水层的一维化有限元方程组装,采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,形成了水下地基场地的整体有限元方程,借助于中心差分法和Newmark法相结合的时步积分法,推导了地震波斜入射下水下地基场地各节点动力响应时程的显式表达式。通过与傅里叶变换得到的理论解和现有的文献进行了对比,初步验证了该文时域算法的有效性和精度,为地震波斜入射下水下地基场地地震动输入提供了一种新型高效的研究方法。     

Weak formulation of Biot's equations in cylindrical coordinates with harmonic expansion in the circumferential direction

A weak symmetric form of Biot's equation in cylindrical coordinates with a spatial Fourier expansion in the circumferential direction is presented. The solid phase displacement and the pore pressure are used as the dependent variables. The original three‐dimensional boundary value problem is here, due to the orthogonality of the harmonic functions and the rotationally symmetric geometry, decomposed into independent two‐dimensional problems, one for each harmonic function. This formulation provides a computationally efficient procedure for vibroacoustic finite element modelling of rotationally symmetric three‐dimensional multilayered structures including porous elastic materials. By numerical simulations, this method is compared with, and verified against, full three‐dimensional Cartesian coordinate system finite element models. Copyright © 2009 John Wiley & Sons, Ltd.     

A finite element solution of acoustic propagation in rigid porous media

This paper deals with the acoustical behaviour of a rigid porous material. A finite element method to compute both the response to an harmonic excitation and the free vibrations of a three‐dimensional finite multilayer system consisting of a free fluid and a rigid porous material is considered. The finite element used is the lowest order face element introduced by Raviart and Thomas, that eliminates the spurious or circulation modes with no physical meaning. For the porous medium a Darcy's like model and the Allard–Champoux model are taken into account. The numerical results show that the finite element method allows us to compute the response curve for the coupled system and the complex eigenfrequencies. Some of them have a small imaginary part but there are also overdamped modes. Copyright © 2004 John Wiley & Sons, Ltd.