露天大型射电望远镜在日常运行中，不可避免地受到周围环境的影响，如太阳辐射的作用，这种影响随着太阳位置的变化而变化。天线钢结构表面温度载荷常成为结构设计的控制载荷，在天线结构设计中需要掌握日照作用下最不利的温度场分布和热变形情况；同样在天线的日常运行过程中，了解天线结构特别是面板的热分布和热变形对天线的高效运行有重要的意义，而影响日照条件下温度的因素很多，如太阳辐射、与大气的对流和天线自身的遮挡等。为此，借鉴文^[1]关于上海65 m射电望远镜太阳辐射作用的研究成果，并根据新疆南山站25 m天线的具体情况，以ASHRAE晴空模型为基础，结合北京晴空模型系数回归方程^[2]，研究了新疆南山站的太阳辐射情况，并借助有限元分析软件ANSYS对南山25 m天线面板的热特性进行了模拟分析，得到天线面板的温度场分布和热变形。

依据ASHRAE晴空模型理论，结合北京晴空模型中关于太阳辐射强度拟合得到的随日序数变化的系数^[2]，可以计算瞬时太阳辐射强度，根据此模型，天线接收的太阳辐射由直接辐射、天空散射辐射和地面反射辐射组成^[3]。

(1)直射辐射强度

对于平面上的构件，其垂直入射的直射辐射强度：

(2)散射辐射强度

水平面上的构件接收的散射辐射强度：

(3)反射辐射强度

反射到构件表面上的辐射量：

综上，入射到某一非垂直构件表面的太阳总辐射为

热对流，指流体各部分由于发生相对位移而引起的热量交换现象。构件表面的对流换热可参考牛顿换热定律^[3]，即

构件表面得到的净长波辐射可以表示为

如图 1，太阳光线入射到天线主反射面，记X₀=[x₀,y₀,z₀]^T为天线反射面最外圈边沿上的一点，X_d=[x_d,y_d,z_d]^T为太阳光线入射方向(由太阳的位置决定)，并设X=[x,y,z]^T为此入射光线上的一点，求得该入射光线和反射面的交点X′即可知道天线面板的遮挡情况。

图 1 太阳入射到天线示意图 Fig. 1 Illustration of incidence of sunlight onto the dish surface of a radio antenna

图选项

由起点发出的光线方程的矩阵表示为

考虑方程(16)在Δ>0的情况下，可以求得入射光线与反射面之间的阴影分界线。利用以上讨论的算法，对南山25 m天线主反射面在夏至日的阴影情况进行了分析和计算，结果如图 2，图中用两种颜色表示的范围为阴影区，可以看到从7点到10点以及19点到21点有不同程度的阴影。通过对面板自身遮挡的讨论，可以更加详细地讨论太阳辐射对天线面板的影响。

图 2 天线面板遮挡情况 Fig. 2 Sunlight coverage patterns on the dish surface under the influence of the screening of the dish itself

图选项

完成了25 m天线整体有限元模型的建立，主要包括座架体、反射面支架、反射面板等，其中面板采用板壳单元，支架采用梁单元，座架结构采用实体梁单元，材料参数包括设置密度、弹性模量、泊松比、热膨胀系数等^[6]。如图 3左图是25 m天线整体结构有限元模型，其中有636个节点，600个壳单元，1 734个梁单元以及232个质量块单元。但本文重点讨论天线面板热特性，其有限元模型如图 3右图，面板主要由8圈组成，其中内3圈由中心体支撑，内圈最少为16块面板，外5圈由32条辐射梁支撑，外圈最多为64块面板，总共有352块面板，面板几何排列如图 3。在建模中考虑面板采用铝合金材料，其厚度假定为3 mm，梁和其他单元使用钢材料，属性如表 1。

图 3 天线有限元模型 Fig. 3 A finite-element model of the antenna

图选项

首先讨论了南山6月21日太阳的位置如图 4，太阳在7点左右升起，升起方位为70°左右，在14点左右太阳高度角达到最大，在21点之后落下，落下时方位角为300°左右。

图 4 夏至日太阳位置图 Fig. 4 Directions to the Sun in June 21st

图选项

同时也分析了352块面板的等效方位角和倾斜角，然后结合上文讨论的方法，可以算得该天各块面板在不同时刻受到的太阳辐射热流密度，作为热载荷加载到有限元模型上，设置相关的参数，并以7点的热流和温度作为初始值。根据长期观测数据，取初始温度为7 ℃，每隔1小时进行一次热流密度加载，在ANSYS中进行瞬态热分析，最后得到天线面板上从太阳升起到落下这段时间内的温度分布，图 5是从早上8点到晚上19点之间天线面板的温度分布图。

图 5 天线面板温度场分布 Fig. 5 Changing temperature distributions on the dish surface in June 21st

图选项

可以看到天线面板上的温度分布，呈现随太阳位置的变化而变化的规律。在8点钟各块面板的温度基本保持在8.2 ℃左右，接近初始温度；然后随着太阳的升起，面板上温度值逐渐增大，温度梯度和温度分布的位置也随着发生变化；在15点左右面板局部温度达到最大值45.3 ℃。

图 6是各个时刻瞬态温度分布曲线，图中所示一天中各个时刻温度最大值从8点到15点一直处于上升状态，在15点达到最大值45.3 ℃，主要是在这段时间内接收的太阳辐射一直在增强，之后到21点这个时间段内逐渐下降。各时刻温差的最大值可达12 ℃多，同样出现在15点左右，因此天线面板在该时刻左右的变形最为严重，对天线指向的影响最大；在整个太阳照射的时间段内面板上的最大温差变化不是很大。

图 6 一天中各个时刻平均温度变化 Fig. 6 Hourly averaged temperature changes in a day

图选项

针对大型天线温度效应的研究是一个极其复杂的课题，本文针对面板日照温度特性研究主要以面板内表面为主，只考虑了天线反射面内表面的日照效应，没有考虑外表面，但是由于面板厚度比较小，产生的影响不是很明显，这方面工作有待深入研究。文^[7]就国外相关领域的工作做过一些综述性的研究，证实了实测和仿真相结合的研究方法的有效性，通过实测数据验证仿真的工作是完全可行的，限于本文主要工作是仿真分析，实测方面有待以后工作中不断完善。

本文主要利用常用的ASHRAE晴空模型，并结合北京晴空模型的系数回归方程，对新疆南山25 m天线受太阳辐射情况进行了讨论，结合天线面板自身遮挡的算法，对夏至日这天照射到天线各个面板上的热流密度进行了讨论和计算；借助有限元分析软件建立了25 m天线的有限元模型并对天线面板内反射面的温度分布做了详细的分析，可以看到以铝合金为材料的面板，其热变化和热分布随着太阳的照射有较大的变化，该天线在15点左右温度分布达到最大值，此时会有较大的面板变形，对天线的指向影响最大。后续工作中将结合温度分布讨论面板的变形情况，以及具体地讨论对天线指向造成的影响，同时结合实测对仿真结果进行验证，从而更好地了解太阳辐射对天线的影响。