Ứng dụng mạng lưới lò xo vào xử lý ổn định toàn khối nền đất yếu tại đồng bằng sông cửu long

Tạ chí Vật liệ Ứng ụng mạng lưới lị xo xử l ổn định toàn gu ễn hế nh đồng ông ửu Long àX ựng số 21 3 hối đất ếu Khoa Xây Dựng, Trường Đại Học Sư Phạm Kỹ Thuật Hồ Chí Minh TỪ KHO Xử lý đất yếu Gia cố nông T M TẮT Nội dung bài báo này nhằm nghiên cứu hiệu phương pháp gia cố nông đất yếu khu vực đồng bằng sông Cửu Long thông qua phương pháp mô số Phương pháp trộn nông xi măng đất Xi măng trộn đất đề xuất và hiệu đánh giá thông qua mô hình LSMnS; là mô hình Mạng lưới lò xo trường rỗng với hay nhiều thành phần rắn khác Tính đúng đắn LSMnS kiểm tra dựa Môi trường rỗng LSM LSMnS phát triển từ mô hình mạng lưới lò xo đàn hồi (LSM) để tính toán các tính chất học môi so sánh với các cơng cụ và mơ hình sẵn có, sai số lớn nhỏ 1,6 % Mô hình này áp dụng để tính các mô đun đàn hồi hữu hiệu số mẫu đất với tỷ lệ xi măng trộn khác nhằm đánh giá tăng cường tính chất các mẫu trộn này Kết cho thấy mô đun cứng và cắt hữu hiệu tăng tới 47 % và 50 % hàm lượng xi măng mức 13 % KEYWORDS ABSTRACT Soil improvement delta by modelisation The shallow soil cement mixing method was proposed and its performance can be Cement soil mixing method Porous media Lattice Spring Models LSMnS The purpose of this paper is to study the effectiveness of various soft soil improvement methods in the Mekong evaluated by lattice spring model for many solids (LSMnS) that was extended from the basic model (LSM) to determinate effective mechanical properties of porous media with one or many solid components This model was validated by comparisons with the other existing models, the maximal difference was always less than 1.6 % Various soil cement mixed samples with different cement proportions were simulated to investigate their different physical properties in order to evaluate the shallow soil cement mixing method It is demonstrated that the effective bulk modulus Ke and shear modulus Ge are increased by 47 % and 50 %, respectively, when the cement proportion is equal to 13 % iới thiệu Vấn đề xử lý đất yếu xuất từ lâu và nghiên mét Độ rỗng đất giảm xuống đất nén chặt thêm nên từ tính chất cọc cát cải thiện Phương pháp đầm chặt đất bằng đầm rung sâu đề xuất cứu nhiều giới Hiện có nhiều phương pháp các nhà Moseley và Kirsch (1993), thép đưa xuống đến độ sâu thông công bố các tạp chí Điều nói lên tầm quan trọng thẳng đứng làm cho đất đầm chặt thép triển sở hạ tầng không Việt Nam mà còn giới toàn phần lớn cho phép, độ lún lệch lớn cho phép; sức chịu tải cát, giếng cát phổ biến ngày Ống thép cắm vào đất làm cho lỏng khoa học đề xuất để xử lý đất yếu ngành xây dựng và giao nghiên cứu này nhằm giải các vấn đề dân sinh phát Tsytovich và cộng (1974) đề xuất phương pháp sử dụng cọc tảng ép chặt lại Đất nén chặt thêm quá trình tạo lỗ khoan Còn nước đất thì bị nén cho thoát ngoài lỗ khoan quy định và máy rung tạo dao động theo phương ngang và phương rút dần lên Độ sâu và khoảng cách đầm phụ thuộc vào các yếu tố: độ lún yêu cầu; sức kháng cắt cần đạt và mức tăng sức kháng chống hoá Phương pháp trộn vôi và xi măng sâu với cách trộn khô trình bày nghiên cứu Chida (1982) Cọc vôi thường nhồi cát Thi công kiểu này đơn giản và bằng vật liệu tương dùng để xử lý, nén chặt các lớp đất yếu nhưbùn, sét và sét pha đệm cát, các loại móng bê tơng, các loại chun dùng để xử lý đất chặt, đường kính cọc vôi tăng lên 20 % làm cho đất xung quanh nén đối rẻ tiền từ sạn sỏi cát thô có chi phí tương đối thấp khác Sản phẩm này chuyên dùng để gia cố cho phần đất yếu dày trạng thái dẻo nhão Việc sử dụng cọc vơi có những tác dụng sau đầm chặt lại ngoài toả nhiệt lượng lớn làm cho nước lỗ rỗng Tác giả liên hệ: Ntanh hcmute.edu.vn Nhận ngày 21/04/2021, giải trình ngày 10/05/2021, chấp nhận đăng 12/07/2021 JOMC 73 Tạ chí Vật liệ bốc làm giảm độ ẩm và tăng nhanh quá trình nén chặt Sau xử lý bằng cọc vôi đất cải thiện đáng kể: độ ẩm đất giảm đến %, lực dính tăng lên khoảng 1,5 đến lần Việc chế tạo cọc đất – ximăng giống cọc đất – vôi, hàm lượng ximăng từ % đến 15 % và kết cho thấy sức kháng xuyên đất tăng lên từ đến lần so với chưa gia cố Trong phương pháp điện thấm (Eggestad, 1983), điện cực đưa vào đất và nối với nguồn điện chiều thì nước lỗ rỗng từ cực dương tới cực âm Do hiệu ứng vật lý này, khu vực gần cực dương có áp lực nước lỗ rỗng thấp mở rộng dần theo thời gian và độ bền đất tăng hệ cố kết đất còn độ nén lún đất thì giảm đi) Trong đất nén lún, điện thấm đẩy nước đến cực âm và nước hút nên không quay trở lại cực dương Sự cố kết đất xảy tương ứng với thể tích đất hút Các phương pháp này tỏ hiệu và áp dụng mạnh mẽ thực tế Tuy nhiên, vẫn còn nhiều khó khăn tồn các địa àX ựng số 21 3 chủ yếu là sức căng bao D BOX truyền vào vật liệu độn bên trong; áp lực truyền xuống lớp đất xấu phía là không đáng kể Tuy nhiên phương pháp này gặp khó khăn việc thi cơng đại trà túi D BOX là không phổ biến dẫn đến khả áp dụng rộng rãi là hạn chế Trong nghiên cứu Nguyễn Sỹ Hùng và Vương Hoàng Thạch (2020) bàn các phương pháp gia cố nông đất phổ biến đồng bằng sông Cửu Long phương pháp gia cố bằng cọc cừ tràm, cọc tre; phương pháp gia cố bằng cọc đá chẻ, cọc bê tông cốt thép, cọc xi măng đất trộn sâu và trộn nông Theo kết nghiên cứu này, phương pháp trộn đất chỗ với phương pháp trộn nông là phương pháp hữu hiệu để nhằm cải thiện các đặc tính kỹ thuật và môi trường lớp đất mềm bị ô nhiễm Từ các kết thấy rằng, phương pháp trộn nơng xi măng đất để xử lý đất yếu, tạo ổn định toàn khối là giải pháp phù hợp với khu vực đồng bằng sơng Cửu Long Do đó, đề tài này chất khác có hiệu khác với phương pháp tập trung nghiên cứu hiệu phương pháp trộn nông xi măng cứu và ứng dụng Việt Nam, nhiên, khác địa chất chính là sét nhão bùn sét yếu Tuy nhiên, điểm đề tài là Hầu hết các phương pháp nêu nghiên nên cần có điều chỉnh để đem lại hiệu tốt điển hình phương pháp cọc cát (Nguyễn Đình Đức, 2009, Nguyễn Uyên, 2011), cọc xi măng đất (TCVN 385 2006, TCCS 05: 2010/VKHTLVN), bấc đất áp dụng cho địa chất điển hình khu vực, có thành phần khơng sâu vào nghiên cứu thực nghiệm, mà đánh giá hiệu phương pháp bằng cách xác định tính chất đất sau cải tạo bằng phương pháp mô số Các mẫu đất trộn xi măng với các tỷ thấm đứng (PVD) (TCVN 9355 2012), công nghệ khoan phụt cao áp lệ khác số hóa, sau mơ bằng phương pháp mạng phương pháp nào loại địa chất công trình là nhiệm vụ độ cứng hữu hiệu (Ke), mô đun cắt hữu hiệu (Ge), mô đun Young, Thời gian gần đây, số phương pháp đại khác đất Phương pháp gia cố nông nhắc đến nhiều trước đây, bài (Nguyễn uốc Dũng và cộng sự, 2005)… Trong thực tế, việc lựa chọn khó khăn kỹ thuật và kinh tế lưới lò xo (Lattice spring model) để tính toán các đặc tính mẫu hệ số Poisson nhằm đánh giá hiệu việc trộn xi măng mẫu các nhà khoa học nước nghiên cứu và đề xuất để phù hợp với báo này trình bày nghiên cứu phát triển mô hình LSM, LSMnS và áp cứu Nguyễn Minh Đức và cộng (2018, 2019, 2020) Việc xác định các tính chất học môi trường rỗng trở điều kiện đặc thù khu vực đồng bằng sông Cửu Long Các nghiên rằng, đất điển hình đồng bằng sơng Cửu Long có thành phần chủ yếu là đất sét bùn sét Tác giả đề xuất sử dụng đệm cát kết hợp với sử dụng vải địa kỹ thuật nhằm làm tăng khả thoát nước dụng chúng vào mô các mẫu đất trộn xi măng với tỉ lệ khác thành chủ đề chú ý, có ứng dụng lớn các mảng học, vật lý, vật liệu và tính toán động đất Các vật liệu rỗng phát triển và dùng rộng rãi nhiều ngành công đất tăng cường khả cố kết đất giải pháp nghiệp có ngành xây dựng lại hiệu tốt để cải tạo đất san lấp bằng sét yếu bùn nạo rỗng phức tạp và khó khăn, nhiều phương pháp sử dụng nông thôn đồng bằng sông Cửu Long, nhiên việc áp dụng với công 2007), Buxton và cộng (2005), (Boutin và Auriault, 1990), phương để gia cố đất yếu Các kết rằng phương pháp này đem vét lòng sông phù hợp với các loại công trình đường giao thông trình xây dựng dân sinh, cơng trình móng nơng còn đem lại nhiều khó khăn, chưa thực hiệu Nguyễn Sỹ Hùng và Hoàng Anh (2020) đề xuất phương pháp sử dụng túi đất D box tương tự viện nghiên cứu Nhật Bản để gia cố đất yếu bằng cát san lấp khu vực tỉnh An Giang Kết nghiên cứu rằng Việc sử dụng móng bao D BOX là giải pháp tốt cho nhà thấp tầng (có quy mơ tầng) Sức chịu tải móng Trước đây, việc tính toán các tính chất học môi trường phương pháp đồng (homogenization method), (Malinouskaya, pháp số học (Nemat Nasser và Iwakuma, 1982), (Torquato, 1998) Cohen (2004)… Tuy nhiên các phương pháp này có nhiều hạn chế khó khăn áp dụng cho các mơi trường có thành phần phức tạp kích thước quá lớn Sau đó, với tiến khoa học công nghệ, mạnh mẽ vượt bậc các máy tính, nhiều công cụ số đề xuất và sử dụng để giải vấn đề ví dụ công cụ giải số FMD Malinouskaya (Malinouskaya, 2007) … JOMC 74 Tạ chí Vật liệ Gần đây, mô hình mạng lưới lò xo (LSM) phát triển mở rộng cho môi trường rỗng Tuy nhiên, mô hình này sử Bài báo này giới thiệu các tính chất LSM mô hình, các loại lò xo, phàn tử đàn hồi, các điều kiện biên… và mở rộng thành mô hình LSMnS để áp dụng cho môi trường rỗng với nhiều pha rắn khác Sau đó, mơ hình LSMnS áp dụng để tính toán tenso độ cứng hữu hiệu và các mô đun đàn hồi (mô đun cứng, cắt…) số môi trường rỗng giả định và thực tế Các kết này so sánh với các phương pháp số học trước (Malinouskaya, 2007), (Nemat Nasser và Iwakuma, 1982), (Torquato, 1998) Cohen (2004) và 3 (5) Nếu lò xo góc hợp lò xo đơn b và b’, biến dạng góc tính bằng dụng cho mơi trường rỗng với pha rắn các môi trường thực tế thường có nhiều pha rắn và rỗng kết hợp với ựng số 21 b u −u u= m n Padzniakou, (2012) sử dụng các máy tính mạnh mẽ để mô và tính toán các môi trường đàn hồi kích thước lớn., sau àX    (b,b ) =  (b ) −  (b) (6) Trong khối lập phương bản, độ dài tất các lò xo đơn bằng nửa độ dài ban đầu Do đó, hệ số đàn hồi lò xo  (b) = 2 (b) với  (b ) là hệ số tương ứng với độ dài ban đầu Lúc đó, tổng lượng khối lập phương đơn vị xác định bằng công thức (7)   =  (b) | u(b) |2 +   (b,b ) ( (b,b ) )2 E 2b b,b với  (b) là hệ số lò xo góc với kết thực nghiệm [7] Các phép so sánh này cho thấy khác biệt nhỏ ( %) Do đó, mơ hình LSMnS cho thấy khả và chính xác Mơ h nh mạng lưới lò xo B ic L ttic S ring Mo 2.1 3D LSM, lò xo đơn, lị xo góc l (a) (b) (c) (d) Phương pháp mạng lưới lò xo dùng để mô mơi trường đàn hồi, sau mở rộng để tính toán các môi trường rỗng (Padzniakou, 2012) Trong mô hình 3D, thành phần rắn thể lưới lập phương với 18 lò xo (Hình 1a) Các vecto tọa độ ứng với nút lân cận loại (màu xanh, Hình 1a) sau c (b) (1) = ( a,0,0),(0,  a,0),(0,0,  a), b = 1, ,6 12 nút lân cận loại xác định bằng các vecto sau (màu đỏ, Hình 1a) c(b) =(a, a,0),(0, a, a),(a,0, a), b =7, ,18 (2) nh Mô hình 3D LSM với loại lò xo a) 18 lò xo đơn b) Lò xo góc π/4 c) Lò xo góc π/3 với a là kích thước cạnh khối lập phương (màu đỏ, Hình 1b) Vecto đơn vị trương ứng theo phương xác định (b) n(b) = c | c(b) | (3) Một vài nút mạng xem nằm trung tâm khối lập phương (Hình 1b), đó, các vecto tương ứng là (b) r(b) = c (4) Trong LSM, có loại lò xo tính toán là lò xo đơn nối hai nút mạng gần và và lò xo góc tạo hai lò xo đơn mặt phẳng Trong đơn vị lưới lập phương, có 24 lò xo góc loại π/3 và 24 lò xo góc loại π/4 thể Hình 1c,d Năng lượng đơn vị là hàm phức tạp biểu Năng lượng tất lò xo đơn khối đơn vị xác định Ec =   (b) | u(b) |2 2b (8) Theo Pazniakou (2012), độ dịch chuyển tương đối u(b) tính thơng qua tenso ứng suất sau (b) | r(b) |  n(b) ui(b) = = ij r j ij j Với chuyển vị bé, chúng ta có (10) | u(b) |= u(b)  n(b) =| rb | ij n(jb)ni(b) Do đó, lượng c từ công thức (8) thể viết lại sau Ec 1= =  (b) | u(b) |2  |rb | ni(b) n(jb) | n(b) n(b)ij diễn qua chuyển vị và biến dạng góc Nếu các điểm nút lò xo k l kl 2 2b b đơn là m và n thì chuyển vị tính theo (9) (11) Và mật độ lượng đàn hồi khối đơn vị JOMC 75 Tạ chí Vật liệ (12) Ec 1 b (b) (b) (b) (b) W= c V= [ |r | ni n j | nk nl ]ij kl a b Từ đó, ta xác định tenso độ cứng đơn tương ứng với các lực tập trung (13) C (cf ) =  |rb | ni(b) n(jb) n(b) n(b) ijkl a3 b k l Với giá trị tính sau  (b) , các hệ số tenso độ cứng 3= =  cf cf cf Cxxxx ; Cxxyy ; Cxyxy = a a a (14) Tương tự, từ mật độ lượng lò xo góc khối đơn vị, các hệ số tenso độ cứng góc tương ứng xác định sau 8 ; C (af ) = 4 (af ) = (af ) 6 Cxxxx xxyy − ; Cxxxx = a3 a a3 vị cần chú ý thêm hệ số nhớt àX  ựng số 21 3 nhằm giảm rung động mô (Buxton và cộng sự, 2005) Thuật toán Velvet thể sau a (t ) vi (t +  t ) = vi (t ) + i  t 2 ui (t +  t ) = ui (t ) + vi (t +  t ) t Fi (t +  t )  t) − vi (t +  t ) (t += mi (20) a (t +  t ) t vi (t +  t ) = vi (t +  t ) + i 2 Theo Padzniakou (2012), (Nguyen, 2015), giá trị tối ưu  có (15) thể dung là 0,5 và mức độ vĩ mơ (macroscopic level), tính chất đàn Từ đó, ta xác định tenso độ cứng tổng từ các tenso C( eff ) Khi quá trình đạt cân bằng, chúng ta nhận tenso ứng ( eff ) suất trung bình và C tính dựa công thức (eff ) (21) =  σ C : ε σ và ε là các giá trị trung bình ứng suất và biến đơn và góc tương ứng là 3 + 8 ; C  − 4 ; C  + 6 (16) Cxxxx = xxyy = xxxx = a a3 a a3 a a3 Do đó, các tính chất mơi trường xác định sau  − 4  − 4 ;  =  6 5 ;  =a2 s = s a+ 3; K = a a3 a a 4 + 4 a2 Với Poisson  s và  s (17) là các hệ số Lamé, K mô đun cứng và 2.2 Xác định ứng s ất hữ hiệ th ật toán lặ Vel et  hệ số xton cộng hồi môi trường rỗng thể tenso độ cứng hữu hiệu dạng cấp độ vĩ mô 2.3 Một số ấn đề ới mơ hình LSM bản, điề iện bi n Trên thực tế, các thí nghiệm mô môi trường đàn hồi và môi trường rỗng bằng mô hình LSM với các tính chất nêu mang đến hai vấn đề Theo Pazdniakou, lực tác động các lò xo đơn lên các nút mạng i thể = Fi   ((un − ui ).cˆi,n )cˆi,n n (18) (a) cˆ i , n là vecto chuẩn hóa nối các nút i và n Và lực tác động lò xo góc b i b’ lên các nút cuối b, b’ và nút góc i tính bằng  F(b) =  (b,b ) (b,b) n(b)  n (b,b) (b) | n(b)  n || c |   F(b ) =  (b,b ) (b,b)  n(b )  n (b,b) (b)  | n(b )  n || c | (19)  = F(i) (F(b) + F(b ) ) Để xác định tenso độ cứng hữu hiệu môi trường rỗng bằng mô hình LSM, điều kiện biên không gian theo chu kỳ (spatial periodic boundary condition) cần áp dụng (xem chi tiết (Malinouskaya, 2007), (Padzniakou và Adler, 2013)) Sau đó, thuật toán Vevelt áp dụng để xác định phát triển trường chuyển (b) nh Mô kéo đơn giản theo phương x vật liệu đàn hồi đẳng hướng a) Kết theo lý thuyết đàn hồi b) Kết mô trực tiếp bằng LSM Vấn đề đầu tiên là hệ số Poisson, với cơng thức (17), dao động từ tới 0,25 hệ số β/α biến đổi từ tới ∞ Tuy nhiên với số vật liệu thực tế, giá trị lớn 0,25 JOMC 76 Tạ chí Vật liệ Vấn đề thứ hai là điều kiện biên, ta thấy với công thức (17), giá trị các hệ số đàn hồi α và β là giống cho lò xo đơn và góc Điều này đúng chúng nằm lòng môi trường; àX ựng số 21 3 Ứng xử đàn hồi môi trường rỗng thể qua công thức (21) Để xác định tenso độ cứng hữu hiệu C( eff ) mẫu rỗng, chúng thực tế thì nhiều lò xo nằm bề mặt bề mặt tiếp xúc ta dùng mô lập trình hệ ngôn ngữ Fortran với mô kéo đơn giản theo phương x khối lập phương làm bằng vật liệu x, , cứng giống vị trí chúng Theo lý thuyết đàn hồi, khối kéo và thí nghiệm cắt đươc tiến hành để xác định lỗ rỗng (Ladd và cộng sự, 1997) Để chứng minh điều đó, mơ đàn hồi đẳng hướng (Hình 2) tiến hành xem các lò xo có độ lập phương biến dạng thành hình hộp chữ nhật Hình 2a Tuy hình LSM và PTĐH bao gồm mô kéo đơn giản theo phương và mô phòng cắt theo các mặt phẳng x , x và Nếu môi trường xem đẳng hướng, cần thí nghiệm C( eff ) nhiên, kết mô trực tiếp với LSM theo điều kiện cho mặt phẳng biến dạng Hình 2b Rõ ràng, hệ số đàn hồi các lò xo nằm bề mặt cần phải khác biệt Để giải vấn đề này, (Ladd và cộng sự, 1997), (Ladd và cộng sự, 1997b) đề xuất ý tưởng phần tử đàn hồi (elastic element) Một phần tử đàn hồi (PTĐH) là khối lập phương thiết lập nút mạng liên kết với bằng các lò xo (Hình 3) và là phần tử nhỏ môi trường đàn hồi Phụ thuộc vào số PTĐH liền kề, nh Mô kéo xác định hệ số lò xo thuộc vài phần tử khác Ví dụ, lò xo đơn hệ số dạng [1 0] thuộc 1, 2, hay phần tử khác nhau; lúc này độ cứng chúng lần lượt là α/4, α/2, 3α/4 và α Còn lò xo đơn dạng [1 0] thuộc hoăc phần tử và độ cứng tương ứng là α/2 α Tương tự với lò xo góc π/4, hệ số là β/2 β ứng với trường hợp thuộc phần tử; riêng loại π/3 thuộc phần tử nên hệ số ln là β Rõ ràng lúc này, các lò xo các vị trí khác Padzniakou, (2012) rằng, PTĐH giải xác định C xyxy Với trường hợp tổng quát, mơ cần tiến có dạng C vấn đề hệ số Poisson nêu bằng cách gán các giá trị âm cho hệ Cxxxx và C xxyy và cắt hành để xác định giá trị tenso độ cứng hữu hiệu mơi trường có các hệ số đàn hồi khác nh Phần tử đàn hồi (elastic element) LSM b ( eff )  Cxxxx   Cxxyy  Cxxzz =     Cxxyy C yyyy C yyzz 0 Cxxzz C yyzz Czzzz 0 0 0 C yzyz 0 0 0 Cxzxz         Cxyxy  0 0 (22) h t triển LSM thành LSMnS cho môi trường rỗng gồm nhiều thành hần rắn h c nh u Mô hình LSM dùng rộng rãi để mô mơi trường rỗng, từ có nhiều áp dụng học, vật liệu, tính toán động đất  lớn 0,25 Điều này là hợp lí (Padzniakou, 2012), (Malinouskaya, 2007), (Buxton và cộng sự, 2005) bù vào lực tác dụng gây giá trị âm β (Ladd và cộng sự, môi trường rỗng gồm pha rắn và lỗ rỗng Điều này gây nhiều khó số đàn hồi β, lúc PTĐH, giữa hai nút cuối lò xo góc, ln có lò xo đơn và 1997), (Ladd và cộng sự, 1997b) Với PTĐH, điều kiện biên lực tác dụng áp dụng Lực phân bố lên PTĐH quy đổi nút mạng thuộc mặt phẳng tác dụng và tương tự hệ số đàn hồi, độ lớn lực phụ thuộc vào số lượng phần tử mà nút mạng thuộc vào (Padzniakou, 2012) 2.4 Mơ hỏng tính tốn tenso độ cứng hữ hiệ môi trường rỗng Tuy nhiên, hạn chế lớn là phạm vi áp dụng bị giới hạn với khăn ứng dụng vì hầu hết môi trường, vật liệu rỗng thực tế thường có ít hai pha rắn kết hợp lỗ rỗng Do đó, mơ hình mở rộng LSMnS phát triển để mơ mơi trường rỗng với nhiều pha rắn nhằm khắc phục hạn chế Mô hình LSMnS phát triển dựa LSM với PTĐH nên những tính chất nêu mục vẫn giữ lại JOMC 77 Tạ chí Vật liệ 3.1 Dạng hình học môi trường àX ựng số 21 3 là mô đun Young, ν hệ số Poisson thành phần rắn Do đó, ta có các hệ số đàn hồi 1 , 1 ứng với pha rắn và  ,  với pha Dạng hình học môi trường nhiều pha rắn rời rạc rắn Một lò xo mạng lưới lúc này thuộc PTĐH loại xử lí bằng kỹ thuật chụp cắt lớp m CT với bề dày lớp bằng kích thước lúc này khơng phụ thuộc vào số lượng mà còn phụ thuộc vào hóa mơi trường pha rắn (Nguyen, 2015) Mẫu rỗng a PTĐH (thường vài làm m ) Lúc này theo phương, mẫu chia Ncx  Ncy  Ncz khối lập phương có kích thước a và đặc trưng hàm số biểu diễn pha (phase function) i, , Với môi trường nhiều pha rắn, pha rắn gán tương ứng số (1,2 ) Trong phạm vi bài báo này, chúng ta xét mơi trường có pha rắn và pha rỗng, lúc đó, hàm số pha trở thành Z(i,j,k)=  khối thuộc pha rỗng   1 khối thuộc pha rắn  khối thuộc pha rắn   (23) Lúc này, để tính toán tác động giữa các pha rắn với và giữa loại hai loại Hình Do đó, hệ số đàn hồi loại PTĐH mà thuộc Lúc này hệ số các lò xo đơn [1 0], [1 0], lò xo góc π/4 hay π/3 xác định sau     + = [1 0] 4     + =  /4 2     + a = [1 0] 2 (25)    +  =  /3 1 2 với η1 và η2 là số lượng PTĐH loại và mà thuộc 3.3 Điề iện bi n th ật toán Trong giới hạn bài báo, chúng ta xét đến môi trường rỗng tuần pha rắn và pha rỗng, các PTĐH tương tự LSM sử dụng hoàn theo khơng gian và hoàn toàn đại diện ô rắn tương ứng với PTĐH loại (màu xanh dương, Hình 5b) với môi trường, biến dạng vỹ mô tương ứng với mô đơn giản Dựa vào hàm số pha trước đó, các khối lập phương sở thuộc pha nút mạng là đỉnh, và pha rắn tương ứng với PTĐH loại (xanh lá) Sự tác động giữa các pha thể qua các lò xo nằm mặt tiếp xúc, chúng làm việc theo hai loại PTĐH mạng sở (Padzniakou, 2012) Để tính toán độ cứng hữu hiệu mục 2.4 áp dụng cho ô mạng sở để xác định ứng xử ô mạng toàn mơi trường (Padzniakou và Adler, 2013) Lúc đó, biến dạng và chuyển vị các mặt đối diện phụ thuộc lẫn và tuân thủ điều kiện tuần hoàn không gian (Burla và cộng sự, 2009) Các điều kiện phân bổ lực phân bố các nút mạng hoàn toàn tương tự mô hình LSM nêu Trong trường hợp tổng quát, ba mô kéo và ba mô cắt theo ba phương là cần thiết để xác định C( eff ) Và môi trường xem là đẳng hướng ( eff ) thì cần thí nghiệm kéo và cắt đơn giản để xác định C thành phần Với mô hình LSMnS, thuật toán Velvet áp dụng hoàn (a) toàn giống LSM các mô xây dựng ngôn ngữ (b) nh Mô hình mạng lưới với pha rắn a) Dạng hình học tương ứng với hàm số pha b) loại PTĐH ứng với pha rắn (xanh dương), pha rắn (xanh lá), màu cam là mặt tiếp xúc 3.2 ệ số đàn hồi lò xo ới nhiề Tuy nhiên, thời gian các mô chủ yếu dành cho vòng lặp LSMnS không quá chênh lệch với LSM Thực tế, với mô hình LSMnS, chúng ta sử dụng tiêu chuẩn hỗ trợ OpenMP (Open Rất rõ ràng rằng độ cứng hữu hiệu môi trường phụ thuộc vào phân bố các thành phần rắn đó, và độ cứng các lò xo nằm các vị trí khác mạng lưới khác Theo lý thuyết đàn hồi và cơng thức (17), ta xác định số đàn hồi lò xo đơn (α) và góc (β) PTĐH tương ứng với pha rắn sau  để xác định chính xác hệ số đàn hồi lò xo ô mạng sở Velvet và bước xác dịnh này nằm ngoài vòng lặp nên thời gian mơ rắn a3E(4 −1) aE  = 5(1− 2 ) 20(2 −1)(1+ ) lập trình Fortran Sự phức tạp LSMnS là phải dùng nhiều thời gian MultiProcesing) Fortran để giảm bớt thời gian mô Chuẩn này giúp mô chạy song song máy tính có nhiều nhân Theo qua tắc Amdahl (Che và Nguyen, 2014), với LSMnS, dùng OMP và chạy song song máy tính nhân, tốc độ chạy nhanh gấp 5,2 lần và 10,4 máy 48 nhân (Nguyen, 2015) Ứng ụng mô hỏng ố ố môi trường rỗng 4.1 Môi trường ới hai rắn có cùng hệ số oisson (24) JOMC 78 Tạ chí Vật liệ àX ựng số 21 3 (Nemat Nasser và Iwakuma, 1982), (Torquato, 1998), Cohen (2004) đề xuất số phương pháp dự báo độ cứng hữu hiêu các môi trường không đồng (Nemat Nasser và Iwakuma, 1982) sử dụng phương pháp cộng hưởng đàn hồi, (Torquato, 1998), Cohen (2004) dùng các khai triển mở rộng để tính toán các các hệ số học hữu hiệu cho môi trường dạng khối lập phương 3D với các khối cầu vật liệu (Hình 6) (a) nh Môi trường pha rắn dạng khối lập phương với các khối cầu vật liệu Pha rắn suốt, pha rắn màu xám Dựa vào các kết đó, chúng ta thực mơ với mơi trường có mạng sở Hình 6, tạo hai pha rắn có các tính chất học thỏa mãn các điều kiện quy định (Nemat Nasser và Iwakuma, 1982), (Torquato, 1998) là = = 0.3 và G / G = 2 với G là mô đun cắt chất rắn Do môi trường trường hợp (b) này là đẳng hướng, chúng ta cần thực mô kéo và cắt theo phương x để xác định các thành phần tenso độ cứng hữu hiệu (22) Từ đó, mơ đun cứng hữu hiêu, mô đun cắt hữu hiệu đầu tiên và thứ hai xác định sau = Ke Cxxxx + 2Cxxyy * Cxxxx − Cxxyy ; Ge C= = yzyz ; Ge (26) Các kết có từ mơ LSMnS cho các mẫu có kích thước và tỉ lệ thể tích pha rắn  = Vran / Vtong khác so sánh với các kết (Nemat Nasser và Iwakuma, 1982), (Torquato, 1998)Cohen (2004) các Hình 7.a, 7.b và 7.c Ta thấy rằng các kết mô và tính toán gần nhau, sai số lớn là 1,6 % Lưu ý rằng các so sánh thực với nhiều kích thước khác môi trường và kết không đem lại nhiều khác biệt; rõ ràng điều này độ chính xác cao LSMnS (c) nh So sánh giữa kết mô bằng LSMnS và (Nemat Nasser và Iwakuma, 1982), (Torquato, 1998)Cohen (2004) (a) Mô Ke / K (b) Mô đun cắt hữu hiệu đầu tiên (c) Mô đun cắt hữu hiệu thứ hai G* / G Ge / G e 2 đun cứng hữu hiệu 4.2 Mơi trường ới hai rắn bất ì JOMC 79 Tạ chí Vật liệ àX ựng số 21 3 LSMnS chứng minh khả làm việc hiệu với mơi trường có pha rắn có hệ số Poisson, tiếp theo, chúng ta kiểm tra với hai pha rắn bất kì bằng cách so sánh với kết công cụ tính toán FMD đề xuất Malinouskaya (Malinouskaya, 2007) FMD phát triển dựa phương pháp thể tích hữu hạn với bài toán rời rạc không gian bậc Môi trường FMD tạo bằng các lưới tứ diện; dựa vào số lượng tứ diện ô mạng sở, FMD có loại lưới mạng là SCT6, SCT24, TUT, SUT Theo (Malinouskaya, 2007), SCT24 là loại dễ mô và đem lại kết chính xác nên bài báo này, chúng ta so sánh kết LSMnS với FMD dùng lưới tứ diện SCT24 (b) (a) (b) nh a) Môi trường tạo công cụ FMD (Malinouskaya, 2007) b) Tứ diện SCT24 ô mạng sở Môi trường các so sánh này có dạng (c) Hình với các tỷ lệ thể tích ρ thay đổi từ 0,05 đến 0,5 Hai pha rắn lúc nh So sánh kết mô với LSMnS (đen) và SCT24 (đỏ) a) này không bị ràng buộc so sánh trước, trường hợp này, K e b) Ge c) Ge* chúng có các mô đun đàn hồi Bảng Bảng Hệ số đàn hồi hai pha rắn dùng so sánh Mô đun đàn Mô đun Mô đun Mô đun Chất rắn 30 20 49,091 hồi Chất rắn cứng 100 cắt 50 Các kết mô so sánh Hình Sự khác Hệ số Young Poisson 128,57 0,2857 0,2273 giữa LSMnS và SCT24 là nhỏ, với tỉ lệ thể tích ρ = 0,1, sai số cho mô đun cứng, mô đun cắt thứ và thứ hai 0,1 %; sai số lớn là 0,8 % với ρ = 0,5 cho G* Nhờ vào các kết này, ta có e thể nói LSMnS làm việc tốt với các mơi trường với hai pha rắn có các mô đun đàn hồi khác ụng cho c c mẫu xi măng trộn đất h ẩn bị mẫ Các mẫu mơ số hóa từ các mẫu trộn phòng thí nghiệm, quy trình chuẩn bị mẫu bao gồm xác định tỉ lệ phù hợp, trộn mẫu thực tế số hóa mẫu để có mẫu đảm bảo đầu vào mô hình mạng lưới lò xo Việc lựa chọn tỷ lệ xi măng với đất ảnh hưởng lớn đến tính chất hỗn hợp vật liệu xi măng đất và giá thành công trình Tỷ lệ xi (a) măng với đất tính theo % khối lượng xi măng so với khối lượng đất khô theo phương pháp thí nghiệm phòng xác định sức kháng nén xi măng đất (TCVN:9403 2012) Tỷ lệ xi măng với đất thích hợp JOMC 80 Tạ chí Vật liệ àX ựng số 21 3 thay đổi theo loại đất và có giá trị biến đổi phạm vi định Theo Thân Văn Vân (2009), lượng xi măng trộn vào là % 15 % trọng lượng khô đất cần gia cố lượng xi măng từ 180 250 kg/m3 đất gia cố Thông thường hàm lượng hạt sét đất yếu tăng thì lượng xi măng yêu cầu tăng Dựa vào thực tế đất khu vực đồng bằng sông Cửu Long chủ yếu là sét và bùn sét yếu (Nguyễn Minh Đức và cộng sự, 2018), nghiên cứu này sử dụng các mẫu với tỉ lệ khác từ %, %, 13 % 18 % và 21 % để đánh giá hiệu phương pháp trộn xi măng đất và ảnh hưởng tỷ lệ xi măng tới cường độ kháng nén hỗn hợp trộn nh Mẫu F21 thực tế và mặt cắt ngang (tỉ lệ lớn hơn), có lẫn thành phần khác: màu xanh là XM, màu đỏ là cát, màu vàng là sét Bảng Thành phần cấp phối các mẫu đất trộn xi măng Mẫu F0 F7 F9 F13 F18 F21 Tỷ lệ xi măng (%) 13 18 21 ết ả mô hỏng Trong mục này, chúng ta áp dụng LSMnS để mô các mẫu với tỉ lệ cấp phối khác trình bày Bảng Sau xử lí bằng kỹ thuật chụp hình ảnh và Direct Samples, chúng ta có Số hóa mẫ mẫu dạng lập phương chia làm 480  480  480 ô mạng lập uy trình xử lý số hóa mẫu từ mẫu trộn thực tế trình bày nghiên cứu Padzniakou (2012) và Nguyễn Thế Anh (2015) thông qua các bước chụp cắt lớp mẫu (m CT), phân tích thông qua công cụ Direct samples (Nguyen K T, 2013) Mẫu chụp cắt lớp nhiều lần, thơng qua xác định phân bố các thành phần lớp cắt, các thành phần rắn khác gán các giá trị 0, 1, 2… từ phương sở Hình 10 và Hình 11 Các mô đun đàn hồi hạt sét và xi măng sử dụng mô lấy theo các kết nghiên cứu Han (1986), Arns (1996), Nguyễn Uyên (2009), Wang và cộng (2015), cụ thể cho Bảng Bảng Đặc tính vật lý các hạt rắn mô K (GPa) G (GPa) ρ (g/cm3) Sét 37 44 2,65 XM 17,5 13,1 2,56 tạo nên mẫu số hóa trình bày mục Trong nghiên cứu Nguyễn Thế Anh (2015), các mẫu chụp m CT với bề dày lớp chụp khoảng 5,7 nm; rõ ràng, công nghệ này đòi hỏi chi phí cao so với giới hạn đề tài, đó, tác giả sử dụng công nghệ chụp ảnh các lớp mẫu, dựa vào màu ảnh và xử lý qua công cụ Direct samples (Nguyen Bảng toán mô hình mạng lưới lò xo áp dụng thành Mẫu Ke (GPa) Ge (GPa) F0 19,62 20,8307 F7 23,56 21,56907376 F9 24,84 25,3881931 F13 28,96 31,25033126 F18 33,32 38,4717207 F21 35,54 39,60876769 K T, 2013) để có các mẫu cần dùng cho mơ Do giới hạn tính phần rắn và tỉ lệ các vật chất khác mẫu bé (Hình 10, Hình 11) nên các mẫu này xem gồm thành phần sét, xi măng và lỗ rỗng mẫu Biện pháp này đơn giản hơn, chi phí thấp nhiên làm giảm độ thực tế mẫu số hóa Kết mơ mơ đun đàn hồi hữu hiệu Ke và Ge mẫu Các kết mô là các thành phần tenso ứng suất hữu hiệu C( eff ) ngoài sử dụng để tính các mô đun đàn hồi hữu hiệu K e nh Mẫu đất thành phần rắn F18 số hóa và mặt cắt ngang nó, màu xanh là XM và Ge thể Bảng và Hình 12 Kết này tương đồng với các kết nghiên cứu Thân Văn Văn (2009), Nguyễn Sỹ Hùng và Vương Hoàng Thạch (2020) cho thấy các tính chất đất JOMC 81 Tạ chí Vật liệ gia cường tăng tỷ lệ hàm lượng xi măng trộn đất àX ựng số 21 3 Kết mô còn thiếu liên hệ với kết thí nghiệm Tại hàm lượng xi măng 13 % (gần tương ứng 340 kg/m ), giá trị Ke thực tế cường độ kháng nén mẫu có chênh lệch tương càng tăng thì mô đun đàn hồi hữu hiệu mẫu càng tăng, nhiên liên hệ tương quan giữa hệ số đàn hồi hữu hiệu và cường độ kháng nén tăng 47 % Ge tăng 50 % so với mẫu không trộn Hàm lượng sau hàm lượng vượt qua 13 % thì mức độ tăng các tính chất giảm lại đối lớn giữa mơ và thực tế Điều này các cơng thức cơng bố có biên thiên quá lớn mô hình chưa xét đến tương tác hóa học giữa xi măng quá trình đơng kết và hạt sét hoăc quá trình số hóa mẫu quá thô không sử dụng công nghệ m CT Điều này là giới hạn bài báo, cần nghiên cứu mở rộng để đánh giá chính xác kết mô so sánh với thực tế ài liệu th m hảo nh Mối quan hệ giữa tỉ lệ xi măng trộn đất và Mô đun đàn hồi hữu hiệu Ke và Ge Các kết liên hệ với kết nén thực tế phòng thí nghiệm [1] C.H Arns, “The influence of morphology on physical properties of reservoir [2] C Boutin and J L Auriault, “Dynamic behaviour of porous media saturated [3] tồn nhiều vấn đề cần giải tính toán sơ cho thấy kết kháng nén mẫu đất qua mô và thực tế chưa đồng với Vấn đề này mơ hình chưa xét đến tương tác [4] hóa học giữa xi măng quá trình đông kết và hạt sét hoăc quá trình số hóa mẫu quá thơ khơng sử dụng công nghệ m CT Điều này là giới hạn đề tài, cần nghiên cứu mở rộng để đánh giá chính xác kết mô so sánh với thực tế ết luận mở rộng [5] [6] [7] thấy phương pháp trộn nông xi măng đất là những biện pháp [8] Việc sử dụng các mô hình số để mô các loại vật liệu và kết cấu trở thành xu hướng và sử dụng rộng rãi Nghiên cứu này đề xuất và phát triển công cụ hữu ích là mơ hình mạng lưới lò xo đàn hồi Mô hình LSM hay mở rộng LSMnS dễ dàng áp dụng cho các mơi trường rỗng với hay nhiều thành phần rắn để xác định các tính học chúng Việc mô các mẫu đất trộn với tỉ lệ xi măng khác cho thấy hàm lượng xi măng càng tăng thì các tính chất đất càng tốt Tại hàm lượng 13 %, giá trị Ke và Ge tăng lần lượt tương ứng là 47 % và 50 %, nhiên sau mức độ tăng giảm dần theo hàm lượng xi măng thêm vào Bài báo đề xuất hàm lượng hợp lý là 13 % đến 15 % Sci., vol 28, no 11, pp 1157–1181, 1990 R Burla, A Kumar, and B Sankar, “Implicit boundary method for determination of effective properties of composite microstructures,” nt Soli s Str ct., vol 46, Mar 2009, doi: 10.1016/j.ijsolstr.2009.02.003 G A Buxton, R Verberg, D Jasnow, and A C Balazs, “Newtonian fluid meets an elastic solid: Coupling lattice Boltzmann and lattice spring [9] h s e , vol 71, no 5, p 056707, May 2005, doi: 10.1103/PhysRevE.71.056707 H Che and M Nguyen, “Amdahl’s law for multithreaded multicore processors,” arallel Distrib om t., vol 10, no 74, pp 3056–3069, 2014, doi: 10.1016/j.jpdc.2014.06.012 I Cohen, “Simple Algebraic Approximations for The Effective Elastic Moduli of a Cubic Array of Spheres,” Mech h s Soli s S L DS, Kết nghiên cứu cho thấy địa hình đất yếu khu vực phù hợp, triển khai rộng rãi và chi phí phù hợp by a viscoelastic fluid Application to bituminous concretes,” nt ng models,” đồng bằng sông Cửu Long là có thành phần chủ yếu là sét bùn sét yếu Nhiều phương pháp xử lý đề xuất và nghiên cứu này cho rocks,” Ph.D.thesis, University of technology, Aachen, 1996 vol 52, pp 10.1016/j.jmps.2004.02.008 2167–2183, Sep M 2004, S doi: Chida, S, “Development of dry jet mixing methods Public Works Research Institute,” Japan: Ministry of Construction, 1982, pp 29–35 A Eggestad, “Improvement of cohesive soils” State of the art report, roc o the th SM , Main Session A, 1983, pp.991 1007 S Kazemian and B Huat, “Assessment of stabilization methods for soft soils by admixtures,” nt on Sci Soc 10.1109/CSSR.2010.5773714 es SS , 2010, doi: [10] D.Han, “Effects of porosity and clay content on acoustic properties of sandstones and unconsolidated sediments,” Ph.D thesis, University of Stanford, 1986 [11] S Horpibulsuk, N Miura, H Koga, and T S Nagaraj, “Analysis of strength development in deep mixing: a field study,” roc nst m ro , vol 8, no 10.1680/grim.2004.8.2.59 2, pp 59–68, i ng Jan 2004, ro n doi: [12] A J C Ladd, J H Kinney, and T M Breunig, “Deformation and failure in cellular materials,” h s e , vol 55, no 3, pp 3271–3275, Mar 1997, doi: 10.1103/PhysRevE.55.3271 JOMC 82 Tạ chí Vật liệ [13] A J C Ladd and J H Kinney, “Elastic constants of cellular structures,” h s Stat Mech ts l., vol 240, no 1, pp 349–360, Jun 1997, doi: đất yếu,” Tạ hí 66–69, Sep 2006 hoa ọc àX ựng số 21 3 ỹ Th ật Thủ Lợi Và Môi Trường, no 26, pp 10.1016/S0378 4371(97)00158 [34] Tsytovich N., Berezantsev V., Dalmatov B and Abelev M, o n ation soils M Eng Thesis, Asian Institute of Technology, Bangkok, Thailand, 1989 [35] TCCS 05:2010/VKHTLVN, “Hướng dẫn sử dụng cọc xi măng đất thi công [14] K H Law, “Strength & deformation characteristics of cement treated clay,” [15] G A Lorenzo and D T Bergado, “Fundamental Characteristics of Cement Admixed Clay in Deep Mixing,” Mater i ng., vol 18, no 2, pp 161– 174, Apr 2006, doi: 10.1061/(ASCE)0899 1561(2006)18:2(161) [16] L Malinouskaya, “Propagation des ondes acoustiques dans les milieux heterogeneous,” Ph.D thesis, Université Pierre et Marie Curie, Paris VI, 2007 [17] N Miura, B T Bergado and A Sakai Improvement of soft marine clays by special admixtures using dry and wet jet mixing method roc th S Bangkok, 1987, 8:46 8:35 [18] M P Moseley and K Kirsch, Eds., , ro n im ro ement, 2nd ed London ; New York: Spon Press, 2004 an s bstr ct res Moscow: Mir publisher, 1974 theo phương pháp Jet grouting để xử lý đất yếu và chống thấm công trình bằng đất.” Viện Khoa học Thủy lợi Việt Nam [36] TCXDVN 385:2006, “Gia cố đất yếu bằng trụ đất xi măng.” [37] S Torquato, “Effective stiffness tensor of composite media : II Applications to isotropic dispersions,” Mech h s Soli s, vol 46, no 8, pp 1411– [38] 1440, Aug 1998, doi: 10.1016/S0022 5096(97)00083 Wang et al., “A new dynamic testing method for elastic, shear modulus and Poisson’s ratio of concrete,” onstr il Mater., vol 100, pp 129– 135, Dec 2015, doi: 10.1016/j.conbuildmat.2015.09.060 [19] S Nemat Nasser, and T H M Iwakuma On composite with periodic structure, Mech.Mater 1, 1982, 239 267 [20] A Nguyen The, “Acoustic wave in porous media Numerical study of wave propagation in porous media with one or many mineral components Applications to real Fontainebleau and STATOIL samples,” Ph.D thesis, Université Pierre et Marie Curie, Paris VI, 2015 [21] K T Nguyen, “Direct samples,” Ph.D thesis, Universit e Pierre et Marie Curie, 2013 [22] Đ Nguyễn Đình, “Nghiên cứu giải pháp cọc cát để gia cố tầng đất yếu khu vực thành phố Hải Phòng.” Đề tài nghiên cứu khoa học trường Đại học dân lập Hải Phòng, 2009 [23] B Nguyễn Ngọc, B Lê Thị Thanh, and P Vũ Đình, Đất x dựng, 2005 [24] U Nguyễn, Xử lý đất ế x [25] D Nguyễn ựng – Địa chất ựng Hà Nội: Nhà xuất Xây cơng trình ỹ th ật cải tạo đất x ựng Hà Nội: NXB Xây dựng, 2009 uốc, A Phùng Vĩnh, and H Nguyễn uốc, ông nghệ hoan hụt cao xử lý đất ế Nhà xuất Nông nghiệp, 2005 [26] Đ Nguyễn Minh, Đ Nguyễn Văn, and A Nguyễn Thế, “Nghiên cứu cường độ đất sét yếu gia cường vải địa kỹ thuật điều kiện thí nghiệm CBR trường,” Tạ hí X Dựng, no 617, pp 192–196, 2019 [27] Đ Nguyễn Minh and L Lê Đức, “Nghiên cứu cải tạo đất bùn nạo vét lòng sông thay cát công tác san lấp mặt bằng tỉnh An Giang,” Tạ X Dựng, pp 59–63, Nov 2018 hí [28] Đ Nguyễn Minh and L Trần Ngọc, “Nghiên cứu ứng xử lún cơng trình móng nơng cát san lấp và đất sét yếu bão hòa tỉnh An Giang,” Tạ hí X Dựng, pp 53–58, 2018 [29] H Nguyễn Sỹ and A Hoàng, “Nghiên cứu thực nghiệm giải pháp gia cố đất yếu có cát san lấp bằng túi đất D BOX,” Tạ hí X Dựng, 2020 [30] H Nguyễn Sỹ and T Vương Hoàng, “Gia cường nơng đất yếu có cát san lấp bằng cọc xi măng đất,” Tạ hí X Dựng, 2020 [31] Pazdniakou, A, “Lattice model in porous media studies,” Ph.D thesis, Université Pierre et Marie, Paris VI, 2012 [32] A Pazdniakou and P M Adler, “Dynamic permeability of porous media by the lattice Boltzmann method,” ater eso r., vol 62, pp 292–302, Dec 2013, doi: 10.1016/j.advwatres.2013.06.001 [33] Thân Văn Văn, “Lựa chọn tỷ lệ xi măng với đất chế tạo cọc xử lý JOMC 83 ... xác định các tính chất học môi trường rỗng trở điều kiện đặc thù khu vực đồng bằng sông Cửu Long Các nghiên rằng, đất điển hình đồng bằng sơng Cửu Long có thành phần chủ yếu là đất. .. kết thấy rằng, phương pháp trộn nông xi măng đất để xử lý đất yếu, tạo ổn định toàn khối là giải pháp phù hợp với khu vực đồng bằng sơng Cửu Long Do đó, đề tài này chất khác có hiệu khác... cho thấy khả và chính xác Mơ h nh mạng lưới lò xo B ic L ttic S ring Mo 2.1 3D LSM, lị xo đơn, lị xo góc l (a) (b) (c) (d) Phương pháp mạng lưới lò xo dùng để mô môi trường đàn hồi, sau