第 卷 第 期
年 月
空 间 科 学 学 报 ,
绳系卫星释放及工作态动力学分析 ‘
顾晓勤
中山学院机械系 中
摘 要
考虑绳系卫星系统绳索质量的影响 , 导出三维空间动力学方程 , 运用 极小值
理论求出释放过程中张力优化控制规律 得到最佳释放速度变化曲线 对子星工作阶段进行
动力学分析 , 讨论系统质心沿椭圆轨道运动与姿态运动的祸合 文中附有算例
关键词 绳系卫星一动力学一优化控制一振动
引言
绳系卫星的概念最早由 提出 , 此系统由母星 一绳索 一子星组成 , 美国的 “双
子星座 号 ” 飞船进行了第一次空间绳实验 通过改变绳索的长度 , 可精确地控制微重
力水平 , 完成各种微重力实验 绳系卫星还可用于发电
、 轨道之间货物转运
、 太空垃圾的
收集 、 子星的在轨推进和上大气层研究等
等研究子星与母星间的相互祸合作用 , 文献 提出绳系卫星系统质心沿
椭圆轨道运动时子星释放
, 并讨论质心沿圆轨道运动时绳系卫星系统的摆动运动 ,
等考虑了空间残余大气对绳系卫星系统运动稳定性的影响 网 等研究了通
过长摆动绳传递到子星的速度增量以改变子星轨道 , 对绳索无电机驱动释放中不确定性
进行了灵敏度分析 叫 , 并提出使子星入轨高度最大的绳射技术方案 , 但为简化起见 , 忽
略不计绳索质量 , 并在轨道平面内建立动力学控制方程 , 加以讨论
于绍华等讨论了绳系卫星系统二维平面运动和常规动力学 , 顾晓勤等讨论了绳系卫
星绳索横向振动的纵向优化控制方法 , 本文导出了绳系卫星系统三维空间动力学方程 ,
并且考虑了绳索质量的影响 , 运用优化理论求出张力优化控制规律 , 得到最佳释放速度
对于系统处于绳长恒定的工作状态 , 得到了轨道平面内、 外摆角及绳索张力变化规律 , 讨
论了系统平面天平动 , 研究了系统质心沿椭圆轨道运动时真近点角与姿态角的关系 为阐
明本文所述方法 , 文中附有数值模拟示例
动力学方程
设地球质心为 , 重力场为球对称 以 。 为原点建立惯性坐标系 一 ‘ ,
沿黄道面与赤道平面交线指向春分点 , 轴沿地球极轴 , 轴由右手法则确定
‘ 轴
轨道
国家自然科学基金项 目
一 一 收到原稿 , 一 一 收到修定稿
期 顾晓勤 绳系卫星释放及工作态动力学分析
面倾角 云、 升交点赤径 刀 和轨道角 如图 所示 母星
、
子星 、 绳索 组成的
系统质心 口 , 建立轨道坐标系 口一 , 其中 轴沿当地垂直轴指向天顶 , 轴沿当
地水平轴指向系统质心飞行方向, 轴垂直于轨道平面 设释放 、 回收过程及子星工作
中绳索保持直线 , 母星 、 子星为质点 , 系统主轴坐标系 口一 那 , 轴沿 , 系统在轨道
面内的摆角 沪以及绕 , 轴转动的面外摆角 口定义如图 所示 一 与 一 之间
的方向余弦矩阵为 图
⋯
价 口
一 价
价 尽
价 口
价
价 口
一 口
尽
沿 、 , 、 轴及 、
、
轴单位矢量分别为 落, , 及 葱。
、 。 、 杨 设地心 。
到绳系卫星系统质心 矢径 , , 质量微元 相对于 点的矢径为 蕊 , 微元相对地
心 。 的矢径为
衬 「
。 , 、二
,
⋯
。 ,
户 一 一 毋了 十 户 凡
·
了一
图 坐标系及绳系卫星位形 图 转角 叭 口的定义
价 口
设绳索释放 、 回收控制器安装在母星中 , 绳索密度 。 , 放出绳长 为时间 艺的函数
质量微元 上作用着万有引力 一
拼
,
‘ , 引力参数 拼 ‘ 令
‘ 一 , , ‘ 一 场 , 一 艺
乞
二‘, 由系统质心定义计算可得到
十
,
,
了 乙 一 了 去 一 二 乙
匕
系统角速度在主轴坐标系投影列阵 。 一 哟 口, 口, 哟 用, 系统相对于
空 间 科 学 学 报 卷
质心 的动能 和势能 分别为
二 。 。 饥 一 户 了 ‘ , 、「, 。 、 一 。
一 二
—
一 二 ” 、 、口 寸 甲 ‘ 西 月 寸 , 尸 】
了过
一 上旦竺些燮些
其中
「饥 饥 爪 一 二 ‘ 艺‘ ,
、
了 了矛、 一
”
上 “ ‘ 了
封 气‘少一
—
州广 暇二 一 下二二 尸叫 ‘
、 义 」
设系统质心 沿半径为 , 的圆轨道运动 , 角速度 户二 俩万声
, 绳系卫星系统的几何尺度远
小于轨道运动几何尺度 , 姿态运动与轨道运动之间祸合作用极其微小 , 故在公式 一
推导中忽略其相互影响 记绳索张力为 约, 将式 代入 方程 , 整理得到以轨
道面内外摆角 诚
、
斑 , 绳长 为变量的绳系卫星系统动力学控制方程
其中
毋“ , 珊‘
‘卜 ‘一月户 ” 铸卜 ” ·价一价一 。
,
口卜韶‘
‘
,户‘ , ” 铸, 一 “ ”‘ “ “一
,
竺业黯匹卫‘。卜 里止气群卫些。, 。 十
生 竺竺丝业丝卫全些丝上竺内 ‘ 一二
、 ‘卜 二竺
目
了竺黔‘竺鱼‘
·
张力优化控制
以绳索张力为控制量 , 要求 令 价 , 铸 , 口 , 户 , , , , 记
, 劝 , 外 , 二 , 卯 , , ‘ , 由式 得到系统状态方程
全 二 ,
、 一 ‘二 ”一 箫‘
‘
, ” 毋, 一 ”‘·‘。。·‘,
。一珊‘
‘ ”‘, 一 ” 毋,’ 一 ‘””‘·“一“
,
。 。 一 。 功 口一 夕 夕 必 口 户内 , 一
户 一 户
取性能指标
,
一 百“ ‘
, 甲 一
’ 、
百」
。 “ ‘ 万忠 ‘“‘,
加权矩阵
,
, , ⋯ , , 鳍 ,
期 顾晓勤 绳系卫星释放及工作态动力学分析
、“ 、 、 、
,
、 二 一 , 、 。
构适 函数 月 , 生
’
, 一 百工 ‘ 艺万 人‘ , 里
‘
, 佯随间量 一 ‘ , 入 , “
‘
, 入 ‘
令 。 二 。 一 户 , 运用 , 极小值原理可得到协态方程
“ 一“‘ ‘ ”一 ‘一 号”
·‘·
, ‘ ‘· 十 丛签碧入
,
‘ 一“ 济一 ‘ ‘ 贪‘‘卜 户‘二月 ” 踌卜‘
·
, “ 鱼罕旦
,
‘ 一“ ”一 ‘ 裂男立 ” 功, ”一钊‘ 一 口一 ‘ 鱼篆黔
,
“ 一“ ”一 ‘ ” 毋,‘二“一 “ 入贪‘‘ , “ 箫
,
“ 一“‘‘ , 贪一 贪
’ ‘‘ , 入 ” 毋, ‘ 户, 备 。 。 一 , 一 。‘卜
。‘ 一 竺竺丝竺竺卫业丝业亘竺壑丝鱼、 。
,
入 一【 久 户
一 二 ‘‘ , 贪‘ ” 场, ‘
·户,一 ‘
·
式中 凡一 鲁
, 、
及终了状态 入 ,
户 , 一 一 户
释放初始状态 二 。
二 。 功。, 铸。 , 尽。, 户。 , , 。 ,
入‘ 了 二 , 二 、 , ,
由控制方程得到 入 , 联立状态方程
、
求出绳张力控制规律
, , ⋯ ,
协态方程 , 考虑边界条件 , 采用数值法
释放速度优化控制
以释放速度 为控制量 , 令 二一 冲‘
, 踌 , 口 , , 户 , 由式 得到系统状态方
程
念 二 , , , 如 卯
式中 , 外 , 卯 的定义见式 记 扩 , 二 , 外 , , 即 取性能指标
二
。
百‘ ’ ‘ ’‘ ’ 艺
, 一
、
⋯
、
百 。‘ ’ ‘ ‘ ”’‘ ’ ‘ ‘,
加权矩阵 ’ 鳍 , , , ,
· 一 会二一 。 一 ‘ ‘
’ ‘ 留 ‘
, , , 构造 函数 二 , ,
, , 其中 。 乘子向量 入 入 , 入 , “ , 入 运
用优化理论求绳索释放速度控制规律 , 得到协态方程
空 间 科 学 学 报 卷
二 一石叻 入 乡 。
、
叻一 万入 口 功 月
,
一“ 必一 入 ‘ 贪‘‘ , 一 户‘二口 ” 必,入 ‘
· 口,
尽 沪
一 气又, 丁 十‘ 卢口
一 入 沙 沪“ 汐 价 。 尽,
“应一 “ ” 踌,‘二月一 入 入贪‘‘
,工入
·
入
·
入
。
入了少、
控制方程
入 沪 入 口
释放初始状态 二 。 及终了状态 入 了
二
入 ,
功。 , 功。 , 口。 , 尽。
、二 , , 乞 , , ,
联立正则方程 、 和控制方程 , 考虑边界条件式 , 采用数值法求出绳速控制
规律 , 一次积分可得绳长变化曲线
工作阶段
当子星释放结束后 , 系统开始绕平衡位置作空间天平动 由式 得到系统动力学方
程及张力计算式
功 一 口户沙 必 乡 功。 沪
口 沙 踌“ 。 功乡 口。 , 口
,
,
价 口一 乡 沙 铸“ 口 户
—
、。了
‘
由上式可以看到 绳系卫星系统在轨道平面内
、
外摆角 拭约 、 侧 , 与系统静止工作状态
初始条件即释放终了时刻状态量 二
沪 、 口 的函数 一旦 帅了, 功 ,
内
、
外摆角及绳索张力便可确定
当系统在轨道平面内摆动时 ,
勺 有关 , 与绳索长度无关 绳索张力 为 及
街 , 街 , 确定 , 绳系卫星系统静止工作状态轨道平面
即做平面天平动 积分为
、 , ” ‘· 功
方程 第一式的
几 勺
‘
系统相对轨道坐标系作姿态运动时的相对动能与势能之和由释放终了状态确定 , 即 二
、 。 、 。 、
百 · 价 万夕‘ 、 ‘ ‘ 沪
·
变量曲线族
记 二
几 了
式 给出相平面 帅 , 句 中以 为参
、一 士 。 一 昌六 护
期 顾晓勤 绳系卫星释放及工作态动力学分析
消去时间变量后整理得到 自治形式的一阶微分方程
沪
功
沙 沪
功
上述方程的奇点对应系统的相对平衡位置 对 价。 沪。 , 由于绳系卫星系统主惯量
矩 几 份 几 勺 人 份 , 奇点类型为中心 , 即稳定平衡 系统做微幅天平动的周期
“开 一 八下瓜牙弃若了 , 摆幅 功。 、 鲜 于兰
甲
一 一 “
当绳系卫星系统质心沿椭圆轨道运动时 , 质心运动与姿态运动祸合 设近地点幅角
, 真近点角 , 轨道角 。 将式 代入 方程 , 得到轨道平面内动力
学方程
必十 纂 价一
其中 城
、
创 均为真近点角 , 的函数
亡 ”‘卜福‘ 一“ ,
为半轴参数 , 为偏心率 将式 自变量 ‘置换为真近点角 , 得到周期变系数的非
齐次非线性方程
价
圣
价
—
, 二万 十 了 , 了一一一一一一一 二二叮
, 艺又 以
价
入
当系统做微幅摆动时 , 将上式线性化并做变量置换 梦 。 沪, 得到非齐次 方
程
写
二下二
仃宝
。
了一下一 一一芍犷 夕
十 三 口
对 《 的椭圆轨道 , 将上述方程的解展开成偏心率 的幂级数 , 求解得到 夕 。 一
护 十 ⋯, 回代后得到
一 、 。
、
月 一 、 , 吐 , 山 、
甲 一 一一一二不下一甲叮丁下一 —一
。气
十 口
真近点角 为时间 的函数 , 系统质心椭圆运动周期 二 一探
, · 为长半轴
数值模拟
算例 设母星质量 , 子星 , 绳索线密度 二 一 ,
。 , 系统质心沿 的圆轨道运动 绳索全部释放 , 释放终了 行 ,
衍 一 一“ “
, 街 一 , 舟二 。 得到绳系卫星系统轨道平面内
、
外摆角角速度 拭
、
斑约变化规律如图 所示 绳索张力变化曲线如图 所示 , 图 为
面内、 外摆角 拭约
、
斑 变化规律
空 间 科 学 学 报 卷
之一的舀
一
、甲︵。︶扮﹃劣三瞥代
一
入入入入 八八八 八八八八
厂厂 人 产产
了
‘‘ 、、
狄狄
丫 犷犷
又又又
忿
图 轨道平面内外摆角角速度变化规律
,,,,,,
、、、、 、、、
丁丁丁丁 ‘
廿廿
下下
工工
洲洲阿阿
卜卜卜卜卜
亡
图 绳索张力变化曲线
乃乃乃乃 八八八
了了了 从从 从从从
八八八
从从从 刀刀 人人
代代代代 、、 又又沪沪
甘甘甘甘 粉粉粉 勺””
山
一
人。︶。一助口侧
图 轨道平面内外摆角变化曲线
︵。︶出留侧
月月月口口巨巨口口巨巨口口口口口口口口口口
回回回回回回回回回回回回
洲洲洲洲洲洲洲洲洲瓜瓜瓜瓜
尸尸尸尸尸尸尸尸尸尸尸尸
才才才才
勿勿勿口口口口口口口口口口口口口口日日口口
一、、
丫丫丫丫丫丫丫
又又
︵。︶月留代
亡
图 真近点角曲线
升
图 轨道平面内摆角变化规律
期 顾晓勤 绳系卫星释放及工作态动力学分析
算例 系统质量参数及绳长与上例相同 质心沿椭圆轨道运动 , ,
已 设过近地点时刻 , 得真近点角 约变化规律如图 所示 , 周期
限于篇幅略去 曲线图 轨道平面内绳系卫星系统绕平衡位置摆角 拭 曲线如图
所示
结论
本文在不忽略绳索质量的假设下 , 导出绳系卫星系统三维空间动力学运动方程 利用
极小值理论得到释放子星中张力优化控制规律 , 求出释放速度最佳控制规律
对于系统处于绳长恒定的工作状态 , 进行动力学分析 , 得到轨道平面内、 外摆角及绳索张
力变化规律 , 讨论系统平面天平动 , 研究系统质心轨道运动与姿态运动的祸合 , 得到系统
质心作 运动时真近点角与姿态角的关系 数值计算表明本文方法的有效性 本文
进一步工作是分析绳系卫星系统子星回收动力学与控制技术
参 考 文 献
〕 , 、。
云 , , 一
【」 , 叮 ‘ , ,
一
【」 ,
“ £ , 、 , 夕。 , , 一
」 , 。。 ‘ ,
, 一
【 一
伪
·
£, , 一
于绍华 , 刘强 , 杨林娜 绳系卫星系统二维平面运动和常规动力学
顾晓勤 , 谭朝阳 绳系卫星横向振动的控制方法 空间科学学报 ,
宇航学报 , , 一
, 一
刘延柱 航天器姿态动力学 北京 国防工业出版社 , 一
司刁
空 间 科 学 学 报 卷
那时二 。‘ , 叼 。 , 叼
,
,
,
,
, , ,