(Luận văn thạc sĩ) phân tích ứng xử dầm bê tông cốt cứng bằng phương pháp phần tử hữu hạn

i LỜI CẢM ƠN Thực tế cho thấy, thành công gắn liền với hỗ trợ, giúp đỡ người xung quanh cho dù giúp đỡ hay nhiều, trực tiếp hay gián tiếp Trong suốt thời gian từ bắt đầu làm luận văn đến nay, nhận quan tâm, bảo, giúp đỡ thầy cơ, gia đình đờng nghiệp Với lịng biết ơn vô sâu sắc, học viên xin gửi lời cảm ơn chân thành đến quý Thầy Cô Trường Đại học Kinh Tế Công Nghiệp Long An dùng tri thức tâm huyết để truyền đạt cho tơi vốn kiến thức q báu suốt thời gian học tập tại trường Học viên xin tỏ lịng kính trọng biết ơn sâu sắc thầy giáo PGS.TS Trương Tích Thiện trực tiếp tận tình hướng dẫn cung cấp tài liệu thông tin khoa học cần thiết cho luận văn Nhờ có lời hướng dẫn, dạy bảo đó, luận văn tơi hồn thành cách suất sắc Xin chân thành cảm ơn Lãnh đạo Trường Đại học Kinh Tế Công Nghiệp Long An, khoa Kiến trúc xây dựng tạo điều kiện cho tơi hồn thành tốt cơng việc nghiên cứu khoa học Bài luận văn thực hiện 06 tháng Vì vốn kiến thức có hạn, nên khơng tránh khỏi thiếu sót, tơi mong nhận ý kiến đóng góp q Thầy Cơ anh (chị) học viên lớp luận văn hoàn thiện HỌC VIÊN THỰC HIỆN Bạch Thanh Vũ Luan van ii BẢN CAM KẾT Ngoài kết tham khảo từ cơng trình khác ghi luận văn, xin cam kết luận văn tơi thực hiện luận văn nộp tại Trường Đại học Kinh tế Công nghiệp Long An Tôi xin cam đoan rằng: Số liệu kết nghiên cứu luận văn hoàn toàn trung thực chưa sử dụng cơng bố cơng trình khác Mọi giúp đỡ cho việc thực hiện luận văn cảm ơn thơng tin trích dẫn luận văn ghi rõ nguồn gốc HỌC VIÊN THỰC HIỆN Bạch Thanh Vũ Luan van iii TĨM TẮT Luận văn sử dụng chương trình ANSYS – chương trình tính tốn xây dựng dựa phương pháp phần tử hữu hạn – để phân tích ứng xử dầm bê tông cốt cứng Luận văn thực hiện hai tốn, kết tính tốn có so sánh với kết từ thí nghiệm ABSTRACT In this thesis, ANSYS – the engineering simulation software based on the finite element method – is used to analyze behaviour of reinforced concrete beam Two model is analyzed, results from ANSYS is also compared with result of experiment Luan van iv MỤC LỤC LỜI CẢM ƠN i BẢN CAM KẾT ii TÓM TẮT iii ABSTRACT iii MỤC LỤC iv DANH MỤC BẢNG vi DANH MỤC HÌNH ẢNH .vii DANH MỤC CÁC TỪ VIẾT TẮT x Chương 1.TỔNG QUAN 1.1 Giới thiệu chung kết cấu liên hợp thép – bê tông 1.2 Tình hình nghiên cứu đề tài 1.2.1 Nghiên cứu nước 1.2.2 Nghiên cứu nước 1.3 Tính cấp thiết đề tài 1.4 Mục tiêu nghiên cứu 1.5 Đối tượng phạm vi nghiên cứu 1.6 Phương pháp nghiên cứu 1.7 Lợi ích đề tài 1.7.1 Lợi ích khoa học 1.7.2 Lợi ích thực tiễn Chương 2.CƠ SỞ LÝ THUYẾT 2.1 Giới thiệu vê vật liệu bê tông 2.1.1 Cường độ chịu nén 2.1.2 Quan hệ ứng suất -biến dạng cho bê tông chịu tải trọng nén trục 2.1.3 Quan hệ ứng suất - biến dạng bê tông thường 12 2.1.4 Đường cong ứng suất - biến dạng bê tông chịu kéo 15 2.1.5 Các đặc trưng bê tông 16 2.1.6 Quan hệ ứng suất - biến dạng cho bê tông chịu tải trọng hai trục 18 2.2 Các mơ hình vật liệu bê tơng 20 Luan van v 2.2.1 Mơ hình Hognestad (1951) 20 2.2.2 Mô hình Mendis (2001) 22 2.3 Các đặc trưng tính thép 24 2.4 Mơ hình vật liệu thép 25 2.5 Phương pháp phần tử hữu hạn 26 2.5.1 Giới thiệu 26 2.5.2 Xây dựng mơ hình phần tử hữu hạn ANSYS 28 2.6 Phân tích kết cấu chương trình ANSYS 34 Chương PHÂN TÍCH ỨNG XỬ DẦM BÊ TƠNG CỐT CỨNG BẰNG CHƯƠNG TRÌNH ANSYS 36 3.1 Mơ hình tốn dầm 36 3.1.1 Mô hình hình học 36 3.1.2 Thông số vật liệu 37 3.1.3 Mơ hình PTHH 41 3.1.4 Tải trọng – điều kiện biên 42 3.1.5 Kết phân tích 43 3.2 Khảo sát ảnh hưởng thép hình đến khả chịu lực dầm 49 3.3 Nhận xét 50 Chương 4.KẾT LUẬN VÀ KIẾN NGHỊ 51 4.1 Kết luận 51 4.1.1 Những nội dung làm 51 4.1.2 Những nội dung cịn thiếu sót 51 4.2 Kiến nghị 52 TÀI LIỆU THAM KHẢO 53 PHỤ LỤC CODE ANSYS 55 Luan van vi DANH MỤC BẢNG Bảng 2.1 Giá trị α dùng để tính đổi kết thử cường độ mẫu chuẩn Bảng 3.1 Thơng số kích thước thiết diện ngang cường độ vật liệu dầm FEB [2] 37 Bảng 3.2 Bảng cấp phối bê tông [2] 37 Bảng 3.3 Thông số vật liệu bê tông [2] 38 Bảng 3.4 Thông số vật liệu thép 39 Bảng 3.5 Thông số vật liệu thép 40 Bảng 3.6 Kết độ võng dầm tính ANSYS với mức tải khác 43 Bảng 3.7 Bảng thông số ds bf 49 Luan van vii DANH MỤC HÌNH ẢNH Hình 1.1 Các dạng tiết diện dầm – cột liên hợp thép bê tông [3] Hình 2.1 Máy thí nghiệm nén mẫu trụ tròn Hình 2.2 Mẫu thử dùng cho thí nghiệm nén Hình 2.3 Mơ hình thí nghiệm nén Hình 2.4 Quan hệ ứng suất biến dạng bê tông chịu kéo nén trục 11 Hình 2.5 Hệ số giãn nở cho bê tông chịu tải trọng nén trục 12 Hình 2.6 (a) Quan hệ ứng suất - biến dạng; (b) Ảnh hưởng tốc độ gia tải 13 Hình 2.7 Thí nghiệm kiểm sốt theo mức độ gia tải mức độ biến dạng 14 Hình 2.8 Đường cong ứng suất-biến dạng bê tơng chịu kéo 16 Hình 2.9 Thí nghiệm khối trụ chẻ 16 Hình 2.10 Phương pháp xác định loại mô đun đàn hồi 17 Hình 2.11 Sự phân bố ứng suất hai trục 18 Hình 2.12 Độ bền dạng phá hỏng bê tơng chịu ứng suất hai trục 19 Hình 2.13 Mơ hình ứng suất biến dạng Hognestad Todeschini 21 Hình 2.14 Mơ hình ứng suất - biến dạng Mander 21 Hình 2.15 Mơ hình Mendis 23 Hình 2.16 Quan hệ ứng suất - biến dạng thép 25 Hình 2.17 Quan hệ ứng suất - biến dạng thép ASTM A615 Grade60 26 Hình 2.18 Phần tử khối bê tơng Solid65 29 Hình 2.19 Các mơ hình cốt thép bê tơng cốt thép 29 Hình 2.20 Mơ hình phần tử LINK180 ANSYS 30 Hình 2.21 Mơ hình phần tử SHELL181 ANSYS 31 Hình 2.22 Đường cong ứng suất- biến dạng bê tông 34 Luan van viii Hình 2.23 Đường cong ứng suất- biến dạng thép 34 Hình 3.1 Mơ hình thí nghiệm dầm [2] 36 Hình 3.2 Mặt cắt ngang dầm FEB [2] 37 Hình 3.3 Đường cong ứng suất biến dạng vật liệu bê tông 300 [2] 38 Hình 3.4 Đường cong ứng suất biến dạng thép thu từ thí nghiệm kéo đơn trục [2] 39 Hình 3.5 Mơ hình đàn – dẻo song tuyến tính cho thép 40 Hình 3.6 Đường cong ứng suất biến dạng thép thu từ thí nghiệm kéo đơn trục [2] 41 Hình 3.7 Mơ hình ½ dầm ANSYS 42 Hình 3.8 Mơ hình PTHH dầm 42 Hình 3.9 Mơ hình PTHH cốt thép thép hình dầm 42 Hình 3.10 Hình ảnh minh họa điều kiện biên tải trọng 43 Hình 3.11 Đờ thị biểu diễn kết độ võng dầm theo giá trị lực tăng dần [2] 44 Hình 3.12 Đồ thị biểu diễn kết độ võng dầm tính tốn từ ANSYS theo giá trị lực tăng dần 44 Hình 3.13 Kết độ võng dầm với mức tải 5000N 45 Hình 3.14 Kết độ võng dầm với mức tải 10000N 45 Hình 3.15 Kết độ võng dầm với mức tải 13300N 46 Hình 3.16 Trường ứng suất tương đương von-Mises cốt thép tại mức tải 13.3kN 47 Hình 3.17 Trường ứng suất tương đương von-Mises thép hình tại mức tải 13.3kN 47 Hình 3.18 Vết rạn xuất hiện bê tông tại mức tải 13.3kN 48 Hình 3.19 So sánh kết độ võng dầm ANSYS thí nghiệm 49 Luan van ix Hình 3.20 Kết độ võng dầm tương ứng với trường hợp thép hình chữ I có kích thước khác 50 Luan van x DANH MỤC CÁC TỪ VIẾT TẮT Tên đầy đủ PP PTHH Ý nghĩa Phương pháp phần tử hữu hạn FEB Fully Enscaped beam FEM Finite Element Method E Module đàn hồi  Hệ số Poisson fy Cường độ tính tốn thép theo giới hạn chảy vật liệu G Mô đun trượt  , , Các ứng suất fC Độ bền nén đơn trục bê tông Bề mặt phá hủy xác định từ năm S  ah thông số f t , f c , f cb , f1 , f Trạng thái ứng suất thủy tĩnh điều kiện môi trường xung quanh Luan van 43 Tải trọng: theo [2], tải tăng từ từ với bước tải 5kN Do để so sánh với kết thí nghiệm, tải đặt vào mơ hình ANSYS tăng từ từ với bước tải 5kN Tại bề mặt đặt tải, ta đặt tải 100kN/diện tích đặt tải thiết lập cho ANSYS giải với 20 bước giải (mỗi bước giải lực tăng lên tương đương 5kN) Hình 3.10 Hình ảnh minh họa điều kiện biên tải trọng 3.1.5 Kết phân tích 3.1.5.1 Kết độ võng dầm Kết độ võng dầm tính từ chương trình ANSYS với mức tải tăng dần trình bày bảng 3.6 hình 3.13 đến 3.15 Theo [2], tại mức tải 13300N, dầm bắt đầu xuất hiện dấu hiệu rạn nứt (hình 3.11) Bảng 3.6 Kết độ võng dầm tính ANSYS với mức tải khác Tải trọng (kN) 5000 10000 13300 15000 20000 25000 30000 Độ võng dầm (mm) 0.16 0.42 0.97 1.24 1.66 2.21 3.02 Luan van 44 Hình 3.11 Đờ thị biểu diễn kết độ võng dầm theo giá trị lực tăng dần [2] 35 30 Tải (kN) 25 20 15 ANSYS 10 0 0.5 1.5 2.5 3.5 Độ võng dầm (mm) Hình 3.12 Đồ thị biểu diễn kết độ võng dầm tính tốn từ ANSYS theo giá trị lực tăng dần Luan van 45 Hình 3.13 Kết độ võng dầm với mức tải 5000N Hình 3.14 Kết độ võng dầm với mức tải 10000N Luan van 46 Hình 3.15 Kết độ võng dầm với mức tải 13300N 3.1.5.2 Đánh giá độ bền dầm Điểm mạnh chương trình ANSYS dễ dàng mơ trường ứng suất kết cấu Kết ứng suất von-Mises cốt thép thép hình với mức tải 13.3kN thể hiện hình 3.15 3.16 Các kết cho thấy ứng suất tương đương cực đại thép nhỏ giới hạn bền cốt thép lẫn thép hình nằm vùng đàn hời Bê tơng vật liệu dịn, ta cần tính tốn ứng suất tương đương Mohr cực đại để đánh giá độ bền bê tông Trong ANSYS, xuất kết ứng suất thứ thứ kết hợp với cơng thức (3.1) để tìm ứng suất tương đương Morh cực đại mô hình Tại mức tải 13.3kN, ứng suất uất tương đương Morh cực đại mơ hình có giá trị 3.14MPa Giá trị vượt giới hạn phá hủy kéo đơn trục bê tông lúc xuất hiện vết rạn nứt Điều phù hợp với kết thí nghiệm Luan van 47 f mohr  f1   ft  f  fc  (3.1) Hình 3.16 Trường ứng suất tương đương von-Mises cốt thép tại mức tải 13.3kN Hình 3.17 Trường ứng suất tương đương von-Mises thép hình tại mức tải 13.3kN Luan van 48 Hình 3.18 Vết rạn xuất hiện bê tông tại mức tải 13.3kN 3.1.5.3 So sánh kết ANSYS thí nghiệm Kết so sánh độ võng dầm trình bày hình 3.16 Kết ANSYS thí nghiệm phù hợp tại mức tải 13.3kN, 15kN 30kN, cịn tại mức tải khác sai số hai phương pháp lớn Sai số mơ hình ANSYS, bỏ qua trượt thép bê tơng Vì thực tế, xảy trượt thép bê tơng độ cứng uốn dầm giảm đáng kể Ngoài ra, lý tưởng hóa mơ hình vật liệu ảnh hưởng đến kết tính tốn Theo tiêu chuẩn ACI 318M-08, mơ đun đàn hồi E c bê tông chọn theo công thức Ec   0.33 1 Ecm Khi khai báo mơ hình vật liệu bê tơng ANSYS, học viên chọn mô đun đàn hồi E c bê tông giá trị Ecm Đây nguyên nhân làm giảm độ tin cậy kết mô Luan van 49 35 30 25 Tải (kN) 20 ANSYS 15 EXP 10 0 0.5 1.5 2.5 3.5 Độ võng dầm (mm) Hình 3.19 So sánh kết độ võng dầm ANSYS thí nghiệm 3.2 Khảo sát ảnh hưởng thép hình đến khả chịu lực dầm Mơ hình dầm FEB tốn mục 3.1 sử dụng lại, mơ hình vật liệu giữ cũ có thay đổi kích thước ds bf Học viên đề xuất hai phương án (A2 A3) thay đổi kích thước ds bf Bảng 3.7 Bảng thơng số ds bf Mơ hình ds (mm) bf (mm) A 200 146 A2 262 200 A3 150 146 Kết độ võng dầm 03 trường hợp trình bày hình 3.20 Kết thể hiện rõ độ cứng uốn dầm tăng lên tăng kích thước ds bf Luan van 50 35 30 Tải (kN) 25 20 A3 15 A2 10 A 0 Độ võng dầm (mm) Hình 3.20 Kết độ võng dầm tương ứng với trường hợp thép hình chữ I có kích thước khác 3.3 Nhận xét Luận văn phân tích ứng xử mơ hình dầm tác động tải tĩnh, giá trị tải tĩnh chia thành nhiều bước lực để tính Điều chưa phù hợp với điều kiện thí nghiệm thí nghiệm tải tăng theo thời gian Đây nguyên nhân gây sai số kết theo ANSYS kết thí nghiệm Ngồi ra, việc lựa chọn giá trị module đàn hồi cho bê tông vấn đề cần lưu ý, khoảng giá trị cho phép lớn Khi tăng kích thước mặt cắt thép hình độ võng dầm giảm đi, độ cứng uốn dầm tăng lên, điều phù hợp với thực tế Luan van 51 Chương KẾT LUẬN VÀ KIẾN NGHỊ 4.1 Kết luận Kỹ thuật công nghệ ngày phát triển nhanh, yêu cầu người điều kiện sống ngày cao Các toán kỹ thuật đổi liên tục ngày phức tạp hơn.Do đó, nghiên cứu tìm hiểu phương pháp, cơng cụ chương trình tính tốn hiện đại cần thiết Hịa tinh thần đó, học viên lựa chọn đề tài ứng dụng chương trình ANSYS để phân tích ứng xử dầm bê tơng cốt cứng 4.1.1 Những nội dung làm Luận văn hoàn thành yêu cầu đặt đề cương Học viên tìm hiểu đượccác mơ hình vật liệu bê tông, thép phần tử dùng để mơ hình thép, bê tơng ANSYS cách thức mơ hình, cách thức kết nối phần tử với Luận văn ứng dụng phần mềm ANSYS để xây dựng mơ hình dầm bê tơng cốt cứng phân tích ứng xử dầm tác động tải tĩnh Các kết tính tốn từ Ansys kiểm chứng với kết thực nghiệm.Các toán luận văn thực hiện dựa giả thuyết tải tĩnh thí nghiệm thực tế tải trọng tăng theo thời gian Sự khác biệt dẫn đến sai số đáng kể kết tính theo ANSYS kết thí nghiệm Ảnh hưởng kích thước thép hình đến khả chịu lực dầm khảo sát luận văn Kết tính tốn cho thấy tăng kích thước thép hình độ võng dầm giảm Kết hoàn toàn phù hợp với thực tế 4.1.2 Những nội dung cịn thiếu sót - Các tồn luận văn dừng lại phân tích tải tĩnh, chưa xét đến toán tải theo thời gian - Chưa tìm mơ hình vật liệu bê tơng tốt cho toán luận văn Luan van 52 - Tuy luận văn có phân tích dự đốn xuất vết rạn mơ hình dầm, nhiên kết sai khác nhiều so với thực nghiệm 4.2 Kiến nghị Với kết đạt luận văn, tác giả xin đưa vài kiến nghị: - Thực giải toán nhiều lần để tìm cách tính mơ đun đàn hồi bê tơng hợp lý - Khi phân tích kết cấu Bê tông cốt thép, để đơn giản trình tính tốn, ta bỏ qua trượt bê tông cốt thép - Một lưu ý sử dụng ANSYS để phân tích tốn BTCT: không nên đặt tải trực tiếp lên bê tông Luan van 53 TÀI LIỆU THAM KHẢO [1] P V Hội, Kết cấu liên hợp thép bêtông dùng nhà cao tầng, Hà Nội: Nhà xuất Khoa học Kỹ thuật, 2006 [2] Huỳnh Phúc Linh, Hồ Hữu Chỉnh Độ bền độ võng nứt dầm liên hợp có thép hình nằm hồn tồn bê tơng cốt thép Kỷ yếu Hội nghị Khoa học Công nghệ lần thứ 12 Trường đại học Bách Khoa – ĐHQG HCM, 10/2011 [3] E 4, Design of composite steel and concrete structures: Part 1-1: General rules and rules, UK, 2005 [4] Đinh Văn Thuật, Phạm Văn Hội, “Giải pháp kết cấu liên hợp thép bê tông cho nhà nhiều tầng Việt Nam” [5] Yassin A.Y.M., Nethercot D.A., Cross-sectional properties complex composite beams, Engineering structures 29, 2007 [6] TCXDVN 338:2005, Kết cấu thép - Tiêu chuẩn thiết kế [7] EC4: ENV 1994-1-1: Eurocode 4: Design of Composite Steel and Concrete Structures Part 1.1: General rules and rules for buildings [8] Uy B., Liew J.Y.R., "Chap 51," in Composite steel-concrete structures, CRC Press, 2003, p 19 [9] "ACI Institute, ACI 318M-08 Building code requirements for structural concrete and commentary, American Concrete Institute, Farmington Hills, Michigan, 2008." [10] "ANSYS Inc - Theory reference for ANSYS and ANSYS workbench release 17, 2017." [11] Jarmai K., Farkas J., Liszkai T., "Optimum design and comparison of hollow flange beams,," in Proceedings of the Eighth International Symposium, Singapore,, 1998 [12] E.J Sapountzakis, G.C Tsiatas, "Elastic flexural buckling analysis of composite beams of variable cross-section by BEM," Engineering Structures, vol 29, p 675–681, 2007 Luan van 54 [13] N SERES, BEHAVIOUR AND RESISTANCE OF CONCRETE ENCASED EMBOSSMENTS IN COMPOSITE FLOORS, Budapest: PhD Thesis, 2012 [14] M AL-Bdoor, Behavior of Composite Columns Subjected to Lateral Cyclic Loading, Master thesis, 2013 [15] K D Tsavdaridis, "Experimental and computational study of the vertical shear behaviour of partially encased perforated steel beams," Engineering structures, vol 56, pp 805-822 Luan van 55 PHỤ LỤC CODE ANSYS /prep7 b=250 h=400 l=4000/2 b1=52 bf=146 tf=8 tw=6 d=100 lf=1400 dt=14 area=153.9380 s1=10 s2=5 s3=2 s4=50 Ad=3.14*(dt**2)/4 blc4,0,0,b,h ADELE,1 es=1 ! MODEL K,5,198,55 K,6,52,55 K,7,52,345 K,8,198,345 LANG,1,5,90 LANG,2,5,90 LANG,1,6,90 LANG,4,6,90 LANG,4,7,90 LANG,3,7,90 LANG,3,8,90 LANG,7,8,90 K,17,52,100 K,18,52,300 K,19,198,100 K,20,198,300 K,21,125,100 K,22,125,300 A,28,5,9,30 A,5,10,2,9 L,18,22 L,22,20 L,17,21 L,21,19 L,22,21 K,31,52,100,-l !K,32,125,200,-4000 !L,31,32 L,31,17 LANG,13,18,90 LANG,26,17,90 LANG,7,19,90 LANG,30,20,90 !MATERTIAL ET,1,SOLID65 !khai bao vat lieu be tong R,1,0,0,0,0,0,0 RMORE,0,0,0,0,0 L,7,8 LANG,34,22,90 L,5,6 LANG,37,21,90 LANG,17,27,90 LANG,9,28,90 A,4,14,7,13 A,7,27,29,14 A,8,15,29,27 A,16,3,15,8 A,7,13,23,18 A,27,7,18,22 A,8,27,22,20 A,16,8,20,26 A,18,23,24,17 A,18,22,21,17 A,22,20,19,21 A,20,26,25,19 A,24,17,6,12 A,17,21,28,6 A,21,19,5,28 A,19,25,10,5 A,12,6,11,1 A,6,11,30,28 Luan van MP,EX,1,16700 MP,PRXY,1,0.25 TB,CONCR,1,1,9 TBDATA,1,0.3,0.8,1.06 5,19.728 TB,MISO,1,1,19,0 TBTEMP,22 TBPT,,0.0001,1.67 TBPT,,0.0002,3.34 TBPT,,0.0003,5.01 TBPT,,0.0004,6.68 TBPT,,0.0005,8.35 TBPT,,0.0006,10.02 TBPT,,0.0007,11.69 TBPT,,0.0008,13.36 TBPT,,0.0009,15.03 TBPT,,0.001,15.353 TBPT,,0.0011,16.281 TBPT,,0.0012,17.099a TBPT,,0.0013,17.806 TBPT,,0.0014,18.402 TBPT,,0.0015,18.888 TBPT,,0.0016,19.264 TBPT,,0.0017,19.529 TBPT,,0.0018,19.683 56 TBPT,,0.0019,19.728 ! khai bao vat lieu thep hinh ET,2,181 MP,EX,2,2e5 MP,PRXY,2,0.3 sectype,1,shell !be day ban canh secdata,tf sectype,2,shell !be day ban bung secdata,tw TB,BISO,2,1,2 TBDATA,,460,2100 !khai bao vat lieu thep doc ET,3,link180 MP,EX,3,2e5 MP,PRXY,3,0.3 KEYOPT,1,3,1 sectype,3,link SECDATA,area !KHAI BAO PLANE 182 ET,4,PLANE182 !Khai bao vat support ET,5,solid185 MP,EX,3,2e5 MP,PRXY,3,0.3 !MESHING type,4 MAT,3 lsel,s,,,21,24,1 lesize,all,,,s2 lsel,all lesize,25,,,s1 lieu lesize,46,,,s1 lesize,47,,,s1 lesize,30,,,s1 lesize,26,,,s1 lesize,34,,,s2 lesize,35,,,s2 lesize,40,,,s2 lesize,42,,,s2 lesize,17,,,s2 lesize,38,,,s2 lesize,37,,,s2 lesize,9,,,s2 lesize,42,,,s2 lesize,15,,,s3 lesize,14,,,s3 lesize,27,,,s3 lesize,29,,,s3 lesize,12,,,s3 lesize,1,,,s3 lesize,3,,,s3 lesize,20,,,s3 lesize,33,,,s3 lesize,31,,,s3 lesize,8,,,s3 lesize,5,,,s3 lesize,4,,,s3 lesize,16,,,s3 lesize,41,,,s3 lesize,18,,,s3 lesize,19,,,s3 lesize,13,,,s3 lesize,44,,,s3 lesize,36,,,s3 lesize,45,,,s3 lesize,32,,,s3 lesize,28,,,s3 lesize,48,,,s3 lesize,39,,,s3 lesize,49,,,s3 lesize,7,,,s3 Luan van lesize,11,,,s3 lesize,10,,,s3 lesize,43,,,s3 lesize,6,,,s3 lesize,2,,,s3 amesh,all ASEL,S,LOC,Z,0 CM,ALTEM,AREA type,1 mat,1 esize,s4 vdrag,ALL,,,,,,50 CMSEL,S,ALTEM ACLEAR,ALL ALLSEL K,62,52,300 K,63,125,300 k,64,198,300 k,65,52,100 k,66,125,100 k,67,198,100 K,68,52,300,-l K,69,125,300,-l k,70,198,300,-l k,71,52,100,-l k,72,125,100,-l k,73,198,100,-l type,2 mat,2 secnum,1 A,62,63,69,68 A,63,64,70,69 A,65,66,72,71 A,66,67,73,72 lesize,130,,,s2 lesize,134,,,s2 lesize,141,,,s2 lesize,137,,,s2 57 lesize,139,,,s2 lesize,143,,,s2 lesize,132,,,s2 lesize,136,,,s2 lesize,133,,,l/s4 lesize,131,,,l/s4 lesize,135,,,l/s4 lesize,140,,,l/s4 lesize,138,,,l/s4 lesize,142,,,l/s4 asel,s,,,90,93,1 amesh,all allsel secnum,2 A,63,66,72,69 lesize,144,,,s1 lesize,145,,,s1 asel,s,,,94 amesh,94 !amesh,all allsel !THEP DOC K,74,52,55 K,75,198,55 K,76,52,345 K,77,198,345 K,78,52,55,-l K,79,198,55,-l K,80,52,345,-l K,81,198,345,-l L,74,78 L,75,79 L,76,80 L,77,81 type,3 mat,3 secnum,3 lsel,s,,,146,149,1 lesize,all,,,l/s4 lmesh,all allsell ! Mo hinh va chia luoi support k,82,0,400,-1350 k,83,250,400,-1350 k,84,0,400,-1450 k,85,250,400,-1450 k,86,0,0,-150 k,87,250,0,-150 k,88,0,0,-250 k,89,250,0,-250 a,82,83,85,84 a,86,87,89,88 type,5 mat,3 VEXT,95, , ,0,30,0,,,, VEXT,96, , ,0,-30,0,,,, wpoff,52 wprot,,,90 wpoff,,,73 FLST,2,2,6,ORDE,2 FITEM,2,21 FITEM,2,-22 VSBW,P51X TYPE, MAT, REAL, ESYS, SECNUM, , SFA,118,1,PRES,1.2 SFA,108,1,PRES,1.2 DA,122,UX, DA,122,Uy, !DA,122,Uz DA,132,UX, DA,132,Uy, !DA,132,Uz DA,142,UX, DA,142,Uy, !DA,142,Uz DA,137,UX, DA,137,Uy, !DA,137,Uz asel,s,loc,z,-2000 !Da,all,ux da,all,uz alls /sol CNVTOL,U,1e6,1e-6,0,, CNVTOL,,-1,1 AUTOTS,0 NSUBST,30 OUTRES,ERASE OUTRES,ALL,ALL TIME,1 ncnv,1, solve /POST1 !* /EFACET,1 PLNSOL, U,Y, 0,1.0 !Load NLdiag,nrre,stat,1 LSCLEAR,ALL SFA,113,1,PRES,1.2 SFA,128,1,PRES,1.2 Luan van ... cốt cứngcịn ít, đa phần tập trung vào phương pháp nghiên cứu theo tiêu chuẩn tính nước ngồi Vì vậy, luận văn này, tác giả chọn đề tài ? ?Phân tích ứng xử dầm bê tông cốt cứng phương pháp phần tử hữu. .. tông cốt cứng (PEB) Phương pháp phần tử hữu hạn (PTHH) dùng để phân tích nứt cho mơ hình, kết kiểm chứng thực nghiệm - Eva Maria Regy, Basil Sabu (2017) nghiên cứu ứng xử dầm bê tông cốt cứng tác... văn - Tìm kiếm tài liệu lý thuyết báo khoa học phân tích ứng xử dầm bê tông cốt cứng, làm sở so sánh với kết phân tích từ phần mềm - Sử dụng chương trình ANSYS để mơ hình phân tích dạng ứng xử