西太平洋锚泊潜标异常沉降的仿真分析.docxVIP

西太平洋锚泊潜标异常沉降的仿真分析 现阶段，很难将数字信息有效地传输到水环境。 迄今为止,数学模型已经能够逐步模拟出极为近似锚泊线的柱状模型,不论是否采用弹性存在的方法,都能对其张力进行分析 深海潜标系统多为单点系泊 1 深度海洋系统的动态模拟原理 1.1 主浮体及浮体的装置 锚泊系统实际全长为5 870 m,布放水深6 070 m。该潜标系统主要包括1个主浮体、7组玻璃浮球组、2组温深盐链、深海海流观测单元、1组声学释放单元以及1组锚系单元,所有单元通过10 mm凯夫拉缆绳串联布放。主要搭载的设备如下:主浮体搭载2台75 k Hz声学多普勒剖面测量仪(ADCP),其流速观测范围为水面至水下1 000 m;在温盐链的设计中,水深200 m至500 m均匀布放四台温盐深测量仪(SBE37CTD),水深250 m至1 000 m布放10台温仪(SBE56);深海海流情况通过位于1 800 m、3 000 m和5 800 m的3个单点海流计观测深海的流速结构特征。潜标系统结构示意图如图1所示。 1.2 缆系缆系的流场—潜标动力学分析方法 当整个锚系系统在有流海域布放时,水流冲击力、重力锚阻力、锚系重力和浮力会使整个锚系系统进行动态飘移,整个缆系必将在三维空间中沿两个水平的方向横摇ΔX和ΔY,在铅直方向上产生一定的沉降ΔZ 其中,ρ为流体密度,C 缆体分析则是将柔性体绳分成有限个刚性有限元,使整个模拟系统形成以多体动力学理论为基础的集中质量–弹簧模型,其锚系缆绳上第i个质量的运动方程如式(2)所示 其中,m 1.3 缆绳受力分析 模拟仿真自容式潜标系统在水下工作时,将对系泊与环境进行如下假设,以便简化分析: 1)缆绳模型由n个微分缆绳单元组成,缆绳单元为刚直杆体,两个微分单元之间进行铰接,传递受力。同时由于n足够大,铰接缆绳模型的净重以及张力弯曲性质等都无限接近实际缆绳模型。 2)在全部深度上,海流方向都需分为u和v两个水平方向和一个垂直方向w,作用在任何一个元素上的三维模型方向一致,流体阻力作用在锚泊缆绳与水流速度向量的平面内。 3)假定系泊连接缆绳为柔性结构,不能传递弯矩。 4)永远假定锚系底部自由度为0。 假设中n可以无限大以达成精密度极高的模拟计算,但由于模型边界条件的理想化、回收数据量有限、忽略柔性体传递弯矩、水下震荡现象的作用等,可能导致模拟计算结果与实际情况存在一定差异。 根据假设绘制出第n个缆绳单元的受力分析图如图3所示,列出受力平衡姿态下计算公式。 式中,T 式中,C 将F 1.4 缆绳的精密度 据玻璃浮球和主浮体位置将锚系缆绳分成8段进行等比例建模,由于不同深度存在不同流速,进而对于缆绳的精密度划分要求有明显区别,3 000 m以下海洋潜流影响极小的缆系受力单一且稳定,适当降低精密度可减少结构的复杂性。其中凯夫拉缆绳的参考水下密度为1.05 g/cm 缆绳建模完成后在结构相应位置固定等比例大小的CTD、单点海流计、ORE等锚系悬挂设备,并施加同等重力。 2 影响潜标流场压力的因素 2.1 垂向海流对锚泊缆系的影响 潜标锚系系统于2017年4月布放于日本东南部西太平洋(30°N,146°E),于2018年6月进行回收。潜标回收后获得ADCP以及SBE37CTD,提取设备升降情况和海洋要素,搭载双75 k Hz ADCP的主浮体设备的深度变化情况,经过正交分解后的U、V和W三种方向时间-深度-海流关系图如图4、图5、图6。 仅通过本次对垂向海流数据的处理可知布放期间海流方向铅直向下,依据公式(3)—公式(9)计算可知,在0～1 000 m内垂向海流对锚泊缆系的下压流阻范围为0～10 N,由于整个锚泊系统总净浮力约6 600 N,故垂向海流对于缆绳整体姿态影响程度较小。 根据图中主浮体深度曲线可之,4月19日锚泊系统下放至姿态稳定后,主浮体位置与预设水下400 m基本相同。锚泊系统布放期间有3个沉降峰值,分别为(1)2017年7月5日下潜至水下577 m,实际下沉177 m,受东北方向海流影响;(2)2017年8月28日下潜至水下635 m,实际下沉235 m,受西南方向海流影响;(3)2018年4月27日下潜至水下829 m,实际下沉429 m,受西南方向海流影响。 2.2 海水密度对流阻情况 全布放时间段200 m至500 m SBE37CTD解析数据结果图如图7所示,由盐度和水温与深度对比曲线可以得知环境水温均低于17℃,代入相应海水盐度公式。 其中,S为盐度;?为海水比重;T为水温。计算可得200 m至500 m系泊大深度沉降时海水密度变化范围为1 023.5 kg/m 受海水密度影响的因素包括:a)根据公式(1)、(3)、(4)可知,海水密度会对流阻情况产生影响,流阻系数取1.2时,根据实际数据可知在2018