Tổng quan luận án

Phân tích ổn định kết cấu và dao động phi tuyến là hai bài toán cơ bản trong lĩnh vực kỹ thuật công trình và cơ học ứng dụng. Trong thực tế tính toán thiết kế, các cấu kiện thanh, cột chịu nén và các hệ cơ học chịu tải trọng động thường biểu hiện các đặc tính phi tuyến phức tạp. Mặc dù bài toán ổn định của cột có tiết diện không đổi đã được nghiên cứu tương đối toàn diện, bài toán ổn định của cột có tiết diện thay đổi (dạng bậc thang, biến thiên liên tục theo hàm số mũ, hàm lũy thừa hoặc hàm đa thức) vẫn gặp nhiều khó khăn trong việc tìm kiếm lời giải giải tích đơn giản và chính xác.

Bên cạnh đó, các phương pháp số dư trọng số, tiêu biểu là phương pháp Galerkin, đóng vai trò nền tảng trong việc rời rạc hóa các phương trình vi phân chi phối. Tuy nhiên, độ chính xác của nghiệm xấp xỉ bậc nhất trong phương pháp Galerkin phụ thuộc chặt chẽ vào tính trực giao giữa phần dư phương trình và hàm thử thông qua toán tử lấy trung bình. Việc sử dụng toán tử lấy trung bình toàn cục thông thường trong phương pháp Galerkin kinh điển (GCA) hoặc phương pháp tuyến tính hóa tương đương (ELM) truyền thống thường làm gia tăng sai số khi bậc phi tuyến của hệ tăng cao hoặc khi cấu kiện có tiết diện biến thiên phức tạp. Khoảng trống nghiên cứu được đặt ra là cần cải tiến tính chất trực giao bằng cách biến đổi toán tử lấy trung bình, kết hợp các tiếp cận đối ngẫu và trung bình cục bộ có trọng số để nâng cao độ chính xác của nghiệm xấp xỉ bậc nhất mà không làm phức tạp hóa quy trình tính toán.

Mục tiêu nghiên cứu của luận án được xác định cụ thể gồm 4 nội dung:

  1. Phát triển tính chất trực giao của phương pháp số dư trọng số (cụ thể là phương pháp Galerkin) với trung bình thông thường và áp dụng tiêu chuẩn đối ngẫu có trọng số của phương pháp tuyến tính hóa tương đương cho bài toán ổn định đàn hồi.
  2. Giải quyết bài toán chọn một hàm trọng số cụ thể trong số lớp các hàm trọng số một tham số và phát triển tính chất trực giao bằng cách đề xuất xây dựng một trung bình cục bộ có trọng số mới ứng dụng vào phương pháp Galerkin và kết hợp với phương pháp bình phương tối thiểu khảo sát bài toán ổn định với thanh có tiết diện không đổi và có tiết diện thay đổi.
  3. Đề xuất một công cụ tính toán thay thế mới và hiệu quả để tính toán kỹ thuật trong việc thiết kế các hệ kết cấu với các tiết diện thay đổi.
  4. Xây dựng quy trình, phát triển tính chất trực giao thông qua phương pháp tuyến tính hóa đối ngẫu một trọng số cho các hệ động tiền định phi tuyến và áp dụng để phân tích tần số của dao động tự do phi tuyến một số trường hợp điển hình và so sánh kết quả một số phương pháp gần đúng để đánh giá hiệu quả của phương pháp đề xuất.

Đối tượng nghiên cứu của luận án là tính chất trực giao trong các phương pháp xấp xỉ gần đúng, tập trung vào phương pháp Galerkin, phương pháp tuyến tính hóa tương đương đối ngẫu và kỹ thuật lấy trung bình cục bộ có trọng số. Phạm vi nghiên cứu bao gồm:

  • Về bài toán ổn định: Khảo sát bài toán mất ổn định đàn hồi của cột Euler với 4 điều kiện biên cơ bản (hai đầu bản lề P-P, hai đầu ngàm C-C, một đầu ngàm một đầu tự do C-F, một đầu ngàm một đầu bản lề C-P) trong hai trường hợp tiết diện không đổi và tiết diện thay đổi (theo quy luật hàm số mũ, hàm lũy thừa, biến thiên tuyến tính và bậc hai).
  • Về bài toán dao động: Khảo sát các hệ cơ học một bậc tự do phi tuyến tiền định ở trạng thái dao động tự do, bao gồm hệ có lực phục hồi bậc phân số, dao động điều hòa Duffing, dao động phi tuyến có khả năng mở rộng hữu hạn và dao động kiểu Duffing.

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

Lịch sử nghiên cứu về bài toán ổn định đàn hồi được khởi xướng từ các thí nghiệm của Piter Musschenbroek (1729) và các phân tích toán học của Leonhard Euler về tải trọng tới hạn của thanh chịu nén dọc trục. Các phát triển lý thuyết và thực nghiệm tiếp theo trong thế kỷ 18 và 19 gắn liền với các công trình của Lagrange, Fr. Jasinski và M. Tổng quan lịch sử chi tiết về ổn định đàn hồi được tổng kết bởi Hoff (1956) và Timoshenko (1953). Các nghiên cứu về mất ổn định sau phân nhánh được thiết lập bởi Koiter (1945), Sewell, Thompson và Pope. Về mặt phân loại cơ học, Libove phân tích hiện tượng mất ổn định chuyển vị hữu hạn, trong khi Bolotin (1964) nghiên cứu ổn định động lực học dưới tải trọng biến thiên điều hòa theo thời gian. Các nghiệm giải tích chính xác về mất ổn định đàn hồi cho nhiều kết cấu cột, dầm, vòm, bản, vỏ có tiết diện biến thiên được Wang CM, Wang CY và Reddy tổng hợp có hệ thống.

Về các phương pháp gần đúng giải phương trình vi phân, phương pháp biến phân Rayleigh - Ritz dựa trên việc cực trị hóa phiếm hàm năng lượng. Phương pháp số dư trọng số (WRM) được khởi xướng bởi Bubnov (1913), Galerkin (1915) và được Duncan (1937) tổng thuật toàn diện. Mối quan hệ giữa phương pháp số dư trọng số và phương pháp biến phân được Finlayson và Scriven, cùng Leipholz phân tích rõ ràng. Trong phân tích động lực học và ổn định, phương pháp Galerkin được áp dụng và mở rộng trong nhiều công trình: Barati và Zenkour phân tích dao động vỏ composite nano xốp; Sofiyev nghiên cứu vỏ trụ FGM; Ding và cộng sự khảo sát sự hội tụ của phương pháp Galerkin cho dầm trên nền phi tuyến; Elishakoff và Zingales minh họa sự hội tụ chuỗi vô hạn; Laura và Cortinez khảo sát lũy thừa không nguyên trong hàm thử; Elishakoff và Damien Boutur chỉ ra các hạn chế của phương pháp Galerkin truyền thống khi áp dụng cho cột có tiết diện thay đổi giật cấp liên quan đến hàm Delta Dirac.

Trong lĩnh vực dao động phi tuyến, phương pháp tuyến tính hóa tương đương (ELM) bắt nguồn từ nghiên cứu của Kryloff và Bogoliubov (1937) cho hệ tiền định và được Caughey (1963) phát triển cho hệ ngẫu nhiên thông qua tiêu chuẩn cực tiểu hóa sai số bình phương trung bình. Hàng loạt công trình tổng quan của Spanos, Roberts, Socha và Soong, Falsone và Ricciardi, Elishakoff, Proppe, Crandall đã ghi nhận sự phát triển mạnh mẽ của ELM ngẫu nhiên. Đối với hệ tiền định, các tiêu chuẩn tuyến tính hóa được đóng góp bởi Adelberg và Denman, Sinha và Srinivasan, Agrwal và Denman, Chattopadhyay và Chakraborty, Beléndez và cộng sự. Các phương pháp tiệm cận và bán giải tích khác gồm phương pháp nhiễu loạn đồng vị HPM của He (1999), phương pháp cân bằng điều hòa HBM của Mickens, Lim và Wu, Febbo, cùng các tổng thuật của Cveticanin và Bayat.

Để khắc phục nhược điểm giảm độ chính xác khi tính phi tuyến tăng, Nguyễn Đông Anh (2015) đề xuất phương pháp lấy trung bình có trọng số thay thế trung bình cổ điển. Tiếp đó, khái niệm tiếp cận đối ngẫu trong tuyến tính hóa tương đương (kết hợp phép thay thế lượt đi từ hệ phi tuyến sang hệ tuyến tính và phép thay thế lượt về từ hệ tuyến tính sang hệ phi tuyến) được Nguyễn Đông Anh và các cộng sự phát triển thành tiêu chuẩn đối ngẫu có trọng số (WDC). Tiêu chuẩn này được Nguyễn Ngọc Linh áp dụng phân tích đáp ứng mô-men bậc hai cho dao động ngẫu nhiên.

Vị trí của luận án: Luận án kế thừa ý tưởng của phương pháp tuyến tính hóa tương đương đối ngẫu và tiêu chuẩn đối ngẫu có trọng số (WDC), chuyển hướng ứng dụng từ động lực học ngẫu nhiên sang bài toán giá trị biên ổn định đàn hồi và dao động tiền định phi tuyến. Luận án giải quyết khoảng trống trong việc sửa đổi toán tử lấy trung bình của phương pháp Galerkin thông qua kỹ thuật lấy trung bình cục bộ có trọng số (WLA), hình thành hai phương pháp tính toán mới là GWLA (Galerkin with weighted local averaging) và SGWLA (Simplified Galerkin with weighted local averaging), cho phép tính toán tường minh tải trọng tới hạn của cột có tiết diện biến thiên với độ chính xác cao dựa trên hàm thử đơn giản bậc nhất.

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

Khung lý thuyết và khái niệm sử dụng

  1. Lý thuyết ổn định đàn hồi và bài toán trị riêng: Ổn định đàn hồi được khảo sát qua phương trình vi phân uốn dọc của thanh. Xét trạng thái cân bằng phân nhánh, bài toán cân bằng đàn hồi được đưa về bài toán tìm giá trị riêng nhỏ nhất (tương ứng với tải trọng tới hạn $P_{cr}$) và hàm riêng (dạng mất ổn định).
  2. Phương pháp số dư trọng số và tính chất trực giao: Với toán tử vi phân $\mathcal{L}(u) = 0$ trên miền $\Gamma$, nghiệm xấp xỉ $\tilde{u} = \sum_{i=1}^n a_i g_i(x)$ tạo ra phần dư $R_\Gamma$. Điều kiện triệt tiêu phần dư được biểu diễn qua tích phân trực giao: $$\int_\Gamma W_i R_\Gamma d\Gamma = 0$$ Trong phương pháp Galerkin, hàm trọng số chính là đạo hàm của hàm xấp xỉ theo các bậc tự do: $W_i = \partial\tilde{u}/\partial a_i = g_i(x)$.
  3. Ý tưởng tiếp cận đối ngẫu và tiêu chuẩn đối ngẫu có trọng số (WDC): Dựa trên bài toán tối ưu hóa đa mục tiêu (MOO) theo phương pháp tổng có trọng số, quá trình thay thế tương đương giữa hàm phi tuyến và hàm tuyến tính được thiết lập đồng thời theo hai chiều:
    • Chiều đi: thay thế hàm phi tuyến bằng hàm tuyến tính với sai số $e_{df}$.
    • Chiều về: thay thế hàm tuyến tính bằng hàm phi tuyến với sai số $e_{dr}$. Hàm mục tiêu tổ hợp vô hướng được lập thông qua trọng số chuẩn hóa $p \in [0, 1]$: $$WDC = p \cdot E[e_{df}^2] + (1-p) \cdot E[e_{dr}^2] \to \min$$
  4. Kỹ thuật lấy trung bình cục bộ có trọng số (WLA): Thay vì tích phân toàn cục trên toàn bộ chiều dài miền $[0, l]$, toán tử lấy trung bình được thực hiện cục bộ tại một tọa độ không thứ nguyên $r \in [0, 1]$ với hàm cục bộ $r f(r)$, sau đó lấy tích phân có trọng số để thu được giá trị tương đương đại diện.

Phương pháp nghiên cứu thực tế

  • Phương pháp giải tích và hình học giải tích: Sử dụng các phép biến đổi toán học giải tích để thiết lập biểu thức thế năng, phương trình vi phân phần dư, dẫn xuất các công thức tính hệ số tương đương và tần số dao động. Mô hình hóa hình học giải tích ý tưởng đối ngẫu thông qua các phép biến đổi hình học (như phép thay thế đối ngẫu giữa hình tròn và hình vuông).
  • Phương pháp số và công cụ tính toán: Luận án sử dụng các lệnh tính toán số và vẽ đồ thị trên phần mềm MATLAB, đồng thời lập trình các thuật toán giải tích biến phân và giải hệ phương trình đại số phi tuyến trên phần mềm Mathematica.

Nguồn tư liệu và dữ liệu kiểm chứng

Nguồn dữ liệu đối chứng bao gồm các nghiệm giải tích chính xác đã được công bố trong các công trình kinh điển của Leonhard Euler, S.P. Timoshenko, C.M. Wang, C.Y. Wang và J.N. Reddy cho bài toán ổn định cột; nghiệm giải tích và chuỗi tiệm cận của Cveticanin, Mickens, He, Febbo cho các hệ dao động phi tuyến.

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

Chương 1: Tổng quan về phương pháp giải bài toán ổn định đàn hồi và dao động phi tuyến hệ một bậc tự do

Chương 1 hệ thống hóa các khái niệm nền tảng về ổn định và dao động:

  • Các mô hình mất ổn định: Phân tích ba dạng mất ổn định cơ bản:
    • Mất ổn định loại 1 (mất ổn định phân nhánh hay cổ điển): Kết cấu chuyển từ trạng thái cân bằng ổn định sang trạng thái cân bằng không ổn định lân cận khi tải trọng đạt giá trị tới hạn $P_{cr}$.
    • Mất ổn định loại 2 (snapthrough buckling): Xuất hiện bước nhảy đột ngột từ trạng thái cân bằng này sang trạng thái cân bằng khác không liền kề có chuyển vị lớn hơn.
    • Mất ổn định nhiễu loạn hữu hạn (chuyển vị hữu hạn): Kết cấu mất độ cứng nghiêm trọng, trạng thái mất ổn định chỉ có thể duy trì khi giảm tải trọng.
  • Bốn bài toán ổn định đàn hồi cột Euler tiết diện không đổi:
    • Thanh hai đầu bản lề (P-P): $P_{cr} = \pi^2 EI / l^2 \approx 9,8696 EI/l^2$.
    • Thanh hai đầu ngàm (C-C): $P_{cr} = 4\pi^2 EI / l^2 \approx 39,4784 EI/l^2$.
    • Thanh một đầu ngàm, một đầu tự do (C-F): $P_{cr} = \pi^2 EI / (4l^2) \approx 2,4674 EI/l^2$.
    • Thanh một đầu ngàm, một đầu bản lề (C-P): $P_{cr} \approx 20,1906 EI/l^2$.
  • Mô hình dao động phi tuyến: Thiết lập phương trình vi phân tổng quát của hệ một bậc tự do: $$\ddot{x} + 2h\dot{x} + \omega_0^2 x + g(x, \dot{x}) = f(t)$$ Điểm qua 10 hệ dao động phi tuyến điển hình: hệ Lutes Sarkani, Van der Pol, hệ có cản phi tuyến bậc ba, hệ đàn hồi phi tuyến theo quy luật mũ, hệ Duffing (với đường cong thế năng giếng đơn và giếng đôi), hệ Atlik-Utku, dao động có lực phục hồi bậc phân số, dao động điều hòa Duffing, dao động phi tuyến có khả năng mở rộng hữu hạn và dao động kiểu Duffing.
  • So sánh các phương pháp số dư trọng số (WRM): Khảo sát bài toán thanh đàn hồi tiết diện không đổi chịu lực phân bố tuyến tính $q = cx$ và lực dọc $P=0$ ($A=1, E=1, c=1, L_T=1$). Kết quả so sánh chuyển vị $u$ và đạo hàm $u'$ được tổng hợp tại Bảng 1.1 trong luận án:
Vị trí khảo sát Nghiệm chính xác Phương pháp ấn định điểm Phương pháp miền phụ và bình phương tối thiểu Phương pháp ấn định điểm bình phương tối thiểu Phương pháp Galerkin
$u$ tại $x = L_T/2$ 0,2292 0,1250 0,1875 0,2500 0,2292
$u$ tại $x = L_T$ 0,3333 0,1667 0,2500 0,3333 0,3333
$u$ tại $x = 0$ 0,5000 0,3333 0,5000 0,6667 0,5833
$u'$ tại $x = L_T/2$ 0,3750 0,1667 0,2500 0,3333 0,3333
$u'$ tại $x = L_T$ 0 0 0 0 0,0833
  • Đánh giá phương pháp Galerkin: Ứng dụng giải bài toán cột C-P với 7 hàm thử chuyển vị $W(x)$ khác nhau thỏa mãn điều kiện biên, cho thấy sai số của lực tới hạn $P_{Galerkin}$ so với nghiệm Euler ($20,1906$) dao động từ $4,26%$ đến $134,05%$. Điều này khẳng định độ chính xác của phương pháp Galerkin phụ thuộc rất lớn vào dạng hàm thử và toán tử trực giao.

Chương 2: Phương pháp tuyến tính hóa tương đương

Chương 2 trình bày có hệ thống về cơ sở giải tích của phương pháp tuyến tính hóa tương đương (ELM) và bước phát triển sang tiêu chuẩn đối ngẫu:

  • Tuyến tính hóa tương đương cho hệ một bậc tự do: Thay thế thành phần phi tuyến $g(x, \dot{x})$ bằng các thành phần tuyến tính tương đương $b\dot{x} + kx$. Hệ số cản $b$ và độ cứng $k$ được xác định thông qua việc cực tiểu hóa sai số bình phương trung bình $e_1 = E[(g(x, \dot{x}) - b\dot{x} - kx)^2] \to \min$, dẫn đến: $$b = \frac{E[\dot{x}g]}{E[\dot{x}^2]}, \quad k = \frac{E[xg]}{E[x^2]}$$ Hệ phương trình được đóng kín bằng phương trình Fokker - Planck - Kolmogorov (FPK) với đáp ứng bình phương trung bình: $$\langle x^2 \rangle = \frac{\sigma^2}{2(2h+b)(\omega_0^2+k)}, \quad \langle \dot{x}^2 \rangle = \frac{\sigma^2}{2(2h+b)}$$
  • Các biến thể tiêu chuẩn tuyến tính hóa: Phân tích tiêu chuẩn cực tiểu sai số thế năng, tiêu chuẩn tuyến tính hóa có điều chỉnh và tiêu chuẩn tuyến tính hóa từng phần.
  • Tiêu chuẩn đối ngẫu và tiêu chuẩn đối ngẫu có trọng số (WDC): Trình bày bản chất hình học và giải tích của phép thay thế hai chiều. Sử dụng hệ số tương quan bình phương $r^2$ để phản ánh mức độ phi tuyến giữa hàm gốc và hàm thay thế, làm cơ sở toán học để thiết lập trọng số kết nối hai hàm mục tiêu trong tối ưu hóa đa mục tiêu.

Chương 3: Phát triển tính chất trực giao áp dụng trong phân tích bài toán ổn định

Chương 3 là một trong hai chương đóng góp học thuật trọng tâm của luận án, tập trung phát triển tính chất trực giao cho bài toán ổn định đàn hồi:

  • Áp dụng tiêu chuẩn đối ngẫu có trọng số (WDC) cho bài toán ổn định: Khảo sát bài toán xác định lực tới hạn cho 4 dạng liên kết cột (P-P, C-C, C-F, C-P). Khảo sát sự thay đổi của lực tới hạn $P_{crdn}$ khi thay đổi hàm thử và các hàm trọng số $p_i(r)$.
  • Xây dựng kỹ thuật lấy trung bình cục bộ có trọng số (WLA): Đề xuất toán tử lấy trung bình cục bộ tại điểm $r$ với hàm $r f(r)$ trên đoạn $[0, 1]$.
  • Phương pháp GWLA và SGWLA:
    • Tích hợp kỹ thuật WLA vào phương pháp Galerkin tạo thành phương pháp GWLA.
    • Đơn giản hóa quy trình tính toán để hình thành phương pháp SGWLA, cho phép xác định trực tiếp tải trọng tới hạn chuẩn hóa $P_{SGWLA}^{const}$ mà chỉ cần dùng một hàm thử đơn giản.
  • Khảo sát mất ổn định của cột Euler có tiết diện không đổi: So sánh kết quả tải trọng tới hạn giữa phương pháp Galerkin thông thường ($P_{GCA}^{const}$), GWLA ($P_{GWLA}^{const}$), SGWLA ($P_{SGWLA}^{const}$) với nghiệm chính xác ($P_{exact}^{const}$) trên cả 4 sơ đồ liên kết cột (P-P, C-P, C-C, C-F). Kết quả cho thấy GWLA và SGWLA giảm thiểu đáng kể sai số so với GCA khi sử dụng cùng một hàm thử bậc thấp.
  • Khảo sát mất ổn định của cột Euler có tiết diện thay đổi:
    • Đưa ra quy trình chuyển đổi cột có tiết diện biến thiên thành cột tương đương có tiết diện không đổi thông qua các hệ số độ cứng tương đương $k_{CA}^{var}$ và $k_{WLA}^{var}$.
    • Xét hai quy luật biến thiên tiết diện chính: mô-men quán tính biến thiên theo hàm số mũ ($k_{CA}^{exp}, k_{WLA}^{exp}$) và theo hàm lũy thừa ($k_{CA}^{pow}, k_{WLA}^{pow}$).
    • Khảo sát chi tiết cho trường hợp tiết diện biến thiên tuyến tính ($a=1$) và biến thiên bậc hai ($a=2$).
    • Kết quả tính toán chứng minh tải trọng tới hạn $P_{GWLA}^{exp, WLA}$ và $P_{GWLA}^{pow, WLA}$ hội tụ sát với nghiệm chính xác $P_{exact}^{exp}$ và $P_{exact}^{pow}$ hơn rõ rệt so với phương pháp GCA truyền thống.

Chương 4: Phát triển tính chất trực giao áp dụng trong phân tích bài toán dao động phi tuyến

Chương 4 mở rộng việc phát triển tính chất trực giao và phương pháp tuyến tính hóa đối ngẫu sang phân tích tần số dao động tự do của các hệ động lực học phi tuyến tiền định:

  • Quy trình tuyến tính hóa đối ngẫu cho hệ tiền định: Thiết lập tiêu chuẩn đối ngẫu có trọng số cho phương trình vi phân phi tuyến thuần túy và phi tuyến bậc cao. Lựa chọn trọng số $p_n$ và $1-p_n$ liên hệ với hệ số tương quan $r^2$.
  • Phân tích số và kiểm chứng trên 4 bài toán dao động tự do phi tuyến:
    1. Dao động với lực phục hồi bậc phân số ($\ddot{x} + c^2 x |x|^{n-1} = 0$, $c > 0, n > 1$): Xác định tần số xấp xỉ và đánh giá sai số theo các bậc phân số khác nhau.
    2. Dao động điều hòa Duffing ($\ddot{x} + \frac{x^3}{1+x^2} = 0$): Khảo sát tần số dao động với lực phục hồi hữu tỉ chưa tối giản.
    3. Dao động phi tuyến có khả năng mở rộng hữu hạn ($\ddot{x} + \frac{x}{1-x^2} = 0$): Khảo sát hệ phi tuyến mạnh có lực phục hồi dạng phân thức tối giản.
    4. Dao động kiểu Duffing ($\ddot{x} + \omega_0^2 x + \gamma x^{2n+1} = 0$): Khảo sát với các bậc phi tuyến lẻ $n = 1, 2, 3$ với các giá trị biên độ $a$, tần số riêng tuyến tính $\omega_0 = 1$ và hệ số phi tuyến $\gamma$ thay đổi.
  • Đánh giá độ chính xác: So sánh tần số xấp xỉ thu được từ tiêu chuẩn đối ngẫu có trọng số ($\omega_{d2}$) với nghiệm chính xác ($\omega_e$), tiêu chuẩn kinh điển ($\omega_c$), tiêu chuẩn trung bình có trọng số ($\omega_{aw}$), phương pháp tiệm cận ($\omega_{as}$), phương pháp tham số mở rộng ($\omega_{PEM}$), phương pháp cân bằng điều hòa bậc thấp nhất ($\omega_{lhb}$) và phương pháp cân bằng điều hòa tuyến tính ($\omega_{hb}$). Kết quả cho thấy phương pháp đề xuất duy trì sai số nhỏ và ổn định ngay cả khi biên độ dao động hoặc hệ số phi tuyến $\gamma$ tăng cao.

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

Luận án tự xác định các đóng góp mới về mặt khoa học và kỹ thuật bao gồm:

Về mặt lý luận

  1. Phát triển thành công tính chất trực giao trong phương pháp số dư trọng số (cụ thể là phương pháp Galerkin) bằng việc thay thế toán tử lấy trung bình toàn cục thông thường bằng phép lấy trung bình cục bộ có trọng số (WLA).
  2. Thiết lập cơ sở toán học cho phương pháp Galerkin với phép lấy trung bình cục bộ có trọng số (GWLA) và phương pháp Galerkin đơn giản hóa (SGWLA), giải quyết bài toán lựa chọn hàm trọng số tối ưu trong lớp hàm một tham số thông qua hệ số tương quan $r^2$.
  3. Mở rộng phạm vi ứng dụng của phương pháp tuyến tính hóa tương đương đối ngẫu và tiêu chuẩn đối ngẫu có trọng số (WDC) từ các bài toán ngẫu nhiên sang bài toán giá trị biên ổn định đàn hồi và bài toán dao động tự do phi tuyến tiền định.

Về mặt thực tiễn và ứng dụng tính toán

  1. Đề xuất một công cụ tính toán thay thế giải tích gần đúng mới, cho phép xác định tải trọng tới hạn của cột Euler có tiết diện không đổi và tiết diện thay đổi (hàm mũ, hàm lũy thừa, bậc nhất, bậc hai) với độ chính xác cao chỉ bằng việc sử dụng một hàm thử đơn giản bậc nhất, tránh được việc phải giải các hệ phương trình vi phân phức tạp hoặc chia nhỏ phần tử hữu hạn.
  2. Xây dựng quy trình tuyến tính hóa đối ngẫu từng bước áp dụng để xác định tần số dao động tự do cho các hệ phi tuyến phức tạp (lực phục hồi phân số, hàm hữu tỉ, phi tuyến bậc cao), cho kết quả giải tích tường minh với độ chính xác vượt trội so với phương pháp tuyến tính hóa kinh điển.

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

Dựa trên phạm vi và nội dung được trình bày trong văn bản luận án, các hạn chế và định hướng nghiên cứu tiếp theo bao gồm:

  • Hạn chế:
    • Luận án tập trung chủ yếu vào các cấu kiện dạng thanh/cột một chiều (1D) và các hệ cơ học một bậc tự do (1-DOF) chịu kích động tiền định hoặc trạng thái tự do.
    • Chưa mở rộng khảo sát cho các kết cấu hai chiều phức tạp hơn như tấm, vỏ mỏng có tiết diện thay đổi hoặc hệ nhiều bậc tự do (M-DOF) liên tục chịu tải trọng ngẫu nhiên phi dừng.
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng phát triển tính chất trực giao và kỹ thuật WLA cho bài toán ổn định và dao động của kết cấu tấm và vỏ mỏng đàn hồi có cơ tính biến thiên (FGM) hoặc vật liệu composite phân lớp.
    • Ứng dụng phương pháp tuyến tính hóa đối ngẫu có trọng số cho các hệ động lực học phi tuyến nhiều bậc tự do và hệ liên tục chịu kích động ngẫu nhiên phức tạp.
    • Nghiên cứu kết hợp kỹ thuật GWLA/SGWLA vào thuật toán của phương pháp phần tử hữu hạn (FEM) để phát triển các phần tử dầm - cột đặc biệt có tiết diện biến thiên liên tục.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong lĩnh vực Cơ kỹ thuật, Kỹ thuật Cơ khí và Kỹ thuật Xây dựng - Công trình:

  • Đối với nghiên cứu sinh, học viên cao học và nhà nghiên cứu cơ học: Luận án cung cấp khung phương pháp luận hoàn chỉnh về kỹ thuật lấy trung bình cục bộ có trọng số (WLA) và phương pháp tuyến tính hóa đối ngẫu (WDC). Các phân tích so sánh chi tiết ở Chương 1, Chương 3 và Chương 4 là tài liệu tham khảo hữu ích cho việc phát triển các phương pháp giải tích xấp xỉ mới.
  • Đối với giảng viên và người giảng dạy: Cung cấp hệ thống ví dụ minh họa chi tiết, các bảng so sánh sai số và phân tích bản chất cơ học của 4 bài toán biên cột Euler, các mô hình mất ổn định phân nhánh, snapthrough và 10 mô hình dao động phi tuyến điển hình.
  • Đối với kỹ sư tính toán kết cấu công trình: Cung cấp các công thức giải tích đơn giản hóa (SGWLA) để ước lượng nhanh tải trọng tới hạn của các cột có tiết diện thay đổi trong thiết kế kết cấu xây dựng, cầu đường và tháp công trình mà không đòi hỏi phần mềm số phức tạp.

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

1. Phương pháp GWLA và SGWLA khác biệt như thế nào so với phương pháp Galerkin thông thường (GCA)?

Phương pháp Galerkin thông thường (GCA) sử dụng toán tử lấy trung bình tích phân toàn cục trên toàn miền khảo sát, khiến sai số phụ thuộc chặt chẽ vào việc lựa chọn hàm thử; nếu dùng hàm thử đơn giản, sai số thường lớn. Phương pháp GWLA thay thế toán tử toàn cục bằng phép lấy trung bình cục bộ có trọng số (WLA) thông qua hàm cục bộ $r f(r)$. Phương pháp SGWLA là phiên bản đơn giản hóa của GWLA, cho phép tính toán trực tiếp tải trọng tới hạn với độ chính xác cao chỉ bằng một hàm thử bậc thấp duy nhất.

2. Tiêu chuẩn đối ngẫu có trọng số (WDC) xử lý bài toán tuyến tính hóa tương đương theo nguyên lý nào?

WDC tiếp cận bài toán tuyến tính hóa tương đương dưới dạng bài toán tối ưu hóa đa mục tiêu (MOO). Thay vì chỉ xét phép thay thế một chiều từ hệ phi tuyến sang hệ tuyến tính (chiều đi với sai số $e_{df}$), WDC kết hợp đồng thời phép thay thế ngược lại từ hệ tuyến tính sang hệ phi tuyến (chiều về với sai số $e_{dr}$). Hai hàm mục tiêu sai số bình phương trung bình được tổ hợp vô hướng thông qua trọng số chuẩn hóa $p$, liên kết chặt chẽ với hệ số tương quan bình phương $r^2$.

3. Luận án khảo sát bài toán ổn định cột Euler với những dạng tiết diện thay đổi nào?

Luận án khảo sát cột Euler có mô-men quán tính thay đổi theo hai quy luật giải tích chính: quy luật hàm số mũ và quy luật hàm lũy thừa. Cụ thể, các trường hợp tiết diện biến thiên tuyến tính ($a=1$) và biến thiên bậc hai ($a=2$) được giải tích chi tiết trên cả 4 dạng điều kiện biên (P-P, C-C, C-F, C-P).

4. Bốn hệ dao động phi tuyến tiền định được phân tích trong Chương 4 gồm những hệ nào?

Bốn bài toán dao động tự do phi tuyến tiền định được khảo sát bao gồm:

  1. Dao động với lực phục hồi bậc phân số: $\ddot{x} + c^2 x |x|^{n-1} = 0$.
  2. Dao động điều hòa Duffing có lực phục hồi là hàm hữu tỉ chưa tối giản: $\ddot{x} + \frac{x^3}{1+x^2} = 0$.
  3. Dao động phi tuyến có khả năng mở rộng hữu hạn với lực phục hồi là hàm hữu tỉ tối giản: $\ddot{x} + \frac{x}{1-x^2} = 0$.
  4. Dao động kiểu Duffing với lực đàn hồi phi tuyến bậc lẻ: $\ddot{x} + \omega_0^2 x + \gamma x^{2n+1} = 0$ ($n = 1, 2, 3$).

Kết luận

Luận án tiến sĩ của tác giả Trần Tuấn Long đã nghiên cứu phát triển thành công tính chất trực giao trong phương pháp số dư trọng số và phương pháp tuyến tính hóa tương đương thông qua việc tích hợp tiêu chuẩn đối ngẫu có trọng số (WDC) và kỹ thuật lấy trung bình cục bộ có trọng số (WLA). Các phương pháp mới được đề xuất (GWLA, SGWLA và quy trình tuyến tính hóa đối ngẫu tiền định) đã giải quyết hiệu quả bài toán xác định tải trọng tới hạn của cột đàn hồi có tiết diện biến thiên phức tạp và phân tích tần số dao động tự do của các hệ động lực học phi tuyến. Kết quả nghiên cứu cung cấp một công cụ giải tích gần đúng có độ chính xác cao, đóng góp hoàn thiện lý thuyết xấp xỉ trong cơ học kết cấu và kỹ thuật công trình.