Header Page of 258 -i - Lời cảm ơn Tác giả bày tỏ lòng biết ơn sâu sắc tới tập thể thầy hướng dẫn khoa học, PGS.TS Đinh Quang Cường GS.TS Phạm Khắc Hùng, tận tâm hướng dẫn giúp đỡ tác giả để hoàn thành luận án Ngoài kiến thức khoa học quý báu, thầy động viên, quan tâm hỗ trợ để tác giả vượt qua nhiều thời điểm khó khăn trình thực luận án Tác giả xin đặc biệt cảm ơn GS TS Phạm Khắc Hùng cho phép tác giả vận dụng phần sáng chế để giải vấn đề luận án Tác giả chân thành cảm ơn đồng nghiệp Viện Xây dựng Công trình Biển, cán Khoa Sau Đại học trường Đại học Xây Dựng đóng góp ý kiến chuyên môn tạo điều kiện tốt để tác giả hoàn thành luận án Tác giả cảm ơn gia đình yêu quý mình, đặc biệt vợ, cha mẹ hai bên nội ngoại tin tưởng, khích lệ, cảm thông cho tác giả năm tháng làm luận án Tác giả Mai Hồng Quân Footer Page of 258 Header Page of 258 -ii - Lời cam đoan Tôi xin cam đoan công trình nghiên cứu riêng Các số liệu, kết nêu luận án trung thực chưa công bố công trình khác Ngày … tháng… năm 2014 Nghiên cứu sinh Mai Hồng Quân Footer Page of 258 Header Page of 258 -iii - Mục lục Lời cảm ơn i Danh mục chữ viết tắt ký hiệu ix Danh mục hình vẽ xii Danh mục bảng biểu xiv Mở đầu Ý nghĩa khoa học thực tiễn luận án Ý nghĩa khoa học: Ý nghĩa thực tiễn luận án: CHƯƠNG 1: TỔNG QUAN VỀ CÁC PHƯƠNG PHÁP ĐÁNH GIÁ AN TOÀN KẾT CẤU CHÂN ĐẾ CÔNG TRÌNH BIỂN CỐ ĐỊNH BẰNG THÉP (JACKET) 1.1 Quá trình phát triển xây dựng công trình biển cố định thép 1.1.1 Khái quát công trình biển cố định thép 1.1.2 Các tải trọng tác động lên công trình biển cố định 1.1.3 Yêu cầu thiết kế thi công 1.1.4 Quá trình phát triển xây dựng công trình biển cố định thép giới 1.1.5 Tình hình ứng dụng triển vọng phát triển loại công trình biển cố định thép để khai thác dầu khí Việt Nam 1.2 Tình hình nghiên cứu ứng dụng phương pháp đánh giá an toàn kết cấu jacket tiêu chuẩn hành 1.2.1 Các phương pháp đánh giá an toàn sử dụng tiêu chuẩn hành 1.2.1.1 Phương pháp đánh giá theo trạng thái giới hạn 1.2.1.2 Phương pháp đánh giá theo độ tin cậy 1.2.2 Nhận xét phương pháp đánh giá an toàn sử dụng tiêu chuẩn hành 10 1.2.3 Tình hình nghiên cứu nước quốc tế 10 1.2.3.1 Tình hình nghiên cứu giới 11 1.2.3.2 Tình hình nghiên cứu nước 11 1.3 Đặt vấn đề nghiên cứu luận án 12 Footer Page of 258 Header Page of 258 -iv - 1.3.1.1 Nguyên lý tổng quát để đánh giá an toàn loại kết cấu công trình biển theo sáng chế GS Phạm Khắc Hùng 12 1.3.2 Đặt vấn đề nghiên cứu 13 1.3.3 Nhiệm vụ nghiên cứu luận án 13 1.4 Các giả thiết giới hạn nghiên cứu luận án 14 1.5 Kết luận chương 15 CHƯƠNG 2: ĐÁNH GIÁ AN TOÀN CỦA KẾT CẤU JACKET CÔNG TRÌNH BIỂN CỐ ĐỊNH DỰA TRÊN TÍNH TOÁN BỀN VÀ MỎI TRUYỀN THỐNG 2.1 Mở đầu 16 2.1.1 Các trạng thái tác động sóng lên công trình biển 16 2.1.2 Đánh giá an toàn kết cấu theo tiêu chuẩn hành 17 2.1.2.1 Đánh giá theo điều kiện bền 17 2.1.2.2 Đánh giá an toàn theo điều kiện mỏi 17 2.2 Mô tả chuyển động sóng biển bề mặt 18 2.2.1 Mô tả sóng theo quan điểm tiền định 18 2.2.1.1 Các lý thuyết sóng 18 2.2.1.2 Miền áp dụng lý thuyết sóng 19 2.2.2 Mô tả sóng theo quan điểm ngẫu nhiên 19 2.2.2.1 Mặt cắt (profile) sóng ngẫu nhiên 19 2.2.2.2 Phổ lượng sóng 20 2.2.2.3 Các phổ sóng thông dụng thiết kế kết cấu công trình biển 21 2.2.2.4 Phổ vận tốc gia tốc phần tử nước sóng ngẫu nhiên 22 2.3.Tải trọng sóng tác dụng lên phần tử mảnh kết cấu jacket 22 2.3.1 Tải trọng sóng tiền định 22 2.3.2 Tải trọng sóng ngẫu nhiên 23 2.4 Đánh giá an toàn kết cấu jacket theo điều kiện bền truyền thống 24 2.4.1 Đánh giá an toàn kết cấu dựa mô hình sóng tiền định 24 2.4.1.1 Xác định phản ứng động kết cấu theo mô hình tiền định 24 Footer Page of 258 Header Page of 258 -v - 2.4.1.2 Kiểm tra bền kết cấu theo mô hình tiền định 25 2.4.2 Đánh giá an toàn kết cấu dựa mô hình sóng ngẫu nhiên 25 2.4.2.1 Phương pháp phổ 25 2.4.2.2 Phương pháp giải miền thời gian 26 2.5 Đánh giá an toàn kết cấu jacket theo điều kiện mỏi truyền thống 28 2.5.1 Đánh giá an toàn mỏi kết cấu dựa mô hình sóng tiền định 28 2.5.1.1 Tính toán mỏi theo phương pháp tổn thất tích luỹ 28 2.5.1.2 Tính toán tổn thất mỏi theo mô hình sóng tiền định 29 2.5.2 Đánh giá an toàn mỏi kết cấu dựa mô hình sóng ngẫu nhiên 30 2.5.2.1 Ứng suất ngẫu nhiên điểm nóng 30 2.5.2.2 Xác định tổn thất mỏi trung bình điểm nóng trạng thái biển ngắn hạn phương pháp phổ 31 2.5.2.3 Tuổi thọ mỏi trung bình điểm nóng kết cấu jacket 33 2.6 Kết luận chương 33 CHƯƠNG 3: PHƯƠNG PHÁP LUẬN ĐÁNH GIÁ MỨC SUY GIẢM ĐỘ TIN CẬY THEO THỜI GIAN CỦA KẾT CẤU JACKET CÁC CÔNG TRÌNH BIỂN CỐ ĐỊNH 3.1 Mở đầu 35 3.1.1 Dạng tổng quát độ tin cậy theo điều kiện bền truyền thống 35 3.1.2 Dạng tổng quát độ tin cậy theo điều kiện mỏi truyền thống 35 3.2 Dạng tổng quát đánh giá suy giảm ĐTC tổng thể kết cấu jacket 35 3.2.1 Dạng tổng quát độ độ tin cậy dựa điều kiện bền mở rộng 35 3.2.2 Dạng tổng quát độ độ tin cậy theo điều kiện mỏi mở rộng 36 3.3 Xác định độ tin cậy bền kết cấu jacket trạng thái biển ngắn hạn cực đại 37 3.3.1 Ứng suất ngẫu nhiên kết cấu 37 3.3.2 Độ tin cậy bền kết cấu jacket ứng suất có phổ dải hẹp 38 3.3.3 Độ tin cậy bền kết cấu jacket ứng suất có phổ dải rộng 38 Footer Page of 258 Header Page of 258 -vi - 3.4 Xác định độ tin cậy theo điều kiện mỏi điểm nóng kết cấu jacket phụ thuộc vào thời gian khai thác 40 3.4.1 Xác định kỳ vọng phương sai tổn thất mỏi trạng thái biển ngắn hạn điểm nóng 40 3.4.1.1 Biểu diễn ứng suất ngẫu nhiên điểm nóng miền thời gian 40 3.4.1.2 Xác định số lượng chu trình ứng suất kỹ thuật đếm dòng mưa 43 3.4.1.3 Xác định kỳ vọng phương sai tổn thất mỏi điểm nóng trạng thái biển ngắn hạn 44 3.4.2 Xây dựng hàm phân phối xác suất tổn thất mỏi điểm nóng năm……… 48 3.4.2.1 Mật độ xác suất tổn thất mỏi điểm nóng năm 48 3.4.2.2 Hàm phân phối xác suất tổn thất mỏi điểm nóng năm 48 3.4.2.3 Đánh giá độ tin cậy mỏi điểm nóng năm………………… 49 3.5 Xác định độ tin cậy mỏi điểm nóng thời điểm T(năm) 50 3.6 Xác định ĐTC kết cấu điểm xét dựa điều kiện bền mở rộng 51 3.6.1 Độ tin cậy ứng với trường hợp ứng suất kết cấu có phổ dải hẹp 51 3.6.2 Độ tin cậy ứng với trường hợp ứng suất kết cấu có phổ dải rộng 51 3.7 Xác định ĐTC kết cấu điểm xét dựa điều kiện mỏi mở rộng 51 3.7.1 Kỳ vọng tỷ số tổn thất mỏi mở rộng 51 3.7.2 Phương sai tỷ số tổn thất mỏi mở rộng 52 3.7.3 Độ tin cậy tính theo điều kiện mỏi mở rộng 52 3.8 Đánh giá suy giảm ĐTC theo thời gian KCCĐ jacket dựa ĐTC thực tế (tổng thể) KC điểm xét 52 3.9 Đánh giá mức độ suy giảm khả chịu tải điều kiện biển cực đại theo thời gian khai thác công trình 53 3.10 Sơ đồ thuật toán đánh giá an toàn kết cấu jacket theo phương pháp truyền thống theo phương pháp luận Luận án 55 3.11 Kết luận chương 60 Footer Page of 258 Header Page of 258 -vii - CHƯƠNG 4: VÍ DỤ ỨNG DỤNG 4.1 Mở đầu 61 4.2 Các số liệu đầu vào sử dụng ví dụ 61 4.2.1 Số liệu công trình 61 4.2.2 Số liệu môi trường 62 4.3 Các phần mềm máy tính sử dụng ví dụ 64 4.3.1 Các phần mềm thương mại 64 4.3.2 Phần mềm tự lập “ RFCAL” 65 4.4 Kết tính độ tin cậy theo điều kiện bền truyền thống 65 4.4.1 Kết tính nội lực ngẫu nhiên kết cấu 65 4.4.2 Kiểm tra bền phần tử 67 4.5 Độ tin cậy theo điều kiện mỏi truyền thống kết cấu jacket 67 4.5.1 Tính toán tổn thất mỏi 68 4.5.1.1 Đầu vào tính mỏi 68 4.5.1.2 Tính toán ứng suất điểm nóng 69 4.5.1.3 Tính toán tổn thất mỏi điểm nóng 71 4.5.2 Kết tính toán độ tin cậy theo điều kiện phá hủy mỏi truyền thống 73 4.6 Đánh giá suy giảm độ tin cậy tổng thể kết cấu jacket 73 4.6.1 Độ tin cậy điểm đặc trưng kết cấu bắt đầu khai thác 73 4.6.2 Đánh giá suy giảm độ tin cậy khả chịu tải kết cấu trình khai thác 74 4.6.2.1 Độ tin cậy theo điều kiện bền mở rộng 74 4.6.2.2 Kết tính độ tin cậy theo điều kiện mỏi mở rộng 75 4.6.2.3 Nghiên cứu bổ sung với trường hợp biến đổi khí hậu bất thường: 77 4.6.2.4 Kết tính toán độ tin cậy tổng thể công trình 79 4.7 Đánh giá mức độ suy giảm khả chịu tải điều kiện biển cực đại theo thời gian khai thác công trình điểm xét 80 4.8 Kết luận chương 81 Footer Page of 258 Header Page of 258 -viii - PHẦN KẾT LUẬN Những kết đạt 83 Những đóng góp luận án 83 Kiến nghị luận án 84 Hướng nghiên cứu phát triển luận án 84 PHỤ LỤC: Chương trình tính mỏi ngẫu nhiên RFCAL ………… ……….…….90 Footer Page of 258 Header Page of 258 -ix - Danh mục chữ viết tắt ký hiệu CTB Công trình biển CTB CĐ TTB ĐLNN Công trình biển cố định Trạng thái biển Đại lượng ngẫu nhiên QTNN ĐTC Quá trình ngẫu nhiên Độ tin cậy Ứng suất kết cấu ߪ ሾߪ ሿ Ứng suất cho phép ܴ Khả chịu lực vật liệu ߛ௠ ULS WSD LRFD Hệ số an toàn vật liệu ܴ௞ ߛ௜ FLS ‫ܦ‬ ሾ‫ ܦ‬ሿ Cường độ vật liệu Trạng thái ứng suất cực hạn Thiết kế theo ứng suất cho phép Phương pháp thiết kế theo hệ số tải trọng hệ số cường độ Hệ số tải trọng thứ i Trạng thái giới hạn phá hủy mỏi Tỷ số tổn thất mỏi Tỷ số tổn thất mỏi cho phép DT f (DT) Tỷ số tổn thất mỏi tích luỹ thời gian T Mật độ xác suất DT ߛ am N C, m Hệ số an toàn mỏi Chiều sâu vết nứt số chu trình ứng suất theo đường cong S-N C, m thông số theo quy tắc Paris ∆௄ ܲ, ሾܲሿ Hệ số cường độ ứng suất điểm nóng bị phá hủy mỏi ߚ௧௞ , ሾߚ ሿ Chỉ số độ tin cậy tính toán số độ tin cậy cho phép ܲ௙ , ൣܲ௙ ൧ Độ tin cậy độ tin cậy cho phép kết cấu Xác suất phá hủy xác suất phá hủy cho phép ܲ஻ Độ tin cậy theo điều kiện bền ܲ௠் Độ tin cậy theo điều kiện mỏi thời điểm T năm ܲ௠ Footer Page of 258 Độ tin cậy theo điều kiện mỏi Header Page 10 of 258 -x - ܲ஻்௢௧ Độ tin cậy theo điều kiện bền tổng thể ߟ(௧) Mặt cắt bề mặt sóng biển (Profile sóng ) Tz Tc H, Hs, Hmax ܵఎఎ (߱) Chu kỳ trung bình cắt không Chu kỳ trung bình đỉnh sóng cực đại (Hmax ) Chiều cao sóng, chiều cao sóng đáng kể, chiều cao sóng cực đại ܲ௠்௢௧ ߝ Độ tin cậy theo điều kiện mỏi mở rộng Phổ sóng Thông số bề rộng dải phổ Sv v (ω ) x x Phổ vận tốc phần tử nước theo phương x phương y Sv v ( ω ) y y Sa a ( ω ) x x Phổ gia tốc phần tử nước theo phương x phương y S a a (ω ) y y F (t) Tải trọng sóng tác động lên kết cấu tính theo công thức Morison FD(t) FI(t) Thành phần lực cản tải trọng sóng Thành phần quán tính tải trọng sóng Mật độ nước biển ߩ CD Hệ số cản vận tốc ‫ݒ‬, ‫ݒ‬ሶ Hệ số nước kèm Hệ số quán tính Vận tốc gia tốc phần tử nước ‫ݑ‬, ‫ݑ‬ሶ , ‫ݑ‬ሷ Chuyển vị, vận tốc gia tốc vật cản ω ω1 Ma trận độ cứng, véc tơ chuyển vị nút véc tơ tải trọng quy nút Hệ số khuếch đại động Tần số dao động sóng Tần số dạng dao động riêng thứ kết cấu CM CI ܵிி (߱ ) Phổ tải trọng sóng lên phần tử kết cấu ߜ௩௫ K,U(t), F(t) Kđ Độ lệch chuẩn thành phần vận tốc ‫ݒ‬௫ T1 ε Φ(n*n) Footer Page 10 of 258 Chu kỳ dạng dao động riêng thứ Hệ số cản kết cấu dao động Ma trận dạng dao động riêng kết cấu Header Page 104 of 258 -89 - 19 Bureau Veritas, MCS, “ARIANE-3D V6.3 – User’s Guide” 20 Barltrop N.D.P and Adams, A.J, “Dynamics of fixed marine structures”, 3rd Edition Butterworth Heinemann/MTD Ltd, 1991 21 Bathe K.J, “Finite Element Procedures in Engineering Analysis”, 1996 22 Bea, R.G , “Reliability-based Requalification Criteria for Offshore Platforms” , Proc 12th OMAE Vol 11: 351-361 New York: ASME, 1993 23 Bea, R.G and Young, C.N., “Loading and Capacity Effects on Platform Performance in Extreme Storm Waves and Earthquakes”, Proc 25th OTC: OTC 7140 Houston: Offshore Technology Conf 1993 24 Cornell, C.A., “Structural Reliability - Some contributions to Offshore Technology” , Proc 27th OTC: 535-542, OTC 7753 Houston: Offshore Technology Conf, 1995 25 DNV OSS 101, “Rules for Classification of the Offshore Drilling and Support Units”, 2012 26 DNV OS 201, “Structure Design of Offshore Units”, 2012 27 DNV, “Guidline for offshore structure reliability-General”, DNV report No 952018, 1996 28 DNV, “Guidline for offshore structure reliability-Application to jacket Platform”, DNV report No 95-3203, 1996 29 D Benasciutti and R Tovo, “Spectral methods for lifetime prediction under wide-band stationary random processes ”, University of Ferrara, Italy, 2002 30 HSE “A review of Reliability considerations for Fixed Offshore platform”, John Wiley & Son, Inc, 2000 31 Energo Engineering, “Assessment of fixed offshore platform performance in hurricanes katrina and rita, final report”, Inc.3100 Wilcrest Drive, Suite 240, Houston, Texas 77042, 2007 32 Guedes Soares & Das, “Analysis and Design of Marine Structures”, Taylor & Francis Group, London, ISBN 978-0-415-54934-9, 2009 Footer Page 104 of 258 Header Page 105 of 258 -90 - 33 Gerhard Ersdal, “Assessment of existing offshore structures for life extension”, Doctorial Thesis, Faculty of Science and Technology, Department of Mechanical and Structural Engineering and Material Science, 2005 34 Heideman, L.C & Weaver T.O., “ Static Wave Force Procedure for Platform Design”, Proc Civil Engng in the Oceans: 496-517.College station, Texas: ASCE, 1992 35 ISO 19901-6, Specific requirements for offshore structures – Part 6: Marine Operations, 1st Edition December 2009 36 ISO 19901-2, Specific requirements for offshore structures – Part 2: Seismic Design Procedures and Criteria, 1st Edition November 2004 37 Isaac Elishakoff, “ Probability Methods in the theory of Structures”, John Willey & Son, USA, 1983 38 James F.Wilson , “Dynamics of offshore structures”, John Wiley & Son, Inc, 2003 39 James F Wilson, “Dynamics of offshore structures” ISBN 0-471-26467-9, 1984 40 Krieger, W.F et al., “Process for Assessment of Existing Platforms to Determine their Fitness for Purpose”, Proc 20th OTC: OTC7482 Houston: Offshore Technology Conference, 1994 41 LS Etube, “Fatigue and Fracture Mechanics of Offshore Structures”, Professional Engineering Publishing Limited, London and Bury St Edmunds, UK, 2001 42 Lloyd, J.R & Karsan, “ Development of a Re liability-based Alternative to API RP2A.Proc 20th”, OTC Vol.4:593-600, OTC 5882 Houston: Offshore Technology Conf, 1988 43 Moan, T., “Review of Probabilistic Inspection Analysis Methods”, Offshore Technology Report OTO 1999061, Health and Safety Executive, UK, 1999 44 Moan, T., Johannesen, J.M and Vårdal, O.T, “Probabilistic Inspection Planning of jacket Structures”, Proc 31st OTC: OTC 10848 Houston:Offshore Tecnhology Conf, 1999 45 Minoo H.Patel, “Dynamics of offshore structures”, Bulter&Tanner, Froem, Somerset, 1989 Footer Page 105 of 258 Header Page 106 of 258 -91 - 46 Moan, T., “Target Levels for Structural Reliability and Risk Analysis of Offshore Structures”, C Guedes Soares (ed), Rotterdam: A.A Balkema, 1998 47 Michael E McCormick, “Ocean Engineering Mechanics”, United States Naval Academy, 1998 48 N Haritos, “Introduction to the Analysis and Design of Offshore Structures– An Overview”, The University of Melbourne, Australia, 2007 49 Moan, T., “Safety of Offshore Structures”, Proc 4th ICASP: 41-85 Bologna: Pitagora Editrice, 1983 50 NORSOK Standard N-001, Integrity of Offshore Structures, 7th Edition June 2010 51 NORSOK Standard N-004, Design of Steel Structures, 2nd Edition October 2004 52 NORSOK Standard N-006 , Assessment of Structure Integrity for Existing Offshore Load-bearing Structures, 1st Edition March 2009 53 Nallayarasu, “Offshore structures – Analysis and Design”, Indian Institute of Technology Madras, 2000 54 Naser SHABAKHTY, “Durable Reliability of Jack-up Platforms- The impact of Fatigue, Fracture and Effect of Extreme Environmental Loads on the Structural Reliability”, Technische Universiteit Delft, The Netherlands, 2004 55 Palle Thoft-Christensen, Michael J.Baker, “Structural Reliability Theory and Its Applications”, Springer-Verlag Berlin Heidelberg New York, 1982 56 Subrata Chakrabarti, “Handbook of Offshore Engineering”, Elsevier,Ltd, 2005 57 Subrata K Chakrabarti, “Handbook of offshore engineering”, Volume I, Offshore Structure Analysis, Inc., Plainfield, llinois, USA, 2005 58 Subrata K Chakrabarti, “Handbook of offshore engineering”, Volume II, Offshore Structure Analysis, Inc., Plainfield, llinois, USA, 2005 59 SESAM User manual, “Profast— Probabilistic Fatigue Analysis”, Det Norske Veritas, Norway, 2002 60 SACS Version 5.1 Manual report Footer Page 106 of 258 Header Page 107 of 258 -92 - 61 Singiresu S Rao, “Reliability- Based design”, McGraw Hill, Inc, 1992 62 Turan Dirlik, “Application of computer in fatigue analysis”, Doctor thesis, University of Warwick Coventry, England, 1985 63 Yu L., Zhou M., Das P K., “Fatigue reliability of bulk carrier – A case study”, Advanced Structural Reliability Analysis Network, Edinburgh, 2010 64 Yu, L., Das, P K., & Zheng, Y L, “A response surface approach to fatigue reliability of ship structures”, Ship and Offshore Structures, 4(3), 253-259, 2009 65 Yu, L., & Das, P K., “Fatigue design assessment based on pseudo-excitation method”, International Maritime Association of Mediterranean (IMAM), 2007 66 Yu, L., Das, P K., & Barltrop, N D, “A new look at the effect of bandwidth and non-normality on fatigue damage”, Fatigue & Fracture of Engineering Material & Structures, 27(1), 51-58, 2004 67 Yong Bai, “Marine Structural design”, ELSEVIER, 2003 68 Yu, L., Das, P K., & Zheng, Y L., “Stepwise response surface method and its application to reliability analysis of ship structures”, Journal of Offshore Mechanics and Arctic Engineering, 124(4), 226-230, 2002 69 Yu, L., Das, P K., & Barltrop, N D., “Importance sampling method with kernel density estimate – A general approach”, Proceedings of the 8th International Conference on Structural Safety and Reliability (ICOSSAR 2001), California, 2001 Footer Page 107 of 258 Header Page 108 of 258 -93 - PHỤ LỤC: CHƯƠNG TRÌNH TÍNH MỎI “RFCAL” Phần Giao diện chương trình RFCAL Chương trình RFCAL lập ngôn ngữ Visual Basic để tính toán tổ n thất mỏi cho trạng thái biển năm cho nhiều năm Chương trình tự động tạo số liệu kết nối với SACS 5.3 để tính toán nội lực Một số giao diện chương trình sau: Hình Giao diện chương trình, bước tính toán Footer Page 108 of 258 Hình Đọc số liệu phổ sóng tạo file liệu đầu vào cho SACS 5.3 Hình Đọc, lọc tách tính Hình Tính toán tổn thất mỏi cho toán ứng suất phần tử năm Header Page 109 of 258 -94 - Phần Mã nguồn chương trình RFCAL Attribute VB_Name = "Module1" 'Option Explicit Global Xmax, Ymax, Xmin, Ymin Type EventMouseType Event As Byte End Type Type CommandType CommandSt As String Event As Byte End Type Global Command As CommandType Global EventMouse As EventMouseType Global Tyle Sub Init() Xmax = 10 Ymax = 10 Xmin = -10 Ymin = -10 EventMouse.Event = 0# Command.CommandSt = "" Command.Event = 0# End Sub Sub QScale(ByVal XSmin, ByVal YSmin, ByVal XSmax, ByVal YSmax) Dim dx, dy Dim dwx, dwy Dim Tg Dim TlX, TlY Dim Xc, Yc If Xmax < Xmin Then Tg = Xmax Xmax = Xmin Xmin = Tg End If If Ymax < Ymin Then Footer Page 109 of 258 Header Page 110 of 258 -95 - Tg = Ymax Ymax = Ymin Ymin = Tg End If Xmin = XSmin Ymin = YSmin Xmax = XSmax Ymax = YSmax dx = (Xmax - Xmin) dy = (Ymax - Ymin) Xc = (Xmax + Xmin) / 2# Yc = (Ymax + Ymin) / 2# dwx = Form1.Picture1.Width dwy = Form1.Picture1.Height TlX = dx / dwx TlY = dy / dwy If (TlX = 0) Or (TlY = 0) Then MsgBox ("Loi ty le TlX=0 hoac TlY=0") Exit Sub End If If TlX > TlY Then Tyle = TlX Else Tyle = TlY End If Xmin = Xc - dwx * Tyle / Xmax = Xc + dwx * Tyle / Ymin = Yc - dwy * Tyle / Ymax = Yc + dwy * Tyle / Form1.Picture1.Scale (Xmin, Ymax)-(Xmax, Ymin) End Sub Footer Page 110 of 258 Header Page 111 of 258 -96 - Sub QLine(ByVal X1, ByVal Y1, ByVal X2, ByVal Y2) Form1.Picture1.Line (X1, Y1)-(X2, Y2) End Sub Sub QCircle(ByVal Xc, ByVal Yc, ByVal R) Form1.Picture1.Circle (Xc, Yc), R End Sub Sub QPset(ByVal Xp, ByVal Yp) Form1.Picture1.PSet (Xp, Yp) End Sub Sub QSetColor(ByVal Color) Form1.Picture1.ForeColor = Color End Sub Sub QSetBackColor(ByVal Color) Form1.Picture1.BackColor = Color End Sub Function TachSo(ByVal St As String, ByVal StF1 As String, ByVal StF2 As String) Dim Vt1, Vt2 Dim StTg Vt1 = InStr(1, St, StF1) Vt1 = Vt1 + Len(StF1) Vt2 = InStr(Vt1, St, StF2) StTg = Mid(St, Vt1, Vt2 - Vt1) TachSo = Val(StTg) End Function Function DocSo(Path) Dim FileNumber, FileNumberFXA, FileNumberOut Dim FileName, FileNameOut Dim ii, jj Dim St, St1, St2, St3 Footer Page 111 of 258 FileNumberMXA, FileNumberMYA, Header Page 112 of 258 -97 - Dim StTg(0 To 2) 'Dim MaxSeet Dim X, FXA, MXA, MYA Dim Vt1, Vt2 Dim L FileNumber = FileNumber + FileNumberOut = FileNumber + 100 On Error GoTo 10 Open Path + FileName For Input As #FileNumber On Error Resume Next MaxSeet = L=0 Do FileNameOut = Str(MaxSeet) FileNameOut = Mid(FileNameOut, 2, Len(FileNameOut) - 1) Do FileNameOut = "0" + FileNameOut Loop Until Len(FileNameOut) > For ii = To FileNumberOut = FileNumberOut + Select Case ii Case 1: St1 = ".FXA" Case 2: St1 = ".MXA" Case 3: St1 = ".MYA" End Select Open Path + FileNameOut + St1 For Output As #FileNumberOut Do Input #FileNumber, St L=L+1 Print #FileNumberOut, St If EOF(FileNumber) Then Footer Page 112 of 258 Header Page 113 of 258 -98 - Close (FileNumberOut) Close (FileNumber) DocSo = MaxSeet Exit Function End If Loop Until (InStr(1, St, "") 0) 'Or (EOF(FileNumber)) Close (FileNumberOut) Next ii MaxSeet = MaxSeet + Loop Until (EOF(FileNumber)) 'Or (St = "") Close (FileNumber) DocSo = MaxSeet - 10: Close (FileNumber) MaxSeet = Exit Function End Function Sub DocLai(Path, MaxSeet) Dim FileNumber, FileNumberFXA, FileNumberOut Dim FileName, FileNameOut Dim ii, jj Dim St, St1, St2, St3 Dim StTg(0 To 2) 'Dim MaxSeet FileNumberMXA, FileNumberMYA, Dim X, FXA, MXA, MYA Dim Vt1, Vt2 Dim L Dim D, t, A, W 'Bien de tinh moi Dim Ni, di, Tongdi '================================================== 'Dac trung hinh hoc D = 609 t = 12.7 Footer Page 113 of 258 Header Page 114 of 258 -99 - A = 3.14 * (D ^ - (D - * t) ^ 2) / W = 3.14 / 32 * (D ^ - (D - * t) ^ 4) / D '================================================== FileNumber = FileNumber + FileNumberOut = FileNumber + 1000 For ii = To MaxSeet FileNameOut = Str(ii) FileNameOut = Mid(FileNameOut, 2, Len(FileNameOut) - 1) Do FileNameOut = "0" + FileNameOut Loop Until Len(FileNameOut) > FileNumber = FileNumber + Open Path + FileNameOut + ".FXA" For Input As #FileNumber Open Path + FileNameOut + ".MXA" For Input As #FileNumber + Open Path + FileNameOut + ".MYA" For Input As #FileNumber + 'Mo file de ghi ung suat Open Path + FileNameOut + ".S" For Output As #FileNumber + 'Print #FileNumber + 3, "P" For jj = To Input #FileNumber + jj, St Next jj 'Print #FileNumber + 3, St For jj = To Input #FileNumber + jj, St Next jj 'Print #FileNumber + 3, St For jj = To Input #FileNumber + jj, St Next jj 'Print #FileNumber + 3, St Do Footer Page 114 of 258 Header Page 115 of 258 -100 - For jj = To Input #FileNumber + jj, StTg(jj) If InStr(1, StTg(0), "") = Then Vt1 = InStr(1, StTg(jj), "y=""") Vt2 = InStr(Vt1 + 3, StTg(jj), """") St1 = Mid(StTg(jj), Vt1 + 3, Vt2 - Vt1 - 3) Select Case jj Case 0: FXA = Val(St1) Case 1: MXA = Val(St1) Case 2: MYA = Val(St1) End Select End If Next jj If InStr(1, StTg(0), "") = Then Vt1 = InStr(1, StTg(0), "x=""") Vt2 = InStr(Vt1 + 3, StTg(0), """") St1 = Mid(StTg(0), Vt1 + 3, Vt2 - Vt1 - 2) X = Val(St1) '================================================== 'Tinh ung suat S = FXA * 1000 / A + MXA * 1000000 / W + MYA * 1000000 / W '================================================== St = Format(S, "0.0000") Print #FileNumber + 3, St End If Loop Until (EOF(FileNumber)) Or (InStr(1, StTg(0), "") 0) Close (FileNumber) ' 'Tinh moi FileNumber = FileNumber + Open Path + FileNameOut + ".S" For Input As #FileNumber Close (FileNumber) ' - Footer Page 115 of 258 Header Page 116 of 258 -101 - Close (FileNumber + 1) Close (FileNumber + 2) Close (FileNumber + 3) Kill Path + FileNameOut + ".FXA" Kill Path + FileNameOut + ".MXA" Kill Path + FileNameOut + ".MYA" Next ii End Sub Sub Doc() Dim FileName, FileNumber Dim ii Dim P(1 To 100) Dim MaxFile, MaxSeet Dim St '================================================================ ===== FileNumber = Open FileName For Input As #FileNumber MaxFile = MaxFile + Do Input #FileNumber, St MaxSeet = DocSo(St) If MaxSeet > Then Call DocLai(St, MaxSeet) MaxFile = MaxFile + End If Loop Until EOF(FileNumber) Close (FileNumber) MsgBox ("Xong! Tong so file da conver: " + Str(MaxFile)) End Sub Sub TinhMoi() Dim FileNumber, FileNumberP, FileNumberO Dim Path Dim FileName, FileNameP, FileNameO Footer Page 116 of 258 Header Page 117 of 258 -102 - Dim St, St1, St2, StP Dim Row Dim Huong, Hs, ChuKy, Seed Dim Smax, Smin, S, Ni, di, Tongdi Dim Vt1 k = 1E+15 m=3 FileNumberP = FileNumberP + FileNumber = FileNumberP + 100 Open FileNameP For Input As FileNumberP Row = Do Line Input #FileNumberP, StP 'Tach so lieu Huong = Val(TachSo(StP, "huong", "\")) Hs = Val(TachSo(StP, "Hs=", "\")) ChuKy = Val(TachSo(StP, "To=", "\")) Seed = Val(TachSo(StP, "\0", ".")) Row = Row + Sheet2.Cells(Row, 2) = Huong Sheet2.Cells(Row, 3) = Hs Sheet2.Cells(Row, 4) = ChuKy Sheet2.Cells(Row, 5) = Seed 'Doc file ung suat FileName = StP FileNumber = FileNumber + Open FileName For Input As FileNumber Tongdi = Do Line Input #FileNumber, St Vt1 = InStr(1, St, ",") Footer Page 117 of 258 Header Page 118 of 258 St1 = Mid(St, 1, Vt1 - 1) St2 = Mid(St, Vt1 + 1, Len(St) - Vt1) Smax = Val(St1) Smin = Val(St2) S = Abs(Smax - Smin) Ni = k * S ^ (-m) di = / Ni Tongdi = Tongdi + di Loop Until EOF(FileNumber) Close #FileNumber Sheet2.Cells(Row, 6) = Tongdi Loop Until EOF(FileNumberP) Close (FileNumberP) End Sub Footer Page 118 of 258 -103 - ... CHƯƠNG 3: PHƯƠNG PHÁP LUẬN ĐÁNH GIÁ MỨC SUY GIẢM ĐỘ TIN CẬY THEO THỜI GIAN CỦA KẾT CẤU JACKET CÁC CÔNG TRÌNH BIỂN CỐ ĐỊNH 3.1 Mở đầu 35 3.1.1 Dạng tổng quát độ tin cậy theo điều... cứu suy giảm độ tin cậy kết cấu công trình biển cố định thép theo điều kiện bền thời điểm kết cấu chịu tổn thất mỏi - Nhiệm vụ 3: nghiên cứu xác định độ tin cậy kết cấu công trình biển cố định. .. truyền thống Chương 3: Phương pháp luận đánh giá mức suy giảm độ tin cậy theo thời gian kết cấu jacket công trình biển cố định Chương 4: Ví dụ ứng dụng Kết luận luận án Phụ lục Luận án bao gồm 90