基于ansys船舶轴系的振动校核计算

《基于ansys船舶轴系的振动校核计算》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、基于ANSYS船舶轴系的振动校核计算1概述　　船舶轴系是由推力轴、中间轴、艉轴、推力轴承、滑动轴承、联轴节、螺旋桨等组成的复杂系统，在船舶运行过程中，它会发生弯曲振动现象，对船舶正常运行产生不利影响。船舶轴系振动有三种类型：由旋转轴不平衡引起的横向振动，可以是垂直方向的，也可以是水平方向的，会造成艉管密封漏水或漏油，轴承座松动，甚至破裂;由螺旋桨推力不均匀引起的纵向振动，情况严重时可以造成推力轴承敲击，曲柄箱破裂，有齿轮传动时，还会损坏齿轮;此外，从主机通过轴系传递功率至螺旋桨造成轴段来回摆动，各轴段间的扭角不相同，从而

2、产生扭转振动，破坏的结果是轴系断裂，有齿轮传动时，会造成齿轮敲击。因此，在船舶设计过程中，有必要对船舶轴系进行振动校合计算。　　对于轴系这样的复杂结构，运用有限元方法进行振动计算具有明显的优越性。本文针对上海交通大学和某造船厂共同设计开发的46000吨集装箱船，应用ANSYS有限元软件6.0版本对其传动轴系进行振动校合计算，为进一步的设计提供参考。ANSYS是美国ANSYS公司开发的大型通用有限元分析软件，它具有结构静力分析、结构动力分析、瞬态分析、模态分析、流体动力学分析、电磁场分析等多种功能。本文即是利用ANSYS软

3、件的模态分析功能，完成对船舶轴系这一复杂结构的建模和有限元分析。实践证明，这种方法可以有效的提高工作效率，缩短分析周期，对工程实际是非常有效的。　　2轴系计算的有限元模型　　进行校合计算的46000吨集装箱船，采用的是瓦西兰公司的32缸柴油发动机组，发动机输出法兰通过齿轮箱变速后，和中间轴连接，中间轴和艉轴之间有联轴节。中间轴长3.68m，外径0.4m，无轴承支承。艉轴长5.3m，外径0.48m，前后分别有两个轴承，前轴承宽0.48m，后轴承宽1.08m，轴承刚度由轴承说明书给出。中间轴和艉轴中都布置有润滑系统。螺旋桨是

4、变距螺旋桨，总重14500kg。根据实际需要，只需对船舶轴系的自由振动情况进行校合计算，不考虑受迫振动情况。所以在轴系的有限元建模中，只保留从齿轮箱输出法兰到螺旋桨部分的轴系。　　根据轴系的实际结构，建模过程中进行了以下简化：　　对轴系本体部分采用BEAM188梁单元模拟。BEAM188单元是三维梁单元，每个节点具有六个自由度：UX、UY、UZ、ROTX、ROTY、ROTZ，可以满足各种振动计算的要求。设置不同的梁截面，可以模拟不同直径的轴结构。考虑到润滑系统的布置，这里都设置为内径100外径不同的环形截面。　　对弹性支

5、承的轴承部分采用COMBINE14弹簧单元模拟。COMBINE14通常是一维线性弹簧单元，可以分别有三个方向的自由度UX、UY、UZ，只沿弹簧方向传递力。由于轴承有一定的宽度，可以有力矩作用，所以考虑在轴承部分的每个节点上都设置弹簧单元，来模拟力矩对轴承的影响。由于是一维弹簧单元，所以考虑在轴的水平和垂直方向分别设置两个弹簧，来分别模拟轴承部分在Y向和Z向的弹性。所以最后是在轴承部分的每个节点上有两个弹簧单元，弹簧单元一端直接连接在轴的节点上，一端设置为固定端。　　在轴系和齿轮箱法兰的连接处，考虑存在弹性连接，所以在纵向

6、上设置一个弹簧单元来模拟纵向的弹性连接，弹簧的刚度由经验数据给出。在水平和垂直方向上也设置两个弹簧，来模拟齿轮箱法兰对轴系的支承作用。　　对联轴节部分，为了计算方便将其同样简化为梁单元，梁单元的内径不变，只是将梁单元的外径适当放大，来模拟这部分的强度。　　对螺旋桨部分，将艉轴部分适当延长来模拟螺旋桨部分的长度，将螺旋桨的质量加上附水质量(变距桨按30%的螺旋桨干质量计算)简化为集中质量，集中质量直接加在螺旋桨的几何中心位置。　　经过以上简化处理，可以建立轴系的有限元计算模型，见图1。轴系共有节点63个，其中方向节点27个

7、，BEAM188梁单元27个，采用了5种不同的截面形状，COMBINE14弹簧单元15个，MASS21质量单元1个。　　材料的弹性常数为：弹性模量E=2.1x1011N/m2，泊松比μ=0.3，密度ρ=7.8x103g/m3。　　　　图1船舶轴系的有限元计算模型　　3轴系横向振动的计算　　轴系横向弯曲振动计算中，假设轴承的刚度在各个方向上是相同的，轴系在水平和垂直方向上的振动是相同的，所以只计算垂直方向的振动。ANSYS模态分析中，BEAM188单元只保留UY、ROTZ自由度，其他自由度都去掉。模态分析后可以得出各阶固有

8、频率，各节点的相对位移值、转角值，各单元的弯矩值、剪力值。如果在模态分析的结果上，作垂直方向上的谐响应分析，就可以得到各阶模态对应的模态参与因子。　　横向振动的固有频率见表1。第一阶固有频率14.286Hz下的参数值见表2。前两阶的计算结果图示如下，见图2-9。　　表1横向振动的固有频率　　　　表2横振频率f=12.