• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 21
    • 下载费用:30 金币  

    重庆时时彩后二组选复式: 直升机典型旋转部件的空间多体运动仿真分析方法.pdf

    关 键 词:
    直升机 典型 旋转 部件 空间 运动 仿真 分析 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201510822574.1

    申请日:

    2015.11.24

    公开号:

    CN106777438A

    公开日:

    2017.05.31

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/50申请日:20151124|||公开
    IPC分类号: G06F17/50 主分类号: G06F17/50
    申请人: 中国直升机设计研究所
    发明人: 张树桢; 孙东红; 陈浩; 李明强; 赖凌云; 付双检; 王国胜; 宋芹芹; 马小艳; 范学伟; 沈安澜; 王文涛
    地址: 333001 江西省景德镇市航空路6号
    优先权:
    专利代理机构: 中国航空专利中心 11008 代理人: 郭平
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201510822574.1

    授权公告号:

    |||

    法律状态公告日:

    2017.06.23|||2017.05.31

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明属于结构动力学分析范畴,特别是针对直升机典型旋转部件的空间运动分析。通过多体动力学的建模理论,对直升机典型旋转部件进行建模分析,主要包括旋翼和减速器系统。其动力学方程可应用于分析各类直升机旋转结构的空间运动问题,可在直升机旋翼系统设计的各阶段进行动力学特性的相关研究,得到可靠的运动轨迹和作用力结果。其需要的各参数既可以通过有限元计算方法给出,也可以通过实际测量得到。该计算分析可用于初步计算,给出初步的旋转部件空间运动情况及桨根处的作用力大小,由于数值积分方法的引入,尤其适用于求解单片或多片桨叶失效,以及其他类的故障条件下的运动分析,直接为设计方法、故障类型的判断及故障条件下的运动情况提供指导性的理论依据。

    权利要求书

    1.直升机典型旋转部件的多体运动仿真分析方法,其特征是,包括以下步
    骤:
    第一步,建立仿真分析模型如下:
    原点O为桨叶根部传动轴的轴心处,建立公共基坐标系o-x-y,设桨
    叶质心处为原点oi,建立桨叶质心坐标系oi-xi-yi,设传动轴与桨叶连接
    点为球铰o′,建立坐标系o′-x′-y′,第i片桨叶相对其球铰原点坐标的转角
    为φi;由于桨叶相对传动轴对称布置,且实际条件下的各片桨叶参数基本一致,
    设定桨叶的质量为m,桨叶长度为l,按照材料属性均匀分布的假设,得到转动
    惯量为J=1/3ml2;设定传动轴处传递至球铰处的稳定转动角速度为ω,得
    到球铰o′的约束方程为

    其雅可比矩阵表达式为

    独立约束方程个数为3,则拉格朗日乘子σi1,σi2,σi3即为3个,即
    <mrow> <msub> <mi>&sigma;</mi> <mi>i</mi> </msub> <mo>=</mo> <msup> <mfenced open = '(' close = ')'> <mtable> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>3</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
    其重力加速度的坐标阵为
    <mrow> <mi>g</mi> <mo>=</mo> <msup> <mfenced open = '(' close = ')'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>g</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    则其作用力项为
    <mrow> <msub> <mi>z</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mi>m</mi> <msup> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mi>g</mi> </mrow> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
    设铰点处的刚度项为k,则其作用力为
    <mrow> <msub> <mi>z</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <msup> <mfenced open = '(' close = ')'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>k&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
    桨叶的动力学方程的简略表达式为
    <mrow> <mi>M</mi> <mover> <mi>X</mi> <mo>&CenterDot;&CenterDot;</mo> </mover> <mo>+</mo> <msub> <mi>&Phi;</mi> <mrow> <mi>q</mi> <mi>i</mi> </mrow> </msub> <msub> <mi>&sigma;</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>F</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
    其具体表达式为
    <mrow> <mtable> <mtr> <mtd> <mrow> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <mi>m</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>m</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> <msup> <mi>ml</mi> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <msub> <mover> <mi>x</mi> <mo>&CenterDot;&CenterDot;</mo> </mover> <mi>l</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>y</mi> <mo>&CenterDot;&CenterDot;</mo> </mover> <mi>l</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mover> <mi>&phi;</mi> <mo>&CenterDot;&CenterDot;</mo> </mover> <mi>l</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>lsin&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>lcos&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>3</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mfenced open = '[' close = ']'> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>m</mi> <mi>g</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>k&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
    由约束方程可得如下的加速度约束方程

    将式(9)代入式(8),得到其拉格朗日乘子,从而计算得到拉格朗日乘子的表
    达式为
    <mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mn>3</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>l</mi> <mi>m</mi> <msubsup> <mover> <mi>&phi;</mi> <mo>&CenterDot;</mo> </mover> <mi>i</mi> <mn>2</mn> </msubsup> <msub> <mi>cos&phi;</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>lsin&phi;</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>i</mi> </msub> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mn>3</mn> <mo>(</mo> <mrow> <msup> <mi>l</mi> <mn>2</mn> </msup> <msup> <mi>m&omega;</mi> <mn>2</mn> </msup> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&omega;</mi> <mi>t</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>l</mi> <mi>m</mi> <msubsup> <mover> <mi>&phi;</mi> <mo>&CenterDot;</mo> </mover> <mi>i</mi> <mn>2</mn> </msubsup> <msub> <mi>sin&phi;</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>m</mi> <mi>g</mi> <mo>+</mo> <msub> <mi>lcos&phi;</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>i</mi> </msub> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>/</mo> <mn>3</mn> <mo>(</mo> <mrow> <msup> <mi>l</mi> <mn>2</mn> </msup> <msup> <mi>m&omega;</mi> <mn>2</mn> </msup> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&omega;</mi> <mi>t</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>1</mn> <mo>/</mo> <mn>3</mn> <mrow> <mo>(</mo> <mrow> <msup> <mi>l</mi> <mn>2</mn> </msup> <msup> <mi>m&omega;</mi> <mn>2</mn> </msup> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&omega;</mi> <mi>t</mi> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&phi;</mi> <mi>i</mi> </msub> <mi>k</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
    由于笛卡尔坐标系,桨叶坐标系相对铰点坐标系的方向余弦矩阵为
    <mrow> <msub> <mi>A</mi> <mi>i</mi> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>cos&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>sin&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>sin&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>cos&phi;</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
    且将式(2)进行相应处理,得到其相对雅克比表达式为

    由多体动力学的相关理论得到作用力的表达式为
    <mrow> <msub> <mi>F</mi> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>&Phi;</mi> <mo>&OverBar;</mo> </mover> <mrow> <mi>q</mi> <mi>i</mi> </mrow> </msub> <msub> <mi>A</mi> <mi>i</mi> </msub> <msubsup> <mi>C</mi> <mi>i</mi> <mi>p</mi> </msubsup> <mo>)</mo> </mrow> <mi>T</mi> </msup> <msub> <mi>&sigma;</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
    <mrow> <msub> <mi>M</mi> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>&Phi;</mi> <mo>&OverBar;</mo> </mover> <mrow> <mi>q</mi> <mi>i</mi> </mrow> </msub> <msub> <mi>IA</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mover> <mi>&Phi;</mi> <mo>&OverBar;</mo> </mover> <mrow> <mi>&phi;</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mi>T</mi> </msup> <msub> <mi>&sigma;</mi> <mrow> <mi>i</mi> <mi>r</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
    其中I为单位阵。同时,通过力Fi,Mi和刚度k的关系f=kx,得到各桨
    叶的位移和相对球铰的转角。
    由于传动轴的输入角速度相对稳定,同理得到各个桨叶和传动轴的相对转角
    和作用力。
    对于空间旋翼模型,在平面模型的基础上增加z向作用力和相对转角θ,
    得到球铰o′的约束方程,其方程式的推导如平面模型一致,并代入减速器的质
    量mz和转动惯量Jzx,Jzy,Jzz,同样根据计算或实测得到球铰处的刚度k,
    具体为绕铰点坐标的3个转动刚度kθx,kθy,kθz;
    第二步,通过测量和有限元计算得到所需的各参数数值,包括单片桨叶的质
    量m,单片桨叶绕重心的转动惯量J,传动轴中心到铰点的距离l,桨叶质量中
    心相对传动轴坐标系的坐标(xi,yi),桨叶间的夹角由减速器传递至传动
    轴的角速度ω,减速器的质量mz,绕自身质心坐标系的坐标Jzx,Jzy,Jzz,
    等效计算出球铰连接位置处的扭转刚度kθx,kθy,kθz及相应的阻尼系数;
    第三步,将第二步测量获得桨叶的动力学参数输入第一步的模型中,得到式
    (7)中各矩阵的表达式;
    第四步,对得到的方程进行常微分方程组的计算求解;可采用龙格库塔变步
    长积分方法,求解铰点处的位移和速度结果。并在此基础上,对铰点处的作用力
    进行求解,得到多体动力学分析中较为关心的逆问题求解,即已知运动结果的条
    件下,求解系统的输入力,从而得到桨叶和减速器的运动情况。

    说明书

    直升机典型旋转部件的空间多体运动仿真分析方法

    技术领域

    本发明属于结构动力学分析范畴,特别是针对直升机典型旋转部件的空间运
    动分析。

    背景技术

    直升机的空间姿态保持完全依靠旋翼所提供的作用力。目前,机体的动力学
    分析基本以有限元方法为主,通过结构离散化对系统进行处理。旋翼由多片桨叶
    构成,单片桨叶的静刚度、疲劳及振动等问题均可以一定程度以有限元方法进行
    相关分析,其结果可信度较高。

    随研究问题全面性以及实际情况复杂性的提高,有限元方法在实际应用中难
    以对典型旋转部件的运动情况进行研究,主要由于旋转部件的转速高、质量大、
    部件相对位移显著,因此必须采用多体动力学的思路对旋转部件进行建模,以分
    析其空间运动情况,满足对旋转部件动平衡、局部失效的运动情况以及进一步分
    析结构的正逆多体动力学问题。

    发明内容

    本发明的目的是:

    利用多体动力学的仿真思路,建立旋转部件的空间运动方程,以此为基础,
    分析桨叶和减速器整体的空间运动情况,解决如桨叶故障等工况条件下的直升机
    旋转部件动力学分析问题。

    本发明的技术方案是:

    直升机典型旋转部件的多体运动仿真分析方法,包括以下步骤:

    第一步,建立仿真分析模型如下:

    设桨叶根部传动轴的轴心为原点O,建立公共基坐标系o-x-y,

    设桨叶质心处为原点oi,建立桨叶质心坐标系oi-xi-yi,

    设传动轴与桨叶连接点为球铰o′,建立坐标系o′-x′-y′,在球铰坐标系
    中,第i片桨叶相对其球铰原点坐标的转角为φi。

    设定每片桨叶的质量为m,桨叶长度为l,按照材料属性均匀分布的假设,
    得到转动惯量为J=1/3ml2,设定传动轴输出的转动角速度为ω,得到球铰
    o′处的约束方程Φ为


    则雅可比矩阵表达式为


    其独立约束方程个数为3,则拉格朗日乘子σi1,σi2,σi3即为3个,即

    σi=(σi1 σi2 σi3)T (3)

    其重力加速度的坐标阵为

    g=(0 -g 0)T (4)

    则其作用力项为

    zi1=(0 -mg 0)T (5)

    设球铰处的刚度项为k,则其作用力为

    zi2=(0 0 -kφi)T (6)

    桨叶的动力学方程的简略表达式为


    式(7)的具体表达式为


    其中为球铰处的x向、y向和转角的加速度量。由约束方程可
    得如下的加速度约束方程


    将式(9)代入式(8),得到其拉格朗日乘子,从而计算得到拉格朗日乘子的表
    达式为


    由于笛卡尔坐标系,桨叶坐标系相对铰点坐标系的方向余弦矩阵为


    且整理式(2),得到其相对雅克比表达式和为


    由多体动力学的相关理论得到作用力的表达式为



    其中I为单位阵。同时,通过力Fi,Mi和刚度k的关系f=kx,得到各桨
    叶的位移和相对球铰的转角。

    由于传动轴的输入角速度相对稳定,同理得到各个桨叶和传动轴的相对转角
    和作用力。

    对于空间旋翼模型,在平面模型的基础上增加z向作用力和相对转角θ,
    得到球铰o′的约束方程,其方程式的推导如平面模型一致,并代入减速器的质
    量mz和转动惯量Jzx,Jzy,Jzz,同样根据计算或实测得到球铰处的刚度k,
    具体为绕铰点坐标的3个转动刚度kθx,kθy,kθz。

    第二步,通过测量和有限元计算得到所需的各参数数值,包括单片桨叶的质
    量m,单片桨叶绕重心的转动惯量J,传动轴中心到铰点的距离l,桨叶质量中
    心相对传动轴坐标系的坐标(xi,yi),桨叶间的夹角由减速器传递至传动
    轴的角速度ω,减速器的质量mz,绕自身质心坐标系的坐标Jzx,Jzy,Jzz,
    等效计算出球铰连接位置处的扭转刚度kθx,kθy,kθz及相应的阻尼系数。

    第三步,将第二步测量获得桨叶的动力学参数输入第一步的模型中,得到式
    (7)中各矩阵的表达式。

    第四步,对得到的方程进行常微分方程组的计算求解,得到桨叶相对传动轴
    的转角。对该常微分方程组的求解,可采用龙格库塔变步长积分方法,求解铰点
    处的位移和速度结果。并在此基础上,对铰点处的作用力进行求解,得到多体动
    力学分析中较为关心的逆问题求解,即已知运动结果的条件下,求解系统的输入
    力。

    本发明的有益效果

    直升机典型旋转部件的动力学方程可应用于求解各类直升机旋转结构的多
    体动力学分析??稍谥鄙砩杓频母鹘锥谓卸ρ匦缘南喙胤治?,得到
    可靠的运动轨迹和作用力结果。其需要的各参数既可以通过有限元计算方法给
    出,也可以通过实际测量得到。该计算分析可用于初步计算,给出初步的旋转部
    件空间运动情况及桨根处的作用力大小,由于数值积分方法的引入,尤其适用于
    求解单片或多片桨叶失效,以及其他类的故障条件下的运动分析,直接为设计方
    法、故障类型的判断及故障条件下的运动情况给出指导性的理论依据。

    附图说明

    图1为本发明所述平面单片桨叶模型示意图,其中,1为传动轴及位置,2
    为球铰及位置,3为桨叶及其质心位置。

    图2为本发明所述平面旋翼模型示意图。

    图3为本发明所述空间旋翼模型示意图,其中,1为本说明书中,编号为1
    号的桨叶,2为本说明书中,编号为2号的桨叶,3为本说明书中,编
    号为3号的桨叶,4为本说明书中,编号为4号的桨叶,5为减速器。

    图4为本发明所述摆振弯矩的时间历程。

    图5为本发明所述摆振角度的时间历程。

    图6为本发明所述挥舞弯矩及角度的时间历程。

    图7为本发明所述扭转弯矩及角度的时间历程。

    图8为本发明所述脱离时刻为0.5s,尾减质心的运动轨迹。

    图9为本发明所述脱离时刻为0.5s,尾减质心的运动轨迹投影。

    图10为本发明所述脱离时刻为0.6s,尾减质心的运动轨迹。

    图11为本发明所述脱离时刻为0.6s,尾减质心的运动轨迹投影。

    具体实施方式

    为表明该仿真分析方法的实际应用价值,对单片桨叶缺失条件下的旋转部件
    系统进行仿真分析,且采用时间延续的计算方法,模拟该系统在某时刻自机身脱
    离后的空间运动轨迹和桨叶的摆振、挥舞和扭转状态,并在Matlab软件中作图,
    以对该特殊条件下的运动情况进行分析。

    直升机典型旋转部件的多体运动仿真分析方法,包括以下步骤:

    第一步,建立仿真分析模型如下:

    桨叶、传动轴和减速器系统的空间构型如图3所示,设定1号,2号为一组,
    3号,4号为一组,1、2号和3、4号之间沿传动轴方向的间距为0.2m,3,4号
    桨叶靠近传动轴顶端,且1号和4号桨叶的角度为60°,设定相关仿真参数与约
    束条件:

    1)桨叶无气动载荷且完全刚性;

    2)外侧3号桨叶缺失;

    3)旋转部件系统在未脱离机体前,空间固定,即无地面相对速度;

    4)在桨叶根部和传动轴之间附加万向铰接和等效刚度,模拟桨根的连接形
    式;

    5)在某一时刻去除减速器根部约束及其输出的驱动角速度,以模拟由于离
    心力非对称,导致减速器部件自机身脱离的运动状态,分为0.5s和0.6s
    时刻两种情况。

    计算减速器正常约束至脱离机身条件下,1、2和4号桨叶的摆振、挥舞和
    扭转时间历程,并设定脱离时刻为0.5s和0.6s以作对比分析。

    根据图1所示,原点O为桨叶根部传动轴的轴心处,建立公共基坐标系
    o-x-y,设桨叶质心处为原点oi,建立桨叶质心坐标系oi-xi-yi,设
    传动轴与桨叶连接点为球铰o′,建立坐标系o′-x′-y′,第i片桨叶相对其球
    铰原点坐标的转角为φi。由于桨叶相对传动轴对称布置,且实际条件下的各片
    桨叶参数基本一致,设定桨叶的质量为m,桨叶长度为l,按照材料属性均匀分
    布的假设,得到转动惯量。设定传动轴处传递至球铰处的稳定转动角速度为ω,
    得到球铰o′的约束方程为


    其雅可比矩阵表达式为


    独立约束方程个数为3,拉格朗日乘子即为3个,即

    σi=(σi1 σi2 σi3)T (3)

    其重力加速度的坐标阵为

    g=(0 -g 0)T (4)

    则其作用力项为

    zi1=m(0 -g 0)T (5)

    设铰点处的刚度项为k,则其作用力为

    zi2=(0 0 -kφi)T (6)

    桨叶的动力学方程的简略表达式为


    其具体表达式为


    由约束方程可得如下的加速度约束方程


    将式(9)代入式(8),得到其拉格朗日乘子,从而计算得到拉格朗日乘子的表
    达式为


    由于笛卡尔坐标系,桨叶坐标系相对铰点坐标系的方向余弦矩阵为


    且将式(2)进行相应处理,得到其相对雅克比表达式为


    由多体动力学的相关理论得到作用力的表达式为



    同时得到各桨叶的位移及相对球铰的转角。对于空间旋翼模型,在平面模型
    的基础上增加z向作用力和相对转角θ,得到球铰o′的约束方程,其方程式
    的推导如平面模型一致。

    第二步,根据实测和计算得到的桨叶的质量为3kg,桨叶绕自身重心的转动
    惯量为0.4kg·m2,传动轴中心到铰点的距离0.3m,,桨叶间的夹角60°传动
    轴的角速度为125.7rad/s,减速器的质量为40kg,绕自身重心的转动惯量为
    0.79kg·m2,0.79kg·m2,0.67kg·m2,桨叶根部的等效摆振刚度为
    2.86×104 N·m/deg,挥舞刚度为600N·m/deg,扭转刚度为150N·m/deg。

    第三步,将第二步测量获得桨叶的动力学参数输入第一步的模型中,得到式
    (7)中各矩阵的表达式。

    第四步,对得到的方程进行常微分方程组的计算求解。采用龙格库塔变步长
    积分方法,首先求解系统自静止状态下,当减速器固定,传动轴输入定频转速后,
    桨叶的运动情况;其次当系统稳定后,令减速器脱离,对桨叶-传动轴-减速器组
    成空间运动系统进行积分运算,求解桨叶的摆振、挥舞、扭转和减速器质心的空
    间运动情况;最后将得到的解进行Matlab出图,分析其摆振、挥舞和扭转的运
    动情况,以对该类问题进行总结,并依据质心的运动情况,分析其运动模式。

    由于微分方程组的复杂性,其理论解析较为困难,同时由于多体动力学分析
    问题的复杂程度较高,需设定不同的初始条件,以对系统进行相应的动力学分析,
    采用数值积分方法不仅有利于求解定常的动力学运动情况,也有利于求解非定常
    条件下的运动学分析。

    图4,5中,从零位移,零速度的初始条件进行仿真分析。桨叶自静止,受
    到传动轴的驱动,从而发生摆振,由于阻尼的存在,该瞬态运动历经0.05s消除。
    脱离条件下,内侧及外侧尾桨叶的摆振角度稳定在±0.3°以内。

    图6中,在减速器固定条件下,由于未施加气动载荷,因此无挥舞发生。在
    脱离条件下,外侧编号1的挥舞弯矩及角度在0.508s之前大于内侧,随后发生
    交替。

    图7中,在减速器固定条件下,由于未施加气动载荷,所以扭转弯矩及扭转
    角不存在。脱离条件下,扭转弯矩和扭转角并不一一对应,同其他运动存在复杂
    的耦合。

    将减速器的质心处设置为参考点,对质心进行运动轨迹分析??占湓硕旒?br />分别如图8、10所示,在航向、侧向及垂向的投影如图9、11所示。其中图8、
    9描绘了0.5s时刻减速器自机身脱离后的运动状态,图10、11描绘了0.6s时
    刻减速器自机身脱离后的运动状态。

    在单片桨叶缺失条件下,重点分析减速器自机身脱离后,剩余桨叶及减速器
    系统的运动特性。通过多体动力学仿真分析表明:

    1)桨叶的挥舞运动十分显著。在减速器脱离的短时间内,桨叶即产生了较
    大的挥舞角,扭转其次,摆振运动较不明显;

    2)桨叶和减速器系统的运动轨迹,受脱离瞬时桨叶离心力的合力指向影响
    突出。在不同的脱离时刻,桨叶和减速器系统的运动轨迹存在显著的差异,该差
    异也表明了运动的复杂性和仿真的真实性;

    3)桨叶的摆振、挥舞和扭转运动的计算结果具有普遍适应性。不同脱离时
    刻,尾桨叶相对传动轴连接处的摆振、挥舞及扭转运动只存在和不同脱离时刻相
    关的时滞,表明计算结论具有普遍性。

    关于本文
    本文标题:直升机典型旋转部件的空间多体运动仿真分析方法.pdf
    链接地址://www.4mum.com.cn/p-6021044.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 重庆时时杀一码 哪里买云南时时走势图 中巨奖的征兆 北京pk10下载 骰子猜大小赌博软件 买平肖出特肖算不算百度知道 老时时20110701 百乐门怎么玩 功夫时时计划软件免费版 金沙手机投注 北京pk10公式大全论坛 福彩3d猜大小怎么玩 大小单双稳赚倍投法 11选5任5稳赚投注技巧 3d组六6码遗漏 赛车怎么买稳赚不赔