Điểm osculation của sóng rayleigh trong một số mô hình

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - DOÃN THU HƢƠNG ĐIỂM OSCULATION CỦA SÓNG RAYLEIGH TRONG MỘT SỐ MÔ HÌNH LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội - 2015 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - DOÃN THU HƢƠNG ĐIỂM OSCULATION CỦA SÓNG RAYLEIGH TRONG MỘT SỐ MÔ HÌNH Mã số: 60 44 01 07 Chuyên ngành: Cơ học vật thể rắn LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC TS TRẦN THANH TUẤN Hà Nội - 2015 LỜI CẢM ƠN Lời đầu tiên, em xin gửi lời cảm ơn chân thành tới thầy giáo hướng dẫn TS Trần Thanh Tuấn, người giao đề tài quan tâm, tận tình hướng dẫn em suốt trình thực luận văn Em xin bày tỏ lòng biết ơn chân thành tới nhóm Seminar môn Cơ học PGS TS Phạm Chí Vĩnh chủ trì, toàn thể thầy cô giáo khoa Toán - Cơ - Tin học, trường Đại học Khoa Học Tự Nhiên - Đại Học Quốc Gia Hà Nội đãdạy bảo, cung cấp kiến thức bổ ích cho em suốt trình học tập nghiên cứu khoa Em xin cảm ơn thầy, cô giáo, cán Phòng Sau đại học, Trường Đại học Khoa học Tự nhiên - Đại Học QuốcGia Hà Nội tạo điều kiện thuận lợi trình thực luận văn Nhân dịp này, em cảm ơn gia đình, bạn bè động viên, tạo điều kiện cho em suốt trình học tập thực luận văn Hà Nội, ngày tháng năm 2015 Học viên Doãn Thu Hƣơng Mục Lục Mục Lục Lời mở đầu Chương 1: Phương trình tán sắc sóng mặt Rayleigh truyền đàn hồi trực hướng 1.1 Các phương trình truyền sóng 1.2 Trường hợp có hai mặt tự 1.3 Trường hợp có mặt tự do, mặt bị ngàm 10 Chương Các công thức xác định điểm tiếp xúc 13 2.1 Trường hợp có hai mặt tự 13 2.2 Trường hợp có mặt đáy bị ngàm 15 Trường hợp: C Trường hợp C  A   C 1 A 2  A 17 20 2.3 Tính trơn đường cong phổ vận tốc điểm tiếp xúc 23 Chương Trường hợp đẳng hướng ví dụ minh họa số 31 3.1 Tấm có hai biên tự 31 3.2 Trường hợp có mặt tự do, mặt đáy ngàm 33 3.3 Ví dụ minh họa số tập nghiệm điểm tiếp xúc S1 , S S3 35 Kết luận 39 Tài liệu tham khảo 40 Các công trình khoa học công bố 41 Lời mở đầu Phương trình tán sắc sóng mặt Rayleigh mô hình khác thường dẫn phương trình siêu việt dạng ẩn phụ thuộc vào hai biến vận tốc truyền sóng tần số sóng với tham số vật liệu mô hình Trong việc giải số tìm nghiệm phương trình tán sắc này, tần số sóng thường cho trước vận tốc truyền sóng tìm phương pháp số khác Nói chung, với giá trị tần số sóng, có nhiều nghiệm vận tốc nghiệm vận tốc ứng với mode truyền sóng khác sóng mặt Rayleigh Khi nghiệm vận tốc truyền sóng tìm với giá trị khác tần số sóng tranh miêu tả phụ thuộc chúng gọi đường cong phổ mode truyền sóng Thông thường đường cong phổ nằm xen kẽ Tuy nhiên số trường hợp đặc biệt giá trị tham số mô hình, tồn cặp đường cong (ứng với mode khác nhau) tiến gần “tiếp xúc” với Các điểm tiếp xúc điểm thuộc hai mode khác toán truyền sóng Rayleigh chúng điểm tương ứng với nghiệm bội phương trình tán sắc Có nhiều thuật ngữ tiếng Anh cho điểm đặc biệt “osculation points” hay “avoided crossing points” luận văn sử dụng thuật ngữ “điểm tiếp xúc” Những điểm tiếp xúc xuất toán truyền sóng Rayleigh mà xuất nhiều toán thuộc lĩnh vực khác vật lý lượng tử, vật lý chất rắn, học, với nhiều thuật ngữ khác (xem Kausel cộng sự, 2015, với tài liệu tham khảo báo) Nói chung điểm tiếp xúc nghiệm bội toán giá trị riêng tương ứng với lĩnh vực trên, chúng có số tính chất đặc biệt Trong lĩnh vực địa chấn, cụ thể phương pháp tỷ số H/V-là phương pháp liên quan đến sóng mặt Rayleigh, tính chất đặc biệt đường cong tỷ số H/V phát điểm tiếp xúc Đó điểm tiếp xúc, đường cong có điểm cực đại chuyển thành điểm không (xem Trần Thanh Tuấn, 2009) Do điểm cực đại điểm không hai điểm quan trọng phương pháp tỷ số H/V nên điểm tiếp xúc tập đường cong phổ vận tốc sóng mặt Rayleigh cần nghiên cứu Trong lĩnh vực địa chấn, điểm tiếp xúc quan sát thấy từ lâu (ví dụ Sezawa Kanai, 1935) công trình nghiên cứu lý thuyết điểm Theo Kausel cộng (2015) nói điểm tiếp xúc lĩnh vực địa chấn đề cập rõ ràng sách Levshin (1973) sau đề cập nhắc đến số công trình Forbriger (2006) Liu cộng (2009) Gần đây, số kết giải tích điểm tiếp xúc sóng Rayleigh đàn hồi, cụ thể công thức xác định điểm tiếp xúc, công bố Trần Thanh Tuấn (2009) bổ sung Kausel cộng (2015) Tuy nhiên công thức tìm cho trường hợp đàn hồi đẳng hướng Nội dung luận văn cao học tìm công thức điểm tiếp xúc sóng Rayleigh với điều kiện biên khác làm từ vật liệu trực hướng Hơn nữa, tính chất trơn phổ đường cong vận tốc điểm tiếp xúc khảo sát Luận văn phần mở đầu kết luận có chương Nội dung chương tìm phương trình tán sắc sóng Rayleigh truyền trong trường hợp có hai biên tự trường hợp có biên tự biên ngàm Chương khảo sát phương trình tán sắc tìm để tìm công thức xác định điểm tiếp xúc khảo sát tính trơn phổ đường cong vận tốc điểm tiếp xúc Chương trình bày kết nhận trường hợp đẳng hướng minh họa vài kết ví dụ số Chƣơng 1: Phƣơng trình tán sắc sóng mặt Rayleigh truyền đàn hồi trực hƣớng Chương sử dụng phương pháp truyền thống để tìm phương trình tán sắc sóng mặt Rayleigh truyền trực hướng Đầu tiên, phương trình trạng thái phương trình chuyển động trình bày lại theo sách chuyên khảo Sau đó, tùy vào điều kiện biên tấm, phương trình tán sắc sóng Rayleigh thiết lập Các phương trình tán sắc sử dụng việc nghiên cứu điểm tiếp xúc chương 1.1 Các phƣơng trình truyền sóng Xét toán trực hướng có độ dày h thông số vật liệu c 1 , c , c 2 , c 6 Sóng mặt Rayleigh truyền mặt phẳng theo trục x1 trùng với hướng tắt dần theo trục phẳng Trục O x nằm có phương trình x đáy có phương trình  h x2 x2  vuông góc với mặt mặt Do toán truyền sóng Rayleigh biến dạng phẳng nên trường chuyển dịch có dạng u i  u i ( x1 , x , t ) , ( i  1, ) , u ( x , x , t )  , (1.1) t thời gian Mối liên hệ ứng suất chuyển dịch cho (ví dụ xem Ting, 1996)  1  c 1 u ,1  c u ,  22  c u ,1  c 2 u , (1.2)   c 6 ( u ,  u ,1 ) dấu phẩy đạo hàm theo biến không gian Trong trường hợp không xét đến trọng lực phương trình chuyển động sóng Rayleigh có dạng  1 ,1   ,   u1 ,  ,1   Giả sử sóng lan truyền theo phương 22 ,2 x1 (1.3)   u2 với vận tốc c số sóng k , hàm chuyển dịch biểu diễn dạng ui  U i  x2  e ik ( x1  c t ) , ( i  1, ) (1.4) Thay dạng hàm chuyển dịch vào phương trình chuyển động (1.3) sau sử dụng phương trình trạng thái (1.2), ta thu hệ phương trình vi phân chuyển động U i ( x ) Giải hệ ta có nghiệm tổng quát hàm chuyển dịch có dạng (xem Phạm Chí Vĩnh Ogden, 2004) u  B1e u   B1e  B2e k b1 x  k b1 x  B 3e   1B2e k b1 x  k b1 x  B4e   B3e số tích phân B i ( i  1, ) k b3 x  k b3 x k b3 x   3B4e (1.5)  k b3 x nghiệm phương trình b1 , b c 2 c 6 b   ( c  c 6 )  c 2 ( X  c 1 )  c 6 ( X  c 6 )  b  ( c 1  X )( c 6  X )    với X  c 2 Chú ý phương trình trùng phương chung có bốn nghiệm phức b1 , b  b1 và  b b1 thực dương bi số thực âm, bi bi nói thực phức b3 chúng Nghĩa là, trường hợp chọn số phức có phần thực dương Nếu b (1.6) b i ( i  1, 3) số thực dương, phức, bi bi số số ảo có phần ảo dương Trong phương trình (1.5), ta ký hiệu  k   i  k  (U (1.7) / U 1)k với k  b k ( c1  c 6 ) c1  X  c 6 b k  c 22bk  c 66  X ( c1  c 6 ) b k , ( k  1, ) (1.8) Sử dụng đại lượng không thứ nguyên e1  c1 c 66 , e2  c 22 c 66 , e3  c1 , x  c 66 X (1.9) c 66 phương trình (1.6) có dạng e b   ( e  1)  e ( x  e1 )  ( x  1)  b  ( e1  x )(1  x )    2 (1.10) , ( k  1, ) (1.11) (1.8) có dạng k  b k ( e  1) e2bk   x e1  x  b k  ( e  1) b k Theo công thức Viet ta có: ( e  1)  e ( x  e1 )  ( x  1) S ( x )  b1  b   2 P ( x )  b1  b  2 , e2 ( e1  x )(1  x ) (1.12) e2 Các số hạng công thức hàm chuyển dịch (1.5) tương ứng với bốn thành phần sóng gồm hai sóng lên hai sóng xuống sóng qP qSV Phương trình tán sắc để xác định vận tốc truyền sóng c phụ thuộc vào tần số xác định từ điều kiện biên Trong phần chương này, hai trường hợp biên xem xét Đó trường hợp có hai mặt biên tự trường hợp có mặt tự mặt bị ngàm Hai trường hợp trường hợp tới hạn mô hình đặt bán không gian Trường hợp đầu trường hợp tới hạn bán không gian có độ cứng nhỏ, trường hợp sau trường hợp bán không gian có độ cứng lớn so với độ cứng 1.2 Trƣờng hợp có hai mặt tự Từ điều kiện tự ứng suất mặt mặt ta có  12 ( )   22  12 ( h )   (h)  22 (0 )  (1.13) Sử dụng công thức chuyển dịch (1.5) ứng suất (1.2) vào điều kiện biên thu hệ phương trình đại số số tích phân B , B , B , B dạng ma trận sau: M  [ B1 , B , B , B ]  T ma trận M M 1 (1.14) có dạng c 6 b1   c 2  b1   b  c 6 b1 e   b1  c 2  b1 e c 6 b1 c 6 b3 c 6 b3 c 2  b1 c 2 b3 c 2 b3 c 6 b1 e   b1 c 2  b1 e   b1 c 6 b3 e  b3 c 2 b3 e  b3 c 6 b3 e   b3 c 2 b3 e   b3       (1.15) với   k h Để hệ phương trình có nghiệm không tầm thường định thức tương ứng ma trận phải Từ ta thu phương trình tán sắc sóng mặt Rayleigh sau B0  B0 c o s h (  b1 ) c o s h (  b )  B0 B0 s in h (  b1 ) s in h (  b )  (1.16) B  b  ( S e  e  x )(1  x )  e x b1    (1.17) B  b1  ( S e  e  x )(1  x )  e x b    với S biểu diễn (1.12) Khi biểu diễn thông qua tham số tấm, phương trình tán sắc (1.16) có dạng c o s h (  b1 ) c o s h (  b )  B s in h (  b1 ) s in h (  b ) b1 1 (1.18) b3 với ( S e  e  x ) (1  x ) S  e x P S  e x ( S e  e  x )(1  x ) P B  P 2 ( S e  e  x ) (1  x )  e x P  e x S ( S e  e  x )(1  x ) S 2 (1.19) cho phương trình (1.12) 1.3 Trƣờng hợp có mặt tự do, mặt dƣới bị ngàm Từ điều kiện tự ứng suất mặt điều kiện ngàm mặt ta có u1 (0 )  u (0 )  ,  12 ( h )   (h)  22 (1.20) Tương tự trường hợp hai biên tự do, sử dụng công thức chuyển dịch (1.5) ứng suất (1.2) vào điều kiện biên thu hệ phương trình đại số số tích phân B , B , B , B dạng ma trận sau: Từ suy     x a1 , xa   s ig n  x  x  a1 t  , x   4t  , x  a1 a1 a1 a1   e2 P ( xa )     M ( xa ) e2  (2.72) với P ( xa )  Trong trường hợp này, thay hàm A(x) C ( x ) ta có A ( x a  e1    e2  1 e3  , M ( xa )  x  xa xa (2.73) vào biểu thức đạo hàm )   1, C ( x a )  1, A x ( x a )  C x ( x a )  1 Do từ phương trình (2.64) ta có F x  S     t1 (1  t )    t    t (1  t )   x x  x =   t1 t    t     t1 t  t  (2.74) x x x    t ( t1 t  1)     t  t1 t    S  x x   Do t1 t ( S )   Tương tự ta có: F   S     t  t1 t    S     Xét tập hợp nghiệm S2 (2.75) phƣơng trình(2.55) Trong trường hợp tập nghiệm này, vận tốc truyền sóng e1 e  e xa  e2  e3 2 , (2.76) ta có t3 ( S )  t3 ( a , x a )   2 Tương tự phần ta có (2.77)      ( , x a ) (  a , x a )  lim  a 2   a xa  xa  lim  a t3 ( , x a )  2   a ( e  1) 2 2 Từ phương trình (2.77), biểu thức không xác định điểm (  a  /  Tại lân cận (  a ta có khai triển Taylor hàm , xa ) (2.78) M ( xa ) 2 , xa ) có dạng t3 ( , x ) theo hai biến có dạng t3 ( , x )     t3  (  a , x a )( x  x a )   t 2 x 2   (    a )    t3 ( , x a )   ( t3 )  ( a , x a ) (   a )  o    a 2 2 2  (2.79) Ta có   t   S    ta n h   b3    ta n h    b3   b3   ta n h    (2.80)  t   t  b  t b ( x a ,  a ) 2 Ở ta sử dụng kết phương trình (2.77) Từ suy  Để tính    x    ( a , x a ) 2 t3  S  b3  S (  a , x a )  s ig n (    a ) 2 , ta sử dụng công thức   ( a , x )  x  x a ( e  1)  2  x  t3 ( a , x )     x  x  xa  xa M ( xa ) e2  ( x  xa )  2 với 2 x  xa t3 ( a , x ) 2  2 x  xa (2.81)  ( a , x ) (  a , x a )  lim a2 xác định từ (2.65) (2.66) Tại lân cận (  a , xa ) ( e  e3 ) e2 P ( x ) 2 M (x) (2.82) , từ khai triển Taylor (2.79) ta có t3 ( a , x )    t3 2  x (  a , x a )( x  x a )  o ( x  x a ) 2 2 (2.83) với t  x ( a , x a )  t3 ( a , x a ) a 2 2  b3  x ( x a ) (2.84) Do    x ( xa  xa ) (  a , x a )  s ig n ( x  x a ) 2 ( e  1) 2 2 t  (S )   2 2 x   ( e  e3 ) e P ( x a ) M ( xa ) 2 (2.85) P ( xa )   e1   e 2 e3   e3 2 Để tính toán đạo hàm riêng hàm lân cận nghiệm F ( , x ) e (  e  e  e )  e1 e ( e  1) A( xa )  A ( xa )  ' ; C ( x a )  1, ( e  e )(1  e ) 2 ta có 2 e3 ( e  e3 ) S2 (2.86) C ( xa )  ' ; 2 e3 Thay giá trị vào phương trình (2.64) ta nhận F x ( S )  A ( x a )  t1 (1  t )   S  x   t (1  t )  t   S  ,  x (2.87)  S2  (2.88) S2   F  ( S )  A ( x a )  t1 (1  t )    2 Đối với tập nghiệm S3 ,  t1 (1  t )  t    đạo hàm riêng hàm  F tính hoàn toàn tương tự với ý  a thay  a Tính chất đạo hàm vận tốc truyền sóng tần số điểm tiếp xúc Đối với ba tập nghiệm điểm tiếp xúc    x S1 , S S ta có ( a , x a )  s ig n ( x  x a ) R ( a , x a ) i i i i (2.89) i     ( a , x a )  s ig n (   a ) Q ( a , x a ), i i i i i  i  1,  (2.90) đó, hàm R ( a , x a ) i i Q (  a , x a ) ứng i với tập nghiệm cho i phương trình (2.72) (2.81), (2.85) Như vậy, ta thấy đạo hàm riêng hàm  ( , x ) s ig n (    a ) không liên tục điểm tiếp xúc hàm s ig n ( x  x a ) gián đoạn điểm Điều làm cho đạo hàm toàn phần vận tốc truyền sóng tần số bị gián đoạn Tuy nhiên, thấy đạo hàm riêng bên trái hàm phải hàm   ( , x )  ( , x ) đạo hàm riêng bên  Và từ công thức (2.61) (2.62), ta thấy điểm tiếp xúc, đạo hàm toàn phần bên trái đường cong vận tốc tần số mode đối xứng đạo hàm toàn phần bên phải mode phản đối xứng ngược lại Chính điều làm cho nhánh bên trái mode đối xứng dường nối liền với nhánh bên phải mode phản đối xứng Điều xảy với hai nhánh lại Hiện tượng quan sát tính toán số Sezawa Kanai (1935) Hiện tượng lý giải từ ý nghĩa vật lý toán truyền sóng sau Do vận tốc nhóm tính từ đạo hàm toàn phần vận tốc pha có dạng (ví dụ xem Achenbach (1973)) vg  d  c k dk dc (2.91) dk vận tốc nhóm đặc trưng cho vận tốc truyền lượng mode Hơn nữa, lượng đại lượng vật lý nên vận tốc thay đổi cách đột ngột Chính vậy, điểm tiếp xúc, nhánh mode đối xứng phải nối cách trơn với nhánh khác mode phản đối xứng Nếu điều không xảy ra, vận tốc truyền lượng bị gián đoạn điểm tiếp xúc Đối với điểm tiếp xúc thuộc trường hợp 3, điểm tiếp xúc loại này, tính trơn đường cong vận tốc chưa xác định Bước đầu biết lim  2 t3  t3  ( , x ) t3  t31 2  T ( , x ) t  t lim 2 2 t3  t3 t3  t31 2   Do đạo hàm đường cong vận tốc tần số ban đầu có dạng (2.92)       xác định công cụ giải tích tốt hơn, ví dụ sử dụng định lý L’Hospital Tuy nhiên, biểu thức t 321 t 322 có biểu diễn phức tạp nên luận văn không nghiên cứu tính trơn điểm tiếp xúc loại Chƣơng Trƣờng hợp đẳng hƣớng ví dụ minh họa số Trong trường hợp đẳng hướng, hệ số vật liệu vô hướng (1.9) có dạng (Phạm Chí Vĩnh Nguyễn Thị Khánh Linh, 2012) e1  e  /  , e  /   với     (3.1) đặc trưng cho số vật liệu đẳng hướng   vận tốc sóng dọc ngang tấm, nghiệm b1 , b phương trình đặc trưng (1.10) có dạng b1    x , b   x (3.2) 3.1 Tấm có hai biên tự Xét phương trình tán sắc trường hợp tự (1.16), trường hợp đẳng hướng ta có ( x  )  (  x )( x  1) B  8( x  2) 1 x 1 x (3.3) Khi đó, phương trình (1.16) đưa dạng phương trình tán sắc lớp đẳng hướng tự trình bày khóa luận tốt nghiệp Doãn Thu Hương (2011) Khi phương trình tán sắc (1.16) biểu diễn tách thành hai nhánh đối xứng phản đối xứng phương trình (2.6), ta có t1  B  B  t3 t1 t3  B0  B0  (3.4) Bằng cách thay biểu thức B (3.3) vào biểu diễn ý mối liên hệ ta n h ( x )   i ta n ( ix ) , ta nhận phương trình mode đối xứng phản đối xứng sau ta n (  x  1 / ) ta n ( x  1 / )   x 1  x 1 (x  2) (3.5) ta n ( ta n (  x  1 / ) x  1 / ) (x  2)   x 1  x 1 (3.6) Các phương trình trình bày nhiều sách chuyên khảo ví dụ Achenback (1973) Trong trường hợp có hai biên tự do,do điều kiện điểm tiếp xúc phức tạp nên ta xét trường hợp đẳng hướng có tham số vật liệu e1  e  , e  Đây trường hợp hệ số    /  Theo Phạm Chí Vĩnh Ogden (2004), môi trường đẳng hướng, để phương trình Rayleigh có nhiều nghiệm thực hệ số  phải thỏa mãn   Giá trị   ví dụ minh họa thỏa mãn điều kiện Với giá trị tham số trên, phương trình Rayleigh có nghiệm thực Nghiệm nhỏ x R  , hai nghiệm lại x  x  4.6 Nghiệm x R vận tốc truyền sóng Rayleigh bán không gian làm vật liệu Hai nghiệm lại vận tốc truyền sóng Rayleigh điểm tiếp xúc Tần số điểm tiếp xúc xác định từ phương trình (2.13) Hình 1: Các đường cong contour phổ vận tốc sóng Rayleigh tự điểm tiếp xúc Hình vẽ đường cong contour phương trình tán sắc (1.16) với tham số chọn Trong vùng tần số vận tốc khảo sát hình vẽ, quan sát điểm tiếp xúc Hai điểm đánh dấu hình vuông điểm tiếp xúc có vận tốc x bốn điểm đánh dấu tròn điểm tiếp xúc có vận tốc x1 Tần số điểm tiếp xúc xác định từ phương trình (2.13) 3.2 Trƣờng hợp có mặt tự do, mặt đáy ngàm Trong trường họp bị ngàm đáy,thay thông số vật liệu đẳng hướng (3.1) vào phương trình (1.25) ta có x  4x  A   4( x  2) (   1) x  x (   )  C   4( x  2)  x 1 x (3.7) Thay biểu thức hàm A ( x ) C ( x ) vào phương trình tán sắc trường hợp ngàm (1.23) ta dễ dàng nhận lại phương trình tán sắc sóng Rayleigh mô hình lớp có đáy bị ngàm trình bày phương trình (2.14) luận án tiến sỹ Trần Thanh Tuấn (2009) Đối với lớp nghiệm thứ điểm tiếp xúc nên  S1 (2.38), e3   / Do đó, từ công thức (2.28) ta có ca xa     (3.8)  ký hiệu vận tốc sóng ngang truyền lớp Tần số điểm tiếp xúc có dạng (từ phương trình (2.38)) a   (  m ) 2 ( m  ,1, , ) (3.9) Do   kh   a  2 f a h 1 ca  f a1 : fa h   ca  (  m)  (  m) (3.10) Kết trùng với kết nhận Trần Thanh Tuấn (2009) (xem phương trình 4.11) Đối với lớp nghiệm S2 , trường hợp đẳng hướng ta có \ e1 e  e 4 xa  e  e3 2 Từ điều kiện xa   4  (3.11) ta có điều kiện tham số vật liệu lớp để tồn lớp nghiệm điểm tiếp xúc   / Kết nhận Trần Thanh Tuấn (2009) Kausel cộng (2015) Tần số điểm tiếp xúc 4  a   Các ràng buộc R2 p  xác định từ phương trình (2.55) có dạng S2 (  q  ) 4  (1   ) , p  1, , , q  ,1, , (3.12) từ phương trình (2.56) số vật liệu lớp trường hợp   với điều kiện /   1/2 2p (3.13)  2(q  p ) với ý p  q  dẫn đến q   p  q  (3.14) Chú ý rằng, trường hợp đẳng hướng, phương trình xác định tần số điểm tiếp xúc (2.46) đưa hai phương trình (4.12) (4.13) Trần Thanh Tuấn (2009) Tuy nhiên, cách tìm tần số điểm tiếp xúc ràng buộc tham số vật liệu khác so với cách thực luận văn Về mặt công thức tìm luận văn điểm tiếp xúc (2.50) gọn so với công thức (4.15) (4.16) Trần Thanh Tuấn (2009), mặt toán học công thức tương đương với Trong Trần Thanh Tuấn (2009), phương trình (2.46), thay lấy tổ hợp tổng hiệu phương trình, phương trình chia cho Tương tự vậy, tần số điểm tiếp xúc thuộc tập nghiệm S3 xác 4   r 4   , r  1, , , s  ,1, ,   s    (1   ) 2  (3.15) định từ phương trình (2.55) có dạng a  Các ràng buộc R từ phương trình (2.57)của số vật liệu lớp trường hợp   với điều kiện /   1/2  s  r  S3 (3.16)  2(r  s) với ý r  s  dẫn đến Tập nghiệm 2s s  (3.17) trường hợp đẳng hướng không tìm Trần Thanh Tuấn (2009), nhiên bổ sung Kausel cộng (2015) 3.3 Ví dụ minh họa số tập nghiệm điểm tiếp xúc Để khảo sát số tập nghiệm có giá trị thỏa mãn S1 S1 , S S3 phương trình (2.38) ta chọn tham số e1  , e  , e  Hình 2: Các đường cong tán sắc với điểm tiếp xúc S Hình vẽ biểu diễn đường contour phương trình tán sắc (1.23) với thông số Các đường contour bao gồm các mode đối xứng phản đối xứng Bốn đường tròn nhỏ biểu thị cho bốn điểm tiếp xúc thuộc lớp nghiệm S tính toán từ phương trình (2.38) ta có x a  6.3 3 ; e a  1.2 , 3.7 , 6.2 , 8.7 (3.18) Để minh họa số tập nghiệm S2 ta chọn tham số phương trình ràng buộc R thỏa mãn p  q  Khi từ điều kiện ràng buộc (2.57) ta có e e ( e1  1) ( e  e )( e  1) Có nhiều tập giá trị  (3.19) mãn đẳng thức với điều kiện thêm e , e , e thỏa e1 , e , e  , e1 e  e  , e1  e  2 (3.20) Các điều kiện xuất phát từ điều kiện lượng biến dạng xác định dương vật liệu điều kiện x a  Giả sử ta chọn e1  , e  , e  truyền sóng điểm tiếp xúc Khi từ công thức (2.55) ta có vận tốc x a  / tần số a    /   Các kết minh họa Hình vẽ 3a Nếu ta chọn tham số khác có giá trị xa  /  a  /   4.3 e1  / , e  , e  Ta có Các kết minh họa Hình vẽ 3b Trong ví dụ minh họa ta chọn hệ số p q công thức (2.55) có giá trị p  1, q  Khi đó, tham số vật liệu phải thỏa mãn điều kiện sau để tồn lớp nghiệm điểm tiếp xúc S2 e e ( e1  1) ( e  e )( e  1) 2  21 25 (3.21) (a) (b) Hình 3: Các điểm tiếp xúc thuộc nhóm nghiệm Giả sử ta chọn hai giá trị tham số e1  1 / 8 , e2  , e3  / S2 với p  q  e1  1 / 4 , e  1, e  / Cả hai giá trị cho vận tốc điểm tiếp xúc (theo công thức (2.55)) x a  8 Và số thứ cho tần số (thay q  thay p  (2.55))  a  2.4 5 số thứ hai cho tần số  a  3.4 2 (a) Hình 4: Các điểm tiếp xúc thuộc nhóm nghiệm S2 với p  1, q  Một số kết số điểm tiếp xúc thuộc tập nghiệm S2 thường xảy mode bậc cao Và nhiều kết số có tập nghiệm S tồn điểm tiếp xúc xảy mode (mode 0), mặt toán học điều chưa chứng minh Kết quan trong địa vật lý nói chung thiết bị đo đạc đo tín hiệu mode mang phần lớn lượng sóng mặt Rayleigh Kết luận Luận văn khảo sát toán truyền sóng mặt Rayleigh trực hướng chịu hai điều kiện biên khác Phương trình tán sắc sóng Rayleigh trường hợp điều kiện biên nhận phương pháp truyền thống Các phương trình tán sắc sử dụng để khảo sát điểm tiếp xúc sóng mặt Rayleigh trường hợp Các kết đạt luận văn là: - - - - Nhận phương trình tán sắc sóng mặt Rayleigh truyền trực hướng hai trường hợp điều kiện biên: có hai mặt tự có mặt tự do, mặt ngàm Các công thức xác định điểm tiếp xúc trường hợp xác định sử dụng ý tưởng từ phương pháp lý thuyết tia Tolstoy Usdin (1953) Đã khảo sát đạo hàm đường cong phổ vận tốc sóng mặt Rayleigh điểm tiếp xúc Kết điểm tiếp xúc, đường cong phổ vận tốc tính trơn Đạo hàm chúng điểm không liên tục Tuy nhiên, đạo hàm trái đường cong phổ vận tốc điểm tiếp xúc mode đối xứng với đạo hàm phải tương ứng mode phản đối xứng Tương tự vậy, đạo hàm phải đường cong phổ vận tốc điểm tiếp xúc mode đối xứng với đạo hàm trái tương ứng mode phản đối xứng Trong trường hợp đẳng hướng, kết nhận luận văn đưa kết nhận tác giả khác Đã khảo sát số số trường hợp nghiệm điểm tiếp xúc Các kết đạt luận văn có ý nghĩa khoa học Tài liệu tham khảo 10 11 12 13 Doãn Thu Hương (2011) , “Sóng Rayleigh Lam truyền môi trường môi trường không theo phương z” Khóa luận tốt nghiệp ngành học, Trường ĐH Khoa Học Tự Nhiên Achenbach, J D "Waves in elastic solids." Nord Holland, Amsterdam (1973) Forbriger, Thomas "Einige Gedanken zu: Oskulationen von Dispersionskurven, Entartung und Hybridisierung von Moden." (2006) Kausel, Eduardo, Peter Malischewsky, and João Barbosa "Osculations of spectral lines in a layered medium." Wave Motion 56 (2015): 22-42 Levshin, A L "Surface and channel seismic waves." Nauka, Moscow (1973) Liu, Xue‐Feng, You‐Hua Fan, and Xiao‐Fei Chen "Research on the Cross of the Dispersion Curves of Rayleigh Waves and Multi‐ModesCoupling Phenomenon." Chinese Journal of Geophysics 52.5 (2009): 994-1002 Sezawa, Katsutada, and Kiyoshi Kanai "Discontinuity in Dispersion Curves of Rayleigh-Waves." Proceedings of the Imperial Academy 11.1 (1935): 1314 Ting.T.C.T (1996) Anisotropic Elasticity: Theory and Applications, Oxford Unversity Press NewYork Tolstoy, Ivan, and Eugene Usdin "Dispersive properties of stratified elastic and liquid media: A ray theory." Geophysics 18.4 (1953): 844-870 Tran Thanh Tuan "The ellipticity (H/V-ratio) of Rayleigh surface waves." PhD diss., 2009 Vinh, Pham Chi, and Nguyen Thi Khanh Linh "An approximate secular equation of Rayleigh waves propagating in an orthotropic elastic half-space coated by a thin orthotropic elastic layer." Wave motion 49.7 (2012): 681-689 Vinh, Pham Chi, and R W Ogden "On formulas for the Rayleigh wave speed."Wave Motion 39.3 (2004): 191-197 Vinh, Pham Chi, and R W Ogden "Formulas for the Rayleigh wave speed in orthotropic elastic solids." Archives of Mechanics 56.3 (2004): 247-265 Các công trình khoa học công bố Trần Thanh Tuấn, Peter Malischewsky, Doãn Thu Hương (2013) Tính chất tỷ số H/V điểm osculation mô hình lớp có đáy bị ngàm Tuyển tập Hội nghị Khoa học toàn quốc Cơ học Vật rắn biến dạng lần thứ XI Thành phố Hồ Chí Minh, 7-9/11/2013, p.1275-1282 \ ... vận tốc truyền sóng tần số sóng với tham số vật liệu mô hình Trong việc giải số tìm nghiệm phương trình tán sắc này, tần số sóng thường cho trước vận tốc truyền sóng tìm phương pháp số khác Nói... NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - DOÃN THU HƢƠNG ĐIỂM OSCULATION CỦA SÓNG RAYLEIGH TRONG MỘT SỐ MÔ HÌNH Mã số: 60 44 01 07 Chuyên ngành: Cơ học vật thể rắn LUẬN VĂN THẠC SĨ KHOA... chung, với giá trị tần số sóng, có nhiều nghiệm vận tốc nghiệm vận tốc ứng với mode truyền sóng khác sóng mặt Rayleigh Khi nghiệm vận tốc truyền sóng tìm với giá trị khác tần số sóng tranh miêu tả