WO2018086160A1 - 基于粗糙表面的直齿轮三维接触刚度计算方法 - Google Patents

基于粗糙表面的直齿轮三维接触刚度计算方法 Download PDF

Info

Publication number
WO2018086160A1
WO2018086160A1 PCT/CN2016/107233 CN2016107233W WO2018086160A1 WO 2018086160 A1 WO2018086160 A1 WO 2018086160A1 CN 2016107233 W CN2016107233 W CN 2016107233W WO 2018086160 A1 WO2018086160 A1 WO 2018086160A1
Authority
WO
WIPO (PCT)
Prior art keywords
contact
stiffness
gear
micro
convex
Prior art date
Application number
PCT/CN2016/107233
Other languages
English (en)
French (fr)
Inventor
刘志峰
张涛
赵永胜
毕舒心
Original Assignee
北京工业大学
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 北京工业大学 filed Critical 北京工业大学
Priority to US15/798,320 priority Critical patent/US10088386B2/en
Publication of WO2018086160A1 publication Critical patent/WO2018086160A1/zh

Links

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

Definitions

  • the invention belongs to the field of gear mechanics analysis and dynamics research, and particularly relates to a method for calculating a three-dimensional contact stiffness of a spur gear based on a rough surface.
  • the internal excitation of the gear pair is generated inside the system during the gear pair meshing process, and is also the focus and difficulty in the research of gear system dynamics.
  • the internal excitation of the gear system is mainly divided into three types: time-varying stiffness excitation, time-varying error excitation and meshing impact excitation.
  • the time-varying stiffness excitation is mainly composed of three parts: the bending stiffness, shear stiffness and contact stiffness of the gear teeth.
  • the current application of the IOS6336 meshing stiffness calculation formula and the stiffness excitation of the gear are equivalent to Fourier expansion.
  • the calculation of contact stiffness is based on the Hertz contact theory of smooth surfaces, so it is not possible to accurately simulate the deformation of the contact state of the actual tooth surface. Therefore, in order to establish a more effective gear dynamics model, it is necessary to study the calculation method of three-dimensional contact stiffness of spur gear based on rough surface.
  • the object of the present invention is to provide a calculation method for three-dimensional contact stiffness of a spur gear based on a rough surface, which can effectively and accurately determine the accuracy of the gear contact stiffness model, and fully considers factors such as assembly error and manufacturing error.
  • the problem of uneven distribution of the tooth engagement load is to provide a calculation method for three-dimensional contact stiffness of a spur gear based on a rough surface, which can effectively and accurately determine the accuracy of the gear contact stiffness model, and fully considers factors such as assembly error and manufacturing error.
  • S1 uses a hexahedral mesh, requiring each square mesh area A 0 on the tooth surface to be equal;
  • S2 calculates the smooth tooth surface contact pressure distribution based on the finite element method, and extracts the contact zone node pressure P i ;
  • the S1 specifically includes the following steps:
  • the step 2) specifically includes the following steps:
  • the step 3) specifically includes the following steps:
  • the mesh used in the finite element analysis is already small enough to treat the extracted contact zone node pressure P i as the contact pressure in the square mesh on the tooth surface, so the contact pressure of the different square mesh on the tooth surface is expressed. for:
  • n is the frequency index
  • z is the height of the rough surface profile
  • D is the fractal dimension
  • ⁇ is the size parameter of the spectral density
  • G is the roughness parameter.
  • x represents the sampling length coordinate. The same material and processing method are used to obtain the gear sample block, and the coordinates of the tooth surface topography are obtained on the three-dimensional shape measuring instrument. Then the power spectral density method is used to fit the WM function to obtain two parameters of D and G.
  • the contact of the two rough surfaces in the contact model of the gear pair is simplified to the contact of a rigid surface with the rough surface.
  • the dotted line indicates that the neutral surface indicates the contact position of the ideal smooth gear pair, and the portion where the rough surface above the neutral surface protrudes is defined as a microprotrusion.
  • the rough surface is composed of a myriad of circular microprotrusions of various sizes, assuming that each microprotrusion is separated from each other and the interaction is negligible.
  • a' is the cross-sectional area where the micro-convex body intersects the rigid surface before denaturation
  • a the contact area of the micro-convex body after deformation with the rigid surface is defined as the true contact area.
  • the deformation of the micro-convex is divided into elastic and plastic deformation, and the critical cross-sectional area that determines the elastic or plastic deformation of the micro-convex is determined by the material properties of the gear pair:
  • the lower corners e and p represent elastic and plastic deformation, respectively.
  • H represents the hardness of the softer material
  • H 2.8Y
  • R and ⁇ represent the curvature and normal deformation of the single microprotrusion, respectively.
  • the total elastic deformation force F E and the total plastic deformation force F P of a single cube mesh are expressed as:
  • the step 3) further includes the following steps:
  • the deformation is divided into two stages: full elastic deformation and complete plastic deformation, wherein the plastic deformation stage has a stiffness of zero.
  • the single microprotrusion normal contact stiffness where complete elastic deformation occurs can be expressed as
  • the normal contact stiffness in a single square mesh is obtained by integrating in the fully elastic region.
  • ⁇ denotes the coefficient of static friction
  • ⁇ 1 , ⁇ 2 , G 1 , G 2 represent the Poisson's ratio and shear modulus of the two materials, respectively
  • t and f represent the tangential and normal loads of a single micro-convex, respectively
  • r represents the true contact radius of a single micro-convex
  • a indicates the true contact area.
  • T ⁇ b A r
  • ⁇ b the shear strength of the softer material
  • a r the actual contact area of the square grid
  • the step 4) specifically includes the following steps:
  • the tangential contact stiffness K 2 of the contact tooth surface is
  • ⁇ represents the pressure angle of the gear.
  • the three-dimensional contact stiffness calculation method of the spur gear based on the rough surface proposed by the invention further corrects the gear contact stiffness calculation model, solves the disadvantages of the smooth contact according to the Hertz theory in the original calculation process, and reveals the roughness parameter during the rough tooth surface contact process.
  • the influence law and calculation method of the stiffness characteristics of the gear are based on the finite element three-dimensional contact stiffness calculation model, which is more accurate than the original two-dimensional calculation model or calculation formula, and takes into account the influence of errors in the actual assembly and manufacturing process of the gear, so that the gear contact stiffness is more accurate, The basis of gear dynamics analysis.
  • the established spur gear contact stiffness model contains roughness-related parameters for the first time, and indirectly describes the relationship between machining mode and gear contact stiffness.
  • the fractal theory is applied to the different positive direction grids of the contact area for the first time.
  • the three-dimensional contact stiffness of the gears is obtained, and the gear pair assembly is avoided.
  • the influence of the eccentric load caused by the manufacturing error on the meshing stiffness is avoided.
  • S1, S2, and S3 are equally applicable to gears other than spur gears.
  • Figure 1 is a simplified contact model diagram of a gear pair.
  • Figure 2 is a diagram of a single microprotrusion contact area.
  • Figure 3 is a schematic diagram of the spur gear contact stiffness transformation.
  • a method for calculating a three-dimensional contact stiffness of a spur gear based on a rough surface comprising the following steps:
  • the step 1) specifically includes the following steps:
  • the step 2) specifically includes the following steps:
  • the step 3) specifically includes the following steps:
  • the mesh used in the finite element analysis is small enough to consider the extracted contact zone node pressure P i as the contact pressure in the square mesh on the tooth surface, so the contact pressure of the different square mesh on the tooth surface It can be expressed as:
  • the outline of the rough flank can be represented by the W-M function as:
  • n is the frequency index
  • z is the height of the rough surface profile
  • D is the fractal dimension
  • ⁇ is the size parameter of the spectral density
  • G is the roughness parameter
  • x Indicates the sampling length coordinate.
  • the contact of the two rough surfaces in the contact model of the gear pair is simplified to the contact of a rigid surface with the rough surface, as shown in FIG.
  • the dotted line is the neutral surface indicating the contact position of the ideal smooth gear pair, and the portion where the rough surface above the neutral surface is convex is defined as a micro-convex, as shown in FIG.
  • Rough surface is made by A plurality of circular microprotrusions of different sizes are formed, assuming that each microprotrusion is separated from each other, and the interaction is negligible.
  • a' is the cross-sectional area where the micro-convex body intersects the rigid surface before denaturation
  • a the contact area of the micro-convex body after deformation with the rigid surface is defined as the true contact area.
  • the deformation of the micro-convex is divided into elastic and plastic deformation, and the critical cross-sectional area that determines the elastic or plastic deformation of the micro-convex is determined by the material properties of the gear pair:
  • the lower corners e and p represent elastic and plastic deformation, respectively.
  • E 1 E 2 v 1 v 2 represents the elastic modulus and Poisson's ratio of the two contact tooth faces, respectively
  • H represents the hardness of the softer material
  • H 2.8Y
  • R and ⁇ respectively represent a single microconvex The curvature of the body and the amount of normal deformation.
  • the total elastic deformation force F E and the total plastic deformation force F P of a single cube mesh can be expressed as:
  • the step 3) specifically includes the following steps:
  • the deformation is divided into two stages: full elastic deformation and complete plastic deformation, wherein the plastic deformation stage has a stiffness of zero.
  • the single microprotrusion normal contact stiffness where complete elastic deformation occurs can be expressed as
  • the normal contact stiffness in a single square mesh is obtained by integrating in the fully elastic region.
  • ⁇ denotes the coefficient of static friction
  • ⁇ 1 ⁇ 2 G 1 G 2 represents the Poisson's ratio and shear modulus of the two materials, respectively
  • t and f represent the tangential and normal loads of a single micro-convex, respectively
  • r represents the true contact area radius of a single micro-convex and Where a represents the true contact area.
  • T ⁇ b A r
  • ⁇ b the shear strength of the softer material
  • a r the actual contact area of the square grid
  • the step 4) specifically includes the following steps:
  • the tangential contact stiffness K 2 of the contact tooth surface is
  • the contact stiffness of the tooth surface can be expressed as
  • ⁇ represents the pressure angle of the gear.

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

Abstract

一种基于粗糙表面的直齿轮三维接触刚度计算方法，采用六面体网格划分，要求齿面上每个正方形网格面积相等；基于有限元法计算光滑齿面接触压力分布，提取接触区节点压强；基于分形理论计算单个正方形网格法向接触刚度，切向接触刚度；计算齿面接触刚度。本方法进一步精确了齿轮接触刚度计算模型，解决了原先计算过程中按照赫兹理论光滑接触的弊端，揭示了粗糙齿面接触过程中粗糙度参数对于齿轮刚度特性的影响规律和计算方法。同时基于有限元的三维接触刚度计算模型，比原先二维计算模型或者计算公式更精确，并且考虑齿轮实际装配和制造过程中误差的影响，从而使齿轮接触刚度更加准确，为齿轮动力学分析奠定了基础。

Description

基于粗糙表面的直齿轮三维接触刚度计算方法 技术领域
本发明属于齿轮力学分析与动力学研究领域，具体涉及一种基于粗糙表面的直齿轮三维接触刚度计算方法。
背景技术
齿轮副内部激励是在齿轮副啮合过程中在系统内部产生的，也是齿轮系统动力学研究中的重点和难点。齿轮系统的内部激励主要分为三种：时变刚度激励、时变误差激励和啮合冲击激励。对于时变刚度激励主要由三部分构成：轮齿的弯曲刚度、剪切刚度和接触刚度，目前应用较为广泛的是IOS6336啮合刚度计算公式以及将齿轮的刚度激励等效为傅里叶展式。接触刚度的计算采用的是基于光滑表面的赫兹接触理论进行推导的，因此不能精确模拟实际轮齿粗糙表面接触状态的变形。所以为了建立更有效的齿轮动力学模型，有必要深入研究基于粗糙表面的直齿轮三维接触刚度的计算方法。
发明内容
本发明的目的是提供一种基于粗糙表面的直齿轮三维接触刚度的计算方法，该方法既能有效的精确齿轮接触刚度模型的准确性，同时又充分考虑到由于装配误差、制造误差等因素下的轮齿啮合载荷分布不均的问题。
为了实现上述目的，本发明采用的技术方案是：
S1采用六面体网格划分，要求齿面上每个正方形网格面积A0相等；
S2基于有限元法计算光滑齿面接触压力分布，提取接触区节点压强Pi；
S3基于分形理论计算单个正方形网格法向接触刚度KN，切向接触刚度KT；
S4计算齿面接触刚度K。
所述S1具体包括以下步骤：
齿轮对采用六面体网格划分，约束六面体网格的边长为lmm，因此得到齿轮 对的接触面上均是标准的正方形网格，网格面积A0＝l2mm2。
所述步骤2)具体包括以下步骤：
将网格划分好的齿轮对导入有限元仿真软件中，设置静力学分析环境，通过定义每个齿轮的材料属性、接触对、固定被动齿轮的中心孔，然后在主动齿轮轴线上施加扭矩T，以载荷步的方式得到当前齿轮特定接触面的应力和应变云图，此时提取特定接触面上接触区节点压强Pi。
所述步骤3)具体包括以下步骤：
在有限元分析中使用的网格已经足够小，将提取到的接触区节点压强Pi看作是齿面上正方形网格内的接触压强，因此齿面上不同的正方形网格的接触压力表示为：
Fi＝PiA0         (1)
由于机械加工的齿面具有分形特征，存在局部与整体自相似的特性。粗糙齿面的轮廓由W-M函数表示为：
Figure PCTCN2016107233-appb-000001
其中，n为频率指数，n＝0和nmax是最低及最高截止频率对应的序列，z表示粗糙表面轮廓高度，D表示分形维数，γ表示谱密度的尺寸参数，G表示粗糙度参数，x表示采样长度坐标。采用相同的材料和加工方法得到齿轮样品方块，在三维形貌测量仪上获得齿面形貌点坐标，再采用功率谱密度方法拟合W-M函数得到D和G两个参数。
在齿轮对的接触模型中两个粗糙表面的接触被简化为一个刚性表面与粗糙表面的接触。虚线是中性面表示理想光滑齿轮对的接触位置，将中性面的上方的粗糙表面凸出的部分定义为微凸体。粗糙表面是由无数多个大小不一的圆形微凸体组成，假设每个微凸体之间彼此分立，相互作用忽略不计。当粗糙表面与刚性表面相互接触时，不同微凸体在压力的作用下发生弹性或塑性形变，则微凸体横截面积a′的分布规律满足：
Figure PCTCN2016107233-appb-000002
式中，a′为微凸体变性前与刚性表面相交的横截面积，a表示微凸体变形后与刚性面的接触面积定义为真实接触面积，当微凸体发生弹性变形时，a′＝2a,当发生塑性形变时，a′＝a，D表示分形维数，a′l表示最大横街面积，
Figure PCTCN2016107233-appb-000003
表示域拓展因子，可由超越方程(4)求得：
[ψ(2-D)/2-(1+ψ-D/2)-(2-D)/D]/[(2-D)/D]＝1       (4)
根据微凸体横截面积的不同，将微凸体的变形分为弹性和塑性变形，则决定微凸体发生弹性或塑性变形的临界横截面积是由齿轮对的材料属性决定的：
Figure PCTCN2016107233-appb-000004
H表示较软材料的硬度，H＝2.8Y，Y表示屈服强度值；E*表示等效弹性模量,
Figure PCTCN2016107233-appb-000005
E1、E2、v1、v2分别表示两接触齿面的弹性模量和泊松比。如果微凸体的横截面积a′>a′c，这发生弹性变形，若横截面积a′≤a′c，这发生塑性变形。粗糙表面单个微凸体法向变形量由W-M函数中峰-谷间幅值来表示，
δ＝23-DG(D-1)(lnγ)1/2π(D-2)/2(a′)(2-D)/2      (6)
单个微凸体曲率R为，
Figure PCTCN2016107233-appb-000006
对于单个微凸体的弹性或塑性形变，其法向载荷f与横截面积a′满足如下关系，
Figure PCTCN2016107233-appb-000007
fp＝Ha′         (9)
式中，下角标e和p分别代表弹性和塑性变形，
Figure PCTCN2016107233-appb-000008
表示等效弹性模量，H表示较软材料的硬度，H＝2.8Y，R和δ分别表示单个微凸体的曲率和法向变形量。
单个正方体网格总的弹性变形力FE和总塑性变形力FP表示为：
Figure PCTCN2016107233-appb-000009
Figure PCTCN2016107233-appb-000010
式中符号的意义与上文相同，由此得到单个正方形网格内总的接触压力,
Fi＝FE+FP           (12)
所述步骤3)还包括以下步骤：
根据式(1)和(12)可以得到单个正方形网格当中最大横截面积a′l。
对于单个微凸体，其变形分为完全弹性变形及完全塑性变形两个个阶段，其中塑性变形阶段刚度为0。根据刚度定义，发生完全弹性变形的单个微凸体法向接触刚度可表示为，
Figure PCTCN2016107233-appb-000011
因此通过在完全弹性区域进行积分，得到单个正方形网格内法向接触刚度
Figure PCTCN2016107233-appb-000012
研究表明单个微凸体法切向变形量δt表示为
Figure PCTCN2016107233-appb-000013
式中μ表示静摩擦系数，G′表示结合面的等效剪切模量且满足1/G′＝(2-ν1)/G1+(2-ν2)/G2，其中ν1、ν2、G1、G2分别表示两材料的泊松比和剪切模量，t和f分别表示单个微凸体的切向和法向载荷，r表示单个微凸体的真实接触区域半径且
Figure PCTCN2016107233-appb-000014
a表示真实接触面积。
则单个微凸体的切向刚度可以表示为，
Figure PCTCN2016107233-appb-000015
对于单个微凸体切向载荷t与法向载荷f满足如下关系，
Figure PCTCN2016107233-appb-000016
T＝τbAr，τb表示较软材料的剪切强度，Ar表示正方形网格的实际接触面积
Figure PCTCN2016107233-appb-000017
则单个正方形网格的切向接触刚度KT可以表示为，
Figure PCTCN2016107233-appb-000018
所述步骤4)具体包括以下步骤：
设特定齿面上有m个正方形网格，则接触齿面的法向接触刚度K1表示为，
Figure PCTCN2016107233-appb-000019
接触齿面的切向接触刚度K2为
Figure PCTCN2016107233-appb-000020
则齿面的接触刚度可以表示为
K＝K1cosα+K2sinα          (22)
式中α表示齿轮的压力角。
相对于现有技术，本发明的有益效果为：
本发明提出的基于粗糙表面的直齿轮三维接触刚度计算方法，进一步精确了齿轮接触刚度计算模型，解决了原先计算过程中按照赫兹理论光滑接触的弊端，揭示了粗糙齿面接触过程中粗糙度参数对于齿轮刚度特性的影响规律和计算方法。同时本发明是基于有限元的三维接触刚度计算模型，比原先的二维计算模型或者计算公式更精确，并且考虑到齿轮实际装配和制造过程中误差的影响，从而使齿轮接触刚度更加准确，为齿轮动力学分析奠定了基础。
建立的直齿轮接触刚度模型中首次包含了粗糙度有关参数，间接阐述了加工方式与齿轮接触刚度的关系。
通过特殊有限元网格划分方法，首次将分形理论分别应用于接触区不同正方向网格中，通过对不同接触区内网格的累加，得到了齿轮三维接触刚度，避免了因齿轮对因装配或制造误差引起的偏载对啮合刚度的影响。
采用的技术方案中S1、S2、S3对除了直齿轮以外的其他类型齿轮同样适用。
附图说明
图1是齿轮对简化接触模型图。
图2是单个微凸体接触区域图。
图3是直齿轮接触刚度转化示意图。
图4是所述技术方案的计算流程图。
具体实施方式
下面结合附图1-4和实施例对本发明作进一步说明。
一种基于粗糙表面的直齿轮三维接触刚度计算方法，包括以下步骤：
1)齿轮对特定网格划分，要求各网格面积A0相等；
2)基于有限元法计算光滑齿面接触压力分布，提取接触区节点压强Pi；
3)基于分形理论计算单个正方形网格法向接触刚度KN，切向接触刚度KT；
4)计算齿面接触刚度K。
所述步骤1)具体包括以下步骤：
齿轮对采用六面体网格划分，约束六面体网格的边长为2mm，因此可以得到齿轮对的接触面上均是标准的正方形网格，网格面积A0＝4mm2。
所述步骤2)具体包括以下步骤：
将网格划分好的齿轮对导入有限元仿真软件中，设置静力学分析环境，通过定义每个齿轮的材料属性，接触对，固定被动齿轮的中心孔，然后在主动齿轮轴线上施加扭矩T，以载荷步的方式可以得到当前齿轮特定接触面的应力和应变云图，此时可以提取特定接触面上接触区节点压强Pi。
所述步骤3)具体包括以下步骤：
在有限元分析中使用的网格已经足够小，可以将提取到的接触区节点压强Pi看作是齿面上正方形网格内的接触压强，因此齿面上不同的正方形网格的接触压力可以表示为：
Fi＝PiA0            (1)
由于机械加工的齿面具有分形特征，存在局部与整体自相似的特性。粗糙齿面的轮廓可以由W-M函数表示为：
Figure PCTCN2016107233-appb-000021
其中n为频率指数，n＝0和nmax是最低及最高截止频率对应的序列，z表示粗糙表面轮廓高度，D表示分形维数，γ表示谱密度的尺寸参数，G表示粗糙度参数，x表示采样长度坐标。采用相同的材料和加工方法得到齿轮样品方块，在三维形貌测量仪上获得齿面形貌点坐标，再采用功率谱密度方法拟合W-M函数得到D和G两个参数。
在齿轮对的接触模型中两个粗糙表面的接触被简化为一个刚性表面与粗糙表面的接触，如图1所示。虚线是中性面表示理想光滑齿轮对的接触位置，将中性面的上方的粗糙表面凸出的部分定义为微凸体，如图2所示。粗糙表面是由无 数多个大小不一的圆形微凸体组成，假设每个微凸体之间彼此分立，相互作用忽略不计。当粗糙表面与刚性表面相互接触时，不同微凸体在压力的作用下发生弹性或塑性形变，则微凸体横截面积a′的分布规律满足：
Figure PCTCN2016107233-appb-000022
式中，a′为微凸体变性前与刚性表面相交的横截面积，a表示微凸体变形后与刚性面的接触面积定义为真实接触面积，当微凸体发生弹性变形时，a′＝2a,当发生塑性形变时，a′＝a，D表示分形维数，a′l表示最大横街面积，
Figure PCTCN2016107233-appb-000023
表示域拓展因子，可由超越方程(4)求得：
[ψ(2-D)/2-(1+ψ-D/2)-(2-D)/D]/[(2-D)/D]＝1     (4)
根据微凸体横截面积的不同，将微凸体的变形分为弹性和塑性变形，则决定微凸体发生弹性或塑性变形的临界横截面积是由齿轮对的材料属性决定的：
Figure PCTCN2016107233-appb-000024
H表示较软材料的硬度，H＝2.8Y，Y表示屈服强度值；E*表示等效弹性模量,
Figure PCTCN2016107233-appb-000025
E1 E2 v1 v2分别表示两接触齿面的弹性模量和泊松比。如果微凸体的横截面积a′>a′c，这发生弹性变形，若横截面积a′≤a′c，这发生塑性变形。粗糙表面单个微凸体法向变形量可由W-M函数中峰-谷间幅值来表示，
δ＝23-DG(D-1)(lnγ)1/2π(D-2)/2(a′)(2-D)/2      (6)
单个微凸体曲率R为，
Figure PCTCN2016107233-appb-000026
对于单个微凸体的弹性或塑性形变，其法向载荷f与横截面积a′满足如下关系，
Figure PCTCN2016107233-appb-000027
fp＝Ha′          (9)
式中，下角标e和p分别代表弹性和塑性变形，
Figure PCTCN2016107233-appb-000028
表示等效弹性模量，E1 E2 v1 v2分别表示两接触齿面的弹性模量和泊松比，H表示较软材料的硬度，H＝2.8Y，R和δ分别表示单个微凸体的曲率和法向变形量。
单个正方体网格总的弹性变形力FE和总塑性变形力FP可以表示为：
Figure PCTCN2016107233-appb-000029
Figure PCTCN2016107233-appb-000030
式中符号的意义与上文相同，由此得到单个正方形网格内总的接触压力,
Fi＝FE+FP         (12)
所述步骤3)具体包括以下步骤：
根据式(1)和(12)可以得到单个正方形网格当中最大横截面积a′l。
对于单个微凸体，其变形分为完全弹性变形及完全塑性变形两个个阶段，其中塑性变形阶段刚度为0。根据刚度定义，发生完全弹性变形的单个微凸体法向接触刚度可表示为，
Figure PCTCN2016107233-appb-000031
因此通过在完全弹性区域进行积分，得到单个正方形网格内法向接触刚度
Figure PCTCN2016107233-appb-000032
研究表明单个微凸体法切向变形量δt可表示为
Figure PCTCN2016107233-appb-000033
式中μ表示静摩擦系数，G′表示结合面的等效剪切模量且满足1/G′＝(2-ν1)/G1+(2-ν2)/G2，其中ν1 ν2 G1 G2分别表示两材料的泊松比和剪切模量，t和f分别表示单个微凸体的切向和法向载荷，r表示单个微凸体的真实接触区域半径且
Figure PCTCN2016107233-appb-000034
式中a表示真实接触面积。
则单个微凸体的切向刚度可以表示为，
Figure PCTCN2016107233-appb-000035
对于单个微凸体切向载荷t与法向载荷f满足如下关系，
Figure PCTCN2016107233-appb-000036
其中T＝τbAr，τb表示较软材料的剪切强度，Ar表示正方形网格的实际接触面积
Figure PCTCN2016107233-appb-000037
则单个正方形网格的切向接触刚度KT可以表示为，
Figure PCTCN2016107233-appb-000038
所述步骤4)具体包括以下步骤：
设特定齿面上有m个正方形网格，则接触齿面的法向接触刚度K1表示为，
Figure PCTCN2016107233-appb-000039
接触齿面的切向接触刚度K2为
Figure PCTCN2016107233-appb-000040
根据图3可得齿面的接触刚度可以表示为
K＝K1cosα+K2sinα            (22)
式中α表示齿轮的压力角。

Claims (6)

  1. 基于粗糙表面的直齿轮三维接触刚度计算方法，其特征在于：本方法包括如下步骤，
    S1采用六面体网格划分，要求齿面上每个正方形网格面积A0相等；
    S2基于有限元法计算光滑齿面接触压力分布，提取接触区节点压强Pi；
    S3基于分形理论计算单个正方形网格法向接触刚度KN，切向接触刚度KT；
    S4计算齿面接触刚度K。
  2. 根据权利要求1所述的基于粗糙表面的直齿轮三维接触刚度计算方法，其特征在于：所述S1具体包括以下步骤：
    齿轮对采用六面体网格划分，约束六面体网格的边长为lmm，因此得到齿轮对的接触面上均是标准的正方形网格，网格面积A0＝l2mm2。
  3. 根据权利要求1所述的基于粗糙表面的直齿轮三维接触刚度计算方法，其特征在于：所述步骤2)具体包括以下步骤：
    将网格划分好的齿轮对导入有限元仿真软件中，设置静力学分析环境，通过定义每个齿轮的材料属性、接触对、固定被动齿轮的中心孔，然后在主动齿轮轴线上施加扭矩T，以载荷步的方式得到当前齿轮特定接触面的应力和应变云图，此时提取特定接触面上接触区节点压强Pi。
  4. 根据权利要求1所述的基于粗糙表面的直齿轮三维接触刚度计算方法，其特征在于：所述步骤3)具体包括以下步骤：
    在有限元分析中使用的网格已经足够小，将提取到的接触区节点压强Pi看作是齿面上正方形网格内的接触压强，因此齿面上不同的正方形网格的接触压力表示为：
    Fi＝PiA0            (1)
    由于机械加工的齿面具有分形特征，存在局部与整体自相似的特性；粗糙齿面的轮廓由W-M函数表示为：
    Figure PCTCN2016107233-appb-100001
    其中，n为频率指数，n＝0和nmax是最低及最高截止频率对应的序列，z表示粗糙表面轮廓高度，D表示分形维数，γ表示谱密度的尺寸参数，G表示粗糙度参数，x表示采样长度坐标；采用相同的材料和加工方法得到齿轮样品方块，在三维形貌测量仪上获得齿面形貌点坐标，再采用功率谱密度方法拟合W-M函数得到D和G两个参数；
    在齿轮对的接触模型中两个粗糙表面的接触被简化为一个刚性表面与粗糙表面的接触；虚线是中性面表示理想光滑齿轮对的接触位置，将中性面的上方的粗糙表面凸出的部分定义为微凸体；粗糙表面是由无数多个大小不一的圆形微凸体组成，假设每个微凸体之间彼此分立，相互作用忽略不计；当粗糙表面与刚性表面相互接触时，不同微凸体在压力的作用下发生弹性或塑性形变，则微凸体横截面积a′的分布规律满足：
    Figure PCTCN2016107233-appb-100002
    式中，a′为微凸体变性前与刚性表面相交的横截面积，a表示微凸体变形后与刚性面的接触面积定义为真实接触面积，当微凸体发生弹性变形时，a′＝2a,当发生塑性形变时，a′＝a，D表示分形维数，a′l表示最大横街面积，
    Figure PCTCN2016107233-appb-100003
    表示域拓展因子，可由超越方程(4)求得：
    [ψ(2-D)/2-(1+ψ-D/2)-(2-D)/D]/[(2-D)/D]＝1        (4)
    根据微凸体横截面积的不同，将微凸体的变形分为弹性和塑性变形，则决定微凸体发生弹性或塑性变形的临界横截面积是由齿轮对的材料属性决定的：
    Figure PCTCN2016107233-appb-100004
    H表示较软材料的硬度，H＝2.8Y，Y表示屈服强度值；E*表示等效弹性模量,
    Figure PCTCN2016107233-appb-100005
    E1、E2、v1、v2分别表示两接触齿面的弹性模量和泊松比；如果微
    凸体的横截面积a′>a′c，这发生弹性变形，若横截面积a′≤a′c，这发生塑性变形；粗糙表面单个微凸体法向变形量由W-M函数中峰-谷间幅值来表示，
    Figure PCTCN2016107233-appb-100006
    单个微凸体曲率R为，
    Figure PCTCN2016107233-appb-100007
    对于单个微凸体的弹性或塑性形变，其法向载荷f与横截面积a′满足如下关系，
    Figure PCTCN2016107233-appb-100008
    fp＝Ha′            (9)
    式中，下角标e和p分别代表弹性和塑性变形，
    Figure PCTCN2016107233-appb-100009
    表示等效弹性模量，H表示较软材料的硬度，H＝2.8Y，R和δ分别表示单个微凸体的曲率和法向变形量；
    单个正方体网格总的弹性变形力FE和总塑性变形力FP表示为：
    Figure PCTCN2016107233-appb-100010
    Figure PCTCN2016107233-appb-100011
    由此得到单个正方形网格内总的接触压力,
    Fi＝FE+FP              (12)
    。
  5. 根据权利要求4所述的基于粗糙表面的直齿轮三维接触刚度计算方法，其特征在于：所述步骤3)还包括以下步骤：
    根据式(1)和(12)可以得到单个正方形网格当中最大横截面积a′l；
    对于单个微凸体，其变形分为完全弹性变形及完全塑性变形两个个阶段，其中塑性变形阶段刚度为0；根据刚度定义，发生完全弹性变形的单个微凸体法向接触刚度可表示为，
    Figure PCTCN2016107233-appb-100012
    因此通过在完全弹性区域进行积分，得到单个正方形网格内法向接触刚度
    Figure PCTCN2016107233-appb-100013
    研究表明单个微凸体法切向变形量δt表示为
    Figure PCTCN2016107233-appb-100014
    式中μ表示静摩擦系数，G′表示结合面的等效剪切模量且满足1/G′＝(2-ν1)/G1+(2-ν2)/G2，其中ν1、ν2、G1、G2分别表示两材料的泊松比和剪切模量，t和f分别表示单个微凸体的切向和法向载荷，r表示单个微凸体的真实接触区域半径且
    Figure PCTCN2016107233-appb-100015
    a表示真实接触面积；
    则单个微凸体的切向刚度可以表示为，
    Figure PCTCN2016107233-appb-100016
    对于单个微凸体切向载荷t与法向载荷f满足如下关系，
    Figure PCTCN2016107233-appb-100017
    T＝τbAr，τb表示较软材料的剪切强度，Ar表示正方形网格的实际接触面积
    Figure PCTCN2016107233-appb-100018
    则单个正方形网格的切向接触刚度KT表示为，
    Figure PCTCN2016107233-appb-100019
  6. 根据权利要求1所述的基于粗糙表面的直齿轮三维接触刚度计算方法，其特征在于：所述步骤4)具体包括以下步骤：
    设特定齿面上有m个正方形网格，则接触齿面的法向接触刚度K1表示为，
    Figure PCTCN2016107233-appb-100020
    接触齿面的切向接触刚度K2为
    Figure PCTCN2016107233-appb-100021
    则齿面的接触刚度可以表示为
    K＝K1cosα+K2sinα           (22)
    式中α表示齿轮的压力角。
PCT/CN2016/107233 2016-11-09 2016-11-25 基于粗糙表面的直齿轮三维接触刚度计算方法 WO2018086160A1 (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US15/798,320 US10088386B2 (en) 2016-11-09 2017-10-30 Device and method for measuring three-dimensional contact stiffness of spur gear based on rough surface

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610985604.5 2016-11-09
CN201610985604.5A CN106844818B (zh) 2016-11-09 2016-11-09 基于粗糙表面的直齿轮三维接触刚度计算方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/798,320 Continuation-In-Part US10088386B2 (en) 2016-11-09 2017-10-30 Device and method for measuring three-dimensional contact stiffness of spur gear based on rough surface

Publications (1)

Publication Number Publication Date
WO2018086160A1 true WO2018086160A1 (zh) 2018-05-17

Family

ID=59146316

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2016/107233 WO2018086160A1 (zh) 2016-11-09 2016-11-25 基于粗糙表面的直齿轮三维接触刚度计算方法

Country Status (2)

Country Link
CN (1) CN106844818B (zh)
WO (1) WO2018086160A1 (zh)

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108875198A (zh) * 2018-06-13 2018-11-23 株洲时代新材料科技股份有限公司 一种基于流体有限元的大变形弹性体产品静刚度计算方法
CN108959790A (zh) * 2018-07-12 2018-12-07 中南大学 一种控制齿轮单双齿啮合刚度平缓性的方法
CN109271688A (zh) * 2018-09-03 2019-01-25 武汉理工大学 高速齿轮啮合传动温度场的精确计算方法
CN109271711A (zh) * 2018-09-25 2019-01-25 重庆大学 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法
CN109284544A (zh) * 2018-09-04 2019-01-29 中南大学 螺旋锥齿轮形性协同制造参数驱动混合反调及控制方法
CN109408857A (zh) * 2018-09-04 2019-03-01 中南大学 螺旋锥齿轮形性协同制造的智能参数驱动模块化设计方法
CN109446655A (zh) * 2018-10-30 2019-03-08 东北大学 一种双分形结合面法向接触刚度的三维分形预测方法
CN109492307A (zh) * 2018-11-12 2019-03-19 中南大学 一种弧齿锥齿轮齿面载荷接触性能参数的数值计算方法
CN109783840A (zh) * 2018-11-30 2019-05-21 华南理工大学 一种直齿圆柱齿轮内啮合齿轮副时变啮合刚度的计算方法
CN110532669A (zh) * 2019-08-26 2019-12-03 西安建筑科技大学 一种用于机械结合面接触刚度建模的方法
CN110688717A (zh) * 2019-09-25 2020-01-14 西北工业大学 齿轮齿条结构的可靠性分析方法及装置
CN110705161A (zh) * 2019-09-27 2020-01-17 湘潭大学 考虑摩擦热分布的机械密封不对中状态热仿真分析方法
CN110826273A (zh) * 2019-10-28 2020-02-21 长安大学 一种考虑浮动特性的行星传动多体齿轮承载接触特性分析方法
CN110940587A (zh) * 2019-12-17 2020-03-31 北京工业大学 一种基于多尺度理论的床身-地基结合面接触刚度计算方法
CN111209705A (zh) * 2020-01-15 2020-05-29 同济大学 一种基于有限元的玻璃升降器三维柔性装配公差预测方法
CN111353234A (zh) * 2020-03-11 2020-06-30 西安建筑科技大学 一种考虑微凸体相互作用的粗糙表面接触模型构建方法及系统
CN111382481A (zh) * 2020-03-11 2020-07-07 西安建筑科技大学 一种考虑微凸体相互作用的结合面法向刚度建模方法及系统
CN112257193A (zh) * 2020-09-27 2021-01-22 沈阳建筑大学 一种曲面织构的表征方法
CN112417593A (zh) * 2020-11-18 2021-02-26 西北工业大学 一种航空发动机圆弧端齿连接装配安装角度的优化方法
CN112434431A (zh) * 2020-11-27 2021-03-02 西安交通大学 一种内啮合直齿轮磨损下的啮合刚度的建模方法
CN112434432A (zh) * 2020-11-27 2021-03-02 西安交通大学 一种外啮合直齿轮磨损下的啮合刚度的建模方法
CN113010978A (zh) * 2021-02-09 2021-06-22 中国航发哈尔滨东安发动机有限公司 一种基于动态仿真的航空直齿轮修形方法
CN113076606A (zh) * 2021-03-24 2021-07-06 西北工业大学 考虑接触变形的航空管路扩口接头泄漏率计算方法及系统
CN113139313A (zh) * 2021-04-25 2021-07-20 重庆大学 基于有限元理论的等效轻量化结构分析方法
CN113239430A (zh) * 2021-04-29 2021-08-10 西安交通大学 一种输电钢管塔用三环板节点环板作用力计算方法
CN113343383A (zh) * 2021-06-01 2021-09-03 江苏科技大学 一种机器人关节谐波减速器动态传动误差预测方法
CN113392544A (zh) * 2021-05-28 2021-09-14 东北林业大学 一种基于变形协调理论的行星螺纹滚柱轴承接触载荷计算方法
CN113704910A (zh) * 2021-08-20 2021-11-26 武汉大学 一种微米级粗糙度螺栓结合面真实接触率的计算方法
CN114297822A (zh) * 2021-11-30 2022-04-08 国网辽宁省电力有限公司锦州供电公司 高压开关柜梅花触头触指压力及温度计算方法
CN114354104A (zh) * 2021-12-03 2022-04-15 北京工业大学 一种测量余弦波规则形貌结合面切向刚度的装置及方法
CN115438446A (zh) * 2022-07-01 2022-12-06 上海交通大学 一种考虑宏微观形貌特征的界面切向接触刚度预测方法
CN117094200A (zh) * 2023-10-17 2023-11-21 安徽大学 一种考虑不对中误差的齿轮时变啮合刚度计算方法
CN117634057A (zh) * 2023-10-18 2024-03-01 南京航空航天大学 一种含剥落故障的弧齿锥齿轮时变啮合刚度计算方法
CN115438446B (zh) * 2022-07-01 2024-05-28 上海交通大学 一种考虑宏微观形貌特征的界面切向接触刚度预测方法

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609231B (zh) * 2017-08-24 2019-08-30 中南大学 一种蜗杆磨削齿轮表面微观形貌仿真方法和系统
CN107545127B (zh) * 2017-10-13 2021-01-05 北京工业大学 一种考虑接触的工业机器人关节刚度建模方法
CN108168492A (zh) * 2017-12-17 2018-06-15 北京工业大学 一种基于三维微观表面形貌的接触特性分析方法
CN109059838B (zh) * 2018-08-28 2020-08-11 沈阳建筑大学 一种具有表面微织构的粗糙表面参数表征方法
CN109187008B (zh) * 2018-09-13 2020-11-03 中国航发湖南动力机械研究所 齿轮毂加扭装置
CN109472071B (zh) * 2018-10-30 2023-06-16 东北大学 一种结合面静摩擦系数三维分形预测方法
CN109829224B (zh) * 2019-01-25 2020-08-11 北京科技大学 一种计算粗糙弹性界面切向接触刚度的方法
CN109887079B (zh) * 2019-03-21 2022-09-20 中南大学 螺旋锥齿轮三维建模方法
CN110188508B (zh) * 2019-06-19 2023-06-20 江苏太平洋齿轮传动有限公司 基于实际载荷的齿面接触斑点齿轮的设计方法
CN111999054B (zh) * 2020-08-17 2021-10-22 湖南大学 一种考虑实际齿面弹塑性接触的圆柱直齿磨损模拟方法
CN113204890A (zh) * 2021-05-21 2021-08-03 北京工业大学 一种基于齿面形貌参数的渐开线齿廓改进方程设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003194754A (ja) * 2001-12-28 2003-07-09 Toyota Motor Corp 熱処理シミュレーションによる歪の解析方法
CN104573196A (zh) * 2014-12-18 2015-04-29 西安交通大学 一种斜齿圆柱齿轮时变啮合刚度解析计算方法
CN105224744A (zh) * 2015-09-29 2016-01-06 西安交通大学 一种基于啮合刚度的剥落齿轮啮合模型
CN105302991A (zh) * 2015-11-17 2016-02-03 天津百利机械装备研究院有限公司 一种渐开线圆柱斜齿轮机构参数化分析方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003194754A (ja) * 2001-12-28 2003-07-09 Toyota Motor Corp 熱処理シミュレーションによる歪の解析方法
CN104573196A (zh) * 2014-12-18 2015-04-29 西安交通大学 一种斜齿圆柱齿轮时变啮合刚度解析计算方法
CN105224744A (zh) * 2015-09-29 2016-01-06 西安交通大学 一种基于啮合刚度的剥落齿轮啮合模型
CN105302991A (zh) * 2015-11-17 2016-02-03 天津百利机械装备研究院有限公司 一种渐开线圆柱斜齿轮机构参数化分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Surface Morphology and Stiffness Model Based on Fractal Theory of Inner Bevel Planetary Gear", SCIENCE -ENGINEERING (B), CHINA MASTER'S THESES FULL-TEXT DATABASE, 15 June 2014 (2014-06-15), ISSN: 1674-0246 *

Cited By (56)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108875198A (zh) * 2018-06-13 2018-11-23 株洲时代新材料科技股份有限公司 一种基于流体有限元的大变形弹性体产品静刚度计算方法
CN108875198B (zh) * 2018-06-13 2022-12-02 株洲时代新材料科技股份有限公司 一种基于流体有限元的大变形弹性体产品静刚度计算方法
CN108959790A (zh) * 2018-07-12 2018-12-07 中南大学 一种控制齿轮单双齿啮合刚度平缓性的方法
CN109271688A (zh) * 2018-09-03 2019-01-25 武汉理工大学 高速齿轮啮合传动温度场的精确计算方法
CN109271688B (zh) * 2018-09-03 2022-11-18 武汉理工大学 高速齿轮啮合传动温度场的精确计算方法
CN109284544A (zh) * 2018-09-04 2019-01-29 中南大学 螺旋锥齿轮形性协同制造参数驱动混合反调及控制方法
CN109408857A (zh) * 2018-09-04 2019-03-01 中南大学 螺旋锥齿轮形性协同制造的智能参数驱动模块化设计方法
CN109408857B (zh) * 2018-09-04 2023-03-24 中南大学 螺旋锥齿轮形性协同制造的智能参数驱动模块化设计方法
CN109284544B (zh) * 2018-09-04 2023-03-31 中南大学 螺旋锥齿轮形性协同制造参数驱动混合反调及控制方法
CN109271711A (zh) * 2018-09-25 2019-01-25 重庆大学 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法
CN109446655A (zh) * 2018-10-30 2019-03-08 东北大学 一种双分形结合面法向接触刚度的三维分形预测方法
CN109492307B (zh) * 2018-11-12 2023-04-07 中南大学 一种弧齿锥齿轮齿面载荷接触性能参数的数值计算方法
CN109492307A (zh) * 2018-11-12 2019-03-19 中南大学 一种弧齿锥齿轮齿面载荷接触性能参数的数值计算方法
CN109783840B (zh) * 2018-11-30 2022-10-25 华南理工大学 一种直齿圆柱齿轮内啮合齿轮副时变啮合刚度的计算方法
CN109783840A (zh) * 2018-11-30 2019-05-21 华南理工大学 一种直齿圆柱齿轮内啮合齿轮副时变啮合刚度的计算方法
CN110532669B (zh) * 2019-08-26 2022-09-13 西安建筑科技大学 一种用于机械结合面接触刚度建模的方法
CN110532669A (zh) * 2019-08-26 2019-12-03 西安建筑科技大学 一种用于机械结合面接触刚度建模的方法
CN110688717A (zh) * 2019-09-25 2020-01-14 西北工业大学 齿轮齿条结构的可靠性分析方法及装置
CN110705161A (zh) * 2019-09-27 2020-01-17 湘潭大学 考虑摩擦热分布的机械密封不对中状态热仿真分析方法
CN110826273B (zh) * 2019-10-28 2024-02-02 长安大学 一种考虑浮动特性的行星传动多体齿轮承载接触特性分析方法
CN110826273A (zh) * 2019-10-28 2020-02-21 长安大学 一种考虑浮动特性的行星传动多体齿轮承载接触特性分析方法
CN110940587A (zh) * 2019-12-17 2020-03-31 北京工业大学 一种基于多尺度理论的床身-地基结合面接触刚度计算方法
CN111209705B (zh) * 2020-01-15 2024-04-26 同济大学 一种基于有限元的玻璃升降器三维柔性装配公差预测方法
CN111209705A (zh) * 2020-01-15 2020-05-29 同济大学 一种基于有限元的玻璃升降器三维柔性装配公差预测方法
CN111382481B (zh) * 2020-03-11 2024-02-13 深圳万知达科技有限公司 一种考虑微凸体相互作用的结合面法向刚度建模方法及系统
CN111353234B (zh) * 2020-03-11 2023-03-28 西安建筑科技大学 一种考虑微凸体相互作用的粗糙表面接触模型构建方法及系统
CN111382481A (zh) * 2020-03-11 2020-07-07 西安建筑科技大学 一种考虑微凸体相互作用的结合面法向刚度建模方法及系统
CN111353234A (zh) * 2020-03-11 2020-06-30 西安建筑科技大学 一种考虑微凸体相互作用的粗糙表面接触模型构建方法及系统
CN112257193B (zh) * 2020-09-27 2023-09-26 沈阳建筑大学 一种曲面织构的表征方法
CN112257193A (zh) * 2020-09-27 2021-01-22 沈阳建筑大学 一种曲面织构的表征方法
CN112417593B (zh) * 2020-11-18 2024-02-23 西北工业大学 一种航空发动机圆弧端齿连接装配安装角度的优化方法
CN112417593A (zh) * 2020-11-18 2021-02-26 西北工业大学 一种航空发动机圆弧端齿连接装配安装角度的优化方法
CN112434432A (zh) * 2020-11-27 2021-03-02 西安交通大学 一种外啮合直齿轮磨损下的啮合刚度的建模方法
CN112434431B (zh) * 2020-11-27 2024-04-05 西安交通大学 一种内啮合直齿轮磨损下的啮合刚度的建模方法
CN112434432B (zh) * 2020-11-27 2024-04-05 西安交通大学 一种外啮合直齿轮磨损下的啮合刚度的建模方法
CN112434431A (zh) * 2020-11-27 2021-03-02 西安交通大学 一种内啮合直齿轮磨损下的啮合刚度的建模方法
CN113010978A (zh) * 2021-02-09 2021-06-22 中国航发哈尔滨东安发动机有限公司 一种基于动态仿真的航空直齿轮修形方法
CN113076606A (zh) * 2021-03-24 2021-07-06 西北工业大学 考虑接触变形的航空管路扩口接头泄漏率计算方法及系统
CN113076606B (zh) * 2021-03-24 2024-03-26 西北工业大学 考虑接触变形的航空管路扩口接头泄漏率计算方法及系统
CN113139313A (zh) * 2021-04-25 2021-07-20 重庆大学 基于有限元理论的等效轻量化结构分析方法
CN113139313B (zh) * 2021-04-25 2023-03-14 重庆大学 基于有限元理论的等效轻量化结构分析方法
CN113239430A (zh) * 2021-04-29 2021-08-10 西安交通大学 一种输电钢管塔用三环板节点环板作用力计算方法
CN113239430B (zh) * 2021-04-29 2022-12-09 西安交通大学 一种输电钢管塔用三环板节点环板作用力计算方法
CN113392544A (zh) * 2021-05-28 2021-09-14 东北林业大学 一种基于变形协调理论的行星螺纹滚柱轴承接触载荷计算方法
CN113392544B (zh) * 2021-05-28 2022-08-26 东北林业大学 一种基于变形协调理论的行星螺纹滚柱轴承接触载荷计算方法
CN113343383A (zh) * 2021-06-01 2021-09-03 江苏科技大学 一种机器人关节谐波减速器动态传动误差预测方法
CN113704910A (zh) * 2021-08-20 2021-11-26 武汉大学 一种微米级粗糙度螺栓结合面真实接触率的计算方法
CN113704910B (zh) * 2021-08-20 2024-04-30 武汉大学 一种微米级粗糙度螺栓结合面真实接触率的计算方法
CN114297822B (zh) * 2021-11-30 2024-04-16 国网辽宁省电力有限公司锦州供电公司 高压开关柜梅花触头触指压力及温度计算方法
CN114297822A (zh) * 2021-11-30 2022-04-08 国网辽宁省电力有限公司锦州供电公司 高压开关柜梅花触头触指压力及温度计算方法
CN114354104A (zh) * 2021-12-03 2022-04-15 北京工业大学 一种测量余弦波规则形貌结合面切向刚度的装置及方法
CN115438446A (zh) * 2022-07-01 2022-12-06 上海交通大学 一种考虑宏微观形貌特征的界面切向接触刚度预测方法
CN115438446B (zh) * 2022-07-01 2024-05-28 上海交通大学 一种考虑宏微观形貌特征的界面切向接触刚度预测方法
CN117094200B (zh) * 2023-10-17 2024-01-16 安徽大学 一种考虑不对中误差的齿轮时变啮合刚度计算方法
CN117094200A (zh) * 2023-10-17 2023-11-21 安徽大学 一种考虑不对中误差的齿轮时变啮合刚度计算方法
CN117634057A (zh) * 2023-10-18 2024-03-01 南京航空航天大学 一种含剥落故障的弧齿锥齿轮时变啮合刚度计算方法

Also Published As

Publication number Publication date
CN106844818A (zh) 2017-06-13
CN106844818B (zh) 2018-03-30

Similar Documents

Publication Publication Date Title
WO2018086160A1 (zh) 基于粗糙表面的直齿轮三维接触刚度计算方法
US10088386B2 (en) Device and method for measuring three-dimensional contact stiffness of spur gear based on rough surface
CN106709207B (zh) 一种考虑粗糙表面微凸体相互作用影响的确定受载结合部法向接触刚度的方法
CN107545127B (zh) 一种考虑接触的工业机器人关节刚度建模方法
CN104077440B (zh) 一种基于曲面拟合的确定结合面接触面积及刚度的方法
CN106529038B (zh) 一种从微-宏尺度模型识别螺栓连接结合部切向阻尼特性的方法
Marimuthu et al. Investigation of load carrying capacity of asymmetric high contact ratio spur gear based on load sharing using direct gear design approach
CN110188508A (zh) 基于实际载荷的齿面接触斑点齿轮的设计方法
CN107704657A (zh) 一种由残余应力释放引起的板类零件翘曲变形的计算方法
CN109472071B (zh) 一种结合面静摩擦系数三维分形预测方法
Bruzzone et al. 2D nonlinear and non-Hertzian gear teeth deflection model for static transmission error calculation
CN109829224B (zh) 一种计算粗糙弹性界面切向接触刚度的方法
CN108846189B (zh) 一种齿轮副啮合特性分析方法
Yang et al. Numerical simulation of 3D rough surfaces and analysis of interfacial contact characteristics
CN106777463A (zh) 一种基于三维修正分形理论的单面接触刀柄‑主轴结合部刚度分析方法
Urzica et al. Simulation of the non-gaussian roughness with specified values for the high order moments
CN116187066A (zh) 一种适用于磨削、喷丸表面的分层叠加重构建模方法
Rodríguez-Arós et al. A bending–stretching model in adhesive contact for elastic rods obtained by using asymptotic methods
CN115659611A (zh) 一种非高斯分布下装配结合面接触热阻分析方法
CN110532669B (zh) 一种用于机械结合面接触刚度建模的方法
Yang et al. Fractal characteristic evaluation and interpolation reconstruction for surface topography of drilled composite hole wall
Damtie et al. Contact stress analysis of involute spur gear by Finite Element Method (FEM)
Grigorenko et al. On the equilibrium of nonthin cylindrical shells with a dent
CN113704910B (zh) 一种微米级粗糙度螺栓结合面真实接触率的计算方法
CN112214845A (zh) 一种基于新的点蚀模型的齿轮时变啮合刚度计算方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16921454

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16921454

Country of ref document: EP

Kind code of ref document: A1

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 08.10.2019)

122 Ep: pct application non-entry in european phase

Ref document number: 16921454

Country of ref document: EP

Kind code of ref document: A1