Phân tích ổn định mái đập đá đổ bằng phương pháp phần tử rời rạc

Nghiên cứu này chỉ xét đến ảnh hưởng của trọng lực lên đập, kết quả nghiên cứu thể hiện miền thể tích của những phần tử có chuyển vị lớn, xu hướng bứt phá ra khỏi kết cấu mái đập.

BÀI BÁO KHOA HỌC PHÂN TÍCH ỔN ĐỊNH MÁI ĐẬP ĐÁ ĐỔ BẰNG PHƯƠNG PHÁP PHẦN TỬ RỜI RẠC Nguyễn Thanh Hải1 Tóm tắt: Hiện việc tính tốn hệ số mái cho loại đập đá đổ nhiều hạn chế nghiên cứu Việt Nam Trong báo này, sử dụng phương pháp phần tử rời rạc DEM để phân tích tương tác phần tử đá bên thân đập, ảnh hưởng đến ổn định mái đập Trong mô hình chiều (2D), phần tử đá mơ phần tử hồn tồn cứng có hình dạng đa giác Các hạt không bị biến dạng Nghiên cứu xét đến ảnh hưởng trọng lực lên đập, kết nghiên cứu thể miền thể tích phần tử có chuyển vị lớn, xu hướng bứt phá khỏi kết cấu mái đập Bằng phương pháp DEM, xác định hệ số ổn định mái đập cho phần thể tích, rõ vị trí trượt mái kết cấu đập Kết nghiên cứu thể trường vận tốc phần trượt mái Từ khóa:Vật liệu rời rạc, ổn định mái đập, hệ số ổn định mái, đập đá đổ, phương pháp phần tử rời rạc ĐẶT VẤN ĐỀ* Vấn đề ổn định mái dốc, mái đập đá đổ nghiên cứu rộng rãi nay, cịn nhiều tốn chưa thể giải triệt để Trong bối cảnh đó, ổn định mái dốc cơng trình xây dựng ổn định mái dốc đường, ổn định mái cho đập giữ nước,… quan tâm lớn Chính có nhiều tham số đặc trưng gây ảnh hưởng đến việc ổn định mái đập hệ số mái, hệ số ma sát, đặc tính vật liệu, mực nước ngầm, mưa…(L Nansheng et al., 2015, S Keming, M.R Bagale, 2012 , N.M Pinyol et al., 2008) gây bất lợi cho cơng trình Việc ứng dụng phương pháp số để giải đặt nhiều vấn đề chưa làm rõ Trong phương pháp số thực việc tính tốn ổn định mái dốc, phải kể đến phương pháp cân giới hạn (limit equilibrium method - LEM) sử dụng phổ biến Hệ số ổn định – the safety factor phương pháp định nghĩa tỷ số lực chống trượt lực trượt mặt trượt tới hạn Trong phương pháp số, phương pháp phần tử hữu hạn (Finite element method - FEM) Trường Đại học Bách Khoa, Đại học Đà Nẵng, Việt Nam xem xét xác định hệ số ổn định mái dốc (D.V Griffiths, P.A Lane, 1999) Cung trượt xuất xảy biến dạng lớn liên kết, lúc ứng suất cắt nhỏ, cường độ chịu cắt lớn, kết ổn định mái phụ thuộc lớn đến liệu đầu vào khai báo tính toán Một phương pháp khác nghiên cứu ổn định mái dốc có phương pháp phần tử rời rạc (Discrete element method - DEM)(Y Guan et al., 2017, Y Lu et al., 2018) Khi xét đến ổn định mái dốc, xem kết cấu tổ hợp phần tử rời rạc liên kết với (K.J Chang, A Taboada, 2009) tổ hợp thành phần rời rạc liên tục (Y Lu, 2018) Một số nghiên cứu sử dụng kết hợp phương pháp cân giới hạn phương pháp phần tử rời rạc (DEM)(E Li et al., 2014, W.J Xu et al., 2019, Y Guan et al., 2017) để phân tích ổn định mái dốc nhận kết khả quan so sánh độc lập với LEM hay DEM Công trình nghiên cứu (W Zheng et al., 2014) phân tích mối tương quan tính liên tục khơng liên tục để đánh giá ổn định mái dốc đồ thị Đặc điểm phương pháp phần tử rời rạc có xét đến chi tiết bên vật liệu, liên kết phần tử cấu thành vật liệu Điều cho phép nghiên cứu phạm vi cấu trúc vật KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 67 (12/2019) 71 liệu (F Radjai, F Dubois 2011,F Radjai and V Richefeu 2009 ,S Nezamabadi et al 2017) Do vậy, việc xem xét vật liệu bị phá hủy phần tử liên kết tạm thời hay liên kết bị phá bỏ hồn tồn Chính điều gây trượt kết cấu hay phá hủy phần kết cấu Sự phản ánh chi tiết diễn biến liên kết phần tử phương pháp phần tử rời rạc cho phép rõ q trình, vị trí liên kết hay phá hủy liên kết Trong báo này, sử dụng phương pháp phần tử rời rạc để mô cho trường hợp ổn định mái dốc đập đá đổ cấu thành từ phần tử đa giác (xem mục 2, mục 3) Mô kết cấu đập để có phân tích ổn định mái đập, thể trường vận tốc (xem mục 4) Chúng tơi có kết luận mục với kết (a) phản ánh ổn định mái dốc có định hướng cho nghiên cứu PHƯƠNG PHÁP SỐ Phương pháp phần tử rời rạc DEM cho phép mơ tốn môi trường phần tử rời rạc đá, cát, bột,… (F Radjai, F Dubois 2011,F Radjai and V Richefeu 2009) Trong phương pháp này, phương trình chuyển động phần tử hồn tồn cứng tích hợp cách tính đến ràng buộc động học tương tác tiếp xúc Những tương tác đặc trưng thông số: hệ số ma sát, hệ số đàn hồi va chạm, lực dính phần tử Tại tương tác xác định giá trị lực pháp tuyến có xét đến hệ số ma sát, lực dính (T.H Nguyen et al 2017,S Nezamabadi et al 2017,K Krabbenhoft et al 2012) (b) Hình Hai loại liên kết phần tử đa giác (a) liên kết điểm - cạnh (bên trái); cạnh-cạnh (bên phải); (b) xác định véc tơ đơn vị lực liên kết Xét hệ thống nhiều phần tử, tương tác hệ thống mô tả đại lượng vô hướng thể khả liên kết phần tử hình Lúc hệ thống tồn mạng lưới tiếp xúc tương tác phần tử Các điểm tiếp xúc xem hoạt động điểm tiếp xúc nhận giá trị lực khác không Ngược lại, điểm tiếp xúc không hoạt động hai phần tử không tồn tương tác lực, khơng góp phần vào việc truyền ràng buộc Trong nghiên cứu này, quan tâm đến phần tử có hình dạng đa giác, xãy tương tác phần tử quan niệm tính hình a,b Xác định liên kết hai phần tử đa giác, liên kết đơn giản điểm cạnh lực pháp tuyến vng góc với cạnh Ngược lại liên kết cạnh - cạnh, lúc xuất véc tơ đơn vị giá trị lực pháp tuyến 72 Trong trường hợp này, hai giá trị quy đổi thành giá trị đại diện cho mối liên kết (D.H Nguyen et al 2018) Ngoài ra, trường hợp liên kết hai phần tử điểm - điểm, trường hợp xảy Nếu trường hợp tiếp xúc điểm - điểm tồn tại, (D.H Nguyen et al 2015, 2018) đề xuất xem chúng liên kết đơn giản phần tử liên kết điểm - cạnh cạnh - cạnh Hình giới thiệu hệ thống lực tập hợp phần tử đa giác cứng Lực tương tác phần tử thể đoạn thẳng nối phần tử Độ đậm, nhạt đoạn thẳng tỷ lệ thuận với giá trị độ lớn lực tương tác Đặc trưng cường độ cắt (shear strength) vật liệu môi trường rời rạc góc ma sát vật liệu tensơ ứng suất giai đoạn biến dạng phải tính tốn từ hệ KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 67 (12/2019) thống liên kết giá trị lực Để xác định ten sơ ứng suất, sử dụng mômen lực Mi phần tử (particle i)(J J Moreau, 1997): định từ giá trị trung bình q/p trạng thái định vật liệu (J Mitchell and K Soga, 2005) đại diện cho cường độ cắt nội vật liệu (1) Trong đó, vectơ lực tác dụng lên phần tử i điểm tiếp xúc c, khoảng cách từ tâm phần tử i đến điểm đặt lực Trường hợp xuất nhiều điểm tiếp xúc phần tử i tính tổng tất điểm tiếp xúc c phần tử lân cận với phần tử i Như mơmen lực tập hợp phần tử hệ thống tổng số mômen lực phần tử thành phần (3) Xét quan hệ ứng suất biến dạng Coulomb tốn ứng suất phẳng ta xác định quan hệ: (4) Trong thành phần ứng suất tiếp ứng suất pháp, C lực dính kết đơn vị phần tử Xét cho phương hợp với mặt phẳng ngang góc thành phần ứng suất tiếp ứng suất pháp thể sau: (5) (6) Hình Hệ thống lực tập hợp phần tử đa giác cứng Độ đậm đường phần tử tỷ lệ thuận với độ lớn lực Xét trong hệ thống có n phần tử miền thể tích V, tensơ ứng suất xác định theo công thức(J J Moreau, 1997): (2) Trong đó, khoảng cách nối tâm phần tử tiếp xúc Từ đây, xác định ứng suất trung bình Như vậy, áp dụng phương pháp phần tử rời rạc môi trường vật liệu rời cho phép xác định thành phần ứng suất, thành phần lực chi tiết bên vật liệu Đây điều kiện để vận dụng phương pháp phần tử rời rạc phân tích ổn định mái dốc, taluy, mái đập,… sử dụng vật liệu rời ÁP DỤNG TÍNH TỐN ỔN ĐỊNH MÁI ĐẬP ĐÁ ĐỔ Trong báo thực mô hình tốn đập đá đổ gồm phần tử đa giác dựa tỷ lệ cấp phối hàm lượng hạt bảng theo TCVN10777:2015 “Cơng trình thủy lợi - Đập đá đổ mặt bê tông – yêu cầu thiết kế” Bảng Tỷ lệ hàm lượng hạt mơ hình STT Hàm lượng hạt 5% 20% 75% Kích thước hạt d = [0.05 – 0.075]mm d = [0.075 – 5]mm d = [5 – 800]mm ứng suất lệch , ứng suất tốn phẳng chiều Góc ma sát φ vật liệu rời xác Phần tử đá mô mơ hình hạt có hình dạng đa giác với số cạnh từ 5-8 cạnh Phần tử hạt có hình dạng đơn giản đa giác đều, kích thước khác Các phần tử đá KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 67 (12/2019) 73 mô báo hạt hồn tồn cứng, khơng bị biến dạng, có khối lượng riêng 2,5(g/cm3) Xét đến ma sát phần tử, mơ hình chúng tơi mặc định hệ số ma sát  phần tử 0,5 Góc ma sát vật liệu xác định theo công thức (3) cho vật liệu rời rạc Đối với vật liệu phần tử đá có dạng đa giác, D.H Nguyen et al.,chỉ rõ tỷ số q/p trạng thái tới hạn nằm phạm vi từ 0,4 đến 0,5, có nghĩa góc ma sát từ 23° đến 30°(D.H Nguyen, 2015) Giá trị tương đối hợp lý cho trường hợp vật liệu đá rời rạc Trong mơ hình chúng tơi đề xuất nghiên cứu cho trường hợp góc phi 25° Gia tốc trọng trường mơ hình chọn 9,81 (m/s2) Biên trái, biên phải đáy nhóm (cluster) cấu thành từ phần tử cứng, khơng biến dạng, khơng có chuyển vị Mặt đập mặt thống (a) (b) Hình Mơ hình mơ có hình dạng mái đập đá đổ có hệ số mái m = 1,2: (a) Kích thước chi tiết phần kết cấu đập; (b) Mơ hình kết cấu đập cấu thành từ phần tử đá Như hình giới thiệu mơ hình mái đập có hệ số mái m = 1,2, có vật liệu phần tử đá hồn tồn cứng hình 2b Hình 2a thể chi tiết kích thước chiều dài, cao, rộng cạnh Độ rỗng thân đập xét chiếm khoảng 17% thể tích kết cấu đập Một kết cấu mái dốc đảm bảo ổn định có hệ số K tỷ lệ thành phần chống trượt thành phần gây trượt không nhỏ hình Hình Miền thể tích V chứa tập hợp phần tử phần thể tích thành phần Ai có khả ổn định ổn định toàn mái đập Trong báo này, xác định phạm vi việc liên kết giới hạn tổng số phần tử có 74 khả trượt khỏi kết cấu ban đầu so với tổng số phần tử lưu lại kết cấu đập Xét toán tập hợp phần tử miền thể tích V có xu hướng liên kết so với phần thể tích cịn lại kết cấu mái đập, lúc từ phương trình (4) hệ số ổn định K phần thể tích V xác định bởi: (7) Nếu K < mái dốc có khả xảy trượt mái miền thể tích V, ngược lại K ≥ mái đập đảm bảo ổn định Khi K lớn khả ổn định lớn KẾT QUẢ VÀ BÌNH LUẬN Sử dụng ngơn ngữ lập trình C++, thực mơ hình số cho tốn ổn định mái đập tác dụng trọng lực, xét trường hợp không chịu ảnh hưởng lực mao dẫn, áp lực nước, Do vậy, phần tử kết cấu đập tương tác với nhau, trình dịch chuyển, chuyển vị gây ổn định mái KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 67 (12/2019) Safety Factor 2.0 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 Hình Mơ ổn định kết cấu đập thực mô thời điểm 10.000 bước tính tốn Hình thể khối trượt hay cung trượt trình mơ sau 10.000 bước tính tốn Kết cho thấy có khối trượt, xét từ phải sang trái, độ đậm màu giảm dần Để xác định giá trị Ki, miền thể tích Vi có tập hợp phần tử có xu hướng liên kết với kết cấu, tính tốn cơng thức: ổn định K = [0,8 - 1,0] (phần tử miền V2 tổng phần tử phần thể tích A1 A2 ) Miền thể tích V3 cấu thành từ phần tử phần thể tích A1 , A2 A3 có hệ số ổn định K = [1,0 - 1,2] Miền thể tích V4 tập hợp phần tử cịn lại thân đập có gam màu sáng Như vậy, với kết nhận phần tử vùng V1 ,V2 có khả trượt khỏi mặt đập Ngược lại xét phần tử thuộc từ miền V3 hay V4 đảm bảo ổn định Xảy khả trượt V1 ,V2 phần tử miền V1 chuyển vị lớn tác động trượt phần tử gần nó, kéo theo phần tử miền thể tích V2 trượt khỏi mặt đập đá đổ Tuy nhiên tác dụng trọng lực, hệ số ma sát phần tử, độ gồ ghề bề mặt phần tử nên làm giảm khả trượt hay tăng khả chống trượt xét miền V3 V4 (8) Trong nghiên cứu này, mô ổn định đập cấu thành từ phần tử đá, nên giá trị lực dính phần tử đá giả thiết nhỏ, bỏ qua Tập hợp phần tử miền Vi định nghĩa tập hợp phần tử có chuyển vị phạm vi gần nhau, chênh lệch phần tử có chuyển vị lớn với nhỏ miền Vi không vượt 1,5 lần Do có việc xuất nhiều Vi cho phép thấy rõ vị trí trượt cục Lưu ý đây, tập hợp phần tử tính tốn miền thể tích V phạm vi cấu thành từ phần tử thể tích thành phần A1, 2,…,n tức V = A1+A2+…+An Từ kết hình 4, nhận thấy có vùng màu khác nhau, định nghĩa miền diện tích V1 , V2 , V3, V4 từ phải sang trái hay mức độ thang màu giảm dần từ đậm sang nhạt Trong miền V1 (ngoài cùng, bên phải) có hệ số ổn định K = [0,6 - 0,8] tập hợp phần tử vùng màu đậm (phần tử miền V1 phần tử phần thể tích A1 ) Tiếp miền V2 vùng màu sáng V1 có hệ số Hình Trường vận tốc phần tử mơ hình thời điểm 10.000 bước tính tốn Chúng ta nhận thấy trường vận tốc phần tử mơ hình hình 5, phản ánh kết hình Khi phần tử phía ngồi (bên phải) V1 V2 nhận giá trị lớn, giảm dần tiến vào lõi đập hay phía trái đập Như trình bày mục 2, phần tử đá không tiếp xúc hay nói cách khác khơng có tương tác, lúc khơng xuất lực phần tử Do vậy, mái đập có ổn định, xảy q trình trượt mái liên kết phần tử khối thể tích kết cấu mái Lúc này, phần thể tích chứa phần tử đá có xu hướng trượt phần thể tích khác tác dụng trọng lực hay nói cách khác chuyển vị phần tử tăng lớn Phần thể tích lại xem ổn định, trường vận tốc gần nhỏ KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 67 (12/2019) 75 KẾT LUẬN Trong báo tác giả sử dụng phương pháp phần tử rời rạc - DEM để mô mối liên kết phần tử đa giác lập trình ngơn ngữ C++, đại diện cho phần tử đá thân đập đá đổ, xem xét ảnh hưởng đặc trưng đến việc ổn đỉnh mái Mái đập ổn định ảnh hưởng trọng lực ứng với phần thể tích trượt khác với hệ số mái đập đá đổ 1,2 Các miền thể tích chứa phần tử rời ổn định, trượt bề mặt mái đập V1, V2 liên kết với khối thể tích cịn lại đập Sự chuyển vị phần tử làm cho phân bố ứng suất mái đập thay đổi hay giảm xuống liên kết, giảm tương tác với khối phần tử Miền thể tích V3 thời điểm tính tốn xác định ổn định K = [1,0 - 1,2], nhiên miền có khả ổn định phần tử miền thể tích V2 trượt xuống, kéo theo phần tử tiếp giáp A2 A3 thời gian sau Trường vận tốc phần tử thể rõ, cho thấy chuyển vị phần tử ổn định, trượt khỏi kết cấu mái Ưu điểm mơ hình phân tích ổn định DEM nghiên cứu chi tiết nội bên vật liệu rời như: (i) độ rỗng, độ chặt kết cấu; (ii) hình dạng xếp phần tử; (iii) xác định liên kết hay phá hủy liên kết phần tử bên kết cấu So với phương pháp phân tích ổn định phổ biến nay, mơ hình báo mơ tả trình ổn định từ bắt đầu có chuyển vị phần tử hay liên kết phần tử kết cấu Xác định vị trí xác liên kết kết cấu để có biện pháp khắc phục, hạn chế ổn định trượt mái đập Hơn nữa, sử dụng mơ hình phân tích, đề xuất hình dạng, kích thước, cấp phối đá nhằm giảm độ rỗng, tăng độ chặt kết cấu Mơ hình cịn tồn vấn đề chưa xem xét ảnh hưởng mao dẫn, bỏ qua giá trị lực dính, đa dạng hình dạng hình học phần tử, số lượng phần tử cịn Mơ hình chưa xét ảnh hưởng nước ngầm dòng thấm Bên cạnh đó, hạn chế phương pháp DEM khối lượng tính tốn lớn, tốn nhiều thời gian số lượng phần tử mơ hình tăng lên TÀI LIỆU THAM KHẢO D.H Nguyen, E.Azéma, F.Radjai, and P.Sornay (2014), “Effect of size polydispersity versus particle shape in dense granular media”, Phys Rev E 90, 012202 D.H Nguyen, E Azema, P Sornay, and F Radjai (2015),“Effects of shape and size polydispersity on strength properties of granular materials”, Phys Rev E 91, 032203 D.H Nguyen, E Azéma, P Sornay, F Radjai (2018), “Rheology of granular materials composed of crushable particles”, Eur Phys J E, 41: 50 D.V Griffiths, P.A Lane (1999), “Slope stability analysis by finite elements”, Géotechnique, V49.3, 387-403 E Li, X.Zhuang, W.Zheng, Y Cai (2014), “Effect of graph generation on slope stability analysis based on graph theory” Journal of Rock Mechanics and Geotechnical Engineering, 6, 380–386 F Radjai and V Richefeu (2009), “Contact dynamics as a nonsmooth discrete element method”, Mechanics of Materials, 41:715–728, 2.1 F Radjai, F Dubois (2011), “Discrete-element modeling of granular materials”, Wiley-ISTE J Mitchell and K Soga (2005), “Fundamentals of soil behavior”, Wiley, New York J J Moreau (1997), “Numerical investigation of shear zones in granular materials”, HLRZ-Workshop on Friction, Arching, Contact Dynamics, Singapore pp.233-247 K Krabbenhoft, A.V Lyamin, J Huang, M Vicente da Silva (2012), “Granular contact dynamics using mathematical programming methods”, Computers and Geotechnics, 43,165-176 76 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 67 (12/2019) K.J Chang, A Taboada, (2009), “Discrete element simulation of the Jiufengershan rock-and-soil avalanche triggered by the 1999 Chi-Chi earthquake, Taiwan”, Journal of geophysical research, VOL 114, F03003 L.Nansheng, T Bo, and X Lihui (2015), “Slope stability analysis of earth-rockfill Dams using MGA and UST”, Journal of Computational Engineering, Volume 2015, Article ID 895142, 10 pages N.M Pinyol, E.E Alonso, S Olivella (2008), “Rapid drawdown in slopes and embankments”, Water resources research, V 44, W00D03 S.Keming, M.R.Bagale (2012), “Study on dam slope stability under the condition of rainfall”, International journal of scientific & Technology research, Volum 1, Issue S Nezamabadi, T.H Nguyen, J.Y Delenne, and F Radjai (2017), “Modeling soft granular materials”, Granular Matter, 19.8 T.H Nguyen, S Nezamabadi, J.Y Delenne, and F Radjai (2017), “Compaction of granular materials com- posed of deformable particles”, EPJ Web of Conferences, 140:05013 W.J Xu, S Wang, M Bilal (2019), “LEM-DEM coupling for slope stability analysis”, Sci China Tech, 62, https://doi.org/10.1007/s11431-018-9387-2 W Zheng, X.Zhuang, D.D.Tannant,Y.Cai, S.Nunoo (2014), “Unified continuum/discontinuum modeling framework for slope stability assessment”, Engineering Geology, 179, 90–101 Y Guan, X Liu, E Wang, and S Wang (2017), “The stability analysis method of the cohesive granular slope on the basis of graphtheory”, Materials,10, 240; doi:10.3390/ma10030240 Y Lu, Y Tan, X Li (2018), “Stability analyses on slopes of clay-rock mixtures using discrete element method”, Engineering Geology, V 244, 116-124 Abstract: STUDY ON SLOPE STABILITY OF ROCKFILL DAMS BY DISCRETE ELEMENT MOTHOD (DEM) The calculation of slope stability of rockfill dams is still limited in studies in Vietnam Rockfill dams are being constructed using experience studies or similar constructions In this paper, we use a discrete element method DEM to analyze the interactions between the rock elements within the dam body in twodimension Rock elements are simulated by rigid particles that have a polygonal shape These rock elements are not deformed In this paper, oneonly takes into account the influence of gravity on the dam We determined the safety factor of slope and the velocity field by DEM method Keywords: Granular material, slope stability, safety factor, Rockfill dam, discrete element method Ngày nhận bài: 10/10/2019 Ngày chấp nhận đăng: 30/12/2019 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 67 (12/2019) 77 ... này, sử dụng phương pháp phần tử rời rạc để mô cho trường hợp ổn định mái dốc đập đá đổ cấu thành từ phần tử đa giác (xem mục 2, mục 3) Mơ kết cấu đập để có phân tích ổn định mái đập, thể trường... phương pháp phần tử rời rạc môi trường vật liệu rời cho phép xác định thành phần ứng suất, thành phần lực chi tiết bên vật liệu Đây điều kiện để vận dụng phương pháp phần tử rời rạc phân tích ổn định. .. khơng xuất lực phần tử Do vậy, mái đập có ổn định, xảy q trình trượt mái liên kết phần tử khối thể tích kết cấu mái Lúc này, phần thể tích chứa phần tử đá có xu hướng trượt phần thể tích khác tác