Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Một phương pháp đánh giá mức độ an toàn của kết cấu khung chịu tải trọng động theo lý thuyết tập mờ" do nghiên cứu sinh Lê Công Duy thực hiện dưới sự hướng dẫn khoa học của GS. Lê Xuân Huỳnh, bảo vệ năm 2014 tại Trường Đại học Xây dựng Hà Nội (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.08).

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

Trong tính toán thiết kế kết cấu công trình, các đại lượng đầu vào liên quan đến đặc trưng vật liệu, hình học và tải trọng thường chứa đựng các thông tin ngẫu nhiên, không rõ ràng hoặc không thể chính xác hóa—gọi chung là các đại lượng không chắc chắn. Các phương pháp tính toán truyền thống dựa trên ứng suất cho phép hoặc trạng thái giới hạn (mô hình tiền định) chưa phản ánh toàn diện sự làm việc thực tế của kết cấu do chưa xét đầy đủ ảnh hưởng của các yếu tố không chắc chắn. Hệ quả là trong thực tế, nhiều công trình vẫn xảy ra hư hỏng dù được thiết kế với hệ số an toàn tương đối lớn.

Để mô tả các đại lượng không chắc chắn, các nhà nghiên cứu thường sử dụng số khoảng, đại lượng ngẫu nhiên, số mờ hoặc đại lượng ngẫu nhiên - mờ:

  • Mô hình ngẫu nhiên: Đánh giá kết cấu dựa trên lý thuyết độ tin cậy và lý thuyết xác suất thống kê. Mô hình này đòi hỏi tập dữ liệu thống kê đủ lớn để xác định hàm mật độ phân bố xác suất; tuy nhiên, trong thực tế xây dựng, nhiều trường hợp số liệu khảo sát bị thiếu hụt hoặc không tuân theo các quy luật phân phối chuẩn tắc.
  • Mô hình mờ: Sử dụng lý thuyết tập mờ (Fuzzy Set Theory) để xử lý các thông tin nhậm nhằng (vague), mờ (fuzzy) khi thiếu dữ liệu định lượng. Tại thời điểm nghiên cứu, việc ứng dụng lý thuyết tập mờ trong phân tích trạng thái kết cấu và đánh giá mức độ an toàn của khung nhà chịu tải trọng động vẫn là hướng tiếp cận mới, số lượng công trình nghiên cứu hoàn chỉnh còn hạn chế.

Khoảng trống nghiên cứu mà luận án tập trung giải quyết là: Xây dựng thuật toán giải hệ phương trình phương pháp phần tử hữu hạn (PTHH) có tham số mờ và thiết lập công thức đánh giá độ tin cậy mờ (mức độ an toàn) cho kết cấu khung nhiều tầng chịu tải trọng động khi các tham số đầu vào (độ cản, vật liệu, biên độ và tần số tải trọng) được mô tả dưới dạng số mờ tam giác.

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

Mục tiêu tổng quát của luận án bao gồm:

  1. Nghiên cứu ứng dụng lý thuyết tập mờ, đề xuất cách giải phương trình cơ bản của phương pháp PTHH có tham số đầu vào mờ.
  2. Vận dụng, triển khai và chứng minh công thức đánh giá độ tin cậy mờ của kết cấu (Phương pháp tỷ số diện tích).
  3. Áp dụng tính toán và đánh giá mức độ an toàn đối với khung phẳng nhiều tầng chịu tải trọng động trong trường hợp các yếu tố tác động (độ cản kết cấu, đặc trưng vật liệu, đặc trưng biên độ và tần số của tải trọng động) được xét dưới dạng các số mờ tam giác.

Nội dung nghiên cứu cụ thể được triển khai theo 5 nhiệm vụ:

  • Nhiệm vụ 1: Vận dụng cơ sở lý thuyết tập mờ, các thuật toán hỗ trợ tính toán số mờ kết hợp phần mềm Maple.13.
  • Nhiệm vụ 2: Phân tích các phương pháp đánh giá mức độ an toàn cho kết cấu theo quan điểm ngẫu nhiên của lý thuyết xác suất và quan điểm mờ của lý thuyết tập mờ; từ đó triển khai và chứng minh công thức đánh giá mức độ an toàn theo lý thuyết tập mờ.
  • Nhiệm vụ 3: Nghiên cứu phương pháp PTHH có tham số mờ, đề xuất cách giải phương trình đại số tuyến tính có tham số mờ; ứng dụng phân tích kết cấu chịu tải trọng tĩnh với tham số đầu vào dạng số mờ tam giác.
  • Nhiệm vụ 4: Nghiên cứu mô hình tính kết cấu khung chịu tải trọng động, thiết lập phương trình dao động cho kết cấu khung phẳng nhiều tầng chịu tải trọng động có tham số đầu vào mờ.
  • Nhiệm vụ 5: Ứng dụng công thức đã triển khai để đánh giá mức độ an toàn cho kết cấu khung khi các yếu tố đầu vào (vật liệu, khối lượng, độ cản, tải trọng động) là các số mờ tam giác; tính toán và so sánh kết quả với một số phương pháp đánh giá khác.

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

  • Đối tượng nghiên cứu: Mức độ an toàn của kết cấu khung phẳng nhiều tầng chịu tải trọng động có tham số đầu vào không chắc chắn dạng mờ.
  • Phạm vi và giả thiết học thuật:
    • Kết cấu khung phẳng làm việc trong giai đoạn đàn hồi tuyến tính.
    • Chân cột liên kết ngàm cứng trên nền móng tuyệt đối (chưa xét đến liên kết đàn hồi tại chân công trình và tương tác động lực học giữa đất nền và công trình).
    • Sàn và dầm xem như tuyệt đối cứng trong mặt phẳng của nó; bỏ qua biến dạng dọc trục của dầm và cột; mỗi tầng xem như một bậc tự do chuyển vị ngang của khối lượng tập trung đặt tại mức sàn.
    • Tải trọng động tác dụng là hàm theo thời gian $t$, xét hai dạng chính: tải trọng điều hòa có chu kỳ hình sin và tải trọng xung hình chữ nhật.
    • Các biến đại lượng mờ đầu vào được giả thiết không tương quan.
    • Bài toán áp dụng ở mức độ cơ học kết cấu có số bậc tự do không lớn, tồn tại nghiệm giải tích đóng đối với các biến đầu vào mờ.

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

Luận án hệ thống hóa tiến trình phát triển của lý thuyết an toàn và độ tin cậy kết cấu trên thế giới và tại Việt Nam qua ba mô hình tính toán: tiền định, ngẫu nhiên và mờ.

Tiến trình mô hình hóa an toàn kết cấu:
[Mô hình Tiền định] (Ứng suất cho phép / Trạng thái giới hạn)
[Mô hình Ngẫu nhiên] (Xác suất thống kê / Lý thuyết độ tin cậy)
[Mô hình Mờ / Ngẫu nhiên - Mờ] (Lý thuyết tập mờ / Số học mờ)

Các nghiên cứu quốc tế

  • Giai đoạn khởi đầu của lý thuyết độ tin cậy: Ra đời từ thập niên 1930 trong kỹ thuật điện tử, máy tính, hàng không vũ trụ. Các công trình đầu tiên của Mayer và Khoialép thể hiện độ tin cậy qua ứng suất cho phép và hệ số an toàn. Khotsialov bắt đầu đưa thống kê toán học vào cơ học kết cấu, đặt nền móng cho phương pháp trạng thái giới hạn. Đầu thập niên 1950 tại Mỹ, Freudenthal và Pugsley phát triển bài toán độ tin cậy kết cấu chịu tải trọng tĩnh và mô hình thống kê tuổi thọ kết cấu chịu tải trọng động.
  • Chuẩn hóa thiết kế theo độ tin cậy: Tiêu chuẩn quốc tế ISO 2394:1998; tiêu chuẩn Trung Quốc GB 50153-92 (JB50153-92), JGJ125-99. Ủy ban Liên hiệp về An toàn Kết cấu (JCSS) phân chia tính toán độ tin cậy thành 3 mức: Mức 1 (phương pháp bán xác suất qua hệ số thành phần), Mức 2 (xấp xỉ lặp số, xấp xỉ tuyến tính hóa Taylor), Mức 3 (tính toán giải tích chính xác hàm mật độ phá hoại toàn phần).
  • Sự hình thành và phát triển của lý thuyết tập mờ: Lotfi Zadeh (1965, Đại học California) đề xuất lý thuyết tập mờ và logic mờ. Các ứng dụng kỹ thuật tiêu biểu: Mamdani (1977, điều khiển học), Fordon & Bezdek (1979, chẩn đoán y học), Cao & Chen (1983, dự báo khí tượng).
  • Ứng dụng tập mờ trong độ tin cậy kết cấu xây dựng:
    • Blockley (1979) [107]: Sử dụng tập mờ đánh giá an toàn công trình đang khai thác.
    • Brown (1979) [103]: Đề xuất phép đo độ tin cậy mờ cho công trình dựa trên tỷ số diện tích hàm thuộc.
    • Yao (1979) [89]: Tổng kết nội dung và ứng dụng lý thuyết tập mờ trong kỹ thuật công trình.
    • Shiraishi & Furuta (1989) [91]: Đánh giá độ tin cậy kết cấu dựa trên giao thoa giữa tải trọng ngẫu nhiên và khả năng chịu lực mờ.
    • Bernd Möller, Michael Beer và cộng sự (2004, 2008) [83], [84]: Nghiên cứu tính không chắc chắn ngẫu nhiên - mờ, phương pháp độ tin cậy bậc nhất mờ (FFORM) và thuật toán tối ưu mức-$\alpha$.
    • Các nghiên cứu khác: Mô hình hồi quy mờ tam giác [92]; phương pháp thực nghiệm trong cơ học kết cấu mờ [94]; mô hình giao thoa mờ tổng quát của Wood & Lodwick [96]; phương pháp nhãn mờ xây dựng hàm trạng thái giới hạn [90]; phân tích độ tin cậy hệ thống [104].

Các nghiên cứu tại Việt Nam

  • Về lý thuyết độ tin cậy ngẫu nhiên: Các nhà khoa học Phạm Khắc Hùng, Nguyễn Văn Phó, Phan Văn Khôi, Lê Xuân Huỳnh, Nguyễn Hữu Lộc, Lê Ngọc Thạch, Nguyễn Xuân Chính, Nguyễn Viết Trung... Luận án tiến sĩ của tác giả tại tài liệu [22] là công trình đầu tiên tại Việt Nam phân tích dao động và độ tin cậy của hệ thanh trực giao chịu tải trọng động ngẫu nhiên.
  • Về ứng dụng lý thuyết tập mờ trong kết cấu:
    • Các công bố của Nguyễn Văn Phó, Lê Xuân Huỳnh, Lê Ngọc Thạch, Bùi Đức Chính, Nguyễn Đình Xân [5], [23], [25], [27], [28], [44], [45], [46] về cơ học thông tin mờ, chẩn đoán vết nứt dầm bê tông cốt thép, khung nhà chịu gió bão, dầm liên kết nửa cứng.
    • Tiêu chuẩn xây dựng Việt Nam [65] ban hành chỉ dẫn đánh giá mức độ nguy hiểm của kết cấu nhà dựa trên nguyên lý tập mờ.
    • Luận án tiến sĩ của Bùi Đức Chính (2006) [78]: "Tính toán độ tin cậy của kết cấu theo lý thuyết tập mờ" là luận án tiến sĩ đầu tiên tại Việt Nam ứng dụng mô hình mờ đánh giá kết cấu chịu tải trọng tĩnh bằng phân tích hồi quy mờ.

Vị trí và điểm giải quyết của luận án

Luận án kế thừa mô hình đo độ tin cậy mờ của Brown [103] và các thuật toán số học mờ của Möller [84], khắc phục hạn chế của các phương pháp trước đó:

  1. Giải quyết tình trạng bất cập của mô hình lai ngẫu nhiên - mờ (vốn tích hợp hai hàm không cùng độ đo xác suất và mức độ thuộc).
  2. Mở rộng bài toán phân tích từ tải trọng tĩnh sang bài toán dao động động lực học của khung nhiều tầng có kể đến tham số mờ ở cả đầu vào và đầu ra.
  3. Cung cấp thuật giải thực hành giải hệ phương trình đại số tuyến tính mờ bằng công cụ tính toán ký hiệu trên phần mềm Maple.13.

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

Cơ sở lý thuyết tập mờ và số học mờ

  • Định nghĩa tập mờ: Tập mờ $\tilde{A}$ xác định trên tập nền kinh điển $X$ là tập hợp các cặp giá trị: $$\tilde{A} = \left{ (x, \mu_A(x)) \mid x \in X \right}$$ trong đó ánh xạ $\mu_A: X \to [0, 1]$ biểu diễn hàm thuộc (membership function), xác định mức độ thuộc của phần tử $x$ vào tập $\tilde{A}$.
  • Các đặc trưng hình học của tập mờ:
    • Độ cao: $h = \text{height}(\tilde{A}) = \sup_{x \in X} \mu_A(x)$. Nếu $h = 1$ thì tập mờ là chuẩn tắc (normal); nếu $h < 1$ là không chuẩn tắc.
    • Miền xác định: $\text{supp}(\tilde{A}) = {x \in X \mid \mu_A(x) > 0}$.
    • Miền tin cậy (Lõi): $\text{core}(\tilde{A}) = A_1 = {x \in X \mid \mu_A(x) = 1}$.
    • Lát cắt-$\alpha$ (Tập mức-$\alpha$): $A_\alpha = {x \in X \mid \mu_A(x) \ge \alpha}$ với $\alpha \in [0, 1]$.
  • Số mờ tam giác: Một số mờ $\tilde{A} = (a_1, a_2, a_3)$ được xác định bởi cận dưới $a_1$, giá trị trung tâm $a_2$ (có $\mu_A(a_2) = 1$) và cận trên $a_3$. Hàm thuộc có dạng: $$\mu_A(x) = \begin{cases} 0, & x < a_1 \ \frac{x - a_1}{a_2 - a_1}, & a_1 \le x \le a_2 \ \frac{a_3 - x}{a_3 - a_2}, & a_2 \le x \le a_3 \ 0, & x > a_3 \end{cases}$$ Với lát cắt mức-$\alpha$, khoảng biến thiên tương ứng là: $A_\alpha = [a_1 + \alpha(a_2 - a_1), a_3 - \alpha(a_3 - a_2)]$.
Đại lượng mờ Ký hiệu toán học Ý nghĩa trong bài toán kết cấu
Độ cứng nút khung $\tilde{k} = (k_1, k_2, k_3)$ Mô tả liên kết từ khớp ($k=0$) đến ngàm ($k=1$) hoặc liên kết nửa cứng
Mô đun đàn hồi $\tilde{E} = (E_1, E_2, E_3)$ Phản ánh sự phân tán tính chất cơ lý của bê tông, cốt thép
Độ cản nhớt $\tilde{c} = (c_1, c_2, c_3)$ Mô tả mức độ tiêu tán năng lượng dao động của hệ kết cấu
Tải trọng động $\tilde{F}(t) = \tilde{F}_0 \sin(\tilde{\omega} t)$ Phản ánh sự bất định của biên độ và tần số kích động cưỡng bức
  • Toán tử số học mờ:
    • Phân tích khoảng cổ điển: Xác định các phép tính cộng, trừ, nhân, chia trên các khoảng đóng $A = [a_1, a_2]$ và $B = [b_1, b_2]$. Phép cộng: $[a_1+b_1, a_2+b_2]$; Phép trừ: $[a_1-b_2, a_2-b_1]$; Phép nhân và chia theo quy tắc cực trị biên.
    • Nguyên lý mở rộng Zadeh và thuật toán min-max: Biến đổi hàm số qua các biến mờ $\mu_B(y) = \max_{\vec{x}: f(\vec{x})=y} \min(\mu_{A_1}(x_1), \dots, \mu_{A_n}(x_n))$. Nhược điểm là khối lượng rời rạc hóa trục hoành quá lớn khi số biến tăng.
    • Thuật toán tối ưu mức-$\alpha$ (Möller): Rời rạc hóa trục tung thành các mức $\alpha_k \in [0, 1]$. Tại mỗi mức $\alpha_k$, tìm cận dưới $y_{j, \alpha_k}^l$ và cận trên $y_{j, \alpha_k}^r$ thông qua việc giải hai bài toán tối ưu phi tuyến có ràng buộc: $$y_j = f_j(x_1, \dots, x_n) \to \text{Min / Max} \quad \text{với } (x_1, \dots, x_n) \in X_{\alpha_k}$$

Phương pháp phương trình vi phân dao động kết cấu

Theo nguyên lý D'Alembert, phương trình cân bằng động lực học tại tầng thứ $k$ có dạng: $$F_k^{qt} + F_k^{đh} + F_k^c = F_k(t)$$ Hệ phương trình vi phân dao động của hệ khung phẳng $n$ bậc tự do biểu diễn dưới dạng ma trận: $$[M]{\ddot{x}(t)} + [C]{\dot{x}(t)} + [K]{x(t)} = {F(t)}$$ trong đó $[M]$ là ma trận khối lượng tập trung đường chéo, $[C]$ là ma trận độ cản nhớt kết cấu, $[K]$ là ma trận độ cứng ngang của hệ khung, ${x(t)}, {\dot{x}(t)}, {\ddot{x}(t)}$ lần lượt là véc tơ chuyển vị, vận tốc và gia tốc ngang tương đối tại các mức sàn.

Công cụ và nguồn dữ liệu nghiên cứu

  • Công cụ tính toán số: Tác giả lập trình quy trình tính toán trên phần mềm đại số máy tính Maple.13, kết hợp tính toán ký hiệu (symbolic computation) và thuật toán quét lát cắt-$\alpha$ để giải hệ phương trình đại số tuyến tính mờ cũng như tìm đáp ứng dao động.
  • Nguồn dữ liệu: Dữ liệu thông số vật liệu, hình học, tải trọng động và độ cản được thu thập từ các tiêu chuẩn thiết kế hiện hành và các bài toán mẫu trong tài liệu cơ học kết cấu chuyên khảo.

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 phân tích cơ sở lý luận và tổng quan các mô hình đánh giá an toàn công trình, phân loại chi tiết các phương pháp hiện hành thành hai nhóm chính:

Các mô hình đánh giá an toàn kết cấu:
  1. Mô hình tiền định: Xác định mức độ an toàn qua tỷ số $n = R/Q$ hoặc hiệu số $M = R - Q$. Kết cấu xem là an toàn 100% khi $n > 1$ ($M > 0$) và mất an toàn khi $n < 1$ ($M < 0$). Mô hình này không lượng hóa được rủi ro hư hỏng do không xét đến sự biến động phân tán của các tham số thiết kế.
  2. Mô hình ngẫu nhiên: Xem hiệu ứng tải trọng $Q$ và khả năng chịu lực $R$ là các biến ngẫu nhiên tuân theo phân phối xác suất (thường là phân phối chuẩn). Khoảng an toàn $M = R - Q$ có giá trị trung bình $\bar{M} = \bar{R} - \bar{Q}$ và độ lệch chuẩn $\sigma_M = \sqrt{\sigma_R^2 + \sigma_Q^2}$. Xác suất an toàn $P_s = \text{Prob}(M > 0)$ và xác suất phá hoại $P_f = \text{Prob}(M \le 0) = 1 - P_s$. Chỉ số độ tin cậy định nghĩa là $\beta = \bar{M}/\sigma_M$. Khi hàm trạng thái giới hạn phi tuyến $M = f(X_1, \dots, X_n)$, phương pháp tuyến tính hóa Taylor bậc nhất quanh điểm kỳ vọng được sử dụng.
  3. Mô hình mờ: Khảo sát hai nhóm mô hình:
    • Nhóm giao thoa ngẫu nhiên - mờ ([78], [91], [93]): Tính xác suất hỏng mờ $\tilde{P}_f = \int \mu(x) f(x) dx$. Phương pháp này vướng nhược điểm không đồng nhất về mặt toán học do tích phân hai hàm khác nhau về độ đo (diện tích dưới đường mật độ xác suất bằng 1, trong khi tích phân hàm thuộc khác 1).
    • Nhóm giao thoa mờ - mờ ([25], [83], [96], [103], [107]): Đánh giá trực tiếp trên hai tập mờ $\tilde{Q}$ và $\tilde{R}$ hoặc tập hiệu mờ $\tilde{M} = \tilde{R} - \tilde{Q}$.
  4. Mô hình khung phẳng nhiều tầng: Thiết lập mô hình console nhiều bậc tự do tương đương, phương trình vi phân dao động matricial, phân tích tính không chắc chắn của độ cứng nút khung $\tilde{k} \in [0, 1]$, mô đun đàn hồi $\tilde{E}$, hệ số cản $\tilde{c}$, tải trọng gió động và động đất.

Chương 2: Một số phép toán của lý thuyết tập mờ và một cách giải thực hành hệ phương trình đại số tuyến tính mờ

Chương 2 trình bày có hệ thống các thuật toán số học mờ và đề xuất quy trình giải hệ phương trình đại số tuyến tính mờ $[\tilde{A}]{\tilde{x}} = {\tilde{b}}$ trong phương pháp phần tử hữu hạn:

Quy trình giải hệ phương trình PTHH mờ [Ã]{x̃} = {b̃}:
1. Nhập ma trận độ cứng mờ [Ã] và véc tơ tải trọng mờ {b̃} (Số mờ tam giác)
2. Áp dụng lát cắt-α: Chuyển các phần tử mờ về dạng khoảng đóng [a_ij^l(α), a_ij^r(α)]
3. Thiết lập biểu thức giải tích đóng của nghiệm x_i(α) trên phần mềm Maple.13
4. Thực hiện thuật toán tối ưu mức-α: Quét α từ 0 đến 1 để tìm [x_i^l(α), x_i^r(α)]
5. Tổng hợp hàm thuộc nghiệm đầu ra (Nội lực mờ, chuyển vị mờ dạng phi tuyến)
  • Quy tắc tính toán số học mờ: So sánh chi tiết giữa phương pháp phân tích khoảng, thuật toán min-max và thuật toán tối ưu mức-$\alpha$. Phân tích khoảng cho thấy tính chất phân phối không hoàn toàn ($A \times (B + C) \subseteq A \times B + A \times C$), do đó nếu thực hiện phép toán mờ quá sớm sẽ làm tăng độ phân tán và sai số của nghiệm đầu ra.
  • Thuật toán tối ưu mức-$\alpha$: Cho phép tìm chính xác đường bao cực tiểu và cực đại của nghiệm $x_j(\alpha)$ tại từng mức $\alpha_k$ mà không làm khuếch đại độ rộng mờ.
  • Ứng dụng tính toán tĩnh kết cấu: Luận án áp dụng giải bài toán thanh phẳng một chiều và khung phẳng chịu tải trọng tĩnh với các tham số đầu vào (kích thước hình học, mô đun đàn hồi $E$, tải trọng tập trung và phân bố) là số mờ tam giác. Kết quả đầu ra thu được là chuyển vị mờ tại các nút và nội lực mờ (mô men uốn, lực cắt, lực dọc) tại các tiết diện nguy hiểm.

Chương 3: Phương pháp tỷ số diện tích đánh giá mức độ an toàn kết cấu theo lý thuyết tập mờ

Chương 3 là phần đóng góp cốt lõi về mặt lý thuyết của luận án, tập trung phát triển và chứng minh tường minh công thức đánh giá mức độ an toàn theo phương pháp tỷ số diện tích:

Phương pháp đánh giá mờ Công thức tính toán Đánh giá ưu / nhược điểm
Giao thoa ngẫu nhiên - mờ [91], [78] $\tilde{P}f = \int{-\infty}^\infty \mu(x) f(x) dx$ Hai hàm không cùng độ đo; kết quả mang tính xấp xỉ
Lát cắt-$\alpha$ [107] $FP = h/2; \quad SP = 1 - h/2$ Chỉ xét tung độ giao điểm $h$, bỏ qua bề rộng đáy phần giao
Tỷ số giao hội [25] $FP = \Delta / (\tilde{R} + \tilde{Q}); \quad SP = 1 - FP$ Chưa tính trực tiếp $SP$; khó áp dụng khi miền giao bất thường
Độ tin cậy bậc nhất mờ (FFORM) [83] So sánh $\tilde{\beta} = [\beta_l, \beta_c, \beta_u]$ với $\beta_0$ Khối lượng tính toán rất lớn; đòi hỏi tập số liệu thống kê đầy đủ
Phương pháp tỷ số diện tích (Luận án đề xuất) $FRe = \frac{\int_{z \ge 0} \mu_{\tilde{Z}}(z) dz}{\int_{-\infty}^\infty \mu_{\tilde{Z}}(z) dz}$ Tính trực tiếp độ tin cậy từ tập hiệu $\tilde{Z} = \tilde{R} - \tilde{Q}$; xét toàn diện diện tích hàm thuộc
  • Triển khai toán học phương pháp tỷ số diện tích: Xét tập hiệu an toàn mờ $\tilde{Z} = \tilde{R} - \tilde{Q}$ với $\tilde{R}$ là khả năng chịu lực mờ và $\tilde{Q}$ là hiệu ứng tải trọng mờ:
    • Trường hợp tổng quát: Hàm thuộc $\mu_{\tilde{Z}}(z)$ có dạng bất kỳ, độ tin cậy mờ $FRe$ (Mức độ an toàn $SP$) được xác định bằng tỷ số giữa diện tích miền an toàn ($z \ge 0$) và tổng diện tích của miền xác định của $\tilde{Z}$: $$SP = FRe = \frac{S^+}{S^+ + S^-} = \frac{\int_{0}^{+\infty} \mu_{\tilde{Z}}(z) dz}{\int_{-\infty}^{+\infty} \mu_{\tilde{Z}}(z) dz}$$ Mức độ phá hoại (độ không tin cậy mờ $FP$): $$FP = \frac{S^-}{S^+ + S^-} = \frac{\int_{-\infty}^{0} \mu_{\tilde{Z}}(z) dz}{\int_{-\infty}^{+\infty} \mu_{\tilde{Z}}(z) dz} = 1 - SP$$
    • Trường hợp hàm thuộc tam giác: Khi $\tilde{Z} = (z_1, z_2, z_3)$ là số mờ tam giác, diện tích các phần tam giác và hình thang âm/dương được tính toán giải tích chính xác qua các tọa độ đỉnh $(z_1, z_2, z_3)$ cắt trục tung tại $z = 0$.
  • Phương pháp xây dựng hàm thuộc cho các tham số đầu vào: Trình bày phương pháp chuyên gia và phương pháp trực quan đơn giản hóa để xây dựng hàm thuộc mờ tam giác cho biên độ tải trọng động, tần số kích thích và hệ số cản nhớt kết cấu.

Chương 4: Mô hình tính kết cấu khung và đánh giá mức độ an toàn khung phẳng nhiều tầng chịu tải trọng động có tham số mờ

Chương 4 thực hiện các bài toán ứng dụng số nhằm kiểm chứng phương pháp đề xuất:

  • Thiết lập mô hình dao động khung nhiều tầng: Xét khung phẳng nhiều tầng chịu tải trọng động điều hòa $F_k(t) = P_{0k} \sin(\theta t)$ hoặc tải trọng xung. Các tham số đầu vào được mờ hóa gồm:
    • Khối lượng tầng: $\tilde{m}k = (m{k1}, m_{k2}, m_{k3})$.
    • Độ cứng tầng (phụ thuộc mô đun $\tilde{E}$ và liên kết nút mờ): $\tilde{K}_k$.
    • Ma trận cản nhớt: $[\tilde{C}] = \alpha [\tilde{M}] + \beta [\tilde{K}]$ (cản Rayleigh mờ).
    • Biên độ tải trọng $\tilde{P}_0$ và tần số kích động $\tilde{\theta}$.
  • Giải bài toán dao động mờ trên Maple.13:
    • Tách phương trình vi phân dao động thành hệ phương trình đại số thông qua phương pháp khai triển theo các dạng dao động riêng (modal analysis).
    • Áp dụng thuật toán tối ưu mức-$\alpha$ để giải nghiệm chuyển vị động cực đại mờ $\tilde{x}{max}$ và nội lực động cực đại mờ $\tilde{M}{max}, \tilde{Q}_{max}$.
  • Đánh giá mức độ an toàn và so sánh kết quả:
    • Thiết lập tập khả năng chịu lực cho phép $\tilde{R}$ (chuyển vị ngang đỉnh cho phép $[\tilde{f}]$ hoặc mô men giới hạn $[\tilde{M}_u]$).
    • Xác định tập hiệu mờ $\tilde{Z} = \tilde{R} - \tilde{Q}_{max}$.
    • Tính mức độ an toàn $SP$ và mức độ phá hoại $FP$ theo phương pháp tỷ số diện tích.
    • So sánh đối chiếu kết quả tính toán với các phương pháp: Phương pháp lát cắt-$\alpha$ [107], Phương pháp tỷ số giao hội [25] và Phương pháp giao thoa mờ tổng quát [96]. Kết quả so sánh cho thấy phương pháp tỷ số diện tích phản ánh sát sự làm việc thực tế, loại bỏ được các giả thiết gán ép tỷ lệ và xử lý ổn định cả trong các trường hợp hàm thuộc giao nhau bất thường.

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

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

  1. Hoàn thiện phương pháp tỷ số diện tích: Vận dụng, triển khai và chứng minh công thức toán học giải tích xác định mức độ an toàn (độ tin cậy mờ) $SP$ và mức độ phá hoại $FP$ của kết cấu dựa trên tỷ số diện tích hàm thuộc của tập hiệu mờ $\tilde{Z} = \tilde{R} - \tilde{Q}$. Công thức được chứng minh tổng quát cho trường hợp hàm thuộc bất kỳ và trường hợp cụ thể với hàm thuộc dạng tam giác.
  2. Quy trình giải hệ phương trình PTHH có tham số mờ: Đề xuất quy trình tính toán kết hợp phương pháp lát cắt-$\alpha$, thuật toán tối ưu mức-$\alpha$ và công cụ tính toán ký hiệu trên phần mềm Maple.13 để giải hệ phương trình đại số tuyến tính mờ, khắc phục hiện tượng tăng độ phân tán giả tạo của nghiệm trong số học mờ cổ điển.
  3. Mô hình động lực học kết cấu mờ toàn diện: Xây dựng thành công phương trình dao động và thuật giải đáp ứng động lực học của hệ khung phẳng nhiều tầng khi xét đồng thời các yếu tố mờ ở đầu vào gồm: đặc trưng vật liệu, độ cứng nút, độ cản nhớt, khối lượng tập trung và đặc trưng kích động của tải trọng động.

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

  1. Cung cấp một công cụ tính toán và đánh giá an toàn cho các kỹ sư thiết kế công trình khi đối mặt với điều kiện thiếu hụt số liệu thống kê hoặc thông tin đầu vào mang tính không chắc chắn.
  2. Cho phép người thiết kế đánh giá định lượng mức độ an toàn của khung nhà nhiều tầng chịu tải động (như rung chấn máy móc, gió động, xung lực), từ đó điều chỉnh các thông số thiết kế hình học và vật liệu nhằm đảm bảo công trình thỏa mãn tiêu chuẩn kỹ thuật và điều kiện tiện nghi sử dụng.
  3. Xây dựng các module chương trình tính toán số trên Maple.13 trong phần phụ lục, có thể sử dụng trực tiếp để giải các bài toán cơ học kết cấu có tham số mờ.

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

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

  • Luận án chỉ giới hạn khảo sát kết cấu khung phẳng làm việc trong miền đàn hồi tuyến tính; chưa xét đến ứng xử phi tuyến vật liệu và phi tuyến hình học.
  • Chưa xét đến tính tương tác động lực học giữa công trình và nền đất (Soil-Structure Interaction), giả thiết chân cột ngàm cứng tuyệt đối trên nền móng.
  • Các biến mờ đầu vào được giả thiết độc lập, chưa xét đến tính tương quan phụ thuộc giữa các tham số mờ.
  • Quy mô bài toán thực nghiệm số áp dụng ở mức độ kết cấu có số bậc tự do không lớn, đòi hỏi hàm truyền hoặc đáp ứng có dạng giải tích đóng để tối ưu hóa mức-$\alpha$.

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

  1. Mở rộng bài toán phân tích đánh giá mức độ an toàn cho hệ kết cấu không gian chịu tải trọng động ngẫu nhiên - mờ phức tạp (như tải trọng động đất, gió bão bất thường).
  2. Nghiên cứu bài toán phi tuyến động lực học có xét đến biến dạng dẻo và liên kết đàn hồi phi tuyến tại nút khung cũng như chân cột tiếp xúc nền đất.
  3. Phát triển thuật toán giải PTHH mờ cho các hệ kết cấu có số bậc tự do lớn bằng cách tích hợp thuật toán tối ưu hóa tiến hóa (Evolutionary Algorithms / Genetic Algorithms) với các phần mềm phần tử hữu hạn chuyên dụng.

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 thực tiễn:

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên ngành Xây dựng công trình:
    • Cung cấp cơ sở toán học chặt chẽ về lý thuyết tập mờ, số học mờ và phương pháp giải phương trình PTHH mờ (Chương 2).
    • Cung cấp phương pháp luận chi tiết và các chứng minh giải tích về đánh giá độ tin cậy mờ (Chương 3).
  • Đối với kỹ sư thiết kế và chuyên gia kiểm định công trình:
    • Cung cấp phương pháp đánh giá mức độ an toàn khung nhà nhiều tầng chịu tải trọng động trong điều kiện thông tin khảo sát không đầy đủ hoặc có độ nhậm nhằng cao (Chương 4).
    • Cung cấp các đoạn mã nguồn và quy trình tính toán trên phần mềm Maple.13 tại phần Phụ lục để áp dụng vào các bài toán tính toán kết cấu thực tế.

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

1. Sự khác biệt cốt lõi giữa mô hình đánh giá an toàn ngẫu nhiên và mô hình mờ trong luận án là gì?

Mô hình ngẫu nhiên giả định các đại lượng không chắc chắn là biến ngẫu nhiên tuân theo quy luật phân phối xác suất khách quan và đòi hỏi mẫu số liệu thống kê đủ lớn để xác định hàm mật độ $f(x)$. Ngược lại, mô hình mờ mô tả sự không chắc chắn do thiếu thông tin, thông tin không rõ ràng hoặc mang tính định tính thông qua hàm thuộc $\mu(x) \in [0, 1]$, không bắt buộc diện tích dưới đường cong hàm thuộc phải bằng đơn vị và không đòi hỏi quy luật phân phối thống kê chuẩn tắc.

2. Vì sao luận án không sử dụng nhóm phương pháp giao thoa ngẫu nhiên - mờ của các nghiên cứu trước đó?

Nhóm phương pháp giao thoa ngẫu nhiên - mờ (như công thức của Shiraishi & Furuta [91] hay Bùi Đức Chính [78]) tích phân tích số giữa hàm thuộc mờ $\mu(x)$ và hàm mật độ xác suất $f(x)$. Hai hàm này có bản chất và độ đo hoàn toàn khác nhau (diện tích của $f(x)$ với trục hoành bằng 1, trong khi diện tích của $\mu(x)$ khác 1), dẫn đến kết quả tính toán mang tính quy ước và gượng ép về mặt toán học. Luận án khắc phục điều này bằng cách chuyển hoàn toàn bài toán sang không gian mờ (mô hình mờ - mờ) qua tập hiệu $\tilde{Z} = \tilde{R} - \tilde{Q}$.

3. Phương pháp tỷ số diện tích được luận án chứng minh và tính toán như thế nào?

Phương pháp tỷ số diện tích định nghĩa mức độ an toàn $SP$ (chỉ số độ tin cậy mờ $FRe$) bằng tỷ số giữa diện tích phần dương $S^+$ (miền $z \ge 0$, tương ứng khả năng chịu lực lớn hơn hiệu ứng tải trọng) chia cho tổng diện tích $S^+ + S^-$ của đồ thị hàm thuộc tập hiệu $\tilde{Z} = \tilde{R} - \tilde{Q}$. Mức độ phá hoại được tính bằng $FP = S^- / (S^+ + S^-) = 1 - SP$. Đối với số mờ tam giác, các diện tích này được tính toán chính xác bằng công thức hình học giải tích.

4. Thuật toán tối ưu mức-$\alpha$ của Bernd Möller có ưu điểm gì so với thuật toán min-max truyền thống?

Thuật toán min-max thực hiện chia lưới rời rạc trên trục hoành (miền giá trị), dẫn đến số lượng tổ hợp tính toán bùng nổ và dễ gây sai số khuếch đại khoảng mờ. Thuật toán tối ưu mức-$\alpha$ chia lát cắt trên trục tung (mức độ thuộc $\alpha \in [0, 1]$) và giải hai bài toán quy hoạch tối ưu Min/Max tại mỗi lát cắt, giúp xác định chính xác cận dưới và cận trên của biến đầu ra mà không làm phân kỳ nghiệm.

5. Kết cấu khung khảo sát trong luận án được áp dụng những giả thiết động lực học nào?

Khung phẳng nhiều tầng được mô hình hóa thành hệ console có hữu hạn bậc tự do với các giả thiết: vật liệu làm việc trong giai đoạn đàn hồi tuyến tính; dầm và sàn tuyệt đối cứng trong mặt phẳng ngang; bỏ qua biến dạng dọc trục của cột và dầm; mỗi tầng có một bậc tự do chuyển vị ngang của khối lượng tập trung đặt tại mức sàn; liên kết chân cột ngàm cứng trên nền tuyệt đối.

Kết luận

Luận án của nghiên cứu sinh Lê Công Duy đã giải quyết trọn vẹn mục tiêu nghiên cứu về việc đánh giá mức độ an toàn của kết cấu khung phẳng chịu tải trọng động trong điều kiện các tham số đầu vào là số mờ. Tác giả đã hoàn thiện và chứng minh toán học tường minh phương pháp tỷ số diện tích, đồng thời đề xuất quy trình giải hệ phương trình PTHH mờ trên phần mềm Maple.13 dựa trên thuật toán tối ưu mức-$\alpha$. Kết quả nghiên cứu đóng góp cơ sở lý luận cho lĩnh vực cơ học kết cấu mờ và cung cấp công cụ tính toán tin cậy cho công tác thiết kế, kiểm định công trình dân dụng và công nghiệp.