Chương I. Tổng quan Phân tích kết cấu là quá trình xác định ứng xử của hệ kết cấu khi chịu các dạng tải trọng. Phân tích đàn hồi tuyến tính giả thuyết bỏ qua ảnh hưởng của chuyển vị đến ứng xử của kết cấu và do đó quan hệ ứng suất – biến dạng là tuyến tính. Phân tích này thường đơn giản với khối lượng tính toán ít.
Dạng phân tích này hiện đang được áp dụng phổ biến để thiết kế kết cấu với việc kể đến tác động phi tuyến hình học và vật liệu một cách gián tiếp thông qua các công thức thiết kế hoặc các hệ số đơn giản nào đó được đề xuất trong các tiêu chuẩn. Tuy nhiên, dạng phân tích này chưa phản ánh đúng bản chất chịu lực thật của kết cấu. Ngược lại, trong bài toán phân tích phi tuyến quan hệ tải trọng – chuyển vị là phi tuyến, do đó cần phải sử dụng các thuật toán giải lặp để phân tích (vì kết cấu đã bị biến đổi về hình học và tính chất vật liệu cũng đã thay đổi). Do việc phân tích phải trải qua nhiều bước lặp và ma trận độ cứng luôn được cập nhật sau mỗi bước gia tải nên thời gian và khối lượng tính toán của bài toán phân tích phi tuyến sẽ lớn hơn nhiều so với bài toán phân tích đàn hồi tuyến tính.
Một phân tích phi tuyến cho khung thép cần kể đến các yếu tố chính sau: phi tuyến hình học và phi tuyến vật liệu. Phân tích phi tuyến hình học có kể đến ảnh hưởng do sự biến đổi hình học và sự phân bố ứng suất dư ban đầu trong cấu kiện, do đó ma trận độ cứng sẽ có thêm các ẩn số chuyển vị so với ma trận độ cứng thông thường. Nếu trong phân tích tuyến tính thì lời giải có thể tìm trực tiếp thì trong phân tích phi tuyến hình học lời giải phải dùng đến phương pháp gia tải từng bước do có sự biến đổi về mặt hình học của kết cấu sau mỗi bước tải. Phân tích phi tuyến vật liệu là phân tích có kể đến ứng xử phi đàn hồi của vật liệu.
Có hai phương pháp thường được sử dụng khi phân tích phi tuiyến vật liệu là phương pháp khớp dẻo và phương pháp vùng dẻo. 11 Luan van Phương pháp khớp dẻo (plastic hinge) còn gọi là phương pháp dầm-cột là mô hình đơn giản, dễ sử dụng và phổ biến nhất. Trong phương pháp khớp dẻo giả thiết sự chảy dẻo chỉ xảy ra trong một vùng nhỏ ở hai đầu phần tử, phần còn lại được giả thuyết vẫn còn đàn hồi. Phương pháp vùng dẻo (plastic zone), hay còn gọi là phương pháp dẻo phân bố (distributed plasticity), là phương pháp phần tử hữu hạn dựa trên việc chia cấu kiện thành nhiều phần tử dọc theo chiều dài và chia mặt cắt ngang tiết diện thành nhiều thớ.
Phương pháp này có thể: mô phỏng sự lan truyền dẻo qua mặt cắt ngang và dọc theo chiều dài cấu kiện. Tuy nhiên khối lượng tính toán và lưu trữ trong phương pháp này là khá lớn. Do vậy, phương pháp này thường chỉ được dùng trong nghiên cứu để kiểm tra độ tin cậy các phương pháp phân tích khác. Tình hình nghiên cứu Với tính hiệu quả về mặt tính toán, phương pháp dầm-cột đã được nghiên cứu sâu rộng trong phân tích khung thép chịu tải trọng tĩnh và động.
Phương pháp này dựa vào việc mô phỏng cấu kiện bằng việc chia cấu kiện thành một hay hai phần tử. Lui EM và Chen WF (1986) [14] phân tích ứng xử của khung thép phẳng dùng phương pháp khớp dẻo. Ứng xử phi tuyến của liên kết được mô phỏng bằng hàm mũ và có kể đến sự gia tải và dỡ tải của liên kết. Hsieh SH và Deierlein GG (1991) [9] phân tích phi tuyến khung không gian có liên kết nửa cứng.
Ứng xử phi tuyến vật liệu được kể đến bởi việc sử dụng phương pháp khớp dẻo có ma trận giảm dẻo dựa trên mặt dẻo ba tham số để mô phỏng sự chảy dẻo của mặt cắt ngang do tác động của lực dọc trục và mômen uốn theo hai phương. Chan SL và Zhou ZH (1994) [6] đã phát triển một phần tử mới dùng đa thức bậc 5 cho hàm chuyển vị của cấu kiện dầm-cột chịu tải đầu mút cho phân tích phi tuyến khung thép phẳng chịu tải trọng tĩnh. Liew JYR, Chen WF, Chen H (2000) [13] đã phát triển phương pháp khớp dẻo hiệu chỉnh dùng hai mặt chảy dẻo đồng dạng cho phép mô phỏng sự chảy dẻo dần 12 Luan van dần của đầu mút phần tử thay vì sự chảy dẻo đột ngột của thường thấy trong phân tích khớp dẻo đơn giản. Kim SE và Choi SH (2001) [11] trình bày một phương pháp phân tích nâng cao khung thép không gian có xét đến các yếu tố phi tuyến hình học, vật liệu và liên kết bằng cách dùng hàm ổn định và phương pháp khớp dẻo.
Vu (2006) [21] đề xuất một phần tử dầm-cột mới trong phân tích chuyển vị lớn khung thép phẳng đàn hồi sử dụng công thức tổng Lagrange dựa vào hàm nội suy chuyển vị. Ngo-Huu C, Kim SE và Oh JR (2008) [15] đề xuất phương pháp khớp dẻo thớ có chiều dài khớp thớ bằng không để phân tích phi tuyến vật liệu và dùng hàm ổn định truyền thống để phân tích phi tuyến hình học cho phần tử dầm-cột của khung thép không gian chịu tải tĩnh. Sau đó, Tai TH và Kim SE (2011) [24] cũng dùng phương pháp trên để phân tích ứng xử động của hệ khung thép không gian. Tuy nhiên, ở phương pháp trên, quan hệ lực dọc và chuyển vị bỏ qua ảnh hưởng của góc xoay hai đầu phần tử.
Ngo-Huu C, Kim SE (2009) [16] đã phát triển một phần tử dầm-cột khớp thớ phi tuyến cho mô phỏng khung thép không gian chịu tải tĩnh. Tác động phi đàn hồi được mô phỏng dựa vào phương pháp khớp dẻo hiệu chỉnh, cấu kiện được chia thành ba phần tử gồm hai phần tử khớp thớ hai đầu có chiều dài hữu hạn và một phần tử đàn hồi ở giữa. Hàm ổn định truyền thống có được từ lời giải giải tích của cấu kiện dầm-cột chịu lực dọc trục và mô-men uốn ở hai đầu được sử dụng để mô phỏng ứng xử bậc hai của phần tử đàn hồi ở giữa. Lực dọc bỏ qua ảnh hưởng của góc xoay hai đầu phần tử.
Chiorean CG (2009) [8] đã đề xuất một phương pháp dầm cột mới cho phân tích phi tuyến khung thép không gian có liên kết nửa cứng. Quan hệ lực – biến dạng phi đàn hồi phi tuyến và hàm ổn định được dùng để mô phỏng tác động phi tuyến vật liệu và hình học. 13 Luan van Chin-Long Lee và Filip C. Flippou (2009) [7] đề xuất một phần tử dầm-cột sử dụng phương pháp khớp thớ với chiều dài khớp thớ thay đổi (Spreading Inelastic Zone Element – SIZE) để phân tích ứng xử kết cấu dưới tác dụng của tải lặp.
Aslam Kassimali và Juan J. Garcilazo (2010) [1] phân tích phi tuyến hình học khung thép phẳng đàn hồi có xét đến ảnh hưởng của nhiệt độ. Phương pháp đề xuất sử dụng lý thuyết dầm-cột có xét đến hàm ổn định và hàm hiệu ứng cung được trình bày trước đó bởi Oran (1973) [18] và Kassimali (1976) [19]. Trong nghiên cứu này, quan hệ lực dọc và góc xoay hai đầu phần tử được kể đến thông qua các hàm hiệu ứng cung được thiết lập dựa trên các hàm ổn định truyền thống.
Lyon (2010) [22] đề xuất phần tử đồng xoay mới kết hợp lý thuyết khớp dẻo để áp dụng cho phân tích phi tuyến hình học và vật liệu cho khung thép. Phần tử đồng xoay được phát triển có ưu điểm là chỉ cần mô phỏng một phần tử cho một cấu kiện mà vẫn đạt độ chính xác cao, tuy nhiên, phần tử khớp dẻo được đề xuất vẫn là khớp dẻo cứng và sự chảy dẻo dần dần của khớp dẻo vẫn chưa được kể đến. Thanh-Nam Le, Jean-Marc Battini và Mohammed Hjiaj (2011) [25] đề xuất phần tử dầm đồng xoay trong phân tích động khung thép phẳng đàn hồi. Nghiên cứu sử dụng hàm chuyển vị xấp xỉ đa thức bậc 3 cho phần tử dầm nên quan hệ giữa mô- men và góc xoay chưa xét đến ảnh hưởng của lực dọc.
Nguyen Dinh Kien (2012) [17] đã trình bày công thức đồng xoay cho phần tử dầm-cột Timoshenko trong phân tích chuyển vị lớn của khung thép phẳng đàn hồi. Trong đó, ma trận độ cứng được xây dựng bằng phương pháp năng lượng, các chuyển vị được xác định dựa trên hàm dạng đa thức bậc 3 và thuật toán chiều dài cung được sử dụng trong bài toán lặp gia tăng. Bradford (2012) [10] đề xuất phần tử dầm-cột bậc 4 trong phân tích phi tuyến hình học cho khung thép đàn hồi. Ma trận độ cứng phần tử dầm- cột được xây dựng từ hàm năng lượng.
Tính cấp thiết của đề tài Thông thường, khi hệ kết cấu ứng xử phi tuyến, phương pháp phần tử hữu hạn (PTHH) được sử dụng để phân tích. Phương pháp này chia nhỏ một cấu kiện thành nhiều phần tử con, mức độ chính xác phụ thuộc vào số lượng phần tử con được chia. Do việc phân tích phải qua nhiều bước lặp và phải cập nhật lại ma trận độ cứng kết cấu sau mỗi bước gia tải nên khối lượng tính toán và dữ liệu lưu trữ của bài toán phân tích phi tuyến theo phương pháp này sẽ rất lớn. Việc giảm khối lượng tính toán mà vẫn đảm bảo độ chính xác ứng xử phi tuyến của hệ kết cấu là cần thiết và có tính thực tiễn cao.
Trong phương pháp phần tử đồng xoay, yếu tố phi tuyến hình học do sự tương tác giữa lực dọc và mô-men uốn được tính đến bằng việc cải tiến các hàm ổn định từ lời giải của phương trình vi phân cân bằng chịu tải đầu mút. Yếu tố phi tuyến vật liệu được được xét đến bằng cách sử dụng mô hình khớp dẻo hiệu chỉnh. Ưu điểm của việc sử dụng phần tử đồng xoay – khớp dẻo này là chỉ cần sử dụng một hoặc hai phần tử con trên một cấu kiện là có thể mô phỏng khá chính xác ứng xử phi tuyến của kết cấu, do đó hiệu quả tính toán sẽ cao hơn so với phương pháp PTHH truyền thống. Mục tiêu của đề tài • Phát triển phần tử dầm-cột đồng xoay kết hợp với các khớp dẻo hiệu chỉnh ở hai đầu phần tử để phân tích ứng xử phi tuyến hình học và phi tuyến vật liệu cho khung thép phẳng.
• Phát triển một chương trình phân tích tin cậy và hiệu quả cho phân tích ứng xử phi tuyến của khung thép phẳng. Đối tượng và phạm vi nghiên cứu Đối tượng nghiên cứu: • Khung thép phi đàn hồi chịu tải trọng tĩnh. 15 Luan van Phạm vi nghiên cứu: • Khung thép phẳng có liên kết cứng.