技术领域
[0001] 本发明属于齿轮评估技术领域，尤其涉及一种基于齿面残余应力与接触应力的复合应力场构建方法。
相关背景技术
[0002] 齿轮的耐久性分析大都建立在精确的应力分析上，即计算齿轮的应力场，但是现有的应力场构建通常仅仅考虑齿轮啮合过程中的接触应力，而没有考虑到齿轮加工过程中出现的残余应力，导致建立的应力场不够精准。具体的国际上广泛应用的AGMA与ISO齿轮接触疲劳强度计算标准均基于Hertz接触理论，而Hertz接触理论通过把两弹性体的接触问题等效为一个当量弹性体和一个刚性平面的接触，对复杂载荷下齿面接触应力计算精度较低。工程实际中没有对齿轮表面与亚表面的残余应力场进行考虑，不能准确反映齿面残余应力与接触应力的复合应力状态。因此得到的应力场与实际情况差别较大。
[0003] 其中，残余应力是指在没有外加载荷的作用下，存在于结构内部的应力分布，会对结构的力学性能与疲劳寿命产生较大影响。其中，齿轮加工过程中的渗碳淬火‑磨削‑喷丸工艺，会在齿轮表面与亚表面产生残余应力，而齿轮在啮合过程中会受到残余应力与接触应力的复合作用。
具体实施方式
[0027] 下面结合附图以及实施例对本发明进行更为详细的说明。
[0028] 步骤一、螺旋锥齿轮几何参数如表1所示，构建齿轮齿坯模型。结合产形面方程与啮合方程，计算螺旋锥齿轮齿面离散点坐标。输入齿轮几何参数，通过MATLAB编程计算齿轮各齿面离散点坐标。导入CATIA平台生成齿面离散点，采用三次样条曲线拟合齿面离散点。拟合齿面样条曲线，并建立齿根曲面和齿根过渡曲面，形成整个齿槽曲面。
[0029] 表1螺旋锥齿轮几何参数
[0030]
[0031] 根据螺旋锥齿轮的面角、节角、根角、外锥距及齿面宽画出草图，通过旋转命令建立齿轮齿坯。基本建立的齿槽曲面，采用曲面特征分割命令，获得螺旋锥齿轮的齿槽实体，再通过旋转阵列及布尔操作可获得完整的齿轮几何模型。
[0032] 步骤二、齿轮材料参数见表2。基于CATIA平台将构建的螺旋锥齿轮切割成单齿模型；为了合理布置网格密度，通过特征曲线构建分割曲面，再通过分割命令将单齿模型分割为10个区域。
[0033] 表2齿轮材料属性
[0034]
[0035] 将分割后的单齿模型导入ABAQUS平台中，确认各分割区域能够进行结构化网格划分。网格单元类型选择线性缩减积分单元C3D8R。对不同分割区域分别进行结构化网格划分，为提高计算精度，对齿面接触表面与亚表面区域进行细化网格，构建齿轮接触分析有限元模型见图1。
[0036] 步骤三、基于X射线衍射法与电解腐蚀剥层法，测量经渗碳淬火‑磨削‑喷丸工艺齿轮试件的表面与亚表面残余应力，得到距齿面不同深度的残余应力分布曲线如图2所示。
[0037] 步骤四、基于齿轮几何模型，将齿面离散成21×41个节点。结合小样试件与齿轮试件残余应力测试数据，插值拟合各测量点的正应力与切应力，得到在齿面坐标系下齿面上任一点的残余应力。
[0038] 步骤五、结合各齿面节点的方向余弦矩阵，将齿面坐标系下的残余应力转换到全局坐标系中。有限元接触分析中，只能施加全局坐标系下的残余应力场，故需要将齿面坐标系下的残余应力转换到全局坐标系中。结合各节点的方向余弦矩阵，通过下式得到在全局坐标系下的齿面节点的残余应力。
[0039] σn2＝[βn]‑1[σn1]([βn]Τ)‑1
[0040] 步骤六、结合三维有限元接触分析模型，导入齿面节点残余应力信息，采用VISUAL STUDIO编程语言对各齿面节点施加不同的残余应力场
[0041] 步骤七、设置主动轮转速为2000r/min，从动轮转矩为1200N·m，各相互作用的摩擦系数为0.11。进行含残余应力场的齿轮接触分析，构建齿轮表面与亚表面的复合应力场，见图3。分别绘制齿面接触应力场与齿面复合应力场下齿面节点的正应力历程曲线见图4。
[0042] 上述仅为本发明的一个具体导向实施方式，但本发明的设计构思并不局限于此，凡利用此构思对本发明进行非实质性的改动，均应属于侵犯本发明的保护范围的行为。
