CN108595765A - 风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法 - Google Patents
风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法 Download PDFInfo
- Publication number
- CN108595765A CN108595765A CN201810255895.1A CN201810255895A CN108595765A CN 108595765 A CN108595765 A CN 108595765A CN 201810255895 A CN201810255895 A CN 201810255895A CN 108595765 A CN108595765 A CN 108595765A
- Authority
- CN
- China
- Prior art keywords
- roller
- coordinate system
- inner ring
- bearing
- raceway
- 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
Links
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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Rolling Contact Bearings (AREA)
Abstract
本发明公开了一种风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,该方法主要是基于坐标向量模拟双列圆锥滚子轴承内圈及滚动体的变形,具体是通过坐标变换来反映滚道与滚子的接触变形以及滚子与挡边的接触变形,并考虑挡边变形、轴承游隙以及滚子修形的影响,再通过建立滚道与滚子静力学平衡方程得到轴承滚道载荷分布,而后在此基础上对轴承寿命进行求解。本发明方法可高效、准确、快捷地计算出轴承载荷分布及轴承寿命,具有实际推广价值。
Description
技术领域
本发明涉及风力发电机组部件强度分析的技术领域,尤其是指一种风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法。
背景技术
双列圆锥滚子轴承是大型风力发电机组传动系统的重要零部件,在半直驱风机、直驱风机的主轴轴承中广泛应用。轴承载荷工况复杂多变,运行过程中轴承内部的受力情况复杂,载荷分布计算困难。
鉴于此,提出一种基于坐标向量模拟双列圆锥轴承内圈及滚子变形的方法,通过坐标转换来反映滚道与滚子的接触变形,在模型中考虑了挡边模型、轴承游隙等的影响,能够更精确求解轴承滚道载荷分布。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提出了一种风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,可高效、准确、快捷地计算轴承载荷分布及轴承寿命。
为实现上述目的,本发明所提供的技术方案为:风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,该方法主要是基于坐标向量模拟双列圆锥滚子轴承内圈及滚动体的变形,具体是通过坐标变换来反映滚道与滚子的接触变形以及滚子与挡边的接触变形,并考虑挡边变形、轴承游隙以及滚子修形的影响,再通过建立滚道与滚子静力学平衡方程得到轴承滚道载荷分布,而后在此基础上对轴承寿命进行求解;其包括以下步骤:
1)通过坐标向量模拟轴承在变形前和变形后的滚子与滚道、滚子与挡边的位置向量,然后再通过坐标转换将接触点向量坐标从滚子坐标系与内圈坐标系转换到系统坐标系下,建立变形协调方程;
2)根据变形协调方程,建立滚子静力学平衡方程和轴承外圈静力学平衡方程,并依据Newton-Raphson方程进行迭代求解滚道、挡边精确的载荷分布;
3)根据轴承滚道的载荷分布求解双列圆锥滚子轴承的基本额定寿命和修正寿命,而在求解过程中采用切片法对轴承滚道进行切片,分别求解每一个切片的寿命,然后积分得到整个轴承的寿命。
所述步骤1)包括以下具体步骤:
1.1)以轴承中心为原点建立系统坐标系Oxyz,以内圈中心为原点建立内圈坐标系Ohxhyhzh,以滚子中心为原点建立滚子坐标系Owhxwhywhzwh,h代表第几列;
1.2)确定变形量
1.2.1)内圈变形
在内圈坐标系下,在外载荷作用下内圈的位移向量为:
δh={δhx,δhy,δhz}T (1)
内圈偏转角向量为:
γh={0,γhy,γhz}T (2)
考虑轴承游隙后内圈中心在3个方向的坐标为:
δ′hx=δhx±Gr/2 (3)
δ′hy=δhyγhz[er(δhx±Gr/2)] (4)
δ′hz=δhyγhy[er(δhx±Gr/2)] (5)
式中:δhx,δhy,δhz分别为不考虑游隙时内圈坐标系中心的位移;γhy,γhz分别为内圈中心绕yh,zh轴的偏转角;δ′hx,δ′hy,δ′hz为考虑径向游隙时内圈坐标系中心的位移;er为内圈坐标系中心到系统坐标系中心的距离;Gr为径向游隙;当h=1时,取+,当h=2时,取-;
1.2.2)滚子变形
在滚子坐标系下,滚子的位移向量为:
Uwh={uwhx,0,uwhz}T (6)
滚子旋转角向量为:
式中:uwhx,uwhz分别为滚子在xwh,zwh方向的位移;为滚子绕ywh方向的旋转角;
1.3)确定变形前接触点位置
1.3.1)内圈坐标系
内、外圈滚道与滚子接触点在滚子坐标系下zwh轴方向的坐标为:
式中:Dpw为滚子组节圆直径;xwhi,xwhe分别为套圈内、外滚道与滚子接触点在滚子坐标系下xwh轴的坐标;ε为滚子半锥角;α,β分别为外、内圈滚道半锥角;γ为轴承的接触角,c为法向间隙,当h=1时,c=-1,当h=2时,c=1;l为滚子长度;
内、外圈滚道与滚子接触点在内圈坐标系下的位置向量为:
式中:e为滚子坐标系中心与内圈坐标系中心的轴向距离;rhix,rhiy,rhiz为滚子与内滚道接触点在内圈坐标系下的坐标;rhex,rhey,rhez为滚子与外滚道接触点在内圈坐标系下3个方向的坐标;ψ为滚子方位角,在系统坐标系下z轴负方向对应的滚子为第1个滚子,其方位角为0,第j个滚子的方位角为Z为单列滚子的数量;
滚子的端部采用球基面,滚子与挡边的接触点在内圈坐标系下的位置向量为:
式中:rhfx,rhfy,rhfz为滚子与挡边接触点在内圈坐标系下的3个方向的坐标;Rs为球基面半径;θ为内圈挡边与轴承中心的连线与x方向的夹角;Dw为滚子大端直径,λ为滚子球基面对应的半角;
1.3.2)滚子坐标系
滚子与内、外滚道接触点在滚子坐标系下的位置向量为:
滚子与挡边接触点在滚子坐标系下的位置向量为:
式中:swhfx,swhfy,swhfz为滚子与挡边接触点在滚子坐标系下3个方向的坐标;
1.4)确定变形后接触点的位置
1.4.1)内圈坐标系
在内圈发生接触变形后,内圈滚道与滚子接触点在内圈坐标系下的位置向量为:
式中:为内圈滚道与滚子接触点在内圈坐标系下3个方向的坐标;当h=1时,取+,当h=2时,取-;
由于外滚道固定,外滚道与滚子接触点在内圈坐标系上的位置向量与变形前一致,即
滚子与挡边接触点在内圈坐标系下的位置向量为:
式中:为滚子与挡边接触点的位置向量在内圈坐标系下3个方向的坐标;
1.4.2)滚子坐标系
在滚子发生变形后滚子与内、外滚道的接触点在滚子坐标系下的位置向量为:
滚子与挡边接触点在滚子坐标系下的位置向量为:
1.5)变形协调方程
通过坐标转换将接触点从滚子坐标系与内圈坐标系转换到系统坐标系下,则
ew=e+er
式中:ew为滚子坐标系中心到系统坐标系中心的距离,当h=1时取+,当h=2时取取-;
滚子与内滚道、外滚道及挡边的法向接触变形为:
式中:nhi nhe 分别为内圈、外圈、挡边的单位法向量;当h=1时取+,当h=2时取取-。
在步骤2)中,所述滚子静力学平衡方程的具体情况如下:
对单个滚子进行受力分析,其平衡方程为:
Qi+Qe+Fc=0 (28)
Ti+Te+Mg=0 (29)
Mg=JrωGωZsinγ (35)
式中:Ki,Ke分别为内、外圈接触刚度;Qi,Qe分别为滚子与内、外圈的接触处载荷;Ti,Te分别为内、外圈滚道作用在滚子上的力矩;Fc为滚子的离心力;Mg为滚子的陀螺力矩;m为单个圆锥滚子质量;RG为滚子质心的转动半径;ωG为滚子公转角速度;Jr为滚子转动惯量;ω为滚子自转角速度;δhi(xwhi,ψ)滚子与内滚道的法向接触变形,δhe(xwhe,ψ)滚子与外滚道的法向接触变形,δhf(ψ)滚子与挡边的法向接触变形;xwhi,xwhe分别为套圈内、外滚道与滚子接触点在滚子坐标系下xwh轴的坐标;ψ为滚子方位角;γ为轴承的接触角;
所述外圈静力学平衡方程的具体情况如下:
对外圈进行受力分析,其平衡方程为:
式中:Fx,Fy,Fz,My,Mz为轴承中心处五个方向的载荷;qeh为外圈与滚子接触点位置单位长度的接触载荷;当h=1时取+,当h=2时取取-;
静力学平衡方程求解过程如下:
根据上述建立的滚子静平衡方程与外圈静力学平衡方程进行求解,其求解过程是:1)输入轴承主要结构参数:l,Rs,Dw,θ,α,β,Z,e,er,Gr;2)输入工况参数:Fx,Fy,Fz,My,Mz;3)给出滚子变形初始值及内圈变形初始值;4)通过公式(28)-(42)进行迭代计算得到滚道载荷分布;
其中,l为滚子长度;Rs为球基面半径;Dw为滚子大端直径;θ为内圈挡边与轴承中心的连线与x方向的夹角;α,β分别为外、内圈滚道半锥角;Z为单列滚子的数量;e为滚子坐标系中心与内圈坐标系中心的轴向距离;er为内圈坐标系中心到系统坐标系中心的距离;Gr为径向游隙。
在步骤3)中,根据静力学平衡方程得到轴承载荷分布情况,进而可求解滚道寿命,具体如下:
将轴承滚道切分为n等份分,每份切片的寿命为:
式中:QCiht为第t个切片内滚道的基本额定载荷;QCeht为第t个切片内滚道的基本额定载荷;QEiht为第t个切片内滚道的当量载荷;QEeht为第t个切片外滚道的当量载荷;
每份切片的基本额定载荷为:
式中:bm为额定寿命修正系数;λs为修正滚子边缘载荷及应力集中而引入的修正系数;Δl为每个切片的长度;Dt为第t个切片的滚子直径;Dpw为滚子组节圆直径;γ为轴承接触角;γ*为轴承结构系数;
内圈即旋转的套圈的滚道每份切片的当量载荷为:
外圈即静止的套圈的滚道每份切片的当量载荷为:
单列轴承内、外圈寿命为:
单列轴承内、外圈综合寿命为:
L10m=L10im+L10em (50)
整个轴承基本额定寿命为:
整个轴承修正寿命为:
La10DRTRB=a1aisoL10DRTRB (52)
上式中:ψ为滚子方位角,qih为内圈单位长度接触载荷;qeh为外圈单位长度接触载荷;a1为可靠性系数;aiso为润滑油系数。
本发明与现有技术相比,具有如下优点与有益效果:
1、传动双列圆锥滚子轴承计算未考虑挡边变形、轴承径向游隙、滚子修形等问题,而本方法考虑了挡边变形、轴承游隙以及滚子修形对轴承滚道载荷分布的影响,得到更加精确的滚道载荷分布。
2、本方法通过坐标变换来反映滚道与滚子的接触变形,滚子与挡边的接触变形,能够更加精确的分析接触变形。
3、本方法对轴承进行切片,通过切变能更加真实的模拟滚子延滚子长度方向的载荷分布。
附图说明
图1为本发明的双列圆锥滚子轴承结构示意图。
图2为本发明的滚子与挡边接触示意图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
本实施例所提供的风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,针对于传统双列圆锥滚子轴承内部载荷计算模型中未考虑挡边变形、轴承径向游隙、滚子修形等问题,提出了一种基于坐标向量模拟双列圆锥滚子轴承内圈及滚动体变形的方法,具体是通过坐标变换来反映滚道与滚子的接触变形以及滚子与挡边的接触变形,并考虑挡边变形、轴承游隙以及滚子修形的影响,再通过建立滚道与滚子静力学平衡方程得到轴承滚道载荷分布,而后在此基础上对轴承寿命进行求解;其具体情况如下:
1)通过坐标向量模拟轴承在变形前和变形后的滚子与滚道、滚子与挡边的位置向量,然后再通过坐标转换将接触点向量坐标从滚子坐标系与内圈坐标系转换到系统坐标系下,建立变形协调方程;其具体包括以下步骤:
1.1)风力发电机用双列圆锥滚子轴承结构如图1所示,以轴承中心为原点建立系统坐标系Oxyz,以内圈中心为原点建立内圈坐标系Ohxhyhzh,以滚子中心为原点建立滚子坐标系Owhxwhywhzwh,h代表第几列;
1.2)确定变形量
1.2.1)内圈变形
在内圈坐标系下,在外载荷作用下内圈的位移向量为:
δh={δhx,δhy,δhz}T (1)
内圈偏转角向量为:
γh={0,γhy,γhz}T (2)
考虑轴承游隙后内圈中心在3个方向的坐标为:
δ′hx=δhx±Gr/2 (3)
δ′hy=δhyγhz[er(δhx±Gr/2)] (4)
δ′hz=δhyγhy[er(δhx±Gr/2)] (5)
式中:δhx,δhy,δhz分别为不考虑游隙时内圈坐标系中心的位移;γhy,γhz分别为内圈中心绕yh,zh轴的偏转角;δ′hx,δ′hy,δ′hz为考虑径向游隙时内圈坐标系中心的位移;er为内圈坐标系中心到系统坐标系中心的距离;Gr为径向游隙;当h=1时,取+,当h=2时,取-;
1.2.2)滚子变形
在滚子坐标系下,滚子的位移向量为:
Uwh={uwhx,0,uwhz}T (6)
滚子旋转角向量为:
式中:uwhx,uwhz分别为滚子在xwh,zwh方向的位移;为滚子绕ywh方向的旋转角;
1.3)确定变形前接触点位置
1.3.1)内圈坐标系
内、外圈滚道与滚子接触点在滚子坐标系下zwh轴方向的坐标为:
式中:Dpw为滚子组节圆直径;xwhi,xwhe分别为套圈内、外滚道与滚子接触点在滚子坐标系下xwh轴的坐标;ε为滚子半锥角;α,β分别为外、内圈滚道半锥角;γ为轴承的接触角,c为法向间隙,当h=1时,c=-1,当h=2时,c=1;l为滚子长度;
内、外圈滚道与滚子接触点在内圈坐标系下的位置向量为:
式中:e为滚子坐标系中心与内圈坐标系中心的轴向距离;rhix,rhiy,rhiz为滚子与内滚道接触点在内圈坐标系下的坐标;rhex,rhey,rhez为滚子与外滚道接触点在内圈坐标系下3个方向的坐标;ψ为滚子方位角,在系统坐标系下z轴负方向对应的滚子为第1个滚子,其方位角为0,第j个滚子的方位角为Z为单列滚子的数量;
滚子与挡边接触示意图如图2所示,滚子的端部采用球基面,滚子与挡边的接触点在内圈坐标系下的位置向量为:
式中:rhfx,rhfy,rhfz为滚子与挡边接触点在内圈坐标系下的3个方向的坐标;Rs为球基面半径;θ为内圈挡边与轴承中心的连线与x方向的夹角;Dw为滚子大端直径,λ为滚子球基面对应的半角;
1.3.2)滚子坐标系
滚子与内、外滚道接触点在滚子坐标系下的位置向量为:
滚子与挡边接触点在滚子坐标系下的位置向量为:
式中:swhfx,swhfy,swhfz为滚子与挡边接触点在滚子坐标系下3个方向的坐标;
1.4)确定变形后接触点的位置
1.4.1)内圈坐标系
在内圈发生接触变形后,内圈滚道与滚子接触点在内圈坐标系下的位置向量为:
式中:为内圈滚道与滚子接触点在内圈坐标系下3个方向的坐标;当h=1时,取+,当h=2时,取-;
由于外滚道固定,外滚道与滚子接触点在内圈坐标系上的位置向量与变形前一致,即
滚子与挡边接触点在内圈坐标系下的位置向量为:
式中:为滚子与挡边接触点的位置向量在内圈坐标系下3个方向的坐标;
1.4.2)滚子坐标系
在滚子发生变形后滚子与内、外滚道的接触点在滚子坐标系下的位置向量为:
滚子与挡边接触点在滚子坐标系下的位置向量为:
1.5)变形协调方程
通过坐标转换将接触点从滚子坐标系与内圈坐标系转换到系统坐标系下,则
ew=e+er
式中:ew为滚子坐标系中心到系统坐标系中心的距离,当h=1时取+,当h=2时取取-;
滚子与内滚道、外滚道及挡边的法向接触变形为:
式中:nhi nhe 分别为内圈、外圈、挡边的单位法向量;当h=1时取+,当h=2时取取-。
2)根据变形协调方程,建立滚子静力学平衡方程和轴承外圈静力学平衡方程,并依据Newton-Raphson方程进行迭代求解滚道、挡边精确的载荷分布;其中,所述滚子静力学平衡方程的具体情况如下:
对单个滚子进行受力分析,其平衡方程为:
Qi+Qe+Fc=0 (28)
Ti+Te+Mg=0 (29)
Mg=JrωGωZsinγ (35)
式中:Ki,Ke分别为内、外圈接触刚度;Qi,Qe分别为滚子与内、外圈的接触处载荷;Ti,Te分别为内、外圈滚道作用在滚子上的力矩;Fc为滚子的离心力;Mg为滚子的陀螺力矩;m为单个圆锥滚子质量;RG为滚子质心的转动半径;ωG为滚子公转角速度;Jr为滚子转动惯量;ω为滚子自转角速度;δhi(xwhi,ψ)滚子与内滚道的法向接触变形,δhe(xwhe,ψ)滚子与外滚道的法向接触变形,δhf(ψ)滚子与挡边的法向接触变形;xwhi,xwhe分别为套圈内、外滚道与滚子接触点在滚子坐标系下xwh轴的坐标;ψ为滚子方位角;γ为轴承的接触角;
所述外圈静力学平衡方程的具体情况如下:
对外圈进行受力分析,其平衡方程为:
式中:Fx,Fy,Fz,My,Mz为轴承中心处五个方向的载荷;qeh为外圈与滚子接触点位置单位长度的接触载荷;当h=1时取+,当h=2时取取-;
静力学平衡方程求解过程如下:
根据上述建立的滚子静平衡方程与外圈静力学平衡方程进行求解,其求解过程是:1)输入轴承主要结构参数:l,Rs,Dw,θ,α,β,Z,e,er,Gr;2)输入工况参数:Fx,Fy,Fz,My,Mz;3)给出滚子变形初始值及内圈变形初始值;4)通过公式(28)-(42)进行迭代计算得到滚道载荷分布;
其中,l为滚子长度;Rs为球基面半径;Dw为滚子大端直径;θ为内圈挡边与轴承中心的连线与x方向的夹角;α,β分别为外、内圈滚道半锥角;Z为单列滚子的数量;e为滚子坐标系中心与内圈坐标系中心的轴向距离;er为内圈坐标系中心到系统坐标系中心的距离;Gr为径向游隙。
3)根据轴承滚道的载荷分布求解双列圆锥滚子轴承的基本额定寿命和修正寿命,而在求解过程中采用切片法对轴承滚道进行切片,分别求解每一个切片的寿命,然后积分得到整个轴承的寿命;其中,根据静力学平衡方程得到轴承载荷分布情况,可求解滚道寿命,具体如下:
将轴承滚道切分为n等份分,每份切片的寿命为:
式中:QCiht为第t个切片内滚道的基本额定载荷;QCeht为第t个切片内滚道的基本额定载荷;QEiht为第t个切片内滚道的当量载荷;QEeht为第t个切片外滚道的当量载荷;
每份切片的基本额定载荷为:
式中:bm为额定寿命修正系数,双列圆锥滚子轴承取551.2;λs为修正滚子边缘载荷及应力集中而引入的修正系数(对于端部切片取0.61,对于中间切片取1);Δl为每个切片的长度;Dt为第t个切片的滚子直径;Dpw为滚子组节圆直径;γ为轴承接触角;γ*为轴承结构系数;
内圈即旋转的套圈的滚道每份切片的当量载荷为:
外圈即静止的套圈的滚道每份切片的当量载荷为:
单列轴承内、外圈寿命为:
单列轴承内、外圈综合寿命为:
L10m=L10im+L10em (50)
整个轴承基本额定寿命为:
整个轴承修正寿命为:
La10DRTRB=a1aisoL10DRTRB (52)
上式中:ψ为滚子方位角,qih为内圈单位长度接触载荷;qeh为外圈单位长度接触载荷;a1为可靠性系数;aiso为润滑油系数。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (4)
1.风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,其特征在于:该方法主要是基于坐标向量模拟双列圆锥滚子轴承内圈及滚动体的变形,具体是通过坐标变换来反映滚道与滚子的接触变形以及滚子与挡边的接触变形,并考虑挡边变形、轴承游隙以及滚子修形的影响,再通过建立滚道与滚子静力学平衡方程得到轴承滚道载荷分布,而后在此基础上对轴承寿命进行求解;其包括以下步骤:
1)通过坐标向量模拟轴承在变形前和变形后的滚子与滚道、滚子与挡边的位置向量,然后再通过坐标转换将接触点向量坐标从滚子坐标系与内圈坐标系转换到系统坐标系下,建立变形协调方程;
2)根据变形协调方程,建立滚子静力学平衡方程和轴承外圈静力学平衡方程,并依据Newton-Raphson方程进行迭代求解滚道、挡边精确的载荷分布;
3)根据轴承滚道的载荷分布求解双列圆锥滚子轴承的基本额定寿命和修正寿命,而在求解过程中采用切片法对轴承滚道进行切片,分别求解每一个切片的寿命,然后积分得到整个轴承的寿命。
2.根据权利要求1所述的风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,其特征在于:所述步骤1)包括以下具体步骤:
1.1)以轴承中心为原点建立系统坐标系Oxyz,以内圈中心为原点建立内圈坐标系Ohxhyhzh,以滚子中心为原点建立滚子坐标系Owhxwhywhzwh,h代表第几列;
1.2)确定变形量
1.2.1)内圈变形
在内圈坐标系下,在外载荷作用下内圈的位移向量为:
δh={δhx,δhy,δhz}T (1)
内圈偏转角向量为:
γh={0,γhy,γhz}T (2)
考虑轴承游隙后内圈中心在3个方向的坐标为:
δ′hx=δhx±Gr/2 (3)
δ′hy=δhyγhz[er(δhx±Gr/2)] (4)
δ′hz=δhyγhy[er(δhx±Gr/2)] (5)
式中:δhx,δhy,δhz分别为不考虑游隙时内圈坐标系中心的位移;γhy,γhz分别为内圈中心绕yh,zh轴的偏转角;δ′hx,δ′hy,δ′hz为考虑径向游隙时内圈坐标系中心的位移;er为内圈坐标系中心到系统坐标系中心的距离;Gr为径向游隙;当h=1时,取+,当h=2时,取-;
1.2.2)滚子变形
在滚子坐标系下,滚子的位移向量为:
Uwh={uwhx,0,uwhz}T (6)
滚子旋转角向量为:
式中:uwhx,uwhz分别为滚子在xwh,zwh方向的位移;为滚子绕ywh方向的旋转角;
1.3)确定变形前接触点位置
1.3.1)内圈坐标系
内、外圈滚道与滚子接触点在滚子坐标系下zwh轴方向的坐标为:
式中:Dpw为滚子组节圆直径;xwhi,xwhe分别为套圈内、外滚道与滚子接触点在滚子坐标系下xwh轴的坐标;ε为滚子半锥角;α,β分别为外、内圈滚道半锥角;γ为轴承的接触角,c为法向间隙,当h=1时,c=-1,当h=2时,c=1;l为滚子长度;
内、外圈滚道与滚子接触点在内圈坐标系下的位置向量为:
式中:e为滚子坐标系中心与内圈坐标系中心的轴向距离;rhix,rhiy,rhiz为滚子与内滚道接触点在内圈坐标系下的坐标;rhex,rhey,rhez为滚子与外滚道接触点在内圈坐标系下3个方向的坐标;ψ为滚子方位角,在系统坐标系下z轴负方向对应的滚子为第1个滚子,其方位角为0,第j个滚子的方位角为Z为单列滚子的数量;
滚子的端部采用球基面,滚子与挡边的接触点在内圈坐标系下的位置向量为:
式中:rhfx,rhfy,rhfz为滚子与挡边接触点在内圈坐标系下的3个方向的坐标;Rs为球基面半径;θ为内圈挡边与轴承中心的连线与x方向的夹角;Dw为滚子大端直径,λ为滚子球基面对应的半角;
1.3.2)滚子坐标系
滚子与内、外滚道接触点在滚子坐标系下的位置向量为:
滚子与挡边接触点在滚子坐标系下的位置向量为:
式中:swhfx,swhfy,swhfz为滚子与挡边接触点在滚子坐标系下3个方向的坐标;
1.4)确定变形后接触点的位置
1.4.1)内圈坐标系
在内圈发生接触变形后,内圈滚道与滚子接触点在内圈坐标系下的位置向量为:
式中:为内圈滚道与滚子接触点在内圈坐标系下3个方向的坐标;当h=1时,取+,当h=2时,取-;
由于外滚道固定,外滚道与滚子接触点在内圈坐标系上的位置向量与变形前一致,即
滚子与挡边接触点在内圈坐标系下的位置向量为:
式中:为滚子与挡边接触点的位置向量在内圈坐标系下3个方向的坐标;
1.4.2)滚子坐标系
在滚子发生变形后滚子与内、外滚道的接触点在滚子坐标系下的位置向量为:
滚子与挡边接触点在滚子坐标系下的位置向量为:
1.5)变形协调方程
通过坐标转换将接触点从滚子坐标系与内圈坐标系转换到系统坐标系下,则
ew=e+er
式中:ew为滚子坐标系中心到系统坐标系中心的距离,当h=1时取+,当h=2时取取-;
滚子与内滚道、外滚道及挡边的法向接触变形为:
式中:nhi nhe 分别为内圈、外圈、挡边的单位法向量;当h=1时取+,当h=2时取取-。
3.根据权利要求1所述的风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,其特征在于:在步骤2)中,所述滚子静力学平衡方程的具体情况如下:
对单个滚子进行受力分析,其平衡方程为:
Qi+Qe+Fc=0 (28)
Ti+Te+Mg=0 (29)
Mg=JrωGωZsinγ (35)
式中:Ki,Ke分别为内、外圈接触刚度;Qi,Qe分别为滚子与内、外圈的接触处载荷;Ti,Te分别为内、外圈滚道作用在滚子上的力矩;Fc为滚子的离心力;Mg为滚子的陀螺力矩;m为单个圆锥滚子质量;RG为滚子质心的转动半径;ωG为滚子公转角速度;Jr为滚子转动惯量;ω为滚子自转角速度;δhi(xwhi,ψ)滚子与内滚道的法向接触变形,δhe(xwhe,ψ)滚子与外滚道的法向接触变形,δhf(ψ)滚子与挡边的法向接触变形;xwhi,xwhe分别为套圈内、外滚道与滚子接触点在滚子坐标系下xwh轴的坐标;ψ为滚子方位角;γ为轴承的接触角;
所述外圈静力学平衡方程的具体情况如下:
对外圈进行受力分析,其平衡方程为:
式中:Fx,Fy,Fz,My,Mz为轴承中心处五个方向的载荷;qeh为外圈与滚子接触点位置单位长度的接触载荷;当h=1时取+,当h=2时取取-;
静力学平衡方程求解过程如下:
根据上述建立的滚子静平衡方程与外圈静力学平衡方程进行求解,其求解过程是:1)输入轴承主要结构参数:l,Rs,Dw,θ,α,β,Z,e,er,Gr;2)输入工况参数:Fx,Fy,Fz,My,Mz;3)给出滚子变形初始值及内圈变形初始值;4)通过公式(28)-(42)进行迭代计算得到滚道载荷分布;
其中,l为滚子长度;Rs为球基面半径;Dw为滚子大端直径;θ为内圈挡边与轴承中心的连线与x方向的夹角;α,β分别为外、内圈滚道半锥角;Z为单列滚子的数量;e为滚子坐标系中心与内圈坐标系中心的轴向距离;er为内圈坐标系中心到系统坐标系中心的距离;Gr为径向游隙。
4.根据权利要求1所述的风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法,其特征在于:在步骤3)中,根据静力学平衡方程得到轴承载荷分布情况,进而可求解滚道寿命,具体如下:
将轴承滚道切分为n等份分,每份切片的寿命为:
式中:QCiht为第t个切片内滚道的基本额定载荷;QCeht为第t个切片内滚道的基本额定载荷;QEiht为第t个切片内滚道的当量载荷;QEeht为第t个切片外滚道的当量载荷;
每份切片的基本额定载荷为:
式中:bm为额定寿命修正系数;λs为修正滚子边缘载荷及应力集中而引入的修正系数;Δl为每个切片的长度;Dt为第t个切片的滚子直径;Dpw为滚子组节圆直径;γ为轴承接触角;γ*为轴承结构系数;
内圈即旋转的套圈的滚道每份切片的当量载荷为:
外圈即静止的套圈的滚道每份切片的当量载荷为:
单列轴承内、外圈寿命为:
单列轴承内、外圈综合寿命为:
L10m=L10im+L10em (50)
整个轴承基本额定寿命为:
整个轴承修正寿命为:
La10DRTRB=a1aisoL10DRTRB (52)
上式中:ψ为滚子方位角,qih为内圈单位长度接触载荷;qeh为外圈单位长度接触载荷;a1为可靠性系数;aiso为润滑油系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810255895.1A CN108595765B (zh) | 2018-03-27 | 2018-03-27 | 风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810255895.1A CN108595765B (zh) | 2018-03-27 | 2018-03-27 | 风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108595765A true CN108595765A (zh) | 2018-09-28 |
CN108595765B CN108595765B (zh) | 2021-12-28 |
Family
ID=63623767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810255895.1A Active CN108595765B (zh) | 2018-03-27 | 2018-03-27 | 风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108595765B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109580223B (zh) * | 2018-12-14 | 2020-07-28 | 河南科技大学 | 一种考虑芯轴变形影响的组合式轴连轴承载荷分析方法 |
CN111914370A (zh) * | 2020-08-07 | 2020-11-10 | 上海理工大学 | 一种电驱动系统高速轴承可靠性试验载荷谱的编制方法 |
WO2022033358A1 (zh) * | 2020-08-12 | 2022-02-17 | 上海理工大学 | 基于多轴随机道路载荷下轮毂轴承服役寿命快速评估方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006105304A (ja) * | 2004-10-06 | 2006-04-20 | Ntn Corp | 車輪用軸受装置 |
CN104636596A (zh) * | 2014-12-26 | 2015-05-20 | 中国北方车辆研究所 | 一种特定载荷下圆柱滚子轴承非对称修形方法 |
CN106649982A (zh) * | 2016-11-08 | 2017-05-10 | 大连工业大学 | 风力发电机大锥角圆锥主轴承摩擦力矩计算方法 |
CN107563081A (zh) * | 2017-09-12 | 2018-01-09 | 大连交通大学 | 高速空心圆柱滚子轴承的参数优化方法 |
CN107563082A (zh) * | 2017-09-12 | 2018-01-09 | 大连交通大学 | 基于圆柱滚子轴承接触变形量和载荷分布的轴承参数优化方法 |
CN107657122A (zh) * | 2017-09-29 | 2018-02-02 | 无锡三立轴承股份有限公司 | 机床主轴圆柱滚子轴承设计方法 |
CN107688716A (zh) * | 2017-09-12 | 2018-02-13 | 大连交通大学 | 基于载荷分布和疲劳寿命的空心圆柱滚子轴承参数优化方法 |
-
2018
- 2018-03-27 CN CN201810255895.1A patent/CN108595765B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006105304A (ja) * | 2004-10-06 | 2006-04-20 | Ntn Corp | 車輪用軸受装置 |
CN104636596A (zh) * | 2014-12-26 | 2015-05-20 | 中国北方车辆研究所 | 一种特定载荷下圆柱滚子轴承非对称修形方法 |
CN106649982A (zh) * | 2016-11-08 | 2017-05-10 | 大连工业大学 | 风力发电机大锥角圆锥主轴承摩擦力矩计算方法 |
CN107563081A (zh) * | 2017-09-12 | 2018-01-09 | 大连交通大学 | 高速空心圆柱滚子轴承的参数优化方法 |
CN107563082A (zh) * | 2017-09-12 | 2018-01-09 | 大连交通大学 | 基于圆柱滚子轴承接触变形量和载荷分布的轴承参数优化方法 |
CN107688716A (zh) * | 2017-09-12 | 2018-02-13 | 大连交通大学 | 基于载荷分布和疲劳寿命的空心圆柱滚子轴承参数优化方法 |
CN107657122A (zh) * | 2017-09-29 | 2018-02-02 | 无锡三立轴承股份有限公司 | 机床主轴圆柱滚子轴承设计方法 |
Non-Patent Citations (2)
Title |
---|
李震等: "双列圆锥滚子轴承拟静力学分析", 《哈尔滨工程大学学报》 * |
艾思源: "滚动轴承温度分布及润滑性能研究", 《工程科技Ⅱ辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109580223B (zh) * | 2018-12-14 | 2020-07-28 | 河南科技大学 | 一种考虑芯轴变形影响的组合式轴连轴承载荷分析方法 |
CN111914370A (zh) * | 2020-08-07 | 2020-11-10 | 上海理工大学 | 一种电驱动系统高速轴承可靠性试验载荷谱的编制方法 |
WO2022033358A1 (zh) * | 2020-08-12 | 2022-02-17 | 上海理工大学 | 基于多轴随机道路载荷下轮毂轴承服役寿命快速评估方法 |
US11835419B2 (en) | 2020-08-12 | 2023-12-05 | University Of Shanghai For Science And Technology | Method for rapid assessment of service life of hub bearing based on multiaxial random road loading |
Also Published As
Publication number | Publication date |
---|---|
CN108595765B (zh) | 2021-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108595765A (zh) | 风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法 | |
US11835419B2 (en) | Method for rapid assessment of service life of hub bearing based on multiaxial random road loading | |
CN107729597A (zh) | 一种主轴轴承滚道校核工具 | |
CN107688716B (zh) | 基于载荷分布和疲劳寿命的空心圆柱滚子轴承参数优化方法 | |
CN105550403B (zh) | 一种风机外购件强度校核工具 | |
CN106649982A (zh) | 风力发电机大锥角圆锥主轴承摩擦力矩计算方法 | |
CN102201033B (zh) | 航空多转子耦合系统动力学的分析方法 | |
CN103256908B (zh) | 一种变桨轴承径向游隙的确定方法 | |
Stammler et al. | Friction torque of wind-turbine pitch bearings–comparison of experimental results with available models | |
CN109753723B (zh) | 一种向心滚动轴承疲劳寿命计算方法 | |
CN110059408A (zh) | 谐波减速器中柔性薄壁轴承疲劳寿命计算方法 | |
CN103256299A (zh) | 一种变桨轴承结构参数的确定方法 | |
Florescu et al. | Research on increasing the performance of wind power plants for sustainable development | |
Chen et al. | Contact stress and deformation of blade bearing in wind turbine | |
CN112861271B (zh) | 风电机组主轴承力学特性的计算方法、系统及介质 | |
CN107563081B (zh) | 高速空心圆柱滚子轴承的参数优化方法 | |
CN107590356A (zh) | 一种风力发电机组主轴轴承自动选型方法及存储设备 | |
Stammler et al. | Load distribution in a roller-type rotor blade bearing | |
CN101476981A (zh) | 一种确定高速滚珠轴承负荷分布的方法 | |
CN113704981A (zh) | 温升过程中高速轴承时变动力学行为分析方法 | |
CN114139294A (zh) | 一种行星滚柱丝杠副传动效率计算方法 | |
CN113468691A (zh) | 风力发电机组的传动轴系的设计方法 | |
CN114357846B (zh) | 一种风电齿轮箱轴承故障识别方法及系统 | |
Tan et al. | Investigation of dynamic characteristics of planetary gear stage in wind turbine considering voltage dip | |
CN109580223A (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 |