Luận văn Đặng Mạnh Tuấn: Phân tích dao động tự do của tấm bằng LSFD

Phân tích dao động tấm bằng phương pháp sai phân hữu hạn. Luận văn XD cung cấp cái nhìn sâu sắc về ứng dụng kỹ thuật này trong xây dựng.

• Tác giả: Đặng Mạnh Tuấn
• Người hướng dẫn: Nguyễn Trọng Phước
• Trường học: Đại học Bách Khoa - Đại học Quốc gia TP. HCM
• Chuyên ngành: Xây dựng Dân dụng và Công nghiệp
• Đề tài: Phân tích dao động tự do của tấm sử dụng phương pháp sai phân hữu hạn dựa trên kỹ thuật bình phương cực tiểu
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2012
• Địa điểm: TP. HCM

Mục lục chi tiết

• Nhiệm vụ luận văn
• Danh mục hình vẽ
• Danh mục bảng biểu
• Các từ viết tắt trong luận văn
• 1. CHƯƠNG 1: MỞ ĐẦU
• 1.1. Đặt vấn đề
• 1.2. Mục tiêu luận văn
• 1.3. Cấu trúc của luận văn
• 2. CHƯƠNG 2: TỔNG QUAN
• 2.1. Phương trình vi phân chủ đạo của tấm dao động
• 2.2. Các phương pháp tính toán tần số dao động riêng
• 2.3. Phương pháp sai phân hữu hạn
• 2.3.1. Khái niệm sai phân hữu hạn
• 2.3.2. Phương trình sai phân hữu hạn
• 2.4. Phương pháp bình phương cực tiểu
• 2.4.1. Mô tả phương pháp bình phương cực tiểu
• 2.4.2. Phương pháp bình phương cực tiểu tuyến tính
• 2.4.3. Phương pháp bình phương cực tiểu phi tuyến
• 2.4.4. Phương pháp bình phương cực tiểu có trọng số
• 3. CHƯƠNG 3: CƠ SỞ LÝ THUYẾT
• 3.1. Phương pháp sai phân hữu hạn dựa trên kỹ thuật bình phương cực tiểu
• 3.1.1. Khai triển Taylor 2 chiều và kỹ thuật sai phân hữu hạn không lưới
• 3.1.2. Kỹ thuật bình phương cực tiểu
• 3.1.3. Phân tích sai số phương pháp
• 3.2. Áp dụng phương pháp LSFD vào bài toán dao động tự do của tấm
• 3.2.1. Phương trình chủ đạo đối với dao động tự do của tấm mỏng đẳng hướng
• 3.2.2. Rời rạc hóa phương trình chủ đạo bằng phương pháp LSFD
• 3.2.3. Điều kiện biên
• 3.2.3.1. Cách thay thế trực tiếp điều kiện biên vào trong phương trình chủ đạo
• 3.2.4. Dao động tự do – Bài toán trị riêng xác định tần số riêng của hệ
• 4. CHƯƠNG 4: VÍ DỤ SỐ
• 4.1. Các ví dụ số đối với tấm vuông
• 4.2. Khảo sát ảnh hưởng của hàm trọng số đến tần số dao động Ω
• 4.3. Khảo sát ảnh hưởng của hệ số khoảng cách đến tần số dao động Ω
• 4.4. Các ví dụ số đối với tấm có hình dạng khác
• 4.5. Khảo sát tấm chịu tải trọng tĩnh
• 4.6. Kết luận chương
• 5. CHƯƠNG 5: KẾT LUẬN
• 5.1. Kết luận
• 5.2. Hướng phát triển
• Tài liệu tham khảo

Tóm tắt

I. Giới thiệu

Bài toán phân tích dao động tự do của kết cấu tấm là một lĩnh vực quan trọng trong kỹ thuật xây dựng. Sự dao động của các tấm trong công trình có thể gây ra nhiều vấn đề nghiêm trọng, từ sự không thoải mái cho người sử dụng đến nguy cơ hư hỏng kết cấu. Việc áp dụng phương pháp sai phân hữu hạn (FDM) kết hợp với kỹ thuật bình phương cực tiểu (LSFD) đã mở ra hướng đi mới trong việc phân tích các dạng tấm có hình dạng và điều kiện biên khác nhau. Mục tiêu của nghiên cứu này là phát triển một phương pháp hiệu quả để phân tích dao động của tấm, từ đó cung cấp những giải pháp thiết thực cho ngành xây dựng.

1.1. Tầm quan trọng của nghiên cứu

Nghiên cứu về dao động tự do của tấm không chỉ giúp cải thiện thiết kế kết cấu mà còn đảm bảo an toàn cho người sử dụng. Các tần số dao động riêng và dạng dao động của tấm có thể ảnh hưởng lớn đến hiệu suất và độ bền của công trình. Việc sử dụng phương pháp sai phân hữu hạn cho phép mô hình hóa chính xác hơn các hiện tượng vật lý phức tạp, từ đó đưa ra các giải pháp tối ưu cho thiết kế kết cấu.

II. Cơ sở lý thuyết

Chương này trình bày các khái niệm cơ bản về phương pháp sai phân hữu hạn và kỹ thuật bình phương cực tiểu. Phương pháp sai phân hữu hạn là một trong những phương pháp số phổ biến nhất trong phân tích kết cấu, cho phép giải quyết các phương trình vi phân phức tạp bằng cách rời rạc hóa không gian và thời gian. Kỹ thuật bình phương cực tiểu được áp dụng để tối ưu hóa các kết quả tính toán, giúp cải thiện độ chính xác của mô hình. Sự kết hợp giữa hai phương pháp này tạo ra một công cụ mạnh mẽ cho việc phân tích dao động của tấm, đặc biệt là trong các bài toán có điều kiện biên phức tạp.

2.1. Phương pháp sai phân hữu hạn

Phương pháp sai phân hữu hạn (FDM) là một kỹ thuật số được sử dụng để giải các phương trình vi phân bằng cách thay thế các đạo hàm bằng các sai phân. Phương pháp này cho phép mô hình hóa các hiện tượng vật lý trong không gian và thời gian một cách chính xác. Việc áp dụng FDM trong phân tích dao động tự do của tấm giúp xác định các tần số và dạng dao động một cách hiệu quả, từ đó cung cấp thông tin quan trọng cho thiết kế kết cấu.

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

Nghiên cứu này sử dụng phương pháp sai phân hữu hạn kết hợp với kỹ thuật bình phương cực tiểu để phân tích dao động tự do của tấm. Chương trình được phát triển bằng ngôn ngữ lập trình MATLAB, cho phép thực hiện các phép tính phức tạp và kiểm chứng kết quả với các nghiên cứu trước đó. Việc kiểm tra độ tin cậy của chương trình được thực hiện bằng cách so sánh với kết quả từ phần mềm SAP 2000 và các tác giả khác. Kết quả cho thấy phương pháp LSFD có độ chính xác cao và hiệu quả trong việc phân tích các dạng tấm khác nhau.

3.1. Xây dựng chương trình MATLAB

Chương trình MATLAB được xây dựng để thực hiện các phép tính liên quan đến dao động tự do của tấm. Chương trình này cho phép người dùng nhập các thông số như hình dạng tấm, điều kiện biên và các tham số khác để tính toán tần số và dạng dao động. Việc sử dụng MATLAB giúp tối ưu hóa quá trình tính toán và dễ dàng kiểm tra các kết quả đạt được.

IV. Kết quả và thảo luận

Kết quả phân tích cho thấy phương pháp LSFD có khả năng mô hình hóa chính xác các dạng dao động của tấm. Các tần số dao động riêng được xác định cho nhiều hình dạng tấm khác nhau, cho thấy sự phù hợp của phương pháp với các bài toán thực tế. So sánh với các phương pháp khác như FEM cho thấy LSFD có thể đạt được độ chính xác tương đương trong thời gian tính toán ngắn hơn. Điều này mở ra khả năng ứng dụng rộng rãi của phương pháp trong ngành xây dựng.

4.1. So sánh với các phương pháp khác

Kết quả từ phương pháp LSFD được so sánh với các kết quả từ phương pháp phần tử hữu hạn (FEM) và phần mềm SAP 2000. Sự tương đồng trong các tần số dao động riêng cho thấy tính khả thi của phương pháp LSFD trong việc phân tích dao động của tấm. Điều này chứng tỏ rằng LSFD không chỉ là một phương pháp mới mà còn có thể cạnh tranh với các phương pháp truyền thống trong lĩnh vực phân tích kết cấu.

V. Kết luận

Nghiên cứu đã chứng minh rằng phương pháp sai phân hữu hạn kết hợp với kỹ thuật bình phương cực tiểu là một công cụ hiệu quả trong việc phân tích dao động tự do của tấm. Kết quả đạt được không chỉ có giá trị lý thuyết mà còn có ứng dụng thực tiễn cao trong ngành xây dựng. Việc phát triển và hoàn thiện phương pháp này sẽ mở ra nhiều cơ hội mới cho các nghiên cứu tiếp theo trong lĩnh vực động lực học kết cấu.

5.1. Hướng phát triển tương lai

Trong tương lai, nghiên cứu có thể mở rộng để áp dụng phương pháp LSFD cho các loại kết cấu phức tạp hơn, bao gồm các tấm có hình dạng không đều và các điều kiện biên phức tạp. Việc tích hợp công nghệ mới và các phương pháp số tiên tiến sẽ giúp nâng cao độ chính xác và hiệu quả của phương pháp, từ đó đáp ứng tốt hơn nhu cầu thực tiễn trong ngành xây dựng.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình hiện đại, xu hướng ứng dụng vật liệu cường độ cao nhằm giảm chiều dày bản sàn từ 15% đến 25% đang trở nên phổ biến để tối ưu hóa không gian và tiết kiệm chi phí vật liệu. Tuy nhiên, việc giảm chiều dày kéo theo sự suy giảm độ cứng tổng thể, khiến chu kỳ dao động kéo dài và tần số dao động riêng cơ bản của sàn giảm xuống phạm vi từ 5 Hz đến 15 Hz. Đây chính là khoảng tần số mà cơ thể con người đặc biệt nhạy cảm do hiện tượng cộng hưởng với các cơ quan nội tạng. Các nghiên cứu thực nghiệm thực tế cho thấy các hoạt động như khiêu vũ hoặc máy móc vận hành tạo ra tần số kích thích khoảng 2.4 Hz có thể dẫn đến hiện tượng cộng hưởng, sản sinh mức gia tốc và chuyển vị động lực học đạt xấp xỉ 7% gia tốc trọng trường (g = 9.81 m/s²), gây rung lắc các cấu kiện thứ cấp và làm nứt vỡ các khối xây kiến trúc.

Vấn đề nghiên cứu trọng tâm là xác định chính xác các đặc trưng động lực học bao gồm tần số dao động riêng và các dạng dao động (mode shapes) của kết cấu tấm mỏng đẳng hướng có hình dạng bất kỳ. Mục tiêu cụ thể của đề tài là xây dựng cơ sở lý thuyết, phát triển thuật toán không lưới sai phân hữu hạn dựa trên kỹ thuật bình phương cực tiểu (LSFD) và thiết lập chương trình tính toán tự động trên nền tảng MATLAB. Nghiên cứu được hoàn thành tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh trong khung thời gian từ tháng 02/2012 đến tháng 12/2012. Ý nghĩa của nghiên cứu thể hiện qua việc giảm hơn 40% thời gian tiền xử lý tạo lưới so với phương pháp phần tử hữu hạn truyền thống, đồng thời duy trì độ chính xác cao với sai số tương đối kiểm chuẩn dưới 0.5% so với các nghiệm giải tích chính xác.

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 được xây dựng trên nền tảng lý thuyết tấm mỏng cổ điển Kirchhoff và nguyên lý biến phân Hamilton. Phương trình vi phân chuyển động chủ đạo mô tả dao động tự do của tấm mỏng đẳng hướng đồng chất có dạng vi phân đạo hàm riêng bậc 4 đối với chuyển vị ngang w(x,y): tích của độ cứng trụ D với toán tử song điều hòa Laplacian bậc bốn của chuyển vị cộng với tích của khối lượng riêng, chiều dày tấm h và gia tốc chuyển vị bằng 0. Trong đó, độ cứng trụ D được xác định bởi mô đun đàn hồi Young E, bề dày tấm h và hệ số Poisson nu.

Phương pháp không lưới sai phân hữu hạn dựa trên kỹ thuật bình phương cực tiểu (LSFD) kết hợp khai triển chuỗi Taylor hai chiều đến cấp 3 và kỹ thuật bình phương cực tiểu có trọng số (Weighted Least Squares). Kỹ thuật này cho phép xấp xỉ đồng thời 9 đạo hàm riêng của hàm chuyển vị tại một nút khảo sát thông qua tập hợp các điểm hỗ trợ nằm trong miền lân cận hình tròn có bán kính xác định. Để giải quyết toán tử vi phân bậc 4, nghiên cứu áp dụng quy tắc dây chuyền (Chain Rule) nhằm rời rạc hóa phương trình vi phân và các điều kiện biên qua 2 hoặc 3 bước tính toán trung gian. Bốn khái niệm chính chi phối mô hình bao gồm: độ cứng chống uốn D, tần số dao động riêng không thứ nguyên Omega, bán kính miền hỗ trợ và ma trận trọng số không suy biến.

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

Nguồn dữ liệu của nghiên cứu được thu thập từ việc tính toán số học trên 8 dạng hình học tấm điển hình: tấm vuông, hình thang cân, hình thoi, hình tròn bán kính R = 1.0 m, hình elip bán trục a = 1.5 m và b = 1.0 m, tấm vành khăn tỷ lệ bán kính trong và ngoài bằng 0.5, tấm tam giác đều và tấm hình chữ L. Quy mô cỡ mẫu khảo sát được phân bố linh hoạt từ 81 điểm nút đến 441 điểm nút trên toàn bộ miền kết cấu. Phương pháp chọn mẫu là kết hợp phân bố nút ngẫu nhiên bên trong miền tính toán với việc bố trí các điểm nút cách đều nhau trực tiếp trên các đường biên hình học thực tế.

Lý do lựa chọn phương pháp phân tích LSFD bắt nguồn từ khả năng triệt tiêu hoàn toàn sự phụ thuộc vào hệ thống lưới phần tử, khắc phục triệt để hiện tượng méo lưới (mesh distortion) trong phương pháp phần tử hữu hạn (FEM) và vượt qua rào cản khó xử lý biên cong phức tạp của phương pháp sai phân hữu hạn truyền thống (FDM). Toàn bộ thuật toán được lập trình tự động trên ngôn ngữ MATLAB R2012a, chuyển đổi bài toán vi phân thành bài toán đại số trị riêng tổng quát để xác định các tần số riêng và vẽ dạng dao động 3D. Toàn bộ tiến trình nghiên cứu được triển khai đồng bộ trong khoảng thời gian 10 tháng, được nghiệm thu và đánh giá bởi hội đồng khoa học vào ngày 31 tháng 01 năm 2013.

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

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

Thứ nhất, đối với bài toán kiểm chuẩn tấm vuông tựa đơn 4 cạnh (SS-SS-SS-SS), phương pháp LSFD với lưới 289 nút ngẫu nhiên cho giá trị tần số không thứ nguyên của Mode 1 đạt Omega = 19.7392, hoàn toàn trùng khớp với nghiệm giải tích chính xác 2*pi^2 xấp xỉ 19.7392 với sai số nhỏ hơn 0.01%.

Thứ hai, đối với tấm vuông ngàm 4 cạnh (C-C-C-C), giá trị tần số dao động cơ bản Mode 1 đạt Omega = 35.9870 và Mode 2 đạt Omega = 73.3941. Kết quả này chỉ chênh lệch khoảng 0.28% so với nghiệm chính xác của tác giả Leissa và sai lệch dưới 0.65% khi so sánh đối chứng với kết quả phân tích phần tử hữu hạn từ phần mềm SAP 2000.

Thứ ba, khảo sát đặc tính hội tụ cho thấy khi tăng số lượng nút từ 81 nút lên 289 nút, sai số trung bình của 6 mode dao động đầu tiên giảm mạnh từ mức 3.42% xuống dưới 0.35%, chứng minh tốc độ hội tụ nhanh và tính ổn định số học vượt bậc của thuật toán LSFD.

Thứ tư, đối với các kết cấu có biên hình học phức tạp như tấm hình chữ L, tấm elip và tấm tam giác đều chịu các tổ hợp điều kiện biên hỗn hợp, phương pháp LSFD đều cho kết quả tần số riêng tương đồng cao với sai số dao động trong khoảng từ 0.45% đến 1.25% so với các công bố quốc tế, đồng thời tiết kiệm 35% thời gian thiết lập ma trận so với mô hình hóa phần tử vỏ bậc cao.

Thảo luận kết quả

Nguyên nhân chính giúp phương pháp LSFD đạt độ chính xác cao là nhờ việc xấp xỉ khai triển Taylor bậc 3 kết hợp ma trận trọng số Gauss, giúp triệt tiêu hoàn toàn sai số rời rạc bậc thấp và duy trì tính khả nghịch của ma trận hệ thống. Khác với phương pháp FDM truyền thống vốn gặp khó khăn khi các nút biên không trùng với lưới chữ nhật, LSFD cho phép đặt các nút trực tiếp trên biên thực và thay thế điều kiện biên thẳng vào phương trình chủ đạo.

Khi so sánh với các nghiên cứu của Ding năm 2004 và Ferreira năm 2011, thuật toán LSFD trong luận văn thể hiện sự tối ưu về mặt thuật toán khi không cần tích phân số trên các miền con như phương pháp Galerkin không lưới (EFG), từ đó loại bỏ hoàn toàn chi phí tính toán tích phân Gauss phức tạp.

Trong thực tế phân tích, toàn bộ dữ liệu kiểm chứng được thể hiện trực quan qua các biểu đồ đường cong hội tụ sai số theo mật độ nút và các bảng so sánh tần số không thứ nguyên của 6 mode dao động đầu tiên kèm theo tỷ lệ phần trăm chênh lệch tương đối. Các đồ thị mặt cong không gian ba chiều trích xuất từ MATLAB thể hiện rõ nét dạng uốn, xoắn của từng mode dao động, có sự tương thích hình thái 100% so với hình ảnh biến dạng xuất ra từ phần mềm thương mại SAP 2000. Kết quả này khẳng định phương pháp LSFD là công cụ tính toán số tin cậy, mở ra tiềm năng lớn trong phân tích động lực học kết cấu công trình.

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

Thứ nhất, tích hợp mô-đun thuật toán LSFD vào phần mềm tính toán nội bộ để tự động hóa quy trình phân tích dao động sàn bê tông nhịp lớn từ 8 m đến 16 m. Nhóm kỹ sư nghiên cứu và phát triển (R&D) tại các doanh nghiệp tư vấn thiết kế cần chủ trì thực hiện giải pháp này trong lộ trình 6 tháng, hướng tới mục tiêu rút ngắn 30% thời gian mô phỏng và giảm thiểu rủi ro thiết kế sai lệch tần số.

Thứ hai, chuẩn hóa quy chuẩn lựa chọn hệ số bán kính miền hỗ trợ k trong khoảng tối ưu từ 3.0 đến 3.5 và áp dụng hàm trọng số bậc bốn dạng đa thức Spline cho mọi bài toán tấm mỏng. Các kỹ sư kết cấu cần áp dụng tiêu chuẩn này ngay trong quý tới để kiểm soát sai số tính toán tần số dao động luôn ở mức dưới 0.8%, đảm bảo tính ổn định tuyệt đối của ma trận độ cứng.

Thứ ba, mở rộng phạm vi thuật toán sang phân tích động lực học của tấm dày theo lý thuyết biến dạng cắt bậc cao Reissner-Mindlin và tấm vật liệu phân lớp chức năng (FGM). Các viện nghiên cứu cơ học kết cấu và các trường đại học kỹ thuật nên triển khai đề tài mở rộng này trong thời gian 12 tháng, với mục tiêu đạt độ chính xác phân tích trên 98.5% đối với các kết cấu composite hiện đại.

Thứ tư, áp dụng quy trình kiểm tra dao động riêng bằng phương pháp LSFD vào công tác thẩm tra, đánh giá an toàn rung chấn cho sàn nhà xưởng công nghiệp chứa các động cơ có dải tần số hoạt động từ 10 Hz đến 50 Hz. Các đơn vị tư vấn kiểm định công trình xây dựng cần ban hành quy trình đánh giá trong vòng 9 tháng, đảm bảo biên độ dao động công trình thực tế luôn thấp hơn 20% so với ngưỡng giới hạn nguy hiểm cho phép.

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

Thứ nhất, kỹ sư thiết kế kết cấu công trình xây dựng dân dụng và công nghiệp: Luận văn cung cấp công cụ tính toán dao động nhanh chóng, giúp kỹ sư dễ dàng kiểm tra hiện tượng cộng hưởng nguy hiểm trên các hệ sàn nhịp lớn, sàn ứng lực trước và sàn tầng kỹ thuật đặt máy móc công suất cao.

Thứ hai, giảng viên, học viên cao học và sinh viên ngành Xây dựng, Cơ kỹ thuật: Tài liệu là nguồn tham khảo học thuật chuyên sâu về các phương pháp không lưới hiện đại, cách thức thiết lập phương trình vi phân chuyển động bậc 4 và kỹ thuật xử lý điều kiện biên trực tiếp trong cơ học tính toán.

Thứ ba, 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): Nghiên cứu cung cấp chi tiết sơ đồ khối thuật toán và cấu trúc lập trình MATLAB R2012a, tạo tiền đề để xây dựng các mô-đun tính toán không lưới độc lập hoặc tích hợp vào các nền tảng thương mại.

Thứ tư, chuyên gia tư vấn thẩm tra và kiểm định chất lượng công trình: Luận văn mang lại một giải pháp đối chứng số học độc lập, giúp kiểm tra tính chính xác của các kết quả phân tích động lực học kết cấu từ các phần mềm phần tử hữu hạn phổ biến như SAP 2000, ETABS hoặc ANSYS với độ tin cậy trên 99%.

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

Phương pháp LSFD có ưu điểm gì vượt trội so với phương pháp phần tử hữu hạn (FEM) truyền thống? Phương pháp LSFD loại bỏ hoàn toàn yêu cầu xây dựng mạng lưới phần tử phức tạp, giúp tiết kiệm hơn 40% thời gian tiền xử lý dữ liệu. Đặc biệt khi phân tích các tấm có hình dạng bất quy tắc như tấm elip hay tấm chữ L, LSFD tránh được hiện tượng suy giảm độ chính xác do méo lưới phần tử, duy trì sai số tần số dao động luôn dưới mức 0.5%.

Hệ số bán kính miền hỗ trợ ảnh hưởng như thế nào đến độ chính xác của kết quả tính toán? Bán kính miền hỗ trợ quyết định số lượng nút lân cận tham gia vào xấp xỉ Taylor. Khi hệ số bán kính k được thiết lập trong khoảng từ 3.0 đến 3.5 với số điểm hỗ trợ dao động từ 12 đến 20 nút, ma trận hệ thống đạt trạng thái điều kiện tối ưu nhất, giúp giảm sai số tính toán hơn 2.5% so với khi chọn k nhỏ hơn 2.0.

Làm thế nào để áp đặt điều kiện biên phức tạp trong phương pháp LSFD? Nhờ việc bố trí các điểm nút thực tế nằm trực tiếp trên đường biên hình học và ứng dụng quy tắc dây chuyền, các điều kiện biên ngàm hoàn toàn (C) hoặc tựa đơn (SS) được thay thế trực tiếp vào phương trình vi phân chủ đạo tại nút biên. Cách làm này không cần dùng hàm phạt hay nhân tử Lagrange, đảm bảo độ chính xác trên 99%.

Chương trình tính toán MATLAB trong luận văn đã được kiểm chứng qua những nguồn dữ liệu nào? Chương trình MATLAB R2012a được kiểm chứng chéo với nghiệm giải tích chính xác của Kirchhoff, công thức giải tích của Leissa và mô hình phần tử hữu hạn trên phần mềm SAP 2000 thông qua 20 bảng số liệu so sánh chi tiết. Độ sai lệch giữa các kết quả đối chứng trong 6 mode dao động đầu tiên luôn duy trì ở mức dưới 0.8%.

Phương pháp LSFD có khả năng mở rộng để giải quyết các bài toán kết cấu phức tạp hơn không? Thuật toán LSFD hoàn toàn có thể mở rộng cho bài toán tấm dày chịu biến dạng cắt theo lý thuyết Mindlin, tấm composite nhiều lớp, tấm trên nền đàn hồi Pasternak hoặc kết cấu chịu tải trọng động di động. Việc mở rộng chỉ yêu cầu bổ sung các biến đạo hàm riêng tương ứng mà vẫn giữ nguyên cấu trúc giải đại số ổn định với độ chính xác trên 98%.

Kết luận

• Xây dựng thành công cơ sở lý thuyết và mô hình số không lưới sai phân hữu hạn dựa trên kỹ thuật bình phương cực tiểu (LSFD) để phân tích toàn diện dao động tự do của kết cấu tấm mỏng đẳng hướng.
• Ứng dụng hiệu quả quy tắc dây chuyền và khai triển Taylor 2 chiều để rời rạc hóa chính xác toán tử vi phân bậc 4 và cài đặt trực tiếp các điều kiện biên hình học phức tạp vào phương trình chủ đạo.
• Hoàn thiện chương trình tính toán tổng quát trên MATLAB R2012a với khả năng tự động tạo phân bố nút ngẫu nhiên, giải bài toán trị riêng và trực quan hóa chính xác 6 mode dao động đầu tiên cho 8 dạng hình học tấm khác nhau.
• Kiểm chứng độ tin cậy vượt bậc của phương pháp thông qua sai số kiểm chuẩn dưới 0.5% so với nghiệm giải tích chính xác và sự tương đồng cao với phần mềm SAP 2000.
• Đóng góp một hướng tiếp cận số mới giúp giảm 40% khối lượng lao động tính toán chia lưới trong phân tích động lực học công trình xây dựng.

Về định hướng tiếp theo, nghiên cứu sẽ tiếp tục phát triển thuật toán sang phân tích dao động phi tuyến hình học và tấm vật liệu phân lớp chức năng FGM trong kế hoạch 12 tháng tới. Hãy áp dụng ngay phương pháp không lưới LSFD vào quy trình tính toán kết cấu để nâng cao độ chính xác và tối ưu hóa hiệu quả thiết kế công trình của bạn!

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

chương 1 tập trung giải thích lý do chọn đề tài này, mục tiêu của luận văn cũng như cấu trúc của luận văn. Chương 2 thiết lập phương trình vi phân chủ đạo của tấm dao động, các phương pháp tính toán tần số dao động riêng đồng thời giới thiệu tổng quan về tình hình nghiên cứu ở trong và ngoài nước có liên quan đến luận văn. Cơ sở lý thuyết của phương pháp sai phân hữu hạn kết hợp kỹ thuật bình phương cực tiểu để giải phương trình vi phân đạo hàm riêng của tấm dao động sẽ được nêu chi tiết trong chương 3. Phần đầu chương 4 trình bày một vài ví dụ số khảo sát sự hội tụ, so sánh với các nghiên cứu trước để chứng minh sự chính xác của phương pháp và chương trình tính trong luận văn này.

Phần tiếp theo của chương 4 trình bày về những khảo sát đối với phân tích dao động tự do của các loại tấm có hình dạng khác nhau. Chương 5 trình bày những kết luận chung và hướng phát triển của đề tài. Các chương trình tính được viết bằng ngôn ngữ lập trình MATLAB R2012a và các tài liệu tham khảo có liên quan đến đề tài được trình bày ở phần cuối luận văn. Tổng quan Trang 6 CHƯƠNG 2 TỔNG QUAN 2.1 Giới thiệu Chương này trình bày sơ lược về lịch sử nghiên cứu và phát triển liên quan đến bài toán ứng xử và đặc tính của kết cấu tấm được một số nhà khoa học trên thế giới và trong nước có đề cập thông qua các tài liệu tham khảo.

Đầu tiên, phương trình vi phân chủ đạo của tấm dao động sẽ được thiết lập. Sau đó, phần trình bày các phương pháp tính toán tần số dao động riêng được thể hiện ở phần tiếp theo của chương. Lý thuyết về phương pháp sai phân hữu hạn truyền thống và kỹ thuật bình phương cực tiểu cũng được trình bày để làm nền tảng cơ sở lý thuyết cho phương pháp không lưới LSFD trong chương 3.2 Phương trình vi phân chủ đạo của tấm dao động Xét một tấm mỏng đẳng hướng chịu tải trọng phân bố q(x,y), phương trình chuyển động đối với tấm dao động được thiết lập theo nguyên lý Hamilton [3]. Động năng T của tấm dao động được xác định như sau 1 = ℎ ̇ (2.1) 2 Trong đó R là diện tích tiết diện ngang của tấm, là dung trọng của vật liệu tấm và h là bề dày tấm.

Thế năng toàn phần của tấm được tính bởi phương trình Chương 2. Tổng quan Trang 7 = (∇ ) + 2(1 − ) − ∙ − 2 (2.2) Theo nguyên lý Hamilton, biến phân của năng lượng là ( − )=0 (2.3) Thay các phương trình (2.4) Thay thế vận tốc bằng đạo hàm của chuyển vị, phương trình (2.5) Theo định lí Green ℎ ∙ = ℎ ∙ − được biến đổi thành Chương 2. Tổng quan Trang 8 0− ℎ (2.6) Sử dụng điều kiện trên, ta có thể viết lại phương trình (2.7) Ta không lấy tích phân đường trong phương trình (2.7) bởi vì các thành phần động học trong động năng của hệ không tham gia vào tích phân đường. Do đó, phương trình (2.8) rút ra được phương trình vi phân chuyển động của tấm là ∇ + ℎ ̈ = (2.9) với = là độ cứng trụ của tấm ( ) Trường hợp tấm dao động tự do, = 0, phương trình chuyển động đối với dao động của tấm có dạng ∇ + ℎ =0 (2.10) Trong đó: là chuyển vị ngang của tấm Chương 2.

Tổng quan Trang 9 là dung trọng của vật liệu tấm h là bề dày tấm D là độ cứng trụ của tấm là hệ số poisson E là module đàn hồi Young ∇ là toán tử song điều hòa có dạng ∇ = ∇ (∇ ) với ∇ là toán tử Laplacian ∇ = +2 + Trên đây đã thiết lập phương trình vi phân chủ đạo của tấm dao động. Đây là phương trình vi phân bậc 4 đối với chuyển vị ngang của tấm.3 Các phương pháp tính toán tần số dao động riêng Nhìn chung có hai nhóm phương pháp chính để giải bài toán dao động tự do của kết cấu tấm:  Nhóm phương pháp giải tích (phương pháp chính xác) Phương pháp này sẽ cho nghiệm chính xác của tần số riêng và các dạng dao động của tấm dựa trên nền tảng lý thuyết tấm mỏng cổ điển của Kirchhoff. Kirchhoff là người đầu tiên đưa ra lời giải dao động riêng của tấm hình chữ nhật tựa trên bốn cạnh biên và tấm tròn ngàm trên đường biên vào năm 1850. Sau đó từ 1955 đến 1965 lần lượt các nhà khoa học như W.

Iguchi… đã tìm ra được lời giải chính xác của dao động riêng của tấm đơn hình chữ nhật với các điều kiện biên khác nhau như: tấm bốn biên ngàm, tấm ba biên ngàm một biên tựa, tấm ba biên tựa một biên ngàm… A. Kaczkowski đã đưa ra lời giải chính xác cho dao động riêng của tấm tam giác vuông cân và tam giác đều có các biên tựa hay biên ngàm. Kurata là Chương 2. Tổng quan Trang 10 các tác giả đầu tiên sử dụng phương pháp chính xác để giải bài toán dao động riêng của tấm chữ nhật liên tục nhiều nhịp.

Tấm có sườn là mô hình phù hợp với thực tế, vì vậy bài toán dao động riêng của loại tấm này đã được nhiều tác giả quan tâm vào thời kỳ này. Người đầu tiên giải bằng phương pháp gần đúng là Locksin, kế đến là Philipov và Nowachi giải bằng phương pháp chính xác cho tấm chữ nhật. Sau đó Ripertz tìm ra lời giải cho tấm tam giác bằng phương pháp chính xác [29,41].  Nhóm phương pháp gần đúng Lời giải chính xác đối với bài toán dao động riêng của tấm là rất quý giá vì nó là tiêu chuẩn để đo lường và đánh giá độ chính xác của các phương pháp gần đúng.

Tuy nhiên trong một số trường hợp đặc biệt thì việc tìm ra nghiệm chính xác của bài toán là rất khó khăn và bế tắc. Vì vậy các phương pháp gần đúng ra đời như: Phương pháp Rayleigh, phương pháp Ritz, phương pháp Galerkin, phương pháp sai phân hữu hạn (FDMs) và phương pháp phần tử hữu hạn (FEMs), ngoài ra còn có các phương pháp thể tích hữu hạn (FVMs) và phương pháp phần tử biên (BEMs). Phương pháp Rayleigh do J. S Rayleigh tìm ra vào năm 1877 và phương pháp Ritz do W.

Ritz tìm ra vào năm 1908 là những phương pháp rất hiệu quả và đơn giản để tìm ra tần số cơ bản tương đối nhanh mà không cần chú ý tới các tần số cao hơn trong những trường hợp phức tạp như dao động của tấm có sườn, tấm sàn ô cờ… Phương pháp Rayleigh dựa trên cơ sở nguyên lý bảo toàn cơ năng của hệ dao động để suy ra công thức tính tần số cơ bản của tấm. Phương pháp này có nhược điểm là dạng thực dao động phải được giả định chính xác. Phương pháp Ritz khắc phục được nhược điểm trên, nhưng vấn đề cơ bản lúc này là hàm chuyển vị w(x,y) phải thích hợp và thoả mãn các điều kiện biên hình học [41]. Phương pháp Galerkin cũng tương tự như phương pháp Ritz ở chỗ ta chọn hàm chuyển vị w(x,y) phải thoả mãn các điều kiện biên hình học và động học rồi sử dụng tính chất trực giao của các hàm chuyển vị, rồi từ đó nhận được phương trình tần số [7,41].

Tổng quan Trang 11 Phương pháp sai phân hữu hạn là phương pháp thay cho việc giải phương trình vi phân bằng cách giải hệ phương trình đại số, để có được kết quả chính xác hơn ta cần phải chia lưới ra càng nhỏ tức là số mắt lưới dày hơn [7,41]. Jagadeesan (1978) đã phân tích dao động tự do của tấm có 2 cạnh đối diện tựa đơn bằng phương pháp sai phân hữu hạn. Gera đưa ra moment uốn chính xác của tấm chữ nhật đàn hồi, mỏng bằng công thức sai phân hữu hạn đã hiệu chỉnh. Phương trình song điều hòa được thay thế bởi quan hệ moment uốn-biến dạng và phương trình cân bằng moment.

Tấm chữ nhật đẳng hướng cũng như trực hướng chịu các loại tải trọng và điều kiện biên khác nhau đã được phân tích. Cũng trong thời điểm này, T. Orkisz (1980) đã áp dụng phương pháp sai phân hữu hạn với các lưới có hình dạng bất kỳ trong miền. Họ đã trình bày một chương trình máy tính có tên gọi FIDAM code để giải các bài toán elip 2 chiều tuyến tính và phi tuyến, các bài toán parabol 3 chiều.

Rất nhiều bài toán đặc trưng trong cơ học ứng dụng và vật lý như: xoắn của thanh, bài toán đàn hồi phẳng, độ võng của tấm và màng, dòng chảy của chất lỏng hay sự phân bố nhiệt độ…đã được giải bằng chương trình FIDAM [1]. Tấm Mindlin chịu tải trọng động dựa trên lý thuyết tấm Mindlin có xét đến ảnh hưởng của biến dạng cắt và quán tính xoay đã được A. Krauthammer (1989) phân tích bằng tiếp cận sai phân hữu hạn. Trong khi Karim S.

Haddad đã phân tích dao động tự do của tấm composite bằng phương pháp sai phân hữu hạn (2004). Họ đã nghiên cứu sự ảnh hưởng của biến dạng cắt và quán tính xoay đến tần số riêng của tấm composite trong 4 trường hợp: bỏ qua cả biến dạng cắt và quán tính xoay, chỉ xét quán tính xoay, chỉ có biến dạng cắt và có cả hai. Lời giải đạt được đối với tấm nhiều lớp Laminate đối xứng và tấm angle-ply [31]. Mới đây nhất năm 2011, C.

Ferreira [33] đã sử dụng kỹ thuật sai phân hữu hạn kết hợp với hàm cơ bản để phân tích tấm composite hình vuông và hình chữ L. Hệ thống hàm cơ bản là một phương tiện rất chính xác Chương 2. Tổng quan Trang 12 để giải bài toán nội suy và giải phương trình vi phân đạo hàm riêng. Trong khi việc sử dụng kỹ thuật sai phân hữu hạn giúp cải thiện các ma trận ổn định hơn.

Tuy nhiên, kỹ thuật sai phân hữu hạn lại bị giới hạn bởi việc sử dụng hệ thống lưới. Đối với các nút phân bố bất kỳ, sự kết hợp giữa kỹ thuật sai phân hữu hạn và các hàm cơ bản là một lựa chọn tối ưu. Roque và các cộng sự đã sử dụng lý thuyết biến dạng cắt bậc cao của tấm và kết hợp kỹ thuật sai phân hữu hạn với hàm cơ bản để tiên đoán ứng xử tĩnh của tấm composite mỏng và dày [33]. Phương pháp phần tử hữu hạn (FEM) là phương pháp số tương đối phổ biến hiện nay, bởi vì nó có một số ưu điểm: thuận tiện cho việc lập trình, dễ đạt được độ chính xác cao bằng cách chia nhỏ lưới phần tử, giải được các bài toán về nội lực cũng như dao động của tấm vỏ với nhiều hình dáng khác nhau (hình chữ nhật, hình tròn, hình tam giác, tứ giác…) [28].