空間站柔性太陽(yáng)電池翼對(duì)日跟蹤控制

趙 真，王 碧，錢志源，陳國(guó)平

（1.南京航空航天大學(xué)機(jī)械結(jié)構(gòu)力學(xué)及控制國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇南京 210016；2.上海宇航系統(tǒng)工程研究所，上海 201109）

0 引言

我國(guó)載人航天工程已進(jìn)入空間站研制階段，計(jì)劃在2022 年前后全面建成?？臻g站采用三艙“T 型”布局，分別為核心艙、實(shí)驗(yàn)艙Ⅰ和實(shí)驗(yàn)艙Ⅱ。在核心艙兩端分別留有載人飛船和貨運(yùn)飛船對(duì)接口，如圖1所示［1］。

圖1 中國(guó)空間站構(gòu)型Fig.1 Schematic diagram of Chinese Space Station（CSS）

為保證我國(guó)空間站15 a 在軌運(yùn)行期間能源平衡，在實(shí)驗(yàn)艙Ⅰ、Ⅱ末端各配置1 套Alpha 對(duì)日定向裝置和2 套太陽(yáng)電池翼，如圖2 所示。為適應(yīng)空間站90 min 每圈的軌道周期，Alpha 對(duì)日定向裝置持續(xù)驅(qū)動(dòng)太陽(yáng)電池翼和末端桁架組成的柔性負(fù)載以0.065（°）/s 的角速度對(duì)日跟蹤，將太陽(yáng)能轉(zhuǎn)化為電能。

圖2 空間站太陽(yáng)電池翼驅(qū)動(dòng)控制系統(tǒng)Fig.2 Schematic diagram of the solar cell wing drive control system of CSS

由于柔性負(fù)載的大尺寸（翼展大于55 m，面積約270 m2）、大慣量（3.4×105kg·m2）、超低頻（基頻小于0.07 Hz）、模態(tài)密集（2 Hz 以內(nèi)頻率超過(guò)600 階）等特點(diǎn)，為柔性負(fù)載地面模態(tài)試驗(yàn)、數(shù)學(xué)建模、模型正確性估計(jì)等帶來(lái)難度。Alpha 對(duì)日定向裝置長(zhǎng)期在軌工作過(guò)程中的性能退化對(duì)長(zhǎng)期保持低速穩(wěn)定驅(qū)動(dòng)控制帶來(lái)挑戰(zhàn)。因此，突破空間站大柔性太陽(yáng)電池翼高穩(wěn)定驅(qū)動(dòng)控制技術(shù)，對(duì)于保障空間站能源安全和平穩(wěn)運(yùn)行至關(guān)重要。

以國(guó)際空間站（International Space Station，ISS）為例，與Alpha 對(duì)日定向裝置輸出端連接的末端桁架上配置4 個(gè)太陽(yáng) 電池翼，如圖3 所示［2］，翼展超過(guò)70 m，面積552 m2，1 階頻率0.06 Hz，轉(zhuǎn)動(dòng)慣量5.4×105kg·m2。為實(shí)現(xiàn)柔性負(fù)載穩(wěn)定驅(qū)動(dòng)，采用帶4 階Butterworth 低通濾波器的速度、位置雙環(huán)系統(tǒng)，控制參數(shù)的確定通過(guò)構(gòu)造與兩環(huán)控制帶寬相關(guān)的受限優(yōu)化函數(shù)（Constrained Optimization Scheme）進(jìn) 行綜合［3-4］。

圖3 國(guó)際空間站Alpha 對(duì)日定向裝置與太陽(yáng)電池翼Fig.3 Solar Alpha rotary joint（SARJ）and solar cell wings of ISS

隨著現(xiàn)代航天器柔性附件（太陽(yáng)翼、天線、桁架等）尺寸越來(lái)越大，對(duì)于大柔性附件高穩(wěn)定驅(qū)動(dòng)控制系統(tǒng)的研究，得到了學(xué)者們重視［5］。除國(guó)際空間站上的應(yīng)用外，BONG 等［6］在帶有大柔性太陽(yáng)翼的哈勃望遠(yuǎn)鏡（Hubble Space Telescope，HST）指向控制中，提出采用帶陷阱濾波的魯棒控制方案，并與經(jīng)典PID控制進(jìn)行對(duì)比，認(rèn)為該方法在確保平臺(tái)速度穩(wěn)定度和撓性振動(dòng)抑制上有較好的效果。YOSHIRO 等［7］在日本工程試驗(yàn)衛(wèi)星八號(hào)（Japan Engineering Test Satellite Ⅷ，ETS-Ⅷ）上，為解決柔性天線大角度跟蹤指向控制，采用線性插值規(guī)劃控制（Linearly Interpolated Gain Scheduling Controller）進(jìn)行運(yùn)動(dòng)擬合，也得到了較好的速度穩(wěn)定度效果。針對(duì)線性控制器對(duì)大柔性體非線性系統(tǒng)調(diào)節(jié)能力不足的問(wèn)題，學(xué)者們考慮將輸入成型控制、魯棒控制［8］、變結(jié)構(gòu)控制、自適應(yīng)控制等應(yīng)用于大柔性體伺服控制，以滿足系統(tǒng)動(dòng)、靜態(tài)性能指標(biāo)［9］。

本文以空間站大柔性太陽(yáng)電池翼高穩(wěn)定對(duì)日跟蹤控制為研究對(duì)象，推導(dǎo)柔性太陽(yáng)電池翼動(dòng)力學(xué)模型和開展有限元分析，建立Alpha 對(duì)日定向裝置機(jī)電系統(tǒng)模型和非線性傳動(dòng)系統(tǒng)模型；針對(duì)柔性負(fù)載低頻和傳動(dòng)系統(tǒng)非線性對(duì)控制性能的影響，設(shè)計(jì)帶變速規(guī)劃和柔性振動(dòng)抑制的速度、位置雙環(huán)控制系統(tǒng)；并通過(guò)仿真校驗(yàn)，證明控制方案的有效性。

1 控制系統(tǒng)建模

1.1 柔性太陽(yáng)電池翼動(dòng)力學(xué)建模

參考我國(guó)空間站Alpha 對(duì)日定向裝置和太陽(yáng)電池翼的相對(duì)布局關(guān)系，采用混合坐標(biāo)法和有限元法對(duì)柔性負(fù)載（末端桁架和兩側(cè)太陽(yáng)電池翼的組合體）進(jìn)行離散化，并用Lagrange 方程建立柔性動(dòng)力學(xué)模型，得到太陽(yáng)電池翼驅(qū)動(dòng)力矩與轉(zhuǎn)速的狀態(tài)方程［10］。

建模中引入空間站本體坐標(biāo)系，末端桁架固連隨動(dòng)坐標(biāo)系{α}。{α}原點(diǎn)Oα在Alpha 對(duì)日定向裝置與末端桁架連接法蘭面幾何中心，rb為空間站本體坐標(biāo)系原點(diǎn)Ob到末端桁架坐標(biāo)系{α}原點(diǎn)Oα的矢徑，rα，j為桁架和太陽(yáng)電池翼組成的柔性體上任意點(diǎn)j相對(duì)于坐標(biāo)系{α}的矢徑，δα，j為j點(diǎn)的變形位移，Rb，j為j點(diǎn)相對(duì)于坐標(biāo)系的位移矢量，具體描述如圖4 所示。

圖4 單自由度轉(zhuǎn)動(dòng)柔性負(fù)載坐標(biāo)系定義Fig.4 Definition of the coordinate system of single-degree-of-freedom rotating flexible body

矢徑Rb，j及其導(dǎo)數(shù)的表達(dá)式為

式中：ωα=(ωα，x，ωα，y，ωα，z)為{α}系相對(duì)于系的角速度矢量在系中的描述，其中，ωα，y為Alpha對(duì)日定向裝置轉(zhuǎn)動(dòng)角速度。

假設(shè)Cb2α為由系到{α}系的變換矩陣：

式中：θα為Alpha 對(duì)日定向裝置轉(zhuǎn)過(guò)的角位移（yα軸方向）。

柔性負(fù)載上任意點(diǎn)j的速度公式（2），轉(zhuǎn)換到隨動(dòng)系{α}系下描述如下：

假設(shè)向量δα表示柔性負(fù)載上所有節(jié)點(diǎn)的線彈性變形，可用柔性負(fù)載模態(tài)分析得到的基向量的線性組合表示：

式中：Φα=(Φα，1，Φα，2，…，Φα，j，…，Φα，q)表示柔性負(fù)載的振型基?向量陣；ηα=(ηα，1，ηα，2，…，ηα，j，…，ηα，q)T表示?柔性負(fù)載的模態(tài)坐標(biāo)，其中，q為柔性負(fù)載的總節(jié)點(diǎn)數(shù)。

經(jīng)運(yùn)動(dòng)學(xué)描述，可得系統(tǒng)動(dòng)能T和勢(shì)能V表達(dá)式為式中：mα，j為柔性負(fù)載上的j點(diǎn)集中質(zhì)量；Κα為剛度矩陣為柔性負(fù)載的固有頻率矩陣（n×n階）。

Jα為柔性負(fù)載相對(duì)坐標(biāo)系{α}原點(diǎn)Oα的轉(zhuǎn)動(dòng)慣量矩陣：

Fα為柔性負(fù)載振動(dòng)對(duì)其轉(zhuǎn)動(dòng)的耦合系數(shù)矩陣：

系統(tǒng)Lagrange 函數(shù)為

將式（6）、式（7）代入式（10），建立第2 類Lagrange 方程如下：

描述系統(tǒng)繞Alpha 對(duì)日定向裝置轉(zhuǎn)動(dòng)和彈性振動(dòng)的系統(tǒng)動(dòng)力學(xué)方程為

式中：Tα為Alpha 對(duì)日定向裝置驅(qū)動(dòng)力矩；ζα=diag(ζα，1，ζα，2，…，ζα，n)為柔性負(fù)載各階模態(tài)的阻尼比對(duì)角陣；Ωα為柔性負(fù)載圓頻率對(duì)角陣Ωα=diag(ωα，1，ωα，2，ωα，3，…，ωα，n)，n為模態(tài)階數(shù)，有

將式（12）、式（13）用狀態(tài)方程描述，輸入為Alpha 對(duì)日定向裝置驅(qū)動(dòng)力矩Tα，y，輸出為柔性負(fù)載轉(zhuǎn)動(dòng)角加速度αα，（y可由αα，y積分得到角速度ωα，y）。狀態(tài)方程的狀態(tài)量為柔性負(fù)載的模態(tài)坐標(biāo)及其變化率，可通過(guò)對(duì)模態(tài)坐標(biāo)動(dòng)態(tài)變化評(píng)估各階頻率在振動(dòng)響應(yīng)中所占權(quán)重和演化規(guī)律，并結(jié)合模態(tài)基向量得到每1 個(gè)離散點(diǎn)的位移量。

狀態(tài)方程形式如下：

式中：A、B、C分別為2n×2n、2n×1、1×2n的系數(shù)矩陣；D為實(shí)數(shù)系數(shù)；n為柔性負(fù)載模態(tài)分析頻率截?cái)嗪蟮谋Ａ綦A數(shù)。則有

由式（15）的第1 式調(diào)整后為

將式（16）代入式（15）第2 式，可得狀態(tài)方程（14）的第1 式。推導(dǎo)過(guò)程如下：

又式（15）的第2 式代入第1 式，調(diào)整后可得狀態(tài)方程（14）的第2 式如下：

得Alpha 對(duì)日定向裝置驅(qū)動(dòng)力矩與柔性負(fù)載轉(zhuǎn)動(dòng)角加速度（輸出量積分得到角速度）的狀態(tài)方程為

1.2 太陽(yáng)電池翼模態(tài)分析與截?cái)?/h3>

柔性體狀態(tài)方程如式（19）所示，為得到方程系數(shù)，需要對(duì)太陽(yáng)電池翼進(jìn)行模態(tài)分析，獲取其各階頻率和振型。如空間站太陽(yáng)電池翼這樣尺寸大、基頻低、模態(tài)密集的結(jié)構(gòu)，直接采用模態(tài)截?cái)嗳匀粚?dǎo)致模型階數(shù)較高。先根據(jù)伺服控制帶寬進(jìn)行初次截?cái)?，再依?jù)慣性完備性準(zhǔn)則和模態(tài)有效質(zhì)量進(jìn)行模態(tài)篩選［11］。

根據(jù)末端桁架和太陽(yáng)電池翼的設(shè)計(jì)方案，用MSC.Nastran 建立非線性有限元模型，如圖5 所示。設(shè)轉(zhuǎn)軸正方向?yàn)閥α、原點(diǎn)為Alpha，對(duì)日定向裝置與末端桁架連接法蘭的幾何中心。有限元模型節(jié)點(diǎn)41 976個(gè)，各類單元41 465個(gè)，各類多點(diǎn)約束連接237個(gè)。采用非線性動(dòng)力學(xué)求解器SOL 106對(duì)柔性負(fù)載進(jìn)行模態(tài)分析，其中前3階主模態(tài)如圖6～圖8所示。

圖5 太陽(yáng)電池翼有限元模型Fig.5 Finite element model of solar cell wings

圖6 太陽(yáng)電池翼第1 階主模態(tài)兩對(duì)稱振型（0.066 9 Hz）Fig.6 The first main mode of solar cell wings（0.066 9 Hz）

圖7 太陽(yáng)電池翼第2 階主模態(tài)兩對(duì)稱振型（0.067 2 Hz）Fig.7 The second main mode of solar cell wings（0.067 2 Hz）

圖8 太陽(yáng)電池翼第3 階主模態(tài)兩對(duì)稱振型（0.123 0 Hz）Fig.8 The third main mode of solar cell wings（0.123 0 Hz）

對(duì)1.0 Hz和2.0 Hz截止頻率下的模態(tài)階數(shù)和模態(tài)有效質(zhì)量進(jìn)行匯總，見表1。

表1 不同截止頻率下的模態(tài)階數(shù)和模態(tài)有效質(zhì)量Tab.1 Modal orders and effective mass at different cutoff frequencies

由于柔性負(fù)載在Alpha對(duì)日定向裝置驅(qū)動(dòng)下實(shí)現(xiàn)繞坐標(biāo)系{α}的yα軸轉(zhuǎn)動(dòng)，在驅(qū)動(dòng)力矩Tα，y作用下，柔性負(fù)載彈性變形主要體現(xiàn)在zα軸方向的平動(dòng)和繞yα軸的轉(zhuǎn)動(dòng)。而模態(tài)截?cái)嘀?.0 Hz或1.0 Hz，其RY向模態(tài)有效質(zhì)量已經(jīng)大于98%，保留的動(dòng)力學(xué)特性可為工程接受。頻率2.0 Hz內(nèi)的模態(tài)有638階，需要根據(jù)模態(tài)有效質(zhì)量進(jìn)一步完備性篩選，減少系統(tǒng)階數(shù)。

1.3 Alpha 對(duì)日定向裝置機(jī)電系統(tǒng)建模

Alpha 對(duì)日定向裝置采用永磁同步電機(jī)配合高精度角度傳感器，實(shí)現(xiàn)寬調(diào)速范圍的多級(jí)閉環(huán)驅(qū)動(dòng)控制。取永磁體基波磁場(chǎng)軸線（轉(zhuǎn)子N極）為d軸，q軸順著旋轉(zhuǎn)方向超前d軸90°電角度。電機(jī)模型遵循如下假設(shè)：1）忽略磁路飽和、磁滯和渦流影響，磁路線性，無(wú)阻尼；2）定子繞組三相對(duì)稱，各繞組軸線在空間上互差120°；3）電機(jī)電勢(shì)正弦，定子電流在氣隙中只產(chǎn)生正弦分布磁勢(shì)，忽略磁場(chǎng)高次諧波［12］。

永磁同步電機(jī)動(dòng)力學(xué)模型為

式中：uq、ud為交軸、直軸電壓；iq、id為交軸、直軸電流；Lq、Ld為電機(jī)直軸、交軸同步電感，Lq=Ld；ψf為轉(zhuǎn)子永磁體的勵(lì)磁磁鏈；ωi為轉(zhuǎn)子電角速度（轉(zhuǎn)子機(jī)械角速度ωm=ωi/Pn）；Rs為定子電阻；Pn為定子繞組極對(duì)數(shù)；Jm為電機(jī)機(jī)械轉(zhuǎn)動(dòng)慣量；Te為電磁轉(zhuǎn)矩；Bm為黏性摩擦系數(shù)；Tm為Alpha 對(duì)日定向裝置電機(jī)端輸出力矩，通過(guò)傳動(dòng)系統(tǒng)放大，并克服傳動(dòng)摩擦力矩TD后，為Alpha對(duì)日定向裝置驅(qū)動(dòng)力矩Tα，對(duì)應(yīng)式（19）中的u。

1.4 傳動(dòng)系統(tǒng)建模

1.4.1 間隙模型

Alpha 對(duì)日定向裝置的多級(jí)減速器會(huì)帶來(lái)設(shè)計(jì)中不可忽略的傳動(dòng)間隙，在控制系統(tǒng)設(shè)計(jì)中需克服間隙引起“死區(qū)”非線性。Alpha 對(duì)日定向裝置減速器裝配間隙等效模型如圖9 所示。

圖9 傳動(dòng)間隙等效模型Fig.9 Equivalent model of transmission gap

圖中：θretarder為減速器折算到柔性負(fù)載端的間隙；Tm、ωm為Alpha 對(duì)日定向裝置電機(jī)端力矩和轉(zhuǎn)動(dòng)角速度；TL、ωf為Alpha 對(duì)日定向裝置輸出端驅(qū)動(dòng)力矩和轉(zhuǎn)動(dòng)角速度。

當(dāng)傳動(dòng)軸彈性形變角度θmf滿足切換條件|θmf| ≤θretarder/2 時(shí)，Alpha 對(duì)日定向裝置輸出端懸空：

式中：輸出驅(qū)動(dòng)力矩TL為零；電機(jī)端負(fù)載力矩Tm為傳動(dòng)系統(tǒng)摩擦力矩TD，反映傳動(dòng)摩擦特性；i為減速比。

當(dāng)間隙消除，即|θmf|＞θretarder/2 時(shí)，

式中：Alpha 對(duì)日定向裝置驅(qū)動(dòng)力矩TL為凈扭轉(zhuǎn)形變力矩。

1.4.2 摩擦模型

傳動(dòng)系統(tǒng)中的摩擦力矩會(huì)降低傳動(dòng)效率，并且在低速轉(zhuǎn)動(dòng)時(shí)由于動(dòng)、靜摩擦力矩的切換導(dǎo)致“爬行現(xiàn)象”。但摩擦力矩的存在也可以有效減緩柔性負(fù)載振動(dòng)對(duì)Alpha 對(duì)日定向裝置電機(jī)的反作用力矩，實(shí)現(xiàn)振動(dòng)隔離，保護(hù)驅(qū)動(dòng)電機(jī)。

結(jié)合試驗(yàn)數(shù)據(jù)，可擬合出Alpha 對(duì)日定向裝置的傳動(dòng)摩擦特性。采用動(dòng)態(tài)LuGre 模型擬合，模擬靜、動(dòng)摩擦間的切換。摩擦力矩TD數(shù)學(xué)表達(dá)式如下：

式中：ωf=ωα，y為Alpha 對(duì)日定向裝置轉(zhuǎn)動(dòng)角速度；z為平均變形量；σ0為剛性系數(shù)；σ1為滑動(dòng)阻尼系數(shù)；σ2為黏性摩擦系數(shù)；TD_dyn為庫(kù)倫摩擦力矩；TD_sta為最大靜摩擦力矩；θ˙s為Stribeck 特征速度［13］。

假設(shè)傳動(dòng)系統(tǒng)靜摩擦力矩為80 N·m，動(dòng)摩擦力矩為60 N·m，構(gòu)造摩擦力矩TD是轉(zhuǎn)動(dòng)角速度ωf的函數(shù)，LuGre 摩擦曲線如圖10 所示。

圖10 LuGre 摩擦模型曲線Fig.10 LuGre friction model curve

2 穩(wěn)定跟蹤控制方案設(shè)計(jì)

為提高驅(qū)動(dòng)穩(wěn)定性，采用帶運(yùn)動(dòng)規(guī)劃和振動(dòng)抑制的速度、位置雙閉環(huán)伺服方案?？刂葡到y(tǒng)框圖如圖11 所示。

圖11 對(duì)日定向裝置驅(qū)動(dòng)控制框圖Fig.11 Drive control block diagram of the SARJ

2.1 變速規(guī)劃方案

提高啟動(dòng)和制動(dòng)等變速過(guò)程的平穩(wěn)性，降低柔性體低頻擾動(dòng)載荷，采用Heaviside 五次樣條變速規(guī)劃策略，實(shí)現(xiàn)1、2 階導(dǎo)數(shù)連續(xù)，且角度和角速度規(guī)劃匹配。在調(diào)速過(guò)程中，通過(guò)Heaviside 變速規(guī)劃，也可克服傳動(dòng)間隙影響，表達(dá)式如下：

2.2 位置環(huán)和速度環(huán)方案

位置環(huán)、速度環(huán)調(diào)節(jié)器為消除靜差，提高跟蹤精度，引入積分環(huán)節(jié)。但傳動(dòng)間隙θretarder會(huì)導(dǎo)致末端反饋測(cè)量偏差較大，造成速度環(huán)調(diào)節(jié)器積分累積超過(guò)限幅閾值，引發(fā)系統(tǒng)振蕩。因此，在位置環(huán)和速度環(huán)都引入積分分離控制器。

當(dāng)被控量與規(guī)劃值偏差較大時(shí)，取消積分作用，以免積分導(dǎo)致系統(tǒng)穩(wěn)定性降低；當(dāng)被控量接近規(guī)劃值時(shí)，引入積分控制，消除靜差，提高跟蹤精度［14］。

積分分離控制算法表示為

式中：Δt為采樣時(shí)間；β為積分項(xiàng)的開關(guān)系數(shù)，

2.3 振動(dòng)抑制方案

太陽(yáng)電池翼受到外界擾動(dòng)時(shí)，會(huì)激發(fā)其彈性低頻振動(dòng)，且持續(xù)時(shí)間較長(zhǎng)。為快速抑制柔性太陽(yáng)電池翼受外界擾動(dòng)載荷引起的振動(dòng)，在速度調(diào)節(jié)器前，設(shè)計(jì)特定頻率的陷阱濾波器（Notch Filter）：

通過(guò)設(shè)置陷阱濾波器參數(shù)ωz、ωp、ζp、ζz，使濾波器作用區(qū)域限制在太陽(yáng)電池翼基頻附近。參考模態(tài)分析結(jié)果，柔性負(fù)載彈性主要體現(xiàn)在1 階模態(tài)上（0.066 9 Hz），取ωz=0.42，ωp=0.377，ζp=0.70，ζz=0.02。陷阱濾波器傳遞函數(shù)Bode 圖如圖12所示。

圖12 陷阱濾波器幅頻、相頻曲線Fig.12 Phase frequency and amplitude frequency curves of the notch filter

3 仿真校驗(yàn)

空間站運(yùn)行在380～410 km 之間的軌道高度，為確保柔性太陽(yáng)翼陣面與太陽(yáng)光入射方向垂直，常見工況為：Alpha 對(duì)日定向裝置加電啟動(dòng)，從停轉(zhuǎn)狀態(tài)加速到0.300（°）/s，開始對(duì)太陽(yáng)進(jìn)行搜索捕獲，捕獲太陽(yáng)后Alpha 對(duì)日定向裝置轉(zhuǎn)入空間站軌道角速度0.065（°）/s 持續(xù)對(duì)日跟蹤模式。

假設(shè)，柔性負(fù)載轉(zhuǎn)動(dòng)慣量為

主模態(tài)頻率fa=（0.066 9，0.067 2，0.123 0，0.180 0，0.181 7，…）Hz。

柔性負(fù)載振動(dòng)對(duì)其轉(zhuǎn)動(dòng)的耦合系數(shù)矩陣為

得到反映柔性負(fù)載傳遞函數(shù)Bode 圖（輸入驅(qū)動(dòng)、輸出角速度）如圖13 所示。圖中，藍(lán)色實(shí)線為其幅頻特性、相頻特性曲線，紅色虛線為簡(jiǎn)化負(fù)載為剛體時(shí)的動(dòng)態(tài)特性，兩者差異明顯。由于給出的是Alpha 對(duì)日定向裝置轉(zhuǎn)動(dòng)方向的特性，由Bode 圖可知0.067、0.123、0.180、0.882、0.956 Hz 等頻率對(duì)其動(dòng)力學(xué)振動(dòng)特性起主導(dǎo)作用。

圖13 柔性負(fù)載動(dòng)力學(xué)特性Bode 圖Fig.13 Bode diagram of flexible body transfer function

設(shè)Alpha 對(duì)日定向裝置永磁同步電機(jī)、傳動(dòng)系統(tǒng)的參數(shù)見表2［15］。

表2 Alpha 對(duì)日定向裝置參數(shù)Tab.2 Input parameters of the SARJ

設(shè)計(jì)控制參數(shù)，位置環(huán)控制帶寬0.027 8 Hz；速度環(huán)控制帶寬0.039 5 Hz。仿真中柔性負(fù)載的啟動(dòng)規(guī)劃時(shí)間180 s，分析結(jié)果如圖14～圖19所示。

圖14 伺服角速度與規(guī)劃指令角速度Fig.14 Servo angular velocity and planning angular velocity

圖15 伺服角速度與規(guī)劃指令角速度偏差Fig.15 Deviation between the servo angular velocity and the planning angular velocity

圖16 伺服角度與規(guī)劃指令角度Fig.16 Servo angle and planning angle of the SARJ

圖17 伺服角度與規(guī)劃指令角度偏差Fig.17 Deviation between the servo angle and the planning angle

圖18 Alpha 對(duì)日定向裝置輸出力矩Fig.18 Output torque of the SARJ

圖19 前3 階主頻率模態(tài)坐標(biāo)Fig.19 Modal coordinates of the first third order frequencies

從圖14～圖17 可知，在啟動(dòng)規(guī)劃過(guò)程中受傳動(dòng)間隙和傳動(dòng)摩擦影響，存在速度跟蹤偏差，180 s 規(guī)劃完成時(shí)已經(jīng)實(shí)現(xiàn)穩(wěn)定跟蹤。跟蹤精度優(yōu)于0.3°，速度偏 差±0.005（° ）/s，穩(wěn)定度優(yōu)于 7%＠0.065（°）/s。

從圖18 和圖19 可知，Alpha 對(duì)日定向裝置的啟動(dòng)和變速過(guò)程中驅(qū)動(dòng)力矩不大于30 N·m。穩(wěn)速運(yùn)行期間，驅(qū)動(dòng)力矩不大于5 N·m。驅(qū)動(dòng)過(guò)程中激發(fā)柔性負(fù)載第1、3 階模態(tài)，當(dāng)驅(qū)動(dòng)處于穩(wěn)速轉(zhuǎn)動(dòng)后，低階柔性振動(dòng)也得到抑制。

4 結(jié)束語(yǔ)

本文研究空間站柔性太陽(yáng)電池翼高穩(wěn)定對(duì)日跟蹤控制問(wèn)題，推導(dǎo)了柔性太陽(yáng)電池翼驅(qū)動(dòng)力矩與轉(zhuǎn)速狀態(tài)方程，建立考慮傳動(dòng)間隙和傳動(dòng)摩擦的Alpha 對(duì)日定向裝置模型，提出一種帶運(yùn)動(dòng)規(guī)劃和振動(dòng)抑制的速度、位置雙閉環(huán)伺服控制方法。文中給出的仿真算例說(shuō)明了該方法在實(shí)現(xiàn)柔性負(fù)載穩(wěn)定驅(qū)動(dòng)的有效性。文中所述方案控制性能好于國(guó)際空間站Alpha 對(duì)日跟蹤控制設(shè)計(jì)指標(biāo)：穩(wěn)態(tài)跟蹤偏差不大于 2°，速度穩(wěn) 定度優(yōu)于10%＠0.065（°）/s［3-4］。