Tự động hóa thiết kế tối ưu kết cấu thép dạng dàn bằng phương pháp phần tử hữu hạn có kể đến các ràng buộc theo tiêu chuẩn việt nam

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc Lập – Tự Do – Hạnh Phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên : NGUYỄN HỮU THỊNH Ngày,tháng, năm sinh : 04/08/1975 Chuyên ngành : XÂY DỰNG DÂN DỤNG & CÔNG NGHIỆP Khóa : 11 (Năm học 2000 – 2003) Phái : Nam Nơi sinh : Quảng Nam Ms ngành : 23.04.10 I TÊN ĐỀ TÀI : TỰ ĐỘNG HÓA THIẾT KẾ TỐI ƯU KẾT CẤU THÉP DẠNG DÀN BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN CÓ KỂ ĐẾN CÁC RÀNG BUỘC THEO TIÊU CHUẨN VIỆT NAM II NHIỆM VỤ VÀ NỘI DUNG : Nhiệm vụ : Lập chương trình tự động hóa thiết kế tối ưu kết cấu thép dạng dàn (dàn phẳng, dàn không gian) thông qua việc nghiên cứu ứng dụng tóan tối ưu hóa sử dụng thuật tóan đơn hình việc phân tích kết cấu phương pháp phần tử hữu hạn Nội dung : - Chương - Chương - Chương - Chương hình tương tích - Chương Khung Không Gian - Chương - Chương - Chương Dàn - Chương : : : : Phần Tổng Quan Bài Toán Tối Ưu Và Các Phương Pháp Thiết Kế Tối Ưu Phương Pháp Quy Hoạch Tuyến Tính Và Thuật Toán Đơn Hình Cơ sở bước phân tích chung phương pháp phần tử hữu hạn sử dụng mô : Thiết Lập Chương Trình Máy Tính Cho Việc Tự Động Hóa Phân Phân Tích Kết Cấu : Các Yêu Cầu Chung Cho Việc Tính Toán Thiết Kế Kết Cấu Dàn Thép : Thiết Lập Bài Toán Thiết Kế Tối Ưu Cho Kết Cấu Dàn Thép Không Gian : p Dụng Số Vào Chương Trình Tự Động Hóa Thiết Kế Tối Ưu Kết Cấu Thép Dạng : Kết Luận Kiến Nghị III NGÀY GIAO NHIỆM VỤ : 25/11/2002 IV NGÀY HÒAN THÀNH NHIỆM VỤ : 07/06/2003 V HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN : Tiến sỹ BÙI CÔNG THÀNH VI HỌ VÀ TÊN CÁN BỘ NHẬN XÉT : VII HỌ VÀ TÊN CÁN BỘ NHẬN XÉT : Cán Bộ Hướng Dẫn Cán Bộ Nhận Xét Cán Bộ Nhận Xét Ts.BÙI CÔNG THÀNH Nội dung đề cương Luận Văn Thạc Só Hội Đồng Chuyên Ngành thông qua TRƯỞNG PHÒNG ĐÀO TẠO – SAU ĐẠI HỌC Ngày……tháng… năm 2003 CHỦ NHIỆM NGÀNH Lời Cảm Ơn Những kiến thức quy báu mà em gặt hái sau trình dài học tập, nghiên cứu hòan thành luận văn nhờ dạy bảo tận tình thầy, cô ban giảng huấn lớp cao học ngành Xây Dựng Dân Dụng -– Công Nghiệp khóa 11 Công lao to lớn Thầy Cô, em xin thầm khắc ghi lòng tỏ lòng biết ơn sâu sắc nhất, kính trọng Em xin chân thành cảm ơn Thầy Tiến sỹ Bùi Công Thành -– Phó chủ nhiệm khoa Kỹ thuật Xây Dựng, tận tình hướng dẫn em suốt thời gian thực luận văn Chính kiến thức mà em có từ truyền đạt tận tụy Thầy tài liệu chuyên ngành Thầy giúp em hòan thành luậa văn Một lần nữa, em xin bày tỏ lòng biết ơn sâu sắc thành kính với Thầy Em xin chân thành cảm ơn Thầy Phó giáo sư Tiến sỹ Chu Quốc Thắng, truyền đạt kiến thức quý báu, đặc biệt kiến thức Phương Pháp Phần Tử Hữu Hạn, định hướng Thầy việc giúp chúng em tìm hiểu phần mềm lập trình Matlab Em xin chân thành cảm ơn Thầy Tiến sỹ Đỗ Kiến Quốc, thầy Tiến Sỹ Nguyễn Văn Hiệp, thầy Giáo sư Tiến sỹ Hòang Văn Tân, cô Tiến sỹ Nguyễn Thị Hiền Lương, cô Thạc sỹ Đỗ Thị Xuân Lan,…đã truyền đạt cho chúng em nhiều kiến thức quý giá phục vụ cho công tác sau Em xin cảm ơn thầy, cô Phòng Quản Lý Sau Đại Học cung cấp cho chúng em thông tin cần thiết quan trọng suốt trình học tập nghiên cứu Xin cảm ơn người bạn chân tình giúp đỡ hòan thành luận văn Cuối cùng, Xin cảm ơn ba mẹ, anh chị em gia đình luôn quan tâm, giúp đỡ, động viên suốt trình học tập nghiên cứu Kính chúc Thầy, Cô, Ba Mẹ, Anh Chị , bạn sức khỏe hạnh phúc ! TÓM TẮT LUẬN VĂN I TÊN ĐỀ TÀI : Tự động hóa thiết kế tối ưu kết cấu thép dạng dàn phương pháp phần tử hữu hạn có kể đến ràng buộc theo tiêu chuẩn Việt Nam II TÓM TẮT : Cùng với phát triển vượt bậc ngành khoa học lý thuyết, khoa học kỹ thuật ứng dụng, đặc biệt phát triển ngành Công Nghệ Tin Học, việc áp dụng phương pháp tối ưu ngành kinh tế, kỹ thuật trở nên phổ biến Trong lónh vực thiết kế công trình, việc áp dụng phương pháp tối ưu mang lại hiệu rõ rệt giá trị vốn đầu tư thường lớn Ở đây, khuôn khổ đề tài Luận văn Thạc sỹ, em xin trình bày tóm tắt nội dung trình nghiên cứu lập chương trình máy tính cho việc tự động hóa thiết kế tối ưu kết cấu thép dạng dàn phẳng, dàn không gian có kể đến điều kiện theo TCVN thông qua việc áp dụng phương pháp PTHH để phân tích kết cấu, thuật tóan đơn hình kết hợp với vòng lặp để tìm nghiệm tối ưu Bài tóan thiết kế tối ưu kết cấu thép dạng dàn với hàm mục tiêu trọng lượng thân tòan dàn Các biến thiết kế diện tích tiết diện dàn (ở chọn lọai tiết diện thép ống) Các điều kiện ràng buộc cần thỏa mãn bao gồm : ràng buộc điều kiện bền, ràng buộc điều kiện ổn định Euler, ràng buộc điều kiện chuyển vị, ràng buộc điều kiện kiến trúc, ràng buộc điều kiện độ mãnh giới hạn, điều kiện ràng buộc khác trình thiết lập tóan tối ưu Do ràng buộc ứng suất, chuyển vị, ổn định hàm phi tuyến nên tóan thuộc dạng quy họach phi tuyến Để giải tóan này, ta sử dụng phương pháp chuỗi chương trình tuyến tính nhằm đưa tóan quy họach phi tuyến dạng tóan quy họach tuyến tính, sau sử dụng thuật tóan đơn hình đẻ giải Do đề tài nghiên cứu mang tính tự động hóa nên vấn đề lập chương trình máy tính tự động bao gồm nội dung sau : tự động hóa phân tích kết cấu PPPTHH; tự động hóa thiết lập ràng buộc, hàm mục tiêu; tự động hóa tuyến tính hóa hàm phi tuyến; tự động hóa tìm lời giải tối ưu thuất tóan đơn hình; tự động hóa xử lý kết quả; tự động hóa xuất kết nội lực, chuyển vị, phản lực gối tựa ứng với kết biến thiết kế sau Quá trình giải tóan qui họach phi tuyến dựa chuỗi chương trình tuyến tính đòi hỏi ta thực nhiều vòng lặp mà kết vòng lặp giá trị khởi đầu cho vòng lặp điều kiện hội tụ đạt tới Quá trình lập chương trình máy tính tự động hóa thiết kế tối dựa ngôn ngữ lập trình Matlab 5.3 Trong bao gồm chương trình chương trình con, tất chương trình thể dạng file.M Sau thiết lập chương trình máy tính phục vụ tóan thiết kế tối ưu, tác giả kiểm tra độ xác chương trình ví dụ số cụ thể ứng với sơ đồ kết cấu dàn phẳng, dàn không gian ứng với nhiều trường hợp khác nhau, sau tác giả đưa nhận xét kết luận Luận Văn bao gồm chương : Chương 1: Phần Tổng Quan Chương 2: Bài Toán Tối Ưu Và Các Phương Pháp Thiết Kế Tối Ưu Chương 3: Phương Pháp Quy Hoạch Tuyến Tính Và Thuật Toán Đơn Hình Chương 4: Cơ Sở Và Các Bước Phân Tích Chung Của Phương Pháp Phần Tử Hữu Hạn Sử Dụng Mô Hình Tương Tích Chương 5: Thiết Lập Chương Trình Máy Tính Cho Việc Tự Động Hóa Phân Phân Tích Kết Cấu Khung Không Gian Bằng Phương Pháp Phần Tử Hữu Hạn Chương 6: Các Yêu Cầu Chung Cho Việc Tính Toán Thiết Kế Kết Cấu Dàn Thép Chương 7: Thiết Lập Bài Toán Thiết Kế Tối Ưu Cho Kết Cấu Dàn Thép Không Gian Chương 8: p Dụng Số Vào Chương Trình Tự Động Hóa Thiết Kế Tối Ưu Kết Cấu Chương 9: Kết Luận Và Kiến Nghị ….………….…………… ….Thuyết Minh Luận Văn Thạc sỹ Ngành XDDD_CN…… ….………………… …… MỤC LỤC Nhiệm vụ luận văn Thạc Só Lời Cám Ơn Tóm Tắt Luận Văn Chương 1: PHẦN TỔNG QUAN………………………………………………………………………………………….1 1.1 Sơ lược lịch sư ûquá trình phát triển thiết kế tối ưu……………………………………….1 1.2 Ứng dụng phương pháp tối ưu ngành kinh tế kỹ thuật………………………….2 1.3 Ứng dụng thiết kế tối ưu kết cấu giới………………………………………………….2 1.3.1 Lónh vực tính toán kết cấu thép……………………………………………………….2 1.3.2 Lónh vực tính toán kết cấu bê tông cốt thép………………………………………3 1.4 Những đặc điểm thiết kế kết cấu thép nước ta ………………………3 Chương 2: BÀI TOÁN TỐI ƯU VÀ CÁC PHƯƠNG PHÁP THIẾT KẾ TỐI ƯU…………………………………4 2.1 Bài toán tối ưu…………………………………………………………………………………………….4 2.1.1 Một số thuật ngữ……………………………………………………………………………4 2.1.2 Phát biểu toán học toán tối ưu…………………………………………….5 2.1.3 Quá trình lập tối ưu……………………………………………………………………… 2.1.4 Sự tồn tính nghiệm tối ưu………………………………… 2.1.5 Thuận lợi việc sử dụng tối ưu hóa phương pháp số …………….7 2.1.6 Giới hạn tính tối ưu phương pháp số……………………………….…7 2.2 Các phương pháp thiết kế tối ưu……………………………………………………………… ….8 2.2.1- Phương pháp cổ điển …………………………………………………………….………8 2.2.2- Phương pháp tiêu chuẩn tối ………………………………………………………… 2.2.3- Phương pháp quy họach tóan học…………………………………………………….8 2.3- Chọn phương pháp tính toán tối ưu ……………………………………………………………….8 2.3.1- Phương pháp quy hoạch toán học…………………………………………………….8 2.3.2 Các phương pháp giải toán quy hoạch phi tuyến có ràng buộc…….11 2.3.3 Phương pháp dựa chuỗi chương trình tuyến tính …………………11 2.3.4 Giới thiệu phần mềm Matlab 5.3 ………………………………………………… 12 Chương 3: PHƯƠNG PHÁP QUY HOẠCH TUYẾN TÍNH VÀ THUẬT TOÁN ĐƠN HÌNH……………….13 3.1 thuật toán đơn hình từ vựng (ý nghóa giải tích)………………………………….…….13 3.2 Tìm từ vựng xuất phát ……………………………………………………………………………… 15 3.3- Trường hợp mục tiêu không giới nội ………………………………………………………… 15 3.4- Sự nghiệm tối ưu ……………………………………………………………………… 15 3.5- Dạng bảng phương pháp đơn hình ……………………………………………………… 16 3.6 Hiện tượng suy biến ……………………………………………………………………………………16 3.7 Ý nghóa hình học quy hoạch tuyến tính phương pháp đơn hình ……………………16 3.7.1 Điểm cực biên nghiệm sở …………………………………………………….16 3.7.2 Trường hợp tập ràng buộc có dạng tắc……………………………………18 3.7.3 Sự tồn đỉnh đỉnh tối ưu ……………………………………………………… 18 3.7.4 Hiện tượng suy biến …………………………………………………………………… 19 3.7.5 Giải thích hình học phương pháp đơn hình ………………………………………19 Nguyễn HữuT hịnh…………………………………… …TpHCM2003…………………… ………………………… ….………….…………… ….Thuyết Minh Luận Văn Thạc sỹ Ngành XDDD_CN…… ….………………… …… Chương 4: CƠ SỞ VÀ CÁC BƯỚC PHÂN TÍCH CHUNG CỦA PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN SỬ DỤNG MÔ HÌNH TƯƠNG TÍCH ……………………………………………………………………………………….19 4.1 Cơ sở nghiên cứu …………………………………………………………………………………… 19 4.1.1 Nguyên lý công ảo (nguyên lý công khả dó) …………………………………….19 4.1.2 Nguyên lý toàn phần dừng ……………………………………………….21 4.2 Phương pháp phần tử hữu hạn mô hình phân tích kết cấu…………… ………22 4.3 Trình tự phân tích toán kết cấu theo phương pháp phần tử hữu hạn ……………27 4.4 Tính gradient hay độ nhạy (sensitivity) ……………………………………….…………………28 Chương 5: THIẾT LẬP CHƯƠNG TRÌNH MÁY TÍNH CHO VIỆC TỰ ĐỘNG HÓA PHÂN PHÂN TÍCH KẾT CẤU KHUNG KHÔNG GIAN BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN……………………………29 5.1 số lý thuyết việc phân tích kết cấu khung không gian ………….29 5.1.1 ký hiệu ………………………………………………………………………………….29 5.1.2 Xác định ma trận độ cứng phần tử hệ toạ độ địa phương ………….30 5.1.3 Xác định ma trận độ cứng phần tử hệ toạ độ tổng thể …………… 30 5.1.4 Xác định ma trận chuyển đổi toạ độ [Gel] ………………………………… ….31 5.1.5 Xác định ma trận độ cứng tổng thể kết cấu [K] ………………………….31 5.1.6 Xác định véc tơ tải tổng thể {P} ………………………………………………………31 5.1.7 Xác định véc tơ chuyển vị nút tổng thể {U}……………………………………….31 5.2 Ma trận chuyển đổi tọa độ phần tử …………………………………………………………… 31 5.2.1 Góc xoay tọa độ không gian ba chiều …………………………………….32 5.2.2 Ma trận chuyển đổi tọa độ …………………………………………………………….33 5.3 Áp đặt điều kiện biên (Quá trình rút gọn ma trận) …………………………………………38 5.4 Lập chương trình máy tính tự động hóa việc phân tích kết cấu khung không gian (frame3d.m) ……………………………………………………………………………………………………………….38 5.4.1 Sau sơ đồ khối chương trình …………………………………………………38 5.4.2 Các file.m sử dụng việc tự động hóa phân tích kết cấu khung không gian …………………………………………………………………………………………………………………40 5.5 p dụng chương trình FRAME3D.m ………………………………………………………………46 5.5.1 gIới thiệu chương trình FRAME3D.m ………………………………………………46 5.5.2 Sử dụng chương trình FRAME3D ……………………………………………………46 Chương 6: CÁC YÊU CẦU CHUNG CHO VIỆC TÍNH TOÁN THIẾT KẾ KẾT CẤU DÀN THÉP ………48 6.1 Các điều kiện tính toán thiết kế kết cấu dàn thép …………………………………….48 6.1.1 Điều kiện bền ………………………………………………………………………………48 6.1.2 Điều kiện ổn định …………………………………………………………………………48 6.1.3 Điều kiện biến dạng (chuyển vị) …………………………………………………….49 6.1.4 Điều kiện độ mảnh giới hạn dàn λmax ………………………….50 Chương 7: THIẾT LẬP BÀI TOÁN THIẾT KẾ TỐI ƯU CHO KẾT CẤU DÀN THÉP KHÔNG GIAN ….50 7.1 Chọn hệ thống đơn vị …………………………………………………………………………………50 7.2 Chọn vật liệu thép …………………………………………………………………………………….50 7.3 Chọn hình dáng tiết diện ngang dàn ………………………………………………… 50 7.4 Xác định biến thiết kế hàm mục tiêu …………………………………………………50 Nguyễn HữuT hịnh…………………………………… …TpHCM2003…………………… ………………………… ….………….…………… ….Thuyết Minh Luận Văn Thạc sỹ Ngành XDDD_CN…… ….………………… …… 7.4.1 Xác định biến thiết kế ………………………………………………………… 51 7.4.2 Xác định hàm mục tiêu ……………………………………………………………… 51 7.5 Xác định ràng buộc cho toán thiết kế tối ưu …………………………………….51 7.5.1 Ràng buộc ứng suất ……………………………………………………………… 51 7.5.2.Ràng buộc ổn định theo công thức Euler …………………………………….51 7.5.3 Ràng buộc điều kiện để áp dụng công thức Euler ………………………52 7.5.4 Ràng buộc điều kiện chuyển vị …………………………………………………52 7.5.5 Ràng buộc độ mảnh giới hạn dàn …………………… 53 7.5.6 Điều kiện giới hạn dịch chuyển (movelimit) phương pháp giải toán qui hoạch phi tuyến dựa chuỗi chương trình tuyến tính …………………53 7.6 Chọn giá trị biến thiết kế ban đầu { X io } …………………………………………………….54 7.7 xác định hệ số {a ji }, {b j }, {ci } chuỗi hàm mục tiêu, ràng buộc từ việc tuyến tính hóa quanh điểm { X io } theo công thức Taylor ……………………………………… 54 7.8 Lập chương trình máy tính phục vụ cho trình tự động hóa thiết kế tối ưu ………………………………………………………………………………………………………………………………….55 7.9 Sơ đồ khối cho chương trình tự động hoá thiết kế tối ưu (MAINPROGRAMME.M) ……………………………………………………………………………………………… 56 Chương 8: ÁP DỤNG SỐ VÀO CHƯƠNG TRÌNH TỰ ĐỘNG HÓA THIẾT KẾ TỐI ƯU KẾT CẤU THÉP DẠNG DÀN …………………………………………………………………………………………………………………58 8.1 Bài tóan 1: Dàn Phẳng Thanh……………………………………………………………………58 8.2 Bài tóan : Dàn Không Gian Thanh ………………………………………………………….66 8.3.Bài Tóan : Dàn Phẳng Thanh ………………………………………………………………….79 8.4 Bài Tóan : Dàn Phẳng Thanh ………………………………………………………………….82 8.5 Bài Tóan : Dàn Không Gian 24 Thanh ……………………………………………………… 85 8.6 Bài Tóan : Dàn Không Gian 27 Thanh ……………………………………………………… 90 Chương 9: KẾT LUẬN VÀ KIẾN NGHỊ …………………………………………………………………………… 98 9.1 kết luận ………………………………………………………………………………………………… 98 9.2 kiến nghị ………………………………………………………………………………………………….98 Tài liệu tham khảo …………………………………………………………………………………………………….99 Nguyễn HữuT hịnh…………………………………… …TpHCM2003…………………… ………………………… ….………….…………… ….Thuyết Minh Luận Văn Thạc sỹ Ngành XDDD_CN…… ….………………… …… CHƯƠNG 1:PHẦN TỔNG QUAN 1.1 Sơ lược lịch sư ûquá trình phát triển thiết kế tối ưu: [11] Việc phân tích kết cấu tối ưu kết cấu nhiều tình xem hai khía cạnh mục đích Vì lý mà chúng phát triển theo thời gian Trong nhiều thời điểm khác nhau, lý thuyết phân tích kết cấu thường quan tâm nhiều Do đó, nhiều đóng góp quan trọng lónh vực lý thuyết phân tích kết cấu theo thực tiễn tạo nhiều khám phá quan trọng lónh vực tối ưu kết cấu Một cách thích hợp để nói công việc Galileo, người tiên phong việc cố gắng mô tả trạng thái ứng suất kết cấu dầm tác dụng tải trọng gây uốn, cho dù kết không xác Ôâng tham gia vào vấn đề tìm kiếm hình dáng tối ưu cho kết cấu dầm sâu có tiết diện thay đổi (Variable depth beam) Công việc nghiên cứu ban đầu phát triển Johan Bernoulli Những nghiên cứu ông cho ứng suất phẳng tiết diện ngang dầm tuân theo định luật Hook từ nảy sinh ứng suất dọc trục tương tự lý thuyết Galileo Từ quan điểm phân tích kết cấu dầm, khám phá quan trọng có tính chất phổ biến Parent vào năm 1708, 1710 Theo đó, ông thừa nhận tồn ứng suất không dọc trục đưa khái niệm trục trung hòa (Neutral Axis) Bằng cách sử dụng lý thuyết này, ông đề xuất vấn đề thiết kế tối ưu cho dầm có ứng suất phân bố chịu tác dụng tải trọng di chuyển Sau đó, Lagrange, người công việc nghiên cứu trước Parent, ôâng thừa nhận lại khái niệm trục trung hòa.Và vào năm 1770-1773, ông tiến hành thiết kế cột tác dụng tảt trọng dọc trục ông phát hình dáng đường cong đàn hồi (the shape of the elastic curve) Từ nghiên cứu này, ông tiến hành thiết kế trọng lượng nhỏ cột có hình dạng đối xứng trục, kết ông thiết kế tối ưu để dẫn đến tiết diện không đổi cột hình trụ tròn không trọng lượng thân giảm dần Tuy nhhiên, vấn đề giải đắn Clausen vào năm 1851 Nghiên cứu lónh vực phát triển bình thường, đến cuối kỷ 19 có đóng góp quan trọng Levy Maxwell Sự đóng góp Levy liên quan đến vấn đề thiết kế kết cấu dàn, vòm với ứng suất phân bố Ông ta kết cấu dạng dàn tam giác Lattice tác dụng tải trọng không đổi phải tónh định để đạt tối ưu Một vài năm trước năm 1869, Maxwell đưa quan điểm lý thuyết kết cấu tối ưu Vấn đề liên quan đến kết cấu dạng dàn có liên kết khớp xoay Pin-Jionted Frameworks chịu tác dụng tải trọng đơn ông đưa định đề sau: ∑ltiσtI ∑lcjσcj = constant, lti, lcj chiều dài phần thử chịu kéo, nén tương ứng ; σtI, σcj ứng suất cho phép phần thử chịu kéo, nén tương ứng Vào năm 1904, A.G.M Michell dựa quan điểm Maxwell để phát triển lý thuyết trình bày thiết kế tối ưu kết cấu dạng dàn có liên kết khớp xoay Pin-Jionted Frameworks chịu ràng buộc ứng suất với trường hợp tải trọng đơn Năm 1943, Smith Cox đưa nguyên lý thiết kế trọng lượng nhỏ Nguyễn Hữu Thịnh…………………………………… …TpHCM2003…………………… ………………………… ….………….…………… ….Thuyết Minh Luận Văn Thạc sỹ Ngành XDDD_CN…… ….………………… …… cho kết cấu chịu nén sau ứng dụng cho kết cấu panel cứng (stiffened panels) Năm 1944, Zahorki nghiên cứu thiết kế tối ưu trọng lượng cho kết cấu dạng Stringer panels Sau Leggett, Hopkins Wittrick tham giam vào thiết kế kết cấu dạng Sandwich panels chịu ràng buộc ổn định Từ hình thành lý thuyết thiết kế tối ưu trọng lượng chịu ràng buộc ổn định Sau này, với phát triển vượt bậc khoa học lý thuyết: toán học, vật lý học,…,khoa học thực nghiệm, với đời máy tính kỹ thuật số, từ phương pháp tối ưu thuật toán tương ứng đời Chính tạo đà cho thiết kế tối ưu phát triển nhanh chóng 1.2 Ứng dụng phương pháp tối ưu ngành kinh tế kỹ thuật: [9], [11], [12] Trong năm gần đây, phương pháp tối ưu hóa ngày áp dụng sâu rộng hiệu vào nghành kinh tế, kỹ thuật, công nghệ thông tin vào ngành khoa học công nghệ khác - Lựa chọn vật liệu để luyện gang, thép hợp kim… ngành luyện kim - Giải vấn đề điện khí hóa vùng có phân bố nguồn lượng khác ngành lượng - Chọn lưu lượng dầu khai thác giếng công nghiệp dầu mỏ - Phân bố hợp lý diện tích trồng trọt ngành nông nghiệp 1.3 ỨNG DỤNG THIẾT KẾ TỐI ƯU KẾT CẤU TRÊN THẾ GIỚI: [9],[10], [11], [12] Trong thiết kế kết cấu công trình nhiều nơi giới người ta áp dụng phương pháp tối ưu hoá vào toán thiết kế toán như: kết cấu thép, kết cấu bê tông cốt thép Quá trình tính toán phức tạp khối lượng tính toán lớn Tuy nhiên, với phát triển máy tính phần mềm lập trình ngày việc giải vấn đề nêu có nhiều thuận lợi 1.3.1 Lónh vực tính toán kết cấu thép: -Thiết kế tối ưu kết cấu dạng dàn có sơ đồ hình học cho trước (Optimal Design of Trusses of Specified) với biến thiết kế Xi diện tích tiết diện ngang dàn, hàm mục tiêu Z trọng lượng daøn n Z = ρ ∑ Li X i i =1 Trong : ρ : trọng lượng riêng dàn Li : chiều dài dàn Xi : tiết diện dàn Trong trình tính toán để Z đạt giá trị nhỏ nhất, biến Xi phải xác định cho kết cấu dàn thỏa mãn yêu cầu ràng buộc gồm : ứng suất kéo nén (compressive stress, tensile stress), chuyển vị (displacement), độ ổn định Euler (Euler Nguyễn Hữu Thịnh…………………………………… …TpHCM2003…………………… ………………………… ….………….…………… ….Thuyết Minh Luận Văn Thạc sỹ Ngành XDDD_CN…… ….………………… …… buckling)… - Thiết kế tối ưu trọng lượng kết cấu dạng dàn vớiù sơ đồ hình học không cho trước (Optimal Design of Trusses including Configuration Variables), biến thiết kế bao gồm: tiết diện ngang dàn, biến liên quan đến hình dạng hình học, vị tướng kết cấu dàn - Thiết kế tối ưu dầm hộp thép, dầm tổ hợp thép với hàm mục tiêu trọng lượng dầm với biến kích thước cánh bụng dầm Các ràng buộc cần thỏa mãn ứng suất độ ổn định - Thiết kế tối ưu hình dáng phận khí liên quan đến khí động lực học (Aerodynamic Shape Optimization): thiết kế tối ưu kết cấu thép cho đuôi trực thăng, cánh máy bay, vỏ máy bay… - Thiết kế tối ưu dàn thép anten 55 phân tử, dàn thép 200 phân tử 1.3.2 Lónh vực tính toán kết cấu bê tông cốt thép: - Ứng dụng thiết kế tối ưu kết cấu BTCT thường (Reinforced Concrete Structure), BTCT dự ứng lực (Prestressed Concrete Structure) - Thiết kế tối ưu cho dầm, khung BTCT thường - Thiết kế tối ưu cho sàn BTCT thường - Thiết kế tối ưu cho dầm BTCT dự ứng lực - Thiết kế tối ưu cho sàn BTCT dự ứng lực Ở nước ta, phương pháp tối ưu hoá nghiên cứu áp dụng nghành xây dựng Tuy nhiên, công tác thiết kế kết cấu công trình phương pháp tối ưu hoá chưa nghiên cứu sâu để áp dụng cách thức phổ biến 1.4 Những đặc điểm thiết kế kết cấu thép nước ta : Trong thiết kế kết cấu thép nước ta đặc biệt kết cấu khung,dàn trình tự thiết kế thường theo bước sau : Chọn trước kích thước tiết diện cấu kiện theo kinh nghiệm có kết hợp với số điều kiện cấu tạo yêu cầu kiến trúc Tính toán nội lực giai đoạn đàn hồi (Dùng phần mềm phổ biến phân tích kết cấu phương pháp PTHH giải trực tiếp phương pháp học kết cấu) Thiết kế lại tiết diện ngang thép theo nội lực tính toán tiết diện chọn theo trạng thái giới hạn Kiểm tra điều kiện cấu tạo, giới hạn chuẩn bị, độ võng,v.v… Khi không thỏa điều kiện yêu cầu chọn lại kích thước tiết diện tính toán lại theo trình tự * Vấn đề cần nghiên cứu: Với bước thiết kế kết cấu thép trình bày đây, ta nhận thấy có nhiều phương án thiết kế thỏa mãn yêu cầu làm việc bền vững công trình người thiết kế cần tìm phương án có phương án thiết kế Vì vậy, người thiết kế không khẳng định phương án có tốt phương án khác hay không Trong số nơi giới người ta áp dụng phương pháp tối ưu Nguyễn Hữu Thịnh…………………………………… …TpHCM2003…………………… ………………………… … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… 14 15 16 17 18 19 20 21 22 23 24 25 26 27 11 10 eps eps eps 12 10 eps eps eps 10 10 eps eps eps 11 10 eps eps eps 12 10 eps eps eps 10 eps eps eps 10 eps eps eps 10 eps eps eps 10 eps eps eps 10 eps eps eps 10 eps eps eps 10 11 10 eps eps eps 11 12 10 eps eps eps 10 12 10 eps eps eps 0]; SupJnt = [1 1 1 1 1]; SupJnt4 = [1 1 0 1 0 1 0 0]; lamdamax=150;%Do manh cho phep cua theo phuong bat ky rad_scale=0.9;%Ty so ban kinh va ngoai cua tiet dien ong L0=zeros(nemax,1); var_cp=zeros(nemax,1); for i=1:nemax x_coords=[JntCoors(MbData(i,2),2),JntCoors(MbData(i,3),2)]; y_coords=[JntCoors(MbData(i,2),3),JntCoors(MbData(i,3),3)]; z_coords=[JntCoors(MbData(i,2),4),JntCoors(MbData(i,3),4)]; ex_bar = [x_coords(2)-x_coords(1) y_coords(2)-y_coords(1) z_coords(2)-z_coords(1)]; L0(i,1)= sqrt(ex_bar*ex_bar')*0.9; %Chieu dai tinh toan cua end for i=1:nemax var_cp(i,1)=(4*L0(i,1)^2*(1-rad_scale^2))/(lamdamax^2*(1+rad_scale^2)); end JntLoad JntLoad4 = = [9 5000 -3000 0]; [9 5000 -3000 0 0]; Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 124 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… XUẤT KẾT QUAÛ : Flect Shape of 27 Space Truss Structure 900 26 25 800 27 14 700 17 16 Y 600 15 18 13 23 500 22 24 400 11 10 300 20 19 200 21 100 12 200 200 Z X XX = 4.271787845303866e+000 4.271787845303866e+000 9.005132887309401e+000 8.543575690607732e+000 1.281536353591160e+001 8.543575690607732e+000 4.271787845303866e+000 4.271787845303866e+000 6.244482379896856e+000 8.543575690607732e+000 1.281536353591160e+001 8.543575690607732e+000 4.271787845303866e+000 4.271787845303866e+000 4.271787845303866e+000 8.543575690607732e+000 1.281536353591160e+001 8.543575690607732e+000 4.271787845303866e+000 8.543575690607732e+000 4.936696782315726e+000 4.271787845303866e+000 8.543575690607732e+000 4.271787845303866e+000 4.271787845303866e+000 Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 125 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… 8.543575690607732e+000 4.271787845303866e+000 vonglap = 15 chenhlech = Columns through 1.378890190569208e+000 1.105117740995423e+000 Columns through 0 0 8.881784197001252e-016 0 0 0 8.881784197001252e-016 0 0 Columns through Columns 10 through 12 Columns 13 through 15 Columns 16 through 18 Columns 19 through 21 Columns 22 through 24 Columns 25 through 27 optimal_rangbuoc = Columns through 12 63 75 109 110 112 113 114 115 116 118 119 120 Columns 13 through 24 121 122 123 124 125 126 127 128 130 131 132 133 Columns 25 through 28 134 135 199 201 Nguyeãn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 126 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… object_value = 7.403409419805855e+004 optimal_rad = 2.675858926489799e+000 2.675858926489799e+000 3.885111063713306e+000 3.784235984838984e+000 4.634723614567045e+000 3.784235984838984e+000 2.675858926489799e+000 2.675858926489799e+000 3.235240702033668e+000 3.784235984838984e+000 4.634723614567045e+000 3.784235984838984e+000 2.675858926489799e+000 2.675858926489799e+000 2.675858926489799e+000 3.784235984838984e+000 4.634723614567045e+000 3.784235984838984e+000 2.675858926489799e+000 3.784235984838984e+000 2.876580963473171e+000 2.675858926489799e+000 3.784235984838984e+000 2.675858926489799e+000 2.675858926489799e+000 3.784235984838984e+000 2.675858926489799e+000 vonglap = 15 chenhlechmax = 1.000000000000000e-006 XX = 4.271787845303866e+000 4.271787845303866e+000 9.005132887309401e+000 8.543575690607732e+000 1.281536353591160e+001 8.543575690607732e+000 4.271787845303866e+000 4.271787845303866e+000 6.244482379896856e+000 8.543575690607732e+000 1.281536353591160e+001 8.543575690607732e+000 4.271787845303866e+000 4.271787845303866e+000 4.271787845303866e+000 8.543575690607732e+000 1.281536353591160e+001 Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 127 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… 8.543575690607732e+000 4.271787845303866e+000 8.543575690607732e+000 4.936696782315726e+000 4.271787845303866e+000 8.543575690607732e+000 4.271787845303866e+000 4.271787845303866e+000 8.543575690607732e+000 4.271787845303866e+000 timecpustop = 4.543001600000000e+004 » ******************************************************************************* * Space Frame Analysis Program FRAME3D * * written in MATLAB by Nguyen Huu Thinh * * Thanh Pho Ho Chi Minh thang 03 nam 2003 * ******************************************************************************* Problem title : Analysis of a Space Truss Structure Units in daN & cm GENERAL STRUCTURE DATA Number of joints : 12 Number of members : 27 Number of materials : Number of support joints : Number of loaded joints : JOINT COORDINATES Joint No X-coord 0.000 0.000 300.000 0.000 0.000 300.000 0.000 0.000 300.000 10 0.000 11 0.000 12 300.000 MATERIAL DATA Material no MEMBER DATA Member e1 1 2 3 4 5 10 11 12 Y-coord 0.000 0.000 0.000 300.000 300.000 300.000 600.000 600.000 600.000 900.000 900.000 900.000 Z-coord 0.000 300.000 0.000 0.000 300.000 0.000 0.000 300.000 0.000 0.000 300.000 0.000 Modulus of elasticity 2.00e+006 e2 6 9 Area 4.27e+000 4.27e+000 9.01e+000 8.54e+000 1.28e+001 8.54e+000 4.27e+000 4.27e+000 6.24e+000 8.54e+000 1.28e+001 8.54e+000 J 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 Shear Modulus 1.20e+006 Iy 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 Iz 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 Mtl 1 1 1 1 1 1 Angle 0 0 0 0 0 0 Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 128 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 10 11 7 10 11 10 10 11 12 12 6 9 11 12 12 4.27e+000 4.27e+000 4.27e+000 8.54e+000 1.28e+001 8.54e+000 4.27e+000 8.54e+000 4.94e+000 4.27e+000 8.54e+000 4.27e+000 4.27e+000 8.54e+000 4.27e+000 SUPPORT JOINT RESTRAINTS (Fixed:1 Released:0) Joint no X Y Z XX 1 1 1 1 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 YY 0 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 2.22e-016 1 1 1 1 1 1 1 0 0 0 0 0 0 0 ZZ 0 JOINT LOADS Joint X Y 5.000e+003 -3.000e+003 Z 0.000e+000 XX 0.000e+000 YY 0.000e+000 ZZ 0.000e+000 DISPLACEMENTS Joint X 0.000e+000 0.000e+000 0.000e+000 6.168e-001 1.756e-001 4.648e-001 1.449e+000 5.677e-001 1.449e+000 10 2.034e+000 11 1.152e+000 12 2.034e+000 Z 0.000e+000 0.000e+000 0.000e+000 1.756e-001 1.756e-001 4.648e-001 3.511e-001 3.511e-001 1.233e+000 5.267e-001 5.267e-001 1.408e+000 XX 5.241e-004 5.425e-004 1.492e-003 7.163e-004 6.386e-004 1.818e-003 6.031e-004 6.348e-004 1.504e-003 5.676e-004 5.561e-004 3.285e-004 YY -6.290e-004 -1.032e-003 -4.452e-004 -1.373e-003 -1.443e-003 -8.280e-004 -2.899e-003 -2.718e-003 -2.907e-003 -2.936e-003 -3.007e-003 -2.961e-003 ZZ -1.526e-003 -3.344e-004 -1.349e-003 -2.010e-003 -1.247e-003 -1.918e-003 -2.191e-003 -1.860e-003 -2.574e-003 -1.955e-003 -1.953e-003 -1.791e-003 xx yy zz Y 0.000e+000 0.000e+000 0.000e+000 1.756e-001 -1.219e-015 -2.165e-001 1.756e-001 -1.667e-015 -4.087e-001 1.756e-001 -1.655e-015 -4.087e-001 FORCES Member x y z 1 -5.000e+003 1.704e-017 2.071e-018 5.000e+003 -1.704e-017 -2.071e-018 6.612e-016 -2.610e-017 -6.612e-016 -5.952e-016 3.272e-015 1.839e-015 2 3.473e-011 -1.215e-017 3.129e-019 -3.473e-011 1.215e-017 -3.129e-019 3.651e-016 9.533e-017 -3.651e-016 -1.892e-016 -4.728e-016 -3.173e-015 3 1.300e+004 -4.977e-018 6.248e-018 -1.300e+004 4.977e-018 -6.248e-018 3.401e-016 -4.558e-016 -3.401e-016 -1.419e-015 9.463e-017 -1.588e-015 4 2.819e-011 -1.307e-018 6.679e-018 -2.819e-011 1.307e-018 -6.679e-018 -5.924e-016 -2.910e-015 5.924e-016 7.618e-017 -4.593e-016 -9.533e-017 5 -3.638e-012 -3.245e-018 5.601e-018 3.638e-012 3.245e-018 -5.601e-018 1.218e-017 -2.485e-015 -1.218e-017 -4.251e-016 -1.297e-015 -3.892e-016 6 -7.071e+003 -1.736e-017 7.325e-018 7.071e+003 1.736e-017 -7.325e-018 -4.860e-016 -4.491e-016 4.860e-016 -2.658e-015 -3.272e-015 -4.092e-015 Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 129 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… 7 -3.274e-011 3.995e-017 4.410e-018 3.274e-011 -3.995e-017 -4.410e-018 1.355e-015 -8.290e-016 -1.355e-015 -4.939e-016 6.260e-015 5.724e-015 8 1.276e-011 -1.459e-017 3.047e-018 -1.276e-011 1.459e-017 -3.047e-018 1.133e-015 -4.626e-016 -1.133e-015 -4.514e-016 -1.280e-015 -3.097e-015 9 8.000e+003 6.133e-017 -5.323e-017 -8.000e+003 -6.133e-017 5.323e-017 1.847e-015 -1.847e-015 7.521e-015 8.448e-015 1.017e-014 8.227e-015 10 2.910e-011 -1.054e-018 -7.458e-018 -2.910e-011 1.054e-018 7.458e-018 2.271e-016 -1.946e-016 -2.271e-016 3.359e-015 -1.865e-016 -2.607e-016 11 -3.638e-012 8.241e-019 3.182e-018 3.638e-012 -8.241e-019 -3.182e-018 1.924e-016 -2.579e-015 -1.924e-016 9.249e-016 -4.660e-016 8.942e-016 12 -7.071e+003 2.051e-018 -5.724e-018 7.071e+003 -2.051e-018 5.724e-018 3.313e-016 2.933e-015 -3.313e-016 -5.042e-016 1.026e-015 -1.554e-016 13 10 -9.095e-012 -7.437e-018 6.943e-021 9.095e-012 7.437e-018 -6.943e-021 3.349e-017 -5.360e-017 -3.349e-017 5.151e-017 -1.465e-015 -7.665e-016 14 11 -3.490e-013 2.416e-018 6.041e-019 3.490e-013 -2.416e-018 -6.041e-019 2.565e-016 -2.072e-016 -2.565e-016 2.593e-017 5.002e-016 2.245e-016 15 12 1.819e-012 -1.391e-017 1.961e-017 -1.819e-012 1.391e-017 -1.961e-017 4.774e-017 -4.681e-015 -4.774e-017 -1.201e-015 -3.246e-015 -9.275e-016 16 10 1.273e-011 -4.727e-019 -3.168e-018 -1.273e-011 4.727e-019 3.168e-018 5.468e-017 -5.468e-017 4.403e-016 9.040e-016 -2.994e-017 -1.706e-016 17 11 -1.455e-011 1.796e-018 -5.837e-018 1.455e-011 -1.796e-018 5.837e-018 -4.350e-016 4.350e-016 8.994e-016 2.133e-015 2.689e-016 6.645e-016 18 12 -1.455e-011 -1.286e-018 -2.687e-018 1.455e-011 1.286e-018 2.687e-018 1.496e-016 -1.496e-016 4.129e-016 7.271e-016 -6.916e-016 1.462e-016 19 -2.001e-011 -5.460e-018 -3.698e-018 2.001e-011 5.460e-018 3.698e-018 -6.781e-016 6.781e-016 6.574e-016 4.519e-016 -9.341e-016 -7.039e-016 20 0.000e+000 7.699e-018 5.069e-018 0.000e+000 -7.699e-018 -5.069e-018 -8.216e-016 -1.719e-015 8.216e-016 -4.318e-016 1.257e-015 2.009e-015 21 5.000e+003 -3.889e-017 8.084e-018 -5.000e+003 3.889e-017 -8.084e-018 -9.781e-016 -2.020e-015 9.781e-016 -4.052e-016 -5.971e-015 -5.697e-015 22 -1.637e-011 -1.997e-018 -7.719e-018 1.637e-011 1.997e-018 7.719e-018 -2.939e-016 2.939e-016 8.907e-016 1.425e-015 -2.526e-016 -3.464e-016 23 1.455e-011 4.493e-018 -3.732e-018 -1.455e-011 -4.493e-018 3.732e-018 -7.032e-016 7.032e-016 9.897e-016 5.937e-016 8.381e-016 1.068e-015 Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 130 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… 24 1.455e-011 -2.576e-017 -2.120e-018 -1.455e-011 2.576e-017 2.120e-018 -8.003e-016 8.003e-016 3.308e-016 3.053e-016 -3.297e-015 -4.431e-015 25 10 11 -1.091e-011 1.385e-018 1.948e-018 1.091e-011 -1.385e-018 -1.948e-018 -1.714e-018 -1.876e-016 1.714e-018 -3.968e-016 1.908e-016 2.249e-016 26 11 12 0.000e+000 -1.413e-018 1.339e-018 0.000e+000 1.413e-018 -1.339e-018 1.732e-016 -3.322e-016 -1.732e-016 -2.360e-016 -2.515e-016 -3.479e-016 27 10 12 0.000e+000 4.410e-018 5.871e-019 0.000e+000 -4.410e-018 -5.871e-019 2.123e-016 -5.154e-017 -2.123e-016 -1.246e-016 4.182e-016 9.049e-016 REACTIONS Joint X Y Z -5.000e+003 -1.000e+004 9.396e-018 1.883e-017 5.474e-011 -2.001e-011 1.819e-012 1.300e+004 -1.819e-012 XX YY 0.000e+000 4.930e-031 1.972e-031 -5.916e-031 -7.889e-031 3.944e-031 ZZ -1.578e-030 -7.889e-031 -2.761e-030 PHẦN : CÁCH SỬ DỤNG CHƯƠNG TRÌNH TỰ ĐỘNG HÓA THIẾT KẾ TỐI ƯU KẾT CẤU THÉP DẠNG DÀN (MAINPROGRAMME.M) Bước1 - : Nhập liệu tóan (Opt_Data.M): Nhập tựa đề cho tóan : heading Nhập đơn vị tính : units Nhập kích thước tóan : StrSize Nhập tọa độ nút : JntCoors Nhập thông số vật liệu : MtlData Nhập cường độ kéo nén cho phép thép : R_cp Nhập hệ số liên quan đến ràng buộc Movelimit : beta1 Nhập giá trị chuyển vị nút cho phép : max_displacement Nhập biến thiết kế diện tích tiết diện dạng symbolics : de_var Nhập giá trị diện tích tiết diện : X01 Nhập liệu phần tử cho tóan dàn : MbData Nhập liệu phần tử cho việc phân tích nội lực từ kết tối ưu (Opt_Frame3d_Anls.M) : MbData4 Nhập hệ số điều kiện làm việc : heso_lamviec Nhập liên kết gối tựa cho tóan dàn : SupJnt Nhập liên kết gối tựa cho việc dàn cho tóan phân tích nội lực từ kết tối ưu (Opt_Frame3d_Anls.M) : SupJnt4 Nhập giá trị độ mảnh cho phép : lamdamax Nhập tỷ số bán kính bán kính ngòai : rad_scale Nhập tải trọng cho tóan dàn :JntLoad Nhập tải trọng cho việc dàn cho tóan phân tích nội lực từ kết tối ưu (Opt_Frame3d_Anls.M) : JntLoad4 Sau lưu liệu thành tên file : Opt_Data.M Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 131 … ……….…………………………….……….….Phụ Lục Luận Văn Thạc sỹ Ngành XDDD_CN………… ……… ….………………… …… Ghi Chú : Khi dùng chương trình Frame3d.M , Opt_Frame3d_Anls.M để phân tích cho tóan dàn không gian, dàn phẳng ta phải gán giá trị mô men quán tính chống uốn, mô men quán tính chống xoắn giá trị eps tức xấp xỉ (eps : giá trị vô bé tức xấp xỉ Matlab5.3), đồng thời khai báo liên gối tựa, véc tơ tải cho phù hợp tóan cụ thể Bước2 : Thao tác - Vào hình làm việc Matlab gõ : Mainprogramme, nhấn enter - Sau thực thao tác theo yêu cầu chương trình hình liên quan đến việc xuất kết dạng file.text, file hình ảnh biến dạng kết cấu ứng với kết tối ưu Nguyễn Hữu Thịnh………………………………………………… …TpHCM2003……………………………………………………………… 132 SƠ ĐỒ KHỐI NHẬP DỮ LIỆU BÀI TOÁN THIẾT KẾ TỐI ƯU Start Chọn sơ đồ kết cấu dàn cho việc thiết kế tối ưu Đơn vị tính (units) Chuyển vị cho phép (Max_displ acement) Kích thước toán (Strsize) Nhập liệu phần tử (MbData) (MbData4) Tọa độ nút kết cấu (Jntcoors) Hệ số làm việc (heso_lam viec) Điều kiện biên (SupJnt) (SupJnt4) Đặc trưng vật liệu thép (MtlData) Độ mảnh giới hạn (lamda_max) Cường độ tính toán thép (R) Tỷ số bán kính tiết diện ống (rad_scale) Tải trọng nút (Jntload) (Jntload4) Hệ số Movelimit (beta1) Giá trị biến thiết kế nhập vào (X01) SƠ ĐỒ KHỐI THIẾT KẾ SƠ BỘ ĐỂ CHỌN GIÁ TRỊ BIẾN THIẾT KẾ BAN ĐẦU Start Nhập liệu toán thiết kế tối ưu Opt_Data.m Phân tích toán ứng với giá trị biến thiết kế Nhập vào X01 X i0 = σ i01 01 X i = A1 [σ ] i = : ne max L20i (1 − γ ) = A2 πE (1 + γ ) i = : ne max X i0 ≥ [σ ] RL20i (1 − γ ) = A3 πE (1 + γ ) i = : ne max X i0 ≤ Num_Truss_Anls m 4πL20i (1 − γ ) X ≥ = A4 λmax (1 + γ ) i i = : ne max XX i = {max( A1, A2, A4, A5) + A3} X i0 = XX i X ≥ u 01 j X 01 j = A5 [u ] j = : 3nojnt j SÔ ĐỒ KHỐI XÁC ĐỊNH CÁC MA TRẬN CÁC HỆ SỐ TỪ CHUỖI CÁC CHƯƠNG TRÌNH TUYẾN TÍNH Start Nhập liệu toán thiết kế tối ưu Opt_Data.m Phân tích toán ứngvới biến thiết kế X = { X i }i =1:ne max Constrain_Func.m Xác định giá trị số {value _ cp j } Xác định hàm ràng buộc f j ({ X i }) j = : ne max + 6noj j = : ne max + 6nojnt Xác định ma trận hệ số {a}, {b}, {c} Từ chuỗi phương trình tuyến tính hóa quanh điểm { X i0 } a ji = f j' ( X i0 )  ne max b = − f ({ X }) + f j' ( X i0 ) X i0 + value _ cp j ∑ j j i  j =1  ' ci = Z ( X i )  i = : 8ne max + 6nojnt  j = : ne max + 6nojnt   SƠ ĐỒ CHƯƠNG TRÌNH CHÍNH TỰ ĐỘNG HOÁ THIẾ T KẾ TỐI ƯU KẾT CẤU THÉP DẠNG DÀN Start Nhập liệu toán thiết kế tối ưu Tính tóan thiết kế sơ chọn giá trị biến thiết kế ban đầu Opt_Data.m X0_Design.m XX = { XX i }i =1:ne max X = XX Xác định ma trận hệ số chuỗi chương trình tuyến tính {a ji },{b j },{ci } Công thức khai triển tuyến tính Taylor Opt_Frame3d_Anls m Giải tìm nghiệm tối ưu thuật tóan đơn hình OptimalSimplex.M XX Chenhlechmax No Yes Hội tụ Xử lý kết SƠ ĐỒ CHƯƠNG TRÌNH CHÍNH TỰ ĐỘNG HOÁ THIẾ T KẾ TỐI ƯU KẾT CẤU THÉP DẠNG DÀN (tiếp theo) Xử lý kết Dữ liệu tóan sau Phân tích kết cấu ứùng với giá trị biến thiết kế sau Opt_Frame3d_Anls m Xuất liệu dạng file.text Tên file nhập vào từ bàn phím Xuất liệu hình ảnh sơ đồ chuyển vị kết cấu Nội lực nút phần tử tọa độ địa phương Các thành phần chuyển vị nút tọa độ tổng thể Các thành phần phản lực gối tựa SƠ ĐỒ CHƯƠNG TRÌNH CHÍNH TỰ ĐỘNG HOÁ PHÂN TÍCH KẾT CẤU KHÔNG GIAN FRAME3D.M Start Nhập liệu toán Phân tích kết cấu Data.m Tính tóan ma trận độ cứng phần tử tọa độ tổng thể Beam3d.m Tính tóan ma trận chứa thứ tự bậc tự phần tử Fdmatrix.m Tính tóan ma trận độ cứng tổng thể kết cấu Assemk.m Tính tóan véc tơ tải nút tổng thể tác dụng tòan kết cấu Assemp.m Giải hệ, tìm ma trận chuyển vị nút tổng thể Tính nội lực, phản lực gối tựa,… Xuất liệu dạng file.text Tên file nhập vào từ bàn phím Xuất liệu hình ảnh sơ đồ chuyển vị kết cấu ... , bạn sức khỏe hạnh phúc ! TÓM TẮT LUẬN VĂN I TÊN ĐỀ TÀI : Tự động hóa thiết kế tối ưu kết cấu thép dạng dàn phương pháp phần tử hữu hạn có kể đến ràng buộc theo tiêu chuẩn Việt Nam II TÓM TẮT... việc thiết kế kết cấu theo tiêu chuẩn thiết kế họ Chính tác giả nhận thấy việc áp dụng phương pháp thiết kế tối ưu vào trình thiết kế kết cấu thép theo tiêu chuẩn Việt Nam cần thiết Hàm mục tiêu. .. hóa phân tích kết cấu PPPTHH; tự động hóa thiết lập ràng buộc, hàm mục tiêu; tự động hóa tuyến tính hóa hàm phi tuyến; tự động hóa tìm lời giải tối ưu thuất tóan đơn hình; tự động hóa xử lý kết