CN113627046A - 一种风机偏航制动盘建模与极限强度校核方法 - Google Patents
一种风机偏航制动盘建模与极限强度校核方法 Download PDFInfo
- Publication number
- CN113627046A CN113627046A CN202110783542.0A CN202110783542A CN113627046A CN 113627046 A CN113627046 A CN 113627046A CN 202110783542 A CN202110783542 A CN 202110783542A CN 113627046 A CN113627046 A CN 113627046A
- Authority
- CN
- China
- Prior art keywords
- yaw
- yaw brake
- brake disc
- load
- modeling
- 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
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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Braking Arrangements (AREA)
Abstract
本发明涉及一种风机偏航制动盘建模与极限强度校核方法,建立包含风机偏航制动载荷传递路径的几何模型,进行网格划分,进行接触对定义,施加边界条件与各极限工况载荷,通过在塔顶法兰中心位置建立两个主节点,将一个主节点与主机架前端法兰通过刚性梁连接,加载3个力及2个力矩载荷;另一主节点与偏航制动卡钳通过刚性梁连接,加载不超过偏航制动卡钳最大制动力矩的Mz载荷,超过部分加载到与主机架相连的主节点,将载荷传递至主机架与偏航制动卡钳,最后传递至偏航制动盘。本发明既考虑了偏航轴承的非线性传递力关系又考虑了偏航制动盘与偏航制动卡钳之间接触力的传递,符合偏航制动盘实际受力关系,能够实现风机偏航制动盘极限强度的准确校核。
Description
技术领域
本发明涉及风机航制动盘技术领域,尤其涉及一种风机偏航制动盘建模与极限强度校核方法。
背景技术
随着社会经济的快速发展,人们对能源的需求也在迅猛增长。而风能作为一种清洁无污染的可再生能源,越来越受到人们的重视。由于国家政策的大力支持,近几年风电机组的国产化程度逐渐提高。风电机组功率的逐渐增大使得机组内部结构件也越来越大,结构件的强度性能也愈发重要.
偏航制动盘是风机组中受力情况较为复杂、可靠性要求较高的关键部件之一,承受着来自风轮总成与机舱总成传递过来的周期性载荷和随机载荷的综合作用,其强度性能直接关系到整个风电机组的安全性与可靠性。
现在大型风机组偏航制动盘主要采用有限元方法进行计算.在现有的偏航制动盘强度性能分析技术中,如申请公布号为CN106547963 A,申请公布日为2017年03月29日的发明专利申请公开了一种风机偏航制动盘极限强度的校核方法及装置,该方法采用有限元软件建立包含制动盘在内的实体模型,并对其进行强度分析。首先建立包含风机偏航载荷传递路径的几何模型,将其导入有限元分析软件进行网格划分;其次在偏航轴承滚球建模方面采用只承受压力不承受拉力的杆单元单元模拟,而在偏航制动盘与偏航制动卡钳的接触关系处理上,则采用了绑定接触的方法模拟;再通过在塔顶法兰中心位置建立主节点,将主节点与主机架前端法兰通过刚性梁单元连接起来;最后在塔顶法兰中心位置施加偏航载荷,通过刚性梁单元将载荷传递至主机架,经过偏航制动卡钳传递至偏航制动盘,以实现对偏航制动盘受力关系的模拟。
现有的偏航制动盘强度计算方法在偏航制动盘与偏航制动卡钳的接触关系上,主要采用绑定接触的方法模拟,而实际上两者之间是摩擦接触的关系,当采用绑定接触时会强制两者同步运动,无法体现两者之间可以发生分离或滑动的真实接触关系,这样会导致偏航制动盘的最终计算结果偏大。同时模拟偏航轴承滚球时普遍采用了只承受压力的杆单元link180建模的方法,虽然能体现出偏航轴承传力时只承受压力不承受拉力的非线性关系,但该单元仅能赋予确定刚度的实常数来描述滚球刚度,而无法体现出滚球实际刚度的变化非线性,即现有技术在偏航制动盘与偏航制动卡钳连接关系及偏航轴承非线性模拟方面考虑的不够全面,不符合实际工作情况,导致偏航制动盘的强度校核不够准确。
发明内容
为解决目前杆单元模拟偏航轴承滚球和绑定接触模拟偏航制动盘及偏航制动卡钳之间关系的方法得到的风机偏航制动盘强度不够准确的问题,本发明提供了一种风机偏航制动盘建模与极限强度校核方法,能够更加真实而准确的模拟制动盘的受力特性。
为达到上述目的,本发明提供了一种风机偏航制动盘建模与极限强度校核方法,包括:
1)对偏航制动盘及其载荷传递路径中的相关部件进行三维几何建模与模型装配,得到三维几何模型;
2)利用有限元前处理软件对所述三维几何模型进行网格划分,对载荷传递路径中的相关部件以及偏航制动盘进行属性定义,得到网格模型;
3)将所述网格模型导入有限元分析软件,并依据相邻部件之间的实际接触关系定义相关接触对,在偏航制动盘与偏航制动卡钳之间建立弹簧单元来模拟两者之间接触关系,完成有限元建模;
4)对有限元模型进行边界约束与极限工况加载;
5)对各极限工况下有限元模型进行非线性求解,求解完成后进入有限元分析软件后处理模块,查看各极限工况下偏航制动盘的应力大小及危险位置,将所述危险位置的部件许用应力值与各极限工况下偏航制动盘的应力大小比较,校核偏航制动盘极限强度,得到极限强度安全系数。
进一步地,所述载荷传递路径相关部件包括:主机架、偏航驱动电机、偏航轴承内圈、偏航轴承外圈、偏航轴承滚球、偏航制动卡钳、偏航制动盘、塔筒顶法兰及塔筒假体部分。
进一步地,所述步骤2)网格划分,包括:
所述主机架、偏航驱动电机和偏航制动盘均采用十节点四面体高阶单元进行网格划分;
所述偏航制动卡钳、偏航轴承内圈、偏航轴承外圈和塔顶法兰均采用二十节点六面体高阶单元进行网格划分;
所述偏航轴承与偏航驱动电机之间的偏航齿轮啮合采用梁单元与杆单元进行建模;
所述塔筒假体部分采用壳单元进行网格划分;
所述偏航轴承滚球采用轴向拉压弹簧单元建模;
所述的杆单元其中通过单元实常数定义设置为仅受压不受拉的特性。
进一步地,所述梁单元和杆单元划分份数均为8~16。
进一步地,所述偏航轴承内圈和所述偏航轴承外圈之间的偏航轴承滚球按照实际个数采用弹簧单元combin39进行建模,每个所述偏航轴承滚球采用六个所述combin39单元辐条轮模拟,六个所述combin39单元相对于滚道水平面分别成±30°、±50°、±70°分布,并采用预定义的力-变形曲线来定义所述偏航轴承滚球承受压力载荷时的刚度变化。
进一步地,所述偏航制动盘分别与8个偏航制动卡钳的上、下两个面发生接触,两者之间的接触面摩擦采用轴向拉压弹簧单元进行建模;每个接触面均采用X、Y、Z 3个方向拉压弹簧来分别模拟偏航制动盘与偏航制动卡钳的X、Y方向横向摩擦力与Z方向的垂向压力;
所述偏航制动盘与偏航制动卡钳之间的作用力共采用24个弹簧单元模拟。
进一步地,每个接触面上的采用在同一位置建立两个节点,两个节点分别位于接触面和目标面上,并在重合的两个节点间建立3个方向的弹簧单元,而弹簧的两端分别通过刚性梁与偏航制动盘及偏航制动卡钳连接;
Z方向的弹簧单元模拟偏航制动卡钳与偏航制动盘之间的垂向压力,具有仅受压不受拉特性,采用预定义的力-变形曲线来定义弹簧单元承受压力载荷时的刚度变化,可根据单个偏航制动卡钳的夹紧力与摩擦片压缩量得到Z方向的力-变形曲线;
X、Y方向的弹簧单元模拟偏航制动卡钳与偏航制动盘之间的横向摩擦力,具有可受压可受拉特性,根据单个偏航制动卡钳的夹紧力及其与偏航制动盘之间的摩擦系数得到横向摩擦力,再采用横向摩擦力与滑移量得到X、Y方向的力-变形曲线。
进一步地,所述步骤3)在进行有限元建模时,所述主机架与偏航制动卡钳及偏航驱动电机之间、偏航轴承与偏航制动盘之间、偏航制动盘与塔顶法兰之间、塔顶法兰与塔筒假体之间的接触关系均设置为绑定接触.
进一步地,所述步骤4)在进行边界约束与加载计算时,包括:
对所述塔筒假体底部进行所有自由度约束;
在所述塔顶法兰中心位置建立两个重合的主节点,并将两个主节点分别与主机架前、后端的塔筒顶法兰及偏航制动卡钳后端面通过刚性梁连接起来,将塔顶中心三个力分量Fx、Fy、Fz与3个力矩分量Mx,My,Mz加载到这两个主节点上。
进一步地,将塔顶中心载荷3个力分量Fx、Fy、Fz与2个力矩Mx、My加载到与主机架相连的主节点上,将小于偏航制动卡钳所能提供的最大制动力矩的载荷分量Mz加载到与偏航制动卡钳相连的主节点上,将超出部分也加载到与主机架相连的主节点上,从而使得塔顶中心载荷通过刚性梁传递至主机架和偏航制动卡钳,最终传递至偏航制动盘,以此来模拟偏航制动盘的受力情况。
进一步地,所述步骤5)在强度结果评估时,查看各极限工况下偏航制动盘应力大小及危险位置,排除偏航制动盘与偏航轴承外圈及塔顶法兰之间的绑定区域,因绑定接触导致该区域结果不准确,将应力最大值与部件许用应力作比较,评估偏航制动盘的最小极限安全系数时需考虑材料安全系数1.1。
本发明的上述技术方案具有如下有益的技术效果:
(1)本发明在有限元建模与强度校核方面既考虑了偏航轴承的非线性传递力关系又考虑了偏航制动盘与偏航制动卡钳之间接触力的传递,符合偏航制动盘实际受力关系,能够实现风机偏航制动盘极限强度的准确校核。
(2)本发明在有限元建模时采用可定义拉压与刚度变化的弹簧单元来模拟偏航轴承滚球,并采用3个方向拉压弹簧来模拟偏航制动盘及偏航制动卡钳之间接触关系,更加全面地体现了偏航轴承本身的非线性传力关系及偏航制动卡钳与偏航制动盘之间的接触传力特性,避免了因偏航制动盘与偏航制动卡钳之间绑定接触造成的偏航制动盘应力过大的问题,从而使得偏航制动盘受力情形模拟得更加真实,同时也降低了由于偏航制动卡钳与偏航制动盘之间采用摩擦接触导致的不收敛难题。
(3)本发明得到的偏航制动盘的应力大小和危险位置准确度更高,符合偏航制动盘实际受力关系,能够实现风机偏航制动盘极限强度的准确校核。
附图说明
图1是本发明风机偏航制动盘建模与极限强度校核方法的流程图;
图2是实施的有限元模型剖面示意图;
图3是偏航驱动装置连接示意图;
图4是偏航轴承滚球部分计算模型示意图;
图5是偏航轴承滚球力-变形曲线示意图;
图6是偏航制动盘与偏航制动卡钳之间弹簧单元示意图;
图7是偏航制动盘与偏航卡钳之间弹簧单元X、Y方向弹簧力-变形曲线示意图;
图8是偏航制动盘与偏航卡钳之间弹簧单元Z方向弹簧力-变形曲线示意图;
图9为偏航制动盘几何装配模型。
图中:1-主机架、2-塔顶法兰、3-塔筒假体、4-与主机架前、后端法兰连接的刚性梁连接、5-与偏航制动卡钳后端面连接的刚性梁连接、6-偏航驱动电机、7-偏航驱动齿轮、8-偏航轴承内圈、9-偏航轴承外圈、10-偏航轴承滚球、11-偏航制动卡钳、12-偏航制动盘、13-偏航制动盘与偏航制动卡钳之间的弹簧接触连接。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面结合具体实施方式并参照附图,对本发明进一步详细说明。应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。
图9为简化的偏航制动盘几何装配模型,除偏航制动卡钳11与偏航制动盘12之间采用弹簧模拟摩擦接触外,其余各相邻部件之间采用绑定接触连接。首先载荷经过主机架1,通过与偏航制动卡钳11之间的绑定接触将载荷传递至偏航制动卡钳11,再通过偏航制动卡钳11与偏航制动盘12之间的摩擦接触将载荷传递至偏航制动盘12,经过偏航制动盘12与塔顶法兰2之间绑定接触将载荷传递至塔顶法兰,同理,通过塔顶法兰2与塔筒3之间的绑定接触将载荷传递至塔筒3,再通过塔筒3将载荷传递到大地中;而超出偏航制动卡钳最大制动力矩的Mz载荷则通过与偏航驱动电机6之间的绑定接触传递至偏航驱动电机6,偏航驱动电机6再投入偏航刹车来抵抗此载荷。
本发明提供了一种风机偏航制动盘建模与极限强度校核方法,结合图1,包括如下步骤:
1)利用三维绘图软件,对所述偏航制动盘及其载荷传递路径相关部件进行三维几何建模与模型装配,得到三维几何模型,其中所述载荷传递路径相关部件包括:主机架1、偏航驱动电机6、偏航轴承内圈8、偏航轴承外圈9、偏航轴承滚球10、偏航制动卡钳11、偏航制动盘12、塔筒顶法兰2及塔筒假体部分3。
2)利用有限元前处理软件对所述三维几何模型进行网格划分,并对各组件主机架1、偏航驱动电机6、偏航轴承内圈8、偏航轴承外圈9、偏航轴承滚球10、偏航制动卡钳11、偏航制动盘12、塔筒顶法兰2、塔筒假体部分3进行单元、材料、实常数、横截面等属性定义,得到网格模型。对偏航制动盘实体模型转化为网格模型定义单元类型。有限元模型剖面示意图,如图2所示。
3)将所述网格模型导入有限元分析软件,并依据相邻部件之间的实际接触关系定义相关接触对,在偏航制动盘与偏航制动卡钳之间建立弹簧单元来模拟两者之间接触关系,完成有限元建模。
4)对所述有限元模型进行边界约束与极限工况加载。
5)对各极限工况下有限元模型进行非线性求解,求解完成后进入有限元分析软件后处理模块,查看各极限工况下偏航制动盘的应力大小及危险位置,将该值与部件许用应力值比较即可校核偏航制动盘极限强度,得到极限强度安全系数。
进一步地,所述步骤2)在进行网格划分时,包括:
所述主机架、偏航驱动电机、偏航制动盘均采用十节点四面体高阶单元进行网格划分;
所述偏航制动卡钳、偏航轴承内圈/外圈、塔顶法兰均采用二十节点六面体高阶单元进行网格划分;
所述的塔筒假体部分采用壳单元进行网格划分。
偏航轴承滚球部分计算模型示意图,如图4所示。偏航轴承滚球力-变形曲线示意图如图5所示,所述偏航轴承与偏航驱动电机之间的偏航齿轮啮合采用梁单元与杆单元进行建模,以实现主机架与偏航轴承外圈之间的载荷传递。
进一步地,所述梁单元和杆单元划分份数均为8~16。
进一步地,通过单元实常数定义所述杆单元仅受压不受拉特性。
进一步地,所述偏航轴承内圈和所述偏航轴承外圈之间的偏航轴承滚球按照实际个数采用弹簧单元combin39进行建模,每一个偏航轴承滚球采用一圈6根弹簧单元combin39轮辐来模拟,且六个弹簧单元combin39相对于滚道水平面分别成±30°、±50°、±70°分布。
进一步地,通过定义实常数即预定义的力-变形曲线赋予弹簧单元,来实现对偏航轴承滚球只受压不受拉及其承受压力载荷时的刚度变化等非线性属性的模拟。
进一步地,所述步骤3)在进行接触对定义时,包括:
所述主机架1与偏航制动卡钳11及偏航驱动电机6之间、偏航轴承与偏航制动盘12之间、偏航制动盘12与塔顶法兰2之间、塔顶法兰2与塔筒假体3之间的接触关系均设置为绑定接触。
所述的偏航轴承滚球10采用轴向拉压弹簧单元建模。
所述偏航制动盘12与偏航制动卡钳11之间的接触采用3个分布于X、Y、Z方向的拉压弹簧单元combin39进行建模,以此来模拟偏航制动盘与偏航制动卡钳之间的作用力,即X、Y方向的横向摩擦力与Z方向的垂直压力。图6为偏航制动盘与偏航制动卡钳之间弹簧单元示意图。图7是偏航制动盘与偏航卡钳之间弹簧单元X、Y方向弹簧力-变形曲线示意图;图8是偏航制动盘与偏航卡钳之间弹簧单元Z方向弹簧力-变形曲线示意图。
进一步地,通过实常数来定义,通过定义三个方向的实常数即力-变形曲线赋予对应方向的弹簧单元,来模拟偏航制动盘与偏航制动卡钳之间的横向摩擦力与垂直压力作用。
进一步地,所述步骤4)在进行边界约束与加载计算时,包括:
对所述塔筒假体3底部进行所有自由度约束。
在所述塔顶法兰2中心位置建立两个主节点,两个节点分别位于接触面和目标面上,并将两个主节点分别与主机架前端法兰及偏航制动卡钳后端面通过刚性梁4连接起来,将各极限工况载荷即塔顶中心三个力分量Fx、Fy、Fz与3个力矩分量Mx,My,Mz加载到这两个主节点上。
进一步地,将塔顶中心载荷3个力分量Fx、Fy、Fz与2个力矩Mx、My加载到与与主机架1相连的主节点上,将小于偏航制动卡钳11所能提供的最大制动力矩的载荷分量Mz加载到与偏航制动卡钳11相连的主节点上,将超出部分也加载到与主机架1相连的主节点上,从而塔顶中心载荷通过刚性梁传递至主机架1和偏航制动卡钳11,最终传递至偏航制动盘,以模拟偏航制动盘的受力情况。
进一步地,所述步骤5)在评估偏航制动盘计算结果时,因绑定接触导致绑定区域结果不准确予以忽略,同时根据GL2010规范,评估偏航制动盘的最小极限强度安全系数时需考虑材料安全系数1.1。
现有技术中关于偏航制动盘建模与极限强度校核方法,大多采用杆单元模拟偏航轴承滚球,并设置绑定接触模拟偏航制动盘及偏航制动卡钳之间的接触关系,载荷全部加载到塔顶中心的主节点,并通过刚性梁与主机架前端法兰连接将载荷由主机架传递至偏航制动盘,以此来评估偏航制动盘强度性能。
本发明首先建立包含风机偏航制动载荷传递路径的几何模型,然后将该模型导入至有限元前处理软件进行网格划分,得到网格模型,再将其导入有限元分析软件进行接触对定义,施加边界条件与各极限工况载荷,通过在塔顶法兰中心位置建立两个主节点,将一个主节点与主机架前端法兰通过刚性梁连接起来,在该位置加载3个力及2个力矩载荷,另一个主节点与偏航制动卡钳通过刚性梁连接,在该主节点加载不超过偏航制动卡钳最大制动力矩的Mz载荷,超过部分加载到与主机架相连的主节点上,以此将载荷传递至主机架与偏航制动卡钳,最后传递至偏航制动盘,以实现对偏航制动盘受力关系的准确模拟。本发明在有限元建模时采用可定义拉压与刚度变化的弹簧单元来模拟偏航轴承滚球,并采用3个方向拉压弹簧来模拟偏航制动盘及偏航制动卡钳之间接触关系,更加全面地体现了偏航轴承本身的非线性传力关系及偏航制动卡钳与偏航制动盘之间的接触传力特性,避免了因偏航制动盘与偏航制动卡钳之间绑定接触造成的偏航制动盘应力过大的问题,从而使得偏航制动盘受力情形模拟得更加真实,同时也降低了由于偏航制动卡钳与偏航制动盘之间采用摩擦接触导致的不收敛难题。通过上述方法,本发明得到的偏航制动盘的应力大小和危险位置准确度更高,符合偏航制动盘实际受力关系,能够实现风机偏航制动盘极限强度的准确校核。
应当理解的是,本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理,而不构成对本发明的限制。因此,在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。此外,本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。
Claims (10)
1.一种风机偏航制动盘建模与极限强度校核方法,其特征在于,该方法包括:
1)对偏航制动盘及其载荷传递路径中的相关部件进行三维几何建模与模型装配,得到三维几何模型;
2)利用有限元前处理软件对所述三维几何模型进行网格划分,对载荷传递路径中的相关部件以及偏航制动盘进行属性定义,得到网格模型;
3)将所述网格模型导入有限元分析软件,并依据相邻部件之间的实际接触关系定义相关接触对,在偏航制动盘与偏航制动卡钳之间建立弹簧单元来模拟两者之间接触关系,完成有限元建模;
4)对有限元模型进行边界约束与极限工况加载;
5)对各极限工况下有限元模型进行非线性求解,求解完成后进入有限元分析软件后处理模块,查看各极限工况下偏航制动盘的应力大小及危险位置,将所述危险位置的部件许用应力值与各极限工况下偏航制动盘的应力大小比较,校核偏航制动盘极限强度,得到极限强度安全系数。
2.根据权利要求1所述的风机偏航制动盘建模与极限强度校核方法,其特征在于,所述载荷传递路径相关部件包括:主机架、偏航驱动电机、偏航轴承内圈、偏航轴承外圈、偏航轴承滚球、偏航制动卡钳、偏航制动盘、塔筒顶法兰及塔筒假体部分。
3.根据权利要求2所述的风机偏航制动盘建模与极限强度校核方法,其特征在于,所述步骤2)网格划分,包括:
所述主机架、偏航驱动电机和偏航制动盘均采用十节点四面体高阶单元进行网格划分;
所述偏航制动卡钳、偏航轴承内圈、偏航轴承外圈和塔顶法兰均采用二十节点六面体高阶单元进行网格划分;
所述偏航轴承与偏航驱动电机之间的偏航齿轮啮合采用梁单元与杆单元进行建模;
所述塔筒假体部分采用壳单元进行网格划分;
所述偏航轴承滚球采用轴向拉压弹簧单元建模;
所述的杆单元其中通过单元实常数定义设置为仅受压不受拉的特性。
进一步地,所述梁单元和杆单元划分份数均为8~16。
4.根据权利要求3所述风机偏航制动盘建模与极限强度校核方法,其特征在于,所述偏航轴承内圈和所述偏航轴承外圈之间的偏航轴承滚球按照实际个数采用弹簧单元combin39进行建模,每个所述偏航轴承滚球采用六个所述combin39单元辐条轮模拟,六个所述combin39单元相对于滚道水平面分别成±30°、±50°、±70°分布,并采用预定义的力-变形曲线来定义所述偏航轴承滚球承受压力载荷时的刚度变化。
5.根据权利要求3所述风机偏航制动盘建模与极限强度校核方法,其特征在于,所述偏航制动盘分别与8个偏航制动卡钳的上、下两个面发生接触,两者之间的接触面摩擦采用轴向拉压弹簧单元进行建模;每个接触面均采用X、Y、Z 3个方向拉压弹簧来分别模拟偏航制动盘与偏航制动卡钳的X、Y方向横向摩擦力与Z方向的垂向压力;
所述偏航制动盘与偏航制动卡钳之间的作用力共采用24个弹簧单元模拟。
6.根据权利要求5所述风机偏航制动盘建模与极限强度校核方法,其特征在于,每个接触面上的采用在同一位置建立两个节点,两个节点分别位于接触面和目标面上,并在重合的两个节点间建立3个方向的弹簧单元,而弹簧的两端分别通过刚性梁与偏航制动盘及偏航制动卡钳连接;
Z方向的弹簧单元模拟偏航制动卡钳与偏航制动盘之间的垂向压力,具有仅受压不受拉特性,采用预定义的力-变形曲线来定义弹簧单元承受压力载荷时的刚度变化,可根据单个偏航制动卡钳的夹紧力与摩擦片压缩量得到Z方向的力-变形曲线;
X、Y方向的弹簧单元模拟偏航制动卡钳与偏航制动盘之间的横向摩擦力,具有可受压可受拉特性,根据单个偏航制动卡钳的夹紧力及其与偏航制动盘之间的摩擦系数得到横向摩擦力,再采用横向摩擦力与滑移量得到X、Y方向的力-变形曲线。
7.根据权利要求2所述风机偏航制动盘建模与极限强度校核方法,其特征在于,所述步骤3)在进行有限元建模时,所述主机架与偏航制动卡钳及偏航驱动电机之间、偏航轴承与偏航制动盘之间、偏航制动盘与塔顶法兰之间、塔顶法兰与塔筒假体之间的接触关系均设置为绑定接触。
8.根据权利要求2所述风机偏航制动盘建模与极限强度校核方法,其特征在于,所述步骤4)在进行边界约束与加载计算时,包括:
对所述塔筒假体底部进行所有自由度约束;
在所述塔顶法兰中心位置建立两个重合的主节点,并将两个主节点分别与主机架前、后端的塔筒顶法兰及偏航制动卡钳后端面通过刚性梁连接起来,将塔顶中心三个力分量Fx、Fy、Fz与3个力矩分量Mx,My,Mz加载到这两个主节点上。
9.根据权利要求8所述风机偏航制动盘建模与极限强度校核方法,其特征在于,将塔顶中心载荷3个力分量Fx、Fy、Fz与2个力矩Mx、My加载到与主机架相连的主节点上,将小于偏航制动卡钳所能提供的最大制动力矩的载荷分量Mz加载到与偏航制动卡钳相连的主节点上,将超出部分也加载到与主机架相连的主节点上,从而使得塔顶中心载荷通过刚性梁传递至主机架和偏航制动卡钳,最终传递至偏航制动盘,以此来模拟偏航制动盘的受力情况。
10.根据权利要求2所述风机偏航制动盘建模与极限强度校核方法,其特征在于,所述步骤5)在强度结果评估时,查看各极限工况下偏航制动盘应力大小及危险位置,排除偏航制动盘与偏航轴承外圈及塔顶法兰之间的绑定区域,因绑定接触导致该区域结果不准确,将应力最大值与部件许用应力作比较,评估偏航制动盘的最小极限安全系数时需考虑材料安全系数1.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110783542.0A CN113627046B (zh) | 2021-07-12 | 2021-07-12 | 一种风机偏航制动盘建模与极限强度校核方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110783542.0A CN113627046B (zh) | 2021-07-12 | 2021-07-12 | 一种风机偏航制动盘建模与极限强度校核方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113627046A true CN113627046A (zh) | 2021-11-09 |
CN113627046B CN113627046B (zh) | 2023-08-08 |
Family
ID=78379518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110783542.0A Active CN113627046B (zh) | 2021-07-12 | 2021-07-12 | 一种风机偏航制动盘建模与极限强度校核方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113627046B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002070917A1 (en) * | 2001-03-01 | 2002-09-12 | Novatech Measurements Ltd | Braking force measurement |
JP2003337083A (ja) * | 2002-03-15 | 2003-11-28 | Ricoh Co Ltd | 回転駆動機構の解析方法、記憶媒体及び回転駆動機構の解析システム |
WO2010066163A1 (zh) * | 2008-12-11 | 2010-06-17 | 三一电气有限责任公司 | 风机偏航试验机构及风机整机试验台 |
US20100236875A1 (en) * | 2009-03-20 | 2010-09-23 | Bendix Spicer Foundation Brake Llc | Air Disk Brake Caliper Pre-Stressing Method and Pre-Stressed Caliper Apparatus |
FR2959717A1 (fr) * | 2010-05-05 | 2011-11-11 | Peugeot Citroen Automobiles Sa | Dispositif et procede pour la supervision de l'etat d'endommagement d'un disque de frein |
CN103440386A (zh) * | 2013-09-06 | 2013-12-11 | 南车株洲电力机车研究所有限公司 | 风力发电机偏航系统驱动力矩计算系统及其方法 |
DE102013013382A1 (de) * | 2013-08-10 | 2014-02-27 | Daimler Ag | Verfahren zur Fahrwerksoptimierung |
CN108547740A (zh) * | 2018-03-21 | 2018-09-18 | 内蒙古久和能源装备有限公司 | 一种风力发电机组主机架强度有限元计算方法 |
CN109492275A (zh) * | 2018-10-23 | 2019-03-19 | 许昌许继风电科技有限公司 | 一种风机主轴制动器强度校核方法及系统 |
JP2020056747A (ja) * | 2018-10-04 | 2020-04-09 | 公益財団法人鉄道総合技術研究所 | 鉄道車両用車輪の踏面損傷推定方法 |
CN111444647A (zh) * | 2020-03-13 | 2020-07-24 | 许昌许继风电科技有限公司 | 一种风机变桨轴承连接螺栓建模与强度校核方法 |
CN111680363A (zh) * | 2020-06-09 | 2020-09-18 | 南方天合底盘系统有限公司 | 一种用于浮动式鼓式制动器行车制动效能因数的计算方法 |
CN112240352A (zh) * | 2019-07-16 | 2021-01-19 | 广州汽车集团股份有限公司 | 一种轮毂轴承的有限元分析方法及其系统 |
WO2021114943A1 (zh) * | 2019-12-11 | 2021-06-17 | 湘电风能有限公司 | 风力发电机组偏航制动系统及其控制方法 |
-
2021
- 2021-07-12 CN CN202110783542.0A patent/CN113627046B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002070917A1 (en) * | 2001-03-01 | 2002-09-12 | Novatech Measurements Ltd | Braking force measurement |
JP2003337083A (ja) * | 2002-03-15 | 2003-11-28 | Ricoh Co Ltd | 回転駆動機構の解析方法、記憶媒体及び回転駆動機構の解析システム |
WO2010066163A1 (zh) * | 2008-12-11 | 2010-06-17 | 三一电气有限责任公司 | 风机偏航试验机构及风机整机试验台 |
US20100236875A1 (en) * | 2009-03-20 | 2010-09-23 | Bendix Spicer Foundation Brake Llc | Air Disk Brake Caliper Pre-Stressing Method and Pre-Stressed Caliper Apparatus |
FR2959717A1 (fr) * | 2010-05-05 | 2011-11-11 | Peugeot Citroen Automobiles Sa | Dispositif et procede pour la supervision de l'etat d'endommagement d'un disque de frein |
DE102013013382A1 (de) * | 2013-08-10 | 2014-02-27 | Daimler Ag | Verfahren zur Fahrwerksoptimierung |
CN103440386A (zh) * | 2013-09-06 | 2013-12-11 | 南车株洲电力机车研究所有限公司 | 风力发电机偏航系统驱动力矩计算系统及其方法 |
CN108547740A (zh) * | 2018-03-21 | 2018-09-18 | 内蒙古久和能源装备有限公司 | 一种风力发电机组主机架强度有限元计算方法 |
JP2020056747A (ja) * | 2018-10-04 | 2020-04-09 | 公益財団法人鉄道総合技術研究所 | 鉄道車両用車輪の踏面損傷推定方法 |
CN109492275A (zh) * | 2018-10-23 | 2019-03-19 | 许昌许继风电科技有限公司 | 一种风机主轴制动器强度校核方法及系统 |
CN112240352A (zh) * | 2019-07-16 | 2021-01-19 | 广州汽车集团股份有限公司 | 一种轮毂轴承的有限元分析方法及其系统 |
WO2021114943A1 (zh) * | 2019-12-11 | 2021-06-17 | 湘电风能有限公司 | 风力发电机组偏航制动系统及其控制方法 |
CN111444647A (zh) * | 2020-03-13 | 2020-07-24 | 许昌许继风电科技有限公司 | 一种风机变桨轴承连接螺栓建模与强度校核方法 |
CN111680363A (zh) * | 2020-06-09 | 2020-09-18 | 南方天合底盘系统有限公司 | 一种用于浮动式鼓式制动器行车制动效能因数的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113627046B (zh) | 2023-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109299532B (zh) | 一种风机主机架与后机架连接螺栓强度校核方法及系统 | |
CN107145663B (zh) | 车轮多目标优化设计方法 | |
CN111444647A (zh) | 一种风机变桨轴承连接螺栓建模与强度校核方法 | |
CN102945307B (zh) | 汽车底盘关键结构件结构优化设计方法 | |
CN105022868B (zh) | 风机偏航轴承连接螺栓极限强度和疲劳强度的计算方法 | |
De et al. | Finite element analysis on squeal-noise in railway applications | |
CN102880764B (zh) | 一种风力发电机组前机架强度的计算方法 | |
CN104965950B (zh) | 风机后机架极限强度和疲劳强度的计算方法 | |
CN102609578B (zh) | 一种汽车钢板弹簧多体模型的构建方法 | |
CN106547963B (zh) | 一种风机偏航制动盘极限强度的校核方法及装置 | |
Xiao et al. | Mathematical modeling and dynamic analysis of planetary gears system with time-varying parameters | |
CN113627046B (zh) | 一种风机偏航制动盘建模与极限强度校核方法 | |
CN105466371A (zh) | 测量飞机起落架轮轴端位置的装置及测量方法 | |
Dukalski et al. | Analysis of the influence of assembly electric motors in wheels on behaviour of vehicle rear suspension system | |
CN106886631A (zh) | 一种风机主轴性能的计算方法 | |
CN109086465B (zh) | 一种风机轮毂与主轴连接螺栓的强度校核方法 | |
Leo et al. | Drop test simulation for an aircraft landing gear via multi-body approach | |
Kanchwala et al. | Pitch reduction and traction enhancement of an EV by real-time brake biasing and in-wheel motor torque control | |
CN106289783A (zh) | 复现风力机6自由度载荷的液压加载试验台及液压加载策略 | |
CN108052753B (zh) | 基于doe与mbd的基础制动装置缓解性能优化方法 | |
CN108215692B (zh) | 自动引导车的悬架以及确定悬架参数的方法 | |
Chung et al. | Development of a systematic process for hydraulic brake system design and for hot judder characteristic estimation | |
CN112560162A (zh) | 一种动力总成悬置路谱载荷的缩减方法及装置 | |
CN110232218B (zh) | 一种风力发电机组叶片锁定装置强度的计算方法及装置 | |
CN116579213A (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 |