当前位置：首页 > news >正文

探索十字型声子晶体案例：从原理到代码实现

news 2026/5/12 5:23:31

十字型声子晶体案例

在材料科学与声学领域，声子晶体是一种极具潜力的人造周期性复合材料，它能够调控弹性波或声波的传播，就像半导体对电子的调控一样。今天咱们就来深入探究一下十字型声子晶体这个有趣的案例。

十字型声子晶体的原理基础

声子晶体依靠其周期性结构产生带隙，当弹性波或声波的频率落在这个带隙范围内时，波的传播就会被抑制。十字型声子晶体，因其独特的十字状单元结构而得名。这种特殊的结构赋予了它一些不同于常规声子晶体的声学特性。

想象一下，十字型结构在空间中以一定周期排列，不同介质（比如固体和空气）组成的十字结构会对弹性波产生散射和干涉。这些散射波相互作用，在某些频率范围内就会导致波的抵消，从而形成带隙。

代码实现与分析

下面咱们用Python结合有限元方法的相关库来模拟十字型声子晶体的声学特性。这里我们借助meep库，它是一个功能强大的用于电磁学和声学模拟的开源软件包。

import meep as mp import numpy as np import matplotlib.pyplot as plt # 定义晶格常数 a = 1.0 # 定义单元结构尺寸 rod_width = 0.2 * a rod_length = 0.8 * a # 定义晶格 sx = 1.0 * a sy = 1.0 * a cell_size = mp.Vector3(sx, sy) # 定义材料 mat_rod = mp.Medium(epsilon=12) # 定义几何结构 geometry = [ mp.Block( material=mat_rod, size=mp.Vector3(rod_width, rod_length), center=mp.Vector3(0, 0.5 * (sy - rod_length)) ), mp.Block( material=mat_rod, size=mp.Vector3(rod_length, rod_width), center=mp.Vector3(0.5 * (sx - rod_length), 0) ) ] # 定义边界条件 pml_layers = [mp.PML(0.3)] # 定义频率范围 fmin = 0.2 fmax = 0.5 nfreq = 50 # 定义模拟对象 sim = mp.Simulation( cell_size=cell_size, geometry=geometry, boundary_layers=pml_layers, resolution=50 ) # 定义频域源 src = mp.Source( src=mp.GaussianSource((fmin + fmax) / 2, fwidth=fmax - fmin), component=mp.Ex, center=mp.Vector3() ) sim.append_near2far( (fmin + fmax) / 2, 0, 1, mp.Near2FarRegion(mp.Vector3(-0.5 * sx, 0), size=mp.Vector3(0, sy)), mp.Near2FarRegion(mp.Vector3(0.5 * sx, 0), size=mp.Vector3(0, sy)), mp.Near2FarRegion(mp.Vector3(0, -0.5 * sy), size=mp.Vector3(sx, 0)), mp.Near2FarRegion(mp.Vector3(0, 0.5 * sy), size=mp.Vector3(sx, 0)) ) sim.run( mp.at_beginning(mp.output_epsilon), mp.to_appended("efield", mp.at_every(0.1, mp.output_efield)), until_after_sources=mp.stop_when_fields_decayed(50, mp.Ex, mp.Vector3(), 1e-3) ) # 获取传输系数 flux_freqs = np.linspace(fmin, fmax, nfreq) flux = sim.get_fluxes() # 绘制传输谱 plt.plot(flux_freqs, flux) plt.xlabel('Frequency (1/a)') plt.ylabel('Transmission') plt.title('Transmission Spectrum of Cross - shaped Phononic Crystal') plt.grid(True) plt.show()

代码分析：

参数定义部分：首先我们定义了晶格常数a，它是整个结构周期性的基础。接着定义了十字型单元结构中棒的宽度rodwidth和长度rodlength，以及晶格的尺寸sx和sy。这里mat_rod定义了构成十字结构的材料特性，这里假设其相对介电常数为12，在声学模拟中可类比为材料的某种弹性参数。
几何结构构建：通过mp.Block来构建十字型结构。这里有两个mp.Block，一个水平放置，一个垂直放置，它们的中心位置和尺寸都根据之前定义的参数来确定。
模拟设置：pml_layers定义了完美匹配层（PML），用于吸收边界处的波，避免反射干扰模拟结果。sim这个对象整合了晶格尺寸、几何结构、边界条件和分辨率等信息。
源与监测：定义了一个高斯脉冲源src，它位于结构中心，用于激发弹性波。sim.append_near2far部分设置了近场到远场的转换区域，以便后续获取传输系数。
模拟运行与结果处理：sim.run启动模拟，这里设置了在开始时输出介电常数分布，并且每隔0.1个时间步输出电场分布，直到场强衰减到一定程度停止。最后获取不同频率下的传输系数，并绘制传输谱。通过这个传输谱，我们就能直观地看到在哪些频率范围内，十字型声子晶体对波的传输起到抑制或增强作用，也就是分析带隙的分布情况。