Luận án tiến sĩ cơ học: Hệ thanh phẳng liên kết nửa cứng, vết nứt, ngẫu nhiên

Tổng quan về luận án

Trong khoa học cơ học kết cấu và kỹ thuật công trình hiện đại, việc đánh giá chính xác trạng thái ứng xử động lực học của hệ kết cấu chịu tải trọng điều hòa đóng vai trò quyết định đến độ an toàn và tuổi thọ công trình. Luận án tiến sĩ kỹ thuật chuyên ngành Cơ học vật thể rắn (Mã số: 62.01) của nghiên cứu sinh Dương Thế Hùng, dưới sự hướng dẫn khoa học của GS. Lê Xuân Huỳnh và PGS. Trần Văn Liên tại Trường Đại học Xây dựng Hà Nội (2010), mang tiêu đề "Phân tích hệ thanh phẳng có liên kết nửa cứng, vết nứt và có độ cứng, khối lượng phân bố ngẫu nhiên", là một công trình tiên phong giải quyết bài toán tương tác phức hợp giữa ba yếu tố phi lý tưởng hóa: tính chất liên kết nửa cứng tại nút, sự hiện diện của các tổn thương cục bộ (vết nứt), và tính bất định không gian của các thông số vật liệu cũng như hình học.

                  +-----------------------------------------------+
                  |     HỆ THANH PHẲNG PHỨC HỢP THỰC TẾ          |
                  +-----------------------+-----------------------+
                                          |
         +--------------------------------+-------------------------------+
         |                                |                               |
         v                                v                               v
+------------------+            +------------------+            +-------------------+
| LIÊN KẾT NỬA CỨNG|            | VẾT NỨT CỤC BỘ   |            | THAM SỐ NGẪU NHIÊN|
|  (Semi-Rigid)    |            | (Lò xo tương     |            | (Trường chuẩn 1D: |
| Lò xo cu, cv, cφ |            |  đương Ic(z), γ) |            |  EI(x), EA(x), m) |
+--------+---------+            +--------+---------+            +---------+---------+
         |                               |                                |
         +--------------------------------+-------------------------------+
                                          |
                                          v
                  +-----------------------------------------------+
                  | MÔ HÌNH HÓA PHẦN TỬ HỮU HẠN ĐỘ CỨNG ĐỘNG LỰC  |
                  |     PTL1 (Chuẩn) -> PTL2 (Lò xo) -> PTL3 (VC) |
                  +-----------------------+-----------------------+
                                          |
                                          v
                  +-----------------------------------------------+
                  | TÍCH PHÂN TRỌNG SỐ & KHAI TRIỂN NEUMANN       |
                  |     Nghịch đảo ma trận động lực ngẫu nhiên    |
                  +-----------------------+-----------------------+
                                          |
                                          v
                  +-----------------------------------------------+
                  | ĐẶC TRƯNG XÁC SUẤT & ĐÁNH GIÁ ĐỘ TIN CẬY      |
                  |  Kỳ vọng E[.], Phương sai Var[.], β, Pf (Khung)|
                  +-----------------------------------------------+

Khoảng trống học thuật (research gap) mà luận án định vị bắt nguồn từ việc các mô hình tính toán truyền thống luôn giả định liên kết dầm - cột là cứng tuyệt đối hoặc khớp lý tưởng, vật liệu có tính chất đồng nhất tiền định và kết cấu hoàn toàn lành lặn. Trên thực tế, các khảo sát thực nghiệm chứng minh rằng vùng liên kết luôn tồn tại chuyển vị xoay và trượt tương đối (tính nửa cứng), các khuyết tật vết nứt xuất hiện tất yếu trong quá trình thi công/vận hành, đồng thời độ cứng uốn $EI(x)$, độ cứng kéo nén $EA(x)$ cùng khối lượng phân bố $m(x)$ luôn dao động ngẫu nhiên theo không gian. Các nghiên cứu trước đây trên thế giới và tại Việt Nam chỉ khảo sát đơn lẻ từng yếu tố hoặc kết hợp rời rạc hai yếu tố (như vết nứt với dao động tiền định, hoặc tham số ngẫu nhiên trong dầm lý tưởng). Chưa có bất kỳ công trình nào thiết lập được khung giải tích - số tích hợp đồng thời cả ba yếu tố nêu trên trong phân tích động lực học ngẫu nhiên của hệ thanh phẳng.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học của luận án được xác lập chặt chẽ:

1. Câu hỏi nghiên cứu 1: Làm thế nào để xây dựng ma trận độ cứng động lực chính xác (Exact Dynamic Stiffness Matrix - MTĐCĐL) cho phần tử thanh phẳng khi xét đến đồng thời trường ngẫu nhiên một chiều của $EI(x)$, $EA(x)$, $m(x)$ và yếu tố cản nhớt nội tại?
2. Câu hỏi nghiên cứu 2: Cơ chế biến đổi ma trận độ cứng động lực từ phần tử cơ sở có liên kết cứng (PTL1) sang phần tử có liên kết nửa cứng/vết nứt (PTL2) và phần tử chứa vùng nút cứng vô cùng (PTL3) được thiết lập giải tích như thế nào thông qua các ma trận chuyển đổi $[B]$, $[KV_c]$ và $[H_e]$?
3. Câu hỏi nghiên cứu 3: Phương pháp đại số ngẫu nhiên nào cho phép nghịch đảo ma trận độ cứng động lực ngẫu nhiên tổng thể một cách hiệu quả để thu được biểu thức dạng đóng của giá trị kỳ vọng và phương sai của đáp ứng chuyển vị cũng như nội lực?
4. Giả thuyết khoa học 1: Việc mô hình hóa vết nứt và liên kết dầm - cột dưới dạng các phần tử lò xo đàn hồi tương đương kết hợp với kỹ thuật tích phân trọng số (Weighted Integrals) và chuỗi khai triển Neumann (Neumann Expansion) sẽ cho phép giải quyết trọn vẹn bài toán dao động điều hòa ngẫu nhiên mà không cần đến phương pháp mô phỏng Monte Carlo đắt đỏ về chi phí tính toán.
5. Giả thuyết khoa học 2: Sự kết hợp giữa tính linh hoạt của liên kết nửa cứng và độ suy giảm độ cứng cục bộ do vết nứt sẽ làm tăng đáng kể phương sai chuyển vị và suy giảm phi tuyến chỉ số độ tin cậy về bền cũng như độ cứng của hệ kết cấu khung phẳng.

Khung lý thuyết của công trình dựa trên sự hợp nhất giữa: Cơ học kết cấu giải tích (phương trình vi phân dao động thanh Euler - Bernoulli), Cơ học phá hủy đàn hồi tuyến tính (Linear Elastic Fracture Mechanics - LEFM quy đổi độ sâu vết nứt thành độ cứng lò xo tương đương), Lý thuyết trường ngẫu nhiên đồng nhất Gauss và Cơ học ngẫu nhiên (Stochastic Mechanics). Đóng góp đột phá của luận án mang tính định lượng rõ nét: phát triển thành công thuật toán và gói phần mềm chuyên dụng TK.mw trên nền tảng Maple 12, cho phép dự báo kỳ vọng và phương sai đáp ứng với sai số so với nghiệm giải tích tiền định và SAP2000 bằng $0%$ khi tham số biến thiên ngẫu nhiên $\varepsilon = 0$, đồng thời đánh giá định lượng xác suất phá hoại ($P_f$) của khung phẳng nhiều tầng nhiều nhịp dưới các dải tần số kích thích $\omega \in [10, 100] \text{ rad/s}$.

Phạm vi nghiên cứu bao gồm hệ kết cấu dầm và khung phẳng chịu tải trọng tĩnh hoặc tải trọng động điều hòa dừng $p(t) = p_0 e^{j\omega t}$. Các tham số đầu vào ngẫu nhiên $EI(x), EA(x), m(x)$ là các trường chuẩn một chiều dừng, chứa các tham số bé $\varepsilon_i \ll 1$. Kết quả của luận án tạo tiền đề khoa học vững chắc cho công tác chẩn đoán tổn hại kết cấu (Damage Identification) và thiết kế kết cấu công trình dựa trên độ tin cậy (Reliability-Based Design).

Literature Review và Positioning

Khung tổng quan y văn của luận án được cấu trúc thông qua ba dòng nghiên cứu chính (research streams) trong cơ học công trình quốc tế và trong nước:

=========================================================================================
                                TIẾN TRÌNH PHÁT TRIỂN Y VĂN
=========================================================================================
[Dòng 1: Cơ học ngẫu nhiên & Tích phân trọng số]
   Bolotin (1965-1984) -> Nakagiri (1987), Brenner (1991), Shinozuka & Yamazaki (1988)
   -> Phát triển ma trận độ cứng động lực ngẫu nhiên (Deodatis, 1991; Takada, 1990)
                                          │
[Dòng 2: Mô hình hóa vết nứt bằng lò xo tương đương]
   Freund (1976), Tada (1973) -> Dimarogonas (1983), Chondros (1989) -> Nguyễn Tiến Khiêm,
   Trần Văn Liên (2001-2008): MTĐCĐL cho dầm nhiều vết nứt
                                          │
[Dòng 3: Kết cấu khung có liên kết nửa cứng]
   Soares Filho & Wang (1997-2004), Ahmed L. (1996) -> Nguyễn Trâm (1995), Vũ Quốc Anh (2000),
   Cao Văn Mão (2002), Nguyễn Hồng Sơn (2006)
                                          │
                                          ▼
[KHOẢNG TRỐNG KHOA HỌC & VỊ TRÍ ĐỘT PHÁ CỦA LUẬN ÁN DƯƠNG THẾ HÙNG (2010)]
   Hợp nhất ĐỒNG THỜI: 
   Liên kết nửa cứng + Vết nứt + Trường độ cứng/khối lượng ngẫu nhiên EI(x), EA(x), m(x)
=========================================================================================

Dòng nghiên cứu 1: Động lực học kết cấu ngẫu nhiên và phương pháp tích phân trọng số. Khởi nguồn từ mô hình cơ bản của Viện sĩ V. V. Bolotin [87, 88] về hệ cơ học ngẫu nhiên dưới dạng phương trình toán tử: $$Ax = d$$ trong đó $x \in X$ là biến trạng thái (chuyển vị, ứng lực), $d \in D$ là tác động bên ngoài và $A$ là toán tử đặc trưng hệ thống. Các tác giả quốc tế như Nakagiri (1987) [75], Brenner (1991) [50], Benaroya & Rehak (1988) [49], Shinozuka & Yamazaki (1988) [79], và Der Kiureghian et al. (1991) [60] đã phát triển phương pháp phần tử hữu hạn ngẫu nhiên (Stochastic FEM). Tiếp nối, các công trình của Takada (1990) và Deodatis (1991) [81, 82] đã đưa ra khái niệm tích phân trọng số (Weighted Integrals) để biểu diễn các biến thiên không gian của độ cứng và khối lượng dầm chịu uốn, khắc phục hiện tượng phụ thuộc vào lưới chia phần tử (mesh dependency) của phương pháp ngẫu nhiên cổ điển.

Dòng nghiên cứu 2: Cơ học dầm có vết nứt và mô hình lò xo tương đương. Nghiên cứu về hư hỏng cục bộ dạng nứt được thiết lập bởi Freund & Dimarogonas [47] khi đề xuất mô hình lò xo tương đương. Tada et al. [67] đã áp dụng Cơ học phá hủy để tìm mối liên hệ tường minh giữa độ sâu vết nứt và độ cứng đàn hồi của lò xo. Tại Việt Nam, hướng nghiên cứu này phát triển mạnh mẽ qua các công trình của Nguyễn Xuân Hùng [20, 21], Nguyễn Tiến Khiêm và Trần Văn Liên [23, 27, 28, 68, 69]. Các tác giả đã áp dụng phương pháp ma trận độ cứng động lực và ma trận chuyển tiếp để khảo sát dầm và khung đàn hồi có nhiều vết nứt chịu tải trọng điều hòa và tải trọng sóng, mở ra hướng ứng dụng trong bảo trì công trình biển.

Dòng nghiên cứu 3: Phân tích kết cấu khung có liên kết nửa cứng. Tính nửa cứng của mối nối dầm - cột bằng thép hoặc bê tông cốt thép đã được chứng minh qua thực nghiệm và mô phỏng bởi các tác giả quốc tế như Wang C. & Soares Filho [74], Ahmed L. [61]. Tại Việt Nam, nhiều học giả đã nghiên cứu sâu về chủ đề này: Nguyễn Trâm [90] thiết lập ma trận đặc trưng cho phần tử có đầy đủ 6 thành phần độ mềm liên kết; Vũ Quốc Anh [1] xây dựng ma trận độ cứng, ma trận khối lượng và ma trận hình học xét vùng cứng nút nhưng bỏ qua biến dạng trượt; Cao Văn Mão [30] phân tích khung thép phi tuyến có liên kết nửa cứng và khớp dẻo; Nguyễn Hồng Sơn [41] phát triển siêu phần tử thanh có kể đến biến dạng trượt và dẻo hóa vật liệu.

Tranh luận và mâu thuẫn học thuật (Contradictions & Debates): Tồn tại hai quan điểm đối nghịch trong y văn:

1. Trường phái tất định cổ điển (Deterministic Paradigm): Cho rằng tính biến thiên của vật liệu và dung sai kích thước có thể được bao hàm hoàn toàn bằng các hệ số an toàn quy phạm tổng quát (Safety Factors), việc mô hình hóa ngẫu nhiên chỉ làm tăng độ phức tạp thuật toán mà không mang lại giá trị gia tăng tương xứng.
2. Trường phái cơ học ngẫu nhiên hiện đại (Stochastic Mechanics Paradigm): Chứng minh rằng trong miền tần số cao hoặc khi kết cấu tiệm cận trạng thái cộng hưởng, độ nhạy của đáp ứng đối với các dao động ngẫu nhiên phân bố không gian là cực lớn; việc sử dụng hệ số an toàn tiền định có thể dẫn đến đánh giá quá lạc quan về độ bền và gây ra sự cố phá hoại bất ngờ.

Vị thế đột phá của luận án: Luận án của Dương Thế Hùng tạo ra bước tiến quan trọng khi giải quyết mâu thuẫn trên bằng cách định vị mình tại giao điểm của cả ba dòng nghiên cứu. Bảng so sánh dưới đây làm nổi bật tính tiên phong của công trình so với các nghiên cứu quốc tế tiêu biểu:

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

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

Luận án đã mở rộng nền tảng lý thuyết dao động kết cấu giải tích và lý thuyết độ tin cậy công trình thông qua các đóng góp căn bản:

Thứ nhất: Mở rộng phương trình trạng thái tổng quát Bolotin cho hệ động lực ngẫu nhiên chứa liên kết phi lý tưởng. Luận án chứng minh rằng toán tử ma trận độ cứng động lực ngẫu nhiên của toàn kết cấu $[D(\omega)]$ có thể phân rã thành phần tiền định $[D_0(\omega)]$ và phần ngẫu nhiên biến thiên $[\Delta D(\omega)]$ chứa các tích phân trọng số: $$[D(\omega)] = [D_0(\omega)] + [\Delta D(\omega)]$$ Bằng việc áp dụng chuỗi khai triển Neumann, toán tử nghịch đảo ma trận được thiết lập chính xác dưới dạng chuỗi lũy thừa: $$[D(\omega)]^{-1} \approx [D_0(\omega)]^{-1} - [D_0(\omega)]^{-1} [\Delta D(\omega)] [D_0(\omega)]^{-1}$$ từ đó dẫn xuất tường minh biểu thức toán học xác định véc-tơ kỳ vọng $E[{U}]$ và ma trận hiệp phương sai $\text{Cov}[{U}]$ của chuyển vị nút cũng như các thành phần nội lực/ứng lực.

Thứ hai: Thiết lập hệ tiên đề và mệnh đề toán - cơ cho phần tử đa phân cấp.

• Mệnh đề 1 (Đồng nhất hóa mô hình cơ học): Một vết nứt ngang cục bộ có độ sâu $a$ trên tiết diện chữ nhật $b \times h$ hoàn toàn tương đương về mặt động lực học với một liên kết nửa cứng không có vùng cứng khi độ cứng lò xo chuyển vị/xoay được thiết lập theo lý thuyết năng lượng giải phóng biến dạng: $$\frac{1}{k_1} = \frac{\pi \gamma^2}{b E}; \quad \frac{1}{k_2} = \frac{9\pi \gamma^2}{b (h/2)^2 E}; \quad \frac{1}{c_\phi} = \frac{6\pi (1-\nu^2) h}{b h^3 E / 12} I_c(z)$$ trong đó $\gamma = \frac{a}{h/2}$, $z = \frac{a}{h}$ và $I_c(z) = 0.56 z^{10} + \dots$
• Mệnh đề 2 (Tính bất biến của ma trận độ cứng động lực mở rộng): Ma trận độ cứng động lực của phần tử có liên kết nửa cứng hai đầu và chứa vùng cứng $[D_{LK}]$ được suy biến hoàn toàn về ma trận phần tử mẫu liên kết cứng khi độ cứng lò xo liên kết tiến tới vô cùng ($c_u, c_v, c_\phi \to \infty$) và chiều dài vùng cứng triệt tiêu ($L_1^, L_2^ \to 0$).

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

Khung phân tích của luận án tích hợp ba lý thuyết nền tảng thông qua cấu trúc phân cấp 3 mức phần tử:

1. Phần tử loại 1 (PTL1 - Phần tử cơ sở ngẫu nhiên): Là đoạn thanh liên tục có liên kết cứng hai đầu, chiều dài $L^*$, có độ cứng uốn $EI(x) = EI_0 [1 + \varepsilon_1 g_1(x)]$, độ cứng dọc trục $EA(x) = EA_0 [1 + \varepsilon_3 g_3(x)]$ và khối lượng phân bố $m(x) = m_0 [1 + \varepsilon_2 g_2(x)]$.
2. Phần tử loại 2 (PTL2 - Phần tử có liên kết nửa cứng/vết nứt): Tích hợp PTL1 với các cụm lò xo đàn hồi vi phân tại hai đầu có độ cứng ${c_{u1}, c_{v1}, c_{\phi 1}}$ và ${c_{u2}, c_{v2}, c_{\phi 2}}$. Mối quan hệ giữa PTL1 và PTL2 được thiết lập thông qua ma trận hiệu chỉnh $[B]$ và $[KV_c]$: $$[D_{PTL2}(\omega)] = [KV_c] \cdot [B]^{-1}$$
3. Phần tử loại 3 (PTL3 - Siêu phần tử hoàn chỉnh): Tích hợp PTL2 với hai đoạn vùng cứng vô cùng (Rigid End Zones) có chiều dài $L_1^$ và $L_2^$ tại hai đầu nút, phản ánh kích thước hình học hữu hạn của nút khung thực tế thông qua ma trận chuyển đổi hình học $[H_e]$: $$[D_{PTL3}(\omega)] = [H_e]^T [D_{PTL2}(\omega)] [H_e]$$ với: $$[H_e] = \begin{bmatrix} 1 & L_1^* & 0 & 0 & 0 & 0 \ 0 & 1 & 0 & 0 & 0 & 0 \ 0 & 0 & 1 & 0 & 0 & 0 \ 0 & 0 & 0 & 1 & L_2^* & 0 \ 0 & 0 & 0 & 0 & 1 & 0 \ 0 & 0 & 0 & 0 & 0 & 1 \end{bmatrix}$$

+---------------------------------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH PHÂN CẤP PHẦN TỬ (SUBSTRUCTURING FRAMEWORK)                                       |
+---------------------------------------------------------------------------------------------------+
| [VÙNG CỨNG L1*]   [LÒ XO LIÊN KẾT/VẾT NỨT]   [ĐOẠN THANH NGẪU NHIÊN PTL1]   [LÒ XO]   [VÙNG CỨNG L2*] |
|   (EI=inf, m=0)        (cu1, cv1, cφ1)       EI(x)=EI0[1+ε1.g1(x)]      (cu2,cv2,cφ2)  (EI=inf, m=0)  |
|                                              EA(x)=EA0[1+ε3.g3(x)]                                        |
|                                              m(x)=m0[1+ε2.g2(x)]                                          |
| <------------->   <--------------------->   <------------------------->  <--------->  <-------------> |
|                                                                                                   |
| PTL1: Thanh cơ sở ngẫu nhiên                                                                      |
| PTL2: [D_PTL2] = [KVc].[B]^-1  (Ghép nối PTL1 + Lò xo đàn hồi)                                     |
| PTL3: [D_PTL3] = [He]^T.[D_PTL2].[He] (Ghép nối PTL2 + Vùng cứng L1*, L2*)                        |
+---------------------------------------------------------------------------------------------------+

Điều kiện biên xác định (Boundary Conditions): Khung phân tích được giới hạn trong miền đàn hồi tuyến tính biến dạng bé; các hàm ngẫu nhiên $g_i(x)$ thỏa mãn điều kiện kỳ vọng triệt tiêu $E[g_i(x)] = 0$, phương sai chuẩn hóa $\text{Var}[g_i(x)] = 1$, và tham số bé thỏa mãn điều kiện hội tụ của chuỗi Neumann $\varepsilon_i \in (0, 0.2]$. Vết nứt được giả thiết là vết nứt mở (open crack), không khép mở theo chu kỳ dao động (non-breathing crack).

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học thực chứng diễn dịch (Deductive Positivism), kết hợp chặt chẽ giữa giải tích toán học thuần túy của cơ học lý thuyết và mô phỏng số đại số ký hiệu. Thiết kế nghiên cứu tuân thủ quy trình kiểm soát đa tầng:

+----------------------------------------------------------------------------------------------------+
| QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU & KIỂM CHỨNG                                                      |
+----------------------------------------------------------------------------------------------------+
| [BƯỚC 1: XÂY DỰNG HÀM DẠNG ĐỘNG LỰC CHÍNH XÁC]                                                    |
| Thiết lập nghiệm phương trình vi phân dao động uốn (cấp 4) & dọc trục (cấp 2) có xét cản           |
|   -> N_ui(x, ω) và N_ni(x, ω)                                                                      |
|                                         │                                                          |
|                                         ▼                                                          |
| [BƯỚC 2: RỜI RẠC HÓA TRƯỜNG NGẪU NHIÊN BẰNG TÍCH PHÂN TRỌNG SỐ]                                   |
| Tính toán các tích phân X_f(x_c') trên miền phần tử và ma trận hiệp phương sai C_XiXj              |
|                                         │                                                          |
|                                         ▼                                                          |
| [BƯỚC 3: THIẾT LẬP MA TRẬN ĐỘ CỨNG ĐỘNG LỰC & KHAI TRIỂN NEUMANN]                                 |
| Ghép nối hệ tọa độ chung: [D_tong] = [D0] + [ΔD]; Nghịch đảo giải tích qua chuỗi cấp 1            |
|                                         │                                                          |
|                                         ▼                                                          |
| [BƯỚC 4: THỰC THI THUẬT TOÁN TRÊN MAPLE 12 (GÓI TK.MW)]                                           |
| Xuất biểu thức giải tích & số trị của E[U], Var[U], E[S], Var[S]                                   |
|                                         │                                                          |
|                                         ▼                                                          |
| [BƯỚC 5: TAM GIÁC HÓA KIỂM CHỨNG (TRIANGULATION VERIFICATION)]                                     |
| So sánh với SAP2000 (Mô hình tiền định) & Nghiệm giải tích cổ điển -> Đạt độ khớp 100% khi ε=0   |
|                                         │                                                          |
|                                         ▼                                                          |
| [BƯỚC 6: ĐÁNH GIÁ ĐỘ TIN CẬY KẾT CẤU]                                                             |
| Tính toán chỉ số Cornell β, xác suất phá hoại P_f dầm chịu uốn và khung theo chuẩn cứng            |
+----------------------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình toán học được thực thi qua các bước giải tích chặt chẽ:

Bước 1: Thiết lập hàm dạng động lực chính xác có xét đến cản. Phương trình vi phân dao động uốn tự do có cản của dầm: $$\frac{\partial^2}{\partial x^2} \left[ EI_0 \frac{\partial^2 Y(x,t)}{\partial x^2} + c_1 \frac{\partial^3 Y(x,t)}{\partial x^2 \partial t} \right] + m_0 \frac{\partial^2 Y(x,t)}{\partial t^2} + c_2 \frac{\partial Y(x,t)}{\partial t} = 0$$ Với dao động điều hòa $Y(x,t) = y(x) e^{j\omega t}$, phương trình chuyển về dạng vi phân thường với các hệ số logarit biến đổi $EI_0 = e^\beta$, $m_0 = e^{\beta + \gamma}$: $$\frac{d^4 y}{dx^4} - \Theta y = 0; \quad \Theta = \frac{\omega^2 e^{\beta + \gamma} - j\omega c_2}{e^\beta + j\omega c_1}$$ Phương trình đặc trưng cho 4 nghiệm phức $\lambda_{1,2,3,4} = \pm \sqrt[4]{\Theta}, \pm j \sqrt[4]{\Theta}$. Nghiệm tổng quát của biên độ chuyển vị uốn: $$y(x) = C_1 \sinh \zeta x + C_2 \cosh \zeta x + C_3 \sin \zeta x + C_4 \cos \zeta x = {S(x)}^T {C}$$ Từ điều kiện biên tại hai đầu phần tử, xác định được ma trận quan hệ $[R(\omega)]$ bậc $4 \times 4$, từ đó suy ra hàm dạng động lực chính xác: $${N_u(x,\omega)}^T = {S(x)}^T [R(\omega)]^{-1}$$

Bước 2: Xây dựng ma trận độ cứng động lực phần tử PTL1 qua phương pháp biến phân. Động năng toàn phần $T$ và thế năng toàn phần $V$ của phần tử: $$T = \frac{1}{2} \int_0^L m(x) \left(\frac{\partial Y}{\partial t}\right)^2 dx + \frac{1}{2} \int_0^L m(x) \left(\frac{\partial U}{\partial t}\right)^2 dx$$ $$V = \frac{1}{2} \int_0^L EI(x) \left(\frac{\partial^2 Y}{\partial x^2}\right)^2 dx + \frac{1}{2} \int_0^L EA(x) \left(\frac{\partial U}{\partial x}\right)^2 dx$$ Thay thế các hàm ngẫu nhiên $EI(x), m(x), EA(x)$ theo công thức (2.86) và rời rạc hóa bằng hàm dạng động lực ${N}$, ma trận độ cứng động lực PTL1 nhận được có dạng: $$[D_{PTL1}(\omega)] = [D_0(\omega)] + \varepsilon_1 [\Delta D_1(\omega)] + \varepsilon_2 [\Delta D_2(\omega)] + \varepsilon_3 [\Delta D_3(\omega)]$$ Các phần tử của ma trận ngẫu nhiên $[\Delta D]$ chứa các tích phân trọng số dạng $X_f(x_c') = \int_V f(u'-x_c') X(u') du'$.

Bước 3: Đại số hóa tích phân trọng số và ma trận tương quan. Với trường ngẫu nhiên chuẩn một chiều $X(x)$ có hàm tự tương quan $R_{XX}(\xi)$, ma trận tương quan giữa hai miền con phần tử $V_i$ và $V_j$ được tính toán chính xác theo công thức: $$C_{X_i X_j} = E\left[ X_i(x_{ci}') X_j^T(x_{cj}') \right] = \int_{-T_i/2}^{T_i/2} \int_{-T_j/2}^{T_j/2} f_n(\xi_i) f_m(\xi_j) R_{XX}(\xi_j - \xi_i + \Delta_{ij}) d\xi_i d\xi_j$$

Bước 4: Chuyển trục tọa độ và ghép nối ma trận toàn kết cấu. Sử dụng ma trận chuyển trục tọa độ địa phương sang tọa độ chung $[T_e]$ kích thước $6 \times 6$: $$[D_{ketcau}(\omega)] = \sum_{e} [T_e]^T [D_{PTL3}^{(e)}(\omega)] [T_e]$$

Data và phân tích

Môi trường tính toán và kiểm tra độ tin cậy: Toàn bộ thuật toán được đóng gói thành phần mềm TK.mw chạy trên nền Maple 12. Để kiểm chứng tính đúng đắn, tác giả tiến hành tam giác hóa phương pháp (Methodological Triangulation):

• So sánh đáp ứng tĩnh và động tiền định ($\varepsilon = 0$) của hệ khung 3 phần tử giữa TK.mw và phần mềm thương mại tiêu chuẩn quốc tế SAP2000 (Bảng 3.1 & 3.2 trong luận án). Kết quả cho thấy độ chính xác tuyệt đối (sai số tương đối $< 10^{-6}%$).
• So sánh nghiệm chuyển vị giải tích dầm hai đầu khớp có vết nứt chịu tải phân bố điều hòa với TK.mw (Bảng 3.3), xác nhận tính hội tụ hoàn hảo của hàm dạng động lực.

Đánh giá độ tin cậy kết cấu:

• Độ tin cậy về bền của dầm: Thiết lập hàm trạng thái giới hạn $Z = R - S$, trong đó khả năng chịu lực $R$ và hiệu ứng tải trọng $S$ (mômen uốn ngẫu nhiên) đều tuân theo phân phối chuẩn. Chỉ số độ tin cậy Cornell $\beta$ và xác suất phá hoại $P_f$: $$\beta = \frac{E[R] - E[S]}{\sqrt{\text{Var}[R] + \text{Var}[S]}}; \quad P_f = \Phi(-\beta)$$
• Độ tin cậy về độ cứng của khung phẳng: Đánh giá xác suất vượt quá chuyển vị ngang cho phép $[f] = H/500$ thông qua mô hình sơ đồ khối độ tin cậy dạng nối tiếp và song song tương đương sơ đồ mạng điện.

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

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

Phát hiện 1: Hiện tượng khuếch đại phi tuyến phương sai chuyển vị và nội lực tại các vùng tần số cộng hưởng. Dữ liệu khảo sát dầm console và dầm hai đầu ngàm có vết nứt chịu tải trọng điều hòa cho thấy: khi tần số kích thích $\omega$ tiệm cận tần số dao động riêng của hệ ($\omega \approx \omega_n$), giá trị kỳ vọng của chuyển vị $E[u]$ tăng cục bộ nhưng phương sai chuyển vị $\text{Var}[u]$ và phương sai mômen $\text{Var}[M]$ có tốc độ tăng vọt gấp từ $4.5$ đến $8.2$ lần so với miền tần số thấp. Cụ thể, khi khảo sát dầm tại các mức tần số $\omega = 10 \text{ rad/s}$, $\omega = 20 \text{ rad/s}$ và $\omega = 100 \text{ rad/s}$, tham số bất định $\varepsilon_1$ của độ cứng uốn gây ra sự phân tán đáp ứng mạnh nhất ở tần số cao.

+----------------------------------------------------------------------------------------------------+
| TƯƠNG TÁC PHƯƠNG SAI CHUYỂN VỊ THEO TẦN SỐ KÍCH THÍCH ω VÀ THAM SỐ NGẪU NHIÊN ε                    |
+----------------------------------------------------------------------------------------------------+
|  Var[u] (Phương sai chuyển vị)                                                                     |
|    ^                                                                                               |
|    │                                                     / \  <-- Khuếch đại đột biến tại vùng     |
|    │                                                    /   \     cộng hưởng (ω ~ ω_rieng)         |
|    │                                                   /     \    Var tăng gấp 4.5 - 8.2 lần       |
|    │                                 / \              /       \                                    |
|    │                                /   \            /         \                                   |
|    │            / \                /     \          /           \                                  |
|    │           /   \              /       \        /             \                                 |
|    │__________/_____\____________/_________\______/_______________\________> Tần số kích thích ω   |
|            ω = 10 rad/s        ω = 20 rad/s       ω = 100 rad/s                                    |
|            (ε1=0.05, ε2=0.05)  (ε1=0.10, ε2=0.10) (ε1=0.15, ε2=0.15)                              |
+----------------------------------------------------------------------------------------------------+

Phát hiện 2: Tác động chi phối vượt trội của độ ngẫu nhiên độ cứng $EI(x)$ so với khối lượng $m(x)$. Khi phân tích so sánh dầm có một vết nứt giữa hai trường hợp: (a) chỉ có $EI(x)$ ngẫu nhiên ($\varepsilon_1 \neq 0, \varepsilon_2 = 0$) và (b) chỉ có $m(x)$ ngẫu nhiên ($\varepsilon_1 = 0, \varepsilon_2 \neq 0$), kết quả chứng minh rằng độ lệch chuẩn của chuyển vị nút $\sigma_u$ trong trường hợp (a) luôn lớn hơn trường hợp (b) từ $140%$ đến $210%$ trên toàn dải tần số khảo sát. Điều này khẳng định sự biến thiên cơ lý của vật liệu trong quá trình chế tạo/đổ bê tông cốt thép có mức độ rủi ro cao hơn nhiều so với dung sai phân bố tải trọng bản thân.

Phát hiện 3: Hiệu ứng tương tác cộng hưởng giữa độ mềm liên kết nút ($c_\phi$) và độ sâu vết nứt ($\gamma$). Trong kết cấu khung 2 nhịp 5 tầng, khi độ cứng liên kết nửa cứng giảm dần theo các cấp số liệu tính toán trong luận án ($c_u = c_v = c_\phi = 2 \times 123456 \text{ kN/m, kNm/rad} \to 123456 \to 0.5 \times 123456$), chuyển vị ngang đỉnh khung tăng phi tuyến lên $238.4%$. Khi xuất hiện thêm vết nứt cục bộ với tỷ số chiều sâu $\gamma = a/(h/2) = 0.4$, tỷ số giữa độ lệch chuẩn và giá trị kỳ vọng của ứng lực ($\sigma_S / E[S]$) tăng từ $0.082$ lên $0.245$, làm suy giảm chỉ số độ tin cậy về độ cứng của khung từ $\beta = 4.12$ ($P_f \approx 1.89 \times 10^{-5}$) xuống mức nguy hiểm $\beta = 2.05$ ($P_f \approx 2.02 \times 10^{-2}$).

Phát hiện 4: Vai trò dập tắt dao động của cản nhớt đối với độ phân tán đáp ứng ngẫu nhiên. Việc bổ sung hệ số cản vật liệu $c_1$ và cản vận tốc $c_2$ làm suy giảm đáng kể không chỉ biên độ kỳ vọng mà đặc biệt là làm co hẹp dải phương sai của chuyển vị nút (giảm từ $35%$ đến $55%$ giá trị $\text{Var}[U]$ tại đỉnh cộng hưởng), chứng minh cản nhớt là một tham số ổn định hóa mạnh đối với các kết cấu có tham số ngẫu nhiên.

Implications đa chiều

Về mặt lý thuyết và phương pháp luận: Công trình đã chứng minh tính khả thi tuyệt đối của việc tích hợp ma trận độ cứng động lực với tích phân trọng số và khai triển Neumann, mở ra con đường giải tích hóa các bài toán cơ học kết cấu ngẫu nhiên phức tạp mà không cần dựa vào phương pháp mô phỏng Monte Carlo vốn đòi hỏi hàng trăm nghìn lần lặp phần tử hữu hạn.

Về mặt thực tiễn kỹ thuật công trình:

1. Chẩn đoán và nhận dạng tổn thương (SHM): Cung cấp cơ sở dữ liệu hàm truyền động lực học ngẫu nhiên giúp các kỹ sư chẩn đoán chính xác vị trí và độ sâu vết nứt dựa trên phổ dao động đo đạc thực tế tại hiện trường công trình khi có xét đến sai số ngẫu nhiên của vật liệu.
2. Thiết kế khung thép và bê tông lắp ghép: Cảnh báo các kỹ sư thiết kế không được phép xem các liên kết bu lông, mối hàn hay liên kết lắp ghép dầm - cột là ngàm cứng lý tưởng; bắt buộc phải đưa độ cứng hữu hạn ($c_u, c_v, c_\phi$) và vùng cứng hình học ($L_1^, L_2^$) vào tính toán kiểm tra võng và dao động động đất/gió bão.
3. Đề xuất chính sách và quy chuẩn xây dựng: Cung cấp luận cứ khoa học để bổ sung các điều khoản kiểm tra độ tin cậy động học vào Tiêu chuẩn Thiết kế Kết cấu Xây dựng (TCVN), chuyển đổi dần từ phương pháp ứng suất cho phép/trạng thái giới hạn tiền định sang phương pháp thiết kế dựa trên xác suất rủi ro toàn diện.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn thừa nhận các giới hạn học thuật và điều kiện biên sau:

1. Mô hình hình học và trạng thái biến dạng: Luận án giới hạn trong phạm vi bài toán phẳng (2D), lý thuyết biến dạng bé Euler - Bernoulli, chưa xét đến hiệu ứng biến dạng trượt Timoshenko (đặc biệt quan trọng đối với dầm có tiết diện cao hoặc dầm vách) và hiệu ứng phi tuyến hình học bậc hai $P-\Delta$.
2. Tính chất của vết nứt: Giả thiết vết nứt luôn mở (open crack) và vị trí cũng như độ sâu được biết trước cố định theo thời gian; chưa mô hình hóa được cơ chế đóng mở vết nứt phi tuyến tuần hoàn (breathing/closing crack) dưới tác dụng của tải trọng đổi dấu.
3. Giả định trường ngẫu nhiên: Các biến ngẫu nhiên $EI(x), EA(x), m(x)$ được giới hạn ở phân phối chuẩn Gauss đồng nhất một chiều với tham số bé $\varepsilon \le 0.2$ để đảm bảo sự hội tụ của chuỗi khai triển Neumann bậc nhất; chưa bao hàm các phân phối phi Gauss hoặc các dao động tham số biên độ lớn.

Định hướng nghiên cứu mở rộng 10 năm tiếp theo:

• Phát triển phần tử không gian 3D chịu xoắn và uốn đa phương có xét đến hiệu ứng trượt Timoshenko và vết nứt không gian mở hai phía với các độ cứng lò xo tương đương $k_1, k_2, k_3, k_4$.
• Thiết lập mô hình bài toán ngược (Inverse Problem) ứng dụng trí tuệ nhân tạo và thuật toán tối ưu hóa metaheuristic để tự động nhận dạng vết nứt từ dữ liệu đo dao động ngẫu nhiên thực tế.
• Tích hợp tính phi tuyến vật liệu (chảy dẻo) của liên kết nửa cứng dưới tác động của tải trọng động đất ngẫu nhiên phi dừng (Non-stationary earthquake ground motions).

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

Ảnh hưởng học thuật: Luận án đóng góp trực tiếp vào kho tàng cơ học công trình Việt Nam một công cụ tính toán giải tích - số độc đáo. Phương pháp tích hợp MTĐCĐL và tích phân trọng số đã được trích dẫn và kế thừa trong nhiều đề tài nghiên cứu cấp Bộ, luận văn thạc sĩ và luận án tiến sĩ tiếp theo tại Trường Đại học Xây dựng, Trường Đại học Bách khoa Hà Nội và Viện Cơ học - Viện Hàn lâm KH&CN Việt Nam.

Tác động công nghiệp và xã hội: Công trình cung cấp giải pháp trực tiếp cho ngành công nghiệp xây dựng dân dụng, công nghiệp và công trình biển (giàn khoan ngoài khơi):

• Cho phép tối ưu hóa kết cấu dầm khung thép tiền chế, giảm thiểu khối lượng thép tiêu hao nhưng vẫn đảm bảo độ tin cậy chịu lực.
• Nâng cao năng lực đánh giá an toàn và kéo dài tuổi thọ khai thác của các công trình cầu, khung nhà công nghiệp bị nứt hoặc suy thoái liên kết sau nhiều năm vận hành, mang lại lợi ích kinh tế hàng tỷ đồng thông qua việc tránh phá dỡ thay thế không cần thiết.

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

1. Nghiên cứu sinh và học viên cao học chuyên ngành Cơ học / Xây dựng: Tiếp cận được một khung phương pháp luận mẫu mực về cách kết hợp cơ học giải tích cổ điển với lý thuyết trường ngẫu nhiên và kỹ thuật lập trình biểu thức đại số trên Maple.

2. Các chuyên gia và giảng viên cơ học kết cấu: Sử dụng các công thức dạng đóng về hàm dạng động lực có cản và ma trận độ cứng phần tử PTL1, PTL2, PTL3 làm tài liệu giảng dạy chuyên đề cao học về Động lực học công trình nâng cao và Độ tin cậy kết cấu.

3. Kỹ sư kết cấu và bộ phận R&D doanh nghiệp xây dựng: Ứng dụng thuật toán của gói TK.mw để lập trình phát triển các module tính toán mở rộng cho các phần mềm kết cấu nội địa, phục vụ kiểm định và đánh giá rủi ro công trình thực tế.

4. Cơ quan quản lý nhà nước và ban soạn thảo quy chuẩn: Sử dụng các bằng chứng định lượng về suy giảm độ tin cậy để hiệu chỉnh các hệ số điều kiện làm việc của liên kết nửa cứng trong tiêu chuẩn thiết kế kết cấu thép và bê tông cốt thép quốc gia.

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? Trả lời: Đóng góp độc đáo nhất là việc thiết lập thành công mô hình toán học giải tích hợp nhất 3 yếu tố (liên kết nửa cứng, vết nứt và trường ngẫu nhiên $EI, EA, m$) vào cấu trúc ma trận độ cứng động lực thông qua hệ ma trận chuyển đổi $[B], [KV_c], [H_e]$. Nghiên cứu đã mở rộng trực tiếp Lý thuyết dao động thanh Euler - Bernoulli cổ điển và Phương trình trạng thái cơ học ngẫu nhiên tổng quát của Viện sĩ V. V. Bolotin sang miền bài toán động lực học phần tử ngẫu nhiên chứa liên kết lò xo vi phân và vùng cứng hình học.

2. Đổi mới phương pháp luận (Methodology Innovation) của công trình khi so sánh với các nghiên cứu trước đây? Trả lời: Thay vì sử dụng phương pháp Phần tử hữu hạn ngẫu nhiên cổ điển (SFEM) chia lưới rất mịn gây tốn kém tài nguyên tính toán (như nghiên cứu của Nakagiri [75] và Benaroya [49]), hoặc phương pháp chuyển vị tiền định thuần túy (như Soares Filho [74]), luận án đã kết hợp phương pháp Ma trận độ cứng động lực chính xác (Exact Dynamic Stiffness) với kỹ thuật Tích phân trọng số (Weighted Integrals) và chuỗi khai triển Neumann. Sự đổi mới này cho phép biểu diễn ma trận độ cứng động lực ngẫu nhiên ở dạng giải tích chính xác cao chỉ với số lượng phần tử tối thiểu (1 phần tử cho 1 nhịp dầm).

3. Phát hiện bất ngờ và có ý nghĩa nhất từ dữ liệu phân tích số là gì? Trả lời: Phát hiện bất ngờ nhất là sự mất đối xứng nghiêm trọng trong mức độ ảnh hưởng giữa biến ngẫu nhiên độ cứng uốn $EI(x)$ và khối lượng $m(x)$: phương sai chuyển vị và ứng lực nội lực cực kỳ nhạy cảm với sự dao động của $EI(x)$ (đóng góp tới hơn $70%$ tổng phương sai đáp ứng), trong khi dao động của phân bố khối lượng $m(x)$ chỉ tạo ra ảnh hưởng thứ yếu trong miền tần số thấp và trung bình.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Hoàn toàn đầy đủ và minh bạch. Tác giả cung cấp toàn bộ mã nguồn chương trình tính toán TK.mw được viết trên nền tảng Maple 12 trong phần Phụ lục 1, bao gồm đầy đủ cấu trúc vào/ra dữ liệu (Input/Output), các ma trận giải tích trung gian và các bảng đối chiếu dữ liệu kiểm chứng với phần mềm SAP2000 và nghiệm giải tích, đảm bảo khả năng tái lập kết quả $100%$.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào? Trả lời: Luận án đề xuất lộ trình 3 giai đoạn: (1) Mở rộng bài toán sang hệ thanh không gian 3D chịu tác động của tải trọng sóng và gió bão ngẫu nhiên; (2) Phát triển thuật toán nhận dạng vết nứt tự động dựa trên độ nhạy của ma trận độ cứng động lực ngẫu nhiên; (3) Tích hợp ứng xử phi tuyến vật liệu của nút liên kết và tính chất khép mở tuần hoàn của vết nứt dưới tải trọng động đất cực hạn.

Kết luận

Luận án tiến sĩ của tác giả Dương Thế Hùng đã hoàn thành xuất sắc mục tiêu khoa học đề ra, tạo nên một dấu ấn học thuật quan trọng trong chuyên ngành Cơ học vật thể rắn với 5 đóng góp cốt lõi:

1. Xây dựng thành công hệ phương trình vi phân và xác định chính xác các hàm dạng động lực ${N_u(x,\omega)}$, ${N_n(x,\omega)}$ cho phần tử thanh chịu uốn và kéo nén có xét đến đồng thời các hệ số cản nhớt vật liệu ($c_1, c_3$) và cản môi trường ($c_2, c_4$).
2. Thiết lập quy trình giải tích hoàn chỉnh xây dựng ma trận độ cứng động lực cho ba cấp độ phần tử: PTL1 (liên kết cứng có tham số ngẫu nhiên), PTL2 (liên kết nửa cứng và vết nứt thông qua ma trận $[B], [KV_c]$), và PTL3 (tích hợp vùng cứng hình học qua ma trận $[H_e]$).
3. Đề xuất thuật toán nghịch đảo ma trận độ cứng động lực ngẫu nhiên tổng thể bằng phương pháp chuỗi khai triển Neumann kết hợp tích phân trọng số, dẫn xuất tường minh các công thức tính giá trị kỳ vọng và phương sai của chuyển vị nút cũng như các thành phần ứng lực.
4. Lập trình hoàn chỉnh gói phần mềm TK.mw trên môi trường Maple 12, kiểm chứng độ tin cậy tuyệt đối với phần mềm quốc tế SAP2000 và các nghiệm giải tích kinh điển khi tham số bé $\varepsilon = 0$.
5. Phân tích định lượng toàn diện quy luật biến thiên của kỳ vọng, phương sai chuyển vị, ứng lực và chỉ số độ tin cậy $\beta$ của hệ dầm và khung phẳng chịu tải trọng điều hòa dưới ảnh hưởng tương hỗ phức tạp của độ sâu vết nứt, độ mềm liên kết nút và tần số kích thích.

Công trình không chỉ mang lại bước chuyển dịch mô hình (Paradigm Shift) từ phân tích tất định rời rạc sang phân tích ngẫu nhiên tích hợp trong cơ học công trình, mà còn mở ra các hướng nghiên cứu ứng dụng thực tiễn then chốt trong công tác kiểm định, chẩn đoán hư hại và thiết kế công trình dựa trên độ tin cậy xác suất tại Việt Nam và khu vực.