CN104298814B - 一种基于参数误差累积的齿轮系统性能可靠度计算方法 - Google Patents
一种基于参数误差累积的齿轮系统性能可靠度计算方法 Download PDFInfo
- Publication number
- CN104298814B CN104298814B CN201410490573.7A CN201410490573A CN104298814B CN 104298814 B CN104298814 B CN 104298814B CN 201410490573 A CN201410490573 A CN 201410490573A CN 104298814 B CN104298814 B CN 104298814B
- Authority
- CN
- China
- Prior art keywords
- parameter
- gear train
- stress
- gear
- performance reliability
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种基于参数误差累积的齿轮系统性能可靠度计算方法,其步骤如下:一:装配参数类型确定;二:敏感参数类型确定;三:敏感参数响应面分析;四:齿轮系统性能可靠度计算。本发明弥补了传统性能可靠度计算方法在对装配累积误差影响考虑不足的缺点,通过对齿轮系统模型建立坐标并进行模型偏移和旋转处理,使模型装配形式接近实际情况,实现了模型装配敏感参数在其误差范围内的连续变换;针对各组设计点的应力仿真由软件在控制了其他仿真条件后自动进行,保证了条件的一致性同时也提高了计算效率;通过对响应面的分析确定齿轮系统不同装配下的安全因子分布,由此来确定性能可靠度,减少了仿真次数,提高了效率。
Description
技术领域
本发明提供一种基于参数误差累积的齿轮系统性能可靠度计算方法,本计算方法主要考虑齿轮系统在多种装配参数误差累计影响下的性能可靠度。属于机械产品可靠度设计分析技术领域。
背景技术
产品的装配偏差在实际加工装配时必然存在,装配偏差是产品性能可靠度的一个重要的影响因素。目前,机械产品的性能可靠度计算,很少考虑由于装配偏差产生的影响。针对产品装配误差累积的性能可靠度计算方法也比较少。在齿轮性能可靠度计算方法的研究,还没有综合考虑齿轮中心距偏差和齿轮轴线平行度误差的性能可靠度计算方法。
发明内容
(1)目的:本发明的目的是提供一种基于参数误差累积的齿轮系统性能可靠度计算方法,基于装配参数在各自公差范围所产生的误差积累,综合考虑装配参数之间的耦合关系,通过多次有限元仿真分析建立响应面及对响应面进行抽样,拟合基于材料屈服极限的安全因子分布,得到多参数综合作用下齿轮系统的性能可靠度。
(2)技术方案:
本发明是通过以下技术方案实现的,首先对齿轮系统进行分析,得出齿轮系统的具体结构尺寸信息,根据公差与配合使用手册,确定齿轮装配时存在偏差的参数类型,并确定参数的公差范围;其次,通过传统的应力分析确定齿轮系统受力薄弱位置,对齿轮系统处于极限公差装配的模型进行应力分析,通过敏感性分析确定敏感参数;针对敏感参数进行抽样仿真,形成安全因子的响应面,并对响应面进行抽样得到安全因子的分布,计算产品的性能可靠度。
本发明是一种基于参数误差累积的齿轮系统性能可靠度计算方法,其步骤如下:
步骤一:装配参数类型确定
对齿轮系统进行结构分析,测量得出齿轮系统的具体结构尺寸信息,根据公差与配合使用手册,计算确定齿轮装配时存在偏差的参数类型,并确定存在偏差参数的公差范围;
步骤二:敏感参数类型确定
对齿轮系统中存在的装配参数类型利用计算机仿真软件进行单偏差参数应力仿真分析,并分别与无偏差状态下的应力仿真分析结果进行比较,确定每种参数类型对齿轮系统影响的敏感程度。并对齿轮系统中存在的装配偏差类型利用计算机仿真软件进行组合双偏差参数应力仿真分析,并分别与无偏差状态下的应力仿真分析结果进行比较,确定参数类型组合对齿轮系统的影响;
综合以上结果将敏感程度高的参数类型作为齿轮系统的性能可靠度仿真的输入参数;
步骤三:敏感参数响应面分析
选择Latin Hypercube Sampling Design(LHS)超拉丁法为设计点抽样方法,设定抽取样本设计点的数量,利用计算机对各设计点进行自动应力仿真,得出齿轮系统模型在各组设计点所对应的装配偏差状态下的应力输出参数(主要为安全因子);设置响应面类型,利用各组设计点的应力仿真结果生成响应面;
步骤四:齿轮系统性能可靠度计算
在所生成的响应面上,利用Latin Hypercube Sampling Design(LHS)超拉丁法进行样本点的抽样,得到样本点应力输出参数的分布,将满足产品设计要求的安全系数百分比,作为齿轮系统的性能可靠度计算结果。
其中,在步骤一中所述的“对齿轮系统进行分析,得出齿轮系统的具体结构尺寸信息”,其作法如下:
1)对齿轮系统的各组成齿轮进行基本尺寸的测量;
2)结合所得的尺寸测量数据,对齿轮系统的装配形式进行分析;
3)综合尺寸数据及装配形式确定齿轮系统的结构尺寸信息。
其中,在步骤二中所述的“计算机仿真软件”,是指ANSYS有限元分析软件。
其中,在步骤二中所述的“敏感参数类型确定”,主要包括四个步骤:
1)针对齿轮系统模型进行坐标系的建立,结合装配参数类型及坐标系对模型进行偏移旋转的设置;
2)施加分析条件,对无偏差齿轮系统模型进行应力仿真分析,提取出啮合区域齿根处的应力仿真结果;
3)对齿轮系统模型进行单个装配参数的修改,保持分析条件不变,对修改后的模型进行应力仿真分析,提取出啮合区域齿根处的应力仿真结果;
4)对齿轮系统装配参数两两组合后对模型进行装配参数进行的修改,保持分析条件不变,对修改后的模型进行应力仿真分析,提取出啮合区域齿根处的应力仿真结果。
其中,在步骤二中所述的“应力仿真分析”是指应用计算机仿真软件(ANSYS)对齿轮系统模型进行静力仿真,得到模型正常工作状态下的应力状态,并对所得结果进行分析。
其中,在步骤三中所述的“应力输出参数”是指通过对齿轮系统模型进行应力仿真模拟后得到的,齿轮啮合区域齿根处最大等效应力以及根据材料属性确定的齿轮安全因子。
其中,在步骤三中所述的“安全因子”指的是材料的屈服极限与计算所得的最大等效应力之比。
其中,在步骤四中所述的“性能可靠度”是指在允许装配误差范围内齿轮系统实际的安全因子在设计安全系数范围内的百分比。
(3)优点及功效:本发明提供一种基于参数误差累积的齿轮系统性能可靠度计算方法,具有以下优点:
1)本发明基于装配参数在各自公差范围所产生的误差累积,综合考虑装配参数之间的耦合关系,弥补了传统性能可靠度计算方法在对装配累积误差影响考虑不足的缺点。
2)本发明通过对齿轮系统模型建立坐标并进行模型偏移和旋转处理,使模型装配形式接近实际情况,实现了模型装配敏感参数在其误差范围内的连续变换;
3)本发明针对参数所在的公差范围及分布情况选择的样本点抽取方法,所得的样本点能反映出参数分布类型的真实情况,各组样本设计点更为合理。
4)本发明针对各组设计点的应力仿真由软件在在控制了其他仿真条件后自动进行,保证了条件的一致性同时也提高了计算效率;
5)本发明通过对响应面的分析确定齿轮系统不同装配下的安全因子分布,由此来确定性能可靠度,减少了仿真次数,提高了效率。
附图说明
图1为本发明所述方法流程框图
具体实施方式
以下将结合附图和实施例对本发明做进一步的详及细说明。
本发明是一种基于参数误差累积的齿轮系统性能可靠度计算方法,见图1所示,该方法具体步骤如下:
步骤一:装配参数类型确定
对齿轮系统进行分析,得出齿轮系统的具体结构尺寸信息,根据公差与配合使用手册,确定齿轮装配时存在偏差的参数类型,并确定参数的公差范围。
根据齿轮对的实物及互换性的分析,一对齿轮的装配偏差参数包括齿轮中心距偏差及齿轮轴线平行度误差。结合齿轮系统结构尺寸和装配尺寸,查找互换性手册后易得出:齿轮对的中心距极限偏差为0.02mm,齿轮轴线平行度为:轴线平行度△Y公差为0.027mm,轴线平行度△X公差为0.054mm。利用反三角函数,将齿轮轴线平行度公差转化为角度。齿轮中心距偏差,齿轮轴线平行度误差分别对应参数fa,fy,fx。具体参数情况如下表1所示:
表1输入参数列表
步骤二:敏感参数类型确定
对齿轮系统中存在的装配参数类型进行单因素仿真分析,并分别与无偏差状态下的应力仿真分析结果进行比较,确定每种参数类型对齿轮系统影响的敏感程度。并对齿轮系统中存在的装配偏差类型进行组合双因素仿真分析,并分别与无偏差状态下的应力仿真分析结果进行比较,确定参数类型组合对齿轮系统的影响。综合以上结果将敏感程度高的参数类型作为齿轮系统的性能可靠度仿真的输入参数。
利用步骤一种中的结果,以齿轮A的中心为坐标原点,以由齿轮A中心指向齿轮B的中心为H方向,以齿轮轴的轴线方向为Z方向,第三方向为V。以建立的坐标轴为基础,将沿H轴负方向偏移距离设为fa,单位mm;将绕H的旋转角度设为fy,单位为°;将绕V的旋转角度设为fx,单位为°。以fa,fy,fx均为0时的模型为齿轮的标准(无偏差)装配模型。
首先对无偏差装配模型利用Ansys进行静力仿真,得出应力分析结果。在仿真过程中,将两齿轮啮合区域设置为摩擦接触,对齿轮分别添加旋转副;将啮合部分齿面的网格进行细化;在齿轮A上添加700N·M,在齿轮B上施加固定约束,将齿根处屈服极限设为1200MPa用以计算安全因子。
分别对fa,fy,fx三个参数进行单独修改,实现模型以不同的单个装配偏差的装配,保持上述应力分析条件不变,分别对修改后模型进行静力仿真,针对每个参数变化,共得出六组应力参数结果。将六组参数结果依次与无偏差装配模型应力参数结果进行比较。经比较得出模型在三个参数变换后,静力仿真后所得应力分析结果均有较大的改变,因此确定此三个参数为影响齿轮性能可靠度的敏感参数。
将三个参数fa,fy,fx两两组合并进行修改,实现模型以不同的组合装配偏差的装配,保持上述应力分析条件不变,分别对修改后模型进行静力仿真,针对每组组合参数的变化,共得出六组应力参数结果。将六组参数结果依次与无偏差装配模型应力分析结果进行比较。经比较得出模型在三个参数进行两两组合变换后,静力仿真后所得应力参数结果仍有较大的改变,因此确定此三个参数对齿轮对性能可靠度的影响并未彼此抵消。具体结果如下表2所示:
表2仿真应力汇总
步骤三:敏感参数响应面分析
选择Latin Hypercube Sampling Design(LHS)超拉丁法为设计点抽样方法,设定抽取样本设计点的数量,利用计算机对各设计点进行自动应力仿真,得出齿轮系统模型在各组设计点所对应的装配偏差状态下的应力输出参数(主要为安全因子)。设置响应面类型,利用各组设计点的应力仿真结果生成响应面。
在参数的分布范围内利用超拉丁法进行30次随机样本的抽取,,然后利用Ansys进行样本自动应力仿真,得出30组样本的啮合区域齿根处和齿面处最大等效应力及对应安全因子。利用样本的应力仿真结果,将响应面类型设置为Kriging,进行响应面生成。
步骤四:性能可靠度计算
选择Latin Hypercube Sampling Design(LHS)超拉丁法为设计点抽样方法,设定抽取样本设计点的数量,利用计算机对各设计点进行自动应力仿真,得出模型在各组设计点所对应的装配偏差状态下的应力输出参数。设置响应面类型,利用各组设计点的应力仿真结果生成响应面。
利用超拉丁法,针对生成的响应面进行10000次样本点的抽取。针对样本点的抽取结果,对样本点的应力分布进行分析。在齿根处的样本最大等效应力分布中输入应力极限1200MPa,可得样本的应力均在应力极限之内;在安全因子的分布中输入安全系数1.5,在安全系数范围内的百分比(性能可靠度)为96.1%。
Claims (6)
1.一种基于参数误差累积的齿轮系统性能可靠度计算方法,其特于征在:其步骤如下:
步骤一:装配参数类型确定
对齿轮系统进行结构分析,测量得出齿轮系统的具体结构尺寸信息,根据公差与配合使用手册,计算确定齿轮装配时存在偏差的参数类型,并确定存在偏差参数的公差范围;
步骤二:敏感参数类型确定
对齿轮系统中存在的装配偏差参数类型利用计算机仿真软件进行单偏差参数应力仿真分析,并分别与无偏差状态下的应力仿真分析结果进行比较,确定每种参数类型对齿轮系统影响的敏感程度;并对齿轮系统中存在的装配偏差参数类型利用计算机仿真软件进行组合双偏差参数应力仿真分析,并分别与无偏差状态下的应力仿真分析结果进行比较,确定参数类型组合对齿轮系统的影响;
综合以上结果将敏感程度高的参数类型作为齿轮系统的性能可靠度仿真的输入参数;
步骤三:敏感参数响应面分析
选择Latin Hypercube Sampling Design即LHS超拉丁法为设计点抽样方法,设定抽取样本设计点的数量,利用计算机对各设计点进行自动应力仿真,得出齿轮系统模型在各组设计点所对应的装配偏差状态下的应力输出参数;设置响应面类型,利用各组设计点的应力仿真结果生成响应面;
步骤四:齿轮系统性能可靠度计算
在所生成的响应面上,利用Latin Hypercube Sampling Design即LHS超拉丁法进行样本点的抽样,得到样本点应力输出参数的分布,将满足产品设计要求的安全系数百分比,作为齿轮系统的性能可靠度计算结果;
其中,在步骤三中所述的应力输出参数,是指通过对齿轮系统模型进行应力仿真模拟后得到的,齿轮啮合区域齿根处最大等效应力以及根据材料属性确定的齿轮安全因子;
其中,在步骤四中所述的性能可靠度,是指在允许装配误差范围内齿轮系统实际的安全因子在设计安全系数范围内的百分比。
2.根据权利要求1所述的一种基于参数误差累积的齿轮系统性能可靠度计算方法,其特于征在:在步骤一中所述的对齿轮系统进行结构分析,测量得出齿轮系统的具体结构尺寸信息,其作法如下:
1)对齿轮系统的各组成齿轮进行基本尺寸的测量;
2)结合所得的尺寸测量数据,对齿轮系统的装配形式进行分析;
3)综合尺寸数据及装配形式确定齿轮系统的结构尺寸信息。
3.根据权利要求1所述的一种基于参数误差累积的齿轮系统性能可靠度计算方法,其特于征在:在步骤二中所述的计算机仿真软件,是指ANSYS有限元分析软件。
4.根据权利要求1所述的一种基于参数误差累积的齿轮系统性能可靠度计算方法,其特于征在:在步骤二中所述的敏感参数类型确定,包括四个步骤:
1)针对齿轮系统模型进行坐标系的建立,结合装配参数类型及坐标系对模型进行偏移旋转的设置;
2)施加分析条件,对无偏差齿轮系统模型进行应力仿真分析,提取出啮合区域齿根处的应力仿真结果;
3)对齿轮系统模型进行单个装配参数的修改,保持分析条件不变,对修改后的模型进行应力仿真分析,提取出啮合区域齿根处的应力仿真结果;
4)对齿轮系统装配参数两两组合后对模型进行两个装配参数的修改,保持分析条件不变,对修改后的模型进行应力仿真分析,提取出啮合区域齿根处的应力仿真结果。
5.根据权利要求1所述的一种基于参数误差累积的齿轮系统性能可靠度计算方法,其特于征在:在步骤二中所述的应力仿真分析,是指应用计算机仿真软件ANSYS对齿轮系统模型进行静力仿真,得到模型正常工作状态下的应力状态,并对所得结果进行分析。
6.根据权利要求1所述的一种基于参数误差累积的齿轮系统性能可靠度计算方法,其特于征在:在步骤三中所述的安全因子,指的是材料的屈服极限与计算所得的最大等效应力之比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410490573.7A CN104298814B (zh) | 2014-09-23 | 2014-09-23 | 一种基于参数误差累积的齿轮系统性能可靠度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410490573.7A CN104298814B (zh) | 2014-09-23 | 2014-09-23 | 一种基于参数误差累积的齿轮系统性能可靠度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104298814A CN104298814A (zh) | 2015-01-21 |
CN104298814B true CN104298814B (zh) | 2017-05-17 |
Family
ID=52318538
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410490573.7A Active CN104298814B (zh) | 2014-09-23 | 2014-09-23 | 一种基于参数误差累积的齿轮系统性能可靠度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104298814B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105404738B (zh) * | 2015-11-17 | 2018-10-02 | 天津百利机械装备研究院有限公司 | 一种用于指导齿轮应力检测的齿根应力分析方法 |
CN106599448B (zh) * | 2016-12-12 | 2019-07-05 | 北京航空航天大学 | 一种基于动态可靠度的齿轮系统公差优化计算方法 |
CN109409028B (zh) * | 2019-01-09 | 2022-04-19 | 北京航空航天大学 | 基于确信可靠度的齿轮可靠性分析方法 |
CN110704940B (zh) * | 2019-08-30 | 2023-08-18 | 中国航空工业集团公司沈阳飞机设计研究所 | 飞机发动机安装交点位置偏差处理方法及发动机推力销 |
CN113051677A (zh) * | 2021-03-16 | 2021-06-29 | 南京航空航天大学 | 一种基于ansys-workbench的转子系统圆弧端齿连接结构参数化建模方法 |
CN113239491B (zh) * | 2021-05-17 | 2022-04-01 | 湖南南方宇航高精传动有限公司 | 一种风电齿轮箱中箱体加强筋的多参数优化设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393864A (zh) * | 2011-06-28 | 2012-03-28 | 电子科技大学 | 一种基于故障物理的航天器用谐波齿轮的可靠性优化方法 |
CN104008240A (zh) * | 2014-05-22 | 2014-08-27 | 北京航空航天大学 | 一种在轨空间柔性齿轮机构动态耦合时变故障率分析方法 |
-
2014
- 2014-09-23 CN CN201410490573.7A patent/CN104298814B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393864A (zh) * | 2011-06-28 | 2012-03-28 | 电子科技大学 | 一种基于故障物理的航天器用谐波齿轮的可靠性优化方法 |
CN104008240A (zh) * | 2014-05-22 | 2014-08-27 | 北京航空航天大学 | 一种在轨空间柔性齿轮机构动态耦合时变故障率分析方法 |
Non-Patent Citations (4)
Title |
---|
"基于响应面法齿轮啮合传动可靠性灵敏度分析";李昌 等;《航空动力学报》;20110331;第26卷(第3期);第711-715页 * |
"考虑安装与制造误差的齿轮动态接触仿真";佟操 等;《东北大学学报(自然科学版)》;20140731;第35卷(第7期);第996-1000页 * |
"考虑装配误差的轴系结构静力学对比分析";王宇宁 等;《中国工程机械学报》;20140228;第12卷(第1期);第13-17、47页 * |
"蒙特卡洛模拟及其在公差设计中的应用";袁贵星 等;《天津科技大学学报》;20080630;第23卷(第2期);第60-64页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104298814A (zh) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104298814B (zh) | 一种基于参数误差累积的齿轮系统性能可靠度计算方法 | |
Xie et al. | Variable selection methods for probabilistic load forecasting: Empirical evidence from seven states of the united states | |
CN103267507B (zh) | 基于有限元分析提取机械结构平面的平面度误差的方法 | |
CN104915483B (zh) | 一种地震作用下建筑物基础稳定性验算方法 | |
CN109918833B (zh) | 一种数值模拟可信度的定量分析方法 | |
CN106055827A (zh) | 一种油藏数值模拟参数敏感性分析装置及方法 | |
丁力 et al. | Design and experiment of seed plate of corn air suction seed metering device based on DEM-CFD | |
CN105868496A (zh) | 一种面向装配的矩形平面形状误差评价参数确定方法 | |
CN103366056A (zh) | 基于有限元分析的机械零件圆柱度误差的测量方法 | |
CN105631135B (zh) | 基于abaqus基坑开挖快速建模分析系统及其应用方法 | |
CN104036085B (zh) | 用于cae仿真分析的复杂型面汽车声学部件的厚度统计方法 | |
CN108709511A (zh) | 一种汽轮机全实缸状态下通流间隙检测方法 | |
CN104036118B (zh) | 一种电力系统并行化轨迹灵敏度获取方法 | |
CN102495964A (zh) | 一种二维多工位装配成功率计算方法 | |
CN103366095B (zh) | 一种基于坐标变换的最小二乘拟合信号处理方法 | |
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 | |
CN109388868A (zh) | 一种考虑钉孔接触及钉弯曲的钉载计算方法 | |
CN100552722C (zh) | 三角网格上多尺度曲率特征提取方法 | |
Kosnikov | Increasing the information capacity of the interface of the control system for multyparameter objects | |
Grzesik | A variational principle for the Milne problem linear extrapolation length | |
CN110717244A (zh) | 基于平均偏离度算法的数据信任度分析计算机仿真方法 | |
Bi et al. | Reliability Analysis of Cantilever Beam based on Hybrid Stochastic Finite Element Method | |
CN104771186B (zh) | 一种用于硬膜下血肿体积计算的应用系统 | |
CN104636557A (zh) | 计算机辅助工程cae模拟仿真曲轴模锻工艺的方法 | |
Yuan | Research on Manufacturing Process Planning Technology of Aero-Engine Blades Based on Particle Swarm Optimization |
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 |