Preliminary Study on Storage Life Distribution of Semiconductor Device Based on Weibull Distribution
-
摘要: 对有机电致发光二极管（Organic Light-Emitting Diode，OLED）微型显示器件进行90℃、80℃、70℃的高温贮存试验，获得产品的失效数据。基于威布尔分布模型，采用最小二乘法进行参数估计，对失效数据分析，获得OLED微型显示器件失效分布函数。应用经典可靠性理论，计算产品在90℃、80℃、70℃的特征寿命、可靠寿命及平均故障间隔时间（Mean Time Between Failure，MTBF）。采用Arrhenius模型，依据90℃、80℃、70℃的贮存特征寿命，获得常温下产品的贮存特征寿命。分析结果表明，该方法合理、简便、有效，数据结果可以进一步应用到推导产品常温贮存寿命。Abstract: In this study, a semiconductor device was tested for high temperature storage at 90℃, 80℃ and 70℃, and the failure data is obtained. Based on the Weibull distribution model, parameter estimation was carried out by the least square method. The failure distribution function of the semiconductor device was obtained. And the classical reliability theory was applied to calculate the characteristic life, reliable life and MTBF of the product at 90℃, 80℃ and 70℃. Using the Arrhenius model, the storage characteristic life of the semiconductor device at room temperature was obtained, according to the storage characteristic life of 90℃, 80℃ and 70℃. The results show that the method is reasonable, simple and effective, and the results can be used to derive the normal temperature storage life.
-
0. 引言
随着遥感器技术的进步，空间遥感器相机向着多谱段、大幅宽、实时传输等方向不断发展。获取多个谱段的景物信息对资源普查、矿物勘测、农作物监测等具有极高的应用前景[1]。多光谱遥感相机具有较宽谱段的成像能力，美国于2008年发射Geoeye-1高分辨率卫星，运行于高度681 km，轨道倾角98°的太阳同步轨道。其全色谱段地面像元分辨率为0.41 m，包括4个经典波段（红、绿、蓝和近红外）多光谱的地面像元分辨率为1.65 m[2]。国内于2020年成功发射资源三号03星，运行于轨道高度506 km的太阳同步轨道，地面分辨率下视2.1 m，前后视2.5 m，多光谱5.8 m，覆盖宽度单载荷大于50 km，立体成像大于45 km[3]。
本文所介绍的空间相机运行于近地圆轨道，相比于太阳同步轨道，全寿命周期内β角（轨道面与太阳光矢量之间的夹角）变化很大，在－58.4°~58.4°之间，相机安装于卫星舱体外，造成相机每个面都有可能受照[4]。相机在轨可实现全天时工作，工作期间姿态多变且伴随着大角度的俯仰与侧摆，轨道外热流非常复杂，热流扰动大。相机拥有可见光与红外双成像通道，在外热流变化复杂的空间环境下，相机需要对以下方面开展散热设计：可见光焦面需要在较低的温度水平下工作，以减小暗电流和热噪声对成像质量的影响[5]；低温红外焦面所配备的压缩制冷机功耗大，且工作时间长。为保证其在轨期间正常进行制冷工作，需将其产生的热量快速排散。
为满足该空间相机内热源的散热要求，本文对该相机所受到的外热流进行分析并根据空间环境进行了散热设计。针对热设计开展热仿真计算，验证了相机散热设计的正确性与合理性。
1. 相机结构及外热流分析
1.1 相机概述及散热设计分析
该相机采用同轴三反光学系统，主体构型示意图与主要内热源分布情况见图 1，内热源由可见光组件、红外组件和制冷机组件等组成。
![]()
图 1 相机主要结构及内热源分布示意图Figure 1. Schematic diagram of the main structure of the camera and the distribution of internal heat sources相机内热源中红外组件与制冷机组件需要长期工作。其中制冷机组件包含制冷机热端部分与制冷机压缩机部分，压缩机部分安装在相机本体外，固定在卫星舱板上。相机工作热耗219 W，待机热耗为110 W，温控指标如表 1所示。
表 1 内热源温控指标Table 1. Temperature index of the internal heat source℃ Component Operating temperature Visible components -5~40 Infrared components -25~55 Refrigerator compressor -15~30 refrigerator hot end ≤20 1.2 外热流分析
该相机运行于近地圆轨道，β角在－58.4°~58.4°之间变化。相机在轨运行姿态为：阳照区非工作时段－Z对日，工作期间调整为+Z对地，阴影区长期+Z对地，+Z向为相机光线入射方向的反向，+X向为速度矢量方向，+Y向由右手坐标系确定。其中，在阳照区成像过程中相机入光口方向需要与太阳光矢量方向成最小30°规避角，避免太阳光入射相机入光口影响成像。
根据相机飞行轨道参数分析轨道外热流大小和变化规律，是确定相机散热面方位和定义极端工况的重要输入条件[6]。利用热分析软件计算出相机在冬至时刻运行姿态下随β角变化的平均到达外热流，如图 2所示。
相机+Z面、+X面和－X面平均到达外热流变化趋势相同，都随着β角绝对值增大呈下降的趋势。－Y面平均到达外热流随β角增大呈下降趋势，+Y面外热流变化与－Y面相反，±Y面平均到达外热流以β角0°呈对称趋势，且变化较大。
相机成像过程中，由于规避角影响，使相机±X和±Y面受到太阳直射，相较于正常姿态其外热流会发生较大变化。由于相机在不同方向摆动时，－X面与+X面，－Y面与+Y面到达外热流变化趋势相同，故以单向摆动为例，在极限β角（0°与58.4°）下计算±X与±Y面平均到达外热流，如图 3、图 4。
±Y面在β角为58.4°时受阳光直射时间长，平均到达外热流大，相较于β角为0°下的±X面，热环境恶劣，不利于内热源的散热设计。而±X面虽然均会受到阳光直射，但在相机在轨运行中存在单一侧面长时间背阳的情况，相较于±Y面，到达外热流小，使内热源散热压力较小。根据相机结构布局与到达外热流情况，将相机散热面设置在±X侧。
2. 相机散热设计
针对以上相机情况，结合相机所处空间环境与内热源特性，总结以下散热设计难点：
相机所处空间环境外热流变化复杂。在相机处于工作姿态时，除入光口外相机各面均会受到阳光直射。由于相机在轨期间姿态多变，使其受到的外热流变化复杂。
相机具有可见光与红外双成像通道，结构复杂，制冷机压缩机部分位于相机本体外，内热源位置分散。部分内热源在轨期间处于长期工作状态，发热量大，且相机所处空间环境复杂，不利于大功率热源散热。
针对以上相机散热设计的难点，同时兼顾相机内热源的热排散与热补偿，采用被动热控和主动热控相结合的热设计方案。
2.1 散热系统设计
为保证相机内热源产生的热量有效排散，需要对其进行被动热控设计。由于无法在相机本体上开设散热面，热设计中额外设置了辐射冷板进行散热。辐射冷板表面喷涂KS-ZA白漆来增大散热效果，内热源通过槽道热管将热量传导至辐射冷板，并最终辐射至外部空间。冷板与卫星舱体通过聚酰亚胺隔热垫与钛合金螺钉进行隔热安装，冷板与热管，热管与热管工装面之间填充导热填料，增大热量传导效率[6]。
相机各部分内热源工作模式与工作指标各不相同，需要设置多个辐射冷板进行散热设计。可见光组件热耗65 W，红外组件总热耗24 W，二者温度指标相近，且可见光组件单轨工作60 s，散热压力较小，使用一块辐射冷板对可见光组件与红外组件同时散热。可见光组件通过两根热管与辐射冷板相连，三部分红外组件各通过一根热管与辐射冷板相连，辐射冷板设置在相机－X侧。制冷机压缩机部分长期工作，其待机功耗为38.5 W，成像工作时功耗为45.5 W。使用一块辐射冷板对压缩机散热，压缩机通过一根热管与辐射冷板相连，辐射冷板设置在相机+X侧。制冷机热端部分是相机内部最大热源，其待机功耗为71.5 W，成像工作时功耗为84.5 W。由于制冷机热端长期工作，且功耗大，而相机在成像过程中±X面均有可能受到太阳直射，单一辐射冷板不利于热量排散，并且由于单侧设置辐射冷板受到的外热流大，使辐射冷板面积过大，不利于相机的结构设计。采用双辐射冷板耦合散热的方案，可以利用成像过程中相机一侧背阳，到达外热流较小的辐射冷板提升制冷机热端的散热效率，同时减小了辐射冷板面积，有利于相机的结构设计。辐射冷板耦合设计如图 5所示，制冷机热端通过两根热管将热量分别引至±X侧辐射冷板，两侧的辐射冷板通过一根U型热管连接。
![]()
图 5 制冷机热端散热方案示意图Figure 5. Schematic diagram of the heat dissipation for refrigerator hotend为使短期工作热源与长期工作热源均满足温度指标与合理规划热控资源，需要精确计算辐射冷板面积并代入仿真计算进行校核，最终确定辐射冷板的面积[6]。根据外热流与辐射冷板涂层属性，计算辐射冷板吸收热量Q1为：
$$ Q_1=Q_{\text {内 }}+\alpha \times\left(q_1+q_2\right) \times A+\varepsilon \times q_3 \times A $$ (1) 式中：Q内为内热源功耗；α为辐射冷板太阳吸收率；q1为太阳直射的到达外热流密度；q2为阳光反照的到达外热流密度；ε为辐射冷板表面发射率；q3为地球红外的到达外热流密度；A为辐射冷板面积。
辐射冷板散出的热量Q2为：
$$ Q_2=A \times \varepsilon \times \sigma \times T^4 $$ (2) 式中：σ为斯特藩-玻尔兹曼常数，其值为5.67×10－8 W/(m2·K4)[7]。
当热平衡时，有：
$$ Q_1=Q_2 $$ (3) 通过仿真计算得出，可见光组件与红外组件辐射冷板面积为0.19 m2，制冷机压缩机辐射冷板面积为0.26 m2，制冷机热端辐射冷板两面各0.25 m2，总面积0.5 m2。
2.2 热补偿设计
主动热控是在被动热控的基础上进一步提高并维持内热源温度水平的热设计方法，采用补偿加热的方式控制内热源温度，使其在低温工况下或待机状态下也满足温控指标[8-9]。补偿加热回路采用趋势控温法进行闭环控温，当温度低于温度指标时，回路连接，加热片工作，当温度达到温度指标时，回路断开，加热片停止工作。加热回路设置在各个辐射冷板上，保证内热源温度指标的同时也可以确保热管处于正常工作温区内。加热片使用聚酰亚胺薄膜型加热片，热敏电阻测温，由相机下位机采集控温。
2.3 传热路径
可见光组件内热源包含两部分，其中一部分通过其安装框架将热量传导至电箱基板，通过安装在电箱基板上的热管将热量传导至辐射冷板，最终辐射至外部空间。另一部分通过一个铝扩热板将热量集中传导至电箱外壳，通过热管将热量传导至辐射冷板并辐射至外部空间。红外组件上的发热元件与红外组件电箱外壳直接接触，接触界面填充导热填料增强导热，电箱外壳通过热管将热量传导至辐射冷板，最终将热量辐射至外部空间。制冷机热端产生的热耗直接由安装在其表面的两根热管传导至相机两侧辐射冷板，两侧的辐射冷板同时将热量辐射至外部空间。制冷机压缩机安装在相机本体外，通过铝合金安装架固定在卫星舱板上。压缩机产生的热耗传导至压缩机上压盖，接触面之间填充导热填料增强导热，通过安装在上压盖上的热管传导至辐射冷板，并辐射至外部空间。相机各内热源传热路径如图 6所示，图中点划线为传导导热路径，实线为热管传热，虚线为辐射传热路径。
导热热阻分为材料自身的传导热阻与接触面之间的接触传热热阻，计算公式分别为：
$$ {R_{\text{c}}} = \frac{\delta }{{{A_{\text{c}}} \cdot \lambda }} $$ (4) $$ {R_{{\text{ct}}}} = \frac{1}{{K \cdot {A_{{\text{ct}}}}}} $$ (5) 式中：Rc为材料自身的传导热阻；δ为热量传导方向的长度；Ac为横截面面积；λ为材料导热率；Rct为接触传热热阻；K为接触传热系数；Act为接触面积[10]。
3. 热仿真计算
3.1 建立热仿真模型
根据上述热设计方案，利用有限元软件对相机进行建模，如图 7所示。壳单元划分，单元厚度按照等效厚度计算，采用热耦合的方式对模型简化，共建立了209个热耦合，划分了16437个单元。
3.2 工况
仿真计算中根据卫星太阳位置、β角和散热面涂层表面属性定义高温与低温两个极端工况，如表 2所示。
表 2 相机工况定义Table 2. Definition of camera working conditionsSolar position β Heat sink surface coating surface properties High temperature working condition Winter solstice 0° α/ε=0.23/0.92（End period） Low temperature working condition Summer Solstice 60° α/ε=0.13/0.92（Initial period） 3.3 计算结果
根据工况定义，利用仿真软件计算相机在各工况下的温度结果，单轨结束时热源温度分布云图与热源温度变化曲线如图 8、9所示。
![]()
图 8 高温工况温度分布云图与温度变化曲线Figure 8. High temperature distribution clouds and temperature variation curves![]()
图 9 低温工况温度分布云图与温度变化曲线Figure 9. Low temperature distribution clouds and temperature variation curves从图 8、9可以看出，相机在温度稳定后高温工况下可见光组件温度为－1.1℃~12.9℃，红外组件温度为3℃~10.5℃，制冷机热端温度为－5℃~0.3℃，制冷机压缩机温度为－1℃~21.3℃。相机在低温工况下可见光组件温度为－1.9℃~9℃，红外组件温度为1.7℃~6.7℃，制冷机热端温度为－12℃~－10.8℃，制冷机压缩机温度为－11.3℃~5.4℃。均满足热控指标要求。
4. 结论
本文为解决相机在外热流复杂的空间环境下相机大功率热源的散热问题，通过分析相机各个侧面到达外热流情况与相机内热源工作模式，采用在相机两侧设置辐射冷板的散热方式，并使用热管将两侧辐射冷板耦合的方案，增大热量的排散效率，减小了辐射冷板的面积。最终，通过热仿真计算验证了设计的合理性。仿真结果表面：可见光组件温度为－1.9℃~12.9℃，红外组件温度为1.7℃~10.5℃，制冷机热端温度为－12℃~0.3℃，制冷机压缩机温度为－11.3℃~21.3℃，均满足温控指标要求，解决了复杂外热流下相机大功率热源的散热问题。
-
[1] 王忆锋, 唐利斌, 岳清. OLED器件寿命衰退模型的MATLAB分析计算[J].电子器件, 2010, 33(4): 412-415. http://www.cnki.com.cn/Article/CJFDTotal-DZQJ201004003.htm WANG Yifeng, TANG Libing, YUE Qing. MATLAB Analysis and Computation of OLED Lifetime Degradation Model[J]. Electron Devices, 2010, 33(4): 412-415. http://www.cnki.com.cn/Article/CJFDTotal-DZQJ201004003.htm
[2] 张方晖, 席俭飞, 王秀峰.硫系玻璃薄膜封装层对OLED寿命的影响[J].半导体光电, 2010, 31(1): 30-33. http://www.cnki.com.cn/Article/CJFDTotal-BDTG201001008.htm ZHANG Fanghui, XI Jianfei, WANG Xiufeng. Influence of Sulfide Glass Thin Film Encapsulation Layer on Lifetime of OLEDs[J]. Semiconductor Optoelectronics, 2010, 31(1): 30-33. http://www.cnki.com.cn/Article/CJFDTotal-BDTG201001008.htm
[3] Eisenbrand F, Karrenbauer A, XU C. Algorithms for longer OLED lifetime[J]. Lecture Notes in Computer. Science, 2007, 4525: 338-351. DOI: 10.1007/978-3-540-72845-0_26
[4] HAN C W, KIM H K, PANG H S, et al. Dual plate OLED display(DOD) embedded with white OLED[J]. J. Display Technol., 2009, 5(12): 541-545 DOI: 10.1109/JDT.2009.2024009
[5] NElSON W B. Accelerated Life Testing-step-stress Model Sand Data Analysis[J]. IEEE Transactions on Reliability, 1980, 29(2): 103-108.
[6] 陈循, 张春华, 汪亚顺, 等.加速寿命试验技术与应用[M].北京:国防工业出版社, 2013. CHEN Xun, ZHANG Chunhua, WANG Yashun, et al. Technology And Application of Accelerated Life Test[M]. Beijing: National Defence Industry Press, 2013.
[7] ELSAYED A. Reliability Engineering[M]. 2nd (SECOND EDITION), 1984.
[8] 陈光宇, 张文, 张小民.威布尔分布下系统全寿命周期成本建模与决策[J].系统工程理论与实践, 2016, 36(11): 2933-2940. http://www.cnki.com.cn/Article/CJFDTOTAL-XTLL201611020.htm CHEN Guangyu, ZHANG Chunhua, WANG Yashun. Modeling and decision-making of system lifecycle costs under Weibull distribution[J]. Systems Engineering-Theory & Practice, 2016, 36(11): 2933-2940. http://www.cnki.com.cn/Article/CJFDTOTAL-XTLL201611020.htm
[9] ROLINK. Infrared Focal Plane Array Storage Life Assessment By Accelerated Aging[J]. Quality and Reliability Engineering International, 1998, 14(6): 425-432. DOI: 10.1002/(SICI)1099-1638(199811/12)14:6<425::AID-QRE226>3.0.CO;2-J
[10] 蒋仁言.可靠性模型及应用[M].北京:机械工业出版社, 1999. JIANG Renyan. Reliability Model and Application[M]. Beijing: China Machine Press, 1999.
[11] 赵宇.可靠性数据分析[M].北京:国防工业出版社, 2014. ZHAO Yu. Reliability Data Analysis[M]. Beijing: National Defence Industry Press, 2014.
下载:
计量
- 文章访问数: 297
- HTML全文浏览量: 125
- PDF下载量: 30
