Tìm nghiệm giải tích gần đúng của bài toán tĩnh và động lực tấm composite lớp trên nền đàn hồi theo mô hình vinkler

MỞ ĐẦU
Vật liệu composite ngày nay được sử dụng rộng rãi khi thiết kế chế tạo
những kết cấu hàng không, tên lửa, vũ trụ, tàu thuyền. Composite được ứng dụng
ngày càng nhiều trong những lĩnh vực khác nhau của ngành chế tạo máy và nền
kinh tế quốc dân.
Composite được ứng dụng và phát triển như vậy vì chúng rất nhẹ và bền.
Để có thể thiết kế tối ưu vật liệu và các kết cấu composite, cần thiết phải hiểu rõ
bản chất và những quy luật ứng xử cơ học khá phức tạp của loại vật liệu này.
Trong thực tế thường gặp các kết cấu đặt tiếp xúc trên bề mặt một môi
trường hoặc vật thể đàn hồi khác, ví dụ dầm móng đặt trên nền đất, cầu phao phà
đặt trên mặt nước. Bài toán xác định nội lực và chuyển vị của kết cấu trên nền
đàn hồi là dạng bài toán siêu tĩnh, trong đó phản lực nền là một hệ lực phân bố
liên tục trên bề mặt tiếp xúc, phụ thuộc vào biến dạng của kết cấu cũng như quan
niệm về mô hình nền. Trong luận văn này ta sử dụng mô hình đơn giản, thường
dùng trong kỹ thuật là mô hình Vinkler. Theo đó, cường độ phản lực của nền tại
một điểm tỷ lệ thuận với độ lún của nền tại điểm đó. Nếu kí hiệu p là áp suất
phản lực, y là độ lún, K là hệ số nền thì p = Ky. Thứ nguyên của hệ số nền là
[Lực/(chiều dài)3].
Dao động là một hiện tượng phổ biến trong tự nhiên và trong kỹ thuật. Các
máy, các phương tiện giao thông vận tải, các tòa nhà cao tầng, các cầu, các mạch
điện là các hệ dao động trong kỹ thuật. Nghiên cứu về dao động ngày nay trở
thành bộ phận không thể thiếu được cho tất cả các kết cấu, công trình.
Trong [1], [6] đã nghiên cứu bài toán dao động của vỏ trụ và vỏ thoải
composite có gân gia cường. Dao động phi tuyến của tấm composite lớp có gân
gia cường được tính toán trong [2]. Trong [5] đã nghiên cứu bài toán phi tuyến,
đưa ra các hệ thức tính toán tĩnh và động cho vỏ thoải composite hai độ cong bất
kỳ. Trong [7] đã tính toán dao động vỏ thoải composite.

1

Mục đích của luận văn là tìm nghiệm giải tích gần đúng của bài toán
tĩnh và động lực tấm composite lớp trên nền đàn hồi theo mô hình Vinkler.
Từ hệ phương trình cân bằng đã sử dụng hàm ứng suất và phương pháp BubnovGalerkin để nhận được phương trình dao động phi tuyến của tấm. Lời giải số tìm
được theo phương pháp bước lặp và sơ đồ tính toán Newmark, đã xem xét quan
hệ tần số - biên độ dao động phi tuyến, ảnh hưởng của hệ số nền và tần số dao
động ngoại lực đến lời giải bài toán động lực của tấm. Báo cáo đã sử dụng phương
pháp phần tử hữu hạn tính toán độ võng tấm composite lớp trên nền đàn hồi. Đã so
sánh kết quả thu được theo hai phương pháp giải tích và phần tử hữu hạn.
Luận văn gồm 3 chương:
Chương 1. Các hệ thức cơ sở của tấm composite lớp trên nền đàn hồi.
Chương 2. Phương pháp phần tử hữu hạn cho tấm composite lớp trên nền đàn hồi.
Chương 3. Tính toán số cho tấm composite lớp trên nền đàn hồi.
Qua đây, em xin được gửi lời cảm ơn chân thành đến thầy Vũ Đỗ Long người
đã tận tình giúp đỡ em trong suốt quá trình thực hiện luận văn. Em cũng xin bày tỏ
lòng biết ơn đến các thầy cô trong bộ môn Cơ học và các thầy cô trong khoa Toán –
Cơ – Tin học đã trang bị kiến thức giúp em hoàn thành luận văn này.
Các kết quả chính của luận văn đã được trình bày và thảo luận tại Hội nghị
khoa học toàn quốc “Cơ học vật rắn biến dạng lần thứ 10”. Tác giả đã nhận được
những góp ý bổ ích từ các thành viên của Hội nghị. Tuy nhiên do bước đầu tiếp cận
nghiên cứu khoa học về lĩnh vực vật liệu composite, chắc chắn luận văn không
tránh khỏi những thiếu sót. Vì vậy, tác giả rất mong tiếp tục nhận được những đánh
giá và góp ý của các thầy cô giáo và các bạn để luận văn được hoàn thiện hơn.
Hà Nội, tháng 11 năm 2010
Học viên
Nguyễn Thị Huệ

2


CHƯƠNG 1

CÁC HỆ THỨC CƠ SỞ CỦA TẤM COMPOSITE LỚP TRÊN
NỀN ĐÀN HỒI
1.1. Phương trình tổng quát của tấm composite lớp trên nền đàn hồi
1.1.1. Mối liên hệ chuyển dịch – biến dạng của tấm composite lớp
Xét một tấm composite lớp có x1 , x2 là các trục tọa độ nằm trong mặt
phẳng giữa theo các cạnh, còn x3 ≡ z hướng theo phương pháp tuyến với mặt
giữa (Hình 1).

x1

x3

x2

0
Hình 1

Theo lý thuyết của Kirchhoff – Love mối liên hệ phi tuyến dịch chuyển –
biến dạng của tấm:

ε11 = ε10 + zφ1
ε 22 = ε 20 + zφ2
γ 12 = ε 60 + zφ6

3


Trong đó:
2

ε10

∂u 1  ∂w 
=
+ 
÷
∂x1 2  ∂x1 

2

∂v 1  ∂w 
=
+ 
÷
∂x2 2  ∂x2 
∂u ∂v ∂w ∂w
ε 60 =
+
+
∂x2 ∂x1 ∂x1 ∂x2
∂ 2w
φ1 = − 2
∂x1
∂2w
φ2 = − 2
∂x2
∂2w
φ6 = −
∂x1∂x2

ε 20

(1.1)

Còn u , v, w là chuyển vị của phương ngang, phương dọc và độ võng của
các điểm giữa thuộc mặt phẳng của tấm; ε10 , ε 20 , ε 60 là các biến dạng tại mặt giữa;

φ1,φ2 ,φ6 là các biến thiên độ cong của tấm. Chúng thỏa mãn phương trình tương
thích biến dạng:
2

∂ 2ε10 ∂ 2ε 20 ∂ 2ε 60  ∂ 2 w  ∂ 2 w ∂ 2 w
+ 2 −
=
÷ − 2
∂x1∂x2  ∂x1∂x2 
∂x22
∂x1
∂x1 ∂x22

(1.2)

1.1.2. Quan hệ ứng suất biến dạng của tấm composite lớp
Sử dụng giả thiết Kirchhoff có thể bỏ qua thành phần ứng suất vuông góc
với mặt giữa: σ 13 = σ 23 = σ 33 = 0
Liên hệ ứng suất biến dạng của lớp composite thứ k của tấm là [5]:
σ 11 
 
σ 22 
σ 

 12 

(k )

 Q11 Q12 Q16 


=  Q12 Q22 Q26 


Q 16 Q26 Q66 

4

(k )

 ε1 
 
ε 2 
ε 
 6

(k )

(1.3)


Trong đó ký hiệu các thành phần biến dạng trong mặt phẳng lớp thứ k:
ε11 ≡ ε1 , ε 22 ≡ ε 2 , ε 6 ≡ ε12
Trường hợp phương của sợi lệch một góc θ với trục x1 của tấm, thay ma

k
k
trận Qij bằng ma trận Qijk . Trong đó Qijk tính qua Qij theo công thức [4]:

Q11 = Q11 cos 4 θ + 2 ( Q12 + 2Q66 ) sin 2 θ cos 2 θ + Q22 sin 4 θ

(

Q12 = ( Q11 + Q22 − 4Q66 ) sin 2 θ cos 2 θ + Q12 sin 4 θ + cos 4 θ

)

Q16 = ( Q11 − Q12 − 2Q66 ) sin θ cos3 θ + ( Q12 − Q22 + 2Q66 ) sin 3 θ cosθ
Q22 = Q11 sin 4 θ + 2 ( Q12 + 2Q66 ) sin 2 θ cos 2 θ + Q22 cos 4 θ

(1.4)

Q26 = ( Q11 − Q12 − 2Q66 ) sin 3 θ cosθ + ( Q12 − Q22 + 2Q66 ) sin θ cos3 θ

(

Q66 = ( Q11 + Q22 − 2Q12 − 2Q66 ) sin 2 θ cos 2 θ + Q66 sin 4 θ + cos 4 θ

)

Biểu thức các hằng số độ cứng qua các mô đun đàn hồi trong hệ trục chính
như sau:
Q11 =

E1

E 2 ,
1 − 2 ν12
E1

Q12 =

Q16 = 0

E2
ν
E2 2 12 ,
1 − ν 12
E1

Q26 = 0

Q22 =

E2
E 2
1 − 2 ν12
E1

(1.5)

Q66 = G12

trong đó: E1 , E2 là các môđun đàn hồi của tấm theo phương trục chính của lớp
vật liệu composite; ν 12 là hệ số Poisson của vật liệu, G12 là môđun trượt trong
hệ trục chính của lớp vật liệu.

Lực pháp, lực tiếp, mômen uốn, mômen xoắn được xác định theo công thức:
h

N1 =
−

∫h σ11dz

h

N2 =
−

h

2

M1 =
−

2

2

∫h

2

h

σ 22dz

M2 =
−

2

∫h

2

∫h

h

σ 11zdz

N6 =
−

2

∫h

h

σ 22 zdz

M6 =
−

2

5

2

2

2

∫h

σ 12dz

2

(1.6)

σ 12 zdz


N12 = N 21 = N 6 , M12 = M 21 = M 6

Ở đây:

Thay (1.3) vào (1.6) ta được:
4

h

N1 = ∑

2

∫

k =1 − h
4

=∑

h

k =1 −

h

=

2

(

2

∫h

( k)

Q11 ε1 + Q12

( k)

(

2

2

)

( k)

ε 2 + Q16 ε 6 dz

)

(

)

(

Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ
1
1
12
2
2

16
6
6
 11

0
0
0


∫h Q11 ( ε1 + zφ1 ) + Q12 ( ε 2 + zφ2 ) + Q16 ( ε 6 + zφ6 ) dz +
4

+

( 1)

h

4

( 1)

(

)

( 1)

(

)

(

∫h

Q ( 3) ε 0 + zφ + Q ( 3) ε 0 + zφ + Q ( 3) ε 0 + zφ
1
1
12
2
2
16
6
6
 11

+
−

−h

+
−

∫h

4
4

2

(

)

(

)

(

) dz +

(

)

(

)

(

) dz

Q ( 4 ) ε 0 + zφ + Q ( 4 ) ε 0 + zφ + Q ( 4 ) ε 0 + zφ
1
1

12
2
2
16
6
6
 11

= A11ε10 + A12ε 20 + A16ε 60 + B11φ1 + B12φ2 + B16φ6
4

h

N2 = ∑

2

∫

k =1 − h

(

N6 = ∑

h

k =1 −

2

∫h

( k)

)

( k)

A22ε 20 + A26ε 60 + B12φ1 + B22φ2 + B26φ6

(Q

16

2

( k)

Q12 ε1 + Q22 ε 2 + Q26 ε 6 dz

2
0
A12ε1 +
4

) dz +

∫

Q ( 2) ε 0 + zφ + Q ( 2 ) ε 0 + zφ + Q ( 2 ) ε 0 + zφ
1
1
12
2
2
16
6
6
 11

0
0

=

) dz

( k)

( k)

)

( k)

ε1 + Q26 ε 2 + Q66 ε 6 dz

= A16ε10 + A26ε 20 + A66ε 60 + B16φ1 + B26φ2 + B66φ6

6


h

4

M1 = ∑

k =1 −

=∑

k =1 −

h

=

∫h

h

4

2

∫
h

4
h

2

(Q

11

2

2

∫h

( k)

( k)

(

2

)

( k)

ε1 + Q12 ε 2 + Q16 ε 6 zdz

)

(

)

(

Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ
1
1
12
2
2
16
6
6
 11

(

)

(

)

(

Q ( 1) ε 0 + zφ + Q ( 1) ε 0 + zφ + Q ( 1) ε 0 + zφ
1

1
12
2
2
16
6
6
 11
4

(

)

(

)

(

) zdz +

( 2)
( 2) 0
( 2)
+ ∫ Q11 ε10 + zφ1 + Q12
ε 2 + zφ2 + Q16 ε 60 + zφ6

0
0

+
−

∫h

(

+
−

∫h

(

)

(

) zdz +
)

Q ( 3) ε 0 + zφ + Q ( 3) ε 0 + zφ + Q ( 3) ε 0 + zφ zdz +
1
1
12
2
2
16
6

6 
 11


4

−h

)

) zdz

4

2

(

)

(

)

(

)

Q ( 4 ) ε 0 + zφ + Q ( 4 ) ε 0 + zφ + Q ( 4 ) ε 0 + zφ zdz
1

1
12
2
2
16
6
6 
 11


= B11ε10 + B12ε 20 + B16ε 60 + D11φ1 + D12φ1 + D16φ6
4

h

M2 = ∑

k =1 −

2

∫h

2

(

)

(

)

(

) zdz

(

) zdz

Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ
1
1
22
2
2
26
6
6
 12

= B12ε10 + B22ε 20 + B26ε 60 + D12φ1 + D22φ1 + D26φ6
4

M6 = ∑

h

k =1 −

2

∫h

2

(

)

(

)

Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ + Q ( k ) ε 0 + zφ
1
1
26
2
2
66
6
6
 16

= B16ε10 + B26ε 20 + B66ε 60 + D16φ1 + D26φ1 + D66φ6

7

N zk +1

( Aij , Bij , Dij ) = ∑ ∫ Qij ( k ) ( 1, z, z 2 ) dz

Với:

k =1 zk

(i, j = 1, 2, 6)

(1.7)

Mối quan hệ giữa lực pháp, lực tiếp, momen uốn, momen xoắn và biến
dạng, biến thiên độ cong được viết dưới dạng ma trận như sau:
 N1   A11
N  A
 2   12
 N 6   A16
 =
 M1   B11
 M 2   B12
  
 M 6   B16

A12
A22
A26
B12
B22

A16 B11 B12
A26 B12 B22
A66 B16 B26
B16 D11 D12
B26 D12 D22

B26

B66

D16

D26

0
B16  ε1 
÷ 
B26 ÷ε 20 
B66 ÷ε 0 
÷ 6 
D16 ÷ φ 
1
D26 ÷φ 
÷ 2 
D66 ÷
 φ6 

(1.8)

với: N = 4 là số lớp của tấm;  Aij  ,  Bij  ,  Dij  lần lượt là ma trận độ cứng dãn
nén, độ cứng tương tác dãn – uốn – xoắn và độ cứng uốn của tấm composite lớp.
Giả thiết tấm xếp lớp đối xứng qua mặt giữa ta có  Bij  = 0 và xem các đại
lượng A16 , A26 , D16 , D26 là nhỏ có thể bỏ qua.
Khai triển (1.8) ta được biểu thức lực màng của tấm composite lớp:
N1 = A11ε10 + A12ε 20
N 2 = A12ε10 + A22ε 20

(1.9)

N 6 = A66ε 60
Giải ngược lại suy ra:

ε10 =
Trong đó:


1 
A12
1 
A12 
0
N
−
N
ε
=
N
−
N1 ÷,

,

÷

1
2
2
2
*
A 22
A
E1* 
E

11

2

E1*

2
2
A11 A22 − A12
A11 A22 − A12
*
=
, E2 =
, G* = A66
A22
A11

8

ε 60 =

1
N6
G*


Và mômen trong của tấm composite lớp được tính theo công thức:
 ∂ 2w
∂ 2w 
M1 = D11φ1 + D12φ2 = − D1  2 + µ2 2 ÷
∂x2 
 ∂x1
 ∂ 2w
∂ 2w 
M 2 = D12φ1 + D22φ2 = − D2  2 + µ1 2 ÷
∂x1 
 ∂x2
2
∂ w
M 6 = D66φ6 = −2 Dk
∂x1∂x2

(1.10)

Trong đó:
D1 = D11 ,

µ2 =

D12
,
D11

D2 = D22 ,

µ1 =

D12
,
D22

Dk = D66

µ1 µ2
=
D1 D2

D3 = 2 Dk + D1µ 2 = 2 Dk + D2 µ1
1.1.3. Phương trình chuyển động của tấm composite lớp trên nền đàn hồi
Phương trình chuyển động của tấm composite lớp trên nền đàn hồi với hệ
số nền K theo mô hình Vinkler được viết như sau [5]:
∂N1 ∂N 6
∂ 2u
∂3w
+
= J O 2 − J1

∂x1 ∂x2
∂t
∂x1∂t 2
∂N 6 ∂N 2
∂ 2v
∂ 3w
+
= J O 2 − J1
∂x1 ∂x2
∂t
∂x2∂t 2

(1.11)
∂ M6 ∂ M2
∂ M1
∂  ∂w
∂w  ∂  ∂w
∂w 
+
2
+
+
N
+
N
+
N
+
N


÷

÷
1
6
6
2
∂x1∂x2
∂x1  ∂x1
∂x2  ∂x2 
∂x1
∂x2 
∂x12
∂x22
 ∂ 3u
 ∂ 4w
∂2w
∂ 3v 
∂ 4w 
= J O 2 + J1 
+
− J 2  2 2 + 2 2 ÷− q (t ) + Kw
2
2÷
∂t
∂
x
∂
t
∂

x
∂
t
∂x2 ∂t 
2
 1

 ∂x1 ∂t
2

2

2

N zk +1

Trong đó J i được xác định theo công thức: J i = ∑

∫
k =1

zk

ρ ( k ) z i dz

ρ ( k ) là mật độ khối lượng của lớp thứ k, q(t ) là lực phân bố

9

1.2. Nghiệm của bài toán
Giả thiết lực ngang q(t ) phân bố đều và mật độ khối lượng của lớp thứ k là
hằng số. Khi đó ta có:
N zk +1

J1 = ∑

∫
k =1

ρ ( k ) zdz = 0

zk

Theo Volmir [8] thì các số hạng quán tính trong hai phương trình đầu của
(1.11) có thể bỏ qua. Do vậy phương trình chuyển động có dạng:
∂N1 ∂N 6
+
=0
∂x1 ∂x2

(1.12)

∂N 6 ∂N 2
+
=0
∂x1 ∂x2
∂ 2 M1
∂ 2M 6 ∂2M 2
∂  ∂w

∂w 
+2
+
+
+ N6
 N1
÷+
2
2
∂x1∂x2
∂x1  ∂x1
∂x2 
∂x1
∂x2
 ∂ 4w
∂  ∂w
∂w 
∂ 2w
∂ 4w 
+
+ N2
 N6
÷ = J O 2 − J 2  2 2 + 2 2 ÷− q (t ) + Kw
∂x2 
∂x1
∂x2 
∂t
∂x2 ∂t 
 ∂x1 ∂t

(1.13)

(1.14)

Phương trình (1.12), (1.13) thỏa mãn khi đưa vào hàm ứng suất φ dạng:
∂ 2ϕ
N1 = hσ11 = 2
∂x2
∂ 2ϕ
N 2 = hσ 22 = 2
∂x1
∂ 2ϕ
N 6 = hσ 12 = −
∂x1∂x2

10

(1.15)


Thay (1.9) và (1.15) vào phương trình tương thích biến dạng ta được:
4
∂ 2ε10 1  ∂ 4ϕ
* ∂ ϕ 
=
−
ν

÷
1

∂x22
E1*  ∂x24
∂x12∂x22 
4
∂ 2ε 20 1  ∂ 4ϕ
* ∂ ϕ 
= *  4 −ν 2 2 2 ÷
∂x12
E2  ∂x1
∂x1 ∂x2 

∂ε 60
1 ∂ 3ϕ
=− * 2
∂x1
G ∂x1 ∂x2
∂ 2ε 60
1 ∂ 4ϕ
=− * 2 2
∂x1∂x2
G ∂x1 ∂x2
∂ 2ε10 ∂ 2ε 20 ∂ 2ε 60
1 ∂ 4ϕ 1 ∂ 4ϕ  1 ν 1* ν 2*  ∂ 4ϕ
+ 2 −
=
+
+ −
−
Suy ra:
÷

∂x1∂x2 E1* ∂x24 E2* ∂x14  G E1* E2*  ∂x12∂x22
∂x22
∂x1
Vì

ν1* ν*2
=
E1* E2*

nên

1 ν 1* ν 2* 1
ν1*
−
−
= −2 *
G E1* E2* G
E1

Do đó phương trình tương thích biến dạng trở thành:
2

1 ∂ 4ϕ  1
ν1*  ∂ 4ϕ
1 ∂ 4ϕ  ∂ 2 w  ∂ 2 w ∂ 2 w
+
−2 *÷ 2 2 + * 4 =
÷ − 2
E2* ∂x14  G*
E1  ∂x1 ∂x2 E1 ∂x2  ∂x1∂x2 

∂x1 ∂x22
Thay (1.10) và (1.15) vào (1.14) ta được:
 ∂4w
∂ 2 M1
∂4w 
=
−
D
+
µ
÷
1
2
4
∂x12
∂x12∂x22 
 ∂x1
 ∂4w
∂2M 2
∂4w 
= − D2  4 + µ1 2 2 ÷
∂x22
∂x1 ∂x2 
 ∂x2
∂ 2M 6
∂ 4w
2
= −4 Dk 2 2
∂x1∂x2
∂x1 ∂x2

11

(1.16)


∂  ∂w
∂w  ∂  ∂ 2ϕ ∂w
∂ 2ϕ ∂w 
N
+
N
=
+

÷=
 1
÷
6
∂x1  ∂x1
∂x2  ∂x1  ∂x22 ∂x1 ∂x1∂x2 ∂x2 
∂ 3ϕ ∂w ∂ 2ϕ ∂ 2 w
∂ 3ϕ ∂w
∂ 2ϕ ∂ 2 w
=
+
−
−
∂x1∂x22 ∂x1 ∂x22 ∂x12 ∂x12∂x2 ∂x2 ∂x1∂x2 ∂x1∂x2
∂  ∂w

∂w  ∂  ∂ 2ϕ ∂w ∂ 2ϕ ∂w 
+ N2
+
−
÷=
 N6
÷=
∂x2 
∂x1
∂x2  ∂x2  ∂x1∂x2 ∂x1 ∂x12 ∂x2 
∂ 3ϕ ∂w
∂ 2ϕ ∂ 2 w
∂ 3ϕ ∂w ∂ 2ϕ ∂ 2 w
=−
−
+
+
∂x1∂x22 ∂x1 ∂x1∂x2 ∂x1∂x2 ∂x12∂x2 ∂x2 ∂x12 ∂x22
Phương trình (1.14) trở thành:
∂ 4w
∂4w
∂ 4w
∂4w
∂4w
∂ 2ϕ ∂ 2 w
− D1 4 − D1µ2 2 2 − D2 4 − D2 µ1 2 2 − 4D k 2 2 + 2
+
∂x1
∂x1 ∂x2
∂x2

∂x1 ∂x2
∂x1 ∂x2 ∂x1 ∂x22
 ∂ 4w
∂ 2ϕ ∂ 2 w
∂ 2ϕ ∂ 2 w
∂ 2w
∂ 4w 
+ 2
−2
= J O 2 − J 2  2 2 + 2 2 ÷− q (t ) + Kw
∂x1∂x2 ∂x1∂x2
∂x2 ∂x12
∂t
∂x2 ∂t 
 ∂x1 ∂t
hay
∂ 4w
∂ 4 w  ∂ 2ϕ ∂ 2 w ∂ 2ϕ ∂ 2 w
∂ 2ϕ ∂ 2 w 
∂ 2w
D1 4 + D2 4 −  2
+
−2
÷+ J O 2
∂x1∂x2 ∂x1∂x2 
∂x1
∂x2  ∂x1 ∂x22 ∂x22 ∂x12
∂t
 ∂4w
∂4w 

∂4w
− J 2  2 2 + 2 2 ÷− q (t ) + Kw + ( D1µ2 + D2 µ1 + 4D k ) 2 2 = 0
∂x2 ∂t 
∂x1 ∂x2
 ∂x1 ∂t
Do đó (1.14) được viết dưới dạng
∂ 4w
∂ 4 w  ∂ 2ϕ ∂ 2 w ∂ 2ϕ ∂ 2 w
∂ 2ϕ ∂ 2 w 
∂ 2w
D1 4 + D2 4 −  2
+
−2
÷+ J O 2
∂x1∂x2 ∂x1∂x2 
∂x1
∂x2  ∂x1 ∂x22 ∂x22 ∂x12
∂t
(1.17)
 ∂4w
∂4w 
∂ 4w
− J 2  2 2 + 2 2 ÷− q (t ) + Kw + 2D3 2 2 = 0
∂x2 ∂t 
∂x1 ∂x2
 ∂x1 ∂t
Cuối cùng bài toán bao gồm hai phương trình đạo hàm riêng phi tuyến (1.16),
(1.17).

12

Điều kiện biên tựa bản lề thỏa mãn nếu chọn hàm độ võng dạng:
w = f ( t ) sin

π x1 π x2
sin
a
b

(1.18)

trong đó a và b là độ dài các cạnh của tấm, f ( t ) là độ võng cực đại.
Thế (1.18) vào vế phải của phương trình (1.16) ta có
2

4
 ∂2w 
2 π
2 π x1
2 π x2
=
f
c
os
c
os

÷
a

b
a 2b 2
 ∂x1∂x2 
4
∂2w ∂ 2w
πx
πx
2 π
= f 2 2 sin 2 1 sin 2 2
2
2
a
b
∂x1 ∂x2
ab

Suy ra
2

 ∂2w  ∂2w ∂ 2w
=

÷ − 2
2
∂
x
∂
x
∂
x

∂
x
1
2
 1 2
4
1 π 
2π x1 
2π x2  
2π x1 
2π x2  
= f2
1
+
c
os
1
+
c
os
−
1
−
c
os
1
−
c
os


÷
÷

÷
÷
2 a 2b 2 
a 
b  
a 
b  
1 π4 
2π x1
2π x2 
= f2
c
os
+
c
os

÷
2 a 2b 2 
a
b 
Vậy (1.16) trở thành:
4
1 ∂ 4ϕ  1
ν1*  ∂ 4ϕ
1 ∂ 4ϕ
2π x1

2π x2 

21 π
+ * −2 * ÷ 2 2 + * 4 = f
cos
+ cos
÷
*
4
2 2
2a b 
a
b 
E2 ∂x1  G
E1  ∂x1 ∂x2 E1 ∂x2

Khi đó nghiệm của phương trình sẽ là:

ϕ = Acos

2π x1
2π x2
2π x1
2π x2
+ Bcos
+ C sin
sin
a
b
a

b

13

(1.19)


Thế (1.19) vào (1.16) ta được:
4

4

4

4

∂4w
2π x1
π x1 π x2
 2π 
π 
=
A
c
os
+
C
sin
sin


÷

÷
a
a
a
a
b
∂x14


 
∂4w
2π x2
πx
πx
 2π 
π 
= B
+ C  ÷ sin 1 sin 2
÷ cos
4
b
a
b
∂x2
 b 
b
∂ 4ϕ
π4

π x1 π x2
=
C
sin
sin
a
b
∂x12∂x22
a 2b 2
Suy ra
1
E2*

4

4

2π x1 1  2π 
2π x2
 2π 
A
+ * B
+
÷ cos
÷ cos
a
a
b
b
E





1

4
 1  π 4  1
ν1*  π 4
1  π   π x1
πx
+C  *  ÷ +  * − 2 * ÷ 2 2 + *  ÷  sin
sin 2
a
b
E1  a b
E1  b  
 E2  a   G

= f2

1 π4 
2π x1
2π x2 
c
os
+
c
os


÷
2 a 2b 2 
a
b 

Do đó
E2* f 2a 2
A=
32b 2

E1* f 2b 2
B=
32a 2

C =0

(1.20)

Thế (1.18), (1.19) vào (1.17) ta được:
D  πx
πx
πx
πx
 D 2D
W = π 4 f  41 + 2 32 + 42 ÷sin 1 sin 2 + Kf (t )sin 1 sin 2 +
a
b
a
b
a b

b 
a
4
π 
2π x1
2π x2  π x1
πx
− f (t ) 2 2  4 Acos
+ 4 Bcos
sin 2 +
÷sin
a
b 
a
b
ab 

1  d 2 f
πx
πx
 1
+  J O + J 2π 2  2 + 2 ÷ 2 sin 1 sin 2 − q(t ) = 0
a
b
b   dt
a

Đây là phương trình cân bằng của tấm composite lớp trên nền đàn hồi.

14

(1.21)


1.2.1. Bài toán tĩnh
d2 f
= 0 , lực phân bố đều q(t ) = q0 = const
dt 2
Khi đó phương trình cân bằng của tấm trên nền đàn hồi có dạng:
Ta có

D  πx
πx
πx
πx
 D 2D
W = π 4 f  41 + 2 32 + 42 ÷sin 1 sin 2 + Kf (t )sin 1 sin 2 +
a
b
a
b
a b
b 
a
4
4π 
2π x1
2π x2  π x1
πx
− f (t ) 2 2  Acos

+ Bcos
sin 2 − q0 = 0
÷sin
a
b 
a
b
ab 
Áp dụng phương pháp Bubnov – Galerkin:
ab

∫∫ W sin
00

π x1 π x2
sin
dx1dx2 = 0
a
b

ta được
ab

ab

 4  D1 2D3 D2 
 2 π x1 2 π x2
π x1 π x2
π
f

+
+
+
Kf
sin
sin
dx
dx
−
q
sin
sin
dx1dx2

÷
1
2
0
∫∫   a 4 a 2b2 b4  
∫∫
a
b
a
b
00
00
ab

π4 
2π x1

2π x2  2 π x1 2 π x2
− ∫∫ 4 f 2 2  Acos
+ Bcos
sin
dx1dx2 = 0
÷sin
a
b
a
b
a
b


00

(1.22)

Sau một số phép biến đổi và sử dụng (1.20) khi đó phương trình (1.22) trở thành:
16q
( K + m1 ) f + m3 f 3 − 20 = 0
(1.23)
π
trong đó
D 
π 4  E2* E1* 
 D 2D
m1 = π 4  41 + 2 32 + 42 ÷, m3 =
 4 + 4÷
16

a
a
b
b
a 


b
Đây là phương trình cân bằng xác định độ võng cực đại theo phương pháp lý
thuyết bản vỏ mỏng trong trường hợp không có dao động.
1.2.2. Bài toán động
d2 f
≠ 0 , lực phân bố đều q(t ) = q0 sin ( ωt )
dt 2
Khi đó phương trình cân bằng của tấm trên nền đàn hồi có dạng:

Ta có

15


D  πx
πx
πx
πx
 D 2D
W = π 4 f  41 + 2 32 + 42 ÷sin 1 sin 2 + Kf (t )sin 1 sin 2 +
a
b
a

b
a b
b 
a
π4 
2π x1
2π x2  π x1
πx
− f (t ) 2 2  4 Acos
+ 4 Bcos
sin 2 +
÷sin
a
b 
a
b
ab 
2

1  d f
πx
πx
 1
+  J O + J 2π 2  2 + 2 ÷ 2 sin 1 sin 2 − q(t ) = 0
a
b
b   dt
a

Áp dụng phương pháp Bubnov – Galerkin:

ab

∫∫ W sin
00

π x1 π x2
sin
dx1dx2 = 0
a
b

ta được
ab



∫∫ π
00
ab

4


D 
πx
πx
 D 2D
f  41 + 2 32 + 42 ÷+ Kf  sin 2 1 sin 2 2 dx1dx2 +
a
b

a b
b 
a



1  d2 f
πx
πx
2 1
+ ∫∫  J O + J 2π  2 + 2 ÷÷ 2 sin 2 1 sin 2 2 dx1dx2 +
a
b
b   dt
a
0 0
ab

πx
πx
− ∫∫ q(t )sin 1 sin 2 dx1dx2 +
a
b
00

(1.24)

ab

π4 

2π x1
2π x2  2 π x1 2 π x2
− ∫∫ 4 f 2 2  Acos
+ Bcos
sin
dx1dx2 = 0
÷sin
a
b
a
b
a
b


00
Sau một số phép biến đổi và sử dụng (1.20) khi đó phương trình (1.24) trở
thành:
d2 f
16q(t )
m 2 + ( K + m1 ) f + m3 f 3 −
=0
dt
π2

(1.25)

trong đó:
2D3 D2 
1

1 
π 4  E2* E1* 
4  D1
m = J O + J 2π  2 + 2 ÷, m1 = π  4 + 2 2 + 4 ÷, m3 =
+

÷
16  b 4 a 4 
b 
a b
b 
a
a
2

Đây là phương trình dao động của tấm composite trên nền đàn hồi xác định
độ võng cực đại theo phương pháp lý thuyết bản vỏ mỏng.

CHƯƠNG 2

16


PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN CHO TẤM
COMPOSITE LỚP TRÊN NỀN ĐÀN HỒI
2.1. Khái niệm về phương pháp phần tử hữu hạn
Phương pháp phần tử hữu hạn (PP PTHH) là một phương pháp số đặc biệt
có hiệu quả để tìm dạng gần đúng của một hàm chưa biết trong miền xác định V
của nó. Tuy nhiên PP PTHH không tìm dạng xấp xỉ của hàm cần tìm trên toàn
miền V mà chỉ trong từng miền con Ve (phần tử) thuộc miền xác định V. Do đó

phương pháp này rất thích hợp với hàng loạt bài toán vật lý và kỹ thuật trong đó
hàm cần tìm được xác định trên những miền phức tạp gồm nhiều vùng nhỏ có
đặc tính hình học, vật lý khác nhau chịu những điều kiện biên khác nhau.
Trong PP PTHH miền V được chia thành một số hữu hạn các miền con, gọi
là phần tử. Các phần tử này được nối kết với nhau tại các điểm định trước trên
biên phần tử, gọi là nút. Trong phạm vi mỗi phần tử đại lượng cần tìm được lấy
xấp xỉ trong một hàm đơn giản gọi là các hàm xấp xỉ. Và các hàm xấp xỉ này
được biểu diễn qua các giá trị của hàm (và có khi cả các giá trị của đạo hàm của
nó) tại các điểm nút trên phần tử. Các giá trị này được gọi là các bậc tự do của
phần tử và được xem là ẩn số cần tìm của bài toán.
2.2. Tính toán tấm composite lớp trên nền đàn hồi theo phương pháp phần
tử hữu hạn
Bài toán: Xét tấm composite chữ nhật N lớp trên nền đàn hồi có kích thước
(a x b), độ dày h, chịu lực phân bố đều q(x,y) = q = const.
Bước 1: Rời rạc hóa miền vật thể.
Chia tấm thành Ne phần tử hình chữ nhật Ve (e = 1…Ne) có kích thước bằng
nhau. Chọn hệ tọa độ Đề-các có x, y là các trục tọa độ nằm trong mặt phẳng giữa
theo các cạnh, trục z hướng theo phương pháp tuyến với mặt giữa.

17


Giả thiết tấm có các lớp được xếp đối xứng qua mặt trung bình. Khi đó xét
điểm M(x,y) nằm cách mặt trung hòa khoảng z. Khi đó chuyển vị tại M theo
phương Ox, Oy tương ứng là: u = − z

∂w
∂w
, v = −z
∂y

∂x

Tại mỗi điểm M(x,y) trong phần tử e giả sử có các thành phần chuyển vị là
{w(x, y), θx(x, y), θy(x, y)}.
Mỗi phần tử hình chữ nhật có bốn nút, mỗi nút có ba bậc tự do. Do đó phần
tử có mười hai bậc tự do nên việc chọn hàm xấp xỉ phải đảm bảo có mười hai
tham số.
Bước 2: Chọn hàm xấp xỉ
Chọn hàm độ võng dưới dạng đa thức gồm mười hai số hạng [3]:
w( M ) = w( x, y ) = a1 + a2 x + a3 y + a4 x 2 + a5 xy + a6 y 2 + a7 x 3 + a8 x 2 y + a9 xy 2
+ a10 y 3 + a11x 3 y + a12 xy 3
Góc xoay của tấm tại M trên mặt giữa theo các phương Ox, Oy tương ứng:
∂w
= a2 + 2a 4 x + a5 y + 3a 7 x 2 + 2a 8 xy + a9 y 2 + 3a11x 2 y + a12 y 3
∂x
∂w
θy =
= a3 + a 5 x + 2a6 y + a 8 x 2 + 2a 9 xy + 3a10 y 2 + a11x 3 + 3a12 xy 2
∂y

θx =

Chuyển dịch của phần tử
a
2
2
3
2
2
y 3 x3 y xy 3   1 

 w  1 x y x xy y x x y xy
  a2 
  
2
2
2
3
{ u} e = θ x  = 0 1 0 2 y 0 3x 2xy y 0 3x y y   
... 
2
2
3
2 
θ y  0 0
1
0
x
2
y
0
x
2x
y
3
y
x
3x
y
  
 a 

 12 
=  F ( x, y )  { a} e
( 2.1)
Thay tọa độ các nút i, j, k, l đã biết vào phương trình (2.1) ta được véctơ
chuyển vị nút :

18


{ q} e

 F ( xi , yi ) 
 F (x , y )
j
j 
=
{ a} = [ A] { a} e
 F ( xk , yk )  e


 F ( xl , yl ) 

với

 F ( xi , yi ) 
 F (x , y )
[ A] =  j j 
F ( xk , yk )



 F ( xl , yl ) 

{ a} e = [ A] −1 { q} e

Suy ra :

(2.2)

Từ (2.1) và (2.2) có :

{ u} e = [ F ] [ A] { q} e
−1

 N1 
=  N 2  { q} e =  N ( x, y )  { q} e
 
 N 3 

trong đó: [A] là ma trận tọa độ nút phần tử, [N(x,y)] là ma trận hàm dạng
Bước 3 : Biểu diễn biến dạng và ứng suất qua chuyển vị nút
Trạng thái biến dạng phần tử tính theo hệ thức Cauchy:
∂u
∂ 2w
εx =
= − z 2 = zφ x
∂x
∂x
∂v
∂ 2w
ε y = = − z 2 = zφ y

∂y
∂y
∂u ∂v
∂ 2w
γ xy =
+
= −2z
= zφxy
∂y ∂x
∂x∂y
Liên hệ ma trận biến dạng và ma trận độ cong như sau:
ε 
φ 
x
 
 x
{ ε} =  ε y  = z φy  = z[ φ ]
 
 
γ
xy
 
φxy 

19

(2.3)


với

 ∂2w 
 ∂ 2 N1 
 2 

2 
 ∂x 
 ∂x 
 ∂ 2 w 
 ∂ 2 N 
[ φ ] = −  2  = −  21 { q} e = − [ B ] { q} e
 ∂y 
 ∂y 
 2∂ 2 w 
 2∂ 2 N 
1



 ∂x∂y 
 ∂x∂y 

(2.4)

Liên hệ ứng suất - biến dạng được biểu diễn theo định luật Hooke:
σ 
ε 
x
 
 x

{ σ } = σ y  = Qij   ε y  = Qij  { ε }
 
 
τ xy 
γ xy 
Công thức xác định mômen uốn, mômen xoắn là:
h

Mx =
−

2

∫h

h

zσ x dz , M y =
−

2

2

∫h

h

zσ y dz , M xy = M yx =

( k)

Suy ra

 M  h σ 
2  x
 x
 M y  = ∫ z σ y 

 −h 2  
 M xy 
τ xy 

h

dz =
−

2

∫h

Qij 
 

( k)

∫h

zτ xy dz

(2.5)

{ φ} z 2dz = [ D ] { φ}

(2.6)

−

2

2

2

2

(do { φ} = { φ ( x, y ) } )
trong đó [ D ] là ma trận độ cứng uốn của tấm composite lớp:
 D11
[ D ] =  D12
 D16

D12
D22
D26

h
D16 
2

N = 4 zk +1
( k) 2
( k)

D26  =  Dij  = ∫ Qij  z dz = ∑ ∫ Qij  z 2dz
k =1 zk
−h
D66 
2

20


Các phần tử Dij được tính như sau:
h

D11 =

∫h Q11
h

D12 =

2

∫h Q12

h

z 2dz + ∫ Q12

2

( 1)

z 2dz +

h

( 1)

( 1)

2

4
2

∫
h

4

z dz +

( 2)

h /4

∫ Q16

0
h
4

( 1)

0
h
4

Q66 z dz + ∫ Q66
2

0

∫h

Q12

∫h

Q16

∫h

Q22

z 2dz +

z dz +
−

( 2)

( 3)

∫h

z dz +
−

∫h

Q11

z 2 dz +

Q12

∫h

Q16

∫h

Q22

Q26
Q66

( 3)

∫h

2

z dz +

4

z 2dz +

∫h

z 2dz
z 2dz

( 4)

z 2dz

2
4

Q26

−
2
−h

4
−

4

( 4)

2
−h
4

−h

z 2dz

( 4)

∫h

−

( 3)

( 4)

2
−h
4

4

0

2

z 2dz +

( 3)

4

2
−h
4

−

0

2

z 2 dz +

∫h

−

4
0

−

( 2)

( 3)

4
0

z 2dz +

2

−

4
0

−

( 2)

( 3)

Q11 z dz +

z 2dz +

( 2)

−h

∫h

−

0
h
4

4
2

2

0

4

2

0

−

∫ Q26 z dz + ∫ Q26

h

D66 =

( 1)

( 2)

2
∫ Q22 z dz + ∫ Q22

h

D26 =

0
h
4

4

∫h Q16
h

D22 =

2

4

z dz + ∫ Q11
2

4

h

D16 =

( 1)

h

Q66

( 4)

( 4)

z 2dz
z 2dz

2

Công thức tính Qij tại từng lớp được tính như sau:
Lớp 1 và lớp 4 đối xứng nhau (góc xoay θ1 =
Q11

π
( 1)
( 4)
) nên Qij = Qij . Cụ thể:
4

( 1)

= Q11cos 4θ1 + Q22 sin 4 θ1 + 2 ( Q12 + 2Q66 ) sin 2 θ1cos 2θ1

( 1)

= ( Q11 + Q22 − 4Q66 ) sin 2 θ1cos 2θ1 + Q12 sin 4 θ1 + cos 4θ1

( 1)

= ( Q11 − Q12 − 2Q66 ) sin θ1cos3θ1 + ( Q12 − Q22 + 2Q66 ) sin 3 θ1cosθ1

( 1)

= Q11 sin 4 θ1 + Q22cos 4θ1 + 2 ( Q12 + 2Q66 ) sin 2 θ1cos 2θ1

( 1)

= ( Q11 − Q12 − 2Q66 ) sin 3 θ1cosθ1 + ( Q12 − Q22 + 2Q66 ) sin θ1cos3θ1

( 1)

= ( Q11 + Q22 − 2Q12 − 2Q66 ) sin 2 θ1cos 2θ1 + Q66 sin 4 θ1 + cos 4θ1

Q12
Q16
Q22
Q26
Q66

(

(

21

)

)


Lớp 2 và lớp 3 đối xứng nhau (góc xoay θ 2 = −
Q11

π
( 2)
( 3)
) nên Qij = Qij . Cụ thể:
4

( 2)

= Q11cos 4θ 2 + Q22 sin 4 θ 2 + 2 ( Q12 + 2Q66 ) sin 2 θ 2cos 2θ 2

( 2)

= ( Q11 + Q22 − 4Q66 ) sin 2 θ 2cos 2θ 2 + Q12 sin 4 θ 2 + cos 4θ 2

( 2)

= ( Q11 − Q12 − 2Q66 ) sin θ 2cos3θ 2 + ( Q12 − Q22 + 2Q66 ) sin 3 θ 2cosθ 2

( 2)

= Q11 sin 4 θ 2 + Q22cos 4θ 2 + 2 ( Q12 + 2Q66 ) sin 2 θ 2cos 2θ 2

( 2)

= ( Q11 − Q12 − 2Q66 ) sin 3 θ 2cosθ 2 + ( Q12 − Q22 + 2Q66 ) sin θ 2cos3θ 2

( 2)

= ( Q11 + Q22 − 2Q12 − 2Q66 ) sin 2 θ 2cos 2θ 2 + Q66 sin 4 θ 2 + cos 4θ 2

Q12
Q16
Q22
Q26
Q66

(

(

)

)

trong đó Qij được xác định theo công thức (1.5).

M 
 x
Ký hiệu : { Σ} =  M y  . Từ (2.6) ta được: { Σ} = [ D ] { φ }


 M xy 

(2.7)

Bước 4: Thiết lập thế năng toàn phần, ma trận độ cứng, véctơ tải của phần
tử.
Thế năng biến dạng của phần tử tấm:
U etam =

1
ε xσ x + ε yσ y + γ xyτ xy dV
2 ∫∫∫
V

(

)

e

Thay (2.3) vào hệ thức trên ta được :
U etam =

1
zφ xσ x + zφ yσ y + φxyτ xy dxdydz

2 ∫∫∫
V

(

)

e

h
h
h


2
2
2
1
= ∫∫ φx dxdy ∫ σ x zdz + ∫∫ φ y dxdy ∫ σ y zdz + ∫∫ φxy dxdy ∫ τ xy zdz ÷
÷
2  F
÷
F
F
−h
−h
−h

2
2

2


22


Theo (2.5) ta được :
U etam =

T
1
1
M
φ
+
M
φ
+
M
φ
d
xd
y
=
φ
{
}
x x
y y
xy xy

∫∫ { Σ} dxdy
2 ∫∫
2
F
F

(

)

Từ (2.4) và (2.7) suy ra
U etam =

T

1
1 T
T
φ
D
φ
d
xd
y
{
}
{
}
[
]

=
q
B
D
B
d
xd
y
{
}

÷
{ q} e
[
]
[
]
[
]
e  ∫∫
÷
2 ∫∫
2
F
F


U etam =

1

{ q} Te [ K ] tam
{ q} e
e
2

trong đó [ K ] e

tam



T
=  ∫∫ [ B ] [ D ] [ B ] dxdy ÷ được gọi là ma trận độ cứng phần tử

÷
F


Thế năng của nền đàn hồi với hệ số nền K là :
1
Kw2 ( x, y ) dxdy
∫∫
2F
T
1
= K ∫∫  w ( x, y )   w ( x, y )  dxdy
2 F
1
T
T

= K ∫∫ { q} e [ N1 ] [ N1 ] { q} e dxdy
2 F

1
T
T
= K { q} e  ∫∫ [ N1 ] [ N1 ] dxdy ÷{ q} e

÷
2
F


1 T
T
= { q} e  K ∫∫ [ N1 ] [ N1 ] dxdy ÷{ q} e

÷
2
 F

1 T
nen
= { q} e [ K ] e { q} e
2

U enen =
U enen
U enen
U enen

U enen
U enen
với

[ K ] enen = K ∫∫ [ N1 ] T [ N1 ] dxdy
F

Vậy thế năng của tấm composite lớp trên nền đàn hồi sẽ là:

23


U e = U etam + U enen =
với [ K ] e = [ K ] e

tam

+ [ K]e

nen

1 T
{ q} e [ K ] e { q} e
2

được gọi là ma trận độ cứng phần tử tấm composite

lớp trên nền đàn hồi.
Công của ngoại lực do q gây ra trên độ võng của phần tử e là:
Ae = ∫∫ ( q ( x, y ) dF ) w ( x, y )

F

Ae = ∫∫ q ( x, y )  N1 ( x, y )  { q} e dxdy
F

Ae = { q} e

T

Ae = { q}
trong đó

T
e

T

∫∫ q ( x, y )  N1 ( x, y )  dxdy
F

{ P} e
T

{ P} e = ∫∫  N1 ( x, y )  q ( x, y ) dxdy
F

là véctơ tải phần tử của tấm

composite trên nền đàn hồi.
Bước 5: Ghép nối phần tử, tìm ma trận độ cứng tổng thể và véctơ tải tổng

thể.
Giả sử { q} là véctơ chuyển vị nút tổng thể , còn { q} e (e = 1...Ne) là chuyển
vị nút của một phần tử tương ứng. Vậy giữa hai đại lượng này phải có liên hệ:

{ q} e = [ L ] e { q} , trong đó [ L ] e

gọi là ma trận định vị của phần tử có các phần tử

cho bởi quy tắc sau:
0
Lij = 
1

neáu qei ≠ q j
neáu qei = q j

Thế năng toàn phần của tấm composite lớp trên nền đàn hồi là:
Ne

T
1 T

∏e = ∑  { q} e [ K ] e { q} e − { q} e { P} e ÷

e=1  2

24


Ne

T
T
T
1 T

∏e = ∑  { q} [ L ] e [ K ] e [ L ] e { q} − { q} [ L ] e { P} e ÷

e=1  2
1 T
T
∏e = { q} [ K ] { q} − { q} { P}
2

Ne

Trong đó :  K  = ∑

e=1

T
 L  e  K  e  L  e ,

Ne

T

{ P} = ∑  L  e  P  e
e=1

tương ứng là ma trận độ cứng tổng thể, véctơ tải tổng thể của tấm composite lớp
trên nền đàn hồi.
Bước 6: Áp điều kiện biên tìm chuyển vị nút chưa biết
Điều kiện cực trị thế năng toàn phần là
∂Π
= 0 ∀i = 1, N e
∂qi
∂ 1 T
T

 { q} [ K ] { q} − { q} { P} ÷ = 0 ∀i = 1, N e
∂qi  2

Hệ trên trở thành :

[ K ] { q} = { P}

(2.8)

Đây chính là hệ phương trình xác định chuyển vị nút.
Áp điều kiện biên :
• Tấm khớp bản lề tại bốn cạnh:
 Trên biên có: W = 0
 Trên cạnh y = 0, y = b thì θx = 0

(2.9)

 Trên cạnh x = 0, x = a thì θy = 0
• Tấm có hai cạnh ngàm, hai cạnh khớp bản lề
 Trên biên có W = 0

 Tại hai cạnh ngàm: θx = θy = 0
 Tại hai cạnh khớp bản lề: θx = 0 hoặc θy = 0

25

(2.10)