一种风机叶片建模方法与流程

文档序号:11251109阅读:3526来源:国知局
一种风机叶片建模方法与流程
本发明涉及参数化的复杂曲面建模,尤其涉及一种风机叶片建模方法。
背景技术
:风机叶片是风力发电机组中的核心部件,其成本占整个风机组成本15%-20%,并且风机叶片采用合理的结构、先进的材料和科学的工艺才能使其能承担风力、自重和离心力等载荷给予的各种弯矩与拉力,这样才能保证整个风电机组的整体性能和发电效率。而建立风机叶片模型并对其进行载荷分析,存在以下问题:(1)叶片结构非常复杂,建模耗费大量的时间和精力,而简化的模型精度一般。(2)对三维模型的有限元模型的建立,修改和网格划分上需要花费大量的时间。(3)ansys中包含的命令更多的偏向于建模和有限元分析,但是在数值计算和矩阵操作方面较为不足。(4)叶片模型的网格化在数值模拟过程中非常耗时,同时网格质量情况决定着数值模拟结果的精确程度甚至决定数值模拟是否收敛,而叶片曲面复杂,每个截面的厚度都不同,所以无论采用自由网格划分还是映射网格划分,都不能得到质量高的网格。因此,本领域的技术人员致力于开发一种风机叶片建模方法,其采用参数化建模方法对复杂曲面的风机叶片进行三维建模和有限元建模和计算,为设计研发阶段提供结构强度校核和力学性能评估,同时本发明提供了解决上述问题的风机叶片建模方法:(1)在建立三维模型时,采用分段建模的方法,然后再合成一体,可以减少建模的复杂性和时间长度。(2)本发明采用hypermesh软件,其是一个高性能的有限元前后处理器,提供高度交互及可视化的环境进行仿真分析工作。与其他的有限元前后处理器比较,hypermesh的图形用户界面易于学习,特别是它支持直接输入已有的三维cad几何模型(ug、pro/e和catia等)已有的有限元模型,并且导入的效率和模型质量都很高,可以大大减少很多重复性的工作。(3)本发明采用matlab和ansys共同进行有限元分析和优化,由于matlab是功能强大的科学计算软件,具有丰富的矩阵操作命令和数值计算函数,可以方便和快速的使用优化算法实现对象设计参数的优化,因此利用matlab软件计算得到实际三维坐标点,并导入ansys中进行参数化建模。(4)本发明采用的专业软件hypermesh。在cae(compueraidedengineering)领域中,hypermesh最突出的特点是具有强大的有限元网格前处理和后处理,并且能直接与cad几何模型和有限元模型对接,以减少建模的重复工作,高质量的自动网格划分,极大地简化复杂几何的有限元建模过程。技术实现要素:有鉴于现有技术的上述缺陷,本发明所要解决的技术问题是提供一种采用参数化建模方法对复杂曲面的风机叶片进行三维建模和有限元建模的风机叶片建模方法。为实现上述目的,本发明提供了一种风机叶片建模方法,包括:将制造数据处理成适用于三维建模软件的叶片几何参数;将所述叶片几何参数导入到所述三维建模软件中以建立叶片三维模型;将所述叶片三维模型进行网格化处理以建立叶片网格模型;将所述叶片网格模型进行铺层设计以建立复合材料叶片有限元参数化模型。进一步地,还包括:在所述复合材料叶片有限元参数化模型上,进行载荷的施加和约束以计算相关力学参数。进一步地,所述制造数据处理包括采用matlab编程方式将离散的数据点处理成所述三维建模软件能够识别的数据格式。进一步地,所述三维建模软件采用catia软件。进一步地,所述将所述叶片几何参数导入到所述三维建模软件中以建立叶片三维模型包括以下步骤:将所述叶片几何参数导入到所述三维建模软件中,并且利用曲面造型功能建立叶片外形轮廓;根据所述叶片外形轮廓将叶片分为至少2个叶片段,并且对每个叶片段分别建立叶片段模型;将所有叶片段模型合并成一个完整的叶片模型。进一步地,所述网格化处理包括采用hypermesh软件对所述叶片的三维模型进行网格化处理。进一步地,所述网格化处理包括以下步骤:根据叶片结构将所述叶片的三维模型分组;对每个组分别进行网格化处理。进一步地,所述叶片结构包括主梁、腹板、前缘和后缘。进一步地,采用apdl语言结合matlab在ansys中建立复合材料叶片有限元参数化模型。进一步地,所述铺层设计包括沿叶片轴向进行分段铺层,直至叶尖位置。进一步地,所述分段铺层包括以下步骤:沿叶片轴向将叶片分段成多个子叶片段;对每一子叶片段,获取所述子叶片段的铺层坐标向量,根据所获取坐标向量对应的坐标位置依次对该叶片段进行多层的铺层。进一步地,所述力学参数包括最大挥舞、最小挥舞、最大摆振和最小摆振。根据本发明提供的风机叶片建模方法,采用编程手段进行自动三维建模和自动有限元建模,采用分段建模思路和分段铺层方法大大提高了建模效率和确保了复合材料铺层建模的准确性,最后计算所得结果与实验测试结果精度可以控制在10%以内。以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明,以充分地了解本发明的目的、特征和效果。附图说明图1是本发明的优选实施例的风机叶片建模方法示意图;图2是截取的叶片分段模型1结构图;图3是截取的叶片分段模型2结构图;图4是截取的叶片分段模型3结构图;图5是截取的叶片分段模型4结构图;图6是叶片整体模型结构图;图7是叶片精确模型结构图;图8是叶片部分精确模型图;图9是叶片部分网格模型图;图10是铺层apdl代码的流程图;图11是叶片剖面结构图;图12是载荷位置示意图;下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。具体实施方式下面结合附图和某型号风力发电机叶片对本发明作进一步说明,采用编写的程序对其进行三维自动建模和有限元铺层建模,最后计算其静力学和模态,与实验测试结果的对比。几何建模工具采用三维建模软件catiap3v5r20;有限元网格划分采用hypermeshv12.0工具;计算软件采用ansys;在生成apdl代码过程中,采用matlab编程工具。第一步:采用matlab制作数据进行处理(制作数据是从已有的叶片获取结构数据并进行进一步优化后的数据),生成catia三维建模软件所识别的数据格式,即使用matlab将离散的数据点处理成可以导入catia的数据格式。第二步:将已处理的数据文件导入三维建模软件catia中,使用几何坐标点自动插值生成曲面的功能建立叶片的外形模型。具体操作如下:(2.1)在catia安装目录下找到\intel_a\code\command\gsdpointsplineloftfromexel.xls文件。(2.2)将叶片几何数据写入gsdpointsplineloftfromexel.xls文件中(2.3)单击菜单“工具”→“宏”,运行fueil1.main后输入3,即可放样生成曲面。其中放样是将一个二维形体对象作为沿某个路径的剖面转换成其相对应的三维模型。例如“宏”的运行方式是多点成线,多线成面,具体是选多个沿叶片长度方向的截面,然后将各截面曲线离散的多个点输入,通过样条曲线生成曲线,然后将个条曲线连接形成叶片外形轮廓,这个叶片外形轮廓是简化的叶片外形,下面需要建立精确的叶片外形三维模型。但是,由于叶片形状不规则,包含的坐标点数据太多,一次性整体建模所需电脑运行时间较长,因此为了建模方便快速,将叶片分为多段(需要分的段数可以根据电脑性能和时间来考量,一般分为几十段),然后分别同时每个叶片段的建模再合并成一个模型的方式来缩短建模时间,具体采用一下步骤:将叶片分成多个叶片段,并且对每个叶片段分别建立叶片段模型。例如,如图2-5所示,从叶片根部到尖部截取出4个有代表性结构的叶片分段模型。将所有叶片段模型合并成一个完成的叶片模型。例如,如图6所示。更加优选的方案,对叶片模型结构包括外形和内部结构(例如主梁1、2、腹板3、4、前缘5和后缘6)进行进一步的构建和划分,以建立叶片精确三维模型,例如,如图7和图8所示。第三步:将叶片的外形模型导入到hypermesh中后,将叶片的各个部分切分出来,主要是将主梁、腹板、前缘以及后缘进行切分,以得到叶片精确的几何外形。例如,如图11所示,将叶片模型按照叶片结构进行分组,即分别将主梁1、2、腹板3、4、前缘5和后缘6划分成各个组,然后对每个叶片结构(即每个组别)进行模态分析前期网格化处理,也是为后面铺层进行预处理。根据叶片形状不规则,曲面复杂,需要在尺寸较大的部分将网格划分的稀一些,以减少计算量,而在尺寸小及形状变化剧烈的部分则将网格划分得密一些,以保证网格划分的质量。因此运用hypermesh自动网格划分功能实现高质量的网格划分,例如,如图9所示的叶片网格模型图,该模型根据叶片形状选取了多个节点,并且根据这些节点划分了多个高精度壳单元。第四步:根据叶片铺层方案图纸,对大梁与腹板布置和叶片主体进行有限元建模叶片铺层,采用matlab编程的方式编写铺层apdl代码,即沿着叶片半径方向分段生成铺层代码,直至叶尖位置。由于风机叶片受力复杂、结构特殊、有精度要求等特点,故采用合理的铺层设计,可使结构设计更为合理,主要包括:(1)叶片各截面的构造;(2)各截面铺层材料的选择和厚度的计算,所以通过优化设计减轻重量以更好地发挥结构的效能。对复合材料叶片铺层进行优化设计,即在满足某种约束条件使各项力学参数达到最优。在铺层材料的选择,例如叶片的最外层为蒙皮,蒙皮主要的功能提供叶片的气动外形,同时承担部分弯曲载荷和大部分剪切载荷。蒙皮的层状结构包括胶衣层、玻纤毡增强层、强度层等,胶衣层提供光滑的气动表面,以提高叶片的气动性能,玻纤毡增强层提供缓冲,强度层为蒙皮的承载层,获取以上复合材料的参数,进行有效计算,达到铺层设计的需求。如图10所示,铺层apdl代码的流程包括:(4.1)沿叶片轴向将叶片分段成多个子叶片段,例如分为子叶片段1、子叶片段2、……、子叶片段n。(4.2)获取子叶片段i的铺层坐标向量,对子叶片段i进行m层的铺层,其中i为循环变量参数,i=1,2,…,n,例如当i的值为2时,获取子叶片段2的铺层坐标向量,然后根据获取的子叶片段2的铺层坐标向量相应的坐标位置调用第1层铺层函数以实现第一层的铺层,调用第2层铺层函数以实现第二的铺层,……,调用第m层铺层函数以实现第二的铺层,然后i自动累加,进入到下一个子叶片段的,重复循环上面铺层的流程,即获取子叶片段3的铺层坐标向量,然后根据获取的子叶片段3的铺层坐标向量相应的坐标位置调用第1层铺层函数以实现第一层的铺层,调用第2层铺层函数以实现第二的铺层,……,调用第m层铺层函数以实现第n的铺层,依次循环,直至子叶片段n的m层铺层完成,达到叶片的尖点位置,整个叶片模型的铺层完成。更好的优选方案,对子叶片段i进行m层铺层时,由于叶片不同结构所使用的铺层材料也不同,所以可以根据叶片结构来划分不同的区域进行分别铺层,例如,如图11所示,将叶片段i划分为主梁1、2、腹板3、4、前缘5和后缘6部分,分别对这些部分进行m层铺层操作。另外,可以利用matlab软件计算得到实际三维坐标点,并导入ansys中进行参数化建模。每隔一个单位长度选取1个截面。由于叶尖主要起降低噪声的作用,对叶片整体结构强度影响较小,在建立有限元模型时去叶尖部分,采用shell99和shell91两种3d壳体单元模拟叶片结构。而matlab和ansys的联合优化:将matlab和apdl语言结合进行二次开发,基于叶片三维形状拟合的方法,将matlab计算所得的叶片空间三维数据点按照特定的格式输出保存,然后用apdl语言编制程序宏文件与matlab建立数据传递机制,建立参数化模型,例如首先通过matlab自带ga函数调用其优化工具箱中遗传算法完成算法优化,然后通过文本文档作为中介进行matlab和ansys的数据传递,最后matlab通过system命令调用ansys批处理模式完成适应度计算。第五步:根据所提供的载荷位置和载荷大小进行约束和载荷的施加。即可计算其力学行为,获取所需的力学参数,其中力学参数包括最大挥舞、最小挥舞、最大摆振和最小摆振等。例如,如图12所示,以不同加载方式对叶片9个位置11、12、13、14、15、16、17、18、19进行施压,以测量最大挥舞、最小挥舞,其测量结果如下:从测量结果看出沿叶片轴向到尖部的最大挥舞、最小挥舞都是逐渐增大,直至尖部最大。载荷位置最大挥舞最小挥舞位置111270.61481.9位置122541.12963.7位置132811.74445.6位置145082.35927.5位置156352.97409.3位置167623.48891.2位置17889410373位置181016511855位置191143513337第六步:计算结果与实验对比研究如下表所示,是计算结果与实测数据对比,其最后模态分析后的结果和实验测试结果之间误差控制在10%以内。以上详细描述了本发明的较佳具体实施例。应当理解,本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此,凡本
技术领域
中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案,皆应在由权利要求书所确定的保护范围内。当前第1页12
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1