毕业设计跟聚焦超声相关,主要是论文复现,做一个仿真和一个实验,顺便记录一下过程。
建模操作
新建模型
物理场选择压力声学和电流,研究选择稳态。
全局定义
参数定义
阶跃函数
单位设置
在几何1的设置窗口,将单位设置为mm
模型绘制
换能器
换能器的设计需要考虑发射信号的频率和功率;超声波传播介质的衰减和吸收;当然还有换能器本身的位置和尺寸。
- 近场距离 N,计算公式为: N = D 2 f / 4 c N=D^2f/4c N=D2f/4c。D是换能器直径,f是频率,c是介质中声速。
- 焦距F,即换能器与相当于目标区的聚焦点之间的距离。
- 场深或聚焦区,它表示 -6dB 信号的幅度与最高幅度之间的差距,计算公式为: F Z = 2 F 2 N + F 2 F_Z=\frac{2F^2}{N+\frac{F}{2}} FZ=N+2F2F2
论文中使用半径35mm,焦距70mm的碗状换能器,即换能器元件的曲率半径等于焦距70mm。由于碗状换能器无法直接绘制,我们可以用矩形和圆取交集。
首先绘制一个半径70mm的圆,z轴坐标70mm,绘制一个矩形,宽35mm,高9.38mm。然后使用布尔操作和分割,选择交集。
这里矩形的高需要用勾股定理来计算,已知圆的半径和弦长,求出矩形的高为9.38mm。
生物组织
接下来绘制生物组织,注意这里需要绘制完美匹配层(PML)。
绘制一个宽48.6mm,高75.5mm,z轴坐标9.38mm的矩形,并在右侧和顶面绘制5mm的层。
如果想要焦区位于组织的正中心,z轴坐标的计算需要考虑焦距和组织的长度,
z
=
70
−
35
=
35
m
m
z=70-35=35mm
z=70−35=35mm。
但是论文中的组织和换能器是直接贴合在一起的,因此z轴坐标应该是换能器的厚度
z
=
9.38
m
m
z=9.38mm
z=9.38mm。此时,焦区中心坐标
z
=
35
+
9.38
=
44.38
m
m
z=35+9.38=44.38mm
z=35+9.38=44.38mm。
椭圆形聚焦区
绘制一个a轴7.5mm,b轴1.5mm,扇形角180°,旋转270°的椭圆,中心坐标 z = 35 + 9.38 = 44.38 m m z=35+9.38=44.38mm z=35+9.38=44.38mm。
形成联合体
完美匹配层
在组件->定义中加入完美匹配层。
几何类型选择圆柱形,坐标拉伸选择有理数。
压力声学
指定水的物理属性
从声压级参考压力列表中选择使用水的参考压力,
C
r
e
f
C_{ref}
Cref输入1483[m/s]。
压力声学1
在压力声学模型栏中,选择用户定义的衰减,衰减系数输入alpha_water。
压力声学2
在压力声学模型栏中,选择生物组织包含的域,并选择用户定义的衰减,衰减系数输入alpha_tissue。
添加材料
水域
在内置材料中选择Water,liquid。
生物组织
新建空材料,选择生物组织所包含的域,并输入对应参数(相对介电常数未知)。