弹性力学的重构核粒子边界无单元法.docVIP

弹性力学的重构核粒子边界无单元法第55卷第7期2006年7月 1000--329012006155(07)13215--08 物理 ACTAPHYSICASINICA Vo1.55,No.7,July,2006 ④2006Chin.Phys.Soc. 弹性力学的重构核粒子边界无单元法* 秦义校程玉民 1)(上海大学上海市应用数学和力学研究所,上海200072) 2)(太原科技大学,太原030024) (2006年1月17日收到;2006年2月28日收到修改稿) 将重构核粒子法(RKPM)和边界积分方程方法结合,提出了一种新的边界积分方程无网格方法——重构核粒 子边界无单元法(RKP.BEFM).对弹性力学问题,推导了其重构核粒子边界无单元法的公式,研究其数值积分方案, 建立了重构核粒子边界无单元法离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位移和应力积分 公式.重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点的精确值,所以本方法具有 更高的精度.最后给出了数值算例,验证了本方法的有效性和正确性. 关键词:重构核粒子法,弹性力学,边界无单元法 PACC:0200,4630C 1.引言 起源于1977年的光滑粒子法,标志着无网格方 法的产生….1992年以来,无网格方法迅速发展,其 研究和应用领域不断拓展.目前,无网格方法形成形 函数的方法主要有光滑粒子法,移动最小二乘法,单 位分解法,径向基函数法,重构核粒子法和复变量移 动最tl~Z.乘法等.与有限元法一样,无网格方法一般 采用变分原理或加权残数法建立未知量求解方 程引,而局部边界积分方程法是另外一种形式的 全域无网格方法,其在整个求解域的每个离散点为 中心的子域构造形函数,在每个子域建立局部边界 积分方程,最终建立整体域节点未知量求解方 程.1995年,I_ju提出了重构核粒子法 (reproducingkernelparticlemethod,RKPM)…’,该方 法在光滑粒子法中引入核函数修正项,使理论上精 确重构有限域的近似函数成为可能.基于重构核粒 子法建立的无网格方法已广泛应用于力学分析的各 领域,成为无网格方法的重要分支. 1997年,Mukherjee发表了在求解域的边界配置 节点,采用边界积分方程形成求解方程组的边界点 法(boundarynodemethod,BNM),这是边界积分方程 的无网格方法的开端¨.到目前为止形成的边界 积分方程无网格方法还有边界无单元法(boundary element-freemethod,BEFM)¨’J,杂交边界点法 (hybridboundarynodemethod,HBNM)¨,边界点插值 法(boundarypointinterg,olationmethod,BPIM)和 边界云团法(boundarycloudmethod,BCM)加等,其中 边界无单元法作为边界积分方程直接解法,避免了 处理本质边界条件的困难. 本文采用重构核粒子法构造形函数,提出了重 构核粒子边界无单元法(reproducingkernelparticle boundaryelementfreemethod,RKP.BEFM).对弹性力 学问题,推导了其重构核粒子边界无单元法的公式, 研究其数值积分方案,建立了重构核粒子边界无单 元法离散化边界积分方程,并给出了重构核粒子边 界无单元法的内点位移和应力积分公式.重构核粒 子法形成的形函数具有重构核函数的光滑性,且能 再现多项式在插值点的精确值,所以本文方法具有 更高的精度.作为直接方法,重构核粒子边界无单元 法既具有边界元法半解析,精度高和降维等优点,又 具有无网格方法的优点,也可方便地施加边界条件. 最后通过数值算例分析,验证了本文方法的有效性 和正确性. *国家自然科学基金(批准号和上海市重点学科建设项目(批准号:Y0103)资助的课题 ‘E—mail:ymeheng@sh163.net;ymcheng@staff.shu.edu.cn 物理55卷 2.重构核粒子法的形函数 2.1.重构核粒子法的核函数 重构核粒子法中,未知函数u()的近似函数 u.()可表示为 u.()=I.(;—s)u(s)ds,(1) 其中.(;—s)为改进的核函数 . (;—s)=C(;—s).(—s),(2) 也称为重构核,.(—s)为权函数,对有限域问题 一 般写成 — s)=(),(3) C(;—s)称作修正函数 C(x;—s)=∑6()(—s)’ 毫//(—s)西(),(4) 这里 //(—s):[1,—s,(—s),…,(—s)], (5) 西()=[b.(),b(),…,b()],(6) 其中b()是满足重构条件的函数. 把(1)式用Taylor级数展开,