喷油嘴微小孔磨流流加工数值毕业论文毕业设计

喷油嘴微小孔磨流流加工数值毕业 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 毕业设计 编号 本科生毕业论文 喷油嘴微小孔磨粒流加工数值模拟 Micro-hole of Nozzle Abrasive Flow Machining Numerical Simulation 毕业设计(论文)原创性声明和使用授权说明 原创性声明 本人郑重承诺:所呈交的毕业设计,论文,～是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知～除文中特别加以标注和致谢的地方外～不包含其他人或组织已经发表或公布过的研究成果～也不包含我为获得 及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体～均已在文中作了明确的说明并表示了谢意。 作 者 签 名: 日 期: 指导教师签名: 日 期: 使用授权说明 本人完全了解 大学关于收集、保存、使用毕业设计,论文,的规定～即:按照学校要求提交毕业设计,论文,的印刷本和电子版本,学校有权保存毕业设计,论文,的印刷本和电子版～并提供目录检索与阅览服务,学校可以采用影印、缩印、数字化或其它复制手段保存论文,在不以赢利为目的前提下～学校可以公布论文的部分或全部内容。 作者签名: 日 期: 学位论文原创性声明 本人郑重声明:所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外～本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体～均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。 作者签名: 日期: 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定～同意学校保留并向国家有关部门或机构送交论文的复印件和电子版～允许论文被查阅和借阅。本人授权 大学可以将本学位论文的全部或部分内容编入有关数据库进行检索～可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 涉密论文按学校规定处理。 作者签名: 日期: 年 月 日 导师签名: 日期: 年 月 日 指导教师评阅书 指导教师评价: 一、撰写,设计,过程 1、学生在论文,设计,过程中的治学态度、工作精神 ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、学生掌握专业知识、技能的扎实程度 ? 优 ? 良 ? 中 ? 及格 ? 不及格 3、学生综合运用所学知识和专业技能分析和解决问题的能力 ? 优 ? 良 ? 中 ? 及格 ? 不及格 4、研究方法的科学性,技术线路的可行性, 设计方案 关于薪酬设计方案通用技术作品设计方案停车场设计方案多媒体教室设计方案农贸市场设计方案 的合理性 ? 优 ? 良 ? 中 ? 及格 ? 不及格 5、完成毕业论文,设计,期间的出勤情况 ? 优 ? 良 ? 中 ? 及格 ? 不及格 二、论文,设计,质量 1、论文,设计,的整体结构是否符合撰写规范, ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、是否完成指定的论文,设计,任务,包括装订及附件,, ? 优 ? 良 ? 中 ? 及格 ? 不及格 三、论文,设计,水平 1、论文,设计,的理论意义或对解决实际问题的指导意义 ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、论文的观念是否有新意,设计是否有创意, ? 优 ? 良 ? 中 ? 及格 ? 不及格 3、论文,设计说明书,所体现的整体水平 ? 优 ? 良 ? 中 ? 及格 ? 不及格 建议成绩:? 优 ? 良 ? 中 ? 及格 ? 不及格 (在所选等级前的?内画“?”) 指导教师: (签名) 单位: (盖章) 年 月 日 评阅教师评阅书 评阅教师评价: 一、论文,设计,质量 1、论文,设计,的整体结构是否符合撰写规范, ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、是否完成指定的论文,设计,任务,包括装订及附件,, ? 优 ? 良 ? 中 ? 及格 ? 不及格 二、论文,设计,水平 1、论文,设计,的理论意义或对解决实际问题的指导意义 ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、论文的观念是否有新意,设计是否有创意, ? 优 ? 良 ? 中 ? 及格 ? 不及格 3、论文,设计说明书,所体现的整体水平 ? 优 ? 良 ? 中 ? 及格 ? 不及格 建议成绩:? 优 ? 良 ? 中 ? 及格 ? 不及格 (在所选等级前的?内画“?”) 评阅教师: (签名) 单位: (盖章) 年 月 日 教研室(或答辩小组)及教学系意见 教研室(或答辩小组)评价: 一、答辩过程 1、毕业论文,设计,的基本要点和见解的叙述情况 ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、对答辩问题的反应、理解、表达情况 ? 优 ? 良 ? 中 ? 及格 ? 不及格 3、学生答辩过程中的精神状态 ? 优 ? 良 ? 中 ? 及格 ? 不及格 二、论文,设计,质量 1、论文,设计,的整体结构是否符合撰写规范, ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、是否完成指定的论文,设计,任务,包括装订及附件,, ? 优 ? 良 ? 中 ? 及格 ? 不及格 三、论文,设计,水平 1、论文,设计,的理论意义或对解决实际问题的指导意义 ? 优 ? 良 ? 中 ? 及格 ? 不及格 2、论文的观念是否有新意,设计是否有创意, ? 优 ? 良 ? 中 ? 及格 ? 不及格 3、论文,设计说明书,所体现的整体水平 ? 优 ? 良 ? 中 ? 及格 ? 不及格 评定成绩:? 优 ? 良 ? 中 ? 及格 ? 不及格 (在所选等级前的?内画“?”) 教研室主任(或答辩小组组长): (签名) 年 月 日 教学系意见: 系主任: (签名) 年 月 日 摘 要 喷嘴嘴作为电喷发动机的关键部件，它的工作好坏将严重影响发动机的性能。随着人们对于产品小型化和加工过程微型化的追求，磨料流加工技术随之应运而生。本文介绍磨料流加工的基本工作原理及其特点，利用GAMBIT进行模型创建及网格划分工作，利用FLUENT软件对喷油嘴零件进行数值模拟，通过数值 模拟获得磨粒流加工喷油嘴的动态压强、静态压强、速度等仿真曲线，并对其进行了分析与总结。通过对零件进行数值模拟技术，我们可以检验设计中的失误，提高设计的效率和减少设计成本，实现零件的最优化设计。 关键词:磨料流工 FLUENT GAMBIT 数值模拟 I Abstract Nozzle is the key component to EFI engine ,its quality of work will seriously affect the engine performance. With the product miniaturization and the pursuit of miniaturization process, it come into being along Abrasive Flow Machining technology. This article describes the AFM basic working principle and features, we use GAMBIT mesh to create and work, simulate the nozzle part in fluent, by numerical simulation we can obtain the simulation curve which is about the dynamic pressure. static pressure. velocitye and so on ,and get its analysis and summary. Through the numerical simulation of Parts, we can test the design errors, improve design efficiency and reduce design costs, and gets the optimal design of the parts at last. ;;;Abrasive Flow Machining FLUENT GAMBIT SimulationKeywords: I 目录 摘 要 ...................................................... 1 ABSTRACT ................................................... III 目录 ........................................................ IV 第1章 绪论 .................................................. 1 1.1引言 .......................................................... 1 1.2 研究磨粒流的目的和意义 ........................................ 4 1.3 磨粒流技术国内外的发展现状 ................................... 6 1.4 小结 .......................................................... 9 第2章 FLUENT软件............................................ 10 2.1 GAMBIT软件介绍 ............................................. 10 2.2 FLUENT软件概述 ............................................. 13 2.3 小结 ........................................................ 15 第三章 数学模型和计算方法 ..................................... 16 3.1数学模型 ..................................................... 16 3.1.1 控制方程数学模型 ........................................ 16 3.1.2 流体相数学模型 .......................................... 17 3.1.3 颗粒相碳化硅数学模型 .................................... 17 3.1.4 湍流数值模拟数学模型 .................................... 18 3.2 计算方法 .................................................... 18 3.2.1 有限差分法 .............................................. 19 3.2.2 有限元素法 .............................................. 19 3.2.3 有限体积法 .............................................. 19 3.3 FLUENT求解方法的选择 ....................................... 19 3.4 小结 ........................................................ 20 第四章 湍流模型与边界条件 ..................................... 21 4.1 湍流模型 .................................................... 21 4.2 边界条件 .................................................... 23 4.3 小结 ........................................................ 25 第五章 数值模拟算法与多相流的选择 .............................. 26 5.1 求解器的选择 ................................................ 26 5.2 离散格式的选择 .............................................. 26 5.3 压力速度耦合方法的选择 ...................................... 26 5.4 多相流仿真的选择 ............................................ 26 5.5 小结 ........................................................ 27 第六章 喷油嘴数值模拟 ........................................ 28 I 6.1 喷油嘴二维数值模拟 ........................................... 28 6.2 小结 ......................................................... 33 结 论 ..................................................... 34 参考文献 .................................................... 35 致谢 ........................................................ 36 I 第1章 绪论 1.1引言 尽管磨料流加工技术已有50余年的研究和应用历史，而且国内外的一些研究工作者从20世纪70年代以来一直在不懈地从事着该技术的研究工作，但该技术的研究和应用基本上还停留于实验研究方面，或者说仍处于研究和应用的初级阶段，由于其加工过程只能由工程技术人员主动地进行控制，所以目前主要还是适于单件加工生产，很难推广应用于大规模生产。这主要是因为目前还缺乏对磨料流加工过程的综合复杂特性的支持。而磨料流叫那个所用夹具的设计是决定磨料流加工效果甚至是决定磨料流加工成败的一个关键因素。如夹具上磨料流的进出口位置、大小、方向、数量等的设计确定，对于磨料流在工件行腔或者夹具型腔的压力与流速分布具有重要影响，从而影响其加工的质量和效率。这样的一些设计问题目前基本完全依赖于设计人员个人的实践经验，其结构设计是否合理甚至是否达到最优设计，尚不能在设计阶段作出可靠的分心评价，而只能通过实际 【】1加工，根据加工的结果对其结构设计进行可行性评价。实际中大多数的设计尤其是复杂的设计很难保证一次设计的成功率，往往要经过多次的实验和反复修改，这不仅严重影响了加工生产周期的缩短，而且在无法进行结构设计修改的情况下容易造成夹具或模具的报废，从而造成很大的浪费，增加生产成本，延误市场竞争的战机。 从现代市场竞争的观点来看，也就是制约着企业对市场需求反应的灵敏性，即从硬件技术方面制约先进制造技术的理念在一个企业的实施。实际中需要用磨料流进行表面抛光加工的工件的表面形状和尺寸各异，且差异很大，如何根据其材料、形状和尺寸大小以及表面质量初始状态等来合理地选用磨料流中磨料类型、粒度及其粘结剂、润滑剂等的合理配比以及磨料流加工的合理压力、流速、往复行程次数等目前尚不能从理论上给出一个合理的答案，主要还是根据经验和实验来决定。这些问题正是磨料流加工技术目前难以应用于大批量生产以及推广应用的主要困难和技术障碍，可以相信这些问题如能较好给予研究解决，一定能够大大推动磨料流加工技术的实际使用，从而产生巨大的经济和社会效益。 要深入理解和控制磨料流加工过程，就必须更多地了解金属去除与表面生成的机理，过程模型、加工能力和控制方法，而要做到这些并解决上述存在的主要问题，就必须一方面借助于磨料流加工的传统实验研究方法、大量从事磨料流加工机理和工艺规律的实验研究和数据积累，同时在此基础上充分利用现代数值计算和计算机虚拟仿真技术，开着磨料流加工虚拟技术的研究，从理论上对磨料流加工的进行虚拟仿真，实现虚拟制造。利用磨料流加工虚拟技术可在设计阶段实 1 现对磨料流加工夹具设计 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 的可行性评价以至实现其结构性和磨料流加工过 【】2程的优化和监控。 随着产品的不断精密化，人们对产品表面质量及功能性边缘的要求不断提高，传统的加工工艺已无法达到其精度要求。在加工过程中，工件的加工稳定性如同公差一样要尽量完美，并且有着同样程度的要求。精密零件制造中的最终精加工是一种劳动强度大而不易控制的过程，它在全部制造成本中所占的比重非常高。磨粒流加工技术是一种能够保证精度、效率、经济的自动化光整加工方法，是解决精密零件最终精加工的一种有效手段。它是以一定的压力强迫含磨料的粘弹性物质性体，称其为柔性磨料或粘弹性磨料，通过被加工表面，利用其中磨粒的刮削作用去除工件表面微观不平材料而达到对工件表面光整加工的目的。磨粒流加工技术是国外20世纪70年代以来开始推广应用的一种先进光整加工技术，宇航用的液压阀体孔道多，孔道相互交叉，台阶孔，交叉孔内的毛刺很难去除，柔性磨体加工技术就是为解决宇航用液压阀体孔道内的毛刺去除问题从美国发 【】3展起来的。这项技术出现后很快在液压、模具、航空、纺织机械、汽轮机、齿轮等机械行业中获得应用。后来，日本、西欧和前苏联都引进和采用了磨粒流加工技术。 我们国家曾将磨料流加工技术列为七五攻关技术项目之一进行了研究，取得了一定的实验研究成果，在20世纪80年代初成功的研究开发了磨料流加工机床、加工所用的柔性磨料及加工技术，成果推广应用于多个企业。 磨粒流工艺的发明，可以说创立了内外表面边角处理的新概念，使得金属表面处理更具特色。其更大优势在于，随着现代制造技术微型化，密集化，越来越多的零件结构呈现出封闭的无规则的同时又要求高质量表面的各种型腔、通道等，磨粒流工艺在这些部位的微量加工和抛光方面可谓独树一帜。不但如此，磨粒流工艺还可以成功地对各种材料进行微量研磨加工，包括较软的有色金属直至坚韧的镍合金、陶瓷和硬质合金等制作的零件，非常适应现代制造材料和制造技 【】4术的发展。 磨粒流加工是利用磨粒流中的磨砂充作无数的切削刀具(自锐性)，以其坚硬的锋利的棱角对工件表面进行反复切削，从而达到一定的加工目的。在工艺实施中，通常采用两个相对的磨粒缸使磨粒在零件和夹具所形成的通道中来回挤动。磨削作用就产生在流体受到限制的部位，即挤压部位。当磨粒均匀而渐进地对通道表面或边角进行工作时，产生去毛刺、抛光及倒角的作用。 磨粒流加工工艺中有3个重要环节。 (1)挤压研磨机床。固定工件和夹具，在一定的压力作用下，使磨粒通过被 加工表面，达到研磨，去毛刺，倒角及抛光的目的(图1-1)。由机床控制挤出压 2 力，压力范围从7一 224 kg / cm。 (2)磨粒。由柔性的半固态载体和一定量磨砂拌制而成，有不同粘度、流变、磨砂粒度和密度。最常用的磨砂是碳化硅。根据被加工材料，还可选择立方氮化硼、氧化铝和金钢砂。砂粒尺寸为0.005^' 1.5 mma高粘度磨粒可用来对零件的壁面和大通道进行均匀研磨。低粘度磨粒可用来对零件边角倒圆和小通道的研磨【】5。 磨粒的粘度，挤压压力和通道的大小决定了磨粒的流速，影响到研磨量，磨削均匀性和边角倒圆大小。 (3)夹具。使零件定位。并引导磨粒通达被加工部位。堵住不需要加工的部位。有些零件的磨粒流加工不需要夹具辅助，如模具等。有些加工仅需简单夹具(图1-2)。大批量零件生产所用的夹具，要设计得易于装、卸、清洗，通常须安装在 【】6分度台上。这样的夹具，一次可加工许多零件。 图1-1 磨粒流的加工原理图和去毛刺 图1-2 不同粘度的磨粒流 3 1.2 研究磨粒流的目的和意义 当今材料、计算机、精密零件的不断发展，将带动制造技术整体推进。特别以航天、航空工业为标志的高精技术密集行业，出现了许多应用新材料，新技术设计的各种零件精度更高，重量更轻，结构更复杂。而制造这些零件的高端工艺也给我们带来了挑战，使我们的光整工艺不断完善。尤其在航天领域，几乎所有大型飞机制造公司的发动机厂和维修厂都使用磨粒流工艺加工各种各样的零件。除了作为一种抛光手段，磨粒流工艺还可以对那些形状公差、质量要求极其严格的零件进行微量磨削加工。主要应用场合分布在以下几个领域:调整叶片、叶轮、油嘴、化油器和定子的气流阻力，增大涡流数值，去除激光和电火花加工所产生的重铸层，改善压缩机和涡轮零件翼形表面状态，去除毛刺，对零件的内外边角进行倒圃，齿轮去毛刺和超精抛光，对旧零件讲行修招.去除积碳.以及改善表面的招休完雄性。具体应用介绍如下。 (1)调整喷油嘴小孔流量。以下是喷油嘴图片。图7为电火花加工后小孔状况; [7]图8为同一小孔经微孔磨粒流加工后的状况。可以看出光滑的孔边缘和均匀的表面(图1-3)。 图1-3 粗加工和磨粒流加工的表面对比 喷油嘴经微孔磨粒流加工，减小了流量散差，使流量更为精确。孔的出口保持锐利，而在内径处产生均匀的圆角，并使表面光滑，轮廓完整。进口处的光洁和圆角，使小孔流量增大，寿命提高，这些经微孔磨粒流处理过的油嘴，其流量散差可控制在1%以内。 (2)改善疲劳强度。微孔磨粒流均匀倒圆了油道与油腔的交叉处，改善了疲劳强度，如图9所示。喷油器体、高压泵壳体、增压器、阀和针阀体等柴油机零件经过磨粒流加工后，在高压孔交叉处产生大的圆角((0.1mm到0.5mm以上)，改善了高脉冲产生的疲劳强度。柴油机中的高压油管，要不断经受很高压力的冲击，达1000Ba以上。因此在一些应力集中的地方会产生疲劳损坏。磨粒流方法 [8]光整表面，去除疵点，倒圆锐角，大大增加了零件的可靠性，延长了使用寿命。 (3)提高性能。汽车进气管，手工抛光其内表面时，只能先切割开，抛好以 4 后再焊接起来。而用磨粒流加工方法，不需要切割打开，就可以使磨粒挤压通过所有的管道(图1-4)。磨粒流除了使进气管道的壁面得到抛光外，还使管道内部空间增大，气流量增高。经空气推力测试，抛光后推力增加17%到23%。研磨量的多少，很大程度上取决于零件浇铸表面的粗糙度。 图1-4 磨粒流加工处理管道 磨粒流工艺己广泛用于汽车零件的精加工:进排气管、进气门、增压腔、喷油器、喷油嘴、泵、二冲程、四冲程汽缸头、涡轮壳体、涡轮叶片、花键、油头和齿轮等,下图(图1-5)就是一个涡轮发动机中涡轮抛光的例子。 各种机械结构中的传动齿轮，包括涡轮、涡杆和伞齿轮都可以用磨粒流加工。可以提高齿轮使用寿命，减少传动嗓音。 图1-5 磨粒流涡轮抛光前后对比 铸件可直接在专门的生产性磨粒流系统上抛光。这个磨粒流系统每小时抛光3U个汽缸头。计算机控制所有加工参数。旋转台面上有两个工位。当一个部件在加工时，另一个工位做装卸、准备工作。， (4)用于模具抛光。模具需要光洁的表面，以便于成型、脱模，确保产品的表面质量。这种要求很高的工序，传统上是由熟练的技术工人进行手工研磨，其不稳定性是不可避免的。磨粒流的出现使通道、型腔的高精度成为现实。磨粒流使 [9]模具表面光洁度提高而更为可靠并降低劳动成本。 由电火花加工而成的三通道铝挤出模，在磨粒流机床上经过约六分钟的加工，其表面粗糙度从Re 0.4提高到Ra 0.1 。这类模具使用细磨粒磨削，不会改 5 变工作带形状，每个齿廓都经抛光，能保持均衡的挤出速度。 电火花成形、线切割、铣削以及磨削加工的零件表面，经磨粒流加工，其表面粗糙度可提高10个数量级，该图是模具表面经电火花加工后的金相图，每一张图的状态经过2分钟抛光，其粗糙度等级可提高10个数量级。 磨粒流加工可用于抛光各种模具:挤出模、成型模、拉丝模、锻压模和冷锻模具等。 (5)磨粒流加工技巧。磨粒流加工是一种特殊的表面处理工艺，自然有其与众不同的特点。如果充分注意这些特点并加以利用，将获得最佳效果。比如磨粒流无论粘度如何，在夏天都会变稀，要达到相同的目的必须适当增加挤动次数，当然如果加工参数选择过高，将导致磨粒发热，使其切削性下降。从小孔径流量理论来说，狭窄空间流量快，故在工艺设计时要注意这一因素，例如对大径或盲孔的磨粒流加工最好设计专用喷嘴，如图1-6所示，以此来形成狭窄空间。这里需要强调的是去毛刺，抛光，倒圆角用磨粒流不一样，建议不要使用同一种磨粒流完成不同工艺过程。此外，还要倡导绿色制造，加强磨粒收集，进行环保处理。对于完成磨粒流加工的零件内残存的磨粒，只需经过液体浸泡，震动即可清除， 【】10然后再清洗，吹干。 图1-6 磨粒流加工技巧应用 1.3 磨粒流技术国内外的发展现状 1. 国外的发展现状 由德国制造的 Perfect Finish GmbH 磨粒流流体动力研磨系统，主要运用在 【】12航天及汽车工业,有着复杂几何图形合金含量较高的部件抛光及去毛剌。这种全自动研磨技术可以替代以前各种费时的人工去毛剌和抛光工序，人工操作通常无法达到长久一致的表面质量保证。现在，磨粒流体动力研磨技朮作为一种加工方法,适用于在内外部都需要高质量表面的产品。打破以传统手工研磨抛光工序: 适用于在内外部都需要高质量表面的产品例如:涡轮机内部零件/航天/汽车/各类 【】13精密工件。挤压工业:平面/分流/多孔/精细复杂模具。药用业/纺织业/液压/压缩/气动工件。 6 图1-7德国Perfect Finish Gmbh磨粒流流体动力研磨系统 SPKS挤压研磨流体抛光机械为微精处理机械，对于凹陷面与弯曲孔道等通常刀、磨具达不到的复杂形状优为有效，该技术打破了传统的手工研磨抛光工序，使微孔、多孔、长孔、弯孔、异形孔的工件抛光研磨便利、轻松，特别是在气体、液体类的导通管内进行镜面抛光，使研磨痕和流体通过方向一致，有效的提高模具或工件的性能、质量、光洁度，达到镜面等级，同时延长模具及工件的使用寿 【】15命，更能提升产品的品质和产能。 图1-8 SPKS 挤压研磨流体抛光机械 2. 国内的发展现状 我国很早就引进了磨粒流光整技术，并用于铝型材模具的抛光。引进的铝型材生产线上一般都配有磨粒流设备通称挤压研磨设备。国内有几家研究单位已将该工艺用于不同类型零件的光整加工，在夹具设计及加工控制方面积累了一些经验。已在航空、航天、汽车、纺机、模具等领域取得了可观的经济效益。但由于对这项光整技术宣传不够，很多人至今还不了解，在一些领域里处于空白状态。另外，由于该技术在国内应用尚不广泛，经验也不足，工艺上也不成熟，自动化程度不高，未能达到一定的生产规模，磨粒流加工的优势尚未充分发挥出来。 北京航空工艺研究所经过几年的探索和研究，现在已取得了可惜的进展。自 7 行开发的磨粒流介质，性能达到了国际先进水平，可以替代进口产品。 图1-9中航工业北京航空制造工程研究所 近年来，Fletcher等研究了磨料流加工中应用的高分子聚合物的热特性和流变性，认为介质的流变性对磨料流加工的成败具有重要的作用。Davies和Fletcher研究了几种配料的流变性与其相应加工参数之间的关系，结果表明粘度和磨料的比例都会影响温度和介质通过工件的压力下降，在磨料加工过程中温度是影响介 【】16质粘度的一个重要因素。Williams和Rajurkar的研究表明，介质的粘度和挤压力主要决定着表面的粗糙度和材料去除率，表面粗糙度精度的改善主要发生在磨料介质的前几个挤压往复进程中，并提出了估算动态有效切削磨粒数目的方法和每个行程中磨粒损伤的计算方法。他们还提出了多孔抛光中金属出去分布的实验方法与定量分析方法，发现用磨粒流加工一个具有中心孔和 四个外围孔的工件时，中心孔的金属去除率比外围孔的金属出去率高30%。Williams和Rajurkar研究探索了磨料流加工过程特性的一些方面和表面的特性化以及过程建模，研究了工艺输入参数对工艺性能参数的影响，利用DDS随机建模与分析技术研究了磨料流加工表面;磨料流加工表面轮廓模型的格林 函数 excel方差函数excelsd函数已知函数     2 f x m x mx m      2 1 4 2拉格朗日函数pdf函数公式下载 揭示了其特性形状是双指数的叠加。Williams等还研究提出了基于监控策略和磨料流加工的声发性特性的磨料流加工声发射在线监控和自适应控制系统，但是这些研究工作仅仅考虑了一部分过程参数而忽略了其他一些关键参数。磨料流加工去除了传统以及特种加工对工件表面的影响，使表面更加均匀一致;磨料流加工与磨削加工有很多相似之处。 在国内王纯、杨建明和王洁针对传统的磨料流加工在磨料介质流速增大的情 8 况下容易出现剪切变稀现象从而易失效的不足，研究开发出了磨料流振动抛光机床和相应的加工技术，从而有利于较大幅度提高抛光效率，并提出了磨料流加工流动新的边界条件假设，用简易实验验证了这种假设。汤勇等对磨料流加工存在的壁画滑动现象进行了实验研究，结果表明:磨料平均速度存在着临界值，在于平均速度时存在壁面滑动现象，同时壁滑速度随磨料平均速度增加而增大;磨料粘度升高，会使平均速度临界减少，而壁滑速度增加的程度却增大，存在壁滑是实现磨料流加工的前提条件。 1.4 小结 本节我们回顾了磨料流加工的形成及其加工特点，研究了其目的和意义以及国内外发展现状，使我们对其有了一个大致而清晰的了解。 9 第2章 Fluent软件 2.1 GAMBIT软件介绍 GAMBIT是面向CFD的几何建模和网格生成软件，是目前CFD分析中最优秀的前置处理器，它包括先进的几何建模和网格划分方法。用户既可以在GAMBIT中直接建立点、线、面、体的几何模型，也可以从PRO/E、UG、CATIA、SOLIDWORKS、ANSYS等主流的CAD/CAE软件中导入创建好的实体与网格。GAMBIT与CAD软件的接口和功能强大的布尔运算能力可使用户方便的建立复杂几何模型。借助其功能灵活、完全集成和易于操作的界面，GAMBIT软件可以显著减少CFD应用中前置处理的时间。复杂的模型可直接采用GAMBIT固有几何模块生成，或由CAD/CAE构型系统导入。高度自动化的网格生成工具保证其最佳质量的网格生成，如结构化的、非结构化的、多块的、混合的网格。GAMBIT可以生产FLUENT、FIDAP、POLYFLOW等求解器所需要的网格，它提供的非结构化的网格生成程序，对相对复杂的几何结构网格生成非常有效。FLUENT还可根据计算结果调整网格，这种网格的自适应能力对于精确求解有较大梯度的流场有很实际的作用。由于网格自适应和调整只是在需要加密的流动区域里实施，而非整个流场，因此可以节约计算时间。 GAMBIT是为了帮助分析者和设计者建立并网格化计算流体力学(CFD)模型和其它科学应用而设计的一个软件，GAMBIT通过它的用户界面(GUI)来接受用户的输入。GAMBIT GUI具有简单而又直接的做出建立模型、网格化模型、指定模型区域大小等功能。计算机网格根据性质的不同，可以分为结构化网格、非结构化网格和混合网格。对于二维平面系统模型而言，四边形网格属于结构网格，而三角形网格则属于非结构网格。三维立体系统模型则较复杂，有四面体、五面体、六面体和楔形等等立体结构网格。其中，六面立方体或六面长方体皆属于结构网格，其他则皆属于非结构网格。 因为流动域的几何模型比较复杂，喷油嘴零件的计算网格主要采用非结构化的三角形、四边形的面网格以及非结构化的四面体、五面体和六面体的体网格，在靠近壁面的流动区域采用了棱柱型的边界层网格，以便准确地模拟流动和换热情况。图2-1给出了GAMBIT中常用的网格类型。 相对四边形的面网格和六面体的体网格而言，在流动域内生成三角形的面网格和四面体的体网格更容易一些，而且网格数相对也少。但采用四边形面网格可以允许网格单元的变形更大一些，而三角形、四面体网格单元会因为网格的变形过大而增加网格的扭曲率，从而影响计算结果的收敛性和准确性。需要指出的是如果采用非结构化的三角形、四面体网格模型来求解，需要采用二阶求解精度。 10 为了提高求解精度，防止数值离散，对于相对简单的几何结构最好还是采用四边形、六面体网格来划分。而对于复杂的几何结构采用三角形、四面体网格来划分，此类网格容易生成网格单元并且可以减少网格数，比较适用于几何结构复杂的喷油嘴喷孔结构。喷油嘴二维及三维网格划分如图2 -2所示。 图2-1 GAMBIT常用网格类型 当影响网格大小的因素较多，主要有以下几个方面:几何实体的大小是决定网格尺寸的主要因素;理想的三角形网格是等边的，面网格的尺寸将影响体网格的大小进而影响边界层的求解;相邻结构的网格尺寸也会影响网格的大小，若与凹坑处相邻边界面的尺寸较大，那么凹坑处生成的网格只有一个体网格，这样该处的网格质量就会变差。无论采用何种方式来划分计算网格，网格的类型、尺度和质量都会直接影响到CFD分析结果的准确性和稳定性。 喷油嘴二维网格划分 喷油嘴三维网格划分 图2-2 喷油嘴模型网格划分 磨粒流加工过程属于复杂的湍流流动，复杂的湍流流动在传输平均动量和其它标量过程中起着主要作用，如果要实现对湍流流动的正确模拟，那么对靠近壁面流动区域划分计算网格是有要求的。由于湍流的平均流动和脉动之间存在着很强的相互作用，湍流的数值模拟计算结果比层流流动更容易受到网格的影响。为 11 了保证计算结果的准确，对于靠近壁面区域湍流的求解最好选择在平均流动变化快而且平均应力较大的边界层所在的区域。在计算网格的处理上，在流动域靠近壁面部分要采用加附面层的方法，在流体流动区域的表面网格和内部网格之间再划分出一个特定区域，用以保证壁面附近的湍流计算的准确性。 磨粒流加工过程属于复杂的湍流流动，复杂的湍流流动在传输平均动量和其它标量过程中起着主要作用，如果要实现对湍流流动的正确模拟，那么对靠近壁面流动区域划分计算网格是有要求的。由于湍流的平均流动和脉动之间存在着很强的相互作用，湍流的数值模拟计算结果比层流流动更容易受到网格的影响。为了保证计算结果的准确，对于靠近壁面区域湍流的求解最好选择在平均流动变化快而且平均应力较大的边界层所在的区域。在计算网格的处理上，在流动域靠近壁面部分要采用加附面层的方法，在流体流动区域的表面网格和内部网格之间再划分出一个特定区域，用以保证壁面附近的湍流计算的准确性。喷油嘴模型的网格划分如图2-3、2-4所示。 图2-3 喷油嘴二维网格划分 可用GAMBIT软件生成混合型四面体网格，在靠近壁面的边界层内生成一层很薄的三棱柱网格，在流体内部区域生成四面体网格。与在整个流动域全部生成四面体网格相比，此种方法在靠近壁面流动区域为棱柱型混合网格，可以更好的模拟流体接近壁面区域的流场的情况，计算结果也更为准确可靠。 图2-4 喷油嘴三维网格划分 需要指出的是，本章网格划分采用指定网格间距(interval size)模式，支路网格密度大于干路网格密度。当网格划分成功后，将创建成功的模型导出保存为*mesh文件，启动FLUENT读取此文件，选择恰当的仿真参数即可进行数值分析。 12 2.2 FLUENT软件概述 FLUENT是目前处于世界领先地位的商业CFD软件包之一，最初由FLUENT Inc.公司发行。2006年2月ANSYS Inc.公司收购FLUENT Inc.公司后成为全球最大的CAE软件公司。FLUENT6.3.26就是由ANSYS Inc.公司发布的新版本，CFD软件通常都包含3个主要功能部分:前处理器、求解器、后处理器。 FLUENT是一个用于模拟和分析复杂几何区域内的流体流动与传热现象的专用软件。FLUENT提供了灵活的网络热性，可以支持多种网络。用户可以自由选择使用结构化或者非结构化网络来划分复杂的几何区域，例如针对二维问题支持三角形网络或者四边形网络;针对三维问题支持四面体.六面体.凌锥.多面体网络;同时也支持混合网络。用户也可以利用FLUENT提供的网络自适应特性在求解过程中根据所获得的计算结果来优化网格。 FLUENT是使用C语言开发的，支持并行计算，支持UNIX和WINDOWS等多平台，采用用户.多服务器的结构，能够在安装不同操作系统的工作站和服务器之间协同完成同一个任务。FLUENT通过菜单界面和用户进行交互，用户可 【】11以通过多窗口的方式随时观察计算的进步和计算结果。 FLUENT软件包主要由GAMBIT、Tgrid、Fliters、FLUENT几部分组成。 (1)前处理器。包括GAMBIT、Tgrid和Fliters。其中GAMBIT是由FLUENT Inc.公司自主开发的专用CFD前置处理器，用于模拟对象的几何建模以及网格组成。Tgrid是一个附加前置处理器，它可以从GAMBIT或其他CAD/CAE软件包中读入所生成的模拟对象的几何结构，从现有的边界网格开始生成由三角形、四面体、混合网格组成的体网格。Fliters实际上就是其他CAD/CAE软件包，例如ANSYS、CGNS等于FLUENT之间的接口，通过接口可以将由其它CAD/CAE软件包所生成的面网格或体网格读入到FLUENT。 (2)求解器。它是CFD软件包的核心，FLUENT实际上是一个求解器，FLUENT6.3.26是一个基于非结构化网格的通用求解器，支持并行计算，分单精度和双精度两种。一旦所生成的网格读入到FLUENT中，所有剩下的操作都可以在FLUENT里面完成，其中包括设置边界条件、定义材料性质、执行求解、根据计算结果优化网格、对计算结果进行后处理等。 (3)后处理器。FLUENT本身就附带有强大的后处理功能，有云图、等值线图、矢量图、剖面图、XY散点图、粒子轨迹图、动画等多种方式显示、存储和输出计算结果，可以平移、缩放、旋转、镜像图像，也可以将计算结果导出到其他CFD、FEM软件或其他后处理软件中，例如 Tecplor。 FLUENT的强大的求解功能使其广泛应用于国防、航空航天、机器制造、汽车、船泊、兵器、电子、铁道、石油天然气、材料工程等领域。 利用FLUENT进行流体流动与传热的模拟的计算流程如图2-5所示。首先利 13 用GAMBIT进行流动区域几何形状的构建、边界类型以及网格的生成，并输出 用于FLUENT求解器计算的格式，然后利用FLUENT求解器对流动区域进行求 解计算，并进行计算结果的后处理。 图2-5 基本结构程序示意图 FLUENT程序的用途 1 采用三角形、四边形、四面体、六面体及其混合网格计算二维和三位流动 问题。计算过程中，网络的可以自适应。 2 可压缩与不可压缩流动问题 3 稳态和瞬态流动问题 4 无粘流，层流及湍流问题 5 牛顿流体以及非牛顿流体 6对流换热问题(包括自然对流和混合对流) 7 导热于对流热耦合问题 8 辐射换热 9惯性坐标系和非惯性坐标系下的流动问题模拟 10多运动坐标系下的流动问题 11 化学组成混合与反应 12 可以处理热量、质量、动量和化学组成的源项 13用lagrangian轨道模型模拟稀疏相 14 多孔介质流动 15 一维风扇、热交换器性能计算 16 两相流问题 14 17 复杂表面形状下的自由面流动 2.3 小结 在进行磨粒流数值分析前，必须首先应对被分析零件进行模型的创建与网格的划分工作，网格的质量将影响数值分析的精度。利用CFD前处理软件GAMBIT来进行模型的创建与网格划分工作，这样求解的结果就可以直接导入FLUENT软件中进行求解，最后达到对问题进行处理和分析的目的。 15 第三章 数学模型和计算方法 3.1数学模型 任何流体运动的规律都是以质量守恒定律、动量守恒定律和能量守恒定律为基础，这些定律可以由数学方程组来描述，如欧拉方程、N-S方程。采用数值计算方法，通过计算机求解这些控制流体流动的数学方程，进而研究流体的运动规律，这样的学科就是计算流体力学CFD。通过数值模拟，我们可以得到极其复杂的流场内各个位置上的基本物理量，如速度、压力、温度、浓度等物理量的分布云图，从中可以分析流体的运动规律，在一定程度上可以代替实验并节约资金。 磨粒流由碳化硅颗粒和流体介质按一定比例混合而成的，文中将流体和颗粒都看作连续介质，可据连续介质理论的质量守恒和动量守恒定律建立流体相和颗粒相的数学模型。 3.1.1 控制方程数学模型 对于所有的流动性问题，FLUENT都遵循质量守恒、动量守恒和能量守恒方程。 (1)质量守恒方程 ,(),u,,i (3-1) ，,Sm,,txi 该方程是质量守恒的总形式，可以适合可压和不可压流动。源项是稀疏相Sm增加到连续相中的质量，如液体蒸发变成气体或者质量源项。 对于二维轴对称几何条件，连续方程可以变成: ,,,,,,,()()uvv，，，,S (3-2) m,,,txrr 式中，x是轴向坐标，r是径向坐标，u和v分别是轴向和径向速度分量。 (2)动量守恒方程 在惯性坐标系下，i方向的动量守恒方程为: ,,(),,uu,()u,,pijiji (3-3) ，,,，，，gF,ii,,,,txxcjij ，，,uuu,,2jil,,,,,式中，p是静压，是应力张量，定义为:，(),，,,,ijijij,,,xxx3,,jil，， FF是两相之间作用的重力体积力和其他体积力，还可以包括其他模型源项或者ii 用户定义源项。 16 (3)能量守恒方程 能量守恒定律又称为热力学第一定律，能量方程可由下式来表达。 222,,,,,,,()()()uTvTTqTTT,,,,eff (3-4) ，，,，，，()222,,,,,,xyzCCxyzpp 式中，q是热源，为比热容，为有效导热系数。 C,peff 3.1.2 流体相数学模型 磨粒流加工时的流动属于湍流流动，这里选用柱坐标下轴对称物理模型，采k,,用 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 二方程湍流模型建立封闭的数学模型。 ,U1,(yV)连续性方程: (3-5) ，,0,xy,r 轴向动量方程: 2,,,,,,,1UPU1V,U,yUV,,U(,)1(,) ，,，，(y,)(,)(,,)，,(,)eee2,,,,,,,,xy,y,x,xyyyxxxyyx(3-6) 径向动量方程: 2,,,,,()1()1,,UVyVV2,V,,,,,,VPUV1 e，,，(),()()()yy,,,,，，,e2eee,,,,,xyyxxyy,,,,,,yyxyyyyy (3-7) 湍流动能k方程: ,,()1(),,UkyVk,,,,,,kk1tt [()][()],,,,，，，，,yP，,p,,,,xxyyy,,,,xyykk (3-8) 湍流动能耗散率方程: , ,,,,,,,,,1,,()1(),,,,UyVtt yCPC[][()]()，，，，,,,,,，,12Pxxyyyk,,,,,,,,xyy,, (3-9) ,,,,UVVUV2222其中: p,，，，，,{2[()()()]()}p1,,,,xyyyx 以上各式中，U、V分别为流体介质的轴向和径向平均速度;x、y分别为流 ,,体介质的轴向和径向坐标;为流体介质密度;P为流体介质静压;为动力粘 2,,,,,度系数;有效粘性系数=+;涡黏系数;、C、C、、,,,,Ck/Cet12k,,t, 是经验常数，在进行仿真运算时选择默认值即可。 3.1.3 颗粒相碳化硅数学模型 根据液固两相流理论，将碳化硅颗粒看作离散相，忽略颗粒与颗粒之间的作 用力，不考虑颗粒对磨料介质的影响，可建立离散相数学模型: 17 dUC3Re, (3-10) pDp,,()UU2epdtd4,pp 式中:U、U分别为航空煤油和碳化硅磨粒速度分量，d为碳化硅颗粒直径，epp C为阻力系数。 D 3.1.4 湍流数值模拟数学模型 流体实验表明，当雷诺数小于某一临界值时，流动是平滑的，相邻的流体层彼此有序的流动，这种流动称作层流;当雷诺数大于临界值时，会出现一系列复杂的变化，最终导致流动特性的本质性变化，流动呈无序的混乱状态。这时，即使是边界条件保持不变，流体也是不稳定的，速度等流动特性都随之变化，这种状态称为湍流。 一般认为无论湍流运动多么复杂，非稳态的连续方程和Navier-Stoke方程对于湍流的瞬时运动仍然是适用的。在此为了考虑不可压缩流动中脉动的影响，目前广泛采用的方法是时间平均法，即把湍流运动看作由两个流动叠加而成，一是时间平均流动，二是瞬时脉动流动。湍流流动的基本方程如下: ,,divu()0连续方程: (3-11) ，,,t, Stoke方程): 动量方程(Navier- 2，，,,,,,()()()(),,,,,uuuvuw,,,,, divuudivugradu()()，,,,，,,,,,,txxyz,,,,,,,，， 2，，,,,,,()()()(),,,,,vuvvvw,,,,, divvudivugradv()()，,,,，,,,,,,tyxyz,,,,,,,，， 2，，,,,,,()()()(),,,,,wuwvww,,,,, divwudivugradw()()，,,,，,,,,,,tzxyz,,,,,,,，， (3-12) 湍流流动是一种高度非线性的复杂流动，目前的湍流数值模拟方法可分为直接模拟和非直接模拟。所谓直接数值模拟就是直接求解瞬时湍流控制方程(4.12)和(4.13)。而非直接数值模拟就是不直接计算湍流的脉动特性，而是设法对湍流作某种程度的近似和简化处理。 3.2 计算方法 目前商业CFD软件的求解过程一般首先设定好边界条件和适当的物理参数，然后选择软件内假定好的物性模型，并经过适当的程序离散化成有限维度的近似问题，再经由计算机运算上述设定的条件，并将其近似解求出，最后由数学分析来确保所解出来的数值近似解的收敛性与稳定性。不同软件所用的数值方法存在 18 差异，目前主要分为有限差分法、有限元素法、有限体积法。 3.2.1 有限差分法 有限差分法是数值解法中最经典的算法。该方法是将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。有限差分法发展较早，比较成熟，较多用于求解抛物型问题，用该方法求解边界条件复杂、尤其是椭圆型问题不如有限元法或有限体积法有效。 3.2.2 有限元素法 有限元素法主要用来解决结构力学、流体动力学、热传与电磁场等工程上的问题。有限元素法一般用于解决复杂边界问题，该方法把偏微分方程的正解假想落在某一个适当的函数区间内，可以说是将其转换成每个元素区域的常微分方程或代数方程，最后用数值方法求解每个元素区域内的方程式。有限元素法的基础是极值原理和插值求解，该方法吸收了有限差分法中离散处理的内核，选择逼近函数并对其进行合理的积分。有限元素法适用于几何及物理条件比较复杂的问题，并且便于程序的标准化，对椭圆型方程问题更具有实用性。 3.2.3 有限体积法 有限体积法的基本思想是将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积，将待解的微分方程对每一个控制体积积分，便可解出一组离散方程，其中的未知数便是网格点上的因变量的数值。从积分区域的选取方面看来，有限体积法属于加权剩余法中的子区域法。从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。有限体积法运用于随时间变化而移动的不连续面时，其采用的方法是全隐式法。全隐式法可以允许任何时间差的变化，根据情况的不同，还应考虑最大时间间隔限制问题。对于暂态问题，其时间差距必须小到瞬间差;对于稳态问题，则需要较大的时间间隔才可以快速的达到稳定状态。本文的数值模拟计算方法采用有限体积法。 3.3 FLUENT求解方法的选择 FLUENT提供三种不同的解格式:非耦合求解 ( Segregated )、耦合隐式求解 ( Coupled Implicit )、耦合显式求解 ( Coupled Explicit )。 非耦合求解方法主要用于不可压缩或压缩性不强的流体流动。耦合求解则可以用在高速可压缩流动。FLUENT默认设置是非耦合求解，但对于高速可压流动，有强的体积力、浮力或离心力的流动，求解问题时网格要比较密，建议采用耦合隐式求解方法，可以耦合求解能量和动量方程，能比较快的得到收敛解。如果必须要耦合求解，但是计算机的内存较小，这时候可以考虑用耦合显式解法器求解问题，该解法器也耦合了动量，能量及组分方程，但内存占用却比隐式求解方法小，但收敛时间比较长。 19 这里需要指出的是非耦合求解的一些模型在耦合求解解法器里并不都有，耦合解法器没有的模型包括:多相流模型，混合分数/PDF燃烧模型，预混燃烧模型，污染物生成模型，相变模型，Rosseland辐射模型，确定质量流率的周期性流动模型及周期性换热模型等。磨粒流加工数值模拟属于压缩性不强的流体多相流模拟，故选择非耦合求解计算。 3.4 小结 本节介绍了数学模型和计算方法，我们知道了要对流体运动进行分析，都必须以质量守恒定律，动量守恒定律，动能守恒定律为基础，然后在这些方程组的基础上，进行流体运动的求解，进而研究流体的运动规律。通过数值模拟，我们可以得出流体在各个位置的云图，分心运动规律从而代替实验节约资金。并通过对FLUENGT的求解方法进行了阐述，对FLUENGT的正确求解有了深入了解。 20 第四章 湍流模型与边界条件 磨粒流加工过程的流动状态属于湍流加工，若获得正确的数值模拟结果，应首先设置正确的边界条件，下面简要介绍一下FLUENT软件所能处理的湍流流动及边界条件的设置问题。 4.1 湍流模型 湍流出现在速度发生变动的地方，这种波动使得流体介质之间相互交换动量、能量和浓度变化，引起了数量的波动。由于波动是小尺度且高频率的，因此在实际工程计算中直接利用计算机模拟对计算机的要求很高。实际上瞬时控制方程可能在时间上、空间上是均匀的，或者可以人为的改变尺度，修改后的方程可耗费较少的计算机内存。不幸的是没有一个湍流模型对于所有的问题是通用的，选择模型时主要依靠以下几点:流体是否可压、建立特殊的可行的问题、精度的要求、计算机的运算能力、时间的限制。为了选择最好的模型，需要了解不同模型的适用范围和限制，下面介绍几种常见的湍流模型。 k,,(1)标准模型 此模型是最简单的完整湍流模型，该模型是两个方程模型，需要求解湍动能及其耗散率两个变量，速度和长度尺度。湍动能输运方程是通过精确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的。 k,,标准模型需要求解湍流能及其耗散率方程，湍流能输运方程是通过精确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得 k,,到的。该模型假设流动为完全湍流，分子粘性的影响可以忽略。因此，标准模型只适合完全湍流的流动过程模拟。 k,,标准模型湍流能K和耗散率方程如下图所示: , ，，,,,Dkk,,t (4-1) GGY,，，，,,,,,,,,,,kbMDtxx,,,iki,,，， 2，，,,,D,,,,,,t (4-2) CGCGC,，，，,,,,，，,,,,,,,132kbDtxxkk,,,iki,,，， GG在上述方程中，表示由于平均速度梯度引起的湍流能产生，是由于福bk Y利影响的湍流能产生，可压速湍流脉动膨胀对总的耗散率的影响，湍流粘性M 2k,,系数=。 Ct,, 21 在FLUENT中，作为默认值常数，=1.44，=1.92，=0.09，湍流能CCC1,2,,k与耗散率的湍流普朗特数分别是=1.0，=1.3。 ,,,k, k,,(2) RNG 模型 k,,k,,RNG 模型来源于严格的技术统计，它和标准模型相似，但是有如下改进: (1) RNG模型在方程中加了一个条件，有效地改善了精度。 , (2) 考虑到了湍流漩涡，提高了这方面的精度。 k,,(3) RNG理论为湍流Prandtl数提供了一个解析公式，然而标准模式 使用的是用户提供的常数。 k,,(4) 标准模型是一种高雷诺数的模型，RNG理论提供了一个考虑低雷 诺数流动粘性的解析公式。 k,,这些公式的效用依靠正确的对待近壁区域RNG 模型是从暂态N-S方程 k,,中推导得出的。RNG 模型方程如下: ,,,,yk()()(),,,,,,kkuGGYS，,，，,,， ikeffkbMktxx,,,,ijxj 2,,,,,,, (4-3) ()()()()uCGCGCS，,，，,，,,,,,,,ieffkb,,,,,132txxxkk,,,,ijj G是由层流速度梯度而产生的湍流动能，G是由浮力而产生的湍流动能，YkbM由于在可压缩湍流中，过渡的扩散产生的波动，C，C，C，是常量，a和a是k123ke方程和e方程的湍流系数，S和S是用户定义的。 ke k,,在此模型中，利用RNG理论求得速度参数C等于0.0845，与标准准模型u k,,k,,的0.09很接近。这些特点使得RNG 模型比标准模型在更广泛的流动中有更高的可信度和精度。 k,,(3)带旋流修正的模型 k,,k,,带旋流修正的模型是近期出现的，比起标准模型来有两个主要的不同点: k,,(1) 带旋流修正的模型为湍流粘性增加了一个湍流公式。 (2) 为耗散率增加了新的传输方程，这个方程来源于一个为层流速度波动 而作的精确方程，模型要确保在雷诺压力中有数学约束。 k,,带旋流修正的模型直接的优势在于对平板和圆柱射流的发散比率有更精确的预测，并且它对旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。此模型对旋转均匀剪切流、包括喷射的自由流、混合流、管道和边界流、圆柱射流及分离流均有效。 22 (4) 标准 k-ω模型 标准k-ω模型是基于Wilcox k-ω模型，它是为考虑低雷诺数、可压缩性和剪切流传播而修改的。Wilcox k-ω模型预测了自由剪切流传播速率，像尾流、混合流动、平板绕流、圆柱绕流和放射状喷射，因而可以应用于墙壁束缚流动和自由剪切流动。标准k-ω模型是一种经验模型，是基于湍流能量方程和扩散速率方程。 由于k-ω模型已经修改多年，k方程和ω方程都增加了项，这样增加了模型的精度。标准k-ω模型的方程如下: ,,,,k()()(),, kkuGYS，,,，,，ikkkK,,,,txxxijj ,,,,,()()()uGYS，,,，,， (4-4) ,,,,i,,,,txxx,,,,ijj 在方程中，G是由层流速度梯度而产生的湍流动能，G是由ω方程产生的，Tkωk和T表明了k和ω的扩散率，Y和Y由于扩散产生的湍流，S和S是用户定义的。 ωkωke (6) The Spalart-Allmaras 模型 Spalart-Allmaras模型是相对简单的单方程模型，只需求解湍流粘性的输运方程，不需要求解当地剪切层厚度的长度尺度。该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果，在透平机械湍流模拟方面也有较好结果。由于没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合，如平板射流问题、从有壁面影响流动突然变化到自由剪切流、流场尺 Spalart-Allmaras度变化明显等场合。 模型常用于航空领域，来墙壁束缚流动，并且已经获得了很好的模拟结果。需要注意的是，Spalart-Allmaras 模型是一种新出现的模型，现在还不能断定它适用于所有的复杂的工程流体. k,,对于磨粒流加工过程的数值模拟，可选用标准模型或带旋流修正的k,,模型来进行。 4.2 边界条件 恰当的选择边界条件对使用CFD软件FLUENT进行磨粒流数值模拟的准确性是至关重要的，FLUENT提供很多边界条件设置流体进入或者流出求解域。主要有速度入口边界条件、压力入口边界条件、质量入口边界条件、压力出口边界条件、压力远场边界条件、自由出流边界条件、进口通风边界条件、进口风扇边界条件、出口通风边界条件、排气扇边界条件、对称边界条件、周期性边界条件和固壁边界条件等十余种边界条件。现简要介绍其中几种边界条件。 (1)速度入口边界条件 23 该边界条件常用于不可压缩流动问题，对可压缩问题不合适，需给定入口边界上的速度及其其他相关标量值。 (2)压力入口边界条件 该边界条件对可压缩与不可压缩问题都适用，边界条件需要输入总压、总温、流动方向、静压、湍流量、辐射参数、化学组分质量分数、第二相体积分数等参数。 (3) 质量入口边界条件 该边界条件主要用于可压缩流动。对于不可压缩流动，由于密度微常数，可选用速度入口边界条件。质量入口条件包括两种:质量流量和质量通量。质量流量是单位时间内通过总面积的质量，质量通量是单位时间单位面积内通过的质量。如果是二维轴对称问题，质量流量是单位时间单位内通过的质量。若想达到给定的流量，FLUENT会调解局部进口总压来调解速度，这样使得计算速度变慢。因此，如果压力入口边界条件和质量入口边界条件都合适时，应优先选用压力入口边界条件。 (4) 压力出口边界条件 该边界条件只能用于模拟亚音速流动，如果当地速度已经超过音速，该压力在计算过程中就不再适用了。在参数设置上，需给定表压强。该边界条件可以处理出口有回流的问题，对于有回流的出口，该边界条件比自由出流边界条件更容易收敛。 (5) 压力远场边界条件 该边界条件可以用于可压缩流动，常用于理性气体定律计算密度问题，如果知道来流的静压和马赫数，就可以用压力远场边界条件来模拟该类问题，为了满足无穷远压力边界条件，需要把边界放到关心的区域足够远的地方。 (6)自由出流边界条件 该边界条件常用于出流边界上的压力或速度均为未知的情形，此类边界条件的特点是不需要给出出口边界条件。但对于包含压力进口条件、可压缩流动问题、有密度变化的非稳定流动问题，此边界条件不适用。 (7)对称边界边界条件 该边界条件适用于流动及传热场是对称的情形，适用于粘性流中运动边界处理。在对称轴或对称面上，既无质量的交换，也无热量等其他物理量的交换。在对称轴或对称平面上，所有物理量在其垂直方向上的梯度为零，垂直于边界的速度分量也为零。计算过程中不需要给定任何参数，只需确定合理的对称位置即可。 (8)固壁边界边界条件 该边界条件用于壁面有平移运动或者旋转运动的情况，可以指定壁面切向速度分量，也可以指定壁面切应力从而模拟壁面滑移。根据流动情况，可以计算壁 24 面切应力和流体换热情况。对于粘性流动问题，FLUENT默认设置是壁面无滑移条件，利用标准函数壁面法可处理近壁面处的流动问题。 通过本文所设计的磨粒流加工装备可实时监测磨粒流加工不同零件的入口压力和出口压力，并且当地速度属于亚音速流动,因此磨粒流加工数值模拟的边界条件设置为压力入口和压力出口边界条件。 4.3 小结 本章介绍了湍流模型和边界条件、以及边界条件各自的使用范围和条件，有利于数值模拟获得准确值。 25 第五章 数值模拟算法与多相流的选择 对于管道、高压喷嘴等微小孔结构中的复杂流体问题, 可以利用FLUENT来求解,该软件提供了湍流方程,可模拟湍流的流动状态。利用FLUENT流体分析软件对共轨管及喷油嘴内部磨粒流的流动形态进行数值模拟，可形象直观地获得磨粒流加工特性。 5.1 求解器的选择 FLUENT软件提供了四种求解器，分别为二维单精度解法器、三维单精度解法器、二维双精度解法器、三维双精度解法器。由于磨粒流加工过程属于湍流流动，对二维模型进行数值模拟时，我们利用二维双精度解法器进行求解，对三维模型进行数值模拟时，我们利用三维双精度解法器进行数值模拟。 5.2 离散格式的选择 FLUENT可以选择每一个控制方程的对流项的离散格式，主要包含一阶迎风格式、二阶迎风格式和QUICK格式。当流动和网格成一条线时，选用一阶迎风离散格式。对于三角形和四面体网格，流动和网格不成一条线，此时选用二阶迎风格式可获取高精度的结果。对于有复杂流动的四边形或六面体网格，如果使用二阶离散格式，可获得较好的仿真结果。 5.3 压力速度耦合方法的选择 在分离求解器中，FLUENT提供了SIMPLE、SIMPLEC以及PISO三种压力速度耦合算法。定常计算一般选用SIMPLE或SIMPLEC方法，对于过渡计算推荐使用PISO方法。PISO方法还可以用于高度倾斜网格的定常状态计算和过渡计算。需要注意的是压力速度耦合只用于分离求解器，对于耦合求解器不可以使用。 对于没有附加模型激活的层流流动的简单问题，收敛性已经被压力速度耦合所限制，通常可以用SIMPLEC算法很快得到收敛解。对于包含湍流或附加物理模型的复杂流动，只要有压力速度耦合做限制，SIMPLEC方法就可以提供其收敛性。 5.4 多相流仿真的选择 混合模型是一种简化的多相流模型，它用于模拟各相有不同速度的多相流，但是假定了在短空间尺度上局部的平衡，也用于模拟有强烈耦合的各向同性多相流和 26 各相以相同速度运动的多相流。 本文所进行的磨粒流加工属于多相流问题，在数值模拟工程中选用二阶迎风格式、SIMPLE压力速度耦合算法、多相流混合模型进行仿真。 5.5 小结 本文介绍了数值模拟的算法和多相流的选择，这样对于复杂流体的计算有了充分的依据。 27 第六章 喷油嘴数值模拟 6.1 喷油嘴二维数值模拟 喷油嘴同样作为发动机系统的重要工作部件，喷油嘴内部的微小孔抛光及去毛刺问题一直是有待解决与提高的问题。在进行数值模拟的过程中，同样假设颗粒相与介质相之间不发生溶解或结晶等化学过程;颗粒状固体磨料与介质具有相同的压力作用;颗粒状固体磨料与介质分别满足动量守恒和能量守恒方程;颗粒状固体磨料与介质之间的相互作用通过阻力系数来实现的。 根据企业提供的数据，我们选用微小孔直径为0.16mm的喷油嘴来进行数值模拟。由于喷油嘴上的五个微小喷油孔是均匀分布的，在不影响分析结果的前提下，将其简化成二维轴对称几何结构。 当参数设置完成后，采用SIMPLE算法求解两相流动方程，经过多次计算，对于喷油嘴二维数值模拟大约经过277步迭代就可得到理想的精度。计算结束后，获得如下图6-1所示的残差监测变化曲线。 图6-1 残差监视变化曲线 从残差曲线可以看出，在迭代277次左右开始收敛，说明模型结构设计和求解参数设置比较合理，能得到满意的收敛解。 需要指出的是，由于喷油嘴内腔结构相对复杂，所以在建立喷油嘴二维物理 28 模型时，首先在AUTOCAD软件中建立二维模型，之后将此模型保存为.sat格式输出。在GAMBIT中将其导入后进行面域及网格的划分，并设置边界条件。当完成边界条件设置后，保存为.mesh文件，再将此文件载入FLUENT中即可进行数值模拟。喷油嘴磨粒流数值模拟参数与共轨管磨粒流数值模拟参数设置相同。图6-2为不同进出口压力条件下的静态压强云图。 a)进口压力2 Mpa，出口1Mpa b)进口压力3 Mpa，出口2Mpa c)进口压力5 Mpa，出口3Mpa d)进口压力7 Mpa，出口5Mpa e)进口压力8Mpa，出口5Mpa f)进口压力9Mpa，出口6Mpa 图6-2 静态压强云图 从数值模拟结果来看，静态压强从入口到出口呈梯度递减趋势，入口静态压强最大，出口静态压强最小，符合喷油嘴磨粒流加工数值模拟的参数设置。 29 a)进口压力2 Mpa，出口1Mpa b)进口压力3 Mpa，出口2Mpa c)进口压力5 Mpa，出口3Mpa d)进口压力7 Mpa，出口5Mpa e)进口压力8Mpa，出口5Mpa f)进口压力9Mpa，出口6Mpa 图6-3 动态压强云图 从动态压强仿真结果来看，在喷油嘴微小孔处的动态压强远大于喷油嘴内腔的动态压强，磨粒流在喷油嘴微小孔处的运动更为激烈，加工效率将越高，能更好的保证喷油嘴微小孔处的加工精度。由于喷油嘴整个内腔与喷油嘴喷孔的直径与长度的相差较大，单位面积上的作用力将有很大的差别，作用于喷油嘴微小孔处的作用力更大，这样更加有利于喷油嘴喷孔处的精抛加工。图6-4所示的喷油嘴速度云图更能体现上述观点。 30 a)进口压力2 Mpa，出口1Mpa b)进口压力3 Mpa，出口2Mpa c)进口压力5 Mpa，出口3Mpa d)进口压力7 Mpa，出口5Mpa e)进口压力8Mpa，出口5Mpa f)进口压力9Mpa，出口6Mpa 图6-4速度云图 通过比较不同加工条件下的速度云图可以发现，磨粒流在喷油嘴微小孔处的速度最大，在喷油嘴内腔的壁面速度最小，仅存在很小的相对滑移，这样更加有利于磨粒流对喷油嘴微小孔的孔壁进行抛光。根据经验可知，运动速度越大，磨粒流中的磨粒与孔壁的接触机会就越大，金属表面的去除量也越大，就更有利于磨粒流对微小孔结构的精加工。需要说明的是，随着进口压力的不断增加，磨粒流加工入口边界条件的动压设置不断增大。从速度云图中还可以看出，在相同进出口压力差的条件下，动压设置越大，相对应的速度越小。图6-5给出了喷油嘴磨粒流加工数值模拟的有效粘度云图。 31 a)进口压力2 Mpa，出口1Mpa b)进口压力3 Mpa，出口2Mpa c)进口压力5 Mpa，出口3Mpa d)进口压力7 Mpa，出口5Mpa e)进口压力8Mpa，出口5Mpa f)进口压力9Mpa，出口6Mpa 图6-5 有效粘度云图 通过观察喷油嘴磨粒流加工有效粘度云图可以知道，在喷油嘴微小孔处的有效粘度最大，内腔孔壁的有效粘度次之，内腔中心处的有效粘度最小。有效粘度越大，磨粒流加工作用越明显，越有利于微小孔流道的表面精加工。经过观察还可以知道，流体的有效粘度只于压差有关，与磨粒流的初始动能无关，相同的压差会得到等值的有效粘度。 32 6.2 小结 通过分析磨粒流加工喷油嘴微小孔流道的数值模拟结果来看，利用磨粒流加工喷油嘴确实可以在喷孔处获得理想的加工效果，提高喷油嘴的喷射性能，进而提高发动机的燃油性能及工作效率。而且利用磨粒流加工喷油嘴零件是一种行之有效的加工方法。利用磨粒流加工将会在喷油嘴喷孔处获得理想的加工精度，获得极好的喷射效率，使柴油在发动机内获得最佳的燃烧性能。 33 结 论 通过对磨粒流加工原理和运动规律的分析，探讨喷油嘴微小孔磨粒流加工特 k,,性。以流体力学软件FLUENT为平台，采用固液两相Mixture湍流模型对磨粒流加工过程中磨粒流的流动形态进行数值模拟。 根据磨粒流加工特点和喷油嘴零件结构建立二维仿真模型，应用FLUENT软件对磨粒流加工过程进行数值模拟。分析了静态压强、动态压强、湍流速度、以及湍流粘度对磨粒流加工的影响。 从模拟结果可知，通过改变入口压力可获得不同的磨粒流加工速度。对不同工艺的磨粒流加工作了比较，从模拟结果可以看出，随着压力的增大，喷油嘴微小孔出口处的磨粒流速度也随之增大，湍流粘度降低，磨粒切削能力增强，从而提高磨粒流加工效率。在其他条件不变的情况下，压力越大，磨粒的切削能力越强，磨粒流对微小孔流道的加工效果越好。为了获得较好的圆角和理想的孔壁粗糙度，应适当增加入口压力，选择较小运动粘度的磨粒流介质。 根据速度分布图可以看出，在孔交叉处的速度较大。由于磨粒流的挤压研磨作用，会在此区域产生较大的圆角，能够有效降低高脉冲产生的疲劳应力，明显改善零件的疲劳强度。用该加工方法能使表面光整，增强零件可靠性，延长使用寿命。 通过分析不同条件下的流场内部压力、速度、粘度等因素的数值模拟图对磨粒流加工的参数化有一定意义，为磨粒流加工参数的优化选择提供了理论依据。 34 参考文献 [1] 毕建忠，李继亮，磨料流加工技术及应用[J ]，工程设计与应用研究，2001(3);7~11。[2] 罗华.磨粒流加工在生产中的发.技术报告，2004。 [3] 艾冬梅,贾志新. 小孔加工技术发展现状[J ] . 机械工程师,2000 (1) :8-10。 [4] 贾宝贤,王振龙,赵万生. 基于特种加工的微小孔加工技术[J ] . 电加工与模具,2005 (2) :1-5。 [5] 聂先桥，杨建民，磨料流加工磨料流动边界的确定[J ].机械制造自动化，2001，446(49) 9-11。 [6] 汪文友，粘性磨料流加工机理的探究，磨料磨具与磨削，2001，(3):17~21。 [7] 戴根村.磨料流技术去除交叉孔毛刺[J ].航天制造技术，2002(1);25~28。 [8] 李元宗.磨料流加工[J ] ，电子技术，2002(4);10~12 [9] 杨贵铭.磨粒流挤压技术研究应用[J ].航空制造工程，2005(11);12~13 [10]李融 磨料流加工机器应用范围[J ]，机床，1999(5);21~22。 [11] FLUENT入门与进阶教程 北京理工大学出版。 [12] L.J. Rhoades，Abrasive flow machining，Manuf .Eng.2002(101); 75-78。 [13] R.K.Jain，V.K.Jain，Modeling of material removal and surface roughess in abrasive flow machining process，Int.J.Mach.Tools Manuf2007(39);1903-1923。 [14] R.E.Williams，K.P.Rajurkar，Stochastic modeling and analysis of abrasive flow machining， ASME J.Eng.Ind.2003(114); 74-81. [15] P.Koshy，V.K.Jain，G.KLal，Stochastic simulation approach to modeling diamond wheel topography，Int.J.Mach.Tools Manuf.2005(37); 751-761。 [16] H.Yoshikawa，T.Sata，Simulated grinding process by Monte Carlo method，Ann.CIRP 1999 (16 ):297-302. 35 致谢 这次毕业论文是在导师的亲切关怀和悉心指导下完成的，首先我要对导师在毕业设计过程中给予的谆谆教诲致以深深的谢意～ 从毕业设计选题到设计完成，在短短的几个月的毕业设计期间，李老师给予了我耐心的指导和细心关怀。正是在李老师耐心的指导和关怀下我才不会在设计的过程中迷失方向，失去前进的动力。李老师严谨的治学态度，求真务实的工作作风，注重学生的思维，拓展学习知识面以及谦和的为人风范对我的学习产生了极大的影响，同时这些对我今后的工作和学习也具有深远的意义。还要感谢在这次毕业设计中帮助我的学长，在做毕业设计的困难时期给予我很大的帮助。 最后，感谢关心我支持我的朋友，感谢学校领导、老师们，感谢你们给予我的帮助和关怀，感谢长春理工大学以及学院给予我四年良好的学习环境，谢谢～ 36 学位论文原创性声明 本人郑重声明:所呈交的学位论文～是本人在导师的指导下进行的研究工作所取得的成果。尽我所知～除文中已经特别注明引用的内容和致谢的地方外～本论文不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出重要贡献的个人和集体～均已在文中以明确方式注明并表示感谢。本人完全意识到本声明的法律结果由本人承担。 学位论文作者,本人签名,: 年 月 日 学位论文出版授权书 本人及导师完全同意《中国博士学位论文全文数据库出版章程》、《中国优秀硕士学位论文全文数据库出版章程》(以下简称“章程”)～愿意将本人的学位论文提交“中国学术期刊,光盘版,电子杂志社”在《中国博士学位论文全文数据库》、《中国优秀硕士学位论文全文数据库》中全文发表和以电子、网络形式公开出版～并同意编入CNKI《中国知识资源总库》～在《中国博硕士学位论文评价数据库》中使用和在互联网上传播～同意按“章程”规定享受相关权益。 论文密级: ?公开 ?保密(___年__月至__年__月)(保密的学位论文在解密后应遵守此协议) 作者签名:_______ 导师签名:_______ _______年_____月_____日 _______年_____月_____日 37 独 创 声 明 本人郑重声明:所呈交的毕业设计(论文)，是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议。尽我所知，除文中已经注明引用的内容外，本设计(论文)不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。 本声明的法律后果由本人承担。 作者签名: 二〇一〇年九月二十日 毕业设计(论文)使用授权声明 本人完全了解滨州学院关于收集、保存、使用毕业设计(论文)的规定。 本人愿意按照学校要求提交学位论文的印刷本和电子版，同意学校保存学位论文的印刷本和电子版，或采用影印、数字化或其它复制手段保存设计(论文);同意学校在不以营利为目的的前提下，建立目录检索与阅览服务系统，公布设计(论文)的部分或全部内容，允许他人依法合理使用。 (保密论文在解密后遵守此规定) 作者签名: 二〇一〇年九月二十日 38 致 谢 时间飞逝～大学的学习生活很快就要过去～在这四年的学习生活中～收获了很多～而这些成绩的取得是和一直关心帮助我的人分不开的。 首先非常感谢学校开设这个课题～为本人日后从事计算机方面的工作提供了经验～奠定了基础。本次毕业设计大概持续了半年～现在终于到结尾了。本次毕业设计是对我大学四年学习下来最好的检验。经过这次毕业设计～我的能力有了很大的提高～比如操作能力、分析问题的能力、合作精神、严谨的工作作风等方方面面都有很大的进步。这期间凝聚了很多人的心血～在此我表示由衷的感谢。没有他们的帮助～我将无法顺利完成这次设计。 首先～我要特别感谢我的知道郭谦功老师对我的悉心指导～在我的论文书写及设计过程中给了我大量的帮助和指导～为我理清了设计思路和操作方法～并对我所做的课题提出了有效的改进方案。郭谦功老师渊博的知识、严谨的作风和诲人不倦的态度给我留下了深刻的印象。从他身上～我学到了许多能受益终生的东西。再次对周巍老师表示衷心的感谢。 其次～我要感谢大学四年中所有的任课老师和辅导员在学习期间对我的严格要求～感谢他们对我学习上和生活上的帮助～使我了解了许多专业知识和为人的道理～能够在今后的生活道路上有继续奋斗的力量。 另外～我还要感谢大学四年和我一起走过的同学朋友对我的关心与支持～与他们一起学习、生活～让我在大学期间生活的很充实～给我留下了很多难忘的回忆。 最后～我要感谢我的父母对我的关系和理解～如果没有他们在我的学习生涯中的无私奉献和默默支持～我将无法顺利完成今天的学业。 四年的大学生活就快走入尾声～我们的校园生活就要划上句号～心中是无尽的难舍与眷恋。从这里走出～对我的人生来说～将是踏上一个新的征程～要把所学的知识应用到实际工作中去。 回首四年～取得了些许成绩～生活中有快乐也有艰辛。感谢老师四年来对我孜孜不倦的教诲～对我成长的关心和爱护。 学友情深～情同兄妹。四年的风风雨雨～我们一同走过～充满着关爱～给我留下了值得珍藏的最美好的记忆。 在我的十几年求学历程里～离不开父母的鼓励和支持～是他们辛勤的劳作～无私的付出～为我创造良好的学习条件～我才能顺利完成完成学业～感激他们一直以来对我的抚养与培育。 最后～我要特别感谢我的导师赵达睿老师、和研究生助教熊伟丽老师。是他们在我毕业的最后关头给了我们巨大的帮助与鼓励～给了我很多解决问题的思路～在此表示衷心的感激。老师们认真负责的工作态度～严谨的治学精神和深厚的理论水平都使我收益匪浅。他无论在理论上还是在实践中～都给与我很大的帮助～使我得到不少的提高这对于我以后的工作和学习都有一种巨大的帮助～感谢他耐心的辅导。在论文的撰写过程中老师们给予我很大的帮助～帮助解决了不少的难点～使得论文能够及时完成～这里一并表示真诚的感谢。 39 学位论文原创性声明 本人郑重声明:所呈交的学位论文～是本人在导师的指导下进行的研究工作所取得的成果。尽我所知～除文中已经特别注明引用的内容和致谢的地方外～本论文不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出重要贡献的个人和集体～均已在文中以明确方式注明并表示感谢。本人完全意识到本声明的法律结果由本人承担。 学位论文作者,本人签名,: 年 月 日 学位论文出版授权书 本人及导师完全同意《中国博士学位论文全文数据库出版章程》、《中国优秀硕士学位论文全文数据库出版章程》(以下简称“章程”)～愿意将本人的学位论文提交“中国学术期刊,光盘版,电子杂志社”在《中国博士学位论文全文数据库》、《中国优秀硕士学位论文全文数据库》中全文发表和以电子、网络形式公开出版～并同意编入CNKI《中国知识资源总库》～在《中国博硕士学位论文评价数据库》中使用和在互联网上传播～同意按“章程”规定享受相关权益。 论文密级: ?公开 ?保密(___年__月至__年__月)(保密的学位论文在解密后应遵守此协议) 作者签名:_______ 导师签名:_______ _______年_____月_____日 _______年_____月_____日 40 独 创 声 明 本人郑重声明:所呈交的毕业设计(论文)，是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议。尽我所知，除文中已经注明引用的内容外，本设计(论文)不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。 本声明的法律后果由本人承担。 作者签名: 二〇一〇年九月二十日 毕业设计(论文)使用授权声明 本人完全了解滨州学院关于收集、保存、使用毕业设计(论文)的规定。 本人愿意按照学校要求提交学位论文的印刷本和电子版，同意学校保存学位论文的印刷本和电子版，或采用影印、数字化或其它复制手段保存设计(论文);同意学校在不以营利为目的的前提下，建立目录检索与阅览服务系统，公布设计(论文)的部分或全部内容，允许他人依法合理使用。 (保密论文在解密后遵守此规定) 作者签名: 二〇一〇年九月二十日 41 致 谢 时间飞逝～大学的学习生活很快就要过去～在这四年的学习生活中～收获了很多～而这些成绩的取得是和一直关心帮助我的人分不开的。 首先非常感谢学校开设这个课题～为本人日后从事计算机方面的工作提供了经验～奠定了基础。本次毕业设计大概持续了半年～现在终于到结尾了。本次毕业设计是对我大学四年学习下来最好的检验。经过这次毕业设计～我的能力有了很大的提高～比如操作能力、分析问题的能力、合作精神、严谨的工作作风等方方面面都有很大的进步。这期间凝聚了很多人的心血～在此我表示由衷的感谢。没有他们的帮助～我将无法顺利完成这次设计。 首先～我要特别感谢我的知道郭谦功老师对我的悉心指导～在我的论文书写及设计过程中给了我大量的帮助和指导～为我理清了设计思路和操作方法～并对我所做的课题提出了有效的改进方案。郭谦功老师渊博的知识、严谨的作风和诲人不倦的态度给我留下了深刻的印象。从他身上～我学到了许多能受益终生的东西。再次对周巍老师表示衷心的感谢。 其次～我要感谢大学四年中所有的任课老师和辅导员在学习期间对我的严格要求～感谢他们对我学习上和生活上的帮助～使我了解了许多专业知识和为人的道理～能够在今后的生活道路上有继续奋斗的力量。 另外～我还要感谢大学四年和我一起走过的同学朋友对我的关心与支持～与他们一起学习、生活～让我在大学期间生活的很充实～给我留下了很多难忘的回忆。 最后～我要感谢我的父母对我的关系和理解～如果没有他们在我的学习生涯中的无私奉献和默默支持～我将无法顺利完成今天的学业。 四年的大学生活就快走入尾声～我们的校园生活就要划上句号～心中是无尽的难舍与眷恋。从这里走出～对我的人生来说～将是踏上一个新的征程～要把所学的知识应用到实际工作中去。 回首四年～取得了些许成绩～生活中有快乐也有艰辛。感谢老师四年来对我孜孜不倦的教诲～对我成长的关心和爱护。 学友情深～情同兄妹。四年的风风雨雨～我们一同走过～充满着关爱～给我留下了值得珍藏的最美好的记忆。 在我的十几年求学历程里～离不开父母的鼓励和支持～是他们辛勤的劳作～无私的付出～为我创造良好的学习条件～我才能顺利完成完成学业～感激他们一直以来对我的抚养与培育。 最后～我要特别感谢我的导师赵达睿老师、和研究生助教熊伟丽老师。是他们在我毕业的最后关头给了我们巨大的帮助与鼓励～给了我很多解决问题的思路～在此表示衷心的感激。老师们认真负责的工作态度～严谨的治学精神和深厚的理论水平都使我收益匪浅。他无论在理论上还是在实践中～都给与我很大的帮助～使我得到不少的提高这对于我以后的工作和学习都有一种巨大的帮助～感谢他耐心的辅导。在论文的撰写过程中老师们给予我很大的帮助～帮助解决了不少的难点～使得论文能够及时完成～这里一并表示真诚的感谢。 内部资料， 42 请勿外传～ 43