Nghiên cứu điều kiện địa kỹ thuật phục vụ thiết kế và thi công hố đào sâu

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ XÂY DỰNG VIỆN KHOA HỌC CÔNG NGHỆ XÂY DỰNG NGUYỄN TRƯỜNG HUY NGHIÊN CỨU ĐIỀU KIỆN ĐỊA KỸ THUẬT PHỤC VỤ THIẾT KẾ VÀ THI CÔNG HỐ ĐÀO SÂU LUẬN ÁN TIẾN SĨ CHUYÊN NGÀNH: KỸ THUẬT XÂY DỰNG CƠNG TRÌNH NGẦM MÃ SỐ: 62.58.02.04 NGƯỜI HƯỚNG DẪN KHOA HỌC: TS TRỊNH VIỆT CƯỜNG PGS TS ĐOÀN THỀ TƯỜNG HÀ NỘI - 2015 i LỜI CAM ĐOAN Tôi xin cam đoan cơng trình khoa học tơi thực Các kết quả, số liệu luận án trung thực chưa công bố cơng trình khoa học khác Tác giả hồn tồn chịu trách nhiệm tính xác thực nguyên luận án Tác giả luận án Nguyễn Trường Huy ii LỜI CẢM ƠN Tác giả luận án xin bày tỏ lòng biết ơn sâu sắc trân trọng đến TS Trịnh Việt Cường PGS TS Đoàn Thế Tường hai thầy hướng dẫn trực tiếp tận tình bảo, hướng dẫn giúp đỡ tác giả suốt q trình thực hồn thành luận án Tác giả xin chân thành cảm ơn đến Lãnh đạo Viện Khoa học công nghệ xây dựng, Viện Thông tin Đào tạo Tiêu chuẩn hóa, Viện chuyên ngành Địa kỹ thuật tạo điều kiện thuận lợi, đóng góp ý kiến quý báu cho tác giả trình nghiên cứu Tác giả bày tỏ lòng biết ơn đến TS Trần Thương Bình anh chị em phịng thí nghiệm Địa kỹ thuật Trường Đại Học Kiến trúc Hà Nội, Phịng thí nghiệm Trường Đại học Xây dựng tạo điều kiện thuận lợi, giúp đỡ, động viên, trình thực nghiên cứu thực nghiệm luận án Tác giả xin bày tỏ lòng biết ơn chân thành đến đơn vị nơi tác giả công tác Trường Đại học Kiến trúc Hà Nội, tạo điều kiện thuận lợi, giúp đỡ, động viên để tác giả yên tâm tập trung nghiên cứu để hồn thành luận án Để hồn thành luận án tác giả động viên, ủng hộ, chia sẻ kịp thời từ gia đình lúc khó khăn nhất, tác giả xin bày tỏ lòng biết ơn chia sẻ thành cơng có thân đến gia đình Cuối cùng, tác giả xin chân thành cảm ơn đến bạn bè, đồng nghiệp giúp đỡ, động viên, ủng hộ, chia sẻ q trình tác giả hồn thành luận án iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC HÌNH VẼ VÀ ĐỒ THỊ xii DANH MỤC BẢNG BIỂU xi DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT vii MỞ ĐẦU 1 Tính cấp thiết đề tài Đối tượng phạm vi nghiên cứu Mục tiêu Nội dung nghiên cứu luận án Cách tiếp cận phương pháp nghiên cứu Những luận điểm bảo vệ Những điểm khoa học Ý nghĩa khoa học thực tiễn Cơ sở tài liệu CHƯƠNG TỔNG QUAN CÁC VẤN ĐỀ NGHIÊN CỨU ĐÁNH GIÁ ĐIỀU KIỆN ĐỊA KỸ THUẬT HỐ ĐÀO SÂU 1.1 THỰC TRẠNG NGHIÊN CỨU LIÊN QUAN ĐẾN ĐỊA KỸ THUẬT PHỤC VỤ THIẾT KẾ HỐ ĐÀO SÂU: 1.1.1 Đặc điểm phạm vi nghiên cứu hố đào sâu 1.1.2 Các tượng địa kỹ thuật xảy thi công hố đào sâu 11 1.1.3 Các yếu tố ảnh hưởng đến dịch chuyển đất xung quanh hố đào sâu 12 1.1.4 Trạng thái ứng suất đất [38] 13 1.1.5 Tính tốn, thiết kế thi công hố đào sâu 15 1.1.6 Đánh giá ưu điểm, nhược điểm số mơ hình đất 17 1.2 CÁC NGHIÊN CỨU LIÊN QUAN ĐẾN KHẢO SÁT CUNG CẤP SỐ LIỆU ĐỊA KỸ THUẬT PHỤC VỤ THIẾT KẾ THI CÔNG HỐ ĐÀO SÂU 22 iv 1.2.1 Các phương pháp thí nghiệm địa kỹ thuật 27 1.2.2 Các thông số đặc trưng đất toán hố đào sâu 30 1.3 KẾT LUẬN CHƯƠNG 1: 35 CHƯƠNG THÍ NGHIỆM XÁC ĐỊNH THÔNG SỐ ĐẶC TRƯNG CỦA ĐẤT NỀN TRÊN MÁY BA TRỤC GIẢM ỨNG SUẤT NGANG 36 2.1 MỞ ĐẦU 36 2.2 CÁC PHƯƠNG PHÁP THÍ NGHIỆM HIỆN NAY XÁC ĐỊNH THÔNG SỐ ĐẦU VÀO ĐỂ THIẾT KẾ THI CÔNG HỐ ĐÀO SÂU 37 2.2.1 Thông số địa kỹ thuật để thiết kế thi công hố đào sâu 37 2.2.2 Thí nghiệm nén ba trục: 38 2.3 MƠ HÌNH THÍ NGHIỆM BA TRỤC GIẢM ỨNG SUẤT NGANG 39 2.3.1 Nội dung phương pháp thí nghiệm 40 2.3.2 Thiết bị, dụng cụ thí nghiệm 41 2.3.3 Quy trình thí nghiệm 42 2.3.4 Dỡ tải 44 2.3.5 Tính tốn, báo cáo kết thí nghiệm 45 2.3.7 Xử lý số liệu thí nghiệm 46 2.3.8 Kết thí nghiệm 47 2.4 NHẬN XÉT CHƯƠNG 55 CHƯƠNG MƠ HÌNH LADE CẢI TIẾN CHO BÀI TOÁN HỐ ĐÀO SÂU 58 3.1 ỨNG XỬ CỦA ĐẤT NỀN 58 3.2 MƠ HÌNH ĐẤT NỀN LADE 58 3.2.1 Giới thiệu mơ hình Lade 58 3.2.2 Lý lựa chọn mơ hình Lade 58 3.2.3 Mơ hình Lade 59 3.3 MƠ HÌNH ĐẤT NỀN LADE CẢI TIẾN 65 3.3.1 Cơ sở để xây dựng mơ hình Lade cải tiến 65 3.3.2 Đặc trưng đàn hồi 66 3.3.3 Mặt phá hoại 68 v 3.3.4 Mặt biến dạng dẻo 70 3.3.5 Mặt chảy dẻo 74 3.3.6 Xác định đặc trưng mơ hình Lade cải tiến 79 3.3.7 Đánh giá độ tin cậy mơ hình Lade cải tiến 80 3.4 PHẦN MỀM TÍNH TỐN THEO MƠ HÌNH ĐẤT NỀN LỰA CHỌN 84 CHƯƠNG ÁP DỤNG TÍNH TỐN HỐ ĐÀO SÂU THEO MƠ HÌNH LADE CẢI TIẾN 89 4.1 NỘI DUNG TÍNH TỐN 89 4.2 MƠ HÌNH HĨA HỐ ĐÀO SÂU 90 4.2.1 Cơng trình Bệnh viện 108 91 4.2.2 Cơng trình Hoabinh Green City Minh Khai 98 4.2.3 Cơng trình đập thủy điện số .104 KẾT LUẬN VÀ KIẾN NGHỊ 112 KẾT LUẬN 112 KIẾN NGHỊ 113 CÁC CƠNG TRÌNH KHOA HỌC TÁC GIẢ ĐÃ CƠNG BỐ .114 DANH MỤC TÀI LIỆU THAM KHẢO 115 PHỤ LỤC 120 PHỤ LỤC A CHẾ TẠO THIẾT BỊ THÍ NGHIỆM 120 A.1 Buồng ba trục 120 A.2 Thiết bị dùng để tạo trì áp lực nước 121 A.3 Thiết bị tạo lực dọc trục 122 A.4 Thiết bị giảm áp lực buồng .123 A.5 Thiết bị đo biến biến dạng, đo áp suất 123 A.6 Thiết bị ghi chép số liệu 128 A.7 Các phụ kiện khác 129 A.8 Kết thí nghiệm 129 PHỤ LỤC B: PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN .143 B.1 Thành phần ứng suất đất 143 vi B.2 Thành phần biến dạng đất 144 B.3 Quan hệ ứng suất biến dạng đất 145 B.4 Công dẻo 145 B.5 Ứng xử tái bền 146 B.6 Ma trận đàn dẻo .148 B.8 Phần tử dàn 151 B.9 Phần tử dầm 152 B.10 Phần tử tam giác đẳng tham số .154 B.11 Phần tử tiếp xúc .157 B.12 Các dạng phân tích 158 B.13 Giải hệ phương trình 159 B.14 Phương pháp giải toán phi tuyến .160 B.15 Mô tả q trình thi cơng 162 B.17 Mô đường cong ứng suất - biến dạng thí nghiệm nén ba trục theo mơ hình Lade cải tiến 165 PHỤ LỤC C: MÃ NGUỒN PHẦN MỀM LADEDEEP .167 vii DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT A, B Thông số áp lực nước lỗ rỗng A, Rf, Kl, n, G, F Thông số mơ hình A, b Hệ số phương trình tuyến tính, dùng để xác định giá trị hệ số Poisson, xác định theo kết thí nghiệm Α Hệ số xác định theo tính chất vật liệu a Khoảng cách từ mặt đất đến điểm đặt Ea [A] Ma trận biến trường chuyển vị [A]e Ma trận tọa độ cụ thể nút Bt Mơ đun biến dạng thể tích tiếp tuyến [B] Ma trận biến đổi chuyển vị nút biến dạng tỷ đối CD Thí nghiệm nén ba trục cố kết nước CU Thí nghiệm nén ba trục cố kết khơng nước c Lực dính đất [Dp] Ma trận dẻo [D] Ma trận đàn hồi đất tốn biến dạng phẳng E Mơđun đàn hồi đất E0 Mô đun biến dạng đất Ea Tổng áp lực đất chủ động xác định từ hiệu áp lực đất bị động áp lực đất chủ động tồn chiều sâu trường Ei Mơđun đàn hồi ban đầu Et Môđun đàn hồi tiếp tuyến E20, E25 Môđun đàn hồi cát tuyến e Hệ số rỗng đất Fx, Fy Thành phần thể tích tác dụng theo phương x, y đơn vị thể tích áp lực [F]e Lực nút tương đương phần tử hữu hạn (PTHH) [F]s Véctơ lực nút tương đương kết cấu rời rạc G Mô đun đàn hồi trượt viii G rsief Độ cứng kháng cắt chuẩn ban đầu HĐS Hố đào sâu K Hệ số theo phương ngang K0 Hệ số áp lực đất tĩnh Ka Hệ số áp lực đất chủ động Kp Hệ số áp lực đất bị động K0,NC Hệ số áp lực đất trạng thái cố kết bình thường K0,OC Hệ số áp lực đất tĩnh đất cố kết Kb,m Thơng số mơ hình [K]e Ma trận độ cứng PTHH [K]s Ma trận độ cứng tổng thể kết cấu ksi Mô đun phản lực độ sâu xét ks Mô đun phản lực l Khoảng cách từ mặt đất đến điểm có áp lực zero m Số mũ độ cứng N Giá trị thể tích riêng đất cố kết bình thường [N] Ma trận đạo hàm riêng trường chuyển vị PTHH OCR Hệ số cố kết đất P Tải trọng gây chuyển vị δ p Áp lực đất thực tác dụng lên tường po Áp lực đất tĩnh lên tường pa Áp lực khí σc Ứng suất tiền cố kết đất Pref Sức bền chuẩn Px, Py Các thành phần ngoại lực theo phương x, y tác dụng đơn vị diện tích mặt ngòai vật thể đàn hồi Rf Độ huy động phá hoại SK Độ cố kết t Độ dày phân tử ix t0 Độ sâu ngàm cần thiết UU Thí nghiệm nén ba trục khơng cố kết khơng nước [U]s Véctơ chuyển vị nút kết cấu rời rạc hóa ux Chuyển vị điểm theo phương x uy Chuyển vị điểm theo phương y [u] Véctơ chuyển vị thành phần điểm PTHH [u]s Véctơ thành phần chuyển vị nút PTHH v Thể tích đặc trưng v0 Thể tích đặc trưng ban đầu đất ΔU0 Số gia mật độ lượng biến dạng Δσ Độ lệch ứng suất ΔW Số gia công học Δs Chuyển vị [α] Ma trận thông số αi hàm chuyển vị δ Dịch chuyển ngang tường δW Công học δu Chuyển vị hệ [δε] Véctơ biến dạng [δu] Véctơ chuyển vị εl Biến dạng đứng ( biến dạng dọc trục) εe Biến dạng đàn hồi εp Biến dạng dẻo εv Biến dạng thể tích [ε] Véctơ biến dạng tỷ đối γn Trọng lượng riêng đất lớp thứ n γw Trọng lượng riêng nước φ Góc ma sát đất, góc ma sát tiếp xúc φn Góc ma sát đất lớp thứ n 204 // If i1o3>so.somal^[so.plstg^[(i-1)*so.nstg+h]].pa then si:=so.somal^[so.plstg^[(i-1)*so.nstg+h]].si-so.somal^[so.plstg^[(i1)*so.nstg+h]].dsi*ln(i1o3/so.somal^[so.plstg^[(i-1)*so.nstg+h]].pa)/ln(10) Else si:=so.somal^[so.plstg^[(i-1)*so.nstg+h]].si; If si>so.somal^[so.plstg^[(i-1)*so.nstg+h]].simax then si:=so.somal^[so.plstg^[(i-1)*so.nstg+h]].simax; If siso.somal^[so.plstg^[(i-1)*so.nstg+h]].fimax then fi:=so.somal^[so.plstg^[(i-1)*so.nstg+h]].fimax; If fiso.somal^[so.plstg^[(i-1)*so.nstg+h]].pa then si:=so.somal^[so.plstg^[(i-1)*so.nstg+h]].si-so.somal^[so.plstg^[(i1)*so.nstg+h]].dsi*ln(i1o3/so.somal^[so.plstg^[(i-1)*so.nstg+h]].pa)/ln(10) Else si:=so.somal^[so.plstg^[(i-1)*so.nstg+h]].si; If si>so.somal^[so.plstg^[(i-1)*so.nstg+h]].simax then si:=so.somal^[so.plstg^[(i-1)*so.nstg+h]].simax; If si0 then Begin k:=0; For j:=1 to so.naf If so.af^[j].sc=h then If so.af^[j].num=i then k:=j; // If k=0 then Begin For j:=1 to Begin coor2[j,1]:=so.jnt^[so.anc^[i].jn[j]].x; coor2[j,2]:=so.jnt^[so.anc^[i].jn[j]].y; End; cosin2(cs2,coor2); //D For j:=1 to For k:=1 to glodis2[(j-1)*2+k]:=ddis^[(so.anc^[i].jn[j]-1)*3+k]; //D'=CS*D For j:=1 to Begin locdis2[j]:=0; For k:=1 to locdis2[j]:=locdis2[j]+cs2[j,k]*glodis2[k]; End; // F'=K*D' ell:=sqrt(sqr(coor2[2,1]-coor2[1,1])+sqr(coor2[2,2]-coor2[1,2])); stiffmatanchor(kmi2,so.ancmal^[so.ancstg^[(i-1)*so.nstg+h]].ae,ell,so.ancmal^[so.ancstg^[(i1)*so.nstg+h]].lsp); For j:=1 to Begin locload2[j]:=0; For k:=1 to locload2[j]:=locload2[j]+kmi2[j,k]*locdis2[k]; End; 212 so.ancifs^[(h-1)*so.nanc+i]:=so.ancifs^[(h-1)*so.nanc+i]+locload2[3]; //No compression If so.ancifs^[(h-1)*so.nanc+i]so.ancmal^[so.ancstg^[(i-1)*so.nstg+h]].np then so.ancifs^[(h-1)*so.nanc+i]:=so.ancmal^[so.ancstg^[(i-1)*so.nstg+h]].np; // locload2[1]:=-so.ancifs^[(h-1)*so.nanc+i]; locload2[2]:=0; locload2[3]:=so.ancifs^[(h-1)*so.nanc+i]; locload2[4]:=0; //F=CST*F' For j:=1 to Begin gloload2[j]:=0; For k:=1 to gloload2[j]:=gloload2[j]+cs2[k,j]*locload2[k]; End; //Global load vector For j:=1 to Begin For k:=1 to If so.cdis^[(so.anc^[i].jn[j]-1)*3+k]=0 then dload^[(so.anc^[i].jn[j]-1)*3+k]:=dload^[(so.anc^[i].jn[j]-1)*3+k]+gloload2[(j-1)*2+k]; End; End// Else Begin so.ancifs^[(h-1)*so.nanc+i]:=(step/so.ictrl.nstep)*so.af^[k].f; //Maximum force If so.ancifs^[(h-1)*so.nanc+i]>so.ancmal^[so.ancstg^[(i-1)*so.nstg+h]].np then so.ancifs^[(h-1)*so.nanc+i]:=so.ancmal^[so.ancstg^[(i-1)*so.nstg+h]].np; End; End; // //Beam-Column For i:=1 to so.nbc If so.bcstg^[(i-1)*so.nstg+h]>0 then Begin For j:=1 to Begin coor2[j,1]:=so.jnt^[so.bc^[i].jn[j]].x; coor2[j,2]:=so.jnt^[so.bc^[i].jn[j]].y; End; cosin3(cs3,coor2); //D For j:=1 to For k:=1 to glodis3[(j-1)*3+k]:=ddis^[(so.bc^[i].jn[j]-1)*3+k]; //D'=CS*D 213 For j:=1 to Begin locdis3[j]:=0; For k:=1 to locdis3[j]:=locdis3[j]+cs3[j,k]*glodis3[k]; End; // F'=K*D' ell:=sqrt(sqr(coor2[2,1]-coor2[1,1])+sqr(coor2[2,2]-coor2[1,2])); axx:=ax(so.spro^[so.bcstg^[(i-1)*so.nstg+h]].code,so.spro^[so.bcstg^[(i1)*so.nstg+h]].yd,so.spro^[so.bcstg^[(i-1)*so.nstg+h]].zd,so.spro^[so.bcstg^[(i1)*so.nstg+h]].yb,so.spro^[so.bcstg^[(i-1)*so.nstg+h]].zb,so.spro^[so.bcstg^[(i-1)*so.nstg+h]].d); izz:=iz(so.spro^[so.bcstg^[(i-1)*so.nstg+h]].code,so.spro^[so.bcstg^[(i1)*so.nstg+h]].yd,so.spro^[so.bcstg^[(i-1)*so.nstg+h]].zd,so.spro^[so.bcstg^[(i1)*so.nstg+h]].yb,so.spro^[so.bcstg^[(i-1)*so.nstg+h]].zb,so.spro^[so.bcstg^[(i-1)*so.nstg+h]].d); stiffmatbeamcolumn(kmi3,axx,izz,so.smal^[so.spro^[so.bcstg^[(i1)*so.nstg+h]].nm].e,so.smal^[so.spro^[so.bcstg^[(i-1)*so.nstg+h]].nm].p,ell); For j:=1 to Begin locload3[j]:=0; For k:=1 to locload3[j]:=locload3[j]+kmi3[j,k]*locdis3[k]; End; // For j:=1 to so.bcifs^[(h-1)*so.nbc+i].vl[j]:=so.bcifs^[(h-1)*so.nbc+i].vl[j]+locload3[j]; // //Place procedure for nonlinear beam here //so.BeamNonlinearForce(i,h,fmax,fmin); // For j:=1 to locload3[j]:=so.bcifs^[(h-1)*so.nbc+i].vl[j]; //F=CST*F' For j:=1 to Begin gloload3[j]:=0; For k:=1 to gloload3[j]:=gloload3[j]+cs3[k,j]*locload3[k]; End; //Global load vector For j:=1 to Begin For k:=1 to If so.cdis^[(so.bc^[i].jn[j]-1)*3+k]=0 then dload^[(so.bc^[i].jn[j]-1)*3+k]:=dload^[(so.bc^[i].jn[j]-1)*3+k]+gloload3[(j-1)*3+k]; End; End; //Bar For i:=1 to so.nbr If so.brstg^[(i-1)*so.nstg+h]>0 then Begin For j:=1 to 214 Begin coor2[j,1]:=so.jnt^[so.br^[i].jn[j]].x; coor2[j,2]:=so.jnt^[so.br^[i].jn[j]].y; End; cosin2(cs2,coor2); //D For j:=1 to For k:=1 to glodis2[(j-1)*2+k]:=ddis^[(so.br^[i].jn[j]-1)*3+k]; //D'=CS*D For j:=1 to Begin locdis2[j]:=0; For k:=1 to locdis2[j]:=locdis2[j]+cs2[j,k]*glodis2[k]; End; // F'=K*D' ell:=sqrt(sqr(coor2[2,1]-coor2[1,1])+sqr(coor2[2,2]-coor2[1,2])); axx:=ax(so.spro^[so.brstg^[(i-1)*so.nstg+h]].code,so.spro^[so.brstg^[(i1)*so.nstg+h]].yd,so.spro^[so.brstg^[(i-1)*so.nstg+h]].zd,so.spro^[so.brstg^[(i1)*so.nstg+h]].yb,so.spro^[so.brstg^[(i-1)*so.nstg+h]].zb,so.spro^[so.brstg^[(i-1)*so.nstg+h]].d); stiffmatbar(kmi2,so.smal^[so.spro^[so.brstg^[(i-1)*so.nstg+h]].nm].e,axx,ell,1); For j:=1 to Begin locload2[j]:=0; For k:=1 to locload2[j]:=locload2[j]+kmi2[j,k]*locdis2[k]; End; so.brifs^[(h-1)*so.nbr+i]:=so.brifs^[(h-1)*so.nbr+i]+locload2[3]; //Minimum force If so.brifs^[(h-1)*so.nbr+i](1E+15) then so.brifs^[(h-1)*so.nbr+i]:=(1E+15); // locload2[1]:=-so.brifs^[(h-1)*so.nbr+i]; locload2[2]:=0; locload2[3]:=so.brifs^[(h-1)*so.nbr+i]; locload2[4]:=0; //F=CST*F' For j:=1 to Begin gloload2[j]:=0; For k:=1 to gloload2[j]:=gloload2[j]+cs2[k,j]*locload2[k]; End; //Global load vector For j:=1 to Begin For k:=1 to 215 If so.cdis^[(so.br^[i].jn[j]-1)*3+k]=0 then dload^[(so.br^[i].jn[j]-1)*3+k]:=dload^[(so.br^[i].jn[j]-1)*3+k]+gloload2[(j-1)*2+k]; End; End; //Spring For i:=1 to so.nspr If so.sprstg^[(i-1)*so.nstg+h]>0 then Begin For j:=1 to Begin coor2[j,1]:=so.jnt^[so.spr^[i].jn[j]].x; coor2[j,2]:=so.jnt^[so.spr^[i].jn[j]].y; End; cosin3(cs3,coor2); //D For j:=1 to For k:=1 to glodis3[(j-1)*3+k]:=ddis^[(so.spr^[i].jn[j]-1)*3+k]; //D'=CS*D For j:=1 to Begin locdis3[j]:=0; For k:=1 to locdis3[j]:=locdis3[j]+cs3[j,k]*glodis3[k]; End; // F'=K*D' stiffmatspring(kmi3,so.sprmal^[so.sprstg^[(i-1)*so.nstg+h]].k,so.sprmal^[so.sprstg^[(i1)*so.nstg+h]].lsp); For j:=1 to Begin locload3[j]:=0; For k:=1 to locload3[j]:=locload3[j]+kmi3[j,k]*locdis3[k]; End; // For j:=1 to so.sprifs^[(h-1)*so.nspr+i].vl[j]:=so.sprifs^[(h-1)*so.nspr+i].vl[j]+locload3[j]; // so.SpringNonlinearForce(i,h); // For j:=1 to locload3[j]:=so.sprifs^[(h-1)*so.nspr+i].vl[j]; //F=CST*F' For j:=1 to Begin gloload3[j]:=0; For k:=1 to gloload3[j]:=gloload3[j]+cs3[k,j]*locload3[k]; End; //Global load vector For j:=1 to 216 Begin For k:=1 to If so.cdis^[(so.spr^[i].jn[j]-1)*3+k]=0 then dload^[(so.spr^[i].jn[j]-1)*3+k]:=dload^[(so.spr^[i].jn[j]-1)*3+k]+gloload3[(j-1)*3+k]; End; End; //Tinh so gia chuyen vi k:=0; For i:=1 to so.njnt For j:=1 to If so.cdis^[(i-1)*3+j]=0 then Begin k:=k+1; so.pm^[k]:=preload^[(h-1)*3*so.njnt+(i-1)*3+j]+(step/so.ictrl.nstep)*(load^[(h-1)*3*so.njnt+(i1)*3+j]-preload^[(h-1)*3*so.njnt+(i-1)*3+j])-dload^[(i-1)*3+j]; End; // If so.nblock=1 then Begin so.chobk1_1(1); so.chobk2_1(1); End; If so.nblock>1 then Begin so.chobk1(1); so.chobk2(1); End; // k:=0; For i:=1 to so.njnt For j:=1 to Begin If so.cdis^[(i-1)*3+j]=0 then Begin k:=k+1; totd^[(i-1)*3+j]:=totd^[(i-1)*3+j]+so.pm^[k]; ddis^[(i-1)*3+j]:=so.pm^[k]; End Else ddis^[(i-1)*3+j]:=0; End; //Kiem tra hoi tu told:=0; For i:=1 to so.njnt told:=told+sqrt(sqr(ddis^[(i-1)*3+1])+sqr(ddis^[(i-1)*3+2])); // ttold:=0; For i:=1 to so.njnt ttold:=ttold+sqrt(sqr(totd^[(i-1)*3+1])+sqr(totd^[(i-1)*3+2])); // told:=told/ttold; 217 // Until (toldso.ictrl.toldis); // //Keep last stage balance load For i:=1 to 3*so.njnt preload^[(h-1)*3*so.njnt+i]:=dload^[i]; End;//h>1 // For i:=1 to so.njnt For j:=1 to so.dis^[(h-1)*so.njnt*3+(so.renj^[i]-1)*3+j]:=totd^[(i-1)*3+j]; // Freemem(so.pm,so.neq*sizeof(real)); // frmssi2d.runmode:=1; End Else Begin ListBox.Items.Add('Error: Please check the boundary conditions !');num:=num+1;ListBox.ItemIndex:=num; End; // If so.nblock=1 then Freemem(so.km,so.nkm*sizeof(real)); so.nkm:=0; Freemem(so.block,so.nblock*sizeof(BlockInf)); so.nblock:=0; Freemem(so.cline,so.neq*sizeof(cblock)); so.neq:=0; // Until (h=so.nstg) or (err1=1) or (iter=so.ictrl.niter); // so.DelJCon; // For i:=1 to so.nblock Begin If FileExists(so.MatrixDirectory+'\km'+IntToStr(i)+'.mat') then DeleteFile(so.MatrixDirectory+'\km'+IntToStr(i)+'.mat'); End; //Xoa thu muc ma tran If DirectoryExists(so.MatrixDirectory) then If DeleteDir(so.MatrixDirectory)=true then Begin //Do nothing End; // Freemem(ddis,3*so.njnt*sizeof(real)); // Freemem(dload,3*so.njnt*sizeof(real)); 218 Freemem(totd,3*so.njnt*sizeof(real)); // Freemem(load,so.njnt*sizeof(real)); Freemem(preload,3*so.njnt*sizeof(real)); Freemem(strload,so.njnt*sizeof(real)); Freemem(slload,so.njnt*sizeof(real)); // Freemem(so.cdis,so.njnt*3*sizeof(integer)); Freemem(so.cdel,so.njnt*sizeof(integer)); // If err1=1 then Begin so.delresult; End; // so.rerenamejnt; //Change Sign For i:=1 to so.nstg*so.nbc Begin so.bcifs^[i].vl[1]:=-so.bcifs^[i].vl[1]; so.bcifs^[i].vl[3]:=-so.bcifs^[i].vl[3]; so.bcifs^[i].vl[5]:=-so.bcifs^[i].vl[5]; End; // ListBox.Items.Add('Finish');num:=num+1;ListBox.ItemIndex:=num; CloseFile(fxt); End; ... VẤN ĐỀ NGHIÊN CỨU ĐÁNH GIÁ ĐIỀU KIỆN ĐỊA KỸ THUẬT HỐ ĐÀO SÂU 1.1 THỰC TRẠNG NGHIÊN CỨU LIÊN QUAN ĐẾN ĐỊA KỸ THUẬT PHỤC VỤ THI? ??T KẾ HỐ ĐÀO SÂU: 1.1.1 Đặc điểm phạm vi nghiên cứu hố đào sâu Trong... NGHIÊN CỨU ĐÁNH GIÁ ĐIỀU KIỆN ĐỊA KỸ THUẬT HỐ ĐÀO SÂU 1.1 THỰC TRẠNG NGHIÊN CỨU LIÊN QUAN ĐẾN ĐỊA KỸ THUẬT PHỤC VỤ THI? ??T KẾ HỐ ĐÀO SÂU: 1.1.1 Đặc điểm phạm vi nghiên cứu hố đào. .. SÁT CUNG CẤP SỐ LIỆU ĐỊA KỸ THUẬT PHỤC VỤ THI? ??T KẾ THI CÔNG HỐ ĐÀO SÂU Để thi? ??t kế kết cấu chắn giữ HĐS cần có số liệu cung cấp từ khảo sát địa kỹ thuật, sau lựa chọn kiểu kết cấu chắn giữ tính