CN104850689A - 一种基于固定网格技术的流固耦合计算方法 - Google Patents
一种基于固定网格技术的流固耦合计算方法 Download PDFInfo
- Publication number
- CN104850689A CN104850689A CN201510217427.1A CN201510217427A CN104850689A CN 104850689 A CN104850689 A CN 104850689A CN 201510217427 A CN201510217427 A CN 201510217427A CN 104850689 A CN104850689 A CN 104850689A
- Authority
- CN
- China
- Prior art keywords
- solid
- fluid
- node
- flow field
- file
- 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
Links
- 239000007787 solid Substances 0.000 title claims abstract description 136
- 238000005516 engineering process Methods 0.000 title claims abstract description 18
- 238000010168 coupling process Methods 0.000 title abstract description 7
- 230000008878 coupling Effects 0.000 title abstract description 6
- 238000005859 coupling reaction Methods 0.000 title abstract description 6
- 238000004364 calculation method Methods 0.000 title abstract description 4
- 239000012530 fluid Substances 0.000 claims abstract description 61
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000006073 displacement reaction Methods 0.000 claims abstract description 28
- 238000004088 simulation Methods 0.000 claims abstract description 15
- 239000013598 vector Substances 0.000 claims description 23
- 239000011159 matrix material Substances 0.000 claims description 19
- 230000001133 acceleration Effects 0.000 claims description 8
- 238000013016 damping Methods 0.000 claims description 4
- 238000007654 immersion Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 238000012805 post-processing Methods 0.000 abstract 1
- 230000008569 process Effects 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000010429 evolutionary process Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于固定网格技术的流固耦合计算方法,属于流固耦合模拟技术领域。本发明包括步骤:A、调用网格划分模块,采用两套网格:流体和固体组成的整个物理区域,固体区域;B、调用流场计算CFD模块,更新流场状态变量,同时获得流体作用于固体边界节点上的力;C、调用固体计算CSD模块,更新固体状态变量,同时获取固体边界节点上的位移和速度;D、调用结果输出模块,将固体状态变量和流场信息输出到文件,供后处理使用;E、判断是否结束计算。本发明成功避免使用动网格模块,大大节省了网格更新所需要的计算资源;克服了传统的基于移动网格技术时复杂的内部网格更新与边界网格控制运算,不容易准确计算流固界面作用力的缺点。
Description
技术领域
本发明涉及一种基于固定网格技术的流固耦合计算方法,属于流固耦合模拟技术领域。
背景技术
流固耦合问题的难点之一是流体和固体使用不同的数学描述框架。通常,流体运动使用欧拉描述,而固体运动使用拉格朗日描述:传统的流固耦合方法,流体区域的网格随固体边界的运动而动态更新,网格运算和网格质量控制过程异常复杂,得到的边界层附近的速度分布误差较大。为此,本发明提出基于固定网格技术的流固耦合计算方法,避免了使用动网格计算模块,计算过程简单,易于程序设计。
发明内容
本发明提供了一种基于固定网格技术的流固耦合计算方法,以用于克服传统基于移动网格技术在计算固体与流体耦合作用上的不足。
本发明的技术方案是:一种基于固定网格技术的流固耦合计算方法,所述方法的具体步骤如下:
Step1、调用网格划分模块,采用两套网格:流场区域,固体区域;其中流场区域在欧拉描述下采用笛卡尔网格离散,固体区域在拉格朗日描述下使用适体有限元网格离散;
Step2、调用流场CFD(计算流体动力学)计算模块,采用基于分步投影的浸入边界方法,求解不可压缩粘性牛顿流体的流动控制方程,更新流场状态变量,同时获得流体作用于固体边界节点上的力;
Step3、调用固体CSD(计算固体动力学)计算模块,采用有限单元方法,求解固体动力学控制方程,更新固体状态变量,同时获取固体边界节点上的位移和速度;
Step4、将固体状态变量和流场信息输出到文件,供后处理使用;
Step5、判断是否结束计算:
如果△t·n<T,则进入下一时间步,继续执行步骤Step2、Step3和Step4;
如果△t·n≥T,则结束整个计算;
其中△t为时间步长,T为要求计算的总物理时间,n为时间步数。
所述流场区域包括流体和固体所占据的空间区域,固体区域可包括刚体区域和弹性体区域。
所述步骤Step1中,流场区域在欧拉描述下采用笛卡尔网格划分,其网格单元中心或节点上的流场变量称为欧拉变量,并将网格节点坐标信息xj输出到文件fcor.txt;固体区域在拉格朗日描述下使用适体有限元网格划分,其网格节点上的变量称为拉格朗日变量,相应网格节点坐标信息X输出到文件scor.txt,其中固体边界节点坐标信息可用表示。
所述步骤Step1中,流场区域在欧拉描述下采用笛卡尔网格划分,其网格单元中心或节点上的流场变量称为欧拉变量,并将网格节点坐标信息xj输出到文件fcor.txt;固体区域在拉格朗日描述下使用适体有限元网格划分,其网格节点上的变量称为拉格朗日变量,相应网格节点坐标信息X输出到文件scor.txt,其中固体边界节点坐标信息可用表示。
所述步骤Step2中,通过基于分步投影浸入边界方法的CFD计算,实现两大功能:
一是获得流体作用于固体边界节点上的力:
通过流场区域预测速度u′(xj,t)和近似光滑函数得到的固体边界拉格朗日节点上的速度应该等于给定的固体边界节点的自然速度来实现固体边界节点上力密度的求解,并将结果输出到文件sfor.txt,供步骤Step3中CSD计算模块使用;其中,为固体边界节点坐标的解释,t为时间;
二是更新流场状态变量:
流场区域速度校正值Δu(xj,t)为:
式中,Δsi为第i段固体边界的面积,h为流体网格间距,M为固体边界节点数,N为流体欧拉网格节点数,Cji为信息转换矩阵,定义如下:
式中,分别为固体边界节点坐标在x,y,z方向的分量,xj,yj,zj分别为欧拉网格节点坐标xj在x,y,z方向的分量;函数φ可表示为:
式中,r为函数φ的自变量;
流场区域速度可由下式更新,
u(xj,t)=u′(xj,t)+Δu(xj,t) (4)
并将结果输出到文件fvel.txt,供步骤Step4中的后处理使用。
所述步骤Step3中,通过有限单元方法,求解离散后的固体动力学控制方程,实现两大功能,更新固体状态变量,同时获取固体边界节点上的位移和速度;
一是获取固体边界节点上的位移和速度:
采用有限单元数值离散方法,可获得固体的一组离散非线性方程组,写成矩阵形式有:
式中,Ms为质量矩阵,Cs阻尼矩阵,Ks为刚度矩阵,Qs为外力矢量,包括流体作用于固体边界的力密度和其他外力矢量,和分别代表固体t时刻的加速度矢量,速度矢量和位移矢量;
式(5)求解在时间推进上采用纽马克法,可获得t+Δt时刻有限元节点新的位移矢量速度矢量和加速度矢量提取t+Δt时刻边界节点上的位移矢量速度矢量令并将结果输出到文件sdisvel.txt文件,供步骤Step2中CFD计算模块使用;
二是更新固体状态变量:
通过对式(5)的求解,获得下一时刻有限元节点新的位移矢量通过固体本构方程,可进一步获得固体各节点的应力应变等固体其他状态变量的值,并将结果输出到文件sstress.txt,供步骤Step4中的后处理使用。
本发明的有益效果是:
1、避免使用动网格技术,大量节省计算资源:传统的基于移动网格技术的流固耦合方法,需要借助动网格技术,而对于复杂的流固耦合问题,大幅度的固体运动往往导致流体区域网格更新的失败,而本发明正是弥补这一重要缺陷,在固体与流体耦合作用过程中成功避免使用动网格技术。
2、即使近年来也有少量成功应用动网格技术来执行流固耦合计算,但也很难保证耦合界面信息交换时能量守恒和边界层网格质量有效控制等关键问题,大大降低了模拟精度。本发明采用适宜的近似光滑函数,利用流固界面一致条件,确保流体与固体界面作用力平衡,保证耦合计算的有效性。
3、本发明使用的计算流程物理概念清晰,计算过程简单,易于程序设计,便于处理复杂的流固耦合问题。
附图说明
图1为本发明中的流程图;
图2为本发明中算例几何尺寸和整个系统计算区域的示意图;
图3为本发明中计算得到的悬臂梁自由端的位移随时间的演化过程。
具体实施方式
实施例1:如图1-3所示,一种基于固定网格技术的流固耦合计算方法,所述方法的具体步骤如下:
Step1、调用网格划分模块,采用两套网格:流场区域,固体区域;其中流场区域在欧拉描述下采用笛卡尔网格离散,固体区域在拉格朗日描述下使用适体有限元网格离散;
Step2、调用流场CFD(计算流体动力学)计算模块,采用基于分步投影的浸入边界方法,求解不可压缩粘性牛顿流体的流动控制方程,更新流场状态变量,同时获得流体作用于固体边界节点上的力;
Step3、调用固体CSD(计算固体动力学)计算模块,采用有限单元方法,求解固体动力学控制方程,更新固体状态变量,同时获取固体边界节点上的位移和速度;
Step4、将固体状态变量和流场信息输出到文件,供后处理使用;
Step5、判断是否结束计算:
如果△t·n<T,则进入下一时间步,继续执行步骤Step2、Step3和Step4;
如果△t·n≥T,则结束整个计算;
其中△t为时间步长,T为要求计算的总物理时间,n为时间步数。
所述流场区域包括流体和固体所占据的空间区域,固体区域可包括刚体区域和弹性体区域。
所述步骤Step1中,流场区域在欧拉描述下采用笛卡尔网格划分,其网格单元中心或节点上的流场变量称为欧拉变量,并将网格节点坐标信息xj输出到文件fcor.txt;固体区域在拉格朗日描述下使用适体有限元网格划分,其网格节点上的变量称为拉格朗日变量,相应网格节点坐标信息X输出到文件scor.txt,其中固体边界节点坐标信息可用表示。
所述步骤Step1中,流场区域在欧拉描述下采用笛卡尔网格划分,其网格单元中心或节点上的流场变量称为欧拉变量,并将网格节点坐标信息xj输出到文件fcor.txt;固体区域在拉格朗日描述下使用适体有限元网格划分,其网格节点上的变量称为拉格朗日变量,相应网格节点坐标信息X输出到文件scor.txt,其中固体边界节点坐标信息可用表示。
所述步骤Step2中,通过基于分步投影浸入边界方法的CFD计算,实现两大功能:
一是获得流体作用于固体边界节点上的力:
通过流场区域预测速度u′(xj,t)和近似光滑函数得到的固体边界拉格朗日节点上的速度应该等于给定的固体边界节点的自然速度来实现固体边界节点上力密度的求解,并将结果输出到文件sfor.txt,供步骤Step3中CSD计算模块使用;其中,为固体边界节点坐标的解释,t为时间;
二是更新流场状态变量:
流场区域速度校正值Δu(xj,t)为:
式中,Δsi为第i段固体边界的面积,h为流体网格间距,M为固体边界节点数,N为流体欧拉网格节点数,Cji为信息转换矩阵,定义如下:
式中,分别为固体边界节点坐标在x,y,z方向的分量,xj,yj,zj分别为欧拉网格节点坐标xj在x,y,z方向的分量;函数φ可表示为:
式中,r为函数φ的自变量;
流场区域速度可由下式更新,
u(xj,t)=u′(xj,t)+Δu(xj,t) (4)
并将结果输出到文件fvel.txt,供步骤Step4中的后处理使用。
所述步骤Step3中,通过有限单元方法,求解离散后的固体动力学控制方程,实现两大功能,更新固体状态变量,同时获取固体边界节点上的位移和速度;
一是获取固体边界节点上的位移和速度:
采用有限单元数值离散方法,可获得固体的一组离散非线性方程组,写成矩阵形式有:
式中,Ms为质量矩阵,Cs阻尼矩阵,Ks为刚度矩阵,Qs为外力矢量,包括流体作用于固体边界的力密度和其他外力矢量,和分别代表固体t时刻的加速度矢量,速度矢量和位移矢量;
式(5)求解在时间推进上采用纽马克法,可获得t+Δt时刻有限元节点新的位移矢量速度矢量和加速度矢量提取t+Δt时刻边界节点上的位移矢量速度矢量令并将结果输出到文件sdisvel.txt文件,供步骤Step2中CFD计算模块使用;
二是更新固体状态变量:
通过对式(5)的求解,获得下一时刻有限元节点新的位移矢量通过固体本构方程,可进一步获得固体各节点的应力应变等固体其他状态变量的值,并将结果输出到文件sstress.txt,供步骤Step4中的后处理使用。
实施例2:如图1-3所示,一种基于固定网格技术的流固耦合计算方法,所述方法的具体步骤如下:
一静止刚性方柱浸没在不可压缩流场中,左边进口流动假设为均匀流场,速度大小为vx=v0,vy=0。方柱正后方中心附着一个弹性悬臂梁(固体密度ρs=100Kg/m3,弹性模量E=25Gpa,泊松比ν=0.35)。在静止时梁轴线平行于远场流动方向,进口速度v0=51.3cm/s,雷诺数Re=ρfDv0/μf=333,式中D=1cm为刚性方柱的边长,柔性悬臂梁高度为0.06D。
S1:网格划分
流场区域(包括流体和固体区域)在欧拉描述下采用笛卡尔网格划分,网格为均匀四边形网格,网格间距为h=0.01D,并将单元坐标信息输出到文件fcor.txt,固体区域(此算例包括刚体区域和弹性体区域,即刚性方柱和弹性悬臂梁),在拉格朗日描述下使用有限元四边形网格划分,间距为0.01D,相应网格节点信息文件输出到scor.txt文件。
S2:调用流场CFD计算模块,更新流场状态变量,同时获得流体作用于固体边界节点上的力;
通过流场区域预测速度u′(xj,t)和近似光滑函数得到的固体边界拉格朗日节点上的速度应该等于给定的固体边界节点的自然速度。对刚性方柱,其固体边界节点的自然速度对弹性悬臂梁边界节点的自然速度即等于固体CSD计算得到的固体边界的速度,在初始时刻,可取从而实现固体边界节点上力密度的求解,并将结果输出到文件sfor.txt,供下一步骤S3中CSD计算模块使用。同时,更新流场速度,并将结果输出到文件fvel.txt,供步骤S4中的后处理使用。
S3:调用固体CSD计算模块,更新固体状态变量,同时获取固体边界节点上的位移和速度
采用有限单元数值离散方法,可获得固体的一组离散非线性方程组,写成矩阵形式有:式中,Ms为质量矩阵,Cs阻尼矩阵,Ks为刚度矩阵,Qs为外力矢量,包括流体作用于固体边界的力密度和其他外力矢量,和分别代表固体t时刻的加速度矢量,速度矢量和位移矢量。
在时间推进上采用纽马克法,可获得t+Δt时刻有限元节点新的位移矢量速度矢量和加速度矢量提取t+Δt时刻边界节点上的位移矢量速度矢量令并将结果输出到文件sdisvel.txt文件,供步骤S2中CFD计算模块使用。
S4:调用结果输出模块,将作用在固体边界上的力以及流场信息输出到文件,供后处理软件读取显示。
S5:时间推进
在一个时间步内计算完成后,转入下一个时间步,重复上述步骤S2-S4,直到满足计算时间要求停止计算。本例取时间步长Δt=0.001s,计算总时间T=12s,共计时间步n=12000。
图3为悬臂梁自由端的位移随时间的演化过程。从图可见,弹性悬臂梁在流体动力作用下经历了一个长期的不稳定振动过程,约8s后振动趋于稳定的周期振动,梁自由端的振幅约1.4cm左右。
上面结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (5)
1.一种基于固定网格技术的流固耦合计算方法,其特征在于:所述方法的具体步骤如下:
Step1、调用网格划分模块,采用两套网格:流场区域,固体区域;其中流场区域在欧拉描述下采用笛卡尔网格离散,固体区域在拉格朗日描述下使用适体有限元网格离散;
Step2、调用流场CFD计算模块,采用基于分步投影的浸入边界方法,求解不可压缩粘性牛顿流体的流动控制方程,更新流场状态变量,同时获得流体作用于固体边界节点上的力;
Step3、调用固体CSD计算模块,采用有限单元方法,求解固体动力学控制方程,更新固体状态变量,同时获取固体边界节点上的位移和速度;
Step4、将固体状态变量和流场信息输出到文件,供后处理使用;
Step5、判断是否结束计算:
如果△t·n<T,则进入下一时间步,继续执行步骤Step2、Step3和Step4;
如果△t·n≥T,则结束整个计算;
其中△t为时间步长,T为要求计算的总物理时间,n为时间步数。
2.根据权利要求1所述的基于固定网格技术的流固耦合计算方法,其特征在于:所述流场区域包括流体和固体所占据的空间区域,固体区域可包括刚体区域和弹性体区域。
3.根据权利要求1所述的基于固定网格技术的流固耦合计算方法,其特征在于:所述步骤Step1中,流场区域在欧拉描述下采用笛卡尔网格划分,其网格单元中心或节点上的流场变量称为欧拉变量,并将网格节点坐标信息xj输出到文件fcor.txt;固体区域在拉格朗日描述下使用适体有限元网格划分,其网格节点上的变量称为拉格朗日变量,相应网格节点坐标信息X输出到文件scor.txt,其中固体边界节点坐标信息可用表示。
4.根据权利要求1所述的基于固定网格技术的流固耦合计算方法,其特征在于:所述步骤Step2中,通过基于分步投影浸入边界方法的CFD计算,实现两大功能:
一是获得流体作用于固体边界节点上的力:
通过流场区域预测速度u′(xj,t)和近似光滑函数得到的固体边界拉格朗日节点上的速度应该等于给定的固体边界节点的自然速度来实现固体边界节点上力密度的求解,并将结果输出到文件sfor.txt,供步骤Step3中CSD计算模块使用;其中,为固体边界节点坐标的解释,t为时间;
二是更新流场状态变量:
流场区域速度校正值Δu(xj,t)为:
式中,Δsi为第i段固体边界的面积,h为流体网格间距,M为固体边界节点数,N为流体欧拉网格节点数,Cji为信息转换矩阵,定义如下:
式中,分别为固体边界节点坐标在x,y,z方向的分量,xj,yj,zj分别为欧拉网格节点坐标xj在x,y,z方向的分量;函数φ可表示为:
式中,r为函数φ的自变量;
流场区域速度可由下式更新,
u(xj,t)=u′(xj,t)+Δu(xj,t) (4)
并将结果输出到文件fvel.txt,供步骤Step4中的后处理使用。
5.根据权利要求1所述的基于固定网格技术的流固耦合计算方法,其特征在于:所述步骤Step3中,通过有限单元方法,求解离散后的固体动力学控制方程,实现两大功能,更新固体状态变量,同时获取固体边界节点上的位移和速度;
一是获取固体边界节点上的位移和速度:
采用有限单元数值离散方法,可获得固体的一组离散非线性方程组,写成矩阵形式有:
式中,Ms为质量矩阵,Cs阻尼矩阵,Ks为刚度矩阵,Qs为外力矢量,包括流体作用于固体边界的力密度和其他外力矢量,和分别代表固体t时刻的加速度矢量,速度矢量和位移矢量;
式(5)求解在时间推进上采用纽马克法,可获得t+Δt时刻有限元节点新的位移矢量速度矢量和加速度矢量提取t+Δt时刻边界节点上的位移矢量速度矢量令并将结果输出到文件sdisvel.txt文件,供步骤Step2中CFD计算模块使用;
二是更新固体状态变量:
通过对式(5)的求解,获得下一时刻有限元节点新的位移矢量通过固体本构方程,可进一步获得固体各节点的应力应变等固体其他状态变量的值,并将结果输出到文件sstress.txt,供步骤Step4中的后处理使用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510217427.1A CN104850689B (zh) | 2015-04-30 | 2015-04-30 | 一种基于固定网格技术的流固耦合计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510217427.1A CN104850689B (zh) | 2015-04-30 | 2015-04-30 | 一种基于固定网格技术的流固耦合计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104850689A true CN104850689A (zh) | 2015-08-19 |
CN104850689B CN104850689B (zh) | 2019-01-04 |
Family
ID=53850332
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510217427.1A Active CN104850689B (zh) | 2015-04-30 | 2015-04-30 | 一种基于固定网格技术的流固耦合计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104850689B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105760602A (zh) * | 2015-12-30 | 2016-07-13 | 南京航空航天大学 | 一种有限体积加权基本无振荡格式的全流场数值模拟方法 |
CN106126797A (zh) * | 2016-06-20 | 2016-11-16 | 华中科技大学 | 一种水下探测装置无动力运动预报自动计算方法 |
CN106202782A (zh) * | 2016-07-20 | 2016-12-07 | 昆明理工大学 | 一种水轮机导叶主动旋转输入功率的数值计算方法 |
CN106407600A (zh) * | 2016-10-18 | 2017-02-15 | 南京航空航天大学 | 一种基于流‑固耦合及代理模型的入水仿真分析方法 |
CN107895093A (zh) * | 2017-12-08 | 2018-04-10 | 兰州理工大学 | 一种风力机叶片流固耦合模态设计方法 |
CN108038329A (zh) * | 2017-12-25 | 2018-05-15 | 大连理工大学 | 一种微胶囊自愈合材料损伤和愈合机理的数值仿真方法 |
CN108763741A (zh) * | 2018-05-28 | 2018-11-06 | 青岛科技大学 | 一种液压软管流固耦合数值预测方法 |
CN110633556A (zh) * | 2019-10-24 | 2019-12-31 | 南京航空航天大学 | 一种陶瓷基复合材料流固耦合响应计算方法 |
CN110717285A (zh) * | 2019-09-11 | 2020-01-21 | 哈尔滨工程大学 | 一种大幅六自由度运动的流固耦合模拟方法 |
CN111062154A (zh) * | 2019-11-18 | 2020-04-24 | 华北水利水电大学 | 一种南水北调渠系建筑物流固耦合计算方法 |
CN111950174A (zh) * | 2020-07-08 | 2020-11-17 | 上海交通大学 | 流固耦合的计算方法、装置及电子设备 |
CN114757082A (zh) * | 2022-03-07 | 2022-07-15 | 同济大学 | 一种基于拉格朗日-欧拉稳定配点的流固耦合计算方法 |
CN114840889A (zh) * | 2022-07-04 | 2022-08-02 | 海南浙江大学研究院 | 一种网箱网衣受力易损部位检测系统 |
CN116341419A (zh) * | 2023-05-17 | 2023-06-27 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种流固耦合的数值确定方法、系统及电子设备 |
CN117172161A (zh) * | 2023-11-03 | 2023-12-05 | 北京大学 | 一种流场模拟方法、装置、计算机设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2204565A1 (en) * | 2008-12-31 | 2010-07-07 | Caterpillar, Inc. | Engine generator set control |
CN103778326A (zh) * | 2014-01-09 | 2014-05-07 | 昆明理工大学 | 一种基于预测刚体与流体耦合作用的浸入边界反馈力方法 |
CN103970989A (zh) * | 2014-04-15 | 2014-08-06 | 昆明理工大学 | 一种基于流固界面一致条件的浸入边界流场计算方法 |
-
2015
- 2015-04-30 CN CN201510217427.1A patent/CN104850689B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2204565A1 (en) * | 2008-12-31 | 2010-07-07 | Caterpillar, Inc. | Engine generator set control |
CN103778326A (zh) * | 2014-01-09 | 2014-05-07 | 昆明理工大学 | 一种基于预测刚体与流体耦合作用的浸入边界反馈力方法 |
CN103970989A (zh) * | 2014-04-15 | 2014-08-06 | 昆明理工大学 | 一种基于流固界面一致条件的浸入边界流场计算方法 |
Non-Patent Citations (1)
Title |
---|
王文全: "薄壁结构流固耦合数值模拟及计算方法研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105760602A (zh) * | 2015-12-30 | 2016-07-13 | 南京航空航天大学 | 一种有限体积加权基本无振荡格式的全流场数值模拟方法 |
CN106126797A (zh) * | 2016-06-20 | 2016-11-16 | 华中科技大学 | 一种水下探测装置无动力运动预报自动计算方法 |
CN106202782A (zh) * | 2016-07-20 | 2016-12-07 | 昆明理工大学 | 一种水轮机导叶主动旋转输入功率的数值计算方法 |
CN106202782B (zh) * | 2016-07-20 | 2019-11-08 | 昆明理工大学 | 一种水轮机导叶主动旋转输入功率的数值计算方法 |
CN106407600A (zh) * | 2016-10-18 | 2017-02-15 | 南京航空航天大学 | 一种基于流‑固耦合及代理模型的入水仿真分析方法 |
CN106407600B (zh) * | 2016-10-18 | 2019-11-15 | 南京航空航天大学 | 一种基于流-固耦合及代理模型的入水仿真分析方法 |
CN107895093A (zh) * | 2017-12-08 | 2018-04-10 | 兰州理工大学 | 一种风力机叶片流固耦合模态设计方法 |
CN108038329A (zh) * | 2017-12-25 | 2018-05-15 | 大连理工大学 | 一种微胶囊自愈合材料损伤和愈合机理的数值仿真方法 |
CN108038329B (zh) * | 2017-12-25 | 2021-02-26 | 大连理工大学 | 一种微胶囊自愈合材料损伤和愈合机理的数值仿真方法 |
CN108763741A (zh) * | 2018-05-28 | 2018-11-06 | 青岛科技大学 | 一种液压软管流固耦合数值预测方法 |
CN110717285A (zh) * | 2019-09-11 | 2020-01-21 | 哈尔滨工程大学 | 一种大幅六自由度运动的流固耦合模拟方法 |
CN110633556B (zh) * | 2019-10-24 | 2020-05-26 | 南京航空航天大学 | 一种陶瓷基复合材料流固耦合响应计算方法 |
CN110633556A (zh) * | 2019-10-24 | 2019-12-31 | 南京航空航天大学 | 一种陶瓷基复合材料流固耦合响应计算方法 |
CN111062154A (zh) * | 2019-11-18 | 2020-04-24 | 华北水利水电大学 | 一种南水北调渠系建筑物流固耦合计算方法 |
CN111950174A (zh) * | 2020-07-08 | 2020-11-17 | 上海交通大学 | 流固耦合的计算方法、装置及电子设备 |
CN111950174B (zh) * | 2020-07-08 | 2024-05-14 | 上海交通大学 | 流固耦合的计算方法、装置及电子设备 |
CN114757082A (zh) * | 2022-03-07 | 2022-07-15 | 同济大学 | 一种基于拉格朗日-欧拉稳定配点的流固耦合计算方法 |
CN114757082B (zh) * | 2022-03-07 | 2024-04-12 | 同济大学 | 一种基于拉格朗日-欧拉稳定配点的流固耦合计算方法 |
CN114840889A (zh) * | 2022-07-04 | 2022-08-02 | 海南浙江大学研究院 | 一种网箱网衣受力易损部位检测系统 |
CN116341419A (zh) * | 2023-05-17 | 2023-06-27 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种流固耦合的数值确定方法、系统及电子设备 |
CN116341419B (zh) * | 2023-05-17 | 2023-08-01 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种流固耦合的数值确定方法、系统及电子设备 |
CN117172161A (zh) * | 2023-11-03 | 2023-12-05 | 北京大学 | 一种流场模拟方法、装置、计算机设备及存储介质 |
CN117172161B (zh) * | 2023-11-03 | 2024-02-02 | 北京大学 | 一种流场模拟方法、装置、计算机设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104850689B (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104850689A (zh) | 一种基于固定网格技术的流固耦合计算方法 | |
CN110298105B (zh) | 饱和多孔介质大变形分析的ccpdi-impm方法 | |
CN108920768B (zh) | 一种针对弹性薄壁结构的流固耦合方法 | |
CN101908088A (zh) | 一种基于时域双向迭代的叶轮机叶片颤振应力预测方法 | |
CN109446581B (zh) | 一种波浪作用下浮体的水动力响应的测量方法及系统 | |
CN103778326A (zh) | 一种基于预测刚体与流体耦合作用的浸入边界反馈力方法 | |
CN103853921A (zh) | 一种大变形超弹性结构流激振动特性预测方法 | |
Suleiman | Text Book on Dynamic Relaxation Method | |
Chen et al. | Nonlinear aeroelastic analysis of an airfoil-store system with a freeplay by precise integration method | |
CN112949065B (zh) | 模拟层状岩体力学行为的双尺度方法、装置、存储介质及设备 | |
CN109977475A (zh) | 一种用于三维复杂管路流固耦合计算的动网格更新方法 | |
CN112001109A (zh) | 再生核粒子算法实现结构冲击动力学仿真方法 | |
CN106354954A (zh) | 一种基于叠层基函数的三维力学模态仿真模拟方法 | |
CN113792440B (zh) | 计算浮式结构物在非稳态载荷作用下结构力响应的方法 | |
CN114638143B (zh) | 一种适用于模拟水弹性问题的耦合数值计算方法 | |
CN106021186B (zh) | 一种高效求解大规模非线性随机结构系统状态的多尺度迭代方法 | |
CN105224770B (zh) | 高速轻载机构非线性动态系统结构拓扑优化方法 | |
Lv et al. | Multiscale method for geometrical nonlinear analysis of fluid actuated cellular structures with arbitrary polygonal microstructures | |
CN104866712A (zh) | 一种新的投影浸入边界法的速度校正方法 | |
CN104809301A (zh) | 一种反映非线性刚柔混合连接特性的时域子结构方法 | |
Muffo et al. | Interface Velocity Consistency in time-accurate flow simulations on dynamic meshes | |
CN103761402A (zh) | 一种确定结构多尺度动力分析中子模型尺寸的方法 | |
Zorkipli et al. | Simulation of aeroelastic system with aerodynamic nonlinearity | |
Chae et al. | Efficient Coupled Time Integration Methods for Transonic Aeroelastic Analysis | |
CN105335565A (zh) | 高速轻载机构非线性动态系统结构进化优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |