多段铰接式扬矿硬管系统建模与动力学分析.

1、多段较接式扬矿硬管系统建模与动力学分析徐博文（黑龙江科技大学，矿业工程学院，采矿 13-1 , 2013020658【24】） 摘 要：为了减弱深海采矿系统中扬矿管道螺纹连接处的应力集中与弯曲应力，并减小振动在管壁内的传递，采用球较接头对扬矿硬管进行连接，替代原有的螺纹连接。基于柔性多体系统动力学理论，建立多段较接式扬矿 硬管系统的运动学模型，分别对四、六级海况下扬矿系统的横向运动进行数值模拟。根据流体力学相似理论， 建立1?100扬矿系统实验系统并通过实验验证数学模型。研究结果表明：六级海况下较接式扬矿系统水下第1根管所受最大弯矩为 322 kN m,刚性连接系统的最大弯矩为1.890 MN

2、-m；较接式扬矿系统中，第 1根管的弯曲变形量最大，为4.79 mm；应用柔性多体系统动力学理论来分析较接式扬矿系统是可行的；减小系统弯矩， 较接式扬矿系统明显优于固接式扬矿系统；弯曲变形对系统运动影响很小。关键词：扬矿系统；较接；柔性多体系统动力学；横向运动；数值模拟Modeling and dynamics analysis of multiple articulated liftingsystemAbstract: To reduce stress concentration and bending stress on the screw joint, and decrease vibr

3、ation transmission through the pipes, spherical joint was used to connect lifting pipes to replace screwed connection in lifting pipes of deep-sea mining system. Based on dynamics of flexible multi-body systems, the mechanics model of articulated lifting system was established. Under the four-grade

4、and six-grade oceanic conditions, dynamic responses of lifting system weresimulated. Based on the fluid mechanics and similarity theory, experiment model of lifting system with the scale of 1 ?100was established to verify mechanics model. The results show that the maximum moment of flexion is 322 kN

5、 m on thefirst pipe under six-grade condition in the articulated system. The maximum moment of flexion is 1.890 MN m in thefixed one. The maximum bending deformation is 4.79 mm on the center of first pipe in articulated system. It is feasible toanalyze articulated lifting system by applying theory o

6、f flexible multi-body dynamics. Articulated lifting system is betterthan fixed one on reducing moment of flexion clearly. And deformation has a little effect on the motion of the system.Key words: lifting system; articulate; flexible multi-body dynamics; lateral movement; numerical simulation深海采矿系统中

7、，扬矿系统主要由供矿装置、柔性软管、中间舱、提升泵、主提升管道等组成1扬矿硬管系统由多个管道组成，为了便于装卸，相互间采用螺纹连接0。整个扬矿系统在海水中承受着由采矿船的纵摇、横摇、升沉以及洋流、海浪等的作用，呈现出复杂的动力学特性3。很多研究者对扬矿系统的运动特性进行了研究，如：P. Bar-Avi 4用解析法分析了等直径扬矿硬管在涡流作用下的横向振动；K. H.Kim等回采用模态叠加法将该系统把扬矿硬管考虑成连续的受拉伸的梁；凌胜 等6采用伽辽金（Galerkin）方法对带有集中质量的阶梯式扬矿硬管的横向偏移进行建模，发现刚性连 接的扬矿系统承受很大的弯曲应力，连接部位容易出现应力集中、应

8、力腐蚀，影响管道系统的强度。 为了减弱管道螺纹连接处的应力集中与弯曲应力，并减小振动在管壁内的传递，采用球较接头对扬矿硬管进行连接，替代原有的刚性连接，只需对系统的 2个关键位置进行分析，就可以说明扬 矿系统的运动规律。在扬矿系统中，泵与中间舱的形状不同于扬矿硬管，在仿真过程中作为关键位 置进行描述。在此，本文作者基于柔性多体系统动力学理论，建立多段较接式扬矿系统的数学模型， 研究分析扬矿系统的运动特性。1研究方法与动力学建模图1所示为球钱接头结构，采用上、下接头连接件和卡环相连接的形式7,螺纹部分采用API(AmericanPetroleum Institute)标准网，用螺纹脂密封，可有效

9、防止扬矿硬管中物质的外漏。将这种球钱接头应用在扬矿硬管之间，扬矿系统可看作“刚柔耦合管线”。应用拉格朗日方程10对钱接式扬矿系统进行动力学建模。扬矿硬管嵯圈蝶母螺检 甘环卜接头球较接头管端 外球传球公端连接件连接件图1球较接头结构Fig,l Spherical joint stmctiue of lifting systemi.i波流联合作用的液动力海流和波浪水质点的速度产生拖曳力，波浪水质点的加速度产生惯性力。根据实际海况求得海流速度和波浪水质点的速度和加速度，假设波浪为线性波，采用莫里森方程(Morison)计算流体作用于扬矿系统的载荷iio在我国勘定的开辟区，洋面流速为0.772 m/s

10、,海底流速为0.15 m/s；根据文献12给出海深y处的海流速度公式: TOC o 1-5 h z “八(1 000 + V V-、，广二0,15 + 0.622 x - 口(1)cI. 1 000 J根据莫里森方程，管体所受的波浪载荷为：耳=g C/(4口 + 口 + 鹿 |& +jlD2 o -、小下产仙生。44 T式中：CD为拖曳力系数，根据我国海港水文的规定选取CD=1.2; CM为惯性力系数，根据横截面形状，选取CM=2； Uwn和wn U .分别为水质点的法向速度和法向加速度，与海浪波高、波长 及周期有关。扬矿硬管的变形柔性多体系统含有柔性部件，因而其变形不可忽略。采用钱接接头连接

11、扬矿硬管，相邻两接头 之间的管单元可以视为简支梁。在海流海浪联合液动力作用下，管体截面上有弯矩也有剪力。扬矿 硬管单元长为11 m,外径为0.219 m,长度远大于外径，所以，剪力对弯曲变形的影响可以忽略。假设扬矿硬管中点作用有一正弦力Fsin(wt ),则挠曲线方程10为:,、2FI)sin(1 v(a t) =-EIti4.7LY I . 3兀丫 1 . 5 兀Y sinsin+ smL 34 L 54 L从式(3)可以看出，高阶模态的影响迅速衰减，分析系统弯曲变形时， 扬矿硬管的振型仅取第 1阶模态即可得出近似值。附加质量在流体中运动的物体，由于物体的动能和周围流体的动能都随时间而变化，

12、周围流体的附加质量在此时将对物体产生附加作用力13 o在深海采矿系统中，扬矿系统的运动会产生一个散射速度势，使周围的附加流体质量改变运动状态，即改变了速度场的分布，使扬矿系统受到1个附加水动力载荷作用。假设波面是理想、不可压缩的，波的运动是无旋的。当物体在不可压缩理想流体中运动时，其 附加质量只与物体的几何形状有关，而与流场的运动状态无关。扬矿硬管管体单位长度附加质量为 入=Pnr 2(其中，p为海水密度，r为管体半径)。数学建模相邻2段扬矿硬管之间的约束为球钱形式。每段管与中心线的夹角为e i(i=1,2,，n-1, n),如图2所示。整个系统在洋流、海浪及重力、浮力作用下(暂不考虑船的运动

13、对其在轴向的激励)，并且假设在z方向没有运动，则扬矿系统可看作是在二维空间内受相互约束的运动系统，绕O点在垂直平面内摆动。选取 。为广义坐标，在 x-y平面内建立拉格朗日动力学方程。图2扬矿系统的简化模型Fig*2 Simplified model of lifting system对于1 km深海采矿海试系统钱接式扬矿系统，其每段管长为1,质量为m ,较接接头质量为 M在300m和900 m处分别加装提升泵和中间舱。提升泵和中间舱在水中的质量分别为6.077 0 t和2.170 5 to为避免扬矿硬管在泵启动时发生旋转，泵与扬矿硬管之间采用螺纹连接；中间舱与扬矿硬管采用“十”字型连接装置。初

14、始条件为：e i(0)=0 , 0)0( = i 9 .；环境参数：4级海况，其中波高为2.5 m,周期为8 s； 6级海况，其中波高为 4 m,周期为10 s。I,则管的动能为:T 1 , 4 =5(加+i2 .(厩)2 + 彳(R)2 +打一17=1在重力场中，因为每根管都要绕质心在垂直平面内摆动，设转动惯量为FL11 Jl2Z Z 2/ 2 4% 85(4 % ) + COS (d -4)22r -ir COS h74Fit1C0S(n?)(7)2jt8Ein5fl a +；=i(4)式(4)中：右边最后3项为扬矿硬管变形对系统动能的贡献;Fsin wt为系统所受的液动力；w为海浪的圆频

15、率；EI为扬矿硬管材料的抗弯刚度。(EZ)2JTS接头的动能为:fl-1 n-1Z(M)2+ZZ2/2q%cos(4-%)0 (5)r=l系统的总势能可以表示为:fl 1 M4=-mglcos0j - nig -cosw -A/g/2cos 0j=i(6)根据扬矿系统载荷分析和初始条件，利用拉格朗日方程，建立微分方程组。cLcL=Qi 式中：i=1,2,，n； Qi为作用在第i个扬矿硬管上的作用外力之和；L为拉格朗日函数。方程组用矩阵的形式表示为:式中：A和C为系数矩阵；A为nxn阶方阵；C为nxi阶列阵。矩阵中的元素根据以上能量分析及初始条件来确定，是接头和扬矿硬管的质量及海水的附加质量、材

16、料的抗弯刚度等参数的组合表达式；n为硬管总数，每段扬矿硬管的长度为11 m,整个系统共需要 82段。在四、六级海况下，系统顶端的采矿船所受的激励分别为1 877.81sin (2nt/8)和3 810.17sin (2nt/10)。根据以上分析，运用 MATLAB 软件分别对四、六级海况下，泵及中间舱的横向运动进行数值模拟， 结果见图3和图4。从图3可知，在四级海况下，泵和中间舱的最大横向位移分别为8.0 m和12.6 m从图4可知，六级海况下，泵和中间舱的最大横向位移分别为9.7 m和13.5 m,曲线都呈周期变化。曲线在初期有一段滞后，这是因为运动传递到泵和中间舱需要一段时间。时间S1一泵

17、：2一中间舱图3四级海几下扬矿系统潢向位移计算维紧Fig-8 Computed results of latern 1 displacement of lifriiig system under fbvu-giade conditionif-1 2F 20406080100120时间总 1一泵；2 中间舱图4省级海况下扬矿系统横向位移计算结果Fig.4 Computed results of lateral displacement oflifting system tuidei six-grade condition2实验验证根据流体力学相似理论及量纲分析14,建立1?100扬矿系统实验模型

18、。通过正弦机构产生的运动模拟波浪，并带动扬矿硬管模型运动，进而模拟海洋环境下扬矿硬管的力学特性。实验系统主要包括：正弦机构（用来模拟采矿船运动）；调速系统；数字采集系统。一串较接连接 的扬矿硬管模型是实验的研究对象，由正弦机构带动。调速系统由三相异步电机和变频器组成，主要设备和仪器包电机型号 Y112M-4,功率为4 kW,额定转速为1 500 r/min ；减速器为双轴型一 级减速，传动比为 11,机型号为A型4级。变频器型号为： CIMR-F7A。数字采集系统由数字采集 卡及同步摄像头组成，数字采集卡型号为中佳T1500,摄像头型号为 DV-2410CDW。实验在装满水的有机玻璃圆筒内进行

19、。圆筒直径为1 m,总高度为10 mo运动周期通过调节变频器的输出频率来改变。运动位移通过数字采集卡与摄像头采集。摄像头和数字采集卡与计算机相 连，在计算机上对数据进行分析15 o从电动机传递到正弦机构的转速为：148050 x 11xH*2.7H其中：“1 480”为电机的同步转速；“11”为减速器传动比；H为变频器显示的数值。对钱接式扬矿系统在四、六级海况下的运动进行实验模拟。在四级海况下，扬矿系统的实验曲线见图5。六级海况下扬矿系统的实验曲线见图6。从图5可知，泵和中间舱的最大横向位移分别为7.2 m和11.1m。由图6可知，泵和中间舱的最大横向位移分别为9m和13 m。对比图3与图5,

20、图4与图6可知，数值模拟与实验结果间存在一定的偏差。其原因大致为：实验采用的介质为水，与实际海水有所差别；系统的数学模型进行了简化，并且计算机模拟时的计算步 长不能精确确定；处理数据时存在人为误差。但总体来说，计算结果与实验模拟趋势相同，大体相 符。说明所建模型是正确的。1612-12111J048121620时间人 1泵；2中间舱图5 四级海况下扬矿系统横向位移实脸结案Fig.5 Expeiiinental results of lateral displacement oflifting system luider fom-grade condition时间/s1一泵：2中间舱图6六级海况

21、下扬矿系统横向位移实验结果Fig,6 Experimental resiilTs of lateiTil displacemeiiT oflifting svsTem imder six-srade conditionU 痴J在六级海况下，水下第1根管承受整个系统的重力，并且液动力在此处最大，所以，其承受弯矩也是最大的。分析比较刚性连接系统与钱接连接系统在同一海况下运动时所受的弯矩见图7。从图7(a)可知：刚性连接的扬矿系统，在相同海况下，水下第 1根管的最大弯矩为 1.890 MN m,且弯 矩只在1.875 MN-m附近变化；而钱接连接的扬矿系统，水下第1根管的弯矩最大值为322 kN m

22、(见图7(b),对应的弯曲应力为312 MPa,这对于扬矿硬管的材料和结构，是较安全的。由此可知，钱(Exwsx接式连接可以减小系统所受弯矩，有利于海上实际作业。0)刚性连接杨犷 为味:(b) 接连接拘IT拿洗S 7 在左投海况t不闫扬行第电第1根首的直姆曲热 ig.7 %fccn匠cit of Hexioii on the first pipe in dierem liftinz syreni= nziiercxeaD：c c口口出口口口图8所示为扬矿硬管单元在某一时刻的弯曲变形量。变形最大的是第1根管，变形量为4.79 mm,其他各点变形量为0.51.5 mm ,变形量很小，说明弯曲变形对

23、系统运动影响不大。204060部 J钟讯L效图8 六就海况下揭#段管单元中点弯曲变舟曲线Fia.8 Bending deformation cime of pipes under six-zrade condition3结论a.基于柔性多体动力学理论，考虑海底运动中产生的附加质量，建立多段较接式扬矿硬管系统 的运动学模型，并对四、六级海况下扬矿系统的运动进行数值模拟。在四级海况下，泵和中间舱的 最大横向位移分别为 8.0 m和12.6 m；在六级海况下，泵和中间舱的最大横向位移分别为9.7 m和13.5 m。b.根据相似理论建立了实验装置，并对四、六级海况下扬矿系统的横向运动进行实验模拟。实

24、验结果为：在四级海况下，泵和中间舱的横向位移分别为7.2m和11.1 m；在六级海况下，泵和中间舱的最大横向位移分别为9 m和13 mo计算结果和实验结果变化趋势相同，除去误差因素，2种结果基本相符，说明应用柔性多体系统动力学理论分析钱接式扬矿系统是可行的。c.在六级海况下，扬矿系统的水下第 1根管所受的最大弯矩为322 kN m；对于刚性连接的扬矿系统，其第1根管的最大弯矩为 1.890 MN m,说明钱接式连接可以减小系统所受的弯矩；第 1根 管的弯曲变形量最大，为4.79 mm,其他各点变形量为0.51.5 mm,对系统的运动影响不大。此结论可为研究1 km较接式连接扬矿系统进行海上开采

25、（试马歙）提供理论依据。f-y 、. 4 卜参考文献：Chung J S, Whitney A K. Flow-induced moment and lift for acircular cylinder with cable attachmentJ. International Journal of Offshore and Polar Engineering, 1993, 3(4): 280287.Chung J S, Cheng B R. MSE and FEM modeling of thrusts toelastic joints of long vertical pipe in 3

26、-D nonlinear motionsJ.International Journal of Offshore and Polar Engineering, 1999,9(2): 117125.CHENG Bao-rong, HOU Zhi-chao, ZHENG Zhao-chang.Improved discrete element method for linear and nonlinear dynamic analysis of 3-D beam structuresJ. Tsinghua Scienceand Technology, 1997, 2(2): 52535.Bar-Av

27、i P. Dynamic responses of risers conveying fluidJ.Journal of Offshore Mechanics and Arctic Engineering,2000(122): 188-193.Kim K H, Choi H S, Seong W J, et al. Dynamic positioningcontrol for the buffer of a deep-sea mining system using FIRfilterC/Proceeding of the Eleventh International Offshore andP

28、olar Engineering Conference. Stavanger, Norway, 2001: 422.6凌 胜，肖林京，申会华，等.深海采矿开采系统运动状态和动态特性影响因素分析研究J.中国工程科学，2002, 4(3):78-83.LING Sheng, XIAO Lin-jing, SHEN Yan-hua, et al. A study of the factors influencing the kinematic condition and dynamic characteristics of deep seabed mining systemsJ. Engineeri

29、ng Science, 2002, 4(3):-783.FENG Ya-li, ZHANG Wen-ming, FENG Fu-zhang. Design and strength analysis of a spherical connector for lifting subsystem in deep sea mining systemJ. China Ocean Engineering, 2006,20(4): 60573.8王涮，张汝忻，邹家祥，等.API圆螺纹套管接头应力场分布试验J.北京科技大学学报，2002, 24(12): 555- 558.WANG Li, ZHANG Ru

30、-xin, ZOU Jia-xiang, et al. Experimental research of stress field distribution on API round thread casing connectionJ. Journal of University of Science and TechnologyBeijing, 2002, 24(12): 555558.9冯雅丽，冯福璋，张文明，等.扬矿球钱接头的螺纹滑脱及脂密封性能分析J.润滑与密封，2006(12): 12&131.FENG Ya-li, FENG Fu-zhang, ZHANG Wen-ming, et

31、 al. Thread off and thread compound performance analysis for the spheric joint of deep sea bed mining lifting systemJ. Lubrication Engineering, 2006(12): 128-131.10陆佑方.柔性多体系统动力学M.北京：高等教育出版社,1993.LU You-fang. Dynamics of flexible multibody systemsM.Beijing: Higher Education Press, 1993.11竺艳蓉.海洋工程波浪力学M.天津：天津大学出版社,1991.ZHU Yan-rong. Ocean engineering wave mechanicsM. Tianjin: Tianjin University Press, 1991.Aso K, Kan K, Doki H. Experim