I. Hướng dẫn toàn diện tính toán bền cho mái vòm composite
Luận văn thạc sĩ của Nguyễn Hải Bằng (2012) trình bày một phương pháp tính toán bền cho mái vòm làm hoàn toàn bằng vật liệu composite có dạng lượn sóng dưới tác dụng của tải trọng gió. Nghiên cứu này có ý nghĩa thực tiễn to lớn, đặc biệt trong bối cảnh các công trình xây dựng ven biển và trong môi trường ăn mòn đòi hỏi vật liệu có độ bền cao và khả năng chống oxy hóa. Vật liệu composite cốt sợi (FRP), như composite cốt sợi thủy tinh (GFRP) hay cốt sợi carbon (CFRP), nổi lên như một giải pháp thay thế ưu việt cho các vật liệu truyền thống như bê tông cốt thép hay kim loại. Các kết cấu vỏ mỏng dạng vòm composite không chỉ nhẹ, dễ thi công mà còn có khả năng chịu lực tốt. Tuy nhiên, việc phân tích kết cấu của chúng, đặc biệt với hình dạng lượn sóng phức tạp, là một thách thức. Nghiên cứu tập trung vào việc thiết lập một mô hình toán học đáng tin cậy để dự đoán trạng thái ứng suất và biến dạng của kết cấu, đảm bảo ổn định kết cấu khi chịu các tải trọng tác dụng khắc nghiệt. Mục tiêu chính là xây dựng một công cụ tính toán thuận tiện, cho phép các kỹ sư xác định độ dày tối thiểu của vỏ mái vòm hoặc giới hạn sức gió mà kết cấu có thể chịu đựng, từ đó tối ưu hóa thiết kế và đảm bảo an toàn cho công trình. Luận văn đã thành công trong việc xây dựng mô hình lý thuyết, áp dụng các phương pháp giải tiên tiến và phát triển phần mềm tính toán chuyên dụng.
1.1. Giới thiệu về vật liệu composite cốt sợi FRP trong xây dựng
Việc sử dụng vật liệu composite cốt sợi (FRP) trong các công trình dân dụng và quốc phòng đang trở thành một xu hướng tất yếu. Các vật liệu này, chẳng hạn như composite nền epoxy cốt sợi thủy tinh, sở hữu những đặc tính vượt trội. Chúng có độ bền vật liệu cao, khối lượng riêng thấp, và đặc biệt là khả năng kháng hóa chất và chống ăn mòn tuyệt vời. Điều này làm cho chúng trở thành lựa chọn lý tưởng cho các công trình trong môi trường biển, hải đảo hoặc các khu công nghiệp có môi trường axit. Luận văn nhấn mạnh, "khi thiết kế và thi công các công trình trong môi trường biển có muối mặn cần một loại vật liệu và kết cấu không có cốt thép". Cơ học vật liệu composite cho thấy, với cấu trúc vật liệu dị hướng, chúng ta có thể tùy chỉnh các đặc tính cơ học theo các phương chịu lực chính, mang lại hiệu quả kinh tế và kỹ thuật cao.
1.2. Tầm quan trọng của việc phân tích kết cấu mái vòm lượn sóng
Mái vòm dạng tấm lượn sóng được sử dụng để tăng cường độ cứng và khả năng chịu lực của kết cấu mà không làm tăng đáng kể trọng lượng. Tuy nhiên, hình dạng hình học phức tạp này đòi hỏi một phương pháp phân tích kết cấu chi tiết và chính xác. Việc tính toán độ bền và ổn định của các kết cấu không gian này còn chưa được chuẩn hóa đầy đủ, đặc biệt khi chúng được làm hoàn toàn từ vật liệu composite. Một phân tích kết cấu chính xác sẽ giúp dự đoán các vị trí tập trung ứng suất, độ võng tối đa và nguy cơ mất ổn định, từ đó đưa ra các giải pháp thiết kế an toàn và tiết kiệm vật liệu. Nghiên cứu này giải quyết trực tiếp khoảng trống đó bằng cách cung cấp một quy trình tính toán hoàn chỉnh.
II. Thách thức lớn trong tính toán ổn định kết cấu mái vòm
Việc tính toán bền cho mái vòm composite lượn sóng phải đối mặt với nhiều thách thức kỹ thuật. Đầu tiên, đặc tính vật liệu dị hướng của composite làm cho các phương trình cơ học trở nên phức tạp hơn nhiều so với vật liệu đẳng hướng truyền thống. Các hằng số đàn hồi thay đổi theo phương, đòi hỏi phải áp dụng lý thuyết đàn hồi cho vật liệu trực hướng. Thứ hai, hình học lượn sóng của vỏ làm tăng độ khó trong việc thiết lập các phương trình cân bằng và điều kiện biên. Các lý thuyết vỏ mỏng cổ điển cần được điều chỉnh để phản ánh đúng sự thay đổi về độ cứng theo các phương khác nhau. Thách thức lớn nhất là kết hợp chính xác ảnh hưởng của hình học lượn sóng và tính dị hướng của vật liệu. Ngoài ra, việc xác định tải trọng tác dụng của gió lên một bề mặt cong và không đồng nhất cũng là một bài toán phức tạp. Luận văn đã chỉ ra sự cần thiết phải tuân thủ tiêu chuẩn thiết kế TCVN 2737-1995 để phân vùng áp lực gió một cách hợp lý. Việc giải quyết đồng thời các vấn đề này đòi hỏi một phương pháp luận chặt chẽ, kết hợp giữa lý thuyết cơ học vững chắc và các công cụ giải số hiệu quả để đảm bảo độ tin cậy kết cấu.
2.1. Phân tích ảnh hưởng của tải trọng tác dụng từ môi trường
Tải trọng gió là một trong những tải trọng tác dụng nguy hiểm và khó lường nhất đối với các kết cấu mái che có khẩu độ lớn như mái vòm. Theo tiêu chuẩn thiết kế TCVN 2737–1995, áp lực gió không phân bố đều trên bề mặt vòm mà thay đổi tùy thuộc vào vị trí và hình dạng khí động học của công trình. Luận văn đã áp dụng tiêu chuẩn này để chia mái vòm thành các vùng chịu áp lực hút và đẩy khác nhau, tạo ra một sơ đồ tải trọng phức tạp. Việc mô hình hóa chính xác sơ đồ tải trọng này là bước đầu tiên và quan trọng nhất để đảm bảo kết quả phân tích kết cấu phản ánh đúng thực tế, từ đó đánh giá chính xác ổn định kết cấu.
2.2. Hạn chế của lý thuyết cổ điển và sự phức tạp của vật liệu
Lý thuyết màng (lý thuyết phi mô men) thường được áp dụng cho các kết cấu vỏ mỏng, tuy nhiên nó không đủ để phân tích trạng thái uốn cục bộ tại các vị trí liên kết hoặc nơi có sự thay đổi đột ngột về hình học. Đối với vỏ lượn sóng, sự thay đổi độ cứng do hình dạng yêu cầu phải xem xét cả trạng thái uốn và kéo-nén. Hơn nữa, cơ học vật liệu composite với các đặc tính như vật liệu dị hướng đòi hỏi các hệ thức quan hệ giữa ứng suất và biến dạng phức tạp hơn định luật Hooke cho vật liệu đẳng hướng. Luận văn đã vượt qua những hạn chế này bằng cách sử dụng một lý thuyết vỏ cải tiến và các phương pháp giải số mạnh mẽ, thay vì chỉ dựa vào các công thức giải tích đơn giản.
III. Phương pháp Seydel Tối ưu tính toán mái vòm lượn sóng
Để giải quyết sự phức tạp của hình học lượn sóng, luận văn của Nguyễn Hải Bằng đã áp dụng một cách sáng tạo phương pháp Seydel. Đây là một phương pháp kinh điển và hiệu quả để đưa bài toán tính toán một kết cấu vỏ mỏng có gân hoặc lượn sóng về bài toán tính toán một vỏ trụ trơn tương đương. Cốt lõi của phương pháp này là thay thế các thông số hình học phức tạp bằng cách điều chỉnh các hằng số đàn hồi của vật liệu. Vỏ lượn sóng ban đầu được xem như một vỏ trơn làm bằng vật liệu dị hướng với các mô đun đàn hồi hiệu quả mới. Các mô đun này (E'₁, E'₂, G') được tính toán dựa trên các thông số của sóng như chu kỳ, biên độ (H) và chiều dài cung. Cách tiếp cận này giúp đơn giản hóa đáng kể các phương trình vi phân đạo hàm riêng mô tả trạng thái cân bằng của vỏ. Thay vì phải giải quyết một bài toán với hình học biến đổi liên tục, người nghiên cứu có thể làm việc với một hệ phương trình trên một miền hình học đơn giản hơn. Luận văn trích dẫn: "công trình [5] đã sử dụng lý thuyết của Seydel để đưa dạng vỏ có dạng lượn sóng về tương đương với vỏ ở dạng không lượn sóng với các hằng số đàn hồi bị thay đổi".
3.1. Nguyên lý chuyển đổi kết cấu vỏ mỏng lượn sóng về vỏ tương đương
Nguyên lý cơ bản của phương pháp Seydel là bảo toàn độ cứng uốn và độ cứng kéo-nén của kết cấu vỏ mỏng trong các phương chính. Hình dạng lượn sóng làm tăng độ cứng của vỏ theo phương vuông góc với đường sinh một cách đáng kể, trong khi độ cứng dọc theo đường sinh gần như không đổi. Phương pháp Seydel quy đổi sự gia tăng độ cứng hình học này thành sự gia tăng của mô đun đàn hồi vật liệu. Về bản chất, một vỏ composite lượn sóng được mô hình hóa thành một vỏ composite trơn nhưng có tính vật liệu dị hướng rõ rệt hơn, với các mô đun đàn hồi mới được tính toán theo công thức (2-18) của luận văn. Điều này cho phép áp dụng các lý thuyết tính toán vỏ trụ đã được phát triển trước đó.
3.2. Xác định các hằng số đàn hồi mới cho tấm lượn sóng
Việc xác định các mô đun đàn hồi hiệu quả là bước then chốt trong phương pháp Seydel. Luận văn đã trình bày chi tiết công thức tính toán các độ cứng trụ mới D'₁, D'₂, D'₃ và mô đun trượt G' dựa trên các thông số của tấm lượn sóng như chu kỳ sóng (λ) và chiều cao sóng (H). Các công thức này được xây dựng từ việc phân tích mô men quán tính của mặt cắt ngang lượn sóng. Sau khi có các độ cứng mới, các mô đun đàn hồi hiệu quả E'₁ và E'₂ được suy ra. Quá trình này biến đổi một bài toán phức tạp về hình học thành một bài toán dễ giải quyết hơn về cơ học vật liệu composite dị hướng, mở đường cho việc áp dụng các phương pháp giải số tiếp theo.
IV. Bí quyết giải bài toán bền bằng phương pháp Ritz Galerkin
Sau khi đã đơn giản hóa bài toán hình học bằng phương pháp Seydel, luận văn tiếp tục giải hệ phương trình vi phân bằng cách kết hợp phương pháp Ritz và phương pháp Bubnov-Galerkin. Đây là những kỹ thuật biến phân mạnh mẽ, thường được sử dụng trong phân tích kết cấu và là nền tảng của phương pháp phần tử hữu hạn (FEM). Phương pháp Ritz yêu cầu tìm nghiệm gần đúng dưới dạng một chuỗi các hàm được chọn trước (hàm cơ sở) thỏa mãn các điều kiện biên của bài toán. Luận văn đã đề xuất các dạng hàm nghiệm cho độ võng (w) và hàm ứng suất (φ) dưới dạng chuỗi hàm lượng giác và đa thức. Việc lựa chọn các hàm cơ sở này là cực kỳ quan trọng, chúng phải "thỏa mãn tất cả các điều kiện biên, trong đó điều kiện biên tại hai cạnh x = ± a tự do là khó nhất". Sau đó, phương pháp Bubnov-Galerkin được áp dụng để xác định các hệ số chưa biết trong chuỗi hàm. Nguyên lý của phương pháp này là làm cho sai số (phần dư) của phương trình vi phân trực giao với tất cả các hàm cơ sở. Quá trình này dẫn đến một hệ phương trình đại số tuyến tính, từ đó có thể giải để tìm các hệ số và xác định được trường ứng suất và biến dạng trong toàn bộ mái vòm.
4.1. Thiết lập hệ phương trình cho ứng suất và biến dạng
Bài toán được mô tả bởi một hệ hai phương trình vi phân đạo hàm riêng bậc bốn, một phương trình tương thích biến dạng (liên quan đến hàm ứng suất φ) và một phương trình cân bằng (liên quan đến độ võng w). Việc đưa vào hàm ứng suất φ giúp thỏa mãn tự động hai phương trình cân bằng trong mặt phẳng. Do đó, bài toán quy về việc tìm hai hàm φ và w. Hệ phương trình này phản ánh mối liên hệ phức tạp giữa các thành phần nội lực, mô men, ứng suất và biến dạng của kết cấu vỏ mỏng dưới tác dụng của tải trọng ngoài.
4.2. Tìm nghiệm qua chuỗi hàm và điều kiện biên của mái vòm
Thành công của phương pháp Ritz phụ thuộc rất nhiều vào việc chọn dạng nghiệm. Luận văn đã chọn các chuỗi hàm (3-22) một cách khéo léo để thỏa mãn các điều kiện biên phức tạp: ngàm cứng tại hai biên dọc đường sinh (y = ±b) và tự do tại hai đầu vòm (x = ±a). Các điều kiện ngàm yêu cầu chuyển vị và góc xoay bằng không. Các điều kiện biên tự do yêu cầu nội lực và mô men pháp tuyến bằng không. Bằng cách áp dụng phương pháp Bubnov-Galerkin, một hệ phương trình đại số tuyến tính được thiết lập. Việc giải hệ này, dù phức tạp, có thể thực hiện được bằng các công cụ tính toán số, chẳng hạn như phần mềm do chính tác giả xây dựng.
V. Kết quả mô phỏng số và phân tích độ bền vật liệu mái vòm
Chương 4 của luận văn tập trung vào việc trình bày các kết quả tính toán cụ thể, sử dụng một phần mềm mô phỏng được tác giả tự xây dựng trên nền tảng Visual Basic. Phần mềm này cho phép nhập các thông số hình học (chiều dài L, bán kính R, độ dày h, biên độ sóng H) và thông số vật liệu (E, ν) để tiến hành mô phỏng số. Kết quả đầu ra bao gồm các bảng giá trị và biểu đồ trực quan về độ võng, mô men uốn (My), và các thành phần ứng suất (σx, σy, σ₁₂) tại một lưới các điểm trên bề mặt mái vòm. Một trong những kết quả quan trọng nhất là việc xác định cường độ ứng suất tương đương (σu) theo thuyết bền năng lượng biến đổi hình dạng. Giá trị này được so sánh với giới hạn bền của vật liệu để đánh giá độ tin cậy kết cấu. Nghiên cứu chỉ ra rằng "vị trí các giá trị này đạt Max hoặc min là không thay đổi" khi chỉ có cường độ gió hoặc độ dày h thay đổi. Cụ thể, độ võng lớn nhất xuất hiện ở gần đỉnh vòm, phía sau hướng gió. Mô men uốn My đạt cực đại tại biên ngàm, trong khi cường độ ứng suất lớn nhất lại nằm ở một vị trí khác. Những phát hiện này cung cấp cái nhìn sâu sắc về hành vi chịu lực của kết cấu.
5.1. Phân tích cường độ ứng suất và mô men uốn qua phần mềm
Phần mềm tính toán đóng vai trò như một công cụ mô phỏng số hiệu quả, tương tự như các phần mềm mô phỏng thương mại (ANSYS, ABAQUS) nhưng được chuyên biệt hóa cho bài toán này. Kết quả từ Bảng 4.6 cho thấy, với độ dày h = 2 cm và các thông số cho trước, cường độ ứng suất cực đại vẫn nằm trong giới hạn đàn hồi cho phép của vật liệu. Biểu đồ mô men uốn (Hình 4.2) chỉ rõ sự tập trung mô men tại các đường biên ngàm, đây là thông tin quan trọng cho việc thiết kế các liên kết. Việc phân tích kết cấu bằng công cụ này giúp người thiết kế nhanh chóng đánh giá được mức độ an toàn của mái vòm.
5.2. Mối quan hệ giữa độ dày vỏ tải trọng và độ tin cậy kết cấu
Mục tiêu cốt lõi của đề tài là thiết lập mối quan hệ giữa các tham số thiết kế. Kết quả tính toán cho phép xây dựng các biểu đồ hoặc bảng tra, thể hiện sự phụ thuộc của ứng suất lớn nhất vào độ dày vỏ (h) và áp lực gió (P₀). Từ đó, kỹ sư có thể trả lời hai câu hỏi thiết kế quan trọng: (1) Với một cấp gió bão cho trước, độ dày vỏ tối thiểu là bao nhiêu để đảm bảo an toàn? (2) Với một mái vòm có độ dày cho trước, nó có thể chịu được sức gió tối đa là bao nhiêu? Việc trả lời được hai câu hỏi này chính là cơ sở để đảm bảo độ tin cậy kết cấu và tối ưu hóa hình học cũng như chi phí vật liệu.
VI. Ứng dụng thực tiễn và tương lai của kết cấu không gian
Kết quả nghiên cứu của luận văn "Tính toán bền cho mái vòm làm hoàn toàn bằng vật liệu Composite có dạng lượn sóng" mang lại giá trị ứng dụng thực tiễn cao. Phương pháp và công cụ tính toán được phát triển có thể áp dụng trực tiếp vào việc thiết kế các công trình kết cấu không gian như nhà chờ, mái che sân vận động, kho chứa hàng tại các khu vực ven biển, hải đảo hoặc môi trường công nghiệp khắc nghiệt. Việc sử dụng vật liệu composite giúp giải quyết triệt để vấn đề ăn mòn, giảm chi phí bảo trì và kéo dài tuổi thọ công trình. Phần mềm tính toán đi kèm là một công cụ hữu ích cho các kỹ sư, giúp rút ngắn thời gian thiết kế và đưa ra các quyết định dựa trên cơ sở khoa học vững chắc. Hướng phát triển trong tương lai của lĩnh vực này có thể bao gồm việc tối ưu hóa hình học của dạng lượn sóng (ví dụ: hình dạng tối ưu không phải sin mà là dạng khác) để tăng tối đa độ cứng. Ngoài ra, có thể mở rộng mô hình để xem xét các yếu tố phức tạp hơn như phân tích dao động của kết cấu dưới tải trọng gió động, hoặc phân tích bài toán phi tuyến vật liệu và phi tuyến hình học khi kết cấu chịu tải trọng lớn.
6.1. Khả năng ứng dụng thực tiễn trong các công trình đặc thù
Các công trình quốc phòng và dân sinh tại các vùng biển đảo luôn yêu cầu vật liệu và kết cấu đặc biệt. Mái vòm composite lượn sóng là một giải pháp lý tưởng cho các trạm radar, nhà trú ẩn, kho tàng... nhờ khả năng chống chịu môi trường muối mặn và gió bão. Kết quả của luận văn cung cấp một cơ sở lý thuyết tin cậy để các đơn vị thiết kế và thi công có thể tự tin áp dụng loại vật liệu và kết cấu tiên tiến này, đảm bảo độ tin cậy kết cấu trong những điều kiện vận hành khó khăn nhất. Việc tính toán được độ dày cần thiết giúp tiết kiệm vật liệu mà vẫn đảm bảo an toàn tuyệt đối.
6.2. Hướng phát triển cho việc tối ưu hóa hình học kết cấu vòm
Nghiên cứu này đã mở ra nhiều hướng phát triển tiềm năng. Một trong số đó là bài toán tối ưu hóa hình học. Thay vì sử dụng dạng sóng hình sin cố định, các thuật toán tối ưu hóa có thể được áp dụng để tìm ra biên dạng sóng tối ưu nhằm đạt được độ cứng lớn nhất với cùng một lượng vật liệu. Ngoài ra, việc kết hợp phương pháp phần tử hữu hạn (FEM) với mô hình giải tích đã xây dựng có thể giúp kiểm chứng và hiệu chỉnh kết quả, đồng thời cho phép phân tích các dạng mái vòm có hình dạng tổng thể phức tạp hơn. Nghiên cứu về phân tích dao động và ổn định động lực học cũng là một hướng đi cần thiết để hiểu đầy đủ hành vi của các kết cấu không gian nhẹ dưới tác động của gió giật.
