Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

储液容器内液体晃荡的非线性动力学分析

李文盛 赵友清 贾善坡 王凯 谭继可

王江波, 丁俊升, 王晓东, 杜忠华, 高光发. 粗骨料粒径对混凝土动态压缩行为的影响研究[J]. 爆炸与冲击, 2022, 42(2): 023101. doi: 10.11883/bzycj-2021-0147
引用本文: 李文盛, 赵友清, 贾善坡, 王凯, 谭继可. 储液容器内液体晃荡的非线性动力学分析[J]. 爆炸与冲击, 2014, 34(1): 87-92.
WANG Jiangbo, DING Junsheng, WANG Xiaodong, DU Zhonghua, GAO Guangfa. Effect of coarse aggregate size on the dynamic compression behavior of concrete[J]. Explosion And Shock Waves, 2022, 42(2): 023101. doi: 10.11883/bzycj-2021-0147
Citation: Li Wen-sheng, Zhao You-qing, Jia Shan-po, Wang Kai, Tan Ji-ke. Numerical analysis on liquid sloshing in storage container by nonlinear dynamics method[J]. Explosion And Shock Waves, 2014, 34(1): 87-92.

储液容器内液体晃荡的非线性动力学分析

基金项目: 国家自然科学基金项目(50978033);湖北省教育厅科学研究计划项目(Q20101301)
详细信息
    作者简介:

    李文盛(1968—), 男, 副教授, 硕士生导师

    通讯作者:

    Zhao You-qing, z_uqing@163.com

  • 中图分类号: O351.3

Numerical analysis on liquid sloshing in storage container by nonlinear dynamics method

Funds: Supported by the National Natural Science Foundation of China (50978033)
  • 摘要: 基于非线性波动理论模型,求解储液容器内液体晃动的固有频率、模态及动力学响应问题。流体使用us-up状态方程,利用ABAQUS软件的自适应网格技术,建立储液容器液体晃动数学模型,通过施加水平简谐激励得到液体晃动的固有频率和模态,并与解析解对比,验证了该方法的准确性与可行性。然后,分析了矩形储液容器在多种激励作用下液体非线性晃动响应特性。
  • 混凝土材料是一种由水泥砂浆基体与粗骨料颗粒组成的应变率敏感的工程复合材料,其动态冲击下的力学行为与准静态状况下有着显著的区别。随着应变率的增加,混凝土材料的抗拉强度与抗压强度都有明显的增加,这种率效应是由于水泥基体的微裂纹扩展和粘弹性特性引起的[1]。SHPB装置是研究混凝土材料动态力学性能最常用的试验装置[2],可以获得其应变率在101~103 s−1范围内的冲击力学行为。粗骨料通常占混凝土混合体积40%以上,在混凝土的力学性能和物理特性方面起着重要的作用[3]。事实上,与传统的金属材料不同,在材料动态力学性能测试试件的尺度上,粗骨料的影响是不可忽视的[4-6]

    当前,土木工程领域的国内外学者分别对大尺寸混凝土试件静态力学性能和准细观尺寸效应开展了一些试验研究,研究表明:不同受力状态、不同截面形状、不同强度的混凝土准静态力学性能的尺寸效应非常明显,且该效应与试件尺寸、试件形状和材料强度耦合[7-8]。同时,混凝土准静态三点弯曲试验、准静态压缩和拉伸试验表明,粗骨料尺寸和含量对混凝土准静态力学性能、破坏能量和裂纹扩展形态皆具有非常明显的影响[9-11]。Hao等[12]对最大骨料粒径为4和8 mm的混凝土和砂浆试件进行了SHPB试验,试验结果表明:含骨料的混凝土试件动态增长因子总高于砂浆试件,并且随着骨料粒径的增大,混凝土的动态增长因子减小。Kim等[13]对最大骨料粒径为13、19和25 mm的砂浆和混凝土进行了SHPB试验发现:纯率动态增长因子随着应变率的增加而增加,但不同最大粗骨料尺寸的总体结果对动态增长因子没有显示出很强的影响。Li等[14]对骨料最大粒径为7.5 mm的方形和圆柱形试件进行了SHPB试验,结果表明圆柱体试件的动态抗压强度高于立方体试件,即试件形状对抗压强度也有一定的影响。Grote等[15]进行了骨料对平板冲击试验混凝土材料动态增长因子的影响,并对粗骨料体积率从0%增加到42%进行数值模拟研究,结果表明混凝土试件的动态强度得到了提高。Cadoni等[4]对含有粗骨料最大尺寸为25 mm的混凝土试样进行Hopkinson钢筋束拉伸试验,发现骨料的含量越小,试件的动态拉伸强度和能量吸收能力越高。

    动态增长系数(动态增长因子,动态抗压强度与准静态抗压强度之比)常用来衡量混凝土材料的应变率效应。研究人员提出了各种经验公式来描述混凝土材料的动态增长因子与应变率的关系[1, 15-20],其中使用最广泛的是CEB-FIP方程[19]。CEB模型表明,在应变率为30 s−1时,动态增长因子和对数应变率(lg˙ε)曲线中存在斜率的突变。Bischoff等[20]根据不同文献的研究结果对动态增长因子与应变率的关系进行了综合评述,认为在动态冲击下(约10 s−1),混凝土的抗压强度比准静态压缩强度高出85%~100%,但试验结果有很大的离散性。高光发等[21-23]对近30年相关文献中混凝土的动态压缩和拉伸试验数据进行总结,得到了混凝土材料在不同应变率区间的强度动态增长因子预测表达式。

    综上所述,发现混凝土动态力学性能的研究主要集中在应变率效应方面,而粗骨料粒径对混凝土动态冲击性能的研究不够充分。基于此,本文中对粗骨料平均粒径对混凝土抗压强度及动态增长因子的影响进行研究,对3种不同平均骨料粒径的混凝土及砂浆进行SHPB试验,并定义无量纲骨料尺寸¯D¯D=da/ds,粗骨料平均粒径da与试件直径ds之比)用来比较粗骨料平均粒径对混凝土动态抗压强度动态增长因子的影响。

    试验材料为中强度砂浆和混凝土,试验原材料包括硅酸盐水泥、河砂、粗骨料、矿粉、水和减水剂,不同粗骨料粒径的混凝土与砂浆相比,除了增加相同体积分数的粗骨料以外,其它材料及配比完全一样。另外,为了保证强度达到设计要求以及试件的准确性,对粗骨料进行了水洗和烘干的预处理,减小其含泥量。由于本试验中研究不同粒径对混凝土力学性能的影响,故对粗骨料进行筛选,得到4~8 mm、10~14 mm和22~26 mm这3个范围的粗骨料尺寸,可认为粗骨料等效平均粒径分别为6、12和24 mm。

    试验中加工了砂浆强度型号为M60的砂浆试件,在此基础上制备出3种不同粗骨料尺寸的C60混凝土,配合比如表1所示,需要说明的是,这里的C60标号与传统的标号意义不一致,这里是指M60砂浆添加粗骨料后的混凝土,其目的是为了研究粗骨料对混凝土力学性能的影响。同时选用3种粗骨料,其粗骨料平均粒径分别为6、12和24 mm,每种粒径的质量占比都相同。混凝土试件的命名规则为砂浆等级以及加入粗骨料的平均粒径,比如对于C60-G6混凝土为M60的砂浆和平均粒径为6 mm的粗骨料均匀混合加工而成。

    表  1  不同材料试样的配比
    Table  1.  Mix proportion of different grades of mortar and concrete
    试件w(水泥)/(kg·m−3w(砂子)/(kg·m−3w(粗骨料)/(kg·m−3w(矿粉)(kg·m−3w(水)/(kg·m−3w(减水剂)/(kg·m−3
    M607141000222232.22
    C60547 767897171711.70
    下载: 导出CSV 
    | 显示表格

    按照表1的配合比制备尺寸为1 000 mm×1 000 mm×400 mm混凝土和砂浆靶板,根据规定进行常规养护时间28 d,然后通过钻孔取芯切割的方法获取试件,对试件的两端面进行打磨,使试件的两端面平行度小于0.2°,打磨后的试件直径为68 mm,准静态和动态压缩试件长径比分别为2.0和0.5。

    利用MTS试验机进行了砂浆和不同粒径混凝土试件的准静态压缩试验,应变率为10−4 s−1,试件的侧面中间对称贴有一组应变片,用于记录试验中试件的应变变化,试件放置于压头和底座之间,如图1所示。试验中,加载力p通过试验机的数据采集软件获得,试件的应力σ=p/As可通过公式σ=p/As计算,式中As为试件的横截面积。

    图  1  准静态压缩试验
    Figure  1.  Quasi-static compression test

    图2为准静态压缩条件下不同类型试件的应力-应变曲线。需要说明的是,每种试件共获得3条有效的曲线,准静态压缩的屈服强度应为其3个有效数据的平均值,为了展示材料的应力应变曲线,选取应力峰值在中间的曲线为典型的应力-应变曲线。

    图  2  准静态压缩下的应力-应变曲线
    Figure  2.  Quasi-static compressive stress-strain curves
    1.3.1   SHPB试验装置

    试验采用直径为80 mm的SHPB装置,如图3所示。该装置由发射器、撞击杆、入射杆、透射杆、能量吸收装置和数据采集系统组成,杆的几何参数和材料参数见表2。根据一维弹性波理论,试件与杆件的接触面上会产生透射和反射现象。入射脉冲和反射脉冲由粘贴在入射杆中间位置的一对轴对称的应变片记录,而透射脉冲由在透射杆上另一对轴对称粘贴的应变片测量。根据入射/透射杆上应变片测量到的入射、反射和透射应变信号,可通过下列公式计算出试样的应变率˙ε、应变ε以及应力σ

    图  3  SHPB试验装置
    Figure  3.  Schematic of the SHPB test system
    表  2  SHPB装置中部件的主要参数
    Table  2.  Specifications of the SHPB experimental system
    部件材质几何参数物理参数
    直径/mm长度/mmEb/GPaρb/(g·cm−3)νbcb/(m·s−1)
    撞击杆40Cr钢80500/10002107.850.225210
    入射杆40Cr钢8060002107.850.225210
    透射杆40Cr钢8040002107.850.225210
    下载: 导出CSV 
    | 显示表格
    {˙ε(t)=2cbLsεr(t)ε(t)=2cbLst0εr(t)dtσ(t)=AbAsEbεt(t) (1)

    式中:Ab为杆的横截面积;AsLs分别为试件的横截面积和长度;cbEb为杆的弹性波速和杨氏模量;εr(t)εt(t)为反射波和透射波。

    1.3.2   双脉冲整形器使用

    混凝土的动态压缩破坏应变通常小于1%,因此需要缓慢升高的入射应力使试件内部应力均匀。然而,常规SHPB试验(无脉冲整形器)产生的入射应力迅速增大,会导致混凝土试件在达到峰值应力之前发生破坏。为此,在SHPB试验中采用了脉冲整形技术,以保证试件处于动态应力平衡和恒应变率加载状态。在本研究中,将直径为8~12 mm和24~40 mm的橡胶和T2紫铜放置在入射杆的撞击端作为脉冲整形器,如图4(a)所示。试件放置于入射杆与透射杆之间(图4(b)),撞击杆冲击时,较小的橡胶首先变形,消除了入射脉冲中的高频振荡波,保证了试件端面与加载杆的充分接触。然后,较大的紫铜片对入射脉冲进行整形,产生加载脉冲,使试样失效。通过改变脉冲整形器的厚度和直径,可以调整入射脉冲的形状。

    图  4  试验时整形器与试件位置图
    Figure  4.  Position of the pulse shaper and the specimen in the experiment
    1.3.3   应变率的确定

    混凝土试件的SHPB试验中整个加载区域内很难实现恒定的应变率,而关于应变率的取法也众说纷纭。有的文献对整个加载过程的应变率求平均值认为是代表应变率[14],然而,这种不能表征试件破坏时的真实应变率,可能导致得到的应变率大于实际情况。Kim等[13]将试件应力峰值达到最大点对应的应变率定义为代表应变率。但是这种方法只能代表试件失效时的应变率,并不能完全表征试件在整个压缩过程中的应变率。本文中,代表应变率采用应变率平台段初始至试件达到最大应力之间范围的平均值,并对此进行线性拟合,如图5所示。这种应变率的选取方法可以表征试验过程中试件应力上升至破坏的整个过程[24]

    图  5  典型试件的试验结果
    Figure  5.  Typical test results of specimens
    1.3.4   动态应力平衡

    众所周知,利用SHPB装置进行混凝土材料的动态力学性能测试时,必须保证试件在加载过程中的恒应变率变形和轴向的动态应力平衡。这是由于在SHPB动态压缩试验时,试件的应力、应变和应变率计算都是基于一维应力波和轴向动态应力平衡假设。试件两端的应力可计算获得弹性关系:

    σ1=AbAs(σi+σr) (2)
    σ2=AbAsσt (3)

    式中:σ1σ2分别为试件前和后的应力;σiσrσt分别为入射、反射和透射应力。则平均应力σave可表示为:

    σave=12(σ1+σ2) (4)

    所以,为了表征试件在破坏时是否达到动态应力平衡状态,一些学者定义了一个动态平衡因子R[25, 26]用来表示试件前后应力的差异,定义为:

    R=|f1f2fave| (5)

    式中:f1f2fave分别为试件的前后及平均压缩强度,即为σ1σ2的最大值。一般认为,当试件两端应力差与平均应力之比小于5%时可视为动态应力平衡[13]。同时,混凝土类材料SHPB试验相对较准确的应变率上限应满足:

    ˙εεscsnLs (6)

    式中:n为应力波在试件轴向的反射次数,对于脆性材料,n=4。本试验所有试件的长度为34 mm,密度为2.35 g/cm3。根据图2可知,杨氏模量最小25.5 GPa,则试件的波速根据cs=Es/Esρsρs可计算为3294 m/s,破坏应变假设为0.5%,则可以计算出最大应变率应该小于121.1/s。图6为试验应变率为55和150 s−1时的试验结果,在试件应力峰值初,应变率为55 s−1R为3.2%,试验过程中试件两端的应力大致相等,而对于应变率为150 s−1R为15.7%,超过动态应力平衡规定的5%,并且试件两端的应力-时间曲线重复性差,故可以认为当试验的应变率大于上限时,试验过程中难以保证试件内部的动态应力平衡。

    图  6  试件的动态应力平衡状态
    Figure  6.  Dynamic stress equilibrium in the specimens

    图7为M60、C60-G6、C60-G12和C60-G24材料在不同应变率下典型的应力-应变曲线,并与准静态压缩试验下相应的应力-应变曲线进行了比较。从图中可以看出,4种材料的应力应变行为具有明显的应变率效应,即随着应变率的增加,抗压强度逐渐增大;同时,4种材料的应力-应变曲线呈现出相似的变化趋势,即初期应力与应变呈线性关系,随后为非线性阶段,当应力到达峰值以后急剧减小,使得试件破坏。应力应变曲线的下降阶段说明了峰值应力后的应变软化现象,试件开始出现不同程度的破碎。砂浆和混凝土试件的应力-应变曲线在动态和准静态压缩的初始阶段重合,这说明4种材料的弹性模量没有应变率效应。

    图  7  不同应变率下试件的动态应力-应变曲线
    Figure  7.  Dynamic stress-strain curves of specimens at various strain rates

    关于混凝土材料强度的应变率效应有不同的解释。Bischoff等[20]和Hao等[27]认为混凝土中的自由水在动态压缩时会产生Stefan效应,即内部孔洞中自由水的运动会阻碍裂纹的扩展。Cotsovos等[28]认为这是SHPB试验中试件端面摩擦效应引起的横向惯性效应,并不是真正的应变率效应。Ma等[29]认为在动态压缩试验时,试件内部裂纹的扩展速度远低于应力波的加载速度,裂缝并不只是沿薄弱区域产生和扩展,同时还会穿过粗骨料颗粒。短时间内试件内部的冲击能量累积增大,导致其强度增大。

    不同平均骨料粒径的混凝土和砂浆的动态抗压强度与应变率的关系如图8(b)所示,很明显,抗压强度随着应变率的增加而增加。图8(b)还显示,在相同的应变率范围内,C60-G12混凝土的动态抗压强度最高,与准静态压缩状况有很大的差别(图8(a))。这说明当加入了粗骨料后,混凝土的准静态压缩强度会随着骨料平均粒径的不同而出现差异,但总体都低于砂浆的抗压强度。然而,加入粗骨料后会增强混凝土的动态抗压强度。

    图  8  不同试件压缩强度与应变率的关系图
    Figure  8.  Relationship between the compressive strength and the strain rate of different specimens

    动态增长因子通常定义为动态压缩强度σd与准静态强度σqs之比:

    χ=σdσqs (7)

    现在已经广泛用于评价应变率效应对混凝土类材料的影响。研究表明:类混凝土材料的动态增长因子与应变率的对数呈线性关系。基于这些结果,学者们[16-171930-32]给出了混凝土材料抗压强度动态增长因子一些表达式,下面列出了几种被广泛应用的模型。

    CEB方程[19]根据应变率、准静态应变率和准静态压缩强度给出了计算动态增长因子的经验公式:

    χ=σdσqs={(˙ε/˙ε0s)kα˙ε30s1β(˙ε/˙ε0s)1/3˙ε30s1 (8)

    式中:˙ε0s=3×105s1k为常数,σd是与试件准静态压缩强度σd和参考压缩强度σrs(σrs=10MPa)有关的常量,其中,α=1/(5+9σd/σrs)lgβ=6.16α2k=1.026。

    Tedesco等[17]根据动态压缩试验给出了混凝土动态抗压强度的动态增长因子表达式:

    χ={0.000695lg˙ε+1.0581.0˙ε63.1s10.758lg˙ε0.2892.5˙ε63.1s1 (9)

    式中:从低应变率敏感到高应变率敏感的过渡应变率为63.1 s−1,比CEB模型中的要高出一些。

    Zhou等[16]基于岩石材料的SHPB试验结果拟合了动态增长因子与应变率对数的一个经验公式:

    χ={0.0225lg˙ε+1.12˙ε10s10.2713(lg˙ε)2+0.3563(lg˙ε)+1.227510s1˙ε2000s1 (10)

    Hartmann等[30]基于文献中混凝土材料的动态压缩试验结果得到了动态增长因子经验公式:

    χ=0.5(˙ε/˙ε0)0.13+0.90˙ε0=1.0s1 (11)

    Lee等[31]提出了一种混凝土抗压强度动态增长因子与纯应变速率的经验公式,此公式中消除了试验中横向惯性效应对动态增长因子的影响:

    χ=(˙ε/105)0.0147 (12)

    Al-Salloum等[32]基于最小二乘法得到了抗压强度动态增长因子与应变率的多项式分布:

    χ=(3.54˙ε+430.6)/(˙ε+447.3) (13)

    使用文献中模型的经验公式做出的曲线并和本文中试验所得的动态增长因子一同绘制在图9中。图9证实上面的公式都不能很好地描述M60砂浆和不同粗骨料平均粒径混凝土的动态增长因子演变,这些模型一部分低估或高估了试件的动态增长因子,一部分对于临界转换应变率的定义不适用本文中的试验结果。究其原因,认为现阶段文献中动态增长因子的经验公式大多数都是基于混凝土类材料的强度来进行拟合的,并没有区分砂浆、骨料粒径以及含量等对动态增长因子的细致影响。因此,当进行研究时,会发现利用上面提到的模型会与本文中的结果有一定的差异。

    图  9  试验动态增长因子与不同模型对比
    Figure  9.  Comparison of DIFs with different models

    Bischoff等[20]建议混凝土材料的动态增长因子和应变率可以表达为:

    χ=Alg(˙ε/˙εs)+1 (14)

    式中:˙εs=30s1为参考应变率,参数A可认为是不同粒径粗骨料混凝土的应变率强化因子。为此,对本M60、C60-G6、C60-G12和C60-G24材料的动态增长因子利用式(14)进行了线性拟合,相应的曲线如图10所示,具体拟合结果如表3所示。同时,对试验数据线性拟合的拟合度R2进行了评价,其值接近1,表明线性拟合结果能够对试验数据进行高精度的预测。从拟合的结果可以看出,在应变率大于30 s-1时,4种材料的试验动态增长因子由式(14)中对数曲线可以很好的拟合,可以用来预测不同骨料粒径混凝土在高应变率下的力学行为。需要说明的是,试验结果中混凝土材料的动态增长因子可能受动静态压缩试件不同长径比的影响。

    图  10  不同骨料粒径的混凝土动态增长因子与应变率关系
    Figure  10.  Fitting curves of the dynamic increase factors and the strain rate for different aggregate sizes
    表  3  不同试样的动态增长因子在公式中的拟合结果
    Table  3.  Fitting results of dynamic increace factors of different samples in the formula
    试样A˙εs/s-1R2
    M600.94300.998
    C60-G61.45300.996
    C60-G122.13300.995
    C60-G241.08300.996
    下载: 导出CSV 
    | 显示表格

    在2.2节中,基于式(14)对不同骨料的混凝土和砂浆的动态增长因子进行了拟合,为了更直观表示粗骨料平均粒径对混凝土动态增长因子的影响,本文中定义了粗骨料无量纲尺寸¯D=da/dadsds(粗骨料平均粒径da与试件直径ds之比)用来评价骨料粒径对混凝土动态抗压强度的影响。图11为粗骨料无量纲尺寸¯D与混凝土应变率强化因子A的关系图(¯D=0表示砂浆材料)。从图中可以看出,3种粗骨料无量纲尺寸的混凝土材料应变率强化因子分别为1.45、2.13和1.08,均大于砂浆材料的0.94,这说明粗骨料对混凝土类材料的动态增长因子有强化的作用。在混凝土试件测试这个尺度上,粗骨料尺寸的影响不可忽视。随着粗骨料无量纲尺寸¯D的增大,应变率强化因子A呈现先增大后减小的趋势。这是由于在动态压缩条件下,混凝土材料中裂纹扩展时会穿过粗骨料,而粗骨料的强度较大,破坏时耗散的能量就越多,导致混凝土的应变率强化因子要高于砂浆。随着粗骨料粒径的逐渐增大,裂缝扩展时的路径变长,导致其消耗的能量就越多,从而表现出其动态增长因子增大[33];当粗骨料粒径进一步增大时,使得混凝土材料在准细观尺度上的不均匀性,从而导致宏观上应变率强化因子的降低。

    图  11  无量纲粗骨料尺寸与应变率强化因子的关系
    Figure  11.  Relationship between the dimensionless coarse aggregate size and the strain rate strengthening factor

    利用80 mm口径的SHPB装置对不同粗骨料平均粒径(6、12和24 mm)的混凝土和M60砂浆进行了动态压缩试验,并且使用了双脉冲整形器来保证试验中的恒应变率加载和动态应力平衡,研究了应变率效应对混凝土和砂浆材料动态力学性能的影响,分析了粗骨料平均粒径对混凝土动态增长因子公式中应变率强化系数的影响规律,主要得到以下结论。

    (1)在准静态压缩条件下,M60砂浆的抗压强度高于C60混凝土的抗压强度,且随着粗骨料平均粒径的增加,混凝土的抗压强度呈现出先增大后恒定的趋势;

    (2)动态加载下,不同粗骨料粒径混凝土和对应砂浆材料的抗压强度都有明显的应变率效应,在本文的动态应变率范围内(40~110 s−1),C60-G12混凝土的动态抗压强度高于其他类型的混凝土和砂浆;

    (3)不同粗骨料粒径混凝土的动态增长因子均高于砂浆材料,表明粗骨料对混凝土材料的动态抗压强度具有增强作用,定义了粗骨料无量纲尺寸¯D来评价粗骨料粒径的尺寸效应,即随着¯D的增加,应变率强化系数A呈现出先增大后减小的趋势。

  • 图  1  非线性波动理论模型

    Figure  1.  The 2D model considering nonlinear wave theory

    图  2  储液系统晃动波高

    Figure  2.  Surface wave displacement for the tank-liquid system

    图  3  不同谐频下自由液面点B的波面响应

    Figure  3.  Free surface elevation of poit B for harmonic frequencies

    图  4  第1阶频率作用下液体晃动波高曲线

    Figure  4.  Variation in time of the surface wave in the first sloshing mode

    图  5  第1阶液体晃动模态图和液体晃动位移矢量图

    Figure  5.  Liquid shapes corresponding to the first sloshing mode and displacement vector plot

    图  6  不同幅值激励下液体晃动的波高曲线

    Figure  6.  Free surface elevation of liquid for different amplitudes

    图  7  EI Centro地震波和液面点AB液体晃动的波高曲线

    Figure  7.  EI Centro seismic wave and free surface displacement curve at points A and B

    表  1  液体自由晃动频率

    Table  1.   Frequencies of liquid sloshing

    n ωnum/(rad·s-1) ωana/(rad·s-1) εω/%
    1 3.20 3.20 0
    2 4.51 4.52 0.221 2
    3 5.60 5.53 1.265 8
    下载: 导出CSV
  • [1] Hasheminejad S M, Aghabeigi M. Sloshing characteristics in half-full horizontal elliptical tanks with vertical baffles[J]. Applied Mathematical Modelling, 2012, 36(1): 57-71. doi: 10.1016/j.apm.2011.02.026
    [2] Xue Mi-an, Lin Peng-zhi. Numerical study of ring baffle effects on reducing violent liquid sloshing[J]. Computers & Fluids, 2011, 52: 116-129. http://www.sciencedirect.com/science/article/pii/S0045793011002921
    [3] Amsden A A, Harlow F H. A simplified MAC technique for incompressible fluid flow calculations[J]. Journal of Computational Physics, 1970, 6(2): 322-325. doi: 10.1016/0021-9991(70)90029-X
    [4] Armenio V. An improved MAC method(SIMAC)for unsteady high-reynolds free surface flows[J]. International Journal for Numerical Methods in Fluids, 1997, 24(2): 185-214. doi: 10.1002/(SICI)1097-0363(19970130)24:2<185::AID-FLD487>3.0.CO;2-Q
    [5] 包光伟, 王振强, 张天翔.火箭推进剂液体晃动关机响应的数值仿真[J].宇航学报, 2002, 23(2): 84-88. doi: 10.3321/j.issn:1000-1328.2002.02.018

    Bao Guang-wei, Wang Zhen-qiang, Zhang Tian-xiang. Numerical simulation of slosh of the propellant fuel during the period the shut-down of the rocket[J]. Jourual of Astronautics, 2002, 23(2): 84-88. doi: 10.3321/j.issn:1000-1328.2002.02.018
    [6] Jung J H, Yoon H S, Lee C Y, et al. Effect of the vertical baffle height on the liquid sloshing in a three-dimensional rectangular tank[J]. Ocean Engineering, 2012, 44: 79-89. doi: 10.1016/j.oceaneng.2012.01.034
    [7] 刘永涛, 马宁, 顾解忡.各种激励作用下液舱内液体晃荡的计算与分析[J].船海工程, 2009, 38(5): 7-12. doi: 10.3963/j.issn.1671-7953.2009.05.002

    Liu Yong-tao, Ma Ning, Gu Xie-chong. Calculation and analysis of liquid sloshing loads in tanks under different kinds of stimulations[J]. Ship and Ocean Engineering, 2009, 38(5): 7-12. doi: 10.3963/j.issn.1671-7953.2009.05.002
    [8] 周宏, 李俊峰, 王天舒.用ALE有限元模拟的网格更新方法[J].力学学报, 2008, 40(2): 267-272. http://www.cqvip.com/qk/91029x/200802/26789636.html

    Zhou Hong, Li Jun-feng, Wang Tian-shu. Mesh update algorithm in ale finite method within free surface flow[J]. Chinese Journal of Theoretical and Applied Mechanics, 2008, 40(2): 267-272. http://www.cqvip.com/qk/91029x/200802/26789636.html
    [9] Sygulski R. Boundary element analysis of liquid sloshing in baffled tanks[J]. Engineering Analysis with Boundary Elements, 2011, 35(8): 978-983. doi: 10.1016/j.enganabound.2011.03.001
    [10] Wang C Z, Khoo B C. Finite element analysis of two-dimensional nonlinear sloshing problems in random excitations[J]. Ocean Engineering, 2005, 32(2): 107-133. doi: 10.1016/j.oceaneng.2004.08.001
    [11] ABAQUS/explicit user's manual: Version 6.4[M]. Rhode Island: Hibbit, Karlsson and Sorensen Inc, 2002.
    [12] ABAQUS/theory user's manual: Version 6.4[M]. Rhode Island: Hibbit, Karlsson and Sorensen Inc, 2002.
    [13] 梅强中.水波动力学[M].北京: 科学出版社, 1984.
  • 期刊类型引用(7)

    1. 王高辉,彭涛,卢文波,任心语. 卵石和碎石混凝土的冲击动态力学性能分析. 武汉大学学报(工学版). 2025(01): 11-19 . 百度学术
    2. 张湘茹,程月华,吴昊. 基于3D细观模型的混凝土动态压缩行为分析. 爆炸与冲击. 2024(02): 55-71 . 本站查看
    3. 卢俊平,王亚,刘磊,张成良,徐泽辉. 粉煤灰高性能混凝土动态压缩性能试验研究. 混凝土. 2024(08): 55-62 . 百度学术
    4. 周玄,王伯通,武一丁,陆文成,马铭辉,余毅磊,高光发. 霍普金森杆实验方法中材料弹性阶段杨氏模量及其曲线准确性分析. 爆炸与冲击. 2024(09): 130-143 . 本站查看
    5. 汪凯,燕远岭,赵哲,张杯,李致坤. 界面过渡区与骨料特征对混凝土强度及变形影响的数值模拟研究. 硅酸盐通报. 2023(04): 1298-1308+1322 . 百度学术
    6. 李贤,宁麟,陈亮亮,李炀,邓小江. 基于核磁共振技术分析骨料粒径对水泥砂浆微观结构的影响. 硅酸盐通报. 2023(12): 4216-4223 . 百度学术
    7. 亓晓鹏,张杰,赵婷婷,王志勇,王志华. 考虑骨料级配的混凝土靶板接触爆炸破坏模式. 兵工学报. 2023(12): 3641-3653 . 百度学术

    其他类型引用(9)

  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  721
  • HTML全文浏览量:  325
  • PDF下载量:  9
  • 被引次数: 16
出版历程
  • 收稿日期:  2012-06-08
  • 修回日期:  2012-10-15
  • 刊出日期:  2014-01-25

目录

/

返回文章
返回