phân tích độ tin cậy của kết cấu tấm gia cường sử dụng phương pháp bề mặt đáp ứng cho hàm trạng thái giới hạn phi tuyến bậc cao

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỞ THÀNH PHỐ HỒ CHÍ MINH PHẠM NGUYỄN DUY HẢI PHÂN TÍCH ĐỘ TIN CẬY CỦA KẾT CẤU TẤM GIA CƯỜNG SỬ DỤNG PHƯƠNG PHÁP BỀ MẶT ĐÁP ỨNG CHO HÀM TRẠNG THÁI GIỚI HẠN PHI TUYẾN BẬC CAO Chuyên ngành : Xây dựng Công trình dân dụng Công nghiệp Mã số chuyên ngành : 60 58 02 08 LUẬN VĂN THẠC SỸ XÂY DỰNG Người hướng dẫn khoa học: TS NGUYỄN TRẦN CHÂN TP Hồ Chí Minh, Năm 2015 iii Tóm tắt Luận văn nhằm phân tích độ tin cậy kết cấu có gân gia cư ờng phương pháp bề mặt đáp ứng kết hợp với phương pháp mô Monte-Carlo Biến ngẫu nhiên chọn số mô-đun đàn hồi, tải trọng tác dụng bề dày Hàm trạng thái giới hạn chuyển vị giới hạn kết cấu Thuật toán phân tích độ tin cậy sử dụng luận văn gồm bước sau: Bước 1: Giải toán ứng xử tĩnh chuyển vị có gân gia cường, ứng xử phân tích ứng xử Mindlin, sử dụng phần tử tam giác CS-DSG3 Bước 2: Xấp xỉ hàm trạng thái giới hạn chuyển vị kết cấu phương pháp bề mặt đáp ứng Bước 3: Đánh giá độ tin cậy phương pháp Monte-Carlo với hàm trạng thái giới hạn chuyển vị giới hạn kết cấu Kết phân tích ứng xử tĩnh kết cấu có gân gia cường sử dụng phần tử tam giác CS-DSG3 lập trình ngôn ngữ Matlab so sánh kết tính toán ứng xử tĩnh phần mềm Ansys Kết phân tích độ tin cậy kết cấu có gân gia cường sử dụng phương pháp RSM-MCS so sánh với phương pháp FORM Từ khóa: Phương pháp phần tử hữu hạn (FEM), phương pháp đánh giá độ tin cậy MonteCarlo-Simulation (MCS), phương pháp bề mặt đáp ứng (RSM), phương pháp phần tử hữu hạn trơn (SFEM) iv Mục lục Lời cam đoan .i Lời cảm ơn .ii Tóm tắt iii Danh mục hình vẽ vii Danh mục bảng biểu ix Danh mục ký hiệu x Chữ viết tắt x Các hàm x Ma trận véc-tơ xi Các ký hiệu xiii Chương GIỚI THIỆU TỔNG QUAN 1.1 Giới thiệu 1.1.1 Tổng quan có gân gia cường 1.1.2 Tổng quan phương pháp phần tử hữu hạn (FEM) phần tử tam giác CS-DSG3 1.1.3 Tổng quan phương pháp bề mặt đáp ứng (RSM) 1.1.4 Tổng quan phương pháp đánh giá độ tin cậy 1.2 Tình hình nghiên cứu 1.2.1 Các công trình nghiên cứu nước 1.2.2 Các công trình nghiên cứu nước 1.3 Mục tiêu nghiên cứu 1.4 Phạm vi nghiên cứu 1.5 Cấu trúc luận văn v Chương CƠ SỞ LÝ THUYẾT 11 2.1 Lý thuyết dày Mindlin-Reissner 11 2.1.1 Khái niệm 11 2.1.2 Các giả thuyết 12 2.1.3 Qui ước dấu 12 2.1.4 Trường chuyển vị, ứng suất, biến dạng 15 2.1.5 Năng lượng biến dạng đàn hồi 18 2.2 Lý thuyết dầm Timoshenko 19 2.2.1 Trường chuyển vị, ứng suất, biến dạng 19 2.2.2 Năng lượng biến dạng đàn hồi dầm 23 2.3 Lý thuyết Reissner-Mindlin có gân gia cường 26 2.4 Phương pháp PTHH có gân gia cường 28 2.5 Phân tích ứng xử Mindlin sử dụng phần tử CS-DSG3 32 2.5.1 Tóm tắt phương pháp DSG3 32 2.5.2 Tóm tắt phương pháp CS-DSG3 35 2.5.3 Hệ phương trình có gân gia cường 36 2.6 Lý thuyết phương pháp bề mặt đáp ứng (Response Surface Method) 37 2.6.1 Phương pháp RSM cho xấp xỉ hàm trạng thái giới hạn 37 2.6.2 Phương pháp bề mặt đáp ứng (RSM) 39 2.7 Lý thuyết phân tích độ tin cậy 41 2.7.1 Phương pháp MCS (Monte Carlo Simulation) 41 2.7.2 Phương pháp RSM kết hợp với phương pháp Monte Carlo để đánh giá độ tin cậy 46 Chương CÁC KẾT QUẢ SỐ 47 3.1 Phân tích ứng xử tĩnh kết cấu có gân gia cường 47 3.1.1 Bài toán 47 3.1.1.1 Bài toán tĩnh học gia cường 48 vi 3.1.1.2 Bài toán tĩnh học có gia cường 50 3.1.2 Bài toán 53 3.1.2.1 Bài toán tĩnh học gia cường 54 3.1.2.2 Bài toán tĩnh học có gia cường 56 3.2 Phân tích độ tin cậy 59 3.2.1 Phân tích độ tin cậy hàm giải tích 59 3.2.1.1 Bài toán 59 3.2.1.2 Bài toán 61 3.2.2 So sánh kết tính toán FEM RSM-FEM 62 3.2.3 Phân tích độ tin cậy kết cấu có gân gia cường 64 3.2.3.1 Bài toán 64 3.2.3.2 Bài toán 67 Chương KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 70 4.1 Kết luận 70 4.2 Các vấn đề tồn hướng phát triển đề tài 71 Tài liệu tham khảo 72 Phụ lục 75 Phụ lục: Một số đoạn mã lập trình 75 vii Danh mục hình vẽ Hình 1.1: Tấm có gân gia cường dọc theo hướng chịu tải Hình 1.2: Tấm gia cường ứng dụng giao thông hay nhà dân dụng Hình 1.3: Tấm có gân gia cường ứng dụng kết cấu mái siêu thị, trạm xăng dầu Hình 1.4: Tấm có gân gia cường ứng dụng kết cấu bể chứa Hình 2.1: Tấm Mindlin-Reissner 11 Hình 2.2: Quy ước dấu 13 Hình 2.3: a) Quy ước dấu mặt phẳng Oxz; b) Quy ước dấu mặt phẳng Oyz 14 Hình 2.4: Biến dạng 15 Hình 2.5: Dầm hệ tọa độ địa phương Orsz 20 Hình 2.6: Quy ước dấu cho dầm 21 Hình 2.7: Biến dạng dầm 21 Hình 2.8: Đổi biến dầm 25 Hình 2.9: Tấm có gân gia cường theo hai phương x y 27 Hình 2.10: Tấm dầm rời rạc hóa tập hợp điểm nút 28 Hình 2.11: Phần tử tam giác điểm nút hệ tọa độ Descartes hệ tọa độ tự nhiên 32 Hình 2.12: Ba tam giác đ ược tạo từ tam giác 123 phương pháp CS-DSG3 35 Hình 2.13: a) bề mặt phản ứng f(x) độc lập b) đường đồng mức f(x) độc lập 38 Hình 2.14: a) bề mặt phản ứng f(x) có tương tác hai biến b) đường đồng mức f(x) có tương tác hai biến 38 viii Hình 2.15: a) bề mặt phản ứng f(x) bậc hai b) đường đồng mức f(x) bậc hai 39 Hình 2.16: Sơ đồ phương pháp Monte-Carlo 42 Hình 2.17: Sơ đồ biến ngẫu nhiên trục tọa độ số thực R 42 Hình 3.1: Tấm có gân gia cường 48 Hình 3.2: Chuyển vị gân gia cường cho toán (CS-DSG3) 49 Hình 3.3: Chuyển vị gân gia cường cho toán (Ansys) 49 Hình 3.4: Chuyển vị có gân gia cường cho toán (CS-DSG3) 50 Hình 3.5: Chuyển vị có gân gia cường cho toán (Ansys) 51 Hình 3.6: Quá trình hội tụ độ võng lớn cho toán 52 Hình 3.7: Tấm có hai gân gia cường 53 Hình 3.8: Chuyển vị gân gia cường cho toán (CS-DSG3) 55 Hình 3.9: Chuyển vị gân gia cường cho toán (Ansys) 55 Hình 3.10: Chuyển vị có gân gia cường cho toán (CS-DSG3) 57 Hình 3.11: Chuyển vị có gân gia cường cho toán (Ansys) 57 Hình 3.12: Quá trình hội tụ độ võng lớn cho toán 58 Hình 3.13: Hàm g=0 59 Hình 3.14: Hàm g hàm xấp xỉ g = 60 Hình 3.15: Biểu đồ so sánh xác suất phá hủy hàm chuyển vị cho toán 65 Hình 3.16: Xác suất phá hủy tăng dần tỉ lệ thay đổi biến P cho toán 66 Hình 3.17: Biểu đồ so sánh xác suất phá hủy hàm chuyển vị cho toán 68 Hình 3.18: Xác suất phá hủy tăng dần tỉ lệ thay đổi biến P cho toán 69 ix Danh mục bảng biểu Bảng 3.1: Dữ liệu toán 47 Bảng 3.2: Kết chuyển vị lớn gia cường cho toán 48 Bảng 3.3: Kết chuyển vị lớn có gia cường cho toán 50 Bảng 3.4: Dữ liệu toán 53 Bảng 3.5: Kết chuyển vị lớn gia cường cho toán 54 Bảng 3.6: Kết chuyển vị lớn có gia cường cho toán 56 Bảng 3.7: Kết xác suất phá hủy hàm g g cho toán 60 Bảng 3.8: Kết xác suất phá hủy hàm g g cho toán 61 Bảng 3.9: So sánh kết chuyển vị lớn FEM RSM-FEM có gân 62Error! Bookmark not defined Bảng 3.10: So sánh kết chuyển vị lớn FEM RSM-FEM có gân62 Bảng 3.11: Kết xác suất phá hủy không gia cường cho toán 64 Bảng 3.12: Kết xác suất phá hủy có gân gia cường cho toán 64 Bảng 3.13: Kết xác suất phá hủy VP thay đổi biến P cho toán 65 Bảng 3.14: Kết xác suất phá hủy không gia cường cho toán 67 Bảng 3.15: Kết xác suất phá hủy có gân gia cường cho toán 67 Bảng 3.16: Kết xác suất phá hủy VP thay đổi biến P cho toán 68 x Danh mục ký hiệu Chữ viết tắt FEM Finite Element Method - Phương pháp phần tử hữu hạn DSG3 Discrete Shear Gap- Phương pháp rời rạc độ lệch trượt áp dụng cho phần tử tam giác điểm nút CS-FEM Cell-based Smoothed Finite Element Method – Phương pháp phần tử hữu hạn trơn dựa phần tử CS-DSG3 Cell-based Smoothed Discrete Shear Gap – Phương pháp rời rạc độ lệch trượt trơn hóa dựa phần tử FSDT First- Order Shear Deformable Theory – Lý thuyết biến dạng cắt bậc FORM First Order Reliability Method– Phương pháp đánh giá độ tin cậy bậc SORM Second Order Reliability Method–Phương pháp đánh giá độ tin cậy bậc hai PTHH Phần tử hữu hạn RSM Response Surface Method–Phương pháp bề mặt đáp ứng MCS Monte-Carlo Simulation–Thuật toán Monte-Carlo ASF Assumed Strain Field– Phương pháp giả định trường biến dạng TBT Timoshenko Beam Theory–Lý thuyết dầm Timoshenko CDF Cumulative Distribution Function–Hàm phân bố xác suất Các hàm Φ (.) Hàm tích lũy Gaussian tiêu chuẩn g ( x) Hàm trạng thái giới hạn g (x) Hàm đa thức FX ( x) Hàm xác suất dùng cho hai biến ngẫu nhiên liên tục rời rạc Φ −1 Nghịch đảo hàm phân phối tích lũy xác suất dạng chuẩn xi TM Hàm đa thức Chebyshev loại Ma trận véc-tơ x Véc-tơ biến thiết kế y Véc-tơ biến ngẫu nhiên u Véc-tơ biến chuẩn hóa u (pe ) Véc-tơ hàm chuyển vị phần tử u b( e ) Véc-tơ hàm chuyển vị phần tử dầm εb Véc-tơ biến dạng uốn εm Véc-tơ biến dạng màng γ Véc-tơ biến dạng cắt σb Véc-tơ ứng suất uốn σm Véc-tơ ứng suất màng σs Véc-tơ ứng suất cắt εb Véc-tơ biến dạng dầm σb Véc-tơ ứng suất dầm d (pe ) Véc-tơ chuyển vị nút phần tử tam giác điểm nút d (pie ) Véc-tơ chuyển vị nút phần tử tam giác điểm nút nút thứ i d b( e ) Véc-tơ chuyển vị nút phần tử điểm nút d bi( e ) Véc-tơ chuyển vị nút phần tử điểm nút nút thứ i dp Véc-tơ chuyển vị nút tổng thể db Véc-tơ chuyển vị nút tổng thể dầm d Véc-tơ chuyển vị nút tổng thể có gân gia cường f Véc-tơ tải trọng tác dụng 60 Hình 3.14: Hàm g hàm xấp xỉ g = Nhận xét: ta nhận thấy đồ thị hàm trùng với nhau, nói xấp xỉ hàm số g phương pháp RSM cho ta kết tương đối xác Tiếp theo, tác giả tiến hành phân tích độ tin cậy hàm số g phương pháp Monte-Carlo FORM hàm số g phương pháp Monte-Carlo, với biến ngẫu nhiên trình bày Bảng 3.7: Kết xác suất phá hủy hàm g g cho toán Phương pháp MCS (N=10 ) RSM-MCS FORM Xác suất phá hủy Pf 2.46x10-3 2.39x10-3 3.90 x10-2 Thời gian tính toán (s) 35s 21s 1s 61 Bảng 3.7 cho thấy thời gian tính toán xác suất phá hủy hàm trạng thái giới hạn g dùng phương pháp FORM giải nhanh kết có sai số định Kết xác suất phá hủy dùng phương pháp RSM-MCS cho kết gần giống với phương pháp MCS 3.2.1.2 Bài toán Xét hàm số g= 0.5(x-1) 2y- y2+z3-0.2sin(xy)+1.149 Phân tích độ tin cậy hàm số này, với x, y z biến ngẫu nhiên phân bố khoảng [0 1] Biết phá hủy hàm g hàm g [...]... dụng phương pháp bề mặt đáp ứng (RSM) kết hợp với phương pháp Monte-Carlo, trong đó phương pháp RSM được sử dụng để xấp xỉ hàm trạng thái giới hạn và phương pháp Monte-Carlo được sử dụng để phân tích độ tin cậy cho kết cấu tấm gia cường theo các biến ngẫu nhiên Phương pháp kết hợp này sẽ khắc phục được những nhược điểm của các phương pháp trên 1.1.1 Tổng quan về tấm có gân gia cường So với kết cấu tấm. .. luận văn này, tác giả sẽ sử dụng phương pháp bề mặt đáp ứng (RSM) nhằm xấp xỉ hàm trạng thái giới hạn, đồng thời kết hợp với phương pháp Monte-Carlo để đánh giá độ tin cậy của kết cấu tấm có gân gia cường 3 Bài toán phân tích ứng xử kết cấu tấm có gân gia cường được thực hiện bằng nhiều phương pháp khác nhau như phương pháp giải tích, phương pháp số Tuy nhiên phương pháp giải tích chỉ giải quyết những... đồng thời cũng phân tích độ tin cậy của loại kết cấu này khi sử dụng cho sân bay nổi trên biển 8 + Lr.wally, Orisamolu, K-Tma [4] đã sử dụng phương pháp phần tử hữu hạn ngẫu nhiên để phân tích độ tin cậy của kết cấu tấm có gân gia cường + S Peng-Cheng, H Dade và W Zongmu [5] đã sử dụng hàm B-spline như hàm tọa độ để phân tích ứng xử tĩnh học, động học và ổn định của tấm thép có gân gia cường + T P Holopainen... đã sử dụng phương pháp PTHH để phân tích dao động tự do của tấm có gân gia cường lệch tâm, trong đó tác giả đã xem xét ứng xử của phần tử tấm theo lý thuyết tấm Reissener-Mindlin + M Mukhopadhyay [7] đã sử dụng phương pháp sai phân hữu hạn để phân tích dao động và ổn định tấm có gân gia cường Từ kết quả tham khảo của các bài báo trên ta thấy việc phân tích độ tin cậy của kết cấu tấm có gân gia cường. .. thế giới, bài toán phân tích độ tin cậy của kết cấu tấm có gân gia cường được thực hiện bằng nhiều phương pháp khác nhau Một số công trình nghiên cứu tiêu biểu được liệt kê như sau: + He Jian, Chen Xiaoyan [2] đã phân tích độ tin cậy của kết cấu tấm có gân gia cường dựa trên phương pháp entropy cực đại + Okada, và các cộng sự [3] đã sử dụng phương pháp phần tử hữu hạn phi tuyến để tính toán tấm gia cường, ... như xây dựng dân dụng, tàu biển, hàng không và kỹ thuật cơ khí Cho đến nay, nhiều mô hình kết cấu tấm khác nhau đã ra đời như kết cấu tấm gấp, tấm lượn sóng, tấm có gân gia cường, trong đó tấm có gân gia cường là kết cấu tấm được gia cường bởi các gân và được sử dụng khá rộng rãi trong lĩnh vực xây dựng Mặc dù tấm gia cường có tính ứng dụng cao nhưng việc phân tích ứng xử của kết cấu này lại tương... FORM do hàm trạng thái giới hạn được xấp xỉ thành đa thức bậc hai nhưng việc tính toán lại tương đối phức tạp Để khắc phục những nhược điểm của các phương pháp trên, cần có những phương pháp xấp xỉ hàm trạng thái giới hạn thông minh hơn như phương pháp bề mặt đáp ứng (RSM) hay phương pháp mạng thần kinh nhân tạo (ANN) và sau đó kết hợp với phương pháp Monte-Carlo để đánh giá độ tin cậy của kết cấu nhằm... cậy tấm có dầm gia cường Hàm trạng thái giới hạn là giới hạn chuyển vị, còn các biến ngẫu nhiên được chọn là hằng số mô-đun đàn hồi, tải trọng tác dụng và chiều dày tấm Nghiên cứu được thực hiện bao gồm các công việc sau: • Nắm vững lý thuyết và lập trình code Matlab bằng phần tử CS-DSG3 để phân tích ứng xử kết cấu tấm gia cường; • Xấp xỉ hàm trạng thái giới hạn bằng phương pháp bề mặt đáp ứng; • Phân. .. Phân tích độ tin cậy của tấm gia cường bằng phương pháp Monte-Carlo với hàm trạng thái giới hạn bị ràng buộc về chuyển vị 1.4 Phạm vi nghiên cứu Nội dung nghiên cứu của luận văn được thực hiện trong phạm vi sau: • Phân tích ứng xử tĩnh kết cấu tấm có gân gia cường • Ứng xử của vật liệu là đàn hồi tuyến tính • Biến dạng của kết cấu được coi là biến dạng nhỏ • Biến ngẫu nhiên là mô-đun đàn hồi, tải tác dụng. .. không an toàn Do đó, các phương pháp phân tích độ tin cậy được đề xuất nhằm tính toán xác suất phá hủy hay độ an toàn của kết cấu Có một số phương pháp đánh giá độ tin cậy phổ biến hiện nay gồm: • Phương pháp mô phỏng Monte-Carlo (Monte Carlo Simulation – MCS) • Phương pháp đánh giá độ tin cậy bậc nhất (First Order Reliability Method– FORM) • Phương pháp đánh giá độ tin cậy bậc hai (Second Order Reliability ... CS-DSG3 để phân tích ứng xử kết cấu gia cường; • Xấp xỉ hàm trạng thái giới hạn phương pháp bề mặt đáp ứng; • Phân tích độ tin cậy gia cường phương pháp Monte-Carlo với hàm trạng thái giới hạn bị... nhằm phân tích độ tin cậy kết cấu gia cường sử dụng phương pháp bề mặt đáp ứng (RSM) kết hợp với phương pháp Monte-Carlo, phương pháp RSM sử dụng để xấp xỉ hàm trạng thái giới hạn phương pháp. .. dụng phần tử tam giác CS-DSG3 Bước 2: Xấp xỉ hàm trạng thái giới hạn chuyển vị kết cấu phương pháp bề mặt đáp ứng Bước 3: Đánh giá độ tin cậy phương pháp Monte-Carlo với hàm trạng thái giới hạn