Đánh giá hư hại của khung bê tông cốt thép chịu động đất có xét đến tương tác giữa đất nền và kết cấu

ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA LÊ VĂN HÙNG ĐÁNH GIÁ HƯ HẠI CỦA KHUNG BÊ TƠNG CỐT THÉP CHỊU ĐỘNG ĐẤT CĨ XÉT ĐẾN TƯƠNG TÁC GIỮA ĐẤT NỀN VÀ KẾT CẤU Chuyên ngành: Kỹ thuật xây dựng cơng trình dân dụng cơng nghiệp Mã số: 60 58 02 08 LUẬN VĂN THẠC SĨ TP HỒ CHÍ MINH, tháng 06 năm 2018 Cơng trình hoàn thành tại: Trường Đại Học Bách Khoa – ĐHQG - HCM Cán hướng dẫn khoa học: TS Cao Văn Vui…………………………………… Cán chấm nhận xét 1: ………………………………………………………… Cán chấm nhận xét 2: ………………………………………………………… Luận văn thạc sĩ bảo vệ Trường Đại học Bách Khoa, ĐHQG Tp.HCM, ngày….tháng… năm 2018 Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm: CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG ĐẠI HỌC QUỐC GIA TP.HCM CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM TRƯỜNG ĐẠI HỌC BÁCH KHOA Độc Lập - Tự Do - Hạnh Phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: LÊ VĂN HÙNG MSHV: 1570159 Ngày, tháng, năm sinh: 28/08/1987 Nơi sinh: Hà Tĩnh Chuyên ngành: Kỹ thuật xây dựng dân dụng công nghiệp Mã số: 60 58 02 08 I TÊN ĐỀ TÀI: Đánh giá hư hại khung bê tơng cốt thép chịu động đất có xét đến tương tác đất kết cấu NHIỆM VỤ VÀ NỘI DUNG: Nghiên cứu mơ hình tương tác kết cấu đất - Soil Structure Interaction models - SSI) Từ đó, lựa chọn mơ hình để sử dụng luận văn Nghiên cứu mơ hình đánh giá mức độ hư hại khung bê tơng cốt thép (BTCT) Mơ hình khung BTCT có sử dụng phần tử phi tuyến ứng xử trễ Phân tích phi tuyến theo lịch sử thời gian khung BTCT chịu động đất, có khơng có xét SSI Từ kết phân tích phi tuyến theo lịch sử thời gian, số hư hại (Damage Index –DI) tính tốn cho khung BTCT có khơng có SSI chịu trận động đất Ảnh hưởng SSI đánh giá cách so sánh số DI khung BTCT có xét SSI số DI khung BTCT khơng có SSI V NGÀY GIAO NHIỆM VỤ: 04 - 09- 2017 VI NGÀY HOÀN THÀNH NHIỆM VỤ: 18 – 06 - 2018 VII CÁN BỘ HƯỚNG DẪN: TS Cao Văn Vui Tp HCM, ngày tháng năm 2018 CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM BỘ MÔN ĐÀO TẠO (Họ tên chữ ký) (Họ tên chữ ký) TS Cao Văn Vui TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG (Họ tên chữ ký) LỜI CẢM ƠN Luận văn thạc sĩ Xây dựng cơng trình dân dụng cơng nghiệp nằm hệ thống luận cuối khóa nhằm trang bị cho Học viên cao học khả tự nghiên cứu, biết cách giải vấn đề cụ thể đặt thực tế xây dựng… Đó trách nhiệm niềm tự hào học viên cao học Để hoàn thành luận văn này, cố gắng nỗ lực thân, nhận giúp đỡ nhiều từ tập thể cá nhân Tơi xin ghi nhận tỏ lịng biết ơn đến tập thể cá nhân dành cho tơi giúp đỡ q báu Đầu tiên tơi xin bày tỏ lịng biết ơn sâu sắc đến thầy TS Cao Văn Vui Thầy đưa gợi ý để hình thành nên ý tưởng đề tài Thầy góp ý cho tơi nhiều cách nhận định đắn vấn đề nghiên cứu, cách tiếp cận nghiên cứu hiệu Tôi xin chân thành cảm ơn quý Thầy Cô Khoa Kỹ Thuật Xây Dựng, trường Đại học Bách Khoa Tp.HCM truyền dạy kiến thức quý giá cho tơi, kiến thức khơng thể thiếu đường nghiên cứu khoa học nghiệp sau Luận văn thạc sĩ hoàn thành thời gian quy định với nỗ lực thân, nhiên khơng thể khơng có thiếu sót Kính mong q Thầy Cơ dẫn thêm để tơi bổ sung kiến thức hồn thiện thân Xin trân trọng cảm ơn Tp HCM, ngày 18 tháng 06 năm 2018 Lê Văn Hùng TÓM TẮT LUẬN VĂN THẠC SĨ Trong tiêu chuẩn thiết kế động đất nay, tương tác đất kết cấu (Soil Structure Interaction – SSI) thường bỏ qua SSI xem có lợi Tuy nhiên, bỏ qua ảnh hưởng SSI không phản ánh sát trình làm việc thực tế kết cấu chịu động đất Mục đích nghiên cứu nhằm đánh giá ảnh hưởng SSI đến mức độ hư hại khung bê tông cốt thép (BTCT) chịu động đất với cường độ khác Để đạt mục đích này, khung BTCT tầng tầng có khơng có SSI mơ hình hóa sử dụng phần tử phi tuyến ứng xữ trễ Mơ hình khung kiểm chứng cách so sánh với kết phân tích tác giả khác Sau kiểm chứng phân tích lịch sử thời gian khung chịu động đất với cường độ động đất khác Kết phân tích lịch sử thời gian sử dụng để tính số hư hại (Damage Index - DI) Kết phân tích hư hại khung có SSI so sánh với khung khơng có SSI Kết so sánh cho thấy xét đến SSI mức độ hư hại khung tăng lên so với khung không xét đến SSI, SSI có hiệu ứng bất lợi cho kết cấu Do vậy, đánh giá cơng trình cũ hay thiết kế cơng trình nên cân nhắc đưa SSI vào phân tích kết cấu Kết nghiên cứu cho thấy ảnh hưởng SSI giảm vận tốc sóng cắt (Vs ) tăng ABSTRACT In current seismic design codes, the effect of soil – structure interaction (SSI) is neglected because the effect of SSI is conceptually beneficial However, neglecting the effect of SSI may lead to inaccurate response of structures under seismic excitations This study aims at assessing the effect of SSI on the damage of reinforced concrete structures subjected to different seismic intensities To achieve this goal, 4-storey and 8-storey reinforced concrete frames with and without SSI are modeled using hysteretic nonlinear elements Models of these frames are compared with the results conducted by other researchers, showing good approximations After verifying, time history analyses of these frames subjected to different earthquake intensities are conducted The data of time history analyses are then used for damage analyses, using a damage model The damage states of these frames with and without SSI are compared The results of comparisons show that the effect of SSI increases the damage of these frames, showing the detrimental effect of SSI Thus, the effect of SSI should be considered in evaluating existing structures or designing new structures The results also show that the effect of SSI decreases with increase of shear velocity Vs LỜI CAM ĐOAN Tôi xin cam đoan cơng việc tơi thực hướng dẫn thầy TS Cao Văn Vui Các kết luận văn thật chưa công bố nghiên cứu khác Tôi xin chịu trách nhiệm cơng việc thực Tp HCM, ngày 18 tháng 06 năm 2018 Lê Văn Hùng Mục lục CHƯƠNG 1: MỞ ĐẦU 1.1 Lý chọn đề tài 1.2 Mục đích nghiên cứu 1.3 Phạm vi nghiên cứu 1.4 Phương pháp nghiên cứu 1.5 Ý nghĩa nghiên cứu 1.6 Tính cấp thiết đề tài 1.7 Cấu trúc luận văn CHƯƠNG 2: TỔNG QUAN 2.1 Sơ lược động đất 2.2 Tương tác đất kết cấu (SSI) 2.3 Các cơng trình nghiên cứu giới 2.4 Các cơng trình nghiên cứu nước 17 2.5 Kết luận chương 19 CHƯƠNG 3: CƠ SỞ LÝ THUYẾT 20 3.1 Mơ hình SSI 20 3.2 Quan hệ ứng suất - biến dạng vật liệu 25 3.2.1 Quan hệ ứng suất biến dạng bê tông 25 3.2.2 Quan hệ ứng suất biến dạng thép 28 3.3 Mối quan hệ mơ men - góc xoay, ứng xử trễ phân tích phi tuyến 29 3.3.1 Quan hệ mơ men - góc xoay 29 3.3.2 Ứng xử trễ cấu kiện BTCT 30 3.4 Mơ hình phân tích phi tuyến 31 3.5 Mơ hình hư hại 32 CHƯƠNG 4: PHÂN TÍCH HƯ HẠI CỦA KHUNG BTCT 38 4.1 Khung tầng 38 4.1.1 Mô tả khung tầng 38 4.1.2 Mơ hình kết cấu 39 4.1.3 Mơ hình SSI 41 4.1.4 Kiểm chứng mơ hình 42 4.1.5 Phân tích phi tuyến theo thời gian 45 4.1.6 Phân tích hư hại 47 4.1.7 Kết thảo luận 47 4.2 Khung tầng 64 4.2.1 Mô tả khung tầng 64 4.2.2 Mơ hình kết cấu 65 4.2.3 Mơ hình SSI 68 4.2.4 Kiểm chứng mơ hình sử dụng phân tích đẩy dần (pushover) 68 4.2.5 Phân tích phi tuyến theo thời gian 69 4.2.6 Phân tích hư hại 71 4.2.6 Kết thảo luận 86 4.2.7 Phân tích hư hại khung với vận tốc sóng cắt khác 89 CHƯƠNG 5: KẾT LUẬN VÀ KIẾN NGHỊ 100 5.1 Kết luận 100 5.2 Kiến nghị 100 CƠNG TRÌNH CƠNG BỐ 101 TÀI LIỆU THAM KHẢO 108 PHỤ LỤC…………………………………………………………………… 116 [MRLink104]=[MRLinkType19]; [MRLink105]=[MRLinkType19]; [MRLink106]=[MRLinkType19]; [MRLink107]=[MRLinkType19]; [MRLink108]=[MRLinkType19]; [MRLink109]=[MRLinkType19]; [MRLink110]=[MRLinkType19]; [MRLink111]=[MRLinkType19]; [MRLink112]=[MRLinkType19]; [MRAllLink]= [ [MRLink1] [MRLink2] [MRLink3] [MRLink4] [MRLink5] [MRLink6] [MRLink7] [MRLink8] [MRLink9] [MRLink10] [MRLink11] [MRLink12] [MRLink13] [MRLink14] [MRLink15] [MRLink16] [MRLink17] [MRLink18] [MRLink19] [MRLink20] [MRLink21] [MRLink22] [MRLink23] [MRLink24] [MRLink25] [MRLink26] [MRLink27] [MRLink28] [MRLink29] [MRLink30] [MRLink31] [MRLink32] [MRLink33] [MRLink34] [MRLink35] [MRLink36] [MRLink37] [MRLink38] [MRLink39] [MRLink40] [MRLink41] [MRLink42] [MRLink43] [MRLink44] [MRLink45] [MRLink46] [MRLink47] [MRLink48] [MRLink49] [MRLink50] [MRLink51] [MRLink52] [MRLink53] [MRLink54] [MRLink55] [MRLink56] [MRLink57] [MRLink58] [MRLink59] [MRLink60] [MRLink61] [MRLink62] [MRLink63] [MRLink64] [MRLink65] [MRLink66] [MRLink67] [MRLink68] [MRLink69] [MRLink70] [MRLink71] [MRLink72] [MRLink73] [MRLink74] [MRLink75] [MRLink76] [MRLink77] [MRLink78] [MRLink79] [MRLink80] [MRLink81] [MRLink82] [MRLink83] [MRLink84] [MRLink85] [MRLink86] [MRLink87] [MRLink88] [MRLink89] [MRLink90] [MRLink91] [MRLink92] [MRLink93] [MRLink94] [MRLink95] [MRLink96] [MRLink97] [MRLink98] [MRLink99] [MRLink100] [MRLink101] [MRLink102] [MRLink103] [MRLink104] [MRLink105] [MRLink106] [MRLink107] [MRLink108] [MRLink109] [MRLink110] [MRLink111] [MRLink112] ]; %% Xu li ket qua Moment-Rotation tu SAP2000 analysis % Time and Rotation NoOTS=501; %Number of Output Time Steps (taken from SAP) + 1; NOF1=1; %Can duoi Link duoc ve hinh NOF2=4; %Can tren Link dc ve hinh [tRAll]=xlsread('DeformationsNgam.xlsx','Element Deformations - Links'); [tRLink1]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink2]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink3]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink4]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink5]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink6]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink7]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink8]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink9]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink10]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink11]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink12]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink13]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink14]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink15]=tRAll(1:NoOTS,[1 [tRLink16]=tRAll(1:NoOTS,[1 [tRLink17]=tRAll(1:NoOTS,[1 [tRLink18]=tRAll(1:NoOTS,[1 [tRLink19]=tRAll(1:NoOTS,[1 [tRLink20]=tRAll(1:NoOTS,[1 [tRLink21]=tRAll(1:NoOTS,[1 [tRLink22]=tRAll(1:NoOTS,[1 [tRLink23]=tRAll(1:NoOTS,[1 [tRLink24]=tRAll(1:NoOTS,[1 [tRLink25]=tRAll(1:NoOTS,[1 [tRLink26]=tRAll(1:NoOTS,[1 [tRLink27]=tRAll(1:NoOTS,[1 [tRLink28]=tRAll(1:NoOTS,[1 [tRLink29]=tRAll(1:NoOTS,[1 [tRLink30]=tRAll(1:NoOTS,[1 [tRLink31]=tRAll(1:NoOTS,[1 [tRLink32]=tRAll(1:NoOTS,[1 [tRLink33]=tRAll(1:NoOTS,[1 [tRLink34]=tRAll(1:NoOTS,[1 [tRLink35]=tRAll(1:NoOTS,[1 [tRLink36]=tRAll(1:NoOTS,[1 [tRLink37]=tRAll(1:NoOTS,[1 [tRLink38]=tRAll(1:NoOTS,[1 [tRLink39]=tRAll(1:NoOTS,[1 [tRLink40]=tRAll(1:NoOTS,[1 [tRLink41]=tRAll(1:NoOTS,[1 [tRLink42]=tRAll(1:NoOTS,[1 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; [tRLink43]=tRAll(1:NoOTS,[1 [tRLink44]=tRAll(1:NoOTS,[1 [tRLink45]=tRAll(1:NoOTS,[1 [tRLink46]=tRAll(1:NoOTS,[1 [tRLink47]=tRAll(1:NoOTS,[1 [tRLink48]=tRAll(1:NoOTS,[1 [tRLink49]=tRAll(1:NoOTS,[1 [tRLink50]=tRAll(1:NoOTS,[1 [tRLink51]=tRAll(1:NoOTS,[1 [tRLink52]=tRAll(1:NoOTS,[1 [tRLink53]=tRAll(1:NoOTS,[1 [tRLink54]=tRAll(1:NoOTS,[1 [tRLink55]=tRAll(1:NoOTS,[1 [tRLink56]=tRAll(1:NoOTS,[1 [tRLink57]=tRAll(1:NoOTS,[1 [tRLink58]=tRAll(1:NoOTS,[1 [tRLink59]=tRAll(1:NoOTS,[1 [tRLink60]=tRAll(1:NoOTS,[1 [tRLink61]=tRAll(1:NoOTS,[1 [tRLink62]=tRAll(1:NoOTS,[1 [tRLink63]=tRAll(1:NoOTS,[1 [tRLink64]=tRAll(1:NoOTS,[1 [tRLink65]=tRAll(1:NoOTS,[1 [tRLink66]=tRAll(1:NoOTS,[1 [tRLink67]=tRAll(1:NoOTS,[1 [tRLink68]=tRAll(1:NoOTS,[1 [tRLink69]=tRAll(1:NoOTS,[1 [tRLink70]=tRAll(1:NoOTS,[1 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); 7]); tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; tRAll(1:NoOTS,:)=[]; [tRLink71]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink72]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink73]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink74]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink75]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink76]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink77]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink78]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink79]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink80]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink81]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink82]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink83]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink84]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink85]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink86]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink87]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink88]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink89]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink90]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink91]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink92]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink93]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink94]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink95]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink96]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink97]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink98]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink99]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink100]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink101]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink102]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink103]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink104]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink105]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink106]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink107]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink108]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink109]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink110]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink111]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; [tRLink112]=tRAll(1:NoOTS,[1 7]); tRAll(1:NoOTS,:)=[]; %Moment [MAllLink0]=xlsread('ForcesNgam.xlsx','Element Forces - Links'); [m,n]=size(MAllLink0); MAllLink=MAllLink0; for i=1:(m-1) if MAllLink0(m-i,1)== MAllLink0(m-i+1,1) MAllLink(m-i+1,:)=[]; end end clear MAllLink0 [MLink1]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink2]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink3]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink4]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink5]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink6]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink7]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink8]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink9]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink10]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink11]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink12]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink13]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink14]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink15]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink16]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink17]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink18]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink19]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink20]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink21]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink22]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink23]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink24]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink25]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink26]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink27]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink28]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink29]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink30]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink31]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink32]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink33]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink34]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink35]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink36]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink37]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink38]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink39]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink40]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink41]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink42]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink43]=MAllLink(1:NoOTS,7); [MLink44]=MAllLink(1:NoOTS,7); [MLink45]=MAllLink(1:NoOTS,7); [MLink46]=MAllLink(1:NoOTS,7); [MLink47]=MAllLink(1:NoOTS,7); [MLink48]=MAllLink(1:NoOTS,7); [MLink49]=MAllLink(1:NoOTS,7); [MLink50]=MAllLink(1:NoOTS,7); [MLink51]=MAllLink(1:NoOTS,7); [MLink52]=MAllLink(1:NoOTS,7); [MLink53]=MAllLink(1:NoOTS,7); [MLink54]=MAllLink(1:NoOTS,7); [MLink55]=MAllLink(1:NoOTS,7); [MLink56]=MAllLink(1:NoOTS,7); [MLink57]=MAllLink(1:NoOTS,7); [MLink58]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; MAllLink(1:NoOTS,:)=[]; [MLink59]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink60]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink61]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink62]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink63]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink64]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink65]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink66]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink67]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink68]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink69]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink70]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink71]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink72]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink73]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink74]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink75]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink76]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink77]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink78]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink79]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink80]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink81]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink82]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink83]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink84]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink85]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink86]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink87]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink88]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink89]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink90]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink91]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink92]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink93]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink94]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink95]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink96]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink97]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink98]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink99]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink100]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink101]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink102]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink103]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink104]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink105]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink106]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink107]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink108]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink109]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink110]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink111]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [MLink112]=MAllLink(1:NoOTS,7); MAllLink(1:NoOTS,:)=[]; [Link1]=[tRLink1 MLink1]; [Link2]=[tRLink2 MLink2]; [Link3]=[tRLink3 MLink3]; [Link4]=[tRLink4 MLink4]; [Link5]=[tRLink5 MLink5]; [Link6]=[tRLink6 MLink6]; [Link7]=[tRLink7 MLink7]; [Link8]=[tRLink8 MLink8]; [Link9]=[tRLink9 MLink9]; [Link10]=[tRLink10 MLink10]; [Link11]=[tRLink11 MLink11]; [Link12]=[tRLink12 MLink12]; [Link13]=[tRLink13 MLink13]; [Link14]=[tRLink14 MLink14]; [Link15]=[tRLink15 MLink15]; [Link16]=[tRLink16 MLink16]; [Link17]=[tRLink17 MLink17]; [Link18]=[tRLink18 MLink18]; [Link19]=[tRLink19 MLink19]; [Link20]=[tRLink20 MLink20]; [Link21]=[tRLink21 MLink21]; [Link22]=[tRLink22 MLink22]; [Link23]=[tRLink23 MLink23]; [Link24]=[tRLink24 MLink24]; [Link25]=[tRLink25 MLink25]; [Link26]=[tRLink26 MLink26]; [Link27]=[tRLink27 MLink27]; [Link28]=[tRLink28 MLink28]; [Link29]=[tRLink29 MLink29]; [Link30]=[tRLink30 MLink30]; [Link31]=[tRLink31 MLink31]; [Link32]=[tRLink32 MLink32]; [Link33]=[tRLink33 MLink33]; [Link34]=[tRLink34 MLink34]; [Link35]=[tRLink35 MLink35]; [Link36]=[tRLink36 MLink36]; [Link37]=[tRLink37 MLink37]; [Link38]=[tRLink38 MLink38]; [Link39]=[tRLink39 MLink39]; [Link40]=[tRLink40 MLink40]; [Link41]=[tRLink41 MLink41]; [Link42]=[tRLink42 MLink42]; [Link43]=[tRLink43 [Link44]=[tRLink44 [Link45]=[tRLink45 [Link46]=[tRLink46 [Link47]=[tRLink47 [Link48]=[tRLink48 [Link49]=[tRLink49 [Link50]=[tRLink50 [Link51]=[tRLink51 [Link52]=[tRLink52 [Link53]=[tRLink53 [Link54]=[tRLink54 [Link55]=[tRLink55 [Link56]=[tRLink56 [Link57]=[tRLink57 [Link58]=[tRLink58 MLink43]; MLink44]; MLink45]; MLink46]; MLink47]; MLink48]; MLink49]; MLink50]; MLink51]; MLink52]; MLink53]; MLink54]; MLink55]; MLink56]; MLink57]; MLink58]; [Link59]=[tRLink59 MLink59]; [Link60]=[tRLink60 MLink60]; [Link61]=[tRLink61 MLink61]; [Link62]=[tRLink62 MLink62]; [Link63]=[tRLink63 MLink63]; [Link64]=[tRLink64 MLink64]; [Link65]=[tRLink65 MLink65]; [Link66]=[tRLink66 MLink66]; [Link67]=[tRLink67 MLink67]; [Link68]=[tRLink68 MLink68]; [Link69]=[tRLink69 MLink69]; [Link70]=[tRLink70 MLink70]; [Link71]=[tRLink71 MLink71]; [Link72]=[tRLink72 MLink72]; [Link73]=[tRLink73 MLink73]; [Link74]=[tRLink74 MLink74]; [Link75]=[tRLink75 MLink75]; [Link76]=[tRLink76 MLink76]; [Link77]=[tRLink77 MLink77]; [Link78]=[tRLink78 MLink78]; [Link79]=[tRLink79 MLink79]; [Link80]=[tRLink80 MLink80]; [Link81]=[tRLink81 MLink81]; [Link82]=[tRLink82 MLink82]; [Link83]=[tRLink83 MLink83]; [Link84]=[tRLink84 MLink84]; [Link85]=[tRLink85 MLink85]; [Link86]=[tRLink86 MLink86]; [Link87]=[tRLink87 MLink87]; [Link88]=[tRLink88 MLink88]; [Link89]=[tRLink89 MLink89]; [Link90]=[tRLink90 MLink90]; [Link91]=[tRLink91 MLink91]; [Link92]=[tRLink92 MLink92]; [Link93]=[tRLink93 MLink93]; [Link94]=[tRLink94 MLink94]; [Link95]=[tRLink95 MLink95]; [Link96]=[tRLink96 MLink96]; [Link97]=[tRLink97 MLink97]; [Link98]=[tRLink98 MLink98]; [Link99]=[tRLink99 MLink99]; [Link100]=[tRLink100 MLink100]; [Link101]=[tRLink101 MLink101]; [Link102]=[tRLink102 MLink102]; [Link103]=[tRLink103 MLink103]; [Link104]=[tRLink104 MLink104]; [Link105]=[tRLink105 MLink105]; [Link106]=[tRLink106 MLink106]; [Link107]=[tRLink107 MLink107]; [Link108]=[tRLink108 MLink108]; [Link109]=[tRLink109 MLink109]; [Link110]=[tRLink110 MLink100]; [Link111]=[tRLink111 MLink111]; [Link112]=[tRLink112 MLink112]; clear MLink1 MLink2 MLink3 MLink4 MLink5 MLink6 MLink7 MLink8 MLink10; clear MLink11 MLink12 MLink13 MLink14 MLink15 MLink16 MLink17 MLink19 MLink20; clear MLink21 MLink22 MLink23 MLink24 MLink25 MLink26 MLink27 MLink29 MLink30; clear MLink31 MLink32 MLink33 MLink34 MLink35 MLink36 MLink37 MLink39 MLink40; clear MLink41 MLink42 MLink43 MLink44 MLink45 MLink46 MLink47 MLink49 MLink50; clear MLink51 MLink52 MLink53 MLink54 MLink55 MLink56 MLink57 MLink59 MLink60; clear MLink61 MLink62 MLink63 MLink64 MLink65 MLink66 MLink67 MLink69 MLink70; clear MLink71 MLink72 MLink73 MLink74 MLink75 MLink76 MLink77 MLink79 MLink80; clear MLink81 MLink82 MLink83 MLink84; MLink9 MLink18 MLink28 MLink38 MLink48 MLink58 MLink68 MLink78 clear tRLink1 tRLink2 tRLink3 tRLink4 tRLink5 tRLink6 tRLink7 tRLink8 tRLink9 tRLink10; clear tRLink11 tRLink12 tRLink13 tRLink14 tRLink15 tRLink16 tRLink17 tRLink18 tRLink19 tRLink20; clear tRLink21 tRLink22 tRLink23 tRLink24 tRLink25 tRLink26 tRLink27 tRLink28 tRLink29 tRLink30; clear tRLink31 tRLink32 tRLink33 tRLink34 tRLink35 tRLink36 tRLink37 tRLink38 tRLink39 tRLink40; clear tRLink41 tRLink42 tRLink43 tRLink44 tRLink45 tRLink46 tRLink47 tRLink48 tRLink49 tRLink50; clear tRLink51 tRLink52 tRLink53 tRLink54 tRLink55 tRLink56 tRLink57 tRLink58 tRLink59 tRLink60; clear tRLink61 tRLink62 tRLink63 tRLink64 tRLink65 tRLink66 tRLink67 tRLink68 tRLink69 tRLink70; clear tRLink71 tRLink72 tRLink73 tRLink74 tRLink75 tRLink76 tRLink77 tRLink78 tRLink79 tRLink80; clear tRLink81 tRLink82 tRLink83 tRLink84; [AllLink]= [ [Link1] [Link2] [Link3] [Link4] [Link5] [Link6] [Link7] [Link8] [Link9] [Link10] [Link11] [Link12] [Link13] [Link14] [Link15] [Link16] [Link17] [Link18] [Link19] [Link20] [Link21] [Link22] [Link23] [Link24] [Link25] [Link26] [Link27] [Link28] [Link29] [Link30] [Link31] [Link32] [Link33] [Link34] [Link35] [Link36] [Link37] [Link38] [Link39] [Link40] [Link41] [Link42] [Link43] [Link44] [Link45] [Link46] [Link47] [Link48] [Link49] [Link50] [Link51] [Link52] [Link53] [Link54] [Link55] [Link56] [Link57] [Link58] [Link59] [Link60] [Link61] [Link62] [Link63] [Link64] [Link65] [Link66] [Link67] [Link68] [Link69] [Link70] [Link71] [Link72] [Link73] [Link74] [Link75] [Link76] [Link77] [Link78] [Link79] [Link80] [Link81] [Link82] [Link83] [Link84] [Link85] [Link86] [Link87] [Link88] [Link89] [Link90] [Link91] [Link92] [Link93] [Link94] [Link95] [Link96] [Link97] [Link98] [Link99] [Link100] [Link101] [Link102] [Link103] [Link104] [Link105] [Link106] [Link107] [Link108] [Link109] [Link110] [Link111] [Link112] ]; %% Phan tinh toan Damage Index: [m2,n2]=size(AllLink); DIandEhofLink=zeros(n2/3,4); DIPAandEhofLink=zeros(n2/3,4); %%Plot the damage state of frame: MaxRAllLink=zeros(n2/3,1); for link=1:n2/3 MaxRAllLink(link)=max(AllLink(:,3*link-1)); end %Damage states: 1-No damage (teta=NOF1) & (link= 0.05 Link = num2str(AverageDI(i),2) ; if mod(LinkCo(i,3) , 3000) == % dam if (LinkCo(i,2) < 2000) || ( (5000 < LinkCo(i,2)) && (LinkCo(i,2) < 6000) ) || ( (10000