Báo Cáo Nghiên Cứu: Phân Tích Phi Tuyến Dầm Thép – Bê Tông Liên Hợp Có Xét Đến Ảnh Hưởng Của Tương Tác Bán Phần Bằng Phương Pháp Vùng Dẻo

Tóm tắt nghiên cứu

Nghiên cứu tập trung giải quyết bài toán cơ bản nhưng phức tạp trong kỹ thuật công trình: Hành vi phi tuyến thực tế của dầm thép – bê tông liên hợp (Steel-Concrete Composite Beam – SCCB) diễn biến như thế nào khi đồng thời chịu ảnh hưởng của tương tác bán phần tại bề mặt tiếp xúc, sự lan truyền vùng dẻo trên tiết diện, ứng suất dư và hiệu ứng phi tuyến hình học?

Để trả lời câu hỏi này, tác giả đã phát triển phần tử dầm liên hợp không gian 8 bậc tự do (8-DOF) dựa trên nguyên lý thế năng toàn phần dừng kết hợp phương pháp Rayleigh – Ritz. Mô hình tích hợp đồng thời:

  • Mô hình phi tuyến vật liệu của bê tông (Karayannis), thép đàn – dẻo lý tưởng và quan hệ phi tuyến lực cắt – trượt của liên kết chống cắt (Ollgaard).
  • Kỹ thuật chia thớ vi phân (fiber approach) với 66 thớ thép và 64 thớ bê tông nhằm theo dõi sự lan truyền dẻo và dịch chuyển trục trung hòa của lõi đàn hồi.
  • Thuật toán phi tuyến gia tăng Euler kết hợp kỹ thuật điều khiển công hằng (Constant Work Control), được lập trình tự động hóa hoàn toàn bằng ngôn ngữ C++ thông qua chương trình SCCB.

Kết quả nghiên cứu chứng minh rằng hiện tượng trượt tại bề mặt tiếp xúc làm suy giảm đáng kể độ cứng tổng thể và khả năng chịu tải cực hạn ($P_u$) của dầm so với giả thiết tương tác toàn phần. Mô hình tính toán đã được thẩm tra và cho kết quả trùng khớp cao với các dữ liệu thực nghiệm kinh điển quốc tế (Chapman & Balakrishnan, Haremza, Ansourian), mở ra giải pháp mô phỏng số chính xác và tiệm cận hành vi thực của kết cấu liên hợp.

Bối cảnh và tầm quan trọng

Kết cấu liên hợp thép – bê tông (SCCB) là giải pháp kết cấu ưu việt được ứng dụng rộng rãi trong các công trình nhà cao tầng, cầu nhịp lớn và công trình ngầm nhờ tận dụng tối đa khả năng chịu nén của bê tông và khả năng chịu kéo vượt trội của thép. Hiệu quả làm việc đồng thời giữa hai vật liệu phụ thuộc trực tiếp vào các liên kết chống cắt (shear connectors) trên bề mặt tiếp xúc.

+-------------------------------------------------------------+
|                     BẢN SÀN BÊ TÔNG                         |
+-------------------------------------------------------------+
       ||                 ||                 ||    <-- Liên kết chống cắt
=============================================================
|                                                           |
|                       DẦM THÉP HÌNH                       |
|                                                           |
=============================================================

Tuy nhiên, trong thực tế thiết kế và tính toán kỹ thuật hiện nay tồn tại những khoảng trống nghiên cứu (research gaps) đáng kể:

  1. Giả thiết giản lược quá mức: Các tiêu chuẩn thiết kế hiện hành (như TCVN, Eurocode 4, AISC, BS 5400) thường giả thiết liên kết là hoàn hảo (tương tác toàn phần – full interaction) hoặc sử dụng các công thức quy đổi tiết diện tương đương tuyến tính. Cách tiếp cận này bỏ qua độ trượt cục bộ tại mặt tiếp xúc (partial interaction/slip), dẫn đến việc đánh giá sai lệch độ võng và độ cứng của hệ kết cấu.
  2. Hạn chế của mô hình khớp dẻo tập trung: Phần lớn các nghiên cứu trong nước trước đây chỉ xét phi tuyến vật liệu tại hai đầu phần tử theo mô hình khớp dẻo (plastic hinge) hoặc khớp dẻo hiệu chỉnh. Các mô hình này không thể mô phỏng được sự lan truyền dẻo dần dần (spread of plasticity) qua chiều cao mặt cắt ngang và dọc theo chiều dài dầm, đồng thời bỏ qua sự dịch chuyển của trục trung hòa lõi đàn hồi.
  3. Bỏ qua ứng suất dư: Quá trình cán nóng và hàn dầm thép hình luôn để lại ứng suất dư (residual stress) rất lớn, làm thay đổi thời điểm bắt đầu chảy dẻo của cấu kiện ngoài miền đàn hồi nhưng thường bị loại bỏ trong tính toán lý thuyết.

Xuất phát từ yêu cầu cấp thiết nhằm đánh giá chính xác khả năng chịu tải thực tế, tính năng làm việc của kết cấu ngoài miền đàn hồi và phòng ngừa sự cố sụp đổ công trình, nghiên cứu này mang tính thời sự cao, đặt nền móng cho việc áp dụng phân tích phi tuyến bậc cao vào thực tiễn kỹ thuật xây dựng tại Việt Nam.

Methodology và approach

Nghiên cứu tiếp cận bài toán bằng phương pháp phần tử hữu hạn phi tuyến nâng cao kết hợp giải thuật số hiện đại, với quy trình chặt chẽ:

1. Xây dựng phần tử hữu hạn 8 bậc tự do (8-DOF)

Phần tử dầm liên hợp được xây dựng dựa trên lý thuyết dầm Euler – Bernoulli có xét đến chuyển vị trượt tương đối dọc trục tại mặt tiếp xúc theo mô hình động học Newmark: $$\delta(x) = u_c(x) - u_s(x) + h_t v'(x)$$ Trong đó $u_c(x)$ và $u_s(x)$ lần lượt là chuyển vị dọc trục của bản bê tông và dầm thép, $h_t$ là khoảng cách giữa hai trục trọng tâm. Hàm chuyển vị dọc trục $u(x)$ được mô tả qua hàm dạng tuyến tính, trong khi độ võng $v(x)$ được biểu diễn bằng các hàm bậc ba Hermite.

2. Mô hình vật liệu phi tuyến và phân chia thớ vi phân

  • Bê tông: Áp dụng mô hình ứng suất – biến dạng phi tuyến của Karayannis (1994), bỏ qua khả năng chịu kéo: $$\sigma_c = f_c' \left[ 2\left(\frac{\varepsilon_c}{\varepsilon_0}\right) - \left(\frac{\varepsilon_c}{\varepsilon_0}\right)^2 \right]$$
  • Thép dầm và thép gia cường: Giả thiết quan hệ ứng suất – biến dạng đàn – dẻo tuyệt đối (không xét tái bền).
  • Liên kết chống cắt: Sử dụng quan hệ lực cắt – trượt thực nghiệm của Ollgaard (1971): $$P = P_{\max} (1 - e^{-\beta s})^\alpha$$
  • Ứng suất dư: Tích hợp các biểu đồ phân bố ứng suất dư tiêu chuẩn của Vogel (1985) và Lehigh Notes (1965).
  • Kỹ thuật chia thớ (Fiber Element): Mặt cắt ngang dầm thép được chia thành 66 phần tử diện tích nhỏ, bản bê tông được chia thành 64 phần tử vi phân. Kỹ thuật này cho phép tính toán chính xác trạng thái ứng suất – biến dạng của từng thớ sợi, xác định sự phát triển của vùng dẻo và tính toán tọa độ trục trung hòa dịch chuyển ($d_{zec}, d_{zes}$) ở từng bước gia tải.

3. Thiết lập phương trình cân bằng năng lượng

Ứng dụng ten xơ biến dạng Green trong hệ tọa độ Lagrangian để bao quát hiệu ứng phi tuyến hình học (chuyển vị lớn, biến dạng bé). Toàn bộ phương trình thế năng biến dạng tổng cộng $\Pi$ (gồm năng lượng biến dạng dầm thép, bê tông, cốt thép, năng lượng trượt liên kết và thế năng tải trọng ngoài) được thiết lập. Áp dụng điều kiện thế năng toàn phần dừng: $$\frac{\partial \Pi}{\partial d_i} = 0 \quad (i = 1 \dots 8)$$ Từ đó suy dẫn ra ma trận độ cứng cát tuyến $[K]$ và ma trận độ cứng tiếp tuyến gia tăng $[K_T]$.

4. Thuật toán phi tuyến và chương trình SCCB

Để giải hệ phương trình phi tuyến ngoài miền đàn hồi và vượt qua các điểm tới hạn (limit points), nghiên cứu kết hợp thuật toán gia tăng Euler với kỹ thuật điều khiển công hằng (Constant Work Control): $$\Delta W_n = {\Delta d_n}^T {\Delta r_n} = \text{const}$$ Kỹ thuật này tự động điều chỉnh hệ số gia tăng tải trọng $\Delta \lambda_n$ nghịch biến với độ cứng tiếp tuyến tức thời của kết cấu, đảm bảo sai số được phân bố đồng đều giữa các bước lặp. Toàn bộ thuật toán được tối ưu hóa và lập trình trên C++, hình thành phần mềm chuyên dụng SCCB có tốc độ xử lý vượt trội và khả năng đọc/xuất dữ liệu linh hoạt.

Phát hiện chính

Nghiên cứu đã thực hiện các khảo sát mô phỏng số chi tiết và đối chiếu trực tiếp với các bộ dữ liệu thí nghiệm thực tế nổi tiếng thế giới. Các phát hiện trọng tâm gồm:

Bộ dữ liệu thực nghiệm Cấu kiện khảo sát Hiện tượng phân tích Kết quả đối chứng từ SCCB
Chapman & Balakrishnan (1964) Dầm E1, U4 Tương tác toàn phần vs. Tương tác bán phần; Độ trượt tiếp xúc Đường cong $P - \Delta$ và phân bố trượt $s(x)$ trùng khớp xuất sắc với thực nghiệm.
Haremza (2009) Dầm A3 Các mức độ tương tác: Toàn phần, bán phần, không tương tác Mô phỏng chính xác sự sụt giảm độ cứng uốn khi mức độ liên kết giảm dần.
Ansourian (1981) Dầm CTB1, CTB4 Tải trọng tới hạn $P_u$, sự phát triển biến dạng dẻo sâu Sai số tải trọng tới hạn $P_u$ dự báo dưới 3% so với số liệu đo đạc thực tế. 
Tải trọng (P)
  ^
  |                     /--- Tương tác toàn phần (Full Interaction)
  |                   /
  |                 /---\  <-- Tương tác bán phần (Partial Interaction - Thực tế)
  |               /      \
  |             /         \
  |           /            \
  |         /               \--- Không tương tác (No Interaction)
  |       /
  |      /
  +---------------------------------------------> Chuyển vị / Độ võng

1. Suy giảm độ cứng và tải trọng tới hạn do tương tác bán phần

Khi số lượng chốt liên kết chống cắt giảm (tương tác bán phần), độ trượt tương đối tại mặt phân giới thép – bê tông tăng lên rõ rệt dọc theo chiều dài nhịp. Sự xuất hiện của biến dạng trượt làm giảm độ cứng chống uốn tổng thể của dầm, khiến độ võng tăng từ 20% đến 50% ở cùng một cấp tải trọng làm việc, đồng thời làm suy giảm tải trọng giới hạn cực hạn ($P_u$) so với tính toán theo lý thuyết tương tác toàn phần.

2. Sự dịch chuyển trục trung hòa và lan truyền vùng dẻo

Khác với giả thiết truyền thống coi trục trung hòa là cố định, phương pháp vùng dẻo chỉ ra rằng khi các thớ ngoài cùng của dầm thép hoặc bê tông bắt đầu chảy dẻo/nén vỡ, vùng đàn hồi bị thu hẹp dần. Trục trung hòa của phần lõi đàn hồi liên tục dịch chuyển lên phía trên hoặc xuống dưới tùy thuộc vào trạng thái chịu lực của tiết diện. Việc bỏ qua sự dịch chuyển này trong các phương pháp khớp dẻo truyền thống sẽ dẫn đến sai số lớn khi đánh giá khả năng tiêu tán năng lượng của kết cấu.

3. Tác động bất lợi của ứng suất dư

Ứng suất dư nén tại mép cánh dầm thép hình (theo mô hình Vogel/Lehigh) làm cho các thớ thép biên đạt tới giới hạn chảy dẻo $\sigma_y$ sớm hơn đáng kể so với dầm lý tưởng không có ứng suất dư. Điều này dẫn đến hiện tượng kết cấu chuyển sang giai đoạn làm việc phi tuyến sớm hơn, làm mềm kết cấu ở giai đoạn tải trọng sử dụng.

4. Tính ổn định và hội tụ của thuật toán điều khiển công hằng

Việc kết hợp thuật toán Euler với điều khiển công hằng ($\Delta W = \text{const}$) chứng minh tính ưu việt vượt trội trong việc theo dõi chính xác đường cong cân bằng sau điểm cực hạn (post-buckling/softening branch) mà không gặp hiện tượng phân kỳ toán học như các phương pháp gia tải theo lực thông thường.

Đóng góp khoa học

Đóng góp về mặt lý thuyết

  • Thiết lập hoàn chỉnh hệ phương trình thế năng toàn phần dừng phi đàn hồi cho phần tử dầm liên hợp 8 bậc tự do, đồng thời bao quát 4 yếu tố phi tuyến phức hợp: phi tuyến hình học bậc hai, phi tuyến vật liệu đa thành phần, phi tuyến liên kết tiếp xúc và ảnh hưởng của ứng suất dư ban đầu.
  • Cung cấp cơ sở giải tích chuẩn xác để mô tả quy luật dịch chuyển trục trung hòa của dầm liên hợp trong suốt tiến trình biến dạng dẻo.

Đóng góp về mặt phương pháp luận

  • Chuẩn hóa phương pháp phân tích vùng dẻo (Plastic Zone Method) thông qua việc phân chia lưới thớ vi phân (fiber section model), khắc phục triệt để các nhược điểm cố hữu của phương pháp khớp dẻo tập trung.
  • Ứng dụng thành công thuật toán điều khiển công hằng giúp tối ưu hóa số bước lặp và kiểm soát sai số tích lũy trong bài toán phi tuyến kết cấu phức tạp.

Đóng góp về mặt ứng dụng thực tiễn

  • Phát triển thành công phần mềm SCCB bằng ngôn ngữ C++ với hiệu năng xử lý cao, giao diện dữ liệu tệp văn bản (*.txt) tinh gọn, dễ dàng tích hợp vào các quy trình nghiên cứu, thiết kế và kiểm định kết cấu thực tế.
  • Cung cấp căn cứ khoa học định lượng giúp các kỹ sư thiết kế tối ưu hóa số lượng chốt chống cắt, dự báo chính xác độ võng thực tế và đánh giá độ an toàn kết cấu dầm liên hợp trong công trình dân dụng và cầu đường.

Đối tượng quan tâm

Nghiên cứu này mang lại giá trị thiết thực cho nhiều nhóm đối tượng trong ngành xây dựng và cơ học công trình:

  • Nhà nghiên cứu và Giảng viên chuyên ngành Xây dựng – Cầu đường: Nguồn tài liệu học thuật chuyên sâu về giải thuật phần tử hữu hạn phi tuyến, mô hình cơ học vi phân thớ và kỹ thuật giải lặp số.
  • Kỹ sư kết cấu và Chuyên gia tư vấn thiết kế: Công cụ hỗ trợ đánh giá ứng xử thực của dầm liên hợp khi thiết kế theo tiêu chuẩn Eurocode 4 hoặc AISC, đặc biệt trong các trường hợp bố trí liên kết chống cắt bán phần nhằm tiết kiệm chi phí thi công.
  • Các đơn vị thẩm tra và kiểm định an toàn công trình: Cơ sở tính toán tin cậy để đánh giá khả năng chịu tải còn lại của các công trình dầm sàn liên hợp hiện hữu chịu tải trọng nặng hoặc công trình có dấu hiệu hư hỏng liên kết bề mặt.
  • Cơ quan biên soạn quy chuẩn, tiêu chuẩn kỹ thuật: Dữ liệu tham khảo quan trọng phục vụ công tác hoàn thiện và cập nhật các tiêu chuẩn thiết kế kết cấu liên hợp tại Việt Nam theo hướng tiếp cận phân tích phi tuyến bậc cao.

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

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

Trả lời: Phát hiện then chốt là việc chỉ ra rằng độ trượt bán phần tại mặt tiếp xúc kết hợp với ứng suất dư làm suy giảm đáng kể độ cứng chống uốn và đẩy nhanh tiến trình chảy dẻo qua mặt cắt ngang. Các phương pháp tính toán tuyến tính hoặc khớp dẻo truyền thống đã đánh giá thấp độ võng thực tế và không phản ánh đúng cơ chế phá hoại của dầm liên hợp.

2. Phương pháp vùng dẻo (Plastic Zone Method) có gì vượt trội so với phương pháp khớp dẻo (Plastic Hinge Method)?

Trả lời: Phương pháp khớp dẻo chỉ giả định biến dạng dẻo tập trung cục bộ tại các nút đầu phần tử và coi phần còn lại hoàn toàn đàn hồi. Ngược lại, phương pháp vùng dẻo chia tiết diện thành hàng chục thớ nhỏ (fiber elements), cho phép theo dõi trực tiếp sự lan truyền dẻo từng điểm qua chiều cao tiết diện và dọc chiều dài cấu kiện, đồng thời tính toán chính xác sự dịch chuyển trục trung hòa của lõi đàn hồi.

3. Kết quả của chương trình SCCB có thể tổng quát hóa cho mọi loại dầm không?

Trả lời: Có. Thuật toán và phần tử 8-DOF của chương trình SCCB được xây dựng trên các nguyên lý cơ học tổng quát (Lagrangian, Rayleigh – Ritz), cho phép áp dụng hiệu quả cho nhiều loại sơ đồ kết cấu dầm đơn giản, dầm liên tục, dầm có tiết diện thép hình và bản bê tông khác nhau với các điều kiện biên đa dạng.

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

Trả lời: Hướng phát triển tự nhiên của đề tài bao gồm: mở rộng mô hình cho bài toán khung không gian 3D, xét đến ảnh hưởng của tải trọng động/động đất, hiện tượng từ biến – co ngót dài hạn của bê tông, hiệu ứng nứt cục bộ và ứng xử của kết cấu liên hợp trong điều kiện nhiệt độ cao (hỏa hoạn).

5. Ứng dụng thực tế của phần mềm SCCB là gì?

Trả lời: Phần mềm SCCB có thể được sử dụng trực tiếp để thẩm tra độ võng phi tuyến, xác định tải trọng phá hủy tới hạn của dầm thép – bê tông liên hợp trong các dự án cầu vượt, sàn liên hợp nhà cao tầng, giúp tối ưu hóa số lượng chốt chống cắt mà vẫn đảm bảo an toàn tuyệt đối.

Kết luận

Đề tài NCKH "Phân tích phi tuyến dầm thép – bê tông liên hợp có xét đến ảnh hưởng của tương tác bán phần bằng phương pháp vùng dẻo" do ThS. Lê Phương Bình chủ nhiệm đã hoàn thành xuất sắc các mục tiêu đề ra. Nghiên cứu đã xây dựng thành công phần tử hữu hạn 8 bậc tự do tiên tiến, làm chủ thuật toán điều khiển công hằng giải bài toán phi tuyến phức hợp và đóng gói hoàn chỉnh phần mềm chuyên dụng SCCB.

Kết quả nghiên cứu không chỉ làm sáng tỏ bản chất cơ học thực tế của kết cấu dầm liên hợp ngoài miền đàn hồi mà còn cung cấp một công cụ phân tích số đắc lực, mở ra tiềm năng ứng dụng to lớn trong việc tối ưu hóa thiết kế và nâng cao mức độ an toàn cho các công trình xây dựng hiện đại.