学位论文 > 优秀研究生学位论文题录展示
无网格法及其在任意形状接触界面金属成形中的应用研究
作 者: 苏建民
导 师: 王耘;胡树根;宋小文
学 校: 浙江大学
专 业: 机械设计及理论
关键词: 伽辽金无网格法 弹塑性 任意形状接触界面 质点积分 投影子域
分类号: TG306
类 型: 硕士论文
年 份: 2010年
下 载: 54次
引 用: 0次
阅 读: 论文下载
内容摘要
金属塑性成形技术在机械制造业中具有重要地位。数值模拟技术提高了塑性成形过程的可控化水平,是优化工艺参数和模具结构的有力工具,对保证工件质量,减少材料消耗和开发时间,降低制造成本有重要的意义。有限元法是应用最为广泛和有效的数值模拟技术。但迄今为止,有限元法模拟金属塑性成形过程时,存在应力场精度达不到工程要求和网格容易发生畸变等诸多瓶颈性问题。这些困难的根本原因在于有限元法中单元网格的局限性。无网格法不需要产生网格,避免了有限元法中单元网格带来的一系列困难。它不仅能保证计算的精度,还能降低计算的难度,因此将无网格法用于分析金属塑性成形过程有着重要的理论和实际意义。考虑到任意形状接触界面是金属塑性成形中十分常见的接触情形,因此本文将无网格法用于分析任意形状接触界面的金属塑性成形问题。计算效率是无网格法的一个关键性制约因素,而积分运算和接触界面处理算法是影响无网格法效率的两大重要原因,因此论文从这两方面入手,着重研究提高无网格法计算效率的相关算法。论文选择基于移动最小二乘近似函数的伽辽金无网格法作为数值分析方法,采用弹塑性材料模型作为计算的材料模型。在任意形状接触界面的处理方面,开展了模具型腔的描述方法、接触对搜索算法、法向接触力和切向摩擦力的计算方法等工作。在此基础上,研究了基于更新拉格朗日描述的弹塑性材料变形过程控制求解方程及其相应矩阵的定义式,并使用基于中心差分格式的动力显式求解算法求解该方程,最终实现了伽辽金无网格法在任意形状接触界面金属塑性成形过程的数值模拟。论文进行的创新性工作有以下两点:1将质点积分法用于任意形状接触界面的金属成形过程分析,在保证精度的前提下实现了高效率的积分运算。质点积分法是一种高效稳定的积分方法,本文基于对任意形状接触界面处理算法的深入研究,将其用于对一个使用球头形凸模挤压块状铜质工件的塑性变形过程进行数值分析。将计算结果与基于ANSYS/LS-DYNA环境的有限元法得到的结果进行比较,说明了本文提出方案的正确性和可行性。2在系统研究主从面法和一体化算法的基础上,提出了一种新的接触对快速搜索算法一投影子域法。将模具型腔表面的主接触块和工件表面的从接触点均沿压力机运动方向投影,采用子域分割思想,用矩形子域分割接触界面投影区域,实现从接触点和主接触块快速匹配判定运算,进而实现测试对和接触对的确定。算例表明子域投影法具有处理速度快、实施方便等特点。
|
全文目录
摘要 4-5 Abstract 5-9 第一章 绪论 9-24 1.1 引言 9-10 1.2 有限元法分析金属塑性成形过程的现状 10-14 1.2.1 板料成形过程有限元模拟及其存在的问题 11-13 1.2.2 体积成形过程有限元模拟及其存在的问题 13-14 1.3 无网格法分析金属塑性成形过程的现状 14-21 1.3.1 配点型无网格法在金属塑性成形中的应用 15-16 1.3.2 伽辽金型无网格法在金属塑性成形中的应用 16-21 1.4 论文主要工作内容 21-24 第二章 伽辽金无网格法研究及其关键技术分析 24-41 2.1 紧支近似函数 24-33 2.1.1 移动最小二乘近似函数的构造 25-27 2.1.2 权函数的选择及影响域大小的确定 27-33 2.1.3 移动最小二乘近似函数的一致性探讨 33 2.2 伽辽金无网格法 33-41 2.2.1 弹性力学问题的伽辽金无网格法 34-35 2.2.2 数值积分方法分析及其选用 35-38 2.2.3 几种本质边界条件的处理方法 38-40 2.2.4 伽辽金型无网格法程序流程图 40-41 第三章 有限应变弹塑性材料的无网格法研究 41-51 3.1 物体构形和运动的描述方法 41-43 3.2 弹塑性材料模型的本构关系 43-47 3.2.1 三种应变张量及其应用 43-44 3.2.2 速率梯度张量及其应用 44-45 3.2.3 三种应力张量及其应用 45-46 3.2.4 有限应变弹塑性材料的本构关系方程 46-47 3.3 弹塑性问题无网格法方程的构建 47-51 3.3.1 动力虚功率方程 48-49 3.3.2 无网格控制方程 49 3.3.3 动力显式求解算法的应用 49-50 3.3.4 显式算法的稳定性分析 50 3.3.5 集中质量矩阵和集中阻尼矩阵的应用 50-51 第四章 质点积分法分析复杂界面的接触问题 51-70 4.1 接触问题的一般描述 51-53 4.2 模具型腔描述方法 53-54 4.3 两种接触搜索算法 54-60 4.3.1 主从面法 54-58 4.3.2 一体化算法 58-60 4.4 接触力的计算方法 60-62 4.5 质点积分法 62-64 4.6 数值算例及结果分析 64-70 第五章 投影子域法及其应用研究 70-87 5.1 投影矩形区域及其子域的定义 71-72 5.2 接触点的整数坐标及其子域编号 72-74 5.3 主接触块的相交子域及其编号计算 74-76 5.4 测试对和接触对的确定 76-77 5.5 算法实现流程及分析 77-80 5.6 数值算例及结果分析 80-87 第六章 结论和展望 87-90 6.1 本文主要工作 87-88 6.2 本文主要创新点和结论 88-89 6.3 不足和展望 89-90 附录 算法中的类及其实现程序 90-101 参考文献 101-106 本人在攻读硕士学位期间的科研工作及成果 106-107 致谢 107
|
相似论文
- 带填充墙框架结构非线性有限元分析,TU323.5
- 土的三剪统一帽子模型研究,TU411.7
- 圆柱形粗糙峰的塑性接触模型研究,TG84
- 混凝土防渗墙加固的土石坝有限元分析研究,TV543.8
- 金属粉末的粘弹塑性本构方程的研究,O175
- 基于非局部摩擦的钢丝滚道球轴承接触仿真与实验研究,TH133.3
- 地震作用下钢筋混凝土结构弹塑性分析,TU375
- 沉入式大圆筒防波堤稳定性的三维弹塑性静动力数值分析,U656.2
- 钢管混凝土框架—钢板组合剪力墙结构体系抗震性能分析,TU398.9
- 相邻结构体系的主动和半主动控制研究,TU352.11
- 人字形偏心支撑钢框架结构影响系数,TU391
- 基于蠕变特性土质边坡稳定性的研究,TU43
- 基于最佳侧移刚度分布的多高层钢框架结构抗震设计方法,TU973.13
- 硬盘磁头悬架窝点与舌尖的微动接触力学分析及仿真研究,TP333.35
- 强震作用下大跨度拱形立体桁架结构动力强度破坏研究,TU399
- 楼板开洞的多高层抗震弹塑性能研究,TU973.31
- 振动摆动辗压成形体积效应机理研究及计算机模拟,TG376
- 框架结构静力与动力弹塑性抗震性能分析,TU352.11
- 塑性铰长度对平面框架滞回耗能计算影响分析,TU313
- 隔爆电气设备外壳的设计与有限元分析,TD684
中图分类: > 工业技术 > 金属学与金属工艺 > 金属压力加工 > 一般性问题 > 压力加工工艺
© 2012 www.xueweilunwen.com
|