山地轨道运载机牵引钢丝绳空间结构分析与几何建模

欧阳玉平，王天玉，孙 晗， 孙旭东,2，洪添胜，黄志平，姜小刚,2

(1. 华东交通大学 机电与车辆工程学院，南昌 330013；2. 载运工具与装备教育部重点实验室(华东交通大学)，南昌 330013；3. 华南农业大学 工程学院，广州 510642；4. 广东振声科技股份有限公司，广东 梅州 514700)

随着我国水果产业的壮大，实现果园运输机械化不仅仅是农业机械化的一个重要组成部分，同时也是水果产业不断发展的重要标志. 传统的果园种植缺少科学合理规划，立地条件差，果树大多生长在难以形成较完善交通运输网，甚至是陡峭的梯田上. 这就导致了常规的水果采运车难以在该环境下推广使用，直接给成熟果子的运输造成了较大的困难[1]. 目前，我国90%的柑橘种植在山地，土地条件差，山地果园机械化水平低，生产主要靠人力完成，劳动强度非常大. 特别是山地果园肥料和果品运送方面主要靠肩挑或背篓背负，必须用机械化代替艰苦的劳动. 柑橘产业要持续发展，农机农艺相融合是唯一出路. 未来强壮的“果农”就是农机农艺融合[2]. 山地果园钢丝绳牵引式双轨运输机是一种由钢丝绳牵引，载物滑车沿着具有一定坡度的轨道行驶的运载工具，轨道由2条相互平行的圆管及辅助横梁焊接而成，轨道坡度集中在10°～40°之间[3]. 国内外已研制出多种山地果园轨道运输机械，山地果园牵引式双轨运输机是其中一种代表机型[4]. 该运输机缓解了山地地区果农的劳动强度，有效实现了果品及其他农资省力化运载需求[5]. 由于运输机作业环境复杂，运载果品时若频繁振荡，易造成钢丝绳振荡冲击损伤，因此有必要对该类运输机的牵引系统的失效行为及钢丝绳受力情况进行研究[6].

钢丝绳是一种金属结构，在许多工业工程、海洋工程和土木工程领域已经使用了很长一段时间[7-8]，在吊桥、起重机、建筑电梯和矿井提升机等领域也有其应用[9-10]. 钢丝绳经常被使用，因为它们具有较高的抗拉强度和弯曲灵活性[11]，而且与其他结构相比，它们具有更高的强度和更长的寿命[12]. 在理论推导、实验研究和有限元模拟这三种途径中，几何模型是研究钢丝绳力学模型及失效行为的研究基础. 通过精确的数学表达式，建立钢丝绳的精确几何模型，从而进行相应的力学分析[13]. 王等[14]提出了钢丝绳和股线中心线的数学表示及其半径的计算方法, 给出了钢丝绳和股线缠绕关系的空间参数几何方程，并用经验值或迭代法以外的方法计算钢丝和钢绞线的半径是合理的. I.I. Argatov 等[15]建立了基于Archard-Kragelsky磨损定律的微动磨损数学模型，并给出了磨损定律参数校准的示例. W. Ma等[16]首先推导出每根导线的质心轴(单、双或三螺旋)的坐标方程，然后利用MatlabTM生成形成该轴的控制点坐标的数据文件，最后根据得到的基螺旋的坐标方程，用递推法建立了高水平螺旋. 吴等[17]利用Frenet框架和微分几何建立了钢丝绳二次螺旋的数学模型, 并在此基础上推导出二次螺旋线的几何参数，利用Love弹性细杆理论建立钢丝绳的等效力学模型，推导出钢丝绳的等效弹性模量和等效剪切模量的计算公式. Ivan Argatov[18]建立了简单螺旋钢丝绳股的精细离散数学模型，通过泊松比和局部接触变形(钢丝压扁)研究了钢绞线横向收缩的影响.

以上研究都对不同环境下的钢丝绳给出了精准的空间几何表达式，目的都是为了更好地研究钢丝绳力学性能及其失效行为，但针对山地轨道运载机牵引钢丝绳的空间几何表达式目前还没有明确的定论. 因此，本文为了后续更深入研究山地轨道运载机牵引钢丝绳力学模型的建立与损伤机理的研究，对适用于山地轨道运载机牵引钢丝绳的空间几何结构进行数学公式推导，同时推导了曲率及挠率方程及其连续变化规律，这样有利于减少在不同参数和操作条件下进行昂贵的试验测试的需要.

1 运输机钢丝绳选型及捻制特征

1.1 运输机钢丝绳选型

山地轨道牵引式双轨运输机属于新式农用运输机，钢丝绳应用特点与矿山或港口起重机等有较大差别，需对钢丝绳绳径进行重新计算，避免绳径选择不适，造成浪费或强度不够. 运输机钢丝绳作业时，钢丝绳在卷筒、约束轮及滑轮等约束下沿弯曲轨道转弯运行，应用工况与农用索道类似，选用圆股线接触右交互捻(SZ)，型号为6×19+FC类钢丝绳[19]. 在前面研究中建立钢丝绳牵引系统虚拟样机模型[4]，简化示意图如图1所示，图1中A点为载物滑车最低滑行位置，B点为载物滑车最高滑行位置.

1—卷筒; 2—载物滑车; 3—托辊; 4—横梁; 5—轨道; 6—载物滑车; 7—钢丝绳; 8—滑轮

钢丝绳最大有效长度与整绳最小破断拉力应满足式(1)、式(2)[20]，当载物滑车运行至A点时，钢丝绳与卷筒相接处O点承受最大的静拉力如式(3)所示，钢丝绳质量如式(4)[20]所示，整绳最大静拉力如式(5)所示：

L=L1+L2+L3,

(1)

(2)

Fmax-Mg(sinα+f1cosα)+mgL(sinα+f2cosα)，

(3)

m=KgD2,

(4)

Fmax=Mg(sinα+f1cosα)+KgD2gL(sinα+f2cosα).

(5)

式中：Fmax为整绳最大静拉力，kN；Fmin为钢丝绳最小破断拉力，kN；Ka为钢丝绳许用安全系数；L1为A、B两点间距离，m；L2为B点至顶端滑轮距离，m；L3为卷筒至顶端滑轮距离，m；M为载物滑车及所载物资总质量，kg；g为重力加速度，m/s2；α为轨道倾角，°；f1为轮轨间滚动摩擦系数；f2为钢丝绳与托辊间摩擦系数；m为钢丝绳质量，kg/m；D为钢丝绳公称直径，mm；Kg为钢丝绳重量系数，kg/100 m·mm2.

钢丝绳最小破断拉力Fn、钢丝绳破断拉力总和Fmin[21]为

Fn=KminD2Et,

(6)

Fmin=KhFn,

(7)

式中:Et为钢丝绳公称抗拉强度，MPa；Kmin为钢丝绳最小破断拉力系数；Kh为破断拉力换算系数.

将式(6)代入式(7)得

Fmin=KhKminD2Et,

(8)

将式(5)及式(8)代入式(2)得如下关系式：

Mg(sinα+f1cosα)+KgD2gL(sinα+f2cosα)≤

(9)

(10)

对于6×19+FC钢丝绳，f1取0.05，f2取0.1[19]；重量系数Kg取0.36，最小破断拉力系数Kmin取0.338，破断拉力换算系数Kh取1.214，g取9.8 m/s2；Ka取7，α取40°，M取500 kg，L取200 m；Et取1 670 MPa[20]；于是

6.03 mm.

查阅机械设计手册[19]，选取钢丝绳的直径为7.7 mm.

1.2 钢丝绳捻制特征

6×19+FC类圆股右交互捻钢丝绳由侧丝绕中心丝左旋捻制成单股，再由6侧股绕绳芯右旋捻制而成，侧股侧丝由内层丝及外层丝组成. 所采用绳芯材质为纤维，因此该类钢丝绳只有螺旋股，无直股. 钢丝绳捻制示意图如图2所示，每股由1条中心丝，6条内层丝与12条外层丝捻制而成. 为便于阐述钢丝间空间位置关系，对股及丝进行编号，以初始截面最左端为第1股，逆时针依次编为第1、2、3、4、5、6股；侧股中心丝逆时针依次编为1、2、3、4、5、6丝；侧股内层丝以每股中心线与绳芯中心线连线逆时针的第一钢丝中心为第01丝，逆时针依次编为第02、03、04、05、06丝；外层丝以每股中心线与绳芯中心线连线逆时针的第一钢丝中心为第07丝，逆时针依次编为第08、09、10、11、12、13、14、15、16、17、18丝.

(a)钢丝绳整体示意图

(b)钢丝绳截面示意图

将整绳解剖后依次展开，如图3所示，β0、β1、β2分别表示股在绳中、内层丝在股中、外层丝在股中的捻角；r0、r1、r2分别表示股在绳中螺旋半径(捻制半径)、内层丝在股中螺旋半径、外层丝在股中螺旋半径；S、S0、S1、S2分别表示整绳、侧股(中心丝)、侧股内层丝、侧股外层丝长度；αa、αb、αc分别表示钢丝绳在捻制过程中，股在绳中转过的角度、内层丝在股中转过的角度、外层丝在股中转过的角度.

图3 钢丝绳展开示意图

由图3几何位置可知：

AB=r0αa，BC=r1αb，BD=r2αc，

(11)

由钢丝绳已知数据和几何关系可得

r1=0.50 mm，r2=0.97 mm，r0=2.63 mm，

β0=18.30°，β1=7.16°，β2=13.10°.

结合钢丝绳已求结构参数，型号为6×19+FC钢丝绳的关键几何参数，如表1所示.

表1 运输机钢丝绳几何参数

2 钢丝绳螺旋方程

2.1 直线段钢丝绳空间螺旋方程

直线段钢丝绳六侧股一次螺旋线的初始捻制角依次记为θ1、θ2、θ3、θ4、θ5及θ6，以逆时针为正，初始捻制角分别为0°、60°、120°、180°、240°、300°；侧股中6条内层丝二次螺旋线的初始捻制角θ01～θ06分别为45°、105°、165°、225°、285°、345°；侧股中12条外层丝二次螺旋线的初始捻制角θ07～θ18分别为0°、30°、60°、90°、120°、150°、180°、210°、240°、270°、300°、330°.

本文以第1股及第1股第18丝为研究对象，建立直线段钢丝绳空间曲线模型，如图4所示. 图4中侧股中心丝中心线上点A用向量OA表示，侧丝中心线上点B用OB表示，点A处的切向量、法向量、从法向量依次记为t、n、b. 设向量AB为Frenet标架n-b-t内的平面矢量，则侧丝空间二次螺旋线可由向量OB=OA+AB表示.

令向量OA在x-y-z标架上为OAx，向量AB在n-b-t标架上为ABn，由图4钢丝的位置关系及钢丝绳已知参数可得

(12)

(13)

图4 钢丝绳第1股中心丝及其第18丝空间曲线模型

侧股缠绕于绳芯转换矩阵[n,b,t]表达式为

(14)

为求向量OB在坐标系x-y-z的表达式，ABn需转换为x-y-z标架表达式ABx. 由空间坐标变换原理可知，向量AB在坐标系x-y-z中表达式可表示为

ABx=[n,b,t]×ABn.

(15)

将式(13)及(14)代入式(15)可得

(16)

由式(12)及(16)可知，向量OB在坐标系x-y-z中表达式为

(17)

令

(18)

由式(12)可得，直线段钢丝绳第i股中心丝方程表达式：

(19)

将式(18)带入(17)可得直线段钢丝绳第i股内层第j丝空间方程表达式与直线段钢丝绳第i股外层丝第k丝方程表达式：

(20)

(21)

2.2 弯曲段钢丝绳空间螺旋线方程

运输机钢丝绳弯曲段体现在钢丝绳与约束轮、滑轮及卷筒接触，如图5所示. 为研究钢丝绳弯曲段空间螺旋方程，以钢丝绳绕过滑轮为例进行说明. 绳与轮接触时，绳中心轴线呈弯曲圆弧状态，绳中各股股芯中心线绕滑轮圆弧呈空间一次螺旋线，侧股侧丝中心线绕圆弧呈二次螺旋线，如图6所示.

(a)钢丝绳与约束轮接触 (b)钢丝绳与滑轮接触(c)钢丝绳与卷筒接触

1—钢丝绳中心路径; 2—钢丝绳单股中心路径; 3—钢丝绳单股侧丝中心路径

为构建弯曲段钢丝绳空间螺旋方程，建立绳轮接触空间曲线模型，如图7所示. 以动点C替代直线段钢丝绳坐标系原点O，建立弯曲段钢丝绳空间直角坐标系LMN. 图7中一次螺旋线绕弧线旋转，二次螺旋线绕一次螺旋线旋转，弯曲状态侧丝中心线的二次螺旋线向量O′B的表达式为

O′B=O′C+CA+AB.

(22)

1—侧股螺旋半径(整绳螺旋半径)对应圆; 2—侧股侧丝中心线；3—侧股侧丝螺旋半径画出的圆; 4—侧股股芯中心线

图7中，将x-y-z标架上的一次螺旋线与二次螺旋线分别视为建立在新标架L-M-N弯曲一次螺旋线与二次螺旋线，则向量O′C在标架L-M-N的表达式O′Cl为

(23)

标架x-y-z向标架L-M-N转换，对坐标系LMN的转换矩阵为

(24)

标架n-b-t对坐标系LMN的转换矩阵[N,B,T]为

[N,B,T]=[X,Y,Z]×[n,b,t]=

(25)

ABn转换为坐标系LMN上的表达式ABl为

ABl=[N,B,T]×ABn=

(26)

由图7可知，弯曲的一次螺旋线O′Al的向量方程表达式为

O′Al=O′Cl+CAl，

(27)

(28)

CAx转换为坐标系LMN上的表达式CAl为

CAl=[X,Y,Z]×CAx=

在广州，这样的优秀教师比比皆是，他们的故事书之不尽。华阳小学的曹定金老师为了放飞学生的梦想，独创了“第30号课室”；文德路小学的陈秀茹老师把一个不幸的孩子送的罐子作为自己的能量罐，让爱充盈课堂；开发区二小的余雪云老师为了纠正学生的不良习惯，给学生的心灵“解套”，连续4年对学生做思想工作，用不伤害孩子尊严的方法纠正了孩子的坏毛病……正是这些好教师，把教育梦想写成教育故事，在教育现场书写教育传奇，支撑起了广州好学校、广州好教育的高楼大厦。

(29)

O′Al=O′Cl+CAl=

(30)

O′Bl=O′Cl+CAl+ABl=

(31)

图8为钢丝绳侧股芯丝与滑轮接触示意图，式(31)中含有3个未知量，为求解式(31)关于αa的表达式，需求出图8中各变量的相互关系.

图8 钢丝绳侧股芯丝与滑轮接触示意图

在弯曲段钢丝绳中，αd与αa关系为

(32)

令

(33)

将式(18)及式(33)代入式(30)可得弯曲段钢丝绳第i股中心丝方程表达式为

(34)

将式(18)及式(33)代入式(31)可得弯曲段钢丝绳第i股内层第j丝方程为

(35)

将式(18)及式(33)代入式(31)可得弯曲段钢丝绳第i股外层第k丝方程为

(36)

3 钢丝绳曲率和挠率求解

钢丝绳空间螺旋方程、曲率和挠率变化规律是研究整绳三维实体建模及力学模型的基础，曲率是表示曲线在某点的弯曲程度，挠率是表示空间曲线的扭曲程度，钢丝绳空间螺旋曲线上任意点均有确定的曲率和挠率，整绳曲率及挠率具有一定变化规律.

钢丝绳曲率k、挠率τ与空间螺旋曲线方程f(αa)的关系如下[22]：

(37)

(38)

式中:f(αa)为钢丝绳空间螺旋曲线方程，x、y、z为对应的极角αa的函数. 由式(37)及式(38)，可分别求出直线段、弯曲段钢丝绳各钢丝的曲率与挠率空间表达式，但由于钢丝绳空间螺旋曲线方程较为复杂，求三次导后更为复杂，表示较困难. 应用MATLAB软件语言符号功能对钢丝绳的曲率与挠率进行求解，在求解弯曲段钢丝的曲率及挠率时，设定滑轮直径为180 mm. 钢丝绳曲率及挠率在MATLAB中的计算流程如图9所示.

将钢丝绳直立及弯曲段的二维曲线函数图像导入Photoshop进行合并，得到直立与弯曲段钢丝绳曲率与挠率关系图，如图10及11所示. 在图10中，直线段钢丝绳的曲率与挠率呈明显的周期性变化规律，变化周期均为180°；挠率峰值均出现于曲率的波谷处，曲率处于波峰时，挠率变化较为平缓；曲率与挠率的极值较明显. 在图11中，弯曲段钢丝绳曲率与挠率同样具有明显的周期性变化规律，变化周期均为360°；曲率处于波峰时，挠率发生突变；曲率处于波谷时，挠率变化率最大；曲率和挠率的极值大小明显.

图9 钢丝绳曲率及挠率在MATLAB中计算流程

图10 直线段钢丝绳外层丝曲率与挠率变化规律

图11 弯曲段钢丝绳外层丝曲率与挠率变化规律

4 钢丝绳三维建模

上文推导了钢丝绳的空间数学方程、曲率和挠率表达式，分析了钢丝绳在直立及弯曲状态，侧股中侧丝的曲率及挠率的连续变化规律. 为更直观了解该类钢丝绳的结构，通过已建立钢丝绳空间数学模型，运用MATLAB参数化功能及Solidworks软件曲线、扫描和实体阵列等功能建立钢丝绳三维实体模型.

应用SolidWorks曲线及扫描功能建立钢丝绳侧股芯丝，通过Matlab获取钢丝绳侧股侧丝的二次螺旋线参考点坐标，制成txt文件，如表2为第1股第18丝参考点系列坐标；然后运用Solidworks“XYZ曲线”功能，在参数输入设置中导入已制好的坐标文件；最后应用SolidWorks扫描、装配及实体阵列建立钢丝绳整绳三维几何模型. 图12～14分别为单丝、丝与丝缠绕、直立段及绳轮处钢丝绳三维几何模型.

(a)侧股芯丝 (b)侧股芯丝与内层丝缠绕 (c)侧股芯丝、内层丝与外层丝缠绕

图13 直立状态钢丝绳整绳三维实体模型

图14 弯曲段钢丝绳三维实体模型

经过参数设计及三维建模最后得出第1股第18丝参考点系列坐标表格如表2所示。

表2 第1股第18丝参考点系列坐标

5 结 论

1)在确定运输机钢丝绳型号及直径基础上，对钢丝绳股、丝进行编号，运用微分几何学及空间坐标变换原理，推导了6×19+FC类钢丝绳在直立及弯曲状态的空间螺旋方程，并推导了该类钢丝绳侧股外层钢丝曲率及挠率的连续变化规律，为钢丝绳在各种载荷下的工作状态的数值分析和实验分析提供理论依据.

2)图10结果表明：直线段钢丝绳的曲率与挠率呈明显的周期性变化规律，变化周期均为180°；挠率峰值均出现于曲率的波谷处，曲率处于波峰时，挠率变化较为平缓；曲率与挠率的极值较明显. 图11结果表明：弯曲段钢丝绳曲率与挠率同样具有明显的周期性变化规律，变化周期均为360°；曲率处于波峰时，挠率发生突变；曲率处于波谷时，挠率变化率最大；曲率和挠率的极值大小明显.

3)根据所推导的钢丝绳空间螺旋方程、曲率及挠率表达式，运用MATLAB参数化及SolidWorks曲线、扫描及装配等功能，绘制了运输机钢丝绳直立及弯曲状态的三维几何模型，并得到运输机钢丝绳股、丝参考点系列坐标，为后续建立钢丝绳力学模型及分析结构参数对钢丝绳损伤的影响奠定基础.