首页 / 专利库 / 物理 / 临界状态 / 一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法

一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法

阅读:451发布:2021-07-21

专利汇可以提供一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法专利检索,专利查询,专利分析的服务。并且本 发明 提供一种基于I-V曲线与 能量 平衡的 小卫星 太阳阵寿命预测方法,针对当前小卫星太阳阵寿命预测存在的损伤规律空间环境影响因素较少,及较大程度上局限于 单体 电池 的寿命预测而无法解决太阳阵整体寿命预测等问题,本发明基于空间环境模拟试验,在考虑日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹 角 、 温度 、 太阳 辐射 等重要空间环境因素影响下太阳阵整体寿命 预测模型 ,进而部分解决考虑太阳阵重要寿命影响因素的寿命预测及适用于不同批次、不同型号太阳阵整体寿命预测的通用性问题。所述I-V曲线为太阳单体电池的 电流 与 电压 曲线。本发明以试验为 基础 ,具有较好的模型通用性以及更强的工程实用性。,下面是一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法专利的具体信息内容。

1.一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法,其特征在于:该方法通过如下步骤实现:
步骤一、日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹确定;
根据轨道高度、降交点地方时和预测起始时间,计算每天的日地距离因子、轨道周期Te、轨道地影时间、每轨太阳光线与太阳阵法线的夹角的变化规律,得到随时间变化的定量数据,用于后续的太阳阵I-V曲线和能量平衡分析;
步骤二、太阳阵I-V曲线模型确定;
根据太阳阵特性构建太阳阵计算模型,同时考虑太阳入射角、辐照衰减、日地因子、损失因子因素的影响,计算不同季节、不同轨道条件和不同工况下太阳阵输出电压、输出电流,以表征太阳阵输出功率实时及长期变化情况;
以标准状态的I-V曲线特征点为参数,考虑多种环境因素对太阳阵寿命的影响,计算太阳阵的输出特性;利用公式(Equ.1)太阳阵I-V曲线的计算机解析模型,得到不同条件下的太阳阵的I-V特性曲线;该模型在光照强度小于2个太阳常数时,有很高的精确性;太阳同步轨道小卫星的光照情况满足这一条件:
式中:
I——太阳阵输出电流,单位为A;
Isc'——太阳阵短路电流,典型参数或实测值,单位为A;
C1——公式系数1;
V——太阳阵输出电压,单位为V;
C2——公式系数2;
Vov'——太阳阵开路电压,典型参数或实测值,单位为V;
Imp'——太阳阵最佳工作点输出电流,典型参数或实测值,单位为A;
Vmp'——太阳阵最佳工作点输出电压,典型参数或实测值,单位为V;
太阳阵开路电压和最佳工作点输出电压计算模型如下:
式中:
Vov——单体太阳电池开路电压,单位为V;
Vmp——单体太阳电池最佳工作点电压,单位为V;
βVBOL——单体太阳电池寿命初期电压温度系数,单位为V/℃;
KVRAD——太阳阵开路电压辐照衰降因子;
T——太阳阵温度,单位为℃;
太阳阵短路电流和最佳工作点电流计算模型如下:
ISC——单体太阳电池短路电流,单位为A;
Imp——单体太阳电池最佳工作点电流,单位为A;
αI—单体太阳电池电流温度系数,单位为A/℃;
θ(t)—— 一圈轨道内太阳光线与太阳阵法线方向的夹角,单位为度;
T——太阳阵温度,单位为℃;
KIRAD________太阳阵短路电流辐照衰降因子;
Frd——日地距离因子;
利用“太阳阵开路电压及短路电流辐照衰降因子计算模型”预测LEO轨道辐射环境对卫星太阳电池输出参数衰减的影响,在该模型中Isc即为KIRAD,Vov即为KVRAD;
A.模型输入参数定义如下:
电池类型:单结GaAs太阳电池;石英玻璃盖片厚度:120μm;轨道高度:
300km~3000km;倾角:只针对99°;时间单位:月;
B.模型输出参数定义如下:
最大输出功率Pmax、短路电流Isc、开路电压Vov,其输出形式:给出Pmax、Isc和Vov经过m个月后,Pmax、Isc和Vov为初始值的百分比,即给出Pmax、Isc和Vov关于时间month的函数;
以下为该太阳阵开路电压及短路电流辐照衰降因子的计算模型:
不同轨道高度位移损伤剂量计算如下:x为轨道高度,month为在轨月数,y为计算得到的位移损伤剂量;
当300km<=x<=600km时,计算公式为:
2 3 4 5
y=14(A0+A1·x+A2·x+A3·x+A4·x+A5·x)·month (Equ.4)
其中,A0=-5.72637E6,A1=69074.68933,A2=-329.19032,
A3=0.77634,A4=-9.13546E-4,A5=4.49106E-7
当600km2 3 4
y=14(A0+A1·x+A2·x+A3·x+A4·x)·month (Equ.5)
其 中,A0=-5.80893E7,A1=321272.30685,A2=-663.23216,A3=0.59526,A4=-1.77968E-4
当1000km2 3 4 5
y=14(A0+A1·x+A2·x+A3·x+A4·x+A5·x)·month (Equ.6)
其中,A0=5.01219E8,A1=-1.76649E6,A2=2453.54778,
A3=-1.65135,A4=5.32602E-4,A5=-5.18233E-8
GaAs/Ge太阳电池的Pmax、Isc和Voc的计算模型为:
最大输出功率衰减,即Pmax的计算模型:
Pmax=1.0-C×log10(1+(y/Dx)) (Equ.7)其中,C=0.242,Dx=3.47e9,y为计算得到的位移损伤剂量;
短路电流衰减,即Isc的计算模型:
KIRAD=Isc=1.0-C×log10(1+(y/Dx)) (Equ.8)其中,C=0.213,Dx=8.3e19
开路电压衰减,即Voc的计算模型:
KVRAD=Vov=1.0-C×log10(1+(y/Dx)) (Equ.9)其中,C=0.07,Dx=1.8e9
步骤三、太阳阵能量平衡计算模型确定;
在进行能量平衡计算时,根据在轨数据或地面提供的数据对能量平衡的临界状态进行实时监控,如果太阳阵提供能量的多余电量Qresidual(c)由正值转变为小于等于零的任意值,则表明太阳阵已处于严重损伤状态,且Qresidual(c)的计算式为:
其中:
Qresidual(c)——在轨第c圈太阳阵可提供的多余电量,单位为C;
te——阴影期时间,单位为s;
ISA(c)——在轨第c圈方阵电流箝位点电流值,单位为A;
Iload_mean(c)——光照期负载电流Iload(A),其为在轨第c圈负载电流每周期的平均值,单位为A;
Id(t)——阴影期,蓄电池放电电流,单位为A;
根据指定时期太阳阵I-V曲线的方程,在给出相应的光照区母线电压Vs_bus及太阳阵隔离二极管和供电线缆压降之和Vs_dioline时,得到该指定时期I-V曲线上太阳阵工作电压箝位点Vop1处的电流值Is_op1;由能量平衡计算可知,该指定时期的太阳阵提供能量的多余电量Qresidual(c)可表示为:
式中:
Is_opl—指点时期I-V曲线上太阳阵工作电压箝位点Vop1处的电流值,单位为A;
Is_load_mean——指定时期光照区负载电流/在轨所有负载电流数据的平均值,单位为A;
Id(t)——指定周期蓄电池在阴影区的放电电流值,单位为A;
太阳阵工作电压点输出功率计算模型如下:
PSA(t)=Vbus(t)Iop1(t), ISA(t)=Iop1(t)
进一步可得:
Ps_op1(c)=Vs_bus(c)Is_op1(c)
Vs_op1(c)=Vs_bus(c)+Vs_dioline (Equ.12)
式中:
t――第c圈轨道周期内时刻,0PSA——太阳阵输出功率,单位为W;
Ps_OP1(c)——第c圈太阳阵输出为箝位点Vs_op1时输出功率,单位为W;
Is_opI(c)——第c圈I-V曲线上太阳阵工作电压箝位点Vs_op1时电流值,单位为A;
Vs_bus——光照区母线电压,单位为V;
Vs_dioline——太阳阵隔离二极管和供电线缆压降之和,单位为V;
对于在轨电源系统,由于电源控制器使得光照期母线电压始终保持为定值,因而,可以认为:在电源控制器正常工作的前提下,光照期母线电压始终不变;同时,在得到Vs_dioline值后,可够通过上述建立的I-V曲线得到指定时期箝位点电流值Is_op1(c),以用于能量平衡计算;
如果系统预测得到某指定时期Qresidual(c)=0,则说明此时太阳阵已到寿;
其中,进行蓄电池组放电电流计算时,蓄电池组的放电电流取决于蓄电池组的放电功率、放电调节器效率、蓄电池组供电线路损耗因子、电池组电压因素;
阴影区,蓄电池组放电电流为:
式中:
t――第c圈轨道周期内时刻,0Iload(t)――在轨负载电流需求随时间变化的函数;
ηBDR――放电调节器效率;
nline――蓄电池组供电线路损耗因子;
Vbat(t)――蓄电池组放电电压;基于在轨第c圈轨道蓄电池放电初压和放电终压值,可以近似地认为Vbat(t)由放电初压至放电终压线性变化;
Vbus:放电时母线电压,单位为V;
其中:阴影区太阳阵ISA(t)电流为零,Vbat(t)是一个被积函数,该函数是由蓄电池在轨放电初压和放电终压确定的线性函数。
2.根据权利要求1所述的一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法,其特征在于:所述的太阳阵温度通过太阳阵温度模型计算如下:
太阳阵温度随卫星进出影状态的变化而变化,地影区,太阳阵温度逐渐下降,直至降至出影前的最低温度;光照期,太阳阵温度从出影后迅速上升,直至达到光照期的温度平衡点,此后温度保持不变直至卫星进入下一轨道圈的地影期,周而复始;
太阳阵温度变化的简化模型如下:
在地影期内,太阳阵温度从光照期的最高平衡温度线性下降至地影期最低温度,出影后,太阳阵温度在8分钟内从地影期最低温度上升至60℃,在20分钟内从60℃上升至光照期最高平衡温度,直至下次进影;
光照期的最高平衡温度、地影期最低温度的默认值分别为:光照期的最高平衡温度Tx;地影期最低温度TSAE。

说明书全文

一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方

技术领域

[0001] 本发明属于小卫星太阳阵预测技术领域,具体涉及一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法

背景技术

[0002] 小卫星应用于除载人航天外的所有空间领域,包括在遥感、通信、导航、技术验证、空间科学等发挥了重要的作用。电源系统是卫星的关键服务系统,负责产生、存储和为卫星整个寿命期间为整星用电负载提供稳定的不间断的能源,其供电能、电源品质直接影响卫星的工作状态、可靠性与使用寿命。对电源系统控制的性能、可靠性有很高的要求。
[0003] 由于不同卫星运行轨道有其明显的环境特点,导致不同轨道运行的卫星系统性能退化规律和特征有所区别。小卫星的运行轨道大部分为中低轨道,受复杂空间环境影响,如:太阳压力、电离层、带电粒子等,太阳入射在一年中的变化范围大,小卫星进出影频繁,高低温交替变换等,太阳阵间歇性工作次数多,工作时间短,工作电流大,高低温冲击剧烈。这些因素不可避免的将对卫星的性能与寿命产生影响。最终在内部和外部因素的综合作用下,导致系统失效。图1给出了,由于小卫星各子系统失效而导致小卫星整体失效所占的比例情况。由此可知,电源系统失效占据大部分比例,且太阳阵失效是电源系统失效的最主要因素。
[0004] 然而,由于小卫星属于不可修复产品,且受质量和尺寸的限制,不能采取冗余部件的方法以提高小卫星太阳阵的可靠性,使得小卫星太阳阵寿命预测对于小卫星的设计、生产、使用起着重要指导作用。
[0005] 目前应用于空间环境条件下小卫星太阳阵寿命预测的方法可归纳为:(一)、紫外加速寿命试验法:建立模拟空间环境下的紫外加速寿命试验装置,对太阳电池进行紫外加速寿命试验技术研究,获得太阳电池开路电压随着紫外辐照时间的变化数据。通过试验数据处理,获得太阳电池开路电压随着等效紫外辐照时间的衰减规律,采用加严判据理论,研究紫外辐射对太阳电池的损伤规律,预测紫外辐照环境下的电池寿命;(二)、热应变极值法:以试验测试为基础分析多次热循环后太阳电池板的热应变演变规律,提出以热应变极大值(或残余热应变)作为多层胶接结构的损伤参量,建立预测太阳单体电池结构寿命的数学模型等方法。其中,紫外加速寿命试验法仅仅考虑了紫外这一空间环境条件,且较其它寿命影响因素,紫外的影响相对较弱,单一开展紫外辐射对太阳阵的损伤规律研究不能很好的揭露太阳阵在复杂空间环境条件下的损伤规律;热应变极值法从通过试验建立了单体电池的数据模型,对于经试验测试,且仅考虑单体电池在温度影响下的寿命预测是有效的,对于不同批次、不同型号等太阳单体需要重新建立寿命预测模型,且无法解决太阳阵整体寿命预测问题。

发明内容

[0006] 针对当前小卫星太阳阵寿命预测存在的问题,本发明提出一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法,在空间环境模拟试验的基础上,构建在考虑日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹角、温度、太阳辐射等重要空间环境因素影响下的太阳阵整体寿命预测模型,进而部分解决考虑太阳阵重要寿命影响因素的寿命预测及太阳阵整体寿命预测的通用性问题。所述I-V曲线为太阳单体电池的电流与电压曲线。
[0007] 本发明采用的技术方案为:一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法,该方法通过如下步骤实现:
[0008] 步骤一、日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹角确定;
[0009] 根据轨道高度、降交点地方时和预测起始时间,计算每天的日地距离因子、轨道周期Te、轨道地影时间、每轨太阳光线与太阳阵法线的夹角的变化规律,得到随时间变化的定量数据,用于后续的太阳阵I-V曲线和能量平衡分析;
[0010] 步骤二、太阳阵I-V曲线模型确定;
[0011] 根据太阳阵特性构建太阳阵计算模型,同时考虑太阳入射角、辐照衰减、日地因子、损失因子因素的影响,计算不同季节、不同轨道条件和不同工况下太阳阵输出电压、输出电流,以表征太阳阵输出功率实时及长期变化情况;
[0012] 以标准状态的I-V曲线特征点为参数,考虑多种环境因素对太阳阵寿命的影响,计算太阳阵的输出特性;利用公式(Equ.1)太阳阵I-V曲线的计算机解析模型,得到不同条件下的太阳阵的I-V特性曲线;该模型在光照强度小于2个太阳常数时,有很高的精确性;太阳同步轨道小卫星的光照情况满足这一条件:
[0013]
[0014] 式中:
[0015] I——输出电流,单位为A;
[0016] Isc'——太阳阵短路电流,典型参数或实测值,单位为A;
[0017] C1——公式系数1;
[0018] V——太阳阵输出电压,单位为V;
[0019] C2——公式系数2;
[0020] Vov'——太阳阵开路电压,典型参数或实测值,单位为V;
[0021] Imp'——太阳阵最佳工作点输出电流,典型参数或实测值,单位为A;
[0022] Vmp'——太阳阵最佳工作点输出电压,典型参数或实测值,单位为V;
[0023] 太阳阵开路电压和最佳工作点输出电压计算模型如下:
[0024]
[0025] 式中:
[0026] Vov——单体太阳电池开路电压,单位为V;
[0027] Vmp——单体太阳电池最佳工作点电压,单位为V;
[0028] βVBOL——单体太阳电池寿命初期电压温度系数,单位为V/℃;
[0029] KVRAD——太阳阵开路电压辐照衰降因子;
[0030] T——阳阵温度,单位为℃;
[0031] 太阳阵短路电流和最佳工作点电流计算模型如下:
[0032]
[0033] ISC——单体太阳电池短路电流,单位为A;
[0034] Imp——单体太阳电池最佳工作点电流,单位为A;
[0035] αI—单体太阳电池电流温度系数,单位为A/℃;
[0036] θ(t)—— 一圈轨道内太阳光线与太阳阵法线方向的夹角,单位为度;
[0037] T—太阳阵温度,单位为℃;
[0038] KIRAD—太阳阵短路电流辐照衰降因子;
[0039] Frd——日地距离因子;
[0040] 利用“太阳阵开路电压及短路电流辐照衰降因子计算模型”预测LEO轨道辐射环境对卫星太阳电池输出参数衰减的影响,在该模型中Isc即为KIRAD,Vov即为KVRAD;
[0041] A.模型输入参数定义如下:
[0042] 电池类型:单结GaAs太阳电池;石英玻璃盖片厚度:120μm;轨道高度:300km~3000km;倾角:只针对99°;时间单位:月;
[0043] B.模型输出参数定义如下:
[0044] 最大输出功率Pmax、短路电流Isc、开路电压Vov,其输出形式:给出Pmax、Isc和Vov经过m个月后,Pmax、Isc和Vov为初始值的百分比,即给出Pmax、Isc和Vov关于时间month的函数;
[0045] 以下为该太阳阵开路电压及短路电流辐照衰降因子的计算模型:
[0046] 不同轨道高度位移损伤剂量计算如下:x为轨道高度,month为在轨月数,y为计算得到的位移损伤剂量;
[0047] 当300km<=x<=600km时,计算公式为:
[0048] y=14(A0+A1·x+A2·x2+A3·x3+A4·x4+A5·x5)·month (Equ.4)
[0049] 其中,A0=-5.72637E6,A1=69074.68933,A2=-329.19032,
[0050] A3=0.77634,A4=-9.13546E-4,A5=4.49106E-7
[0051] 当600km
[0052] y=14(A0+A1·x+A2·x2+A3·x3+A4·x4)·month (Equ.5)
[0053] 其 中,A0=-5.80893E7,A1=321272.30685,A2=-663.23216,A3=0.59526,A4=-1.77968E-4当1000km
[0054] y=14(A0+A1·x+A2·x2+A3·x3+A4·x4+A5·x5)·month (Equ.6)
[0055] 其中,A0=5.01219E8,A1=-1.76649E6,A2=2453.54778,
[0056] A3=-1.65135,A4=5.32602E-4,A5=-5.18233E-8
[0057] GaAs/Ge太阳电池的Pmax、Isc和Voc的计算模型为:
[0058] 最大输出功率衰减,即Pmax的计算模型:
[0059] Pmax=1.0-C×log10(1+(y/Dx)) (Equ.7)[0060] 其中,C=0.242,Dx=3.47e9,y为计算得到的位移损伤剂量;
[0061] 短路电流衰减,即Isc的计算模型:
[0062] KIRAD=Isc=1.0-C×log10(1+(y/Dx)) (Equ.8)[0063] 其中,C=0.213,Dx=8.3e19
[0064] 开路电压衰减,即Voc的计算模型:
[0065] KVRAD=Vov=1.0-C×log10(1+(y/Dx)) (Equ.9)[0066] 其中,C=0.07,Dx=1.8e9
[0067] 步骤三、太阳阵能量平衡计算模型确定;
[0068] 在进行能量平衡计算时,根据在轨数据或地面提供的数据对能量平衡的临界状态进行实时监控,如果太阳阵提供能量的多余电量Qresidual(c)由正值转变为零,则表明太阳阵已处于严重损伤状态,且Qresidual(c)的计算式为:
[0069]
[0070] 其中:
[0071] Qresidual(c)——在轨第c圈太阳阵可提供的多余电量,单位为C;
[0072] te——阴影期时间,单位为s;
[0073] ISA(c)——在轨第c圈方阵电流箝位点电流值,单位为A;
[0074] Iload_mean(c)——光照期负载电流Iload(A),其为在轨第c圈负载电流每周期的平均值,单位为A;
[0075] Id(t)——阴影期,蓄电池放电电流,单位为A;
[0076] 根据指定时期太阳阵I-V曲线的方程,在给出相应的光照区母线电压VS_bus及太阳阵隔离二极管和供电线缆压降之和Vs_dioline时,得到该指定时期I-V曲线上太阳阵工作电压箝位点Vop1处的电流值Is_op1;由能量平衡计算可知,该指定时期的太阳阵提供能量的多余电量Qresidual(c)可表示为:
[0077]
[0078] 式中:
[0079] Is_op1——指点时期I-V曲线上太阳阵工作电压箝位点Vop1处的电流值,单位为A;
[0080] Is_load_mean——指定时期光照区负载电流/在轨所有负载电流数据的平均值,单位为A;
[0081] Id(t)——指定周期蓄电池在阴影区的放电电流值,单位为A;
[0082] 太阳阵工作电压点输出功率计算模型如下;
[0083] PsA(t)=Vbus(t)Iop1(t), Isa(t)=Iop1(t)
[0084] 进一步可得:
[0085] Ps_op1(c)=Vs_bus(c)Is_op1(c)
[0086] Vs_op1(c)=Vs_bus(c)+Vs_dioline (Equ.12)[0087] 式中:
[0088] t――第c圈轨道周期内时刻,从0
[0089] PSA——太阳阵输出功率,单位为W;
[0090] Ps_op1(c)——第c圈太阳阵输出为箝位点Vs_op1时输出功率,单位为W;
[0091] Is_op1(c)——第c圈I-V曲线上太阳阵工作电压箝位点Vs_op1时电流值,单位为A;
[0092] Vs_bus——光照区母线电压,单位为V;
[0093] Vs_dioline——太阳阵隔离二极管和供电线缆压降之和,单位为V;
[0094] 对于在轨电源系统,由于电源控制器使得光照期母线电压始终保持为定值,因而,可以认为:在电源控制器正常工作的前提下,光照期母线电压始终不变;同时,在得到Vs_dioline值后,可够通过上述建立的I-V曲线得到指定时期箝位点电流值Is_op(c),以用于能量平衡计算;
[0095] 如果系统预测得到某指定时期Qresidual(c)=0,则说明此时太阳阵已到寿;
[0096] 其中,进行蓄电池组放电电流计算时,蓄电池组的放电电流取决于蓄电池组的放电功率、放电调节器效率、蓄电池组供电线路损耗因子、电池组电压因素;
[0097] 阴影区,蓄电池组放电电流为:
[0098]
[0099] 式中:
[0100] t――第c圈轨道周期内时刻,从0
[0101] Iload(t)――在轨负载电流需求随时间变化的函数; ηBDR――放电调节器效率;
[0102] nine――蓄电池组供电线路损耗因子;
[0103] Vbat(t)――蓄电池组放电电压;基于在轨第c圈轨道蓄电池放电初压和放电终压值,可以近似地认为Vbat(t)由放电初压至放电终压线性变化;
[0104] Vbus:放电时母线电压,单位为V;
[0105] 其中:阴影区太阳阵ISA(t)电流为零,Vbat(t)是一个被积函数,该函数是由蓄电池在轨放电初压和放电终压确定的线性函数。
[0106] 其中,所述的太阳阵温度通过太阳阵温度模型计算如下:
[0107] 太阳阵温度随卫星进出影状态的变化而变化,地影区,太阳阵温度逐渐下降,直至降至出影前的最低温度;光照期,太阳阵温度从出影后迅速上升,直至达到光照期的温度平衡点,此后温度保持不变直至卫星进入下一轨道圈的地影期,周而复始;
[0108] 太阳阵温度变化的简化模型如下:
[0109] 在地影期内,太阳阵温度从光照期的最高平衡温度线性下降至地影期最低温度,出影后,太阳阵温度在8分钟内从地影期最低温度上升至60℃,在20分钟内从60℃上升至光照期最高平衡温度,直至下次进影;
[0110] 光照期的最高平衡温度、地影期最低温度的默认值分别为:
[0111] 光照期的最高平衡温度TSAS;地影期最低温度TSAE。
[0112] 本发明的优点为:
[0113] (1)、多寿命影响因素综合:本发明以试验为基础,综合考虑日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹角、温度、太阳辐射等重要空间环境因素影响下的太阳阵整体寿命预测模型;
[0114] (2)、模型通用性:对于太阳同步轨道小卫星太阳阵,包括:Si及GaAs各类型太阳电池阵的整体寿命预测,避免了只能对特定单体电池进行寿命预测的局限;
[0115] (3)、工程实用性:本申请的太阳阵寿命预测方法属于物理模型与数据驱动相结合的方法。在建模过程中采取了部分简化方式,同时,模型所需要的数据基于现有的小卫星所采集的参数,容易获取,简化了寿命预测模型的复杂度及数据获取困难等问题,进而具有更强的工程实用性。附图说明
[0116] 图1为各子系统失效导致小卫星失效比列分配关系图;
[0117] 图2为太阳阵寿命预测流程图
[0118] 图3为小卫星太阳电池特征参数衰减计算流程图;
[0119] 图4为太阳阵年日地距离因子变化规律曲线图;
[0120] 图5为太阳入射角年变化规律曲线图;
[0121] 图6为HY-1B小卫星太阳阵I-V曲线图;
[0122] 图7为太阳阵提供多余电量随时间变化曲线图。

具体实施方式

[0123] 下面结合附图和实施例对本发明进行详细说明。
[0124] 一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法,该方法通过如下步骤实现:
[0125] 步骤一、日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹角确定;
[0126] 根据轨道高度、降交点地方时和预测起始时间,调用STK软件,计算每天的日地距离因子、轨道周期Te、轨道地影时间、每轨太阳光线与太阳阵法线的夹角的变化规律,得到随时间变化的定量数据,用于后续的太阳阵I-V曲线和能量平衡分析。
[0127] 步骤二、太阳阵I-V曲线模型确定;
[0128] 根据太阳阵特性构建太阳阵计算模型,同时考虑太阳入射角、辐照衰减、日地因子、损失因子等因素的影响,计算不同季节、不同轨道条件和不同工况下太阳阵输出电压、输出电流,以表征太阳阵输出功率实时及长期变化情况。
[0129] 以标准状态的I-V曲线特征点为参数,考虑多种环境因素对太阳阵寿命的影响,计算太阳阵的输出特性。利用Equ.1太阳阵I-V曲线的计算机解析模型,可以得到不同条件下的太阳阵的I-V特性曲线。该模型在光照强度小于2个太阳常数时,有很高的精确性。太阳同步轨道小卫星的光照情况满足这一条件:
[0130]
[0131] 式中:
[0132] I——太阳阵输出电流,单位为A;
[0133] ISC′——太阳阵短路电流,典型参数或实测值,单位为A;
[0134] C1——公式系数1;
[0135] V——太阳阵输出电压,单位为V;
[0136] C2——公式系数2;
[0137] Vov'_____阵开路电压,典型参数或实测值,单位为V;
[0138] Imp′——太阳阵最佳工作点输出电流,典型参数或实测值,单位为A;
[0139] Vmp'——太阳阵最佳工作点输出电压,典型参数或实测值,单位为V。
[0140] (1)太阳阵典型特征点参数计算模型;
[0141] (1.1)太阳阵开路电压和最佳工作点输出电压:
[0142]
[0143] 式中:
[0144] Vov——单体太阳电池开路电压(QJ 1019-1995规定的AM0,25℃),单位为V;
[0145] Vmp——单体太阳电池最佳工作点电压(QJ 1019-1995规定的AM0,25℃),单位为V;
[0146] βVBOL——单体太阳电池寿命初期电压温度系数(太阳电池的温度改变1℃时,其输出电压的变化值),单位为V/℃;
[0147] KVRAD——太阳阵开路电压辐照衰降因子;
[0148] T——太阳阵温度,单位为℃。
[0149] (1.2)太阳阵短路电流和最佳工作点电流计算模型:
[0150]
[0151] ISC——单体太阳电池短路电流(QJ 1019-1995规定的AM025℃条件),单位为A;
[0152] Imp——单体太阳电池最佳工作点电流(QJ 1019-1995规定的AM0,25℃条件),单位为A;
[0153] αI——单体太阳电池电流温度系数(太阳电池的温度改变1℃时,其输出电流的变化值),单位为A/℃;
[0154] θ(t)——一圈轨道内太阳光线与太阳阵法线方向的夹角,单位为度;
[0155] T——太阳阵温度,单位为℃;
[0156] KIRAD——太阳阵短路电流辐照衰降因子;
[0157] Frd——日地距离因子。
[0158] (2)太阳阵开路电压及短路电流辐照衰降因子计算模型;
[0159] 本模型用于预测LEO轨道辐射环境对卫星太阳电池输出参数衰减的影响,模型中Isc即为KIRAD,Vov即为KVRAD。
[0160] A.模型输入参数定义如下:
[0161] 电池类型:单结GaAs太阳电池;石英玻璃盖片厚度:120μm;轨道高度:300km~3000km;倾角:只针对99°;时间单位:月;
[0162] B.模型输出参数定义如下:
[0163] 最大输出功率Pmax、短路电流Isc、开路电压Vov(输出形式:给出Pmax、Isc和Vov经过m个月后,Pmax、Isc和Vov为初始值的百分比,即给出Pmax、Isc和Vov关于时间month的函数);
[0164] 以下为该计算模型的计算式:
[0165] (2.1)不同轨道高度位移损伤剂量计算
[0166] x为轨道高度,month为在轨月数,y为计算得到的位移损伤剂量(质子位移损伤剂量)。
[0167] 当300km<=x<=600km时,计算公式为:
[0168] y=14(A0+A1·x+A2·x2+A3·x3+A4·x4+A5·x5)·month (Equ.4)
[0169] 其中,A0=-5.72637E6,A1=69074.68933,A2=-329.19032,
[0170] A3=0.77634,A4=-9.13546E-4,A5=4.49106E-7
[0171] 当600km
[0172] y=14(A0+A1·x+A2·x2+A3·x3+A4·x4)·month (Equ.5)
[0173] 其 中,A0=-5.80893E7,A1=321272.30685,A2=-663.23216,A3=0.59526,A4=-1.77968E-4当1000km
[0174] y=14(A0+A1·x+A2·x2+A3·x3A4·x4+A5·x5)·month (Equ.6)
[0175] 其中,A0=5.01219E8,A1=-1.76649E6,A2=2453.54778,
[0176] A3=-1.65135,A4=5.32602E-4,A5=-5.18233E-8
[0177] (2.2)GaAs/Ge太阳电池的Pmax、Isc和Voc的计算模型
[0178] 最大输出功率衰减(Pmax的计算模型):
[0179] Pmax=1.0-C×log10(1+(y/Dx)) (Equ.7)[0180] 其中,C=0.242,Dx=3.47e9,y为计算得到的位移损伤剂量。
[0181] 短路电流衰减(Isc的计算模型):
[0182] KIRAD=Isc=1.0-C×log10(1+(y/Dx)) (Equ.8)[0183] 其中,C=0.213,Dx=8.3e19
[0184] 开路电压衰减(Voc的计算模型):
[0185] KVRAD=Vov=1.0-C×log10(1+(y/Dx)) (Equ.9)[0186] 其中,C=0.07,Dx=1.8e9
[0187] (3)太阳阵温度模型;
[0188] 太阳阵温度随卫星进出影状态的变化而变化。地影区,太阳阵温度逐渐下降,直至降至出影前的最低温度;光照期,太阳阵温度从出影后迅速上升,直至达到光照期的温度平衡点,此后温度保持不变直至卫星进入下一轨道圈的地影期,周而复始。
[0189] 太阳阵温度变化的简化模型如下:
[0190] 在地影期内,太阳阵温度从光照期的最高平衡温度线性下降至地影期最低温度,出影后,太阳阵温度在8分钟内从地影期最低温度上升至60℃,在20分钟内从60℃上升至光照期最高平衡温度,直至下次进影。
[0191] 光照期的最高平衡温度、地影期最低温度的默认值分别为:
[0192] 光照期的最高平衡温度TSAS;地影期最低温度TSAE。
[0193] 步骤三、太阳阵能量平衡计算模型;
[0194] (1)能量平衡计算;
[0195] 根据在轨数据或地面提供的数据对能量平衡的临界状态进行实时监控,如果太阳阵提供能量的多余电量Qesidual(c)由正值转变为零,则表明太阳阵已处于严重损伤状态,且Qesidual(c)的计算式为:
[0196]
[0197] 其中:
[0198] Qresidual(c)——在轨第c圈太阳阵可提供的多余电量,单位为C;
[0199] te——阴影期时间,单位为s;
[0200] ISA(c)——在轨第c圈方阵电流箝位点电流值,单位为A;
[0201] Iload_mean(c)——光照期负载电流Iload(A)(在轨第c圈负载电流每周期的平均值),单位为A;
[0202] Id(t)——阴影期,蓄电池放电电流,单位为A;
[0203] 由‘太阳阵工作电压点输出功率计算模型’可知,根据指定时期太阳阵I-V曲线的方程,在给出相应的光照区母线电压Vs_bus及太阳阵隔离二极管和供电线缆压降之和Vs_dioline时,可以得到该指定时期I-V曲线上太阳阵工作电压箝位点Vop1处的电流值Is_opl。由能量平衡计算可知,该指定时期的太阳阵提供能量的多余电量Qesidual(c)可表示为:
[0204]
[0205] 式中:
[0206] Is_op1——指点时期I-V曲线上太阳阵工作电压箝位点Vop1处的电流值;
[0207] Is_load_mean——指定时期光照区负载电流/在轨所有负载电流数据的平均值;
[0208] Id(t)——指定周期蓄电池在阴影区的放电电流值;
[0209] 如果系统预测得到某指定时期Qresidual(c)=0,则说明此时太阳阵已到寿。
[0210] (2)太阳阵工作电压点输出功率计算模型;
[0211] PSA(t)=Vbus(t)Iop1(t), ISA(t)=Iop1(t)
[0212] 进一步可得:
[0213] Ps_OP1(c)=Vs_bus(c)Is_op1(c)
[0214] Vs_op1(c)=Vs_bus(c)+Vs_dioline (Equ.12)
[0215] 式中:
[0216] t――第c圈轨道周期内时刻,从0
[0217] PSA——太阳阵输出功率,单位为W;
[0218] Ps_OP1(c)——第c圈太阳阵输出为箝位点Vs_op1时输出功率,单位为W;
[0219] Is_opl(c)——第c圈I-V曲线上太阳阵工作电压箝位点Vs_op1时电流值,单位为A;
[0220] Vs_bus——光照区母线电压,单位为V;
[0221] Vs_dioline——太阳阵隔离二极管和供电线缆压降之和,单位为V;
[0222] 对于在轨电源系统,由于电源控制器使得光照期母线电压始终保持为定值,因而,可以认为:在电源控制器正常工作的前提下,光照期母线电压始终不变;同时,在得到Vs_dioline值后,可够通过上述建立的I-V曲线得到指定时期箝位点电流值Is_op1(c),以用于能量平衡计算。
[0223] (3)蓄电池放电电流计算;
[0224] 蓄电池组的放电电流取决于蓄电池组的放电功率、放电调节器效率、蓄电池组供电线路损耗因子、电池组电压等因素。
[0225] 阴影区,蓄电池组放电电流为:
[0226]
[0227] 式中:
[0228] t――第c圈轨道周期内时刻,从0
[0229] Iload(t)――在轨负载电流需求随时间变化的函数;
[0230] ηBDR――放电调节器效率
[0231] ηline――蓄电池组供电线路损耗因子
[0232] Vbat(t)――蓄电池组放电电压;基于在轨第c圈轨道蓄电池放电初压和放电终压值,可以近似地认为Vbat(t)由放电初压至放电终压线性变化;
[0233] Vbus:放电时母线电压,单位为V;
[0234] 其中:阴影区太阳阵ISA(t)电流为零,Vbat(t)是一个被积函数,该函数是由蓄电池在轨放电初压和放电终压确定的线性函数。本专利仅考虑太阳阵的寿命预测问题,即假定蓄电池及电源控制器工作状态一切正常或蓄电池及电源控制器性能衰退对太阳阵不产生影响,因而,此处蓄电池在轨放电初压和放电终压可用定值处理。
[0235] 实施例1如下:
[0236] 本发明是一种基于I-V曲线与能量平衡的小卫星太阳阵寿命预测方法,所述的寿命预测方法是把太阳阵性能衰退过程视为能量平衡的衰退过程,从而,在综合考虑了小卫星运行空间环境因素(包括:日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹角、温度、太阳辐射等)影响下,建立具有较好通用性的太阳阵综合寿命预测模型。图2所示为本发明的寿命预测方法的总体流程图,具体实施步骤如下:
[0237] 步骤一、日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹角确定;
[0238] 根据轨道高度、降交点地方时和预测起始时间,调用STK软件,计算日地距离因子、轨道周期、轨道地影时间、轨道太阳光线与太阳阵法线的夹角的变化规律,得到随时间变化的定量数据,用于后续的太阳阵I-V曲线和能量平衡分析。
[0239] 步骤二、太阳阵I-V曲线模型;
[0240] (1)太阳阵开路电压及短路电流辐照衰降因子计算模型;
[0241] 图3给出的是LEO轨道辐射环境对太阳电池特征参数衰减影响的定量计算流程。结合图3,首先,根据小卫星轨道高度计算位移损伤剂量随时间的变化函数y=f(month)(其中,相应参数参见Equ.4~Equ.6)。
[0242] 其次,利用公式Equ.8-9分别计算太阳阵开路电压辐照衰降因子随时间的变化函数KVRAD=gv(month)与太阳阵短路电流辐照衰降因子随时间的变化函数KIRAD=gI(month)。
[0243] (2)太阳阵典型特征点参数计算模型;
[0244] 对于计算公式Equ.2-3,其中,Vov,Vop,βVBOL,ISC,Imp,αI,Np Ns皆为已知参数数据;θ(t),Frd通过上述步骤已计算得到;KIRAD,KVRAD为month的函数;由太阳阵温度模型获得光照期太阳阵的温度数据。把上述参数数据及函数带入Equ.2可得太阳阵开路电压和最佳工作点输出电压、太阳阵短路电流和最佳工作点电流,可表示如下:
[0245]
[0246] 把公式Equ.14带入Equ.1,分别计算公式系数C1,C2,进而可建立太阳阵I-V曲线随时间变化的规律,可表示为:
[0247] I=V(month) (Equ.15)
[0248] 将I-V曲线上的三个特征点即开路电压点(Voc')、短路电流点(ISC')和最大功率点(Imp',Vmp')用不同条件下的光照、温度及辐照损失系数进行修正后,代入此解析表达式,即可得到不同条件下的太阳阵的I-V特性曲线。
[0249] 步骤三、太阳阵能量平衡计算模型;
[0250] 由Vs_bus,Vs_dioline值通过Equ.12计算Vs_op1值,由公式Equ.15计算得到相应的Is_op1值,进而建立Is_opl与时间month之间的函数关系;
[0251] 对于Equ.13中的参数ηBDR,ηline,Vbus为已知量,在阴影区ISA(t)为零,Iload(t)利用已有小卫星在轨负载电流数据的均值做近似处理,即Iload(t):=Iload_mean(c),其中,‘:=’表示‘定义为’,ISA(t)=0A。
[0252] 蓄电池放电电流计算:阴影期,蓄电池放电电流Id(t)公式中的Vbat(t)作为一个被积函数。该函数是由蓄电池在轨放电初压和放电终压确定的线性函数,且其放电初压恒定,放电终压由已有小卫星蓄电池放电终压均值近似。
[0253] 对于Equ.11Te,te已知,Vs_op1及 积分项通过上述描述可求,IS_load_mean(c)利用已有小卫星光照期负载电流均值近似。把上述参数,带入Equ.11,由于IS_load_mean(c)×(Te-te)及 近似为常数项,而Is_op1(c)由Equ.12及Equ.15确定,因
而,Is_op1(c)可表示为:Is_op1(c)=Is_op1(month),最终,得到QS_residual(c)=QS_residual(month),即得到太阳阵提供能量的多余电量随时间变化的函数,从而可确定太阳阵寿命。
[0254] 实施例2如下:
[0255] 本实施例以我国航天HY-1B小卫星太阳阵为对象,通过本实施例的详细阐述,进一步说明本发明的实施过程及工程应用过程。
[0256] 步骤一、日地距离因子、轨道地影时间、太阳光线与太阳阵法线的夹角确定;
[0257] 小卫星轨道高度645km、降交点地方时15:00PM±30min——根据轨道高度、降交点地方时和预测起始时间,调用STK软件,计算日地距离因子Frd、轨道周期、轨道地影时间、轨道太阳光线与太阳阵法线的夹角θ(t)的变化规律,得到随时间变化的定量数据,用于后续的太阳阵输出功率计算和能量平衡分析。
[0258] Frd=[1033 1.023 1.008 0.991 0.977 0.968 0.968 0.976 0.991 1.008 1.0241.033](一年中每个月的日地距离因子,如图4所示);轨道周期:Te=100.8min;轨道地影时间:te=33.5217min分钟;太阳入射角θ(t)具有周期性,如图5所示为太阳入射角年变化规律曲线。
[0259] 步骤二、太阳阵I-V曲线模型;
[0260] (1)太阳阵开路电压及短路电流辐照衰降因子计算模型;
[0261] 由于HY-1B小卫星的轨道高度为645km,依据公式Equ.5计算位移损伤剂量y随时间month变化规律:y=f(month)。利用公式Equ.8-9分别计算太阳阵开路电压辐照衰降因子随时间的变化函数KVRAD=gv(month)与太阳阵短路电流辐照衰降因子随时间的变化函数KIRAD=gI(month)。
[0262] (2)太阳阵典型特征点参数计算模型;
[0263] 对于计算公式Equ.2-3,其中,Vov=2.65V,Vmp=2.32V,βVBOL=-6.7mV/°C,ISC2
=0.396A,Imp=0.375A,aI=0.014mA/cm°C,Np=114,Ns=18,T=82°C(由太阳阵
温度模型可知,光照期,太阳阵温度可迅速上升至平衡温度,因而,此处以最高平衡温度TSAS=82°C作为HY-1B小卫星太阳阵温度);θ(t),Frd由步骤一给出,且KVRAD=gv(month),KIRAD=gI(month)。进而,分别建立形如Equ.14的表达式,把得到的表达式带入Equ.1,并最终建立HY-1B小卫星太阳阵I-V曲线随时间变化的规律,图6为HY-1B小卫星太阳阵寿命
周期内某一给定时刻的I-V曲线。
[0264] 步骤三、太阳阵能量平衡计算模型;
[0265] 由Vs_bus=29.3V,Vs_dioine=2.0 V值通过Equ.12计算Vs_op1=31.3V值,由公式Equ.15计算得到相应的Is_op1值,进而建立Is_op1与时间month之间的函数关系Is_opl=Is_opl(month);
[0266] 对于Equ.13中的参数ηBOR=0.92,ηline=0.9673,Vbus=28.6V为已知量,Iload(t),在阴影区IsA(t)=0A,经计算得: Iload(t):=Iload_mean(c)=9.0A,其中,‘:=’表示‘定义为’。
[0267] 阴影期,蓄电池放电电流Id(t)公式中的Vbat(t)作为一个被积函数。该函数是由蓄电池在轨放电初压和放电终压确定的线性函数,且其放电初压恒定,放电终压由已有小卫星蓄电池放电终压均值近似,且放电初压为1.383 V*18(蓄电池组串联单体蓄电池个数)=24.8940V,放电终压均值21.6V。
[0268] 把Te=100.8min, te=33.5217min,及通过计算得到的Vs_op1及 IS_load_mean(c)=11.4A带入Equ.11,由于IS_load_mean(c)×(Te-te)及 近似为常数
项,而Is_op1(c)由Equ.12及Equ.15确定,因而,Is_op1(c)可表示为:Is_op1(c)=Is_op1(month),由此,Equ.11转化为QS_residual(c)=QS_residual(month),即得到太阳阵提供能量的多余电量随时间变化的函数,图7所示为多余电量随时间的变化曲线。由图可知,多余电量随着日地距离因子,太阳入射角等周期影响因素的综合作用下,呈现明显的周期性,且多余电量持续减少,与平线相交处即为太阳阵使用寿命终止位置
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈