专利名称：一种进行卫星表面带电分析的方法
技术领域：
本发明涉及一种进行卫星表面带电分析的软件，具体地，涉及一种对地球同步轨 道卫星高压太阳电池阵表面带电情况进行仿真分析的软件。
背景技术：
国外已经开发出了多款用于卫星表面带电分析的软件，但由于技术封锁，我们只 能开发自有知识产权的卫星表面带电分析软件，目前该工作在国内还没有取得突破。地球同步轨道卫星高压太阳电池阵的表面带电问题具有一定代表性，主要表现在 带电环境和卫星太阳电池阵工作达到上百伏。目前对卫星带电无法从卫整星角度进行分 析，多根据郎缪尔探针原理进行简单的估计，不仅针对性不强而且误差较大。根据目前的方 法，通常是给出一个单一的充电电位数值，无法给出更详细的分析结果。
发明内容
因此，本发明解决的问题就是克服现有技术的不足，提供了一种较全面的卫星带 电分析软件，可以在针对不同光照条件、空间等离子体环境、卫星主体材料和高压太阳电池 阵材料，通过数值计算，以图表形式给出卫星主体结构、高压太阳帆板和太阳电池阵玻璃盖 片的电子产生率计算结果，卫星表面电位或充电电流随时间的变化结果，卫星太阳电池阵 表面带电电位的二维分布结果，并生成详细的技术报告。整个参数的设置均通过软件图形用户界面完成，可以将已有的空间环境参数、卫 星主体材料参数和高压太阳电池阵材料参数集成到软件中，并可以根据要求自行设定上述 参数。本发明的技术解决方案是通过建立电流收集模型，根据卫星表面带电电流平衡 方程，通过编程完成数值计算，得到平衡电位。在确定太阳电池阵表面电位分布的基础上， 结合空间环境特性和材料特性，分析在典型磁暴和亚暴条件下电池阵与空间等离子体之间 放电的可能性，也就是要确定各种空间环境条件下的放电阈值，最后实现了计算结果的可 视化。本发明的原理是确定卫星表面电位的理论基础是电流平衡方程，即流到卫星的 所有电流总和为零时，卫星应该处于什么样的电位。因此，首先要解决的问题是确定各种电 流的解析表达式。在GEO轨道，卫星收集电流的理论基础是稳定的、未磁化等离子体的厚鞘层近似， 这是因为电子和离子的回旋半径即德拜长度大于卫星尺寸，而卫星速度小于周围等离子体 的热速度。卫星表面带电的物理过程是非常复杂的。如果假定等离子体符合麦克斯韦分布， 则地球同步轨道处的表面带电公式变得十分简单。所有卫星表面充电的物理过程都是一个 电流平衡的过程，即在恒稳态下，所有电流的总和为零。卫星表面的净电流包含多个分量，并且计算时需要考虑多种因素对它的影响。只有当卫星处于恒稳态时，其表面的电势才能达到均衡状态，而这种恒稳态在地球同步轨道 是很难遇到的。卫星表面电势达到均衡状态时，各种电流的总和将为零。电流总和为零时 可以得到卫星处于恒稳态时的表面电位。通过考察该电位的大小即可评估卫星是否存在放 电危险。本发明与现有技术相比的优点在于在本发明中，通过图形用户界面实现了参数 设置和计算结果的的可视化，并能在局部材料参数发生改变的情况下快速完成仿真计算， 自动生成图表形势的详细技术报告，计算结果与EQUIP0T比较误差不超过30%。
图1为卫星平衡电位及充电时间算法流程图充电初始电位取VO，初始时间取0，计算初始平衡电流JO，判断初始平衡电流正负 号来决定迭代步长电位的正负，即判断是向正电位方向还是负电位方向充电，确定初始迭 代电位步长Δν，用来减少迭代次数。判断平衡电流Jv和初始充电电流JO同号的目的是确 保整个充电向同一个方向进行充电，异号表明充电电位超过平衡电位，此时需要调整充电 步长，调整步长流程为黄色标示。计算充电时间t和平衡电位Vs，电位平衡的条件为迭代电位步长小于充电电位 的 0. 001。参数说明V:充电电位ΔΥ 充电电位迭代算法所取电位步长，在程序中根据条件在不断的调整Δ T 总充电时间At 卫星表面充电AV对应的充电时间Δ Jv:充电电位为V时的平衡电流AVs:平衡电位
具体实施例方式研究卫星表面充电的理论基础是电流平衡方程，通过确定各种电流的解析表达 式，来求解卫星表面的电位分布。在地球同步轨道上，与卫星表面相互作用的空间等离子体 可以近似为厚鞘层，这是由于空间等离子体中的电子和离子的回旋半径即德拜长度大于卫 星尺寸，而卫星速度小于周围等离子体的热速度。建立电流收集模型卫星表面收集的电子电流、离子电流、光效应电流、次级电流、后向散射电流等等 效电流有可能在卫星上的某个电位上达到动态平衡。此时，卫星在整个表面上的电流代数 和为零，即达到电流平衡的状态。电流方程可以表示成
权利要求
1. 一种进行卫星表面带电分析的方法，其中包括以下步骤 1)建立电流收集模型卫星表面收集的电子电流、离子电流、光效应电流、次级电流、后向散射电流等等效电流有可能在卫星上的某个电位上达到动态平衡。此时，卫星在整个表面上的电流代数和为零，即达到电流平衡的状态。电流方程可以表示成 dVCA — = Jnet (V) 二 Λ (V) - U (η + L· (η + Λ, (V) + Jbe (V) + Jph (V) + Jc (V)] = 0ατ (1)其中V 卫星表面充电电位 t 卫星表面充电时间 jnet 到达卫星表面的净电流 je 电子电流 ji 质子电流jse 入射电子产生的次级电子电流 jsi 入射质子产生的次级电子电流 jbe 后向散射电子电流 jph 光电流 jc 电导电流CA 表面单位面积相对周围空间等离子体的充电电容电流方向以电子电流方向为正； 确定平衡状态的是表面电位V，求解(1)式找到一个这样的电位，使得净电流为零； 根据卫星表面带电电流平衡方程，通过编程完成数值计算，其中，充电初始电位取V0， 初始时间取0，根据电流平衡方程，计算初始平衡电流J0，通过判断初始平衡电流正负号来 决定迭代步长电位的正负，确定初始迭代电位步长。判断平衡电流Jv和初始充电电流JO 同号的目的是确保整个充电向同一个方向进行充电，异号表明充电电位超过平衡电位，通 过调整充电步长，不断对这个过程进行迭代，当满足迭代电位步长的条件时，获得充电时间 t和平衡电位Vs。在确定太阳电池阵表面电位分布的基础上，主要通过计算所采取的环境、材料特性参 数，获得卫星星体和电池阵与空间等离子体之间相互作用产生的充电电位V与充电时间t 的关系，并依据两者之间的击穿阈值电压，判断它们之间放电的可能性； 最后，将计算结果可视化，方便进行分析。
全文摘要
本发明是一种用于进行卫星表面带电分析的方法，集中解决地球同步轨道卫星高电压太阳电池阵表面带电问题。利用软件图形用户界面设定卫星轨道环境参数、卫星主体材料参数、卫星高压太阳电池阵材料参数，可以实现对卫星表面带电的仿真分析，并自动生成技术报告。
文档编号G06F17/50GK102117345SQ20091021715
公开日2011年7月6日 申请日期2009年12月31日 优先权日2009年12月31日
发明者冯伟泉, 刘业楠, 徐炎林, 王志浩 申请人:北京卫星环境工程研究所