瞬态热传导问题的无网格局部PetrovGalerkin法.PDF

第２７卷 第２期 桂 林 工 学 院 学 报 Ｖｏｌ２７Ｎｏ２ ２００７年５月 ＪｏｕｒｎａｌｏｆＧｕｉｌｉｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ Ｍａｙ　２００７  文章编号：１００６－５４４Ｘ（２００７）０２－０２９４－０４ 瞬态热传导问题的无网格局部 Ｐｅｔｒｏｖ－Ｇａｌｅｒｋｉｎ法 １，２ ３ １，４ ４ 王虎奇 ，陈莘莘 ，王晓峰 ，刘光焰 （１广西大学 土木建筑工程学院，南宁　５３０００４；２广西工学院机械工程系，广西 柳州　５４５００６； ３湖南工业大学 包装与印刷学院，湖南 株洲　４１２００８；４桂林工学院土木工程系，广西 桂林　５４１００４） 摘　要：基于加权余量法和采用移动最小二乘近似函数作为试函数，研究了无网格局部Ｐｅｔｒｏｖ－ Ｇａｌｅｒｋｉｎ方法在瞬态热传导问题中的应用．阐述了该方法应用于瞬态热传导问题的过程和基本理 论，并编制了相应的计算程序进行计算．算例表明该方法用来求解瞬态热传导问题是有效的． 关键词：无网格法；局部Ｐｅｔｒｏｖ－Ｇａｌｅｒｋｉｎ法；移动最小二乘近似函数；瞬态热传导问题 中图分类号：Ｏ２４１８２；Ｏ４１４１　　 　　　　　　文献标志码：Ａ 场．最后通过算例说明无网格局部 Ｐｅｔｒｏｖ－Ｇａｌｅｒ ０　引言 ｋｉｎ法在瞬态热传导问题中的应用是完全可行的． 计算力学的发展面临着许多难以处理的问题， 在此基础上，可将无网格局部 Ｐｅｔｒｏｖ－Ｇａｌｅｒｋｉｎ法 如具有非常大变形的冲压成型问题、裂纹动态扩 用于瞬态温度应力的计算． 展问题和流固耦合等问题．近年来兴起的无网格 １　移动最小二乘法 法由于只需节点而不需单元或网格信息，前处理 特别简单，能够很方便地用于裂纹开裂、非线性 考虑点ｘ的邻域（子域） ，它位于全域 内， Ω Ω ｘ 大变形等问题的数值计算，具有许多传统的数值 为了近似函数ｕ在子域 的分布，在有限个随机分 Ω ｘ 方法无法比拟的优点，从而受到了国际计算力学 布的节点ｘ（ｉ＝１，２，…，ｎ）上，函数 ｕ的近似式 ｉ 界的高度重视． ｈ ｕ（ｘ）， ｘ ，可以定义为  ∈Ω ｘ ［１，２］ 无网格局部 Ｐｅｔｒｏｖ－Ｇａｌｅｒｋｉｎ法 是一种真 ｍ ｈ Ｔ 　　ｕ（ｘ）＝ ｐ（ｘ）ａ（ｘ）＝ｐ（ｘ）ａ（ｘ）． （１） 正的无