三维矩形弹性液舱内液体晃荡数值模拟探究.pdfVIP

第42卷第5期 船海工程 V01．42No．5 2013年10月 SHIP＆0CEANENGINEERING 0ct．2013 DOI：10．3963／j．issn．1671-7953．2013．05．027 三维矩形弹性液舱内液体晃荡数值模拟研究 陈星，蒋梅荣 (大连理工大学海岸和近海工程国家重点实验室，辽宁大连116023) 摘要：为了解三维弹性液舱内液体晃荡情况，采用ADINA模拟不同激励频率下液舱内液体的晃荡，将 所得的结果与解析解进行对比，得到不同载液率下弹性液舱内液体晃荡特点，得到液舱载液率、外界激励等参 数变化对液体自由表面运动及晃荡荷载的影响。 关键词：液体晃荡；流固耦合；弹性液舱；ADINA 中图分类号：U674．13 文献标志码：A 文章编号：1671-7953(2013)05-0099-06 液体晃荡是极为复杂的一种流体运动形式， 制方程——流固耦合系统方程中，在单一时间步 它具有高度的非线性和随机性。在以往的研究 内，对流固耦合系统方程所有变量同时进行求解。 中，液体晃荡问题大部分的研究工作主要集中在 采用位移一压力格式描述流固耦合系统。 刚性液舱数值模拟方面¨引，然而实际工程中所使 流固耦合系统方程为 用的液舱都有一定的弹性，假定液舱是刚性结构， (1) FX三一墨，ds(xs’*0 忽略了液舱结构的变形，同时也忽略了液舱变形 oF。[x。，下，(墨)]1 与内部流体的相互作用。近年来有部分学者对弹 式中：卜耦合系统的解向量，X=(X，，x。)。 性液舱的液体晃荡问题进行了一些研究H引，对于 其中：x，，墨——定义在流体和结构节点上的解向 弹性液舱晃荡，必须研究晃荡过程中带自由液面 量。 流体的非线性动力学问题及其结构变形的相互耦 因此流固耦合面上结构的位移和流体应力分 合问题。 别表示为d。=d。(鼍)，T，=Tf(墨)。一，只分别是 流体与结构相应的有限元方程。 1 数值方法 对于单独的流体或结构模型来说，其线性方 假设矩形液舱为线弹性结构，液体为不可压 程如下。 缩的粘性流体。流固耦合接触面上，流体和结构 1)流体方程。 不直接脱离，即流体与结构在接触面上任意时刻 B，=一F；；Ff[群，Add：+(1一A。)d!。1] 法向位移都相等。 (2) 采用有限元法对流体域连续方程、运动方程、 2)结构方程。 质量守恒方程等进行离散求解；采用隐式积分法 B。=一F：兰F。[x：，A，T：+(1一A，)T；。1] 对结构域几何方程、运动方程、物理方程等进行积 (3) 分求解∞J。 将流体方程和结构方程放在一个方程组中 ADINA在求解流固耦合问题时有直接