Bài giảng Phương pháp phần tử hữu hạn - Bài 4 (Hệ thanh - Phần 1/3: Các hàm dạng và các ma trận phần tử) - Nguyễn Xuân Thành

Nội dung text: Bài giảng Phương pháp phần tử hữu hạn - Bài 4 (Hệ thanh - Phần 1/3: Các hàm dạng và các ma trận phần tử) - Nguyễn Xuân Thành

1. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương BÀI GIẢNG 4 (hệ thanh - phần 1/3 các hàm dạng và các ma trận phần tử) Nguyễn Xuân Thành tkris1004@nuce.edu.vn Bộ môn Cơ học Kết cấu Trường Đại học Xây dựng Ngày 29 tháng 8 năm 2013
2. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương NỘI DUNG CHÍNH 1 Hàm dạng Hàm chuyển vị cho phần tử tuyến tính Một quy trình xác định hàm dạng Hàm dạng của một số phần tử thanh tuyến tính 2 Ma trận độ cứng Quy trình chung Phần tử thanh chịu kéo/nén đúng tâm Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu 3 Véc-tơ lực nút tương đương Quy trình chung Phần tử thanh chịu kéo/nén và uốn
3. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương NỘI DUNG CHÍNH 1 Hàm dạng Hàm chuyển vị cho phần tử tuyến tính Một quy trình xác định hàm dạng Hàm dạng của một số phần tử thanh tuyến tính 2 Ma trận độ cứng Quy trình chung Phần tử thanh chịu kéo/nén đúng tâm Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu 3 Véc-tơ lực nút tương đương Quy trình chung Phần tử thanh chịu kéo/nén và uốn
4. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm chuyển vị cho phần tử tuyến tính ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ ∙ Phần tử tuyến tính Khi chỉ có các điểm nút đặt ở đỉnh của phần tử ⇒ phần tử "tuyến tính"* Đơn giản cho việc thiết lập công thức các ma trận phần tử Trong các tính toán yêu cầu độ chính xác cao, ít dùng!
5. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm chuyển vị cho phần tử tuyến tính Hàm chuyển vị Thường được chọn ở dạng đa thức Trường hợp phần tử tuyến tính, chọn các hàm chuyển vị là các đa thức có số thành phần đúng bằng số chuyển vị tại nút của phần tử . Với phần tử thanh tuyến tính, bài toán 1 chiều: ( ) = 1 + 2 Với phần tử 2D tam giác tuyến tính: ( , ) = 1 + 2 + 3 ( , ) = 4 + 5 + 6
6. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm chuyển vị cho phần tử tuyến tính Hàm chuyển vị Với phần tử 2D hình chữ nhật tuyến tính: ( , ) = 1 + 2 + 3 + 4 ( , ) = 5 + 6 + 7 + 8 Với phần tử 3D hình chóp tuyến tính: ( , , ) = 1 + 2 + 3 + 4 ( , , ) = 5 + 6 + 7 + 8 ( , , ) = 9 + 10 + 11 + 12
7. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm chuyển vị cho phần tử tuyến tính Hàm chuyển vị Với phần tử 3D hình hộp tuyến tính: ( , , ) = 1 + 2 + 3 + 4 + 5 + ··· + 8 ( , , ) = 9 + 10 + 11 + 12 + 13 + ··· + 16 ( , , ) = 17 + 18 + 19 + 20 + 21 + ··· + 24 Nhận xét: Điều kiện tương thích của trường chuyển vị trên các biên chung của các phần tử hữu hạn đòi hỏi số các thành phần của mỗi hàm chuyển vị phải bằng số nút của của phần tử hữu hạn và ngược lại. Chú ý! Sự có mặt của số hạng kiểu như hoặc , nhưng dạng hình học lại đặc biệt, nên
8. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Một quy trình xác định hàm dạng Hàm dạng là hàm thể hiện mối liên hệ giữa chuyển vị tại một điểm bất kỳ bên trong phần tử với các chuyển vị tại nút của phần tử. Như vậy, để xác định hàm dạng, ta xuất phát từ hàm chuyển vị, kết hợp với điều kiện biên động học của phần tử (tức là các chuyển vị tại nút của phần tử). Có thể viết trường chuyển vị (tập hợp các hàm chuyển vị) như sau: u( , , ) = A( , , )α trong đó u là véc-tơ các chuyển vị thành phần của điểm ( , , ) bên trong phần tử; A( , , ) là ma trận chứa các hàm cơ sở 푖( , , ); và α là véc-tơ các hệ số 푖 trong các hàm chuyển vị.
9. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Một quy trình xác định hàm dạng Biểu thức trên nghiệm đúng tại các điểm, trong đó có các điểm nút của phần tử, mà tại đó: u( , , )|( , , )tại nút = q Như vậy*: −1 (︁ )︁ −1 α = A( , , )|( , , )tại nút q = A푒 q Thay vào biểu thức của trường chuyển vị, có: −1 u( , , ) = A( , , )A푒 q ⏟ ⏞ N( , , )
10. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm dạng của phần tử thanh chịu kéo/nén đúng tâm 푞1 푖 , 푗 푞4 ∙ ∙ 퐿 Hàm chuyển vị: [︂ ]︂ 1 [︀ ]︀ ( ) = 1 + 2 = A( ) với A( ) = 1 2 Xét điều kiện biên tại các đầu thanh: [︂푞 ]︂ [︂A(0)]︂ [︂ ]︂ [︂1 0]︂ [︂ ]︂ q = 1 = 1 = 1 푞4 A(퐿) 2 1 퐿 2 ⏟ ⏞ A푒
11. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm dạng của phần tử thanh chịu kéo/nén đúng tâm Từ đó, rút ra: [︂1 0]︂−1 [︁(︁ )︁ ]︁ N( ) = A( )A−1 = [︀1 ]︀ = 1 − 푒 1 퐿 퐿 퐿
12. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm dạng của phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu 푞2 푞5 푞 푞 3 , 6 ∙ ∙ 푖 푗 퐿 [︀ ]︀ Có bốn chuyển vị tại nút q = 푞2 푞3 푞5 푞6 Hàm chuyển vị: ⎡ ⎤ 1 2 3 ⎢ . ⎥ ( ) = 1 + 2 + 3 + 4 = A( ) ⎣ . ⎦ 4 với A( ) = [︀1 2 3]︀
13. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm dạng của phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu Xét điều kiện biên tại các đầu thanh: ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ ⎡ ⎤ 푞2 A(0) 1 1 0 0 0 1 ⎢푞3⎥ ⎢A, (0)⎥ ⎢ 2⎥ ⎢0 1 0 0 ⎥ ⎢ 2⎥ q = ⎢ ⎥ = ⎢ ⎥ ⎢ ⎥ = ⎢ 2 3 ⎥ ⎢ ⎥ ⎣푞5⎦ ⎣ A(퐿) ⎦ ⎣ 3⎦ ⎣1 퐿 퐿 퐿 ⎦ ⎣ 3⎦ 2 푞6 A, (퐿) 4 0 1 2퐿 3퐿 4 ⏟ ⏞ A푒 Từ đó, rút ra: ⎡1 0 0 0 ⎤−1 −1 [︀ 2 3]︀ ⎢0 1 0 0 ⎥ N( ) = A( )A = 1 ⎢ 2 3 ⎥ 푒 ⎣1 퐿 퐿 퐿 ⎦ 0 1 2퐿 3퐿2
14. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Hàm dạng của phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu Triển khai cụ thể: ⎡ 3 2 2 3 ⎤ 1 − + ⎢ 퐿2 퐿3 ⎥ ⎢ ⎥ ⎢ 2 2 3 ⎥ ⎢ − + ⎥ ⎢ 퐿 퐿2 ⎥ N( ) = ⎢ ⎥ ⎢ 3 2 2 3 ⎥ ⎢ − ⎥ ⎢ 2 3 ⎥ ⎢ 퐿 퐿 ⎥ ⎣ 2 3 ⎦ − + 퐿 퐿2 Các trường hợp phần tử khác Thực hiện theo cách tương tự trên.
15. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương NỘI DUNG CHÍNH 1 Hàm dạng Hàm chuyển vị cho phần tử tuyến tính Một quy trình xác định hàm dạng Hàm dạng của một số phần tử thanh tuyến tính 2 Ma trận độ cứng Quy trình chung Phần tử thanh chịu kéo/nén đúng tâm Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu 3 Véc-tơ lực nút tương đương Quy trình chung Phần tử thanh chịu kéo/nén và uốn
16. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Quy trình chung Từ quan hệ chuyển vị - biến dạng và các hàm dạng đã được xác định ở bước trên, ta đi tìm ma trận B liên hệ giữa biến dạng với chuyển vị tại nút của phần tử B = ∇N Từ quan hệ ứng suất - biến dạng, ta xác định ma trận E. Trong bài toán thanh, ma trận đàn hồi E có dạng E = Thay vào công thức xác định ma trận độ cứng phần tử K ∫︁ K = B EB Ω Ω
17. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu kéo/nén đúng tâm 푞1 푖 , 푗 푞4 ∙ ∙ 퐿 Do thành phần biến dạng duy nhất trong phần tử này chỉ là biến dạng dọc trục 휀 và có liên hệ với chuyển vị dọc trục ( ) theo 휀 = ( )/ nên ma trận B được xác định như sau: [︁(︁ )︁ ]︁ [︂ 1 1 ]︂ B = N = 1 − = − 퐿 퐿 퐿 퐿
18. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu kéo/nén đúng tâm Trong bài toán kéo/nén đơn, đúng tâm, không có các ứng suất trước và biến dạng trước, thì quan hệ giữa ứng suất và biến dạng tuân theo Định luật Hooke quen thuộc sau: 휎 = 휀 Như vậy E = , từ đó dẫn đến: 퐿 ⎡ 1 ⎤ ∫︁ ∫︁ − [︂ ]︂ 퐿 1 1 K = B B Ω = ⎢ ⎥ − ⎣ 1 ⎦ 퐿 퐿 Ω 0 퐿 [︂ 1 −1]︂ hay K = 퐿 −1 1
19. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu 푞2 푞5 푞 푞 3 , 6 ∙ ∙ 푖 푗 퐿 Xét dầm Euler-Bernoulli, bỏ qua thành phần biến dạng trượt, lúc này biến dạng 휀 ở tiết diện có tọa độ , tại một phân tố cách trục trung hòa đoạn bằng như sau: 휀 ( , ) = = − , ( ) 휌
20. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu Như vậy, ta có ma trận B được xác định như sau: ⎡ 3 2 2 3 ⎤ 1 − + ⎢ 퐿2 퐿3 ⎥ ⎢ ⎥ ⎢ 2 2 3 ⎥ 2 2 ⎢ − + ⎥ ⎢ 퐿 퐿2 ⎥ B = − N = − ⎢ ⎥ 2 2 ⎢ 3 2 2 3 ⎥ ⎢ − ⎥ ⎢ 2 3 ⎥ ⎢ 퐿 퐿 ⎥ ⎣ 2 3 ⎦ − + 퐿 퐿2
21. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu Cụ thể: ⎡ 3 12 ⎤ − 2 + 3 ⎢ 퐿 퐿 ⎥ ⎢ 4 6 ⎥ ⎢ − + ⎥ ⎢ 2 ⎥ ⎢ 퐿 퐿 ⎥ B = − ⎢ ⎥ ⎢ 6 12 ⎥ ⎢ 2 − 3 ⎥ ⎢ 퐿 퐿 ⎥ ⎣ 2 6 ⎦ − + 퐿 퐿2 Dễ dàng xác định ma trận E =
22. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu Thay vào biểu thức xác định ma trận độ cứng phần tử K Thành phần ở hàng 푖 cột 푗 với (푖, 푗 = 1, , 4) của ma trận độ cứng này được xác định như sau: 퐿 ⎛ ⎞ ∫︁ ∫︁ 2 퐾푖푗 = ⎝( 푖), ( 푗), ⎠ 0 trong đó 푖( ) và 푗( ) tương ứng là thành phần thứ 푖 và thứ 푗 của véc-tơ các hàm dạng N
23. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu ∫︁ Dễ dàng nhận thấy 2 là mô-men quán tính của mặt cắt ngang xung quanh trục Thay vào trên và tiến hành các phép lấy đạo hàm và tích phân cần thiết, ta có kết quả cuối cùng sau: ⎡12 6퐿 −12 6퐿 ⎤ 2 2 ⎢ 4퐿 −6퐿 2퐿 ⎥ K = ⎢ ⎥ 퐿3 ⎣ 12 −6퐿⎦ đx 4퐿2 Các trường hợp phần tử khác Thực hiện theo cách tương tự trên.
24. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương NỘI DUNG CHÍNH 1 Hàm dạng Hàm chuyển vị cho phần tử tuyến tính Một quy trình xác định hàm dạng Hàm dạng của một số phần tử thanh tuyến tính 2 Ma trận độ cứng Quy trình chung Phần tử thanh chịu kéo/nén đúng tâm Phần tử thanh chịu uốn ngang phẳng, nút cứng hai đầu 3 Véc-tơ lực nút tương đương Quy trình chung Phần tử thanh chịu kéo/nén và uốn
25. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Quy trình chung Thay ma trận các hàm dạng N và ma trận B liên hệ giữa biến dạng với chuyển vị tại nút phần tử đã tìm được ở trên, cùng với trường tải trọng đã được cho trước vào công thức xác định lực nút tương đương R ∫︁ ∫︁ R푖 = B Eε0 Ω − B 휎0 Ω ∫︁ ∫︁ + N f Ω + N t 푆 + P Nhận xét: Không cần sử dụng đến ma trận B nếu nguyên nhân tác dụng không bao gồm ứng suất trước và biến dạng trước.
26. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu kéo/nén và uốn 푅2 ( ) 푅5 푅 푅 푅 3 , , 6 푅 1 ∙ ∙ 4 푖 푗 ( ) 퐿 Do phần tử chịu đồng thời cả kéo/nén và uốn, nên ta xét trường hợp tổng quát là mỗi đầu nút phần tử có 3 chuyển vị tại nút. Như vậy, véc-tơ chuyển vị tại nút tổng cộng có 6 thành phần [︀ ]︀ q = 푞1 푞2 . . . 푞6
27. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu kéo/nén và uốn Tương ứng với nó là véc-tơ lực nút tương đương cũng có 6 thành phần [︀ ]︀ R = 푅1 푅2 . . . 푅6 Chuyển vị tại một điểm bất kỳ có tọa độ : [︂ ( )]︂ u = ( ) Ma trận các hàm dạng sau khi kể đến tất cả các chuyển vị tại nút: [︂ ( ) 0 0 ( ) 0 0 ]︂ N = 1 4 0 2( ) 3( ) 0 5( ) 6( )
28. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu kéo/nén và uốn Trường tải trọng tác dụng: [︂ ( )]︂ p = ( ) Do không xét đến các nguyên nhân ứng suất ban đầu và biến dạng ban đầu, công thức xác định lực nút tương đương sẽ chỉ là: ∫︁ R = N p Ω Ω
29. Hàm dạng Ma trận độ cứng Véc-tơ lực nút tương đương Phần tử thanh chịu kéo/nén và uốn Cụ thể các thành phần của véc-tơ tải trọng như sau: 퐿 ∫︁ 푅1 = 1( ) ( ) 0 퐿 ∫︁ 푅2 = 2( ) ( ) 0 và cứ như vậy Một số trường hợp đặc biệt của tải trọng