191
1
2
3
4
5
1
0
0 0
0
0
0
0
0 1
150
0
0
2
0 0
0
0
1 0
0
0
0 0
S
c
S
s
s
S
c
c
S
s
s
S
=
−
với
sin ,
cos
s
c
θ
θ
=
=
và các ẩn cần tìm là
i
S
(
1,.., 5)
i
=
. Hãy viết chương trình Matlab ñể
tính các lực này, khi cho trước góc
θ . Chạy
chương trình với
56
o
θ
=
.
Triển khai trong Matlab như sau:
function Vidu92
disp('Vao goc theta ()so duong n: ')
theta = input('theta = ')
c = cosd(theta); s = sind(theta);
A = [c 1 0 0 0;
0 s 0 0 1;
0 0 2*s 0 0;
0 -c c 1 0;
0 s s 0 0];
b = [0; 0; 150; 0; 0];
x=A^-1*b;
disp('Ung luc trong cac thanh: ')
x
Chạy chương trình và nhập theta = 56 (ñộ), ta nhận ñược kết quả, ñó là ứng lực
trong các thanh.
x =
161.7802
-90.4663
90.4663
-101.1763
75.0000
Hệ thanh siêu tĩnh
ðối với các hệ thanh siêu tĩnh, do có các liên kết thừa, các phương trình cân bằng
nút không ñủ ñể ta giải tìm ra ứng lực trong các thanh. ðể giải ñược bài toán này
chúng ta cần phải chú ý thêm các ñiều kiện ràng buộc, các ñiều kiện liên kết. Các
ñiều kiện này sẽ thể hiện sự phụ thuộc lẫn nhau biến dạng dài của các thanh. Dưới
ñây là một ví dụ.
θ
1
S
2
θ
1
θ
1
θ
1
P=150N
S
3
S
3
S
4
S
4
S
2
S
1
S
1
S
5
S
5
Hình 9-2
