CN113742859A - 一种基于齿面残余应力与接触应力的复合应力场构建方法 - Google Patents
一种基于齿面残余应力与接触应力的复合应力场构建方法 Download PDFInfo
- Publication number
- CN113742859A CN113742859A CN202110929491.8A CN202110929491A CN113742859A CN 113742859 A CN113742859 A CN 113742859A CN 202110929491 A CN202110929491 A CN 202110929491A CN 113742859 A CN113742859 A CN 113742859A
- Authority
- CN
- China
- Prior art keywords
- stress
- tooth surface
- gear
- residual stress
- tooth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 20
- 238000010276 construction Methods 0.000 title claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 16
- 238000012360 testing method Methods 0.000 claims abstract description 16
- 238000002441 X-ray diffraction Methods 0.000 claims abstract description 6
- 238000009826 distribution Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 5
- 238000005260 corrosion Methods 0.000 claims description 3
- 230000007797 corrosion Effects 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 7
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 230000006872 improvement Effects 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 238000005422 blasting Methods 0.000 description 2
- 238000005255 carburizing Methods 0.000 description 2
- 229920001971 elastomer Polymers 0.000 description 2
- 239000000806 elastomer Substances 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 210000003781 tooth socket Anatomy 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种基于齿面残余应力与接触应力的复合应力场构建方法。本发明考虑了加工过程中引入的残余应力,采用X射线衍射残余应力测试,构建齿面残余应力场;基于齿轮有限元接触分析,计算齿面接触应力场;进行含残余应力的齿轮接触分析,构建齿面复合应力场;得到的齿轮应力场更接近与实际状态,使得对于齿轮强度和耐久度的评估更加准确。
Description
技术领域
本发明属于齿轮评估技术领域,尤其涉及一种基于齿面残余应力与接触应力的复合应力场构建方法。
背景技术
齿轮的耐久性分析大都建立在精确的应力分析上,即计算齿轮的应力场,但是现有的应力场构建通常仅仅考虑齿轮啮合过程中的接触应力,而没有考虑到齿轮加工过程中出现的残余应力,导致建立的应力场不够精准。具体的国际上广泛应用的AGMA与ISO齿轮接触疲劳强度计算标准均基于Hertz接触理论,而Hertz接触理论通过把两弹性体的接触问题等效为一个当量弹性体和一个刚性平面的接触,对复杂载荷下齿面接触应力计算精度较低。工程实际中没有对齿轮表面与亚表面的残余应力场进行考虑,不能准确反映齿面残余应力与接触应力的复合应力状态。因此得到的应力场与实际情况差别较大。
其中,残余应力是指在没有外加载荷的作用下,存在于结构内部的应力分布,会对结构的力学性能与疲劳寿命产生较大影响。其中,齿轮加工过程中的渗碳淬火-磨削-喷丸工艺,会在齿轮表面与亚表面产生残余应力,而齿轮在啮合过程中会受到残余应力与接触应力的复合作用。
发明内容
为解决上述问题,本发明公开一种基于齿面残余应力与接触应力的复合应力场构建方法,本发明考虑了加工过程中引入的残余应力,采用X射线衍射残余应力测试,构建齿面残余应力场;基于齿轮有限元接触分析,计算齿面接触应力场;进行含残余应力的齿轮接触分析,构建齿面复合应力场;得到的齿轮应力场更接近与实际状态,使得对于齿轮强度和耐久度的评估更加准确。
一种基于齿面残余应力与接触应力的复合应力场构建方法,包括以下步骤:
步骤一、输入齿轮几何参数,建立齿轮三维几何模型;
步骤二、输入齿轮材料参数,导入齿轮三维几何模型,对齿面接触区域进行网格细化,构建齿轮三维有限元接触分析模型;
步骤三、得到齿轮的试验样本残余应力沿齿厚方向分布曲线;
步骤四、基于齿轮几何模型,将齿面离散成齿面节点,结合试验样本残余应力沿齿厚方向分布曲线,插值拟合各测量点的正应力与切应力,得到在齿面坐标系下齿面上各齿面节点的残余应力;
步骤五、结合各齿面节点的方向余弦矩阵,将齿面坐标系下的残余应力转换到全局坐标系中,得到在全局坐标系下的各齿面节点的残余应力:
σn2=[βn]-1[σn1]([βn]Τ)-1 (1.1)
其中,σn2为全局坐标系下第n个齿面节点的残余应力,σn1为齿面坐标系下第n个齿面节点的残余应力,βn为第n个齿面节点的方向余弦,n表示节点编号;T表示矩阵转置;
步骤六、结合齿轮三维有限元接触分析模型,导入全局坐标系下的各齿面节点的残余应力,施加齿面残余应力场;
步骤七、设置边界条件与相互作用,进行含残余应力场的齿轮接触分析,构建齿轮表面与亚表面的复合应力场。
进一步的改进,所述步骤一中,输入齿轮几何参数,采用CATIA平台建立齿轮三维几何模型;所述齿轮几何参数包括齿数、大端模数、轴交角、中点螺旋角、法向压力角、齿宽和外锥距。
进一步的改进,所述步骤二中,采用ABAQUS平台构建齿轮三维接触分析模型,其中齿轮材料参数包括材料类型、弹性模量、泊松比和密度。
进一步的改进,所述步骤三中,采用X射线衍射残余应力测试与电解腐蚀剥层,得到齿轮的试验样本残余应力沿齿厚方向分布曲线。
进一步的改进,所述步骤六中,使用VISUAL STUDIO编程语言对各齿面节点施加不同的残余应力场。
进一步的改进,所述边界条件包括主动轮的转速与从动轮所受转矩,所述相互作用指两个啮合齿轮之间的摩擦系数。
本发明的优点:
基于X射线衍射残余应力测试,得到齿轮表面与亚表面的残余应力分布,构建的齿面残余应力场更能反映实际工作情况。构建含残余应力场的三维有限元接触分析模型,可更准确的反映齿面残余应力与接触应力的复合应力状态,使得对于齿轮强度和耐久度的评估更加准确。
附图说明
图1齿轮有限元接触分析模型;
图2齿面残余应力测试值;
图3齿面复合应力场;
图4齿面节点正应力历程曲线。
具体实施方式
下面结合附图以及实施例对本发明进行更为详细的说明。
步骤一、螺旋锥齿轮几何参数如表1所示,构建齿轮齿坯模型。结合产形面方程与啮合方程,计算螺旋锥齿轮齿面离散点坐标。输入齿轮几何参数,通过MATLAB编程计算齿轮各齿面离散点坐标。导入CATIA平台生成齿面离散点,采用三次样条曲线拟合齿面离散点。拟合齿面样条曲线,并建立齿根曲面和齿根过渡曲面,形成整个齿槽曲面。
表1螺旋锥齿轮几何参数
根据螺旋锥齿轮的面角、节角、根角、外锥距及齿面宽画出草图,通过旋转命令建立齿轮齿坯。基本建立的齿槽曲面,采用曲面特征分割命令,获得螺旋锥齿轮的齿槽实体,再通过旋转阵列及布尔操作可获得完整的齿轮几何模型。
步骤二、齿轮材料参数见表2。基于CATIA平台将构建的螺旋锥齿轮切割成单齿模型;为了合理布置网格密度,通过特征曲线构建分割曲面,再通过分割命令将单齿模型分割为10个区域。
表2齿轮材料属性
将分割后的单齿模型导入ABAQUS平台中,确认各分割区域能够进行结构化网格划分。网格单元类型选择线性缩减积分单元C3D8R。对不同分割区域分别进行结构化网格划分,为提高计算精度,对齿面接触表面与亚表面区域进行细化网格,构建齿轮接触分析有限元模型见图1。
步骤三、基于X射线衍射法与电解腐蚀剥层法,测量经渗碳淬火-磨削-喷丸工艺齿轮试件的表面与亚表面残余应力,得到距齿面不同深度的残余应力分布曲线如图2所示。
步骤四、基于齿轮几何模型,将齿面离散成21×41个节点。结合小样试件与齿轮试件残余应力测试数据,插值拟合各测量点的正应力与切应力,得到在齿面坐标系下齿面上任一点的残余应力。
步骤五、结合各齿面节点的方向余弦矩阵,将齿面坐标系下的残余应力转换到全局坐标系中。有限元接触分析中,只能施加全局坐标系下的残余应力场,故需要将齿面坐标系下的残余应力转换到全局坐标系中。结合各节点的方向余弦矩阵,通过下式得到在全局坐标系下的齿面节点的残余应力。
σn2=[βn]-1[σn1]([βn]Τ)-1
步骤六、结合三维有限元接触分析模型,导入齿面节点残余应力信息,采用VISUALSTUDIO编程语言对各齿面节点施加不同的残余应力场
步骤七、设置主动轮转速为2000r/min,从动轮转矩为1200N·m,各相互作用的摩擦系数为0.11。进行含残余应力场的齿轮接触分析,构建齿轮表面与亚表面的复合应力场,见图3。分别绘制齿面接触应力场与齿面复合应力场下齿面节点的正应力历程曲线见图4。
上述仅为本发明的一个具体导向实施方式,但本发明的设计构思并不局限于此,凡利用此构思对本发明进行非实质性的改动,均应属于侵犯本发明的保护范围的行为。
Claims (6)
1.一种基于齿面残余应力与接触应力的复合应力场构建方法,其特征在于,包括以下步骤:
步骤一、输入齿轮几何参数,建立齿轮三维几何模型;
步骤二、输入齿轮材料参数,导入齿轮三维几何模型,对齿面接触区域进行网格细化,构建齿轮三维有限元接触分析模型;
步骤三、得到齿轮的试验样本残余应力沿齿厚方向分布曲线;
步骤四、基于齿轮几何模型,将齿面离散成齿面节点,结合试验样本残余应力沿齿厚方向分布曲线,插值拟合各测量点的正应力与切应力,得到在齿面坐标系下齿面上各齿面节点的残余应力;
步骤五、结合各齿面节点的方向余弦矩阵,将齿面坐标系下的残余应力转换到全局坐标系中,得到在全局坐标系下的各齿面节点的残余应力:
σn2=[βn]-1[σn1]([βn]Τ)-1 (0.1)
其中,σn2为全局坐标系下第n个齿面节点的残余应力,σn1为齿面坐标系下第n个齿面节点的残余应力,βn为第n个齿面节点的方向余弦,n表示节点编号;T表示矩阵转置;
步骤六、结合齿轮三维有限元接触分析模型,导入全局坐标系下的各齿面节点的残余应力,施加齿面残余应力场;
步骤七、设置边界条件与相互作用,进行含残余应力场的齿轮接触分析,构建齿轮表面与亚表面的复合应力场。
2.如如权利要求1所述的基于齿面残余应力与接触应力的复合应力场构建方法,其特征在于,所述步骤一中,输入齿轮几何参数,采用CATIA平台建立齿轮三维几何模型;所述齿轮几何参数包括齿数、大端模数、轴交角、中点螺旋角、法向压力角、齿宽和外锥距。
3.如如权利要求1所述的基于齿面残余应力与接触应力的复合应力场构建方法,其特征在于,所述步骤二中,采用ABAQUS平台构建齿轮三维接触分析模型,其中齿轮材料参数包括材料类型、弹性模量、泊松比和密度。
4.如如权利要求1所述的基于齿面残余应力与接触应力的复合应力场构建方法,其特征在于,所述步骤三中,采用X射线衍射残余应力测试与电解腐蚀剥层,得到齿轮的试验样本残余应力沿齿厚方向分布曲线。
5.如如权利要求1所述的基于齿面残余应力与接触应力的复合应力场构建方法,其特征在于,所述步骤六中,使用VISUAL STUDIO编程语言对各齿面节点施加不同的残余应力场。
6.如如权利要求1所述的基于齿面残余应力与接触应力的复合应力场构建方法,其特征在于,所述边界条件包括主动轮的转速与从动轮所受转矩,所述相互作用指两个啮合齿轮之间的摩擦系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110929491.8A CN113742859B (zh) | 2021-08-13 | 2021-08-13 | 一种基于齿面残余应力与接触应力的复合应力场构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110929491.8A CN113742859B (zh) | 2021-08-13 | 2021-08-13 | 一种基于齿面残余应力与接触应力的复合应力场构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113742859A true CN113742859A (zh) | 2021-12-03 |
CN113742859B CN113742859B (zh) | 2024-03-08 |
Family
ID=78731013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110929491.8A Active CN113742859B (zh) | 2021-08-13 | 2021-08-13 | 一种基于齿面残余应力与接触应力的复合应力场构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113742859B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003194637A (ja) * | 2001-12-26 | 2003-07-09 | Toshiba Corp | 有限要素法による残留応力解析方法 |
US20080141782A1 (en) * | 2004-12-16 | 2008-06-19 | Frontics, Inc. | Evaluating Method of the Residual Stress Determining Method Using the Continuous Indentation Method |
JP2009236822A (ja) * | 2008-03-28 | 2009-10-15 | Komatsu Ltd | 歯車の圧縮残留応力評価方法および歯車設計方法 |
CN106503350A (zh) * | 2016-10-25 | 2017-03-15 | 北京航空航天大学 | 一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的设计方法 |
CN106649971A (zh) * | 2016-10-25 | 2017-05-10 | 北京航空航天大学 | 一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 |
CN107335848A (zh) * | 2017-06-20 | 2017-11-10 | 西北工业大学 | 三维铣削残余应力预测方法 |
CN109271711A (zh) * | 2018-09-25 | 2019-01-25 | 重庆大学 | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 |
CN110008594A (zh) * | 2019-04-08 | 2019-07-12 | 长安大学 | 一种行星传动人字齿轮有限元网格自动建模与装配方法 |
CN112069618A (zh) * | 2020-09-03 | 2020-12-11 | 上海理工大学 | 解决齿轮接触分析仿真的方法 |
CN112948989A (zh) * | 2021-01-26 | 2021-06-11 | 浙江大学 | 一种数据驱动的齿面强化方法 |
-
2021
- 2021-08-13 CN CN202110929491.8A patent/CN113742859B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003194637A (ja) * | 2001-12-26 | 2003-07-09 | Toshiba Corp | 有限要素法による残留応力解析方法 |
US20080141782A1 (en) * | 2004-12-16 | 2008-06-19 | Frontics, Inc. | Evaluating Method of the Residual Stress Determining Method Using the Continuous Indentation Method |
JP2009236822A (ja) * | 2008-03-28 | 2009-10-15 | Komatsu Ltd | 歯車の圧縮残留応力評価方法および歯車設計方法 |
CN106503350A (zh) * | 2016-10-25 | 2017-03-15 | 北京航空航天大学 | 一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的设计方法 |
CN106649971A (zh) * | 2016-10-25 | 2017-05-10 | 北京航空航天大学 | 一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 |
CN107335848A (zh) * | 2017-06-20 | 2017-11-10 | 西北工业大学 | 三维铣削残余应力预测方法 |
CN109271711A (zh) * | 2018-09-25 | 2019-01-25 | 重庆大学 | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 |
CN110008594A (zh) * | 2019-04-08 | 2019-07-12 | 长安大学 | 一种行星传动人字齿轮有限元网格自动建模与装配方法 |
CN112069618A (zh) * | 2020-09-03 | 2020-12-11 | 上海理工大学 | 解决齿轮接触分析仿真的方法 |
CN112948989A (zh) * | 2021-01-26 | 2021-06-11 | 浙江大学 | 一种数据驱动的齿面强化方法 |
Non-Patent Citations (2)
Title |
---|
刘晓龙;高玉魁;刘蕴韬;陈东风;: "孔挤压强化残余应力场的三维有限元模拟和实验研究", 航空材料学报, no. 02, 1 April 2011 (2011-04-01), pages 24 - 27 * |
刘竹丽;陈赟;王祝新;: "齿面残余应力对齿轮轮齿弯曲疲劳寿命的影响分析", 郑州大学学报(工学版), no. 03, 28 June 2020 (2020-06-28), pages 53 - 56 * |
Also Published As
Publication number | Publication date |
---|---|
CN113742859B (zh) | 2024-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021721B (zh) | 一种渗碳圆柱齿轮参数化仿真分析方法及cae系统 | |
CN110008594B (zh) | 一种行星传动人字齿轮有限元网格自动建模与装配方法 | |
Jyothirmai et al. | A finite element approach to bending, contact and fatigue stress distribution in helical gear systems | |
CN107436982A (zh) | 考虑基体刚度修正的剥落斜齿轮副的啮合特性分析方法 | |
CN109271711B (zh) | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 | |
CN104408241A (zh) | 一种修形圆柱齿轮的有限元网格自动生成方法 | |
CN111666643A (zh) | 一种用于确定复杂齿面接触性能的方法 | |
Zhang et al. | Design, meshing characteristics and stress analysis of cylindrical gears with curvilinear tooth profile | |
CN113742859B (zh) | 一种基于齿面残余应力与接触应力的复合应力场构建方法 | |
Bodzás | Computer aided designing and modelling of x-zero gear drive | |
CN115026354A (zh) | 一种复杂齿形的车齿刀具逆向包络设计方法 | |
CN113192180B (zh) | 一种基于插齿加工原理的椭圆齿轮参数化精确建模方法 | |
CN111075920B (zh) | 基于fft和润滑影响的rv减速器摆线针轮残余应力求解方法 | |
Batsch | Rapid prototyping and tooth contact analysis of eccentric cycloid gear mesh | |
Liu et al. | Meshing analysis and optimization for plane-generated enveloping toroid hourglass worm drive | |
Singh et al. | Comparative study of stress analysis of helical gear using AGMA Standards and FEM | |
CN109190324B (zh) | 一种弧齿锥齿轮含误差的齿面载荷接触变形数值计算方法 | |
CN113010978A (zh) | 一种基于动态仿真的航空直齿轮修形方法 | |
CN113704907A (zh) | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 | |
Maeda | Research on simplification of face gear tooth surface design | |
CN110188475B (zh) | 一种基于刚轮参数的工作状态下的柔轮测量与评价方法 | |
Wu et al. | The aeronautics face-gear NC hobbing machining technology | |
Murthy et al. | Stress analysis of helical gear by FEM techniques with variation in face width And helix angle | |
CN113868729B (zh) | 一种基于共轭齿廓曲线作图法的精确圆柱蜗轮建模方法 | |
Shaker | Optimization of tooth-root profile for maximum load-carrying capacity: spur and bevel gears |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |