收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

接触碰撞问题的非线性无网格方法及在金属成型分析中的应用

SIDIBE Kalilou(西迪贝 卡里鲁)  
【摘要】:第一章 引言 无网格方法是近年来兴起的一种新的数值方法。在无网格方法中,研究的问 题域由一系列的离散点组成,在进行仿真计算时,只需节点离散数据,而不 需要单元信息,只采用基于节点的近似,可以完全抛开网格,克服了有限元 法由于有网格的存在而难以方便处理大变形畸变问题、裂纹扩展问题等缺 陷。 在无网格方法中,典型的、具有代表性的有以下几种,它们是:(1)光滑质 点流体动力学法(Smooth Particle Hydrodynamics,SPH);(2)模糊单元法(Diffuse Element Method,DEM);(3)无网格伽辽金法(Element-Free Galerkin(EFG)Method); (4)无网格局部伽辽金法(Meshless Local Petrov-Galerkin(MLPG)Method,);(5)再 生核质点法(Reproducing Kernel Particle Method,RKPM);(6)单位分解法(Partition of Unity Method,PUM);(7)点插值法PIM(Point Interpolation Method)等等。 再生核质点RKPM法是Liu等人基于再生核(reproducing kernel)思想及小波 变换理论提出的一种无网格方法。采用窗口函数和傅立叶变换建立新的形函 数,由于窗口函数可以平移、缩放,可以应用于弹性、塑性和动力问题。本文 正是利用RKPM方法分析接触碰撞问题。 接触碰撞仿真问题的研究已经有很长的历史了。早期的接触体都是假定一 个为刚体另一个为简单的弹性体,分析主要集中于总接触力的计算。牛顿第三 定律及库伦定律用于计算接触界面的相互作用力。随着数值科学和工程应用的 发展,更精确的接触算法也应运而生。本文提出一种新的从点到面的接触算法 来处理无网格接触问题。接触力的计算用的是罚函数法。 积分格式在数值计算中非常重要,它关系到计算精度和CPU耗时。无网格 方法最大的缺点是计算时间长。本文提出一种新的积分法则:降阶高斯积分法 则,可以大大降低计算时间。 博士学位论文:摘要 在数值仿真中,公式及算法的验证也是非常重要的。本文通过提出一个具 有一般意义的E七n chn址甘ktest算法来验证本文中运用的降阶高斯积分方案及接 触算法在实际应用中的正确性及准确性,二维和三维问题都得到成功的验证。 第2章RRrM法的基本原理 RKpM方法是再生核质点法(斑犷oduc吨K治n岭1 Part记leM剧助d)的缩写。 同SPH方法一样,RKPM方法对变量域的近似也是采用积分的形式 ur(,)一工u(x)砚(x一y)方。不同于sPH法的是:吩M方法~‘函数中多出 了一个修正函数C(x,x一y),即,砚(l一y)=C仁,x一y)丸(x一y),修正函数是 通过施加再生条件获得的,可以精确的再生多项式。电(x一y)即K淤n姆l函数。 由于在构造形函数时施加了再生条件,RKPM法满足一致连续性,这一点是 SPH法缺乏的。因此RK卫M法可以看成是利用修正函数恢复连续性的SPH法。 类似其他无网格法,RKPM形函数利用离散质点建立插值函数。RKPM方法中 的原理及公式的推导在Liu现K,C址n工s,Lis.,Belytscbko等的论文中做了较 系统的论叙。在他们的工作的基础上,本章对K即叱l函数汽(l一y)及其支持 域,待定系数br(x)的计算等内容都做了详细的介绍和推导。在本文中,x表示 一点在t时刻的位置矢量,X表示材料粒子的初始位置矢量。 2.1 Ke.el函数及其支持域 问题域由一系列离散点卜:,xZ,…,x,}组成,x,是粒子I的位置矢量,N’P 是粒子的总数量。无网格法的一个共同的特征就是它们都含有权函数,也就是 K治n粗l函数。任何一种权函数都有一个集中域,在这个域里的函数值是非零 值,而集中域以外的函数值为零。这个集中域即前面所讲的支持域或影响域。 这里所讲的权函数在小波理论里也称窗函数。 120 通常使用的支持域为圆形(球形)或矩形(立方形),如下图所示 (a)圆形 助矩形 图2.1支持域 上图中的支持域。;为问题域。的一个子域。 根据Monagh川(1 982),权函数必须满足以下条件: a .bC‘d (2(2(2(2 1.牙(s)0 2.平(s)==. 在。:中(正值性) 不在。,中(紧凑性) J甲(s)‘一1,(合一性) O 不(s)为单调递减的函数,(递减性) 牙(s)峥6(s),当峥0时 (2 .le) 8(s)即de、函数,:={巨二兰业 a即支持域的宽度。 常见的权函数有: 指数权函数二(s)一介一(’/a”,_黔‘1 口,当sl (2 .2a) 一{ 号一‘·’+‘·’,当·‘全 警一‘一‘一警·’,当合二‘, (2 .2b) 当s1 121 博士学位论文:摘要 四次方样条权函数二(·)·}i一‘·’+8·’一,·‘,肖·‘, LO,兰sl (2 .Ze) 2.2 RKPM中的Keruel近似形式 SpH方法是最早基于K图凹1近似的无网格方法。u(l)的近似形式是通过对 K日比吧1函数的积分得到的,即: u‘(x)一仲。(x一y)u份)方 口 ‘,____、1_,,}!x一y}}、止_,二‘__,。。__小二‘,~‘ ,七气x一y)二一尸又一J刀万J二1毋i幽资汉,W林F乙苗工、甲划花K胭如哭。 (2 .3) aa (2.3减的离散形式如下 召尸 u‘(x)二艺.。(x一x,)u(x,)成 I二t (2,4) 方程俘Ic)的积分形式能够保证零阶连续,然而恤wK氰1卯sc)和 氏加毗如氰1更场脂出它并不能保证2.4式的离散形式连续,因此也无法保证问 题的收敛性。为了解决这个问题,玩等人(l99阮)在K朗.1


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 娄路亮,曾攀,方刚;无网格方法及其在体积成形中的应用[J];塑性工程学报;2001年03期
2 王卫东,赵国群,栾贻国;材料加工数值模拟的新进展——无网格方法[J];材料导报;2002年09期
3 曾清红,卢德唐;含有启动压力梯度的渗流问题及其无网格解法[J];计算力学学报;2005年04期
4 王辉;许君风;赵旭升;;功能梯度梁的无网格分析[J];中原工学院学报;2008年04期
5 苏华,王光远,张学鸿;细长杆与圆孔孔壁接触碰撞问题的动力间隙元法[J];地震工程与工程振动;1996年01期
6 樊成;栾茂田;黎勇;杨庆;;Kriging插值无网格方法及其在力学边值问题中的应用[J];岩石力学与工程学报;2007年01期
7 韩光;樊成;王大国;;具有过点插值性质的无网格法及其应用[J];辽宁工程技术大学学报;2007年S2期
8 蔡山;李治平;张明友;;使用无网格方法计算有动边界的试井典型曲线研究[J];油气井测试;2008年01期
9 吴欣;赵国群;徐淑波;管延锦;;刚塑性无网格伽辽金方法近似方案研究[J];塑性工程学报;2008年03期
10 赵群;;基于无网格方法的钢结构构件优化模拟分析[J];武汉工程大学学报;2009年05期
11 熊渊博;龙述尧;李光耀;;层合板分析的无网格局部Petrov-Galerkin方法[J];复合材料学报;2005年06期
12 吴欣;赵国群;路平;管延锦;;基于一致径向基函数的金属塑性成形过程无网格方法分析[J];机械工程学报;2008年02期
13 樊成;王大国;;Kriging插值无网格法形函数性质[J];辽宁工程技术大学学报(自然科学版);2009年04期
14 赵国群;吴欣;王卫东;;刚塑性无网格方法中体积闭锁问题的缓解算法[J];机械工程学报;2006年06期
15 吴欣;赵国群;王卫东;;任意边界形状的二维刚塑性成形过程的无网格分析方法[J];计算力学学报;2007年03期
16 李雷;周毅英;谢水生;贺睿;;应变梯度形变理论的无网格数值方法研究[J];塑性工程学报;2009年04期
17 赵国群;吴欣;管延锦;;刚塑性无网格方法中精度影响因素的研究[J];塑性工程学报;2006年03期
18 张小华;欧阳洁;孔倩;;聚合物流动的多尺度模拟[J];化工学报;2007年08期
19 陈海龙;陈摇娟;初文华;位摇莎;;无网格方法在爆轰数值模拟中的应用[J];哈尔滨工业大学学报;2010年11期
20 曾亿山,卢德唐,晏忠良,曾清红;无网格方法在油藏工程中的应用初探[J];油气井测试;2003年05期
中国重要会议论文全文数据库 前10条
1 吴兴;李志鹏;王昌生;;流固耦合问题的非边界适应方法及无网格方法的研究进展[A];第二十三届全国水动力学研讨会暨第十届全国水动力学学术会议文集[C];2011年
2 陈勇;黄建明;陆明万;;地震波传播计算中的无网格方法[A];1999年中国地球物理学会年刊——中国地球物理学会第十五届年会论文集[C];1999年
3 倪国喜;王瑞利;林忠;;无网格方法中的粒子初始均匀分布[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
4 刘高联;李孝伟;;基于当地笛卡儿架构的无网格方法[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
5 雷雨冰;袁亚雄;;无网格数值方法研究及其在流体力学中的应用[A];加入WTO和中国科技与可持续发展——挑战与机遇、责任和对策(下册)[C];2002年
6 郑云英;杨骁;;饱和黏弹性多孔介质动力学问题的无网格方法[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
7 王瑞利;倪国喜;林忠;;关于无网格方法中撒点算法[A];中国工程物理研究院科技年报(2005)[C];2005年
8 潘小飞;张雄;陆明万;;配点型无网格方法的研究[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
9 王东东;寇勇;金科;王省哲;;基于无网格方法的铁磁梁式板磁弹性屈曲与弯曲分析[A];中国力学学会学术大会'2009论文摘要集[C];2009年
10 刘方玉;蔡山;雷霆;;使用无网格方法计算有动边界的渗流问题[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
中国博士学位论文全文数据库 前10条
1 SIDIBE Kalilou(西迪贝 卡里鲁);接触碰撞问题的非线性无网格方法及在金属成型分析中的应用[D];湖南大学;2005年
2 李坤;基于无网格方法的声学问题数值模拟研究[D];华中科技大学;2011年
3 任红萍;插值型无网格方法研究[D];上海大学;2010年
4 李九红;复变量无网格方法及其应用研究[D];西安理工大学;2004年
5 秦伶俐;无网格方法在结构计算上的改进[D];中国农业大学;2005年
6 蔡永昌;流形方法的理论与应用研究[D];重庆大学;2001年
7 郑兴;SPH方法改进研究及其在自由面流动问题中的应用[D];哈尔滨工程大学;2010年
8 张雷;基于无网格法的连铸结晶器内坯壳凝固与热应力的数值模拟[D];清华大学;2005年
9 熊渊博;Kirchhoff板问题的无网格局部Petrov-Galerkin方法研究[D];湖南大学;2005年
10 戴保东;改进的无网格局部边界积分方程方法研究[D];上海大学;2006年
中国硕士学位论文全文数据库 前10条
1 郑兴;光滑质点流体动力学(SPH)算法研究[D];哈尔滨工程大学;2005年
2 王晓光;弹性及弹塑性问题的无网格局部边界元法[D];浙江工业大学;2003年
3 丁慧芬;对流扩散方程的径向基函数法[D];辽宁师范大学;2008年
4 黄娟;改进广义移动最小二乘近似的无网格方法[D];苏州大学;2007年
5 马修彦;弹性问题无网格局部边界元法及其应用研究[D];浙江工业大学;2003年
6 苏芳;基于局部多顶式近似空间的单位分解方法及其误差分析[D];湘潭大学;2005年
7 孙阳光;自然边界元的无网格方法和拟小波方法研究[D];华中科技大学;2005年
8 刘建军;基于无网格方法的板壳大变形数值模拟关键技术研究[D];山东理工大学;2012年
9 张小华;无网格方法在计算流体力学中的应用研究[D];西北工业大学;2006年
10 李武;无网格数值积分方法的理论研究[D];安徽理工大学;2007年
中国重要报纸全文数据库 前10条
1 田涛;西方人眼中的中国法律(之十)[N];法制日报;2007年
2 本报记者  管宏业;“奔腾”撞出高质量[N];中国商报;2006年
3 王乾;投资7亿改善人居环境[N];黑龙江经济报;2008年
4 记者 李宏见习记者 鲍蔓华;宁波转移支付1.8亿元支持域内欠发达地区建设[N];人民政协报;2008年
5 本报记者 冯竞 吴晋娜;让科技更好融入经济社会发展[N];科技日报;2011年
6 本报记者  盛若蔚;钟志华:中国创造 时不我待[N];人民日报;2006年
7 通讯员 钟吉福;给卫星和飞船建造舒适的“家”[N];科技日报;2002年
8 记者 谢涛;通用汽车中国科研基金成立[N];中国证券报;2000年
9 严华;下一代网络面面观(优势篇)[N];人民邮电;2001年
10 记者 俞安;掀起跑项目上项目建设项目热潮[N];抚顺日报;2007年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978