冰柱冲击问题的数值仿真分析.docxVIP

? ? 冰柱冲击问题的数值仿真分析 ? ? 王 超， 杨 波， 张 媛， 郭春雨， 叶礼裕 (哈尔滨工程大学 船舶工程学院，哈尔滨 150001) 冰的撞击过程和随之引起的动态载荷是关系极地科学领域和冰区科学问题的重要研究内容之一，例如极地船舶的设计和建造、极区科考探索以及资源开采、航空和星际空间任务过程遭遇环境.因此研究冰在撞击过程的破坏过程以及动态载荷变化，一方面能够为上述实际问题提供安全作业依据，另一方面有助于为各种有与冰发生碰撞可能的结构物提供规范设计指导. 目前研究冰冲击破坏过程有试验研究和数值研究方法.试验方法能够揭示冰真实撞击过程的破坏机理，然而冰的冲击试验并没有标准的规范章程，因此冰模型的制备具有不统一性且很难定量把控试验过程的对比工况参数.可靠的数值模型计算方法可以弥补试验中的缺陷，研究定量参数对冰冲击过程的影响，基于商用LS-DYNA软件中网格方法的冰冲击过程预报是一种最为常见的数值方法[1-4].然而软件库中的本构模型并不能完全描述海冰复杂的力学性能，后期学者们在用户自定义模块改进了冰的本构模型[5-7]，增加温度、应变率及孔隙率等因素的影响.有学者基于黏聚单元模型构建了海冰的均化弹塑性本构，采用LS-DYNA计算了海洋结构物与平整冰相互作用的冰载荷[8-9]，然而传统的网格方法在处理断裂问题方面仍旧具有一定的局限性.后来，无网格粒子法逐渐被应用于冰的冲击破坏问题，例如光滑流体动力学(SPH)[10]和近场动力学 (PD)[11]方法. 近场动力学方法是近二十年新兴的一种非局部、无网格粒子、积分形式的计算方法，在断裂问题域和非断裂问题域上都可以连续求解[12-14].该方法最早由Silling提出并建立推导了其控制方程，给出了基本的数值实现程序案例[15].PD方法主要有两种表现形式，经推导和比较，两种表达方式大同小异且具有相类似的计算精度[16].PD发展至今主要有键型、常规状态型和非常规状态型3种形式，目前已成功应用于材料的弹性、塑性、热力学等多种物理问题计算分析领域[17-18].然而PD应用于冰计算领域刚刚起步，多数文献采用基于键型的弹性模型来模拟冰的破坏特性，忽略了泊松比的真实状态[19-20].在键型PD中，物体中两点之间的力被假定为位移和这些点的初始位置的函数.由于忽略了其他键变形对材料点对相互作用的影响，该公式存在如下缺点：泊松比的限制(对于二维平面应力情况，将泊松比固定为1/3，对于平面应变和三维情况，将泊松比固定为1/4)以及仅经历体积应变的多孔材料的塑性建模能力[21]. 本文的主要目的是建立基于状态型PD方法的冰的弹性本构模型，消除键型对泊松比的限制，并研究和讨论冰本身物理参数、数值预报方法参数的敏感性特征.首先通过冰柱冲击试验验证了数值模型的准确性，分析了泊松比对冰体裂纹扩展的影响，然后进行了数值计算工况的变参数预报分析，即参数敏感性分析. 1 理论基础 1.1 常规状态型近场动力学 图1 物质点x与其临近物质点x′相互作用Fig.1 Interaction of a material point x with its neighboring points x′ 考虑到力矢量状态，基于OSB-PD的积分微分方程描述为[23] t′[x,t]〈x-x′〉]dVx′+b(x,t) (1) (2) (3) (4) (5) 式中：V′、V分别为材料点x和x′的体积. OSB-PD的最终表达式可以表示为 (6) 式中：a、b、d分别为PD专属参数；x′-x为变形前材料点相对位置；y′-y为变形后材料点相对位置；Λ为PD附加参数；θ和θ′分别为材料点x和x′处的膨胀变形；s为材料点之间变形伸长量(简称键). 1.2 刚性碰撞载荷 冰物质点和碰撞体有相对位移时，海冰物质点将与刚性物体表面发生接触.接触以后，物质点会渗入到刚体内，如图2所示.为反映真实的物理情况，渗入到刚体内的物质点需要进行位置再分配，将该物质点分配到临近处的刚体表面，如图2(c)所示. 图2 物质点的重新分配Fig.2 Relocation of material particle 新分配物质点k的坐标可以通过下式来计算： (7) (8) (9) (10) (11) 1.3 失效准则 当键处于拉伸状态，s为正值.由于s不依赖于变形方向，故这种材料是各向同性的.为在本构模型中引入材料破坏的概念，使用一个比较简单判定条件，即某一时刻键伸长率超过一定值时，认为物质点之间的键永久发生断裂.从此刻起，此对物质点之间力密度为0.该条件称为近场动力学的破坏准则. 设μ是历史变形判断标准，当键伸长率达到破坏准则，键发生断裂，μ=0；反之，μ=1，可定义为 (12) 式中：s0为键破坏时键伸长率的极限值，可称为极限伸长率.对于典型的标准