BÀI GIẢNG MATLAB - Trang 200

198

Trong ví dụ này ta cần chú ý ñến các lực, ngẫu lực tập trung, gối mềm và gối cứng:

Lực tập trung tác dụng tại vị trí

iF

suy ra

(

)

/

q iF

F h

=

.

Mômen tập trung tại vị trí

iM

ñược thay bằng hai lực ngược chiều cùng trị số tác

dụng tại hai nút (trước và sau)

2

2

(

1)

/ 2 ,

(

1)

/ 2

q iM

M

h

q iM

M

h

= −

+

=

Gối cứng

(

)

0

i iGc

w

=

=

, tức là dòng thứ i của ma trận A chỉ có phần tử

( , )

1,

a i i

=

( , )

0,

,

a i j

i

j

=

( )

0

b i

=

.

Gối mềm, coi lực lò xo,

lx

i

F

cw

= −

, như là lực tập trung tác dụng tại nút

( / )

iGm

iGm

q

c h w

= −

,

Như thế phương trình với chỉ số

iGm

ñược viết lại thành

4

2

1

1

2

4

(6

)

4

0

i

i

i

i

i

h c

w

w

w

w

w

EI h

+

+

+

+

=

(gối mềm tại ñây)

Thực hiện trong Matab với số ñiểm chia

3001

n

=

, số ñoạn chia là 3000, thu ñược

kết quả như trên hình 9-11.

Nếu cho hệ số cứng lò xo tăng lên rất lớn (

15

10

c

=

N/m) thì ta có thể coi như gối

cứng và thu ñược kết quả như trên hình 9-12.

0

1

2

3

4

5

- 5

0

5

x 10

-4

Do thi do vong cua dam

w

(x

)

[m

]

0

1

2

3

4

5

- 100

0

100

Bieu do momen

M

[

N

m

]

0

1

2

3

4

5

- 200

0

200

Bieu do luc cat

x [m]

Q

[

N

]

Hình 9-11. Kết quả tính toán (xem Beam_SPHH_Vd3.m)