果酱 应用数学学报 1687 - 0042 1110 - 757 x Hindawi出版公司 817680 10.1155 / 2010/817680 817680 研究文章 弹性波在预应力介质中的传播 辛格 英德尔 1 马丹 迪·库马尔 2 古普塔 Manish 3. Jaiani 乔治 1 数学系 J. V. M. G. R. R. (P.G)学院 Charkhi Dadri,哈里亚纳邦127306 印度 jvmgrr.org 2 数学系 纺织与科学技术研究所 Bhiwani 127021 印度 titsbhiwani.org 3. 电子与仪器学系 纺织与科学技术研究所 Bhiwani 127021 印度 titsbhiwani.org 2010 09 12 2010 2010 13 06 2010 29 10 2010 2010 版权所有©2010 Inder Singh等人。 这是一篇在知识共享署名许可下发布的开放存取的文章,它允许在任何媒体上无限制地使用、传播和复制,只要原始作品被适当地引用。

在外力存在的情况下,导出了均匀预应力介质的动力方程的三维解。由一般解得到二维平面波的解,并证明存在两种类型的平面波,即准p波和准sv波。推导了这些波的慢度面和视速度的解析、数值和图形表达式。

1.介绍

事实上,地球是预应力介质,由于许多物理原因,即重力变化、缓慢蠕变过程和温度变化等。因此,由于这些问题在找矿和地震预测中的应用,引起了地震学家的极大兴趣。

为了研究弹性波在无限大预应力固体中的传播,Sidhu和Singh [ 1 2],诺里斯[ 3., Day等[ 4]采用了比奥运动方程的一种简单得多的形式[ 5].早期研究人员考虑的介质是均匀的预应力与增量弹性系数具有正交各向异性由于初始应力的法向分量。

在本文中,我们研究了三维传播问题 P 讨论了均匀预应力介质的SV波和SH波。本文所采用的预应力介质模型比以前的研究人员所采用的更为普遍。

由不同初始应力条件下的一般解得到二维平面波解。推导了慢度曲面图。

2.基本方程

考虑一个均匀的预应力固体。因此,预应力状态由六个分量定义,即: 年代 11 年代 22 年代 33 年代 12 年代 21 年代 31 年代 13 , 年代 23 年代 32 .让所有的应力分量都是 x y z ) .初始应力状态引入各向异性,因此即使对于由两个拉梅常数定义的初始各向同性固体 λ μ ,增量弹性系数的个数 B j ) 总是大于2。让 X Y Z ) 分别为力沿坐标轴的分量,其中 X Y Z 都是常数。预应力固体在外力作用下的动力学方程的一般形式由Biot给出[ 552页): 年代 11 x + 年代 12 y + 年代 31 z + ρ Δ x + ρ ω z Y - ω y Z ) - ρ e X - e x x 年代 11 x - e y y 年代 12 y - e z z 年代 31 z - e y z 年代 31 y + 年代 12 z ) - e z x 年代 11 z + 年代 31 x ) - e x y 年代 12 x + 年代 11 y ) + 年代 11 - 年代 22 ) W z y - 2 年代 12 W z x + 年代 12 W x z + 年代 33 - 年代 11 ) W y z + 2 年代 12 W y x - 年代 31 W x y + 年代 23 W y y - W z z ) ρ 2 u 1 t 2 另外两个方程是通过循环排列得到的 x y z , 1, 2, 3和 X Y Z ρ 为密度 u 1 u 2 u 3. ) 沿轴的位移分量。 Δ X 是否假定增量体力的分量满足方程 Δ X u j X j 0 e e x x + e y y + e z z 年代 j 为预应力分量,假定满足平衡方程。 年代 j j + ρ X 0 与初始应变有关 ϵ j 由胡克定律 年代 j λ ϵ δ j + 2 μ ϵ j 增量的压力 年代 j 应该与增量应变线性相关吗 e j 通过增量弹性系数 B j 年代 11 B 11 e x x + B 12 e y y + B 13 e z z 年代 22 B 21 e x x + B 22 e y y + B 23 e z z 年代 33 B 31 e x x + B 32 e y y + B 33 e z z 年代 23 2 1 e y z 年代 31 2 2 e z x 年代 23 2 3. e x y 我们假设 u 1 u 2 u 3. ) u v w ) X 1 X 2 X 3. ) X Y Z ) e x x u x e y y v y e z z w z e x y 1 2 u y + v x ) e y z 1 2 v z + w y ) e z x 1 2 w x + u z ) W x 1 2 w y - v z ) W y 1 2 u z - w x ) W z 1 2 v x - u y ) 用( 2.2), ( 2.5)和( 2.6) ( 2.1),我们有 B 11 2 u x 2 + 3. - P 1 2 ) 2 u y 2 + 2 + P 3. 2 ) 2 u z 2 - 年代 12 2 v x 2 - 1 2 年代 12 2 v z 2 - 年代 31 2 w x 2 - 1 2 年代 31 2 w y 2 + 年代 12 2 u x y + 年代 23 2 u y z + 年代 31 2 u x z + B 12 + 3. + P 1 2 ) 2 v x y + 1 2 年代 31 2 v y z - 1 2 年代 23 2 v x z + B 13 + 2 - P 3. 2 ) 2 w x z + 1 2 年代 12 2 w y z - 1 2 年代 23 2 w x y - 年代 11 x + ρ X ) u x - 1 2 年代 12 x + 年代 11 y ) u y - 1 2 年代 31 x + 年代 11 z + ρ Z ) u z - 1 2 年代 12 x + 年代 11 y - ρ Y ) v x + 年代 12 y + ρ X ) v y - 1 2 年代 31 y + 年代 12 z ) v z - 1 2 年代 31 x + 年代 11 z - ρ Z ) w x - 1 2 年代 31 y + 年代 12 z ) w y - 年代 31 z + ρ X ) w z ρ 2 u t 2 另外两个方程由( 2.7)按循环排列 x y z , 1, 2, 3, X Y Z u v w

在这里 P 1 年代 11 - 年代 22 P 2 年代 22 - 年代 33 P 3. 年代 33 - 年代 11

3.波的传播

对于沿方向余弦所指定的方向传播的弹性波 l n ) 沿着坐标轴,我们取 u U 经验值 P ) v V 经验值 P ) w W 经验值 P ) U V W ) 振幅因子沿轴和 P 阶段的因素: P k c t - l x + y + n z ) 在哪里 c 为相速度, k 是波数。

把( 3.1)和( 3.2) ( 2.7),我们得到 Ω 1 + Ω 1 + ρ c 2 ) U + 1 + 1 ) V + K 1 + K 1 ) W 0 K 2 + K 2 ) U + Ω 2 + Ω 2 + ρ c 2 ) V + 2 + 2 ) W 0 3. + 3. ) U + K 3. + K 3. ) V + Ω 3. + Ω 3. + ρ c 2 ) W 0 在哪里 Ω 1 - B 11 l 2 + 3. - P 1 2 ) 2 + 2 + P 3. 2 ) n 2 + 年代 12 l + 年代 23 n + 年代 31 n l 1 年代 12 l 2 + 年代 12 2 ) n 2 - B 12 + 3. + P 1 2 ) l - 年代 31 2 ) n + 年代 23 2 ) l n K 1 年代 31 l 2 + 年代 31 2 ) 2 - B 13 + 2 - P 3. 2 ) l n - 年代 12 2 ) n + 年代 23 2 ) l n Ω 1 k 年代 11 x + ρ X ) l + 1 2 年代 12 x + 年代 11 y + ρ Y ) + 1 2 年代 31 x + 年代 11 z + ρ Z ) n 1 k 1 2 年代 12 x + 年代 11 y - ρ Y ) l - 年代 12 y + ρ X ) + 1 2 年代 31 y + 年代 12 z ) n K 1 k 1 2 年代 31 x + 年代 11 z - ρ Z ) l + 1 2 年代 31 y + 年代 12 z ) + 1 2 年代 31 z + ρ X ) n 的显式表达式 Ω 2 Ω 3. ) 2 3. ) , 2 3. ) Ω 2 Ω 3. ) 2 3. ) , 2 3. ) ,由( 3.4) - ( 3.9),通过循环排列 x y z , 1, 2, 3, l n X Y Z ,分别。

设定系数的行列式 U V W ( 3.3= 0,化简后,我们得到三次方程 ρ c 2 ρ 3. c 6 + ρ 2 c 4 一个 + B + D ) + ρ c 2 一个 B + B D + D 一个 - F J - E - H G ) + 一个 B D + E F G + H J - 一个 F J - E D - H G B ) 0 在哪里 一个 Ω 1 + Ω 1 E 1 + 1 H K 1 + K 1 B Ω 2 + Ω 2 F 2 + 2 K 2 + K 2 D Ω 3. + Ω 3. G 3. + 3. J K 3. + K 3.

4.特殊情况

有两个特殊情况可以立即处理。

4.1.三维预应力介质 以下4.4.1。波沿唯一轴的传播

n 1 l 0 在( 3.3) - ( 3.8),在表示 Ω 2 Ω 3. ) 2 3. ) , 2 3. ) Ω 2 Ω 3. ) 2 3. ) , 2 3. ) ,我们得到 Ω 11 - 2 + P 3. 2 ) Ω 21 - 1 - P 2 2 ) Ω 31 - B 33 11 年代 12 2 21 0 31 年代 31 K 11 0 K 21 年代 12 2 K 31 年代 23 Ω 11 2 k 年代 31 x + 年代 11 z + ρ z ) Ω 21 2 k 年代 23 y + 年代 22 z + ρ z ) Ω 31 k 年代 33 x + ρ z ) 11 2 k 年代 31 y + 年代 12 z ) 21 - k 年代 23 z + ρ Y ) 31 k 年代 33 x + 年代 31 z - ρ X ) K 11 k 年代 31 z + ρ X ) K 21 2 k 年代 23 x + 年代 12 z ) K 31 2 k 年代 33 y + 年代 23 z + ρ Y ) 方程( 3.10)以形式 ρ 3. c 6 + ρ 2 c 4 一个 1 + B 1 + D 1 ) + ρ c 2 一个 1 B 1 + B 1 D 1 + D 1 一个 1 - F 1 J 1 - E 1 1 - H 1 G 1 ) + 一个 1 B 1 D 1 + E 1 F 1 G 1 + H 1 1 J 1 - 一个 1 F 1 J 1 - E 1 1 D 1 - H 1 G 1 B 1 ) 0 在哪里 一个 1 Ω 11 + Ω 11 B 1 Ω 21 + Ω 21 D 1 Ω 31 + Ω 31 E 1 11 + 11 F 1 21 G 1 31 + 31 H 1 K 11 1 K 21 + K 21 J 1 K 31 + K 31

4.1.2。在没有身体力量的情况下

X Y Z 0 年代 11 年代 22 年代 33 年代 12 等等都是常数,那么( 4.3.在…的帮助下成为… 4.2)和( 4.4), ρ 3. c 6 + ρ 2 c 4 Ω 11 + Ω 21 + Ω 31 ) + ρ c 2 Ω 11 Ω 21 + Ω 21 Ω 31 + Ω 31 Ω 11 - 11 K 21 ) + Ω 11 Ω 21 Ω 31 - 11 K 21 Ω 31 ) 0

4.1.3。预应力由法向分量定义

在介质中只有正常的预应力分量,即预应力 年代 12 年代 31 年代 23 0 在( 4.1), ( 4.5)以形式 Ω 11 + ρ c 2 ) Ω 21 + ρ c 2 ) Ω 31 + ρ c 2 ) 0 化简后,我们得到三个值 c 2 c 2 2 + P 3. / 2 ρ c 2 1 - P 2 / 2 ρ c 2 B 33 ρ

4.2.二维预应力介质 4.2.1。准备面波的解决方案

在这里,我们考虑平面波在 x y 平面垂直于 z 设在;把 n 0 在( 3.4) - ( 3.9),在表示 Ω 2 Ω 3. ) 等等,我们从( 3.3): Ω 111 + Ω 111 + ρ c 2 ) U + 111 + 111 ) V 0 K 211 + K 211 ) U + Ω 211 + Ω 211 + ρ c 2 ) V 0 在哪里 Ω 111 - B 11 l 2 + 3. - P 1 2 ) 2 + 年代 12 l 111 年代 12 l 2 - B 12 + 3. + P 1 2 ) l Ω 211 - 3. + P 1 2 ) l 2 + B 22 2 + 年代 12 l K 211 年代 12 2 - B 12 + 3. - P 1 2 ) n l Ω 111 k 年代 11 x + ρ X ) l + 1 2 年代 12 x + 年代 11 y + ρ Y ) Ω 211 k 1 2 年代 22 x + 年代 12 y + ρ X ) l + 年代 22 y + ρ Y ) 111 k 1 2 年代 12 x + 年代 11 y - ρ Y ) l - 年代 12 y + ρ X ) K 211 k 年代 12 x + ρ Y ) l + 1 2 年代 22 x + 年代 12 y + ρ X ) 方程( 4.8)有非平凡解时 ρ 2 c 4 + Ω 111 + Ω 211 + Ω 111 + Ω 211 ) ρ c 2 - 111 + 111 ) K 211 + K 211 ) 0 它是二次方程 ρ c 2 ,它有两个值 c 2 对应于准sv波和准sv波 P 波。

在没有身体力量的情况下 年代 11 年代 22 ,等为常量,则从( 4.8)和( 4.9),我们得到 Ω 111 + ρ c 2 ) U + 111 V 0 K 211 U + Ω 211 + ρ c 2 ) V 0 如果 年代 12 年代 23 年代 31 0 ,然后( 4.11)以形式 Ω 1111 + ρ c 2 ) U + 1111 V 0 K 2111 U + Ω 2111 + ρ c 2 ) V 0 在哪里 Ω 1111 - B 11 l 2 + 3. - P 1 2 ) 2 1111 - B 12 + 3. + P 1 2 ) l Ω 2111 - 3. + P 1 2 ) l 2 + B 22 2 K 2111 - B 21 + 3. - P 1 2 ) l 齐次集合( 4.12) U V 什么时候有非平凡解 | Ω 1111 + ρ c 2 ) 1111 K 2111 Ω 2111 + ρ c 2 ) | 0 这个二次方程 ρ c 2 可以解决得到吗 2 ρ c 2 B 11 + 3. + P 1 2 ) l 2 + B 21 + 3. - P 1 2 ) 2 ± B 11 - 3. - P 1 2 ) l 2 + B 22 - 3. + P 1 2 ) 2 2 + 4 B 21 + 3. - P 1 2 ) B 12 + 3. + P 1 2 ) l 2 2 因此,在预应力介质的二维模型中,一般存在两类平面波,即准平面波 P 波和准sv波的相速度分别对应于( 4.15).

4.2.2。平面波在正交异性介质中的传播

考虑一个均匀的预应力弹性实体。材料在有限应变下是各向同性的,在正交各向异性对称下则是各向异性的。初始应力的主方向与弹性对称方向和坐标轴重合。设均匀初始应力状态有主应力 年代 11 年代 22 , 年代 33 .我们进一步假设 年代 22 年代 33 年代 11 年代 22 是常数。的主应力 年代 33 没有明确地加入到运动方程中。然而,它的影响间接地包含在增量弹性系数的值中。我们把 l θ 因为 θ ;( 4.15)可以写成 2 ρ c 2 θ ) B 11 + 3. + P 1 2 ) 2 θ + B 21 + 3. - P 1 2 ) 因为 2 θ ± B 11 - 3. - P 1 2 ) 2 θ + B 22 - 3. + P 1 2 ) 因为 2 θ 2 + 一个 在哪里 4 B 12 + 3. + P 1 / 2 ) B 12 + 3. + P 1 / 2 ) 2 θ 因为 2 θ ,在那里 B 21 - P 1 B 12 C P θ ) C SV θ ) 是…的价值 c 在( 4.16),对应于准-的速度 P 分别为波和准sv波。因此,这些表达式与Sidhu和Singh得到的表达式一致[ 1为准速度 P 波和准sv波。

5.数值计算与讨论

这里我们考虑一个初始应力介质的模型: B 11 λ + 2 μ + P 1 B 12 λ + P 1 B 22 λ + 2 μ 3. μ 使用( 5.1) ( 4.17),我们得到一个无量纲形式的速度方程为 C ̂ P θ ) 1 2 δ + 3. + p ) + δ + 1 + p ) 2 + δ p δ + 1 + p ) 2 θ 1 / 2 C ̂ SV θ ) 1 2 δ + 3. + p ) - δ + 1 + p ) 2 + δ p δ + 1 + p ) 2 θ 1 / 2 在哪里 δ λ μ p P 1 2 μ β 2 μ ρ 准-的视速度 P 波和准sv波可由 5.2), C P 一个 θ ) C ̂ P θ ) θ C SV 一个 θ ) C ̂ SV θ ) θ 无量纲慢度的数值 1 / C ̂ P θ ) 1 / C ̂ SV θ ) ) 是从( 5.2)假设 δ 1 的不同值 p 不同 - 0.8 0 8 对于不同的值 θ 不同 0 ° 90 ° .数据 1 2给出了准的无量纲慢度的变化 P 准sv波的入射角分别为−0.8、−0.6、0.2、0.4、0.6、0.8。

6.结论

研究表明,准 P 波和准sv波受介质中存在的初始应力和传播方向的影响很大。

承认

作者感谢格鲁吉亚第比利斯州立大学的George Jaiani博士提出的宝贵意见和鼓励编写论文。

Sidhu r S。 辛格 美国J。 对“弹性波在自由面上的初应力反射:P和SV运动”的评论 美国声学学会杂志 1983 74 1640 1644 Sidhu r S。 辛格 美国J。 P波和SV波在预应力弹性半空间自由表面的反射 美国声学学会学报 1984 76 2 594 598 2 - s2.0 - 0021639240 Norries a . N。 平面波在预应力弹性介质中的传播 美国声学学会杂志 1983 47 1642 1643 一天 年代。 N。 达塔 一个。 P波和S波的传播和P波在正态初始应力下在介质中的反射 地球物理研究通讯 1985 23 1 28 毕奥 m·A。 增量变形力学 1965 美国纽约 约翰·威利父子公司