共查询到20条相似文献，搜索用时 46 毫秒
1.
This paper outlines the development as well as implementation of a numerical procedure for coupled finite element analysis of dynamic problems in geomechanics, particularly those involving large deformations and soil-structure interaction. The procedure is based on Biot’s theory for the dynamic behaviour of saturated porous media. The nonlinear behaviour of the solid phase of the soil is represented by either the Mohr Coulomb or Modified Cam Clay material model. The interface between soil and structure is modelled by the so-called node-to-segment contact method. The contact algorithm uses a penalty approach to enforce constraints and to prevent rigid body interpenetration. Moreover, the contact algorithm utilises a smooth discretisation of the contact surfaces to decrease numerical oscillations. An Arbitrary Lagrangian–Eulerian (ALE) scheme preserves the quality and topology of the finite element mesh throughout the numerical simulation. The generalised-α method is used to integrate the governing equations of motion in the time domain. Some aspects of the numerical procedure are validated by solving two benchmark problems. Subsequently, dynamic soil behaviour including the development of excess pore-water pressure due to the fast installation of a single pile and the penetration of a free falling torpedo anchor are studied. The numerical results indicate the robustness and applicability of the proposed method. Typical distributions of the predicted excess pore-water pressures generated due to the dynamic penetration of an object into a saturated soil are presented, revealing higher magnitudes of pore pressure at the face of the penetrometer and lower values along the shaft. A smooth discretisation of the contact interface between soil and structure is found to be a crucial factor to avoid severe oscillations in the predicted dynamic response of the soil. 相似文献
2.
《Computers and Geotechnics》1988,5(3):175-195
A method of analysis is presented for problems in which the deformation of the individual blocks play a significant role. The blocks are modeled as single quadrilateral elements and a constitutive model has been presented for computing the contact forces. In order to illustrate the influence of the deformation of the individual blocks two examples have been presented. The results of the analysis with deformable blocks and with rigid blocks are compared. These examples clearly demonstrate the importance of the deformation of blocks in the class of problems represented by these two examples. 相似文献
3.
Simulation of frictional contact between soils and rigid or deformable structure in the framework of smoothed particle hydrodynamics (SPH) is presented in this study. Two algorithms are implemented into the SPH code to describe contact behavior, where the contact forces are calculated using the law of conservation of momentum based on ideal plastic collision or using the criteria of partial penetrating. In both algorithms, the problem of boundary deficiency inherited from SPH is properly handled so that the particles located at contact boundary can have precise acceleration, which is critical for contact detection. And the movement and rotation of the rigid structure are taken into account so that it is easy to simulate the process of pile driving or movement of a retaining wall in geotechnical engineering analysis. Furthermore, the capability of modeling deformability of a structure during frictional contact simulations broadens the fields of SPH application. In contrast to previous work dealing with contact in SPH, which usually use particle‐to‐particle contact or ignoring sliding between particles and solid structure, the method proposed here is more efficient and accurate, and it is suitable to simulate interaction between soft materials and rigid or deformable structures, which are very common in geotechnical engineering. A number of numerical tests are carried out to verify the accuracy and stability of the proposed algorithms, and their results are compared with analytical solutions or results from finite element method analysis. Good agreement obtained from these comparisons suggests that the proposed algorithms are robust and can be applied to extend the capability of SPH in solving geotechnical problems. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
4.
将基于圆化多边形离散单元法与有限元方法结合,提出一种可变形圆化多边形离散单元法。此法对块体离散元进行圆化处理,可较好地表征不规则块体外形,又保留了颗粒离散元计算高效的优势。在求解接触力时,消除了角点处法向奇异等问题,同时增强计算的稳定性和简化接触判断。同时对切向接触力计算模型进行修正,使得接触力计算效率得到提高。此法突破了圆化多边形刚体假设的限制,可以精确计算任意形状不规则离散单元之间的相互作用,对单元的运动和变形进行模拟。通过超静定梁冲击试验、不规则块体单轴压缩试验和料斗流动“卡阻”试验3个数值模拟算例,论证此法可以有效地捕捉单元的碰撞、分离和变形等空间运动和自身特性以及其细观力学表征。 相似文献
5.
The paper examines the axisymmetric problem of the indentation of a poroelastic halfspace that is reinforced with an inextensible permeable/impermeable membrane located at a finite depth by a rigid indenter. The constitutive behavior of the poroelastic halfspace is described by the three-dimensional theory of poroelasticity proposed by M.A. Biot. The contact conditions between the indenter and the poroelastic halfspace are varied to accommodate both adhesive/frictionless contact and impermeable/permeable conditions. The formulation of the mixed boundary value problems uses the stress function approaches applicable to semi-infinite domains. Successive applications of Laplace and Hankel integral transforms are used to reduce the mixed boundary value problems to sets of coupled Fredholm integral equations of the second kind. These integral equations are solved using numerical approaches, applicable both for the solution of the systems of coupled equations and for Laplace transform inversion, to examine the time-dependent displacement of the rigid indenter. The analytical-numerical estimates for the time-dependent displacements of the rigid indenter are compared with results obtained using a finite element approach. 相似文献
6.
Many low‐order displacement‐based finite elements with exact integration are not suitable for estimating collapse loads of undrained geotechnical problems, especially for axisymmetric cases. As a result, higher‐order elements have to be used for these situations. In this technical note, the enhanced assumed strain (EAS) finite element method proposed by Simo and Rifai for elasticity problems are extended to plasticity problems to determine collapse loads. The numerical results for the problem of a smooth rigid surface footing on a deep purely cohesive undrained soil layer are given. It is demonstrated that the four‐noded quadrilateral EAS finite element is capable of estimating the collapse loads accurately for both undrained plane strain and axisymmetric problems. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
7.
数值流形法（NMM）在接触处理过程包括：接触方式判断,包括点-点和点-线接触;开-闭迭代,确定块体系统的约束状态;以及接触传递,将有效的接触对传递到下一个时步。过程稍显繁琐和耗时,且与真实的物理接触状态有异。由Munjiza所提出的NBS（no binary search）接触检测算法,将单元映射到规则格子中,以链表结构将其有效地连接在一起,只在单元所在格子以及周围格子内部进行接触判断,接触检测效率大为提高,计算量仅随单元数线性增长,内存需求也很低。在计算接触力时以所定义的势为媒介,用重叠面积来衡量接触力的大小,属于分布式接触力,更接近于实际,避免了原NMM接触处理过程的繁琐。因此,将其作为一种平行的接触处理方法引入到NMM中,并以算例验证,证实了算法的可行性。 相似文献
8.
模拟三维裂纹问题的扩展有限元法 总被引：4，自引：1，他引：3
扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其计算网格与不连续面相互独立,因此模拟移动不连续面时无需对网格进行重新剖分。给出了模拟三维裂纹问题的扩展有限元法。在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性。用两个水平集函数表示裂纹。采用线性互补法求解裂纹面非线性接触条件,不需要迭代,提高了计算效率。采用两点位移外推法计算裂纹前缘应力强度因子。给出了3个三维弹性静力问题算例,其结果显示了所提方法能获得高精度的应力强度因子,并能有效地处理裂纹面间的接触问题,同时表明扩展有限元结合线性互补法求解不连续问题具有较好的前景。 相似文献
9.
This paper describes a three-dimensional discrete element method of analysis for granular materials. Individual particles are modelled as rigid, six-sided solids with six degrees of freedom. Various types of contacts are possible among the three-dimensional discrete elements. Models for these contact types are described. Computational aspects of the method, such as techniques for contact detection and data structure, are addressed. Selective results of three examples that demonstrate the application of the method to granular flow problems are presented and discussed. 相似文献
10.
11.
The numerical simulation of soil-pile interaction problems, by means of full 3D finite element models, involves a large number of degrees of freedom (DOF) and difficulties during the mesh generation process. In order to reduce the unknowns and simplify and properly analyze such class of geotechnical problems, the so-called embedded beam elements (EBE) have recently been developed. In a preceding contribution of the authors, an improved EBE formulation, which brings into play the soil-pile interaction surface, was proposed with the aim to localize material plasticity in the soil surrounding the pile. This embedded beam model couples two different finite elements, each described by distinct kinematics (ie, solid and beam). The coupling is incorporated in the formulation by means of kinematical constrains established over the solid and beam displacement fields on the interaction surface. One of the main advantages of the embedded elements is that the addition of beams structural members immersed within the 3D soil model does not represent a constraint for the solid mesh, which can be adopted independently from the beam mesh. In this paper, the lateral loading of pile groups is studied by means of the proposed EBE approach with elasto-plastic interfaces. In order to represent a rigid cap, a master node and a special set of kinematical restrictions are incorporated into the formulation. The paper presents results obtained by means of the present formulation compared against other well-established analysis methods and test results published in the literature, for both elastic and elasto-plastic cases. 相似文献
12.
A smoothed particle hydrodynamics (SPH) framework for three-dimensional dynamic soil-multibody interaction modeling is presented, where both soils and rigid bodies are discretized using SPH particles. In the framework, soils are modeled using the Drucker-Prager model, while rigid bodies are considered with a multibody dynamics solver. A hybrid contact method suitable for three-dimensional simulations is developed to model the soil-body and body-body frictionless and frictional contacts, where contact forces are calculated based on ideal plastic collision and the unit normal/tangential vectors of the actual surface. Owing to its simplicity in contact detection and accuracy in contact force calculation, the hybrid contact method can be easily incorporated into SPH. Furthermore, graphics processing unit (GPU) parallelization is utilized to improve efficiency. The presented numerical framework and the hybrid contact method are validated using several examples. Numerical results are compared with analytical solutions and results from the literature. Furthermore, two three-dimensional simulations involving dynamic soil-multibody interaction are included to demonstrate the application. 相似文献
13.
14.
Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method 
下载免费PDF全文
下载免费PDF全文 A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
15.
In geotechnical engineering, numerical analysis of pile capacity is often performed in such a way that piles are modeled using only the geometry of their final position in the ground and simply loaded to failure. In these analyses, the stress changes caused by the pile installation are neglected, irrespective of the installation method. For displacement piles, which are either pushed or hammered into the ground, such an approach is a very crude simplification. To model the entire installation process of displacement piles a number of additional nonlinear effects need to be considered. As the soil adjacent to the pile is displaced significantly, small deformation theory is no longer applicable and a large deformation finite element formulation is required. In addition, the continuously changing interface between the pile and the soil has to be considered. Recently, large deformation frictional contact has been used to model the pile installation and cone penetration processes. However, one significant limitation of the analysis was the use of linear elements, which have proven to be less accurate than higher order elements for nonlinear materials such as soils.
This paper presents a large deformation frictional contact formulation which can be coupled consistently with quadratic solid elements. The formulation uses the so-called mortar-type discretisation of the contact surfaces. The performance of this contact discretisation technique is demonstrated by accurately predicting the stress transfer between the pile and the soil surfaces. 相似文献
16.
17.
土-结构-流体动力相互作用的实时耦联动力试验 总被引：3，自引：0，他引：3
针对振动台试验中无限地基难以模拟和数值分析中流-固耦合作用难以计算两个难题,将最近发展的实时耦联动力试验方法引入土-结构-流体动力相互作用问题的研究。以一个渡槽结构为例,其中渡槽-水体作为物理子结构,采用振动台进行物理试验,而无限地基作为数值子结构,采用集总参数模型进行数值模拟。两个子结构之间实时交换数据,联合评估整个耦合体系的动力响应。试验结果和有限元数值模拟结果吻合良好,表明该试验方法具有较高精度。对不同特性地基土进行的试验对比分析结果表明：对于软土地基,考虑土-结构相互作用(SSI)的结构反应幅值明显减小,周期延长;随着地基土变硬,SSI效应逐渐变弱,结构反应最终收敛至刚性地基解。 相似文献
18.
19.
Coupling of solid deformation and pore pressure for undrained deformation—a discrete element method approach 下载免费PDF全文
下载免费PDF全文 This paper presents a numerical scheme for fluid‐particle coupling that uses the discrete element method by taking into consideration solid deformation and pore pressure generation. A new water particle element is introduced to calculate pore water pressure due to porosity changes. The water particle element has the same size and shape as the solid element and experiences the same amount of deformation. On the basis of the effective stress principle at the element contact, the total force is equal to the sum of the force transmitted through the solid element contact and the water particle force due to pore water pressure. Analytical solutions of traditional soil mechanics problems, such as isotropic compression and consolidated triaxial undrained test, are used to quantitatively validate the proposed model. The numerical results show good agreement between the model and the analytical solutions. The model therefore provides an effective method to calculate pore pressure in a porous medium in discrete modeling. 相似文献
20.
依据空气单元模拟排水孔的思路,在渗流场任意可能逸出边界外布置一层渗透性较强的单元,与其他实体单元一起进行渗流计算,可以实现该边界的自动求解。与现行方法相比,不仅省去了逸出边界迭代计算的麻烦,而且避免了由于逸出点定位不准确而可能引起的渗流计算不合理或收敛困难等缺陷。算例分析表明,空气单元法的计算精度与空气单元的相对渗透系数R有关,与厚度L无关。当选取了合适的R值后,可以较好地逼近现行方法。对于因引入空气单元而造成的强弱透水介质界面处数值拟真性差的现象,还提出了加入接触面薄层单元的改进方法。方法改进后,计算精度可进一步提高。最后,以小湾水电站22#坝段渗流场计算为例,说明该法的实用性和可行性。 相似文献