CN104361220B - 一种考虑热固耦合的齿轮系时变应力场计算方法 - Google Patents
一种考虑热固耦合的齿轮系时变应力场计算方法 Download PDFInfo
- Publication number
- CN104361220B CN104361220B CN201410602264.4A CN201410602264A CN104361220B CN 104361220 B CN104361220 B CN 104361220B CN 201410602264 A CN201410602264 A CN 201410602264A CN 104361220 B CN104361220 B CN 104361220B
- Authority
- CN
- China
- Prior art keywords
- gear
- overbar
- heat production
- gear train
- field
- 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.)
- Active
Links
Landscapes
- Gears, Cams (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种考虑热固耦合的齿轮系时变应力场计算方法,它有七大步骤:一、确定齿轮系的啮合产热模型;二、确定经验模型的各项参数;三、计算产热及分配产热功率;四、传热分析计算稳态温度场;五、稳态温度场解析表达;六、离散齿轮系的一个旋转周期及热固耦合静力计算;七、插值计算一个循环周期内的应力场。本发明可以计算齿轮系的稳态温度场,并可以进行热固耦合下的应力场计算。它提供了一种将复杂物理过程简化组合的思路,对于其他力学仿真具有重要的借鉴意义。
Description
技术领域
本发明涉及一种考虑热固耦合的齿轮系时变应力场计算方法,它是用与齿轮系耐久性分析的应力分析方法,本方法主要考虑了两个问题,第一是齿轮系的摩擦生热带来的热固耦合问题;第二是考虑了齿轮系啮合传递过程中的应力时变问题。属于机械产品耐久性设计分析技术领域。
背景技术
机械产品的耐久性分析大都建立在精确的应力分析基础上,可以说应力分析的准确与否直接影响耐久性的分析精度与可信性。目前,齿轮系的应力分析大都局限于静力分析,由于齿轮系在啮合过程中还要受到自身啮合产热的影响(尤其对于重载齿轮系更加明显),所以在温度场对应力场的影响是不可忽略的。另一方面传统的强度校核过程,取危险截面进行强度校核即可,但耐久性分析要求得到工作剖面上完整的时变应力场。
本发明提供了一种考虑热固耦合的齿轮系时变应力场计算方法,首先在传统的计算静应力分布的基础上,以齿轮系统为研究对象,综合考虑齿轮啮合产热及向环境散热过程得到稳态温度场的影响,得到热固耦合情况下的静力分析结果。其次考虑齿轮啮合过程的周期性,离散齿轮啮合过程,以离散点为样本点进行热固耦合计算,最后以样本点作为基本数据进行插值计算,得到齿轮啮合周期的考虑了热固耦合的时变的应力场。
发明内容
本发明的目的是提供一种考虑热固耦合的齿轮系时变应力场计算方法。
本发明是通过以下技术方案实现的,首先选择齿轮系的产热经验模型,然后对齿轮系的进行具体分析,得到齿轮系的产热相关参数。另一方面根据齿轮系的实际工况和润滑状态得到另一部分产热相关参数。计算齿轮系的产热功率后,将所得产热功率进行分配,得到每个齿轮的产热功率。然后对齿轮系进行Fluent流体传热仿真,得到齿轮系的稳态温度场。将温度场作为载荷施加到静力分析上,得到热固耦合情况下的应力结果。在静力分析上要同时施加转速,以便得到考虑离心力情况下的计算结果。将齿轮的一个选择周期进行离散化,每一个离散点作为一个样本点分别进行静力分析。然后对样本点的计算结果进行插值计算得到整个旋转周期内的时变应力场。
本发明是一种考虑热固耦合的齿轮系时变应力场计算方法,该方法具体步骤如下:
步骤一:确定齿轮系的啮合产热模型
齿轮系啮合产热有多种模型,目前最为常见的是根据油膜模型推倒的理论模型和在实践中总结出的经验模型。由于经验模型参数简单精度良好所以本发明选择Andersonand Loewenthal经验模型。
步骤二:确定经验模型的各项参数
Anderson and Loewenthal经验模型的参数分为三部分:齿轮对形状参数、润滑状态参数和工况参数。所以要确定这些参数需要对齿轮对、润滑状态和工况进行分析。需要确定的参数有八个,分别是齿宽、传递功率、主动齿轮转速、压力角、主动齿轮齿数、从动轮齿数、润滑油动力粘度以及模数。
步骤三:计算产热及分配产热功率
根据产热模型分别计算每一对齿轮的产热功率。假设一对齿轮的啮合产热过程中热量是均匀传导到两个齿轮上。把齿轮对的产热功率均分到两个齿轮上,得到每一个齿轮的产热功率。
步骤四:传热分析计算稳态温度场
齿轮系大都包括润滑和降温系统,所以齿轮系在产热时还包括传热过程并最终达到一个稳态温度场。本发明采用Fluent有限元软件对传热进行分析。
步骤五:稳态温度场解析表达
在Fluent热传导分析步骤中,对齿轮的形状进行了相应的简化,简化成了相应的圆盘,最后得到的稳态温度场是圆盘半径的函数。通过对稳态温度场的插值计算得到相应的解析的温度场的表达式,并以此作为整个齿轮的温度分布函数。
步骤六:离散齿轮系的一个旋转周期及热固耦合静力计算
齿轮系的应力场具有周期性,其周期性是本齿轮的一个旋转周期。所以计算一个旋转周期内的应力场即可。把一个旋转周期离散为若干样本点分别进行计算。步骤五得到的稳态温度场作为载荷施加到静力分析模型上,同时应该施加转速载荷以便考虑齿轮系离心力的作用。进行应力仿真得到具体啮合状况下的应力场分布。
步骤七:插值计算一个循环周期内的应力场
以上计算的到的应力场都是样本点的计算结果,将样本点的计算结果进行插值计算,即可得到循环周期内的时变应力场。
其中,步骤三、步骤五以及步骤七所述的操作均通过MATLAB编程计算完成。
其中,步骤六通过ANSYS有限元软件进行计算,热固耦合应力仿真计算的步骤包括:设置接触方式,模型网格划分,施加扭矩、温度载荷、旋转载荷及约束,应力结果分析。
本发明是一种考虑热固耦合的齿轮系时变应力场计算方法,具有以下优点:
1.本发明可以计算齿轮系的稳态温度场,并可以进行热固耦合下的应力场计算。
2.本发明通过离散旋转周期,计算样本点的应力场,并对样本点结果进行插值,可得到一个周期内的应力场变化情况,通过静态的应力场得到了时变的应力场。
3.本发明是通过计算样本点的应力场得到时变的应力场的,由于样本点的计算是静力计算,其计算速度快、精度高,而且计算方法易于普通工程师掌握。
4.本发明提供了一种将复杂物理过程简化组合的思路,对于其他力学仿真具有重要的借鉴意义。一种考虑热固耦合的齿轮系时变应力场计算方法。
附图说明
图1为本发明的流程框图
具体实施方式
下面将结合具体的某齿轮系考虑热固耦合的时变应力场计算过程对本发明做进一步详细说明。
本发明是一种考虑热固耦合的齿轮系时变应力场计算方法,见图1所示,其具体实施步骤如下:
步骤一:齿轮系啮合产热模型的确定
齿轮系啮合产热有多种模型模型,目前最为常见的是根据油膜模型推倒的理论模型和在实践中总结出的经验模型。由于经验模型参数简单精度良好所以本发明选择Anderson and Loewenthal经验模型。Anderson and Loewenthal法将齿轮的功率分为三部分:滑动、滚动和风阻损失。由于风阻损失较小,忽略风阻损失。滑动损失及滚动损失分别由以下公式确定:
式中,为齿轮的平均法向载荷;和分别为齿轮的平均滑动速度和平均滚动速度;f为摩擦系数;为平均油膜厚度;B,ε和β(直齿轮中此参数为零)分别为齿宽、重合度和基圆螺旋角。
式中,Tp,np分别为齿轮传递转矩和主动齿轮转速;Dp,α分别为齿轮的分度圆直径和端面压力角;ls为齿轮啮合线长度;z1、z2分别为主、被动齿轮齿数;μ为润滑油的动力粘度;Req为当量接触半径。
齿轮啮合总功率损失为:
步骤二:确定步骤一所确定的经验模型的各项参数
Anderson and Loewenthal经验模型的参数分为三部分:齿轮对形状参数、润滑状态参数和工况参数。所以要确定这些参数需要对齿轮对、润滑状态和工况进行分析。需要确定的参数有八个分别是齿宽、传递功率、主动齿轮转速、压力角、主动齿轮齿数、从动轮齿数、润滑油动力粘度以及模数。各项参数确定需要结合齿轮及润滑油相应的国标、行标或实验值。
取一对齿轮啮合,其参数如下表。
步骤三:计算产热及分配产热功率
根据产热模型分别计算每一对齿轮的产热功率。假设一对齿轮的啮合产热过程中热量是均匀传导到两个齿轮上。把齿轮对的产热功率均分到两个齿轮上,得到每一个齿轮的产热功率。根据Anderson and Loewenthal经验模型编制MATLAB程序,将步骤二的参数代入可得产热单元1的产热功率为28.6KW,产热单元2的产热功率为28.4KW。把产热单元的产热功率均分到产热单元中两个齿轮上,并进行热功率叠加,得1齿轮产热功率为14.3KW,2齿轮的产热功率为28.5KW,3齿轮的产热功率为14.2KW。
步骤四:传热分析计算稳态温度场
齿轮系大都包括润滑和降温系统,所以齿轮系在产热时还包括传热过程并最终达到一个稳态温度场。本发明采用Fluent有限元软件对传热进行分析。
步骤五:稳态温度场解析表达
在Fluent热传导分析步骤中,对齿轮的形状进行了相应的简化,简化成了相应的圆盘,最后得到的稳态温度场是圆盘半径的函数。通过对稳态温度场的多项式拟合计算得到相应的解析的温度场的表达式,并以此作为整个齿轮的温度分布函数。通过对步骤四得到的温度场进行多项式拟合得到温度场分别为,T(R1)=2018×R1,T(R2)=3457×R2,T(R3)=1989×R3。
步骤六:离散齿轮系的一个旋转周期及热固耦合静力计算、
齿轮系的应力场具有周期性,其周期性是本齿轮的一个旋转周期。所以计算一个旋转周期内的应力场即可。把一个旋转周期离散为若干样本点分别进行计算。步骤五得到的温度场作为载荷施加到静力分析模型上,同时应该施加转速载荷以便考虑齿轮系离心力的作用。进行应力仿真得到具体啮合状况下的应力场分布。
本例只计算1和2齿轮一个齿从进入啮合到打开的过程。离散点角间距为两度。扭矩为700N.M,温度场取步骤五的计算结果,转速分别为8015r/min和6338r/min。
计算结果取前后齿根处的平均应力。计算结果如下:
步骤七:插值计算一个循环周期内的应力场
以上计算的到的应力场都是样本点的计算结果,将样本点的计算结果进行插值计算,及可得到循环周期内的时变应力场。
Claims (1)
1.一种考虑热固耦合的齿轮系时变应力场计算方法,其特征在于:该方法具体步骤如下:
步骤一:齿轮系啮合产热模型的确定
选择Anderson and Loewenthal经验模型;Anderson and Loewenthal法将齿轮的功率分为三部分:滑动、滚动和风阻损失;由于风阻损失较小,忽略风阻损失;滑动损失及滚动损失分别由以下公式确定:
式中,为齿轮的平均法向载荷;和分别为齿轮的平均滑动速度和平均滚动速度;f为摩擦系数;为平均油膜厚度;B,ε和β分别为齿宽、重合度和基圆螺旋角;
式中,Tp,np分别为齿轮传递转矩和主动齿轮转速;Dp,α分别为齿轮的分度圆直径和端面压力角;ls为齿轮啮合线长度;z1、z2分别为主、被动齿轮齿数;μ为润滑油的动力粘度;Req为当量接触半径;
齿轮啮合总功率损失为:
步骤二:确定步骤一所确定的经验模型的各项参数
Anderson and Loewenthal经验模型的参数分为三部分:齿轮对形状参数、润滑状态参数和工况参数;所以要确定这些参数需要对齿轮对、润滑状态和工况进行分析;需要确定的参数有八个分别是齿宽、传递功率、主动齿轮转速、压力角、主动齿轮齿数、从动轮齿数、润滑油动力粘度以及模数;各项参数确定需要结合齿轮及润滑油相应的国标、行标或实验值;
取一对齿轮啮合,其参数如下:
产热单元1中,包含齿轮编号1/2,齿宽为0.009,传递功率为3,主动齿轮转速为8015,压力角为0.349,主动齿轮齿数为34,从动轮齿数为43,从动轮齿数为0.0029,模数为2.54;产热单元2中,包含齿轮编号2/3,齿宽为0.009,传递功率为3,主动齿轮转速为6338,压力角为0.349,主动齿轮齿数为43,从动轮齿数为34,从动轮齿数为0.0029,模数为2.54;
步骤三:计算产热及分配产热功率
根据产热模型分别计算每一对齿轮的产热功率;假设一对齿轮的啮合产热过程中热量是均匀传导到两个齿轮上;把齿轮对的产热功率均分到两个齿轮上,得到每一个齿轮的产热功率;根据Anderson and Loewenthal经验模型编制MATLAB程序,将步骤二的参数代入可得产热单元1的产热功率为28.6KW,产热单元2的产热功率为28.4KW;把产热单元的产热功率均分到产热单元中两个齿轮上,并进行热功率叠加,得1齿轮产热功率为14.3KW,2齿轮的产热功率为28.5KW,3齿轮的产热功率为14.2KW;
步骤四:传热分析计算稳态温度场
齿轮系包括润滑和降温系统,所以齿轮系在产热时还包括传热过程并最终达到一个稳态温度场;采用Fluent有限元软件对传热进行分析;
步骤五:稳态温度场解析表达
在Fluent热传导分析步骤中,对齿轮的形状进行了相应的简化,简化成了相应的圆盘,最后得到的稳态温度场是圆盘半径的函数;通过对稳态温度场的多项式拟合计算得到相应的解析的温度场的表达式,并以此作为整个齿轮的温度分布函数;通过对步骤四得到的温度场进行多项式拟合得到温度场分别为,T(R1)=2018×R1,T(R2)=3457×R2,T(R3)=1989×R3;
步骤六:离散齿轮系的一个旋转周期及热固耦合静力计算
齿轮系的应力场具有周期性,其周期性是齿轮的一个旋转周期;所以计算一个旋转周期内的应力场即可;把一个旋转周期离散为若干样本点分别进行计算;步骤五得到的温度场作为载荷施加到静力分析模型上,同时应该施加转速载荷以便考虑齿轮系离心力的作用;进行应力仿真得到具体啮合状况下的应力场分布;
步骤七:插值计算一个循环周期内的应力场
以上计算的到的应力场都是样本点的计算结果,将样本点的计算结果进行插值计算,及可得到循环周期内的时变应力场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410602264.4A CN104361220B (zh) | 2014-10-31 | 2014-10-31 | 一种考虑热固耦合的齿轮系时变应力场计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410602264.4A CN104361220B (zh) | 2014-10-31 | 2014-10-31 | 一种考虑热固耦合的齿轮系时变应力场计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104361220A CN104361220A (zh) | 2015-02-18 |
CN104361220B true CN104361220B (zh) | 2017-06-23 |
Family
ID=52528479
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410602264.4A Active CN104361220B (zh) | 2014-10-31 | 2014-10-31 | 一种考虑热固耦合的齿轮系时变应力场计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104361220B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104776191B (zh) * | 2015-04-02 | 2017-02-22 | 清华大学 | 一种飞溅润滑建模过程中齿轮模型的处理方法 |
CN105136322A (zh) * | 2015-09-15 | 2015-12-09 | 安徽精科检测技术有限公司 | 一种齿轮齿根温度测量方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08211946A (ja) * | 1995-02-08 | 1996-08-20 | Hitachi Ltd | 温度制御方法及び装置 |
-
2014
- 2014-10-31 CN CN201410602264.4A patent/CN104361220B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08211946A (ja) * | 1995-02-08 | 1996-08-20 | Hitachi Ltd | 温度制御方法及び装置 |
Non-Patent Citations (2)
Title |
---|
基于热网络法的直升机主减速器热分析;杨攀;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20140715(第7期);第7-8、26-29页 * |
齿轮齿根动应力分析及其结构优化设计;戴进;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20081215(第12期);第52-64页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104361220A (zh) | 2015-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | A mesh stiffness method using slice coupling for spur gear pairs with misalignment and lead crown relief | |
Mastrone et al. | Oil distribution and churning losses of gearboxes: Experimental and numerical analysis | |
Li et al. | Thermal analysis of helical gear transmission system considering machining and installation error | |
Zhan et al. | A CAD-FEM-QSA integration technique for determining the time-varying meshing stiffness of gear pairs | |
CN108052760B (zh) | 一种齿轮副非线性动力学计算方法 | |
Saxena et al. | Effect of shaft misalignment and friction force on time varying mesh stiffness of spur gear pair | |
CN110147624A (zh) | 一种基于载荷谱的齿轮接触疲劳寿命预测方法 | |
CN107133421A (zh) | 一种湿式双离合器自动变速箱齿轮敲击仿真方法及系统 | |
Wang et al. | Investigation into the meshing friction heat generation and transient thermal characteristics of spiral bevel gears | |
Mao et al. | Polymer gear surface thermal wear and its performance prediction | |
CN104502096A (zh) | 齿轮动态传动误差计算方法 | |
Zhang et al. | Determination of surface temperature rise with the coupled thermo-elasto-hydrodynamic analysis of spiral bevel gears | |
Ruzek et al. | On windage losses in high-speed pinion-gear pairs | |
Concli et al. | Influence of lubricant temperature, lubricant level and rotational speed on the churning power losses in an industrial planetary speed reducer: computational and experimental study | |
CN104361220B (zh) | 一种考虑热固耦合的齿轮系时变应力场计算方法 | |
Dai et al. | Determination of an optimal oil jet nozzle layout for helical gear lubrication: Mathematical modeling, numerical simulation, and experimental validation | |
CN103726888A (zh) | 一种汽轮机转子低周疲劳在线监控方法 | |
Lukac et al. | Transmission system power flow model | |
Li et al. | Effects of the helix angle, the friction coefficient and mechanical errors on unsteady-state temperature field of helical gear and thermal sensitivity analysis | |
CN105787149A (zh) | 一种由弧齿锥齿轮传动系统轴上功率谱向齿面应力谱精确转换的方法 | |
Sun et al. | Analytical models for thermal deformation and mesh stiffness of spur gears under steady temperature field | |
Li et al. | Analysis of thermal characteristic of spur/helical gear transmission | |
Guo et al. | A new algorithm to solve meshing-in impact considering the measured pitch error and to investigate its influence on the dynamic characteristics of a gear system | |
Liu et al. | Investigation on contact and bending stress of face-hobbed and face-milled hypoid gear | |
CN203365034U (zh) | 一种带齿轮转盘轴承装配及摩擦力矩测量机 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |