2025《激光热弹效应基本理论概述》2600字.docxVIP

激光热弹效应基本理论概述

目录

TOC\o1-3\h\u4078激光热弹效应基本理论概述 1

270991热弹性理论 1

246352热弹耦合方程的有限元求解理论基础 2

276663热弹耦合方程的有限元模型的建立和验证 3

1热弹性理论

当激光源对钢轨材料进行辐照时，激光能量会被钢轨表面瞬间吸收，温度升高，进而在表面产生热膨胀，如图1.1所示。在温度升高的过程中，钢轨材料表面产生的温度场是非均匀的，这种不均匀会导致应力、应变在钢轨中产生，进而产生超声波。通常，Gauss分布能较好地解释激光能量在空间上的呈现，结合圆柱坐标系，对各向同性的钢轨材料开展热弹分析，其中，热扩散方程为：

(1.1)

其中：为激光源在t时刻的温度，k为钢轨材料的热扩散系数，c为钢轨材料热容量，为和材料密度。

图2SEQ图2\*ARABIC1激光源对钢轨材料辐照示意图

钢轨材料在被激光源辐照后，会在钢轨表面产生热膨胀效应。在此过程中，钢轨表面加载的外部热源随时间和空间也都呈现高斯分布，激光超声波模型就是以这种热辐照效应作为边界条件从而建立的。在二维坐标系中，设钢轨的上表面厚度为d，径向宽度为R，进一步可以得到钢轨表面所施加的热流边界条件：

(1.2)

其中，为材料的表面吸收率；为激光源的峰值能量密度；激光束空间分布函数和时间分布函数的具体表达式为：

(1.3)

(1.4)

其中，为激光点源半经；为激光源上升时间。钢轨材料的初始温度值为300k。

激光热弹机制下，钢轨中产生的超声波信号满足Navier-Stokes方程：

(1.5)

式中，代表瞬态位移向量；代表瞬态体力源；及分别为材料的常数；代表钢轨热弹效应的耦合系数，公式为：

(1.6)

式中，表示热膨胀系数。另外，钢轨材料需要满足相应的自由边界条件，即上表面（Z=0）和下表面（Z=h）为零：

(1.7)

式中，n为与钢轨面垂直的单位向量；σ和I分别为应力张量和单位张量。

在满足上下面自由边界条件后，还需要约束模型两侧的边界位移，同样地，设置边界位移等于零。另外，模型的起始位移和起始速度应符合：

(1.8)

2热弹耦合方程的有限元求解理论基础

在热弹性机理下，激光光源与钢轨表面热传导的有限元公式为：

(1.9)

式中，和分别表示热传导矩阵和温度矩阵；和分别表示热流矢量和热源矢量；和分别表示温度变化率以及热容量矩阵。

不考虑阻尼影响，在各向同性的弹性介质中，超声波传播控制方程的有限元表达式为：

(1.10)

和分别为刚度矩阵、质量矩阵；、和分别为位移矢量、加速度矢量以及外界作用力矩阵。

在热弹机制下，钢轨表面热载荷力矢量为：

(1.11)

(1.12)

式中，[B]表示个体单元应变矩阵，[D]表示钢轨材料的参数矩阵，{?0}为热应变矢量，{Tref}表示在热分析过程中的参考温度矩阵。

由热膨胀导致的钢轨表面变化的位移信号与时间变化间存在关系曲线，该曲线可以通过对式（1.12）求解得到。求解的过程需要采用Newmark时间积分法进行时间的离散化求解。

综上，钢轨表面的位移和式（1.12）的一次导数为：

(1.13)

(1.14)

式中，Δt为时间步长，γ和α分别为Newmark算法中的数值积分，该值关系着整个算法求解的稳定性与计算的精度。一般而言，γ=0.25和α=0.5时是常见的平均加速度法，此时，积分无条件稳定，t+Δt的运动方程为：

(1.15)

结合式（1.13）可得：

(1.16)

最后，将式1.15代入式1.14中，根据，及计算得到的计算公式：

(1.17)

综上，则可得到热弹机制下，激光束在钢轨表面引起的位移变化。

3热弹耦合方程的有限元模型的建立和验证

在进行有限元模型建立之前，需确定钢轨仿真模型的单元网格划分尺寸及时间步长，这是影响限元数值分析结果精度的关键步骤。

在有限元仿真分析过程中，通常都是根据建模对象的形状自由选择相应的网络形状，其中三角形网格和四边形网格是常用的网格形状。在进行热弹波求解的过程中，一般要求单位网格大小比超声波传播波长的1/4小，这是为了能够