关 键 词：

YX15 225 900 彩板彩钢瓦压瓦 成型 设计

资源描述：

YX15-225-900彩板彩钢瓦压瓦成型工机设计,YX15,225,900,彩板彩钢瓦压瓦,成型,设计

内容简介：

攀枝花学院Panzhihua University本科毕业设计（论文）文献综述院 （系）： 机电工程学院 专 业： 机械设计制造及其自动化 班 级： 03级机制一班 学生姓名： 肖 洋 学 号: 200310621040 2007 年 4 月 25 日本科生毕业设计（论文）文献综述评价表毕业设计（论文）题目拔叉80-08的加工工艺及夹具设计综述名称拔叉80-08的加工工艺及夹具设计综述评阅教师姓名职称评 价 项 目优良合格不合格综述结构01文献综述结构完整、符合格式规范综述内容02能准确如实地阐述参考文献作者的论点和实验结果03文字通顺、精练、可读性和实用性强04反映题目所在知识领域内的新动态、新趋势、新水平、新原理、新技术等参考文献05中、英文参考文献的类型和数量符合规定要求，格式符合规范06围绕所选毕业设计（论文）题目搜集文献成绩综合评语： 评阅教师（签字）： 年 月 日文献综述：拔叉80-08的加工工艺及夹具设计1 前言随着科学技术的发展，各种新材料、新工艺和新技术不断涌现，机械制造工艺正向着高质量、高生产率和低成本方向发展。各种新工艺的出现，已突破传统的依靠机械能、切削力进行切削加工的范畴，可以加工各种难加工材料、复杂的型面和某些具有特殊要求的零件。数控机床的问世，提高了更新频率的小批量零件和形状复杂的零件加工的生产率及加工精度。特别是计算方法和计算机技术的迅速发展，极大地推动了机械加工工艺的进步，使工艺过程的自动化达到了一个新的阶段。“工欲善其事，必先利其器。”工具是人类文明进步的标志。自20世纪末期以来，现代制造技术与机械制造工艺自动化都有了长足的发展。但工具（含夹具、刀具、量具与辅具等）在不断的革新中，其功能仍然十分显著。机床夹具对零件加工的质量、生产率和产品成本都有着直接的影响。因此，无论在传统制造还是现代制造系统中，夹具都是重要的工艺装备。2 夹具的发展史 夹具在其发展的200多年历史中，大致经历了三个阶段：第一阶段，夹具在工件加工、制造的各工序中作为基本的夹持装置，发挥着夹固工件的最基本功用。随着军工生产及内燃机，汽车工业的不断发展，夹具逐渐在规模生产中发挥出其高效率及稳定加工质量的优越性，各类定位、夹紧装置的结构也日趋完善，夹具逐步发展成为机床工件工艺装备工艺系统中相当重要的组成部分。这是夹具发展的第二阶段。这一阶段，夹具发展的主要特点是高效率。在现代化生产的今天，各类高效率，自动化夹具在高效，高精度及适应性方面，已有了相当大的提高。随着电子技术，数控技术的发展，现代夹具的自动化和高适应性，已经使夹具与机床逐渐融为一体，使得中，小批量生产的生产效率逐步趋近于专业化的大批量生产的水平。这是夹具发展的第三个阶段，这一阶段，夹具的主要特点是高精度，高适应性。可以预见，夹具在不一个阶段的主要发展趋势将是逐步提高智能化水平。3 机床夹具机床夹具是零件在机床上加工时，用以装夹工件（和引导刀具）的一种工艺装备。其作用是正确确定工件与刀具之间的相对位置，并将工件固定地夹紧。使用夹具可以有效地保证工件的加工质量，提高劳动生产率，扩大机床的工艺范围和减轻劳动强度。3.1 机床夹具的分类若根据夹具的使用特点来划分，则有通用夹具、专用夹具、组合夹具、通用可调整夹具和成组夹具等。若按使用夹具的机床类型来划分，则为车床夹具、铣床夹具、钻床夹具、镗床夹具、磨床夹具、齿轮加工夹具、组合机床夹具等。若按夹具夹紧动力源的不同来划分，则有手动夹具、气动夹具、液压夹具、气液夹具、电动夹具、电磁夹具和真空夹具等。3.2 夹具的组成通常夹具由三大部分组成：3.2.1 定位装置定位装置包括定位元件极其组合，其作用是确定工件在夹具中的位置。常用的定位元件有支承钉、定位销、V型块等。3.2.2 夹紧装置夹紧装置的作用是将工件压紧夹牢，保证工件在定位时所占据的位置，在加工过程中不因受外力（切削力、重力、惯性力）等作用而产生位移，同时可以减轻或防止振动，它通常由夹紧元件（夹爪、压板等）、传动机构（如杠杆、斜楔等）和动力装置（气缸、液压缸等）组成。3.2.3 夹具体夹具体是用于连接夹具上所有的元件和装置，使其成为一个整体的基础件，它还用来与机床的有关连接，以确定夹具在机床上的位置。3.3 机床夹具的作用3.3.1 提高劳动生产率依靠夹具所设置的专门定位元件和高效夹紧装置，可以快速而准确地完成工件在加工工位上的定位和夹紧，省去对工件的逐个找正调整的装夹过程，大大缩短了每一工件的装夹辅助工时，这对于大批量生产的工件，尤其是对外形轮廓较复杂，不易找正的工件，选用高效夹具来安装，其效用就更大。3.3.2 保证工件的加工精度，稳定整批工件的加工质量夹具的设计和应用注重于解决工件的可靠定位和稳定装夹，可使用同一批工件的安装结果高度统一，稳定的装夹使各工件间的加工条件差异性大为减小，所以，采用夹具可以在保证工件加工精度的基础上极大地稳定整批工件的加工质量。3.3.3 改善工人劳动条件采用夹具后，使工件的装卸方便而快捷，减轻了工人的劳动强度。对较重的箱壳及大型工件，设计夹具时往往同时要考虑工件的工位间滚动输送装置及本工序的初始定位。另外，为保证工人生产安全，夹具设计还要考虑必要的防护封闭装置，加之以现代化的气动、液动及自动化机械手等装备的使用，无疑使工人的劳动条件大为改善，劳动强度大大降低。目前，现代化的加工自动线已经达到无人管理程度，即输送、安装、加工、卸料全部自动化，使得操作工人可以从繁重的劳动中彻底解脱出来。3.3.4 降低对操作工人的技术等级要求夹具的应用使得工件的装夹，操作大为简化，使得一些生产技术不熟练的技工可能胜任原来只能由熟练技术工人才能完成的复杂工件的精确装夹工作。因此，自动、高效夹具的实际应用，可以相应地降低对操作工人的装夹技术要求。4 总结一项优秀的夹具结构设计，往往可以使得生产效率大幅度提高，并使产品的加工质量得到极大地稳定。尤其是那些外形轮廓结构较复杂的，不规则的拔叉类，杆类工件，几乎各道工序都离不开专门设计的高效率夹具。目前，中等生产规模的机械加工生产企业，其夹具的设计，制造工作量，占新产品工艺准备工作量的50%80%。生产设计阶段，对夹具的选择和设计工作的重视程度，丝毫也不压于对机床设备及各类工艺参数的慎重选择。夹具的设计，制造和生产过程中对夹具的正确使用，维护和调整，对产品生产的优劣起着举足轻重的作用。 参 考 文 献1 李庆余，张 佳机械制造装备设计M 机械工业出版社，200382 濮良贵，纪名刚机械设计M 高等教育出版，20013 陈作模，孙 恒机械原理M高等教育出版，20024 王先逵机械制造工艺学M机械工业出版社，20005 冯 道机械零件切削加工工艺与技术标准实用手册M 安徽文化音像出版社，20036 徐鸿本机床夹具设计手册M 辽宁科学技术出版社，2003107 胡建新机床夹具M 中国劳动社会保障出版社，20015郑州科技学院毕业设计（论文）任务书题目 YX15-225-900彩板瓦成型工艺及专机设计 专业 机制 学号 姓名 徐 浩 主要内容、基本要求、时间安排、主要参考资料等：主要内容1：成型花纹展开图一张A1，2：成型机器装配图2张A03：辊轮轴零件装配图1张A04：机架装配图1.5张 A15：设计计算说明书一本（不少于2万字）5：设计计算说明书不低于2万字，翻译不低于3000汉字,文献综述不低于3000字。基本要求：1：原材料为彩色或镀锌板卷，宽度1000mm，厚度0.5mm-1mm 2：波纹瓦成型高度为29mm,间距125mm，倾角45，顶宽24mm(参考零件图)3：成型最小长度不小于1000mm, 最大长度不限，4：生产速度不小于0.5m/s，5：成型后不能产生翘曲扭转等现象，6：成型后表面涂层不能有划伤。时间安排第七学期17周到19周 前期准备13 周 查阅相关文献，写出文献综述，开题报告，翻译一篇有关的外文资料，到工厂进行调研、实习；46 周 拟定总体方案，绘制花型展开图。79 周 进行设备的总体方案设计；画出设计草图。 绘制滚压成型机装配图和零件图；1013周 整理设计资料、参考文献，写出完整详细的设计说明书。14 周 准备答辩主要参考资料： 机械设计（手册），机械原理，压力加工手册，滚压成型完 成 期 限： 指导教师签名： 专业负责人签名： 填 表 日 期： 1、工作原理彩板瓦成型装备又称辊式冷弯成型设备，是指在正常工作温度条件下，一定宽度的钢带或彩卷经过一组水平排列的的轧辊，在摩擦力的带动下前行并逐步的变形，达到符合要求的截面形状和尺寸长度的一种冷弯设备。彩钢成型设备加工流程一般分为三部分:辊压轴部分,压型部分,剪切部分。彩钢板的冷弯过程在工艺上近似于弯曲工艺，同时又与轧制工艺有所区别，属于特种生产方法。由于钢材在冷加工变形过程中会提高钢材的屈服强度，一般可以提高10-15，有的可以提高20，因此用这种逐步变形方法制造出来的彩板瓦，质量轻，强度高的特点。彩板瓦设备成型过程如图所示，当彩卷或彩钢板依次通过连续的辊轮时，在摩擦力的作用下，随着辊轮的回转向前进，同时逐渐进行弯曲变形，从而获得所要求的截面形状。另外在压型时，由于挤压，金属会发生变形，弯曲部位的金属厚度变薄，减薄量一般很小。在压型前后会认为钢板的厚度不变。这种彩板瓦成型装备适合大批量的连续的工业化生产，其产品的质量经久耐用，加工成本较低，生产效率比较高。图成型过程成卷连续生产的特点是在冷弯成型过程中可以轧制较大的弯曲深度或弯曲角。因为成卷连续生产分为变形区和成型区，卷料在开始成型时才开始进入成型辊中，弯曲时产生的的缺陷也仅仅会发生在料头和料尾，随后可以剪除掉，并且在喂料时不会存在喂料困难问题。同时连续生产效率高，不限制卷料的长度，生产极为方便。但是连续生产一般较为复杂，需要对不同形状的产品配置不同的轧辊、剪切机。而且剪切机剪刃加工不易，费用高。同时剪切机还需要控制系统进行控制。单件生产的特点是即适用于自动化生产，也可以用于人工生产。特别是用于人工生产时，设备较为简单，操作方便，而且适用于多种型号，各种产品。但是单件生产效率低，不能连续进行操作。根据生产的彩板瓦轧制的深度深，变形角度大，且连续轧制效率高的特点。本文设计的彩板瓦设备成型工艺为成卷连续生产。其生产工艺过程如下图轧辊侧面示意图1-矫直辊 2-压型辊 3-成型辊 4-检测装置 5-机架特点1，成型过程是摩擦力带着卷料连续进入辊轮并不断的向前运行的，有利于提高生产效率。2，成型过程中卷料的长度是无限制的。3，在压型过程中辊轮只与卷料局部接触，有利于减少耗能。4，在辊压部分分为变形区和成型区，卷料在成型之前就与辊轮接触并发生变形，因而变形区较长。5，由于卷料的厚度薄，变形阻力小，对设备的要求较低而且结构简单，易于操作。6，在设备连续运行时，成型工艺状态是稳定的，耗能、各个机构和辊轮中产生的作用力及作用力方向是不变的。7，在成型过程中，彩卷的横截面的厚度和面积基本不变，而截面形状发生显著变化，其高度增大，宽度减少，长度几乎不变。4 彩板瓦设备用来将彩卷经过多道辊轮逐渐变形，整形达到生产要求的形状和尺寸的设备。一般包括压辊，支撑定位压辊位置的压辊机架，成型机架以及传动系统。一般机械传动方式很多。主要有三种：1，带传动。带传动是一种挠性传动。基本组成零件为带轮和传动带。带传动具有结构简单、传动平稳、价格低廉和缓冲吸振等特点。通过传动带和带轮的摩擦来传递运动和动力，属于摩擦型传动。2，链传动。链传动是一种挠性传动。基本零件组成为链轮和链条。通过链轮的轮齿与链条的啮合来传递运动和动力。3，齿轮传动。与摩擦型的带传动相比，链传动无弹性滑动和整体打滑现象，因而能够保持准确的平均传动比，传动效率较高，一般在95左右。而且链条油不用像带那样张的很紧，作用于轴上的径向压力很小。与齿轮相比，链传动的制造和安装精度要求较低，成本也低。同时，在远距离的传动时，链传动在结构上比齿轮轻松的多。通过比较，本设计选定链传动为本次设计的传动方式。1、工作原理彩板瓦成型装备又称辊式冷弯成型设备，是指在正常工作温度条件下，一定宽度的钢带或彩卷经过一组水平排列的的轧辊，在摩擦力的带动下前行并逐步的变形，达到符合要求的截面形状和尺寸长度的一种冷弯设备。彩钢成型设备加工流程一般分为三部分:辊压轴部分,压型部分,剪切部分。彩钢板的冷弯过程在工艺上近似于弯曲工艺，同时又与轧制工艺有所区别，属于特种生产方法。由于钢材在冷加工变形过程中会提高钢材的屈服强度，一般可以提高10-15，有的可以提高20，因此用这种逐步变形方法制造出来的彩板瓦，质量轻，强度高的特点。彩板瓦设备成型过程如图所示，当彩卷或彩钢板依次通过连续的辊轮时，在摩擦力的作用下，随着辊轮的回转向前进，同时逐渐进行弯曲变形，从而获得所要求的截面形状。另外在压型时，由于挤压，金属会发生变形，弯曲部位的金属厚度变薄，减薄量一般很小。在压型前后会认为钢板的厚度不变。这种彩板瓦成型装备适合大批量的连续的工业化生产，其产品的质量经久耐用，加工成本较低，生产效率比较高。图成型过程成卷连续生产的特点是在冷弯成型过程中可以轧制较大的弯曲深度或弯曲角。因为成卷连续生产分为变形区和成型区，卷料在开始成型时才开始进入成型辊中，弯曲时产生的的缺陷也仅仅会发生在料头和料尾，随后可以剪除掉，并且在喂料时不会存在喂料困难问题。同时连续生产效率高，不限制卷料的长度，生产极为方便。但是连续生产一般较为复杂，需要对不同形状的产品配置不同的轧辊、剪切机。而且剪切机剪刃加工不易，费用高。同时剪切机还需要控制系统进行控制。单件生产的特点是即适用于自动化生产，也可以用于人工生产。特别是用于人工生产时，设备较为简单，操作方便，而且适用于多种型号，各种产品。但是单件生产效率低，不能连续进行操作。根据生产的彩板瓦轧制的深度深，变形角度大，且连续轧制效率高的特点。本文设计的彩板瓦设备成型工艺为成卷连续生产。其生产工艺过程如下图轧辊侧面示意图1-矫直辊 2-压型辊 3-成型辊 4-检测装置 5-机架特点1，成型过程是摩擦力带着卷料连续进入辊轮并不断的向前运行的，有利于提高生产效率。2，成型过程中卷料的长度是无限制的。3，在压型过程中辊轮只与卷料局部接触，有利于减少耗能。4，在辊压部分分为变形区和成型区，卷料在成型之前就与辊轮接触并发生变形，因而变形区较长。5，由于卷料的厚度薄，变形阻力小，对设备的要求较低而且结构简单，易于操作。6，在设备连续运行时，成型工艺状态是稳定的，耗能、各个机构和辊轮中产生的作用力及作用力方向是不变的。7，在成型过程中，彩卷的横截面的厚度和面积基本不变，而截面形状发生显著变化，其高度增大，宽度减少，长度几乎不变。4 彩板瓦设备用来将彩卷经过多道辊轮逐渐变形，整形达到生产要求的形状和尺寸的设备。一般包括压辊，支撑定位压辊位置的压辊机架，成型机架以及传动系统。一般机械传动方式很多。主要有三种：1，带传动。带传动是一种挠性传动。基本组成零件为带轮和传动带。带传动具有结构简单、传动平稳、价格低廉和缓冲吸振等特点。通过传动带和带轮的摩擦来传递运动和动力，属于摩擦型传动。2，链传动。链传动是一种挠性传动。基本零件组成为链轮和链条。通过链轮的轮齿与链条的啮合来传递运动和动力。3，齿轮传动。与摩擦型的带传动相比，链传动无弹性滑动和整体打滑现象，因而能够保持准确的平均传动比，传动效率较高，一般在95左右。而且链条油不用像带那样张的很紧，作用于轴上的径向压力很小。与齿轮相比，链传动的制造和安装精度要求较低，成本也低。同时，在远距离的传动时，链传动在结构上比齿轮轻松的多。通过比较，本设计选定链传动为本次设计的传动方式。立体光照成型的注塑模具工艺的综合模拟摘要功能性零部件都需要设计验证测试，车间试验，客户评价，以及生产计划。在小批量生产零件的时候，通过消除多重步骤，建立了有快速成型形成的注塑模具，这种方法可以保证缩短时间和节约成本。这种潜在的一体化由快速成型形成注塑模具的方法已经被多次证明是可行的。无论是模具设计还是注塑成型的过程中，缺少的是对如何修改这个模具材料和快速成型制造过程的影响有最根本的认识。此外，数字模拟技术现在已经成为模具设计工程师和工艺工程师开注塑模具的有用的工具。但目前所有的做常规注塑模具的模拟包已经不再适合这种新型的注塑模具，这主要是因为模具材料的成本变化很大。在本文中，以完成特定的数字模拟注塑液塑造成快速成型模具的综合方法已经发明出来了，而且还建立了相应的模拟系统。通过实验结果表明，目前这个方法非常适合处理快速成型模具中的问题。关键词注塑成型,数字模拟,快速成型引言在注塑成型中，聚合物熔体在高温和高压下进入模具中。因此，模具的材料需要有足够的热性能和机械性能来经受高温和高压的塑造循环。许多研究的焦点都是直接有快速成型形成注塑模具的过程。在生产小批量零件的时候，通过消除多重步骤，直接由快速成型形成的注塑模具可以保证缩短时间和节约成本。这种潜在的有快速成型形成注塑模具的方法已经被证明成功了。快速成型模具在性能上是有别与传统的金属模具。主要差异是导热性能和弹性模量（刚性）。举例来说，在立体光照成型模具中的聚合物的导热率小于铝制的工具的千分之一。在用快速成型技术来制造铸模时，整个模具设计和注塑成型工艺参数都需要修改和优化，传统的方法是改变彻底的刀具材料不过，目前还没有对如何修改这个模具材料的方法有根本的了解在当前的模具中，仅仅改变一些材料的性能是不能得到一个合理的结果的。同样，使用传统方法的时候，实际生产的零件也会有出先次品。因此，研究出一个快速成型过程，材料和注塑模具之间的互动关系是非常火急的。这样就可以确定模具设计标准和快速模具的注塑的技术。此外，计算机模拟是一种预测模塑件的质量的有效的方法。目前，商用仿真软件包已经成为模具设计师和工艺工程师在注塑过程中例行性的工具。不幸的是，目前常规注塑成型的模拟程序已经不再适用于这个快速成型模具，因为它极大的需要不同的刀具材料。例如，利用现在的仿真软件在铝和立体光照模具之间做个实验比较一下，虽然铝模具模拟植的部分失真是合理的，但是结果是不可以接受的，因为误差超过了百分之五十。在注塑成型中，失真主要是由于塑料零件的收缩和翘曲，模具也是一样的。对于通常模具，失真的主要因素是塑料件的收缩和翘曲，这个在目前的模拟中能测试准确。但是对于快速成型模具，潜在的失真会更多，在当前的测试中，其中就会有些失真会被忽视。例如，用一个简单的三步骤模拟分析模具变形的时候，就会出现很多偏差。在本文中，基于以上分析，一个新的快速成型模具的仿真系统已经开发出来了。拟议制度着重于预测部分失真，主要是用与预测快速成型模具的缺陷。先进的仿真系统可以用于预测快速成型模具设计和工艺是否最合理。我们的仿真系统已经被我们的实验证明是没有错误的。虽然有很多材料可以用于快速成型技术，但是我们还是专注于利用立体光照模具的技术来制造聚合物模具立体光照成型的过程是利用激光能量一层一层建立零件的部分。使用立体光照则可以体现出双方在快速成型工业的商业优势，而且在以后也可以生产出准确的，高品质的零部件。直到最近，立体光照主要是用于建立物理模型，为了检查视觉效果，仅仅只利用了它的一点点功能。不过，新一代的立体光照的光改善了立体化，机械性能，热学性能，所以它可以更好的应用于实际的模具中。2 综合仿真的成型过程2.1 方法 为了在注塑成型过程中模拟立体光照模具的功能，反复的试验中得到了一个方法。不同的软件组已经开发出来了，而且也已经做到了这一点。主要的假设是，温度和负载边界条件造成立体光照模具的扭曲，仿真步骤如下：部分几何模型则作为一个实体模型，这将通过流量分析软件包被翻译到一个文件中。模拟光聚合物模具中熔融体填充的过程，然后输出温度和压力的资料。在前一步获得了热负荷和边界条件，然后对光模具进行结构分析，其中失真的计算是在该注塑过程中进行的。如果模具的扭曲收敛了，那么直接进行下一步否则，扭曲的型腔（改动扭曲后的型腔的尺寸）返回第二个步骤，以熔体形式模拟注入扭曲的模具中。然后注射成型零件的收缩和翘曲模拟就开始应用了，算出该成型零件最终的扭曲部分上述的模拟流动中，基本上是三个仿真模块。2.2充型模拟的熔体2.2.1数字建模 计算机仿真技术已经能成功的预测到在极其复杂的几何形状下的填充情况。然而，目前大多数字模拟是基于一种混合有限元和有限差的中性平面上的。模拟软件包的应用过程基于这一模型说明图。然而，不同与系统中模具设计中的表面实体模型，这里所谓的中性平面（如图所示，图）是一个假想的在中间型腔中有距离和方向的一个平面，这个平面可能会在应用的过程中带来很大的不便。举例来说，模具表面常用于目前的快速成型系统中（通常是格式），所以当用模拟软件包的时候，第二次建模是不可避免的。那是因为模型在快速成型系统和仿真系统中是不一样的。考虑到这些缺点，在模拟系统中，型腔的表面将以基准面来引入，而不是中性平面。根据以往的调查，流量和温度场的方程式可以写为：X,Y是中性平面坐标系中的两个平面，是高度坐标，是，方向上的速度，是整体的平均厚度，, ,CP (T), K(T)分别表示聚合物的粘性，密度，周期热，热导率。图 是中性平面的模拟程序是维表面模型，是中性平面模型，是网状的平面模型，是最后的模拟结果此外，在高度方向上的边界条件的误差可以表示为：正如图中的中表示，TW 是恒壁温度.结合方程和方程，表明了u, v, T, P在坐标上面应该是对称的，因此在上半个高度中的平均u, v应该和整个高度中的平均u, v是一样的。根据这个特点，我们可以把整个型腔在上下高度上分为两个部分，正如图中的第一部分和第二部分。同时，型腔（如图）表面产生的三角有限元将替代了中性平面（如图）。因此，在高度方向上的有限元误差仅仅限于型腔表面，正如图所示，高度上的误差将从到。这是中性平面上的单一性。此外，从图到图，坐标也随之改变了。为了配合上述调整，方程仍是用方程。然而，原来的边界条件高度方向则改写为：与此同时，为了保持在同一坐标（）上的两部分能够流动，那么更多的边界条件必须满足。下标I和II则分别代表第一部分和第二部分的参数Cm-I 和Cm-II 则表示在填充阶段中分开的两个表面上的自由移动的熔融线。应该指出的是，方程与和方程与不同，和在数字模拟过程中将变的更难，主要原因是以下几点：同一个断层的表面都已经都已经有着特殊的网格，这将导致同一层上的独特的格局因此，在比较两个熔接口的时候，应该计算出各自的u, v, T, P。因为两个部分都有各自的流道通向节点和节点（如图所示）在同一段中，有可能两个都充满，也有可能一个满，一个空这两个情况应该分开处理，应该平均流动，使后者也分配到流动。这意味着在前线熔合处出现一点点小的误差是可以允许的通过控制时间和选择更好的位置来控制前线熔合节点。每个流场的边界都扩张到熔线前线，所以核查方程是否准确是相当重要的。鉴于上述分析，在同一个节点处的物理参数应该加以比较和调整。所以在进行模拟之前，描述同一节点有限元的信息应该准备好，也就是说，匹配的原理应该先预备好。图 表明表面模型中的中性平面的高度方向上的边界条件2.2.2数字模拟 压力场在建模中，粘度 是由于熔提的剪切速率，温度和压力引起的性能剪切变稀后，这就代表一个跨越式的模式，例如：其中对应于幂律指数，的特点是在在牛顿和幂律渐近极限之间的剪应力过渡区。无论在温度还是压力指数上，0(T, P)都可以有合理的表示，详情如下：方程11和12构成了一个五个常数，可以代表粘度，而且通过粘度的剪切速率的计算可以得到：根据上述情况，通过方程14，我们可以推断出一下充气压力方程：其中S是由计算出来的。运用伽辽金方法，对压力的有限元方程推导为：其中l是所有要素的的导线，包括节点N，而且其中i和j代表此处的N节点的数目，的计算方法如下：其中代表三角有限元，而代表有限元中的压力。 温度场中，为了确定高度方向上的误差，应该在模具表面上分为一层一层的三角有限元的网格。左边的能量方程4可以表示为：其中代表每一层N节点上的温度。热传导的计算方法是：其中l是所有要素，包括节点N，而且i和j分别代表此处的N节点个数。对流项的计算方法是：当是粘性热时，计算方法是：把方程1720带入方程4，温度方程变为：2.3 模具结构分析 结构分析的目的是预测在填充过程中，模具由于热和机械压力而产生的变形。这个模型是基于一个三维热边界元法。边界元法是比较适合这个应用的，因为只有变形的模具表面才有这样的信息。此外，边界元法有一个优点，那就是在计算变形的模具的时候，它的计算是不会白费的。 模具在所受载荷超过弹性范围的时候会产生应力。因此，在决定模具变形的时候，模具材料是一个基准。模具的热性能和力学性能是各向同性的，而且温度也是独立的。 尽管这个过程是循环的，但是相同时间的温度和热流都是可以用于计算模具变形的通常情况下，在模具里面每个瞬间温度都局限于型腔的表面和喷嘴的顶端。在观察距离的时候，瞬间的衰减变化是很微笑的，小于毫米这说明在模具的喷嘴处的变形是很小的，因此，忽略这个影响也是合理的稳态温度场满足拉普拉斯方程2T = 0的边界条件。至于机械边界条件，型腔表面受到熔体的压力，模具的表面会连接到工作台上的，而其他的外部表面将会假设是自由的.热边界的推导方程是大家都知道的，这是由于：其中uk, pk和分别是位移，牵引力和温度。, 是代表材料的膨胀系数和泊松比。Ulk是在方向上基本的位移。在一个三维空间中，各向同性弹性区域中，由一个单元产生的负荷主要集中在xl方向上，它是以下面的形式产生的：其中lk是Kronecker三角函数，是该模具材料的剪切模量。Plk的基本收缩都是在模具表面的每个节点处测量的，可以表示为：整个将分散在模具的表面上，转变为方程：其中n是指在这个区域上的表面成分。把恰当的线性函数代入方程，得到的线性边界方程就是模具的方程这个方程适用于每个离散的模具表面，从而组合成线性方程组，其中是节点的总数。每个节点有八个相关数量，三个位移组成部分，三个牵引组成部分，还有温度和热流量。在稳态热模型中，每个节点处的温度和磁场是已知的，余下的个量中，三个必须是已知的。此外，在若干个节点处的位移值的方程必须消除刚体运动和刚体自转的奇异系统。由此产生的系统方程式是一个集合起来的综合矩阵，它可以为有限元方法求解。基于方程的注塑假设，下面将给出元件的应力和应变：该偏元件的应力和应变分别是：用类似的方法可以预测在回火玻璃中的残余应力了。以积分的形式在平面上分析粘性和弹性结构关系时，可以表示为以下公式：其中G1是材料的的剪切模量。扩张的应变的情况如下：其中是材料体积的弹性模量，和的定义是：如果(t) = 0，那么方程到方程的结果则为：同样的，利用方程到方程消除应变xx(z, t)，得到：利用拉普拉斯变化方程，辅助系数R()由下面的方程得出：利用上述方程，并简化在模具中的应力和应变的形式，那么注塑中残余的应力在冷却阶段中，由下面的方程获得：方程可以通过梯形正交被解决。由于材料的时间在快速的变化，所以需要一个准数控程序来检测。辅助模量是检测数控梯形的规则。关于翘曲分析，节点位移和曲率将以壳单元表达为：其中 k 单元刚度矩阵，Be是衍生算子矩阵，d是位移，re是负载单元，可以由下面的方程得出：使用完整的三维有限元分析法的好处就是可以准确知道翘曲的结果。但是，当零件的形状很复杂的时候，它也是相当麻烦的。在本文中，在壳体理论基础上介绍了一种二维有限元分析方法。这种方法被大量使用是因为大多数注塑模具的零件都有一些部分几何的厚度远远小于其他部分。因此，那些部分则可以被作为一个集会的单元来预测翘曲。每三个节点壳单元组合成一个恒应变三角单元和一个离散克希霍夫三角元，如图所示，因此翘曲可以分为平面伸展变形和板弯曲变形。并相应的以单元刚度矩阵来描述翘曲的拉伸刚度矩阵和弯曲刚度矩阵。图 a-c是壳单元在局部坐标系统里的变形分解a是平面伸展元素，b是平面弯曲元素，c是壳单元三 实验验证 对提出的模型进行了评定和发展，最后核查是非常重要的。从模型模拟中得到的扭曲数据将和文献中的立体光照模具数据比较。如图所示，有一个注塑尺寸36 36 6毫米和实验数据中是相同的。薄壁和加强筋的厚度都是1.5毫米，这个注塑材料是聚丙烯。注塑机的型号是ARGURYHydronica320-210-750，它的工艺参数是，熔解温度是度，模具温度是度，注塑压力是.帕，保压时间是秒，冷却时间是秒。立体光照模具材料使用杜邦SOMOSTM树脂，能抵御高达度的高温。如上所述，热传导是区分立体光照模具和传统模具的一个重要因素。模具中的热量转移会产生温度的不均匀分布，所以导致了成型零件的翘曲立体光照成型模具的周期是可以预测的。以高的热传导率金属为背面做的薄壳立体光照模具将会增加自身的热传导率。图 模型腔图 不同的热传导率下，在方向上的扭曲失真比较实验值，三步走和常规都是指最后的实验结果常规是指实验中最好的结果三步走步骤的模拟过程分别与传统的注塑成型相似图 在不同的热传导率下，在方向上的扭曲失真比较 图 在不同热传导率下，在方向上扭曲失真比较图 不同热传导率下各个捻度变量的比较对于这个部分，扭曲包括三个方向上的位移和捻度（两个最初的平行边的夹角的误差）如图到图，实验结果表明，这些数值也包括通过传统注塑模具模拟系统预测的扭曲值和报道中的三步骤。结论本文介绍了一个综合模拟的快速成型模具的方法，并且建立了相应的仿真系统。为了验证这个系统，实验还进行了快速焊接立体光照成型模具。很明显，立体光照模具也会出现传统的注塑模具模拟软件一样的故障假设由于注射中的温度和负载荷引起了扭曲那么用三步骤完成的话，结果也会出现比较多的误差。不过更先进的模型会使结果更接近与实验。立体光照模具改进了热传导率极大的增加了零件质量由于温度比压力（负载）对模具的影响更大，所以改进立体光照模具的热传导率可以更显著的提高零件质量。无论零件多么复杂，快速成型技术可以使人们造型更快，更便捷，更便宜在快速成型稳步发展的基础上，快速制造也将随之而来，并且需要更多的精确工具来确定工艺过程的参数现行的模拟工具不能满足研究者研究模具相对的变化。正如本文中所述，对于一个综合模型来说，要预测最后零件质量是相当重要的。在不久的将来，我们期待看到通过快速成型扩展到快速模具制造的模拟程序。参考文献1 Wang KK (1980) System approach to injection molding process.Polym-Plast Technol Eng 14(1):7593.2 Shelesh-Nezhad K, Siores E (1997) Intelligent system for plastic injectionmolding process design. J Mater Process Technol 63(13):458462.3 Aluru R, Keefe M, Advani S (2001) Simulation of injection moldinginto rapid-prototyped molds. Rapid Prototyping J 7(1):4251.4 Shen SF (1984) Simulation of polymeric flows in the injection moldingprocess. Int J Numer Methods Fluids 4(2):171184.5 Agassant JF, Alles H, Philipon S, Vincent M (1988) Experimental andtheoretical study of the injection molding of thermoplastic materials.Polym Eng Sci 28(7):460468.6 Chiang HH, Hieber CA, Wang KK (1991) A unified simulation of thefilling and post-filling stages in injection molding. Part I: formulation.Polym Eng Sci 31(2):116124.7 Zhou H, Li D (2001) A numerical simulation of the filling stagein injection molding based on a surface model. Adv Polym Technol20(2):125131.8 Himasekhar K, Lottey J, Wang KK (1992) CAE of mold cooling in injectionmolding using a three-dimensional numerical simulation. J EngInd Trans ASME 114(2):213221.9 Tang LQ, Pochiraju K, Chassapis C, Manoochehri S (1998) Computeraidedoptimization approach for the design of injection mold coolingsystems. J Mech Des, Trans ASME 120(2):165174.10 Rizzo FJ, Shippy DJ (1977) An advanced boundary integral equationmethod for three-dimensional thermoelasticity. Int J Numer MethodsEng 11:17531768.11 Hartmann F (1980) Computing the C-matrix in non-smooth boundarypoints. In: New developments in boundary element methods, CML Publications,Southampton, pp 367379.12 Chen X, Lama YC, Li DQ (2000) Analysis of thermal residual stress inplastic injection molding. J Mater Process Technol 101(1):275280.13 Lee EH, Rogers TG (1960) Solution of viscoelastic stress analysisproblems using measured creep or relaxation function. J Appl Mech30(1):127134.14 Li Y (1997) Studies in direct tooling using stereolithography. Dissertation,University of Delaware, Newark, DE.IntegratedIntegratedIntegratedIntegrated simulationsimulationsimulationsimulation ofofofof thethethethe injectioninjectioninjectioninjection moldingmoldingmoldingmolding processprocessprocessprocesswithwithwithwith stereolithographystereolithographystereolithographystereolithography moldsmoldsmoldsmoldsAbstractAbstractAbstractAbstractFunctional parts are needed for design verification testing,field trials,customer evaluation, and production plan ning. By eliminating multiple steps, thecreationofthe injec tion mold directly by a rapid prototyping (RP) process holds thebest promise of reducing the time and cost needed to mold low-volume quantities ofparts. The potential of this integra tion of injection molding with RP has beendemonstrated many times. Whatismissingisthe fundamental understanding of howthe modifications to the mold material and RP manufacturing process impact both themold design and the injection mold ing process. In addition, numerical simulationtechniques have now become helpful tools of mold designers and process engi neersfor traditional injection molding. Butallcurrent simulation packages for conventionalinjection molding are no longer ap plicable to this new typeofinjection molds,mainly because the propertyofthe mold material changes greatly.Inthis paper, anintegrated approach to accomplish a numerical simulation of in jection molding intorapid-prototyped moldsisestablished and a corresponding simulation systemisdeveloped. Comparisonswithexperimental results are employed for verification,which show that the present schemeiswellsuited to handle RP fabri catedstereolithography (SL) molds.KeywordsKeywordsKeywordsKeywordsInjection moldingNumerical simulationRapid prototyping1 1 1 1 IntroductionIntroductionIntroductionIntroductionIn injection molding, the polymer melt at high temperatureisinjected into themold under high pressure 1. Thus, the mold material needs to have thermal andmechanical properties capa bleofwithstanding the temperatures and pressures ofthe mold ing cycle. The focus of many studies has been to create theinjection mold directly by a rapid prototyping (RP) process. By eliminatingmultiple steps, this method of tooling holds the best promise of reducing the time andcost needed to create low-volume quantities of parts in a production material. ThepotentialofintegratinginjectionmoldingwithRPtechnologieshasbeendemonstrated many times. The properties of RP molds are very different from thoseof traditional metal molds. The key differ ences are the properties of thermalconductivity and elastic mod ulus (rigidity). For example, the polymers used inRP-fabricated stereolithography (SL) molds have a thermal conductivity thatislessthan one thousandth that of an aluminum tool. In using RP technologies to createmolds, the entire mold design and injection-molding process parameters need to bemodified and optimized from traditional methodologies due to the completelydifferent tool material. However, thereisstillnota fundamen tal understanding ofhow the modifications to the mold tooling method and material impact both the molddesign and the injec tion molding process parameters. One cannot obtain reasonableresultsbysimply changing a few material properties in current models. Also, usingtraditional approaches when making actual parts may be generating sub-optimalresults. So thereisa dire need to study the interaction between the rapid tooling (RT)pro cess and material and injection molding, so as to establish the mold designcriteria and techniques for an RT-oriented injection molding process.In addition, computer simulationisaneffective approach for predicting thequality of moldedparts. Commerciallyavailablesimulation packages of thetraditional injection molding process have now become routine toolsofthe molddesigner and pro cess engineer 2. Unfortunately, current simulation programs forconventional injection molding arenolonger applicable to RP molds, because of thedramatically dissimilar tool material. For instance, in using the existing simulationsoftware with alu minum and SL molds and comparing with experimental results,though the simulation values of part distortion are reasonable for the aluminum mold,results are unacceptable, with the error exceeding 50%. The distortion duringinjection moldingisdue to shrinkage and warpage of the plastic part, aswellas themold. For ordinarily molds, the main factoristhe shrinkage and warpage of theplastic part, whichismodeled accurately in cur rent simulations. But for RP molds,the distortion of the mold has potentially more influence, which have been neglectedin current models. For instance, 3 used a simple three-step simulation process toconsider the mold distortion, which had too much deviation.In this paper, based on the above analysis, a new simula tion system for RPmoldsisdeveloped. The proposed system focuses on predicting part distortion, whichisdominating defect in RP-molded parts. The developed simulationcanbe applied asan evaluation tool for RP mold design and process opti mization. Our simulationsystemisverifiedbyan experimental example.Although many materials are available for use in RP tech nologies, weconcentrateonusing stereolithography (SL), the original RP technology, to createpolymer molds. The SL pro cess uses photopolymer and laser energy to build a partlayerbylayer. Using SL takes advantage of both the commercial domi nanceofSLin the RP industry and the subsequent expertise base that has been developed forcreating accurate, high-quality parts.Untilrecently, SL was primarily used to createphysical models for visual inspection and form-fitstudieswithvery limitedfunc tional applications. However,thenewer generationstereolitho graphicphotopolymers have improved dimensional,mechanical and thermal propertiesmakingitpossible to use them for actual functional molds.2 2 2 2 IntegratedIntegratedIntegratedIntegrated simulationsimulationsimulationsimulation ofofofof thethethethe moldingmoldingmoldingmolding processprocessprocessprocess2.1 MethodologyIn order to simulate the use of an SL mold in the injection molding process, aniterative methodisproposed. Different soft ware modules have been developed andused to accomplish this task. The main assumptionisthat temperature and loadbound ary conditions cause significant distortions in the SL mold. The simulationsteps are as follows:1The part geometryismodeled as a solid model, whichistranslated to afilereadable by theflow analysis package.2Simulate the mold-fillingprocess of the melt into a pho topolymer mold,whichwilloutput the resulting temperature and pressure profiles.3Structural analysisisthen performed on the photopolymer mold modelusing the thermal and load boundary conditions obtained from the previous step,which calculates the distor tion that the mold undergo during the injection process.4Ifthe distortion of the mold converges, move to the next step. Otherwise,the distorted mold cavityisthen modeled (changes in the dimensions of the cavityafter distortion), and returns to the second step to simulate the melt injection into thedistorted mold.5The shrinkage and warpage simulation of the injection molded partisthenapplied, which calculates thefinaldistor tions of the molded part.In above simulationflow, there are three basic simulation mod ules.2.2Filling simulationof themelt2.2.1 Mathematical modelingIn order to simulate the use of an SL mold in the injection molding process, aniterative methodisproposed. Different software modules have been developed andused to accomplish this task. The main assumptionisthat temperature and loadboundary conditions cause significant distortionsinthe SL mold. The simulation stepsare as follows:1. The part geometryismodeled as a solid model, whichistranslated to a filereadable by the flow analysis package.2. Simulate the mold-filling process of the melt into a photopolymer mold, whichwilloutput the resulting temperature and pressure profiles.3. Structural analysisisthen performedonthe photopolymer mold model usingthe thermal and load boundary conditions obtained from the previous step, whichcalculates the distortion that the mold undergo during the injection process.4.Ifthe distortion of the mold converges, move to the next step. Otherwise, thedistorted mold cavityisthen modeled (changesinthe dimensions of the cavity afterdistortion), and returns to the second step to simulate the melt injection into thedistorted mold.5. The shrinkage and warpage simulationofthe injection molded partisthenapplied, which calculates the final distortionsofthe molded part.In above simulation flow, there are three basic simulation modules.2.2 Filling simulation ofthe melt2.2.1 Mathematical modelingComputer simulation techniques have had success in predictingfillingbehaviorin extremely complicated geometries. However, most of the current numericalimplementationisbasedona hybrid finite-element/finite-difference solution with themiddleplane model. The application processofsimulation packages basedonthismodelisillustrated in Fig. 2-1. However, unlike the surface/solidmodel inmold-design CAD systems, the so-called middle-plane (as shown in Fig. 2-1b)isanimaginary arbitrary planar geometry at the middle of the cavity in the gap-wisedirection, which should bring about great inconvenience in applications. For example,surface models are commonly used in current RP systems (generally STL file format),so secondary modelingisunavoidable when using simulation packages because themodels in the RP and simulation systems are different. Considering these defects, thesurface model of the cavityisintroduced as datum planes in the simulation, instead ofthe middle-plane.According to the previous investigations 46, fillinggoverning equations for theflow and temperature field can be written as:wherex, yare the planar coordinates in the middle-plane, andzisthe gap-wisecoordinate;u, v,ware the velocity componentsinthex, y, zdirections;u, vare theaverage whole-gap thicknesses; and, ,CP(T), K(T)represent viscosity, density,specific heat and thermal conductivity of polymer melt, respectively.Fig.2-1Fig.2-1Fig.2-1Fig.2-1 a a a a d. d. d. d. Schematic procedure of thesimulation with middle-plane model. a a a aThe3-D surfacemodelb b b bThemiddle-plane model c c c c Themeshed middle-plane modeld d d dThedisplay of thesimulation resultIn addition, boundary conditions in the gap-wise direction can be defined as:whereTWisthe constantwalltemperature (shown in Fig. 2a).Combining Eqs. 14 with Eqs. 56,itfollows that the distributions of theu, v, T,Patzcoordinates should be symmetrical, with the mirror axis beingz= 0, andconsequently theu, vaveraged in half-gap thicknessisequal to that averaged inwholegap thickness. Basedonthis characteristic, we can divide the whole cavity intotwo equal parts in the gap-wise direction, as described byPartIandPartIIin Fig. 2b.At the same time, triangular finite elements are generatedinthe surface(s) of thecavity(atz= 0 in Fig. 2b), insteadofthe middle-plane(atz= 0 in Fig. 2a).Accordingly, finite-difference increments in the gapwise direction are employed onlyin the inside of the surface(s)(wallto middle/center-line), which, in Fig. 2b, meansfromz= 0 toz=b. Thisissingle-sided instead of two-sided with respect to themiddle-plane (i.e. from the middle-line to two walls).Inaddition, the coordinatesystemischanged from Fig. 2a toFig.2b to alter the finite-element/finite-differencescheme, as shown in Fig. 2b. With the above adjustment, governing equations are stillEqs. 14. However, the original boundary conditionsinthe gapwise direction arerewritten as:Meanwhile, additional boundary conditions must be employed atz=bin orderto keep the flows at the juncture of the two parts at the same section coordinate 7:where subscripts I,IIrepresent the parametersofPartIandPartII, respectively,and Cm-I and Cm-II indicate the moving free melt-fronts of the surfaces of thedivided two parts in the filling stage.Itshould be noted that, unlike conditions Eqs. 7 and 8, ensuring conditions Eqs.9 and 10 are upheld in numerical implementations becomes more difficult due to thefollowing reasons:1. The surfaces at the same section have been meshed respectively, which leadsto a distinctive pattern of finite elements at the same section. Thus, an interpolationoperation should be employed foru, v, T, Pduring the comparison between the twoparts at the juncture.2. Because the two parts have respective flow fields with respect to the nodes atpoint A and point C (as shown in Fig. 2b) at the same section,itispossible to haveeither both filled or one filled (and one empty). These two cases should be handledseparately, averaging the operation for the former, whereas assigning operation for thelatter.3.Itfollows that a small difference between the melt-frontsispermissible. Thatallowance can be implementedbytime allowance control or preferable locationallowance control of the melt-front nodes.4. The boundaries of the flow field expandbyeach melt-front advancement, soitisnecessary to check the condition Eq. 10 after each change in the melt-front.5. In view of above-mentioned analysis, the physical parameters at the nodes ofthe same section should be compared and adjusted, so the information describingfinite elements of the same section should be prepared before simulation, that is, thematching operation among the elements should be preformed.Fig.Fig.Fig.Fig. 2a,b.2a,b.2a,b.2a,b. Illustrative of boundary conditionsinthe gap-wise direction a a a aof themiddle-planemodelb b b bof thesurfacemodel2.2.2 Numerical implementationPressure field.In modeling viscosity, whichisa functionofshear rate,temperature and pressureofmelt, the shear-thinning behavior can bewellrepresentedby a cross-type model such as:wherencorresponds to the power-law index, and*characterizes the shearstress level of the transition region between the Newtonian and power-law asymptoticlimits. In terms ofanArrhenius-type temperature sensitivity and exponential pressure dependence,0(T, P)can be represented with reasonable accuracy as follows:Equations 11 and 12 constitute a five-constant(n,* ,B,Tb,)representationfor viscosity. The shear rate for viscosity calculationisobtainedby:Based on the above, we can infer the following filling pressure equation from thegoverning Eqs. 14:whereSiscalculatedbyS=b0/(bz)2dz. Applying the Galerkin method, thepressure finite-element equationisdeduced as:wherel_ traversesallelements, including nodeN, and whereIandjrepresent thelocal node number in elementl_ corresponding to the node number N andN_ in thewhole, respectively. TheD(l_)ijiscalculated as follows:whereA(l_)represents triangular finite elements, andL(l_)iisthe pressure trialfunction in finite elements.Temperature field.To determine the temperature profile across the gap, eachtriangular finite element at the surfaceisfurther divided intoNZlayers for thefinite-difference grid.The leftitemofthe energy equation (Eq. 4)canbe expressed as:whereTN, j,trepresents the temperature of thejlayerofnodeNat timet. Theheat conductionitemiscalculatedby:whereltraversesallelements, including nodeN, andiandjrepresent the localnode number in elementlcorresponding to the node numberNandN_ in the whole,respectively.The heat convectionitemiscalculatedby:For viscous heat,itfollowsthat:Substituting Eqs. 1720 into the energy equation (Eq. 4), the temperatureequation becomes:2.3 Structural analysis ofthemoldThe purpose of structural analysisisto predict the deformation occurring in thephotopolymer mold due to the thermal and mechanical loads of the filling process.This modelisbased on a three-dimensional thermoelastic boundary element method(BEM). The BEMisideally suited for this application becauseonlythe deformationof the mold surfacesisof interest. Moreover, the BEMhasan advantage over othertechniques in that computing effortisnot wasted on calculating deformation withinthe mold.The stresses resulting from the process loads arewellwithin the elastic rangeofthe mold material. Therefore, the mold deformation modelisbasedona thermoelasticformulation. The thermal and mechanical properties of the mold are assumed to beisotropic and temperature independent.Although the processiscyclic, time-averaged values of temperature and heatflux are used for calculating the mold deformation. Typically, transient temperaturevariations within a mold have been restricted to regions local to the cavity surface andthe nozzletip8. The transients decay sharply with distance from the cavity surfaceand generally little variationisobserved beyond distances as small as 2.5 mm. Thissuggests that the contribution from the transients to the deformation at the mold blockinterfaceissmall, and thereforeitisreasonable to neglect the transient effects. Thesteadystatetemperaturefieldsatisfies Laplaces equation2T=0 andthetime-averaged boundary conditions. The boundary conditions on the mold surfacesare describedindetail by Tang et al. 9. As for the mechanical boundary conditions,the cavity surfaceissubjected to the melt pressure, the surfaces of the mold connectedto the worktable are fixed in space, and other external surfaces are assumed to bestress free.The derivation of the thermoelastic boundary integral formulationiswellknown10.Itisgivenby:whereuk,pkandTare the displacement, traction and temperature,representthe thermal expansion coefficient and Poissons ratio of the material, andr=|yx|.clk(x)isthe surface coefficient which dependsonthe local geometry atx, theorientation of the coordinate frame and Poissons ratio for the domain 11. Thefundamental displacementulkat a pointyin thexkdirection, in a three-dimensionalinfinite isotropic elastic domain, results from a unit load concentrated at a pointxacting in thexldirection andisof the form:wherelkisthe Kronecker delta function andisthe shear modulus of the moldmaterial.The fundamental tractionplk, measured at the pointyon a surface with unitnormaln n n n,is:Discretizing the surface of the mold into atotalofNelements transforms Eq. 22to:wherenrefers to thenthsurface elementonthe domain.Substituting the appropriate linear shape functions into Eq. 25, the linearboundary element formulation for the mold deformation modelisobtained. Theequationisapplied at each node on the discretized mold surface, thus giving a systemof 3Nlinear equations, whereNisthetotalnumber of nodes. Each node has eightassociated quantities: three components of displacement, three components of traction,a temperature and a heat flux. The steady state thermal model supplies temperatureand flux values as known quantities for each node, and of the remaining six quantities,three must be specified. Moreover, the displacement values specified at a certainnumber of nodes must eliminate the possibility of a rigid-body motion or rigid-bodyrotation to ensure a non-singular system of equations. The resulting system ofequationsisassembled into a integrated matrix, whichissolved withaniterativesolver.2.4 Shrinkage and warpage simulation ofthemoldedpartInternal stresses in injection-molded components are the principal cause ofshrinkage and warpage. These residual stresses are mainly frozen-in thermal stressesdue to inhomogeneous cooling, when surface layers stiffen sooner than the coreregion, as in free quenching. Based on the assumption of the linear thermo-elastic andlinearthermo-viscoelasticcompressiblebehaviorofthepolymericmaterials,shrinkage and warpage are obtained implicitly using displacement formulations, andthe governing equationscanbe solved numerically using a finite element method.With the basic assumptionsofinjection molding 12, the components of stressand strain are givenby:The deviatoric components of stress and strain, respectively, are given byUsing a similar approach developedbyLee and Rogers 13 for predicting theresidual stresses in the tempering of glass, an integral form of the viscoelasticconstitutive relationshipsisused, and the in-plane stresses can be related to the strainsby the following equation:WhereG1isthe relaxation shear modulus of the material. The dilatationalstresses can be related to the strain as follows:WhereKisthe relaxation bulk modulus of the material, and the definition ofandis:If(t)=0, applying Eq. 27 to Eq. 29 resultsin:Similarly, applying Eq. 31 to Eq. 28 and eliminating strainxx(z, t)resultsin:Employing a Laplace transform to Eq. 32, the auxiliary modulusR()isgivenby:Using the above constitutive equation (Eq. 33) and simplified forms of thestresses and strains in the mold, the formulation of the residual stress of the injectionmolded part during the cooling stageisobtainby:Equation 34 can be solved through the application of trapezoidal quadrature. Dueto the rapid initial change in the material time, a quasi-numerical procedureisemployed for evaluating the integral item. The auxiliary modulusisevaluatednumerically by the trapezoidal rule.For warpage analysis, nodal displacements and curvatures for shell elements areexpressed as:wherek isthe element stiffness matrix,Be isthe derivative operator matrix,disthe displacements, andreisthe element load vector which can be evaluatedby:The use of afullthree-dimensional FEM analysiscanachieve accurate warpageresults, however,itiscumbersome when the shape of the partisvery complicated. Inthis paper, a twodimensional FEM method, basedonshell theory, was used becausemost injection-molded parts have a sheet-like geometry in which the thicknessismuch smaller than the other dimensions of the part. Therefore, the part can beregarded asanassembly of flat elements to predict warpage. Each three-nodeshellelementisa combinationofa constant strain triangular element (CST) and a discreteKirchhoff triangular element (DKT), as shown in Fig. 3. Thus, the warpage can beseparatedintoplane-stretchingdeformationoftheCSTandplate-bendingdeformation of the DKT, and correspondingly, the element stiffness matrix todescribe warpagecanalso be divided into the stretching-stiffness matrix andbending-stiffness matrix.Fig.Fig.Fig.Fig. 3a3a3a3a c. c. c. c. Deformation decomposition of shell elementinthelocal coordinate system. a a a aIn-planestretching elementb b b bPlate-bending element c c c cShell element3 3 3 3 ExperimentalExperimentalExperimentalExperimental validationvalidationvalidationvalidationTo assess the usefulnessofthe proposed model and developed program,verificationisimportant. The distortions obtained from the simulation model arecompared to the ones from SL injection molding experiments whose dataispresentedin the literature 8. A common injection molded partwiththe dimensions of 36366mmisconsidered in the experiment, as shown in Fig. 4. The thickness dimensions ofthe thin walls and rib are both 1.5 mm; and polypropylene was used as the injectionmaterial. The injection machine was a production level ARGURY Hydronica320-210-750 with the following process parameters: a melt temperatureof250 C; anambient temperatureof30 C;aninjection pressure of 13.79 MPa; an injection timeof 3s;and a cooling time of 48 s. The SL material used, Dupont SOMOSTM 6110resin,hasthe ability to resist temperatures ofupto 300 C temperatures. Asmentioned above, thermal conductivity of the moldisa major factor that differentiatesbetweenanSL and a traditional mold. Poor heat transfer in the mold would produce anon-uniform temperaturedistribution,thus causingwarpagethat distortsthecompleted parts. For an SL mold, a longer cycle time would be expected. The methodof using a thin shell SL mold backed with a higher thermal conductivity metal(aluminum) was selected to increase thermal conductivityofthe SL mold.Fig.Fig.Fig.Fig. 4. 4. 4. 4. Experimental cavity modelFig.Fig.Fig.Fig. 5. 5. 5. 5. A comparison of the distortion variationinthe X direction for different thermalconductivity; where “Experimental”, “present”, “three-step”, and “conventional” mean the resultsof the experimental, the presented simulation, the three-step simulation process and theconventional injection molding simulation, respectively.Fig.Fig.Fig.Fig. 6. 6. 6. 6. Com