基于自然单元法的极限上限分析.pdfVIP

第33卷第1期 固体力学学报 V01．33No．1 CHINESE OFSOLIDMECHANICS 2012年2月 JOURNAL 2012 February 基于自然单元法的极限上限分析‘ 周书涛1 刘应华p+ 陈莘莘2 (1清华大学航天航空学院工程力学系，北京，100084) (2湖南工业大学土木工程学院，株洲，412008) 摘 要 自然单元法是一种基于离散点集的Voronoi图和Delaunay三角化几何信息，以自然邻近插值为试函 数的新型数值方法．相对于一般无网格法中常采用的移动最小二乘近似而言，自然邻近插值不涉及到复杂的矩阵 求逆运算，更不需要任何人为的参数，可以提高计算效率．采用该方法构造的形函数满足Delta函数的性质，可以像 有限元一样准确地施加边界条件，可以方便处理场函数及其导数的不连续性的问题．论文将自然单元法应用到极 限上限分析中，编制了相应的计算程序，通过极限分析的几个经典算例进行了验证，同时采用类似于分片应力磨平 的方式，编制相应的磨平程序，由计算点上的塑性耗散功外推得到了节点上的塑性耗散功的值，从而画出了极限状 态下结构的塑性耗散功的分布云图．计算结果表明采用自然单元法求解极限上限分析具有稳定性好，精度高，收敛 快等优点． 关键词 自然单元法，自然邻近插值，极限上限分析，Sibson插值 更不需要任何人为的参数，使得该方法的计算效率 0 引言 比一般的无网格法有较大的提高．本文根据极限分 析的上限定理，给出了基于自然单元法的极限分析 在结构的设计与安全评估中，结构的塑性极限 上限求解格式，最后通过极限分析的几个经典算例 分析是一种评估结构承载能力的有效方法之一，通 的计算和对比分析，验证了自然单元法应用于极限 过确定各类工程结构的极限载荷，可以为工程设计 分析中的有效性和合理性． 和安全评定提供准确可靠的理论依据．目前极限分 析方面的研究都是围绕着发展有效的极限分析数值 1 自然单元法 方法而展开，传统极限分析的数值方法都是基于网 格的，如有限元法，边界元法等[1。4]．近年来兴起的无 网格方法受到了国际计算力学界的广泛重视．虽然 年首先提出，并用自然邻近插值作为试函数来求解 该方法目前在严格的数学论证，计算效率和实际工 偏微分方程．自然临近插值是基于离散节点的 程应用方面还不能与成熟的有限元相比，但是该方 法是采用基于点的近似，可以部分或彻底消除网格， 方案，按照插值基函数的不同，常见的应用于自然单 从而能有效的避免对于依赖网格而带来的种种困 难．采用不同的加权残量法和近似函数，无网格法可 以分为很多种[53．自然单元法邸1是近十多年来发展 等不同的方案．这三种不同插值方案各有利弊， 起来的求解偏微分方程的一种数值方法，这是一种 基于离散节点的Voronoi图和Delaunay三角化，以 只在某些特殊的变分原理中采用；Sibson插值方案 自然邻近插值作为试函数，以Galerkin法作为加权 残量法而得到的新型数值方法．该方法以Delaunay 三角形或四面体为背景网格进行数值积分，由于该 形函数在节点处是Co连续，在除节点以外的区域是 插值方法构造简单，不涉及到复杂的矩阵求逆运算，