储箱内液体晃动动力学分析及防晃结构优化.docx

该【储箱内液体晃动动力学分析及防晃结构优化 】是由【老狐狸】上传分享，beplayapp体育下载一共【10】页，该beplayapp体育下载可以免费在线阅读，需要了解更多关于【储箱内液体晃动动力学分析及防晃结构优化 】的内容，可以使用beplayapp体育下载的站内搜索功能，选择自己适合的beplayapp体育下载，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此beplayapp体育下载到您的设备，方便您编辑和打印。储箱内液体晃动动力学分析及防晃构造优化刘富;童明波;陈建平【摘要】基于光滑粒子流体动力学(SPH)方法对无阻尼板和装有阻尼板的矩形储箱在加速度鼓励作用下,储箱内液体的晃动与冲击进展了三维数值模拟,将测试点的计算压力及液体晃动模拟与试验进展了比照,:,基于iSIGHT优化平台,通过拉丁方试验设计方法设计出样本点,构造出液体晃动的二次响应面模型,承受序列二次规划方法对防晃挡板的安装位置及尺寸进展优化,在保证液体晃动重心满足约束条件的根底上,防晃板的重量最轻.%SPHexecutes3-urateandsuperioratsolvingstrongnon-,samplepointsaredesignedwithLatindesigningmethod,.【期刊名称】《计算机应用与软件》【年(卷),期】2023(028)012【总页数】4页(P202-205)【关键词】光滑粒子流体动力学;储箱;晃动;重心;优化【作者】刘富;童明波;陈建平【作者单位】南京航空航天大学飞行器先进设计技术国防重点学科试验室江苏南京210016;南京航空航天大学飞行器先进设计技术国防重点学科试验室江苏南京210016;南京航空航天大学构造工程与力学系江苏南京210016【正文语种】中文【中图分类】TP391引言储箱内带自由液面的液体晃动不仅会对贮箱壁产生冲击载荷,引起构造的疲乏破坏。对于载液比重大的载液系统,如航天器、大型民机以及油罐车等,在运动过程中,短时间内液体的大幅晃动引发的液体重心变化还会对这些载液系统的运动轨迹或姿态产生不利的影响。无网格粒子形式以及拉格朗日性质的光滑粒子流体动力学(SPH)方法[1]作为一种相对而言比较的数值计算方法,解决了传统的基于网格法的数值方法在求解液体大幅晃动问题时消灭的大变形导致网格发生畸变的问题。自Monaghan[2]首次将SPH方法应用于自由外表水波的流淌模拟以来,SPH方法开头渐渐应用在工程实际中[3,4],其计算精度、稳定性及自适用性都已经到达实际工程应用的允许范围。近年来,国内学者承受SPH方法对自由外表流淌的数值模拟取得了肯定的争论成果[5,6]。带自由液面储箱的防晃设计,主要是通过设计合理的储箱构造外形以及利用构造内部阻尼装置来提高晃动阻尼。对于外形固定的储箱构造,则主要通过试验方法或者数值方法争论防晃挡板的外形(环形、圆形及矩形等)、尺寸、安装位置及安装方式(横向或纵向)等参数对晃动阻尼的影响[7,8]。为了验证SPH方法求解液体大幅晃动问题的准确性与优越性,基于冲击碰撞动力学分析程序PAM-CRASH,承受SPH方法进展三维矩形储箱内液体晃动与冲击的数值模拟,并与文献中的试验进展比照。目前,国内外对于液体晃动过程中液体重心问题的争论较少,本文推导了液体晃动重心的表达式,基于iSIGHT优化平台,对矩形储箱内防晃挡板的安装位置及尺寸进展优化,在保证液体晃动重心满足约束条件的根底上,防晃挡板的重量最轻。SPH数值方法掌握方程连续介质流体动力学N-S方程可分别由式(1)和式(2)描述:式中,上标α和β表示坐标方向及求和重复指标。v表示粒子的速度矢量,σ为总应力张量,F为体力。通过应用SPH方法的核近似法和粒子近似法将N-S方程进行空间离散化,得到一系列与时间相关的常微分方程,表达式如式(3)和式(4)所示,最终通过对时间的积分进展求解。式中,vij表示粒子i和粒子j的相对速度,Wij为粒子j对粒子i产生影响的光滑函数。SPH方法关键问题为了准确模拟流体动力学的问题,防止在冲击域内的求解结果产生非物理震荡,必需对SPH算法进展一些特别的处理,如添加使用最为广泛的Monaghan型的人工粘度∏ij[9]。对于不行压缩流,可以通过利用Monaghan提出的人工压缩率[10]将其视为可压缩流。SPH数值计算不考虑储箱内液体外表所受张力,分别对Craig试验[11]中无阻尼板和装有阻尼板的矩形储箱在平动加速度鼓励作用下,储箱内液体的晃动与冲击进展数值计算,并与试验结果及CFD计算结果[11]进展比照。试验中在矩形储箱的底面和侧壁布置两个压力测试点P1和P2,储箱的长、宽、高分别为500mm、400mm、400mm,储箱的主视图及压力测试点布置如图1所示。液深为储箱高度一半,储箱固定于一辆轻型卡车上,车由静止加速再减速至静止状态,小车的实测加速度曲线及滤波后用于数值计算的曲线如图2所示,表1给出了主要的材料属性参数。图1矩形储箱示意图在有无阻尼板两个工况的计算中,粒子数一样,均到达175760,。光滑核函数取Monaghan[12]提出的分段三次样条函数,三次样条光滑函数具有计算稳定且精度高的优点,大多文献中的SPH方法均承受此光滑函数。人工粘性系数α和β取值均为1,选取Monaghan提出的人工压缩率方法。计算在HP工作站进展,承受6个CPU并行计算,均耗时28小时左右。表1材料属性属性名称数值矩形贮箱及阻尼板密度/(kg/m3)7830矩形贮箱及阻尼板杨氏模量/(Pa)×/(kg/m3)998液体体积模量/(Pa)×109液体粘性系数/(N·s/m2),两种数值计算方法对晃动的模拟与试验的比照,SPH计算较真实地模拟出了试验中自由液面随着储箱加速度变化,从安静到大幅晃动的过程,并且从图中可以看出本文SPH方法对晃动的模拟与试验结果吻合度在某些时刻比CFD方法要好一些。图3液体晃动的数值模拟与试验比照图4给出了试验、CFD方法及SPH方法所得的测试点P1处的压力时间历程曲线。从图4(a)中可以看到5s之后测试点的压力所呈现的双峰特征,这是液体晃动过程中迟滞特性的一种典型表现。其中前一个峰值是由于液体对储箱壁的冲击而产生的高脉冲压力,时间短,具有比较明显的局部特性。这可以从图2中的加速度曲线分析出:5s-10s这段时间内加速度变化大,储箱壁主要承受水跃或者行进波晃动现象引起的瞬时冲击压力。后一个峰值也是由于冲击产生的脉冲压力,但是时间相对较长,而且冲击压力峰值变小,这主要是由于快速并且连续的液体压力作用在没有被完全漂浮的储箱壁,压力特性主要表现为连续性,而不像前一个峰值中所明显表现出的瞬时性。当给储箱加装阻尼板后,液体对储箱壁的冲击力度降低,由于阻尼板的面积较大,阻尼板的安装肯定程度上可以等效为增加了储箱的充液比,在大充液比状况下,迟滞现象消逝,因此图4(b)中没有明显的双峰特性,且压力峰值减小,晃动现象更多的表现为驻波运动。另外,在无阻尼板工况下,当加速度减至0并根本趋于稳定之后,P1处的压力值开头在初始静态压力值处四周振荡,并呈现出有阻尼的渐渐衰减的趋势;在有阻尼板工况下,P1处压力值则很快收敛于初始静态压力值。图4SPH数值计算与试验及CFD数值计算压力曲线比照通过将图4中的SPH计算结果与试验及CFD计算结果比照可以觉察,SPH方法计算所得到的曲线较CFD方法不管在变化趋势、压力作用周期还是压力峰值方面都与试验吻合更好。SPH方法所计算的零时刻的压力均为0,主要是由于粒子受到重力场的作用从零时刻开头,处于液面下的测试点,其四周粒子感受上面粒子的重力对这些粒子的作用需要一个短暂的时间,因此四周这些粒子对测试点产生的静压力也需要一个过程,随着时间的推动,压力测试点就会受到重力场作用下的液体静压作用以及液体冲击带来的动压。图4中的计算峰值稍低于试验所测峰值,这是由于前面所提及的对小车的加速度鼓励曲线进展了滤波,滤去了一局部波峰的原因。对于8s-10s时间段,两张图中的计算曲线与试验结果吻合均不是很抱负,从数值计算所得曲线的斜率可以看出压力变化比较突然,这可能是由于试验过程中低频加速度鼓励测量误差所致[11]。该时间段加速度测量值变化程度要比真实状况猛烈,导致数值计算中的液体冲击幅度较大,而在真实试验中,此时间段内储箱壁承受的压力主要为持续性的压力,因而压力曲线变化比较平缓。防晃挡板优化由于副油箱内的液体由一系列在空间运动的、附有质量的粒子模拟,因此液体在任一时刻t的重心可以通过如下表达式求得:式中,xt、yt和zt分别是t时刻的液体重心坐标,坐标原点为晃动中心处。xi、yi和zi分别为粒子i在该时刻的坐标,mi为粒子i的质量,n为粒子总数。图5为防晃挡板优化模型示意图。水平阻尼板关于Y轴对称,位于储箱高度1/2处,垂直挡板位于Y轴上。其主要设计参数为挡板的中心位置及挡板的宽度。液体在晃动过程中沿X轴方向的重心变化峰值XCG1和XCG2示意如图6所示,其中CG1表示重心点位于平衡位置左侧,即XCG1为负值,CG2表示重心点位于平衡位置右侧,XCG2为正值。图5防晃挡板优化模型图6晃动液体重心示意图基于iSIGHT优化平台,承受拉丁方试验设计方法设计出40组试验样本点,建立液体晃动的二次响应面近似模型,承受序列二次规划算法对防晃挡板的中心位置、尺寸进展优化。优化目标挡板重量最轻,即挡板的宽度之和U=U1+2U4最小。约束条件X方向的重心变化XCG1≥–,XCG2≤,U2-≥0,U2+-≤0,U3-≥0,U3+-≤0。单位:m。≤U1≤,≤U2≤,≤U3≤,≤U4≤,单位:m。建立的二次响应面模型,XCG1和XCG2的数学表达式分别如式(6)和式(7)所示,单位为m。XCG1二次响应面模型对真实样本点的靠近程度指数为:复相关系数R2=,RAdj2=。XCG2二次响应面模型对真实样本点的靠近程度指数为:复相关系数R2=,RAdj2=。从系数可以看出二次响应面多项式对样本点拟合很好。承受序列二次规划优化算法,各个设计变量及优化结果的迭代收敛曲线如图7所示,经过迭代求得的优化结果为:优化变量U1=,U2=,U3=,U4=;优化目标U=U1+2U4=;约束XCG1=–,XCG2=。承受SPH数值计算方法求得XCG1=–,XCG2=,%%,说明所构造响应面模型的计算结果与程序计算结果拟合很好。计算模型的初始值及优化值如表2所示。从表中可以看出在满足重心约束的条件下,优化后的防晃挡板重量比初始重量大大减轻。图7防晃挡板优化迭代收敛曲线表2防晃挡板优化前后参数值模型参数设计变量约束目标U1/mU2/mU3/mU4/mXCG1/mXCG2/mU/--,(XCG2)(XCG1)。在12s左右后,随着储箱运行速度的减慢,液体趋于安静,重心的变化渐渐收敛于0。图8矩形贮箱液体晃动X方向重心变化时间历程曲线3结语承受SPH方法,对矩形贮箱在平动加速度鼓励作用下自由液面的晃动与冲击进展了三维数值模拟,并基于iSIGHT优化平台对储箱内防晃挡板的安装位置和尺寸进展优化,可以得出如下结论:通过与文献中试验数据的比照,觉察计算所得到的压力曲线在变化趋势、压力作用周期和峰值方面都与试验吻合较好,SPH方法进展带自由液面的液体晃动计算是合理准确的。文中算例承受的SPH方法计算结果较CFD方法更接近于试验结果;整个晃动过程模拟也与试验相近,并且比较真实地模拟出了自由液面大幅晃动下承受有限差分法、有限元方法等难以模拟的碎波现象,显示了SPH方法在求解带有强非线性现象的液体大幅晃动问题方面的优越性。在一样的储箱外形和工况条件下,阻尼板的安装一方面可以显著转变液体的晃动特性,另一方面降低了测试点的压力,适宜的阻尼板安装对液体的防晃具有重要的作用。通过对防晃挡板的宽度和挡板中心安装位置进展优化,在满足液体晃动重心的约束下,挡板重量大大减轻。对于载液系统的防晃设计供给了很好的参考依据。参考文献[1]GingoldRA,:theoryandapplicationtonon-sphericalstars[J].MonthlyNoticesofRoyalAstronomicalSociety,1977,181:375-389.[2][J].putationalPhysics,1994,110(2):399-406.[3]VuystTD,VignjevicR,[J].InternationalJournalofImpactEngineering,2023,31(8):1054-1064.[4]ShaoSD,ChangMJ,GrahamDI,[J].CoastalEngineering,2023,53(9):723-735.[5]陈正云,朱仁庆,[J].船海工程,2023,37(2):44-47.[6]杜小弢,吴卫,龚凯,[J].水动力学争论与进展,2023,21(5):579-586.[7]BiswalKC,BhattacharyyaSK,-linearsloshinginpartiallyliquidfilledcontainerswithbaffles[J].,2023,68(3):317-337.[8]MalekiA,[J].JournalofSoundandVibration,2023,311(1-2):372-385.[9][J].putationalPhysics,1989,82(1):1-15.[10][J].putationalPhysics,1994,110(2):399-406.[11]CraigKJ,[J].StructuralandMultidisciplinaryOptimization,2023,33(1):71-87.[12]MonaghanJJ,[J].AstronomyandAstrophysics,1985,149(1):135-143.