openEMS电磁仿真在工业无线仪表中的应用：从天线设计到信号完整性

📶 4G Cat1 无线仪表

工作频段：B1/B3/B5/B8/B34/B38/B39/B40/B41（约 700MHz ~ 2.7GHz），典型带宽 10MHz/20MHz。

在这个频段，天线尺寸和PCB走线尺寸都接近信号波长（以B3为例，λ≈16cm），天线匹配网络和射频走线的阻抗控制直接决定通信质量。

# 4G Cat1 模块射频前端仿真示例（Python + openEMS）
# 仿真目标：1.6GHz微带线 + π型匹配网络

import openEMS
import CSXCAD
import numpy as np

# 仿真频段（LTE B3: 1710-1785MHz 上行）
f_low = 1.71e9
f_high = 1.785e9
f_center = (f_low + f_high) / 2  # 1.7475 GHz

# 创建仿真
sim = openEMS.OpenEMS(f_start=f_low, f_stop=f_high)
sim.kernels = {'n_threads': 8}  # 8线程并行

# PCB叠层参数（4层板，FR4）
stackup = {
    'L1_top': {'material': 'copper', 'thick': 0.035},  # 顶层铜
    'L1_diel': {'epsR': 4.4, 'thick': 0.2},            # FR4介质
    'L2_gnd': {'material': 'copper', 'thick': 0.035},  # 地平面
}

# 微带线设计（50Ω特征阻抗）
# W = 0.3mm, H = 0.2mm (FR4), 目标 Z0 ≈ 50Ω
width = 0.3e-3    # 线宽 0.3mm
height = 0.2e-3   # 介质厚度 0.2mm
length = 10e-3    # 走线长度 10mm

# 创建微带线几何
msline = CSXCAD.MicrostripLine(
    L=length, W=width, H=height,
    EpsR=4.4, MuR=1,
    Z0=50  # 目标特征阻抗
)

# 端口设置（波导端口激励）
port1 = openEMS.WaveguidePort('L1_top', 'microstrip',
                               L=length, W=width,
                               f_center=f_center, f_range=[f_low, f_high])

# 匹配网络（电感-电容-电感 π型）
# 目标：在 1.7475GHz 将 50Ω 匹配到天线阻抗 Z_ant
# L1 = 3.9nH, C = 2.7pF, L2 = 3.9nH
L_match = CSXCAD.Inductor(L=3.9e-9)   # 匹配电感
C_match = CSXCAD.Capacitor(C=2.7e-12)  # 匹配电容

# 运行仿真
sim.add_geometry([msline, L_match, C_match, port1])
results = sim.run()

# 提取S参数
S11 = results['S11']  # 回波损耗
S21 = results['S21']  # 插入损耗

print(f"在 {f_center/1e9:.2f} GHz:")
print(f"  S11 = {S11:.2f} dB (目标 < -10dB)")
print(f"  S21 = {S21:.2f} dB (目标 > -1dB)")

关键仿真指标：

📡 LoRa 无线仪表

工作频段：470MHz ~ 510MHz（中国LoRa频段），调制方式为扩频（Chirp Spread Spectrum），带宽 125/250/500kHz。

相比4G Cat1，LoRa频率更低（λ≈60cm），天线尺寸更大，但PCB走线的阻抗控制和低噪声放大器输入匹配仍然是设计难点。

# LoRa 470MHz 缝隙天线仿真（openEMS）
# 仿真目标：470MHz 缝隙天线，PCB尺寸 30mm x 15mm

import openEMS
import CSXCAD
import numpy as np

f_center = 470e6   # 470 MHz
f_range = [460e6, 480e6]

sim = openEMS.OpenEMS(f_start=f_range[0], f_stop=f_range[1])

# PCB substrate (Rogers RO4003, 低损耗板材)
substrate = CSXCAD.Box(
    name='substrate',
    size=[30e-3, 15e-3, 0.8e-3],  # 30mm x 15mm x 0.8mm
    material={'epsR': 3.38, 'kappa': 0.0001}  # Rogers RO4003C
)

# 地平面开缝（缝隙天线）
slot = CSXCAD.Box(
    name='slot',
    size=[20e-3, 1.5e-3],  # 缝隙 20mm x 1.5mm
    material='air',
    layer='L1_top'
)

# 馈电点（偏移缝隙中心）
feed = CSXCAD.Feed(
    name='feed',
    x=10e-3, y=0, z=0.8e-3,  # 缝隙中点
    f_center=f_center
)

# 近场-远场变换（计算天线方向图）
nf2ff = openEMS.NF2FF(
    center=[15e-3, 7.5e-3, 5e-3],
    size=[50e-3, 50e-3, 50e-3],  # NF2FF 球半径
    f_center=f_center
)

# 运行仿真
sim.add_geometry([substrate, slot, feed, nf2ff])
results = sim.run()

# 提取天线参数
gain_dbi = nf2ff.calc_gain()  # 增益 (dBi)
efficiency = nf2ff.calc_efficiency()  # 辐射效率
pattern = nf2ff.calc_pattern_2d()  # 2D 方向图

print(f"470MHz 缝隙天线仿真结果:")
print(f"  增益: {gain_dbi:.1f} dBi (目标 > -2 dBi)")
print(f"  效率: {efficiency*100:.1f}% (目标 > 50%)")

LoRa仪表电磁仿真重点：

• PCB尺寸与波长的关系（30mm ≈ λ/4）
• 馈电网络与匹配（50Ω → 天线阻抗）
• 地平面完整性对方向图的影响
• 金属外壳对天线的遮挡效应

📡 雷达液位计

工作频段：26GHz / 80GHz（微波雷达 FMCW），测量原理是发射微波信号到被测介质表面，接收反射信号，通过飞行时间（ToF）计算距离。

这是微波频段最高频的应用，λ≈11mm（26GHz）或 3.75mm（80GHz），天线设计的微小误差都会显著影响波束和测量精度。

# 雷达液位计 26GHz 喇叭天线仿真（openEMS）
# 仿真目标：26GHz 波导喇叭天线，增益 > 20dBi

import openEMS
import CSXCAD
import numpy as np

f_center = 26e9      # 26 GHz (K波段)
f_range = [24e9, 28e9]

sim = openEMS.OpenEMS(f_start=f_range[0], f_stop=f_range[1])

# 标准矩形波导 WR-34 (尺寸: 86.36mm x 43.18mm)
# 但液位计天线是小型化喇叭，这里用锥形喇叭
flange_length = 30e-3   # 喇叭长度 30mm
flange_width = 20e-3    # E面宽度 20mm
flange_height = 15e-3   # H面高度 15mm

# 锥形喇叭几何
horn = CSXCAD.PyramidalHorn(
    L=flange_length,
    W=flange_width,
    H=flange_height,
    throat_length=5e-3
)

# 介质透镜（可选，提升方向性）
lens = CSXCAD.Box(
    name='lens',
    size=[flange_width, flange_height, 3e-3],
    material={'epsR': 2.1, 'kappa': 0},  # PTFE 透镜 epsR=2.1
    position=[0, 0, flange_length]
)

# 波导馈电
waveguide = CSXCAD.RectangularWaveguide(
    a=3.1e-3, b=1.55e-3,  # WR-34 波导内径
    L=10e-3
)

# 近场-远场变换
nf2ff = openEMS.NF2FF(
    center=[0, 0, flange_length + 10e-3],
    distance='auto',
    f_center=f_center
)

sim.add_geometry([horn, lens, waveguide, nf2ff])
results = sim.run()

# 分析方向图
theta_range = np.linspace(-90, 90, 361)  # 方位面
phi_range = np.linspace(-90, 90, 181)   # 俯仰面
pattern_3d = nf2ff.calc_pattern_3d(theta_range, phi_range)

# 计算 -3dB 波束宽度
hp_theta = nf2ff.calc_half_power_beamwidth('theta')  # E面
hp_phi = nf2ff.calc_half_power_beamwidth('phi')        # H面

print(f"26GHz 喇叭天线仿真结果:")
print(f"  增益: {nf2ff.calc_gain():.1f} dBi")
print(f"  E面波束宽度: {hp_theta:.1f}° (目标 < 20°)")
print(f"  H面波束宽度: {hp_phi:.1f}° (目标 < 25°)")
print(f"  副瓣电平: {nf2ff.calc_sidelobe_level():.1f} dB (目标 < -20dB)")

# 远场方向图（用于评估雷达分辨率）
radar_range_res = nf2ff.calc_range_resolution(
    c=3e8,         # 光速
    bandwidth=1e9   # FMCW 带宽 1GHz
)
print(f"  理论距离分辨率: {radar_range_res:.2f} m (带宽1GHz时约15cm)")

雷达液位计仿真核心指标：

⚡ 模拟前端PCB（信号完整性）

仪表的模拟前端（AFE）处理传感器的小信号（μV~mV级），对噪声极为敏感。PCB上微弱的射频干扰（RFI）会直接叠加在测量信号上，导致测量误差。

# 模拟电路PCB 噪声耦合仿真（openEMS）
# 仿真目标：检测相邻走线间的串扰（Crosstalk）强度

import openEMS
import CSXCAD
import numpy as np

f_range = [0, 1e9]  # 仿真DC~1GHz，看宽带噪声耦合

sim = openEMS.OpenEMS(f_start=f_range[0], f_stop=f_range[1])

# FR4 PCB叠层（4层板）
stackup = [
    {'name': 'L1_cu', 'thick': 0.035, 'material': 'copper'},
    {'name': 'diel1', 'thick': 0.2, 'material': {'epsR': 4.4, 'kappa': 0}},
    {'name': 'L2_gnd', 'thick': 0.035, 'material': 'copper'},
    {'name': 'diel2', 'thick': 1.2, 'material': {'epsR': 4.4, 'kappa': 0}},
    {'name': 'L3_pwr', 'thick': 0.035, 'material': 'copper'},
    {'name': 'diel3', 'thick': 0.2, 'material': {'epsR': 4.4, 'kappa': 0}},
    {'name': 'L4_cu', 'thick': 0.035, 'material': 'copper'},
]

# 受干扰走线（模拟小信号输入，靠近传感器）
victim = CSXCAD.MicrostripLine(
    name='victim',
    L=50e-3, W=0.2e-3, H=0.2e-3,
    EpsR=4.4, Z0=50,
    layer='L1_cu'
)

# 干扰源走线（数字开关信号，频率100MHz）
aggressor = CSXCAD.MicrostripLine(
    name='aggressor',
    L=50e-3, W=0.2e-3, H=0.2e-3,
    EpsR=4.4, Z0=50,
    layer='L1_cu',
    y_offset=1e-3  # 与 victim 间距 1mm
)

# 端口
port_victim = openEMS.WaveguidePort('victim', f_range=f_range)
port_aggressor = openEMS.WaveguidePort('aggressor', f_range=f_range)

# 端口阻抗分析（提取S参数）
results = sim.run()
S_vv = results['S11']  # 受害网络自反射
S_va = results['S21']  # 从干扰源到受害网络的耦合（串扰）

print(f"串扰仿真结果（间距1mm，长度50mm）:")
print(f"  NEXT (Near-End Crosstalk): {S_va[0]:.1f} dB @ 100MHz")
print(f"  FEXT (Far-End Crosstalk): {S_va[1]:.1f} dB @ 100MHz")

# 如果串扰过大，建议增加走线间距或加地孔屏蔽
improvement = sim.sweep({'spacing': [1e-3, 2e-3, 3e-3, 5e-3]})
print(f"  间距优化建议: ≥ {improvement['min_spacing']*1000:.1f} mm")

模拟电路PCB仿真重点：

• 近端串扰（NEXT）：受害网络接收端测得的耦合噪声
• 远端串扰（FEXT）：受害网络驱动端测得的耦合噪声
• 地 bounce：数字开关噪声通过地平面耦合到模拟地
• 电源完整性（PI）：PDN（电源分配网络）阻抗仿真