Tạp chí Khoa học Cơng nghệ Xây dựng, NUCE 2021 15 (3V): 1–15 PHÂN TÍCH ĐẶC TRƯNG DAO ĐỘNG CỦA DẦM FGM THEO LÝ THUYẾT DẦM EULER-BERNOULLI BẰNG TIẾP CẬN GIẢI TÍCH Nguyễn Văn Longa,∗, Trần Minh Túa , Trần Hữu Quốca a Khoa Xây dựng Dân dụng Công nghiệp, Trường Đại học Xây dựng, 55 đường Giải Phóng, quận Hai Bà Trưng, Hà Nội, Việt Nam Nhận ngày 12/05/2021, Sửa xong 26/06/2021, Chấp nhận đăng 28/06/2021 Tóm tắt Bài báo phân tích dao động riêng đáp ứng động dầm vật liệu có tính biến thiên (FGM) tác dụng tải trọng phân bố vng góc với bề mặt dầm Với cách chọn hệ tọa độ quy chiếu qua mặt trung hòa, hệ thức quan hệ phương trình chuyển động dầm FGM thiết lập sở lý thuyết dầm Euler–Bernoulli Lời giải giải tích dạng hiển xây dựng cho trường hợp dầm FGM có liên kết hai đầu khớp Ví dụ kiểm chứng thực qua so sánh với công bố tác giả khác sử dụng hệ quy chiếu gắn với mặt trung bình Ảnh hưởng tham số vật liệu, kích thước hình học, tải trọng cưỡng lên đặc trưng dao động dầm khảo sát cụ thể qua ví dụ số Từ khố: phân tích dao động; dầm FGM; mặt trung hịa; lý thuyết dầm Euler-Bernoulli VIBRATIONAL CHARACTERISTIC OF FGM BEAM BASED ON EULER-BERNOULLI BEAM THEORY BY USING ANALYTICAL APPROACH Abstract In this paper, free vibration and transient analysis of FGM beam under transverse distributed load is presented The constitutive relations and governing equations are obtained based on Euler-Bernoulli beam theory including the neutral surface position concept The analytical direct solution is proposed for simply supported FGM beam The validated examples have been conducted by comparison with those of other authors using a reference coordinate system coinciding with the middle surface The effects of material, geometric parameters, excitation on vibrational characteristics of the FGM beams are investigated through numerical examples Keywords: vibration analysis; functionally graded beam; neutral surface position; Euler-Bernoulli beam https://doi.org/10.31814/stce.nuce2021-15(3V)-01 © 2021 Trường Đại học Xây dựng (NUCE) Giới thiệu Trong khoa học vật liệu, vật liệu có tính biến thiên (FGM) loại composite hệ có tính biến thiên liên tục dọc theo phương chiều dày kết cấu Loại vật liệu coi loại vật liệu composite không đồng tạo thành từ hỗn hợp thành phần gốm (ceramic) kim loại (metal) [1], tỷ phần thể tích vật liệu thành phần thay đổi trơn theo phương chiều dày Điều góp phần tránh tập trung ứng suất gây gián đoạn pha vật liệu vật liệu composite lớp hay composite cốt sợi Do tận dụng đặc tính kháng nhiệt kháng ăn mịn gốm, kết hợp với độ bền dẻo kim loại, vật liệu FGM sử dụng rộng rãi nhiều ứng dụng kỹ thuật, đặc biệt cho kết cấu môi trường nhiệt độ cao hàng không vũ trụ, điện hạt nhân, khí, giao thơng vận tải ∗ Tác giả đại diện Địa e-mail: longnv@nuce.edu.vn (Long, N V.) Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Các phân tích ứng xử tĩnh động kết cấu dầm, vật liệu FGM tập trung nghiên cứu năm gần Nhiều cơng trình nghiên cứu đặc trưng dao động ứng xử uốn kết cấu dầm FGM thực Trong số này, Phuong cs [2] xây dựng nghiệm Navier phân tích ứng xử uốn dầm FGM có lỗ rỗng vi mơ sử dụng mơ hình dầm Timoshenko Aydogdu Taskin [3] thiết lập lời giải giải tích dạng nghiệm Navier phân tích dao động tự dầm đơn giản FGM sử dụng lý thuyết biến dạng cắt bậc cao lý thuyết dầm cổ điển Li [4] đề xuất cách tiếp cận hợp nhất, xây dựng lời giải giải tích phân tích uốn dao động riêng dầm FGM sử dụng hai mơ hình dầm Timoshenko Euler-Bernoulli Sina cs [5] phân tích dao động riêng dầm FGM với điều kiện biên khác nhau, sử dụng lời giải giải tích thiết lp t mụ hỡnh tm bc nht Sáimsáek v Kocatăurk [6] phân tích dao động tự đáp ứng động dầm FGM liên kết hai đầu khớp tác dụng tải trọng tập trung di động điều hòa, dựa lý thuyết dầm Euler–Bernoulli khai triển nghiệm chuyển vị dạng đa thức Sử dụng lý thuyết dầm cổ điển, dầm bậc lý thuyết dầm bậc cao, S¸ims¸ek [7] phân tích dao động tự dầm FGM với điều kiện biên khác nhau; đó, thành phần chuyển vị khai triển theo chuỗi đa thức Cũng sử dụng khai triển thành phần chuyển vị dạng chuỗi đa thức, S¸ims¸ek [8] phân tích dao động dầm FGM tác dụng tải trọng di động theo số lý thuyết dầm bao gồm: lý thuyết dầm cổ điển, bậc bậc ba Thai Vo [9] thiết lập nghiệm Navier cho tốn phân tích uốn dao động riêng dầm FGM liên kết hai đầu khớp Áp dụng phương pháp Rayleigh–Ritz, Pradhan Chakraverty [10] phân tích tần số dao động riêng dầm FGM với liên kết khác theo lý thuyết dầm cổ điển dầm bậc Mashat cs [11] phân tích dao động tự dầm FGM với điều kiện biên khác sử dụng lý thuyết biến dạng cắt mơ hình phần tử hữu hạn chiều Su Banerjee [12] áp dụng thuật toán Wittrick–Williams phát triển phương pháp độ cứng động lực phân tích dao động riêng dầm FGM theo lý thuyết dầm Timoshenko với liên kết khác Le cs [13] phân tích đáp ứng động dầm FGM nhiều nhịp tác dụng tải trọng di động phương pháp phần tử hữu hạn sử dụng lý thuyết dầm bậc Hầu hết nghiên cứu tập trung phân tích tần số dao động riêng đáp ứng động dầm, với cách tiếp cận thông thường - hệ tọa độ quy chiếu chọn trùng với mặt phẳng trung bình Với vật liệu FGM, đặc trưng học vật liệu thay đổi theo chiều dày kết cấu, mặt trung bình hình học thường khơng trùng với mặt trung hòa, tồn tương tác màng-uốn quan hệ nội lực - biến dạng [14] Trong nghiên cứu FGM, Zhang [15, 16] tương tác loại bỏ mặt phẳng tham chiếu lựa chọn qua mặt trung hịa Ý tưởng sau số tác giả khác áp dụng cho dầm FGM, chẳng hạn Yaghoobi Fereidoon [17], Larbi cs [18], Li cs [19], Zhang [20] Có thể thấy rằng, với công bố theo cách tiếp cận này, với dầm FGM, theo hiểu biết tác giả, chưa có nghiên cứu đưa biểu thức hiển để khảo sát đáp ứng động Để góp phần làm phong phú thêm nghiên cứu lý thuyết ứng xử học dầm FGM, báo đưa cách tiếp cận giải tích nhằm phân tích dao động riêng đáp ứng động dầm FGM Với mục đích đưa lời giải dạng hiển phục vụ cho tính tốn, thiết kế sơ bộ, lý thuyết dầm Euler-Bernoulli với hệ tọa độ quy chiếu qua mặt trung hòa sử dụng Các phương trình chủ đạo dầm thiết lập sở nguyên lý Hamilton, có kể đến vị trí mặt trung hịa nên có dạng đơn giản Tuy nhiên để thu nghiệm giải tích dạng hiển tần số dao động riêng đáp ứng động qua việc giải phương trình vi phân dao động cưỡng bức, dầm FGM với liên kết khớp hai đầu, môi trường không cản đối tượng chọn Các kết số thực nhằm đánh giá ảnh hưởng tham số vật liệu, tham số hình học, tải trọng cưỡng đến tần số dao động riêng đáp ứng động dầm FGM Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Mơ hình hóa dầm vật liệu FGM Vật liệu FGM cấu thành từ hai vật liệu thành phần gốm kim loại Các đặc trưng học giả thiết biến đổi trơn dọc chiều cao dầm Tỷ lệ thể tích thành phần vật liệu thường biến đổi theo quy luật hàm lũy thừa (P-FGM), hàm mũ (E-FGM) hàm Sigmoid (S-FGM) [21] Trong vật liệu E-FGM thường sử dụng phân tích phá hủy, vật liệu P-FGM thường sử dụng phân tích tĩnh, dao động ổn định Vật liệu S-FGM thường sử dụng với kết cấu FGM nhiều lớp để tránh tập trung ứng suất bề mặt phân cách lớp Do vậy, nghiên cứu này, dầm vật liệu P-FGM, với chiều dài L, mặt ngang chữ nhật với bề rộng b, chiều cao h Hình xem xét ặ trung bình mặt ặ trung hịa củaủ dầmầ FGM Hình 1.ịVị trí mặt ấ ệ ụ ủ ậ ệ ô đun đàn hồ ối lượ Các tính chất hiệu dụng vật liệu FGM bao gồm mô đun đàn hồi E, khối lượng riêng ρ, ρ đượ ể ễ dướ biểu diễn dạng [18]: p zth + C ) (1) P(zth) = Pm + (P − P + c m h với Pm Pc tương ứng tính gốm kim loại; p ≥ số tỷ lệ thể tích; C khoảng tương ứng tính củ ố ≥ ố ỷ ệ ể cách mặt trung hòa mặt trung bình Trong nghiên cứu này, để đơn giản tính tốn, hệ số ả ặ tọa độ chiều cao ặ dầm, ν = const ứu này, để Poisson coi không thay ữđổi theo ệ ốFGM Poisson coi tọa khơng thay độ hịa ề mặt Dođơn mô giả đun đàn hồi E vật liệu thayđược đổi theo độ chiều caođổdầm, mặtọatrung ầ ν trung bình hình học dầm thường khơng trùng (khi p 0) Trên Hình 1, trục tọa độ x, z tương ứng với phương trục chiều cao dầm, tọa độ điểm hệ tọa độ qua mặt Do mô đun đàn hồ ủ ậ ệ FGM thay đổ ọa độ ề ầ ặ trung bình mặt trung hịa biểu diễn tương ứng ztb zth Khi đó, vị trí mặt trung hòa ặ ọ ủ ầm thườ ≠ 0) dầm FGM xác định từ điều kiện [18, 22]: ụ ọa độ tương ứ ới phương ụ ề ầ ọa độ h/2 điể ấ ỳ ệ ọa độ qua mặ ặt trung hòa đượ ể ễ ztb E(ztb )dztb h/2 tương ứ ặ−h/2 ủ ầm FGM xác đị z z Khi đó, vị h (Eđược c − Em ) p (ztb − C) E(ztb )dztb = ⇒ C = = (2) điề ệ h/2 (p + 2) (pEm + Ec ) −h/2 E(ztb )dztb −h/2 Các phương trình lý thuyết dầm Euler-Bernoulli Sử dụng hệ quy chiếu gắn với mặt trung hòa, trường chuyển vị theo lý thuyết dầm Euler-Bernoulli [23]: ∂w0 (x, t) (3a) u(x, zth , t) = −zth phương trình bả ủ ế ầ ∂x ụ ệ ế ắ ặt trung hòa, trườ , ) ể ị ế ầ Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng (3b) w(x, zth , t) = w0 (x, t) w0 độ võng điểm mặt trung hòa theo phương trục zth ; t biến thời gian Các thành phần biến dạng xác định thông qua thành phần chuyển vị: ∂u ∂2 w0 = −zth ∂x ∂x εx = γ xzth = (4a) ∂w ∂u ∂w0 ∂w0 + − =0 = ∂x ∂zth ∂x ∂x (4b) Như vậy, lý thuyết dầm Euler-Bernoulli, biến dạng góc: γ xzth = Ứng suất pháp liên hệ với biến dạng theo định luật Hooke: σ x = E(zth )ε x = −E(zth )zth ∂2 w0 ∂x2 (5) Với lý thuyết dầm Euler-Bernoulli, ta định nghĩa thành phần nội lực dầm bao gồm: h/2−C σ x dA = b Nx = A σ x dzth (6a) zth σ x dzth (6b) −h/2−C h/2−C zth σ x dA = b Mx = A −h/2−C Thay (5) vào (6), ta được: h/2−C N x = −b ∂2 w ∂2 w0 E(zth )zth dzth = − b ∂x ∂x M x = −b zth E(zth )dzth = −B11 ∂2 w0 ∂x2 (7a) z2th E(zth )dzth = −D11 ∂2 w0 ∂x2 (7b) −h/2−C −h/2−C h/2−C h/2−C ∂2 w ∂2 w0 z2th E(zth ) dzth = − b ∂x ∂x −h/2−C h/2−C −h/2−C Các số độ cứng dầm (7) xác định bởi: h/2−C zth E(zth )dzth = B11 = b (8a) −h/2−C h/2−C z2th E(zth )dzth D11 = b (8b) −h/2−C Có thể thấy rằng, với việc sử dụng hệ tọa độ qua mặt trung hòa, B11 = 0; dẫn đến lực dọc N x = 0; biểu thức xác định thành phần nội lực dầm FGM trở nên đơn giản, tương tự với dầm đẳng hướng Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Phương trình chuyển động cho dầm xây dựng dựa nguyên lý Hamilton [24, 25]: T (9) (δU + δV − δK) dt 0= δU biến phân biến dạng đàn hồi dầm, δV biến phân tải trọng δK biến phân động Biến phân biến dạng đàn hồi dầm: L L σ x δε x dAdx = − δU = A 0 L ∂2 δw0 zth σ x dAdx = − ∂x2 A Mx ∂2 δw0 dx ∂x2 (10) Biến phân tải trọng uốn: L δV = − (11) qδw0 dx Biến phân động xác định bởi: L L ρ(zns ) (˙uδ˙u + wδ ˙ w) ˙ dAdx = δK = A I0 w˙ δw˙ + I2 ∂w˙ ∂δw˙ dx ∂x ∂x (12) Dấu (.) biến chuyển vị thể đạo hàm theo thời gian t Các mơ men qn tính khối lượng định nghĩa (12) tính theo: h/2−C ρ(zth )dA = b I0 = A ρ(zth )dzth (13a) z2th ρ(zth )dzth (13b) −h/2−C h/2−C z2th ρ(zth )dA I2 = A =b −h/2−C Thay biểu thức δU, δV δK từ (10), (11) (12) vào (9), tiến hành tích phân phần, ta được: = V xz δw0 L ∂δw0 − Mx ∂x L L − ∂2 M x wă ău0 + I + q I w ă I w0 dx 0 x x2 x2 (14) wă M x I0 uă + I2 l lc ct hiệu dụng ∂x ∂x Cho hệ số biến phân độ võng chiều dài dầm không, phương trình chuyển động dầm thu được: ∂2 M x wă + q = I w ă I (15) 0 ∂x2 ∂x2 đó: V xz = Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Các tham số điều kiện biên: chuyển vị, lực rút từ (14): ∂w0 , Mx ∂x (w0 , V xz ) ; (16) Thay (7) vào (15), ta phng trỡnh chuyn ng theo chuyn v: D11 wă w q + I w ă − I =0 0 ∂x4 ∂x2 (17) Đây phương trình chủ đạo để phân tích đáp ứng động dầm FGM theo lý thuyết dầm EulerBernoulli Tiếp theo, báo tập trung phân tích dao động riêng dao động cưỡng dầm FGM với liên kết hai đầu khớp (SS) Các biểu thức điều kiện biên dầm bao gồm x = x = L: (18a) w0 = M x = hay ∂2 w0 =0 ∂x2 (18b) Phân tích dao động riêng Trong phân tích dao động riêng, bỏ qua tải trọng uốn: q = 0, phương trình chuyn ng ca dm FGM cú dng: wă w (19) D11 + I0 wă − I2 = ∂x ∂x Nghiệm độ võng giả thiết dạng: w0 (x, t) = W(x) sin (ωt + ϕ) (20) ω tần số dao động riêng dầm, ϕ pha ban đầu Thay (20) vào (19), ta được: D11 d4 W 2d W + I ω − I0 ω W = dx4 dx2 (21) Đây phương trình vi phân cấp tuyến tính nhất; phương trình đặc trưng: D11 k4 + I2 ω2 k2 − I0 ω2 = (22) Bằng cách đặt s = k2 , phương trình (22) đưa dạng: (23) D11 s2 + I2 ω2 s − I0 ω2 = Phương trình bậc hai (23) có hai nghiệm trái dấu: −I2 ω2 + s1 = I22 ω4 + 4D11 I0 ω2 2D11 = λ2 > (24a) Long, N V., cs / Tạp chí Khoa học Công nghệ Xây dựng −I2 ω2 − I22 ω4 + 4D11 I0 ω2 s2 = 2D11 (24b) = −µ2 < Do đó, nghiệm phương trình (22): −I2 ω2 + k1,2 = ± I22 ω4 + 4D11 I0 ω2 2D11 I2 ω2 + k3,4 = ±i I22 ω4 + 4D11 I0 ω2 2D11 = ±λ (25a) = ±iµ (25b) Từ đó, suy nghiệm phương trình (21) có dạng: W(x) = C1 sin µx + C2 cos µx + C3 sinh λx + C4 cosh λx (26) Các số C1 , C2 , C3 , C4 phụ thuộc vào điều kiện biên (18); ba bốn số tần số dao động riêng ω xác định thơng qua bốn phương trình điều kiện biên Điều kiện để dầm dao động tự số C1 , C2 , C3 , C4 không đồng thời không Với trường hợp dầm liên kết hai đầu khớp (SS); bốn phương trình điều kiện biên cho kết quả: C2 = C3 = C4 = (27a) C1 sin µL = (27b) Để số C1 , C2 , C3 , C4 không đồng thời khơng, suy C1 sin µL = ⇒ µL = mπ; 0; từ suy ra: m = 1, 2, 3, (28) Từ đây, ta thu tần số dao động riêng: ωm = m2 π2 L2 D11 / I0 + I2 m2 π2 L2 (29) Tần số dao động tần số nhỏ số tần số dao động riêng thu được; rõ ràng với dầm liên kết hai đầu khớp, m = 1: ωcb = {ωm } = ω1 = π2 L2 D11 / I0 + I2 π2 L2 (30) Hàm độ võng dầm đó: M Wm sin W(x) = m=1 mπx L M ⇒ Wm sin w0 (x, t) = m=1 mπx sin (ωm t + ϕm ) L (31) Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Phân tích dao động cưỡng Xét dầm FGM liên kết hai đầu khớp (SS), tác dụng tải trọng phân bố q(t) = Q sin Ωt vng góc với trục dầm Phương trình chuyển động dm (17) c vit li thnh: D11 wă w0 + I w ă I = Q sin Ωt 0 ∂x4 ∂x2 (32) Nghiệm tổng quát phương trình vi phân (32) bao gồm: nghiệm phương trình (khi khơng có vế phải) nghiệm riêng phương trình có vế phải Nghiệm phương trình phương trình (32) phương trình vi phân dao động tự do, có dạng: M mπx h sin (ωm t + ϕm ) (33) Wmh sin w0 (x, t) = L m=1 đó: Wmh ϕm tương ứng biên độ pha ban đầu dao động tự do; chúng phụ thuộc vào điều kiện ban đầu toán Để tìm nghiệm riêng phương trình (32), lưu ý rằng, tải trọng phân bố Q khai triển theo biến không gian x dạng chuỗi lượng giác tương tự độ võng: M Qm sin Q(x) = m=1 mπx L (34) mπx dx L (35) Với tải trọng bất kỳ: L Qm = L Q(x) sin Đối với cho tải trọng phân bố Q = Q0 , hệ số Qm xác định bởi: m = 2, 4, 6, 2Q0 Qm = (1 − cos mπ) ⇒ Qm = 4Q mπ m = 1, 3, 5, Qm = mπ (36) Thay (34) vào (32) ta tìm nghiệm riêng phương trình là: M p w0 = mπx Fm sin sin Ωt L − Ω2 ω2 m=1 m (37) Qm m π2 I0 + I2 L Từ đó, ta có nghiệm tổng quát độ võng: đó: Fm = M p Wmh sin (ωm t + ϕm ) + w0 = wh0 + w0 = m=1 mπx Fm sin Ωt sin L − Ω2 ω2m (38) Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Trong báo này, tác giả sử dụng điều kiện ban đầu toán là, t = 0: w0 = w˙ = 0, ∀x ∈ [0, L] Do đó, ta thu được: ϕ =0 m Ω/ωm ∀m (39) h Fm Wm = 2 Ω − ωm Thay kết thu (39) vào (38) ta được: M Ω/ωm sin mπx sin ωm t + w0 (x, t) = sin Ωt Fm L ωm − Ω Ω2 − ω2m m=1,3,5, (40) Phương trình (40) sử dụng để phân tích đáp ứng động dầm FGM liên kết hai đầu khớp Từ đây, ta tìm độ võng dầm (x = L/2): M Ω/ωm mπ sin (41) sin ωm t + Fm w0 (L/2, t) = sin Ωt 2 2 ωm − Ω Ω − ωm m=1,3,5, Kết số thảo luận Xét dầm vật liệu P-FGM (Al/Al2 O3 ), liên kết hai đầu khớp; chiều cao tiết diện dầm h = 0,1 m, bề rộng tiết diện dầm b = 0,05 m, chiều dài dầm L Mô đun đàn hồi, hệ số Poisson, khối lượng riêng kim loại (Al): Em = 70 GPa, νm = 0,3, ρm = 2702 kg/m3 ; gốm (Al2 O3 ): Ec = 380 GPa, νm = 0.3, ρm = 3800 kg/m3 6.1 Ví dụ kiểm chứng Xét dầm tiết diện chữ nhật, liên kết hai đầu khớp, vật liệu P-FGM (Al/Al2 O3 ) Bảng trình bày L2 ρm /Em cho dầm hai kết tính tốn tần số dao động riêng không thứ nguyên, ω ¯ i = ωi h Bảng Kiểm chứng tần số dao động riêng không thứ nguyên, ω ¯ dầm FGM liên kết hai đầu khớp L/h 20 m Nguồn p 0,5 10 Thai Vo [9] Bài báo 5,3953 5,3953 4,5931 4,5932 4,1484 4,1485 3,7793 3,7796 3,4921 3,4923 Thai Vo [9] Bài báo 20,6187 20,6187 17,5415 17,5452 15,7982 15,8063 14,3260 14,3400 13,2376 13,2521 Thai Vo [9] Bài báo 43,3483 43,3483 36,8308 36,8640 33,0278 33,0996 29,7458 29,8650 27,4752 27,6041 Thai Vo [9] Bài báo 5,4777 5,4777 4,6641 4,6641 4,2163 4,2163 3,8472 3,8472 3,5547 3,5547 Thai Vo [9] Bài báo 21,8438 21,8438 18,5987 18,5987 16,8100 16,8100 15,3334 15,3335 14,1676 14,1677 Thai Vo [9] Bài báo 48,8999 48,8999 41,6328 41,6330 37,6173 37,6178 34,2954 34,2961 31,6883 31,6891 Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng trường hợp kích thước dầm L/h = 5; 20 áp dụng cho dạng dao động với số tỷ lệ thể tích p = 0; 0,5; 1; 2; 10 Nghiệm giải tích báo so sánh với kết tác giả Thai Vo [9] sử dụng lý thuyết dầm Euler-Bernoulli (CBT) với dạng nghiệm Navier; nhiên tính tốn hệ trục gắn với mặt trung bình Từ kết tính tốn kiểm chứng trên, thấy lời giải giải tích chương trình máy tính sử dụng báo có độ tin cậy (sai lệch lớn L/h = 5, p = 10 m = 0,4691%) 6.2 Khảo sát tần số dao động riêng dầm FGM Xét dầm chữ nhật, liên kết hai đầu khớp, vật liệu P-FGM dao động tự môi trường hồn tồn lý tưởng (khơng cản) Hình thể dạng dao động (biến thiên độ võng w0 theo chiều dài dầm) dầm FGM với dạng dao động với, tương ứng với m = 1; 2; Có thể thấy rằng, giá trị m số nửa bước sóng hình sin đường đàn hồi; m tăng, tần số dao động riêng dầm tăng Tần số dao động ứng với trường hợp m = (đường đàn hồi có dạng nửa bước sóng hình sin) có trị số nhỏ nhất, nhỏ nhiều so với dạng dao động ủ dầm FGM ωbả Hình 3.ế Biến thiên ế ầ tầnốsố ầcơbảbản, ố ủ cb theo p L/h ầ Hình Ba dạngạng daodao động riêngđộ dầm đầ củađầ ạngđộ dao FGM tương ứng với m = 1; 2; ủ ầ tương ứ ủ ầ tương ứ ầ số dao ủ FGM ầ ủ với ỷ ỉtích, ệ ố ểỷp ỷ thước ố ỷdầm bảndầm TầnTần số dao động ωcb Tần số động dao động cơcủa ầ cácớchỉ sốớtỷỉ lệốthể ệ tỷểsố kích ố Biến thiên theo hai L/hkích tính tốn Bảngbày thướ ầ trình đượcbày tínhnhư tốn trình Bảngtần số dao ế động ế ủ ầ ủ ầ tham sốkích thướ thểầhiện Hìnhtính tốn trình bày Bảng ố dao động bả ố đượ ể ệ dao động bả p tăng (đồngốnghĩa đượ ệ hàm lượng kim loại vật liệu Các ốkết chocơ thấy, với việcể tăng Bảng ố dao ố cấu ỷ ệgiảm) ể Với FGM), khiầ L/h tăngđộng (dầmcơ dàibảra), tần số dao cơủ bảnầ ωcb giảm ớ(độ cứngỉ kết ủ ầ ỉ ố ỷ ệ ể Bảng ầ ố dao động bả tỷ số kích thước dầm L/h, tần số ω giảm nhanh p tăng khoảng [0, 2], sau ỷ ố kíchcbthướ ầ ỷ ố kích thướ ầ giảm chậm lại gần không thay đổi p tiến tới 10 Đồng thời, với số tỷ lệ thể tích p, tần số ωcb giảm nhanh L/h tăng khoảng [10, 20], sau giảm chậm lại dầm trở nên dài (L/h tiến tới 50) 6.3 Khảo sát ảnh hưởng tham số vật liệu, kích thước hình học, tải trọng uốn lên đáp ứng dao động dầm FGM Trong phần này, báo xem xét dầm chữ nhật, liên kết hai đầu khớp, vật liệu P-FGM (Al/Al2 O3 1); dầm dao động môi trường không cản tác dụng tải trọng cưỡng phân bố q(t) = Q0 sin Ωt 10 ế ậ ả ệ ế ấ ả ặ ấ tăng (đồng nghĩa vớ ệc tăng hàm lượ ệc tăngtăng (dầ (đồng nghĩa ầ vớ ố cơtăng bả hàm lượ ảm (độ Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Bảng Tần số dao động bản, ωcb [rad/s] dầm FGM với số tỷ lệ thể tích, p tỷ số kích thước dầm L/h khác p L/h 0,5 10 10 20 30 40 50 2837,5 2403,0 2165,1 1967,9 1905,8 1865,8 1830,4 1807,5 711,5 602,6 543,1 493,8 478,2 468,2 459,3 453,5 316,4 268,0 241,5 219,6 212,7 208,2 204,3 201,7 178,0 ầ 150,8 135,9 123,6 119,7 117,2 114,9 113,5 ầ ữ ậ ết hai đầ ậ ầ ầ ữ ậ ết hai đầ ậ ệ ệ 113,9 96,5 ầm dao 87,0 79,1 76,6 75,0 động ản dướ ụ73,6 ụủ ảủ72,6ả ầm dao động môimôi trườtrườ ản dướ ọng cưỡ ứ ứ ọng cưỡ ố đềố đề Khảo sát ảnh hưởng số số hạng khai triển chuỗi lượng giác, M Khảo hưởng số hạng số hạng chuỗi lượng Khảo sátsát ảnhảnh hưởng củacủa số số khaikhai triểntriển chuỗi lượng giácgiác Xét dầm FGM với tỷ số kích thước dầm L/h = 50 số tỷ lệ thể tích p = 2; hai tần số dao ố kích thướ ỷ ểệ ể tác dụng tải ầ ầdầm là: ω ớ= ỷ 79,08 ốỷ kích thướ ỉ ố ỉdao ỷ ốệđộng động riêng rad/s ωầ3 =ầ710,65 rad/s Dầm ầ ố dao độ ủ ầ ω ω ầ phân ố dao độ q(t) = sin ủ Ωtầ (kN/m,ωs) ω ầ ầ trọng bố đều: dướ ụ ụ ủ ủả ảọ ọ độđộ dướ ố đềố đề (b) Ω = ΩΩ== 600 rad/s (b)(b) (a)(a) = 30 30 (a) ΩΩ =Ω rad/s = rad/s 30 rad/s Hình So hưởng sánh ảnhcủa củahạng sốtrong hạng khai triển khai triển chuỗi lượng giác Mlên đáp SoSo sánh ảnh số số chuỗi lượng giác sánh ảnh hưởnghưởng số sốsốhạng khai triển chuỗi lượng giác lên đáp lên đáp ứng động dầm FGM ứngứng động dầm FGM động dầm FGM Hình khảo sát đường cong tải trọng - độ võng dầm FGM với hai trường hợp tần số lực cưỡng ảo ảo sátsát đườđườ ả ảọ ọ độ độ ủ ủầ ầ ới hai ợ ợ ớitrườ haikhai trườtriển khảo sát bao gồm: (a) Ω = 30 rad/s, (b) Ω = 600 rad/s Số số hạng chuỗi ầ ầ ố ốựcực cưỡcưỡ ức ức đượđượ ả ả ồ Ω Ω Ω Ω ố ố lượng giác nghiệm độ võng lấy với giá trị: M = 1, M = 10 M = 100 Các kết đồ ố ạthấy, tần số lựcể cưỡngỗibức lượgần với tần sốệm độ võng đượ ấ dầm, ị khai thị cho daoệm động bảnđượ ố ể ỗi lượ độ võng ấớ sốớ số hạng ị ế cho ả đồ đáp ị ứng động ấ gần ầnhưốtrùng ực cưỡ ứ Tuy ầ nhiên, triển M = 1, M = 10 M = 100 kết khớp ế ả đồ ị ấ ầ ố ực cưỡ ứ ầ số lựcbả cưỡng trường ầ ố hợp dao tần động ủ ầ gầnốvớiố tần số dao động riêng ể dầm, lấy ầ ố dao động bả ủ ầ ố ố ể số hạng trongếkhai ảtriển 1) đáp ứng trùng động bảo độ hội tụ (sai lệch lớn đáp (M ứng= độ ần khớdầm chưa đảm y nhiên, trườ ế khiả đáp ứng độ ần trùng khớ y nhiên, trườ nhấtợvề biên độ độ võng M = so với M = 10 M = 100 lên đến 14,5%) ầ ố ực cưỡ ứ ầ ầ ố dao độ ế ủ ầ ế ỉ ấ ợ ầ ố ực cưỡ ứ ầ ầ ố dao độ ế ủ ầ ế ỉ ấ ộ ố ển (M = 1) đáp ứng độ ủ ầm chưa đả ảo độ ộ Khảo sát số tỷển lệ thể vàđáp tỷ sốứng kíchđộ thước ủ dầm,ầm L/hchưa đả ộ ảnh ố hưởng (M tích, = 1)p,thì ảo độ ộ ụ ệ ấ ề biên độ độ Xét ụ dầm ệFGMớdưới ấtác ềdụng biêncủa độ tải độ trọng cưỡng phân ớbố q(t) = sin 30t (kN/m, s) đế bao gồm đường cong độ võng-thời gian dầm FGM với số tỷ lệ thể tích, p, Hình đế ả ảnh hưở ủ ỉ ố ỷ ệ ể ỷ ố kích thướ ầ ả ảnh hưở ủ ỉ ố ỷ ệ ể 11 ỷ ố kích thướ ầ ầm FGM dướ ụ ủ ả ọng cưỡ ứ ố đề ầm FGM dướ ụ ủ ả ọng cưỡ ứ ố đề ồm đườ độ ủ ầ Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng tỷ số kích thước dầm, L/h khác Các kết đồ thị cho thấy, tăng p L/h làm tăng biên độ độ võng dầm Cụ thể là: p = 10, dầm có biên độ độ võng lớn gấp 4,19 lần so với p = 0; L/h = 50, dầm có biên độ độ võng tăng 59,2 lần so với L/h = 20 Ảnh hưở ệểđápể hưở ủ tỷủlệỉthểốỉ tích ỷố pệỷ lên Hình Ảnh Ảnh hưởng số ứng động dầm FGM lênlên đápđáp ứngứng độ độ ủ ủ ầ ầ ả ả Ảnh hưởng tỷ số kích hưởng tỷkích số kích thước Hình Ảnh Ảnh hưởng củacủa tỷcủa số thước L/hthước lên đáp ứng động dầm FGM động dầm FGM lênlên đápđáp ứngứng động củacủa dầm FGM hưở ủ ủbiên biên ảnhảnh hưở độ độ ực ực cưỡcưỡ ứ ứ Khảo sát ảnh hưởng biên độ lực cưỡng bức, Q Hình khảo sát đáp ứng động dầm FGM (L/h = 50, p = 2, ω1 = 79,08 rad/s) tác dụng tải trọng cưỡng phân bố q(t) = Q sin 30t (kN/m, s); ba trường hợp biên độ lực cưỡng xem xéthưở bao gồm:ủ Q =ỉ0,5; Có thể thấy biên độ số lựckích cưỡng Q Ảnh ố 1ỷ vàệ (kN/m) ể Ảnhrằng, hưởng tỷ thước tăng ảnh hưởng lên biên độ dao động, vận tốc độ võng, chúng tăng tuyến tính theo biên độ lực lên đáp ứng độ ủ ầ lên đáp ứng động dầm FGM cưỡng Đường cong quỹ đạo pha trường hợp tần số lực cưỡng xa so với tần số bản,ảtrở nênảnh hỗnhưở loạn, gồmủ nhiều cong biênđường độ ực cưỡgiao cắt ứ Đáp ứngứng động củacủa dầm FGM táctác dụng củacủa tải tải trọng uốnuốn Đáp động dầm FGM dụng trọng vớivới cáccác biên độ độ củacủa lựclực cưỡng bứcbức khác biên cưỡng khác ảo ảo sátsát đápđáp ứngứng độ độ ủ ủ ầ ầ ω ω dướ cưỡcưỡ ứ ứ ố đề (kN/m, s); s); ba ba trườ dướ ụ ụ ủ ủả ả ọngọng ố đề (kN/m, trườ ợ ợ n độ cưỡcưỡ ức ức đượđượ ồ ể ể n độủ ủực ực ấ ấ ằng, khikhi biên độ độ ực ực cưỡcưỡ ứ ứ tăng ảnhảnh hưởng lênlên biên độ độ daodao độ độ ậ ậố ố ằng, biên tăng hưởng biên Đáp ứng động dầm FGM tác dụng trọng uốn qs) vớiQcác ủa độ võng, chúng tăng tuyến tínhtính theo biên độ ực cưỡ Đườ ỹ đạ võng, chúng tăng tuyến theo biên độ ực ức ỹ đạđộ Hình độ Đáp ứng động dầm FGM tác dụng tải trọng uốncưỡ q = tải Qức sin 30tĐườ (kN/m, biên lực cưỡng khác nhauớ phapha trườ ố ốủcáccủa cưỡ ứlực bả bả ở ỗ ỗ biên độ trườ ợ ợ ầ ầvới ủực ực cưỡcủa ứ cưỡng ớầ khác ầ ố ốnhau ạ ồ ều ều đườ ắ ắ đườ ứng ầ ω s) lên tượng Hình khảo ảo sát sát ảnh đáp hưởng củađộ biên độủlực cưỡng bức: q(t) = Q sin 70t (kN/m, phách FGM (L/h 50,ủa p biên = 2, 79,08 thấy rằng, trường ả ủảdầm hưở độωđộ cưỡcưỡrad/s) ứ ứDễ dàng(kN/m, 1ực ảnh hưở dướ điều hòa ụ ả ảnh ọng cưỡ =ủa ứbiên ố= ực đề s);trong ba trườ tượng điềđiề ức đượ ủ ủ ầ ầ 12 ω ω ễ ễể phách ợ ệnnện độtượng ủ phách ực cưỡ ấ ấkhi ằng, ực ực cưỡ ứ ứbiên ầ ầđộ dao ố bả ằng, ốủ ủảnh cưỡ lên ớầ ầđộ ố ấ ằng, biên độ trườ ựctrườ cưỡ ợ ợứầ ầ ố tăng hưởng ậ bảố đườ ỹ đạ ần ần nhưnhư đườ , dao độỹđộđạ ỹ đạtăng tuyến làtheo đườ độ ực ạcưỡ ắ ức ắố ốĐườ , dao ủađườ độ võng, chúng tính biên ủ ủtrong ầmầm FGM tương daodao điề độ độ ớtăng FGM tương pha trườ ợ tựầtự ốđộng ủđộng ựcđiề cưỡ iên ứ iên ầtăngố ếcơếbả ề ềỗ Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng hợp tần số lực cưỡng gần với tần số bản, đường cong quỹ đạo pha gần đường cong dạng xoắn ốc khép kín, dao động dầm FGM tương tự dao động điều hòa Biên độ phách tăng tuyến tính theo chiều tăng biên độ lực cưỡng bức, chiều dài phách khơng đổi Nhìn chung,ện khitượ tượng phách biênđộđộ độ võng võng tăng kể kể so với trường hợp tầnợsố ảy xảy ra, ra, cáccác biên độ tăngđáng đáng ới trườ lực cưỡng xa tần số dao động dầm ầ ố ực cưỡ ứ ầ ố dao động bả ủ ầ ện tượ ố ực cưỡ ầ ảy ra, biên độ độ võng tăng đáng kể ầ ố dao động bả ủ ầ ứ ới trườ ợ hưở biên biên ực cưỡtới hiệnứtượng phách ện tượng Hình So sánh ảnh ảnh hưởng độ lựcđộ cưỡng điều hòaphách dầmđiề FGM ủ ả ảnh hưở ủ ầ ố ực cưỡ ầ ức, Ω lên tượng phách điề Khảo sát ảnh hưởng tần số lực cưỡng bức, Ω lên tượng phách điều hịa Hình mô tả đáp ứng động dầm FGM: (L/h = 50, p = 2, ω1 = 79,08 rad/s) tác dụng hưở bố biên ựcΩtcưỡ trường ện tượng phách tải trọng cưỡng ảnh phân q(t) =độ sin (kN/m, ứ s); hai hợp tần số củađiề lực cưỡng ủ ầ xem xét bao gồm: Ω = 65 70 (rad/s) Rõ ràng là, tần số lực cưỡng gần vớiảtần số bản, biên độủphách dài pháchức, đềuΩtăng ảnh hưở ầ vàốchiều ực cưỡ lên theo tượng phách điề ả h hưở ủ ầ ố ực cưỡ ủ ầ ứ ện tượng phách điề ả đáp ứng độ ủ ầ ω dướ ụ ủ ả ọng cưỡ ứ ố đề ) = 2sinΩ (kN/m, s); hai trườ ợ ầHìnhố9 So ủ sánh ực ảnh Ωớ phách ảcưỡ hhưởng hưở ức ủ sốầlực ốcưỡng ực cưỡtới ồhiệnứ tượng ện tượng điề đượ tần điều hòaphách dầm FGM ủ ầ ầ ố ủ ực cưỡ ứ ầ ầ ố bản, biên độ ề phách tăng theo ả đáp ứng độ ủ ầ ω dướ ụ ủ ả ọng cưỡ ứ ố đề ) = 2sinΩ (kN/m, s); hai trườ ợ ầ ố ủ ực cưỡ ức đượ Ω 13 Kết luận ầ ố ủ ực cưỡ ứ ầ ầ ố bản, biên độ ề Trên sở theo ế ầ ọ ệ ọa độ ếu qua mặ phách đềucơtăng trung hồn, báo thiế ập phương trình chủ đạo để phân tích đáp ứng độ Long, N V., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Kết luận Trên sở lý thuyết dầm Euler-Bernoulli chọn hệ tọa độ quy chiếu qua mặt trung hoàn, báo thiết lập phương trình chủ đạo để phân tích đáp ứng động dầm FGM chịu tải trọng uốn dạng đơn giản, tương tự dầm đẳng hướng Biểu thức giải tích dạng hiển tần số dao động riêng đáp ứng động dầm FGM liên kết khớp hai đầu thiết lập cho thấy ưu điểm việc sử dụng hệ trục tọa độ qua mặt trung hịa để tính tốn Chương trình tính tự viết Matlab nghiệm giải tích thu chứng tỏ độ tin cậy qua so sánh với kết tác giả khác sử dụng nghiệm Navier tính tốn mặt trung bình Các khảo sát số thực cho phép đánh giá ảnh hưởng tham số hình học, vật liệu đến đáp ứng động dầm FGM liên kết khớp hai đầu Các kết nhận cho thấy: - Tần số dao động dầm ln ứng với dạng dao động gồm nửa bước sóng (m = 1); - Việc lựa chọn lấy đầy đủ số số hạng khai triển nghiệm độ võng cần thiết để khảo sát đáp ứng động; - Khi tăng p, tỷ số kích thước L/h làm giảm độ cứng kết cấu dầm, dẫn đến giảm tần số dao động tự độ võng tăng; - Biên độ lực cưỡng tỷ lệ với độ võng dầm phân tích đáp ứng động tuyến tính; - Khi tần số lực cưỡng gần với tần số dao động tự do, tượng phách điều hòa xảy với dầm Biên độ lực cưỡng tăng làm tăng biên độ phách; tần số lực cưỡng ảnh hưởng tới biên độ chiều dài phách Tài liệu tham khảo [1] Koizumi, M (1993) The concept of FGM [2] Phuong, N T B., Tu, T M., Phuong, H T., Long, N V (2019) Bending analysis of functionally graded beam with porosities resting on elastic foundation based on neutral surface position Journal of Science and Technology in Civil Engineering (STCE) - NUCE, 13(1):33–45 [3] Aydogdu, M., Taskin, V (2007) Free vibration analysis of functionally graded beams with simply supported edges Materials & Design, 28(5):1651–1656 [4] Li, X.-F (2008) A unified approach for analyzing static and dynamic behaviors of functionally graded Timoshenko and Euler–Bernoulli beams Journal of Sound and Vibration, 318(4-5):1210–1229 [5] Sina, S A., Navazi, H M., Haddadpour, H (2009) An analytical method for free vibration analysis of functionally graded beams Materials & Design, 30(3):741747 [6] Sáimsáek, M., Kocatăurk, T (2009) Free and forced vibration of a functionally graded beam subjected to a concentrated moving harmonic load Composite Structures, 90(4):465–473 [7] S¸ims¸ek, M (2010) Fundamental frequency analysis of functionally graded beams by using different higher-order beam theories Nuclear Engineering and Design, 240(4):697–705 [8] S¸ims¸ek, M (2010) Vibration analysis of a functionally graded beam under a moving mass by using different beam theories Composite Structures, 92(4):904–917 [9] Thai, H.-T., Vo, T P (2012) Bending and free vibration of functionally graded beams using various higher-order shear deformation beam theories International Journal of Mechanical Sciences, 62(1):57– 66 [10] Pradhan, K K., Chakraverty, S (2013) Free vibration of Euler and Timoshenko functionally graded beams by Rayleigh–Ritz method Composites Part B: Engineering, 51:175–184 [11] Mashat, D S., Carrera, E., Zenkour, A M., Khateeb, S A A., Filippi, M (2014) Free vibration of FGM layered beams by various theories and finite elements Composites Part B: Engineering, 59:269–278 [12] Su, H., Banerjee, J (2015) Development of dynamic stiffness method for free vibration of functionally graded Timoshenko beams Computers & Structures, 147:107–116 14 Long, N V., cs / Tạp chí Khoa học Công nghệ Xây dựng [13] Le, T H., Gan, B S., Trinh, T H., Nguyen, D K (2014) Finite element analysis of multi-span functionally graded beams under a moving harmonic load Mechanical Engineering Journal, 1(3):CM0013– CM0013 [14] Bellifa, H., Benrahou, K H., Hadji, L., Houari, M S A., Tounsi, A (2015) Bending and free vibration analysis of functionally graded plates using a simple shear deformation theory and the concept the neutral surface position Journal of the Brazilian Society of Mechanical Sciences and Engineering, 38(1):265– 275 [15] Zhang, D.-G (2013) Modeling and analysis of FGM rectangular plates based on physical neutral surface and high order shear deformation theory International Journal of Mechanical Sciences, 68:92–104 [16] Zhang, D.-G., Zhou, Y.-H (2008) A theoretical analysis of FGM thin plates based on physical neutral surface Computational Materials Science, 44(2):716–720 [17] Yaghoobi, H., Fereidoon, A (2010) Influence of neutral surface position on deflection of functionally graded beam under uniformly distributed load World Applied Sciences Journal, 10(3):337–341 [18] Larbi, L O., Kaci, A., Houari, M S A., Tounsi, A (2013) An Efficient Shear Deformation Beam Theory Based on Neutral Surface Position for Bending and Free Vibration of Functionally Graded Beams# Mechanics Based Design of Structures and Machines, 41(4):421–433 [19] Li, J., Wang, G., Guan, Y., Zhao, G., Lin, J., Naceur, H., Coutellier, D (2021) Meshless analysis of bidirectional functionally graded beam structures based on physical neutral surface Composite Structures, 259:113502 [20] Zhang, D.-G (2013) Nonlinear bending analysis of FGM beams based on physical neutral surface and high order shear deformation theory Composite Structures, 100:121–126 [21] Chi, S.-H., Chung, Y.-L (2006) Mechanical behavior of functionally graded material plates under transverse load—Part I: Analysis International Journal of Solids and Structures, 43(13):3657–3674 [22] Atmane, H A., Tounsi, A., Bernard, F (2015) Effect of thickness stretching and porosity on mechanical response of a functionally graded beams resting on elastic foundations International Journal of Mechanics and Materials in Design, 13(1):71–84 [23] Reddy, J N (2006) Theory and Analysis of Elastic Plates and Shells CRC Press [24] Reddy, J N (2017) Energy principles and variational methods in applied mechanics John Wiley & Sons [25] Dym, C L., Shames, I H (1973) Solid mechanics Springer 15 ... Chakraverty [10] phân tích tần số dao động riêng dầm FGM với liên kết khác theo lý thuyết dầm cổ điển dầm bậc Mashat cs [11] phân tích dao động tự dầm FGM với điều kiện biên khác sử dụng lý thuyết biến... sử dụng lý thuyết biến dạng cắt bậc cao lý thuyết dầm cổ điển Li [4] đề xuất cách tiếp cận hợp nhất, xây dựng lời giải giải tích phân tích uốn dao động riêng dầm FGM sử dụng hai mơ hình dầm Timoshenko... pháp độ cứng động lực phân tích dao động riêng dầm FGM theo lý thuyết dầm Timoshenko với liên kết khác Le cs [13] phân tích đáp ứng động dầm FGM nhiều nhịp tác dụng tải trọng di động phương pháp