风机叶片覆冰的机组载荷的计算方法与流程

本发明实施例涉及风电
技术领域：
，尤其涉及一种风机叶片覆冰的机组载荷的计算方法。
背景技术：
：目前，风力发电机组叶片覆冰的机组载荷仿真计算一般采用gl2010规范中给出的相关方法，gl2010规范采用的覆冰质量模型如图1所示。依据叶尖外形插值计算叶尖最小结冰弦长cmin，进而可以进一步计算出叶片覆冰的最大线密度μe。叶片覆冰的最大线密度μe的计算公式如下所示：μe＝ρekcmin(cmin+cmax)其中，冰密度ρe＝700kg/m3；系数k＝0.00675+0.3exp(-0.32r)，cmax为最大弦长。从图1所示的覆冰质量模型可知，由风轮中心o至风轮半径中间距离s的覆冰线密度由零线性增加至最大线密度值μe，之后保持最大线密度μe值至叶尖b。根据风机叶片的这种覆冰质量模型结合冰密度可以计算出风机叶片的覆冰质量。然而，依据gl2010规范中的叶片覆冰质量计算方法得出的覆冰质量偏保守，导致冰载极限、疲劳载荷偏大，且未考虑覆冰对叶片气动性能的影响，计算结果不准确。iec61400-1标准第四版提出了一种叶片覆冰质量计算的新方法。iec61400-1标准第四版采用的覆冰质量模型如图2所示。各截面处叶片的覆冰线密度m(r)由风轮中心o至叶尖b线性增加至m。各截面处叶片的覆冰线密度m(r)的计算公式如下所示：m(r)＝0.125×c85×r其中，c85为85％叶轮半径位置处的弦长，r为叶片各截面距轮毂中心o的距离。另外，由于覆冰后叶片的气动性能衰减，iec61400-1标准第四版还明确了覆冰翼型气动性能参数(升阻力系数)修正方法——惩罚因子。并且，不同载荷分析类型，例如极限载荷和疲劳载荷的叶片覆冰质量的考虑方式不同，如下表所示：分析类型叶片1叶片2叶片3惩罚因子极限载荷2×m(r)2×m(r)-是疲劳载荷m(r)m(r)m(r)×0.5是iec61400-1标准第四版提出的叶片覆冰质量计算方法及覆冰翼型气动性能参数(升阻力系数)修正方法更加合理。但是，由于bladed载荷仿真计算软件中默认覆冰质量模型采用gl2010规范的设定，所以，目前无法在bladed载荷仿真计算软件中实现iec61400-1标准第四版的覆冰的建模，进而无法开展准确合理的冰载工况假定及计算分析工作。技术实现要素：本发明实施例的目的在于提供一种风机叶片覆冰的机组载荷的计算方法，能够更加精确地计算出风机叶片覆冰后的机组载荷。本发明实施例的一个方面提供一种风机叶片覆冰的机组载荷的计算方法，采用基于风机的第一标准的载荷仿真计算软件，所述第一标准以风机叶片的第一覆冰质量模型建模。所述方法包括：基于第二标准来修正所述载荷仿真计算软件中所述风机叶片的翼型，其中，所述第二标准以风机叶片的第二覆冰质量模型建模；基于所述第一标准计算所述风机叶片的第一覆冰质量矩；基于所述第二标准计算所述风机叶片的第二覆冰质量矩；找出计算出的所述第一覆冰质量矩与所述第二覆冰质量矩相等时的所述风机叶片的最小叶尖结冰弦长；将找出的所述最小叶尖结冰弦长输入到所述载荷仿真计算软件中，基于所述第二标准在所述载荷仿真计算软件中建立所述风机叶片的覆冰模型；以及基于所述第二标准的覆冰工况通过所述载荷仿真计算软件来计算获得所述风机叶片覆冰的机组载荷。进一步地，所述载荷仿真计算软件包括bladed软件。进一步地，所述风机的第一标准包括gl2010规范，所述风机的第二标准包括iec61400-1标准第四版。进一步地，所述基于第二标准来修正所述载荷仿真计算软件中所述风机叶片的翼型包括：基于所述第二标准来修正所述载荷仿真计算软件中所述风机叶片的气动性能参数。进一步地，所述基于所述第二标准来修正所述载荷仿真计算软件中所述风机叶片的气动性能参数包括：使用惩罚因子来修正所述风机叶片的气动性能参数。进一步地，所述使用惩罚因子来修正所述风机叶片的气动性能参数包括：最小攻角至失速攻角范围内的攻角，分别使用升力系数惩罚因子和阻力系数惩罚因子来获得所述风机叶片修正后的升力系数和阻力系数；及在-180度至180度范围内除最小攻角至失速攻角范围之外的攻角，使用外推方法来获得所述风机叶片修正后的升力系数和阻力系数。进一步地，所述基于所述第一标准计算所述风机叶片的第一覆冰质量矩包括：基于所述第一覆冰质量模型中叶片的第一覆冰线密度来计算所述叶片的所述第一覆冰质量矩。进一步地，所述第二覆冰质量矩包括第二极限单叶片覆冰质量矩，所述基于第二标准计算所述风机叶片的第二覆冰质量矩包括：基于所述第二覆冰质量模型中叶片的第二覆冰线密度的2倍来计算所述风机叶片的所述第二极限单叶片覆冰质量矩。进一步地，所述找出计算出的所述第一覆冰质量矩与所述第二覆冰质量矩相等时的所述风机叶片的最小叶尖结冰弦长包括：找出计算出的所述第一覆冰质量矩与所述第二极限单叶片覆冰质量矩相等时的所述风机叶片的最小叶尖结冰弦长，所述建立所述风机叶片的覆冰模型包括：建立所述风机叶片的极限覆冰模型，所述计算获得所述风机叶片的覆冰载荷包括：计算获得所述风机叶片的极限覆冰的机组载荷。进一步地，所述第二覆冰质量矩包括第二疲劳单叶片覆冰质量矩，所述基于第二标准计算所述风机叶片的第二覆冰质量矩包括：基于所述第二覆冰模型中叶片的第二覆冰线密度的0.5倍来计算所述风机叶片的所述第二疲劳单叶片覆冰质量矩。进一步地，所述找出计算出的所述第一覆冰质量矩与所述第二覆冰质量矩相等时的所述风机叶片的最小叶尖结冰弦长包括：找出计算出的所述第一覆冰质量矩与所述第二疲劳单叶片覆冰质量矩相等时的所述风机叶片的最小叶尖结冰弦长，所述建立所述风机叶片的覆冰模型包括：建立所述风机叶片的疲劳覆冰模型，所述计算获得所述风机叶片的覆冰载荷包括：计算获得所述风机叶片的疲劳覆冰的机组载荷。进一步地，所述建立所述风机叶片的覆冰模型包括：在所述载荷仿真计算软件中将风机叶片的原线密度增加所述第二覆冰线密度的0.5倍。本发明实施例的风机叶片覆冰的机组载荷的计算方法采用基于风机的第一标准的载荷仿真计算软件，实现了基于第二标准的叶片覆冰质量计算方法及基于第二标准的风机叶片的翼型修正方法在该载荷仿真计算软件中的建模。从而，能够更加精确地计算出风机叶片覆冰后的机组载荷。附图说明图1为gl2010规范中采用的覆冰质量模型的示意图；图2为iec61400-1标准第四版中采用的覆冰质量模型的示意图；图3为本发明一个实施例的风机叶片覆冰的机组载荷的计算方法的流程图；图4为一种风机叶片的示意图。具体实施方式这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施例并不代表与本发明相一致的所有实施例。相反，它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置的例子。在本发明实施例使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。除非另作定义，本发明实施例使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本发明说明书以及权利要求书中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。同样，“一个”或者“一”等类似词语也不表示数量限制，而是表示存在至少一个。“多个”或者“若干”表示两个及两个以上。除非另行指出，“前部”、“后部”、“下部”和/或“上部”等类似词语只是为了便于说明，而并非限于一个位置或者一种空间定向。“包括”或者“包含”等类似词语意指出现在“包括”或者“包含”前面的元件或者物件涵盖出现在“包括”或者“包含”后面列举的元件或者物件及其等同，并不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而且可以包括电性的连接，不管是直接的还是间接的。在本发明说明书和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。本发明实施例提出了一种风机叶片覆冰的机组载荷的计算方法，该计算方法采用基于风机的第一标准的载荷仿真计算软件，实现了基于第二标准的叶片覆冰质量计算方法及基于第二标准的风机叶片的翼型修正方法在该载荷仿真计算软件中的建模。其中，第一标准以风机叶片的第一覆冰质量模型建模，第二标准以风机叶片的第二覆冰质量模型建模。图1揭示了本发明一个实施例的风机叶片覆冰的机组载荷的计算方法的流程图。如图1所示，本发明一个实施例的风机叶片覆冰的机组载荷的计算方法可以包括步骤s11-s16。在步骤s11中，基于第二标准来修正载荷仿真计算软件中风机叶片的翼型。在步骤s12中，基于第一标准计算风机叶片的第一覆冰质量矩。在步骤s13中，基于第二标准计算风机叶片的第二覆冰质量矩。在步骤s14中，找出计算出的第一覆冰质量矩与第二覆冰质量矩相等时的风机叶片的最小叶尖结冰弦长tipchord。在步骤s15中，将找出的最小叶尖结冰弦长tipchord输入到该载荷仿真计算软件中，基于第二标准在该载荷仿真计算软件中建立风机叶片的覆冰模型。在步骤s16中，基于第二标准的覆冰工况通过该载荷仿真计算软件来计算获得风机叶片覆冰的机组载荷。如图3所示，本发明实施例的风机叶片覆冰的机组载荷的计算方法包括以上多个列举的步骤，但是，本发明实施例的风机叶片覆冰的机组载荷的计算方法还可以在列举的步骤之前、之后以及之间包括额外的步骤。在一些实施例中，列举的步骤中的一个或多个步骤可以不同于图3所示的顺序执行。在一个实施例中，本发明实施例的载荷仿真计算软件可以包括bladed软件，风机的第一标准可以包括gl2010规范，风机的第二标准可以包括iec61400-1标准第四版。其中，gl2010规范采用图1所示的第一覆冰质量模型，iec61400-1标准第四版采用图2所示的第二覆冰质量模型。然而，本发明实施例并不局限于此。在一些实施例中，步骤s11的基于第二标准来修正载荷仿真计算软件中风机叶片的翼型包括：基于第二标准来修正载荷仿真计算软件中风机叶片的气动性能参数。例如，在第二标准为iec61400-1标准第四版的实施例中，可以使用惩罚因子来修正风机叶片的气动性能参数。以下将详细介绍如何基于iec61400-1标准第四版使用惩罚因子来修正风机叶片的气动性能参数。在最小攻角α0至失速攻角αs范围内的攻角，即α0～αs，可以分别使用升力系数惩罚因子和阻力系数惩罚因子来获得风机叶片修正后的升力系数和阻力系数。惩罚因子公式如下所示：cl，pen(α)＝-0.0014α2-0.0017α+0.9509cd，pen(α)＝0.0191α+3.1151其中，cl，pen为升力系数惩罚因子，cd，pen为阻力系数惩罚因子，α为攻角。则风机叶片修正后的升力系数和阻力系数即覆冰翼型升阻力系数可以分别表示为：cl，iced＝cl，pen×cl，cleancd，iced＝cd，pen×cd，clean其中，cl，iced为覆冰翼型升力系数，cl，clean为翼型原升力系数；cd，iced为覆冰翼型阻力系数，cd，clean为翼型原阻力系数。因此，在α0～αs范围的攻角，可以根据风机叶片翼型的原升力系数cl，clean和原阻力系数cd，clean分别与升力系数惩罚因子cl，pen和阻力系数惩罚因子cd，pen的乘积来得到风机叶片修正后(即叶片覆冰后)的升力系数cl，iced和阻力系数cd，iced。在-180度至180度范围内除最小攻角α0至失速攻角αs范围之外的攻角，可以使用外推方法，例如viterna方法外推来获得风机叶片修正后的升力系数和阻力系数，弯矩系数cm不做改动。具体如下所述。在αs～90°攻角范围内，风机叶片修正后的升力系数cl，iced和阻力系数cd，iced可由下式得出：cl，iced＝cdmaxsin(2α)/2+a2cos2(α)/sin(α)cd，iced＝cdmaxsin2(α)+b2cos(α)其中，cdmax表示攻角α＝90°时的阻力系数，cdmax可以根据如下公式得出：cdmax＝min(1.11+0.018ar，2.01)其中，ar为叶片展弦比，代表风轮半径与80％风轮半径处弦长的比值。a2＝(cls-cdmaxsin(αs)cos(αs))sin(αs)/cos2(αs)b2＝(cds-cdmaxsin2(αs))/cos(αs)其中，cls为失速攻角αs处对应的升力系数；cds为失速攻角αs处对应的阻力系数。对于90°～180°-αs、-180°+αs～-αs攻角范围，翼型的升力系数可以通过缩放和镜像αs～90°的升力系数曲线获得。对于不对称翼型，缩放系数cladj为0.7。攻角为±180°时，则升力系数为0。其它升力系数曲线所缺部分(-180°～-180°+αs、-αs～α0、180°-αs～180°)，可以通过线性插值获得。翼型的阻力系数不缩放，可以直接镜像获得。因此，在-180°～-180°+αs攻角范围内，风机叶片修正后的升力系数cl，iced和阻力系数cd，iced如下所示：cl，iced＝(180+α)/αsclscladjcd，iced＝cdmaxsin2(180+α)+b2cos(180+α)在-180°+αs～-90°攻角范围内，风机叶片修正后的升力系数cl，iced和阻力系数cd，iced如下所示：cl，iced＝cladj(cdmaxsin(2(α+180))/2+a2cos2(180+α)/sin(180+α))cd，iced＝cdmaxsin2(180+α)+b2cos(180+α)在-90°～-αs攻角范围内，风机叶片修正后的升力系数cl，iced和阻力系数cd，iced如下所示：cl，iced＝cladj(-cdmaxsin(2(-α))/2-a2cos2(-α)/sin(-α))cd，iced＝cdmaxsin2(-α)+b2cos(-α)在-αs～α0攻角范围内，风机叶片修正后的升力系数cl，iced和阻力系数cd，iced如下所示：cl，iced＝-clscladj+(α+αs)/(α0+αs)(clscladj+cl0)cd，iced＝cd0+(-α+α0)(cds–cd0)/(αs+α0)其中，cl0为最小攻角α0处对应的升力系数；cd0为最小攻角α0处对应的阻力系数。cl0和cd0可以由上面提到的惩罚因子修正方法得出。在90°～180°-αs攻角范围内，风机叶片修正后的升力系数cl，iced和阻力系数cd，iced如下所示：cl，iced＝cladj(-cdmaxsin(2(180-α))/2-a2cos2(180-α)/sin(180-α))cd，iced＝cdmaxsin2(180-α)+b2cos(180-α)在180°-αs～180°攻角范围内，风机叶片修正后的升力系数cl，iced和阻力系数cd，iced如下所示：cl，iced＝(α-180)/αsclscladjcd，iced＝cdmaxsin2(α-180)+b2cos(α-180)在一些实施例中，步骤s12的基于第一标准计算风机叶片的第一覆冰质量矩包括：基于第一覆冰质量模型中叶片的第一覆冰线密度来计算叶片的第一覆冰质量矩。以下将以第一标准为gl2010规范为例来详细介绍如何基于采用图1所示的第一覆冰质量模型的gl2010规范来计算叶片的第一覆冰质量矩。(a1)叶轮半径：r＝0.5d，其中d为叶轮直径；(b1)叶轮半径一半位置：s＝0.5r；(c1)叶尖最小结冰弦长：cmin＝tipchord，tipchord表示叶尖最小结冰弦长，此值与叶尖形状有关，所以，在叶片选定后，叶尖形状就不变，此值即为常数，其中tipchord为bladed软件界面中输入的参数；(d1)叶片的第一覆冰线密度即覆冰单位长度质量μe如下所示：μe＝ρekcmin(cmin+cmax)其中：冰密度ρe＝700kg/m3；系数k＝0.00675+0.3exp(-0.32r)，cmax为最大弦长；(e1)叶片各截面距轮毂中心o的距离数组：rr＝r+l，其中，r为叶片各截面距叶根a的距离数组，l为叶根a距轮毂中心o的距离；(f1)将位置s放入数组rr，并对数组rr从小到大排列，新排列的数组为rrr；例如，以某叶片模型为例说明：r可从叶片模型中获取，r为叶片各截面距叶根a的距离数组，为一系列离散数值，如：0、0.4、1.4、2.4、3.4……70.67、71.69，l为固定数值，从模型中读取为1.668，因为rr＝r+l，所以rr为1.668、2.068、3.068、4.068、5.068……72.338、73.348，s＝0.5r＝0.5×73.348＝36.674，然后，将s与数组rr合并为：1.668、2.068、3.068、4.068、5.068……72.338、73.348，36.674，排序得到新数组rrr：1.668、2.068、3.068、4.068、5.068……36.674(s)……72.338、73.348。(g1)加入位置s后，叶片各截面距叶根a的距离数组q＝rrr-l；(h1)数组rrr中叶根a位置到位置s的覆冰单位长度质量数组为：μe1＝rrr1μe/rrr1(s)图4揭示了一种叶片的示意图。如图4所示，图4中的1至n分别表示叶片的n个截面。在上式中，rrr1为位置1到位置s的叶片各截面距轮毂中心o的距离数组，其中，位置1代表叶根a位置，即数组rrr中的第一个数；rrr1(s)为位置s距轮毂中心o的距离；(i1)数组rrr中位置s+1到叶尖b位置的覆冰单位长度质量数组为：μe2＝μe(j1)叶根a位置到叶尖b位置覆冰单位长度质量数组为：mgl＝[μe1，μe2](k1)覆冰单位长度质量矩数组为：igl＝mgl×q(l1)igl是叶片各截面的覆冰质量矩，为一系列的离散数值。由于叶片是一个连续的实物，因此，要想通过覆冰质量矩一系列的离散值igl来获得连续的覆冰质量矩，需要对igl进行积分。叶片的第一覆冰质量矩可以表示为：igl＝∫igldq因此，通过上述步骤，可以获得叶片的第一覆冰质量矩igl。在一些实施例中，第二覆冰质量矩包括第二极限单叶片覆冰质量矩。因此，步骤s13的基于第二标准计算风机叶片的第二覆冰质量矩包括：基于第二覆冰质量模型中叶片的第二覆冰线密度的2倍来计算风机叶片的第二极限单叶片覆冰质量矩。以下将以第二标准为iec61400-1标准第四版为例来详细介绍如何基于采用图2所示的第二覆冰质量模型的iec61400-1标准第四版来计算叶片的第二极限单叶片覆冰质量矩。(a2)叶轮半径：r＝0.5d，其中d为叶轮直径；(b2)85％叶轮半径位置：c85_p＝0.85r；(c2)叶片各截面距轮毂中心o的距离数组：rr＝r+l，其中r为叶片各截面距叶根a的距离数组，l为叶根a距轮毂中心o的距离；(d2)85％叶轮半径位置处弦长：c85＝(c85_p–rr(t))(c(t+1)–c(t))/(rr(t+1)-rr(t))+c(t)其中，rr(t)代表85％叶轮半径位置前的位置，rr(t+1)代表85％叶轮半径位置后的位置，c(t)代表85％叶轮半径位置前的位置对应的弦长，c(t+1)代表85％叶轮半径位置后的位置对应的弦长；(e2)叶片的第二覆冰线密度即覆冰单位长度质量m1如下所示：m1＝2×0.125×c85×r(f2)叶根a位置到叶尖b位置的覆冰单位长度质量数组为：miec1＝rrm1/rr(tip)其中，rr(tip)为叶尖b距轮毂中心o的距离；(g2)覆冰单位长度质量矩数组为：iiec1＝miec1×r(h2)iiec1是叶片各截面的覆冰质量矩，为一系列的离散数值。由于叶片是一个连续的实物，因此，要想通过覆冰质量矩一系列的离散值iiec1来获得连续的覆冰质量矩，需要对iiec1进行积分。因此，第二极限单叶片覆冰质量矩可以表示为：iiec1＝∫iiec1dr因此，通过上述步骤，可以获得叶片的第二极限单叶片覆冰质量矩iiec1。在计算极限工况时，步骤s14的找出计算出的第一覆冰质量矩与第二覆冰质量矩相等时的风机叶片的最小叶尖结冰弦长包括：找出计算出的第一覆冰质量矩igl与第二极限单叶片覆冰质量矩iiec1相等时的风机叶片的最小叶尖结冰弦长tipchord。相应地，步骤s15中的建立风机叶片的覆冰模型包括：建立风机叶片的极限覆冰模型。步骤s16中的计算获得风机叶片的覆冰载荷包括：计算获得风机叶片的极限覆冰的机组载荷。具体地，将找出的第一覆冰质量矩igl与第二极限单叶片覆冰质量矩iiec1相等时的最小叶尖结冰弦长tipchord值输入到bladed软件中，勾选两叶片结冰，输入冰密度ρe＝700kg/m3。保存模型，完成极限载荷分析覆冰建模。依据建立的极限覆冰bladed模型，按照iec61400-1标准第四版的覆冰工况要求进行工况假定及bladed仿真计算，并对载荷计算结果进行后处理。从而可以最终获得风机叶片的极限覆冰的机组载荷。在另一些实施例中，第二覆冰质量矩包括第二疲劳单叶片覆冰质量矩。因此，步骤s13的基于第二标准计算风机叶片的第二覆冰质量矩包括：基于第二覆冰模型中叶片的第二覆冰线密度的0.5倍来计算风机叶片的第二疲劳单叶片覆冰质量矩。以下将以第二标准为iec61400-1标准第四版为例来详细介绍如何基于采用图2所示的第二覆冰质量模型的iec61400-1标准第四版来计算叶片的第二疲劳单叶片覆冰质量矩。(a3)叶轮半径：：r＝0.5d，其中d为叶轮直径；(b3)85％叶轮半径位置：c85_p＝0.85r；(c3)叶片各截面距轮毂中心o的距离数组：rr＝r+l，其中r为叶片各截面距叶根a的距离数组，l为叶根a距轮毂中心o的距离；(d3)85％叶轮半径位置处弦长：c85＝(c85_p–rr(t))(c(t+1)–c(t))/(rr(t+1)-rr(t))+c(t)其中，rr(t)代表85％叶轮半径位置前的位置，rr(t+1)代表85％叶轮半径位置后的位置，c(t)代表85％叶轮半径位置前的位置对应的弦长，c(t+1)代表85％叶轮半径位置后的位置对应的弦长；(e3)叶片的第二覆冰线密度即覆冰单位长度质量m2如下所示：m2＝0.5×0.125×c85×r(f3)叶根a位置到叶尖b位置的覆冰单位长度质量数组为：miec2＝rrm2/rr(tip)其中，rr(tip)为叶尖b距轮毂中心o的距离；(g3)覆冰单位长度质量矩数组为：iiec2＝miec2×r(h3)第二疲劳单叶片覆冰质量矩可以表示为：iiec2＝∫iiec2dr因此，通过上述步骤，可以获得叶片的第二疲劳单叶片覆冰质量矩iiec2。在疲劳工况下，步骤s14的找出计算出的第一覆冰质量矩与第二覆冰质量矩相等时的风机叶片的最小叶尖结冰弦长包括：找出计算出的第一覆冰质量矩igl与第二疲劳单叶片覆冰质量矩iiec2相等时的风机叶片的最小叶尖结冰弦长tipchord。相应地，步骤s15中的建立风机叶片的覆冰模型包括：建立风机叶片的疲劳覆冰模型，在步骤s15中的建立风机叶片的覆冰模型还包括：在载荷仿真计算软件中将风机叶片的原线密度增加第二覆冰线密度的0.5倍。原线密度指叶片自身的线密度，即结冰前的线密度。步骤s16中的计算获得风机叶片的覆冰载荷包括：计算获得风机叶片的疲劳覆冰的机组载荷。具体地，将找出的第一覆冰质量矩igl与第二疲劳单叶片覆冰质量矩iiec2相等时的最小叶尖结冰弦长tipchord值输入到bladed软件中，勾选两叶片结冰，输入冰密度ρe＝700kg/m3。之后，修改bladed软件模型中的原线密度，在之前的基础上增加第二覆冰线密度的0.5倍。保存模型，完成疲劳载荷分析覆冰建模。依据建立的疲劳覆冰bladed模型，按照iec61400-1标准第四版的覆冰工况要求进行工况假定及bladed仿真计算，并对载荷计算结果进行后处理。从而可以最终获得风机叶片的疲劳覆冰的机组载荷。本发明实施例基于iec61400-1标准第四版中叶片覆冰质量计算方法及覆冰后叶片气动性能参数(升阻力系数)修正方法，搭建了合理且有效的覆冰工况bladed模型，填补了按照iec61400-1标准第四版计算分析叶片覆冰载荷对各部件影响的空白。本发明实施例的风机叶片覆冰的机组载荷的计算方法实现了基于iec61400-1标准第四版中叶片覆冰质量计算方法及覆冰翼型气动性能参数(升阻力系数)修正方法在bladed软件中的建模，避免了gl2010规范中覆冰质量偏保守导致的冰载极限、疲劳载荷偏大问题，解决了依据iec61400-1标准第四版无法进行冰载计算及分析的问题。以上对本发明实施例所提供的风机叶片覆冰的机组载荷的计算方法进行了详细的介绍。本文中应用了具体个例对本发明实施例的风机叶片覆冰的机组载荷的计算方法进行了阐述，以上实施例的说明只是用于帮助理解本发明的核心思想，并不用以限制本发明。应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明的精神和原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也均应落入本发明所附权利要求书的保护范围内。当前第1页12