转向杆加工工艺分析钻削大端Φ14锁紧孔的立式钻床夹具设计【含三维PROE】【含5张CAD图纸+PDF图】

喜欢这套资料就充值下载吧。资源目录里展示的都可在线预览哦。下载后都有，请放心下载，文件全都包含在内，【有疑问咨询QQ:414951605 或 1304139763】 江西科技学院本科生毕业设计（论文）任务书（工科及部分理科专业使用）题 目： 转向杆加工工艺分析与夹具设计 学 院： 专 业： 班 级： 学 号： 学生姓名： 陈焕林 起讫日期： 指导教师： 孙淑梅 职称： 高级工程师 学院院长： 审核日期： 说 明1. 毕业论文任务书由指导教师填写，并经学院或专业学科组审定，下达到学生。2. 进度表由学生填写，每两周交指导教师签署审查意见，并作为毕业论文工作检查的主要依据。3. 学生根据指导教师下达的任务书独立完成开题报告，1个月内提交给指导教师批阅。4. 本任务书在毕业设计（论文）完成后，与论文一起交指导教师，作为论文评阅和毕业设计（论文）答辩的主要档案资料，是学士学位论文成册的主要内容之一。一、毕业设计（论文）的要求和内容（包括原始数据、技术要求、工作要求）1、任务及要求：根据给定的要求，分析研究转向杆的加工工艺，明确工件的结构、材料，产量等，遵照六点定位规则确定工件的定位方案，正确设计钻削大端14锁紧孔的立式钻床夹具，定位元件、导引元件、夹紧元件和夹具体。绘制夹具装配图和相关零件图。转向杆材料35#钢，中等批量生产。2、原始参数：转向杆的工作简图如下：二、毕业设计（论文）图纸内容及张数1、编制零件加工工艺过程卡 ；2、运用Pro/E和AutoCAD绘制夹具装配图及相关零件工作图； 3、编制设计计算说明书 。三、毕业设计（论文）实物内容及要求 无 四、毕业设计（论文）进度计划序号各阶段工作内容起讫日期备 注五、主要参考资料六、毕业设计（论文）进度表（本表每两周由学生填写一次，交指导教师签署审查意见）第一、二周时间: （ 月 日 至 月 日）学生主要工作：指导教师审核：年 月 日第三、四周时间: （ 月 日 至 月 日）学生主要工作：指导教师审核：年 月 日第五、六周时间: （ 月 日 至 月 日）学生主要工作：指导教师审核：年 月 日第七、八周时间: （ 月 日 至 月 日）学生主要工作：指导教师审核：年 月 日第九、十周时间: （ 月 日 至 月 日）学生主要工作：指导教师审核：年 月 日第十一、十二周时间: （ 月 日 至 月 日）学生主要工作：指导教师审核：年 月 日第十三、十四周时间: （ 月 日 至 月 日）学生主要工作：指导教师审核：年 月 日七、其他（学生提交）1开题报告1份 2论文1份（5000字以上） 指 导 教 师： 学 院 院 长： 学生开始执行任务书日期： 学 生 姓 名： 送交毕业设计（论文）日期： 夹具夹紧力的优化及对工件定位精度的影响B.Li 和 S.N.Mellkote布什伍德拉夫机械工程学院，佐治亚理工学院，格鲁吉亚，美国研究所由于夹紧和加工，在工件和夹具的接触部位会产生局部弹性变形，使工件尺寸发生变化，进而影响工件的最终加工质量。这种效应可通过最小化夹具设计优化，夹紧力是一个重要的设计变量，可以得到优化，以减少工件的位移。本文提出了一种确定多夹紧夹具受到准静态加工部位的最佳夹紧力的新方法。该方法采用弹性接触力学模型代表夹具与工件接触，并涉及制定和解决方案的多目标优化模型的约束。夹紧力的最优化对工件定位精度的影响通过3-2-1式铣夹具的例子进行了分析。关键词：弹性 接触 模型 夹具 夹紧力 优化 前言 定位和夹紧的工件加工中的两个关键因素。要实现夹具的这些功能，需将工件定位到一个合适的基准上并夹紧，采用的夹紧力必须足够大，以抑制工件在加工过程中产生的移动。然而，过度的夹紧力可诱导工件产生更大的弹性变形 ，这会影响它的位置精度，并反过来影响零件质量。所以有必要确定最佳夹紧力，来减小由于弹性变形对工件的定位误差，同时满足加工的要求。在夹具分析和综合领域上的研究人员使用了有限元模型的方法或刚体模型的方法。大量的工作都以有限元方法为基础被报道参考文献1-8。随着得墨忒耳8，这种方法的限制是需要较大的模型和计算成本。同时，多数的有限元基础研究人员一直重点关注的夹具布局优化和夹紧力的优化还没有得到充分讨论，也有少数的研究人员通过对刚性模型9-11对夹紧力进行了优化，刚型模型几乎被近似为一个规则完整的形状。得墨忒耳12，13用螺钉理论解决的最低夹紧力，总的问题是制定一个线性规划，其目的是尽量减少在每个定位点调整夹紧力强度的法线接触力。接触摩擦力的影响被忽视，因为它较法线接触力相对较小，由于这种方法是基于刚体假设，独特的三维夹具可以处理超过6个自由度的装夹，复和倪14也提出迭代搜索方法，通过假设已知摩擦力的方向来推导计算最小夹紧力，该刚体分析的主要限制因素是当出现六个以上的接触力是使其静力不确定，因此，这种方法无法确定工件移位的唯一性。 这种限制可以通过计算夹具工件系统15的弹性来克服，对于一个相对严格的工件，该夹具在机械加工工件的位置会受夹具点的局部弹性变形的强烈影响。Hockenberger和得墨忒耳16使用经验的接触力变形的关系（称为元功能），解决由于夹紧和准静态加工力工件刚体位移。同一作者还考察了加工工件夹具位移对设计参数的影响17。桂 18 等 通过工件的夹紧力的优化定位精度弹性接触模型对报告做了改善，然而，他们没有处理计算夹具与工件的接触刚度的方法，此外，其算法的应用没有讨论机械加工刀具路径负载有限序列。李和Melkote 19和乌尔塔多和Melkote 20用接触力学解决由于在加载夹具夹紧点弹性变形产生的接触力和工件的位移，他们还使用此方法制定了优化方法夹具布局21和夹紧力22。但是，关于multiclamp系统及其对工件精度影响的夹紧力的优化并没有在这些文件中提到 。本文提出了一种新的算法，确定了multiclamp夹具工件系统受到准静态加载的最佳夹紧力为基础的弹性方法。该法旨在尽量减少影响由于工件夹紧位移和加工荷载通过系统优化夹紧力的一部分定位精度。接触力学模型，用于确定接触力和位移，然后再用做夹紧力优化，这个问题被作为多目标约束优化问题提出和解决。通过两个例子分析工件夹紧力的优化对定位精度的影响，例子涉及的铣削夹具3-2-1布局。1 夹具工件联系模型 11 模型假设该加工夹具由L定位器和带有球形端的c形夹组成。工件和夹具接触的地方是线性的弹性接触，其他地方完全刚性。工件夹具系统由于夹紧和加工受到准静态负载。夹紧力可假定为在加工过程中保持不变，这个假设是有效的，在对液压或气动夹具使用。在实际中，夹具工件接触区域是弹性分布，然而，这种模式的发展，假设总触刚度（见图1）第i夹具接触力局部变形如下： (1) 其中（j=x，y，z）表示，在当地子坐标系切线和法线方向的接触刚度第 19 页 共 15 页图1 弹簧夹具工件接触模型。 表示在第i个接触处的坐标系（j=x，y，z）是对应沿着xyz方向的弹性变形，分别 （j= x，y，z）的代表和切向力接触 ，法线力接触。12 工件夹具的接触刚度模型集中遵守一个球形尖端定位，夹具和工件的接触并不是线性的，因为接触半径与随法线力呈非线性变化 23。由于法线力接触变形作用于半径和平面工件表面之间，这可从封闭赫兹的办法解决缩进一个球体弹性半空间的问题。对于这个问题， 是法线的变形，在文献23 第93页中给出如下： （2）其中式中 和是工件和夹具的弹性模量，、分别是工件和材料的泊松比。切向变形沿着和切线方向）硅业切力距有以下形式文献23第217页 （3）其中、 分别是工件和夹具剪切模量一个合理的接触刚度的线性可以近似从最小二乘获得适合式 （2），这就产生了以下线性化接触刚度值：在计算上述的线性近似， （4） (5)正常的力被假定为从0到1000N，且最小二乘拟合相应的R2值认定是0.94。2夹紧力优化 我们的目标是确定最优夹紧力，将尽量减少由于工件刚体运动过程中，局部的夹紧和加工负荷引起的弹性变形，同时保持在准静态加工过程中夹具工件系统平衡，工件的位移减少，从而减少定位误差。实现这个目标是通过制定一个多目标约束优化问题的问题，如下描述。2.1 目标函数配方工件旋转，由于部队轮换往往是相当小17的工件定位误差假设为确定其刚体翻译基本上，其中 、和 是 沿，和三个正交组件（见图2）。图2 工件刚体平移和旋转工件的定位误差归于装夹力，然后可以在该刚体位移的范数计算如下： （6）其中表示一个向量二级标准。 但是作用在工件的夹紧力会影响定位误差。当多个夹紧力作用于工件，由此产生的夹紧力为,有如下形式： （7）其中夹紧力是矢量，夹紧力的方向矩阵，是夹紧力是矢量的方向余弦，、和 是第i个夹紧点夹紧力在、和方向上的向量角度（i=1、2、3.,C）。在这个文件中，由于接触区变形造成的工件的定位误差，被假定为受的作用力是法线的，接触的摩擦力相对较小，并在进行分析时忽略了加紧力对工件的定位误差的影响。意指正常接触刚度比，是通过（i=1，2L）和最小的所有定位器正常刚度相乘，并假设工件、取决于、的方向，各自的等效接触刚度可有下式计算得出（见图3），工件刚体运动，归于夹紧行动现在可以写成： (8)工件有位移，因此，定位误差的减小可以通过尽量减少产生的夹紧力向量 范数。因此，第一个目标函数可以写为：最小化 （9）要注意，加权因素是与等效接触刚度成正比的在、和 方向上。通过使用最低总能量互补参考文献15，23的原则求解弹性力学接触问题得出A的组成部分是唯一确定的，这保证了夹紧力和相应的定位反应是“真正的”解决方案，对接触问题和产生的“真正”刚体位移，而且工件保持在静态平衡，通过夹紧力的随时调整。因此，总能量最小化的形式为补充的夹紧力优化的第二个目标函数，并给出：最小化 （10）其中代表机构的弹性变形应变能互补，代表由外部力量和力矩配合完成，是遵守对角矩阵的， 和是所有接触力的载体。如图3 加权系数计算确定的基础内蒙古科技大学本科生毕业设计（外文翻译）2.2 摩擦和静态平衡约束在（10）式优化的目标受到一定的限制和约束，他们中最重要的是在每个接触处的静摩擦力约束。库仑摩擦力的法律规定（是静态摩擦系数）,这方面的一个非线性约束和线性化版本可以使用，并且19有： （11）假设准静态载荷，工件的静力平衡由下列力和力矩平衡方程确保（向量形式）： (12)其中包括在法线和切线方向的力和力矩的机械加工力和工件重量。2.3界接触力由于夹具工件接触是单侧面的，法线的接触力只能被压缩。这通过以下的的约束表（i=1，2,L+C） （13）它假设在工件上的法线力是确定的，此外，在一个法线的接触压力不能超过压工件材料的屈服强度（）。这个约束可写为： （i=1，2，,L+C） (14) 如果是在第i个工件夹具的接触处的接触面积，完整的夹紧力优化模型，可以写成：最小化 (15)3模型算法求解式（15）多目标优化问题可以通过求解约束24。这种方法将确定的目标作为首要职能之一，并将其转换成一个约束对。该补充（）的主要目的是处理功能，并由此得到夹紧力（）作为约束的加权范数最小化。对为主要目标的选择，确保选中一套独特可行的夹紧力，因此，工件夹具系统驱动到一个稳定的状态（即最低能量状态），此状态也表示有最小的夹紧力下的加权范数。 的约束转换涉及到一个指定的加权范数小于或等于，其中是 的约束，假设最初所有夹紧力不明确，要确定一个合适的。在定位和夹紧点的接触力的计算只考虑第一个目标函数（即）。虽然有这样的接触力，并不一定产生最低的夹紧力，这是一个“真正的”可行的解决弹性力学问题办法，可完全抑制工件在夹具中的位置。这些夹紧力的加权系数，通过计算并作为初始值与比较，因此，夹紧力式（15）的优化问题可改写为: 最小化 （16）由： (11)(14) 得。类似的算法寻找一个方程根的二分法来确定最低的上的约束, 通过尽可能降低上限，由此产生的最小夹紧力的加权范数。 迭代次数K，终止搜索取决于所需的预测精度和，有参考文献15： （17）其中表示上限的功能，完整的算法在如图4中给出。 图4 夹紧力的优化算法（在示例1中使用）。图5 该算法在示例2使用4 加工过程中的夹紧力的优化及测定上一节介绍的算法可用于确定单负载作用于工件的载体的最佳夹紧力，然而，刀具路径随磨削量和切割点的不断变化而变化。因此，相应的夹紧力和最佳的加工负荷获得将由图4算法获得，这大大增加了计算负担，并要求为选择的夹紧力提供标准， 将获得满意和适宜的整个刀具轨迹 ，用保守的办法来解决下面将被讨论的问题，考虑一个有限的数目（例如m）沿相应的刀具路径设置的产生m个最佳夹紧力，选择记为， ， ,在每个采样点，考虑以下四个最坏加工负荷向量： (18)、和表示在、和方向上的最大值，、和上的数字1，2，3分别代替对应的和另外两个正交切削分力，而且有：虽然4个最坏情况加工负荷向量不会在工件加工的同一时刻出现，但在每次常规的进给速度中，刀具旋转一次出现一次，负载向量引入的误差可忽略。因此，在这项工作中，四个载体负载适用于同一位置，（但不是同时）对工件进行的采样 ，夹紧力的优化算法图4，对应于每个采样点计算最佳的夹紧力。夹紧力的最佳形式有： (i=1,2,m) (j=x,y z,r) (19)其中是最佳夹紧力的四个情况下的加工负荷载体，(C=1，2，C)是每个相应的夹具在第i个样本点和第j负荷情况下力的大小。是计算每个负载点之后的结果，一套简单的“最佳”夹紧力必须从所有的样本点和装载条件里发现，并在所有的最佳夹紧力中选择。这是通过在所有负载情况和采样点排序，并选择夹紧点的最高值的最佳的夹紧力，见于式 （20）： （k=1，2，C） （20）只要这些具备，就得到一套优化的夹紧力，验证这些力，以确保工件夹具系统的静态平衡。否则，会出现更多采样点和重复上述程序。在这种方式中，可为整个刀具路径确定“最佳”夹紧力 ，图5总结了刚才所描述的算法。请注意，虽然这种方法是保守的，它提供了一个确定的夹紧力，最大限度地减少工件的定位误差的一套系统方法。5影响工件的定位精度它的兴趣在于最早提出了评价夹紧力的算法对工件的定位精度的影响。工件首先放在与夹具接触的基板上，然后夹紧力使工件接触到夹具，因此，局部变形发生在每个工件夹具接触处，使工件在夹具上移位和旋转。随后，准静态加工负荷应用造成工件在夹具的移位。工件刚体运动的定义是由它在、和方向上的移位和自转（见图2），如前所述，工件刚体位移产生于在每个夹紧处的局部变形，假设为相对于工件的质量中心的第i个位置矢量定位点，坐标变换定理可以用来表达在工件的位移，以及工件自转如下: (21)其中表示旋转矩阵,描述当地在第i帧相联系的全球坐标系和是一个旋转矩阵确定工件相对于全球的坐标系的定位坐标系。假设夹具夹紧工件旋转，由于旋转很小，故也可近似为： （22） 方程（21）现在可以改写为： （23）其中是经方程（21）重新编排后变换得到的矩阵式，是夹紧和加工导致的工件刚体运动矢量。工件与夹具单方面接触性质意味着工件与夹具接触处没有拉力的可能。因此，在第i装夹点接触力可能与的关系如下： （24）其中是在第i个接触点由于夹紧和加工负荷造成的变形，意味着净压缩变形，而负数则代表拉伸变形; 是表示在本地坐标系第i个接触刚度矩阵，是单位向量. 在这项研究中假定液压/气动夹具，根据对外加工负荷，故在法线方向的夹紧力的强度保持不变，因此，必须对方程（24）的夹紧点进行修改为： （25）其中是在第i个夹紧点的夹紧力，让表示一个对外加工力量和载体的61矢量。并结合方程（23）（25）与静态平衡方程，得到下面的方程组： （26）其中，其中表示相乘。由于夹紧和加工工件刚体移动，q可通过求解式（26）得到。工件的定位误差向量， （见图6），现在可以计算如下： （27） 其中是考虑工件中心加工点的位置向量，且 6模拟工作 较早前提出的算法是用来确定最佳夹紧力及其对两例工件精度的影响例如：1适用于工件单点力。2应用于工件负载准静态铣削序列 如左图7 工件夹具配置中使用的模拟研究 工件夹具定位联系; 、和全球坐标系。 3-2-1夹具图7所示，是用来定位并控制7075 - T6铝合金（127毫米127毫米38.1毫米）的柱状块。假定为球形布局倾斜硬钢定位器/夹具在表1中给出。工件夹具材料的摩擦静电对系数为0.25。使用伊利诺伊大学开发EMSIM程序参考文献26 对加工瞬时铣削力条件进行了计算，如表2给出例（1），应用工件在点（109.2毫米，25.4毫米，34.3毫米）瞬时加工力，图4中表3和表4列出了初级夹紧力和最佳夹紧力的算法 。该算法如图5所示 ，一个25.4毫米铣槽使用EMSIM进行了数值模拟，以减少起步（0.0毫米，25.4毫米，34.3毫米）和结束时（127.0毫米，25.4毫米，34.3毫米）四种情况下加工负荷载体，（见图8）。模拟计算铣削力数据在表5中给出。图8最终铣削过程模拟例如2。表6中5个坐标列出了为模拟抽样调查点。最佳夹紧力是用前面讨论过的排序算法计算每个采样点和负载载体最后的夹紧力和负载。7结果与讨论例如算法1的绘制最佳夹紧力收敛图9，图9对于固定夹紧装置在图示例假设（见图7），由此得到的夹紧力加权范数有如下形式：.结果表明，最佳夹紧力所述加工条件下有比初步夹紧力强度低得多的加权范数，最初的夹紧力是通过减少工件的夹具系统补充能量算法获得。由于夹紧力和负载造成的工件的定位误差，如表7。结果表明工件旋转小，加工点减少错误从13.1到14.6不等。在这种情况下，所有加工条件改善不是很大，因为从最初通过互补势能确定的最小化的夹紧力值已接近最佳夹紧力。图5算法是用第二例在一个序列应用于铣削负载到工件，他应用于工件铣削负载一个序列。最佳的夹紧力，对应列表6每个样本点，随着最后的最佳夹紧力，在每个采样点的加权范数和最优的初始夹紧力绘图10，在每个采样点的加权范数的，和绘制。结果表明，由于每个组成部分是各相应的最大夹紧力，它具有最高的加权范数。如图10所示，如果在每个夹紧点最大组成部分是用于确定初步夹紧力，则夹紧力需相应设置，有比相当大的加权范数。故是一个完整的刀具路径改进方案。上述模拟结果表明，该方法可用于优化夹紧力相对于初始夹紧力的强度，这种做法将减少所造成的夹紧力的加权范数，因此将提高工件的定位精度。图108结论该文件提出了关于确定多钳夹具，工件受准静态加载系统的优化加工夹紧力的新方法。夹紧力的优化算法是基于接触力学的夹具与工件系统模型，并寻求尽量减少应用到所造成的工件夹紧力的加权范数，得出工件的定位误差。该整体模型，制定一个双目标约束优化问题，使用-约束的方法解决。该算法通过两个模拟表明，涉及3-2-1型，二夹铣夹具的例子。今后的工作将解决在动态负载存在夹具与工件在系统的优化，其中惯性，刚度和阻尼效应在确定工件夹具系统的响应特性具有重要作用。9参考资料：1、J. D. Lee 和L. S. Haynes .柔性夹具系统的有限元分析交易美国ASME，工程杂志工业 ：134-139页。2、W. Cai, S. J. Hu 和J. X. Yuan .“柔性钣金夹具：原理，算法和模拟”，交易美国ASME，制造科学与工程杂志 ：1996 318-324页。3、P. Chandra, S. M. Athavale, R. E. DeVor 和S. G. Kapoor.“负载对表面平整度的影响”工件夹具制造科学研讨会论文集1996，第一卷：146-152页。4、R. J. Menassa 和V. R. DeVries.“适用于选拔夹具设计与优化方法，美国ASME工业工程杂志：113 、 412-414，1991。5、A. J. C. Trappey, C. Su 和J. Hou.计算机辅助夹具分析中的应用有限元分析和数学优化模型， 1995 ASME程序，MED： 777-787页。6、 S. N. Melkote, S. M. Athavale, R. E. DeVor, S. G. Kapoor 和J. Burkey .“基于加工过程仿真的加工装置作用力系统研究”， NAMRI/SME：207214页, 19957、“考虑工件夹具，夹具接触相互作用布局优化模拟的结果” 341-346，1998。 8、E. C. DeMeter. 快速支持布局优化，国际机床制造， 硕士论文 1998。9、Y.-C. Chou, V. Chandru, M. M. Barash .加工夹具机械构造的数学算法：分析和合成，美国ASME，工程学报工业“：1989 299-306页。10、S. H. Lee 和 M. R. Cutkosky. 具有摩擦性的夹具规划 美国ASME，工业工程学报：1991，320327页。11、S. Jeng, L. Chen 和W. Chieng.“最小夹紧力分析”，国际机床制造，硕士论文 1995年。12、E. C. DeMeter.加工夹具的性能的最小最大负荷标准 美国ASME，工业工程杂志 ：199413、E. C. DeMeter .加工夹具最大负荷的性能优化模型 美国ASME，工业工程杂志 1995。14、JH复和AYC倪.“核查和工件夹持的夹具设计”方案优化，设计和制造，4，硕士论文： 307-318，1994。15、T. H. Richards、埃利斯 霍伍德.1977,应力能量方法分析，1977。16、M. J. Hockenberger and E. C. DeMeter. 对工件准静态分析功能位移在加工夹具的应用程序，制造科学杂志与工程： 325331页, 1996。Int J Adv Manuf Technol (2001) 17:104113 2001 Springer-Verlag London Limited Fixture Clamping Force Optimisation and its Impact on Workpiece Location Accuracy B. Li and S. N. Melkote George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Georgia, USA Workpiece motion arising from localised elastic deformation at fixtureworkpiece contacts owing to clamping and machining forces is known to affect significantly the workpiece location accuracy and, hence, the final part quality. This effect can be minimised through fixture design optimisation. The clamping force is a critical design variable that can be optimised to reduce the workpiece motion. This paper presents a new method for determining the optimum clamping forces for a multiple clamp fixture subjected to quasi-static machining forces. The method uses elastic contact mechanics models to represent the fixtureworkpiece contact and involves the formulation and solution of a multi-objective constrained optimisation model. The impact of clamping force optimisation on workpiece location accuracy is analysed through examples involving a 32-1 type milling fixture. Keywords: Elastic contact modelling; Fixture clamping force; Optimisation 1. Introduction The location and immobilisation of the workpiece are two critical factors in machining. A machining fixture achieves these functions by locating the workpiece with respect to a suitable datum, and clamping the workpiece against it. The clamping force applied must be large enough to restrain the workpiece motion completely during machining. However, excessive clamping force can induce unacceptable level of workpiece elastic distortion, which will adversely affect its location and, in turn, the part quality. Hence, it is necessary to determine the optimum clamping forces that minimise the workpiece location error due to elastic deformation while satisfying the total restraint requirement. Previous researchers in the fixture analysis and synthesis area have used the finite-element (FE) modelling approach or Correspondence and offprint requests to: Dr S. N. Melkote, George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0405, USA. E-mail: shreyes.melkoteme.gatech.edu the rigid-body modelling approach. Extensive work based on the FE approach has been reported 18. With the exception of DeMeter 8, a common limitation of this approach is the large model size and computation cost. Also, most of the FE- based research has focused on fixture layout optimisation, and clamping force optimisation has not been addressed adequately. Several researchers have addressed fixture clamping force optimisation based on the rigid-body model 911. The rigid body modelling approach treats the fixture-element and work- piece as perfectly rigid solids. DeMeter 12, 13 used screw theory to solve for the minimum clamping force. The overall problem was formulated as a linear program whose objective was to minimise the normal contact force at each locating point by adjusting the clamping force intensity. The effect of the contact friction force was neglected because of its relatively small magnitude compared with the normal contact force. Since this approach is based on the rigid body assumption, it can uniquely only handle 3D fixturing schemes that involve no more than 6 unknowns. Fuh and Nee 14 also presented an iterative search-based method that computes the minimum clamping force by assuming that the friction force directions are known a priori. The primary limitation of the rigid-body analysis is that it is statically indeterminate when more than six contact forces are unknown. As a result, workpiece displace- ments cannot be determined uniquely by this method. This limitation may be overcome by accounting for the elasticity of the fixtureworkpiece system 15. For a relatively rigid workpiece, the location of the workpiece in the machining fixture is strongly influenced by the localised elastic defor- mation at the fixturing points. Hockenberger and DeMeter 16 used empirical contact force-deformation relations (called meta- functions) to solve for the workpiece rigid-body displacements due to clamping and quasi-static machining forces. The same authors also investigated the effect of machining fixture design parameters on workpiece displacement 17. Gui et al 18 reported an elastic contact model for improving workpiece location accuracy through optimisation of the clamping force. However, they did not address methods for calculating the fixtureworkpiece contact stiffness. In addition, the application of their algorithm for a sequence of machining loads rep- resenting a finite tool path was not discussed. Li and Melkote 19 and Hurtado and Melkote 20 used contact mechanics to Fixture Clamping Force Optimisation 105 solve for the contact forces and workpiece displacement pro- duced by the elastic deformation at the fixturing points owing to clamping loads. They also developed methods for optimising the fixture layout 21 and clamping force using this method 22. However, clamping force optimisation for a multiclamp system and its impact on workpiece accuracy were not covered in these papers. This paper presents a new algorithm based on the contact elasticity method for determining the optimum clamping forces for a multiclamp fixtureworkpiece system subjected to quasi- static loads. The method seeks to minimise the impact of workpiece motion due to clamping and machining loads on the part location accuracy by systematically optimising the clamping forces. A contact mechanics model is used to deter- mine a set of contact forces and displacements, which are then used for the clamping force optimisation. The complete prob- lem is formulated and solved as a multi-objective constrained optimisation problem. The impact of clamping force optimis- ation on workpiece location accuracy is analysed via two examples involving a 32-1 fixture layout for a milling oper- ation. 2. FixtureWorkpiece Contact Modelling 2.1 Modelling Assumptions The machining fixture consists of L locators and C clamps with spherical tips. The workpiece and fixture materials are linearly elastic in the contact region, and perfectly rigid else- where. The workpiecefixture system is subjected to quasi- static loads due to clamping and machining. The clamping force is assumed to be constant during machining. This assumption is valid when hydraulic or pneumatic clamps are used. In reality, the elasticity of the fixtureworkpiece contact region is distributed. However, in this model development, lumped contact stiffness is assumed (see Fig. 1). Therefore, the contact force and localised deformation at the ith fixturing point can be related as follows: F i j = k i j d i j (1) where k i j (j = x,y,z) denotes the contact stiffness in the tangential and normal directions of the local x i ,y i ,z i coordinate frame, d i j Fig. 1. A lumped-spring fixtureworkpiece contact model. x i , y i , z i , denote the local coordinate frame at the ith contact. (j = x,y,z) are the corresponding localised elastic deformations along the x i ,y i , and z i axes, respectively, F i j (j = x,j,z) represents the local contact force components with F i x and F i y being the local x i and y i components of the tangential force, and F i z the normal force. 2.2 WorkpieceFixture Contact Stiffness Model The lumped compliance at a spherical tip locator/clamp and workpiece contact is not linear because the contact radius varies nonlinearly with the normal force 23. The contact deformation due to the normal force P i acting between a spherical tipped fixture element of radius R i and a planar workpiece surface can be obtained from the closed-form Hertz- ian solution to the problem of a sphere indenting an elastic half-space. For this problem, the normal deformation D i n is given as 23, p. 93: D i n = S 9(P i ) 2 16R i (E*) 2 D 1/3 (2) where 1 E* = 1 - n 2 w E w + 1 - n 2 f E f E w and E f are Youngs moduli for the workpiece and fixture materials, respectively, and n w and n f are Poisson ratios for the workpiece and fixture materials, respectively. The tangential deformation D i t (= D i tx or D i ty in the local x i and y i tangential directions, respectively) due to a tangential force Q i (= Q i x or Q i y ) has the following form 23, p. 217: D ti t = Q i 8a i S 2 - n f G f + 2 - n w G w D (3) where a i = S 3P i R i 4 S 1 - n f E f + 1 - n w E w DD 1/3 and G w and G f are shear moduli for the workpiece and fixture materials, respectively. A reasonable linear approximation of the contact stiffness can be obtained from a least-squares fit to Eq. (2). This yields the following linearised contact stiffness values: k i z = 8.82 S 16R i (E*) 2 9 D 1/3 (4) k i x = k i y = 4 E* S 2 - n j G f + 2 - n w G w D - 1 k i z (5) In deriving the above linear approximation, the normal force P i was assumed to vary from 0 to 1000 N, and the correspond- ing R 2 value of the least-squares fit was found to be 0.94. 3. Clamping Force Optimisation The goal is to determine the set of optimal clamping forces that will minimise the workpiece rigid-body motion due to 106 B. Li and S. N. Melkote localised elastic deformation induced by the clamping and machining loads, while maintaining the fixtureworkpiece sys- tem in quasi-static equilibrium during machining. Minimisation of the workpiece motion will, in turn, reduce the location error. This goal is achieved by formulating the problem as a multi- objective constrained optimisation problem, as described next. 3.1 Objective Function Formulation Since the workpiece rotation due to fixturing forces is often quite small 17 the workpiece location error is assumed to be determined largely by its rigid-body translation Dd w = DX w DY w DZ w T , where DX w , DY w , and DZ w are the three orthogonal components of Dd w along the X g , Y g , and Z g axes (see Fig. 2). The workpiece location error due to the fixturing forces can then be calculated in terms of the L 2 norm of the rigid-body displacement as follows: iDd w i = (DX w ) 2 + (DY w ) 2 + (DZ w ) 2 ) (6) where ii denotes the L 2 norm of a vector. In particular, the resultant clamping force acting on the workpiece will adversely affect the location error. When mul- tiple clamping forces are applied to the workpiece, the resultant clamping force, P R C = P R X P R y P R Z T , has the form: P R C = R C P C (7) where P C = P L+1 .P L+C T is the clamping force vector, R C = n L+1 .n L+C T is the clamping force direction matrix, n L+i = cosa L+i cosb L+i cosg L+i T is the clamping force direction cosine vector, and a L+i , b L+i , and g L+i are angles made by the clamping force vector at the ith clamping point with respect to the X g , Y g , Z g coordinate axes (i = 1,2,. . .,C). In this paper, the workpiece location error due to contact region deformation is assumed to be influenced only by the normal force acting at the locatorworkpiece contacts. The frictional force at the contacts is relatively small and is neg- lected when analysing the impact of the clamping force on the workpiece location error. Denoting the ratio of the normal contact stiffness, k i z , to the smallest normal stiffness among all locators, k s z ,byj i (i = 1,. . .,L), and assuming that the workpiece rests on N X , N Y , and N Z number of locators oriented in the X g , Fig. 2. Workpiece rigid body translation and rotation. Y g , and Z g directions, the equivalent contact stiffness in the X g , Y g , and Z g directions can be calculated as k s zSO N X i=1 j iD , k s zSO N Y i=1 j iD , and k s zSO N Z i=1 j iD respectively (see Fig. 3). The workpiece rigid-body motion, Dd w , due to clamping action can now be written as: Dd w = 3 P R X k s zSO N X i=1 j iD P R Y k s zSO N Y i=1 j iD P R Z k s z SO N Z i=1 j iD 4 T (8) The workpiece motion, and hence the location error can be reduced by minimising the weighted L 2 norm of the resultant clamping force vector. Therefore, the first objective function can be written as: Minimize iP R C i w = ! 11 P R X O N X i=1 j i 2 2 + 1 P R Y O N Y i=1 j i 2 2 + 1 P R Z O N Z i=1 j i 2 2 2 (9) Note that the weighting factors are proportional to the equival- ent contact stiffnesses in the X g , Y g , and Z g directions. The components of P R C are uniquely determined by solving the contact elasticity problem using the principle of minimum total complementary energy 15, 23. This ensures that the clamping forces and the corresponding locator reactions are “true” solutions to the contact problem and yield “true” rigid- body displacements, and that the workpiece is kept in static equilibrium by the clamping forces at all times. Therefore, the minimisation of the total complementary energy forms the second objective function for the clamping force optimisation and is given by: Minimise (U* - W*) = 1 2 FO L+C i=1 (F i x ) 2 k i x + O L+C i=1 (F i y ) 2 k i y + O L+C i=1 (F i z ) 2 k i z G (10) = .l T Ql Fig. 3. The basis for the determination of the weighting factor for the L 2 norm calculation. Fixture Clamping Force Optimisation 107 where U* represents the complementary strain energy of the elastically deformed bodies, W* represents the complementary work done by the external force and moments, Q = diag c 1 x c 1 y c 1 z .c L+C x c L+C y c L+C z is the diagonal contact compliance matrix, c i j = (k i j ) - 1 , and l = F 1 x F 1 y F 1 z .F L+C x F L+C y F L+C z T is the vector of all contact forces. 3.2 Friction and Static Equilibrium Constraints The optimisation objective in Eq. (10) is subject to certain constraints and bounds. Foremost among them is the static friction constraint at each contact. Coulombs friction law states that (F i x ) 2 +(F i y ) 2 ) #m i s F i z (m i s is the static friction coefficient). A conservative and linearised version of this nonlinear con- straint can be used and is given by 19: uF i x u + uF i y u #m i s F i z (11) Since quasi-static loads are assumed, the static equilibrium of the workpiece is ensured by including the following force and moment equilibrium equations (in vector form): O F = 0 (12) O M = 0 where the forces and moments consist of the machining forces, workpiece weight and the contact forces in the normal and tangential directions. 3.3 Bounds Since the fixtureworkpiece contact is strictly unilateral, the normal contact force, P i , can only be compressive. This is expressed by the following bound on P i : P i $ 0(i = 1, . . ., L + C) (13) where it is assumed that normal forces directed into the workpiece are positive. In addition, the normal compressive stress at a contact cannot exceed the compressive yield strength (S y ) of the workpiece material. This upper bound is written as: P i # S y A i (i = 1, . . .,L+C) (14) where A i is the contact area at the ith workpiecefixture con- tact. The complete clamping force optimisation model can now be written as: Minimize f = H f 1 f 2 J = H .l T Ql iP R C i w J (15) subject to: (11)(14). 4. Algorithm for Model Solution The multi-objective optimisation problem in Eq. (15) can be solved by the e-constraint method 24. This method identifies one of the objective functions as primary, and converts the other into a constraint. In this work, the minimisation of the complementary energy (f 1 ) is treated as the primary objective function, and the weighted L 2 norm of the resultant clamping force (f 2 ) is treated as a constraint. The choice of f 1 as the primary objective ensures that a unique set of feasible clamping forces is selected. As a result, the workpiecefixture system is driven to a stable state (i.e. the minimum energy state) that also has the smallest weighted L 2 norm for the resultant clamping force. The conversion of f 2 into a constraint involves specifying the weighted L 2 norm to be less than or equal to e, where e is an upper bound on f 2 . To determine a suitable e,itis initially assumed that all clamping forces are unknown. The contact forces at the locating and clamping points are computed by considering only the first objective function (i.e. f 1 ). While this set of contact forces does not necessarily yield the lowest clamping forces, it is a “true” feasible solution to the contact elasticity problem that can completely restrain the workpiece in the fixture. The weighted L 2 norm of these clamping forces is computed and taken as the initial value of e. Therefore, the clamping force optimisation problem in Eq. (15) can be rewritten as: Minimize f 1 = .l T Ql (16) subject to: iP R C i w $e, (11)(14). An algorithm similar to the bisection method for finding roots of an equation is used to determine the lowest upper bound for iP R C i w . By decreasing the upper bound e as much as possible, the minimum weighted L 2 norm of the resultant clamping force is obtained. The number of iterations, K, needed to terminate the search depends on the required prediction accuracy d and ueu, and is given by 25: K = F log 2 S ueu d DG (17) where I denotes the ceiling function. The complete algorithm is given in Fig. 4. 5. Determination of Optimum Clamping Forces During Machining The algorithm presented in the previous section can be used to determine the optimum clamping force for a single load vector applied to the workpiece. However, during milling the magnitude and point of cutting force application changes continuously along the tool path. Therefore, an infinite set of optimum clamping forces corresponding to the infinite set of machining loads will be obtained with the algorithm of Fig. 4. This substantially increases the computational burden and calls for a criterion/procedure for selecting a single set of clamping forces that will be satisfactory and optimum for the entire tool path. A conservative approach to addressing these issues is discussed next. Consider a finite number (say m) of sample points along the tool path yielding m corresponding sets of optimum clamp- ing forces denoted as P 1 opt , P 2 opt ,.,P m opt . At each sampling 108 B. Li and S. N. Melkote Fig. 4. Clamping force optimisation algorithm (used in example 1). point, the following four worst-case machining load vectors are considered: F X max = F max X F 1 Y F 1 Z T F Y max = F 2 X F max Y F 2 Z T F Z max = F 3 X F 3 Y F max Z T (18) F r max = F 4 X F 4 Y F 4 Z T where F max X , F max Y , and F max Z are the maximum X g , Y g , and Z g components of the machining force, the superscripts 1, 2, 3 of F X , F Y , and F Z stand for the other two orthogonal machining force components corresponding to F max X , F max Y , and F max Z , respectively, and iF r max i = max(F X ) 2 +(F Y ) 2 +(F Z ) 2 ). Although the four worst-case machining load vectors will not act on the workpiece at the same instant, they will occur once per cutter revolution. At conventional feedrates, the error introduced by applying the load vectors at the same point would be negligible. Therefore, in this work, the four load vectors are applied at the same location (but not simultaneously) on the workpiece corresponding to the sam- pling instant. The clamping force optimisation algorithm of Fig. 4 is then used to calculate the optimum clamping forces corresponding to each sampling point. The optimum clamping forces have the form: P i jmax = C i 1j C i 2j .C i Cj T (i = 1, . . .,m)(j = x,y,z,r) (19) where P i jmax is the vector of optimum clamping forces for the four worst-case machining load vectors, and C i kj (k = 1,. . .,C) is the force magnitude at each clamp corresponding to the ith sample point and the jth load scenario. After P i jmax is computed for each load application point, a single set of “optimum” clamping forces must be selected from all of the optimum clamping forces found for each clamp from all the sample points and loading conditions. This is done by sorting the optimum clamping force magnitudes at a clamping point for all load scenarios and sample points and selecting the maximum value, C max k , as given in Eq. (20): C max k # C i kj (k =