偏置抛物面卫星天线反射器固化变形分析

牛雪娟; 刘江雨; 李辰阳; NIU Xuejuan; LIU Jiangyu; LI Chenyang

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

偏置抛物面卫星天线反射器固化变形分析 PDF

- ORCID：
牛雪娟 ^1,2
✉
- ORCID：
刘江雨 ¹
- ORCID：
李辰阳 ¹

1. 天津工业大学机械工程学院，天津 300387； 2. 天津市现代机电装备技术重点实验室，天津 300387

中图分类号： TB332

最近更新：2022-12-23

DOI：10.12044/j.issn.1007-2330.2022.06.009

摘要

偏置抛物面卫星天线反射器由碳纤维预浸料铺制而成，在固化过程中通常会出现热变形、固化收缩变形及机械变形等，且过大的变形会导致反射器反射信号产生偏差。为了有效预测天线反射器在热压罐固化过程中的变形，本文采用线性黏弹性本构模型表征纤维增强复合材料在固化过程中的应力-应变关系。采用顺序热-力耦合分析法，对天线反射器内部温度场、固化度场及位移场的分布进行建模，并对铺层结构及降温速率对天线反射器变形的影响开展研究。结果表明：相较于纯0°铺层结构，正交和准各向同性铺层天线反射器z向变形量减小了51%以上，且降低降温速率有助于提高天线反射器的固化型面变形。

关键词

天线反射器; 碳纤维预浸料; 固化变形; 顺序热-力耦合

0 引言

偏置抛物面卫星天线反射器的工作原理是由一点发射信号通过反射面进行反射，形成平行信号射向接受卫星，制备过程中过大的变形会使天线反射器的信号产生偏转，因此需要对天线反射器的变形进行合理预测并制定相应的优化方案。

在碳纤维制品的固化变形研究中，相较于试验分析仿真分析因其灵活性高、时间短、成本低等优点被广泛应用，目前国内外已有不少针对固化变形的仿真分析研究。K.K.YEONG等^［

1］试验测定了3501-6树脂的一系列黏弹性行为数据，为复合材料黏弹性本构模型研究奠定了良好的基础。闵荣等^{［参考文献 2

百度学术}2］根据黏弹性本构模型对C型的复合材料结构件进行了仿真与实验对比，验证了黏弹性本构模型预测的有效性。X.Y.LI等^{［参考文献 3

百度学术}3］研究了在热传导过程中不同材质的模具对复合材料结构件温度梯度的影响。元振毅等^{［参考文献 4

百度学术}4］通过引入层间剪切层的概念，模拟了模具与构件之间相互作用对变形的影响，使固化变形预测更加符合实际。祝君军等^{［参考文献 5

百度学术}5］对比了阴模模具和阳模模具对V型复合材料结构件固化变形的影响，得出使用阳模模具会比使用阴模模具变形要大的结论。而针对天线反射器的研究如下：邢思远等^{［参考文献 6

百度学术}6］模拟了偏置抛物面卫星天线反射器在服役过程中不同工况温度下的热变形；江文剑^{［参考文献 7

百度学术}7］建立了反射器模具的补偿公式，使模具在反射器固化成型的时候恰好变形到理论型面。

上述研究中的模型大多为平板型、L型以及C型等简单形状结构件的固化变形。对于偏置抛物面天线反射器的研究也多为服役过程中的热变形研究，而针对天线反射器在制造过程中固化变形的研究较少。本文采用顺序热-力耦合的方式模拟天线反射器的固化变形，通过施加不同的边界条件模拟了天线反射器在热压罐中加热固化和脱模之后两种状态，考虑了脱模之后残余应力的释放对变形的影响，为天线反射器变形的预测以及试验提供了可靠的依据。

1 理论模型

1.1 反射器的热传导-固化模型

碳纤维复合材料在固化过程中的二阶偏微分热传导控制方程表达式^［

8］为：

λ_{x} \frac{\partial^{2} T}{\partial x^{2}} + λ_{y} \frac{\partial^{2} T}{\partial y^{2}} + λ_{z} \frac{\partial^{2} T}{\partial z^{2}} + Q = ρ_{c} c_{c} \frac{\partial T}{\partial t}

(1)

式中，λ_x，λ_y，λ_z分别为复合材料在x，y，z方向上的热传导系数；T为固化工艺温度； $ρ_{c}$ 为密度； $c_{c}$ 为比热容；Q为t时刻复合材料内部产热量，其表达式为：

Q = ρ_{r} V_{r} H_{r} \frac{d α}{d t}

(2)

式中， $ρ_{r}$ 为树脂密度； $V_{r}$ 为树脂体积分数； $H_{r}$ 为单位质量树脂的产热量； $α$ 为固化度； $\frac{d α}{d t}$ 为固化速率。

1.2 反射器固化过程中的应力-应变模型

树脂为高分子聚合物材料，在固化过程中表现出明显的黏弹性行为，在目前的研究中通常采用广义Maxwell模型^［

9］描述树脂材料的力学松弛行为。广义Maxwell模型如图1所示，广义Maxwell模型由n个Maxwell单元和一个弹性模量为

E_{\infty}

的弹簧并联组成。其中，n个Maxwell单元模拟了聚合物运动单元的多重性；最后一个弹簧为平衡模量，表征聚合物材料由于受到网络间交联点的限制，使应力松弛至某一平衡值就不再变化这一特征。对于广义Maxwell模型，在恒定初始应变作用下，其弹性模量表达式为：

E (t) = E_{\infty} + \sum_{m = 1}^{n} w_{m} (E_{0} - E_{\infty}) e^{- \frac{t}{τ_{m}}}

(3)

式中， $E_{0}$ 为广义Maxwell模型未发生松弛时的模量； $E_{\infty}$ 为广义Maxwell模型完全松弛时的模量； $w_{m}$ 为第m个Maxwell单元的权重系数； $τ_{m} = η_{m} / E_{m}$ 为第m个Maxwell单元的松弛时间。

图1 广义Maxwell模型

Fig.1 Generalized Maxwell model

在固化过程中，由于温度与固化度会随时间变化，因此需要考虑温度与固化度对材料黏弹性行为的影响。线性黏弹性材料的本构方程^［

2，10］表示为：

σ_{i} (t) = \int_{0}^{t} C_{i j} (a_{c}, T_{c}, ξ^{t} - ξ^{'}) \frac{\partial ε_{j}^{e f f}}{\partial τ} d τ

(4)

式中， $ξ^{t}$ 和 $ξ^{'}$ 为折算时间； $a_{c}$ ， $T_{c}$ 分别为参考固化度和参考温度； $ε_{j}^{e f f}$ 为有效应变增量；C_ij为刚度矩阵，其表达式为：

C_{i j} = C_{i j}^{\infty} + \sum_{m = 1}^{n - 1} W_{m} (C_{i j}^{0} - C_{i j}^{\infty}) e^{- \frac{ξ}{τ_{m}}}

(5)

由于公式（4）不能直接在有限元软件中执行，因此将其转换为应力增量表达形式：

Δ σ_{i}^{t + Δ t} = σ_{i}^{t + Δ t} - σ_{i}^{t} = C_{i j}^{*} Δ ε_{j}^{t + Δ t} + Δ σ_{i}^{R}

(6)

式中，

C_{i j}^{*} = C_{i j}^{\infty} + \sum_{m = 1}^{n - 1} {\bar{C}}_{i j} W_{m} \frac{a_{T} τ_{m}}{Δ t} (1 - e^{- \frac{Δ ξ^{t + Δ t}}{τ_{m}} - 1})

(7)

Δ σ_{i}^{R} = \sum_{m = 1}^{n - 1} S_{i m}^{t} e^{- \frac{Δ ξ^{t + Δ t}}{τ_{m}} - 1}

(8)

式中， $S_{i m}^{t}$ 为历史状态变量，其表达式为：

S_{i m}^{t} = e^{- \frac{Δ t}{τ_{m}}} S_{i m}^{t - Δ t} + {\bar{C}}_{i j} W_{m} \frac{a_{T} τ_{m} Δ ε_{j}^{t}}{Δ t} (1 - e^{- \frac{Δ ξ^{t}}{τ_{m}}})

(9)

2 天线反射器

2.1 有限元模型建立

本文所用的偏置抛物面卫星天线反射器模型^［

11］由一抛物线绕旋转轴一周所形成的抛物面与一圆柱面相截而成。抛物面的母线方程为

Y ² = 4 \times F \times X

，焦点坐标为（100，0）。圆柱面直径为350 mm，与回旋抛物面中心的偏置距离为219.76 mm。所截的反射器长轴为520 mm，短轴为350 mm。对天线反射器进行有限元建模时，将天线反射器模型导入ABAQUS前处理模块，划分网格后创建为网格部件。然后将网格部件进行实体网格偏移，每层厚度为0.125 mm，共16层，28 320个单元。模型中使用离散坐标系定义每一层的材料方向，其中材料法向轴方向垂直于模型中各个单元，主轴方向定义为x轴方向。主轴沿法向轴旋转不同角度即可得到不同的铺层方向。图2为反射器几何尺寸及45°铺层结构的网格模型示意图，其中x’为铺层在该点处的纤维方向，z’为其局部曲面法向。

图2 反射器几何尺寸及45°铺层网格模型

Fig.2 Reflector geometry and 45°grid model

2.2 固化变形相关的用户子程序设计

在进行天线反射器固化变形仿真时，温度边界条件、材料内部生热及黏弹性本构模型都具有时变特性。仅靠Abaqus软件自身并不能直接定义相关参数，需用Fortran语言编写对应的用户子程序来实现。

在进行热传导-固化分析时，使用DISP子程序定义固化工艺温度。使用FILM子程序定义对流换热边界条件，主要定义对流换热系数和环境温度SINK值。其中，对流换热系数H（1）=38 W/（m²/K）^［

4］，SINK值与工艺温度值一致。使用USDFLD子程序定义固化度，并通过状态变量STATEV与定义内部生热的子程序HETVAl交换数据。使用HETVAL子程序定义固化动力学方程中的固化速率、固化度和材料的内部生热。

天线反射器在热压罐中固化成型时，随着温度载荷升高固化度逐步增加，复合材料力学本构关系不断变化，具有典型的黏弹性特性。本文依据式（6）~式（9），通过UMAT子程序为不断变化的本构模型提供应力增量更新和雅可比矩阵 $\partial σ / \partial ε$ 。并使用UEXPAN子程序定义复合材料在3个方向上的化学收缩应变和热应变特性。本文中预浸料型号为AS4/3501-6，其在未发生松弛时的材料参数见表1。

表1 AS4/3501-6预浸料的未松弛材料参数^{［参考文献 2

百度学术}2，12］

Tab.1 Unrelaxed material parameters of AS4/3501-6 composite ^{［参考文献 2

百度学术}2，12］

弹性模量/GPa			剪切模量/GPa			泊松比			线胀系数/10^-6℃^-1			化学收缩系数/10^-6℃^-1
纵向	横向		面内	横向		面内	横向		纵向	横向		纵向	横向
126	8.30		4.10	2.80		0.25	0.35		0.50	35.30		-167	-8 810

2.3 基于顺序热-力耦合的固化变形模拟

本文采用顺序热-力耦合的方法模拟天线反射器的固化变形，分别通过热传导-固化模块和固化变形两个模块进行分析。本文先以非对称铺层方式（［0°/45°/-45°/90°］₄）为例进行分析，材料选用研究较为完善的AS4/3501-6预浸料，关于3501-6树脂固化动力学所需的各参数值请参考文献［

2］。在热传导-固化分析时使用热传递分析步，单元类型为DC3D8，固化工艺温度变化见图3。起始温度为293 K，以2.5 K/min的速率升温至393 K后保温60 min；然后以3.8 K/min的速率升温至450 K后保温120 min；最后以-2.42 K/min的速率降温至293 K。

图3 固化工艺中温度曲线

Fig.3 Temperature curve in curing process

数值模拟时，对天线反射器下表面施加温度载荷，其余面定义热对流交换条件，所施加的边界条件见图4。在分析固化变形模块时，共设置两个静力通用分析步：一个用于分析步反射器在固化过程中的变形；另一个用于模拟脱模后残余应力释放引起的变形。由于天线反射器为偏置抛物面结构，曲面的变曲率特征使其在固化过程中相对于模具的位移量小，因此，在脱模之前将天线反射器的底面固定，上表面施加0.6 MPa的压力。脱模之前的边界条件见图5，脱模之后固定模型下表面长轴与短轴交点。

图4 热传导-固化分析步边界条件

Fig.4 Thermal conduction-curing boundary conditions

图5 脱模之前边界条件

Fig.5 Boundary conditions before demolding

3 结果分析

3.1 天线反射器温度与固化度分析

热传导-固化模块求解完成后，得到在固化过程中任意时刻反射器的温度、固化度分布。第1 013 s的天线反射器温度场分布云图如图6所示。可知，天线反射器各处的温度分布较为均匀，没有出现明显的温度梯度。图7为天线反射器模型中心点的温度与固化度随时间变化的曲线。从固化度曲线中可以看出，固化度在两个升温阶段固化速率较大，在达到第二个保温阶段时天线反射器的固化基本完成。由温度曲线可知，由于反射器厚度较小，反射器内部的温度与所施加的温度载荷基本一致。

图6 温度场分布图（t=1 013 s）

Fig.6 Temperature field distribution map （t=1 013 s）

图7 固化度与温度曲线

Fig.7 Curing degree and temperature curve

3.2 天线反射器固化变形分析

在分析反射器固化变形模块时，反射器模型的节点编号、网格数量需要与热传导-固化模块的反射器模型完全一致，这样便于将热传导-固化模块得到的温度场后处理结果odb文件作为温度载荷导入固化变形模块中的预定义场。固化变形模块求解后得出天线反射器脱模前后z方向位移云图如图8所示。可知，反射器在脱模之前由于受到压强载荷及固化收缩的影响，使反射器在z轴产生压缩变形；在脱模之后反射器出现翘曲变形，且沿反射器中心成对角分布，最大变形量为2.3 mm左右。可以看出，非对称铺层结构的天线反射器的固化变形量较大。

图8 反射器z方向变形图

Fig.8 Reflector z-direction deformation diagram

（a）　脱模前　　　　　　　（b）　脱模后

图9展示了天线反射器上表面短轴端点A在z方向应变值随时间的变化。在第Ⅰ阶段，树脂处于流动状态，反射器受热膨胀使应变到达5 mm左右。第Ⅱ阶段处于保温阶段，由于树脂的固化收缩，使应变值减小到1 mm左右。第Ⅲ阶段中，热膨胀与固化收缩共同作用，但是固化收缩在这一阶段作用更为显著。因此使应变曲线呈现出先略微增大然后迅速减小的趋势。在第Ⅳ阶段中，降温使反射器产生热收缩应变，应变值进一步迅速降低。最后，在脱模之后反射器底面的约束解除，残余应力释放，使应变值产生阶跃性变化，z方向最终的应变值为-12 mm。本文研究得出的脱模之前应变曲线走势与文献［

13］中基本一致，并且本文考虑了脱模之后残余应力释放对变形的影响。综上，本文所建立的天线反射器有限元模型具有较高的可靠性。

图9 反射器上表面短轴端点A在z方向应变曲线

Fig.9 Strain curve in z-direction at point A on the upper surface of the reflector

3.3 影响变形的因素分析

为了分析不同铺层结构和固化工艺对天线反射器变形的影响，本文对纯0°铺层（［0°］₁₆）、正交铺层（［0°/90°］_4s）和准各向同性铺层（［0°/45°/-45°/90°］_2s）三种对称铺层结构进行对比分析。不同铺层方式的天线反射器脱模之后的位移如图10所示。

（a）［0°］₁₆铺层

（b）［0°/90°］_4s铺层

（c）［0°/45°/-45°/90°］_2s铺层

图10 不同铺层方式在x、y、z方向上的变形图

Fig.10 Deformation diagram in x， y and z directions for different layering methods

由图10可得，纯0°铺层的天线反射器变形量在x方向上最小。但在y方向上，纯0°铺层的反射器相较于正交对称铺层的反射器发生了严重的固化收缩变形，最大变形量在0.7 mm左右。在z方向上，纯0°铺层的反射器变形主要发生在短轴两侧，正交对称铺层的反射器变形位置主要发生在长、短轴端点附近。两种正交对称铺层结构中，准各向同性铺层的反射器变形量相对较小，但是由于±45°铺层的原因，使z方向最大变形位置沿±45°方向发生了偏移。3种不同的铺层结构中，正交对称铺层结构比纯0°铺层结构在z方向上的最大变形减小了51%以上。造成这种现象的原因是固化过程中构件会产生内应力，而正交铺层结构在厚度方向上对称，层间产生的应力能够被部分抵消。

为分析降温速率对变形的影响，以正交铺层天线反射器为例，采用3种降温速率进行对比。图11为3种降温速率下的反射器上表面短轴端点A在z方向的应变曲线。可知，当降温速率分别为-2.4、-2、-1 K/min时，对应的应变值依次减小，所以适当降低降温速率有助于减小天线反射器的固化变形。

图11 不同降温速率下z方向应变曲线图

Fig.11 Strain curve in z-direction under different cooling rates

4 结论

（1）本文将整个制备工艺分为脱模前、脱模后两个工序，并分别设置不同的边界条件和载荷类型，模拟了天线反射器在固化过程中及脱模后的变形行为，从而使天线反射器的变形预测更加符合实际。

（2）对比分析不同铺层结构对固化变形的影响得出，采用正交对称铺层结构的天线反射器在z方向上变形量比纯0°铺层结构的天线反射器变形量减小了51%以上。

（3）对比3种不同的降温速率得出，降温速率越小，z方向上的应变值越小。因此，降低降温速率有助于减小天线反射器的固化变形。

参考文献

YEONG K K， SCOTT R W. Stress relaxation behavior of 3501-6 epoxy resin during cure［J］. Polymer Engineering & Science，1996，36（23）： 2852-2862. [百度学术]

闵荣，元振毅，王永军，等.基于黏弹性本构模型的热固性树脂基复合材料固化变形数值仿真模型［J］.复合材料学报，2017，34（10）：2254-2262. [百度学术]

MIN Rong，YUAN Zhenyi，WANG Yongjun，et al. Numerical simulation for curing deformation of resin matrix thermosetting composites using viscoelastic constitutive model［J］. Acta Materiae Compositae Sinica，2017，34（10）：2254-2262. [百度学术]

LI Xiaoyang， WANG Jihui， LI Shuxin， et al. Cure-induced temperature gradient in laminated composite plate： Numerical simulation and experimental measurement［J］. Composite Structures，2020，253（6）：112822. [百度学术]

元振毅，王永军，王俊彪，等.基于模具-制件相互作用的复合材料制件固化变形数值模型［J］.复合材料学报，2016，33（04）：902-909. [百度学术]

YUAN Zhenyi，WANG Yongjun，WANG Junbiao，et al. Numerical model on curing deformation of composite part based on tool-part interaction［J］. Acta Materiae Compositae Sinica， 2016，33（04）： 902-909. [百度学术]

祝君军，文琼华，罗辑，等.模具形式对V型结构复合材料固化变形的影响［J］.宇航材料工艺，2018，48（05）：44-48. [百度学术]

ZHU Junjun， WEN Qionghua， LUO Ji，et al.Effect of mold forming method on curing deformation of V-shaped composites［J］.Aerospace Materials and Technology，2018，48（05）：44-48. [百度学术]

邢思远，刘洪新，彭志刚，等.高精度复合材料天线反射器结构与热变形仿真优化分析［J］.玻璃钢/复合材料，2018（08）：66-73. [百度学术]

XING Siyuan，LIU Hongxin，PENG Zhigang，et al.Simulation and optimization structure and thermal deformation of high precision composite antenna reflector［J］.Fiber Reinforced Plastics /Composites，2018（08）：66-73. [百度学术]

江文剑，史明，杨亚宁.一种高精确度偏置反射面天线模具设计［J］.太赫兹科学与电子信息学报，2018，16（04）：653-657. [百度学术]

JIANG Wenjian， SHI Ming， YANG Yaning.A high-precision mould design method of offset-reflector antenna［J］. Journal of Terahertz Science & Electronic Information Technology，2018，16（04）：653-657. [百度学术]

BOGETTI T A， GILLESPIE J W.Process-induced stress and deformation in thick-section thermoset composite laminates［J］.Journal of Composite Materials，1992，26（5）： 626-660. [百度学术]

王子川，林福严，刘培伦，等.根据松弛模量求近似松弛谱的方法研究［J］.煤炭技术，2021，40（01）：156-158. [百度学术]

WANG Zichuan， LIN Fuyan， LIU Peilun，et al.Research on solving method of approximate spectrum based on relaxation modulus［J］.Coal Technology，2021，40（01）：156-158. [百度学术]

乔巍，姚卫星，黄杰.考虑模具和黏弹性的复合材料固化变形数值模型［J］.复合材料科学与工程，2020（10）：67-74. [百度学术]

QIAO Wei， YAO Weixing， HUANG Jie.Numerical model of process-induced deformation in composite parts considering mould and material viscoelasticity［J］.Composites Science & Engineering，2020（10）：67-74. [百度学术]

陆炳荣.偏置抛物面的几个几何设计计算问题［J］.电子机械工程，1999（01）：50-53. [百度学术]

LU Bingrong. Some problems to be solved in geometric design and calculation of the offset paraboloid［J］.Electronic Machinery Engineering，1999（01）：50-53. [百度学术]

SCOTT R. White，YEONG K. Kim. Process-induced stress relaxation analysis of AS4/3501-6 laminate［J］.Journal of Reinforced Plastics and Composites，1997，16（1）：2-16. [百度学术]

刘哲.热固性树脂基复合材料固化变形数值模拟分析［D］.哈尔滨：哈尔滨工业大学，2014：45-46. [百度学术]

LIU Zhe. Numerical simulation of curing deformation of thermosetting resin matrix composites［D］. Harbin：Harbin Institute of Technology，2014：45-46. [百度学术]

You are the visitor

Beijing Public Network Security：11010602060147号京ICP备06019415号-2

Address：P.O.Box 9200-73, Beijing 100076, P.R.China

Postcode：100076 Fax：010-68383237

ServiceTel：010-68383269 E-mail：703@china.com

About Journal

Editorial Board

Guide for Authors