Phương pháp MPMM trong phân tích động lực học kết cấu tấm mindlin chịu tải trọng điều hòa di động

ĐẠI HỌC QUỐC GIA TP. HỒ CHÍ MINH
TRƯỜNG ĐẠI HỌC BÁCH KHOA
--------------------

TRẦN TRUNG HIẾU

PHƯƠNG PHÁP MPMM TRONG PHÂN TÍCH ĐỘNG LỰC
HỌC KẾT CẤU TẤM MINDLIN CHỊU TẢI TRỌNG ĐIỀU
HÒA DI ĐỘNG
Chun ngành : Kỹ thuật xây dựng cơng trình dân dụng và công nghiệp
Mã số ngành

: 60580208

LUẬN VĂN THẠC SĨ

TP. HỒ CHÍ MINH, tháng 01 năm 2019


CƠNG TRÌNH ĐƯỢC HỒN THÀNH TẠI
TRƯỜNG ĐẠI HỌC BÁCH KHOA
ĐẠI HỌC QUỐC GIA TP. HỒ CHÍ MINH

Cán bộ hướng dẫn khoa học: PGS. TS. Lương Văn Hải

Cán bộ chấm nhận xét 1: PGS. TS. Nguyễn Trung Kiên

Cán bộ chấm nhận xét 2: TS. Nguyễn Tấn Cường

Luận văn thạc sĩ được bảo vệ tại Trường Đại học Bách Khoa, ĐHQG Tp. HCM
ngày 13 tháng 01 năm 2019.

Thành phần Hội đồng đánh giá Luận văn thạc sĩ gồm:
1. PGS. TS. Bùi Công Thành

Chủ tịch hội đồng

2. TS. Nguyễn Hồng Ân

Thư ký

3. PGS. TS. Nguyễn Trung Kiên

Thành viên. Phản biện 1

4. TS. Nguyễn Tấn Cường

Thành viên. Phản biện 2

5. PGS. TS. Nguyễn Trọng Phước

Thành viên

Xác nhận của Chủ tịch Hội đồng đánh giá LV và Trưởng Khoa quản lý chuyên
ngành sau khi Luận văn đã được sửa chữa (nếu có).
CHỦ TỊCH HỘI ĐỒNG

TRƯỞNG KHOA
KỸ THUẬT XÂY DỰNG


i

ĐẠI HỌC QUỐC GIA TP.HCM

CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM

TRƯỜNG ĐẠI HỌC BÁCH KHOA

Độc Lập - Tự Do - Hạnh Phúc

NHIỆM VỤ LUẬN VĂN THẠC SĨ
Họ và tên học viên: TRẦN TRUNG HIẾU

MSHV: 7141170

Ngày, tháng, năm sinh: 03/02/1991

Nơi sinh: Phú Yên

Chuyên ngành: Kỹ thuật xây dựng CT DD & CN

Mã số: 60580208

I. TÊN ĐỀ TÀI: Phương pháp MPMM trong phân tích động lực học kết cấu
tấm Mindlin chịu tải trọng điều hòa di động
II. NHIỆM VỤ VÀ NỘI DUNG
1. Thiết lập phương trình năng lượng của tấm Mindlin, mơ hình nền đường nhiều
lớp và thiết lập cơng thức ma trận kết cấu tấm Mindlin trên nền nhiều lớp sử dụng
phương pháp phần tử nhiều lớp tấm chuyển động MPMM
2. Phát triển thuật tốn, sử dụng ngơn ngữ lập trình Matlab để giải hệ phương trình
động của bài tốn.

3. Kiểm tra độ tin cậy của chương trình tính bằng cách so sánh kết quả của Luận
văn với kết quả của các phương pháp khác.
4. Tiến hành thực hiện các ví dụ số nhằm khảo sát ảnh hưởng của các nhân tố quan
trọng đến ứng xử động của kết cấu tấm, từ đó rút ra kết luận và kiến nghị.
5. NGÀY GIAO NHIỆM VỤ

: 03/08/2018

6. NGÀY HOÀN THÀNH NHIỆM VỤ : 03/12/2018
7. HỌ VÀ TÊN CÁN BỘ HƯỚNG DẪN: PGS. TS. Lương Văn Hải
Tp. HCM, ngày 03 tháng 12 năm 2018

CÁN BỘ HƯỚNG DẪN

CHỦ TỊCH HỘI ĐỒNG NGÀNH

PGS. TS. Lương Văn Hải
TRƯỞNG KHOA KỸ THUẬT XÂY DỰNG


ii

LỜI CẢM ƠN
Luận văn thạc sĩ Kỹ thuật xây dựng cơng trình dân dụng và cơng nghiệp nằm
trong hệ thống bài luận cuối khóa nhằm trang bị cho Học viên cao học khả năng tự
nghiên cứu, biết cách giải quyết những vấn đề cụ thể đặt ra trong thực tế xây dựng…
Đó là trách nhiệm và niềm tự hào của mỗi học viên cao học.
Để hoàn thành Luận văn này, ngoài sự cố gắng và nỗ lực của bản thân, tôi đã
nhận được sự giúp đỡ từ nhiều tập thể và các cá nhân. Tôi xin ghi nhận và tỏ lòng
biết ơn đến các tập thể và các cá nhân đã dành cho tơi sự giúp đỡ q báu đó.

Đầu tiên tơi xin bày tỏ lịng biết ơn sâu sắc đến Thầy PGS. TS. Lương Văn
Hải. Thầy đã đưa ra gợi ý đầu tiên để hình thành nên ý tưởng của đề tài và Thầy góp
ý cho tơi rất nhiều về cách nhận định đúng đắn trong những vấn đề nghiên cứu, cũng
như cách tiếp cận nghiên cứu hiệu quả.
Tôi xin chân thành cảm ơn quý Thầy Cô Khoa Kỹ Thuật Xây dựng, trường
Đại học Bách Khoa TP. HCM đã truyền dạy những kiến thức q giá cho tơi, đó cũng
là những kiến thức không thể thiếu trên con đường nghiên cứu khoa học và sự nghiệp
của tôi sau này.
Tôi cũng xin gửi lời cảm ơn đến NCS. Cao Tấn Ngọc Thân, NCS. Nguyễn
Xuân Vũ đã giúp đỡ tôi rất nhiều trong quá trình thực hiện Luận văn này.
Luận văn thạc sĩ đã hoàn thành trong thời gian quy định với sự nỗ lực của bản
thân, tuy nhiên không thể khơng có những thiếu sót. Kính mong q Thầy Cơ chỉ dẫn
thêm để tôi bổ sung những kiến thức và hồn thiện bản thân mình hơn.
Xin trân trọng cảm ơn.

Tp. HCM, ngày 03 tháng 12 năm 2018

Trần Trung Hiếu


iii

TÓM TẮT LUẬN VĂN THẠC SĨ
Nhiều nghiên cứu đã được thực hiện gần đây để dự đoán chuyển vị và ứng suất
của nền đường cao tốc, nền đường băng... Trong các nghiên cứu đó, tải trọng tác dụng
bởi xe chạy trên đường cao tốc, máy bay chuyển động trên đường băng thường được
mơ hình là tải trọng hằng số. Trong thực tế, khi xe cộ chạy trên đường, tải trọng tác
dụng lên bề mặt đường sẽ có cường độ thay đổi theo thời gian do độ nhám của bề mặt
đường và hệ thống động cơ của xe. Do đó, việc mơ phỏng xe chạy trên đường là tải
trọng điều hịa di động sẽ phản ánh đúng với thực tế hơn so với khi xem xe cộ là tải

trọng hằng số di động. Và kết cấu đường băng, nền đường thường được cấu tạo nhiều
lớp bao gồm: lớp bê tông, lớp nhựa đường, lớp xi măng đá, đặt trên nền đất. Trong
Luận văn này, phương pháp nhiều lớp tấm chuyển động MPMM (Multi-Layer Plate
Moving Method) được đề xuất để giải quyết bài toán tấm dày trên nền nhiều lớp chịu
tải trọng điều hòa di động. Phương pháp MPMM được đề xuất dựa trên phương pháp
phần tử chuyển động MEM (Moving Element Method). Các phương trình chuyển
động của tấm, các ma trận kết cấu được thiết lập trong hệ tọa độ tương đối chuyển
động cùng vận tốc với vận tốc của lực di chuyển. Do đó, phương pháp này sẽ tránh
được việc cập nhật véctơ tải trọng tương ứng với mơ hình tấm. Ảnh hưởng của các
thông số đến ứng xử động của tấm Mindlin trên nền nhiều lớp như hệ số độ cứng nền,
hệ số độ cản nền, hệ số độ cứng và hệ số độ cản lớp liên kết giữa hai tấm, cũng như
module đàn hồi, chiều dày của hai lớp tấm và ảnh hưởng của tần số tải trọng điều hòa
di động, sự lệch pha và khoảng cách giữa hai tải trọng điều hòa di động cùng tác dụng
lên tấm được khảo sát. Các kết quả nghiên cứu trong Luận văn hy vọng có thể là một
trong những tài liệu tham khảo hữu ích trong việc nghiên cứu tấm Mindlin trên nền
nhiều lớp chịu tải trọng điều hòa di động.


iv

LỜI CAM ĐOAN
Tôi xin cam đoan đây là công việc do chính tơi thực hiện dưới sự hướng dẫn
của Thầy PGS. TS. Lương Văn Hải.
Các kết quả trong Luận văn là đúng sự thật và chưa được công bố ở các nghiên
cứu khác.
Tôi xin chịu trách nhiệm về công việc thực hiện của mình.

Tp. HCM, ngày 03 tháng 12 năm 2018

Trần Trung Hiếu

v

MỤC LỤC
NHIỆM VỤ LUẬN VĂN THẠC SĨ ......................................................................... i
LỜI CẢM ƠN ......................................................................................................... ii
TÓM TẮT LUẬN VĂN THẠC SĨ ......................................................................... iii
LỜI CAM ĐOAN ................................................................................................... iv
MỤC LỤC .............................................................................................................. v
DANH MỤC CÁC HÌNH VẼ .............................................................................. viii
DANH MỤC CÁC BẢNG BIỂU .......................................................................... xii
MỘT SỐ KÝ HIỆU VIẾT TẮT ........................................................................... xiv
CHƯƠNG 1 TỔNG QUAN ................................................................................... 1
1.1. Giới thiệu.................................................................................................... 1
1.2. Tình hình nghiên cứu và tính cấp thiết của đề tài ........................................ 2
1.2.1.

Các cơng trình nghiên cứu ngồi nước ........................................... 3

1.2.2.

Các cơng trình nghiên cứu trong nước ............................................ 5

1.3. Mục tiêu và hướng nghiên cứu .................................................................... 6
1.4. Cấu trúc của Luận văn ................................................................................ 7
CHƯƠNG 2 CƠ SỞ LÝ THUYẾT ......................................................................... 8
2.1. Lý thuyết tấm Mindlin ................................................................................ 8
2.1.1.

Giới thiệu tổng quát........................................................................ 8

2.1.2.

Lý thuyết tấm có kể đến biến dạng trượt của Reissner-Mindlin ...... 9

2.1.3.

Biến dạng của tấm và mối quan hệ giữa biến dạng – chuyển vị .... 11

2.1.4.

Biến dạng của tấm và mối quan hệ giữa ứng suất – biến dạng ...... 12

2.1.5.

Phương trình năng lượng của tấm ................................................. 14

2.2. Phần tử đẳng tham số ................................................................................ 15
2.2.1.

Khái niệm phần tử đẳng tham số .................................................. 15

2.2.2.

Phép tích phân số - Phép cầu phương Gauss ................................. 18

2.3. Mô hình nền đường nhiều lớp ................................................................... 19
2.3.1.

Hệ số độ cứng đàn hồi K s ........................................................... 21


vi

2.3.2.

Hệ số cản Cs ................................................................................ 22

2.4. Phương pháp phần tử nhiều lớp tấm chuyển động MPMM........................ 23
2.4.1.

Thiết lập công thức ma trận kết cấu tấm Mindlin trên nền nhiều
lớp sử dụng phương pháp phần tử nhiều lớp tấm chuyển động
MPMM ........................................................................................ 24

2.4.2.

Tải trọng điều hòa di động............................................................ 34

2.4.3.

Tần số tới hạn và vận tốc tới hạn trong bài toán tấm Mindlin chịu
tác dụng của tải trọng điều hòa di động ........................................ 35

2.5. Phương pháp Newmark............................................................................. 36
2.6. Thuật toán sử dụng trong luận văn ............................................................ 38
2.6.1.

Thông số đầu vào ......................................................................... 38

2.6.2.

Giải bài toán theo dạng chuyển vị ................................................ 40

2.6.3.

Giải bài toán theo dạng gia tốc ..................................................... 40

2.6.4.

Độ ổn định và hội tụ theo phương pháp Newmark........................ 40

2.7. Lưu đồ tính tốn ....................................................................................... 41
CHƯƠNG 3

KẾT QUẢ PHÂN TÍCH SỐ .......................................................... 44

3.1. Kiểm chứng chương trình Matlab ............................................................. 47
3.1.1.

Bài tốn 1a: Phân tích ứng xử của tấm Mindlin trên nền nhiều
lớp khi chịu tác dụng của tải trọng tĩnh khi xem tấm xi măng đá
và nền là cứng vơ cùng ................................................................. 47

3.1.2.

Bài tốn 1b: Phân tích ứng xử của tấm Mindlin trên nền nhiều
lớp chịu tác dụng của tải trọng hằng số di động ............................ 50

3.2. Phân tích động lực học tấm Mindlin trên nền nhiều lớp chịu tác dụng của
tải trọng điều hịa di động ......................................................................... 52
3.2.1.

Bài tốn 2: Khảo sát sự hội tụ của bài toán ................................... 52

3.2.2.

Bài toán 3: Khảo sát ảnh hưởng của tỉ số độ cứng giữa nền và
lớp liên kết giữa hai tấm đối với chuyển vị tấm ............................ 56


vii

3.2.3.

Bài toán 4: Khảo sát ảnh hưởng của tỉ số độ cản giữa nền và lớp
liên kết giữa hai tấm đối với chuyển vị tấm .................................. 64

3.2.4.

Bài toán 5: Khảo sát ảnh hưởng của tỉ số module đàn hồi giữa
hai tấm đối với chuyển vị tấm ...................................................... 70

3.2.5.

Bài toán 6: Khảo sát ảnh hưởng của tỉ số chiều dày giữa hai tấm
đối với chuyển vị tấm ................................................................... 76

3.2.6.

Bài toán 7: Khảo sát ảnh hưởng của biên độ của tải trọng điều hòa
di động Po đối với chuyển vị tấm ................................................ 82

3.2.7.

Bài toán 8: Khảo sát ảnh hưởng của vận tốc V của tải trọng điều
hòa di động đối với chuyển vị tấm................................................ 84

3.2.8.

Bài toán 9: Khảo sát ảnh hưởng của tần số của tải trọng điều hòa
di động ωo đối với chuyển vị tấm ............................................... 88

3.2.9.

Bài toán 10: Khảo sát ảnh hưởng của độ lệch pha  của 02 tải
trọng điều hòa di động cùng tác dụng lên tấm đối với chuyển vị
tấm ............................................................................................... 90

3.2.10. Bài toán 11: Khảo sát ảnh hưởng của khoảng cách a giữa 02 tải
trọng điều hòa di động cùng tác dụng lên tấm đối với chuyển vị
tấm ............................................................................................... 93
CHƯƠNG 4

KẾT LUẬN VÀ KIẾN NGHỊ ....................................................... 98

4.1. Kết luận .................................................................................................... 98
4.2. Kiến nghị ................................................................................................ 101
TÀI LIỆU THAM KHẢO ................................................................................... 102

LÝ LỊCH TRÍCH NGANG ................................................................................. 108


viii

DANH MỤC CÁC HÌNH VẼ
Hình 1.1. Mặt cắt nền đường băng nhiều lớp ......................................................... 1
Hình 1.2. Mơ hình tải trọng cố định và phần tử tấm chuyển động (MPMM) .......... 2
Hình 2.1. Mơ hình động học của kết cấu tấm theo lý thuyết Kirchhoff .................. 9
Hình 2.2. Mơ hình động học của kết cấu tấm theo lý thuyết Mindlin ................... 11
Hình 2.3. Quy ước chiều dương của chuyển vị  và hai chuyển vị xoay  x ,  y
của tấm Mindlin ................................................................................. 11
Hình 2.4. Phần tử tứ giác Q9 trong hệ tọa độ địa phương .................................... 16
Hình 2.5. Phần tử tứ giác Q9 trong hệ tọa độ tự nhiên ......................................... 16
Hình 2.6. Mặt cắt chi tiết nền đường băng nhiều lớp ............................................ 20
Hình 2.7. Mơ hình tấm Mindlin trên nền nhiều lớp .............................................. 21
Hình 2.8. Nền với nhiều lớp đất khác nhau .......................................................... 23
Hình 2.9. Phần tử tứ giác 9 nút, 2 lớp trong phương pháp MPMM………………25
Hình 2.10. Lưu đồ tính tốn………………………………………………………..42
Hình 3.1. Mơ hình kiểm chứng tấm Mindlin trên nền nhiều lớp chịu tải trọng
tĩnh khi xem tấm xi măng đá và nền là cứng vơ cùng.......................... 47
Hình 3.2.

Sự hội tụ chuyển vị lớn nhất w c tại tâm tấm bê tơng ......................... 48

Hình 3.3.

Mơ hình kiểm chứng tấm Mindlin trên nền nhiều lớp chịu tải trọng
hằng số di động .................................................................................. 50

Hình 3.4.

Sự chênh lệch của chuyển vị lớn nhất tại tâm tấm BT theo các bước
thời gian t theo các phương pháp ..................................................... 51

Hình 3.5.

Sự hội tụ của chuyển vị tại tâm tấm BT theo các bước thời gian t .... 52

Hình 3.6.

Sự hội tụ của chuyển vị tại tâm tấm XMĐ theo các bước thời gian t . 53

Hình 3.7.

Sự hội tụ chuyển vị lớn nhất tại tâm tấm bê tông và tâm tấm xi măng
đá trong trường hợp bước thời gian t  0.001s ............................... 54

Hình 3.8.

Biên độ dao động của tâm tấm bê tơng theo thời gian t ..................... 55

Hình 3.9.

Biên độ dao động của tâm tấm xi măng đá theo thời gian t ................ 55


ix

Hình 3.10. Chuyển vị lớn nhất w c tại tâm tấm bê tông và w g tại tâm tấm xi

măng đá khi tỉ số độ cứng k g / kc thay đổi ........................................ 56
Hình 3.11. So sánh chuyển vị tại tâm tấm BT khi tỉ số độ cứng k g / k c thay đổi ... 57
Hình 3.12. So sánh chuyển vị tại tâm tấm XMĐ khi tỉ số độ cứng k g / k c thay đổi 58
Hình 3.13. So sánh biên độ chuyển vị tại tâm tấm bê tông khi tỉ số độ cứng

k g / k c thay đổi ................................................................................. 59
Hình 3.14. So sánh biên độ chuyển vị tại tâm tấm XMĐ khi tỉ số độ cứng

k g / k c thay đổi .................................................................................. 59
Hình 3.15. Chuyển vị lớn nhất w c tại tâm tấm bê tông và w g tại tâm tấm xi
măng đá khi tỉ số độ cứng kc / k g thay đổi ......................................... 60
Hình 3.16. So sánh chuyển vị tại tâm tấm BT khi tỉ số độ cứng kc / k g thay đổi ... 61
Hình 3.17. So sánh chuyển vị tại tâm tấm XMĐ khi tỉ số độ cứng kc / k g thay đổi 62
Hình 3.18. So sánh biên độ chuyển vị tại tâm tấm bê tông khi tỉ số độ cứng

kc / k g thay đổi .................................................................................... 63
Hình 3.19. So sánh biên độ chuyển vị tại tâm tấm XMĐ khi tỉ số độ cứng

kc / k g thay đổi .................................................................................. 63
Hình 3.20. Chuyển vị lớn nhất w c tại tâm tấm bê tông và w g tại tâm tấm xi
măng đá khi tỉ số độ cản cg / cc thay đổi ............................................. 65
Hình 3.21. So sánh chuyển vị tại tâm tấm BT khi tỉ số độ cứng cg / cc thay đổi .... 66
Hình 3.22. So sánh chuyển vị tại tâm tấm XMĐ khi tỉ số độ cứng cg / cc thay đổi. 66
Hình 3.23. So sánh biên độ chuyển vị tại tâm tấm bê tông khi tỉ số độ cứng

cg / c c thay đổi ................................................................................... 67
Hình 3.24. So sánh biên độ chuyển vị tại tâm tấm XMĐ khi tỉ số độ cứng

cg / c c thay đổi .................................................................................. 67

x

Hình 3.25. Chuyển vị lớn nhất w c tại tâm tấm bê tông và w g tại tâm tấm xi
măng đá khi tỉ số độ cản cc / cg thay đổi ........................................... 68
Hình 3.26. So sánh biên độ chuyển vị tại tâm tấm bê tông khi tỉ số độ cứng

cc / c g thay đổi ................................................................................... 69
Hình 3.27. So sánh biên độ chuyển vị tại tâm tấm XMĐ khi tỉ số độ cứng

cc / c g thay đổi .................................................................................. 70
Hình 3.28. Chuyển vị lớn nhất w c tại tâm tấm bê tông và w g tại tâm tấm xi
măng đá khi tỉ số module đàn hồi Eg / Ec của hai tấm thay đổi ......... 71
Hình 3.29. So sánh chuyển vị tại tâm tấm BT khi tỉ số module Eg / Ec thay đổi ... 72
Hình 3.30. So sánh chuyển vị tại tâm tấm XMĐ khi tỉ số module Eg / Ec thay đổi 73
Hình 3.31. Chuyển vị lớn nhất w c tại tâm tấm bê tông và w g tại tâm tấm xi
măng đá khi tỉ số module đàn hồi Ec / Eg của hai tấm thay đổi ........ 74
Hình 3.32. So sánh chuyển vị tại tâm tấm BT khi tỉ số module Ec / Eg thay đổi ... 75
Hình 3.33. So sánh chuyển vị tại tâm tấm XMĐ khi tỉ số module Ec / Eg thay đổi 75
Hình 3.34. Chuyển vị lớn nhất w c tại tâm tấm bê tông và w g tại tâm tấm xi
măng đá khi tỉ số chiều dày hg / hc của hai tấm thay đổi .................... 77
Hình 3.35. So sánh chuyển vị tại tâm tấm BT khi tỉ số chiều dày hg / hc thay đổi . 78
Hình 3.36. So sánh chuyển vị tại tâm tấm XMĐ khi tỉ số chiều dày hg / hc
thay đổi ............................................................................................. 78
Hình 3.37. Chuyển vị lớn nhất tại tâm tấm bê tông và tại tâm tấm xi măng đá
khi tỉ số chiều dày hc / hg của hai tấm thay đổi .................................. 79
Hình 3.38. So sánh chuyển vị tại tâm tấm BT khi tỉ số chiều dày hc / hg thay đổi . 80
Hình 3.39. So sánh chuyển vị tại tâm tấm XMĐ khi tỉ số chiều dày hc / hg
thay đổi ............................................................................................. 81

xi

Hình 3.40. So sánh chuyển vị tại tâm tấm BT khi giá trị biên độ tải trọng điều hòa
di động Po thay đổi ........................................................................... 82
Hình 3.41. So sánh chuyển vị tại tâm tấm XMĐ khi giá trị biên độ tải trọng điều
hịa di động Po thay đổi..................................................................... 83
Hình 3.42. Chuyển vị lớn nhất tại tâm tấm bê tông theo vận tốc ứng với

cc  1.75 106 , cg  3.79  105 ......................................................... 85
Hình 3.43. Chuyển vị lớn nhất tại tâm tấm xi măng đá theo vận tốc ứng với

cc  1.75 106 , cg  3.79  105 .......................................................... 85
Hình 3.44. Chuyển vị lớn nhất tại tâm tấm bê tông theo vận tốc ứng với

cc  c g  1.75  103 ........................................................................... 86
Hình 3.45. Chuyển vị lớn nhất tại tâm tấm xi măng đá theo vận tốc ứng với

cc  c g  1.75  103 ........................................................................... 87
Hình 3.46. Chuyển vị theo tần số lực kích thích ứng với vận tốc di chuyển V  0 88
Hình 3.47. Hệ số động (theo chuyển vị tại tâm tấm bê tông w c ) theo tần số lực
ứng với các trường hợp nền có hệ số cản khác nhau .......................... 89
Hình 3.48. Hệ số động (theo chuyển vị tại tâm tấm xi măng đá w g ) theo tần số
lực ứng với các trường hợp nền có hệ số cản khác nhau .................... 90
Hình 3.49. So sánh chuyển vị tại tâm tấm BT khi pha của tải trọng thay đổi ......... 91
Hình 3.50. So sánh chuyển vị tại tâm tấm XMĐ khi pha của tải trọng thay đổi ..... 92
Hình 3.51. Chuyển vị lớn nhất của tấm bê tông và tấm xi măng đá theo pha của
tải trọng (độ) ..................................................................................... 93
Hình 3.52. Chuyển vị tại tâm tấm bê tơng w c theo pha của tải trọng ứng với
các trường hợp khoảng cách giữa hai tải trọng ................................... 96

Hình 3.53. Chuyển vị tại tâm tấm xi măng đá w g theo pha của tải trọng ứng với
các trường hợp khoảng cách giữa hai tải trọng ................................... 97


xii

DANH MỤC CÁC BẢNG BIỂU
Bảng 2.1. Tọa độ và trọng số trong phép cầu phương Gauss ................................. 19
Bảng 2.2. Thông số nền đường băng nhiều lớp ..................................................... 20
Bảng 2.3. Thông số của tải trọng .......................................................................... 38
Bảng 2.4. Thông số của tấm bê tông ..................................................................... 39
Bảng 2.5. Thông số liên kết giữa hai tấm .............................................................. 39
Bảng 2.6. Thông số của tấm xi măng đá ............................................................... 39
Bảng 2.7. Thông số nền đất .................................................................................. 39
Bảng 3.1. Thông số của tải trọng .......................................................................... 44
Bảng 3.2. Thông số của tấm bê tông ..................................................................... 45
Bảng 3.3. Thông số liên kết giữa hai tấm .............................................................. 45
Bảng 3.4. Thông số của tấm xi măng đá ............................................................... 45
Bảng 3.5. Thông số nền đất .................................................................................. 45
Bảng 3.6. Sự hội tụ chuyển vị lớn nhất w c (mm) tại tâm tấm bê tông ................. 49
Bảng 3.7. Sai số (%) chuyển vị lớn nhất w c tại tâm tấm bê tông của các
phương pháp với lưới chia 60x60 so với nghiệm giải tích ..................... 49
Bảng 3.8. Sự chênh lệch của chuyển vị lớn nhất w c tại tâm tấm BT theo các
bước thời gian t theo các phương pháp .............................................. 51
Bảng 3.9. Sự hội tụ chuyển vị lớn nhất w c tại tâm tấm bê tông ............................ 53
Bảng 3.10. Sự hội tụ chuyển vị lớn nhất w g tại tâm tấm xi măng đá .................... 53
Bảng 3.11. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số độ cứng k g / kc thay đổi ......... 57
Bảng 3.12. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số độ cứng kc / k g thay đổi ......... 60

Bảng 3.13. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số độ cản cg / cc thay đổi ............. 65


xiii

Bảng 3.14. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số độ cản cc / cg thay đổi ............ 69
Bảng 3.15. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số module đàn hồi Eg / Ec của
hai tấm thay đổi .................................................................................... 72
Bảng 3.16. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số module đàn hồi Ec / Eg của
hai tấm thay đổi .................................................................................... 74
Bảng 3.17. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số chiều dày hg / hc của hai tấm
thay đổi................................................................................................. 77
Bảng 3.18. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi tỉ số chiều dày hc / hg của hai tấm
thay đổi................................................................................................. 80
Bảng 3.19. So sánh (%) chênh lệch chuyển vị lớn nhất w c tại tâm tấm bê tông
và w g tại tâm tấm xi măng đá khi giá trị lực di động Po thay đổi ........ 83
Bảng 3.20. Chuyển vị tại tâm tấm BT và tâm tấm XMĐ ứng với các giá trị pha
của tải trọng .......................................................................................... 92
Bảng 3.21. Chuyển vị (mm) của tâm tấm bê tông wc theo pha của tải trọng ứng
với các trường hợp khoảng cách giữa hai tải trọng ................................ 94
Bảng 3.22. Chuyển vị (mm) của tâm tấm xi măng đá w g theo pha của tải trọng
ứng với các trường hợp khoảng cách giữa hai tải trọng ......................... 95
Bảng 3.23. Sai số chuyển vị (%) giữa hai trường hợp pha   0o và   180o ứng
với các trường hợp khoảng cách giữa hai tải trọng ................................ 96

xiv

MỘT SỐ KÝ HIỆU VIẾT TẮT
Chữ viết tắt
MPMM

Phương pháp phần tử nhiều lớp tấm chuyển động (Multi-Layer
Plate Moving Method)

FEM

Phương pháp phần tử hữu hạn (Finite Element Method)

MEM

Phương pháp phần tử chuyển động (Moving Element Method)

FTM

Phương pháp biến đổi Fourier (Fourier Transform Method)

SIM

Structural Impedance Method

DQM

Differential Quadrature Method

DSC

Discrete Singular Convolution

HDQ

Harmonic Differential Quadrature

EEM

Eigenfunction Expansion Method

BEM

Phần tử biên (Boundary Element Method)

Q9

Phần tử tứ giác 9 nút (Quadrilateral Nine-Node Element)

DOF

Bậc tự do (Degree of Freedom)

BT

Bê tông

XMĐ

Xi măng đá

Ma trận và vec tơ
u

Véctơ chuyển vị tại một điểm bất kỳ của tấm

uc

Véctơ chuyển vị tại một điểm bất kỳ của tấm bê tông

ug

Véctơ chuyển vị tại một điểm bất kỳ của tấm xi măng đá

κ

Véctơ độ cong

d

Véctơ chuyển vị nút của phần tử nhiều lớp tấm chuyển động

γ

Ma trận biến dạng cắt

M

Ma trận khối lượng tổng thể

K

Ma trận độ cứng tổng thể

C

Ma trận cản tổng thể


xv

P

Véctơ tải tổng thể

Me

Ma trận khối lượng phần tử nhiều lớp tấm chuyển động

Ce

Ma trận cản phần tử nhiều lớp tấm chuyển động

Ke

Ma trận độ cứng phần tử nhiều lớp tấm chuyển động

Mec

Ma trận khối lượng phần tử tấm bê tông

Cec

Ma trận cản phần tử tấm bê tông

Kec

Ma trận độ cứng phần tử tấm bê tông

Meg

Ma trận khối lượng phần tử tấm xi măng đá

Ceg

Ma trận cản phần tử tấm xi măng đá

Keg

Ma trận độ cứng phần tử tấm xi măng đá

Meff

Ma trận khối lượng hiệu dụng

Peff

Ma trận tải trọng hiệu dụng

Keff

Ma trận độ cứng hiệu dụng

Ký hiệu

L

Chiều dài tấm theo phương x

B

Chiều dài tấm theo phương y

Lc

Chiều dài tấm bê tông theo phương x

Bc

Chiều dài tấm bê tông theo phương y

Lg

Chiều dài tấm xi măng đá theo phương x

Bg

Chiều dài tấm xi măng đá theo phương y

E

Module đàn hồi của tấm

Ec

Module đàn hồi của tấm bê tông

Eg

Module đàn hồi của tấm xi măng đá

G

Module chống cắt đàn hồi của tấm

Gc

Module chống cắt đàn hồi của tấm bê tông


xvi

Gg

Module chống cắt đàn hồi của tấm xi măng đá



Hệ số Poisson của tấm

c

Hệ số Poisson của tấm bê tông

g

Hệ số Poisson của tấm xi măng đá



Trọng lượng riêng của tấm

c

Trọng lượng riêng của tấm bê tông

g

Trọng lượng riêng của tấm xi măng đá

h

Chiều dày tấm

hc

Chiều dày tấm bê tông

hg

Chiều dày tấm xi măng đá

kc

Hệ số độ cứng liên kết giữa tấm bê tông và tấm xi măng đá

cc

Hệ số độ cản liên kết giữa tấm bê tông và tấm xi măng đá

kg

Hệ số độ cứng nền đất

cg

Hệ số độ cản nền đất

x

Góc xoay của tấm quay quanh trục y

y

Góc xoay của tấm quay quanh trục x

s
u , v, w

Hệ số hiệu chỉnh cắt

uc , vc , wc

Chuyển vị của tấm bê tông theo phương x , y và z

u g , v g , wg

Chuyển vị của tấm xi măng đá theo phương x , y và z

V

Vận tốc của tải trọng điều hòa di động



Tần số của tải trọng điều hòa di động



Pha của tải trọng điều hòa di động

P0

Biên độ của tải trọng điều hòa đi động

Chuyển vị của tấm theo phương x , y và z

Tổng quan

1

CHƯƠNG 1
TỔNG QUAN
1.1.

Giới thiệu
Nhiều nghiên cứu đã được thực hiện gần đây để dự đoán chuyển vị và ứng suất

của nền đường cao tốc, nền đường băng, …Lớp bề mặt của áo đường (lớp bê tông
hoặc lớp bê tông nhựa asphalt) thường được mơ hình dựa trên lý thuyết tấm mỏng
của Kirchhoff và các lớp phía dưới thường được mơ hình là nền đàn hồi, thơng dụng
nhất là nền đàn hồi kiểu Winkler. Tải trọng tác dụng bởi xe chạy, máy bay chuyển
động trên đường cao tốc,… thường được mơ hình là tải trọng tĩnh (tải trọng là hằng
số) trong các trường hợp đó.
Trong thực tế, khi xe cộ chạy trên đường, tải trọng tác dụng lên bề mặt đường
sẽ có cường độ thay đổi theo thời gian do độ nhám của bề mặt đường và hệ thống
động cơ của xe theo Kim và cộng sự (2004) [1]. Do đó xem xe chạy trên đường là tải
trọng điều hịa di động sẽ phản ánh đúng với thực tế hơn so với khi xem xe cộ là tải
trọng hằng số di động.
Trong đó, trong kết cấu đường băng, nền đường thường được cấu tạo nhiều
lớp bao gồm: lớp bê tông, lớp nhựa đường, lớp xi măng đá, đặt trên nền đất được thể
hiện như trong Hình 1.1. theo Wu và cộng sự (2014) [2].

Hình 1.1. Mặt cắt nền đường băng nhiều lớp


Tổng quan

2

Trong Luận văn này, phương pháp nhiều lớp tấm chuyển động MPMM (MultiLayer Plate Moving Method) được sử dụng để giải quyết bài toán tấm trên nền nhiều
lớp chịu tải trọng điều hịa di động và có xét đến sự tương tác giữa các lớp với nhau.
Phương pháp MPMM được đề xuất dựa trên phương pháp phần tử chuyển động
MEM. Các phương trình chuyển động của tấm, các ma trận kết cấu được thiết lập
trong hệ tọa độ tương đối chuyển động cùng vận tốc với vận tốc của lực di chuyển.
Do đó, phương pháp này sẽ tránh được việc cập nhật véctơ tải trọng tương ứng với
mơ hình tấm.
Trong Luận văn, giả thuyết lớp bê tông và lớp xi măng đá như là hai tấm dày
(tấm Mindlin), còn liên kết giữa hai tấm bằng lớp nhựa đường được mơ hình hóa
thành hệ số độ cứng đàn hồi và hệ số độ cản tương tác giữa hai tấm. Các hệ số độ
cứng đàn hồi và hệ số độ cản được tính tốn dựa trên mơ hình của Richart và Lysmer,
có xét đến hệ số điều chỉnh theo mơ hình Whitman theo Shambhu (2009) [3] được
thể hiện như trong Hình 1.2.

Hình 1.2. Mơ hình tải trọng cố định và phần tử tấm chuyển động (MPMM)
1.2.

Tình hình nghiên cứu và tính cấp thiết của đề tài
Cùng với sự phát triển của khoa học kỹ thuật, nhiều phương pháp số được thiết

lập để tính tốn và phân tích ứng xử động của kết cấu tấm chịu tải trọng động.


Tổng quan

3

1.2.1. Các cơng trình nghiên cứu ngồi nước
Bài tốn phân tích ứng xử của kết cấu tấm chịu tác dụng của tải trọng di động
đã được nghiên cứu từ rất lâu. Mathews (1958) [4] đã giải quyết bài toán dầm có
chiều dài vơ hạn trên nền đàn hồi chịu tải trọng di chuyển tùy ý bằng phương pháp
biển đổi FTM (Fourier Transform Method). Phương pháp FTM thực chất là một
phương pháp miền tần số, có thể cho lời giải chính xác nhưng gặp bế tắc khi bài tốn
phức tạp, nhiều bậc tự do, khi tải trọng tác động có xét đến sự thay đổi của gia tốc.
Michael và Edward (1989) [5] đã giải quyết bài toán tấm Kirchoff với điều kiện biên
bất kỳ sử dụng phương pháp SIM (Structural Impedance Method). Liew và cộng sự
(1996) [6] đã giải quyết bài toán tấm Mindlin trên nền Winkler bằng phương pháp
DQM (Differential Quadrature Method) với điều kiện biên tựa đơn, tự do và ngàm.
Gbadeyan và Oni (1995) [7] đã thực hiện phân tích ứng xử động của dầm và tấm chữ
nhật chịu tác động của tải trọng chuyển động dựa trên phương pháp Struble’s hiệu
chỉnh. Tiếp đó, Gbadeyan và Dada (2006) [8] thực hiện phân tích ứng xử động của
tấm chữ nhật Mindlin chịu vật thể chuyển động có khối lượng phân bố đều. Kim và
Rosset (1998) [9] đã nghiên cứu đến trạng thái ứng xử của một tấm vô hạn trên nền
đàn hồi chịu tải trọng chuyển động điều hòa khơng đổi. Kim (2004) [1] đã phân tích
mất ổn định và dao động của tấm Kirchoff trên nền đàn nhớt Winkler dưới tác dụng
của tải trọng động bằng phương pháp biến đổi Fourier Transform. Huang và
Thambiratnam (2001) [10] đã sử dụng phương pháp dải hữu hạn để phân tích ứng xử
của tấm trên nền đàn hồi. Sun (2003) [11] đã xây dựng lời giải giải tích cho tấm
Kirchhoff trên nền đàn nhớt chịu tải điều hòa bằng chuỗi Fourier. Sau đó, Sun (2005)
[12] đã phân tích tấm Kirchoff trên nền đàn nhớt chịu tải trọng điều hòa sử dụng hàm
Green. Civalek và Acar (2007) [13] đã áp dụng phương pháp DSC (Discrete Singular
Convolution) để giải quyết bài toán tấm Mindlin chịu uốn trên nền đàn hồi hai thông
số. Tiếp đến, Civalek (2007) [14] đã phân tích phi tuyến tấm Kirchoff trên nền
Winkler-Pasternak bằng phương pháp kết hợp DSC-HDQ (Harmonic Differential
Quadrature). Javad và cộng sự (2013) [15] đã sử dụng phương pháp EEM
(Eigenfunction Expansion Method) để giải quyết bài toán ổn định và ứng xử động
của tấm Mindlin dưới tác động của tải trọng di động.

Tổng quan

4

Phương pháp phần tử hữu hạn truyền thống FEM (Finite Element Method) đã
được sử dụng để giải quyết nhiều bài toán phức tạp. Jong-Shyong và cộng sự (1987)
[16] đã thực hiện phân tích ứng xử động của tấm chịu tải trọng động bằng phương
pháp phần tử hữu hạn FEM bằng việc rời rạc hóa kết cấu, sử dụng phần tử đẳng tham
số tứ giác kết hợp với tích phân Newmark để giải quyết bài toán động. Mức độ ảnh
hưởng của độ lệch tâm, vận tốc và gia tốc ban đầu và chiều dài của tấm là các yếu tố
quan trọng trong nghiên cứu trên. Musharraf và cộng sự (1991) [17] đã phân tích
phản ứng xử động tấm Mindlin trên nền đàn nhớt dưới tác động của tải trọng di
chuyển. Pan và Atluri (1995) [18] đã giải bài toán ứng xử động trong đường băng
bằng cách sử dụng phương pháp kết hợp giữa phần tử hữu hạn và phần tử biên
FEM/BEM (Boundary Element Method). Xing và Liu (2009) [19] đã trình bày
phương pháp giải quyết bài tốn dao động của tấm chữ nhật Mindlin. Lời giải của bài
toán rất hữu ích cho việc thiết kế các thông số ban đầu cũng như phân tích nó trong
thực tiễn. Li và cộng sự (2013) [20] đã phân tích ứng xử động của tấm chữ nhật nền
đàn nhớt dưới tác dụng của vật di chuyển với vận tốc thay đổi. Phương pháp phần tử
hữu hạn FEM đưa ra lời giải bằng cách rời rạc hóa phần tử tấm thành các phần tử hữu
hạn. Các thành phần chuyển vị, nội lực được tính tốn dựa trên các hàm dạng và
chuyển vị của các nút của phần tử. Phương pháp phần tử hữu hạn FEM gặp khó khăn
khi tải trọng tiến đến gần biên của miền hữu hạn phần tử và di chuyển vượt ra ngồi
biên, ngồi ra phương pháp này u cầu phải ln cập nhật vị trí véctơ tải trọng. Do
đó việc giải quyết bài toán sẽ tốn nhiều thời gian và chi phí.
Để giảm bớt những khó khăn khi giải quyết bài toán chịu tải di động, Koh và
cộng sự (2003) [21] đã sử dụng phương pháp phần tử chuyển động MEM trong việc
khảo sát ứng xử động của tàu cao tốc. Sau đó, Koh và cộng sự (2007) [22] đã khảo

sát đến ứng xử động của nền bán không gian đàn hồi sử dụng phương pháp MEM.
Xu và cộng sự (2009) [23] sử dụng phương pháp MEM để phân tích dao động ngẫu
nhiên của tấm Kirchhoff trên nền Kelvin chịu tải trọng di động sử dụng phần tử 2-D
chuyển động. Lei và Wang (2013) [24] đã đề xuất một cách tiếp cận mới tên là phần
tử khung chuyển động cho đường ray, dựa trên phần tử xe và phần tử nền để đánh giá
ứng xử động của tàu và hệ thống nền ba lớp. Ang và cộng sự (2013) [25] đã sử dụng


Tổng quan

5

MEM để khảo sát đến ứng xử của tàu cao tốc trong khoảng thời gian tăng tốc và giảm
tốc và khảo sát hiện tượng nảy bánh xe của tàu cao tốc. Tran và cộng sự (2013) [26]
đã phân tích động lực học của tàu cao tốc trên nền hai thơng số. Tiếp đó, Tran và
cộng sự (2014) [27] tiếp tục sử dụng MEM để phân tích ứng xử động của tàu cao tốc
chuyển động với vận tốc không đều. Gần đây, Tran và cộng sự (2016) [28] đã phân
tích ứng xử động của tàu cao tốc khi phanh đột ngột. Cao và cộng sự (2018) [29] đã
phân tích ứng xử động của mơ hình đường ray xe lửa cao tốc ba chiều sử dụng phương
pháp MEM. Luong và cộng sự (2018) [30] đã phân tích ứng xử tĩnh và động của tấm
Mindlin trên nền nhớt đàn hồi bằng phương pháp MEM. Cao và cộng sự (2018) [31]
đã dùng phương pháp MEM để phân tích ứng xử động của tấm composite trên nền
Pasternak chịu tải trọng di động.
Trong kết cấu đường băng, nền đường được cấu tạo nhiều lớp bao gồm: Lớp
bê tông, lớp nhựa đường, lớp xi măng đá, tiếp đến là nền đất. Wu và cộng sự (2014)
[2] đã khảo sát ứng xử động của tấm Kirchoff trên nền nhiều lớp dưới tác động của
tải trọng di chuyển. Tuy nhiên, hiện nay chưa có bất cứ nghiên cứu nào về kết cấu
tấm Mindlin trên nền nhiều lớp chịu tải trọng điều hịa di động mà có xét đến tương
tác giữa các lớp với nhau.
1.2.2. Các cơng trình nghiên cứu trong nước

Một số Luận văn cao học ngành xây dựng cơng trình dân dụng và cơng nghiệp
tại trường Đại học Bách Khoa Thành phố Hồ Chí Minh cũng đã giải quyết một số bài
toán tải trọng chuyển động đối với dầm và tấm.
Cường (2011) [32] đã phân tích dao động của tấm Mindlin trên nền đàn nhớt
xét đến khối lượng của vật chuyển động sử dụng phương pháp phần tử hữu hạn FEM.
Duy (2013) [33] phân tích ứng xử động tàu cao tốc có xét đến độ cong thanh ray và
tương tác đất nền. Anh (2013) [34] đã phân tích động lực học tàu cao tốc có xét đến
độ nảy bánh xe và tương tác với đất nền. Hải và cộng sự (2013) [35] đã phân tích ứng
xử tàu cao tốc có xét đến độ cong thanh ray và tương tác với đất nền sử dụng MEM.
Nhựt (2014) [36] đã phân tích ứng xử động của tàu cao tốc bằng mơ hình 3-D sử dụng
phương pháp phần tử chuyển động. Thu (2014) [37] đã phân tích động lực học tàu
cao tốc sử dụng phương pháp phần tử nhiều lớp dầm chuyển động có xét đến tương


Tổng quan

6

tác đất nền. Nhi (2014) [38] đã phân tích động lực học tấm Mindlin trên nền đàn nhớt
chịu tải trọng di động sử dụng phương pháp phần tử 2-D chuyển động. An (2014)
[39] đã phân tích bài tốn tương tác tấm Mindlin trên nền đàn nhớt được gia cường
Top-base chịu tải trọng di động. Xuyên (2015) [40] đã phân tích động lực học tấm
trên nền Top-base chịu tải trọng di động sử dụng phần tử tấm chuyển động. Minh
(2016) [41] đã phân tích động lực học tấm Mindlin trên nền nhiều lớp chịu tải trọng
động. Thân và cộng sự (2016) [42] đã phân tích ứng xử động của tấm lớp composite
lớp trên nền đàn nhớt chịu tải trọng di chuyển sử dụng phương pháp phần tử chuyển
động. Thân và cộng sự [43] đã dùng phương pháp tấm trên nền nhiều lớp để phân
tích nền đường chịu tải trọng hằng số di động .Thân và cộng sự (2018) [44] tiếp tục
phân tích ứng xử động kết cấu dầm trên nền Pasternak dưới tác dụng của tải trọng
chuyển động không đều sử dụng phương pháp phần tử chuyển động cải tiến.

Tuy nhiên, trong các nghiên cứu này, bài toán tải trọng có cường độ thay đổi
theo thời gian tác dụng lên nền có nhiều lớp thì chưa được xét đến. Do đó, để tiếp nối
các nghiên cứu về kết cấu tấm, Luận văn sẽ phân tích tấm Mindlin chịu tải trọng điều
hòa di động sử dụng phương pháp tấm nhiều lớp chuyển động (MPMM) và giải
phương trình chuyển động bằng phương pháp Newmark. Từ đó rút ra các kết luận
quan trọng và đề xuất các giải pháp áp dụng trong thực tế.
1.3.

Mục tiêu và hướng nghiên cứu
Mục tiêu của Luận văn nhằm sử dụng phương pháp phần tử nhiều lớp tấm

chuyển động MPMM (Multi-Layer Plate Moving Method) để phân tích động lực học
tấm Mindlin trên nền nhiều lớp chịu tải trọng điều hòa di động. Các vấn đề nghiên
cứu cụ thể trong phạm vi Luận văn bao gồm:
 Thiết lập phương trình năng lượng của tấm Mindlin, mơ hình nền đường
nhiều lớp và thiết lập công thức ma trận kết cấu tấm Mindlin trên nền nhiều
lớp sử dụng phương pháp phần tử nhiều lớp tấm chuyển động MPMM.
 Sử dụng phương pháp Newmark để giải quyết phương trình chuyển động.
 Phát triển thuật tốn, sử dụng ngơn ngữ lập trình Matlab để giải hệ phương
trình động của bài tốn.


Tổng quan

7

 Kiểm tra độ tin cậy của chương trình tính bằng cách so sánh kết quả của
Luận văn với kết quả của các phương pháp khác.
 Tiến hành thực hiện các ví dụ số nhằm khảo sát ảnh hưởng của các nhân tố
quan trọng đến ứng xử động của kết cấu tấm, từ đó rút ra kết luận và kiến nghị.

1.4.

Cấu trúc của Luận văn
Nội dung trong Luận văn được trình bày như sau:
 Chương 1: Giới thiệu tổng quan về tải trọng điều hịa di động, tình hình
nghiên cứu của các tác giả trong và ngoài nước về các loại tấm chịu tải trọng
di động, cũng như mục tiêu và hướng nghiên cứu của đề tài.
 Chương 2: Trình bày tổng quan lý thuyết tấm Mindlin, phần tử đẳng tham
số, mơ hình nền đường nhiều lớp và phương pháp phần tử nhiều lớp tấm
chuyển động MPMM, phương pháp Newmark dạng gia tốc trung bình để giải
quyết phương trình chuyển động.
 Chương 3: Kiểm tra độ tin cậy của chương trình tính bằng cách so sánh với
các phương pháp khác. Trình bày các kết quả phân tích số được tính tốn bằng
ngơn ngữ lập trình Matlab.
 Chương 4: Đưa ra một số kết luận quan trọng đạt được trong Luận văn và
kiến nghị hướng phát triển của đề tài trong tương lai.
 Tài liệu tham khảo: trích dẫn các tài liệu liên quan phục vụ cho mục đích
nghiên cứu của đề tài.