我们采用北京大学力学系研制的SAP5P软件，分析计算麻花钻切削时的变形情况。

1. 等参元的选择

   SAP5P程序提供了多种单元类型。由于麻花钻为三维实体结构，故选用八类节点数可变的三维等参单元，每个单元根据需要取8～21个节点进行计算。

2. 单元的划分

   麻花钻切削锥部为主要受力体，载荷主要作用在主切削刃上，因此，划分单元时可将切削锥部划分得细一些，而导向部分可划分得粗一些，并从锥部与导向部分连接处开始逐渐由密到疏。麻花钻的单元划分如图3所示。由于切削锥部近似为圆锥体，而横刃仅为一条直线，因此必然存在单元间的过渡连接问题。如图4所示，在锥部中间构造一节点(节点28)，过此节点与横刃、两后刀面及两沟槽构造一六面体，将此六面体划分为8个单元，在锥部与导向部分连接处作一横截面，在此横截面与六面体之间构造8个单元，这样切削锥部共划分为16个单元。
   

   图3 麻花钻单元划分图

   图4 切削锥部单元划分图

   图5 主切削刃载荷分布图

   导向部分为带两条螺旋沟槽的圆柱体。沿轴向截取不同高度的实体段，每段划分8个20节点单元。根据计算机容量和麻花钻切削长度选择相应的实体段数及单元数。

3. 节点数学模型的建立

   根据麻花钻切削锥部的实际结构计算节点坐标。

   导向部分外径基本相同，只是实体段上层相对于下层沿螺旋线旋转了一定角度。设任一单元下层节点坐标为X(i)，Y(i)，Z(i)；上层节点坐标为X(j)，Y(j)，Z(j)，则

   X(j)=X(i)cosf - Y(i)sinf

   Y(j)=X(i)cosf + Y(i)sinf

   Z(j)=Z(i)+h
   其中，f＝(h/r)tgw，h为单元厚度(即实体段高度)，r为钻头外圆处半径，w为螺旋角，�b为单元下层节点沿螺旋线旋转的角度。

4. 载荷的简化及计算

   麻花钻的受力情况较复杂，因此可对载荷作一些简化。下面分两种加载方式进行计算。

   1. 载荷以集中力的形式作用于主切削刃最外端。假设麻花钻只受扭矩作用，设扭矩为M(N·mm)，麻花钻最大外径为D(mm)，则作用在主切削刃最外端的力F为

      F＝M/D

   2. 载荷按实际情况作用于主切削刃相应节点上，如图5所示。

      设P为单位半径长度载荷，则
      

      图6 有限元计算流程图

      图7 麻花钻变形图

      切向力 Pt＝9．81HBf/2
      法向力 Pn＝(0.5～1.0)Pt
      径向力 Pr＝Pncosf
      轴向力 Pa＝Pnsinf

   其中，HB为被加工工件材料的硬度，f为钻削进给量(mm/r)。则主切削刃相应节点所分配的载荷近似为

   P1＝PR/6

   P2＝4PR/6

   P3＝PR/6

5. 约束条件

   采用SAP5P程序对麻花钻进行分析计算时，只对切削部分及导向部分进行计算，忽略了柄部的影响。因此，在整体分析中，应对柄部与导向部分连接处的节点加上固定端约束条件，限制其六个方向的自由度；对其它节点则给以X，Y，Z三个方向的平移自由度，限制X，Y，Z三个方向的转动自由度。

6. 前处理通用接口程序

   我们用Borland C语言编写了前处理通用接口程序，用户只需输入麻花钻的主要几何参数，通过运行前处理程序，即可自动划分网格，自动生成供SAP5P运行的数据文件。有限元计算流程如图6所示。

三、计算结果

选取如表1所示的三组麻花钻参数进行计算。弹性模量E＝2.15×105/mm2，剪切弹性模量G＝8.05×104N/mm2，波桑系数µ＝0.25，载荷计算取HB＝200，f＝2.5mm/r，Pn＝0.75Pt。

将三组参数输入前处理程序，生成供SAP5P运行的数据文件，进入SAP5P进行计算。麻花钻变形图见图7。由输出位移及变形图可知，最大变形发生在主切削刃最外端，即节点32(56)处。

以直径D＝20mm的麻花钻为例，分别改变芯厚及螺旋角的大小，其它参数不变，进入SAP5P进行计算，结果如表2所示。图8所示为D0-d变化曲线，图9所示为D0-w变化曲线。

四、结束语

本文以麻花钻为例，采用有限元法分析了螺旋刃刀具主要几何参数的变化对扭转刚度的影响。从计算结果可以看出，麻花钻几何参数不同，其变形量及刚度值也不同；随着钻芯直径及螺旋角的增大，变形量逐渐减小，刚度值逐渐增大，从而可优选出使变形量最小、刚度最大的刀具几何参数，从而提高刀具的切削性能