Một số vấn đề của lý thuyết kết cấu vỏ mỏng đàn hồi

Quá trình phát triển lý thuyết vỏ gắn liền với sự tiến hóa của toán học giải tích và cơ học môi trường liên tục. Giai đoạn sơ khai dựa trên các định đề hình học và vật lý uốn tấm phẳng của nhà cơ học Germain. Việc giải trực tiếp hệ phương trình vi phân đạo hàm riêng với hệ số biến thiên từng gặp nhiều khó khăn toán học. Để áp dụng vào thực tế kỹ thuật, các nhà khoa học đã xây dựng nhiều mô hình đơn giản hóa dựa trên thực nghiệm. Nhiều nhánh lý thuyết chuyên biệt đã ra đời như lý thuyết hiệu ứng biên, lý thuyết bán mômen của vỏ trụ, lý thuyết vỏ thoải và lý thuyết phi mômen. Các giả thiết giản lược giúp phân lập trạng thái ứng lực chính, biến các phương trình phức tạp thành các hệ giải được. Bước chuyển mình này đã thúc đẩy ứng dụng kết cấu mái vòm rộng rãi trong các công trình kiến trúc quy mô lớn.

Kết cấu vỏ mỏng đàn hồi có phạm vi ứng dụng rộng rãi trong nhiều ngành công nghiệp kỹ thuật cao. Trong xây dựng dân dụng và công nghiệp, vỏ mỏng được sử dụng làm mái vòm sân vận động, nhà triển lãm, tháp làm mát và silô chứa vật liệu rời. Khả năng bao che không gian lớn mà không cần cột chống trung gian là ưu điểm nổi bật. Trong ngành công nghiệp dầu khí và hóa chất, vỏ mỏng định hình nên các bồn chứa áp lực cao, đường ống dẫn khí và bình chứa khí hóa lỏng. Trong ngành kỹ thuật hàng không và vũ trụ, cấu trúc vỏ mỏng hợp kim hoặc composite tạo nên thân máy bay, vỏ tên lửa và tàu không gian. Ngành công nghiệp hàng hải ứng dụng vỏ mỏng để thiết kế thân tàu thủy và vỏ tàu ngầm chịu áp lực nước sâu. Cấu trúc mỏng nhẹ giúp tối ưu trọng lượng bản thân và tăng cường khả năng chịu tải tổng thể.

Sự xuất hiện của phương pháp phần tử hữu hạn và công nghệ máy tính điện tử tạo nên bước ngoặt quyết định cho ngành cơ học kết cấu vỏ. Các phương pháp giải tích cổ điển chỉ giải quyết được các dạng hình học cơ bản với điều kiện biên lý tưởng. Trong thực tế, các công trình vỏ có hình dạng bất kỳ, chiều dày thay đổi và chịu tải trọng phức tạp. Phương pháp phần tử hữu hạn cho phép rời rạc hóa kết cấu vỏ thành các phần tử phẳng hoặc cong nhỏ. Máy tính số hóa các ma trận độ cứng và giải hệ phương trình đại số quy mô hàng triệu ẩn số. Việc mô phỏng chính xác trường ứng suất, biến dạng và dự báo vị trí tập trung ứng suất nguy hiểm trở nên dễ dàng. Quá trình thiết kế kết cấu vỏ mỏng hiện nay đạt độ chính xác cao, rút ngắn thời gian tính toán và đảm bảo an toàn tuyệt đối cho công trình.

Trong cơ học kết cấu, vỏ mỏng được phân loại chủ yếu dựa vào dấu của độ cong Gauss trên mặt trung hòa. Vỏ có độ cong Gauss dương bao gồm mặt cầu, mặt ellipsoid và paraboloid elliptic. Loại vỏ này có độ cứng không gian rất cao và khả năng chịu tải trọng phân bố cực tốt. Vỏ có độ cong Gauss âm tiêu biểu là mặt hyperboloid một tầng, paraboloid hyperbolic dạng yên ngựa. Dạng vỏ này thường được áp dụng cho tháp làm mát hoặc mái che kiến trúc phức tạp. Vỏ có độ cong Gauss bằng không gồm các mặt có thể trải phẳng như vỏ trụ và vỏ nón. Mỗi nhóm hình học sở hữu cơ chế truyền lực đặc thù. Việc phân loại chuẩn xác giúp kỹ sư lựa chọn mô hình tính toán giải tích hoặc số phù hợp với từng dạng kết cấu thực tế.

Phương trình mặt trung hòa của vỏ mỏng thường được biểu diễn dưới dạng tham số véctơ bán kính theo hai tọa độ đường cong. Khi thiết lập hệ trục tọa độ vi phân trên mặt vỏ, hệ tọa độ cong trực giao được ưu tiên sử dụng. Trong hệ tọa độ này, hai họ đường cong tọa độ cắt nhau vuông góc tại mọi điểm trên mặt trung hòa. Thông thường, các trục tọa độ được chọn trùng với các đường cong chính của mặt cong. Lựa chọn này giúp loại bỏ thành phần độ cong xoắn, làm các biểu thức vi phân trở nên tinh gọn. Các hệ số Lame đo lường tỷ lệ biến thiên của độ dài cung vi phân theo từng tọa độ tham số. Việc xác định chính xác các hệ số Lame và đạo hàm của véctơ pháp tuyến là bước chuẩn bị quan trọng để xây dựng quan hệ biến dạng - chuyển vị.

Hình học vi phân mặt cong sử dụng hai dạng toàn phương cơ bản để mô tả thuộc tính bề mặt vỏ mỏng. Dạng toàn phương thứ nhất xác định bình phương khoảng cách giữa hai điểm lân cận vô cùng bé trên mặt trung hòa. Các hệ số của dạng toàn phương thứ nhất đo lường độ dài, góc giữa các đường cong và diện tích vi phân. Dạng toàn phương thứ hai xác định khoảng cách từ một điểm trên mặt cong đến mặt phẳng tiếp tuyến lân cận. Nó phản ánh độ cong pháp tuyến và độ uốn của mặt vỏ theo các phương khác nhau. Tỷ số giữa dạng toàn phương thứ hai và dạng toàn phương thứ nhất cho phép tính trực tiếp độ cong pháp tuyến tại điểm khảo sát. Hai dạng toàn phương này cung cấp cơ sở toán học để tính biến dạng kéo nén và biến dạng uốn trong lý thuyết vỏ mỏng.

Để đơn giản hóa bài toán cơ học không gian ba chiều về bài toán hai chiều trên mặt trung hòa, lý thuyết vỏ mỏng áp dụng hệ giả thiết Kirchhoff - Love. Giả thiết thứ nhất khẳng định đoạn thẳng thẳng góc với mặt trung hòa trước biến dạng vẫn giữ nguyên độ thẳng và vuông góc với mặt trung hòa sau biến dạng. Giả thiết này triệt tiêu biến dạng trượt ngang theo chiều dày. Giả thiết thứ hai coi ứng suất pháp tuyến theo phương vuông góc với mặt trung hòa là rất nhỏ và có thể bỏ qua so với ứng suất trong mặt phẳng. Giả thiết thứ ba giả định chuyển vị của các điểm trong vỏ là bé so với chiều dày vỏ, cho phép áp dụng phương trình biến dạng tuyến tính. Hệ giả thiết này giữ vai trò kinh điển, giúp thiết lập các công thức giải tích trực quan và chính xác cho vỏ mỏng.

Xét một phân tố vỏ vi phân vô cùng bé được tách ra từ mặt trung hòa dưới tác dụng của tải trọng ngoài. Các thành phần ứng lực và mômen tác dụng lên các mặt cắt phân tố phải thỏa mãn các định luật tĩnh học. Chiếu toàn bộ hệ lực lên ba trục tọa độ đường cong thu được ba phương trình vi phân cân bằng lực. Lấy tổng mômen đối với ba trục tọa độ cho ba phương trình cân bằng mômen. Trong đó, hai phương trình mômen uốn cho phép biểu diễn lực cắt ngang theo đạo hàm của các mômen. Phương trình mômen quay quanh trục pháp tuyến biểu diễn sự đối xứng của các lực cắt tiếp tuyến. Kết hợp các phương trình tĩnh học tạo ra hệ phương trình vi phân đạo hàm riêng mô tả cân bằng của phân tố vỏ mỏng đàn hồi.

Có hai hướng giải quyết chính cho bài toán biên của lý thuyết vỏ mỏng đàn hồi. Hướng thứ nhất là phương pháp chuyển vị, chọn ba thành phần chuyển vị của mặt trung hòa làm ẩn số cơ bản. Thay quan hệ biến dạng - chuyển vị và định luật vật lý vào phương trình cân bằng để thu được hệ phương trình vi phân đối với chuyển vị. Hướng thứ hai là phương pháp nội lực và hàm ứng lực, đưa các hàm thế để thỏa mãn tự động phương trình cân bằng. Ngoài ra, nguyên lý tương tự tĩnh học - hình học của Goldenveizer và Lurie chỉ ra sự tương đồng toán học giữa các ten-xơ nội lực và ten-xơ biến dạng. Nguyên lý tương tự này giúp giảm một nửa số lượng biến số độc lập, tạo điều kiện thuận lợi cho việc giải tích toán học hệ phương trình vỏ mỏng.

Để trạng thái phi mômen tồn tại thực sự trong kết cấu vỏ mỏng, cần thỏa mãn đồng thời bốn điều kiện cơ bản. Thứ nhất, mặt trung hòa của vỏ phải là mặt cong liên tục và trơn, không có điểm gãy khúc hoặc thay đổi độ cong đột ngột. Thứ hai, chiều dày của vỏ phải không đổi hoặc biến thiên rất từ từ theo các phương tọa độ. Thứ ba, tải trọng tác dụng lên mặt vỏ phải là tải phân bố liên tục, không có lực tập trung hoặc mômen uốn cục bộ. Thứ tư, điều kiện liên kết tại các mép biên phải cho phép vỏ tự do biến dạng theo phương pháp tuyến hoặc lực đỡ biên phải trùng với tiếp tuyến của mặt vỏ. Nếu một trong các điều kiện trên bị vi phạm, mômen uốn sẽ xuất hiện và phá vỡ trạng thái phi mômen thuần túy.

Vỏ tròn xoay là dạng kết cấu phổ biến nhất trong thực tế nhờ tính đối xứng trục hình học. Khi áp dụng lý thuyết phi mômen cho vỏ tròn xoay, hệ phương trình cân bằng vi phân đơn giản hóa đáng kể. Phương trình cân bằng hình chiếu theo phương pháp tuyến dẫn đến phương trình Laplace liên hệ giữa lực pháp tuyến kinh tuyến và lực pháp tuyến vĩ tuyến với áp lực mặt ngoài. Phương trình cân bằng theo phương tiếp tuyến kinh tuyến cho phép tích phân trực tiếp để tìm lực dọc kinh tuyến. Bằng cách chiếu toàn bộ tải trọng phần vỏ phía trên mặt cắt vĩ tuyến, lực kinh tuyến được xác định hoàn toàn thông qua biểu thức giải tích tường minh. Sau khi có lực kinh tuyến, lực vĩ tuyến được suy ra trực tiếp từ phương trình Laplace mà không cần tích phân phức tạp.

Sau khi xác định được các thành phần lực màng, bài toán tiếp theo là tính toán trường chuyển vị đàn hồi của vỏ tròn xoay. Dưới tác dụng của tải trọng đối xứng trục, chuyển vị xoay vòng quanh trục triệt tiêu, chỉ còn lại chuyển vị tiếp tuyến theo đường sinh và chuyển vị pháp tuyến. Sử dụng định luật Hooke cho trạng thái ứng suất phẳng để tính các thành phần biến dạng dài kinh tuyến và vĩ tuyến. Từ quan hệ hình học giữa biến dạng và chuyển vị, chuyển vị pháp tuyến được tính trực tiếp từ biến dạng vĩ tuyến và bán kính cong. Chuyển vị tiếp tuyến thu được bằng cách giải một phương trình vi phân thường bậc nhất với điều kiện biên cố định tại gối tựa. Việc xác định trường chuyển vị giúp kiểm tra độ cứng và đánh giá khả năng tương thích biến dạng tại các liên kết biên.

Hiệu ứng biên là hiện tượng cơ học đặc trưng và nổi bật nhất trong lý thuyết mômen của vỏ mỏng. Khi mép biên của vỏ bị ngàm cứng hoặc chịu lực cắt cục bộ, mômen uốn phát sinh rất lớn tại biên. Tuy nhiên, do độ cứng uốn của vỏ nhỏ, sóng ứng suất uốn này suy giảm rất nhanh theo hàm mũ khi tiến sâu vào trong lòng vỏ. Hiện tượng dao động tắt dần này được gọi là hiệu ứng biên cục bộ. Vùng chịu tác động của hiệu ứng biên chỉ chiếm một dải hẹp có chiều rộng tỷ lệ với căn bậc hai tích giữa bán kính cong và chiều dày vỏ. Bên ngoài vùng này, kết cấu hoàn toàn làm việc theo trạng thái phi mômen. Hiểu rõ quy luật hiệu ứng biên giúp kỹ sư tăng cường chiều dày cục bộ tại gối tựa mà không làm tăng trọng lượng toàn công trình.

Đối với vỏ tròn xoay chịu tải trọng đối xứng trục, lý thuyết mômen dẫn đến một hệ phương trình vi phân bậc bốn phức tạp. Nhà cơ học Meissner đã có đóng góp đột phá khi biến đổi hệ phương trình này thành hai phương trình vi phân bậc hai liên kết. Hai hàm ẩn cơ bản được chọn là góc xoay tiếp tuyến của pháp tuyến và đại lượng lực cắt biến đổi. Hệ phương trình Meissner giải quyết tường minh sự phân bố mômen uốn, lực cắt và lực màng dọc theo đường sinh vỏ. Nghiệm của phương trình Meissner thường biểu diễn dưới dạng các hàm Bessel, hàm Kelvin hoặc hàm mũ lượng giác phức tắt dần. Việc áp dụng các hàm giải tích này cho phép tính toán chính xác ứng suất uốn cực đại tại mép ngàm của vòm cầu, vỏ nón và vỏ trụ tròn xoay.

Trong các công trình xây dựng thực tế, hình học vỏ mỏng thường không đồng nhất và chịu tải trọng bất kỳ. Các phương pháp giải tích cổ điển bộc lộ nhiều hạn chế khi giải các bài toán có hình học phức tạp. Việc triển khai các thuật toán số trên máy tính giải quyết hoàn toàn bài toán kết cấu vỏ mỏng tổng quát. Phương pháp phần tử hữu hạn sử dụng các phần tử vỏ phẳng nhiều bậc tự do hoặc phần tử vỏ cong đẳng tham số. Ma trận độ cứng của toàn hệ thống được lắp ghép tự động và giải thông qua các thuật toán đại số ma trận lớn. Các phần mềm tính toán hiện đại xuất ra trường ứng suất tương đương von Mises, biểu đồ mômen và biến dạng 3D trực quan. Đây là công cụ đắc lực hỗ trợ kỹ sư tối ưu hóa hình dáng và đảm bảo độ an toàn chịu lực của công trình.