本申请涉及风电技术领域，尤其涉及一种风力发电机的极限载荷的确定方法、装置及可读存储介质。
背景技术：
风是没有公害的能源之一，而且取之不尽，用之不竭。对于缺水、缺燃料和交通不便的沿海岛屿、草原牧区、山区和高原地带，因地制宜地利用风力发电，非常适合，大有可为。风力发电是指利用风力发电机把风的动能转为电能。
技术实现要素：
本申请提供一种改进的风力发电机的极限载荷的确定方法、装置及可读存储介质。
根据本申请实施例的一个方面，提供一种风力发电机的极限载荷的确定方法，包括：
获得所述风力发电机的基准弯矩；及
根据关于所述基准弯矩的极限弯矩修正函数，确定表征所述极限载荷的极限弯矩，所述极限弯矩修正函数至少是通过最小二乘法拟合得到的。
根据本申请实施例的另一个方面，提供一种风力发电机的极限载荷的确定装置，包括一个或多个处理器，用于实现上述的确定方法。
根据本申请实施例的另一个方面，提供一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，实现上述的确定方法。
本申请一些实施例根据风力发电机的基准弯矩及关于基准弯矩的极限弯矩修正函数，来确定表征风力发电机极限载荷的极限弯矩，可快速确定风力发电机的极限载荷；并且，极限弯矩修正函数至少是通过最小二乘法拟合得到的，利用该函数可提升确定的极限弯矩的精确度。
附图说明
图1所示为本申请一实施例中在叶根坐标系下风力发电机的极限弯矩的示意图；
图2所示为本申请一实施例中在轮毂旋转坐标系下风力发电机的极限弯矩的示意图；
图3所示为本申请一实施例中在轮毂固定坐标系下风力发电机的极限弯矩的示意图；
图4所示为本申请一实施例中在塔顶坐标系下风力发电机的极限弯矩的示意图；
图5所示为本申请一实施例中在塔底坐标系下风力发电机的极限弯矩的示意图；
图6所示为本申请一实施例中风力发电机的极限载荷的确定方法的流程图；
图7所示为图6所示确定方法的获得基准弯矩步骤的一实施例子流程图；
图8所示为本申请一实施例中叶片转动时叶素扫出的环面示意图；
图9所示为图6所示确定方法的确定极限弯矩步骤的子流程图；
图10所示为本申请一实施例中方位角的示意图；
图11所示为本申请一实施例中风力发电机的侧视图；
图12所示为通过仿真得到的绕x轴的叶根极限弯矩与通过本申请提供的方法得到的绕x轴的叶根极限弯矩的对比示意图；
图13所示为通过仿真得到的绕y轴的叶根极限弯矩与通过本申请提供的方法得到的绕y轴的叶根极限弯矩的对比示意图；
图14所示为通过仿真得到的绕z轴的叶根极限弯矩与通过本申请提供的方法得到的绕z轴的叶根极限弯矩的对比示意图；
图15所示为通过仿真得到的绕x轴的轮毂旋转极限弯矩与通过本申请提供的方法得到的绕x轴的轮毂旋转极限弯矩的对比示意图；
图16所示为通过仿真得到的绕y轴的轮毂旋转极限弯矩与通过本申请提供的方法得到的绕y轴的轮毂旋转极限弯矩的对比示意图；
图17所示为通过仿真得到的绕z轴的轮毂旋转极限弯矩与通过本申请提供的方法得到的绕z轴的轮毂旋转极限弯矩的对比示意图；
图18所示为通过仿真得到的绕x轴的轮毂固定极限弯矩与通过本申请提供的方法得到的绕x轴的轮毂固定极限弯矩的对比示意图；
图19所示为通过仿真得到的绕y轴的轮毂固定极限弯矩与通过本申请提供的方法得到的绕y轴的轮毂固定极限弯矩的对比示意图；
图20所示为通过仿真得到的绕z轴的轮毂固定极限弯矩与通过本申请提供的方法得到的绕z轴的轮毂固定极限弯矩的对比示意图；
图21所示为通过仿真得到的绕x轴的塔顶极限弯矩与通过本申请提供的方法得到的绕x轴的塔顶极限弯矩的对比示意图；
图22所示为通过仿真得到的绕y轴的塔顶极限弯矩与通过本申请提供的方法得到的绕y轴的塔顶极限弯矩的对比示意图；
图23所示为通过仿真得到的绕z轴的塔顶极限弯矩与通过本申请提供的方法得到的绕z轴的塔顶极限弯矩的对比示意图；
图24所示为通过仿真得到的绕x轴的塔底极限弯矩与通过本申请提供的方法得到的绕x轴的塔底极限弯矩的对比示意图；
图25所示为通过仿真得到的绕y轴的塔底极限弯矩与通过本申请提供的方法得到的绕y轴的塔底极限弯矩的对比示意图；
图26所示为通过仿真得到的绕z轴的塔底极限弯矩与通过本申请提供的方法得到的绕z轴的塔底极限弯矩的对比示意图；
图27所示为本申请一实施例中风力发电机的极限载荷的确定装置的框图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
风力发电机的载荷包括稳态载荷、周期载荷、瞬态载荷、随机载荷、共振诱导载荷、极限载荷、疲劳载荷等。其中极限载荷对于风力发电机的结构部件的强度设计至关重要，是保证风力发电机正常运转的关键因素之一。
本申请实施例的风力发电机的极限载荷的确定方法，包括：获得所述风力发电机的基准弯矩；及根据关于所述基准弯矩的极限弯矩修正函数，确定表征所述极限载荷的极限弯矩，所述极限弯矩修正函数至少是通过最小二乘法拟合得到的。
下面结合附图，对本申请提供的风力发电机的极限载荷的确定方法、装置及可读存储介质进行详细说明。在不冲突的情况下，下述的实施例及实施方式中的特征可以相互组合。
图1至图5是不同坐标系下风力发电机的极限弯矩的示意图。参见图1，风力发电机包括叶片101、机舱102、轮毂103及塔架104。塔架104从基础(未图示)向上延伸，机舱102安装在塔架104的顶端，轮毂103安装在机舱102的一端，叶片101安装在轮毂103上。图示实施例中，风力发电机包括三个叶片101，在其他实施例中，叶片101的数量可不为三个。风力发电机的极限弯矩可包括在叶根坐标系下绕x轴的叶根极限弯矩mxbr_1、在叶根坐标系下绕y轴的叶根极限弯矩mybr_1、以及在叶根坐标系下绕z轴的叶根极限弯矩mzbr_1。叶根1011为叶片101与轮毂103连接的根部。
参见图2，风力发电机的极限弯矩还可包括在轮毂旋转坐标系下绕x轴的轮毂旋转极限弯矩mxhr_1、在轮毂旋转坐标系下绕y轴的轮毂旋转极限弯矩myhr_1、以及在轮毂旋转坐标系下绕z轴的轮毂旋转极限弯矩mzhr_1。
参见图3，风力发电机的极限弯矩还可包括在轮毂固定坐标系下绕x轴的轮毂固定极限弯矩mxhs_1。在轮毂固定坐标系下绕y轴的轮毂固定极限弯矩myhs_1、在轮毂固定坐标系下绕z轴的轮毂固定极限弯矩mzhs_1。
参见图4，风力发电机的极限弯矩还可包括在塔顶坐标系下绕x轴的塔顶极限弯矩mxtt_1、在塔顶坐标系下绕y轴的塔顶极限弯矩mytt_1、以及在塔顶坐标系下绕z轴的塔顶极限弯矩mztt_1。
参见图5，风力发电机的极限弯矩还可包括在塔底坐标系下绕x轴的塔底极限弯矩mxtb_1、在塔底坐标系下绕y轴的塔底极限弯矩mytb_1、以及在塔底坐标系下绕z轴的塔底极限弯矩mztb_1。
图6所示为风力发电机的极限载荷的确定方法100的一个实施例的流程图。参见图6，所述确定方法100包括如下步骤110和步骤120。
在步骤110中，获得风力发电机的基准弯矩。
在一些实施例中，基准弯矩和极限弯矩为绕同一坐标轴的不同弯矩，基准弯矩可以为瞬态弯矩或稳态弯矩。稳态弯矩可表征稳态载荷，稳态载荷指的是均匀风速下的载荷。瞬态弯矩可表征瞬态载荷，瞬态载荷可以是湍流等作用下的载荷。极限弯矩是瞬态弯矩的极限值。在一些实施例中，一个或多个坐标轴下的极限弯矩可以作为其他一个或多个坐标轴下的基准弯矩。同一坐标轴指同一坐标系下的同一坐标轴，不同坐标轴指同一坐标系下的不同坐标轴和不同坐标系下的坐标轴。
在步骤120中，根据关于所述基准弯矩的极限弯矩修正函数，确定表征所述极限载荷的极限弯矩，所述极限弯矩修正函数至少是通过最小二乘法拟合得到的。
极限弯矩可以是上文所述的任一极限弯矩。在一些实施例中，至少两个不同坐标轴下的极限弯矩可以根据不同的基准弯矩和/或不同的极限弯矩修正函数来确定。将步骤110获得的基准弯矩代入极限弯矩修正函数，可以计算得到极限弯矩。通过最小二乘法拟合得到极限弯矩修正函数，可以通过最小二乘法线性拟合或通过最小二乘法非线性拟合。
根据风力发电机的基准弯矩及关于基准弯矩的极限弯矩修正函数，来确定表征风力发电机极限载荷的极限弯矩，可快速确定风力发电机的极限载荷；并且，极限弯矩修正函数至少是通过最小二乘法拟合得到的，利用该函数可提升确定的极限弯矩的精确度。下文将对如何确定上文所述的多个极限弯矩进行逐一详细说明。并且，相对于采用数值仿真的手段来确定极限载荷的方案来说，采用本申请提供的确定方法可快速确定风力发电机的极限载荷，节约计算资源。
在一个示例性实施例中，风力发电机的极限弯矩包括绕y轴的叶根极限弯矩mybr_1，如图1所示。下面将介绍确定绕y轴的叶根极限弯矩mybr_1的过程。
在一些实施例中，获得基准弯矩的步骤110包括：获得风力发电机在叶根坐标系下绕y轴的叶根基准弯矩mybr。确定极限弯矩的步骤120包括：根据对应的所述极限弯矩修正函数，确定绕y轴的叶根基准弯矩mybr对应的极限弯矩。其中，在绕y轴的叶根基准弯矩mybr对应的极限弯矩为绕y轴的叶根极限弯矩mybr_1。极限弯矩修正函数为关于绕y轴的叶根基准弯矩mybr的绕y轴的叶根极限弯矩mybr_1的极限弯矩修正函数。
在一些实施例中，基准弯矩与极限弯矩为绕同一坐标轴的弯矩，绕y轴的叶根基准弯矩mybr和绕y轴的叶根极限弯矩mybr_1为在叶根坐标系下绕y轴的弯矩。在一些实施例中，基准弯矩为瞬态弯矩，绕y轴的叶根基准弯矩mybr为挥舞瞬态弯矩。在一些实施例中，可以根据动量理论推算出瞬态弯矩。参见图7，获得基准弯矩的步骤110包括如下步骤111及步骤112。在一些实施例中，获得绕y轴的叶根基准弯矩mybr的步骤可以包括步骤111及步骤112。
在步骤111中，确定稳态弯矩。在一些实施例中，确定绕y轴的叶根稳态弯矩mybr_2。
在一些实施例中，可以根据风轮半径r、额定功率p、功率系数cp和推力系数ct，确定绕y轴的叶根稳态弯矩mybr_2。可以根据叶素动量理论推算出关于风轮半径r、额定功率p、功率系数cp和推力系数ct的稳态弯矩函数。动量理论用来描述作用在风轮旋转平面上的推力与来流速度之间的关系，说明风轮究竟能从动能中转换成多少机械能。
在一些实施例中，根据动量理论，作用在风轮旋转平面的推力的表达式如下表达式(1)：
式中，t表示作用在风轮旋转平面的推力；ct表示推力系数；ρ表示空气密度；r表示风力发电机的风轮半径；u为来流速度。其中，对于理想水平轴风力发电机而言，推力系数ct＝8/9。
参见图8，风力发电机运转时，叶片101的叶素1014转动时形成一个环面1012，风轮旋转平面1013可看作是由一系列环面1012构成的。叶素1014为叶片101在风轮任意半径r处的一个基本单元，是由风轮任意半径r处翼型剖面延伸一小段厚度dr而形成的。风力发电机的单个叶片的绕y轴的叶根稳态弯矩mybr_2的表达式可如下表达式(2)：
式中，b表示叶片个数。
表达式(2)中的积分式进行积分，且结合表达式(1)，可以得到如下表达式(3)：
结合表达式(1)，将功率系数cp和风轮半径r作为特征参数对表达式(3)进行简化，可得到如下表达式(4)：
式中，p表示风力发电机的额定功率。依据统计数据，推力系数ct的值可以在0.41～0.84之间，功率系数cp的值可以在0.1～0.593之间。在一个实施例中，可取ct＝0.574，cp＝0.46。根据表达式(4)可以计算得到稳态弯矩，可以计算得到绕y轴的叶根稳态弯矩mybr_2。
在步骤112中，根据稳态弯矩，确定瞬态弯矩，作为基准弯矩。
推力t的表达式(1)可用来确定稳态弯矩，其中u为均匀风速。稳态弯矩是均匀风速下的弯矩。实际中由于重力、惯性力、来流速度、横风、风剪切、偏航、湍流、变桨等因素的影响，会产生瞬态弯矩，瞬态弯矩的极限值通常大于稳态弯矩的值。可以根据绕y轴的叶根稳态弯矩mybr_2确定在叶根坐标系下绕y轴的叶根的挥舞瞬态弯矩，作为绕y轴的叶根基准弯矩mybr。
在一些实施例中，瞬态弯矩等于稳态弯矩和放大系数的乘积，如表达式(5)：
mybr＝mybr_2cstd-ybr(5)
式中，cstd-ybr表示放大系数。大量的仿真结果表明，挥舞瞬态弯矩的极限值的大小一般为对应的稳态弯矩的大小的两倍以上，因此cstd-ybr的值一般大于或等于2。在一个实施例中，cstd-ybr的值可以取2.6。绕y轴的叶根基准弯矩mybr(挥舞瞬态弯矩)通过表达式(5)计算得到。
如此确定基准弯矩，进一步根据确定的基准弯矩和极限弯矩修正函数，确定极限弯矩。在一些实施例中，极限弯矩修正函数包括关于误差修正函数和所述基准弯矩的函数。所述误差修正函数至少是通过最小二乘法拟合得到的。极限弯矩修正函数是关于误差修正函数和基准弯矩的函数。误差修正函数包括拟合系数，误差修正函数确定后，拟合系数为定值。可以通过最小二乘法拟合得到误差修正函数中的拟合系数，确定误差修正函数，利用误差修正函数建立极限弯矩修正函数。在一些实施例中，极限弯矩修正函数可以包括基准弯矩与误差修正函数相除的关系。误差修正函数的值可以表征修正系数，极限弯矩可以等于基准弯矩除以修正系数。在另一些实施例中，极限弯矩修正函数可以包括基准弯矩与误差修正函数相乘的关系。极限弯矩可以为基准弯矩与修正系数的乘积。如此利用误差修正函数对基准弯矩进行修正，得到极限弯矩。
在一些实施例中，误差修正函数至少是利用仿真数据通过最小二乘法拟合得到的。误差修正函数的构造可以实际风机设计中的大量仿真数据为基础。可以利用已有的仿真数据通过最小二乘法拟合得到拟合系数，从而确定误差修正函数。仿真数据的准确度较高，采用仿真数据可使得确定的极限弯矩修正函数比较精确，进而可提升确定出的极限弯矩的精确度。在其他一些实施例中，可以利用其他途径得到的已有数据，确定误差修正函数。
在一些实施例中，绕y轴的叶根极限弯矩mybr_1的极限弯矩修正函数是关于误差修正函数和绕y轴的叶根基准弯矩mybr的函数，该误差修正函数至少是通过最小二乘法拟合得到的。在一些实施例中，上文提到的其他一个或多个极限弯矩的极限弯矩修正函数也可以是关于对应的误差修正函数和对应的基准弯矩的函数，且误差修正函数也至少是通过最小二乘法拟合得到的，例如，绕x轴的叶根极限弯矩mxbr_1、绕x轴的轮毂旋转极限弯矩mxhr_1。两个或更多个不同坐标轴下的极限弯矩的误差修正函数可以不同。
在一些实施例中，所述误差修正函数至少是通过最小二乘法非线性拟合得到的，可以得到较准确的误差修正函数。在一些实施例中，误差修正函数可以是采用曲面函数形式，通过对曲面函数进行最小二乘法拟合得到的。
在一些实施例中，误差修正函数是关于风轮半径和额定功率的函数。对于一个或多个极限弯矩，例如绕y轴的叶根极限弯矩mybr_1，风轮半径和额定功率对误差修正函数的值的影响较大，根据风轮半径和额定功率确定误差修正函数的值，可以对基准弯矩进行较准确的修正，得到较准确的极限弯矩。
在一些实施例中，误差修正函数至少是利用所述风力发电机的风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合得到的。在一些实施例中，风轮半径参考数据和额定功率参考数据可以为仿真数据，可以使用已有的仿真数据。风轮半径参考数据可以为仿真中所用的风轮半径的数据，额定功率参考数据可以为仿真中使用的额定功率的数据。至少利用所述风力发电机的风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合得到误差修正函数中的拟合系数。在一些实施例中，通过构造自变量包括风轮半径的二次项和额定功率的二次项的曲面函数，对曲面函数进行最小二乘法非线性拟合，可以得到更准确的误差修正函数。
在一些实施例中，误差修正函数至少是利用误差修正函数值的参考数据、风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合的方法得到的。误差修正函数值的参考数据可以根据仿真数据计算得到。在一些实施例中，误差修正函数值的参考数据可以根据基准弯矩参考数据和极限弯矩参考数据获得。误差修正函数是利用基准弯矩参考数据和极限弯矩参考数据、风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合的方法得到的。参考数据可以是仿真数据，也可以是根据仿真数据计算得到的数据。例如基准弯矩参考数据可以是依据仿真数据，利用稳态弯矩表达式(4)及表达式(5)计算得到的数据，风轮半径参考数据、额定功率参考数据和极限弯矩参考数据可以为仿真数据。
在一些实施例中，可以构造误差修正函数的表达式，如下表达式(6)所示：
fybr＝p00+p10x+p01y+p20x2+p11xy+p02y2(6)
式中，fybr表示误差修正函数的值，p00、p10、p01、p20、p11、p02表示拟合系数，x表示风力发电机的风轮半径，y表示风力发电机的额定功率。表达式(6)为误差修正函数关于风轮半径和额定功率的表达式。绕y轴的叶根极限弯矩mybr_1对应的误差修正函数的表达式可以为表达式(6)。
误差修正函数关于基准弯矩mybr和极限弯矩mybr_1的表达式为如下表达式(7)：
绕y轴的叶根极限弯矩mybr_1对应的误差修正函数的值可以为绕y轴的叶根基准弯矩mybr和绕y轴的叶根极限弯矩mybr_1的比值。
利用表达式(6)和(7)，根据拟合算法，得到如下表达式(8)所示的拟合误差函数：
式中，fybri行表示第i组数据中误差修正函数的值，通过表达式(7)计算得到；xi表示第i组数据中风轮半径的值，yi表示第i组数据中额定功率的值。
对函数s1的拟合系数p00、p10、p01、p20、p11、p02分别求偏导，并令其等于0，得到如下方程(9)：
在一些实施例中，获取n(n为大于1的整数)组参考数据。每组参考数据包括基准弯矩参考数据和极限弯矩参考数据，代入表达式(7)，计算得到误差修正函数值的参考数据，作为表达式(9)中的fybri的值。
每组参考数据还包括风轮半径参考数据和额定功率参考数据。将误差修正函数值的参考数据、风轮半径参考数据和额定功率参考数据代入方程(9)。通过求解方程(9)，可得到拟合系数p00、p10、p01、p20、p11、p02的值。将拟合系数p00、p10、p01、p20、p11、p02的值代入表达式(6)，即可得到误差修正函数关于风轮半径及额定功率的表达式。
绕y轴的叶根极限弯矩mybr_1对应的误差修正函数可以通过上述过程确定。上述采用最小二乘法拟合进行误差修正。利用已有的基准弯矩和极限弯矩的仿真数据，通过最小二乘法拟合进行误差修正。误差修正函数是关于风轮半径和功率的非线性函数，通过最小二乘法非线性拟合。以风轮半径和功率为变量，进行二次曲面拟合。
在一些实施例中，参见图9，确定极限弯矩的步骤120包括如下步骤121及步骤122。
在步骤121中，根据所述风力发电机的风轮半径和额定功率，确定所述误差修正函数的值。
将实际的风力发电机的风轮半径和额定功率代入误差修正函数，如表达式(6)，计算得到误差修正函数的值。
在步骤122中，根据确定的所述误差修正函数的值和所述基准弯矩，确定所述极限弯矩。
在一些实施例中，根据步骤121确定的所述误差修正函数的值和绕y轴的叶根基准弯矩mybr，确定绕y轴的叶根极限弯矩mybr_1。具体地，将误差修正函数的值及绕y轴的叶根基准弯矩mybr的值代入对应的极限弯矩修正函数，可得到绕y轴的叶根极限弯矩mybr_1的值。
在一些实施例中，绕y轴的叶根极限弯矩mybr_1的极限弯矩修正函数如表达式(10)所示：
在一个示例性实施例中，风力发电机的极限弯矩包括绕x轴的叶根极限弯矩mxbr_1。下面将详细介绍确定绕x轴的叶根极限弯矩mxbr_1的过程。
在一些实施例中，获得基准弯矩的步骤110包括：获得风力发电机在叶根坐标系下绕x轴的叶根基准弯矩mxbr。确定极限弯矩的步骤120包括：根据对应的极限弯矩修正函数，确定绕x轴的叶根基准弯矩mxbr对应的极限弯矩。其中，绕x轴的叶根基准弯矩mxbr对应的极限弯矩为绕x轴的叶根极限弯矩mxbr_1。绕x轴的叶根极限弯矩mxbr_1的极限弯矩修正函数为关于绕x轴的叶根基准弯矩mxbr的极限弯矩修正函数。绕x轴的叶根基准弯矩mxbr和绕x轴的叶根极限弯矩mxbr_1为绕同一坐标轴的不同弯矩。
在一些实施例中，绕x轴的叶根基准弯矩mxbr可以为瞬态弯矩，为在叶根坐标系下绕x轴的瞬态弯矩。再次参见图7，获得基准弯矩的步骤110包括如下步骤111及步骤112。
在步骤111中，确定稳态弯矩。确定绕x轴的叶根稳态弯矩mxbr_2。
分别确定平均转矩mxbr_21及重力引起的重力弯矩mxbr_22。
在一些实施例中，单个叶片的平均转矩mxbr_21可采用如下表达式(11)计算：
式中，p表示风力发电机的额定功率，ω表示风轮旋转速度，r表示风力发电机的风轮半径，λ表示叶尖速比，u表示来流风速，ρ表示空气密度，cp表示功率系数(此估算中忽略了传递过程中电损耗)。在一些实施例中，叶尖速比λ的值可取为10，功率系数cp的值可取为0.46。
在一些实施例中，假设叶片的中心位置与叶根的距离为r/3，单个叶片的重力弯矩mxbr_22可采用如下表达式(12)计算：
式中，massbl表示单个叶片的质量，r表示风力发电机的风轮半径。单个叶片的质量massbl可通过如下表达式(13)计算得到：
massbl＝cd1r3+cd2r2+cd3r+cd4+cd5p(13)
式中，cd1、cd2、cd3、cd4及cd5表示系数。在一些实施例中，系数cd1、cd2、cd3、cd4及cd5可利用参考数据拟合得到，参考数据包括风轮半径r参考数据、功率p参考数据和单个叶片的质量参考数据。
之后，根据单个叶片的平均转矩mxbr_21及重力弯矩mxbr_22，确定单个叶片的绕x轴的叶根稳态弯矩mxbr_2。单个叶片的绕x轴的叶根稳态弯矩mxbr_2为平均转矩mxbr_21及重力引起的重力弯矩mxbr_22之和，表达式如下表达式(14)：
mxbr_2＝mxbr_21+mxbr_22(14)
在本实施例中，根据绕x轴的叶根稳态弯矩mxbr_2确定摆阵瞬态弯矩，作为绕x轴的叶根基准弯矩mxbr。绕x轴的叶根稳态弯矩mxbr_2是均匀风速下的弯矩，实际中摆阵瞬态弯矩会受到湍流、变桨、偏航等因素的影响，摆阵瞬态弯矩的极限值一般比绕x轴的叶根稳态弯矩的值大。在一些实施例中，摆阵瞬态弯矩(可以为绕x轴的叶根基准弯矩mxbr)等于绕x轴的叶根稳态弯矩mxbr_2与放大系数的乘积，如表达式(15)：
mxbr＝mxbr_2cstd-xbr(15)
式中，cstd-xbr表示放大系数。大量的仿真结果表明，摆阵瞬态弯矩的大小一般为绕x轴的叶根稳态弯矩mxbr_2的大小的两倍以上，因此cstd-xbr的值一般大于或等于2。在一些实施例中，cstd-xbr的值可以取2.6。
在一些实施例中，绕x轴的叶根极限弯矩mxbr_1的极限弯矩修正函数包括关于误差修正函数和基准弯矩的函数。可以构造绕x轴的叶根极限弯矩mxbr_1的误差修正函数的表达式，可如下表达式(16)所示：
fxbr＝p21+p22x+p23y+p24x2+p25xy+p26y2(16)
式中，fxbr表示误差修正函数的值，p21、p22、p23、p24、p25、p26表示拟合系数，x表示风力发电机的风轮半径，y表示风力发电机的额定功率。
绕x轴的叶根极限弯矩mxbr_1对应的误差修正函数至少是通过最小二乘法拟合得到的。通过最小二乘法拟合得到拟合系数，从而得到误差修正函数。在一些实施例中，拟合系数p21、p22、p23、p24、p25、p26可利用参考数据拟合得到，参考数据包括基准弯矩参考数据、绕x轴的叶根极限弯矩mxbr_1的参考数据、风轮半径参考数据和额定功率参考数据。在一些实施例中，基准弯矩参考数据为摆阵瞬态弯矩参考数据。参考数据可以是已有的仿真数据，或者是根据仿真数据计算得到的数据。绕x轴的叶根极限弯矩mxbr_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
拟合系数p21、p22、p23、p24、p25、p26的求解过程与表达式(6)的误差修正函数fybr中拟合系数p00、p10、p01、p20、p11、p02的求解过程类似。绕x轴的叶根极限弯矩mxbr_1对应的误差修正函数至少是通过最小二乘法非线性拟合得到的。绕x轴的叶根极限弯矩mxbr_1对应的误差修正函数至少是利用所述风力发电机的风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合的方法得到的。具体求解过程参见上文所述，不再进行赘述。将p21、p22、p23、p24、p25、p26的值代入表达式(16)，得到误差修正函数关于风轮半径及额定功率的表达式。
在一些实施例中，确定极限弯矩的步骤120包括：根据风力发电机的风轮半径和额定功率，确定误差修正函数的值；且根据确定的误差修正函数的值和基准弯矩，确定极限弯矩。根据风轮半径和额定功率，确定绕x轴的叶根极限弯矩mxbr_1对应的误差修正函数的值。将实际的风力发电机的风轮半径和额定功率代入误差修正函数的表达式(16)，得到误差修正函数的值。
根据确定的误差修正函数的值和绕x轴的叶根基准弯矩mxbr，确定绕x轴的叶根极限弯矩mxbr_1。具体地，将误差修正函数的值及绕x轴的叶根基准弯矩mxbr的值代入绕x轴的叶根极限弯矩mxbr_1的极限弯矩修正函数，得到绕x轴的叶根极限弯矩mxbr_1的值。在一些实施例中，关于绕x轴的叶根极限弯矩mxbr_1的极限弯矩修正函数如表达式(17)所示：
将误差修正函数的值和绕x轴的叶根基准弯矩mxbr的值代入表达式(17)，得到绕x轴的叶根极限弯矩mxbr_1的值。
在一个示例性实施例中，风力发电机的极限弯矩包括绕z轴的叶根极限弯矩mzbr_1。下面将介绍确定绕z轴的叶根极限弯矩mzbr_1的过程。
在一些实施例中，绕z轴的叶根极限弯矩mzbr_1对应的基准弯矩包括绕y轴的叶根极限弯矩mybr_1。确定绕z轴的叶根极限弯矩mzbr_1的步骤包括：根据关于绕y轴的叶根极限弯矩mybr_1的绕z轴的叶根极限弯矩mzbr_1的极限弯矩修正函数，确定绕z轴的叶根极限弯矩mzbr_1。
在一些实施例中，绕z轴的叶根极限弯矩mzbr_1的极限弯矩修正函数至少是利用最小二乘法非线性拟合得到的。在一些实施例中，绕z轴的叶根极限弯矩mzbr_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。在一些实施例中，极限弯矩修正函数至少是利用样条最小二乘法拟合得到的。极限弯矩修正函数例如可采用三次分段样条函数进行拟合得到，对于每一段区间，均满足极限弯矩修正函数。绕z轴的叶根极限弯矩mzbr_1的极限弯矩修正函数至少是利用样条最小二乘法拟合得到的。
在一些实施例中，关于绕y轴的叶根极限弯矩mybr_1的绕z轴的叶根极限弯矩mzbr_1的极限弯矩修正函数如表达式(18)所示：
式中，c1、c2、c3及c4分别表示系数，分段系数为系数c1、c2、c3及c4组成的系数矩阵。
系数c1、c2、c3及c4可以利用参考数据通过样条最小二乘法拟合得到。参考数据包括绕y轴的叶根极限弯矩参考数据和绕z轴的叶根极限弯矩参考数据。参考数据可以是仿真数据。可以采用三次分段样条函数进行拟合得到系数c1、c2、c3及c4。系数c1、c2、c3及c4确定后，代入表达式(18)，得到绕z轴的叶根极限弯矩mzbr_1的极限弯矩修正函数。
表达式(18)的极限弯矩修正函数确定后，将实际的绕y轴的叶根极限弯矩mybr_1代入关系式(18)，得到绕z轴的叶根极限弯矩mzbr_1。绕y轴的叶根极限弯矩mybr_1可以通过上文所述的方法得到，之后可以计算得到绕z轴的叶根极限弯矩mzbr_1。
在另一些实施例中，绕x轴的叶根稳态弯矩可以作为基准弯矩，可以通过误差修正函数修正湍流等动态因素的影响。
在一个示例性实施例中，风力发电机的极限弯矩包括在轮毂旋转坐标系下绕x轴的轮毂旋转极限弯矩mxhr_1。下面将介绍确定绕x轴的轮毂旋转极限弯矩mxhr_1的过程。
在一些实施例中，获得基准弯矩的步骤110，包括：获得风力发电机在轮毂旋转坐标系下绕x轴的轮毂旋转基准弯矩mxhr；确定极限弯矩的步骤120，包括：根据对应的极限弯矩修正函数，确定绕x轴的轮毂旋转基准弯矩mxhr对应的极限弯矩。其中，绕x轴的轮毂旋转基准弯矩mxhr对应的极限弯矩为绕x轴的轮毂旋转极限弯矩mxhr_1。极限弯矩为绕x轴的轮毂旋转极限弯矩mxhr_1与绕x轴的轮毂旋转基准弯矩mxhr为绕同一坐标轴的不同弯矩。
在一些实施例中，绕x轴的轮毂旋转基准弯矩mxhr为瞬态弯矩，为在轮毂旋转坐标系下绕x轴的轮毂旋转瞬态弯矩。在一些实施例中，再次参见图7，获得基准弯矩的步骤110包括如下步骤111及步骤112。
在步骤111中，确定稳态弯矩。
确定在轮毂旋转坐标系下绕x轴的轮毂旋转稳态弯矩mxhr_2。在一些实施例中，通过如下表达式(19)和表达式(20)来确定旋转x轴稳态弯矩mxhr_2：
式中，ω表示风轮旋转速度，p表示风力发电机的额定功率，λ表示叶尖速比，ρ表示空气密度，r表示风力发电机的风轮半径，cp表示功率系数(此估算中忽略了传递过程中电损耗)。
根据旋转x轴稳态弯矩mxhr_2，确定绕x轴的轮毂旋转瞬态弯矩，作为绕x轴的轮毂旋转基准弯矩mxhr。绕x轴的轮毂旋转稳态弯矩mxhr_2表征的是稳态载荷，实际中由于湍流、变桨、偏航等因素的影响，会使得瞬态弯矩的值大于稳态弯矩的值。在一些实施例中，可采用如下表达式(21)确定绕x轴的轮毂旋转瞬态弯矩：
mxhr＝mxhr_2cstd-xhr(21)
式中，cstd-xhr表示放大系数。在一些实施例中，根据仿真结果，cstd-xhr的值可以取1.8。
在一些实施例中，绕x轴的轮毂旋转极限弯矩mxhr_1的极限弯矩修正函数包括关于对应的误差修正函数和对应的基准弯矩的函数。在一些实施例中，可以构造误差修正函数的表达式，如下表达式(22)所示：
fxhr＝p31+p32x+p33y+p34x2+p35xy+p36y2(22)
式中，fxhr表示误差修正函数的值，p31、p32、p33、p34、p35、p36表示拟合系数，x表示风力发电机的风轮半径，y表示额定功率。
在一些实施例中，绕x轴的轮毂旋转极限弯矩mxhr_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。在一些实施例中，绕x轴的轮毂旋转极限弯矩mxhr_1的极限弯矩修正函数包括的误差修正函数至少是通过最小二乘法拟合得到的，进而得到极限弯矩修正函数。误差修正函数的拟合系数p31、p32、p33、p34、p35、p36可通过最小二乘法拟合得到，进而得到误差修正函数。在一些实施例中，误差修正函数至少是通过最小二乘法非线性拟合得到的。误差修正函数的拟合系数p31、p32、p33、p34、p35、p36可通过最小二乘法非线性拟合得到。在一个实施例中，误差修正函数可以是采用最小二乘法的二次曲面拟合得到的。误差修正函数的拟合系数p31、p32、p33、p34、p35、p36可采用最小二乘法的二次曲面拟合得到。
在一些实施例中，拟合系数p31、p32、p33、p34、p35、p36可利用参考数据拟合得到。参考数据可以是已有的仿真数据或根据已有的仿真数据计算得到的数据。该误差修正函数可以至少是利用仿真数据通过最小二乘法拟合得到的。参考数据可以包括绕x轴的轮毂旋转基准弯矩mxhr的参考数据、绕x轴的轮毂旋转极限弯矩mxhr_1的参考数据、风轮半径参考数据和额定功率参考数据。其中，绕x轴的轮毂旋转基准弯矩mxhr的参考数据可以是依据仿真数据，利用表达式(19)、(20)及表达式(21)计算得到，绕x轴的轮毂旋转极限弯矩mxhr_1的参考数据、风轮半径参考数据和额定功率参考数据可以是仿真数据。误差修正函数至少是利用风力发电机的风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合的方法得到的。利用绕x轴的轮毂旋转基准弯矩mxhr的参考数据、绕x轴的轮毂旋转极限弯矩mxhr_1的参考数据、风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合得到拟合系数p31、p32、p33、p34、p35、p36。拟合系数p31、p32、p33、p34、p35、p36的求解过程与上文所述的误差修正函数fybr中拟合系数p00、p10、p01、p20、p11、p02的求解过程类似，不再进行赘述。
拟合系数p31、p32、p33、p34、p35、p36确定后，将拟合系数p31、p32、p33、p34、p35、p36的值代入表达式(22)，得到关于风轮半径及额定功率的误差修正函数。
在一个实施例中，再次参见图9，确定绕x轴的轮毂旋转极限弯矩mxhr_1的步骤包括如下步骤121及步骤122。
在步骤121中，根据风轮半径和额定功率，确定误差修正函数的值。
将风轮半径和额定功率代入误差修正函数的表达式(21)中，可确定误差修正函数的值。
在步骤122中，根据确定的误差修正函数的值和基准弯矩，确定极限弯矩。
根据确定的误差修正函数的值和绕x轴的轮毂旋转基准弯矩mxhr，确定绕x轴的轮毂旋转极限弯矩mxhr_1。具体地，将误差修正函数的值及绕x轴的轮毂旋转基准弯矩mxhr的值，代入对应的极限弯矩修正函数，可得到绕x轴的轮毂旋转极限弯矩mxhr_1的值。
在一些实施例中，关于绕x轴的轮毂旋转基准弯矩mxhr的绕x轴的轮毂旋转极限弯矩mxhr_1的极限弯矩修正函数如表达式(23)所示：
在另一些实施例中，旋转x轴稳态弯矩可以作为基准弯矩，可以通过误差修正函数修正湍流等动态因素的影响。
在一个示例性实施例中，风力发电机的极限弯矩包括绕y轴的轮毂旋转极限弯矩myhr_1。下面将介绍确定绕y轴的轮毂旋转极限弯矩myhr_1的过程。
在一些实施例中，基准弯矩包括绕y轴的叶根极限弯矩mybr_1；确定极限弯矩的步骤120包括：根据关于绕y轴的叶根极限弯矩mybr_1的极限弯矩修正函数，确定绕y轴的轮毂旋转极限弯矩myhr_1。
在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。绕y轴的轮毂旋转极限弯矩myhr_1的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕y轴的轮毂旋转极限弯矩myhr_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，关于绕y轴的叶根极限弯矩mybr_1的绕y轴的轮毂旋转极限弯矩myhr_1的力学关系如下表达式(24)：
myhr_1＝mybr_1+fxrl(24)
式中，fx表示叶根沿着x轴方向的力，rl表示叶根长度。
依据表达式(24)，极限弯矩myhr_1的极限弯矩修正函数可表示成如下表达式(25)：
myhr_1＝k1mybr_1+b1(25)
式中，k1及b1表示拟合系数。
绕y轴的轮毂旋转极限弯矩myhr_1的极限弯矩修正函数是关于绕y轴的叶根极限弯矩mybr_1的线性函数。拟合系数k1及b1可以通过最小二乘法线性拟合得到。依据拟合算法，得到如下表达式(26)所示的拟合误差函数：
式中，myhr_1i表示第i组风机载荷数据中轮毂的绕y轴的轮毂旋转极限弯矩，mybr_1i表示第i组数据中叶根的绕y轴的叶根极限弯矩，k1及b1表示拟合系数。
对拟合误差函数s2的拟合系数k1、b1分别求偏导，并令其等于0，得到如下方程(27)：
在一些实施例中，获取n(n为大于1的整数)组参考数据。每组参考数据包括绕y轴的叶根极限弯矩mybr_1i的参考数据和绕y轴的轮毂旋转极限弯矩myhr_1i的参考数据，代入表达式(27)。参考数据可以为已有的仿真数据。利用仿真数据通过最小二乘法拟合得到拟合系数k1及b1。求解方程(27)，可得到拟合系数k1及b1。将k1及b1的值代入表达式(25)，得到绕y轴的轮毂旋转极限弯矩mybr_1的极限弯矩修正函数。
将实际的绕y轴的叶根极限弯矩mybr_1的值代入表达式(25)，得到绕y轴的轮毂旋转极限弯矩myhr_1。绕y轴的叶根极限弯矩mybr_1可以通过上文所述的方法得到，之后可以计算得到绕y轴的轮毂旋转极限弯矩myhr_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕z轴的轮毂旋转极限弯矩mzhr_1。下面将介绍确定绕z轴的轮毂旋转极限弯矩mzhr_1的过程。
在一些实施例中，基准弯矩包括绕x轴的叶根极限弯矩mxbr_1及在叶根坐标系下绕z轴的叶根极限弯矩mzbr_1；确定极限弯矩的步骤120包括：根据关于绕x轴的叶根极限弯矩mxbr_1和绕z轴的叶根极限弯矩mzbr_1的极限弯矩修正函数，确定绕z轴的轮毂旋转极限弯矩mzhr_1。
在一些实施例中，绕z轴的轮毂旋转极限弯矩的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，绕z轴的轮毂旋转极限弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。在一些实施例中，关于绕x轴的叶根极限弯矩mxbr_1和绕z轴的叶根极限弯矩mzbr_1的绕z轴的轮毂旋转极限弯矩的力学关系如下表达式(28)：
mzhr_1＝mzbr_1cosθ+mxbr_1sinθ+fxrlsinθ(28)
式中，fx表示叶根坐标系下沿着x轴方向的力，rl表示叶根长度，θ为锥角(也即是风力发电机的轮毂的中心轴与水平方向的夹角)。
依据表达式(28)，极限弯矩mzhr_1的极限弯矩修正函数可表示成如下表达式(29)：
mzhr_1＝k2(mzbr_1cosθ+mxbr_1sinθ)+b2(29)
式中，k2及b2表示拟合系数。
在一些实施例中，θ可为5°。拟合系数k2及b2可利用参考数据来确定，参考数据包括绕x轴的叶根极限弯矩参考数据、绕z轴的叶根极限弯矩参考数据及绕z轴的轮毂旋转极限弯矩参考数据。参考数据可以是仿真数据。拟合系数k2及b2的求解过程与拟合系数k1、b1的求解过程类似，不再进行赘述。
将实际的绕x轴的叶根极限弯矩mxbr_1和绕z轴的叶根极限弯矩mzbr_1的值代入表达式(29)，得到绕z轴的轮毂旋转极限弯矩mzhr_1。绕x轴的叶根极限弯矩mxbr_1和绕z轴的叶根极限弯矩mzbr_1可以通过上文所述的方法得到，之后可以计算得到绕z轴的轮毂旋转极限弯矩mzhr_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕x轴的轮毂固定极限弯矩mxhs_1。下面将介绍确定绕x轴的轮毂固定极限弯矩mxhs_1的过程。
在一些实施例中，在轮毂固定坐标系下的x轴与在轮毂旋转坐标系下的x轴重合，绕x轴的轮毂固定基准弯矩mxhs与绕x轴的轮毂旋转基准弯矩mxhr相等，且绕x轴的轮毂固定极限弯矩mxhs_1与绕x轴的轮毂旋转极限弯矩mxhs_1相等。
在一些实施例中，基准弯矩包括绕x轴的轮毂旋转极限弯矩mxhr_1，确定极限弯矩的步骤120，包括：根据极限弯矩修正函数，确定绕x轴的轮毂固定极限弯矩mxhs_1。
关于所述绕x轴的轮毂固定极限弯矩mxhs_1的极限弯矩修正函数可如下表达式(29)所示：
mxhs_1＝fxhs×mxhr_1(29)
在一个实施例中，fxhs＝1，也即是绕x轴的轮毂固定极限弯矩mxhs_1与绕x轴的轮毂旋转极限弯矩mxhr_1相等。
在一个示例性实施例中，风力发电机的极限弯矩包括绕y轴的轮毂固定极限弯矩myhs_1。下面将介绍确定绕y轴的轮毂固定极限弯矩myhs_1的过程。
在一些实施例中，所述基准弯矩包括绕y轴的轮毂旋转极限弯矩myhr_1和绕z轴的轮毂旋转极限弯矩mzhr_1；所述确定极限弯矩的步骤120，包括：根据极限弯矩修正函数，确定绕y轴的轮毂固定极限弯矩myhs_1。
在一些实施例中，所述极限弯矩修正函数至少是利用最小二乘法线性拟合得到的。绕y轴的轮毂固定极限弯矩的极限弯矩修正函数至少是利用最小二乘法线性拟合得到的。在一些实施例中，所述极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕y轴的轮毂固定极限弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。在一些实施例中，关于绕y轴的轮毂旋转极限弯矩myhr_1和绕z轴的轮毂旋转极限弯矩mzhr_1的绕y轴的轮毂固定极限弯矩的力学关系如下表达式(31)：
myhs_1＝myhr_1cosβ-mzhr_1sinβ(31)
对表达式(31)取极值，可得到如下表达式(32)：
式中，ψ为系数，依据表达式(32)，极限弯矩myhs_1的极限弯矩修正函数可表示成如下表达式(33)：
式中，k3及b3表示拟合系数。
在一些实施例中，拟合系数k3及b3可利用参考数据来确定，参考数据包括绕y轴的轮毂旋转极限弯矩参考数据、绕z轴的轮毂旋转极限弯矩参考数据和绕y轴的轮毂固定极限弯矩参考数据。参考数据可以是仿真数据。拟合系数k3及b3的求解过程与上文所述的拟合误差函数s2中拟合系数k1、b1的求解过程类似，不再进行赘述。通过最小二乘法线性拟合得到拟合系数k3及b3。
将实际的绕y轴的轮毂旋转极限弯矩myhr_1的值和绕z轴的轮毂旋转极限弯矩mzhr_1的值代入表达式(33)，得到绕y轴的轮毂固定极限弯矩myhs_1。绕y轴的轮毂旋转极限弯矩myhr_1和绕z轴的轮毂旋转极限弯矩mzhr_1可以通过上文所述的方法得到，之后可以计算得到绕y轴的轮毂固定极限弯矩myhs_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕z轴的轮毂固定极限弯矩mzhs_1。下面将介绍确定绕z轴的轮毂固定极限弯矩mzhs_1的过程。
在一些实施例中，所述基准弯矩包括绕y轴的轮毂旋转极限弯矩myhr_1和绕z轴的轮毂旋转极限弯矩mzhr_1；确定极限弯矩的步骤120，包括：根据极限弯矩修正函数，确定绕z轴的轮毂固定极限弯矩mzhs_1。
在一些实施例中，所述极限弯矩修正函数至少是利用最小二乘法线性拟合得到的。绕z轴的轮毂固定极限弯矩的极限弯矩修正函数至少是利用最小二乘法线性拟合得到的。在一些实施例中，所述极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕z轴的轮毂固定极限弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。在一些实施例中，关于绕y轴的轮毂旋转极限弯矩myhr_1和绕z轴的轮毂旋转极限弯矩mzhr_1的绕z轴的轮毂固定极限弯矩的力学关系如下表达式：
mzhs_1＝mzhr_1cosβ-myhr_1sinβ(34)
对表达式(33)取极值，可得到如下表达式(35)：
式中，为系数，依据表达式(35)，极限弯矩mzhs_1的极限弯矩修正函数可表示成如下表达式(36)：
式中，k4及b4表示拟合系数。
在一些实施例中，拟合系数k4及b4可利用参考数据来确定，参考数据包括绕y轴的轮毂旋转极限弯矩参考数据、绕z轴的轮毂旋转极限弯矩参考数据和绕z轴的轮毂固定极限弯矩参考数据。参考数据可以是仿真数据。拟合系数k4及b4的求解过程与拟合误差函数s2中拟合系数k1、b1的求解过程类似，不再进行赘述。通过最小二乘法线性拟合得到拟合系数k4及b4。
将实际的绕y轴的轮毂旋转极限弯矩myhr_1的值和绕z轴的轮毂旋转极限弯矩mzhr_1的值代入表达式(36)，得到绕z轴的轮毂固定极限弯矩mzhs_1。绕y轴的轮毂旋转极限弯矩myhr_1和绕z轴的轮毂旋转极限弯矩mzhr_1可以通过上文所述的方法得到，之后可以计算得到绕z轴的轮毂固定极限弯矩mzhs_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕y轴的塔顶极限弯矩mytt_1。下面将介绍确定绕y轴的塔顶极限弯矩mytt_1的过程。
在一些实施例中，基准弯矩包括作用在风轮旋转平面上的推力产生的绕y轴的第一塔顶稳态弯矩mytt_thrust、风轮重力产生的绕y轴的第二塔顶稳态弯矩mytt_g、以及机舱和风轮的惯性力产生的绕y轴的第三塔顶稳态弯矩mytt_a；确定极限弯矩的步骤120，包括：根据关于所述第一塔顶稳态弯矩mytt_thrust、所述第二塔顶稳态弯矩mytt_g和所述第三塔顶稳态弯矩mytt_a的所述极限弯矩修正函数，确定绕y轴的塔顶极限弯矩mytt_1。
在一些实施例中，所述极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。绕y轴的塔顶极限弯矩mytt_1的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，所述极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕y轴的塔顶极限弯矩mytt_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，第一塔顶稳态弯矩mytt_thrust可采用如下表达式(37)得到：
参见图11，表达式(37)中的hhfs表示机舱中心线距离塔顶的垂向距离，ct表示推力系数，ρ表示空气密度，r表示风力发电机的风轮半径，u表示来流速度。
在一些实施例中，第二塔稳态限弯矩mytt_g可采用如下表达式(38)得到：
mytt_g＝-massrotorglohg(38)
再次参见图11，表达式(38)中的lohg表示风轮与塔架中心线之间的纵向距离(也即是在塔顶坐标系下的x轴方向上的距离)，massrotor表示风轮质量，g为重力加速度。在一些实施例中，未考虑机舱在塔架中心线前面部分的质量。
在一些实施例中，第三塔稳态限弯矩mytt_a可采用如下表达式(39)得到：
mytt_a＝-(massrotor+massnacelle)ahhfs(39)
式中，massnacelle表示机舱质量，a表示机舱加速度，a的值可取为2m/s2，massrotor表示风轮质量。
在一些实施例中，可构造绕y轴的塔顶极限弯矩mytt_1的极限弯矩修正函数如下表达(40)及(41)：
mytt_temp＝a00+a10mytt_thrust+a01mytt_g(40)
mytt_1＝b00+b10mytt_a+b01mytt_temp(41)
式中，mytt_temp表示极限弯矩的中间变量，设置中间变量可以方便计算；a00、a10、a01、b00、b10、b01表示拟合系数。
在一些实施例中，拟合系数a00、a10、a01、b00、b10、b01可利用参考数据得到，参考数据根据仿真数据，利用稳态弯矩表达式(37)、(38)和(39)计算得到的第一塔顶稳态弯矩参考数据、第二塔顶稳态弯矩参考数据及第三塔顶稳态弯矩参考数据。绕y轴的塔顶极限弯矩参考数据可以是仿真数据。拟合系数a00、a10、a01、b00、b10、b01的求解过程与上文所述的误差修正函数fybr中拟合系数p00、p10、p01、p20、p11、p02的求解过程类似，不再进行赘述。
求解出a00、a10、a01、b00、b10、b01后，将第一塔稳态限弯矩mytt_thrust、第二塔顶稳态弯矩mytt_g和第三塔顶稳态弯矩mytt_a的值代入表达式(40)和表达式(41)，得到绕y轴的塔顶极限弯矩mytt_1。
在另一些实施例中，可以直接构造绕y轴的塔顶极限弯矩mytt_1关于第一塔顶稳态弯矩mytt_thrust、第二塔顶稳态弯矩mytt_g、第三塔顶稳态弯矩mytt_a的极限弯矩修正函数。
在一个示例性实施例中，风力发电机的极限弯矩包括绕z轴的塔顶极限弯矩mztt_1。下面将介绍确定绕z轴的塔顶极限弯矩mztt_1的过程。
在一些实施例中，所述基准弯矩包括绕z轴的轮毂固定极限弯矩mzhs_1；所述确定极限弯矩的步骤120，包括：根据关于所述绕z轴的轮毂固定极限弯矩mzhs_1的所述极限弯矩修正函数，确定绕z轴的塔顶极限弯矩mztt_1。
在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。绕z轴的塔顶极限弯矩mztt_1的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕z轴的塔顶极限弯矩mztt_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，关于绕z轴的轮毂固定极限弯矩mzhs_1的绕z轴的塔顶极限弯矩mztt_1的极限弯矩修正函数如下表达式(42)所示：
mztt_1＝k5mzhs_1+b5(42)
式中，k5及b5表示拟合系数。
在一些实施例中，拟合系数k5及b5可利用参考数据通过最小二乘法拟合得到。拟合系数k5及b5可以通过最小二乘法线性拟合得到。参考数据包括绕z轴的轮毂固定极限弯矩参考数据和绕z轴的塔顶极限弯矩参考数据。参考数据可以是仿真数据。拟合系数k5及b5的求解过程与拟合误差函数s2中拟合系数k1、b1的求解过程类似，不再进行赘述。
将实际的绕z轴的轮毂固定极限弯矩mzhs_1的值代入表达式(42)，可确定绕z轴的塔顶极限弯矩mztt_1。绕z轴的轮毂固定极限弯矩mzhs_1可以通过上文所述的方法得到，之后可以计算得到绕z轴的塔顶极限弯矩mztt_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕x轴的塔顶极限弯矩mxtt_1。下面将介绍确定绕x轴的塔顶极限弯矩mxtt_1的过程。
在一些实施例中，所述基准弯矩包括绕x轴的轮毂固定极限弯矩mxhs_1；所述确定极限弯矩的步骤120，包括：根据关于所述绕x轴的轮毂固定极限弯矩mxhs_1的所述极限弯矩修正函数，确定绕x轴的塔顶极限弯矩mxtt_1。
在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。绕x轴的塔顶极限弯矩mxtt_1的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕x轴的塔顶极限弯矩mxtt_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，关于绕x轴的轮毂固定极限弯矩mxhs_1的绕x轴的塔顶极限弯矩mxtt_1的极限弯矩修正函数如下表达式(43)所示：
mxtt_1＝k6mxhs_1+b6(43)
式中，k6及b6表示拟合系数。
在一些实施例中，拟合系数k6及b6可利用参考数据通过最小二乘法拟合得到。拟合系数k6及b6可以通过最小二乘法线性拟合得到。参考数据包括绕x轴的轮毂固定极限弯矩参考数据和绕x轴的塔顶极限弯矩参考数据。参考数据可以是仿真数据。拟合系数k6及b6的求解过程与拟合误差函数s2中拟合系数k1、b1的求解过程类似，不再进行赘述。
将实际的绕x轴的轮毂固定极限弯矩mxhs_1的值代入表达式(43)，可确定绕x轴的塔顶极限弯矩mxtt_1。绕x轴的轮毂固定极限弯矩mxhs_1可以通过上文所述的方法得到，之后可以计算得到绕x轴的塔顶极限弯矩mxtt_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕y轴的塔底极限弯矩mytb_1。下面将介绍确定绕y轴的塔底极限弯矩mytb_1的过程。
在一些实施例中，所述基准弯矩包括在塔底坐标系下绕y轴的塔底基准弯矩mytb；所述确定极限弯矩的步骤120，包括：根据关于所述绕y轴的塔底基准弯矩mytb的所述极限弯矩修正函数，确定绕y轴的塔底极限弯矩mytb_1。绕y轴的塔底基准弯矩mytb和绕y轴的塔底极限弯矩mytb_1为绕同一坐标轴的弯矩。
在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。绕y轴的塔底极限弯矩mytb_1的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕y轴的塔底极限弯矩mytb_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，绕y轴的塔底基准弯矩mytb的表达式如下表达式(44)：
mytb＝cstd-ytbmytb_2＝cstd-ytb(thtower-massrotorglohg)(44)
式中，t表示均匀风速作用下的推力，massrotor表示风轮质量，lohg表示风轮与塔架中心线之间的纵向距离，g为重力加速度(未考虑机舱在塔架中心线前面部分的质量)；cstd-ytb表示放大系数，类似于上文所述的瞬态弯矩相对于稳态弯矩的放大系数；htower表示塔架高度。
其中，推力t可根据动量理论，采用如下表达式(45)：
式中，ρ表示空气密度，r表示风力发电机的风轮半径，u表示来流风速，可以为均匀风速，ct表示推力系数。来流风速u可采用如下表达式(45)：
式中，p表示风力发电机的额定功率，cp表示功率系数。
结合(44)、(45)及(46)，可得到如下表达式(47)：
至少根据风轮半径、来流风速、塔架高度、风轮质量、风轮到塔架中心的纵向距离和放大系数，计算绕y轴的塔底基准弯矩mytb。利用表达式(47)可计算得到绕y轴的塔底基准弯矩mytb。
在一些实施例中，关于绕y轴的塔底基准弯矩mytb的绕y轴的塔底极限弯矩mytb_1的极限弯矩修正函数如下表达式(48)所示：
mytb_1＝k7mytb+b7(48)
式中，k7及b7表示拟合系数。
在一些实施例中，拟合系数k7及b7可利用参考数据通过最小二乘法拟合得到。绕y轴的塔底极限弯矩mytb_1与绕y轴的塔底基准弯矩mytb为线性关系，拟合系数k7及b7可通过最小二乘法线性拟合得到。参考数据包括绕y轴的塔底基准弯矩参考数据和绕y轴的塔底极限弯矩参考数据。参考数据可以是仿真数据，或者是根据仿真数据计算得到的数据。绕y轴的塔底基准弯矩参考数据可以是依据仿真数据，利用稳态弯矩计算值乘以放大系数得到，绕y轴的塔底极限弯矩参考数据可以是仿真数据。拟合系数k7及b7的求解过程与拟合误差函数s2中拟合系数k1、b1的求解过程类似，不再进行赘述。
将根据表达式(47)计算得到的绕y轴的塔底基准弯矩mytb的值代入表达式(48)，可确定绕y轴的塔底极限弯矩mytb_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕x轴的塔底极限弯矩mxtb_1。下面将介绍确定绕x轴的塔底极限弯矩mxtb_1的过程。
在一些实施例中，所述基准弯矩包括绕y轴的塔底极限弯矩mytb_1；所述确定极限弯矩的步骤120，包括：根据关于所述绕y轴的塔底极限弯矩mytb_1的所述极限弯矩修正函数，确定绕x轴的塔底极限弯矩mxtb_1。
在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。绕x轴的塔底极限弯矩mxtb_1的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕x轴的塔底极限弯矩mxtb_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，关于绕y轴的塔底极限弯矩mytb_1的绕x轴的塔底极限弯矩mxtb_1的极限弯矩修正函数如下统计表达式(49)所示：
mxtb_1＝k8mytb_1+b8(49)
式中，k8及b8表示拟合系数。
在一些实施例中，拟合系数k8及b8可利用参考数据通过最小二乘法拟合得到。绕x轴的塔底极限弯矩mxtb_1和绕y轴的塔底极限弯矩mytb_1为线性关系，拟合系数k8及b8可通过最小二乘法线性拟合得到。参考数据包括绕y轴的塔底极限弯矩参考数据和绕x轴的塔底极限弯矩参考数据。参考数据可以是仿真数据。拟合系数k8及b8的求解过程与上文拟合误差函数s2中拟合系数k1、b1的求解过程类似，不再进行赘述。
将实际的绕y轴的塔底极限弯矩mytb_1的值代入表达式(49)，可确定绕x轴的塔底极限弯矩mxtb_1。绕y轴的塔底极限弯矩mytb_1可以通过上文所述的方法得到，之后可以计算得到绕x轴的塔底极限弯矩mxtb_1。
在一个示例性实施例中，风力发电机的极限弯矩包括绕z轴的塔底极限弯矩mztb_1。下面将介绍确定绕z轴的塔底极限弯矩mztb_1的过程。
在一些实施例中，所述基准弯矩包括绕z轴的塔顶极限弯矩mztt_1；确定极限弯矩的步骤120，包括：根据关于所述绕z轴的塔顶极限弯矩mztt_1的所述极限弯矩修正函数，确定绕z轴的塔底极限弯矩mztb_1。
在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。绕z轴的塔底极限弯矩mztb_1的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，关于基准弯矩的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。绕z轴的塔底极限弯矩mztb_1的极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，关于绕z轴的塔顶极限弯矩mztt_1的绕z轴的塔底极限弯矩mztb_1的极限弯矩修正函数如下表达式(50)所示：
mztb_1＝k9mztt_1+b9(50)
式中，k9及b9表示拟合系数。
在一些实施例中，拟合系数k9及b9可利用参考数据通过最小二乘法拟合得到。绕z轴的塔底极限弯矩mztb_1和绕z轴的塔顶极限弯矩mztt_1是线性关系，拟合系数k9及b9可通过最小二乘法线性拟合得到。参考数据包括绕z轴的塔顶极限弯矩参考数据和绕z轴的塔底极限弯矩参考数据。参考数据可以是仿真数据。拟合系数k9及b9的求解过程与上文所述的拟合误差函数s2中拟合系数k1、b1的求解过程类似，不再进行赘述。
将实际的绕z轴的塔顶极限弯矩mztt_1的值代入表达式(50)，可确定绕z轴的塔底极限弯矩mztb_1。绕z轴的塔顶极限弯矩mztt_1可以通过上文所述的方法得到，之后可以计算得到绕z轴的塔底极限弯矩mztb_1。
综上所述，上文所述任一关于基准弯矩的极限弯矩修正函数可以通过最小二乘法拟合得到。在一些示例性实施例中，所述获得所述风力发电机的基准弯矩，包括：获得所述风力发电机的在叶根坐标系下绕y轴的叶根基准弯矩mybr、绕x轴的叶根基准弯矩mxb和绕x轴的轮毂旋转基准弯矩mxhr中的至少一个。所述根据关于所述基准弯矩的极限弯矩修正函数，确定表征所述极限载荷的极限弯矩，包括：根据对应的所述极限弯矩修正函数，确定所述绕y轴的叶根基准弯矩mybr、所述绕x轴的叶根基准弯矩mxbr和所述绕x轴的轮毂旋转基准弯矩mxhr中的至少一个对应的极限弯矩。根据关于绕y轴的叶根基准弯矩mybr的极限弯矩修正函数，确定绕y轴的叶根极限弯矩mybr_1；根据关于绕x轴的叶根基准弯矩mxbr的极限弯矩修正函数，确定绕x轴的叶根极限弯矩mxbr_1；根据关于绕x轴的轮毂旋转基准弯矩mxhr的极限弯矩修正函数，确定绕x轴的轮毂旋转极限弯矩mxhr_1。该些基准弯矩和对应的极限弯矩为绕同一坐标轴的弯矩。
对于绕y轴的叶根极限弯矩mybr_1、绕x轴的叶根极限弯矩mxbr_1和绕x轴的轮毂旋转极限弯矩mxhr_1，在一些实施例中，对应的极限弯矩修正函数包括关于误差修正函数和所述基准弯矩的函数，误差修正函数至少是通过最小二乘法拟合得到的。在一些实施例中，误差修正函数至少是通过最小二乘法非线性拟合得到的。在一些实施例中，所述误差修正函数至少是利用所述风力发电机的风轮半径参考数据和额定功率参考数据，通过最小二乘法非线性拟合得到的。可以根据所述风力发电机的风轮半径和额定功率，确定所述误差修正函数的值；及根据确定的所述误差修正函数的值和所述基准弯矩，确定所述极限弯矩。在一些实施例中，所述极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
对于绕y轴的叶根极限弯矩mybr_1、绕x轴的叶根极限弯矩mxbr_1和绕x轴的轮毂旋转极限弯矩mxhr_1，在一些实施例中，可以确定稳态弯矩；及根据所述稳态弯矩，确定瞬态弯矩，作为所述基准弯矩。
在一个示例性实施例中，所述基准弯矩包括绕x轴的轮毂旋转极限弯矩mxhr_1、绕y轴的轮毂旋转极限弯矩myhr_1和绕z轴的轮毂旋转极限弯矩mzhr_1。所述根据关于所述基准弯矩的极限弯矩修正函数，确定表征所述极限载荷的极限弯矩，包括：根据所述极限弯矩修正函数，确定绕x轴的轮毂固定极限弯矩mxhs_1、绕y轴的轮毂固定极限弯矩myhs_1和绕z轴的轮毂固定极限弯矩mzhs_1中的至少一个。
对于绕x轴的轮毂固定极限弯矩mxhs_1、绕y轴的轮毂固定极限弯矩myhs_1和/或绕z轴的轮毂固定极限弯矩mzhs_1，在一些实施例中，所述极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，所述极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
在一些实施例中，对于绕y轴的轮毂旋转极限弯矩myhr_1、绕z轴的轮毂旋转极限弯矩mzhr_1、绕y轴的塔顶极限弯矩mytt_1、绕z轴的塔顶极限弯矩mztt_1、绕x轴的塔顶极限弯矩mxtt_1、绕y轴的塔底极限弯矩mytb_1、绕x轴的塔底极限弯矩mxtb_1和/或绕z轴的塔底极限弯矩mztb_1对应的极限弯矩修正函数至少是通过最小二乘法线性拟合得到的。在一些实施例中，极限弯矩修正函数至少是利用仿真数据通过最小二乘法拟合得到的。
为了验证本申请提供的方法确定的极限弯矩的准确度，选择十款不同型号风力发电机，分别采用仿真实验得到极限弯矩，同时采用本申请实施例提供的方法确定极限弯矩。两种方法确定极限弯矩时，输入参数包括风力发电机的额定功率、风轮直径、塔架高度等，得到的极限弯矩包括绕x轴的叶根极限弯矩、绕y轴的叶根极限弯矩以及绕z轴的叶根极限弯矩，在轮毂旋转坐标系下绕x轴的轮毂旋转极限弯矩、绕y轴的轮毂旋转极限弯矩、以及绕z轴的轮毂旋转极限弯矩，绕x轴的轮毂固定极限弯矩、绕y轴的轮毂固定极限弯矩、以及绕z轴的轮毂固定极限弯矩，绕x轴的塔顶极限弯矩、绕y轴的塔顶极限弯矩、以及绕z轴的塔顶极限弯矩，绕x轴的塔底极限弯矩、绕y轴的塔底基准弯矩、以及绕z轴的塔底极限弯矩。
将通过上述两种方法得到的极限弯矩数据进行对比，对比结果如图12至图26所示。图12至图26中，横坐标表示风力发电机的编号，纵坐标表示极限弯矩。其中，图中所示的极限弯矩的值为进行归一化处理后的值。在一些实施例中，在归一化处理时，所有风力发电机的极限弯矩均除以第一组风力发电机对应的极限载荷。
图12所示为通过仿真得到的绕x轴的叶根极限弯矩与通过本申请提供的方法得到的绕x轴的叶根极限弯矩的对比示意图，其中mxbr_s表示通过仿真得到的绕x轴的叶根极限弯矩，mxbr_1表示采用本申请提供的方法得到的绕x轴的叶根极限弯矩。由图12可知，通过仿真得到的绕x轴的叶根极限弯矩、与本申请提供的方法得到的绕x轴的叶根极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕x轴的叶根极限弯矩的准确度较高。
图13所示为通过仿真得到的绕y轴的叶根极限弯矩与通过本申请提供的方法得到的绕y轴的叶根极限弯矩的对比示意图，其中mybr_s表示通过仿真得到的绕y轴的叶根极限弯矩，mybr_1表示采用本申请提供的方法得到的绕y轴的叶根极限弯矩。由图13可知，通过仿真得到的绕y轴的叶根极限弯矩、与本申请提供的方法得到的绕y轴的叶根极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕y轴的叶根极限弯矩的准确度较高。
图14所示为通过仿真得到的绕z轴的叶根极限弯矩与通过本申请提供的方法得到的绕z轴的叶根极限弯矩的对比示意图，其中mzbr_s表示通过仿真得到的绕z轴的叶根极限弯矩，mzbr_1表示采用本申请提供的方法得到的绕z轴的叶根极限弯矩。由图14可知，通过仿真得到的绕z轴的叶根极限弯矩、与本申请提供的方法得到的绕z轴的叶根极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕z轴的叶根极限弯矩的准确度较高。
图15所示为通过仿真得到的绕x轴的轮毂旋转极限弯矩与通过本申请提供的方法得到的绕x轴的轮毂旋转极限弯矩的对比示意图，其中mxhr_s表示通过仿真得到的绕x轴的轮毂旋转极限弯矩，mxhr_1表示采用本申请提供的方法得到的绕x轴的轮毂旋转极限弯矩。由图15可知，通过仿真得到的绕x轴的轮毂旋转极限弯矩、与本申请提供的方法得到的绕x轴的轮毂旋转极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕x轴的轮毂旋转极限弯矩的准确度较高。
图16所示为通过仿真得到的绕y轴的轮毂旋转极限弯矩与通过本申请提供的方法得到的绕y轴的轮毂旋转极限弯矩的对比示意图，其中myhr_s表示通过仿真得到的绕y轴的轮毂旋转极限弯矩，myhr_1表示采用本申请提供的方法得到的绕y轴的轮毂旋转极限弯矩。由图16可知，通过仿真得到的绕y轴的轮毂旋转极限弯矩、与本申请提供的方法得到的绕y轴的轮毂旋转极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕y轴的轮毂旋转极限弯矩的准确度较高。
图17所示为通过仿真得到的绕z轴的轮毂旋转极限弯矩与通过本申请提供的方法得到的绕z轴的轮毂旋转极限弯矩的对比示意图，其中mzhr_s表示通过仿真得到的绕z轴的轮毂旋转极限弯矩，mzhr_1表示采用本申请提供的方法得到的绕z轴的轮毂旋转极限弯矩。由图17可知，通过仿真得到的绕z轴的轮毂旋转极限弯矩、与本申请提供的方法得到的绕z轴的轮毂旋转极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕z轴的轮毂旋转极限弯矩的准确度较高。
图18所示为通过仿真得到的绕x轴的轮毂固定极限弯矩与通过本申请提供的方法得到的绕x轴的轮毂固定极限弯矩的对比示意图，其中mxhs_s表示通过仿真得到的绕x轴的轮毂固定极限弯矩，mxhs_1表示采用本申请提供的方法得到的绕x轴的轮毂固定极限弯矩。由图18可知，通过仿真得到的绕x轴的轮毂固定极限弯矩、与本申请提供的方法得到的绕x轴的轮毂固定极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕x轴的轮毂固定极限弯矩的准确度较高。
图19所示为通过仿真得到的绕y轴的轮毂固定极限弯矩与通过本申请提供的方法得到的绕y轴的轮毂固定极限弯矩的对比示意图，其中myhs_s表示通过仿真得到的绕y轴的轮毂固定极限弯矩，myhs_1表示采用本申请提供的方法得到的绕y轴的轮毂固定极限弯矩。由图19可知，通过仿真得到的绕y轴的轮毂固定极限弯矩、与本申请提供的方法得到的绕y轴的轮毂固定极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕y轴的轮毂固定极限弯矩的准确度较高。
图20所示为通过仿真得到的绕z轴的轮毂固定极限弯矩与通过本申请提供的方法得到的绕z轴的轮毂固定极限弯矩的对比示意图，其中mzhs_s表示通过仿真得到的绕z轴的轮毂固定极限弯矩，mzhs_1表示通过本申请提供的方法得到的绕z轴的轮毂固定极限弯矩。由图20可知，通过仿真得到的绕z轴的轮毂固定极限弯矩、与本申请提供的方法得到的绕z轴的轮毂固定极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕z轴的轮毂固定极限弯矩的准确度较高。
图21所示为通过仿真得到的绕x轴的塔顶极限弯矩与通过本申请提供的方法得到的绕x轴的塔顶极限弯矩的对比示意图，其中mxtt_s表示通过仿真得到的绕x轴的塔顶极限弯矩，mxtt_1表示采用本申请提供的方法得到的绕x轴的塔顶极限弯矩。由图21可知，通过仿真得到的绕x轴的塔顶极限弯矩、与本申请提供的方法得到的绕x轴的塔顶极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕x轴的塔顶极限弯矩的准确度较高。
图22所示为通过仿真得到的的绕y轴的塔顶极限弯矩与通过本申请提供的方法得到的绕y轴的塔顶极限弯矩的对比示意图，其中mytt_s表示通过仿真得到的绕y轴的塔顶极限弯矩，mytt_1表示采用本申请提供的方法得到的绕y轴的塔顶极限弯矩。由图22可知，通过仿真得到的绕y轴的塔顶极限弯矩、与本申请提供的方法得到的绕y轴的塔顶极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕y轴的塔顶极限弯矩的准确度较高。
图23所示为通过仿真得到绕z轴的塔顶极限弯矩与通过本申请提供的方法得到的绕z轴的塔顶极限弯矩的对比示意图，其中mztt_s表示通过仿真得到的绕z轴的塔顶极限弯矩，mztt_1表示通过本申请提供的方法得到的绕z轴的塔顶极限弯矩。由图23可知，通过仿真得到的绕z轴的塔顶极限弯矩、与本申请提供的方法得到的绕z轴的塔顶极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕z轴的塔顶极限弯矩的准确度较高。
图24所示为通过仿真得到的绕x轴的塔底极限弯矩与通过本申请提供的方法得到的绕x轴的塔底极限弯矩的对比示意图，其中mxtb_s表示通过仿真得到的绕x轴的塔底极限弯矩，mxtb_1表示采用本申请提供的方法得到的绕x轴的塔底极限弯矩。由图24可知，通过仿真得到的绕x轴的塔底极限弯矩、与本申请提供的方法得到的绕x轴的塔底极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕x轴的塔底极限弯矩的准确度较高。
图25所示为通过仿真得到的绕y轴的塔底极限弯矩与通过本申请提供的方法得到的绕y轴的塔底极限弯矩的对比示意图，其中mytb_s表示通过仿真得到的绕y轴的塔底极限弯矩，mytb_1表示采用本申请提供的方法得到的绕y轴的塔底极限弯矩。由图25可知，通过仿真得到的绕y轴的塔底极限弯矩、与本申请提供的方法得到的绕y轴的塔底极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕y轴的塔底极限弯矩的准确度较高。
图26所示为通过仿真得到的绕z轴的塔底极限弯矩与通过本申请提供的方法得到的绕z轴的塔底极限弯矩的对比示意图，其中mztb_s表示通过仿真得到的绕z轴的塔底极限弯矩，mztb_1表示通过本申请提供的方法得到的绕z轴的塔底极限弯矩。由图26可知，通过仿真得到的绕z轴的塔底极限弯矩、与本申请提供的方法得到的绕z轴的塔底极限弯矩吻合度较好，说明本申请实施例提供的方法确定的绕z轴的塔底极限弯矩的准确度较高。
综上可知，本申请实施例提供的方法得到的极限弯矩的准确度较高，对于风力发电机的设计具有重要的意义。
本申请实施例还提供了一种风力发电机的极限载荷的确定装置400。参见图27，确定装置400包括一个或多个处理器401，用于实现上述任一实施例所述的确定方法。在一些实施例中，确定装置400可以包括计算机可读存储介质404，计算机可读存储介质可以存储有可被处理器401调用的程序，可以包括非易失性存储介质。在一些实施例中，确定装置400可以包括内存403和接口402。在一些实施例中，确定装置400还可以根据实际应用包括其他硬件。
本申请实施例还提供了一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，实现上述任一实施例所述的确定方法。
本申请可采用在一个或多个其中包含有程序代码的存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。计算机可读存储介质包括永久性和非永久性、可移动和非可移动媒体，可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机可读存储介质的例子包括但不限于：相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。
需要说明的是，在本文中，诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个......”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上对本申请实施例所提供的方法和装置进行了详细介绍，本文中应用了具体个例对本申请的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本申请的方法及其核心思想；同时，对于本领域的一般技术人员，依据本申请的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本申请的限制。