CN107451376B - 一种偏航变桨轴承滚道校核工具 - Google Patents
一种偏航变桨轴承滚道校核工具 Download PDFInfo
- Publication number
- CN107451376B CN107451376B CN201710753617.4A CN201710753617A CN107451376B CN 107451376 B CN107451376 B CN 107451376B CN 201710753617 A CN201710753617 A CN 201710753617A CN 107451376 B CN107451376 B CN 107451376B
- Authority
- CN
- China
- Prior art keywords
- load
- bearing
- rolling element
- row
- contact
- 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
- 230000003068 static effect Effects 0.000 claims abstract description 69
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000011156 evaluation Methods 0.000 claims abstract description 29
- 238000001228 spectrum Methods 0.000 claims abstract description 20
- 238000005096 rolling process Methods 0.000 claims description 114
- 238000004364 calculation method Methods 0.000 claims description 45
- 229910000831 Steel Inorganic materials 0.000 claims description 26
- 239000010959 steel Substances 0.000 claims description 26
- 238000005452 bending Methods 0.000 claims description 12
- 238000002940 Newton-Raphson method Methods 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000008676 import Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000009472 formulation Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims 3
- 238000012795 verification Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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/04—Ageing analysis or optimisation against ageing
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Geometry (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Computational Mathematics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Primary Health Care (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Rolling Contact Bearings (AREA)
- Coloring Foods And Improving Nutritive Qualities (AREA)
Abstract
本发明公开了一种偏航变桨轴承滚道校核工具,直接利用LDD载荷谱进行疲劳寿命评估,对时序载荷进等效处理,形式与LDD载荷谱形式一致;工具可根据经验公式或静力学方程法求解轴承静强度,也可通过积分公式法或静力学方程法求解轴承疲劳寿命,同时可求解最大动态剪切应力及硬化层深度。工具可以根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命;工具可以根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命。本发明工具可非常有效缩短整机开始设计周期,提高风机可靠性。
Description
技术领域
本发明涉及风力发电机组部件强度分析领域,尤其是指一种偏航变桨轴承滚道校核工具。
背景技术
业内习知,偏航变桨轴承作为风力发电机组的回转支承部件,其结构的可靠性决定风机运行的情况,偏航变桨轴承的结构形式一般为单排四点接触球轴承、双排四点接触球轴承,目前每一家轴承供应商的计算方法、评价指标不一,这就为后期的载荷评估带来麻烦,因此有必要针对偏航变桨轴承建立一个集成不同算法、不同评价标准的工具。
发明内容
本发明的目的在于克服现有偏航变桨轴承滚道校核计算过程中的不足,提供了一种高效、快捷、系统的偏航变桨轴承滚道校核工具,该工具集成了不同的计算方法及评价指标,能够满足不同供应商计算轴承的要求,可非常有效缩短整机开始设计周期,更好服务于风力发电机设计、认证和评估,提高风机发电机组设计的可靠性,节约成本。
为实现上述目的,本发明所提供的技术方案为:一种偏航变桨轴承滚道校核工具,所述工具为基于Matlab开发的偏航变桨轴承滚道校核软件,主要用于风力发电机偏航、变桨轴承初期选型、认证以及后期的结构强度评估;
所述工具直接利用LDD载荷谱进行疲劳寿命评估,同时能够对时序载荷进等效处理,将每一个工况等效为一组载荷,等效的载荷谱的形式与LDD载荷谱形式一致;
偏航变桨轴承通常采用单排四点接触球轴承或双排四点接触球轴承的结构形式,所述工具能够根据经验公式或静力学方程法求解轴承静强度,也能够通过积分公式法或静力学方程法求解轴承疲劳寿命,同时求解最大动态剪切应力及硬化层深度;
所述工具针对不同轴承供应商的轴承寿命提出以下寿命评估方法:
对于单排四点接触球轴承,根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命;
对于双排四点接触球轴承,根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命;
其中,按转数进行评估轴承寿命时应以载荷谱中实际的转数为评价指标;按时间进行评估轴承寿命时陆上风电应以20年的运行时间为评价指标,海上风电应以25年的运行时间为评价指标。
所述时序载荷的处理情况如下:
时序载荷:
对每一个时序载荷工况的力、弯矩、速度通过下面方程1.1-1.3处理得到等效载荷工况表;
式1.1-1.3中,F为轴向力、径向力等效载荷,Fi为时序载荷中第i步的力,M为合弯矩等效载荷,Mi为时序载荷中第i步的弯矩,n为等效转速,ni为时序载荷中第i步的转速;
将每一个疲劳工况的发生时间按下述方法求解:
对于发生次数的工况:
T20/25=L×T1×0.05×b (1.4)
对于发生时间的工况:
T20/25=6×L×T1×0.05×b (1.5)
式1.4-1.5中,L为风机的寿命要求,对于海上工况为25年,陆上工况20年,T1为每个工况发生的次数或小时数,b为每个工况发生的步数,0.05表示步长,T20/25为20年或25年每个疲劳工况发生的时间;
对于LDD载荷则可直接导入工具进行寿命计算;
所述静力学方程的情况如下:
1)四点接触球轴承的受载后几何形变如下:
坐标系采用GL规范中叶根坐标系,假设轴承周边结构具有足够的刚性,外滚道沟渠率中心位置坐标系为xo,yo,内滚道沟渠率中心位置坐标系为xi,yi,轴承外圈固定,内圈相对于外圈产生轴向位移δa,径向位移δr,扭角θ;坐标系位置也由xi,yi转化xi',yi'再到xi”,yi”最后到xi”',yi”',沟渠率中心距离也由最初的MN变为MN”';
对于双排滚子轴承,由排距d与转角θ对内外滚道产生0.5dθ变形;
对于单排轴承,不存在排距的影响;
2)双排四点接触球轴承的力学模型与负载分布计算
四个接触对在位置角ψ处的接触力分别表示为Q1ψ、Q2ψ、Q3ψ、Q4ψ,根据Hertz点接触理论,法向接触载荷Qiψ和接触变形δiψ的关系如下:
式中,在位置角ψj处,内圈受到轴向载荷Fa、径向载荷Fr、倾覆力矩M以及钢球对内滚道的接触载荷Qiψ的作用;Kn为滚动体与内外圈总的负荷变形常数。
以内圈与滚动体为研究对象,内圈在外部载荷和所有滚动体载荷的作用下处于平衡状态,内圈的力学平衡方程为:
其中,α1ψ、α2ψ、α3ψ、α4ψ分别为受载后四个接触对在位置角ψ处的接触角;
上述方程通过Newton Raphson方法进行迭代求解,得到滚道精确的载荷分布;
3)单排四点接触球轴承的力学模型与负载分布计算
单排四点接触球轴承的静力学方程求解原理与上述双排四点接触球轴承的计算原理相同;
两个接触对在位置角ψ处的接触力分别表示为Q1ψ、Q2ψ,受载后两个接触对在位置角ψ处的接触角为α1ψ、α2ψ;
故平衡方程为:
所述静强度计算的情况具体如下:
1)积分公式法的静强度根据下列公式求解
单排轴承:
双排轴承:
根据赫兹接触理论得:
式3.1-3.4中,Qmax为滚动体最大载荷,Fr为径向载荷,Fa为轴向载荷,M为合弯矩,Z为钢球总数,α为初始接触角,dm为节圆直径,a为接触椭圆短半轴,b为接触椭圆的长半轴,Smin为最小安全系数,[σ]为许用接触应力,σmax为最大接触应力;
2)静力学方程的静强度求解
根据方程2.2-2.7,通过Newton Raphson方法即可求出最大滚动体的载荷Qmax,其余的计算过程与积分公式法一致;
所述疲劳寿命的计算情况如下:
1)滚道的基本额定滚动体负荷:
其中,内滚道取上符号,外滚道取下符号;
式4.1-4.3中,Qci为内滚道的基本额定动负荷,Qce为外滚道的基本额定动负荷,fi为内滚道沟曲率半径系数,fo为外滚道沟曲率半径系数,γ为无量纲几何参数,为Dw为滚子直径,dm为滚道中心径,Z为钢球总个数;
2)积分公式求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
2.1)当量滚动体负荷
相对负荷方向旋转套圈的当量滚动体负荷积分为:
相对负荷方向静止套圈的当量滚动体负荷积分为:
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev1=Qev2=QmaxJ1(ε) (4.6)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu1=Qeu2=QmaxJ2(ε) (4.7)
式4.4-4.7中,ε为载荷系数,ψ为位置角,Qmax为滚动体最大载荷,根据式3.1计算;
2.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev1为:
相对负荷静止套圈主推力滚道的额定寿命L10eu1为:
相对负荷旋转套圈主推力滚道的额定寿命L10ev2为:
相对负荷静止套圈副推力滚道的额定寿命L10eu2为:
3)积分公式求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
3.1)当量滚动体负荷
对于双排四点接触球轴承当量滚动体负荷积分如下:
相对负荷方向旋转套圈的当量滚动体负荷积分为:
相对负荷方向静止套圈的当量滚动体负荷积分为:
主推力滚道最大滚动体负荷为:
副推力滚道最大滚动体负荷为:
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev13=Qmax1J1 (4.16)
Qev24=Qmax2J1 (4.17)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu13=Qmax1J2 (4.18)
Qeu24=Qmax2J2 (4.19)
3.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev13为:
相对负荷静止套圈主推力滚道的额定寿命L10eu13为:
相对负荷旋转套圈主推力滚道的额定寿命L10ev24为:
相对负荷静止套圈副推力滚道的额定寿命L10eu24为:
4)静力学方程求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
4.1)当量滚动体负荷
相对旋转套圈主推力滚道的当量滚动体负荷Qev1为:
相对静止套圈主推力滚道的当量滚动体负荷Qeu1为:
相对旋转套圈反推力滚道的当量滚动体负荷Qev2为:
相对静止套圈反推力滚道的当量滚动体负荷Qeu2为:
其中,Q1ψ是作用于主推力滚道各个滚动体上的负荷,Q2ψ是作用于反推力滚道各个滚动体上的负荷;
4.2)接触滚道的寿命,按百万转计算
接触滚道的寿命由方程4.8-4.11求解;
5)静力学方程求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
5.1)当量滚动体负荷
主推力滚道(接触对1、接触对3)的当量滚动体负荷为:
反推力滚道的当量滚动体负荷为:
其中,Q1ψ、Q3ψ是作用于主推力滚道各个滚动体上的负荷,Q2ψ、Q4ψ是作用于反推力滚道各个滚动体上的负荷,Qev1、Qev2、Qev3、Qev4分别为旋转套圈对应的四个接触位置的当量滚动体负荷,Qeu1、Qeu2、Qeu3、Qeu4分别为静止套圈对应的四个接触位置的当量滚动体负荷。
5.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
相对负荷静止套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
相对负荷旋转套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
相对负荷静止套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
所述寿命评估标准的情况如下:
1)积分公式求解单排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.1)
最大受载滚道的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)]-0.9×106 (6.2)
2)积分公式求解双排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)+L10ev24 (-10/9)+L10eu24 (-10/9)]-0.9×106 (6.3)
单个轴承最大受载下的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)]-0.9×106 (6.4)
最大受载滚道的运行转数
L10=L10ev13×106 (6.5)
3)静力学方程求解单排四点接触球轴承寿命,其方法与式子6.1-6.2相同;
4)静力学方程求解双排四点接触球轴承寿命
整个轴承的运行转数为:
单个轴承最大受载下的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.7)
最大受载滚道的运行转数为:
L10=[L10ei1 (-10/9)+L10ee1 (-10/9)]-0.9×106 (6.8)
5)按时间
根据运行时间评估轴承寿命,需要在上述方程的基础上按下式进行求解
其中,n为等效转速;
所述最大动态剪切应力及硬化层深度的情况如下:
式7.1-7.3中,τ0为最大动态剪切应力,Q为接触力,t为滚道与钢球的接触尺寸的结构参数,z0为硬化层深度。
所述工具的参数输入与结果输出的情况如下:
工具界面分为五部分:轴承输入参数、计算方法、载荷导入、计算结果以及结果输出;
在轴承输入参数中需要输入钢球直径、轴承节圆直径、内滚道沟曲率半径系数、外滚道沟曲率半径系数、初始接触角、单排钢球数目、列数、双排滚子间距、轴承游隙、钢球弹性模量、内套圈弹性模量、外套圈弹性模量、内套圈泊松比、外套圈泊松比、滚道硬度、许用接触应力;
在计算方法中选择积分公式或Newton Raphson方法;
在载荷导入中导入包含极限载荷与LDD载荷谱的文件,在载荷输入方法中选择LDD或时序载荷TS,当选择LDD时,无需选择时序路径,当选择TS时需要在载荷导入中选择时序的文件夹;
在计算结果中显示最终计算结果,包括钢球受到最大接触力、钢球受到最大接触应力、最小静安全系数、最大动态剪切应力、硬化层深度、计算寿命;
在结果输出中显示详细的计算结果报告的存放路径。
本发明与现有技术相比,具有如下优点与有益效果:
1、从功能角度:本工具可对时序载荷等效处理,也可直接导入LDD载荷谱,可根据经验公式或静力学方程法求解静强度,可分别通过积分公式法或静力学方程法求解轴承疲劳寿命,同时可求解最大动态剪切应力及硬化层深度。
2、从用户角度:本工具根据目前各大主流轴承供应商内控标准,提出不同寿命评估方法,对于单排四点接触球轴承可分为单排轴承的运行转数,最大受载滚道的运行转数,单排轴承的运行时间,最大受载滚道的运行时间。对于双排四点接触球轴承:双排轴承的运行转数、单排轴承的运行转数,最大受载滚道的运行转数,双排轴承的运行时间,单排轴承的运行时间,最大受载滚道的运行时间。按转数进行评估时应以载荷谱中实际的转数为评价指标。按时间进行评估时应以20年的运行时间为评价指标(对于海上风电25年),全面满足用户的需求。
3、从评估角度:大大缩短了整机产品的设计周期,原先供应商校核结果需要1个月(含报告),采用本系统理想情况下可缩短至1天且可自动输出计算报告。
附图说明
图1为本发明的轴承受载后几何形变示意图。
图2为本发明的轴承套圈上的受力图。
图3为本发明的软件主界面图。
图4为本发明的偏航变桨轴承滚道校核模块界面图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
本实施例所提供的偏航变桨轴承滚道校核工具,为基于Matlab开发的偏航变桨轴承滚道校核软件,主要用于风力发电机偏航、变桨轴承初期选型、认证以及后期的结构强度评估,其主界面如图3所示。所述工具可直接利用LDD载荷谱进行疲劳寿命评估,同时可对时序载荷进等效处理,将每一个工况等效为一组载荷,等效的载荷谱的形式与LDD载荷谱形式一致。偏航变桨轴承一般采用单排四点接触球轴承或双排四点接触球轴承的结构形式,工具可根据经验公式或静力学方程法求解轴承静强度,也可分别通过积分公式法或静力学方程法求解轴承疲劳寿命,同时可求解最大动态剪切应力及硬化层深度。所述工具针对不同轴承供应商的轴承寿命提出以下寿命评估方法:
对于单排四点接触球轴承:可以根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命。
对于双排四点接触球轴承:可以根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命。
其中,按转数进行评估轴承寿命时应以载荷谱中实际的转数为评价指标;按时间进行评估轴承寿命时陆上风电应以20年的运行时间为评价指标,海上风电应以25年的运行时间为评价指标。
以下为本实施例上述偏航变桨轴承滚道校核工具的具体情况:
所述载荷谱的处理的情况如下:
时序载荷:
对每一个时序载荷工况的力、弯矩、速度通过方程1.1-1.3处理得到等效载荷工况表。
式1.1-1.3中,F为轴向力、径向力等效载荷,Fi为时序载荷中第i步的力,M为合弯矩等效载荷,Mi为时序载荷中第i步的弯矩,n为等效转速,ni为时序载荷中第i步的转速。
在20年(海上风电25年)运行时间内,将每一个疲劳工况的发生时间按下述方法求解:
对于发生次数的工况:
T20/25=L×T1×0.05×b (1.4)
对于发生时间的工况:
T20/25=6×L×T1×0.05×b (1.5)
式1.4-1.5中,L为风机的寿命要求,对于海上工况为25年,陆上工况为20年,T1为每个工况发生的次数或小时数,b为每个工况发生的步数,0.05表示步长,T20/25为20年或25年每个疲劳工况发生的时间。
对于LDD载荷可直接导入工具进行寿命计算。
所述静力学方程的情况如下:
静力学方程求解以双排四点接触球轴承为例进行分析。
1)四点接触球轴承的受载后几何形变如下
如图1所示,坐标系采用GL规范中叶根坐标系,假设轴承周边结构具有足够的刚性。外滚道沟渠率中心位置坐标系为xo,yo,内滚道沟渠率中心位置坐标系为xi,yi,轴承外圈固定,内圈相对于外圈产生轴向位移δa,径向位移δr,扭角θ。坐标系位置也由xi,yi转化xi',yi'再到xi”,yi”最后到xi”',yi”',沟渠率中心距离也由最初的MN变为MN”'。
对于双排滚子轴承,由排距d与转角θ对内外滚道产生0.5dθ变形。
对于单排轴承,不存在排距的影响;
2)双排四点接触球轴承的力学模型与负载分布计算
如图2所示,四个接触对在在位置角ψ处的接触力分别表示为Q1ψ、Q2ψ、Q3ψ、Q4ψ,根据Hertz点接触理论,法向接触载荷Qiψ和接触变形δiψ的关系如下:
式中,在位置角ψj处,内圈受到轴向载荷Fa、径向载荷Fr、倾覆力矩M以及钢球对内滚道的接触载荷Qiψ的作用;Kn为滚动体与内外圈总的负荷变形常数。
以内圈与滚动体为研究对象,内圈在外部载荷和所有滚动体载荷的作用下处于平衡状态,内圈的力学平衡方程为:
其中,α1ψ、α2ψ、α3ψ、α4ψ分别为受载后四个接触对在位置角ψ处的接触角。
上述方程通过Newton Raphson方法进行迭代求解,得到滚道精确的载荷分布。
3)单排四点接触球轴承的力学模型与负载分布计算
单排四点接触球轴承的静力学方程求解原理与上述双排四点接触球轴承的计算原理相同。
两个接触对在位置角ψ处的接触力分别表示为Q1ψ、Q2ψ,受载后两个接触对在位置角ψ处的接触角为α1ψ、α2ψ。
故平衡方程为:
所述静强度计算的情况如下:
1)积分公式法的静强度根据下列公式求解
单排轴承:
双排轴承:
根据赫兹接触理论得:
式3.1-3.4中,Qmax为滚动体最大载荷,Fr为径向载荷,Fa为轴向载荷,M为合弯矩,Z为钢球总数,α为初始接触角,dm为节圆直径,a为接触椭圆短半轴,b接触椭圆的长半轴,Smin为最小安全系数,[σ]为许用接触应力,σmax最大接触应力。
2)静力学方程的静强度求解
根据方程2.2-2.7,通过Newton Raphson方法即可求出最大滚动体的载荷Qmax,其余的计算过程与积分公式法一致。
所述疲劳寿命计算情况如下:
1)滚道的基本额定滚动体负荷为:
其中,内滚道取上符号,外滚道取下符号。
式4.1-4.3中,Qci为内滚道的基本额定动负荷,Qce为外滚道的基本额定动负荷,fi为内滚道沟曲率半径系数,fo为外滚道沟曲率半径系数,γ为无量纲几何参数,为Dw为滚子直径,dm为滚道中心径,Z为钢球总个数。
2)积分公式求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命:
2.1)当量滚动体负荷
相对负荷方向旋转套圈的当量滚动体负荷积分为:
相对负荷方向静止套圈的当量滚动体负荷积分为:
相对负荷方向旋转套圈(接触对1和接触对2的内滚道)的当量滚动体负荷为:
Qev1=Qev2=QmaxJ1(ε) (4.6)
相对负荷方向静止套圈(接触对1和接触对2的外滚道)的当量滚动体负荷为:
Qeu1=Qeu2=QmaxJ2(ε) (4.7)
式4.4-4.7中,ε为载荷系数,ψ为方位角,Qmax为滚动体最大载荷,可根据3.1式计算。
2.2)接触滚道的寿命(百万转)
相对负荷旋转套圈主推力滚道(接触对1的内圈)的额定寿命L10ev1为:
相对负荷静止套圈主推力滚道(接触对1的外圈)的额定寿命L10eu1为:
相对负荷旋转套圈主推力滚道(接触对2的内圈)的额定寿命L10ev2为:
相对负荷静止套圈副推力滚道(接触对2的外圈)的额定寿命L10eu2为:
3)积分公式求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
3.1)当量滚动体负荷
对于双排四点接触球轴承当量滚动体负荷积分如下。
相对负荷方向旋转套圈的当量滚动体负荷积分
相对负荷方向静止套圈的当量滚动体负荷积分
主推力滚道最大滚动体负荷为:
副推力滚道最大滚动体负荷为:
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev13=Qmax1J1 (4.16)
Qev24=Qmax2J1 (4.17)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu13=Qmax1J2 (4.18)
Qeu24=Qmax2J2 (4.19)
3.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev13为:
相对负荷静止套圈主推力滚道的额定寿命L10eu13为:
相对负荷旋转套圈主推力滚道的额定寿命L10ev24为:
相对负荷静止套圈副推力滚道的额定寿命L10eu24为:
4)静力学方程求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
4.1)当量滚动体负荷
相对旋转套圈主推力滚道的当量滚动体负荷Qev1为:
相对静止套圈主推力滚道的当量滚动体负荷Qeu1为:
相对旋转套圈反推力滚道的当量滚动体负荷Qev2为:
相对静止套圈反推力滚道的当量滚动体负荷Qeu2为:
其中,Q1ψ是作用于主推力滚道各个滚动体上的负荷,Q2ψ是作用于反推力滚道各个滚动体上的负荷;
4.2)接触滚道的寿命,按百万转计算
接触滚道的寿命由方程4.8-4.11求解;
5)静力学方程求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
5.1)当量滚动体负荷
主推力滚道(接触对1、接触对3)的当量滚动体负荷为:
反推力滚道的当量滚动体负荷为:
其中,Q1ψ、Q3ψ是作用于主推力滚道各个滚动体上的负荷,Q2ψ、Q4ψ是作用于反推力滚道各个滚动体上的负荷,Qev1、Qev2、Qev3、Qev4分别为旋转套圈对应的四个接触位置的当量滚动体负荷,Qeu1、Qeu2、Qeu3、Qeu4分别为静止套圈对应的四个接触位置的当量滚动体负荷。
5.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
相对负荷静止套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
相对负荷旋转套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
相对负荷静止套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
所述寿命评估标准的情况如下:
根据不同的评估标准,寿命计算结果的最终处理也不一样。
1)积分公式求解单排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.1)
最大受载滚道的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)]-0.9×106 (6.2)
2)积分公式求解双排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)+L10ev24 (-10/9)+L10eu24 (-10/9)]-0.9×106 (6.3)
单个轴承(受载最大)的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)]-0.9×106 (6.4)
最大受载滚道的运行转数为:
L10=L10ev13×106 (6.5)
3)静力学方程求解单排四点接触球轴承寿命
方法与式子6.1-6.2一致。
4)静力学方程求解双排四点接触球轴承寿命
整个轴承的运行转数为:
单个轴承(受载最大)的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.7)
最大受载滚道的运行转数为:
L10=[L10ei1 (-10/9)+L10ee1 (-10/9)]-0.9×106 (6.8)
5)按时间
如根据运行时间评估轴承寿命,需要在上述方程的基础上按下式进行求解
其中,n为等效转速。
所述最大动态剪切应力及硬化层深度的情况如下:
式7.1-7.3中,τ0为最大动态剪切应力,Q为接触力,t为滚道与钢球的接触尺寸的结构参数,z0为硬化层深度。
如图4所示,所述工具的参数输入与结果输出的情况如下:
工具界面分为五部分:轴承输入参数、计算方法、载荷导入、计算结果以及结果输出。
在轴承输入参数中需要输入钢球直径、轴承节圆直径、内滚道沟曲率半径系数、外滚道沟曲率半径系数、初始接触角、单排钢球数目、列数、双排滚子间距、轴承游隙、钢球弹性模量、内套圈弹性模量、外套圈弹性模量、内套圈泊松比、外套圈泊松比、滚道硬度、许用接触应力。
在计算方法中选择积分公式或Newton Raphson方法;
在载荷导入中导入包含极限载荷与LDD载荷谱的文件,在载荷输入方法中可以选择LDD或TS(时序载荷),当选择LDD时,无需选择时序路径,当选择TS时需要在载荷导入中选择时序的文件夹。
在计算结果中显示最终计算结果,其中包括钢球受到最大接触力、钢球受到最大接触应力、最小静安全系数、最大动态剪切应力、硬化层深度、计算寿命。
在结果输出中显示详细的计算结果报告的存放路径。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (3)
1.一种偏航变桨轴承滚道校核工具,其特征在于:所述工具为基于Matlab开发的偏航变桨轴承滚道校核软件,主要用于风力发电机偏航、变桨轴承初期选型、认证以及后期的结构强度评估;
所述工具直接利用LDD载荷谱进行疲劳寿命评估,同时能够对时序载荷进等效处理,将每一个工况等效为一组载荷,等效的载荷谱的形式与LDD载荷谱形式一致;
偏航变桨轴承通常采用单排四点接触球轴承或双排四点接触球轴承的结构形式,所述工具能够根据经验公式或静力学方程法求解轴承静强度,也能够通过积分公式法或静力学方程法求解轴承疲劳寿命,同时求解最大动态剪切应力及硬化层深度;
所述工具针对不同轴承供应商的轴承寿命提出以下寿命评估方法:
对于单排四点接触球轴承,根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命;
对于双排四点接触球轴承,根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命;
其中,按转数进行评估轴承寿命时应以载荷谱中实际的转数为评价指标;按时间进行评估轴承寿命时陆上风电应以20年的运行时间为评价指标,海上风电应以25年的运行时间为评价指标。
2.根据权利要求1所述的一种偏航变桨轴承滚道校核工具,其特征在于:
所述时序载荷的处理情况如下:
时序载荷:
对每一个时序载荷工况的力、弯矩、速度通过下面方程1.1-1.3处理得到等效载荷工况表;
式1.1-1.3中,F为轴向力、径向力等效载荷,Fi为时序载荷中第i步的力,M为合弯矩等效载荷,Mi为时序载荷中第i步的弯矩,n为等效转速,ni为时序载荷中第i步的转速;
将每一个疲劳工况的发生时间按下述方法求解:
对于发生次数的工况:
T20/25=L×T1×0.05×b (1.4)
对于发生时间的工况:
T20/25=6×L×T1×0.05×b (1.5)
式1.4-1.5中,L为风机的寿命要求,对于海上工况为25年,陆上工况20年,T1为每个工况发生的次数或小时数,b为每个工况发生的步数,0.05表示步长,T20/25为20年或25年每个疲劳工况发生的时间;
对于LDD载荷则可直接导入工具进行寿命计算;
所述静力学方程的情况如下:
1)四点接触球轴承的受载后几何形变如下:
坐标系采用GL规范中叶根坐标系,假设轴承周边结构具有足够的刚性,外滚道沟渠率中心位置坐标系为xo,yo,内滚道沟渠率中心位置坐标系为xi,yi,轴承外圈固定,内圈相对于外圈产生轴向位移δa,径向位移δr,扭角θ;坐标系位置也由xi,yi转化xi',yi'再到xi”,yi”最后到xi”',yi”',沟渠率中心距离也由最初的MN变为MN”';
对于双排滚子轴承,由排距d与转角θ对内外滚道产生0.5dθ变形;
对于单排轴承,不存在排距的影响;
2)双排四点接触球轴承的力学模型与负载分布计算
以内圈与滚动体为研究对象,内圈在外部载荷和所有滚动体载荷的作用下处于平衡状态,内圈的力学平衡方程为:
上述方程通过Newton Raphson方法进行迭代求解,得到滚道精确的载荷分布;
3)单排四点接触球轴承的力学模型与负载分布计算
单排四点接触球轴承的静力学方程求解原理与上述双排四点接触球轴承的计算原理相同;
故平衡方程为:
所述静强度计算的情况具体如下:
1)积分公式法的静强度根据下列公式求解
单排轴承:
双排轴承:
根据赫兹接触理论得:
式3.1-3.4中,Qmax为滚动体最大载荷,Fr为径向载荷,Fa为轴向载荷,M为合弯矩,Z为钢球总数,α为初始接触角,dm为节圆直径,a为接触椭圆短半轴,b为接触椭圆的长半轴,Smin为最小安全系数,[σ]为许用接触应力,σmax为最大接触应力;
2)静力学方程的静强度求解
根据方程2.2-2.7,通过Newton Raphson方法即可求出最大滚动体的载荷Qmax,其余的计算过程与积分公式法一致;
所述疲劳寿命的计算情况如下:
1)滚道的基本额定滚动体负荷:
其中,内滚道取上符号,外滚道取下符号;
式4.1-4.3中,Qci为内滚道的基本额定动负荷,Qce为外滚道的基本额定动负荷,fi为内滚道沟曲率半径系数,fo为外滚道沟曲率半径系数,γ为无量纲几何参数,Dw为滚子直径,dm为节圆直径,Z为钢球总个数;
2)积分公式求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
2.1)当量滚动体负荷
相对负荷方向旋转套圈的当量滚动体负荷积分为:
相对负荷方向静止套圈的当量滚动体负荷积分为:
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev1=Qev2=QmaxJ1(ε) (4.6)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu1=Qeu2=QmaxJ2(ε) (4.7)
式4.4-4.7中,ε为载荷系数,ψ为位置角,Qmax为滚动体最大载荷,根据式3.1计算;
2.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev1为:
相对负荷静止套圈主推力滚道的额定寿命L10eu1为:
相对负荷旋转套圈主推力滚道的额定寿命L10ev2为:
相对负荷静止套圈副推力滚道的额定寿命L10eu2为:
3)积分公式求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
3.1)当量滚动体负荷
对于双排四点接触球轴承当量滚动体负荷积分如下:
相对负荷方向旋转套圈的当量滚动体负荷积分为:
相对负荷方向静止套圈的当量滚动体负荷积分为:
主推力滚道最大滚动体负荷为:
副推力滚道最大滚动体负荷为:
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev13=Qmax1J1 (4.16)
Qev24=Qmax2J1 (4.17)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu13=Qmax1J2 (4.18)
Qeu24=Qmax2J2 (4.19)
3.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev13为:
相对负荷静止套圈主推力滚道的额定寿命L10eu13为:
相对负荷旋转套圈主推力滚道的额定寿命L10ev24为:
相对负荷静止套圈副推力滚道的额定寿命L10eu24为:
4)静力学方程求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
4.1)当量滚动体负荷
相对旋转套圈主推力滚道的当量滚动体负荷Qev1为:
相对静止套圈主推力滚道的当量滚动体负荷Qeu1为:
相对旋转套圈反推力滚道的当量滚动体负荷Qev2为:
相对静止套圈反推力滚道的当量滚动体负荷Qeu2为:
4.2)接触滚道的寿命,按百万转计算
接触滚道的寿命由方程4.8-4.11求解;
5)静力学方程求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
5.1)当量滚动体负荷
主推力滚道的当量滚动体负荷为:
反推力滚道的当量滚动体负荷为:
其中,是作用于主推力滚道各个滚动体上的负荷,是作用于反推力滚道各个滚动体上的负荷,Qev1、Qev2、Qev3、Qev4分别为旋转套圈对应的四个接触位置的当量滚动体负荷,Qeu1、Qeu2、Qeu3、Qeu4分别为静止套圈对应的四个接触位置的当量滚动体负荷;
5.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命分别为:
相对负荷静止套圈主推力滚道的额定寿命分别为:
相对负荷旋转套圈副推力滚道的额定寿命分别为:
相对负荷静止套圈副推力滚道的额定寿命分别为:
所述寿命评估标准的情况如下:
1)积分公式求解单排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.1)
最大受载滚道的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)]-0.9×106 (6.2)
2)积分公式求解双排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)+L10ev24 (-10/9)+L10eu24 (-10/9)]-0.9×106 (6.3)
单个轴承最大受载下的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)]-0.9×106 (6.4)
最大受载滚道的运行转数为:
L10=L10ev13×106 (6.5)
3)静力学方程求解单排四点接触球轴承寿命,其方法与方程6.1-6.2相同;
4)静力学方程求解双排四点接触球轴承寿命
整个轴承的运行转数为:
单个轴承最大受载下的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.7)
最大受载滚道的运行转数为:
L10=[L10ei1 (-10/9)+L10ee1 (-10/9)]-0.9×106 (6.8)
5)按时间
根据运行时间评估轴承寿命,需要在上述方程的基础上按下式进行求解
其中,n为等效转速;
所述最大动态剪切应力及硬化层深度的情况如下:
式7.1-7.3中,τ0为最大动态剪切应力,Q为接触力,t为滚道与钢球的接触尺寸的结构参数,z0为硬化层深度。
3.根据权利要求1所述的一种偏航变桨轴承滚道校核工具,其特征在于:所述工具的参数输入与结果输出的情况如下:
工具界面分为五部分:轴承输入参数、计算方法、载荷导入、计算结果以及结果输出;
在轴承输入参数中需要输入钢球直径、轴承节圆直径、内滚道沟曲率半径系数、外滚道沟曲率半径系数、初始接触角、单排钢球数目、列数、双排滚子间距、轴承游隙、钢球弹性模量、内套圈弹性模量、外套圈弹性模量、内套圈泊松比、外套圈泊松比、滚道硬度、许用接触应力;
在计算方法中选择积分公式或Newton Raphson方法;
在载荷导入中导入包含极限载荷与LDD载荷谱的文件,在载荷输入方法中选择LDD或时序载荷TS,当选择LDD时,无需选择时序路径,当选择TS时需要在载荷导入中选择时序的文件夹;
在计算结果中显示最终计算结果,包括钢球受到最大接触力、钢球受到最大接触应力、最小静安全系数、最大动态剪切应力、硬化层深度、计算寿命;
在结果输出中显示详细的计算结果报告的存放路径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710753617.4A CN107451376B (zh) | 2017-08-29 | 2017-08-29 | 一种偏航变桨轴承滚道校核工具 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710753617.4A CN107451376B (zh) | 2017-08-29 | 2017-08-29 | 一种偏航变桨轴承滚道校核工具 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107451376A CN107451376A (zh) | 2017-12-08 |
CN107451376B true CN107451376B (zh) | 2020-07-14 |
Family
ID=60494219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710753617.4A Active CN107451376B (zh) | 2017-08-29 | 2017-08-29 | 一种偏航变桨轴承滚道校核工具 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107451376B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103256908A (zh) * | 2013-05-22 | 2013-08-21 | 河南科技大学 | 一种变桨轴承径向游隙的确定方法 |
CN103324805A (zh) * | 2013-07-04 | 2013-09-25 | 沈阳华创风能有限公司 | 风机底架极限强度校核的计算方法 |
CN105550403A (zh) * | 2015-12-04 | 2016-05-04 | 广东明阳风电产业集团有限公司 | 一种风机外购件强度校核工具 |
-
2017
- 2017-08-29 CN CN201710753617.4A patent/CN107451376B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103256908A (zh) * | 2013-05-22 | 2013-08-21 | 河南科技大学 | 一种变桨轴承径向游隙的确定方法 |
CN103324805A (zh) * | 2013-07-04 | 2013-09-25 | 沈阳华创风能有限公司 | 风机底架极限强度校核的计算方法 |
CN105550403A (zh) * | 2015-12-04 | 2016-05-04 | 广东明阳风电产业集团有限公司 | 一种风机外购件强度校核工具 |
Also Published As
Publication number | Publication date |
---|---|
CN107451376A (zh) | 2017-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107729597B (zh) | 一种主轴轴承滚道校核工具 | |
CN107704652B (zh) | 一种风力发电机组轴承刚度计算工具 | |
CN111896255B (zh) | 基于多轴随机道路载荷下轮毂轴承服役寿命快速评估方法 | |
EP2532904A2 (en) | Structure of slewing ring bearing | |
CN102819635B (zh) | 特大型双排四点接触球轴承静力学模型建立方法 | |
CN105550403B (zh) | 一种风机外购件强度校核工具 | |
WO2007061439A3 (en) | Wind turbine pitch bearing and method | |
CN109753723B (zh) | 一种向心滚动轴承疲劳寿命计算方法 | |
Stammler et al. | Friction torque of wind-turbine pitch bearings–comparison of experimental results with available models | |
CN108843760A (zh) | 变速装置 | |
CN103256908B (zh) | 一种变桨轴承径向游隙的确定方法 | |
CN108692938B (zh) | 一种获取滚动轴承寿命的方法 | |
CN103674546B (zh) | 风电转盘轴承加速疲劳寿命试验方法 | |
CN112861271B (zh) | 风电机组主轴承力学特性的计算方法、系统及介质 | |
CN107451376B (zh) | 一种偏航变桨轴承滚道校核工具 | |
CN103174740A (zh) | 一种四点接触球轴承沟曲率半径系数的设计方法 | |
CN107563081B (zh) | 高速空心圆柱滚子轴承的参数优化方法 | |
RU2626800C1 (ru) | Способ комплектования шарикоподшипников | |
CN106934165A (zh) | 一种交叉圆锥滚子轴承静刚度的快速提取方法 | |
Rabréau et al. | Influence of bearing kinematics hypotheses on ball bearing heat generation | |
CN111859575A (zh) | 一种单列四点接触球转盘轴承额定动静载荷计算方法 | |
Sun et al. | Rolling-element bearings in China: From ancient times to the 20th century | |
CN106640559A (zh) | 非对称双排四点接触球轴承非对称角的设置方法以及装置 | |
CN211288484U (zh) | 一种s形定位回转径结构的偏航轴承 | |
CN211174993U (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201216 Address after: No.2, Luohu West Road, Houpu village committee, Jieshi Town, Lufeng City, Shanwei City, Guangdong Province Patentee after: Shanwei Mingyang New Energy Technology Co.,Ltd. Address before: 528437 No. 22 Torch Road, Torch Development Zone, Zhongshan City, Guangdong Province Patentee before: MING YANG SMART ENERGY GROUP Co.,Ltd. |
|
TR01 | Transfer of patent right |