Winkler地基上厚板分析的自然单元法.pdfVIP

第29卷第2期 固体力学学报 V01．29No．2 CHINESE 0FSOLIDMECHANICS 2008年6月 JOURNAL 2008 June Winkler地基上厚板分析的自然单元法。 曾祥勇p’ 朱爱军2 邓安福1 (1重庆大学土木工程学院，重庆，400045)(2贵州大学土木建筑工程学院，贵阳，550003) 摘 要 地基上厚板弯曲挠度的自然单元法求解控制方程，并进行了相应的程序实现，最后通过算例分析，表明了该文方法 的可行性和有效性． 关键词 自然单元法，Laplace插值函数，自然邻接点，Mindlin板，无网格法 张建辉，邓安福[14]首次应用EFGM计算分析筏板 0 引言 基础和弹性地基板，取得较好的成果．张伟星，庞辉 等[1朝应用无单元Galerkin法研究分析了弹性半空 自然单元法作为目前众多无网格方法中一种独 间地基板的弯曲问题．方电新，李卧东，王元汉[1钉用 具特色的方法，近年来受到了国内外研究人员的广 无单元Galerkin法计算了平板弯曲问题．熊渊博 elementmethod) 泛关注[1-z]．自然单元法(natural 的概念由Braun与Sambridge[33于1995年在国际弹性地基板分析． 著名刊物Nature上首次提出并对这种新型的数值 与其它基于代数方法构造形函数的无网格方法 方法进行了成功的应用．该方法是一种基于给定结 (如EFGM等)相比，本文自然单元法形函数构造基 点的Voronoi图及Delaunay三角形剖分图，以自然于自身独特的几何形式，其形函数数学表达式均为 邻接点插值函数作为试函数的Galerkin方法．由于简单的加减乘除计算，未涉及到相对复杂的矩阵运 试函数采用了这种独具特色的几何构造方法，因此 自然单元法相比其它无网格法具有简明高效的突出 三角形划分来搜索计算点的自然邻接点，不像 优点[4]，深入开展自然单元法在工程领域中的应用 EFGM等方法，除了需要搜索积分点影响域内的各 研究更具实用价值．目前国外开展这方面的研究代 场点及计算各场点的权重外，还要进行大量的求矩 表学者有Braun、Sambridge、Sukumar、Belikov、阵逆阵的计算，因此求解同一问题时计算时间较 Gonzalez等[3。]，在国内王兆清对自然单元法的进展EFGM等方法有显著的减少，在计算效率方面具有 进行了总结、朱怀球等亦对自然单元法进行了相应 突出的优越性． 的研究[8-lo]． 本文针对厚板弯曲问题，在自然临接点Laplace 本文将自然单元法用于分析地基板的受弯计算 插值形函数的基础上，求出了其导函数，并基于 问题，属于板问题的无网格法受力分析研究领域，在 该领域，国内外已有一些学者开展了研究，如国外 曲挠度的无网格自然单元法求解控制方程，编制了 Krysl和Belytschko[11]发展了一种无单元Galerkin相应的程序进行计算，通过算例的自然单元法解与 法(EFGM)并对薄板的静力问题进行了分析．Ouat—文献=18]级数解的比较，表明了本文方法用于厚板问 ouati和Johnson[123也采用了EFGM对Kirchhoff题受力分析的可行性和有效性． G X 板进行了分析．LiuR和ChenL[”3也采用了无 同时，由于Mindlin板理论同时亦适用于薄板 网格法对薄板的静力和动力问