遗传算法在水平轴洋流机叶片优化设计中的应用余
在全球气候变暖的背景下,随着能源短缺和经济高速发展,可再生能源发电技术更需进一步的飞跃,不仅要求处理温室效应而且要求保护自然环境以及与自然环境和谐共存。世界各国可持续发展的主要方向转向开发和利用绿色可再生能源。潮流能作为一种对环境几乎不会产生污染的可再生能源,蕴藏量丰富,而且具备相对成熟的利用技术,开发潜力极为巨大[1-2]。
叶片的性能直接影响着洋流机的性能,也是洋流机捕获洋流能的唯一部件,其制造成本也占整个发电机组的百分之二十左右[2]。叶片设计是洋流能利用效率的优劣和洋流机运行安全以及发电机组成本控制的关键因素。目前在国内外洋流机叶片设计大多数是采用了风力机的叶片设计方法[3],其中最为常用方法是与相关的约束模型结合的叶素-动量理论,如Wilson模型[4-5]和Glauert模型[6]。其中,Wilson在Glauert研究的基础上考虑了升阻比和叶尖损失的影响提出的以单叶素截面输出最大功率为设计目标的叶片设计方法应用最为广泛[8-9]。该设计方法速度快,仅考虑了单个叶素的气动性能而忽略了叶素之间的相互影响。
洋流机工作环境和水泵相似,因此需要考虑洋流机叶片的汽蚀问题。汽蚀现象是在流体流动过程中,当某一局部区域的压力不大于水温相对应的汽化压力时,流体就会在该局部区域发生汽化[9-10]。汽化压力被定义为液体发生汽化时的压力,和液体种类、温度有关[11]。汽化产生的气泡顺流进入高压区发生破裂,引发周围液体高频碰撞从而导致材料受到破坏[12-13]。
2018年12月4日,第11届主题为“车联世界 数字运营 出行智能”TC汽车互联网大会在上海隆重拉开帷幕。在本次大会上,华为重磅宣布,车联网平台云服务将在今年12月底上线公测,这是华为在车联网发展中的又一里程碑。
本文基于洋流机叶片全局优化的目的,并且考虑叶轮汽蚀现象,采用遗传算法优化设计60 kW的洋流机叶片,与文献[2]中的Wilson设计叶片进行对比,分析这两种设计叶片弦长和扭角以及额定工况和偏工况下的水动性能,从而验证遗传算法在洋流机叶片设计应用中的有效性。
船舶污染物接收处置能力的组成分为两部分。一是船舶污染物接收处置相关资源的能力(即“硬件能力”),二是将已有资源投入到船舶污染物接收处置活动中使其充分发挥作用的资源运用能力(即“软件能力”)。
1 叶轮汽蚀
在洋流机叶片设计过程中需要考虑汽蚀问题。在叶轮旋转时发生汽蚀,升力系数下降,阻力系数上升,没有足够的升力,输出功率受到直接的影响[12-13]。汽蚀的发生是由翼型最小压力系数Cpmin所决定。翼型最小压力系数Cpmin定义为
由叶素理论推得轴向力F和转矩M的表达式分别为
(1)
式中:p为流体压力,Pa;pLmin为最小局部压力,Pa;ρ为流体密度,kg·m-3;U为来流速度,m·s-1。
汽蚀系数σ定义为
(2)
式中,pvap为汽化压力,Pa。
在优化设计叶片过程中,定义K为避免汽蚀的安全系数。通过设置K小于1,使得叶片形状上具有更高抗气蚀能力。为了避免汽蚀,必须保证方程(3)成立,即
Kσ+Cpmin≥0
(3)
2 叶素-动量设计理论
2.1 动量理论
动量理论引入轴向和径向诱导因子的概念,假设叶轮旋转平面为单元流管致动盘模型,分析致动盘前后压力变化情况,根据伯努利方程和质量守恒定理可得到制动盘所受到载荷与诱导因子之间的关系。动量理论轴流模型如图1所示,图中:U1为尾流速度,m·s-1;p+、p-分别为盘前、盘后压力,Pa;a为轴向诱导因子。
图1 动量理论轴流模型Fig.1 Axial flow model based on momentum theory
由动量理论推得轴向力F和转矩M表达式分别为
dF=4πρa(1-a)U2rdr
许飞卖力地搓洗着我的衣服,他累得满头是汗,扭头问我,陈胖子,你过的是什么日子,你妈的铲子没把你练成个人型么?
(4)
dM=2πρΩ(1-a)Ur3dr
(5)
式中:Ω为水轮转速,r·min-1;r为半径,m。
2.2 叶素理论
叶素理论假设叶片沿展向分为多个片段,每一个片段视为一段叶素,各个叶素相对独立互不影响。考虑叶片旋转效应,结合翼型气动特性分析每段叶素受力并且沿翼展方向积分,求得叶片上所受的力,联立动量理论求解每一段叶素处的诱导因子。叶素速度三角形及受力分析如图2所示,图中:W为相对速度,m·s-1;α、β、φ分别为攻角、扭角和入流角,(°);b为切向诱导因子;FL、FD分别为叶素升力和叶素阻力,N;Fy=FDsin φ+FLcos φ,Fx=FLsin φ-FDcos φ。
图2 叶素速度受力图Fig.2 Velocity and forces of the blade-element
相对速度W为
主动调查取证。强化与公安机关以及生态环境等行政执法部门的联动,建立检察机关提前介入环境损害调查机制,充分发挥各自职能优势,共同开展现场勘查等调查工作。通过询问证人、现场勘查、暗中走访等方式进行明察暗访,听取群众意见,调查环境受损第一手资料。强化科技引领,积极探索运用便携式污染检测设备、无人机航拍,借助卫星遥感等新技术进行调查,获取、保全技术性证据,提高线索比对和取证固证水平。
(6)
[1] Edward Wong, “Chinese Military Seeks to Extend Its Naval Power”, New York Times, April 23, 2010.
(7)
(8)
式中:c为叶素截面弦长,m;N为叶片数;CL为升力系数;CD为升力系数。
2.3 叶素动量理论
(2) 偏工况下,在额定流速以及低于额定流速时优化设计叶片能量利用系数更大,平均提高了4.6%。同时表明,在额定流速以下优化设计叶片可以更好地利用洋流能。
(9)
(10)
3 遗传算法优化模型
在宏观意义上,遗传算法是一种仿生算法,它是模仿生物在自然环境中的进化和遗传过程而形成的一种自适应全局优化概率搜索算法。通过对原始群体进行交叉、复制以及突变等遗传学操作,并以目标函数为适应度,不断更新种群从而得到最优个体[14]。最早由Holland于1975年提出这种不断更新迭代得到最优个体的优化算法[15]。遗传算法因其易实现性和通用性,使得其在工程问题中得到广泛应用[16-23]。
3.1 参数化曲线
叶片设计主要是为了获得沿叶片展向扭角和弦长的分布规律,因此每段叶素的弦长和扭角是遗传算法优化控制变量。为减少控制变量以加快计算速度,通过Bezier曲线定义弦长和扭角的分布规律。在Bezier曲线中,设定控制点为Qj(j=0、1、2、…、n),则n次Bezier曲线Q(t)的方程[23-25]为
(11)
Bj,n(t)=Cj,ntj(1-t)n-j
(12)
(13)
式中,t∈[0,1]。
扭角控制方法与弦长相同。以弦长控制为例,弦长沿展向给定5个控制点,通过Bezier曲线原理得到弦长沿叶片展向分布规律。通过改变控制点纵向位置来控制弦长的变化,以达到改变每段叶素弦长的目的。弦长控制点如表1所示。Bezier曲线如图3所示,其中R为叶轮半径,m。
表1 弦长控制点Tab.1 Control point of the chord
相对半径r/R0.200000.260860.434290.693711.00000弦长c/m0.6420.4510.3320.2230.100
图3 弦长控制点和Bezier曲线Fig.3 Chord control points and the Bezier curve
3.2 优化目标
优化目标为额定流速(1.8 m·s-1)能量利用系数最大。因此,优化目标为max Cp
“意之所之玄之又玄”即是气要追求的境界。而此意是道之意,是远之意,也是琴之意。这就像是投石击水,水波泛起了层层涟漪。琴者的演奏虽然已经停止但是琴韵悠远,琴心已经翻越了千山万山,来到了一片广阔之地。因而是“求之弦中有不足,求之弦外则有余也”。
(14)
式中,Cp为能量利用系数。
3.3 约束条件
在叶片径向位置定义控制点,弦长和扭角变化曲线也就是控制点所定义的Bezier参数化曲线。本文对于弦长和扭角分别定义5个控制点,共10个控制点,以弦长为例(如图3所示)。需要考虑洋流机在水流环境下的汽蚀现象,并对汽蚀安全系数Cf进行约束。控制点和避免汽蚀安全系数约束方程为
(15)
式中:cmax、cmin分别为最大弦长和最小弦长,m;ccpi(i=1、2、3、4、5)为各弦长控制点;tmax、tmin分别为最大扭角和最小扭角,(°);tcpi为各扭角控制点。
4 结果分析
叶片优化设计参数和遗传算法参数如表2所示。
陈山利率二人避开日本人的目光,踏着齐腰深的草丛,抢到日军前面,凭借三棵大树掩护,发现是一股逃兵。待敌人走近,陈山利突然现身,将冲锋枪对准日兵,大喊着:“缴枪不杀!”
表2 主要设计参数及遗传算法参数Tab.2 Main design parameters of the blade and its genetic algorithm configuration parameters
参数名称数值参数名称数值额定功率P/kW60遗传操作种群数100叶轮半径R/m3.5遗传突变比例0.04轮毂半径r1/m0.7遗传交叉比例0.7额定流速Ur/(m·s-1)1.8最大遗传代数150叶片数N3
叶片优化具体过程如图4所示。经过104代遗传优化得到了优化设计叶片,其与Wilson设计叶片的弦长和扭角比较如图5所示。
图4 优化设计遗传算法算流程Fig.4 Flowchart of the optimization design according to genetic algorithm
图5 优化叶片与Wilson设计叶片的弦长和扭角Fig.5 Comparison of chord and twist between the optimized blade and the blade based on Wilson method
如图5(a)所示,优化设计叶片沿径向弦长均小于Wilson设计叶片弦长,平均减小约16.2%。减小了叶片弦长从而很大程度上减少了叶片质量,从而减轻了叶片受到的重力载荷,不仅叶片的运行安全得到了保障,而且节省了材料,从而降低了发电成本。
(1) 优化设计叶片与Wilson设计叶片相比弦长平均减小了16.2%,叶根处扭角减小,很大程度上减轻了叶片的质量,降低了叶片扭曲程度;
学生上课使用手机的内因 统计结果见表5。总体上,学生上课使用手机的三种主要内因所占比例均达到三成以上,但学生上课使用手机更倾向于个人自控能力与缺乏学习动力、抵抗手机诱惑能力弱。
在实际操作中,他却并没有这样做。2016年以来,三都县在建的千万元以上项目有127个,但与脱贫攻坚有关的只有41个。他把精力和资金都集中到与脱贫攻坚工作无直接关系的“养生谷”“千神广场”等“高大上”的综合开发项目上了。都江镇曾是梁嘉庚对口帮扶的乡镇,但他却很少进村,就算去了也只是带着商人看项目,在村委会开个会就走了。
作为当代的小学数学教师需要具备一定的创新能力,在实际的教学中需要保留传统教学中的一些优点,同时根据现在的学生的认知水平将一些先进的教学思想有效地融合到自己的课堂中。通过全新的信息技术手段不断地对自己的教学形式进行创新,通过教学模式的创新,促使学生在学习数学中得到全面发展。比如说在学习人教版五年级数学下册有关长方体和正方体的知识时,因为小学生的年龄特点,他们的抽象思维能力和空间能力还相对较弱,那么这时候我们就可以融入信息技术手段,通过直观的方式展示这些立体图形,方便学生进行理解和学习。让学生的学习效率得到大大的提升。
图6 优化叶片和Wilson设计叶片的能量利用系数随来流速度的变化Fig.6 Relationship between the energy coefficient and flow velocity for the optimized blade and the blade based on Wilson method
5 结 论
考虑叶轮旋转产生汽蚀问题,并以叶素动量理论为基础,采用遗传算法建立叶片优化设计模型,设计了60 kW洋流机叶片与Wilson设计叶片相比较,主要得到以下结论:
两种设计叶片在定速定桨距工况下不同洋流速度下的能量利用系数如图6所示。从图6中可看出,优化设计来流流速小于额定流速偏工况下,叶片能量利用系数提高很多,最大提高了8.6%;其在额定流速度时,叶片能量利用系数由0.452增加到0.461,相对增加了2.0%,平均增加了4.6%。因此,优化设计叶片具有更强的捕获洋流能的能力,可以更好地利用洋流能,达到优化设计目的。
如图5(b)所示,优化设计叶片叶根处扭角小于Wilson设计叶片,减小了叶轮启动时所需力矩。更小的扭角一定程度上减少了叶根弯矩,保障了叶片的运行安全。
叶素-动量理论是由动量理论和叶素理论耦合而成[2-7]。结合动量理论和叶素理论推导的轴向推力和轴向力矩可以得到诱导因子和弦长关系,即式(4)、(5)分别等于式(7)、(8),则可以得到诱导因子和弦长关系为
参考文献:
[1] 马舜.水平轴潮流能转换系统能量转换率及功率控制研究[D].杭州:浙江大学,2011.
[2] 徐学涛.水平轴潮流能发电机械叶片设计及特性分析[D].杭州:浙江大学,2012.
[3] BATTEN W M J,BAHAJ A S,MOLLAND A F,et al.Hydrodynamics of marine current turbines[J].Renewable Energy,2006,31(2):249-256.
“研讨与练习”题是课文教学后,学生全面理解课文的补充途径。教师能从学生的答题情况了解学生的学习情况以及学生之间存在的个体差异,从而反思教学过程;“研讨与练习”题也能提高学生解读课文的能力,集中学生的学习精力、活跃学生的思维。因此,“研讨与练习”是教学环节中必不可少的内容,是提高学生语文学习能力的关键点。根据现代散文“研讨与练习”题的特点,加上语文学科核心素养要求以及“知识与技能、过程与方法、情感态度与价值观”的目标要求,可以采用以下策略实施课后习题教学:
[4] 张玉良,李仁年,杨从新.水平轴风力机的设计与流场特性数值预测[J].兰州理工大学学报,2007,33(2):54-57.
从2014年9月新生开学开始,外籍教师进入大学英语课堂,在14级医学院进行视听说一体化试行。鉴于新生进校英语基础差、参差不齐,尤其是听力、口语几乎是零基础。这一教学改革深受学生们的欢迎,她们期盼外教能提高她们学习英语的兴趣,能听懂并且运用口语进行简单的日常对话、会话等。教师可以通过音、视频等多种媒介将语言、文化、社交等多方面知识展示给学生,内容多为影视片段赏析、音乐欣赏、情景模拟等。学生的兴趣提高了,注意力更加集中,学习效率也得到了提高。同时,在相关内容的赏析学习过后马上进入口语训练阶段,学生经过适当练习,很快能将听力练习中学到的新的表达方式应用于口语交际,能力提高迅速。
[5] 李春,叶舟,高伟,等.现代大型风力机设计原理[M].上海:上海科学技术出版社,2013.
(3) 工程处治方案。项目区山体高陡,多种形式的崩塌频发而冲击能量大,加之工程工期紧张,考虑到线路以路基形式通过,结合项目区坡脚存在较宽的河流阶地等有利条件,决定在坡脚路基部位设置装配式棚洞被动防护方案。即棚洞结构的桩柱及纵梁采用现浇型式,而棚洞顶板采用预制装配的方式,以快速完成工程施工。放弃了工程施工速度较慢、全部现浇和对地基整体强度要求较高的明洞防护方案。棚洞上部根据危岩落石的冲击能量,设置厚约4 m碎石土缓冲层,并在其中设置厚70 cm的EPS缓冲层,以进一步提高棚洞的抗冲击能力。
[6] 张仲柱,王会社,赵晓路,等.水平轴风力机叶片气动性能研究[J].工程热物理学报,2007,28(5):781-783.
[7] 李朝斌,戴韧.水平轴风力机失速模型及对比验证[J].能源研究与信息,2010,26(4):225-229.
[8] 王俭超.水平轴潮流水轮机叶片设计和模型试验研究[D].青岛:中国海洋大学,2011.
[9] 何川,郭立君.泵与风机[M].北京:中国电力出版社,2008.
[10] BRUJAN E A,IKEDA T,MATSUMOTO I Y.On the pressure of cavitation bubbles[J].Experimental Thermal and Fluid Science,2008,32(5):1188-1191.
[11] FRANC J P,MICHEL J M,et al.Fundamentals of cavitation[M].Netherlands:Springer,2005.
[12] BAL simulation on horizontal axis marine current turbines[J].Renewable Energy,2015,80:15-25.
[13] TRUDEAU M G.Structural and hydrodynamic design optimization enhancements with application to marine hydrokinetic turbine blades[D].Pennsylvania:The Pennsylvania State University,2011.
[14] SIVANANDAM S N,DEEPA S N.Introduction to genetic algorithms[M].Berlin:Springer,2008.
[15] 葛继科,邱玉辉,吴春明,等.遗传算法研究综述[J].计算机应用研究,2008,25(10):2911-2916.
[16] 马永杰,云文霞.遗传算法研究进展[J].计算机应用研究,2012,29(4):1201-1206.
[17] 张宇.大型风力机叶片的振动分析与优化设计[D].沈阳:沈阳工业大学,2013.
[18] 罗文,王莉娜.基于小波分解与遗传算法和支持向量机的短期风速预测[J].太阳能学报,2012,33(8):1327-1333.
[19] HOLLAND J H.Adaption in natural and artificial systems[M].Ann Arbor,MI:The University of Michigan Press,1975.
[20] QIN X Y,CHEN L G,SUN F R,et al.Optimization for a steam turbine stage efficiency using a genetic algorithm[J].Applied Thermal Engineering,2003,23(18):2307-2316.
[21] CHEN Y,LI H,JIN K,et al.Wind farm layout optimization using genetic algorithm with different hub height wind turbines[J].Energy Conversion and Management,2013,70:56-65.
[22] CHAQUET J M,CARMONA E J,CORRAL R.Using genetic algorithms to improve the thermodynamic efficiency of gas turbines designed by traditional methods[J].Applied Soft Computing,2012,12(11):3627-3635.
[23] 施法中.计算机辅助几何设计与非均匀有理B样条[M].北京:北京航空航天大学出版社,1994.
[24] 卓文涛,季锃钏,陈二云,等.翼型气动性能与噪声的综合优化设计方法[J].动力工程学报,2012,32(6):481-486.
[25] 马华,刘峰,任春丽.Bezier曲线的计算机实现[J].西安电子科技大学学报:自然科学版,2002,29(4):565-568.