CN115903475A - 基于空间太阳能电站地面验证系统对日定向模型控制方法 - Google Patents

基于空间太阳能电站地面验证系统对日定向模型控制方法 Download PDF

Info

Publication number
CN115903475A
CN115903475A CN202211316894.6A CN202211316894A CN115903475A CN 115903475 A CN115903475 A CN 115903475A CN 202211316894 A CN202211316894 A CN 202211316894A CN 115903475 A CN115903475 A CN 115903475A
Authority
CN
China
Prior art keywords
cable
sun
structure platform
aerial structure
angle
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.)
Granted
Application number
CN202211316894.6A
Other languages
English (en)
Other versions
CN115903475B (zh
Inventor
杜莹春
张逸群
段宝岩
陈光达
赵泽
樊冠恒
纪祥飞
王东旭
李欣童
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN202211316894.6A priority Critical patent/CN115903475B/zh
Publication of CN115903475A publication Critical patent/CN115903475A/zh
Application granted granted Critical
Publication of CN115903475B publication Critical patent/CN115903475B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种基于空间太阳能电站地面验证系统对日定向模型控制方法,包括如下步骤:步骤一:根据静力学平衡方程,建立空中结构平台的三索对日定向物理模型;步骤二:选择空中结构平台的初始平衡位置以及平衡姿态;步骤三:根据对日定向要求,计算实时索长,并且调整空中结构平台的姿态;步骤四:根据北斗差分定位装置反馈的数据,进行索长误差补偿调整,实现闭环控制;本发明针对空间太阳能电站地面演示验证系统中空中部分对日定向要求,为其测试与控制这一子系统,轨道运行这一子链路提供的适用于该地面演示验证系统研制方案的三索对日定向模型及其控制方法,实现了实时对日定向的目标,达到精度要求下的闭环控制能力。

Description

基于空间太阳能电站地面验证系统对日定向模型控制方法
技术领域
本发明属于空间太阳能电站地面演示验证系统的对日定向姿态控制技术领域,涉及一种基于空间太阳能电站地面验证系统对日定向模型控制方法。
背景技术
空间太阳能电站自上世纪被首次提出以来,世界各国相继研发了多种方案。为了填补国际上对空间太阳能电站全链路、全系统的工程可行性验证的研究空白,段宝岩院士团队以欧米伽(OMEGA)空间太阳能电站创新方案为理论基础,进行了空间太阳能电站全链路、全系统地面演示验证系统的研制。其中,空中部分包含了轨道运行、聚光、射频微波传能等链路环节,以及光能收集与转换、发射天线这些子系统。空中部分需要完成对日定向功能。
现有的建立在地面,对于指向有要求的系统多采用伺服机构进行指向控制,例如雷达的伺服控制系统,可以根据目标方位,进行相应的俯仰角和方位角的快速调整,但是其俯仰角度调整受限,并不能适用于空间太阳能电站地面验证系统中的对日定向要求。
针对这一问题,设计提出了一种采用索拉控制系统来进行姿态调整,以此来达到对日定向的目的。
发明内容
本发明的目的是提供一种基于空间太阳能电站地面验证系统对日定向模型控制方法,采用该方法能够实现针对空间太阳能电站地面验证系统中空中部分的对日定向的要求。
本发明所采用的技术方案是,基于空间太阳能电站地面验证系统对日定向模型控制方法,具体包括如下步骤:
步骤1,根据静力学平衡方程,建立空中结构平台的三索对日定向物理模型;
步骤2,选择空中结构平台的初始平衡位置以及平衡姿态;
步骤3,根据对日定向要求,计算实时索长,并且调整空中结构平台的姿态;
步骤4,根据北斗差分定位装置反馈的数据,进行索长误差补偿调整,实现闭环控制。
本发明的特点还在于:
步骤1的具体过程为:
步骤1.1,建立坐标系:
以三座吊塔构成的等边三角形几何中心点为坐标原点O,正东方向为X轴,正北方向为Y轴,根据右手定则确定Z轴,建立OXYZ坐标系,确定点坐标A0、B0、C0为已知量,分别表示三个吊塔的顶点坐标;A1、B1、C1三点坐标为已知量,分别表示对应三个卷扬机的三个出索口坐标;A、B、C三点坐标,表示了每一个对日定向姿态下的,空中结构平台上面三个索吊点的位置坐标:
Figure BDA0003909724760000021
Figure BDA0003909724760000031
Figure BDA0003909724760000032
以三个卷扬机的出索口A1、B1、C1为拉索起始端点,则I号卷扬机调控的索长范围表示为A1→A0→A,II号卷扬机调控的索长范围表示为B1→B0→B,3号卷扬机调控的索长范围表示为C1→C0→C,其中,在利用静力学平衡方程构建空中结构平台的三索对日定向物理模型中,仅对索吊点至吊塔顶点部分作为索长求解量,出索口至吊塔顶点的索长为固定量,即定义三根索长L1、L2、L3分别表示为:
Figure BDA0003909724760000033
步骤1.2,确定未知量:
对于整个模型,设置未知量分别是三根索力F1、F2、F3以及三根索长L1、L2、L3,由于三个吊塔顶点A0、B0、C0为已知量,而三个索吊点A、B、C的坐标在姿态调整过程中要随时变化,所以将每次对日定向状态下索长的求解,转化为空中结构平台上面三个索吊点A(xA,yA,zA)、B(xB,yB,zB)、C(xC,yC,zC)的坐标分量求解,因此,空间太阳能电站地面演示验证系统的三索对日定向模型包含12个未知量,建立统一方程形式:
qx=b (5);
其中,未知量为:
x=(xF,xp) (6);
其中,xF=(F1,F2,F3)T,xp=(xA,yA,zA,xB,yB,zB,xC,yC,zC)T,q为系数矩阵,b为常数约束向量;
步骤1.3,建立力平衡和力矩平衡约束关系:
将每一个对日定向状态,看做是静力平衡状态,空中结构平台要满足静力平衡,以及力矩平衡,由于空中结构平台为中心对称结构,重力均匀分布,质心位置与几何中心重合;
对空中结构平台进行静力平衡方程分析如下:
Figure BDA0003909724760000041
矩阵表示为:
Figure BDA0003909724760000042
系数矩阵为:
Figure BDA0003909724760000043
静力平衡方程表示为:
B·xF=[0,0,G]T (10);
对空中结构平台进行力矩平衡方程建立:
Figure BDA0003909724760000051
系数矩阵为:
Figure BDA0003909724760000052
力矩平衡方程可以写为:
W·xP=0 (13);
其中,F1、F2、F3分别对应三个索吊点处的索力大小,G表示空中结构平台的自身重力,α、β、γ分别表示索力在坐标轴X、Y、Z上的方向角,下标1、2、3分别表示索力F1、F2、F3处的方向角;
步骤1.4,建立几何约束关系:
几何约束包括了空中结构平台对日定向的指向姿态约束以及三个吊点之间的固有几何属性约束;
空中结构平台的质心在固定高度h下工作,并且满足空中结构平台的法向量
Figure BDA0003909724760000053
指向太阳,用
Figure BDA0003909724760000054
表示太阳高度角,即太阳光入射线与地平面的夹角,用θ表示太阳方位角,即太阳光入射线在地平面的投影线与正南方向的夹角;
则根据高度约束以及角度约束可以得到一组几何约束方程:
Figure BDA0003909724760000061
其中,其中γn表示法向量
Figure BDA0003909724760000062
与Z坐标轴的夹角,γn和高度角互为余角;αn表示法向量
Figure BDA0003909724760000063
在XOY平面上与X坐标轴的夹角,αn与方位角相等;
提取系数矩阵R:
Figure BDA0003909724760000064
则有:
Figure BDA0003909724760000065
其中,Ri、R'i表示向量xP与角度计算相关的系数;
对于空中结构平台的三个吊点,是一个等边三角形的三个顶点,所以索吊点A到索吊点B的距离、索吊点A到索吊点C的距离以及索吊点B到索吊点C的距离互相保持一致,设计边长为l,则几何约束方程表示为:
Figure BDA0003909724760000066
提取系数矩阵E:
Figure BDA0003909724760000067
其中,Ei、E'i、E″i分别表示向量xP与距离l计算相关的系数;
则有:
E·xP=[l,l,l]T (19);
步骤1.5,建立空中结构平台的三索对日定向物理模型:
联立(10)、(13)、(16)、(19)四个方程,可以得到:
Figure BDA0003909724760000071
结合方程(5),构建统一方程形式:
Figure BDA0003909724760000072
Figure BDA0003909724760000073
x=(F1,F2,F3,xA,xB,xC,yA,yB,yC,zA,zB,zC)T (23);
设空中结构平台质心:
Figure BDA0003909724760000074
结合方程(4)、(24),将空中结构平台在三个坐标轴X、Y、Z方向上的姿态角ηx、ηy、ηz,索长L1、L2、L3,以及空中结构平台质心坐标,这些最终输出量与未知量x之间的关系矩阵写作U,则有:
[L,Kcxyz]T=Ux (25);
联立方程(5)、(25)可得空中结构平台的三索对日定向物理模型:
[L,Kcxyz]T=Uq-1b (26);
其中,L=[L1,L2,L3]。
步骤2中,设定空中结构平台的初始位姿如下:
Figure BDA0003909724760000081
整个空中结构平台处于水平姿态,空中结构质心位置高度固定为h0、三根索长长度相等分别为L0
步骤3的具体过程为:
步骤3.1,实时计算太阳高度角和方位角:
Figure BDA0003909724760000082
Figure BDA0003909724760000083
其中,
Figure BDA0003909724760000084
表示太阳高度角,θ表示太阳方位角,φ表示当地的地理纬度,δ表示太阳赤纬,t表示时角;
步骤3.2,根据步骤3.1建立的空中结构平台的三索对日定向物理模型,将太阳高度角
Figure BDA0003909724760000085
和方位角θ输入进去,计算出此时的理论空中结构姿态数据,包括空中结构平台质心位置Kc,在三个坐标轴X、Y、Z方向上的姿态角ηx、ηy、ηz以及三根索长L=[L1 0,L2 0,L3 0];
步骤3.3,根据步骤3.2得出的一组索长和步骤2得到的一组初始状态下的索长数据,进行作差计算得到本次对日定向状态下的三根索分别对应的索长调整量:
Figure BDA0003909724760000086
步骤3.4,对进行空中结构平台的三索对日定向物理模型进行修正,具体为:
当考虑热膨胀影响的时候,需要将绳子原长按照从出索口开始计算,即考虑温差为ΔT,热膨胀系数为δ,那么由此得到的变形量ΔLT表示为:
Figure BDA0003909724760000091
钢索受其中张力的影响,会发生弹性变形,贯穿全段,变形量ΔLE为:
Figure BDA0003909724760000092
考虑卷扬机回差,经过实际测量三个卷扬机回差值均为ΔLh
根据步骤3.3所述,在得到的三根索长调整量基础上,进行索长调整量修正:
Figure BDA0003909724760000093
即在步骤3.3的基础上,再分别对应调整三根索长的调整量为ΔL1'、ΔL2'、ΔL3'。
步骤4的具体过程为:
步骤4.1,根据步骤3.2所述,由输入的太阳高度角
Figure BDA0003909724760000094
以及太阳方位角θ,经过空中结构平台的三索对日定向物理模型计算,得到三个坐标轴X、Y、Z方向上的理论姿态角ηx、ηy、ηz,再由北斗差分定位装置反馈的数据,得到空中结构平台在该时刻下的实测姿态角ηx'、ηy'、ηz',判断是否满足精度要求|ηx'-ηx|≤Δx,|ηy'-ηy|≤Δy,|ηz'-ηz|≤Δz,其中Δx、Δy、Δz分别对应目标姿态角的误差容许值;若满足,则回到步骤3.2,进行下一时刻对日定向姿态调整,若不满足,则进行步骤4.2补偿调整;
步骤4.2,由北斗差分定位装置反馈的数据,得到了一组实测姿态角数据ηx'、ηy'、ηz',再联立方程(21)、(22)、(26),代入姿态角实测值ηx'、ηy'、ηz',便可以计算出,此时刻下,三根钢索的实际长度La=[L1 a,L2 a,L3 a];再根据步骤3.2所述,得到的理论三根钢索长度L=[L1 0,L2 0,L3 0],两组数据进行作差计算得到三根钢索分别对应的索长补偿量为:
Figure BDA0003909724760000101
根据步骤3.4所述,在索长修正调整后的基础上,进一步进行索长补偿调整。然后返回到步骤4.1。
本发明的有益效果如下:
1.建立了空中结构平台的三索对日定向物理模型,实现了不同时刻的对日定向姿态下,计算出三根索长调整量的能力;
2.基于迭代计算,闭环控制的思想,利用姿态信息反馈装置,实现了误差范围允许内的最佳索长调整策略,即满足了对日定向姿态调整指标,完成了精密控制的要求。
附图说明
图1是本发明基于空间太阳能电站地面验证系统对日定向模型控制方法中建立的空中结构平台的三索对日定向物理模型的主流程图;
图2是本发明基于空间太阳能电站地面验证系统对日定向模型控制方法中建立的空中结构平台的三索对日定向物理模型流程图;
图3是本发明基于空间太阳能电站地面验证系统对日定向模型控制方法中三索对日定向模型系统组成与全局坐标系图;
图4(a)、(b)是本发明基于空间太阳能电站地面验证系统对日定向模型控制方法中空中结构受力分析图;
图5(a)、(b)是本发明基于空间太阳能电站地面验证系统对日定向模型控制方法中指向约束示意图;
图6是本发明基于空间太阳能电站地面验证系统对日定向模型控制方法中计算实时索长且进行姿态调整流程图;
图7是本发明基于空间太阳能电站地面验证系统对日定向模型控制方法中索长闭环控制调整流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明基于空间太阳能电站地面验证系统对日定向模型控制方法,如图1所示,包括吊塔、卷扬机、空中结构平台以及北斗定位装置,其具体工作步骤以及控制方法如下:
吊塔共有三个,用于承重,每个吊塔结构相同,且在塔顶、塔身以及底部塔座均安装有定滑轮,用于固定钢索滑道。
卷扬机共有三个,构成索拉控制系统,根据实时计算得到的索长数据,分别调整对应三根钢索的长度。
空中结构平台,作为被控制对象,与三根钢索相连接,完成对日定向的工作要求。
北斗定位装置共有两个,安装在空间太阳能电站地面演示验证系统的空中结构平台上面,用于实时监测空中部分的姿态角度,反馈数据,进行闭环控制。
本发明基于空间太阳能电站地面验证系统对日定向模型控制方法,包括如下步骤:
步骤1,根据静力学平衡方程,建立空中结构平台的三索对日定向物理模型;
如图2所示,步骤1的具体过程为:
步骤1.1,建立坐标系:
如图3所示,以三座吊塔构成的等边三角形几何中心点为坐标原点O,正东方向为X轴,正北方向为Y轴,根据右手定则确定Z轴,建立OXYZ坐标系。确定点坐标A0、B0、C0为已知量,分别表示三个吊塔的顶点坐标;A1、B1、C1三点坐标为已知量,分别表示对应三个卷扬机的三个出索口坐标;A、B、C三点坐标,表示了每一个对日定向姿态下的,空中结构平台上面三个索吊点的位置坐标。
Figure BDA0003909724760000121
Figure BDA0003909724760000131
Figure BDA0003909724760000132
以三个卷扬机的出索口A1、B1、C1为拉索起始端点,则I号卷扬机调控的索长范围表示为A1→A0→A,II号卷扬机调控的索长范围表示为B1→B0→B,III号卷扬机调控的索长范围表示为C1→C0→C。其中,在利用静力学平衡方程构建空中结构平台的三索对日定向物理模型中,仅对索吊点至吊塔顶点部分作为索长求解量,出索口至吊塔顶点的索长为固定量。即定义三根索长L1、L2、L3分别表示为:
Figure BDA0003909724760000133
步骤1.2,确定未知量:
对于整个模型,设置未知量分别是三根索力F1、F2、F3以及三根索长L1、L2、L3。由于三个吊塔顶点A0、B0、C0为已知量,而三个索吊点A、B、C的坐标在姿态调整过程中要随时变化,
所以将每次对日定向状态下索长的求解,转化为空中结构平台上面三个索吊点A(xA,yA,zA)、B(xB,yB,zB)、C(xC,yC,zC)的坐标分量求解。
因此,空间太阳能电站地面演示验证系统的三索对日定向模型包含12个未知量,建立统一方程形式:
qx=b (5);
其中,未知量为:
x=(xF,xp) (6);
其中xF=(F1,F2,F3)T,xp=(xA,yA,zA,xB,yB,zB,xC,yC,zC)T,q为系数矩阵,b为常数约束向量。
步骤1.3,建立力平衡和力矩平衡约束关系,具体为:
如图4(a)、(b)所示(图4(a)为三根索力与空中结构平台自身重力G的分布状态示意图,图4(b)为三根所里在坐标轴上的方向角示意图),可以将每一个对日定向状态,看做是静力平衡状态,空中结构平台要满足静力平衡,以及力矩平衡。由于空中结构平台为中心对称结构,重力均匀分布,质心位置与几何中心重合。
对空中结构平台进行静力平衡方程分析如下:
Figure BDA0003909724760000141
矩阵表示为:
Figure BDA0003909724760000142
系数矩阵为:
Figure BDA0003909724760000143
静力平衡方程可以写为:
B·xF=[0,0,G]T (10);
对空中结构平台进行力矩平衡方程建立:
Figure BDA0003909724760000151
系数矩阵为:
Figure BDA0003909724760000152
力矩平衡方程可以写为:
W·xP=0 (13);
其中,F1、F2、F3分别对应三个索吊点处的索力大小,G表示空中结构平台的自身重力,α、β、γ分别表示索力在坐标轴X、Y、Z上的方向角,下标1、2、3分别表示索力F1、F2、F3处的方向角。
步骤1.4,建立几何约束关系,具体为:
几何约束包括了空中结构平台对日定向的指向姿态约束以及三个吊点之间的固有几何属性约束。
如图5(a)、5(b)所示,空中结构平台的质心需要在固定高度h下工作,并且需要满足空中结构平台的法向量
Figure BDA0003909724760000153
指向太阳,用
Figure BDA0003909724760000154
表示太阳高度角,即太阳光入射线与地平面的夹角,用θ表示太阳方位角,参见图5(a),即太阳光入射线在地平面的投影线与正南方向的夹角,图5(b)为三个索吊点A、B、C的坐标与法向量
Figure BDA0003909724760000155
的关系图。
则根据高度约束以及角度约束可以得到一组几何约束方程:
Figure BDA0003909724760000161
其中,其中γn表示法向量
Figure BDA0003909724760000162
与Z坐标轴的夹角,它和高度角互为余角;αn表示法向量
Figure BDA0003909724760000163
在XOY平面上与X坐标轴的夹角,它与方位角相等。
提取系数矩阵R:
Figure BDA0003909724760000164
则有:
Figure BDA0003909724760000165
其中,Ri、R'i表示向量xP与角度计算相关的系数。
对于空中结构平台的三个吊点,它们是一个等边三角形的三个顶点,所以索吊点A到索吊点B的距离、索吊点A到索吊点C的距离以及索吊点B到索吊点C的距离互相保持一致,设计边长为l,则几何约束方程可以写为:
Figure BDA0003909724760000166
提取系数矩阵E:
Figure BDA0003909724760000167
其中,Ei、E'i、E″i分别表示向量xP与距离l计算相关的系数。
则有:
E·xP=[l,l,l]T (19);
步骤1.5,建立空中结构平台的三索对日定向物理模型,具体为:
联立(10)、(13)、(16)、(19)等四个方程,可以得到:
Figure BDA0003909724760000171
结合方程(5),构建统一方程形式:
Figure BDA0003909724760000172
Figure BDA0003909724760000173
x=(F1,F2,F3,xA,xB,xC,yA,yB,yC,zA,zB,zC)T (23);
设空中结构平台质心:
Figure BDA0003909724760000174
结合方程(4)、(24),将空中结构平台在三个坐标轴X、Y、Z方向上的姿态角ηx、ηy、ηz,索长L1、L2、L3,以及空中结构平台质心坐标,这些最终输出量与未知量x之间的关系矩阵写作U,则有:
[L,Kcxyz]T=Ux (25);
联立方程(5)、(25)可得空中结构平台的三索对日定向物理模型:
[L,Kcxyz]T=Uq-1b (26);
其中,L=[L1,L2,L3]。
步骤2,选择空中结构平台的初始平衡位置以及平衡姿态;
步骤2中,设定空中结构平台的初始位姿:
Figure BDA0003909724760000181
整个空中结构平台处于水平姿态,空中结构质心位置高度固定为h0、三根索长长度相等分别为L0
步骤3,根据对日定向要求,计算实时索长,并且调整空中结构平台的姿态;
如图6所示,步骤3的具体过程为:
步骤3.1,实时计算计算太阳高度角和方位角:
Figure BDA0003909724760000182
Figure BDA0003909724760000183
其中,
Figure BDA0003909724760000184
表示太阳高度角,θ表示太阳方位角,φ表示当地的地理纬度,δ表示太阳赤纬,t表示时角。
步骤3.2,根据步骤1建立的空中结构平台的三索对日定向物理模型,将太阳高度角
Figure BDA0003909724760000185
和方位角θ输入进去,计算出此时的理论空中结构姿态数据,包括空中结构平台质心位置Kc,在三个坐标轴X、Y、Z方向上的姿态角ηx、ηy、ηz以及三根索长L=[L1 0,L2 0,L3 0]。
步骤3.3,根据步骤3.2得出的一组索长,和步骤二里面得到的一组初始状态下的索长数据,进行作差计算得到本次对日定向状态下的三根索分别对应的索长调整量。
Figure BDA0003909724760000191
步骤3.4,综合考虑误差来源,包括热膨胀、弹性变形、回差等因素,进行空中结构平台的三索对日定向物理模型修正。
当考虑热膨胀影响的时候,需要将绳子原长按照从出索口开始计算,即考虑温差为ΔT,热膨胀系数为δ,那么由此得到的变形量ΔLT表示为:
Figure BDA0003909724760000192
钢索受其中张力的影响,会发生弹性变形,贯穿全段,变形量ΔLE为:
Figure BDA0003909724760000193
考虑卷扬机回差,经过实际测量三个卷扬机回差值均为ΔLh
根据步骤3.3所述,在得到的三根索长调整量基础上,进行索长调整量修正:
Figure BDA0003909724760000194
即在步骤3.3的基础上,再分别对应调整三根索长的调整量为ΔL1'、ΔL2'、ΔL3'。
步骤4,根据北斗差分定位装置反馈的数据,进行索长误差补偿调整,实现闭环控制。
如图7所示,步骤4的具体过程为:
步骤4.1,根据步骤3.2所述,由输入的太阳高度角
Figure BDA0003909724760000201
以及太阳方位角θ,经过空中结构平台的三索对日定向物理模型计算,得到三个坐标轴X、Y、Z方向上的理论姿态角ηx、ηy、ηz,再由北斗差分定位装置反馈的数据,得到空中结构平台在该时刻下的实测姿态角ηx'、ηy'、ηz',判断是否满足精度要求|ηx'-ηx|≤Δx,|ηy'-ηy|≤Δy,|ηz'-ηz|≤Δz,其中Δx、Δy、Δz分别对应目标姿态角的误差容许值。若满足,则回到步骤3.2,进行下一时刻对日定向姿态调整。若不满足,则进行下一步补偿调整。
步骤4.2,由北斗差分定位装置反馈的数据,得到了一组实测姿态角数据ηx'、ηy'、ηz',再联立方程(21)、(22)、(26),代入姿态角实测值ηx'、ηy'、ηz',便可以计算出,此时刻下,三根钢索的实际长度La=[L1 a,L2 a,L3 a]。再根据步骤3.2所述,得到的理论三根钢索长度L=[L1 0,L2 0,L3 0],两组数据进行作差计算得到三根钢索分别对应的索长补偿量为:
Figure BDA0003909724760000202
根据步骤3.4所述,在其索长修正调整后的基础上,进一步进行索长补偿调整。然后返回到步骤4.1。

Claims (5)

1.基于空间太阳能电站地面验证系统对日定向模型控制方法,其特征在于,具体包括如下步骤:
步骤1,根据静力学平衡方程,建立空中结构平台的三索对日定向物理模型;
步骤2,选择空中结构平台的初始平衡位置以及平衡姿态;
步骤3,根据对日定向要求,计算实时索长,并且调整空中结构平台的姿态;
步骤4,根据北斗差分定位装置反馈的数据,进行索长误差补偿调整,实现闭环控制。
2.根据权利要求1所述的基于空间太阳能电站地面验证系统对日定向模型控制方法,其特征在于,所述步骤1的具体过程为:
步骤1.1,建立坐标系:
以三座吊塔构成的等边三角形几何中心点为坐标原点O,正东方向为X轴,正北方向为Y轴,根据右手定则确定Z轴,建立OXYZ坐标系,确定点坐标A0、B0、C0为已知量,分别表示三个吊塔的顶点坐标;A1、B1、C1三点坐标为已知量,分别表示对应三个卷扬机的三个出索口坐标;A、B、C三点坐标,表示了每一个对日定向姿态下的,空中结构平台上面三个索吊点的位置坐标:
Figure FDA0003909724750000011
Figure FDA0003909724750000021
Figure FDA0003909724750000022
以三个卷扬机的出索口A1、B1、C1为拉索起始端点,则I号卷扬机调控的索长范围表示为A1→A0→A,II号卷扬机调控的索长范围表示为B1→B0→B,III号卷扬机调控的索长范围表示为C1→C0→C,其中,在利用静力学平衡方程构建空中结构平台的三索对日定向物理模型中,仅对索吊点至吊塔顶点部分作为索长求解量,出索口至吊塔顶点的索长为固定量,即定义三根索长L1、L2、L3分别表示为:
Figure FDA0003909724750000023
步骤1.2,确定未知量:
对于整个模型,设置未知量分别是三根索力F1、F2、F3以及三根索长L1、L2、L3,由于三个吊塔顶点A0、B0、C0为已知量,而三个索吊点A、B、C的坐标在姿态调整过程中要随时变化,所以将每次对日定向状态下索长的求解,转化为空中结构平台上面三个索吊点A(xA,yA,zA)、B(xB,yB,zB)、C(xC,yC,zC)的坐标分量求解,因此,空间太阳能电站地面演示验证系统的三索对日定向模型包含12个未知量,建立统一方程形式:
qx=b (5);
其中,未知量为:
x=(xF,xp) (6);
其中,xF=(F1,F2,F3)T,xp=(xA,yA,zA,xB,yB,zB,xC,yC,zC)T,q为系数矩阵,b为常数约束向量;
步骤1.3,建立力平衡和力矩平衡约束关系:
将每一个对日定向状态,看做是静力平衡状态,空中结构平台要满足静力平衡,以及力矩平衡,由于空中结构平台为中心对称结构,重力均匀分布,质心位置与几何中心重合;
对空中结构平台进行静力平衡方程分析如下:
Figure FDA0003909724750000031
矩阵表示为:
Figure FDA0003909724750000032
系数矩阵为:
Figure FDA0003909724750000033
静力平衡方程表示为:
B·xF=[0,0,G]T (10);
对空中结构平台进行力矩平衡方程建立:
Figure FDA0003909724750000041
系数矩阵为:
Figure FDA0003909724750000042
力矩平衡方程可以写为:
W·xP=0 (13);
其中,F1、F2、F3分别对应三个索吊点处的索力大小,G表示空中结构平台的自身重力,α、β、γ分别表示索力在坐标轴X、Y、Z上的方向角,下标1、2、3分别表示索力F1、F2、F3处的方向角;
步骤1.4,建立几何约束关系:
几何约束包括了空中结构平台对日定向的指向姿态约束以及三个吊点之间的固有几何属性约束;
空中结构平台的质心在固定高度h下工作,并且满足空中结构平台的法向量
Figure FDA0003909724750000043
指向太阳,用
Figure FDA0003909724750000044
表示太阳高度角,即太阳光入射线与地平面的夹角,用θ表示太阳方位角,即太阳光入射线在地平面的投影线与正南方向的夹角;
则根据高度约束以及角度约束可以得到一组几何约束方程:
Figure FDA0003909724750000051
其中,其中γn表示法向量
Figure FDA0003909724750000052
与Z坐标轴的夹角,γn和高度角互为余角;αn表示法向量
Figure FDA0003909724750000053
在XOY平面上与X坐标轴的夹角,αn与方位角相等;
提取系数矩阵R:
Figure FDA0003909724750000054
则有:
Figure FDA0003909724750000055
其中,Ri、R'i表示向量xP与角度计算相关的系数;
对于空中结构平台的三个吊点,是一个等边三角形的三个顶点,所以索吊点A到索吊点B的距离、索吊点A到索吊点C的距离以及索吊点B到索吊点C的距离互相保持一致,设计边长为l,则几何约束方程表示为:
Figure FDA0003909724750000056
提取系数矩阵E:
Figure FDA0003909724750000057
其中,Ei、E'i、E”i分别表示向量xP与距离l计算相关的系数;则有:
E·xP=[l,l,l]T (19);
步骤1.5,建立空中结构平台的三索对日定向物理模型:
联立(10)、(13)、(16)、(19)四个方程,可以得到:
Figure FDA0003909724750000061
结合方程(5),构建统一方程形式:
Figure FDA0003909724750000062
Figure FDA0003909724750000063
x=(F1,F2,F3,xA,xB,xC,yA,yB,yC,zA,zB,zC)T (23);
设空中结构平台质心:
Figure FDA0003909724750000064
结合方程(4)、(24),将空中结构平台在三个坐标轴X、Y、Z方向上的姿态角ηx、ηy、ηz,索长L1、L2、L3,以及空中结构平台质心坐标,这些最终输出量与未知量x之间的关系矩阵写作U,则有:
[L,Kcxyz]T=Ux (25);
联立方程(5)、(25)可得空中结构平台的三索对日定向物理模型:
[L,Kcxyz]T=Uq-1b (26);
其中,L=[L1,L2,L3]。
3.根据权利要求2所述的基于空间太阳能电站地面验证系统对日定向模型控制方法,其特征在于:所述步骤2中,设定空中结构平台的初始位姿如下:
Figure FDA0003909724750000071
整个空中结构平台处于水平姿态,空中结构质心位置高度固定为h0、三根索长长度相等分别为L0
4.根据权利要求3所述的基于空间太阳能电站地面验证系统对日定向模型控制方法,其特征在于:所述步骤3的具体过程为:
步骤3.1,实时计算太阳高度角和方位角:
Figure FDA0003909724750000072
Figure FDA0003909724750000073
其中,
Figure FDA0003909724750000074
表示太阳高度角,θ表示太阳方位角,φ表示当地的地理纬度,δ表示太阳赤纬,t表示时角;
步骤3.2,根据步骤3.1建立的空中结构平台的三索对日定向物理模型,将太阳高度角
Figure FDA0003909724750000075
和方位角θ输入进去,计算出此时的理论空中结构姿态数据,包括空中结构平台质心位置Kc,在三个坐标轴X、Y、Z方向上的姿态角ηx、ηy、ηz以及三根索长L=[L1 0,L2 0,L3 0];
步骤3.3,根据步骤3.2得出的一组索长和步骤2得到的一组初始状态下的索长数据,进行作差计算得到本次对日定向状态下的三根索分别对应的索长调整量:
Figure FDA0003909724750000081
步骤3.4,对进行空中结构平台的三索对日定向物理模型进行修正,具体为:
当考虑热膨胀影响的时候,需要将绳子原长按照从出索口开始计算,即考虑温差为ΔT,热膨胀系数为δ,那么由此得到的变形量ΔLT表示为:
Figure FDA0003909724750000082
钢索受其中张力的影响,会发生弹性变形,贯穿全段,变形量ΔLE为:
Figure FDA0003909724750000083
考虑卷扬机回差,经过实际测量三个卷扬机回差值均为ΔLh
根据步骤3.3所述,在得到的三根索长调整量基础上,进行索长调整量修正:
Figure FDA0003909724750000084
即在步骤3.3的基础上,再分别对应调整三根索长的调整量为ΔL1'、ΔL2'、ΔL3'。
5.根据权利要求4所述的基于空间太阳能电站地面验证系统对日定向模型控制方法,其特征在于:所述步骤4的具体过程为:
步骤4.1,根据步骤3.2所述,由输入的太阳高度角
Figure FDA0003909724750000091
以及太阳方位角θ,经过空中结构平台的三索对日定向物理模型计算,得到三个坐标轴X、Y、Z方向上的理论姿态角ηx、ηy、ηz,再由北斗差分定位装置反馈的数据,得到空中结构平台在该时刻下的实测姿态角ηx'、ηy'、ηz',判断是否满足精度要求|ηx'-ηx|≤Δx,|ηy'-ηy|≤Δy,|ηz'-ηz|≤Δz,其中Δx、Δy、Δz分别对应目标姿态角的误差容许值;若满足,则回到步骤3.2,进行下一时刻对日定向姿态调整,若不满足,则进行步骤4.2补偿调整;
步骤4.2,由北斗差分定位装置反馈的数据,得到了一组实测姿态角数据ηx'、ηy'、ηz',再联立方程(21)、(22)、(26),代入姿态角实测值ηx'、ηy'、ηz',便可以计算出,此时刻下,三根钢索的实际长度La=[L1 a,L2 a,L3 a];再根据步骤3.2所述,得到的理论三根钢索长度L=[L1 0,L2 0,L3 0],两组数据进行作差计算得到三根钢索分别对应的索长补偿量为:
Figure FDA0003909724750000092
根据步骤3.4所述,在索长修正调整后的基础上,进一步进行索长补偿调整,然后返回到步骤4.1。
CN202211316894.6A 2022-10-26 2022-10-26 基于空间太阳能电站地面验证系统对日定向模型控制方法 Active CN115903475B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211316894.6A CN115903475B (zh) 2022-10-26 2022-10-26 基于空间太阳能电站地面验证系统对日定向模型控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211316894.6A CN115903475B (zh) 2022-10-26 2022-10-26 基于空间太阳能电站地面验证系统对日定向模型控制方法

Publications (2)

Publication Number Publication Date
CN115903475A true CN115903475A (zh) 2023-04-04
CN115903475B CN115903475B (zh) 2023-09-15

Family

ID=86486345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211316894.6A Active CN115903475B (zh) 2022-10-26 2022-10-26 基于空间太阳能电站地面验证系统对日定向模型控制方法

Country Status (1)

Country Link
CN (1) CN115903475B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
UA72511U (ru) * 2011-12-26 2012-08-27 Виктор Васильевич Ожинский Способ контроля солнечной ориентации фотоэлектрических преобразований космических аппаратов
CN105905317A (zh) * 2016-06-07 2016-08-31 湖北航天技术研究院总体设计所 一种卫星对日定向控制系统及其控制方法
CN106155074A (zh) * 2016-08-15 2016-11-23 上海航天控制技术研究所 一种保证星地链路的卫星三轴对日定向控制方法
CN111007865A (zh) * 2019-12-18 2020-04-14 中国人民解放军国防科技大学 以对日指向偏差为约束的卫星平稳对地定向方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
UA72511U (ru) * 2011-12-26 2012-08-27 Виктор Васильевич Ожинский Способ контроля солнечной ориентации фотоэлектрических преобразований космических аппаратов
CN105905317A (zh) * 2016-06-07 2016-08-31 湖北航天技术研究院总体设计所 一种卫星对日定向控制系统及其控制方法
CN106155074A (zh) * 2016-08-15 2016-11-23 上海航天控制技术研究所 一种保证星地链路的卫星三轴对日定向控制方法
CN111007865A (zh) * 2019-12-18 2020-04-14 中国人民解放军国防科技大学 以对日指向偏差为约束的卫星平稳对地定向方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
方元坤 等: "分布式遥感编队多星协同观测中的姿态控制", 光学精密工程, vol. 27, no. 1, pages 58 - 68 *

Also Published As

Publication number Publication date
CN115903475B (zh) 2023-09-15

Similar Documents

Publication Publication Date Title
CN110308746B (zh) 一种适用于三轴转台式测控天线的星体标校方法
CN106985139B (zh) 基于扩展状态观测与补偿的空间机器人自抗扰协调控制方法
Ren et al. Underactuated control and analysis of single blade installation using a jackup installation vessel and active tugger line force control
CN110027678A (zh) 一种基于海浪主动补偿的混联登乘机构运动规划方法
CN111259097A (zh) 一种应用于光伏产业无人机巡检的精细化航点校核方法
CN105015800B (zh) 航天器舱段在地面上的自动化装配系统
CN102568034B (zh) 空间光学遥感器对实际地物成像的计算机仿真系统
Jiang et al. Optimization and analysis on cable net structure supporting the reflector of the large radio telescope FAST
CN112769466B (zh) 低轨卫星星座构型保持方法
CN113641182B (zh) 星间激光通信系统的高精度瞄准指向方法及系统
CN111830845A (zh) 一种无人机风机叶片自动巡检仿真系统及方法
CN114018161A (zh) 一种用于测控系统的高精度距离标校方法
CN108897239A (zh) 一种航天器两级姿态控制模拟系统
CN107589431A (zh) 一种提高机载光电系统目标定位精度的校靶方法
CN105511500B (zh) 一种空间目标动态光学特性地面模拟中目标和光源控制模拟系统及方法
US7938364B2 (en) Apparatus and method for making optimal use of a photovoltaic array on an airborne power system
CN109543227B (zh) 大型反射面天线的风致指向误差实时测量方法、雷达天线
CN115903475B (zh) 基于空间太阳能电站地面验证系统对日定向模型控制方法
CN114221111A (zh) 一种机载卫通天线安装误差动态修正方法
US20240118694A1 (en) Intelligent course planning method and controller for unmanned surface vehicle
CN113674391A (zh) 一种升降轨成像双星地物表观辐亮度差异计算方法
CN112208798A (zh) 一种绕飞编队高码速率星间链路切换方法及系统
CN112015201B (zh) 一种基于预测校正的四旋翼飞行器位置控制方法
CN109213190A (zh) 一种日心悬浮轨道上电动帆航天器编队飞行协同控制方法
Qiu et al. Robust hybrid Controller design based on feedback linearization and µ synthesis for UAV

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
GR01 Patent grant
GR01 Patent grant