Ứng xử của đất nền dưới đáy móng công trình biển trọng lực chịu tác động của tổ hợp các dạng tải trọng (luận án (theses))

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẶNG XUÂN TRƯỜNG ỨNG XỬ CỦA ĐẤT NỀN DƯỚI ĐÁY MĨNG CƠNG TRÌNH BIỂN TRỌNG LỰC CHỊU TÁC ĐỘNG CỦA TỔ HỢP CÁC DẠNG TẢI TRỌNG LUẬN ÁN TIẾN SĨ KĨ THUẬT TP HỒ CHÍ MINH NĂM 2019 ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẶNG XUÂN TRƯỜNG LUẬN ÁN TIẾN SĨ KĨ THUẬT ỨNG XỬ CỦA ĐẤT NỀN DƯỚI ĐÁY MĨNG CƠNG TRÌNH BIỂN TRỌNG LỰC CHỊU TÁC ĐỘNG CỦA TỔ HỢP CÁC DẠNG TẢI TRỌNG Chuyên ngành: KỸ THUẬT ĐỊA CHẤT Mã số chuyên ngành: 62.52.05.01 Phản biện độc lập 1: PGS.TS Vũ Chí Hiếu Phản biện độc lập 2: PGS.TS Tạ Đức Thịnh Phản biện 1: PGS.TS Lê Văn Nam Phản biện 2: PGS.TS Đoàn Thế Tường Phản biện 3: PGS.TS Nguyễn Minh Tâm NGƯỜI HƯỚNG DẪN KHOA HỌC PGS.TS Đậu Văn Ngọ PGS.TS Nguyễn Hữu Bảng LỜI CAM ĐOAN Tác giả xin cam đoan cơng trình nghiên cứu thân tác giả Các kết nghiên cứu kết luận luận án trung thực, không chép từ nguồn hình thức nào.Việc tham khảo nguồn tài liệu thực trích dẫn ghi nguồn tài liệu tham khảo quy định Tác giả luận án Chữ ký Đặng Xn Trường TĨM TẮT LUẬN ÁN Cơng trình biển trọng lực hai dạng chủ yếu áp dụng xây dựng cơng trình biển cố định, bao gồm khơi ven bờ Tải trọng tác động lên cơng trình vấn đề móng cịn nhiều thách thức nhà thiết kế Hơn nữa, chi phí xây dựng dạng cơng trình thường lớn Chính vậy, việc xác định sở liệu đầu vào phương pháp tính tốn nhằm đảm bảo tính ổn định lâu dài mà tiết kiệm chi phí đầu tư vấn đề thực phải quan tâm Luận án thực với mong muốn giải vài vấn đề số nhiều vấn đề cần phải quan tâm nêu Các nội dung nghiên cứu cụ thể sau:  Phần thứ luận án xác định dạng tải trọng tác động vào cơng trình biển trọng lực Trong lý thuyết hồi quy quy hoạch thực nghiệm áp dụng để xác định thông số ngẫu nhiên tải trọng mà cụ thể thơng số sóng biển Trong phần này, tác giả xây dựng tốn xác định áp lực sóng cực đại tác động lên trụ cơng trình biển trọng lực kích thước lớn có tiết diện thay đổi đặt vùng nước sâu  Phần thứ hai toán xác định biến dạng đất phân tích theo học vật rắn biến dạng lý thuyết đàn hồi bán không gian vô hạn Xây dựng phần mềm máy tính để thiết lập miền đồng mức chuyển vị đất đáy móng cơng trình biển trọng lực chịu tác động đồng thời dạng tải trọng gây Các nghiên cứu luận án tác giả xây dựng phần mềm tính tốn nhằm tự động hóa q trình tính tốn cần thay đổi thông số đầu vào Các kết nghiên cứu luận án công bố tạp chí uy tín nước nước (02 báo Quốc tế, 01 báo nước báo khác có liên quan) ABSTRACT Gravity based structure (GBS) is one of two main types used in the construction of fixed marine structures, including offshore and inshore The impact of load on the structure and foundation issues are still challenging for designers Moreover, the cost of construction on these types of structures is usually very large Therefore, the determination of the input database and the method of calculation to ensure long-term stability and save the cost of investment are the real issues to pay attention This dissertation is emplemented with the desire to solve the several problems of the many issues that need attention as mentioned above Specific research contents are as follows:  The part 1: The thesis determines the types of load to impact on gravity based structures In which the regression theory and experimental planning are applied to determine the random parameters of the load, namely the wave parameters In this Chapter, the author has prepared the calculation to determine the maximum wave pressure acting on a large gravity based structures with varied cross section located in deep water  The part 2: This is the main content of the thesis It said that the vertical load impacted on the bottom of the foundation to transform from the uniform load to cross section on the entire bottom of the foundation The problem of determining the deformation of the soil is analyzed by solid and elastic theory in the semiinfinite space Computer program establish in the contour domain of the soil substrate under the bottom of the gravity based structures subjected to simultaneous effects of various types of load The main content of thesis is prepared by the author aim to automate the calculation process when the input parameters are changed The research results in the dissertation are also published in prestigious domestic and international journal (two international and one national journal) LỜI CẢM ƠN Lời đầu tiên, xin gửi lời cảm ơn chân thành tới Ban giám hiệu, Phòng đào tạo Sau đại học, Khoa Kỹ thuật Địa chất Dầu khí, Bộ mơn Địa kỹ thuật trường đại học Bách khoa – ĐHQG TP Hồ Chí Minh tạo điều kiện giúp đỡ tơi q trình học tập nghiên cứu trường Tôi xin bày tỏ lòng biết ơn sâu sắc đến NGƯT.PGS.TS Đậu Văn Ngọ PGS.TS Nguyễn Hữu Bảng (Trường Đại học Mỏ Địa chất) tận tình giúp đỡ, hướng dẫn động viên thực luận án Tôi xin gửi lời cảm ơn chân thành đến TS Tạ Quốc Dũng (Trưởng khoa Kỹ thuật Kỹ thuật Địa chất Dầu khí), PGS.TS Nguyễn Việt Kỳ (Chủ nhiệm Bộ môn Địa Kỹ thuật), PGS.TS Trần Văn Xuân (Chủ nhiệm Bộ mơn Địa chất – Dầu khí) Q thầy khoa Kỹ thuật Địa chất Dầu khí nói chung, mơn Địa Kỹ thuật nói riêng tạo điều kiện, giúp đỡ động viên tơi hồn thành luận án Xin chân thành cảm ơn Đảng uỷ - Ban Giám hiệu khoa Kỹ thuật xây dựng trường đại học Giao thông vận tải TP Hồ Chí Minh (nơi tơi cơng tác) tạo điều kiện thuận lợi cho học tập nghiên cứu Cuối cùng, tơi xin bày tỏ lịng biết ơn sâu sắc đến tất người thân gia đình động viên khích lệ tơi hồn thành nhiệm vụ Thành phố Hồ Chí Minh, năm 2019 Tác giả luận án Đặng Xuân Trường MỤC LỤC DANH MỤC CÁC HÌNH ẢNH DANH MỤC BẢNG BIỂU 12 DANH MỤC CÁC TỪ VIẾT TẮT .13 MỞ ĐẦU .14 Tính cấp thiết đề tài nghiên cứu 14 Đối tượng phạm vi nghiên cứu 14 Mục đích luận án .15 Những luận điểm bảo vệ 16 Nhiệm vụ luận án 16 Nội dung nghiên cứu luận án 17 Phương pháp nghiên cứu luận án 17 Những điểm mặt khoa học luận án .18 Ý nghĩa khoa học thực tiễn luận án 18 10 Cơ sở tài liệu luận án 19 11 Cấu trúc luận án 19 CHƯƠNG NGHIÊN CỨU TỔNG QUAN VỀ CƠNG TRÌNH BIỂN TRỌNG LỰC VÀ ĐẶC ĐIỂM MÔI TRƯỜNG BIỂN KHU VỰC NGHIÊN CỨU 20 1.1 Tổng quan cơng trình biển trọng lực 20 1.1.1 Cơng trình biển trọng lực giới 20 1.1.2 Phương pháp xây dựng 28 1.1.3 Cơng trình biển trọng lực Việt Nam 30 1.1.4 Nhận xét đề xuất 32 1.2 Đặc điểm khu vực nghiên cứu 33 1.2.1 Giới thiệu Bể Nam Côn Sơn .33 1.2.2 Các yếu tố cấu trúc kiến tạo 33 1.2.3 Lịch sử phát triển địa chất 34 1.2.4 Địa tầng, trầm tích môi trường .36 1.3 Đặc điểm môi trường dạng tải trọng tác động 44 1.3.1 Các tình để xác định tải trọng tác động lên CTBTL 44 1.3.2 Các loại tải trọng tác động lên CTBTL .44 1.3.3 Tổ hợp tải trọng 46 1.3.4 Xác định tải trọng tác động lên CTBTL 47 1.4 Tổng quan ứng xử đất đáy móng CTBTL .51 1.4.1 Tình hình nghiên cứu ứng xử đất đáy móng cơng trình biển trọng lực giới 51 1.4.2 Tình hình nghiên cứu ứng xử đất đáy móng cơng trình biển trọng lực Việt Nam .53 1.5 Nhận xét cuối chương 53 CHƯƠNG PHÂN TÍCH THƠNG SỐ SĨNG BẰNG LÝ THUYẾT HỒI QUY VÀ QUY HOẠCH THỰC NGHIỆM 55 2.1 Tải trọng sóng lý thuyết áp dụng 55 2.1.1 Lý thuyết sóng điều hịa 55 2.1.2 Lý thuyết sóng ngẫu nhiên 56 2.2 Đề xuất cách xác định thơng số sóng biển .57 2.2.1 Cách xác định thơng số sóng biển 57 2.2.2 Phương pháp xác định tải trọng sóng cực đại tác động lên trụ CTBTL có tiết diện thay đổi vùng nước sâu 58 2.3 Lý thuyết hồi quy .58 2.4 Tính hệ số hàm hồi quy .62 2.5 Quy hoạch thực nghiệm 64 2.6 Ý nghĩa hàm hồi quy 67 2.7 Áp dụng phương pháp hàm hồi quy để tính chiều cao sóng biển 68 2.7.1 Tổng hợp số liệu khảo sát 68 2.7.2 Kết tính tốn .69 2.8 Tính tốn tải trọng sóng tác dụng lên trụ cơng trình biển trọng lực có tiết diện thay đổi vùng nước sâu 70 2.8.1 Đặc điểm tải trọng sóng tác dụng lên trụ CTBTL 70 2.8.2 Tổng quan tình hình nghiên cứu tải trọng sóng tác dụng lên trụ 71 2.8.3 Nhận xét đề xuất 72 2.8.4 Cơ sở lý thuyết tốn sóng tác động lên trụ 73 2.8.5 Áp dụng cho mơ hình nghiên cứu .74 CHƯƠNG ÁP DỤNG LÝ THUYẾT ĐÀN HỒI VÀ CƠ HỌC VẬT RẮN BIẾN DẠNG ĐỂ XÁC ĐỊNH CHUYỂN VỊ CỦA NỀN ĐẤT DƯỚI ĐÁY MĨNG CƠNG TRÌNH BIỂN TRỌNG LỰC .79 3.1 Các vấn đề đất 79 3.2 Móng cơng trình biển trọng lực 79 3.3 Sức chịu tải móng chịu tải trọng nén tâm 80 3.4 Móng trọng lực chịu tải trọng tổng quát 82 3.5 Chuyển vị đất bán không gian vô hạn .83 3.5.1 Chuyển vị bán không gian vô hạn tải trọng tập trung gây 83 3.5.2 Thuật tốn để lập chương trình máy tính 87 3.5.3 Tính chuyển vị theo phương r với tải trọng phân bố q - tâm 88 3.5.4 Tính chuyển vị theo phương r với trường hợp tổng quát chịu tải trọng phân bố q momen M 90 3.5.5 Đề xuất áp dụng tính chuyển vị đất 91 CHƯƠNG NGHIÊN CỨU ỨNG XỬ CỦA ĐẤT NỀN DƯỚI ĐÁY MĨNG CƠNG TRÌNH BIỂN TRỌNG LỰC TẠI KHU VỰC MỎ ĐẠI HÙNG .92 4.1 Điều kiện địa chất cơng trình 92 4.1.1 Đánh giá điều kiện địa chất khu vực mỏ Đại Hùng 92 4.1.2 Xác định tải trọng tác động .97 4.1.3 Kết tính tốn chương trình máy tính 98 4.1.4 Nhận xét 100 4.2 Bài toán so sánh – Tính tốn sức chịu tải theo Cơ học đất .100 4.2.1 Áp lực đáy móng 100 4.2.2 Sức chịu tải 101 4.2.3 Tính lún cho cơng trình 101 4.2.4 Nhận xét 107 4.3 Mô hình tốn PLAXIS 108 4.3.1 Các thông số vật liệu nhập vào toán 108 4.3.2 Mơ hình tốn .108 4.3.3 Kết tính tốn .110 4.3.4 Phân tích ổn định so sánh độ lún vùng khác khối móng 111 4.4 Nhận xét kết 112 KẾT LUẬN VÀ KIẾN NGHỊ .113  Kết luận 113  Kiến nghị 114 DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ 115 TÀI LIỆU THAM KHẢO 117 PHỤ LỤC 120 PHỤ LỤC 125 PHỤ LỤC 129 PHỤ LỤC CHƯƠNG TRÌNH TÍNH HÀM HỒI QUY Private Sub Form_Load() Show Rem CHUONG TRINH TIM HAM HOI QUY Static a(4, 4), ad(4, 4), tich(10, 10), p(10), p1(4, 4) As Double Static km(4, 4), b(10), xo(10), y(10), dao(4, 4), xt(10, 10) As Double Static nhan1(4, 8), nhan2(10, 10), hesob(10), x(8, 4), xxdd(1000), xxra(1000) As Double Show Scale (-320, 240)-(320, -240) Rem BAI TOAN CUA TRUONG, DUA TRUC TIEP SO LIEU VAO y(i) y(1) = 7.07: y(2) = 6.7: y(3) = 4.87: y(4) = 5.7: y(5) = 6.27: y(6) = 6.33: y(7) = 7.2: y(8) = 6.7 Rem dua ve dang da chuan hoa thi nghiem Rem So lieu vao ma tran X dong = 8: cot = x(1, 2) = 10.7: x(1, 3) = 8.9: x(1, 4) = 1# x(2, 2) = 11.3: x(2, 3) = 6.5: x(2, 4) = 2# x(3, 2) = 12#: x(3, 3) = 3.6: x(3, 4) = 3# x(4, 2) = 12.6: x(4, 3) = 5#: x(4, 4) = 4# x(5, 2) = 13.4: x(5, 3) = 6.1: x(5, 4) = 5# x(6, 2) = 4.1: x(6, 3) = 5.9: x(6, 4) = 6# x(7, 2) = 14.5: x(7, 3) = 6.7: x(7, 4) = 7# x(8, 2) = 11.3: x(8, 3) = 6.5: x(8, 4) = 8# Rem cot 1=1 For i = To dong For j = To x(i, j) = Next j, i Print Print " So lieu vao ma tran X" Print For i = To dong For j = To cot Print " x("; i; ","; j; ") ="; x(i, j); Next j Print: Next i Print Print " So lieu vao vecto y" For i = To dong Print "y("; i; ")="; y(i) 'Print "y" & i & " = "; y(i) Next i Rem chuyen vi ma tran x For i = To dong For j = To cot xt(j, i) = x(i, j) Next j, i 'Print For i = To cot 120 For j = To dong 'Print xt(i, j); Next j 'Print Next i For i = To cot For j = To cot a(i, j) = For d = To dong a(i, j) = a(i, j) + xt(i, d) * x(d, j) Next d, j, i 'Print For i = To cot For j = To cot 'Print a(i, j); Next j 'Print Next i Rem cho ma tran don vi p1(1, 1) = 1: p1(1, 2) = 0: p1(1, 3) = 0: p1(1, 4) = p1(2, 1) = 0: p1(2, 2) = 1: p1(2, 3) = 0: p1(2, 4) = p1(3, 1) = 0: p1(3, 2) = 0: p1(3, 3) = 1: p1(3, 4) = p1(4, 1) = 0: p1(4, 2) = 0: p1(4, 3) = 0: p1(4, 4) = For kkk = To cot For i = To cot For j = To cot km(i, j) = a(i, j) p(i) = p1(kkk, i) Next j, i Rem GIAI HE PHUONG TRINH DAI SO TUYEN TINH mhe = cot For k = To mhe b(k) = p(k): Rem doi don vi Tan, Tan.m Next k k=2 s=1 n=1 101 If (km(n, n) = 0!) Then GoTo 301 601 de = km(n, n) For i = k To mhe da = km(i, n) / de b(i) = b(i) - b(n) * da For j = s To mhe km(i, j) = km(i, j) - km(n, j) * da Next j, i k=k+1 s=s+1 n=n+1 If (k = mhe + 1) Then GoTo 901 GoTo 101 301 n1 = n + 201 For j = To mhe ac = km(n, j) 121 km(n, j) = km(n1, j) km(n1, j) = ac Next j ac = b(n) b(n) = b(n1) b(n1) = ac n1 = n1 + If (km(n, n) = 0!) Then GoTo 9001 GoTo 9101 9001 If (n1 = mhe + 1) Then GoTo 901 9101 If (km(n, n) = 0) Then GoTo 201 GoTo 601 901 For i = To mhe If (km(i, i) = 0) Then GoTo 701 Next i xo(mhe) = b(mhe) / km(mhe, mhe) k1 = mhe k2 = mhe - For i = To k2 k4 = mhe - i r = 0! For j = k1 To mhe r = r + xo(j) * km(k4, j) Next j xo(k4) = (b(k4) - r) / km(k4, k4) k1 = k1 - Next i Rem in vec to nghiem xo(i) For i = To mhe 'Print " xo(" & i & ") = "; xo(i) Next i 'Print GoTo 150 701 uu$ = " dinh thuc bang khong" Print , uu$ Print , " He phuong trinh vo nghiem" GoTo 151 150 ' Print 151 'Print For j = To mhe ad(j, kkk) = xo(j) Next j For j = To mhe 'Print ad(j, kkk) Next j Next kkk 'Print " ma tran nghich dao (Xt*X)-1" 'Print For i = To cot For j = To cot 'Print " " & Format(ad(i, j), "#######.00000"); Next j 'Print 122 Next i For i = To cot For j = To cot tich(i, j) = For d = To cot tich(i, j) = tich(i, j) + a(i, d) * ad(d, j) Next d, j, i 'Print 'Print " In tich ma tran (Xt*X)-1 va matran (Xt*X phai bang don vi " For j = To cot For i = To cot 'Print " " & Format(tich(j, i), "##.0"); Next i 'Print Next j Rem NHAN ADAO VOI XT DE THU XEM CO BANG MA TRAN DON VI KHONG 'Print "In tich ma tran (Xt*X)-1 va ma tran Xt" For i = To cot For j = To dong nhan1(i, j) = For d = To cot nhan1(i, j) = nhan1(i, j) + ad(i, d) * xt(d, j) Next d, j, i 'Print For ik = To cot For jk = To dong 'Print " " & Format(nhan1(ik, jk), "0.00000"); Next jk 'Print Next ik 'Print For i = To cot hesob(i) = For d = To dong hesob(i) = hesob(i) + nhan1(i, d) * y(d) Next d, i Print Rem ham hoi quy la chieu cao song H=y Rem Chon diem T=6 (sec); huong=3 ; tháng=6 , ta co: Hmax=4,9 m voi T=6 sec Print " In he so b cua ham hoi quy " For i = To cot Print "b" & i & "="; Format(hesob(i), " 0.0000") & " 1/m" Next i Rem bai cua truong Print Print " Ham hoi quy" Print "y = b(1) + b(2) * x(2)+ b(3) * x(3) + b(4) * b(4)" yy = hesob(1) + hesob(2) * + hesob(3) * 3# + hesob(4) * Print " T = "; & " sec" Print " H = " & Format(yy, " 0.00") & " met" 'Print "y = b(1) + b(2) * + b(3) * + b(4) * 6" Rem DOAN TINH SO SONG k t = 6: 't la chu ky 123 omega = * 3.1416 / t Print "omega ="; omega; " (1/m)" buoc = kk0 = 0.1117 ee = 0.0001 'Print "khi kh=1 thi k = " & Format(kk0, " 0.0000") 'Print " da nhan k= " & Format(kk0, " 0.0000") 100 h = 15# kh = kk0 * h t1 = Exp(kh) t2 = / t1 ssh = (t1 - t2) / cch = (t1 + t2) / tth = ssh / cch kk = omega * omega / (9.81 * tth) 'Print "buoc "; buoc; ": kh = " & Format(kh, " 0.000"); " th = " & Format(tth, " 0.0000"); " k = " & Format(kk, " 0.0000") If Abs(kk - kk0) < ee Then GoTo 199 kk0 = kk buoc = buoc + 'Print " da nhan k= " & Format(kk, " 0.0000") GoTo 100 199 nghiem = kk Print "lamda ="; Format(2 * 3.1416 / nghiem, " 0.0000") & " m" 'Print "lamda ="; * 3.1416 / nghiem Rem KHI VUNG NUOC CO DO SAU h LON nghia la kh >3.14 thi lay k = omega * omega / 9.81 Print kk = omega * omega / 9.81 Print " Voi vung nuoc sau lay k = " & Format(kk, " 0.0000") Print " k= " & Format(kk, " 0.0000") & " 1/m" Rem NUOC SAU THI VAN TOC GIAM NHANH THEO DO SAU omega = Sqr(9.81 * k) Print " omega= " & Format(omega, " 0.0000") End Sub 124 da tinh PHỤ LỤC CHƯƠNG TRÌNH TÍNH TẢI TRỌNG SĨNG Private Sub HScrtanso_Change() 'Private Sub Form_Load() Rem giai phuong trinh bac Rem CHUONG TRINH VE DAO DONG KHUNG CONG TRINH BIEN Static a(4, 4), ad(4, 4), tich(10, 10), p(10), p1(4, 4) As Double Static km(4, 4), b(10), xo(10), y(10), dao(4, 4), xt(10, 10) As Double Static nhan1(4, 8), nhan2(10, 10), hesob(10), x(8, 4), xxdd(1000), xxra(1000) As Double Show Scale (-320, 240)-(320, -240) Rem bai cua truong Rem DOAN TINH SO SONG k t = 6: 't la chu ky omega = * 3.1416 / t 'Print "omega ="; omega; " (1/m)" buoc = kk0 = 0.1117 ee = 0.0001 'Print "khi kh=1 thi k = " & Format(kk0, " 0.0000") 'List1.AddItem " da nhan k= " & Format(kk0, " 0.0000") 100 h = 15# kh = kk0 * h t1 = Exp(kh) t2 = / t1 ssh = (t1 - t2) / cch = (t1 + t2) / tth = ssh / cch kk = omega * omega / (9.81 * tth) 'Print "buoc "; buoc; ": kh = " & Format(kh, " 0.000"); " th = " & Format(tth, " 0.0000"); " k = " & Format(kk, " 0.0000") If Abs(kk - kk0) < ee Then GoTo 199 kk0 = kk buoc = buoc + 'List1.AddItem " da nhan k= " & Format(kk, " 0.0000") GoTo 100 199 nghiem = kk 'Print "lamda ="; * 3.1416 / nghiem Rem KHI VUNG NUOC CO DO SAU h LON nghia la kh >3.14 thi lay 'omega * omega =gk 'k = omega * omega / 9.81 Print kkk = omega * omega / 9.81 'Print " Voi vung nuoc sau lay k = " & Format(kkk, " 0.0000") List1.AddItem " Voi vung nuoc sau lay " List1.AddItem " k= " & Format(kkk, " 0.0000") & " 1/m" 125 da tinh Rem NUOC SAU THI VAN TOC GIAM NHANH THEO DO SAU k = 0.1118 ro = 1.024: cvt = 1: cqt = 2: d1 = 30#: d2 = 15#: h = 110#: h0 = 4.9 Dim fx(200), mx(200) ro = 1.024: cvt = 1: cqt = 2: d1 = 30#: d2 = 15#: h = 110#: h0 = 4.9 dentad = (d1 - d2) / tank1 = * (h + h0 / 2) / (d1 - d2) tank2 = / tank1 'List2.AddItem " k= " & Format(k, " 0.0000") & " 1/m" 'List2.AddItem " tang1= " & Format(tank1, " 0.0000") 'List2.AddItem " tang2= " & Format(tank2, " 0.0000") omega = Sqr(9.81 * k) 'List2.AddItem " omega= " & Format(omega, " 0.0000") Dim omet(100), hopfx(100), hopmx(100) hmax = h + (h0 / 2) Texttanso.Text = HScrtanso.Value Rem VE DONG VE DONG VE DONG VE DONG VE DONG Cls: Rem SAU CLS CO CHAM THI LUU LAI DO THI TAN SO j = HScrtanso.Value: 'Mo j 'For j = To 15:' dong for Rem VE TINH VE TINH VE TINH VE TINH VE TINH VE TINH VE TINH 'j = HScrtanso.Value:'dong j 'For j = To 15: 'mo for omet(j) = omega * (6 + j * 0.2) List2.AddItem " " & " omega(" & j & ") = " & Format(omet(j), "0.00") 'eta = (h0 / 2) * Cos(omet(j)): ' o dinh song Cos(omega)=0 zy = h + (h0 / 2) * Cos(omet(j)): ' tinh o dinh song, o day bien y=0 u1 = Exp(k * h) u2 = Exp(-k * h) ss = (u1 - u2) / hs = (omega * h0 / 2) / ss 'List2.AddItem " hsvx= " & Format(hs, " 0.000000000") hsax = (omega * omega * h0 / 2) / ss 'List2.AddItem " hsax= " & Format(hsax, " 0.000000000") For i = To hmax dtt = (h + (h0 / 2) * Cos(omet(j)) - i) * tank2 d = d2 + * dtt r1 = Exp(k * i) r2 = / r1 cc = (r1 + r2) / vx = hs * cc * Cos(omet(j)) ax = hs * omega * cc * Sin(omet(j)) Rem in vx de so sanh voi vung nuoc nong gan bo 'List2.AddItem " vx= " & i & " " & Format(vx, " 0.0000") & " ax= " & Format(ax, " 0.0000") fx(i) = 0.5 * ro * cvt * d * Abs(vx) * vx + ro * cqt * (3.1416 * d) * (3.1416 * d) * ax / mx(i) = fx(i) * (i - 0.5) 'List2.AddItem " " & (i) & " D=" & Format(d, "0.0") & " fx=" & Format(fx(i), " 0.0000") Scale (-320, 240)-(320, -240) xphai = 30 Line (-xphai, 0)-(-xphai, 115) Rem ve vx de so sanh voi vung nuoc nong gan bo 'PSet (vx * - xphai + 20, i), QBColor(15) 126 Rem VE CHO TRUONG HOP MAX la j=7 If j = Then GoTo 222 DrawWidth = Rem ve tinh tinh tinh 'PSet (fx(i) / 50 - xphai, i), QBColor(15): 'mo lenh ve Rem ve dong ve dong PSet (-fx(i) / 50 - xphai, i), QBColor(15): 'mo lenh ve 'Line (-fx(i) / 50 - xphai, i)-(0 - xphai, i) GoTo 333 222 DrawWidth = PSet (-fx(i) / 50 - xphai, i), QBColor(12) DrawWidth = 'Line (-fx(i) / 50 - xphai, i)-(0 - xphai, i), QBColor(11) 333 Next i hopfx(j) = hopmx(j) = For i = To hmax hopfx(j) = hopfx(j) + fx(i) hopmx(j) = hopmx(j) + mx(i) Next i DrawWidth = Dim omegao(111) omegao(j) = omet(j) PSet (omet(j) * 25 - 50 - * xphai, hopfx(j) / 700), QBColor(11) Rem doi kN ve Tan (chia cho 10) List1.AddItem " " & " omega(" & j & ") = " & Format(omet(j), "0.00") List1.AddItem " Fx =" & Format(hopfx(j) / 10, "0.00") & " Tan, " List1.AddItem " mz =" & Format(hopmx(j) / 10, "0.00") & " Tm" tay = hopmx(j) / hopfx(j) List1.AddItem " b =" & Format(tay, "0.00") & " m" Dim xdr(100), yd2(100) List1.AddItem "============================= " ome = omet(7) Rem tinh cho truong hop Fx cuc dai For j = To omet(j) = ome 'eta = (h0 / 2) * Cos(omet(j)): ' o dinh song Cos(omega)=0 zy = h + (h0 / 2) * Cos(omet(j)): ' tinh o dinh song, o day bien y=0 u1 = Exp(k * h) u2 = Exp(-k * h) ss = (u1 - u2) / hs = (omega * h0 / 2) / ss 'List2.AddItem " heso= " & Format(hs, " 0.000000000") For i = To hmax dtt = (h + (h0 / 2) * Cos(omet(j)) - i) * tank2 d = d2 + * dtt r1 = Exp(k * i) r2 = / r1 127 cc = (r1 + r2) / vx = hs * cc * Cos(omet(j)) ax = hs * omega * cc * Sin(omet(j)) 'List2.AddItem " vx= " & Format(vx, " 0.0000") & " ax= " & Format(ax, " 0.0000") fx(i) = 0.5 * ro * cvt * d * Abs(vx) * vx + ro * cqt * (3.1416 * d) * (3.1416 * d) * ax / mx(i) = fx(i) * (i - 0.5) List2.AddItem " y= " & i & ", D=" & Format(d, "0.0") & ", fx=" & Format(fx(i) / 10, " 0.0000") Scale (-320, 240)-(320, -240) xphai = 30 DrawWidth = 'PSet ((-fx(i)) / 50 - xphai, i), QBColor(12) DrawWidth = Line (-fx(i) / 50 - xphai, i)-(0 - xphai, i) Next i hopfx(j) = hopmx(j) = For i = To hmax hopfx(j) = hopfx(j) + fx(i) hopmx(j) = hopmx(j) + mx(i) Next i Rem doi kN ve Tan (chia cho 10) 'List1.AddItem " Fx =" & Format(hopfx(j) / 10, "0.00") & " Format(omet(j), "0.00") 'List1.AddItem " mz =" & Format(hopmx(j) / 10, "0.00") & " Tm" tay = hopmx(j) / hopfx(j) 'List1.AddItem " b =" & Format(tay, "0.00") & " m" Next j 'List1.AddItem " ome =" & Format(ome, "0.00") DrawWidth = PSet (0 - xphai, tay), QBColor(12) List2.AddItem "============================= " DrawWidth = End Sub 128 Tan " & " omegat= " & PHỤ LỤC CHƯƠNG TRÌNH VẼ MIỀN ĐỒNG MỨC CỦA CHUYỂN VỊ Private Sub Command1_Click() End End Sub Private Sub Form1_Click() End End Sub Private Sub Hscrdd_Change() Scale (-320, 240)-(320, -240) Show 'Cls Rem CHI VE DEN d=60 (voi d=60 thi chuyen vi=0 tren truc dung) DE LAY DAU CUA CHUYEN VI Rem NEU d=61 THi CHUYEN Vi BAT DAU DOi DAU Dim vte(1000), uer(1000), urq(1000), tx1(1000), tx2(1000), dx1(1000), dy1(1000), dx2(1000), dy2(1000) As Double Dim r(1000), rr(1000), ttx(1000), tty(1000), x0(1000), y0(1000) As Double Dim urp(1000), urch(1000), urt(1000), urpp(1000), xima(900) Dim tetaq(1000), tetap(1000) Open "d:\cv18.txt" For Output As #1 Textdd.Text = HScrdd.Value jd = HScrdd.Value + 0: 'Gan j de viet goc bang dd = jd: Rem Dieu kien duong kinh d