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

技术编号:5240825 阅读:357 留言:0更新日期:2012-04-11 18:40
本发明专利技术公布了一种织物流固耦合模型的动态流场混合分析方法,方法采用计算步骤为:(1)建立织物三维网格模型;(2)基于LS-DYNA的ALE法对织物进行动态分析,得到织物动态变化情况;(3)将瞬时织物外形作为流场边界条件,并且划分贴体网格;(4)基于FLUENT进行流场分析。本发明专利技术首次利用LS-DYNA的ALE法(ArbitraryLagrangian-Eulerian)获得织物动态形状作为FLUENT边界条件分析流场情况,有效解决了ALE法流场描述不精确和柔性织物大变形下动网格流场计算困难的问题。

【技术实现步骤摘要】

本专利技术涉及折叠织物流固耦合中流场的动态分析,首次采用LS-DYNA/FLUENT混 合分析的方法可以获得织物动态变化时更为精确的流场情况。解决了 ALE法流场描述不精 确和柔性织物大变形下的动网格流场计算难题。
技术介绍
高速流体作用下的织物运动,其流场变化剧烈,必须考虑流场结构的影响,流固耦 合方法成为织物工作过程研究的重要手段。由于织物工作过程属于典型的流固耦合、强非 线性时变系统,其结构分析属于几何非线性与材料非线性并存的瞬间大变形结构动力学问 题,流场研究则属于柔性折叠体大变形下的湍流流动问题,耦合研究非常困难。在国内,已 有不少学者对气囊的工作过程进行了研究,取得了一定的成果,但基于流场结构耦合方法 对复杂折叠状态下的织物工作过程研究较少,尤其是流场的动态变化规律更是研究中的盲 点。LS-DYNA有限元分析软件是一款优良的显式动态分析软件,在许多流固耦合问题上得到 了较好的应用,但由于其流场模型的局限性,难以准确描述复杂的流场细节。而FLUENT软 件则是优秀的流场分析软件,但是织物在高速流体作用下的运动是个随机过程,很难用函 数来描述其各个节点的运动轨迹,因而采用FLUENT的动网格计算较难实现。对于柔性织物 流固耦合模型的动态流场分析无论基于哪种软件,都无法克服其局限性。
技术实现思路
本专利技术目的是针对现有技术存在的缺陷,解决织物流固耦合动态流场分析问题, 提供一种物流固耦合模型的动态流场混合分析方法。本专利技术为实现上述目的,采用如下技术方案本专利技术,包括以下步骤 第一步对织物建立三维模型使用三角形单元对织物划分网格;如果织物需要折叠,使用直接折叠法或者初始矩阵 法完成织物的折叠;动态流场使用六面体网格,织物使用壳结构单元,可以穿插于动态流场 网格中;第二步采用ALE法对织物进行充气过程计算对织物建立FSI模型,控制方程如下流场控制方程由质量方程、动量方程和能量方程组成,分别为权利要求1. 一种,其特征在于包括以下步骤 第一步对织物建立三维模型使用三角形单元对织物划分网格;如果织物需要折叠,使用直接折叠法或者初始矩阵 法完成织物的折叠;动态流场使用六面体网格,织物使用壳结构单元穿插于动态流场网格 中;第二步采用ALE法对织物进行充气过程计算对织物建立FSI模型,控制方程如下流场控制方程由质量方程、动量方程和能量方程组成,分别为dp dv= dp-f- = —P~--Mil- !dt Qsci Qx^, Qyl,^ , QE QEi^uphA-pwS,式中巧表示物质速度;Wf表示相对速度1 |=兮-吟;吟表示网格的速度Oe表示应力张丨ισβ = -P5i + M巧J + Vjj) ; ~表示单位体积力;各表示Kronecker S ~函数; 网格控制方程为m m 1 m為表示拉格朗日坐标A表示欧拉坐标;Wf表示相对速度; 结构控制方程为Aiw^Cw^Kw = P其中Μ、C、K分别表示单元质量、阻尼模量和弹性模量;F表示膜单元所受合力; 采用显式的松散耦合方法进行计算在起始步将流固耦合界面的形状作为流场的边界 进行流场求解,计算出流固耦合界面上的流体力,然后把求得的流体力当作耦合面上的荷 载,计算气囊的结构变形,得到新的流固耦合界面位置;不断重复上述交互过程,直至达到 预订求解时间或者问题收敛,得到织物外形动态变化过程; 第三步导出数据,建立流场网格导出第二步所述织物外形动态变化过程中的织物充气过程各单元瞬时的外形及运动 速度,以及流场入口压力曲线;对各瞬时外形的气囊建立囊内流场贴体网格;囊衣流场边 界为无滑移边界条件,即囊衣各单元的运动速度为物面边界条件;充气边界设定为压力入 Π ;第四步基于FLUENT进行流场计算基于FLUENT软件对流场采用大涡模拟或雷诺平均N-S方法进行计算,直至结果收敛。全文摘要本专利技术公布了一种,方法采用计算步骤为(1)建立织物三维网格模型;(2)基于LS-DYNA的ALE法对织物进行动态分析,得到织物动态变化情况;(3)将瞬时织物外形作为流场边界条件,并且划分贴体网格;(4)基于FLUENT进行流场分析。本专利技术首次利用LS-DYNA的ALE法(ArbitraryLagrangian-Eulerian)获得织物动态形状作为FLUENT边界条件分析流场情况,有效解决了ALE法流场描述不精确和柔性织物大变形下动网格流场计算困难的问题。文档编号G06F17/50GK102004821SQ201010536440公开日2011年4月6日 申请日期2010年11月9日 优先权日2010年11月9日专利技术者余莉, 刘雄, 呼政魁, 张绳, 程涵 申请人:南京航空航天大学本文档来自技高网...

【技术保护点】
一种织物流固耦合模型的动态流场混合分析方法,其特征在于包括以下步骤:第一步:对织物建立三维模型使用三角形单元对织物划分网格;如果织物需要折叠,使用直接折叠法或者初始矩阵法完成织物的折叠;动态流场使用六面体网格,织物使用壳结构单元穿插于动态流场网格中;第二步:采用ALE法对织物进行充气过程计算对织物建立FSI模型,控制方程如下:流场控制方程由质量方程、动量方程和能量方程组成,分别为:***式中v↓[i]表示物质速度;w↓[i]表示相对速度w↓[i]=v↓[i]-u↓[i];u↓[i]表示网格的速度;σ↓[ij]表示应力张量σ↓[ij]=-pδ↓[*]+μ(v↓[ij]+v↓[ji]);b↓[i]表示单位体积力;δ↓[ij]表示Kroneckerδ-函数;网格控制方程为:***X↓[i]表示拉格朗日坐标;x↓[i]表示欧拉坐标;w↓[i]表示相对速度;结构控制方程为:M*+C*+Kw=F其中M、C、K分别表示单元质量、阻尼模量和弹性模量;F表示膜单元所受合力;采用显式的松散耦合方法进行计算:在起始步将流固耦合界面的形状作为流场的边界进行流场求解,计算出流固耦合界面上的流体力,然后把求得的流体力当作耦合面上的荷载,计算气囊的结构变形,得到新的流固耦合界面位置;不断重复上述交互过程,直至达到预订求解时间或者问题收敛,得到织物外形动态变化过程;第三步:导出数据,建立流场网格导出第二步所述织物外形动态变化过程中的织物充气过程各单元瞬时的外形及运动速度,以及流场入口压力曲线;对各瞬时外形的气囊建立囊内流场贴体网格;囊衣流场边界为无滑移边界条件,即囊衣各单元的运动速度为物面边界条件;充气边界设定为压力入口;第四步:基于FLUENT进行流场计算基于FLUENT软件对流场采用大涡模拟或雷诺平均N-S方法进行计算,直至结果收敛。...

【技术特征摘要】

【专利技术属性】
技术研发人员:余莉程涵刘雄张绳呼政魁
申请(专利权)人:南京航空航天大学
类型:发明
国别省市:84[中国|南京]

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1