Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Kết cấu khung liên hợp thép - bê tông là giải pháp kết cấu kết hợp giữa thép hình và bê tông cốt thép nhằm tận dụng các ưu điểm chịu lực của cả hai loại vật liệu, giúp nâng cao khả năng chịu lực, độ ổn định, khả năng kháng chấn và chống cháy so với kết cấu thép hoặc bê tông đơn thuần.

Trong thực tế tính toán thiết kế công trình, phương pháp truyền thống thường được áp dụng qua hai bước riêng biệt:

  1. Phân tích đàn hồi tuyến tính kết hợp nguyên lý cộng tác dụng để xác định nội lực và chuyển vị của toàn hệ.
  2. Kiểm tra khả năng chịu lực, ứng suất giới hạn và độ ổn định của từng cấu kiện đơn lẻ dựa trên hệ số chiều dài tính toán $K$.

Quy trình thiết kế hai bước này bộc lộ nhiều hạn chế: không phản ánh đúng quan hệ phi tuyến thực tế giữa tải trọng và chuyển vị, chưa mô tả được ứng xử phi tuyến của vật liệu khi kết cấu chuyển sang giai đoạn chảy dẻo mà chưa bị phá hủy, đồng thời dẫn đến việc sử dụng vật liệu chưa tối ưu.

Để khắc phục hạn chế trên, các phương pháp phân tích phi tuyến vật liệu đã được phát triển với hai hướng tiếp cận chính:

  • Phương pháp vùng dẻo (Plastic Zone Method / Phương pháp thớ): Rời rạc hóa cấu kiện thành nhiều phần tử con và chia nhỏ tiết diện thành các thớ để theo dõi quá trình phát triển dẻo. Phương pháp này có độ chính xác cao nhưng khối lượng tính toán rất lớn (thời gian tính toán gấp hàng trăm lần so với phương pháp khớp dẻo theo Ziemian), chỉ khả thi với các kết cấu đơn giản, khó áp dụng cho khung công trình thực tế.
  • Phương pháp khớp dẻo (Plastic Hinge Method): Giả thiết chiều dài vùng dẻo bằng 0 ($l_h = 0$), khớp dẻo hình thành tập trung tại hai đầu phần tử khi nội lực đạt đến mặt chảy dẻo hoàn toàn (như mặt chảy Orbison 1982 hoặc AISC-LRFD 1994). Mô hình này chuyển đột ngột từ trạng thái đàn hồi sang chảy dẻo hoàn toàn, bỏ qua quá trình chảy dẻo lan truyền từ từ trên tiết diện và dọc theo chiều dài phần tử, không xét đến ảnh hưởng của hình dạng tiết diện và góc xoay dẻo.

Tại Việt Nam, các nghiên cứu về kết cấu liên hợp chủ yếu tập trung vào cấu kiện riêng lẻ hoặc liên kết, chưa có tiêu chuẩn thiết kế riêng và còn thiếu các phương pháp phân tích dẻo tổng thể cho toàn hệ kết cấu khung cột thép dầm liên hợp có xét đến sự lan truyền biến dạng dẻo. Do đó, việc xây dựng phương pháp phân tích dẻo hiệu quả cho hệ kết cấu này là vấn đề cần thiết.

Mục tiêu nghiên cứu

Luận án xác định 4 mục tiêu nghiên cứu cụ thể:

  1. Xây dựng đường quan hệ mô men – độ cong ($M-\phi$) của tiết diện dầm liên hợp có xét đến sự làm việc dẻo của vật liệu nhằm phản ánh đúng ứng xử thực tế của dầm liên hợp chịu tải trọng.
  2. Xây dựng phương trình mặt giới hạn đàn hồi, mặt chảy dẻo trung gian và mặt chảy dẻo hoàn toàn (mặt phá hoại) của tiết diện cột thép chữ I chịu nén uốn hai phương ($P-M_z-M_y-\theta$) bằng phương pháp giải tích.
  3. Xây dựng phương pháp phần tử hữu hạn (PTHH) với phần tử thanh đa điểm dẻo để phân tích phi tuyến hệ kết cấu khung cột thép dầm liên hợp, có kể đến sự lan truyền vùng biến dạng dẻo dọc trục thanh và trên tiết diện.
  4. Xây dựng chương trình máy tính ứng dụng (SPH) để thực hiện phân tích phi tuyến hệ kết cấu khung cột thép dầm liên hợp dưới tác dụng của tải trọng tĩnh.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Phân tích phi tuyến kết cấu khung cột thép dầm liên hợp chịu tải trọng tĩnh có xét đến sự làm việc dẻo của vật liệu.
  • Phạm vi nghiên cứu:
    • Kết cấu dầm đơn lẻ, dầm liên tục và hệ kết cấu khung phẳng gồm cột thép hình chữ I và dầm liên hợp thép - bê tông cốt thép.
    • Tải trọng tác dụng là tải trọng tĩnh, tăng dần đơn điệu và không đảo chiều.
    • Mô hình vật liệu: Thép sử dụng mô hình đàn dẻo lý tưởng (không xét giai đoạn củng cố/tái bền); bê tông chịu nén theo mô hình Kent - Park (1973) và chịu kéo theo mô hình Vebo - Ghali (1977).
    • Không xét ảnh hưởng của biến dạng trượt cắt trong cấu kiện.
    • Không xét hiện tượng mất ổn định cục bộ của bản thành/bản cánh và mất ổn định ngoài mặt phẳng (uốn - xoắn).
    • Không xét phi tuyến hình học ($P-\Delta$, $P-\delta$) trong quá trình thiết lập phương trình cân bằng.

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

Các hướng nghiên cứu trước

Luận án đã tổng quan các công trình nghiên cứu trong và ngoài nước liên quan đến các khía cạnh chính:

  • Phương pháp phân tích và thiết kế kết cấu: Phương pháp thiết kế theo ứng suất cho phép (ASD), thiết kế dẻo (PD) và thiết kế theo hệ số tải trọng và sức kháng (LRFD) theo các tiêu chuẩn AISC-LRFD (1994), ANSI/AISC 360-10, Eurocode 3 (EC3), Eurocode 4 (EC4), ASCE.
  • Phân tích khớp dẻo và khớp dẻo hiệu chỉnh: Các nghiên cứu của Neal (1956), Hodge (1959), Orbison, Mc Guire và Abel (1982), Liew và Chen (1995), White (1993), Chen và Sohal (1995), Chen, Kim và Choi (1996). Tại Việt Nam có các công trình của Vũ Quốc Anh và Nghiêm Mạnh Hiến (2000, phần mềm KCW V1.0), Bùi Công Thành và Chu Việt Cường (2012), Đoàn Ngọc Tịnh Nghiêm (2016), Nguyễn Văn Tú, Võ Thanh Lương.
  • Phân tích vùng dẻo và phương pháp thớ (Fiber Element): Các công trình của Ziemian (1993), Chan và Chui (2000), Izzuddin và Elnashai (1993), Filip và Fabio (1995), Ngo Huu Cuong, Kim và Ryul (2008, 2010), Đoàn Ngọc Tịnh Nghiêm và Ngô Hữu Cường, Phạm Thị Thanh Thủy.
  • Mặt chảy dẻo của tiết diện: Mô hình mặt chảy dẻo trơn 3 chiều của Orbison (1982), mặt chảy dẻo phụ thuộc hình dạng tiết diện của Duan và Chen (1990), mặt chảy dẻo hai mặt phẳng của AISC-LRFD (1994), mô hình của Balling và Lyon (2000).
Phương pháp phân tích Đặc điểm mô hình Ưu điểm Nhược điểm
Khớp dẻo truyền thống Chiều dài vùng dẻo $l_h = 0$, khớp xuất hiện tại 2 đầu phần tử Tính toán nhanh, thuật toán đơn giản Bỏ qua sự lan truyền dẻo trên tiết diện và dọc trục; không xét góc xoay dẻo
Khớp dẻo hiệu chỉnh Dùng hàm suy giảm độ cứng tại 2 đầu cấu kiện Chính xác hơn khớp dẻo truyền thống Chủ yếu áp dụng cho khung thép phẳng, ít nghiên cứu cho dầm liên hợp
Vùng dẻo (Chia thớ) Chia cấu kiện thành $n$ phần tử, chia tiết diện thành nhiều thớ Phản ánh chính xác ứng xử phi đàn hồi Khối lượng tính toán rất lớn, thời gian phân tích kéo dài gấp hàng trăm lần
Đề xuất của luận án Phần tử đa điểm dẻo, hàm độ cứng bậc 3 $EI(x)$, 3 mặt chảy dẻo giải tích Phản ánh lan truyền dẻo liên tục, thời gian tính toán hợp lý Giới hạn trong bài toán phẳng, chưa xét phi tuyến hình học

Khoảng trống nghiên cứu luận án giải quyết

  1. Khắc phục sự gián đoạn của mô hình khớp dẻo truyền thống bằng việc xây dựng quan hệ $M-\phi$ giải tích liên tục cho dầm liên hợp và hệ thống mặt chảy dẻo 3 trạng thái (đàn hồi, trung gian, phá hoại) cho cột thép chữ I.
  2. Giải quyết mâu thuẫn giữa độ chính xác và khối lượng tính toán của phương pháp vùng dẻo thông qua việc đề xuất phần tử thanh dầm - cột đa điểm dẻo có hàm độ cứng chống uốn $EI(x)$ biến thiên liên tục dạng bậc 3 dọc theo trục phần tử.

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

Khung lý thuyết và mô hình vật liệu

  • Mô hình vật liệu thép: Mô hình ứng suất – biến dạng ($\sigma - \varepsilon$) đàn dẻo lý tưởng theo Eurocode 3, được xác định bởi môđun đàn hồi $E_s$ và giới hạn chảy $f_y$, bỏ qua vùng củng cố.
  • Mô hình vật liệu bê tông chịu nén: Mô hình Kent và Park (1973) có kể đến ảnh hưởng kiềm chế nở hông của cốt thép đai. Quan hệ ứng suất – biến dạng $\sigma_c - \varepsilon_c$ gồm nhánh parabol tăng cường độ đến giá trị $K \cdot f_c'$ và nhánh tuyến tính suy giảm đến biến dạng cực hạn $\varepsilon_u$: $$\sigma_c = K f_c' \left[ \frac{2\varepsilon_c}{\varepsilon_0} - \left(\frac{\varepsilon_c}{\varepsilon_0}\right)^2 \right]$$
  • Mô hình vật liệu bê tông chịu kéo: Mô hình của Vebo và Ghali (1977), mô tả sự làm việc đàn hồi tuyến tính trước khi nứt và suy giảm ứng suất sau khi nứt dựa trên cường độ chịu kéo $f_{ct}$ và môđun đàn hồi $E_c$.
  • Lý thuyết dầm: Áp dụng lý thuyết dầm Euler - Bernoulli (giả thiết tiết diện phẳng trước và sau khi biến dạng vẫn phẳng và vuông góc với trục trung hòa).

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

  • Phương pháp giải tích: Thiết lập phương trình tích phân ứng suất trên tiết diện dầm liên hợp để rút ra quan hệ mô men – độ cong ($M-\phi$) cho cả vùng mô men dương và âm. Thiết lập hệ phương trình giải tích xác định mặt giới hạn đàn hồi, mặt chảy dẻo trung gian và mặt chảy dẻo hoàn toàn của tiết diện cột chữ I chịu nén uốn hai phương ($P-M_z-M_y-\theta$).
  • Phương pháp phần tử hữu hạn: Phát triển phần tử thanh 2 nút với $n$ điểm kiểm tra dẻo bên trong. Sử dụng hàm xấp xỉ độ cứng chống uốn tiếp tuyến $EI(x)$ dạng đa thức bậc 3 để tích phân giải tích ma trận độ cứng phần tử và véc-tơ tải trọng nút quy đổi.
  • Thuật toán giải phi tuyến: Áp dụng phương pháp gia tải từng bước kết hợp thuật toán lặp Newton - Raphson và Newton - Raphson cải tiến để giải hệ phương trình cân bằng phi tuyến toàn hệ: $$[K_t] {\Delta U} = {\Delta R} + {R_{ub}}$$
  • Công cụ lập trình: Xây dựng phần mềm SPH (Steel Plastic Hinge) trên môi trường Delphi XE8.
  • Nguồn dữ liệu kiểm chứng: Các kết quả thực nghiệm dầm liên hợp liên tục (CTB1, CTB2) của Ansourian (1981), các bài toán khung portal và khung 3 tầng 2 nhịp của Li, Chen cùng các phần mềm thương mại (Abaqus với phần tử khối C3D8R).

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

Chương 1: Tổng quan về vấn đề nghiên cứu

Chương 1 trình bày tổng quan về kết cấu khung liên hợp thép - bê tông, xu hướng thiết kế kết cấu hiện đại, các mức độ phân tích phi tuyến (đàn hồi bậc 1, đàn hồi bậc 2, đàn dẻo bậc 1, đàn dẻo bậc 2), các mô hình vật liệu thép và bê tông, cùng các phương pháp tính toán kết cấu khung khi xuất hiện biến dạng dẻo.

Chương này phân tích chi tiết ưu nhược điểm của các dạng phần tử trong phương pháp phần tử hữu hạn:

  • Phần tử khớp dẻo lý tưởng tại hai đầu;
  • Phần tử có liên kết lò xo phi tuyến;
  • Phần tử có chiều dài vùng dẻo tại hai đầu;
  • Phần tử chia thớ dọc theo chiều dài;
  • Phần tử khối 3D (C3D8R trong Abaqus).

Từ việc đánh giá các mô hình mặt chảy dẻo hiện hành (Orbison, Duan - Chen, AISC-LRFD, Balling), chương 1 khẳng định sự cần thiết phải xây dựng một mô hình phần tử mới có khả năng mô tả sự lan truyền biến dạng dẻo liên tục trên tiết diện và dọc trục phần tử mà không làm tăng quá mức khối lượng tính toán.

Chương 2: Xây dựng quan hệ mô men - độ cong của tiết diện dầm liên hợp và mặt chảy dẻo của tiết diện cột thép

Chương 2 tập trung giải quyết bài toán cơ học ở cấp độ tiết diện cho dầm liên hợp và cột thép chữ I bằng phương pháp giải tích:

  1. Khảo sát quá trình chảy dẻo của tiết diện: Khảo sát sự phân bố ứng suất – biến dạng qua 3 giai đoạn: đàn hồi, đàn dẻo và chảy dẻo hoàn toàn.
  2. Xây dựng quan hệ $M-\phi$ của tiết diện dầm liên hợp:
    • Trường hợp mô men dương ($M > 0$): Bản sàn bê tông chịu nén, dầm thép chữ I chịu kéo và một phần chịu nén, cốt thép sàn chịu nén. Luận án phân chia các trường hợp vị trí trục trung hòa dẻo (nằm trong bản bê tông, nằm trong bản cánh trên hoặc nằm trong bản bụng dầm thép) để thiết lập công thức tính mô men chống uốn tương ứng với từng cấp độ cong $\phi$.
    • Trường hợp mô men âm ($M < 0$): Bản bê tông chịu kéo (bị nứt và suy giảm khả năng chịu lực theo mô hình Vebo - Ghali), cốt thép sàn dọc chịu kéo, dầm thép chịu nén ở bản cánh dưới và kéo ở bản cánh trên.
    • Kết quả giải tích được so sánh và kiểm chứng với phương pháp chia thớ trên tiết diện I400, cho thấy sự trùng khớp về đường cong quan hệ $M-\phi$.
  3. Xây dựng mặt chảy dẻo của tiết diện cột thép chữ I chịu nén uốn hai phương:
    • Thiết lập phương trình giải tích cho Mặt giới hạn đàn hồi $(P - M_{ze0} - M_{ye0} - \theta)$: Xác định thời điểm thớ ngoài cùng của tiết diện bắt đầu đạt giới hạn chảy $f_y$.
    • Thiết lập phương trình giải tích cho Mặt chảy dẻo hoàn toàn $(P - M_z - M_y - \theta)$: Xác định trạng thái toàn bộ tiết diện đạt ứng suất chảy dẻo $f_y$ (ứng với góc nghiêng của trục trung hòa $\theta$).
    • Thiết lập phương trình Mặt chảy dẻo trung gian: Mô tả sự suy giảm độ cứng tiếp tuyến và biểu diễn giai đoạn chuyển tiếp đàn dẻo của tiết diện cột thép W14x426 dưới tác dụng đồng thời của lực dọc $P$ và mô men uốn hai trục $M_z, M_y$.
+-------------------------------------------------------------------------+
|                  MẶT CHẢY DẺO CỘT THÉP CHỮ I (P-Mz-My-θ)                |
+-------------------------------------------------------------------------+
|  Trạng thái 1: Mặt giới hạn đàn hồi (Elastic Limit Surface)             |
|  - Thớ ngoài cùng bắt đầu đạt giới hạn chảy fy                          |
|  - Tiết diện làm việc hoàn toàn đàn hồi với độ cứng ban đầu EI0         |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|  Trạng thái 2: Mặt chảy dẻo trung gian (Intermediate Yield Surface)     |
|  - Vùng dẻo lan truyền từ thớ biên vào trục trung hòa                   |
|  - Độ cứng tiếp tuyến EIt suy giảm dần theo tỷ lệ diện tích chảy dẻo    |
+-------------------------------------------------------------------------+
+-------------------------------------------------------------------------+
|  Trạng thái 3: Mặt chảy dẻo hoàn toàn (Full Plastic Surface / Phá hoại) |
|  - Toàn bộ tiết diện đạt ứng suất dẻo fy                                |
|  - Hình thành khớp dẻo hoàn toàn, không tiếp nhận thêm mô men uốn       |
+-------------------------------------------------------------------------+

Chương 3: Phương pháp PTHH phân tích dẻo kết cấu khung cột thép dầm liên hợp xét đến sự lan truyền vùng biến dạng dẻo của phần tử

Chương 3 trình bày việc thiết lập công thức phần tử hữu hạn mới phục vụ phân tích hệ kết cấu:

  1. Mô hình phần tử thanh đa điểm dẻo: Đề xuất phần tử dầm và cột gồm 2 nút ở 2 đầu (mỗi nút 3 bậc tự do trong bài toán phẳng: $u, v, \theta$) và $n$ điểm kiểm tra trạng thái dẻo được bố trí dọc theo chiều dài phần tử $L$.
  2. Hàm biến thiên độ cứng chống uốn $EI(x)$ dạng bậc 3: Nhằm mô tả sự suy giảm độ cứng không đồng đều khi vùng dẻo phát triển dọc trục thanh, độ cứng uốn tiếp tuyến $EI(x)$ được xấp xỉ bằng đa thức bậc 3: $$EI(x) = a_0 + a_1 x + a_2 x^2 + a_3 x^3$$ Các hệ số $a_0, a_1, a_2, a_3$ được xác định thông qua giá trị độ cứng tiếp tuyến $EI_i$ tại các điểm kiểm tra dẻo $x_i$ (thu được từ quan hệ $M-\phi$ ở dầm hoặc mặt chảy dẻo ở cột).
  3. Thiết lập ma trận độ cứng phần tử $[k_e]$: Dựa trên nguyên lý năng lượng bù biến dạng $U^*$ và biểu thức độ cứng biến thiên, các hệ số của ma trận độ cứng tiếp tuyến phần tử $[k_{dp}]$ được tích phân giải tích trực tiếp, loại bỏ việc phải chia nhỏ thanh thành nhiều phần tử con.
  4. **Véc-tơ tải trọng nút quy đổi ${f}:$$ Thiết lập công thức xác định tải trọng nút quy đổi cho phần tử thanh có độ cứng biến thiên khi chịu tải trọng phân bố đều $q$ và tải trọng tập trung $P_y$.
  5. Phương trình cân bằng toàn hệ: Thiết lập ma trận độ cứng toàn hệ $[K]$ và véc-tơ tải trọng nút ${R}$ thông qua ma trận chuyển trục tọa độ $[T]$.

Chương 4: Xây dựng chương trình phân tích dẻo và khảo sát một số bài toán

Chương 4 trình bày thuật toán, cấu trúc chương trình SPH và tiến hành kiểm chứng qua 4 bài toán kết cấu cụ thể:

  1. Thuật toán chương trình SPH: Chương trình được cấu trúc thành các mô-đun: nhập dữ liệu, tính toán đặc trưng tiết diện ($M-\phi$ và mặt chảy dẻo), lắp ghép ma trận độ cứng, giải lặp Newton - Raphson cập nhật trạng thái dẻo, kiểm tra điều kiện hội tụ sai số lực/chuyển vị ($\varepsilon \le \varepsilon_{tol}$), xác định hệ số tải trọng giới hạn $\lambda_p$ và tỷ lệ chảy dẻo của các tiết diện.
  2. Bài toán 1 - Dầm đơn giản liên hợp thép - bê tông: Nhịp dầm chịu tải trọng tập trung $P = 100\text{ kN}$ tại giữa nhịp. Kết quả phân tích chỉ ra quan hệ tải trọng - chuyển vị, sự suy giảm độ cứng tiếp tuyến $EI_t/EI_{max}$ và quá trình hình thành khớp dẻo tại giữa nhịp.
  3. Bài toán 2 - Dầm liên tục liên hợp 2 nhịp (Mẫu CTB1 và CTB2 theo thí nghiệm Ansourian 1981):
    • So sánh giá trị tải trọng giới hạn và mô men dẻo $M_p$ tính từ chương trình SPH với số liệu đo đạc thực nghiệm của Ansourian.
    • Kết quả mô phỏng thể hiện chính xác thứ tự hình thành khớp dẻo: khớp dẻo thứ nhất xuất hiện tại gối tựa giữa (chịu mô men âm), sau đó phân phối lại nội lực và hình thành các khớp dẻo tiếp theo tại giữa nhịp cho đến khi hình thành cơ cấu phá hoại dẻo.
  4. Bài toán 3 - Khung Portal liên hợp 1 tầng 1 nhịp: Khung gồm 2 cột thép hình và xà ngang dầm liên hợp chịu tải trọng đứng và tải trọng ngang. Chương trình xác định được đường cong quan hệ tải trọng - chuyển vị ngang đỉnh khung, thứ tự xuất hiện các khớp dẻo tại chân cột và nách khung.
  5. Bài toán 4 - Khung phẳng liên hợp 3 tầng 2 nhịp (Mô hình Li): Khảo sát khung liên hợp chịu tải trọng phức hợp, đánh giá tỷ lệ chảy dẻo của từng tiết diện dầm/cột qua từng bước tăng tải, xác định hệ số tải trọng giới hạn $\lambda_p$ và cơ cấu phá hoại dẻo của toàn hệ.
Bài toán khảo sát Dạng kết cấu Loại tải trọng Dữ liệu đối chứng Kết quả kiểm chứng chính
Dầm đơn giản Dầm thép chữ I kết hợp bản bê tông Tải tập trung giữa nhịp Lời giải giải tích & chia thớ Xác định suy giảm $EI_t/EI_{max}$, khớp dẻo tại giữa nhịp
Dầm liên tục CTB1, CTB2 Dầm liên hợp 2 nhịp liên tục Tải tập trung tại các nhịp Thí nghiệm Ansourian (1981) Khớp dẻo xuất hiện tại gối giữa trước, sau đó đến nhịp; giá trị $M_p$ khớp thực nghiệm
Khung Portal 1 tầng 1 nhịp (cột thép, dầm liên hợp) Tải trọng đứng và tải ngang Nghiên cứu của các tác giả quốc tế Xác định chuyển vị ngang đỉnh, thứ tự hình thành khớp dẻo tại cột và dầm
Khung 3 tầng 2 nhịp Khung phẳng liên hợp nhiều tầng Tải trọng nút và tải phân bố Mô hình khung của Li Mô tả tỷ lệ chảy dẻo từng bước tải, xác định hệ số tải trọng giới hạn $\lambda_p$

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

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

  1. Xây dựng hoàn chỉnh phương pháp giải tích để xác định quan hệ mô men – độ cong ($M-\phi$) của tiết diện dầm liên hợp thép – bê tông cốt thép cho cả vùng mô men dương và mô men âm, có xét đến mô hình phi tuyến của bê tông chịu nén (Kent - Park) và bê tông chịu kéo (Vebo - Ghali).
  2. Thiết lập hệ phương trình giải tích xác định ba mặt trạng thái của tiết diện cột thép chữ I chịu nén uốn hai phương ($P-M_z-M_y-\theta$): mặt giới hạn đàn hồi, mặt chảy dẻo trung gian và mặt chảy dẻo hoàn toàn, cho phép mô tả liên tục quá trình phát triển dẻo trên tiết diện cột phụ thuộc vào góc xoay dẻo.
  3. Đề xuất mô hình phần tử thanh dầm - cột đa điểm dẻo với quy luật biến thiên độ cứng chống uốn $EI(x)$ dạng đa thức bậc 3 dọc theo trục thanh. Mô hình này tích phân giải tích ma trận độ cứng phần tử mà không cần chia nhỏ cấu kiện thành nhiều phần tử con, phản ánh được sự lan truyền biến dạng dẻo trong toàn hệ kết cấu.

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

  1. Xây dựng và hoàn thiện bộ chương trình tính toán SPH bằng ngôn ngữ Delphi XE8, cung cấp công cụ phân tích phi tuyến vật liệu cho kết cấu dầm và khung phẳng liên hợp thép - bê tông dưới tác dụng của tải trọng tĩnh.
  2. Chương trình cho phép xác định trực tiếp đường cong ứng xử tải trọng – chuyển vị, hệ số tải trọng giới hạn $\lambda_p$, tỷ lệ phần trăm chảy dẻo của các tiết diện, vị trí và thứ tự hình thành các khớp dẻo trong công trình.
  3. Giúp rút ngắn thời gian tính toán đáng kể so với việc sử dụng phương pháp phần tử khối 3D hoặc phương pháp chia thớ truyền thống, tạo cơ sở khoa học hỗ trợ công tác thiết kế phân tích trực tiếp kết cấu liên hợp.

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

Hạn chế của luận án

  • Mô hình phân tích mới được thiết lập và kiểm chứng cho hệ kết cấu phẳng (2D), chưa mở rộng cho hệ khung không gian 3D tổng quát.
  • Chưa xét đến ảnh hưởng đồng thời của phi tuyến hình học (hiệu ứng $P-\Delta$, $P-\delta$ và chuyển vị lớn).
  • Vật liệu thép chưa xét đến giai đoạn củng cố (tái bền) sau khi chảy dẻo; bỏ qua hiện tượng trượt tương đối tại mặt tiếp xúc giữa bản sàn bê tông và dầm thép (giả thiết liên kết chống cắt hoàn toàn).
  • Tải trọng khảo sát mới giới hạn ở tải trọng tĩnh tác dụng đơn điệu, chưa xét tải trọng lặp, tải trọng động hoặc động đất.

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

  • Phát triển phần tử thanh đa điểm dẻo cho bài toán khung không gian 3D chịu tải trọng phức hợp.
  • Tích hợp đồng thời phi tuyến vật liệu và phi tuyến hình học vào ma trận độ cứng của phần tử thanh đa điểm dẻo.
  • Mở rộng mô hình có xét đến sự làm việc mềm của liên kết dầm - cột (liên kết nửa cứng) và độ trượt dẻo của liên kết neo chống cắt thép - bê tông.
  • Nghiên cứu ứng xử của kết cấu khung liên hợp dưới tác dụng của tải trọng động lực học và tải trọng động đất có xét đến suy giảm độ cứng theo chu kỳ.

Giá trị tham khảo

Luận án là tài liệu tham khảo cho các đối tượng và lĩnh vực sau:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Kỹ thuật Xây dựng công trình dân dụng và công nghiệp: Cung cấp cơ sở lý thuyết về giải tích dẻo, phương pháp thiết lập quan hệ $M-\phi$ cho cấu kiện liên hợp và kỹ thuật xây dựng ma trận độ cứng phần tử hữu hạn phi tuyến với độ cứng biến thiên.
  • Kỹ sư kết cấu và chuyên gia thiết kế công trình: Tham khảo quy trình phân tích phi tuyến một bước, đánh giá khả năng chịu tải thực tế và cơ cấu phá hoại dẻo của hệ khung dầm liên hợp cột thép.
  • Các nhóm nghiên cứu phát triển phần mềm tính toán kết cấu: Sử dụng thuật toán phần tử thanh đa điểm dẻo, sơ đồ khối và mã nguồn chương trình SPH (trong phần phụ lục) để tích hợp vào các phần mềm phân tích kết cấu chuyên dụng.

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

1. Mô hình vật liệu thép và bê tông nào được sử dụng để xây dựng quan hệ $M-\phi$ trong luận án?

Luận án sử dụng mô hình ứng suất – biến dạng đàn dẻo lý tưởng theo Eurocode 3 cho vật liệu thép (không kể giai đoạn củng cố). Đối với bê tông, mô hình chịu nén của Kent và Park (1973) được sử dụng để kể đến hiệu ứng kiềm chế nở hông của cốt thép đai; mô hình chịu kéo của Vebo và Ghali (1977) được áp dụng để mô tả giai đoạn bê tông nứt và suy giảm ứng suất chịu kéo.

2. Mặt chảy dẻo của cột thép chữ I trong luận án có điểm gì khác biệt so với các mặt chảy dẻo truyền thống như Orbison hay AISC-LRFD?

Các mặt chảy dẻo truyền thống (Orbison 1982, AISC-LRFD 1994) chỉ mô tả trạng thái chảy dẻo hoàn toàn, coi tiết diện chuyển đột ngột từ đàn hồi sang dẻo hoàn toàn và không phụ thuộc vào góc xoay dẻo hay hình dạng tiết diện. Luận án đã xây dựng giải tích hệ thống gồm 3 mặt: mặt giới hạn đàn hồi, mặt chảy dẻo trung gian và mặt chảy dẻo hoàn toàn ($P-M_z-M_y-\theta$), giúp mô tả liên tục quá trình lan truyền dẻo từ thớ ngoài cùng vào trong tiết diện theo từng cấp tải trọng.

3. Phần tử thanh đa điểm dẻo trong luận án được thiết lập trên nguyên lý nào?

Phần tử thanh đa điểm dẻo gồm 2 nút ở hai đầu và $n$ điểm kiểm tra trạng thái dẻo bố trí dọc chiều dài phần tử. Thay vì chia nhỏ cấu kiện thành nhiều phần tử con, luận án sử dụng đa thức bậc 3 để biểu diễn quy luật biến thiên của độ cứng chống uốn tiếp tuyến $EI(x) = a_0 + a_1 x + a_2 x^2 + a_3 x^3$ dọc theo trục thanh. Từ đó, ma trận độ cứng tiếp tuyến của phần tử được tích phân giải tích trực tiếp dựa trên năng lượng bù biến dạng.

4. Chương trình SPH được xây dựng bằng công cụ nào và độ tin cậy được kiểm chứng ra sao?

Chương trình SPH (Steel Plastic Hinge) được lập trình bằng ngôn ngữ Delphi XE8. Độ tin cậy của chương trình được kiểm chứng thông qua việc so sánh kết quả tính toán với nghiệm giải tích, phương pháp chia thớ (Fiber Element), phần mềm Abaqus (sử dụng phần tử khối 3D C3D8R), và số liệu thí nghiệm thực tế trên dầm liên hợp liên tục CTB1, CTB2 của Ansourian (1981).

5. Vì sao phương pháp đề xuất lại tiết kiệm thời gian tính toán hơn phương pháp vùng dẻo (chia thớ)?

Phương pháp vùng dẻo đòi hỏi phải chia cấu kiện thành nhiều phần tử nhỏ dọc chiều dài và chia mỗi tiết diện thành hàng chục đến hàng trăm thớ nhỏ, dẫn đến số bậc tự do và số phép tính tăng vọt. Phương pháp của luận án giữ nguyên mô hình 1 cấu kiện là 1 phần tử (2 nút đầu mút) và tích phân giải tích ma trận độ cứng dựa trên hàm $EI(x)$ bậc 3, qua đó giảm khối lượng tính toán nhưng vẫn theo dõi được sự lan truyền dẻo dọc trục.

Kết luận

Luận án đã giải quyết bài toán phân tích dẻo cho kết cấu khung phẳng cột thép dầm liên hợp chịu tải trọng tĩnh bằng việc xây dựng quan hệ $M-\phi$ giải tích cho dầm liên hợp, hệ thống ba mặt chảy dẻo cho cột thép chữ I và mô hình phần tử thanh đa điểm dẻo với hàm độ cứng bậc 3. Chương trình tính toán SPH lập trình trên Delphi XE8 đã được kiểm chứng với các kết quả thực nghiệm và lý thuyết quốc tế, phản ánh chính xác quá trình lan truyền biến dạng dẻo và tải trọng giới hạn của kết cấu. Công trình đóng góp giải pháp tính toán hiệu quả, kết nối giữa phương pháp khớp dẻo đơn giản và phương pháp vùng dẻo phức tạp trong kỹ thuật kết cấu công trình.