CN110032769B - 一种通风机叶片流固耦合模型 - Google Patents
一种通风机叶片流固耦合模型 Download PDFInfo
- Publication number
- CN110032769B CN110032769B CN201910200932.3A CN201910200932A CN110032769B CN 110032769 B CN110032769 B CN 110032769B CN 201910200932 A CN201910200932 A CN 201910200932A CN 110032769 B CN110032769 B CN 110032769B
- Authority
- CN
- China
- Prior art keywords
- fluid
- equation
- blade
- solid
- domain
- 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.)
- Active
Links
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]
Abstract
本发明公开了一种通风机叶片流固耦合数值分析方法,该方法的通过分别建立流体域、固体域的动力学模型和基本方程、边界条件以及湍流模型,将方程代入流固耦合模型流程图。在耦合边界上固体的振动引起流体的运动,而流体运动时的压力又反作用于固体,形成两者间的耦合作用,建立叶片流固耦合问题的双向耦合模型。本发明与已有的叶片流固耦合模型相比,更加趋近于实际,不单纯将流体的作用处理为阻尼作用,把流体域和固体域视为一个完整的工程问题。
Description
技术领域
本发明涉及的是一种通风机叶片流固耦合模型,具体而言,本发明是有关于通风机叶片动力学分析建模,属于流固耦合力学领域。
背景技术
通风机因其风压稳定、风量大、效率高、噪声小等特点,广泛应用于工厂、矿井、隧道、冷却塔、船舶和建筑物等场所。叶片是通风机的关键部件,承受离心力、流体动力、振动、温差、介质等的综合作用。
通风机的工作介质是气体,因此风机内部的各种流固耦合现象对故障的产生和发展有很大的影响,增大了准确诊断通风机叶片故障的难度。在通风机的各种故障中,叶片失效占很大比例,主要故障类型为:叶片裂纹故障和叶片折断故障。其中失速颤振是风机最为多发的颤振故障,其起因就是由于叶片本身震荡与其激发流场非定常性之间的交互作用。这种气动力和叶片振动特性的相互作用常会导致通风机的重大事故,也给企业生产带来极大的经济损失。
叶片颤振及故障是典型的流固耦合问题,由于其属于边缘科学涉及诸多学科分支,对其的研究还不是很深入。一般将机械结构视作固体结构,它们的振动分析与流体动力学问题分析都有各自的求解理论与方法。当固体与流体相互耦合时,在耦合边界上固体的振动引起流体的运动,而流体运动时的压力又反作用于固体,形成两者间的耦合作用。传统意义上的流固耦合模型求解是将叶片振动视为主振体,而把流体的作用处理成研究阻尼,再根据叶片振动的条件求流体运动。另一种方法是将固体与流体分别视作一个系统内的两大子结构,流体的动力学方程与固体的动力学方程单独建立,在流体与固体相接边界上,引进相互作用的力作为广义坐标,得到较为简便的流固耦合系统振动问题求解模型。以上模型大部分为单向耦合模型,由分析可知,建立趋近叶片实际工况、双向耦合的流固耦合模型具有重要意义。
发明内容
本发明能建立趋近与实际工况的流固耦合模型,实现双向耦合。
为实现上述目的。本发明采用如下的技术方案。具体的实现步骤如下所述:
步骤1:建立固体域控制方程,通风机叶片的基本方程为两个自由度的振动方程,分为叶片的垂直运动变化量α和扭转运动的变化角度β,其无量纲的形式的简化方程如下:
其中:Sβ和Iβ为叶片的无量纲静矩和转矩,Cl为升力系数,Cm为转矩系数,M∞为来流马赫数,ωα和ωβ分别是两个方向上的固有频率,μ为质量比,U为来流速度,KS叶片刚度系数矩阵,MS是叶片的质量矩阵,β0是叶片的静态攻角,FS是叶片的力矩阵。
步骤2:建立流体域控制方程,假设流体是均匀、小扰动的流动,建立流体域的运动方程和连续条件方程为
其中:P是叶片工作中的气动压力,u,v,w分别为流体单元沿x,y,z三个自由度的位移分量,ρ是流体的密度,Kf流体的压缩模量;
步骤3:湍流方程采用的是低雷诺数双方程模型q-ω,计算中的方程模型为式(3);
其中:μt是流体域流体的黏性,是流体域的速度分量,p是流体域气体的压力值,q是流体域种气体作的功,k是流体域中湍流动能,Re是特征雷诺数,J是叶片工作中的湍流速度场的应变率,ω是紊流动能的耗散率,lω是流体域当前单元到最近固体面的距离,θ,C1,C2,C3,Cμ为封闭常数。
步骤4:边界耦合条件
在流体域和固体域的交界面,要同时满足运动学和动力学的条件,法向速度和法向虚位移一至;
步骤5:耦合算法
对建立的数学模型分别进行离散化,利用伽辽金法对方程(1)、(2)化简为有限元方程。将方程各式带入耦合算法流程图,建立叶片得双向流固耦合模型。
本发明对于叶片流固耦合分析模型具有以下显著特点
1、叶片的流固耦合模型采用双向耦合计算。传统的单向耦合计算仅仅将叶片振动视为主振体,而把流体的作用处理成研究阻尼,再根据叶片振动的条件求流体运动。双向耦合将流体域和固体域视为一个完整的工程问题。
2、叶片的固体域模型不采用悬臂梁结构,将其运动方程视为两个自由度的振动方程。
附图说明
图1为本发明的流固耦合模型流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细描述。
步骤1:建立固体域控制方程,通风机叶片的基本方程为两个自由度的振动方程,分为叶片的垂直运动变化量α和扭转运动的变化角度β,其无量纲的形式的简化方程如下:
其中:Sβ和Iβ为叶片的无量纲静矩和转矩,Cl为升力系数,Cm为转矩系数,M∞为来流马赫数,ωα和ωβ分别是两个方向上的固有频率,μ为质量比,U为来流速度,KS叶片刚度系数矩阵,MS是叶片的质量矩阵,β0是叶片的静态攻角,FS是叶片的力矩阵。
步骤2:建立流体域控制方程,假设流体是均匀、小扰动的流动,建立流体域的运动方程和连续条件方程为
其中:P是叶片工作中的气动压力,u,v,w分别为流体单元沿x,y,z三个自由度的位移分量,ρ是流体的密度,Kf流体的压缩模量;
步骤3:湍流方程采用的是低雷诺数双方程模型q-ω,计算中的方程模型为式(3);
其中:μt是流体域流体的黏性,是流体域的速度分量,p是流体域气体的压力值,q是流体域种气体作的功,k是流体域中湍流动能,Re是特征雷诺数,J是叶片工作中的湍流速度场的应变率,ω是紊流动能的耗散率,lω是流体域当前单元到最近固体面的距离,θ,C1,C2,C3,Cμ为封闭常数。
步骤4:边界耦合条件
在流体域和固体域的交界面,要同时满足运动学和动力学的条件,法向速度和法向虚位移一至;
步骤5:耦合算法
对建立的数学模型分别进行离散化,利用伽辽金法对方程(1)、(2)化简为有限元方程。
2、根据权利要求1所述的一种通风机叶片流固耦合数值分析方法,其特征在于:
本方法实施的具体步骤如下:
步骤1中采用4阶龙格—库塔法求解非线性微分方程(1);
其中r1=α,r2=β,Δt是时间步长。
步骤2中,流体域中,将控制方程无量纲化可得非线性微分方程,同时为了边界条件得处理方便,将方程转化到任意坐标系中,并且进行线性化,化简得到方程式(5);利用隐式得时间推进有限元差分方法克朗可—尼科尔森方法求解,牛顿-迭代使每一个时间步内收敛;
将基本方程各式代入叶片流固耦合算法中,实现研究叶片振动问题的流固耦合数值分析方法。
Claims (2)
1.一种通风机叶片流固耦合数值分析方法,其特征在于:该分析方法包含以下几个步骤,
步骤1:建立固体域控制方程,通风机叶片的基本方程为两个自由度的振动方程,分为叶片的垂直运动变化量α和扭转运动的变化角度β,其无量纲的形式的简化方程如下:
其中:Sβ和Iβ为叶片的无量纲静矩和转矩,Cl为升力系数,Cm为转矩系数,M∞为来流马赫数,ωα和ωβ分别是两个方向上的固有频率,μ为质量比,U为来流速度,KS叶片刚度系数矩阵,MS是叶片的质量矩阵,β0是叶片的静态攻角,FS是叶片的力矩阵;
步骤2:建立流体域控制方程,假设流体是均匀、小扰动的流动,建立流体域的运动方程和连续条件方程为
其中:P是叶片工作中的气动压力,u,v,w分别为流体单元沿x,y,z三个自由度的位移分量,ρ是流体的密度,Kf流体的压缩模量;
步骤3:湍流方程采用的是低雷诺数双方程模型q-ω,计算中的方程模型为式(3);
其中:μt是流体域流体的黏性,是流体域的速度分量,q是流体域种气体作的功,k是流体域中湍流动能,Re是特征雷诺数,J是叶片工作中的湍流速度场的应变率,ω是紊流动能的耗散率,lω是流体域当前单元到最近固体面的距离,θ,C1,C2,C3,Cμ为封闭常数;
步骤4:边界耦合条件
在流体域和固体域的交界面,要同时满足运动学和动力学的条件,法向速度和法向虚位移一至;
其中,nf是流体边界的单位法向量,νSn,νfn分别是流体域、固体域的法向速度;在这里,下角标s,f分别代表叶片的固体域和流体域参数;
步骤5:耦合算法
对建立的数学模型分别进行离散化,利用伽辽金法对方程(1)、(2)化简为有限元方程。
2.根据权利要求1所述的一种通风机叶片流固耦合数值分析方法,其特征在于:
本方法实施的具体步骤如下:
步骤1中采用4阶龙格—库塔法求解非线性微分方程(1);
其中r1=α,r2=β,Δt是时间步长;
步骤2中,流体域中,将控制方程无量纲化可得非线性微分方程,同时为了边界条件得处理方便,将方程转化到任意坐标系中,并且进行线性化,化简得到方程式(5);利用隐式得时间推进有限元差分方法克朗可—尼科尔森方法求解,牛顿-迭代使每一个时间步内收敛;
将基本方程各式代入叶片流固耦合算法中,实现研究叶片振动问题的流固耦合数值分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910200932.3A CN110032769B (zh) | 2019-03-18 | 2019-03-18 | 一种通风机叶片流固耦合模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910200932.3A CN110032769B (zh) | 2019-03-18 | 2019-03-18 | 一种通风机叶片流固耦合模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110032769A CN110032769A (zh) | 2019-07-19 |
CN110032769B true CN110032769B (zh) | 2023-05-02 |
Family
ID=67236189
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910200932.3A Active CN110032769B (zh) | 2019-03-18 | 2019-03-18 | 一种通风机叶片流固耦合模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110032769B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111854070B (zh) * | 2020-06-09 | 2022-03-08 | 海信(山东)空调有限公司 | 一种空调及风量控制方法 |
CN112380747B (zh) * | 2020-11-13 | 2022-05-24 | 中国电力工程顾问集团西南电力设计院有限公司 | 一种双曲线型钢结构冷却塔设计方法 |
CN113591416A (zh) * | 2021-08-10 | 2021-11-02 | 西安西热节能技术有限公司 | 一种适用于汽轮机叶片流热固耦合分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882177A (zh) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | 叶轮机械计入叶间相位角的气弹稳定性流固耦合预测方法 |
CN107506562A (zh) * | 2017-09-29 | 2017-12-22 | 西安科技大学 | 一种水润滑橡胶轴承双向热流固耦合计算方法 |
CN107895093A (zh) * | 2017-12-08 | 2018-04-10 | 兰州理工大学 | 一种风力机叶片流固耦合模态设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10698980B2 (en) * | 2015-10-10 | 2020-06-30 | John Michael Snider, SR. | Methods for constructing surfaces for optimizing fluid flow |
-
2019
- 2019-03-18 CN CN201910200932.3A patent/CN110032769B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882177A (zh) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | 叶轮机械计入叶间相位角的气弹稳定性流固耦合预测方法 |
CN107506562A (zh) * | 2017-09-29 | 2017-12-22 | 西安科技大学 | 一种水润滑橡胶轴承双向热流固耦合计算方法 |
CN107895093A (zh) * | 2017-12-08 | 2018-04-10 | 兰州理工大学 | 一种风力机叶片流固耦合模态设计方法 |
Non-Patent Citations (2)
Title |
---|
复合材料风力机叶片流固耦合分析方法研究;胡芳琳等;《太阳能学报》;20141028(第10期);第1917-1922页 * |
通风机故障叶片的模态分析;付胜 等;《矿山机械》;20150430;第25-28页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110032769A (zh) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110032769B (zh) | 一种通风机叶片流固耦合模型 | |
CN107895069B (zh) | 一种基于复合材料结构的流固耦合数值预测方法 | |
Smith et al. | CFD-based analysis of nonlinear aeroelastic behavior of high-aspect ratio wings | |
CN110287643B (zh) | 空腔流致振动与流致噪声耦合特性的风洞实验模拟方法 | |
CN110162826B (zh) | 薄壁结构热气动弹性动响应分析方法 | |
CN105134482B (zh) | 大型智能风机叶片系统灰色组合建模与优化振动控制的方法 | |
CN109840349A (zh) | 一种固定翼飞机阵风响应建模分析方法 | |
CN110889169A (zh) | 基于多体系统传递矩阵法的舵面系统非线性颤振模型建模方法 | |
Zare et al. | Airborne dust-induced performance degradation in NREL phase VI wind turbine: a numerical study | |
CN111723438A (zh) | 一种消除点阵夹芯板结构热屈曲及抑制非线性颤振的方法 | |
Firouz-Abadi et al. | Analysis of non-linear aeroelastic response of a supersonic thick fin with plunging, pinching and flapping free-plays | |
Amini et al. | Piezoelectric energy harvesting from vertical piezoelectric beams in the horizontal fluid flows | |
Virdis et al. | A lumped parameter model for the analysis of dynamic effects in Wells turbines | |
CN108287948B (zh) | 一种基于空腔多场耦合方程的边界条件建立方法 | |
Castorrini et al. | Numerical study on the passive control of the aeroelastic response in large axial fans | |
Kabir et al. | Aerodynamic Analysis on Double Wedge Airfoil at Different Mach Numbers with Varying Angle of Attacks Using Computational Fluid Dynamics | |
Lahooti et al. | Thick strip method for efficient large-eddy simulations of flexible wings in stall | |
Rajaomazava III et al. | Numerical analysis of hydrofoil dynamics by using a fluid-structure interaction approach | |
Bekhti et al. | Numerical analysis of dynamic stall on wind turbine airfoils | |
Salman et al. | Structural Nonlinearities and their Impact on the Fidelity of Critical Steady Maneuver Loads and Trimming Con figuration of Very Flexible Airframes | |
Ramaiah et al. | Aero-elastic analysis of stiffened composite wing structure | |
Qi et al. | Unsteady aerodynamics of offshore floating wind turbines using free vortex wake model | |
Kim et al. | Hydrodynamic analysis of ducted marine turbine: Study on the duct and turbine efficiency | |
Karthik et al. | Modeling and simulation studies on aeroelasticity of an airfoil with simulink | |
Guangzhi et al. | Transient aerodynamic characteristics of a moving train on bridge in tornado environment |
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 |