Biên Hấp Thụ Tương Hỗ Miền Thời Gian

Luận án tiến sĩ phân tích time domain reciprocal absorbing boundaries, xây dựng cơ sở lý luận, kiểm chứng thực nghiệm, đóng góp tri thức mới cho ngành.

• Tác giả: Cuong Tan Nguyen
• Người hướng dẫn: John L. Tassoulas, Supervisor
• Trường học: The University Of Texas At Austin
• Chuyên ngành: Engineering
• Đề tài: Time-domain Reciprocal Absorbing Boundaries
• Loại tài liệu: Dissertation
• Năm xuất bản: 2017
• Địa điểm: Austin

Mục lục chi tiết

• Acknowledgments
• Abstract
• List of Tables
• List of Figures
• 1. CHƯƠNG 1: Reciprocal absorbing boundary condition for the time-domain numerical analysis of wave motion in unbounded layered media
• 1.1. Notation and definitions
• 1.2. Time-domain reciprocity theorems for elastodynamics
• 1.3. Derivation of reciprocal absorbing boundary condition
• 1.4. Time-stepping algorithm
• 1.5. Reciprocal absorbing boundary condition (RABC)
• 1.5.1. Semi-infinite beam on simple supports
• 1.5.2. Semi-infinite beam on flexible supports
• 1.5.3. Semi-infinite rod
• 1.5.4. Stratum in plane strain
• 1.5.5. Stratum in antiplane shear
• 2. CHƯƠNG 2: Application of reciprocal absorbing boundary condition to transient analysis of acoustic wave propagation
• 2.1. Reciprocity theorems for acoustics
• 2.1.1. The linear wave equation for velocity potential
• 2.1.2. Acoustic reciprocity theorems
• 2.1.2.1. Reciprocity in the frequency domain
• 2.1.2.2. Reciprocity in the time domain
• 2.2. Reciprocal absorbing boundary condition
• 2.2.1. Finite elements for acoustic wave propagation
• 2.2.1.1. The semi-discrete form
• 2.2.1.2. Time-stepping algorithm
• 2.2.2. Derivation of reciprocal absorbing boundary condition
• 2.2.2.1. The original semi-infinite problem
• 2.2.2.2. The construction of reciprocal absorbing boundary condition
• 2.2.2.2.1. One-dimensional problem
• 2.2.2.2.2. Two-dimensional problem
• 3. CHƯƠNG 3: Reciprocal absorbing boundary condition with perfectly matched discrete layers for SH waves in a layered half-space
• 3.1. Reciprocal absorbing boundary condition
• 3.2. Derivation of reciprocal absorbing boundary condition
• 3.3. Time-stepping algorithm
• 3.4. Reciprocal absorbing boundary condition
• 3.5. Perfectly Matched Discrete Layers
• 3.5.1. Finite element discretization
• 3.5.2. Reciprocal absorbing boundary condition with perfectly matched discrete layers
• 3.5.2.1. A two-layer stratum truncated by RABC
• 3.5.2.2. A homogeneous layer truncated by PMDLs
• 3.5.2.3. A homogeneous half-space truncated by RABC and PMDLs
• 3.6. A layered half-space truncated by RABC & PMDLs
• 4. CHƯƠNG 4: Reciprocal absorbing boundary condition with perfectly matched discrete layers for SV-P waves in a layered half-space
• 4.1. Reciprocal absorbing boundary condition for a stratum
• 4.1.1. Derivation on the basis of the displacement-based reciprocity theorem
• 4.1.2. Derivation based on the dynamic stiffness matrix
• 4.2. Derivation based on the acceleration unit-impulse response matrix
• 4.3. Convolution quadrature and time-stepping scheme
• 4.4. The construction of reciprocal absorbing boundary condition (RABC)
• 4.5. Computing initial values
• 4.6. Perfectly matched discrete layers for SV-P waves
• 4.6.1. Finite element discretization
• 4.7. Reciprocal absorbing boundary condition with perfectly matched discrete layers
• 4.7.1. RABC for a two-layer stratum
• 4.7.2. A homogeneous half-space
• 4.7.3. A layered half-space
• 5. CHƯƠNG 5: Transient analysis of full-space unbounded domains by reciprocal absorbing boundaries
• 5.1. Reciprocal absorbing boundary condition (RABC)
• 5.1.1. Acceleration based reciprocity theorem
• 5.1.2. The construction of Reciprocal Absorbing Boundary Condition
• 5.1.3. Perfectly matched discrete layers (PMDLs)
• 5.1.3.1. Semidiscrete equation of motion
• 5.1.3.1.1. Anti-plane shear
• 5.1.3.2. Steady-state phase for harmonic forces
• 5.2. Reciprocal absorbing boundary condition with perfectly matched discrete layers
• Bibliography
• Vita

Tóm tắt

I. Giới Thiệu Chung Về Biên Hấp Thụ Tương Hỗ Miền Thời Gian

Bài toán truyền sóng trong môi trường không giới hạn đặt ra nhiều thách thức trong mô phỏng số. Trong thực tế, khu vực quan tâm thường chỉ là một phần nhỏ, gọi là "near field", của một không gian rộng lớn, thậm chí vô hạn. Để giải quyết vấn đề này, biên hấp thụ được đưa vào để mô phỏng sự hấp thụ sóng vào miền bên ngoài. Các biên hấp thụ này thể hiện mức độ chính xác khác nhau của "far field" khi nhìn từ "near field", nơi thường chứa các yếu tố bất thường và phi tuyến. Do đó, việc nghiên cứu và phát triển các phương pháp biên hấp thụ hiệu quả là vô cùng quan trọng. Các biên hấp thụ cho phân tích miền thời gian của truyền sóng trong môi trường không giới hạn có thể được phân loại thành hai loại: cục bộ và toàn cục. Nghiên cứu này tập trung vào việc phát triển các biên hấp thụ "toàn cục" mới bằng cách sử dụng các định lý tương hỗ, như tài liệu gốc đề cập.

1.1. Tại sao cần biên hấp thụ tương hỗ cho mô phỏng truyền sóng

Việc sử dụng biên hấp thụ tương hỗ (RABC) cho phép mô phỏng chính xác hơn truyền sóng trong môi trường không giới hạn. Thay vì mô phỏng toàn bộ không gian vô hạn, chỉ cần mô phỏng một miền hữu hạn với RABC đặt tại biên. RABC mô phỏng chính xác ảnh hưởng của miền vô hạn bên ngoài lên miền hữu hạn bên trong, loại bỏ hoặc giảm thiểu sự phản xạ sóng giả tạo tại biên, đảm bảo kết quả mô phỏng chính xác hơn.

1.2. Ứng dụng của biên hấp thụ trong các lĩnh vực kỹ thuật nào

Biên hấp thụ có ứng dụng rộng rãi trong nhiều lĩnh vực kỹ thuật, bao gồm địa chấn học và kỹ thuật động đất (dự đoán chuyển động mặt đất mạnh, tương tác đất-công trình động), các thí nghiệm không phá hủy (đánh giá tình trạng đường cao tốc và đường băng sân bay, phát hiện và mô tả các khuyết tật trong các cấu trúc dài như cọc và đường ống, và xác định tính chất đất từ các phép đo chuyển động sóng tại các công trường xây dựng), âm học (tán xạ sóng âm trong không khí và nước), hàng không (dòng chảy xung quanh cánh máy bay), kỹ thuật điện (tán sắc sóng điện từ), thủy động lực học (khúc xạ sóng nước, tương tác cấu trúc-chất lỏng, vận chuyển chất ô nhiễm), khí tượng học (dự báo thời tiết toàn cầu và cục bộ), hải dương học (lập bản đồ đáy biển), thăm dò địa vật lý (di cư địa chấn)...

II. Thách Thức Khi Mô Phỏng Truyền Sóng Trong Môi Trường Không Giới Hạn

Mô phỏng truyền sóng trong môi trường không giới hạn đối mặt với nhiều khó khăn. Việc tính toán trực tiếp trên miền vô hạn là bất khả thi. Việc cắt miền tính toán và áp dụng các điều kiện biên đơn giản thường dẫn đến phản xạ sóng giả tại biên, làm sai lệch kết quả mô phỏng. Đòi hỏi các phương pháp tính toán phải hiệu quả và chính xác để mô phỏng chính xác sự lan truyền của sóng mà không bị ảnh hưởng bởi biên giới nhân tạo. Các phương pháp như FDTD, FEM, MOM đòi hỏi tài nguyên tính toán lớn và có thể trở nên không hiệu quả đối với các bài toán lớn và phức tạp. Việc phát triển các biên hấp thụ tương hỗ hiệu quả là rất quan trọng để giải quyết những thách thức này.

2.1. Tại sao phản xạ sóng giả tại biên lại là vấn đề lớn

Phản xạ sóng giả tại biên có thể gây ra sai lệch nghiêm trọng trong kết quả mô phỏng. Sóng phản xạ từ biên sẽ tương tác với sóng tới, tạo ra các hiện tượng giao thoa và nhiễu xạ không mong muốn. Điều này đặc biệt quan trọng trong các ứng dụng như thiết kế anten, nơi mà sự chính xác của mô phỏng là yếu tố then chốt.

2.2. Hạn chế của các phương pháp tính toán điện từ trường truyền thống

Các phương pháp tính toán điện từ trường truyền thống như FDTD, FEM, và MOM có thể tốn kém về mặt tính toán và khó áp dụng cho các bài toán lớn và phức tạp. Ví dụ, phương pháp FDTD yêu cầu lưới tính toán mịn để đảm bảo độ chính xác, dẫn đến số lượng phần tử lớn và thời gian tính toán kéo dài. Các phương pháp khác như MOM có thể gặp khó khăn trong việc xử lý các cấu trúc có hình dạng phức tạp.

III. Phương Pháp Biên Hấp Thụ Tương Hỗ Cho Bài Toán Truyền Sóng

Nghiên cứu này đề xuất một phương pháp mới dựa trên định lý tương hỗ để xây dựng biên hấp thụ cho bài toán truyền sóng trong miền thời gian. Ý tưởng chính là lợi dụng tính tương hỗ giữa các miền vô hạn chỉ khác nhau về vị trí biên. Điều này cho phép xây dựng biên hấp thụ một cách trực tiếp trong miền thời gian mà không cần chuyển đổi qua lại giữa miền tần số. Phương pháp này đã được áp dụng thành công cho nhiều bài toán khác nhau, bao gồm dầm bán vô hạn, thanh bán vô hạn, và lớp vật liệu có lớp trong trạng thái biến dạng phẳng và trượt ngang.

3.1. Định lý tương hỗ được áp dụng như thế nào trong phương pháp này

Định lý tương hỗ cho phép thiết lập mối quan hệ giữa hai trạng thái của hệ thống. Trong trường hợp này, hai trạng thái khác nhau chỉ ở vị trí biên. Bằng cách áp dụng định lý tương hỗ, có thể thiết lập một điều kiện biên mà mô phỏng chính xác ảnh hưởng của miền vô hạn bên ngoài lên miền hữu hạn bên trong.

3.2. Ưu điểm của phương pháp miền thời gian so với miền tần số

Phương pháp miền thời gian cho phép mô phỏng trực tiếp sự lan truyền của sóng theo thời gian, phù hợp với các bài toán có nguồn kích thích thay đổi theo thời gian. Hơn nữa, phương pháp miền thời gian thường dễ dàng hơn trong việc xử lý các vật liệu phi tuyến và các hiện tượng phức tạp khác so với phương pháp miền tần số.

IV. Kết Hợp RABC và PMDLs Cho Môi Trường Bán Vô Hạn Phân Lớp

Nghiên cứu này tiếp tục mở rộng việc sử dụng Biên Hấp Thụ Tương Hỗ (RABC) bằng cách kết hợp nó với Perfectly Matched Discrete Layers (PMDLs) để giải quyết các bài toán truyền sóng SH (Shear Horizontal) trong môi trường bán vô hạn phân lớp. Sự kết hợp này cho phép mở rộng phạm vi ứng dụng của RABC cho các môi trường bán vô hạn và tránh sử dụng các phần tử hấp thụ góc. PMDLs là một kỹ thuật được sử dụng để hấp thụ sóng tại biên của miền tính toán mà không gây ra phản xạ, kết hợp với RABC giúp tăng cường hiệu quả và độ chính xác của mô phỏng. Sự kết hợp này mang lại kết quả tốt trong việc mô phỏng truyền sóng trong các môi trường phức tạp.

4.1. Tại sao cần kết hợp RABC với Perfectly Matched Discrete Layers PMDLs

Việc kết hợp RABC với PMDLs giúp cải thiện đáng kể hiệu quả và độ chính xác của mô phỏng. PMDLs có khả năng hấp thụ sóng tốt hơn ở các góc tới lớn, trong khi RABC cung cấp một điều kiện biên chính xác hơn cho miền tính toán. Sự kết hợp này giúp giảm thiểu phản xạ sóng giả và tăng độ ổn định của phương pháp số.

4.2. Ứng dụng của RABC và PMDLs trong kỹ thuật địa chấn

Kỹ thuật địa chấn có thể ứng dụng RABC và PMDLs, đặc biệt là trong việc mô phỏng sự lan truyền của sóng địa chấn trong lòng đất. Các lớp đất khác nhau có thể được mô hình hóa như các lớp phân lớp, và RABC cùng PMDLs có thể được sử dụng để mô phỏng chính xác ảnh hưởng của miền vô hạn bên dưới lên miền hữu hạn bên trên. Điều này giúp cải thiện độ chính xác của việc dự đoán chuyển động mặt đất trong các trận động đất.

V. Phân Tích Chuyển Động Tạm Thời Trong Miền Không Giới Hạn Bằng RABC

Nghiên cứu này cũng trình bày về việc ứng dụng Biên Hấp Thụ Tương Hỗ (RABC) để phân tích chuyển động tạm thời trong các miền không giới hạn. Định lý tương hỗ dựa trên gia tốc được sử dụng để xây dựng RABC. Perfectly Matched Discrete Layers (PMDLs) cũng được sử dụng để tăng cường khả năng hấp thụ sóng. Kết quả cho thấy phương pháp này có thể mô phỏng chính xác sự lan truyền của sóng trong các miền không giới hạn, mang lại những hiểu biết sâu sắc về hành vi của sóng trong các môi trường phức tạp. Phương pháp này mở ra tiềm năng cho việc nghiên cứu các hiện tượng truyền sóng trong nhiều lĩnh vực khác nhau.

5.1. Tại sao sử dụng định lý tương hỗ dựa trên gia tốc

Định lý tương hỗ dựa trên gia tốc có thể đơn giản hóa việc xây dựng RABC trong một số trường hợp. Gia tốc là một đại lượng vật lý quan trọng trong các bài toán động lực học, và việc sử dụng nó trong định lý tương hỗ có thể dẫn đến các công thức đơn giản và hiệu quả hơn.

5.2. So sánh hiệu quả của RABC và PMDLs trong phân tích chuyển động tạm thời

RABC và PMDLs đều có những ưu điểm và nhược điểm riêng. RABC cung cấp một điều kiện biên chính xác hơn, trong khi PMDLs có khả năng hấp thụ sóng tốt hơn ở các góc tới lớn. Việc kết hợp cả hai phương pháp có thể mang lại kết quả tốt nhất trong phân tích chuyển động tạm thời.

VI. Kết Luận Tiềm Năng và Hướng Phát Triển Của Biên Hấp Thụ Tương Hỗ

Nghiên cứu về Biên Hấp Thụ Tương Hỗ (RABC) đã chứng minh tiềm năng to lớn trong việc giải quyết các bài toán truyền sóng trong môi trường không giới hạn. Phương pháp này cung cấp một giải pháp hiệu quả và chính xác để mô phỏng sự lan truyền của sóng, đồng thời giảm thiểu các ảnh hưởng tiêu cực do biên giới nhân tạo. Việc kết hợp RABC với các kỹ thuật khác như PMDLs mở ra những hướng phát triển mới cho việc nghiên cứu các hiện tượng truyền sóng phức tạp trong nhiều lĩnh vực khoa học và kỹ thuật. Trong tương lai, nghiên cứu có thể tập trung vào việc cải thiện hiệu suất tính toán, mở rộng phạm vi ứng dụng, và phát triển các phương pháp RABC mới cho các loại sóng khác nhau.

6.1. Hướng phát triển của RABC trong tương lai

Hướng phát triển của RABC trong tương lai có thể bao gồm việc cải thiện hiệu suất tính toán bằng cách sử dụng các thuật toán song song và các kỹ thuật tối ưu hóa. Ngoài ra, việc mở rộng phạm vi ứng dụng của RABC cho các loại sóng khác nhau, chẳng hạn như sóng điện từ và sóng âm, cũng là một hướng nghiên cứu quan trọng. Cuối cùng, việc phát triển các phương pháp RABC mới dựa trên các nguyên lý vật lý khác nhau có thể dẫn đến những đột phá trong lĩnh vực này.

6.2. Tại sao Biên Hấp Thụ Tương Hỗ là một công cụ quan trọng cho nghiên cứu

Biên Hấp Thụ Tương Hỗ (RABC) cung cấp một công cụ mạnh mẽ cho các nhà nghiên cứu và kỹ sư trong nhiều lĩnh vực. Nó cho phép mô phỏng chính xác các hiện tượng truyền sóng trong môi trường không giới hạn, giúp hiểu rõ hơn về hành vi của sóng và thiết kế các hệ thống hiệu quả hơn. Ví dụ, trong kỹ thuật viễn thông, RABC có thể được sử dụng để thiết kế anten hiệu quả hơn. Trong kỹ thuật địa chấn, RABC có thể được sử dụng để dự đoán chuyển động mặt đất trong các trận động đất. Trong y học, RABC có thể được sử dụng để mô phỏng sự lan truyền của sóng siêu âm trong cơ thể con người.

Nội dung chính

Tổng quan luận án

Luận án tiến sĩ "Time-domain Reciprocal Absorbing Boundaries" của nghiên cứu sinh Cuong Tan Nguyen, thực hiện dưới sự hướng dẫn của Giáo sư John L. Tassoulas tại Đại học Texas tại Austin (The University of Texas at Austin, bảo vệ tháng 12 năm 2017), tập trung giải quyết bài toán mô phỏng số sự lan truyền sóng trong các môi trường vô hạn (unbounded media). Trong cơ học tính toán ứng dụng cho kỹ thuật địa chấn, tương tác động lực học đất - kết cấu (soil-structure interaction), kiểm tra không phá hủy và âm học, miền nghiên cứu hữu hạn (near field) thường chỉ là một phần nhỏ kết nối với một môi trường vô hạn rộng lớn (far field). Để tính toán số trị, miền vô hạn bắt buộc phải được cắt tỉa bằng các điều kiện biên nhân tạo nhằm hấp thụ hoặc cho phép sóng phản xạ truyền qua mà không dội ngược trở lại.

Khoảng trống nghiên cứu được tác giả xác định xuất phát từ giới hạn của các phương pháp hiện hữu. Các điều kiện biên truyền nhất quán (Consistent Transmitting Boundary - CTB) phát triển rất hoàn thiện trong miền tần số (frequency domain), nhưng việc chuyển đổi sang miền thời gian (time domain) đòi hỏi các phép biến đổi Fourier xuôi và ngược cồng kềnh, gặp trở ngại nghiêm trọng khi hệ thống xuất hiện các yếu tố phi tuyến tính ở miền gần (như đất phi tuyến hoặc tiếp xúc trượt ma sát tại đáy móng). Mặt khác, phương pháp phương trình tích phân biên trong miền thời gian đòi hỏi phải biết trước nghiệm cơ bản (Green function) giải tích dạng tường minh—vốn cực kỳ phức tạp hoặc bất khả thi đối với các môi trường phân lớp. Luận án giải quyết khoảng trống này bằng cách đề xuất một điều kiện biên hấp thụ toàn cục (global absorbing boundary condition) trực tiếp trong miền thời gian dựa trên các định lý tương hỗ đàn hồi động học, loại bỏ hoàn toàn yêu cầu phải có hàm Green giải tích từ trước.

Mục tiêu và câu hỏi nghiên cứu của luận án được triển khai cụ thể qua các nhiệm vụ:

1. Thiết lập cơ sở lý thuyết cho định lý tương hỗ dựa trên gia tốc (acceleration-based reciprocity theorem) trực tiếp trong miền thời gian cho môi trường đàn hồi tuyến tính và môi trường âm học.
2. Xây dựng công thức đại số và thuật toán tích bước thời gian cho Điều kiện biên hấp thụ tương hỗ (Reciprocal Absorbing Boundary Condition - RABC) kết hợp cùng phương pháp phần tử hữu hạn.
3. Ứng dụng RABC vào bài toán sóng âm trong ống dẫn sóng (acoustic waveguides) sử dụng hàm thế vận tốc.
4. Phát triển mô hình lai ghép RABC với các lớp rời rạc phối hợp hoàn hảo (Perfect Matched Discrete Layers - PMDLs) nhằm xử lý sóng cắt ngang SH và sóng kết hợp SV-P trong bán không gian phân lớp (layered half-space), loại bỏ nhu cầu sử dụng các phần tử góc (corner absorbing elements).
5. Mở rộng mô hình RABC-PMDL cho miền toàn không gian (full-space) chịu các chế độ kích động dao động dọc trục, xuyên tâm, xoắn và dao động ngang.

Đối tượng nghiên cứu của luận án là trường sóng cơ học đàn hồi động học (elastodynamic wave fields), bao gồm sóng cắt ngoài mặt phẳng (anti-plane shear / SH), sóng biến dạng phẳng (plane strain / SV-P) và sóng âm tuyến tính. Phạm vi không gian nghiên cứu bao gồm: thanh và dầm bán vô hạn, tầng đất phân lớp trên nền cứng (layered strata), bán không gian phân lớp (layered half-space) và toàn không gian hai chiều (full-space). Phạm vi thời gian tập trung vào phân tích phản ứng quá độ (transient analysis) dưới tác động của các tải trọng thay đổi theo thời gian (hàm bước nhảy Heaviside, xung tam giác, xung chữ nhật, sóng con Ricker) và giai đoạn dao động điều hòa xác lập (steady-state phase).

Tổng quan tài liệu và vị trí của luận án

Luận án điểm lại hệ thống các công trình nghiên cứu nền tảng về định lý tương hỗ và điều kiện biên hấp thụ trong cơ học sóng:

Về các định lý tương hỗ động học đàn hồi, tác giả viện dẫn công trình khởi đầu của Graffi (1963) [16] thiết lập quan hệ tương hỗ dạng tích chập thời gian cho vật thể đàn hồi tuyến tính, đồng nhất và đẳng hướng. Wheeler và Sternberg (1968) [53] mở rộng định lý này cho các miền vô hạn với khái niệm trạng thái đàn hồi động học có quá khứ tĩnh lặng (quiescent past). De Hoop (1966) [3] tổng quát hóa định lý tương hỗ tích chập thời gian cho môi trường bất đồng nhất, dị hướng và nhớt đàn hồi, trong khi Achenbach (1975) [1] hệ thống hóa các quan hệ tương hỗ trong trạng thái điều hòa thời gian.

Về biên truyền sóng nhất quán (Consistent Transmitting Boundary - CTB), Lysmer và Drake (1971) [31], Lysmer và Waas (1972) [33] đã đề xuất mô hình trong miền tần số cho tầng đất phân lớp chịu sóng cắt ngoài mặt phẳng thông qua việc giải bài toán trị riêng theo phương thẳng đứng. Waas (1972) [52] tiếp tục hoàn thiện CTB cho sóng biến dạng phẳng và sóng đối xứng trục. Kausel và Roësset (1975, 1977) [25, 26] mở rộng CTB cho các chuyển động sóng trụ tổng quát và phát triển các phần tử siêu cấp (hyperelements). Tassoulas và Kausel (1983) [44] mở rộng phần tử siêu cấp cho miền có điều kiện biên không thuần nhất. Kausel và Peek (1982) [24] thiết lập hệ hàm Green bán rời rạc cho môi trường phân lớp. Ngoài ra, Lotfi và cộng sự (1987) [30] đã mở rộng CTB cho hệ đập nước - hồ chứa, và Bougacha cùng Tassoulas (1991) [7] áp dụng cho môi trường phân lớp rỗng hai pha.

Về hạn chế của các điều kiện biên cục bộ (local absorbing boundary conditions), tác giả trích dẫn nghiên cứu của Ditkowski và Sever (2007) [13], chỉ ra rằng các điều kiện biên cục bộ không thể hiện được quán tính và độ cứng thực tế mà miền ngoài truyền cho miền trong.

Khoảng trống luận án lựa chọn giải quyết nằm ở việc các công thức CTB kinh điển gắn chặt với miền tần số, trong khi các phương pháp biên tích phân miền thời gian lại lệ thuộc vào hàm Green giải tích. Luận án định vị một giải pháp mới: sử dụng tính tương hỗ giữa hai miền vô hạn chỉ khác nhau bởi sự tịnh tiến vị trí biên thẳng đứng. Nhờ đó, phương pháp RABC cho phép mô phỏng trực tiếp trong miền thời gian mà hàm Green được tự động tích hợp trong quá trình giải, không cần xác định trước bằng công thức giải tích.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án dựa trên lý thuyết đàn hồi động học tuyến tính (linear elastodynamics). Trạng thái đàn hồi động học được định nghĩa theo Wheeler và Sternberg: một cặp có thứ tự $\mathcal{F} = [\mathbf{u}, \boldsymbol{\sigma}]$ gồm trường chuyển vị $\mathbf{u}(\mathbf{x}, t)$ và tensor ứng suất $\boldsymbol{\sigma}(\mathbf{x}, t)$ thỏa mãn:

• Phương trình chuyển động: $\nabla \cdot \boldsymbol{\sigma} + \mathbf{f} = \rho \ddot{\mathbf{u}}$
• Định luật cấu hooks liên hệ ứng suất - biến dạng: $\boldsymbol{\sigma} = \rho(c_d^2 - 2c_s^2)(\nabla \cdot \mathbf{u})\mathbf{I} + 2\rho c_s^2 \boldsymbol{\epsilon}$ trong đó $\rho$ là khối lượng riêng, $c_d$ là vận tốc sóng nén (dilatational wave velocity), $c_s$ là vận tốc sóng cắt (shear wave velocity), $\mathbf{f}$ là mật độ lực khối, và $\boldsymbol{\epsilon} = \frac{1}{2}(\nabla \mathbf{u} + (\nabla \mathbf{u})^T)$ là tensor biến dạng nhỏ. Trạng thái có quá khứ tĩnh lặng thỏa mãn $\mathbf{u} = \mathbf{0}$ khi $t \le 0$, ký hiệu $\mathcal{F} \in \mathcal{E}_0(\mathbf{f}, \rho, c_d, c_s; \bar{R})$.

Khung phân tích chính là định lý tương hỗ dựa trên gia tốc (acceleration-based reciprocity theorem) do tác giả chứng minh bằng cách áp dụng biến đổi Laplace cho hai trạng thái $\mathcal{F}^I$ và $\mathcal{F}^E$, sau đó lấy biến đổi Laplace ngược thông qua định lý tích chập: $$\int_{\Gamma} t_j^I * a_j^E , d\Gamma + \int_R f_j^I * a_j^E , dR + \rho \int_R [u_j^I(\mathbf{x}, 0^+) \dot{a}j^E + v_j^I(\mathbf{x}, 0^+) a_j^E] , dR = \int{\Gamma} t_j^E * a_j^I , d\Gamma + \int_R f_j^E * a_j^I , dR + \rho \int_R [u_j^E(\mathbf{x}, 0^+) \dot{a}j^I + v_j^E(\mathbf{x}, 0^+) a_j^I] , dR$$ trong đó $t_j = \sigma{ij} n_i$ là lực mặt, $a_j = \ddot{u}_j$ là gia tốc, và dấu $*$ biểu diễn tích phân tích chập thời gian: $$\phi(t) * \psi(t) = \int_0^t \phi(t - \tau) \psi(\tau) , d\tau$$

Phương pháp nghiên cứu thực tế là mô phỏng số định lượng, sử dụng các công cụ và kỹ thuật:

• Rời rạc hóa không gian: Phương pháp phần tử hữu hạn (FEM).
• Tính toán tích chập thời gian: Sử dụng quy tắc hình thang (trapezoidal rule) với bước thời gian rời rạc $\Delta t$, $t_k = (k-1)\Delta t$.
• Tích phân bước thời gian: Sử dụng thuật toán Newmark ở chế độ gia tốc trung bình không đổi (constant-average acceleration scheme) với các tham số $\beta = 1/4$ và $\gamma = 1/2$, đảm bảo tính ổn định vô điều kiện cho hệ phương trình động lực học.
• Xác định giá trị ban đầu: Sử dụng ma trận độ cứng động học $S_\infty(\omega)$ tại thời điểm $t = 0^+$ thu được từ bài toán trị riêng của biên truyền sóng nhất quán (CTB).
• Mô hình cản vật liệu (damping): Giả định ma trận cản cục bộ dạng cản Rayleigh tỷ lệ với ma trận khối lượng và độ cứng: $[\mathbf{C}]^e = \alpha_0 [\mathbf{K}]^e + \alpha_1 [\mathbf{M}]^e$, với $\alpha_0 = \frac{2\omega_1 \omega_2}{\omega_1 + \omega_2} \xi$ và $\alpha_1 = \frac{2}{\omega_1 + \omega_2} \xi$, trong đó $\xi$ là tỷ số cản và $\omega_1, \omega_2$ là các tần số biên chặn vùng quan tâm.
• Kỹ thuật lớp phối hợp rời rạc (PMDLs): Rời rạc hóa các lớp mở rộng với tọa độ phức để triệt tiêu năng lượng phản xạ ở các biên sâu và biên góc trong bán không gian và toàn không gian.

Nội dung chính theo từng chương

Chapter 1: Introduction

Chương mở đầu hệ thống hóa vai trò của các điều kiện biên hấp thụ trong việc tính toán số đối với miền vô hạn trong kỹ thuật công trình, địa chấn, thủy động lực học, khí tượng học và âm học. Tác giả phân loại ranh giới hấp thụ thành nhóm cục bộ (local) và toàn cục (global/non-local). Chương 1 đặt tiền đề cho toàn bộ luận án: xây dựng điều kiện biên toàn cục trong miền thời gian dựa trên nhận định hình học rằng tính tương hỗ luôn tồn tại giữa hai miền vô hạn đồng nhất chỉ khác nhau bởi vị trí của đường biên phân cắt. Cấu trúc 5 bài báo khoa học tương ứng với các chương tiếp theo cũng được phác thảo chi tiết.

Chapter 2: Reciprocal absorbing boundary condition for the time-domain numerical analysis of wave motion in unbounded layered media

Chương này công bố công thức nền tảng của RABC áp dụng cho môi trường phân lớp. Tác giả xem xét hai hệ thống vô hạn $I$ và $E$ giống hệt nhau, phân chia ranh giới thành biên gần thẳng đứng ($\Gamma_0$), biên xa vô tận ($\Gamma_\infty$), bề mặt tự do ứng suất thuần nhất ($\Gamma_N$) và đáy ngàm cứng ($\Gamma_D$). Bằng cách triệt tiêu các tích phân biên trên $\Gamma_\infty$, $\Gamma_N$ và $\Gamma_D$, quan hệ tương hỗ rút gọn chỉ còn tích phân tích chập trên biên gần: $$\int_{\Gamma_0^I} \tau_j^I(\mathbf{x}, t) * a_j^E(\mathbf{x}, t) , d\Gamma = \int_{\Gamma_0^E} \tau_j^E(\mathbf{x}, t) * a_j^I(\mathbf{x}, t) , d\Gamma$$

Rời rạc hóa miền giữa $\Gamma_0^I$ và $\Gamma_0^E$ thành các phần tử hữu hạn với $n$ bậc tự do. Áp dụng $n$ trường hợp tải trọng xung bước nhảy đơn vị $\mathbf{P}_I^j = \mathbf{e}_j$ ($1 \le j \le n$) tại biên trong. Áp dụng định lý tương hỗ giữa bài toán thứ $j$ và bài toán thứ $l$ tại biên ngoài dẫn đến: $$(\mathbf{P}_I^j)^T * \mathbf{a}_E^l = (\mathbf{P}_E^l)^T * \mathbf{a}_I^j \quad (1 \le j, l \le n)$$

Thiết lập hệ phương trình đại số gồm $2n^2$ phương trình cân bằng phần tử hữu hạn và $n^2$ phương trình tương hỗ để giải chính xác $3n^2$ ẩn số gia tốc $\mathbf{a}_I, \mathbf{a}_E$ và lực biên $\mathbf{P}_E$.

Tác giả giải bài toán trị riêng để xác định ma trận độ cứng động học $S_\infty(\omega)$ nhằm tìm nghiệm ban đầu cho:

• Dầm bán vô hạn trên gối tựa đơn (nghiệm phương trình đặc trưng $s_2 X^2 + 2s_1 X + s_2 = 0$).
• Dầm bán vô hạn trên gối tựa đàn hồi liên tục với độ cứng lò xo $k_s$ (phương trình bậc hai theo $X = \cosh K$).
• Thanh bán vô hạn chịu lực dọc trục.
• Tầng đất phân lớp trong trạng thái biến dạng phẳng (giải bài toán trị riêng bậc hai $(k^2 \mathbf{A} + ik\mathbf{B} + \mathbf{G} - \omega^2 \mathbf{M})\boldsymbol{\Delta} = \mathbf{0}$ cho $4n$ trị riêng và chọn một nửa số mốt thỏa mãn điều kiện bức xạ sóng).
• Tầng đất phân lớp chịu cắt ngoài mặt phẳng (antiplane shear với phương trình $(k^2 \mathbf{A} + \mathbf{G} - \omega^2 \mathbf{M})\boldsymbol{\Delta} = \mathbf{0}$).

Kết quả tính toán kiểm chứng số cho thấy:

• Đối với dầm bán vô hạn và thanh bán vô hạn chịu lực bước nhảy đơn vị, độ xoay và chuyển vị tại nút 1, 2 tính bằng RABC trùng khít với nghiệm lưới mở rộng (extended mesh) 100 phần tử.
• Phân tích sai số theo chuẩn $L_2$ khẳng định sai số tiến về 0 khi bước thời gian $\Delta t \to 0$ với bậc hội tụ bậc hai (second-order rate of convergence), hoàn toàn tương thích với sơ đồ Newmark và quy tắc hình thang.
• Các thí nghiệm số với tải trọng xung tam giác đối xứng, xung chữ nhật và sóng con Ricker trên tầng đất biến dạng phẳng với các hệ số Poisson $\nu = 0.33, 0.40, 0.45$ đều thể hiện tính ổn định và triệt tiêu sóng phản xạ hoàn toàn.

Chapter 3: Application of reciprocal absorbing boundary condition to transient analysis of acoustic wave propagation

Chương 3 mở rộng RABC sang phân tích quá độ sóng âm tuyến tính trong ống dẫn sóng (waveguides). Phương trình sóng âm tuyến tính được biểu diễn qua thế vận tốc $\phi$: $$\nabla^2 \phi = \frac{1}{c^2} \ddot{\phi}$$ với vận tốc hạt $\mathbf{v} = \nabla \phi$ và áp suất âm học $p = -\rho \dot{\phi}$. Định lý tương hỗ âm học trong miền thời gian được thiết lập giữa hai trạng thái vận tốc âm và áp suất.

Luận án kiểm tra bài toán âm học 1D với nguồn bước nhảy Heaviside và nguồn Ricker wavelet, xác nhận tính ổn định và tốc độ hội tụ bậc hai. Trong bài toán 2D, tác giả mô hình hóa một hốc tròn (circular cavity) nằm trong một kênh dẫn vô hạn (infinite channel). Lưới phần tử hữu hạn mô phỏng trường sóng từ thời điểm $t = 5\text{ s}$ đến $t = 50\text{ s}$. Kết quả cho thấy các mặt sóng bức xạ truyền ra ngoài êm thuận tại các nút quan sát A, B, C mà không ghi nhận bất kỳ sự phản xạ nhân tạo nào dội ngược từ các biên cắt RABC.

Chapter 4: Reciprocal absorbing boundary condition with perfectly matched discrete layers for SH waves in a layered half-space

Chương 4 giải quyết bài toán truyền sóng cắt ngoài mặt phẳng (SH waves) trong bán không gian phân lớp. Do bán không gian không có đáy cứng cố định mà trải rộng xuống độ sâu vô hạn, việc chỉ đặt RABC ở biên đứng là chưa đủ. Tác giả phân chia bán không gian thành lớp bề mặt và bán không gian bên dưới, kết hợp RABC tại các biên đứng với các lớp rời rạc phối hợp hoàn hảo (PMDLs) tại đáy miền tính toán.

Ưu điểm nổi bật của cấu trúc lai ghép này là loại bỏ hoàn toàn việc phải thiết lập các phần tử hấp thụ góc (corner elements) phức tạp—vốn là điểm nghẽn toán học trong các mô hình biên hấp thụ truyền thống. Các thử nghiệm tính toán bao gồm:

• Tầng đất hai lớp cắt bởi RABC chịu lực bước nhảy tại nút 1 và nút 2.
• Tầng đồng nhất cắt bằng PMDLs chịu lực sóng con Ricker tại điểm A.
• Bán không gian đồng nhất và bán không gian phân lớp cắt bởi tổ hợp RABC và PMDLs.

Các ảnh chụp mặt sóng (snapshots) tại các thời điểm $t = 1\text{ s}, 2\text{ s}, 3\text{ s}, 4\text{ s}$ chứng minh năng lượng sóng SH đi qua biên đáy PMDLs và biên bên RABC một cách triệt để.

Chapter 5: Reciprocal absorbing boundary condition with perfectly matched discrete layers for SV-P waves in a layered half-space

Chương 5 mở rộng bài toán sang hệ sóng kết hợp phức tạp hơn nhiều: sóng dọc và sóng cắt trong mặt phẳng (SV-P waves) thuộc bài toán biến dạng phẳng trong bán không gian phân lớp. Tác giả dẫn xuất RABC theo ba cách tiếp cận khác nhau:

1. Dựa trên định lý tương hỗ chuyển vị truyền thống.
2. Dựa trên ma trận độ cứng động học.
3. Dựa trên ma trận đáp ứng xung đơn vị gia tốc (acceleration unit-impulse response matrix).

Hệ thống PMDLs cho sóng SV-P được cấu trúc thành: phần tử thường trong miền tính toán, phần tử cạnh $\text{PMDL}_x$ (cắt biên ngang), phần tử cạnh $\text{PMDL}z$ (cắt biên đứng) và phần tử góc $\text{PMDL}{xz}$ (xử lý vùng góc đáy), đi kèm quy tắc tích phân điểm Gauss đặc thù.

Phát hiện quan trọng của chương này là việc kiểm tra tính ổn định lâu dài (long-term stability):

• Khi chỉ sử dụng thuần túy PMDLs, phản ứng của tầng đất đồng nhất và bán không gian phân lớp xuất hiện hiện tượng mất ổn định (instability/drift) khi thời gian phân tích kéo dài.
• Khi kết hợp RABC cho biên đứng và một số lượng nhỏ các lớp PMDLs (chỉ cần 3 lớp PMDL) ở đáy, hệ thống duy trì được độ chính xác cao và triệt tiêu hoàn toàn sự mất ổn định trong phản ứng dài hạn.
• Đánh giá chỉ số điều kiện ma trận $\kappa(\mathbf{M})$ và sai số tương đối theo chuẩn $L_2$ khẳng định mô hình RABC-PMDL có tính ổn định số vượt trội so với các phương pháp biên hấp thụ cục bộ hoặc PML truyền thống.

Chapter 6: Transient analysis of full-space unbounded domains by reciprocal absorbing boundaries

Chương 6 mở rộng phương pháp RABC kết hợp PMDLs cho miền toàn không gian hai chiều vô hạn (full-space) ngoại vi của một hốc/đĩa tròn bán kính $R$. Miền tính toán được chia tách và áp dụng biên hấp thụ cho 4 chế độ dao động riêng biệt của đĩa tròn:

1. Dao động dọc trục (axial / out-of-plane / vertical direction).
2. Dao động xuyên tâm (radial direction).
3. Dao động xoắn (torsional / circumferential direction).
4. Dao động ngang trong mặt phẳng (horizontal / lateral direction).

Tác giả thiết lập bộ thông số vật liệu quy chuẩn: khối lượng riêng $\rho = 1$, mô đun trượt $\mu = 1$, hệ số Poisson $\nu = 1/3$. Bộ tham số này cho vận tốc sóng cắt $c_s = \sqrt{\mu/\rho} = 1\text{ m/s}$ và vận tốc sóng nén $c_d = \sqrt{(\lambda + 2\mu)/\rho} = 2\text{ m/s}$.

Tác giả kiểm tra đáp ứng quá độ dưới tác dụng của tải trọng sóng con Ricker với hai chu kỳ đặc trưng $t_s = 3/\pi$ và $t_s = 3/2\pi$, đồng thời khảo sát phản ứng dao động điều hòa xác lập ở các tần số $f = 0.01\text{ Hz}, 1\text{ Hz}, 5\text{ Hz}$ (tần số không thứ nguyên $a_0 = \pi$).

Bảng sai số tương đối chuẩn $L_2$ tại các điểm quan sát mẫu giữa nghiệm RABC-PMDLs và nghiệm lưới mở rộng (extended mesh) khẳng định độ chính xác cao của phương pháp trên cả 4 chế độ dao động:

Chế độ dao động	Tải trọng kiểm tra	Vận tốc sóng liên quan	Phạm vi sai số tương đối chuẩn $L_2$
Dọc trục (Vertical/Axial)	Ricker wavelet ($t_s = 3/\pi, 3/2\pi$) & Điều hòa ($f = 0.01, 2.005\text{ s}$)	Sóng cắt ($c_s = 1\text{ m/s}$)	Nhỏ, hội tụ hoàn toàn với lưới mở rộng
Xuyên tâm (Radial)	Ricker wavelet & Điều hòa ($f = 0.01\text{ Hz}, 1\text{ Hz}, 5\text{ Hz}$)	Sóng nén ($c_d = 2\text{ m/s}$)	Dưới mức sai số số trị thông thường
Xoắn (Torsional)	Ricker wavelet & Điều hòa ($f = 0.01\text{ Hz}, 1\text{ Hz}, 5\text{ Hz}$)	Sóng cắt thuần túy ($c_s = 1\text{ m/s}$)	Trùng khít nghiệm giải tích/lưới lớn
Ngang (Horizontal)	Ricker wavelet & Điều hòa ($f = 0.01\text{ Hz}, 1\text{ Hz}, 5\text{ Hz}$)	Kết hợp sóng nén và sóng cắt	Độ chính xác cao, không phản xạ biên

Kết quả và những đóng góp mới

Luận án tự xác định các đóng góp mới về mặt lý luận khoa học và thực tiễn kỹ thuật tính toán:

Về mặt lý luận:

• Thiết lập thành công dạng thức toán học của Định lý tương hỗ dựa trên gia tốc (acceleration-based reciprocity theorem) trực tiếp trong miền thời gian cho các bài toán đàn hồi động học và âm học.
• Đưa ra nguyên lý biên hấp thụ tương hỗ: chứng minh rằng tính tương hỗ tồn tại giữa hai miền vô hạn đồng nhất có biên tịnh tiến cho phép trích xuất đáp ứng động lực học của miền ngoài mà không cần sử dụng hàm Green giải tích.

Về mặt thực tiễn tính toán:

• Xây dựng một phương pháp biên hấp thụ toàn cục trực tiếp trong miền thời gian (RABC) có tốc độ hội tụ bậc hai ổn định, liên kết hoàn toàn tự nhiên với thuật toán Newmark và phương pháp phần tử hữu hạn.
• Đề xuất giải pháp lai ghép RABC với PMDLs: RABC chịu trách nhiệm cắt biên theo phương ngang (biên đứng), trong khi PMDLs xử lý biên đáy vô hạn. Mô hình lai này khắc phục triệt để hiện tượng mất ổn định dài hạn (long-term instability) thường gặp của các lớp PML/PMDL thuần túy, đồng thời loại bỏ sự phụ thuộc vào các phần tử góc hấp thụ phức tạp.
• Cung cấp giải pháp tính toán hiệu quả cho các bài toán tương tác đất - kết cấu trong miền thời gian, mở đường cho việc phân tích các bài toán có tính phi tuyến cục bộ trong miền gần.

Về kiến nghị giải pháp:

• Đối với các bài toán nền phân lớp trên bán không gian vô hạn, tác giả khuyến nghị sử dụng mô hình kết hợp RABC ở biên bên và tối thiểu 3 lớp rời rạc PMDL ở biên đáy để đạt sự tối ưu giữa chi phí bộ nhớ tính toán và độ ổn định lâu dài của nghiệm số.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu được xác định từ phạm vi giả định trong văn bản:

• Phương pháp RABC dựa trên giả định hai miền vô hạn $I$ và $E$ phải hoàn toàn đồng nhất về mặt hình học và đặc tính cơ lý vật liệu theo phương truyền sóng tịnh tiến. Do đó, phương pháp chưa áp dụng trực tiếp được cho các môi trường có sự thay đổi hình học lớp bất quy tắc hoặc vật liệu biến thiên ngẫu nhiên theo phương ngang của miền xa.
• Quá trình khởi tạo tính toán tại thời điểm $t = 0^+$ vẫn cần sử dụng ma trận độ cứng động học $S_\infty(\omega)$ suy biến từ việc giải bài toán trị riêng của miền tần số.

Hướng nghiên cứu tiếp theo:

• Mở rộng thuật toán RABC kết hợp PMDLs cho các bài toán không gian ba chiều (3D) hoàn chỉnh với hình học móng và kết cấu phức tạp.
• Ứng dụng RABC để khảo sát các bài toán có phi tuyến vật liệu mạnh ở miền gần, ví dụ như sự hình thành vùng dẻo trong nền đất hoặc hiện tượng tách - trượt tiếp xúc giữa móng công trình và môi trường đất dưới tác động của động đất mạnh.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các nhóm đối tượng:

• Nghiên cứu sinh và nhà nghiên cứu cơ học tính toán: Cung cấp phương pháp luận chặt chẽ về việc khai thác định lý tương hỗ động học đàn hồi và kỹ thuật tích phân tích chập số trị trong miền thời gian.
• Kỹ sư và chuyên gia địa kỹ thuật, động lực học công trình và kỹ thuật địa chấn: Hữu ích cho việc phát triển các mô hình tương tác động lực học giữa công trình và nền đất vô hạn (đặc biệt là công trình ngầm, cọc, móng máy và đập bê tông trọng lực).
• Kỹ sư phát triển phần mềm mô phỏng phần tử hữu hạn: Bản mô tả thuật toán ở Chương 2 (cách thức thiết lập và giải hệ phương trình $3n^2$ ẩn số) và kỹ thuật ghép nối phần tử giãn nở phức PMDLs ở Chương 4 và Chương 5 cung cấp các chỉ dẫn cụ thể để lập trình mô đun biên hấp thụ miền thời gian.

Câu hỏi thường gặp

1. Điểm khác biệt cốt lõi giữa RABC trong luận án này với phương pháp phương trình tích phân biên (BIEM) truyền thống là gì? Cả hai phương pháp đều có thể xuất phát từ định lý tương hỗ, nhưng BIEM đòi hỏi phải biết trước nghiệm cơ bản (hàm Green) giải tích trong miền thời gian—điều rất khó thiết lập cho môi trường phân lớp. Trong khi đó, RABC tận dụng tính tương hỗ giữa hai miền giống hệt nhau chỉ khác vị trí biên, qua đó hàm Green tự động được xác định trong quá trình tính toán số mà không cần nghiệm giải tích trước.

2. Thuật toán tích phân thời gian và quy tắc tính tích phân tích chập nào được sử dụng trong luận án? Luận án sử dụng phương pháp Newmark với sơ đồ gia tốc trung bình không đổi ($\beta = 1/4, \gamma = 1/2$) để tích bước thời gian cho phương trình chuyển động, và sử dụng quy tắc hình thang (trapezoidal rule) để tính toán các tích phân tích chập thời gian, mang lại tốc độ hội tụ bậc hai theo chuẩn sai số $L_2$.

3. Tại sao tác giả lại kết hợp RABC với PMDLs trong các bài toán bán không gian và toàn không gian? RABC được xây dựng tối ưu cho các biên đứng của tầng phân lớp nhưng gặp khó khăn khi xử lý biên đáy mở rộng vô hạn. PMDLs xử lý tốt biên đáy nhưng dễ bị mất ổn định dài hạn và đòi hỏi phần tử góc phức tạp. Sự kết hợp giữa RABC ở biên đứng và PMDLs ở biên đáy giúp loại bỏ hoàn toàn phần tử góc và khắc phục triệt để hiện tượng trôi mất ổn định theo thời gian của PMDL thuần túy.

4. Trong bài toán toàn không gian ở Chương 6, các thông số vật liệu chuẩn được thiết lập ra sao? Tác giả thiết lập các thông số cơ lý chuẩn hóa gồm: khối lượng riêng $\rho = 1$, mô đun trượt $\mu = 1$, và hệ số Poisson $\nu = 1/3$. Hệ thông số này cho vận tốc sóng cắt $c_s = 1\text{ m/s}$ và vận tốc sóng nén $c_d = 2\text{ m/s}$.

5. Ma trận cản vật liệu của đất được tính toán theo mô hình nào? Ma trận cản vật liệu được mô hình hóa theo dạng cản Rayleigh cục bộ tỷ lệ với ma trận khối lượng và độ cứng: $[\mathbf{C}]^e = \alpha_0 [\mathbf{K}]^e + \alpha_1 [\mathbf{M}]^e$. Các hệ số $\alpha_0, \alpha_1$ được xác định dựa trên tỷ số cản quy định $\xi$ cùng hai tần số $\omega_1, \omega_2$ bao quanh dải tần số quan tâm.

Kết luận

Luận án của Cuong Tan Nguyen đã giải quyết thành công bài toán biên hấp thụ toàn cục trong miền thời gian cho môi trường vô hạn thông qua việc đề xuất Điều kiện biên hấp thụ tương hỗ (RABC) dựa trên định lý tương hỗ gia tốc. Phương pháp loại bỏ sự phụ thuộc vào hàm Green giải tích, đồng thời khi kết hợp cùng kỹ thuật PMDLs đã xử lý hiệu quả hiện tượng sóng bức xạ trong các môi trường phức tạp từ tầng đất phân lớp, bán không gian đến toàn không gian hai chiều. Các kết quả kiểm chứng số khẳng định RABC đạt độ chính xác cao, hội tụ bậc hai và duy trì tính ổn định số bền vững trong các phân tích động lực học quá độ dài hạn.