Microsoft Word LA Sauphanbien docx B� GIÁO D�C VÀ �ÀO T�O TR��NG ��I H�C XÂY D�NG MAI H NG QUÂN NGHIÊN C U S� SUY GI�M �� TIN C�Y THEO TH�I GIAN C A CHÂN �� CÔNG TRÌNH BI�N C� ��NH B�NG THÉP DO �NH H�[.]
B TR GIÁO D C VÀ ÀO T O NG I H C XÂY D NG MAI H NG QUÂN NGHIÊN C U S C A CHÂN SUY GI M TIN C Y THEO TH I GIAN CƠNG TRÌNH BI N C NH H NH B NG THÉP DO NG C A T N TH T M I LU N ÁN TI N S K THU T HÀ N I 03-2014 B TR GIÁO D C VÀ ÀO T O NG I H C XÂY D NG MAI H NG QUÂN NGHIÊN C U S C A CHÂN SUY GI M TIN C Y THEO TH I GIAN CƠNG TRÌNH BI N C NH H NH B NG THÉP DO NG C A T N TH T M I LU N ÁN TI N S K THU T Chuyên ngành: Xây d ng cơng trình bi n Mã s : 62 58 02 02-2 (mã s m i: 62580203) T p th hư ng d n khoa h c Th y h ng d n s 1: PG.TS inh Quang C Th y h ng d n s 2: GS TS Ph m Kh c Hùng HÀ N I 03-2014 ng -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 PGS.TS d n giúp inh Quang C tác gi báu, th y ã ng d n khoa h c, ng GS.TS Ph m Kh c Hùng, ã t n tâm h ng hoàn thành lu n án Ngoài ki n th c khoa h c quý ng viên, quan tâm h tr tác gi v t qua th i i m khó kh n trình th c hi n lu n án Tác gi xin c nhi u c bi t c m n GS TS Ph m Kh c Hùng ã cho phép tác gi v n d ng m t ph n sáng ch c a gi i quy t v n lu n án Tác gi chân thành c m n Bi n, cán b Khoa Sau i h c tr v chuyên môn t o i u ki n t t nh t ng nghi p Vi n Xây d ng Cơng trình ng i h c Xây D ng ã óng góp ý ki n tác gi hoàn thành lu n án Tác gi c m n gia ình yêu quý c a mình, c bi t i v i v , cha m hai bên n i ngo i ã tin t !ng, khích l , c m thông cho tác gi nh"ng n m tháng làm lu n án Tác gi Mai H ng Quân -ii - L i cam oan Tôi xin cam oan ây cơng trình nghiên c u c a riêng Các s li u, k t qu nêu lu n án trung th c ch a t#ng c công b b t c công trình khác Ngày … tháng… n m 2014 Nghiên c u sinh Mai H ng Quân -iii - M cl 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 c a lu n án Ý ngh%a khoa h c: Ý ngh%a th c ti&n c a lu n án: CH NG 1: T NG QUAN V K T C U CHÂN CÁC PH NG PHÁP 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 1.1.1 Khái qt v cơng trình bi n c 1.1.2 Các t i tr ng tác ÁNH GIÁ AN TOÀN nh b ng thép 'nh b(ng thép ng lên cơng trình bi n c 'nh 1.1.3 Yêu c u c b n v 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 b(ng thép th 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 thép 'nh b(ng 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 hi n hành 1.2.1 Các ph ng pháp ánh giá an toàn s) d ng tiêu chu*n hi 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 1.2.2 Nh n xét v ph tin c y ng pháp ánh giá an toàn s) d ng tiêu chu*n hi 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 th gi i 11 1.2.3.2 Tình hình nghiên c u n 1.3 tv n c 11 nghiên c u c!a lu"n án 12 -iv - 1.3.1.1 Nguyên lý t+ng quát ánh giá an toàn c a lo i k t c u cơng trình bi n theo sáng ch c a GS Ph m Kh c Hùng 12 1.3.2 tv n nghiên c u 13 1.3.3 Nhi m v nghiên c u c a 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 c!a ch ơng 15 CH NG 2: BI N C ÁNH GIÁ AN TỒN C$A K T C U JACKET CƠNG TRÌNH NH D%A TRÊN TÍNH TỐ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 c a k t c u theo tiêu chu*n hi 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) c a sóng ng u nhiên 19 2.2.2.2 Ph+ n ng l ng c a 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 c a 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 c!a 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 c!a 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 c a 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 c a k t c u theo mơ hình ti n 'nh 24 -v - 2.4.1.2 Ki m tra b n c a k t c u theo mơ hình ti n 'nh 25 2.4.2 ánh giá an toàn c a 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 c!a 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 v m i c a k t c u d a mơ hình sóng ti n 'nh 28 2.5.1.1 Tính tốn m i theo ph ng pháp t+n th t tích lu, 28 2.5.1.2 Tính tố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 v m i c a 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 t i i m nóng 30 2.5.2.2 Xác 'nh t+n th t m i trung bình c a i m nóng tr ng thái bi n ng n h n b(ng ph ng pháp ph+ 31 2.5.2.3 Tu+i th m i trung bình c a i m nóng k t c u jacket 33 2.6 K t lu"n c!a ch ơng 33 CH NG 3: PH NG PHÁP LU-N ÁNH GIÁ M.C SUY GI/M TIN C-Y THEO TH1I 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 t2ng quát ánh giá s suy gi m TC t2ng th c!a 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 v b n c!a k t c u jacket tr ng thái bi n ng n h nc c i 37 3.3.1 -ng su t ng u nhiên k t c u 37 3.3.2 tin c y v b n c a k t c u jacket ng su t có ph+ d i h p 38 3.3.3 tin c y v b n c a k t c u jacket ng su t có ph+ d i r ng 38 -vi - 3.4 Xác nh ) tin c"y theo i u ki n m,i t i m)t i m nóng c!a 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 c a t+n th t m i tr ng thái bi n ng n h n t i m t i m nóng 40 3.4.1.1 Bi u di&n ng su t ng u nhiên t i 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 b(ng k, thu t 3.4.1.3 Xác 'nh k v ng ph m dòng m a 43 ng sai c a t+n th t m i t i i m nóng m t 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 c a t+n th t m i t i i m nóng m t n m……… 48 3.4.2.1 M t xác su t c a t+n th t m i t i i m nóng m t n m 48 3.4.2.2 Hàm phân ph i xác su t c a t+n th t m i t i i m nóng m t n m 48 3.4.2.3 ánh giá tin c y v m i t i i m nóng n m………………… 49 3.5 Xác nh ) tin c"y v m,i t i i m nóng ' th i i m T(n3m) 50 3.6 Xác nh TC c!a k t c u t i 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 c!a k t c u t i i m xét d a i u ki n m,i m' r)ng 51 3.7.1 K v ng c a t/ s t+n th t m i m! r ng 51 3.7.2 Ph 3.7.3 3.8 ng sai c a t/ s t+n th t m i m! r ng 52 tin c y tính theo i u ki n m i m! r ng 52 ánh giá s suy gi m TC theo th i gian c!a KCC jacket d a TC th c t (t2ng th ) c!a KC t i i m xét 52 3.9 ánh giá m c ) suy gi m kh n3ng ch u t i c!a 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 tốn ánh giá an tồn c!a k t c u jacket theo ph ơng pháp truy n th ng theo ph ơng pháp lu"n c!a Lu"n án 55 3.11 K t lu"n c!a ch ơng 60 -vii - CH NG 4: VÍ D4 NG D4NG 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 v cơng trình 61 4.2.2 S li u v 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 qu tính tin c y theo i u ki n b n truy n th ng 65 4.4.1 K t qu tính n i l c ng u nhiên k t c u 65 4.4.2 Ki m tra b n c a ph n t) 67 4.5 ) tin c"y theo i u ki n m,i truy n th ng c!a k t c u jacket 67 4.5.1 Tính tố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 tốn ng su t i m nóng 69 4.5.1.3 Tính tốn t+n th t m i t i i m nóng 71 4.5.2 K t qu tính tốn tin c y theo i u ki n phá h y m i truy n th ng 73 4.6 ánh giá s suy gi m ) tin c"y t2ng th c!a k t c u jacket 73 4.6.1 tin c y t i i m 4.6.2 ánh giá s suy gi m c tr ng c a k t c u b t u khai thác 73 tin c y kh n ng ch'u t i c a 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 qu 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 4.6.2.4 K t qu tính tốn 4.7 ng h p bi n +i khí h u b t th ng: 77 tin c y t+ng th c a cơng trình 79 ánh giá m c ) suy gi m kh n3ng ch u t i c!a i u ki n bi n c c i theo th i gian khai thác cơng trình t i i m xét 80 4.8 K t lu"n c!a ch ơng 81 -viii - PH5N K T LU-N Nh"ng k t qu ã t c 83 Nh"ng óng góp m i c a lu n án 83 Ki n ngh' c a lu n án 84 H ng nghiên c u phát tri n c a lu n án 84 PH4 L4C: Ch ng trình tính m i ng u nhiên RFCAL ………… ……….…….90 -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 -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 -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 -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 -93 - PH4 L4C: CH NG TRÌNH TÍNH M&I “RFCAL” Ph(n Giao di n ch ơng trình RFCAL Ch ng trình RFCAL c l p b(ng ngôn ng" Visual Basic th t m i cho tr ng thái bi n n m cho nhi u n m Ch t ng t o s li u k t n i v i SACS 5.3 M t s giao di n c a ch ng trình c3ng tính tốn n i l c ng trình nh sau: Hình Giao di n c a ch ng trình, tu n t b c tính tốn Hình tính tốn t+ n c, l c tách tính tốn ng su t c a ph n t) Hình c s li u ph+ sóng t o file d" li u u vào cho SACS 5.3 Hình Tính tốn t+n th t m i cho 1n m -94 - Ph(n Mã ngu n c!a 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 -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 -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 FileNumberMXA, FileNumberMYA, -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 -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 -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 -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) ' - -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 -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, ",") -103 - 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