原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第25卷第4期 计算力学学报 V01．25。No．4
Chinese of Mechanics 2008
2008年8月 JournalComputational August
自然单元法在Winkler地基薄板计算中的应用
曾祥勇“， 朱爱军2， 邓安福1
(1．重庆大学土木工程学院．重庆400044；2．贵州大学土木建筑工程学院．贵阳550003)
摘要：自然单元法是一种基于Voronoi图及Delaunay三角形剖分图，以自然邻接点插值为试函数的一种无网
格数值方法．本文以目前该方法中自然邻接点的Laplace插值形函数为基础，求出了其一阶及二阶导函数，建立
了winkler地基上薄板弯曲挠度的自然单元法求解控制方程，并编制了相应的计算程序。通过算例分析表明了
本文方法的可行性和有效性。
关键词：自然单元法lLaplace插值函数；自然邻接点；薄板弯曲；无网格法；Winkler地基
中图分类号：TU470 文献标识码：A
对要求域内导数连续的C1问题(例如薄板弯曲问
1 引 言
题)，在自然邻接点Laplace插值形函数的基础上，
目前求解计算力学问题的经典数值方法主要 求出了其一阶及二阶导函数，建立了Winkler弹性
有有限元法、有限差分法和边界元法等。由于有限 地基上薄板弯曲挠度的无网格自然单元法求解控
元法、有限差分法依赖于求解区域的网格划分，使 制方程，并编制了相应的程序进行计算，通过算例
其在裂纹扩展、大变形和移动相边界等问题的处理
上，必须进行网格重构，耗费大量的人力和物力。
为克服这一困难，一些新的数值方法被研究出来， 法)解的比较，表明了本文方法的可行性和有效性。
如无网格Galerkin法等各种无网格方法和主要应
2 自然单元法的插值形函数
用于岩石力学的流形方法[1矗]。自然单元法即为
目前众多无网格法(亦称无单元法)中的一种。 2．1 自然单元法的基本理论[．1
在二维欧几里德空间R2上，考虑一系列离散
Braun与Sambridge1995年在Nature上发表文
Element
章，第一次提出了自然单元法(Natural
Voronoi图)由二维平面上一系列闭合凸面多边形
Method)的概念并将这种新型的数值方法用于偏
微分方程的求解[3]。自然单元法基于给定节点的 或无界的子区域组成，每一个节点九，对应的子区
域T，内，任意一个属于1’，区域的点至n-(最近的
Voronoi图及Delaunay三角形剖分图，是一种以
邻接点)的距离均小于该点至其他节点九』(J≠J)
自然邻接点插值函数作为试函数的Galerkin方
法。目前国外开展该方法研究的代表学者有 的距离。以上定义可用式(1)表示为
Braun，Sambridge，N．Sukumar和Belikov等E3-53，
国内王兆清对自然单元法的进展进行了总结，蔡永 (1)
昌和朱怀球等亦对自然单元法进行了研究[1—7。。 由图1可以看出，构成Voronoi图的各子区域
目前关于自然单元法的研究大多为将其应用 丁，均由不同节点间的垂直平分线分割而成。
于只要求场函数本身连续的C0问题中[3可]，本文针
文档评论（0）