CN105808792B - 一种贮箱晃动质量的数值计算方法 - Google Patents
一种贮箱晃动质量的数值计算方法 Download PDFInfo
- Publication number
- CN105808792B CN105808792B CN201410835095.9A CN201410835095A CN105808792B CN 105808792 B CN105808792 B CN 105808792B CN 201410835095 A CN201410835095 A CN 201410835095A CN 105808792 B CN105808792 B CN 105808792B
- Authority
- CN
- China
- Prior art keywords
- tank
- liquid
- slosh
- mode
- excitation
- 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
Abstract
本发明属于液体动力学技术领域,具体涉及一种贮箱晃动质量的数值计算方法,目的在于解决现有技术难以获取晃动质量参数的问题。该方法包括建立贮箱结构有限元模型、建立贮箱液体有限元模型、施加重力加速度、获取贮箱液体晃动的各阶模态和振型和反应谱分析步骤。本发明利用数值计算中的优势将复杂的晃动力学平衡振动微分方程进一步简化为单质量块单自由度的力学平衡方程,避免了实际晃动试验中需要多点激励以及拟合求解的麻烦和误差,只需要通过单点激励和简单计算就能获得晃动质量参数。利用振型参与因子和反应谱方法避免了时域激励响应的稳定时间长的问题。在提出参数方面只需提取晃动的合力,在数值提取上具有更高的精度。
Description
技术领域
本发明属于液体动力学技术领域,具体涉及一种贮箱晃动质量的数值计算方法。
背景技术
贮箱晃动问题在工程界的研究始于二十世纪三十年代地震对大型水闸、油箱、水库以及水塔等大型贮水设备的影响。针对大型火箭贮箱液体晃动问题也早在二十世纪六、七十年代就开展了大量的研究工作。
基于线性假设的贮箱内液体晃动理论可以说已经非常成熟,工程上常采用弹簧-质量-阻尼或是阻尼单摆等效力学模型来描述液体晃动基本上可以满足工程应用要求。贮箱内液体晃动的阻尼单摆等效力学模型如图1所示。
其中xi(t)为晃动激励位移;是第i阶晃动质量响应位移;mi为第i阶晃动质量;m0为晃动过程中相对静止质量;li为第i阶晃动等效摆长;ci为第i阶晃动阻尼系数;hi为第i阶晃动等效摆悬挂点高度。
虽然线性晃动理论已经非常成熟,但是在贮箱获取等效力学模型的参数尤其是贮箱晃动质量时,采用的理论以及工程计算方法存在一定的局限性,而试验研究程序比较复杂、研究费用高、周期长。因此,国内许多学者一直在寻找一种简单易行的方法来补充或是取代现有的模型参数获取方法。
在所有的等效力学模型的参数的获取中,晃动质量参数的获取是最困难的。
一些学者发展了晃动参数的工程计算方法,在一定的程度上突破工程计算方法在贮箱形状的局限性,但是对于具有复杂结构的贮箱要准确地获得晃动质量参数仍然十分困难。
随着数值计算方法的发展,有限体积的VOF方法以及液体有限单元应用的成熟,许多专家学者开始使用有限元法、有限体积法进行虚拟晃动试验,但是该数值方法中使用时域强迫激励存在谐振响应稳定时间长、频率干涉以及拟合误差等问题,因此通过该方法获取晃动质量参数时仍然比较困难且精度不高。
工程计算与数值计算相结合的方法]综合工程计算方法的快捷以及有数值计算方法的应用范围广的优势,因此得到国内许多学者的青睐。但是通过从数值计算中提取特征函数再结合工程计算方法时不仅需要自编处理程序,而且还需要面对诸多繁琐计算公式,因此在工程应用中要使用该方法获得晃动质量参数过程不仅十分繁琐,而且对实施操作人员理论水平要求很高。
发明内容
本发明的目的在于解决现有技术难以获取晃动质量参数的问题,提供一种贮箱晃动质量的数值计算方法。
本发明是这样实现的:
一种贮箱晃动质量的数值计算方法,包括如下步骤:
第一步:建立贮箱结构有限元模型;
第二步:建立贮箱液体有限元模型;
第三步:施加重力加速度;
第四步:获取贮箱液体晃动的各阶模态和振型;
第五步:反应谱分析;
第六步:提取贮箱受到的激励加速度方向的合力;
第七步:获得分析晃动阶次的晃动质量参数。
如上所述的建立贮箱结构有限元模型步骤,提取贮箱的内表面,用壳单元给贮箱划分有限单元网格,降低质量密度至液体密度的0.1%,在数值上虚化贮箱结构的质量效应,仅提供给液体晃动边界;约束贮箱底端顶点平动和转动6个自由度。
如上所述的建立贮箱液体有限元模型步骤,贮箱液体使用势流体单元模拟,与结构耦合的位置使用共节点、不拟合的力传递耦合方法;赋予势流体单元体积模量以及密度属性,液体表面设置为自由表面边界。
如上所述的获取贮箱液体晃动的各阶模态和振型步骤,使用利用振型参与因子和反应谱法计算所得到的前三阶晃动模态,屏蔽液体二、三阶晃动模态,解耦各模态之间的相互影响,仅留下一阶的液体晃动模态;
得到晃动力学方程为:
式中,i——自然数序列,1、2、3…;ω——外部激励频率;——贮箱加速度;F(ω)——贮箱晃动力;ωi——第i阶晃动频率。
如上所述的三阶晃动模态均为非对称模态。
如上所述的反应谱分析步骤,使用正弦交变场加速激励对需要分析的模态进行激励,设定正弦交变场加速度角频率与分析阶次晃动模态频率之比小于0.01,得到如下公式:
如上所述的反应谱分析步骤,设定正弦交变场加速度频率为0.001HZ,幅值为单位1。
如上所述的提取贮箱受到的激励加速度方向的合力步骤,在第五步得到的结果中提取贮箱模型在激励方向的支反力,贮箱模型在激励方向的支反力幅值等于贮箱所受到液体晃动的晃动合力F(ω)的幅值。
如上所述的获得分析晃动阶次的晃动质量参数步骤,在提取贮箱模型在激励方向的支反力后,根据公式(2)计算获得分析阶次的晃动质量参数。
本发明的有益效果在于:
本发明采用建立贮箱结构有限元模型、建立贮箱液体有限元模型、施加重力加速度、获取贮箱液体晃动的各阶模态和振型、反应谱分析、提取贮箱受到的激励加速度方向的合力和获得分析晃动阶次的晃动质量参数步骤,利用数值计算中的优势将复杂的晃动力学平衡振动微分方程进一步简化为单质量块单自由度的力学平衡方程,从而避免了实际晃动试验中需要多点激励以及拟合求解的麻烦和误差,只需要通过单点激励和简单计算就能获得晃动质量参数。利用振型参与因子和反应谱方法避免了时域激励响应的稳定时间长的问题。而在提出参数方面只需要提取晃动的合力,而不是一些数值计算方法中采用的晃动波高,因此在数值提取上具有更高的精度。在整个计算过程中不需要提取贮箱结构的型函数,只需要考虑贮箱外形以及箱内装置的耦合边界影响,因此该方法可以推广应用到任意形状以及装有各种液体管理装置的贮箱液体晃动分析。
附图说明
图1是贮箱内液体晃动的等效力学模型的示意图;
图2是本发明的一种贮箱晃动质量的数值计算方法的流程图;
图3是采用本发明的方法计算得到的一阶晃动质量曲线。
具体实施方式
下面结合附图和具体实施例对本发明的一种贮箱晃动质量的数值计算方法进行描述:
一种贮箱晃动质量的数值计算方法,包括如下步骤:
第一步:建立贮箱结构有限元模型;
提取贮箱的内表面,用壳单元给贮箱划分有限单元网格。降低质量密度至液体密度的0.1%,在数值上虚化贮箱结构的质量效应,仅提供给液体晃动边界。约束贮箱底端顶点平动和转动6个自由度。本例中,对于直径为一米的圆筒平底贮箱,计算贮箱内液体晃动质量,设定壳单元厚度为10cm,材料密度设定为1kg/m3,弹性模量设定为2.1x1012。
第二步:建立贮箱液体有限元模型;
贮箱液体使用势流体单元模拟,与结构耦合的位置使用共节点、不拟合的力传递耦合方法。赋予势流体单元体积模量以及密度属性,液体表面设置为自由表面边界。本例中设定势流体单元密度为1000kg/m3,体积模量为2.56x109。
第三步:施加重力加速度;
整个计算域施加重力加速度,本例中设定重力加速度为9.8m/s2。
第四步:获取贮箱液体晃动的各阶模态和振型;
使用利用振型参与因子和反应谱法计算所得到的前三阶晃动模态,三阶晃动模态均为非对称模态。屏蔽不需要分析的液体二、三阶晃动模态,解耦各模态之间的相互影响,仅留下一阶的液体晃动模态。其中,前三阶振型如图2所示。在屏蔽不需要分析的二、三阶液体晃动模态后,晃动力学方程可以简化为
式中,i——自然数序列,为1、2、3…;ω——外部激励频率;——贮箱加速度;F(ω)——贮箱晃动力;ωi——第i阶晃动频率。
第五步:反应谱分析;
使用正弦交变场加速激励对需要分析的模态进行激励,设定正弦交变场加速度角频率与分析阶次晃动模态频率之比小于0.01。在本实施例中,设定正弦交变场加速度频率为0.001HZ,幅值为单位1,即的幅值为1,从而将公式(1)简化为:
第六步:提取贮箱受到的激励加速度方向的合力;
在第五步得到的结果中提取贮箱模型在激励方向的支反力,由于贮箱模型底端顶点被固支,因此贮箱模型在激励方向的支反力幅值就等于贮箱所受到液体晃动的晃动合力F(ω)的幅值。
第七步:获得分析晃动阶次的晃动质量参数;
在提取贮箱模型在激励方向的支反力后,根据公式(2)计算获得分析阶次的晃动质量参数。
为验证本方法的可行性和准确性,使用该方法计算了圆柱贮箱的不同液高时的晃动质量。其对比结果如图3所示,从图中可以看出,计算与理论值的误差在不同液高时均小于1%,远远高于试验要求5%的误差要求,因此证明该方法在计算晃动质量时是具有非常高的精度。
本方法的可以通过adina、ansys、nastran、abaqus等有限元软件实现。
本发明采用建立贮箱结构有限元模型、建立贮箱液体有限元模型、施加重力加速度、获取贮箱液体晃动的各阶模态和振型、反应谱分析、提取贮箱受到的激励加速度方向的合力和获得分析晃动阶次的晃动质量参数步骤,利用数值计算中的优势将复杂的晃动力学平衡振动微分方程进一步简化为单质量块单自由度的力学平衡方程,从而避免了实际晃动试验中需要多点激励以及拟合求解的麻烦和误差,只需要通过单点激励和简单计算就能获得晃动质量参数。利用振型参与因子和反应谱方法避免了时域激励响应的稳定时间长的问题。而在提出参数方面只需要提取晃动的合力,而不是一些数值计算方法中采用的晃动波高,因此在数值提取上具有更高的精度。在整个计算过程中不需要提取贮箱结构的型函数,只需要考虑贮箱外形以及箱内装置的耦合边界影响,因此该方法可以推广应用到任意形状以及装有各种液体管理装置的贮箱液体晃动分析。
Claims (6)
1.一种贮箱晃动质量的数值计算方法,包括如下步骤:
第一步:建立贮箱结构有限元模型;
第二步:建立贮箱液体有限元模型;
第三步:施加重力加速度;
第四步:获取贮箱液体晃动的各阶模态和振型;
第五步:反应谱分析;
第六步:提取贮箱受到的激励加速度方向的合力;
第七步:获得分析晃动阶次的晃动质量参数;
所述的第五步:反应谱分析,使用正弦交变场加速激励对需要分析的模态进行激励,设定正弦交变场加速度角频率与分析阶次晃动模态频率之比小于0.01,得到如下公式:
式中:F(ω)为贮箱晃动力;为贮箱加速度;mi为第i阶晃动质量;
所述的第五步:反应谱分析,设定正弦交变场加速度频率为0.001HZ,幅值为单位1;
所述的第六步:提取贮箱受到的激励加速度方向的合力,在第五步得到的结果中提取贮箱模型在激励方向的支反力,贮箱模型在激励方向的支反力幅值等于贮箱所受到液体晃动的晃动合力F(ω)的幅值。
2.根据权利要求1所述的一种贮箱晃动质量的数值计算方法,其特征在于:所述的建立贮箱结构有限元模型步骤,提取贮箱的内表面,用壳单元给贮箱划分有限单元网格,降低质量密度至液体密度的0.1%,在数值上虚化贮箱结构的质量效应,仅提供给液体晃动边界;约束贮箱底端顶点平动和转动6个自由度。
3.根据权利要求1所述的一种贮箱晃动质量的数值计算方法,其特征在于:所述的建立贮箱液体有限元模型步骤,贮箱液体使用势流体单元模拟,与结构耦合的位置使用共节点、不拟合的力传递耦合方法;赋予势流体单元体积模量以及密度属性,液体表面设置为自由表面边界。
4.根据权利要求1所述的一种贮箱晃动质量的数值计算方法,其特征在于:所述的获取贮箱液体晃动的各阶模态和振型步骤,使用利用振型参与因子和反应谱法计算所得到的前三阶晃动模态,屏蔽液体二、三阶晃动模态,解耦各模态之间的相互影响,仅留下一阶的液体晃动模态;得到晃动力学方程为:
式中,i——自然数序列,1、2、3…;ω——外部激励频率;——贮箱加速度;F(ω)——贮箱晃动力;ωi——第i阶晃动频率。
5.根据权利要求4所述的一种贮箱晃动质量的数值计算方法,其特征在于:所述的三阶晃动模态均为非对称模态。
6.根据权利要求1所述的一种贮箱晃动质量的数值计算方法,其特征在于:所述的获得分析晃动阶次的晃动质量参数步骤,在提取贮箱模型在激励方向的支反力后,根据公式(2)计算获得分析阶次的晃动质量参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410835095.9A CN105808792B (zh) | 2014-12-29 | 2014-12-29 | 一种贮箱晃动质量的数值计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410835095.9A CN105808792B (zh) | 2014-12-29 | 2014-12-29 | 一种贮箱晃动质量的数值计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105808792A CN105808792A (zh) | 2016-07-27 |
CN105808792B true CN105808792B (zh) | 2019-04-19 |
Family
ID=56979854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410835095.9A Active CN105808792B (zh) | 2014-12-29 | 2014-12-29 | 一种贮箱晃动质量的数值计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105808792B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106383964B (zh) * | 2016-10-10 | 2019-07-12 | 北京宇航系统工程研究所 | 一种悬挂式充液贮箱动力学建模方法 |
CN110222355A (zh) * | 2019-04-11 | 2019-09-10 | 北京宇航系统工程研究所 | 一种考虑热力耦合的低温推进剂贮箱防晃板分析方法 |
CN111397825B (zh) * | 2020-03-03 | 2022-12-16 | 东方电气(广州)重型机器有限公司 | 容器内的液体晃动效应确定方法、系统及存储介质 |
CN111487867B (zh) * | 2020-04-16 | 2021-09-28 | 浙江大学 | 一种液体晃动模型参考自适应抑制运动控制方法 |
CN113252131B (zh) * | 2021-04-06 | 2022-10-18 | 上海宇航系统工程研究所 | 一种运载火箭上升段贮箱内推进剂余量的测量系统与方法 |
CN113553661B (zh) * | 2021-06-30 | 2022-09-16 | 航天科工火箭技术有限公司 | 获取火箭推进剂晃动参数的方法及计算机可读存储介质 |
CN114428119B (zh) * | 2022-01-20 | 2023-11-21 | 重庆大学 | 一种用于反演各向异性特征的复合材料弹性常数的方法 |
CN115358121B (zh) * | 2022-08-12 | 2024-02-09 | 上海核工程研究设计院股份有限公司 | 一种斜底直立环形水箱流体晃动等效力学模型的建模方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102991729A (zh) * | 2012-11-30 | 2013-03-27 | 北京控制工程研究所 | 一种轻质网式表面张力贮箱 |
CN103455644A (zh) * | 2012-05-31 | 2013-12-18 | 北京宇航系统工程研究所 | 一种基于耦合质量的液体推进剂动力学模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006163625A (ja) * | 2004-12-03 | 2006-06-22 | Osaka Univ | 有限要素解析方法、有限要素解析装置、及びコンピュータプログラム |
JP2013008223A (ja) * | 2011-06-24 | 2013-01-10 | Toyota Motor Corp | スポット溶接部応力算出方法 |
-
2014
- 2014-12-29 CN CN201410835095.9A patent/CN105808792B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103455644A (zh) * | 2012-05-31 | 2013-12-18 | 北京宇航系统工程研究所 | 一种基于耦合质量的液体推进剂动力学模拟方法 |
CN102991729A (zh) * | 2012-11-30 | 2013-03-27 | 北京控制工程研究所 | 一种轻质网式表面张力贮箱 |
Non-Patent Citations (1)
Title |
---|
模态分析在液体晃动等效力学模型建模中的应用;郭经经等;《上海航天》;20101231(第6期);第11-15页 |
Also Published As
Publication number | Publication date |
---|---|
CN105808792A (zh) | 2016-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105808792B (zh) | 一种贮箱晃动质量的数值计算方法 | |
Xue et al. | Numerical study of ring baffle effects on reducing violent liquid sloshing | |
Wang et al. | Finite element analysis of two-dimensional nonlinear sloshing problems in random excitations | |
Windt et al. | Evaluation of the overset grid method for control studies of wave energy converters in OpenFOAM numerical wave tanks | |
Zhang et al. | Apply MPS method to simulate liquid sloshing in LNG tank | |
CN106934185A (zh) | 一种弹性介质的流固耦合多尺度流动模拟方法 | |
Jin et al. | Viscous effects on liquid sloshing under external excitations | |
CN104317985A (zh) | 一种基于界带有限元和拉格朗日坐标的流体仿真方法 | |
Li et al. | Numerical studies on sloshing in rectangular tanks using a tree-based adaptive solver and experimental validation | |
Matioc | An explicit solution for deep water waves with Coriolis effects | |
Mandal et al. | Pressure based Eulerian approach for investigation of sloshing in rectangular water tank | |
Malik et al. | Transient numerical simulations for hydrodynamic derivatives predictions of an axisymmetric submersible vehicle | |
Peñas-López et al. | The history effect on bubble growth and dissolution. Part 2. Experiments and simulations of a spherical bubble attached to a horizontal flat plate | |
CN114510887B (zh) | 一种基于代理模型的海洋条件下大容器液位快速预测方法 | |
Huang et al. | Time-domain simulation of tank sloshing pressure and experimental validation | |
Li et al. | An approximate analytical solution of sloshing frequencies for a liquid in various shape aqueducts | |
Larsen | Impact loads on circular cylinders | |
Yao et al. | An improved SPH-FEM coupling approach for modeling fluid–structure interaction problems | |
Yeylaghi et al. | Sph modeling of hydrodynamic loads on a point absorber wave energy converter hull | |
Zhang et al. | MPS-FEM coupled method for interaction between sloshing flow and elastic structure in rolling tanks | |
Davis et al. | Systematic identification of drag coefficients for a heaving wave follower | |
Bašić et al. | Analytical and numerical computation of added mass in ship vibration analysis | |
Jendzelovsky et al. | Numerical modeling of cylindrical tank and compare with experiment | |
Ozdemir et al. | FSI methods for seismic analysis of sloshing tank problems | |
Eswaran et al. | Numerical simulation of 2D and 3D sloshing waves in a regularly and randomly excited container |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |