6.1 waveform命令
功能:在模型中指定用于源的波形。
命令格式为:#waveform: str1 f1 f2 str2
str1:为波形的类型,它可以是
gaussian:是一个高斯波形。
gaussiandot:高斯波形的一阶导数。
gaussiandotnorm:高斯波形的一阶导数的归一化。
gaussiandotdot:高斯波形的二阶导数。
gaussiandotdotnorm:高斯波形的归一化二阶导数。
ricker:高斯波形的负的、归一化的二阶导数。
sine:正弦波形的一个单周期。
contsine :一个连续的正弦波形。为了避免在计算中引入噪声,对正弦波的第一个周期(斜坡
激励)调制波形的幅度。
f1:为波形最大幅度的缩放(对于#赫兹偶极子,单位为安培,对于#电压源或#传输线,单位为伏
特)。
f2:为波形的频率(赫兹)。
str2:用于将波形分配给源的标识符
例如,要指定一个振幅为1,中心频率为1.2 GHz的高斯波形。则指令为:
#waveform: gaussian 1 1.2e9 my_gauss_pulse.
注:用于创建波形的函数可以在工具部分找到
6.2 excitation_file命令
功能:指定一个ASCII文件,该文件包含幅度值列,这些幅度值列指定可以与模型中的源一起
使用的自定义波形形状。每列的第一行必须以标识符字符串开头,该字符串将用作每个波形的名
称。如果振幅值少于将要执行的迭代次数,则振幅值序列的末尾将用零值填充,直至迭代次数。如
果指定的幅度值超出需要,则忽略。
命令格式为:#excitation_file: file1
file1:为与输入文件在同一目录下包含指定波形的文件的名称,或者file1可以是包含指定波形的文
件的完整路径(允许您指定任何位置)。
例如,要指定包含两个自定义波形形状的文件my_waves.txt。则指令为:
#excitation_file: my_waves.txt
my_waves.txt文件的内容将采用如下形式:
然后使用my_pulse1自定义波形形状。
例如,一个z极化赫兹偶极子源。则指令为:
#hertzian_dipole: z 0.5 0.5 0.5 my_pulse1
注:使用自定义波形激励时,不需要使用#waveform命令,只需要使用#excitation_file命令和
任何将与自定义波形激励一起使用的源。
6.3 hertzian_dipole命令
功能:在电场位置指定电流密度项最简单的激励,通常称为加性或软源。
式中𝐽𝑠为电流密度,𝐼为电流,∆𝑙为无穷小电偶极子的长度,∆x,∆y,∆z为网格的空间分辨
率。根据指定的极化,∆𝑙被设置为∆x,∆y或∆z
命令格式为:#hertzian_dipole: c1 f1 f2 f3 str1 [f4 f5]
c1:为光源的偏振,可以是x、y或z。
f1、f2、f3:为模型中源的坐标(x、y、z)。
f4、f5:为可选参数。f4是启动源的时间延迟。f5是去除源的时间。
如果时间窗口长于源去除时间,则源将在源去除时间后停止。如果源移除时间长于时间窗
口,则源将rx间窗口内处于活动状态。如果省略f4 f5,则源将在时间窗口开始时启动,在时间窗口
结束时停止。
str1:为与源一起使用的波形的标识符。
例如,使用单位振幅的x极化赫兹偶极子和600 MHz中心频率的里克波形。则指令为:
#waveform: ricker 1 600e6 my_ricker_pulse
#hertzian_dipole: x 0.05 0.05 0.05 my_ricker_pulse
当#hertzian_dipole在2d模拟中使用时,它在模拟的无限方向上充当电流的线源。
6.4 magnetic_dipole命令
功能:模拟一个无穷小的磁偶极子。这通常被称为添加剂或软源。
命令格式为:#magnetic_dipole: c1 f1 f2 f3 str1 [f4 f5]
c1:为光源的偏振,可以是x、y或z。
f1、f2、f3:为模型中源的坐标(x、y、z)。
f4、f5:为可选参数。f4是启动源的时间延迟。f5是去除源的时间。
如果时间窗口长于源去除时间,则源将在源去除时间后停止。如果源移除时间长于时间窗
口,则源将rx间窗口内处于活动状态。如果省略f4 f5,则源将在时间窗口开始时启动,在时间窗口
结束时停止。
str1:为与源一起使用的波形的标识符。
6.5 voltage_source命令
功能:在电场位置引入电压源。如果它的电阻为零,即指定电场分量的时间变化是规定的,
它可以是硬源,或者如果它的电阻不为零,它可以作为电阻电压源。当模型中包含天线的物理特性
时,对天线的激励是有用的。
命令格式为:#voltage_source: c1 f1 f2 f3 f4 str1 [f5 f6]
c1:为光源的偏振,可以是x、y或z。
f1、f2、f3:为模型中源的坐标(x、y、z)。
f4、f5:为可选参数。f4是启动源的时间延迟。f5是去除源的时间。
如果时间窗口长于源去除时间,则源将在源去除时间后停止。如果源移除时间长于时间窗
口,则源将rx间窗口内处于活动状态。如果省略f4 f5,则源将在时间窗口开始时启动,在时间窗口
结束时停止。
str1:为与源一起使用的波形的标识符。
例如,一个y方向,内阻为50欧姆,幅度为5,中心频率为1.2 GHz高斯波形的电压源。则指令为:
#waveform: gaussian 5 1.2e9my_gauss_pulse
#voltage_source: y 0.05 0.05 0.05 50 my_gauss_pulse.
6.6 transmission_line:命令
功能:在电场位置引入一维传输线模型。传输线的指定电阻可以大于零,小于自由空间的阻
抗(376.73欧姆)。当模型中包含天线的物理特性时,对天线的激励是有用的。
命令格式为:#transmission_line: c1 f1 f2 f3 f4 str1 [f5 f6]
c1:为传输线的极化,可以是x、y或z。
f1、f2、f3:为模型中传输线的坐标(x、y、z)
f4为:传输线源特性电阻,单位为欧姆。可以是大于零且小于自由空间阻抗(376.73欧姆)的值。
f5、f6:为可选参数。f5是启动传输线励磁的时间延迟。f6为消除传输线励磁的时间。如果时间窗
长于输电线拆除时间的激励,则在输电线拆除时间的激励之后,输电线的激励将停止。如果传输线
的激励解除时间长于时间窗,则传输线的激励将在整个时间窗内处于活动状态。如果省略f5 f6,则
传输线的励磁将在时间窗开始时开始,在时间窗结束时停止。
str1:为与源一起使用的波形的标识符。
传输线中电压和电流值的时间历史记录保存到输出文件中。这些都记录在输出文件部分中。
这些参数对于计算天线的特性有用。
例如,输入阻抗或s参数。gprMax包含一个Python模块(在工具包中)来帮助您查看使用传输线
馈送的天线模型的输入阻抗和s11参数。
例如,一个z向传输线源,其电阻为75欧姆,幅度为5,中心频率为1.2 GHz高斯波形。则指令
为:
#waveform: gaussian 5 1.2e9 my_gauss_pulse
#transmission_line: z 0.05 0.05 0.05 75 my_gauss_pulse
6.7 rx命令
功能:在模型中引入输出点。在这些位置上,电场和磁场分量的值在模型迭代次数上将被保
存到文件中。
命令格式为:#rx: f1 f2 f3 [str1 str2]
f1、f2、f3:为模型中接收机的坐标(x,y,z)。
str1:为接收方的标识符。
str2:为接收器的输出列表。它可以是Ex、Ey、Ez、Hx、Hy、Hz、Ix、Iy或Iz中的任何选项。
注:当没有可选参数str1和str2时,所有电场和磁场分量将随接收点输出。
6.8 rx_array命令
功能:在模型中定义多个输出点。
命令格式为:#rx_array: f1 f2 f3 f4 f5 f6 f7 f8 f9
f1、f2、f3:为输出线/面/体积的左下(x,y,z)坐标,
f4、f5、f6:为输出线/面/体积的右上(x,y,z)坐标。
f7、f8、f9:为(x,y,z)增量,它定义了每个方向上输出点的数量。
f7、f8、f9可以设置为零,以防止在特定方向上出现任何输出点。否则,f7的最小值为∆x,f8
的最小值为∆y,f9的最小值为∆z。
6.9 src_steps和rx_steps命令
功能:在模型运行之间移动所有简单源(#src_steps)或所有接收器(#rx_steps)的位置。
命令的语法为:#src_steps: f1 f2 f3
命令的语法为:#rx_steps: f1 f2 f3
f1、f2、f3:为增量(x,y,z)来移动所有简单源(#hertzian_dipole或#magnetic_dipole)或所有接收器
(使用#rx或#rx_array命令创建)。
6.10 geometry_view命令
功能:输出关于模型几何的文件信息。
命令格式为:#geometry_view: f1 f2 f3 f4 f5 f6 f7 f8 f9 file1 c1
f1、f2、f3:为几何视图体积的左下角(x、y、z)坐标,以米为单位。
f4、f5、f6:为几何视图体积的右上角(x,y,z)坐标,以米为单位。
f7、f8、f9:为几何视图的空间离散,以米为单位。
通常这些参数与模型的空间离散化相同,也可以改得更小。
file1:为文件的文件名,几何视图将存储在与输入文件相同的目录中。
c1:n (normal)或f (fine),它指定是以每个单元为基础(n)还是以每个单元边缘为基)输出几何信息。
fine模式应该保留用于查看占用小体积的几何结构的详细部分,因为使用此模式可以生成具有大文
件大小的几何文件。
注:只检查模型的几何形状时,可以使用可选的命令行参数--geometry-only运行gprMax。构
建模型并生成任何几何视图文件,但不会运行模拟。
6.11 snapshot命令
功能:在给定的时间瞬间获得关于模型体积内电磁场的信息。
命令格式为:#snapshot: f1 f2 f3 f4 f5 f6 f7 f8 f9 f10 file1
#snapshot: f1 f2 f3 f4 f5 f6 f7 f8 f9 i1 file1
f1、f2、f3:为快照体积的左下角(x,y,z)坐标,单位为米。
f4、f5、f6:为快照体积的右上角(x,y,z)坐标,单位为米。
f7、f8、f9:为快照的空间离散,单位为米。
f10、i1:为以秒为单位的时间(浮点数)或迭代数(整数),表示快照将被拍摄的时间点。
file1:为要保存快照的文件名。快照文件会自动存储在输入文件名后面加上' _snaps '的目录中。对于多个模型运行,每个模型运行将有自己的目录,即' _snaps1 ', ' snaps2 '等…
例如,在模拟时间为3纳秒时保存模型中电磁场的快照。则指令为:
#snapshot: 0 0 0 1 1 1 0.1 0.1 0.1 3e-9 snap1
可以利用Python脚本创建一系列快照。
例如,要以0.1ns的间隔从0.1ns开始到3ns创建30个快照,请在输入文件中使用以下代码片
段。相应替换x1 y1 z1 x2 y2 z2 dx dy dz
#python:
for i in range(1, 31):
print('#snapshot: x1 y1 z1 x2 y2 z2 dx dy dz {} snapshot{}'.format((i/10)*1e-9,i))
#end_python: