-
滚动体与沟道接触为有限长线接触,超出了Hertz理论的范畴。有限长接触问题的基本方程为[9]:
$$ \iint_\varOmega {p\left( {x,y} \right)}{\mathrm{d}}x{\mathrm{d}}y = Q $$ (1) $$ \frac{1}{{\text{π} {E'}}}\iint_\varOmega {\frac{{p\left( {{x'},{y'}} \right){\mathrm{d}}{x'}{\mathrm{d}}{y'}}}{{\sqrt {{{\left( {x - {x'}} \right)}^2} + {{\left( {y - {y'}} \right)}^2}} }}} = \delta - z\left( {x,y} \right) $$ (2) 式中:
Ω ——接触区域;
p ——接触应力(MPa);
Q ——滚动体与沟道的作用力(N);
$ {E^{'}} $ ——当量弹性模量(MPa);
$ \delta $ ——滚动体与沟道接触中心弹性趋近量(mm);
z ——滚动体与沟道接触表面的初始距离(mm)。
针对滚子与滚道线接触问题,可以采用罗继伟[21]的方法,如图2所示,将接触区域沿滚动体轴向(y轴)分割为n个单元,假定切片区域j内的接触应力pj沿母线方向符合均匀分布,沿x轴为赫兹分布:
$$ {p_j}\left( {x,y} \right) = {p_{{\mathrm{o}}j}}\sqrt {1 - {{\left( {\frac{x}{{a{}_j}}} \right)}^2}} $$ (3) 式中:
poj ——接触切片j中心处的最大接触应力(MPa);
2aj、2hj ——接触切片j的长边和短边(mm)。
切片j产生的接触载荷为:
$$ {q_j} = 2{h_j}\int_{ - {a_j}}^{{a_j}} {{p_{{\mathrm{o}}j}}} \sqrt {1 - {{\left( {x/{a_j}} \right)}^2}} {\mathrm{d}}x = \text{π}{a_j}{h_j}{p_{{\mathrm{o}}j}} $$ (4) 柔度系数为:
$$ {D_{ij}} = \int\limits_{ - {a_j}}^{{a_j}} {\int\limits_{{y_j} - {h_j}}^{{y_j} + {h_j}} {\frac{{\sqrt {1 - {{\left( {{x'}/{a_j}} \right)}^2}} }}{{\sqrt {{x'}^2 + {{\left( {{y_i} - {y_j} - {y'}} \right)}^2}} }}} } {\mathrm{d}}{x'}{\mathrm{d}}{y'} $$ (5) 将式(1)、式(2)进行离散化处,可以得到:
$$ Q = \text{π} \sum\limits_{j = 1}^n {a{}_j{h_j}} {p_{{\mathrm{o}}j}} $$ (6) $$ \frac{1}{{\text{π} {E'}}}\sum\limits_{j = 1}^n {{D_{ij}}} {p_{{\mathrm{o}}j}} = \delta - {Z_i}\left( {{y_i}} \right)\;\; \left( {{{i = 1,2,}} \cdots n} \right)$$ (7) 式中:
Zi(yi)——滚动体第i个切片处与沟道表面的初始间距(mm)。
式(6)和式(7)构成了滚子与滚道接触力学模型的方程组,运用数值解法即可得到各切片单元接触应力poj。利用式(3)和式(4)进一步计算得到整个接触区域的应力分布和每个切片的接触载荷。
-
风电主轴承运行工况复杂多变,在进行轴承滚道疲劳寿命计算时,需要考虑润滑、污染、倾斜、轴承工作游隙等因素的影响。所以采用ISO/TS
16281 标准[22]介绍的方法,进行轴承修正参考额定寿命的计算。对应于90%可靠度的基本参考额定寿命为:
$$ {L_{10{\mathrm{r}}}} = {\left\{ {\sum\limits_{k = 1}^{{n_{\mathrm{s}}}} {\left[ {{{\left( {\dfrac{{{q_{{{k}}{\mathrm{ci}}}}}}{{{q_{{{k}}{\mathrm{ei}}}}}}} \right)}^{ - 9/2}} + {{\left( {\dfrac{{{q_{{{k}}{\mathrm{ce}}}}}}{{{q_{{{k}}{\mathrm{ee}}}}}}} \right)}^{ - 9/2}}} \right]} } \right\}^{ - 8/9}} $$ (8) 在基本参考额定寿命的基础上,进行修正得到更符合轴承实际运行的修正参考额定寿命:
$$ \begin{split} &{L_{n{\mathrm{mr}}}} = \\&{a_1}{\left( {\sum\limits_{k = 1}^{{n_{\mathrm{s}}}} {\left\{ {{{\left[ {{a_{{\mathrm{ISO}}}}\left( {\frac{{{e_{\mathrm{C}}}{C_{{\mathrm{ur}}}}}}{{{P_{k{\mathrm{s}}}}}},\kappa } \right)} \right]}^{ - 9/8}}\left[ {{{\left( {\frac{{{q_{k{\mathrm{ci}}}}}}{{{q_{k{\mathrm{ei}}}}}}} \right)}^{ - 9/2}} + {{\left( {\frac{{{q_{k{\mathrm{ce}}}}}}{{{q_{k{\mathrm{ee}}}}}}} \right)}^{ - 9/2}}} \right]} \right\}} } \right)^{ - 8/9}} \end{split} $$ (9) 式中:
$ {a}_{1} $ ——可靠度修正系数;
$ {a}_{\mathrm{I}\mathrm{S}\mathrm{O}} $ ——修正系数;
$ {e}_{\mathrm{C}} $ ——污染系数;
$ {C}_{\mathrm{u}\mathrm{r}} $ ——径向载荷疲劳极限(N);
$ \kappa $ ——润滑剂黏度比;
$ {P}_{{k}\mathrm{s}} $ ——轴承第k个切片的当量动载荷(N);
$ {q}_{{k}\mathrm{c}\mathrm{i}} $、$ {q}_{{k}\mathrm{c}\mathrm{e}} $——内圈和外圈某一轴承切片处的基本额定动载荷(N);
$ {q}_{{k}\mathrm{e}\mathrm{i}} $、$ {q}_{{k}\mathrm{e}\mathrm{e}} $——内圈和外圈某一切片处的当量动载荷(N)。
对于多工况载荷谱,可以分别计算每个工况下的轴承滚道寿命,在运用GL规范介绍方法计算合并修正参考额定寿命为:
$$ {L_{10{\mathrm{mr}}}} = \dfrac{{\displaystyle \sum {{q_i}} }}{{\displaystyle \sum {\dfrac{{{q_i}}}{{{L_{10{\mathrm{mr}}i}}}}} }} $$ (10) 式中:
qi ——第i个载荷等级的分时(h);
L10mri ——第i个载荷等级的轴承修正参考额定寿命(h)。
-
本文利用精益传动设计软件,进行建模和分析计算。精益传动设计软件可以进行传动系统二维和三维建模,能涵盖整个传动产品设计流程。吸收了国际设计标准规范,不仅可以实现齿轮、轴承、轴等传动部件的分析计算,还可以进行系统NVH、变形、模态以及转子动力学等关键技术的动力学分析计算。支持多任务、多数据快速并行算法,是一款多物理场耦合的大型传动系统设计分析软件。主要面向汽车、风电、船舶、燃机等行业的国产传动系统仿真计算工业软件。
-
以某风电机组主轴系为例,利用精益传动设计软件,建立了轮毂、主轴、上风向轴承和下风向轴承一体化的分析模型,如图3所示。
其中轮毂和主轴设置好材料属性柔性化处理后导入精益传动设计软件,上下风向轴承可以根据轴承参数在软件中自定义设定添加,材料属性和轴承参数分别如表1和表2所示。轮毂中心施加风机叶轮中心载荷,软件中设置上下风向单列圆锥滚子轴承总轴向游隙,进行仿真分析。
表 1 材料属性
Table 1. Material properties
名称 弹性模量/GPa 泊松比 密度/(kg·m−3) 主轴 200 0.30 7 850 轮毂 169 0.28 7 300 表 2 轴承参数
Table 2. Bearing parameters
参数 上风向轴承 下风向轴承 内径/mm 1 500 1 350 外径/mm 1 965 1 740 宽度/mm 220 190 重量/kg 1 751 1 153 -
主轴转速为8.8 r/min,采用油润滑,运动粘度为320 mm2/s,应用环境为正常清洁度。通过分析模型的计算,得到了不同轴向游隙下,上下风向轴承的合并修正参考额定寿命、承载区域以及最大接触应力的变化趋势。图4~图6横坐标负数为正游隙,表示轴承间隙;正数为负游隙,表示轴承预紧。
风电机组疲劳载荷分析形式主要有3种:实时的时间序列、载荷持续分布(LDD)和等效疲劳载荷[23]。主轴轴承滚道疲劳寿命一般采用载荷持续分布进行计算[24],载荷包含力矩和力在各方向上的分量,即Mx,My,Mz,Fx,Fy,Fz,以及时间和转速。
图4为在轮毂中心LDD疲劳载荷下,不同轴向游隙下,轴承滚道合并修正参考额定寿命的变化趋势。随着轴向游隙的减小,上风向轴承滚道寿命先增加后减小,下风向轴承滚道寿命逐渐减小;当轴向游隙为−0.6 mm,上风向轴承和下风向轴承滚道寿命接近;当轴向游隙值小于−0.74 mm时,下风向轴承滚道寿命低于要求设计寿命7 300 d;当轴向游隙为−0.36 mm时,上风向轴承滚道寿命达到最大。因此,上风向轴承和下风向轴承同时满足大于要求设计寿命和轴向工作游隙为负游隙的条件下,合适轴向工作游隙范围为−0.74 mm~0,配对轴承寿命最优轴向工作游隙为−0.36 mm。
图5为不同轴向游隙下,上风向轴承和下风向轴承承载区的变化。由图可知,随着轴向游隙的减小,上风向轴承和下风向轴承承载区逐渐增大。下风向轴承承载区始终大于180°,有超过一半滚子个数在承载区内。轴向游隙大于−0.3 mm时,上风向轴承承载区小于180°,小于一半滚子个数位于承载区。
为了减少单个滚子的接触应力,风电轴承往往要求至少有一半的滚子位于承载区。上风向轴承和下风向轴承同时满足至少有一半滚子位于承载区的轴向工作游隙范围是−1~−0.3 mm。
综上所述,同时考虑轴向游隙与寿命、承载区的关系,可以得到配对的上风向轴承和下风向轴承理想的轴向工作游隙范围为−0.74~−0.3 mm。
图6为轮毂中心极限载荷下,不同轴向游隙下,上风向轴承和下风向轴承最大接触应力变化。由图可知,随着轴向游隙减小,上风向轴承最大接触应力逐渐减小,下风向轴承最大接触应力逐渐增大;轴向工作游隙−0.74~−0.3 mm范围内,上下风向轴承最大接触应力均小于最大接触应力要求值2 500 MPa。因此,上述轴向工作游隙范围内,上下风向轴承满足极限的设计要求,进一步证明上述轴向工作游隙是合理的。
Analysis on Operating Clearance of Matched Single Row Tapered Roller Bearings for Wind Turbine Main Shaft
-
摘要:
目的 单列圆锥滚子轴承双支撑的主轴系结构在风电机组中得到了广泛的应用,确定合适的轴承工作游隙,对轴承和主轴系的可靠运行至关重要。 方法 以某中速永磁风电机组主轴系为研究对象,建立了轮毂、主轴、上风向轴承和下风向轴承一体化的动力学分析模型;计算得到了不同轴向游隙下,轴承滚道寿命、承载区域以及最大接触应力的变化趋势。 结果 确定了配对单列圆锥滚子轴承理想的轴向工作游隙范围为−0.74~−0.3 mm。 结论 研究结果对轴承和主轴系的可靠运行具有重大意义,为配对单列圆锥滚子轴承轴向工作游隙的确认提供了方法和参考。 Abstract:Introduction The main shafting structure with double supports of single row tapered roller bearings has been widely used in wind turbines, and determining the appropriate bearing operating clearance is crucial for the reliable operation of bearings and the main shafting. Method A dynamics analysis model was established focusing on the main shafting of a medium speed permanent magnet wind turbine, incorporating the hub, main shaft, and both upwind and downwind bearing. The variation trends of bearing raceway life, load area, and maximum contact stress under different axial clearances were calculated. Result The ideal axial operating clearance range of matched single row tapered roller bearings has been identified as −0.74 ~ −0.3 mm. Conclusion The research results have significant implications for the reliable operation of bearings and main shafting, providing a method and reference for determining the axial operating clearance of matched single row tapered roller bearings. -
表 1 材料属性
Tab. 1. Material properties
名称 弹性模量/GPa 泊松比 密度/(kg·m−3) 主轴 200 0.30 7 850 轮毂 169 0.28 7 300 表 2 轴承参数
Tab. 2. Bearing parameters
参数 上风向轴承 下风向轴承 内径/mm 1 500 1 350 外径/mm 1 965 1 740 宽度/mm 220 190 重量/kg 1 751 1 153 -
[1] 陈嘉豪, 尹梓炜, 朱容宽, 等. 漂浮式海上风机平台阻尼结构设计与研究 [J]. 南方能源建设, 2021, 8(3): 18-25. DOI: 10.16516/j.gedi.issn2095-8676.2021.03.003. CHEN J H, YIN Z W, ZHU R K, et al. Design and research on the foundation damping structure of semi-submersible floating offshore wind turbines [J]. Southern energy construction, 2021, 8(3): 18-25. DOI: 10.16516/j.gedi.issn2095-8676.2021.03.003. [2] 施伟, 郑侃, 任年鑫. 南海海况下半潜浮式风机在故障工况下的动力学响应分析 [J]. 南方能源建设, 2018, 5(4): 12-20. DOI: 10.16516/j.gedi.issn2095-8676.2018.04.002. SHI W, ZHENG K, REN N X. Dynamic analysis of semi-type floating offshore wind turbine with failure conditions under metocean conditions in South China Sea [J]. Southern energy construction, 2018, 5(4): 12-20. DOI: 10.16516/j.gedi.issn2095-8676.2018.04.002. [3] 王慧. 大功率海上风电机组关键部位轴承选型与技术分析 [J]. 风能, 2013(12): 92-95. DOI: 10.3969/j.issn.1674-9219.2013.12.043. WANG H. Analysis of bearings selections and technical features for key positions in large capacity offshore wind turbines [J]. Wind energy, 2013(12): 92-95. DOI: 10.3969/j.issn.1674-9219.2013.12.043. [4] 张利慧, 李殊瑶, 李晓波, 等. 基于时频域特征参数的风电机组滚动轴承故障诊断方法及应用 [J]. 内蒙古电力技术, 2024, 42(3): 13-19. DOI: 10.19929/j.cnki.nmgdljs.2024.0035. ZHANG L H, LI S Y, LI X B, et al. Diagnostic method of wind turbine rolling bearing failure based on time-frequency domain characteristic parameters and its application [J]. Inner Mongolia electric power, 2024, 42(3): 13-19. DOI: 10.19929/j.cnki.nmgdljs.2024.0035. [5] 陈方述, 阳雪兵, 石峰. 永磁直驱风力发电机组弯头形机舱优化分析 [J]. 南方能源建设, 2023, 10(1): 154-159. DOI: 10.16516/j.gedi.issn2095-8676.2023.01.020. CHEN F S, YANG X B, SHI F. Optimization analysis of permanent-magnet direct-drive wind turbine elbow shape nacelle [J]. Southern energy construction, 2023, 10(1): 154-159. DOI: 10.16516/j.gedi.issn2095-8676.2023.01.020. [6] 申屠东华, 何先照, 卢江跃, 等. 大型风电机组传动系统支撑方式特性分析 [J]. 装备制造技术, 2019(1): 129-133. DOI: 10.3969/j.issn.1672-545X.2019.01.034. SHENTU D H, HE X Z, LU J Y, et al. Analysis of characteristics of transmission system support of large wind turbine [J]. Equipment manufacturing technology, 2019(1): 129-133. DOI: 10.3969/j.issn.1672-545X.2019.01.034. [7] 杜静, 江术新, 黄强, 等. 风力发电机主轴轴承载荷分布及接触应力分析 [J]. 机械设计, 2014, 31(4): 64-67. DOI: 10.13841/j.cnki.jxsj.2014.04.014. DU J, JIANG S X, HUANG Q, et al. Analysis of load distribution and contact stress of wind turbine main shaft bearing [J]. Journal of machine design, 2014, 31(4): 64-67. DOI: 10.13841/j.cnki.jxsj.2014.04.014. [8] 付大鹏, 褚加瑞. 基于Romax的风电机组主轴轴承接触应力分析 [J]. 东北电力大学学报, 2018, 38(6): 62-67. DOI: 10.3969/j.issn.1005-2992.2018.06.011. FU D P, CHU J R. Analysis on contact stress for the wind turbines main shaft bearing based on Romax [J]. Journal of Northeast Electric Power University, 2018, 38(6): 62-67. DOI: 10.3969/j.issn.1005-2992.2018.06.011. [9] 杨剑飞, 李秀珍, 欧惠宇, 等. 基于有限长接触的风电机组主轴承接触应力分析 [J]. 轴承, 2020(2): 1-6. DOI: 10.19533/j.issn1000-3762.2020.02.001. YANG J F, LI X Z, OU H Y, et al. Analysis on contact stress of main bearings for wind turbines based on finite length contact [J]. Bearing, 2020(2): 1-6. DOI: 10.19533/j.issn1000-3762.2020.02.001. [10] 牛宝禛, 李伦, 李济顺, 等. 基于多因素的风电主轴轴承疲劳寿命分析 [J]. 轴承, 2022(8): 9-14,18. DOI: 10.19533/j.issn1000-3762.2022.08.002. NIU B Z, LI L, LI J S, et al. Analysis on fatigue life of main shaft bearings in wind turbines based on multiple factors [J]. Bearing, 2022(8): 9-14,18. DOI: 10.19533/j.issn1000-3762.2022.08.002. [11] 姜宏伟. 过盈配合下风电机组主轴轴承寿命研究 [J]. 重型机械, 2021(1): 73-76. DOI: 10.3969/j.issn.1001-196X.2021.01.013. JIANG H W. Research on wind turbine main bearing life considering the influence of interference fit [J]. Heavy machinery, 2021(1): 73-76. DOI: 10.3969/j.issn.1001-196X.2021.01.013. [12] 李云峰, 范雨晴, 王高峰, 等. 风电主轴双支承圆锥滚子轴承疲劳寿命计算 [J]. 轴承, 2023(6): 26-35. DOI: 10.19533/j.issn1000-3762.2023.06.004. LI Y F, FAN Y Q, WANG G F, et al. Fatigue life calculation of tapered roller bearings for doubly supporting wind turbine main shaft [J]. Bearing, 2023(6): 26-35. DOI: 10.19533/j.issn1000-3762.2023.06.004. [13] 董营, 吕超, 赵登利, 等. 风机主轴用单列圆锥滚子轴承安装游隙控制工艺 [J]. 轴承, 2019(8): 19-21,29. DOI: 10.19533/j.issn1000-3762.2019.08.005. DONG Y, LÜ C, ZHAO D L, et al. Control process for assembly clearance of single row tapered roller bearings for wind turbine mainshaft [J]. Bearing, 2019(8): 19-21,29. DOI: 10.19533/j.issn1000-3762.2019.08.005. [14] 肖文, 王忠强, 裴世源, 等. 大兆瓦风电主轴双列圆锥滚子轴承的承载接触机理 [J]. 机械设计与制造, 2021(6): 90-94. DOI: 10.3969/j.issn.1001-3997.2021.06.021. XIAO W, WANG Z Q, PEI S Y, et al. Contact mechanism of double-row tapered roller bearings of main shaft of wind turbine under loading condition [J]. Machinery design & manufacture, 2021(6): 90-94. DOI: 10.3969/j.issn.1001-3997.2021.06.021. [15] ZHENG J Y, JI J C, YIN S, et al. Fatigue life analysis of double-row tapered roller bearing in a modern wind turbine under oscillating external load and speed [J]. Proceedings of the institution of mechanical engineers, part C: journal of mechanical engineering science, 2020, 234(15): 3116-3130. DOI: 10.1177/0954406220911966. [16] 章滔, 阳雪兵, 何录忠, 等. 风电机组回转支承滚子修形曲线设计方法研究 [J]. 机械强度, 2021, 43(1): 183-190. DOI: 10.16579/j.issn.1001.9669.2021.01.027. ZHANG T, YANG X B, HE L Z, et al. Research on the design method of roller modification curve for wind turbine [J]. Journal of mechanical strength, 2021, 43(1): 183-190. DOI: 10.16579/j.issn.1001.9669.2021.01.027. [17] 吕润楠, 郝旭, 于长鑫, 等. 双列圆锥滚子轴承滚道表面波纹度对轴承振动特性的影响研究 [J]. 振动与冲击, 2022, 41(20): 126-132. DOI: 10.13465/j.cnki.jvs.2022.20.016. LÜ R N, HAO X, YU C X, et al. Effects of vibration characteristics of a double row tapered roller bearing with raceway surface waviness [J]. Journal of vibration and shock, 2022, 41(20): 126-132. DOI: 10.13465/j.cnki.jvs.2022.20.016. [18] 孟高强. 浅析配对单列圆锥滚子轴承的应用 [J]. 装备制造技术, 2016(12): 176-178. DOI: 10.3969/j.issn.1672-545X.2016.12.058. MENG G Q. The application of matching single row tapered roller bearing is analyzed [J]. Equipment manufacturing technology, 2016(12): 176-178. DOI: 10.3969/j.issn.1672-545X.2016.12.058. [19] 陈建斌. 圆锥滚子轴承轴向游隙对接触应力影响分析 [D]. 南昌: 华东交通大学, 2018. CHEN J B. Influence of axial clearance of tapered roller bearings on contact stress [D]. Nanchang: East China Jiaotong University, 2018. [20] 彭智. 圆锥滚子轴承接触力学特性分析及优化设计 [D]. 西安: 长安大学, 2021. DOI: 10.26976/d.cnki.gchau.2021.001840. PENG Z. Tapered roller bearing contact mechanics characteristics analysis and optimization design [D]. Xi'an: Chang'an University, 2021. DOI: 10.26976/d.cnki.gchau.2021.001840. [21] 罗继伟, 罗天宇. 滚动轴承分析计算与应用 [M]. 北京: 机械工业出版社, 2009. LUO J W, LUO T Y. Analysis and application of rolling bearing [M]. Beijing: China Machine Press, 2009. [22] ISO. Rolling bearings-methods for calculating the modified reference rating life for universally loaded bearings: ISO/TS 16281: 2008 [S]. Geneva: International Organization for Standardization, 2008. [23] 许移庆. 风电机组疲劳载荷影响因素分析(三) [J]. 风能, 2013(10): 82-86. DOI: 10.3969/j.issn.1674-9219.2013.10.016. XU Y Q. Analysis of impact factors on wind turbine fatigue loads (Ⅲ) [J]. Wind energy, 2013(10): 82-86. DOI: 10.3969/j.issn.1674-9219.2013.10.016. [24] 李英昌, 刘桂然, 员一泽. 风力发电机主轴轴承加速寿命试验研究 [J]. 机械, 2023, 50(3): 14-20. DOI: 10.3969/j.issn.1006-0316.2023.03.003. LI Y C, LIU G R, YUAN Y Z. Study on accelerated life test of main shaft bearing of wind turbine [J]. Machinery, 2023, 50(3): 14-20. DOI: 10.3969/j.issn.1006-0316.2023.03.003.