基于ANSYSLS—DYNA不同初速度下船—船碰撞分析.docVIP

基于ANSYSLS—DYNA不同初速度下船—船碰撞分析 　　摘要：在船舶航行过程中，船舶碰撞而产生的海难事故较为多发，因此研究船舶碰撞是十分必要的。本文基于ANSYS/LS-DYNA14.5应用软件，采用显示有限元求解方法，以两艘295TEU集装箱船相撞为研究对象。主要介绍了船体结构建模，并分析了船-船碰撞速度，能量转换和损伤情况。讨论不同初速度下，被撞击船舷侧的损伤，从而为日后的船舶结构设计、评估和维护提供理论支持和实践参考。 　　关键词：集装箱船；碰撞；非线性有限元；数值模拟 　　Analysis of Ship-Ship Collision at Seriesof speedbased on ANSYS/LS-DYNA 　　Fu Yuliang，Jiang Guohe，Zhao Kaiqi 　　（Merchant Marine College，Shanghai University，Shanghai 201306，China） 　　ABSTRACT：During the voyage of the ship，more collision shipwreck generated，it is necessary to study the ship collision.Therefore，according to explicit finite element method，based on ANSYS/LS-DYNA software，the collision processes is two 295 TEU container ships.It is mainly introduced the structure modeling，obtain and analysis of the speed of ship-ship collision，energy conversion and damage.Discuss the influence of different initial velocity to the struck ship.Thus，it can provide theoretic support and practical reference to ship structure design，evaluation and maintenance. 　　KEYWORDS：Container ship；Collision；Nonlinear finite element；Numerical simulation 　　引言 　　船舶碰撞是海上事故中比较严重的情形，可能导致船舶的结构变形、损伤、货物泄露等，严重时会危及到生命和财产安全。因此，无论是从安全上、经济上，海上从环境保护上来看，研究并提高船舶的抗碰撞性能都具有重要的意义[1]。本文选取了295TEU集装箱船相互碰撞为例，基于ANSYS/LS-DYNA软件进行分析。该295TEU集装箱船总长120.38m，两柱间长113.00m，型宽19.0m，型深8.0m，方形系数CB为0.730，采用横骨架和纵骨架组合式结构，双壳双底，双层底高1300mm。 　　1.船舶碰撞数值仿真的基本理论 　　1.1显示有限元求解方法 　　在关于动力学分析过程中，通常有两种方法可以计算，一种是隐式算法（Newmark法），另一种是显式算法（中心差分法）。其中，显式算法适合于求解瞬间动态问题[2]。在有限元方法中，碰撞的运动方程可以???示为： 　　[M]{α}+[C]{v}+[K]{d}={Fex} 　　式中，[M]为质量矩阵；{α}为加速度向量；[C]为阻尼矩阵；{v}为速度向量；[K]为刚度矩阵；{d}为位移向量；{Fex}为包括碰撞力在内的外力矢量。在ANSYS/LS-DYNA动力分析过程中采用中心差分法时，无需求解联立的方程组，也不必对矩阵分解或者求逆。在求解的过程中存在一个重要的问题就是选取合适的时间步长。由于此法的稳定性，为保证其精度，选择的时间步长不宜超过稳定的临界时间步长。在实际问题中，选取临界步长的时候用最小有限元网格特征长度和应力波速之比，来作为近似结果[3]，即为： 　　Δt≤Δter=min（Lε/c） 　　1.2沙漏控制 　　单个积分点的实体单元和壳单元容易形成零能模式，主要表现为产生一种自然震荡并且比所有结构响应的同期要短很多，网格变形呈现锯齿状，易产生沙漏变形。必须有效控制分析中可能出现的沙漏变形，以确保分析的正确性。当发生沙漏现象时，容易在计算过程中造成沙漏能量损失，这会对系统的总能量产生一定的影响。当显式动态分析使用缩减积分单元时，沙漏能量不能超过内能的10%即被认为是可以接受的分析结果。 　　2.船舶碰撞有