Header Page 1 of 148.

BỘ GIÁO DỤC VÀ ĐÀO TẠO

BỘ QUỐC PHÒNG

HỌC VIỆN KỸ THUẬT QUÂN SỰ

PHẠM VĂN ĐẠT

PHÂN TÍCH KẾT CẤU DÀN CHỊU TẢI TRỌNG
TĨNH THEO SƠ ĐỒ BIẾN DẠNG

LUẬN ÁN TIẾN SĨ KỸ THUẬT

HÀ NỘI - NĂM 2015
Footer Page 1 of 148.


Header Page 2 of 148.

BỘ GIÁO DỤC VÀ ĐÀO TẠO

BỘ QUỐC PHÒNG

HỌC VIỆN KỸ THUẬT QUÂN SỰ

PHẠM VĂN ĐẠT

PHÂN TÍCH KẾT CẤU DÀN CHỊU TẢI TRỌNG
TĨNH THEO SƠ ĐỒ BIẾN DẠNG


Phạm Văn Đạt

Footer Page 3 of 148.


Header Page 4 of 148.

II

LỜI CẢM ƠN
Tác giả xin bày tỏ lòng biết ơn vô hạn ñến GS. TSKH. Hà Huy Cương
nhà khoa học chân chính. Thầy ñã luôn chỉ bảo, ñộng viên tận tình hướng dẫn
giúp ñỡ tác giả nâng cao kiến thức khoa học ñể hoàn thành luận án này.
Tác giả xin trân trọng cảm ơn PGS. TS. Nguyễn Phương Thành ñã giúp
ñỡ cho tác giả nhiều chỉ dẫn khoa học và tạo mọi ñiều kiện thuận lợi ñể tác
giả hoàn thành luận án này.
Tác giả xin trân trọng cảm ơn các thầy cô, các bạn ñồng nghiệp trong Bộ
môn Sức bền – Cơ kết cấu, trong Khoa Xây dựng và Khoa Tại chức, Trường
Đại học Kiến trúc Hà Nội ñã luôn quan tâm, giúp ñỡ và ñộng viên tác giả
hoàn thành luận án.
Tác giả xin trân trọng cảm ơn Ban giám ñốc, Khoa sau ñại học, Viện các
công trình ñặc biệt và tập thể các thầy cô trong Bộ môn Cơ sở kỹ thuật công
trình Học viện Kỹ thuật Quân sự ñã giúp ñỡ tác giả trong quá trình học tập
nghiên cứu tại Học viện.
Tác giả xin cảm ơn ñối với người thân trong Gia ñình ñã ñộng viên,
khích lệ và không ngại vất vả trong công việc gia ñình ñể tác giả yên tâm
hoàn thành luận án.
Tác giả luận án


pháp nguyên lý cực trị Gauss
2.1.1 Phân tích tuyến tính kết cấu dàn theo cách thứ nhất
2.1.2 Phân tích tuyến tính kết cấu dàn theo cách thứ hai
2.1.3 Phân tích phi tuyến hình học kết cấu dàn theo cách thứ nhất
2.1.4 Phân tích phi tuyến hình học kết cấu dàn theo cách thứ hai
2.2 Phương pháp xác ñịnh các thành phần chuyển vị tại các nút dàn và
nội lực trong các thanh dàn ñối với bài toán dàn phi tuyến hình học
2.3 Một số kết quả nghiên cứu bài toán kết cấu dàn
2.3.1 Tính toán dàn theo cách thứ nhất
2.3.2 Tính toán dàn theo cách thứ hai
2.3.3 Ảnh hưởng của thông số vật liệu ñến ñộ chênh lệch kết quả phân
tích nội lực trong các thanh dàn giữa PTTT và PTPTHH
2.3.4 Ảnh hưởng của giá trị tải trọng tác dụng ñến ñộ chênh lệch kết
quả nội lực trong các thanh dàn giữa PTTT và PTPTHH
2.4 Kết luận chương

Footer Page 5 of 148.

Trang
I
II
III
VI
VIII
XII
XIV
1
5
7
17

3.2.1 Tính toán dàn vòm phẳng tĩnh ñịnh trong, siêu tĩnh ngoài
3.2.2 Ảnh hưởng ñộ thoải của dàn vòm phẳng tĩnh ñịnh trong, siêu tĩnh
ngoài ñến PTCL chuyển vị, nội lực giữa PTPTHH và PTTT
3.3 Phân tích phi tuyến hình học dàn vòm phẳng siêu tĩnh trong, tĩnh
ñịnh ngoài
3.3.1 Tính toán dàn vòm phẳng siêu tĩnh trong, tĩnh ñịnh ngoài
3.3.2 Ảnh hưởng ñộ thoải của dàn vòm phẳng siêu tĩnh trong, tĩnh ñịnh
ngoài ñến PTCL chuyển vị, nội lực giữa PTPTHH và PTTT
3.4 Phân tích phi tuyến hình học dàn vòm phẳng siêu tĩnh trong và siêu
tĩnh ngoài
3.4.1 Tính toán dàn vòm phẳng siêu tĩnh trong và siêu tĩnh ngoài
3.4.2 Ảnh hưởng ñộ thoải của dàn vòm phẳng siêu tĩnh trong và siêu
tĩnh ngoài ñến PTCL chuyển vị, nội lực giữa PTPTHH và PTTT
3.5 Kết luận chương
CHƯƠNG 4
PHÂN TÍCH PHI TUYẾN HÌNH HỌC KẾT CẤU DÀN KHÔNG
GIAN
4.1 Phân tích phi tuyến hình học dàn cầu không gian một lớp
4.1.1 Tính toán phi tuyến hình học dàn Kiewitt 8
4.1.2 So sánh kết quả tính toán chuyển vị, nội lực giữa PTTT và
PTPTHH
4.1.3 Ảnh hưởng ñộ thoải của dàn cầu không gian K8 ñến PTCL
chuyển vị, nội lực giữa PTTT và PTPTHH
4.2 Phân tích phi tuyến hình học kết cấu dàn vòm không gian một lớp
4.2.1 Tính toán dàn vòm không gian một lớp loại 1
4.2.2 So sánh kết quả tính toán chuyển vị, nội lực giữa PTTT và
PTPTHH
4.2.3 Ảnh hưởng ñộ thoải của dàn vòm không gian một lớp loại 1 ñến
PTCL chuyển vị, nội lực giữa PTTT và PTPTHH
4.3 Kết luận chương

V

CHƯƠNG 5
TÍNH TOÁN ỔN ĐỊNH PHI TUYẾN HÌNH HỌC KẾT CẤU
DÀN VÒM PHẲNG
5.1 Phương pháp chuyển vị cưỡng bức ñể xác ñịnh tải trọng hạn trong
bài toán nén dọc trục
5.1.1 Bài toán ổn ñịnh thanh chịu nén
5.1.2 Phương pháp chuyển vị cưỡng bức
5.1.3 Phương pháp phần tử hữu hạn ñể xác ñịnh tải trọng tới hạn thanh
hai ñầu khớp chịu nén dọc trục
5.2 Phương pháp xác ñịnh tải trọng tới hạn lên kết cấu dàn có kể ñến
tính phi tuyến hình học
5.3 Xác ñịnh tải trọng tới hạn lên dàn vòm phẳng tĩnh ñịnh chịu tải
trọng thẳng ñứng tại nút dàn vòm
5.3.1 Ví dụ phân tích
5.3.2 Nghiên cứu ảnh hưởng ñộ thoải của dàn vòm phẳng tĩnh ñịnh
ñến giá trị tải trọng tới hạn tác dụng lên dàn
5.4 Tính toán ổn ñịnh dàn vòm phẳng tĩnh ñịnh trong, siêu tĩnh ngoài
chịu tải trọng thẳng ñứng tại nút dàn
5.4.1 Ví dụ phân tích
5.4.2 Ảnh hưởng ñộ thoải của dàn vòm phẳng tĩnh ñịnh trong, siêu tĩnh
ngoài ñến giá trị tải trọng tới hạn tác dụng lên dàn vòm
5.5 Tính toán ổn ñịnh dàn vòm phẳng siêu tĩnh trong, tĩnh ñịnh ngoài
chịu tải trọng thẳng ñứng tại nút dàn
5.5.1 Ví dụ phân tích
5.5.2 Ảnh hưởng ñộ thoải của dàn vòm phẳng siêu tĩnh trong, tĩnh ñịnh
ngoài ñến giá trị tải trọng tới hạn tác dụng lên dàn vòm
5.6 Tính toán ổn ñịnh dàn vòm phẳng siêu tĩnh trong và siêu tĩnh ngoài
chịu tải trọng thẳng ñứng tại nút dàn vòm

136
140
141
142


Header Page 8 of 148.

VI

DANH MỤC CÁC BẢNG
Số hiệu
Bảng 2.1
Bảng 2.2
Bảng 2.3
Bảng 2.4
Bảng 2.5
Bảng 3.1
Bảng 3.2
Bảng 3.3
Bảng 3.4
Bảng 3.5
Bảng 3.6
Bảng 3.7
Bảng 3.8
Bảng 3.9
Bảng 3.10
Bảng 3.11
Bảng 3.12
Bảng 3.13

Kết quả chuyển vị theo hai phương tại các nút dàn
66
Kết quả so sánh nội lực trong các thanh dàn giữa PTTT
67
và PTPTHH
PTCL nội lực trong các thanh dàn ứng với các giá trị
70
k=f/l khác nhau
Kết quả chuyển vị theo hai phương tại các nút dàn
73
Kết quả so sánh nội lực trong các thanh dàn giữa PTTT
74
và PTPTHH
PTCL nội lực trong các thanh dàn ứng với các giá trị
78
k=f/l khác nhau
Kết quả chuyển vị theo hai phương tại các nút dàn
82
Kết quả so sánh nội lực trong các thanh dàn giữa PTTT
83
và PTPTHH
PTCL nội lực trong các thanh dàn ứng với các giá trị
86
k=f/l khác nhau
Kết quả PTCL nội lực trong các thanh giữa PTTT và
97
PTPTHH của dàn cầu Kiewitt 8 ứng với các giá trị k=f/l
khác nhau
Kết quả PTCL nội lực trong các thanh giữa PTTT và 106
PTPTHH của dàn vòm không gian một lớp loại 1 ứng


VIII

DANH MỤC CÁC HÌNH VẼ
Số hiệu

Nội dung hình vẽ

Trang

Hình 1.1

Sân vận ñộng Astrodome

6

Hình 1.2

Nhà thi ñấu Superdome

6

Hình 1.3

Nhà thi ñấu Nagoya Dome

6

Hình 1.4


Sơ ñồ chuyển vị của nút thanh trong hệ phẳng PTTT

26

Hình 2.4

28

Hình 2.7

Sơ ñồ chuyển vị của nút thanh trong hệ không gian
PTTT
Sơ ñồ chuyển vị của nút thanh trong hệ phẳng
PTPTHH
Sơ ñồ chuyển vị của nút thanh trong hệ không gian
PTPTHH
Sơ ñồ khối chương trình.

Hình 2.8

Dàn ví dụ 2.2

37

Hình 2.9

Hình dạng dàn trước và sau khi biến dạng ví dụ 2.2

41


Hình 2.15

Nội lực thanh 1, 2

48

Hình 2.16

Nội lực thanh 3, 4

48

Hình 2.17

Nội lực thanh 5

48

Hình 2.5
Hình 2.6

Footer Page 10 of 148.

31
33
36


Header Page 11 of 148.


50

Hình 2.23

Hình dạng dàn

50

Hình 3.1

54

Hình 3.2
Hình 3.3
Hinh 3.4
Hình 3.5
Hình 3.6
Hình 3.7
Hinh 3.8
Hình 3.9
Hình 3.10

Dàn vòm phẳng tĩnh ñịnh chịu tải trọng thẳng ñứng tại
các nút dàn
Vị trí các nút dàn vòm
Hình dạng dàn trước và sau khi biến dạng
Chuyển vị theo phương trục x
Chuyển vị theo phương trục y
Dàn vòm phẳng tĩnh ñịnh trong, siêu tĩnh ngoài
Hình dạng dàn trước và sau khi biến dạng

80

Hình 3.15

Hình dạng dàn trước và sau khi biến dạng

81

Hinh 3.16
Hình 3.17
Hình 4.1

Chuyển vị theo phương trục x
Chuyển vị theo phương trục y
Một số dạng kết cấu dàn cầu không gian

85
85
90

Hình 4.2

Dàn cầu không gian K8

91

Hình 4.3

Số hiệu nút của dàn K8



Hình 4.8

Hình dạng kết cấu dàn trước và sau biến dạng khi
k=1/8
Chuyển vị của các nút dàn theo phương x

95

Hình 4.9

Chuyển vị của các nút dàn theo phương y

95

Hình 4.10

Chuyển vị của các nút dàn theo phương z

95

Hình 4.11

Biểu ñồ so sánh nội lực

95

Hình 4.12

Biểu ñồ chuyển vị theo phương x


101

Hình 4.18

Kết quả các thành phần chuyển vị tại các nút dàn (cm)

103

Hình 4.19

Kết quả nội lực trong các thanh dàn (kN)

103

Hình 4.20

104

Hình 4.21

Hình dạng kết cấu dàn trước và sau biến dạng khi
k=1/3
Chuyển vị của nút dàn theo phương x

105

Hình 4.22

Chuyển vị của nút dàn theo phương y


108

Hình 5.1

Thanh ñầu hai ñầu khớp chịu nén ñúng tâm

112

Hình 5.2

Phần tử dầm

114

Hình 5.3

Ổn ñịnh thanh hai ñầu khớp

115

Hình 5.4

Ví dụ 5.1

121

Footer Page 12 of 148.




Header Page 14 of 148.

XII

DANH MỤC CÁC KÝ HIỆU
Ký hiệu

Đại lượng

A

Diện tích mặt cắt ngang của thanh dàn

B

Bề rộng của dàn vòm không gian 1 lớp

c

Hàm ràng buộc dạng bất ñẳng thức

C

Số liên kết nối ñất

ceq

Hàm ràng buộc dạng ñẳng thức



z

Tổng hình chiếu sai số theo phương x
Tổng hình chiếu sai số theo phương y
Tổng hình chiếu sai số theo phương z

G

Mô ñun ñàn hồi trượt

g

Hàm ràng buộc

h

Chiều cao dàn

I

Ma trận ñơn vị

k

Độ thoải của dàn

[K]

Ma trận ñộ cứng kết cấu

Mô men

N

Lực dọc trong thanh

P

Tải trọng tác dụng

Pth

Tải trọng tới hạn

Q

Lực cắt

r

Bán kính cong của dàn

u

Thành phần chuyển vị theo phương x

v

Thành phần chuyển vị theo phương y


Biến dạng dài tỉ ñối

λ

Thừa số Largrage

χ

Biến dạng uốn

θ1

Góc xoay tại nút thứ nhất của phần tử thanh

θ2

Góc xoay tại nút thứ hai của phần tử thanh

Footer Page 15 of 148.


Header Page 16 of 148.

XIV

DANH MỤC CÁC CHỮ VIẾT TẮT
PTCL
PTPTHH
PTTT


mại, xưởng sửa chữa bảo dưỡng máy bay v.v…
Trước ñây việc tính toán phân tích nội lực cho kết cấu dàn thường ñược
thực hiện tính toán bằng thủ công với các phương pháp ñơn giản như: Phương
pháp tách mắt, Phương pháp mặt cắt ñơn giản, Phương pháp mặt cắt phối hợp,
Phương pháp họa ñồ - Giản ñồ Maxwell-Cremona v.v… Hiện nay do sự phát
triển của công nghệ tin học ñiện tử nên việc tính toán ñơn giản và thuận tiện
hơn rất nhiều nhờ các phần mềm phân tích tính toán ứng dụng ñược viết dựa
theo phương pháp phần tử hữu hạn như phần mềm Sap, Etabs v.v…, ñặc biệt
các phần mềm này có thể phân tích tính toán với các kết cấu siêu tĩnh bậc cao.
Tuy nhiên khi áp dụng các phương pháp này ñể tính toán kết cấu dàn thực tế

Footer Page 17 of 148.


Header Page 18 of 148.

2

thường giả thiết chuyển vị nút dàn là bé, tức là bỏ qua sự thay ñổi góc của các
trục thanh dàn giữa trước và sau khi dàn biến dạng. Do kết cấu dàn ngày càng
mỏng, vượt khẩu ñộ lớn và vật liệu có ñộ bền cao, nên khi chịu lực làm cho
góc của các trục thanh dàn giữa trước và sau khi dàn biến dạng thay ñổi.
Chính vì thế, kết quả phân tích tuyến tính kết cấu dàn hiện nay là chưa sát với
sự làm việc thực tế của kết cấu do chưa kể ñến ảnh hưởng của sự thay ñổi góc
các thanh dàn trong quá trình kết cấu biến dạng. Ngoài ra, khi tính toán ổn

ñịnh cho kết cấu dàn hiện nay thường mới chỉ nghiên cứu tính toán ổn ñịnh
tổng thể hoặc tính toán ổn ñịnh cục bộ tuyến tính cho kết cấu dàn. Với lý do
trên, luận án nghiên cứu với ñề tài: “Phân tích kết cấu dàn chịu tải trọng



Phương pháp nghiên cứu:
Dựa trên phương pháp nguyên lý cực trị Gauss do GS.TSKH. Hà Huy
Cương ñề xuất và kết hợp với các phương pháp quy hoạch toán học. Phương
pháp nguyên lý cực trị Gauss có cách nhìn ñơn giản ñể phân tích các bài toán
kết cấu và người kỹ sư có thể dễ dàng áp dụng ñể tính toán trong các bài toán
kết cấu thực tế.

Phạm vi nghiên cứu:
Giới thiệu phương pháp mới ñể tính toán nội lực, chuyển vị và ổn ñịnh
cục bộ của dàn vòm có xét ñến tính phi tuyến hình học và vật liệu làm việc
trong giai ñoạn ñàn hồi chịu tác dụng của tải trọng tĩnh tại các nút dàn.

Ý nghĩa khoa học của luận án:
Xét ñược tính phi tuyến hình học của kết cấu dàn khi phân tích nội lực,
chuyển vị và ổn ñịnh chịu tải trọng tĩnh là vấn ñề rất khoa học và có ý nghĩa
thực tiễn.

Bố cục của luận án
Ngoài phần mở ñầu, phần kết luận, danh mục các công trình nghiên cứu

của tác giả liên quan ñến nội dung luận án, tài liệu tham khảo và phụ lục; nội
dung chính của luận án ñược bố cục trong 5 chương:
- Chương 1 Tổng quan về phân tích tính toán kết cấu dàn: Trình bày các

ñặc ñiểm, ứng dụng của kết cấu dàn cũng như tổng quan về tình hình nghiên
cứu kết cấu dàn ở trong nước và trên thế giới. Đồng thời trong chương này
trình bày một số khái niệm về toán quy hoạch và ổn ñịnh kết cấu dàn. Cuối
chương tác giả ñưa ra các vấn ñề cụ thể giải quyết của luận án.
- Chương 2 Cơ sở lý thuyết phân tích phi tuyến hình học kết cấu dàn:



Header Page 21 of 148.

5

CHƯƠNG 1
TỔNG QUAN VỀ PHÂN TÍCH TÍNH TOÁN KẾT CẤU DÀN
1.1 Đặc ñiểm và ứng dụng của kết cấu dàn
Kết cấu dàn là kết cấu ñược tạo thành từ các thanh liên kết với nhau tại
các nút dàn, nút dàn phải nằm tại giao ñiểm của các trục thanh. Khi lực chỉ

ñặt tại nút thì các thanh dàn chủ yếu làm việc chịu kéo hoặc nén, do ñó ta có
thể coi các nút dàn là khớp. Trong phụ lục 5 của luận án tác giả cũng ñã phân
tích cùng một kết cấu dàn cho hai trường hợp (khi coi nút dàn là khớp và khi
coi nút dàn là liên kết cứng), kết quả quả cho thấy lực dọc trong các thanh dàn
khi coi nút dàn là khớp lý tưởng thường nhỏ hơn khi coi nút dàn là liên kết
cứng khoảng 5%. Do kết cấu dàn khi chịu lực, các thanh chủ yếu chỉ chịu kéo
hoặc nén nên tận dụng hết ñược khả năng làm việc của vật liệu. Vì vậy kết
cấu dàn là kết cấu tiết kiệm vật liệu và về phương diện kiến trúc có thể tạo

ñược nhiều hình dáng khác nhau, nên kết cấu dàn ñược sử dụng nhiều trong
các công trình cầu, dàn khoan, cột truyền tải ñiện và làm kết cấu mái che cho
các công trình nhà thi ñấu, sân vận ñộng, nhà hát, sân bay v.v…[31], [33],
[51], [64], [65].
Kết cấu dàn ñầu tiên trên thế giới ñược xây dựng năm 1863 là công
trình Schwedler Dome tại Berlin do kỹ sư Schwedler người Đức thiết kế, có
dạng kết cấu vòm ñược tạo bởi các lưới ô tam giác và vượt ñược khẩu ñộ 30m
[62]. Đến năm 1889 tại Pari Pháp xây dựng tháp Eiffel nằm cạnh sông Seine
có chiều cao 325 m trở thành biểu tượng của kinh ñô ánh sáng. Năm 1898 tại

Nagoya Dome có sức chứa 40.500 người với kích thước khẩu ñộ trên 180m
(hình 1.3).

Footer Page 22 of 148.


Header Page 23 of 148.

7

Năm 2007 Trung Quốc ñã xây dựng nhà hát lớn tại Bắc Kinh dạng hình
Elipsoid, với kích thước một chiều 144m và một chiều 212m. Chiều cao của
công trình 46m và công trình có sức chứa 5.452 người (hình1.4).
Ngoài ứng dụng làm kết cấu
cho các công trình nhịp lớn như ñã
kể trên, kết cấu dàn còn có tác dụng
giảm chấn cho các kết cấu công
trình chịu ñộng ñất. Khi có ñộng ñất

Vïng tiªu t¸n n¨ng l−îng

xẩy ra thì trên kết cấu dàn STMFs
(Special Truss moment frames) xuất

Hình 1.5 Kết cấu STMFs

hiện các vị trí biến dạng dẻo (vùng tiêu tán năng lượng) như hình 1.5, làm
tăng khả năng giảm chấn cho công trình [28].
Ngoài ra, do cách tính ñơn giản của dàn nên có thể dùng sơ ñồ dàn ảo


có ý nghĩa thực tiễn. Qua các tài liệu tham khảo có thể thấy, hiện nay khi
phân tích kết cấu dàn không gian thường sử dụng một trong các phương pháp
sau: phương pháp cân bằng nút [14], [66]; phương pháp mặt cắt ñơn giản [14],
[66]; phương pháp mặt cắt phối hợp [14], [66]; phương pháp phần tử hữu hạn
[12], [29], [66]; phương pháp sai phân hữu hạn [66]; phương pháp mô hình
tấm mỏng tương ñương [66].

Phân tích phi tuyến kết cấu dàn: Kết cấu dàn thường nhẹ, ngày càng mỏng
và vượt khẩu ñộ lớn nên khi phân tích kết cấu dàn ñã có một số nghiên cứu
xem xét ảnh hưởng của tính phi tuyến hình học cũng như phi tuyến vật liệu

ñến sự làm việc của kết cấu dàn:
- Năm 2001 khi trình bày phân tích phi tuyến hình học kết cấu dàn
Carlos A.Fellippa sử dụng phương pháp phần tử hữu hạn với phần tử hữu hạn
không tuyến tính, biến dạng dài tỉ ñối của thanh dàn ñược xác ñịnh :
(s) 2
ij

(0) 2
ij
(0) 2
ij

(l ) − (l )
ε=
2 (l )

(1.1)

Dựa theo nguyên lý thế năng biến dạng cực tiểu sẽ xây dựng ñược


ổn ñịnh tổng thể phi tuyến hình học của dàn Mises chịu tải trọng thẳng ñứng.
Trong nghiên cứu này, tác giả ñã sử dụng phương pháp cân bằng nút ñể thiết
lập ñược ñường cân bằng cho bài toán [34].
- Năm 2012 M.Greco và các cộng sự ñã nghiên cứu phân tích phi tuyến
hình học của kết cấu dàn theo hai cách: cách thứ nhất là xây dựng theo vị trí
nút dựa trên nguyên lý công ảo; cách thứ hai phân biến dạng của kết cấu dàn
ra làm hai thành phần là biến dạng thể tích tương ñối và biến dạng quay cứng
xung quanh ba trục tọa ñộ. Trong cả hai cách của các tác giả là cuối cùng ñưa
về dạng các phương trình cân bằng phi tuyến, ñể giải các phương trình này
các tác giả ñã sử dụng phương pháp lặp chiều dài cung ñể giải [38].
Ngoài các nghiên cứu phân tích phi tuyến hình học vừa trình bày trên
còn một số các nghiên cứu khác như [27], [37], [45], [46] v.v… và trong hầu
hết các nghiên cứu này các tác giả thường sử dụng phương pháp phần tử hữu

Footer Page 25 of 148.