结构物受水下爆兴波感化的数值分析.pdfVIP

数解决时域初值问题；从Stoker[9]对直接时域方法进行讨论，研究人员就开始了关于时域计算的某些 基本原理和计算方法的研究；Wehausent]o]重新表达初值问题，并应用含时间因子的时域格林函数导 出了相应的积分方程。 本文首先基于一个简单的气泡模型的假定，用四阶龙格．库塔法求解气泡半径控制方程从而得到 气泡运动方程，然后通过求解拉普拉斯方程得到气泡脉动引起的自由水面兴波情况，之后用求解出 的水波作为入射势，用边界元的方法分别求解波浪与固定结构物和漂浮结构物的作用问题。 1理论和数值模型 1．I最简单的气泡模型 由于时间尺度上的原因，水下爆炸兴起的水波不大可能是冲击波引起的，而是由气泡引起的。 S一25 即使这样，相对于自由水面波浪的固有周期(大约在o．5 s)气泡的持续时间仍然很短。因此，气 泡兴起的波浪的特点是相当于水域受到初始瞬息扰动而产生的。图2显示的是气泡产生的水波可能 过程。 。 o o Bubble Bubblecdapdns ～一d舶n翻她 Wave Wave produdng propels 图2气泡兴波的可能过程 真实情况的气泡在坍塌末期会有射流形成㈣，为非球对称运动而且气泡在脉动过程会有迁移运 动‘培l，气泡的运动过程十分复杂，不过本文讨论的是深水爆炸下气泡脉动引起的兴波，气泡可认为 处于无限域中瞄J，故可以假设气泡模型：‘ (1)气泡是球形的，中心位于空间一直不动； (2)气泡的埋深足够深，因此自由表面对气泡的影响不计； (3)只考虑气泡对自由表面的影响。 在这样的假设下，我们可以把气泡当成一个脉动的源。于是，这个问题就变为在脉动源强为Q(力 的作用下的瞬变波问题。这个问题分别在气泡动力学(可以求Q(f))和水波理论中得到解决。 根据气泡动力学，假设：流体理想，不可压缩；气泡的埋深为d。计算坐标系如图3所示。 则存在速度势饥满足： V2西。=0 流体区域内 (1) 一十—1 7lU■。trlHII— IZ- 孕：_dR 在气泡壁面上 (2) 一=一 lVo小专。 在无穷远 (3) 一p警一詈fV吼[2-pgd：e14万∥-2在气泡壁面 (4) 土g ● ·d ，， ● 6…_Q 图3计算模型及坐标系 其中尺是气泡的半径。将气泡的脉动用源汇来代替，即①铲Q(D／r，那么由条件(2)可得： Q=一妇2 (5) 将(5玳入方程(4)可以得到半径R的控制方程： 妇+23_k2=鲥+笠pf，k鱼R13，J (6) 公式(6)中，Po和凰是气泡初始时的压力和半径，它们是由炸药量来决定的。 1．2气泡半径控制方程 求解上述简单气泡模型可以归结为求解以下初值问题： j}=F(足应)=百gd+旦Rpfk鱼R13)。一三2壁R