金属材料表面的声弹性理论研究.docxVIP

金属材料表面的声弹性理论研究 0 声弹性法检测应力系统实验研究 声音弹性法的强度检测具有快速、损伤和材料适应性强的特点。它可以检测物体的内部和外部压力，这对设计应用具有很大的潜在价值。该方法通过测量超声波在被检测对象中的传播速度的变化来确定应力值。由于应力引起的超声波速度变化很小,对声波传播时间的分辨率要求达到纳秒级,再加上对传感器与表面耦合因素的影响研究不够,目前声弹性法检测应力总体上还处在实验研究阶段。 利用表面波测量物体表面的工作应力和残余应力方面,近年来的研究开始活跃起来,如Husson等人研究了在一定声程内表面波速度的变化和应力之间的关系。Hirao曾研究了各向同性弹性体中的形变与表面波传播速度的一般关系,但不便于直接进行应力检测。M.Duquennoy从弹性波理论出发研究建立了表面波传播相位变化与应力的线性关系,但实验研究还不够。 国内相关的实验研究目前还处于起步阶段,由于仪器精度不够,目前只对声弹性效果较显著的航空有机玻璃的表面应力进行实验研究,而对工程结构中普遍使用的金属材料的声弹性实验研究还没有文献报道。本文采用5GHz采样频率的采集系统对金属材料的表面波应力检测做了详细的实验研究,简化了声弹性公式,得到了相对传播时间差和应力水平之间的表达式,并分析了实际检测中各种因素对实验结果的影响。 1 初始位移v0 小幅度表面波在受应力作用的各向同性弹性体表面上传播时,表面波的速度与无应力作用时有所不同。具有质点速度v的表面波,在沿着应力场为T的弹性体表面传播了一定的距离后,相位的变化为δ?(ω)。 δ?(ω)=-ω4Ρ∫VGdV(1)δ?(ω)=?ω4P∫VGdV(1) 在(1)式中ω表示圆频率;V表示表面波在弹性体内的传播区域。G为材料的弹性二阶(λ,μ)、三阶(l,m,n)弹性系数和表面波传播区域的位移梯度组成的函数。P表示一段时间内表面波垂直于传播方向上一个波长的宽度内的平均能量流,P可以表达为表面波的速度以及弹性体密度等组成的函数。 Ρ=ωρ0V02(1/V02+ks22ks-2Κ2/V20+Κ4ks2ks+k1+Κ22/V02+Κ42ks22k1)(2)P=ωρ0V02(1/V02+ks22ks?2K2/V20+K4ks2ks+k1+K22/V02+K42ks22k1)(2) 其中: ks=(1V02-1VS2)1/2k1=(1V02-1Vl2)1/2Κ2=2ksk1[1/V02+ks2]Κ4=2[1+V02ks2]ks=(1V02?1VS2)1/2k1=(1V02?1Vl2)1/2K2=2ksk1[1/V02+ks2]K4=2[1+V02ks2] 上式中,V0,VS和Vl分别为表面波、横波和纵波在无应力作用的弹性体中的传播速度;ρ0表示无应力状态下传播介质的密度。 假设表面波沿着图1所示的a3方向传播,a3和a2组成一个纵向的平面,且为表面波的偏振平面。 表面波只在弹性体的表面一个波长左右的深度传播,沿着a3方向声程为L0,且表面应力可以近似视为常数。假设弹性体是各向同性的,表面只有沿着a3轴方向应力σ33作用,应力、应变关系满足胡克定律。则初始位移b3、b1的梯度可以表达如下: ?b3?a3=λ+μμ(3λ+2μ)σ33(3)?b1?a1=?b2?a2=-λ2μ(3λ+2μ)σ33(4)?b3?a3=λ+μμ(3λ+2μ)σ33(3)?b1?a1=?b2?a2=?λ2μ(3λ+2μ)σ33(4) 则(1)式的积分形式可以简化为: δ?33(ω)=-L0ω4Ρ∫∞0αi//Fi(α2,ω)σ33da2(5)δ?33(ω)=?L0ω4P∫∞0αi//Fi(α2,ω)σ33da2(5) 下标i∈{1,…,5}表示求和。δ?αβ(ω)表示沿着αβ方向的单向应力对沿着aα方向传播的表面波的相位的改变量。δ?α(β+γ)则为αβ、αγ方向的平面应力对沿着aα方向传播的表面波的相位改变量。 Fi(a2,ω)=ω2{fi1e-2ωksa2+fi2e-2ωk1a2+fi3e-ω(k1+k2)a2}(6)Fi(a2,ω)=ω2{fi1e?2ωksa2+fi2e?2ωk1a2+fi3e?ω(k1+k2)a2}(6) fij为常数,i∈{1,...,5},j∈{1,2,3}, 又 δ?33(ω)2πω=Δt33(7)δ?33(ω)2πω=Δt33(7) Δt33=t-t0为在声程L0上的表面波相对传播时间差。 由(5)式可得: σ33=8πωΡΔt33/[-L0ω∫∞0αi//Fi(a2,ω)da2](8)σ33=8πωPΔt33/[?L0ω∫∞0αi//Fi(a2,ω)da2](8) 同理对于a1方向作用的单向应力而表面波沿着a3方向传播时: σ11=8πωΡΔt31/[-L0ω∫∞0αi⊥Fi(a2,ω)da](9)σ1