船舶与海洋平台结构冰载荷的高性能扩展多面体离散元方法

船舶与海洋平台结构的冰载荷是寒区海洋工程结构物设计中的关键参数,而离散元方法是有效计算结构冰载荷的重要手段. 本文采用基于闵可夫斯基和原理的扩展多面体离散元方法模拟船舶与海洋平台结构的相互作用过程. 其中,构造扩展多面体的近似包络函数并建立了基于优化模型的快速接触搜索算法;考虑单元间粘结作用的刚度软化过程建标识码元间的粘结-破碎模型. 同时,发展了 CPU-GPU 协同异构环境下的高性能并行算法. 为分析海冰与海洋结构作用中的冰载荷,采用ISO标准验证了扩展多面体离散元分析结构冰载荷的准确性. 采用离散元方法计算了船舶结构的冰载荷,研究了船舶结构表明的线载荷分布特点,并采用船舶结构冰阻力经验公式验证了计算结果的合理性. 采用离散元方法计算了平整冰区与多桩腿平台结构的相互作用,分析各桩腿上的冰载荷特点. 针对碎冰区的海冰管理过程,采用离散元方法分析了船舶结构绕行过程中的船舶和海洋平台结构冰载荷. 本文方法可有效应用于海洋结构冰载荷分析,能为极地船舶与海洋平台结构的设计和安全运行提供科学的分析手段.      

基于扩展多面体的离散单元法及其作用于圆桩的冰载荷计算

对于具有复杂几何形态的多面体单元,线性接触模型不能准确地计算不同接触模式下的作用力,且接触变形和作用力方向也不易判断.基于闵可夫斯基和(Minkowski sum)方法的扩展多面体单元能够准确描述非规则颗粒单元的几何形态,并可精确计算单元间的接触碰撞作用.该方法具有接触判断简单、计算效率高的特点.它将基本多面体和扩展球体相叠加以形成具有光滑棱边和角点的扩展多面体单元.考虑扩展多面体单元相互作用过程中角点、棱边和平面之间的不同接触模式,发展了相应的非线性黏弹性接触模型. 该接触模型将不同接触模型下的法向刚度统一表述为单元接触中接触点处等效曲率半径的函数;黏滞力和切向弹性力接触模型则借鉴球体单元非线性接触模型的处理方法. 为检验扩展多面体的可靠性,对碎冰区冰块对圆桩结构的冰载荷进行了离散元分析. 采用沃洛诺伊(Voronoi)切割算法获得了碎冰的初始随机分布状态,并考虑了海冰在运动过程中的海水浮力和拖曳力.计算表明该扩展多面体单元可描述海冰在海流拖曳下的运动过程以及圆桩结构的动冰力特性.在此基础上进一步分析了冰速和冰块尺寸对圆桩冰力的影响,并确定了冰力在圆桩上的分布规律. 最后,讨论了目前扩展多面体单元在计算冰载荷方面的局限性和改进方法.     

基于多介质、多尺度离散元方法的冰载荷数值冰水池

极地船舶与海洋工程结构冰载荷的确定是其结构抗冰设计、冰区安全运行和结构完整性管理的重要研究内容. 当前快速发展的高性能计算技术和多介质、多尺度数值方法为准确、高效地计算结构冰载荷提供了有效的途径, 其中以离散元方法为代表的数值方法取得了出色的研究成果. 为此, 本文针对目前极地船舶与海洋工程结构对冰载荷及力学响应的工程需求, 同时考虑国内外对海冰、工程结构与流体相互耦合的多介质、多尺度数值方法研究现状, 对极地船舶与海洋工程数值冰水池的概念、框架、开发技术以及基于离散元方法的软件实现与工程应用进行了论述. 数值冰水池在船舶与海洋工程结构冰载荷确定方面具有可靠性、经济性、快速性、扩展性和情景化等显著优势. 本文工作借鉴数值水池的研究思路, 以典型船舶和海洋平台结构冰载荷及结构力学响应的离散元计算为例, 探讨了数值冰水池研究的可行性和工程应用前景, 阐述其与理论分析、现场测量和模型试验研究相结合的必要性. 以上研究有益于中国在极地船舶与海洋工程领域形成具有独立知识产权的数值计算分析平台, 对中国极地海洋强国的战略实施具有很好的启发和指导意义.      

船舶在碎冰区航行的离散元模型及冰载荷分析

采用离散元模型对碎冰区浮冰与船舶结构的相互作用进行了数值研究。碎冰由三维圆盘单元构成,并考虑其在海流作用下的浮力、拖曳力和附加质量。船体结构由一系列三角形单元组合构造。通过海冰与船体单元间的接触判断和接触力计算,确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰况（冰速、冰厚、冰块尺寸和密集度）以及航速条件下,海冰对船体的动力作用过程进行了数值分析,对比分析了以上因素对船体冰载荷的影响,可为冰区船舶的安全运行和结构设计提供一定的借鉴作用。     

平整冰和碎冰区单桩锥体风电基础结构冰载荷离散元分析

在渤海和北黄海海域,冰载荷是影响风电基础结构振动响应及疲劳寿命的关键因素.对于单桩风电柔性结构而言,在水线处安装抗冰锥可有效降低冰载荷,保障海上风电安全运行.为分析海冰对锥体风电结构的冰载荷,基于扩展多面体离散元方法模拟了平整冰和碎冰对锥体风电结构的作用过程,分析了不同锥径和海冰密集度下的锥体风电结构冰载荷及冰激结构振动响应.计算结果表明,锥体结构冰载荷具有显著的随机性和周期性,平整冰对锥体结构冰载荷要明显高于碎冰;平整冰对上锥体的冰载荷大于下锥体冰载荷,而碎冰对上锥的冰载荷却小于下锥体冰载荷;所提方法可有效模拟海冰对锥体风电结构的作用力.     

极地船舶冰区操纵性能的六自由度离散元分析

极地船舶操纵破冰性能是破冰船设计建造过程中需重点考虑的问题。为分析极地船舶操纵性能,本文发展了冰区船舶六自由度操纵破冰运动模型,采用扩展多面体离散元模拟海冰,舵桨操纵模型提供破冰力。开展平整冰区定常直航模拟计算,并与Lindqvist船体冰阻力经验公式展开对比验证;开展雪龙号敞水35°舵角操纵性仿真,并与实船试航结果进行对比。在此基础上,开展冰厚及舵角影响下的船体结构冰载荷及破冰轨迹计算,模拟得到操纵破冰航行中船体垂荡、横摇及纵摇运动时程;最后分析操纵破冰船体线载荷分布与直航破冰下的差异。     

冰脊压剪试验及其对直立结构冰载荷的离散元分析

为更准确地分析冰脊与海洋结构的相互作用,本文提出一种基于粘结球体单元的冰脊离散元方法,其中龙骨与脊帆由不同尺寸的单元随机排列而成,整体呈各向同性.采用该方法模拟了一系列龙骨压剪试验,结果表明此方法能较好地表现龙骨的力学特性,且弹性模量、粘聚强度、内摩擦角等计算参数不仅影响龙骨的压剪强度,还影响其破坏模式.最后计算了冰脊...     

基于离散元-有限元模型的冰激锥体海洋平台结构振动分析

冰荷载是影响海洋平台结构及其上部设备作业安全的重要因素。针对海洋平台与海冰相互作用的动力过程,采用具有粘接-破碎性能的离散元方法(DEM)对海冰的破碎特性进行分析,采用有限元方法(FEM)对海洋平台结构的动力响应进行计算;同时考虑海冰与海洋平台结构的耦合作用,将海冰的离散元方法与海洋平台结构的有限元方法相结合,建立了冰激海洋平台结构振动的DEM-FEM模型,并由此计算海冰作用下海洋平台结构的振动响应以及冰荷载特性。     

冰区航行中船舶结构冰载荷的现场测量与反演方法研究进展

冰载荷是影响极地船舶航行安全的重要环境因素,而对船舶结构的现场监测是获取冰载荷的可靠途径.鉴于船-冰相互作用的复杂性, 目前还难以直接测量冰载荷,一般通过结构应变、六自由度运动参数等船舶结构局部或总体响应的测量数据间接反演冰载荷.根据冰载荷的作用范围,本文将船舶结构冰载荷现场监测划分为局部冰载荷现场监测与总体冰载荷现场监测两大类.对国内外18艘极地船舶冰载荷现场测量试验的开展时间、试验海域、测量方案等信息进行了系统的总结和分析.从基本原理、适用范围、应用现状和发展前景等方面全面地介绍了船舶结构冰载荷反演的影响系数矩阵法、支持向量机法、格林函数法、运动参数法和功能关系法,并重点分析了\"MV Timofey Guzhenko\"极地穿梭油轮与\"IBRV Araon\"破冰考察船的冰载荷测量结果.在此基础上对船体局部冰压、冰力峰值、冰载荷概率分布和冰激振动加速度等相关研究进展进行了深入的讨论.最后从测量技术、反演方法、冰载荷特性等方面剖析了当前船舶结构冰载荷现场监测中存在的问题,并探讨了相应的研究方向.本文对国内外极地船舶冰载荷现场测量与反演方法的论述可为后续研究与工程应用提供科学参考,从而更好地促进我国极地船舶的抗冰结构设计与冰区航行技术的发展.      

水位变化对正倒锥体冰载荷影响的离散元分析

在海冰与锥体海洋结构的相互作用过程中,潮汐水位变化时海冰作用于锥体结构的位置改变对冰载荷具有显著影响.本文采用具有粘结破碎功能的离散元方法计算海冰与锥体作用的破坏过程.同时考虑海冰上下表面温度差异对海冰强度的影响,将离散元计算冰载荷及海冰破坏模式与渤海现场实测数据进行对比验证.离散元结果表明,海冰与正锥和倒锥碰撞时均发生弯曲破坏,且冰载荷均随水线处锥径的增大而增大.在水线处锥径相同的情况下,正锥冰载荷大于倒锥冰载荷,而正锥作用下海冰的断裂长度则较小.基于离散元计算结果和渤海现场观测资料分析了海冰与正锥、倒锥作用时冰载荷和断裂长度差异的主要原因.海冰与正倒锥交界线处作用时,一般发生弯曲破坏.当冰层中心高度与正倒锥交界线的高度相同时,海冰才会发生局部挤压破碎,但冰荷载并没有明显升高.由此可见,倒锥体结构可有效降低冰载荷从而具有较好的抗冰性能.以上研究表明离散元方法可确定海冰与锥体结构作用时的海冰破碎规律和冰载荷特性,为海洋工程结构的抗冰设计提供参考依据.     

冰块尺寸对碎冰区船舶结构冰阻力影响的离散元分析

船舶在碎冰区航行时的冰阻力会影响船舶的安全性和航行性能。本文以雪龙号破冰船为研究对象,采用三维离散元方法对碎冰环境下船体结构冰阻力进行数值模拟。基于随机生成算法对不规则形状的碎冰区进行参数化建模以形成复杂的碎冰初始场,同时考虑海水对船体结构的作用力建立船体碎冰阻力的计算模型,并将离散元计算结果与碎冰区DuBrovin冰阻力计算公式进行对比。在此基础上将冰块尺寸进行无量纲化,研究不同冰厚和航速下冰块尺寸对碎冰阻力的影响规律。以平整冰区船体冰阻力的Lindqvist公式为参照,得到碎冰向平整冰转化的临界尺寸。该浮冰临界尺寸可作为识别碎冰区与平整冰的判据。本文研究工作可为船舶在极地冰区的安全航行提供一定的参考。     

海冰与自升式海洋平台相互作用GPU离散元模拟

在海冰与自升式海洋平台结构的相互作用过程中,冰载荷是影响平台结构振动响应和疲劳寿命的重要因素. 采用具有粘接-破碎效应的离散元模型,可对海冰与自升式海洋平台结构作用中的海冰破碎特征及相应冰载荷进行数值分析. 针对自升式海洋平台的多桩腿结构特性及其冰载荷离散元分析的大规模计算需求,建立了基于GPU 的并行算法并开发了相应的计算程序. 为实现离散元分析的高效计算,采用网格排序方法创建单元邻居列表,以快速确定海冰单元间及其与平台结构间的接触模式和作用力. 此外,还发展了球体单元与圆柱形结构在不同接触形式下的计算模型. 为检验该离散元模型的有效性,对渤海锥体海洋平台结构的作用过程进行了计算,并与现场实测冰力数据进行了对比验证. 在此基础上对多桩腿自升式平台结构的冰载荷进行了离散元分析,获得了海冰的破坏特性,确定了不同桩腿上的冰力时程. 该模型可进一步应用于不同类型海洋结构的冰载荷分析,为冰区海洋平台的结构设计和现役平台结构的疲劳分析提供参考依据.     

基于Voronoi切割算法的碎冰区构造及离散元分析

离散元方法广泛应用于海冰,特别是碎冰区的动力过程及其对海洋结构作用过程的数值模拟。为构造碎冰区中的冰块几何特性,基于二维Voronoi图方法对计算域进行随机切割以生成碎冰区中冰块的几何形态,并采用球体单元对每个碎冰块单元进行填充,从而确定碎冰区的初始分布场。在采用Voronoi图进行碎冰区构造时,可对冰块尺寸、几何形态和密集度等海冰参数进行设定。为确定冰块的不同几何规则度,综合采用排斥法和扰动法以定量地控制碎冰块几何形态从完全随机分布到规则分布的连续变换。为分析不同几何规则度下碎冰块的几何特性概率分布规律,对计算域内冰块的面积和边数等参数进行统计分析,从而可更合理地参数化控制初始冰场中碎冰块的几何特性。在此基础上,本文基于粘接-破碎的球体离散元方法对不同冰况下锥体结构的冰荷载进行了数值计算,讨论分析了碎冰区的海冰密集度、冰块面积和几何规则度对冰载荷的影响。     

船-冰作用中海冰典型破坏模式的离散元数值模拟

在破冰船破冰过程中,冰排主要表现为挤压与弯曲两种破坏模式.本文基于黏结离散单元法对船-冰作用中的这两种典型破坏模式进行数值模拟.海冰离散元数值试样采用随机排布方式生成,采用单轴压缩试验与三点弯曲试验相结合的方式标定模型中的细观参数.将船-冰碰撞中的挤压和弯曲作用方式简化为直立或倾斜 平板与海冰的作用模式,构建挤压与弯曲...     

爆炸载荷下海洋平台夹层板结构研究

海洋平台在服役期间,常会遭到油气泄漏而引起的爆炸冲击破坏.夹层板结构是近代发展起来的一种比较先进的结构形式,具有比强度高、比刚度大等特点.以某海洋平台燃油舱围壁为研究对象,设计了三种夹层板结构,分别为U型夹层板、V型夹层板和蜂窝型夹层板.基于此,使用MSC.Dytran对不同舱壁结构在爆炸载荷下的动态响应进行数值仿真计算,并与传统平板的计算结果进行比较分析,从而获得防爆效果较好的夹层板防爆舱壁结构.为模拟海洋平台燃油舱油气泄漏导致爆炸产生的爆炸冲击波对平台多舱室结构的破坏作用,采用了能够考虑耦合面破裂的多欧拉域的流固耦合算法.     

极地海冰动力过程的多面体离散元方法

极地海冰影响局部海域的通航条件乃至地质尺度上的气候预测,由于基于连续介质理论的海冰数值模式如网格方法难以应对浮冰尺度的模拟,近年来离散元方法越来越多用于海冰动力学领域。本文使用基于表面网格划分的多面体描述浮冰,通过使用GJK(Gilbert-Johnson-Keerthi)算法以进行初步接触判断并提高计算效率,采用基于能量守恒接触理论确定单元间接触力,由此建立了极地海冰多面体离散元方法。将该方法应用于极地海冰动力过程模拟,并考虑Hopkins提出的冰脊形成导致的海冰间塑性接触力。首先通过两个浮冰间碰撞模拟以验证接触力计算的正确性,再通过对单向均匀风场下浮冰挤压的模拟验证冰脊生成的模拟能力,最后研究冻结强度参数对剪切风场下冰盖断裂和冰脊生成的影响。结果表明,极地海冰多面体离散元方法能够合理地反映地质尺度海冰的动力过程,并且可有效地模拟冰盖线性运动引起的冰隙和冰脊等动力特征。     

碎冰区海冰与船舶结构相互作用的离散元分析

采用离散元模型对碎冰区海冰与船舶结构的相互作用进行了数值模拟。碎冰由三维圆盘单元构成,并考虑了其在海流作用下的浮力、拖曳力、附加质量;采用单元组合的方式构造了船体的简易模型。通过海冰与船体单元间的接触判断和接触力计算确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰速、冰厚、冰块尺寸、密集度条件下海冰与船体的动力作用过程进行了数值计算,并对比分析了不同海冰参数下冰对船体作用力的幅值。结果表明：冰对船体的作用力随冰厚、流速、密集度、冰块尺寸的增加而增大。本文工作可为冰区船舶的安全运行和结构设计提供参考。     

不同载荷频率下有砟铁路道床累积沉降特性的离散元分析

为研究有砟道床在高频载荷作用下的累积沉降规律,建立工程尺度下枕木-有砟道床的离散元计算模型。考虑道砟的非规则真实几何形态,利用不可破碎的簇颗粒模型近似构造。在10 Hz～100 Hz载荷频率区间内,对比分析枕木位移以及道砟颗粒在垂向和横向上的运动分布规律,研究分析道床在循环正弦载荷作用下的动力特性。数值结果表明,在载荷频率不超过60 Hz时,道床累积沉降随着载荷频率的提高而逐渐增大;当载荷频率超过60 Hz时,道床累积沉降随着载荷频率的提高而急剧增大,并在载荷频率为80 Hz时达到最大值,而后道床累积沉降量随着频率的提高反而逐渐减小。此外,分析对比不同载荷频率下道床速度场分布及典型位置道砟的加速度等物理量,发现80 Hz的载荷频率下道床离散元模型的振动最强烈。     

船舶与海洋平台撞击的荷载模拟

以某导管架结构的X斜撑在导管架安装阶段受到大吨位起重铺管船撞击后的损伤检测结果为依据。对碰撞过程进行了数值模拟,按照构件模型从简单到复杂的模拟思路得到了船舶对海洋平台结构的撞击力。建立了船舶与海洋平台碰撞系统的力学模型,对导管架平台受损构件进行了等效静力强度计算分析。为了考虑在船舶冲撞下的整体平台结构动力效应的影响,在分析时考虑了碰撞过程的水动力效应,桩-土-平台结构的相互作用,采用非线性弹簧来模拟受损构件局部凹陷损伤,在船舶以不同的速度与平台发生碰撞的情况下,模拟出平台受撞击构件的永久凹陷变形与碰撞接触时间的曲线．反演分析计算了船舶对平台结构碰撞的最大撞击力。为确定受到撞击后导管架平台结构构件和管节点的损伤程度,提出合理可行的修理加固方案提供了分析依据。     

碎冰区海冰与船舶结构相互作用的离散元分析

采用离散元模型对碎冰区海冰与船舶结构的相互作用进行了数值模拟。碎冰由三维圆盘单元构成,并考虑了其在海流作用下的浮力、拖曳力、附加质量;采用单元组合的方式构造了船体的简易模型。通过海冰与船体单元间的接触判断和接触力计算确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰速、冰厚、冰块尺寸、密集度条件下海冰与船体的动力作用过程进行了数值计算,并对比分析了不同海冰参数下冰对船体作用力的幅值。结果表明:冰对船体的作用力随冰厚、流速、密集度、冰块尺寸的增加而增大。本文工作可为冰区船舶的安全运行和结构设计提供参考。