CN109726358A - 一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法 - Google Patents
一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法 Download PDFInfo
- Publication number
- CN109726358A CN109726358A CN201711035311.1A CN201711035311A CN109726358A CN 109726358 A CN109726358 A CN 109726358A CN 201711035311 A CN201711035311 A CN 201711035311A CN 109726358 A CN109726358 A CN 109726358A
- Authority
- CN
- China
- Prior art keywords
- jet pipe
- actuator
- coordinate system
- coordinate
- axis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 13
- 230000009466 transformation Effects 0.000 title claims abstract description 12
- 230000008878 coupling Effects 0.000 title description 5
- 238000010168 coupling process Methods 0.000 title description 5
- 238000005859 coupling reaction Methods 0.000 title description 5
- 230000008859 change Effects 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims abstract description 8
- 238000006243 chemical reaction Methods 0.000 claims abstract description 7
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000004044 response Effects 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 4
- 239000007921 spray Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000567 combustion gas Substances 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明属于适量耦合计算技术领域,具体涉及一种实现推力矢量控制用伺服系统的高精度与快速响应设计的基于三维坐标变换的前摆心喷管铰链耦合解耦算法;包括两台伺服伺服作动器均成90°配置,定义静坐标系O‑XYZ以喷管摆心为原点,固连于喷管摆心;X轴、Y轴分别为喷管2个单摆方向;动坐标系在初始状态下与静坐标系O‑X″Y″Z″一致,与喷管固连,随着喷管姿态的变化而变化,从喷管尾部看,静坐标系和动坐标系Z(Z″)轴垂直纸面向外;零位时,两个坐标系原点重合,坐标轴方向相平行;本发明包括以下步骤:步骤1.根据姿态角得到坐标系转换矩阵;步骤2.计算上述姿态角状态下作动器的位置;步骤3.得到作动器活塞杆的伸缩量。
Description
技术领域
本发明属于适量耦合计算技术领域,具体涉及一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法。
背景技术
推力矢量控制技术是通过改变火箭发动机或其他动力装置产生的燃气流方向来控制箭体的飞行方向和姿态角。固体火箭发动机推力矢量控制对伺服系统的控制精度、响应速度以及可靠性都有极高的要求,同时火箭控制系统对产品自身的重量、体积都有严格的要求。因此,在体积与重量等约束条件严格的条件下,实现伺服系统的精确控制和快速响应控制是一个关键的技术问题。
全轴摆动喷管实现摆动的驱动装置是2路互成90°夹角作动器的伺服机构,通常当一个方向伺服作动器运动时,另一个方向作动器由于存在机械限制,使喷管除了主要向摆动作动器伸缩方向运动之外,也将在另一个作动器方向产生移动,这就是铰链耦合效应。伺服机构的控制方案分为理想状态和非理想状态,理想状态是2路作动器不发生牵连,每路作动器伸缩时正负摆角对称,否则为非理想状态,对于导弹常采用的后摆心结构形式的喷管,作动器的牵连作用影响作用比较小,一般情况下可忽略。而对于前摆心深潜入式的喷管,由于喷管尺寸、结构等方面的限制,在推力向量控制过程中不能忽略2个方向的上相互牵连,为提高控制精度,需要考虑牵连效应对喷管摆动的影响,把非设计状态控制在一定范围内。
本发明提出一种推力矢量铰链耦合解耦控制用方法,对俯仰和偏航2个方向上作动器控制指令进行解耦计算。
发明内容
本发明的目的是,针对现有技术不足,提供一种实现推力矢量控制用伺服系统的高精度与快速响应设计的基于三维坐标变换的前摆心喷管铰链耦合解耦算法。
本发明的技术方案是:
一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法,包括两台伺服伺服作动器均成90°配置,定义静坐标系O-XYZ以喷管摆心为原点,固连于喷管摆心;X轴、Y轴分别为喷管2个单摆方向;动坐标系在初始状态下与静坐标系O-X″Y″Z″一致,与喷管固连,随着喷管姿态的变化而变化,从喷管尾部看,静坐标系和动坐标系Z(Z″)轴垂直纸面向外;零位时,两个坐标系原点重合,坐标轴方向相平行;
本发明包括以下步骤:
步骤1.根据姿态角得到坐标系转换矩阵;
步骤2.计算上述姿态角状态下作动器的位置;
步骤3.得到作动器活塞杆的伸缩量。
所述步骤一中
若姿态角为α,β为喷管在2个方向上的单摆角度,坐标系转换矩阵为:
所述步骤二中作喷管摆动时,作动器下支点会随着喷管运动,作动器上支点始终保持不动;作动器1、2的上支点坐标分别为(c1,c2,c3)、(d1,d2,d3)
作动器1、2下支点在喷管运动前,相对于动坐标为(a1,a2,a3),(b1,b2,b3),定义下支点计算坐标为:
作动器下支点在喷管运动后,相对于静坐标系使用向量表示为:
则作动器1、2方向分别摆动角度α,β之后两作动器的长度为:
作动器活塞杆的伸缩量,即位移可由作动器的上、下支点之间的距离减去作动器初始长l0来确定为:
Δli=li-l0(i=1,2)
本发明的有益效果是:
(1)实现前摆芯喷管两通道完全解耦控制,降低喷管负载力矩要求,实现前摆芯喷管在运载型号中的应用;
(2)解析法求解,无数值计算误差,精度高;
(3)采用数字伺服控制器易于实现,计算量小,运算速度快,不影响伺服闭环控制特性。
附图说明
图1是系统坐标系定义图;
图2是系统坐标系含作动器定义图;
图3为本发明给出的方法所使用的系统坐标系和旋转角度绕Y’轴旋转时的情况定义图;
图4为本发明给出的方法所使用的系统坐标系和旋转角度绕X轴旋转时的情况定义图。
具体实施方式
下面结合附图与实施例对本发明进行进一步的介绍:
一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法,包括两台伺服伺服作动器均成90°配置,定义静坐标系O-XYZ以喷管摆心为原点,固连于喷管摆心;X轴、Y轴分别为喷管2个单摆方向;动坐标系在初始状态下与静坐标系O-X″Y″Z″一致,与喷管固连,随着喷管姿态的变化而变化,从喷管尾部看,静坐标系和动坐标系Z(Z″)轴垂直纸面向外;零位时,两个坐标系原点重合,坐标轴方向相平行;
本发明包括以下步骤:
步骤1.根据姿态角得到坐标系转换矩阵;
步骤2.计算上述姿态角状态下作动器的位置;
步骤3.得到作动器活塞杆的伸缩量。
所述步骤一中
若姿态角为α,β为喷管在2个方向上的单摆角度,坐标系转换矩阵为:
所述步骤二中作喷管摆动时,作动器下支点会随着喷管运动,作动器上支点始终保持不动;作动器1、2的上支点坐标分别为(c1,c2,c3)、(d1,d2,d3)
作动器1、2下支点在喷管运动前,相对于动坐标为(a1,a2,a3),(b1,b2,b3),定义下支点计算坐标为:
作动器下支点在喷管运动后,相对于静坐标系使用向量表示为:
则作动器1、2方向分别摆动角度α,β之后两作动器的长度为:
作动器活塞杆的伸缩量,即位移可由作动器的上、下支点之间的距离减去作动器初始长l0来确定为:
Δli=li-l0(i=1,2)
实施例
如图1和图2所示,两台伺服作动器均成90°配置,A、B分别为作动器1、2下支点,C、D分别为作动器1、2上支点。静坐标系O-XYZ以喷管摆心为原点,固连于喷管摆心;X轴、Y轴分别为喷管2个单摆方向。
假设不考虑喷管形变,则作动器下支点与喷管摆心的相对位置是固定的,即在坐标系O-XYZ中下支点的坐标是唯一确定的。当控制系统的指令为(α,β),合成摆动可以分解为两步完成:先沿作动器1运动方向摆动角度α,到达中间过渡位置,再沿作动器2运动方向摆动角度β,得到O-X″Y″Z″。
所以变换矩阵为:
其中α,β的符号由右手定则确定。
作喷管摆动时,作动器下支点会随着喷管运动,作动器上支点始终保持不动。作动器1、2的上支点坐标分别为(c1,c2,c3)、(d1,d2,d3)
作动器1、2下支点在喷管运动前,相对于动坐标为(a1,a2,a3),(b1,b2,b3),定义下支点计算坐标为:
作动器下支点在喷管运动后,相对于静坐标系使用向量表示为:
则作动器1、2方向分别摆动角度α,β之后两作动器的长度为:
作动器活塞杆的伸缩量(即位移)可由作动器的上、下支点之间的距离减去作动器初始长l0来确定为:
Δli=li-l0(i=1,2)。
Claims (4)
1.一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法,其特征于:两台伺服伺服作动器均成90°配置,定义静坐标系O-XYZ以喷管摆心为原点,固连于喷管摆心;X轴、Y轴分别为喷管2个单摆方向;动坐标系在初始状态下与静坐标系O-X″Y″Z″一致,与喷管固连,随着喷管姿态的变化而变化,从喷管尾部看,静坐标系和动坐标系Z(Z″)轴垂直纸面向外;零位时,两个坐标系原点重合,坐标轴方向相平行;
本发明包括以下步骤:
步骤1.根据姿态角得到坐标系转换矩阵;
步骤2.计算上述姿态角状态下作动器的位置;
步骤3.得到作动器活塞杆的伸缩量。
2.如权利要求1所述的一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法,其特征在于:所述步骤一中
若姿态角为α,β为喷管在2个方向上的单摆角度,坐标系转换矩阵为:
。
3.如权利要求1所述的一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法,其特征在于:所述步骤二中作喷管摆动时,作动器下支点会随着喷管运动,作动器上支点始终保持不动;作动器1、2的上支点坐标分别为(c1,c2,c3)、(d1,d2,d3)
作动器1、2下支点在喷管运动前,相对于动坐标为(a1,a2,a3),(b1,b2,b3),定义下支点计算坐标为:
作动器下支点在喷管运动后,相对于静坐标系使用向量表示为:
则作动器1、2方向分别摆动角度α,β之后两作动器的长度为:
。
4.如权利要求1所述的一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法,其特征在于:作动器活塞杆的伸缩量,即位移可由作动器的上、下支点之间的距离减去作动器初始长l0来确定为:
Δli=li-l0 (i=1,2) 。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711035311.1A CN109726358A (zh) | 2017-10-30 | 2017-10-30 | 一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711035311.1A CN109726358A (zh) | 2017-10-30 | 2017-10-30 | 一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109726358A true CN109726358A (zh) | 2019-05-07 |
Family
ID=66291510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711035311.1A Pending CN109726358A (zh) | 2017-10-30 | 2017-10-30 | 一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109726358A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111412083A (zh) * | 2020-04-03 | 2020-07-14 | 北京中科宇航技术有限公司 | 一种确定伺服作动器长度变化的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007225365A (ja) * | 2006-02-22 | 2007-09-06 | Fujitsu Ltd | 姿勢角検出方法および姿勢角検出装置 |
CN104764401A (zh) * | 2015-03-11 | 2015-07-08 | 湖北航天技术研究院总体设计所 | 一种发动机柔性喷管摆角及摆心测量方法 |
CN105116906A (zh) * | 2015-07-17 | 2015-12-02 | 中国空间技术研究院 | 基于向量理论的航天器变轨发动机干扰力矩计算方法 |
-
2017
- 2017-10-30 CN CN201711035311.1A patent/CN109726358A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007225365A (ja) * | 2006-02-22 | 2007-09-06 | Fujitsu Ltd | 姿勢角検出方法および姿勢角検出装置 |
CN104764401A (zh) * | 2015-03-11 | 2015-07-08 | 湖北航天技术研究院总体设计所 | 一种发动机柔性喷管摆角及摆心测量方法 |
CN105116906A (zh) * | 2015-07-17 | 2015-12-02 | 中国空间技术研究院 | 基于向量理论的航天器变轨发动机干扰力矩计算方法 |
Non-Patent Citations (1)
Title |
---|
赵亮 等: "基于牵连补偿算法的喷管摆动控制系统设计", 《固体火箭技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111412083A (zh) * | 2020-04-03 | 2020-07-14 | 北京中科宇航技术有限公司 | 一种确定伺服作动器长度变化的方法 |
CN111412083B (zh) * | 2020-04-03 | 2021-09-14 | 北京中科宇航技术有限公司 | 一种确定伺服作动器长度变化的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110471456B (zh) | 高超声速飞行器俯冲段制导、姿控、变形一体化控制方法 | |
JP4742329B2 (ja) | 浮遊移動体の制御システム | |
CN108363305B (zh) | 基于主动干扰补偿的战术导弹鲁棒过载自驾仪设计方法 | |
CN101804862B (zh) | 一种无人机推力变向装置及其控制方法 | |
CN105865272B (zh) | 一种用于半捷联制导导弹的一体化控制方法 | |
CN105425812B (zh) | 一种基于双模型下的无人机自动着舰轨迹控制方法 | |
TWI448969B (zh) | Three - axis dynamic simulation platform system and its control method | |
CN105116914B (zh) | 一种平流层飞艇解析模型预测路径跟踪控制方法 | |
Jiang et al. | Improving tracking accuracy of a novel 3-DOF redundant planar parallel kinematic machine | |
CN115406312B (zh) | 考虑视场角和舵机延时约束的导弹制导控制一体化方法 | |
CN110641738B (zh) | 一种空间五自由度自由飞行机械臂的轨迹跟踪控制方法 | |
CN104881553A (zh) | 单滑块滚喷模式变质心飞行器模型及其结构布局参数的设计方法 | |
CN109726358A (zh) | 一种基于三维坐标变换的前摆心喷管铰链耦合解耦算法 | |
CN113919081A (zh) | 一种考虑惯性耦合的柔性飞行动力学建模与分析方法 | |
CN103291493B (zh) | 运动解耦的轴对称矢量喷管调节机构 | |
CN110823016B (zh) | 一种转捩研究用高精度三维空间制导方法 | |
CN112327892B (zh) | 一种auv误差受限的抗扰控制方法 | |
CN111547176B (zh) | 自平衡机器人控制方法、系统、自平衡机器人及介质 | |
CN117234086A (zh) | 一种全垫升气垫船航速航向指定时间快速控制方法 | |
CN109571484B (zh) | 一种用于自组装任务的空间机械臂柔顺装配控制方法 | |
RU2392186C2 (ru) | Способ управления двухдвигательным самолетом и система для его осуществления | |
CN107831782B (zh) | 一种系留无人机近地回收位置控制器设计方法 | |
Meng et al. | Force-sensorless contact force control of an aerial manipulator system | |
Blouin et al. | Trajectory optimization of a small airship in a moving fluid | |
Li et al. | Path-following control for multi-axle car-like wheeled mobile robot with nonholonomic constraint |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190507 |
|
RJ01 | Rejection of invention patent application after publication |