CT 1 I. ĐẶT VẤN ĐỀ Trong bài báo này, liên kết dị hướng được hiểu là liên kết có phản lực liên kết thay đổi theo chiều cũng như độ lớn của chuyển vị của điểm liên kết. Hình 1 mô tả sự làm việc của cọc trong nền đất. Để tính toán nội lực và chuyển vị của cọc tiến hành rời rạc hoá cọc thành n phần tử thanh, qui đổi áp lực pháp tuyến thành liên kết dị hướng theo phương vuông góc với cọc tại các nút và lực ma sát thành bên cọc thành các liên kết dị hướng bên trái và bên phải cọc tại các nút (hình 1b). Áp lực pháp tuyến (kí hiệu N) quan hệ với chuyển vị ngang là một đường cong trên hình 1c, rời rạc hoá quan hệ này thành các đoạn thẳng có các độ cứng k i (trên hình rời rạc hoá thành 2 TÍNH CỌC ĐƠN TRONG NỀN ĐẤT VỚI MÔ HÌNH LIÊN KẾT DỊ HƯỚNG TS. LƯƠNG XUÂN BÍNH ThS. ĐỖ XUÂN QUÝ Bộ môn Sức bền Vật liệu Khoa Công trình Trường Đại học Giao thông Vận tải Tóm tắt: Cho đến nay, khi phân tích ứng xử của cọc, chúng ta thường mô hình hoá liên kết giữa cọc và nền đất bởi các lò xo đàn hồi đẳng hướng theo hai chiều (ứng xử của lò xo theo chiều kéo và chiều nén như nhau). Tuy nhiên, khi cọc bị biến dạng, áp lực pháp tuyến và lực ma sát của nền đất tác dụng lên mặt bên cọc là rất phức tạp. Ứng xử của liên kết giữa cọc và nền là phi tuyến. Trên những vùng cọc không tiếp xúc với nền, liên kết giữa cọc và nền không làm việc. Ngoài ra chuyển vị theo phương thẳng đứng của cọc cũng ảnh hưởng đến lực ma sát trên thành bên. Do đó, việc mô hình hoá liên kết giữa cọc và nền bởi các liên kết đàn hồi đẳng hướng theo hai chiều là không thoả đáng. Trong bài báo này, các tác giả đề xuất mô hình liên kết dị hướng (ứng xử của liên kết là phi tuyến) cho cả áp lực pháp tuyến và lực ma sát giữa cọc và nền đất, đồng thời có xét đến ảnh hưởng giữa chuyển vị dọc của cọc đến lực ma sát, từ đó xây dựng thuật toán và chương trình tính cọc trên máy tính bằng phương pháp phần tử hữu hạn (PP PTHH). Summary: Up to now, when analyzing the response of the pile, we usually modelize the interaction between the pile and soil media with the elastic springs, of which the response in both compression and tension directions are the same. However, when the pile is deformed, the normal pressure and the friction of soil media acting on the surface of pile are quite complicated. The pile-soil interaction is non-linear. On the region, which the pile does not contact the soil media, the connection between pile and soil is not effective. Besides, the longitudinal displacement also affects the lateral friction. Therefore, modeling the interaction between pile and soil as isotropic springs is not adequate. In this paper, the authors would like to propose anisotropic interaction (pile-soil response is non-linear) for both normal pressure and friction between soil and pile, simultaneously, and take the effect of longitudinal displacement on the friction into account. Then, the algorism and computer program to analyze the single pile with the finite element method is built. T CT 1 đoạn). Mô hình hoá áp lực pháp tuyến thành liên kết theo phương ngang cọc loại liên kết nhiều lần tuyến tính. Lực ma sát bên thành cọc (kí hiệu F ms ) quan hệ với chuyển vị dọc của cọc là đường cong trên hình 1d, rời rạc hoá thành các đoạn thẳng có độ cứng k ’ i . Mô hình hoá lực ma sát này thành liên kết theo phương dọc của cọc loại liên kết nhiều lần tuyến tính. Lực ma sát thành cọc là một hàm của chuyển vị dọc và độ cứng k ’ i , trong khi đó độ cứng này lại phụ thuộc vào áp lực ngang của đất vào thành cọc (áp lực ngang mà lớn thì hệ số này cũng lớn). Áp lực ngang cọc lại phụ thuộc vào chuyển vị ngang cọc. Như vậy lực ma sát thành cọc là một hàm của chuyển vị ngang và chuyển vị dọc. Như vậy: ),(fF )(fN docngangms ngang (1) Nên phương trình cơ bản của bài toán trở thành phương trình phi tuyến. a) N M m F ms H M P b) H M P N nga c) O k 1 k 2 N 0 k 3 a1 a2 a3 N a 3 N a2 N a1 F ms d d) O k k F a1 F a2 F a3 F a4 a a a a k Mô hình làm vi ệc CT 1 Hình 1. Cọc và mô hình tính II. MÔ HÌNH TOÁN CỦA LIÊN KẾT Áp lực ngang bên thành trái cọc (theo hình 1c): nan n j jajaja kkNN )()( 2 111 (2) vì khi chuyển vị bằng 0 áp lực ngang N = N 0 nên 2211201 )( kkNN aaaa (3) Áp lực ngang bên thành phải cọc cũng được tính theo công thức như đối với áp lực ngang thành trái cọc, với chú ý là: ' và chiều của áp lực N’ ngược chiều áp lực N. Lực ma sát thành trái cọc (theo hình 1d) : nan n j jajajams kkFF ')(')( 2 111 (4) vì khi chuyển vị dọc bằng 0 lực ma sát cũng bằng 0 nên 112221 ')(' kkF aaaa (5) trong đó các hệ số k’ i phụ thuộc vào áp lực ngang của cọc Nak ii ' (6) với a i là các hằng số (a 1 = a 3 ). Lực ma sát thành phải cọc: tính tương tự như lực ma sát ở thành trái cọc. Ngoài áp lực ngang, lực ma sát hai bên trái và phải cọc thì việc dời lực ma sát hai bên thành cọc về trục thành còn sinh ra hai mô men uốn ngược chiều nhau. Mô men tổng hợp được tính bằng: 2 )'( d FFM msms (7) trong đó d là kích thước của cọc. III. PHƯƠNG TRÌNH CƠ BẢN CỦA PP PTHH CHO KẾT CẤU CÓ LIÊN KẾT DỊ HƯỚNG Phương trình cơ bản của phương pháp phần tử hữu hạn cho kết cấu thông thường là hệ phương trình đại số tuyến tính được viết dưới dạng như sau: PK (8) trong đó: K - ma trận độ cứng của kết cấu, - véc tơ chuyển vị nút của kết cấu, P - véc tơ tải trọng nút. Với kết cấu có liên kết dị hướng, véc tơ tải trọng nút ngoài các tải trọng tác dụng trên kết cấu còn có các phản lực của các liên kết dị hướng (ta coi các phản lực ứng với các liên kết dị hướng như là các ngoại lực). Tách véc tơ tải trọng thành hai véc tơ: 1) véc tơ tải trọng gồm các T CT 1 ngoại lực tại các nút như véc tơ P trong (6); và 2) véc phản lực liên kết dị hướng kí hiệu là N . Khi đó phương trình cơ bản của phương pháp phần tử hữu hạn cho kết cấu có liên kết dị hướng có dạng như sau: PNK (9) trong đó các thành phần của véc tơ phản lực liên kết N được xác định theo (2) (4). Cần chú ý rằng véc tơ phản lực liên kết N là hàm phi tuyến của các chuyển vị nút , nên (9) là hệ phương trình phi tuyến. Để tìm được lời giải cho kết cấu có liên kết dị hướng, ta phải giải hệ phương trình phi tuyến. Đó chính là sự khác biệt giữa bài toán tính kết cấu dị hướng và kết cấu thông thường bằng phương pháp phần tử hữu hạn. IV. GIẢI HỆ PHƯƠNG TRÌNH PHI TUYẾN BẰNG PHƯƠNG PHÁP NEWTON Hệ phương trình (9) được viết dưới dạng khai triển như sau: j iijiji f 0 (10) trong đó f i là hàm số của n biến i . Khai triển Taylor với hàm một biến f(x) tại x = a có dạng: )( ! 2 )( ))(()()( 2 '' ' ax af axafafxf (11) Giá trị hàm số tại x = a + x được tính: ! 2 )( )()()( 2 '' ' x af xafafxaf (12) Nếu như x đủ nhỏ, ta có thể bỏ qua các vô cùng bé bậc cao kết hợp với điều kiện f(x) = 0, ta được: 0)()()( ' xafafxaf (13) Nếu f là một hàm nhiều biến f = f(x) = f(x 1 ,x 2 , , x n ), tại x 0 = T n xxx 00 2 0 1 , ,, công thức (13) trở thành: 0 )( )()( )( 0 2 2 0 1 1 0 0 n n x x xf x x xf x x xf xf (14) Kết hợp các phương trình (10) và (14) với chú ý rằng x 0 = 0 ta thu được hệ phương trình sau: CT 1 )( )( )()( )( )( )()( )( )( )()( )( )( )()( 0 0 2 2 0 1 1 0 0 3 0 3 2 2 0 3 1 1 0 3 0 2 0 2 2 2 0 2 1 1 0 2 0 1 0 1 2 2 0 1 1 1 0 1 nn n nnn n n n n n n f fff f fff f fff f fff (15) Do đó, ta có trình tự các bước giải lặp như sau. Bước 1. Cho véc tơ chuyển vị nút một giá trị ban đầu. 0 ii có thể lấy 0 i Bước 2. Tính giá trị các hàm số f i và các đạo hàm riêng phần của chúng, j i f )( 0 tại 0 . Giải hệ phương trình (15) tìm được véc tơ số gia chuyển vị nút, i . Bước 3. Tính lại véc tơ chuyển vị nút i 1 = i 0 + i . Bước 4. Kiểm tra điều kiện dừng chương trình. Sai số: )( 1 2 i f Kiểm tra sai số: Nếu thì dùng chương trình lấy nghiệm 1 ii . Nếu tiếp tục thực hiện lặp với 10 ii . Nếu số lần lặp nhiều quá thì cho dừng chương trình, lúc này phép lặp không hội tụ. V. VÍ DỤ TÍNH TOÁN VÀ ĐÁNH GIÁ KẾT QUẢ Tính cọc dài 30 m làm bằng vật liệu có E = 4000 kN/cm 2 , mặt cắt ngang cọc có J = 213333 cm 4 , F = 1600 cm 2 , kích thước cọc d = 40 cm. Cọc chịu tác dụng của tải trọng ngang H = 100 kN và tải trọng đứng P = 1000 kN. Rời rạc cọc thành 30 phần tử, mỗi phần tử có chiều dài 1m. Liên kết của cọc với đất: liên kết ngang: k 1 = 10 kN/cm, k 2 = 100 kN/cm, a1 = -10 cm, a2 = -4 cm, N 0 = -100 kN; liên kết dọc: a 1 = 10 1/cm, a 2 = 100 1/cm, a1 = -10 cm, a2 = -5 cm ( k’ = a.N). Kết quả tính toán được thể hiện trong hình 2 và bảng 1. Tính toán cọc theo mô hình liên kết dị hướng cho kết quả khác so với mô hình liên kết Winkler, trong ví dụ này sự sai khác được thể hiện theo bảng 1. Bảng 1. So sánh kết quả với mô hình nền Winkler Mô hình liên kết Chuyển vị dọc lớn nhất (cm) Mô men lớn nhất (kNcm) Độ võng lớn nhất (cm) Winker 0.0308605 12108.377 1.2873678 Dị hướng 0.0368346 14713.806 1.0119647 T CT 1 Sai lệch (%) 19 21.6 21 VI. KẾT LUẬN VÀ KIẾN NGHỊ Các tác giả đã xây dựng được cơ sở lý thuyết, thuật toán cũng như chương trình tự động hóa tính toán cọc có liên kết dị hướng bằng pháp phần tử hữu hạn. Phân tích so sánh cho thấy mức độ sai số khi không xét đến sự làm việc của liên kết dị hướng là đáng kể. Do đó, kiến nghị đưa mô hình liên kết dị hướng vào tính toán cọc trong nền đất. Hình 2. Biểu đồ nội lực và chuyển vị của cọc Tài liệu tham khảo [1] I. M. Rabinovich: "Giáo trình Cơ học kết cấu", Gosstroiizđat Matxcơva, 1954. Biểu đồ lực dọc Biểu đồ chuyển vị dọc trục Biểu đồ mô men Biểu đồ độ võng CT 1 [2] Nguyễn Văn Hợi, Cao Chu Quang: "Tính công trình ngầm có xét đến liên kết tiếp xúc một chiều giữa kết cấu và môi trường đất đá theo phương pháp quy hoạch toàn phương", Tuyển tập công trình Hội nghị Khoa học Toàn quốc Cơ học vật rắn biến dạng lần thứ 7, Tập 1, Nhà xuất bản ĐHQG Hà Nội, 2004. [3] Vũ Đình Lai, Nguyễn Xuân Lựu: "Tính kết cấu có liên kết dị hướng", Tuyển tập công trình Hội nghị Khoa học Toàn quốc Cơ học vật rắn biến dạng lần thứ 8, Tập 1, Nhà xuất bản Khoa học Tự nhiên và Công nghệ, 2006. [4] Lương Xuân Bính, Nguyễn Xuân Lựu, Đỗ Xuân Quý: “Tính kết cấu có liên kết dị hướng bằng phương pháp phần tử hữu hạn”, Tuyển tập Công trình Hội Nghị Cơ học Toàn quốc lần thứ 8, Nhà xuất bản Khoa học Tự nhiên và Công nghệ, 2007 . Trong bài báo này, liên kết dị hướng được hiểu là liên kết có phản lực liên kết thay đổi theo chiều cũng như độ lớn của chuyển vị của điểm liên kết. Hình 1 mô tả sự làm việc của cọc trong nền. -5 cm ( k’ = a.N). Kết quả tính toán được thể hiện trong hình 2 và bảng 1. Tính toán cọc theo mô hình liên kết dị hướng cho kết quả khác so với mô hình liên kết Winkler, trong ví dụ này sự. lực ma sát của nền đất tác dụng lên mặt bên cọc là rất phức tạp. Ứng xử của liên kết giữa cọc và nền là phi tuyến. Trên những vùng cọc không tiếp xúc với nền, liên kết giữa cọc và nền không làm