基于离散元方法的竖直群桩冰荷载分析

在高纬度海域或高寒地区的江河中,碎冰对海洋平台桩腿和桥墩等直立群桩结构产生很大的冲击作用,并在一定条件下造成不同程度的破坏。依据碎冰在自然条件下的离散分布特性,采用Voronoi切割算法构造随机分布和具有非规则几何形态的碎冰初始状态,在此基础上采用若干个具有黏结-破碎功能的球形颗粒构造具有一定厚度、尺寸和形状等物理参数的河冰离散单元模型,将竖直群桩简化为刚体结构,从而对碎冰的运移及其对群桩的冲击作用进行数值模拟。在此基础上可确定碎冰撞击下的群桩冰荷载,并分析冰块尺寸、冰速和冰厚对桩结构所受冰荷载的影响。     

海冰与海洋平台碰撞分析中的可视化方法

针对基于离散元方法的海冰与海洋平台结构碰撞分析系统,对运行于CPU GPU异构高性能计算平台上的大规模粒子模拟进行可视化方法的研究.使用OpenGL完成不同规模的平整冰、浮冰、冰脊和不规则形态海冰的呈现,不同形状海洋平台结构与受力网格的绘制,海冰粒子的运动轨迹与速度的表达,海浪效果的模拟以及碰撞过程的动画演示等.通过定义合理的交互模式与接口,在一定程度上融合GPU加速的离散元计算与后处理显示,基于多进程管道通信、多线程并行输出等实现在GPU加速计算的同时实时显示粒子模拟计算结果.该方法在帮助研究者随时掌握程序执行状况的同时,大幅减少后续传输、处理和存储的数据.     

冰区航行船舶冰阻力及六自由度运动响应的离散元分析EI北大核心CSCD

冰阻力是冰区船舶航行过程中的重要影响因素,对船舶航行性能及航行安全带来严峻挑战。该文采用扩展多面体离散元方法建立冰船动力作用过程的三维离散元方法,基于船舶推进功率与推进力间近似关系,将螺旋桨推力、冰载荷、舵力及水动力等载荷分开考虑,开展恒功率破冰船六自由度非线性操纵直航冰阻力及运动响应计算。为验证该离散元方法的可靠性,对比了碎冰区DUBROVIN及平整冰区LINDQVIST冰阻力公式计算结果。基于船舶在不同主机推进功率及复杂冰况下的直航冰阻力计算结果,展开对船舶破冰航速的影响因素分析。在此基础上,给出了船舶垂荡、横摇、纵摇等冰区船舶运动响应模拟结果和极地船舶运动响应预报结果,最后对碎冰及平整冰区船舶冰阻力及运动响应展开了对比分析。该方法可有效计算船舶冰阻力,其模拟结果可为冰区船舶运动响应及航行安全预警提供重要参考。     

船体结构冰载荷反演方法及试验验证

冰载荷是影响船舶冰区航行期间结构安全的重要环境载荷。船舶的冰压监测通常采用应变传感器,合理地布放传感器是识别冰载荷的基础。通过对比船体外板结构试验中的冲击载荷和不同测试方案下的应变信号,确定了最佳应变传感器布放方案;采用Green核函数方法建立了船体外板结构应变冲击载荷间的响应关系,并对采集信号在噪声影响下反演的不适定性进行了分析;采用Tikhonov正则化方法克服了载荷反演过程中出现的数值不稳定问题;最后将试验中的响应用到载荷识别分析中,反演的载荷可以较为准确地反映冲击载荷的时域特征并且载荷识别精度良好。     

返回舱着陆姿态和动力响应的DEM-FEM耦合分析

针对返回舱安全着陆问题,采用陆面颗粒离散元模型（discrete element model, DEM）与返回舱壳体有限元模型（finite element method, FEM）相耦合的方法,对返回舱回收着陆过程的动力特性进行数值分析,并与试验结果进行对比,验证计算结果的正确性。分析讨论垂直和倾斜2种着陆姿态的返回舱结构应力分布,从冲击能量的角度讨论陆面颗粒对其缓冲作用。对影响返回舱回弹响应的主要因素进行分析并提出相应的改进方法。结果认为：陆面颗粒是返回舱着陆冲击能量吸收的主体,且倾斜着陆时陆面颗粒吸能占比更加显著;当陆面颗粒之间的摩擦因数较小或返回舱质量较大时,返回舱着陆的回弹速度相对较小。     

失效测点影响下极地船舶结构冰载荷的有效识别方法

冰载荷是极地船舶在冰区航行时受到的一种极端环境载荷,而对船体结构的现场监测是研究冰载荷的重要途径。一般通过在船体结构上安装应变传感器直接测量冰载荷引起的应变响应;采用影响系数矩阵法根据应变响应反演识别冰载荷。然而当应变传感器无法正常工作时,识别结果的准确性将难以得到保障。通过对典型极地船舶舷侧板架结构的有限元分析,研究了测点失效对识别结果的影响。基于对“天恩号”多用途冰级船实测应变数据以及典型极地船舶舷侧板架结构有限元应变数据的深入分析,确定了测点应变的空间分布规律,并进一步提出了失效测点影响下基于最小二乘拟合的冰载荷识别方法。在此基础上,较为准确地识别出7种典型工况下的冰载荷,大幅降低了识别误差,最终验证了该方法的有效性。     

渤海海冰的重叠长度分析

基于渤海海冰现场观测资料,分析了重叠冰的基本特征;综合考虑海冰重叠过程中的屈曲破坏和挤压破坏模式,建立了计算一维重叠冰和二维指状冰重叠长度的力学模型,并依此对不同冰厚条件下的海冰重叠长度进行了计算。结果表明：海冰重叠长度主要由海冰屈曲破坏条件来决定;在相同冰厚条件下,指状冰的重叠长度约为一维重叠冰的1／2左右。     

渤海冰脊分析及其设计参数

在冰脊实测资料的基础上,分析了压力冰脊的几何特征;综合考虑冰脊形成过程中的弯曲破坏和屈曲破坏模式,建立了冰脊帆高的力学模型;给出了渤海不同海区不同重视期的冰脊设计参数。以上工作可用于渤海海水数值模拟以及近海工程设计中的冰脊参数计算。     

渤海海冰动力学的质点网格法数值模拟

为准确计算渤海海冰的动力过程，本文将质点网格方法(Particle—in—cell，简称PIC)，引入到渤海海冰数值模拟中。该方法首先在欧拉坐标系下对海冰动力方程和连续方程进行差分求解，并插值出各网格内海冰质点的速度和密集度，然后在拉格朗日坐标系下对海冰质点的位移和各网格的平均冰厚进行确定。PIC方法将欧拉法与拉格朗日法相结合，避免了欧拉坐标下有限差分法的数值扩散和拉格朗日坐标下光滑质点流体动力学计算量大的缺点。在对渤海海冰动力过程的数值模拟中采用了Hibler的粘塑性本构模型，并考虑了海冰热力作用过程。利用PIC方法对辽东湾海冰进行了48h数值模拟，结果表明：该方法可成功地处理海冰流变过程，精确计算出冰缘线位置和海冰分布状况，对海冰厚度和密集度的计算精度均优于有限差分法。PIC方法在精确模拟海冰的重叠和堆积过程，以及为冰区油气作业提供详实的海冰信息方面都有良好的应用前景。     

面向船体结构动冰载荷监测的线性形函数识别方法研究EI北大核心CSCD

船体结构冰载荷监测系统中的载荷识别算法是对船舶在冰区安全运行进行实时评估的关键环节,而常用的影响系数矩阵法尚未考虑冰载荷的动力效应。在时域内将动冰载荷离散为诸多时间元,在每个时间元内利用线性形函数的组合形式逼近构造动冰载荷;通过结构的动力响应分析得到形函数的响应矩阵,并以此建立船体结构冰载荷识别的正问题。利用共轭梯度最小二乘迭代型正则化算法和终止迭代准则对冰载荷识别问题中的不适定性进行控制。该数值算例和试验验证均表明,该冰载荷识别模型具有良好的识别精度、求解速度及稳定性。