CN110263492A - 一种谐波减速器双圆弧齿廓扭转刚度计算方法 - Google Patents

一种谐波减速器双圆弧齿廓扭转刚度计算方法 Download PDF

Info

Publication number
CN110263492A
CN110263492A CN201910633304.4A CN201910633304A CN110263492A CN 110263492 A CN110263492 A CN 110263492A CN 201910633304 A CN201910633304 A CN 201910633304A CN 110263492 A CN110263492 A CN 110263492A
Authority
CN
China
Prior art keywords
indicate
gear
stiffness
speed reducer
circular arc
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
Application number
CN201910633304.4A
Other languages
English (en)
Other versions
CN110263492B (zh
Inventor
刘志峰
张涛
杨聪彬
张彩霞
胡秋实
纪新宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201910633304.4A priority Critical patent/CN110263492B/zh
Publication of CN110263492A publication Critical patent/CN110263492A/zh
Application granted granted Critical
Publication of CN110263492B publication Critical patent/CN110263492B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

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)
  • Retarders (AREA)
  • Gears, Cams (AREA)

Abstract

本发明公开了一种谐波减速器双圆弧齿廓扭转刚度计算方法,步骤一:计算谐波减速器双圆弧齿廓单齿啮合刚度此时双圆弧接触齿对被视为悬臂梁模型;步骤二:计算谐波减速器双圆弧齿廓多齿啮合刚度;步骤三:计算谐波减速器双圆弧齿廓整机扭转刚度K;本发明通过分析谐波减速器刚柔轮共轭包络理论,通通过分析谐波减速器扭转刚度的组成,研究双圆弧齿廓单齿啮合刚度、啮合齿数以及其他结构件的扭转刚度串联模型,为谐波减速器扭转刚度和传动精度的提高奠定了理论基础。

Description

一种谐波减速器双圆弧齿廓扭转刚度计算方法
技术领域
本发明涉及谐波减速器的设计与制造领域技术领域,特别是涉及一种谐波减速器双圆弧齿廓扭转刚度计算模型。
背景技术
谐波减速器是机器人关节的核心元件,谐波减速器柔轮和刚轮的啮合属于大变形条件下小模数多齿啮合。谐波减速器的扭转刚度的定义为在输入端固定的条件下,输出端负载扭矩和转角之间的关系曲线。目前谐波减速器双圆弧齿廓的整机扭转刚度模型多采用实验和有限元仿真的方式获得,实验的方式对于设备和工装的要求较多,每个谐波减速器都进行实验研究过于繁琐,而将谐波减速器整机放入有限元模型中进行计算,有限元网格过度,对于服务器的性能要求较高。因此,本文采用一种基于有限元仿真和数值计算相结合的方式,提出一种谐波减速器双圆弧齿廓扭转刚度计算模型。
发明内容
本发明目的是:为提高谐波减速器传动精度,通过分析谐波减速器扭转刚度的组成,研究双圆弧齿廓单齿啮合刚度、啮合齿数以及其他结构件的扭转刚度串联模型,提出一种谐波减速器双圆弧齿廓扭转刚度计算模型。
本发明所采取的技术方案是:
一种谐波减速器双圆弧齿廓扭转刚度计算方法,该方法包括如下步骤:
步骤一:计算谐波减速器双圆弧齿廓单齿啮合刚度此时双圆弧接触齿对被视为悬臂梁模型:
式中,角标1和2分别表示柔轮和刚轮,Kb表示轮齿的弯曲刚度,Ks表示轮齿的剪切刚度,Ka表示轮齿的轴向压缩刚度,Kf表示轮齿的基础形变刚度,Kh表示柔轮和刚轮齿廓的接触刚度。
步骤1.1:计算柔轮和刚轮的弯曲刚度
步骤1.2:计算剪切刚度
步骤1.3:计算轴向压缩刚度
步骤1.4:计算基础形变刚度
式中,d表示啮合点与齿根圆的水平距离,E表示材料的弹性模量G表示材料的剪切模量,h表示啮合点与齿廓中心轴线的距离,x表示啮合点与积分位置的水平距离,Ix表示积分位置的截面惯性矩,Ax表示积分位置的截面面积,α表示啮合点的压力角,W表示齿宽,Sf表示齿根圆半径,uf表示啮合点延长线与齿廓中轴线的交点与齿根圆的水平距离。A*、B*、C*和D*表示的是基础形变系数,可以通过如下矩阵计算:
式中,θf表示的是单齿廓中心角的一半,hfi表示的是壁厚系数;
式中,ρ表示柔轮中性层变形之后的矢径,在不同波发生器转角位置的齿根圆半径是不同的,δ表示柔轮壁厚,rw表示刚轮外圆半径,rf2表示刚轮齿根圆半径。
步骤1.5:计算柔轮和刚轮的接触刚度
刚轮和柔轮的接触属于粗糙表面的接触,接触刚度Kh通过提取啮合点的接触压力和接触面积采用分形理论获得。
步骤1.6:计算柔轮和刚轮的啮合点的压力角α
谐波减速器双圆弧齿廓在设计过程中啮合压力角处于实时变化的过程,因此需要计算任意啮合点的瞬时压力角,压力角表示的是接触压力和瞬时速度的夹角。
步骤二:计算谐波减速器双圆弧齿廓多齿啮合刚度:
步骤2.1:计算瞬时接触齿数z及对应的波发生器转角
步骤2.2:计算波发生器转角条件下的瞬时作用力臂di和单齿啮合刚度Kdi
步骤2.3:计算刚柔轮齿对的等效啮合刚度:
步骤三:计算谐波减速器双圆弧齿廓整机扭转刚度K:
K=1/(1/Kin+1/Kout+1/Kb+1/KH+1/Kc) (1-9)
式中,Kin和Kout分别表示输入和输出轴的扭转刚度,Kb表示的波发生器的扭转刚度,KH表示柔轮筒臂的扭转刚度,Kc表示刚柔轮齿的啮合刚度。
本发明具有的优点和积极效果是:
本发明通过分析谐波减速器刚柔轮共轭包络理论,通通过分析谐波减速器扭转刚度的组成,研究双圆弧齿廓单齿啮合刚度、啮合齿数以及其他结构件的扭转刚度串联模型,为谐波减速器扭转刚度和传动精度的提高奠定了理论基础。
附图说明
图1谐波减速器双圆弧齿廓单齿刚度模型图;
图2谐波减速器单齿基础形变刚度示意图;
图3谐波减速器多齿啮合力和力臂示意图;
图4谐波减速器整机结构及刚度组成示意图;
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
一种谐波减速器双圆弧齿廓扭转刚度计算模型,该模型包括如下步骤:
步骤一:计算谐波减速器双圆弧齿廓单齿啮合刚度此时双圆弧接触齿对被视为悬臂梁模型:
式中,角标1和2分别表示柔轮和刚轮,Kb表示轮齿的弯曲刚度,Ks表示轮齿的剪切刚度,Ka表示轮齿的轴向压缩刚度,Kf表示轮齿的基础形变刚度,Kh表示柔轮和刚轮齿廓的接触刚度。
步骤1.1:计算柔轮和刚轮的弯曲刚度
步骤1.2:计算剪切刚度
步骤1.3:计算轴向压缩刚度
步骤1.4:计算基础形变刚度
式中,如图1所示,d表示啮合点与齿根圆的水平距离,E表示材料的弹性模量G表示材料的剪切模量,h表示啮合点与齿廓中心轴线的距离,x表示啮合点与积分位置的水平距离,Ix表示积分位置的截面惯性矩,Ax表示积分位置的截面面积,α表示啮合点的压力角,W表示齿宽,Sf表示齿根圆半径,uf表示啮合点延长线与齿廓中轴线的交点与齿根圆的水平距离。A*、B*、C*和D*表示的是基础形变系数,可以通过如下矩阵计算:
式中,如图2所示,θf表示的是单齿廓中心角的一半,hfi表示的是壁厚系数;
式中,ρ表示柔轮中性层变形之后的矢径,在不同波发生器转角位置的齿根圆半径是不同的,δ表示柔轮壁厚,rw表示刚轮外圆半径,rf2表示刚轮齿根圆半径。
步骤1.5:计算柔轮和刚轮的接触刚度
刚轮和柔轮的接触属于粗糙表面的接触,接触刚度Kh通过提取啮合点的接触压力和接触面积采用分形理论获得。
步骤1.6:计算柔轮和刚轮的啮合点的压力角α
谐波减速器双圆弧齿廓在设计过程中啮合压力角处于实时变化的过程,因此需要计算任意啮合点的瞬时压力角,压力角表示的是接触压力和瞬时速度的夹角。
步骤二:计算谐波减速器双圆弧齿廓多齿啮合刚度:
步骤2.1:计算瞬时接触齿数z及对应的波发生器转角
步骤2.2:计算波发生器转角条件下的瞬时作用力臂di和单齿啮合刚度Kdi
步骤2.3:如图3所示,计算刚柔轮齿对的等效啮合刚度:
步骤三:如图4所示,计算谐波减速器双圆弧齿廓整机扭转刚度K:
K=1/(1/Kin+1/Kout+1/Kb+1/KH+1/Kc) (1-9)
式中,Kin和Kout分别表示输入和输出轴的扭转刚度,Kb表示的波发生器的扭转刚度,KH表示柔轮筒臂的扭转刚度,Kc表示刚柔轮齿的啮合刚度。
本发明具有的优点和积极效果是:
本发明通过分析谐波减速器刚柔轮共轭包络理论,通通过分析谐波减速器扭转刚度的组成,研究双圆弧齿廓单齿啮合刚度、啮合齿数以及其他结构件的扭转刚度串联模型,为谐波减速器扭转刚度和传动精度的提高奠定了理论基础。

Claims (1)

1.一种谐波减速器双圆弧齿廓扭转刚度计算方法,其特征在于:该方法包括如下步骤,
步骤一:计算谐波减速器双圆弧齿廓单齿啮合刚度此时双圆弧接触齿对被视为悬臂梁模型:
式中,角标1和2分别表示柔轮和刚轮,Kb表示轮齿的弯曲刚度,Ks表示轮齿的剪切刚度,Ka表示轮齿的轴向压缩刚度,Kf表示轮齿的基础形变刚度,Kh表示柔轮和刚轮齿廓的接触刚度;
步骤1.1:计算柔轮和刚轮的弯曲刚度
步骤1.2:计算剪切刚度
步骤1.3:计算轴向压缩刚度
步骤1.4:计算基础形变刚度
式中,d表示啮合点与齿根圆的水平距离,E表示材料的弹性模量G表示材料的剪切模量,h表示啮合点与齿廓中心轴线的距离,x表示啮合点与积分位置的水平距离,Ix表示积分位置的截面惯性矩,Ax表示积分位置的截面面积,α表示啮合点的压力角,W表示齿宽,Sf表示齿根圆半径,uf表示啮合点延长线与齿廓中轴线的交点与齿根圆的水平距离;A*、B*、C*和D*表示的是基础形变系数,可以通过如下矩阵计算:
式中,θf表示的是单齿廓中心角的一半,hfi表示的是壁厚系数;
式中,ρ表示柔轮中性层变形之后的矢径,在不同波发生器转角位置的齿根圆半径是不同的,δ表示柔轮壁厚,rw表示刚轮外圆半径,rf2表示刚轮齿根圆半径;
步骤1.5:计算柔轮和刚轮的接触刚度
刚轮和柔轮的接触属于粗糙表面的接触,接触刚度Kh通过提取啮合点的接触压力和接触面积采用分形理论获得;
步骤1.6:计算柔轮和刚轮的啮合点的压力角α
谐波减速器双圆弧齿廓在设计过程中啮合压力角处于实时变化的过程,因此需要计算任意啮合点的瞬时压力角,压力角表示的是接触压力和瞬时速度的夹角;
步骤二:计算谐波减速器双圆弧齿廓多齿啮合刚度:
步骤2.1:计算瞬时接触齿数z及对应的波发生器转角
步骤2.2:计算波发生器转角条件下的瞬时作用力臂di和单齿啮合刚度Kdi
步骤2.3:计算刚柔轮齿对的等效啮合刚度:
步骤三:计算谐波减速器双圆弧齿廓整机扭转刚度K:
K=1/(1/Kin+1/Kout+1/Kb+1/KH+1/Kc) (1-9)
式中,Kin和Kout分别表示输入和输出轴的扭转刚度,Kb表示的波发生器的扭转刚度,KH表示柔轮筒臂的扭转刚度,Kc表示刚柔轮齿的啮合刚度。
CN201910633304.4A 2019-07-15 2019-07-15 一种谐波减速器双圆弧齿廓扭转刚度计算方法 Active CN110263492B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910633304.4A CN110263492B (zh) 2019-07-15 2019-07-15 一种谐波减速器双圆弧齿廓扭转刚度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910633304.4A CN110263492B (zh) 2019-07-15 2019-07-15 一种谐波减速器双圆弧齿廓扭转刚度计算方法

Publications (2)

Publication Number Publication Date
CN110263492A true CN110263492A (zh) 2019-09-20
CN110263492B CN110263492B (zh) 2023-05-02

Family

ID=67926055

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910633304.4A Active CN110263492B (zh) 2019-07-15 2019-07-15 一种谐波减速器双圆弧齿廓扭转刚度计算方法

Country Status (1)

Country Link
CN (1) CN110263492B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113722891A (zh) * 2021-08-10 2021-11-30 北京工业大学 一种谐波减速器传动系统扭转刚度拟合方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104074948A (zh) * 2014-07-02 2014-10-01 天津工业大学 具有公切线型双圆弧齿廓的杯形谐波齿轮及其齿廓设计方法
CN105299151A (zh) * 2014-06-11 2016-02-03 南京理工大学 一种谐波齿轮减速器
CN108533715A (zh) * 2018-06-28 2018-09-14 西安交通大学 一种用于谐波齿轮传动的双向共轭齿形设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105299151A (zh) * 2014-06-11 2016-02-03 南京理工大学 一种谐波齿轮减速器
CN104074948A (zh) * 2014-07-02 2014-10-01 天津工业大学 具有公切线型双圆弧齿廓的杯形谐波齿轮及其齿廓设计方法
CN108533715A (zh) * 2018-06-28 2018-09-14 西安交通大学 一种用于谐波齿轮传动的双向共轭齿形设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIU ZHIFENG: "Experimental Studies on Torsional Stiffness of Cycloid Gear Based on Machining Parameters of Tooth Surfaces" *
辛洪兵;: "双圆弧谐波齿轮传动基本齿廓设计" *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113722891A (zh) * 2021-08-10 2021-11-30 北京工业大学 一种谐波减速器传动系统扭转刚度拟合方法

Also Published As

Publication number Publication date
CN110263492B (zh) 2023-05-02

Similar Documents

Publication Publication Date Title
Zheng et al. Design and manufacture of new type of non-circular cylindrical gear generated by face-milling method
CN104820756B (zh) 一种考虑延长啮合的裂纹齿轮转子系统动力参数确定方法
CN109630652B (zh) 一种三圆弧谐波齿轮插齿刀及其齿廓设计方法
CN104408220B (zh) 一种改进的轮齿加载接触分析方法
CN105299151B (zh) 一种谐波齿轮减速器
CN110334460A (zh) 圆柱齿轮啮合刚度计算方法
CN109271715B (zh) 一种基于梯度下降法的谐波减速器柔轮结构集成优化方法
CN107944174B (zh) 一种圆柱齿轮齿向载荷分布系数获取方法
CN108416120B (zh) 一种直齿圆柱齿轮双齿啮合区载荷分配率的确定方法
CN104573196A (zh) 一种斜齿圆柱齿轮时变啮合刚度解析计算方法
CN110245417B (zh) 一种谐波减速器双圆弧齿形啮合点法线斜率计算方法
CN104462637A (zh) 齿轮齿条式变比转向器齿轮副齿条齿廓的设计方法
CN110263492A (zh) 一种谐波减速器双圆弧齿廓扭转刚度计算方法
CN106055825A (zh) 一种基于cad的建模实现方法及系统
CN104455211A (zh) 高阶变性傅里叶非圆齿轮副的设计方法
CN110516384A (zh) 一种多输入多输出齿轮传动系统动力学建模方法
ITBO20090714A1 (it) Ruota dentata con profilo atto ad ingranare con semi-incapsulamento in un'apparecchiatura idraulica ad ingranaggi
Zeng et al. Accurate modelling and transient meshing analysis of involute spur gear based on the principle of gear shaping
Chen et al. The investigation of elongation of the neutral line in harmonic drive
CN102615653A (zh) 双圆弧谐波机器人关节
CN113486476B (zh) 一种磨削双圆弧谐波减速器刚轮插刀的砂轮齿廓设计方法
CN113868755A (zh) 基于基节误差的斜齿轮啮合刚度计算方法
Liang et al. Calculation method of sliding ratios for conjugate-curve gear pair and its application
Mote et al. Design and FEM analysis of helical gear
Mu et al. Dynamic Contact Analysis for RV Reducer

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