开题准备

流程手续

文献调研

1. Codalike Multiple Scattering of Elastic Waves in Dense Granular Media

2. Probing intermittency and reversibility in a dense granular suspension under shear using multiply scattered ultrasound

3. Ultrasound Propagation in Externally Stressed Granular Media

4. Monte Carlo Simulations of Ultrasound Scattering and Absorption in Finite-Size Heterogeneous Materials

5. Time reversal of ultrasound in granular media

6. Transitional nonlinear elastic behaviour in dense granular media

工作日志

10 月

2023.10.13

完成了声学探测装置的组装, 通过网线连接到了实验室的电脑里. 初步掌握了配套软件的使用方法. 由于最终目的是要将声学探测和 CT 探测结合起来, 所以两者的同步率也非常重要. 幸运的是, 清远公司在这个设备上做了时间戳处理, 时间同步是通过网络时钟确定的. 但是有关 CT 设备在时间戳方面的设置还不清楚(准确地说, 是对 CT 装置本身都还不清楚!), 以后组装的时候需要对这方面进行考虑.

初步拟定了一个工作计划:

2023.10.18

完成了任务书和开题报告的交稿, 估计后面还会继续和毕设的助教在开题报告的格式上面扯皮(很奇怪, 在 Latex 已经广泛应用的今天, 为什么不直接做一个模板, 大家直接往里面填内容就完事了呢?). 和师兄交流了一下, 下一步应该是设计和搭建实验装置, 我现在有两件事情需要做:

2023.10.25

今天师兄告诉我声学功率放大器到了实验室, 我准备今天把装置取到实验室, 结合之前的超声波探测器一起测试一下装置搭建的可能性.

评估内容你自己设定，可以参考一下下面这些：

1. 被动接收模式：实验室环境噪声的强度水平是多少分贝，人说话的声音是什么水平，颗粒物质运动的声音是什么水平，可以用采集软件自带的功能做FFT，看看频谱特征频率怎么样等等。

2. 主动接收模式：用两个传感器相对放置，一个从信号发生器接收信号，另一个采集。可以看看输入信号分别为正弦波，方波，脉冲波等等的时候，通过空气介质或者颗粒介质接收到的信号分别是什么波形、两个传感器表面间距对于采集信号的影响、使用功率放大器与否的影响等等。

因为是试用品所以只能用到 10.31, 所以需要充分利用这段时间. 注意保存配套零件以便到时候再寄还.

和老师的讨论后的结论是先做一些简单的实验, 所以设计环切装置的事情可以先搁置. 但是, 学习 SolidWorks 使用, 设计方形或者圆筒形的颗粒容器也是需要继续的.

重新连接了 RAEM 软件, 里面有关 IP 地址的设定让人头昏脑胀, 后续会写一个专门的教程. 如果明白有关 TCP 的知识应该会好很多, 不至于像今晚这样盲人摸象一样匹配, 但是物理系的学生应当秉承 “拿来主义”, 不求甚解; 能用就行.

明年会和实验室里的师兄协商一下信号发生器的使用时间. 许家兆师兄申请的样机其实连接还算比较简单, 关键是我对于电路知之甚少, 只知道信号放大可以是 $0-20\text{dB}$, 步长为 $2\text{dB}$, 至于 $\text{Rin}$ 和 $\text{Rout}$ 应该是和信号源的内阻有关, 但是更细一点的内幕就不太清楚了.

今天的收获主要还是在自学了怎么从 0 开始配置超声探测装置, 前面师兄说的是只要会用就行, 但是今天到了实验室发现之前的配置全部都失效了, 所以也是看着说明书一步步配置完了. 还有就是实验室的电脑只有 1080P, 看起来有些糊, 趁手程度当然还是比不上自己的电脑. 也许后面会和师兄一样弄个向日葵连接.

2023.10.30

今天参与了组会, 虽然自己的实验进度不堪入目, 但是还是得到了很多建设性的意见.

我的大致思路是, 方盒中的颗粒会产生背景辐射, 并且还会产生非常明显的特征频率. 并且, 人声, 摇动颗粒等行为均不会对该特征频率产生任何影响.

但是有一个问题: 没有考虑在空气中进行裸测, 无法排除是不是实验室本身就有着这种特征频率的信号. 为了避免这种可能, 需要补测空气中的信号.

另外, 邢义师兄表示, 空气其实是超声波的不良导体, 应该考虑在水或者其它介质中进行实验.

在许家兆师兄提前看过我的展示 PPT 并且听过我的讲述之后, 给出了新的实验思路:

一、被动式：

1. 不放在颗粒中测量直接环境噪声（43dB会不会有点高了，我记得上次设置40dB的门限的时候是接收不到信号的）。
2. 测量搅动颗粒的声信号（即把探头放在装有盒子里和颗粒直接接触，然后用棍子搅动颗粒模拟剪切行为，观察接收到的信号变化）

二、主动式：

1. 评估对信号的转换和接收：信号发生器用固有频率150Hz输入不同波形（正弦波、方波、脉冲波）的信号，两个探头直接接触相对放置（放在自带的夹具里磁吸在一起），看接收信号的波形是否和输入一致；
2. 评估测量的频率范围：选定一种波形，将输入信号频率从60kHz到400kHz范围内调节，观察接收信号的变化（强度和波形）；
3. 评估空气的衰减：选定一种波形，将两个探头的距离从0到20cm之间等距调节，观察接收到的信号的变化（主要看强度）；
4. 功率放大器测试：两个探头直接接触，看功放的放大倍率和放大效果。以及使用功放与否，能探测到信号的传感器之间的极限距离的差别。给出功放是否必要使用的结论。
5. 静态颗粒加压模拟测试：用一个小袋子装一些颗粒，用夹具隔着颗粒磁吸在一起模拟加压情况。可以竖直放置，这样可以在上面压金属块调节压强。可以作为预实验。参数你自己根据颗粒尺寸设置。给出声学装置测量颗粒声学是否好用的结论。

今天还在罗茹丹师姐的帮助下熟悉了 1 号楼 3 楼的 3D 打印实验室, 初步学会了颗粒的制备流程, 后续应该会写一份颗粒制备的教程. 这里先简要记录一下:

在电脑上导入颗粒的模型(这个需要提前在 SolidWorks 上面做好), 利用复制操作做好颗粒阵列, 注意在颗粒之间留出一定的间距.

维护 3d 打印机. 注意支撑材料(SUP)和打印材料(ABS)的提前更换. 瓶盖需要轻微扭松以确保气压平衡. 注意清洁陶瓷板(打印平台), 打印机本身非常容易积灰. 注意清理废料袋.

打印完成后, 放入冰箱进行冷冻, 等待成型后分别取出含有支撑材料的颗粒板, 并且在这个过程中完成陶瓷板的清洁.

将颗粒取入放到具有较大尺寸和孔径的筛网容器中, 放置在加热干燥室中, 等待支撑材料熔化流出. 等待完成后转移到具有较小尺寸和孔径的筛网容器中.

超声清洗. 可以提前烧好开水减少加热时间. 注意控制温度为 $58-62$ ℃, 原则上不应该超过 $60$ ℃, 但是温度绝不能低于 $58$ ℃.

清洗的具体过程将是明天中午的时间讲解.

2023.10.31

粗糙记录一下听到的有关清洗步骤的内容:

在将颗粒转移到大孔径大尺寸筛网上并放置于烘烤箱中时, 可以在此时一并清洗陶瓷板. 此时可以使用无水乙醇将陶瓷板上的所有有机材料一并洗去.

提前用热水壶烧开水, 配合冷水混合配置出 $58-62$ ℃ 的水温.

第一次和第二次的清洗需要加入洗洁精, 用量约为 $5\text{ml}$, 计算在按压上约为 $15-18$ 次. 倒入时可能会有部分液体粘于容器壁, 可以使用超声波清洗机里的热水荡洗溶入.

注意干燥箱里的熔化出的支撑材料需要及时清理, 避免其溢出烘烤箱.

超声清洗装置加热的设定温度为 $64$ ℃, 但是实际操作过程中应当配合使用酒精温度计.

倒入提前配好的热水, 注意水线约为筛网容器容深减去约 $3\text{cm}$. 戴好手套, 用手将取自干燥箱的颗粒尽快放入已经盛好水的筛网中, 开启超声进行清洗.

一般设定的一轮清洗需要进行约 $10$ 次玻璃棒搅拌. 可以配合记录笔记提醒自己完成完整次数的操作.

如果清洗过程中水的温度下降较快, 可以将盖合上. 在这个过程中可以继续使用热水壶烧制开水, 以备下一轮的清洗使用. 一般是 $2-3$ 壶热水倒入热水桶中并且倒入冷水中和, 过程中可以使用酒精温度计测量温度.

等待 $30$ 分钟后, 一轮清洗完成, 开闸放水, 将洗涤残余支撑材料的废水导入红色水桶中. 注意此时的温度仍然比较高, 可以选择将其放置变凉倒入洗手池或者开冷水倒入.

注意放水过程中, 需要注意不要让颗粒露出液面接触空气(具体原因未知), 接近露出时闭闸, 然后按照上述的洗洁精加入等步骤进行第二轮清洗.

等待下午进行后续清洗步骤的讲解.

下午因为临时调度没有办法帮上清洗颗粒的忙

托赵西梅老师借了 3 号楼的声学装置:

信号发生器 * 1, 示波器 * 1, 声速测量仪器(两个压电陶瓷做的换能器, 上面有刻度尺) * 1, T 型 BNC 接口, BNC 线材 * 4(接口上三根外加额外的一根)

算是重温了一下本科物理实验, 声速测量有

1. 共振法

共振频率可以通过调节发射频率, 观察接收到的信号是否出现极大值, 然后按照此法依次调节发射信号的十位, 个位等数位的具体数值, 由此确定发射频率 $f$.

因为两个压电陶瓷实际上是等同的, 所以接收端的探测器在接收到发射过来的信号时也会反射信号至发射端, 当两者间距为半波长的整数倍时就会形成共振响应, 体现在接收端上就是形成信号峰的极大值.

所以记录下峰为极大值的位置, 使用逐差法就可以确认(半)波长 $\lambda$.

根据以上数据就可以得到声速的测量值 $v = f \lambda$.

2. 相位法

将输入信号作为 X 轴, 接受信号作为 Y 轴, 在示波器上就会出现各种李萨如图形. 形成一三象限的直线即为同相, 二四象限的直线即为反相. 记录出现(反)直线的位置, 使用逐差法计算得到(半)波长 $\lambda$. 计算出声速的方法与上述公式相同.

3. 时差法

我猜测论文里经常提到的 Time of Flight (TOF) 就是这个方法. 之前还在使用 RAEM1 的时候我就在考虑怎么确定发射时间, 现在有了一个 T 头的 BNC 接线就了然了.

要测定声速还是要用其良导体才好, 这里在导师千叮咛万嘱咐之后在实验桌上加了水. 设定一个合适频率的脉冲波, 当然还要注意设定示波器合适的 Scale, 要让发射的脉冲信号(比如黄色信号)和接受信号(比如紫色/粉色信号)同屏出现, 使用 Cursor 就可以轻松计算出两者的时间差(一般几十到一百多微秒).

这里还要注意的是, 这样的时间差 $\Delta t_{i}$ 和探头距离 $L_{i}$ 不能直接相除, 因为还需要考虑声电/电声效应的固有时间, 所以正确的做法是拟合两者出现的直线, 然后将其斜率作为采用值.

$$ L = \color{red}{v}\times (\Delta t - t_{0}) $$

11 月

2023.11.1

昨天晚上的实验相当于是用空气和水做了一次练手, 现在需要考虑的是:

1. 如果是在颗粒物质中会怎样?

从这里就需要涉及弹性力学相关的知识了. 加上前面的在学 Edwards 统计方法, 我觉得我也应该安排时间来学习一下热力学的知识了, 虽然书是已经买了放到了工位, 但是除了预录取面试前夕的抱佛脚, 目前因为杂七杂八的事情都还没有好好学过.

如果要真正研究应力等因素对有效介质的影响, 没有系统学过弹性介质肯定是不行的.

另外, 原本用作空气和水的容器比较大, 用于颗粒物质的话, $5\text{mm}$ 所需要的数量会很大, 肯定不是借用许家兆师兄的备用颗粒能够做到的, 怎么借调到这么多的颗粒也是一个问题.

2. 如何将 RAEM1 的这一套装置实现出声学仪器类似的效果?

广州的这个声学装置的主要优点可能也就在可以方便地做一个 FFT 进行频率分析, 但是说实话可能没有什么大用. 比如前面提到的 TOF, 如何实现在 RAEM1 上还是一个悬而未决的问题.

另外, 尚不清楚 RAEM1 和示波器之间的适配性, 理论上来说只是将声学装置中的压电陶瓷探头换为了 RAEM1 的探头, 但是能否实现, 如何实现等问题都亟待解决.

2023.11.03

今天稍微整理了一下这一周以来的实验情况, 也和导师交流了一下, 发现了之前做实验没有考虑到的一些问题. 现在趁着还有记忆在这里陈列一下:

$\sqrt{}$ 凡是都要考虑空白对照, 比如之前摇晃颗粒的时候确实记录了对应过程的信号, 但是没有考虑到 无任何处理 的颗粒介质中的信号, 所以这个实验回头需要补测.

$\sqrt{}$ 不要太过依赖软件的计算功能(比如 FFT), 尽可能自己写好脚本来处理数据.

$\sqrt{}$ 和技术再次确认如何稳定使用连续采集而非包络采集;

? 为什么 RAEM1 和示波器上观察到的图像不一样. e.g. 正弦波和带有锯齿的正弦波?

? 为什么 RAEM1 上没有观察到对应的倍频效应?

$\sqrt{}$ 确认单位 dB 的计算方式, 由此确定信号之间的差异倍数是如何得到的.

$\sqrt{}$ G150 探头的含义是共振频率为 $150\text{ kHz}$, 但是为什么实验测得的共振频率为 $129.2\text{ kHz}$? 是不是漏测数据了?

$\sqrt{}$ 和技术再次确认探头的信号发射的方向性的问题, 理论上超声信号在空气中的传播应当方向性很好.

$\sqrt{}$ 没有考虑到信号发生器的共振频率的问题, 找相关人员或者资料确认对此频率数值是否有研究的必要, 或者还是说可以忽略这个问题.

$\sqrt{}$ 用相对粗糙的压强施加方法来制备一个致密堆积然后测量, 至少要先观察到和 Jia 论文里相似的波形, 即一个相对快, 幅值大的相干弹道波和一条相对慢, 幅值小, 广谱的散射波.

接下来的一段时间应该是处理开题答辩的事情. 我已经受够繁文缛节了.jpg

2023.11.06

完成了开题答辩, 因为太经典力学(不涉及任何量子效应), 太实验(不涉及任何模拟计算), 答辩专家们反而无计可施了, 除了郑杭老师问了问实验室的配置, 没有其他问题可问.

9 月的时候和导师面谈, 说是 “比较独立的方向”, 现在我信了.

完成了这些形式主义的东西之后应该就是继续回到实验, 准备把上面的几个问题一一解决. 另外今天下单了网口拓展坞, 明天就能到货, 希望可以顺利迁移程序到新电脑上.

1. dB 的定义方式

首先需要明确的是 以 $10$ 为底, 并且是功率相比.

因此我们也就可以理解为

$$ N_{\text{dB}} = 10\lg{\frac{p_{i}}{p_{o}}} $$

所以功率放大器增益 $x\text{dB}$ 就是指的

$$ N_{\text{dB}} = 10\lg{\frac{U_{i}^{2}/R_{i}}{U_{o}^{2}/R_{o}}} $$

2. Matlab 实现 FFT

目前已知 Origin Pro 也具有非常方便的 FFT 功能. 但是为了方便地进行数据分析, 最后还是要落实到自己写代码上.

另外确认了一下, RAEM1 的采样频率应当是 $2\text{ MHz}$, 也就是 $2000\text{ kHz}$. 而借来的示波器的采样频率是 $1\text{ GSa/s}=1\text{ GHz}$. 这意味着, 我们的探测设备很有可能难以满足后续的实验需求.

但是我查阅了一下 Stick-Slip 的文章, 里面的信号采样频率也就是 $2\text{ MHz}$, 明天应该再去找导师去确认一下.

clear;clc;

%%% This code was created to read the dB.csv created by
%%% RAEM1. 

%%% 1. Read the data in the *.csv file

data = readmatrix("dB.csv");
[Length, Width] = size(data);
Amp = data(4:Length, 1);
num = Length - 3;
Freq = num / 5 * 10^3;
delta_t = 1/Freq;

tspan = (0:num - 1) / Freq;
fspan = (0:num - 1) * (Freq / num);

%%% 2. Draw the time-domain signal

figure(1)
subplot(2,1,1)
plot(tspan*10^3,Amp)
xlabel("Time(ms)")
ylabel("Amplitude(V)")
subplot(2,1,2)

ref_Vol = 1e-5;
%%% RAEM1 的分贝值是如何计算出来的?
Amp_dB = 20 * log10(abs(Amp)/ref_Vol);
plot(tspan*10^3,Amp_dB)
xlabel("Time(ms)")
ylabel("Amplitude(dB)")

%%% 3. FFT. 
%%% Remember that the sampling rate is 2 MHz, so the time step
%%% equals to 1 us.

DFT = abs(fft(Amp));
half_num = floor(num / 2);
figure(2)
plot(fspan(1:half_num)/10^3, DFT(1:half_num))
xlabel("Frequency(kHz)");
ylabel("Amp");

3. 信号发生器的共振频率

如果探头的共振频率和信号发生器的压电陶瓷共振频率不同, 确实可能会影响到实验结果. 但是信号发生器是否具有这样一个特质: 当我们选择到一个频率 $f_{0}$ (比如 $129.2\text{ kHz}$)的时候发射, 信号发生器会将自己的共振频率选定为 $f_{0}$?

还有一个问题, 信号发生器本身到底有没有所谓的一个压电陶瓷存在? 我怀疑其本身只是调制了电流而并没有一个真正的压电陶瓷.

查阅了一下, 并不是所有的信号发生器都是使用的压电陶瓷作为振荡器, 这取决于具体的设计和需求. 明天应该看看这个信号发生器的型号(顺便也把实验室里其余的信号发生器信号也记一下), 看下都各自是什么工作原理, 再做进一步打算.

压电控制器 & 信号发生器.

4. 超声辐射的方向性

超声波甚至可以用来测速. 这是不是说明, 这个 G150 探头产生的信号强度其实并不太高? 否则为什么只是隔了几厘米, 信号强度就衰减到和背景噪音强度差不多的程度了.

超声波的频率越高, 方向性会越尖锐, 但随之而来的是衰减也越大.

2023.11.07

今天不仅完成了软件的迁移, 还学会了如何操纵 RAEM1 的连续采集方式. 而且, 许家兆师兄带来的金属片很好地充当了临时活塞头的效果, 把各种金属(师兄们加工过的金属块, 还有一大袋的金属球珠), 形成的致密颗粒堆积的传声效果确实比之前的松散堆积好多了.

今天通过一些简单的实验也发现了不少问题, 现在趁着记忆还在简单梳理一下:

1. 连续采集虽然声称是 0 间隔进行区间采集, 但是细看 “全波形” 图像会发现实际上还是有空隙的. 如何确保所需的数据完整的被采集到一个区间内将是一个需要研究的问题.
2. 关于设置信号频率的问题. 首先需要明白: 为什么要测共振频率? 在测水的声速进行练手的时候, 需要注意到我们实际上并不是使用的水中的共振频率, 而是简单的设置了 $100\text{ kHz}$, 然后观察接收信号出现的位置和时间.

实验装置的需求整理:

1. 能够方便地安装和卸下 G150 探头, 并且能够通过滑尺等方式测量探头之间的距离.

2. 能够方便地控制施加应力, 目前的容器施加压强的方式比较定性, 无法精确控制.

还有关于 RAEM1 和示波器之间的协调关系. 示波器是不是本身在外部有 BNC 线材接入的情况下就会出现一个误导性的电压抬升和波动?

具体来说, 就是把频率调为 $100\text{ Hz}$ 的时候, 就可以看到类正弦的接收信号在 $t = 0$ 处会有一个相对高峰, 当时还以为真的测到信号了. 然而, 无论是将信号发生器的输出关掉, 还是干脆直接把信号发生器上面的 BNC 接口直接拔下, 这个疑似相干弹道波反而还更清晰了(在没有任何信号被激发的情况下, 接收端竟然接收到了信号!)这无疑是非常荒谬的.

因为在参考的论文中设定的脉冲是 $2\mu\text{s}$, 我理解的是这个数值就是指的脉冲的宽度(Width), 所以频率就在 $10-1000\text{Hz}$ 里面选择, 这种思路和测水中声速是一脉相承的. 并且, 确实可以在紫色信号看到有意义的峰, 而非散乱的信号杂音. 但是, 并没有看到所谓的相干弹道波, 这一点对于我们的散射理论反而是更重要的.

目前的话, 我认为应该重读一下当时参考的声学探测论文, 还有就是以周志刚的这篇博士论文为纲学习颗粒物质. 之前给的夏成杰师兄的文章固然好, 可惜并不是从声学角度给出的, 如果我要从声学入手颗粒物质检测, 那么跟着周的文章学自然是最轻松的.

昨天的开题报告, 实际上已经是部分同学的结题报告, 不由得感叹人生的参差. 悟已往之不谏, 知来者之可追, 在这一年抓紧时间熟悉课题组的工作内容, 更应该是现在应该去上心的事情.

11.13

记录一下周四和导师交流产生的悬而未决的问题.

1. 如何解决 RAEM1 只有单通道的问题. (涉及测量时间差从而测定声速的问题)

2. 实验室的装置, 应当先斩后奏, 买了再报销; 不能一直被动地等.

3. 重新调研信号和颗粒的参数, 需要将对应的参数匹配到自己的装置中.

4. 如何观察到相干弹道波的现象.

目前我的研究进展的痛点:

1. 急需一个能够满足自己研究需求的装置, 这个装置需要:

• 回顾 Jia 和 Zhou 的论文, 总结他们的装置参数和对应的实验图像, 从而做出自己的实验装置设计和图像预期;
• 使用 SolidWorks 进行设计;

2. 声学和通用的颗粒知识学习应当齐头并进, 不可偏废.

记录一下新的装置到了之后产生的效果:

1. 信号发生器能够产生 N-Cycle 的信号, 基本上可以用于复现 Jia 的脉冲信号;
2. 示波器比起之前借用的示波器更高档一些, 但是仍然是噪音波动还会有残留效应, 导致在观察波形(包含波前起点)的时候比较困难, 需要用到信号放大器把发生器的信号再次放大在可以看到比较明显的现象.

不难观察到, 随着循环数的增加, 接收信号得到的幅值也会更高. 但是, 我们不能无限制地增加循环数, 因为这样就会导致信号的重合. Jia 文章设定为 $10$, 但是在实际的操作过程中, 可以将循环数设定到 $15$ 左右.

通过示波器的 Cursor 获取时间差(~ $200\text{ us}$), 再通过刻度尺测量两个探头之间的距离(~ $5\text{ cm}$), 最后计算出的声速约为 $200-300\text{ m/s}$, 和 Zhou 论文里的数量级大致相仿, 具体的差异可能还要看颗粒参数和边界条件.

我的想法是, RAEM1 的单通道性能其实比较废物, 我目前也没有想到什么除了再买一个 RAEM1 然后接入到 SWAE 以外的方法. 显然这是很浪费的; 并且我有理由认为, 既然 RAEM1 的采集性能也就是 $2\text{ MHz}$, 那么 SWAE 的数据处理性能上限可能也就是在这个数量级(因为从单次采集, 即 $5\text{ ms}$, 得到的数据数量刚好就是根据 $2\text{MHz}$ 来的), 所以我认为另一条路径就是, 设法将示波器和主机连接, 直接输出其波形数据. 但是这样肯定就比不上 SWAE 处理数据的方便程度了(主要是即时性, 虽然一些数据功能可以通过自己编写 Matlab 来实现, 但是终究是比较麻烦的)

另外, 目前是已经可以得到波在颗粒介质里的频谱展宽现象, 但是还没有掌握具体的分析方法(因为需要具体情况具体分析, 不能直接引用 Jia 的结论), 所以最后还是要落回到细读原文的研究方法上.

今天准备复现之前的 10 循环发生信号, 结果发现这个信号发生器比想象的要难用很多, 一直都很难调出和之前一样的效果. 导师这边又在催进度, 好像都默认周末就是用来加班加点做实验的. 虽然说的确在声学实验上投入了不少钱(新的发生器, 示波器还有专门的电脑), 但是这些都是工作上的投入, 和计算组的同学相比我觉得我都算是勤奋的了. 新的装置又不可能一夜之间就做出来了, 前面的技术和知识铺垫哪一个都是需要花费时间的. 现在如果要说尽快有进展, 那就只能换颗粒材料, 看下玻璃珠和钢珠都是什么情形, 但是这个本质上是没有在技术上获得新进步的, 所有的测量操作都是一样的. 总结下来就是一句话: 欲速则不达.

11.14

今天又额外举行了一次组会, 虽然说是定的下午两点半, 但是实际上还是延误了十来分钟才开始的. 今天开的组会也是重量级, 一次开了接近五个小时. 下了楼四餐早关门了, 就去全家买了个三明治, 感觉还不错.

先简单记录一下今天参与组会的感受.

1. 要避免像王鑫师兄那样, 一次性考虑太多的因素. 比如, 在做粮仓效应的时候考虑尺度, 摩擦等众多因素. 导师的意见是: 我们设计实验是为了解决问题而不是提出更多的问题, 要成为实验学者而非理论学者.

2. 积极参与到组会中的思维训练中, 导师随时都有可能抽奖到你来回答一个科研思路问题. 这就要求听组会的过程中要精神集中, 清楚各个师兄师姐都在干什么.

再简单说一下今天做的实验. 今天主要是考虑了循环数对于响应信号的影响, 目前结果是已经出来了, 但是如何对结果进行分析, 仍然不得要领. 另外, 速度时已经可以测出来了, 应该是比 ABS 的塑料小球要快一些.

昨天做了一个小的固定时间测试实验, 即把两个探头使用夹具相吸固定测时差, 结果显示的确存在约有 $1\mu\text{s}$ 的时间差, 而测速度的时候的时间差也是~ $200\mu\text{s}$, 虽然说相比之下确实小很多, 但是将其考虑进去也会更加科学.

缺少信号放大器带来的效果也很快显现了, 虽然信号发生器已经尽可能调到最大幅值(10 Vpp), 但是反应在接收信号上, 对比度就显得有点小了.

下一步的打算是:

1. 想办法初步控制压强(应力);

2. 学会使用信号发生器的任意波形功能, 模仿判断飞行时间选取的小实验.

11.15

论文重读和问题研究:

1. Sound pulse broadening in stressed granular media

为了获得相干压缩波在颗粒介质中的脉冲响应，我们可以用源传感器发送的信号 S 对通过介质传输的信号进行解卷积(deconvolve). 这里的解卷积的原理和用途是什么?

该文章中的发射信号并不是之前设定的方波脉冲, 而是自己设定的波形.

飞行时间的定义方法涉及参考点的选取: 峰值*10%($t_{0}$), 峰值($t_{1}$), 首次过零($t_{2}$).

参考信号是在参考介质中测定的. 实验发现, 用 $t_{1}$ 测得的速度几乎与探头间距无关. 在正式实验中, 显然这也是一个值得探究的点, 可以稍微做一点数据出来说明选取的合理性.

可以通过重标定(rescale)检查响应信号, 即 $A/A_{1}, t/t_{1}$.

还有极其复杂的一维展宽推导和数值模拟

2. Ultrasound Propagation in Externally Stressed Granular Media

目前已经能够掌握粗糙的飞行时间速度测量, 如何从简单的波形上确定群速度? 如何统一两者在颗粒介质中的描述思想?

除了飞行时间速度以外, 还有周志刚提到的反射声波法, 多普勒效应法, 共振柱法, 后面应当再独自调研一下, 了解为主

颗粒的尺寸将会显著影响相干波的形状. 可见, 在颗粒尺度较大的情况下信号被严重散射了, 应当考虑使用尺寸更小的颗粒, 这样观察到的相干波应当会产生一个对比度极高的峰.

注意要区分两种相干波, $E$ 波和 $R$ 波. 后者是前者的反射波.

• 信号参数: $500\text{kHz}$ 为中心的宽带频谱($20-1000\text{ kHz}$), 持续时间 $2\text{ ms}$, 振子振幅为 $10\text{ nm}$(需要用干涉仪才能够测定).

• 颗粒参数: 都是多分散, 分为三组: $0.2-0.3, 0.4-0.8, 1.5\pm 0.15\text{ mm}$

飞行时间确定是首次达到峰值的 $5$%. 有效介质速度为 $v_{\text{eff}} = 1070\pm 30\text{m/s}$.

注意: 探测器尺寸也会影响到信号, 减小尺寸会增强散射波(S)的振幅. S 对构型高度敏感而 E 则复现程度高, 用 $\Delta P/P=(P’-P)/P$ 来评估这一参数. 这说明 E 是一个自平均信号. E 信号才是定义声速的信号.

声学共振频率的确定: $\nu\approx v_{\text{eff}}/d$.

3. Codalike Multiple Scattering of Elastic Waves in Dense Granular Media

需要关心力链关联长度和颗粒介质中的波长之间的关系. 波长减小到晶粒大小的时候, 力链产生的散射效应将会非常显著.

• 信号参数: $20\mu\text{s}$, 10-Cycle, 中心频率 $500\text{ kHz}$.

• 颗粒参数: $kd = \omega d/v\approx 2.9, v = 750\text{ m/s}$. 应力 $P = 750\text{ kPa}$, 厚度为 $L = 11.4\text{ mm}$. 直径 $d = 0.6-0.8\text{ mm}$.

• 容器参数: 直径 $D = 30\text{ mm}$, 高度为 $7-15\text{ mm}$.

$E_{p}$ 幅值较小.

需要考虑: 如何做到的去除噪声? 如何做到系综平均?

传输的时域信号是相干分量(低频)+剧烈散射的非相干分量(高频), 能够在最新的示波器上观察到吗?

注意: 相干波存在两种, 一种是压缩波($E_{P}$), 一种是剪切波($E_{S}$), 只有 $E_{p}$ 才是自平均的, 所以它对构型不敏感; 后者使用大尺度探头更容易探测到.

为了对散射波进行研究, 可以通过数字滤波去除低频的相干波.

如何测定ABS球和钢珠的共振频率?

4. Ultrasound propagation in disordered granular media

相干弹道波表现的是有效介质特征, 杂乱的散射波表现的是介质不均匀力场.

• 颗粒参数: $400\mu\text{m}-5\text{ mm}$. 实例数据为 $0.4-0.8\text{ mm}$.

• 容器参数: $30\text{ mm}$,

• 信号参数: 以 $500\text{ kHz}$ 为中心的一个或多个周期的输入脉冲. 持续时间为 $2\mu\text{s}$ 的单周期脉冲波, 传播距离约为 $L = 11\text{ mm}$.

表现出明显的相干弹道波特征. 随着颗粒尺度增大, E 和 S 都表现出低频特性.

$\lambda = v_{\text{eff}}/f_{c}$ 与颗粒尺度 $d$ 相当时, 将出现频率截断现象(低通特性).

如何研究多重散射声波的传输特性? 使用 10-Cycle 脉冲, $f = 500\text{ kHz}$.

11.16

今天听了邢义师兄的博士学位预答辩, 足足有三个小时之长. 可见如果没有充分的理解和研究, 面对老师的质询答不出来的话, 场面将会极其尴尬. 所以功夫在平时, 做好每一步的调研和思考, 想常人之难想, 方能天衣无缝, 滴水不漏.

在听的过程中也记录了一下问题, 目前尚不能理解:

1. Edwards 系综解释中, 我的理解是态数目为体积固定为 $V$ 时可能的构型数目, 这个计数应当是不涉及具体的外部驱动过程的, 为什么说需要将循环剪切和振动两种驱动类型分开看待?

2. 如何考虑微重力下, 外部施加等重力压强, 这个可以等价于普通重力场下的颗粒体系吗? 我认为是不可以的, 因为需要考虑粮仓效应或者说 Jassen Effect, 按照反作用力去理解, 均匀外加压强显然不是同一回事.

另外在发呆过程中也思考了几个问题:

1. 如何去除空气中的超声杂声? 为什么会存在这样的信号? 这些莫名其妙的峰可以通过使用简单的挡板来去除吗? 还是说需要用上数字滤波器?

2. 为什么自己的 FFT 算法得出的结果和 SWAE 算出来的结果差这么多? ($\text{Amp%}$)

3. 根据简单的波方程知识, 可以知道波速, 频率, 波长之间存在关系:

$$ \lambda = v/f $$

前面可以粗测得到等效声速 $v_{\text{eff}} \approx 220\text{ m/s}$, 所以如果是声波被剧烈散射, 就应当取颗粒直径 $d = 5\text{ mm}$ 作为波长, 计算得到的频率就是 $44\text{ kHz}$; 如果要观察到明显的相干弹道波, 就应该让频率再降低, 使得长波极限被采用, 也就是~ $4.4\text{ kHz}$.

问题在于, 感觉测定的声速有点偏小, 不知道是不是就是这么多. 可以回去再测测钢珠的, 感觉数量级都差不多.

首先是, 有关于飞行时间参考点选取的问题, 还需要设计实验来验证; 其次, 虽然 Jia 声称参考信号就是用的 “中心频率为 $500\text{kHz}$”, 但是并不清楚这个是怎么做到的, 反而是那个参考点选取的实验, 虽然波形看上去有点奇怪, 但是至少是可预计的未来时间内可以自己设计波形, 然后使用任意波形发生器激发做到的.

另外, 虽然 Jia 设定的信号幅值很小($100\text{ mV}$), 但是传播距离也很小($11\text{ mm}$). 我想, 目前涉及的不只是颗粒材质的问题, 更应该考虑频率和颗粒尺度之间的匹配关系. 目前颗粒尺度还有待考量, 但是频率尚且是可以自己控制的, 可以先做一点实验对比一下.

11.20

今天进行了一次巨长无比的组会, 虽然今天没有抽问环节, 但是仍然感觉身心俱疲(14:40-18:30, 19:30-23:00, 接近七个半小时), 主要还是最后的哲学分享环节十分折磨, 但愿下次没有这种事情发生了.

说回正事, 今天主要是讲了一下之前的 PPT 的总结, 导师没有给出什么新的建议, 主要还是要和贾小平老师及时联系, 沟通一下装置可能存在的问题. 另外就是要把三篇文件看了, 一方面是再次熟悉一下大课题的方法处理, 另一方面就是学习一下 Jia 的实验处理方法.

11.21

今天解决了一个工程学上的难点, 即在颗粒层数较多或者颗粒直径较小的情况下, 很难在已经倒入颗粒之后重新把探头埋进去. 现在想起来可以用双面胶将其固定在容器壁, 这样就可以很轻松的重复制备颗粒堆积而无需考虑探头位置的变化. 唯一需要考虑的可能就是多出来的双面胶也会固定一些颗粒从而造成浪费.

今天在没有改变容器构型的情况下, 多用了两种颗粒堆积, 但是仍然没有观察到明显的弹道波. 难道是对 $\lambda$ 的计算有误解吗?

现在重新整理一下实验做下来仍然没有搞清楚的点:

1. 经常提到的 “以 $500\text{kHz}$ 为中心的宽带频谱” 到底是怎么实现的?

2. 外部施加压强的影响应该需要超过重力带来的影响才可以. 但是实际上, 即使没有加上外加重物, 也足以观察到现象. 需要设计新的容器;

3. 施加压强实际上会极高(~$\text{MPa}$), 所以面积就不能太大, 所以设想的容器不能够把接线长度引入, 而是必须要挖孔;

4. 等效波长 $\lambda_{\text{eff}} = v_{\text{eff}}/\nu$ 里面的 $\nu$ 到底是指的什么?

5. 颗粒尺度已经来到了 $2\text{ mm}$, 再尝试把距离减小试试看?

6. 新的容器制作必须提上日程.

11.27

上周五的时候老师没有把关键问题问出来, 今天组会的时候才问了几个没答上来的问题. 我也发现了之前实验中存在的一些思路上的问题.

1. 声学探测的原理是什么? 我们关心什么物理量的演化?

2. Jia 为什么要选定一个 $500\text{kHz}$ 为中心频率? 它和颗粒的直径有着怎样的联系?

3. 不要把目光局限在静态堆积上. 尝试回答声学能够解决什么问题(关心波形, 能量)? Jia 的声学方法是什么, 解决了什么问题?

更多的信息可能需要明天看完贾老师的五小时演讲才会知道怎么做.

让我不满的是, 我本身就是以复现论文结果为主的, 让我突然成为一个声学专家是怎么回事? 还说所谓考试和他无关, 那我怎么应对明年的一系列 GPA 课程还有博资考呢?

感觉有点不太现实, 突然把一个完全很大的领域交给我一个人, 一周的时间来精通, 搞笑呢

12 月

12.05

昨天简单的介绍一下自己对于贾老师研究的理解情况, 但是邢义师兄对于我的未来学习安排有些担心. 仔细想想自己对于未来比较宏观的学习安排还是有些不太清晰, 昨天说的测定响应函数也只不过是和老师交流之后产生的新想法, 但是我对于接下来应该怎么做, 具体去做些什么, 其实并没有真正的想法.

在这里简单列举一下应该做的事情:

1. $\sqrt{}$ 熟悉探头的响应情况, 规划响应函数的测定方案;

2. $\sqrt{}$ 学习 声学原理 的声强计原理, 以及非线性声学的处理思路;

3. $\sqrt{}$ 自己写一个滤波函数出来, 学会如何使用 high pass 和 low pass 函数, 关于频率的截断是如何做到的.

4. 继续学习 SolidWorks 操作, 制作一个简单的圆筒壁容器, 方便直接使用贾老师用的扩散函数.

5. 系统学习声学的扩散理论, 清楚这个方程背后的具体原理;

6. 继续学习贾老师的暑校演讲, 争取这周把后面的 $2/3$ 看了;

7. $\sqrt{}$ 地震学笔记, 继续阅读学习;

8. $\sqrt{}$ 问技术探头的制造原理; 学习压电效应, 搞清楚是否在线性区域内进行工作.

另外吐槽一句, 昨天收到了这样一份邮件:

你好。

由于你在报考24级直博生预录取阶段，填写的意向导师王宇杰老师，近期即将离职交大，所以你2024年9月入学后，需要更换一名意向导师，并根据学院政策，在你的第一学年结束时，再确定正式导师。

特此通知，收到本邮件还请回复一下，以便我确认你已知晓上述信息。

如有疑问，请随时与我联系，谢谢。

合着就我一个人不知道老师都要离职了. 在实验室呆了两个月, 直接人都没了可还行. 总之后面少不了多的麻烦.

今天花了不少时间在个人网页的重新整理上, 实验上的安排估计不多, 所以还是放在阅读和学习上. 明天起个早班来把实验做一下.

12.06

今天花了不少时间在声强计原理和非线性声学上, 估计这一周就可以学完. 今天看累了, 用剩余的时间写了一个滤波函数, 检验了一下还不错. 看到频域上有个 $130\text{ kHz}$ 的峰就感觉火大, 明天可以做一点粗糙而对实验检验一下自己的一些想法.

顺便今天收到了新的回信:

你好。

经请示主管领导，你的疑问答复如下：

1、王老师调离交大，等你的第一学年结束时需要正式确定导师的时候，王老师不会在导师库里，你只能在导师库里选择导师。至于你博士工作是不是还是和王老师继续合作，由你选择的导师和王老师之间确定。

2、第一学年结束时，选择的导师就是你的正式导师，该导师需要符合所有的关于导师招生人数的限额要求。

3、如（1）所述，王老师调离交大，选择正式导师时候，你无法选择王老师作为系统里面的正式导师。

特此回复，望知晓，谢谢。

总之就是很麻, 不知道去找谁, 问题就在于这个傻逼学校只能让你选有名额的导师, 那当时预推免的时候又拿出一个名额, 有点搞笑了. 为什么要让学生来承担后果呢?

12.11

今天测定了脉冲波和白噪声下的响应效率函数, 果然效果很差. 但是客观来说, 这是探头的工作范围本来就在 $60-400\text{ kHz}$ 所致.

下一步的学习安排:

1. 学习贾老师暑校 (2/3), 学习剪切带如何用声波探测;
2. 自己写一个 Correlation 函数判定;

$$ R(i,j) = \frac{C_{i,j}}{\sqrt{C_{i,i}\cdot C_{j,j}}},\quad C_{i,j}(t) = \int_{0}^{T}S_{i}(t)\cdot S_{j}(t)\mathrm{d}t $$

3. 整理 huwei 的实验仪器参数, 提出仪器需求;
4. 尝试拟合自己的实验曲线, 得到扩散特征时间 $\tau_{\alpha}$;
5. *搞清楚扩散方程的原理;
6. *寻找一下专业的校准超声探头的技术手册, 学习其校准方法;
7. 重新设计自己的曲线包络算法, 得到更理想的数据, 尝试得到衰减特征时间 $\tau_{\alpha}$;
8. 重新学习如何测定相速度 $v_{\varphi}$.

12.12

尝试写出关联函数. 比较奇怪的是, 离散积分得到的值是 $[0, T]$ 上整个区间的, 应该只有一个值, 为什么能够得到一段 $R(i,j)$ 的时域分布呢? 没有搞懂是这个是如何定义出来的.

前面响应函数的测定, 在用连续波激励的时候存在一个问题, 就是只考虑了倍数关系, 没有考虑声阻抗实际上是一个复数, 所以真正计算的时候是需要处理相位信息的. 但是实际上, 我们可以做一个简单的心算:

$$ ka = \frac{2\pi}{\lambda_{E}}\frac{d}{2} = \frac{\pi f_{c}}{v_{E}} =\frac{\pi\cdot 2\times 10^{-2}\cdot 129\times 10^3}{10^3}\approx 7\gg 1 $$

而声阻抗表达式

$$ Z_{S} \approx \rho_{0}c_{0}S\frac{(ka)^{2}}{2} + i\rho_{0}c_{0}S\frac{8ka}{3\pi} $$

所以应当采取短波极限, 也就是大约为 $2$ 左右. 相位畸变几乎可以忽略不计. 感觉可以不用深究有关于声阻抗造成的相位畸变.

感觉要做的事情一时间太多, 稍微整理一下思路.

12.18

工位办公室竟然养了一只猫.

12.19

昨天没有开组会, 倒是留出了一天的有效时间来阅读文献. 讲真, 如果组会频率下降, 就不用每周都要单独分配一天来做 PPT, 效率还提高了不少.

现在简单整理一下这周的任务计划:

12.20

今天开了周一没有开的组会, 这里总结一下除了上面的清单要做的事情:

12.21

询问了一下赵西梅老师有关于声学校正的看法, 没有得到太多有用的信息.

比较麻烦的是 BNC 三通转接头也被还回去了, 现在只能紧急订购几个, 在到货之前没办法做实验了.

12.25

1. $\sqrt{}$ 重构有关非线性响应的代码, 要求识别并且排除信噪比较低的数据, 自动从合格数据开始分析.

现在的代码已经可以自主选择非线性分析起点了, 可以说是智能程度大大增加了. 另外, 需要注意的是, 虽然在做实验的时候可以把 $\delta V_{\text{input}}$ 取小一点, 但是这对于改善信号的信噪比, 也就是说有效数据起点没有任何用处. 所以真正需要的是那些在小振幅就可以很高效地传播信号, 自己用的超精细颗粒/精细颗粒无意间让弹性模量变得很小, 这显然是不利于实验正常进行的.

也就是说竟然无意间达到了 “开局即巅峰” 的境界, 玻璃珠 + 90khz 效果最好. 既然如此, 可以再尝试一下 2mm 铁质颗粒或者其他颗粒. 或者加大一下压强, 效果应该也会明显一些.

明天把买的转接头拿到实验室重新搭一下线路, 然后就是仔细读一下有关声学润滑和诱导雪崩的文章, 之前看了但是没有看懂, 看懂了就可以设计新一轮实验了.

12.26

今天拿到了转接头, 一看发现型号根本不对, 怪不得那么便宜, 可谓是便宜没好货.

1. 重新购买转接头. 非线性效应只能先将就着做;

2. 联系许家兆师兄看如何购买信号放大器;

3. 重新阅读润滑和雪崩的文章, 考虑如何设计参数;

12.27

今天阅读了高时间分辨率的文章, 里面谈及了有关主滑移和微滑移的频谱相似性. 直觉告诉我, 这东西可能和时域相关性是等价的, 但是还没有使用傅里叶变换进行严格证明.

记录一下明天要做的事情:

12.28

我们自己的探头只能检测到压缩波吗? 能够检测到剪切波吗?

现在也算是做了一点实验, 也学了一点声学方法, 现在简单总结一下. 这里就省去引言在内的无关紧要的废话, 只谈文章脉络和发展思路. 已经做过的就简单一些, 还没有做过的就补充实验设计的方法和思路.

12.29

和老师进行了交流.

1. 调研各个厂家的探头

2. 询问技术有没有更硬件的控制方式

3. 更深入地阅读有关探测剪切带的问题, jia 是如何做出关于剪切带从无到有产生的推测

4. 询问代学亮师兄的剪切盒子, 配置好有关剪切力传感器等物件, 尝试进行同步性分析

5. 咨询一下技术有关于探头的有效探查频域在物理上是如何实现的

和王鑫师兄聊天后得到建议与经验

1. 参考课表进行提前选课(软物质, 博资考需要注意高量和热统)

2. 多找袁后非师兄, 邢义师兄, 曾志坤师兄交流, 询问一下想法是否可行, 看下有什么参考意见. 不要偷懒

3. 优先推进自己的项目, 有时间再看看小组学习的知识;

4. 推荐阅读 granular media 和 组里的基金本子, 师兄师姐的毕业论文也是值得参考的对象.

研究致密堆积的时候, 需要有效利用相似度定义波前进行 时间平均. 这在处理非线性效应的谐波分量时也非常好用.

高斯轮廓拟合, 得到波形图. 这样处理的意义是什么? 因为有效信号是低频波所以需要进行高斯模糊, 去除杂乱的高频成分?

1 月

1.3

研究一下有关于传感器在和剪切盒耦合的时候应该怎么处理. 那个耦合剂能起到什么作用? 老师的意见是可能会起到低通的作用, 实际上会是什么作用?

1.4

今天在代学亮师兄的帮助下熟悉了剪切盒的使用, 也成功做出一条应力-应变曲线, 但是目前还存在不少问题.

首先记录一下实验仪器的参数

1. 电机

最大活动范围是 $150\text{ mm}$, 剪切速度通过 $\text{pps}$ 定义. 粗测了一下 $10\text{ pps}$ 对应的剪切速度 $\approx 37\mu\text{m/s}$.

贾老师的剪切速度可以达到 $0.6\mu\text{m/s}$, 这还是相对于颗粒尺度为 $\sim 700\mu\text{m}$. 我们相对剪切速度为贾老师的 $10$ 倍左右.

可以使用 jog 或者保存的配置文件来调整电机的参数. 一个移动单位对应 $0.1\text{ mm}$.

2. 力学传感器

最快的采样频率为 $1/\tau = 1/0.0005 = 2000\text{ Hz}$. 虽然这个采样频率已经算是比较高, 但是对于 $10\text{ MHz}$ 还是相形见绌. 师兄还告诉我, 这个采样频率可能对应的时间间隔相对于现实世界的时间流逝而言未必准确. 即使如此, 在放大之后观察其曲线, 仍然会有一定的粗糙感.

3. 剪切盒

目前的问题:

• 缺少较短的连杆&套筒, 还需要采购一下多通道采集卡. 三头 BNC 转接头也需要尽快就位;

• 无法观察高度变换. 可能需要类似于限位器的装置, 能够保持表面平整的同时还可以具有高度测量的作用.

• 目前仍然欠缺长时间稳定的声学测量手段.

1.5

恍然大悟. 高斯调制的原理原来如此简单.

回到调制信号上, 假设中心频率为 $f = 150\text{ kHz}$, 定义

$$ S(t) = \sin{(2\pi ft)}\cdot G(t) = \sin{(2\pi ft)}\cdot \frac{1}{\sqrt{2\pi\sigma^{2}}}e^{-\frac{(t-\mu)^{2}}{2\sigma^{2}}} $$

$\mu$ 一般与实际实验中的循环数有关, 即一般都是取其一半所对应的时间

$S(t)$ 进行傅里叶变换:

$$ S(\omega) = \int_{-\infty}^{+\infty}S(t)e^{-i\omega t}\mathrm{d}t $$

记高斯函数 $G(t)$ 的傅里叶变换为 $G(\omega)$, 正弦函数的傅里叶变换为 $\delta(\omega - 2\pi f) + \delta(\omega + 2\pi f)$. 根据卷积定理有

$$ S(\omega) = G(\omega)\cdot(\delta(\omega - 2\pi f) + \delta(\omega + 2\pi f)) = G(\omega - 2\pi f) + G(\omega + 2\pi f) $$

这就不难理解了: 在高斯分布上乘 $f$ 的单色波, 等效于将高斯分布在频域上移动至 $2\pi f$ 处.

目前的两种方案:

1. 监测法. 记录整个剪切过程 (约 $1000\text{ s}$) 的声学信号. 目前还在探索如何使用 RAEM1 做到这一点.

2. 探测法. 等待后续剪切盒的设计和使用. 可能需要涉及对容器顶/底进行挖孔测试.

以 $1\text{ Hz}$ 或者 $2\text{ Hz}$ 频率发射试探信号, 记录这些信号.

计划:

1. 力学信号处理

• 如何在数学上定义滞滑事件, 对滞滑事件进行标记与统计;
• 记录事件之间的时间间隔, 并且对其进行统计;
• $\sqrt{}$ 学习时间窗口配合加权平均(考虑分别使用抛物线权重和高斯权重)对曲线进行平滑处理;

只是看不太出来各个方法的优劣是什么. 只有一点是确定的: 力学传感器的力学精度只有 $0.1\text{ N}$, 如果更精细一些就好了.

• $\sqrt{}$ 双轴法绘制实验图像. $x$ 轴: 时间/位移(Displacement); $y$ 轴: 应力/相似性参数;
• $\sqrt{}$ 学习如何对局部图像进行放大的绘图方法.

2. 连续声学采集方法

• 再次尝试使用 RAEM1 采集规定时间的信号, 并且保证高可重复性

为什么设定采集 $10\text{ s}$ 的数据, 表格生成出来却对应了 $25\text{ s}$ 或者 $7\text{ s}$ 的数据点?

3. 物资采购

• 连杆和套筒
• 多通道采集卡
• 三母头的 BNC 信号转接头

虽然现在 $\Gamma_{i,i+1}(t)$ 还没有办法做出来, 但是现在作图工具已经备齐. 虽然这个放大工具好用, 可惜对于双 $y$ 轴的放大效果还是不行, 需要后续手动调整. 目前把四个轴都占用上, 反而触发 bug 使其放大失败失败(没错, 失败两次!), 从而做出还不错的效果.

后续把 $\Gamma_{i,i+1}(t)$ 做出来之后, 这种示意图可能就需要自己手动写代码控制了.

1.7

记录一下组会上的一些质疑和建议

1. 目前对于 micro- 的定义仍然不算清晰, $0.1\text{ N}$ 可能包含了较多的噪声, 使得记录的 micro 事件数量偏多; 后续应该使用另一种直方图绘制方法, 将这些 $0.1\text{ N}$ 的信号标记出来.

2. 阅读学习老师推荐的 Sounds of Failure: Passive Acoustic Measurements of Excited Vibrational Modes.

3. 思考如何构造更好的剪切盒装置;

4. 测试更大的剪切形变($\sim 4\text{ cm}$)

5. 记录一下剪切盒的尺寸参数.

1.8

当探头铺设到剪切盒底部时, 没有记录到有效的声学信号.

忘记先和师兄商量剪切盒能不能开洞就已经执行了, 真是脑子短路了.

PS. 后面还得知误拆了卢海洋师兄的剪切加压装置, 也是非常抱歉. 在开学后的组会上私下再次道了歉, 希望之后没有再次事件.

3 月

3.3

简单写一下可能的仪器采购清单

先抄写一下参考文章所用的 NI 仪器.

3.4

恢复手感用的实验. 加入放大器的因素.

最大加入放大器后的不畸变电压值: 150 V. 可能的问题在于: 无法比较自由的调整放大电压间隔.

SURUGA DSCONTROL-WIN

Force Recorder Standard

软件迁移到了自己的电脑上, 但是因为缺线所以还没有实际做出电机的效果.

3.7

把装置的图纸画出来了, 并且已经和制作装置的师傅联系上了, 比较蛋疼的一点是真正制作的时候才知道 板材实际上只有 $5,6,8,10,12\text{ mm}$ 的尺寸. 还好师傅说不用重新画图, 只需要稍微改一下厚度就可以投入生产了.

还有需要注意的一点是, 开 M4 螺纹孔实际上并不需要将螺纹孔打通, 只需要在对应位置留出一定的空间就可以了. 公差在画图的时候也是需要考虑的, 师傅帮忙考虑的是活塞留出半径 $\pm 0.25\text{ mm}$, 支撑柱留出半径 $\pm 0.5\text{ mm}$. 到此就算是完成 DGM 装置的设计和交付了, 明天就会收到报价.

另外还和导师讨论了有关多通道采集卡的事情, 最后得到的共识就是贵的实际上不是采集卡(BNC-2110)而是机箱(PXle1084), 机箱将数据上传到电脑上, 重量级价格(58000+), 而老师听说的价格是十万, 这一点我尚不清楚到底是怎么回事, 应该再具体整理一下各文献中所使用仪器的具体价格. 在 NI 方面, 除了索要报价外, 最好是看看能不能申请到样机; 声学探头则是主要看一下我们自己用的几千块钱的探头和 Micro-II Express-8 的探头性能差距到底在哪里.

简单写一下要具体做的实验系列

3.18

1. 测量 face-to-face time

peak-to-peak: 14.2 $\mu$ s wave front-to-wave front: 4.08 $\mu$ s

3.19

face-2-face time, 信号参数: 100kHz, 10Vpp, 4循环

2. 考察不同参考点(设第一个峰值为 $A_{1}$, 分别取 $A_{1}/10, A_{1}, 0$) 对应的事件为参考点, 研究对应声速的变化, 确定最佳声速参考点选取

参考文章: Sound pulse broadening in stressed granular media, 10.1103/PhysRevE.91.022205

实验条件

装置自重: 1031.2 g

容器参数: 内径 9cm, 容器限高 12.5 cm

顶部活塞自重: $\sim 260$ g.

参考文章不同距离取点: 10

第一组实验, 使用 5mm + 6mm 双分散钢珠, 测试声学信号是否明显

第二组实验, 使用 5mm + 6mm 双分散玻璃珠, 可以观察到显著信号.

偶然理解了 Tektronix 示波器的时间平均功能, 得到了无法想象的完美曲线, 我也算是理解到贾小平老师那么干净的图像究竟是怎么得到的了.

那么现在再根据示波器的双通道同步采集, 再加上事后对于波形信号文件的处理, 得到该实验想要的曲线将会易如反掌.

设定示波器的平均次数: 128

发射信号参数设定: 30 kHz, 5 循环 10 Vpp(10dB 增益), 发射间隔 5ms 通过 TekScope 软件导出了示波器的数据, 目前得到的都是静态数据, 不知道后续会不会有什么多通道采集卡的替代方法.

受压的玻璃珠堆积对应的飞行时间法测声速图像.

使用重物的质量: 5.4kg 对应的压强:

$$ \frac{Mg}{\pi r^{2}} = \frac{(5.4+0.26)\times 9.8}{3.14\times (4.5\times 10^{-2})^{2}} = 8719\text{ Pa} $$

通过呕心沥血写出来的波前和波峰搜索算法投入使用, 最后得到了不同参考点随距离变化的图像.

由此可知在较短约化长度的情况下, 声速测量误差较大, 排除一开始的误差后, 得到的图像和贾小平教授的结论基本相近.

更详细的数据可能需要换用更小更多的颗粒, 并且需要配备相应的前置放大器.

3.24

3. 声速与外部应力的关系

第一组.

2mm 钢珠 - 10 mm.

信号参数: 30kHz-5循环正弦波列, 10Vpp(+10dB), 触发间隔 10ms.

依次称量放上去的金属块质量, 每次的应力是之前放上的金属块质量之和

然而，人们观察到，在堵塞球体堆积中，低 $P$ 值时声速为 $V\propto P^{1/4}$，高 $P$ 值时声速为 $V\propto P^{1/6}$（我们将证明本研究中的实验仍处于 $V\propto P^{1/6}$ 阶段）

$P^{\alpha}: P\uparrow, \alpha\downarrow, p = 1/\text{slope} = \frac{1}{\alpha}\uparrow$

在自己的实验中, 施加的压力应当始终都是较小的.

3.26

注意多次测量取平均值, 并且加上标准差计算而来的误差棒

4. 相速度测量

令传播信号为

$$ s(x,t) = A_{0}e^{-\alpha x}e^{i(\omega x/{V_{\varphi}(\omega)} - \omega t)} $$

频域上的衰减函数 $\alpha(\omega)$

$$ \alpha(f) = \frac{1}{x_{2} - x_{1}}\ln{\frac{S_{1}(f)}{S_{2}(f)}} $$

频域上的相速度

$$ V_{\varphi}(f) = 2\pi f\frac{\Delta x}{\varphi_{2}(f) - \varphi_{1}(f)} $$

之前为了测量飞行时间存的文件实际上并不太适用于计算相速度.

波形文件需要采集的时间长度尽可能完整(包括相干波和散射波); 设计算法, 只选取相干波而筛去散射波;

使用 到达时间-距离 进行拟合, 这种方法比较接近传统的声速测量方法.

Elastic wave propagation in confined granular systems PHYSICAL REVIEW E 72, 021301 (2005)

模仿实验的数据:

可以看到虽然选取的参考点不同, 但是在距离上的传播是线性的, 归根到底是相干波中具有不同的频率成分, 不同的频率具有不同的波速. 为了更准确全面地描述这个 Sound Pulse Broadening in Granular Media 的过程, 最好的方法仍然是测量信号的相速度.

$$ \frac{400\text{ m/s}}{30\text{ kHz}} \approx 13.3\text{ mm} \sim 6d $$

基本可以使用长波近似.

这种速度和飞行时间速度(T.O.F. Velocity) 的比较:

可见在测量距离大于一定情况下可以使用飞行时间法进行近似

模仿绘制声脉冲展宽的图像, 以及学习处理其归一化

注意这些声波重叠处理的时候, 归一处理的类型不同: 同一距离/厚度, 不同应力/压强(这里为了测相速度在两个距离 $x_{1}$ 和 $x_{2}$ 都测了, 但是施加的应力是相同的) 同一距离/厚度, 不同应力/压强

模仿的 Rescale 图像

• 时间和振幅均归一化

• 时间对齐于峰, 振幅归一化

使用归一化宽度来描述在逐渐变厚的介质中的扩散行为

$$ W\sim C_{W}\left(\frac{L}{d}\right)^{-\frac{1}{2}} $$

进行的模仿实验表明出的指数非常接近($-0.47493\approx -0.5$)

需要写出一个适用于示波器的响应函数修正. 示波器的采样频率可以达到很高($\gg 2\text{ MHz}$), 按照 fft 得到的频域来一个一个试响应振幅显然是不现实的.

有一点是可以明确的: 关心的频域一般都是在 $[0,1000]\text{ kHz}$. 有一定的原因是之前的采集频率是 $2\text{ MHz}$, 但是也要根据实际的探头性能来, 目前几乎没有什么探头能够在高于 1000 kHz 的频域还能有较好的响应的.

L = []? M = [1,2,3,4,5,6,7,8,9] ensemble_size = 7

有关 TekScope 的使用小技巧

将网线连接示波器和电脑之后, 启动 TekScope, 然后再开机示波器, 此时示波器会根据 DHCP 自动分配 ip 地址, 将其填入 Host Address 即可秒连.

有关示波器和信号放大器的使用小技巧

制备颗粒堆积系综的时候, 由于使用的平均采样, 示波器 CH1 的源信号有时候会迅速衰弱. 此时调整一下信号放大器的 Gain 值即可令其回正.

相速度测量

相速度测量很难准确的原因在于

$$ \Delta\varphi(\omega) = \phi(\omega) + 2N(\omega)\pi $$

显然, 这里的 $N(\omega)$ 是很有可能彼此之间毫无关系的, 我们也并不能通过简单的加减 $2\pi$ 重置区间来得到正确的数值. 上面的文章中介绍了一种基于互相关方法粗略估计声速的方法. 关键词是 “连续小波变换”

理解归一化宽度 $W = (t_{1} - t_{0}) / t_{1}$ 的物理含义.

傅里叶变换的两个特性

$$ g(t) = f(t - t_{1}) \Rightarrow \widetilde{g}(\omega) = e^{i\omega t_{1}}\widetilde{f}(\omega) $$

$$ g(t) = f(at) \Rightarrow \widetilde{g}(\omega) = \frac{1}{|a|}\widetilde{f}\left(\frac{\omega}{a}\right) $$

$$ \begin{aligned} a(x,t) &= \frac{1}{2\pi}\int_{-\infty}^{+\infty}\widetilde{a}(\omega)e^{-i\omega t}\mathrm{d}\omega = \frac{1}{2\pi}\int_{-\infty}^{+\infty}A(\omega)e^{ik(\omega)x}e^{-i\omega t}\mathrm{d}\omega\\ &= \frac{1}{2\pi}\int_{-\infty}^{+\infty}A(\omega)e^{i[k^{\prime}(\omega) + ik^{\prime\prime}(\omega)]x}e^{-i\omega t}\mathrm{d}\omega \end{aligned} $$

对于类 Dirac 的脉冲波, 有 $A(\omega) = \delta(\omega-\omega_{0})$

1. 一维质量球链($k^{\prime\prime} = 0$)

对于一维质量弹簧链, 我们已经知道色散关系

$$ \omega^{2} = \frac{4K_{s}}{M}\sin^{2}{\frac{ka}{2}} $$

使用长波极限得到波速度

$$ \omega^{2} = \frac{K_{s}}{M}k^{2}a^{2} \Rightarrow V = \frac{K_{s}}{M}a $$

具体到一维质量球, 即有 $a = 2R$, 从而得到原文中的色散关系式

$$ k = \frac{1}{R}\sin^{-1}\left(\frac{\omega R}{V}\right) $$

其中已经知道 $\arcsin{x}$ 的麦克劳林级数

$$ \sin^{-1}{x} = x + \frac{1}{2}\frac{x^{3}}{3} + \frac{1\cdot 3}{2\cdot 4}\frac{x^{5}}{5} + \cdots $$

采取长波极限 $V\rightarrow \infty$

$$ k = \frac{\omega}{V} + \frac{1}{6}\left(\frac{\omega}{V}\right)^{3}R^2 $$

所以波数项 $\text{exp}(ikx)$:

$$ \text{exp }i\left[\frac{\omega}{V}x + \frac{1}{6}\left(\frac{\omega}{V}\right)^{3}R^{2}x\right] $$

使用 $t_{1} = \frac{x}{V}$ 和 $\omega_{1} = \left(\frac{2V^3}{R^{2}x}\right)^{\frac{1}{3}}$ 简化表达式:

$$ e^{i\omega t_{1}}e^{i\left[\frac{1}{3}\left(\frac{\omega}{\omega_{1}}\right)^{3}\right]} $$

第一项可以利用傅里叶变换的时间平移性质处理. 第二项对应于 Airy 函数傅里叶变换.

$$ \mathcal{F}(\text{Ai})(\omega) = \frac{1}{2\pi}\int_{-\infty}^{+\infty}\text{Ai}(t)e^{- i\omega t}\mathrm{d}t = e^{i\frac{1}{3}(\omega)^{3}} $$

只是使用了 $\omega_{1}$ 进行了 Rescale 操作.

所以信号

$$ a(x,t) = \omega_{1}\text{Ai}[\omega_{1}(t_{1} - t)] = \omega_{1}\text{Ai}\left[\omega_{1} t_{1}\frac{t_{1} - t}{t_{1}}\right] $$

而 $(t_{1} - t_{0})w_{1}\approx \frac{\pi}{2}$

$$ W = \frac{t_{1} - t_{0}}{t_{1}} \sim \frac{1}{\omega_{1}t_{1}} = \left(\frac{R^{2}}{2x^{2}}\right)^{\frac{1}{3}}\propto x^{-\frac{2}{3}} $$

取 $x = L$ 即得到我们熟悉的结论. 显然 $-\frac{2}{3}$ 离 $-\frac{1}{2}$ 还是有一定的差距的, 这是因为没有考虑 $k^{\prime\prime}(\omega)$ 导致的.

2. 频域上的衰减 $\alpha(\omega)$

是不是可以使用求相速度的方法, 能量/功率耗散图来求解这个问题?

存在衰减

$$ \alpha = k^{\prime\prime}(\omega) = \left[\sigma_{K}^{2}L_{c}^{n-1}\right]\left(\frac{\omega}{V}\right)^{n} $$

时域信号的傅里叶变换为

$$ \widetilde{a}(\omega) = e^{ik(\omega)x} = e^{ik^{\prime }x}e^{-k^{\prime\prime }x} = e^{i\omega t_{1}}e^{-\left(\frac{\omega}{\omega_{1}}\right)^{n}} $$

其中

$$ t_{1} = \frac{x}{V},\quad \omega_{1} = \left(\frac{V^{n}}{\sigma_{K}^{2}L_{c}^{n-1}x}\right)^{\frac{1}{n}} $$

令 $s(t)$ 是 $\text{exp }(-\omega^{n})$ 的逆傅里叶变换. 则时间信号写作

$$ a(x,t) = \omega_{1}s[\omega_{1}(t-t_{1})] $$

所以

$$ W\sim\frac{1}{\omega_{1}t_{1}} = (\sigma_{K})^{\frac{2}{n}}\left(\frac{L_{c}}{x}\right)^{\frac{n-1}{n}} $$

在一维随机分层介质中, $n=2$; 代入 $x = L$ 即有

$$ W\propto L^{-\frac{1}{2}} $$

3. 更一般的情形(3D, $n$)

$$ \alpha\sim\sigma_{K}^{2}L_{c}^{n-1}k^{n}\sim\omega^{n} $$

频散能量图公式推导.

$$ \Gamma_{i,j}^{t}(\tau) = \int_{-\infty}^{+\infty}S_{i}(t)S_{j}(t+\tau)\mathrm{d}t $$

$$ \begin{aligned} \mathcal{F}[\Gamma_{i,j}^{t}(\tau)] &= \int_{-\infty}^{+\infty}e^{-i\omega\tau}\int_{-\infty}^{+\infty}S_{i}(t)S_{j}(t+\tau)\mathrm{d}t\mathrm{d}\tau\\ &= \int_{-\infty}^{+\infty}S_{i}(t)\int_{-\infty}^{+\infty}S_{j}(\tau+t)e^{-i\omega\tau}\mathrm{d}\tau\mathrm{d}t\\ &\overset{u = \tau + t}{=} \int_{-\infty}^{+\infty}S_{i}(t)\int_{-\infty}^{+\infty}S_{j}(u)e^{-i\omega(u-t)}\mathrm{d}u\mathrm{d}t\\ &= \int_{-\infty}^{+\infty}S_{i}(t)e^{i\omega t}\int_{-\infty}^{+\infty}S_{j}(u)e^{-i\omega u}\mathrm{d}u\mathrm{d}t\\ &=\mathcal{F}[S_{i}(t)]^{\star}\mathcal{F}[S_{j}(t+\tau)] \end{aligned} $$

傅里叶变换的平移性质:

$$ \mathcal{F}[S_{j}(t+\tau)] = e^{i\omega\tau}\mathcal{F}[S_{j}(t)] $$

所以有

$$ \mathcal{F}[\Gamma_{i,j}^{t}(\tau)] = \mathcal{F}[S_{i}(t)]^{\star}e^{i\omega\tau}\mathcal{F}[S_{j}(t)] = \mathcal{F}[S_{i}(t)]^{\star}\mathcal{F}[S_{j}(t)]e^{i\omega t} $$

或者

$$ \mathcal{F}[\Gamma_{i,j}^{t}(\tau)] = S_{i}(\omega)^{\star}S_{j}(\omega)e^{i\omega t} $$

具体在编程上, 就是取 $\delta(\omega_{n})$ 与 $\mathcal{F}[\Gamma_{i,j}^{t}(\tau)]$ 相乘, 然后进行逆傅里叶变换得到时域分布

$$ \mathcal{F}^{-1}\left\{\delta(\omega_{n})\mathcal{F}[\Gamma_{i,j}^{t}(\tau)]\right\} = \mathcal{F}^{-1}\left\{S_{i}(\omega_{n})^{\star}S_{j}(\omega_{n})e^{i\omega_{n} t}\delta(\omega_{n})\right\} $$

我们已经知道, 频域上的一个 $\delta(\omega)$ 函数, 对应于时域上的三角函数:

$$ \begin{aligned} \mathcal{F}^{-1}\left\{\delta(\omega_{n})\mathcal{F}[\Gamma_{i,j}^{t}(\tau)]\right\} &= \int_{-\infty}^{+\infty}S_{i}(\omega_{n})^{\star}S_{j}(\omega_{n})e^{i\omega_{n} t}\delta(\omega_{n})e^{i\omega\tau}\mathrm{d}\omega\\ &= S_{i}(\omega_{n})^{\star}S_{j}(\omega_{n})e^{i\omega_{n} (t+\tau)} \end{aligned} $$

这表明当

$$ \omega_{n}(\tau+t) = 2N\pi\Rightarrow \tau = \frac{2N\pi}{\omega_{n}} - t $$

时, 信号 $S_{j}(t)$ 角频率为 $\omega_{n}$ 的分量通过平移时间 $\tau$, 可以与 $S_{i}(t)$ 角频率为 $\omega_{n}$ 的分量完全重合(在功率谱上表现为 $\Gamma_{i,j}=1$), 也就是说相速度即为

$$ V_{\varphi}(\omega_{n}) = \frac{\Delta x}{\tau(N)} $$

4 月

04.08

很不容易借到了 Nonlinear Acoustics 教科书, 是 1997 年发行的, 能够在今天再看到着实不易. 很巧的是, 看到了前一位借书的同学是来自密院的郝家兴同学, 并且里面还附带了有关复印备注信息的纸条, 而且借书就在 4 月 5 日, 没有过多久.

既然在网上都没有找到什么这本书的扫描版, 我想我也应该去复制一份扫描存储, 在 DOI 高度发达的今天竟然还有书籍没有电子版, 实在是不能让人忍受的行为.

04.15

成功报名游泳培训班, 这下不用寄人篱下了, 可以安然度过本科最后的时光.

这次组会没有上去讲, 因为上一周因为杂七杂八的事情也没怎么做实验, 这周要把上周的进度补回来.

零碎的想法

1. 怎么使用一个声学探头, 以目前的信号发生器来实现混频($f_{1}\leq f_{2}$)?

2. 基波和谐波的衰减测定方法: 需要用到多距离拟合 $\alpha_{3f}\neq \alpha_{3*f}$

紧锣密鼓地水论文。去除了一下模板里面关于用法部分的水分，发现虽然从昨天到今天也写了一部分，但是字数只有 4890/5078，这显然是令人沮丧的，准备继续细化一下论文内容。

04.18

你麻痹不是才期中吗，我怎么感觉我要死了。

感觉诸事不顺…买的插线板今天终于到货了，但是发现之前买的地面走线盒大小又不合适了，需要现买，到货又需要一段时间，实验无法进行。我在想要不然直接通电先把数据拿到得了

毕业论文。看规范是一共要写到两万字，现在最多也就才写到四分之一。

NI 多通道采集卡。今天和经销商聊了一下并且带到老师那里去简单交流了一下，确认了大致的需求。虽然财务应该还是大师姐来管，但是有关力学传感器和一部分仪器的研判估计还是少不了我的份。

游泳。傻逼游泳馆只接受校园卡换手带，极限转点找宣嘉乐借了卡，但是明天晚上有软凝聚态的课也补不了进度。麻。

b 站一位大佬过了地灵殿 LNN，可喜可贺。我的 lnb 依旧搁浅中…

04.22

为了测衰减做出来的效果并不好, 可能的原因:

1. 探头响应太低了;
2. 施加应力太小. 目前是 $9\text{ kPa}$. 继续加金属块.

方案一

在低频域范围, 用小 $\delta f$, G150 细测一遍钢珠衰减.

方案二

用 W800 测.

除此以外, 还有改变振幅测可逆性的实验, 注意信号放大器最大能放大到 25 V, 这一点还需要注意.

1. 测量不同初始振幅, 各 Gain 对应的振幅值.

$$ \text{dB} = 10\lg{\frac{V}{V_{0}}} $$

$$ 2\text{dB} = 10\lg{\frac{V_{2}}{V_{0}}} - 10\lg{\frac{V_{1}}{V_{0}}} = 10\lg{\frac{V_{2}}{V_{1}}}\Rightarrow \lg\frac{V_{2}}{V_{1}} = 0.2\Rightarrow V_{2} = 10^{0.2}V_{1} $$

所以 Gain dB 对应于

$$ V_{n} = V_{0} 10^{\text{Gain}/10} $$

最大振幅为 25V, 若要对应于 20dB Gain, 则

$$ 25 = V_{0} 10^2 \Rightarrow V_{0} = 0.25\text{V} $$

若初始振幅为 $10\text{Vpp}$, 则对应的 Gain 为

$$ 25 = 10^{\text{Gain}/ 10}\Rightarrow \text{Gain}\approx 13.97 $$

实测下来 Gain 数还会比这个更少, 所以还是需要具体打表确认.

2. 源振幅与颗粒结构之间的关系.

这次要彻底模仿 Jia 的振幅时域 protocol, 从而更好地复现结果.

3. 相速度测量.

想办法废物利用一下 G150 的数据. 浪费好好的周末时间, 还是多少出点数据吧.

4. 扩散模型.

2024.04.23

通过 $\Delta V = 1\text{ Vpp}$, $\Delta\text{Gain}=2\text{ dB}$ 打表测了振幅. 可能需要写一个小程序输出一下从小到大顺序的源振幅的 $(V_{\text{AFG}},\text{Gain})$ 组合.

我想模仿的图像

“严谨设计” 的实验做出来的效果反而不如之前的好, 有点幽默了.

$$ v = f\lambda\Rightarrow \lambda_{\text{eff}} = \frac{V_{\text{tof}}}{f_{c}} = \frac{300\text{ m/s}}{30\text{kHz}} = 10\text{ mm} \approx 3d $$

共振峰可能不止 30 khz 一个. 具体实验的时候可以大概扫一下. 31.3

成功地把 kickstart 部署到电脑上了. 可以自由地控制信号发生器.

04.25

重新做一下 irreversible 的实验。

再用类似的protocol做一遍, 效果还不好就算了.

$f_g = 47\times 10^3\text{ Hz}, f_s = 31.2\text{MHz}$

实验流程：

1. 制备玻璃珠的随机密堆积。
2. 粗略扫频，寻找共振峰。
3. kickstart 连接示波器和发生器。 （1）调整发生器为 10-Cycle 正弦脉冲波，幅值初始为 2.0 Vpp（间隔为 0.5 或者 0.25）。 （2）调整示波器设置：双通道，采样长度200k，平均采集128，探头衰减x1.

颗粒介质对 protocol, 即制备历史具有记忆性. 1-2-1-3-1-4 和1-2-3-4 是完全不同的激励方式. 两者对于相空间探索的方式完全不同, 造成谐波分析对于前者而言很难成功.

后续做一个连续激励的, 即复现简陋装置的实验.

f = 47k

似乎声学并不一定是软化.

周志刚介绍平均场模型的附录内容有问题。多了一个莫名其妙的 $v_0$。在我自己的论文里重写了一遍.

声学黏度和动力学黏度(爱因斯坦关系)是不是等效的呢？

可能不能使用共振峰信号。这个信号对于颗粒介质的诱导效果最大。

再做一组连续的 protocol, 要是还是没效果那我就真没辙了, 用简陋实验的数据写论文算了.

9 kpa, 32 khz, 512 次平均采集, 5-循环正弦波列, 触发间隔 20ms.

循环数都是越大越好吗?

先记 10 个 1VPP 的点, 再开始 1-10, delta = 0.25 的点, 最后再做一点弛豫/aging 的点.

不知道这个波包分离,观察到的是回波(echo)还是压缩波与剪切波的分离. 作为压缩波的探头, 是否能够探测到剪切波的存在? 这个现象似乎还是因为减少了循环数才观察到的.

5 月

毕业论文完成得差不多了, 就剩一点结尾. 喜提新课题, 完全没有接触过, 论文成堆地要去读. 但是我的作业和论文还有不少没做, 还有傻逼游泳测试.