双相介质半空间内椭圆夹杂对透射SH波的散射

齐辉; 龚曲; 曾庆友

doi:10.11883/bzycj-2016-0149

双相介质半空间内椭圆夹杂对透射SH波的散射

doi: 10.11883/bzycj-2016-0149

哈尔滨工程大学航天与建筑工程学院, 黑龙江哈尔滨 150001

基金项目:

黑龙江省自然科学基金项目 A201404

详细信息

作者简介:
齐辉(1963-), 男, 博士, 教授, 博导, qihui205@sina.com

中图分类号: O347.3;O343.1
计量
- 文章访问数: 5319
- HTML全文浏览量: 2021
- PDF下载量: 134
- 被引次数: 3
出版历程
- 收稿日期: 2016-05-13
- 修回日期: 2016-07-08
- 刊出日期: 2018-05-25

Scattering of transmission SH-wave by elliptic inclusion in bi-material half-pace

College of Aerospace and Civil Engineering, Harbin Engineering University, Harbin 150001, Heilongjiang, China

摘要

摘要: 为探究双相介质弹性半空间内椭圆弹性夹杂对透射SH波的散射问题，主要采用Green函数法、复变函数法、保角映射法和极坐标移动技术。首先，引入复变量并在复平面上运用保角映射的方法将椭圆边界映射为单位圆边界；然后，将双相介质沿垂直边界剖开分成两个四分之一空间，在剖分面上作用附加力系使SH波在垂直边界上满足位移和应力连续的条件，并构造四分之一空间内点源荷载作用下的Green函数位移场；进而，利用"契合"的思想在垂直边界上建立定解积分方程组，并利用SH波衰减的性质进行有限项截断来求解未知附加力系。最后，通过具体算例得出在不同参数情况下椭圆夹杂周边动应力集中因子分布情况。结果得知，SH波的入射角度和频率以及介质的性质对夹杂周边动应力集中分布有一定影响。
- 散射 /
- 透射SH波 /
- 双相介质 /
- 椭圆夹杂 /
- 动应力集中因子
Abstract: In this paper we studied the scattering of the transmission SH-wave form by elliptic elastic inclusions in a bi-material half-space using mainly the Green function, the complex function, the conformal mapping and the polar coordinate. Firstly, by introducing the complex variables and using the conformal mapping, we mapped the elliptic boundary into a unit circle. Next, by dividing the double-phase medium along the vertical boundaries into two quarter-spaces and adding a force system on the subdivision surface to make the SH wave satisfy the displacement and stress continuity conditions in the vertical boundary, we constructed a quarter-space Green function of displacement field. Then, we established an additional force system of integro differential equations in the vertical interface using the "conjunction" technique, and solved the unknown additional force system using the effective truncation. Finally, we presented and discussed the distribution of the dynamic stress concentration factor around the elliptic inclusion under the conditions of different parameters. The results show that the SH-wave incident angle and frequency, and the media both have influences on the distribution of the dynamic stress concentration factor.
- scattering /
- transmission SH-wave /
- bi-material /
- elliptic inclusion /
- dynamic stress concentration factor

HTML全文

弹性波散射理论一直以来都是弹性动力学中重要研究课题之一，巧妙地运用了一些数学物理方法解析地求解了一些复杂的波动问题，对地震工程、岩土工程及地下结构工程等相关技术的研究与应用有着重要价值。对弹性波在全空间中传播时遇缺陷发生散射的研究已日趋完善，相关的研究已扩展到半空间、四分之一空间等更复杂的情况，更多的界面模型被涉及^[1-8]。双相介质半空间中的缺陷对SH波的散射问题是在近几年才备受科研人员的重视，而且他们基本上是讨论入射SH波与缺陷在同一个介质中的情况^[2-6]，对于缺陷与入射波处于不同介质中的研究还非常少。本文中采用Green函数法、复变函数法、保角映射法、“镜像”法、极坐标移动技术以及“契合”的思想解析地求解双相介质弹性半空间内椭圆弹性夹杂对透射SH波的散射问题。并通过具体的算例得出在不同的入射角、SH波频率和介质性质的情况下椭圆夹杂周边环向动应力集中分布情况，以期获得一些具体的理论结果为相关科研及工程实际应用提供参考。

1. 问题模型

如图 1所示，由介质Ⅰ和介质Ⅱ组成的双相介质半空间内有一个椭圆夹杂，椭圆夹杂为介质Ⅲ，这3种介质均为连续、均匀、各向同性的弹性介质。椭圆夹杂长半轴和短半轴长度分别为a和b，夹杂中心到垂直边界B_V的距离为h，到水平边界B_H距离为d，建立如图 1所示直角坐标系xOy和x″O″y″。SH波从介质Ⅱ中入射，遇垂直边界B_V产生透射SH波进入介质Ⅰ中。主要研究在不同参数条件下介质Ⅰ中椭圆夹杂周边在透射SH波作用下产生的环向动应力集中情况，并对结果进行分析，突出反映透射SH波的危害性。

图 1 理论模型

Figure 1. Theoretical model

下载: 全尺寸图片幻灯片

2. 控制方程

二维平面SH波位移函数W(x, y)与时间的依赖关系为e^-iωt, 满足控制方程：

$\frac{{{\partial ^2}W}}{{\partial {x^2}}} + \frac{{{\partial ^2}W}}{{\partial {y^2}}} + {k^2}W = 0$

(1)

式中: $k = \frac{\omega }{{{c_{\rm{s}}}}}$ ， ${c_{\rm{s}}} = \sqrt {\frac{\mu }{\rho }}$ ，k为波数，ω为位移函数圆频率，c_s为介质的剪切波速，μ为介质的剪切模量，ρ为介质密度。引入复变量z=x+iy，z=x-iy，位移函数W(z, z)在复平面(z, z)上控制方程的表达形式为：

$\frac{{{\partial ^2}W}}{{\partial z\partial \bar z}} + \frac{1}{4}{k^2}W = 0$

(2)

引入保角映射函数：

$z = \omega \left( \eta \right) = R\left( {\eta + \frac{m}{\eta }} \right)$

(3)

式中: $\eta = \mathit{R}{\mathit{e}^{{\rm{i}}\theta }}, R = \frac{{a + b}}{2}, m = \frac{{a-b}}{{a + b}}$ ，a、b分别为椭圆长半轴和短半轴长。

控制方程在映射平面(η, η)上可以表示为：

$\frac{1}{{\omega '\left( \eta \right)\overline {\omega '\left( \eta \right)} }}\frac{{{\partial ^2}W}}{{\partial \eta \partial \bar \eta }} + \frac{1}{4}{k^2}W = 0$

(4)

与式(4)相应的应力表达式为：

${\tau _{rz}} = \frac{\mu }{{R\left| {\omega '\left( \eta \right)} \right|}}\left( {\eta \frac{{\partial W}}{{\partial \eta }} + \bar \eta \frac{{\partial W}}{{\partial \bar \eta }}} \right),\;\;\;\;\;\;\;\;{\tau _{\theta z}} = \frac{{{\rm{i}}\mu }}{{R\left| {\omega '\left( \eta \right)} \right|}}\left( {\eta \frac{{\partial W}}{{\partial \eta }} - \bar \eta \frac{{\partial W}}{{\partial \bar \eta }}} \right)$

(5)

3. Green函数

本文中在求解实际问题之前先构造问题的Green函数，其中Green函数Ⅰ具体为含椭圆弹性夹杂的四分之一空间在垂直边界上任一位置处的出平面点源荷载作用下产生的位移场，四分之一空间为介质Ⅰ，椭圆夹杂为介质Ⅲ，如图 2所示。采用“虚设点源”法构造Green函数入射波位移场表达式如下：

${G^{\left( {{\rm{in}}} \right)}}\left( {\eta ,\bar \eta } \right) = \frac{{\rm{i}}}{{2{\mu _1}}}H_0^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - \omega \left( {{\eta _0}} \right)} \right|} \right) + \frac{{\rm{i}}}{{2{\mu _1}}}H_0^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - \overline {\omega \left( {{\eta _0}} \right)} + 2{\rm{i}}d} \right|} \right)$

(6)

图 2 虚设点源模型

Figure 2. Model of dummy point source loads

下载: 全尺寸图片幻灯片

式中:i为虚数单位，H₀⁽¹⁾为零阶第一类Hankel函数，ω(η₀)为点源荷载的位置矢量。构造四分之一空间内椭圆弹性夹杂产生的散射波位移场表达式^[4]为：

${G^{\left( {\rm{s}} \right)}}\left( {\eta ,\bar \eta } \right) = \sum\limits_{n = - \infty }^\infty {{A_n}\sum\limits_{j = 1}^4 {S_n^{\left( j \right)}} }$

(7)

式中：

$S_n^{\left( 1 \right)} = H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right)} \right|} \right){\left[ {\frac{{\omega \left( \eta \right)}}{{\left| {\omega \left( \eta \right)} \right|}}} \right]^n}，$

$S_n^{\left( 2 \right)} = H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - 2h} \right|} \right){\left[ {\frac{{\omega \left( \eta \right) - 2h}}{{\left| {\omega \left( \eta \right) - 2h} \right|}}} \right]^{ - n}}，$

$S_n^{\left( 3 \right)} = {\left( { - 1} \right)^n}H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - 2{\rm{i}}d} \right|} \right){\left[ {\frac{{\omega \left( \eta \right) - 2{\rm{i}}d}}{{\left| {\omega \left( \eta \right) - 2{\rm{i}}d} \right|}}} \right]^{ - n}}，$

$S_n^{\left( 4 \right)} = {\left( { - 1} \right)^n}H_n^{\left( 1 \right)}\left( {{k_1}\left| {\omega \left( \eta \right) - 2{\rm{i}}d - 2h} \right|} \right){\left[ {\frac{{\omega \left( \eta \right) - 2{\rm{i}}d - 2h}}{{\left| {\omega \left( \eta \right) - 2{\rm{i}}d - 2h} \right|}}} \right]^n}。$

椭圆弹性夹杂内部驻波表达式如下：

${G^{\left( {\rm{t}} \right)}}\left( {\eta ,\bar \eta } \right) = \sum\limits_{{\rm{n}} = - \infty }^\infty {{{\rm{B}}_n}{J_n}\left( {{k_3}\left| {\omega \left( \eta \right)} \right|} \right){{\left[ {\frac{{\omega \left( \eta \right)}}{{\left| {\omega \left( \eta \right)} \right|}}} \right]}^n}}$

(8)

然后依据椭圆夹杂边界上应力和位移连续条件可建立如下方程组：

$\left\{ \begin{array}{l} {G^{\left( {{\rm{in}}} \right)}}\left( {\eta ,\bar \eta } \right) + {G^{\left( {\rm{s}} \right)}}\left( {\eta ,\bar \eta } \right) = {G^{\left( {\rm{t}} \right)}}\left( {\eta ,\bar \eta } \right)\\ \tau _{rz}^{\left( {{\rm{in}}} \right)}\left( {\eta ,\bar \eta } \right) + \tau _{rz}^{\left( {\rm{s}} \right)}\left( {\eta ,\bar \eta } \right) = \tau _{rz}^{\left( {\rm{t}} \right)}\left( {\eta ,\bar \eta } \right) \end{array} \right.$

(9)

在方程两边同时乘以e^-imθ，然后在(-π, π)区间内积分并进行有限项截断求出系数A_n和B_n，本文中取n=8。Green函数Ⅰ为:

${G_1}\left( {\eta ,\bar \eta } \right) = G_{\rm{I}}^{\left( {{\rm{in}}} \right)}\left( {\eta ,\bar \eta } \right) + G_{\rm{I}}^{\left( {\rm{s}} \right)}\left( {\eta ,\bar \eta } \right)$

(10)

Green函数Ⅱ为四分之一空间中无椭圆夹杂时的位移场，即:

${G_2}\left( {\eta ,\bar \eta } \right) = G_{{\rm{II}}}^{\left( {{\rm{in}}} \right)}\left( {\eta ,\bar \eta } \right)$

(11)

4. 入射波、反射波、透射波和散射波

采用“镜像”法将半空间双相介质问题转化为全空间双相介质问题，入射波、反射波和透射波可分别表示为：

$\begin{array}{l} {W^{\left( {{\rm{in}}} \right)}} = {W_0}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _0}}}} \right]} \right\} + \\ {W_0}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h - 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h + 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _0}}}} \right]} \right\} \end{array}$

(12)

$\begin{array}{l} {W^{\left( {{\rm{re}}} \right)}} = {W_2}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h} \right){{\rm{e}}^{{\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h} \right){{\rm{e}}^{{\rm{i}}{\alpha _0}}}} \right]} \right\} + \\ {W_2}\exp \left\{ {\frac{{{\rm{i}}{k_2}}}{2}\left[ {\left( {\omega \left( \eta \right) + h - 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _0}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h + 2{\rm{i}}d} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _0}}}} \right]} \right\} \end{array}$

(13)

$\begin{array}{l} {W^{\left( {{\rm{tr}}} \right)}} = {W_4}\exp \left\{ {\frac{{{\rm{i}}{k_1}}}{2}\left[ {\left( {\omega \left( \eta \right) + h} \right){{\rm{e}}^{ - {\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _4}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h} \right){{\rm{e}}^{ - {\rm{i}}{\alpha _4}}}} \right]} \right\} + \\ {W_4}\exp \left\{ {\frac{{{\rm{i}}{k_1}}}{2}\left[ {\left( {\omega \left( \eta \right) + h - 2{\rm{i}}d} \right){{\rm{e}}^{{\rm{i}}\left( {{\rm{ \mathsf{ π} }} - {\alpha _4}} \right)}} - \left( {\overline {\omega \left( \eta \right)} + h + 2{\rm{i}}d} \right){{\rm{e}}^{{\rm{i}}{\alpha _4}}}} \right]} \right\} \end{array}$

(14)

运用Snell定律可以得到入射波、反射波和透射波位移幅值之间的关系式如下：

${W_2} = {W_0}\frac{{\cos {\alpha _0} - \frac{{{c_2}}}{{{c_1}}}\frac{{{\mu _1}}}{{{\mu _2}}}\cos {\alpha _4}}}{{\cos {\alpha _0} + \frac{{{c_2}}}{{{c_1}}}\frac{{{\mu _1}}}{{{\mu _2}}}\cos {\alpha _4}}},\;\;\;\;\;\;\;\;\;\;{W_4} = {W_0}\frac{{2\cos {\alpha _0}}}{{\cos {\alpha _0} + \frac{{{c_2}}}{{{c_1}}}\frac{{{\mu _1}}}{{{\mu _2}}}\cos {\alpha _4}}}$

(15)

由于入射波的作用，同样在介质Ⅰ中会产生散射波，在椭圆夹杂内部会产生驻波，其过程及相关系数的求解参考构造Green函数的过程，在此不再叙述。

5. 定解积分方程组

如图 3所示，本文中将双相介质沿垂直界面剖开，在剖分面上作用连续分布的入平面和出平面点源荷载，分别表示为f₁(r″₀, θ″₀)和f₂(r″₀, θ″₀)，在垂直边界上满足应力连续条件，即:

$\tau _{\theta ''z''}^{\left( {\rm{I}} \right)}\sin {{\theta ''}_0} + {f_1}\left( {{{r''}_0} + {{\theta ''}_0}} \right) = \tau _{\theta ''z''}^{\left( {{\rm{II}}} \right)}\sin {{\theta ''}_0} + {f_2}\left( {{{r''}_0} + {{\theta ''}_0}} \right)$

(16)

图 3 剖分面模型

Figure 3. Model of cutaway interface

下载: 全尺寸图片幻灯片

介质Ⅰ中存在透射波和散射波，即:

${W^{\left( {\rm{I}} \right)}} = {W^{\left( {{\rm{tr}}} \right)}} + {W^{\left( {\rm{s}} \right)}},\;\;\;\;\;\;\;\tau _{\theta z}^{\left( {\rm{I}} \right)} = \tau _{\theta z}^{\left( {{\rm{tr}}} \right)}$

(17)

介质Ⅱ中存在入射波和反射波，即:

${W^{\left( {{\rm{II}}} \right)}} = {W^{\left( {{\rm{in}}} \right)}} + {W^{\left( {{\rm{re}}} \right)}},\;\;\;\;\;\;\;\tau _{\theta z}^{\left( {{\rm{II}}} \right)} = \tau _{\theta z}^{\left( {{\rm{in}}} \right)} + \tau _{\theta z}^{\left( {{\rm{re}}} \right)}$

(18)

然后利用“契合”的思想将介质Ⅰ和介质Ⅱ“契合”在一起，在垂直界面处τ_θz⁽ⁱⁿ⁾+τ_θz^(re)=τ_θz^(tr)，于是可得f₁(r″₀, θ″₀)=f₂(r″₀, θ″₀)。根据前面构造的Green函数，通过积分的方法可以得到外力系f₁(r″₀, θ″₀)与f₂(r″₀, θ″₀)相应的位移表达式，即:

${W^{\left( {{f_1}} \right)}}\left( {r'',\theta ''} \right) = \int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_1}\left( {r'',\theta '',{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}} + \\ \int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_1}\left( {r'',\theta '',{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}}$

(19)

${W^{\left( {{f_2}} \right)}}\left( {r'',\theta ''} \right) = - \int_0^h {{f_2}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_2}\left( {r'',\theta '',{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}} - \\ \int_0^\infty {{f_2}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){G_2}\left( {r'',\theta '',{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right){\rm{d}}{{r''}_0}}$

(20)

利用垂直边界处的关系：W^(Ⅰ)+W^(f₁)=W^(Ⅱ)+W^(f₂)，W⁽ⁱⁿ⁾+W^(re)=W^(tr)可得到如下定解积分方程组：

$\begin{array}{*{20}{c}} {\int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} + }\\ {\int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} = \\ - {W^{\left( {\rm{s}} \right)}}\left| {_{{{\theta ''}_0} = \frac{{\rm{ \mathsf{ π} }}}{2}}} \right.} \end{array}$

(21)

$\begin{array}{*{20}{c}} {\int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} + }\\ {\int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)\left[ {{G_1}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right) + {G_2}\left( {r'',\frac{{\rm{ \mathsf{ π} }}}{2},{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)} \right]{\rm{d}}{{r''}_0}} = \\ - {W^{\left( {\rm{s}} \right)}}\left| {_{{{\theta ''}_0} = - \frac{{\rm{ \mathsf{ π} }}}{2}}} \right.} \end{array}$

(22)

式中:G₁与G₂分别为之前构造的在介质Ⅰ和介质Ⅱ中的Green函数，利用散射波的衰减性并进行有限项截断求得未知附加力系。

6. 动应力集中因子

本文中主要是运用弹性波的散射理论来研究弹性波的绕射与动应力集中问题。环向动应力集中因子可写成如下形式：

$\tau _{\theta z}^ * = \left| {\frac{\tau }{{{\tau _0}}}} \right|$

(23)

式中:τ₀=μ₂k₂W₄为半空间透射SH波应力的最大幅值，τ为椭圆夹杂周边环向总应力。

在椭圆夹杂周边环向应力表达式：

${\tau _{\theta z}} = \tau _{\theta z}^{\left( {{\rm{tr}}} \right)} + \tau _{\theta z}^{\left( {\rm{s}} \right)} + \int_0^h {{f_1}\left( {{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)\frac{{{\mu _1}}}{r}\frac{{\partial {G_1}\left( {r,\theta ,{{r''}_0},\frac{{\rm{ \mathsf{ π} }}}{2}} \right)}}{{\partial \theta }}{\rm{d}}{{r''}_0}} + \\ \int_0^\infty {{f_1}\left( {{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)\frac{{{\mu _1}}}{r}\frac{{\partial {G_1}\left( {r,\theta ,{{r''}_0}, - \frac{{\rm{ \mathsf{ π} }}}{2}} \right)}}{{\partial \theta }}{\rm{d}}{{r''}_0}}$

(24)

7. 算例及结果分析

通过引入量纲一参数并对其赋值得到了一些具体算例的结果，给出了椭圆夹杂周边环向动应力集中因子在不同参数情况下的分布情况。取量纲一参数 ${\mu _{21}} = \frac{{{\mu _2}}}{{{\mu _1}}},{\mu _{31}} = \frac{{\mu 3}}{{{\mu _1}}},{k_{21}} = \frac{{{k_2}}}{{{k_1}}},{k_{31}} = \frac{{k3}}{{{k_1}}}$ ；透射波量纲一波数为k₁b；椭圆夹杂位置坐标量纲一参数为 $\frac{d}{b}$ 和 $\frac{h}{b}$ 。

图 4所示为椭圆夹杂周边环向动应力集中因子τ_θz^*随不同入射角α₀的分布情况。此处取k₂₁=1.0，k₃₁=1.0，μ₂₁=2.0，μ₃₁=0.5， $\frac{a}{b}$ =1.25, $\frac{h}{b} = \frac{d}{b} = 12.0$ 。可以看出，当SH波水平向入射时，椭圆夹杂周边环向动应力集中因子明显大于其他入射角度时的相应值，且当入射角α₀=45°时的相应值最小，α₀=30°和α₀=60°时的相应值较为接近且大小处于居中位置。当透射波波数k₁b=0.1时，即在“准静态”情况下，SH波水平入射时椭圆夹杂周边环向动应力集中因子最大值是入射角α₀=45°时相应最大值的1.7倍，比其他透射波波数情况下相应值大。由此可以得出，当SH波水平入射产生透射波时，椭圆夹杂周边环向动应力集中情况最严重，且不同的透射波波数对动应力集中的分布有一定的影响。

图 4 SH波以不同角度入射时动应力集中因子的分布情况

Figure 4. Distribution of dynamic stress concentration factor with different incident angles disturbed by SH wave

下载: 全尺寸图片幻灯片

图 5所示为SH波水平入射产生透射波时，椭圆夹杂周边环向动应力集中因子τ_θz^*随不同波数比k₂₁的分布情况。其中量纲一参数k₃₁=1.2，μ₃₁=1.0，μ₂₁=0.2， $\frac{a}{b}$ =1.25， $\frac{h}{b}$ =2.0， $\frac{d}{b}$ =12.0。可以看出，当透射波波数k₁b=0.1时，即在“准静态”情况下，波数比k₂₁对夹杂周边环向动应力集中因子的分布几乎没有影响，在其他透射波波数情况下，夹杂周边环向动应力集中因子最大值随k₂₁增大而略微减小，整体上变化并不大。当透射波波数发生变化时，夹杂周边环向动应力集中因子最大值的位置发生了明显变化。由此可以得出，垂直界面右侧介质的性质对椭圆夹杂周边动应力集中因子的分布影响较小，但是透射波频率对椭圆夹杂周边动应力集中因子最大值的位置影响较大。

图 5 SH波水平向入射时动应力集中因子随k₂₁的分布情况

Figure 5. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with k₂₁ disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

图 6所示为SH波水平向入射产生透射波时，椭圆弹性夹杂周边环向动应力集中因子τ_θz^*随波数比k₃₁变化的分布情况。其中量纲一参数k₂₁=1.0，μ₂₁=2.0，μ₃₁=1.0， $\frac{a}{b}$ =1.25， $\frac{h}{b} = \frac{d}{b} = 12.0$ 。可以看出，当透射波波数k₁b=0.1时，即在“准静态”情况下，波数比k₃₁对夹杂周边环向动应力集中因子的分布影响较小。在中频和高频透射波情况下，夹杂周边环向动应力集中因子变化非常明显。当k₁b=0.5，k₃₁=4.0时，环向动应力集中因子τ_θz^*达到了极值|τ_θz^*|=7.63，比k₁b=0.5，k₃₁=0.5时的相应极值增大了3.06倍。由此可得知，在中频和高频透射波情况下，波数比k₃₁对椭圆夹杂周边环向动应力集中因子影响非常大，且极值的位置会随透射波频率的变化而发生改变。

图 6 SH波水平向入射时动应力集中因子随k₃₁的分布情况

Figure 6. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with k₃₁ disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

图 7所示为SH波水平入射产生透射波时，椭圆弹性夹杂周边在θ=90°处环向动应力集中因子τ_θz^*随 $\frac{h}{b}$ 变化的分布情况。其中量纲一参k₂₁=1.0，k₃₁=1.0，μ₂₁=2.0，μ₃₁=2.0， $\frac{a}{b}$ =1.25， $\frac{d}{b}$ =12.0。可以看出，椭圆夹杂周边在θ=90°处环向动应力集中因子τ_θz^*随 $\frac{h}{b}$ 的增大而呈周期性变化，透射波频率越大，周期越小，且振幅呈衰减趋势，衰减到一定程度后逐渐趋于稳定。当透射波波数k₁b=0.1时，即在“准静态”情况下接近为一条直线。可以得知，椭圆弹性夹杂距离垂直边界一定距离后，该距离对夹杂周边环向动应力集中因子的分布影响可以忽略，且这个距离值随透射波频率增大而增大。

图 7 SH波水平向入射时动应力集中因子随h/b的分布情况

Figure 7. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with h/b disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

8. 结论

SH波入射的角度、介质的性质、透射波的频率以及椭圆夹杂与垂直边界的距离对椭圆夹杂周边环向动应力集中因子的分布均有不同程度的影响。当SH波水平入射产生透射波时，椭圆夹杂周边环向动应力集中程度较大。在一定范围内，垂直界面右侧的介质性质对动应力集中因子的极值影响较小，但在不同透射波频率下该极值的位置会发生明显变化。在中频和高频透射波情况下，椭圆夹杂的性质对其周边动应力集中的分布影响较大，且在一定条件下动应力集中程度会非常严重。当椭圆夹杂距离垂直边界一定距离后，该距离的影响可以忽略不计。总之，在实际工程中应该重视双相介质中透射波对结构可能造成的不利影响。

图 1 理论模型

Figure 1. Theoretical model

下载: 全尺寸图片幻灯片

图 2 虚设点源模型

Figure 2. Model of dummy point source loads

下载: 全尺寸图片幻灯片

图 3 剖分面模型

Figure 3. Model of cutaway interface

下载: 全尺寸图片幻灯片

图 4 SH波以不同角度入射时动应力集中因子的分布情况

Figure 4. Distribution of dynamic stress concentration factor with different incident angles disturbed by SH wave

下载: 全尺寸图片幻灯片

图 5 SH波水平向入射时动应力集中因子随k₂₁的分布情况

Figure 5. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with k₂₁ disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

图 6 SH波水平向入射时动应力集中因子随k₃₁的分布情况

Figure 6. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with k₃₁ disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

图 7 SH波水平向入射时动应力集中因子随h/b的分布情况

Figure 7. Distribution of dynamic stress concentration factor around the elliptic inclusion edge with h/b disturbed by SH wave horizontally

下载: 全尺寸图片幻灯片

参考文献(8)

[1]	齐辉, 蔡立明, 罗广龙, 等.双相介质半空间中圆孔对透射SH波的稳态分析[J].爆炸与冲击, 2015, 35(4):591-598. doi: 10.11883/1001-1455(2015)04-0591-08 QI Hui, CAI Liming, LUO Guanglong. Steady state analysis for circular cavity impacted by transmitted SH wave in a bi-material half space[J]. Explosion and Shock Waves, 2015, 35(4):591-598. doi: 10.11883/1001-1455(2015)04-0591-08
[2]	齐辉, 折勇, 李宏亮, 等.SH波入射时垂直半空间中双相介质界面附近圆孔的动力分析[J].爆炸与冲击, 2009, 29(1):73-79. doi: 10.11883/1001-1455(2009)01-0073-07 QI Hui, SHI Yong, LI Hongliang, et al.Dynamic analysis for scattering of SH-wave by circular cavities near bi-material interfaces in a vertical half-space[J]. Explosion and Shock Waves, 2009, 29(1):73-79. doi: 10.11883/1001-1455(2009)01-0073-07
[3]	刘殿魁, 林宏.SH波对双相介质界面附近圆形孔洞的散射[J].固体力学学报, 2003, 24(2):197-204. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gtlxxb200302010 LIU Diankui, LIN Hong. Scattering of SH-waves by circular cavities near bimaterial interface[J]. Acta Mechanica Solida Sinica, 2003, 24(2):197-204. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gtlxxb200302010
[4]	齐辉, 丁晓浩, 赵元博.双相介质弹性半空间垂直界面附近椭圆形夹杂对SH波的散射[J].振动与冲击, 2015, 34(24):76-81. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_zdycj201524013 QI Hui, DING Xiaohao, ZHAO Yuanbo. Scattering of SH-wave by an elliptic inclusion near vertical interface in a bi-material half-space[J]. Journal of Vibration and Shock, 2015, 34(24):76-81. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_zdycj201524013
[5]	齐辉, 杨杰.SH波入射双相介质半空间浅埋任意位置圆形夹杂的动力分析[J].工程力学, 2012, 29(7):320-327. doi: 10.6052/j.issn.1000-4750.2010.09.0680 QI Hui, YANG Jie. Dynamic analysis for shallowly buried circular inclusions of arbitrary positions impacted by SH-wave in bi-material half space[J]. Engineering Mechanics, 2012, 29(7):320-327. doi: 10.6052/j.issn.1000-4750.2010.09.0680
[6]	杨杰. 半空间垂直界面裂纹及圆夹杂对SH波的散射[D]. 哈尔滨: 哈尔滨工程大学, 2012: 22-34. http://cdmd.cnki.com.cn/Article/CDMD-10217-1012518460.htm YANG Jie. Scattering of SH-wave by vertical interface crack and circular inclusions in half space[D]. Harbin: Harbin Engineering University, 2012: 22-34. http://cdmd.cnki.com.cn/Article/CDMD-10217-1012518460.htm
[7]	杨在林, 许华南, 黑宝平.SH波上方垂直入射时界面附近椭圆夹杂与裂纹的动态响应[J].岩土力学, 2013, 34(8):2378-2384. http://www.cnki.com.cn/Article/CJFDTotal-YTLX201308054.htm YANG Zailin, XU Huanan, HEI Baoping. Dynamic response of elliptical inclusion and crack neat interface under vertically incident SH-wave from above[J]. Rock and Soil Mechanics, 2013, 34(8):2378-2384. http://www.cnki.com.cn/Article/CJFDTotal-YTLX201308054.htm
[8]	杨在林, 许华南, 黑宝平.半空间椭圆夹杂与裂纹对SH波的散射[J].振动与冲击, 2013, 32(11):56-61. doi: 10.3969/j.issn.1000-3835.2013.11.012 YANG Zailin, XU Huanan, HEI Baoping. Interaction of elliptical and crack under incident SH-wave in a half-space[J]. Journal of Vibration and Shock, 2013, 32(11):56-61. doi: 10.3969/j.issn.1000-3835.2013.11.012

施引文献

期刊类型引用(2)

1.	屈恩相，张景颢，王伟业，胡兴森，王浩宁，徐静怡，李聪. 镜像法求解SH波散射问题的应用综述. 安徽建筑. 2024(07): 115-117 . 百度学术
2.	屈恩相，齐辉，郭晶，杨杰，郑易. “分区”与“契合”思想求解半空间含凸起地形对SH波散射问题的研究进展. 黑龙江工业学院学报(综合版). 2022(04): 74-82 . 百度学术

其他类型引用(1)

资源附件(0)

访问统计

图(7)

计量

文章访问数: 5319
HTML全文浏览量: 2021
PDF下载量: 134
被引次数: 3

1. 问题模型
2. 控制方程
3. Green函数
4. 入射波、反射波、透射波和散射波
5. 定解积分方程组
6. 动应力集中因子
7. 算例及结果分析
8. 结论

双相介质半空间内椭圆夹杂对透射SH波的散射

doi: 10.11883/bzycj-2016-0149

作者简介: 齐辉(1963-), 男, 博士, 教授, 博导, qihui205@sina.com

计量

出版历程