Header Page of 148 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 Footer Page of 148 Header Page of 148 i LỜI CAM ĐOAN Tôi xin cam đoan công trình khoa học 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ả hoàn toà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 Footer Page of 148 Header Page of 148 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 trình thực hoà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 để hoàn thành luận án Để hoà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ẻ trình tác giả hoàn thành luận án Footer Page of 148 Header Page of 148 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 toá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 Footer Page of 148 Header Page of 148 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 toá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 Footer Page of 148 Header Page of 148 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 TOÁN THEO MÔ HÌNH ĐẤT NỀN LỰA CHỌN 84 CHƯƠNG ÁP DỤNG TÍNH TOÁN HỐ ĐÀO SÂU THEO MÔ HÌNH LADE CẢI TIẾN 89 4.1 NỘI DUNG TÍNH TOÁ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 Footer Page of 148 Header Page of 148 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ả 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 Footer Page of 148 Header Page of 148 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 thoát nước CU Thí nghiệm nén ba trục cố kết không thoát nước c Lực dính đất [Dp] Ma trận dẻo [D] Ma trận đàn hồi đất toá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 toà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 Footer Page of 148 Header Page of 148 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ử Footer Page of 148 Header Page 10 of 148 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 thoát 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 Footer Page 10 of 148 Header Page 219 of 148 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; Footer Page 226 of 148 Header Page 227 of 148 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 Footer Page 227 of 148 Header Page 228 of 148 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 Footer Page 228 of 148 Header Page 229 of 148 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 Footer Page 229 of 148 Header Page 230 of 148 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 Footer Page 230 of 148 Header Page 231 of 148 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; Footer Page 231 of 148 Header Page 232 of 148 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)); Footer Page 232 of 148 Header Page 233 of 148 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; Footer Page 233 of 148 ... hng n s lm vic ca h o h) nh hng s dng cỏc phng phỏp thi cụng, quy trỡnh thi cụng, trỡnh t thi cụng Trong thit k thi cụng HS cú nhiu ni dung nh thit k h tng chn, chng y tri, h mc nc ngm, h chng... bi ỏp lc hụng gim dn vi cựng ln, tc l 1= - 1.1.5 Tớnh toỏn, thit k thi cụng h o sõu Mc ớch ca cụng tỏc thit k thi cụng HS l thit k thi cụng h kt cu chn gi thnh h o, m bo cỏc yờu cu sau [10]:... TO THIT B TH NGHIM 120 A.1 Bung ba trc 120 A.2 Thit b dựng to v trỡ ỏp lc nc 121 A.3 Thit b to lc dc trc 122 A.4 Thit b gim ỏp lc bung .123 A.5 Thit