人人文库网 > 毕业设计 > 毕业论文 > 微细电火花加工数值模拟【25500字】【优秀机械毕业设计论文】

数值模拟与影响铝电解槽磁热耦合问题.doc

微细电火花加工数值模拟【25500字】【优秀机械毕业设计论文】

资源目录

微细电火花加工数值模拟【】【优秀机械毕业设计论文】.rar

翻译

毕业论文翻译封面.doc---(点击预览)

数值模拟与影响铝电解槽磁热耦合问题.doc---(点击预览)

sdarticle (4).pdf---(点击预览)

压缩包内文档预览：

编号：606115 类型：共享资源大小：5.62MB 格式：RAR 上传时间：2016-02-28 上传人：木*** IP属地：江苏

50
积分

版权申诉 word格式文档无特别注明外均可编辑修改；预览文档经过压缩，下载后原文更清晰！ 立即下载

关键词：: 微细电火花加工数值模拟摹拟优秀优良机械毕业设计论文

资源描述：

文档包括:
说明书一份,51页,25500字左右.
任务书一份.
开题报告一份.
翻译一份.

摘要:
微细电火花加工是一种微细特种加工工艺方法，在微机械以及微细结构加工领域具有举足轻重的地位。由于RC微能电源简单可靠，在微细电火花加工中得到了大量的应用。为了深入对微细电火花加工机理的认识，进一步提高微细电火花加工能力和效率，本文进行了相关的数值模拟和试验研究，在理论研究上和工程应用上取得了一些结论。
本文主要研究内容如下：
1．针对RC微能电源，建立了极间电场物理模型，采用数值模拟方法，探讨了电极形状、加工间隙等对电场分布的影响，掌握了相关的电场分布规律。
2．根据传热学理论，建立了微细电火花加工工件表面的热物理模型，采用数值模拟方法，考虑热物性参数、对流换热、相变潜热等随温度变化的影响因素，掌握了工件表面温度场分布，得出了微细电火花加工极间温度场的分布规律。
3．在上述研究基础上，采用微细电解方法在线制作柱状钨电极，合理选择工艺参数，进行了微细电火花加工试验，并分析了电参数对加工尺度和效率的影响，获得了三组应用于不同加工场合下优化的电参数。

关键词：微细电火花加工，加工能力，效率，数值模拟，电场，温度场，电参数

指导老师签名：

ABSTRACT:
Micro EDM(electric discharge machining) is a kind of special machining technique for micro machining .It has the pivotal status in the domain of micro machine as well as the micro structure processing . To raise the processing scale and the efficiency, this article has carried on the related numerical simulation and experimental study, and has made some progress in the fundamental research as well as the engineering application.
Main content of this article is as follows:
1. Using the numerical simulation method, electric field distribution influenced by the electrode geometry and the machining polarly has been discussed. Also, the related electric field distributed rule of the micro EDM has been got.
2. According to heat transfer theory, the thermal physical model of work piece’s surface in micro EDM has been established, using the numerical simulation method, considering some factors influenced along with the changing of temperature such as the thermal natural parameter, the heat convection, the latent heat and so on, the temperature field distribution of the work piece’s surface has been grasped, and distributed rule of the interelectrode temperature field in micro EDM has been obtained.
3. Basing the research foundation above, using the micro electric chemical machining to manufacture columnar tungsten electrode on-site, choosing the reasonable technological parameter, micro EDM experiments have been carried on. Also, the influence to the working accuracy and the efficiency caused by electrical parameter has been analyzed. What’s more, three groups optimized electrical parameter have been obtained to apply in the different manufacturing places.
Key words: Micro EDM, Processing capacity, Efficiency, Numerical simulation, Electric field, Temperature field, Electrical parameter

第一章绪论 1
1.1 引言 1
1.2 数值模拟在电火花蚀除机理研究中的现状 2
1.3 微细电火花加工的现状和发展 5
1.3.1 微细电火花微能电源的介绍 5
1.3.2 微细电火花加工的研究现状和发展 6
1.4 课题背景意义 9
1.5 本文主要内容及章节安排 10
第二章微细电火花加工数值模拟 12
2.1 基于微能电源电火花加工极间电场分析 12
2.1.1 概述及物理模型 12
2.1.2 电火花电场数学模型 13
2.1.3 电火花电场的求解 14
2.1.4 ANSYS电场的仿真 16
2.1.5 影响电场分布因素的探究 20
2.2 放电状态下的极间的物理过程 22
2.2.1 温度场模型的建立 22
2.2.2 温度场模型的求解 24
2.2.3 温度场的仿真 27
2.3 本章小结 32
第三章微细电火花加工试验研究 33
3.1 试验设备—微细电火花机床介绍 33
3.2 试验电极的制作 34
3.3 常用材料性质 35
3.4 极性及电参数对微细电火花加工的影响 36
3.4.1极性对微细电火花加工的影响 36
3.4.2 电容对微细电火花加工的影响 39
3.4.3 空载电压对微细电火花加工的影响 40
3.4.4 脉宽对微细电火花加工的影响 42
3.4.5 脉间对微细电火花加工的影响 45
3.4.6 试验结果及讨论 47
3.5 本章小结 50
第四章总结和展望 46
4.1 全文总结 46
4.2 研究展望 47
参考文献 48
致谢 51

毕业设计（论文）任务书

I、毕业设计(论文)题目：

微细电火花加工数值模拟

II、毕业设计(论文)使用的原始资料(数据)及设计技术要求：

本课题旨在对微细电火花加工进行数值模拟。电火花加工是利用工具电极与工件间
火花放电对工件进行逐层蚀除来实现零件的加工。由于电火花加工受多种因素综合影响、
制约着整个放电过程的进行。利用有限元分析软件，可对加工过程进行初步模拟分析。
研究主要要求如下：
（1）加工材料为不锈钢，厚度为200μm，加工电极为直径50μm的钨电极。
（2）加工间隙为5μm，加工孔径为150μm，加工速度为30μm/min。
（3）加工介质为火花油，其介电常数为2.05。
（4）进行电场分析，加工电压为0-10V。
（5）进行加工温度场分析，设置放电脉宽为50μm、100μm。

III、毕业设计(论文)工作内容及完成时间：
1. 查阅相关资料，外文资料翻译（6000字符以上），撰写开题报告
第1周—第2周
2．建立电场分析的数学和物理模型第3周—第5周
②电极不同形态对电场间电势及电场分布的影响；
③极性对电场分布的影响。
2011.03.21～ 2周
3．分析加工电场分布，电极状态对电场影响第6周—第9周
4．建立温度场分析的数学和物理模型第10周—第11周
5．分析脉宽、加工时间等对温度场的影响第12周—第13周
6．撰写毕业论文第14周—第16周
7. 毕业设计审查、毕业答辩第17周

Ⅳ 、主要参考资料：
[1] 赵万生，王振龙，郭东明.国外特种加工技术的最新进展[J]，电加工，1999，（5）：9-12
[2] 余承业.特种加工新技术[M]，北京，国防工业出版社，1995：28
[3] 基于RC微能电源的微细电火花加工数值模拟和试验研究.南京航空航天大学硕士学位论文，2008
[4] Shankar P. Analysis of spark profiles during EDM process[J], Machining Science and Technology,
1997,(1):195~217
[5] 孔祥谦.有限元法在传热学中的应用[M]，北京，科学出版社，1998：103

航空与机械工程系机械设计制造及其自动化专业类 0781052 班
学生（签名）：

填写日期： 2011 年 01 月 03 日

指导教师：

助理指导教师(并指出所负责的部分)：

机械设计制造及其自动化系（室）主任：贺红林

附注:任务书应该附在已完成的毕业设计说明书首页。

微细电火花加工数值模拟

内容简介：: 铝电解槽热磁耦合问题数值模拟 Y. , M. J. 桑联邦理工学院，分析与科学计算研究所，瑞士洛桑 1015 第 8 站 2006 年 12 月 27 日初稿； 2008 年 2 月 4 日修订； 2008 年 2 月 8 日被收录； 2008 年 2 月 29 日可网上搜索摘要本文运用一系列偏微分方程对铝电解槽的热磁耦合行为进行了数值模拟。热模型被认为是一个由焦耳效应引起的非线性对流扩散热方程组成的两阶段史蒂芬问题。该磁流体动力领域的主导是纳维进化组合（切尔诺夫）用于获取电解槽仿真壁架的温度和凝固层剖面的稳态解。利用有限元方法的数值算法来获取流体速度，电势，磁感应和温度。同时也利用了迭代算法和三维数值模拟结果。 2008年爱思唯尔公司保留所有权利。关键词：铝电解；切尔诺夫组合；热方程；磁流体力学；壁架；凝固 1. 绪论本文研究了由电解槽热磁耦合作用模型引起的相位变化问题。在一个利用霍尔属部分是由三氧化二铝电解融化在熔融冰晶石材质的槽中制造而成的 1。该电解槽中产生了多种现象，图 1为一个横截面示意图。电解槽中稳定的电流通过铝液在阳极和阴极棒之间产生。送到槽中的电流产生重要的磁场，该磁场连同电解槽中流通的电流共同产生一个维持这两种导电液体耦合运动的拉氏力量作用域。电解槽中会产生磁流体动力学相互作用。另一方面，由槽体中的电阻率引起焦耳效应，热源也随之产生。冰冻层冰冻层电解液阳极块铝液阴极层图 1. 铝电解槽横截面所谓的壁架在固态槽壁层上建立。这些壁架能够避免电解槽侧壁腐蚀性点解，并降低电解槽热损耗（见 2第 23页）。此外，它的轮廓严重影响磁流体动力稳定性，引起铝液和槽体接触面振荡，降低电流效率。因此最佳的层剖面是电解槽侧壁设计的目标之一。冶炼槽内的热凝固问题已经被几个专家解决 3据我们所知，在热磁耦合上，该问题一直没有被重视。本文的目的是解决类似的耦合问题。我们期待，关于该问题的详细资料记入在萨法论文集 6。数学上，该问题解决了偏微分方程组、麦克斯韦方程组和纳维中，偏微分方程组包含由焦耳效应引起的热方程，麦克斯韦方程组以导电率作为温度的函数。铝液和和槽体之间的接触面是未知的。壁架被认为是电绝缘体，热模型是静止的两阶段史蒂芬问题。本文大纲如下：第 2节介绍物理模型，第 3节给出算法，第 4节得出数值计算结果。 2. 模型为了介绍该模型，我们首先描述一些几何和物理量。概括描述几何图形定义如图 1所示。下面介绍物理符号： 21 ：流体和固态层， 21 ：电极，：表示电解槽的域另外，我们定义如下接触面： 21 ：铝液和槽体之间的自由接触面，未知，， i 1,2， 21 ：电极的外边界。我们必须处理的未知物理场列举如下：流体动力场： u ： i 1,2，（固态层中 u 为 0）， p ：压力。电磁场： b ：磁感应场， e ：电场， j ：电流密度。热场： H ：总热能，：温度。材料属性定义如下：：密度， b与：槽体内、外导电率，：流体粘度， 0：空隙导磁率， k ：热导率，热容，：潜热。物理假设该模型需要以下基本假设： 1. 各流体不相融，不可压缩，并且遵守牛顿定律。 2 在每个域内，i = 1,2，各流体遵守静态纳维 3 电磁场满足静态麦克斯韦方程组，此外，欧姆定律应该在整个电解槽内有效。 4 槽外的电流密度已知（即阴极棒内的电流）。 5 导电率是液体和电极部分的温度的函数。 6 粘度，密度和比热容 7流体1和固体2的体积为给定值（质量守恒）。 8电解槽中的流产生的焦耳效应提供唯一的热源。 9忽略化学反应的影响 7，马朗戈尼效应 8,9，表面张力以及气流的存在。流体动力学问题在这一部分，我们考虑温度场和电磁场 j ，并且磁感应场 b 为已知。我们选择的用 y,x ,yh:y ,的参数化形式表示铝液和槽体间的未知接触面，其中 D 通常是一个与铝阴极界面的参数化相对应的矩形区域。考虑到 h ，我们用 2 ,i,和的相互关系。根据假定 7得出如下关系： 1D yh 其中1 单位法线指向为 1。我们定义水动力场的标准方程组如下： 1 i v, ，（ 1） 1 （ 2），（ 3）其中 3,2,1,21 ,T 。这里 (.,.)通常是方程（ 1） -（ 3）对应于第 1和第 2条假设。我们通过引进包含流体的域于任意场 w ， w 的跳跃，即 al u m in iu mb a t h 域中，u 和 p 具体为： u （ 4） 0 ，（ 5） 0D2 。（ 6） 2,1 的流体部分只是所有域固前的一个子域。为了解决一个固定域们使用包含惩罚函数的“虚拟域”方法。之后会定义液体和固体中的速度和压力。我们将术语斯托克斯方程，度函数的固相组分。函数 K 由科泽尼定律给定： 32 21，其中， P 为平均孔径， C 为通过实验确定的常量（见 10）。修改方程（ 1）为 1 i v, （ 7）液相状态下 K 为 0，上述公式简化为一般的纳维相对于其他状态，糊状区内 K 可能很大，并且上述方程模拟了达西定律：。当 1们得到并且固态区内 u 为 0。最终得到流体动力学问题已知 j ， b 和，求得 u ， p 和 h ，并且满足以下条件 1 i v, ，（ 8） 1 （ 9），（ 10） u （ 11） 0 ，（ 12） 0D2 ，（ 13） 1D yh （ 14）电磁问题我们假定速度场 u 及温度场已知。根据法拉第定律，我们令 e 为 0，电场由e 给定，其中表示在中计算得到的电势场。我们仍然用 u 表示速度的连续延拓（在中为 0），同时考虑到安培定律：令及欧姆定律：在中，因此，我们给出电守恒定律： bu 。我们用n表示运算 ,n ，其中 n 是的外部单位法线。关于电势，我们介绍以下边界条件： n ， 2jn ， 1 其中0为电流 j 的一个函数，磁感应强度萨伐尔关系求得： 30，其中0 电磁问题 u 和 h 已知，通过以下公式求得， b 和 j bu ，（ 15） n ，（ 16） 2jn ，（ 17） 1，（ 18），（ 19） 30。（ 20）温度问题我们假定流体动力场 u 和电磁场 j 已知。我们正在寻求的稳定解将在这里作为一个随时间变化的热方程被求得。因此，在本小节我们介绍了进化的热模型。在对流扩散问题中，凝固前（固液分界线）的位置事先未知，因此需要作为解法的一部分来确定。该问题被广泛地称为“ 史蒂芬问题”，并且是高度非线性的。为了克服有关史蒂芬接触面条件的非线性困难，我们定义一个焓函数，它表示每单位体积的物质的总热量。焓可以用温度，潜热和固相分数： 1 。（ 21）由于焓 H 为单调函数，我们引进函数，由以下关系定义： H 。（ 22）函数 H 是由在 H, 列表中的值经过数学处理（插值法）计算得出的，上述列表值与由方程 21 得出的反比关系式 1H 相对应。通过这种关系，我们可以将该问题表示成有关温度和焓的史蒂芬问题，形式如下： ,d ，（ 23） H ，（ 24）该问题是一个非线性对流扩散系统。表达式 ,u 用表示 u 的数积，为仅由焦耳效应提供的热源。表达方式如下： 2 。（ 25）考虑到分配性，温度及可以用标准数值技术来解决相变问题。温度遵守罗宾边界条件：，（ 26）其中n是在指向的单位外法线的方向导数，是由空间和温度决定热传递系数，a是外界温度。热传递是由于对流和辐射。而辐射是由如下表达式间接考虑： C/321 其中1c，2 HH 0 的初始条件是假定的。对于一个给定的表示积分时间的标量值 T ，表达式如下： T ,0a n d,0 。温度问题 u ， h 和 j 已知，通过以下公式求得和 H ： in,d ，（ 27），（ 28），（ 29） 0fo r, （ 30）完整问题我们刚刚描述了流体力学，电磁和热力学的问题。在每个具体域中，我们都假定其他域已知。我们要解决的问题就是找到同时满足上述三个问题的速度 u ，压力 p ，电势，焓 H 和温度；函数 H ，， ss 均给定，并且已知常量，0，，PC，a，1V，和。 3. 数学方法上述的数学问题的数值解是基于一个迭代过程的，在迭代过程中，我们交替进行未知的三种域的计算：流体域电磁域热域在本节中我们提出了含用有限元法进行空间离散的全局“伪进化”算法被用来解决三域耦合问题。流体动力场计算流体动力学问题通过迭代求解。在每一个求解步骤中，我们首先解决接触面无正常受力平衡条件的固定形状问题，然后利用非平衡法向力更新接触面位置。解决替代应用问题的两个步骤如下：第一步：通过给定的几何结构 h ，并且考虑到接触条件，解决流体动力学问题如下： , ， 0, ， t 为 h 的切向量，该问题之后会通过弱公式化简单表达。第二步：更新接触面位置，以便验证接触面所受法向力的平衡，通过以下表达式选择 : ： s , , 这里我们用示通过以下条件求得的常量： D 。利用迭代法计算 1 1 1 1m 和 1上述函数的函数值通过之前的迭代 m 求得。设定 ,和 ,，并且通过下列表达式定义求解步骤 i v, 1111 ，（ 31），（ 32） mm ，（ 33） 0,1 ， t 为 h 的切向量，（ 34） ,11 ，（ 35） f C 11 ，（ 36） 0 （ 37）注意到这种算法的停止条件是基于 mm 1 的规范估计 1H ，它必须小于公差。电磁场计算磁感应强度 b 直接取决于电流 j ，间接取决于位势场，对于一个已知速度场 1求解 m 的步骤中，我们用迭代法确定。迭代法中，利用公式（ 15）、边界条件（ 16） -（ 18）和值计算 1m ，利用公式（ 19）相继计算 1随后，我们利用毕奥萨伐尔定律求得作为 1函数的1值。求解步骤 1 ，（ 38） o n 01，（ 39） 201 ，（ 40） 11 0 m，（ 41），（ 42） 3101 3。（ 43）停止条件是基于 1的规范估计 2L ，它必须小于公差。热场计算如前所述，我们用伪进化类型作为收敛于温度问题（ 27） -（ 30）稳定解的数学均值。在运用（ 27） -（ 30）时，利用半隐式法离散化得到： 1 ，（ 44）其中， m 和在时 H ，和的值；为离散化的时间间隔。为了关闭系统（ 44），我们利用切诺夫方法，即： 11 ，（ 45）其中是正松弛参数。通过在（ 44）中替换（ 45），我们得到一个计算时间1： 11 ，（ 46）只要满足下列条件，该方法就是稳定的 11： Rs 在上面的式子中， 1m 不是在时间 1糊状区温度的良值，其良值可从 1求得。为了避免可能出现的误解，用表示前者。在时间离散化形式中，我们假定已知在时间区间内的另外用以下方法计算求解步骤 m ， 1 1m ： d i v 11 ，（ 47） 1a1 ，（ 48） 11 ，（ 49）。（ 50）伽辽金构想利用伽辽金构想对三套方程组方法适用于在一个四面体网格上进行一阶分段多项式有限元法。图示用于计算的四面体网格利用经典的稳定有限元法1P（见 12）对纳维题值求解，用含一阶分段多项式的标准有限元方法对位势问题值求解。该小节中，我们把重点放在与热问题相应的有限元方法上。考虑到局部沛克莱数的大小（在本例中大约为 1000），我们使用定法（流线迎风彼得罗夫 13。定义有限元空间 ;0 ，其中表示网格重叠。与方程集（ 47） -（ 50）对应的有限元由下列条件给定：找到 3, H 满足 ,dd,dx ，（ 51） 11 ，（ 52） r ，（ 53）其中的插值，公式如下： 0 00 ，（ 54）图 2. 域四面体网格 1 11K，（ 55）其中 K 的大小，变量我们分别用 mm 和表示时间间隔 m 内的导热系数和传热系数。注。很明显，焓 H 在先验未知的凝固前有一个突增，但我们可以用连续函数子空间内的1 。因此，该近似值表示了在一个精确焓值突增的小区间内的焓值差。我们注意到，即使建立了很好的离散问题，逼近收敛于 H 的焓近似值仅在 4. 数值计算结果我们使用广义最小余数法解决流体动力学问题温度问题矩阵系统。另一方面，由于与电场问题关的该矩阵系统是对称且正定的，我们利用代数多重网格法利用计算机（奔腾（ R） 4，行磁 10小时内可获得全局算法的收敛性。电势计算相关结果见图内温度分布。壁架模型见图我们可清晰观察到固化前一个小区间内由定法得出的相关数值扩散的影响。值得指出的是，这张图片与代表速度场的图片（图 6）具有一致性，特别是壁架大的位置与各场中数值小的位置也是一致的。很容易观察到在液体分数小于 1的部分域内，速度场的数值计算结果符合达西定律。这表明，在流体动力学模型中，利用含有速度场惩罚值的“辅助域”方法是有效的。用多个元素对流体层沿深度进行离散化，以便更准确地表示流体动力场。另外如前所述，界面铝保证了垂直受力平衡（见第步）。液体深度误差的主要误差值是 6。在铝们取得最高的对流效果，最佳壁架厚度见图保证总散热量与内部产热量的一致性是很重要的。为了保证结果与槽稳态条件相一致，关键是保证实现上述条件。总散热量等于各部分 .,1, （其中 N 表示槽中电流穿过的各部分数量）产生的焦耳热之和：总散热量。总散热量相当于槽边界处的对流耗散： a 。图 3. 槽内电势结果图 4. 槽内温度结果焦耳总散热量的数值误差为该误差值相当于槽切片内使用有限元分析代码进行电热计算的误差，见 14。图 5. 壁架内液体分数图 6. 固化前速度场（上图）和固化后速度场（下图），平均值 =s 图 7. 壁架模型和电热计算的液体分数使用另一种方法求得壁架模型。该方法依赖于无速度场的电热计算。利用金属槽体内的人工导热系数进行对流效应模拟。对比图得对称壁架如图而更易观察固体壁架结构内速度场的影响。 5. 结论这项工作是在前人研究基础上的延伸（ 15， 16），目的在于引出判断霍尔在上面提到的参考文献中，这些标准是从磁流体动力学线性方程组定态解的频数分析得到的。它产生于 15 和 16中进行的数值研究，这些标准的稳定性在很大程度上依赖于求得的定态解的精确度。精确度不仅取决于正确的数值方法，也取决于电解槽特征表述与模型的妥善性。本文中，壁架模型和速度场的温度分度影响已考虑在内。通过观察分析上述结果，得到以下结论：流体动力场的影响是一个决定电解槽热行为的重要因素。从图图可以看到，速度场对壁架形状有很大的影响。忽略流体应力张量对壁架的侵蚀作用完成的计算。该问题应在将来解决。虽然处理多域相互作用问题有很大难度、几何条件也很复杂，切尔诺夫方法在解决热磁力耦合问题上很稳定。最后，我们注意到，多个作者进行了热膨胀影响下的铝电解槽钢壳热 19在这些研究中从未有热磁耦合的计算。本文中所示温度场计算对于之后由萨法等人（见 22）为展示热对流对槽体结构力学的影响以及速度场和钢壳机械变形的相关性而进行的弹性热计算而言是富有成效的。致谢作者衷心感谢瑞士国家科学基金会和加拿大铝业公司的支持。参考文献 1 P. J. H. A. J. 的电解霍尔 He三编。，铝，出版社， 2001。 2 K. H. 杜塞尔多夫生产霍尔 1986. 3 V. M. J. 示范热电和磁流体的铝 500 电解数学模型细胞在：第 42届会议的拿大温哥华， 2003。 4 L. 自由边界的存在性问题，在热电造型的铝溶液电解槽 , 欧洲应用数学学报 14 (2003) 201216. 5

温馨提示:
1: 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

人人文库网所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。