Xác định thành phần động của tải trọng gió dựa trên phân tích động lực học

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỞ THÀNH PHỐ HỒ CHÍ MINH - - HỒNG TRỌNG KHẢM XÁC ĐỊNH THÀNH PHẦN ĐỘNG CỦA TẢI TRỌNG GIÓ DỰA TRÊN PHÂN TÍCH ĐỘNG LỰC HỌC LUẬN VĂN THẠC SỸ XÂY DỰNG CƠNG TRÌNH DÂN DỤNG VÀ CƠNG NGHIỆP Tai Lieu Chat Luong Thành Phố Hồ Chí Minh tháng 03 năm 2016 NHẬN XÉT VỀ VIỆC HOÀN THÀNH LUẬN VĂN Xác nhận Học Viên: Hoàng Trọng Khảm, LỚP: XD3, Khóa 2013-2016 Tên đề tài luận văn: Xác Định Thành Phần Động Của Tải Trọng Gió Dựa Trên Phân Tích Động Lực Học Đủ điều kiện bảo vệ luận văn TPHCM, Ngày tháng năm 2016 TS NGUYỄN TRỌNG PHƯỚC LỜI CAM ĐOAN Tôi xin cam đoan Luận văn tơi tìm hiểu phát triển dựa vào tài liệu tham khảo trích dẫn Các thông tin, kết nghiên cứu người khác sử dụng luận văn điều trích dẫn quy định Kết Luận văn hoàn toàn trung thực chưa công bố nghiên cứu khác Tơi xin chịu hồn tồn trách nhiệm cơng việc thực Tác giả luận văn Hồng Trọng Khảm i 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 trường Đại Học Mở TPHCM nằm hệ thống luận cuối khóa nhằm trang bị cho học viên có đề tài nghiên cứu khoa học bổ ích, thực tiễn kiến thức chuyên sâu trước trường Để hoàn thành luận văn này, cố gắng nhận nhiều giúp đỡ từ thầy hướng dẫn, bạn bè người thân Tôi chân thành cảm ơn giúp đỡ quan tâm Đầu tiên tơi xin gửi lời cảm ơn đến thầy hướng dẫn TS Nguyễn Trọng Phước giúp tơi định hướng có nhìn tổng qt đề tài Trong trình làm luận văn Thầy hướng dẫn nhiệt tình cung cấp tài liệu q báu giúp tơi hồn thành luận văn thời gian sớm Tôi xin chân thành cảm ơn thầy trưởng khoa Xây Dựng trường Đại Học Mở TPHCM TS Trần Tuấn Anh quan tâm định hướng tơi đề tài có ứng dụng thực tiễn Luận văn thạc sĩ hoàn thành thời gian quy định với cố gắng thân khơng thể tránh sai sót Kính mong quý thầy cô đưa nhận xét giúp tơi có kiến thức nhìn đắn Xin trân trọng cảm ơn! ii TÓM TẮT Trong tính tốn thiết kế khung nhà cao tầng cao 40m phải xét đến thành phần động tải trọng gió Hiện Việt Nam để tính tốn thành phần gió động người thiết kế vào tiêu chuẩn Việt Nam TCVN 229:1999 TCVN 2737:1995 Tuy nhiên tiêu chuẩn hướng dẫn sử dụng cơng thức tính tốn tĩnh học tương đương gần để xác định lực thành phần động tải trọng gió gây mà khơng xét đến yếu tố xung vận tốc gió thay đổi theo thời gian ảnh hưởng đến làm việc khung nhà Vì có nhiều trường hợp xung lực gió tác động lên hệ khung nhà cao tầng nguy hiểm so với phương pháp tính tốn theo TCVN Do việc tính tốn nghiên cứu thành phần gió động theo lý thuyết Động lực học việc làm cần thiết, tài liệu tham khảo cho người thiết kế có nhìn sâu vấn đề Luận văn đưa giả thiết chuỗi số liệu áp lực gió tác dụng lên hệ kết cấu khung nhà nhiều tầng Kết tính tốn cho thấy vài trường hợp phương pháp tính tốn thành phần động tải trọng gió theo TCVN chưa phải thiên an toàn iii MỤC LỤC CHƯƠNG GIỚI THIỆU 1.1 Đặt vấn đề 1.2 Mục tiêu luận văn 1.3 Phương pháp thực 1.4 Phạm vi luận văn 1.5 Cấu trúc luận văn CHƯƠNG TỔNG QUAN 2.1 Giới thiệu 2.2 Khái niệm gió 2.3 Cách xác định tải trọng gió theo TCVN 2.4 Thực nghiệm Hầm gió 10 2.5 Tình hình nghiên cứu liên quan luận văn 12 2.6 Kết luận 19 CHƯƠNG CƠ SỞ LÝ THUYẾT 20 3.1 Giới thiệu 20 3.2 Thiết lập sơ đồ tính tốn đặc trưng động lực học 20 3.3 Mơ tả gió tải trọng dạng xung 26 3.4 Mơ tả gió tải trọng theo chuỗi thời gian 29 3.5 Phương pháp số giải toán 32 3.6 Phương trình chuyển động 32 3.7 Thuật toán Newmark 35 3.8 Xây dựng phổ phản ứng 37 3.9 Lưu đồ tính tốn 39 3.10 Kết luận chương 40 CHƯƠNG KẾT QUẢ SỐ 41 4.1 Giới thiệu 41 4.2 Bài toán lý thuyết cho hệ kết cấu chịu tải trọng gió dạng xung 46 4.3 Bài toán thực tế cho hệ kết cấu chịu tải trọng gió dạng xung 55 4.4 Bài toán lý thuyết cho hệ kết cấu với tải trọng gió dạng chuỗi thời gian 62 iv 4.5 Bài toán thực tế cho hệ kết cấu với tải trọng gió dạng chuỗi thời gian 65 4.6 Kết luận 70 CHƯƠNG KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 71 5.1 Kết luận 71 5.2 Hướng phát triển 72 Tài liệu tham khảo 73 Phụ lục 76 v DANH MỤC HÌNH VẼ Hình 2.1 Lực gió tĩnh theo TCVN Hình 2.2 Chi tiết cánh quạt mơ hình thực nghiệm Hầm gió 10 Hình 2.3 Mơ hình máy bay A380 thực nghiệm Hầm gió 11 Hình 3.1.Sơ đồ tính tốn hệ công xôn 21 Hình 3.2 Ma trận khối lượng thu gọn 22 Hình 3.3 Độ cứng tầng 23 Hình 3.4 Độ cứng tầng 23 Hình 3.5 Tải trọng xung hình sin 26 Hình 3.6 Tải trọng gió xung đơi nửa hình sin 27 Hình 3.7 Tải trọng gió xung ba nửa hình sin 27 Hình 3.8 Tải trọng xung hình chữ nhật 28 Hình 3.9 Tải trọng xung hình tam giác 29 Hình 3.10 Mơ vận tốc gió cấp 9-10 theo thời gian 30 Hình 3.11 Vận tốc gió theo thời gian độ cao 13; 33; 70; 160 feet 31 Hình 3.12 Sơ đồ thiết lập phương trình chuyển động 33 Hình 3.13 Kết cấu chịu tải trọng gió thay đổi theo thời gian 33 Hình 3.14 Sơ đồ lực tác dụng lên bậc tự 34 Hình 3.15 Hệ kết cấu chịu tải trọng nửa hình sin 37 Hình 3.16 Phổ phản ứng 38 Hình 4.1 Mơ hình cơng xơn tịa nhà 20 tầng 41 Hình 4.2 Sơ đồ tính tốn độ cứng 42 Hình 4.3 Sơ đồ mặt cao ốc văn phòng 45 Hình 4.4 Mơ hình khung không gian 45 Hình 4.5 Biểu đồ hệ số Kd tải trọng gió nửa hình sin 49 Hình 4.6 Tải trọng gió xung đơi nửa hình sin 50 vi Hình 4.7 Biểu đồ hệ số Kd tải trọng gió nửa hình sin 50 Hình 4.8 Tải trọng gió xung ba nửa hình sin 51 Hình 4.9 Biểu đồ hệ số Kd tải trọng gió nửa hình sin 51 Hình 4.10 Xung hình chữ nhật 51 Hình 4.11 Biểu đồ hệ số Kd xung hình chữ nhật 52 Hình 4.12 Tải trọng xung hình tam giác 53 Hình 4.13 Biểu đồ hệ số Kd xung hình tam giác 55 Hình 4.14 Phổ phản ứng hệ kết cấu 55 Hình 4.15 Mặt chung cư 56 Hình 4.16 Chuyển vị thành phần tĩnh gây 58 Hình 4.17 Chuyển vị thành phần tĩnh động gây 60 Hình 4.18 Biểu đồ mơment chân cột C1 60 Hình 4.19 Biểu đồ chuyển vị vị trí chân Cột C1 61 Hình 4.20 Biểu đồ mơment vị trí chân cột C1 40 Hình 4.21 Chuyển vị đỉnh cột C1 65 Hình 4.22 Chuyển vị đỉnh cột C1 66 vii DANH MỤC BẢNG BIỂU Bảng 2.1 Giá trị áp lực gió W0 Bảng 2.2 Hệ số k Bảng 4.1 Tải trọng gió tĩnh 48 Bảng 4.2 Thành phần động tải trọng gió 48 Bảng 4.3 Bảng so sánh hệ số động 56 Bảng 4.4 Thành phần tĩnh tải trọng gió cơng trình chung cư 57 Bảng 4.5 Chu kì tần số 58 Bảng 4.6 Bảng so sánh chuyển vị, môment, hệ số động 62 Bảng 4.7 Tổng lực gió tĩnh theo phương X 62 Bảng 4.8 Thành phần tĩnh tải trọng gió 63 Bảng 4.9 Thành phần động tải trọng gió 63 Bảng 4.10 So sánh hệ số động 64 Bảng 4.11 Thành phần tĩnh tải trọng gió cơng trình chung cư 65 Bảng 4.12 So sánh chuyển vị, môment, hệ số động 68 viii Phụ lục 0 0 0 0 0 -1 -1 0 0 0; 0 0 0 0 0 0 -1 -1 0 0 0; 0 0 0 0 0 0 -1 -1 0 0; 0 0 0 0 0 0 0 -1 -1 0 0; 0 0 0 0 0 0 0 -1 -1 0; 0 0 0 0 0 0 0 0 -1 -1 0; 0 0 0 0 0 0 0 0 -1 -1; 0 0 0 0 0 0 0 0 0 -1 1]; % kN/m, ma tran cung anpha0=0.19; anpha1=0.0098; C=anpha0*M+anpha1*K; % Phan tich Mode [V,D]=eig(K,M); [W,k]=sort(diag(D)); V=V(:,k) ; Factor=diag(V'*M*V); Phi=V*inv(sqrt(diag(Factor))); Omega=diag(sqrt(Phi'*K*Phi)); % Chon tan so dao dong dau tien hinh dang Mode tuong ung lambda=1; w = diag(Omega(1:lambda)); phi = Phi(:,1:lambda) ; umaxT=[]; num=1; vectorx=[]; Xác định thành phần động tải trọng gió dựa phân tích Động lực học 104 Phụ lục for j=0:.2:3 vectorx(num)=j; T=2.46;%s T1=j*T; dt=0.001; %s buoc thoi gian t=0:dt:T1; % N=length(t); vecto thoi gian % so buoc thoi gian %************************************************ alfa=0.5; beta=0.25; %dieu kien ban dau u{1}=zeros(20,1); ud{1}=zeros(20,1); udd{1}=zeros(20,1); %************************************************ for i=1:N omega1=pi/T1; %rad/s tan so luc kich thich p=[128.299248 173.01438 145.1108736 157.6247904 166.0938048 178.228512 182.9528316 187.6455504 192.3382692 195.608952 198.7374312 201.8659104 204.9943896 208.1228688 211.251348 213.4950048 215.5806576 217.6663104 219.5386578 221.363604]' ; %KN bien luc kich thich PP{i}=p*(1-(i*dt/T1)); end %PP NEWMARK Meff=M+C*alfa*dt+beta*K*dt^2; Xác định thành phần động tải trọng gió dựa phân tích Động lực học 105 Phụ lục MMM=inv(Meff); NNN=inv(K); for i=1:N-1 Peff{i}=PP{i+1}-K*u{i}-(C+K*dt)*ud{i}C*alfa*dt*udd{i}-beta*K*dt^2*udd{i}; udd{i+1}=MMM*Peff{i}; ud{i+1}=ud{i}+alfa*dt*udd{i}+alfa*dt*udd{i+1}; u{i+1}=u{i}+ud{i}*dt+beta*dt^2*udd{i}+beta*dt^2*udd{i+1}; end uT=[]; for i=1:N uT(:,i)=u{i}; end for i=1:20 umax(i)=max(abs(uT(i,:))); end umaxT(num)=umax(1,20); num=num+1; end [ut]=NNN*p; kd= umaxT./ut(20,1); plot(vectorx,kd,'-ok') grid on Phổ phản ứng Xác định thành phần động tải trọng gió dựa phân tích Động lực học 106 Phụ lục clear; clc; % tham so dong luc hoc M =900*diag([1 1 1 1 1 1 1 1 1 1 1]) ; %T mass K=10^6*[2 -1 0 0 0 0 0 0 0 0 0; -1 -1 0 0 0 0 0 0 0 0 0; -1 -1 0 0 0 0 0 0 0 0; 0 -1 -1 0 0 0 0 0 0 0 0; 0 -1 -1 0 0 0 0 0 0 0; 0 0 -1 -1 0 0 0 0 0 0 0; 0 0 -1 -1 0 0 0 0 0 0; 0 0 0 -1 -1 0 0 0 0 0 0; 0 0 0 -1 -1 0 0 0 0 0; 0 0 0 0 -1 -1 0 0 0 0 0; 0 0 0 0 -1 -1 0 0 0 0; 0 0 0 0 0 -1 -1 0 0 0 0; 0 0 0 0 0 -1 -1 0 0 0; 0 0 0 0 0 0 -1 -1 0 0 0; 0 0 0 0 0 0 -1 -1 0 0; 0 0 0 0 0 0 0 -1 -1 0 0; 0 0 0 0 0 0 0 -1 -1 0; 0 0 0 0 0 0 0 0 -1 -1 0; 0 0 0 0 0 0 0 0 -1 -1; 0 0 0 0 0 0 0 0 0 -1 1]; %kN/m Ma tra cung % Phan tich mode anpha0=0.19; Xác định thành phần động tải trọng gió dựa phân tích Động lực học 107 Phụ lục anpha1=0.0098; C=anpha0*M+anpha1*K; [V,D]=eig(K,M); [W,k]=sort(diag(D)); V=V(:,k) ; Factor=diag(V'*M*V); Phi=V*inv(sqrt(diag(Factor))); Omega=diag(sqrt(Phi'*K*Phi)); % Chon tan so dao dong dau tien hinh dang Mode tuong ung lambda=1; w = diag(Omega(1:lambda)); phi = Phi(:,1:lambda) ; umaxT=[]; num=1; vectorx=[]; for j=0:.2:3 vectorx(num)=j; T=2.46;%s T1=j*T; dt=0.001; %s t=0:dt:T1; N=length(t); % buoc thoi gian vecto thoi gian % so buoc thoi gian %************************************************ alfa=0.5; beta=0.25; %Dieu kien ban dau Xác định thành phần động tải trọng gió dựa phân tích Động lực học 108 Phụ lục u{1}=zeros(20,1); ud{1}=zeros(20,1); udd{1}=zeros(20,1); %************************************************ for i=1:N omega1=pi/T1; %rad/s tan so luc kich thich p=[128.299248 145.1108736 157.6247904 166.0938048 173.01438 178.228512 182.9528316 187.6455504 192.3382692 195.608952 198.7374312 201.8659104 204.9943896 208.1228688 211.251348 213.4950048 215.5806576 217.6663104 219.5386578 221.363604]'; %Bien luc kich thich PP{i}=p*(1-(i*dt/T1));%xung hinh tam giac end % PP NEWMARK Meff=M+beta*K*dt^2; MMM=inv(Meff); NNN=inv(K); for i=1:N-1 Peff{i}=PP{i+1}-K*u{i}-(C+K*dt)*ud{i}C*alfa*dt*udd{i}-beta*K*dt^2*udd{i}; udd{i+1}=MMM*Peff{i}; ud{i+1}=ud{i}+alfa*dt*udd{i}+alfa*dt*udd{i+1}; u{i+1}=u{i}+ud{i}*dt+beta*dt^2*udd{i}+beta*dt^2*udd{i+1}; end uT=[]; for i=1:N Xác định thành phần động tải trọng gió dựa phân tích Động lực học 109 Phụ lục uT(:,i)=u{i}; end for i=1:20 umax(i)=max(abs(uT(i,:))); end % [ut]=NNN*p; umaxT(num)=umax(1,20); num=num+1; end [ut]=NNN*p; kd= umaxT./ut(20,1); hold on plot(vectorx,kd,'-pk') grid on %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % umaxT=[]; num=1; vectorx=[]; for j=0:.2:3 vectorx(num)=j; T=2.46;%s T1=j*T; dt=0.001; %s t=0:dt:T1; N=length(t); % buoc thoi gian vecto thoi gian %ío buoc thoi gian Xác định thành phần động tải trọng gió dựa phân tích Động lực học 110 Phụ lục %************************************************ alfa=0.5; beta=0.25; %Dieu kien ban dau u{1}=zeros(20,1); ud{1}=zeros(20,1); udd{1}=zeros(20,1); %************************************************ for i=1:N omega1=pi/T1; p=[128.299248 %rad/s forcing frequency 145.1108736 157.6247904 166.0938048 173.01438 178.228512 182.9528316 187.6455504 192.3382692 195.608952 198.7374312 201.8659104 204.9943896 208.1228688 211.251348 213.4950048 215.5806576 217.6663104 219.5386578 221.363604]' ; %KN bien luc kich thich PP{i}=p*sin(omega1*i*dt); %xung hinh sin end % PP NEWMARK Meff=M+beta*K*dt^2; MMM=inv(Meff); NNN=inv(K); for i=1:N-1 Peff{i}=PP{i+1}-K*u{i}-(C+K*dt)*ud{i}C*alfa*dt*udd{i}-beta*K*dt^2*udd{i}; udd{i+1}=MMM*Peff{i}; ud{i+1}=ud{i}+alfa*dt*udd{i}+alfa*dt*udd{i+1}; Xác định thành phần động tải trọng gió dựa phân tích Động lực học 111 Phụ lục u{i+1}=u{i}+ud{i}*dt+beta*dt^2*udd{i}+beta*dt^2*udd{i+1}; end uT=[]; for i=1:N uT(:,i)=u{i}; end for i=1:20 umax(i)=max(abs(uT(i,:))); end umaxT(num)=umax(1,20); num=num+1; end [ut]=NNN*p; kd= umaxT./ut(20,1); plot(vectorx,kd,'-*k') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % umaxT=[]; num=1; vectorx=[]; for j=0:.2:3 vectorx(num)=j; T=2.46;%s T1=j*T; dt=0.001; %s buoc thoi gian Xác định thành phần động tải trọng gió dựa phân tích Động lực học 112 Phụ lục t=0:dt:T1; % N=length(t); vecto thoi gian % so buoc thoi gian %************************************************ alfa=0.5; beta=0.25; %Dieu kien ban dau u{1}=zeros(20,1); ud{1}=zeros(20,1); udd{1}=zeros(20,1); %************************************************ for i=1:N omega1=pi/T1; p=[128.299248 %rad/s forcing frequency 145.1108736 157.6247904 166.0938048 173.01438 178.228512 182.9528316 187.6455504 192.3382692 195.608952 198.7374312 201.8659104 204.9943896 208.1228688 211.251348 213.4950048 215.5806576 217.6663104 219.5386578 221.363604]' ; %KN bien luc kich thich PP{i}=p; %xung hinh chu nhat end % PP NEWMARK Meff=M+beta*K*dt^2; MMM=inv(Meff); NNN=inv(K); for i=1:N-1 Peff{i}=PP{i+1}-K*u{i}-(C+K*dt)*ud{i}C*alfa*dt*udd{i}-beta*K*dt^2*udd{i}; udd{i+1}=MMM*Peff{i}; Xác định thành phần động tải trọng gió dựa phân tích Động lực học 113 Phụ lục ud{i+1}=ud{i}+alfa*dt*udd{i}+alfa*dt*udd{i+1}; u{i+1}=u{i}+ud{i}*dt+beta*dt^2*udd{i}+beta*dt^2*udd{i+1}; end uT=[]; for i=1:N uT(:,i)=u{i}; end for i=1:20 umax(i)=max(abs(uT(i,:))); end umaxT(num)=umax(1,20); num=num+1; end [ut]=NNN*p; kd= umaxT./ut(20,1); plot(vectorx,kd,'-dk') Tai gio dang chuoi thoi gian clear; clc; %Cac tham so dong luc hoc M =900*diag([1 1 1 1 1 1 1 1 1 1 1]) ; %T, Ma tran khoi luong K=10^6*[2 -1 0 0 0 0 0 0 0 0 0; -1 -1 0 0 0 0 0 0 0 0 0; -1 -1 0 0 0 0 0 0 0 0; Xác định thành phần động tải trọng gió dựa phân tích Động lực học 114 Phụ lục 0 -1 -1 0 0 0 0 0 0 0 0; 0 -1 -1 0 0 0 0 0 0 0; 0 0 -1 -1 0 0 0 0 0 0 0; 0 0 -1 -1 0 0 0 0 0 0; 0 0 0 -1 -1 0 0 0 0 0 0; 0 0 0 -1 -1 0 0 0 0 0; 0 0 0 0 -1 -1 0 0 0 0 0; 0 0 0 0 -1 -1 0 0 0 0; 0 0 0 0 0 -1 -1 0 0 0 0; 0 0 0 0 0 -1 -1 0 0 0; 0 0 0 0 0 0 -1 -1 0 0 0; 0 0 0 0 0 0 -1 -1 0 0; 0 0 0 0 0 0 0 -1 -1 0 0; 0 0 0 0 0 0 0 -1 -1 0; 0 0 0 0 0 0 0 0 -1 -1 0; 0 0 0 0 0 0 0 0 -1 -1; 0 0 0 0 0 0 0 0 0 -1 1]; %kN/m, Ma tran cung anpha0=0.19; anpha1=0.0098; C=anpha0*M+anpha1*K; p=[48.88876608 65.9275848 55.29488026 60.06334118 63.29048141 67.91444352 69.71465794 71.50283078 73.29100363 74.53730592 75.72942115 76.92153638 78.11365162 79.30576685 80.49788208 81.35283341 82.1475769 82.94232038 83.65578329 84.35118384]'; % Phan tich Mode Xác định thành phần động tải trọng gió dựa phân tích Động lực học 115 Phụ lục [V,D]=eig(K,M); [W,k]=sort(diag(D)); V=V(:,k) ; Factor=diag(V'*M*V); Phi=V*inv(sqrt(diag(Factor))); Omega=diag(sqrt(Phi'*K*Phi)); % Chon tan so dao dong dau tien hinh dang Mode tuong ung lambda=1; w = diag(Omega(1:lambda)); phi = Phi(:,1:lambda) ; load Ewing.txt -ascii data=Ewing; NNN=inv(K); umaxT=[]; num=1; vectorx=[]; a=zeros(20,1); b=zeros(20,1); c=zeros(20,1); [ut]=NNN*p; T2=90; dt=0.01; %s t=0:dt:T2; N=length(t); buoc thoi gian % vecto thoi gian % so buoc thoi gian %************************************************ alfa=0.5; Xác định thành phần động tải trọng gió dựa phân tích Động lực học 116 Phụ lục beta=0.25; %Dieu kien ban dau u{1}=a; ud{1}=b; udd{1}=c; %************************************************ for i=1:N PP{i}=data(i,:)'; end % PP NEWMARK Meff=M+beta*K*dt^2; MMM=inv(Meff); for i=1:N-1 Peff{i}=PP{i+1}-K*u{i}-(C+K*dt)*ud{i}C*alfa*dt*udd{i}-beta*K*dt^2*udd{i}; udd{i+1}=MMM*Peff{i}; ud{i+1}=ud{i}+alfa*dt*udd{i}+alfa*dt*udd{i+1}; u{i+1}=u{i}+ud{i}*dt+beta*dt^2*udd{i}+beta*dt^2*udd{i+1}; end uT=[]; for i=1:N uT(:,i)=u{i}; end % -///// for i=1:20 Xác định thành phần động tải trọng gió dựa phân tích Động lực học 117 Phụ lục umax(i)=max(abs(uT(i,:))); end % -///// for i=1:20 kd(i)= umax(i)./ut(i,1); end kd' kd(20) %*************************************************** %VE DO THI uT kham=uT'; stepti=size(uT); uTT=uT(20,:)'; for i=1:stepti(2) xii(i)=i./100; end plot(xii,uT(20,:)) grid on Xác định thành phần động tải trọng gió dựa phân tích Động lực học 118