TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP HÀ NỘI HANOI UNIVERSITY OF INDUSTRY Tập san SINH VIÊN NGHIÊN CỨU KHOA HỌC Số 14 2024 291
NGHIÊN CỨU ẢNH HƯỞNG CỦA BIÊN DẠNG MẶT ĐƯỜNG
ĐẾN DAO ĐỘNG Ô TÔ
THE STUDY OF THE INFLUENCE OF ROAD SURFACE PROFILE
ON VEHICLE VIBRATION
Trần Huy Hoàng1,*, Nghiêm Ngọc Cường1, Nguyễn Văn Duy1,
Nguyễn Minh Đức1, Nguyễn Mạnh Trường1, Lê Duy Long2
1Lớp KTOT 02- K16, Trường Cơ khí - Ô tô, Trường Đại học Công nghiệp Hà Nội
2Trường Cơ khí - Ô tô, Trường Đại học Công nghiệp Hà Nội
*Email: hoanghuytrank7pl@gmail.com
TÓM TẮT
Đề tài xây dựng mô hình cấu trúc không gian xe ô tô con khảo sát bao gồm các thành phần trong hệ thống treo. Từ mô
hình đã xây dựng được tiến hành mô phỏng mấp mô mặt đường từ đó phân tích ảnh hưởng của các thông số này tới độ êm
dịu khi chuyển động của xe ô tô, đánh giá thời gian dập tắt dao động của hệ thống treo phỏng biên dạng mấp
mặt đường thực tế. Việc đưa biên dạng đường thực tế vào trong phương trình dao động giúp chúng ta được kết quả
chính xác hơn trong quá trình nghiên cứu. Kết quả đạt được sẽ được đem so sánh với kết quả thực nghiệm đo thực tế từ đó
làm cơ sở thông số đầu vào cho bài toán dao động ô tô.
Từ khóa: Biên dạng mặt đường, Mô phỏng, Matlab-Simulink, Dao động ô tô.
ABSTRACT
The project aims to construct a structural model of a passenger car's spatial configuration, including components of the
suspension system. From the constructed model, the road surface roughness is simulated to analyze the impact of these
parameters on the ride comfort of the car, evaluate the damping time of the suspension system, and simulate the actual
road surface roughness profile. Incorporating the real road profile into the oscillation equations allows us to achieve more
accurate results in the research process. The obtained results will be compared with experimental measurements from real-
world data, thereby serving as the input parameters for the automotive oscillation problem.
Keywords: Road profile, Simulation, Matlab-Simulink, Vehicle dynamics.
1. MỞ ĐẦU
Ô loại phương tiện vận tải vai trò hết sức quan
trọng trong xã hội hiện nay. Hiện nay, ô vẫn đang được
sử dụng rộng rãi trong tất cả các lĩnh vực kinh tế đời sống
xã hội. Cùng với sự phát triển mạnh mẽ của khoa học công
nghệ, công nghiệp ô trên thế giới đặc biệt tại Việt Nam
đã và đang những bước phát triển từng ngày với yêu cầu
ngày càng cao, ngành ô hiện nay đã cho ra đời nhiều mẫu
cải tiến để phục vụ cho các nhu cầu vận chuyển. Việc
nghiên cứu đánh giá các yếu tố ảnh hưởng tới quá trình
dao động ô giúp chúng ta thêm thuyết để từ đó phát
triển đề tài. Khi ô tô chuyển động, ảnh hưởng của dao động
là một trong các yếu tố quan trọng ảnh hưởng tới quá trình
vận hành, đặc biệt là tính êm dịu của ô tô. Một trong những
hệ thống liên quan trực tiếp tới hệ thống dao động của ô tô
là hệ thống treo, hệ thống treo trong quá trình xe di chuyển
phải chịu những tác động từ mấp mô mặt đường do đó việc
đánh giá ảnh hưởng của mấp mặt đường tới dao động
của ô tô là điều cần thiết nhằm đảm bảo độ êm dịu của ô tô.
Việc xây dựng mô hình dao động, hình toán
phỏng hệ thống treo trên ô tô cần thiết để khảo sát đánh
giá chất lượng dao động của ô tô. Do đó, đề tài nghiên cứu
khoa học sẽ định hướng nghiên cứu và đánh giá ảnh hưởng
của mấp mặt đường tới dao động ô tô. Trên những
sở đó, nhóm tác giả quyết định nghiên cứu ảnh hưởng của
biên dạng mặt đường đến dao động ô tô.
Nhóm nghiên cứu ứng dụng phần mềm Matlab-
Simulink trong việc thiết kế, phỏng không gian của xe
khảo sát. Các kết quả nghiên cứu này có thể được coi là cơ
sở quan trọng cần thiết cho c nghiên cứu tiếp theo
trong lĩnh vực dao động ô tô.
2. NỘI DUNG NGHIÊN CỨU
Xây dựng hình cấu trúc không gian xe ô tô con khảo
sát bao gồm các thành phần trong hệ thống treo. Từ hình
TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP HÀ NỘI HANOI UNIVERSITY OF INDUSTRY Tập san SINH VIÊN NGHIÊN CỨU KHOA HỌC Số 14 2024 292của xe khảo sát, nhóm tác giả đã xây dựng được tiến hành
phỏng mấp mặt đường từ đó phân tích ảnh hưởng
của các thông số này tới độ êm dịu khi chuyển động của xe
ô đánh giá thời gian dập tắt dao động của hệ thống
treo. Với nguyên lý đó, nhóm nghiên cứu đã mô phỏng dựa
trên các biên dạng mặt đường thực tế gồm: biên dạng mặt
đường hình sin, biên dạng mặt đường hình xung, biên dạng
mặt đường ngẫu nhiên theo tiêu chuẩn ISO 8608 khảo
sát cho xe chạy qua trong các dải vận tốc khác nhau để thu
thập đồ thị khảo sát dao động xe ô tô [4].
2.1. Các dạng mấp mô mặt đường
Chuyển động của ô trên bề mặt đường không bằng
phẳng sẽ phát sinh các dao động của các khối lượng phần
treo và khối lượng phần không được treo của ô tô. Độ mấp
của bề mặt đường nguồn kích thích chính cho ô tô
dao động. Khi nghiên cứu hình dao động của ô cần
thiết phải tả toán học biến dạng bề mặt đường sẽ tham
gia vào phải của hệ phương trình vi phân tả chuyển động
dao động của hệ Điều kiện đường trong thực tế sử dụng ô
tô rất đa dạng. Ảnh hưởng của chúng tới dao động của ô tô
được xác định bởi kích thước hình học, hình dạng đặc
tính thay đổi của chúng. Tuỳ theo chiều dài của mấp mô
hoặc chiều cao của nó thể phân ra c nhóm đặc trưng
khác nhau của biên dạng bề mặt đường, thể phân thành
ba loại chính như sau: loại 1: mấp chiều dài ngắn, c
động của chúng lên các va đập (tác động xung). loại 2: mấp
dạng hàm điều hoà (hình sin). loại 3: mấp ngẫu
nhiên.
Hình 1. Các dạng mấp mô mặt đường loại 1
Trong nghiên cứu dao động của ô rất nhiều
hình dao động của nhiều tác giả khác nhau, y thuộc vào
mục tiêu nghiên cứu. Trong nghiên cứu này, nhóm tác gi
chỉ đưa ra 3 loại hình đặc trưng với mục đích một
cách nhìn tổng quan về phương pháp xây dựng hình
gồm mô hình dao động ¼, hình dao động ½ và mô hình
không gian.
2.2. Cơ sở lý thuyết chung về dao động
Dao động được tả theo khái niệm sau: dao động
một quá trình trong đó một đại lượng vật (như hóa học,
sinh học,...) thay đổi theo thời gian và có một đặc điểm nào
đó lặp lại ít nhất một lần. Ô tô được xem như là một hệ dao
động cơ học, gồm nhiều khối lượng được liên kết với nhau
tương tác chặt chẽ với mặt đường dạng phức tạp.
Trong thuyết ô tô, khi nghiên cứu về dao động của nó,
thường giả định rằng khối lượng của ô được tập trung,
được chia thành hai phần: phần được treo phần không
được treo. Dao động của phần được treo phức tạp
thể được xem như sự kết hợp của một số dao động đơn
giản. Trong trường hợp y, phần được treo thể thực hiện
đồng thời ba loại di chuyển tịnh tiến theo các trục tọa độ O,
O, O ba loại di chuyển góc xung quanh các trục này.
Trong đó, di chuyển tịnh tiến theo trục O di chuyển góc
xung quanh trục O (gọi lắc dọc), O (gọi lắc ngang)
thường có ảnh hưởng lớn nhất đến sêm ái và an toàn trong
việc di chuyển của ô tô.
Các nguồn gây ra dao động gồm: mấp mặt đường,
điều kiện tự nhiên, lực quán tính,… Chuyển động của ô
trên mặt đường không bằng phẳng sẽ sinh ra các dao động
của khối lượng phần không được treo khối lượng phần
được treo của ô tô. Độ mấp của mặt đường nguồn
kích thích chính cho ô tô dao động.
2.3. Xây dựng mô hình cấu trúc ô tô
Nhóm tác giả đã chọn xe Toyota Vios 2016 để tiến hành
xây dựng mô hình cấu trúc không gian:
Hình 2. Hình ảnh xe Toyota Vios 2016
Ô con khối lượng kích thước nhỏ. Ô được
thiết kế khối lượng phân bố đối xứng qua mặt phẳng dọc.
Thân xe được phân thành phần được treo trước sau,
tương ứng các cầu xe trước sau. Các githiết trong
quá trình mô tả động lực học được sử dụng: mô hình động
lực học không gian, mô hình ô tô con có trục đối xứng theo
trục dọc của xe, thân xe coi như một tấm phẳng có khối
lượng ma đặt tại trọng tâm thân xe. Trong bài toán xét dao
động của xe khi chuyển động thẳng, thân xe ba chuyển
động bao gồm: chuyển động tịnh tiến theo trục thẳng đứng
z hai góc xoay ơng ứng trục x y, bốn bánh xe với
các thông số độ cứng và hệ số giảm chấn. Từ qtrình phân
tích ta xây dựng hình không gian cho xe để tiến hành
mô phỏng:
TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP HÀ NỘI HANOI UNIVERSITY OF INDUSTRY Tập san SINH VIÊN NGHIÊN CỨU KHOA HỌC Số 14 2024 293
Hình 3. Mô hình không gian xe con tiến hành mô phỏng
Từ mô hình đã xây dựng áp dụng định luật hai Newton,
nhóm tác giả xây dựng bảy phương trình toán học tả dao
động của xe gồm thân xe, góc lắc ngang thân xe, góc lắc
dọc thân xe và bốn bánh xe như sau:
- Phương trình:
ξ
=−  +
= − .ξ + +
. ξ + +
+  .( ξ)
= −ξ+  ξ(+  )
+   + 
- Phương trình:
ξ
=−  +
= − .ξ +
. ξ +
+  .( ξ)
= −ξ+ + ξ(+  )
+ + + 
- Phương trình:
ξ
=−  +
= − .ξ
. ξ 
+  .( ξ)
= −ξ+ ++ ξ(+  )
+ +++ 
- Phương trình:
ξ
=−  +
= − .ξ +
. ξ +
+  .( ξ)
= −ξ+ + ξ(+  )
+ ++ 
- Phương trình:
= + + +
+  + + +
=.ξ+ +
+ .ξ ++
+ .ξ +
+ . ξ+
+ .ξ
+ .ξ
+ .ξ+
+ .ξ+
= −.(+ + + )+ ξ+ξ+ ξ
+ξ+ (+)
+(+)
(+++)+ξ+ξ
+ξ+ξ
+(+)+(
+)
- Phương trình
= −  +
+ + +  .
= −.ξ++
ξ++
+. ξ+
+ξ+
+. ξ
+ξ
. ξ+
ξ+.
=(+)ξ+ξ+ξ
ξ(+++)
(+)
+(+)ξ
+ξ+ξξ
(+++)(
+)
- Phương trình:
= ( + + +).1
+( + + +).2
=−.ξ++
+ξ++
+. ξ+
+ξ+.
+.ξ
+ξ
+. ξ+
+ξ+.
TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP HÀ NỘI HANOI UNIVERSITY OF INDUSTRY Tập san SINH VIÊN NGHIÊN CỨU KHOA HỌC Số 14 2024 294=(+)ξξ
+ξ+ξ
(+)
(
+
+
+
)
+(+)ξ
ξ+ξ+ξ
(+)
(
+
+
+
)
2.4. Xây dựng mô hình xe trong Simulink
Nhóm tác giả xây dựng phương trình trạng thái của xe:
=+
= + 
Trong đó:
= [   0 0 0] vetor
đầu vào
=
là vetor trạng thái
=
là vetor đầu ra
Với [s] = 
Từ phương trình trạng thái, nhóm nghiên cứu tiến hành
xây dựng mô hình mô phỏng xe trên Simulink:
Hình 4. Mô hình xe trên Simulink
Mô hình xe trong Simulink có thể được mô tả một cách
ngắn gọn như sau: khối mặt đường sẽ tạo ra các tín hiệu
mấp thể hiện cho biên dạng mặt đường. Khối lực kích
thích mặt đường sẽ nhận tín hiệu mấp được tạo ra t
khối mặt đường, tác dụng lực kích thích n bốn bánh xe với
bảy tín hiệu đầu vào được khối Mux gộp lại thành một tín
hiệu để làm tín hiệu đầu vào cho phương trình trạng thái.
Khối kết quả sẽ nhận hai mươi mốt tín hiệu đầu ra ứng với
các giá trcủa sự dịch chuyển, vận tốc gia tốc của thân
xe và bốn bánh bánh xe.
3. KẾT QUẢ
Sau khi nhóm tác giả tiến hành mô phỏng và đưa ra kết
quả dịch chuyển của bốn bánh xe, thân xe và các góc xoay
thể thấy dạng mặt đường hình sin đặc biệt nhất bởi
dạng dịch chuyển thẳng đứng của thân xe tại trọng tâm
một hàm điều hòa. Sau khi qua mấp mặt đường, thân xe
vẫn dịch chuyển với một dạng biên độ điều hòa.
Dưới tác động của biên dạng mặt đường hình sin, nhóm
tác giả đã phỏng, đưa ra kết quả như hình bên dưới
nhận thấy các giá trị này thay đổi trong một khoảng thời
gian ban đầu sau đó giá trịc gia tốc này trở lại giá trị điều
hòa. Giá trgia tốc tại các vị trí tiếp xúc của thân xe thay
đổi trong giai đoạn ban đầu khi tín hiệu mặt đường hình sin
tác động vào xe sau thời gian đó các giá trị gia tốc cũng trở
lại giá trị điều hòa.
Hình 5. Gia tốc của thân xe khi xe đi qua mấp mô hình sin
Hình 6. Gia tốc thân xe khi xe đi qua mấp mô hình xung
Hình 7. Gia tốc thân xe khi xe đi qua mấp mô ngẫu nhiên
TRƯỜNG ĐẠI HỌC CÔNG NGHIỆP HÀ NỘI HANOI UNIVERSITY OF INDUSTRY Tập san SINH VIÊN NGHIÊN CỨU KHOA HỌC Số 14 2024 295Để xác định mức độ dao động của xe, sử dụng đồ thị gia
tốc của THÂN XE để khảo sát. Qua ba loại mấp mặt
đường, để giảm độ cao của thân xe khi va chạm với mấp
mô, ta cần di chuyển qua mấp với tốc độ 20 km/h để
đảm bảo thân xe ổn định, nếu di chuyển qua những mấp mô,
đặc biệt là hình xung với tốc độ lớn hơn 20 km/h hệ thống
treo của xe Toyota Vios 2016 tiến hành khảo sát sẽ không
thể dập tắt dao động trong khoảng thời gian ngắn hơn 0,6
giây độ dịch chuyển biên độ cao hơn sẽ làm xe mất
khả năng cân bằng, gây ảnh hưởng tới hệ thống treo và độ
an toàn của người lái xe.
4. KẾT LUẬN
Từ nghiên cứu, mô phỏng các hệ phương trình vi phân,
khảo sát dao động của ô Toyota Vios 2016 trên miền thời
gian cho ba dạng mấp mô mặt đường là mấp mô biên dạng
hình sin, mấp biên dạng hình xung mấp biên dạng
ngẫu nhiên. Đưa ra được tốc độ tối ưu để hệ thống treo của
xe khảo sát có thể dập tắt dao động một cách tốt nhất, đảm
bảo độ êm dịu khi xe di chuyển. Các kết quả khảo sát dao
động trên miền thời gian đối với ô tô Toyota Vios 2016 khi
đi qua các mấp đơn dạng sin, hình xung đi trên đường
dạng hàm điều hòa cho thấy với các thông số của hệ thống
treo hiện trên xe đảm bảo được các chỉ tiêu về êm dịu
chuyển động đối với xe con.
TÀI LIỆU THAM KHẢO
[1]. Võ Văn Hường, Nguyễn Tiến Dũng, Dương Ngọc Khánh, Đàm Hoàng Phúc, 2014. Động lực học ô tô. NXB Giáo
dục Việt Nam.
[2]. Nguyễn Phùng Quang. 2005. Matlab & Simulink. NXB Khoa học và Kỹ thuật, Hà Nội.
[3]. Đặng Việt Hà, 2010. Nghiên cứu ảnh hưởng của một số thông số đến độ êm dịu chuyển động của ô tô khách được
đóng mới tại Việt Nam, luận án tiến sỹ, Đại học Giao Thông Vận Tải.
[4]. ISO 5982:2001, Mechanical vibration and shock - Range of idealized values to characterize seated-body
biodynamic response under vertical vibration.