激光声表面波：早期龋齿与牙隐裂精准诊断的新曙光

激光声表面波：早期龋齿与牙隐裂精准诊断的新曙光一、绪论1.1研究背景口腔健康是全身健康的重要组成部分，世界卫生组织将口腔疾病列为全球重点防治的非传染性疾病之一。在各类口腔疾病中，早期龋齿和牙隐裂较为常见，严重影响患者的生活质量。龋齿，俗称蛀牙，是一种细菌性口腔疾病，主要特征为牙齿硬组织的慢性进行性破坏。初期通常表现为牙齿表面的白垩色斑，随着病情发展，会逐渐形成龋洞。相关数据显示，全球范围内龋齿的发病率一直居高不下。《第四次全国口腔健康流行病学调查报告》指出，我国5岁儿童乳牙龋患率为70.9%，12岁儿童恒牙龋患率为34.5%，35-44岁中年人群龋患率为89.0%，65-74岁老年人群龋患率更是高达98.0%。早期龋齿若未得到及时诊断和治疗，病变会持续发展，细菌进一步侵蚀牙齿深层组织，引发牙髓炎、根尖周炎等并发症，导致剧烈疼痛，严重时甚至会造成牙齿丧失，影响咀嚼功能和营养摄入，进而对全身健康产生负面影响。牙隐裂则是指牙冠表面出现的非生理性细小裂纹，这些裂纹通常难以用肉眼直接察觉。牙隐裂的发生与多种因素有关，如咬硬物、磨牙及咬牙习惯、温度作用、牙体修复材料、创伤性牙合力、牙齿发育缺陷、牙尖斜度以及牙齿增龄性变化等。由于牙隐裂早期症状不明显，患者可能仅在咬硬物或遇到冷热刺激时感到短暂不适，容易被忽视。但随着时间推移，隐裂纹会逐渐加深，一旦裂缝扩展至牙髓腔，细菌侵入就会引发牙髓炎、牙髓坏死等问题，导致牙齿持续疼痛，影响正常生活。若裂纹进一步发展，还可能致使牙齿劈裂，最终不得不拔除患牙，给患者带来极大的痛苦和经济负担。早期诊断对于早期龋齿和牙隐裂的治疗及预防病情恶化至关重要。在龋齿早期，病变仅局限于牙齿表面的釉质层，此时进行干预治疗，操作简单，治疗效果好，患者痛苦小，费用也相对较低。例如，通过专业的诊断方法发现早期龋齿后，可以采用非创伤性修复治疗技术，使用专门的牙科材料对龋损部位进行填充，就能有效阻止龋齿的进一步发展。而对于牙隐裂，早期发现后采取调磨牙尖、降低咬合等保守治疗措施，能够减轻患牙的负担，防止裂纹加深，避免牙齿出现更严重的病变，最大程度地保留牙齿的功能和结构完整性。然而，传统的口腔疾病诊断方法在早期龋齿和牙隐裂的检测上存在一定的局限性。例如，常规的口腔视诊和探诊对于早期龋齿和细微的牙隐裂难以准确判断；X线检查虽然能发现一些明显的龋齿和牙齿病变，但对于早期龋齿和浅表的牙隐裂敏感度较低，容易出现漏诊；染色法、光纤光源透照法等虽然在一定程度上有助于牙隐裂的诊断，但也存在操作复杂、准确性有限等问题。因此，寻找一种更加准确、高效、便捷的早期诊断方法，对于提高早期龋齿和牙隐裂的治疗效果、改善患者的口腔健康状况具有重要意义。1.2研究意义激光声表面波技术作为一种新兴的无损检测技术，在早期龋齿和牙隐裂的诊断与评估中展现出独特的优势，具有重要的研究意义。在诊断准确性方面，传统诊断方法对于早期龋齿和细微牙隐裂的检测存在局限性，容易导致漏诊和误诊。激光声表面波技术利用激光脉冲照射牙齿表面，激发出声表面波，通过分析声表面波的传播特性，如速度、频率、振幅等变化，能够精确识别牙齿内部的微小结构变化和缺陷。相关研究表明，在早期龋齿检测中，激光声表面波技术能够检测到牙齿釉质表面下深度小于0.1mm的早期脱矿病变，相比传统X线检查，其对早期龋齿的检出率提高了30%以上。对于牙隐裂，该技术可以检测到宽度小于0.05mm的细微裂纹，显著提高了牙隐裂的早期诊断率，为及时治疗提供了有力保障。从患者体验角度来看，激光声表面波技术属于非接触式检测，避免了传统探诊等方法对牙齿和口腔组织的直接接触刺激，大大减少了患者在检查过程中的不适感和痛苦。同时，该技术检测过程快速，一般单个牙齿的检测时间在数秒内即可完成，相较于一些复杂耗时的传统检查方法，如牙髓活力测试等，能有效缩短患者的就诊时间，提高就医效率，尤其对于儿童、老年人以及对口腔检查较为敏感的患者群体具有重要意义。在临床治疗决策方面，激光声表面波技术能够提供详细的牙齿病变信息，帮助医生更准确地评估早期龋齿和牙隐裂的严重程度、病变范围及发展趋势。医生可以根据这些精确的诊断结果制定个性化的治疗方案，避免过度治疗或治疗不足的情况发生。例如，对于早期龋齿，根据激光声表面波技术检测到的病变深度和范围，医生可以选择合适的治疗方法，如非创伤性修复治疗、预防性树脂充填等；对于牙隐裂，根据裂纹的深度和长度，决定是采取调合、全冠修复还是根管治疗等不同的治疗措施，从而提高治疗效果，最大程度地保留牙齿的功能和结构完整性，降低患者后续治疗的复杂性和成本。激光声表面波技术的研究与应用，还能推动口腔医学诊断技术的创新发展，为口腔疾病的早期诊断和精准治疗提供新的思路和方法，具有重要的科学价值和临床应用前景，有望在未来口腔医疗领域发挥重要作用，改善广大患者的口腔健康状况。1.3国内外研究现状1.3.1激光超声无损检测发展概述激光超声无损检测技术起源于20世纪60年代，R.M.White于1963年首次发现脉冲激光辐射固体材料表面时，会在样品表面激发出低于激光频率的声表面波，这一发现为激光超声技术的发展奠定了基础。此后，众多学者围绕这一现象展开了深入研究，推动了激光超声无损检测技术的逐步发展。20世纪70年代中期，激光超声技术作为一种新兴的无损检测技术开始崭露头角。该技术利用Q开关脉冲激光器发出的激光束照射被测物体，激发出超声波，通过检测超声波的特性来获取物体内部的信息。在这一时期，研究主要集中在激光超声的激发机理和基本检测方法上。随着研究的不断深入，激光超声检测技术逐渐形成了电学检测法和光学检测法两大类。电学检测法包括接触式和非接触式两种类型，接触式主要利用压电换能器将超声信号转化为电信号，具有较高的灵敏度和适中的价格，在超声检测中被广泛应用，但检测时需要与被测物体表面接触或距离很近；非接触式检测方法则包含电容换能器以及电磁-声换能器。光学检测法包含非干涉法和干涉法，非干涉法使用光反射技术、光偏转技术以及光衍射技术等，将超声波信号调制到光强信号中，通过光电检测器进行检测；干涉法则借助声波引起的光束频率和相位调制来实现测量，常见的干涉仪有外差干涉仪以及共焦F-P干涉仪等。到了20世纪80年代，激光超声无损检测技术在理论和实验方面都取得了重要进展。研究者们对激光与材料相互作用的机理有了更深入的理解，进一步完善了激光超声的激发和检测理论。同时，随着计算机技术和信号处理技术的飞速发展，激光超声检测系统的性能得到了显著提升，能够更准确地检测和分析超声信号。在应用方面，激光超声无损检测技术开始在航空航天、材料科学等领域得到初步应用，用于检测金属材料、复合材料等的内部缺陷和结构完整性。进入21世纪，激光超声无损检测技术迎来了更快速的发展。新型激光源和探测器的不断涌现，使得激光超声检测的灵敏度、分辨率和检测速度得到了大幅提高。例如，超短脉冲激光器的应用，能够产生更短的激光脉冲，从而提高了超声信号的时间分辨率，使得对微小缺陷的检测成为可能。此外，多模态检测技术的发展，将激光超声与其他无损检测技术（如红外热成像、X射线检测等）相结合，实现了对物体更全面、更准确的检测。在应用领域，激光超声无损检测技术不仅在航空航天、汽车制造、电力能源等传统工业领域得到广泛应用，还在生物医学、文物保护等新兴领域展现出巨大的应用潜力。在生物医学领域，激光超声技术可用于生物组织的无损检测和成像，为疾病的早期诊断提供了新的手段；在文物保护领域，该技术能够对文物的内部结构和损伤情况进行无损检测，为文物的修复和保护提供重要依据。近年来，随着人工智能、大数据等新兴技术的不断发展，激光超声无损检测技术也开始与这些技术深度融合。通过建立超声信号与缺陷特征之间的深度学习模型，能够实现对缺陷的自动识别和分类，提高检测效率和准确性。同时，利用大数据分析技术，可以对大量的检测数据进行处理和分析，挖掘数据背后的潜在信息，为设备的维护和管理提供决策支持。1.3.2早期龋齿和牙隐裂诊断的研究进展早期龋齿的传统诊断方法主要包括视诊、探诊和X线诊断。视诊是通过肉眼观察牙齿表面的颜色、形态等变化来判断是否存在龋齿，这种方法简单易行，但对于早期龋齿，尤其是牙齿表面仅有白垩色斑而无明显龋洞的情况，很难准确判断。探诊则是利用牙科探针探测牙齿表面的硬度和粗糙程度，以发现潜在的龋损，但对于早期龋齿，由于病变部位较浅，探针可能无法准确感知，容易出现漏诊。X线诊断是目前临床上常用的龋齿诊断方法之一，它能够通过X线影像观察牙齿内部的结构，发现牙齿硬组织的脱矿和龋洞。然而，X线检查对于早期龋齿的敏感度较低，尤其是对于牙齿邻面和颊舌面的早期龋损，容易漏诊。研究表明，X线检查对早期龋齿的检出率约为50%-70%。牙隐裂的传统诊断方法主要有视诊、染色法、咬诊和X线检查。视诊时，医生需要仔细观察牙齿表面是否有裂纹，但由于牙隐裂的裂纹通常非常细微，肉眼很难直接发现，尤其是在牙齿的咬合面和邻面，漏诊率较高。染色法是将特定染料（如碘酊、亚甲蓝等）涂抹在牙齿表面，使裂纹染色后更容易观察。然而，这种方法对于较浅的裂纹或被牙菌斑、牙结石覆盖的裂纹，染色效果不佳，影响诊断准确性。咬诊是让患者咬合特定物品（如棉签、牙线等），通过观察患者是否出现疼痛或不适来判断隐裂部位及程度，但该方法主观性较强，不同患者的感受差异较大，且对于早期牙隐裂，疼痛症状可能不明显，容易误诊。X线检查对于牙隐裂的诊断也存在一定局限性，由于牙隐裂的裂纹方向和深度不同，X线影像上可能无法清晰显示裂纹，导致漏诊。研究显示，X线检查对牙隐裂的诊断准确率约为40%-60%。为了克服传统诊断方法的局限性，近年来，新兴诊断技术不断涌现。在早期龋齿诊断方面，激光荧光诊断技术得到了广泛研究和应用。该技术利用激光激发牙齿表面的荧光物质，通过分析荧光强度和波长的变化来判断牙齿是否存在龋损。研究表明，激光荧光诊断技术对早期龋齿的检出率可达到80%以上，显著高于传统诊断方法。此外，光学相干断层扫描（OCT）技术也逐渐应用于早期龋齿的诊断。OCT技术能够对牙齿进行高分辨率的断层成像，清晰显示牙齿内部的结构和病变情况，对于早期龋齿的诊断具有较高的准确性。在牙隐裂诊断方面，口腔显微镜检查法能够放大牙齿表面的细节，使医生更容易观察到细微的裂纹，提高了牙隐裂的诊断准确率。同时，锥形束CT（CBCT）技术也为牙隐裂的诊断提供了更准确的影像学依据。CBCT能够进行三维成像，全面展示牙齿的结构和裂纹的走向、深度等信息，有效弥补了传统X线检查的不足。1.3.3牙齿的超声检测研究进展超声检测技术在牙齿结构和病变检测中有着广泛的应用。在牙齿结构检测方面，超声可以用于测量牙齿的厚度、密度等参数，帮助了解牙齿的生理状态。研究表明，通过超声测量可以准确获取牙齿釉质和牙本质的厚度，为口腔医学研究和临床诊断提供重要数据。在龋齿检测方面，超声检测能够通过分析超声信号在牙齿中的传播特性，判断牙齿内部是否存在脱矿和龋洞。一些研究利用超声的反射和散射特性，成功检测出了牙齿内部的早期龋损，为早期龋齿的诊断提供了新的方法。在牙隐裂检测方面，超声检测可以通过检测裂纹对超声传播的影响，来判断牙隐裂的存在和程度。有研究通过实验发现，当超声遇到牙隐裂时，会发生反射、散射和衰减等现象，通过分析这些变化可以有效检测牙隐裂。然而，现有牙齿超声检测研究在检测精度和适用范围等方面仍存在不足。在检测精度方面，虽然超声检测能够发现一些牙齿病变，但对于早期龋齿和细微牙隐裂的检测，其精度仍有待提高。早期龋齿的病变范围较小，细微牙隐裂的裂纹宽度极窄，现有的超声检测技术难以准确检测到这些微小变化。在适用范围方面，超声检测受到牙齿位置、形态以及周围组织的影响较大。例如，对于位于口腔深部或被其他牙齿遮挡的牙齿，超声检测的准确性会受到影响。此外，牙齿周围的软组织和骨骼也会对超声信号产生干扰，增加了检测的难度。同时，目前的超声检测设备在操作的便捷性和检测速度方面也有待进一步改进，以满足临床快速、准确诊断的需求。1.4本文主要工作本文聚焦激光声表面波在早期龋齿和牙隐裂诊断与评估中的应用，开展了一系列深入研究。首先，全面梳理了激光超声无损检测技术的发展脉络，详细阐述其在早期龋齿和牙隐裂诊断领域的研究现状，同时深入分析了现有牙齿超声检测研究在精度、适用范围等方面的不足。在理论研究方面，深入探究激光与牙齿组织相互作用的物理机制，明确激光声表面波在牙齿中的激发原理与传播特性。通过建立理论模型，详细分析声表面波的速度、频率、振幅等参数在牙齿病变情况下的变化规律，为后续的实验研究提供坚实的理论基础。在实验研究环节，精心搭建激光声表面波检测实验平台。该平台配备高能量脉冲激光器、高灵敏度激光干涉仪以及高精度位移控制系统等先进设备，确保实验数据的准确性和可靠性。针对早期龋齿，利用人工制备的龋齿样本以及临床收集的早期龋齿病例，开展大量实验检测。通过分析激光声表面波在正常牙齿和早期龋齿牙齿表面传播时的信号差异，建立早期龋齿的特征信号数据库，探索基于激光声表面波技术的早期龋齿定量诊断方法。对于牙隐裂，采用在健康牙齿上制造模拟裂纹以及获取临床牙隐裂患牙的方式，进行检测实验。深入研究激光声表面波在遇到牙隐裂时的反射、散射和衰减等现象，分析裂纹深度、长度和方向对声表面波传播特性的影响，建立牙隐裂的检测和评估指标体系。此外，将激光声表面波技术与其他无损检测技术（如激光荧光、光学相干断层扫描等）进行对比研究，通过实验验证其在早期龋齿和牙隐裂诊断中的优势和不足。同时，结合人工智能算法，对激光声表面波检测得到的大量数据进行分析处理，建立基于深度学习的早期龋齿和牙隐裂诊断模型，提高诊断的准确性和效率。最后，总结研究成果，分析激光声表面波技术在实际临床应用中可能面临的问题和挑战，并对未来的研究方向进行展望，为推动激光声表面波技术在口腔疾病诊断领域的实际应用提供参考依据。二、激光声表面波的激发检测机理及其频散特征2.1激光在固体中激发超声的热弹机制2.1.1激光的吸收及其辐照效应当激光作用于固体材料表面时，其与固体材料的相互作用是一个复杂的物理过程。从微观层面来看，激光是由大量的光子组成，这些光子具有特定的能量和动量。当激光光子入射到固体材料表面时，会与材料中的原子、分子等微观粒子发生相互作用。对于金属材料而言，其内部存在大量的自由电子，激光光子与自由电子相互作用，将能量传递给自由电子，使自由电子获得较高的动能。自由电子在材料内部运动过程中，会与晶格原子发生碰撞，将动能传递给晶格原子，从而使晶格原子的振动加剧，宏观上表现为材料对激光能量的吸收。对于非金属材料，其原子通过共价键、离子键等化学键结合在一起，激光光子与材料中的电子云相互作用，激发电子跃迁到更高的能级，在这个过程中，光子的能量被材料吸收。根据能量守恒定律，入射到材料表面的激光能量E_0一部分被反射，记为E_R，一部分被材料吸收，记为E_A，还有极少部分可能会透过材料（对于不透明材料，透过部分可忽略不计），即E_0=E_R+E_A。反射系数R=\frac{E_R}{E_0}，吸收系数\alpha=\frac{E_A}{E_0}，且R+\alpha=1。材料的反射系数和吸收系数与材料的性质、表面状态以及激光的波长、入射角等因素密切相关。例如，金属材料对激光的反射率通常较高，在某些情况下，金属对特定波长激光的反射率可达90%以上，这意味着其对激光的吸收系数相对较低。而一些非金属材料，如陶瓷、玻璃等，对激光的吸收特性则与它们的化学成分、晶体结构等因素有关。材料吸收激光能量后，会产生一系列的辐照效应。其中，最为直接的效应是温度升高。由于激光能量被材料吸收，转化为热能，使得材料内部的微观粒子热运动加剧，温度迅速上升。根据热传导理论，材料内部的温度分布会随着时间和空间发生变化。在激光辐照的初始阶段，材料表面温度急剧升高，形成一个高温区域，随着时间的推移，热量会逐渐向材料内部扩散。在这个过程中，材料内部会产生温度梯度，导致热应力的产生。热应力是由于材料内部不同部位温度变化不一致，热胀冷缩程度不同而产生的应力。当热应力超过材料的屈服强度时，材料会发生塑性变形；若热应力进一步增大，超过材料的断裂强度，材料则可能会出现裂纹甚至破裂。此外，材料吸收激光能量后的膨胀也是一个重要的辐照效应。随着温度的升高，材料原子间的间距增大，宏观上表现为材料的膨胀。对于各向同性材料，其在各个方向上的膨胀程度相同；而对于各向异性材料，如某些晶体材料，其在不同晶向的膨胀系数不同，会导致材料在不同方向上的膨胀程度存在差异。这种膨胀效应在激光激发超声的过程中起着关键作用，它是热弹机制下产生超声的重要原因之一。在激光脉冲辐照材料的极短时间内，材料表面因吸收激光能量而迅速膨胀，这种快速的膨胀会在材料内部产生应力波，进而激发超声。2.1.2热弹机制下激光辐照材料的温度场分布在热弹机制下，研究激光辐照材料的温度场分布对于理解超声的激发过程至关重要。根据热传导理论，对于各向同性的均匀固体材料，在不考虑内热源（除激光热源外）和材料内部热辐射的情况下，其热传导方程可表示为：\frac{\partialT}{\partialt}=a\left(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}}\right)其中，T(x,y,z,t)表示材料在位置(x,y,z)处、时刻t的温度；a=\frac{k}{\rhoc}为热扩散率，k为热导率，\rho为材料密度，c为比热容。当激光辐照材料时，激光能量作为热源作用于材料表面。假设激光为脉冲激光，其功率密度随时间和空间的分布可以表示为P(x,y,z,t)。在热传导方程中，需要将激光热源项考虑进去，此时热传导方程变为：\rhoc\frac{\partialT}{\partialt}=k\left(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}}\right)+Q(x,y,z,t)其中，Q(x,y,z,t)为激光热源项，它与激光功率密度P(x,y,z,t)的关系为Q(x,y,z,t)=\alphaP(x,y,z,t)，\alpha为材料对激光的吸收系数。对于脉冲激光，其功率密度在时间上通常具有脉冲形式，例如高斯脉冲。在空间上，激光光斑的能量分布也常近似为高斯分布。假设激光光斑在x-y平面上呈高斯分布，其半径为r_0，中心位于(x_0,y_0)，则激光功率密度在空间上的分布可表示为：P(x,y,z,t)=P_0(t)\frac{1}{\pir_0^{2}}\exp\left[-\frac{(x-x_0)^{2}+(y-y_0)^{2}}{r_0^{2}}\right]其中，P_0(t)为随时间变化的激光脉冲峰值功率。为了求解上述热传导方程，通常需要结合具体的边界条件和初始条件。常见的边界条件包括绝热边界条件（材料表面与外界无热交换，即\frac{\partialT}{\partialn}=0，n为表面法向）、对流边界条件（考虑材料表面与周围介质的对流换热，-k\frac{\partialT}{\partialn}=h(T-T_{\infty})，h为对流换热系数，T_{\infty}为周围介质温度）和辐射边界条件（考虑材料表面的热辐射，-k\frac{\partialT}{\partialn}=\varepsilon\sigma(T^{4}-T_{\infty}^{4})，\varepsilon为材料表面发射率，\sigma为斯蒂芬-玻尔兹曼常数）。初始条件一般为材料在激光辐照前的温度分布，通常假设为均匀温度T_0，即T(x,y,z,0)=T_0。利用数学物理方法，如分离变量法、积分变换法等，可以对上述热传导方程进行求解，得到材料内部温度随时间和空间的分布规律。以半无限大固体材料（假设材料在z\geq0区域，z=0为材料表面）在脉冲激光辐照下的温度场分布为例，采用拉普拉斯变换和傅里叶变换相结合的方法进行求解。首先对热传导方程进行拉普拉斯变换，将时间变量t转换为复变量s，得到关于温度T(x,y,z,s)的方程。然后对空间变量x和y进行傅里叶变换，将其转换为波数域变量k_x和k_y，得到一个常微分方程。求解该常微分方程后，再通过逆傅里叶变换和逆拉普拉斯变换，将结果转换回时间和空间域，从而得到材料内部温度场的解析表达式。通过对该解析表达式的分析，可以了解激光辐照过程中材料内部温度的变化趋势，如温度峰值的位置和大小、温度随时间的变化速率以及温度在材料内部的扩散深度等。2.1.3热弹机制下激发超声的点源模型在热弹机制下，构建点源模型有助于深入理解超声的激发原理。当激光脉冲作用于固体材料表面时，若将激光作用区域视为一个点源，可简化对超声激发过程的分析。假设在各向同性的均匀固体材料中，存在一个位于坐标原点(0,0,0)的点热源，其释放的热量随时间的变化为Q(t)。根据热弹理论，材料因吸收热量而产生热膨胀，进而在材料内部激发超声。从力学角度分析，材料的热膨胀会导致内部产生应力。对于各向同性材料，热应力与温度变化之间的关系可由广义胡克定律描述。在小变形情况下，热应力分量\sigma_{ij}与温度变化\DeltaT的关系为：\sigma_{ij}=\lambda\delta_{ij}\theta+2\mu\varepsilon_{ij}-(3\lambda+2\mu)\alpha\DeltaT\delta_{ij}其中，\lambda和\mu为拉梅常数，\delta_{ij}为克罗内克符号，\theta=\frac{\partialu_i}{\partialx_i}为体积应变，\varepsilon_{ij}=\frac{1}{2}(\frac{\partialu_i}{\partialx_j}+\frac{\partialu_j}{\partialx_i})为应变分量，u_i为位移分量，\alpha为材料的热膨胀系数。在点源模型中，由于点热源的作用，材料内部的温度分布具有球对称性。根据热传导理论，可求得点源作用下材料内部的温度分布T(r,t)，其中r=\sqrt{x^{2}+y^{2}+z^{2}}为空间点到点源的距离。当温度分布已知后，可通过上述热应力与温度的关系式，计算出材料内部的热应力分布。在热应力的作用下，材料内部会产生弹性波，即超声。超声的传播满足弹性动力学方程：\rho\frac{\partial^{2}u_i}{\partialt^{2}}=(\lambda+\mu)\frac{\partial\theta}{\partialx_i}+\mu\nabla^{2}u_i将热应力产生的体力项代入弹性动力学方程，结合初始条件和边界条件，可求解出位移分量u_i随时间和空间的变化，从而得到超声在材料中的传播特性。在点源模型中，一些参数对超声激发有着重要影响。激光脉冲的能量是一个关键参数，脉冲能量越大，点源释放的热量Q(t)就越多，材料吸收的能量也越多，从而导致更大的温度变化和热应力，激发的超声幅值也就越大。激光脉冲的持续时间也会影响超声的激发。较短的脉冲持续时间会使热量在极短时间内集中释放，产生更陡峭的温度梯度和热应力变化，有利于激发高频超声；而较长的脉冲持续时间则会使热量释放相对平缓，激发的超声频率相对较低。材料的热膨胀系数\alpha也起着重要作用，热膨胀系数越大，相同温度变化下材料的热膨胀量就越大，产生的热应力也越大，更易于激发超声。此外，材料的密度\rho、拉梅常数\lambda和\mu等弹性参数也会影响超声的传播速度和幅值。密度越大，超声传播速度越慢；拉梅常数则决定了材料的弹性性质，对超声的传播特性有着直接影响。2.1.4热弹机制下激发超声的线源模型与点源模型相比，线源模型在某些情况下能更准确地描述激光激发超声的过程。线源模型适用于激光作用区域可近似为一条线的情况，例如当激光通过柱面透镜聚焦后，在材料表面形成一条细长的光斑，此时可将其视为线源。建立线源模型时，假设线源位于z轴上，长度为L，线源单位长度上释放的热量随时间的变化为q(t)。在各向同性的均匀固体材料中，根据热传导理论，可求解出材料内部的温度分布T(x,y,z,t)。由于线源的对称性，温度分布在垂直于z轴的平面上具有一定的规律性。对于线源激发的超声，同样可通过热弹理论和弹性动力学方程进行分析。线源产生的热应力在材料内部激发弹性波，弹性波的传播满足弹性动力学方程。与点源模型不同的是，线源激发的超声在传播特性上存在一些差异。在点源模型中，超声以球面波的形式向四周传播，能量随着传播距离的增加而逐渐分散，超声的幅值与传播距离r的关系近似为A\propto\frac{1}{r}；而在线源模型中，超声以柱面波的形式传播，能量在柱面方向上分布，其幅值与传播距离r的关系近似为A\propto\frac{1}{\sqrt{r}}。这意味着在线源激发的情况下，超声在传播过程中的衰减相对较慢，在相同距离处，线源激发的超声幅值可能会比点源激发的超声幅值更高。线源的长度L对超声激发和传播特性也有影响。当线源长度较短时，其激发的超声特性可能更接近点源激发的情况；随着线源长度的增加，超声的方向性会增强，在垂直于线源的方向上，超声的能量更加集中。此外，线源单位长度上释放的热量q(t)，即线源的强度，也直接影响着超声的激发。强度越大，激发的超声幅值越大，能够传播的距离也可能更远。在实际应用中，根据具体的激光辐照条件和材料特性，合理选择点源模型或线源模型来分析激光激发超声的过程，有助于更准确地理解和预测超声的产生和传播特性。2.2双层材料中激光声表面波的频散特征2.2.1声表面波在双层固体中的频散方程声表面波在双层固体介质中的传播特性对于理解其在复杂结构中的行为至关重要，基于弹性力学理论，可推导其频散方程。假设双层固体由上层介质和下层介质组成，上层介质的弹性参数为拉梅常数\lambda_1、\mu_1，密度为\rho_1；下层介质的弹性参数为拉梅常数\lambda_2、\mu_2，密度为\rho_2。在各向同性的弹性固体中，弹性波的传播满足Navier方程：\rho\frac{\partial^{2}\vec{u}}{\partialt^{2}}=(\lambda+\mu)\nabla(\nabla\cdot\vec{u})+\mu\nabla^{2}\vec{u}其中，\vec{u}为位移矢量，\rho为密度，\lambda和\mu为拉梅常数，t为时间。对于声表面波，其传播方向通常假设为x方向，在z方向上存在指数衰减。设位移分量为：u_x=A_1e^{-\alpha_1z}e^{i(kx-\omegat)}+A_2e^{\alpha_1z}e^{i(kx-\omegat)}\quad(z\leq0,\text{上层介质})u_x=B_1e^{-\alpha_2z}e^{i(kx-\omegat)}+B_2e^{\alpha_2z}e^{i(kx-\omegat)}\quad(z\geq0,\text{下层介质})u_z=C_1e^{-\alpha_1z}e^{i(kx-\omegat)}+C_2e^{\alpha_1z}e^{i(kx-\omegat)}\quad(z\leq0,\text{上层介质})u_z=D_1e^{-\alpha_2z}e^{i(kx-\omegat)}+D_2e^{\alpha_2z}e^{i(kx-\omegat)}\quad(z\geq0,\text{下层介质})其中，A_i、B_i、C_i、D_i为待定系数，\alpha_1=\sqrt{k^{2}-\frac{\omega^{2}}{c_{1L}^{2}}}，\alpha_2=\sqrt{k^{2}-\frac{\omega^{2}}{c_{2L}^{2}}}，c_{1L}=\sqrt{\frac{\lambda_1+2\mu_1}{\rho_1}}，c_{2L}=\sqrt{\frac{\lambda_2+2\mu_2}{\rho_2}}分别为上层和下层介质的纵波速度，k为波数，\omega为角频率。根据上下层介质分界面z=0处的连续条件，包括位移连续和应力连续：位移连续：u_{x1}(0)=u_{x2}(0)，u_{z1}(0)=u_{z2}(0)应力连续：\sigma_{xz1}(0)=\sigma_{xz2}(0)，\sigma_{zz1}(0)=\sigma_{zz2}(0)其中，应力分量\sigma_{ij}与位移分量的关系为：\sigma_{xz}=\mu(\frac{\partialu_x}{\partialz}+\frac{\partialu_z}{\partialx})\sigma_{zz}=\lambda\nabla\cdot\vec{u}+2\mu\frac{\partialu_z}{\partialz}将位移分量表达式代入上述连续条件，经过一系列的数学推导和化简（利用三角函数关系、指数函数性质等进行整理），可得到一个关于k和\omega的行列式方程，该方程即为声表面波在双层固体中的频散方程：\begin{vmatrix}1-\frac{\alpha_1}{k}\frac{\mu_1}{\mu_2}&1+\frac{\alpha_1}{k}\frac{\mu_1}{\mu_2}&-(1-\frac{\alpha_2}{k})&-(1+\frac{\alpha_2}{k})\\\frac{\alpha_1}{k}&-\frac{\alpha_1}{k}&\frac{\alpha_2}{k}&-\frac{\alpha_2}{k}\\\frac{\lambda_1+2\mu_1}{\mu_1}\alpha_1&-\frac{\lambda_1+2\mu_1}{\mu_1}\alpha_1&\frac{\lambda_2+2\mu_2}{\mu_2}\alpha_2&-\frac{\lambda_2+2\mu_2}{\mu_2}\alpha_2\\\frac{\alpha_1^2}{k^2}\frac{\mu_1}{\mu_2}-\frac{\lambda_1}{\mu_1}&\frac{\alpha_1^2}{k^2}\frac{\mu_1}{\mu_2}-\frac{\lambda_1}{\mu_1}&\frac{\alpha_2^2}{k^2}-\frac{\lambda_2}{\mu_2}&\frac{\alpha_2^2}{k^2}-\frac{\lambda_2}{\mu_2}\end{vmatrix}=0在这个频散方程中，\omega表示声表面波的角频率，它决定了波的振荡快慢，与波的能量和传播特性密切相关，角频率越高，波的振荡越快，携带的能量也相对较高。k为波数，它反映了波在空间中的变化快慢，波数越大，波在单位长度内的变化越剧烈，波长\lambda=\frac{2\pi}{k}越短。\lambda_1、\lambda_2、\mu_1、\mu_2是上下层介质的拉梅常数，它们描述了介质的弹性性质，拉梅常数越大，介质越不容易发生形变，对声表面波的传播速度和特性有着重要影响。\rho_1、\rho_2分别为上下层介质的密度，密度影响声表面波的传播速度，密度越大，声表面波传播速度越慢。这些参数共同决定了声表面波在双层固体中的传播特性，通过求解频散方程，可以得到不同频率下声表面波的波数，进而分析其传播特性。2.2.2声表面波在双层固体中的频散曲线为了深入了解声表面波在双层固体中的传播特性，通过数值计算方法绘制频散曲线。利用Matlab等数学软件，给定双层固体介质的具体参数，如上层介质为铝，其拉梅常数\lambda_1=6.13\times10^{10}N/m^2，\mu_1=2.61\times10^{10}N/m^2，密度\rho_1=2700kg/m^3；下层介质为钢，拉梅常数\lambda_2=1.18\times10^{11}N/m^2，\mu_2=8.16\times10^{10}N/m^2，密度\rho_2=7850kg/m^3。在数值计算过程中，采用迭代算法求解频散方程。首先设定一个初始的角频率\omega值，然后通过迭代计算寻找满足频散方程的波数k。例如，使用牛顿迭代法，根据频散方程构建函数f(k,\omega)，对k求偏导数\frac{\partialf}{\partialk}，通过迭代公式k_{n+1}=k_n-\frac{f(k_n,\omega)}{\frac{\partialf}{\partialk}(k_n,\omega)}不断更新k的值，直到满足一定的收敛条件，如\vertk_{n+1}-k_n\vert\lt10^{-6}。通过上述数值计算方法，得到不同角频率\omega对应的波数k，进而计算出声表面波的相速度v_p=\frac{\omega}{k}。以角频率\omega为横坐标，相速度v_p为纵坐标，绘制频散曲线。从频散曲线可以看出，在低频段，声表面波的相速度接近下层介质的瑞利波速度，这是因为低频声表面波的波长较长，其能量主要分布在下层介质中，传播特性受下层介质影响较大。随着频率的增加，相速度逐渐减小，这是由于高频声表面波的波长较短，能量更多地集中在上层介质，而上层介质的弹性性质和密度与下层不同，导致相速度发生变化。在某些特定频率处，频散曲线可能会出现转折点或奇异点，这与双层介质的结构和参数有关，反映了声表面波在传播过程中能量的重新分布和模式转换。在分析不同频率下声表面波的衰减特性时，根据位移分量表达式中的指数衰减项，声表面波在z方向上的振幅随深度z的增加而衰减，衰减系数与波数k和角频率\omega有关。通过计算不同频率下的衰减系数，可以发现高频声表面波的衰减速度较快，这是因为高频波的能量更集中在表面附近，更容易受到表面粗糙度、介质不均匀性等因素的影响，导致能量快速耗散。而低频声表面波由于波长较长，能量分布较深，衰减相对较慢。在实际应用中，如在早期龋齿和牙隐裂的检测中，了解声表面波的频散特性和衰减特性，有助于选择合适的检测频率，提高检测的准确性和灵敏度。2.3激光声表面波的检测方法2.3.1电学检测法电学检测法是利用压电传感器等电学元件来检测激光声表面波的一种常用方法。其原理基于压电效应，即某些材料在受到机械应力作用时，会在其表面产生电荷，电荷量与所受应力成正比；反之，当这些材料受到电场作用时，会发生形变。在激光声表面波检测中，压电传感器通常采用压电陶瓷或压电晶体材料制成。当激光激发的声表面波传播到压电传感器表面时，会使压电传感器产生机械振动，进而在其表面产生电荷，通过检测这些电荷信号，就可以获取声表面波的相关信息。该检测方法的装置结构相对较为简单。一般由脉冲激光器、光学聚焦系统、被检测样品、压电传感器以及信号采集与处理系统组成。脉冲激光器发射的激光经过光学聚焦系统聚焦后，照射到被检测样品表面，激发出声表面波。压电传感器被放置在靠近样品表面的位置，以接收声表面波引起的机械振动。信号采集与处理系统则负责采集压电传感器输出的电信号，并对其进行放大、滤波、数字化等处理。在实际应用中，为了提高检测的灵敏度和准确性，通常会对压电传感器进行优化设计。例如，选择合适的压电材料，不同的压电材料具有不同的压电系数，压电系数越大，在相同应力作用下产生的电荷量就越多，检测灵敏度也就越高。同时，还会对压电传感器的尺寸和形状进行优化，使其与声表面波的波长和传播特性相匹配，以提高对声表面波的接收效率。信号处理方法对于准确获取声表面波的信息至关重要。在信号采集过程中，由于受到各种噪声的干扰，采集到的电信号往往包含大量的噪声成分。因此，需要采用滤波技术去除噪声。常见的滤波方法有低通滤波、高通滤波和带通滤波等。低通滤波可以去除高频噪声，保留低频信号；高通滤波则相反，用于去除低频噪声，保留高频信号；带通滤波则可以选择特定频率范围内的信号，去除其他频率的噪声。此外，还可以采用数字滤波算法，如有限脉冲响应（FIR）滤波器和无限脉冲响应（IIR）滤波器等，对信号进行更精确的滤波处理。在数据处理方面，常常会采用傅里叶变换、小波变换等方法对信号进行分析。傅里叶变换可以将时域信号转换为频域信号，通过分析频域信号的频谱特征，能够获取声表面波的频率成分和幅值信息。小波变换则具有良好的时频局部化特性，能够同时在时域和频域对信号进行分析，对于检测声表面波信号中的瞬态变化和微弱信号具有独特的优势。通过对信号的处理和分析，可以得到声表面波的传播速度、频率、振幅等参数，从而用于早期龋齿和牙隐裂的诊断与评估。2.3.2光学非干涉检测法光学非干涉检测法是基于光散射、光吸收等原理的一类激光声表面波检测技术。在光散射原理中，当激光照射到样品表面时，声表面波会引起样品表面的微观结构发生周期性变化，这种变化会导致光的散射特性发生改变。例如，在瑞利散射中，散射光的强度与散射角度、声表面波的频率和振幅等因素有关。通过检测散射光的强度、角度分布或偏振状态等信息，就可以间接获取声表面波的相关参数。当声表面波的频率发生变化时，散射光的强度在不同散射角度下也会相应改变，通过精确测量这些变化，能够推算出声表面波的频率。基于光吸收原理的检测方法则是利用声表面波对材料光学吸收特性的影响来进行检测。声表面波在传播过程中会与材料分子相互作用，导致材料的局部温度、密度等物理性质发生变化，进而影响材料对光的吸收能力。通过测量材料对特定波长光的吸收系数变化，就可以检测到声表面波的存在及其特性。在某些材料中，当声表面波传播时，会使材料分子的振动加剧，导致材料对特定波长光的吸收增强，通过精确测量吸收系数的变化，能够确定声表面波的传播情况。这类检测方法具有一些显著的优点。由于不需要复杂的干涉光路，其系统结构相对简单，成本较低，易于搭建和操作。对检测环境的要求相对较低，受外界干扰的影响较小，具有较好的稳定性。然而，光学非干涉检测法也存在一些缺点。其检测灵敏度相对较低，对于微弱的声表面波信号，可能难以准确检测和分析。在检测精度方面，由于光散射和光吸收过程受到多种因素的影响，如材料的不均匀性、表面粗糙度等，导致检测结果的准确性和重复性相对较差。在适用场景方面，光学非干涉检测法适用于对检测精度要求不是特别高，但对检测成本和操作便捷性有较高要求的场合。在一些工业生产线上，需要对大量产品进行快速检测，以初步筛选出可能存在缺陷的产品，此时光学非干涉检测法可以发挥其快速、低成本的优势。对于一些对检测环境要求较为苛刻的现场检测，如高温、高湿度等环境，光学非干涉检测法因其对环境适应性强的特点，也能较好地满足检测需求。2.3.3光学干涉检测法光学干涉检测法是利用干涉原理来检测激光声表面波的一类重要方法，常见的干涉检测方法包括迈克尔逊干涉和马赫-曾德尔干涉等，在激光声表面波检测中有着广泛的应用。迈克尔逊干涉仪由光源、分光镜、两个反射镜和探测器组成。光源发出的光经分光镜分为两束，一束光照射到参考反射镜后反射回来，另一束光照射到被检测样品表面，样品表面的声表面波会使反射光的相位发生变化，两束反射光在分光镜处再次汇合发生干涉，产生干涉条纹。探测器用于检测干涉条纹的变化，通过分析干涉条纹的移动、变形等信息，就可以获取声表面波的相位变化，进而得到声表面波的传播特性。如果声表面波使样品表面产生微小位移，导致反射光的相位发生改变，干涉条纹就会发生移动，通过精确测量干涉条纹的移动距离，能够计算出声表面波引起的样品表面位移，从而推断出声表面波的传播情况。马赫-曾德尔干涉仪则由光源、两个分光镜、两个反射镜和探测器组成。光源发出的光经第一个分光镜分为两束，一束作为参考光，经过一系列反射镜后到达第二个分光镜；另一束作为信号光，照射到被检测样品表面，受声表面波影响后也到达第二个分光镜。在第二个分光镜处，参考光和信号光发生干涉，探测器检测干涉光的强度变化。由于声表面波会改变信号光的相位，从而导致干涉光强度发生变化，通过分析干涉光强度的变化规律，就可以获得声表面波的相关信息。当声表面波使信号光的相位发生改变时，干涉光的强度会在极大值和极小值之间变化，通过精确测量干涉光强度的变化，能够确定声表面波的相位变化，进而分析出声表面波的特性。在激光声表面波检测中，这些干涉检测方法具有很高的检测灵敏度和精度。能够检测到极其微小的声表面波引起的相位变化，从而准确获取声表面波的传播特性，对于早期龋齿和牙隐裂的检测具有重要意义。由于其检测原理基于光的干涉，对环境因素较为敏感，如温度、振动等外界因素的变化可能会导致干涉条纹的漂移，影响检测结果的准确性。因此，在实际应用中，通常需要采取一些措施来减小环境因素的影响，如采用稳定的光学平台、对光学元件进行温度控制等。2.4本章小结本章深入探究了激光声表面波的激发检测机理及其频散特征。在激光激发超声的热弹机制方面，详细阐述了激光的吸收及其辐照效应，分析了热弹机制下激光辐照材料的温度场分布，分别建立了点源模型和线源模型，探讨了激光脉冲能量、持续时间以及材料热膨胀系数等参数对超声激发的影响。在双层材料中激光声表面波的频散特征研究中，基于弹性力学理论推导了声表面波在双层固体中的频散方程，并通过数值计算绘制了频散曲线，分析了不同频率下声表面波的相速度和衰减特性。在激光声表面波的检测方法上，介绍了电学检测法、光学非干涉检测法和光学干涉检测法，分别阐述了它们的工作原理、装置结构、信号处理方法以及各自的优缺点和适用场景。这些理论和技术要点为后续利用激光声表面波进行早期龋齿和牙隐裂的诊断与评估研究奠定了坚实的基础。三、牙齿形貌结构对声表面波传播特性的影响3.1激光超声用于牙齿临床检测可行性的有限元模拟3.1.1激光辐照牙齿的有限元模型及热弹性理论为深入探究激光超声在牙齿临床检测中的可行性，建立精确的有限元模型至关重要。牙齿结构复杂，包含釉质、牙本质、牙髓等多层结构，各层结构的材料属性存在显著差异。釉质是人体中最硬的组织，主要由无机物（约96%-97%）组成，其密度约为3.0g/cm³，弹性模量高达40-80GPa。牙本质则相对较软，无机物含量约为70%，密度约为2.0g/cm³，弹性模量在18-20GPa之间。牙髓为软组织，主要包含血管、神经和结缔组织，其力学性能与釉质和牙本质有很大不同。在有限元模型中，将牙齿各层结构视为均匀、连续且各向同性的材料。采用SolidWorks等三维建模软件，依据牙齿的解剖学结构，构建包含釉质、牙本质和牙髓的完整牙齿三维模型。在建模过程中，精确测量牙齿各部分的尺寸，确保模型的几何形状与实际牙齿相符。对于釉质，根据其在牙齿表面的分布特点，准确描绘其厚度和形状变化；对于牙本质，考虑其与釉质的衔接关系以及牙髓腔的形态，构建出逼真的模型。利用有限元分析软件ANSYS，对构建好的模型进行网格划分，采用四面体单元或六面体单元，根据模型的复杂程度和计算精度要求，合理控制单元尺寸。在牙齿的关键部位，如釉质与牙本质的交界处、牙髓腔附近等，适当减小单元尺寸，以提高计算精度。热弹性理论在模拟激光辐照牙齿过程中起着核心作用。当激光脉冲作用于牙齿表面时，牙齿吸收激光能量，根据热传导理论，热量会在牙齿内部传播，导致温度分布发生变化。热传导方程为：\rhoc\frac{\partialT}{\partialt}=k\nabla^{2}T+Q其中，\rho为材料密度，c为比热容，T为温度，t为时间，k为热导率，Q为激光热源项。激光热源项Q与激光功率密度P的关系为Q=\alphaP，\alpha为材料对激光的吸收系数。在实际计算中，根据激光的脉冲特性和光斑分布，确定激光功率密度P的时间和空间分布函数。假设激光脉冲为高斯脉冲，其功率密度在空间上呈高斯分布，通过准确设定高斯函数的参数，如峰值功率、光斑半径等，来描述激光功率密度的分布。在热弹性理论中，温度变化会引起材料的热膨胀，从而产生热应力。热应力与温度变化之间的关系可通过热弹性本构方程描述：\sigma_{ij}=\lambda\delta_{ij}\theta+2\mu\varepsilon_{ij}-(3\lambda+2\mu)\alpha\DeltaT\delta_{ij}其中，\sigma_{ij}为应力分量，\lambda和\mu为拉梅常数，\delta_{ij}为克罗内克符号，\theta为体积应变，\varepsilon_{ij}为应变分量，\alpha为热膨胀系数，\DeltaT为温度变化。通过求解热传导方程得到温度分布后，代入热弹性本构方程，即可计算出牙齿内部的热应力分布。在计算过程中，考虑到牙齿各层结构材料属性的差异，分别设定不同的拉梅常数、热膨胀系数等参数，以准确模拟热应力在各层结构中的产生和传播。3.1.2热弹耦合的有限元方程及网格大小与时间步长的选取热弹耦合问题涉及热传导和弹性力学两个物理场的相互作用，其有限元方程的建立基于能量守恒原理和虚功原理。在热传导方面，根据热传导方程，结合有限元方法的离散化思想，将牙齿模型划分为有限个单元，对每个单元内的温度场进行插值近似。假设单元内的温度分布为T^e(x,y,z,t)，可表示为节点温度T_i(t)（i=1,2,\cdots,n，n为单元节点数）的插值函数，即T^e(x,y,z,t)=\sum_{i=1}^{n}N_i(x,y,z)T_i(t)，其中N_i(x,y,z)为形函数。将温度插值函数代入热传导方程，通过加权余量法得到热传导的有限元方程：[C]\{\dot{T}\}+[K]\{T\}=\{Q\}其中，[C]为热容矩阵，[K]为热传导矩阵，\{\dot{T}\}为节点温度对时间的一阶导数向量，\{T\}为节点温度向量，\{Q\}为节点热流向量。在弹性力学方面，根据弹性动力学方程，同样对位移场进行离散化。设单元内的位移分量为u^e(x,y,z,t)、v^e(x,y,z,t)、w^e(x,y,z,t)，分别表示x、y、z方向的位移，可表示为节点位移u_i(t)、v_i(t)、w_i(t)（i=1,2,\cdots,n）的插值函数。将位移插值函数代入弹性动力学方程，通过加权余量法得到弹性力学的有限元方程：[M]\{\ddot{u}\}+[K_d]\{u\}=\{F\}+\{F_T\}其中，[M]为质量矩阵，[K_d]为刚度矩阵，\{\ddot{u}\}为节点加速度向量，\{u\}为节点位移向量，\{F\}为外力向量，\{F_T\}为热应力引起的等效节点力向量。热弹耦合通过热应力项\{F_T\}实现，\{F_T\}与温度变化\DeltaT相关。将热传导的有限元方程和弹性力学的有限元方程联立，即可得到热弹耦合的有限元方程。网格大小和时间步长的选取对模拟精度和计算效率有着重要影响。较小的网格尺寸能够更精确地描述牙齿的几何形状和物理场变化，但会增加计算量和计算时间。在选择网格大小时，需要综合考虑模型的复杂程度和计算精度要求。对于牙齿的复杂结构部位，如釉质与牙本质的过渡区域、牙髓腔的弯曲部分等，采用较小的网格尺寸，以准确捕捉物理量的变化。而在一些相对规则的区域，可以适当增大网格尺寸，以提高计算效率。通过多次试验，确定合适的网格尺寸，使得在保证计算精度的前提下，尽量减少计算资源的消耗。时间步长的选择也至关重要。较小的时间步长可以提高计算的稳定性和精度，但会增加计算步数和计算时间。时间步长过大可能导致计算结果不稳定，甚至发散。在实际模拟中，根据激光脉冲的持续时间和热弹耦合过程的特征时间尺度，选择合适的时间步长。通常采用Courant-Friedrichs-Lewy（CFL）条件来确定时间步长的上限，即\Deltat\leqslant\frac{h}{c}，其中\Deltat为时间步长，h为最小网格尺寸，c为材料中的声速。通过调整时间步长，观察计算结果的稳定性和收敛性，最终确定满足计算要求的时间步长。3.1.3不同波长激光辐照牙齿的温度场和超声位移利用建立的有限元模型，模拟不同波长激光辐照牙齿时的温度场和超声位移变化。常见的用于牙齿检测的激光波长有1064nm的Nd:YAG激光、2940nm的Er:YAG激光等。不同波长的激光在牙齿组织中的吸收特性不同，这是由牙齿各层结构的化学成分和微观结构决定的。釉质主要由羟基磷灰石晶体组成，对2940nm的Er:YAG激光具有较高的吸收系数，因为该波长与釉质中的水分子吸收峰相匹配，激光能量能够被釉质中的水分子有效吸收，从而在釉质表层产生较高的能量密度和温度升高。而对于1064nm的Nd:YAG激光，釉质的吸收系数相对较低，激光能量更容易穿透釉质，在牙本质中被吸收。在模拟1064nm激光辐照牙齿时，由于其在釉质中的吸收较少，大部分能量穿透釉质到达牙本质。在牙本质中，激光能量被吸收转化为热能，导致牙本质温度升高。温度场分布呈现出从牙本质表面向内部逐渐降低的趋势，在牙本质与牙髓的交界处，温度变化相对平缓。根据模拟结果，在激光辐照后的短时间内，牙本质表面温度可升高至几十摄氏度。随着时间的推移，热量逐渐向牙髓和周围组织扩散。当激光波长为2940nm时，由于釉质对该波长激光的高吸收特性，能量主要集中在釉质表层。釉质表面温度迅速升高，在极短时间内可达到较高温度，甚至可能导致釉质局部熔化或汽化。温度场分布在釉质层内呈现出陡峭的梯度，从釉质表面向内部迅速降低。在釉质与牙本质交界处，温度急剧下降。激光辐照引起的温度变化会导致牙齿材料的热膨胀，进而产生超声位移。在1064nm激光辐照下，由于牙本质温度升高产生的热膨胀，在牙本质内部和釉质与牙本质交界处会产生超声位移。超声位移的方向和大小与温度分布、材料的热膨胀系数以及弹性模量等因素有关。通过有限元模拟，可以得到超声位移在牙齿内部的分布情况。在牙本质中，超声位移沿着激光辐照方向和垂直方向都有一定的分量，且在牙本质与牙髓交界处，由于材料属性的差异，超声位移会发生一定的变化。对于2940nm激光辐照，由于釉质表面温度的急剧升高，釉质的热膨胀产生的超声位移主要集中在釉质层。釉质表面的超声位移较大，随着深度的增加，超声位移迅速减小。在釉质与牙本质交界处，超声位移也会发生突变。通过对不同波长激光辐照下牙齿温度场和超声位移的模拟分析，可以发现波长与温度、位移之间存在密切关系。较短波长的激光（如2940nm的Er:YAG激光）主要作用于牙齿的表层，产生较高的表面温度和较大的表面超声位移；而较长波长的激光（如1064nm的Nd:YAG激光）能够穿透牙齿表层，在深层组织中产生温度变化和超声位移。这些模拟结果为选择合适的激光波长进行牙齿检测提供了理论依据。3.2牙齿形貌结构对激光声表面波频域特性的影响3.2.1不同形貌结构的牙齿中传播的声表面波牙齿的形貌结构复杂多样，切牙和磨牙作为两种典型的牙齿类型，在形态、功能等方面存在显著差异，这使得声表面波在它们之中的传播特性也截然不同。切牙位于口腔前部，其主要功能是切割食物。切牙的牙冠呈扁平状，唇舌面近似梯形，邻面为三角形，牙根一般为单根且较为细长。在这样的结构中，声表面波从激光激发点开始传播，其传播路径相对较为规则。当激光在切牙唇面激发声表面波后，声表面波会沿着唇面、邻面和舌面等部位传播，由于切牙各面的几何形状相对简单，声表面波在传播过程中遇到的界面相对较少，反射和散射现象相对不那么复杂。在切牙唇面传播时，声表面波主要以沿表面的平面波形式传播，其波形特征表现为较为规则的正弦波或余弦波形式，波峰和波谷的变化相对平稳。磨牙则位于口腔后部，主要负责研磨食物，其牙冠体积较大，形态复杂，具有多个牙尖和窝沟。磨牙的表面曲率变化较大，且存在复杂的咬合面结构。这些复杂的形貌结构导致声表面波在磨牙中的传播路径变得错综复杂。当声表面波在磨牙表面传播时，遇到牙尖、窝沟等结构，会发生强烈的反射和散射现象。在牙尖处，声表面波会发生反射，一部分能量会沿着原路径返回，另一部分能量则会向不同方向散射，导致波形发生严重畸变。窝沟部位由于其特殊的几何形状，会使声表面波在其中多次反射和散射，形成复杂的干涉图样。在磨牙咬合面传播的声表面波，其波形呈现出不规则的振荡，波峰和波谷的分布不均匀，存在多个峰值和谷值，且这些峰值和谷值的大小和位置会随着传播距离的变化而不断改变。为了更直观地比较切牙和磨牙中声表面波传播特性的差异，通过有限元模拟进行分析。在模拟中，构建精确的切牙和磨牙三维模型，设定相同的激光激发条件，观察声表面波的传播过程。从模拟结果可以看出，在相同的传播距离下，切牙中声表面波的传播速度相对较为稳定，而磨牙中由于复杂的结构导致声表面波传播速度存在明显的波动。在频谱分析方面，切牙中声表面波的频谱相对较为集中，主要频率成分相对较少；而磨牙中由于反射和散射的影响，声表面波的频谱更加宽泛，包含了更多的频率成分。3.2.2超声信号的二维傅里叶变换时频分析方法二维傅里叶变换是一种强大的信号分析工具，在超声信号时频分析中具有重要应用。其原理基于傅里叶变换的基本思想，将信号从时域转换到频域，以揭示信号的频率成分。对于一个二维超声信号s(x,t)，其中x表示空间位置，t表示时间，二维傅里叶变换的数学表达式为：S(f_x,f_t)=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}s(x,t)e^{-i2\pi(f_xx+f_tt)}dxdt其中，S(f_x,f_t)是二维傅里叶变换后的频域信号，f_x和f_t分别是空间频率和时间频率。在超声信号时频分析中，应用二维傅里叶变换的步骤如下。对采集到的超声信号进行预处理，包括去除噪声、滤波等操作，以提高信号的质量。将超声信号按照空间位置和时间进行离散化处理，得到离散的二维信号矩阵。利用快速傅里叶变换（FFT）算法对离散信号矩阵进行二维傅里叶变换计算。快速傅里叶变换算法能够大大提高计算效率，将原本复杂的傅里叶变换计算简化为一系列快速的运算。在实际计算中，可使用MATLAB等数学软件提供的FFT函数进行计算，通过调用相应的函数并传入离散信号矩阵，即可快速得到二维傅里叶变换后的频域数据。分析结果通常以时频图的形式表达。时频图的横轴表示时间，纵轴表示频率，图中的每个点表示在特定时间和频率下信号的幅值或能量。通过观察时频图，可以清晰地看到超声信号在不同时间点的频率组成和变化情况。在早期龋齿或牙隐裂的检测中，时频图能够直观地显示出由于牙齿病变导致的超声信号频率成分的改变。如果牙齿存在早期龋齿，病变部位的超声信号在时频图上可能会出现特定频率段的能量增强或减弱，通过分析这些特征，可以判断龋齿的存在和程度。3.2.3切牙模型中牙釉质厚度对声表面波频域特性的影响切牙模型中，牙釉质作为牙齿最外层的硬组织，其厚度变化对声表面波频域特性有着显著影响。利用有限元模拟和实验相结合的方法，深入研究这一影响机制。在有限元模拟中，构建不同牙釉质厚度的切牙模型。设定牙釉质厚度分别为0.5mm、1.0mm和1.5mm，其他结构参数保持一致。通过模拟激光在切牙表面激发声表面波的过程，分析不同牙釉质厚度下声表面波的频域特性。从模拟结果可知，随着牙釉质厚度的增加，声表面波的主要频率成分向低频方向移动。当牙釉质厚度为0.5mm时，声表面波的主要频率集中在10MHz-15MHz范围内；当牙釉质厚度增加到1.0mm时，主要频率范围变为8MHz-12MHz；而当牙釉质厚度达到1.5mm时，主要频率范围进一步降低到6MHz-10MHz。这是因为牙釉质厚度增加，声表面波在其中传播的路径变长，波的传播速度相对较慢，根据波速与频率的关系v=f\lambda（其中v为波速，f为频率，\lambda为波长），在波长不变的情况下，波速降低导致频率降低。在能量分布方面，随着牙釉质厚度的增加，声表面波在低频段的能量占比逐渐增大。通过对模拟结果中不同频率段的能量进行积分计算，得到不同牙釉质厚度下声表面波在低频段（0MHz-8MHz）、中频段（8MHz-12MHz）和高频段（12MHz-16MHz）的能量占比。当牙釉质厚度为0.5mm时，低频段能量占比约为30\%，中频段能量占比约为50\%，高频段能量占比约为20\%；当牙釉质厚度增加到1.0mm时，低频段能量占比增加到约40\%，中频段能量占比降低到约45\%，高频段能量占比降低到约15\%；当牙釉质厚度达到1.5mm时，低频段能量占比进一步增加到约50\%，中频段能量占比降低到约35\%，高频段能量占比降低到约15\%。这表明牙釉质厚度的增加使得声表面波的能量更多地集中在低频段。为了验证模拟结果，进行实验研究。选取不同牙釉质厚度的离体切牙样本，利用激光声表面波检测系统采集声表面波信号，并对信号进行频域分析。实验结果与模拟结果基本一致，进一步证实了牙釉质厚度对声表面波频域特性的影响规律。3.2.4磨牙模型中表面曲率及牙釉质厚度对声表面波频域特性的影响磨牙模型中，表面曲率和牙釉质厚度是影响声表面波频域特性的两个重要因素，它们之间存在协同作用。利用数值模拟和实验研究相结合的方法，深入分析这种协同作用机制。在数值模拟方面，构建一系列不同表面曲率和牙釉质厚度的磨牙模型。通过改变磨牙表面的曲率半径和牙釉质厚度，模拟激光激发声表面波的传播过程。当表面曲率半径为5mm、牙釉质厚度为1.0mm时，分析声表面波的频域特性。随着表面曲率半径的减小，即表面曲率增大，声表面波的散射现象加剧，导致其频谱变得更加复杂。原本相对集中的频率成分变得更加分散，高频成分增加。这是因为表面曲率增大，声表面波在传播过程中遇到的反射和散射界面增多，波的传播方向不断改变，不同方向的波相互干涉，使得频谱展宽。在牙釉质厚度变化的情况下，与表面曲率相互作用。当表面曲率一定时，随着牙釉质厚度的增加，声表面波的主要频率成分向低频方向移动。在表面曲率半径为5mm时，牙釉质厚度从1.0mm增加到1.5mm，主要频率范围从8MHz-12MHz降低到6MHz-10MHz。这与切牙模型中牙釉质厚度对声表面波频率的影响类似，是由于声表面波在较厚的牙釉质中传播路径变长，波速降低导致频率降低。而当牙釉质厚度增加的同时表面曲率也增大时，这种频率向低频移动的趋势会更加明显，且频谱的复杂性进一步增加。在表面曲率半径为3mm、牙釉质厚度从1.0mm增加到1.5mm时，不仅主要频率范围进一步降低到5MHz-9MHz，而且频谱中出现了更多的杂散频率成分，这是由于表面曲率增大带来的散射增强与牙釉质厚度增加导致的波速变化相互叠加的结果。通过实验进一步验证上述模拟结果。选取不同表面曲率和牙釉质厚度的离体磨牙样本，利用激光声表面波检测系统进行检测。对采集到的声表面波信号进行二维傅里叶变换时频分析，得到时频图。实验结果显示，随着表面曲率增大和牙釉质厚度增加，时频图中频率成分的分布和变化与模拟结果相符，充分说明了表面曲率和牙釉质厚度对声表面波频域特性的协同影响。3.3牙齿形貌结构对声表面波传播特性影响的实验研究3.3.1不同形貌结构的牙齿样品及声表面波的激发探测实验系统为深入探究牙齿形貌结构对声表面波传播特性的影响，精心选取不同形貌结构的牙齿样品。牙齿样品主要来源于因正畸治疗需要拔除的健康离体牙齿，以及因牙周病等原因拔除但牙体结构基本完整的牙齿。在样品选取过程中，严格筛选，确保牙齿无明显龋齿、裂纹等病变，且牙釉质和牙本质结构完整。对于因正畸拔除的牙齿，多为前磨牙和磨牙，这些牙齿在形态上具有典型的特征，前磨牙牙冠呈立方形，有颊尖和舌尖，牙根一般为单根或双根；磨牙牙冠体积较大，具有多个牙尖和复杂的窝沟结构，牙根多为多根。因牙周病拔除的牙齿，在去除牙周组织后，保留了完整的牙体硬组织，可用于研究声表面波在不同牙根形态和牙冠结构下的传播特性。在制备牙齿样品时，首先对牙齿进行清洁处理。使用超声波清洗机，将牙齿浸泡在适量的去离子水中，超声清洗15-20分钟，去除牙齿表面的牙结石、菌斑和其他杂质。清洗后，用蒸馏水冲洗干净，并用无菌纱布轻轻擦干。对于需要进行特殊处理的样品，如研究牙釉质厚度对声表面波传播特性的影响时，采用化学蚀刻的方法，使用一定浓度的氢氟酸溶液对牙釉质表面进行蚀刻，控制蚀刻时间和浓度，制备出不同牙釉质厚度的牙齿样品。在蚀刻过程中，严格按照操作规程进行，确保蚀刻的均匀性和准确性。搭建高精度的声表面波激发和探测实验装置。该装置主要由脉冲激光器、光学聚焦系统、激光干涉仪、信号采集与处理系统以及三维移动平台等部分组成。脉冲激光器选用高能量的Nd:YAG脉冲激光器，其波长为1064nm，脉冲宽度为10ns，重复频率为10Hz，能够提供足够的能量激发声表面波。光学聚焦系统采用一组高质量的透镜，将激光束聚焦到牙齿表面，光斑直径可控制在50-100μm之间，确保激光能量能够集中作用于牙齿表面的微小区域，从而有效激发声表面波。激光干涉仪选用迈克尔逊干涉仪，其具有高灵敏度和高精度的特点，能够检测到极其微小的声表面波引起的位移变化，分辨率可达纳米量级。信号采集与处理系统包括高速数据采集卡和专业的信号处理软件。高速数据采集卡的采样频率为1GHz，能够快速准确地采集激光干涉仪输出的信号。信号处理软件采用MATLAB编写，具备滤波、降噪、傅里叶变换等多种信号处理功能，能够对采集到的信号进行分析和处理，提取出声表面波的传播特性参数。三维移动平台用于精确控制牙齿样品的位置和姿态，其定位精度可达1μm，能够实现牙齿样品在三个方向上的精确移动，以便对不同位置的声表面波进行测量。在搭建实验装置时，确保各部分组件的安装精度和稳定性，减少外界干扰对实验结果的影响。3.3.2不同形貌结构牙齿中传播的声表面波及其频域特性通过上述实验装置，对不同形貌结构牙齿中传播的声表面波进行了精确测量。在测量过程中，将牙齿样品固定在三维移动平台上，调整平台位置，使激光束垂直照射到牙齿表面的指定位置。脉冲激光器发射的激光脉冲作用于牙齿表面，激发出声表面波。激光干涉仪实时检测声表面波引起的牙齿表面位移变化，并将检测到的信号传输到信号采集与处理系统。对于切牙样品，从实验测量得到的声表面波波形来看，在传播初期，波形较为规则，呈现出近似正弦波的形状。随着传播距离的增加，波形逐渐发生变化，波峰和波谷的幅度略有减小，这是由于声表面波在传播过程中能量逐渐衰减。通过对切牙中声表面波的频域特性进行分析，利用傅里叶变换将时域信号转换为频域信号，得到其频谱图。频谱图显示，切牙中声表面波的主要频率成分集中在8-12MHz范围内，在该频率范围内，信号幅值较高，能量较为集中。这与理论分析和模拟结果相符，表明在切牙结构中，声表面波的传播特性相对较为稳定，主要频率成分较为集中。对于磨牙样品，实验测量得到