Slide Động lực học công trình_chương 2

CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.1 Phương trình vi phân tổng quát:
P(t)
mn

mk

M(t)

mi
m2
m1

yk (t)
yi (t)

yn (t)

y2 (t)
y1 (t)

(a)
Xét dao động của khung không trọng lượng mang các
khối lượng tập trung (hình a). Chịu các lực kích thích thay
đổi theo thời gian. Bỏ qua biến dạng dọc của khung, ta
có bài toán dao động có n bậc tự do.


CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.1 Phương trình vi phân tổng quát:
P(t)

mn

mk

M(t)

mi
m2
m1

yk (t)
yi (t)

yn (t)

y2 (t)
y1 (t)

(a)

Zk(t)

P(t)

M(t)

R2(t)

Rk(t)
Zi(t)

Z2(t)

R1(t)

Z1(t)

Ri(t)

Zn(t)

Rn(t)

(b)

Xét tại thời điểm bất kỳ t dưới tác dụng của các lực:
* Lực kích thích: M(t), P(t), q(t).
Zk ( t ) = −mk .yk ( t )
* Lực quán tính:
* Lực cản: Rk(t)


CHƯƠNG 2:
2.1 Phương trình vi phân tổng quát:
P(t)
mn

mk

M(t)

mi
m2
m1

(a)

yk (t)
yi (t)
y2 (t)
y1 (t)

yn (t)

Zk(t)

P(t)

M(t)

R2(t)

Rk(t)
Zi(t)
Z2(t)

R1(t)

Z1(t)

Ri(t)

(b)

Zn(t)

Rn(t)

δk

δi
δ2

δn
Zi=1

δ1

(c)

Gọi δ ki là chuyển vị khối lượng do Z = 1 tác dụng tĩnh
gây ra:
∆ kP(t) chuyển vị khối lượng mk do lực kích thích gây ra.
Xem hệ đàn hồi là tuyến tính, chuyển vị là rất nhỏ:


CHƯƠNG 2:
2.1 Phương trình vi phân tổng quát:
P(t)
mn

mk

M(t)

mi
m2
m1

(a)

yk (t)
yi (t)
y2 (t)
y1 (t)

yn (t)

Zk(t)

P(t)

M(t)

R2(t)

Rk(t)
Zi(t)
Z2(t)

R1(t)

Z1(t)

Ri(t)

(b)

Zn(t)

Rn(t)

δk

δi
δ2

δn
Zi=1

δ1

(c)

Phương trình chuyển vị của các khối lượng:

yk ( t ) = δ k1 [Z1 ( t ) − R1 ( t )]+ δ k 2 [Z2 ( t ) − R2 ( t )]+ ...

+ Zkn [Zn ( t ) − Rn ( t )]+ ∆ kP ( t ); k = 1,2 ,3,...,n

CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.1 Phương trình vi phân tổng quát:

yk ( t ) + δ k1 [ m1 y1( t ) + R1( t )] + δ k 2 [ m2 y2 ( t ) + R2 ( t )] + ... +

+ δ kn [ mn yn ( t ) + Rn ( t )] − ∆ kP ( t ) = 0 ; k = 1, 2 , ..., n.
Đây là phương trình vi phân tổng quát mô tả dao
động cưỡng bức có kể đến lực cản của hệ có bậc tự do
bằng n.


CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.2 Dao động riêng khi không lực cản:
Không kể đến lực kích thích và lực cản. Phương
trình được viết lại như sau:

m1 .δ k1 .y1 ( t ) + m2 .δ k 2 .y2 ( t ) + ... + mn .δ kn .yn ( t ) + yk ( t ) = 0
Nghiệm tổng quát thứ k của phương trình được
biểu thị dưới dạng tổng của các nghiệm riêng:

yk ( t ) =

n

n

i

i =1

∑ yki( t ) = ∑ yki Fi ( t )

yki : các hằng số chưa biết;
Fi(t): hàm số phụ thuộc thời gian t, chưa xác định.


CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.2 Dao động riêng khi không lực cản:
Với một nghiệm riêng thứ i, tại các khối lượng ta có:

y1i ( t ) = y1i Fi ( t ),
y2i ( t ) = y2i Fi ( t ),
..............................,
yni ( t ) = yni Fi ( t ).

y1i y2i
m1

yii

yki

yni

m2 mi m k m n

Điều này chứng tỏ tỷ lệ giữa chuyển vị của các
khối lượng không phụ thuộc vào thời gian t. Đường
cong tạo bởi các tung độ y1i , y2i , … là đường cong
đàn hồi của dầm và là dạng chính thứ i của dao động

riêng.


CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.2 Dao động riêng khi không lực cản:
Đạo hàm nghiệm riêng thứ i và thay vào phương
trình cơ bản, ta thu được:

( m1δ 11 − ui )
m2δ 12
m1δ 21
( m2δ 22 − ui )
D=
...
...
m1δ n1
m 2 δ n2
Trong đó:

...
mnδ 1n
...
mnδ 2 n
=0
...
...
... ( mnδ nn − ui )

1
ui = 2

ωi

Phương trình này được gọi là phương trình tần số
hoặc phương trình thế kỷ. Giải hệ này ta thu được các
giá trị ui,Từ các giá trị này ta tìm được các tần số dao
động riêng ω i (phổ tần số dao động riêng).


CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.3. Dạng chính của dao động riêng:
Đạo hàm nghiệm riêng thứ i và thay vào phương
trình cơ bản, ta thu được:

 ( t ) + y F ( t ) = 0
[ m1δ k1 y1i + m2δ k 2 y 21 + ... + mnδ kn yni ] F
i
ki i
⇒

i ( t )
F
=
Fi ( t ) −

yki
m1δ k1 y1i + m2δ k 2 y2i + ... + mnδ kn yni

Vế trái phụ thuộc vào thời gian t, vế phải chỉ phụ
thuộc vào kết cấu, vị trí và trị số các khối lượng, nên tỷ
số này là một hằng số và bằng -ω i2.

 ( t ) + ω 2 F ( t ) = 0
F
i
i i
[ m1δ k1 y1i + m2δ k 2 y21 + ... + mnδ kn yni ] ω i − yki = 0
2


CHƯƠNG 2: DAO ĐỘNG CỦA HỆ CÓ BẬC TỰ DO HỮU HẠN
2.3. Dạng chính của dao động riêng:

Như vậy đối với hệ có n bậc tự do luôn luôn tìm
được n giá trị tần số dao động riêng. Ứng với mỗi tần số
dao động riêng ω i có một dạng chính của dao động xác
định bằng các chuyển vị y1i, y2i, …, yni của các khối lượng.
Phương trình dao động của khối lượng thứ k với
tần số ω i có dạng:

yki ( t ) = yki ai sin( ωi t + ϕi )
Phương trình dao động tổng quát của khối lượng
thứ k:
n

yk ( t ) =

∑ yki ai sin(ωi t + ϕi )
i =1

Ví dụ 1: Cho dầm như hình vẽ với m1 = m2 = m. Tìm các
tần số dao động riêng?
m1
m2
Phương trình tần số cho
l/3
l/3
l/3
bài toán 2 khối lượng:

δ 12 m2
( δ 11m1 − u )
=0
δ 21m1
( δ 22 m2 − u )

Z1=1

2l / 9

Z2=1
2l / 9

2 =
⇒ u 2 − u( δ 11m1 + δ 22 m2 ) + m1m2 ( δ 11δ 22 − δ 12
) 0

33
4l

3 7l 3

5ml ; δ = δ ml
δ 11
⇒= δu221 ==
; u122 = 21 =
.
243
162EI
EI
486 486
EI EI


Ví dụ 1: Cho dầm như hình vẽ với m1 = m2 = m. Tìm các
tần số dao động riêng?
m1
m2
Tần số dao động
l/3
l/3
l/3
riêng được xác định:
Z1=1

2l / 9

1
EI
ω1 =

= 5 ,69
u1
ml 3
1
EI
ω2 =
= 22
u2
ml 3

Z2=1
2l / 9

m1

m2

m1

m2


Ví dụ 2: Tìm các tần số dao
động riêng và các dạng
dao động riêng chính của
dầm công xôn trên hình vẽ.
Cho biết EI = const.
Giải:

l

l

m2=m

m1=3m

EI

l
Z1=1

2l

( M1 )

Z2=1

( M2 )

Hệ có bậc tự do là 2, dao động không cản, phương
trình tần số dao động có dạng:
δ 12 m2
( δ 11m1 − u )
=0
δ 21m1
( δ 22 m2 − u )
2 =
⇒ u 2 − u( δ 11m1 + δ 22 m2 ) + m1m2 ( δ 11δ 22 − δ 12
) 0

Vẽ các biểu đồ mô men uốn đơn vị Z1 = 1, Z2 = 2. Xác
định các chuyển vị δ 11, δ 12, δ 21, δ 22.


Ví dụ 2: Tìm các tần số dao
động riêng và các dạng
dao động riêng chính của
dầm công xôn trên hình vẽ.
Cho biết EI = const.
l .l 2l
l3
δ 11 = ( M 1 )( M 1 ) =
=
2 EI 3 3 EI

δ 22

2l .2l 2.2l 8l 3
= ( M 2 )( M 2 ) =
=
2 EI 3
3 EI

δ 12

5l 3
= δ 21 = ( M 1 )( M 2 ) =
6 EI

1

EI
⇒ ω1 =
= 0 ,5345
, 1/ s
3
u1
ml

m2=m

m1=3m

EI
l
l

l
Z1=1

2l

( M1 )

Z2=1

( M2 )

; ω 2 = 2 ,5

EI

ml

3

, 1/ s


Xác định các dạng
chính của dao động:
( δ 11m1 − u1 )y11 + δ 12 m2 y21 = 0

δ 21m1 y11 + ( δ 22 m2 − u1 )y21 = 0
Cho y11 =1  y21 = 3
Tương ứng với ω 2, cũng
thực hiện tương tự như
trên, cho y12 = 1 ta sẽ tìm
được dạng chính thứ hai
của dao động riêng chuyển
vị tương ứng tại các khối
lượng.

m1=3m

EI
l

m2=m
l

l

Z1=1

2l

( M1 )

Z2=1

( M2 )

ω1
y11=1

ω2
y12=1

m1

y21=3
m2
m2
y22= -1


m2

m1

2EI

2EI
EI

3m

Ví dụ 3: Tìm các tần số
dao động riêng và các
dạng dao động riêng chính
của khung như hình vẽ.
Cho biết EI = 34,8.104 N.m2,
m = 1000/g.Ns2/m. m1 = 2m, m2
= m.

2,86

EI

2m

2m

Z1=1

5,2

2m

2m

3,12
2,08

Hệ có hai bậc tự do,
Phương trình tần số có dạng:

δ 12 m2
( δ 11m1 − ui )
=0
δ 21m1
( δ 22 m2 − ui )
Vẽ các biểu đồ mô men
uốn đơn vị Z1 = 1 và Z2 = 1.

8,97

0,78

3,90

(M1).1/13

6,24

Z2=1

2,34
9,68
(M2).1/13

0 ,356
)=
EI
0
δ 22 = ( M 2 )( M 2 ) = 0 ,4266
EI

0
δ 11 = ( M 1 )( M 1

δ 12 =

0
δ 21 = ( M 1 )( M 2

)=

0 ,12
EI

m2

m1

2EI

2EI
3m

Để xác định các chuyển
vị δ ik ta tạo các trạng thái
khả dĩ và vẽ các biểu đồ
mô men uốn đơn vị (M1o)
và (M2o) tương ứng trong
hệ cơ bản.
Áp dụng các nhân biểu đồ:

EI

EI

2m

2m

2m

2m

P=1

1

o
(M 1 )

P=1

1

o

(M 2)


1
ω1 =
= 65 ,69 1 / s
u1

ω2 =

1
= 99 ,1 1 / s
u2

Chu kỳ dao động:

m2

m1

2EI

2EI
3m

Thay các giá trị tìm
được vào phương trình tần

số ta thu được:

EI
2m

EI
2m

2m

2m

2π
2π
= 0 ,0956 ; T2 =
= 0 ,0634
T1 =
ω1
ω2

Với ω 1 = 65,69, ta có:

( δ 11m1 − u1 )y11 + δ 12 m2 y21 = 0

δ 21m1 y11 + ( δ 22 m2 − u1 )y21 = 0
Cho y11 = 1  y21 = - 0,6587


m2

m1

2EI

2EI
3m

Dạng chính thứ nhất của
dao động riêng và chuyển vị
tương ứng của các khối
lượng như hình vẽ.

y11( t ) = y11a1 sin( ω1t + ϕ1 )
y11( t ) = a1 sin( 65 ,69t + ϕ1 );
y21( t ) = y21a1 sin( ω1t + ϕ1 )
y21( t ) = -0,6587a1 sin( 65 ,69t + ϕ1 );

EI

EI
2m

2m

2m

2m

m2
m1

m1

y21 = 0,6587


Tương tự với ω 2 = 99,1 1/s:

δ 21m1 y12 + ( δ 22 m2 − u2 )y22 = 0

2EI

2EI
3m

( δ 11m1 − u2 )y12 + δ 12 m2 y22 = 0

m2

m1
EI

EI

2m

2m

2m

2m

Cho y12 = 1  y22 = 3,037
Dạng chính thứ hai của
dao động riêng và chuyển vị
tương ứng của các khối
lượng như hình vẽ.

y12( t ) = y12a2 sin( ω2t + ϕ2 )
y12( t ) = a1 sin( 99 ,10t + ϕ2 );
y22( t ) = y22a2 sin( ω2t + ϕ2 )
y22( t ) = 3,037a2 sin( 99 ,10t + ϕ2 );

m1
y12 = 1

y21 = 0,6587

m1
y12 = 1

y22 = 3,037


Phương trình dao động tổng quát của các khối lượng:

y1 ( t ) = y11 a1 sin( ω1t + ϕ1 ) + y12 a2 sin( ω 2 t + ϕ 2 ) =
= a1 sin( 65 ,69.t + ϕ1 ) + a2 sin( 99 ,1.t + ϕ 2 )
y2 ( t ) = y21 a1 sin( ω1t + ϕ1 ) + y22 a2 sin( ω 2 t + ϕ 2 ) =
= -0,6587 a1 sin( 65 ,69.t + ϕ1 ) + 3,037 a2 sin( 99 ,1.t + ϕ 2 )
Các đại lượng a1, ϕ 1, a2, ϕ 2, được xác định theo

điều kiện ban đầu của dao động ở thời điểm t = 0.

y1 ( t ) = y1 ( 0 ); y 1 ( t ) = v1 ( 0 );
y2 ( t ) = y2 ( 0 ); y 2 ( t ) = v2( 0 )


* Cách sử dụng tính đối xứng của hệ trong dao động:
Đối với những hệ đối xứng mang các khối lượng có
giá trị và vị trí được bố trí đối xứng, hệ sẽ dao động
tương ứng với hai loại dao động chính sau:
•
Dạng dao động đối xứng tương ứng với các lực
quán tính tác dụng đối xứng.
•
Dạng dao động phản xứng tương ứng với các lực
quán tính tác dụng phản xứng


1) Biện pháp biến đổi hệ về sơ đồ nửa hệ tương đương:
m1

m2

l

m1

m1

l

l

m2 / 2

Đối xứng
m1
m3

m2

m1

m4

m6

m4
m5

l

l

m2//2

m1
m3

m3

m6

m6

m4

m1
l

l

Phản xứng
m2//2

m1
m3

m4

m6

l

m2 / 2

m5

l

2) Biện pháp sử dụng chuyển vị kép:
• Dạng dao động đối xứng:
Biểu thị chuyển vị tại
cặp khối lượng có vị trí
đối xứng theo chuyển vị
kép:

Yk (t ) = 2 yk (t ), k = 1,2...( n − 1)
Khối lượng mn không
có chuyển vị kép. yn(t).
Các cặp khối lượng
đối xứng phát sinh các
cặp lực quán tính đối
xứng

y1

yk

m1

mk

yn-1

mn-1

yn

yn-1

mn

mn-1

l

y1

m1

mk
l

Z1=1

Z1=1

Zk=1

yk

(M 1 )
Zk=1

(M k)
Zn=1

(M n)

Gọi δ ik (i ≠ k, k ≠ n):
chuyển vị đơn vị ứng
với vị trí và phương của
các cặp lực quán tính
đối xứng Zi(t) do cặp lực
đối xứng Zk = 1 tác dụng
tĩnh tạ mk.
Gọi δ im : chuyển vị đơn
vị ứng với vị trí và
phương của các cặp lực
quán tính đối xứng Zi(t)
do lực Zn = 1 tác dụng
tĩnh tạ vị trí mn .

y1

yk

yn-1

yn

yn-1

l

y1

l

Z1=1

Z1=1

Zk=1

yk

(M 1 )
Zk=1

(M k)
Zn=1

(M n)