Tổng quan luận án

Chuyên đề nghiên cứu "Một số vấn đề của lý thuyết kết cấu vỏ mỏng đàn hồi" (Chuyên đề tiến sĩ số 2, khối lượng 3 đơn vị học trình) được thực hiện tại Bộ môn Sức bền vật liệu - Khoa Xây dựng, Trường Đại học Xây dựng Hà Nội vào tháng 6 năm 2003, thuộc chuyên ngành Sức bền vật liệu và Cơ học kết cấu (Mã số: 2-01-02), dưới sự hướng dẫn khoa học của TS. Nguyễn Văn Hùng và GS. Nguyễn Văn Phó.

Tính cấp thiết của đề tài xuất phát từ nhu cầu thực tiễn trong thiết kế và xây dựng các công trình dân dụng, công nghiệp, công trình chứa nhiên liệu, hàng hải và hàng không. Kết cấu vỏ mỏng không gian có khả năng chịu lực lớn nhờ tận dụng tối đa tính chất cơ lý của vật liệu và tạo thành hệ thống chịu lực không gian ba chiều có khả năng truyền tải trọng liên hoàn, mở rộng khẩu độ kết cấu. Tuy nhiên, việc tính toán vỏ đòi hỏi giải quyết các hệ phương trình đạo hàm riêng phức tạp với hệ số biến thiên. Khoảng trống nghiên cứu được xác định là nhu cầu hệ thống hóa các phương pháp giải giải tích (lý thuyết phi mômen, lý thuyết mômen, lý thuyết hiệu ứng biên, phương pháp tương tự tĩnh học - hình học) và kết nối với các công cụ giải số hiện đại trên máy tính điện tử.

Mục tiêu và nhiệm vụ nghiên cứu của chuyên đề bao gồm:

  1. Hệ thống hóa cơ sở hình học vi phân của mặt vỏ, các giả thiết cơ bản và thiết lập hệ phương trình vi phân chủ đạo của lý thuyết vỏ mỏng đàn hồi.
  2. Trình bày các hướng tiếp cận giải bài toán đàn hồi: giải theo nội lực, giải theo chuyển vị, phương pháp tương tự tĩnh học - hình học và phương pháp phần tử hữu hạn.
  3. Nghiên cứu sâu lý thuyết phi mômen và lý thuyết mômen (kể đến hiệu ứng biên) đối với vỏ tròn xoay chịu tải trọng đối xứng trục.
  4. Ứng dụng giải quyết cụ thể hai bài toán kết cấu thực tế: bể chứa nước ngầm hình trụ bằng bê tông cốt thép đặt trong đất và bể chứa dầu mỏ (kết cấu liên hợp vỏ cầu - vành đai - vỏ trụ - bản đáy), đồng thời kiểm chứng bằng phần mềm SAP2000 Nonlinear.

Đối tượng nghiên cứu là kết cấu vỏ mỏng đàn hồi tuyến tính ($h/R \le 0,05$), tập trung vào các dạng hình học cơ bản trong công trình kỹ thuật như vỏ trụ tròn, vỏ cầu (chỏm cầu), bản tròn phẳng và vành liên kết. Phạm vi nghiên cứu giới hạn ở bài toán tĩnh đối xứng trục dưới tác dụng của trọng lượng bản thân, tải trọng đất đắp, áp lực thủy tĩnh của nước và áp lực dầu mỏ.

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

Lý thuyết vỏ mỏng đàn hồi hình thành từ thế kỷ 19, bắt nguồn từ các định đề hình học và vật lý của lý thuyết uốn bản do G. đề xuất. Đến đầu thế kỷ 20, Staerman công bố những công trình đầu tiên về vỏ bán cầu bằng cách tích phân phương trình vi phân cân bằng, đưa việc tính toán vòm mái về bài toán dầm trên nền đàn hồi, mở đường cho việc ứng dụng mái vòm trong kiến trúc.

Do việc tích phân hệ phương trình vi phân tổng quát gặp nhiều khó khăn, các nhà khoa học đã phát triển nhiều lý thuyết tính toán đơn giản hóa:

  • Lý thuyết hiệu ứng biên, lý thuyết bán mômen và lý thuyết vỏ thoải: Giúp giải quyết bài toán cho các dạng mặt có độ cong bằng không, mặt tròn xoay, mặt tịnh tiến và mặt bậc hai.
  • Nghiên cứu phi tuyến: K. Ong nghiên cứu sự ổn định của vỏ có độ võng lớn (phi tuyến hình học); Yu. Rabotnov đạt được các kết quả trong bài toán phi tuyến vật lý của lý thuyết vỏ.
  • Các hướng nghiên cứu chuyên sâu khác: V. Bolotin nghiên cứu lĩnh vực đàn hồi khí, động học và ổn định của vỏ; A. Iliouchin khảo sát biến dạng dẻo, ổn định và động học ngoài miền đàn hồi, vỏ đàn nhớt; A.R. Rzhanitsyn phát triển phương pháp tính tấm và vỏ theo trạng thái giới hạn; V.Z. Vlasov đề xuất các phương pháp tính gần đúng vỏ trụ dựa trên việc bỏ qua mômen uốn dọc, mômen xoắn và biến dạng trượt của mặt trung bình; B.B. Novozhilov và A.L. Goldenveizer chứng minh sai số của các giả thiết cơ bản trong lý thuyết vỏ mỏng có cùng bậc với bình phương bề dày tương đối $(h/R)^2$.

Vị trí của chuyên đề nằm ở việc tổng hợp có hệ thống các lý thuyết giải tích kinh điển, làm rõ phạm vi áp dụng của lý thuyết phi mômen và hiệu ứng biên, từ đó xây dựng thuật toán giải cho các hệ vỏ liên hợp phức tạp (gồm nhiều loại kết cấu vỏ ghép nối qua vành biên) và so sánh trực tiếp với phương pháp phần tử hữu hạn số trị.

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

Chuyên đề dựa trên nền tảng lý thuyết đàn hồi tuyến tính và lý thuyết vỏ mỏng cổ điển với ba giả thiết cơ bản (Kirchhoff - Love):

  1. Giả thiết phần tử thẳng (Kirchhoff - Love): Đoạn pháp tuyến thẳng góc với mặt trung bình trước biến dạng vẫn thẳng và vuông góc với mặt trung bình sau biến dạng, chiều dài không đổi: $\epsilon_{\alpha z} = \epsilon_{\beta z} = \epsilon_z = 0$.
  2. Giả thiết ứng suất pháp: Ứng suất pháp theo phương chiều dày $\sigma_z$ rất nhỏ so với các ứng suất khác nên được bỏ qua trong tính toán ($\sigma_z = 0$).
  3. Giả thiết biến dạng bé và đàn hồi tuyến tính: Biến dạng của vỏ là bé, vật liệu làm việc hoàn toàn tuân theo định luật Hooke.

Phương pháp nghiên cứu thực hiện bao gồm:

  • Phương pháp giải tích vi tích phân: Thiết lập và giải hệ phương trình vi phân cân bằng tĩnh học, phương trình hình học và phương trình vật lý.
  • Phương pháp lực và phương pháp chuyển vị: Sử dụng hệ phương trình chính tắc để khử các ẩn số lực thừa hoặc chuyển vị nút tại các vị trí liên kết biên và liên kết giữa các vỏ thành phần.
  • Lý thuyết hiệu ứng biên Meisner - Zimmerman: Sử dụng toán tử Meisner và các hàm suy giảm Zimmerman ($\eta_1, \eta_2, \eta_3, \eta_4$) để giải bài toán uốn cục bộ tại mép vỏ.
  • Phương pháp tương tự tĩnh học - hình học (Goldenveizer): Thiết lập hàm liên tục $\Phi_1, \Phi_2, \Phi_3$ để chuyển đổi tương đương giữa hệ phương trình liên tục biến dạng và phương trình cân bằng nội lực.
  • Mô hình dầm trên nền đàn hồi: Ứng dụng mô hình nền Winkler với hệ số nền $k$ để giải bản đáy bể chứa nước và dầu.
  • Phương pháp phần tử hữu hạn (FEM): Sử dụng phần mềm SAP2000 Nonlinear Version 7.42 để phân tích số trị 3D, chia lưới phần tử vỏ (shell elements), giải các trường hợp tải trọng và tổ hợp tải trọng nhằm đối chiếu với kết quả giải tích.

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

Chương I: Một số kiến thức cơ bản về hình học của vỏ mỏng

Chương I trình bày các khái niệm hình học vi phân mô tả mặt trung bình của vỏ:

  • Đặc trưng hình học: Vỏ là vật thể không gian 3 chiều có bề dày $h$ nhỏ hơn rất nhiều so với hai kích thước còn lại. Mặt nằm cách đều hai mặt giới hạn là mặt trung bình. Phân loại theo bề dày tương đối: $h/R \le 0,05$ là vỏ mỏng, $h/R > 0,05$ là vỏ dày ($R$ là bán kính cong chính nhỏ nhất).
  • Độ cong Gauss ($k = k_1 k_2$):
    • Vỏ có độ cong Gauss dương ($k > 0$): mặt cầu, mặt ellipxoid.
    • Vỏ có độ cong Gauss bằng không ($k = 0$): mặt trụ, mặt nón.
    • Vỏ có độ cong Gauss âm ($k < 0$): mặt hypecbolic.
  • Hệ tọa độ cong và dạng toàn phương:
    • Dạng phương trình mặt véc-tơ: $\vec{r} = \vec{r}(\alpha, \beta) = x(\alpha, \beta)\vec{i} + y(\alpha, \beta)\vec{j} + z(\alpha, \beta)\vec{k}$.
    • Dạng toàn phương thứ nhất (xác định vi phân độ dài cung $ds$): $$ds^2 = E d\alpha^2 + 2F d\alpha d\beta + G d\beta^2 = A^2 d\alpha^2 + B^2 d\beta^2$$ Trong đó $A, B$ là các hệ số Lamé biến đổi gia số tọa độ thành gia số chiều dài ($ds_\alpha = A d\alpha$, $ds_\beta = B d\beta$), $C = 0$ khi hệ tọa độ cong trực giao.
    • Dạng toàn phương thứ hai (xác định độ cong tiết diện pháp tuyến $k$): $$k = \frac{L d\alpha^2 + 2M d\alpha d\beta + N d\beta^2}{E d\alpha^2 + 2F d\alpha d\beta + G d\beta^2}$$ Sáu hệ số $E, F, G, L, M, N$ hoàn toàn xác định bề mặt hình học và thỏa mãn điều kiện tương thích Gauss - Codazzi.

Chương II: Hệ phương trình cơ bản của lý thuyết vỏ mỏng và phương pháp giải

Chương II thiết lập hệ 17 phương trình cơ bản giải bài toán vỏ mỏng:

  • Các thành phần ứng lực:
    • Ứng lực màng (nằm trong mặt phẳng tiếp tuyến): lực dọc $N_\alpha, N_\beta$, lực trượt $S_{\alpha\beta}, S_{\beta\alpha}$.
    • Lực uốn xoắn: mômen uốn $M_\alpha, M_\beta$, lực cắt $Q_\alpha, Q_\beta$, mômen xoắn $M_{\alpha\beta}, M_{\beta\alpha}$.
    • Đối với vỏ mỏng, tỉ số $z/R \ll 1$, lực trượt và mômen xoắn có tính đối ứng: $S_{\alpha\beta} \approx S_{\beta\alpha} = S$, $M_{\alpha\beta} \approx M_{\beta\alpha} = H$.
  • Hệ thống 17 phương trình chứa 17 ẩn số:
    1. 5 phương trình cân bằng tĩnh học: 3 phương trình hình chiếu lực lên các trục $\alpha, \beta, z$ và 2 phương trình cân bằng mômen (phương trình $\sum M_z = 0$ đồng nhất tự thỏa mãn).
    2. 6 phương trình hình học (liên hệ biến dạng - chuyển vị): Xác định các biến dạng dài $\epsilon_{\alpha 0}, \epsilon_{\beta 0}$, biến dạng góc $\gamma_{\alpha\beta 0}$, biến dạng độ cong $\chi_\alpha, \chi_\beta$ và biến dạng xoắn $\chi_{\alpha\beta}$ theo các chuyển vị $u_0, v_0, w_0$ của mặt trung bình.
    3. 6 phương trình vật lý: Thiết lập quan hệ giữa 8 ứng lực và các biến dạng dựa trên định luật Hooke và độ cứng uốn hình trụ $D = \frac{Eh^3}{12(1-\nu^2)}$.
  • Đường lối giải bài toán đàn hồi:
    • Giải theo nội lực: Rút các lực cắt $Q_\alpha, Q_\beta$, kết hợp phương trình cân bằng với điều kiện liên tục biến dạng kiểu Beltrami - Michell để thu được hệ 6 phương trình chứa 6 ẩn số nội lực ($N_\alpha, N_\beta, M_\alpha, M_\beta, S, H$).
    • Giải theo chuyển vị: Biểu diễn nội lực qua chuyển vị $u, v, w$, dẫn đến hệ phương trình vi phân toán tử $L_{ij} u_j = f_i$ ($i, j = 1 \div 3$).
    • Phương pháp tương tự tĩnh học - hình học: Dựa trên sự tương đương cấu trúc giữa phương trình cân bằng và phương trình liên tục biến dạng, sử dụng các hàm thế $\Phi_1, \Phi_2, \Phi_3$.
    • Điều kiện biên: Hệ phương trình vi phân bậc 8 đối với mỗi phương tọa độ đòi hỏi 4 điều kiện biên trên mỗi cạnh (tổng cộng 16 điều kiện cho 4 biên: biên ngàm, biên tựa di động, biên tựa trượt pháp tuyến, biên khớp cố định).

Chương III: Một số lý thuyết vỏ

Chương III phân tích chi tiết hai trạng thái ứng suất cơ bản của vỏ tròn xoay:

1. Lý thuyết phi mômen vỏ đàn hồi tròn xoay

  • Giả thiết ứng suất pháp và tiếp phân bố đều theo chiều dày vỏ: $M_\alpha = M_\beta = H = 0$, $Q_\alpha = Q_\beta = 0$. Hệ phương trình vi phân tĩnh định gồm 3 phương trình xác định 3 ứng lực màng $N_\phi, N_\theta, S$.
  • Điều kiện tồn tại trạng thái phi mômen: Mặt vỏ trơn liên tục, bề dày không đổi hoặc biến đổi đều đặn; tải trọng phân bố liên tục, không có lực tập trung; liên kết biên cho phép chuyển vị tự do theo phương pháp tuyến; phản lực biên nằm trong mặt phẳng tiếp tuyến.
  • Với tải trọng đối xứng trục ($S = 0, q_\theta = 0$), lực dọc kinh tuyến $N_\phi$ được tích phân trực tiếp: $$N_\phi = \frac{V_0}{2\pi R_2 \sin^2\phi}$$ Trong đó $V_0$ là tổng hợp lực theo phương thẳng đứng của ngoại lực tác dụng lên phần vỏ cắt ra. Lực vòng $N_\theta$ được tính từ phương trình Laplace: $$\frac{N_\phi}{R_1} + \frac{N_\theta}{R_2} = q_n$$
  • Công thức chuyển vị chân vỏ quy ước ($\Delta_{a0}, \vartheta_{a0}$) được thiết lập cho vỏ trụ đứng và vỏ chỏm cầu dưới tác dụng của trọng lượng bản thân và áp lực chất lỏng.

2. Lý thuyết mômen của vỏ tròn xoay và hiệu ứng biên

  • Khi điều kiện biên cản trở biến dạng màng hoặc có mômen uốn cục bộ, bài toán được giải bằng hệ phương trình vi phân cấp 2 liên hợp Meisner: $$L(V) + \nu V = Eh R_1 \psi + \Phi(\phi)$$ $$L(\psi) - \nu \psi = -\frac{R_1}{D} V$$
  • Hiệu ứng biên: Đối với vỏ không quá thoải ($\phi > 20^\circ \div 30^\circ$), mômen uốn tắt rất nhanh khi ra xa mép biên. Bằng cách đặt biến không thứ nguyên $\xi = t/s$ (với thông số đặc trưng chiều dài suy giảm $s = \frac{\sqrt{R h}}{\sqrt[4]{12(1-\nu^2)}}$), phương trình vi phân uốn đưa về dạng: $$\frac{d^4 M}{d\xi^4} + 4M = 0$$
  • Nghiệm tổng quát suy giảm biểu diễn qua 4 hàm số Zimmerman:
    • $\eta_1(\xi) = e^{-\xi}\cos\xi$
    • $\eta_2(\xi) = e^{-\xi}\sin\xi$
    • $\eta_3(\xi) = e^{-\xi}(\cos\xi + \sin\xi)$
    • $\eta_4(\xi) = e^{-\xi}(\cos\xi - \sin\xi)$
  • Khi $\xi > 3$, ảnh hưởng của mômen uốn biên coi như triệt tiêu. Tác giả thiết lập ma trận độ mềm/độ cứng và các chuyển vị đơn vị quy ước ở chân vỏ tròn xoay có chu vi ngàm ($a_{11} = s$, $a_{12} = a_{21} = -\frac{s^2}{2}\sin\phi_0$, $a_{22} = \frac{s^3}{2}\sin^2\phi_0$).

Chương IV: Giải một số bài toán của kết cấu vỏ

Bài toán 1: Bể chứa nước ngầm hình trụ đặt trong nền đất

  • Thông số kết cấu và vật liệu: Đường kính trong $D = 10\text{ m}$ ($R = 5\text{ m}, R_t = 5,3\text{ m}$), chiều cao thành $H = 6\text{ m}$, bề dày nắp $h_{nắp} = 0,15\text{ m}$, bề dày thành $\delta = 0,20\text{ m}$, bề dày lớp đất đắp trên nắp $h_{đất} = 0,5\text{ m}$. Bê tông mác 300 ($R_n = 1,6\text{ kN/cm}^2, R_k = 0,095\text{ kN/cm}^2$), thép CT5 ($R_a = 27\text{ kN/cm}^2$). Địa chất: trọng lượng riêng đất $\gamma_d = 17\text{ kN/m}^3$, góc ma sát trong $\phi = 32^\circ$, hệ số nền $k = 3$.
  • Tải trọng tính toán: Tải nắp và đất đắp $p = 13,5\text{ kN/m}^2$; áp lực nước đáy bể $p_n = 60\text{ kN/m}^2$; áp lực bên của đất phân bố từ $p_1 = 2,87\text{ kN/m}^2$ (đỉnh) đến $p_2 = 37,3\text{ kN/m}^2$ (đáy).
  • Quy trình giải giải tích:
    • Tính toán tách biệt liên kết nắp - thành (biên trên) và liên kết thành - bản đáy trên nền đàn hồi (biên dưới) bằng phương pháp lực.
    • Giải hệ phương trình chính tắc tìm các phản lực mômen uốn biên và lực đẩy ngang:
      • Tại đỉnh thành: Do trọng lượng nắp + đất đắp: $M_{10} = -45,92\text{ kNm}$, $H_{10} = -58,78\text{ kN}$; do áp lực đất: $M_{10} = -0,067\text{ kNm}$, $H_{10} = -1,206\text{ kN}$.
      • Tại chân thành: Do trọng lượng nắp + thành: $M_{21} = -27,83\text{ kNm}$, $H_{21} = -35,70\text{ kN}$; do áp lực đất: $M_{22} = -1,91\text{ kNm}$, $H_{22} = -17,0\text{ kN}$; do áp lực nước: $M_{23} = 3,04\text{ kNm}$, $H_{23} = -27,32\text{ kN}$.
    • Tính toán cốt thép chịu lực và kiểm tra chống nứt cho từng cấu kiện:
      • Nắp bể: Cốt thép hướng kính bố trí $\phi 16a100$ ($F_a = 16\text{ cm}^2/\text{m}$ tại mép ngàm) và $\phi 12a100$ ($F_a = 8,2\text{ cm}^2/\text{m}$ tại nhịp); thép vòng bố trí $\phi 12a100$.
      • Thành bể: Cốt thép phương đứng bố trí $\phi 12a100$ ($F_a = 9,94\text{ cm}^2/\text{m}$, đoạn trên) và $\phi 12a150$ ($F_a = 5,8\text{ cm}^2/\text{m}$, đoạn dưới); cốt thép vòng bố trí $\phi 12a100$ (đoạn $0 \div 3\text{ m}$, $N_\beta = 283,26\text{ kN/m}$) và $\phi 14a100$ (đoạn $3 \div 6\text{ m}$, $N_\beta = 355,29\text{ kN/m}$), đều thỏa mãn điều kiện không nứt.
      • Bản đáy: Cốt thép phương bán kính bố trí $\phi 12a100$ ($F_a = 8,8\text{ cm}^2/\text{m}$), cốt thép vòng bố trí $\phi 12a100$.

Bảng tổng hợp nội lực và bố trí cốt thép bể chứa nước ngầm (Bài toán 1)

Bộ phận kết cấu Vị trí kiểm tra Mômen uốn $M$ ($\text{kNm/m}$) Lực dọc/vòng $N$ ($\text{kN/m}$) Cốt thép tính toán $F_a$ ($\text{cm}^2/\text{m}$) Bố trí cốt thép thực tế
Nắp tròn Mép ngàm biên $45,92$ - $16,00$ $\phi 16a100$
Giữa nhịp $25,18$ - $8,20$ $\phi 12a100$
Thành trụ Đỉnh thành (thép đứng) $45,92$ $35,78$ $9,94$ $\phi 12a100$
Đáy thành (thép đứng) $29,74$ $68,78$ $5,80$ $\phi 12a150$
Vòng (đoạn cao $0 \div 3\text{ m}$) - $283,26$ $10,50$ $\phi 12a100$
Vòng (đoạn cao $3 \div 6\text{ m}$) - $355,29$ $13,16$ $\phi 14a100$
Bản đáy Mép biên (thép bán kính) $27,70$ $51,70$ $8,80$ $\phi 12a100$
Thép vòng $2,42$ $8,38$ $1,20$ $\phi 10a100$

Bài toán 2: Bể chứa dầu mỏ nắp chỏm cầu, thành trụ và đáy phẳng

  • Thông số kết cấu và vật liệu: Bể đặt trên nền đất, bán kính vỏ trụ $R = 6\text{ m}$, chiều cao thành $H = 4,7\text{ m}$, bề dày thành đỉnh $\delta_1 = 0,10\text{ m}$, bề dày thành đáy $\delta_2 = 0,22\text{ m}$ (trung bình $\delta_0 = 0,16\text{ m}$). Nắp chỏm cầu có bán kính cong $r = 11,1\text{ m}$, góc ôm $\phi_0 = 32^\circ 10'$, bề dày $\delta_1 = 0,10\text{ m}$. Vành đai liên kết nút nắp - thành có tiết diện chữ nhật $b \times h = 30 \times 40\text{ cm}$. Bản đáy phẳng dày $\delta_3 = 0,20\text{ m}$. Bê tông mác 200 ($R_n = 1,15\text{ kN/cm}^2, R_k = 0,075\text{ kN/cm}^2$), thép CT3 ($R_a = 21\text{ kN/cm}^2$). Dầu có trọng lượng riêng $\gamma_{dầu} = 9\text{ kN/m}^3$.
  • Quy trình giải và kết quả:
    • Tại nút liên kết đỉnh (chỏm cầu + vành đai + thành trụ): Áp dụng phương pháp chuyển vị với 2 ẩn số chuyển vị nút: góc xoay $\alpha = 0,1075$ và chuyển vị ngang $u = -0,163$.
    • Xác định nội lực truyền vào các cấu kiện tại nút:
      • Vỏ trụ: $M' = -44,09\text{ kNm}$, $H' = -9,88\text{ kN}$.
      • Vành đai: $M''' = 0,26\text{ kNm}$, $H''' = 34,6\text{ kN}$ (lực kéo vòng $Z = 207,6\text{ kN}$).
      • Nắp vỏ cầu: $M'' = -13,4 - 16,6 = -30,0\text{ kNm}$, $H'' = 16,47\text{ kN}$.
    • Tại liên kết chân thành trụ với bản đáy: Giải theo phương pháp lực, kết hợp mô hình dầm trên nền đàn hồi.
    • Bố trí cốt thép:
      • Nắp chỏm cầu: Cốt thép vòng $\phi 8a100 \div \phi 8a200$; cốt thép phương kinh tuyến $\phi 8a150$ ($F_a = 2,46\text{ cm}^2/\text{m}$).
      • Vành đai liên kết: $F_a = 9,88\text{ cm}^2$, bố trí đai và thép dọc chịu kéo vòng.
      • Thành bể: Thép đứng $\phi 8a200$ (đoạn trên) và $\phi 8a100$ (đoạn dưới); thép vòng $\phi 10a100$ ($F_a = 6,2\text{ cm}^2/\text{m}$, kiểm tra chống nứt $N_{k,gh} = 186\text{ kN} > 139\text{ kN}$).
      • Bản đáy: Thép phương bán kính $\phi 10a100$ ($F_a = 17,5\text{ cm}^2/\text{m}$), thép vòng $\phi 10a200$ ($F_a = 1,11\text{ cm}^2/\text{m}$).

Phụ lục: Kết quả phân tích phần tử hữu hạn (SAP2000 Nonlinear v7.42)

Phụ lục chuyên đề cung cấp toàn bộ kết quả xuất từ phần mềm SAP2000 Nonlinear Version 7.42 cho hai bài toán:

  • Tệp xuất dữ liệu Bài toán 1 (bai toan 1 MKN.OUT): Định nghĩa 357 nút, các điều kiện biên ngàm/tựa tự do, bảng cân bằng lực tổng thể (Global Force Balance) cho các trường hợp tải trọng riêng lẻ: trọng lượng nắp đất đắp (DATNAP), áp lực nước (NUOC), áp lực đất hông (DHONG), tổ hợp nắp và nước (NAPNUOC), và tổ hợp toàn phần (TONG). Các biểu đồ nội lực phần tử vỏ 2D/3D gồm lực màng $F_{11}, F_{22}$, mômen uốn $M_{11}, M_{22}$ và hình dạng biến dạng (Deformed Shape) được in kèm chi tiết.
  • Tệp xuất dữ liệu Bài toán 2 (bai toan 2.OUT): Định nghĩa mô hình nút từ 130 đến 822, giải phóng bậc tự do và phân tích phi tuyến với tải trọng dầu (LOAD1DAU), tải trọng nắp (LOAD2NAP) và tổ hợp (COMB1). Kết quả số trị thể hiện sự trùng khớp về quy luật phân bố và giá trị cực trị của mômen uốn cục bộ tại các vị trí nút giao vành đai và đáy bể so với lời giải giải tích của lý thuyết hiệu ứng biên.

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

Chuyên đề đã tổng kết và đạt được các kết quả cụ thể:

  1. Hệ thống hóa lý thuyết giải tích: Tổng hợp hoàn chỉnh hệ thống 17 phương trình vi phân cơ bản của lý thuyết vỏ mỏng đàn hồi từ nền tảng hình học vi phân và 3 giả thiết Kirchhoff - Love, làm sáng tỏ các mối quan hệ tương hỗ giữa phương pháp lực, phương pháp chuyển vị và phương pháp tương tự tĩnh học - hình học.
  2. Thiết lập quy trình tính toán hiệu ứng biên cho vỏ liên hợp: Xây dựng chi tiết phương pháp giải bán giải tích kết hợp giữa lý thuyết phi mômen và lý thuyết hiệu ứng biên (sử dụng toán tử vi phân Meisner và chuỗi hàm Zimmerman suy giảm nhanh) cho các kết cấu vỏ tròn xoay phức hợp có xét đến vành đai liên kết trung gian.
  3. Giải quyết hoàn chỉnh hai bài toán công trình thực tế: Thực hiện tính toán cụ thể từ khâu xác định tải trọng, lập phương trình chính tắc, tìm nội lực giải tích, kiểm tra điều kiện bền, tính toán diện tích cốt thép đến kiểm tra khả năng chống nứt của bê tông cốt thép cho bể chứa nước ngầm $D = 10\text{ m}$ và bể chứa dầu mỏ $D = 12\text{ m}$.
  4. Ứng dụng và kiểm chứng số trị bằng FEM: Ứng dụng chương trình SAP2000 Nonlinear v7.42 để lập mô hình số 3D, phân tích nội lực và biến dạng của hệ kết cấu vỏ, chứng minh tính tin cậy và sự phù hợp giữa phương pháp giải tích giải gần đúng hiệu ứng biên với phương pháp phần tử hữu hạn.

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

Dựa trên phạm vi nghiên cứu được công bố trong văn bản:

  • Hạn chế: Chuyên đề chỉ tập trung nghiên cứu trạng thái biến dạng bé và đàn hồi tuyến tính đối xứng trục của vật liệu đẳng hướng. Các ảnh hưởng của tính phi tuyến hình học (độ võng lớn), phi tuyến vật lý (trạng thái dẻo, đàn nhớt của bê tông), tải trọng động lực học, tải trọng nhiệt độ và các dạng tải trọng phi đối xứng chưa được khảo sát trong tính toán giải tích.
  • Hướng nghiên cứu tiếp: Mở rộng phương pháp giải cho các dạng vỏ thoải có mặt cong phức tạp chịu tải trọng bất kỳ; nghiên cứu tính toán vỏ ngoài miền đàn hồi và phân tích ổn định động học của kết cấu vỏ mỏng không gian dưới tác động của động đất và tải trọng gió thay đổi.

Giá trị tham khảo

Chuyên đề tiến sĩ là tài liệu tham khảo chuyên khảo dành cho:

  • Nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Cơ học kết cấu, Xây dựng dân dụng và công nghiệp, Công trình thủy, Xây dựng cầu hầm: Tham khảo phương pháp thiết lập phương trình vi phân mặt cong, các giả thiết cơ bản và kỹ thuật giải hệ phương trình vi phân bậc cao bằng hàm Zimmerman.
  • Kỹ sư thiết kế kết cấu công trình: Cung cấp quy trình từng bước tính toán kết cấu vỏ bê tông cốt thép thực tế (bể nước, si-lô, bể chứa dầu mỏ), từ việc tính toán nội lực màng, nội lực uốn cục bộ tại biên đến công thức bố trí thép và kiểm tra bề rộng vết nứt.
  • Giảng viên giảng dạy các môn học: Sức bền vật liệu nâng cao, Lý thuyết đàn hồi, Lý thuyết tấm và vỏ mỏng, và Ứng dụng phương pháp phần tử hữu hạn trong kết cấu công trình.

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

1. Điều kiện hình học để một kết cấu được xem là vỏ mỏng trong công trình là gì?
Theo Mục 1-2-2 (Chương I), một kết cấu được coi là vỏ mỏng khi bề dày tương đối thỏa mãn điều kiện $h/R \le 0,05$ (trong đó $h$ là chiều dày vỏ và $R$ là bán kính cong chính nhỏ nhất). Nếu $h/R > 0,05$, kết cấu được xếp vào loại vỏ dày.

2. Hệ phương trình cơ bản của lý thuyết vỏ mỏng đàn hồi gồm bao nhiêu phương trình và ẩn số?
Theo Mục 2-2-3 (Chương II), hệ thống phương trình đầy đủ gồm 17 phương trình với 17 ẩn số, bao gồm:

  • 5 phương trình cân bằng tĩnh học.
  • 6 phương trình hình học (biến dạng - chuyển vị).
  • 6 phương trình vật lý (định luật Hooke). Các ẩn số gồm 8 thành phần ứng lực ($N_\alpha, N_\beta, S, M_\alpha, M_\beta, H, Q_\alpha, Q_\beta$), 6 thành phần biến dạng mặt trung bình ($\epsilon_{\alpha 0}, \epsilon_{\beta 0}, \gamma_{\alpha\beta 0}, \chi_\alpha, \chi_\beta, \chi_{\alpha\beta}$) và 3 thành phần chuyển vị ($u_0, v_0, w_0$).

3. Tại sao lực trượt và mômen xoắn trong vỏ mỏng được coi là đối ứng?
Trong lý thuyết vỏ tổng quát, $S_{\alpha\beta} \ne S_{\beta\alpha}$ và $M_{\alpha\beta} \ne M_{\beta\alpha}$ do ảnh hưởng của các số hạng tỉ số $z/R$. Tuy nhiên, đối với vỏ mỏng, do tỉ số $z/R \ll 1$ (có thể bỏ qua so với 1), các biểu thức tích phân cho phép rút gọn thành $S_{\alpha\beta} \approx S_{\beta\alpha} = S$ và $M_{\alpha\beta} \approx M_{\beta\alpha} = H$.

4. Hàm số Zimmerman đóng vai trò gì trong việc giải bài toán hiệu ứng biên?
Hàm Zimmerman ($\eta_1, \eta_2, \eta_3, \eta_4$) là các hàm suy giảm dạng tích của hàm mũ và hàm lượng giác ($e^{-\xi}\cos\xi, e^{-\xi}\sin\xi,...$) xuất hiện trong nghiệm của phương trình vi phân uốn mép bậc 4 ($\frac{d^4 M}{d\xi^4} + 4M = 0$). Các hàm này mô tả chính xác sự triệt tiêu nhanh chóng của mômen uốn biên khi đi vào sâu trong lòng vỏ (khi đại lượng không thứ nguyên $\xi = t/s > 3$, mômen uốn tắt hoàn toàn).

5. Trong Bài toán 2, tác giả đã giải quyết liên kết nút giữa nắp chỏm cầu, vành đai và thành trụ bằng phương pháp nào?
Tác giả sử dụng phương pháp chuyển vị (Mục 4-2-1), lập hệ phương trình chính tắc 2 ẩn số chuyển vị nút gồm góc xoay $\alpha$ và chuyển vị ngang $u$. Độ cứng uốn $EJ$ của nắp vỏ cầu được chọn làm chuẩn để tính toán các hệ số độ cứng quy ước của vành đai và thành vỏ trụ, từ đó xác định góc xoay $\alpha = 0,1075$ và chuyển vị ngang $u = -0,163$ để phân phối nội lực cho từng cấu kiện.

Kết luận

Chuyên đề tiến sĩ "Một số vấn đề của lý thuyết kết cấu vỏ mỏng đàn hồi" đã hoàn thành việc hệ thống hóa toàn diện cơ sở toán học và cơ học của lý thuyết vỏ mỏng đàn hồi tuyến tính. Tác giả đã làm rõ mối liên hệ giữa các phương pháp giải giải tích cổ điển và lý thuyết hiệu ứng biên Meisner - Zimmerman, áp dụng thành công vào việc tính toán, thiết kế cốt thép và kiểm tra nứt cho hai công trình bể chứa chất lỏng thực tế. Kết quả phân tích giải tích được đối chiếu và kiểm chứng độ tin cậy thông qua mô hình phần tử hữu hạn trên phần mềm SAP2000 Nonlinear Version 7.42, tạo nên tài liệu khoa học có giá trị học thuật và ứng dụng thực tiễn trong ngành Cơ học kết cấu và Xây dựng công trình.