首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到16条相似文献,搜索用时 78 毫秒
1.
自然单元法研究进展   总被引:15,自引:2,他引:13  
王兆清  冯伟 《力学进展》2004,34(4):437-445
自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.  相似文献   

2.
Delaunay多边形单元的有理函数插值格式   总被引:11,自引:0,他引:11  
王兆清  冯伟 《力学季刊》2004,25(3):375-381
本文提出了基于Delaunay多边形化的多边形单元有理函数插值格式。给出了Delaunay多边形化的概念和Delaunay多边形单元有理函数插值形函数的计算表达式。与Delaunay三角化网格不同,Delaunay多边形化网格形成对区域的唯一剖分。Delaunay多边形单元有理函数插值是以Delaunay多边形的顶点作为插值点,构造的有理函数形式插值。Delaunay多边形单元有理函数插值克服了有限元方法中难以构造边数大于4单元多项式形式位移插值的困难。有理函数插值形函数在多边形单元的内部是无穷次光滑的,在多边形的边界上是线性的。在三角形单元和矩形单元上,有理函数插值分别等价于有限元的三角形面积坐标插值和四边形双线性插值。给出了Delaunay多边形有理函数插值在圆域温度分布插值近似中的两个算例。  相似文献   

3.
高度不规则网格多边形单元的有理函数插值格式   总被引:4,自引:0,他引:4  
借鉴自然邻点插值法,提出了基于高度不规则网格多边形单元的有理函数插值格式一多边形有理函数插值.给出了多边形有理函数插值形函数的计算表达式.该插值格式以多边形的顶点作为插值点,插值形函数为有理函数形式,克服了传统有限方法中构造边数大于4单元多项式形式位移插值的困难.  相似文献   

4.
江涛  章青 《力学与实践》2008,30(4):79-83
基于Lasserre体积算法推导了两种插值方案下自然单元法形函数及其导数 的具体计算方法,特别是对计算点处于某些特殊位置时可能造成计算失败的原因和处理方法 进行了较为深入的研究. 算例结果验证了Sibson与non-Sibson插值形函数在三角形外接圆 的圆周上具有不同的连续性,自然单元法形函数在凸区域的边界结点间是线性变化的,因而 可以方便地施加本质边界条件.  相似文献   

5.
自然单元法在Winkler地基薄板计算中的应用   总被引:3,自引:1,他引:2  
自然单元法是一种基于Voronoi图及Delaunay三角形剖分图,以自然邻接点插值为试函数的一种无网格数值方法.本文以目前该方法中自然邻接点的Laplace插值形函数为基础,求出了其一阶及二阶导函数,建立了Winkler地基上薄板弯曲挠度的自然单元法求解控制方程,并编制了相应的计算程序,通过算例分析表明了本文方法的可行性和有效性.  相似文献   

6.
Winkler地基上厚板分析的自然单元法   总被引:2,自引:0,他引:2  
以自然单元法中自然邻接点的Laplace插值形函数为基础,基于Mindlin厚板理论,建立了Winkler地基上厚板弯曲挠度的自然单元法求解控制方程,并进行了相应的程序实现,最后通过算例分析,表明了该文方法的可行性和有效性.  相似文献   

7.
基于自然单元法的极限上限分析   总被引:2,自引:0,他引:2  
自然单元法是一种基于离散点集的Voronoi图和Delaunay三角化几何信息,以自然邻近插值为试函数的新型数值方法.相对于一般无网格法中常采用的移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,更不需要任何人为的参数,可以提高计算效率.采用该方法构造的形函数满足Delta函数的性质,可以像有限元一样准确地施加边界条件,可以方便处理场函数及其导数的不连续性的问题.论文将自然单元法应用到极限上限分析中,编制了相应的计算程序,通过极限分析的几个经典算例进行了验证,同时采用类似于分片应力磨平的方式,编制相应的磨平程序,由计算点上的塑性耗散功外推得到了节点上的塑性耗散功的值,从而画出了极限状态下结构的塑性耗散功的分布云图.计算结果表明采用自然单元法求解极限上限分析具有稳定性好,精度高,收敛快等优点.  相似文献   

8.
基于局部搜索算法的自然邻接点方法   总被引:9,自引:0,他引:9  
蔡永昌  朱合华 《力学学报》2004,36(5):623-628
自然邻接点方法(NNM)采用自然邻接点形函数进行插值,其插值形函数具有严格定义,且与 有限元形函数一样形式简洁、性能优良,因而避免了EFG法里难以准确施加位移边界条件和 材料不连续条件等诸多主要困难. 但是从形式上看自然邻接点方法仍然属于有网格的方法, 其研究和应用受到了较大的限制. 为了克服这个缺点,对于任意给定的数值积分点,提出了 一种基于局部搜索自然邻接点的寻找算法对NNM进行改进. 改进后的NNM与无单元伽辽金法 (EFG)的插值和求解过程类似,兼具有EFG的真正无网格特性及NNM的便于处理边界和材料 不连续条件等优点. 所得计算结果表明,改进后的NNM的计算精度和计算时间与NNM相当, 是一种比较理想的数值求解方法.  相似文献   

9.
三维弹塑性自然单元法算法实现   总被引:1,自引:0,他引:1  
自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然相邻插值(C∞)的伽辽金法。该方法计算精度与四边形或六面体单元有限元法相当,自然相邻插值函数比其他无网格法插值函数的计算速度快。由于自然相邻插值在凸域的边界上的相邻点之间是严格线性的,所以自然单元法在边界面的处理也相当简单。本文研究了在自然单元法中采用Von.Mises,Mohr-Coulomb和Drucker-Prager屈服准则解决三维弹塑性问题,并编制了相应计算程序,最后通过算例验证算法的正确性。  相似文献   

10.
多边形有限元研究进展   总被引:10,自引:0,他引:10  
王兆清 《力学进展》2006,36(3):344-353
有限元法是数值求解偏微分方程边值问题的重要方法, 采用 不规则多边形单元网格, 可以方便有效地模拟材料的力学性能, 又使得区域网格剖分变得灵 活方便. 特别是对于复杂的几何形状, 多边形单元网格具有更大的优势. 本文对国内外有关 多边形有限元法的最新进展作了初步的总结和评述, 主要以基于位移法的多边形有限元为主. 论述了多边形有限元的发展历史, 给出了多边形单元上的Wachspress插值、Laplace插值和 重心坐标的一些最新研究成果. 与经典有限元法形函数为多项式形式不同, 多边形单元的形 函数为有理函数或者无理函数形式. 多边形单元插值形函数满足线性完备性, 可以再现线性 位移场, 像经典有限元法一样直接施加本质边界条件; 插值函数在多边形的边界上是线性的, 确保不同单元间的自动协调. 不同单元的插值形函数表达公式形式统一, 方便混合单元网格 计算的程序编写. 提出了多边形有限元法今后需要研究的问题.  相似文献   

11.
利用基于局部移动Kriging插值无网格法对层合板自由振动进行了数值分析,基于一阶剪切层合理论导出了层合板振动的控制方程和边界条件,进一步得到了自由振动的离散化特征方程。由于Kriging插值函数具有Kronecker delta函数性质,可以直接施加本质边界条件。通过本文给出的方法,对不同边界条件、不同跨厚比、不同材料参数和铺设角度的层合板的振动频率进行了计算,均得到满意结果。最后用该方法对层合板的铺设角度进行优化设计,得到了与已有文献完全一致的优化结果。数值结果充分表明了无网格Kriging方法分析层合板自由振动问题的有效性和高精确度。  相似文献   

12.
利用模态综合法分析车辆与桥梁之间的相互作用时,合理地构造桥梁的插值振型函数可以大幅提高计算精度.其中,分段三次Hermite插值函数和三次样条插值函数较为常用.为研究二者的异同,以简支梁桥为例分别采用这两种插值函数构造结构梁单元模型的一维插值振型函数和板单元模型的二维插值振型函数.基于以上两类插值振型函数,分析单自由度簧上质量匀速过桥时,桥梁的跨中位移、跨中梁底正应力和轮-桥接触力时程响应.结果表明:无论是一维问题还是二维问题,由三次样条插值法构造的插值振型函数与结构的实际振型较为吻合,计算结果具有较高的收敛性和精度.而要达到相同的精度,分段三次Hermite插值法则须加密单元网格,但其误差仅存在于独立网格内,不会累积放大.  相似文献   

13.
范纪华  章定国 《力学学报》2016,48(4):843-856
研究了基于变形场不同离散方法的柔性机器人动力学建模和仿真问题.针对多杆空间链式柔性机器人系统,采用假设模态法、有限元法、Bezier插值方法和B样条插值方法对柔性杆变形场进行描述,构造统一形式,运用Lagrange方法,结合4×4齐次变换矩阵,在计入柔性杆横向弯曲变形引起的纵向缩短的情况下,推导得到多杆空间柔性机器人动力学方程,并编制基于4种变形场不同离散方法的多杆空间链式柔性机器人仿真软件.通过仿真算例对柔性机器人系统的动力学问题进行研究.仿真结果表明:有限元法的计算效率较低;假设模态法在处理较大变形问题时其精度低于Bezier插值方法和B样条插值方法的精度;作为新的变形体离散方法,Bezier插值方法和B样条插值方法可以有效地描述柔性杆的变形场,并能运用到多杆空间柔性机器人动力学建模中.  相似文献   

14.
A local basis algorithm for searching natural neighbours in Natural Element Method (NEM) is presented for solving the elasticity problems in this paper. Comparison with the global sweep algorithm used in natural element method or Natural Neighbour Method (NNM) for searching natural neighbours, the proposed algorithm is more expedient and convenient in the constructions and computation of natural neighbour interpolations. In the proposed NEM based on local search, the Laplace (non-sibson) interpolations are constructed with respect to the natural neighbour nodes of the given point which have been locally defined. The shape functions from the Laplace approximations have the delta function property and the Laplace interpolants are strictly linear between adjacent nodes, which facilitate imposition of essential boundary conditions and treatment of material discontinuity with ease as it is in the conventional finite element method. The Laplace interpolants derived from the local algorithm and the global algorithm in NEM are identical because of the uniqueness of the Voronoi diagram. Numerical results and convergence studies also show that the present NEM based on local search algorithm possesses the same accuracy and rate of convergence as they are in previous NEM.  相似文献   

15.
基于非结构混合网格的N-S方程求解器和结构柔度影响系数法,发展了一种考虑气动、结构非线性的基于RBF插值技术CFD/CSD耦合分析方法,适用于解决现代大展弦比飞机的非线性静气动弹性问题。该方法采用时间相关法(即求解非定常方程组,用长时间的渐近解趋于定常状态)求解静气弹分析时的定常流动。考虑大展弦比飞机结构变形问题为大变形小应力问题,在利用柔度系数法求解结构方程时,假设每次求解结构方程时应力与应变为线性关系,整体静气弹分析过程为非线性关系,因此每次求解结构方程时要更新柔度影响系数矩阵。在非定常N-S方程每求解一个时间步耦合一次结构有限元分析,由于结构有限元分析的时间相对于气动分析时间是很短的,所以这种方法实际上近似使用了一次求解非定常气动力的时间完成了整个静气动弹性分析的过程。对于气动网格与结构有限元网格不一致性,本文采用径向基函数(RBF)插值方法中的TPS方法进行结构弹性变形和气动载荷插值,采用虚功原理完成气动载荷数据交换。为了节省气弹分析时间,采用动网格方法对气动网格进行更新,本文基于RBF插值方法发展一种适用于混合网格(四面体、三棱柱、金字塔和六面体)变形的动网格方法,可以保证附面层网格的质量与分布从而准确模拟其流动。利用该方法对M6机翼、DLR-F6翼身组合体和某大型客机机翼进行了静气动弹性特性分析,结果验证了本文开发的非线性CFD/CSD耦合分析方法的可行性、精确性和高效性。  相似文献   

16.
粗糙表面测量轮廓的分形插值模拟   总被引:4,自引:1,他引:3  
引入粗糙表面的分形插值理论,提出用分形插值函数模拟粗糙表面轮廓的方法。并以车削,磨削及磨损表面为研究对象,对粗糙表面进行了分形插值模拟,系统地研究了影响其模拟模拟效果的因素和规律。发现可以采用分形插值理论来模拟或表征具有分形特征的粗糙表面。  相似文献   

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

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