Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình dân dụng và công nghiệp, kết cấu thép đóng vai trò then chốt nhờ cường độ chịu lực cao và khả năng vượt nhịp lớn. Dầm thép tiết diện chữ I cán nóng và tổ hợp hàn được ứng dụng phổ biến nhờ hiệu quả chịu uốn tối ưu quanh trục quán tính chính. Tuy nhiên, trong thực tế xây dựng, dầm thường xuyên chịu các hệ tải trọng đặt lệch tâm như tải trọng tường xây biên, dầm phụ liên kết một phía, dầm đỡ ban công hoặc ray cầu trục. Do thuộc dạng tiết diện thành mỏng hở, độ cứng chống xoắn của dầm chữ I rất thấp, dẫn đến hiện tượng xoắn kết hợp uốn phức tạp.

Nghiên cứu tập trung giải quyết bài toán phân tích trạng thái ứng suất và biến dạng của dầm thép chữ I chịu xoắn kiềm chế theo tiêu chuẩn của Viện Kết cấu Thép Hoa Kỳ (AISC). Mục tiêu cụ thể là hệ thống hóa cơ sở lý thuyết thanh thành mỏng, làm rõ sự tương tác giữa xoắn St. Venant và xoắn uốn vênh (warping torsion), đồng thời xây dựng công cụ phần mềm tự động hóa tính toán ứng suất cục bộ. Phạm vi nghiên cứu khảo sát dầm một nhịp dài 12 m, chịu tải trọng phân bố đều 10 kN/m đặt lệch tâm 0,1 m, làm việc trong giới hạn đàn hồi tuyến tính.

Ý nghĩa thực tiễn của công trình thể hiện qua việc lượng hóa chính xác sự gia tăng ứng suất nguy hiểm. Kết quả tính toán chỉ ra rằng hiệu ứng xoắn kiềm chế làm tăng ứng suất pháp tổng cộng tại gối tựa lên tới 98,3% (từ 5,9 kN/cm² lên 11,697 kN/cm²) và làm tăng ứng suất tiếp cục bộ tại bụng dầm lên tới 170% (từ 0,697 kN/cm² lên 1,884 kN/cm²). Việc bỏ qua xoắn kiềm chế trong thiết kế truyền thống có thể gây mất an toàn nghiêm trọng cho công trình.

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

Khung lý thuyết áp dụng

Cơ sở lý thuyết của luận văn dựa trên sự kết hợp giữa lý thuyết xoắn thuần túy của Saint-Venant và lý thuyết thanh thành mỏng tiết diện hở của V.Z. Vlasov. Khi một tiết diện chữ I chịu mô men xoắn ngoại lực, tổng mô men xoắn chống lại biến dạng bao gồm hai thành phần: mô men xoắn St. Venant sinh ra từ ứng suất cắt thuần túy và mô men xoắn uốn sinh ra do sự kiềm chế biến dạng vênh ngoài mặt phẳng của cánh dầm.

Phương trình vi phân tổng quát mô tả góc xoắn của thanh thành mỏng có dạng: $$GJ\theta' - EC_w\theta''' = T(z)$$

Trong đó:

  • $E$ là mô đun đàn hồi của thép ($210000\text{ N/mm}^2$).
  • $G$ là mô đun trượt cắt đàn hồi ($81000\text{ N/mm}^2$ với tỷ số $E/G = 2,6$).
  • $J$ là mô men quán tính xoắn St. Venant ($\approx \sum bt^3/3$).
  • $C_w$ là mô men quán tính quạt kháng vênh ($I_y h^2/4$).
  • $\theta, \theta', \theta'', \theta'''$ lần lượt là góc xoắn và các đạo hàm bậc một, bậc hai, bậc ba của góc xoắn theo chiều dài dầm $z$.

Hệ thống ứng suất trong dầm chữ I chịu xoắn kiềm chế bao gồm:

  1. Ứng suất tiếp do xoắn thuần túy: $\tau_t = G t \theta'$.
  2. Ứng suất pháp kháng vênh: $\sigma_w = -E \omega \theta''$ (với $\omega$ là tọa độ quạt).
  3. Ứng suất tiếp kháng vênh: $\tau_{ws} = -(E S_\omega / t) \theta'''$ (với $S_\omega$ là mô men tĩnh quạt).

Về quy chuẩn thiết kế, nghiên cứu đối chiếu hai phương pháp cốt lõi trong tiêu chuẩn AISC: phương pháp thiết kế theo hệ số tải trọng và độ bền (LRFD) với hệ số $\phi = 0,9$ cho trạng thái chảy dẻo và $\phi_c = 0,85$ cho mất ổn định; cùng phương pháp ứng suất cho phép (ASD) với hệ số an toàn $\Omega = 1,67$.

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

Nghiên cứu sử dụng phương pháp giải tích giải phương trình vi phân bậc cao kết hợp với kỹ thuật lập trình tự động hóa trên môi trường Visual Basic for Applications (VBA) tích hợp trong Microsoft Excel.

Quy trình chọn mẫu và phân tích bài toán:

  • Tập mẫu khảo sát dựa trên 12 sơ đồ chất tải và điều kiện biên điển hình theo chỉ dẫn thiết kế của AISC.
  • Nghiên cứu đi sâu phân tích trường hợp dầm chịu mô men xoắn phân bố đều với hai đầu ngàm kiềm chế vênh hoàn toàn (trường hợp 7).
  • Kích thước tiết diện khảo sát gồm dầm I500x300x20x20 mm, thép kết cấu CT42 có giới hạn chảy $F_y = 245\text{ N/mm}^2$ ($24,5\text{ kN/cm}^2$).

Lý do lựa chọn phương pháp giải tích kết hợp VBA xuất phát từ tính phức tạp của các hàm lượng giác hyperbol ($\sinh$, $\cosh$) trong nghiệm giải tích. Việc phát triển chương trình máy tính giúp số hóa toàn bộ quá trình xác định các đặc trưng hình học quạt, tự động vi phân tìm các giá trị cực trị của $\theta(z)$ và tổ hợp ma trận ứng suất đa trục tại các điểm nút nguy hiểm trên mặt cắt, loại bỏ hoàn toàn các sai số làm tròn khi tính toán thủ công.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phân tích số liệu từ mô hình tính toán dầm nhịp $L = 12\text{ m}$ dưới tác dụng của tải trọng uốn $w_u = 10\text{ kN/m}$ và mô men xoắn phân bố $t_u = 1\text{ kNm/m}$ mang lại các kết quả nổi bật sau:

  1. Gia tăng đột biến của ứng suất pháp tại gối tựa: Tại tiết diện gối ngàm ($z/L = 0$), ứng suất uốn thuần túy chỉ đạt $\sigma_b = 5,9\text{ kN/cm}^2$. Khi tính thêm thành phần ứng suất pháp do xoắn uốn $\sigma_w = 5,797\text{ kN/cm}^2$, ứng suất pháp tổng cộng tại mép ngoài bản cánh vọt lên mức $f_u = 11,697\text{ kN/cm}^2$, tương ứng với mức tăng 98,3%.
  2. Sự biến đổi phân bố ứng suất tiếp dọc trục dầm: Tại vị trí cách gối $z = 0,2L = 2,4\text{ m}$, ứng suất tiếp cắt uốn trên bụng dầm là $\tau_b = 0,697\text{ kN/cm}^2$. Khi cộng thêm ứng suất tiếp xoắn St. Venant $\tau_t = 1,187\text{ kN/cm}^2$, tổng ứng suất tiếp đạt $f_v = 1,884\text{ kN/cm}^2$, tăng 170% so với trường hợp không xét xoắn.
  3. Quy luật phân bố biến dạng xoắn: Góc xoắn đạt giá trị cực đại tại giữa nhịp ($z/L = 0,5$) với $\theta_{max} = 0,0285\text{ rad}$ ($1,63^\circ$), trong khi tại gối ngàm góc xoắn và biến dạng vênh bị kiềm chế triệt để ($\theta = 0, \theta' = 0$).
  4. Đánh giá an toàn theo tiêu chuẩn AISC/LRFD: Đối với thép CT42 ($F_y = 24,5\text{ kN/cm}^2$), cường độ thiết kế danh nghĩa đạt $\phi F_y = 0,9 \times 24,5 = 22,05\text{ kN/cm}^2$. Ứng suất pháp lớn nhất $11,697\text{ kN/cm}^2$ chiếm khoảng 53% khả năng chịu lực, đảm bảo độ bền tổng thể nhưng cho thấy mức độ dự trữ chịu tải giảm đáng kể so với tính toán uốn phẳng đơn thuần.

Thảo luận kết quả

Nguyên nhân cốt lõi khiến ứng suất tăng vọt nằm ở đặc tính hình học của tiết diện hở. Mô men quán tính xoắn St. Venant của dầm chỉ đạt $J = 283\text{ cm}^4$, rất nhỏ so với mô men quán tính kháng uốn $I_{xx} = 50823\text{ cm}^4$. Khi chịu tải xoắn, độ cứng chống xoắn St. Venant ($GJ$) không đủ khả năng triệt tiêu nội lực, buộc kết cấu phải chuyển hóa năng lượng biến dạng sang cơ chế xoắn uốn vênh ($EC_w = 2,1 \times 10^4 \times 5,2 \times 10^6\text{ kN}\cdot\text{cm}^4$).

Đặc trưng tỷ lệ xoắn uốn $L/a = 5,49$ (với $a = 218,3\text{ cm}$) cho thấy hiệu ứng xoắn kiềm chế chiếm ưu thế tuyệt đối trong phạm vi khoảng cách $3a \approx 6,55\text{ m}$ tính từ vị trí gối ngàm. Ứng suất pháp vênh tác động kéo - nén ngược chiều trên hai mép của cùng một bản cánh, làm thay đổi hoàn toàn quy luật phân bố ứng suất cổ điển của Bernoulli.

Dữ liệu phân tích có thể được mô hình hóa trực quan qua bảng tổng hợp ứng suất 4 góc cánh dầm và biểu đồ phân bố ứng suất tiếp 3D dọc theo chu tuyến tiết diện. So sánh với các nghiên cứu giải tích quốc tế cho thấy phương pháp lập trình tự động hóa VBA cho kết quả trùng khớp hoàn toàn với dung sai dưới 0,5%, chứng minh độ tin cậy vượt trội trong việc thay thế các bảng tra thủ công phức tạp.

Đề xuất và khuyến nghị

Nhằm nâng cao chất lượng thiết kế kết cấu thép chịu tải lệch tâm và phòng ngừa sự cố công trình, các đề xuất cụ thể được đưa ra như sau:

  1. Chuẩn hóa quy trình tính toán xoắn kiềm chế trong tiêu chuẩn Việt Nam:

    • Hành động: Bộ Xây dựng và các viện nghiên cứu chuyên ngành cần bổ sung chương mục tính toán cấu kiện thanh thành mỏng chịu uốn xoắn đồng thời dựa trên tiêu chuẩn AISC 360 và AISC Design Guide 9 vào hệ thống tiêu chuẩn thiết kế kết cấu thép quốc gia.
    • Chỉ số mục tiêu: Đạt 100% tài liệu hướng dẫn kỹ thuật có ví dụ tính toán chi tiết về dầm thép chữ I chịu xoắn.
    • Thời gian thực hiện: Giai đoạn 2026–2027.
    • Chủ thể thực hiện: Ban soạn thảo quy chuẩn Bộ Xây dựng.

  2. Áp dụng phần mềm tự động hóa tính toán trong các đơn vị tư vấn:

    • Hành động: Các doanh nghiệp tư vấn thiết kế kết cấu cần triển khai tích hợp các module lập trình giải thuật vi phân góc xoắn bằng VBA hoặc Python vào quy trình tính toán nội bộ.
    • Chỉ số mục tiêu: Rút ngắn 85% thời gian lập bảng tính thủ công và loại bỏ 100% sai sót tính toán do tra biểu đồ.
    • Thời gian thực hiện: Quý 3 năm 2026.
    • Chủ thể thực hiện: Trưởng phòng kết cấu và kỹ sư thiết kế tại các công ty tư vấn xây dựng.

  3. Gia cường cấu tạo kháng vênh tại các nút liên kết nguy hiểm:

    • Hành động: Bố trí hệ sườn tăng cường bản cánh, bản chặn đầu ngàm (end plates) hoặc sườn tăng cường dạng hộp tại vị trí liên kết ngàm và các điểm đặt tải trọng tập trung lệch tâm lớn.
    • Chỉ số mục tiêu: Giảm ít nhất 40% góc xoắn cục bộ và triệt tiêu tới 60% ứng suất pháp vênh tại mép cánh dầm.
    • Thời gian thực hiện: Áp dụng trực tiếp trong giai đoạn thiết kế bản vẽ thi công các dự án mới.
    • Chủ thể thực hiện: Kỹ sư thiết kế và kỹ sư giám sát công trình.

  4. Cập nhật chương trình đào tạo chuyên sâu về cơ học thanh thành mỏng:

    • Hành động: Các trường đại học khối kỹ thuật xây dựng cần tăng cường thời lượng giảng dạy lý thuyết Vlasov và bài toán xoắn kiềm chế trong chương trình đào tạo kỹ sư và thạc sĩ.
    • Chỉ số mục tiêu: 100% học viên tốt nghiệp nắm vững bản chất bimoment và khả năng tổ hợp ứng suất uốn xoắn.
    • Thời gian thực hiện: Niên khóa 2026–2028.
    • Chủ thể thực hiện: Khoa Công trình và Bộ môn Kết cấu Thép - Gỗ tại các trường đại học.

Đối tượng nên tham khảo luận văn

Tài liệu luận văn mang lại giá trị học thuật và thực tiễn sâu sắc cho bốn nhóm đối tượng chính:

  1. Kỹ sư thiết kế kết cấu công trình:

    • Lợi ích: Nắm vững quy trình kiểm tra bền dầm chữ I chịu uốn - xoắn tổ hợp theo AISC LRFD và ASD.
    • Tình huống sử dụng: Thiết kế dầm đỡ tường biên nhà cao tầng khung thép, dầm đỡ ban công lệch tâm hoặc dầm cầu trục nhà xưởng công nghiệp.

  2. Học viên cao học và nghiên cứu sinh ngành Xây dựng:

    • Lợi ích: Tiếp cận hệ thống lý thuyết thanh thành mỏng Vlasov một cách tường minh, kèm theo các thuật toán vi phân chuẩn hóa.
    • Tình huống sử dụng: Làm tài liệu tham khảo nền tảng cho các đề tài nghiên cứu về ổn định phi tuyến, dao động xoắn hoặc phân tích phần tử hữu hạn nâng cao.

  3. Giảng viên và sinh viên các trường đại học kỹ thuật:

    • Lợi ích: Sử dụng làm giáo trình phụ trợ, bài giảng chuyên đề kết cấu thép và cơ học vật rắn biến dạng.
    • Tình huống sử dụng: Giảng dạy đồ án kết cấu thép và hướng dẫn sinh viên làm đồ án tốt nghiệp kỹ sư xây dựng.

  4. Kỹ sư thẩm tra và quản lý chất lượng dự án:

    • Lợi ích: Có công cụ đối soát nhanh chóng, chính xác để phát hiện các trường hợp dầm biên bị thiếu khả năng chịu tải do bỏ qua ứng suất xoắn.
    • Tình huống sử dụng: Thẩm tra hồ sơ thiết kế kỹ thuật và đánh giá an toàn chịu lực của các công trình công nghiệp quy mô lớn.

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

Tại sao dầm thép chữ I lại có khả năng chống xoắn thuần túy rất kém?

Dầm thép chữ I là tiết diện thành mỏng hở. Mô men quán tính xoắn $J$ của tiết diện hở tỷ lệ với bậc ba chiều dày bản thép ($\sum bt^3/3$), dẫn đến giá trị $J = 283\text{ cm}^4$ rất nhỏ so với mô men kháng uốn $I_{xx} = 50823\text{ cm}^4$. Khi chịu tải lệch tâm, cánh dầm dễ bị biến dạng trượt và vênh ngoài mặt phẳng, đòi hỏi phải xét thêm hiệu ứng xoắn kiềm chế.

Sự khác biệt cơ bản giữa xoắn St. Venant và xoắn kiềm chế là gì?

Xoắn St. Venant xảy ra khi mặt cắt ngang được tự do vênh ngoài mặt phẳng, chỉ sinh ra ứng suất tiếp thuần túy $\tau_t$. Ngược lại, xoắn kiềm chế xuất hiện khi các liên kết gối tựa hoặc cấu tạo ngăn cản sự vênh, làm phát sinh thêm ứng suất pháp vênh $\sigma_w = -E\omega\theta''$ và ứng suất tiếp vênh $\tau_{ws} = -(ES_\omega/t)\theta'''$.

Tiêu chuẩn AISC kiểm soát ứng suất xoắn kết hợp theo các trạng thái giới hạn nào?

Theo AISC, cấu kiện được kiểm tra theo hai phương pháp chính: LRFD với điều kiện ứng suất pháp tổ hợp $f_{un} \le \phi F_y = 0,9 F_y$ và ứng suất tiếp $f_{uv} \le \phi (0,6 F_y)$; hoặc ASD với điều kiện $f_n \le 0,6 F_y$ và $f_v \le 0,4 F_y$, đồng thời kiểm toán mất ổn định uốn - xoắn tổng thể.

Tại sao ứng suất pháp lớn nhất lại xuất hiện tại gối ngàm thay vì tại giữa nhịp?

Tại gối ngàm, biến dạng vênh của dầm bị kiềm chế hoàn toàn, khiến đạo hàm bậc hai của góc xoắn $\theta''$ đạt giá trị cực đại ($0,64\text{ rad/cm}^2$). Sự kiềm chế biến dạng này tạo ra mô men xoắn uốn vênh lớn nhất, làm phát sinh ứng suất pháp vênh $\sigma_w = 5,797\text{ kN/cm}^2$ cộng hưởng trực tiếp với ứng suất uốn lớn nhất $\sigma_b = 5,9\text{ kN/cm}^2$.

Ứng dụng ngôn ngữ Visual Basic for Applications (VBA) mang lại lợi ích gì cho việc tính toán?

VBA giúp giải tự động phương trình vi phân góc xoắn chứa các hàm hyperbol phức tạp dọc theo 12 trường hợp biên của AISC. Chương trình tự động hóa xác định các đặc trưng hình học quạt ($C_w, S_\omega$), tính đạo hàm 4 bậc của góc xoắn và lập biểu đồ ứng suất tức thì, giúp tiết kiệm hơn 85% thời gian tính toán.

Kết luận

  • Luận văn đã xây dựng hoàn chỉnh phương pháp tính toán dầm thép chữ I chịu xoắn uốn đồng thời dựa trên sự kết hợp giữa lý thuyết thanh thành mỏng Vlasov và tiêu chuẩn AISC.
  • Kết quả kiểm chứng trên dầm nhịp 12 m cho thấy ứng suất pháp tổng cộng tăng 98,3% (đạt 11,697 kN/cm²) và ứng suất tiếp tăng 170% (đạt 1,884 kN/cm²) dưới tác dụng của tải lệch tâm.
  • Chương trình tính toán phát triển trên nền tảng Visual Basic/Excel đã số hóa thành công các phương trình giải tích phức tạp, cung cấp công cụ tra cứu chính xác cho kỹ sư thực hành.
  • Nghiên cứu mở ra hướng phát triển mở rộng cho các bài toán phân tích phi tuyến hình học, dầm cong chịu xoắn và dầm thép tiết diện thay đổi trong giai đoạn tiếp theo.
  • Hãy chủ động ứng dụng quy trình tính toán xoắn kiềm chế theo AISC và tích hợp các công cụ tự động hóa vào dự án thực tế ngay hôm nay để tối ưu hóa an toàn kết cấu công trình.