Giảm chấn nhà cao tầng do hồ nước mái đặt trên cao su lõi chì

ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐỖ THỊ NGỌC TAM GIẢM CHẤN NHÀ CAO TẦNG DO HỒ NƯỚC MÁI ĐẶT TRÊN CAO SU LÕI CHÌ Chuyên ngành: Xây dựng dân dụng - công nghiệp Mã số ngành: 60 58 20 LUẬN VĂN THẠC SĨ Thành phố Hồ Chí Minh, tháng 09 năm 2011 CƠNG TRÌNH ĐƯỢC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH Cán hướng dẫn khoa học: PGS TS ĐỖ KIẾN QUỐC (Ghi rõ họ, tên, học hàm, học vị, chữ ký) Cán chấm nhận xét 1: (Ghi rõ họ, tên, học hàm, học vị, chữ ký) Cán chấm nhận xét 1: (Ghi rõ họ, tên, học hàm, học vị, chữ ký) Luận văn thạc sĩ bảo vệ tại: HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA TP.HỒ CHÍ MINH Ngày ……….tháng ………năm ……… TRƯỜNG ĐẠI HỌC BÁCH KHOA CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM KHOA KỸ THUẬT XÂY DỰNG Độc lập – Tự – Hạnh Phúc Tp HCM, ngày ……… tháng ……… Năm 2011 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: ĐỖ THỊ NGỌC TAM Giới tính: Nữ Ngày, tháng, năm sinh: 26 - 07 - 1985 Nơi sinh: Sông Bé Chuyên ngành: Xây Dựng Dân Dụng – Công Nghiệp MSHV: 09210206 TÊN ĐỀ TÀI: GIẢM CHẤN NHÀ CAO TẦNG DO HỒ NƯỚC MÁI ĐẶT TRÊN CAO SU LÕI CHÌ NHIỆM VỤ LUẬN VĂN: + Phân tích ứng xử động lực học kết cấu khung có hồ nước đặt cao su lõi chì chịu tải trọng động đất + Xây dựng chương trình tính tốn ngơn ngữ lập trình Matlab để so sánh ứng xử kết cấu kết cấu có khơng có hồ nước mái đặt cao su lõi chì NGÀY GIAO NHIỆM VỤ: 03 - 07 - 2010 NGÀY HOÀN THÀNH NHIỆM VỤ: 29 - 08 - 2011 CÁN BỘ HƯỚNG DẪN: PGS TS ĐỖ KIẾN QUỐC Nội dung đề cương luận văn thạc sĩ Hội Đồng Chuyên Ngành thông qua CÁN BỘ HƯỚNG DẪN (Họ tên chữ ký) PGS TS ĐỖ KIẾN QUỐC CHỦ NHIỆM BỘ MÔN QUẢN LÝ CHUYÊN NGÀNH (Họ tên chữ ký) KHOA QL CHUYÊN NGÀNH (Họ tên chữ ký) LỜI CẢM ƠN Em xin gửi lời cảm ơn chân thành sâu sắc tới PGS TS ĐỖ KIẾN QUỐC Những lời động viên, nhiệt tình hướng dẫn cộng với bảo thầy không giúp em có kiến thức khoa học mà cịn hình thành tác phong làm việc khoa học Nhân dịp này, em gửi lời tri ân sâu sắc tới Thầy Cô trực tiếp truyền đạt kiến thức khoa học, cung cấp kinh nghiệm nghiên cứu hỗ trợ nhiều tài liệu quý báu để em hoàn thành tốt luận văn Những lời cảm ơn ân tình xin gửi tới bạn học viên cao học K2009, chuyên ngành xây dựng Dân dụng – Công nghiệp trường đại học Bách khoa Tp HCM, người bạn gắn bó, động viên giúp đỡ suốt hai năm qua Lời cảm ơn cuối xin gửi tới ba mẹ anh chị gia đình tạo điều kiện tốt mặt, khích lệ tinh thần lúc để em hồn thành luận văn Vì thời gian hồn thành luận văn có hạn nên khơng thể tránh khỏi thiếu sót Mọi góp ý em xin ghi nhận cập nhật thời gian sớm để đề tài hoàn chỉnh Tp Hồ Chí Minh, ngày 29 tháng 08 năm 2011 Đỗ Thị Ngọc Tam Mục lục MỤC LỤC Danh sách hình vẽ iv Danh sách bảng biểu ix Viết tắt xiii Chương TỔNG QUAN 1.1 Giới thiệu thiết bị kháng chấn hệ điều chỉnh khối lượng _TMD 1.2 Các nghiên cứu thực với TMD 1.3 Ứng dụng TMD việc kháng chấn cho nhà cao tầng 1.4 Mục tiêu phạm vi đề tài 1.5 Tóm tắt luận văn 10 Chương GIỚI THIỆU CAO SU LÕI CHÌ 11 2.1 Khái niệm 11 2.2 Quá trình phát triển 11 2.3 Cơng thức tính toán 12 2.4 Ứng dụng thực tiễn 14 Chương PHÂN TÍCH ĐỘNG LỰC HỌC KHUNG BÊ TÔNG CỐT THÉP PHẲNG 17 3.1 Thiết lập trường chuyển vị phần tử khung phẳng 17 3.2 Các giả thiết phân tích động lực học khung phẳng 19 3.3 Phương trình cân dao động 20 3.4 Xác định ma trận tính chất hệ kết cấu 23 3.4.1 Ma trận độ cứng phần tử 23 3.4.2 Ma trận khối lượng tương thích 24 3.4.3 Ma trận cản 26 3.5 Giải hệ phương trình vi phân động lực học phương pháp tích phân số 27 3.5.1 Giới thiệu phương pháp giải 27 3.5.2 Phương pháp tích phân Newmark 28 3.5.3 Các bước tính tốn theo phương pháp gia tốc trung bình 31 i Mục lục 3.6 Phân tích ứng xử động lực học khung bê tơng cốt thép phẳng chịu tải trọng động đất 32 3.6.1 Các trường hợp phân tích động lực học cho tốn khung phẳng 32 3.6.2 Trường hợp phân tích khung phẳng khơng xét phi tuyến hình học 33 3.6.3 Trường hợp phân tích khung phẳng có xét phi tuyến hình học 33 3.7 Sơ đồ thuật tốn phân tích động học khung bê tông cốt thép phẳng chịu tải trọng động đất 35 3.8 Đáp ứng lượng kết cấu có hệ cản bị động TMD 37 Chương CƠ SỞ THIẾT KẾ VÀ GIỚI THIỆU CHƯƠNG TRÌNH PHÂN TÍCH 38 4.1 Cơ sở thiết kế TMD 38 4.2 Tính tốn độ cứng liên kết hồ nước khung nhà 40 4.3 Các trường hợp giả định để tính độ cứng k 40 4.3.1 Sóng đơn sắc 40 4.3.2 Xem liệu động đất dãy sóng hẹp 40 4.3.3 Xem sóng động đất White noise 40 4.4 Chương trình phân tích 42 4.4.1 Các tập tin sử dụng chương trình 42 4.4.2 Phạm vi ứng dụng 43 Chương CÁC VÍ DỤ MINH HỌA 44 5.1 Số liệu trận động đất sử dụng để phân tích kết cấu 44 5.2 Kiểm chứng chương trình LEADRUBBER 50 5.2.1 Kiểm chứng tốn phân tích khung phẳng dao động tự 50 5.2.2 Kiểm chứng toán phân tích khung chịu tải trọng động đất 53 5.2.3 Kiểm chứng tốn phân tích khung có hồ nước mái đặt cao su lõi chì dao động tự 54 5.2.4 Kiểm chứng tốn phân tích khung có hồ nước mái đặt cao su lõi chì chịu tải trọng động đất 56 5.3 Phân tích động lực học kết cấu khung phẳng 10 tầng 58 ii Mục lục 5.3.1 Phân tích khung dao động tự 59 5.3.2 Phân tích khung chịu tải trọng động đất 62 5.4 Phân tích động lực học kết cấu khung phẳng 20 tầng 73 5.4.1 Phân tích khung dao động tự 75 5.4.2 Phân tích khung chịu tải trọng động đất 77 5.5 Xét ảnh hưởng phi tuyến hình học 89 5.5.1 Dao động tự 89 5.5.2 Chịu tải trọng động đất Kobe 90 5.6 Xét ảnh hưởng tỉ số cản liên kết vào kết phân tích 91 5.6.1 Dao động tự 91 5.6.2 Chịu tải trọng động đất Kobe 93 5.7 Xét ảnh hưởng vị trí hồ nước đến khả giảm chấn nhà 95 5.7.1 Dao động tự 95 5.7.2 Chịu tải trọng động đất Kobe 96 Chương KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 97 6.1 Kết luận 97 6.2 Hướng phát triển luận văn 98 Tài liệu tham khảo 99 Phụ lục 102 iii Danh sách hình vẽ DANH SÁCH CÁC HÌNH VẼ Hình 1.1 Điều khiển bị động với Tuned Mass Dampers Hình 1.2 Hệ hấp thu dao động Frahm Hình 1.3 Biến dạng TMD – Thí nghiệm bàn rung Hình 1.4 TMD Centerpoint Tower (Kwok and MacDonald, 1987) Hình 1.5 TMD John Hancock Tower ( Petersen, 1980,1981) Hình 1.6 TMD Citicorp Center ( Petersen, 1980,1981) Hình 1.7 TMD Chiba Port Tower (Ohtake et al., 1992) Hình 1.8 TMD Taipei 101 Hình 2.1 Vật liệu chịu lực đàn hồi 11 Hình 2.2 Cao su lõi chì 12 Hình 2.3 Mơ hình LRB – phản ứng song tuyến tính 13 Hình 2.4 Thiết bị dùng The Toushin 24 Ohmori Building 15 Hình 2.5 Thiết bị dùng Bridgestone Toranomon Building 15 Hình 3.1 Phần tử khung phẳng 17 Hình 3.2 Mơ hình thực tế sơ đồ tính 20 Hình 3.3 Sơ đồ khối tốn phân tích động lực học khung BTCT chịu tải trọng động đất 36 Hình 4.1 Mơ hình đơn giản hóa cho thiết bị TMD 38 Hình 5.1 Đồ thị gia tốc trận động đất Elcentro 1940 46 Hình 5.2 Cường độ lượng phổ trận động đất Elcentro 1940 46 Hình 5.3 Đồ thị gia tốc trận động đất Hachinohe 1968 47 Hình 5.4 Cường độ lượng phổ trận động đất Hachinohe 1968 47 Hình 5.5 Đồ thị gia tốc trận động đất Northridge 1994 48 Hình 5.6 Cường độ lượng phổ trận động đất Northridge 1994 48 Hình 5.7 Đồ thị gia tốc trận động đất Kobe 1995 49 Hình 5.8 Cường độ lượng phổ trận động đất Kobe 1995 49 Hình 5.9 Sơ đồ hình học khung tầng – nhịp 50 iv Danh sách hình vẽ Hình 5.10 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp khơng có hồ nước đặt cao su lõi chì dao động tự theo chương trình LEADRUBBER 52 Hình 5.11 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp khơng có hồ nước đặt cao su lõi chì dao động tự theo chương trình Sap2000 52 Hình 5.12 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp hồ nước đặt cao su lõi chì chịu tải trọng động đất Elcentro theo chương trình LEADRUBBER 53 Hình 5.13 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp khơng có hồ nước đặt cao su lõi chì chịu tải trọng động đất Elcentro theo chương trình Sap2000 54 Hình 5.14 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp có hồ nước mái đặt cao su lõi chì dao động tự theo chương trình LEADRUBBER 55 Hình 5.15 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp có hồ nước mái đặt cao su lõi chì dao động tự theo chương trình Sap2000 56 Hình 5.16 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp có hồ nước mái đặt cao su lõi chì chịu tải trọng động đất Elcentro theo chương trình LEADRUBBER 57 Hình 5.17 Đồ thị chuyển vị ngang nút đỉnh khung tầng – nhịp có hồ nước mái đặt cao su lõi chì chịu tải trọng động đất Elcentro theo chương trình Sap2000 57 Hình 5.18 Sơ đồ hình học khung 10 tầng – nhịp 59 Hình 5.19 Đồ thị so sánh chuyển vị ngang nút đỉnh khung 10 tầng – nhịp dao động tự 60 Hình 5.20 Đồ thị so sánh vận tốc ngang nút đỉnh khung 10 tầng – nhịp dao động tự 60 Hình 5.21 Độ giảm phản ứng kết cấu khung 10 tầng – nhịp dao động tự 61 v Danh sách hình vẽ Hình 5.22 Đồ thị so sánh chuyển vị ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Elcentro ứng với TH1 62 Hình 5.23 Đồ thị so sánh vận tốc ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Elcentro ứng với TH1 62 Hình 5.24 Độ giảm phản ứng kết cấu khung 10 tầng – nhịp chịu tải trọng động đất Elcentro 64 Hình 5.25 Đồ thị so sánh chuyển vị ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Hachinohe ứng với TH1 65 Hình 5.26 Đồ thị so sánh vận tốc ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Hachinohe ứng với TH1 65 Hình 5.27 Độ giảm phản ứng kết cấu khung 10 tầng – nhịp chịu tải trọng động đất Hachinohe 67 Hình 5.28 Đồ thị so sánh chuyển vị ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Northridge ứng với TH1 68 Hình 5.29 Đồ thị so sánh vận tốc ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Northridge ứng với TH1 68 Hình 5.30 Độ giảm phản ứng kết cấu khung 10 tầng – nhịp chịu tải trọng động đất Northridge 70 Hình 5.31 Đồ thị so sánh chuyển vị ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Kobe ứng với TH1 71 Hình 5.32 Đồ thị so sánh vận tốc ngang nút đỉnh khung 10 tầng – nhịp chịu tải trọng động đất Kobe ứng với TH1 71 Hình 5.33 Độ giảm phản ứng kết cấu khung 10 tầng – nhịp chịu tải trọng động đất Kobe 73 Hình 5.34 Sơ đồ hình học khung 20 tầng – nhịp 74 Hình 5.35 Đồ thị so sánh chuyển vị ngang nút đỉnh khung 20 tầng – nhịp dao động tự 75 Hình 5.36 Đồ thị so sánh vận tốc ngang nút đỉnh khung 20 tầng – nhịp dao động tự 75 vi Phụ lục ti(j)=(j-1)*dt; if Loaitaitrong==1 %Truong hop phan tich tai dong dat for i=1:sonut P(3*i-2)=loading(3*i-2,j)+Pn(3*i-2); P(3*i-1)=Pn(3*i-1); P(3*i)=Pn(3*i); end P(3*sonut+1)=loading(3*sonut+1,j)+Pn(3*sonut+1); else P=zeros(1,3*sonut+1); end P=double(P); Nenmatran; %Function xac dinh cac ma tran sau gan dieu kien rang buoc PP=double(PP'); Tanso; %Function phan tich tan so va mode dao dong TichphanNewmark; %Function giai he phuong trinh vi phan dong luc hoc bang phuong phap tich phan Newmark ut=double(ut); vt=double(vt); at=double(at); %******* Xac dinh vecto chuyen vi nut tong the tai thoi diem ti ******* for ij=1:(3*sonut+1) q(ij)=0; end dem=0; for i=1:(sonut) if Ux(i)~=0 %Bac tu duoc phep chuyen vi dem=dem+1; q(3*i-2)=ut(3*dem-2); q(3*i-1)=ut(3*dem-1); q(3*i)=ut(3*dem); v(3*i-2)=vt(3*dem-2); v(3*i-1)=vt(3*dem-1); v(3*i)=vt(3*dem); a(3*i-2)=at(3*dem-2); a(3*i-1)=at(3*dem-1); a(3*i)=at(3*dem); end end q(3*sonut+1)=ut(3*dem+1); v(3*sonut+1)=vt(3*dem+1); a(3*sonut+1)=at(3*dem+1); it=j; Uxdinh(it)=q(3*sotang-2); theo thoi gian Vxdinh(it)=v(3*sotang-2); thoi gian axdinh(it)=a(3*sotang-2); thoi gian %Chuyen vi theo phuong x cua nut dinh %van toc theo phuong x cua nut dinh theo %gia toc theo phuong x cua nut dinh theo if abs(Uxdinh(it))>abs(Uxdinhmax) Uxdinhmax=Uxdinh(it); itdinhmax=it; %Vong lap ung voi chuyen vi lon nhat end 114 Phụ lục if abs(Vxdinh(it))>abs(Vxdinhmax) Vxdinhmax=Vxdinh(it); itvdinhmax=it; %Vong lap ung voi van toc lon nhat end if abs(axdinh(it))>abs(axdinhmax) axdinhmax=axdinh(it); itadinhmax=it; %Vong lap ung voi gia toc lon nhat end if Vxdinhmax==Vxdinh(1); itvdinhmax=0 end if axdinhmax==axdinh(1); itadinhmax=0 end %*** Xac dinh moment tai cac vi tri nut dau va nut cuoi moi phan tu *** Mo=zeros(sophtu,2); qe=zeros(sophtu,6); for i=1:sophtu qe(i,1)=q(3*iFra(i)-2); qe(i,2)=q(3*iFra(i)-1); qe(i,3)=q(3*iFra(i)); qe(i,4)=q(3*jFra(i)-2); qe(i,5)=q(3*jFra(i)-1); qe(i,6)=q(3*jFra(i)); Matrantinhmomen; %Function xay dung ma tran tinh momen tu vect o chuyen vi phan tu Mo(i,1)=Se(1,:)*qe(i,:).'; Mo(i,2)=Se(2,:)*qe(i,:).'; end Mo1t(it)=Mo(1,1); %*** Tinh luc cat *** for i=1:sophtu Q(i,1)=(Mo(i,2)-Mo(i,1))/L(i)+PnQ(i,1)+Pnq(i,1); Q(i,2)=(Mo(i,2)-Mo(i,1))/L(i)-PnQ(i,2)-Pnq(i,2); end Q1t(it)=Q(1,1); %*** Xac dinh luc doc cot *** idau=1; icuoi=1; for i=1:socot N0=0; if iFra(i+1)~=jFra(i) %Chuyen truc cot, tu truc sang 2, icuoi=i; for ii=icuoi:-1:idau N1=0;N2=0; for m=socot+1:sophtu if jFra(ii)==jFra(m) N1=-Q(m,2); elseif jFra(ii)==iFra(m) N2=+Q(m,1); break 115 Phụ lục end end N(ii,2)=N1+N2+N0+(Pn(3*jFra(ii)-1)); N(ii,1)=N(ii,2); N0=N(ii,1); end idau=icuoi+1; i=idau; end end %*** Xac dinh luc doc dam *** N0=0; for i=socot+1:sophtu Q1=0;Q2=0; for k=1:socot if jFra(k)==iFra(i) Q1=-Q(k,2); elseif iFra(k)==iFra(i) Q2=+Q(k,1); break end end N(i,1)=Q1+Q2+N0+(Pn(3*iFra(i)-2)); N(i,2)=N(i,1); N0=N(i,2); Q1=0;Q2=0; if i==sophtu break elseif iFra(i+1)~=jFra(i) N0=0; end end if abs(N(1,1))>abs(Phanlucmax) Phanlucmax=N(1,1); itPhanluc=it; end end %Ket thuc vong lap for theo thoi gian (tai dong dat loading) %************************************************************************** Uxdinhmax itdinhmax Vxdinhmax itvdinhmax axdinhmax itadinhmax Phanlucmax itPhanluc LucCatChan=max(max(Q1t),abs(min(Q1t))) MomenChan=max(max(Mo1t),abs(min(Mo1t))) figure(3); plot(ti,Uxdinh,'LineWidth',1.2); xlabel('Thoi gian t(s)'); ylabel('Chuyen vi ngang tai nut dinh (m)'); title('Do thi chuyen vi ngang nut dinh theo t Uxdinh(t) '); grid on; 116 Phụ lục figure(4); plot(ti,Vxdinh,'LineWidth',1.2); xlabel('Thoi gian t(s)'); ylabel('Van toc ngang tai nut dinh (m/s)'); title('Do thi van toc ngang nut dinh theo t Vxdinh(t) '); grid on; figure(5); plot(ti,axdinh,'LineWidth',1.2); xlabel('Thoi gian t(s)'); ylabel('Gia toc ngang tai nut dinh (m/s2)'); title('Do thi gia toc ngang nut dinh theo t axdinh(t) '); grid on; disp('************************************************************ *******') disp('* Exe1.m *' ) disp('* FUNCTION GIAI TIM CHUYEN VI VA NOI LUC KHUNG *' ) disp('* CO XET P-DELTA *') disp('* Hoc vien : DO THI NGOC TAM *' ) disp('* Huong dan khoa hoc: PGS TS DO KIEN QUOC *') disp('* Ten chuong trinh : LEADRUBBER *') disp('* Ngay hoan : 29/06/2011 *') disp('*******************************************************************' ) disp(' ') disp(' Don vi tinh toan: KN-m-Kg ') disp(' ') disp('*******************************************************************' ) Omega=[]; Omega1=[]; N=zeros(sophtu,2); Matrandocung; %Function xac dinh ma tran cung tong th e khung co/khong xet p-delta Matrankhoiluong; %Function xac dinh ma tran khoi luong Nenmatran; %Function xac dinh cac ma tran sau gan dieu kien rang buoc duoc KK, MM, PP, Uxx Tanso; %Function phan tich tan so va mode dao dong Omega1=Omega(1); %Tan so co ban cua khung theo thoi gian Matrancan; %Function xac dinh ma tran can tong the cua khung CC ut=u0*ones(3*sonutRb+1,1); %Vecto chuyen vi ban dau vt=v0*ones(3*sonutRb+1,1); %Vecto van toc ban dau at=inv(MM)*(PP'-KK*ut-CC*vt) %Vecto gia toc ban dau disp('===================================================================' ) disp('=== CAC HE SO GIAI HE PHUONG TRINH VI PHAN DONG LUC HOC ') disp(' BANG PHUONG PHAP TICH PHAN NEWMARK = ==') disp('===================================================================' ) alpha=0.25; delta=0.5; a0=1/(alpha*dt^2); a1=delta/(alpha*dt); a2=1/(alpha*dt); a3=1/(2*alpha)-1; a4=delta/alpha-1; a5=dt/2*(delta/alpha-2); 117 Phụ lục a6=dt*(1-delta); a7=delta*dt; %************************************************************************** Uxdinh=[]; Uxdinh(1)=ut(3*sotang-2); %Chuyen vi theo phuong x tai nut dinh khung Vxdinh=[]; Vxdinh(1)=vt(3*sotang-2); %Van toc theo phuong x tai nut dinh khung axdinh=[]; axdinh(1)=at(3*sotang-2); %Gia toc theo phuong x tai nut dinh khung Uxdinhmax=Uxdinh(1); Vxdinhmax=Vxdinh(1); axdinhmax=axdinh(1); Phanlucmax=0; %Phan luc lon nhat Q1t=[]; ; Mo1t=[]; Q1t(1)=0; Mo1t(1)=0; ti=zeros(tf/dt+1,1); ti(1)=dt %***************** CAC VONG LAP THEO THOI GIAN ***************************% if Loaitaitrong==1 load Eg.txt ug=Eg; if length(ug)abs(Uxdinhmax) Uxdinhmax=Uxdinh(it); itdinhmax=it; %Vong lap ung voi chuyen vi lon nhat end if abs(Vxdinh(it))>abs(Vxdinhmax) Vxdinhmax=Vxdinh(it); itvdinhmax=it; %Vong lap ung voi van toc lon nhat end if abs(axdinh(it))>abs(axdinhmax) axdinhmax=axdinh(it); 119 Phụ lục itadinhmax=it; %Vong lap ung voi gia toc lon nhat end if Vxdinhmax==Vxdinh(1); itvdinhmax=0 end if axdinhmax==axdinh(1); itadinhmax=0 end %*** Xac dinh moment tai cac vi tri nut dau va nut cuoi moi phan tu *** Mo=zeros(sophtu,2); qe=zeros(sophtu,6); for i=1:sophtu qe(i,1)=q(3*iFra(i)-2); qe(i,2)=q(3*iFra(i)-1); qe(i,3)=q(3*iFra(i)); qe(i,4)=q(3*jFra(i)-2); qe(i,5)=q(3*jFra(i)-1); qe(i,6)=q(3*jFra(i)); Matrantinhmomen; %Function xay dung ma tran tinh momen tu vecto chuyen vi phan tu Mo(i,1)=Se(1,:)*qe(i,:).'; Mo(i,2)=Se(2,:)*qe(i,:).'; end Mo1t(it)=Mo(1,1); %*** Tinh luc cat *** for i=1:sophtu Q(i,1)=(Mo(i,2)-Mo(i,1))/L(i)+PnQ(i,1)+Pnq(i,1); Q(i,2)=(Mo(i,2)-Mo(i,1))/L(i)-PnQ(i,2)-Pnq(i,2); end Q1t(it)=Q(1,1); %*** Xac dinh luc doc cot *** idau=1; icuoi=1; for i=1:socot N0=0; if iFra(i+1)~=jFra(i) %Chuyen truc cot, tu truc sang 2, icuoi=i; for ii=icuoi:-1:idau N1=0;N2=0; for m=socot+1:sophtu if jFra(ii)==jFra(m) N1=-Q(m,2); elseif jFra(ii)==iFra(m) N2=+Q(m,1); break end end N(ii,2)=N1+N2+N0+(Pn(3*jFra(ii)-1)); N(ii,1)=N(ii,2); N0=N(ii,1); end idau=icuoi+1; i=idau; end 120 Phụ lục end %*** Xac dinh luc doc dam *** N0=0; for i=socot+1:sophtu Q1=0;Q2=0; for k=1:socot if jFra(k)==iFra(i) Q1=-Q(k,2); elseif iFra(k)==iFra(i) Q2=+Q(k,1); break end end N(i,1)=Q1+Q2+N0+(Pn(3*iFra(i)-2)); N(i,2)=N(i,1); N0=N(i,2); Q1=0;Q2=0; if i==sophtu break elseif iFra(i+1)~=jFra(i) N0=0; end end if abs(N(1,1))>abs(Phanlucmax) Phanlucmax=N(1,1); itPhanluc=it; end % Kiem tra su hoi tu cua luc doc N cho bai toan phi tuyen hinh hoc, neu % thoa thi vong lap while se ket thuc GiasoN=abs(N-Ns); %Xac dinh gia so giua vong lap ke tiep dem=0; %Gan bien dem for i=1:sophtu for n=1:2 if(GiasoN(i,n)-MtDungsai(i,n))