Giáo Trình Phân Tích Ứng Suất (Lý Trường Thành) – Tài Liệu Học Thuật Cơ Học Vật Rắn Biến Dạng


Tổng quan về giáo trình

Giáo trình Phân Tích Ứng Suất do tác giả Lý Trường Thành biên soạn, được Nhà xuất bản Khoa học Tự nhiên và Công nghệ xuất bản năm 2010 (Chịu trách nhiệm xuất bản: GS. TSKH Nguyễn Khoa Sơn; Biên tập: Trần Phương Đông; Số xuất bản: 295-2010/CXB/026-02/KHTNCN). Tài liệu được hoàn thiện với sự đóng góp chuyên môn từ tập thể giảng viên Bộ môn Sức bền vật liệu – Cơ học kết cấu, Trường Đại học Thủy lợi.

Trong chương trình đào tạo kỹ thuật công trình và cơ học, môn học Phân tích ứng suất là học phần chuyên sâu thuộc ngành Cơ học vật rắn biến dạng, đóng vai trò kế thừa và mở rộng các nguyên lý từ môn Sức bền vật liệu và Cơ học kết cấu. Mục tiêu học tập cốt lõi của giáo trình là trang bị cho người học phương pháp xác định trạng thái ứng suất, trường biến dạng và trường chuyển vị của các vật thể đàn hồi chịu tác dụng của tải trọng ngoài khi vật liệu làm việc trong giới hạn đàn hồi tuyến tính.

Về cấu trúc và phương pháp tiếp cận, giáo trình tích hợp có hệ thống giữa hai phân môn: Lý thuyết đàn hồi ứng dụng (tiếp cận giải tích) và Phương pháp phần tử hữu hạn (tiếp cận số trị). Điểm đặc trưng của tài liệu là mở rộng đối tượng nghiên cứu từ kết cấu dạng thanh (1 chiều) sang các kết cấu dạng tấm phẳng, bản dầm tường, vỏ và khối không gian (2 chiều và 3 chiều). Đồng thời, tài liệu tập trung giải quyết các bài toán cơ học phục vụ trực tiếp cho thiết kế và kiểm tra an toàn công trình xây dựng, đặc biệt là công trình Thủy lợi (thân đập bê tông, tường chắn đất, đường ống áp lực và kết cấu ngầm).

                      HỆ THỐNG MÔN HỌC
           ┌──────────────────────────────────────┐
           │     Cơ học vật rắn biến dạng         │
           └──────────────────┬───────────────────┘
                              │
         ┌────────────────────┴────────────────────┐
         ▼                                         ▼
┌──────────────────┐                     ┌──────────────────┐
│ Sức bền vật liệu │                     │  Cơ học kết cấu  │
│ (Vật thể thanh)  │                     │ (Hệ thanh/khung) │
└────────┬─────────┘                     └────────┬─────────┘
         │                                        │
         └────────────────────┬───────────────────┘
                              ▼
           ┌──────────────────────────────────────┐
           │        PHÂN TÍCH ỨNG SUẤT            │
           │  (Dạng thanh, tấm, vỏ, khối 1D-2D-3D)│
           └──────────────────┬───────────────────┘
                              │
         ┌────────────────────┴────────────────────┐
         ▼                                         ▼
┌─────────────────────────┐               ┌─────────────────────────┐
│  Lý thuyết đàn hồi      │               │ Phương pháp phần tử     │
│  ứng dụng (Giải tích)   │               │ hữu hạn (Mô hình số)    │
└─────────────────────────┘               └─────────────────────────┘

Nội dung kiến thức cốt lõi

Các chương và chủ đề chính

Nội dung giáo trình được tổ chức thành 4 chương chuyên đề chính cùng hệ thống bài tập thực hành và tài liệu tham khảo:

                                  MỤC LỤC TỔNG QUAN
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ MỞ ĐẦU: Phạm vi, đối tượng, giả thiết đàn hồi tuyến tính                                │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ CHƯƠNG 1: Những phương trình cơ bản của bài toán đàn hồi tuyến tính đẳng hướng         │
│   ├─ 1.1 Lý thuyết ứng suất (Navier, Saint-Venant, 3 ứng suất chính, bất biến I1,I2,I3)│
│   ├─ 1.2 Lý thuyết về biến dạng (Cauchy, 6 phương trình liên tục Saint-Venant)         │
│   ├─ 1.3 Quan hệ ứng suất - biến dạng (Định luật Hooke, hằng số Lamé, ma trận [D])     │
│   ├─ 1.4 Phương pháp giải bài toán đàn hồi (Beltrami-Michell, Lamé, nửa ngược)         │
│   └─ 1.5 Nguyên lý công và năng lượng (Công khả dĩ, Castigliano, Reissner-Hellinger)   │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ CHƯƠNG 2: Bài toán phẳng trong hệ toạ độ vuông góc                                     │
│   ├─ 2.1 Ứng suất phẳng & Biến dạng phẳng (Hằng số E*, μ*)                             │
│   ├─ 2.2 Hệ phương trình cơ bản (Navier, Cauchy, Saint-Venant, Hooke dạng phẳng)       │
│   ├─ 2.3 Phương trình song điều hoà & Hàm ứng suất Airy ∇⁴φ = 0                        │
│   └─ 2.4 Các bài toán tiêu biểu (Dầm conson, đập tam giác Lévy, đập chữ nhật,          │
│          dầm tường Filon - Ribiere)                                                    │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ CHƯƠNG 3: Bài toán phẳng trong hệ toạ độ cực                                           │
│   ├─ 3.1 Biến đổi toạ độ và toán tử Laplace dạng cực                                   │
│   ├─ 3.2 Phương trình vi phân cân bằng & Quan hệ hình học dạng cực                     │
│   ├─ 3.3 Hàm Airy trong toạ độ cực & Phương pháp nửa ngược                             │
│   ├─ 3.4 Bài toán đối xứng trục / Không phụ thuộc góc cực θ                            │
│   └─ 3.5 Các bài toán kỹ thuật (Thanh cong Golovin, ống dày Lamé, lỗ khoét tròn,       │
│          nêm phẳng, bán không gian đàn hồi)                                            │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ CHƯƠNG 4: Giải bài toán lý thuyết đàn hồi bằng phương pháp số (Phần tử hữu hạn)        │
│   ├─ 4.1 Rời rạc hoá kết cấu, chuyển vị nút, ma trận độ cứng phần tử [k]e              │
│   ├─ 4.2 Phần tử tam giác ứng suất phẳng (CST), hàm dạng, thế năng toàn phần           │
│   ├─ 4.3 Phần tử thanh kéo nén, uốn phẳng, thanh không gian phức tạp                   │
│   └─ 4.4 Ma trận độ cứng tổng thể [K], xử lý điều kiện biên, tối ưu băng thông         │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ PHỤ LỤC & TÀI LIỆU THAM KHẢO                                                           │
└────────────────────────────────────────────────────────────────────────────────────────┘

Chương 1: Những phương trình cơ bản của bài toán đàn hồi tuyến tính đẳng hướng

Chương mở đầu thiết lập hệ thống công thức nền tảng trong không gian 3 chiều:

  • Lý thuyết ứng suất: Định nghĩa 9 thành phần ứng suất trên phân tố hình hộp vi phân ($3$ ứng suất pháp $\sigma_x, \sigma_y, \sigma_z$ và $6$ ứng suất tiếp $\tau_{xy}, \tau_{yx}, \tau_{yz}, \tau_{zy}, \tau_{xz}, \tau_{zx}$), luật đối ứng ứng suất tiếp $\tau_{ij} = \tau_{ji}$, thiết lập hệ phương trình vi phân cân bằng Navier dạng ma trận $CS + P = 0$. Khảo sát phân tố tứ diện để suy ra điều kiện biên lực mặt nghiêng $P_v = S \cdot L$. Phân tích phương trình đặc trưng bậc 3 để tìm 3 ứng suất chính $\sigma_1 \ge \sigma_2 \ge \sigma_3$ thông qua các bất biến ứng suất $I_1, I_2, I_3$ và tính toán ứng suất tiếp cực trị.
  • Lý thuyết biến dạng: Chuyển vị điểm $(u, v, w)$, thiết lập hệ phương trình hình học Cauchy ($\varepsilon = AU$) liên hệ giữa 6 thành phần biến dạng ($\varepsilon_x, \varepsilon_y, \varepsilon_z, \gamma_{xy}, \gamma_{yz}, \gamma_{zx}$) và các đạo hàm chuyển vị. Dẫn xuất 6 phương trình tương thích liên tục Saint-Venant nhằm đảm bảo tính liên tục của môi trường biến dạng.
  • Định luật Hooke tổng quát: Biểu diễn mối quan hệ vật lý qua môđun đàn hồi $E$, hệ số Poisson $\mu$, môđun trượt $G = \frac{E}{2(1+\mu)}$, hằng số Lamé $\lambda = \frac{\mu E}{(1+\mu)(1-2\mu)}$ và thiết lập ma trận đàn hồi $[D]$.
  • Phương pháp giải và tính duy nhất nghiệm: Trình bày phương pháp giải theo ứng suất (hệ Beltrami-Michell), giải theo chuyển vị (hệ Lamé), phương pháp nửa ngược Saint-Venant; kèm theo chứng minh toán học về tính duy nhất của nghiệm bài toán đàn hồi bằng phương pháp phản chứng.
  • Nguyên lý công và năng lượng: Nguyên lý công khả dĩ Lagrange, công bù khả dĩ, biểu thức thế năng biến dạng $W(u)$, thế năng bù $B(\sigma)$, các công thức vi phân Green, Castigliano, nguyên lý cực tiểu thế năng Lagrange ($\delta J = 0$), nguyên lý thế năng bù Castigliano ($\delta U = 0$) và phiếm hàm biến phân hỗn hợp Reissner-Hellinger ($\delta R = 0$).

Chương 2: Bài toán phẳng trong hệ toạ độ vuông góc

Chương này chuẩn hóa bài toán 3D về bài toán 2D trên mặt phẳng $Oxy$:

  • Phân loại bài toán phẳng:
    • Trạng thái ứng suất phẳng: Vật thể dạng tấm mỏng chiều dày $t$, $\sigma_z = \tau_{xz} = \tau_{yz} = 0$, $\varepsilon_z = -\frac{\mu}{E}(\sigma_x + \sigma_y) \neq 0$.
    • Trạng thái biến dạng phẳng: Vật thể lăng trụ dài (đập nước, hầm ngầm), $w = 0, \varepsilon_z = \gamma_{xz} = \gamma_{yz} = 0$, $\sigma_z = \mu(\sigma_x + \sigma_y)$. Đưa ra công thức chuyển đổi thông số đàn hồi tương đương $E^* = \frac{E}{1-\mu^2}$ và $\mu^* = \frac{\mu}{1-\mu}$.
  • Phương pháp hàm ứng suất Airy: Chuyển hệ phương trình vi phân cân bằng và tương thích về phương trình song điều hòa thuần nhất: $$\nabla^2 \nabla^2 \phi = \frac{\partial^4 \phi}{\partial x^4} + 2\frac{\partial^4 \phi}{\partial x^2 \partial y^2} + \frac{\partial^4 \phi}{\partial y^4} = 0$$
  • Lời giải các bài toán thực tế:
    1. Dầm công sôn: Chịu lực tập trung $P$ ở đầu tự do, đối chiếu độ chính xác với công thức Sức bền vật liệu.
    2. Đập mặt cắt tam giác (Lời giải Lévy): Xét đồng thời áp lực nước tĩnh và trọng lượng bản thân đập, xác định trường ứng suất phân bố theo hàm bậc nhất của toạ độ.
    3. Đập mặt cắt chữ nhật: Giải theo đa thức bậc cao, phân tích sai số tại vùng đỉnh và kiểm chứng thực nghiệm cho nguyên lý Saint-Venant.
    4. Dầm tường (Lời giải Filon – Ribiere): Áp dụng chuỗi lượng giác Fourier biểu diễn hàm Airy để giải bài toán dầm chịu tải trọng biên phân bố phức tạp.

Chương 3: Bài toán phẳng trong hệ toạ độ cực

Khảo sát các kết cấu có biên cong hoặc đối xứng tròn thông qua hệ toạ độ $(r, \theta)$:

  • Thiết lập hệ phương trình vi phân cân bằng, hệ phương trình hình học Cauchy và toán tử Laplace trong toạ độ cực $\nabla^2 = \frac{\partial^2}{\partial r^2} + \frac{1}{r}\frac{\partial}{\partial r} + \frac{1}{r^2}\frac{\partial^2}{\partial \theta^2}$.
  • Biểu diễn nghiệm tổng quát hàm ứng suất Airy cho bài toán đối xứng trục không phụ thuộc góc cực $\theta$: $$\phi(r) = A\ln r + Br^2\ln r + Cr^2 + D$$
  • Dẫn xuất lời giải giải tích cho các kết cấu chuyên dụng:
    • Bài toán thanh cong chịu uốn thuần túy (Lời giải Golovin, 1881).
    • Bài toán ống dày chịu áp lực trong, áp lực ngoài, ống có độ dôi và ống đặt trong môi trường đàn hồi (Bài toán Lamé, 1852).
    • Bài toán tập trung ứng suất quanh lỗ khoét tròn trong tấm chịu kéo.
    • Bài toán nêm phẳng chịu tải trọng tập trung, mômen và tải trọng phân bố; bài toán nền bán không gian chịu lực tập trung trên mặt biên.

Chương 4: Giải bài toán lý thuyết đàn hồi bằng phương pháp số

Trình bày cơ sở toán học và thuật toán của Phương pháp Phần tử Hữu hạn (FEM) theo mô hình chuyển vị:

  • Thiết lập sơ đồ rời rạc hóa kết cấu: Vectơ chuyển vị nút, vectơ tải trọng nút, hàm dạng (shape functions) và ma trận quan hệ biến dạng - chuyển vị $[B]$.
  • Xây dựng ma trận độ cứng phần tử $[k]^e$ thông qua nguyên lý cực tiểu thế năng toàn phần $\Phi^e$:
    • Phần tử tam giác bậc nhất (CST) cho bài toán phẳng 2D.
    • Phần tử thanh chịu kéo nén đúng tâm (1D).
    • Phần tử thanh chịu uốn phẳng và phần tử thanh không gian chịu lực phức tạp.
  • Kỹ thuật tính toán kết cấu tổng thể: Ma trận biến đổi toạ độ, kỹ thuật ghép nối ma trận độ cứng tổng thể $[K]$ và vectơ tải trọng ${F}$ bằng phương pháp số mã, thuật toán khử bậc tự do để xử lý điều kiện biên chuyển vị, kỹ thuật đánh số nút tối ưu nhằm thu hẹp bề rộng băng của ma trận $[K]$ phục vụ giải hệ phương trình tuyến tính trên máy tính số.

Kiến thức nền tảng được xây dựng

Giáo trình xây dựng hệ thống kiến thức cơ học dựa trên 3 nhóm giả thiết vật lý:

  1. Giả thiết môi trường: Vật liệu liên tục, đồng chất và đẳng hướng.
  2. Giả thiết cơ học: Vật liệu đàn hồi tuyến tính hoàn toàn và tuân theo định luật Hooke.
  3. Giả thiết hình học: Chuyển vị và biến dạng của vật thể là vô cùng bé.

Từ hệ giả thiết này, tính tuyến tính của các phương trình vi phân được thiết lập, cho phép áp dụng nguyên lý cộng tác dụng (chồng chất nghiệm) trong phân tích kết cấu.


Kỹ năng phát triển cho người học

Nhóm kỹ năng Nội dung cụ thể trong giáo trình Ứng dụng kỹ thuật
Phân tích giải tích Thiết lập hàm ứng suất Airy $\phi(x,y)$ và $\phi(r,\theta)$; giải hệ phương trình vi phân đạo hàm riêng thỏa mãn điều kiện biên. Tính toán ứng suất chính xác cho đập thủy lợi, dầm tường, ống áp lực.
Đánh giá sai số cơ học Vận dụng nguyên lý Saint-Venant để thay thế điều kiện biên lực tương đương; so sánh định lượng nghiệm Đàn hồi và nghiệm Sức bền vật liệu. Đánh giá phạm vi áp dụng của các công thức thanh cổ điển; kiểm tra hiện tượng tập trung ứng suất.
Mô hình hóa số trị Rời rạc hóa miền liên tục thành phần tử hữu hạn; lập ma trận độ cứng $[k]^e$, thiết lập ma trận toàn hệ $[K]$. Phát triển thuật toán và lập trình phần mềm tính toán kết cấu công trình.

Phương pháp giảng dạy và học tập

              TIẾN TRÌNH TIẾP CẬN SƯ PHẠM
┌────────────────────────────────────────────────────────┐
│ BƯỚC 1: Xây dựng lý thuyết đàn hồi không gian 3D       │
│ (Thiết lập hệ phương trình Navier, Cauchy, Hooke)      │
└──────────────────────────┬─────────────────────────────┘
                           │ Thu hẹp giả thiết
                           ▼
┌────────────────────────────────────────────────────────┐
│ BƯỚC 2: Thu gọn về bài toán phẳng 2D (Descartes & Cực) │
│ (Giải tích chính xác qua hàm song điều hoà Airy)       │
└──────────────────────────┬─────────────────────────────┘
                           │ Mở rộng biên phức tạp
                           ▼
┌────────────────────────────────────────────────────────┐
│ BƯỚC 3: Chuyển giao sang phương pháp phần tử hữu hạn   │
│ (Số hoá ma trận, lập trình giải hệ phương trình)       │
└────────────────────────────────────────────────────────┘

Giáo trình triển khai phương pháp sư phạm diễn dịch kết hợp quy nạp kỹ thuật:

  • Tiến trình logic: Xây dựng khung lý thuyết tổng quát cho không gian 3D tại Chương 1 $\rightarrow$ Rút gọn thành các bài toán phẳng 2D trong hệ toạ độ vuông góc và toạ độ cực tại Chương 2 và Chương 3 $\rightarrow$ Chuyển giao sang phương pháp phần tử hữu hạn tại Chương 4 để giải quyết các kết cấu có dạng hình học và điều kiện biên phức tạp mà nghiệm giải tích không thể xử lý.
  • Hệ thống bài tập mẫu và kiểm tra: Mỗi chương đều tích hợp các ví dụ giải tích chi tiết (như dầm conson, đập tam giác, ống dày) đi kèm chuỗi bài tập tự luyện (ví dụ các bài toán B-2-1 đến B-2-7 ở Chương 2). Hệ thống bài tập tập trung vào việc kiểm tra tính chất hàm Airy $\nabla^4 \phi = 0$, xác định tải trọng biên và tìm phân bố ứng suất cụ thể trên các cạnh của tấm chữ nhật hoặc tam giác.
  • Đánh giá và tự học: Người học có thể tự kiểm tra kết quả thông qua việc đối chiếu giữa nghiệm lý thuyết đàn hồi thu được với nghiệm sức bền vật liệu tương ứng (được phân tích chi tiết trong các bài toán dầm tường và thanh cong), từ đó nắm vững bản chất toán - cơ của từng phương pháp.

Điểm nổi bật và cập nhật

  1. Khắc phục hạn chế của cơ học công trình cổ điển: Khác với môn Sức bền vật liệu và Cơ học kết cấu chỉ giới hạn trên kết cấu thanh và đưa ra nghiệm gần đúng (do không thỏa mãn toàn bộ điều kiện biên), giáo trình cung cấp nghiệm thỏa mãn đầy đủ hệ phương trình vi phân và điều kiện biên, giải quyết được các bài toán ứng suất cục bộ, tập trung ứng suất quanh lỗ khoét và bài toán tiếp xúc.
  2. Cầu nối giữa giải tích và tính toán số: Tài liệu xây dựng hoàn chỉnh cầu nối từ các nguyên lý năng lượng cơ bản (công khả dĩ, cực tiểu thế năng Lagrange, Castigliano) trực tiếp sang phương pháp số (FEM), tạo điều kiện cho người học hiểu rõ bản chất giải thuật bên trong các phần mềm tính toán kết cấu hiện đại.
  3. Định hướng ứng dụng công trình Thủy lợi: Các nghiệm giải tích cổ điển (lời giải Lévy cho đập tam giác, dầm tường Filon - Ribiere, bài toán Lamé cho ống dày) được trình bày gắn liền với các công trình thực tế như đập bê tông dâng nước, tường chắn đất, hầm dẫn nước và kết cấu vỏ chịu áp.

Đối tượng sử dụng giáo trình

  • Sinh viên đại học: Sinh viên từ năm thứ 3 trở lên thuộc các ngành Kỹ thuật Xây dựng, Kỹ thuật Công trình Thủy, Kỹ thuật Công trình Biển, Kỹ thuật Cơ sở hạ tầng, Cơ kỹ thuật.
  • Học viên cao học và nghiên cứu sinh: Sử dụng làm tài liệu nghiên cứu chuyên đề về Cơ học vật rắn biến dạng, Lý thuyết đàn hồi ứng dụng và Cơ học tính toán.
  • Giảng viên và kỹ sư thiết kế: Làm giáo trình giảng dạy môn Phân tích ứng suất, Lý thuyết đàn hồi, hoặc tài liệu tra cứu các công thức nghiệm giải tích chuẩn mực để kiểm chứng (benchmark) mô hình số trong tính toán kết cấu công trình thủy lợi – thủy điện.

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

1. Giáo trình này phù hợp với đối tượng nào?

Giáo trình được thiết kế cho sinh viên năm 3, năm 4, học viên cao học các ngành Kỹ thuật Công trình (Thủy lợi, Xây dựng, Giao thông) và Cơ học, cũng như kỹ sư công trình cần tra cứu phương pháp phân tích ứng suất chính xác.

2. Cần kiến thức nền nào để học giáo trình này?

Người học cần hoàn thành các học phần: Giải tích toán học (phương trình vi phân đạo hàm riêng, chuỗi Fourier), Đại số tuyến tính (ma trận, định thức), Cơ học lý thuyết và Sức bền vật liệu cơ bản.

3. Điểm khác biệt lớn nhất so với giáo trình Sức bền vật liệu là gì?

Sức bền vật liệu chỉ nghiên cứu vật thể dạng thanh và sử dụng các giả thiết đơn giản hóa để tìm nghiệm gần đúng. Giáo trình Phân Tích Ứng Suất nghiên cứu vật thể 1D, 2D, 3D (thanh, tấm, vỏ, khối) bằng hệ phương trình vi phân đàn hồi chặt chẽ, thỏa mãn đầy đủ các điều kiện biên và giải quyết được bài toán ứng suất tập trung, ứng suất cục bộ.

4. Làm sao để tự học giáo trình này hiệu quả?

Người học nên bám sát quy trình: Nắm vững hệ phương trình vi phân cơ bản ở Chương 1 $\rightarrow$ Tự giải tay các bài toán phẳng mẫu ở Chương 2 và Chương 3 bằng hàm Airy $\rightarrow$ Thực hành lập trình thuật toán ma trận phần tử hữu hạn ở Chương 4 trên máy tính.

5. Tài liệu có cung cấp các bài toán kiểm chứng thực tế không?

Có. Giáo trình trình bày chi tiết các bài toán kỹ thuật tiêu biểu đã được các nhà khoa học công bố như: Lời giải Lévy (đập tam giác), lời giải Filon – Ribiere (dầm tường), lời giải Galovin (thanh cong) và lời giải Lamé (ống dày chịu áp).


Kết luận

Giáo trình Phân Tích Ứng Suất của tác giả Lý Trường Thành là tài liệu học thuật chuyên khảo mang tính hệ thống cao trong lĩnh vực Cơ học vật rắn biến dạng. Bằng việc kết hợp chặt chẽ giữa lý thuyết đàn hồi giải tích và phương pháp phần tử hữu hạn, cuốn sách cung cấp đầy đủ công cụ toán học và cơ học để khảo sát ứng suất - biến dạng trong các công trình kỹ thuật.

Lộ trình học tập đề xuất: Bắt đầu từ việc nắm vững hệ phương trình cân bằng – tương thích – vật lý, làm chủ kỹ thuật thiết lập hàm ứng suất Airy cho bài toán phẳng 2D, và hoàn thiện bằng kỹ năng số hóa ma trận phần tử hữu hạn phục vụ tính toán kết cấu trên máy tính. Người học có thể tham khảo thêm danh mục tài liệu chuyên sâu được dẫn tại trang 124 của giáo trình để mở rộng nghiên cứu.