共查询到19条相似文献,搜索用时 179 毫秒
1.
变形多孔介质流固耦合模型及数值模拟研究 总被引:2,自引:0,他引:2
假定骨架、固体颗粒和水均是可压缩的,在此基础上采用两相不混溶流体的理论推导了水的连续方程,通过引入气压恒定这一假定进一步简化为水的非饱和渗流连续方程。基于广义Biot理论给出了固体骨架积分形式的平衡方程,结合非饱和渗流连续方程采用加权残值法推导了流固耦合方程组的有限元列式。通过干燥介质吸水的数值模拟来考察非饱和流固耦合模型的预测能力,数值模拟的结果表明耦合模型可以准确地反映吸水过程的规律。将耦合模型应用于水下大断面隧洞开挖的瞬态分析,可以模拟出开挖引起的EDZ区域孔隙水压力急剧升高、有效应力减小、渗透系数动态变化以及排水对洞室稳定性的影响,计算的结果与国外大型原位实验的一般性观测结论相吻合。 相似文献
2.
在已有研究工作基础上对非饱和多孔介质应变局部化问题进行研究,给出非饱和多孔介质的分析控制方程,其中饱和度与毛细压力关系由实验给出。采取适用于非饱和砂土的改进的广义塑性本构模型对应局部变化过程进行数值模拟,给出了试件应变局部化发展过程以及孔隙压力的变化规律。对初始饱和土中所产生非饱和剪切带进行计算的结果表明,采用非饱和模型较饱和模型将获得更为合理的结果。 相似文献
3.
在基于混合物理论的多孔介质模型的基础上,将固体相视为弹粘塑性体,建立了饱和多孔介质的弹粘塑性模型。模型的基本思想是在无粘弹塑性本构关系中引入-时间参数,使固体骨架具备了粘性效应。利用Galerkin加权残值法推导得到了罚有限元格式,并采用Newmark预估校正法求解率相关饱和多孔介质的非线性有限元动力方程,此算法可以很... 相似文献
4.
基于数值流形方法中覆盖函数的基本思想,构造了适用于饱和多孔介质动力耦合分析的三节点平面流形单元,该单元满足Babuska-Brezzi稳定性准则与Zienkiewicz-Taylor分片试验条件,对于位移和孔隙压力具有不等阶的插值函数,且所有节点上具有相同自由度。用标准Galerkin法和Newmark法将饱和多孔介质动力基本方程在空间和时间上离散,得到饱和多孔介质动力分析的流形元离散的算法公式。数值结果表明,与传统有限元相比在孔隙流体不可压缩且非渗流的条件下,数值流形单元对于压力场的计算具有良好的数值稳定性。 相似文献
5.
Biot饱和多孔介质波动行为的数值模拟在众多工程领域中具有重要的意义和作用,由于固相与液相耦合方程难以解耦,使该问题的数值模拟难度较大。针对饱和多孔介质中部分耦合u-p及全耦合u-p-U方程形式的特征,推导了相应动力耦合控制方程的有限元弱形式,并引入不同耦合形式的饱和多孔介质时域黏性边界,综合利用Comsol Multiphysics提供的偏微分方程应用模式进行二次开发求解,通过一维饱和多孔介质动力响应的解析解和数值解验证了模型求解技术的合理性和可行性,基于u-p-U耦合形式探讨了冲击荷载作用下干砂饱和砂地基动力固结中应力波传播特性。计算结果表明慢纵波对动力固结的影响比较显著,合理的冲击荷载持续时间有利于固结效果的改善。 相似文献
6.
《振动与冲击》2017,(9)
Biot饱和多孔介质波动行为的数值模拟在众多工程领域中具有重要的意义和作用,由于固相与液相耦合方程难以解耦,使该问题的数值模拟难度较大。针对饱和多孔介质中部分耦合u-p及全耦合u-p-U方程形式的特征,推导了相应动力耦合控制方程的有限元弱形式,并引入不同耦合形式的饱和多孔介质时域黏性边界,综合利用Comsol Multiphysics提供的偏微分方程应用模式进行二次开发求解,通过一维饱和多孔介质动力响应的解析解和数值解验证了模型求解技术的合理性和可行性,基于u-p-U耦合形式探讨了冲击荷载作用下干砂饱和砂地基动力固结中应力波传播特性。计算结果表明慢纵波对动力固结的影响比较显著,合理的冲击荷载持续时间有利于固结效果的改善。 相似文献
7.
8.
各向异性双重孔隙介质的应力与油水两相渗流耦合理论模型 总被引:1,自引:0,他引:1
针对天然裂缝性油藏的特性,建立了描述双重孔隙介质中油水两相流体流动特性的流固耦合理论模型。该模型不仅考虑了渗透率的各向异性,而且考虑了岩石固体骨架变形的各向异性。渗流方程是依据双重孔隙的概念建立起来的,而固体骨架变形控制方程则是根据Biot 的等温、线性孔隙弹性理论建立起来的。同时,给出了横向各向同性及结构各向异性、固体材料各向同性时的双重孔隙介质的应力与油水两相渗流耦合理论模型。对该模型进行了简化,并将其简化后模型与单相流的各项同性和各向异性双重孔隙介质流固耦合理论模型进行了比较。 相似文献
9.
10.
11.
Xikui Li Zejia Liu R. W. Lewis Kiichi Suzuki 《International journal for numerical methods in engineering》2003,57(6):875-898
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. 相似文献
12.
K. Mori M. Shiomi K. Osakada 《International journal for numerical methods in engineering》1998,42(5):847-856
The effect of microscopic rotation of powder particles in compaction is included in the rigid-plastic finite element method on the basis of the Cosserat continuum theory. In the Cosserat continuum theory, couple stress induced from the microscopic rotation is introduced, and the equilibrium equations of moment for the couple stress are solved simultaneously with those of force. A yield criterion for the Cosserat porous continuum is proposed by taking the effect of the couple stress into consideration, and constitutive equations for the rigid-plastic porous material are derived from the yield criterion on the basis of the associated flow rule. The equilibrium equations of force and moment for the Cosserat continuum are formulated by the use of the Galerkin method. The effect of microscopic rotation of powder particles in plane-strain closed-die compaction is examined. In addition, the calculated result is compared with that for the conventional continuum without the microscopic rotation. © 1998 John Wiley & Sons, Ltd. 相似文献
13.
14.
The finite element models of Cosserat continuum in two- and three-dimensions are presented. The size effects of a cantilever beam and a micro-rod, the well-posedness, the mesh-independent solutions of the boundary value problems with non-associated elastoplastic and strain softening constitutive behavior, and the progressive failure of the two- and three-dimensional vertical excavations are studied. Numerical results illustrate that the proposed Cosserat continuum models are capable of reflecting the size effects of micro-structures, preserving the well-posedness of the boundary value problem characterized by the strain localization, ensuring mesh-independent solutions, and simulating the entire progressive failure process occurring in engineering structures. 相似文献
15.
In this paper, a numerical model is developed for the fully coupled analysis of deforming porous media containing weak discontinuities which interact with the flow of two immiscible, compressible wetting and non-wetting pore fluids. The governing equations involving the coupled solid skeleton deformation and two-phase fluid flow in partially saturated porous media are derived within the framework of the generalized Biot theory. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three-phase formulation. The other variables are incorporated into the model via the experimentally determined functions that specify the relationship between the hydraulic properties of the porous medium, i.e. saturation, permeability and capillary pressure. The spatial discretization by making use of the extended finite element method (XFEM) and the time domain discretization by employing the generalized Newmark scheme yield the final system of fully coupled non-linear equations, which is solved using an iterative solution procedure. Numerical convergence analysis is carried out to study the approximation error and convergence rate of several enrichment strategies for bimaterial multiphase problems exhibiting a weak discontinuity in the displacement field across the material interface. It is confirmed that the problems which arise in the blending elements can have a significant effect on the accuracy and convergence rate of the solution. 相似文献
16.
Djordje Peri Jianguo Yu D. R. J. Owen 《International journal for numerical methods in engineering》1994,37(8):1351-1379
The a posteriori error estimates based on the post-processing approach are introduced for elastoplastic solids. The standard energy norm error estimate established for linear elliptic problems is generalized here to account for the presence of internal variables through the norm associated with the complementary free energy. This is known to represent a natural metric for the class of elastoplastic problems of evolution. In addition, the intrinsic dissipation functional is utilized as a basis for a complementary a posteriori error estimates. A posteriori error estimates and adaptive refinement techniques are applied to the finite element analysis of a strain localization problem. As a model problem, the constitutive equations describing a generalization of standard J2-elastoplasticity within the Cosserat continuum are used to overcome serious limitations exhibited by classical continuum models in the post-instability region. The proposed a posteriori error estimates are appropriately modified to account for the Cosserat continuum model and linked with adaptive techniques in order to simulate strain localization problems. Superior behaviour of the Cosserat continuum model in comparison to the classical continuum model is demonstrated through the finite element simulation of the localization in a plane strain tensile test for an elastopiastic softening material, resulting in convergent solutions with an h-refinement and almost uniform error distribution in all considered error norms. 相似文献
17.
18.
Nicola A. Nodargi Paolo Bisegna 《International journal for numerical methods in engineering》2019,120(11):1227-1248
A mixed membrane eight-node quadrilateral finite element for the analysis of masonry walls is presented. Assuming that a nonlinear and history-dependent 2D stress-strain constitutive law is used to model masonry material, the element derivation is based on a Hu-Washizu variational statement, involving displacement, strain, and stress fields as primary variables. As the behavior of masonry structures is often characterized by strain localization phenomena, due to strain softening at material level, a discontinuous, piecewise constant interpolation of the strain field is considered at element level, to capture highly nonlinear strain spatial distributions also within finite elements. Newton's method of solution is adopted for the element state determination problem. For avoiding pathological sensitivity to the finite element mesh, a novel algorithm is proposed to perform an integral-type nonlocal regularization of the constitutive equations in the present mixed formulation. By the comparison with competing serendipity displacement-based formulation, numerical simulations prove high performances of the proposed finite element, especially when coarse meshes are adopted. 相似文献
19.
Xikui Li Zejia Liu R. W. Lewis 《International journal for numerical methods in engineering》2005,64(5):667-708
A mixed finite element for coupled thermo‐hydro‐mechanical (THM) analysis in unsaturated porous media is proposed. Displacements, strains, the net stresses for the solid phase; pressures, pressure gradients, Darcy velocities for pore water and pore air phases; temperature, temperature gradients, the total heat flux are interpolated as independent variables. The weak form of the governing equations of coupled THM problems in porous media within the element is given on the basis of the Hu–Washizu three‐filed variational principle. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of the THM constitutive model for unsaturated porous media based on the CAP model. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elasto‐plastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is utilized. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization and the softening behaviours caused by thermal and chemical effects. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献