(Luận văn) 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 lu Xin chân thành cảm ơn Lãnh đạo Trường Đại học Kinh Tế Công Nghiệp Long an 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 va 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 quý 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ũ 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 lu an va n Bạch Thanh Vũ 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 toá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 lu an va n 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 lu an 1.3 Tính cấp thiết đề tài va 1.4 Mục tiêu nghiên cứu n 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 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 lu an 3.1.5 Kết phân tích 43 va 3.2 Khảo sát ảnh hưởng thép hình đến khả chịu lực dầm 49 n 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 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 lu an va n 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 lu an Hình 2.10 Phương pháp xác định loại mô đun đàn hồi 17 va n 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 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 lu Hình 3.10 Hình ảnh minh họa điều kiện biên tải trọng 43 an va Hình 3.11 Đồ thị biểu diễn kết độ võng dầm theo giá trị lực tăng dần 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 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 lu an va n x DANH MỤC CÁC TỪ VIẾT TẮT Tên đầy đủ Ý nghĩa Phương pháp phần tử hữu hạn PP PTHH 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 Trạng thái ứng suất thủy tĩnh điều kiện môi trường xung quanh lu  ah thông số f t , f c , f cb , f1 , f an va n 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) lu an Hình 3.10 Hình ảnh minh họa điều kiện biên tải trọng va 3.1.5 Kết phân tích n 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 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] lu an 35 va 30 n 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 45 lu Hình 3.13 Kết độ võng dầm với mức tải 5000N an va n Hình 3.14 Kết độ võng dầm với mức tải 10000N 46 lu Hình 3.15 Kết độ võng dầm với mức tải 13300N an 3.1.5.2 Đánh giá độ bền dầm va n Đ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 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 lu an 13.3kN va n 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 48 Hình 3.18 Vết rạn xuất hiện bê tông tại mức tải 13.3kN lu 3.1.5.3 So sánh kết ANSYS thí nghiệm an Kết so sánh độ võng dầm trình bày hình 3.16 Kết va n ANSYS thí nghiệm phù hợp tại mức tải 13.3kN, 15kN 30kN, 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ể Ngồ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ô 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 lu an giữ cũ có thay đổi kích thước ds bf Học viên đề xuất hai va phương án (A2 A3) thay đổi kích thước ds bf n 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 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 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 an tĩnh chia thành nhiều bước lực để tính Điều chưa phù hợp với điều va kiện thí nghiệm thí nghiệm tải tăng theo thời gian Đây nguyên nhân 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ế 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 lu phần tử với an va 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 n cứng phân tích ứng xử dầm tác động tải tĩnh Các kết tính toá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 toàn luận văn dừng lại phân tích tải tĩnh, chưa xét đến tốn tải theo thời gian - Chưa tìm mơ hình vật liệu bê tơng tốt cho tốn luận văn 52 - Tuy luận văn có phân tích dự đoán xuất vết rạn mơ hình dầm, nhiên kết cịn 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 tố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 q 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 lu an va n 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 lu [6] TCXDVN 338:2005, Kết cấu thép - Tiêu chuẩn thiết kế an [7] EC4: ENV 1994-1-1: Eurocode 4: Design of Composite Steel and Concrete va n 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 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 lu an va n 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 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 lu 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 TB,CONCR,1,1,9 TBDATA,1,0.3,0.8,1.06 5,19.728 n K,17,52,100 K,18,52,300 K,19,198,100 K,20,198,300 va 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 an 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 MP,EX,1,16700 MP,PRXY,1,0.25 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 CMSEL,S,ALTEM ACLEAR,ALL ALLSEL 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 K,62,52,300 K,63,125,300 k,64,198,300 k,65,52,100 k,66,125,100 k,67,198,100 n !MESHING type,4 MAT,3 lsel,s,,,21,24,1 lesize,all,,,s2 lsel,all lesize,25,,,s1 lieu type,1 mat,1 esize,s4 vdrag,ALL,,,,,,50 va !Khai bao vat support ET,5,solid185 MP,EX,3,2e5 MP,PRXY,3,0.3 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 an !KHAI BAO PLANE 182 ET,4,PLANE182 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 lu !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 lesize,46,,,s1 lesize,47,,,s1 lesize,30,,,s1 lesize,26,,,s1 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 type,3 mat,3 secnum,3 lsel,s,,,146,149,1 lesize,all,,,l/s4 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 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 n L,74,78 L,75,79 L,76,80 L,77,81 k,86,0,0,-150 k,87,250,0,-150 k,88,0,0,-250 k,89,250,0,-250 va K,78,52,55,-l K,79,198,55,-l K,80,52,345,-l K,81,198,345,-l ! 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 an !THEP DOC K,74,52,55 K,75,198,55 K,76,52,345 K,77,198,345 SFA,118,1,PRES,1.2 SFA,108,1,PRES,1.2 lu secnum,2 A,63,66,72,69 lesize,144,,,s1 lesize,145,,,s1 asel,s,,,94 amesh,94 !amesh,all allsel lmesh,all allsell TYPE, MAT, REAL, ESYS, SECNUM, , /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