Tổng quan về luận án

Sự phát triển vượt bậc của kỹ thuật xây dựng hiện đại đòi hỏi các giải pháp kết cấu tối ưu hóa khả năng chịu lực, tiết kiệm vật liệu và nâng cao độ an toàn. Kết cấu liên hợp thép - bê tông cốt thép đã khẳng định ưu thế vượt trội nhờ kết hợp đồng thời cường độ chịu nén cao của bê tông và khả năng chịu kéo, uốn linh hoạt của thép hình. Tuy nhiên, quy trình thiết kế kết cấu công trình truyền thống theo tiêu chuẩn hiện hành (như ASD, PD hay LRFD) vẫn duy trì phương pháp tiếp cận hai bước tách rời: phân tích nội lực bằng lý thuyết đàn hồi tuyến tính kết hợp nguyên lý cộng tác dụng, sau đó kiểm tra khả năng chịu lực của từng cấu kiện đơn lẻ qua hệ số chiều dài tính toán $K$. Cách tiếp cận này bộc lộ khoảng trống nghiên cứu sâu sắc (research gap): không phản ánh được quan hệ phi tuyến thực tế giữa tải trọng và chuyển vị, bỏ qua sự suy giảm độ cứng liên tục và khả năng phân phối lại nội lực khi vật liệu bước vào miền chảy dẻo, gây lãng phí vật liệu đáng kể.

Luận án tiến sĩ kỹ thuật "Phân tích dẻo kết cấu khung cột thép dầm liên hợp chịu tải trọng tĩnh" của nghiên cứu sinh Hoàng Hiếu Nghĩa (chuyên ngành Kỹ thuật Xây dựng Công trình Dân dụng và Công nghiệp, Mã số 62 58 02 08, Trường Đại học Kiến trúc Hà Nội, 2020) dưới sự hướng dẫn khoa học của PGS.TS. Vũ Quốc Anh và PGS.TS. Nghiêm Mạnh Hiến, là công trình tiên phong giải quyết toàn diện bài toán phân tích phi tuyến vật liệu của toàn hệ kết cấu. Công trình đặt ra 4 câu hỏi nghiên cứu và giả thuyết khoa học tương ứng:

  1. Câu hỏi 1: Làm thế nào để thiết lập quy luật suy giảm độ cứng và quan hệ phi tuyến mô men - độ cong ($M-\phi$) của tiết diện dầm liên hợp thép - bê tông chịu uốn âm và uốn dương bằng phương pháp giải tích chính xác? Giả thuyết 1: Mô hình biến dạng phẳng Bernoulli kết hợp hàm ứng suất phi tuyến của bê tông nén (Kent & Park 1973) và bê tông kéo (Vebo & Ghali 1977) cho phép mô tả liên tục quá trình thoái hóa độ cứng $EI_t$ từ đàn hồi đến chảy dẻo hoàn toàn.
  2. Câu hỏi 2: Có thể xây dựng một phương trình mặt chảy dẻo không gian phản ánh được cả trạng thái đàn hồi, giai đoạn chảy dẻo trung gian và phá hoại hoàn toàn cho cột thép hình chữ I chịu nén uốn hai phương ($P-M_z-M_y-\theta$) hay không? Giả thuyết 2: Việc giải tích hóa biểu đồ phân bố ứng suất pháp trên tiết diện chữ I với góc xoay trục trung hòa $\theta$ bất kỳ sẽ xác lập được quỹ đạo mặt giới hạn đàn hồi và mặt chảy dẻo chuyển tiếp.
  3. Câu hỏi 3: Làm sao khắc phục nhược điểm tốn kém tài nguyên tính toán của phương pháp vùng dẻo (chia thớ) nhưng vẫn vượt qua giới hạn cứng nhắc của phương pháp khớp dẻo tập trung ($l_h = 0$)? Giả thuyết 3: Thiết lập phần tử thanh dầm - cột đa điểm dẻo (Multi-Plastic-Point Beam-Column Element) với quy luật suy giảm độ cứng dọc trục dạng hàm đa thức bậc 3 ($EI(x)$) sẽ mô phỏng chính xác sự lan truyền dẻo dọc chiều dài cấu kiện.
  4. Câu hỏi 4: Thuật toán giải phi tuyến gia tăng nào đảm bảo hội tụ ổn định khi phân tích cơ cấu phá hoại dẻo của toàn hệ khung không gian? Giả thuyết 4: Kết hợp thuật giải Euler gia số với phương pháp Newton-Raphson cải tiến trên nền tảng chương trình SPH (ngôn ngữ Delphi XE8) sẽ xác định chính xác hệ số tải trọng giới hạn $\lambda_p$ và trình tự xuất hiện khớp dẻo.

Phạm vi nghiên cứu tập trung vào hệ dầm đơn giản, dầm liên tục liên hợp hai nhịp (đối sánh với thực nghiệm quốc tế của Ansourian 1981) và hệ khung phẳng, khung Portal chịu tải trọng tĩnh đơn điệu, mở ra bước tiến vượt bậc trong việc áp dụng thiết kế trực tiếp một bước (Direct Analysis Method).

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận sự tiến hóa của lý thuyết phân tích phi tuyến kết cấu qua hai trường phái chính: phương pháp vùng dẻo (Plastic Zone Method / Fiber Element Method) và phương pháp khớp dẻo (Plastic Hinge Method). Trường phái vùng dẻo được phát triển bởi các học giả hàng đầu như Ziemian (1993), Izzuddin & Elnashai (1993), Filip & Fabio (1997), Li & Lui (2000), Chen, Kim & Choi (2001), Alvarenga (2009). Phương pháp này chia nhỏ cấu kiện thành $n$ phân đoạn và chia tiết diện thành hàng trăm sợi (thớ) vi mô. Dù đạt độ chính xác cao nhất về cơ học vật rắn biến dạng, y văn khẳng định phương pháp vùng dẻo đòi hỏi khối lượng tài nguyên tính toán khổng lồ. Đúng như luận án trích dẫn nhận định kinh điển: "Phương pháp vùng dẻo là cách tính chính xác để kiểm tra các phương pháp phân tích khác, nhưng phương pháp này phức tạp và đòi hỏi thời gian phân tích lớn (gấp hàng trăm lần tính theo phương pháp khớp dẻo – theo Ziemian [34])." Thậm chí, khi mô hình hóa chi tiết dầm liên hợp bằng phần tử khối 3D C3D8R trên phần mềm Abaqus, Kim & Lee (2014) và Mutlu (2017) phải sử dụng tới 45.514 phần tử và 58.000 điểm nút chỉ cho một cấu kiện đơn lẻ, khiến việc áp dụng vào toàn hệ khung nhà cao tầng trở nên bất khả thi trong thực tế kỹ thuật.

Ở chiều đối cực, trường phái khớp dẻo tập trung ($l_h = 0$) do Orbison (1982), Liew & Chen (1995), White (1993), Chen & Sohal (1995), Chan & Chui (2000) phát triển đã đơn giản hóa bài toán bằng cách giả định khớp dẻo hình thành tức thời tại hai đầu mút phần tử. Tuy nhiên, tranh luận học thuật (academic debate) nảy sinh khi các mặt chảy dẻo chuẩn tắc như Orbison (1982) hay AISC-LRFD (1994) áp đặt trạng thái dẻo hoàn toàn cứng nhắc: $$\left(\frac{P}{P_y}\right)^2 + \left(\frac{M_z}{M_{pz}}\right)^2 + \left(\frac{M_y}{M_{py}}\right)^2 = 1$$ Luận án chỉ rõ nghịch lý cốt lõi: "Mặt hạn chế của các phương pháp phân tích khớp dẻo là sử dụng các mặt chảy dẻo hoàn toàn Orbison 1982, AISC-LRFD 1994…, nó chưa thể hiện được sự chảy dẻo lan truyền từ khi tiết diện ở trạng thái đàn hồi chuyển sang trạng thái chảy dẻo hoàn toàn, sự chảy dẻo không phụ thuộc vào biến dạng dẻo của tiết diện và không phụ thuộc vào hình dáng tiết diện khác nhau."

Tại Việt Nam, các nghiên cứu của Bùi Công Thành & Chu Việt Cường (2014) về khung liên hợp mới dừng lại ở mô hình khớp dẻo suy giảm tại liên kết nửa cứng dầm - cột; Đoàn Ngọc Tịnh Nghiêm (2014) phân tích khung thép phẳng bằng phương pháp đồng xoay; Nguyễn Văn Tú (2014) và Võ Thanh Lương (2009) khảo sát khớp dẻo trên khung thép. Chưa có công trình nào giải quyết trọn vẹn bài toán lan truyền dẻo xuyên suốt cả chiều dày tiết diện và chiều dài cấu kiện dầm liên hợp thép - bê tông. Luận án định vị chính xác vào khoảng trống này: tạo ra điểm giao thoa tối ưu (trade-off optimization) giữa độ chính xác nghiêm ngặt của phương pháp chia thớ và tốc độ xử lý siêu việt của phương pháp phần tử hữu hạn giải tích.

                 KHOẢNG TRỐNG NGHIÊN CỨU (RESEARCH GAP)
                 (Chưa có mô hình cân bằng giữa độ chính xác
                  và tốc độ cho khung dầm liên hợp đa tầng)

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp vào cơ học kết cấu phi tuyến hệ thống lý thuyết giải tích mở rộng, giải quyết triệt để các giới hạn của mô hình đàn dẻo cổ điển:

  1. Mở rộng lý thuyết tương tác dẻo đa trục: Xây dựng thành công phương trình mặt chảy dẻo trung gian và mặt phá hoại $P-M_z-M_y-\theta-\phi_p$ cho cột thép tiết diện chữ I chịu nén uốn lệch tâm hai phương. Khác với mặt chảy Orbison độc lập hình học, mô hình giải tích của luận án tích hợp trực tiếp tỷ lệ kích thước hình học bản cánh ($b_f, t_f$), bản bụng ($h_w, t_w$) và góc nghiêng của trục trung hòa $\theta$, cho phép dự báo sự suy giảm độ cứng tiếp tuyến $EI_t$ chính xác tại mọi mức độ dẻo hóa trung gian.
  2. Lý thuyết ứng xử uốn phi đàn hồi của tiết diện liên hợp: Tích hợp tường minh mô hình ứng xử phi tuyến của bê tông chịu nén có xét đến hiệu ứng kiềm chế nở hông của cốt thép đai theo Kent & Park (1973), mô hình bê tông chịu kéo có xét đến suy giảm ứng suất sau nứt của Vebo & Ghali (1977), cùng mô hình đàn dẻo lý tưởng của thép hình theo Eurocode 3 (2005). Luận án chứng minh quan hệ $M-\phi$ của dầm liên hợp là một đường cong trơn phi tuyến chuyển tiếp liên tục, bác bỏ giả định nhảy bậc từ $M_e$ sang $M_p$ trong các tính toán truyền thống.

Khung phân tích độc đáo

Khung phân tích của luận án thiết lập mối liên kết chặt chẽ giữa ứng xử biến dạng vi mô ở cấp độ vật liệu với phản ứng động học vĩ mô của toàn hệ khung:

  • Tính toán giải tích mô men kháng uốn dẻo: Đối với mô men dương ($M > 0$), trục trung hòa dẻo đi qua bản sàn bê tông hoặc cánh thép; đối với mô men âm ($M < 0$), vùng bê tông chịu kéo bị nứt, sức kháng uốn do dầm thép kết hợp cốt thép sàn chịu kéo đảm nhận: $$M_{pl.Rd} = \int_{A_c} \sigma_c (y) y , dA + \int_{A_s} \sigma_s (y) y , dA + \int_{A_{st}} \sigma_{st} (y) y , dA$$
  • Mặt giới hạn đàn hồi và chảy dẻo trung gian của cột I: Thiết lập hệ phương trình biên phân chia 3 trạng thái ứng xử: miền đàn hồi ($P-M_{ze0}-M_{ye0}$), miền chảy dẻo một phần và miền chảy dẻo hoàn toàn.
  • Điều kiện biên và phạm vi áp dụng: Khung phân tích áp dụng cho các cấu kiện không bị mất ổn định cục bộ bản cánh, bản bụng và bỏ qua biến dạng trượt, bảo đảm tương thích hoàn toàn cho bài toán phân tích phi tuyến bậc nhất vật liệu dưới tác dụng của tải trọng tĩnh không đổi dấu.

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

Thiết kế nghiên cứu

Luận án vận dụng triệt để thế giới quan thực chứng (Positivism Paradigm) kết hợp chủ nghĩa duy thực phê phán (Critical Realism), lấy cơ học giải tích làm nền tảng toán học và mô phỏng số làm công cụ thực thi. Thiết kế nghiên cứu đa tầng bao gồm:

  • Cấp độ 1 (Cross-sectional Level): Khảo sát trạng thái ứng suất - biến dạng trên tiết diện ngang dầm liên hợp và cột thép chữ I.
  • Cấp độ 2 (Element Level): Phát triển phần tử thanh dầm - cột 2D/3D đa điểm dẻo (Multi-Plastic-Point Element) gồm 2 nút quy ước ở hai đầu và $n$ điểm biến dạng dẻo phân bố dọc trục.
  • Cấp độ 3 (System Level): Thiết lập ma trận độ cứng toàn hệ kết cấu $[K]$ và giải phương trình cân bằng phi tuyến gia tăng.

Quy trình nghiên cứu rigorous

Để chuyển tải quy luật suy giảm độ cứng từ cấp độ tiết diện sang cấp độ phần tử, tác giả đã đề xuất một đột phá toán học: xấp xỉ quy luật biến thiên của độ cứng kháng uốn $EI_t(x)$ dọc theo chiều dài phần tử $L$ dưới dạng hàm đa thức bậc 3: $$EI(x) = a_0 + a_1 x + a_2 x^2 + a_3 x^3$$ Dựa trên nguyên lý năng lượng bù biến dạng $U^*$ kết hợp biểu thức nội lực uốn $M(x)$ và lực dọc $N(x)$, các hệ số của ma trận độ cứng phần tử thanh đa điểm dẻo $[k_{dp}]$ được xác định chính xác thông qua các tích phân dạng: $$f_{ij} = \int_{0}^{L} \frac{m_i(x) m_j(x)}{EI(x)} , dx$$ Quy trình này cho phép ma trận độ cứng tự động cập nhật mức độ lan truyền dẻo dọc trục mà không cần phải chia nhỏ cấu kiện thành hàng chục phần tử hữu hạn con.

Data và phân tích

Luận án thực hiện quá trình kiểm chuẩn chéo (cross-validation) thông qua việc xây dựng bộ phần mềm chuyên dụng SPH (Steel Plastic Hinge) lập trình trên ngôn ngữ Embarcadero Delphi XE8. Thuật toán phân tích cấu trúc gồm 3 module liên hoàn:

  1. Module 1: Tự động tạo dữ liệu quan hệ $M-\phi$ cho dầm liên hợp và mặt chảy $P-M_z-M_y-\theta$ cho cột thép.
  2. Module 2: Lắp ghép ma trận độ cứng tiếp tuyến toàn hệ $[K]$ và véc-tơ tải trọng nút quy đổi ${F}$.
  3. Module 3: Giải hệ phương trình phi tuyến bằng thuật toán Euler từng bước kết hợp vòng lặp hiệu chỉnh Newton-Raphson cải tiến với sai số khống chế nghiêm ngặt $\epsilon_{tol} \le 10^{-4}$.

Dữ liệu kiểm chứng bao gồm:

  • Thí nghiệm dầm đơn giản liên hợp khẩu độ $L = 5.0\text{ m}$, thép chữ I kết hợp bản sàn bê tông dày $h_c = 100\text{ mm}$, tải trọng tập trung giữa nhịp $P = 100\text{ kN}$.
  • Hai mẫu dầm liên hợp liên tục 2 nhịp đối xứng CTB1 và CTB2 chịu tải trọng tập trung giữa nhịp theo thực nghiệm quốc tế của Ansourian (1981) với kích thước nhịp $2 \times 4.5\text{ m}$.
  • Khung Portal thép - bê tông 1 tầng 1 nhịp chịu tải trọng ngang và tải trọng đứng.
  • Khung phẳng 3 tầng 2 nhịp chịu tải trọng phân bố đều trên dầm và tải trọng ngang tại các nút sàn.

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả mô phỏng và phân tích số trên chương trình SPH mang lại 5 phát hiện then chốt mang tính đột phá:

  1. Bản chất thoái hóa độ cứng tiếp tuyến: Biểu đồ quan hệ $M-\phi$ của dầm liên hợp minh chứng độ cứng tiếp tuyến $EI_t$ suy giảm phi tuyến mạnh mẽ ngay từ khi thớ bê tông ngoài cùng đạt biến dạng nén $\epsilon_c = 0.002$. Tại thời điểm hình thành khớp dẻo hoàn toàn, tỷ lệ độ cứng tiếp tuyến $EI_t/EI_{\max}$ của dầm suy giảm về dưới $5%$, phản ánh chính xác hiện tượng dẻo hóa từng phần mà các mô hình khớp dẻo lý tưởng bỏ qua hoàn toàn.
  2. Độ chính xác vượt trội so với thực nghiệm quốc tế: Khi kiểm chứng trên dầm liên tục hai nhịp CTB1 và CTB2 của Ansourian (1981), mô men dẻo $M_p$ và tải trọng phá hoại dẻo giới hạn $P_u$ tính từ chương trình SPH đạt sai số cực nhỏ: đối với CTB1 đạt $223.5\text{ kNm}$ (chênh lệch dưới $2.1%$ so với thực nghiệm $228\text{ kNm}$); đối với CTB2 đạt $276.8\text{ kNm}$ (chênh lệch $1.8%$ so với thực nghiệm $272\text{ kNm}$).
  3. Hiện tượng tái phân phối nội lực trong dầm liên tục: Trong dầm liên tục liên hợp, khớp dẻo đầu tiên luôn xuất hiện tại tiết diện gối tựa giữa (chịu mô men âm), sau đó độ cứng tại gối suy giảm, kích hoạt quá trình tái phân phối mô men sang nhịp cho đến khi hình thành khớp dẻo thứ hai tại giữa nhịp, tạo thành cơ cấu phá hoại dẻo hoàn chỉnh.
  4. Hiệu ứng góc xoay dẻo lên mặt chảy dẻo cột chữ I (W14x426): Mặt chảy dẻo trung gian $P-M_z-M_y-\theta-\phi_p$ biến thiên liên tục theo góc xoay dẻo $\phi_p$. Ở trạng thái đàn hồi ban đầu ($\phi_p = 0$), mặt giới hạn đàn hồi thu hẹp hơn $35-40%$ so với mặt phá hoại Orbison (1982), chứng minh rằng việc sử dụng trực tiếp mặt Orbison cho phân tích khớp dẻo thông thường đã đánh giá quá cao khả năng làm việc đàn hồi của cột.
  5. Trình tự xuất hiện khớp dẻo trên khung 3 tầng 2 nhịp: Luận án ghi nhận rõ ràng tiến trình dẻo hóa từng bước: khớp dẻo hình thành đầu tiên tại các đầu dầm liên hợp tầng 1 và tầng 2, tiếp tục lan truyền sang dầm tầng 3, và cuối cùng xuất hiện tại chân cột thép tầng 1 để tạo thành cơ cấu phá hoại tổng thể ở hệ số tải trọng giới hạn $\lambda_p = 1.42$.
Đối tượng kiểm chuẩn Mô hình lý thuyết của Luận án (SPH) Dữ liệu Thực nghiệm / Quốc tế Mức độ sai lệch (%)
Dầm liên tục CTB1 ($M_p$) $223.5\text{ kNm}$ $228.0\text{ kNm}$ (Ansourian 1981) $-1.97%$
Dầm liên tục CTB2 ($M_p$) $276.8\text{ kNm}$ $272.0\text{ kNm}$ (Ansourian 1981) $+1.76%$
Mặt chảy cột thép W14x426 Mặt giải tích $P-M_z-M_y-\theta$ Mô hình Orbison (1982) Phản ánh chính xác vùng trung gian
Mô hình số dầm liên hợp Phần tử dầm đa điểm dẻo (SPH) 45.514 phần tử C3D8R (Abaqus) Tiết kiệm $>98%$ thời gian tính

Implications đa chiều

  • Về mặt học thuật: Luận án đặt nền móng phương pháp luận cho việc phân tích phi tuyến giải tích kết cấu liên hợp, thu hẹp khoảng cách giữa lý thuyết cơ học môi trường liên tục phức tạp và bài toán kết cấu công trình thực dụng.
  • Về mặt phương pháp luận: Thuật toán phần tử thanh đa điểm dẻo với hàm độ cứng bậc 3 $EI(x)$ mở ra hướng nghiên cứu mới, có thể áp dụng rộng rãi cho các dạng kết cấu phức tạp khác như vòm thép liên hợp, kết cấu giàn không gian hoặc kết cấu bê tông cốt thép dự ứng lực.
  • Về mặt thực tiễn kỹ thuật: Luận án cung cấp công cụ phần mềm SPH độc lập, cho phép các kỹ sư kết cấu tại Việt Nam phân tích trực tiếp một bước (Direct Analysis) mà không phụ thuộc vào các phần mềm thương mại đắt tiền, tối ưu hóa từ $12-18%$ khối lượng thép hình trong các dự án nhà nhiều tầng.
  • Về mặt chính sách và tiêu chuẩn: Kết quả nghiên cứu là cơ sở khoa học tin cậy để Bộ Xây dựng tham khảo xây dựng và hoàn thiện Tiêu chuẩn Quốc gia (TCVN) về Thiết kế Kết cấu Liên hợp Thép - Bê tông, vốn đang là khoảng trống tiêu chuẩn lớn tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nghiên cứu (research limitations):

  1. Giới hạn về loại hình tải trọng: Nghiên cứu mới khảo sát tải trọng tĩnh đơn điệu tăng dần, chưa xét đến tác động của tải trọng đổi chiều có chu kỳ (như tải trọng động đất) hoặc tải trọng lặp gây mỏi kết cấu.
  2. Bỏ qua phi tuyến hình học bậc hai: Luận án tập trung chuyên sâu vào phi tuyến vật liệu (bậc nhất đàn dẻo), chưa xét đến hiệu ứng $P-\Delta$ và $P-\delta$ của kết cấu khi chuyển vị lớn.
  3. Giả thiết về mất ổn định cục bộ và trượt liên kết: Nghiên cứu giả định liên kết chống cắt giữa bản bê tông và dầm thép là liên kết hoàn toàn (full shear connection), bỏ qua hiện tượng trượt mềm tương đối tại mặt tiếp xúc và chưa tính đến hiện tượng oằn cục bộ bản cánh/bản bụng thép hình.

Hướng nghiên cứu phát triển tiếp theo bao gồm:

  • Tích hợp phi tuyến hình học (hàm ổn định bậc 2) vào phần tử thanh đa điểm dẻo.
  • Mở rộng mô hình dẻo hóa cho bài toán phân tích động phi tuyến (Nonlinear Time-History Analysis) chịu động đất và cháy nổ.
  • Mô hình hóa phi tuyến liên kết nửa cứng và độ trượt dọc trục của neo chống cắt (Shear Connectors).

Tác động và ảnh hưởng

Luận án mang lại tầm ảnh hưởng sâu rộng trên cả 4 phương diện:

  • Tác động học thuật: Công trình đã công bố 15 bài báo trên các tạp chí khoa học chuyên ngành trong nước, 02 bài báo gửi phản biện tại các tạp chí thuộc danh mục ISI/SCOPUS uy tín, cùng 4 báo cáo khoa học tại các hội thảo chuyên ngành toàn quốc.
  • Đổi mới ngành công nghiệp xây dựng: Cung cấp giải pháp tính toán giải tích giúp các tổng thầu thiết kế và thi công rút ngắn thời gian kiểm toán kết cấu khung liên hợp quy mô lớn từ nhiều ngày (nếu dùng Abaqus) xuống vài phút (bằng chương trình SPH).
  • Ảnh hưởng xã hội và môi trường: Việc khai thác hiệu quả vùng làm việc sau đàn hồi giúp giảm thiểu lượng tiêu thụ thép và xi măng, trực tiếp cắt giảm phát thải khí nhà kính $CO_2$ trong chuỗi sản xuất vật liệu xây dựng.
  • Giá trị quốc tế: Khung phân tích giải tích của luận án hoàn toàn tương thích với các tiêu chuẩn tiên tiến thế giới như Eurocode 4 (EN 1994-1-1) và AISC 360-16.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận bộ khung lý thuyết giải tích hoàn chỉnh về quan hệ $M-\phi$ và mặt chảy dẻo trung gian $P-M_z-M_y-\theta$ làm tiền đề phát triển các đề tài phân tích động lực học công trình.
  • Các chuyên gia và nhà nghiên cứu kết cấu: Sở hữu công thức toán học tường minh về phần tử thanh đa điểm dẻo với hàm độ cứng $EI(x)$ bậc 3 để mở rộng cho các bài toán tối ưu hóa kết cấu phi tuyến.
  • Kỹ sư R&D và thiết kế kết cấu công trình: Được trang bị phần mềm SPH và phương pháp phân tích tiên tiến một bước, giúp nâng cao năng lực cạnh tranh trong thiết kế các công trình công cộng, nhà siêu cao tầng khẩu độ lớn.
  • Cơ quan quản lý nhà nước và ban biên soạn tiêu chuẩn: Có đầy đủ luận cứ khoa học để hoàn thiện hệ thống Tiêu chuẩn Xây dựng Việt Nam về kết cấu liên hợp.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Đóng góp độc đáo nhất là việc thiết lập thành công phương trình mặt chảy dẻo trung gian và mặt phá hoại $P-M_z-M_y-\theta-\phi_p$ bằng phương pháp giải tích cho cột thép chữ I chịu nén uốn hai trục. Luận án đã mở rộng căn bản lý thuyết mặt chảy cổ điển của Orbison (1982) và AISC-LRFD (1994) bằng cách đưa góc xoay trục trung hòa dẻo $\theta$ và thông số biến dạng dẻo $\phi_p$ vào phương trình, phản ánh chính xác sự chuyển tiếp dẻo lan truyền trên tiết diện thực.

2. Đột phá phương pháp luận của luận án vượt trội hơn các nghiên cứu quốc tế trước đây như thế nào? So với phương pháp vùng dẻo chia thớ của Ziemian (1993) hay mô hình 3D Solid C3D8R của Kim & Lee (2014) tiêu tốn hàng chục ngàn phần tử hữu hạn, luận án đề xuất phần tử thanh đa điểm dẻo với hàm độ cứng uốn bậc 3 $EI(x)$. Phương pháp này giảm hơn $98%$ thời gian tính toán nhưng vẫn bảo toàn độ chính xác tương đương thực nghiệm quốc tế (sai số $<2.1%$ so với thí nghiệm Ansourian 1981).

3. Phát hiện bất ngờ nhất trong quá trình phân tích số là gì? Phát hiện bất ngờ nhất là tại thời điểm xuất hiện vết nứt kéo đầu tiên trong bản bê tông của dầm liên hợp, sự phân phối lại nội lực diễn ra với tốc độ phi tuyến tính cực nhanh, làm độ cứng tiếp tuyến $EI_t$ tại gối tựa dầm giảm đột ngột hơn $40%$, kích hoạt sự gia tăng mô men dương ở giữa nhịp sớm hơn nhiều so với dự báo của lý thuyết đàn dẻo lý tưởng truyền thống.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Có. Toàn bộ thuật giải tích phân giải tích, phương trình ma trận độ cứng phần tử $[k_{dp}]$, sơ đồ khối thuật toán phân tích dẻo và mã nguồn (source code) các module tính toán then chốt bằng ngôn ngữ Delphi XE8 đều được công khai chi tiết trong phần phụ lục luận án, bảo đảm khả năng tái lập và kiểm chứng độc lập hoàn hảo.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Lộ trình 10 năm mở rộng khung phân tích SPH theo 3 giai đoạn: (1) Tích hợp phi tuyến hình học chuyển vị lớn $P-\Delta$ trong 2 năm đầu; (2) Mở rộng sang bài toán phân tích động lực học phi tuyến kháng chấn và sụp đổ lũy tiến (progressive collapse) trong 5 năm tiếp theo; (3) Tích hợp trí tuệ nhân tạo và học máy để tối ưu hóa thiết kế kết cấu liên hợp theo độ tin cậy trong giai đoạn 10 năm.

Kết luận

Luận án tiến sĩ của tác giả Hoàng Hiếu Nghĩa đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp mang tính nền tảng:

  1. Xác lập phương pháp giải tích chính xác xây dựng đường cong quan hệ mô men - độ cong ($M-\phi$) của tiết diện dầm liên hợp thép - bê tông, tích hợp đầy đủ ứng xử phi tuyến của vật liệu bê tông chịu nén (Kent & Park 1973) và chịu kéo (Vebo & Ghali 1977).
  2. Xây dựng hệ phương trình giải tích tường minh 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 ($P-M_z-M_y-\theta$) cho cột thép hình chữ I chịu nén uốn hai trục.
  3. Phát triển phần tử thanh dầm - cột đa điểm dẻo đột phá với hàm suy giảm độ cứng dọc trục bậc 3 ($EI(x)$), dung hòa hoàn hảo giữa độ chính xác cơ học và tốc độ tính toán.
  4. Xây dựng hoàn chỉnh chương trình máy tính SPH trên nền Delphi XE8 tích hợp các thuật giải phi tuyến Newton-Raphson cải tiến, tự động hóa phân tích tải trọng giới hạn dẻo $\lambda_p$.
  5. Kiểm chứng độc lập và chính xác hóa mô hình lý thuyết qua các thực nghiệm kinh điển quốc tế của Ansourian (1981), đạt độ tin cậy tuyệt đối với sai số dưới $2.1%$.
  6. Định hình phương pháp thiết kế phân tích trực tiếp một bước cho kết cấu liên hợp tại Việt Nam, mở ra các hướng nghiên cứu sâu rộng về phân tích phi tuyến kết cấu trong tương lai.