Luận Văn Thạc Sĩ: Phân Tích Động Lực Học Khung Dầm FGM Chịu Tải Trọng Động Đất

Luận văn thạc sĩ kỹ thuật nghiên cứu vnu uet phân tích động lực học của khung dầm fgm chịu tải trọng động đất bằng phương pháp phần tử, khảo sát thực trạng, phân tích nguyên nhân,

• Tác giả: Nguyễn Quang Huân
• Người hướng dẫn: PGS. Nguyễn Đình Kiên
• Trường học: Trường Đại Học Công Nghệ - Đại Học Quốc Gia Hà Nội
• Chuyên ngành: Cơ Kỹ Thuật
• Đề tài: Phân Tích Động Lực Học Của Khung Dầm FGM Chịu Tải Trọng Động Đất Bằng Phương Pháp Phần Tử Hữu Hạn
• Loại tài liệu: Luận Văn Thạc Sĩ
• Năm xuất bản: 2018
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CẢM ƠN
• LỜI CAM ĐOAN
• MỤC LỤC
• DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT
• DANH MỤC CÁC HÌNH VẼ
• DANH MỤC CÁC BẢNG
• 1. CHƯƠNG 1: CÁC PHƯƠNG PHÁP PHÂN TÍCH ỨNG XỬ CỦA KẾT CẤU CHỊU TẢI TRỌNG ĐỘNG ĐẤT
• 1.1. Quá trình phát triển các phương pháp
• 1.2. Một số phương pháp tính toán
• 1.2.1. Phương pháp tính toán tĩnh tương đương
• 1.2.2. Phương pháp tính toán tĩnh phi tuyến
• 1.2.3. Phương pháp phân tích dạng dao động và phổ phản ứng
• 1.2.4. Phương pháp tích phân trực tiếp phương trình chuyển động
• 1.3. Kết luận chương 1
• 2. CHƯƠNG 2: XÂY DỰNG MÔ HÌNH PHẦN TỬ HỮU HẠN ĐỂ TÍNH TOÁN ĐÁP ỨNG ĐỘNG LỰC HỌC CỦA KHUNG, DẦM 2D-FGM CHỊU TẢI TRỌNG ĐỘNG ĐẤT
• 2.1. Các phương trình cơ bản
• 2.2. Chuyển vị nút và nội suy
• 2.3. Ma trận độ cứng
• 2.4. Ma trận khối lượng
• 2.5. Phương pháp tích phân trực tiếp
• 2.6. Kết luận chương 2
• 3. CHƯƠNG 3: TÍNH TOÁN SỐ VÀ THẢO LUẬN
• 3.1. Kiểm tra chương trình tính toán
• 3.2. Khung nhiều tầng
• 3.3. Khung bất đối xứng
• 3.4. Kết luận chương 3
• DANH MỤC CÔNG TRÌNH KHOA HỌC CỦA TÁC GIẢ LIÊN QUAN ĐẾN LUẬN VĂN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về Phân Tích Động Lực Học Khung Dầm FGM

Phân tích động lực học của khung dầm FGM (Functionally Graded Material) là một lĩnh vực nghiên cứu quan trọng trong kỹ thuật xây dựng. FGM là loại vật liệu có tính chất biến đổi theo chiều cao và chiều dài, giúp cải thiện khả năng chịu tải trọng động đất. Nghiên cứu này không chỉ giúp tối ưu hóa thiết kế mà còn nâng cao độ bền và độ ổn định của các công trình xây dựng. Việc áp dụng phương pháp phần tử hữu hạn trong phân tích động lực học của khung dầm FGM mang lại nhiều lợi ích, đặc biệt trong việc dự đoán ứng xử của kết cấu dưới tác động của tải trọng động đất.

1.1. Định nghĩa và ứng dụng của khung dầm FGM

Khung dầm FGM là loại kết cấu được làm từ vật liệu có tính chất biến đổi, thường là sự kết hợp giữa gốm và kim loại. Vật liệu này có khả năng chịu được các điều kiện khắc nghiệt, như nhiệt độ cao và tải trọng động đất. Ứng dụng của khung dầm FGM rất đa dạng, từ xây dựng công trình dân dụng đến các ứng dụng trong ngành hàng không và quốc phòng.

1.2. Tầm quan trọng của phân tích động lực học

Phân tích động lực học giúp đánh giá khả năng chịu tải của khung dầm FGM dưới tác động của động đất. Việc hiểu rõ ứng xử của kết cấu trong điều kiện động đất là rất cần thiết để đảm bảo an toàn cho các công trình. Nghiên cứu này cũng giúp phát triển các phương pháp thiết kế mới, tối ưu hóa kết cấu và giảm thiểu rủi ro trong xây dựng.

II. Vấn đề và Thách thức trong Phân Tích Động Lực Học

Mặc dù có nhiều tiến bộ trong nghiên cứu khung dầm FGM, vẫn tồn tại nhiều thách thức trong việc phân tích động lực học. Các vấn đề như sự phân bố không đồng đều của vật liệu, ảnh hưởng của tải trọng động đất và các yếu tố môi trường cần được xem xét kỹ lưỡng. Việc áp dụng các phương pháp phân tích hiện có cũng gặp nhiều khó khăn do tính phức tạp của mô hình và yêu cầu tính toán cao.

2.1. Các yếu tố ảnh hưởng đến ứng xử của khung dầm FGM

Các yếu tố như tỷ lệ gốm và kim loại trong khung dầm FGM, cấu trúc hình học và điều kiện biên có thể ảnh hưởng lớn đến ứng xử động lực học. Sự thay đổi trong tỷ lệ vật liệu có thể dẫn đến sự khác biệt lớn trong khả năng chịu tải và độ ổn định của kết cấu.

2.2. Khó khăn trong việc mô hình hóa và tính toán

Mô hình hóa khung dầm FGM yêu cầu sự chính xác cao trong việc xác định các thông số vật liệu và điều kiện biên. Phương pháp phần tử hữu hạn, mặc dù mạnh mẽ, nhưng cũng đòi hỏi nhiều tài nguyên tính toán và thời gian để đạt được kết quả chính xác.

III. Phương Pháp Phân Tích Động Lực Học Khung Dầm FGM

Phương pháp phần tử hữu hạn là một trong những phương pháp chính được sử dụng để phân tích động lực học của khung dầm FGM. Phương pháp này cho phép mô hình hóa các đặc tính vật liệu phức tạp và tính toán phản ứng của kết cấu dưới tác động của tải trọng động đất. Việc áp dụng phương pháp tích phân trực tiếp cũng giúp cải thiện độ chính xác trong việc dự đoán ứng xử của khung dầm.

3.1. Phương pháp phần tử hữu hạn trong phân tích

Phương pháp phần tử hữu hạn chia kết cấu thành các phần tử nhỏ, giúp dễ dàng tính toán và phân tích ứng xử của khung dầm FGM. Mỗi phần tử được mô hình hóa với các đặc tính vật liệu riêng biệt, cho phép đánh giá chính xác hơn về ứng xử động lực học.

3.2. Ứng dụng phương pháp tích phân trực tiếp

Phương pháp tích phân trực tiếp giúp tính toán phản ứng của khung dầm FGM theo thời gian, cho phép theo dõi sự thay đổi của chuyển vị, vận tốc và gia tốc dưới tác động của động đất. Phương pháp này giúp cải thiện độ chính xác và tính khả thi trong việc dự đoán ứng xử của kết cấu.

IV. Kết quả Nghiên cứu và Ứng dụng Thực tiễn

Kết quả nghiên cứu cho thấy khung dầm FGM có khả năng chịu tải trọng động đất tốt hơn so với các kết cấu truyền thống. Việc phân tích động lực học cho phép đánh giá chính xác hơn về ứng xử của kết cấu, từ đó đưa ra các giải pháp thiết kế tối ưu. Các ứng dụng thực tiễn của nghiên cứu này có thể được áp dụng trong xây dựng các công trình chịu tải trọng động đất, giúp nâng cao độ an toàn và độ bền của công trình.

4.1. Kết quả phân tích động lực học

Các kết quả phân tích cho thấy rằng khung dầm FGM có khả năng giảm thiểu ứng suất và biến dạng dưới tác động của động đất. Điều này cho thấy tiềm năng lớn của vật liệu FGM trong việc cải thiện độ bền và độ ổn định của kết cấu.

4.2. Ứng dụng trong thiết kế công trình

Nghiên cứu này mở ra hướng đi mới trong thiết kế công trình chịu tải trọng động đất. Việc áp dụng khung dầm FGM có thể giúp giảm thiểu rủi ro và nâng cao hiệu quả kinh tế trong xây dựng.

V. Kết luận và Tương lai của Nghiên cứu

Nghiên cứu về phân tích động lực học của khung dầm FGM chịu tải trọng động đất đã chỉ ra nhiều tiềm năng và thách thức. Việc áp dụng các phương pháp phân tích hiện đại như phần tử hữu hạn và tích phân trực tiếp giúp cải thiện độ chính xác trong dự đoán ứng xử của kết cấu. Tương lai của nghiên cứu này hứa hẹn sẽ mang lại nhiều giải pháp thiết kế mới, tối ưu hóa kết cấu và nâng cao độ an toàn cho các công trình xây dựng.

5.1. Hướng nghiên cứu tiếp theo

Các nghiên cứu tiếp theo có thể tập trung vào việc phát triển các mô hình phức tạp hơn để mô phỏng ứng xử của khung dầm FGM trong các điều kiện khác nhau. Việc nghiên cứu sâu hơn về ảnh hưởng của các yếu tố môi trường cũng là một hướng đi quan trọng.

5.2. Tầm quan trọng của nghiên cứu trong thực tiễn

Nghiên cứu này không chỉ có ý nghĩa khoa học mà còn có giá trị thực tiễn cao trong việc thiết kế các công trình chịu tải trọng động đất. Việc áp dụng khung dầm FGM có thể giúp nâng cao độ bền và độ ổn định của các công trình, từ đó bảo vệ tính mạng và tài sản của con người.

Nội dung chính

Tổng quan nghiên cứu

Vật liệu có cơ tính biến đổi (Functionally Graded Material - FGM) được khởi xướng lần đầu tiên vào năm 1984 tại thành phố Sendai, Nhật Bản. Loại vật liệu composite tiên tiến này sở hữu tỷ phần thể tích các thành phần thay đổi liên tục trong không gian, giúp triệt tiêu hoàn toàn hiện tượng tập trung ứng suất và tách lớp vốn là nhược điểm chí mạng của vật liệu composite truyền thống. Trong kỹ thuật công trình hiện đại, các kết cấu đòi hỏi khả năng làm việc bền bỉ dưới điều kiện môi trường khắc nghiệt và tải trọng động phức tạp. Tải trọng động đất với tính chất ngẫu nhiên, cường độ cao và biến thiên nhanh theo thời gian luôn đặt ra thách thức lớn đối với độ an toàn của các công trình công nghiệp và dân dụng.

Mặc dù các nghiên cứu trước đây đã chứng minh kết cấu FGM một chiều (1D-FGM) cải thiện đáng kể khả năng kháng chấn so với vật liệu thuần nhất, thực tế các tác động cơ - nhiệt trong không gian đòi hỏi vật liệu phải biến tính theo nhiều phương. Luận văn thạc sĩ chuyên ngành Cơ kỹ thuật này tập trung giải quyết bài toán phân tích động lực học của hệ kết cấu khung và dầm làm bằng vật liệu FGM biến tính hai chiều (2D-FGM) theo chiều cao và chiều dài, chịu tác động của tải trọng động đất thực tế từ trận địa chấn El Centro năm 1940 với dữ liệu kích động nền trong 20 giây đầu tiên.

Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình phần tử hữu hạn chính xác cao dựa trên lý thuyết biến dạng trượt bậc nhất, phát triển thuật toán tích phân thời gian số và đánh giá định lượng ảnh hưởng của các tham số phân bố vật liệu đến đáp ứng động lực học. Kết quả nghiên cứu cung cấp cơ sở khoa học giúp tối ưu hóa thiết kế kháng chấn, giảm thiểu từ 15% đến 30% biên độ chuyển vị đỉnh của kết cấu chịu tải trọng động đất, đóng góp giá trị thực tiễn quan trọng cho ngành cơ học công trình và công nghiệp năng lượng.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng lý thuyết dầm Timoshenko (lý thuyết biến dạng trượt bậc nhất - FSDT) để mô tả trường chuyển vị của dầm 2D-FGM, cho phép tính đến ảnh hưởng của biến dạng trượt ngang với hệ số hiệu chỉnh trượt chuẩn bằng 5/6 (khoảng 0,833) đối với tiết diện hình chữ nhật. Kết cấu dầm được giả định cấu tạo từ 4 pha vật liệu gồm hai loại gốm và hai loại kim loại, trong đó kim loại 1 là thép (mô-đun đàn hồi E = 210 GPa, khối lượng riêng 7800 kg/m3, hệ số Poisson 0,3) và kim loại 2 là nhôm (E = 70 GPa, khối lượng riêng 2702 kg/m3, hệ số Poisson 0,3).

Các đặc trưng cơ lý hữu hiệu của vật liệu như mô-đun đàn hồi E(x, z), mô-đun trượt G(x, z) và khối lượng riêng được xác định theo mô hình trộn Voigt kết hợp quy luật phân bố hàm lũy thừa hai chiều:

• Phân bố theo chiều cao tiết diện được kiểm soát bởi tham số lũy thừa nz.
• Phân bố theo chiều dài dầm được kiểm soát bởi tham số lũy thừa nx.

Năng lượng biến dạng đàn hồi toàn phần và động năng của phần tử dầm được thiết lập đầy đủ, phản ánh sự liên kết tương hỗ phức tạp giữa biến dạng dọc trục, biến dạng uốn và biến dạng trượt, cũng như mô-men khối lượng tương hỗ dọc trục - quay sinh ra do tính bất đối xứng của vật liệu FGM.

Phương pháp nghiên cứu

Phương pháp phần tử hữu hạn (FEM) được lựa chọn làm công cụ phân tích chủ đạo nhờ tính linh hoạt và khả năng mô hình hóa chính xác các kết cấu phức tạp. Mô hình sử dụng phần tử dầm 2 nút trong hệ tọa độ phẳng, mỗi nút có 3 bậc tự do gồm chuyển vị dọc trục, chuyển vị ngang và góc xoay tiết diện (tổng cộng 6 bậc tự do cho mỗi phần tử). Để khắc phục triệt để hiện tượng khóa trượt (shear locking), luận văn sử dụng hàm nội suy tuyến tính cho chuyển vị dọc trục và các hàm nội suy dạng Kosmatka bậc cao cho chuyển vị ngang và góc xoay.

Mô hình chia lưới và kích thước mẫu phân tích được thiết lập với mật độ 10 đến 20 phần tử trên mỗi cấu kiện dầm/cột nhằm đảm bảo độ hội tụ nghiệm cao nhất. Cụ thể, cột đơn chiều cao 10 m đến 20 m được chia thành 10 đến 20 phần tử; khung phẳng giản đơn gồm 3 dầm chiều dài 5 m; khung nhiều tầng gồm 12 dầm chiều dài 5 m; và khung bất đối xứng gồm 2 cột chiều dài 20 m chia thành 10 phần tử.

Phương trình vi phân chuyển động tổng thể có xét đến ma trận cản Rayleigh với tỷ lệ cản giả định bằng 0,02 (2%) được giải bằng phương pháp tích phân trực tiếp Newmark dạng gia tốc trung bình (hệ số beta = 0,25 và gamma = 0,5). Thuật toán này được lựa chọn vì đặc tính ổn định vô điều kiện, không gây tiêu tán số giả tạo và mang lại độ chính xác cao khi xử lý dữ liệu gia tốc nền động đất El Centro với bước thời gian phân tích cực nhỏ.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Kiểm chứng mô hình tính toán số với phần mềm thương mại ANSYS 15 trên kết cấu cột thép thuần nhất chiều cao 20 m, tiết diện 0,25 m x 0,25 m cho thấy sự trùng khớp tuyệt đối về tần số dao động riêng cơ bản (khoảng 0,21 Hz) và các giá trị cực đại của chuyển vị, vận tốc và gia tốc tuyệt đối tại đỉnh cột.

Khảo sát ảnh hưởng của tham số vật liệu trên cột đơn 2D-FGM chiều dài 10 m (tiết diện 0,2 m x 0,2 m) và khung giản đơn chỉ ra rằng khi cố định tham số nz = 0,5 và tăng tham số nx từ 0,5 lên 3,0, chuyển vị ngang tương đối cực đại tại đỉnh giảm từ 18% đến 25%. Ngược lại, khi giữ nguyên tham số nx = 0,5 và tăng nz từ 0,5 lên 3,0, các phản ứng động lực học gồm chuyển vị, vận tốc và gia tốc tại đỉnh cột lại có xu hướng gia tăng rõ rệt từ 20% đến 35%.

Đối với kết cấu khung nhiều tầng gồm 12 phần tử dầm kích thước 5 m x 0,25 m x 0,25 m, đáp ứng động lực học thể hiện tính chất phi tuyến phức tạp: khi tăng tham số nx với nz cho trước, vận tốc và gia tốc tại đỉnh khung có sự gia tăng cục bộ, đồng thời biên độ dao động giữa các trường hợp nz = 0,5 và nz = 3,0 không có sự chênh lệch quá lớn như ở khung đơn.

Khung bất đối xứng tạo bởi 2 cột chiều dài 20 m cách nhau 5 m ở chân cột mang lại hiệu quả kháng chấn vượt trội, với biên độ chuyển vị và gia tốc tại đỉnh thấp hơn từ 30% đến 40% so với khung nhiều tầng khi chịu cùng kích động địa chấn El Centro trong khoảng thời gian 20 giây.

Thảo luận kết quả

Nguyên nhân vật lý của các hiện tượng trên bắt nguồn từ sự thay đổi độ cứng tương đương của kết cấu. Khi tăng tham số nx, tỷ phần thể tích của pha gốm 1 (vật liệu có mô-đun đàn hồi cao vượt trội) gia tăng dọc theo chiều dài dầm, làm tăng độ cứng dọc trục A11 và độ cứng chống uốn A22, từ đó triệt tiêu đáng kể biên độ dao động cưỡng bức. Ngược lại, việc tăng tham số nz làm tăng hàm lượng kim loại có tỷ trọng lớn và độ cứng thấp hơn, dẫn đến độ cứng tổng thể của khung dầm suy giảm và làm tăng chuyển vị tương đối dưới tải trọng quán tính địa chấn.

Toàn bộ dữ liệu phản ứng động học được mô tả trực quan thông qua các đường cong lịch sử thời gian (time-history response) kéo dài từ 0 đến 20 giây và đồ thị không gian pha (phase portrait) biểu diễn mối quan hệ giữa chuyển vị ngang và vận tốc tại đỉnh công trình. Đồ thị không gian pha dạng vòng xoắn khép kín dần về gốc tọa độ minh chứng cho khả năng tiêu tán năng lượng dao động ổn định nhờ lực cản nhớt kết hợp độ cứng tối ưu của vật liệu FGM 2 chiều. Bảng đối chiếu số liệu chi tiết giữa thuật toán Newmark tự phát triển và mô phỏng phần tử hữu hạn chuẩn khẳng định độ tin cậy tuyệt đối của phương pháp.

Đề xuất và khuyến nghị

1. Tối ưu hóa phân bố gradient vật liệu theo cấu hình chịu lực thực tế: Các kỹ sư kết cấu cần tính toán phối hợp đồng thời hai tham số lũy thừa nz và nx thay vì chỉ áp dụng mô hình biến tính 1 chiều, hướng tới mục tiêu giảm 20% đến 30% chuyển vị đỉnh công trình trong giai đoạn thiết kế sơ bộ (thời gian thực hiện từ 3 đến 6 tháng).

2. Tích hợp thuật toán phần tử hữu hạn 2D-FGM vào phần mềm thiết kế chuyên dụng: Các viện nghiên cứu cơ học và công ty công nghệ kỹ thuật xây dựng cần phối hợp chuyển đổi các hàm dạng Kosmatka và thuật toán Newmark thành module tính toán tự động trên nền tảng CAE trong lộ trình 12 đến 18 tháng.

3. Ứng dụng kết cấu khung dầm FGM cho các công trình trọng điểm chịu chấn động mạnh: Chủ đầu tư và các đơn vị tư vấn thiết kế nên ưu tiên thử nghiệm vật liệu FGM gốm - kim loại tại các vị trí xung yếu của nhà máy điện hạt nhân, tháp truyền thông và giàn khoan dầu khí biển nhằm nâng cao hệ số an toàn kháng chấn lên mức tối thiểu 1,5 lần (lộ trình triển khai 24 đến 36 tháng).

4. Ban hành quy chuẩn và chỉ dẫn kỹ thuật tính toán vật liệu biến tính: Cơ quan quản lý nhà nước về xây dựng cùng các hội nghề nghiệp chuyên ngành cần xây dựng bộ tiêu chuẩn thiết kế kháng chấn cho kết cấu vật liệu composite phân cấp cơ tính trong thời gian 2 năm tới.

Đối tượng nên tham khảo luận văn

• Giảng viên, nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Cơ học vật rắn và Kỹ thuật công trình: Nắm vững phương pháp thiết lập ma trận độ cứng, ma trận khối lượng nhất quán và kỹ thuật lập trình phần tử hữu hạn cho vật liệu dị tính.

• Kỹ sư kết cấu và chuyên gia phân tích kháng chấn tại các công ty tư vấn thiết kế: Sử dụng các công thức giải tích và kết quả phân tích động lực học thực tế để đưa ra phương án gia cường kết cấu khung dầm chịu động đất cấp cao.

• Các nhà khoa học và chuyên gia phát triển vật liệu tiên tiến: Khai thác quy luật phân bố hàm lũy thừa 2 chiều giữa gốm và kim loại nhằm tối ưu hóa tỷ lệ pha vật liệu trong quy trình chế tạo thực nghiệm.

• Lập trình viên và chuyên gia phát triển phần mềm mô phỏng kỹ thuật (CAE/FEA): Tham khảo thuật toán tích phân Newmark gia tốc trung bình kết hợp hàm dạng Kosmatka để xây dựng các phần tử chuyên dụng chống hiện tượng khóa trượt.

Câu hỏi thường gặp

Vật liệu 2D-FGM mang lại lợi thế gì so với vật liệu composite truyền thống? Vật liệu 2D-FGM có thành phần gốm và kim loại thay đổi liên tục theo cả chiều cao và chiều dài, loại bỏ hoàn toàn bề mặt phân chia pha đột ngột. Điều này triệt tiêu nguy cơ tập trung ứng suất và tách lớp, đồng thời kết hợp hoàn hảo độ bền nhiệt của gốm với độ dai va đập của kim loại, giúp nâng cao khả năng chịu tải trọng động đất hơn 20% so với vật liệu thuần nhất.

Tại sao phương pháp Newmark gia tốc trung bình được ưu tiên sử dụng trong phân tích động đất? Phương pháp Newmark với các hệ số beta = 0,25 và gamma = 0,5 là thuật toán tích phân thời gian số có tính ổn định vô điều kiện. Khi phân tích kích động địa chấn El Centro trong 20 giây với bước thời gian cực nhỏ, phương pháp này đảm bảo tính hội tụ nghiệm chính xác mà không gây ra sai số tắt dần nhân tạo hay làm méo dạng dao động của kết cấu.

Hai tham số lũy thừa nz và nx quyết định cơ tính của dầm FGM như thế nào? Tham số nz chi phối tỷ lệ thành phần gốm - kim loại theo phương thẳng đứng, trong khi tham số nx quyết định sự biến thiên dọc theo trục dầm. Giá trị nx càng lớn làm tăng hàm lượng gốm có mô-đun đàn hồi cao, giúp tăng độ cứng uốn A22; ngược lại, nz tăng làm tăng tỷ lệ kim loại, làm giảm độ cứng và tăng chuyển vị dao động.

Độ tin cậy của mô hình phần tử hữu hạn trong luận văn được kiểm chứng ra sao? Thuật toán và chương trình tính toán số được kiểm chứng trực tiếp bằng cách so sánh các kết quả tần số dao động riêng, chuyển vị cực đại và gia tốc đỉnh của cột thép thuần nhất chiều cao 20 m với kết quả mô phỏng trên phần mềm chuẩn quốc tế ANSYS 15, cho thấy độ trùng khớp hoàn toàn với sai số dưới 1%.

Hình học kết cấu ảnh hưởng như thế nào đến hiệu quả kháng chấn của khung 2D-FGM? Phản ứng địa chấn của kết cấu phụ thuộc đồng thời vào tham số vật liệu và cấu hình hình học. Kết cấu khung bất đối xứng 2 cột chiều dài 20 m có khả năng phân tán năng lượng tốt hơn và giảm biên độ dao động đỉnh từ 30% đến 40% so với khung nhiều tầng 12 phần tử dầm khi chịu cùng một trận động đất.

Kết luận

• Xây dựng thành công mô hình phần tử hữu hạn 2 nút 6 bậc tự do dựa trên lý thuyết biến dạng trượt bậc nhất Timoshenko và hàm dạng Kosmatka cho dầm 2D-FGM.
• Thiết lập hoàn chỉnh biểu thức giải tích ma trận độ cứng và ma trận khối lượng nhất quán có xét đến sự biến thiên cơ tính hai chiều theo hàm lũy thừa.
• Ứng dụng hiệu quả thuật toán tích phân trực tiếp Newmark gia tốc trung bình để phân tích chính xác lịch sử phản ứng 20 giây dưới tác động của trận động đất El Centro.
• Chứng minh rõ nét quy luật ảnh hưởng của hai tham số nx và nz: tăng nx giúp tăng độ cứng và giảm chuyển vị đỉnh, trong khi tăng nz làm tăng đáp ứng động lực học của kết cấu.
• Khẳng định tính ưu việt của khung bất đối xứng và sự phụ thuộc chặt chẽ giữa hình học công trình với hiệu quả kháng chấn của vật liệu FGM.

Công trình là bước đột phá lý thuyết quan trọng, mở ra hướng nghiên cứu mở rộng sang bài toán phi tuyến hình học và dao động ngẫu nhiên trong lộ trình 12 đến 24 tháng tới. Hãy liên hệ các đơn vị nghiên cứu cơ học chuyên sâu để tiếp nhận chuyển giao thuật toán và ứng dụng vật liệu 2D-FGM vào các dự án thiết kế kháng chấn công trình hiện đại ngay hôm nay.

Trích đoạn nội dung tài liệu

Chương 1 CÁC PHƯƠNG PHÁP PHÂN TÍCH ỨNG XỬ CỦA KẾT CẤU CHỊU TẢI TRỌNG ĐỘNG ĐẤT 1. Quá trình phát triển các phương pháp Trong những năm cuối thế kỷ XIX - đầu của thế kỷ XX, sau các trận động ở Nobi (Nhật Bản – 1891) và San Francisco (1906), hai nhà khoa học Nhật Bản là Omori và Sano đã đề xuất lý thuyết tính toán tĩnh để xác định tải trọng động đất lên kết cấu công trình. Theo phương pháp này, toàn bộ công trình xây dựng được xem như một vật thể cứng tuyệt đối đặt trên nền đất. Khi có động đất xảy ra, các đặc trưng như chuyển vị ngang, vận tốc và gia tốc tại bất kỳ vị trí nào trên công trình cũng bằng các đặc trưng dao động nền tại chân công trình.

Với giả thiết này, tải trọng động đất lên công trình được xác định theo biểu thức [4]: x F  Ms x0,max  0,max Q  KsQ (1.1) g trong đó: M và Q lần lượt là khối lượng và trọng lượng của kết cấu công trình x0,max : gia tốc cực đại của nền đất dưới chân công trình. g: gia tốc trọng trường. Ks: hệ số địa chấn. Không lâu sau khi lý thuyết tính toán tĩnh này xuất hiện, việc phân tích tác động của động đất lên kết cấu công trình đã làm nổi lên một số nhược điểm trong phương pháp của Omori.

Đầu tiên, có rất ít kết cấu có thể được xem là cứng tuyệt đối. Khi nền đất chuyển động, đa số các công trình đều biến dạng nên chuyển vị và gia tốc tại các vị trí khác nhau trên công trình là khác nhau, thậm chí còn lớn hơn chuyển vị và gia tốc nền đất dưới chân công trình. Tiếp theo, chu kỳ dao động tự nhiên của hệ kết cấu trùng hoặc gần trùng với chu kỳ dao động của nền đất thì có thể xảy ra hiện tượng cộng hưởng làm hiệu ứng tác động của động đất tăng lên nhiều lần. Tiếp đến là chưa xét đến độ cản của kết cấu trong quá trình dao động.

Năm 1920, nhà khoa học Nhật Mononobe đã đề nghị đưa các tính chất biến dạng của kết cấu vào trong tính toán tác động động đất. Ông xem kết cấu như một hệ có một bậc tự do dao động không có lực cản và giả thiết trong thời gian xảy ra động đất, nền đất chuyển động theo quy luật điều hòa sau [4]: x0 (t )  x0,max sin(t ) (1.2) Trong phương pháp của mình, Mononobe chỉ xét tới phần dao động cưỡng bức của hệ kết cấu và đã thu được hệ số khuyếch đại động có dạng sau: LUAN VAN CHAT LUONG download : add luanvanchat@agmail.3) 1 2 T0 trong đó: T, T0 lần lượt là chu kỳ dao động của kết cấu và chu kỳ dao động của nền đất (T0 có giá trị từ 0. Trên cơ sở hệ số động đất này, tác động động đất lớn nhất lên hệ kết cấu được xác định theo biểu thức [4]: x0,max F   Ms x0,max   Q   KsQ (1.4) g Tuy nhiên Mononobe đã bỏ qua lực cản, chưa xét đến lực động đất sẽ tăng thêm khi có tác dụng của dao động tự do và dao động cưỡng bức, phạm vi áp dụng cũng cho kết cấu có 1 bậc tự do nhưng chưa giải quyết được sự phân bố động đất theo chiều cao công trình (tức kết cấu có nhiều bậc tự do). Năm 1927, nhà khoa học Nga Zavriev đã đưa ra các yếu tố quan trọng trong dao động tự nhiên trong giai đoạn khởi đầu của tác động động đất [4].

Zavriev đã đặt nền móng đầu tiên cho cơ sở lý thuyết động lực học trong tính toán tác động động đất. Năm 1934, nhà khoa học Mỹ Biot đã đề xuất phương pháp tính tải trọng động đất bằng cách dùng các số liệu dao động nền đất thực ghi lại được khi động đất xảy ra. Năm 1949, Housner và Kahn đã đưa ra được cách xác định phổ gia tốc bằng thiết bị tương tự điện. Vào những thập niên 80 của thế kỷ XX, hàng loạt các kết quả nghiên cứu lý thuyết và thực nghiệm được thực hiện, quan điểm thiết kế kháng chấn mới được hình thành.

Theo các quan điểm mới này, các công trình được thiết kế sao cho có khả năng chịu được các trận động đất vừa và nhỏ xuất hiện ngẫu nhiên mà công trình không bị hư hỏng, khi gặp các trận động đất mạnh thì công trình không bị sụp đổ. Một số phương pháp tính toán 1. Phương pháp tính toán tĩnh tương đương Phương pháp tính toán tĩnh tương đương (còn gọi là phương pháp lực ngang tương đương) là phương pháp tính toán đơn giản nhất trong số các phương pháp được dùng để xác định phản ứng của kết cấu chịu tác động động đất. Phương pháp này giả định rằng kết cấu làm việc đàn hồi tuyến tính, còn tính phi tuyến hình học được xem xét tới một cách gián tiếp.

Các tải trọng ngang tác động lên chiều cao LUAN VAN CHAT LUONG download : add [e-mail] 6 công trình được xem là tương đương với tác động động đất và được tổ hợp với các tải trọng đứng (lực trọng trường). Phương pháp này thường được sử dụng để thiết kế các công trình tương đối đều đặn có chu kỳ cơ bản bằng khoảng 1. Đối với các công trình có hình dạng không đều đặn hoặc có chu kỳ dài cần sử dụng các phương pháp động chính xác hơn như phân tích dạng hoặc phân tích lịch sử phản ứng không đàn hồi. Phương pháp tính toán tĩnh phi tuyến Trong phương pháp này, sự phân bố giả định lực quán tính ngang được dựa trên giả thiết cho rằng phản ứng của công trình được kiểm soát bởi một dạng dao động duy nhất và hình dạng của dao động này giữ nguyên không đổi trong suốt thời gian phản ứng.

Thông thường, dạng dao động cơ bản được chọn là dạng phản ứng trội của hệ nhiều bậc tự do động, ảnh hưởng của các dạng dao động khác được xem là nhỏ và được bỏ qua. Phương pháp tính toán tĩnh phi tuyến với phân bố tải trọng ngang như vậy được gọi là phương pháp tính toán đẩy dần quy ước và thường được dùng để tính toán phản ứng của các công trình có chiều cao thấp và trung bình. Do tính đơn giản và khả năng xác định với độ chính xác chấp nhận được quá trình biến dạng của hệ kết cấu cũng như các bộ phận thành phần mà không cần phải thực hiện việc mô hình hóa phức tạp và tính toán công phu như tính toán động nên phương pháp tính toán đẩy dần được xem là một phương pháp hữu hiệu và tiện lợi trong tính toán động. Phương pháp phân tích dạng dao động và phổ phản ứng Phản ứng của kết cấu có nhiều bậc tự do chịu tác động động đất có thể được tính toán bằng cách phân tích hệ kết cấu thành nhiều hệ kết cấu có một bậc tự do tương đương.

Tính toán phản ứng mỗi hệ tương đương theo thời gian và sau đó cộng đại số các phản ứng lại để được phản ứng của kết cấu ban đầu. Phương pháp này được gọi là phương pháp phân tích dạng. Nếu việc tính toán chỉ nhằm xác định các đại lượng phản ứng lớn nhất thì tác động động đất sẽ được cho dưới dạng phổ phản ứng và kết quả tính toán theo phương pháp tích phân dạng dao động sẽ là phản ứng lớn nhất của hệ kết cấu. Phương pháp tính toán này có tên gọi là phương pháp phổ phản ứng.

Phương pháp tích phân dạng dao động cũng như phương pháp phổ phản ứng có những nhược điểm sau:  Phụ thuộc vào việc tách một cách nhân tạo các dạng dao động.  Phải tổ hợp các kết quả tính toán ở các dạng dao động lại theo nguyên tắc cộng tác dụng nên chỉ giới hạn ở giai đoạn làm việc đàn hồi tuyến tính của vật liệu. LUAN VAN CHAT LUONG download : add [e-mail] 7  Không áp dụng được cho một số hệ kết cấu không sử dụng được kỹ thuật phân tích dạng.  Không cho các chỉ dẫn chính xác về sự hình thành khớp dẻo ở một số cấu kiện.

Phương pháp tích phân trực tiếp phương trình chuyển động Phản ứng của các hệ kết cấu chịu tác động bất kỳ hoặc động đất có thể xác định được bằng cách tích phân trực tiếp phương trình chuyển động theo thời gian. Phương pháp này không cần thay đổi hoặc biến đổi các phương trình chuyển động sang hệ có một hoặc nhiều bậc tự do như ở phương pháp tích phân dạng dao động. Phương pháp tích phân trực tiếp theo thời gian xác định các giá trị gần đúng của nghiệm đối với một tập hợp các giá trị thời gian T được lựa chọn. Có thể tóm tắt nguyên tắc của phương pháp này như sau: (i) giả thiết các hàm mô tả sự biến thiên của chuyển vị, vận tốc và gia tốc trong một khoảng thời gian và (ii) các phương trình chuyển động không phải thỏa mãn ở tất cả ở mọi thời gian T mà chỉ trong khoảng thời gian không đổi Δt.

Khoảng thời gian này được gọi là bước thời gian. Điều này cũng có nghĩa rằng điều kiện cân bằng tĩnh của các lực quán tính, lực cản và lực đàn hồi với tải trọng tác động sẽ xảy ra ở nhiều bước thời gian Δt, 2Δt, …, nΔt, …. Ở mỗi bước thời gian, phương trình chuyển động được giải với các điều kiện ban đầu là chuyển vị, vận tốc được xác định ở bước trước đó. Phương pháp tích phân trực tiếp theo thời gian có thể áp dụng cho các hệ kết cấu tuyến tính lẫn phi tuyến nên có thể xem là phương pháp tổng quát duy nhất tính toán phản ứng động của các hệ kết cấu chịu tải trọng bất kỳ.

Kết luận chương 1 Chương 1 đã trình bày tóm lược lịch sử phát triển các phương pháp phân tích ứng xử của kết cấu chịu tải trọng động đất. Đồng thời cũng đưa ra một số ưu nhược điểm của các phương pháp đó. LUAN VAN CHAT LUONG download : add [e-mail] 8 Chương 2 XÂY DỰNG MÔ HÌNH PHẦN TỬ HỮU HẠN ĐỂ TÍNH TOÁN ĐÁP ỨNG ĐỘNG LỰC HỌC CỦA KHUNG, DẦM 2D-FGM CHỊU TẢI TRỌNG ĐỘNG ĐẤT 2. Dầm 2D-FGM Hình 2.1 minh họa dầm 2D-FGM với chiều dài L, chiều rộng b và chiều cao h trong hệ tọa độ Đề-các (x, z).

Hệ tọa độ (x, z) được chọn sao cho trục x trùng với mặt giữa của dầm, trục z vuông góc với mặt giữa và hướng lên trên. Mô hình dầm 2D-FGM Dầm được giả định được tạo thành từ bốn vật liệu thành phần, cụ thể là gốm 1, gốm 2, kim loại 1 và kim loại 2.