CN106528991B - 基于Taylor随机有限元对齿轮箱进行优化设计的方法 - Google Patents
基于Taylor随机有限元对齿轮箱进行优化设计的方法 Download PDFInfo
- Publication number
- CN106528991B CN106528991B CN201610955932.0A CN201610955932A CN106528991B CN 106528991 B CN106528991 B CN 106528991B CN 201610955932 A CN201610955932 A CN 201610955932A CN 106528991 B CN106528991 B CN 106528991B
- Authority
- CN
- China
- Prior art keywords
- gear
- variance
- mean value
- allowable
- gearbox
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000005452 bending Methods 0.000 claims abstract description 15
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 230000005540 biological transmission Effects 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 3
- 239000000463 material Substances 0.000 abstract description 2
- 230000007423 decrease Effects 0.000 abstract 2
- 230000002860 competitive effect Effects 0.000 abstract 1
- 239000011159 matrix material Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000002994 raw material Substances 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Gears, Cams (AREA)
Abstract
本发明提出基于Taylor随机有限元对齿轮箱进行优化设计的方法,对齿轮箱进行优化设计。核心方案为:建立优化数学模型,目标函数为各级齿轮传动实际中心距之和,计算齿轮弯曲疲劳强度的约束条件,轴强度的约束条件需要用到Taylor随机有限元,齿轮使用的网格采用二十节点六面体参数单元,轴使用的网格采用轴对称四边形环形单元,采用随机方向法求解优化问题,编写计算机运算程序,最后运行计算机运算程序获得最优解。本发明的优化效果显著,齿轮箱质量下降,体积减小,原材料费用下降,提高了产品质量,使产品更具有竞争力。
Description
技术领域
本发明提出基于Taylor随机有限元对齿轮箱进行优化设计的方法,机械设计、机械优化设计、机械现代设计方法领域。
背景技术
齿轮箱在机床,工程机械,冶金机械,矿山机械,石油机械,农业机械,车辆等领域有广泛的应用。随着计算机技术的发展,产生了机械优化设计方法。机械可靠性设计与优化设计相结合形成了可靠性优化设计,既能定量地预测产品的可靠性,又能使产品的设计参数获得优化解。机械可靠性设计仅能对简单零件进行设计。很多现代的结构***具有很高的结构复杂度。在随机的载荷和工作环境下,先进的数值技术、著名的有限元方法被用来分析结构。绝大多数的应用被限制在确定的载荷和工作环境下,尽管随机和不确定的因素达到相当的程度。随机因素对结构的影响越来越受到国内外学者的重视。随着人类认识的深入,忽略随机性的有限元是不符合实际的。有限元分析要想提高计算精度,必须考虑随机因素的影响。考虑随机因素的有限元称为随机有限元。随机有限元的计算方法主要有直接Monte Carlo法,Taylor展开法,摄动法,Neumann展开法,Neumann-PCG法等。
目前,还没有出现基于Taylor随机有限元对齿轮箱进行优化设计的方案。
发明内容
本发明提出基于Taylor随机有限元对齿轮箱进行优化设计的方法,对齿轮箱进行优化设计,使齿轮箱的重量减轻、提高了产品质量。
为此,本发明的技术方案如下:基于Taylor随机有限元对齿轮箱进行优化设计的方法,包括如下步骤:
(1)计算齿轮弯曲疲劳强度的约束条件,用Taylor随机有限元计算齿轮弯曲应力的均值和方差,轴危险截面应力的均值和方差;齿轮使用的网格采用二十节点六面体参数单元,轴使用的网格采用轴对称四边形环形单元,采用随机方向法求解优化问题,生成有限元模型;求解齿轮弯曲疲劳强度的许用均值和许用方差、轴强度的许用均值和许用方差;
(2)建立齿轮箱的优化数学模型
设计变量为:齿轮模数、齿轮齿数、轴的直径、轴的长度;
约束条件为:齿轮弯曲应力的均值和方差小于齿轮弯曲疲劳强度的许用均值和许用方差、轴危险截面应力的均值和方差小于轴强度的许用均值和许用方差;
目标函数为:各级齿轮传动实际中心距之和;
从而立建立齿轮箱的优化数学模型;
(3)根据齿轮箱的优化数学模型,编写计算机运算程序,最后运行计算机运算程序获得最优解。
有益效果:本发明基于Taylor随机有限元对齿轮箱进行优化设计优化,优化效果显著,优化后齿轮箱质量下降、体积减小,原材料费用下降,提高了产品质量,使产品更具有竞争力。
附图说明
图1是一种需要进行优化设计的齿轮箱结构图。
图2所示Taylor随机有限元计算齿轮弯曲应力的均值和方差的框图。
具体实施方式
通过以下内容,进一步详细描述本发明。
图1是一种需要进行优化设计的齿轮箱结构,有12个齿轮,4根轴,标号1-12表达的是齿轮,标号Ⅰ、Ⅱ、Ⅲ、Ⅳ表达的是轴。
结合图1详细描述本发明,基于Taylor随机有限元对齿轮箱进行优化设计的方法,包括如下步骤:
(1)通过齿轮箱原始设计图纸的参数,利用三维建模软件构建齿轮箱的三维实体模型;
(2)将步骤(1)得到的齿轮箱的三维实体模型导入到有限元软件中,计算齿轮弯曲疲劳强度的约束条件,齿轮使用的网格采用二十节点六面体参数单元,轴使用的网格采用轴对称四边形环形单元,采用随机方向法求解优化问题,生成有限元模型;求解齿轮弯曲疲劳强度的许用均值和许用方差、轴强度的许用均值和许用方差;
所述得到均值和方差函数的详细过程:
Taylor随机有限元
静载下有限元支配方程可以写为[K]{δ}={F}
[K]为整体刚度矩阵,{δ}为各结点位移列阵,{F}为各结点载荷列阵。
机械零件材料性能参数,几何尺寸,受到的载荷被看作正态随机变量n个随机变量a1,a2,…,ai,…,an。
上式对ai求偏导数,得
上式对aj求偏导数,得
将位移在随机变量的均值点处进行展开,并在两边同时取均值,可得
(E{δ}表示{δ}的均值,Cov(ai,aj)为ai和aj的协方差。{δ}的方差可由下式计算
{σ}=[D][B]{δ}
{σ}为应力矩阵,[D]为弹性矩阵,[B]为应变矩阵;
上式对ai求偏导数,得
上式对aj求偏导数,得
将应力在随机变量的均值点处进行展开,并在两边同时取均值,可得
E{σ}表示{σ}的均值,Cov(ai,aj)为ai和aj的协方差。{σ}的方差可由下式计算
(3)建立齿轮箱的优化数学模型
图1中的齿轮箱由12个齿轮和4根轴组成;
为了清楚表达下述函数,把图1中标号Ⅰ、Ⅱ、Ⅲ、Ⅳ表达的是轴用1、2、3、4替代;
设计变量为:x=(m1,z1,z2,m2,z3,z4,m3,z5,z6,m4,z7,z8,m5,z9,z10,m6,z11,z12,b1,b2,b3,b4,d1,l1,d2,l2,d3,l3,d4,l4)T,其中,m为齿轮模数,Z为齿轮齿数,d为轴的直径,l为轴的长度;
目标函数为:各级齿轮传动实际中心距之和,具体为:
约束条件为
其中,为齿轮弯曲应力的均值和方差。为齿轮弯曲疲劳强度的许用均值和许用方差。
其中,为轴危险截面应力的均值和方差,为轴强度的许用均值和许用方差。
mkl≤mk≤mks(k=1,2,…,6)
zkl≤zk≤zks(k=1,2,…,12)
bkl≤bk≤bks(k=1,2,3,4)
dkl≤dk≤dks(k=1,2,3,4)
lkl≤lk≤lks(k=1,2,3,4)
其中,mkl,zkl,bkl,dkl,lkl为设计变量下界值。mks,zks,bks,dks,lks为设计变量上界值;
(4)根据齿轮箱的优化数学模型,编写计算机运算程序,最后运行计算机运算程序获得最优解。
下表1是图1所示齿轮箱的原始设计与优化设计参数比较;
表1设计参数比较
m<sub>1</sub> | m<sub>2</sub> | m<sub>3</sub> | m<sub>4</sub> | m<sub>5</sub> | m<sub>6</sub> | z<sub>1</sub> | z<sub>2</sub> | z<sub>3</sub> | z<sub>4</sub> | z<sub>5</sub> | z<sub>6</sub> | z<sub>7</sub> | z<sub>8</sub> | z<sub>9</sub> | |
原始设计 | 4 | 4 | 4 | 4 | 4 | 4 | 18 | 44 | 27 | 43 | 35 | 35 | 31 | 39 | 25 |
优化设计 | 3 | 3.5 | 3.5 | 3.5 | 4 | 4 | 20 | 40 | 28 | 44 | 36 | 36 | 32 | 40 | 27 |
z<sub>10</sub> | z<sub>11</sub> | z<sub>12</sub> | b<sub>1</sub> | b<sub>2</sub> | b<sub>3</sub> | b<sub>4</sub> | d<sub>1</sub> | l<sub>1</sub> | d<sub>2</sub> | l<sub>2</sub> | d<sub>3</sub> | l<sub>3</sub> | d<sub>4</sub> | l<sub>4</sub> | |
原始设计 | 41 | 19 | 47 | 25 | 25 | 25 | 25 | 50 | 350 | 50 | 280 | 50 | 340 | 65 | 290 |
优化设计 | 40 | 19 | 42 | 18 | 22 | 25 | 25 | 45 | 285 | 46 | 210 | 48 | 286 | 62 | 235 |
从表1可以看出,优化效果十分显著,齿轮箱质量下降,体积减小,原材料费用下降,提高了产品质量。
Claims (1)
1.基于Taylor随机有限元对齿轮箱进行优化设计的方法,包括如下步骤:
(1)通过齿轮箱原始设计图纸的参数,利用三维建模软件构建齿轮箱的三维实体模型;
(2)将步骤(1)得到的齿轮箱的三维实体模型导入到有限元软件中,用Taylor随机有限元计算齿轮弯曲应力的均值和方差,轴危险截面应力的均值和方差;齿轮使用的网格采用二十节点六面体参数单元,轴使用的网格采用轴对称四边形环形单元,采用随机方向法求解优化问题,生成有限元模型;求解齿轮弯曲疲劳强度的许用均值和许用方差、轴强度的许用均值和许用方差;
(3)建立齿轮箱的优化数学模型
设计变量为:齿轮模数、齿轮齿数、轴的直径、轴的长度;
约束条件为:齿轮弯曲应力的均值和方差小于齿轮弯曲疲劳强度的许用均值和许用方差、轴危险截面应力的均值和方差小于轴强度的许用均值和许用方差;
目标函数为:各级齿轮传动实际中心距之和;
从而立建立齿轮箱的优化数学模型;
(4)根据齿轮箱的优化数学模型,编写计算机运算程序,最后运行计算机运算程序获得最优解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610955932.0A CN106528991B (zh) | 2016-10-27 | 2016-10-27 | 基于Taylor随机有限元对齿轮箱进行优化设计的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610955932.0A CN106528991B (zh) | 2016-10-27 | 2016-10-27 | 基于Taylor随机有限元对齿轮箱进行优化设计的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106528991A CN106528991A (zh) | 2017-03-22 |
CN106528991B true CN106528991B (zh) | 2019-07-09 |
Family
ID=58325608
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610955932.0A Active CN106528991B (zh) | 2016-10-27 | 2016-10-27 | 基于Taylor随机有限元对齿轮箱进行优化设计的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106528991B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107992688B (zh) * | 2017-12-07 | 2021-09-24 | 中车大连机车车辆有限公司 | 机车牵引齿轮自动设计设备及自动设计方法 |
CN111090961B (zh) * | 2020-01-17 | 2022-05-20 | 湖北汽车工业学院 | 基于Taylor随机有限元的机械零件模糊可靠度计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003049905A (ja) * | 2001-08-07 | 2003-02-21 | Mitsuboshi Belting Ltd | 有限要素法解析を用いた伝動ベルトの荷重分担予測方法及び装置、並びにプログラム |
CN103177168A (zh) * | 2013-04-24 | 2013-06-26 | 湖南大学 | 一种齿轮齿面摩擦因数的反求方法 |
CN105138858A (zh) * | 2015-09-25 | 2015-12-09 | 南车株洲电力机车研究所有限公司 | 一种基于多体多力学的风力发电机齿轮箱优化设计方法 |
CN105631084A (zh) * | 2015-06-04 | 2016-06-01 | 重庆大学 | 行星齿轮减速齿轮箱箱体轻量化结构方法 |
-
2016
- 2016-10-27 CN CN201610955932.0A patent/CN106528991B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003049905A (ja) * | 2001-08-07 | 2003-02-21 | Mitsuboshi Belting Ltd | 有限要素法解析を用いた伝動ベルトの荷重分担予測方法及び装置、並びにプログラム |
CN103177168A (zh) * | 2013-04-24 | 2013-06-26 | 湖南大学 | 一种齿轮齿面摩擦因数的反求方法 |
CN105631084A (zh) * | 2015-06-04 | 2016-06-01 | 重庆大学 | 行星齿轮减速齿轮箱箱体轻量化结构方法 |
CN105138858A (zh) * | 2015-09-25 | 2015-12-09 | 南车株洲电力机车研究所有限公司 | 一种基于多体多力学的风力发电机齿轮箱优化设计方法 |
Non-Patent Citations (1)
Title |
---|
高速泵齿轮箱的模态分析;杨敏;《化工设备与管道》;20160831;第53卷(第4期);55-58 |
Also Published As
Publication number | Publication date |
---|---|
CN106528991A (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marjanovic et al. | A practical approach to the optimization of gear trains with spur gears | |
Pedersen | Improving bending stress in spur gears using asymmetric gears and shape optimization | |
Shweiki et al. | Transmission Error and strain analysis of lightweight gears by using a hybrid FE-analytical gear contact model | |
Wei | Stresses and deformations in involute spur gears by finite element method | |
CN106407609B (zh) | 基于Neumann随机有限元对齿轮箱进行优化设计的方法 | |
CN106528991B (zh) | 基于Taylor随机有限元对齿轮箱进行优化设计的方法 | |
Krol et al. | 3D-modeling and optimization spindle's node machining centre SVM1F4 | |
Dong et al. | Optimum design of the tooth root profile for improving bending capacity | |
CN111027156B (zh) | 含裂纹齿轮的工业机器人减速器传动精度可靠性分析方法 | |
Chen et al. | Probabilistic design optimization of wind turbine gear transmission system based on dynamic reliability | |
Li et al. | Structural design and optimization of the crossbeam of a computer numerical controlled milling-machine tool using sensitivity theory and NSGA-II algorithm | |
Chavadaki et al. | Finite element analysis of spur gear to find out the optimum root radius | |
CN104102774B (zh) | 一种精密机床滚动直线导轨安装面直线度的设计方法 | |
Zhao et al. | Monomer model: an integrated characterization method of geometrical deviations for assembly accuracy analysis | |
Mao et al. | Numerical calculation method of meshing stiffness for the beveloid gear considering the effect of surface topography | |
İrsel | Effects of modification on the strength–weight ratio of standard bevel gears | |
Chen et al. | Nonlinear mesh stiffness model using slice coupling for straight bevel gear considering axial mesh force component and extended tooth contact | |
Ma et al. | Meshing contact impact properties of circular arc tooth trace cylindrical gear based on rotating knife dish milling process | |
Van Thoan et al. | Optimization design for spur gear with stress-relieving holes | |
Gao et al. | Optimization and experimental research on a new-type short cylindrical cup-shaped harmonic reducer | |
Kursat Celik et al. | Analytical and finite element method based stress analysis of rotary elements: case study for the motion transmission gears of a rotary drum mower | |
CN107256294B (zh) | 一种基于等效刚度法的sss离合器动力学建模方法 | |
Mirandola et al. | Geometry evolution prediction and process settings influence in profiled ring rolling | |
Mote et al. | Design and FEM analysis of helical gear | |
Jani et al. | Design, Modelling and Analysis of Helical Gear Pair using ANSYS and AGMA Standards for Calculating a Bending and Contact Stress on Gear Profiles. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |