Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008 192
NGHIÊN CỨU TÍNH TOÁN KẾT CẤU BẰNG PHẦN MỀM
MATLAP THEO PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN
STUDING ABOUT CONSTRUCTION CALCULATION BY MATLAB
SOFTWARE ON LIMITED ELEMENT METHOD

SVTH: NGUYỄN THÀNH NHẬT
Lớp 03X3B, Trường Đại học Bách khoa, Đại học Đà Nẵng
GVHD:ThS NGUYỄN LAN
Trường Đại học Bách khoa, Đại học Đà Nẵng

TÓM TẮT
Ngày nay,nhiều phương pháp tính số đã và đang phát triển mạnh mẽ và trở thành một công
cụ hữu hiệu không thể thiếu được khi giải quyết các bài toán Khoa học – Kỹ thuật. Trong đó
phương pháp phần tử hữu hạn đã trở thành công nghệ phần mềm phổ biến và hiệu quả .Công
việc này không thể tách rời một công cụ đồng hành là ngôn ngữ lập trình kỹ thuật Matlap.
Nghiên cứu này giới thiệu việc ứng dụng ngôn ngữ Matlap để phân tích cho tải trọng di động
thông qua đường ảnh hưởng, tính toán nội lực của dầm liên tục.
ABSTRACT
Nowaday, number calculating methods were and developing strongly, became a effective tool
can’t short when solving scientis-technich mathematic excercises.in where, limited element
method became a popular and effective process software. This work can’t separated with a
tool parner is matlab technich program languages. This study introduce the language applying
to analysic for move load throught influenced line, interforce caculating of continuous beam.

1. Đặt vấn đề
1.1. Lý do chọn đề tài
Như đã biết, lý thuyết phần tử hữu hạn đã được ứng dụng vào lĩnh vực tính kết cấu

(t), U
’
(t), U(t), F(t): Véc tơ gia tốc, vận tốc, chuyển vị nút và véc tơ tải trọng thay
đổi theo thời gian. Các ma trận độ cứng, khối lượng, ma trận cản đều là ma trận vng đối
xứng, chúng được lắp ghép từ các ma trận tương ứng của từng phần tử trong kết cấu.
– Trường hợp phân tích tĩnh (Static Analysis): F(t)= F
Phương trình (1) trở thành: K. U = F (2)
Giải hệ phương trình (2) tìm tất cả các thành phần chuyển vị tại các nút, sau đó tính nội
lực ứng suất cho từng phần tử.
– Trường hợp phân tích tần số dao động riêng (Eigen value Annalysis):
Khi tải trọng ngồi bằng zero, bỏ qua lực cản của mơi trường lúc đó kết cấu dao dộng
điều hòa chuyển vị của hệ có dạng:
U=U. sin(t) và U
’’
= -U. 
2
. sin(t) (3)
-M.U. 
2
. sin(t) + K. U. sin(t) = {0} (K - 
2.
M). U = {0} (4)
Giải phương trình (4) bằng phương pháp SUBSPACE sẽ cho các giá trị riêng và véc tơ
riêng từ đó tính được các tần số riêng (eigen frequencies) và dạng dao động riêng (mode
shape) tương ứng.
3. Ví dụ phân tích cầu dầm liên tục năm nhịp
3.1. Giới thiệu tổng quan về cầu dầm liên tục năm nhịp
Kết cấu cầu chính loại cầu dầm liên tục năm nhịp 40m+3x70m+40m ...
3.2. Mơ hình phần tử hữu hạn và các yếu tố được xem xét
Kết cấu được mơ hình hố theo phương

ẠP ÂÀÛT CẠC ÂIÃƯU KIÃÛN RNG BÜC
CHO MA TRÁÛN V VẸC TÅ CA HÃÛ
TÊNH CẠC MA TRÁÛN V CẠC VẸC TÅ
CHO MÄÙI PHÁƯN TỈÍ
NHÁÛP VO SÄÚ LIÃÛU BAN ÂÁƯU
ÂÀÛC TRỈNG HÇNH HC, VÁÛT LIÃÛU, TI TRNG
BÀÕT ÂÁƯU
V V IN CẠC KÃÚT QU THEO MONG MÚN
TÊNH CẠC BIÃÚN THỈÏ CÁÚP TỈÌ CHUØN
VË NỤT(MÄMEN,LỈÛC CÀÕT...)
KÃÚT THỤC
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008 194

Hình 4.1: Đường ảnh hưởng momen tại tiết diện cách đầu dầm một đoạn 20m

Hình 4.2: Đường ảnh hưởng lực cắt tại tiết diện cách đầu dầm một đoạn 20m

Hình 4.3: Biểu đồ mômen do tĩnh tải giai đoạn 1 gây ra Hình 4.4: Biểu đồ lực cắt do tĩnh tải giai đoạn 1 gây ra
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008

195

Hình 4.5: Biểu đồ mômen do tĩnh tải giai đoạn 2 gây ra


luận sau:
- Mặt dù việc tính toán nội lực của cầu dầm liên tục còn nhiều vấn đề phải giải quyết
như co ngót, từ biến… Tuy nhiên trên trên cơ sở những bài toán cơ bản đã được giải như bài
toán tính nội lực do hoạt tải , bài toán nội lực do tĩnh tải … sẽ giúp ích rất nhiều cho việc thiết
kế các bài toán cầu. Trong thời gian có hạn việc tính toán không tránh khỏi thiếu xót. Trong
thời gian tới tác giả sẽ cố gắng hoàn thiện và phát triển hơn nữa việc ứng dụng ngôn ngữ lập
trình Matlap để giải quyết hơn nữa các bài toán về chuyên ngành cầu như các bài toán về phân
tích đông, các bài toán về phi tuyến hình học. Kết quả tính toán cầu dầm liên tục bằng lập trình
Matlap có thể sử dụng được để tính toán các bài toán mà thực tế yêu cầu.

TÀI LIỆU THAM KHẢO

[1] Võ Như Cầu (2005), Tính kết cấu theo phương pháp phần tử hữu hạn, NXB Xây dựng,
Hà Nội.
[2] Chu Quốc Thắng (1997), Phương pháp phần tử hữu hạn, NXB Khoa học và Kỹ Thuật, Hà
Nội .
[3] Lều Thọ Trình (2005), Cơ học kết cấu tập 2, NXB Khoa học và Kỹ Thuật, Hà Nội .
[4] Nguyễn Hoài Sơn, Lê Thanh Phong, Mai Đức Đãi (2008), Ứng dụng phương pháp phần
tử hữu hạn trong tính toán kết cấu , NXB Đại học Quốc gia TP. Hồ Chí Minh, TP.Hồ Chí
Minh .
[5] Nguyễn Hoài Sơn, Vũ Như Phan Thiện, Đỗ Thanh Việt (2001), Phương pháp phần tử hữu
hạn với Matlap , NXB Đại học Quốc gia TP. Hồ Chí Minh, TP.Hồ Chí Minh .
[6] Nguyễn Hoàng Hải, Nguyễn Khắc Kiểm, Nguyễn Trung Dũng, Hà Trần Đức (2003), Lập
trình Matlap, NXB Khoa học và Kỹ Thuật, Hà Nội .
[7] Calfem.