CN102004821A - 织物流固耦合模型的动态流场混合分析方法 - Google Patents

织物流固耦合模型的动态流场混合分析方法 Download PDF

Info

Publication number
CN102004821A
CN102004821A CN2010105364400A CN201010536440A CN102004821A CN 102004821 A CN102004821 A CN 102004821A CN 2010105364400 A CN2010105364400 A CN 2010105364400A CN 201010536440 A CN201010536440 A CN 201010536440A CN 102004821 A CN102004821 A CN 102004821A
Authority
CN
China
Prior art keywords
fabric
flow field
dynamic
expression
fluid
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
CN2010105364400A
Other languages
English (en)
Other versions
CN102004821B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN 201010536440 priority Critical patent/CN102004821B/zh
Publication of CN102004821A publication Critical patent/CN102004821A/zh
Application granted granted Critical
Publication of CN102004821B publication Critical patent/CN102004821B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公布了一种织物流固耦合模型的动态流场混合分析方法，方法采用计算步骤为：(1)建立织物三维网格模型；(2)基于LS-DYNA的ALE法对织物进行动态分析，得到织物动态变化情况；(3)将瞬时织物外形作为流场边界条件，并且划分贴体网格；(4)基于FLUENT进行流场分析。本发明首次利用LS-DYNA的ALE法(ArbitraryLagrangian-Eulerian)获得织物动态形状作为FLUENT边界条件分析流场情况，有效解决了ALE法流场描述不精确和柔性织物大变形下动网格流场计算困难的问题。

Description

织物流固耦合模型的动态流场混合分析方法
技术领域
本发明涉及折叠织物流固耦合中流场的动态分析，首次采用LS-DYNA/FLUENT混合分析的方法可以获得织物动态变化时更为精确的流场情况。解决了ALE法流场描述不精确和柔性织物大变形下的动网格流场计算难题。
背景技术
高速流体作用下的织物运动，其流场变化剧烈，必须考虑流场结构的影响，流固耦合方法成为织物工作过程研究的重要手段。由于织物工作过程属于典型的流固耦合、强非线性时变系统，其结构分析属于几何非线性与材料非线性并存的瞬间大变形结构动力学问题，流场研究则属于柔性折叠体大变形下的湍流流动问题，耦合研究非常困难。在国内，已有不少学者对气囊的工作过程进行了研究，取得了一定的成果，但基于流场结构耦合方法对复杂折叠状态下的织物工作过程研究较少，尤其是流场的动态变化规律更是研究中的盲点。LS-DYNA有限元分析软件是一款优良的显式动态分析软件，在许多流固耦合问题上得到了较好的应用，但由于其流场模型的局限性，难以准确描述复杂的流场细节。而FLUENT软件则是优秀的流场分析软件，但是织物在高速流体作用下的运动是个随机过程，很难用函数来描述其各个节点的运动轨迹，因而采用FLUENT的动网格计算较难实现。对于柔性织物流固耦合模型的动态流场分析无论基于哪种软件，都无法克服其局限性。
发明内容
本发明目的是针对现有技术存在的缺陷，解决织物流固耦合动态流场分析问题，提供一种物流固耦合模型的动态流场混合分析方法。
本发明为实现上述目的，采用如下技术方案：
本发明织物流固耦合模型的动态流场混合分析方法，包括以下步骤：
第一步：对织物建立三维模型
使用三角形单元对织物划分网格；如果织物需要折叠，使用直接折叠法或者初始矩阵法完成织物的折叠；动态流场使用六面体网格，织物使用壳结构单元，可以穿插于动态流场网格中；
第二步：采用ALE法对织物进行充气过程计算
对织物建立FSI模型，控制方程如下：
流场控制方程由质量方程、动量方程和能量方程组成，分别为：
Figure 912833DEST_PATH_IMAGE001
，
，
Figure 394630DEST_PATH_IMAGE002
，
式中
Figure 277135DEST_PATH_IMAGE003
表示物质速度；表示相对速度
Figure 103326DEST_PATH_IMAGE005
；
Figure 820222DEST_PATH_IMAGE006
表示网格的速度；
Figure 2010105364400100002DEST_PATH_IMAGE007
表示应力张量
Figure 771997DEST_PATH_IMAGE008
；
Figure 2010105364400100002DEST_PATH_IMAGE009
表示单位体积力；
Figure 775725DEST_PATH_IMAGE010
表示Kronecker
Figure 2010105364400100002DEST_PATH_IMAGE011
-函数；
 网格控制方程为： 
Figure 2010105364400100002DEST_PATH_IMAGE013
表示拉格朗日坐标；
Figure 158613DEST_PATH_IMAGE014
表示欧拉坐标；
Figure 2010105364400100002DEST_PATH_IMAGE015
表示相对速度；
结构控制方程为：
Figure 914080DEST_PATH_IMAGE016
其中M、C、K分别表示单元质量、阻尼模量和弹性模量；F表示膜单元所受合力；
采用显式的松散耦合方法进行计算：在起始步将流固耦合界面的形状作为流场的边界进行流场求解，计算出流固耦合界面上的流体力，然后把求得的流体力当作耦合面上的荷载，计算气囊的结构变形，得到新的流固耦合界面位置；不断重复上述交互过程，直至达到预订求解时间或者问题收敛，得到织物外形动态变化过程；
第三步：导出数据，建立流场网格
导出第二步所述织物外形动态变化过程中的织物充气过程各单元瞬时的外形及运动速度，以及流场入口压力曲线；对各瞬时外形的气囊建立囊内流场贴体网格；囊衣流场边界为无滑移边界条件，即囊衣各单元的运动速度为物面边界条件；充气边界设定为压力入口；
第四步：基于FLUENT进行流场计算
基于FLUENT软件对流场采用大涡模拟或雷诺平均N-S方法进行计算，直至结果收敛。
本分析方法的有益效果是：考虑了流场-结构相互作用的影响，获得了更为准确的动态流场细节，结果可靠。
附图说明
图1气囊初始模型，a.气囊折叠模型，b.折叠仰视图（中心孔为进气口）；
图2 充气口压力曲线；
图3流场网格划分；
图4流场压力云图，左.基于ALE法的LS-DYNA结果，右.混合分析方法结果；a. t=0.02s
时流场压力云图，b. t=0.08s时流场压力云图，c. t=0.16s时流场压力云图；
图5流场速度矢量，左.基于ALE法的LS-DYNA结果，右.混合分析方法结果；a. t=0.02s时内部流场速度矢量，b. t=0.08s时内部流场速度矢量，c. t=0.16s时内部流场速度矢量；
图6气囊展开过程中囊体外形、应力、流场压力、速度，左.外形及应力变化，中.内部流场压力云图，右.内部流场速度矢量；a. t=0.005s，b. t=0.02s，c. t=0.04s，
d. t=0.08s，e. t=0.12s，f. t=0.16s，g. t=0.2s。
具体实施方式
下面以某折叠气囊充气展开过程为例进行流场混合分析
首先，在基于HYPERMESH平台下，建立了LS-DYNA流固耦合分析所需初始充气状态下的结构模型（图1）。气囊厚度1mm，充满气体时长，宽，高分别为1m，1m，0.5m，其基本参数见表1。
采用质量流量为入口边界条件，囊外为定常压力，囊内流场采用无反射边界条件。经过LS-DYNA计算获得了囊体结构动态变化情况。
导出织物充气过程各单元瞬时的外形及运动速度以及充气口压力曲线（图2）。对各瞬时外形的气囊建立囊内流场贴体网格(图3)。假设囊衣流场边界为无滑移边界条件，即囊衣各单元的运动速度为物面边界条件；充气边界设定为压力入口。
基于FLUENT软件，对流场采用标准
Figure DEST_PATH_IMAGE017
两方程模型进行计算，直至结果收敛。
为了验证混合分析方法准确性的，图4，图5为LS-DYNA分析法和混合分析方法所得到的流场结果对比，可以发现两种分析方法所得到的流场基本规律相似，但混合分析方法的压力梯度分明，涡旋更清晰，并且能够捕捉到更细小的流场空间结构，可以更为准确的对充气过程的流场动态特性进行分析。因此采用LS-DYNA/FLUENT混合分析的方法可以获得织物动态变化时更为精细的流场情况（图6）。
 
表1 气囊模型参数
部件 单元类型 弹性模量 泊松比 密度
气囊 BT壳单元 5E8Pa 0.2
Figure 506735DEST_PATH_IMAGE018

Claims (1)

1.一种织物流固耦合模型的动态流场混合分析方法，其特征在于包括以下步骤：
第一步：对织物建立三维模型
使用三角形单元对织物划分网格；如果织物需要折叠，使用直接折叠法或者初始矩阵法完成织物的折叠；动态流场使用六面体网格，织物使用壳结构单元穿插于动态流场网格中；
第二步：采用ALE法对织物进行充气过程计算
对织物建立FSI模型，控制方程如下：
流场控制方程由质量方程、动量方程和能量方程组成，分别为：
Figure 2010105364400100001DEST_PATH_IMAGE002
，
Figure 2010105364400100001DEST_PATH_IMAGE004
，
Figure 2010105364400100001DEST_PATH_IMAGE006
，
式中
Figure DEST_PATH_IMAGE008
表示物质速度；
Figure DEST_PATH_IMAGE010
表示相对速度
Figure DEST_PATH_IMAGE012
；表示网格的速度；
Figure DEST_PATH_IMAGE016
表示应力张量
Figure DEST_PATH_IMAGE018
；表示单位体积力；
Figure DEST_PATH_IMAGE022
表示Kronecker-函数；
 网格控制方程为： 
Figure DEST_PATH_IMAGE026
Figure DEST_PATH_IMAGE028
表示拉格朗日坐标；
Figure DEST_PATH_IMAGE030
表示欧拉坐标；表示相对速度；
结构控制方程为：
Figure DEST_PATH_IMAGE032
其中M、C、K分别表示单元质量、阻尼模量和弹性模量；F表示膜单元所受合力；
采用显式的松散耦合方法进行计算：在起始步将流固耦合界面的形状作为流场的边界进行流场求解，计算出流固耦合界面上的流体力，然后把求得的流体力当作耦合面上的荷载，计算气囊的结构变形，得到新的流固耦合界面位置；不断重复上述交互过程，直至达到预订求解时间或者问题收敛，得到织物外形动态变化过程；
第三步：导出数据，建立流场网格
导出第二步所述织物外形动态变化过程中的织物充气过程各单元瞬时的外形及运动速度，以及流场入口压力曲线；对各瞬时外形的气囊建立囊内流场贴体网格；囊衣流场边界为无滑移边界条件，即囊衣各单元的运动速度为物面边界条件；充气边界设定为压力入口；
第四步：基于FLUENT进行流场计算
基于FLUENT软件对流场采用大涡模拟或雷诺平均N-S方法进行计算，直至结果收敛。
CN 201010536440 2010-11-09 2010-11-09 织物流固耦合模型的动态流场混合分析方法 Expired - Fee Related CN102004821B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010536440 CN102004821B (zh) 2010-11-09 2010-11-09 织物流固耦合模型的动态流场混合分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010536440 CN102004821B (zh) 2010-11-09 2010-11-09 织物流固耦合模型的动态流场混合分析方法

Publications (2)

Publication Number Publication Date
CN102004821A true CN102004821A (zh) 2011-04-06
CN102004821B CN102004821B (zh) 2013-07-24

Family

ID=43812181

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010536440 Expired - Fee Related CN102004821B (zh) 2010-11-09 2010-11-09 织物流固耦合模型的动态流场混合分析方法

Country Status (1)

Country Link
CN (1) CN102004821B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103559346A (zh) * 2013-10-30 2014-02-05 曙光信息产业(北京)有限公司 Ls-dyna的计算任务的实现方法和装置
CN104317985A (zh) * 2014-09-19 2015-01-28 大连理工大学 一种基于界带有限元和拉格朗日坐标的流体仿真方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101236573A (zh) * 2007-05-25 2008-08-06 清华大学 一种柔性翼微型飞行器的流固耦合数值仿真方法
JP2008242516A (ja) * 2007-03-23 2008-10-09 Aichi Prefecture 織物の３次元モデル構築方法及び織物の３次元モデル構築装置
CN101727541A (zh) * 2009-12-07 2010-06-09 南京航空航天大学 伞衣织物透气性计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008242516A (ja) * 2007-03-23 2008-10-09 Aichi Prefecture 織物の３次元モデル構築方法及び織物の３次元モデル構築装置
CN101236573A (zh) * 2007-05-25 2008-08-06 清华大学 一种柔性翼微型飞行器的流固耦合数值仿真方法
CN101727541A (zh) * 2009-12-07 2010-06-09 南京航空航天大学 伞衣织物透气性计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
余莉等: "气囊充气过程流固耦合数值模拟", 《南京航空航天大学学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103559346A (zh) * 2013-10-30 2014-02-05 曙光信息产业(北京)有限公司 Ls-dyna的计算任务的实现方法和装置
CN103559346B (zh) * 2013-10-30 2016-10-05 曙光信息产业(北京)有限公司 Ls-dyna的计算任务的实现方法和装置
CN104317985A (zh) * 2014-09-19 2015-01-28 大连理工大学 一种基于界带有限元和拉格朗日坐标的流体仿真方法
CN104317985B (zh) * 2014-09-19 2017-05-24 大连理工大学 一种基于界带有限元和拉格朗日坐标的流体仿真方法

Also Published As

Publication number Publication date
CN102004821B (zh) 2013-07-24

Similar Documents

Publication Publication Date Title
Lian et al. Membrane wing aerodynamics for micro air vehicles
CN102799730B (zh) 一种燃气轮机风扇叶片反扭过程的预估方法
CN102867094B (zh) 一种移动粒子半隐式算法中自由表面流动模型的构建方法
CN102156785A (zh) 一种动态布料仿真方法
CN103559335B (zh) 一种基于流固耦合的伞降仿真方法
CN104133933A (zh) 一种高超声速飞行器热环境下气动弹性力学特性分析方法
CN103400035B (zh) 一种高可信度快速预测飞行器滚转动导数的方法
CN103544347A (zh) 一种充气展开织物工作过程的有限元分析方法
CN102682146B (zh) 直升机旋翼可压缩旋流场的数值模拟方法
CN103389649B (zh) 一种基于球面拼接网格的飞行器机动运动模拟方法
CN107423511A (zh) 满足无滑边界条件和连续性条件浸入边界隐式迭代求解法
CN109408836A (zh) 利用Boltzmann方程进行流体仿真的方法
Esmaeili et al. Active feedback VIV control of sprung circular cylinder using TDE-iPID control strategy at moderate Reynolds numbers
CN102004821B (zh) 织物流固耦合模型的动态流场混合分析方法
CN101727541A (zh) 伞衣织物透气性计算方法
CN106508020B (zh) 一种可用于工程设计的复杂飞行器气动加热计算方法
CN103544053A (zh) 一种考虑织物透气性的降落伞气动性能分析方法
CN103914602A (zh) 可压缩旋流场的数值模拟方法
CN106773782A (zh) 一种气动伺服弹性混合建模方法
US9390205B2 (en) Vorticity-refinement based numerical method for simulating aircraft wing-tip vortex flows
CN105701308A (zh) 一种基于无网格流固耦合的降落伞分析方法
Hu Simulation and analysis of airbag parameter variation under different gas flows
Bekhti et al. Flap/lead-lag computational investigations on NREL S809 airfoil
CN112926252A (zh) 液体与布料交互现象模拟方法
CN102004820B (zh) 基于fsi的复杂折叠织物逆向建模方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130724

Termination date: 20141109

EXPY Termination of patent right or utility model