下载文档原格式
风力发电机叶片等效载荷计算及载荷谱分析米良;聂国林;程珩【摘要】风力发电机叶片部位通常受到随机变幅载荷的作用,所受随机载荷的随机性和无序性给载荷数据的处理带来了很大的困难.当前的数据处理方法通常是将随机变幅载荷转化为等效恒幅载荷进行分析,但由于其未能考虑低于疲劳极限的载荷对疲劳损伤所产生的影响故而会产生较大误差.针对上述问题,提出一种基于模糊理论的等效载荷计算方法,引入恰当的隶属函数,充分考虑低于疲劳极限的载荷对疲劳损伤所造成的影响,更加符合实际情况,以期提高等效载荷的计算精度.%The blade of wind turbine is usually subjected to random variable amplitudeload,which makes it difficult to process the load data.The current data processing method usually transforms the random load into equivalent constant amplitude load to reduce the data processingcapacity.However,the current method of equivalent load calculation fails to take into account the influence on the fatigue life made by the stress amplhude below the convention fatigue limit.So there is a large error in the equivalent load calculation by the current method.In view of the problems above,It properly considers the effects on fatigue life caused by load stress amplitude below the fatigue limit and presents a method of equivalent toad calculation by introducing the appropriate membership function based on fuzzy theory,which is more close to the actual situation.Thus improves the accuracy of the equivalent load calculation.【期刊名称】《机械设计与制造》【年(卷),期】2017(000)003【总页数】3页(P141-143)【关键词】模糊理论;等效载荷;隶属函数;程序载荷谱【作者】米良;聂国林;程珩【作者单位】太原理工大学新型传感器与智能控制教育部与山西省重点实验室,山西太原030024;太原理工大学新型传感器与智能控制教育部与山西省重点实验室,山西太原030024;太原理工大学新型传感器与智能控制教育部与山西省重点实验室,山西太原030024【正文语种】中文【中图分类】TH16;TK83风力发电机叶片是风电机系统中的重要工作部件，承载了主要的风力载荷，最容易发生疲劳破坏。
螺栓有效载荷计算公式
螺栓有效载荷是指螺栓或螺钉能够承受的最大拉力或剪力。
它是工程设计中非常重要的一个指标，能够保证机械装置的安全运行。
螺栓有效载荷的计算公式是根据螺栓的材料、直径、螺纹规格以及紧固力等参数来确定的。
一般而言，螺栓的有效载荷计算公式可分为拉力和剪力两种情况。
对于拉力情况，螺栓的有效载荷计算公式为：
P = F / A
其中，P代表螺栓的有效载荷，F代表螺栓所承受的拉力，A代表螺栓的截面积。
对于剪力情况，螺栓的有效载荷计算公式为：
P = F / A_s
其中，P代表螺栓的有效载荷，F代表螺栓所承受的剪力，A_s代表螺栓的剪切截面积。
在实际应用中，为了保证螺栓的安全性，通常会对螺栓的有效载荷进行安全系数的调整。
安全系数可以根据具体的工程要求来确定，一般建议在设计时选择适当的安全系数，以确保螺栓的可靠性和稳定性。
螺栓有效载荷的计算公式在工程设计中起着重要的作用，它能够帮助工程师评估螺栓的承载能力，从而选择合适的螺栓规格和数量。
合理的螺栓设计不仅可以提高机械装置的安全性和可靠性，还能够减少材料的浪费，降低成本。
螺栓有效载荷计算公式是工程设计中不可或缺的一部分，它能够帮助工程师评估螺栓的承载能力，确保机械装置的安全运行。
在实际应用中，我们需要根据具体情况选择合适的螺栓规格和数量，并考虑安全系数的影响，以确保螺栓的可靠性和稳定性。
收稿日期:2007-06-20.基金项目:国家863计划资助项目(2006AA5Z429).作者简介:王湘明(1963-),男,辽宁沈阳人,副教授,硕士,主要从事风力发电技术等方面的研究.文章编号:1000-1646(2008)01-0042-04海上风力发电机组塔架海波载荷的分析王湘明,陈 亮,邓 英,王婀娜(沈阳工业大学风能技术研究所,沈阳110023)摘 要:由于海上风力发电机组的外界载荷条件比陆地上的风力发电机组更加复杂,因此在机组设计中要考虑到海上风机载荷条件.在研究海波性能和运动规律的基础上,给出了海水中的塔架载荷的计算公式;针对海上风力发电机组的结构特点,将研究的海波载荷计算方法用于风机的塔架载荷计算,以110MW 风力发电机组为例,对塔架所承受的海波载荷和风载荷进行了计算,经过对比分析可知,风载荷大于海波载荷并为同一数量级的载荷,因此,在设计中更应注意这两种载荷的耦合作用.关 键 词:海上风力发电机组;海波载荷;风轮载荷;海波性能;运动规律中图分类号:T K 81 文献标志码:AAnalysis on w ave loads of off shore wind turbineWAN G Xiang 2ming ,CHEN Liang ,DEN G Y ing ,WAN G E 2nuo(Wind Energy Institute of Technology ,Shenyang University of Technology ,Shenyang 110023,China )Abstract :The load condition of offshore wind turbine is more complicated than that on land ,and must be considered in design.In view of the wave feature and its moving rule ,the formulas of tower loads on offshore wind turbines were presented.Based on the structure feature of offshore wind turbines ,the computation method of the wave loads was applied to calculate the tower loads.The wave loads and wind loads of 110MW wind turbine on tower were obtained.By means of analysis and contrast ,it is concluded that wind loads are heavier than wave loads ,but both are in same order.Therefore ,the coupling action between two loads should be paid a more attention in design.K ey w ords :offshore wind turbine ;wave loads ;wind loads ;wave feature ;moving rule1 海上风力发电机组载荷计算通常海上风力发电机组的外界载荷条件比陆地上的风力发电机组更加复杂,其原因是除了受到通常的风载荷以外,更主要的是受到海波和海浪的影响,而海上风机主要是受风轮载荷的影响,此外还受到海浪和海波的影响,这样海水下面塔架的外界载荷就变得比较复杂.因此在机组设计中,除了要考虑陆上风机的风作用下的21种载荷工况之外,还要考虑到42种海上风机载荷条件[1].如果风机装在渤海湾的某一浅海滩上,如图1所示,设海平面为起始位置0,则塔架在海平面以下的部分为0到(-h ),在海面以上部分为0到h 0.塔架总载荷包括风轮载荷、风载载荷和海波载荷.F wind 表示风轮载荷,F h 0表示风载载荷,F h 表示海波载荷.以风机塔架计算为例,风轮载荷为F wind =12C t ρair v 2s A ・B ・S (1)式中:C t ———扭矩系数;ρair ———空气密度;v s ———风速;A ———风轮扫掠面积;B ———叶片的支数;S ———安全系数. 第30卷第1期2008年2月沈 阳 工 业 大 学 学 报Journal of Shenyang University of T echnologyVol 130No 11Feb 12008图1 海上风力发电机组结构图Fig 11 Structure of offshore wind turbine塔架的风载载荷为F ho =12ρair Aw v 2tw (h 0) (2)式中:A w ———风吹到塔架的投影面积;v tw (h )———海平面上高度为h 的平均风速.则塔架总风载载荷为F w ind +F ho ,用blade 软件对风机25m 风速时功率运行工况进行模拟得到塔架的总风载荷,如图2所示.图2 发电时塔架水平方向风载荷图Fig 12 Wind load of tow er in horizontaldirection under generating根据文献[1],海波载荷为F h =(C m -1)ρπD 24L v t +ρπD 24L v D +12C dρDL v t |v t |(3)式中:F h ———作用在直径为D ,长度为L 的截断面上的正常作用力;ρ———海水密度;C m ———惯性系数;C d ———阻力系数.则塔架所受的最大弯曲应力δ=F wind (h 0+h )+F h 0h 02+h +F hh2π(D 4-D 41)32D(4)式中:D ———塔架外径;D 1———塔架内径.塔架的风载荷计算直接参考机械设计手册[4],海流的载荷计算直接与海波条件相关,计算比较复杂,可按参考文献[2]的方法进行计算.2 海波载荷的计算211 海波载荷的计算公式在时间t 时海波中点(x ,y ,z )处的波粒子速度和加速度矢量分别用v D 和 v D 表示.相同地点的总的海流速度矢量用v c 表示,塔架结构本身的速度和加速度分别用v s 和 v s 表示.则海波相对于塔架结构流动的速度v t 和加速度 v t 分别为v t =v D +v c -v sv t = v D + v s 在塔架受海波载荷的研究中,如果塔架的衍射参数D/λ的值小于012(D 为塔架的直径,λ为海波的波长),塔架在水平方向上实际的海波载荷可用Morison 方程式(3)求得.要计算塔架作用力,可以把塔架近似为相等高度的10个子圆柱体.使用Morison 相对运动方程可计算子圆柱体轴上的作用力.把各个参数代入方程(3)就可以计算出塔架的海波载荷.阻力系数和惯性系数可以参考海上风力发电机组相关文献,也可通过海波实验确定,同时,还要考虑海波的结构动力的总衰减,总的衰减包括波辐射衰减、阻力衰减和结构衰减.212 计算公式系数的确定在方程(3)中,阻力系数C d 和惯性系数C m是未知参数,需要得到海波粒子的速度和加速度的测量值.在参考文献[3-5]中采用很多方法给出了不同C d 和C m 的数值,本文按期望的方差确定系数的准确值[6].在参考文献[3]中给出了确定阻力系数和惯性系数的两种方法.第一种方法是当流体的速度和加速度为零时解出莫里森方程.当速度为零时可以确定惯性系数C m ,当加速度为零时可以确定阻力系数C d .这种方法的缺点是只用到了少量的测量数据,并且34第1期王湘明,等:海上风力发电机组塔架海波载荷的分析 对相移误差比较敏感[4];第二种方法是基于最小二乘法,这样可以使压力的测量值和计算值的平方差ε2最小.为了消除计算误差,本文要用近似的雷诺数将计算值和测量值的数据进行分组,从而确定C m 和C d 的值.213 典型的系数值根据文献[4],阻力系数和惯性系数取值见表1.表1中的雷诺数为R n =v a Du (5)式中:v a ———水流速度的幅值,m/s ;u ———粘性系数,m 2/s .表1中的KC 为K eulegan Carpenter 数KC =v a TD(6)式中,T 为震荡周期.表1 系数值T ab 11 Coeff icient valueKCR n <105C dC m R n >105C d C m≤10112112210115≥10016016210115214 海波速度的计算海波包括深水海波和过渡海波,因此研究海波速度也应该从这两方面考虑.对于深水波周期T 应该满足条件d/(g T 2)>018;对于过渡海波周期T 应该满足条件01002<d/(g T 2)<0108.其中d 为水深,g 为重力加速度.起源波的参数包括波粒子的速度和加速度,在这里用线性波理论来定义这些参数[7].21411 深水区的波粒子速度首先建立一个x 2y 2z 右手坐标系统,用x 2y 座标轴所在的平面表示不受干扰的海平面,z 轴垂直向上,那么参考速度可以表示为v R =015H D ・ωD ・exp {i (ωD t -K D ξ)} (7)式中:ωD =2π/T D ,K D =ω2/g ,ξ=x ・cos μ+y ・sin μ,i =-1;ξ轴———海波沿地平线传播的方向;μ———x 轴和ξ轴之间的夹角.在海平面下坐标为(x ,y ,z )点的合成速度矢量v D 为v D =(v D x ,v D y ,v D z )=(i cos μ,i sin μ,-1)・v R ・exp {K D ・z}这样可以用相位纠正的方法来计算v D x/y/z (x/y/z指x 轴,y 轴或z 轴的方向),也就是实部和虚部的相移为T D /4时,实部在t =T/4时达到波峰.在一个海波周期中最大的速度幅度值为v D x/y/z max=[Re (v D x/y/z )]2+[Im (v D x/y/z )]2(8)21412 过渡海波的波粒子速度用参考速度v R 来定义过渡海波的参考速度为v Rs =v R (K Ds )/sin h (K Ds ・d ) (9)式中,K Ds 为分布方程得解,分布方程为K Ds ・|tanh {K Ds ・d}|=ω2/g (10) 因此,过渡海波在海面下的坐标为(x ,y ,z ),点的速度矢量为v Ds =(v Ds x ,v Ds y ,v Ds z )=(i cos h ξ,cos μ,i cos h ξ,sin μ,-sin ξ)・v Rs式中,ξ=K Ds ・(z +d )最大值为v Ds x/y/z m ax=[R e (v Ds x/y/z )]2+[Im (v Ds x/y/z )]2(11)21413 海波粒子的加速度深水区的波粒子加速度用a D 表示,过渡区波粒子的加速度用a Ds 表示,则a D =i ωD ・v D ,a Ds =i ωD ・v Ds (12)那么加速度的最大值a D x/y/z max =|ωD |・v D x/y/z max a Ds x/y/z max=|ωD |・v Ds x/y/z max(13)3 海波载荷计算举例下面以渤海湾实际数据为例来计算塔架单位长度的海波载荷.取T =9s ,D =8m ,v a =114m/s ,ρ=11022×103kg/m 3,C d =112,C m =210.将这些数据代入公式(3)中,用matlab 计算可得结果如图3所示.图3 塔架水平方向的海波载荷Fig 13 W ave load of tow er in horizontal direction由此看出塔架的海波载荷为周期性变化的载荷.以110MW 风力发电机组为例,取塔架根部外径为412m ,内径为41156m ,塔架海平面以上高44 沈 阳 工 业 大 学 学 报第30卷度为70m,海平面以下高度为10m,则得出塔架的最大应力为3217MPa.塔架海面上的风轮的载荷计算如图2所示,最大为21215kN.海面下的海波载荷如图3所示,最大为100kN.因此,海波载荷加上风的载荷等于塔架的总载荷.从计算结果可知,海波载荷对塔架的载荷影响很大,设计时还应特殊考虑海波极端工况的塔架载荷.4 结 论1)对海浪和海波的运动规律进行了描述,给出了海波载荷的计算公式和计算结果,同时,也给出了塔架风载荷的计算公式和计算结果.2)假设110MW风机装在渤海湾近海的海滩上,按照matlab软件计算结果(见图2),塔架在海平面以下的海波载荷为100kN,且为周期性载荷,当海流速度比较大时海波载荷至少为正常情况时的2倍.同样,根据blade软件计算塔架的风载荷(见图3)为塔架所受阻力,正常发电时的海面以上塔架载荷为120kN,发生最大风速扰动时至少也为2倍以上.所以,在设计时,应重点考虑风流和海流的极端载荷工况.3)采用以上的计算公式,对110MW风机为例计算了海波载荷和风载荷,经过对比分析可知,风载荷与海波载荷是同一数量级的载荷.它们之间的耦合作用较大,将会影响机组的稳定性运行.参考文献(R eference):[1]全国风力机械标准化技术委员会.风力机械标准汇编[M].北京:中国标准出版社,2006.(Wind Power Machine Standard Institute.Wind power machine criterion assembly[M].Beijing:Standard Press of China,2006.)[2]Peeringa J M.Wave loads on offshore wind turbines[J].E leclronic C ommunication Network,2004(4):7-10.[3]Dean R G,Dalrympl R A.Water wave mechanics forengineers and scientists[M].India:World Scientific Publishing Co Ltd,1991.[4]Journée.O ffshore hydromechanics[D].Delft:Delft Universi2ty of T echnology,2001.[5]Sarpkaya T,Isaacson I.Mechanics of wave forces onoffshore structures[J].Electronic Communication Net2 work,2004(4):17-19.[6]Grundlehner G J.Systematic model tests on a harsh en2vironment jack2in elevated conditions[C]//Procee2 dings of5th International Conference on Jack2up Plat2 forms.City University,London,1995.[7]Frost W.Regulations for the certification of offshorewind Energy conversion systems[J].Wind Engineer2 ing,1994(1):11-15.[8]徐灏,邱宣怀,蔡春源,等.机械设计手册[M].北京:机械工业出版社,1992.(XU Hao,Q IU Xuan2huai,CAI Chun2yuan,et al.Ma2 chine design manual[M].Beijing:China Machine Press,1992.)(责任编辑:邓美艳 英文审校:杨俊友)(上接第23页)[5]Bertotti G.Dynamic generalization of the scalar preisachmodel of hysteresis[J].IEEE Trans on Magnetics, 1992,28:2592-2601.[6]Mayergoyz I D.Vector preisach hysteresis models[J].Appl Phys,1988,63(8):2995-3000.[7]Nakata T,Takahashi N,K awase Y.Finite element ana2lysis of magnetic fields taking into account hysteresis characteristics[J].IEEE Trans on Magnetics,1985,21(5):1856-1858.[8]王小平.遗传算法———理论、应用与软件实现[M].西安:西安交通大学出版社,2002.(WAN G Xiao2ping.G enetic algorithms—theories,appli2 cation and software implementation[M].Xi’an:Xi’an Jiaotong University Press,2002.)[9]Kumada M.Development of a model4tesla permanentmagnet[C]//Proceedings of Particle Accelerator Con2 ference.Chicago,2001:3221-3223.(责任编辑:邓美艳 英文审校:杨俊友)54第1期王湘明,等:海上风力发电机组塔架海波载荷的分析 。
兆瓦级风力发电机主要承载螺栓强度的有限元计算沃晓临;马春翔【摘要】兆瓦级风力发电机承受载荷复杂,受风速变化及本身支撑结构的影响较大,载荷有一定的随机性,因此在设计中必须考虑多自由度上多个载荷的叠加作用,传统工程算法一般只能考虑螺栓连接处的单轴受载情况已经不再适应.为了尽量模拟实际受载情况,采用了ANSYS对兆瓦级风力发电机轮毂与主轴的连接螺栓做有限元分析,载荷基于动力学计算软件计算得出的时间序列,对连接螺栓本身,在极限载荷与疲劳载荷下进行了仿真计算.【期刊名称】《电力与能源》【年(卷),期】2011(032)003【总页数】4页(P228-231)【关键词】有限元仿真;螺栓强度计算;多轴疲劳计算【作者】沃晓临;马春翔【作者单位】上海交通大学机械与动力工程学院,上海200240;上海电气风电设备有限公司,上海200030;上海交通大学机械与动力工程学院,上海200240【正文语种】中文【中图分类】TK83连接螺栓作为关系到风机结构安全的关键零件之一,需结合风机载荷特点,采用可靠、准确的计算方法对其校核,以提高风机整体可靠性,保证其在复杂外载荷下20年的工作寿命。
风力机的风轮所受载荷包括风载和重力载荷,受载情况十分复杂。
其中影响风轮风载的因素有垂直扫风面的稳态气流、存在偏航误差的气流、主轴的倾角、风切变、塔影效应、尾流影响等;影响风轮的重力载荷的因素有离心力、偏航动作时的陀螺载荷、风轮制动时的制动载荷、叶片挥舞带来的载荷等。
由于载荷的复杂性,只考虑单一倾覆力矩或单一工作载荷的常用工程算法或准则计算方法,已经不能直接应用于风电行业。
本文将结合风力机设计载荷的特点,依据德国船级社风力机认证导则GL2003的要求,使用ANSYS12.0软件,用数值仿真的方法,对兆瓦级风力发电机组的主轴与轮毂的连接螺栓进行计算。
利用数值仿真的结果,对采用高强度螺栓工程计算标准VDI2230计算所得的计算应力进行修正,以得到更接近实际的计算值。
