极大型空间刚-柔耦合系统多阶模态主动振动抑制方法与流程

1.本发明属于航天器控制技术领域，具体涉及一种多阶模态主动振动抑制方法。


背景技术：

2.随着航天器结构设计趋于大型化和复杂化，大型挠性空间桁架结构应用也更加广泛。由于桁架结构具有尺度大、阻尼小的特点，在轨运行时容易受到干扰产生振动，进而影响航天器的正常工作。为满足稳定性需求，对于空间大尺度刚-柔耦合系统进行振动抑制具有重要意义。
3.对于空间大尺度桁架的振动抑制问题，文献《空间智能桁架的有限时间振动抑制控制》等给出了桁架结构独立模态空间控制方法，但是缺乏刚-柔耦合系统层面的振动控制方法；申请号为cn107194077a的中国发明专利，提出了一种将卫星动力学、振动抑制与姿态控制解析式耦合的计算模型，设计多个控制器分别对刚体和柔性天线进行控制，但模型动力学方程没有统一，无法直观看出系统刚-柔耦合关系，控制器设计较为复杂。


技术实现要素：

4.为了克服现有技术的不足，本发明提供了一种极大型空间刚-柔耦合系统多阶模态主动振动抑制方法，首先建立空间大尺度刚-柔耦合系统动力学模型；再设计模态滤波器；然后设计状态观测器；接下来设计滑模面及趋近律；最终设计终端滑模控制器，完成系统多阶模态振动抑制；本发明在前有发明的基础上，提出一种基于高精度等效动力学模型的终端滑模振动抑制方法，创造性地将刚体姿态信息与模态坐标统一，在完善此类系统动力学模型基础上，可在有限时间内控制系统多阶模态，使系统稳定。
5.本发明解决其技术问题所采用的技术方案包括如下步骤：
6.步骤1：建立空间大尺度刚-柔耦合系统动力学模型；
[0007][0008][0009]
qy＝[q
y1
,q
y2
,q
y3
,......q
yn
]
t
[0010]
qz＝[q
z1
,q
z2
,q
z3
,......q
zn
]
t
[0011]
其中，η是状态变量，m、c、k为系数矩阵，q为广义力矩阵，p为剩余矩阵，qy、qz分别为等效梁y向和z向模态坐标，θ为系统面内姿态角，为系统面外姿态角、γ为系统滚转姿态角，q
yn
表示第n阶y向模态坐标，q
zn
表示第n阶z向模态坐标；n≥1，表示模态阶数；
[0012]
步骤2：设计模态滤波器；
[0013]
ψ
t
＝φ
tm[0014]
其中，ψ表示模态滤波器向量，φ表示振型函数，m表示系统总体质量矩阵；
[0015]
步骤3：构造状态观测器：
[0016][0017][0018]
其中，li，i＝1,2,3为待设计参数，g(.)为针对高频抖振的非线性函数，σ为线性区间长度，a∈(0,1)，sgn(z)为z的符号函数，z为通用参数，计算g(e)用e替换z；e为状态误差；x1表示实际状态变量，相当于空间大尺度刚-柔耦合系统动力学模型中的η，x2表示x1的导数，x3为中间变量；分别表示状态x1、x2、x3的估计值，的估计值，分别表示状态速度的估计值；
[0019]
步骤4：设计滑模面趋近律；
[0020]
设计系统滑模参数为：
[0021][0022]
其中，s表示滑动变量，α1》0，w∈(0,1)，sgn(x)为符号函数；
[0023]
设计系统滑模面趋近律为：
[0024][0025]
其中，为滑动变量的导数，α》1，0《β《1，k1》0，k2》0；
[0026]
步骤5：设计终端滑模控制器，完成系统多阶模态振动抑制；
[0027]
终端滑模控制器为：
[0028][0029][0030][0031]
其中，u为控制输入，ρ为自适应增益，ρ≥1；表示对总干扰的观测误差上界；us表示饱和函数的阈值参数，uf表示等效控制输入。
[0032]
本发明的有益效果如下：
[0033]
1、本发明提出了一种空间大尺度刚柔耦合系统高精度动力学建模方法，创造性的将刚体的姿态信息与模态坐标统一，可以直观描述系统刚-柔耦合规律。
[0034]
2、本发明提出的基于状态观测器的终端滑模控制器，结合动态补偿思想可以实时估计模型总干扰，且在有限时间内完成对系统多阶模态振动抑制。
附图说明
[0035]
图1为本发明方法流程示意图。
[0036]
图2为本发明方法空间大尺度刚-柔耦合系统示意图。
具体实施方式
[0037]
下面结合附图和实施例对本发明进一步说明。
[0038]
本发明的目的是设计一种基于空间大尺度刚-柔耦合系统高精度等效动力学模型的终端滑模振动抑制控制方法，分别设计模态滤波器将系统实际物理坐标转换为模态坐标、状态观测器获得刚体姿态角速度信息和模态加速度信息、终端滑模控制器对各阶模态展开控制。
[0039]
如图1所示，一种极大型空间刚-柔耦合系统多阶模态主动振动抑制方法，包括如下步骤：
[0040]
步骤1：建立空间大尺度刚-柔耦合系统动力学模型；
[0041][0042][0043]
qy＝[q
y1
,q
y2
,q
y3
,......q
yn
]
t
[0044]
qz＝[q
z1
,q
z2
,q
z3
,......q
zn
]
t
[0045]
其中，η是状态变量，m、c、k为系数矩阵，q为广义力矩阵，p为剩余矩阵，qy、qz分别为等效梁y向和z向模态坐标，θ为系统面内姿态角，为系统面外姿态角、γ为系统滚转姿态角，q
yn
表示第n阶y向模态坐标，q
zn
表示第n阶z向模态坐标；n≥1，表示模态阶数；
[0046]
步骤2：设计模态滤波器；
[0047]
ψ
t
＝φ
tm[0048]
其中，ψ表示模态滤波器向量，φ表示振型函数，m表示系统总体质量矩阵；
[0049]
步骤3：构造状态观测器：
[0050][0051][0052]
其中，li，i＝1,2,3为待设计参数，g(.)为针对高频抖振的非线性函数，σ为线性区间长度，a∈(0,1)，sgn(z)为z的符号函数，z为通用参数，计算g(e)用e替换z；e为状态误差；x1表示实际状态变量，相当于空间大尺度刚-柔耦合系统动力学模型中的η，x2表示x1的导数，x3为中间变量；分别表示状态x1、x2、x3的估计值，的估计值，分别表示状态速度的估计值；
[0053]
步骤4：设计滑模面趋近律；
[0054]
设计系统滑模参数为：
[0055][0056]
其中，s表示滑动变量，α1》0，w∈(0,1)，sgn(x)为符号函数；
[0057]
设计系统滑模面趋近律为：
[0058][0059]
其中，为滑动变量的导数，α》1，0《β《1，k1》0，k2》0；
[0060]
步骤5：设计终端滑模控制器，完成系统多阶模态振动抑制；
[0061]
终端滑模控制器为：
[0062][0063][0064][0065]
其中，u为控制输入，ρ为自适应增益，ρ≥1；表示对总干扰的观测误差上界；us表示饱和函数的阈值参数，uf表示等效控制输入。
[0066]
具体实施例：
[0067]
(1)建立空间大尺度刚-柔耦合系统动力学模型：
[0068]
如图2所示，空间大尺度刚柔耦合系统由两侧刚体卫星和中段桁架组成。其中，卫星简化成质点，桁架通过等效梁模型获得结构参数。系统的状态量包括系统轨道面内角θ、面外角滚转角γ、等效梁的模态坐标qy、qz。
[0069]
空间大尺度刚-柔耦合系统动力学方程为：
[0070][0071][0072]
qy＝[q
y1
,q
y2
,q
y3
,......q
yn
]
t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(0.3)
[0073]
qz＝[q
z1
,q
z2
,q
z3
,......q
zn
]
t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(0.4)
[0074]
(2)设计系统模态滤波器
[0075]
根据模态正交性和模态矩阵归一化运算，可以得到模态滤波器为：
[0076]
ψ
t
＝φ
tmꢀꢀ
(0.5)其中，ψ表示模态滤波器向量，φ表示振型函数，m表示系统总体质量矩阵。
[0077]
(3)设计状态观测器
[0078]
令x1＝η，则系统动力学模型为：
[0079]
[0080]
令x3＝-m-1
kx
1-m-1
cx2 m-1
p表示环境干扰和模型摄动组成的系统总干扰，则系统动力学模型可表示为：
[0081][0082]
构造状态观测器为：
[0083][0084][0085]
其中，li(i＝1,2,3)为待设计参数，g(z)为针对高频抖振的非线性函数，σ为线性区间长度，a∈(0,1)，sgn(z)为z的符号函数。
[0086]
(4)设计系统滑模参数为：
[0087][0088]
其中，a1》0，w∈(0,1)，sgn(x)为符号函数。
[0089]
为减少系统抖振，提高滑模面收敛速度，设计系统趋近律为:
[0090][0091]
其中，α》1，0《β《1，k1》0，k2》0
[0092]
(5)终端滑模控制器为：
[0093][0094][0095][0096]
其中，ρ为自适应增益，ρ≥1；表示对总干扰的观测误差上界；us表示饱和函数的阈值参数。
[0097]
根据系统动力学方程，设计李亚普诺夫函数为：其导数为。将公式(1.6)、公式(1.8)和公式(1.12)代入得到：
[0098][0099]
由于故s和d
max
一致且有界。
[0100]
设计李亚普诺夫函数为：其导数为：
[0101][0102]
故系统可以在有限时间内稳定。
[0103]
因此，本发明设计的基于状态观测的终端滑模控制器可以在有限时间内完成对空间大尺度刚-柔耦合系统振动抑制。