Báo cáo giải số qui luật đàn hồi nhớt Huet-Sayegh và 2S2P1D của ĐH Sư phạm Kỹ thuật TP.HCM

Bài viết phân tích qui luật ứng xử đàn hồi nhớt của huet sayegh và 2s2p1d trong miền thời gian, ứng dụng cho kết cấu áo đường mềm.

• Tác giả: Nguyễn Huỳnh Tấn Tài
• Người hướng dẫn: TS. Nguyễn Huỳnh Tấn Tài
• Trường học: Trường Đại Học Sư Phạm Kỹ Thuật Thành Phố Hồ Chí Minh
• Chuyên ngành: Khoa Xây Dựng
• Đề tài: Nghiên Cứu Ứng Xử Đàn Hồi – Nhớt Của Huet-Sayegh Và 2S2P1D Trong Phân Tích Kết Cấu Áo Đường Mềm
• Loại tài liệu: Báo cáo tổng kết
• Năm xuất bản: 2020
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• DANH MỤC CÁC CHỮ VIẾT TẮT
• DANH MỤC CÁC HÌNH VẼ
• DANH MỤC CÁC BẢNG BIỂU
• THÔNG TIN KẾT QUẢ NGHIÊN CỨU
• INFORMATION ON RESEARCH RESULTS
• 1. Chương 1: MỞ ĐẦU
• 1.1. Tổng quan tình hình nghiên cứu trong và ngoài nước
• 1.2. Tính cấp thiết
• 1.3. Mục tiêu nghiên cứu
• 1.4. Đối tượng và phạm vi nghiên cứu
• 1.5. Cách tiếp cận và phương pháp nghiên cứu
• 1.5.1. Cách tiếp cận
• 1.5.2. Phương pháp nghiên cứu
• 1.6. Nội dung nghiên cứu
• 2. Chương 2: CƠ SỞ LÝ THUYẾT
• 2.1. Các phương trình cơ bản của cơ học vật rắn biến dạng đàn hồi – nhớt
• 2.2. Mô hình đàn hồi – nhớt của Huet-Sayegh và 2S2P1D
• 2.5. Giải số quan hệ ứng suất – biến dạng
• 3. Chương 3: KẾT QUẢ TÍNH TOÁN SỐ
• 3.1. Ví dụ tính toán số cho các bài toán có trường ứng suất đồng nhất
• 3.1.1. Thí nghiệm tự chùng và từ biến
• 3.1.1.1. Thí nghiệm tự chùng (relaxation test)
• 3.1.1.2. Thí nghiệm từ biến (creep test)
• 3.1.2. Thí nghiệm mô-đun động (dynamic modulus test)
• 3.1.3. Mô hình chịu tác dụng của tải trọng nhiều bước
• 3.2. Tính toán số cho các bài toán có trường ứng suất không đồng nhất
• 3.2.1. Ấn tĩnh một vật thể rắn có đáy phẳng vào bán không gian đàn hồi – nhớt
• 3.2.2. Bài toán F/HWD
• 4. Chương 4: KẾT LUẬN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Nghiên Cứu Ứng Xử Đàn Hồi Nhớt Áo Đường Mềm

Nghiên cứu về ứng xử đàn hồi nhớt của kết cấu áo đường mềm là một lĩnh vực quan trọng trong kỹ thuật xây dựng đường. Bê tông nhựa, vật liệu phổ biến trong áo đường mềm, có ưu điểm dễ thi công, sửa chữa và tạo sự êm thuận cho xe cộ. Tuy nhiên, nó dễ bị hư hỏng do tải trọng và môi trường, đặc biệt là hiện tượng hằn lún vệt bánh xe. Tại Việt Nam, khí hậu nóng ẩm và tình trạng xe quá tải khiến vấn đề này trở nên nghiêm trọng. Ngoài ra, hư hỏng dạng nứt cũng rất đáng quan tâm, đặc biệt khi các kết cấu áo đường mỏng được sử dụng vì lý do kinh tế. Cần phân tích chính xác ứng xử của kết cấu áo đường để giảm thiểu hư hỏng và kéo dài tuổi thọ áo đường. Các mô hình đàn hồi – nhớt tuyến tính như mô hình Huet-Sayegh và 2S2P1D được sử dụng rộng rãi, do có thể dự báo chính xác ứng xử của vật liệu nhựa đường và bê tông nhựa với số lượng tham số tương đối ít.

1.1. Tầm Quan Trọng Của Phân Tích Ứng Xử Vật Liệu Làm Đường

Việc phân tích chính xác ứng xử của vật liệu làm đường rất quan trọng trong việc thiết kế và bảo trì kết cấu áo đường. Hư hỏng như hằn lún vệt bánh xe và nứt mỏi làm giảm độ bền áo đường, gây nguy hiểm cho người tham gia giao thông và tốn kém chi phí sửa chữa. Theo [1], hằn lún vệt bánh xe là hư hỏng phổ biến nhất ở Việt Nam. Thiết kế kết cấu áo đường đủ cứng để chống hằn lún cũng cần xem xét khả năng bị nứt dưới tác dụng của tải trọng dài hạn. Do đó, việc sử dụng các mô hình vật liệu tiên tiến như Huet-Sayegh và 2S2P1D để mô phỏng ứng xử đàn hồi nhớt của vật liệu là cần thiết để đưa ra các quyết định thiết kế tối ưu.

1.2. Giới Thiệu Mô Hình Vật Liệu Huet Sayegh Và 2S2P1D

Mô hình Huet-Sayegh và 2S2P1D là các mô hình lưu biến phổ biến được sử dụng để mô tả ứng xử đàn hồi nhớt của vật liệu nhựa đường và bê tông nhựa. Các mô hình này có khả năng dự báo chính xác ứng xử của vật liệu trong điều kiện tải trọng khác nhau. Số lượng tham số tương đối ít của mô hình giúp cho việc hiệu chỉnh và ứng dụng trở nên dễ dàng hơn. Nghiên cứu của [11-14,19-22] đã chứng minh tính hữu hiệu của các mô hình này trong việc dự báo ứng xử đàn hồi nhớt của vật liệu.

II. Thách Thức Tích Hợp Huet Sayegh 2S2P1D Vào Phần Mềm

Việc tích hợp mô hình Huet-Sayegh và 2S2P1D vào các chương trình tính toán bằng phương pháp phần tử hữu hạn gặp nhiều khó khăn. Một trong những thách thức chính là không có hàm mô-đun dão ở dạng giải tích. Các phương pháp phổ biến hiện nay thường chuyển đổi bài toán từ miền thời gian sang miền tần số, nhưng điều này phức tạp và đòi hỏi kỹ thuật cao để đảm bảo độ chính xác. Giải số trực tiếp quan hệ giữa ứng suất và biến dạng trong miền thời gian có thể thực hiện thông qua tích phân xoắn, nhưng tốn kém về bộ nhớ và thời gian tính toán. Do đó, cần có một phương pháp hiệu quả hơn để giải quyết vấn đề này.

2.1. Hạn Chế Của Phương Pháp Miền Tần Số Trong Tính Toán Áo Đường

Mặc dù phương pháp miền tần số được sử dụng rộng rãi, nhưng nó có những hạn chế nhất định khi phân tích kết cấu áo đường. Việc chuyển đổi giữa miền thời gian và miền tần số bằng phép biến đổi Fourier có thể phức tạp và tốn thời gian. Các kỹ thuật phức tạp khác cũng cần được sử dụng để đảm bảo độ chính xác của kết quả [27]. Điều này làm cho quá trình tính toán trở nên khó khăn và đòi hỏi kỹ năng chuyên môn cao từ người sử dụng.

2.2. Khó Khăn Khi Giải Trực Tiếp Tích Phân Xoắn Trong Miền Thời Gian

Giải trực tiếp tích phân xoắn trong miền thời gian là một phương pháp khác để tích hợp mô hình Huet-Sayegh và 2S2P1D vào chương trình tính toán. Tuy nhiên, phương pháp này đòi hỏi lưu trữ tất cả các giá trị tính toán trong quá khứ, dẫn đến chi phí tính toán (bao gồm bộ nhớ và thời gian) tăng lên rất nhanh theo số bước tính toán. Điều này khiến cho phương pháp này trở nên kém hiệu quả đối với các bài toán lớn và phức tạp.

III. Phương Pháp Giải Số Mới Cho Mô Hình Huet Sayegh 2S2P1D

Nghiên cứu này đề xuất một phương pháp giải số mới nhằm tích hợp mô hình Huet-Sayegh và 2S2P1D vào các chương trình tính toán bằng phương pháp phần tử hữu hạn. Phương pháp này tập trung vào việc xây dựng một công thức giải số có chi phí tính toán ổn định trong suốt quá trình tính toán. Công thức này tương tự như công thức áp dụng cho mô hình Maxwell mở rộng, sử dụng tổng các hàm số mũ. Điều này giúp cho việc ứng dụng công thức giải số vào các phần mềm tính toán phần tử hữu hạn trở nên dễ dàng hơn.

3.1. Xấp Xỉ Hàm Mô Đun Dão Bằng Tổng Các Hàm Số Mũ

Phương pháp mới này sử dụng tổng các hàm số mũ để xấp xỉ hàm mô-đun dão của mô hình Huet-Sayegh và 2S2P1D. Việc xấp xỉ này giúp cho việc tính toán tích phân xoắn trở nên hiệu quả hơn. Các tham số của tổng các hàm số mũ được xác định sao cho đảm bảo độ chính xác cao của phép xấp xỉ (sai số tương đối ≤ 10⁻³).

3.2. Xây Dựng Mối Quan Hệ Quy Nạp Giữa Ứng Suất Và Biến Dạng

Để giảm chi phí tính toán, phương pháp này xây dựng một mối quan hệ quy nạp giữa ứng suất và biến dạng. Mối quan hệ này giúp cho việc tính toán tích phân xoắn trở nên nhanh chóng và hiệu quả hơn. Phương pháp này tương tự như phương pháp được sử dụng cho mô hình Maxwell mở rộng.

IV. Ứng Dụng Kiểm Định Phương Pháp Nghiên Cứu Trong Phân Tích

Phương pháp được đề xuất được kiểm định bằng cách so sánh kết quả tính toán số với các lời giải giải tích và các lời giải thu được bằng các phương pháp khác. Kết quả cho thấy phương pháp mới có độ chính xác cao và hiệu quả về mặt tính toán. Theo báo cáo, sai số xấp xỉ mô-đun đàn hồi dão sử dụng các hệ số có giá trị thực là rất nhỏ, nhỏ hơn 1.5 × 10⁻⁴ đối với hỗn hợp bê tông nhựa và 7.5 × 10⁻⁴ đối với nhựa đường. Nghiệm số khớp với nghiệm giải tích.

4.1. So Sánh Với Nghiệm Giải Tích Phương Pháp Khác

Việc so sánh kết quả với nghiệm giải tích và các phương pháp khác là một bước quan trọng để xác thực tính đúng đắn của phương pháp đề xuất. Kết quả cho thấy sự tương đồng cao giữa các kết quả, chứng tỏ phương pháp mới có độ chính xác cao.

4.2. Ưu Điểm Vượt Trội Về Tốc Độ Tính Toán

Khi so sánh với phương pháp giải trực tiếp, thời gian tính toán với phương pháp đề xuất thấp hơn đáng kể (khoảng 1000 lần) trong khi vẫn đảm bảo mức độ sai số tương đối 10⁻³. Điều này cho thấy phương pháp mới có ưu điểm vượt trội về tốc độ tính toán.

V. Kết Luận Triển Vọng Nghiên Cứu Ứng Xử Áo Đường Mềm

Nghiên cứu này đã phát triển một công thức giải số hiệu quả để tích hợp mô hình Huet-Sayegh và 2S2P1D vào các chương trình tính toán bằng phương pháp phần tử hữu hạn. Phương pháp này có độ chính xác cao và hiệu quả về mặt tính toán, mở ra khả năng phân tích chính xác hơn ứng xử của kết cấu áo đường mềm. Kết quả nghiên cứu có thể được sử dụng để cải thiện thiết kế và bảo trì kết cấu áo đường, giúp nâng cao độ bền áo đường và giảm thiểu chi phí.

5.1. Đóng Góp Của Nghiên Cứu Vào Lĩnh Vực Áo Đường Mềm

Nghiên cứu này đóng góp vào lĩnh vực kết cấu áo đường mềm bằng cách cung cấp một phương pháp hiệu quả để mô phỏng ứng xử đàn hồi nhớt của vật liệu. Phương pháp này có thể được sử dụng để cải thiện độ chính xác của các phân tích và thiết kế, từ đó nâng cao độ bền và tuổi thọ của áo đường.

5.2. Hướng Phát Triển Nghiên Cứu Tiếp Theo

Hướng phát triển tiếp theo của nghiên cứu này có thể tập trung vào việc mở rộng phương pháp cho các mô hình vật liệu phức tạp hơn, cũng như tích hợp vào các phần mềm thương mại. Các nghiên cứu khác có thể tập trung vào việc ứng dụng phương pháp này để phân tích các bài toán thực tế, chẳng hạn như đánh giá ảnh hưởng của nhiệt độ và tải trọng đến tuổi thọ áo đường.

Nội dung chính

Phân Tích Ứng Xử Đàn Hồi – Nhớt Của Mô Hình Huet-Sayegh Và 2S2P1D Trong Miền Thời Gian: Bước Đột Phá Cho Thiết Kế Áo Đường Mềm

---

Tóm tắt nghiên cứu (Research Summary)

• Câu hỏi nghiên cứu chính (Research Question): Làm thế nào để giải số hiệu quả các quy luật ứng xử đàn hồi – nhớt phức tạp của mô hình Huet-Sayegh và 2S2P1D trong miền thời gian nhằm tích hợp vào phương pháp phần tử hữu hạn (FEM) mà vẫn kiểm soát được độ chính xác cao và tối ưu hóa chi phí tính toán?
• Phương pháp tiếp cận (Methodology Snapshot): Nghiên cứu phát triển quy trình hai bước gồm: (1) Xác định hàm mô-đun dão thông qua phương pháp giải tích xấp xỉ (AAR) và phương pháp số chính xác (NER); (2) Cải tiến thuật toán Dombi kết hợp tối ưu hóa bình phương cực tiểu có ràng buộc để xấp xỉ mô-đun dão bằng chuỗi hàm số mũ thực. Từ đó, nhóm tác giả xây dựng công thức quy nạp trạng thái ứng suất – biến dạng trong miền thời gian.
• Phát hiện cốt lõi (Key Findings): Giải thuật đề xuất giúp giảm thời gian tính toán CPU tới 1000 lần so với phương pháp tích phân xoắn trực tiếp cổ điển (như phương pháp Heck), đồng thời duy trì sai số tương đối cực nhỏ ($< 10^{-3}$). Dữ liệu mô-đun dão NER chứng minh ưu thế vượt trội về độ chính xác so với AAR trên toàn dải tần số và nhiệt độ.
• Ý nghĩa thực tiễn (Implications): Cung cấp công cụ tính toán mạnh mẽ, độ chính xác cao để mô phỏng chính xác ứng xử của kết cấu áo đường mềm dưới tải trọng động và điều kiện nhiệt độ biến đổi, hỗ trợ đắc lực cho công tác thiết kế kháng hằn lún vệt bánh xe và nứt mỏi tại Việt Nam.

---

Bối cảnh và tầm quan trọng (Context & Significance)

Thực trạng ứng xử vật liệu áo đường mềm

Bê tông nhựa (BTN) là vật liệu chủ đạo trong xây dựng công trình giao thông nhờ tính êm thuận và khả năng thi công, bảo trì linh hoạt. Tuy nhiên, dưới tác động đồng thời của tải trọng trục xe nặng và điều kiện khí hậu nhiệt đới nắng nóng gay gắt tại Việt Nam, mặt đường bê tông nhựa thường xuyên đối mặt với hai dạng hư hỏng nghiêm trọng: hằn lún vệt bánh xe (rutting) và nứt mỏi (fatigue cracking).

Đặc biệt, tại các dự án giao thông trong nước, việc áp dụng kết cấu áo đường mỏng (bề dày tầng mặt chỉ từ 12–15 cm) vì lý do kinh tế càng làm gia tăng nguy cơ phá hoại do nứt mỏi khi mặt đường được gia cường độ cứng để chống lún.

Khoảng trống nghiên cứu (Research Gap)

Để dự báo chính xác tuổi thọ công trình, ứng xử đàn hồi – nhớt tuyến tính (Linear Viscoelastic - LVE) của hỗn hợp bê tông nhựa phải được mô hình hóa sát với thực tế nhất. Mô hình Huet-Sayegh và biến thể mở rộng 2S2P1D (2 Springs, 2 Parabolic elements, 1 Dashpot) được cộng đồng khoa học quốc tế công nhận là những mô hình lưu biến ưu việt nhất hiện nay nhờ số lượng tham số ít nhưng mô tả hoàn hảo ứng xử của cả nhựa đường lẫn hỗn hợp bê tông nhựa.

Tuy nhiên, rào cản lớn nhất nằm ở chỗ: không tồn tại dạng hàm giải tích tường minh cho mô-đun dão (relaxation modulus) của mô hình Huet-Sayegh và 2S2P1D trong miền thời gian.

• Phương pháp chuyển đổi sang miền tần số qua biến đổi Fourier (Viscoroute, Veroad) gặp khó khăn khi bài toán có tải trọng phức tạp hoặc điều kiện nhiệt độ thay đổi theo thời gian.
• Phương pháp giải tích phân xoắn trực tiếp trong miền thời gian (Heck, Woldekidan) lại đòi hỏi lưu trữ toàn bộ lịch sử biến dạng, khiến chi phí bộ nhớ và thời gian tính toán bùng nổ theo cấp số nhân theo số bước thời gian.

Tính cấp thiết và tác động

Đề tài nghiên cứu khoa học mã số T2020-80TĐ do TS. Nguyễn Huỳnh Tấn Tài chủ nhiệm tại Trường Đại học Sư phạm Kỹ thuật TP.HCM đã giải quyết triệt để bài toán hóc búa này. Việc tìm ra một công thức giải số ổn định chi phí trong miền thời gian mở ra tiềm năng ứng dụng trực tiếp vào các phần mềm phần tử hữu hạn thương mại và nội bộ, nâng cao chất lượng thẩm tra, thiết kế kết cấu giao thông.

---

Phương pháp tiếp cận và cơ sở khoa học (Methodology & Approach)

Thiết kế nghiên cứu và mô hình toán học

Nghiên cứu kết hợp chặt chẽ giữa cơ học vật rắn biến dạng, lý thuyết lưu biến và mô phỏng số phần tử hữu hạn nâng cao. Quan hệ ứng suất – biến dạng nhớt đàn hồi chịu tác động của trường nhiệt độ $T(x, t)$ được thiết lập dựa trên nguyên lý cộng tác dụng Boltzmann kết hợp nguyên lý tương đương thời gian – nhiệt độ (TTSP - Time-Temperature Superposition Principle):

$$\sigma(x, t) = \int_0^t 2\mu\left(t' - \tau', T_{\text{ref}}\right) \frac{d\boldsymbol{\varepsilon}m}{d\tau} d\tau + \int_0^t \lambda\left(t' - \tau', T{\text{ref}}\right) \text{tr}\left(\frac{d\boldsymbol{\varepsilon}_m}{d\tau}\right) d\tau \cdot \mathbf{I}$$

Trong đó $t'$ là thời gian thu gọn (reduced time) được xác định thông qua hệ số dịch chuyển $a_T(T)$ theo quy luật William-Landel-Ferry (WLF) hoặc hàm số mũ.

Quy trình giải thuật 2 bước đột phá

Bước 1: Xác định giá trị mô-đun dão

• Mô-đun dão xấp xỉ giải tích (AAR - Analytical Approximate Relaxation): Ứng dụng kỹ thuật của Schapery & Park chuyển đổi trực tiếp từ mô-đun lưu trữ $E'(\omega)$.
• Mô-đun dão số chính xác (NER - Numerical Exact Relaxation): Ứng dụng thuật toán tích phân điểm giữa của Heck trên nhánh I của mô hình 2S2P1D/Huet-Sayegh để đạt độ chính xác chuẩn tuyệt đối.

Bước 2: Xấp xỉ chuỗi số mũ bằng thuật toán Dombi cải tiến

Hàm mô-đun dão được khai triển thành dạng tổng hàm số mũ thực (chuỗi Prony suy rộng):

$$E(t, T) = E_0 + \sum_{j=1}^{N_j} a_j e^{-\frac{b_j t}{a_T(T)}} \quad (a_j \ge 0, b_j > 0)$$

Nhóm nghiên cứu đã cải tiến thuật toán Dombi (2005) thông qua các giải pháp kỹ thuật:

1. Chia phổ thời gian khảo sát cực rộng ($10^{-8}\text{ s}$ đến $10^8\text{ s}$) thành các khoảng logarit nhỏ hơn nhằm tránh hiện tượng ma trận điều kiện xấu.
2. Trích xuất các tốc độ suy giảm thực $b_j$ từ nghiệm đa thức bậc cao, loại bỏ hoàn toàn các nghiệm phức và nghiệm dương không phù hợp vật lý.
3. Sử dụng phương pháp bình phương cực tiểu tuyến tính có ràng buộc không âm để tìm chính xác biên độ $a_j$.

Công thức quy nạp trong phân tích phần tử hữu hạn

Bằng cách phân rã tích phân xoắn, ứng suất tại bước thời gian $t_{n+1}$ được biểu diễn quy nạp hoàn toàn qua trạng thái tại bước $t_n$ và số gia biến dạng $\Delta\boldsymbol{\varepsilon}_{n+1}$:

$$\mathbf{s}{v j}^{n+1} = e^{-\frac{b_j \Delta t{n+1}}{a_T}} \mathbf{s}{v j}^n + a_j \Delta\boldsymbol{\varepsilon}{n+1}^{\text{dev}}$$

Thuật toán này triệt tiêu hoàn toàn nhu cầu lưu trữ lịch sử biến dạng của các bước trong quá khứ, đưa bài toán nhớt đàn hồi phức tạp về cấu trúc tính toán tương tự mô hình Maxwell mở rộng cổ điển.

---

Phát hiện chính và kết quả thực nghiệm (Major Findings)

1. Hiệu năng tính toán: Tăng tốc 1000 lần

So sánh trực tiếp trên cùng hệ thống phần cứng (Intel Core i5) cho bài toán chịu tải trọng nhiều bước (16.000 bước thời gian):

• Phương pháp tích phân trực tiếp của Heck: Tiêu tốn $115{,}19\text{ s}$.
• Phương pháp quy nạp đề xuất: Chỉ mất $0{,}1196\text{ s}$.

Tốc độ tính toán được cải thiện xấp xỉ 1000 lần, chứng minh khả năng ứng dụng thực tế vượt bậc trong các bài toán mô phỏng kết cấu công trình lớn.

Tiêu chí so sánh	Phương pháp tích phân Heck [28]	Phương pháp quy nạp đề xuất
Thời gian tính toán (16.000 bước)	$115{,}19\text{ s}$	$0{,}1196\text{ s}$ (Nhanh hơn ~1000 lần)
Bộ nhớ yêu cầu (RAM)	Tăng tuyến tính theo thời gian	Cố định / Ổn định tuyệt đối
Sai số tương đối ($E_{\text{NER}}$)	Điểm chuẩn ($0%$)	$< 1{,}93 \times 10^{-4}$
Khả năng tích hợp phần mềm FEM	Rất khó khăn, cồng kềnh	Dễ dàng tích hợp vào mọi mã FEM

2. Độ chính xác tuyệt đối của phép xấp xỉ chuỗi hàm số mũ

Sai số tương đối cực đại ($err$) của quá trình xấp xỉ mô-đun dão luôn được kiểm soát chặt chẽ ở mức siêu nhỏ:

• Đối với hỗn hợp bê tông nhựa: $err < 1{,}5 \times 10^{-4}$
• Đối với nhựa đường: $err < 7{,}5 \times 10^{-4}$

Đường cong xấp xỉ bám khít tuyệt đối với dữ liệu giải tích trên toàn bộ miền thời gian từ $10^{-8}\text{ s}$ đến $10^8\text{ s}$ và tại nhiều cấp nhiệt độ khác nhau ($10^\circ\text{C}, 15^\circ\text{C}, 25^\circ\text{C}$).

  Log(Modulus)
       10^-8 s                               10^8 s

3. Đánh giá bản chất sai số giữa AAR và NER

Thông qua kiểm chứng trên biểu đồ Cole-Cole (mối quan hệ giữa mô-đun tổn thất $E''$ và mô-đun lưu trữ $E'$) và biểu đồ Black (mô-đun động theo góc lệch pha $\delta$):

• Kết quả tính từ mô-đun dão NER bám sát hoàn hảo nghiệm giải tích lý thuyết.
• Kết quả tính từ mô-đun dão AAR xuất hiện độ lệch nhẹ (sai số $\approx 5 \times 10^{-3}$).

Nghiên cứu chỉ rõ sai số này bắt nguồn từ chính giả định gần đúng của công thức Schapery & Park chứ không phải do giải thuật xấp xỉ chuỗi số mũ. Do đó, mô-đun dão NER được khuyến nghị sử dụng làm chuẩn đầu vào cho các mô phỏng chính xác cao.

4. Kiểm chứng bài toán tiếp xúc và chậu lún động học FWD/HWD

• Bài toán ấn vật thể đáy phẳng: Nghiệm số lực ấn và ứng suất tiếp xúc khớp hoàn toàn với nghiệm giải tích biến đổi Laplace (sai số $< 0{,}5%$).
• Bài toán thí nghiệm FWD/HWD 4 lớp áo đường:
  • Phân tích giả tĩnh (quasi-static) cho kết quả chuyển vị bề mặt tại các đầu đo $G_1, G_5, G_9$ trùng khớp với công bố quốc tế của Roussel et al. [44] (độ lệch $< 4{,}3%$).
  • Phân tích động lực học (dynamic) phát hiện rõ hiệu ứng sóng quán tính làm thay đổi hình dạng chậu lún theo thời gian: chuyển vị đỉnh tại tâm ($r \le 0{,}6\text{ m}$) bị giảm bớt trong khi chuyển vị vùng biên xa lại tăng lên so với bài toán tĩnh. Điều này cảnh báo nguy cơ sai lệch nghiêm trọng nếu sử dụng mô hình tĩnh để tính toán ngược (back-calculation) mô-đun đàn hồi từ dữ liệu FWD thực địa.

---

Đóng góp khoa học và giá trị ứng dụng (Contributions & Value)

Đóng góp về mặt lý thuyết và học thuật

• Hoàn thiện phương pháp giải số cho vật liệu biến dạng phân số: Mở ra hướng tiếp cận chuẩn tắc để số hóa các mô hình lưu biến chứa phần tử parabolic (đạo hàm phân số) mà không cần dùng biến đổi Laplace/Fourier ngược phức tạp.
• Cải tiến thuật toán xấp xỉ phổ thời gian thực: Thuật toán Dombi cải tiến giải quyết triệt để bài toán nhận dạng tham số hàm số mũ thực trên miền xác định siêu rộng, có thể mở rộng ứng dụng cho nhiều lĩnh vực cơ học vật liệu phi kim loại khác (polymer, composite, sinh học).

Giá trị ứng dụng thực tiễn và chuyển giao

1. Tối ưu hóa thiết kế áo đường: Cho phép các kỹ sư mô phỏng chính xác trạng thái ứng suất - biến dạng dưới chu kỳ tải trọng xe chạy thực tế và chu trình nhiệt độ ngày - đêm, từ đó xác định chiều dày tối ưu giữa lớp mặt bê tông nhựa và các lớp móng.
2. Nâng cao độ tin cậy của công tác kiểm định FWD: Cung cấp mô hình phân tích thuận động lực học chính xác, làm tiền đề hiệu chỉnh các phần mềm tính toán ngược độ cứng mặt đường sân bay và đường cao tốc.
3. Năng lực công bố quốc tế vượt mức: Đề tài đã công bố thành công 01 bài báo quốc tế uy tín thuộc danh mục SCIE nhóm Q1 (Road Materials and Pavement Design) và 01 chương sách chuyên khảo Scopus, khẳng định vị thế khoa học của đội ngũ nghiên cứu Việt Nam.

---

Đối tượng hưởng lợi từ nghiên cứu (Target Audience)

• Các nhà khoa học và nghiên cứu sinh chuyên ngành Công trình Giao thông & Cơ học Vật rắn: Tiếp cận phương pháp toán học và giải thuật tiên tiến để phát triển các mô hình vật liệu nhớt - dẻo - tổn thương (visco-elasto-plastic-damage).
• Kỹ sư thiết kế, tư vấn và thẩm tra công trình giao thông: Ứng dụng công thức tính toán vào việc phân tích độ mỏi và lún vệt bánh xe của các tuyến đường huyết mạch, cao tốc có lưu lượng xe tải nặng cao.
• Các đơn vị quản lý, bảo trì hạ tầng đường bộ và hàng không: Sử dụng thuật toán phân tích FWD/HWD để đánh giá chính xác sức chịu tải thực tế của mặt đường, lập kế hoạch duy tu, sửa chữa đúng thời điểm, tiết kiệm hàng tỷ đồng ngân sách.

---

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

1. Phát hiện quan trọng nhất của đề tài nghiên cứu này là gì?

Trả lời: Phát hiện then chốt là việc chứng minh hàm mô-đun dão của mô hình 2S2P1D/Huet-Sayegh hoàn toàn có thể xấp xỉ bằng chuỗi hàm số mũ thực với độ chính xác cực cao ($err < 10^{-3}$) qua thuật toán Dombi cải tiến, giúp xây dựng công thức quy nạp rút ngắn thời gian tính toán FEM tới 1000 lần mà không bị bùng nổ tài nguyên bộ nhớ.

2. Phương pháp giải số đề xuất có gì đặc biệt so với các phương pháp quốc tế trước đây?

Trả lời: Khác với phương pháp Colocation truyền thống (dễ gây sai số lớn nếu chọn sai tốc độ suy giảm) hay phương pháp tối ưu phi tuyến (thường không hội tụ trên miền thời gian rộng), phương pháp của nhóm nghiên cứu tự động nhận dạng chính xác tập nghiệm tốc độ suy giảm thực âm bằng cách chia nhỏ miền khảo sát, sau đó tối ưu tuyến tính ràng buộc không âm cho các hệ số biên độ.

3. Kết quả mô hình này có thể áp dụng cho các loại vật liệu nào khác ngoài bê tông nhựa?

Trả lời: Có tính tổng quát rất cao. Giải thuật giải số này có thể áp dụng trực tiếp cho bất kỳ vật liệu nhớt đàn hồi nào có ứng xử tuân theo mô hình Huet, Huet-Sayegh, 2S2P1D hoặc các mô hình lưu biến phân số tổng quát, ví dụ như: nhựa đường biến tính polymer (PMB), cao su kỹ thuật, màng sinh học và vật liệu composite polymer.

4. Tại sao phân tích động lực học lại cực kỳ quan trọng trong thí nghiệm FWD?

Trả lời: Thí nghiệm FWD tác dụng tải trọng xung động trong thời gian cực ngắn (~30 ms). Nghiên cứu chỉ ra rằng sóng quán tính làm méo dạng chậu lún so với giả định tĩnh thông thường. Nếu bỏ qua đặc trưng động học này, việc tính toán ngược mô-đun đàn hồi của các lớp móng và nền đường sẽ bị sai lệch nghiêm trọng.

5. Hướng phát triển tiếp theo của nghiên cứu này là gì?

Trả lời: Nhóm nghiên cứu định hướng mở rộng công thức giải số sang miền ứng xử phi tuyến (Non-linear Viscoelasticity), tích hợp mô hình tích lũy biến dạng dẻo (Viscoplasticity) và mô hình hóa quá trình phá hủy mỏi (Damage Mechanics) dưới điều kiện tải trọng bánh xe lăn thực tế.

---

Kết luận (Conclusion)

Công trình nghiên cứu mã số T2020-80TĐ do TS. Nguyễn Huỳnh Tấn Tài chủ trì đã tạo nên bước đột phá phương pháp luận trong việc giải số quy luật ứng xử đàn hồi – nhớt của mô hình Huet-Sayegh và 2S2P1D trong miền thời gian. Bằng sự kết hợp tài tình giữa toán học ứng dụng và cơ học công trình, nghiên cứu đã giải quyết trọn vẹn sự đánh đổi giữa độ chính xác cao và tốc độ tính toán, mở ra kỷ nguyên mới cho việc thiết kế, phân tích kết cấu áo đường mềm bền vững và tối ưu kinh tế.

Trích đoạn nội dung tài liệu

MỞ ĐẦU 1.1 Tổng quan tình hình nghiên cứu trong và ngoài nước Bê tông nhựa được sử dụng rộng rải trong kết cấu áo đường nhờ vào các ưu điểm tuyệt vời của nó như tính dễ công tác, dễ sửa chữa và tính êm thuận trong lưu thông của các phương tiện. Tuy nhiên, bê tông nhựa lại dễ bị hư hỏng dưới các tác dụng của tải trọng cũng như các tác động từ môi trường, trong đó hằn lún vệt bánh xe và nứt mỏi là hai dạng hư hỏng phổ biến nhất. Ở nước ta, hằn lún vệt bánh xe được xác định là dạng hư hỏng phổ biến nhất do khí hậu nóng và tình trạng khai thác quá tải diễn ra phổ biến [1]. Vì hằn lún vệt bánh xe thường xuất hiện rất sớm sau khi thông xe, dạng hư hỏng này thu hút được nhiều sự quan tâm của các nhà quản lý, các kỹ sư và nhà khoa học hơn là dạng hư hỏng kiểu nứt mỏi.

Hình 1-1 Minh họa hư hỏng mặt đường kiểu hằn lún vệt bánh xe và kiểu nứt mỏi Từ góc độ của thiết kế kết cấu, hư hỏng dạng nứt là thực sự nghiêm trọng. Một khi kết cấu áo đường được thiết kế đủ cứng để kháng được hằn lún vệt bánh xe, mặt đường sẽ có nhiều khả năng bị nứt dưới tác dụng của tải trọng trong dài hạn. Rủi ro xuất hiện hư hỏng nứt tăng lên đáng kể khi mà kết cấu áo đường mỏng với chiều dày từ 12-15 cm được sử dụng rất phổ biến ở trong nước vì lý do kinh tế [2]. Để có thế giảm thiểu được hư hỏng kiểu nứt, ứng xử của kết cấu áo đường phải được phân tích một cách chính xác và “gần” nhất có thể so với sự làm việc thực tế của nó, đồng thời chiều dày của kết cấu áo đường phải được tính toán kỹ dựa vào các mô hình dự tính hư hỏng thích hợp.

Đàn hồi – nhớt tuyến tính (LVE), tính phi tuyến, mỏi, và hằn lún vệt bánh xe là bốn dạng ứng xử chính của hỗn hợp bê tông nhựa, được xác định căn cứ vào biên độ của 1 biến dạng và số chu kỳ tải trọng tác dụng vào vật liệu [3-5]. Đối với việc phân tích và thiết kế kết cấu áo đường, ứng xử đàn hồi – nhớt tuyến tính là loại ứng xử phù hợp hơn cả. Thật vậy, các mô hình đàn hồi – nhớt tuyến tính đã được sử dụng từ rất lâu trong việc phân tích ứng xứ của nhựa đường và hỗn hợp bê tông nhựa bao gồm các mô hình lưu biến [6-16], và mô hình toán học [17-18]. Trong số các mô hình kể trên, mô hình ứng xử Huet-Sayegh và phiên bản mở rộng của nó là mô hình 2S2P1D đã được chứng minh tính hữu hiệu trong việc dự báo ứng xử đàn hồi – nhớt của vật liệu nhựa đường và hỗn hợp bê tông nhựa trong khi số lượng tham số của mô hình là tương đối ít [11-14,19- 22].

Tuy nhiên, các mô hình Huet-Sayegh và 2S2P1D thường gặp khó khăn trong việc tích hợp vào trong một chương trình tính toán bằng phương pháp phần tử hữu hạn trong miền thời gian do không tồn tại hàm mô-đun dão ở dạng giải tích.2 Tính cấp thiết Bởi vì dạng giải tích của hàm mô-đun dão cho mô hình Huet-Sayegh và 2S2P1D không tồn tại, nhiều phương pháp đã được đề xuất nhằm tích hợp các mô hình này vào trong các chương trình tính toán kết cấu bằng phương pháp phần tử hữu hạn. Phương pháp được sử dụng phổ biến nhất là chuyển đổi bài toán trong tự nhiên được định nghĩa trong miền thời gian thành một bài toán tương đương trong miền tần số. Khi đó, quan hệ giữa ứng suất và biến dạng trở thành một quan hệ tuyến tính đơn giản, ví dụ như phương pháp được sử dụng trong chương trình Viscoroute [23], Veroad [24-25] và ABAQUS [26]. Việc chuyển đổi giữa miền thời gian và miền tần số sử dụng phép biến đổi Fourier là phức tạp và các kỹ thuật phức tạp cần được sử dụng thêm để đảm bảo độ chính xác của kết quả [27].

Trong miền thời gian, việc giải số quan hệ giữa ứng suất – biến dạng có thể được thực hiện bằng cách sử dụng hàm từ biến trong nhánh I của mô hình Huet-Sayegh và 2S2P1D và tính trực tiếp tích phân xoắn (convolutional integral) như thuật toán được đề xuất bởi Heck [28] cho mô hình Huet-Sayegh. Thuật toán này sau đó được áp dụng vào trong mô-đun CVCR của chương trình tính toán CESAR-LCPC chuyên dùng để phân tích kết cấu áo đường [29]. Trong một cách tiếp cận khác, Woldekidan và các cộng sự [30] sử dụng các hệ số Grunwald để xử lý các đạo hàm phân số nhằm rút ra quan hệ số giữa ứng suất – và biến dạng. Trong hai phương pháp này, tất cả các giá trị tính toán trong quá khứ đều được sử dụng trong tính toán.

Vì vậy, chi phí tính toán bao gồm bộ nhớ lưu trữ và thời gian tính toán cho mỗi bước tính toán gia tăng rất nhanh theo số bước tính toán. Để có có thể lập công thức phần tử hữu hạn hiệu quả cho mô hình 2S2P1D, Tiouajni và 2 cộng sự đã đề xuất một thuật toán dùng để xấp xỉ mô-đun đàn hồi phức của mô hình 2S2P1D bằng mô-đun đàn hồi phức của mô hình Kelvin-Voigt mở rộng trong miền tần số [31]. Thuật toán xấp xỉ được xây dựng để xác định n giá trị của độ cứng lò xo Ei với giả thiết rằng tỷ số giữa độ nhớt của cản nhớt i và độ cứng của lò xo Ei , hay còn được gọi là thời gian đặc trưng  i , là biết trước. Thuật toán này tương tự với phương pháp được gọi “colocation method” phát triển bởi Schapery năm 1961 [32] nhằm xấp xỉ hàm mô-đun dão bằng chuổi Prony ở chổ một tập các giá trị của thời gian đặc trưng (characteristic time) được lựa chọn trước.

Như được thảo luận ở các phần sau, việc lựa chọn các giá trị của tốc độ suy giảm (decay rates) để xấp xỉ một cách “sát sao” trong một khoảng rất rộng của thời gian đánh giá là rất khó khăn.3 Mục tiêu nghiên cứu Mục tiêu chung của nghiên cứu này là phát triển một công thức giải số nhằm tích hợp mô hình Huet-Sayegh và 2S2P1D vào trong các chương trình tính toán bằng phương pháp phần tử hữu hạn sao cho chi phí tính toán trong từng bước tính toán là ổn định trong suốt quá trình tính toán. Công thức số tìm kiếm tương tự công thức số áp dụng cho mô hình Maxwell mở rộng ở chổ sử dụng tổng các hàm số mũ. Vì vậy, công thức giải số cho mô hình đàn hồi – nhớt của Huet-Sayegh và 2S2P1D có thể được ứng dụng dễ dàng trong các phần mềm tính toán phần tử hửu hạn trong miền thời gian sẵn có để phân tích ứng xử của kết cấu áo đường. Phương pháp giải quyết vấn đề trong báo cáo này gồm 2 bước.

Bước một là xác định hàm mô-đun dão cho mô hình Huet-Sayegh và 2S2P1D và bước hai là xấp xỉ hàm mô-đun dão bằng tổng của các hàm số mũ với độ chính xác cao để xây dựng công thức quy nạp cho quan hệ ứng suất – biến dạng. Để có thể đạt được mục tiêu nêu trên, các mục tiêu cụ thể trong nghiên cứu này bao gồm:  Xây dựng hàm mô-đun dão bằng phương pháp xấp xỉ giải tích (AAR) cho mô hình Huet-Sayegh và 2S2P1D sử dụng phương pháp của Schapery và Park [30] và thiết lập thuật toán xác định mô-đun dão chính xác bằng phương pháp số (NER) sử dụng thuật toán được đề xuất bởi Heck [28];  Cải tiến phương pháp của Dombi [33] để xác định các tham số của tổng các hàm số mũ trong phép xấp xỉ các giá trị của mô-đun dão nhằm đảm bảo độ mức độ chính xác cao cho phép xấp xỉ (sai số tương đối ≤ 10 3 );  Xây dựng mối quan hệ quy nạp của ứng suất và biến dạng để có thể giảm chi phí tính toán trong phép tính tích phân xoắn; 3  Xác thực phương pháp đề xuất thông qua việc so sánh kết quả tính toán số thu được với các lời giải giải tích và các lời giải thu được bằng các phương pháp khác.4 Đối tượng và phạm vi nghiên cứu Mô hình đàn hồi – nhớt tuyến tính của Huet-Sayegh và 2S2P1D. Hình 1-2 Minh họa mô hình lưu biến của Huet-Sayegh (a) và 2S2P1D (b) 1.5 Cách tiếp cận và phương pháp nghiên cứu 1.1 Cách tiếp cận Lý thuyết kết hợp với thực nghiệm và mô phỏng bằng phương pháp số.2 Phương pháp nghiên cứu - Phân tích ứng xử đàn hồi – nhớt của bê tông nhựa. - Mô hình hóa qui luật ứng xử bằng mô hình Huet-Sayegh và 2S2P1D.

- Xây dựng công thức giải số cho mô hình ứng xử đàn hồi nhớt trong miền thời gian.6 Nội dung nghiên cứu - Nghiên cứu tổng quan về các phương pháp giải số mô hình đàn hồi – nhớt Huet- Sayegh và 2S2P1D. - Xác định các thông số của mô hình ứng xử vật liệu và xấp xỉ các giá trị của mô-đun dão bằng tổng các hàm số mũ. - Xây dựng công thức giải số cho mô hình ứng xử đàn hồi nhớt trong miền thời gian. - Báo cáo kết quả nghiên cứu và xuất bản 4 Chương 2: CƠ SỞ LÝ THUYẾT 2.1 Các phương trình cơ bản của cơ học vật rắn biến dạng đàn hồi – nhớt Gọi σ ( x ,t ) là ten-xơ ứng suất Cauchy, mô tả trạng thái ứng suất tại một điểm x và tại thời điểm t của một vật rắn bất kỳ ở trạng thái cân bằng dưới tác dụng của tải trọng và tác động bên ngoài.

Nếu không có ghi chú khác, qui ước tổng của Einstein được áp dụng trong suốt báo cáo này. Phương trình cân bằng cục bộ tại mọi điểm của vật rắn được biểu thị bằng phương trình (1): div (σ (x ,t ))   0f  0 (1) trong đó 0 là khối lượng riêng của vật rắn, f là lực khối tác dụng lên mọi điểm của vật  ij rắn, và toán tử div được định nghĩa  div(σ ) i . Gọi ε ( x , t ) là ten-xơ tổng biến x j dạng, ε m ( x , t ) là ten-xơ biến dạng do tác động cơ học và ε th ( x , t ) ten-xơ biến dạng do tác động của nhiệt độ tại điểm x và thời gian t. Trong khuôn khổ của giả thiết biến dạng bé, chúng ta có thể viết biến dạng tổng bằng tổng của biến dạng cơ học và biến dạng do nhiệt: ε ( x, t )  ε m ( x, t )  ε th ( x, t ).