基于涡尾迹方法的风力机载荷与响应计算

采用涡尾迹方法计算准定常气动性能,结合Leishman-Beddoes动态失速模型计算叶片的非定常气动力。应用有限元分析软件,分析风力机在不同转速下的风轮模态,塔架的前后、左右振动模态。基于非定常气动力和结构振动模态,针对风力机风轮,建立其结构运动方程。求解出的叶片振动速度加入到非定常气动力的计算中,从而建立了考虑叶片振动的风力机载荷、响应计算模型。分析了Phase VI叶片和某1.5MW风力机各个叶片截面的气动力、振动位移和振动速度随时间的响应曲线,表明该方法能够较好地计算出风力机的气动性能、载荷和响应。     

倾转旋翼过渡状态瞬态响应分析与试验

基于多体动力学分析方法建立倾转旋翼过渡状态瞬态响应分析模型,研究过渡状态下倾转旋翼非线性非定常气弹耦合动力学特性;通过引入倾转过程旋翼尾迹弯曲的影响,修正了直升机旋翼常规动态入流模型。集成非定常动态入流方程与倾转过渡状态的多体动力学方程,建立倾转旋翼过渡状态下时域非定常耦合分析模型。以两片桨叶的跷跷板旋翼为例,分析倾转过渡状态旋翼瞬态挥舞响应及旋翼气动力的时间响应历程。利用旋臂式模型机动飞行试验机进行倾转过渡状态下旋翼瞬态拉力的试验研究,试验与理论计算结果表现出很好的相关性。数值计算和试验结果表明:建立的时域模型能够有效分析倾转旋翼过渡状态的瞬态特性;倾转过渡状态旋翼尾迹弯曲对非定常动态入流的影响是显著的,对于倾转过渡状态动态入流尾迹弯曲的修正是很必要的。      

风力机叶片旋转平面内弯曲对动力学响应的影响

采用变分渐近梁截面分析方法(VABS)精确计算旋转平面内弯曲的复合材料叶片截面特性,耦合空间一维Euler-Bernoulli梁模型,分析了不同弯曲幅度对考虑弯扭耦合特性的风力机叶片动力学性能的影响。并采用非定常叶素动量理论,考虑动态入流和动态失速模型求解叶片的气动载荷,分析了叶片前掠和后掠在相同控制策略时各风速下对气动功率的影响,以及在极端阵风条件下的气动弹性响应情况。结果表明叶片弯曲对截面特性和前几阶的固有频率的影响较小,但在扭转方向具有较大影响。同时在相同控制策略下,叶片后掠会使风力机的功率降低,但能有效降低极端阵风下的载荷波动。     

风力机柔性叶片振动变形对其气动阻尼的影响分析

为了研究大型水平轴风力机柔性叶片的振动变形对其气动阻尼的影响,用叶素动量理论计算了风力机叶片气动力,用有限元法计算了结构动力响应,基于能量损失法建立了风力机的气动阻尼分析模型.重点分析了叶片振动扭角、挥舞倾角、摆振倾角对攻角的影响,考虑了挥舞速度,摆振速度对入流速度、入流角的影响.以某5MW海上风力机为例,分析计算了风力机整个叶片的气动阻尼.结果表明挥舞倾角、扭角对挥舞方向气动阻尼影响较大,会使其减小;摆振倾角对摆振方向气动阻尼的影响较大,会使其增大;振动变形对气动阻尼沿叶片的分布没有影响.     

大型风力机气动弹性响应计算研究

发展了一套快速预测风力机动态响应的分析方法.针对南京航空航天大学设计的兆瓦级大型风力机叶片NH1500,按照工程梁方法对叶片及塔架结构进行简化,生成有限元梁模型.根据动量叶素理论编写气动力程序,为有限元模型提供气动载荷,从而完成时域情况下风力机的响应计算.最后,分析了叶片的几何非线性效应对响应的影响规律.     

基于非线性气动力的飞机飞行载荷计算

为了研究气动力非线性对飞机静气动弹性特性及飞行载荷的影响,发展了一种基于非线性气动力的飞机静气弹及飞行载荷计算方法。通过引入CFD计算结果,实现对线性方法中的气动力影响系数矩阵的非线性修正,采用模态法求解静气弹配平方程得到飞机非线性的气弹变形及飞行载荷。以某翼身组合体构型为计算算例,分别针对刚性模型和柔性模型,基于非线性气动力完成了飞机静气动弹性分析及载荷计算,并与线性结果进行了对比分析。结果表明,非线性方法可对大迎角范围内飞机静气弹变形及飞行载荷做出较合理的预测。     

柔性风轮的动态入流效应研究

基于BEM修正模型和Pitt-Peter模型,推导出非定常流下的动态入流模型.利用四阶Runge-Kutta法求解动态入流模型中的一阶微分方程.考虑大型机组柔性叶片工作过程中发生大尺度变形的特点,加入柔性结构变形对动态入流模型的反馈.分别研究风力机在风剪、风湍流、偏航、叶片桨距角和风轮转速变化过程中的诱导因子的渐变机理.运用MATLAB进行编程,对某5MW风力机进行动态诱导因子计算,建立起诱导因子的动态描述,得到适合柔性叶片翼型的动态入流模型及分析方法.通过对比几种不同模型的计算结果,验证了加柔性结构反馈的动态入流模型的正确性和实用性.     

电控旋翼气弹动力学建模研究

研究了电控旋翼的气弹动力学建模方法.以有限元旋翼气弹分析程序LORA01为基础,首先给出了一种能计及襟翼移轴补偿影响的时域综合非定常气动力模型,并将伺服襟翼作为额外引入的质量体计入结构动力学模型;之后结合有限元法,根据Hamilton原理导出由广义力表示的带襟翼桨叶非线性运动方程,采用基于Newmark方法的隐式数值积分法对桨叶运动方程进行求解.最后,利用主动控制襟翼旋翼的试验数据对模型进行了比较验证,证明了所建模型的正确性.     

新型倾转旋翼机气动弹性稳定性分析模型

针对新型半铰接式、摆振柔软的倾转旋翼机,建立了旋翼/动力舱/机翼气动弹性稳定性分析模型,并导出了系统的动力学方程。分析模型采用了当量铰、弹性约束的刚性桨叶和弹性机翼假设,采用Pitt-Peters动力入流模型计入了非定常气动力的影响。所建模型为进行新型倾转旋翼机大速度飞行时气弹稳定性分析奠定了基础。     

风扇/压气机叶型厚度对颤振特性的影响

采用计算流体力学与结构动力学相结合的方法,数值模拟了三维振荡叶栅非定常黏性流场;通过对叶片表面非定常气动力及其所做非定常气动功的计算分析,采用能量法对叶片颤振与否进行预估判断。针对三维直叶片和风扇转子叶片,通过调整叶片的相对厚度,研究了叶片厚度变化对风扇/压气机颤振特性的影响。此外,通过对振荡叶栅非定常流场结构的研究,发现了叶片吸力面气流分离与叶片振动之间的耦合关系。本文的研究在风扇/压气机设计中,可用于评估最大相对厚度等叶型结构设计参数对气弹稳定性的影响,对叶轮机颤振机理研究具有一定的参考意义。     

涡轮叶片模态实测阻尼比的有限元应用

为了更加精确地通过有限元来模拟涡轮叶片的谐响应特性,剖析了有限元分析中阻尼矩阵的构建方法,并且选定了适用于真实涡轮叶片设计的阻尼矩阵模型.在进行涡轮叶片有限元谐响应计算时,以模态测试得的阻尼比为基本参数,建立构建阻尼矩阵的方法,该方法能处理具有单独共振频率或具有临近共振频率的振动问题.通过模拟涡轮叶片的响应试验,测得了试验系统的阻尼比及模拟涡轮叶片的位移响应.根据测得的系统阻尼比,运用构建阻尼矩阵的方法,对具有单独共振频率特征的模拟涡轮叶片进行了谐响应计算,结果显示模拟涡轮叶片的位移响应与试验结果基本一致.运用该方法进行涡轮叶片的谐响应分析可比较精确地得到其谐响应特性.      

周向非均匀流动引起的转子动力学载荷分析

以轴流涡轮为分析对象,分析了流动沿周向分布不均匀的条件下涡轮转子受到的气动载荷。假设涡轮级内的流动沿周向呈1次谐波形式的余弦函数分布,考虑转子叶片叶顶间隙,利用双耦合激励盘模型描述叶片通道中的流动,采用谐波分析方法求得涡轮转子受到的气动载荷。分析了不同轴向位置处、不同非均匀流动参数的影响,讨论了转子偏心所产生的非均匀流动,给出了非稳态流动条件下的气动载荷。结果表明:1次谐波形式的非均匀流动会产生附加载荷作用于转子,从而影响转子的动力学性能;非稳态不均匀流动会影响附加载荷的大小,但不改变其在动坐标系中的方向。     

风力机叶片表面压力的计算与外场测试分析

研究在外场工况下,对风力机叶片表面压力的测试方法,并将测试结果与CFD计算结果进行比较。为了获取叶片在外场非稳态工况下的压力信息,沿叶片展向选取7个典型段面布置带式压力传感器。在数值计算中,通过数码扫描得到试验翼型的几何形状并建立计算模型,用非压缩的N—S方程和SSTk-w。湍流模型耦合,分别对7个翼型的气动性能进行计算。通过对试验和计算结果的对比发现,因为三维旋转效应的存在,基于动量叶素理论的二维翼型计算常常低估了实际风轮动力的产生,旋转叶轮表面压力分布和二维翼型计算结果明显不同。     

叶尖射流对风力机叶尖流场影响的数值研究

为了设计出更加适合非并网系统的风力机,采用在叶尖加入射流的方法来改变叶尖流场分布.在风力机叶片顶端沿弦长布置3个喷口,采用CFD数值模拟方法,通过改变风力机转速获得原型和带喷口的风力机模型的气动特性以及流场分布.发现在转速低于1 200r/min时,带有安装在不同位置的喷口的风力机功率增长率几乎都为零,射流在这一转速范围内对风力机的气动性能几乎没有影响.而转速高于1 200r/min时,随着转速的增大,喷口位于叶尖中部的风力机的功率增长率快速地增大,射流影响了75%以上叶高的表面的压力分布,在大转速下吸力面低压区范围较大,其叶尖涡涡量低于其他方案中,并且在下游扩散得比其他方案快,改善了风力机下游流场,提高了风力机效率.喷口布置在叶尖前缘时其叶尖涡的局部涡量较原型叶片稍大,降低了风力机功率的输出.喷口布置在叶尖尾缘时基本和原型叶片相同.该结论为设计适用于非并网系统的定桨距变转速风力机提供了基础.     

基于人群搜索算法的Dykas汽轮机叶栅型线重构

针对汽轮机叶片智能优化设计中的叶片参数化建模与型线重构问题,首先基于三次Bezier曲线推导了叶片型线的坐标方程,根据叶片基本几何参数实现了汽轮机叶片型线的参数化表达。在构建的参数化模型的基础上,根据已知的叶型坐标点数据,采用人群搜索算法(SOA)重构实际汽轮机叶片型线,并通过叶型重构算例对SOA和粒子群算法(PSO)进行了对比分析。结果表明：基于Bezier曲线构建的参数化型线光滑性好质量高,所提方法可以任意修改叶型几何参数并再生模型,实用性较强且效率较高。基于SOA优化算法的叶片型线重构方法可以准确重构Dykas汽轮机叶片型线,且SOA重构叶型的方法收敛速度快、收敛结果更稳定,为汽轮机叶片优化设计奠定基础。     

风力机标模非定常数值模拟中的影响因素研究

基于自主研发的“亚跨超 CFD 软件平台”(TRIP3.0)，采用刚性运动网格技术和动态拼接网格技术，开发了针对旋转类机械的非定常求解模块。本文开展了 NREL Phase VI 风力机标模非定常数值模拟中的影响因素研究，影响因素主要包括不同时间步长的影响、不同湍流模型的影响、刚性动网格技术和动态拼接网格技术的影响三个方面。本文数值模拟采用的控制方程为雷诺平均 N-S 方程，采用有限体积法离散控制方程，离散方程的时间方向采用“双时间步”方法求解，空间方向无粘项离散采用 MUSCL-Roe 格式，湍流模型采用 SA 和 SST 湍流模型，并引入多重网格和并行计算技术加速求解。数值模拟结果表明：时间步长、湍流模型、网格方案等因素主要影响风力机叶片吸力面的流动结构，进而影响吸力面的压力分布，而对压力面的流动结构和压力分布基本没有影响；采用刚性运动网格结合 SA 湍流模型所得到的压力分布更接近实验值。     

基于广义力的旋翼振动载荷计算

为描述挥舞、摆振铰和变距轴承的运动引入了三个刚体运动自由度 ,旋翼桨叶通过 5节点15自由度有限单元离散 ,计入了桨叶刚性运动与非线性弹性变形之间的动力学耦合效应。利用曲线坐标系下的本构方程 ,对经典的中等变形梁理论进行了重新的推导。另外 ,采用了Leishman Beddoes非定常和动态失速模型 ,入流由自由尾迹分析获得。导出的旋翼桨叶非线性时变常微分方程以广义力的形式给出。桨叶截面载荷和运动方程在位形空间中同时求解。由本文分析得出的桨叶振动载荷与SA349/ 2小羚羊直升机飞行测试数据吻合程度很好。     

涡轮复杂气冷叶盘结构变形分析模型简化方法

基于高压转子开展高压(HP)涡轮转子叶片叶尖变形分析可提高叶尖间隙的数值模拟精度，而高压涡轮转子叶片由于其复杂的气冷结构，有限元分析网格数量巨大；叶片和轮盘的榫接结构属于非线性分析，也需要足够的计算机时。针对该问题提出了一种复杂气冷叶片的简化方法和榫接结构接触计算简化方法，在不影响计算精度的前提下提高计算效率。采用该方法对典型结构高压涡轮转子进行了变形分析，与采用复杂气冷叶片模型和接触分析方法的变形分析结果进行比较。结果表明:涡轮叶片叶尖最大径向变形相对误差为0.47%，计算机时减少99%，证明简化方法和计算方法的有效性。