Cỏc bước thiết lập chương trỡnh tớnh cho hệ khung cú Ne thanh với R nỳt: Bước 1 Lập bảng đỏnh số thanh nỳt : cú Ne thanh với R nỳt

Một phần của tài liệu Bài giảng phương pháp số (Phương pháp phần tử hữu hạn) (Trang 84 - 86)

- hàm dạng của phần tử thanh chịu kộo , nộn dọc trục

3.Cỏc bước thiết lập chương trỡnh tớnh cho hệ khung cú Ne thanh với R nỳt: Bước 1 Lập bảng đỏnh số thanh nỳt : cú Ne thanh với R nỳt

Bước 1. Lập bảng đỏnh số thanh- nỳt : cú Ne thanh với R nỳt

Nỳt T.T

thanh Đầu : i Cuối : j

1 1 2

... ... ...

k ki kj

... ... ...

Ne Nei Nej

Ở bảng này đỏnh số và nhập bằng tay ( nếu số lượng ớt) hoặc viết một đoạn chương trỡnh để tự động đỏnh số thứ tự cỏc thanh và cỏc nỳt. Số liệu ghi vào mảng được khai bỏo :

Bước 2. Lập bảng nhập cỏc thụng số của cỏc thanh dầm : (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) TT thanh L E G Jx Jy Jz x  mx nx y my ny z mz nz 1 ... k ... Ne

Số liệu ở bảng này được ghi vào mảng :

Thanh(1 to Ne, 1 to 15) as double

Chỳ ý : Cỏc trục chớnh x , y , z ở cỏc thanh cú gốc là nỳt đầu và phải lập thành tam giỏc thuận.

Bước 3. Lập bảng ma trận chỉ số [b] : Thanh – Nỳt - Bậc tự do : Nỳt đầu : i Nỳt cuối : J T.T thanh q6i -5 q6i -4 q6i -3 q6i -2 q6i -1 q6i q6j-5 q6j-4 q6j-3 q6j-2 q6j-1 q6j 1 ... k ... Ne

Số liệu sẽ tự tớnh và ghi vào mảng :

b(1 to Ne , 1 to 12) as integer

Bước 4. Tớnh cỏc giỏ trị của [T](e) được ghi vào mảng: T(1 to 12 , 1 to 12) as double

Bước 5. Tớnh [ K](e) được ghi vào mảng :

Ke(1 to 12, 1 to 12) as double

KeT(1 to 12, 1 to 12) as double

Bước 7. Ghộp [K'](e) vào ma trận cứng tổng thể [K'] được ghi vào mảng: KT ( 1 to N, 1 to N) as double

Ở đõy N = 6 . R tổng số bậc tự do

Bước 8. Tớnh vộc tơ tải phần tử => vộc tơ tải tổng thể

Bước 9. Áp cỏc điều kiện biờn

Chỳ ý rằng sau khi cú được ma trận cứng tổng thể [K'] ( mà định thức của nú = 0) , ta sẽ tớnh cho vộc tơ tải tổng thể trờn cơ sở tớnh cỏc vộc tơ tải phần tử. Nhờ cỏc điều kiện biờn ta sẽ bỏ đi được một số phương trỡnh và một số biến ( số lượng của chỳng bằng nhau), điều đú cú nghĩa là ma trận cứng tổng thể [K'] sẽ bị bỏ đi một số hàng và một số cột và sẽ thành ma trận cứng tổng thể [K'* ] . Ma trận [K'* ] là ma trận vuụng cú cấp nhỏ hơn N và định thức của nú  0 , do vậy sẽ dẫn về một hệ Cramer để giải cỏc biến là cỏc bậc tự do ma ta cần tỡm.

Bước 10. Giải hệ phương trỡnh tớnh cỏc chuyển vị qk ( giỏ trị cỏc bậc tự do theo hệ tọa độ tổng), từ đú tớnh {q}e = [T]e. {q’}e

Bước 11. Tớnh cỏc đại lượng : biến dạng , ứng suất .

Một phần của tài liệu Bài giảng phương pháp số (Phương pháp phần tử hữu hạn) (Trang 84 - 86)