本发明公开了一种基于风荷载的自适应塔筒叶尖净空测量方法，在风机的机舱安装激光雷达，获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度；分别构建叶片的有限元分析模型和塔筒的有限元分析模型，获取任意来流风速下叶片任一位置的形变矢量和塔筒任一位置的形变矢量；计算得到在当前来流风速下叶尖的坐标和激光雷达的坐标；得到激光雷达信号的当前发射角度，获取激光雷达所要旋转的角度；当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，计算得到塔筒叶尖净空。本发明只用一个激光雷达就能够实现任意风荷载下的塔筒叶尖净空的精准测量。
1.一种基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述方法包括步骤：S1、在风机的机舱下方中轴线处安装一激光雷达，并根据塔筒的底部中心点位置和风机的迎风方向确定二维坐标体系；
S2、对风机进行标定，获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度，发射角度指激光雷达信号的发射方向和重力方向的夹角；
S3、根据获取的叶片升力与来流风速的关系、叶片推力与来流风速的关系、叶片的几何信息以及材料信息，构建叶片的有限元分析模型，并获取任意来流风速下叶片任一位置的形变矢量；
S4、根据获取的塔筒的几何信息和材料信息，计算塔筒和机舱连接处接收到的来自风机叶轮的推力，构建塔筒的有限元分析模型，并获取任意来流风速下塔筒任一位置的形变矢量；
S5、根据所述任意来流风速下叶片任一位置的形变矢量，计算得到在当前来流风速下叶尖的坐标，根据所述任意来流风速下塔筒任一位置的形变矢量，计算得到在当前来流风速下激光雷达的坐标；
S6、根据所述当前来流风速下叶尖的坐标和激光雷达的坐标，计算得到在当前来流风速下激光雷达信号的当前发射角度，并获取激光雷达所要旋转的角度，该旋转角度使激光雷达信号对准叶尖；
S7、当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，并计算得到塔筒叶尖净空。
2.如权利要求1所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S1包括：
所述二维坐标体系的坐标原点为塔筒底部中心点， x轴的正方向为风机的叶轮迎风方向，y轴的负方向为重力方向。
3.如权利要求2所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S2包括：
获得无风停机状态下风机的叶尖的初始坐标(xtip0,ytip0)，以及获取激光雷达的初始坐标(xld0,yld0)，计算得到激光雷达信号的初始发射角度θ0为：（1）；
其中，xtip0、ytip0分别为无风停机状态下风机的叶尖的初始x轴坐标和初始y轴坐标，xld0、yld0分别为激光雷达的初始x轴坐标和初始y轴坐标。
4.如权利要求3所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S3包括：
根据获得的风机的叶片升力与来流风速u的关系Cl(u,rc)，计算得到叶片在来流风速u下受到的总升力CL(u)为：
根据获得的叶片推力与来流风速u的关系Ct(u,rc)，叶片在来流风速u下受到的总推力CT(u)为；
其中，rc为风机的叶片某一点到叶根的距离，RC为叶片的叶根到叶尖的长度；
叶片的根部固定要求满足公式（4）的约束条件： （4）；
其中，Δx1(0,u)表示叶片根部在来流风速u下x方向的偏移量，Δy1(0,u)表示叶片根部在来流风速u下y方向的偏移量；
以公式（2）、（3）、（4）所满足的条件作为有限元分析模型的边界条件，构建叶片的有限元分析模型，计算得到任意来流风速u下叶片任一位置的形变矢量X1为： （5）；
其中，Δx1(rc,u)为叶片某一点位置rc在来流风速u下的x方向偏移量，Δy1(rc,u)为叶片某一点位置rc在来流风速u下的y方向偏移量，Gblade叶片的几何信息，Mblade为叶片的材料信息，FEMb表示叶片的有限元分析模型。
5.如权利要求4所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S4包括：
通过查询风机的塔筒的系统设计手册，得到塔筒的几何信息Gtower和材料信息Mtower，计算塔筒和机舱连接处收到来自风机叶轮的推力Ftop： （6）；
其中，CT(u)为单个叶片在来流风速u下受到的总推力；
塔筒的底部固定要求满足公式（7）的约束条件： （7）；
其中，Δx2(0,u)表示塔筒底部在来流风速u下的x方向偏移量，Δy2(0,u)表示塔筒底部在来流风速u下的y方向偏移量；
以公式（6）、（7）所满足的条件作为有限元分析模型的边界条件，构建塔筒的有限元分析模型，计算得到任意来流风速u下塔筒上任一意位置的形变矢量X2为： （8）；
其中，l为塔筒某一点高度，即塔筒某一点到塔筒底部的距离，Δx2(l,u)为塔筒某一点高度l在来流风速u下的x方向偏移量，Δy2(l,u)为塔筒某一点高度l在来流风速u下的y方向偏移量，Gtower塔筒的几何信息，Mtower为塔筒的材料信息，FEMt表示塔筒的有限元分析模型。
6.如权利要求5所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S5包括：
连接风机的主控系统，获取当前来流风速unow，计算当前来流风速下的叶尖坐标Xtip(unow)为：
其中，unow为当前来流风速，xtip(unow)为当前来流风速下叶尖的x轴坐标，ytip(unow)为当前来流风速下叶尖的y轴坐标，RC为叶片的叶根到叶尖的长度，Δx1(Rc, unow)为叶尖在当前来流风速下的x方向偏移量，Δy1(Rc, unow)为叶尖在当前来流风速下的y方向偏移量，xtip0、ytip0分别为无风停机状态下风机的叶尖的初始x轴坐标和初始y轴坐标。
7.如权利要求6所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S5包括：
计算得到当前来流风速下激光雷达坐标为Xld(unow)：（10）；
其中，unow为当前来流风速，xld(unow)为当前来流风速下激光雷达的x轴坐标, yld(unow)为当前来流风速下激光雷达的y轴坐标，H为塔筒的高度，Δx2(H,unow)为塔筒在当前来流风速下的x方向偏移量，Δy2(H,unow)为塔筒在当前来流风速下的y方向偏移量。
8.如权利要求7所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S6包括：
计算得到在当前来流风速下激光雷达信号的当前发射角度θ(unow)，并计算得到激光雷达所需要旋转的角度Δθ；
其中，xtip(unow)为当前来流风速下叶尖的x轴坐标，ytip(unow)为当前来流风速下叶尖的y轴坐标，xld(unow)为当前来流风速下激光雷达的x轴坐标，yld(unow)为当前来流风速下激光雷达的y轴坐标，θ0为激光雷达信号的初始发射角度；
将Δθ输入给激光雷达的伺服电机，以控制激光雷达按照所需要旋转的角度Δθ进行转动，使激光雷达信号对准叶尖。
9.如权利要求8所述的基于风荷载的自适应塔筒叶尖净空测量方法，其特征在于，所述步骤S7包括：
当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离s，并对塔筒进行测量，获取到塔筒根部的距离为h的位置和对应该位置的塔筒横截面半径r的关系为r(h)，计算得到叶片塔筒净空D为： （13）；
其中，θ(unow)为当前来流风速下激光雷达信号的当前发射角度，xld(unow)为当前来流风速下激光雷达的x轴坐标，r(ytip(unow))表示h取值为ytip(unow)时由r(h)得到的塔筒横截面半径，ytip(unow)为当前来流风速下叶尖的y轴坐标。
10.一种基于风荷载的自适应塔筒叶尖净空测量系统，其特征在于，所述系统包括：坐标系模块，用于在风机的机舱下方中轴线处安装一激光雷达，并根据塔筒的底部中心点位置和风机的迎风方向确定二维坐标体系；
第一计算模块，用于对风机进行标定，获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度，发射角度指激光雷达信号的发射方向和重力方向的夹角；
叶片有限元模型模块，用于根据获取的叶片升力与来流风速的关系、叶片推力与来流风速的关系、叶片的几何信息以及材料信息，构建叶片的有限元分析模型，并获取任意来流风速下叶片任一位置的形变矢量；
塔筒有限元模型模块，用于根据获取的塔筒的几何信息和材料信息，计算塔筒和机舱连接处接收到的来自风机叶轮的推力，构建塔筒的有限元分析模型，并获取任意来流风速下塔筒任一位置的形变矢量；
第二计算模块，用于根据所述任意来流风速下叶片任一位置的形变矢量，计算得到在当前来流风速下叶尖的坐标，根据所述任意来流风速下塔筒任一位置的形变矢量，计算得到在当前来流风速下激光雷达的坐标；
旋转模块，用于根据所述当前来流风速下叶尖的坐标和激光雷达的坐标，计算得到在当前来流风速下激光雷达信号的当前发射角度，并获取激光雷达所要旋转的角度，该旋转角度使激光雷达信号对准叶尖；
塔筒叶尖净空模块，当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，并计算得到塔筒叶尖净空。
一种基于风荷载的自适应塔筒叶尖净空测量方法及系统
技术领域
[0001] 本发明涉及长度测量技术领域，尤其涉及一种基于风荷载的自适应塔筒叶尖净空测量方法及系统。
背景技术
[0002] 塔筒叶尖净空是指风力发电机轮毂转动时叶片扫过塔筒时叶尖部位距离塔筒的最小几何距离，叶片迎风受力旋转所形成的旋转曲面形状往往会随风力荷载的变化而变
化，为了避免叶片与塔筒间的运动干涉，在设计时需要通过仿真计算出各种工况下的塔筒
叶尖净空距离，确保风机的运行安全。尽管如此，风力发电机组在运行过程中，由于叶片受
损、传感器故障、控制系统故障或者遭遇极端风况等情况的发生，可能会出现整机振动失
稳，叶片的尖端和塔筒的距离(净空)急剧降低，最终导致叶片与塔筒相碰，造成叶片和塔筒
断裂，该现象又称为“扫塔”。对于风力发电机组而言，一旦发生叶片扫塔则需要更换叶片。
单只叶片的成本较高，同时由于更换叶片期间还会造成风力发电机组停机，因此在更换叶
片期间还会造成发电量损失，这会为风电厂带来较大的经济损失。
[0003] 现有技术中，一般通过照相机采集叶片的图像数据，对图像数据进行分析，得到塔筒叶尖净空，但是照相机会受到采集环境的影响，比如天气不好的时候，不能实时进行采
集。或者通过激光雷达直接测距得到塔筒叶尖净空，但是由于激光雷达安装位置的受限，激
光信号不能照射任意风荷载下的叶尖位置，因此测量得到的塔筒叶尖净空导致不精确。专
利公开号为CN106091941A的风力发电机叶尖塔筒净空的测量方法，该方法在塔筒上安装了
3个2D激光扫描仪以实现对叶片的全方位监测，其缺点是成本过高。专利公开号为
CN113309674A的风力发电机组的净空距离确定方法，该方案中净空距离通过图像处理得
到，净空距离测量不确定度受拍摄机位跟被测风机的相对角度影响很大，而且恶劣天气和
夜间无法测量数据。专利公开号为CN113217283A的风电机组叶片净空控制系统，风机运行
过程中需要根据不同的偏航位置不断的调整拍摄位置，难以保证数据的一致性。
发明内容
[0004] 鉴于此，本发明提供一种基于风荷载的自适应塔筒叶尖净空测量方法及系统，通过计算机辅助技术，只用一个激光雷达就能够实现任意风荷载下的塔筒叶尖净空的精准测
量。
[0005] 为实现上述目的，本发明提供一种基于风荷载的自适应塔筒叶尖净空测量方法，所述方法包括步骤：
[0006] S1、在风机的机舱下方中轴线处安装一激光雷达，并根据塔筒的底部中心点位置和风机的迎风方向确定二维坐标体系；
[0007] S2、对风机进行标定，获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度，发射角度指激光雷达信号的发射方向和重
力方向的夹角；
[0008] S3、根据获取的叶片升力与来流风速的关系、叶片推力与来流风速的关系、叶片的几何信息以及材料信息，构建叶片的有限元分析模型，并获取任意来流风速下叶片任一位
置的形变矢量；
[0009] S4、根据获取的塔筒的几何信息和材料信息，计算塔筒和机舱连接处接收到的来自风机叶轮的推力，构建塔筒的有限元分析模型，并获取任意来流风速下塔筒任一位置的
形变矢量；
[0010] S5、根据所述任意来流风速下叶片任一位置的形变矢量，计算得到在当前来流风速下叶尖的坐标，根据所述任意来流风速下塔筒任一位置的形变矢量，计算得到在当前来
流风速下激光雷达的坐标；
[0011] S6、根据所述当前来流风速下叶尖的坐标和激光雷达的坐标，计算得到在当前来流风速下激光雷达信号的当前发射角度，并获取激光雷达所要旋转的角度，该旋转角度使
激光雷达信号对准叶尖；
[0012] S7、当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，并计算得到塔筒叶尖净空。
[0013] 优选的，所述步骤S1包括：所述二维坐标体系的坐标原点为塔筒底部中心点，x轴的正方向为风机的叶轮迎风方向， y轴的负方向为重力方向。
[0014] 优选的，所述步骤S2包括：
[0015] 获得无风停机状态下风机的叶尖的初始坐标(xtip0,ytip0)，以及获取激光雷达的初始坐标(xld0,yld0)。激光雷达的信号需要指向风机的叶尖，计算得到激光雷达信号的初始发
射角度θ0为：
[0016] （1）；
[0017] 其中，xtip0、ytip0分别为无风停机状态下风机的叶尖的初始x轴坐标和初始y轴坐标，xld0、yld0分别为激光雷达的初始x轴坐标和初始y轴坐标。
[0018] 优选的，所述步骤S3包括：
[0019] 根据获得的风机的叶片升力与来流风速u的关系Cl(u,rc)，计算得到叶片在来流风速u下受到的总升力CL(u)为：
[0020] （2）；
[0021] 根据获得的叶片推力与来流风速u的关系Ct(u,rc)，叶片在来流风速u下受到的总推力CT(u)为：
[0022] （3）；
[0023] 其中，rc为风机的叶片某一点到叶根的距离，RC为叶片的叶根到叶尖的长度；
[0024] 叶片的根部固定要求满足公式（4）的约束条件:
[0025] (4)；
[0026] 其中，Δx1(0,u)表示叶片根部在来流风速u下的x方向偏移量，Δy1(0,u)表示叶片根部在来流风速u下的y方向偏移量；
[0027] 以公式（2）、（3）、（4）所满足的条件作为有限元分析模型的边界条件，构建叶片的有限元分析模型，计算得到任意来流风速u下叶片任一位置的形变矢量X1为：
[0028] （5）；
[0029] 其中，Δx1(rc,u)为叶片某一点位置rc在来流风速u下的x方向偏移量，Δy1(rc,u)为叶片某一点位置rc在来流风速u下的y方向偏移量，Gblade叶片的几何信息，Mblade为叶片的
材料信息，FEMb表示叶片的有限元分析模型。
[0030] 优选的，所述步骤S4包括：
[0031] 通过查询风机的塔筒的系统设计手册，得到塔筒的几何信息Gtower和材料信息Mtower，计算塔筒和机舱连接处收到来自风机叶轮的推力Ftop：
[0032] （6）；
[0033] 其中，CT(u)为单个叶片在来流风速u下受到的总推力；
[0034] 塔筒的底部固定要求满足公式（7）的约束条件：
[0035] （7）；
[0036] 其中，Δx2(0,u)表示塔筒底部在来流风速u下的x方向偏移量，Δy2(0,u)表示塔筒底部在来流风速u下的y方向偏移量；
[0037] 以公式（6）、（7）所满足的条件作为有限元分析模型的边界条件，构建塔筒的有限元分析模型，计算得到任意来流风速u下塔筒上任一意位置的形变矢量X2为：
[0038] （8）；
[0039] 其中，l为塔筒某一点高度，即塔筒某一点到塔筒底部的距离，Δx2(l,u)为塔筒某一点高度l在来流风速u下的x方向偏移量，Δy2(l,u)为塔筒某一点高度l在来流风速u下的
y方向偏移量，Gtower塔筒的几何信息，Mtower为塔筒的材料信息，FEMt表示塔筒的有限元分析
模型。
[0040] 优选的，所述步骤S5包括：
[0041] 获取当前来流风速unow，计算当前来流风速下的叶尖坐标Xtip(unow)为：
[0042] （9）；
[0043] 其中，unow为当前来流风速，xtip(unow)为当前来流风速下叶尖的x轴坐标，ytip(unow)为当前来流风速下叶尖的y轴坐标，RC为叶片的叶根到叶尖的长度，Δx1(Rc, unow)为叶尖在
当前来流风速下的x方向偏移量，Δy1(Rc, unow)为叶尖在当前来流风速下的y方向偏移量，
xtip0、ytip0分别为无风停机状态下风机的叶尖的初始x轴坐标和初始y轴坐标。
[0044] 优选的，所述步骤S5还包括：
[0045] 计算得到当前来流风速下激光雷达坐标为Xld(unow)：
[0046] （10）；
[0047] 其中，unow为当前来流风速，xld(unow)为当前来流风速下激光雷达的x轴坐标, yld(unow)，为当前来流风速下激光雷达的y轴坐标，H为塔筒的高度，Δx2(H,unow)为塔筒在当前
来流风速下的x方向偏移量，Δy2(H,unow)为塔筒在当前来流风速下的y方向偏移量。
[0048] 优选的，所述步骤S6包括：
[0049] 计算得到在当前来流风速下激光雷达信号的当前发射角度θ(unow)，并计算得到激光雷达所需要旋转的角度Δθ；
[0050] （11）；
[0051] （12）；
[0052] 其中，xtip(unow)为当前来流风速下叶尖的x轴坐标，ytip(unow)为当前来流风速下叶尖的y轴坐标，xld(unow)为当前来流风速下激光雷达的x轴坐标, yld(unow)为当前来流风速
下激光雷达的y轴坐标，θ0为激光雷达信号的初始发射角度
[0053] 将Δθ输入给激光雷达的伺服电机，以控制激光雷达按照所需要旋转的角度Δθ进行转动，使得激光雷达信号能够对准叶尖。
[0054] 优选的，所述步骤S7包括：
[0055] 当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离s，并对塔筒进行测量，获取到塔筒根部的距离为h的位置和对应该位置的塔筒横截面半径r的关
系为r(h)，计算得到叶片塔筒净空D为：
[0056] （13）；
[0057] 其中，θ(unow)为当前来流风速下激光雷达信号的当前发射角度，Xld(unow)为当前来流风速下激光雷达x轴坐标，r(ytip(unow))表示h取值为ytip(unow)时，由r(h)得到的塔筒横截
面半径，ytip(unow)为当前来流风速下叶尖的y轴坐标。
[0058] 为实现上述目的，本发明提供一种基于风荷载的自适应塔筒叶尖净空测量系统，所述系统包括：
[0059] 坐标系模块，用于在风机的机舱下方中轴线处安装一激光雷达，并根据塔筒的底部中心点位置和风机的迎风方向确定二维坐标体系；
[0060] 第一计算模块，用于对风机进行标定，获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度，发射角度指激光雷达信号
的发射方向和重力方向的夹角；
[0061] 叶片有限元模型模块，用于根据获取的叶片升力与来流风速的关系、叶片推力与来流风速的关系、叶片的几何信息以及材料信息，构建叶片的有限元分析模型，并获取任意
来流风速下叶片任一位置的形变矢量；
[0062] 塔筒有限元模型模块，用于根据获取的塔筒的几何信息和材料信息，计算塔筒和机舱连接处接收到的来自风机叶轮的推力，构建塔筒的有限元分析模型，并获取任意来流
风速下塔筒任一位置的形变矢量；
[0063] 第二计算模块，用于根据所述任意来流风速下叶片任一位置的形变矢量，计算得到在当前来流风速下叶尖的坐标，根据所述任意来流风速下塔筒任一位置的形变矢量，计
算得到在当前来流风速下激光雷达的坐标；
[0064] 旋转模块，用于根据所述当前来流风速下叶尖的坐标和激光雷达的坐标，计算得到在当前来流风速下激光雷达信号的当前发射角度，并获取激光雷达所要旋转的角度，该
旋转角度使激光雷达信号对准叶尖；
[0065] 塔筒叶尖净空模块，当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，并计算得到塔筒叶尖净空。
[0066] 与现有技术相比，本发明提供的一种基于风荷载的自适应塔筒叶尖净空测量方法及系统，所带来的有益效果为：将激光雷达安装在风机的机舱上，使得对塔筒叶尖净空的测
量不受风机偏航方向的影响；通过计算机辅助技术的方法，根据实时获取的风荷载数据，通
过有限元分析预测塔筒叶尖净空的范围，将激光雷达的角度调整到最佳测量角度，测量得
到更高精度的塔筒叶尖净空，从而实现只用一个激光雷达就能实现任意工况下塔筒叶尖净
空的精准测量；避免了安装多个激光雷达追踪叶轮的需求，节省了用户成本。
附图说明
[0067] 图1是根据本发明的一个实施例的基于风荷载的自适应塔筒叶尖净空测量方法的流程示意图。
[0068] 图2是根据本发明的一具体实施例的二维坐标体系的示意图。
[0069] 图3是根据本发明的一个实施例的基于风荷载的自适应塔筒叶尖净空测量系统的系统示意图。
具体实施方式
[0070] 以下将结合附图所示的具体实施方式对本发明进行详细描述，但这些实施方式并不限制本发明，本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的
变换均包含在本发明的保护范围内。
[0071] 如图1所示的本发明的一个实施例，本发明提供一种基于风荷载的自适应塔筒叶尖净空测量方法，所述方法包括步骤：
[0072] S1、在风机的机舱下方中轴线处安装一激光雷达，并根据塔筒的底部中心点位置和风机的迎风方向确定二维坐标体系；
[0073] S2、对风机进行标定，获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度，发射角度指激光雷达信号的发射方向和重
力方向的夹角；
[0074] S3、根据获取的叶片升力与来流风速的关系、叶片推力与来流风速的关系、叶片的几何信息以及材料信息，构建叶片的有限元分析模型，并获取任意来流风速下叶片任一位
置的形变矢量；
[0075] S4、根据获取的塔筒的几何信息和材料信息，计算塔筒和机舱连接处接收到的来自风机叶轮的推力，构建塔筒的有限元分析模型，并获取任意来流风速下塔筒任一位置的
形变矢量；
[0076] S5、根据所述任意来流风速下叶片任一位置的形变矢量，计算得到在当前来流风速下叶尖的坐标，根据所述任意来流风速下塔筒任一位置的形变矢量，计算得到在当前来
流风速下激光雷达的坐标；
[0077] S6、根据所述当前来流风速下叶尖的坐标和激光雷达的坐标，计算得到在当前来流风速下激光雷达信号的当前发射角度，并获取激光雷达所要旋转的角度，该旋转角度使
激光雷达信号对准叶尖；
[0078] S7、当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，并计算得到塔筒叶尖净空。
[0079] 在风机的机舱下方中轴线处安装一激光雷达，并在该激光雷达上安装伺服电机，用以调整激光雷达的信号的发射角度。机舱始终随着叶轮一起旋转。图2所示的二维坐标体
系，根据风机塔筒的底部中心点位置和风机的迎风方向确定一个二维坐标体系，其中，该二
维坐标体系的坐标原点为塔筒底部中心点，由于叶轮始终迎风，x轴的正方向为风机的叶轮
迎风方向，即来流风向的反方向，y轴的负方向为重力方向。本实施例中以激光雷达为例，除
了激光雷达，超声波雷达、毫米波雷达等也是本发明的保护范围之内。
[0080] 获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度，发射角度指激光雷达信号的发射方向和重力方向的夹角。在无风
停机状态下，对风机系统进行初始标定，获得无风停机状态下风机的叶尖的初始坐标
(xtip0,ytip0)，以及获取激光雷达的初始坐标(xld0,yld0)。激光雷达的信号需要指向风机的叶
尖，计算得到激光雷达信号的初始发射角度θ0为：
[0081] （1）；
[0082] 其中，xtip0、ytip0分别为无风停机状态下风机的叶尖的初始x轴坐标和初始y轴坐标，xld0、yld0分别为激光雷达的初始x轴坐标和初始y轴坐标。
[0083] 根据获取的叶片升力与来流风速的关系、叶片推力与来流风速的关系、叶片的几何信息以及材料信息，构建叶片的有限元分析模型，并获取任意来流风速下叶片任一位置
的形变矢量。通过查询风机叶片的系统设计手册，得到风机的叶片升力与来流风速u的关系
Cl(u,rc)、叶片推力与来流风速u的关系Ct(u,rc)，以及获取叶片的几何信息Gblade和材料信
息Mblade。
[0084] 根据获得的风机的叶片升力与来流风速u的关系Cl(u,rc)，计算得到叶片在来流风速u下受到的总升力CL(u)为：
[0085] （2）；
[0086] 其中，rc为风机的叶片某一点到叶根的距离，RC为叶片的叶根到叶尖的长度；
[0087] 根据获得的叶片推力与来流风速u的关系Ct(u,rc)，叶片在来流风速u下受到的总推力CT(u)为：
[0088] （3）；
[0089] 其中，rc为风机的叶片某一点到叶根的距离，RC为叶片的叶根到叶尖的长度。
[0090] 叶片的根部固定要求满足公式（4）的约束条件:
[0091] (4)；
[0092] 其中，Δx1(0,u)表示叶片根部在来流风速u下的x方向偏移量，Δy1(0,u)表示叶片根部在来流风速u下的y方向偏移量。
[0093] 以公式（2）、（3）、（4）所满足的条件作为有限元分析模型的边界条件，构建叶片的有限元分析模型，并根据叶片的有限元分析模型，计算得到任意来流风速u下叶片任一位置
的形变矢量X1为：
[0094] （5）；
[0095] 其中，Δx1(rc,u)为叶片某一点位置rc在来流风速u下的x方向偏移量，Δy1(rc,u)为叶片某一点位置rc在来流风速u下的y方向偏移量，Gblade叶片的几何信息，Mblade为叶片的
材料信息，FEMb表示叶片的有限元分析模型，可以使用任意有限元软件，如ANSYS等。
[0096] 根据获取的塔筒的几何信息和材料信息，计算塔筒和机舱连接处接收到的来自风机叶轮的推力，构建塔筒的有限元分析模型，并获取任意来流风速下塔筒任一位置的形变
矢量。通过查询风机的塔筒的系统设计手册，得到塔筒的几何信息Gtower和材料信息Mtower，
计算塔筒和机舱连接处收到来自风机叶轮的推力Ftop：
[0097] （6）；
[0098] 其中，CT(u)为单个叶片在来流风速u下受到的总推力，
[0099] 塔筒的底部固定要求满足公式（7）的约束条件：
[0100] （7）；
[0101] 其中，Δx2(0,u)表示塔筒底部在来流风速u下的x方向偏移量，Δy2(0,u)表示塔筒底部在来流风速u下的y方向偏移量。
[0102] 以公式（6）、（7）所满足的条件作为有限元分析模型的边界条件，构建塔筒的有限元分析模型，并根据塔筒的有限元分析模型，计算得到任意来流风速u下塔筒上任一意位置
的形变矢量X2为：
[0103] （8）；
[0104] 其中，l为塔筒某一点到塔筒底部的距离，Δx2(l,u)为塔筒某一点高度l在来流风速u下的x方向偏移量，Δy2(l,u)为塔筒某一点高度l在来流风速u下的y方向偏移量，Gtower
塔筒的几何信息，Mtower为塔筒的材料信息，FEMt表示塔筒的有限元分析模型。
[0105] 根据所述任意来流风速下叶片任一位置的形变矢量，计算得到在当前来流风速下叶尖的坐标，根据所述任意来流风速下塔筒任一位置的形变矢量，计算得到在当前来流风
速下激光雷达的坐标。连接风机的主控系统，获取当前来流风速unow，计算当前来流风速下
的叶尖坐标Xtip(unow)为：
[0106] （9）；
[0107] 其中，unow为当前来流风速，xtip(unow)为当前来流风速下叶尖的x轴坐标，ytip(unow)为当前来流风速下叶尖的y轴坐标，RC为叶片的叶根到叶尖的长度，Δx1(Rc, unow)为叶尖在
当前来流风速下的x方向偏移量，Δy1(Rc, unow)为叶尖在当前来流风速下的y方向偏移量，
xtip0、ytip0分别为无风停机状态下风机的叶尖的初始x轴坐标和初始y轴坐标。
[0108] 根据所述任意来流风速下塔筒任一位置的形变矢量，计算得到在当前来流风速下激光雷达的坐标。由于激光雷达安装在机舱上，所处位置位于塔筒顶部，因此计算得到当前
来流风速下激光雷达坐标为Xld(unow)：
[0109] （10）；
[0110] 其中，unow为当前来流风速，xld(unow)为当前来流风速下激光雷达的x轴坐标, yld(unow)为当前来流风速下激光雷达的y轴坐标，H为塔筒的高度，Δx2(H,unow)为塔筒在当前
来流风速下的x方向偏移量，Δy2(H,unow)为塔筒在当前来流风速下的y方向偏移量。
[0111] 根据所述当前来流风速下叶尖的坐标和激光雷达的坐标，计算得到在当前来流风速下激光雷达信号的当前发射角度，并获取激光雷达所要旋转的角度，该旋转角度使激光
雷达信号对准叶尖。计算得到在当前来流风速下激光雷达信号的当前发射角度θ(unow)，并
计算得到激光雷达所需要旋转的角度Δθ；将Δθ输入给激光雷达的伺服电机，以控制激光
雷达按照所需要旋转的角度Δθ进行转动，使得激光雷达信号能够对准叶尖；
[0112] （11）；
[0113] （12）；
[0114] 其中，xtip(unow)为当前来流风速下叶尖的x轴坐标，ytip(unow)为当前来流风速下叶尖的y轴坐标，xld(unow)为当前来流风速下激光雷达的x轴坐标, yld(unow)为当前来流风速
下激光雷达的y轴坐标，θ0为激光雷达信号的初始发射角度。
[0115] 当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，并计算得到塔筒叶尖净空。当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处
的距离s，并对塔筒进行测量，获取到塔筒根部的距离为h的位置和对应该位置的塔筒横截
面半径r的关系为r(h)，计算得到叶片塔筒净空D为；
[0116] （13）；
[0117] 其中，θ(unow)为当前来流风速下激光雷达信号的当前发射角度，Xld(unow)为当前来流风速下激光雷达x轴坐标，r(ytip(unow))表示h取值为ytip(unow)时，由r(h)得到的塔筒横截
面半径，ytip(unow)为当前来流风速下叶尖的y轴坐标。
[0118] 如图3所示的本发明的一实施例，本发明提供一种基于风荷载的自适应塔筒叶尖净空测量系统，所述系统包括：
[0119] 坐标系模块30，用于在风机的机舱下方中轴线处安装一激光雷达，并根据塔筒的底部中心点位置和风机的迎风方向确定二维坐标体系；
[0120] 第一计算模块31，用于对风机进行标定，获得无风停机状态下风机的叶尖的初始坐标和激光雷达的初始坐标，计算激光雷达信号的初始发射角度，发射角度指激光雷达信
号的发射方向和重力方向的夹角；
[0121] 叶片有限元模型模块32，用于根据获取的叶片升力与来流风速的关系、叶片推力与来流风速的关系、叶片的几何信息以及材料信息，构建叶片的有限元分析模型，并获取任
意来流风速下叶片任一位置的形变矢量；
[0122] 塔筒有限元模型模块33，用于根据获取的塔筒的几何信息和材料信息，计算塔筒和机舱连接处接收到的来自风机叶轮的推力，构建塔筒的有限元分析模型，并获取任意来
流风速下塔筒任一位置的形变矢量；
[0123] 第二计算模块34，用于根据所述任意来流风速下叶片任一位置的形变矢量，计算得到在当前来流风速下叶尖的坐标，根据所述任意来流风速下塔筒任一位置的形变矢量，
计算得到在当前来流风速下激光雷达的坐标；
[0124] 旋转模块35，用于根据所述当前来流风速下叶尖的坐标和激光雷达的坐标，计算得到在当前来流风速下激光雷达信号的当前发射角度，并获取激光雷达所要旋转的角度，
该旋转角度使激光雷达信号对准叶尖；
[0125] 塔筒叶尖净空模块36，当激光雷达信号对准叶尖时测量得到从叶尖到激光雷达信号发射处的距离，并计算得到塔筒叶尖净空。
[0126] 尽管为示例目的，已经公开了本发明的优选实施方式，但是本领域的普通技术人员将意识到，在不脱离由所附的权利要求书公开的本发明的范围和精神的情况下，各种改
进、增加以及取代是可能的。
