Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Phân tích phi tuyến khung thép bán cứng không gian chịu tải trọng tĩnh và động bằng phương pháp vùng dẻo" do nghiên cứu sinh Lê Văn Bình thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh (năm 2024), dưới sự hướng dẫn khoa học của PGS. Ngô Hữu Cường. Công trình thuộc 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ố: 62580208).

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

Trong kỹ thuật xây dựng hiện đại, kết cấu khung thép chịu lực được ứng dụng rộng rãi nhờ cường độ cao, trọng lượng nhẹ, tốc độ thi công nhanh và khả năng tiêu tán năng lượng khi chịu động đất thông qua sự làm việc dẻo dai ngoài miền đàn hồi cùng ứng xử vòng trễ của liên kết dầm - cột. Các tiêu chuẩn thiết kế hiện đại như Eurocode 3 và AISC/ANSI 360-22 đã chuyển dịch sang định hướng thiết kế trực tiếp (Direct Design / Performance-based Design), đòi hỏi các công cụ tính toán phi tuyến có khả năng mô phỏng sát với ứng xử thực tế của hệ kết cấu.

Phần lớn các nghiên cứu và phần mềm hiện hành thường áp dụng phương pháp khớp dẻo (plastic-hinge method) hoặc khớp dẻo hiệu chỉnh (refined plastic-hinge method) do ưu điểm tiết kiệm thời gian tính toán. Tuy nhiên, phương pháp khớp dẻo giả thiết sự chảy dẻo chỉ tập trung tại hai đầu phần tử, không phản ánh chính xác sự lan truyền dẻo dọc theo chiều dài cấu kiện và trên tiết diện ngang khi cấu kiện chịu lực phức tạp, chịu tải trọng phân bố hoặc có tác động của ứng suất dư. Ngược lại, phương pháp vùng dẻo (plastic-zone method) phân chia phần tử thành các điểm tích phân để theo dõi trạng thái ứng suất - biến dạng cục bộ, nhưng các nghiên cứu trước đây chủ yếu dừng lại ở bài toán phẳng, bài toán tĩnh hoặc chưa tích hợp đồng thời các yếu tố phi tuyến liên kết không gian, ứng xử vòng trễ và tải trọng động.

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

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

  1. Xây dựng một phần tử hữu hạn lai (hybrid finite element) mới cho bài toán phân tích phi tuyến khung thép bán cứng không gian chịu tải trọng tĩnh và động, có xét đầy đủ các yếu tố phi tuyến chủ đạo gồm: phi tuyến hình học (chuyển vị lớn), phi tuyến liên kết (độ mềm và ứng xử lặp trễ) và phi tuyến vật liệu (chảy dẻo từng phần, tái bền và ứng suất dư).
  2. Xây dựng chương trình phân tích số tin cậy (đặt tên là PZNASS) viết bằng ngôn ngữ MATLAB, áp dụng các thuật toán giải phi tuyến nâng cao cho hệ phương trình cân bằng tĩnh và phương trình chuyển động động lực học.

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

  • Đối tượng nghiên cứu: Ứng xử phi tuyến của hệ kết cấu khung thép không gian có liên kết bán cứng dưới tác dụng của tải trọng tĩnh và động.
  • Phạm vi nghiên cứu:
    • Kết cấu khung thép dạng phẳng và không gian ba chiều (3D).
    • Tải trọng tác dụng bao gồm tải trọng tĩnh đơn điệu, tải trọng động biến thiên theo thời gian, tải xung và tải trọng động đất (gia tốc nền).
    • Cấu kiện thanh được xây dựng dựa trên lý thuyết dầm Euler-Bernoulli (tiết diện phẳng trước và sau biến dạng, bỏ qua biến dạng trượt).
    • Đánh giá trạng thái dẻo của vật liệu theo tiêu chuẩn chảy dẻo von Mises, có xét ảnh hưởng đồng thời của ứng suất pháp và ứng suất tiếp.
    • Giả thiết biến dạng trong phần tử là bé nhưng chuyển vị và góc xoay của hệ kết cấu có thể rất lớn; không xét đến 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 ngang - xoắn.

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

Tác giả đã hệ thống hóa các công trình nghiên cứu trong và ngoài nước theo ba nhóm định hướng chính:

Các hướng nghiên cứu về phân tích phi tuyến tĩnh khung thép

  • Phương pháp khớp dẻo và khớp dẻo hiệu chỉnh: Lui và Chen (1987) khởi xướng việc đưa lò xo xoay vào mô phỏng liên kết bán cứng. Liew và White (1994, 1995) bổ sung sự chảy dẻo dần trên mặt cắt ngang. Hsieh và Deierlein (1991), Kim và Choi (2001), Nguyen và Kim (2014, 2016) mở rộng phương pháp khớp dẻo cho khung không gian 3D và hiện tượng mất ổn định ngang xoắn. Ngo-Huu cùng cộng sự (2012, 2013) phát triển phương pháp khớp dẻo thớ (fiber plastic hinge) kết hợp hàm ổn định. Doan-Ngoc cùng cộng sự (2019) sử dụng công thức đồng xoay kết hợp thuật toán chiều dài cung (arc-length).
  • Phương pháp vùng dẻo và phần tử hữu hạn: Pi và Trahair (1994, 1995) khảo sát dầm thành mỏng uốn hai phương và xoắn. Foley và Vinnakota (1997), Foley (2001) xây dựng phần tử hữu hạn phẳng theo dõi dẻo hóa qua tiết diện và chiều dài, kết hợp xử lý song song. Teh và Clarke (1998) áp dụng công thức đồng xoay và tiêu chuẩn von Mises. Jiang cùng cộng sự (2002) áp dụng công thức Lagrange cập nhật và kỹ thuật kết hợp phần tử lan truyền dẻo. Chirorean (2009) dùng mô hình Ramberg-Osgood cho mặt cắt. Da Silva cùng cộng sự (2015), Viana cùng cộng sự (2018) nghiên cứu phần tử phi tuyến có xét biến dạng trượt Timoshenko.

Các hướng nghiên cứu về phân tích phi tuyến động khung thép

  • Zhu cùng cộng sự (1995) đề xuất mô hình đường biên (bounding-line) 4 thông số cho liên kết chịu tải lặp.
  • Chan (1994), Chui và Chan (1997) phân tích phi tuyến hình học khung bán cứng sau mất ổn định và đánh giá khả năng tiêu tán năng lượng của liên kết.
  • Al-Bermani và Zhu (1995) áp dụng mô hình tái bền động học mặt biên kết hợp khớp dẻo.
  • Lui và Lopes (1997), Awkar và Lui (1999) khảo sát khung phẳng bán cứng chịu tải xung và động đất theo quy tắc CQC nhưng bỏ qua biến dạng dọc trục và hiệu ứng dỡ tải.
  • Sekulovic cùng cộng sự (2001, 2004) mô phỏng liên kết bán cứng bằng lò xo xoay kết hợp cản nhớt (dashpot) và phương pháp khớp dẻo hiệu chỉnh.
  • Kim cùng cộng sự (2001) dùng phương pháp lịch sử thời gian và mặt dẻo hiệu chỉnh Orbison.
  • Thai và Kim (2011), Nguyen và Kim (2014) mở rộng khớp dẻo thớ và phần tử liên kết 6 lò xo cho bài toán động.
  • Nguyen cùng cộng sự (2014) kiểm chứng phần tử chia thớ với thí nghiệm bàn lắc (shaking table test). Lemonis cùng cộng sự (2016) nghiên cứu cơ chế tiêu tán năng lượng qua vùng dẻo.
  • Các luận án trong nước liên quan: Vũ Quốc Anh (2008), Cao Văn Mão (2010), Nguyễn Hồng Sơn (2011), Nguyễn Quốc Hùng (2014), Nguyễn Hải Quang (2017), Đoàn Ngọc Tịnh Nghiêm (2018).
                            KHUNG TỔNG QUAN TÀI LIỆU
 [Phương pháp Khớp dẻo]       [Phương pháp Vùng dẻo]      [Phân tích Động lực học]
 - Lui & Chen (1987)          - Foley & Vinnakota (1997)   - Zhu et al. (1995)
 - Liew & White (1994, 1995)  - Teh & Clarke (1998)        - Chui & Chan (1997)
 - Kim & Choi (2001)          - Jiang et al. (2002)        - Kim et al. (2001)
 - Ngo-Huu et al. (2012, 2013)- Chirorean (2009)          - Nguyen & Kim (2014)
                       [Khoảng trống nghiên cứu xác định]
  Chưa có mô hình phần tử lai vùng dẻo 3D tích hợp đồng thời: Lagrange cập nhật,
  lò xo bán cứng 2 phương uốn, vòng trễ phi tuyến, thuật toán MGDCM & Newmark-NR.

Khoảng trống luận án giải quyết

Phần lớn các nghiên cứu đã công bố chỉ tập trung vào từng yếu tố phi tuyến riêng lẻ hoặc mô hình kết cấu phẳng. Các mô hình không gian 3D chủ yếu dựa trên giả thiết khớp dẻo đơn giản hóa. Luận án lựa chọn giải quyết khoảng trống: thiết lập phần tử hữu hạn lai không gian hoàn chỉnh theo phương pháp vùng dẻo, tích hợp liên kết bán cứng uốn 2 trục với mô hình ứng xử vòng trễ, giải quyết đồng thời bài toán phi tuyến tĩnh và phân tích lịch sử thời gian phi tuyến động lực học.

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

Cơ sở lý thuyết cơ học và biến dạng

  • Công thức Lagrange cập nhật (Updated Lagrangian): Cấu hình của kết cấu tại thời điểm $C_t$ được dùng làm trạng thái tham chiếu để xác định các đại lượng tại thời điểm tiếp theo $C_{t+\Delta t}$.
  • Phương trình công ảo tuyến tính hóa: $$\int_{^tV} {^tC_{ijkl}} , {t e{kl}} , \delta({t e{ij}}) , d(^tV) + \int_{^tV} {^t\tau_{ij}} , \delta({t\eta{ij}}) , d(^tV) = \delta(^{t+\Delta t}W_{ext}) - \int_{^tV} {^t\tau_{ij}} , \delta({t e{ij}}) , d(^tV)$$ Trong đó: $^tC_{ijkl}$ là ma trận vật liệu; $t e{ij}$ và $t\eta{ij}$ là thành phần tuyến tính và phi tuyến của ten-xơ biến dạng Green-Lagrange; $^t\tau_{ij}$ là ten-xơ ứng suất Cauchy.
  • Quan hệ biến dạng - chuyển vị dầm Euler-Bernoulli: Biến dạng dọc trục $\varepsilon_{xx}$ và hai thành phần biến dạng trượt $\gamma_{xy}, \gamma_{xz}$ tại điểm $N(y_N, z_N)$ được tính từ đạo hàm chuyển vị tâm uốn $u, v, w$ và góc xoắn dọc trục $\theta_x$, có xét thành phần phi tuyến bậc cao: $$\varepsilon_{xx} = u' - y v'' - z w'' + \frac{1}{2}(v'^2 + w'^2) + \frac{1}{2}(y^2 + z^2)\theta_x'^2 - y w'\theta_x' + z v'\theta_x'$$
  • Hàm dạng nội suy: Sử dụng hàm đa thức bậc nhất $H_1, H_2$ cho chuyển vị dọc trục $u$ và xoắn $\theta_x$; hàm Hermite bậc ba $H_3, H_4, H_5, H_6$ cho chuyển vị uốn $v$ và $w$.

Thiết lập ma trận phần tử và kỹ thuật tích phân số

  • Ma trận độ cứng tiếp tuyến phần tử dầm: $K_e = K_0 + K_L$, trong đó $K_0$ là ma trận độ cứng đàn hồi tuyến tính và $K_L$ là ma trận độ cứng hình học phi tuyến.
  • Tích phân số Newton-Cotes: Sử dụng $q$ điểm tích phân dọc theo trục phần tử và lưới $m \times n$ điểm tích phân trên mặt cắt ngang tiết diện. Việc phân bố các điểm nút lẻ trùng khớp với các mép biên và trục trung hòa giúp xác định chính xác quá trình lan truyền dẻo và gán ứng suất dư ban đầu (theo mẫu Vogel/ECCS, Lehigh Notes hoặc Kim 2002).
  • Tiêu chuẩn dẻo von Mises: Trạng thái chảy dẻo của từng điểm tích phân được kiểm tra qua ứng suất tương đương: $$\sigma_{eq} = \sqrt{\sigma_{xx}^2 + 3(\tau_{xy}^2 + \tau_{xz}^2)} \ge \sigma_y$$ Thuật toán quay lui đàn hồi Backward Euler được sử dụng để cập nhật ứng suất phi đàn hồi.

Thiết lập phần tử lai và kỹ thuật rút gọn tĩnh

Để mô phỏng liên kết bán cứng hai phương uốn tại hai đầu dầm, 4 lò xo xoay được ghép nối vào phần tử dầm 12 bậc tự do, tạo thành phần tử mở rộng 16 bậc tự do: $$\mathbf{d} = \left[ {d_a}{4\times 1}, {d_b}{12\times 1} \right]^T$$ Trong đó ${d_a} = [d_1, d_2, d_3, d_4]^T$ là các bậc tự do xoay cục bộ bên trong liên kết lò xo tại nút đầu và nút cuối; ${d_b}$ là 12 bậc tự do ngoài của phần tử dầm-cột.

Hệ phương trình độ cứng 16 bậc tự do: $$\begin{bmatrix} \mathbf{K}{aa} & \mathbf{K}{ab} \ \mathbf{K}{ba} & \mathbf{K}{bb} \end{bmatrix} \begin{Bmatrix} \mathbf{d}_a \ \mathbf{d}_b \end{Bmatrix} = \begin{Bmatrix} \mathbf{r}_a \ \mathbf{r}_b \end{Bmatrix}$$

Áp dụng kỹ thuật rút gọn tĩnh (static condensation) để triệt tiêu các bậc tự do bên trong ${d_a}$: $$\mathbf{d}a = \mathbf{K}{aa}^{-1}\mathbf{r}a - \mathbf{K}{aa}^{-1}\mathbf{K}{ab}\mathbf{d}b$$ Ma trận độ cứng tiếp tuyến rút gọn $\mathbf{K}'$ và véc-tơ nội lực rút gọn $\mathbf{r}'$ của phần tử lai 12 bậc tự do: $$\mathbf{K}' = \mathbf{K}{bb} - \mathbf{K}{ab}^T \mathbf{K}{aa}^{-1} \mathbf{K}{ab}$$ $$\mathbf{r}' = \mathbf{r}b - \mathbf{K}{ab}^T \mathbf{K}_{aa}^{-1} \mathbf{r}_a$$

Các mô hình liên kết bán cứng và ứng xử vòng trễ

  1. Mô hình tuyến tính: $M = R_{ki}\theta_r$, với $R_{ki} = \frac{\gamma}{1-\gamma}\frac{4EI}{L}$.
  2. Mô hình lũy thừa 3 thông số Kishi-Chen: $$M = \frac{R_{ki}\theta_r}{\left[1 + \left(\frac{\theta_r}{\theta_0}\right)^n\right]^{1/n}}, \quad \theta_0 = \frac{M_u}{R_{ki}}$$
  3. Mô hình hàm mũ Chen-Lui: $$M = M_0 + \sum_{j=1}^n C_j \left(1 - e^{-\frac{\theta_r}{2j\alpha}}\right) + R_{kf}\theta_r$$
  4. Mô hình 4 thông số Richard-Abbott: $$M = \frac{(R_{ki} - R_{kp})\theta_r}{\left[1 + \left|\frac{(R_{ki}-R_{kp})\theta_r}{M_0}\right|^n\right]^{1/n}} + R_{kp}\theta_r$$
  5. Mô hình vòng trễ động lực học: Luận án khảo sát ba phương pháp gồm: mô hình tái bền độc lập (independent hardening), mô hình tái bền động học (kinematic hardening) và mô hình mặt biên (bounding surface).
Thông số liên kết theo Chen-Lui Loại A (SWA) Loại B (TSA) Loại C (FEP) Loại D (EEP)
Kiểu cấu tạo liên kết Bản mã liên kết sườn Thép góc bụng và cánh Bản đệm phẳng đỉnh Bản đệm mở rộng
$M_0$ (kNm) 0 0 0 0
$R_{kf}$ (kNm/rad) 0 0 5,905 40,675
$\alpha$ $1{,}527 \times 10^{-4}$ $1{,}527 \times 10^{-4}$ $2{,}206 \times 10^{-4}$ $1{,}400 \times 10^{-4}$
$C_1$ (kNm) $-0{,}214 \times 10^3$ $-0{,}137 \times 10^3$ $-0{,}304 \times 10^3$ $-0{,}428 \times 10^3$
$C_2$ (kNm) $0{,}481 \times 10^3$ $0{,}320 \times 10^3$ $0{,}670 \times 10^3$ $0{,}912 \times 10^3$
$C_3$ (kNm) $-0{,}279 \times 10^3$ $-0{,}192 \times 10^3$ $-0{,}412 \times 10^3$ $-0{,}548 \times 10^3$
$C_4$ (kNm) $0{,}037 \times 10^3$ $0{,}028 \times 10^3$ $0{,}070 \times 10^3$ $0{,}088 \times 10^3$

Động lực học và thuật toán giải

  • Ma trận khối lượng: Kết hợp ma trận khối lượng tương thích $12 \times 12$ của phần tử $\mathbf{M}e$ và ma trận khối lượng tập trung $\mathbf{M}{lump}$: $\mathbf{M}_{sys} = \sum \mathbf{M}e + \mathbf{M}{lump}$.
  • Ma trận cản Rayleigh: $\mathbf{C}_T = \alpha_M \mathbf{M} + \beta_K \mathbf{K}$, với: $$\alpha_M = 2\xi \frac{\omega_1 \omega_2}{\omega_1 + \omega_2}, \quad \beta_K = 2\xi \frac{1}{\omega_1 + \omega_2}$$
  • Giải bài toán phi tuyến tĩnh: Phương pháp điều khiển chuyển vị tổng quát hiệu chỉnh (Modified Generalized Displacement Control Method - MGDCM) kết hợp thông số độ cứng tổng quát (GSP) và tiêu chuẩn hội tụ chuyển vị dư nhỏ nhất, cho phép vượt qua các điểm uốn giới hạn (snap-through) và điểm gập ngược (snap-back).
  • Giải bài toán phi tuyến động: Phương pháp tích phân từng bước Newmark kết hợp vòng lặp Newton-Raphson khử lực không cân bằng tại từng bước thời gian $\Delta t$.

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

Chương 1: Mở đầu

Chương này thiết lập nền tảng khoa học cho toàn bộ luận án:

  • Trình bày tổng quan về kết cấu thép, sự phát triển của tiêu chuẩn thiết kế từ ASD, PD đến LRFD và yêu cầu công cụ phân tích cho phương pháp thiết kế trực tiếp.
  • Phân tích chi tiết ưu điểm và hạn chế của hai phương pháp: phương pháp khớp dẻo (nhanh nhưng thiếu thông tin cục bộ và kém chính xác khi cấu kiện chịu lực phức tạp) và phương pháp vùng dẻo (chính xác, phản ánh đầy đủ quá trình lan truyền dẻo và ứng suất dư).
  • Điểm qua có hệ thống các công trình nghiên cứu phi tuyến tĩnh, động, liên kết bán cứng phẳng và không gian trong nước và quốc tế.
  • Xác định tính cấp thiết, mục tiêu, đối tượng, phạm vi nghiên cứu và phương pháp tiếp cận lý thuyết của đề tài.

Chương 2: Phần tử khung thép bán cứng không gian

Chương 2 trình bày toàn bộ quy trình thiết lập cơ sở toán học cho phần tử lai:

  • Xây dựng quan hệ biến dạng - chuyển vị và ma trận biến dạng $\mathbf{B} = \mathbf{B}0 + \mathbf{B}{NL}$ dựa trên công thức Lagrange cập nhật và lý thuyết Euler-Bernoulli.
  • Thiết lập công thức tích phân số Newton-Cotes xác định ma trận độ cứng đàn hồi $\mathbf{K}_0$, ma trận độ cứng hình học $\mathbf{K}_L$, véc-tơ nội lực $\mathbf{f}$ và tích hợp ứng suất dư.
  • Tích hợp 4 lò xo xoay vào hai đầu phần tử, thiết lập phương trình cân bằng 16 bậc tự do và áp dụng kỹ thuật rút gọn tĩnh để đưa về phần tử hữu hạn lai 12 bậc tự do.
  • Xây dựng ma trận khối lượng tương thích $\mathbf{M}_e$, ma trận cản Rayleigh $\mathbf{C}_T$ và ma trận chuyển trục tọa độ không gian tổng quát $\mathbf{T}_e$ qua các ma trận trung gian $\mathbf{B}_S$ và $\mathbf{R}$.
  • Trình bày chi tiết các mô hình liên kết bán cứng (tuyến tính, Kishi-Chen, Chen-Lui, Richard-Abbott), các mô hình vòng trễ (tái bền độc lập, động học, mặt biên) và thuật toán Backward Euler cập nhật trạng thái dẻo von Mises.

Chương 3: Chương trình phân tích PZNASS

Chương 3 trình bày kiến trúc và giải thuật của bộ phần mềm PZNASS được lập trình trên nền tảng MATLAB:

  • Cấu trúc chương trình phân tích hướng module gồm: module tiền xử lý (nhập thông số hình học, tiết diện, lưới tích phân, mô hình liên kết), module xử lý giải phi tuyến (tĩnh và động) và module hậu xử lý (xuất đường cong cân bằng, phân bố vùng dẻo).
  • Lưu đồ thuật toán phi tuyến tĩnh MGDCM kết hợp tiêu chuẩn hội tụ năng lượng và chuyển vị dư.
  • Lưu đồ thuật toán phi tuyến động Newmark kết hợp Newton-Raphson với việc cập nhật ma trận độ cứng hiệu dụng $\mathbf{K}{eff}$ và véc-tơ tải trọng hiệu dụng $\mathbf{F}{eff}$ tại mỗi bước thời gian.
  • Định nghĩa chi tiết cấu trúc mảng dữ liệu đầu vào và các hàm phân tích chính: hàm tính ma trận độ cứng phần tử lai, hàm kiểm tra trạng thái chảy dẻo từng điểm tích phân, hàm rút gọn tĩnh và hàm chuyển trục tọa độ.

Chương 4: Các ví dụ số kiểm chứng chương trình phân tích đã phát triển

Chương này thực hiện 15 ví dụ số đối chuẩn (benchmark) toàn diện nhằm kiểm chứng độ chính xác và hiệu quả của phần tử lai và chương trình PZNASS.

1. Phân tích phi tuyến tĩnh

  • Khung cổng chữ nhật Vogel: Kiểm tra khả năng mô phỏng liên kết cứng và bán cứng, đối chiếu với nghiệm chuẩn của Vogel; hiển thị biểu đồ phần trăm diện tích chảy dẻo trên tiết diện và dọc chiều dài.
  • Khung phẳng 1 nhịp 3 tầng: Đánh giá quan hệ lực - chuyển vị ngang đỉnh $P - \Delta$, khảo sát vị trí hình thành vùng dẻo tại chân cột và đầu dầm ở mức tải cực hạn $P_u$.
  • Khung phẳng Vogel 2 nhịp 6 tầng: So sánh hệ số tải giới hạn giữa liên kết cứng và các kiểu liên kết bán cứng loại A, B, C, D theo mô hình Chen-Lui.
  • Khung vòm không gian (Framed Dome): Kiểm tra tính ổn định phi tuyến hình học không gian, mô tả hiện tượng bật qua (snap-through) và ảnh hưởng của độ mềm liên kết đến đường cong cân bằng đàn hồi và phi đàn hồi.
  • Khung không gian 1 tầng, 2 tầng và khung Harisson: Khảo sát tương tác uốn hai phương kết hợp nén dọc trục và xoắn, lập biểu đồ phân bố vùng dẻo 3D tại trạng thái tải cực hạn $H_{max}, P_{max}$.
  • Khung không gian Orbison 6 tầng (Orbison's Six-storey Space Frame): Bài toán chuẩn kinh điển trong phân tích phi tuyến không gian. Kết quả PZNASS khớp với các công bố của Orbison, Liew, Jiang và Nguyen-Kim về hệ số tải giới hạn và tiến trình dẻo hóa các tầng.

2. Phân tích phi tuyến động lực học

  • Khung phẳng 1 tầng, 2 tầng chịu tải động: Khảo sát đáp ứng lịch sử thời gian chuyển vị ngang dưới tác dụng của tải xung chữ nhật (RC, NC), tải zigzag và gia tốc nền động đất El Centro (1940), Loma Prieta (1989), San Fernando.
  • Khung phẳng Vogel 2 nhịp 6 tầng chịu tải động: Đánh giá chu kỳ dao động riêng và vòng lặp trễ của các liên kết bán cứng tại các dầm biên.
  • Khung không gian 1 nhịp 2 tầng và khung Campbell: Phân tích đáp ứng 3D dưới tác động động đất, thể hiện lịch sử chuyển vị và biểu đồ chảy dẻo tại các thời điểm đạt đỉnh chuyển vị ($t = 3{,}48\text{ s}$ và $t = 8{,}39\text{ s}$).
  • Khung không gian bán cứng Orbison chịu tải động đất: Khảo sát dao động đàn hồi và phi đàn hồi, vẽ đồ thị vòng lặp trễ mô-men - góc xoay $(M - \theta_r)$ của các liên kết loại C, J, K dưới tác động của trận động đất El Centro và San Fernando.
                           BẢNG TỔNG HỢP CÁC VÍ DỤ SỐ KIỂM CHỨNG (CHƯƠNG 4)

Chương 5: Kết luận và kiến nghị

Tổng kết các kết quả đạt được về mặt lý thuyết, phương pháp số và khả năng ứng dụng thực tiễn của phần tử lai cùng phần mềm PZNASS; đồng thời chỉ ra các điểm còn giới hạn và định hướng các nghiên cứu tiếp theo.

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

Đóng góp về mặt lý luận và khoa học

  1. Thiết lập hoàn chỉnh phần tử hữu hạn lai không gian theo phương pháp vùng dẻo: Xây dựng thành công phần tử dầm không gian dựa trên công thức Lagrange cập nhật, tích hợp lò xo bán cứng hai phương uốn tại hai đầu và áp dụng kỹ thuật rút gọn tĩnh. Đây là mô hình phần tử vùng dẻo không gian đầy đủ xét đồng thời các yếu tố phi tuyến hình học, phi tuyến vật liệu và phi tuyến liên kết.
  2. Mô phỏng chính xác sự lan truyền dẻo và ứng suất dư: Sử dụng mạng lưới tích phân số Newton-Cotes theo ba phương $(q \times m \times n)$ kết hợp tiêu chuẩn von Mises, cho phép theo dõi chi tiết quá trình dẻo hóa từng phần trên tiết diện ngang và dọc theo chiều dài thanh mà không cần dùng hàm ổn định đơn giản hóa hay mặt dẻo giả định.
  3. Mô hình hóa ứng xử vòng trễ trong không gian 3D: Kết hợp các mô hình toán học của liên kết bán cứng (Kishi-Chen, Chen-Lui, Richard-Abbott) với các quy luật vòng trễ (tái bền độc lập, động học, mặt biên) trong phân tích động lực học khung không gian.

Đóng góp về mặt thực tiễn và kỹ thuật tính toán

  1. Phát triển chương trình tính toán PZNASS: Xây dựng thành công gói phần mềm bằng ngôn ngữ MATLAB, cung cấp công cụ phân tích có độ tin cậy cao phục vụ phương pháp thiết kế trực tiếp (Direct Design) theo tiêu chuẩn Eurocode 3 và AISC 360-22.
  2. Cung cấp biểu đồ trực quan về trạng thái giới hạn: Chương trình cho phép trích xuất biểu đồ tỷ lệ phần trăm diện tích chảy dẻo trên từng tiết diện và biểu đồ phân bố vùng dẻo toàn khung tại mọi bước tải tĩnh hoặc thời điểm kích động động đất, cung cấp cơ sở cho việc đánh giá hư hỏng và thiết kế dựa trên hiệu năng (Performance-based Design).
  3. Độ tin cậy được kiểm chứng: 15 ví dụ số đối chuẩn đã khẳng định kết quả từ PZNASS hoàn toàn phù hợp với các nghiệm giải tích, số liệu thực nghiệm và kết quả từ các phần mềm thương mại cũng như công bố quốc tế uy tín.

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

Hạn chế của nghiên cứu

Dựa trên các giả thiết và phạm vi nghiên cứu đã công bố trong luận án:

  • Phần tử được xây dựng dựa trên lý thuyết dầm Euler-Bernoulli, do đó bỏ qua ảnh hưởng của biến dạng trượt; điều này có thể làm giảm độ chính xác khi phân tích các cấu kiện có độ mảnh nhỏ (dầm chuyển, cột ngắn).
  • Chưa xét đến hiện tượng mất ổn định cục bộ của bản cánh, bản bụng thép hình và mất ổn định ngang - xoắn của cấu kiện uốn.
  • Giả thiết biến dạng trong phần tử là biến dạng bé (mặc dù chuyển vị và góc xoay của toàn hệ là lớn).
  • Khối lượng tính toán của phương pháp vùng dẻo lớn hơn nhiều so với phương pháp khớp dẻo, đòi hỏi cấu hình máy tính cao và thời gian xử lý dài đối với các công trình nhiều tầng, nhiều nhịp.

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

  • Mở rộng lý thuyết phần tử sang mô hình dầm Timoshenko hoặc dầm bậc cao để xét đến biến dạng trượt và hiện tượng vênh khi xoắn không đều.
  • Nghiên cứu tích hợp hiện tượng mất ổn định cục bộ và suy giảm độ cứng cục bộ vào mạng lưới điểm tích phân.
  • Tối ưu hóa thuật toán tính toán, ứng dụng kỹ thuật véc-tơ hóa và xử lý song song (parallel computing/GPU) nhằm rút ngắn thời gian phân tích cho các kết cấu quy mô lớn.
  • Mở rộng bài toán phân tích phi tuyến dưới tác động đồng thời của nhiệt độ cao (cháy nổ) và tải trọng cơ học.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu có giá trị học thuật và ứng dụng cao cho các đối tượng:

  • Nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Xây dựng / Cơ học Công trình: Cung cấp tài liệu hoàn chỉnh về lý thuyết phi tuyến kết cấu, công thức Lagrange cập nhật, phương pháp tích phân số Newton-Cotes và các thuật toán giải phương trình phi tuyến nâng cao (MGDCM, Newmark - Newton-Raphson). Chương 2 và Chương 3 là phần đáng tham khảo nhất về mặt toán học và thuật toán.
  • Giảng viên và nhà nghiên cứu: Cung cấp nguồn ví dụ số chuẩn (benchmark data) phong phú ở Chương 4 để kiểm chứng, so sánh khi phát triển các phần tử hữu hạn hoặc thuật toán mới.
  • Kỹ sư kết cấu và chuyên gia thiết kế công trình: Giúp hiểu rõ cơ chế ứng xử thực tế của liên kết dầm - cột bán cứng, cơ chế tiêu tán năng lượng khi chịu động đất và ứng dụng phân tích nâng cao phục vụ thiết kế trực tiếp theo tiêu chuẩn AISC 360-22 và Eurocode 3, cũng như đánh giá an toàn, gia cố sửa chữa công trình khung thép hiện hữu.

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

1. Sự khác biệt cốt lõi giữa phương pháp vùng dẻo trong luận án và phương pháp khớp dẻo truyền thống là gì?

Phương pháp khớp dẻo giả thiết toàn bộ biến dạng dẻo chỉ tập trung cục bộ tại hai đầu phần tử thông qua các khớp dẻo giả định, phần thân dầm hoàn toàn đàn hồi. Ngược lại, phương pháp vùng dẻo trong luận án chia phần tử thành nhiều điểm tích phân Newton-Cotes dọc chiều dài ($q$ điểm) và trên mặt cắt ngang ($m \times n$ điểm). Nhờ đó, phương pháp này theo dõi được chính xác sự phát triển và lan truyền dẻo từng phần trong toàn bộ thể tích cấu kiện, có xét đến ứng suất dư mà không cần giả định trước vị trí hình thành khớp dẻo.

2. Kỹ thuật rút gọn tĩnh (static condensation) được sử dụng để giải quyết vấn đề gì khi thiết lập phần tử lai?

Khi gắn 4 lò xo xoay phi tuyến vào hai đầu phần tử dầm không gian (2 lò xo uốn quanh trục $y$ và 2 lò xo uốn quanh trục $z$), phần tử mở rộng có 16 bậc tự do (gồm 12 bậc tự do chuyển vị ngoài và 4 bậc tự do góc xoay cục bộ bên trong liên kết lò xo). Kỹ thuật rút gọn tĩnh được áp dụng để khử 4 bậc tự do bên trong này, đưa phần tử lai về đúng 12 bậc tự do tương thích với phần tử dầm không gian chuẩn, giúp quá trình lắp ghép ma trận độ cứng tổng thể của toàn hệ kết cấu diễn ra thuận lợi.

3. Thuật toán MGDCM trong chương trình PZNASS có ưu điểm gì so với phương pháp Newton-Raphson tiêu chuẩn khi phân tích tĩnh?

Phương pháp Newton-Raphson tiêu chuẩn kiểm soát theo tải trọng thường bị phân kỳ hoặc mất ổn định số tại các điểm tới hạn (cực trị tải trọng). Thuật toán điều khiển chuyển vị tổng quát hiệu chỉnh (MGDCM) kết hợp thông số độ cứng tổng quát (GSP) và tiêu chuẩn chuyển vị dư nhỏ nhất cho phép chương trình tự động điều chỉnh tham số gia tải và dấu của bước tải, giúp bám sát đường cong cân bằng phi tuyến ngay cả qua các điểm uốn giới hạn (snap-through) và điểm gập ngược (snap-back).

4. Luận án sử dụng mô hình nào để đánh giá trạng thái chảy dẻo của vật liệu thép và xử lý quan hệ ứng suất - biến dạng?

Luận án áp dụng tiêu chuẩn chảy dẻo von Mises kết hợp với thuật toán quay lui đàn hồi Backward Euler cho vật liệu đẳng hướng. Trạng thái dẻo được đánh giá thông qua ứng suất tương đương có xét đồng thời ứng suất pháp $\sigma_{xx}$ và các ứng suất tiếp $\tau_{xy}, \tau_{xz}$. Các mô hình quan hệ ứng suất - biến dạng bao gồm mô hình đàn dẻo lý tưởng và mô hình đàn dẻo có xét tái bền động học, tái bền độc lập hoặc mô hình mặt biên khi chịu tải trọng đổi chiều.

5. Những kiểu liên kết bán cứng nào được khảo sát cụ thể theo mô hình hàm mũ Chen-Lui?

Luận án khảo sát 4 kiểu liên kết bán cứng tiêu chuẩn với các bộ thông số thực nghiệm xác định:

  • Loại A (Single Web Angle - SWA): Liên kết dùng bản mã liên kết sườn bụng dầm.
  • Loại B (Top and Seated Angle - TSA): Liên kết dùng thép góc ở bản cánh trên và dưới.
  • Loại C (Flush End Plate - FEP): Liên kết dùng bản đệm phẳng ngang đỉnh dầm.
  • Loại D (Extended End Plate - EEP): Liên kết dùng bản đệm mở rộng vượt qua cánh dầm.

Kết luận

Luận án của nghiên cứu sinh Lê Văn Bình đã xây dựng thành công phần tử hữu hạn lai không gian và chương trình tính toán PZNASS dựa trên phương pháp vùng dẻo, giải quyết đồng thời các bài toán phi tuyến hình học, phi tuyến vật liệu và phi tuyến liên kết bán cứng dưới tác dụng của tải trọng tĩnh và động. Công trình đã bổ sung cơ sở lý thuyết chặt chẽ cho việc mô phỏng ứng xử lan truyền dẻo và phản ứng vòng trễ tiêu tán năng lượng của khung thép 3D. Với 15 ví dụ số kiểm chứng toàn diện, luận án cung cấp một công cụ phân tích số đáng tin cậy, phục vụ đắc lực cho xu hướng thiết kế trực tiếp và thiết kế công trình dựa trên hiệu năng trong kỹ thuật xây dựng hiện đại.