体外冲击波碎石机(Extracorporeal shock wave lithotripter,ESWL)是治疗结石的主要手段,现阶段对其微观工作机理、动态声场分布等问题的研究尚处于起步阶段。通过对ESWL声场的建模及仿真可以对上述性能的研究带来很大的帮助。应用非线性传播方程建立了ESWL声场的数学模型,并通过有限差分算法结合MATLAB软件的图形模块进行了数值仿真。模型不仅考虑了非线性的影响,还引入了衍射、和热粘性吸收对声波的作用,最大限度地与实际声场所受到的各种影响保持一致。仿真很好地反映出ESWL产生的高强度超声波传播的非线性及一些与线性声波传播显著不同的特性。对形象地了解ESWL声波的传播过程和ESWL反射罩优化设计有很好的参考价值。