基于COSMOSWorks的直齿圆柱齿轮齿根应力分析

    4.1建立一个静态分析项目
   
    在"COSMOSWorks管理程序"里，右键单击"齿轮"图标，在弹出的菜单中选择"算例"，"网格类型"选择为"实体网格"，"类型"选择为"static(静态)"，单击"确定"按钮完成算例项目的创建。
   
    4.2 定义齿轮材料
   
    单击"应用材料到所选零部件"图标，弹出的"材料"属性管理器，在"选择材料来源"栏中选择"自定义"，根据表2输入各项参数，"模型类型"选择为"线性弹性同向性"。
   

    4.3定义约束
   
    单击"约束"图标，弹出的"约束"属性管理器，选择约束"类型"为"固定"，选择图2中的齿轮键槽面为约束面。
   
    4.4施加栽荷
   
    为了在齿顶上施加啮合力，需在齿轮实体上建立一个过齿顶线且与基圆相切的基准面，然后单击"载荷"图标，弹出"力"属性管理器，在"类型"上选择"应用力、力矩"，在"力的面、边线、顶点"上选择齿顶线，在"方向的面、边线、基准面、基准轴"中选择上述基准面，在"沿基准面方向1"中设置力的大小为1882．41N(法向力大小计算见4．2节)，在"沿基准面方向2"和"沿基准面方向3"中不设置力的大小，单击"确定"按钮，完成载荷施加。
   
    4．5划分网格
   
    单击"网格"图标，弹出的"网格"属性管理器。COSMOSWorks可根据零件情况自动划分网格的形状及大小，这里采用四面体单元格，网格大小为8．36mm，误差为0．418mm，共划分33788个单元格，54383个节点。
   
    4.6分析计算
   
    以上参数设置完成后，单击"运行"图标，系统开始分析计算并弹出运行进展显示框。由于COSMOSWorks采用独有的快速有限元技术(FFE)和精确计算技术，其分析计算的速度比较快。静态分析运行成功后，COSMOSWorks将在管理器中生成"结果"文件夹，含有"应力"、"位移"、"应变"和"变形"等图标。
   
    4.7结果显示
   
    轮齿应力分布云图如图4所示，图形的颜色变化反映了齿轮内部应力的分布情况。从图中可以看出，最大应力发生在轮齿受拉侧的齿根处，为97．2N／mm2；而在轮齿的另一侧，齿根处受压，最大压应力为-22．6N／mm2，显然远小于拉应力，故轮齿受拉侧先产生疲劳裂纹。
   

    5齿根弯曲应力的理论计算
   
    5.1齿轮工作转矩的计算
   
    由于齿轮传递功率P=4kW，n1=960r／min，则齿轮工作转矩T1为：T1=9550000P/n1
   
    将已知数据代入式(2)求出T1=39800N·mm。
   
    5．2齿轮啮合时法向力R的计算
   
    根据齿轮传动的平稳性可知，齿轮在啮合过程中，法向力R的大小与方向不变，因而轮齿在齿顶处的法向力可由节点处的法向力求出。
   
    Fn=2T1/d1cosa(3)
   
    式中：a为分度圆压力角。将T1和已知条件代入式(3)，得Fn=1882．41N。
   
    5．3验算轮齿弯曲应力及讨论
   
    查表得Yp=3．75(Z=90，x=0)，ψd=1，为了与COSMOSWorks仿真条件相同，这里不考虑载荷集中和动载荷的影响，即Kβ=1，Kv=1代入式(1)得齿根处最大弯曲应力σF=98．3N／mm2。而COSMOSWorks的仿真值为97.2N／mm2，误差为1.1％。这说明：①仿真结果是正确的；②在齿根弯曲应力的理论计算中，将轮齿假设为悬臂梁，采用300近似法确定齿根处的危险截面是合理的，完全满足机械设计的要求；③由于理论值略大于仿真值，因此理论计算是略偏于安全的。
   
    上述计算与仿真都是假设全部载荷由一对轮齿来承担，即重叠系数εa=1，根据重叠系数计算公式
    

    可知，实际重叠系数占εa=1.75，这表明，在轮齿转过一个基圆齿距的时间里，两对轮齿同时啮合的时间占75％，此时齿根处的最大应力为48.6N／mm2，即97.2/2N／mm2；而一对轮齿啮合的时间占25％，此时齿根处的最大应力为97．2N／mm2，因此在轮齿啮合过程中，齿根处最大应力达97．2N／mm2只占整个啮合时间的1／4。
   
    6结论
   
    通过齿根弯曲应力的有限元分析，为齿轮的强度设计提供了参考依据，也对齿轮的进一步优化设计提供重要的指导意义。由于COSMOSWorks具有操作简单，运算速度快和计算精度高等优点，而且分析的模型与结果与SolidWorks共享一个数据库，从CAD到CAE无需进行转换，因此在机械零件的设计中有着广泛的应用前景，也是机械设计人员必备的技能之一。