Tổng quan luận án

Luận án tiến sĩ cơ học của tác giả Nguyễn Như Hiếu với tiêu đề "Tiêu chuẩn đối ngẫu trong phương pháp tuyến tính hóa tương đương cho hệ phi tuyến nhiều bậc tự do chịu kích động ngẫu nhiên", chuyên ngành Cơ học vật rắn (mã số 62440107), được hoàn thành năm 2018 tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội dưới sự hướng dẫn khoa học của GS. Nguyễn Đông Anh và PGS. Ninh Quang Hải.

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

Phương pháp tuyến tính hóa tương đương là một trong những công cụ giải tích xấp xỉ phổ biến và lâu đời nhất trong tính toán động lực học ngẫu nhiên. Phương pháp này thay thế các phần tử phi tuyến bằng các phần tử tuyến tính tương đương nhằm đơn giản hóa việc giải phương trình vi phân chuyển động dưới tác dụng của kích động ngẫu nhiên. Tuy nhiên, các nghiên cứu tuyến tính hóa tương đương truyền thống (kinh điển) dựa trên tiêu chuẩn cực tiểu hóa sai số bình phương trung bình hầu như chỉ đạt độ chính xác chấp nhận được khi áp dụng cho các hệ có tính phi tuyến yếu. Trong trường hợp hệ có tính phi tuyến mạnh, phương pháp truyền thống dẫn đến sai số rất lớn so với nghiệm chính xác hoặc nghiệm mô phỏng số. Do đó, yêu cầu cải tiến phương pháp để nâng cao độ chính xác khi phân tích đáp ứng của các hệ nhiều bậc tự do có tính phi tuyến mạnh chịu kích động ngẫu nhiên là một vấn đề cấp thiết, tạo nên khoảng trống nghiên cứu mà luận án tập trung giải quyết.

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

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

  1. Tăng được độ chính xác của phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu cho hệ ngẫu nhiên phi tuyến nhiều bậc tự do.
  2. Xây dựng công thức xác định các hệ số tương đương của phương trình tuyến tính hóa theo tiêu chuẩn đối ngẫu.
  3. Áp dụng phương pháp cho việc khảo sát đáp ứng của một số hệ liên tục mà phương trình của nó có thể đưa về dạng phương trình của hệ nhiều bậc tự do.

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

  • Đối tượng nghiên cứu: Tiêu chuẩn đối ngẫu trong phương pháp tuyến tính hóa tương đương áp dụng cho các hệ cơ học nhiều bậc tự do có tính phi tuyến mạnh chịu kích động ngoài ngẫu nhiên; các hàm phi tuyến trong hệ là hàm liên tục của các biến trạng thái (vị trí, vận tốc, gia tốc).
  • Phạm vi nghiên cứu:
    • Về mặt lý thuyết và giải tích: Giới hạn tính toán giải tích chi tiết cho các mô hình hệ tập trung hai bậc tự do (độ cứng phi tuyến bậc ba, độ cứng phi tuyến bậc năm, cản và độ cứng phi tuyến đồng thời) và các bài toán kết cấu hệ liên tục (dầm Euler-Bernoulli tựa giản đơn trên nền đàn hồi chịu tải ngẫu nhiên phân bố, dầm mang khối lượng tập trung) được đưa về hệ phương trình vi phân thường cho các biên độ mode dao động (đơn mode, hai mode và nhiều mode lên đến 5 mode).
    • Về kích động: Kích động ngoài được khảo sát là quá trình ngẫu nhiên Gaussian ồn trắng dừng, có giá trị kỳ vọng toán học bằng không.

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

Các hướng nghiên cứu trước đây được luận án điểm lại

Luận án hệ thống hóa lịch sử phát triển của phương pháp tuyến tính hóa tương đương và các kỹ thuật xấp xỉ trong cơ học ngẫu nhiên thông qua các công trình tiêu biểu:

  • Phương pháp tuyến tính hóa tương đương kinh điển: Được đề xuất đồng thời bởi Booton (1954), Kazakov (1954) và Caughey (1956) trong thập niên 1950 dựa trên tiêu chuẩn cực tiểu hóa sai số bình phương trung bình. Sau đó, phương pháp được mở rộng cho hệ nhiều bậc tự do bởi Foster (1968), Iwan và Yang (1972), Atalik và Utku (1976), cùng các tổng kết chuyên khảo của Roberts và Spanos (1990), Socha (2008). Atalik và Utku (1976) đã chứng minh định lý quan trọng cho phép biểu diễn ma trận hệ số tuyến tính hóa qua kỳ vọng của gradient của véc tơ phi tuyến khi các biến trạng thái có phân bố chuẩn.
  • Các phương pháp giải tích xấp xỉ khác: Phương pháp nhiễu (Crandall 1963, Schuëller 1996), phương pháp mô men (Gard 1988), phương pháp trung bình ngẫu nhiên (Spanos 1992), phương pháp khai triển theo hàm riêng (Wong 1964, Payne 1967, Atkinson 1970). Tuy nhiên, các kỹ thuật này thường chỉ giới hạn ở hệ một bậc tự do và có quy trình giải tích phức tạp.
  • Tuyến tính hóa theo tiêu chuẩn năng lượng: Zhang và đồng nghiệp (1990) đề xuất cực tiểu hóa sai số bình phương trung bình của thế năng; Elishakoff và Zhang (1992) nghiên cứu tiêu chuẩn đẳng thức năng lượng; Elishakoff và Colombi (1993), Falsone và Elishakoff (1994) phát triển tiêu chuẩn dựa trên hàm hao tán năng lượng và hàm trọng; Colajani và Elishakoff (1998) hiệu chỉnh tiêu chuẩn cho các hệ có lực đàn hồi dạng hyperbolic tangent.
  • Tuyến tính hóa phi tham số (TTHPTS): Fujimura và Der Kiureghian (2007) đề xuất phương pháp tuyến tính hóa phi tham số bám sát các điểm thiết kế trong không gian biến chuẩn, tiếp tục được phát triển bởi Luca và Der Kiureghian (2010), Reza và Mohsen (2014, 2016), Broccardo và Armen (2015), Wang và Armen (2016).
  • Tuyến tính hóa tương đương điều chỉnh: Phương pháp điều chỉnh một bước của Nguyễn Đông Anh và Di Paola (1995), và điều chỉnh hai bước của Elishakoff, Andrimasy và Dolley (2009).
  • Ứng dụng cho hệ trễ và kết cấu công trình: Casciati và Faravelli (1993), Hurtado và Barbat (2000), Silva và Ruiz (2000), Spanos và Agathoklis (2013) với mô hình Bouc-Wen; Colangelo (2017) cho tương tác lực dọc và mô men uốn; Günay và Sucuoğlu (2009), Jara và đồng nghiệp (2012), Liu và đồng nghiệp (2014) trong phân tích kháng chấn và cách ly rung động; Pradlwarter và Li (1991), Sun và đồng nghiệp (1998), Dogan và Vaicaitis (1999) cho tấm, vỏ, dầm và khung.
  • Nghiên cứu nghiệm chính xác từ phương trình Fokker-Planck: Andronov (1933), Caughey và Dienes (1961), Crandall (1962, 1964, 2006), Stratonovich (1963) với khái niệm thế năng dừng, Cai và Lin (1988) với phương pháp thế năng dừng tổng quát hóa, Wen-Yao và Tong (1992).
  • Nghiên cứu tại Việt Nam: Nguyễn Đông Anh (Anh và Schiehlen 1999; Anh 2010; Anh và Hieu 2012), Ninh Quang Hải (2000), Lưu Xuân Hùng (2002), Đặng Văn Hiếu (2010), Nguyễn Như Hiếu (2011), Nguyễn Ngọc Linh (2014), Dương Ngọc Hảo (2014), cùng nghiên cứu ứng dụng giảm chấn của Nguyễn Đông Anh và Lã Đức Việt (2007).
  • Tiêu chuẩn đối ngẫu ban đầu: Ý tưởng tiêu chuẩn đối ngẫu được khởi xướng bởi Nguyễn Đông Anh (2010) và nhóm cộng sự (2012a, 2012b, 2012d), nhưng mới chỉ dừng lại ở các hệ dao động phi tuyến một bậc tự do.

Khoảng trống luận án lựa chọn giải quyết

Trong khi tiêu chuẩn đối ngẫu đã khẳng định ưu thế giảm sai số trên các hệ một bậc tự do, lý thuyết này hoàn toàn chưa được thiết lập cho các hệ động lực ngẫu nhiên nhiều bậc tự do và hệ kết cấu liên tục. Luận án lựa chọn giải quyết khoảng trống này bằng cách: thiết lập cơ sở toán học cho tiêu chuẩn đối ngẫu trên hệ $n$ bậc tự do; xây dựng hệ phương trình ma trận đóng kín để tìm đồng thời các ma trận hệ số tuyến tính hóa; mở rộng giải thuật cho dao động uốn của dầm liên tục có/không có nền đàn hồi và dầm mang tải trọng tập trung; đồng thời đánh giá độ chính xác qua việc so sánh đối chứng với nghiệm chính xác từ phương trình Fokker-Planck và nghiệm mô phỏng Monte-Carlo.

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

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

Luận án vận dụng các lý thuyết và công cụ giải tích toán học trong lý thuyết xác suất và cơ học ngẫu nhiên:

  • Lý thuyết quá trình ngẫu nhiên Gaussian: Hàm đặc trưng véc tơ, các đặc trưng mô men bậc nhất và bậc hai, tính chất triệt tiêu của các mô men bậc lẻ, quy tắc phân rã mô men bậc chẵn thành tổ hợp các cặp tương quan, và đẳng thức vi phân cho các biến chuẩn (tính chất 1 đến 6 theo Socha 2008).
  • Quá trình Markov và quá trình Wiener: Quá trình Wiener với các số gia độc lập Gaussian dừng; định nghĩa quá trình ngẫu nhiên ồn trắng như đạo hàm suy rộng của quá trình Wiener trên không gian hàm thử $D(T)$, có hàm tự tương quan dạng Dirac-delta $K(t_1, t_2) = c\delta(t_2 - t_1)$ và mật độ phổ năng lượng không đổi $S(\omega) = c$.
  • Khái niệm hệ phi tuyến yếu và mạnh: Phân loại định lượng theo Babitsky và Krupenin (2001) dựa trên tham số bé $\varepsilon$. Hệ được xem là phi tuyến yếu khi hệ số phi tuyến có độ lớn cỡ $O(\varepsilon)$, và là phi tuyến mạnh khi hệ số phi tuyến có độ lớn hữu hạn không phụ thuộc $\varepsilon$, ví dụ cỡ $O(1)$.
  • Định lý Atalik và Utku (1976): Biểu diễn kỳ vọng tích $E[\mathbf{y} \mathbf{\Phi}^T(\mathbf{y})] = E[\mathbf{y}\mathbf{y}^T] E[\nabla \mathbf{\Phi}^T(\mathbf{y})]$ cho véc tơ Gaussian trung bình không, đóng vai trò bản lề để tuyến tính hóa các lực phi tuyến.
  • Phương trình Fokker-Planck (FP) và thế năng dừng: Dạng phương trình vi phân Parabolic cho hàm mật độ xác suất dừng $f(\mathbf{x})$. Trong trường hợp dòng xác suất $G_i(\mathbf{x})$ triệt tiêu trên toàn không gian $\mathbb{R}^n$ (Stratonovich 1963, Cai và Lin 1988), nghiệm chính xác của hàm mật độ xác suất dừng được xác định qua hàm thế năng dừng $U(\mathbf{x})$: $$f(\mathbf{x}) = \frac{1}{b(\mathbf{x})} \exp\left[ -2 \int_{\mathbf{a}0}^{\mathbf{x}} \frac{1}{b(\mathbf{x})} \sum{i=1}^n A_i(\mathbf{x}) dx_i \right]$$
  • Khung phân tích tiêu chuẩn đối ngẫu: Gồm hai bước thay thế liên tiếp:
    • Bước thuận: Thay thế hàm phi tuyến $\mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}})$ bằng hàm tuyến tính $\mathbf{C}_e \dot{\mathbf{q}} + \mathbf{K}_e \mathbf{q}$, sinh ra véc tơ sai số $\mathbf{e} = \mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}}) - \mathbf{C}_e \dot{\mathbf{q}} - \mathbf{K}_e \mathbf{q}$.
    • Bước đối ngẫu: Thay thế ngược lại tổ hợp tuyến tính $\mathbf{C}_e \dot{\mathbf{q}} + \mathbf{K}_e \mathbf{q}$ bằng một thành phần phi tuyến cùng lớp $\mathbf{D}\mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}})$, sinh ra véc tơ sai số đối ngẫu $\mathbf{\psi} = \mathbf{C}_e \dot{\mathbf{q}} + \mathbf{K}_e \mathbf{q} - \mathbf{D}\mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}})$.
    • Phiếm hàm tối ưu kết hợp: Cực tiểu hóa đồng thời hai sai số với tham số điều chỉnh $\lambda > 0$: $$E[\mathbf{e}^T \mathbf{e}] + \lambda E[\mathbf{\psi}^T \mathbf{\psi}] \to \min_{\mathbf{C}_e, \mathbf{K}_e, \mathbf{D}}$$

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

Luận án kết hợp phương pháp nghiên cứu lý thuyết giải tích và phương pháp mô phỏng số:

  1. Phương pháp giải tích tuyến tính hóa: Thiết lập hệ phương trình ma trận đại số đóng kín bậc cao để xác định đồng thời các ma trận hệ số $\mathbf{C}_e, \mathbf{K}_e, \mathbf{D}$. Kết hợp ma trận hàm đáp ứng tần số $\mathbf{H}(\omega) = [-\omega^2 \mathbf{M} + i\omega(\mathbf{C} + \mathbf{C}_e) + (\mathbf{K} + \mathbf{K}_e)]^{-1}$ và tích phân Duhamel trong miền tần số để xác định ma trận mô men bậc hai của đáp ứng $E[\mathbf{q}\mathbf{q}^T]$ và $E[\dot{\mathbf{q}}\dot{\mathbf{q}}^T]$.
  2. Phương pháp giải tích chính xác: Áp dụng phương trình Fokker-Planck dừng để tìm nghiệm giải tích chính xác của hàm mật độ xác suất và các mô men đáp ứng cho các hệ hai bậc tự do có độ cứng phi tuyến và bài toán dầm liên tục (trường hợp đơn mode, hai mode) làm chuẩn đánh giá sai số.
  3. Phương pháp mô phỏng số Monte-Carlo (MC): Ứng dụng mô phỏng số quá trình ngẫu nhiên đầu vào và tích phân số phương trình chuyển động phi tuyến theo thời gian để tính toán giá trị thống kê bình phương trung bình và ước lượng hàm mật độ xác suất thực nghiệm (pdf), phục vụ đối sánh cho các bài toán không tồn tại nghiệm giải tích chính xác (như hệ có cả cản và độ cứng phi tuyến).

Nguồn dữ liệu nghiên cứu

Dữ liệu nghiên cứu bao gồm các mô hình kết cấu cơ học tiêu chuẩn: hệ 2 bậc tự do phi tuyến bậc ba, bậc năm; hệ 2 bậc tự do có cản phi tuyến; phương trình dao động phi tuyến dầm Euler-Bernoulli liên tục và dầm mang khối lượng tập trung $M_0$ chịu tải trọng ngẫu nhiên dừng.

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

Chương 1. Tổng quan về phương pháp tuyến tính hóa tương đương

Chương 1 thiết lập bức tranh tổng thể về cơ sở lý thuyết và các phương pháp giải dao động ngẫu nhiên phi tuyến:

  • Trình bày lịch sử phát triển của phương pháp tuyến tính hóa tương đương từ các công trình nguyên bản của Booton (1954), Kazakov (1954), Caughey (1956) đến các chuyên khảo của Roberts và Spanos (1990), Socha (2008). Sơ đồ khối tổng quát của phương pháp tuyến tính hóa được mô hình hóa qua hai bước lặp: thiết lập biểu thức hệ số tương đương và giải hệ tuyến tính để cập nhật các đặc trưng đáp ứng cho tới khi hội tụ.
  • Hệ thống hóa các khái niệm giải tích ngẫu nhiên: hàm đặc trưng n chiều $\Phi_{\mathbf{X}}(\mathbf{\theta})$, phân bố chuẩn đa chiều $g_G(\mathbf{x})$, các tính chất mô men bậc chẵn và mô men bậc lẻ triệt tiêu của quá trình Gaussian trung bình không (công thức (1.10)–(1.16)).
  • Trình bày định nghĩa toán học của quá trình Markov, quá trình Wiener và đạo hàm suy rộng tạo nên ồn trắng. Đối với ồn màu, luận án định nghĩa thời gian tương quan $t_c = \frac{1}{R_x(0)}\int_0^\infty |R_x(\tau)| d\tau$ để xác định điều kiện chuyển về xấp xỉ ồn trắng hoặc quá trình dải rộng với mật độ phổ $S_0$.
  • Phân biệt định lượng giữa hệ phi tuyến yếu (tham số phi tuyến nhỏ bậc $O(\varepsilon)$) và hệ phi tuyến mạnh (tham số phi tuyến bậc $O(1)$) trên cơ sở lý thuyết của Babitsky và Krupenin (2001).
  • Trình bày kỹ thuật tuyến tính hóa tương đương thông thường dựa trên tiêu chuẩn sai số bình phương trung bình của Atalik và Utku (1976), dẫn ra công thức hệ số khối lượng, cản và độ cứng tương đương qua đạo hàm riêng của véc tơ phi tuyến: $m_{ij}^e = E[\partial \Phi_i / \partial \ddot{q}j]$, $c{ij}^e = E[\partial \Phi_i / \partial \dot{q}j]$, $k{ij}^e = E[\partial \Phi_i / \partial q_j]$.
  • Khảo sát các tiếp cận theo phiếm hàm năng lượng: tiêu chuẩn bình phương trung bình của sai số thế năng đàn hồi $E[(V(\mathbf{q}) - V_e(\mathbf{q}))^2] \to \min$ (Zhang và cs 1990) và tiêu chuẩn dựa trên hàm hao tán năng lượng $E[s(\dot{\mathbf{x}})(W(\dot{\mathbf{x}}) - W_e(\dot{\mathbf{x}}))^2] \to \min$ với hàm trọng $s(\dot{\mathbf{x}})$ (Elishakoff và Colombi 1993).
  • Phân tích điều kiện khả tích của phương trình Fokker-Planck dừng đa chiều dạng Itô khi véc tơ dòng xác suất $G_i(\mathbf{x})$ triệt tiêu, dẫn ra nghiệm giải tích chính xác của hàm mật độ xác suất theo công thức thế năng dừng (1.20).

Chương 2. Phương pháp tuyến tính hóa tương đương theo tiêu chuẩn đối ngẫu

Chương 2 là chương lý thuyết cốt lõi, mở rộng tiêu chuẩn đối ngẫu cho hệ nhiều bậc tự do:

  • Xét phương trình chuyển động phi tuyến tổng quát: $$\mathbf{M}\ddot{\mathbf{x}} + \mathbf{C}\dot{\mathbf{x}} + \mathbf{K}\mathbf{x} + \mathbf{\Phi}(\mathbf{x}, \dot{\mathbf{x}}) = \mathbf{Q}(t)$$ với véc tơ lực kích động $\mathbf{Q}(t)$ là quá trình Gaussian ồn trắng trung bình không có ma trận mật độ phổ $\mathbf{S}_Q(\omega)$.
  • Hệ tuyến tính hóa tương đương được thiết lập: $$\mathbf{M}\ddot{\mathbf{q}} + (\mathbf{C} + \mathbf{C}_e)\dot{\mathbf{q}} + (\mathbf{K} + \mathbf{K}_e)\mathbf{q} = \mathbf{Q}(t)$$
  • Định nghĩa hai véc tơ sai số: sai số thuận $\mathbf{e} = \mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}}) - \mathbf{C}_e\dot{\mathbf{q}} - \mathbf{K}_e\mathbf{q}$ và sai số đối ngẫu $\mathbf{\psi} = \mathbf{C}e\dot{\mathbf{q}} + \mathbf{K}e\mathbf{q} - \mathbf{D}\mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}})$, trong đó $\mathbf{D} = [d{ij}]{n \times n}$ là ma trận hệ số thay thế trở lại.
  • Thiết lập điều kiện cực trị của phiếm hàm $E[\mathbf{e}^T\mathbf{e}] + \lambda E[\mathbf{\psi}^T\mathbf{\psi}] \to \min$. Tác giả chứng minh rằng bài toán tối ưu trên tách thành $n$ bài toán cực trị độc lập cho từng thành phần thứ $i$ ($i = 1, \dots, n$): $$\frac{\partial \varepsilon_i}{\partial c_{ij}^e} = 0, \quad \frac{\partial \varepsilon_i}{\partial k_{ij}^e} = 0, \quad \frac{\partial \varepsilon_i}{\partial d_{ij}} = 0$$

  • Dẫn ra hệ phương trình đại số ma trận dạng đóng (2.13) kích thước $3n \times 3n$ để giải đồng thời các véc tơ hàng $\mathbf{K}_i^{eT}, \mathbf{C}_i^{eT}, \mathbf{D}_i^T$: $$\begin{bmatrix} (1+\lambda) E[\mathbf{q}\mathbf{q}^T] & (1+\lambda) E[\mathbf{q}\dot{\mathbf{q}}^T] & -\lambda E[\mathbf{q}\mathbf{\Phi}^T] \ (1+\lambda) E[\dot{\mathbf{q}}\mathbf{q}^T] & (1+\lambda) E[\dot{\mathbf{q}}\dot{\mathbf{q}}^T] & -\lambda E[\dot{\mathbf{q}}\mathbf{\Phi}^T] \ -\lambda E[\mathbf{\Phi}\mathbf{q}^T] & -\lambda E[\mathbf{\Phi}\dot{\mathbf{q}}^T] & \lambda E[\mathbf{\Phi}\mathbf{\Phi}^T] \end{bmatrix} \begin{bmatrix} \mathbf{K}_i^{eT} \ \mathbf{C}_i^{eT} \ \mathbf{D}_i^T \end{bmatrix}

    \begin{bmatrix} E[\Phi_i \mathbf{q}] \ E[\Phi_i \dot{\mathbf{q}}] \ \mathbf{0} \end{bmatrix}$$
  • Xây dựng ma trận hàm đáp ứng tần số $\mathbf{H}(\omega) = [-\omega^2 \mathbf{M} + i\omega(\mathbf{C} + \mathbf{C}_e) + (\mathbf{K} + \mathbf{K}e)]^{-1}$ và thiết lập công thức tích phân xác định ma trận hiệp phương sai đáp ứng: $$E[\mathbf{q}\mathbf{q}^T] = \int{-\infty}^\infty \mathbf{H}(\omega) \mathbf{S}Q(\omega) \mathbf{H}^{*T}(\omega) d\omega, \quad E[\dot{\mathbf{q}}\dot{\mathbf{q}}^T] = \int{-\infty}^\infty \omega^2 \mathbf{H}(\omega) \mathbf{S}_Q(\omega) \mathbf{H}^{*T}(\omega) d\omega$$
  • Khảo sát các trường hợp riêng: hệ chỉ có độ cứng phi tuyến, hệ chỉ có cản phi tuyến, và thiết lập công thức giải tích cụ thể cho 3 mô hình hệ hai bậc tự do:
    1. Hệ có độ cứng phi tuyến bậc ba.
    2. Hệ có độ cứng phi tuyến bậc năm.
    3. Mô hình hệ có cản và độ cứng phi tuyến đồng thời.

Chương 3. Phân tích đáp ứng của dầm chịu tải trọng ngẫu nhiên

Chương 3 mở rộng ứng dụng tiêu chuẩn đối ngẫu sang hệ liên tục:

  • Thiết lập phương trình vi phân đạo hàm riêng phi tuyến mô tả dao động uốn của dầm Euler-Bernoulli tựa giản đơn có kể đến hiệu ứng kéo dọc trục do biến dạng lớn và tựa trên nền đàn hồi Winkler hệ số $k_f$: $$\rho A \frac{\partial^2 w}{\partial t^2} + c \frac{\partial w}{\partial t} + E I \frac{\partial^4 w}{\partial x^4} + k_f w - \left[ \frac{E A}{2L} \int_0^L \left(\frac{\partial w}{\partial x}\right)^2 dx \right] \frac{\partial^2 w}{\partial x^2} = q(x, t)$$
  • Sử dụng phương pháp Galerkin, khai triển độ võng theo các hàm dạng riêng $w(x, t) = \sum_{m=1}^M w_m(t) \sin\left(\frac{m\pi x}{L}\right)$ để đưa phương trình liên tục về hệ phương trình vi phân thường phi tuyến bậc ba đối với các biên độ mode $w_m(t)$.
  • Áp dụng tiêu chuẩn đối ngẫu thiết lập hệ số độ cứng tương đương cho từng phương trình mode.
  • Tính toán đáp ứng giải tích cho hai bài toán dầm tựa giản đơn:
    • Trường hợp đơn mode: Xác định nghiệm chính xác từ phương trình Fokker-Planck dừng, so sánh với nghiệm xấp xỉ của ba phương pháp: tuyến tính hóa thông thường (TT), tuyến tính hóa năng lượng (NL) và tiêu chuẩn đối ngẫu (DN).
    • Trường hợp hai mode: Khảo sát hiệu ứng tương tác phi tuyến liên kết giữa mode 1 và mode 2.
  • Mở rộng tiêu chuẩn đối ngẫu cho bài toán dầm mang khối lượng tập trung $M_0$ tại vị trí $d_0$:
    • Bổ sung số hạng quán tính Dirac-delta $M_0 \delta(x - d_0) \frac{\partial^2 w}{\partial t^2}$ vào phương trình vi phân chuyển động.
    • Xây dựng ma trận khối lượng và độ cứng tương đương có tính đến vị trí và độ lớn của tải trọng tập trung.
    • Thiết lập tiêu chuẩn đối ngẫu và tính toán nghiệm xấp xỉ của đáp ứng bình phương trung bình.

Chương 4. Tính toán và mô phỏng số đánh giá độ chính xác của phương pháp

Chương 4 thực hiện các tính toán số nhằm kiểm chứng và đánh giá độ chính xác của phương pháp đối ngẫu đề xuất qua 14 bảng số liệu và 9 đồ thị:

  • Hệ hai bậc tự do có độ cứng phi tuyến bậc ba và bậc năm:
    • Đối sánh đáp ứng bình phương trung bình $E[x_1^2], E[x_2^2]$ giữa nghiệm chính xác (CX), tuyến tính hóa thông thường (TT), tuyến tính hóa năng lượng (NL) và tiêu chuẩn đối ngẫu (DN) khi tham số phi tuyến $\beta_1$ thay đổi (Bảng 4.1 và Bảng 4.2). Kết quả chứng minh khi $\beta_1$ tăng từ nhỏ (phi tuyến yếu) lên lớn (phi tuyến mạnh), sai số của phương pháp đối ngẫu luôn nhỏ hơn rõ rệt so với phương pháp thông thường.
  • Hệ hai bậc tự do có cản và độ cứng phi tuyến đồng thời:
    • So sánh kết quả xấp xỉ của tiêu chuẩn đối ngẫu với nghiệm mô phỏng Monte-Carlo (MC) qua các bảng sai số khi các tham số $\alpha_1, \alpha_2$ thay đổi (Bảng 4.4, 4.5 cho $\alpha_2 = 0$ và Bảng 4.6, 4.7 cho $\alpha_2 = 2$).
    • Dựng đồ thị hàm mật độ xác suất (pdf) của các đáp ứng $x_1, x_2$ trong các trường hợp $\alpha_4 = 0$ và $\alpha_4 = 2$ (Hình 4.1, 4.2, 4.3), khẳng định phân bố xấp xỉ từ phương pháp đối ngẫu bám sát đường cong mô phỏng Monte-Carlo.
  • Nghiệm số cho bài toán dao động dầm liên tục:
    • Bảng 4.8 khảo sát dầm đơn mode với các tham số chuẩn hóa ($\omega_0 = 1, \beta = 1, \gamma = 0.1, \rho A = 1, S_0 = 1$) theo sự biến thiên của đại lượng $1/R$ (nghịch đảo bán kính quán tính, đặc trưng cho độ mảnh của dầm). Kết quả đáp ứng bình phương trung bình $E[w_1^2]$ (đơn vị $\text{m}^2$) được so sánh giữa CX, TT, NL và DN, minh chứng phương pháp DN cho kết quả tiệm cận nghiệm chính xác nhất trên toàn dải $1/R$.
    • Bảng 4.9 đánh giá ảnh hưởng của hệ số độ cứng nền đàn hồi $\tilde{k}_f$ (với $R = 1/200$), cho thấy nền đàn hồi làm giảm biên độ dao động dầm nhưng độ chính xác của DN vẫn vượt trội so với TT và NL.
    • Bảng 4.10 khảo sát trường hợp hai mode dao động theo biến thiên $1/R$.
    • Bảng 4.11 và Hình 4.6 khảo sát trường hợp 5 mode dao động ($M = 5$), định lượng tỉ số đáp ứng bình phương trung bình giữa mode 1 với các mode bậc cao (mode 2, 3, 4, 5) khi $1/R$ thay đổi.
  • Đáp ứng của dầm mang khối lượng tập trung:
    • Bảng 4.12 và các Hình 4.7, 4.8, 4.9 phân tích tỉ số đáp ứng bình phương trung bình phụ thuộc vào độ mảnh dầm $L/h$, tỉ số khối lượng tập trung $M_0/m_0$ (khảo sát tại $M_0/m_0 = 5$ hoặc thay đổi), và vị trí đặt khối lượng tập trung $d_0/L$ (tại $L/h = 200, M_0/m_0 = 5$).
  • Phân tích tham số điều chỉnh $\lambda$:
    • Bảng 4.13 và Bảng 4.14 khảo sát sự thay đổi sai số của $E[x_1^2]$ và $E[x_2^2]$ theo tham số điều chỉnh $\lambda$ khi $\beta_1$ thay đổi, chứng minh tính ổn định của giải thuật và sự phù hợp khi chọn $\lambda = 1$.
Nhóm bài toán kiểm chứng Các thông số khảo sát trong văn bản Đối tượng và phương pháp đối sánh Kết quả rút ra từ các bảng/hình
Hệ 2 bậc tự do độ cứng phi tuyến $\beta_1, \beta_3, \beta_5$; trường hợp bậc ba và bậc năm So sánh DN, TT, NL với nghiệm chính xác CX (Bảng 4.1, 4.2, 4.3) Sai số của DN nhỏ hơn nhiều so với TT khi tham số phi tuyến $\beta_1$ lớn
Hệ 2 bậc tự do có cản phi tuyến $\alpha_1, \alpha_2, \alpha_4$ thay đổi So sánh DN với mô phỏng Monte-Carlo MC (Bảng 4.4 - 4.7, Hình 4.1 - 4.3) DN cho kết quả bám sát MC; hàm mật độ xác suất trùng khít tốt
Dầm tựa giản đơn đơn mode $\omega_0=1, \beta=1, \gamma=0.1, S_0=1$, $1/R$ thay đổi So sánh CX, TT, NL, DN (Bảng 4.8, Hình 4.4) Sai số của DN thấp nhất trong các phương pháp xấp xỉ khi xét biến dạng lớn
Dầm trên nền đàn hồi Nền $\tilde{k}_f$ thay đổi, $R = 1/200$ So sánh giữa các phương pháp xấp xỉ (Bảng 4.9, Hình 4.5) Nền đàn hồi làm giảm đáp ứng; phương pháp DN bảo toàn độ tin cậy cao
Dầm đa mode và khối lượng tập trung $M=5$; độ mảnh $L/h$, $M_0/m_0$, vị trí $d_0/L$ DN so với TT và tuyến tính (Bảng 4.10 - 4.12, Hình 4.6 - 4.9) Đánh giá được tương tác đa mode và ảnh hưởng định lượng của tải tập trung
Khảo sát tham số điều chỉnh $\lambda$ $\lambda$ thay đổi, $\beta_1$ thay đổi Đánh giá độ nhạy sai số của $E[x_1^2], E[x_2^2]$ (Bảng 4.13, 4.14) Việc chọn $\lambda = 1$ mang lại nghiệm cân bằng và ổn định

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

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

  1. Thiết lập hoàn chỉnh tiêu chuẩn đối ngẫu cho hệ nhiều bậc tự do: Luận án đã mở rộng thành công lý thuyết tiêu chuẩn đối ngẫu từ hệ một bậc tự do sang hệ ngẫu nhiên nhiều bậc tự do có tính phi tuyến mạnh (cho cả độ cứng phi tuyến và cản phi tuyến).
  2. Xây dựng hệ phương trình ma trận đóng: Luận án thiết lập được hệ phương trình ma trận đóng kích thước $3n \times 3n$ tổng quát để xác định đồng thời ma trận độ cứng tương đương $\mathbf{K}_e$, ma trận cản tương đương $\mathbf{C}_e$ và ma trận thay thế ngược $\mathbf{D}$, kết hợp với ma trận hàm đáp ứng tần số $\mathbf{H}(\omega)$ trong một quy trình giải lặp tự quán.
  3. Phát triển phương pháp cho hệ liên tục: Mở rộng thành công tiêu chuẩn đối ngẫu cho các bài toán động lực học kết cấu liên tục phi tuyến (dầm Euler-Bernoulli tựa giản đơn có xét hiệu ứng kéo dọc trục, có xét liên kết nền đàn hồi và dầm mang khối lượng tập trung) thông qua phép chiếu Galerkin về hệ các phương trình mode dao động.

Đóng góp mới về mặt thực tiễn và tính toán

  1. Nâng cao độ chính xác dự báo đáp ứng: Các kết quả tính toán số khẳng định rằng trong miền phi tuyến mạnh, phương pháp tiêu chuẩn đối ngẫu cho sai số bình phương trung bình nhỏ hơn đáng kể so với phương pháp tuyến tính hóa tương đương thông thường và phương pháp năng lượng khi so sánh với nghiệm chính xác từ phương trình Fokker-Planck hoặc nghiệm mô phỏng Monte-Carlo.
  2. Cung cấp công thức và bảng số liệu kiểm chuẩn: Luận án cung cấp các công thức giải tích cụ thể cho các hệ hai bậc tự do và hệ dầm, đồng thời xây dựng hệ thống 14 bảng số liệu chi tiết làm cơ sở dữ liệu tham chiếu tin cậy cho các nghiên cứu tiếp theo về cơ học ngẫu nhiên.

Kiến nghị và giải pháp đề xuất

  • Đề xuất áp dụng tham số điều chỉnh chuẩn $\lambda = 1$ trong phiếm hàm đối ngẫu cho các tính toán kỹ thuật chung nhằm bảo đảm sự cân bằng tối ưu giữa bước thay thế thuận và bước thay thế đối ngẫu.
  • Đề xuất ứng dụng quy trình tuyến tính hóa đối ngẫu vào việc tính toán thiết kế và kiểm tra độ an toàn của các kết cấu thanh, dầm chịu tải trọng ngẫu nhiên trong xây dựng và giao thông.

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

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

Dựa trên phạm vi và các nội dung được trình bày trong văn bản:

  • Các tính toán giải tích chi tiết và nghiệm số chủ yếu tập trung vào các hệ tập trung hai bậc tự do và bài toán dầm liên tục đưa về số mode dao động hữu hạn (từ 1 đến 5 mode); chưa mở rộng tính toán số cho các hệ có số bậc tự do rất lớn.
  • Mô hình lực kích động ngẫu nhiên trong luận án mới chỉ xem xét cho kích động ngoài dạng quá trình Gaussian ồn trắng trung bình không ở trạng thái dừng; chưa giải quyết cho các lớp kích động ngẫu nhiên phi Gaussian (non-Gaussian), kích động tham số hay các quá trình ngẫu nhiên không dừng biến đổi theo thời gian.
  • Các hàm phi tuyến khảo sát giới hạn ở các hàm giải tích liên tục (dạng đa thức bậc ba, bậc năm); chưa xét đến các dạng phi tuyến không trơn, phi tuyến gián đoạn hoặc các đặc tính trễ hình học phức tạp (như mô hình trễ Bouc-Wen hay va đập).

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

  • Mở rộng tiêu chuẩn đối ngẫu cho các hệ ngẫu nhiên nhiều bậc tự do có số bậc tự do lớn hơn và các hệ kết cấu không gian phức tạp hơn như tấm, vỏ, khung nhiều tầng.
  • Nghiên cứu tiêu chuẩn đối ngẫu dưới tác dụng của các kích động ngẫu nhiên dải hẹp, ồn màu hoặc kích động không dừng.
  • Nghiên cứu áp dụng tiêu chuẩn đối ngẫu cho các hệ có phi tuyến trễ và các bài toán điều khiển, giảm chấn kết cấu công trình chịu động đất và gió bão.

Giá trị tham khảo

Đối tượng thụ hưởng

  • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Cơ học vật rắn, Cơ học kỹ thuật, Xây dựng công trình, Cơ học công trình và Toán ứng dụng; phục vụ trực tiếp cho việc giảng dạy và nghiên cứu chuyên sâu về dao động ngẫu nhiên phi tuyến.
  • Kỹ sư tính toán kết cấu: Các kỹ sư thiết kế công trình dân dụng, cầu đường, công trình biển cần công cụ tính toán xấp xỉ giải tích nhanh và có độ chính xác cao để dự báo đáp ứng dao động của kết cấu dầm, thanh dưới tác dụng của tải trọng ngẫu nhiên (gió, sóng, động đất).

Các phần có giá trị tham khảo chuyên sâu

  • Chương 2: Toàn bộ phần thiết lập toán học của hệ phương trình ma trận đóng kín (2.13) và giải thuật lặp xác định các ma trận hệ số tuyến tính hóa theo tiêu chuẩn đối ngẫu.
  • Chương 3: Kỹ thuật rời rạc hóa dầm Euler-Bernoulli có xét biến dạng lớn và khối lượng tập trung đưa về các phương trình mode dao động kết hợp tiêu chuẩn đối ngẫu.
  • Chương 4: Tập hợp 14 bảng dữ liệu số chi tiết về đáp ứng bình phương trung bình và sai số tương đối, đóng vai trò là tập dữ liệu chuẩn (benchmark) cho các nghiên cứu phát triển thuật toán xấp xỉ và mô phỏng số.

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

1. Bản chất cốt lõi của "tiêu chuẩn đối ngẫu" trong luận án là gì và khác gì so với phương pháp tuyến tính hóa tương đương thông thường?

Trong phương pháp tuyến tính hóa tương đương thông thường, thành phần phi tuyến $\mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}})$ chỉ được thay thế một chiều bằng thành phần tuyến tính $\mathbf{C}_e\dot{\mathbf{q}} + \mathbf{K}_e\mathbf{q}$ thông qua việc cực tiểu hóa sai số $\mathbf{e}$. Tiêu chuẩn đối ngẫu bổ sung thêm bước thứ hai (bước đối ngẫu): thay thế ngược lại tổ hợp tuyến tính này bằng một thành phần phi tuyến cùng lớp $\mathbf{D}\mathbf{\Phi}(\mathbf{q}, \dot{\mathbf{q}})$ với sai số $\mathbf{\psi}$. Bản chất của tiêu chuẩn đối ngẫu là cực tiểu hóa đồng thời phiếm hàm kết hợp $E[\mathbf{e}^T\mathbf{e}] + \lambda E[\mathbf{\psi}^T\mathbf{\psi}]$, giúp tạo lập sự cân bằng tương hỗ giữa phép thay thế thuận và thay thế nghịch, từ đó khắc phục nhược điểm dự báo sai lệch lớn của phương pháp kinh điển khi hệ có tính phi tuyến mạnh.

2. Tham số điều chỉnh $\lambda$ trong tiêu chuẩn đối ngẫu có vai trò gì và trong luận án được lựa chọn như thế nào?

Tham số điều chỉnh $\lambda$ đóng vai trò là trọng số điều tiết mức độ đóng góp của sai số đối ngẫu so với sai số thuận. Khi $\lambda = 0$, tiêu chuẩn đối ngẫu quay trở về đúng phương pháp tuyến tính hóa kinh điển. Trong các tính toán giải tích và mô phỏng số của luận án, tác giả chọn giá trị chuẩn là $\lambda = 1$ để cân bằng bình đẳng giữa hai bước thay thế. Luận án cũng thực hiện phân tích độ nhạy của $\lambda$ tại Bảng 4.13 và Bảng 4.14, chứng minh rằng kết quả đáp ứng có tính ổn định cao và việc lựa chọn $\lambda = 1$ đem lại sai số nhỏ trên toàn dải biến thiên của tham số phi tuyến.

3. Luận án sử dụng những công cụ đối chuẩn nào để kiểm chứng độ chính xác của phương pháp đề xuất?

Luận án sử dụng hai công cụ đối chuẩn chính:

  1. Nghiệm chính xác (CX): Nhận được từ việc giải giải tích phương trình Fokker-Planck dừng dựa trên lý thuyết thế năng dừng cho các hệ hai bậc tự do có độ cứng phi tuyến và bài toán dầm đơn mode, hai mode.
  2. Nghiệm mô phỏng số Monte-Carlo (MC): Được thực hiện thông qua mô phỏng số các quá trình kích động ngẫu nhiên và giải tích phân số phương trình chuyển động theo thời gian để ước lượng các mô men đáp ứng và hàm mật độ xác suất trong các bài toán không có nghiệm chính xác (như hệ có cả cản và độ cứng phi tuyến).

4. Luận án đã mở rộng tiêu chuẩn đối ngẫu cho bài toán dầm liên tục như thế nào?

Tác giả đã mở rộng tiêu chuẩn đối ngẫu cho dầm Euler-Bernoulli tựa giản đơn có xét đến lực kéo dọc trục do biến dạng lớn và tựa trên nền đàn hồi Winkler $k_f$, cũng như mô hình dầm mang khối lượng tập trung $M_0$ tại vị trí $d_0$. Bằng phương pháp Galerkin, chuyển vị dầm được chiếu lên các hàm dạng riêng $\sin(m\pi x/L)$ để chuyển phương trình đạo hàm riêng liên tục thành hệ phương trình vi phân thường cho các biên độ mode dao động $w_m(t)$, sau đó áp dụng tiêu chuẩn đối ngẫu để xác định các hệ số độ cứng tương đương cho từng mode (đơn mode, hai mode và khảo sát lên đến 5 mode).

5. Sự vượt trội về độ chính xác của tiêu chuẩn đối ngẫu so với các phương pháp khác thể hiện rõ nhất trong điều kiện nào?

Sự vượt trội của tiêu chuẩn đối ngẫu thể hiện rõ nhất khi hệ chuyển từ trạng thái phi tuyến yếu sang trạng thái phi tuyến mạnh (tương ứng với các hệ số phi tuyến $\beta_1, \alpha_1$ lớn hoặc dầm có độ mảnh $1/R$ lớn). Trong miền này, phương pháp thông thường (TT) và phương pháp năng lượng (NL) dẫn đến sai số tăng vọt so với nghiệm chính xác, trong khi phương pháp đối ngẫu (DN) vẫn duy trì sai số ở mức rất thấp và hàm mật độ xác suất xấp xỉ bám sát đường cong mô phỏng Monte-Carlo.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Như Hiếu đã xây dựng thành công cơ sở lý thuyết và thuật toán hoàn chỉnh cho tiêu chuẩn đối ngẫu trong phương pháp tuyến tính hóa tương đương áp dụng cho hệ ngẫu nhiên nhiều bậc tự do và hệ liên tục chịu kích động ngẫu nhiên. Công trình đã giải quyết thành công bài toán đóng kín hệ phương trình xác định đồng thời các ma trận hệ số tuyến tính hóa, đồng thời chứng minh tính hiệu quả vượt trội trong việc giảm thiểu sai số đối với các hệ có tính phi tuyến mạnh qua so sánh đối chứng chặt chẽ với nghiệm chính xác của phương trình Fokker-Planck và nghiệm mô phỏng Monte-Carlo. Kết quả nghiên cứu cung cấp một công cụ giải tích xấp xỉ có độ tin cậy cao, đóng góp thiết thực cho lý thuyết động lực học ngẫu nhiên và thực tiễn tính toán kết cấu công trình chịu tải trọng ngẫu nhiên.