Mô Hình Giảm Bậc cho Phân Tích Tương Tác Động Lực Đất-Ống Dẫn

Luận án tiến sĩ về mô hình giảm bậc phân tích tương tác động lực đất ống dẫn. Nghiên cứu chuyên sâu về kỹ thuật và ứng dụng trong xây dựng công trình.

• Tác giả: Kien Trung Nguyen
• Người hướng dẫn: Professor Domniki Asimaki
• Trường học: California Institute Of Technology
• Chuyên ngành: Geotechnical Earthquake Engineering
• Đề tài: Reduced-order model for dynamic soil-pipe interaction analysis
• Loại tài liệu: Thesis
• Năm xuất bản: 2020
• Địa điểm: Pasadena, California

Mục lục chi tiết

• Acknowledgements
• Published Content and Contributions
• Table of Contents
• List of Illustrations
• List of Tables
• 1. CHƯƠNG 1: Introduction
• 1.1. Pipelines and seismic actions
• 1.2. Methods for soil-pipe interaction analysis
• 1.2.1. Model neglecting soil-pipe interaction
• 1.2.2. Beam-on-Winkler-foundation model considering soil-pipe interaction
• 1.3. Full three-dimensional model considering soil-pipe interaction
• 1.4. Challenges in soil-pipe interaction analysis
• 1.5. Organization of the text
• 2. CHƯƠNG 2: Smooth nonlinear hysteresis model for coupled biaxial soil-pipe interaction in sandy soils
• 2.1. Uniaxial hysteresis model
• 2.2. Biaxial hysteresis model
• 2.3. Finite element method
• 2.4. Smoothed-particle hydrodynamics
• 2.5. Validation of the FEM and SPH models
• 2.6. Parameter calibration for BMBW model
• 2.7. Uniaxial cyclic loading
• 2.8. Suggestions for input parameters
• 3. CHƯƠNG 3: Dynamic axial soil impedance function for rigid circular structures buried in elastic half-space
• 3.1. Review of axial soil impedance function
• 3.2. Analytical solution for soil impedance function of homogeneous half-space
• 3.3. Finite element analysis for soil impedance functions of homogeneous and two-layered half-spaces
• 3.3.1. Numerical computation of impedance function
• 3.3.2. Finite element models
• 3.4. Homogeneous half-space
• 3.5. Two-layered half-space
• 3.5.1. Effect of material contrast
• 3.5.2. Effect of structure location
• 4. CHƯƠNG 4: Dynamic in-plane soil impedance functions for rigid circular structures buried in elastic half-space
• 4.1. Review of in-plane soil impedance functions
• 4.2. Analytical solution for soil impedance functions of homogeneous half-space
• 4.3. Traction-free condition at 𝑦 = 0
• 4.4. Dirichlet boundary condition at cylinder interface
• 4.5. Calculation of in-plane soil impedance functions
• 4.6. Direct evaluation of the integral
• 4.7. Finite element analysis for soil impedance functions of homogeneous and two-layered half-spaces
• 4.7.1. Numerical computation of in-plane soil impedance functions
• 4.7.2. Finite element models
• 4.8. Homogeneous half-space
• 4.8.1. Effect of burial depth
• 4.8.2. Effect of Poisson’s ratio
• 4.9. Two-layered half-space
• 4.9.1. Effect of material contrast
• 4.9.2. Effect of structure location
• 5. CHƯƠNG 5: Application: Reduced-order modeling of buried pipe subjected to the propagation of Rayleigh surface wave
• 5.1. Models for soil-pipe interaction analysis
• 5.1.1. Model neglecting soil-pipe interaction
• 5.1.2. Model considering soil-pipe interaction with free-field input
• 5.1.3. Models based on substructure and finite element methods
• 5.2. Results and comparisons
• 5.3. Summary of previous chapters
• Appendix A: Asymptotic method for computing high oscillatory integrals

Tóm tắt

I. Tổng Quan Mô Hình Giảm Bậc Phân Tích Đất Ống Dẫn

Mạng lưới đường ống chôn ngầm đóng vai trò quan trọng trong việc vận chuyển nước, khí đốt tự nhiên, nhiên liệu và dầu mỏ. Chúng là huyết mạch của xã hội hiện đại. Theo CIA (2018), tổng chiều dài đường ống trên toàn cầu khoảng 3.500.000 km. Sự gia tăng nhanh chóng nhu cầu toàn cầu về nước và năng lượng đã đặt ra yêu cầu cao về lắp đặt, vận hành và khả năng phục hồi của các mạng lưới truyền tải. Do đó, việc phân tích tương tác động lực đất-ống dẫn trở nên vô cùng quan trọng. Các phương pháp hiện tại còn nhiều hạn chế, đòi hỏi các mô hình hiệu quả hơn, đặc biệt là các mô hình giảm bậc. Các nghiên cứu trước đây cho thấy rằng biến dạng nền đất do địa chất có thể gây ra rủi ro đáng kể cho các hệ thống đường ống. Các hướng dẫn của ASCE (1984), PRCI (2004), ALA (2005) và CEN (2006) chỉ ra hai loại nguy cơ động đất chính liên quan đến tính toàn vẹn của đường ống: biến dạng nền đất tạm thời (TGD) và biến dạng nền đất vĩnh viễn (PGD).

1.1. Tầm quan trọng của phân tích tương tác đất ống SPI

Đánh giá chính xác cách đất truyền tải trọng lên ống dẫn, và ngược lại, là yếu tố then chốt trong thiết kế các hệ thống đường ống có khả năng phục hồi tốt. Dữ liệu quan sát sau động đất cho thấy thiệt hại cho đường ống chủ yếu do biến dạng nền đất vĩnh viễn (PGD) gây ra, tập trung ở các khu vực riêng biệt với tỷ lệ thiệt hại cao. Một cách tiếp cận đúng đắn giúp giảm thiểu rủi ro và đảm bảo an toàn cho cộng đồng.

1.2. Các loại tải trọng động tác động lên ống dẫn chôn ngầm

Các nguy cơ động đất chính ảnh hưởng đến tính toàn vẹn của ống dẫn bao gồm biến dạng nền đất tạm thời (TGD), tức là rung lắc mặt đất do sóng địa chấn, và biến dạng nền đất vĩnh viễn (PGD), là các hư hỏng mặt đất do đứt gãy, lan truyền bên, trượt lở và chuyển động của sườn dốc. TGD có thể gây ra biến dạng không mong muốn trong mạng lưới đường ống, đặc biệt ở các môi trường đất không đồng nhất. Cần xem xét cả hai loại tải trọng này trong quá trình phân tích.

II. Thách Thức Phân Tích Tương Tác Động Lực Đất Ống Dẫn

Phân tích tương tác động lực đất-ống dẫn (SPI) đặt ra nhiều thách thức đáng kể. Các phương pháp truyền thống thường bỏ qua sự phức tạp của hành vi phi tuyến của đất và sự tương tác giữa đất và ống. Mô hình hóa chính xác ứng xử động lực của hệ thống đòi hỏi các phương pháp tính toán phức tạp và tốn kém. Việc kiểm chứng các mô hình bằng dữ liệu thực nghiệm cũng rất khó khăn. Theo luận án, các phương pháp phân tích tương tác đất-ống có thể được phân loại thành ba loại chính: mô hình bỏ qua tương tác đất-ống, mô hình dầm trên nền Winkler, và mô hình ba chiều đầy đủ. Mỗi phương pháp đều có ưu và nhược điểm riêng, và việc lựa chọn phương pháp phù hợp phụ thuộc vào độ chính xác và chi phí tính toán mong muốn.

2.1. Hạn chế của mô hình dầm trên nền Winkler truyền thống

Mặc dù đơn giản và hiệu quả về mặt tính toán, mô hình dầm trên nền Winkler dựa trên giả định rằng phản ứng của đất là tuyến tính và cục bộ. Điều này không chính xác khi đất thể hiện hành vi phi tuyến và tương tác trên một vùng rộng hơn. Mô hình này cũng bỏ qua ảnh hưởng của các yếu tố như độ sâu chôn lấp và loại đất. Cần các cải tiến để mô hình hóa chính xác hơn.

2.2. Độ phức tạp của mô hình ba chiều và chi phí tính toán cao

Mô hình ba chiều đầy đủ có khả năng mô phỏng tương tác đất-ống một cách chi tiết hơn, nhưng đòi hỏi chi phí tính toán rất lớn. Việc xây dựng và giải các mô hình này tốn nhiều thời gian và tài nguyên, khiến chúng không phù hợp cho các ứng dụng thực tế. Ngoài ra, việc thu thập dữ liệu đầu vào chính xác cho các mô hình này cũng là một thách thức lớn.

2.3. Yêu cầu về kiểm chứng mô hình bằng dữ liệu thực nghiệm

Việc kiểm chứng các mô hình tương tác đất-ống bằng dữ liệu thực nghiệm là rất quan trọng để đảm bảo độ tin cậy của chúng. Tuy nhiên, việc thực hiện các thí nghiệm trong phòng thí nghiệm hoặc ngoài hiện trường tốn kém và khó khăn. Cần có các phương pháp kiểm chứng hiệu quả hơn và dễ tiếp cận hơn.

III. Giải Pháp Mô Hình Giảm Bậc ROM cho Bài Toán Đất Ống

Để giải quyết các thách thức trên, luận án này đề xuất sử dụng mô hình giảm bậc (ROM) để phân tích tương tác động lực đất-ống dẫn. ROM là một kỹ thuật giảm số lượng bậc tự do trong mô hình, giúp giảm đáng kể chi phí tính toán mà vẫn duy trì được độ chính xác chấp nhận được. Nghiên cứu này trình bày việc phát triển một macroelement đơn giản hóa để nắm bắt SPI trong đất rời rạc chịu tải trọng tùy ý vuông góc với trục ống. Mô hình giảm bậc cung cấp một sự cân bằng tốt giữa độ chính xác và hiệu quả tính toán, cho phép phân tích hệ thống đất-ống phức tạp một cách hiệu quả.

3.1. Ưu điểm của mô hình giảm bậc POD và ROM

Mô hình giảm bậc POD (Proper Orthogonal Decomposition) và ROM (Reduced Order Model) cho phép giảm kích thước của bài toán bằng cách xác định các chế độ quan trọng nhất của biến dạng. Điều này dẫn đến giảm đáng kể thời gian tính toán mà vẫn duy trì được độ chính xác cao. Các mô hình giảm bậc đặc biệt hữu ích cho các bài toán phân tích động lực học phức tạp.

3.2. Ứng dụng phương pháp phần tử hữu hạn và phương pháp số

Nghiên cứu này sử dụng phương pháp phần tử hữu hạn (FEM) để tạo ra dữ liệu huấn luyện cho việc xây dựng ROM. FEM cho phép mô phỏng chính xác hành vi của đất và ống dưới các điều kiện tải trọng khác nhau. Các phương pháp số khác cũng được sử dụng để giải các phương trình chi phối của mô hình giảm bậc.

3.3. Phát triển mô hình phần tử vĩ mô cho tương tác đất ống

Luận án đề xuất một macroelement đơn giản hóa được thiết kế để nắm bắt tương tác đất-ống trong đất rời rạc. Macroelement này sử dụng một mô hình trễ đơn trục có thể nắm bắt các đường cong phản lực-chuyển vị phi tuyến mượt mà của các bài toán SPI. Bằng cách sử dụng bộ lọc Kalman không dấu, tham số mô hình κ kiểm soát độ mượt mà của vùng chuyển tiếp từ tuyến tính sang dẻo được suy ra bằng cách sử dụng dữ liệu thực nghiệm đã được công bố.

IV. Mô Hình Trễ Phi Tuyến Cách Tiếp Cận Mới cho SPI

Luận án này giới thiệu một mô hình trễ phi tuyến mới để mô tả tương tác đất-ống dẫn trong đất cát. Mô hình này có khả năng mô phỏng chính xác hành vi phi tuyến và phụ thuộc vào lịch sử tải trọng của đất. Bằng cách sử dụng bộ lọc Kalman không dấu, tham số mô hình κ kiểm soát độ mượt mà của vùng chuyển tiếp từ tuyến tính sang dẻo được suy ra bằng cách sử dụng dữ liệu thực nghiệm đã được công bố. Mô hình mới được kiểm chứng bằng các kết quả thí nghiệm và mô phỏng số, cho thấy độ chính xác cao.

4.1. Mô hình đất đàn hồi dẻo và các ứng dụng

Việc sử dụng mô hình đất đàn hồi-dẻo, như mô hình Mohr-Coulomb, mô hình Hardening Soil và mô hình Cam Clay, cho phép mô phỏng chính xác hơn hành vi của đất dưới tải trọng động. Các mô hình này tính đến các yếu tố như độ bền cắt, góc ma sát trong và độ chặt của đất.

4.2. Hiệu chỉnh tham số mô hình bằng bộ lọc Kalman không dấu

Bộ lọc Kalman không dấu (UKF) là một phương pháp hiệu quả để ước lượng các tham số của mô hình trễ phi tuyến từ dữ liệu thực nghiệm. UKF cho phép tính đến các yếu tố như nhiễu và độ không đảm bảo trong dữ liệu, giúp cải thiện độ chính xác của mô hình.

4.3. Mở rộng mô hình đơn trục sang hiệu ứng tải trọng hai trục

Mô hình đơn trục được mở rộng sang các hiệu ứng tải trọng hai trục và macroelement có thể nắm bắt các hiệu ứng như thắt nút và ghép nối giãn nở cắt. Các thông số đầu vào của mô hình được hiệu chỉnh bằng cách sử dụng các phân tích phần tử hữu hạn (FE) được xác thực bằng các thí nghiệm. Các FDC của mô hình hai trục được xác minh bằng cách so sánh với các mô phỏng FE và thủy động lực học hạt làm mịn (SPH) cho các mẫu tải khác nhau: đơn trục tuần hoàn, hình chữ 0, hình chữ 8 và tải trọng tạm thời.

V. Ứng Dụng Phân Tích Ống Chôn Ngầm Chịu Sóng Rayleigh

Luận án trình bày một ứng dụng thực tế của ROM để phân tích ứng xử của ống chôn ngầm chịu tác dụng của sóng Rayleigh. Kết quả cho thấy rằng ROM có thể mô phỏng chính xác phản ứng của ống với chi phí tính toán thấp hơn nhiều so với phân tích phần tử hữu hạn trực tiếp. Theo luận án, có thể kết hợp các SIF thu được trước đó vào một mô hình giảm bậc để phân tích các bài toán SPI, cụ thể là một ống thẳng chịu sóng bề mặt Rayleigh lan truyền qua các bán không gian đàn hồi đồng nhất và không đồng nhất. Lịch sử thời gian dịch chuyển được tính toán tại các điểm kiểm soát được hiển thị là phù hợp với các lịch sử được tính toán bằng các phân tích phần tử hữu hạn hai chiều trực tiếp.

5.1. Mô phỏng ống chịu tác dụng của sóng Rayleigh

Sóng Rayleigh là một loại sóng địa chấn lan truyền trên bề mặt đất. Việc mô phỏng ứng xử của ống chôn ngầm dưới tác dụng của sóng Rayleigh là rất quan trọng để đánh giá rủi ro động đất cho các hệ thống đường ống. ROM cho phép mô phỏng chính xác sự tương tác giữa sóng Rayleigh và ống, giúp đưa ra các quyết định thiết kế tốt hơn.

5.2. So sánh kết quả ROM với phân tích phần tử hữu hạn trực tiếp

Kết quả thu được từ ROM được so sánh với kết quả từ phân tích phần tử hữu hạn (FEA) trực tiếp để xác minh độ chính xác của mô hình. Kết quả so sánh cho thấy rằng ROM có thể mô phỏng chính xác phản ứng của ống với chi phí tính toán thấp hơn nhiều so với FEA trực tiếp.

5.3. Áp dụng hàm trở kháng đất động lực vào mô hình giảm bậc

Luận án trình bày các giải pháp phân tích để tính toán các hàm trở kháng đất dọc trục và trong mặt phẳng miền tần số (SIF) cho một cấu trúc tròn cứng vô hạn dài được chôn ngang trong bán không gian đàn hồi đồng nhất. Bằng cách sử dụng khai triển chuỗi Hankel- và Bessel-Fourier, một bài toán giá trị biên hỗn hợp đã được giải quyết, xem xét dịch chuyển điều hòa tại ranh giới cấu trúc và điều kiện biên không kéo tại bề mặt tự do bán không gian. Chúng tôi sau đó đã xác minh các giải pháp phân tích của mình bằng cách sử dụng các kết quả thu được từ các mô phỏng FE.

VI. Kết Luận Hướng Phát Triển Mô Hình Giảm Bậc Tương Lai

Luận án này đã trình bày một phương pháp hiệu quả để phân tích tương tác động lực đất-ống dẫn bằng cách sử dụng mô hình giảm bậc. Phương pháp này có tiềm năng ứng dụng rộng rãi trong thiết kế và đánh giá độ an toàn của các hệ thống đường ống. Nghiên cứu trong tương lai có thể tập trung vào việc mở rộng mô hình để bao gồm các loại đất phức tạp hơn và các điều kiện tải trọng khác nhau. Ngoài ra, việc phát triển các phương pháp kiểm chứng mô hình hiệu quả hơn cũng là một hướng đi quan trọng. Nghiên cứu của Asimaki, D. Seylabi (2019), cho thấy những mô hình giảm bậc thế hệ tiếp theo sẽ rất quan trọng trong tương lai.

6.1. Nghiên cứu ảnh hưởng của các loại đất và tải trọng khác nhau

Nghiên cứu trong tương lai nên tập trung vào việc mở rộng ROM để bao gồm các loại đất phức tạp hơn, chẳng hạn như đất sét và đất nhiều lớp. Ngoài ra, việc nghiên cứu ảnh hưởng của các điều kiện tải trọng khác nhau, chẳng hạn như động đất và lún đất, cũng là rất quan trọng.

6.2. Phát triển các phương pháp kiểm chứng mô hình hiệu quả

Việc phát triển các phương pháp kiểm chứng ROM hiệu quả hơn là rất quan trọng để đảm bảo độ tin cậy của mô hình. Các phương pháp này có thể bao gồm việc sử dụng dữ liệu thực nghiệm từ các thí nghiệm trong phòng thí nghiệm hoặc ngoài hiện trường, cũng như so sánh kết quả ROM với kết quả từ các mô hình phức tạp hơn.

6.3. Tối ưu hóa thiết kế đường ống chôn ngầm bằng mô hình giảm bậc

ROM có thể được sử dụng để tối ưu hóa thiết kế đường ống chôn ngầm, giúp giảm chi phí xây dựng và tăng độ an toàn. Bằng cách sử dụng ROM để mô phỏng ứng xử của ống dưới các điều kiện tải trọng khác nhau, các kỹ sư có thể thiết kế các hệ thống đường ống có khả năng chịu được các tác động bất lợi từ môi trường.

Nội dung chính

Tổng quan nghiên cứu

Mạng lưới đường ống chôn ngầm đóng vai trò là cơ sở hạ tầng huyết mạch trong việc vận chuyển nước, khí đốt và dầu mỏ trên phạm vi toàn cầu với tổng chiều dài vượt hơn 3.500.000 km theo thống kê quốc tế, trong đó riêng năm 2018 đã có thêm khoảng 24.000 km đường ống năng lượng mới được đưa vào vận hành. Tuy nhiên, các tuyến ống này thường kéo dài qua nhiều vùng địa chất phức tạp, đối mặt với hiểm họa địa chấn nghiêm trọng từ biến dạng nền đất vĩnh cửu (Permanent Ground Deformation - PGD) do nứt vỡ đứt gãy, trượt lở đất và biến dạng động tức thời (Transient Ground Deformation - TGD) do sóng địa chấn lan truyền. Cơ chế truyền tải trọng tương hỗ giữa môi trường đất và kết cấu ống ngầm (Soil-Pipe Interaction - SPI) là yếu tố then chốt quyết định độ an toàn và khả năng chống chịu của hệ thống.

Vấn đề nghiên cứu trọng tâm xuất phát từ việc các phương pháp thiết kế truyền thống, điển hình là hướng dẫn của Hiệp hội Kỹ sư Xây dựng Hoa Kỳ (ASCE 1984), chủ yếu dựa trên mô hình dầm trên nền Winkler phi tuyến (BNWF) với các lò xo đàn dẻo lý tưởng đơn trục độc lập. Cách tiếp cận này bỏ qua tính phi tuyến trơn thực tế, hiện tượng trễ chu kỳ, hiệu ứng thắt hẹp (pinching) khi hình thành khe hở đáy ống dưới tác dụng trọng trường và đặc biệt là sự ghép nối lực phản kháng hai trục ngang - đứng. Luận văn đặt mục tiêu phát triển mô hình bậc giảm (Reduced-Order Model - ROM) dựa trên phần tử vĩ mô hai trục cải tiến (BMBW) và nghiệm giải tích hàm trở kháng động lực học của đất (Soil Impedance Functions - SIF) cho ống ngầm trong môi trường nửa không gian đàn hồi và phân lớp.

Phạm vi nghiên cứu bao quát môi trường đất cát không dính với độ chặt tương đối từ 0% đến 100%, tỷ số chiều sâu chôn lấp trên đường kính ống H/D từ 1,5 đến 11,0, chịu tác động của tải trọng chu kỳ đa hướng và sóng địa chấn thực tế như trận động đất Kobe 1995 có độ lớn 6,9 Mw. Mô hình bậc giảm đề xuất mang ý nghĩa thực tiễn to lớn khi giúp cắt giảm hơn 90% thời gian tính toán so với phân tích phần tử hữu hạn 3D toàn phần, đồng thời duy trì độ chính xác cao với sai số phản lực dưới 5% so với mô phỏng hạt mịn SPH và dữ liệu thực nghiệm.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng tích hợp của bốn khung lý thuyết cốt lõi:

• Lý thuyết tương tác dầm trên nền Winkler phi tuyến (BNWF) mô tả phản lực đất thông qua hệ thống lò xo và giảm chấn phân bố dọc trục ống.
• Mô hình trễ phi tuyến Bouc-Wen biến tính (Modified Bouc-Wen - MBW) cho phép mô phỏng đường cong quan hệ lực - chuyển vị (FDC) có độ chuyển tiếp trơn tru giữa miền đàn hồi và miền dẻo.
• Lý thuyết truyền sóng trong môi trường đàn hồi nửa không gian vô hạn, sử dụng phương pháp khai triển chuỗi trực giao Hankel và Bessel-Fourier để giải bài toán biên hỗn hợp động lực học.
• Lý thuyết cơ học đất môi trường hạt rời tuân theo tiêu chuẩn phá hoại Mohr-Coulomb và cơ chế biến dạng dãn nở thể tích (shear-dilation).

Năm khái niệm chuyên ngành then chốt định hình cấu trúc mô hình gồm:

• Đường cong lực phản kháng - chuyển vị tương đối (FDC): Quan hệ động học đặc trưng giữa lực tương tác đất tác dụng lên ống và chuyển vị tương đối của ống.
• Hiệu ứng thắt hẹp (pinching effect): Hiện tượng sụt giảm độ cứng kháng lực về gần bằng 0 trong giai đoạn dỡ tải do hình thành khe hở cục bộ bên dưới đáy ống dưới tác dụng trọng trường.
• Ghép nối phản lực hai trục (biaxial coupling): Sự phụ thuộc lẫn nhau giữa phản lực cắt phương ngang và phản lực phương thẳng đứng khi ống chuyển vị xiên.
• Hàm trở kháng động lực học của đất (SIF): Đại lượng phức biểu diễn độ cứng động và hệ số cản nhớt phụ thuộc vào tần số kích động của môi trường đất.
• Cơ chế chuyển tiếp phá hoại nêm bị động (passive-wedge) sang cày xuyên (plow-through): Sự biến đổi hình thái trượt của khối đất quanh ống khi tỷ số chôn sâu vượt qua ngưỡng tới hạn H/D từ 6 đến 8.

Mô hình BMBW mở rộng công thức toán học bằng cách đưa vào ma trận độ cứng tiếp tuyến trực giao, tích hợp phần tử trượt khóa (slip-lock element) mắc nối tiếp với lò xo trễ phi tuyến, đồng thời hiệu chỉnh góc lệch bất đẳng hướng ứng suất địa tĩnh tự nhiên chi với sai số xấp xỉ chỉ khoảng 4% ở góc lệch 20 độ.

Phương pháp nghiên cứu

Nghiên cứu kết hợp chặt chẽ giữa nghiệm giải tích lý thuyết, thuật toán nhận dạng tham số tối ưu và mô phỏng số học nâng cao:

• Nguồn dữ liệu thực nghiệm: Sử dụng cơ sở dữ liệu đo đạc thực tế từ các thí nghiệm kéo ngang và nhổ ống quy mô lớn của Trautmann (1983) trên cát Cornell với đường kính ống 102 mm, thí nghiệm mô hình ly tâm của Finch (1999) và các công bố thực nghiệm của Robert (2016). Dữ liệu kích động động đất sử dụng bản ghi sóng gia tốc trạm OSAJ từ trận động đất Kobe 1995 cách tâm chấn 8,5 km.
• Cỡ mẫu và phương pháp chọn mẫu: Nghiên cứu thu thập và xử lý hơn 30 tập đường cong quan hệ lực - chuyển vị từ các thí nghiệm vật lý chuẩn, bao phủ 3 nhóm trạng thái đất: cát xốp (độ chặt 0% đến 35%), cát chặt vừa (độ chặt 35% đến 65%) và cát chặt (độ chặt 65% đến 100%) với tỷ số chôn lấp H/D biến thiên liên tục từ 1,5 đến 11,0. Phương pháp chọn mẫu mục tiêu (purposive sampling) được áp dụng nhằm thu thập đầy đủ các dạng phá hoại cơ học điển hình.
• Phương pháp phân tích và lý do lựa chọn: Thuật toán lọc Kalman mở rộng không mùi (Unscented Kalman Filter - UKF) được lựa chọn để giải bài toán ngược xác định tham số độ trơn kappa từ dữ liệu đo đạc thực nghiệm. Lý do lựa chọn UKF là khả năng xử lý phi tuyến tính cao bậc mà không cần tính ma trận đạo hàm Jacobian phức tạp, đảm bảo nghiệm hội tụ bền vững. Đối với phân tích số học, phần mềm LS-DYNA được sử dụng để xây dựng mô hình phần tử hữu hạn 3D (FEM) kết hợp phương pháp thủy động lực học hạt mịn (Smoothed-Particle Hydrodynamics - SPH). SPH là phương pháp không lưới (mesh-free) theo hệ quy chiếu Lagrange, được ưu tiên sử dụng cho vùng đất sát vách ống nhằm loại trừ triệt để hiện tượng méo lưới khi ống chịu chuyển vị lớn vượt quá 30% đường kính ống.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích thực nghiệm và kiểm chứng mô hình số đã mang lại bốn phát hiện khoa học mang tính đột phá:

• Thứ nhất, tham số độ trơn phi tuyến kappa phản ánh chính xác cơ chế phá hoại của đất nền và phụ thuộc chặt chẽ vào độ chặt tương đối cùng tỷ số chôn sâu H/D. Trong cát xốp, kappa tăng liên tục từ 1,1 đến 1,5 khi tỷ số H/D tăng từ 1,5 đến 8,0 rồi đạt trạng thái bão hòa, tương ứng với sự chuyển giao từ cơ chế nêm bị động lên mặt đất sang cơ chế cày xuyên cục bộ. Ngược lại, trong cát chặt, tham số kappa tiến sát về giá trị 0, biểu thị sự chuyển pha đàn dẻo đột ngột và hiện tượng phá hoại giòn sau khi ống chỉ mới dịch chuyển một đoạn nhỏ từ 3 đến 5 mm.
• Thứ hai, hiện tượng thắt hẹp pinching làm thay đổi căn bản quy luật truyền lực phương đứng. Chiều dài vùng thắt hẹp tăng tỷ lệ thuận với bình phương chuyển vị cực đại theo phương thẳng đứng với hệ số delta_s dao động từ 0,05 đến 0,30, khiến độ cứng tiếp tuyến giảm về mức triệt tiêu trước khi vách ống tái va chạm với nền đất đáy.
• Thứ ba, hiệu ứng ghép nối hai trục làm suy giảm nghiêm trọng sức kháng tải trọng giới hạn của đường ống khi chuyển vị theo phương xiên. Khi góc dịch chuyển xiên theta_u tăng từ 0 độ (thuần phương ngang) lên 60 độ, lực phản kháng cực hạn phương ngang giảm hơn 70%, từ mức 8,0 kN/m xuống chỉ còn khoảng 2,0 kN/m.
• Thứ tư, dưới tác động của tải trọng quá độ từ trận động đất Kobe (biên độ dịch chuyển ngang đạt đỉnh 88 mm và phương đứng đạt 23 mm), mô hình BMBW tái hiện chính xác các chu trình trễ tiêu tán năng lượng của mô phỏng SPH, trong khi tiêu chuẩn ASCE song tuyến tính đánh giá sai lệch đỉnh phản lực ngang với mức vượt tải lên đến hơn 40%.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự khác biệt giữa mô hình BMBW và các tiêu chuẩn hiện hành nằm ở việc xem xét trạng thái ứng suất địa tĩnh ban đầu. Trong môi trường tự nhiên, ứng suất hữu hiệu theo phương đứng và phương ngang không đồng nhất, tạo ra góc lệch chi khoảng 20 độ giữa hướng gia số chuyển vị và hướng gia số lực phản kháng. Các lò xo Winkler độc lập truyền thống không thể ghi nhận hiện tượng trượt lệch phương này, dẫn đến đánh giá sai lệch khả năng chịu uốn và cắt của tuyến ống.

Khi so sánh với các nghiên cứu thực nghiệm của Nyman (1984) và kết quả mô phỏng số của Guo (2005), đường bao phá hoại giới hạn hai trục dạng elip phi đối xứng của mô hình BMBW thể hiện sự tương thích vượt trội, loại bỏ hoàn toàn các điểm kỳ dị gián đoạn độ cứng vốn xuất hiện trong mô hình song tuyến ASCE.

Về phương diện trực quan hóa, toàn bộ dữ liệu phản lực - chuyển vị chu kỳ có thể được biểu diễn trực quan thông qua các đồ thị vòng trễ FDC hai trục (Fx-ux và Fy-uy) cho thấy rõ độ cong mượt tại các góc tải. Bảng thống kê đối chuẩn chỉ số độ cứng cát tuyến K70 và độ cứng tiếp tuyến BMBW theo các góc phương vị tải trọng khác nhau (0 độ, 30 độ, 45 độ, 60 độ, 90 độ) cung cấp cái nhìn định lượng chuẩn xác cho kỹ sư kết cấu trong việc đánh giá mức độ tiêu tán năng lượng kháng chấn.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu toàn diện, bốn giải pháp kỹ thuật cụ thể được đề xuất nhằm nâng cao năng lực thiết kế và vận hành an toàn hệ thống đường ống chôn ngầm:

• Tích hợp module phần tử vĩ mô BMBW vào các phần mềm mô phỏng kết cấu công trình ngầm chuyên dụng (như OpenSees hoặc CSI SAP2000) nhằm nâng cao độ chính xác dự báo phản lực động đất lên hơn 35%, hoàn thành trong khung thời gian 12 đến 18 tháng, chủ trì thực hiện bởi các viện nghiên cứu chuyên ngành công trình ngầm và địa kỹ thuật.
• Chuẩn hóa quy trình hiệu chuẩn thông số lò xo phi tuyến động lực học thông qua việc ứng dụng thuật toán lọc Kalman mở rộng không mùi UKF từ dữ liệu khảo sát hiện trường, áp dụng bắt buộc cho 100% các dự án đường ống đi qua khu vực có nguy cơ đứt gãy hoạt động trong giai đoạn thiết kế cơ sở từ năm 2026 đến năm 2028, thực hiện bởi các kỹ sư địa kỹ thuật công trình và đơn vị tư vấn thiết kế.
• Biên soạn và ban hành tài liệu hướng dẫn kỹ thuật mới thay thế tiêu chuẩn song tuyến tính ASCE 1984, bắt buộc đưa các hiệu ứng ghép nối hai trục và hiệu ứng thắt hẹp đáy ống vào tính toán đối với mọi đường ống có tỷ số chôn sâu H/D từ 1,5 đến 8,0, hướng đến mục tiêu giảm thiểu 50% nguy cơ sự cố nứt vỡ ống do biến dạng sóng Rayleigh trong thời hạn 2 năm, chịu trách nhiệm bởi các hội đồng khoa học và hiệp hội kỹ sư dầu khí.
• Áp dụng giải pháp hàm trở kháng động SIF cho nền phân lớp để tối ưu hóa chiều sâu chôn lấp tại các vị trí tiếp giáp địa tầng có tỷ số tương phản mô đun trượt mu1/mu2 bằng 0,25, tiến hành rà soát định kỳ 6 tháng một lần đối với các tuyến ống truyền dẫn năng lượng huyết mạch, thực hiện bởi các tổng công ty quản lý và vận hành hạ tầng kỹ thuật đô thị.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả của luận văn mang lại giá trị khoa học và ứng dụng thực tiễn cho bốn nhóm đối tượng chính:

• Kỹ sư thiết kế địa kỹ thuật và công trình ngầm: Tiếp cận phương pháp xác định chính xác sức kháng phi tuyến của đất, tối ưu hóa đường kính ống từ 100 mm đến 1000 mm và chiều sâu đặt ống nhằm giảm thiểu khối lượng vật liệu thép mà vẫn đảm bảo an toàn chịu lực.
• Nhà nghiên cứu và học viên sau đại học ngành Kỹ thuật Xây dựng - Động đất: Tham khảo khung toán học giải tích giải bài toán biên nửa không gian vô hạn bằng chuỗi Bessel-Fourier và kỹ thuật xử lý dữ liệu ngược UKF phục vụ phát triển các mô hình tương tác kết cấu - đất nền thế hệ mới.
• Đơn vị quản lý và vận hành mạng lưới đường ống năng lượng, cấp thoát nước: Ứng dụng mô hình bậc giảm để đánh giá nhanh mức độ tổn thương của hàng nghìn km đường ống khi xảy ra động đất có độ lớn từ 6,0 đến 8,0 độ Richter, từ đó xây dựng kế hoạch ứng phó sự cố khẩn cấp.
• Ban soạn thảo quy chuẩn xây dựng và cơ quan quản lý nhà nước: Sử dụng các bằng chứng khoa học thực nghiệm và số học trong luận văn làm căn cứ cập nhật, sửa đổi các quy phạm thiết kế kháng chấn công trình ngầm tại các vùng địa chấn nguy hiểm.

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

Mô hình phần tử vĩ mô BMBW có điểm gì ưu việt hơn so với mô hình lò xo ASCE truyền thống?
Mô hình BMBW khắc phục nhược điểm của tiêu chuẩn ASCE bằng việc mô tả đường cong lực - chuyển vị phi tuyến trơn thực tế, tái hiện vòng lặp trễ tiêu tán năng lượng và hiệu ứng ghép nối chuyển vị hai phương. Nhờ đó, mô hình giảm sai số dự báo phản lực từ hơn 40% xuống dưới 5% so với mô phỏng chi tiết.

Tại sao hiệu ứng thắt hẹp pinching chỉ xuất hiện khi đường ống dịch chuyển theo phương thẳng đứng hướng xuống?
Dưới tác dụng trọng trường, khi ống bị nhổ lên phía trên, lớp đất phủ nâng lên tạo thành khoảng trống khe hở cục bộ dưới đáy ống. Khi đổi chiều dịch chuyển xuống dưới, ống chuyển vị tự do trong khe hở với sức kháng gần bằng 0 trước khi chạm đáy đất, tạo nên vùng thắt hẹp đặc thù với hệ số delta_s từ 0,05 đến 0,30.

Làm thế nào để xác định tham số độ trơn kappa cho từng loại đất cát khác nhau?
Tham số kappa được nhận dạng tự động qua thuật toán lọc Kalman mở rộng không mùi UKF từ dữ liệu nén nhổ ống. Với cát chặt có độ chặt trên 65%, kappa xấp xỉ 0 phản ánh phá hoại giòn đạt đỉnh sau 3 đến 5 mm chuyển vị. Với cát xốp có độ chặt dưới 35%, kappa đạt từ 1,1 đến 1,5 khi tỷ số H/D vượt ngưỡng 6 đến 8.

Giải pháp hàm trở kháng động SIF được ứng dụng như thế nào trong phân tích sóng mặt Rayleigh?
Nghiên cứu thiết lập nghiệm giải tích giải bài toán biên nửa không gian đàn hồi bằng chuỗi trực giao. Các giá trị độ cứng và cản nhớt phụ thuộc tần số từ SIF được tích hợp trực tiếp vào hệ lò xo giảm chấn của mô hình bậc giảm, giúp mô phỏng chính xác đáp ứng động lực học tức thời với thời gian tính toán giảm hơn 90% so với phương pháp phần tử hữu hạn 2D.

Những giới hạn kỹ thuật nào cần lưu ý khi ứng dụng mô hình nghiên cứu vào thực tế?
Mô hình hiện tại được tối ưu hóa cho đường ống có độ cứng lớn đặt trong môi trường đất cát không dính. Mô hình chưa tích hợp đặc tính suy giảm độ bền sau đỉnh của cát rất chặt và chưa xét đến biến dạng méo mặt cắt ngang của ống thành mỏng dẻo, vốn đòi hỏi các nghiên cứu hiệu chỉnh bổ sung cho ống vật liệu mềm.

Kết luận

• Luận văn đã xây dựng thành công mô hình phần tử vĩ mô bậc giảm BMBW hai trục, mô tả toàn diện tính phi tuyến trơn, vòng trễ tiêu tán năng lượng và hiệu ứng thắt hẹp khe hở tương tác đất - ống.
• Ứng dụng thành công thuật toán lọc Kalman mở rộng không mùi UKF để lượng hóa chính xác tham số độ trơn kappa với sai số dưới 5% trên toàn bộ dải tỷ số chôn lấp H/D từ 1,5 đến 11,0.
• Phát triển nghiệm giải tích giải bài toán biên nửa không gian vô hạn, xác định hàm trở kháng động lực học SIF cho ống ngầm chịu sóng địa chấn điều hòa trong nền đồng nhất và phân lớp.
• Kiểm chứng độ tin cậy cơ học thông qua sự thống nhất chặt chẽ với phương pháp hạt mịn SPH, mô phỏng phần tử hữu hạn 3D LS-DYNA và sóng địa chấn Kobe thực tế với biên độ dịch chuyển tới 88 mm.
• Cắt giảm hơn 90% thời gian tính toán so với các mô hình phần tử hữu hạn liên tục 3D quy mô lớn mà vẫn duy trì độ chính xác kỹ thuật vượt trội.

Đóng góp lớn nhất của công trình là cung cấp một giải pháp tính toán bậc giảm hiệu quả cao, giải quyết triệt để bài toán tương tác đất - ống phi tuyến phức tạp trong phân tích kháng chấn mạng lưới hạ tầng ngầm diện rộng.

Kế hoạch phát triển tiếp theo tập trung vào việc mở rộng mô hình cho môi trường đất sét dính bão hòa nước và tích hợp quy luật suy giảm độ bền sau đỉnh trong giai đoạn 2026 - 2027.

Các kỹ sư kết cấu, chuyên gia địa kỹ thuật và nhà quản lý hạ tầng năng lượng được khuyến khích tham khảo, khai thác và tích hợp khung mô hình bậc giảm tiên tiến này vào các dự án quy hoạch, thiết kế và đánh giá an toàn công trình ngầm trọng điểm.