CN113868801B - 一种基于全齿面搜索的轮齿接触分析方法 - Google Patents
一种基于全齿面搜索的轮齿接触分析方法 Download PDFInfo
- Publication number
- CN113868801B CN113868801B CN202111167325.5A CN202111167325A CN113868801B CN 113868801 B CN113868801 B CN 113868801B CN 202111167325 A CN202111167325 A CN 202111167325A CN 113868801 B CN113868801 B CN 113868801B
- Authority
- CN
- China
- Prior art keywords
- point
- tooth surface
- tooth
- gear
- coordinate system
- 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
Images
Classifications
-
- 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
-
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Gears, Cams (AREA)
Abstract
本发明涉及一种基于全齿面搜索的轮齿接触分析方法,属于齿轮接触分析技术领域。在建立齿轮真实齿面的基础上,将齿轮真实齿面进行网格划分,计算网格点处小齿轮齿面点与大齿轮齿面点的距离。给定小齿轮转角,调整大齿轮转角,找到小齿轮齿面和大齿轮齿面之间的最小值,将此时的齿面点作为瞬时接触点。给定不同的小齿轮转角,得到不同瞬时接触点,形成接触轨迹。根据齿轮转角和齿面间最小距离计算传动误差,两个齿面间距离值小于滚检实验涂层厚度的区域即为接触印痕。本方法不受特定齿面的限制,且该方法在进行轮齿接触分析时可考虑到齿面误差,及单对齿面同时存在多个啮合点的情况,并能求解同一时刻单对齿面的所有啮合点。
Description
技术领域
本发明属于齿轮接触分析技术领域,具体涉及一种基于全齿面搜索的轮齿接触分析方法。
背景技术
轮齿接触分析方法(Tooth Contact Analysis)用于计算求解齿轮副的接触印痕与静态传动误差,是评价齿轮齿面设计质量、加工质量和啮合性能的重要手段。
Litvin教授提出的基于两刚性光滑齿面位矢和法矢相同的接触分析算法,是目前应用最广泛的传统分析方法,适应于两光滑齿面构成的点接触齿轮副。其核心是根据两齿轮齿面点处具有共同的法矢和位矢,构造出五个非线性方程(其中位矢相同决定3个方程,法矢相同决定2个方程),联立求解确定啮合点位置及齿轮转角。这种方法在求解时,给定小齿轮的位置转角,只能求解一个瞬时接触点或者接触线(线接触情况),无法考虑齿面局部存在干涉的情况,也不能考虑由于加工误差引起的齿面接触不连续等情况。
方宗德教授基于三坐标测量的齿面偏差,经过三次B样条拟合得到法向偏差曲面,通过理论齿面与偏差曲面叠加构造斜齿轮实际齿面,推导了齿面法矢、齿面接触条件与边缘接触条件,实现了对各种实际齿面进行接触分析仿真。该方法可用于各种加工误差和修形条件下的齿面接触分析仿真,为圆柱齿轮实际齿面接触分析提供了方法。但该方法仍然基于了两齿面法矢和位矢相等的原理进行求解,不能进行全齿面搜索。
孙殿柱教授考虑轮齿热变形和齿面加工误差,将曲面造型技术、CAGD理论和经典齿轮啮合原理交叉融合,针对螺旋圆柱齿轮提出了真实齿面瞬时啮合点的两种求解算法:零间隙法和相切法。相切法计算精度高、速度快,但是对非线性方程组的初值选取很敏感。零间隙法求解速度较慢,但对初值条件无要求,而且可检验啮合点之外是否有相交干涉。两种方法可互补使用。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于全齿面搜索的轮齿接触分析方法,以解决目前传统齿轮轮齿接触分析方法无法考虑单对齿面同时存在多个啮合点的问题。
技术方案
一种基于全齿面搜索的轮齿接触分析方法,其特征在于步骤如下:
步骤1:建立随动坐标系S1和S2,所述坐标系S1的坐标轴为x1,y1,z1,所述坐标系S2的坐标轴为x2,y2,z2,;S1与小齿轮刚性连接,其z1轴与小齿轮轴线重合;同样,S2与大齿轮刚性连接,其z2轴与大齿轮轴重合;
小齿轮齿面Σ(1)在随动坐标系S1下的齿面方程为:
同样,大齿轮齿面Σ(2)在随动坐标系S2下的齿面方程为:
步骤2:选取初始位置:设定初始啮合位置为Oa和Ob重合时的位置,其中初始点Oa和Ob分别位于小齿轮和大齿轮齿宽中点且距离齿顶面为1倍模数处;小齿轮和大齿轮的转角和从初始位置记,取小齿轮初始转角为大齿轮初始转角为即默认两齿轮从初始位置开始求解接触点;
步骤3:划分网格:以初始位置小齿轮齿面Σ(1)上初始点Oa为原点,建立与小齿轮刚性连接的随动坐标系Sa,坐标轴为xa,ya,za,坐标系平面xaya与小齿轮齿面Σ(1)相切于原点Oa,轴xa的方向与小齿轮齿向方向相同;在坐标系平面xaya上划分m行n列的网格T(1),网格横向与轴xa平行,纵向与轴ya平行;同理,以Ob为原点,建立与大齿轮刚性连接的随动坐标系Sb,坐标轴为xb,yb,zb,划分网格T(2);定义T(1)或T(2)上位于每个网格几何中心的点为网格点,其余点为非网格点;由网格点或非网格点在齿轮齿面上投影获得的点为投影点;在后面的距离搜索时将以T(1)上的网格点为索引;
为便于后续距离搜索时的坐标转换,在初始位置处,即时,建立与坐标系Sa重合的固定坐标系Sp,坐标轴为xp,yp,zp,同理建立与坐标系Sb重合的固定坐标系Sg,坐标轴为xg,yg,zg;由于初始位置处Oa和Ob重合,由随动坐标系Sa、Sb各轴关系可知坐标系Sp绕轴xp旋转180°即为坐标系Sg;
步骤4:计算齿面点Z坐标
计算小齿轮齿面Σ(1)在坐标系Sa下的齿面方程为:
由此可计算两齿轮齿面Σ(1)或Σ(2)上任意点P在相应随动坐标系Sa或Sb下的Z坐标,可通过任意点P在其周围投影点的Z坐标的插值获得;因此在距离搜索时,只有网格点在齿面上的投影点是通过齿面方程求其Z坐标的,非网格点在齿面上的投影点其Z坐标通过二维插值方法获得;
步骤5:齿面距离计算
设T(1)上的某网格点在小齿轮齿面Σ(1)上的投影点为点点到T(2)平面的投影点记为M,此时点M在网格平面T(2)上通常为非网格点,因此在求点M到大齿轮齿面Σ(2)上的投影点时需要用到前述的插值;为计算点到点的距离,需要将两齿面方程置于同一坐标系下;点在坐标系Sb下的坐标向量可表示为:
其中为坐标系Sa下的矢量方程;Mba是Sa到Sb的坐标转换矩阵,转换方式为:Sa→Sp→Sg→Sb;Ha是小齿轮初始点M(1)到大齿轮初始点M(2)的矢量;ib,jb,kb分别为Sb下三个轴的单位向量;
假设小齿轮初始转动角度为按照传动比计算大齿轮初始旋转角度为计算m行n列的网格T(1)上所有网格点沿za向小齿轮齿面Σ(1)的全部投影点;依据上述齿面距离计算方法,计算小轮齿面Σ(1)上全部投影点到大轮齿面Σ(2)上对应点的距离;寻找其中最小距离,记为Δzba,min;
如果Δzba,min>ε,则说明两齿轮齿面没有互相接触,需要迭代求解瞬时接触点;如果Δzba,min≤ε,则说明两齿轮齿面存在相互接触甚至相互嵌入的情况;ε为给定容差,用于判断两齿面的状态;
步骤6:迭代求解瞬时接触点
如果Δzba,min>ε,则需要做下面调整,否则直接计算传动误差;最小距离处小轮齿面Σ(1)上点为K,大轮齿面Σ(2)上点为Q,记点K和点Q为齿面接近点;
Δzba,min在Q点线速度VQ上的投影设为ΔV,则
式中,QC为点Q到大齿轮轴的最短距离;
之后,保持小齿轮转动角度不变,使大齿轮转动角度为继续搜索齿面接近点,计算新的如此不断迭代调整大齿轮转动角度求解齿面最小距离;当最小距离Δzba,min小于等于给定容差ε时,此时认为齿面已经接触,其对应的齿面接近点为接触点;
当齿面接触上时,根据传动误差的定义计算传动误差;
寻找全部瞬时接触点:为模拟滚检涂层厚度,T(1)上每个网格点在小齿轮齿面Σ(1)上的投影点沿外法线方向移动距离ξ;以新的投影点为索引,再次搜索两齿面间距离Δzba,此次搜索不再迭代调整大齿轮转角;齿面距离Δzba小于等于滚检实验涂层厚度ξ的齿面点都视为接触点;
步骤7:以增量分别沿正向和负向改变小齿轮转角再按照传动比计算大齿轮转角按照步骤5和6所述的齿面距离计算方法和瞬时接触点求解方法,计算求解下一个啮合位置处的齿面接触点,并由此计算传动误差、瞬时接触线或接触区域;当小齿轮转角经历[-2π/N1,2π/N1]后则停止计算,其中N1是小齿轮齿数。
一种计算机***,其特征在于包括:一个或多个处理器,计算机可读存储介质,用于存储一个或多个程序,其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现上述的方法。
有益效果
本发明提出了一种基于全齿面搜索的轮齿接触分析方法,在建立齿轮真实齿面的基础上,将齿轮真实齿面进行网格划分,计算网格点处小齿轮齿面点与大齿轮齿面点的距离。给定小齿轮转角,调整大齿轮转角,找到小齿轮齿面和大齿轮齿面之间的最小值,将此时的齿面点作为瞬时接触点。给定不同的小齿轮转角,得到不同瞬时接触点,形成接触轨迹。根据齿轮转角和齿面间最小距离计算传动误差,两个齿面间距离值小于滚检实验涂层厚度的区域即为接触印痕。
该方法可适用于不同类型齿轮齿面的接触分析,其算法不受特定齿面的限制。该方法提到的齿面是将齿轮加工时产生的齿面误差通过样条拟合后与理论齿面叠加构造而成。因此本方法在进行轮齿接触分析时可考虑到齿面误差,具有更接近实际的结果。此外,本方法可以考虑到单对齿面同时存在多个啮合点的情况,并能求解同一时刻单对齿面的所有啮合点,弥补了现有传统轮齿接触分析方法只能计算一个啮合点的不足。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1是本发明基于全齿面搜索接触分析方法的坐标系S1、S2及初始点Oa和Ob的示意图,图中为表示清楚将两齿轮齿间距增大;
图2是本发明基于全齿面搜索接触分析方法的小齿轮网格划分示意图;
图3是本发明基于全齿面搜索接触分析方法的随动坐标系Sa、Sb及固定坐标系Sp、Sg各轴关系示意图;
图4是本发明基于全齿面搜索接触分析方法的小齿轮齿面上投影点Z坐标计算方法示意图,为方便表示将小齿轮齿面Σ(1)的曲率增大;
图5是本发明基于全齿面搜索接触分析方法由齿面接近点计算ΔV的示意图;
图6是本发明基于全齿面搜索接触分析方法两齿面距离搜索的示意图;
图7是本发明基于全齿面搜索接触分析方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图7所示为本发明的流程图,本发明的基于全齿面搜索轮齿接触分析方法主要包括以下步骤:(1)在两齿面参考点处建立局部坐标系,设定齿轮初始位置,对齿轮齿面进行网格划分,给定两齿轮转角分别为和(2)计算网格点在齿轮齿面上投影点的Z坐标,非网格点处齿轮齿面上投影点以二维插值方式可快速获得Z坐标;(3)沿着方向Zb,求解小齿轮任意齿面点到大齿轮齿面的最短距离并以此寻找齿面接近点;(4)根据最短距离和齿面接近点,迭代调整大齿轮转角使得齿面间最短距离小于给定容差ε,此时认为两齿面相切,并求得瞬时接触点;(5)根据瞬时接触点及两齿轮转角求解传动误差。考虑齿轮副滚动印痕检查时的涂层厚度,再对齿轮齿面网格点进行一次距离搜索,寻找接触区域。
上述步骤具体如下:
1、建立随动坐标系S1(坐标轴为x1,y1,z1)和S2(坐标轴为x2,y2,z2)。S1与小齿轮刚性连接,其z1轴与小齿轮轴线重合。同样,S2与大齿轮刚性连接,其z2轴与大齿轮轴重合。
由三坐标测量获得齿轮齿面偏差数据,经过三次B样条拟合得到法向偏差曲面s1(1)(u(1),v(1)),将理论齿面与法向偏差曲面叠加可构造出齿轮实际齿面。
小齿轮齿面Σ(1)在随动坐标系S1下的齿面方程为:
同样,大齿轮齿面Σ(2)在随动坐标系S2下的齿面方程为:
2、选取初始位置。设定初始啮合位置为Oa和Ob重合时的位置,其中初始点Oa和Ob分别位于小齿轮和大齿轮齿宽中点且距离齿顶面为1倍模数处。小齿轮和大齿轮的转角和从初始位置记,取小齿轮初始转角为大齿轮初始转角为即默认两齿轮从初始位置开始求解接触点。坐标系S1、S2及初始点Oa和Ob如图1所示,图中为表示清楚将两齿轮齿间距增大。
3、划分网格。以初始位置小齿轮齿面Σ(1)上初始点Oa为原点,建立与小齿轮刚性连接的随动坐标系Sa(坐标轴为xa,ya,za),坐标系平面xaya与小齿轮齿面Σ(1)相切于原点Oa,轴xa的方向与小齿轮齿向方向相同。在坐标系平面xaya上划分m行n列的网格T(1),网格横向与轴xa平行,纵向与轴ya平行。同理,以Ob为原点,建立与大齿轮刚性连接的随动坐标系Sb(坐标轴为xb,yb,zb),划分网格T(2)。定义T(1)或T(2)上位于每个网格几何中心的点为网格点,其余点为非网格点;由网格点或非网格点在齿轮齿面上投影获得的点为投影点。在后面的距离搜索时将以T(1)上的网格点为索引。小齿轮网格划分示意图如图2所示。
为便于后续距离搜索时的坐标转换,在初始位置处,即时,建立与坐标系Sa重合的固定坐标系Sp(坐标轴为xp,yp,zp),同理建立与坐标系Sb重合的固定坐标系Sg(坐标轴为xg,yg,zg)。由于初始位置处Oa和Ob重合,由随动坐标系Sa、Sb各轴关系可知坐标系Sp绕轴xp旋转180°即为坐标系Sg,如图3所示。
4、计算齿面点Z坐标。计算小齿轮齿面Σ(1)在坐标系Sa下的齿面方程为:
由此可计算两齿轮齿面Σ(1)或Σ(2)上任意点P在相应随动坐标系Sa或Sb下的Z坐标,可通过任意点P在其周围投影点的Z坐标的插值获得。因此在距离搜索时,只有网格点在齿面上的投影点是通过齿面方程求其Z坐标的,非网格点在齿面上的投影点其Z坐标通过二维插值方法获得。小齿轮齿面上投影点Z坐标的计算方法如图4所示,为方便表示将小齿轮齿面Σ(1)的曲率增大。
5、齿面距离计算。设T(1)上的某网格点在小齿轮齿面Σ(1)上的投影点为点点到T(2)平面的投影点记为M,此时点M在网格平面T(2)上通常为非网格点,因此在求点M到大齿轮齿面Σ(2)上的投影点时需要用到前述的插值。为计算点到点的距离,需要将两齿面方程置于同一坐标系下。点在坐标系Sb下的坐标向量可表示为:
其中为坐标系Sa下的矢量方程;Mba是Sa到Sb的坐标转换矩阵,转换方式为:Sa→Sp→Sg→Sb。Ha是小齿轮初始点M(1)到大齿轮初始点M(2)的矢量。ib,jb,kb分别为Sb下三个轴的单位向量。
假设小齿轮初始转动角度为按照传动比计算大齿轮初始旋转角度为计算m行n列的网格T(1)上所有网格点沿za向小齿轮齿面Σ(1)的全部投影点。依据上述齿面距离计算方法,计算小轮齿面Σ(1)上全部投影点到大轮齿面Σ(2)上对应点的距离。寻找其中最小距离,记为Δzba,min。
如果Δzba,min>ε,则说明两齿轮齿面没有互相接触,需要迭代求解瞬时接触点;如果Δzba,min≤ε,则说明两齿轮齿面存在相互接触甚至相互嵌入的情况。ε为给定容差,用于判断两齿面的状态。
6、迭代求解瞬时接触点。如果Δzba,min>ε,则需要做下面调整,否则直接计算传动误差。最小距离处小轮齿面Σ(1)上点为K,大轮齿面Σ(2)上点为Q,记点K和点Q为齿面接近点,如图5所示。
Δzba,min在Q点线速度VQ上的投影设为ΔV,则
式中,QC为点Q到大齿轮轴的最短距离。
之后,保持小齿轮转动角度不变,使大齿轮转动角度为继续搜索齿面接近点,计算新的如此不断迭代调整大齿轮转动角度求解齿面最小距离。当最小距离Δzba,min小于等于给定容差ε时,此时认为齿面已经接触,其对应的齿面接近点为接触点。
当齿面接触上时,根据传动误差的定义计算传动误差。
寻找全部瞬时接触点:为模拟滚检涂层厚度,T(1)上每个网格点在小齿轮齿面Σ(1)上的投影点沿外法线方向移动距离ξ。以新的投影点为索引,再次搜索两齿面间距离Δzba,此次搜索不再迭代调整大齿轮转角。齿面距离Δzba小于等于滚检实验涂层厚度ξ的齿面点都视为接触点。
7、以增量分别沿正向和负向改变小齿轮转角再按照传动比计算大齿轮转角按照步骤5)和6)所述的齿面距离计算方法和瞬时接触点求解方法,计算求解下一个啮合位置处的齿面接触点,并由此计算传动误差、瞬时接触线或接触区域。当小齿轮转角经历[-2π/N1,2π/N1]后则停止计算,其中N1是小齿轮齿数。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明公开的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。
Claims (2)
1.一种基于全齿面搜索的轮齿接触分析方法,其特征在于步骤如下:
步骤1:建立随动坐标系S1和S2,所述坐标系S1的坐标轴为x1,y1,z1,所述坐标系S2的坐标轴为x2,y2,z2,;S1与小齿轮刚性连接,其z1轴与小齿轮轴线重合;同样,S2与大齿轮刚性连接,其z2轴与大齿轮轴重合;
小齿轮齿面∑(1)在随动坐标系S1下的齿面方程为:
同样,大齿轮齿面∑(2)在随动坐标系S2下的齿面方程为:
步骤2:选取初始位置:设定初始啮合位置为Oa和Ob重合时的位置,其中初始点Oa和Ob分别位于小齿轮和大齿轮齿宽中点且距离齿顶面为1倍模数处;小齿轮和大齿轮的转角和从初始位置记,取小齿轮初始转角为大齿轮初始转角为即默认两齿轮从初始位置开始求解接触点;
步骤3:划分网格:以初始位置小齿轮齿面∑(1)上初始点Oa为原点,建立与小齿轮刚性连接的随动坐标系Sa,坐标轴为xa,ya,za,坐标系平面xaya与小齿轮齿面∑(1)相切于原点Oa,轴xa的方向与小齿轮齿向方向相同;在坐标系平面xaya上划分m行n列的网格T(1),网格横向与轴xa平行,纵向与轴ya平行;同理,以Ob为原点,建立与大齿轮刚性连接的随动坐标系Sb,坐标轴为xb,yb,zb,划分网格T(2);定义T(1)或T(2)上位于每个网格几何中心的点为网格点,其余点为非网格点;由网格点或非网格点在齿轮齿面上投影获得的点为投影点;在后面的距离搜索时将以T(1)上的网格点为索引;
为便于后续距离搜索时的坐标转换,在初始位置处,即时,建立与坐标系Sa重合的固定坐标系Sp,坐标轴为xp,yp,zp,同理建立与坐标系Sb重合的固定坐标系Sg,坐标轴为xg,yg,zg;由于初始位置处Oa和Ob重合,由随动坐标系Sa、Sb各轴关系可知坐标系Sp绕轴xp旋转180°即为坐标系Sg;
步骤4:计算齿面点Z坐标
计算小齿轮齿面∑(1)在坐标系Sa下的齿面方程为:
由此可计算两齿轮齿面∑(1)或∑(2)上任意点P在相应随动坐标系Sa或Sb下的Z坐标,可通过任意点P在其周围投影点的Z坐标的插值获得;因此在距离搜索时,只有网格点在齿面上的投影点是通过齿面方程求其Z坐标的,非网格点在齿面上的投影点其Z坐标通过二维插值方法获得;
步骤5:齿面距离计算
设T(1)上的某网格点在小齿轮齿面∑(1)上的投影点为点点到T(2)平面的投影点记为M,此时点M在网格平面T(2)上通常为非网格点,因此在求点M到大齿轮齿面∑(2)上的投影点时需要用到前述的插值;为计算点到点的距离,需要将两齿面方程置于同一坐标系下;点在坐标系Sb下的坐标向量可表示为:
其中为坐标系Sa下的矢量方程;Mba是Sa到Sb的坐标转换矩阵,转换方式为:Sa→Sp→Sg→Sb;Ha是小齿轮初始点M(1)到大齿轮初始点M(2)的矢量;ib,jb,kb分别为Sb下三个轴的单位向量;
假设小齿轮初始转动角度为按照传动比计算大齿轮初始旋转角度为计算m行n列的网格T(1)上所有网格点沿za向小齿轮齿面∑(1)的全部投影点;依据上述齿面距离计算方法,计算小轮齿面∑(1)上全部投影点到大轮齿面∑(2)上对应点的距离;寻找其中最小距离,记为Δzba,min;
如果Δzba,min>ε,则说明两齿轮齿面没有互相接触,需要迭代求解瞬时接触点;如果Δzba,min≤ε,则说明两齿轮齿面存在相互接触甚至相互嵌入的情况;ε为给定容差,用于判断两齿面的状态;
步骤6:迭代求解瞬时接触点
如果Δzba,min>ε,则需要做下面调整,否则直接计算传动误差;最小距离处小轮齿面∑(1)上点为K,大轮齿面∑(2)上点为Q,记点K和点Q为齿面接近点;
Δzba,min在Q点线速度VQ上的投影设为ΔV,则
式中,QC为点Q到大齿轮轴的最短距离;
之后,保持小齿轮转动角度不变,使大齿轮转动角度为继续搜索齿面接近点,计算新的如此不断迭代调整大齿轮转动角度求解齿面最小距离;当最小距离Δzba,min小于等于给定容差ε时,此时认为齿面已经接触,其对应的齿面接近点为接触点;
当齿面接触上时,根据传动误差的定义计算传动误差;
寻找全部瞬时接触点:为模拟滚检涂层厚度,T(1)上每个网格点在小齿轮齿面∑(1)上的投影点沿外法线方向移动距离ξ;以新的投影点为索引,再次搜索两齿面间距离Δzba,此次搜索不再迭代调整大齿轮转角;齿面距离Δzba小于等于滚检实验涂层厚度ξ的齿面点都视为接触点;
2.一种计算机***,其特征在于包括:一个或多个处理器,计算机可读存储介质,用于存储一个或多个程序,其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现权利要求1所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111167325.5A CN113868801B (zh) | 2021-10-04 | 2021-10-04 | 一种基于全齿面搜索的轮齿接触分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111167325.5A CN113868801B (zh) | 2021-10-04 | 2021-10-04 | 一种基于全齿面搜索的轮齿接触分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113868801A CN113868801A (zh) | 2021-12-31 |
CN113868801B true CN113868801B (zh) | 2022-09-16 |
Family
ID=79001786
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111167325.5A Active CN113868801B (zh) | 2021-10-04 | 2021-10-04 | 一种基于全齿面搜索的轮齿接触分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113868801B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114754698B (zh) * | 2022-04-11 | 2023-08-04 | 重庆大学 | 面齿轮齿面测量点规划及在机测量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109992877A (zh) * | 2019-03-28 | 2019-07-09 | 长安大学 | 一种齿轮副几何接触分析的方法 |
CN111666643A (zh) * | 2020-06-16 | 2020-09-15 | 重庆大学 | 一种用于确定复杂齿面接触性能的方法 |
-
2021
- 2021-10-04 CN CN202111167325.5A patent/CN113868801B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109992877A (zh) * | 2019-03-28 | 2019-07-09 | 长安大学 | 一种齿轮副几何接触分析的方法 |
CN111666643A (zh) * | 2020-06-16 | 2020-09-15 | 重庆大学 | 一种用于确定复杂齿面接触性能的方法 |
Non-Patent Citations (1)
Title |
---|
刚支时弧齿锥齿轮齿面啮合迹的确定;刘光磊等;《机械科学与技术》;20010530(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113868801A (zh) | 2021-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhou et al. | A new closed-form calculation of envelope surface for modeling face gears | |
Lin et al. | Numerical tooth contact analysis of a bevel gear set by using measured tooth geometry data | |
CN113868801B (zh) | 一种基于全齿面搜索的轮齿接触分析方法 | |
CN111666643B (zh) | 一种用于确定复杂齿面接触性能的方法 | |
CN108730480B (zh) | 一种rv减速器的摆线轮及其齿廓逆向主动修形方法 | |
Gonzalez-Perez et al. | Reverse engineering of spiral bevel gear drives reconstructed from point clouds | |
CN109992877B (zh) | 一种齿轮副几何接触分析的方法 | |
CN109376448A (zh) | 一种蜗轮滚齿加工齿面三维建模方法、装置及设备 | |
CN103438180B (zh) | 一种定传动比点接触曲线齿锥齿轮的修形方法 | |
CN109241670B (zh) | 一种圆弧齿廓双弧形螺旋齿轮真实齿面创成方法 | |
CN112257192A (zh) | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 | |
CN107617794A (zh) | 一种人字齿内齿圈加工方法 | |
CN110374870B (zh) | 一种双螺杆压缩机螺杆转子间隙布置方法 | |
Liu et al. | Parameterized High‐Precision Finite Element Modelling Method of 3D Helical Gears with Contact Zone Refinement | |
Stein et al. | 3D involute gear evaluation-part II: deviations-basic algorithms for modern software validation | |
Kolivand et al. | Actual tooth contact analysis of straight bevel gears | |
CN106292531A (zh) | 一种计算加工zn1蜗杆盘状成形刀具廓形边界的算法 | |
Zhou et al. | An alignment angle error compensation method of spiral bevel gear tooth surface measurement based on tooth surface matching | |
CN108984889B (zh) | 一种机械式变比转向器齿轮副齿条齿廓数字设计方法 | |
CN111536924A (zh) | 一种滚子包络环面蜗杆齿面检测方法 | |
Liu et al. | Mathematical models and contact simulations of concave beveloid gears | |
CN106767606A (zh) | 基于装配误差的双螺杆泵转子干涉检测方法 | |
CN109492307A (zh) | 一种弧齿锥齿轮齿面载荷接触性能参数的数值计算方法 | |
CN113446960B (zh) | 一种齿面点云理论分布建模方法及测量方法 | |
CN115186493A (zh) | 一种小模数渐开线圆柱齿轮精确建模方法 |
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 |