共查询到20条相似文献，搜索用时 15 毫秒
1.
大多数的岩土工程问题都会涉及到岩土体-结构之间的相互作用,它们之间界面的合理模拟对计算结果有重要的影响。目前常用的Goodman节理单元、Desai薄层单元等虽然简单好用,但也存在着明显的缺陷。相比之下,作为接触问题计算,可更为合理地模拟接触面的力学行为。针对变形体-刚体这种特殊的接触问题,应用变形体的虚功原理建立了相应的虚功方程,在此基础上得出有限元计算公式,并建立了弹塑性接触问题有限元计算的实现方法。通过对刚性条形基础基底压力分布和刚性挡墙土压力分布等算例的计算,验证了该法的正确性和求解工程问题的有效性。 相似文献
2.
3.
The paper presents some numerical results of the problem on contact of stratified elastic half‐spaces possessing unevenness of their surface layers. The approach of the Jacobi polynomials is used to reduce the singular integral equation of the problem to a system of linear algebraic equations. On the basis of the system's solution the dependence of gaps' geometric characteristics on external load and friction is investigated for two types of initial boundary disturbances, namely a local symmetric recess and a periodic set of recesses at a flat boundary. Numerical results are presented in figures for a few elastic and geometric parameters of joint components. The analysis shows that elastic properties and geometric structure of a laminated half‐space have a significant influence on the shape of interface gaps and their location. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
4.
接触问题是一种常见的非线性问题,如何能够很好地模拟接触面的变形及受力特性,以及实现对变形体间的接触问题的真实模拟是该领域研究的难点问题。基于二维接触问题的实际物理意义,分别在法向和切向建立等价的互补模型。用非线性互补函数（NCP）中的Fischer-Burmeister（FB）函数将互补函数模型转化为非光滑方程组表达,用常规的Newton法求解。同时,基于高斯积分法可以用较少的积分点达到较高的精度,为了进一步提高求解精度,改善不连续的通病,对面-面接触模型在高斯点上对接触面进行处理,可通过调节积分点数目对求解精度进行控制,方法易于理解,实现方便。在此基础上建立二维接触有限元模型,通过一系列工程算例验证该方法的可行性与有效性。结果表明,与ABAQUS有限元的计算结果相比,该方法有着较高的精度,更真实地反映问题的实际。 相似文献
5.
6.
This paper addresses the problem of characterizing the shape of a geological surface on the basis of its principal curvatures.
The surface is assumed to be modeled as a set of adjacent triangles defined by the location of their vertices and a method
is proposed for estimating numerically the principal curvatures at the vertices of the triangles using a local C2 interpolant. Also shown is how principal curvatures can be useful for studying the deformation of a geological surface (with
application to 3D balanced unfolding), and analyzing the folding or faulting of the interface between two adjacent layers. 相似文献
7.
8.
双曲线接触面本构模型能比较真实地模拟堆石材料与混凝土面板间的接触关系,扩展Lagrange乘子法能比较精确地计算接触面间接触状态。在吸取Clough双曲线接触面本构关系应用于无厚度的Goodman单元和有厚度的Desai薄层单元的成功经验的基础上,将Clough双曲线接触面模型引入扩展Lagrange乘子法进行摩擦接触问题的求解。介绍了基于扩展Lagrange乘子法的序列二次规划法（SQP）提法,详细推导了在扩展Lagrange算法的非线性接触计算中引入Clough双曲线接触面本构关系的数值列式和实施步骤。数值算例和工程实例的计算结果证明该方法是成功的,能够比较真实地模拟堆石-混凝土墙（混凝土面板）之间的接触状态。 相似文献
9.
三维变形体离散元法能够自动检索接触关系,并对具有不规则形状的粗粒土颗粒和结构物进行有限差分网格离散,因此,具有模拟离散-连续耦合问题的先天优势。采用随机模拟技术生成粗粒土三维数值试样,基于变形体离散元进行粗粒土与结构接触面特性的数值试验,研究了不同接触面粗糙程度的接触面力学特性,对比了粗粒土与结构物在单剪和直剪状态下的接触面力学特性,从宏观和细观两个层面分析了数值试验结果。结果表明,数值试验能较好地反映粗粒土与结构接触面的力学特性,其剪应力-相对剪切位移曲线与试验结果规律相似;接触面粗糙程度对接触面的强度和变形特性影响较大,其机制在于剪切对试样的扰动程度不同;直剪和单剪状态下试样剪应力-相对剪切位移均为双曲线,单剪试验的初始剪切刚度低于直剪试验,两种试验得到接触面抗剪强度指标比较接近。 相似文献
10.
Osvaldo Manzoli Marcelo Sánchez Michael Maedo Jumanah Hajjat Leonardo J. N. Guimarães 《Acta Geotechnica》2018,13(5):1171-1186
The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory. 相似文献
11.
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. 相似文献
12.
Philippe Leroy Arnault Lassin Mohamed Azaroual Laurent André 《Geochimica et cosmochimica acta》2010,74(19):5427-776
The surface tension of the air/water interface is a phenomenon of particular interest in the water-unsaturated zone of porous media because it influences the contact angle and consequently the capillary water volume. A mechanistic model based on the modified Poisson-Boltzmann equation and the Pitzer theory is described and used to predict, under isothermal and isobaric conditions, the surface tension of 1:1 electrolytes at high salinity. These theories enable the determination of the electrical potential at the air/water interface and the activity coefficient of the ionic species in the bulk pore water, respectively. Hydration free energies of the structure-making and structure-breaking ions that influence the surface tension at high salinity are taken into account. Structure-making ions flee the air/water surface because they can better organize the water dipoles in bulk water than at the interface. Structure-breaking ions are positively adsorbed at the air/water interface because the bulk water can better organize their hydrogen-bonding network without these ions. The resulting surface tension increases and decreases, respectively, compared to the surface tension of pure water. The predictions are in good agreement with the surface tension data of 1:1 electrolytes (NaCl, KCl, HCl, NaNO3, KNO3, HNO3 aqueous solutions) and the optimized parameters depend on the effective electrostatic diameters of cations and on the hydration free energies of the ions at the interface. 相似文献
13.
建立准确合理的接触面本构模型是研究桩基承载力、沉降及变形等问题的前提。考虑桩-土界面初始剪切刚度的深度效应,提出基于统计损伤本构模型的无厚度接触面本构模型,根据剪切试验的数据对改进的接触面本构模型参数进行拟合,得到了不同法向应力下的模拟曲线。利用FLAC3D软件内嵌的FISH语言,对桩-土界面进行改进本构模型的二次开发,提出的无厚度接触面模型合理地模拟出桩-土界面刚度的非线性,成功实现了桩-土剪切试验过程的数值模拟。计算结果表明：提出的无厚度接触面本构模型能够合理地模拟桩-土界面的力学行为,可以描述接触面刚度的非线性变化,也适于程序化计算,拓展了接触面统计损伤本构模型的应用领域,为深入探讨桩基承载变形性状奠定了基础。 相似文献
14.
为了研究堆积层-基岩接触面滑坡的形成机理,以陕西省岚皋县广泛分布的堆积层-基岩接触面滑坡为研究对象,采用野外调查和数理统计的方法,总结了该类型滑坡的基本特征,并以祖师庙滑坡为例,借助室内试验和数值模拟手段,研究了堆积层-基岩接触面滑坡的形成机理。结果表明:1)堆积层-基岩接触面滑坡主要发育在坡度为20°~40°、下伏基岩为碎裂状结构的斜坡,多为浅表层滑坡;2)碎石土和堆积层-基岩接触面的内摩擦角随含水量增大而降低,堆积层-基岩接触面是斜坡变形破坏的软弱结构面;3)细颗粒物质的增加降低了堆积层-基岩接触面的内摩擦角,且细颗粒粒径越小,内摩擦角也越小,在饱和状态下降幅度越大;4)因细颗粒物质在堆积层-基岩接触面处聚集造成的强度降低和因水文地质条件改变造成的滑带土饱和是祖师庙滑坡主要的诱发因素。研究结果对揭示堆积层-基岩接触面滑坡的形成机理和地质灾害防治具有一定的指导意义。 相似文献
15.
将地基视为多层各向同性弹性体,对Euler-Bernoulli梁进行有限单元离散分析,对地基-梁接触面采用边界积分法求解,根据地基-梁接触面的竖向位移协调和光滑接触条件,应用有限元与边界元耦合的方法推导出各向同性成层弹性地基上的Euler-Bernoulli梁的半数值半解析解。基于文中理论,编制了相应的程序,通过与现有文献对比,验证了文中理论的正确性,对比分析了分层地基与等效均匀地基两种模型上的弹性地基梁。研究结果表明,分层地基与等效均匀地基两种模型上弹性地基梁性状差异较大,实际弹性地基梁计算中应采用分层地基模型。 相似文献
16.
基于广义位势理论建立的岩土体材料本构模型以及岩土体材料与结构接触面本构模型原理相通,只是前者是在三轴剪切试验条件下的三维应力空间建模,后者是在单剪试验条件下的二维应力空间建模。单剪试验条件下土与结构的接触面问题可以看作是法向与切向应力空间上的二维问题,其试验结果可以表达成由应力、应变组成的二维矢量。结合接触面力学特性,确定应力空间中的势函数以及塑性状态方程,可以推导出双重势面接触面弹塑性本构方程的一般表达式。进一步取两个势函数分别为法向应力和切向应力,建立简化双重势面接触面弹塑性模型的本构方程,该方程可直接应用于有限元等数值分析。结合试验实例对建模方法的合理性进行验证,模型拟合效果良好。研究结果表明,基于广义塑性位势理论建立接触面本构模型无需推求塑性势函数和屈服函数,可以直接得到弹塑性刚度矩阵,且建模方便。 相似文献
17.
基于无黏性可压缩理想流体介质波动理论和Biot流体饱和多孔介质波动理论,考虑水下饱和土的流固耦合,借助Hankel函数积分变换法(HFITM)给出入射平面P1波在海底洞室周围散射问题的解析解。相比传统研究中的"大圆弧假定",Hankel函数积分变换法可以较好地处理半空间表面边界条件。利用该解析解,计算分析了洞室表面透水条件、入射角度、入射频率、海水水深和饱和土的孔隙率等因素对水-土交界面处水平位移、竖向位移和洞室表面动水压力、环周总应力的影响。结果表明:洞室表面透水条件对水-土交界面处水平位移和竖向位移影响较小;随着斜入射角的增加,水-土交界面处竖向位移减小;随着入射频率的增加,水-土交界面处水平位移随之增加;海水水深为2.5倍SV波的波长时,水-土交界面处水平位移及洞室表面动水压力最大值最大;随着孔隙率的增加,水-土交界面处水平位移、竖向位移和洞室表面环周总应力减小,而洞室表面动水压力随之增加。 相似文献
18.
Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
19.
By means of a semi-analytical FE approach an embedded circular footing under monotonic horizontal and moment loading is studied.
In a non-homogeneous soil whose shear modulus is characterized by a power law variation with depth, horizontal, rocking and
coupled modes of displacement, expressed in terms of influence factors are thoroughly examined. The exponent α that controls
the shape of the stiffness variation with depth is termed shear modulus factor. Surface influence coefficients are considered
for the situations where the interface between the soil and the footing is either perfectly rough or perfectly smooth. First,
results of semi-analytical FE analysis in the case of rough footing are established and compared with those of another numerical
method. Results of comparison show good agreement. Then, for different values of α the surface influence coefficients are
presented for an embedded footing in perfect smooth contact with soil. The metacentre is referred to as the depth at which
there is no coupling between the sliding and the rocking modes of footing deformations. Expressions for location and horizontal
influence coefficient corresponding to this particular depth are developed and their variations with α examined. The paper
finishes by showing the effect of interface conditions on the soil normal stresses developed beneath the embedded circular
footing for the case of loading applied at the footing top. 相似文献
20.
《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. 相似文献