CN111274726B - 一种考虑热效应的天线罩电磁性能分析方法 - Google Patents
一种考虑热效应的天线罩电磁性能分析方法 Download PDFInfo
- Publication number
- CN111274726B CN111274726B CN202010047861.0A CN202010047861A CN111274726B CN 111274726 B CN111274726 B CN 111274726B CN 202010047861 A CN202010047861 A CN 202010047861A CN 111274726 B CN111274726 B CN 111274726B
- Authority
- CN
- China
- Prior art keywords
- finite element
- electromagnetic
- element model
- radome
- medium
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 230000000694 effects Effects 0.000 title claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 50
- 230000005540 biological transmission Effects 0.000 claims abstract description 28
- 238000006073 displacement reaction Methods 0.000 claims abstract description 27
- 239000000463 material Substances 0.000 claims abstract description 23
- 230000003068 static effect Effects 0.000 claims abstract description 9
- 230000035699 permeability Effects 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000005284 excitation Effects 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 9
- 230000009471 action Effects 0.000 claims description 8
- 238000000034 method Methods 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 5
- 239000003989 dielectric material Substances 0.000 claims description 4
- 238000012805 post-processing Methods 0.000 claims description 4
- 229920001971 elastomer Polymers 0.000 claims description 3
- 239000000806 elastomer Substances 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 2
- 230000010287 polarization Effects 0.000 claims description 2
- 230000003938 response to stress Effects 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000005672 electromagnetic field Effects 0.000 description 3
- 238000012938 design process Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005405 multipole Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Landscapes
- Details Of Aerials (AREA)
Abstract
本发明提出了一种考虑热效应的天线罩电磁性能分析方法，该方法主要包括如下步骤：S1：建立天线罩力学计算电磁学计算有限元模型；S2：对S1建立的力学有限元模型添加材料属性、施加温度载荷、边界条件；S3：完成对静力学问题的求解，并提取节点位移；S4：基于S3提取的节点位移重构温度载荷下变形的有限元模型；S5：基于HyperMesh软件对S4中变形的有限元模型进行重构得到电磁学计算有限元模型；S6：完成对于S5中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成温度载荷下天线罩透波性能分析；S7：完成对于S1.2中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成原始状态天线罩透波性能分析。
Description
技术领域
本发明属于计算电磁学结构强度交叉领域，尤其涉及一种考虑热效应的天线罩电磁性能分析方法。
背景技术
随着航空航天技术不断朝着高速化、轻量化发展，各类飞行器面临的飞行工况也越来越复杂。高速飞行的飞行器承担的载荷主要体现在高速飞行所带来复杂的温度载荷。天线罩位于整个飞行器的最前端，直接承担高速飞行的飞行器所带来的温度载荷。同时，天线罩又承担着对于整个飞行器雷达制导等关键电子元器件的保护功能。但是，在飞行器实际的设计过程中，天线罩的强度设计和电磁性能设计往往是两个团队，这就导致在实际的设计环节中出现迭代周期长、迭代成本高等一系列的问题。因此，对天线罩联合温度场和电磁场进行计算和仿真就显得极其重要。在各类文献中，针对电大尺寸天线罩考虑热、电场耦合的分析暂还存在空白，几乎没有公开文献对这一块进行过专业的探讨和系列的研究。因此，针对服役在复杂温度场下的天线罩，开展热电联合仿真计算就显得极其重要，这也成为了现在飞行器设计过程中的一个急需解决的技术难题。
发明内容
发明目的：解决目前在天线罩研发过程研发周期长、迭代成本高、力学和电磁学之间没有交叉学科导致天线罩电性能和力学性能无法同时优化等问题，降低天线罩的设计周期和设计成本。
技术方案：为实现本发明的目的，本发明所采用的技术方案是：一种考虑热效应的天线罩电磁性能分析方法，包括以下步骤：
S1：建立天线罩力学计算有限元模型、电磁学计算未变形有限元模型；
S2：对S1建立的力学有限元模型添加材料属性、施加温度载荷、边界条件；
S3：完成对静力学问题的求解，对计算结果进行后处理，提取节点位移；
S4：基于Python语言完重构S3有限元模型，得到温度载荷下变形的有限元模型；
S5：基于HyperMesh软件对S4中变形的有限元模型进行重构得到电磁学计算有限元模型；
S6：完成对于S5中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成温度载荷下天线罩透波性能分析；
S7：完成对于S1中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成原始状态天线罩透波性能分析，并和S6的分析结果进行对比，获得气动热载荷对于天线罩透波性能的影响。
更进一步地，所述步骤S1建立天线罩力学计算有限元模型、电磁学计算未变形有限元模型，具体包括：
S1.1：基于HyperMesh软件为天线罩的几何结构划分六面体单元有限元网格，初步建立力学计算有限元模型。
S1.2：基于HyperMesh软件将天线罩几何结构的外表面划分为三角形单元有限元网格，建立电磁学计算未变形有限元模型。
S1.3：对步骤S1.1建立的有限元模型设置材料属性。
更进一步地，所述步骤S2对S1建立的力学有限元模型添加材料属性、施加温度载荷、边界条件，具体包括：
S2.1：添加材料属性，并将天线罩风洞试验实验数据中天线罩头部和尾部的实际温度施加于S1.1建立的有限元模型上；
S2.2：将上述所有的模型输出为inp文件，基于Abaqus求解器，对S2.1的有限元模型进行传热分析，获得其所有节点的稳态节点温度；
S2.3：将S2.2中的稳态节点温度作为温度载荷施加于S1.1建立的有限元模型；
S2.4：将S2.3中建立的有限元模型的尾部施加固支的边界条件。
更进一步地，所述步骤S3完成对静力学问题的求解，对计算结果进行后处理，提取节点位移，具体包括：
S3.1：设置传热分析求解类型、求解时间、输出变量和输出参数；
S3.2：将S3.1建立的模型输出为inp文件进行静力学求解，求解方式如下：
为了求得天线罩弹性单元的温度应力，首先要导出热弹性力学的基本方程和边界条件，命弹性体的各点的变温ΔT＝T2-T1，即是后一瞬时的温度减去前一瞬时的温度，ΔT升温时为正，降温时为负，弹性体内各个点的微小长度不受到约束，将发生正应变αΔT，其中，α是天线罩弹性单元弹性体的线热膨胀系数，ΔT为温度变化量，在各向同性材料中，系数α不随方向的改变而改变，随意这种正应变的各个方向都是相同的，因而也就没有剪应变的产生，在α不随温度的变化发生变化时，可以得到天线罩弹性单元内各个点的应变分量：εx＝εy＝εz＝αΔT，γxy＝γzy＝γzx＝0，εx,εy,εz分别是变形体在三维空间中X，Y，Z的正应变，γxy＝γzy＝γzx分别是变形体在三维空间中XY，ZY，ZX平面的剪切应变；
天线罩弹性单元受到外在约束以及内部的各个点约束时并不能发生上述的线性应变，就会产生温度应力，因此可以得到式(1)：
σx,σy,σz分别是变形体在三维空间中X，Y，Z的正应力，τxy,τyz,τxz分别是变形体在三维空间中XY，ZY，ZX平面的剪切应力，E为弹性模量，μ为泊松比；
在只有变温ΔT存在的情况下，温度场的位移和应力响应可以等效为平面问题：
σx＝0,τyz＝0,τzx＝0
因此，式(1)可以表示为：
由上式可以求得温度场下的天线罩弹性单元应变分量和由变温ΔT表示的应力分量的物理方程：
可以得到相似的天线罩弹性单元平面应变问题的热弹性力学的物理方程：
弹性力学平面问题的几何方程和平衡微方程为：
u和v分别是平面内X方向和Y方向位移；
平衡微方程给出的是天线罩弹性单元体内的应力和外力之间的平衡关系，不会由于引起应力的原因不同而改变平衡状态。
求解温度应力时，一般都宜于按照位移求解，推导按照位移求解平面变温度应力的基本途径如下：
按照平面应变问题进行公式等效可以得到：
将上式导入平衡微方程，化简后可以得到：
由无面力应力条件边界条件可知求解应力平面问题的应力边界条件：
l，m为边界外法线方向和X、Y方向的夹角的余弦，位移边界条件仍然为：
us，
vs，
为已知边界平面上X方向和Y方向位移；
将上述位移坐标转化为全局坐标系，即可得到结构在X、Y、Z方向上的位移：ux uyuz。
更进一步地，所述步骤S4基于Python语言完重构S3有限元模型，得到温度载荷下变形的有限元模型，具体包括：
S4.1：基于Python脚本读取S3.2生成的inp文件，获得节点的X、Y、Z坐标值；
S4.2：将S4.1读取的节点坐标值与S3.2获取的X、Y、Z方向上的位移ux uy uz位移值进行算数求和；
S4.3：将S4.2求和之后的节点坐标重新写入S3.2生成的inp文件。
更进一步地，所述步骤S5基于HyperMesh软件对S4中变形的有限元模型进行重构得到电磁学计算有限元模型，具体包括：
S5.1：将S4.3生成的inp文件重新导入HyperMesh软件，通过软件提取天线罩模型的外表面有限元模型；
S5.2：将S5.1提取的外表面有限元模型的类型改为三角形，有限元模型三角形单元周长改为波源波长的五分之一，进行重新绘制，输出新的inp文件。
更进一步地，所述步骤S6完成对于S5中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成温度载荷下天线罩透波性能分析，具体包括：
S6.1：将S5.2生成的inp文件导入FEKO中，设置电磁材料参数；
S6.2：建立垂直极化的振子天线模型，并设置激励方式和大小；
S6.3：设置远场参数、求解截断误差；
S6.4：为获得温度场下天线罩的透波性能，将S6.3生成的模型文件按照下式的计算方式进行计算：
等效传输线理论相比于传统的传输线理论针对不同的材料板的厚度和层数有十分出色的通用性。由不同层数的介质板依次叠组成的天线罩模型相当于是不同的传输矩阵的级联，如图2所示，θi为天线罩内天线产生的电磁波在罩内入射时的入射角，di是第i层介质的厚度，ε0,εi,μ0,μi分别是在空气或真空中介电常数、介质材料的相对介电常数、空气或真空中磁导率、介质材料的相对磁导率，数学表达式为：εi＝εr(1-jtanδ),μi＝μr(1-jtanδm)，εr,μr分别是相对介电常数和相对磁导率，tanδ,tanδm分别是电场和磁场中的损耗角正切值，下标m代表磁场。
在热-电一体化分析中，温度场对于天线罩的影响主要体现在介质厚度、介电常数和损耗角正切上。考虑介电常数和损耗角正切在气动热载荷下的变化，可以得到：
其中，Tt为温度场作用下n级四端网络的网络总联矩阵，下标t代表此矩阵受温度场的影响。下标n代表等效网络的级数。在温度场作用下的第n层介质平板的第i层介质的传输矩阵可以表示为：
其中，下标t代表参数受温度场的影响，
Zti温度场下为第i层变形的介质相对于自由空间的归一化特征阻抗，可以表示为：
θti为电磁波从温度场下第i-1层变形介质透射到第i层变形介质的入射角，并且满足斯涅尔定理
εti、μti、εt(i-1)、μt(i-1)表示电磁波从温度场下第i-1层变形介质透射到第i层变形介质时，第i层的介电常数、磁导率，第i-1层的介电常数、磁导率，下标t代表参数受温度场的影响。λ0为电磁波在真空中的波长，“sh”和“ch”分别为双曲正弦和双曲余弦函数，Z0为电磁波在空气中的阻抗，dti表示在温度场下发生变形的第i层介质厚度；
最终可以得到温度场作用下的天线罩介质平板的传输系数Tt：
At、Bt、Ct、Dt可由网络总联矩阵Tt得到；
最终可以得到天线罩在温度场作用下的功率传输系数：
Ttv＝|Tt|2
下标tv代表功率传输系数受到温度场的影响。
更进一步地，所述步骤S7完成对于S1.2中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成原始状态天线罩透波性能分析，具体包括：
S7.1：将S1.2生成的inp文件导入FEKO中，设置电磁材料参数；
S7.2：建立水平极化的振子天线模型，并设置激励方式和大小；
S7.3：设置远场参数、求解截断误差；
S7.4：将S7.3生成的模型文件提交进行计算，获得原始状态下天线罩的透波性能，具体计算方式如下：
由不同层数的介质板依次叠组成的天线罩模型相当于是不同的传输矩阵的级联，如图2所示，θi为天线罩内天线产生的电磁波在罩内入射时的入射角，di是第i层介质的厚度，εi,μi分别是第i层介质的材料的介电常数、磁导率，数学表达式为：εi＝εr(1-jtanδ),μi＝μr(1-jtanδm)，εr,μr分别是相对介电常数和相对磁导率,tanδ,tanδm分别是电场和磁场中的损耗角正切值，下标i代表介质层数。
对于n层介质板，可以等效为n级四端网络，其网络总联矩阵为S：
其中，第n层介质平板的第i层介质的传输矩阵可以表示为：
其中，
λ0为电磁波在真空中的波长，sh和ch分别为双曲正弦和双曲余弦函数，Zi为第i层介质相对于自由空间的归一化特征阻抗，垂直极化的归一化特征阻抗为
Z0为电磁波在空气中的阻抗，下标0代表电磁波传播介质为真空，θi为电磁波从i-1层介质透射到第i层介质的入射角，并且满足斯涅尔定理
di表示第i层介质的厚度；εi、μi、εi-1、μi-1表示电磁波第i-1层介质透射到第i层介质时，第i层的介电常数、磁导率，第i-1层的介电常数、磁导率；
综上，对于介质平板的传输系数Q可以表示为：
A、B、C、D可由网络总联矩阵S得到；
因此，天线罩的功率传输系数为：Tv＝|Q|2，其中，v代表功率传输系数。
有益效果：与现有技术相比，本发明的技术方案具有以下有益技术效果：
本发明提出一种考虑热效应的天线罩电磁性能分析方法，该方法从电磁学等效四端口网络出发，将天线罩实际工作状态下的温度场作为电磁学分析初始条件，能够在天线罩初期设计阶段统筹考虑其力学领域的温度场和电磁学领域的透波性能。此分析方法能够突破在天线罩初期设计阶段中力学和电磁学的学科界限，极大程度上缩短在天线罩设计初期力学设计团队和电磁学设计团队的迭代周期和相关设计成本。具有较高的工程价值。
附图说明
图1热电耦合分析流程图；
图2多层平板电磁波传输示意图；
图3气动热载荷下天线罩透波性能曲线。
具体实施方式
下面对本发明技术方案进行详细说明，但是本发明的保护范围不局限于所述实施例。
本发明提出一种考虑热效应的天线罩电磁性能分析方法，具体包括以下步骤：
S1：建立天线罩力学计算有限元模型、电磁学计算未变形有限元模型；
S1.1：基于HyperMesh软件为天线罩的几何结构划分六面体单元有限元网格，初步建立力学计算有限元模型。
具体的六面体单元尺寸为：100mm。
S1.2：基于HyperMesh软件将天线罩几何结构的外表面划分为三角形单元有限元网格，建立电磁学计算未变形有限元模型。
具体的三角形单元尺寸为：10mm。
S2：对S1建立的力学有限元模型添加材料属性、施加温度载荷、边界条件；
S2.1：添加材料属性，并将天线罩子头部施加1100℃的热流，尾部时间200℃的热流。
具体的，材料密度：1.060×103(Kg/m3)；材料弹性模量：61.3GPA；泊松比：0.325；热导率：2.35W/m·K；比热容：765J/Kg·K。
S2.2：将上述所有的模型输出为inp文件，基于Abaqus求解器，对S2.1的有限元模型进行传热分析，获得其所有节点的稳态节点温度。
S2.3：将S2.2中的稳态节点温度作为温度载荷施加于S1.1建立的有限元模型。
S2.4：将S2.3见建立的有限元模型的尾部施加固支的边界条件。
S3：完成对静力学问题的求解，对计算结果进行后处理，提取节点位移；
S3.1：求解时间：10s、输出变量和参数：节点位移，单元应力。
S3.2：基于S3.2的静力学求解算法，将S3.1建立的模型输出为inp文件进行静力学求解。
S4：基于Python语言完重构S3有限元模型，得到温度载荷下变形的有限元模型；
S4.1：基于Python脚本读取S3.2生成的inp文件，获得节点的X、Y、Z坐标值。
S4.2：将S4.1读取的节点坐标值与S3.3获取的X、Y、Z方向上的位移：ux uy uz位移值进行算数求和。
S4.3：将S4.2求和之后的节点坐标重新写入S3.2生成的inp文件。
S5：基于HyperMesh软件对S4中变形的有限元模型进行重构得到电磁学计算有限元模型；
S5.1：将S4.3生成的inp文件重新导入HyperMesh软件，通过软件提取天线罩模型的外表面有限元模型。
S5.2：将S5.1提取的外表面有限元模型的类型改为三角形，大小改为波源的五分之一，进行重新绘制，输出新的inp文件。
S6：完成对于S6中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成温度载荷下天线罩透波性能分析；
S6.1：将S5.2生成的inp文件导入FEKO中，设置电磁材料参数。
S6.2：建立垂直极化的振子天线模型，并设置激励方式：电压激励；激励大小为：1V。
S6.3：设置多层快速多极子求解算法、远场参数：天线全局坐标系下的XOZ平面、YOZ平面、求解截断：0.018。
S6.4：基于S6.3的电磁场求解算法，将S6.3生成的模型文件提交进行计算，获得温度场下天线罩的透波性能。
S7：完成对于S1.2中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成原始状态天线罩透波性能分析；
S7.1：将S1.2生成的inp文件导入FEKO中，设置电磁材料参数。
S7.2：建立水平极化的振子天线模型，并设置激励方式：电压激励；激励大小为：1V。
S7.3：设置多层快速多极子求解算法、远场参数：天线全局坐标系下的XOZ平面、YOZ平面、求解截断：0.018。
S7.4：基于S7.3的电磁场求解算法，将S7.3生成的模型文件提交进行计算，获得原始状态下天线罩的透波性能。计算结果见图3所示。
以上示意性的对本发明及其实施方式进行了描述，该描述没有限制性，附图中所示的也只是本发明的实施方式之一，实际的结构和方法并不局限于此。所以，如果本领域的普通技术人员受其启示，在不脱离本发明创造宗旨的情况下，不经创造性的设计出与该技术方案相似的结构方式及实施例，均属于本发明的保护范围。
Claims (7)
1.一种考虑热效应的天线罩电磁性能分析方法，其特征在于，所述方法具体包括如下步骤：
S1：建立天线罩力学计算有限元模型、电磁学计算未变形有限元模型；
S2：对S1建立的力学计算有限元模型添加材料属性、施加温度载荷、边界条件；
S3：完成对静力学问题的求解，对计算结果进行后处理，提取节点位移，具体包括：
S3.1：对S1建立的力学计算有限元模型设置传热分析求解类型、求解时间、输出变量和输出参数；
S3.2：将S3.1完成参数设置的力学计算有限元模型输出为inp文件进行静力学求解，求解方式如下：
(3.1.1)首先导出热弹性力学的基本方程和边界条件，命弹性体的各点的变温ΔT＝T2-T1，即是后一瞬时的温度减去前一瞬时的温度，ΔT升温时为正，降温时为负，温度发生变化时，弹性体的热应变为：αΔT，其中，α是天线罩弹性单元弹性体的线热膨胀系数，ΔT为温度变化量，在α不随温度的变化发生变化时，可以得到天线罩弹性单元内各个点的应变分量：εx＝εy＝εz＝αΔT，γxy＝γyz＝γzx＝0，εx，εy，εz分别是变形体在三维空间中X，Y，Z的正应变，γxy＝γyz＝γzx分别是变形体在三维空间中XY，ZY，ZX平面的剪切应变；
(3.1.2)天线罩弹性单元受到外在约束以及内部的各个点约束时并不能发生线性应变，就会产生温度应力，因此可以得到式(1)：
σx，σy，σz分别是变形体在三维空间中X，Y，Z的正应力，τxy，τyz，τzx分别是变形体在三维空间中XY，ZY，ZX平面的剪切应力，E为弹性模量，μ为泊松比；
在只有变温ΔT存在的情况下，温度场的位移和应力响应可以等效为平面问题：
σx＝0，τyz＝0，τzx＝0
因此，式(1)可以表示为：
由上式可以求得温度场下的天线罩弹性单元应变分量和由变温ΔT表示的应力分量的物理方程：
可以得到相似的天线罩弹性单元平面应变问题的热弹性力学的物理方程：
弹性力学平面问题的几何方程和平衡微方程为：
u和v分别是平面内X方向和Y方向位移；
按照平面应变问题进行公式等效可以得到：
将上式导入平衡微方程，化简后可以得到：
由无面力应力条件边界条件可知求解应力平面问题的应力边界条件：
l，m为边界外法线方向和X、Y方向的夹角的余弦，位移边界条件仍然为：
us，
vs，
为已知边界平面上X方向和Y方向位移；
将上述位移坐标转化为全局坐标系，即可得到结构在X、Y、Z方向上的位移：ux uy uz；
S4：基于Python语言重构S3.1中完成参数设置的力学计算有限元模型，得到温度载荷下变形的有限元模型；
S5：基于HyperMesh软件对S4中变形的有限元模型进行重构得到电磁学计算有限元模型；
S6：完成对于S5中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成温度载荷下天线罩透波性能分析；
S7：完成对于S|中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成原始状态天线罩透波性能分析，并和S6的分析结果进行对比，获得气动热载荷对于天线罩透波性能的影响。
2.根据权利要求1所述的一种考虑热效应的天线罩电磁性能分析方法，其特征在于，所述步骤S1建立天线罩力学计算有限元模型、电磁学计算未变形有限元模型，具体包括：
S1.1：基于HyperMesh软件为天线罩的几何结构划分六面体单元有限元网格，初步建立力学计算有限元模型；
S1.2：基于HyperMesh软件将天线罩几何结构的外表面划分为三角形单元有限元网格，建立电磁学计算未变形有限元模型；
S1.3：对步骤S1.1建立的有限元模型设置材料属性。
3.根据权利要求2所述的一种考虑热效应的天线罩电磁性能分析方法，其特征在于，所述步骤S2对S1建立的力学有限元模型添加材料属性、施加温度载荷、边界条件，具体包括：
S2.1：添加材料属性，并将天线罩风洞试验实验数据中天线罩头部和尾部的实际温度施加于S1.1建立的有限元模型上；
S2.2：将所有的模型输出为inp文件，基于Abaqus求解器，对S2.1的有限元模型进行传热分析，获得其所有节点的稳态节点温度；
S2.3：将S2.2中的稳态节点温度作为温度载荷施加于S1.1建立的有限元模型；
S2.4：将S2.3中建立的有限元模型的尾部施加固支的边界条件。
4.根据权利要求3所述的一种考虑热效应的天线罩电磁性能分析方法，其特征在于，所述步骤S4基于Python语言重构S3有限元模型，得到温度载荷下变形的有限元模型，具体包括：
S4.1：基于Python脚本读取S3.2生成的inp文件，获得节点的X、Y、Z坐标值；
S4.2：将S4.1读取的节点坐标值与S3.2获取的X、Y、Z方向上的位移ux uy uz位移值进行算数求和；
S4.3：将S4.2求和之后的节点坐标重新写入S3.2生成的inp文件。
5.根据权利要求4所述的一种考虑热效应的天线罩电磁性能分析方法，其特征在于，所述步骤S5基于HyperMesh软件对S4中变形的有限元模型进行重构得到电磁学计算有限元模型，具体包括：
S5.1：将S4.3生成的inp文件重新导入HyperMesh软件，通过软件提取天线罩模型的外表面有限元模型；
S5.2：将S5.1提取的外表面有限元模型的类型改为三角形，有限元模型三角形单元周长改为波源波长的五分之一，进行重新绘制，输出新的inp文件。
6.根据权利要求5所述的一种考虑热效应的天线罩电磁性能分析方法，其特征在于，所述步骤S6完成对于S5中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成温度载荷下天线罩透波性能分析，具体包括：
S6.1：将S5.2生成的inp文件导入FEKO中，设置电磁材料参数；
S6.2：在S6.1生成的天线罩电磁学计算有限元模型的基础上，建立垂直极化的振子天线模型，并设置激励方式和大小；
S6.3：对S6.2建立的天线和天线罩电磁学计算有限元模型设置远场参数、求解截断误差；
S6.4：为获得温度场下天线罩的透波性能，对S6.3中远场参数、求解截断误差设置完成的天线和天线罩模型按照下式的计算方式进行计算：
(6.4.1)设θi为天线罩内天线产生的电磁波在罩内入射时的入射角，di是第i层介质的厚度，ε0，εi，μ0，μi分别是在空气或真空中介电常数、介质材料的相对介电常数、空气或真空中磁导率、介质材料的相对磁导率，数学表达式为：εi＝εr(1-j tanδ)，μi＝μr(1-j tanδm)，εr，μr分别是相对介电常数和相对磁导率，tanδ，tanδm分别是电场和磁场中的损耗角正切值，下标m代表磁场；
(6.4.2)考虑介电常数和损耗角正切在气动热载荷下的变化，可以得到：
其中，Tt为温度场作用下n级四端网络的网络总联矩阵，下标t代表此矩阵受温度场的影响，下标n代表等效网络的级数，在温度场作用下的第n层介质平板的第i层介质的传输矩阵可以表示为：
其中，下标t代表参数受温度场的影响，
Zti温度场下为第i层变形的介质相对于自由空间的归一化特征阻抗，可以表示为：
θti为电磁波从温度场下第i-1层变形介质透射到第i层变形介质的入射角，并且满足斯涅尔定理
εti、μti、εt(i-1)、μt(i-1)表示电磁波从温度场下第i-1层变形介质透射到第i层变形介质时，第i层的介电常数、磁导率，第i-1层的介电常数、磁导率，下标t代表参数受温度场的影响，λ0为电磁波在真空中的波长，sh和ch分别为双曲正弦和双曲余弦函数，Z0为电磁波在空气中的阻抗，dti表示在温度场下发生变形的第i层介质厚度；
最终可以得到温度场作用下的天线罩介质平板的传输系数Tt：
At、Bt、Ct、Dt可由网络总联矩阵Tt得到；
最终可以得到天线罩在温度场作用下的功率传输系数：
Ttv＝|Tt|2
下标tv代表功率传输系数受到温度场的影响。
7.根据权利要求6所述的一种考虑热效应的天线罩电磁性能分析方法，其特征在于，所述步骤S7完成对于S1中建立的电磁学计算有限元模型的电磁介质参数设置以及天线建模，完成原始状态天线罩透波性能分析，具体包括：
S7.1：将S1.2生成的inp文件导入FEKO中，设置电磁材料参数；
S7.2：建立水平极化的振子天线模型，并设置激励方式和大小；
S7.3：对S7.2建立的天线模型设置远场参数、求解截断误差；
S7.4：将S7.3中的天线模型提交进行计算，获得原始状态下天线罩的透波性能，具体计算方式如下：
(7.3.1)设置θi为天线罩内天线产生的电磁波在罩内入射时的入射角，di是第i层介质的厚度，εi，μi分别是第i层介质的材料的介电常数、磁导率，表达式为：εi＝εr(1-j tanδ)，μi＝μr(1-j tanδm)，εr，μr分别是相对介电常数和相对磁导率，tanδ，tanδm分别是电场和磁场中的损耗角正切值，下标i代表介质层数；
(7.3.2)对于n层介质板，可以等效为n级四端网络，其网络总联矩阵为S：
其中，第n层介质平板的第i层介质的传输矩阵可以表示为：
其中，
λ0为电磁波在真空中的波长，sh和ch分别为双曲正弦和双曲余弦函数，Zi为第i层介质相对于自由空间的归一化特征阻抗，垂直极化的归一化特征阻抗为
Z0为电磁波在空气中的阻抗，下标0代表电磁波传播介质为真空，θi为电磁波从i-1层介质透射到第i层介质的入射角，并且满足斯涅尔定理
di表示第i层介质的厚度；εi、μi、εi-1、μi-1表示电磁波第i-1层介质透射到第i层介质时，第i层的介电常数、磁导率，第i-1层的介电常数、磁导率；
综上，对于介质平板的传输系数Q可以表示为：
A、B、C、D可由网络总联矩阵S得到；
因此，天线罩的功率传输系数为：Tv＝|Q|2
其中，v代表功率传输系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010047861.0A CN111274726B (zh) | 2020-01-16 | 2020-01-16 | 一种考虑热效应的天线罩电磁性能分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010047861.0A CN111274726B (zh) | 2020-01-16 | 2020-01-16 | 一种考虑热效应的天线罩电磁性能分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111274726A CN111274726A (zh) | 2020-06-12 |
| CN111274726B true CN111274726B (zh) | 2021-03-16 |
Family
ID=71001088
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010047861.0A Active CN111274726B (zh) | 2020-01-16 | 2020-01-16 | 一种考虑热效应的天线罩电磁性能分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111274726B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112199880A (zh) * | 2020-10-27 | 2021-01-08 | 北京机电工程总体设计部 | 一种天线罩参数化及强度快速分析仿真系统 |
| CN114398716A (zh) * | 2021-12-14 | 2022-04-26 | 北京遥感设备研究所 | 一种热力电三场耦合下天线罩电性能全波数值计算方法 |
| CN114492102B (zh) * | 2021-12-22 | 2024-06-11 | 中国航天空气动力技术研究院 | 编织防热材料考虑物性演变的传热及力学特性确定方法 |
| CN116187134B (zh) * | 2023-02-10 | 2023-11-21 | 西安电子科技大学 | 网状天线金属丝网透射反射系数分析方法 |
| CN116090307B (zh) * | 2023-02-15 | 2024-08-09 | 哈尔滨锅炉厂有限责任公司 | 一种基于正交异性本构的pche芯体的有限元建模方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105404751A (zh) * | 2015-12-07 | 2016-03-16 | 中国电子科技集团公司第三十二研究所 | 基于热-力-电磁场网络统一的实现方法 |
| CN106528954A (zh) * | 2016-10-19 | 2017-03-22 | 西安电子科技大学 | 空间网状天线温度载荷作用下机电集成优化设计方法 |
| CN110470916A (zh) * | 2019-08-16 | 2019-11-19 | 西安电子科技大学 | 日照温度场作用下反射面天线电性能快速评估与补偿方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108920831B (zh) * | 2018-06-29 | 2020-11-03 | 西安电子科技大学 | 高速飞行器罩体高温烧蚀对天线电性能影响的快速计算方法 |
-
2020
- 2020-01-16 CN CN202010047861.0A patent/CN111274726B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105404751A (zh) * | 2015-12-07 | 2016-03-16 | 中国电子科技集团公司第三十二研究所 | 基于热-力-电磁场网络统一的实现方法 |
| CN106528954A (zh) * | 2016-10-19 | 2017-03-22 | 西安电子科技大学 | 空间网状天线温度载荷作用下机电集成优化设计方法 |
| CN110470916A (zh) * | 2019-08-16 | 2019-11-19 | 西安电子科技大学 | 日照温度场作用下反射面天线电性能快速评估与补偿方法 |
Non-Patent Citations (3)
| Title |
|---|
| Dynamic design of stiffeners for a typical panel using topology optimization and finite element analysis;Yadong Zhou and 等;《Advances in Mechanical Engineering》;20151231;第1-8页 * |
| 基于有限元_边界元的热_声_振耦合分析;李彦斌;《中国力学大会》;20131231;第121页 * |
| 空间热效应对网状天线反射面面型精度的影响分析;周华楠;《中国硕士学位论文全文数据库》;20161231;第1-73页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111274726A (zh) | 2020-06-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111274726B (zh) | 一种考虑热效应的天线罩电磁性能分析方法 | |
| CN111241736B (zh) | 一种天线罩力热电一体化分析方法 | |
| Karakoti et al. | Nonlinear transient analysis of porous P-FGM and S-FGM sandwich plates and shell panels under blast loading and thermal environment | |
| CN110457790B (zh) | 用于结构变形分析的近场动力学非连续伽辽金有限元方法 | |
| Wang et al. | An interpolating boundary element-free method with nonsingular weight function for two-dimensional potential problems | |
| Mahapatra et al. | Nonlinear flexural analysis of laminated composite panel under hygro-thermo-mechanical loading—a micromechanical approach | |
| Singh et al. | Buckling of laminated composite plates subjected to mechanical and thermal loads using meshless collocations | |
| Shegokar et al. | Stochastic finite element nonlinear free vibration analysis of piezoelectric functionally graded materials beam subjected to thermo-piezoelectric loadings with material uncertainties | |
| Peng et al. | A thick shell model based on reproducing kernel particle method and its application in geometrically nonlinear analysis | |
| Cao et al. | Hybrid fundamental solution based finite element method: theory and applications | |
| Ghannadpour et al. | Application of plate decomposition technique in nonlinear and post-buckling analysis of functionally graded plates containing crack | |
| Ye et al. | Buckling analysis of three-dimensional functionally graded sandwich plates using two-dimensional scaled boundary finite element method | |
| Rahmani et al. | Vibration analysis for anti-symmetric laminated composite plates resting on visco-elastic foundation with temperature effects | |
| Jin et al. | Effect of fiber arrangement on residual thermal stress distributions in a unidirectional composite | |
| Nguyen-Ngoc et al. | Three-dimensional polyhedral finite element method for the analysis of multi-directional functionally graded solid shells | |
| Gupta et al. | Static and dynamic NURBS-based isogeometric analysis of composite plates under hygrothermal environment | |
| Zhou et al. | Evaluation of performance of magneto-electro-elastic sensor subjected to thermal-moisture coupled load via CS-FEM | |
| Li et al. | Coupled thermo-mechanical finite element models with node-dependent kinematics for multi-layered shell structures | |
| Dewangan et al. | Numerical thermoelastic eigenfrequency prediction of damaged layered shell panel with concentric/eccentric cutout and corrugated (TD/TID) properties | |
| CN111274725B (zh) | 一种力载荷下天线罩电磁性能分析方法 | |
| Mishra | Fracture analysis of interfacial crack in piezoelectric bimaterial by XIGA approach using Bézier extraction of NURBS | |
| Girish et al. | Thermomechanical postbuckling analysis of symmetric and antisymmetric composite plates with imperfections | |
| Wang et al. | A strip element method for the transient analysis of symmetric laminated plates | |
| Fukunaga et al. | FEM modeling of adaptive composite structures using a reduced higher-order plate theory via penalty functions | |
| Akhras et al. | Three-dimensional thermal buckling analysis of piezoelectric composite plates using the finite layer method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |