一种基于风力机叶片有限元模型的分布力加载方法技术

本发明专利技术提出一种基于风力机叶片有限元模型的分布力加载方法，涉及风力机叶片结构受力分析的加载方法的技术领域。其特征在于：将全尺寸叶片展向离散为M个叶片剖面，以展向段N的有限元模型为研究对象；根据组成展向段N的叶片n截面和n+1截面的剖面载荷（力、力矩及作用点），按照力和力矩等效原则，基于展向段的壳体SS面和PS面的节点坐标和节点个数，将截面集中力载荷转化为作用在展向段N的壳体主梁区域的分布力载荷。本发明专利技术为风力机叶片结构受力分析的高精准加载方式提供了有效的方法指导，具有极大的工程应用价值。极大的工程应用价值。极大的工程应用价值。

【技术实现步骤摘要】
一种基于风力机叶片有限元模型的分布力加载方法


[0001]本专利技术涉及个基于风力机叶片有限元模型的分布力加载方法，属于风电叶片结构受力分析


技术介绍

[0002]风力机叶片作为风力发电机组的关键捕风结构，对机组的稳定、安全运行具有决定性的影响，研究风力机叶片结构承载能力的一项重要内容就是加载方式。由于环境作用、风况工况条件以及细长梁结构形式等影响，叶片服役期内受疲劳载荷和极限载荷的耦合作用，在叶片设计分析中需要模拟近3000个载荷工况条件才能完成叶片载荷计算。众多载荷形式最终转化为风力机各个部件坐标系下的力和力矩。目前，在风力机叶片设计标准(DNV
‑
GL
‑
0376、或TUV等)中，叶片设计和认证的加载形式多为沿叶片展向离散为几十个剖面后，采用剖面多点约束的集中力加载方法，或在全尺寸叶片静载测试中，也是采用四点
‑
六点的剖面集中力加载方式，这种加载方式无法真实模拟风力机叶片在捕风中受到的作用于叶片外表面的分布力载荷，又容易在集中力加载剖面引起应力集中，降低结构校核精度和准度。基于以上情况，提出一种基于风力机叶片有限元模型的分布力加载方法,是具有重大工程应用价值和工程意义的技术问题，也为叶片结构设计和结构校核提供了一种新的解决方法。

技术实现思路

[0003]本专利技术针对现有技术中的不足，提供一种基于风力机叶片有限元模型的分布力加载方法。
[0004]为实现上述技术目的，本专利技术采取的技术方案为：
[0005]一种基于风力机叶片有限元模型的分布力加载方法，包括以下步骤：
[0006]步骤1、加载基于风力机叶片叶根坐标系，
[0007]步骤2、将风力机叶片展向进行截面划分，以叶根至叶尖方向将叶片划分为m个截面，叶根处截面标记为1，
[0008]步骤3、将叶片展向划分的截面n和截面n+1组成的叶片展向段，标记为N，此时基于叶根坐标系，叶片每个剖面的截面载荷为载荷作用点为其中，M表示剖面的弯矩，F表示的剪力，X、Y、Z分别为载荷作用点在叶根坐标系中的坐标值，
[0009]步骤4、以展向段N为叶片有限元模型研究对象，叶片有限元模型包括SS面壳体、PS面壳体和腹板，展向段N内标记SS面壳体和PS面壳体若干个壳体节点，其中，PS面壳体和SS面壳体主梁区域的节点数目分别记为nump，nums，
[0010]步骤5、通过截面n和n+1，计算展向段N内所有壳体节点的总节点合力以及载
荷作用点坐标偏差，其中，为载荷作用点坐标偏差，X
n
,Y
n
,Z
n
分别表示n截面的X、Y、Z轴载荷作用点坐标，X
n+1
,Y
n+1
,Z
n+1
分别表示n+1截面的X、Y、Z轴载荷作用点坐标，同时标记展向段N内合力作用点Loc
n
为截面n的载荷作用点位置，
[0011]步骤6、计算展向段N内PS面壳体和SS面壳体的所有壳体节点的坐标之和，分别标记为：
[0012]∑NodesX
P
∑NodesY
P
∑NodesZ
P
，
[0013]∑NodesX
S
∑NodesY
S
∑NodesZ
S
，
[0014]其中，∑NodesX
P
为PS面壳体的所有壳体节点的X轴坐标值之和，∑NodesY
P
为PS面壳体的所有壳体节点的Y轴坐标值之和，∑NodesZ
P
为PS面壳体的所有壳体节点的Z轴坐标值之和，∑NodesX
S
为SS面壳体的所有壳体节点的X轴坐标值之和，∑NodesY
S
为SS面壳体的所有壳体节点的Y轴坐标值之和，∑NodesZ
S
为SS面壳体的所有壳体节点的Z轴坐标值之和，
[0015]步骤7、计算展向段N内有限元模型的PS面壳体和SS面壳体的所有壳体节点坐标，分别标记为：
[0016]和
[0017]其中，dxp、dyp、dzp为PS面壳体节点X、Y和Z方向的坐标值，dxs、dys、dzs为SS面壳体节点X、Y和Z方向的坐标值，
[0018]步骤8、将展向段N内所有壳体节点的合力以节点坐标均匀分布的方式，基于力和力矩同时等效的原则，计算展向段N内节点的分布力转换矩阵[ΔC]，其中，ΔM为壳体节点的弯矩合力，ΔF为壳体节点的剪力合力,
[0019]步骤9、计算展向段N内PS面壳体和SS面壳体的每个节点的分布力载荷；
[0020]步骤10、重复步骤5至步骤9的方法，进行下一截面的计算，直至m个截面的每个节点的分布力载荷全部计算完成，根据每个节点的分布力载荷模拟风力机叶片真实风况载荷。
[0021]为优化上述技术方案，本专利技术采取的进一步的改进方式如下：
[0022]步骤1所述的风力机叶片叶根坐标系包括3个弯矩M
X
、M
Y
、M
Z
和3个剪力F
X
、F
Y
、F
Z
。
[0023]步骤5所述的展向段N内所有节点合力的计算方法为：
[0024][0025]其中，M
n
为n截面所有节点的弯矩，M
n+1
为n+1截面所有节点的弯矩，F
n
为n截面所有
节点的剪力，F
n+1
为n+1截面所有节点的剪力，载荷坐标的偏差的计算方法为其中，X
n
,Y
n
,Z
n
分别表示n截面的X、Y、Z轴载荷作用点坐标，X
n+1
,Y
n+1
,Z
n+1
分别表示n+1截面的X、Y、Z轴载荷作用点坐标。
[0026]展向段N内节点合力展开计算为：
[0027]ΔM
X
＝M
Xn
‑
M
Xn+1
‑
(ΔY
×
F
Zn+1
‑
ΔZ
×
F
Yn+1
)
[0028]ΔM
Y
＝M
Yn
‑
M
Yn+1
‑
(ΔZ
×
F
Xn+1
‑
ΔX
×
F
Zn+1
)
[0029]ΔM
Z
＝M
Zn
‑
M
Zn+1
‑
(ΔX
×
F
Yn+1
‑
ΔY
×
F
Xn+1
)
[0030]ΔF
X...

【技术保护点】

【技术特征摘要】
1.一种基于风力机叶片有限元模型的分布力加载方法，其特征在于，包括以下步骤：步骤1、加载基于风力机叶片叶根坐标系，步骤2、将风力机叶片展向进行截面划分，以叶根至叶尖方向将叶片划分为m个截面，叶根处截面标记为1，步骤3、将叶片展向划分的截面n和截面n+1组成的叶片展向段，标记为N，此时基于叶根坐标系，叶片每个剖面的截面载荷为载荷作用点为其中，M表示剖面的弯矩，F表示的剪力，X、Y、Z分别为载荷作用点在叶根坐标系中的坐标值，步骤4、以展向段N为叶片有限元模型研究对象，叶片有限元模型包括SS面壳体、PS面壳体和腹板，展向段N内标记SS面壳体和PS面壳体若干个壳体节点，其中，PS面壳体和SS面壳体主梁区域的节点数目分别记为nump，nums，步骤5、通过截面n和n+1，计算展向段N内所有壳体节点的总节点合力以及载荷作用点坐标偏差，Y
n+1
‑
Y
n
，Z
n+1
‑
Z
n
＞，其中，为载荷作用点坐标偏差，X
n
,Y
n
,Z
n
分别表示n截面的X、Y、Z轴载荷作用点坐标，X
n+1
,Y
n+1
,Z
n+1
分别表示n+1截面的X、Y、Z轴载荷作用点坐标，同时标记展向段N内合力作用点Loc
n
为截面n的载荷作用点位置，步骤6、计算展向段N内PS面壳体和SS面壳体的所有壳体节点的坐标之和，分别标记为：∑NodesX
P NodesY
P NodesZ
P
，∑NodesX
S ∑NodesY
s ∑NodesZ
s
，其中，∑NodesX
P
为PS面壳体的所有壳体节点的X轴坐标值之和，∑NodesY
P
为PS面壳体的所有壳体节点的Y轴坐标值之和，∑NodesZ
P
为PS面壳体的所有壳体节点的Z轴坐标值之和，∑NodesX
S
为SS面壳体的所有壳体节点的X轴坐标值之和，∑NodesY
S
为SS面壳体的所有壳体节点的Y轴坐标值之和，∑NodesZ
S
为SS面壳体的所有壳体节点的Z轴坐标值之和，步骤7、计算展向段N内有限元模型的PS面壳体和SS面壳体的所有壳体节点坐标，分别标记为：和其中，dxp、dyp、dzp为PS面壳体节点X、Y和Z方向的坐标值，dxs、dys、dzs为SS面壳体节点X、Y和Z方向的坐标值，步骤8、将展向段N内所有壳体节点的合力以节点坐标均匀分布的方式，基于力和力矩同时等效的原则，计算展向段N内节点的分布力转换矩阵[ΔC]，其中，ΔM为壳体节点的弯矩合力，ΔF为壳体节点的剪力合力,
步骤9、计算展向段N内PS面壳体和SS面壳体的每个节点的分布力载荷；步骤10、重复步骤5至步骤9的方法，进行下一截面的计算，直至m个截面的每个节点的分布力载荷全部计算完成，根据每个节点的分布力载荷模拟风力机叶片真实风况载荷。2.根据权利要求1所述的一种基于风力机叶片有限元模型的分布力加载方法，其特征在于，所述的步骤1中，风力机叶片叶根坐标系包括3个弯矩M
X
、M
Y
、M
Z
和3个剪力F
X
、F
Y
、F
Z
。3.根据权利要求1所述的一种基于风力机叶片有限元模型的分布力加载方法，其特征在于，所述的步骤5中，展向段N内所有节点合力的计算方法为：其中，M
n
为n截面所有节点的弯矩，M
n+1
为n+1截面所有节点的弯矩，F
n
为n截面所有节点的剪力，F
n+1
为n+1截面所有节点的剪力，载荷坐标的偏差的计算方法为Y
n+1
‑
Y
n
，Z
n+1
‑
Z
n
＞，其中，X
n
,Y
n
,Z
n
分别表示n截面的X、Y、Z轴载荷作用点坐标，X
n+1
,Y
n+1
,Z
n+1<...

【专利技术属性】
技术研发人员：李慧，魏鑫，李丹，周勃，包洪兵，
申请(专利权)人：无锡风电设计研究院有限公司，
类型：发明
国别省市：