Luận án Tiến sĩ Nghiên cứu giải pháp giảm ứng suất nhiệt của bê tông đầm lăn trong xây dựng công trình thủy lợi thủy điện tại Việt Nam

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT VIỆN KHOA HỌC THỦY LỢI VIỆT NAM NGUYỄN MINH VIỆT NGHIÊN CỨU GIẢI PHÁP GIẢM ỨNG SUẤT NHIỆT CỦA BÊ TƠNG ĐẦM LĂN TRONG XÂY DỰNG CƠNG TRÌNH THỦY LỢI THỦY ĐIỆN TẠI VIỆT NAM (Chỉnh sửa sau phản biện kín) LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT VIỆN KHOA HỌC THỦY LỢI VIỆT NAM NGUYỄN MINH VIỆT NGHIÊN CỨU GIẢI PHÁP GIẢM ỨNG SUẤT NHIỆT CỦA BÊ TÔNG ĐẦM LĂN TRONG XÂY DỰNG CÔNG TRÌNH THỦY LỢI THỦY ĐIỆN TẠI VIỆT NAM Chuyên ngành: Kỹ thuật xây dựng cơng trình thủy Mã số: 62.58.02.02 NGƯỜI HƯỚNG DẪN KHOA HỌC PGS.TS Hồng Phó Un GS.TS Phạm Ngọc Khánh HÀ NỘI, NĂM 2017 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 văn trung thực, không chép từ nguồn hình thức Việc tham khảo nguồn tài liệu (nếu có) thực trích dẫn ghi nguồn tài liệu tham khảo quy định Tác giả luận án Nguyễn Minh Việt i LỜI CÁM ƠN Sau thời gian thực hiện, với nỗ lực thân với giúp đỡ tận tình Thầy bạn bè đồng nghiệp, Luận án tiến sĩ: “Nghiên cứu giải pháp giảm ứng suất nhiệt bê tông đầm lăn xây dựng cơng trình thủy lợi thủy điện Việt Nam” hoàn thành Tác giả xin bày tỏ biết ơn sâu sắc đến Ban Giám đốc, Trung tâm đào tạo Hợp tác quốc tế, Ban TCHC Viện Khoa học thủy lợi Việt Nam giúp đỡ tạo điều kiện tốt cho NCS thời gian thực Luận án Xin đặc biệt cám ơn hướng dẫn, giúp đỡ tận tình PGS.TS Hồng Phó Uyên, GS.TS Phạm Ngọc Khánh Các Thầy tạo điều kiện tốt cho NCS trình học tập hoàn thành Luận án Tác giả chân thành cám ơn đồng nghiệp bạn bè nhiệt tình giúp đỡ, tạo điều kiện thuận lợi cho tác giả trình học tập thực Luận án Do lực thân nhiều hạn chế, chắn Luận án khơng tránh khỏi thiếu sót Tác giả kính mong Thầy Cơ bảo, đồng nghiệp đóng góp ý kiến để tác giả hoàn thiện, tiếp tục nghiên cứu phát triển đề tài Hà Nội, ngày tháng 06 năm 2017 Tác giả luận án Nguyễn Minh Việt ii MỤC LỤC DANH MỤC CÁC HÌNH VẼ viii DANH MỤC BẢNG BIỂU xii DANH MỤC CÁC TỪ VIẾT TẮT xiv MỞ ĐẦU 1 Tính cấp thiết đề tài luận án Mục đích nghiên cứu .2 Đối tượng nghiên cứu Phạm vi nghiên cứu Phương pháp nghiên cứu 5.1 Phương pháp lý thuyết 5.2 Phương pháp mơ hình tốn .2 Nội dung nghiên cứu .2 Ý nghĩa khoa học thực tiễn đề tài luận án 7.1 Ý nghĩa khoa học 7.2 Ý nghĩa thực tiễn Cấu trúc luận án CHƯƠNG TỔNG QUAN VỀ BÊ TÔNG ĐẦM LĂN VÀ NHỮNG VẤN ĐỀ CẦN NGHIÊN CỨU 1.1 Bê tông đầm lăn 1.2 Tính học BTĐL 1.2.1 Cường độ kháng nén BTĐL 1.2.2 Cường độ kháng kéo BTĐL 1.2.3 Mô đun đàn hồi BTĐL 1.2.4 Biến dạng BTĐL .6 1.2.4.1 Từ biến BTĐL 1.2.4.2 Co ngót BTĐL 1.2.4.3 Giá trị kéo dãn giới hạn BTĐL 1.2.5 Tính học BTĐL theo thời gian 1.3 Tình hình xây dựng đập BTĐL giới Việt Nam .9 1.3.1 Tình hình xây dựng đập BTĐL giới 1.3.2 Tình hình xây dựng đập BTĐL Việt Nam .11 iii 1.4 Tình hình nghiên cứu nhiệt ứng suất nhiệt đập BTĐL 13 1.4.1 Tình hình nghiên cứu nhiệt ứng suất nhiệt BTĐL giới .13 1.4.2 Tình hình nghiên cứu nhiệt ứng suất nhiệt BTĐL Việt Nam 13 1.5 Vấn đề nứt nhiệt đập BTĐL 16 1.6 Vấn đề cần nghiên cứu đặt luận án 17 1.7 Kết luận chương 17 CHƯƠNG NHIỆT VÀ KHỐNG CHẾ ỨNG SUẤT NHIỆT TRONG ĐẬP BÊ TÔNG ĐẦM LĂN 18 2.1 Đặt vấn đề 18 2.2 Nguồn phát sinh nhiệt BTĐL 18 2.3 Vấn đề trao đổi nhiệt BTĐL 19 2.3.1 Trao đổi nhiệt dẫn nhiệt 19 2.3.2 Trao đổi nhiệt đối lưu 21 2.3.3 Trao đổi nhiệt xạ 22 2.4 Cơ chế nứt nhiệt bê tông khối lớn 23 2.4.1 Nứt bề mặt 23 2.4.2 Nứt xuyên 24 2.5 Yêu cầu khống chế nhiệt cho đập BTĐL 24 2.5.1 Chênh lệch nhiệt độ đáy đập 25 2.5.2 Chênh lệch nhiệt độ lớp 26 2.5.3 Chênh lệch nhiệt độ .28 2.5.3.1 Ảnh hưởng chênh lệch nhiệt độ .28 2.5.3.2 Thực tế khống chế chênh lệch nhiệt độ 29 2.6 Phương pháp giải toán nhiệt 30 2.6.1 Các phương pháp giải toán nhiệt 30 2.6.1.1 Phương pháp giải tích 30 2.6.1.2 Phương pháp toán tử 31 2.6.1.3 Phương pháp gần 31 2.6.1.4 Lựa chọn phương pháp giải toán nhiệt 33 2.6.2 Cơ sở tính tốn nhiệt ứng suất nhiệt theo phương pháp phần tử hữu hạn .33 2.6.2.1 Các giả thiết 33 iv 2.6.2.2 Xác định trường nhiệt độ 34 2.6.2.3 Xác định trường ứng suất .36 2.7 Kết luận chương 42 CHƯƠNG TÍNH TỐN NHIỆT VÀ ỨNG SUẤT NHIỆT TRONG ĐẬP BTĐL THEO PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN VÀ ỨNG DỤNG TRONG PHẦN MỀM ANSYS 43 3.1 Đặt vấn đề 43 3.2 Tính tốn nhiệt thủy hóa vật liệu chất kết dính BTĐL 43 3.3 Tính nhiệt ứng suất nhiệt đập BTĐL q trình thi cơng ngơn ngữ lập trình tham số (APDL) ANSYS 46 3.3.1 Giới thiệu phần mềm ANSYS 46 3.3.2 Giải toán phần mềm ANSYS 47 3.3.3 Cơng phân tích nhiệt phần mềm ANSYS 48 3.3.4 Cơ sở xây dựng tốn tính nhiệt ứng suất nhiệt đập BTĐL .49 3.3.5 Sơ đồ khối tính nhiệt ứng suất nhiệt đập BTTL .50 3.3.6 Xây dựng tốn tính nhiệt ứng suất nhiệt đập BTTL 50 3.3.6.1 Mô tả kết cấu đập 50 3.3.6.2 Mơ hình hóa kết cấu đập BTĐL 51 3.3.6.3 Xây dựng tốn tính nhiệt ứng suất nhiệt đập trọng lực BTĐL 51 3.4 Kiểm nghiệm tính tốn nhiệt cho đập BTĐL Sơn La 52 3.4.1 Giới thiệu cơng trình 52 3.4.2 Các tiêu lý nhiệt bê tông đá .53 3.4.3 Các điều kiện biên nhiệt 54 3.4.4 Tiến độ thi công 55 3.4.5 Kết quan trắc nhiệt thân đập BTĐL Sơn La .57 3.4.6 Mặt cắt kiểm tra 57 3.4.7 Mơ hình tính tốn 58 3.4.8 Kết tính toán 60 3.4.9 Nhận xét .62 3.5 Kết luận Chương .63 v CHƯƠNG BIỆN PHÁP GIẢM NHIỆT VÀ ỨNG SUẤT NHIỆT ĐẬP BÊ TÔNG ĐẦM LĂN PHÙ HỢP VỚI ĐIỀU KIỆN VIỆT NAM 64 4.1 Đặt vấn đề 64 4.2 Cơ sở phân vùng nghiên cứu 65 4.2.1 Điều kiện tự nhiên .65 4.2.2 Điều kiện nhiệt độ độ ẩm môi trường .65 4.2.3 Điều kiện nguồn cung ứng vật liệu PGK 68 4.2.4 Phân vùng nghiên cứu 69 4.3 Các tham số dùng cho nghiên cứu 69 4.3.1 Các tham số đầu vào cố định .69 4.3.2 Tham số khống chế nứt .71 4.4 Nghiên cứu ảnh hưởng nhiệt độ khơng khí đến ứng suất nhiệt đập BTĐL 71 4.4.1 Cơ sở nghiên cứu 71 4.4.2 Khi nhiệt độ khơng khí xem nhiệt độ trung bình năm 71 4.4.3 Khi nhiệt độ khơng khí thay đổi theo thời gian (ngày) .72 4.4.4 Nhận xét .73 4.5 Nghiên cứu ảnh hưởng độ ẩm khơng khí đến ứng suất nhiệt đập BTĐL 73 4.5.1 Cơ sở nghiên cứu 73 4.5.2 Kết tính tốn 74 4.5.3 Nhận xét .75 4.6 Nghiên cứu ảnh hưởng thành phần khoáng vật liệu XM đến ứng suất nhiệt đập BTĐL 75 4.6.1 Cơ sở nghiên cứu 75 4.6.2 Kết tính tốn 76 4.6.3 Nhận xét .81 4.7 Nghiên cứu ảnh hưởng nhiệt độ đổ bê tông đến ứng suất nhiệt đập BTĐL .82 4.7.1 Cơ sở nghiên cứu 82 4.7.2 Kết tính tốn 83 4.7.3 Nhận xét .87 4.8 Nghiên cứu ảnh hưởng hàm lượng PGK đến ứng suất nhiệt đập BTĐL .88 4.8.1 Nghiên cứu ảnh hưởng hàm lượng PGK đến nhiệt ứng suất nhiệt thân đập BTĐL 88 vi 4.8.2 Kiến nghị hàm lượng PGK tổng lượng chất kết dính cho khu vực .90 4.8.3 Nhận xét .92 4.9 Đề xuất giải pháp giảm ứng suất nhiệt BTĐL hợp lý cho khu vực 93 4.9.1 Cơ sở đề xuất giải pháp hợp lý 93 4.9.2 Đề xuất giải pháp giảm ứng suất nhiệt hợp lý cho khu vực 94 4.9.2.1 Đối với khu vực miền núi phía Bắc 94 4.9.2.2 Đối với khu vực Bắc Trung .95 4.9.2.3 Đối với khu vực Nam Trung Tây Nguyên 96 4.10 Kiểm nghiệm tính tốn ứng suất nhiệt đập BTĐL Trung Sơn – Thanh Hóa 98 4.10.1 Giới thiệu cơng trình 98 4.10.1.1 Giới thiệu cơng trình 98 4.10.1.2 Các tiêu lý nhiệt bê tông đá 99 4.10.1.3 Các điều kiện biên nhiệt 100 4.10.1.4 Mặt cắt kiểm tra 101 4.10.2 Kết quan trắc nhiệt 102 4.10.3 Tính tốn nhiệt ứng suất nhiệt đập .104 4.10.3.1 Mơ hình tính tốn .104 4.10.3.2 Kết tính toán 105 4.10.3.3 Nhận xét 106 4.10.4 Giải pháp giảm ứng suất nhiệt đập BTĐL Trung Sơn 107 4.11 Kết luận chương 108 KẾT LUẬN VÀ KIẾN NGHỊ 111 Kết luận 111 Những đóng góp Luận án 112 Kiến nghị .112 DANH MỤC CƠNG TRÌNH ĐÃ CÔNG BỐ 113 TÀI LIỆU THAM KHẢO 114 PHỤ LỤC I: BẢNG TỔNG HỢP KẾT QUẢ THÍ NGHIỆM THÀNH PHẦN KHỐNG CỦA MỘT SỐ LOẠI XI MĂNG ĐIỂN HÌNH 119 PHỤ LỤC II: CÁC CHƯƠNG TRÌNH TÍNH TỐN NHIỆT - ỨNG SUẤT NHIỆT 121 vii DANH MỤC CÁC HÌNH VẼ Hình 1 Đập BTĐL xây dựng giới tính đến năm 2009 [50] 11 Hình 2.1 Ứng suất nhiệt bề mặt khối bê tông 23 Hình 2.2 Biến dạng nhiệt độ & ứng suất kiềm chế khối bê tơng 24 Hình 2.3 Chênh lệch nhiệt độ khối đập 28 Hình 2.4 Các phương pháp để giải tốn nhiệt 30 Hình 2.5 Phần tử phẳng sử dụng tính tốn nhiệt 34 Hình 2.6 Phần tử chuyển vị nút phần tử tứ giác tam giác 38 Hình 3.1 Hệ số ảnh hưởng hàm lượng C3A + C3S đến nhiệt thủy hóa XM 46 Hình 3.2 Phần tử sử dụng cho phân tích nhiệt - ứng suất (PLANE55 – PLANE182) 48 Hình 3 Điều kiện biên nhiệt 49 Hình 3.4 Sơ đồ khối tính tốn nhiệt ứng suất nhiệt đập BTĐL .50 Hình 3.5 Mặt cắt ngang đập trọng lực BTĐL 51 Hình Sự phát triển cường độ kháng nén, kéo theo thời gian BTĐL Sơn La 53 Hình Sự phát triển mô đun đàn hồi theo thời gian BTĐL Sơn La 53 Hình 3.8 Tiến độ thi công đập 56 Hình 3.9 Nhiệt độ quan trắc cao trình +114,0 m .57 Hình 3.10 Mặt cắt tính toán 58 Hình 3.11 Mơ hình hình học đập 59 Hình 3.12 Mơ hình PTHH đập làm việc đồng thời với 59 Hình 3.13 Trường nhiệt độ lớn thân đập .60 Hình 3.14 Diễn biến nhiệt độ cao trình +114,0 m theo thời gian 61 Hình 3.15 Trường ứng suất lớn thân đập 61 Hình 3.16 Ứng suất lớn mặt thượng hạ lưu đập .62 Hình 3.17 So sánh nhiệt độ tính tốn thực đo cao trình +114,0 m .62 Hình 4.1 Các đập lớn đã, xây dựng Việt Nam 65 Hình 4.2 Diễn biến nhiệt độ khu vực miền núi phía Bắc .67 Hình 4.3 Diễn biến nhiệt độ khu vực Bắc Trung 68 Hình 4.4 Diễn biến nhiệt độ khu vực Nam Trung Tây Nguyên 68 Hình 4.5 Ảnh hưởng nhiệt độ khơng khí đến nhiệt độ biên khối đổ 71 viii !Dinh nghia phan tu ET,1,PLANE42 ET,2,SOLID70 !Ve diem dac trung K,1,0,0,0 K,2,BD,0,0 K,3,BDD+BC,H-HD,0 K,4,BDD+BC,H,0 K,5,BC,H,0 K,6,BC,HC,0 K,7,-3,6,0 K,8,-4,6,0 K,9,-10.22,12.22,0 K,10,-TLN,12.22,0 K,11,-TLN,-DSN,0 K,12,BD+HLN,-DSN,0 K,13,BD+HLN,10.3,0 K,14,51.51,10.3,0 K,15,47.21,6,0 K,16,46.21,6,0 A,1,2,3,4,5,6!DAP NUMMRG,ALL WPROT,0,-90,0 *DO,I,0,N_LOP-2 WPOFF,0,0,D_LOP ASBW,ALL *ENDDO ALLSEL WPAVE,0 WPROT,0,90,0 A,1,7,8,9,10,11,12,13,14,15,16,2 !NEN ALLSEL,ALL MP,EX,1,EX_DAM !Vat lieu RCC MP,NUXY,1,NUXY_DAM MP,DENS,1,DENS_DAM MP,KXX,1,KXX_DAM MP,C,1,C_DAM MP,ALPX,1,ALPX_DAM MP,EX,2,EX_F !Vat lieu nen MP,NUXY,2,NUXY_F 122 MP,DENS,2,DENS_F MP,KXX,2,KXX_F MP,C,2,C_F MP,ALPX,2,ALPX_F ASEL,S,LOC,Y,0,H AATT,1,,1,,,,, ASEL,S,LOC,Y,-DSN,0 AATT,2,,1,,,, ALLSEL,ALL ESIZE,H/30,0, MSHKEY,2 MSHAPE,0,2D AMESH,ALL TYPE,2 EXTOPT,ESIZE,1,0, EXTOPT,ACLEAR,1 EXTOPT,ATTR,1,1,1 ESEL,ALL VEXT,ALL, , ,0,0,1,,,, ALLSEL,ALL /PNUM,VOLU,1 /NUMBER,1 VPLOT /PNUM,MAT,1 /REPLOT EPLOT NROPT,FULL ESEL,S,MAT,,2,, NSLE,R,1 IC,ALL,TEMP,T_FOUN ALLSEL ASEL,S,LOC,Y,-DSN ASEL,A,LOC,X,-TLN ASEL,A,LOC,X,BD+HLN ASEL,A,LOC,Z,0 ASEL,A,LOC,Z,1 NSLA,R,1 SF,ALL,HFLUX,0 ALLSEL VSEL,S,MAT,,2 123 ASLV,R,1 ASEL,U,LOC,Y,0 ASEL,U,LOC,Z,0 ASEL,U,LOC,Z,1 ASEL,U,LOC,Y,-DSN ASEL,U,LOC,X,-TLN ASEL,U,LOC,X,HLN+BD NSLA,R,1 SF,ALL,CONV,H_FOUD,T_AIR ALLSEL,ALL ESEL,,MAT,,1 EKILL,ALL NSLE,R,1 IC,ALL,TEMP,T_CONC VSEL,S,MAT,,1 ESLV,R,1 EALIVE,ALL ALLSEL VSEL,S,LOC,Y,0,N_LOP*D_LOP ASLV,R,1 ASEL,U,LOC,Z,0,0 ASEL,U,LOC,Z,1,1 *DO,I,1,N_LOP+1 ASEL,U,LOC,Y,(I-1)*D_LOP,(I-1)*D_LOP *ENDDO NSLA,R,1 !Mien nui phia Bac T=I*(TD_LOP+TN_LOP) ND1=7E-09*(T**4) ND2=-6E-06*(T**3) ND3=0.0012*(T**2) ND4=-0.0114*T NDMT=ND1+ND2+ND3+ND4+14.858 !Bac Trung bo T=I*(TD_LOP+TN_LOP) ND1=7E-09*(T**4) ND2=-7E-06*(T**3) ND3=0.0016*(T**2) ND4=-0.0583*T NDMT=ND1+ND2+ND3+ND4+17.314 124 !Tay nguyen T=I*(TD_LOP+TN_LOP) ND1=1E-10*(T**4) ND2=+2E-07*(T**3) ND3=-0.0002*(T**2) ND4=0.0521*T NDMT=ND1+ND2+ND3+ND4+25.886 SF,ALL,CONV,H_CONC,T_AIR ALLSEL TL=N_LOP*TD_LOP+N_LOP*TN_LOP HE00=0.008*(-1.89*61.4+238.07)*(36.75*log(TL)+101.15)*0.0074*40+0.4937) ESEL,S,MAT,,1 NSLE,R,1 NSEL,S,LOC,Y,0,N_LOP*D_LOP ESLN,R,1 BFE,ALL,HGEN,,HE00, , , ALLSEL,ALL /SOLU ANTYPE,TRANS TIME,TL AUTOTS,0 DELTIM,1 KBC,0 SOLVE SAVE FINISH KEYW,PR_STRUC,1 /PREP7 ETCHG,TTS ESEL,S,MAT,,1 MPCHG,1,ALL, ESEL,S,MAT,,2 MPCHG,2,ALL, ALLSEL,ALL NSEL,S,LOC,Y,-DSN D,ALL, ,0, , , ,ALL, , , , , NSEL,S,LOC,X,-TLN NSEL,A,LOC,X,BD+HLN D,ALL, ,0, , , ,UX, , , , , ALLSEL,ALL 125 /SOL ESEL,,MAT,,1 EKILL,ALL ACEL,0,9.81,0 NSEL,S,LOC,Y,0,N_LOP*D_LOP ESLN,R,1 EALIVE,ALL ALLSEL LDREAD,TEMP, , ,TG, ,DAPBTDL','rst',' ' ALLSEL,ALL NROPT,FULL NLGEOM,ON KBC,0 TIME,TG SOLVE SAVE FINISH 126 2: TÍNH NHIỆT - ỨNG SUẤT NHIỆT ĐẬP BTĐL SƠN LA /FILNAME,DAPBTDL /TITLE, PHAN TICH NHIET-UNG SUAT DAP BE TONG DAM LAN SON LA H=131.1 !CHIEU CAO DAP BD=105 !BE RONG DAY DAP BDD=10 !BE RONG DINH DAP HD=0 !CHIEU CAO CO DAP BC=0 !BE RONG CHAN DAP HC=0 !CHIEU CAO CHAN DAP TLN=150 !PHAM VI NEN THUONG LUU HLN=150 !PHAM VI NEN HA LUU DSN=150 !PHAM VI NEN PHIA DUOI D_LOP=1 !M,CHIEU DAY DOT DO N=H/D_LOP N_LOP=NINT(N) !SO DOT DO TD_LOP=1 !THOI GIAN HOAN THANH DOT DO, NGAY TN_LOP=3 !THOI GIAN NGHI GIUA CAC DOT DO, NGAY TG=N_LOP*(TD_LOP+TN_LOP) !THOI GIAN DO DEN DINH DAP EX_F=1.3E10 NUXY_F=0.28 DENS_F=2600 KXX_F=3.1*24*3600/1000 C_F=0.85 EX_DAM=2.5E10 NUXY_DAM=0.2 DENS_DAM=2400 KXX_DAM=1.5*24*3600/1000 C_DAM=0.756 ALPX_DAM=0.787E-5 T_AIR=23.2+4 !NHIET DO KHONG KHI TRUNG BINH NAM H_CONC=30 !HE SO DOI LUU NHIET BE TONG - KHONG KHI H_FOUD=15 !HE SO DOI LUU NHIET NEN - KHONG KHI T_CONC=20 !NHIET DO BE TONG T_FOUN=20 !NHIET DO NEN /PREP7 ET,1,PLANE77 K,1,0,0,0 K,2,BD,0,0 K,3,BDD+BC,H-HD,0 K,4,BDD+BC,H,0 127 K,5,BC,H,0 K,6,BC,HC,0 K,7,-TLN,0,0 K,8,-TLN,-DSN,0 K,9,BD+HLN,-DSN,0 K,10,BD+HLN,0,0 A,1,2,3,4,5,6 NUMMRG,ALL WPROT,0,-90,0 *DO,I,0,N_LOP-2 WPOFF,0,0,D_LOP ASBW,ALL *ENDDO ALLSEL WPAVE,0 WPROT,0,90,0 A,1,7,8,9,10,2 MP,KXX,1,KXX_DAM MP,C,1,C_DAM MP,DENS,1,DENS_DAM MP,DENS,2,DENS_F MP,KXX,2,KXX_F MP,C,2,C_F ASEL,S,LOC,Y,0,H CM,DAP,AREA AATT,1,,1 ASEL,S,LOC,Y,0,-DSN CM,NEN,AREA AATT,2,,1 ALLSEL,ALL ESIZE,H/30,0, MSHKEY,2 MSHAPE,0,2D AMESH,ALL /SOLU ANTYPE,TRANS NROPT,FULL CMSEL,S,NEN NSLA,S,1 IC,ALL,TEMP,T_FOUN 128 CMSEL,S,DAP NSLA,S,1 IC,ALL,TEMP,T_CONC ALLSEL LSEL,S,LOC,Y,-DSN LSEL,A,LOC,X,-TLN LSEL,A,LOC,X,BD+HLN NSLL,S SF,ALL,HFLUX,0 ALLSEL ASEL,S,MAT,,2 LSLA,S,1 LSEL,U,LOC,Y,-DSN,0 NSLL,S,1 NSEL,U,LOC,X,0.1,BD-0.1 SF,ALL,CONV,H_FOUD,T_AIR ALLSEL ESEL,S,MAT,,1 EKILL,ALL !THI CONG BE TONG !DO BE TONG *DO,I,1,N_LOP ASEL,S,LOC,Y,0,(I-1)*D_LOP ESLA,S,1 NSLE,S,ALL SFDELE,ALL,CONV ALLSEL ASEL,S,LOC,Y,(I-1)*D_LOP,I*D_LOP ESLA,S,1 EALIVE,ALL ALLSEL ASEL,S,LOC,Y,0,I*D_LOP LSLA,S,1 *DO,J,1,I LSEL,U,LOC,Y,(J-1)*D_LOP *ENDDO NSLL,R,1 !NHIET DO KHONG KHI THAY DOI THEO NGAY !NHIET DO KHU VUC MIEN NUI PHIA BAC T=I*(TD_LOP+TN_LOP) 129 ND1=7E-09*(T**4) ND2=-6E-06*(T**3) ND3=0.0012*(T**2) ND4=-0.0114*T T_AIR=ND1+ND2+ND3+ND4+14.858 NHIET DO KHU VUC BAC TRUNG BO T=I*(TD_LOP+TN_LOP) ND1=7E-09*(T**4) ND2=-7E-06*(T**3) ND3=0.0016*(T**2) ND4=-0.0583*T T_AIR=ND1+ND2+ND3+ND4+17.314 !NHIET DO KHU VUC TAY NGUYEN T=I*(TD_LOP+TN_LOP) ND1=1E-10*(T**4) ND2=+2E-07*(T**3) ND3=-0.0002*(T**2) ND4=0.0521*T T_AIR=ND1+ND2+ND3+ND4+25.886 SF,ALL,CONV,H_CONC,T_AIR ALLSEL *DO,K,1,I ASEL,S,LOC,Y,(K-1)*D_LOP,K*D_LOP ESLA,S,1 NSLE,S,ALL HE00=0.008*(-1.89*70+238.07)*(36.75*log(((I-K+1)*TD_LOP+(IK)*TN_LOP))+101.15)*0.534 BF,ALL,HGEN,HE00 *ENDDO ALLSEL TIME,I*TD_LOP+(I-1)*TN_LOP DELTIM,TD_LOP AUTOTS,ON KBC,0 OUTRES,ALL,ALL SOLVE !NGHI GIUA CAC DOT ALLSEL *DO,L,1,I ASEL,S,LOC,Y,(L-1)*D_LOP,L*D_LOP 130 ESLA,S,1 NSLE,S,ALL HE00=0.008*(-1.89*70+238.07)*(36.75*log(((I-L+1)*TD_LOP+(IL+1)*TN_LOP))+101.15)*0.534 BF,ALL,HGEN,HE00 *ENDDO ALLSEL TIME,I*TD_LOP+I*TN_LOP DELTIM,TN_LOP AUTOTS,ON KBC,0 OUTRES,ALL,ALL SOLVE *ENDDO FINISH SAVE 131 3: TÍNH NHIỆT - ỨNG SUẤT NHIỆT ĐẬP BTĐL TRUNG SƠN /FILNAME,DAPBTDL /TITLE, PHAN TICH NHIET-UNG SUAT DAP BE TONG DAM LAN TRUNG SON H=84.5 !CHIEU CAO DAP BD=82.7 !BE RONG DAY DAP BDD=8 !BE RONG DINH DAP HD=4.6 !CHIEU CAO CO DAP BC=22.8 !BE RONG CHAN DAP HC=65 !CHIEU CAO CHAN DAP TLN=80 HLN=80 DSN=80 D_LOP=1 !M,CHIEU DAY LOP DO N=H/D_LOP N_LOP=NINT(N) !SO LOP DO TD_LOP=1 !THOI GIAN HOAN THANH LOP DO, NGAY TN_LOP=3 EX_F=1.3E10 NUXY_F=0.28 DENS_F=2600 KXX_F=3.1*24*3600/1000 C_F=0.85 ALPX_F=1.07E-5 EX_DAM=2.5E10 NUXY_DAM=0.2 DENS_DAM=2400 KXX_DAM=1.5*24*3600/1000 C_DAM=0.756 ALPX_DAM=0.787E-5 T_AIR=23.2+1 !NHIET DO KHONG KHI TRUNG BINH NAM H_CONC=30 !HE SO DOI LUU NHIET BE TONG - KHONG KHI H_FOUD=15 !HE SO DOI LUU NHIET NEN - KHONG KHI T_CONC=20 !NHIET DO BE TONG T_FOUN=20 !NHIET DO NEN /PREP7 !Dinh nghia phan tu ET,1,PLANE42 ET,2,SOLID70 !Ve diem dac trung 132 K,1,0,0,0 K,2,BD,0,0 K,3,BDD+BC,H-HD,0 K,4,BDD+BC,H,0 K,5,BC,H,0 K,6,BC,HC,0 K,7,-TLN,0,0 K,8,-TLN,-DSN,0 K,9,BD+HLN,-DSN,0 K,10,BD+HLN,0,0 A,1,2,3,4,5,6!DAP NUMMRG,ALL WPROT,0,-90,0 *DO,I,0,N_LOP-2 WPOFF,0,0,D_LOP ASBW,ALL *ENDDO ALLSEL WPAVE,0 WPROT,0,90,0 A,1,7,8,9,10,2 !NEN ALLSEL,ALL MP,EX,1,EX_DAM !Vat lieu RCC MP,NUXY,1,NUXY_DAM MP,DENS,1,DENS_DAM MP,KXX,1,KXX_DAM MP,C,1,C_DAM MP,ALPX,1,ALPX_DAM MP,EX,2,EX_F !Vat lieu nen MP,NUXY,2,NUXY_F MP,DENS,2,DENS_F MP,KXX,2,KXX_F MP,C,2,C_F MP,ALPX,2,ALPX_F ASEL,S,LOC,Y,0,H AATT,1,,1,,,,, ASEL,S,LOC,Y,-DSN,0 AATT,2,,1,,,, ALLSEL,ALL ESIZE,H/30,0, 133 MSHKEY,2 MSHAPE,0,2D AMESH,ALL TYPE,2 EXTOPT,ESIZE,1,0, EXTOPT,ACLEAR,1 EXTOPT,ATTR,1,1,1 ESEL,ALL VEXT,ALL, , ,0,0,1,,,, ALLSEL,ALL /PNUM,VOLU,1 /NUMBER,1 VPLOT /PNUM,MAT,1 /REPLOT EPLOT NROPT,FULL ESEL,S,MAT,,2,, NSLE,R,1 IC,ALL,TEMP,T_FOUN ALLSEL ASEL,S,LOC,Y,-DSN ASEL,A,LOC,X,-TLN ASEL,A,LOC,X,BD+HLN ASEL,A,LOC,Z,0 ASEL,A,LOC,Z,1 NSLA,R,1 SF,ALL,HFLUX,0 ALLSEL VSEL,S,MAT,,2 ASLV,R,1 ASEL,U,LOC,Y,0 ASEL,U,LOC,Z,0 ASEL,U,LOC,Z,1 ASEL,U,LOC,Y,-DSN ASEL,U,LOC,X,-TLN ASEL,U,LOC,X,HLN+BD NSLA,R,1 SF,ALL,CONV,H_FOUD,T_AIR ALLSEL,ALL 134 ESEL,,MAT,,1 EKILL,ALL NSLE,R,1 IC,ALL,TEMP,T_CONC VSEL,S,MAT,,1 ESLV,R,1 EALIVE,ALL ALLSEL VSEL,S,LOC,Y,0,N_LOP*D_LOP ASLV,R,1 ASEL,U,LOC,Z,0,0 ASEL,U,LOC,Z,1,1 *DO,I,1,N_LOP+1 ASEL,U,LOC,Y,(I-1)*D_LOP,(I-1)*D_LOP *ENDDO NSLA,R,1 SF,ALL,CONV,H_CONC,T_AIR ALLSEL TL=N_LOP*TD_LOP+N_LOP*TN_LOP HE00=0.008*(-1.89*72.5+238.07)*(36.75*log(TL)+101.15)*0.534 ESEL,S,MAT,,1 NSLE,R,1 NSEL,S,LOC,Y,0,N_LOP*D_LOP ESLN,R,1 BFE,ALL,HGEN,,HE00, , , ALLSEL,ALL /SOLU ANTYPE,TRANS TIME,TL AUTOTS,0 DELTIM,1 KBC,0 SOLVE SAVE FINISH KEYW,PR_STRUC,1 /PREP7 ETCHG,TTS ESEL,S,MAT,,1 MPCHG,1,ALL, 135 ESEL,S,MAT,,2 MPCHG,2,ALL, ALLSEL,ALL NSEL,S,LOC,Y,-DSN D,ALL, ,0, , , ,ALL, , , , , NSEL,S,LOC,X,-TLN NSEL,A,LOC,X,BD+HLN D,ALL, ,0, , , ,UX, , , , , ALLSEL,ALL /SOL ESEL,,MAT,,1 EKILL,ALL ACEL,0,9.81,0 NSEL,S,LOC,Y,0,N_LOP*D_LOP ESLN,R,1 EALIVE,ALL ALLSEL LDREAD,TEMP, , ,TL, ,'DAPBTDL','rst',' ' ALLSEL,ALL NROPT,FULL NLGEOM,ON KBC,0 TIME,TL SOLVE SAVE FINISH QUẢ THÍ NGHIỆM THÀNH PHẦN KHOÁNG CỦA TỪNG LOẠI XI MĂNG 136