5. Nội dung nghiên cứu
2.1.1. PLAXIS 3D Foundation
Plaxis được nghiên cứu tại Đại học Delft vào năm 1987, là phần mềm sử dụng phương pháp phần tử hữu hạn để phân tích khu vực ven sông ở Hà Lan. Mãi đến năm 2001 mới có mô hình 2D và khi đó PLAXIS 3D Tunnel đã được phát hành. Ba năm sau, PLAXIS 3D Foundation được tạo ra, nó được phát tiển nhằm phục vụ tính toán nền móng công trình và dùng trong ngành địa kỹ thuật như tính toán ổn định, độ lún và biến dạng.
18
Khi nói rằng chuyển từ phân tích 2D sang 3D rất khả thi nhưng đồng thời sự phức tạp cũng tăng lên. Điều này cũng chính là tăng khối lượng tính toán. Nếu một tính toán chi tiết chỉ mất vài phút để thực hiện trong phân tích 2D, nhưng tính toán tương tự trong phân tích 3D có thể mất hàng giờ. Mô hình phân tích ảnh hưởng của các giai đoạn thi công hố đào đến cọc bên trong hố đào cũng là một trong các trường hợp đó. PLAXIS 3D Foundation gồm ba phần chính, đó là mô hình (model),
tính toán (calculation) và xuất kết quả (Output).
2.1.2. Mô hình
Trong chế độ mô hình thì hình dạng của mô hình được xây dựng. Biên của các lớp đất và đặc trưng vật liệu được thiết lập. Xây dựng các phần tử như tường và dầm tại các vị trí trong mô hình và đặc trưng tiếp xúc được định nghĩa. Cuối cùng lưới được tạo ra và đạt một độ mịn thích hợp. Trong đó việc lựa chọn mô hình đất là rất quan trọng, và sẽ được trình bày trong phần sau.
2.1.3. Tính toán
Trong chế độ tính toán, một số bước tính toán có thể được tạo ra. Khác nhau trường hợp tải và hình dạng được thiết lập để mô phỏng trình tự xây dựng công trình thực tế. Đối với mỗi bước có thể thiết lập các điều kiện mực nước ngầm khác nhau, các bộ phận công trình có thể được kích hoạt. Hố đào được mô phỏng bằng cách chấm dứt hoạt động của các tập hợp. Các loại tính toán phải được định nghĩa có thể là tính dẻo (Plastic) hoặc cố kết (consolidation). Phân tích cố kết
(consolidation) được sử dụng khi mô hình các ứng xử phụ thuộc vào thời gian như sự phát triển và tiêu tán áp lực nước lỗ rỗng hay khi yêu cầu tính toán độ lún do từ biến. Tính toán tích dẻo (Plastic) được dùng để phân tích biến dạng đàn – dẻo
(elastic-plastic) theo lý thuyết biến dạng nhỏ (Brinkgreve, 2007). Ứng suất và biến dạng được tính toán cho tất cả các nút còn trong trạng thái giới hạn.
2.1.4. Xuất kết quả
Trong phần chính thứ ba của Plaxis là chế độ xuất kết quả tính toán và được dùng xử lý kết quả tính toán. Biến dạng, ứng xuất và áp lực nước lỗ rỗng sẽ được
thể hiện trong mỗi bước tính toán, còn đối với các cấu kiện công trình ta có thể xem được moment uốn và lực cắt.
2.2. Tạo mô hình
Để phân tích phần tử hữu hạn trên phần mềm PLAXIS 3D Foundation thì điều quan trọng đầu tiên là phải tạo mô hình hình học cho bài toán. Mô hình này mô tả cấu trúc của công trình trong không gian 3 chiều được chương trình định nghĩa thông qua các mặt phẳng làm việc và các hình trụ hố khoan địa chất. Mô hình bao gồm các lớp địa tầng, kết cấu của công trình và các loại tải trọng. Mô hình phải đủ lớn để biên bài toán không ảnh hưởng đến kết quả phân tích.
Mặt phẳng làm việc (Work Planes): là các mặt phẳng nằm ngang theo trục x - z tương ứng với một cao độ y.
Điểm và đường thẳng (Geometry line): dùng để tạo mô hình hình học cho bài toán.
Phần tử dầm (Beam): dùng để mô hình cho kết cấu thanh mảnh chịu uốn và lực dọc trục như dầm móng…
Phần tử sàn (Floor): dùng mô phỏng cho kết cấu có chiều dày nhỏ theo phương ngang và chịu uốn như bản móng…
Phần tử tường (Wall): dùng mô hình cho kết cấu có chiều dày nhỏ theo phương đứng và chịu uốn như vách tầng hầm…
Phần tử cọc (Pile): dùng mô hình cho các loại cọc.
Phần tử lò xo (Spring): dùng để gắn kết vào một mặt của kết cấu và khống chế mặt đó so với mặt khác. Phần tử này thường dùng mô phỏng sự làm việc của cọc đơn.
Phần tử biên (Line Fixity): dùng để tạo biên khống chế cho bài toán. Phần tử hố khoan (Borehole): dùng khai báo các lớp địa chất cho bài
20
Hình 2.1 – Yêu cầu tối thiểu của mô hình hố đào (Bakker, 2005)[5]
2.3. Chia lƣới phần tử
Để thực hiện tính toán phần tử hữu hạn,mô hình hình học trong PLAXIS 3D Foundation phải được chia thành các phần tử nhỏ hơn, được gọi là chia lưới phần tử hữu hạn. Mỗi phần tử bao gồm một số lượng nhất định các nút hình thành hệ thống phương trình cho việc tính toán. Số lượng các nút nhiều sẽ kéo theo hệ thống phương trình lớn hơn để máy tính giải quyết . Trong 2D mỗi nút có hai bậc tự do, tức là các nút có thể di chuyển theo phương x và y, xem Hình 4.2
Hình 2.2 – Các phần tử và nút trong một mô hình 2D. Mỗi nút có hai bậc tự do,
được mô tả bởi các mũi tên trong hình nhỏ hơn, (Wiberg, 1974)[4]
Khi mô hình trong không gian ba chiều, mỗi nút có ba bậc tự do, kết quả là sẽ cho một hệ thống phương trình lớn hơn vì thực tế rằng mỗi nút cũng có thể di chuyển theo phương z.
Phần tử Nút
Phân tích phần tử hữu hạn sẽ tiến hành theo ba bước, xem hình 2.3. Bước đầu tiên là phân chia mô hình thành các phần nhỏ hơn bằng cách tạo ra các phần tử (Generate elements) vào mô hình, trong đó mỗi phần tử tương đối dễ để giải từng cái một. Bước tiếp theo là phân tích các phần tử (Element analysis). Bước cuối cùng là phân tích hệ thống (System analysis) nơi mà tất cả các phần tử được kết nối với một hệ thống bằng các điều kiện biên, (Wiberg, 1974).
Hình 2.3 – Các bước phân tích phần tử hữu hạn (Wiberg, 1974)[4]
Phân tích phần tử hữu hạn là một phương pháp tính gần đúng và nguồn lỗi nhiều và thường không thể tránh khỏi (Wiberg, 1974). Số lượng nút trong mô hình có tác động đáng kể đến kết quả tính toán. Số lượng phần tử lớn hơn tạo ra một số lượng lớn các nút cho kết quả chính xác hơn mô hình có ít nút. Khi thiết kế các mô hình phức tạp thì thích hợp có số nút cao hơn so với trường hợp đơn giản hơn.
Các thiết lập mặc định cho kích thước cluster trong PLAXIS 3D Foundation là lưới thô (Coarse mesh). Chia lưới thế này có thể đủ kích thước khi mô hình các trường hợp đơn giản và tính chính xác của tính toán không cần cao. Hệ thống phương trình dễ dàng hơn cho máy tính để giải quyết và thời gian tính toán tương đối ngắn.
Nếu cần độ chính xác cao hơn thì cần sự làm mịn các cluster. Nếu toàn bộ mô hình cần được làm mịn, thì sử dụng chức năng Global coarseness để thay đổi kích thước phần tử cả phương đứng và ngang với khoảng từ rất thô (Very coarse) đến rất mịn (Very fine). Lưu ý rằng nếu chúng ta chia lưới 2D quá mịn sẽ làm tăng phần tử khi tiến hành chia lưới 3D và đồng nghĩa là thời gian tính toán sẽ tăng lên.
22
Để xác định mô hình gồm bao nhiêu nút là đủ, Hannes và Daniel (2010) đã thực hiện mô phỏng với số lượng các nút là tăng dần dần trên cùng một mô hình. Kết quả được so sánh với chuyển vị, xem hình 2.4.
Hình 2.4 – Kết quả chuyển vị với số nút tăng dần trong mô hình 3D, (Hannes và
Daniel, 2010) [4]
2.4. Mô hình ứng xử của đất