8.2-12.4GHz喇叭天线设计

记录个人的一些仿真案例。这里参考了李明洋的《HFSS天线设计》

设计目标

设计一个喇叭天线，其指标要求为：

工作频率范围	8.2-12.4GHz
电压驻波比	≤1.5
标称增益值	20dB
接口	BJ100

各频点增益要求为：

频率(GHz)	8.20	8.58	8.96	9.33	9.71	10.09
增益(dB)	19.48	19.79	20.09	20.36	20.61	20.86
频率(GHz)	10.47	10.85	11.22	11.60	11.98	12.40
增益(dB)	21.07	21.30	21.48	21.67	21.84	22.02

一、波导设计

由接口型号BJ100可知，喇叭天线底部为矩形波导，查阅尺寸表，得知其指标为：


a=22.86, b=10.16, A=25.4, B=12.7mm，仿真过程中，暂不考虑内圆角。
矩形波导长度一般取$5\lambda /4$，在这里，将频率定为区间的中心频率，在这里，将频率定为区间的中心频率。

二、喇叭口径计算

2.1 计算公式

由于矩形喇叭的辐射特型基本上是E面和H面扇形喇叭的结合，当波导中只是单模传播(只有TE波)时，矩形喇叭的口径可以通过远场公式推理出来。
此时，矩形喇叭的内径a1、b1和喇叭长度h_horn的计算公式为：
$$\begin{gathered} a_1^4-a{a}_1^3+\frac{3bG{\lambda }^2}{8\pi {\epsilon }_{ap}}{a}_1=\frac{3G^2{\lambda }^4}{32{\pi }^2{\epsilon }_{ap}^2} \\ h\_horn=\frac{a_1-a}{3\lambda }{a}_1 \\ {b}_1=\frac{1}{2}(b+\sqrt{b^2+8\lambda h_{horn}}) \end{gathered}$$
其中，G为无常纲量的增益值，!表示矩形喇叭的口径效率，进行最佳增益设计时，一般取值为0.5。

2.2 代码实现

基于上述公式，编写matlab代码计算各频率对应的喇叭口径参数如表value.xls所示：

clc, clear

gain_sheet = xlsread('D:\workbook\work_test\gain.xlsx');
[ro, co] = size(gain_sheet);

c = 3.0 * 10^8 * 10^3;   % 单位：mm / s
a = 22.86;  % 矩形波导的宽高
b = 10.16;

for i = 1:co
    freq = gain_sheet(1, i) * 10^9;
    gaindB = gain_sheet(2, i);  % dB增益
    gain = 10^(gaindB / 10);

    lambda = c / freq;

    syms a1
    a1 = double(solve(a1^4 - a * a1^3 + 3 * b * gain * lambda^2 * a1 / (8 * pi * 0.5) == 3 * gain^2 * lambda^4 / (32 * pi^2 * 0.5^2)));
    a1 = a1(2);
    h_horn = (a1 - a) * a1 / (3 * lambda);
    b1 = 0.5 * (b + sqrt(b^2 + 8 * lambda * h_horn));
    
    value(1, i) = a1;
    value(2, i) = h_horn;
    value(3, i) = b1;
end

value

xlswrite('value.xls', value);

在这里，首先以10.47GHz对应的参数为例，进行仿真验证。

三、CST仿真

3.1 波导创建

基于波导参数，建立矩形块，这里矩形块的长宽都为矩形波导的外径A、B，为了简化操作，先不进行挖空：

3.2 喇叭口径创建

1，首先建立一个厚度为0的面，面的长宽为喇叭的外口径值(a1+厚度, b1+厚度)，与波导的距离为h_horn：

2，选中两个面

3，使用Loft操作，将两个面合并起来。

4，使用Add操作，将得到的3个模块进行合并。然后，选中波导尾部和喇叭口部两个面，进行挖空操作。（Shape Tools->Shell）



建模完成，之后就是在矩形波导尾部添加波导端口激励，以及对应频率的远场监视器。并进行仿真

四、结果分析

4.1 VSWR

从图中可以看出，此结构的驻波比已达到要求。

4.2 增益

依次执行：Result Templates->Farfield Result->All Setting->将Plot Mode修改为Realized Gain，通过计算即可得到所有频率对应的最大增益，如表中所示。

频率	实际增益	目标增益	误差
8.2	19.76194	19.48	1.43%
8.58	20.33505	19.79	2.68%
8.96	20.40773	20.09	1.56%
9.33	20.67823	20.36	1.54%
9.71	20.96627	20.61	1.70%
10.09	21.06187	20.86	0.96%
10.47	21.36231	21.07	1.37%
10.85	21.45051	21.3	0.70%
11.22	21.67468	21.48	0.90%
11.6	21.77621	21.67	0.49%
11.98	21.78838	21.84	0.24%
12.4	22.19864	22.02	0.80%

从表中可以看出，该结构在8.58GHz的增益误差为2.68%>2.5%，需要对结构进行优化。

4.3 结构优化

使用后处理模板的优化算法，对a1、h_horn、b1进行优化，各参数的优化区间即为之前计算得到的参数范围。优化完成后，得出新的值为：a1=137.6mm；h_horn=204.11mm；b1=107.8mm。

频率	实际增益	目标增益	误差
8.2	19.68488	19.48	1.052%
8.58	19.84746	19.79	0.290%
8.96	20.25649	20.09	0.829%
9.33	20.44076	20.36	0.397%
9.71	20.6184	20.61	0.041%
10.09	20.97304	20.86	0.542%
10.47	21.04859	21.07	0.102%
10.85	21.3803	21.3	0.377%
11.22	21.43593	21.48	0.205%
11.6	21.7346	21.67	0.298%
11.98	21.80362	21.84	0.167%
12.4	21.85564	22.02	0.746%

优化完成后，各频率的增益值误差都符合要求。