CN113297776B - 一种风力机叶片有限元建模与铺层方法 - Google Patents
一种风力机叶片有限元建模与铺层方法 Download PDFInfo
- Publication number
- CN113297776B CN113297776B CN202110659408.XA CN202110659408A CN113297776B CN 113297776 B CN113297776 B CN 113297776B CN 202110659408 A CN202110659408 A CN 202110659408A CN 113297776 B CN113297776 B CN 113297776B
- Authority
- CN
- China
- Prior art keywords
- blade
- layer
- layering
- finite element
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000000463 material Substances 0.000 claims abstract description 68
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 239000002131 composite material Substances 0.000 claims abstract description 15
- 238000013507 mapping Methods 0.000 claims abstract description 7
- 239000010410 layer Substances 0.000 claims description 78
- 239000000835 fiber Substances 0.000 claims description 28
- 239000011162 core material Substances 0.000 claims description 27
- 240000007182 Ochroma pyramidale Species 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000003014 reinforcing effect Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 7
- 230000010354 integration Effects 0.000 claims description 6
- 239000002356 single layer Substances 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 239000000758 substrate Substances 0.000 claims description 3
- 238000013461 design Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013144 data compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种风力机叶片有限元建模与铺层方法,包括以下步骤:()导出各个翼型截面的数据点;()将数据点导入有限元软件中,连接数据点形成翼型截面线,将各个翼型截面线连接起来,形成单元面并对单元面进行编号,多个单元面互相连接构成叶片初步的三维模型;()计算铺层材料的表观材料特性;()根据叶片的单元面编号进行复合材料铺层,构建完成风力机叶片最终的三维模型;()以单元面为单位对叶片进行映射划分单元网格,进行有限元分析;本发明的风力机叶片有限元建模与铺层方法能够准确计算出截面复合材料的特性,在复杂的三维风力机叶片模型中划分出精细的简单的网格,不仅缩短了数值计算时间,而且还提高了计算精度。
Description
技术领域
本发明涉及一种风力机叶片的设计方法,具体涉及一种风力机叶片有限元建模与铺层方法。
背景技术
目前,国内在风力机叶片的结构设计方面,还没有系统的设计模型和方法,只有针对某一方面的模型,这些模型还无法归纳成一套可靠的系统设计模型。因为风力机叶片截面形状复杂,特别是在翼展方向每个截面的弦长、扭角、翼型均不同,复合纤维铺层机构复杂,其性能随不同的铺层角度、铺层顺序和铺层厚度而变化的,使得叶片设计变得极为复杂和困难。
为了解决上述问题,申请公布号为CN 110298097 A的发明专利申请公布了一种风力发电机组风轮叶片铺层设计方法,该方法如下:计算铺层面积,总结叶片IFF系数与铺层关系式,叶片铺层与重量关系式;进行纤维布铺层优化设计。该方法能够减少铺层层数,减轻叶片质量,简化了设计与校核,但是上述叶片铺层设计方法仍然存在以下问题:
1、设计过程中只划分了腹板,壳体,梁帽,其中材料属性在仿真模拟中数值不够精确。
2、叶片模型网格划分复杂,划分网格精确度不高,叶片的模型的网格化在数值模拟过程中非常耗时,同时网格质量情况决定着数值模拟结果的精确程度,甚至影响数值导致无法收敛。
发明内容
本发明的目的在于克服现有技术的不足,提供一种风力机叶片有限元建模与铺层方法,该方法能够准确计算出每个截面复合材料的刚度与强度,并且在复杂的三维风力机叶片模型中划分出精细整齐的网格,网格简单,不仅缩短了数值计算时间,而且还提高了计算精度。
本发明解决上述技术问题的技术方案是:
一种风力机叶片有限元建模与铺层方法,包括以下步骤:
(1)导出各个翼型截面的数据点;
(2)将数据点导入有限元软件中,连接数据点形成翼型截面线,用线段将各个翼型截面线连接起来,形成单元面并对单元面进行编号,多个单元面互相连接构成叶片初步的三维模型;
(3)计算铺层材料的表观材料特性;
(4)根据叶片的单元面编号进行复合材料铺层,构建完成风力机叶片最终的三维模型;
(5)以单元面为单位对叶片进行映射划分四边形的单元网格,然后进行有限元分析。
优选地,在步骤(1)中,包括以下步骤:
S1-1、通过翼型设计软件插值计算并导出各个翼型截面的二维坐标;
S1-2、通过转换公式,将二维坐标转化为三维坐标;
S1-3、将三维坐标数据保存为TXT文件。
优选地,在S1-2中,翼型坐标变换公式为:
(x1,y1)=(x0,y0)-(X,Y)
x2=x1×c×cos θ-y1×c×sin θ
y2=y1×c×cos θ+x1×c×sin θ
z2=r
式中,(x0,y0)为原始翼型坐标,(X,Y)为变桨轴与弦线的交点坐标,(x2,y2,z2)为坐标变换后的三维翼型坐标,c为各翼型截面的弦长,θ为各翼型截面的扭角,r为翼型截面到叶根的距离;
对于后掠式叶片,要根据积叠线形状进行调整,根据后掠式叶片积叠线公式可知,积叠线设计为二次曲线,计算出各翼型截面沿叶片展向方向的平移量,计算公式为:
x3=x2×cosα+Δx
y3=y2
z3=z2-x2×sinα
式中,Δx为各个翼型截面x方向的平移量,α为各个翼型截面随着积叠线变化的偏移角,a为线性项系数,b为二次项系数,(x3,y3,z3)为后掠叶片坐标变换后的三维翼型坐标。
优选地,在步骤(2)中,包括以下步骤:
S2-1、将保存了翼型数据点信息的TXT文件导入有限元软件中;
S2-2、通过有限元软件命令连接各截面翼型点成为翼型截面线;
S2-3、划分翼型截面为六个区域,分别是前缘、前缘夹芯区、梁帽、后缘夹芯区、后缘增强区以及腹板;
S2-4、通过有限元软件命令用线段连接各翼型截面线形成单元面,多个单元面互相连接构成叶片初步的三维模型,完成风力机叶片初步的三维建模;单元面的编号顺序分别按照前缘、前缘夹芯区、梁帽、后缘夹芯区、后缘增强区以及腹板排列,按照单元面编号区分铺层时的各个区域。
优选地,在步骤(3)中,包括以下步骤:
S3-1、设定铺层材料,铺层材料分别有UNIAX、BIX、TRIAX以及BALSA,其中,UNIAX为单轴向材料,BIX为双轴向材料,TRIAX为三轴向材料,BALSA为芯材料;叶片截面的前缘、前缘夹芯区、梁帽,后缘夹芯区以及后缘增强区是由TRIAX、UNIAX以及BALSA材料按照不同方向堆叠而成,而叶片的腹板是由BIX、BALSA材料堆叠而成;
S3-2、结合经典的合层板理论使用微力学方程,得出薄片UNIAX、BIX、TRIAX以及BALSA材料的表观材料特性,其中,
单向薄板在纤维方向上的杨氏模量计算公式为:
E1=Ef1Vf+EmVm
式中,Ef1为纤维的纵向杨氏模量,Em为为基体的杨氏模量,Vf为纤维的体积分数,Vm为基体的体积分数;
单向薄板在纤维横向上的表观杨氏模量计算如下:
式中,Ef2为纤维的横向杨氏模量;
单向薄板的表观面内剪切模量计算公式为:
式中,Gf12为纤维的面内剪切模量,Gm为基体的剪切模量;
单向薄板的表观面外剪切模量计算公式为:
式中,Gf23为纤维的面内剪切模量;
单向薄板的表观面内泊松比计算公式为:
v12=vf12Vf+vmVm
式中,vf12为纤维的泊松比,vm为基体的泊松比;
单向薄板的表观面外泊松比计算公式为:
单向薄板的表观质量密度计算公式为:
ρ=ρfVf+pmVm
式中,ρf是纤维的质量密度,ρm是矩阵的质量密度;
最后,单向薄板的剩余弹性性质由对称性质决定:
E3=E2
G13=G12
v13=ν12
根据上述公式可以分别计算出薄片UNIAX、BIX、TRIAX与BALSA等材料的表观材料特性。
优选地,在步骤(4)中,包括以下步骤:
S4-1、采用SHELL181单元模拟复合材料叶片的铺层结构;
S4-2、根据前缘、前缘夹芯区、梁帽、后缘夹芯区、后缘增强区以及腹板六个不同区域的铺层材料厚度,层数与角度,使用有限元软件命令语言编写循环代码,以弦线为标准将翼型分为上下两部分进行对称铺层。
优选地,步骤S4-2中,包括以下步骤:
S4-2-1、读取风力机叶片三维模型数据;
S4-2-2、进入前处理模式,建立m个存储复合材料层数的1×n数组,其中n为叶片的总面数,由于材料对称堆叠,因此第一层层数等于第m层,第二层层数等于第m-1层,第三层层数等于第m-2层;由于每一层的单层厚度与总厚度是确定的,因此可以计算出每一层材料的层数;
S4-2-3、读取记录m个层数TXT格式的数据文件,将读取的数据存储在步骤S4-2-2创建的数组中,建立壳单元shell181,单元编号为1;
S4-2-4、写入UNIAX、BIX、TRIAX与BALSA材料的三个方向的弹性模量、剪切模量、泊松比以及质量密度,根据材料输入顺序给材料编号为1、2、3、4:
S4-2-5、定义截面属性,铺层循环开始,先进行外循环,初始化循环参数i=1,步长为1,i==n时循环结束,其中n为叶片的总面数;接着进行内循环,首先,对第一层进行铺层循环,初始化循环参数j=1,步长为1,j=m时循环结束,其中m为第一层层数,循环输入单层厚度,材料编号,铺层角度,积分点数量;然后,对第二层进行铺层循环,初始化循环参数k=1,步长为1,k=o时循环结束,其中o为第二层层数,继续循环输入单层厚度,铺层角度,材料编号,积分点数量;接着,依次类推,对剩下的层数进行铺层循环,由于铺层材料是对称铺层,因此实际上循环代码第一层与第m层相同,第二层与第m-1层相同,第三层与第m-2层相同。
优选地,在步骤(5)中,首先,赋予翼型截面的属性,选择单元编号为1,接着,以单元面为单位通过映射划分四边形的单元网格,最后,通过调整划分的网格数量来控制单元数量,以便进行有限元计算分析。
本发明与现有技术相比具有以下的有益效果:
1、本发明的风力机叶片有限元建模与铺层方法通过合层板原理可以准确的计算出复合材料的刚度和强度。
2、本发明的风力机叶片有限元建模与铺层方法通过划分单元面,而且划分的是四边形网格,避免了大小不均匀的三角形网格;能够在复杂的叶片表面上轻松采用扫掠划分出均匀的单元,在减少网格数量的同时,获得较高的精度,缩短计算时间。
3、本发明直接利用有限元软件内置的命令语言进行建模、四边形面划分网格以及铺层,减少了数据传输时因压缩数据导致的信息缺失问题;按照每个面为单位进行铺层,不仅精度高,而且容易控制角度以及厚度。
附图说明
图1为风力机叶片三维模型连接线段的示意图。
图2为风力机叶片翼型截面区域分布示意图。
图3为有限元软件命令语言铺层代码流程图。
图4为风力机叶片网格划分的示意图。
图5为本发明的风力机叶片有限元建模与铺层方法的工作流程图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
参见图1-图5,本发明的一种风力机叶片有限元建模与铺层方法,包括以下步骤:
(1)导出各个翼型截面的数据点;
(2)将数据点导入有限元软件中,连接数据点形成翼型截面线,用线段将各个翼型截面线连接起来,形成单元面并对单元面进行编号,多个单元面互相连接构成叶片初步的三维模型;
(3)计算铺层材料的表观材料特性;
(4)根据叶片的单元面编号进行复合材料铺层,构建完成风力机叶片最终的三维模型;
(5)以单元面为单位对叶片进行映射划分四边形的单元网格,然后进行有限元分析。
本实施例中,步骤(1)包括以下步骤:
S1-1、通过翼型设计软件插值计算并导出各个翼型截面的二维坐标;
S1-2、通过转换公式,将二维坐标转化为三维坐标;
S1-3、将三维坐标数据保存为TXT文件。
本实施例中,在S1-2中,由于翼型设计软件导出的各个翼型二维数据的弦长都默认设置为1,因此在建立三维模型前要将每一个翼型截面都按各自弦长等比例放大;由于翼型截面的扭角随着叶片展向方向旋转,因此可以将翼型用变换公式数据进行处理,其中,翼型坐标变换公式为:
(x1,y1)=(xo,y0)-(X,Y)
x2=x1×c×cos θ-y1×c×sin θ
y2=y1×c×cos θ+x1×c×sin θ
z2=r
式中,(x0,y0)为原始翼型坐标,(X,Y)为变桨轴与弦线的交点坐标,(x2,y2,z2)为坐标变换后的三维翼型坐标,c为各翼型截面的弦长,θ为各翼型截面的扭角,r为翼型截面到叶根的距离;
对于后掠式叶片,要根据积叠线形状进行调整,根据后掠式叶片积叠线公式可知,积叠线设计为二次曲线,计算出各翼型截面沿叶片展向方向的平移量,计算公式为:
x3=x2×cosα+Δx
y3=y2
z3=z2-x2×sin α
式中,Δx为各个翼型截面x方向的平移量,α为各个翼型截面随着积叠线变化的偏移角,a为线性项系数,b为二次项系数,(x3,y3,z3)为后掠叶片坐标变换后的三维翼型坐标。
本实施例中,参见图2,步骤(2)包括以下步骤:
S2-1、将保存了翼型数据点信息的TXT文件导入有限元软件中;
S2-2、通过有限元软件命令连接各截面翼型点成为翼型截面线;
S2-3、划分翼型截面为六个区域,分别是前缘1、前缘夹芯区2、梁帽3、后缘夹芯区4、后缘增强区5以及腹板6;
S2-4、通过有限元软件命令用线段连接各翼型截面线形成单元面,多个单元面互相连接构成叶片初步的三维模型,完成风力机叶片初步的三维建模;单元面的编号顺序分别按照前缘1、前缘夹芯区2、梁帽3、后缘夹芯区4、后缘增强区5以及腹板6排列,按照单元面编号区分铺层时的各个区域。
本实施例中,参见图2,步骤(3)包括以下步骤:
S3-1、设定铺层材料,铺层材料分别有UNIAX、BIX、TRIAX以及BALSA,其中,UNIAX为单轴向材料,BIX为双轴向材料,TRIAX为三轴向材料,BALSA为芯材料;叶片截面的前缘1、前缘夹芯区2、梁帽3,后缘夹芯区4以及后缘增强区5是由TRIAX、UNIAX以及BALSA材料按照不同方向堆叠而成,而叶片的腹板6是由BIX、BALSA材料堆叠而成;
S3-2、结合经典的合层板理论使用微力学方程,得出薄片UNIAX、BIX、TRIAX以及BALSA材料的表观材料特性,其中,
单向薄板在纤维方向上的杨氏模量计算公式为:
E1=Ef1Vf+EmVm
式中,Ef1为纤维的纵向杨氏模量,Em为为基体的杨氏模量,Vf为纤维的体积分数,Vm为基体的体积分数;
单向薄板在纤维横向上的表观杨氏模量计算如下:
式中,Ef2为纤维的横向杨氏模量;
单向薄板的表观面内剪切模量计算公式为:
式中,Gf12为纤维的面内剪切模量,Gm为基体的剪切模量;
单向薄板的表观面外剪切模量计算公式为:
式中,Gf23为纤维的面内剪切模量;
单向薄板的表观面内泊松比计算公式为:
v12=vf12Vf+vmVm
式中,vf12为纤维的泊松比,vm为基体的泊松比;
单向薄板的表观面外泊松比计算公式为:
单向薄板的表观质量密度计算公式为:
ρ=ρfVf+ρmVm
式中,ρf是纤维的质量密度,ρm是矩阵的质量密度;
最后,单向薄板的剩余弹性性质由对称性质决定:
E3=E2
G13=G12
v13=v12
根据上述公式可以分别计算出薄片UNIAX、BIX、TRIAX与BALSA等材料的表观材料特性。
本实施例中,步骤(4)包括以下步骤:
S4-1、采用SHELL181单元模拟复合材料叶片的铺层结构;
S4-2、根据前缘1、前缘夹芯区2、梁帽3、后缘夹芯区4以及后缘增强区5五个不同区域的铺层材料厚度,层数与角度,使用有限元软件命令语言编写循环代码,以弦线为标准将翼型分为上下两部分进行对称铺层。
本实施例中,步骤S4-2包括以下步骤:
S4-2-1、读取风力机叶片三维模型数据;
S4-2-2、进入前处理模式,建立m个存储复合材料层数的1×n数组,其中n为叶片的总面数,由于材料对称堆叠,因此第一层层数等于第m层,第二层层数等于第m-1层,第三层层数等于第m-2层;由于每一层的单层厚度与总厚度是确定的,因此可以计算出每一层材料的层数;
S4-2-3、读取记录m个层数TXT格式的数据文件,将读取的数据存储在步骤S4-2-2创建的数组中,建立壳单元shell181,单元编号为1;
S4-2-4、写入UNIAX、BIX、TRIAX与BALSA材料的三个方向的弹性模量、剪切模量、泊松比以及质量密度,根据材料输入顺序给材料编号为1、2、3、4;
S4-2-5、定义截面属性,铺层循环开始,先进行外循环,初始化循环参数i=l,步长为1,i=n时循环结束,其中n为叶片的总面数;接着进行内循环,首先,对第一层进行铺层循环,初始化循环参数j=1,步长为1,j=m时循环结束,其中m为第一层层数,循环输入单层厚度,材料编号,铺层角度,积分点数量;然后,对第二层进行铺层循环,初始化循环参数k=1,步长为1,k=o时循环结束,其中o为第二层层数,继续循环输入单层厚度,铺层角度,材料编号,积分点数量;接着,依次类推,对剩下的层数进行铺层循环,由于铺层材料是对称铺层,因此实际上循环代码第一层与第m层相同,第二层与第m-1层相同,第三层与第m-2层相同。
本实施例的步骤(5)包括以下步骤:首先,赋予翼型截面的属性,选择单元编号为1,接着,以单元面为单位通过映射划分四边形的单元网格,最后,通过调整划分的网格数量来控制单元数量,以便进行有限元计算分析。
上述为本发明较佳的实施方式,但本发明的实施方式并不受上述内容的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (6)
1.一种风力机叶片有限元建模与铺层方法,其特征在于,包括以下步骤:
(1)导出各个翼型截面的数据点,包括以下步骤:
S1-1、通过翼型设计软件插值计算并导出各个翼型截面的二维坐标;
S1-2、通过转换公式,将二维坐标转化为三维坐标;
其中,坐标转换公式为:
(x1,y1)=(x0,y0)-(X,Y)
x2=x1×c×cosθ-y1×c×sinθ
y2=y1×c×cosθ+x1×c×sinθ
z2=r
式中,(x0,y0)为原始翼型坐标,(X,Y)为变桨轴与弦线的交点坐标,(x2,y2,z2)为坐标变换后的三维翼型坐标,c为各翼型截面的弦长,θ为各翼型截面的扭角,r为翼型截面到叶根的距离;
对于后掠式叶片,要根据积叠线形状进行调整,根据后掠式叶片积叠线公式可知,积叠线设计为二次曲线,计算出各翼型截面沿叶片展向方向的平移量,计算公式为:
x3=x2×cosα+Δx
y3=y2
z3=z2-x2×sinα
式中,Δx为各个翼型截面x方向的平移量,α为各个翼型截面随着积叠线变化的偏移角,a为线性项系数,b为二次项系数,(x3,y3,z3)为后掠叶片坐标变换后的三维翼型坐标;
S1-3、将三维坐标数据保存为TXT文件;
(2)将数据点导入有限元软件中,连接数据点形成翼型截面线,用线段将各个翼型截面线连接起来,形成单元面并对单元面进行编号,多个单元面互相连接构成叶片初步的三维模型;
(3)计算铺层材料的表观材料特性;
(4)根据叶片的单元面编号进行复合材料铺层,构建完成风力机叶片最终的三维模型;
(5)以单元面为单位对叶片进行映射划分四边形的单元网格,然后进行有限元分析。
2.根据权利要求1所述的一种风力机叶片有限元建模与铺层方法,其特征在于,在步骤(2)中,包括以下步骤:
S2-1、将保存了翼型数据点信息的TXT文件导入有限元软件中;
S2-2、通过有限元软件命令连接各截面翼型点成为翼型截面线;
S2-3、划分翼型截面为六个区域,分别是前缘、前缘夹芯区、梁帽、后缘夹芯区、后缘增强区以及腹板;
S2-4、通过有限元软件命令用线段连接各翼型截面线形成单元面,多个单元面互相连接构成叶片初步的三维模型,完成风力机叶片初步的三维建模;单元面的编号顺序分别按照前缘、前缘夹芯区、梁帽、后缘夹芯区、后缘增强区以及腹板排列,按照单元面编号区分铺层时的各个区域。
3.根据权利要求2所述的一种风力机叶片有限元建模与铺层方法,其特征在于,在步骤(3)中,包括以下步骤:
S3-1、设定铺层材料,铺层材料分别有UNIAX、BIX、TRIAX以及BALSA,其中,UNIAX为单轴向材料,BIX为双轴向材料,TRIAX为三轴向材料,BALSA为芯材料;叶片截面的前缘、前缘夹芯区、梁帽,后缘夹芯区以及后缘增强区是由TRIAX、UNIAX以及BALSA材料按照不同方向堆叠而成,而叶片的腹板是由BIX、BALSA材料堆叠而成;
S3-2、结合经典的合层板理论使用微力学方程,得出薄片UNIAX、BIX、TRIAX以及BALSA材料的表观材料特性,其中,
单向薄板在纤维方向上的杨氏模量计算公式为:
E1=Ef1Vf+EmVm
式中,Ef1为纤维的纵向杨氏模量,Em为基体的杨氏模量,Vf为纤维的体积分数,Vm为基体的体积分数;
单向薄板在纤维横向上的表观杨氏模量计算如下:
式中,Ef2为纤维的横向杨氏模量;
单向薄板的表观面内剪切模量计算公式为:
式中,Gf12为纤维的面内剪切模量,Gm为基体的剪切模量;
单向薄板的表观面外剪切模量计算公式为:
式中,Gf23为纤维的面内剪切模量;
单向薄板的表观面内泊松比计算公式为:
v12=vf12Vf+vmVm
式中,vf12为纤维的泊松比,vm为基体的泊松比;
单向薄板的表观面外泊松比计算公式为:
单向薄板的表观质量密度计算公式为:
ρ=ρfVf+ρmVm
式中,ρf是纤维的质量密度,ρm是矩阵的质量密度;
最后,单向薄板的剩余弹性性质由对称性质决定:
E3=E2
G13=G12
v13=v12
根据上述公式可以分别计算出薄片UNIAX、BIX、TRIAX与BALSA等材料的表观材料特性。
4.根据权利要求3所述的一种风力机叶片有限元建模与铺层方法,其特征在于,在步骤(4)中,包括以下步骤:
S4-1、采用SHELL181单元模拟复合材料叶片的铺层结构;
S4-2、根据前缘、前缘夹芯区、梁帽、后缘夹芯区、后缘增强区与腹板六个不同区域的铺层材料厚度,层数与角度,使用有限元软件命令语言编写循环代码,以弦线为标准将翼型分为上下两部分进行对称铺层。
5.根据权利要求4所述的一种风力机叶片有限元建模与铺层方法,其特征在于,步骤S4-2中,包括以下步骤:
S4-2-1、读取风力机叶片三维模型数据;
S4-2-2、进入前处理模式,建立m个存储复合材料层数的1×n数组,其中n为叶片的总面数,由于材料对称堆叠,因此第一层层数等于第m层,第二层层数等于第m-1层,第三层层数等于第m-2层;由于每一层的单层厚度与总厚度是确定的,因此可以计算出每一层材料的层数;
S4-2-3、读取记录m个层数TXT格式的数据文件,将读取的数据存储在步骤S4-2-2创建的数组中,建立壳单元shell181,单元编号为1;
S4-2-4、写入UNIAX、BIX、TRIAX与BALSA材料的三个方向的弹性模量、剪切模量、泊松比以及质量密度,根据材料输入顺序给材料编号为1、2、3、4;
S4-2-5、定义截面属性,铺层循环开始,先进行外循环,初始化循环参数i=1,步长为1,i=n时循环结束,其中n为叶片的总面数;接着进行内循环,首先,对第一层进行铺层循环,初始化循环参数j=1,步长为1,j=m时循环结束,其中m为第一层层数,循环输入单层厚度、材料编号、铺层角度、积分点数量;然后,对第二层进行铺层循环,初始化循环参数k=1,步长为1,k=o时循环结束,其中o为第二层层数,继续循环输入单层厚度、铺层角度、材料编号、积分点数量;接着,依次类推,对剩下的层数进行铺层循环,由于铺层材料是对称铺层,因此实际上循环代码第一层与第m层相同,第二层与第m-1层相同,第三层与第m-2层相同。
6.根据权利要求5所述的一种风力机叶片有限元建模与铺层方法,其特征在于,在步骤(5)中,包括以下步骤:首先,赋予翼型截面的属性,选择单元编号为1,接着,以单元面为单位通过映射划分单元网格,最后,通过调整划分的网格数量来控制单元数量,以便进行有限元计算分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110659408.XA CN113297776B (zh) | 2021-06-11 | 2021-06-11 | 一种风力机叶片有限元建模与铺层方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110659408.XA CN113297776B (zh) | 2021-06-11 | 2021-06-11 | 一种风力机叶片有限元建模与铺层方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113297776A CN113297776A (zh) | 2021-08-24 |
CN113297776B true CN113297776B (zh) | 2022-09-23 |
Family
ID=77328152
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110659408.XA Expired - Fee Related CN113297776B (zh) | 2021-06-11 | 2021-06-11 | 一种风力机叶片有限元建模与铺层方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113297776B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117352111B (zh) * | 2023-12-06 | 2024-03-08 | 城资泰诺(山东)新材料科技有限公司 | 一种复合材料铺层设计优化方法及系统 |
CN117747033B (zh) * | 2024-02-08 | 2024-04-19 | 北京理工大学 | 复合材料网格结构数字化建模方法和建模装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107169191A (zh) * | 2017-05-10 | 2017-09-15 | 上海电气集团股份有限公司 | 一种风机叶片建模方法 |
CN109583023A (zh) * | 2018-10-29 | 2019-04-05 | 兰州理工大学 | 一种风力机叶片尾缘建模方法 |
CN109977526A (zh) * | 2019-03-21 | 2019-07-05 | 北京航空航天大学 | 一种基于三维cst技术的调整机翼有限元模型的方法 |
-
2021
- 2021-06-11 CN CN202110659408.XA patent/CN113297776B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107169191A (zh) * | 2017-05-10 | 2017-09-15 | 上海电气集团股份有限公司 | 一种风机叶片建模方法 |
CN109583023A (zh) * | 2018-10-29 | 2019-04-05 | 兰州理工大学 | 一种风力机叶片尾缘建模方法 |
CN109977526A (zh) * | 2019-03-21 | 2019-07-05 | 北京航空航天大学 | 一种基于三维cst技术的调整机翼有限元模型的方法 |
Non-Patent Citations (5)
Title |
---|
Long Chen等.Optimization design and simulationanalysis of composite material anti-/deicing component for wind turbine blade.《IOP Conference Series: Earth and Environmental Science》.2019, * |
Peter J. Schubel等.Wind Turbine Blade Design.《Energies 2012》.2012, * |
孟瑞等.基于竞争博弈决策方法的风力机叶片气动与结构优化设计.《安徽工业大学学报(自然科学版)》.2017,(第02期), * |
陈文朴等.对比分析不同材料及铺层对叶片力学性能的影响.《能源工程》.2015,(第05期), * |
黄辉秀.一种快速生成风电叶片有限元网格的方法.《玻璃钢/复合材料》.2016,(第01期), * |
Also Published As
Publication number | Publication date |
---|---|
CN113297776A (zh) | 2021-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113297776B (zh) | 一种风力机叶片有限元建模与铺层方法 | |
CN107067471B (zh) | 一种提高悬垂体模型成型质量的自适应扫描速度方法 | |
WO2022037305A1 (zh) | 一种陶瓷基复合材料铺层预制体优化设计方法 | |
CN111460710B (zh) | 一种基于铺丝轨迹的复材固化变形仿真建模方法 | |
CN107391891A (zh) | 一种基于模型融合方法的大展弦比机翼优化设计方法 | |
CN105046021B (zh) | 非定常气动力最小状态有理近似的非线性优化算法 | |
CN112632835B (zh) | 大长细比豆荚杆的临界屈曲载荷分析模型建模方法 | |
CN115238419A (zh) | 一种基于复合材料的轴流压气机叶片设计方法 | |
CN112507587A (zh) | 一种面向压缩稳定性的变刚度复合材料结构优化设计方法 | |
WO2022193477A1 (zh) | 大拉伸量下具有可调节拉胀特性的剪纸超材料及其设计方法 | |
Pei et al. | Path planning based on ply orientation information for automatic fiber placement on mesh surface | |
CN112149322A (zh) | 一种复合材料加筋结构固化变形的有限元分析方法 | |
CN108363835B (zh) | 一种基于t样条曲面的多孔夹芯结构及应用 | |
CN109583023A (zh) | 一种风力机叶片尾缘建模方法 | |
CN110188468B (zh) | 曲线纤维复合材料翼面结构气动弹性剪裁优化方法及系统 | |
CN110210185B (zh) | 一种高超声速隔离段优化设计方法及系统 | |
Maung et al. | Automated manufacture of optimised shape-adaptive composite hydrofoils with curvilinear fibre paths for improved bend-twist performance | |
CN112528537A (zh) | 一种面向压缩稳定性的变刚度复合材料结构分析方法 | |
Wu et al. | Design and optimization of the variable-density lattice structure based on load paths | |
CN108363828A (zh) | 一种变刚度复合材料的建模方法 | |
CN113496060A (zh) | 一种复合材料叶片气动与结构一体化设计方法 | |
CN116305589A (zh) | 一种直升机桨叶结构降阶分析方法、系统、设备及介质 | |
CN110696367A (zh) | 一种壳类零件变厚度曲面分层、路径规划及工艺控制方法 | |
CN116306088A (zh) | 基于共形参数化的多尺度扫掠曲面结构设计方法及设备 | |
Dillinger et al. | Design, manufacturing and identification of an actively controlled flexible wing for subsonic wind tunnel testing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220923 |