共查询到20条相似文献，搜索用时 15 毫秒
1.
The development of a numerical procedure for the finite element analysis of anchors dynamically penetrating into saturated soils is outlined, highlighting its unique features and capabilities. The mechanical behaviour of saturated porous media is predicted using mixture theory. An algorithm is developed for frictional contact in terms of effective normal stress. The contact formulation is based on a mortar segment-to-segment scheme, which considers the interpolation functions of the contact elements to be of order N, thus overcoming a numerical deficiency of the so-called node-to-segment (NTS) contact algorithm. The nonlinear behaviour of the solid constituent is captured by the Modified Cam Clay soil model. The soil constitutive model is also adapted so as to incorporate the dependence of clay strength on strain rate. An appropriate energy-absorbing boundary is used to eliminate possible wave reflections from the artificial mesh boundaries. To illustrate the use of the proposed computational scheme, simulations of dynamically penetrating anchors are conducted. Results are presented and discussed for the installation phase followed by ‘setup’, i.e., pore pressure dissipation and soil consolidation. The results, in particular, reveal the effects of strain rate on the generation of excess pore pressure, bearing resistance and frictional forces. The setup analyses also illustrate the pattern in which pore pressures are dissipated within the soil domain after installation. Hole closure behind a dynamic projectile is also illustrated by an example. 相似文献
2.
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. 相似文献
3.
《Computers and Geotechnics》2006,33(6-7):341-354
This paper presents a numerical formulation for frictional contact problems associated with pile penetration. The frictional contact at the soil–pile interface is formulated using the theory of hardening/softening plasticity, so that advanced models for the interface can be dealt with. A smooth discretisation of the pile surface is proposed using Bézier polynomials. An automatic load stepping scheme is proposed, which features an error control algorithm and automatic subincrementation of the load increments. The numerical algorithms are then used to analyse the installation process of pushed-in axial piles. It is shown that the smooth discretisation of the pile surface is effective in reducing the oscillation in the predicted pile resistances and the automatic load stepping scheme outperforms the classical Newton–Raphson scheme for this type of problem. 相似文献
4.
In this paper, a numerical procedure based on the finite element method is outlined to investigate pile behaviour in sloping ground, which involves two main steps. First a free-field ground response analysis is carried out using an effective stress based stress path model to obtain the ground displacements, and the degraded soil stiffness and strength over the depth of the soil deposit. Next a dynamic analysis is carried out for the pile. The interaction coefficients and ultimate lateral pressure of soil at the pile–soil interface are calculated using degraded soil stiffness and strength due to build-up of pore pressures, and the soil in the far field is represented by the displacements calculated from the free-field ground response analysis. Pore pressure generation and liquefaction strength of the soil predicted by the stress path model used in the free-field ground response analysis are compared with a series of simple shear tests performed on loose sand with and without an initial static shear stress simulating sloping and level ground conditions, respectively. Also the numerical procedure utilised for the analysis of pile behaviour has been verified using centrifuge data, where soil liquefaction has been observed in laterally spreading sloping ground. It is demonstrated that the new method gives good estimate of pile behaviour, despite its relative simplicity. 相似文献
5.
Y. K. Chow 《国际地质力学数值与分析法杂志》1987,11(1):1-15
A numerical procedure is described for the analysis of vertical deformation of smooth, rigid foundations of arbitrary shape on homogeneous and layered soil media. The contact area at the interface of the foundation and soil medium is approximated by square subdivisions. The response of the system is then obtained from the superposition of the influence of the individual subdivisions. The flexibility influence coefficients are based on equivalent smooth, rigid circular areas with the same contact area as the square subdivisions. For foundations on a homogeneous, isotropic elastic half-space, the flexibility coefficients are given analytically by the integrated forms of the Boussinesq's solution. For a layered soil medium, the flexibility coefficients are determined from an axisymmetric finite element analysis which is essentially two dimensional. Thus, there is no necessity for a full three-dimensional finite element analysis. Comparison with solutions obtained using the integral transform technique for smooth, rigid rectangular foundations on a homogeneous, isotropic elastic half-space shows good agreement. Parametric solutions are presented for the response of rectangular foundations on some ‘typical’ soil profiles. The use of a simplified method to estimate the settlement of rectangular foundations on a layered soil medium by superposing solutions for homogeneous, elastic strata is discussed. 相似文献
6.
非饱和土层一维固结特性分析 总被引：1，自引：0，他引：1
在Fredlund非饱和土的一维固结理论的基础上进行假设,由得出的液相及气相的控制方程、Darcy定律及Fick定律,采用Laplace 变换、逆变换等数学方法得到了大面积均布瞬时加载下表面为透水透气面、底面为不透水和不透气面的非饱和土层一维固结时间域内的超孔隙水压力、超孔隙气压力及土层沉降的解析解;应用典型算例,分析了不同气、水渗透系数比情况下土体超孔隙水压力、超孔隙气压力消散及土层沉降随时间的变化规律以及不同时间超孔隙水压力、超孔隙气压力消散随深度的变化规律。将得出的结果退化成相应的饱和土的解与太沙基饱和土固结理论结果比较,验证了其正确性。 相似文献
7.
The equations governing the dynamic behavior of saturated porous media as well as a finite element spatial discretization of these equations are summarized. A three-parameter time integration scheme called the Hilber–Hughes–Taylor α-method is used together with a predictor/multi-corrector algorithm, instead of the widely used Newmark's method, to integrate the spatially discrete finite element equations. The new time integration scheme possess quadratic accuracy and desirable numerical damping characteristics. The proposed numerical solution and bounding surface plasticity theory to describe the constitutive behaviour of soil have been implemented as the computer code DYSAC2. Predictions made by DYSAC2 code are verified using dynamic centrifuge test results for a clay embankment. Importance of initial state of a soil on its dynamic behaviour is demonstrated. 相似文献
8.
桩承式加筋路堤的现场试验及数值分析 总被引：4，自引：0，他引：4
对一桩体面积置换率为8.7 %的低置换率桩承式加筋路堤进行了现场试验及三维有限元分析。现场主要进行了桩、土荷载分担,孔压、沉降及侧向水平变形等内容的观测。将观测数据与常规设计方法及三维有限元分析结果进行了对比研究,在此基础上对设计方法的适用性进行了分析。研究结果表明,路堤填土的土拱效应造成荷载向桩体转移,这种荷载转移大幅度减小了在软土层中产生的超孔隙水压力。当填土高度大于2.5 m时,土拱效应的应力折减系数可用Russell和Pierpoint或Hewlett和Randolph提出的土拱效应分析方法进行计算,其结果与三维有限元分析也较相符,但在路堤高度较小时,只有Russell和Pierpoint方法与实测结果相接近。路堤施工过程中,实测的水平变形与沉降之比仅为0.2左右,这表明采用桩承式加筋路堤不仅可减小沉降,而且可减小水平向的变形,提高路堤的稳定性。 相似文献
9.
Coupled theory of mixtures for clayey soils 总被引：4，自引：0，他引：4
In this work, elasto-plastic coupled equations are formulated in order to describe the time-dependent deformation of saturated cohesive soils (two-phase state). Formulation of these equations is based on the principle of virtual work and the theory of mixtures for inelastic porous media. The theory of mixtures for a linear elastic porous skeleton was first developed by Biot (Theory of elasticity and consolidation for a porous anisotropic solid, Journal of Applied Physics, 1955, 26, 188–185). An extension of Biot's theory into a nonlinear inelastic media was performed by Prevost (Mechanics of continuous porous media, International Journal of Engineering Science, 1980, 18, 787–800). The saturated soil is considered as a mixture of two deformable media, the solid grains and the water. Each medium is regarded as a continuum and follows its own motion. The flow of pore-water through the voids is assumed to follow Darcy's law. The coupled equations are developed for large deformations with finite strains in an updated Lagrangian reference frame. The coupled behavior of the two-phase materials (soil-water state) is implemented in a finite element program. A modified Cam-clay model is adopted and implemented in the finite element program in order to describe the plastic behavior of clayey soils. Penetration of a piezocone penetrometer in soil is numerically simulated and implemented into a finite element program. The piezocone penetrometer is assumed to be infinitely stiff. The continuous penetration of the cone is simulated by applying an incremental vertical movement of the cone tip boundary. Results of the finite element numerical simulation are compared with experimental measurements conducted at Louisiana State University using the calibration chamber. The numerical simulation is carried out for two cases. In the first case, the interface friction between the soil and the piezocone penetrometer is neglected. In the second case, interface friction is assumed between the soil and the piezocone. The results of the numerical simulations are compared with experimental laboratory measurements. 相似文献
10.
Softening solids are analysed under impact loading using a new numerical method which allows displacement discontinuities to propagate arbitrarily through a finite element mesh. The Dirac‐delta distributions that arise in the strain field of classical continuum theory in the presence of strain softening are interpreted as discontinuities in the displacement field. A new finite element procedure with Heaviside jumps added to the underlying displacement interpolation basis is able to capture displacement jumps independent of the spatial discretisation. The amplitudes of displacement jumps are represented by extra degrees of freedom at existing nodes. Numerical results for mode‐I and mode‐II failure due to impact loading are presented. The numerical results highlight the objectivity of the approach with respect to spatial discretisation under dynamic loading conditions. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
11.
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. 相似文献
12.
在波浪循环荷载作用下饱和软黏土地基出现孔隙水压力升高,并导致不排水强度弱化,严重影响防波堤的承载性能。考虑静偏应力影响,基于最大孔隙水压力发展模型和正常固结软黏土不排水强度公式,推导出软黏土不排水强度随循环荷载作用次数和应力水平变化的动态折减规律。结合软黏土不排水强度动态折减规律和M-C屈服准则,在有限元软件ABAQUS上实现软黏土不排水强度循环弱化分析的数值开发和动力运算过程。运用该动力有限元方法对天津港防波堤地基软黏土的动、静三轴试验进行数值模拟运算。结果表明,最大孔隙水压力发展曲线以及循环荷载作用后不排水强度的数值预测结果与动三轴试验结果吻合良好。另外,动力有限元方法（DFEM）能够表示土体强度在循环荷载作用下的具体弱化过程。 相似文献
13.
Yi--> Tian Wenbing Wu Guosheng Jiang M. Hesham El Naggar Guoxiong Mei Meijuan Xu Rongzhu Liang 《国际地质力学数值与分析法杂志》2020,44(8):1170-1183
In engineering practice, a rapid loading rate can result in ground failure when the strength of soft soils is relatively low, and a multistage loading scheme is always utilized to deal with this situation. Firstly, under a multistage load and the continuous drainage boundary, an analytical solution of excess pore-water pressure and consolidation degree is obtained by virtue of the superposition formula of excess pore-water pressure, and a more general continuous drainage boundary under arbitrary time-dependent load is developed. Then, a comparison with existing analytical solutions is conducted to verify the present solution. A preliminary attempt on applying the continuous drainage boundary into the finite element model is made, and the feasibility of the numerical model for the one-dimensional consolidation under the continuous drainage boundary is verified by comparing the results calculated by FEM with that from present analytical solution. Finally, the consolidation behavior of soil is investigated in detail for different int erface parameters or loading scheme. The results show that, in land reclamation projects, a horizontal drain should be placed close to the boundary with a smaller interface parameter to improve the consolidation efficiency. The degree of consolidation is also related to the applied time-dependent load and interface parameters. 相似文献
14.
This paper presents a finite element model for analysing the behaviour of granular material wrapped with polyethylene bags under vertical compression and cyclic shearing. The simple Mohr–Coulomb model is used to represent the soil behaviour. The polyethylene bag is represented by a linear-elastic–perfect-plastic model. The soil-bag interface is modelled with contact constraints. The main purpose of the numerical analysis is to validate the anticipated performance of soilbags under various loading conditions and hence the effectiveness of soilbags as a method of ground improvement. 相似文献
15.
This is a summing-up report based on 48 papers on the applications of numerical method to soil mechanics. Recently, there is a new development of coupling of finite element with other numerical ones in soil mechanics. The dynamic numerical analyses have found further application to nonlinear stress-strain relationship, effective stress, and three dimensional problems. The suitable contact element or conjunct element could be used for contact conditions in finite element calculation. It is an important direction to combine usual numerical method with probabilistic method, which has been used in evaluations of consolidation effect, slope stability, dam, seepage, deformation calculation, interaction, and the performance analyses of pole and foundation. 相似文献
16.
数值计算方法在土力学中的应用 总被引：2，自引：0，他引：2
本文是对48篇有关数值方法在土力学方面应用的论文进行综述的报告。近年来,在土力学方面有限元和其它数值方法的耦合有了新的发展,动力分析在非线性应力应变关系、三维问题和有效应力分析中得到了进一步的应用,接触单元或连接单元可用来处理有限元计算时的接触条件。数值方法和概率方法的结合是一个值得重视的发展方向。这些方法可应用于估计固结效应、土坡稳定、土石坝、渗流、变形计算、共同作用、桩和地基性状分析等方面。 相似文献
17.
Panique Lazcano Daniel R. Galindo Aires Rubén Patiño Nieto Hernán 《Acta Geotechnica》2022,17(5):1897-1915
The calculation of the long-term dynamic bearing capacity arises from the need to consider the generation of maximum pore-water pressure developed from a cyclic load. Under suitable conditions, a long-term equilibrium situation would be reached, when pore-water pressures stabilized. However, excess pore-water pressure generation can lead to cyclic softening. Consequently, it is necessary to define both the cohesion and the internal friction angle to calculate the dynamic bearing capacity of a foundation in the long term, being necessary to incorporate the influence of the self-weight of soil and therefore the width of the foundation. The present work is based on an analysis of the results of cyclic simple shear tests on soil samples from the port of El Prat in Barcelona. From these experimental data, a pore-water pressure generation formulation was obtained that was implemented in FLAC2D finite difference software. A methodology was developed for the calculation of the maximum cyclic load that a footing can resist before the occurrence of the cyclic softening. The type of soil studied is a contractive cohesive soil, which generates positive pore-water pressures. As a numerical result, design charts have been developed for long-term dynamic bearing capacity calculation and the charts were validated with the application of a real case study.
相似文献18.
A finite element model for pile‐driving analysis is developed and used to investigate the behaviour of pre‐bored piles, which are then driven the last 1.25 or 2.25 m to their final design depth. The study was conducted for the case of saturated clays. The model traces the penetration of the pile into the soil and accommodates for large deformations. The non‐linear behaviour of the clay in this study is predicted using the bounding‐surface‐plasticity model, as applied to isotropic cohesive soils. The details of the 3‐D numerical modelling and computational schemes are presented. A significant difference was observed in the pile displacement during driving, and in the computed soil resistance at the pile tip, particularly at the earliest driving stages. No difference in soil resistance at the soil pile interface along the pile shaft was detected between the pre‐bored piles whether driven 1.25 or 2.25 m. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
19.
Y. K. Chow 《国际地质力学数值与分析法杂志》1987,11(6):621-638
A numerical method of analysis based on elasticity theory is presented for the analysis of axially and laterally loaded pile groups embedded in nonhomogeneous soils. The problem is decomposed into two systems, namely the group piles acted upon by external applied loads and pile–soil interaction forces, and a layered soil continuum acted upon by a system of pile–soil interaction forces at the imaginary positions of the piles. The group piles are discretized into discrete elements while the nonhomogeneous soil behaviour is determined from an economically viable finite element procedure. The load–deformation relationship of the pile group system is then determined by considering the equilibrium of the pile–soil interaction forces, and the compatibility of the pile and soil displacements. The influence of soil nonlinearity can be studied by limiting the soil forces at the pile–soil interface, and redistributing the ‘excess forces’ by an ‘initial stress’ process popular in elasto-plastic finite element analysis. The solutions from this approach are compared with some available published solutions for single piles and pile groups in homogeneous and nonhomogeneous soils. A limited number of field tests on pile groups are studied, and show that, in general, the computed response compares favourably with the field measurements. 相似文献
20.
A non-linear three-dimensional finite element procedur is developed and applied for the analysis of pile group foundations. The numerical procedure allows for elastic, non-linear elastic and elastic-plastic hardening behaviour of sand. In order to include the interaction effects involving relative slip and debonding, the thin-layer interface element is used. The predictions for displacements and loads obtained from the numerical procedure are compared with laboratory model test results of a pile group. Displacements, stresses and forces distribution in various components of the pile group are also examined. Furthermore, the effects of the non-linear soil response and relative motions at the interface are indentified and discussed. 相似文献