CN108304643A - 一种空腔流固耦合方程的无量纲化方法 - Google Patents

一种空腔流固耦合方程的无量纲化方法 Download PDF

Info

Publication number
CN108304643A
CN108304643A CN201810084218.8A CN201810084218A CN108304643A CN 108304643 A CN108304643 A CN 108304643A CN 201810084218 A CN201810084218 A CN 201810084218A CN 108304643 A CN108304643 A CN 108304643A
Authority
CN
China
Prior art keywords
cavity
equation
nondimensional
flow
fluid structurecoupling
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.)
Granted
Application number
CN201810084218.8A
Other languages
English (en)
Other versions
CN108304643B (zh
Inventor
王显圣
杨党国
刘俊
周方奇
施傲
杨野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Original Assignee
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center filed Critical High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority to CN201810084218.8A priority Critical patent/CN108304643B/zh
Publication of CN108304643A publication Critical patent/CN108304643A/zh
Application granted granted Critical
Publication of CN108304643B publication Critical patent/CN108304643B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design

Abstract

本发明公开了一种空腔流固耦合方程的无量纲化方法，包括以下步骤：步骤S1：建立空腔流固耦合问题的关键影响参数系统；步骤S2：建立空腔流固耦合问题的表征参数系统；步骤S3：建立无量纲的空腔流固耦合问题的表征参数系统；步骤S4：建立空腔流固耦合方程和边界条件；步骤S5：建立无量纲的空腔流固耦合方程和无量纲的空腔流动边界远场条件。本发明有利于抓住空腔流固耦合问题的本质和关键参数；有利于研究人员对空腔流固耦合问题的无量纲表征参数的选择；提高空腔流固耦合实验和数值研究的效率；能够促进研究人员对于空腔流固耦合问题的理解。

Description

一种空腔流固耦合方程的无量纲化方法
技术领域
本发明涉及航空航天、气动声学和结构动力学交叉技术领域，具体的说，是一种空腔流固耦合方程的无量纲化方法。
背景技术
空腔类型的结构布局在航空航天飞行器中普遍存在，例如内埋武器舱、飞行器起落架舱等。在高速来流条件下，空腔内易发生流激振荡现象，导致空腔壁板结构产生振动，形成典型的空腔流固耦合问题，导致空腔问题的复杂程度显著增加。空腔流动与结构振动方程成为研究空腔流固耦合问题的重要理论基础，不过，传统的空腔流固耦合方程中涉及的物理参数通常包含量纲，而包含量纲的物理参数之间的大小关系不易确定，空腔流固耦合方程不易进行化简，使研究人员不易对空腔流固耦合规律进行理论建模，同时也增加了数值模拟过程中的数值误差，给空腔流固耦合问题的研究带来困难。
目前，国内外研究机构外对于空腔流场变量的无量纲化方法的研究较多，主要通过组合这些流场变量进行消除量纲操作，不过对于空腔流动与结构振动的流固耦合方程的无量纲化方法研究较少，相关研究能力与开发出工程实用的空腔类飞行器部件之间存在较大差距。传统的空腔物理场表征参数的无量纲化方法主要从表征现象的物理量的量纲出发，通过量纲分析消除物理参数的量纲，这种消除物理参数量纲的方法存在诸多不足，一是传统方法需要研究人员对表征物理现象的物理量具有丰富的先验知识，研究人员利用这些先验知识选择物理量的组合方式，这种组合方式存在较大随意性；二是传统方法没有从表征物理现象基本规律的控制方程出发，无法对空腔流固耦合方程进行无量纲化；三是传统方法对于物理参数无量纲化过程的物理含义不够明确，因此无法得到无量纲量之间的函数关系。
发明内容
本发明的目的在于提供一种空腔流固耦合方程的无量纲化方法，解决现有技术对于空腔流固耦合问题研究能力的不足。
本发明通过下述技术方案实现：一种空腔流固耦合方程的无量纲化方法，具体包括以下步骤：
步骤S1：建立空腔流固耦合问题的关键影响参数系统；
步骤S2：建立空腔流固耦合问题的表征参数系统；具体包括以下步骤：
步骤S21：根据空腔流动与结构振动的基本特性，建立表征空腔流固耦合问题的物理参数系统；
步骤S22：根据表征空腔流固耦合问题的物理参数系统中各物理参数的物理含义，结合气动声学与结构动力学先验知识，选择物理参数的特征尺度；
步骤S3：建立无量纲的空腔流固耦合问题的表征参数系统；
步骤S4：建立空腔流固耦合方程和边界条件；
步骤S5：建立无量纲的空腔流固耦合方程和无量纲的空腔流动边界远场条件。
进一步地，为了更好的实现本发明，所述步骤S1具体是指：分析空腔流动与结构振动的基本规律，建立空腔流固耦合问题的关键影响参数系统；所述关键影响参数系统具体包括表征空腔几何参数、表征空腔结构参数和表征来流条件参数；
所述表征空腔几何参数包括空腔长度L、空腔深度D以及空腔宽度W；
所述表征空腔结构参数包括空腔壁板材料的杨氏模量E、空腔壁板材料的泊松比ν、空腔壁板材料的密度ρw以及空腔壁板的厚度h；
所述表征来流条件参数包括来流介质压力ps、来流介质密度ρs、来流介质温度Ts、来流介质粘性μ、来流介质比热比γ、来流介质热传导系数ks、来流介质的定容比热Cν、来流介质的定压比热Cp、来流介质的气体常数R、来流速度Us、来流速度方向的三个方位角αl(l＝1,2,3)、速度边界层厚度δ、位移边界层厚度δ*以及动量边界层厚度θ。
进一步地，为了更好的实现本发明，所述步骤S2中表征空腔多物理场特性的基本物理参数包括空间位置坐标Xk、时间变量t、空腔流场速度分量uk、空腔流场温度T、空腔流场密度ρ、空腔流场压力p、空腔壁板驱动力pw、空腔壁板切向的面内坐标α、空腔壁板切向的面内坐标β以及空腔壁板中平面离面位移η；
所述特征尺度包括空腔长度L、流动特征频率fc、来流速度Us、来流温度Ts、来流密度ρs、来流速压ρsUs 2/2、来流介质压力ps以及空腔壁板厚度h；
其中，下标k(k＝1,2,3)表示物理量在三个空间方向的分量，下标s表示该物理量对应来流介质的物理参数。
进一步地，为了更好的实现本发明，所述步骤S3具体是指：根据空腔流动与结构振动的表征物理参数对应的特征尺度，对表征物理参数进行无量纲化，获得无量纲的空腔流固耦合问题的表征参数系统；所述无量纲的空腔流固耦合问题的表征参数系统包括空间位置坐标的无量纲参数时间变量的无量纲参数空腔流场速度的无量纲参数空腔流场压力的无量纲参数空腔流场密度的无量纲参数空腔流场温度的无量纲参数、空腔壁板驱动力的无量纲参数空腔壁板两个面内坐标的无量纲参数以及空腔壁板中面位移的无量纲参数
进一步地，为了更好的实现本发明，所述步骤S4具体是指：根据空腔流固耦合的基本规律，建立空腔流动与结构振动的耦合方程和每个物理场对应的边界条件：具体包括以下步骤：
步骤S41：建立空腔流动控制方程及空腔流动边界条件；
步骤S42：建立空腔壁板结构振动控制方程及空腔壁板结构振动边界条件。
进一步地，为了更好的实现本发明，所述步骤S41中空腔流动控制方程为：
其中，式中：
e为来流介质内能；
下标k表示物理量在空间三个方向的分量；
下标i、下标j和下标r为哑标；
所述步骤S41中空腔流动边界条件包括空腔流动边界壁面条件和空腔流动边界远场条件；所述空腔流动边界壁面条件的控制方程为：
其中，式中：
Dw为空腔壁面区域；
下标k为物理量在三个空间坐标方向的分量；
nlk表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
所述空腔流动边界远场条件的控制方程为：
其中，式中：
D∞表示远离空腔壁面的空间区域。
进一步地，为了更好的实现本发明，所述步骤S42中具体包括空腔壁板结构振动控制方程为：
其中，式中：
为空腔壁板的抗弯刚度；
pw为空腔壁板驱动应力；
ρw为空腔壁板材料的密度；
所述空腔壁板结构振动边界条件包括空腔壁板结构壁面条件和空腔壁板结构壁面交界线条件；
所述空腔壁板结构壁面条件的控制方程为：
其中，式中：
下标i为哑标；
nli表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
所述空腔壁板结构壁面交界线条件的控制方程为：
其中，式中：Dln表示空腔壁板之间交界线所在空间区域，Λ表示空腔壁板中平面上垂直于空腔壁板之间交界线的空间方向。
进一步地，为了更好的实现本发明，所述步骤S5具体是指：根据有量纲的空腔流固耦合方程与边界条件，依据无量纲的空腔流固耦合问题的表征参数系统，建立空腔流动与结构振动的无量纲方程和无量纲边界条件，获得表征空腔流固耦合问题的无量纲参数之间的关系；具体包括以下步骤：
步骤S51：建立无量纲的空腔流动控制方程及无量纲的空腔流动边界条件；
步骤S52：无量纲的空腔流动边界条件。
进一步地，为了更好的实现本发明，所述步骤S51中的无量纲的空腔流动控制方程具体是指：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(1)，并且对方程中的项进行化简，得到无量纲的空腔流动控制方程；所述无量纲的空腔流动控制方程为：
其中，式中：
所述步骤S51中的无量纲的空腔流动边界条件，具体是指：所述无量纲的空腔流动边界条件包括无量纲的空腔流动边界壁面条件和无量纲的空腔流动边界远场条件；
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(2)得到无量纲的空腔流动边界壁面条件的控制方程，所述无量纲的空腔流动边界壁面条件的控制方程为：
其中，式中：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(3)中得到无量纲的空腔流动边界远场条件，所述无量纲的空腔流动边界远场条件的控制方程为：
进一步地，为了更好的实现本发明，所述无量纲的空腔流动控制方程具体是指：步骤S52具体包括空腔壁板结构振动控制方程的无量纲化和空腔壁板结构振动边界条件的无量纲化；
所述空腔壁板结构振动控制方程的无量纲化具体是指：将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(4)中，并且对方程中的项进行化简，得到无量纲的空腔流动壁板结构振动控制方程，所述无量纲的空腔壁板结构振动控制方程为：
其中，式中：
所述空腔壁板结构振动边界条件的无量纲化具体是指：无量纲的空腔壁板结构振动边界条件包括无量纲的空腔壁板结构壁面条件和无量纲的空腔壁板结构壁面交界线条件；
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(5)，并且对方程中的项进行化简，得到无量纲的空腔壁板结构壁面条件的控制方程，所述无量纲的空腔壁板结构壁面条件的控制方程为：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(5)，并且对方程中的项进行化简，得到无量纲的空腔壁板结构壁面交界线条件的控制方程，所述无量纲的空腔壁板结构壁面交界线条件的控制方程为：
本发明与现有技术相比，具有以下优点及有益效果：
(1)本发明根据空腔流动与结构振动耦合特性，从流体力学和结构动力学基本原理出发，建立空腔流固耦合方程及边界条件，能够从根源上保证所建立方程的正确性，有利于抓住空腔流固耦合问题的本质和关键参数；
(2)本发明根据空腔流固耦合问题表征物理参数的特征量对空腔流固耦合方程及边界条件进行无量纲化，降低无量纲化过程中不同物理量之间组合的随意性，有利于研究人员对空腔流固耦合问题的无量纲表征参数的选择；
(3)本发明所建立的无量纲化的空腔流固耦合方程及边界条件，减小了空腔流固耦合问题表征物理参数的个数，对于空腔流固耦合问题的实验和数值研究具有积极的指导作用，能够提高空腔流固耦合实验和数值研究的效率；
(4)本发明建立的空腔流固耦合方程及边界条件无量纲化方法是基于空腔流固耦合方程及边界条件的等价变换，研究人员在获得无量纲的空腔流固耦合方程及边界条件的同时也获得无量纲的空腔流固耦合问题表征参数之间的方程关系，能够促进研究人员对于空腔流固耦合问题的理解。
附图说明
图1为本发明中空腔流固耦合方程坐标系示意图。
具体实施方式
下面结合实施例对本发明作进一步地详细说明，但本发明的实施方式不限于此。
实施例1：
本发明通过下述技术方案实现，如图1所示，一种空腔流固耦合方程的无量纲化方法，具体包括以下步骤：
步骤S1：建立空腔流固耦合问题的关键影响参数系统；
步骤S2：建立空腔流固耦合问题的表征参数系统；具体包括以下步骤：
步骤S21：根据空腔流动与结构振动的基本特性，建立表征空腔流固耦合问题的物理参数系统；
步骤S22：根据物理参数系统建立表征空腔流固耦合问题的物理参数系统中各物理参数的物理含义，结合气动声学与结构动力学先验知识，选择物理参数的特征尺度；
步骤S3：建立无量纲的空腔流固耦合问题的表征参数系统；
步骤S4：建立空腔流固耦合方程和边界条件；
步骤S5：建立无量纲的空腔流固耦合方程和无量纲的空腔流动边界远场条件。
本发明根据空腔流动与结构振动耦合特性，从流体力学和结构动力学基本原理出发，建立空腔流固耦合方程及边界条件，能够从根源上保证所建立方程的正确性，有利于抓住空腔流固耦合问题的本质和关键参数；根据空腔流固耦合问题表征物理参数的特征量对空腔流固耦合方程及边界条件进行无量纲化，降低无量纲化过程中不同物理量之间组合的随意性，有利于研究人员对空腔流固耦合问题的无量纲表征参数的选择；所建立的无量纲化的空腔流固耦合方程及边界条件，减小了空腔流固耦合问题表征物理参数的个数，对于空腔流固耦合问题的实验和数值研究具有积极的指导作用，能够提高空腔流固耦合实验和数值研究的效率；建立的空腔流固耦合方程及边界条件无量纲化方法是基于空腔流固耦合方程及边界条件的等价变换，研究人员在获得无量纲的空腔流固耦合方程及边界条件的同时也获得无量纲的空腔流固耦合问题表征参数之间的方程关系，能够促进研究人员对于空腔流固耦合问题的理解。
本实施例的其他部分与上述实施例相同，故不再赘述。
实施例2：
本实施例在上述实施例的基础上做进一步优化，如图1所示，一种空腔流固耦合方程的无量纲化方法，具体包括以下步骤：
步骤S1：建立空腔流固耦合问题的关键影响参数系统；具体是指：分析空腔流动与结构振动的基本规律，建立空腔流固耦合问题的关键影响参数系统；所述关键影响参数系统具体包括表征空腔几何参数、表征空腔结构参数和表征来流条件参数；
所述表征空腔几何参数包括空腔长度L、空腔深度D以及空腔宽度W；
所述表征空腔结构参数包括空腔壁板材料的杨氏模量E、空腔壁板材料的泊松比ν、空腔壁板材料的密度ρw以及空腔壁板的厚度h；
所述表征来流条件参数包括来流介质压力ps、来流介质密度ρs、来流介质温度Ts、来流介质粘性μ、来流介质比热比γ、来流介质热传导系数ks、来流介质的定容比热Cν、来流介质的定压比热Cp、来流介质的气体常数R、来流速度Us、来流速度方向的三个方位角αl(l＝1,2,3)、速度边界层厚度δ、位移边界层厚度δ*以及动量边界层厚度θ。
步骤S2：建立空腔流固耦合问题的表征参数系统；具体包括以下步骤：
步骤S21：根据空腔流动与结构振动的基本特性，建立表征空腔流固耦合问题的物理参数系统；
步骤S22：根据表征空腔流固耦合问题的物理参数系统中各物理参数的物理含义，结合气动声学与结构动力学先验知识，选择物理参数的特征尺度；所述步骤S2中表征空腔多物理场特性的基本物理参数包括空间位置坐标Xk、时间变量t、空腔流场速度分量uk、空腔流场温度T、空腔流场密度ρ、空腔流场p压力、空腔壁板驱动力pw、空腔壁板切向的面内坐标α、空腔壁板切向的面内坐标β以及空腔壁板中平面离面位移η；
所述特征尺度包括空腔长度L、流动特征频率fc、来流速度Us、来流温度Ts、来流密度ρs、来流速压ρsUs 2/2、来流介质压力ps以及空腔壁板厚度h；
其中，式中：下标k(k＝1,2,3)表示物理量在三个空间方向的分量；
下标s表示该物理量对应来流介质的物理参数。
步骤S3：建立无量纲的空腔流固耦合问题的表征参数系统；具体是指：
根据空腔流动与结构振动的表征物理参数对应的特征尺度，对表征物理参数进行无量纲化，获得无量纲的空腔流固耦合问题的表征参数系统；所述无量纲的空腔流固耦合问题的表征参数系统包括空间位置坐标的无量纲参数时间变量的无量纲参数空腔流场速度的无量纲参数空腔流场压力的无量纲参数空腔流场密度的无量纲参数空腔流场温度的无量纲参数、空腔壁板驱动力的无量纲参数空腔壁板两个面内坐标的无量纲参数以及空腔壁板中面位移的无量纲参数
步骤S4：建立空腔流固耦合方程和边界条件；具体是指：根据空腔流固耦合的基本规律，建立空腔流动与结构振动的耦合方程和每个物理场对应的边界条件：具体包括以下步骤：
步骤S41：建立空腔流动控制方程及空腔流动边界条件；
步骤S42：建立空腔壁板结构振动控制方程及空腔壁板结构振动边界条件。
进一步地，为了更好的实现本发明，所述步骤S41中空腔流动控制方程为：
其中，式中：e为来流介质内能；
下标k表示物理量在空间三个方向的分量；
下标i、下标j和下标r为哑标；
所述步骤S41中空腔流动边界条件包括空腔流动边界壁面条件和空腔流动边界远场条件；所述空腔流动边界壁面条件的控制方程为：
其中，式中：Dw为空腔壁面区域；
下标k为物理量在三个空间坐标方向的分量；
nlk表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
所述空腔流动边界远场条件的控制方程为：
其中，式中：
D∞表示远离空腔壁面的空间区域。
进一步地，为了更好的实现本发明，所述步骤S42中具体包括空腔壁板结构振动控制方程为：
其中，式中：为空腔壁板的抗弯刚度；
pw为空腔壁板驱动应力；
ρw为空腔壁板材料的密度；
所述空腔壁板结构振动边界条件包括空腔壁板结构壁面条件和空腔壁板结构壁面交界线条件；
所述空腔壁板结构壁面条件的控制方程为：
其中，式中：下标i为哑标；
nli表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
所述空腔壁板结构壁面交界线条件的控制方程为：
其中，式中：Dln表示空腔壁板之间交界线所在空间区域；
Λ表示空腔壁板中平面上垂直于空腔壁板之间交界线的空间方向。
进一步地，为了更好的实现本发明，步骤S5：建立无量纲的空腔流固耦合方程和无量纲的空腔流动边界远场条件，具体是指：根据有量纲的空腔流固耦合方程与边界条件，依据无量纲的空腔流固耦合问题的表征参数系统，建立空腔流动与结构振动的无量纲方程和无量纲边界条件，获得表征空腔流固耦合问题的无量纲参数之间的关系；具体包括以下步骤：
步骤S51：建立无量纲的空腔流动控制方程及无量纲的空腔流动边界条件；所述无量纲的空腔流动控制方程具体是指：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(1)，并且对方程中的项进行化简，得到无量纲的空腔流动控制方程；所述无量纲的空腔流动控制方程为：
其中，式中：
所述无量纲的空腔流动边界条件，具体是指：所述无量纲的空腔流动边界条件包括无量纲的空腔流动边界壁面条件和无量纲的空腔流动边界远场条件；
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(2)得到无量纲的空腔流动边界壁面条件的控制方程，所述无量纲的空腔流动边界壁面条件的控制方程为：
其中，式中：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(3)中得到无量纲的空腔流动边界远场条件，所述无量纲的空腔流动边界远场条件的控制方程为：
所述无量纲的空腔流动控制方程具体是指：步骤S52具体包括空腔壁板结构振动控制方程的无量纲化和空腔壁板结构振动边界条件的无量纲化；
所述空腔壁板结构振动控制方程的无量纲化具体是指将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(4)中，并且对方程中的项进行化简，得到无量纲的空腔壁板结构振动控制方程，所述无量纲的空腔壁板结构振动控制方程为：
其中，式中：
步骤S52：无量纲的空腔流动边界条件；具体是指：无量纲的空腔壁板结构振动边界条件包括无量纲的空腔壁板结构壁面条件和无量纲的空腔壁板结构壁面交界线条件；
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(5)，并且对方程中的项进行化简，得到无量纲的空腔壁板结构壁面条件的控制方程，所述无量纲的空腔壁板结构壁面条件的控制方程为：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(5)，并且对方程中的项进行化简，得到无量纲的空腔壁板结构壁面交界线条件的控制方程，所述无量纲的空腔壁板结构壁面交界线条件的控制方程为：
本实施例的其他部分与上述实施例相同，故不再赘述。
实施例3：
本实施例为本发明的另外一种表述，如图1所示，一种空腔流固耦合方程的无量纲化方法，构建适用于表征空腔流动与结构振动耦合规律的无量纲方程。
本发明通过下述技术方案实现：
1.建立空腔流固耦合问题的关键影响参数系统：
分析空腔流动与结构振动的基本规律，建立空腔流固耦合问题的关键影响参数系统，如表1-表3所示。
根据这些参数的物理含义，将这些参数分成三类：分别为空腔几何参数、空腔结构参数和来流条件参数，其中，
空腔几何参数如表1所示，空腔结构参数如表2所示，来流条件参数如表3所示。
几何参数 物理含义
L 空腔长度
D 空腔深度
W 空腔宽度
表1
结构参数 物理含义
E 空腔壁板材料的杨氏模量
ν 空腔壁板材料的泊松比
ρw 空腔壁板材料的密度
h 空腔壁板的厚度
表2
表3
2.建立空腔流固耦合问题的表征参数系统：
第一步：根据空腔流动与结构振动的基本特性，建立表征空腔流固耦合问题的物理参数系统，如表4所示。
第二步：根据这些物理参数的物理含义，结合气动声学与结构动力学先验知识，选择这些物理参数的特征尺度，如表4所示。
表4为表征空腔多物理场特性的基本物理量及其特征尺度的符号与含义，其中下标k(k＝1,2,3)表示物理量在三个空间方向的分量，下标s表示该物理量对应来流介质的物理参数：
物理量 物理含义 特征尺度 物理含义
xk 空间位置坐标 L 空腔长度
t 时间变量 fc 流动特征频率
uk 空腔流场速度分量 Us 来流介质速度
T 空腔流场温度 Ts 来流介质温度
ρ 空腔流场密度 ρs 来流介质密度
p 空腔流场压力 ps 来流介质压力
pw 空腔壁板驱动力 ps 来流介质压力
α 空腔壁板切向的面内坐标 L 空腔长度
β 空腔壁板切向的面内坐标 L 空腔长度
η 空腔壁板中平面离面位移 h 空腔壁板厚度
表5
3.建立无量纲的空腔流固耦合问题的表征参数系统：
根据空腔流动与结构振动的表征物理参数对应的特征尺度，对这些物理参数进行无量纲化，所获得的无量纲参数如表5所示；表6为表征空腔多物理场特性的无量纲参数及含义：
表5
4.建立空腔流固耦合方程和边界条件：
根据空腔流固耦合的基本规律，建立空腔流动与结构振动的耦合方程和每个物理场对应的边界条件，建立过程共分为两步。
第一步：建立空腔流动控制方程及边界条件：
1)空腔流动控制方程：
其中，式中：
e为来流介质内能；
下标k表示物理量在空间三个方向的分量；
下标i、下标j和下标r为哑标。
2)空腔流动边界条件：
所述空腔流动边界壁面条件的控制方程为:
其中，式中：
Dw表示空腔壁面区域；
nlk表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
下标k为物理量在三个空间坐标方向的分量；
所述空腔流动边界远场条件的控制方程:
其中，式中：
D∞表示远离空腔壁面的空间区域；
第二步：建立空腔壁板结构振动控制方程及边界条件：
1)空腔壁板结构振动控制方程：
其中，式中：为空腔壁板的抗弯刚度；
pw为空腔壁板受到的驱动应力；
ρw为空腔壁板材料的密度。
2)空腔壁板结构振动边界条件：
空腔壁板结构壁面条件的控制方程:
其中，式中：
下标i为哑标；
nli表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
空腔壁板结构壁面交界线条件的控制方程：
其中，式中：
Dln表示空腔壁板之间交界线所在空间区域；
Λ表示空腔壁板中平面上垂直于空腔壁板之间交界线的空间方向。
5.建立无量纲的空腔流固耦合方程和无量纲的空腔流动边界远场条件：
根据有量纲的空腔流固耦合方程与边界条件，依据表5所示的无量纲参数系统，建立空腔流动与结构振动的无量纲方程和无量纲边界条件，获得表征空腔流固耦合问题的无量纲参数之间的关系，建立过程共分为两步。
第一步：建立无量纲的空腔流动控制方程及无量纲的空腔流动边界条件：
1)空腔流动控制方程无量纲化：
将表5所列无量纲参数代入方程(1)中，并且对方程中的项进行化简，得到无量纲的空腔流动控制方程的方程式:
其中，式中：
2)空腔流动边界条件无量纲化：
将表5所列无量纲参数代入方程(2)和方程(3)中，并且对方程中的项进行化简，得到无量纲的空腔流动边界条件。
无量纲的空腔流动边界壁面条件的控制方程的方程式：
其中：式中：
无量纲的空腔流动边界远场条件的控制方程：
第二步：建立无量纲的空腔壁板结构振动方程及边界条件：
1)空腔壁板结构振动控制方程的无量纲化：
将表5所列无量纲参数代入方程(4)中，并且对方程中的项进行化简，得到无量纲的空腔流动控制方程的方程式：
其中，
2)空腔壁板结构振动边界条件的无量纲化：
将表5所列无量纲参数代入方程(5)和方程(6)中，并且对方程中的项进行化简，得到无量纲的空腔壁板结构振动边界条件。
无量纲的空腔壁板结构壁面条件的控制方程：
无量纲的空腔壁板结构壁面交界线条件的控制方程：
以上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化，均落入本发明的保护范围之内。

Claims (10)

1.一种空腔流固耦合方程的无量纲化方法，其特征在于：具体包括以下步骤：
步骤S1：建立空腔流固耦合问题的关键影响参数系统；
步骤S2：建立空腔流固耦合问题的表征参数系统；具体包括以下步骤：
步骤S21：根据空腔流动与结构振动的基本特性，建立表征空腔流固耦合问题的物理参数系统；
步骤S22：根据表征空腔流固耦合问题的物理参数系统中各物理参数的物理含义，结合气动声学与结构动力学先验知识，选择物理参数的特征尺度；
步骤S3：建立无量纲的空腔流固耦合问题的表征参数系统；
步骤S4：建立空腔流固耦合方程和边界条件；
步骤S5：建立无量纲的空腔流固耦合方程和无量纲的空腔流动边界远场条件。
2.根据权利要求1所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S1具体是指：分析空腔流动与结构振动的基本规律，建立空腔流固耦合问题的关键影响参数系统；所述关键影响参数系统具体包括表征空腔几何参数、表征空腔结构参数和表征来流条件参数；
所述表征空腔几何参数包括空腔长度L、空腔深度D以及空腔宽度W；
所述表征空腔结构参数包括空腔壁板材料的杨氏模量E、空腔壁板材料的泊松比ν、空腔壁板材料的密度ρw以及空腔壁板的厚度h；
所述表征来流条件参数包括来流介质压力ps、来流介质密度ρs、来流介质温度Ts、来流介质粘性μ、来流介质比热比γ、来流介质热传导系数ks、来流介质的定容比热Cv、来流介质的定压比热Cp、来流介质的气体常数R、来流速度Us、来流速度方向的三个方位角αl(l＝1,2,3)、速度边界层厚度δ、位移边界层厚度δ*以及动量边界层厚度θ。
3.根据权利要求2所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S2中表征空腔多物理场特性的基本物理参数包括空间位置坐标Xk、时间变量t、空腔流场速度分量uk、空腔流场温度T、空腔流场密度ρ、空腔流场压力p、空腔壁板驱动力pw、空腔壁板切向的面内坐标α、空腔壁板切向的面内坐标β以及空腔壁板中平面离面位移η；
所述特征尺度包括空腔长度L、流动特征频率fc、来流速度Us、来流温度Ts、来流密度ρs、来流速压ρsUs 2/2、来流介质压力ps以及空腔壁板厚度h；
其中，式中：
下标k(k＝1,2,3)表示物理量在三个空间方向的分量；
下标s表示该物理量对应来流介质的物理参数。
4.根据权利要求3所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S3具体是指：根据空腔流动与结构振动的表征物理参数对应的特征尺度，对表征物理参数进行无量纲化，获得无量纲的空腔流固耦合问题的表征参数系统；所述无量纲的空腔流固耦合问题的表征参数系统包括空间位置坐标的无量纲参数时间变量的无量纲参数空腔流场速度的无量纲参数空腔流场压力的无量纲参数空腔流场密度的无量纲参数空腔流场温度的无量纲参数、空腔壁板驱动力的无量纲参数空腔壁板两个面内坐标的无量纲参数以及空腔壁板中面位移的无量纲参数
5.根据权利要求4所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S4具体是指：根据空腔流固耦合的基本规律，建立空腔流动与结构振动的耦合方程和每个物理场对应的边界条件：具体包括以下步骤：
步骤S41：建立空腔流动控制方程及空腔流动边界条件；
步骤S42：建立空腔壁板结构振动控制方程及空腔壁板结构振动边界条件。
6.根据权利要求5述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S41中空腔流动控制方程为：
其中，式中：e为来流介质内能；
下标k表示物理量在空间三个方向的分量；
下标i、下标j和下标r为哑标；
所述步骤S41中空腔流动边界条件包括空腔流动边界壁面条件和空腔流动边界远场条件；所述空腔流动边界壁面条件的控制方程为：
其中，式中：Dw为空腔壁面区域；
nlk表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
下标k为物理量在三个空间坐标方向的分量；
所述空腔流动边界远场条件的控制方程为：
其中，式中：D∞表示远离空腔壁面的空间区域。
7.根据权利要求6所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S42中具体包括空腔壁板结构振动控制方程为：
其中，式中：为空腔壁板的抗弯刚度；
pw为空腔壁板驱动应力；
ρw为空腔壁板材料的密度；
所述空腔壁板结构振动边界条件包括空腔壁板结构壁面条件和空腔壁板结构壁面交界线条件；
所述空腔壁板结构壁面条件的控制方程为：
其中，式中：下标i为哑标；
nli表示第l个空腔壁板的单位法向量在三个空间坐标方向的分量；
所述空腔壁板结构壁面交界线条件的控制方程为：
其中，式中：Dln表示空腔壁板之间交界线所在空间区域；
Λ表示空腔壁板中平面上垂直于空腔壁板之间交界线的空间方向。
8.根据权利要求7所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S5具体是指：根据有量纲的空腔流固耦合方程与边界条件，依据无量纲的空腔流固耦合问题的表征参数系统，建立空腔流动与结构振动的无量纲方程和无量纲边界条件，获得表征空腔流固耦合问题的无量纲参数之间的关系；具体包括以下步骤：
步骤S51：建立无量纲的空腔流动控制方程及无量纲的空腔流动边界条件；
步骤S52：无量纲的空腔流动边界条件。
9.根据权利要求8所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述步骤S51中的无量纲的空腔流动控制方程具体是指：将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(1)，并且对方程中的项进行化简，得到无量纲的空腔流动控制方程；所述无量纲的空腔流动控制方程为：
其中，式中：
所述步骤S51中无量纲的空腔流动边界条件，具体是指：所述无量纲的空腔流动边界条件包括无量纲的空腔流动边界壁面条件和无量纲的空腔流动边界远场条件；
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(2)得到无量纲的空腔流动边界壁面条件的控制方程，所述无量纲的空腔流动边界壁面条件的控制方程为：
其中，式中：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(3)中得到无量纲的空腔流动边界远场条件，所述无量纲的空腔流动边界远场条件的控制方程为：
10.根据权利要求9所述的一种空腔流固耦合方程的无量纲化方法，其特征在于：所述无量纲的空腔流动控制方程具体是指：步骤S52具体包括空腔壁板结构振动控制方程的无量纲化和空腔壁板结构振动边界条件的无量纲化；
所述空腔壁板结构振动控制方程的无量纲化具体是指将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(4)中，并且对方程中的项进行化简，得到无量纲的空腔壁板结构振动控制方程，所述无量纲的空腔壁板结构振动控制方程为：
其中，式中：
所述空腔壁板结构振动边界条件的无量纲化具体是指：无量纲的空腔壁板结构振动边界条件包括无量纲的空腔壁板结构壁面条件和无量纲的空腔壁板结构壁面交界线条件；
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(5)，并且对方程中的项进行化简，得到无量纲的空腔壁板结构壁面条件的控制方程，所述无量纲的空腔壁板结构壁面条件的控制方程为：
将无量纲的空腔流固耦合问题的表征参数系统所列无量纲参数代入方程(5)，并且对方程中的项进行化简，得到无量纲的空腔壁板结构壁面交界线条件的控制方程，所述无量纲的空腔壁板结构壁面交界线条件的控制方程为：
CN201810084218.8A 2018-01-29 2018-01-29 一种空腔流固耦合方程的无量纲化方法 Active CN108304643B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810084218.8A CN108304643B (zh) 2018-01-29 2018-01-29 一种空腔流固耦合方程的无量纲化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810084218.8A CN108304643B (zh) 2018-01-29 2018-01-29 一种空腔流固耦合方程的无量纲化方法

Publications (2)

Publication Number Publication Date
CN108304643A true CN108304643A (zh) 2018-07-20
CN108304643B CN108304643B (zh) 2021-07-06

Family

ID=62866744

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810084218.8A Active CN108304643B (zh) 2018-01-29 2018-01-29 一种空腔流固耦合方程的无量纲化方法

Country Status (1)

Country Link
CN (1) CN108304643B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112685975A (zh) * 2021-01-05 2021-04-20 中南大学 柔性结构-空腔组合系统流固耦合分析的混合建模方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103439475A (zh) * 2013-08-21 2013-12-11 国家电网公司 具有圆形隧道的准饱和粘弹性土振动响应的检测方法
CN104298869A (zh) * 2014-10-07 2015-01-21 北京理工大学 一种弹性水翼的流固耦合特性数值预测方法
CN105092449A (zh) * 2015-07-14 2015-11-25 山东科技大学 基于水的热流固耦合煤体真三轴剪切渗流实验装置及其实验方法
CN105205269A (zh) * 2015-09-25 2015-12-30 宝鸡石油钢管有限责任公司 一种预测特殊螺纹接头油套管流固耦合特性的分析方法
CN106844914A (zh) * 2017-01-09 2017-06-13 西北工业大学 一种空天飞行器机翼振动响应的快速仿真方法
US10267131B2 (en) * 2012-08-13 2019-04-23 Schlumberger Technology Corporation Competition between transverse and axial hydraulic fractures in horizontal well

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10267131B2 (en) * 2012-08-13 2019-04-23 Schlumberger Technology Corporation Competition between transverse and axial hydraulic fractures in horizontal well
CN103439475A (zh) * 2013-08-21 2013-12-11 国家电网公司 具有圆形隧道的准饱和粘弹性土振动响应的检测方法
CN104298869A (zh) * 2014-10-07 2015-01-21 北京理工大学 一种弹性水翼的流固耦合特性数值预测方法
CN105092449A (zh) * 2015-07-14 2015-11-25 山东科技大学 基于水的热流固耦合煤体真三轴剪切渗流实验装置及其实验方法
CN105205269A (zh) * 2015-09-25 2015-12-30 宝鸡石油钢管有限责任公司 一种预测特殊螺纹接头油套管流固耦合特性的分析方法
CN106844914A (zh) * 2017-01-09 2017-06-13 西北工业大学 一种空天飞行器机翼振动响应的快速仿真方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
S.A.M. MEHRYAN: "Analysis of fluid-solid interaction in MHD natural convection in a square cavity equally partitioned by a vertical flexible membrane", 《JOURNAL OF MAGNETISM AND MAGNETIC MATERIALS》 *
李帅军: "管路系统流固耦合振动无量纲计算方法", 《中国会议》 *
王显圣: "弹性空腔流致噪声／结构振动特性试验", 《航空学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112685975A (zh) * 2021-01-05 2021-04-20 中南大学 柔性结构-空腔组合系统流固耦合分析的混合建模方法

Also Published As

Publication number Publication date
CN108304643B (zh) 2021-07-06

Similar Documents

Publication Publication Date Title
CN107895069A (zh) 一种基于复合材料结构的流固耦合数值预测方法
CN107808065B (zh) 三维复杂外形高速飞行器流-固-热快速计算方法
Jantzen et al. Vortex dynamics around pitching plates
CN110287643A (zh) 空腔流致振动与流致噪声耦合特性的风洞实验模拟方法
Weber et al. 3-D transonic flow in a compressor cascade with shock-induced corner stall
CN111380663B (zh) 一种基于稳定性方法的横流转捩实验数据拓展技术
CN105550467B (zh) 高速风洞斜孔壁壁面流动模型的建立方法
Li et al. On the feedback mechanism in supersonic cavity flows
CN108196443A (zh) 变循环发动机的非线性预测控制设计方法
CN104298886A (zh) 航空发动机旋转部件的结冰3-d数值模拟方法
CN105046021A (zh) 非定常气动力最小状态有理近似的非线性优化算法
CN104899418A (zh) 混压式超声速、高超声速进气道不起动振荡频率预测方法
Tamaki et al. Wall modeling for large-eddy simulation on non-body-conforming Cartesian grids
Liu et al. Pressure gradient effects on the secondary instability of Mack mode disturbances in hypersonic boundary layers
Liang et al. Shock train/glancing shock/boundary layer interaction in a curved isolator with sidewall contraction
CN108304643A (zh) 一种空腔流固耦合方程的无量纲化方法
Elmiligui et al. USM3D Simulations for Third Sonic Boom Workshop
Délery et al. The physics of shock wave/boundary layer interaction control: last lessons learned
CN110104164A (zh) 一种用于跨声速机翼的前加载-吸气组合流动控制方法
Fan et al. Prediction of hypersonic boundary layer transition with variable specific heat on plane flow
CN112985755A (zh) 精确预测腔体流声载荷的附面层相似参数模拟方法
Tang et al. Study of multiple solution phenomenon for hypersonic air inlet
Nianhua et al. Unstructured mesh size control method based on artificial neural network
CN110543677A (zh) 一种涡特征驱动的旋转湍流pans模型
CN110287647A (zh) 一种跨声速压气机平面叶栅激波控制的设计方法

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