风电机组塔筒法兰连接螺栓设计方法及装置与流程

本发明涉及风力发电技术领域，特别是涉及一种风电机组塔筒法兰连接螺栓设计方法及装置。

背景技术：

目前我们设计风力发电机组塔筒法兰连接螺栓都是采用工程算法与经验相结合的办法，在进行此种工程计算时，需要对一些几何参数近似估计，而各部件之间的力的作用关系也是依赖理想的抽象数学模型。这就会产生很大的不确定性，大大降低结构的可靠性。因此，为解决上述问题，需要建立新的数值模拟方法计算螺栓的应力，使得工程人员不再只依靠理论公式，将风力发电机组塔筒法兰连接螺栓的不确定性降至最低。

而目前采用有限元计算塔筒法兰连接螺栓的方法都是先进行评估，然后通过solidworks生成三维模型，然后导出至ansys，接着进行有限元的建模和分析。如果计算一次通过，此种方法没有明显的不足。但是，如果经过多次计算，结果仍然不能够满足设计的需求。就需要频繁地更改结构的几何尺寸，频繁地生成三维模型，并导出模型，再投入计算。这会耗费大量的时间与精力，严重影响项目的推进。同时，在繁琐的重复中，会导致操作人员的失误增加，增加结构失效的可能性。所以，现在急需一种高效的有限元数值模拟方法，以弥补上述方法的缺陷。

技术实现要素：

本发明要解决的技术问题是提供一种风电机组塔筒法兰连接螺栓设计方法及装置，从而以较低的计算代价获得精确的连接螺栓的最大疲劳损伤。

为解决上述技术问题，本发明提供了一种风电机组塔筒法兰连接螺栓设计方法，所述方法包括：在ansys软件中对塔筒、塔筒法兰及连接螺栓分别进行建模，并对连接螺栓施加相应的载荷；求解在施加载荷下，塔筒截面弯矩与螺栓应力的分段线性插值函数；根据所述分段线性插值函数，求解螺栓在预定的载荷时间序列下，每个截面上的最大疲劳损伤。

在一些实施方式中，在求解在施加载荷下，塔筒截面弯矩与螺栓应力的分段线性插值函数之后，根据所述分段线性插值函数，求解螺栓在预定的载荷时间序列下，每个截面上的最大疲劳损伤之前，还包括：利用apdl语言，将所有螺栓在应力敏感位置的应力提取出来；将螺栓的应力导出至matlab进行数据分析，求出最大应力幅值，代入到极限应力的求解公式，求出螺栓的极限应力。

在一些实施方式中，所述极限应力的求解公式如下：

其中，z代表极限应力，mb,h代表在法兰位置的最大mx,y力矩，fz,h代表在法兰位置未因数化fz载荷，rtower,h代表在法兰位置处的塔筒半径，nbolts代表连接的螺栓数量，γf,favorable代表最优载荷的部分安全因数。

在一些实施方式中，在根据所述分段线性插值函数，求解螺栓在预定的载荷时间序列下，每个截面上的最大疲劳损伤之后，还包括：将螺栓的最大疲劳损伤输出至excel表格中，绘制相应的疲劳损伤散点图；结合s-n曲线，获得螺栓的局部疲劳载荷安全系数。

在一些实施方式中，求解在施加载荷下，塔筒截面弯矩与螺栓应力的分段线性插值函数，包括：求解在施加载荷下，塔筒截面弯矩与单个螺栓轴向力的分段线性插值函数；根据塔筒截面弯矩与单个螺栓轴向力的分段线性插值函数，确定塔筒截面弯矩与螺栓应力的分段线性插值函数。

在一些实施方式中，塔筒截面弯矩与螺栓应力的分段线性插值函数如下：

其中，σbolt代表螺栓应力，as代表螺栓应力面积。

在一些实施方式中，根据所述分段线性插值函数，求解螺栓在预定的载荷时间序列下，每个截面上的最大疲劳损伤，包括：根据外部输入参数，生成载荷时间序列；将所述载荷时间序列代入所述分段线性插值函数，求得螺栓每个截面上的最大疲劳损伤。

在一些实施方式中，所述外部输入参数包括：风场等级、湍流系数、轮毂中心高度。

在一些实施方式中，在ansys软件中对塔筒、塔筒法兰及连接螺栓分别进行建模，包括：基于二维cad图纸，归纳出建立塔筒法兰连接螺栓三维模型的点的坐标；在ansys经典模块中读取坐标文件，生成关键点，编辑apdl语言，运用这些关键点生成相应的线、面、体，对体进行网格划分，在与螺母接触的区域，尽量做到网格匹配，直至生成仅包含实体185单元的塔筒法兰模型，删除实体特征；在模型中建立连接螺栓，选取的单元是beam188，建立用来施加载荷的载荷伞，选取的单元是beam188；建立各个部件之间的接触，选取的单元类型是conta174和targe170。

此外，本发明还提供了一种风电机组塔筒法兰连接螺栓设计装置，所述装置包括：一个或多个处理器；存储装置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现根据前文所述的风电机组塔筒法兰连接螺栓设计方法。

采用这样的设计后，本发明至少具有以下优点：

本发明提出的风电机组塔筒法兰连接螺栓设计方法及装置，使用apdl语言对风力发电机组塔筒法兰连接螺栓进行参数化建模，考虑了各部件之间的连接关系，真实反映了螺栓在塔筒截面弯矩作用下的力学性能，给出了螺栓的极限应力和20年工作期限内的疲劳损伤，为塔筒的法兰设计提供了可靠的依据，相比工程算法和传统的有限元计算，本方法可靠性高，计算效率提升明显，对专业的要求不高，即使非专业有限元计算人员，也可在只修改几个参数的前提下，进行计算，充分达到了智能化、精细化的设计能力。

附图说明

上述仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，以下结合附图与具体实施方式对本发明作进一步的详细说明。

图1是本发明实施例提供的风电机组塔筒法兰连接螺栓设计方法的流程图；

图2是本发明实施例提供的风电机组塔筒法兰连接螺栓设计方法的流程图；

图3是本发明实施例提供的塔筒的有限元模型示意图；

图4是本发明实施例提供的塔筒连接结构的示意图；

图5是本发明实施例提供的连接螺栓的空间分布示意图；

图6a是本发明实施例提供的塔筒截面弯矩与螺栓轴向力的分段线性插值函数；

图6b是本发明实施例提供的塔筒截面弯矩与螺栓应力的分段线性插值函数；

图7是本发明实施例提供的应力敏感位置的空间分布示意图；

图8是本发明实施例提供的风电机组塔筒法兰连接螺栓设计装置的结构图。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

图1示出了一种实施方式中风电机组塔筒法兰连接螺栓设计方法的流程图。参见图1，风电机组塔筒法兰连接螺栓设计方法包括：

s11，在ansys软件中对塔筒、塔筒法兰及连接螺栓分别进行建模，并对连接螺栓施加相应的载荷。

s12，求解在施加载荷下，塔筒截面弯矩与螺栓应力的分段线性插值函数。

s13，根据所述分段线性插值函数，求解螺栓在预定的载荷时间序列下，每个截面上的最大疲劳损伤。

采用分段线性插值函数，充分考虑了螺栓的分段取值特性，使得对塔筒截面弯矩与螺栓应力之间的函数关系的建模更为精细化，使得对于最大疲劳损伤的计算更为精细化。

图2示出了另一种实施方式中，风电机组塔筒法兰连接螺栓设计方法的流程图。参见图2，风电机组塔筒法兰连接螺栓设计方法包括：

s201，基于二维cad图纸，归纳出建立塔筒法兰连接螺栓三维模型的点的坐标。

s202，在ansys经典模块中读取坐标文件，生成关键点，编辑apdl语言，运用这些关键点生成相应的线、面、体，对体进行网格划分，在与螺母接触的区域，尽量做到网格匹配，直至生成仅包含实体185单元的塔筒法兰模型，删除实体特征。

图3及图4分别示出了塔筒及法兰连接结构的空间模型。

s203，在模型中建立连接螺栓，选取的单元是beam188，建立用来施加载荷的载荷伞，选取的单元是beam188。

图5示出了连接螺栓的有限元模型。

s204，建立各个部件之间的接触，选取的单元类型是conta174和targe170。

s205，依据相关的螺栓预紧力，对螺栓施加与截面尺寸相对应的预紧力，施加的位置选取在螺栓光杆部分的几何中心，方向沿螺杆的轴向。

s206，在载荷伞中心处施加一定数值的弯矩，该弯矩的绝对值要大于该截面位置的最大弯矩。

s207，依据抽象出来的梁单元理论，将法兰看成一根简支梁，求出塔筒截面弯矩与单个螺栓轴向力的分段线性插值函数，进而获得塔筒截面弯矩与螺栓应力的分段线性插值函数。

塔筒截面弯矩与单个螺栓轴向力的分段线性插值函数如下：

z＝0:f′vs＝fv

z＝zi:f′vs＝fv+p×zi

z＝zii:f′vs＝λ^*×zii

zarb＞＞zii:f′vs＝λ^*×zarb

其中，f′vs代表螺栓载荷，z代表螺栓截面载荷，zarb代表任意截面载荷，足够“捕捉”载荷谱。

为了能够确保弯矩均值总是与弯矩被正确、保守的使用，以下要点需要为连接的“压缩面”而计算：

z＝-zif′vs＝fv+p×(-zi)

z＝-ziif′vs＝2×fv+λ^*×(-zii)

zarb＜＜-ziif′vs＝2×fv+λ^*×zarb

其中，用来疲劳分析的非线性关系是弯矩与螺栓应力之间的关系。

螺栓应力按照如下公式计算：

其中，σbolt代表螺栓应力，as代表螺栓应力面积。

承载了大小等于f′vs的载荷的塔筒中的弯矩按照如下公式计算：

其中，mb,h代表在法兰位置的最大mx,y力矩；rtower,h代表在法兰位置的塔筒半径；nbolts代表连接中螺栓数量。

图6a示出了通过上述过程计算得到的塔筒截面弯矩与螺栓轴向力的分段线性插值函数；图6b示出了通过上述过程计算得到的塔筒截面弯矩与螺栓应力的分段线性插值函数。

s208，凭借apdl语言，将所有螺栓在应力敏感位置的应力提取出来。

s209，将螺栓的应力导出至matlab进行数据分析，求出最大应力幅值，代入到极限应力的求解公式，求出螺栓的极限应力。

贡献了每个螺栓的螺栓应力的壁载荷按照如下公式计算：

其中，mb,h代表在法兰位置的最大mx,y力矩；fz,h代表在法兰位置未因数化fz载荷；rtower,h代表在法兰位置处的塔筒半径；nbolts代表连接的螺栓数量；γf,favorable代表最优载荷的部分安全因数。

应该注意，z代表最高螺栓载荷下的壁载荷，也就是极限载荷。计算极限载荷的意义在于，在后续计算最大疲劳损伤过程中，极限载荷给出了对分段线性插值函数的取值范围。

s210，根据风场等级、湍流系数、轮毂中心高度等外部输入参数，基于bladed载荷生成软件，输出风力发电机组在一年中的载荷时间序列。

s211，将载荷时间序列输入matlab,与已有的螺栓应力进行一定的数据后处理，找出这一组螺栓中每个截面上最大疲劳损伤，乘以20，得到该组螺栓在20年中的最大疲劳损伤。

求取上述每个截面上最大疲劳损伤的过程是分段线性插值函数的应用过程。

s212，将螺栓的最大疲劳损伤输出至excel表格中，绘制相应的疲劳损伤散点图。

s213，结合gl2010风力机认证规范中的s-n曲线，获得螺栓的局部疲劳载荷安全系数。

为了对螺栓的全面疲劳分析，对螺栓圆周上各个点的疲劳损伤进行评估十分必要。选取的螺栓在其圆周上每隔30而被评估。应力在ghbladed中每隔时间间隔t通过如下公式计算：

σ(β,t)＝σaxial(t)+σbending_1(t)·cosβ+σbending_2(t)·sinβ

其中，β＝[0°,30°,60°,…,330°]，σaxial(t)，σbending_1(t)及σbending_2(t)全部为预先计算的量。其中，σaxial(t)为由预紧力和轴向力产生，σbending_1(t)为由力矩mx产生的，σbending_2(t)为由力矩my产生的。

图8示出了本发明风电机组塔筒法兰连接螺栓设计装置的结构图。参见图8，风电机组塔筒法兰连接螺栓设计装置包括：中央处理单元(cpu)801，其可以根据存储在只读存储器(rom)中的程序或者从存储部分808加载到随机访问存储器(ram)803中的程序而执行各种适当的动作和处理。在ram803中，还存储有系统操作所需的各种程序和数据。cpu801、rom802以及ram803通过总线804彼此相连。输入/输出(i/o)接口805也连接至总线804。

以下部件连接至i/o接口805：包括键盘、鼠标等的输入部分806；包括诸如阴极射线管(crt)、液晶显示器(lcd)等以及扬声器等的输出部分807；包括硬盘等的存储部分808；以及包括诸如lan卡、调制解调器等的网络接口卡的通信部分809。通信部分809经由诸如因特网的网络执行通信处理。驱动器810也根据需要连接至i/o接口805。可拆卸介质811，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器810上，以便于从其上读出的计算机程序根据需要被安装入存储部分808。

特别的，根据本发明实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括承载在计算机可读介质上的计算机程序，该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分809从网络上被下载和安装，和/或从可拆卸介质811被安装。在该计算机程序被中央处理单元(cpu)801执行时，执行本发明的方法中限定的上述功能。需要说明的是，本发明的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意结合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、光纤、便携式紧凑磁盘只读存储器(cd-rom)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本发明中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者结合使用。而在本发明中，计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何恰当的介质传输，包括但不限于：无线、电线、光缆、rf等等，或者上述的任意合适的组合。

附图中的流程图和框图，图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连的表示的方框实际上可以基本并行的执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰，均落在本发明的保护范围内。