首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 921 毫秒
1.
无单元伽辽金法及其在瞬态温度场中的应用研究   总被引:2,自引:0,他引:2  
无单元伽辽金法(EFGM)采用移动的最小二乘法构造形函数,和有限元相比,它只需结点信息而不需要单元信息.简述了无单元法的基础理论,推导出瞬态温度场的无单元法计算公式,采用罚函数法引入了第一类边界条件,编制了相应的计算程序.通过应用于经典的瞬态温度场例子,和有限元结果作比较,说明了无单元法具有精度高、前后处理简单等优越性,是一种具有较大发展潜力的新数值计算方法.  相似文献   

2.
王志亮  禹峰  殷宗泽 《岩土力学》2004,25(2):199-202
局部点插值法是一种新型的无单元法,该法构造出的位移和孔隙水压力的形函数具有 函数特性,所以,本质边界条件很易处理。首先,介绍了基本理论,给出求解比奥固结问题的主要公式。其次,对复层软基堆载引起的固结和沉降问题作了计算,其中节点采用均匀分布方式,且为了考虑土体的非线性特性计算程序中耦合了邓肯-张模型。最后,把计算值与有限元法、一维固结理论计算结果进行了对比。结果表明该法精度高,令人满意。  相似文献   

3.
未打穿散体材料桩复合地基固结简化解研究   总被引:1,自引:0,他引:1  
张玉国  孙敏  谢康和 《岩土力学》2014,35(5):1303-1310
利用平均固结度普遍解,将未打穿散体材料桩复合地基转化为等效双层地基,推导出未打穿散体材料桩复合地基固结简化解,给出其平均固结度计算简化公式。根据简化解、现有解和有限单元法编制程序,绘制考虑贯入比、桩体渗透系数、土体扰动区渗透系数和土体扰动区影响因子等因素对固结的影响曲线图,分析了各因素对未打穿散体材料桩复合地基固结性状的影响,并对3种计算方法所求出的平均固结度进行了比较研究,进而探讨了简化解的适用性。研究结果表明,无论是在单面排水还是在双面排水条件下,简化解求出的平均固结度总体上是介于现有解和有限单元法二者之间,尤其在固结过程的“中后期”(即固结度≥50%时)。从实用角度来看,简化解的计算结果具有足够精度,且简便、实用;在各种影响因素中,散体材料桩长度和桩体渗透系数对固结过程的影响最大。  相似文献   

4.
地下水流数值模拟是定量研究地下水水量和水质的重要手段。为揭示有限单元法水头反常现象形成的原因,以三维非均质模型为例,刻画了水头反常现象。储量集中法(Lumped Mass方法)是有限单元法中储释水量计算的替代算法,可以改善水头反常现象。对结点控制区域的均衡性质进行了研究,分析了使用Lumped Mass方法时的水头反常现象,结果表明:结点控制区域的对流量计算受到该结点所在层及其相邻层水力坡度的影响;在抽水初始阶段,抽水井结点水头快速下降,导致相邻层中相邻结点控制区域的侧向流入量为正值,控制区域为储水过程,结点水头上升,表现出水头反常现象。因此,有限单元法中对流量的计算是引起水头反常现象的原因。  相似文献   

5.
无单元法在三维断裂力学中的应用   总被引:5,自引:1,他引:5  
介绍了无单元法在三维断裂力学问题中的应用。无单元法基于滑动最小二乘拟合,只需结点信息就可建立离散模型。采用罚函数法来满足无单元法的本征边界条件。用可视准则来处理裂纹面对高斯点影响域的隔离作用,推导了裂纹两表面的相互作用对整体平衡方程组的贡献方程。利用单点位移公式来计算三维裂纹前缘点的应力强度因子。通过算例验证了应用该法确定三维裂纹前缘应力强度因子的有效性。  相似文献   

6.
裂隙岩体因含有发育程度不同的裂隙、节理和断层等不连续面,致其渗透性具有各向异性、不连续性等特点,因此传统的有限元法对分布密集的裂隙岩体渗流场求解有一定的难度。本文提出了采用无单元Glaerkin法求解有自由面裂隙渗流问题,并推导了无单元法求解渗流场的基本方程和积分格式,给出了应用罚函数法处理渗流边界条件和自由面处理方法。采用IDL语言编制了二维无单元法计算软件LIDAREFM。文中以北京怀柔桥梓镇某裂隙岩体边坡渗流场计算为例,研究了复杂裂隙共同作用下渗流场特性和自由面分布,讨论了不同开度、不同连通程度的裂隙对渗流场的影响。研究结果表明:无单元法可以较好地解决有密集裂隙的岩体渗流场的求解问题,实现了裂隙处结点任意加密以及积分网格的独立布置,避免了对有自由面和裂隙穿越的子域的重新处理,简化了渗流问题的求解过程。  相似文献   

7.
无网格自然单元法在弹塑性分析中的应用   总被引:8,自引:1,他引:7  
无网格自然单元法在构造位移插值函数时不需要单元的信息,只需要结点的信息,在裂纹扩展模拟、材料非线性分析、几何非线性分析以及三维计算等方面具有广阔的应用前景。阐述了将无网格自然单元法应用于结构弹塑性分析的过程和基本理论,给出了其位移插值函数的构造过程,并将其结果与常规的有限单元法进行对比,证明了该方法用于弹塑性分析的优越性。  相似文献   

8.
无网格伽辽金法应用的参数选择及内部边界处理   总被引:8,自引:0,他引:8  
无网格辽金法作为-种新的岩土工程数值计算方法, 该法其只需节点信息的无单元特性, 使其具有计算优势。本文结合固结EFGM刚度矩阵公式, 对不同的计算参数进行计算分析, 找出其影响规律。并采用跳跃函数处理内部边界条件, 计算结果表明, EFGM处理内部场函数不连续是准确的。  相似文献   

9.
EFG法在土体固结中的应用   总被引:4,自引:0,他引:4  
EFGM作为一种新的数值计算方法,具有只需节点信息而无段单元的特性,故在解固结方程方面有很大的优势。在计算分析中,此法容易构造固结方程的EFGM刚度矩阵和处理不同边界条件。对单面排水等条件的计算结果表明,EFGM在解决固结变形问题上,精度较高,处理边界准确。  相似文献   

10.
软黏土固结系数确定方法探讨   总被引:2,自引:0,他引:2  
根据揭普高速原状软土的室内固结试验结果和现场实测结果,对7种固结系数确定方法的原理、精度、使用简便性、优缺点等进行了比较研究。研究结果表明,三点法和时间平方根法是较好的方法,但采用时间平方根法时其固结试验应采用较大的加载比。研究还发现,同一场地单元、同类软土、不同土样的固结系数离散性较大。实际工程中,同一单元地段应进行多组固结试验计算固结系数,剔除较高值者,取余下均值为该地段的固结系数。  相似文献   

11.
基于Terzaghi一维固结理论,分析了考虑半透水边界条件的分数阶导数黏弹性饱和土层在随时间变化的任意荷载作用下一维固结问题。首先,应用Laplace变换联立求解饱和土层一维固结微分方程和分数阶Kelvin-Voigt黏弹性本构方程,推导出有效应力和沉降在Laplace变换域内的解析解,采用Crump方法进行Laplace逆变换,得到了时间域内的半解析解。然后将本文得到的半解析解分别退化为半透水边界条件下基于黏弹性假设的一维固结半解析解和双面透水边界条件下基于分数阶黏弹性假设的一维固结半解析解,结果与已有文献的半解析解相同,验证了本研究所提出解的可靠性。最后通过算例分别考察了半透水边界参数、分数阶黏弹性模型参数和荷载参数对饱和土层固结沉降的影响。研究表明,半透水边界条件参数、分数阶次与黏滞系数主要影响饱和土层固结的发展快慢,而饱和土层的最终沉降量主要受到土层压缩模量的影响;另外,饱和土层的固结规律与外荷载变化规律一致。  相似文献   

12.
This paper presents an exact analytical solution to fully coupled axisymmetric consolidation of a semi‐infinite, transversely isotropic saturated soil subjected to a uniform circular loading at the ground surface. The analysis is under the framework of Biot's general theory of consolidation. First, the governing equations of consolidation are transformed into a set of equivalent partial differential equations with the introduction of two auxiliary variables. These partial differential equations are then solved using Hankel–Laplace integral transforms. Once solutions in the transformed domain have been obtained, the actual solutions in the physical domain for displacements and stress components of the solid matrix, pore‐water pressure and fluid discharge can be finally obtained by direct numerical inversion. The accuracy of the numerical solutions developed is confirmed by comparison with an existing exact solution for an isotropic and saturated soil that is a special case of the more general problem addressed. Numerical analyses are also presented to investigate the influence of the degree of material anisotropy on the consolidation settlement. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

13.
邓俊晔  徐卫亚 《岩土力学》2007,28(9):1841-1847
研究了无单元法中单个影响域内的布点形式对误差的影响,提出了半径权值的概念。在经过充分的数值计算后,给出在半径权值对无单元法插值结果的精度的影响值的大小。研究对优化无单元法的节点布置、减小模拟误差有很大的作用,并且使得影响域的大小不再依靠点数的多少来确定,避免了影响域过大的情况和为了单个影响域内有足够多的插值点而使得整体插值点数目过多的弊端。运用单个区域内的插值点健康度的理论,研究了整个区域布点对无单元法模拟函数--滑动最小二乘法模拟精度的影响,并运用遗传算法对整个布点区域上的插值点进行了优化,给出了区域上的优化值。通过一个实例,得到一个具有普遍意义的插值点优化值与插值点密度的关系,并给出无单元法布点的一般性原则。  相似文献   

14.
The study presents semi-analytical solutions of two-dimensional plane strain consolidation problem in unsaturated soils incorporating the lateral semipermeable drainage boundary by adopting Fourier sine series and Laplace transform. The two-dimensional plane strain consolidation equations in the form of two-order partial differential equations with three variables are firstly converted to two-order partial differential equations with two variables, which are similar to those of one-dimensional consolidation problem. The four-order ordinary differential equations about excess pore-air and excess pore-water pressures are got by applying Laplace transform and the substitution method. Then, the solutions of excess pore pressures and settlement are achieved in the Laplace transform domain. Afterwards, on the basis of Crump's method, the inverse Laplace transform is conducted to obtain the analytical solutions in time domain. The comparison is conducted to verify the exactness of the obtained solutions, and the two-dimensional plane strain consolidation property with the lateral semipermeable drainage boundary is illustrated and discussed. Parametric studies are demonstrated for the excess pore pressures and normalized settlement with the change of the boundary parameters, air-water and lateral-vertical permeability coefficients, and the distance and depth. It can be found that the lateral semipermeable drainage boundary impedes the consolidation rate obviously, and when different investigated parameters are adopted, the consolidation property is similar to each other under the later permeable and semipermeable drainage boundary conditions.  相似文献   

15.
Axisymmetric consolidation of a poroelastic soil layer with a compressible fluid constituent induced by groundwater drawdown was studied based on Biot’s axisymmetric consolidation theory. Laplace and Hankel transforms were employed to solve the governing equation. Explicit analytical solutions are obtained in the Laplace–Hankel transform domain when groundwater drawdown is induced by a constant pumping well. Based on the solutions, numerical computations were performed to study the influences of the compressibility of the fluid constituent on the consolidation behavior of the soil layer.  相似文献   

16.
汪磊  李林忠  徐永福  夏小和  孙德安 《岩土力学》2018,39(11):4142-4148
为描述饱和土体的流变特性,引入分数阶导数Kelvin-Voigt黏弹性模型,采用解析方法对半透水边界下的分数阶黏弹性饱和土一维固结特性进行了研究。分别对骤加恒载下饱和土一维固结微分方程和分数阶Kelvin-Voigt黏弹性本构方程进行Laplace变换,并联立求解得到了双边半透水边界条件下分数阶黏弹性饱和土在Laplace变换域内的解析表达式。通过Crump方法实现Laplace数值反演,得到时间域内的半解析解。将所得到的解分别退化为分数阶黏弹性饱和土一维固结半解析解和双边半透水黏弹性饱和土一维固结半解析解,结果与已有文献半解析解相同,验证了提出的双边半透水边界条件下分数阶黏弹性饱和土一维固结解的可靠性。通过算例考察了半透水边界条件和分数阶黏弹性饱和土参数对一维固结特性的影响。研究表明,双边半透水边界下分数阶黏弹性饱和土一维固结发展过程与半透水边界条件、分数阶次和黏滞系数有关,且土体的压缩模量对饱和土一维固结最终沉降量有显著影响。  相似文献   

17.
This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.  相似文献   

18.
INTRODUCTION Theconsolidationofsaturatedsoilisanimpor tantsubjectingeomechanics.Biot(1941)proposed thethree dimensionalconsolidationtheory,whichis basedontheprincipleofeffectivestress,continuity conditionsandequilibriumequations.Theremarka blesuccessofthistheoryistheanalysisofthetime dependenteffectandcouplingeffectofsoilandpore water.Thus,thetheoryisaclassictheoryinsatu ratedsoilstatistics.Itwasalmostunfeasibletoperformanalysesof thistheorywithoutthedevelopmentofhigh speed computersandnu…  相似文献   

19.
This paper presents an analytical layer element solution to axisymmetric thermal consolidation of multilayered porous thermoelastic media containing a deep buried heat source. By applying the Laplace–Hankel transform to the state variables involved in the basic governing equations of porous thermoelasticity, the analytical layer elements that describe the relationship between the transformed generalized stresses and displacements of a finite layer and a half‐space are derived. The global stiffness matrix equation is obtained by assembling the interrelated layer elements, and the real solutions in the physical domain are achieved by numerical inversion of the Laplace–Hankel transform after obtaining the solutions in the transformed domain. Finally, numerical calculations are performed to demonstrate the accuracy of this method and to investigate the influence of heat source's types, layering, and the porous thermoelastic material parameters on thermal consolidation behavior. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号