Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Phân tích tĩnh và động kết cấu vỏ thoải FGM hai độ cong trong môi trường nhiệt" tập trung giải quyết các bài toán cơ học kết cấu đối với vật liệu có cơ tính biến thiên (Functionally Graded Material – FGM). Vật liệu FGM là hệ vật liệu composite được các nhà khoa học Nhật Bản phát triển từ năm 1984 tại Phòng thí nghiệm Hàng không Quốc gia Nhật Bản, thường gồm hai thành phần chính là gốm và kim loại với tỉ lệ biến thiên liên tục theo không gian. Cấu trúc biến đổi trơn này giúp kết hợp tính năng kháng nhiệt, mô đun đàn hồi cao của gốm với độ bền dẻo của kim loại, đồng thời khắc phục các nhược điểm thường gặp ở vật liệu composite phân lớp truyền thống như hiện tượng tách lớp và tập trung ứng suất tiếp xúc.

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

Trong các lĩnh vực kỹ thuật như hàng không vũ trụ, hạt nhân và công nghiệp ô tô, kết cấu vỏ thường xuyên làm việc trong môi trường nhiệt độ cao và chịu tác dụng đồng thời của tải trọng cơ học lẫn trường nhiệt. Việc mô phỏng số và tính toán chính xác ứng xử tĩnh và động của kết cấu vỏ FGM là tiền đề kỹ thuật phục vụ công tác thiết kế.

Tổng quan các nghiên cứu cho thấy phần lớn các công trình trước đây tập trung vào kết cấu tấm phẳng hoặc các loại vỏ kín (vỏ trụ, vỏ nón). Đối với các dạng vỏ hở như vỏ thoải hai độ cong, panel trụ, panel cầu, các nghiên cứu đa số sử dụng lời giải giải tích phi tuyến và bị giới hạn ở các điều kiện biên đơn giản cùng dạng hình học cố định. Khi phân tích bài toán động trong môi trường nhiệt, các công bố thường chia thành hai hướng tiếp cận:

  • Hướng tiếp cận thứ nhất: Nhiệt độ chỉ làm thay đổi cơ tính vật liệu cấu thành.
  • Hướng tiếp cận thứ hai: Nhiệt độ vừa làm thay đổi cơ tính vật liệu, vừa gây biến dạng tĩnh và hình thành trạng thái ứng suất ban đầu (ứng suất trước) trước khi kết cấu dao động.

Tuy nhiên, trong hướng tiếp cận thứ hai, phần lớn các tác giả áp dụng công thức giả thiết ứng suất nhiệt phân bố đều trên mỗi lớp và không phụ thuộc vào điều kiện biên của kết cấu. Rất ít nghiên cứu giải quyết bài toán tĩnh dưới tải trọng nhiệt để xác định chính xác trường ứng suất ban đầu đưa vào phương trình dao động động lực học, cũng như mở rộng cho các kết cấu vỏ có dạng hình học phức tạp.

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

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

  1. Xây dựng các hệ thức quan hệ, phương trình chủ đạo và thiết lập nghiệm giải tích của bài toán phân tích tĩnh và động của kết cấu vỏ thoải FGM hai độ cong trong môi trường nhiệt độ.
  2. Xây dựng mô hình, thuật toán Phần tử hữu hạn, sử dụng phần tử vỏ 3D suy biến phân tích tĩnh và động của vỏ thoải FGM có xét đến ảnh hưởng của nhiệt độ.
  3. Viết chương trình tính trên nền Matlab để khảo sát ảnh hưởng của các tham số vật liệu, kích thước hình học, điều kiện biên và nhiệt độ đến độ võng, ứng suất và ứng xử động của vỏ thoải FGM hai độ cong.

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

  • Đối tượng nghiên cứu: Kết cấu vỏ thoải FGM hai độ cong có hình chiếu bằng là hình chữ nhật làm việc trong môi trường nhiệt độ.
  • Phạm vi nghiên cứu: Tính toán độ võng, ứng suất, tần số dao động riêng và đáp ứng động của vỏ thoải FGM hai độ cong theo lý thuyết biến dạng cắt bậc nhất (FSDT) khi kết cấu vỏ chịu tải trọng cơ học và tải trọng nhiệt độ ở trạng thái dừng.

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

Luận án đã tổng hợp hệ thống tài liệu nghiên cứu trong và ngoài nước liên quan đến ứng xử cơ nhiệt của kết cấu tấm và vỏ FGM:

Các hướng nghiên cứu quốc tế

  • Phân tích ứng suất cơ nhiệt: Reddy và Cheng (2001) sử dụng phương pháp tiệm cận phân tích tấm FGM và tấm FGM gắn áp điện. Wang và cộng sự (2004) thiết lập lời giải tiệm cận cho tấm mỏng P-FGM chịu sốc nhiệt theo lý thuyết Lord và Shulman. Vel và Batra (2004) sử dụng chuỗi lũy thừa nghiên cứu biến dạng cơ nhiệt tấm FGM bốn biên tựa khớp, so sánh lý thuyết tấm cổ điển (CLPT), biến dạng cắt bậc nhất (FSDT) và bậc ba (TSDT). Reddy và Chin (1998) giải phương trình truyền nhiệt và nhiệt đàn hồi cho ống trụ tròn FGM bằng phần tử hữu hạn. Shao và Wang (2006) phân tích đàn hồi nhiệt 3D cho vỏ trụ FGM chịu tải nhiệt dừng. Ootao và Tanigawa (1999, 2007) khảo sát đáp ứng nhiệt đàn hồi theo thời gian của ống trụ và panel trụ FGM. Brischetto và cộng sự (2008), Zenkour (2006, 2007) sử dụng lý thuyết biến dạng cắt dạng sin (SSDT) và lý thuyết bậc ba cho tấm FGM và sandwich. Tounsi và cộng sự (2013), Swaminathan và Sangeetha (2017), Woo và Meguid (2001), Moosaie và Panahi-Kalus (2014), Matsunaga (2009), Bhandari và Purohit (2014), Kiani và cộng sự (2012) khảo sát các trường nhiệt đều, tuyến tính và phi tuyến trên tấm, vỏ cầu và vỏ hai độ cong.
  • Phân tích dao động trong môi trường nhiệt: Yang và Shen (2001, 2003) nghiên cứu dao động riêng và đáp ứng thời gian của tấm FGM có ứng suất trước trên nền đàn hồi. Kim (2005) dùng lý thuyết TSDT và phương pháp Rayleigh-Ritz khảo sát đặc trưng tần số theo hai quy luật truyền nhiệt phi tuyến. Benachour và cộng sự (2011), Chakraverty và Pradhan (2014), Fazzolari (2014), Sobhy (2013), Li và Lu (2009), Attia và cộng sự (2015), Liew và cộng sự (2003, 2004), Park và Kim (2006), Malekzadeh và cộng sự (2009, 2010, 2011), Jooybar và cộng sự (2016), Pradyumna và Bandyopadhyay (2008), Kandasamy và cộng sự (2014), Sheng và Wang (2008), Kadoli và Ganesan (2006) phân tích dao động và ổn định nhiệt của vỏ nón, vỏ trụ và vỏ hai độ cong.

Các hướng nghiên cứu tại Việt Nam

Tại Việt Nam, các nghiên cứu chuyên sâu về kết cấu composite FGM được phát triển mạnh mẽ qua các nhóm nghiên cứu:

  • Hoàng Văn Tùng (2010) với luận án tiến sĩ về ổn định nhiệt đàn hồi của tấm, vỏ FGM; tiếp tục cùng Nguyễn Đình Đức (2010, 2011) công bố về ổn định phi tuyến của tấm, panel trụ FGM chịu tải cơ nhiệt.
  • Nhóm Nguyễn Đình Đức và cộng sự (2012–2016) công bố hàng loạt nghiên cứu về ứng xử phi tuyến, sau ổn định của tấm và vỏ FGM có gân gia cường trên nền đàn hồi trong môi trường nhiệt.
  • Nhóm Đào Huy Bích và cộng sự (2011–2016) nghiên cứu ứng xử phi tuyến tĩnh, ổn định động lực của vỏ cầu thoải FGM chịu áp lực và nhiệt độ.
  • Nhóm Đào Văn Dũng và cộng sự (2013–2016) giải quyết các bài toán phi tuyến giải tích và bán giải tích cho vỏ FGM có gân gia cường chịu nhiệt.
  • Nhóm Hoàng Xuân Lượng, Nguyễn Thái Chung, Phạm Tiến Đạt (Học viện Kỹ thuật Quân sự) nghiên cứu ứng xử tĩnh và động của kết cấu FGM bằng phương pháp phần tử hữu hạn.

Vị trí và khoảng trống luận án giải quyết

Luận án định vị trọng tâm vào việc thiết lập lời giải giải tích theo cách tiếp cận thứ hai: giải bài toán tĩnh chịu tải nhiệt để xác định chính xác trường ứng suất trước phụ thuộc vào điều kiện biên, sau đó tích hợp vào phương trình dao động động lực học. Đồng thời, luận án phát triển mô hình phần tử hữu hạn sử dụng phần tử vỏ 3D suy biến nhằm khắc phục giới hạn hình học và điều kiện biên của phương pháp giải tích.

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

Mô hình vật liệu và quy luật truyền nhiệt

Luận án nghiên cứu loại vật liệu FGM biến thiên theo quy luật hàm lũy thừa (P-FGM), gồm hai pha gốm (ceramic - $c$) và kim loại (metal - $m$). Tỉ lệ thể tích thành phần gốm $V_c$ biến thiên theo chiều dày $z \in [-h/2, h/2]$ với chỉ số thể tích $p$: $$V_c(z) = \left(\frac{z}{h} + \frac{1}{2}\right)^p, \quad V_m(z) + V_c(z) = 1$$

Các đặc trưng cơ nhiệt hiệu dụng $P(z, T)$ (mô đun Young $E$, hệ số giãn nở nhiệt $\alpha$, hệ số dẫn nhiệt $\kappa$) phụ thuộc tọa độ chiều dày và nhiệt độ theo mô hình Voigt và thực nghiệm Touloukian: $$P(z, T) = P_m(T) + [P_c(T) - P_m(T)] V_c(z)$$ $$P(T) = P_0 (1 + P_{-1}T^{-1} + P_1 T + P_2 T^2 + P_3 T^3)$$ Hệ số Poisson $\nu$ được coi là hằng số ($\nu = \text{const}$).

Trường nhiệt độ $T(z)$ theo chiều dày được khảo sát qua 03 quy luật ở trạng thái dừng:

  1. Nhiệt độ tăng đều: $T(z) = T_m + \Delta T$, với $\Delta T = T_c - T_m$.
  2. Truyền nhiệt tuyến tính: $T(z) = T_m + (T_c - T_m)\left(\frac{z}{h} + \frac{1}{2}\right)$.
  3. Truyền nhiệt phi tuyến: Nghiệm của phương trình truyền nhiệt dừng một chiều $\frac{d}{dz}\left[\kappa(z)\frac{dT}{dz}\right] = 0$: $$T(z) = T_m + (T_c - T_m) \cdot \frac{\eta(z)}{C}$$ với $\eta(z)$ và hằng số $C$ được triển khai dưới dạng chuỗi đa thức lũy thừa phụ thuộc vào tỷ số dẫn nhiệt $\kappa_{cm} = \kappa_c - \kappa_m$ và chỉ số $p$.
Vật liệu thành phần Mô đun Young $E$ (GPa) Hệ số truyền nhiệt $\kappa$ (W/mK) Hệ số giãn nở nhiệt $\alpha$ ($10^{-6}$/K) Khối lượng riêng $\rho$ ($\text{Mg/m}^3$)
Gốm (Zirconia) 100 – 250 2.09 10.0 5.70
Kim loại (Lead alloys/Al) Biến thiên nhiệt độ Biến thiên nhiệt độ Biến thiên nhiệt độ 2.70 – 11.3

Khung lý thuyết biến dạng và phương trình chủ đạo

  • Lý thuyết biến dạng cắt bậc nhất (FSDT - Reissner-Mindlin): Trường chuyển vị của vỏ thoải hai độ cong ($R_x, R_y$ là bán kính cong theo phương $x, y$): $$u(x, y, z, t) = u_0(x, y, t) + z \psi_x(x, y, t)$$ $$v(x, y, z, t) = v_0(x, y, t) + z \psi_y(x, y, t)$$ $$w(x, y, z, t) = w_0(x, y, t)$$ trong đó $u_0, v_0, w_0$ là chuyển vị mặt trung bình; $\psi_x, \psi_y$ là các góc xoay của pháp tuyến; hệ số hiệu chỉnh cắt $k_c = 5/6$.
  • Nguyên lý biến phân:
    • Áp dụng nguyên lý thế năng toàn phần cực tiểu $\delta U_t - \delta W_t = 0$ để thiết lập 05 phương trình vi phân cân bằng bài toán tĩnh.
    • Áp dụng nguyên lý Hamilton $\int_{0}^{T} (\delta \mathcal{T}d - \delta U_d - \delta U_0 + \delta W_d) dt = 0$ để thiết lập hệ phương trình chuyển động bài toán động, trong đó $\delta U_0$ là thế năng biến dạng do ứng suất trước $\sigma{xx}^0, \sigma_{yy}^0, \sigma_{xy}^0$ từ bài toán tĩnh gây ra.

Phương pháp tính toán và công cụ thực hiện

  1. Phương pháp Giải tích:
    • Dùng chuỗi lượng giác kép Navier/Galerkin biểu diễn trường chuyển vị và tải trọng nhằm thỏa mãn các điều kiện biên: Bốn biên tựa khớp (SSSS), Khớp - Ngàm - Khớp - Ngàm (SCSC), Bốn biên ngàm cứng (CCCC).
    • Sử dụng phương pháp Galerkin đưa hệ phương trình vi phân đạo hàm riêng về hệ phương trình đại số tuyến tính $[K]{\Delta^t} = {F_{ch}^t} + {F_{nd}}$ cho bài toán tĩnh và phương trình trị riêng $\det([K] + [K_{ini}] - \omega^2 [M]) = 0$ cho bài toán dao động tự do.
    • Sử dụng phương pháp tích phân từng bước Newmark-$\beta$ kết hợp mô hình cản Rayleigh $[C] = a_1[M] + a_2[K]$ để giải bài toán dao động cưỡng bức.
  2. Phương pháp Phần tử hữu hạn (PTHH):
    • Xây dựng phần tử vỏ 3D suy biến (3D-Degenerated shell element) 8 nút phát triển từ phần tử khối theo Ahmad.
    • Thiết lập 04 hệ tọa độ: Hệ tọa độ tổng thể $xyz$, hệ tọa độ địa phương phần tử $x'y'z'$, hệ tọa độ nút $\vec{V}_1\vec{V}_2\vec{V}_3$ (với véc tơ $\vec{V}_3$ chỉ phương theo chiều dày) và hệ tọa độ tự nhiên $(\xi, \eta, \zeta)$.
  3. Chương trình tính toán: Viết 03 chương trình tính độc lập trên nền Matlab gồm: ShellpanelStatic(GT), ShellpanelVibration(GT) và ShellpanelForcedvibration(GT).

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

Chương 1: Tổng quan các vấn đề liên quan đến nội dung đề tài luận án

Chương 1 trình bày tổng quan cấu trúc, cơ tính và các lĩnh vực ứng dụng của vật liệu FGM trong công nghiệp hàng không vũ trụ (lớp cách nhiệt buồng đốt, mũi tên lửa, ống xả khí), công nghiệp điện hạt nhân (vách lò phản ứng), công nghiệp ô tô (xilanh, má phanh, bánh đà) và y sinh (polyme sinh học trong phẫu thuật, nha khoa).

Chương này hệ thống hóa các mô hình cơ tính biến thiên (Voigt, Mori-Tanaka, Hashin-Shtrikman) và quy luật phân bố thể tích (dạng hàm lũy thừa P-FGM, hàm mũ E-FGM, hàm sigmoid S-FGM). Đồng thời, chương 1 tổng kết các mô hình phân bố nhiệt độ (truyền nhiệt 3D, 1D phi tuyến, 1D tuyến tính, truyền nhiệt đều và dạng đa thức layerwise), kèm theo đánh giá chi tiết các công trình nghiên cứu cơ học tính toán trong và ngoài nước về ứng suất cơ nhiệt, dao động tự do và đáp ứng động của tấm/vỏ FGM.

Chương 2: Phân tích tĩnh và động vỏ thoải FGM hai độ cong trong môi trường nhiệt theo lý thuyết biến dạng cắt bậc nhất bằng tiếp cận giải tích

Chương 2 tập trung xây dựng mô hình toán học giải tích cho vỏ thoải FGM hai độ cong có kích thước chiếu bằng $a \times b$, bán kính cong $R_x, R_y$:

Thành phần bài toán Mô tả toán học và quan hệ giải tích
Quan hệ ứng lực - biến dạng $\begin{bmatrix} {N} \ {M} \ {Q} \end{bmatrix} = \begin{bmatrix} [A] & [B] & [0] \ [B] & [D] & [0] \ [0] & [0] & [A_s] \end{bmatrix} \begin{bmatrix} {\epsilon^0} \ {\kappa} \ {\gamma} \end{bmatrix} - \begin{bmatrix} {N^{nd}} \ {M^{nd}} \ {0} \end{bmatrix}$
Độ cứng tích phân chiều dày $(A_{ij}, B_{ij}, D_{ij}) = \int_{-h/2}^{h/2} Q_{ij}(z, T) (1, z, z^2) dz; \quad A_{44, 55} = k_c \int_{-h/2}^{h/2} Q_{44, 55}(z, T) dz$
Hệ phương trình cân bằng tĩnh Thiết lập 05 phương trình vi phân liên hệ chuyển vị $u_0^t, v_0^t, w_0^t, \psi_x^t, \psi_y^t$ và lực nhiệt $N^{nd}, M^{nd}$.
Lời giải giải tích Galerkin tĩnh $[K]{\Delta_{mn}^t} = {F_{ch}^t} + {F_{nd}} \rightarrow$ Xác định chuyển vị tĩnh và ứng suất ban đầu $\sigma_{xx}^0, \sigma_{yy}^0, \sigma_{xy}^0$.
Ma trận độ cứng hình học ban đầu $[K_{ini}]$ được tính từ tích phân các thành phần lực màng nhiệt $N_{xx}^0, N_{yy}^0, N_{xy}^0$ thu được từ bài toán tĩnh.
Phương trình dao động tự do $([K] + [K_{ini}] - \omega^2 [M]){\Delta_{mn}^0} = {0} \rightarrow$ Giải phương trình $\det([K] + [K_{ini}] - \omega^2 [M]) = 0$.
Phương trình dao động cưỡng bức $[M]{\ddot{\Delta}{mn}^d} + [C]{\dot{\Delta}{mn}^d} + ([K] + [K_{ini}]){\Delta_{mn}^d} = {F_{ch}^d(t)}$ giải bằng thuật toán Newmark-$\beta$.

Chương này hoàn thiện 03 lưu đồ thuật toán và chương trình tính toán tương ứng trên Matlab cho các bài toán phân tích tĩnh, xác định tần số dao động riêng và đáp ứng chuyển vị thời gian.

Chương 3: Phân tích tĩnh và động vỏ FGM trong môi trường nhiệt bằng phương pháp phần tử hữu hạn

Chương 3 mở rộng phạm vi giải quyết cho kết cấu vỏ FGM có bề mặt mô tả bởi hàm toán học tổng quát $z = f(x, y)$ và các điều kiện biên phức tạp bằng phương pháp số:

  • Lựa chọn phần tử vỏ 3D suy biến: Đánh giá ưu nhược điểm so sánh giữa phần tử phẳng, phần tử vỏ cong, phần tử khối và phần tử 3D suy biến. Phần tử 3D suy biến 8 nút được chọn nhờ khả năng mô tả chính xác hình học vỏ thực mà không làm tăng bậc tự do như phần tử khối 3D.
  • Xây dựng hệ phương trình phần tử: Thiết lập quan hệ hình học qua các hàm dạng đẳng tham số $N_i(\xi, \eta)$, véc tơ nút $\vec{V}_{3i}$ theo phương pháp tuyến bề mặt; thiết lập ma trận độ cứng phần tử $[K^e]$, ma trận độ cứng hình học phần tử do nhiệt $[K_g^e]$, ma trận khối lượng phần tử $[M^e]$ và véc tơ lực nút phần tử ${F^e}$.
  • Thuật toán ghép nối và giải hệ: Chuyển đổi giữa các hệ tọa độ (nút $\rightarrow$ địa phương $\rightarrow$ tổng thể), ghép nối ma trận toàn cục và áp dụng các thuật toán giải bài toán tĩnh, bài toán trị riêng tìm tần số dao động riêng và tích phân Newmark cho bài toán động.

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

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

  1. Thiết lập lời giải giải tích hoàn chỉnh theo tiếp cận bài toán có ứng suất trước: Xây dựng thành công hệ phương trình và nghiệm giải tích phân tích tĩnh, dao động tự do và dao động cưỡng bức của vỏ thoải FGM hai độ cong chịu tác dụng đồng thời của tải trọng cơ học và nhiệt độ. Lực màng và ứng suất ban đầu do nhiệt độ được xác định trực tiếp từ lời giải bài toán tĩnh, bảo đảm tính phụ thuộc vào điều kiện biên thực tế của kết cấu chứ không cố định giá trị trên từng lớp.
  2. Ứng dụng hiệu quả phần tử 3D suy biến cho kết cấu vỏ FGM: Xây dựng mô hình và thuật toán phần tử hữu hạn sử dụng phần tử vỏ 3D suy biến cho phép mô phỏng vỏ FGM có hình học mặt cong bất kỳ mô tả bởi hàm toán học $z = f(x, y)$, mở rộng phạm vi nghiên cứu vượt qua các giới hạn của nghiệm giải tích.
  3. Phát triển bộ công cụ tính toán số: Xây dựng bộ chương trình tính toán hoàn chỉnh trên nền phần mềm Matlab để khảo sát định lượng ảnh hưởng của các tham số:
    • Chỉ số thể tích vật liệu $p$.
    • Kích thước hình học ($a/b, a/h, R_x/a, R_y/b$).
    • Các dạng điều kiện biên (SSSS, SCSC, CCCC).
    • Ba quy luật truyền nhiệt theo chiều dày (nhiệt tăng đều, tuyến tính, phi tuyến).
    • Tỉ lệ cản Rayleigh và tỉ số tần số kích thích $\Omega/\omega$.
    • Các đại lượng đầu ra gồm: độ võng, các thành phần ứng suất màng/uốn/cắt, tần số dao động riêng cơ bản và đáp ứng chuyển vị theo thời gian.

Đề xuất và kiến nghị kỹ thuật

  • Kết quả tính toán cung cấp căn cứ định lượng để kỹ sư thiết kế lựa chọn cấu hình vật liệu (tỉ lệ gốm/kim loại), độ dày và bán kính cong phù hợp nhằm tránh hiện tượng cộng hưởng và giảm thiểu ứng suất nhiệt trong kết cấu làm việc ở điều kiện nhiệt độ cao.

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

Các giới hạn của đề tài

  1. Giới hạn mô hình giải tích: Nghiệm giải tích trong Chương 2 chỉ áp dụng được cho các kết cấu vỏ có độ cong không đổi ($R_x, R_y = \text{const}$) như vỏ trụ, vỏ cầu, vỏ yên ngựa và bị giới hạn ở một số điều kiện biên đối xứng thông dụng (SSSS, SCSC, CCCC).
  2. Giả thiết nhiệt đàn hồi: Chưa xét đến quá trình truyền nhiệt không dừng (nhiệt độ phụ thuộc thời gian $\partial T/\partial t \neq 0$); tải trọng nhiệt được quy về bài toán tĩnh dừng.
  3. Mô hình biến dạng: Sử dụng lý thuyết biến dạng cắt bậc nhất FSDT với hệ số hiệu chỉnh cắt $k_c = 5/6$ và giả thiết hệ số Poisson không đổi ($\nu = \text{const}$), chưa áp dụng các lý thuyết biến dạng cắt bậc cao (HSDT) không cần hệ số hiệu chỉnh.
  4. Phạm vi ứng xử: Giới hạn trong bài toán cơ học tuyến tính hình học và vật liệu đàn hồi tuyến tính.

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

  • Mở rộng thuật toán PTHH và giải tích cho bài toán phi tuyến hình học (biến dạng lớn) và phi tuyến vật liệu của vỏ FGM trong môi trường nhiệt.
  • Nghiên cứu bài toán nhiệt động lực học có xét đến trường nhiệt không dừng và sốc nhiệt theo thời gian.
  • Khảo sát kết cấu vỏ FGM có gân gia cường hoặc có gắn các lớp vật liệu thông minh (áp điện, từ giảo).

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn cao cho các nhóm đối tượng:

  • Nghiên cứu sinh và học viên cao học: Tham khảo quy trình thiết lập phương trình vi phân biến phân, giải bài toán cơ học kết cấu bằng phương pháp Galerkin và phương pháp Newmark-$\beta$.
  • Giảng viên và nhà nghiên cứu Cơ học/Xây dựng/Hàng không: Cung cấp phương pháp luận giải quyết bài toán động lực học có xét đến ứng suất trước do trường nhiệt gây ra; khung xây dựng phần tử 3D suy biến cho kết cấu vỏ mặt cong toán học.
  • Kỹ sư tính toán thiết kế kết cấu: Ứng dụng các thuật toán và dữ liệu khảo sát tham số trên nền Matlab để lựa chọn chỉ số vật liệu $p$, bán kính cong $R$ và chiều dày $h$ tối ưu khi chế tạo các tấm vỏ chịu nhiệt trong công nghiệp.

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

1. Luận án lựa chọn cách tiếp cận nào để xét ảnh hưởng của nhiệt độ đến bài toán động?

Luận án lựa chọn cách tiếp cận thứ hai: Nhiệt độ vừa làm thay đổi các đặc trưng cơ tính đàn hồi của vật liệu theo mô hình thực nghiệm Touloukian, vừa làm kết cấu biến dạng tĩnh và xuất hiện trạng thái ứng suất trước ($\sigma_{xx}^0, \sigma_{yy}^0, \sigma_{xy}^0$). Trạng thái ứng suất ban đầu này được xác định từ nghiệm của bài toán tĩnh dưới tác dụng của tải trọng nhiệt, từ đó thiết lập ma trận độ cứng ban đầu $[K_{ini}]$ trong phương trình dao động.

2. Sự khác biệt giữa mô hình giải tích và mô hình phần tử hữu hạn trong luận án là gì?

Mô hình giải tích (Chương 2) sử dụng chuỗi lượng giác kép Navier/Galerkin giải các phương trình vi phân FSDT, cho nghiệm chính xác nhưng chỉ áp dụng được cho vỏ có bán kính cong không đổi ($R_x, R_y = \text{const}$) và các biên thông dụng (SSSS, SCSC, CCCC). Mô hình phần tử hữu hạn (Chương 3) sử dụng phần tử vỏ 3D suy biến 8 nút, cho phép tính toán đối với vỏ có bề mặt cong bất kỳ mô tả bởi hàm $z = f(x, y)$ và điều kiện biên tùy ý.

3. Những quy luật truyền nhiệt nào theo chiều dày vỏ được khảo sát trong công trình?

Luận án khảo sát 03 quy luật truyền nhiệt dừng theo phương chiều dày: (1) Nhiệt độ tăng đều trên toàn chiều dày; (2) Truyền nhiệt tuyến tính từ mặt dưới $T_m$ lên mặt trên $T_c$; (3) Truyền nhiệt phi tuyến nhận được từ việc giải phương trình vi phân truyền nhiệt Fourier một chiều $\frac{d}{dz}[\kappa(z)\frac{dT}{dz}] = 0$.

4. Hệ số hiệu chỉnh cắt và hệ số Poisson được xử lý như thế nào trong các biểu thức tính toán?

Luận án áp dụng lý thuyết biến dạng cắt bậc nhất (FSDT) với hệ số hiệu chỉnh cắt $k_c = 5/6$. Do hệ số Poisson của gốm ($\nu_c$) và kim loại ($\nu_m$) có giá trị xấp xỉ nhau và sự biến đổi theo chiều dày gây ra sai số rất nhỏ, luận án giả thiết hệ số Poisson là hằng số ($\nu = \text{const}$) trong toàn bộ các công thức tính toán.

5. Bộ chương trình tính toán số của luận án được xây dựng trên nền tảng nào và gồm những module nào?

Bộ chương trình được viết trên nền tảng phần mềm Matlab, bao gồm 03 module giải tích chính: ShellpanelStatic(GT) (tính toán chuyển vị tĩnh và ứng suất), ShellpanelVibration(GT) (tính toán tần số và dạng dao động riêng) và ShellpanelForcedvibration(GT) (phân tích đáp ứng động lực học theo thời gian bằng phương pháp Newmark-$\beta$).

Kết luận

Luận án đã giải quyết hoàn chỉnh bài toán phân tích tĩnh, dao động tự do và dao động cưỡng bức của kết cấu vỏ thoải FGM hai độ cong trong môi trường nhiệt độ bằng cả hai phương pháp giải tích và phần tử hữu hạn. Công trình đã luận chứng rõ ràng việc đưa ứng suất trước do nhiệt độ từ bài toán tĩnh vào phương trình chuyển động động lực học, giúp phản ánh chính xác sự phụ thuộc của ứng xử động vào điều kiện biên kết cấu. Các kết quả khảo sát số cùng bộ chương trình tính toán trên Matlab tạo thành nguồn tư liệu khoa học hữu ích phục vụ công tác nghiên cứu, đào tạo và tính toán thiết kế kết cấu vật liệu biến tính trong thực tế kỹ thuật.