Microsoft Word Luan an TS Dao van Dinh 5 9 2014 doc BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC GIAO THÔNG VẬN TẢI ĐÀO VĂN DINH DỰ BÁO TUỔI THỌ SỬ DỤNG CỦA CẦU BÊ TÔNG CỐT THÉP VEN BIỂN VIỆT NAM DO XÂM NHẬP[.]
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC GIAO THÔNG VẬN TẢI ĐÀO VĂN DINH DỰ BÁO TUỔI THỌ SỬ DỤNG CỦA CẦU BÊ TÔNG CỐT THÉP VEN BIỂN VIỆT NAM DO XÂM NHẬP CLO CHUY£N NGµNH : Kü THUậT XÂY DựNG CÔNG TRìNH ĐặC BIệT Mà Số : 62.58.02.06 LN ¸N TIÕN SÜ Kü THT Ngêi Híng dÉn khoa häc: GS.TS PHẠM DUY HỮU PGS.TS BÙI TRỌNG CẦU HÀ NỘI - 2014 i L I CAM ðOAN Tơi xin cam đoan lu n án cơng trình nghiên c u c a cá nhân tơi Các k t qu nêu lu n văn trung th c chưa đư c cơng b cơng trình khác HÀ N I NGÀY THÁNG NĂM 2014 TÁC GI ðÀO VĂN DINH ii L I C M ƠN Lu n án ñư c th c hi n dư i s hư ng d n tr c ti p c a GS.TS Ph m Duy H u PGS.TS Bùi Tr ng C u Tôi xin chân thành c m ơn th y hư ng d n ch d n t n tình đóng góp ý ki n q báu ñ giúp th c hi n lu n án Tôi xin c m ơn quý th y b mơn Cơng Trình Giao Thơng Thành Ph Cơng trình Th y, đ c bi t GS.TS Nguy n Vi t Trung giúp đ tơi r t nhi u trình h c t p nghiên c u Tôi xin trân tr ng c m ơn PGS.TS Tr n Quang Vinh đóng góp ý ki n cho lu n án Tôi xin trân tr ng c m ơn Phòng ðào t o sau ð i h c trư ng ð i h c Giao Thơng V n t i t o ñi u ki n thu n l i cho trình h c t p nghiên c u Tôi xin trân tr ng c m ơn Ban Giám Hi u trư ng ð i h c Giao Thơng V n t i, lãnh đ o khoa Cơng Trình t o u ki n đ tơi ñư c h c t p nghiên c u Cu i bày t" c m ơn ñ#ng nghi p, gia ñình ngư i thân ñã giúp đ tơi q trình h c t p, nghiên c u HÀ N I NGÀY THÁNG NĂM 2014 TÁC GI ðÀO VĂN DINH iii M CL C L i cam ñoan i L i c m ơn ii M$c l$c iii Danh m$c b ng bi u vi Danh m$c hình v% vii Danh m$c ch vi t t&t, ký hi u ix M' ð(U CHƯƠNG T+NG QUAN CÁC NGHIÊN C,U V- TU+I TH S/ D0NG C1A C(U BÊ TÔNG C2T THÉP DO XÂM NH3P CLO TRÊN TH4 GI5I VÀ VI6T NAM 1.1 M7 ñ u 1.2 Các nghiên c u v ñ b n tu8i th s9 d$ng th gi i 1.2.1 Các nghiên c u v lý thuy t ñ b n c a bê tông 1.2.2 Các nghiên c u v ch xâm nh p clo vào bê tông 10 1.2.3 Các nghiên c u ch ăn mòn c t thép 13 1.2.4 Các thí nghi m v s c kháng xâm nh p clo c a bê tông 21 1.2.5 Các nghiên c u v h s khu ch tán 30 1.2.6 Các nghiên c u v th i gian kh7i ñ u ăn mòn th i gian lan truy n ăn mòn, tu8i th s9 d$ng 32 1.3 Các nghiên c u nư c 35 1.4 Nh n xét hư ng nghiên c u c a lu n án 37 1.5 M$c tiêu c a lu n án 38 1.6 N i dung phương pháp nghiên c u 38 1.6.1 N i dung nghiên c u 38 1.6.2 Phương pháp nghiên c u 38 CHƯƠNG NGHIÊN C,U XÁC ð:NH H6 S2 KHU4CH TÁN CLO TRONG BÊ TÔNG 39 2.1 Gi i thi u chung 39 2.2 Thí nghi m xác ñ;nh h s khu ch tán clo bê tơng 40 2.2.1 Thí nghi m th m nhanh clo ASTM C1202 40 2.2.2 K t qu thí nghi m th m nhanh clo ASTM C1202 48 iv 2.2.3 Xây d ng cơng th c xác đ;nh h s khu ch tán clo t< k t qu thí nghi m C1202 49 2.2.4 Áp d$ng phương trình DC1202 tính h s khu ch tán clo t< k t qu thí nghi m C1202 50 2.3 T8ng k t k t qu d báo h s khu ch tán th gi i 51 2.3.1 D báo h s khu ch tán (D28) t< t= l nư c xi măng 52 2.3.2 nh hư7ng c a nhi t ñ ñ n h s khu ch tán 52 2.3.3 nh hư7ng c a th i gian ñ n h s khu ch tán 53 2.3.4 nh hư7ng c a mu i silic ñ n h s khu ch tán 55 2.3.5 nh hư7ng c a ñ >m ñ n h s khu ch tán 57 2.3.6 nh hư7ng c a n t ñ n h s khu ch tán 57 2.3.7 H s khu ch tán bi u ki n 58 2.3.8 Tính h s khu ch tán theo d báo theo công th c kinh nghi m 59 2.4 So sánh k t qu tính tốn th o lu n 60 2.4.1 So sánh k t qu tính h s khu ch tán D 60 2.4.2 Phân tích, nh n xét k t qu tính h s khu ch tán D 63 2.5 K t lu n chương 63 CHƯƠNG XÂY D?NG MƠ HÌNH D? BÁO TU+I TH S/ D0NG C(U BÊ TÔNG C2T THÉP VEN BI@N VI6T NAM DO XÂM NH3P CLO 65 3.1 Gi i thi u chung 65 3.1.1 Khái ni m tu8i th s9 d$ng 66 3.1.2 Tu8i th s9 d$ng theo tác ñ ng c a s xâm nh p clo môi trư ng bi n 67 3.2 Xây d ng mơ hình d báo th i gian kh7i ñ u ăn mòn 71 3.2.1 T8ng quát 71 3.2.2 Các tham s c a mơ hình 71 3.2.3 Xây d ng mơ hình d báo th i gian kh7i đ u ăn mịn 80 3.2.4 Xây d ng chương trình tính th i gian kh7i đ u ăn mòn bAng Mathlab 88 3.3 Xây d ng mơ hình d báo th i gian lan truy n ăn mòn 88 3.3.1 Các v n ñ chung 88 3.3.2 Xem xét mơ hình hi n có 89 v 3.3.3 Mơ hình đ xu t 94 3.4 Sơ đ# thu t tốn tính tu8i th s9 d$ng c a c u bê tông “LifeConBridge” 109 3.5 K t qu tính nh n xét 111 3.5.1 Ki m ch ng k t qu tính c a mơ hình đ xu t 111 3.5.2 K t qu tính tốn thí d$ 112 3.5.3 Nh n xét k t qu : 114 3.6 D báo tu8i th s9 d$ng c a c u bê tông c t thép ven bi n Vi t Nam 115 3.6.1 ð c ñi m khí h u vùng ven bi n Vi t Nam bi n đ8i khí h u 115 3.6.2 D báo tu8i th s9 d$ng c a c u bê tông c t thép ven bi n Vi t Nam 118 3.7 K t lu n chương 120 CHƯƠNG CÁC BI6N PHÁP KÉO DÀI TU+I TH S/ D0NG VÀ THÍ D0 TÍNH TỐN .122 4.1 Các bi n pháp dài tu8i th s9 d$ng 122 4.1.1 V i k t c u m i 122 4.1.2 V i k t c u cũ 127 4.2 Thí d$ tính tốn cho m t s b ph n c u T 127 4.2.1 Các thơng sơ tính tu8i th s9 dung c a c u T 127 4.2.2 K t qu tính tốn tu8i th s9 d$ng 129 4.2.3 Các k t lu n rút t< thí d$ nghiên c u 130 4.3 K t lu n chương 131 K4T LU3N VÀ KI4N NGH: 132 1/ K t lu n 132 2/ Ki n ngh; 134 3/ Hư ng nghiên c u ti p theo 135 DANH M0C CƠNG TRÌNH C1A TÁC GI .136 TÀI LI6U THAM KH O 137 PH0 L0C 145 PH0 L0C 155 vi Danh m c b ng bi u B B B B B B B B B B B B B B B B B B B B B B B B ng 1.1 nh hư7ng c a tham s khác ñ i v i q trình ăn mịn 21 ng 1.2 T8ng k t phương pháp thí nghi m xâm nh p clo 30 ng 2.1 Thành ph n bê tông c a t8 m u 42 ng 2.2 M c ñ th m ion clo 45 ng 2.3 K t qu thí nghi m theo ASTM C1202 48 ng 2.4 K t qu tính DC1202 t< thí nghi m 51 ng 2.5 Các giá tr; m cho lo i bê tông 54 ng 2.6 K t qu D theo d báo công th c kinh nghi m 60 ng 2.7 DC1202, D d báo công th c kinh nghi m 61 ng 2.8 So sánh k t qu tính h s D 62 ng 3.1 T c ñ tích lũy n#ng ñ l n nh t c a clo b m t 73 ng 3.2 Tóm t&t tiêu chu>n xác đ;nh gi i h n t i ña cho phép c a clo v a bê tông d ng l c 78 ng 3.3 nh hư7ng c a canxi nitrit (CNI) ngư ng clo t i h n 79 ng 3.4 K t qu tinh th i gian kh7i ñ u ăn mòn H=100% 111 ng 3.5 K t qu tinh th i gian kh7i đ u ăn mịn H=75% 112 ng 3.6 K t qu tinh th i gian lan truy n ăn mòn theo n t bê tông b o v 112 ng 3.7 K t qu tinh th i gian kh7i ñ u ăn mòn theo tham s 113 ng 3.8 K t qu tính th i gian lan truy n ăn mịn theo tham s (quan m 1: n t hồn tồn bê tơng b o v ) 113 ng 3.9 K t qu tính th i gian lan truy n ăn mịn theo tham s (quan m 2: ăn mòn gây nguy hi m) 114 ng 4.1 K t qu tính tu8i th s9 d$ng, v i gi i pháp k t h p 126 ng 4.2 Thông s k t c u v t li u 128 ng 4.3 Thông s môi trư ng 129 ng 4.4 K t qu tính tu8i th s9 d$ng, toán 1D 129 ng 4.5 K t qu tính tu8i th s9 d$ng, toán 2D 130 B ng 4.6 K t qu tính th i gian lan truy n ăn mòn t2 130 vii Danh m c hình v Hình 1.1: S lư ng hư h i quan sát ñư c Nh t b n (ngu#n d li u Prof Hiroshi Mutsuyoshi 2001) Hình 1.2: Bi u đ# Pourbaix quan h th ñi n c c ñ pH c a h Fe-H2O 15 Hình 1.3 Các ph n ng c c dương c c âm (Beeby) 18 Hình 1.4 Th tích tương đ i c a s n ph>m ăn mòn s&t 20 Hình 1.5 Bi u đ# th hi n hư h i ăn mòn gây n t, v , tách l p 20 Hình 1.6 Sơ đ# thí nghi m theo AASHTO T259 (salt ponding) 21 Hình 1.7 Sơ đ# thí nghi m khu ch tán kh i- Bulk Diffusion Test (NordTest NTBuild 443) 22 Hình 1.8 Sơ đ# thí nghi m AASHTO T277 (ASTM C1202) 23 Hình 1.9 Sơ đ# thí nghi m kB thu t n di 25 Hình 1.10 Sơ đ# thí nghi m n di c a Tang Nilsson 26 Hình 1.11 Sơ đ# thí nghi m n di nhanh (NordTest NTBuild 492) 27 Hình 1.12 Thi t b; đo ñi n tr7 su t m t chi u 29 Hình 1.13 Sơ đ# phương pháp ñi m ño c a Wenner 29 Hình 1.14: Tu8i th s9 d$ng c a k t c u bê tông c t thép: Mơ hình hai giai đo n c a Tuuti 1980 33 Hình 2.1 Sơ đ# bơm hút chân khơng m u th9 C1202 43 Hình 2.2 Sơ đ# đo n tích 44 Hình 2.3 Các nh thí nghi m th m nhanh clo theo ASTM C1202 47 Hình 2.4 nh hư7ng c a t= l W/C h s khu ch tán Clo bê tơng nhi t đ 20oC 52 Hình 2.5 nh hư7ng c a tro bay x ñ i v i h s khu ch tán 55 Hình 2.6 nh hư7ng c a silica fum đ i v i h s khu ch tán Hình 2.7 Bi u đ# k t qu tính D 62 Hình 3.1: ð;nh nghĩa tu8i th s9 d$ng kéo dài tu8i th s9 d$ng 67 Hình 3.2: nh hư7ng c a màng sơn ph b m t 75 Hình 3.3: góc ph n tư c t 2D 84 Hình 3.4: Các bi n góc ph n tư c t 2D 85 Hình 3.5: Áp l c bê tơng hình thành s n ph>m ăn mịn (Mơ hình c a Liu) 91 viii Hình 3.6: Lý tư7ng hóa c a bê tơng b o v m t hình tr$ thành dày: 95 Hình 3.7: Kho ng th i gian t< kh7i đ u ăn mịn thép đ n n t hồn tồn bê tơng b o v t i nguy hi m ch;u l c 96 Hình 3.8: Sơ đ# c lư ng cho m t mát bán kính thép ∆rs2 102 Hình 3.9: M i quan h gi a t c đ ăn mịn đ >m tương đ i v i bê tơng tu8i năm có hàm lư ng ion Cl- 1.8kg/m3 nhi t đ 23oC 105 Hình 3.10 M t ñ dòng ăn mòn v i th i gian Ccl=1.25kg/m3, nhi t ñ 20oC, ñ >m H=75% theo Liu Weyers 107 Hình 3.11 Sơ đ# thu t tốn tính tu8i th s9 d$ng xâm nh p clo 110 Hình 3.12: M c tăng nhi t đ trung bình năm (oC) vào gi a th k= 21 theo k;ch b n phát th i trung bình 116 Hình 3.13: M c tăng nhi t đ trung bình năm (oC) vào cu i th k= 21 theo k;ch b n phát th i trung bình 117 Hình 3.14: Các vùng môi trư ng bi n c a tr$ c u bê tơng 119 Hình 3.15: ð;nh tính v phân b n#ng đ clo b m t 120 Hình 4.1: Quan h gi a chi u dày l p bê tông b o v v i th i gian kh7i ñ u ăn mịn w/c=0.35, nhi t đ 20oC, đ >m H=80% cho vùng khí quy n bi n 123 Hình 4.2: Quan h gi a chi u dày l p bê tông b o v v i th i gian t< kh7i đ u ăn mịn đ n n t cho bê tơng f’c=40MPa,đư ng kính c t thép d=16mm, nhi t đ 20oC, đ >m H=80% cho vùng khí quy n bi n 123 Hình 4.3: Quan h gi a t= l nư c xi măng w/c v i th i gian kh7i đ u ăn mịn (chi u dày l p bê tông b o v L=70mm, nhi t đ 20oC, đ >m H=75% cho vùng khí quy n bi n) 124 Hình 4.4: Quan h gi a t= l mu i si líc (silica fume) v i th i gian kh7i ñ u ăn mòn (chi u dày l p bê tơng b o v L=75mm, nhi t đ 25oC, đ >m H=75% cho vùng khí quy n bi n) 125 ix Danh m c ch vi t t t, ký hi u ∆rs1 Lư ng gi m bán kính c a c t thép t i th i ñi m b&t ñ u n t ∆rs Lư ng gi m bán kính c a c t thép g chèn vào v t n t Ecef Mơ đun đàn h#i có hi u (xét ñ n t< bi n); ψ H s t< bi n c a bê tông; δ0 Chi u dày vùng x p bao quanh c t thép; νc H s poisson c a bê tông νc=0,18-0,20; δcon Chuy n v; hư ng tâm c a bê tông; ∆d S thay đ8i đư ng kính c a c t thép νr H s poisson c a g s&t; δrust Chuy n v; nén hư ng tâm c a g (các s n ph>m ăn mòn) ∆T m c tăng nhi t đ thí nghi m C1202 “LifeConBridge”- Service Life of concrete Bridge- Tu8i th s9 d$ng c a c u bê tơng A Di n tích ti t di n m u th9 BTCT Bê tông c t thép C N#ng ñ clo C(x,t) N#ng ñ clo t i chi u sâu x th i gian t C0 N#ng đ clo ban đ u bê tơng CNI Canxi nitrit Cs N#ng ñ clo t i b m t bê tơng Cs(t) N#ng đ clo t i b m t bê tơng t i th i m t d ðư ng kính c t thép D H s khu ch tán clo bê tông D(t) H s khu ch tán clo bê tông th i ñi m t D(T) H s khu ch tán clo bê tơng nhi t đ T D28 H s khu ch tán clo bê tông tu8i 28 ngày DPC H s khu ch tán clo bê tông thư ng DSF H s khu ch tán clo bê tơng có mu i silic (silica fume) E ði n th áp d$ng 141 40 Mark G Richardson (2004), Fundamentals of durable reinforced concrete, published in the Taylor & Francis e-Library, simultaneously published in the USA and Canada 41 Morinagn, S (1989) “Prediction of Service Lives of Reinforced Concrete Buildings Based on Rate of Corrosion of Reinforcing Steel,S” pecial Report of the Institute of Technology, Skimiza Corporation, Japan 42 Morris, W., Moreno, E.I and Sagues, A A.(1996), “Practical Evaluation of Resistivity of Concrete in TestCylinders using a Wenner Array Probe”, Cement and Concrete Research, Vol 26, No 12, pp.1779-1787 43 M.S Shetty (2005) “ Concrete Technology: Theory and Practice” published by S Chand & Company Pvt.Ltd 7361, Ram Naga, New Delhi 44 Mutsuyoshi, H (2001): ‘Present Situation of durability of posttensioned pc bridges in Japan’, in: Durability of post-tensioning tendons (Taerwe, L ed.), Fédération Internationale du Béton, Lausanne, pp 75-88 45 NCHRP- NATIONAL COOPERATIVE HIGHWAY RESEARCH PROGRAM, REPORT 558-Manual on Service Life of CorrosionDamaged Reinforced Concrete Bridge Superstructure Elements 46 Nielsen A (1985) Durability, pp 200-243 in Beton Bogen, Aalborg Cement Company, Aalborg, Portland 47 Omar S Baghabra Al-Amoudi, Walid A Al-Kutti, Shamsad Ahmad, Mohammed Maslehuddin (2009) “Correlation between compressive strength and certain durability indices of plain and blended cement concretes” Cement & Concrete Composites 31 672–676 48 P Vassie (1984), “Reinforcement corrosion and the durability of concrete bridges”, Proceeding of Institution of Civil Engineers 76, 713– 723 142 49 Purvis RL, Babaei K, Clear KC, Markow MJ (1994), Life-cycle cost analysis for protection and rehabilitation of concrete bridges relative to reinforcement corrosion, Washington, DC: Strategic Highway Research Program, National Research Council 50 Rodriguez, O.G & Hooton, R.D (March 2003), “Influence of cracks on chlo-ride ingress into concrete”, ACI Materials Journal, Vol 100, No 2, p 120-126 51 Schießl P.; Wiens U (1995), “Rapid Determination of Chloride Diffusivity in Concrete with Blending Agents”, Proceedings of RILEM International Workshop on Chloride Penetration into Concrete, St.Rémy-lès-Chevreuse, France, pp 115-125 52 Sharif, A., Loughlin, K F., Azad, A K and Nawaz, C M.,(1999), Determination of the effective diffusion coefficient in concrete via a gas diffusion technique, Proceedings of the international conference on concrete durability and repair technology, Edited by Dhir R K and McCarthy M J., published by Thomas Telford 53 Stanish, K (2000), Predicting the Diffusion Coefficient of Concrete from Mix Parameters, University of Toronto Report 54 Swamy,R.N., Hamada and Laiw,J.C, (Jul 1994) “ A citical evaluation of chloride penetration into concrete in marine environment”, in “ Corrosion and Corrosion Protection of steel in concrete” Proceedings of an International Conference, University of Sheffield, England, 404-419 55 Tang, L and Nilsson, L-O (1992), “Chloride diffusivity in high strength concrete at different ages.” Nordic Concrete Research, pp 162-171 56 Tang Luping, Joost Gulikers (2007) “On the mathematics of timedependent apparent chloride diffusion coefficient in concrete” - Cement and Concrete Research 37, 589–595 57 Technical Research Centre of Finland (VTT) (2003), LIFECON Deliverable D 3.2, Service Life ModelsLife -Cycle Management of 143 Concrete Infrastructures for Improved Sustainability, Project funded by the European Community under the Competitive and Sustainable Growth Programme (2001-2003) 58 TRB (1996), Transverse cracking in newly constructed bridge decks, National Co-operative Highway Research Program Report 380, Transportation Research Board, National Academy Press, Washington, 126 p 59 Tuutti, K (1980) “Service life of structures with regard to corrosion of embedded steel”, Proceedings of the International Conference on Performance of Concrete in Marine Environment, ACI SP-65, pp 223236 60.57 T Liu, R.W Weyers (1998), “Modeling the dynamic corrosion process in chloride contaminated concrete structures”, Cement and Concrete Research 28 (3) 365–379 61 Uhlig, H H (1971), Corrosion and Corrosion Control, Wiley, New York 62 Vu, K A T and Stewart, M G (2000): “Structural reliability of concrete bridges including improved chloride-induced corrosion models”, Structural Safety, Vol 22, pp 313-333 63 Weyers, R.E (1998), “Service life model for concrete structures in chloride laden environments.” ACI Materials Journal, Vol 95 (4), pp 445-453 64 Whiting, D (1981), “Rapid Measurement of the Chloride Permeability of Concrete”, Public Roads, Vol 45,No 3, pp 101-112 65 W López, J.A González (1993), “Influence of the degree of pore saturation on the resistivity of concrete and the corrosion rate of steel reinforcement”, Cement and Concrete Research 23 (2) 368–376 66 Youping Liu (October 21, 1996), Modeling the Time-to-Corrosion Cracking of the Cover Concrete in Chloride Contaminated Reinforced 144 Concrete Structures, Doctor of Philosophy, Virginia Polytechnic Institute and State University 67 Zhang, and Gjorv, O.E (October 1995), “Effect of Ionic Interaction in Migration Testing of Chloride Diffusivity in Concrete”, Cement and Concrete Research, V.25, No.7, pp.1535-1542 145 PH# L#C Ph n m m LifeConBridge BÀI TOÁN 1D VÙNG TH1Y TRI-U % Du bao tuoi tho su dung cau be tong cot thep ven bien Viet Nam xam nhap clo-"Sevice Life Concrete Bridges" % Phuong phap sai phan huu han Crank-Nicholson clear all; % Cac thong so dau vao co ban cua bai toan fc = 50; %Cuong chiu nen cua be tong o 28 MPa fsp = 0.59*sqrt(fc); %Cuong chiu keo uon MPa gammac = 2500; % Trong luong rieng cua be tong kg/m3 Ec = 0.043*(gammac)^(1.5)*sqrt(fc);%Mo-duyn dan hoi cua be tong phi = 2;% He so tu bien cua be tong Ecef =Ec/(1+phi); roth =2.5; % Phan tram mat mat dien tich TD thep toi han db = 16; %Duong kinh cot thep vc = 0.18; % He so poisson n = 2.7; %He so gian no the tich n=2.7-3.0 k1 = 0.5; % muc lap day vao vet nut delta0 = 15; % chieu day lop xi micro - m r0 = db/2+delta0/1000; b = 1000; Lc = 75; %Chieu day cau kien -mm % Chieu day lop be tong bao ve-mm dt = 0.1; % buoc thoi gian tinh bang nam dx =5; %Do lon cua luoi chia theo chieu day - mm nx1 =Lc/dx; % So buoc chia cua khoang cach nx = b/dx; NX = 0.35; %Ty le nuoc/xi mang FA = 0; % phu gia Tro bay SF= 2;% phu gia silica-fume SG=0; %Phu gia xi lo Tref = 293.15; T = 299.05; % Nhiet tuyet doi K D28 = exp(-0.165*SF)*10^(-12.06+2.4*NX); % He so khuech tan o 28 tuoi m = 0.2+0.4*(FA/50+SG/70); % D28 thi nghiem Hc = 0.75; %Hc = 75% am chuan H = 0.80; %Do dam tuong doi cua mtH U = 35000; G = 8.314; fT = exp(U/G*(1/Tref-1/T)); fH = 1/(1+((1-H)/(1-Hc))^4); 146 ucth = 0.15; uth = 0.15*ones(nx1+1,1); % Nong clorua tai thoi diem t = for i = 1:nx+1 x(i) =(i-1)*dx; u(i,1) =0; end k=1; u_new =u; while u_new(nx1+1)