基于有限元的X射线聚焦望远镜结构分析

王赛伊圣振王新穆宝忠

摘要：

针对我国X射线时变与偏振卫星（XTP）项目中高能X射线聚焦望远镜的研制需求，初步设计了嵌套式、类WolterI型X射线望远镜并进行了结构优化。基于X射线望远镜的光学设计理论，给出了嵌套式聚焦望远镜的光学结构参数。利用有限元分析软件ABAQUS，建立了嵌套式X射线聚焦望远镜的有限元模型。采用基础激励法模拟了X射线聚焦望远镜模型的随机振动，得到了望远镜结构的应力分布。通过切换层的方式，降低了望远镜结构的应力，达到了多层膜承受应力的要求。该研究为XTP项目中的X射线望远镜载荷设计提供了参考。

关键词：

X射线聚焦望远镜； 光学设计； 有限元分析； 随机振动

中图分类号： TH 743文献标志码： Adoi： 10.3969/j.issn.10055630.2015.05.013

引言

X射线天文望远镜是观测宇宙中辐射X射线天体的重要设备，探测目标包括黑洞、脉冲星、X射线双星、超新星遗迹和太阳活动区等。自上世纪70年代以后，美国航天局（NASA）、日本航天局（JASA）和欧洲航天局（ESA）先后发射了多个聚焦式X射线望远镜，如Einstein望远镜、Chandra望远镜、ASCA望远镜、XMMNewton望远镜、Nustar望远镜等。目前正在研制中的X射线望远镜项目包括日本的AstroH望远镜、美国Nicer望远镜[7]、欧洲的LOFT望远镜和Athena望远镜等。近年来，以黑洞及中子星等为观测目标的X射线时变与偏振探测卫星（XTP）的相关研究已在国内快速展开，X射线聚焦望远镜是该规划中卫星的重要载荷之一。

为了扩展工作于掠入射条件下的X射线望远镜的集光面积，其光学系统是在WolterI型的基础上，采用多层嵌套的结构形式。由于X射线望远镜载荷在发射升空过程中，将受到火箭发动机工作噪声及气动力激振等因素的影响，会产生随机振动载荷，影响望远镜结构的可靠性。为保证X射线望远镜在发射时随机振动环境下的可靠性，如何设计嵌套集成超薄反射镜的望远镜结构一直是X射线聚焦望远镜研究中的关键问题。在结构设计阶段，有限元方法是开展此类分析的重要手段。如AstroH望远镜设计中，基于有限元模态分析方法得到了望远镜反射镜的振型频率为64 Hz时，反射镜的边缘变形比中心处大；振型频率为118 Hz时，反射镜的中心变形比边缘处大；在Nustar望远镜设计中，通过随机振动有限元模拟分析了望远镜物镜的应力分布，得到了133层的嵌套式结构中最内5层为受力较集中区域且应力较大，其应力值约为2 MPa，从而选择比外层更宽的石墨垫片（1.6 mm）进行装配连接，以增强望远镜内层结构的支撑作用。本论文将结合XTP卫星对聚焦式X射线望远镜的光学性能要求，采用有限元分析软件ABAQUS来研究望远镜模型在受到随机激励下的动力学响应，从而获得望远镜结构的应力分布，为XTP项目中X射线聚焦望远镜载荷的结构设计和装配参数提供支持。

根据科学目标和卫星平台的要求，初步确定以下初始参数：焦距4.5 m；望远镜的外口径由天文观测区决定，针对6 keV的主观测能段，最大工作掠入射角可以达到约0.7°，对应的相对孔径D/F=0.1，因此选择最外层直径为450 mm；最内层口径主要依据热弯法制备高精度反射镜的技术能力，选择Din=100 mm作为初始参数；镜片采用厚度为0.21 mm的D263玻璃。结合式（2）～（4），设计得到望远镜的初始光学结构参数，如表1所示。

2.1嵌套式结构的有限元模型

望远镜的装配体以中心芯轴为支撑，由内向外逐层装配而成。层与层之间以石墨条为连接体，根据光学设计结果，将每层的石墨条加工成对应的口径和倾角。X射线望远镜的初始结构模型如图2所示。内层结构由6个60°的扇区模块构成，外层由12个30°的扇区模块构成。由于大口径反射镜的面形极难控制，根据目前的制备水平，外层反射镜选择了30°扇角。60°扇区模块和30°扇区模块的分界为第53层和54层，类似于133层的Nustar望远镜在第67层分界。外层的30°扇区模块中，每个反射镜由5个石墨条支撑，石墨条间隔扇角为7.5°；内层的60°扇区模块中，每个反射镜同样由5个石墨条支撑，石墨条间隔扇角为15°。反射镜采用肖特（Schott）公司的D263玻璃。芯轴采用5 mm厚度的钛合金材料，外径为100 mm。芯轴、石墨垫片和玻璃反射镜之间使用型号为F131的环氧树脂粘结。材料的关键参数如表2所示。

根据光学设计所得到的每一层反射镜和石墨垫片的曲率半径，使用ABAQUS预处理模块，采用实体建模的方式，建立每一层的反射片和石墨垫片。在ABAQUS软件中网格化（MESH）模块，对几何体使用六面体单元进行有限元网格自动划分，网格类型为C3D8。在有限元的预处理中，为了简化模型并且获得结构的整体应力分布，没有考虑环氧树脂层，将石墨垫片与反射镜、石墨垫片与中心芯轴的网格采用共节点的方式连接。

2.2随机振动分析

在结构动力学分析中，常采用模态分析的方法来确定和评价结构的振动特性，如确定结构的自然频率、振型、振型的参与系数等。在模态分析中，因振动被假定为自由振动，所以只有边界条件起作用，其它载荷对分析结果没有影响。望远镜结构是通过固定芯轴的方式与其他机械结构相连接固定，所以其边界条件是对芯轴的前后端面的自由度进行全约束[13]。

随机振动分析采用模态叠加法[14]，先对结构进行模态分析，采用ABAQUS软件的Lanczos算法，设置分析步为线性摄动步（liner perturbation）中的频率提取分析步（frequency extraction），模态分析嵌套式望远镜模型。提取结构模态信息，所求的模态数量应该覆盖激励频率范围内，故只提取了前30阶固有频率和振型。表3为模态分析的固有频率数据，列出了前5阶固有频率数值，图3为望远镜结构的前两阶振型，可知前两阶振型反应了望远镜扇区整体的径向摆动，望远镜的各阶振型都主要是反应了扇区的摆动和扭曲。

模态分析是随机振动分析的基础，在完成模态分析后才能进行随机振动分析。随机振动也称为功率谱密度（power spectral density，PSD）分析，本分析采用加速度功率谱作为载荷，加速度功率谱见表4，沿望远镜的轴向作用于芯轴。

在望远镜芯轴的前后端面施加径向（即图3模型中的X方向）如表4的加速度功率谱，确定该方向振动下望远镜结构的应力情况，分析得到图4的应力云图。

根据分析结果，图4给出了望远镜模型轴向随机振动响应，其结构的应力主要集中在望远镜从第1层反射镜到53层的反射镜，最大应力出现在第1层石墨垫片和芯轴的连接处，即图4所示A点位置，在53层60°扇区模块和30°扇区模块的分界位置应力值比内层略小，即图4所示的B点位置。图5为两个位置A和B的应力均方值曲线，在300～400 Hz之间，应力值迅速增大且达到了最大值，即A点位置的最大应力值约为50 MPa，B点位置的最大应力值约为40 MPa，说明此结构的随机响应主要以400 Hz频率之内的频率成分为主。望远镜的反射镜表面镀制X射线薄膜，通常薄膜所能承受的应力值约为十几个兆帕，根据分析得到望远镜模型的应力值，尤其是A、B点位置的应力值远大于反射镜表面的多层膜能承受的应力大小，可能会导致环氧树脂、薄膜、反射镜基底间脱离，破坏望远镜的结构。

上述装配方式的望远镜结构在受到随机激励下，内层反射镜承受应力很大，结构的可靠性难以保证。从图3的振型分析可以看出：望远镜结构若只由内层支撑结构的整体，每一个60°扇区的摆动和扭曲都很大，导致结构的内53层随机振动的应力响应很大。在原有结构基本不变的情况下，在53、54、55层将反射镜位置切换，使得相邻独立的60°扇区连接起来，这三层反射镜称之为切换层，如图6所示。模态分析了切换层方式的望远镜结构，如图7所示。

针对切换层方式的望远镜结构，得到其第一阶和第二阶固有频率分别为353.06 Hz、696.6 Hz，比较分析图3和图7得出：每一个60°扇区的摆动和扭转比没有切换层方式的望远镜结构明显变小，位移的偏移幅值减小，且切换层方式的望远镜结构的固有频率提高了很多。

由图8分析可得，望远镜的应力集中在内4层和切换层的反射镜位置处，应力最大为5 MPa，同样的位置，从图9应力曲线可知，在整个频率内，应力值增长缓慢，且A点应力值从最大值50 MPa减小为5 MPa，B点位置应力最大值从40 MPa减小为2.5 MPa，其应力值都在多层膜、玻璃和环氧树脂承受的应力大小范围内，图5与图9对比分析，说明了切换层方式的望远镜结构可以有效地降低望远镜结构在随机激励下的响应应力，明显地增强了望远镜结构的可靠性。

3结论

本论文基于X射线聚焦望远镜的初始光学结构，利用ABAQUS有限元分析软件，建立有限元模型并完成结构的模态及随机振动响应分析，结论如下：

（1） 对于6个独立的60°扇区和12个30°扇区（无切换层）的望远镜结构，内53层的反射镜和石墨垫片承受的应力较大，应力最大处在最内层，达到了50 MPa，超过了多层膜可以承受的应力范围。

（2） 采用切换层方式的望远镜结构，其应力主要集中在内4层和切换层位置，最大应力约为5 MPa，极大地降低了望远镜结构内层反射镜和分扇区53层反射镜受到的应力，可以满足多层膜、环氧树脂和玻璃的应力强度要求。

参考文献：

[1]WEISSKOPF M C.The advanced xray astrophysics facilityan overview[J].Astrophysical Letters and Communications，1987，26：16.

[2]WEISSKOPF M C，TANANBAUM H D，VAN SPEYBROECK L P，et al.Chandra Xray observatory （CXO）：overview[J].SPIE，2000，4012：216.

TANAKA Y，INOUE H，HOLT S S.The Xray astronomy satellite ASCA[J].Publications of the Astronomical Society of Japan，1994，46：L37L41.

ASCHENBACH B.Inorbit performance of the XMMNewton Xray telescopes：images and spectra[J].SPIE，2002，4496：822.

[5]CRAIG W W，AN H J，BLAEDEL K L，et al.Fabrication of the NuSTAR flight optics[J].SPIE，2011，8147：81470H81470H14.

[6]KOSAKA T，IGARASHI T，AWAKI H，et al.Vibration properties of hard xray telescope on board satellite[J].SPIE，2010，7732：77323G77323G9.

[7]GENDREAU K C，ARZOUMAIAN Z，OKAJIMA T.The neutron star interior composition explore R （NICER）：an explorer mission of opportunity for soft xray timing spectroscopy[J].SPIE，2012：8443138443138.

[8]NERONOV A，BOYARSKY A，IAKUBOVSKYI D，et al.Potential of LOFT telescope for the search of dark matter[J].Physical Review D，2013：arXiz：1312.5178.

[9]BAVDAZ M，RANDO N，WILLE E，et al.ESAled ATHENA/IXO optics development status[J].SPIE，2011，8147：81470C.

CHEN S H，MU B Z，MA S，et al.Design of hard xray focusing telescope with a large fieldofview[J]SPIE，2014，9272：92721R92721R11.

[11]余成武，卢欣.有限元分析在航天器产品设计中的应用[J].空间控制技术与应用，2008，34（4）：2832.

[12]刘宏颖，穆宝忠，王占山.WolterI 型 X 射线天文望远镜的光学设计[J].光学仪器，2012，34（6）：3136.

[13]刘芳，孟光.随机振动载荷下电路板组件三维有限元模拟[J].振动与冲击，2012，31（20）：6164.

[14]张文元.ABAQUS 动力学有限元分析指南[M].中国图书出版社，2005.

（编辑：张磊）