大变形复合材料柔性梁刚度计算与试验验证
2017-2-20 来源：中航工业直升机设计研究所 作者：熊绍海 金坤健 陈雪冰
摘要：根据复合材料柔性梁的结构与变形特点，基于Halpin-Tsai等公式对材料模量进行修正，采用桨叶专用软件Blade computer aided desigh(BCAD)对柔性梁的刚度(除扭转刚度)进行计算，扭转刚度采用经验公式进行计算。结果表明该方法不仅能提高柔性梁剖面刚度计算的准确性，而且也为无轴承旋翼／尾桨的动力学、栽荷计算的准确性提供了保证。
关键词：复合材料；无轴承；柔性梁；刚度计算
无轴承旋翼／尾桨通过柔性梁的弹性变形实现旋翼的挥舞、摆振及变距运动[1]，柔性梁还传递旋翼几乎所有的载荷。通过刚度的合理裁剪，可以使柔性梁满足弹性耦合运动协调、强度和动力学等多目标要求，因此，柔性梁设计是无轴承旋翼的核心技术。
复合材料柔性梁的剖面特性是无轴承旋翼设计的重要元素，而柔性梁结构特性分析的基础是梁理论，最经典的是Euler—Bernoulli梁理论，该理论假定变形前的横剖面在变形后仍保持为平面Timoshenko理论在此基础上引入了剖面上的剪切变形影响，对于变形后的横剖面并不要求与弹性轴垂直口3；Bauchau[31利用特征函数对Euler—Bernoulli梁的圣维南解进行修正，提出了一种更加符合实际的梁理论；Giavotto等[41建立了适用于平面内和面外翘曲的各向异性梁理论；为了分析大变形复合材料梁的力学特性，Hodges等引入截面三维翘曲及任意旋转张量概念，采用混合变分法，写出了大变形复合材梁的本构方程[5]，并在此基础上开发了工程应用软件VABS：6|。中国国内向锦武等[7]综合考虑直升机旋翼桨叶梁变形特点以及梁剖面内外翘曲对剖面刚度的影响，编写了计算各向异性材料桨叶剖面刚度的有限元软件。
柔性梁一般都采用单向纤维大梁带和编织布混杂铺层的方式，这种大变形梁的刚度不仅受复合材料的各向异性影响，而且与柔性梁的外形、复合材料的铺层角度也有很大关系，因此，应用上述梁理论和各软件计算得到的柔性梁刚度误差往往超过10％，从而影响到无轴承旋翼的动力学和载荷计算结果。
为了进一步提高柔性梁刚度计算精度，本文基于Halpin—Tsai等公式对材料模量进行修正，采用桨叶专用软件BCAD对柔性梁的刚度进行计算，并与刚度测试结果进行对比。
1.BCAD建模方法和计算流程
BCAD(Blade computer aided design)软件专门用于计算直升机桨叶的结构特性，BCAD设计的基本思想是将三维弹性问题分解为二维线性剖面特性分析和一维非线性分析，其计算流程详见图1，经过A，B，C三个模块的程序的填充，可以定义出整片桨叶(或柔性梁)外形，然后通过模块D对各剖面的铺层进行定义，此定义将引用P模块中的材料库数据(材料中的各项参数可以修改)，通过D模块的计算得到除了扭转刚度以外的其他结构参数，扭转刚度采用G模块另外计算。
图1 BCAD专用软件计算流程
式中：Ei为剖面各铺层的弹性模量
为各铺层积分面积；z为距离弹心的厚度；Y为距离弹心的宽度。
为各铺层积分面积；z为距离弹心的厚度；Y为距离弹心的宽度。 2.柔性梁结构特点与材料模量修正
2．1柔性梁结构特点
柔性梁一般分为5个功能区[8]，如图2所示，根部连接区、根部过渡区、扭转变形区(典型段)、桨叶过渡区和桨叶连接区。柔性梁的挥舞变形一般都布置在根部过渡区，扭转变形和摆振变形布置在扭转变形区。其中过渡区厚度和弦长都会逐渐变化，纤维的铺层角度与展向呈一定夹角，对于这种大变形梁的刚度计算，不能直接用材料的弹性模量E，而要对其模量进行一定的修正。
图2柔性梁典型结构
2．2材料模量修正
柔性梁在结构计算时，其主承力材料单向纤维大梁带的弹性模量E和G和常规桨叶的大梁带模量有区别。以某柔性梁所用的单向纤维大梁带为例，根据外形和铺层角度，对大梁带的模量进行修正(见表1，2)。表1中树脂参数和纤维参数已知，大梁带计算数据引用公式如下
表1 某单向纤维大梁带特性数据
表2某单向纤维大梁带在BCAD中弹性模量的设置
图3 大梁带纤维方向与纵向夹角示意
图4柔性梁剖面纤维铺层角度示意
3.柔性梁刚度试验原理
试验需测试的项目为：挥舞刚度、摆振刚度及扭转刚度‘9。。试验采用应变片粘贴方式，在柔性梁各被测剖面处粘贴挥舞、摆振的应变片。
3．1挥舞弯曲刚度测试
柔性梁中央区域固支，在柔性梁加载位置(0．8R，R为尾桨半径)处，使用专用翼型夹具将其夹紧，通过夹具在柔性梁轴线位置上施加垂向力。柔性梁弦平面水平放置，对于不同的测量剖面，调整旋转工作盘以保证被测剖面的弦线处于水平状态。
在加载点位置，利用砝码块加载，载荷级数为5级。试验时，将两组平行应变片组成全桥，测量各级载荷下的输出微应变值。
3．2摆振弯曲刚度测试
柔性梁中央区域固支，在柔性梁加载位置(o．8R)处，使用专用翼型夹具将其夹紧，通过夹具弦向加载。柔性梁弦平面垂直固定，前缘向上、后缘向下，调整旋转工作盘以保证被测剖面弦平面与水平面垂直。
在加载点位置，利用砝码块加载，载荷级数为5级，试验时，选取4个平行应变片(其中2片受拉，2片受压)组成全桥进行测量。对多次应变测量结果取其平均值，将应变值除以2，按式(14)计算柔性梁各级载荷下的摆振弯曲刚度，然后求刚度平均值，得到对应剖面的摆振刚度值。
3．3扭转弯曲刚度测试
柔性梁扭转刚度测试时，柔性梁桨毂连接端弦平面水平固定，保证扭矩作用剖面处于水平状态。在柔性梁加载位置(o．8R)处施加扭矩，在被测剖面左右20 mm处各布置两个激光位移测量传感器(传感器分布在柔性梁中心线的两侧)，用于测量各级力偶矩作用下的位移值。
4.柔性梁刚度计算与测试结果对比
4．1挥舞和摆振刚度对比
本文对柔性梁刚度进行2种状态的计算，将修正前、后的挥舞刚度计算值与试验值[11]进行对比(见表3)，及修正前、后的摆振刚度计算值与试验值[1妇进行对比(见表4)。
(1)大梁带模量每个剖面都选用E=100 000 MPa。
(2)大梁带模量修正，即按表2所列不同剖面不同模量来计算刚度。
表3修正前、后挥舞刚度计算值和测试值对比
表4修正前、后摆振刚度计算值和测试值对比
从表4数据可以得到：
(1)摆振和挥舞刚度测试平均值与修正前计算值的偏差基本都是负偏差，除个别剖面外，偏差都在一10％以内。
(2)摆振和挥舞刚度测试平均值与修正后计算值的偏差基本都是正偏差，且偏差都在4％以内。模量修正后的计算值与实测值更为吻合。
4．2扭转刚度对比
BCAD软件自带的扭转刚度计算主要适用于薄壁盒腔结构，如常规桨叶的铺层结构。而全叠层铺层结构的复合材料截面扭转刚度计算比较繁琐，计算结果和实际测量的结果符合性不好，因此，柔性梁的扭转刚度采用如下经验公式计算得到[12]
表5 BCAD和经验公式扭转刚度计算值和测试值对比
5.柔性梁刚度对尾桨动力学和载荷的影响
5．1 柔性梁刚度对尾桨动力学的影响
5．1．1 柔性梁刚度对尾桨动特性影响
为了判断柔性梁结构参数对尾桨动特性和气弹稳定性的影响，本文计算了5种不同柔性梁刚度的尾桨动特性和气弹稳定性，这5种结构参数是柔性梁理论刚度以及在此基础上分别将柔性梁刚度增加10％，20％及减少10％，20％计算结果，如表6所示。
由表6可见，柔性梁刚度变动对低阶动力学特性有较大影响，假如柔性梁刚度计算值偏小10％或20％，则摆振1阶频率在大总距下与2n接近甚至重合(图5)，极大地影响了柔性梁的结构设计结果。
表6尾桨动特性结果
5．1．2 柔性梁刚度对气弹稳定性影响
针对摆振刚硬无轴承尾桨，通过气弹稳定性计算，分析尾桨在某使用总距下或某使用尾桨转速下，是否会出现气弹耦合不稳定现象。由于1阶摆
图5 1阶摆振频率随总距变化曲线
振模态的气动阻尼最小，因此1阶摆振最容易出现不稳定现象，所以在气弹稳定性分析中最为关注此模态。
本文计算了某尾桨1阶摆振阻尼随柔性梁的刚度变化趋势(图6)，得知随着摆振刚度的增加1阶摆振阻尼有明显的增加，也就是说增加无轴承尾桨的摆振刚度可以提高尾桨的气弹稳定性。
图6 1阶摆振阻尼随柔性梁刚度变化
5．2柔性梁刚度对尾桨载荷影响
目前，多个软件的计算结果表明[1 3|，桨叶挥舞弯矩的计算结果较为准确，而摆振弯矩和扭转力矩的计算结果与实测结果相差较大，原因在于桨叶的挥舞刚度较低，挥舞固有频率远远低于气动载荷中的高频成分，气弹耦合程度较低。制约摆振弯矩和扭转力矩载荷计算精度的因素是扭转响应的求解。扭转响应直接影响桨叶的有效攻角，而有效攻角决定了桨叶截面的气动力。因此，桨叶和柔性梁的扭转刚度是载荷计算扭转动力学方程中的关键因素。另外，柔性梁扭转刚度计算的准确性还影响到拉杆载荷和尾桨助力器作动筒行程的设计。因此，柔性梁刚度计算的准确性很大程度决定了尾桨载荷的准确性。
6.结论
(1)基于Halpin—Tsai等公式对材料模量进行修正，采用桨叶专用软件BCAD对柔性梁的刚度(除扭转刚度)进行计算，扭转刚度采用经验公式进行计算，计算结果与实测结果吻合度很好。该材料模量修正方法同样适用同类桨叶结构参数计算软件。
(2)柔性梁刚度计算的准确性将对无轴承旋翼／尾桨的动特性和稳定性有较大影响，特别是低阶频率和低阶阻尼：
投稿箱：
如果您有机床行业、企业相关新闻稿件发表，或进行资讯合作，欢迎联系本网编辑部， 邮箱：skjcsc@vip.sina.com
如果您有机床行业、企业相关新闻稿件发表，或进行资讯合作，欢迎联系本网编辑部， 邮箱：skjcsc@vip.sina.com
更多相关信息
- target=_blank> 舰载发射智能靶弹在试验海域飞行控制中的关键技术
- target=_blank> 电动多旋翼无人机螺旋桨的性能计算与分析
- target=_blank> 四轮驱动电动车驱动防滑模糊控制研究
- target=_blank>分布式驱动电动客车驱动防滑控制效果分析
- target=_blank>船用螺旋桨振动特性数值计算与分析
新闻资讯
| 更多
- target=_blank>MB2120B型数控内圆磨床电主轴变频器的改造
- target=_blank>【雷尼绍】RESOLUTE™与MELSERVO-J5强强组合，实现更快速、精确、可靠的运动控制
- target=_blank>NSK开发低摩擦轮毂单元轴承 可增加电动汽车续航里程
- target=_blank>T68卧式镗床进给系统的数字化改造
- target=_blank>AGV+协作机器人在零件数控机床加工上下料中的应用
- target=_blank>2020年7月高端装备制造业、工业机器人，行业运行简述
- target=_blank>浅谈线切割机床中走丝与慢走丝
- target=_blank>ANCA整体PCD铣刀 —— 提高生产效率带来新的发展机遇
- target=_blank>伊斯卡，不止专注于金属加工
- target=_blank>EMAG成功收购Scherer Feinbau（舍勒公司），极大扩展了公司产品范围以及客户群体
