动力外荷载作用下浮式板瞬态水弹性分析.pdfVIP

动力外荷载作用下浮式板 瞬态水弹性分析 邱流潮刘桦 (上海交通大学船舶海洋与建筑工程学院，上海200030) 摘要随着浮式平台用途的不断推广，弹性浮体结构的瞬态动力行为研究最近越来 越受到重视．这类结构水平方向的尺寸一般远大于厚度方向的尺寸而具有相对较小的刚 度，使得其水弹性研究比剐体运动的研究更加重要。为了浮式平台的安全运行，就必须 分析它们在各种荷载运动、飞机起降、以及导弹发射等引起的外荷载作用下的动力特 性。本研究建立了计算任意动力外荷载作用下弹性浮板的瞬态响应的时域有限元方法。 其中，结构和流体的运动都假定为小幅运动．几个数值算例验证了本文方法的有效性． 关键词浮式板，水弹性，瞬态响应，流体一结构相互作用，时域分析，有限元法 1前言 浮体结构(包括浮冰、防浪堤、以及各类浮式平台等)的水弹性研究成果具有重要 的应用价值。其中，浮式平台由于用途更加广泛，因而，浮式弹性结构的瞬态行为研究 最近越来越受到重视。这类结构水平方向的尺寸远远大于厚度方向的尺寸而具有相对较 小的刚度，使得其水弹性研究比刚体运动的研究更加重要。为了浮式平台本身及其搭载 设备能够安全运行，需要对它们在各种荷载运动、飞机起降、以及导弹发射等情况下的 动力特性进行分析。上述这些情况下浮体所遭遇的动力外荷载一般都具有冲击荷载或移 动荷载或者二者兼而有之的特性。因而，冲击荷载和移动荷载作用下浮板瞬态水弹性反 应研究具有重要的意义。 浮体的水弹性闯题是一个典型的流体～结构相互作用问题。由于这类闯题的复杂 性，只有极少数问题可获得闭合解析解，而大多数情况下都需要进行数值求解。与其它 数值方法相比，有限元方法由于具有处理任意几何形状和复杂边界条件的能力而被广泛 采用，本文在浮板的瞬态水弹性分析中，流体和结构均用有限单元进行离散。在海洋工 程中，由于海洋在水平方向是无限延伸的，因而在进行时域有限元分析时，需要首先对 流体求解域进行人工截断以获得有限的计算域．同时，为了使得流体一结构相互作用产 生的外行波不会在人工边界发生伪反射进入计算域而影响计算精度，就需要在人工边界 ·839- 施加适当的辐射边界条件。在众多的辐射边界条件中，最简单和使用最广泛的是 Sommerfeldtl】类型的边界条件，不过这类边界条件主要的问题很难正确计算外行波相速 度。特殊情况下，可以假定到达边界的主波的特性是可以识别的，譬如非耗散重力波。 Odanskit21给出了计算波相速度更一般的方法，不过这种方法要求边界附近的网格尺寸足 够小，否则将得到不正确的相速度．另外，还有一种方法就是在人一I：边界设计阻尼层，这 种方法最先是由IsraeliI”等人提出的，不过阻尼层的设计要求非常小心，特别是对亍三维 问题。本文将采用简单有效的重力波辐射边界条件[4】。 2数学公式 流体一结构相互作用系统如图1所示。弹性板的动力响应用挠度w0，¨0表示，流体 水动力问题的解以动水压力P．∽f)表征，位于弹性薄板下的流体(卢1)和位于弹性薄板外的 流体(扛2)由分界线J隔开。t表示时间。x=0．H：)表示空间坐标，z坐标轴垂直向上的方 向为正方向，：=0的平面定义为液体初始自由表面，坐标原点位于板的中心。通常情况 _卜，浮式平台具有较大的长宽比，同时水平方向的尺寸远远大于厚度方向的尺寸，所以 可以采用薄板进行模拟，弹性薄板的振动方程为(其中上标点表示对时间的导数)： DA2w+n^谛+pfgw+pl+q=0 其中，A=a2／c3x2+a2／易r，D=Eh3／020一，，2)]，仃为流体密度，h，占，岛，y分别表示 板的厚度，弹模，密度以及泊松比，g为重力加速度，g伍Ⅳf)表示作用在板上的时间相关 的载荷． 自由悬浮于液体表面的薄板应满足薄板边缘的弯矩和等效剪力都为零的条件，即 垡+。鸳：o， n!芝：o On。 西2 一8nasz 式中H和s分别表示法向和切向．另外，对于方板，在板的四角还必须满足如下的条 件： 盟：o (3) 舐