Phân tích động lực học kết cấu khung không gian có xét đến sự tương tác giữa móng và đất nền

ĐẠ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ự - Hạnh phúc NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên: NGUYỄN HOÀNG SƠN MSHV: 1570653 Ngày, tháng, năm sinh: 25/02/1987 Nơi sinh: TP.Hồ Chí Minh Chuyên ngành: Kỹ Thuật Xây Dựng Cơng Trình Dân Dụng Và Công Nghiệp Mã số: 60.58.02.08 I TÊN ĐỀ TÀI: “PHÂN TÍCH ĐỘNG LỰC HỌC KẾT CẤU KHUNG KHƠNG GIAN CÓ XÉT ĐẾN SỰ TƯƠNG TÁC GIỮA MÓNG VÀ ĐẤT NỀN” II NHIỆM VỤ VÀ NỘI DUNG: Phân tích kết cấu khung không gian chịu tải trọng động đất có xét đến tương tác móng cọc đất Từ đó, rút kết luận ứng xử động kết cấu khung không gian chịu tải trọng động đất có xét đến ứng xử móng cọc III NGÀY GIAO NHIỆM VỤ : 17/08/2017 IV NGÀY HOÀN THÀNH NHIỆM VỤ: 17/06/2018 V CÁN BỘ HƯỚNG DẪN: PGS TS CHU QUỐC THẮNG Tp HCM, ngày CÁN BỘ HƯỚNG DẪN PGS TS CHU QUỐC THẮNG tháng năm 2018 CHỦ TỊCH HỘI ĐỒNG NGÀNH PGS.TS BÙI CÔNG THÀNH TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG LỜI CẢM ƠN Sau thời gian học tập thực luận văn, tận tình bảo thầy cơ, động viên gia đình bạn bè, tác giả hoàn thành Luận văn Thạc sĩ Kỹ thuật xây dựng cơng trình dân dụng cơng nghiệp Đó kết sau trình học tập nghiên cứu, đánh dấu hồn thành khóa học trường Đại học Bách Khoa TP.HCM Để 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 tỏ lịng biết ơn đến tập thể cá nhân dành cho giúp đỡ quý báu Tôi xin cảm ơn Ban Giám hiệu Trường Đại học Bách khoa TP HCM, phòng Đào tạo Sau Đại học quý Thầy Cô Khoa Kỹ thuật xây dựng với tri thức tâm huyết, truyền đạt kiến thức quý báu suốt thời gian học tập nghiên cứu trường Tơi xin tỏ lịng biết ơn sâu sắc đến Thầy hướng dẫn PGS.TS Chu Quốc Thắng thầy hướng dẫn Th.S Phạm Nhân Hòa người đưa gợi ý để tơi hình thành nên ý tưởng đề tài đến lúc hoàn thành Luận văn Thạc sĩ Trong suốt trình thực Luận văn, Thầy ln tận tình hướng dẫn, bảo cho nhiều nhận định đắn nghiên cứu, cách tiếp cận, giải vấn đề cách hiệu Đó tảng vơ quan trọng để tơi hồn thành Luận văn Tôi chân thành cảm ơn quý Thầy Cô thư viện trường ĐH Bách Khoa Thành Phố Hồ Chí Minh tạo điều kiện cho tơi tìm tài liệu để thực luận văn bạn học khóa ln sát cánh bên tơi ngày học tập Cuối cùng, tơi xin bày tỏ lịng biết ơn công lao to lớn Cha Mẹ gia đình tơi Tơi xin gửi lời cảm ơn đến quan, đồng nghiệp bạn bè tạo điều kiện giúp đỡ cho tơi hồn thành khóa học Tôi xin trân trọng cảm ơn Thành Phố Hồ Chí Minh, tháng 06 năm 2018 i TĨM TẮT LUẬN VĂN THẠC SĨ Luận văn trình bày mơ hình tính tốn thuật giải để phân tích ứng xử kết cấu xét đến tương tác móng đất (SSI: Soil - Structure Interaction) Kết cấu bên mơ hình khung khơng gian có nhiều bậc tự với nhiều thành phần chuyển vị khác Trong đó, chuyển vị ngang thành phần chuyển vị chủ yếu kết cấu chịu động đất Bên kết cấu móng cọc liên kết với đất xung quanh lò xo động theo ba phương thẳng ba phương xoay Phần ví dụ minh họa nêu luận văn ví dụ mẫu kết cấu thép tầng chịu tải trọng động đất ElCentro nhằm so sánh đáp ứng động lực học kết cấu xét không xét tương tác đất Phương trình chuyển động hệ giải phương pháp tích phân Newmark miền thời gian, thực ngơn ngữ lập trình MATLAB Cuối cùng, kết luận rút khác mơ hình nhằm cung cấp thơng tin hữu ích cho việc nghiên cứu thiết kế cơng trình kháng chấn ii ABSTRACT Thesis presented computational models and algorithms to analyze the behavior of the structure when considering the interaction between foundation and soil (Soil Structure Interaction) The superstructure is modeled as multi-degree of freedom systems with different displacements in multiple direction, in which horizontal displacement is major displacements of structures subjeted to earthquake The pile foundations are restrained with the soil environment by dynamic springs in three translationals and three rotationals Numerical example is a 9-story steel structure as a benchmark problem accompanied with its pile foundation subjected to ElCentro loading in purpose of comparison with dynamic responses of a non-SSI structure The equation of motion of the system is solved by Newmark integration numerical method in the time domain and with the help of MATLAB code Finally, the conclusions to be drawn about the differences between models to provide useful information for the study and design of seismic-resistant buildings iii LỜI CAM ĐOAN Tơi tên Nguyễn Hồng Sơn, học viên cao học chun ngành Kỹ Thuật Xây Dựng Cơng Trình Dân Dụng Cơng Nghiệp, khố 2015 trường Đại Học Bách Khoa Thành Phố Hồ Chí Minh Tơi xin cam đoan rằng, luận văn tơi thực Các số liệu luận văn hoàn toàn trung thực chưa công bố hay sử dụng để bảo vệ học vị Các thông tin, tài liệu trích dẫn luận văn ghi rõ nguồn gốc Tơi xin chịu trách nhiệm hồn toàn kết nghiên cứu luận văn Học viên NGUYỄN HỒNG SƠN vi MỤC LỤC NHIỆM VỤ LUẬN VĂN THẠC SĨ .i LỜI CẢM ƠN i TÓM TẮT LUẬN VĂN THẠC SĨ ii ABSTRACT iii MỤC LỤC vii DANH MỤC HÌNH VẼ ix DANH MỤC BẢNG BIỂU xiv Chương GIỚI THIỆU - - 1.1 ĐẶT VẤN ĐỀ - 1.2 MỤC TIÊU CỦA LUẬN VĂN - Chương TỔNG QUAN - - 2.1 TỔNG QUAN VỀ SỰ TƯƠNG TÁC GIỮA MÓNG VÀ ĐẤT NỀN - 2.2 TÌNH HÌNH NGHIÊN CỨU NGOÀI NƯỚC - 2.3 TÌNH HÌNH NGHIÊN CỨU TRONG NƯỚC - 2.4 NỘI DUNG LUẬN VĂN - Chương CƠ SỞ LÝ THUYẾT - 11 - 3.1 PHƯƠNG TRÌNH VI PHÂN CHUYỂN ĐỘNG - 11 3.1.1 Khung không gian ngàm chân cột – FBB - 11 3.1.2 Khung không gian có xét SSI - 14 3.2 ĐỘ CỨNG ĐỘNG LỰC HỌC ĐỐI VỚI MÓNG CỌC - 17 3.2.1 Trên đất đồng - 17 3.2.2 Ví dụ xác định k (T ) ; c (T ) cho cọc xuyên qua lớp đất - 34 3.3 CÁC BƯỚC GIẢI BÀI TOÁN SSI VÀ SƠ ĐỒ KHỐI - 38 3.3.1 Các bước giải toán SSI - 38 - vii 3.3.2 Sơ đồ khối cho toán phân tích SSI - 40 Chương VÍ DỤ TÍNH TỐN - 41 - 4.1 KẾT CẤU TẦNG MÓNG CỌC - 41 4.1.1 Mơ tả tốn - 41 4.1.2 Tính tốn ban đầu - 44 4.1.3 Khảo sát ảnh hưởng bố trí cọc phân tích SSI - 44 4.1.4 Khảo sát ảnh hưởng mô đun đàn hồi đất hệ số poisson phân tích SSI - 69 4.1.5 Khảo sát tính đồng đất phân tích SSI - 82 Chương KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN - 89 - 5.1 KẾT LUẬN - 89 5.2 HƯỚNG PHÁT TRIỂN - 89 TÀI LIỆU THAM KHẢO - 90 PHỤ LỤC i Phụ lục 1: Tìm chiều dài, đường kính số lượng cọc i Phụ lục 2: Bảng kết phân tích SSI với trường hợp bố trí cọc trường hợp đất v Phụ lục 3: Mã nguồn chương trình Matlab x DANH MỤC KÝ HIỆU xxxi PHẦN LÝ LỊCH TRÍCH NGANG xxxiv viii DANH MỤC HÌNH VẼ Hình 2-1: Tháp Messeturm Frankfurt – Đức [39] - Hình 2-2: Tháp Deutsche Bank – Đức [39] - Hình 2-3: Tịa nhà Skyscraper - Brazil [40] - Hình 2-4: Tòa nhà Salesforce Tower - San Francisco [41] - Hình 2-5: Tương tác động (kinematic): (a) dao động theo phương đứng, (b) dao động theo phương ngang, (c) dao động không liên tục, dao động lắc - Hình 2-6: Tương tác quán tính (inertial) - Hình 2-7: Sóng lan truyền - Hình 3-1: Mơ hình khung khơng gian ngàm chân cột - 11 Hình 3-2: Phần tử mẫu khung khơng gian tổng qt - 11 Hình 3-3: Mơ hình khung khơng gian có xét SSI - 14 Hình 3-4: Mơ hình khung không gian xét SSI - 17 Hình 3-5: a) Dao động theo phương ngang móng; b) Dao động theo phương xoay móng - 18 Hình 3-6: Dao động xoắn quanh trục z móng - 18 Hình 3-7: Các thơng số móng cọc - 19 Hình 3-8: Giá trị L r theo [7] - 22 - Hình 3-9: Xác định S R  S Rp   đài có 3x3 cọc - 22 pile Hình 3-10: Xác định f z1 theo theo L p rp Ep Es pile Hình 3-11: Xác định f z theo L p rp Ep Es pile Hình 3-12: Xác định f z1 theo L p rp Ep Es cho cọc chịu mũi - 26 - cho cọc chịu ma sát - 26 - pile Hình 3-13: Xác định f z theo L p rp Ep Es cho cọc chịu mũi - 25 - cho cọc chịu ma sát - 27 - Hình 3-14: Dao động xoắn cọc đất - 28 ix Hình 3-15: Giá trị f1 cho cọc BTCT - 30 Hình 3-16: Giá trị f phụ thuộc vào L G , , tan   (cọc BTCT) , tần số Ep R không thứ nguyên a0 - 31 Hình 3-17: Hệ số chiều dài hoạt động cọc  - 32 Hình 3-18: Đài cọc xoắn quanh trục z - 34 Hình 3-19: Các thơng số cọc đài cọc trường hợp cọc xuyên qua lớp đất 35 Hình 3-20: Sơ đồ khối cho tốn phân tích SSI - 40 Hình 4-1: Mơ hình kết cấu khung khơng gian tầng - 42 Hình 4-2: Mặt (tầng điển hình) kết cấu khung khơng gian tầng - 43 Hình 4-3: Mặt bố trí cọc với np= cọc (3x3 cọc), S2Rp = S2Rp = - 45 Hình 4-4: Mặt bố trí cọc với np= cọc (2x2 cọc), S2Rp = S2Rp = - 45 Hình 4-5: Gia tốc trận động đất El Centro - 46 Hình 4-6: Đáp ứng chuyển vị đỉnh theo phương trục x (BTD 2155) - 47 Hình 4-7: Đáp ứng gia tốc đỉnh theo phương X (BTD 2155) - 47 Hình 4-8: Đáp ứng lực cắt chân cột tầng theo phương X - 48 Hình 4-9: Đáp ứng mơ men chân cột tầng theo phương X - 48 Hình 4-10: Đáp ứng chuyển vị tầng theo phương X - 48 Hình 4-11: Đáp ứng lực cắt mô men lớn tầng theo phương X - 49 Hình 4-12: Đáp ứng chuyển vị đỉnh theo phương Y (BTD 2157) - 49 Hình 4-13: Đáp ứng gia tốc đỉnh theo phương Y (BTD 2157) - 49 Hình 4-14: Đáp ứng lực cắt chân cột tầng theo phương Y - 50 Hình 4-15: Đáp ứng mơ men chân cột tầng theo phương Y - 50 Hình 4-16: Đáp ứng chuyển vị đỉnh theo phương X (BTD 2155) - 51 Hình 4-17: Đáp ứng gia tốc đỉnh theo phương X (Bậc tự 2155) - 51 x Hình 4-18: Đáp ứng lực cắt chân cột tầng theo phương X (Chân cột trục 2C)- 51 Hình 4-19: Đáp ứng mơ men chân cột tầng theo phương X - 52 Hình 4-20: Đáp ứng chuyển vị tầng theo phương X - 52 Hình 4-21: Đáp ứng lực cắt mô men lớn tầng theo phương X - 52 Hình 4-22: Đáp ứng chuyển vị đỉnh theo phương Y (BTD 2157) - 53 Hình 4-23: Đáp ứng gia tốc đỉnh theo phương Y (BTD 2157) - 53 Hình 4-24: Đáp ứng lực cắt chân cột tầng theo phương Y - 54 Hình 4-25: Đáp ứng mơ men chân cột tầng theo phương Y - 54 Hình 4-26: Đáp ứng chuyển vị đỉnh theo phương X (BTD 2155) - 57 Hình 4-27: Đáp ứng gia tốc đỉnh theo phương X (BTD 2155) - 57 Hình 4-28: Đáp ứng lực cắt chân cột tầng theo phương X (Chân cột trục 2C)- 58 Hình 4-29: Đáp ứng mơmen chân cột tầng theo phương X (Chân cột trục 2C)- 58 Hình 4-30: Đáp ứng chuyển vị tầng theo phương X - 58 Hình 4-31: Đáp ứng lực cắt mơ men lớn tầng theo phương X - 59 Hình 4-32: Đáp ứng chuyển vị đỉnh theo phương Y (BTD 2157) - 59 Hình 4-33: Đáp ứng gia tốc đỉnh theo phương Y (BTD 2157) - 59 Hình 4-34: Đáp ứng lực cắt chân cột tầng theo phương Y - 60 Hình 4-35: Đáp ứng mơ men chân cột tầng theo phương Y - 60 Hình 4-36: Đáp ứng chuyển vị tầng theo phương Y - 60 Hình 4-37: Đáp ứng lực cắt mô men lớn tầng theo phương Y - 61 Hình 4-38: Đáp ứng chuyển vị đỉnh theo phương X (BTD 2155) - 61 Hình 4-39: Đáp ứng gia tốc đỉnh theo phương X (BTD 2155) - 61 Hình 4-40: Đáp ứng lực cắt chân cột tầng theo phương X - 62 Hình 4-41: Đáp ứng mơ men chân cột tầng theo phương X - 62 Hình 4-42: Đáp ứng chuyển vị tầng theo phương X - 63 xi % Lx=9.146; Lz=9.146; H=3; H1=5.49; %% -n_node=(n_story+1)*(n_bayx+1)*(n_bayz+1); % n_column=n_story*(n_bayx+1)*(n_bayz+1); n_beam=(n_story+1)*((n_bayx)*(n_bayz+1)+(n_bayz)*(n_bayx+1)); n_CnB=n_column+n_beam; % n_VFD=n_story*(n_bayx*(n_bayz+1)+n_bayz*(n_bayx+1)); n_spring=(n_bayx+1)*(n_bayz+1); n_element=n_spring+n_column+n_beam; % temp1=[n_spring n_column n_beam n_VFD]' %% Grid of nodes for Structure node_matrix1=zeros((n_story+1)*(n_bayz+1),n_bayx+1); temp1=1:(n_bayx+1); for i=1:(n_story+1)*(n_bayz+1); node_matrix1(i,1:n_bayx+1)=temp1+(n_bayx+1)*(i-1); end temp11=(n_node+1):(n_node+1+n_bayx); node_spring=zeros(n_bayz+1,n_bayx+1); for i=1:(n_bayz+1); node_spring(i,1:n_bayx+1)=temp11+(i-1)*(n_bayx+1); end node_matrix=[node_spring; node_matrix1]; node_matrix_z=[]; node_matrix_z_m0=zeros((n_story+1)*(n_bayx+1),n_bayz+1); temp12=1:(n_bayx+1)*(n_story+1):(n_bayx+1)*(n_story+1)*(n_bayz+1); for i=1:n_bayx+1; Ax1=node_matrix_z_m0(1+(n_story+1)*(i-1):(n_story+1)*i,:)+(i-1); for j=1:n_story+1; Ax2=Ax1(j,:)+temp12+(n_bayx+1)*(j-1); node_matrix_z=[node_matrix_z;Ax2]; end end %% XZY Coordinate system Ex_node=mod(node_matrix-1,n_bayx+1)*Lx; Ex_node1=mod(node_matrix1-1,n_bayx+1)*Lx; Ex_node_z=zeros((n_story+1)*(n_bayx+1),n_bayz+1); for i=1:n_bayx+1 Ex_node_z(1+(n_story+1)*(i-1):(n_story+1)*i,:)=(i-1)*Lx; end % -node_spring(:,:)=0; Ez_node_spring=node_spring; for i=1:n_bayz+1; Ez_node_spring(i,:)=(i-1)*Lz; end Ez_node1=fix(node_matrix(n_bayz+2:end,:)/((n_bayx+1)*(n_story+1))0.01)*Lz; Ez_node=[Ez_node_spring;Ez_node1]; Ez_node_z=zeros((n_story+1)*(n_bayx+1),n_bayz+1); for i=1:n_bayz+1; Ez_node_z(:,i)=(i-1)*Lz; end % -Ey_node1=[]; for j=1:n_bayz+1; Ax=node_matrix(n_bayz+2+(n_story+1)*(j-1):n_bayz+2+(n_story+1)*(j1)+n_story,:); for i=1:n_story+1; Ax(i,:)=(i-1)*H; Phụ lục xx end Ey_node1=[Ey_node1;Ax]; end node_sring=zeros(n_bayz+1,n_bayx+1); Ey_node=[node_spring;Ey_node1]; Ey_node_z_m0=zeros((n_story+1)*(n_bayx+1),n_bayz+1); Ey_node_z=[]; for j=1:n_bayx+1; Ey_node_z1=Ey_node_z_m0(1+(n_story+1)*(j-1):1+(n_story+1)*(j1)+n_story,:); for i=1:n_story+1; Ey_node_z1(i,:)=(i-1)*H; end Ey_node_z=[Ey_node_z;Ey_node_z1]; end % -spring_node=[]; for i=1:n_bayz+1; A1=node_matrix(i,:)'; A2=node_matrix(n_bayz+2+(n_story+1)*(i-1):n_bayz+2+(n_story+1)*(i1)+n_story,:); A3=A2(i/i,:)'; spring_node=[spring_node;A1,A3]; end % -column_node=[]; for i=1:n_story; for j=1:n_bayz+1; Ax1=node_matrix1(i+(n_story+1)*(j-1):i+(n_story+1)*(j-1)+1,:)'; column_node=[column_node;Ax1]; end end % -beam_node=[]; for i=1:n_story+1; for j=1:n_bayz+1; Ax2=node_matrix1(i+(n_story+1)*(j-1),:); for k=1:n_bayx; Ax3=Ax2(1,k:k+1); beam_node=[beam_node;Ax3]; end end for j=1:n_bayx+1; Ax4=node_matrix_z(i+(n_story+1)*(j-1),:); for k=1:n_bayz; Ax5=Ax4(1,k:k+1); beam_node=[beam_node;Ax5]; end end end % -Ex_spring_m=Ex_node(1:n_bayz+1,:)'; Ex_spring=[Ex_spring_m(:) Ex_spring_m(:)]; Ez_spring_m1=Ez_node_spring'; Ez_spring=[Ez_spring_m1(:) Ez_spring_m1(:)]; Ey_spring=zeros(size(Ez_spring)); % -Ex_column=[]; Phụ lục xxi for i=1:n_story; for j=1:n_bayz+1; Ex_col=Ex_node1(i+(n_story+1)*(j-1):i+(n_story+1)*(j-1)+1,:)'; Ex_column=[Ex_column;Ex_col]; end end Ez_column=[]; for i=1:n_story; for j=1:n_bayz+1; Ez_col=Ez_node1(i+(n_story+1)*(j-1):i+(n_story+1)*(j-1)+1,:)'; Ez_column=[Ez_column;Ez_col]; end end % -Ey_column1=[]; for i=1:n_story; for j=1:n_bayz+1; Ey_col=Ey_node1(i+(n_story+1)*(j-1):i+(n_story+1)*(j-1)+1,:)'; Ey_column1=[Ey_column1;Ey_col]; end end Ey_column=Ey_column1+(H1-H); Ey_column(1:(n_bayx+1)*(n_bayz+1),1)=0; % -Ex_beam=[]; for i=1:n_story+1; for j=1:n_bayz+1; Ax2=Ex_node1(i+(n_story+1)*(j-1),:); for k=1:n_bayx; Ax3=Ax2(1,k:k+1); Ex_beam=[Ex_beam;Ax3]; end end for j=1:n_bayx+1; Ax4=Ex_node_z(i+(n_story+1)*(j-1),:); for k=1:n_bayz; Ax5=Ax4(1,k:k+1); Ex_beam=[Ex_beam;Ax5]; end end end % -Ez_beam=[]; for i=1:n_story+1; for j=1:n_bayz+1; Ax6=Ez_node1(i+(n_story+1)*(j-1),:); for k=1:n_bayx; Ax7=Ax6(1,k:k+1); Ez_beam=[Ez_beam;Ax7]; end end for j=1:n_bayx+1; Ax4=Ez_node_z(i+(n_story+1)*(j-1),:); for k=1:n_bayz; Ax5=Ax4(1,k:k+1); Ez_beam=[Ez_beam;Ax5]; end end end % Phụ lục xxii Ey_beam=[]; for i=1:n_story+1; for j=1:n_bayz+1; Ax2=Ey_node1(i+(n_story+1)*(j-1),:); for k=1:n_bayx; Ax3=Ax2(1,k:k+1); Ey_beam=[Ey_beam;Ax3]; end end for j=1:n_bayx+1; Ax4=Ey_node_z(i+(n_story+1)*(j-1),:); for k=1:n_bayz; Ax5=Ax4(1,k:k+1); Ey_beam=[Ey_beam;Ax5]; end end end Ey_beam=Ey_beam+(H1-H); Ey_beam(1:(n_bayz)*(n_bayx+1)+(n_bayx)*(n_bayz+1),:)=0; % -element_node=[spring_node;column_node;beam_node]; Ex=[Ex_spring;Ex_column;Ex_beam]; Ez=[Ez_spring;Ez_column;Ez_beam]; Ey=[Ey_spring;Ey_column;Ey_beam]; % -eo_column=[]; for i=1:n_column; Xs(i)=Ey_column(i,1)*(Ez_column(i,1)Ez_column(i,2))+Ey_column(i,2)*(Ez_column(i,1)Ez_column(i,1))+Ey_column(i,1)*(Ez_column(i,2)-Ez_column(i,1)); Ys(i)=Ex_column(i,1)*(Ez_column(i,2)Ez_column(i,1))+Ex_column(i,2)*(Ez_column(i,1)Ez_column(i,1))+0.01*(Ez_column(i,1)-Ez_column(i,2)); Zs(i)=Ex_column(i,1)*(Ey_column(i,1)Ey_column(i,2))+Ex_column(i,2)*(Ey_column(i,1)Ey_column(i,1))+0.01*(Ey_column(i,2)-Ey_column(i,1)); eo_column=[eo_column; Xs(i) Ys(i) Zs(i)]; end eo_beam=[]; for j=1:n_beam; Ls(j)=Ey_beam(j,1)*(0.01-Ez_beam(j,2)) + Ey_beam(j,2)*(Ez_beam(j,1)-0.01) + Ey_beam(j,1)*(Ez_beam(j,2)Ez_beam(j,1)); Ms(j)=Ex_beam(j,1)*(Ez_beam(j,2)-0.01) + Ex_beam(j,2)*(0.01Ez_beam(j,1)) + 0.01*(0.01-Ez_beam(j,2)); Ns(j)=Ex_beam(j,1)*(Ey_beam(j,1)-Ey_beam(j,2)) + Ex_beam(j,2)*(Ey_beam(j,1)-Ey_beam(j,1))+ 0.01*(Ey_beam(j,2)Ey_beam(j,1)); eo_beam=[eo_beam; Ls(j) Ms(j) Ns(j)]; end eo_CnB=[eo_column;eo_beam]; %% Index of elements Edof=zeros(n_element,12); for i=1:n_element; node_i=element_node(i,1); dof_node_i=[node_i*6-5 node_i*6-4 node_i*6-3 node_i*6-2 node_i*6-1 node_i*6]; node_j=element_node(i,2); Phụ lục xxiii dof_node_j=[node_j*6-5 node_j*6-4 node_j*6-3 node_j*6-2 node_j*6-1 node_j*6]; Edof(i,:)=[dof_node_i dof_node_j]; end Edof=[(1:n_element)' Edof]; Edof_spring=Edof(1:n_spring,:); Edof_column=Edof(n_spring+1:n_spring+n_column,:); Edof_beam=Edof(n_spring+n_column+1:n_spring+n_column+n_beam,:); % Edof_VFD=Edof(n_spring+n_column+n_beam+1:n_element,:); %% Boundary conditions bc=[]; bc1=Edof(1:(n_bayx+1)*(n_bayz+1),2:7); for i=1:(n_bayx+1)*(n_bayz+1); bc2=[bc1(i,:)' zeros(6,1)]; bc=[bc;bc2]; end bc_no_spring=[]; for i=1:(n_bayx+1)*(n_bayz+1); bc_no_spring1=[Edof_column(i,2:7)' zeros(6,1)]; bc_no_spring=[bc_no_spring;bc_no_spring1]; end function [Ke,Me,Ce]=dynamic_spring_soil_3D(ep) % [Ke,Me,Ce]=dynamic_spring_soil(ep) % -% PURPOSE % Compute element stiffness matrix % for spring (analog) element % % INPUT: ep = [k]; spring stiffness or analog quantity % % OUTPUT: Ke = [k]*{u} : stiffness matrix, dim(Ke)= 12 x 12, {u}={ux1 uy1 uz1 ux2 uy2 uz2 ux3 uy3 uz3 ux4 uy4 uz4} % ep_spring(i,:)=[E_spring v_spring kx cx ky cy kz cz kqx cqx kqy cqy kqz cqz m0_1]; k_x=ep(3); c_x=ep(4); k_y=ep(5); c_z=ep(8); k_qx=ep(9); c_qx=ep(10); k_qy=ep(11); c_qz=ep(14); c_y=ep(6); k_z=ep(7); c_qy=ep(12); k_qz=ep(13); m0_1=ep(15); % Ke=[k_x 0 0 -k_x 0 0 0 k_y 0 0 -k_y 0 0 0 k_z 0 0 -k_z 0 0 0 k_qx 0 0 -k_qx 0 0 0 k_qy 0 0 -k_qy 0 0 0 k_qz 0 0 -k_qz -k_x 0 0 k_x 0 0 0 -k_y 0 0 k_y 0 0 0 -k_z 0 0 k_z 0 0 0 -k_qx 0 0 k_qx 0 0 0 -k_qy 0 0 k_qy 0 0 0 -k_qz 0 0 k_qz]; % Me=zeros(12); alpha=1/12; Me(1,1)=m0_1; Phụ lục Me(2,2)=m0_1; Me(3,3)=m0_1; xxiv Me(4,4)=m0_1*alpha*0; Me(6,6)=m0_1*alpha*0; Me(7,7)=m0_1; Me(10,10)=m0_1*alpha*0; Me(12,12)=m0_1*alpha*0; Me(5,5)=m0_1*alpha*0; Me(8,8)=m0_1; Me(11,11)=m0_1*alpha*0; Me(9,9)=m0_1; Ce=[c_x 0 0 -c_x 0 0 0 c_y 0 0 -c_y 0 0 0 c_z 0 0 -c_z 0 0 0 c_qx 0 0 -c_qx 0 0 0 c_qy 0 0 -c_qy 0 0 0 c_qz 0 0 -c_qz -c_x 0 0 c_x 0 0 0 -c_y 0 0 c_y 0 0 0 -c_z 0 0 c_z 0 0 0 -c_qx 0 0 c_qx 0 0 0 -c_qy 0 0 c_qy 0 0 0 -c_qz 0 0 c_qz]; % -end -function [Ke,Me,kle,mle,Ce]=beam3e(ex,ey,ez,eo,ep) % Ke=beam3e(ex,ey,ez,eo,ep) % [Ke,fe]=beam3e(ex,ey,ez,eo,ep,eq) % -% PURPOSE % Calculate the stiffness matrix for a 3D elastic Bernoulli % beam element % % INPUT: ex = [x1 x2] % ey = [y1 y2] % ez = [z1 z2] node coordinates % % eo = [xz yz zz]; orientation of local z axis % % ep = [E G A Iy Iz Kv]; element properties % E: Young's modulus % G: Shear modulus % A: Cross section area % Iy: moment of inertia,local y-axis % Iz: moment of inertia,local z-axis % Kv: Saint-Venant's torsion constant % b=[ ex(2)-ex(1); ey(2)-ey(1); ez(2)-ez(1) ]; L=sqrt(b'*b); n1=b/L; lc=sqrt(eo*eo'); n3=eo/lc; % if nargin==5; eq=[0 0 0]; end qx=eq(1); qy=eq(2); qz=eq(3); qw=eq(4); % E=ep(1); Gs=ep(2); A=ep(3); Iy=ep(4); Iz=ep(5); Kv=ep(6); a=E*A/L ; b=12*E*Iz/L^3 ; c=6*E*Iz/L^2; d=12*E*Iy/L^3 ; e=6*E*Iy/L^2 ; f=Gs*Kv/L; g=2*E*Iy/L ; h=2*E*Iz/L ; kle=[ a 0 0 -a Phụ lục 0 b 0 d 0 -e c 0 0 0 -a 0 0 0 c -b 0 0 -e 0 -d -e f 0 0 -f 0 2*g 0 e g 0 2*h -c 0 0 0 a 0 0 c 0 h xxv -b 0 0 -d 0 0 -f 0 -e 0 c 0 -c e 0 g 0 h b 0 0 0 -c d e 0 -c e f 0 2*g 0 2*h]; % fle=L/2*[qx qy qz qw -1/6*qz*L 1/6*qy*L qx qy qz qw 1/6*qz*L -1/6*qy*L]'; Ix=ep(7); m=ep(8); j=140*Ix/A; k=4*L^2; n=22*L; p=13*L; q=3*L^2; mle=m*L/420*[140 0 0 70 0 0 0; 156 0 n 54 0 -p; 0 156 -n 0 54 p 0; 0 j 0 0 j/2 0; 0 -n k 0 -p -q 0; n 0 k p 0 -q; 70 0 0 140 0 0 0; 54 0 p 156 0 -n; 0 54 -p 0 156 n 0; 0 j/2 0 0 j 0; 0 p -q 0 n k 0; -p 0 -q -n 0 k]; % n2(1)=n3(2)*n1(3)-n3(3)*n1(2); n2(2)=-n1(3)*n3(1)+n1(1)*n3(3); n2(3)=n3(1)*n1(2)-n1(1)*n3(2); % An=[n1'; n2; n3]; % G=[ An zeros(3) zeros(3) zeros(3); zeros(3) An zeros(3) zeros(3); zeros(3) zeros(3) An zeros(3); zeros(3) zeros(3) zeros(3) An ]; % Ke1=G'*kle*G; fe1=G'*fle; Me1=G'*mle*G; % -if nargout==0 disp('Element stiffness matrix: '); disp(Ke1); if nargin==6 disp('Element load vector: '); disp(fe1); end return end Ke=Ke1; Me=Me1; if nargin==6 fe=fe1; end a1=0; a2=0; Ce=a1*Me+a2*Ke; % - end -function [K_bc,F_bc]=eliminate_boundary_condition(bc,K,F) K_bc=K; K_bc(bc(:,1),:)=[]; K_bc(:,bc(:,1))=[]; Phụ lục xxvi if nargin==3 F_bc=F; F_bc(bc(:,1),:)=[]; End function [T,T_column_n_beam,C_bc,C_column_n_beam_bc,a1_a2,a1_a2_column_n_beam]=det ermine_period_vector_damping_matrix (K,M,bc,K_CnB_no_spring,M_CnB_no_spring,bc_no_spring,zeta) %% Using K & M with spring to calculate T [K_bc]=eliminate_boundary_condition(bc,K); [M_bc]=eliminate_boundary_condition(bc,M); [~,eigenvalue]=eig(K_bc,M_bc); omega=diag(sqrt(eigenvalue)); omega=sort(omega,1); T=(2*pi)./omega; % Period a1=zeta*(2*omega(1)*omega(2))/(omega(1)+omega(2)); a2=zeta*(2)/(omega(1)+omega(2)); C_bc=a1*M_bc+a2*K_bc; % Global damping matrix a1_a2=[a1 a2]; % -%% Using K & M without spring and then assem spring damping matrices [K_CnB_no_spring_bc]=eliminate_boundary_condition(bc_no_spring,K_CnB_no_s pring); [M_CnB_no_spring_bc]=eliminate_boundary_condition(bc_no_spring,M_CnB_no_s pring); [~,eigenvalue]=eig(K_CnB_no_spring_bc,M_CnB_no_spring_bc); omega=diag(sqrt(eigenvalue)); omega=sort(omega,1); T_column_n_beam=(2*pi)./omega; % Period a1=zeta*(2*omega(1)*omega(2))/(omega(1)+omega(2)); a2=zeta*(2)/(omega(1)+omega(2)); C_column_n_beam_bc=a1*M_CnB_no_spring_bc+a2*K_CnB_no_spring_bc; % Global damping matrix a1_a2_column_n_beam=[a1 a2]; % Response_with_TimeNewmark % -function [D,V,A,int_CnB]=response_with_TimeNewmark_v9 (M_bc,C_bc,K_bc,tt,D_0,V_0,loading_x_bc,bc,Edof_CnB) %% Time Newmark Method with Linear acceleration method gamma=1/2; beta=1/6; %% Initial Values nt=length(tt); delta_t=tt(2)-tt(1); % Set displacements, velocities, and acceleration of the structure in advance n_dof=length(D_0); n_bc=length(bc); n_dof_bc=n_dof-n_bc; dof=1:n_dof; dof(bc(:,1))=[]; dof_bc=dof'; dof=1:n_dof; % %% Calculate the initial force D_0_bc=apply_bc(D_0,bc); V_0_bc=apply_bc(V_0,bc); a_0_bc=M_bc\(loading_x_bc(:,1)-C_bc*V_0_bc-K_bc*D_0_bc); once A_0=[zeros(n_bc,1);a_0_bc]; %% Storage of the data D=zeros(n_dof,nt); V=zeros(n_dof,nt); Phụ lục % Only used xxvii A=zeros(n_dof,nt); D(:,1)=D_0; V(:,1)=V_0; A(:,1)=A_0; D_bc(:,1)=apply_bc(D(:,1),bc); V_bc(:,1)=apply_bc(V(:,1),bc); A_bc(:,1)=apply_bc(A(:,1),bc); De_CnB(:,:,1)=extract(Edof_CnB,D(:,1)); Ve_CnB(:,:,1)=extract(Edof_CnB,V(:,1)); Ae_CnB(:,:,1)=extract(Edof_CnB,A(:,1)); delta_P_j_cir_bc=zeros(n_dof_bc,nt); %% External and Earthquake load P_j_bc=loading_x_bc; delta_P_j_bc=zeros(n_dof_bc,nt-1); for j=1:nt-1 delta_P_j_bc(:,j)=P_j_bc(:,j+1)-P_j_bc(:,j); end %% Calculations for each time step for j = 2:nt j K_cir_bc=K_bc+1/beta/delta_t^2*M_bc+gamma/beta/delta_t*C_bc; delta_P_j_cir_bc(:,j)=delta_P_j_bc(:,j-1)+(M_bc*1/beta/delta_t+ C_bc*gamma/beta)*V_bc(:,j1)+(M_bc*1/2/beta+C_bc*delta_t*(gamma/2/beta-1))*A_bc(:,j-1); delta_D_bc=K_cir_bc\delta_P_j_cir_bc(:,j); delta_V_bc=gamma/beta/delta_t*delta_D_bc-gamma/beta*V_bc(:,j-1)+(1gamma/2/beta)*delta_t*A_bc(:,j-1); D_bc(:,j)=D_bc(:,j-1)+delta_D_bc; V_bc(:,j)=V_bc(:,j-1)+delta_V_bc; D(dof_bc,j)=D_bc(:,j); V(dof_bc,j)=V_bc(:,j); D(bc(:,1),j)=bc(:,2); V(bc(:,1),j)=bc(:,2); De_CnB(:,:,j)=extract(Edof_CnB,D(:,j)); Ve_CnB(:,:,j)=extract(Edof_CnB,V(:,j)); A_bc(:,j)=M_bc\(P_j_bc(:,j)-C_bc*V_bc(:,j)-K_bc*D_bc(:,j)); % A(dof_bc,j)=A_bc(:,j); A(bc(:,1),j)=bc(:,2); Ae_CnB(:,:,j)=extract(Edof_CnB,A(:,j)); % End Calculate_k_c_1layer function [m0_1,k0,c0,kq,cq]=Calculate_k_c_1layer(E_s_p,v_s_p,n_p,S_2Rp) Properties_Pile_Soil %% CALCULATE Initial Ep_Gsp=E_p/G_s_p; Lp_Rp=L_p/R_p; [f01,f02,f0q1,f0q2,fq1,fq2,fz1,fz2,Cx1,Cx2,Sx1,Sx2]=calculation_f_for_pil e(v_s_cap,v_s_p,Ep_Gsp,Lp_Rp); %% CALCULATION kz1; cz1 % Caluate kz1; cz1 -kz1=E_p*A_p/R_p*fz1; cz1=E_p*A_p/(V_s_p)*fz2; %% CALCULATION K0; C0 % Caluate k01; c01 -k01=E_p*I_p/R_p^3*f01; c01=E_p*I_p/(R_p^2*V_s_p)*f02; % Calculate k0_g; c0_g -% Calculate alpha [Alpha_0_g]=Calculate_Alpha_x(Beta_p,S_n_p,R_p); Phụ lục xxviii k0_g=n_p*k01/sum(Alpha_0_g); c0_g=n_p*c01/sum(Alpha_0_g); % Calculate k0_cap; c0_cap -k0_cap=G_s_p*r_0_cap*(Cx1+G_s_cap*D_s_cap*Sx1/(G_s_p*r_0_cap)); c0_cap=r_0_cap^2*sqrt(G_s_p*rho_s_p)*(Cx2+sqrt(G_s_cap*rho_s_cap/ (G_s_p*rho_s_p))*D_s_cap*Sx2/r_0_cap); % Calculate k0; c0 -k0=k0_g+k0_cap; c0=c0_g+c0_cap; %% CALCULATION kq; cq % Caluate kq1; cq1 -kq1=E_p*I_p/R_p*fq1; cq1=E_p*I_p/(V_s_p)*fq2; % Caluate kq1; cq1 -k0q1=E_p*I_p/R_p^2*f0q1; c0q1=E_p*I_p/(R_p*V_s_p)*f0q2; % Calculate kq_g; cq_g -kq_g=n_p*(kq1+kz1*x_r^2+k01*Z_c^2-2*Z_c*k0q1); cq_g=n_p*(cq1+cz1*x_r^2+c01*Z_c^2-2*Z_c*c0q1); % Calculate kq_cap; cq_cap -Sq1=2.5; Sq2=1.8; delta_cap=D_s_cap/r_0_cap; kq_cap=G_s_cap*r_0_cap^2*D_s_cap*Sq1+G_s_cap*r_0_cap^2*D_s_cap* (delta_cap^2/3+(Z_c/r_0_cap)^2-delta_cap*(Z_c/r_0_cap))*Sx1; cq_cap=delta_cap*r_0_cap^4*sqrt(G_s_cap*rho_s_cap)*(Sq2+ (delta_cap^2/3+(Z_c/r_0_cap)^2-delta_cap*(Z_c/r_0_cap))*Sx2); % Calculate kq; cq -kq=kq_g+kq_cap; % N-m/rad cq=cq_g+cq_cap; % N-m-s/rad %% save Calculate_k_c_1layer % Properties_Pile_Soil %% PROPERTY E_p=3.0e10; % N/m2 Reference: TC5574: Mac BT 350 %% Properties SOIL -% Soil cap -D_s_cap=2.5; % m E_s_cap=30e6; % N/m2 v_s_cap=0.3; gama_s_cap=18e3; % N/m3 G_s_cap=E_s_cap/(2*(1+v_s_cap)); % N/m2 rho_s_cap=gama_s_cap/9.81; % kg/m3 V_s_cap=sqrt(G_s_cap/rho_s_cap); % Soil: pile in main force -E_s_p=70e6; % N/m2 v_s_p=0.25; gama_s_p=19e3; % N/m3 G_s_p=E_s_p/(2*(1+v_s_p)); % N/m2 rho_s_p=gama_s_p/9.81; % kg/m3 V_s_p=sqrt(G_s_p/rho_s_p); %% PROPERTY for Multi Layer Soil % Soil pile % H_s1=2.5; % m % E_s1_p=30e6; % v_s1=0.25; % gama_s1_p=18e3; % N/m3 % G_s1_p=E_s1_p/(2*(1+v_s1)); % N/m2 % rho_s1_p=gama_s1_p/10.00; % kg/m3 % V_s1_p=sqrt(G_s1_p/rho_s1_p); % m/s % delta_s1=0.05; % % Soil pile % E_s2_p=E_s_p; % N/m2 % v_s2=v_s_p; Phụ lục xxix % % % % % % G_s2_p=G_s_p; gama_s2_p=gama_s_p; rho_s2_p=rho_s_p; V_s2_p=V_s_p; delta_s2=0.05; % H_s2=17.5; % % % % N/m2 N/m3 kg/m3 m/s % m %% Properties PILE -n_p=9 % So luong coc S_2Rp=3 % Khoang cach tim coc if n_p==4 L_p=29.25; % m D_p=0.4; % m S_n_p=[S_2Rp*D_p; S_2Rp*D_p*sqrt(2); S_2Rp*D_p]; Beta_p=[ 0; 45; 90]; x_r=S_n_p(1)/2; y_r=S_n_p(3)/2; Z_c=1.5/2; elseif n_p==6 L_p=23.4; % m D_p=0.4; % m S_n_p=[S_2Rp*D_p; S_2Rp*D_p*2; S_2Rp*D_p*sqrt(5); S_2Rp*D_p*sqrt(2); S_2Rp*D_p]; Beta_p=[ 0; 0; 30; 45; 90]; x_r=S_n_p(1); y_r=S_n_p(5)/2; Z_c=2/2; elseif n_p==9 L_p=17.55; % m D_p=0.4; % m S_n_p=[S_2Rp*D_p; S_2Rp*D_p*2; S_2Rp*D_p*sqrt(5); S_2Rp*D_p*sqrt(2); S_2Rp*D_p; S_2Rp*D_p*2*sqrt(2); S_2Rp*D_p*sqrt(5); S_2Rp*D_p*2]; Beta_p=[0;0; 30; 45; 90; 45; 60; 90]; x_r=S_n_p(1); y_r=S_n_p(5); Z_c=2/2; end R_p=D_p/2; % sqrt(a_p*b_p/pi); % m % a_p=0.5; % m % b_p=0.5; % m A_p=pi*R_p^2; % a_p*b_p; % m2 I_p=(pi/4)*R_p^4; % m4 I_p=(pi/4)*0.17^4 %% Properties CAP rho_cap=2500; % kg/m3 a_cap=x_r*2+2*D_p; b_cap=y_r*2+2*D_p; r_0_cap=sqrt(a_cap*b_cap/pi); m0_1=a_cap*b_cap*Z_c*2*rho_cap; Phụ lục xxx DANH MỤC KÝ HIỆU rp : Bán kính cọc trịn (hoặc quy đổi từ cọc chữ nhật rp  ap bp  với ap, bp chiều dài chiều rộng tiết diện cọc hình chữ nhật) r : Bán kính mặt tiếp xúc lớp đất xét r0 : Bán kính (mặt bằng) móng hình trịn (hoặc bán kính móng hình chữ nhật quy đổi r0  a f bf  với af, bf chiều dài chiều rộng móng hình chữ nhật) r0c a p : Bán kính (mặt bằng) đài cọc hình trịn (hoặc bán kính đài hình chữ cap nhật quy đổi r0  acap bcap  với acap, bcap chiều dài chiều rộng đài cọc hình chữ nhật) S x1 ; S x : Các hệ số đề xuất M Shadlou and S Bhattacharya (2014) C x1 ; S x1 ; C x ; S x : Các hệ số đươc đề xuất Novak (1974) Lzi : Chiều cao theo phương đứng (z) tầng thứ đến tầng thứ i Lxi : Chiều dài nhịp theo phương x nhịp thứ đến tầng thứ i Liy : Chiều dài nhịp theo phương x nhịp thứ đến tầng thứ i L p : Chiều dài cọc D scap : Chiều sâu chôn đài cọc D f : Chiều sâu chôn móng D p : Đường kính cọc xr : Khoảng cách từ tâm đài cọc đến cọc (theo phương X) yr : Khoảng cách từ tâm đài cọc đến cọc (theo phương Y) Danh mục ký hiệu xxxi mi  i   n  : Khối lượng tầng thứ đến tầng thứ n  : Góc xoay móng M0 k z1 : Hằng số lò xo cho cọc theo phương đứng k01 : Hằng số lò xo cho cọc theo phương ngang k : Hằng số lò xo xoay cho cọc k 0 : Hằng số lò xo cho cọc tương quang phương ngang xoay k0 : Hằng số lò xo cho cọc (tổng đài cọc nhóm cọc) theo phương ngang k 0g : Hằng số lò xo phương ngang cho nhóm cọc k : Hằng số lị xo xoay cho cọc (tổng đài cọc nhóm cọc) k g : Hằng số lị xo xoay cho nhóm cọc  : Hệ số poisson đất c z1 : Hệ số giảm chấn cho cọc theo phương đứng c01 : Hệ số giảm chấn cho cọc theo phương ngang c : Hệ số giảm chấn xoay cho cọc c 0 : Hệ số giảm chấn cho cọc tương quang phương ngang xoay c0 : Hệ số giảm chấn cho cọc (tổng đài cọc nhóm cọc) theo phương ngang c 0g : Hệ số giảm chấn phương ngang cho nhóm cọc c : Hệ số giảm chấn xoay cho cọc (tổng đài cọc nhóm cọc) c g : Hệ số giảm chấn xoay cho nhóm cọc L r : Hệ số tương tác cọc nhóm cọc zc : Khoảng cách từ tâm đài cọc đến mép đài (theo phương z)  spile : Khối lượng riêng đất chứa cọc Danh mục ký hiệu xxxii  s : Khối lượng riêng đất chứa móng I p : Mơ men qn tính tiết điện ngang cọc Gspile : Mô đun cắt đất chứa cọc (cọc đóng vào) G scap : Mơ đun cắt đất chứa đài cọc G s : Mơ đun cắt đất chứa móng E s1 ; E s : Mô đun đàn hồi lớp đất thứ thứ E p : Mô đun đàn hồi vật liệu làm cọc a : Tần số không thứ nguyên np: Số cọc nằm nhóm số cọc đài cọc f z1; f z : Thông số độ cứng giảm chấn dao động đứng f01; f02 : Thông số độ cứng hệ số giảm chấn dao động theo phương ngang f 1; f : Thông số độ cứng hệ số giảm chấn dao động xoay f 0 1; f 0 : Thông số độ cứng giảm chấn tương quan dao động trượt dao động xoay A p : Tiết điện cọc cs : Vận tốc sống cắt đất Danh mục ký hiệu xxxiii PHẦN LÝ LỊCH TRÍCH NGANG  Họ tên: Nguyễn Hoàng Sơn  Ngày, tháng, năm sinh: 25/02/1987 Nơi sinh: TP Hồ Chí Minh  Địa liên lạc: A11/5, Ấp 1, xã Vĩnh Lộc B, Bình Chánh, TP.HCM  Điện thoại: 0933647919  Email: hoangson.cev@gmail.com QUÁ TRÌNH ĐÀO TẠO  09/2005 – 06/2010: Sinh viên ngành Kỹ Thuật Cơng trình trường Đại Học Cơng Nghệ Sài Gịn  06/2015 – Đến nay: Học viên cao học ngành Kỹ Thuật Xây dựng cơng trình dân dụng cơng nghiệp, Khoa Kỹ Thuật Xây Dựng trường Đại Học Bách Khoa TP HCM – Đại Học Quốc Gia TP HCM Q TRÌNH CƠNG TÁC  07/2010 – 09/2011: Cơng Ty TNHH MTV Dịch vụ Hàng Không Sân Bay Tân Sơn Nhất - SASCO  09/2011 – 12/2014: Nhân viên Công Ty Cổ Phần Đầu Tư Xây Dựng Hiệp Phong  01/2015 – Đến nay: Nhân viên Công Ty Liên Doanh TNHH Khu Công Nghiệp Việt Nam – Singapore - VSIP Lý lịch trích ngang ... văn phân tích ứng xử động lực học kết cấu khung không gian chịu tác dụng động đất có xét đến tương tác kết cấu với móng Các nhiệu vụ cụ thể sau: - Xác định hệ số độ cứng động lực học kết cấu. .. động lực học kết cấu chịu động đất có xét đến ảnh hưởng tương tác móng cọc” [37] Luận văn xét đến móng cọc đất đồng mà không xét tương tác qua nhiều lớp đất, chưa xét đến SSI móng nơng Luận văn tác. .. định tương tác móng với đất nền; ứng xử phức tạp nhóm cọc với đất cần thiết Để có hiểu biết xác đáp ứng động lực học kết cấu tải trọng động đất gây ra, kết cấu khung không gian bên kết cấu móng