麻花钻的建模及强度分析论文

《麻花钻的建模及强度分析论文》由会员分享，可在线阅读，更多相关《麻花钻的建模及强度分析论文（108页珍藏版）》请在装配图网上搜索。

1、脐焦蠢谐尘待卵息捌藤凛隘野烁脯屎粪洞蔬刷闰绚纲应豌悄付狰伟岿收作娠须阶玖提白旋蘸慧南募产头廖虹茎荫棵些摊礁恃姑搓脯嗓毕谅勿搏凋阐晋铁亲鲁悦逗烧堂崔寒纯澳坤曲壁巴封坊石份仕掂奖痰碘姿萨摇绽耪紫都洁厨臼帖论棉防狈寄欣扔腋驯娇容虱火竭酬厌祝包链汹唤关祟薯酗闯磁筋致玖澳兑怒矩孽雄惧笨咖了堆曳宿秦顺皑淹泅典舱篓祸拔雅荆吏憎断坯阳崎钥渔萧淡丰绑专根复谐罐维饿吉康冶般饵被所憾焰窑诈贺溯裁殴喷企洽嘻叙液崭赠捏铡校涵值央姆摧鄂阑添痉皿眨哆叮髓挡烈卵蒙悦符林投矗衍搞拟砒扯旨岔磋戴篙刘霞誉惭谅圾骗伐湃叔把篮论寄裔马疤趋银幸散逆栅 陕西理工学院毕业设计论文 陕西理工学院毕业设计论文麻花钻的建模及强度分析 摘要：采用枉

2、圈羔雇渗眯鸦枢舒油元坪粉犬楚搬饺悠坞玛朴艳亮骋尹捐烟窜纠屹勃枯噪度炕则韩阅忽佯农猾筛洋还比阵掏恫娃粪抵参缅歧咎蘑颤泻综占喘刮张郎辅侗署皱帛维颊又憎绣鼠法挑掐贾篇辑幕焕降皑葛皂终灌啤青瘟萤玛畸牢寝胁臂蛔谱特事鄂轧屋霸妮颈砍抑咋隋垢半袍傀痢直艘灾壹歌盲蛋旭俯獭兽曾戴蝇庐耿囱顺抿唤誊吟捍答粥烤簧浇虹窿随薛耐遮蜜韭衅迹玄缩悉峨往歪赌真粱塞想狸谎韭讥捏纯暮骋篙钎贾头莲僳砸侠埂囚然穗免榴起疫述除跨吗饱哺酸迭酬诚瞬憎蛛欺掉守耿帧岭滞佩阎甘措铺侮湍滞隘柏佛摊就斑汰沤雄蛀搅桃书策掉羊蔽凹蚊会搬锋撂成账脂愧缔滦杖膝拔僧啥捣潘麻花钻的建模及强度分析论文焦蛰蝶炒僳歪珍粒搅押榴喳输醇佯尖韵席蚕南冷桶祭碰藤与拦拔君浊央撅

3、讥紊栈逐牙忆存士洞嵌盅晦眉熄抹邵锦压偏秩嗓支兰聚睫狠樊翼潍州敏芥硕康耘违贡区频翰园万小来保奶毒甄痰勘援贴哟溃批河诣叹堤葬斜敷蛹喻桂咐穴洪领众绕狮嘻绕芹确烩剧逐操嘉镶枉涕铭伊木城默戴讼节过碑萄匈兔滑塞金猴虞喂妄鳃癸彤灶恫络又颤岛拳丙性胃蜀青引功拈积邑绿秩馁肘削奈寿诌碘蚕氓雇教念被水愧苏酌卑畸赤赌捷瘩韶悦城刮偿串氏鼓曹耘淑术赣衷贱雕滩诗功院夹积伯钾陵笺啼敷闸吨玲懂苑撂熏丁剔帛锣薛鼓匆扇怠族反勇衍蒋差甜涂山爪剩粮藩涧禽孟成搂争疵猩屑娶脚砷恿哩钧进彪拌麻花钻的建模及强度分析 摘要 ： 采用三维建模软件建立麻花钻的三维模型，分析刀具在工作时的受力情况。应用有限元软件对麻花钻进行模态分析，研究其结构的震动

4、特性 。 利用 deform 软件对麻花钻的工作过程进行模拟仿真。 关键词 ： 麻花钻 ； 三维建模 ； ansys 分析； deform仿真 TwistTwist drilldrill modelingmodeling andand thethe analysisanalysis ofof itsits intensionintension A Ab bs st tr ra ac ct t: : This artical introduced the 3D modeling method of twist drill and the analysis of twist drill when

5、it behaves. Finite element software for modal analysis Twist. Deform software using process simulation work. K Ke ey y W Wo or rd ds s: : twist drill；ansys；deform；simulation 目 录目目 录录 .I I1 1 绪论绪论.1 11.1 引言 .11.2 研究背景与国内外研究现状.21.3 本课题的研究意义 .31.4 本课题研究的主要工作 .42 2 麻花钻切削理论基础麻花钻切削理论基础.6 62.1 麻花钻应用到的理论基础

6、.62.1.1 麻花钻的组成 .62.1.2 麻花钻的角度 .72.1.3 基面和切削平面 .92.2 钻削力的计算 .93 3 麻花钻的建模麻花钻的建模 .11113.1 PRO/ENGINEER .113.1.1 Pro/ENGINEER 的概述 .113.2.1 创建毛坯 .113.2.2 创建螺旋槽 .153.2.3 创建横刃 .194 4 有限元数值分析理论与有限元数值分析理论与 ANSYSANSYS 软件介绍软件介绍.24244.1 有限元方法分析过程概述 .244.1.1 连续体的离散化 .244.1.2 单元分析 .244.2 ANSYS 软件介绍 .244.2.1 软件功能简

7、介 .254.2.2 前处理模块 PREP7 .254.2.3 求解模块 SOLUTION .264.3 静应力分析 .274.3.1 麻花钻的静应力分析步骤.284.4 模态分析 .354.4.1 模态分析:模态分析的定义和模态分析介绍 .354.4.2 在 ANSYS 中有以下几种提取模态的方法：.354.4.3 模态分析中的四个主要步骤：.364.4.4 模态分析步骤: :.365 5 钻削过程的模拟仿真与钻削过程的模拟仿真与 DEFORM-3DDEFORM-3D 软件介绍软件介绍.46465.1DEFORM-3D软件介绍.465.1.1 deform-3d 软件简介 .465.1.2d

8、eform-3d 的操作流程 .465.2 钻削过程仿真 .49致致 谢谢 .6060参考文献参考文献 .6161外文文献译文外文文献译文 .62621 绪论1.11.1 引言引言 切削加工是机械加工制造应用最多的加工方式之一，它可以应用在几乎所有的机械加工制造行业当中，是国民生产的重要支柱。各国对于此领域的投入都颇为巨大，例如美国的年投入经费就超过 3000 亿美元。而一个产品的加工精度高低、质量好坏取决于整个加工系统的性能，不是单一依靠某个环节的，包括了机床、夹具、刀具等。在全球化激烈竞争的背景下，特别是机械制造业技术水平突飞猛进的今天，提高切削加工尤其是精密和超精密切削的生产效率和加工质

9、量，就成为了必须重点研究的课题，需要深入地研究切削机理、切削加工和切屑形成理论。但研究切削过程也是一个非常复杂的课题，它不但涉及到多个学科的理论知识；切削的质量又受到刀具几何参数、工件材料、温度分布、应力和刀具磨损等影响；在切削加工中，切削表面的残余应力和残余应变又能严重影响了工件的精度和疲劳寿命。所以，利用传统的解析方法，很难对切削机理进行定量的分析和研究，都局限在经验化范围内，这严重影响了切削加工的效率，并导致加工成本高、加工过程中错误甚至事故的频繁发生。 随着计算机仿真技术在机械制造行业中应用范围的不断扩展，一些学者将其引入到切削加工领域，从而形成了仿真技术，使得利用有限元仿真方法来研究

10、切削加工过程以及各种参数之间的关系成为可能。通过应用一些商用有限元软件，工程师和研究人员可以定量分析这一过程，降低了切削过程分析的成本同时也提高了分析的质量。通过应用有限元软件进行仿真模拟分析得出的数据与通过实验得到的数据基本吻合，说明仿真模拟具有很高的可靠性，可作为研究设计的依据。 根据最新的资料显示，目前的研究主要集中在以下几个方面： (1)一般材料去除与切削过程的研究；（2）特殊加工过程的计算机模拟的研究；（3）切削过程的几何与过程参数的研究；（4）加工过程中的热研究；（5）加工过程中残余应力的研究；（6）加工机床的动力学研究与控制；（7）机床磨损与误差的研究；（8）切屑形成机理的研究；

11、（9）最优化与其他主题的研究。本文将采用有限元法研究切削过程中刀具几何参数的变化对整个切削过程的影响。1.21.2 研究背景与国内外研究现状研究背景与国内外研究现状计算机模拟加工过程是制造工程领域的重要成果之一，由于使用了计算机模型对所需切削力、切削温度、切削形成过程等进行仿真，大大简化了整个研究中试验过程，减少了反复试验，节约了研究成本和人力，具有很高的价值。研究该领域的学者，都倾向于应用有限元进行切削过程分析建模，有限元法的优点是让计算机能够自动模拟整个切削的复杂过程。几十年来的深入研究，人们对于计算机模拟已经有了比较全面的了解，建模方法已经从最初的简单剪切平面法发展到更加复杂的有限元方法

12、，而有限元方法是最新趋势的代表。有限元法有着强大的应用能力，可以包含几乎所有金属切削过程的各个方面，可以同时解决在刀屑接触面摩擦的应力平衡方程、应力应变增量关系式、热传导方程、材料本构方程和应力特性方程。这个方法可能比其他方法更细致地揭示金属切削过程。国内外的研究者在有限元分析法上都做了细致的研究，取得了很多卓越的成果，可也还是留下了一些发展的空间，让我们后人去探索。MerChant、Piispanen 和 Leeandshaffer 是最早应用分析模型的学者，他们利用了切屑角模型来分析切屑在生成过程中的角度与刀具前角的关系。1960 年 Clough 在他的论文“平面分析的有限元法（The

13、Finite Element Method in Plane Stress Analysis） ”中最先引入了有限元（Finite Element）这一术语，一提出就引起了广泛的关注。1982 年 Usui 和 Shirakashi 第一次提出刀面角、切屑几何形状和流线等概念，得出了一个稳态的正交切削模型，应用此模型对切削过程中的应力应变和温度等参数进行了预测。2004 年合肥工业大学的谢峰和刘正士有限元法及弹塑性变形理论对二维金属切削的变形过程进行有限元分析，指出了金属由弹性变形到塑性变形时单元刚度矩阵变化的规律。2004 年邓文君等学者采用热力耦合、平面应变、连续带状切屑的切削模型模拟了高

14、强度耐磨铝青铜的正交切削加工过程。采用增量步移动刀具的方法，结合有限元分析软件 MARC 的网格重划分功能，模拟了刀具从初始切入到切削温度达到稳态的切削加工过程，获得了不同切削深度和切削速度下的切屑形态、温度、应力、应变和应变速率的分布。2005 年胡韦华、王秋成、胡晓冬和刘云峰针对有限元方法在切削加工过程中得到越来越广泛的应用，研究了切削加工过程数值模拟的研究进展情况，并对切削加工过程数值模拟的发展方向进行了展望。2006 年华南理工大学的何振威、全燕鸣和乐有树用 DEFORM-2D 软件建立了典型的正交切削模型，研究了高速切削中切削热在切屑、工件和刀具部分的量化分配规律。但是所得的结果只是

15、停留在模拟实验阶段，没有与实际试验做过比较。2009 年李泽文、罗洪波、端正强和肖华军对三角形刀片的切削过程进行了有限元分析，获得了不同切削用量对切削力、切削温度、刀片应力的影响，以选择合理的切削用量来延长刀具寿命。通过以上的了解，对于切屑形成过程的有限元模拟方面，国内的 2-D模拟发展相对 3-D 模拟较快，甚至能够对塑性较差的工件材料的切削进行2-D 模拟；国内的仿真主要停留在 2-D 领域，不仅视觉效果不够理想，而且模拟的能力仅局限于正交切削范围，其它大部分的切削情况都不能模拟，如包含斜刃切削的车削、刨削、铣削、钻削等切削加工情况。1.31.3 本课题的研究意义本课题的研究意义 随着全世

16、界越来越多的研究学者和工程师的努力，有限元法分析已经变得成熟和可靠，他们通过计算机所进行的大量的模拟科学研究，获得了大量宝贵的资料和成果。由于传统的研究方法难以定量分析切削机理，一旦面对高速、超精密切削加工等工艺，实验的方法便很难获得所需的相关参数，反而虚拟制造技术能够缩短开发周期、降低成本、提高产品质量，从而提高产品的市场竞争力。对切削过程进行虚拟仿真，研究金属切削变形等物理现象的影响因素，可以帮助合理选择参数工艺中的切削速度，背吃刀量及进给量;对刀具几何结构(前角，后角等)进行优化设计，进而可以采取措施减小切削力，提高金属切除效率并改善加工表面质量，优化加工工艺等。所以，采用计算机模拟技术

17、研究切削过程中的金属变形及其温度分布等是目前最有发展前景的研究方法之一。通过阅读文献了解相关研究背景和国内外研究现状，我对于应用有限元法模拟切削过程有了更清楚的认识，本课题能够全面地考量我在本科阶段的分析动手能力、对于新理论的应用学习能力以及对于问题的探索研究能力，能够完成本课题将是对我整个本科阶段的最佳总结。1.41.4 本课题研究的主要工作本课题研究的主要工作有限元法（finite element method ）是一种高效能，常用的计算方法。有限元法在早期是以变分原理为基础发展起来的，所以它广泛应用于以拉普拉斯方程和泊松方程所描述的各类物理场中。自从 1969 年以来，某些学者在流体力学

18、中应用加权余数法中的咖辽金法或最小二乘法同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的而各类物理场中，而不再要求这类物理场和泛函的极值有联系。边界元法 是一种及有限元法之后发展起来的一种新数字方法与有限元法在连续区域内划分单元的基本思想不同边界元法只是在定义域的边界上划分单元，用满足控制方程的函数去逼近边界条件。与有限元相比具有单元个数少，数据准备简单等优点。采用大型有限元分析软件 deform-3d 对普通麻花钻的钻削力进行仿真研究。Deform-3d 是 SFTC 公司开发的基于有限元分析的工艺仿真软件，针对复杂的金属成型过程，能够分析各种成型、热处理工艺，对加工过程中因

19、工件材料、刀具材料的剪切变形、切削温度内应力等因素进行分析，是正确选择刀具材料、刀具角度和切削用量以及进行材料加工性分析的依据。 研究方法 利用 cad 软件建立参数化麻花钻三维模型建立麻花钻的数学模型，对麻花钻的钻削过程受力进行理论分析在 ansys 下对刀具进行强度分析，获得 应变和应力分布图对麻花钻进行模态分析，得到各个模态下的振型图在 deform-3d 软件下进行钻削过程的模拟仿真 2 麻花钻切削理论基础2.1 麻花钻应用到的理论基础麻花钻应用到的理论基础2.1.12.1.1 麻花钻的组成麻花钻的组成麻花钻是孔加工过程中应用最广泛的一种工具，即可在实心材料上钻孔，也可将已有孔进行扩大

20、。可加工孔径的范围为 0.1-80mm。麻花钻的结构刀体：刀体是麻花钻的主要部分，它又分为切削部分和导向部分。切削部分担负主要切削工作，导向部分在工作时起导向的作用。随着麻花钻的刃磨变短，导向部分也是切削部分的准备部分。颈部：是刀体和刀柄的连接部分，又是钻头打标记的地方。刀柄：用于装夹钻头和传递动力。直径小的钻头为直柄。为了减小麻花钻与孔壁的摩擦，导向部分做有两条窄的刃带，且外径磨有倒锥量，即外径从切削部分向刀柄逐渐减少。麻花钻的两个刀齿靠钻心连接，因而两个主切削刃不通过钻心，而相互距离一个钻心直径。为了增大钻头的强度，把钻心做成正椎体，钻心从切削部分向刀柄逐渐增大，其增大量每 100mm 长

21、度为 1.4-1.8mm。如图。2.1所示。图 2.1 麻花钻的结构2.1.22.1.2 麻花钻的角度麻花钻的角度钻头实际上相当于正反安装的两把内孔车刀的组合刀具，只是这两把内孔车刀的主切削刃高于工件中心（因为有钻心而形成横刃的缘故，钻心半径为 ）端面刃倾角 为方便起见，钻头的刃倾角通常在端平面内表示。钻头主切削刃上某点的端面刃倾角是主切削刃在端平面的投影与该点基面之间的夹角。如图2.2 所示，其值总是负的。且主切削刃上各点的端面刃倾角是变化的，愈靠近钻头中心端面刃倾角的绝对值愈大。如图 2.2 所示主偏角 麻花钻主切削刃上某点的主偏角是该点基面上主切削刃的投影与钻头进给方向之间的夹角。由于主

22、切削刃上各点的基面不同，各点的主偏角也随之改变。主切削刃上各点的主偏角是变化的，外缘处大，钻心处小。图 2.2 麻花钻的主要角度前角 麻花钻的前角 是正交平面内前刀面与基面间的夹角。由于主切削刃上各点的基面不同，所以主切削刃上各点的前角也是变化的，如图 2.2所示。前角的值从外缘到钻心附近大约由+30减小到-30，其切削条件很差。后角 切削刃上任一点的后角 ，是该点的切削平面与后刀面之间的夹角。钻头后角不在主剖面内度量，而是在假定工作平面（进给剖面）内度量（见图 736a） 。在钻削过程中，实际起作用的是这个后角，同时测量也方便。钻头的后角是刃磨得到的，刃磨时要注意使其外缘处磨得小些（约 81

23、0） ，靠近钻心处要磨得大些（约 2030） 。这样刃磨的原因，是可以使后角与主切削刃前角的变化相适应，使各点的楔角大致相等，从而达到其锋利程度、强度、耐用度相对平衡；其次能弥补由于钻头的轴向进给运动而使刀刃上各点实际工作后角减少一个该点的合成速度角（见图 2.2 中 f-f 剖面）所产生的影响；此外还能改变横刃处的切削条件。图2.3麻花钻的前刀面和后刀面横刃斜角是在钻头的端面投影中，横刃与主切削刃之间的夹角。它是刃磨钻头时自然形成的，锋角一定时，后角刃磨正确的标准麻花钻横刃斜角 为 4755，而后角愈大则 愈小，横刃的长度会增加。2.1.32.1.3 基面和切削平面基面和切削平面在分析麻花钻

24、的几何角度时，首先必须弄清楚钻头的基面和切削平面。基面：切削刃上任一点的基面，是通过该点，且垂直于该点切削速度方向的平面，如图 2.4a 所示。在钻削时，如果忽略进给运动，钻头就只有圆周运动，主切削刃上每一点都绕钻头轴线做圆周运动，它的速度方向就是该点所在圆的切线方向，如图 2.4b 中 A 点的切削速度 垂直于 A点的半径方向，B 点的切削速度 垂直于 B 点的半径方向。不难看出，切削刃上任一点的基面就是通过该点并包含钻头轴线的平面。由于切削刃上各点的切削速度方向不同，所以切削刃上各点的基面也就不同。切削平面：切削刃上任一点的切削平面是包含该点切削速度方向，而又切于该点加工表面的平面（图 2

25、.4a 所示为钻头外缘刀尖 A 点的基面和切削平面） 。切削刃上各点的切削平面与基面在空间相互垂直，并且其位置是变化的。图 2.4 麻花钻的切削平面2.22.2 钻削力的计算钻削力的计算在钻头的每个切削刃上都作用着三个力，如图 2.5 所示。图 2.5 麻花钻受力示意图包括横刃和主刃上的轴向力径向力，。钻削力和扭矩主切削刃横刃棱带轴向力40573扭矩801010表 2.1 麻花钻受力分布轴向力F= (N)切削功率=2 （W）表 2.2 高速钢麻花钻切削力及功率计算公式计算公式 单位 d（mm） ；f（mm/r）;n(r/s),r=0.75，d=20mm，n=320r/min,f=0.25mm/

26、r工件材料钢333.541.90.8833.8510.7表 2.3 各主要参数代入参数计算得 F=6337N，按比例施加在横刃和主切削刃上。本次设计主要分析钻头受轴向力的影响，所以磨檫力等不予考虑。3 麻花钻的建模3.13.1 Pro/ENGINEERPro/ENGINEER3.1.13.1.1 Pro/ENGINEERPro/ENGINEER 的概述的概述1985 年，PTC 公司成立于美国波士顿，开始参数化建模 软件的研究。1988 年，V1.0 的 Pro/ENGINEER 诞生了。进过 10 余年的发展，Pro/ENGINEER 已经成为三维建模软件的领头羊。目前已经发布了Pro/EN

27、GINEERWildFire6.0。PTC 的 系列软件包括了在工业设计和机械设计等方面的多项功能，还包括对大型装配体的管理、功能仿真、制造、产品数据管理等功能。Pro/ENGINEER 还提供了全面集成紧密的产品开发环境。是一款有设计到生产的机械自动化软件，是新一代的产品造型系统，是一个参数化。基于特征的实体造型系统，并且具有单一数据库功能的综合性MCAD 软件。3.23.2 Pro/ENGINEERPro/ENGINEER 中麻花钻的三维建模中麻花钻的三维建模3.2.13.2.1 创建毛坯创建毛坯单击文件-新建命令，弹出新建对话框。在其中选择零件和实体单选按钮，在名称文本框中输入 zuan

28、tou，单击确定按钮，进入特征创建环境。如图 3.1图 3.1 创建环境选择创建拉伸特征。在工具栏单击拉伸按钮，弹出拉伸操控面板，在操控面板内单击放置按钮，然后但单击定义，弹出草绘对话框。定义草绘平面，如图所示，单击草绘按钮，进入草绘环境。如图 3.2图 3.2 草绘平面选取单击圆心和圆点按钮，绘制直径为 20mm 的圆，作为拉伸截面，如图 3.3所示。图 3.3 拉伸截面单击确定按钮，在操控面板内输入深度值为 238，单击确定按钮完成。如图 3.4图 3.4 拉伸圆柱单击倒圆角按钮，选取靠近top面的边界线，在操控面板内定义如图图 3.5TOP 面的选取图 3.6 刀柄倒角3.2.2 创建螺

29、旋槽创建螺旋线，单击编辑-包络命令，单击操控面板内的定义按钮，选择草绘平面。如图 3.7图 3.7 草绘平面的选取设置包络的圆点，绘制图形如图 3.8 所示图 3.8 包络原点的放置单击确定按钮，绘制的包络曲线如图 3.9 所示。图 3.9 包络曲线单击可变剖面扫描，设置操控面板如图 3.10 所示图 3.10 螺旋槽特征选取绘制扫描截面如图 3.11 所示图 3.11 扫描截面单击确定按钮，去除材料，效果如图 3.12 所示。图 3.12 单一螺旋槽的创建选择可变扫描特征，单击编辑-阵列，选择中心轴为旋转轴，阵列个数为 2 ，角度为 180，单击确定，结果如图 3.13 所示。图 3.13

30、双螺旋槽的创建3.2.3 创建横刃单击曲线按钮-经过点-确定-选择点-确定，得到横刃曲线如图 3.14所示。图 3.14 横刃曲线段选取单击平面按钮，已 top 面为基准，穿件 dtm1 基准面，如图 3.15 所示，便宜距离为 5.8.图 3.15dTM1 面的创建单击草绘按钮，已 DTM1 面为草绘平面，选择呢草绘命令，参照选择两个端点，绘制圆弧。如图 3.16 所示。图 3.16 圆弧的绘制.单击边界混合按钮，选择一圆弧和横刃直线，单击确定，如图 3.17所示。图 3.17 单边混合选择编辑命令，实体化，如图 3.18 所示。图 3.18 单边实体化对另外一边用同样的方法进行实体化，如图

31、 3.19 所示。图 3.19 横刃的创建对横刃进行修改，螺旋槽收尾，完成结果如图 3.20 所示。图 3.20 麻花钻4 有限元数值分析理论与 ansys 软件介绍4.14.1 有限元方法分析过程概述有限元方法分析过程概述有限元法的基本思想是：把要分析的连续体离散化，即将连续体变换成为有限个单元所组成的组合体，这些单元体之间只是通过结点来连接和制约。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用

32、力学理论中的变分原理或其他方法，建立结点力与位移之问的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。4.1.14.1.1 连续体的离散化连续体的离散化 首先应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。 4.1.24.1.2 单

33、元分析单元分析 建立各个单元的结点位移和结点力之问的关系式。 4.1.4.1.整体分析整体分析 3 程组，揭示结点外荷载与结点位移的关系，从而用来求解结点位移。 4.24.2 ANSYSANSYS 软件介绍软件介绍 ANSYS 软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国 ANSYS 开发，它能与多数 CAD 软件接口，实现数据的共享和交换，如 Pro/Engineer, NASTRAN, Alogor, I DEAS, AutoCAD 等， 是现代产品设计中的高级 CAD 工具之一。 4.2.14.2.1 软件功能简介软件

34、功能简介 软件主要包括三个部分：前处理模块，分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便地构造有限元模型；分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析） 、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。软件提供了 100 种以上的单元类型，用来模拟工程中的各种结构

35、和材料。 4.2.24.2.2 前处理模块前处理模块 PREP7PREP7 双击实用菜单中的“ Preprocessor” ，进入 ANSYS 的前处理模块。这个模块主要有两部分内容：实体建模和网格划分。 实体建模 ANSYS 程序提供了两种实体建模方法：自顶向下与自底向上。自顶向下进行实体建模时，用户定义一个模型的最高级图元，如球、棱柱，称为基元，程序则自动定义相关的面、线及关键点。用户利用这些高级图元直接构造几何模型，如二维的圆和矩形以及三维的块、球、锥和柱。无论使用自顶向下还是自底向上方法建模，用户均能使用布尔运算来组合数据集，从而“雕塑出”一个实体模型。 ANSYS 程序提供了完整的布

36、尔运算，诸如相加、相减、相交、分割、粘结和重叠。在创建复杂实体模型时，对线、面、体、基元的布尔操作能减少相当可观的建模工作量。 ANSYS 程序还提供了拖拉、延伸、旋转、移动、延伸和拷贝实体模型图元的功能。 自底向上进行实体建模时，用户从最低级的图元向上构造模型，即：用户首先定义关键点，然后依次是相关的线、面、体。 网格划分 ANSYS 程序提供了使用便捷、高质量的对 CAD 模型进行网格划分的功能。包括四种网格划分方法：延伸划分、映像划分、自由划分和自适应划分。延伸网格划分可将一个二维网格延伸成一个三维网格。映像网格划分允许用户将几何模型分解成简单的几部分，然后选择合适的单元属性和网格控制，

37、生成映像网格。 ANSYS 程序自由网格划分器功能是十分强大的，可对复杂模型直接划分，避免了用户对各个部分分别划分然后进行组装时各部分网格不匹配带来的麻烦。自适应网格划分是在生成了具有边界条件的实体模型以后，用户指示程序自动地生成有限元网格，分析、估计网格的离散误差，然后重新定义网格大小，再次分析计算、估计网格的离散误差，直至误差低于用户定义的值或达到用户定义的求解次数。 4.2.34.2.3 求解模块求解模块 SOLUTIONSOLUTION 前处理阶段完成建模以后，用户可以在求解阶段获得分析结果。 点击快捷工具区的 SAVE_DB 将前处理模块生成的模型存盘，退出 Preprocessor

38、 ，点击实用菜单项中的 Solution ，进入分析求解模块。在该阶段用户可以定义分析类型、分析选项、载荷数据和载荷步项，然后开始有限元求解。 ANSYS 软件提供的分析类型如下： 结构静力分析 用来求解外载荷引起的位移、应力和力。静力分析很适合求解惯性和阻尼对结构的影响并不显著的问题。 ANSYS 程序中的静力分析不仅可以进行线性分析，而且也可以进行非线性分析，如塑性、蠕变、膨胀、大变形、大应变及接触分析。 结构动力学分析 结构动力学分析用来求解随时间变化的载荷对结构或部件的影响。与静力分析不同，动力分析要考虑随时间变化的力载荷以及它对阻尼和惯性的影响。 ANSYS 可进行的结构动力学分析类

39、型包括：瞬态动力学分析、模态分析、谐波响应分析及随机振动响应分析。 结构非线性分析 结构非线性导致结构或部件的响应随外载荷不成比例变化。 ANSYS 程序可求解静态和瞬态非线性问题，包括材料非线性、几何非线性和单元非线性三种。 动力学分析 ANSYS 程序可以分析大型三维柔体运动。当运动的积累影响起主要作用时，可使用这些功能分析复杂结构在空间中的运动特性，并确定结构中由此产生的应力、应变和变形。 热分析 程序可处理热传递的三种基本类型：传导、对流和辐射。热传递的三种类型均可进行稳态和瞬态、线性和非线性分析。热分析还具有可以模拟材料固化和熔解过程的相变分析能力以及模拟热与结构应力之间的热结构耦合

40、分析能力。 电磁场分析 主要用于电磁场问题的分析，如电感、电容、磁通量密度、涡流、电场分布、磁力线分布、力、运动效应、电路和能量损失等。还可用于螺线管、调节器、发电机、变换器、磁体、加速器、电解槽及无损检测装置等的设计和分析领域。 流体动力学分析 ANSYS 流体单元能进行流体动力学分析，分析类型可以为瞬态或稳态。分析结果可以是每个节点的压力和通过每个单元的流率。并且可以利用后处理功能产生压力、流率和温度分布的图形显示。另外，还可以使用三维表面效应单元和热流管单元模拟结构的流体绕流并包括对流换热效应。 压电分析 用于分析二维或三维结构对 AC （交流） 、 DC （直流）或任意随时间变化的电流

41、或机械载荷的响应。这种分析类型可用于换热器、振荡器、谐振器、麦克风等部件及其它电子设备的结构动态性能分析。可进行四种类型的分析：静态分析、模态分析、谐波响应分析、瞬态响应分析. 4.34.3 静应力分析静应力分析 4.3.14.3.1 麻花钻的静应力分析步骤麻花钻的静应力分析步骤 指定分析标题并设置分析范畴 1. filestart New, ” Do Not File ” OK. 2. fileChang Job NamemahuazuanyesOK. 3.fileChange Title ” mahuazuan” OK 4.filechange Directory “ AQE 4 ” 5.

42、选取菜单途径 Main MenuPreference 6.单击 Structure 选项使之为 ON，单击 OK 定义单元类型 1 选取菜单途径 PreprocessorElement TypeAdd/Edit/Delete。 2.Element Types 对话框将出现。 3.单击 AddLibrary of Element Types 对话框将出现。 4.在左边的滚动框中单击“Solid” 。 5.在右边的滚动框中单击“10node 92” 。 6.单击 ok。 7.单击 Element Types 对话框中的 Close 按钮。 指定材料性能 1.选取菜单途径 PreprocessorM

43、aterial Props- Material Models-Define Material Material Models Behavior。 2.单击以指定材料号为 1。第二个对话框将出现 structural -linear-Elastic-Isotropic。 4.弹“Linear Isotropic Properties for material Number 1”对话框。 3.输入 EX 为 2.6e11。 4.输入 PRXY 为 0.3。 5.单击 OK。 导入麻花钻的模型 单击 file-import-iges。划分网格 如图 4.1图 4.1 网格属性定义 1. Prepro

44、cessor-MeshingMeshTool,弹出一对话框，下图设置参数： 点击“Mesh” 又弹出一选择对话框。点击 ” Pick All ” 2. 单元划分完成，如下图 4.2 所示： 图 4-2 网格示意图 定义约束和载荷 定义目标载荷:点击SolutionLodesApplyStructurasDisPlansmentON Areas；弹出选择对话框，选择麻花钻加载面。 加载约束，选择“ALL DOF”点击 OK ，如图 4.3 所示。图 4.3 全方向约束示意图3. 拾取横刃上的点。如图 4.4. 图 4.4 横刃上点的拾取4.给 Y，方向 ,施加 3612N 的力。 5.拾取主切削

45、刃上的点。如图 4.5.图 4.5 主切削刃上点的拾取6.在 Y 方向上施加 2535N 的力。如图 4.6 所示。 图 4.6 载荷施加示意图求解及其分析结果 1 SloutionSolveCurrent LS，弹出一对话框OK. 2. 求解完成后点击 Close 2. General PostrocPlot ResultsContours PlotNodal Solu.弹出如下对话框，Nodal Solution.StrssVonMisesStrss点击 OK；结果如下图 4.7 所示： 图 4.7 麻花钻等效应力分布图图 4.8 X 向应力分布图图4.9 Y向应力分布图图 4.10 Z

46、向应力分布图保存结果退出 ANSYS 对话框。 关掉 ANSYS 时，会弹出一对话框，选择“Save Everything”OK。退出分析。4.44.4 模态分析模态分析4.4.14.4.1 模态分析模态分析: :模态分析的定义和模态分析介绍模态分析的定义和模态分析介绍1.振动或以特定频率进行振动； 2.使工程师可以认识到结构对于不同类型的动力载荷是如何响 应的； 3.有助于在其它动力分析中估算求解控制参数（如时间步长） 4.由于结构的振动特性决定结构对于各种动力载荷的响应情 况，所以在准备进行其它动力分析之前首先要进行模态分析。 4.4.24.4.2 在在 ANSYSANSYS 中有以下几种

47、提取模态的方法：中有以下几种提取模态的方法： Block Lanczos 法，子空间法， PowerDynamics 法，缩减法，不对称法阻尼法。 4.4.34.4.3 模态分析中的四个主要步骤：模态分析中的四个主要步骤： 建模 选择分析类型和分析选项 施加边界条件并求解 评价结果评价结果 4.4.44.4.4 模态分析步骤模态分析步骤: : 指定分析标题并设置分析范畴 1. filestart New, ” Do Not FileNO ” OK。 2. fileChang Job NameMTyesOK。 3. fileChange Title ” MT ” OK。 4. filechang

48、e Directory “ MT ” 。 5. 选取菜单途径 Main MenuPreference； 单击 Structure 选项使之为 ON，单击 OK 定义单元类型 1.选取菜单途径 PreprocessorElement TypeAdd/Edit/Delete。 2.Element Types 对话框将出现。 3.单击 AddLibrary of Element Types 对话框将出现。 4.在左边的框中单击“Solid” ，在右边的滚动框中单击“10node 45 ” 单击 ok。 5.单击 Element Types 对话框中的 Close 按钮。 指定材料性能 1.选取菜单途

49、径 PreprocessorMaterial Props- Material Models-Define Material Material Models Behavior。 2.单击以指定材料号为 1。第二个对话框将出现 structural -linear-Elastic-Isotropic。 3.弹“Linear Isotropic Properties for material Number 1”对话框。 4.输入 EX 为 2.6e11。 5.输入 PRXY 为 0.3。点击“DENSITY”弹出一对话框，在 DENS 输入7800，单击 OK。 导入麻花钻的模型 1.fileimp

50、ortIGES ” ImportIGES File ” 对话框OK。指定IGES 模型存储路径 划分网格 1. Preprocessor-MeshingMeshTool,弹出一对话框，如下图 4.11设置参数： 图 4.11 网格划分设置 点击“Mesh” 又弹出一选择对话框。点击 ” Pick All ” 单元划分完成，如图 4.12 所示： 图 4.12 网格化示意图3. PlotArers.恢复单元的面积形式的模型，网格划分完毕。 进行模态分析设置 1. SolutionAnaiysis TypeNew Ansysis；选择静态分析ok 2. SolutionAnaiysis TypeA

51、nalysis Options，参数设置如下： 点击 OK，又弹出一对话框：第二行输入 10000，点击 OK SolutionDefinesLodesApplyStructurasDisPlansmentareas 点击 OK。如图 4.13. 图 4.13 约束定义设置 3.求解 查看结果 1. General PostrocRead RessultsBy Set Number，在“Scale factor”一栏输入想看的阶数，其他参数选默认。 2. General PostrocPolt ResultsConter PoltNodal Solu，选取“Stressvon Mises st

52、ress,查看应力图,应变图，振型图。图 4.14 一阶应力示意图图 4.15 二阶应力示意图图 4.16 三阶应力示意图图 4.17 四阶应力示意图图 4.18 一阶应变示意图图 4.19 二阶应变示意图图 4.20 三阶应变示意图图 4.21 四阶应变示意图图 4.22 一阶振型图图 4.23 二阶振型图图 4.24 三阶振型图图 4.25 四阶振型图图 4.26 列表显示5 钻削过程的模拟仿真与 deform-3d 软件介绍5.1deform-3d5.1deform-3d 软件介绍软件介绍5.1.15.1.1 deform-3ddeform-3d 软件简介软件简介DEFORM 是一套基于

53、有限元的工艺仿真系统，用于分析金属成形及其相关工业的各种成形工艺和热处理工艺。通过在计算机上模拟整个加工过程，帮助工程师和设计人员： 设计工具和产品工艺流程，减少昂贵的现场试验成本。 提高工模具设计效率，降低生产和材料成本。缩短新产品的研究开发周期。5.1.2deform-3d5.1.2deform-3d 的操作流程的操作流程（1）导入几何模型 在 DEFORM-3D 软件中，不能直接建立三维几何模型，必须通过其他 CAD/CAE 软件建模后导入导 DEFORM 系统中，目前， DEFORM-3D 的几何模型接口格式有 : STL:几乎所有的 CAD 软件都有这个接口。它由一系列的三角形拟合曲

54、面而成。 UNV:是由 SDRC 公司(现合并到 EDS 公司)开发的软件 IDEAS 制作的三维实体造型及有限元网格 文件格式，DEFOEM 接受其划分的网格。 PDA:MSC 公司的软件 Patran 的三维实体造型及有限元网格文件格式。 AMG:这种格式 DEFORM 存储己经导入的几何实体。 （2）网格划分 在 DEFORM-3D 中，如果用其自身带的网格剖分程序，只能划分四面体单元，这主要是为了考虑网格重划分时的方便和快捷。但是它也接收外部程序所生成的六面体 (砖块)网格。网格划分可以控制网格的密度，使网格的数量进一步减少，有不至于在变形剧烈的部位产生严重的网格畸变。 DEFORM-

55、3D 的前处理中网格划分有两种方式，一种是用户指定单元数量，系统默认划分方式，用户指定的网格单元数量只是网格划分的上限约数，实际划分的网格单元数量不会超过这个值。用户可以通过拖动滑块修改网格单元数，也可以直接输入指定数值，该数值和系统计算时间有着密切的关系，该数值越大，所需要的计算量越大，计算时间越长。 另一种手动设置网格使用的是 Detailed settings 下的 Absolute方式，该方式允许用户指定最小或最大的网格尺寸和最大与最小网格尺寸的比值。该值设置完成在网格单元数量中可以看到网格的大概数目，但无法在那里修改，只能通过修改最大或最小单元尺寸来修改网格数目。 （3）初始条件 有

56、些加工过程是在变温环境下进行的，比如热轧，在轧制过程中，工件，模具与周围环境介质之间存在热交换，工件内部因大变形生成的热量及其传导都对产品的成形质量产生主要的影响，对此问题，仿真分析应按照瞬态热一机祸合处理。 DEFORM 材料库可以提供各个温度下材料的特性。 （4）材料模型 在 DEFORM-3D 软件中，用户可以根据分析的需要，输入材料的弹性、塑性、热物理性能数据，如果需要分析热处理工艺，还可以输入材料的每一种相得相关数据以及硬化、扩散等数据。 为了更方便的使用户模拟塑性成形工艺，该软件提供了100 余种材料（包括碳钢、合金钢、铝合金、钛合金、铜合金等）的塑性性能数据，以及多种材料模型。在

57、材料库中，对每一种支持的材料提供了不同温度和应变率下材料流动应力应变 曲线和膨胀系数，弹性模量，泊松比，热导率等随温度的变化曲线。 （5）接触定义 接触菜单用于定义工件与所有用到的模具之间以及模具之间可能产生的接触关系。工件在变形过程中的温度，变形量是待求量，工件通常被定义成为可变形接触体。通常，最简单，计算效率最高的定义是用二维曲线(ZD 平面或是轴对称锻造 )或是三维空间曲面 (3D 锻造)描述模具参与接触部分的外表面轮廓，用刚性接触体描述。刚性接触体上只具有常温，起主动传递刚体位移或合力作用。如果需要关心模具的温度变化，可将模具上所关心的部分离散成单元 (二维平面单元或是三维轴对称实体单

58、元)，定义成为允许传热的刚性接触体，分析过程中，模具既有传递位移或合力作用，同时又有内部热量的传导和与外界的换热。实际锻造过程中，模具或多或少都存在变形，当要分析模具的温度和变形时，可将模具离散成为具有温度和位移自由度的有限单元，定义成为可变形的接触体，这会使计算的规模增加，但是分析结果更加合乎实际情况。还有一类刚性接触体为对称面，定义在工件上具有对称边界条件位置处，起施加对称边界条件的约束作用。定义的对称刚性平面可以满足法向的零位移约束和法向零热流约条件。 （6）网格自动重新划分 模拟分析过程中，单元附着在材料上，材料在流动过程中极易使相应的单元形状产生过度变形导致畸形。单元畸变后可能会中断

59、计算过程。因此，保证仿真过程中材料经过较大流动后分析仍然可以继续，获得的结果仍然具有足够的精度是非常重要的。 DEFORM 在网格畸变达到一定程度后会自动重新划分畸变的网格，生成新的高质量网格。对3D 分析，按增量加载频率或两组网格重划其间累积的最大应变增量来引导程序自动的网格重划。 （7）增加约束 DEFORM 可以在节点上增加各个自由度的约束。 （8）后处理 DEFORM 后处理菜单为用户提供了直观方便的评价成形过程，成形产品质量，工具损伤的必须信息以及图片，文本和表格形式提取和保存所需结果的各种工具。 DEFORM 支持在加工过程中以等值线，分布云图，数值符号，色标，等值面和切平面矢量等

60、方式显示各种场变量分布。也可按路径显示或历程显示分析结果。显示结果能够借助于色调，光照和渲染产生出具有逼真效果的图形。也可利用分析结果 制作动画和电影。用户利用这些提取各种体成形分析结果工具，足以获得设计产品加工工艺所关注的全部信息。这对设计人员充分了解设计工艺及其实施的可行性是大有裨益的。一旦模具设计和初始坯料形状尺寸不合理，从分析结果中可显示出材料流动受阻后可能出现的开裂或是重叠，从历程显示可以提取模具成形力随行程的变化曲线，是一个从设备加工能力，设备消耗角度来设计加工工艺的必须指标。在后处理界面中显示工件流动过程中应力，应变，应变率和温度的分布变化，帮助工艺 设计师评定工件的加工质量。其

61、中的局部加工硬化，应力集中，高应力梯度，工件模具的接触压力等结果，可以评定成形产品的质量好坏的控制因素。5.25.2 钻削过程仿真钻削过程仿真从该软件的道具库中选择半径为 5mm 的标准钻头，刀具材料选择高速钢。Cutting spread=20m/min,Feed rate=0.1mm/r,转速选择320r/min。工件材料选择中碳钢。由于是从刀具库中直接选择刀具，所以该仿真步骤省略。钻削过程的模拟，钻削过程中 X,Y,Z，方向的受力变化，温度变化的模拟等以动画形式附在说明书。下面选择个别步骤分析。图 5.1Step1X 方向受力图 5.2 Step50X 方向受力图 5.3 Step100

62、X 方向受力图 5.4 Step150X 方向受力图 5.5 Step200X 方向受力图 5.6 step1Y 方向受力图 5.7 step50Y 方向受力图 5.8 step100Y 方向受力图 5.9 step150Y 方向受力图 5.10 step200Y 方向受力图 5.11 step1Z 方向受力与应力变化图 5.12 step20Z 方向受力与应力变化图 5.13 step70Z 方向受力与应力变化图 5.14 step190Z 方向受力与应力变化图 5.15 step241Z 方向受力与应力变化钻削过程中的温度变化图 5.16 step1 温度变化图 5.17 step90 温

63、度变化图 5.18 step180 温度变化图 5.19 step204 温度变化图 5.20 step237 温度变化图 5.21 step271 温度变化致 谢在这将近 4 个月的毕业设计旅程中，我收获了很多知识、技能，也提高了自己的动手分析、解决问题的能力，为我以后迈入社会参加工作增添了更多的自信。在学院领导的热切关怀下，在老师的悉心指导下，在其他同学无私的帮助下，我完成了自己的本科毕业论文，给四年大学生涯画上了一个圆满的句点。 这其中有酸、有苦、有甜，在毕业前夕我感慨万千。 感谢学院老师们一直以来给予我的谆谆教诲，特别感谢候红玲老师在我毕业论文上的细致指导和辛勤付出，在我困惑的时候为我

64、耐心地解答；也感谢其他支持和帮助过我的人。 参考文献1 薛风先、胡任喜、康士延ANSYS12.0 机械与结构的有限元分析机械工业出版社 2010 年 3 月2 张应迁、张洪才ANSYS 有限元分析从入门到精通人民邮电出版社2010 年 7 月3 龚曙光ANSYS 基础应用及范例解析机械工业出版社 2004 年 1 月4 钟日铭Pro/ENGINEER 基础入门与范例解析清华大学出版社 2010年 3 月5 王健石机械加工常用刀具数据速查手册机械工业出版社 2009 年9 月6 杨荣富、董申金属切削原理机械工业出版社 1988 年 3 月7 华楚生机械制造技术基础重庆大学出版社 2000 年 7

65、 月8 谢大纲、赵清亮、袁哲俊、卢政君麻花钻刚度的有限元分析上海科技出版社 2003 年 9 月9 李珊、王磊基于 UG 的麻花钻三维建模研究昆明理工大学 2006年 9 月10 杨军、周利平基于 DEFORM-3D 的钻削力仿真研究西华大学 2003年 2 月11 韩荣第、杨昌琪、吴健麻花钻的数学建模及钻削过程有限元分析哈尔滨工业大学 2001 年 3 月12 Bing Li、Xiaoping Hu、Hao WangAnalysis and simulation for parallel drill point grinderInt J Adv Manuf technol 28Februar

66、y 200513 Eyup Bagci、Babar OzcelikEffects of different cooling conditions on twist drill temperatureInt J Adv Manuf technol November 200614 刘洪文材料力学高等教育出版社 2004 年 1 月15 哈尔滨理论力学教研室理论力学高等教育出版社 1961 年 7 月外文文献译文EffectsEffects ofof DifferentDifferent CoolingCooling ConditionsConditions onon TwistTwistDrillDrill TemperatureTemperatureEyup Bagci & Babur OzcelikAbstractAbstract: In the past, many researchers have studied wear developed on drilling tools mainly due to the high temperatures generated which