v

MỤC LỤC
 TRANG


 Error! Bookmark not defined.
L ii
 Error! Bookmark not defined.
 Error! Bookmark not defined.
 iix
 ixi
CHƯƠNG 1: GIỚI THIỆU Error! Bookmark not defined.
1.1 TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CU 1
1.2 CÁC KẾT QUẢ NGHIÊNCU TRONG VÀ NGOÀI NƯỚC 2
1.2.1 N 2
1. 6
1.3 MỤC TIEU NGHIEN CU CỦA DỀ TAI 8
1.4 NHIỆM VỤ CỦA ĐỀ TÀI VÀ GIỚI HẠN ĐỀ TÀI 8
1.5 PHƯƠNG PHÁP NGHIÊN CU 8
CHƯƠNG 2: CƠ SỞ LÝ THUYẾT 9
2.1  9
2.1.1  9
2.1.2  9
2.1.3  Error! Bookmark not
defined.0
2.1.4  Error! Bookmark not defined.1
2.1.5  Error! Bookmark not defined.2
2.2   . Error! Bookmark not
defined.3
vi

vii

5.1  33
5.2  34
5.3  35
5.3.1  35
5.3.2 Ma trcicacutrúcmtcu 37
5.3.3 Ma trcngtrnglccacutrúcmtcu 38
5.3.4 Ma trcicacáp 39
5.3.5 Ma trcngtrnglccacáp 39
5.3.6 Ma trnkhngnhtcacutrúcmtcu 40
5.3.7 Ma trnkhngcacáp 41
5.3.8  41
5.4  41
5.4.1  41
5.4.2  42
CHƯƠNG 6 PHÂN TÍCH ẢNH HƯỞNG CỦA ĐỘNG ĐẤT ĐỐI VỚI CẦU
CÁP TREO 44
6.1  44
6.2  45
6.2.1   45
6.2.2  46
6.2.3  47
6.3 
 48
6.3.1  48
6.3.2  49
6.4 BÀI TOÁN ÁP DNG CHO NG CN
CU TREO 50
6.4.1  50

DANH MỤC CÁC HÌNH
Hình 1.1  6
Hình 2.1  10
Hình 2.2  10
Hình 2.3  10
Hình 2.4  11
Hình 2.5  11
Hình 2.6  14
Hình 3.1  17
Hình 3.2  17
Hình 3.3  18
Hình 3.4  23
Hình 4.1  26
Hình 4.2  32
Hình 5.1  35
Hình 5.2 
ngang 36
Hình 5.3  43
Hình 6.1 
 45
Hình 6.2  47
Hình 6.3  47
Hình 6.4 n hóa bng cách áp dng mô hình mt
bc t do khi xét mt phn t 48
Hình 6.5 
 50
x

võng Tacoma Narrows, xây d phá hy ch 
vào s dc ct di la mng manh ri
gãy vn. Mu th nghi sau Tacoma, mi công
trình cu lu phi qua khâu th nghim ti trng gió ht sc ngt nghèo.
Trong các phân tích nh, các dng t do ca cu treo
có th c phân thành: chuy ng th ng, chuy ng xon và chuyn
ng hai bên. Và c s dng ngày càng ph bin và m
ng s d  xây dng mô hình lý thuyt.
n t hu hn (Finite Element Method  FEM) là m
hiu qu  gii quyt nhic c. Vic tính toán trên máy tính
phân tích phn t hu hr
phân tích v ng hc ca cu treo. n t hu hn là chia min
kho sát thành mt s hu hn các min, gi là các phn t. Công vic
 c gi là ri rc hóa min kho sát, hay ti phn thu hn,
-2-

c áp di vi cu trúc phc tu treo. Vic tính toán trên máy tính
t qu ng ca cu trúc mà có th 
bc t doc thit k và ng d c tính ca
ng hc, các tn s t nhiên, các dng cng và kh  
ng ca các thành phn khác nhau ca cu treo.
Trong lui hng dn và hc viên s dn
t hu hn (FEM) ving ngôn ng Matlab nhm Nghiên cứu động
lực học của cầu cáp treo nhịp dài

Tóm li: Nghiên cu v v ng lc hc ca c
ln lao không ch c. Vi mong muc
nghiên cu và phát trin các v ng lc hc cu cáp treo  Vit Nam bng
i; ng dn và h tài :”Nghiên cứu động
lực học của cầu cáp treo nhịp dài bằng phương pháp phần tử hữu hạn ”.

ng vi các tn s  chính xác cao, còn t l gim xóc
t chính xác.
Robert JamesWestgate [15] nghiên cu ng ca cu treo v nhi
và ti trng. Khi nhi ca cu trúc ct cu b chy và giãn n
ng thi tn s t m trong ngàyi trng thì tn
s t nhiên gimc bit trong các gi m. T c s i v
dng tn s và khng ph thuc vào s i dòng giao thông trên cu, mà
có th m các dng ng ca cu.
Theo M. Zribi, N. B. Almutairi, M. Abdel-Rohman [16]  linh hot ca
cáp và gim chn trên cu trúc cu treo thp khin cu d b ng do gió và ti
trng. Bài báo nghiên cu vic u khing ca cu cáp treo da
trên ti trng di chuyn trên mt cu vi t i,
cáp gia mt cu và cc s dng thit b thy lcu khin
nhm gi      ng, mà thit k da trên thuyt
m bu cng vng. Kt qu cho thu khin làm vic
tt
-4-

Trong nghiên cu ca J.D. Yau, L. Fryba [17], c
dc treo mt nhp, ng ca dm vi thi gian ph thuc vào
u kin biên  ng ca d c chia thành hai phn: ph   n
nglc hc d. T ng ca sóng
a chn dc chiu dài nhp cu, nhp gia ca dm treo không  v trí quan trng
cng gia tc ci. Các di vi  gia tc ci

 Các nghiên cu v ng ca t i vi cu cáp treo
Thm ha t p rt nhiu
các d liu v s  sp ca c bo v tu treo
i ng ci
thit k cu bng cách quan sát phm vi s dng và m rng ra t nhng nghiên


Hình 1.1 Cầu Gavin Canyon đổ sập năm 1994 dưới ảnh hưởng của động đất

1.2.2  c
m  c ta là có nhiu sông rng, bin ln, vc áp dng
kt cu cu cáp treo là mt trong nh
xây d h tng hic nghiên cu tính toán
kt cu cu cáp treo  c nhiu và luôn là bài toán khó và vic t
ng hóa tính toán càng phc tCác công trình nghiên cu  Vit Nam ch
-6-

yu ch dng li  mc nghiên cu ng dng trên lý thuyt vì thit 2 yu t:
i ca th ng và kinh phí trin khai nghiên cu ng dng. Mt s nghiên
cu  Vi
 ng Buffeting là v quan tâm nht khi thit k
cu treo nhp dài, c Huynh, Nguyn Ngc Trung [4] u
khing lc hc bng cách s dng flaps. c gn  hai phía ca mt
cu c dng theo s chuyng ca mt cc di chuyn,
b mt ct ngang ca mt cu chi mt cách liên tc. Kt qu là
lng lc hi. H s ca h thu khic nghiên cu
sut trong quá trình nghiên cu thc nghim và x lý s liu.
Phan c Huynh, Phm Thanh Hoàác gii pháp thay th cho
tính nh cng và rung lc ca cu treo dài s u khin bng
cách s dng Flaps. m cng trng nh mà không có
 cng b  nh cng hc. B u khin phn hi
thông tin bng cách s dng winglets cho vi ng và gim phn
ng rung lc ca cc nghiên cu thông qua mô hình mt cu hai chiu.
Kt qu cho thy rng t n ng rung lc gim
thông qua chuyng kim soát thích hp ca winglets. Ngoài ra nghiên cc
thc hin trên cu treo 3000m.


1.U
- Da vào lý thuyt v phân tích phân t hu hn ca c tính toán.
- Tham kho tài lic v 
-  kim tra kt qu tính toán, mô phng kt qu.
-9-

2
 LÝ THUYT

2.1 GII THIU V PHN T HU HN
2.1.1 Gii thiu chung
  n t hu hn (Finite Element Method  FEM) là mt
 c bit có hiu qu  tìm dng ga mt 
bit trong min xác nh V ca nó. Tuy nhiên, FEM không tìm dng xp x ca n
hàm trên toàn min V ca kt cu mà ch tìm trong tng min con Ve. Chính vì vy
mà FEM có th áp dng cho rt nhiu bài toán k thut và nht i vi bài toán
kt cn hàm cn tìm có th nh trên các min phc tp vi
nhiu kin biên khác nhau.
i vi FEM mic thay th bi mt s hu hn
các mic gi là phn t. Các phn t này ch c ni vi nhau bi
c trên biên gi là nút. Trong phm vi mi phn tng cn
c xp x theo mt dng phân b  s ca hàm xp x
c gi là các tham s hay các t tng quát. Các tham s này lc biu

Hình 2.2: Phần tử một chiều Phn t bc nht Phn t bc hai Phn t bc ba
Hình 2.3: Phần tử hai chiều
-11- Phn t bc nht Phn t bc hai Phn t bc ba
Hình 2.4: Phần tử ba chiều

Phn t bc nht Phn t bc hai Phn t bc ba
Hình 2.5: Phần tử lăng trụ

2.1.4 Tính toán bng phn t hu hn
Cn t hu hn gc chính áp dng cho tính toán trên
cu treo sau [3]:
- c các d liu vào. Các d liu này gm các thông tin mô t
nút và phn t i phn t), các thông s c ca vt li
hi, h s dn nhi ti trng tác dng và thông tin v
liên kt ca kt cu kin biên).
- c 2: Tính toán ma tr cng phn t k, ma trn khng m 
lc nút phn t f ca mi phn t.
- c 3: Xây dng ma trn khi tng th K, ma trn khng M 
lc nút F chung cho c h (ghép ni phn t).
- u kin liên kt trên biên kt cu bng cách bii
ma trn K, M c nút tng th F.
- c 5: Gi   n t hu h  nh giá tr riêng, 
riêng, tn s t nhiên và các dng
- c 6: t qu, v  th biu din các dng cng.

ng, ví d nguyên lý công kh , FEM cung cp m tng quát mang
-13-

tính trc quan theo quy lut t t yêu cu li vi nhng k t
cu.

2.2 TNG QUAN V NGÔN NG LP TRÌNH
2.2.1Lch s i ca Matlab
Matlab vit tt t  » n v ma trnc Cleve
Moler phát minh vào cui thp niên 19c vit bi ngôn
ng t b phn dùng trong ni b ci hc
Stanford.
t li Matlab bng ngôn ng c xây dng
n phc v cho thit k h thu khin, mô phng, . . .Jack xây
dng MATLAB tr thành mô hình ngôn ng l ma trn (matrix-
based programming language).
t trình thông d
Little kt hp vi Moler và Steve Bangert quy án
i -  [1], [12].

2.2.2 Khái nim v Matlab
Matlab là ngôn ng lp trình thc hành bc s d gii các bài
toán v k thut. Matlab tích hp trong vic tính toán, th hin kt qu, lp
trình D liu cùng vn lp trình si s dng thc hin các
ng dng sau [12]:
- S dng các hàm có sn, cho phép tính các phép toán thông
ng
- Cho phép lp trình to ra nhng ng dng mi
- Cho phép mô phng các mô hình thc t
- Phân tích, kho sát và hin th d liu

PHÂN TÍCH NG THNG CA CU
CÁP TREO

3.1GII THIU
Thm ha ca c       ng m n
nghiên cu ca các nhà khoa hc v cu cáp treo. Mng ca
cng d di ca cu gây ra
không ch bi gió mà còn do ti trng lc h di chuyn ca xe c và
t.
Trong nhng t do ca cc
áp dng tính toán trên máy tính dn t hu hn  xác nh
tn s t nhiên và các dngv chính xác cao. Hu ht các tính toán
c da trên cng hóa sut h thng ca khi
din ca cu trúc thc t.
ng thng ca cu treo áp d
pháp phn t hu h nh tn s t nhiên, dng cng
ng hc ca cu treo [8].

3.2 PHÂN TÍCH CU CÁP TREO B
3.2.1 Gi thiu
Trong quá trình nghiên cu ng ca các nhân t khác nhau vào vic
ng hc ca cu treo có các gi nh tng quát sau [6]:
1- Tt c ng sut trong cnh lut Hooke
2- Ti tru ca cáp không gây ra ng sut trong dm
u kin này có th c mt cách tng quát khi n kt cu.
Tuy nhiên nu cc xây dng theo cách mà ti trng sut
-16-

trong dm cu này có th c cân bng trong v ng lc hc bng
cách nghiên cu tt c các ng sut bên trong liên quan.


(

, ) (3.1)
6- Trong thuyt tuyn tính, v trí trng thái cân b  c gi nh là
chuyn v ng nh
7- Lc i tác dng ca ti trng H(t) (hình 3.1-a) là nh khi so
sánh vi lc i tác dng ca ti trH
w
 lc dc
trng (hình 3.1-b) t là nh khi so sánh
vi ti tr
w
.
8-  cou ca dm cc xem là nh khi so sánh v cong ca
cáp, vì vc b qua khi tính toán.
-17-

(a)
(b)
Hình 3.1 – (a) Cấu hình động học của cấu trúc tăng cng, (b) Cấu hình động học
của tháp đỡ

3.2.2u trúc ma trn
n t hu hn, cu trúc cc gi nh chia
thành mt h thng các phn t ri rc kt ni vi nhau ch ti mt s hu
hm nút và trình bày trong (hình 3.2) [4, 6]

c khi Hermitian c 3.2)
   
   
3
2
2
12
3
2
2
11
3
3
2
2
02
3
3
2
3
01
1112
2323
1
x
L
x
L
xlx
L

có th biu din trong gii hn
ca chuyn v nút
                 
tqxltqxltqxltqxltxv
e 412302211101
, 
(3.3)
T c chun hóa

1
() = 1- / và 
2
() = / (3.4)
Vì v3.3) có th c vii dng ma trn
-19-

     
 
  
tqLLtv
e 2
2
11
2
22
2
11
2
121
,23,,23,,



N
e
e
cgCG
w
e
cg
kK
LLLL
LL
LLLL
LL
L
H
k
1
22
22
433
336336
343
336336
30
(3.6)
Ma tr ci ca cáp
 
   


H
w
L
AE
K

**
(3.7)

 
   
 
T
e
N
e
eN
LLLLfff 12/2/12/2/,

22
1





Ma tr ci ca cu trúc cng
 
 
   

ge
kK
LLLL
LL
LLLL
LL
L
IE
k
1
22
22
3
4626
612612
2646
612612
1



(3.8)
Vi:
)/(12
2
LGIE
vegegege


(3.9)















N
e
e
ege
e
e
mM
LLLL
LL
LLLL
LL
L
rm
LLLL
LL
LLLL




= 0 (3.11-a)
-20-

Ta có th gi thin rng có dng  và thay th 
trình (3.11- :






2
+





= 0 (3.11-b)
Giá tr c t 
   
 
 
0
2
 rMK






=

0

(3.14)
Vi :



là ma trn thc,  là giá tr riêng,



là ma tr ,




Nhân ma trn



1
vào hai v 3c:
(


Các thông s chính
Ký hiu
Giá tr

Chiu dài nhp gia
L2
460
m
Moment quán tính ca
cu trúc cng
Ig2
3629
m
2