随着火力发电效率的提升，火电机组关键部件之一的电站“四大厚壁管道”系统(锅炉给水管道、主蒸汽管道、再热蒸汽热段管道和再热蒸汽冷段管道)的壁厚也在不断增加. 其在长期的服役过程中，由于高温运行、材料蠕变蜕化及应力松弛，极易造成管壁局部减薄，最终导致管体轴向爆裂泄漏等严重事故[1]. 因此厚壁管道的无损检测对于保障中国电力行业快速安全发展具有极为重要的意义. 相对于传统的超声波检测技术,超声导波传播距离远、速度快,在在役管道的无损检测中有良好的应用前景[2-3]. 何存富等[4-5]、李隆涛等[6]研究了管道中周向导波的传播机理和频散特性；Clough等[7]利用EMAT传感器在管道中激励周向导波，研究了其在缺陷检测和大小测量方面的有效性；吴斌等[8]、刘增华等[9]利用斜探头对厚壁管表面缺陷进行了扫描检测和成像处理；何跃斌等[10]针对不同倾斜角度的变径管,分析了周向Lamb波的传播特性,给出了激励信号的最佳周期数；高广健等[11]在二阶微扰近似条件下,采用导波模式展开分析方法研究了圆管结构中周向导波的非线性效应；Han等[12]利用不同入射角度的楔块，激励周向导波对厚壁管道内壁和外壁的缺陷进行了检测；戴翔等[13]将数值模拟与实验相结合，研究了周向导波快速检测厚壁管道的方法，优化了探头激励角度范围.

由上可见，国内外众多学者针对管道，尤其是厚壁管道的检测做出了许多卓有成效的工作. 作为厚壁管道超声导波检测的关键部分，传统的导波传感器由于对声束的宽度没有进行专门的优化，其声束较宽，会导致导波能量发散、检测灵敏度降低[14]. 一些学者对导波声场的优化做了卓有成效的研究，譬如通过对传感器阵列的信号进行幅值加权处理[15]，或在激励时对阵列传感器的每个阵元施加时间延迟，通过聚焦[16]来优化其指向性. 但目前厚壁管专用的窄声束超声导波传感器仍有待研究. 因此本文针对上述问题研制了一种声束窄、能量集中的阵列传感器. 首先选取了合适的检测模态，以及该模态的激励频率与角度；然后通过有限元仿真的方法，对该阵列传感器进行了声压场计算，并对其阵元个数(4～8阵元)和阵元间距(0.5λ～1.0λ)进行了优化. 随后制作了传感器实物并进行了声场指向性测量. 实验结果表明，相较于单阵元传感器，该传感器有效地集中了在厚壁管道周向上所辐射出的声场. 该研究对于后续该类型传感器的进一步优化研制具有指导意义.

2015年1月—2016年12月我科收治椎弓根螺钉固定患者40例，其中男27例，女13例；年龄为53.23±16.72（35～71）岁；BMI指数为18.92±2.56kg/m2；颈椎疾患8例，腰椎疾患32例；单节段固定31例，多节段固定9例；随机分为对照组（常规置钉，20例）和实验组（椎弓根螺钉定位导向置钉一体器，20例），两组患者临床一般资料具有统计学可比性（P＞0.05）。

1 厚壁管道超声导波检测模态的选取

选用外径为299 mm、壁厚为16 mm的厚壁钢管作为研究对象，建立空心圆柱体弹性波问题的数学模型，以Disperse为计算平台，通过数值计算的方法得到厚壁管道中周向导波的频散曲线，如图1所示.

早在人类祖先造出小船在海上航行之前，甚至早在地球古大陆分离之前，死树就已经漂洋过海，开始了它们从森林到海洋的新旅程。而如今，当一棵树倒下进入河流或溪流中，却很少有人会去关注或了解它们。

图1 尺寸为299 mm×16 mm的钢管的频散曲线及波结构 Fig.1 Disperse curves and wave structure of a 299 mm×16 mm steel pipe

图1(a)(b)中的虚曲线为周向导波中的周向类S1模态，图1(c)为该模态在260 kHz时的波结构图. 当频率为260 kHz时，该模态具有以下几个优点，是一种较为理想的厚壁管道检测模态.

1) 从图1(a)可以看出，在众多的导波模态中群速度最快，为4 662 m/s，较易分辨；该频率下该模态的群速度达到了极大值，频散效应较小.

第二，坚持以人为本，实施暖人心工程。始终把以人为本作为思想政治工作的核心，注重人文关怀和心理疏导。领导干部定期深入基层调研，听取员工对企业发展的建议，了解员工的实际需求。建立困难家庭档案和联系点，以春节慰问、夏送清凉、金秋助学、冬送温暖为载体，实行一对一困难帮扶。通过入户家访、节日慰问、大病救助、子女助学等形式，帮助困难职工解决生活困难，对于不能解决的问题，视情况及时向上级反映或做好耐心细致的解释工作，使职工群众感受到组织的关怀和温暖。

3) 从图1(c)可以看出，该模态以周向位移为主，轴向位移为零，因此在传播过程中以周向能量为主，利于周向传播；该模态的位移在整个管壁径向上都有分布，在内外壁位移最大，可检测整个管壁厚度方向上的缺陷，尤其是最容易出现缺陷的内外壁.

2 窄声束阵列传感器的有限元仿真优化

Wooh等[17-18]和He等[19]研究了阵列传感器的阵元个数(N)、阵元间距(p)、阵元宽度(l)与纵波波长(λ)等参数对声场分布的影响，并且指出改变这些参数可以改变阵列传感器在弹性固体中的声场. 在此基础上，本节将通过有限元仿真的方法对窄声束阵列传感器的声场进行优化.

综上，本研究结果显示CIPN与糖尿病之间关系密切，糖尿病可以增加CIPN的发生风险，提示在临床上为患有糖尿病的肿瘤患者选择化疗方案时应仔细权衡利弊，调整最佳的治疗方案，并积极预防CIPN，以提高患者的生活质量。鉴于目前缺乏针对CIPN发生风险与糖尿病相关性的大规模的临床研究，今后仍需要进行更多的前瞻性、大样本临床研究，以进一步论证CIPN发生风险与糖尿病的相关性。

在第1节中确定了所需的检测模态为周向类S1模态，并且从频散曲线中得到其检测频率为260 kHz，因此波长(λ)已经确定；此外，在后文中制作该传感器所采用的压电阵元的宽度为8.5 mm，因此阵元宽度(l)也已经确定. 故本节将讨论不同阵元个数(N)与阵元间距(p)下该传感器的声场及优化.

2) 从图1(b)可以看出，该模态的激励角度为28.4°，与别的模态相差较大. 因此若以28.4°的入射角度在该管道中激励超声导波，激励出的导波模态会相对单一.

图3(a)～(e)为上述仿真计算的指向性结果. 可以看出，阵元个数及阵元间距的变化确实引起了传感器声场的变化. 将图3(a)～(e)中主瓣宽度最窄的指向性曲线绘制在同一幅图中，即可得到N=4～8时的最优指向性曲线，如图3(f)所示. 可以看出，当N=5，p=0.7λ时，窄声束阵列传感器的主瓣最窄，旁瓣和栅瓣也得到了较好的抑制. 因此该参数下窄声束阵列传感器的声场指向性最好.

该模型中厚壁管道的外径为299 mm，壁厚16 mm，材质为钢，密度ρ=7 850 kg/m3, 纵波波速cL=6 144.0 m/s, 横波波速cS=3 094.8 m/s；楔块的材质为有机玻璃，密度ρ=1 290 kg/m3, 纵波波速cL=2 702.2 m/s, 横波波速cS=1 330.5 m/s. 网格类型为自由剖分四面体网格，网格尺寸为其波导中纵波波长的1/10.

窄声束阵列传感器仿真采用的阵元入射角度为28.4°，该角度为第1节中所讨论的260 kHz时可以激励出周向类S1模态的入射角度. 在图2(b)中箭头所指的应力加载面上施加法向的简谐应力来模拟压电阵元产生的应力，频率为260 kHz. 由于后文制作传感器实物时，用到的压电片的长×宽为28.5 mm×8.5 mm，因此此处每个应力加载面的长×宽为14.25 mm×8.5 mm.

苏霍姆林斯基说过：“不了解儿童的心灵，就谈不上教学素质。”可见，要想促使学生良好道德品质的形成，必须首先对学生身心发展特点有精准把握。小学生天性好玩，对各种游戏有浓厚的兴趣，表现欲望强烈。学生的这一特点必须引起教师高度关注，结合所学内容组织趣味游戏，让学生在游戏中获得思想的成长，陶冶他们的道德情操。游戏是内化学生道德情感，培养学生良好个性品质的重要途径。借助趣味游戏开展教学活动已然被认可，但如何进行具体操作又是另外一个问题。

传感器的辐射声场通常采用指向性衡量. 如图2(a)所示，本仿真在传感器的前端，厚壁管道的外表面上设置了一条声压级取值线，用以得到该直线上的声压级. 将得到的声压级归一化后，即可得到该传感器的声场指向性. 仿真中采用的阵元个数为N=4，5，6，7，8，共5个；在每个阵元个数下分别仿真计算阵元间距为p=0.5λ，0.6λ，0.7λ，0.8λ，0.9λ，λ时的声场. 因此该仿真计算共进行了30组.

图2 窄声束阵列传感器声场三维有限元仿真模型 Fig.2 Model of sound field simulation of narrow beam array transducer

由于该传感器需要在管道周向上激励超声导波，在管道轴向上保持较小的传感器尺寸，利于提升其在管道周向上的声束指向性，因此该传感器采用沿管道周向上一字排开的阵列布置形式. 以COMSOL为有限元仿真计算工具建立三维模型. 由于本模型要计算的问题是对称的，因此对模型沿中心进行了切割并采用了对称边界，如图2所示. 该操作在不影响最终结果的前提下减小了50%的计算量，提高了计算效率. 图2以N=5，p=0.8λ时建立的三维模型为例. 需要说明的是，阵元间距指的是2个阵元的中心沿厚壁管外壁周向上的距离.

图3 窄声束阵列传感器指向性的优化结果 Fig.3 Directivity optimization results of NBAT

对比起见，也对相同设置下单阵元传感器的声场进行了仿真计算. 图3(f)中的灰色虚线为其指向性，图4为声压级场的对比图. 可以看出，优化后的窄声束阵列传感器在厚壁管道周向上的声束集中，指向性优于传统单阵元传感器.

4 窄声束阵列传感器制作及指向性测量

4.1 传感器制作

根据在COMSOL中建立的模型，将窄声束阵列传感器的楔块按照1∶1的比例进行实际加工. 楔块的材质为与仿真中参数相同的有机玻璃；压电阵元的长×宽×高为28.5 mm×8.5 mm×8.0 mm，振动方式为厚度伸缩型，在压电片置入槽内以纵波耦合剂耦合于楔块之上. 经阻抗分析仪(安捷伦4294A)实测，其谐振频率为214 kHz，第二谐振频率为266 kHz,如图5所示.

图4 窄声束阵列传感器及单阵元传感器的声场对比 Fig.4 Sound fields of NBAT and single element transducer

图5 窄声束阵列传感器实物及压电阵元阻抗 Fig.5 NBAT photo andimpedance of PZT

4.2 窄声束阵列传感器的指向性测量

图6 指向性测量系统示意图及实物图 Fig.6 Directivity measuring system

指向性测量实验在一根长1.5 m、直径299 mm、壁厚16 mm的厚壁管道上进行，如图6所示. 为了降低边界影响，窄声束阵列传感器放置于管道外壁中央. 由于周向类S1模态以周向位移为主，因此以宽频横波直探头(V151，泛美)作为接收探头，测量其面内位移. 激励信号采用汉宁窗调制的五周期正弦波信号，中心频率260 kHz. 进行指向性测量时，函数信号发生器(AFG3021B, 泰克)产生激励信号，经电压放大器(2340, TEGAM)放大后施加于窄声束传感器之上；横波直探头中心到窄声束阵列传感器前端中心的距离为150 mm，在一个半圆(θ=-90°～90°)上每隔10°测量一组声波信号，共19组；测量得到的信号经由数字示波器(DPO 4054, 泰克)显示采集后输入计算机.

图7 θ=-70°～0°时的时域信号 Fig.7 Time-domain signal of θ=-70°-0°

图7为指向性测量实验的典型时域信号(θ=-70°～0°). 其中图7(a)为当横波直探头处于窄声束阵列传感器正前方(θ= 0°)时测得的时域信号. 第1个波包，也就是虚线框里的模态即为周向类S1模态，经计算该波包的群速度为4 419.5 m/s，与Disperse计算得到的周向类S1模态理论群速度的相对误差为5.2%. 通过比较图7(a)～(i)可以看出，时域信号中周向类S1模态幅值的总体趋势趋于减小. 将所有时域信号中第1个波包的幅值提取出来，归一化后绘制在同一幅图中，即可得到窄声束阵列传感器在厚壁管道周向上激励出周向类S1模态的指向性，如图8所示. 其中红色方圈为实验数据点，红实线为其指向性拟合曲线；为对比本实验也测量了相同设置与参数下单阵元传感器的声场指向性，实验数据为如图8中的绿色圆圈，指向性拟合曲线如绿虚线所示. 可以看出，当θ=-48°～43°时，传统的单阵元传感器的声束较宽，并且在θ=0°两侧出现了2个指向性的极大值. 窄声束阵列传感器的声束较窄，并且θ= 0°附近的声束能量远大于其他范围的声束能量；当θ=-48°～-63°或θ=43°～60°时，窄声束阵列传感器与传统单阵元传感器的声场指向性特性较为相近；当θ<-63°或θ>-60°时，窄声束阵列传感器出现了旁瓣，这是由于窄声束阵列传感器在管道周向上的长度较长，多阵元在靠近管道轴向上的位置产生能量泄露造成的. 该范围内窄声束传感器激励出的声波能量较小，所以旁瓣对该传感器在实际检测中的影响有限.

图8 指向性测量实验结果 Fig.8 Directivity measuring results

若以图8中指向性的最大值1为基准值，分别在衰减-2、-4、-6 dB时对图8中2条指向性拟合曲线进行取值，则可得到表1所示的结果. 从表2可以看出，当衰减-2 dB时，窄声束阵列传感器声束分布的角度(θ1)为传统单阵元传感器(θ2)的48.5%；衰减-4 dB时，该值为79.7%；衰减-6 dB时，该值为86.4%. 这表明相较于传统的单阵元传感器，窄声束阵列传感器有效地集中了在厚壁管道周向上辐射出的声场；并且越靠近θ=0°时声束能量越集中，指向性越好.

第二，旅游危机事件网络舆情传播的各个阶段对旅游危机事件的关注点不尽相同。旅游危机事件本身是网络舆情的信息来源。舆情传播的潜伏期通常以事件始末为开端；在舆情传播的爆发期，公众对旅游危机事件的关注扩大到事件的影响和相关管理部门的应对措施；而在舆情传播的成熟期，旅游危机事件的影响和相关管理部门的应对措施成为公众关注的焦点；在舆情传播的平息期，旅游危机事件的影响和结果成为公众关注的主要内容。

表1 不同衰减下传感器的指向性特性

Table 1 Transducer directivity characteristics with different attenuations

衰减值/dB窄声束阵列传感器传统单阵元传感器起始角度/(°)结束角度/(°)角度范围θ1/(°)起始角度/(°)结束角度/(°)角度范围θ2/(°)θ1·θ-12/%-2-201333-35336848.5-4-333063-41387979.7-6-383876-46428886.4

5 结论

1) 以有限元仿真计算了该传感器4～8阵元、0.5λ～λ阵元间距的声场指向性，仿真结果表明：当阵元个数为5，阵元间距为0.7λ时，窄声束阵列传感器的声场指向性最好.

GAO G J, DENG M X, LI M L, et al. Modal expansion analysis of nonlinear circumferential guided wave propagation in a circular tube[J]. Acta Physica Sinica, 2015, 64(18): 206-215. (in Chinese)

3) 窄声束阵列传感器对于厚壁管道的超声导波检测具有潜在实用价值. 该传感器优化研制方法可对该类型传感器后续的设计与优化提供指导意义.

难治性癫痫患者应尽早进行抗癫痫治疗，丙戊酸钠、奥卡西平单药治疗难治性癫痫取得满意疗效［15］。丙戊酸钠是一种常用的丙戊酸制剂，在20世纪60年代就开始运用于临床。有研究表明丙戊酸钠可以与GABA的合成酶和降解酶产生反应，参与三羧酸循环与钠离子通道，可以增强GABA的活性，可以增加GABA的水平，从而发挥抗癫痫作用［16］。奥氮平可在一定程度上改善了患者的工作记忆功能，促进前额叶皮质多巴胺的释放，从而发挥治疗效果［17］。本研究显示观察组的总有效率98.9%（94/95）高于对照组88.4%（84/95）（P＜0.05）。

参考文献：

干旱是在全球气候变化和人类活动双重驱动下，水资源—经济社会—生态与环境复杂系统极值缺水过程在不同时空尺度上的综合体现。按受旱机制分类，干旱分为：气象干旱、水文干旱、农业干旱和经济社会干旱。按受旱成因分类，干旱分为：资源型干旱，即干旱区降水量稀少，土壤缺水，不足以满足人的生存和经济发展的需求；工程型干旱，即缺少水利工程的开发调度；政策型干旱，即非理性的水资源分配；结构型干旱，即水质水量达不到用水要求；布局型干旱，即追求经济社会快速发展情势下的水土资源不匹配；管理型干旱，即水资源得不到高效利用和有效调配，缺少全局可靠规划与应对策略；复合型干旱，即多种因素综合作用下的干旱。

[1] 侯锡瑞. 电站锅炉压力容器压力管道安全技术[M]. 北京: 中国电力出版社, 2005: 1-4.

[13] 戴翔, 罗更生, 张勇芳, 等. 厚壁管道周向超声导波缺陷检测研究[J]. 实验力学, 2016， 31(5): 635-644.

ZHOU Z G, FENG H W. Progress in research of ultrasonic guided wave testing technique[J]. Nondestructive Testing, 2006, 28(2): 57-63. (in Chinese)

[3] 何存富, 郑明方, 吕炎, 等. 超声导波检测技术的发展、应用与挑战[J]. 仪器仪表学报, 2016, 37(8): 1713-1735.

上述译文中，“定罪量刑”标准翻译为“crime”。而crime仅有犯罪含义，与定罪、量刑这两个原文单词相差甚远，不符合法律英语的简短性与精确性，使译文所表达之意影响了原文辞义表达的连贯性。正确译文应为The criteria of conviction and sentencing

[10] 何跃斌, 陈果, 卢超, 等. 变径管周向导波Lamb波传播特性[J]. 无损检测, 2017, 39(7): 1-7.

[4] 何存富, 李隆涛, 吴斌. 周向超声导波在薄壁管道中的传播研究 [J]. 实验力学, 2002, 17(4) : 419-424.

HE C F, LI L T, WU B. Research on guided circumferential waves in a thin-walled pipe[J]. Journal of Experimental Mechanics, 2002, 17(4): 419-424. (in Chinese)

[5] 何存富, 李隆涛, 吴斌. 空心圆柱体中的周向超声导波[J]. 机械工程学报, 2004, 40(8): 7-12.

HE C F, LI L T, WU B. Guided circumferential waves in hollow cylinders [J]. Chinese Journal of Mechanical Engineering, 2004, 40(8): 7-12. (in Chinese)

[6] 李隆涛, 何存富, 吴斌. 周向超声导波对管道纵向缺陷检测的研究[J]. 声学学报, 2005, 30(4): 343-348.

LI L T, HE C F, WU B. Detection of a longitudinal defect in a pipe based on guided circumferential waves techniques[J]. Acta Acustica, 2005, 30(4): 343-348. (in Chinese)

20 世纪初的民国时期可谓气象万千，政局不定，社会动荡，战争密集，都对图书造成了极大的损毁。而东西方思想交错，文化革新使得图书的捐赠成为图书的保存和整理、大众借阅和共享的一种重要途径。另外，民国时期所捐赠的图书多为精品，孤本、善本不乏其中，更为金钱所无法估量。正是有赖于此，中国近代的文化事业、大众的思想文化才会庚续不断，方兴未艾。民国时期图书馆的图书捐赠活动，是对文化的普及，思想的传播，经典的传承的一种良性循环。

[7] CLOUGH M, FLEMING M, DIXON S. Circumferential guided wave EMAT system for pipeline screening using shear horizontal ultrasound[J]. NDT & E International, 2017, 86: 20-27.

[8] 吴斌, 颉小东, 刘增华, 等. 厚壁管道周向导波检测技术实验研究[J]. 实验力学, 2011, 26(5): 625-631.

WU B, XIE X D, LIU Z H, et al. Experimental investigation on thick wall pipes based on circumferential guided waves inspection technology[J]. Journal of Experimental Mechanics, 2011, 26(5): 625-631. (in Chinese)

[9] 刘增华, 颉小东, 吴斌, 等. 基于连续小波变换的厚壁管道周向导波扫描成像试验研究[J]. 机械工程学报, 2013, 49(2): 14-19.

LIU Z H, XIE X D, WU B, et al. Experimental research on circumferential guided wave scanning imaging of thick wall pipes based on continuous wavelet transform[J]. Journal of Mechanical Engineering, 2013, 49(2): 14-19. (in Chinese)

HE C F, ZHENG M F, LÜ Y, et al. Development, applications and challenges in ultrasonic guided waves testing technology[J]. Chinese Journal of Scientific Instrument, 2016, 37(8): 1713-1735. (in Chinese)

HE Y B, CHEN G, LU C, et al. Propagation characteristics of circumferential guided wave soft aperpipe[J]. Nondestructive Testing, 2017, 39(7): 1-7. (in Chinese)

[11] 高广健, 邓明晰, 李明亮, 等. 圆管结构中周向导波非线性效应的模式展开分析[J]. 物理学报, 2015, 64(18): 206-215.

2) 以仿真结果为基础制作了传感器实物. 建立了一套导波传感器声场的测量系统，对该传感器的声场进行了实验测量. 实验结果表明：相较于传统的单阵元传感器，该传感器有效地集中了在厚壁管道周向上辐射出的声场；并且越靠近θ=0°时声束能量越集中，指向性越好.

[12] HAN Z W, LU C, SONG K, et al. Research on detection experiment for defected pipes hidden behind thick-walled using ultrasonic circumferential guided waves[C]∥Mechatronics and Automation Engineering: Proceedings of the International Conference on Mechatronics and Automation Engineering (ICMAE2016). Hongkong: Hongkong Society of Mechanical Engineers, 2017: 277-283.

[2] 周正干, 冯海伟. 超声导波检测技术的研究进展[J]. 无损检测, 2006, 28(2): 57-63.

DAI X, LUO G S, ZHANG Y F, et al. On the defect detection of thick-walled pipe by using circumferential ultrasonic guided wave[J]. Journal of Experimental Mechanics, 2016， 31(5): 635-644. (in Chinese)

[14] 郭程. 基于横波直探头的管道导波激励特性及损伤检测研究[D]. 北京: 北京工业大学, 2016.

GUO C. Research on excitation characteristics of guided waves and damage detection based on bulk shear wave sensor in pipeline[D]. Beijing: Beijing University of Technology, 2016. (in Chinese)

式中：CRn表示产业内规模最大的前n个企业的行业集中度；Xi表示产业内第i位企业的相关经营数据 （产值、销售额等）；n表示前几位大企业的数量；N表示产业内的企业总数。

在连续使用多次单倍体技术，持续对品种进行选育的过程中，用两个DH系通过杂交之后所形成的基础群体进行第二次的选系，可以显著提升单倍体的加倍频率。

[15] MA Q Z, JIAO J P, HU P, et al. Beam forming of lamb waves for nondestructive testing of plates[J]. Applied Mechanics & Materials, 2014, 490/491: 1512-1516.

《世说新语》中有个小故事：“殷中军虽思虑通长，然于才性偏精。忽言及四本，便若汤池铁城，无可攻之势。”意思是殷浩虽然思维能力一般，但对于才性之学特别精通。如果谈到《四本论》，他就像镇守着汤池铁城，有坚不可摧之势。教学中是不是也存在同样的问题呢：作为教师的我，特长在哪里？教学中如何将其发挥出来？

[16] 付君强. 二维阵列导波聚焦方法及无损检测系统研究[D]. 北京: 北京理工大学, 2016.

FU J Q. Two-dimensional array guided wave focusing method and nondestructive testing system research[D]. Beijing: Beijing Institute of Technology, 2016. (in Chinese)

[17] WOOH S C, SHI Y. Optimum beam steering of linear phased arrays[J]. Wave Motion, 1999, 29(3): 245-265.

[18] WOOH S C, SHI Y. A simulation study of the beam steering characteristics for linear phased arrays[J]. Journal of Nondestructive Evaluation, 1999, 18(2): 39-57.

[19] HE C, ZHANG X, LU Y, et al. Sound fields radiated from a linear phased array based on analytical solutions of transient elastic waves[J]. Wave Motion, 2017,73: 86-95.