Bài viết này phát triển Mô hình độ cứng động lực (DSM) để phân tích dao động tự do của dầm bậc nano bằng vật liệu có cơ tính biến thiên (FGM) trên cơ sở Lý thuyết đàn hồi không cục bộ (NET), gọi là mô hình DSMNET.
Tạp chí Khoa học Cơng nghệ Xây dựng, NUCE 2021 15 (2V): 49–64 PHÂN TÍCH DAO ĐỘNG TỰ DO CỦA DẦM BẬC NANO BẰNG VẬT LIỆU CĨ CƠ TÍNH BIẾN THIÊN SỬ DỤNG MƠ HÌNH ĐỘ CỨNG ĐỘNG LỰC KHƠNG CỤC BỘ Phạm Tuấn Đạta , Trần Văn Liênb,∗, Ngô Trọng Đứcc a Viện Kỹ thuật cơng trình đặc biệt, Học viện Kỹ thuật quân sự, 236 đường Hoàng Quốc Việt, quận Cầu Giấy, Hà Nội, Việt Nam b 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 c Fujita Corporation, 52 đường Lê Đại Hành, quận Hai Bà Trưng, Hà Nội, Việt Nam Nhận ngày 15/03/2021, Sửa xong 05/05/2021, Chấp nhận đăng 06/05/2021 Tóm tắt Bài báo phát triển Mơ hình độ cứng động lực (DSM) để phân tích dao động tự dầm bậc nano vật liệu có tính biến thiên (FGM) sở Lý thuyết đàn hồi khơng cục (NET), gọi mơ hình DSMNET NET có xét đến tham số tỷ lệ chiều dài nên giữ lại hiệu ứng tỷ lệ cho cấu trúc nano xét đến tương tác nguyên tử phân tử không liền kề Đặc trưng vật liệu dầm nano FGM thay đổi phi tuyến theo độ dày dầm Dầm bậc nano mô hình hóa theo lý thuyết dầm Timoshenko phương trình chuyển động rút từ nguyên lý Hamilton DSM sử dụng để thu nghiệm xác phương trình chuyển động có xét đến vị trí thực trục trung hòa với điều kiện biên khác So sánh kết tính tốn DSM-NET với kết công bố khẳng định độ tin cậy mơ hình Từ đó, tác giả tiến hành khảo sát số nhằm đánh giá ảnh hưởng tham số phân bố vật liệu, hình học, khơng cục điều kiện biên dao động tự dầm bậc nano FGM Nghiên cứu áp dụng cho nhiều kết cấu nano FGM khác dầm liên tục hay khung nano nhiều bậc phức tạp Từ khoá: mơ hình độ cứng động lực; lý thuyết đàn hồi khơng cục bộ; vật liệu có tính biến thiên; dầm bậc nano; tần số không thứ nguyên FREE VIBRATION ANALYSIS OF FGM STEPPED NANOBEAMS USING NONLOCAL DYNAMIC STIFFNESS MODEL Abstract This paper presents a nonlocal Dynamic Stiffness Model (DSM) for free vibration analysis of Functionally Graded (FG) stepped nanobeams using the Nonlocal Elastic Theory (NET), called DSM-NET model The NET nanobeam model considers the length scale parameter, which can capture the small scale effect of nano structures considering the interactions of non-adjacent atoms and molecules Material characteristics of FG nanobeams are considered nonlinearly varying throughout the thickness of the beam The nanobeam is modelled according to the Timoshenko beam theory and its equations of motion are derived using Hamilton’s principle The DSM is used to obtain an exact solution of the equation of motion taking into account the neutral axis position with different boundary conditions The DSM-NET is validated by comparing the obtained results with published results Numerical results are presented to show the significance of the material distribution profile, nonlocal effect, and boundary conditions on the free vibration of nanobeams It is shown that the study can be applied to other FG stepped nanobeams as well as more complex of framed nanostructures Keywords: DSM; NET; FGM; stepped nanobeam; nondimensional frequency https://doi.org/10.31814/stce.nuce2021-15(2V)-05 © 2021 Trường Đại học Xây dựng (NUCE) ∗ Tác giả đại diện Địa e-mail: lientv@nuce.edu.vn (Liên, T V.) 49 Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Giới thiệu Vật liệu có tính biến thiên (FGM) [1, 2] vật liệu tổng hợp hệ tạo thành từ hai nhiều vật liệu thành phần với thay đổi liên tục tỷ lệ thành phần theo nhiều hướng FGM sử dụng hệ thống điện micro/nano (MEMS/NEMS) để đạt độ nhạy cao hiệu suất mong muốn Các cấu trúc có kích thước nano dầm, vỏ sử dụng rộng rãi thiết bị NEMS, dầm bậc nano đặc biệt thu hút ngày nhiều ý ứng dụng tiềm khác chúng Lý thuyết đàn hồi cổ điển nghiên cứu đầy đủ xác ứng xử học cấu trúc nano ảnh hưởng hiệu ứng kích thước Do đó, lý thuyết đàn hồi không cục (NET) lần đề xuất Eringen [3] với giả thiết tensor ứng suất điểm không hàm biến dạng mà cịn xét đến biến dạng điểm xung quanh Hiện nay, NET sử dụng rộng rãi để thiết lập phương trình chuyển động cấu trúc nano sử dụng vật liệu đồng [4–6] vật liệu FGM [7] Reddy [8] thiết lập phương trình dao động ổn định dầm nano đồng theo NET cho lý thuyết dầm Euler-Bernoulli, Timoshenko, Reddy Levinson Nhiều tác giả khác phát triển phương pháp giải tích [9–12], phương pháp Rayleight-Ritz [13], phương pháp phần tử hữu hạn (FEM) [14–19], phương pháp biến đổi vi phân [20], phương pháp cầu phương vi phân [21], để nghiên cứu ứng xử uốn, ổn định dao động tự nano từ vật liệu đồng Simsek Yurtcu [22], Rahmani Pedram [23] đồng thời nghiên cứu ứng xử uốn ổn định dầm Timoshenko FGM phương pháp giải tích Ngồi ra, Mechab cs [24] nghiên cứu dao động tự do, Uymaz [25] nghiên cứu dao động cưỡng dầm nano, hai sử dụng lý thuyết biến dạng cắt bậc cao Ebrahimi Salari [26] ứng dụng phương pháp bán giải tích để phân tích dao động dầm nano Euler - Bernoulli FGM có xét đến vị trí thực trục trung hịa Narendar Gopalakrishnan [27] phát triển phương pháp PTHH phổ để khảo sát dao động dầm nano không cục Các nghiệm giải tích tìm có dạng chuỗi lượng giác kép Navier, đó, chúng bị giới hạn cho dầm với điều kiện gối tựa đơn giản Đối với điều kiện biên khác, tác giả áp dụng FEM để phân tích ứng xử uốn dao động tự dầm nano FGM theo lý thuyết dầm Euler - Bernoulli [28, 29], lý thuyết Timoshenko [30–32] Gần đây, tác giả [33] tìm tần số riêng dạng dao động dầm nano với điều kiện biên khác cách sử dụng khái niệm khơng gian trạng thái Vì FEM xây dựng dựa hàm dạng đa thức không phụ thuộc vào tần số nên khơng dùng để tìm tất tần số dạng dao động, đặc biệt bậc cao Để khắc phục điều này, phương pháp độ cứng động lực (DSM) phát triển để thay FEM cách sử dụng hàm dạng phụ thuộc tần số [34–36] Mặc dù không dễ dàng tìm nghiệm giải tích xác cho toán dao động dầm nano tìm chúng cho phép nghiên cứu dao động xác dầm dải tần số tùy ý Adhikari cs [4, 37] thu ma trận độ cứng động lực nano đồng chịu dao động dọc Hàm đáp ứng tần số thu dùng DSM đề xuất cho thấy mật độ phổ cao gần tần số cực đại Gần đây, Taima cs [38] sử dụng DSM để phân tích dao động cho dầm bậc nano đồng theo lý thuyết dầm Euler-Bernoulli Theo tìm hiểu tác giả, chưa có nghiên cứu phát triển DSM cho toán dao động dầm nano FGM theo lý thuyết đàn hồi không cục Trong báo này, tác giả phát triển DSM để phân tích dao động tự dầm bậc nano FGM sở lý thuyết dầm Timoshenko lý thuyết đàn hồi khơng cục bộ, gọi mơ hình DSM-NET So sánh kết thu DSMNET với kết công bố cho thấy độ tin cậy mơ hình đề xuất Từ đó, tác giả nghiên cứu ảnh hưởng tham số không cục bộ, phân bố vật liệu tham số hình học đến ba tần 50 Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng số dao động không thứ nguyên đầu DSM-NET tiên dầmCỦA nanoDẦM với điều kiện biên khác MƠ HÌNH TIMOSHENKO BẰNG VẬT LIỆU FGM Mơ hình DSM-NET liệu FGM Đớidầm với Timoshenko mợt dầm nanobằng FGMvật (Hình 1), đặc tính vật liệu thay đổi theo chiều d sau [1]: n z Trục trung hòa Et Gt t t z 1 h x Eb Gb b b P( z ) Pb Pt Pb h 2 h h z 2 b với P lần lượt ký hiệu môđun đàn hồi E, môđun cắt mật độ khối lượng ; ch t b liên quan đến giá trị tương ứng vật liệu lớp dày lớp zdưới Đối với dầm nano FGM (Hình 1), đặc tính vật liệu thay đổi theo chiều nhưcùng; n ch phần thể tích, z tọa đợ từ mặt phẳng dầm Các chuyển vị dầm Timoshenko là: sau [1]: n ) t − u0 P( xb,)t ) z (+z 1 h0 ) ( x, t ) ; w( x, z, t ) w0 ( x, t ) P(z) =u(Pxb, z+, t(P h (1) đó u0 (x, t), w0 (x, t) hlần lượt h chuyển vị dọc, chuyển vị ngang một điểm trục t − ≤ z ≤ hòa; h0 khoảng cách từ 2trục trung2 hòa đến trục x; góc quay mặt cắt quanh trục y Từ ta nhận biến dạng: với P ký hiệu môđun đàn hồi E, môđun cắt G mật độ khối lượng ρ; số vật u0 /liệu x lớp ( z trên h0 )và /lớp x ; xz cùng; w0 / nx chỉ số phần thể t b liên quan đến giá trị tương ứng xxcủa tích, z tọa độ từ mặt phẳngCác giữaquan củahệ dầm chuyển dầm Timoshenko vật Các lý khơng cụcvịbợ cho dầm nano có thểlà: viết dạng [3]: Hình Dầm FGM Hình Dầm FGM u(x, z, t) = u0 (x, t) − (z − h0 )θ(x, t); 2 xxw(x, z, t) = w0 (x,t)2 xz xx E xx ; xz G xz (2) xvị ngang điểm x trục trung hịa; u0 (x, t), w0 (x, t) chuyển vị dọc, chuyển 2trục x; θ góc quay mặt cắt quanh trục y Từ đó, ta nhận h0 khoảng cách từ trục trung eđến đó hịa a tham sớ không cục bộ; e0 số cho vật liệu; a chiều dài biến dạng: trưng bên dầm 2 Dựa trên−nguyên lý Hamilton: ε xx = ∂u (z − h0 )∂θ/∂x; γ xz = ∂w0 /∂x − θ = −ϕ /∂x (3) T Các quan hệ vật lý không cục cho dầm nano có thểđược viết T U [3]: dtởdạng ∂2 σ xx ∂2 σ xz σ xxT−vൠU lần = Eε ; σ − µ = Gγ (4) xx xz xz dầm: đó lượt là động và ∂x2 ∂x2 L u0 u0 w0 w0 u0 dài đặc alà chiều uvật liệu; µ = e20 a2 tham Tsốkhông bộ; e0là số cho I12 I 22 trưngdx I11 cục t t t t t t t t t t bên dầm 0 w0 u0 T U N M Q Q dx x x x 0 Dựa nguyên lý Hamilton: L (δT − δU)dt = N, M, Q lần lượt lực dọc, mômen uốn lực cắt dầm: (5) N của xx dA ; M z h0 xx dA ; Q xz dA T U động dầm: A L δT = A A ∂u ∂δu ∂w0 ∂δw0 ∂u0 ∂δθ ∂θ ∂δu0 ∂θ ∂δθ 0 I11 + − I12 + dx 2+ I22 ∂t ∂t ∂t A11∂t, A12 , A22 ∂t ∂t ∂t E ( z )∂t1, z h∂t , z h dA ; A33 ∂t G ( z )dA 0 A L A (6) ∂δu0 ∂δθ , I ∂δw ( z ) 1, z h , z h dA 12 Q22 − − M I11 , I+ Qδθ dx ∂x ∂x ∂x A hệ số hiệu chỉnh cắt, =5/6 đới với mặt cắt ngang hình chữ nhật; vị trí trục trung hị 51 I12 có giá trị là: δU = N Đạt, P T., cs / Tạp chí Khoa học Công nghệ Xây dựng N, M, Q lực dọc, mômen uốn lực cắt dầm: N= σ xx dA; M= A (z − h0 ) σ xx dA; Q= A σ xz dA (7) A (A11 , A12 , A22 ) = E(z) 1, z − h0 , (z − h0 )2 dA; A33 = η G(z)dA A A (I11 , I12 , I22 ) = (8) ρ(z) 1, z − h0 , (z − h0 ) dA A η hệ số hiệu chỉnh cắt, η = 5/6 mặt cắt ngang hình chữ nhật; vị trí trục trung hịa h0 I12 có giá trị là: h0 = n(RE − 1)h ; 2(n + 2)(n + RE ) I12 = nbh2 ρb Rρ − RE ; 2(2 + n) n + RE Rρ = ρt ; ρb RE = Et Eb (9) Từ (5), ta nhận phương trình dao động tự dầm nano FGM có dạng: ∂2 u0 ∂2 θ ∂2 u0 ∂2 θ ∂4 u0 ∂4 θ =0 − A − I + I + µ I − I 12 11 12 11 12 ∂x2 ∂x2 ∂t2 ∂t2 ∂x2 ∂t2 ∂x2 ∂t2 ∂2 θ ∂2 θ ∂2 u0 ∂w0 ∂2 u0 ∂4 u0 ∂4 θ A22 − A12 + A33 − θ − I22 + I12 − µI12 2 + µI22 2 = ∂x ∂x ∂x ∂t ∂t ∂x ∂t ∂x ∂t ∂2 w0 ∂2 w0 ∂θ ∂4 w0 − I11 + µI11 2 = A33 − ∂x ∂x2 ∂t ∂x ∂t A11 (10) điều kiện biên tương ứng ∂3 θ ∂θ ∂3 u0 ∂u0 =0 − I − A12 + µ I11 12 ∂x ∂x ∂x∂t2 ∂x∂t2 ∂w0 ∂3 w0 w0 = or A33 − θ + µI11 =0 ∂x ∂x∂t2 ∂u0 ∂θ ∂2 w0 ∂3 u0 ∂3 θ θ0 = or A12 − A22 + µ I11 + I12 =0 − I 22 ∂x ∂x ∂t ∂x∂t2 ∂x∂t2 u0 = or A11 Đặt: (11) ∞ {U, Θ, W} = {u0 (x, t), θ(x, t), w0 (x, t)} e−iωt dt (12) −∞ Phương trình dao động (10) miền tần số có dạng: d2 U d Θ − A − µI ω + I11 ω2 U − I12 ω2 Θ = 12 12 dx2 dx2 dW d2 U d Θ + A − µI ω + A33 − I12 ω2 U + I22 ω2 − A33 Θ = 22 22 2 dx dx dx d2 W dΘ A33 − µI11 ω2 − A33 + I11 ω2 W = dx dx A11 − µI11 ω2 − A12 − µI12 ω2 52 (13) Đạt, P T., cs / Tạp chí Khoa học Công nghệ Xây dựng Đưa vào ma trận vectơ: A11 − µI11 ω2 − A12 − µI12 ω2 ˜ = − A12 − µI12 ω2 A A22 − µI22 ω 0 A33 − µI11 ω2 −I12 ω2 0 I11 ω2 ˜ = −I12 ω2 I22 ω2 − A33 ˜ = 0 A33 ; B C ; −A33 0 I11 ω2 (14) U Θ {z} = W Phương trình (13) viết dạng ma trận: ˜ z A ˜ {z} = {0} ˜ z + C + B (15) Chọn nghiệm phương trình (15) dạng {z0 } = {d} eλx dẫn đến ˜ {d} = {0} ˜ +λ B ˜ + C λ2 A (16) Phương trình (16) có nghiệm ˜ =0 ˜ +λ B ˜ + C det λ2 A (17) Từ ta nhận phương trình bậc ba η = λ2 : η3 + aη2 + bη + c = Ký hiệu η1 , η2 , η3 nghiệm phương trình bậc ba λ1,4 = ±k1 ; λ2,5 = ±k2 ; λ3,6 = ±k3 ; kj = √ η j; j = 1, 2, Nghiệm tổng qt phương trình (15) có dạng {z0 (x, ω)} = (18) d j eλ j x Từ phương trình j=1 phương trình thứ ba (16), ta có: α1C1 α2C2 α3C3 α4C4 α5C5 α6C6 C2 C3 C4 C5 C6 {z0 (x, ω)} = C1 β1C1 β2C2 β3C3 β4C4 β5C5 β6C6 k1 x k2 x k3 x −k1 x −k2 x −k3 x e e · e e e e T (19) {C} = (C1 , , C )T số độc lập α1 = A12 − µI12 ω2 k12 + I12 ω2 A11 − µI11 ω2 k12 + I11 ω2 = α4 ; β1 = A33 k1 = −β4 A33 − µI11 ω2 k12 + I11 ω2 (20) Tương tự, α2 = α5 ; β2 = −β5 ; α3 = α6 ; β3 = −β6 , phương trình (19) trở thành: {U, Θ, W}T = [G(x, ω)] {C} (21) α1 ek1 x α2 ek2 x α3 ek3 x α1 e−k1 x α2 e−k2 x α3 e−k3 x k1 x k2 x k3 x −k1 x −k2 x [G(x, ω)] = e e e e e e−k3 x β1 ek1 x β2 ek2 x β3 ek3 x −β1 e−k1 x −β2 e−k2 x −β3 e−k3 x 53 (22) 1 A A 12 I12 k12 I12 2 2 11 I11 k1 I11 ; 1 A33k1 4 2 A I k I 33 11 11 (20 Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng 2một phần nano 5 ;FGM 3 haichiều Hình 6 2., phương tử2dầm Chuyển vị trình lực đầu phần Tương tự,Xét (19) trởnút thành: 5; ; 3như đó tử có dạng: , ,W G ( x, ) C ˆ e = {U1 ,U U Θ1 , W1 , U2 , Θ2 , W2 }T ; {Pe } = {N1 , M1 , Q1 , N2 , M2 , Q2 }T T (23) (21 k3 x k1x k x 1ek1x 2Ue1k2=x z1 (0, (0, ω); W z3 (0,2 ω) ω); 3e k3 x =2 e 3e Θ1 =z12e kx Uk2 x= z1 (L, ω); Θ = z2 (L, ω); W2 = z3k(L, ω) k3 x k3 x k1x 2x 2 (22 G(xN,1 =)− A11e− µI11ω ∂exz1 − A12 −eµI12ω ∂xze2 x=0; Q1 = e− A33 − µI11ωe ∂xz3 − A33z2 x=0 kM1x1 = − A12k2−x µI12 ω2 k∂3xxz1 − A22 −µI k x k1x22 ω2 ∂ x z2 −µI k2 11 x ω2 z3 x=0e e e e e e 3 z 2 N2 = A11 − µI11 ω ∂ x z1 − A12 − µI12 ω ∂ x z2 ; Q2 = A33 − µI11 ω ∂ x z3 − A 33 x=L x=L 2 M2nano = A12FGM − µI12 ωhai∂ xchiều z1 − A22như − µI22 ω ∂ x z2 µI11 ω z3 vị lực tại đầu nú −Chuyển Xét một phần tử dầm Hình x=L (24) hần tử có dạng: đó y Q2 Q1 N1 L x i M1 N2 j W1 W2 M2 U1 U2 1 2 Hình Nút lựcvànút Timoshenko Hìnhvà Nút lực nút phần phần tử tử dầmdầm Timoshenko U {U , ,W ,U , ,W } Thay ˆ (21) vào (24), ta e 1 ˆe = U T 2 [G(0, ω)] [G(L, ω)] · {C} ; ; Pe {N1 , M1 , Q1 , N , M , Q2 }T {Pe } = −BF (G) x=0 BF (G) x=L với [BF ] tốn tử ma trân A11 − µI11 ω2 ∂ x − A12 − µI12 ω2 ∂ x [BF ] = A12 − µI12 ω2 ∂ x − A22 − µI22 ω2 ∂ x −A33 (25) · {C} −µI11 ω2 A33 − µI11 ω2 ∂ x (26) Khử vectơ số {C} (25) dẫn đến {Pe (ω)} = −BF (G) x=0 BF (G) x=L [G(0, ω)] · [G(L, ω)] −1 ˆe = K ˆ e (ω) · U ˆe · U ˆ e ma trận độ cúng động lực phần tử dầm nano FGM K 54 (27) (23 Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng Đối với kết cấu dầm, khung gồm số phần tử dầm nano FGM trên, cách cân ˆ nội lực nút kết cấu, ta thu ma trận độ cứng động lực K(ω) vectơ lực nút Pˆ ˆ vectơ chuyển vị nút kết cấu, phương trình chuyển động kết cấu toàn kết cấu Gọi U theo DSM ˆ ˆ = Pˆ K(ω) · U (28) Do đó, tần số dao động tự {ω} = {ω1 ω2 ωn } tìm thấy từ phương trình ˆ det[K(ω)] =0 (29) Kết tính tốn thảo luận 3.1 So sánh kết tính tốn với kết cơng bố Trong mục này, kết tính tốn số theo mơ hình DSM-NET đề xuất so sánh với kết công bố để kiểm tra độ tin cậy mô hình So sánh thực tần số không nguyên λ1 = ω1 L2 ρt A/Et I dầm nano đồng không cục đơn giản theo bốn lý thuyết dầm không cục bộ: Euler-Bernoulli (EBT), Timoshenko (TBT), Reddy (RBT) Levinson (LBT) Các kết từ mơ hình DSM-NET đề xuất Bảng phù hợp với kết giải tích Reddy [8] với thơng số không cục khác Bảng So sánh tần số không thứ nguyên λ1 = ω1 L2 ρt A/Et I dầm nano đơn giản với tham số không cục khác (tỷ số L/h = 10) µ (10−18 ) EBT [8] TBT [8] RBT [8] LBT [8] Kết 0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 9,8696 9,6347 9,4159 9,2113 9,0195 8,8392 8,6693 8,5088 8,3569 8,2129 8,0761 9,7454 9,5135 9,2973 9,0953 8,9059 8,7279 8,5601 8,4017 8,2517 8,1095 7,9744 9,7454 9,5135 9,2974 9,0954 8,9060 8,7279 8,5602 8,4017 8,2517 8,1095 7,9744 9,7657 9,5333 9,3168 9,1144 8,9246 8,7462 8,5780 8,4193 8,2690 8,1265 7,9911 9,7451 9,5132 9,2971 9,0951 8,9057 8,7277 8,5599 8,4015 8,2515 8,1093 7,9742 So sánh thứ hai thực cho dầm nano FGM đơn giản với thông số hình học vật liệu là: L = 10 nm, b = h = L/10 = nm ; Et = 70 Gpa, Gt = 26 Gpa, ρt = 2700 kg/m3 ; Eb = 393 Gpa, Gb = 157 Gpa, ρb = 3960 kg/m3 [30] Kết tính tốn cho tần số không thứ nguyên λi = ωi L2 ρt A/Et I, i = 1, 2, trục trung hòa (NA) trục mặt phẳng (MA) theo DSM-NET so sánh với kết Eltaher cs [30] sử dụng FEM với 100 phần tử thể Bảng Có thể thấy kết tốt, đặc biệt tính đến vị trí thực trục trung hịa sai số nhỏ 0,5% Các so sánh kết tính tốn với kết cơng bố cho thấy độ tin cậy mơ hình DSM-NET đề xuất 55 Đạt, P T., cs / Tạp chí Khoa học Công nghệ Xây dựng Bảng So sánh tần số không thứ nguyên dầm nano với kết Eltaher cs µ (10−19 ) Frqs (λi ) MA (Kết quả) NA (Kết quả) MA [30] NA [30] MA (Kết quả) n=0 9,7075 37,0962 78,1547 9,7075 37,0962 78,1547 9,7032 37,0382 77,9135 17,0278 64,1874 136,4469 16,7305 64,2634 136,2956 17,1057 65,6963 139,3094 9,7032 37,0382 77,9135 14,5998 51,8360 110,6544 8,8713 27,7303 46,9034 8,8713 27,7303 46,9034 2,0 8,8674 27,6870 46,7622 16,7366 64,3491 136,6522 17,6526 67,2760 142,5265 15,5611 48,0080 81,9216 15,2894 48,0385 81,7958 15,6323 49,1097 83,6044 8,8674 27,6870 46,7622 13,3416 38,8463 66,5345 8,2196 23,0989 36,6272 8,2196 23,0989 36,6272 4,0 8,2160 23,0628 36,5169 15,2950 48,1026 82,0098 16,1321 50,3002 85,5480 14,4179 39,9980 63,9788 14,1662 40,0154 63,8749 14,4840 40,9076 65,2873 13,6100 52,1595 110,3186 17,5472 67,3039 142,4755 17,6845 67,8750 143,8072 17,5542 67,4009 142,8763 12,4351 38,9595 66,1005 13,5863 42,3969 71,5567 12,4377 38,9905 66,2061 16,0358 50,3114 85,5045 16,1613 50,7383 86,3037 16,0421 50,3839 85,7451 n = 0,5 8,2160 23,0628 36,5169 12,3611 32,3883 51,9778 n=5 14,8669 52,7163 119,2343 n = 10 n=0 13,6072 52,1179 110,1423 n = 0,5 n=5 NA [30] n = 10 n=0 MA [30] n = 0,5 n=5 NA (Kết quả) 11,5216 32,4527 51,6184 12,5883 35,3159 55,8791 11,5241 32,4785 51,7008 n = 10 14,1714 40,0687 64,0420 14,9470 41,9022 66,8071 h1 h1 b 14,8578 41,9086 66,7711 14,9741 42,2642 67,3952 14,8637 41,9690 66,6590 3.2 Dao động tự dầm bậc nano FGM b h h hh L1 b1 h b1 h L1 L1 L L1 L L (a) Chiều rộng (b) Chiều caoL (a) Chiều Chiều rộng bậcChiều cao (a) rợng Hình Dầm nano FGM thay đổi dạng(b) (b) Chiều cao Hình nanovớiFGM thay đổi học dạng bậc Trong mục này, dầm3 bậcDầm nano FGM thơng sớ hình vật liệu sau: L=10nm, Hình Dầm nano 3FGM đổi dạng bậc 3, vt=vb=0.3, n=0.5 [30] b=h=1nm, Et=70Gpa, ; Eb= thay 393Gpa, t=2700kg/m b=3960kg/m đượcnày, nghiên cứu (Hình 4) CácFGM kết quảvới đượccác tínhthơng tốn chosớbahình tần sớhọc khơng ngun Trong mục dầm bậc nano vàthứ vật liệu đầu nhưtiênsau: L=10nm, xét đến trục trung hòa với điều kiện3biên: đơn giản (S-S), L A E I ; i 1, 2,3 có i i bậc nano tt =2700kg/m t FGM với ;các b=h=1nm, Et=70Gpa, Ebthông = 393Gpa, =0.3, [30] b=3960kg/m t=vbsau: Trong mục này, dầm số hìnhhọc vật liệu, vnhư L =n=0.5 10 nm, cơng (Hình xơn (CF),4) hai kết đầu (CC) – gối tựa (C-Scho ) vớiba trường hợpkhông nghiên cứu: nghiên quả tính tốn tần sớ thứ nguyên đầu tiên b = hđược = nm, Et =cứu 70 Gpa, ρt ngàm = Các 2700 kg/m ; Engàm = 393 Gpa, ρ = 3960 kg/m , v = v = 0,3, n = b b t b A1) Dầm nano có chiều cao bậc (h1/h = 0,8) tại vị trí khác nhau: L1/L = 0, 0,1, , 1,0 đó 0,5 [30] nghiên 4) Cácxét kếtđến quảbợtrục tốn tần sốtồn khơng thứ nguyên đầu(S-S), i được A1/LEcứu i(Hình 1,ứng 2,3 có với điều kiện biên: đơn giản =t I0 ;tương h1=0,8h tồn dầm,trung Ltính tươngcho ứngba h1=h bợ dầm 1/L=1hịa i L t L tiên λcông = ω L ρ A/E I; i = 1, 2, có xét đến trục trung hịa với điều kiện biên: đơn giản (S-S), B1) Dầm nano có chiều rợng bậc (b /b = 0,8) tại vị trí khác nhau: L /L = 0, 0,1, , 1,0 i ixôn (CF), t ngàm t hai đầu (CC) ngàm – gốibtựa (C-S ) ứng với 4=1trường hợpdầm nghiên cứu: b1/bđầu = tương =0,8b tồntựa bợ dầm, = tương b tồn bộcứu: công xôn (CF), ngàmđóhai (CC)ứng vàb1ngàm – gối (C-S1/b) với trườngb1hợp nghiên A1) Dầm nanoA2) cóDầm chiều cao bậc (h /h = 0,8) tại vị trí khác nhau: L /L = 0, 0,1, nano có chiều cao 1bậc tại vị trí L1/L = 0,5 với tỷ lệ khác nhau: h1/h 1= 1,0,0,8, , 0,2; , 1,0 đó Dầmhnano có chiều bậcbợ tại vị trí L1/L 0,5 vớitương tỷ lệ khác nhau: 1,0,0,8, , bộ 0,2 dầm 56 1/b = L1/L = tươngB2)ứng trênrợng tồn dầm, L1=/L=1 ứng h1b=h tồn 1=0,8h Hình thể hiện sự thay đổi ba tần số không thứ nguyên đầu tiên dầm bậc nano FGM với B1) Dầm nanocáccóvị chiều rợng bậc (b1/bbợ=và0,8) tại vị trí khác nhau: L1/L = 0, 0,1, , 1,0 trí bậc, tham sớ khơng cục điều kiện biên khác Hình thể hiện sự thay đổi đó b1/b = tương b1tiên =0,8b dầm, b1tại/bvị=trí1Ltương b1tỷ=lệbbước, trêntham tồnsớbợ dầm ba tầnứng sớ đầu dầmtồn nanobợ FGM có bậc với 1/L = 0,5ứng không cục bộ và điều kiện biên khác Từ hình vẽ ta có nhận xét sau: A2) Dầm nano có chiều cao bậc tại vị trí L1/L = 0,5 với tỷ lệ khác nhau: h1/h = 1,0,0,8, , 0,2; a) Tồn tại vị trí bậc dầm FGM có ảnh hưởng lớn đến một tần số không thứ nguyên nhất sớsớ phần thểthể tích và tham sớ sớ khơng cụccục bợ bợ tăng (Hình 6b, 6b, 7b).7b) phần tích khơng tăng (Hình Sự thay sớ khơng ngun sớ chỉphần sớ phần thể tích gâylàra rệt so so b)b)Sự thay đổiđổi tầntần sốtham cơ bảnbản khơng thứthứ ngun thể tích gây rõlàrệtrõhơn b) b)SựSự thay đổi tần số bản không thứ ngun sớ phần thể thể tíchtích gâygây làrarõlàrệt so so thay đổi tần số bản không thứ nguyên số phần rõ rệt tham không vớivới tham số số không cụccục bộ bộ gâygây ra vớivới tham số không cục bộ gây tham số không cục bộ gây c)Khi Khi lệ tăng tần sốbản bản không nguyên với điều kiện biên S-S theo lệbậc bậcbậc tăng lên,lên, cáccác tầnsố sốcơ bản không thứ thứ nguyên ứngứng với điều kiện biên S-S c)c) tỷtỷtỷ lệtỷ tăng lên, tần không thứ nguyên ứng với với điềuđiều kiệnkiện biênbiên S-S theo theotheo c)Khi Khi lệ bậc tăng lên, tần số bản không thứ nguyên ứng S-S một xu hướng tăng tần số điều kiện biên C-F giảm một xu hướng tăng tần số điều kiện biên C-F giảm một xuxu hướng tăng cáccác tầntần số điềuđiều kiệnkiện biên C-FC-F giảm một hướng tăng số giảm Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây biên dựng (a)(a) (b) (b) (a) (b) (a)(a) (b) (b) (c)(c) (d) (d) (c) (d) (c)(c) (d) (d) (c)(c) (d) (d) (c) (c) (d) (d) 9 9 (e) (f) (e) (e) (f) (f) 57 Đạt, (e) (e) P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng (f) (f) (e) (f) (e) (f) (g)(g) (g) (g)(g) (h) (h) (h) (h)(h) (i)(i)(i) (i) (i) (j) (j) (j) (j) (j) 1010 10 10 (k) (k) (k) (l) (l)(l) 4.tần Basốsố tần số không thứ nguyên đầu dầmFGM FGM bậc nano vớitrícác vịtham trí bậc, bậc, tham HìnhHình khơng thứthứ ngun đầu đầu tiên củatiên dầm FGM bậc nano với cácvới vị bậc, số phi địasớ Hình 4.Ba Ba tần khơng ngun tiên dầm bậc nano vị trí tham sớ địa phương điều kiện biên khácnhau phương điều kiện biên khác phiphi địa phương vàvà điều kiện biên khác 58 Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng A1) Dầm nano có chiều cao bậc (h1 /h = 0,8) vị trí khác nhau: L1 /L = 0, 0,1, , 1,0 L1 /L = tương ứng h1 = 0,8h toàn dầm, L1 /L = tương ứng h1 = h tồn dầm B1) Dầm nano có chiều rộng bậc (b1 /b = 0,8) vị trí khác nhau: L1 /L = 0, 0,1, , 1,0 b1 /b = tương ứng b1 = 0,8b toàn dầm, b1 /b = tương ứng b1 = b toàn dầm A2) Dầm nano có chiều cao bậc vị trí L1 /L = 0,5 với tỷ lệ khác nhau: h1 /h = 1,0, 0,8, , 0,2; B2) Dầm nano có chiều rộng bậc vị trí L1 /L = 0,5 với tỷ lệ khác nhau: b1 /b = 1,0, 0,8, , 0,2 Hình thể thay đổi ba tần số không thứ nguyên dầm bậc nano FGM với vị trí bậc, tham số khơng cục điều kiện biên khác Hình thể thay đổi ba tần số dầm nano FGM có bậc vị trí L1 /L = 0,5 với tỷ lệ bước, tham số không cục điều kiện biên khác Từ hình vẽ ta có nhận xét sau: a) Tồn vị trí bậc dầm FGM có ảnh hưởng lớn đến tần số khơng thứ ngun định Vị trí gọi vị trí tới hạn tần số cho trước Các vị trí điểm tới hạn tính từ đầu bên trái dầm nano 0,5L tần số ứng với điều kiện biên S-S (Hình 4(a)) C-F (Hình 4(d)); 0,2L với điều kiện biên C-C (Hình 4(g)) C-S (Hình 4(j)); 0,2L cho tần số thứ hai với điều kiện biên C-F (Hình 4(e)), 0,1L với điều kiện biên C-C (Hình 4(h)) C-S (Hình 4(k)) b) Tần số dầm bậc nano FGM nhạy cảm với vị trí bậc, tỷ lệ bậc điều kiện biên Độ lệch lớn tần số 3,5% (Hình 4(a)) 78,4% (Hình 5(a)) cho điều kiện biên S-S, 11,9% (Hình 4(g)) 32,89% (Hình 5(g)) cho điều kiện biên C-C, 12,99% (Hình 4(j)) 5,1% (Hình 5(j)) cho điều kiện biên C-S, lên đến 33,49% (Hình 4(d)) 243,1% (Hình 5(d)) cho điều kiện biên CF c) Các tần số không thứ nguyên ứng với điều kiện biên C-C cao sau C-S, S-S C-F hầu hết trường hợp d) Sự thay đổi tần số không thứ nguyên độ cao bậc gây vị trí khác rõ rệt so với độ rộng bậc gây e) Khi tham số khơng cục µ tăng, độ lệch lớn tần số không thứ nguyên ứng với điều kiện biên S-S, C-C C-S giảm độ lệch lớn tần số ứng với dầm cơng xơn tăng lên đến 35,2% (Hình 4(d)) 271,22% (Hình 5(d)) Khoảng cách hai tần số liên tiếp giảm xuống tham số không cục tăng f) Khi tỷ lệ bước tăng lên, tần số không thứ nguyên ứng với điều kiện biên S-S, C-C C-S tăng tần số không thứ nguyên ứng với điều kiện biên C-F giảm Đây gọi “nghịch lý không cục bộ”, số tác giả đề cập đến cho trường hợp vật liệu đồng [14, 38–40] Hình cho thấy thay đổi tần số không thứ nguyên dầm bậc nano FGM vị trí bậc dọc dầm với tham số không cục số phần thể tích khác với điều kiện biên S-S (Hình 6(a)) C-F (Hình 6(b)) Hình cho thấy thay đổi tần số không thứ nguyên dầm bậc nano FGM vị trí L1 /L = 0,5 với tỷ lệ bậc, tham số không cục số phần thể tích với điều kiện biên S-S (Hình 7(a)) C-F (Hình 7(b)) Từ đồ thị ta có nhận xét sau: a) Các tần số không thứ nguyên dầm nano bậc ứng với điều kiện biên S-S (và C-C, C-S) tăng lên số phần thể tích tăng tham số khơng cục giảm (Hình 6(a), 7(a)) Ngược lại, tần số không thứ nguyên dầm nano bậc ứng với điều kiện biên C-F tăng số phần thể tích tham số khơng cục tăng (Hình 6(b), 7(b)) b) Sự thay đổi tần số khơng thứ ngun số phần thể tích gây rõ rệt so với tham số không cục gây c) Khi tỷ lệ bậc tăng lên, tần số không thứ nguyên ứng với điều kiện biên S-S theo xu hướng tăng tần số điều kiện biên C-F giảm 59 (k) (k) (k) (k) (l)(l) (l) Hình 4.Ba Ba tần không thứ nguyên đầu tiên FGM bậc nano với vịvị trítrí bậc, tham sớsớ Hình 4.4 tần sốsố không thứ nguyên đầu dầm FGM bậc nano với trí bậc, tham sớ Hình tần sớsớ không thứ nguyên đầu tiên dầm FGM bậc bậc, tham Hình 4.Ba Ba tần khơng thứ ngun đầutiên tiêncủa củadầm dầm FGM bậc nano với vị trí phi địa phương và điều kiện biên khác phi địa phương phi địa phương và điều kiện biên khác phi địa phươngvà vàđiều điềukiện kiệnbiên biênkhác khácnhau Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng (a)(a) (a) (a) (b) (b) (b) (b) (b) (a) (c) (c) (c) (c) (d) (d) (d) (d) (d) (c) 11 11 11 11 (e) (e)(e) (f)(f)(f) 60 (f) (e) (f) dựng Đạt,(e) P T., cs / Tạp chí Khoa học Cơng nghệ Xây (e) (e) (f)(f) (g) (g) (h) (h) (h) (h)(h) (g)(g) (g) (i) (j) (i) (i) (i) (i) (j)(j) (j)(j) 12 1212 12 (k)(k) (k) (l)(l) (l) Hình 5.tần Basốtần sớ khơng thứ ngun đầu tiên dầm FGM bậc L1/L 0,5với tỷlệlệ Hình tần sốkhông không thứ nguyên tiên dầm FGM tại vịvị L1/L = =0,5 Hình 5.5.BaBa thứ nguyên đầuđầu tiên dầm FGM bậcbậc vịtại trí Ltrí1trí /L = 0,5 với tỷvới lệtỷbậc, bậc, tham số phi địa phương và điều kiện biên khác bậc,tham thamsốsốphi phiđịađịa phương điều kiện biên khác phương vàvà điều kiện biên khác 61 (k) (l) Hình Ba tần sớ khơng thứ nguyên đầu tiên dầm FGM bậc tại vị trí L1/L = 0,5 với tỷ lệ Đạt,kiện P biên T., khác cs / Tạp chí Khoa bậc, tham số phi địa phương và điều (a) (a) S-S S-S học Cơng nghệ Xây dựng (b)(b)C-F C-F Hình Tần số không thứ nguyên đầu tiên dầm FGM bậc với vị trí bậc, tham sớ phi địa phương sớ phần thể tích khác với điều kiện biên Hình Tần số khơng thứ nguyên dầm FGM bậc với vị trí bậc, tham số phi địa phương (b) C-F sốbậc phần khác Hình Tần sớ khơng thứ nguyên đầu tiên dầm FGM với cácthể vị trítích bậc, tham sớ phi địa với điều kiện biên phương sớ phần thể tích khác với điều kiện biên 13 (a) S-S 14 (a) (a)S-S S-S (b) (b)C-F C-F Hình Tần sớ khơng thứ nguyên đầu tiên dầm FGM bậc với tỷ lệ bậc, tham sớ phi địa Hình Tần số không 14thứ nguyên dầm FGM bậc với tỷ lệ bậc, tham số phi địa phương phương và sớ phần thể tích khác với điều kiện biên a) S-S, b) C-F số phần thể tích khác nhau4 KẾT với điều kiện biên a) S-S, b) C-F LUẬN Kết luận Trong bài báo này, tác giả phát triển mơ hình DSM để phân tích dao đợng tự dầm bậc nano FGM dựa NET, gọi là mơ hình DSM-NET Dầm bậc nano FGM mơ hình hóa theo lý thuyết dầm Timoshenko với phương trình chuyển đợng rút từ ngun lý Hamilton Dùng DSM-NET, tác giả tìm tần sớ và dạng dao động riêng kết cấu nano FGM dạng dầm, khung có xét đến vị trí thực trục trung hòa cho điều kiện biên khác So sánh kết quả tính tốn DSM-NET với kết quả công bố khẳng định độ tin cậy mơ hình Đây là kết quả mới, nhóm tác giả phát triển thời gian gần Từ đó, tác giả tiến hành khảo sát số nhằm đánh giá ảnh hưởng tham sớ phân bớ vật liệu, hình học, khơng cục bợ và điều kiện biên đối với dao động tự dầm bậc nano FGM Các tác giả vị trí tới hạn tại đó sự thay đổi bậc có ảnh hưởng lớn tới một tần số nhất định, đặc biệt tần số không thứ nguyên bản, và sự thay đổi dầm bậc theo độ cao là rõ rệt dầm bậc theo bề rộng Đồng thời sự thay đổi số phần thể tích gây rõ rệt so với tham số không cục bộ, ảnh hưởng tham số không cục bộ đối với tần số cao lớn nhiều so với tần sớ thấp Ngồi ra, tác giả “nghịch lý không cục bộ” xuất hiện đối với tần số đầu tiên dầm công xôn Nghiên cứu này là bước đầu cho nghiên cứu kết cấu nano FGM khác dầm liên tục, khung nano nhiều bậc phức tạp Trong báo này, tác giả phát triển mơ hình DSM để phân tích dao động tự dầm bậc nano FGM dựa NET, gọi mơ hình DSM-NET Dầm bậc nano FGM mơ hình hóa theo lý thuyết dầm Timoshenko với phương trình chuyển động rút từ nguyên lý Hamilton Dùng DSM-NET, tác giả tìm tần số dạng dao động riêng kết cấu nano FGM dạng dầm, khung có xét đến vị trí thực trục trung hịa cho điều kiện biên khác So sánh kết tính tốn DSM-NET với kết công bố khẳng định độ tin cậy mơ hình Đây kết mới, nhóm tác giả phát triển thời gian gần cảmnhằm ơn Từ đó, tác giả tiến hành khảo sátLờisố đánh giá ảnh hưởng tham số phân bố Nghiên cứu này Trường Đại học Xây dựng tài trợ theo đề tài mã sớ 35-2020/KHXD-TĐ vật liệu, hình học, khơng cục điều kiện biên dao động tự dầm bậc nano FGM Các tác giả vị trí tới hạn thay đổi bậc có ảnh hưởng lớn tới tần số định, đặc biệt tần số không thứ nguyên bản, thay đổi dầm bậc theo độ cao rõ rệt 15 rõ rệt so với dầm bậc theo bề rộng Đồng thời thay đổi số phần thể tích gây tham số khơng cục bộ, ảnh hưởng tham số không cục tần số cao lớn nhiều so 62 Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng với tần số thấp Ngoài ra, tác giả “nghịch lý không cục bộ” xuất tần số dầm công xôn Nghiên cứu bước đầu cho nghiên cứu kết cấu nano FGM khác dầm liên tục, khung nano nhiều bậc phức tạp Lời cảm ơn Nghiên cứu Trường Đại học Xây dựng tài trợ theo đề tài mã số 35-2020/KHXD-TĐ Tài liệu tham khảo [1] Shen, H.-S (2016) Functionally graded materials: nonlinear analysis of plates and shells CRC press [2] Mahamood, R M., Akinlabi, E T (2017) Functionally Graded Materials Springer International Publishing [3] Eringen, A C (2002) Nonlocal continuum field theories Springer Science & Business Media [4] Karliˇci´c, D., Murmu, T., Adhikari, S., McCarthy, M (2015) Non-Local Structural Mechanics John Wiley & Sons, Inc [5] Polizzotto, C (2001) Nonlocal elasticity and related variational principles International Journal of Solids and Structures, 38(42-43):7359–7380 [6] Eltaher, M A., Khater, M E., Emam, S A (2016) A review on nonlocal elastic models for bending, buckling, vibrations, and wave propagation of nanoscale beams Applied Mathematical Modelling, 40 (5-6):4109–4128 [7] Salehipour, H., Shahidi, A R., Nahvi, H (2015) Modified nonlocal elasticity theory for functionally graded materials International Journal of Engineering Science, 90:44–57 [8] Reddy, J N (2007) Nonlocal theories for bending, buckling and vibration of beams International Journal of Engineering Science, 45(2-8):288–307 [9] Wang, C M., Zhang, Y Y., He, X Q (2007) Vibration of nonlocal Timoshenko beams Nanotechnology, 18(10):105401 [10] Li, C., Lim, C W., Yu, J L., Zeng, Q C (2011) Analytical solutions for vibration of simply supported nonlocal nanobeams with an axial force International Journal of Structural Stability and Dynamics, 11 (02):257–271 [11] Aydogdu, M (2009) A general nonlocal beam theory: Its application to nanobeam bending, buckling and vibration Physica E: Low-dimensional Systems and Nanostructures, 41(9):1651–1655 [12] Wang, C M., Kitipornchai, S., Lim, C W., Eisenberger, M (2008) Beam bending solutions based on nonlocal Timoshenko beam theory Journal of Engineering Mechanics, 134(6):475–481 [13] Chakraverty, S., Behera, L (2015) Free vibration of non-uniform nanobeams using Rayleigh–Ritz method Physica E: Low-dimensional Systems and Nanostructures, 67:38–46 [14] Eltaher, M A., Alshorbagy, A E., Mahmoud, F F (2013) Vibration analysis of Euler–Bernoulli nanobeams by using finite element method Applied Mathematical Modelling, 37(7):4787–4797 [15] Eltaher, M A., Mahmoud, F F., Assie, A E., Meletis, E I (2013) Coupling effects of nonlocal and surface energy on vibration analysis of nanobeams Applied Mathematics and Computation, 224:760– 774 [16] Pradhan, S C (2012) Nonlocal finite element analysis and small scale effects of CNTs with Timoshenko beam theory Finite Elements in Analysis and Design, 50:8–20 [17] de Sciarra, F M (2014) Finite element modelling of nonlocal beams Physica E: Low-dimensional Systems and Nanostructures, 59:144–149 [18] Tuna, M., Kirca, M (2017) Bending, buckling and free vibration analysis of Euler-Bernoulli nanobeams using Eringen’s nonlocal integral model via finite element method Composite Structures, 179:269–284 [19] Alotta, G., Failla, G., Zingales, M (2014) Finite element method for a nonlocal Timoshenko beam model Finite Elements in Analysis and Design, 89:77–92 [20] Ebrahimi, F., Nasirzadeh, P (2015) A nonlocal Timoshenko beam theory for vibration analysis of thick nanobeams using differential transform method Journal of Theoretical and Applied Mechanics, 1041 63 Đạt, P T., cs / Tạp chí Khoa học Cơng nghệ Xây dựng [21] Jena, S K., Chakraverty, S (2018) Free vibration analysis of variable cross-section single layered graphene nano-ribbons (SLGNRs) using differential quadrature method Frontiers in Built Environment, [22] S¸ims¸ek, M., Yurtcu, H H (2013) Analytical solutions for bending and buckling of functionally graded nanobeams based on the nonlocal Timoshenko beam theory Composite Structures, 97:378–386 [23] Rahmani, O., Pedram, O (2014) Analysis and modeling the size effect on vibration of functionally graded nanobeams based on nonlocal Timoshenko beam theory International Journal of Engineering Science, 77:55–70 [24] Mechab, I., Meiche, N E., Bernard, F (2016) Free vibration analysis of higher-order shear elasticity nanocomposite beams with consideration of nonlocal elasticity and Poisson effect Journal of Nanomechanics and Micromechanics, 6(3):04016006 [25] Uymaz, B (2013) Forced vibration analysis of functionally graded beams using nonlocal elasticity Composite Structures, 105:227–239 [26] Ebrahimi, F., Salari, E (2015) A semi-analytical method for vibrational and buckling analysis of functionally graded nanobeams considering the physical neutral axis position CMES: Comput Model Eng Sci, 105(2):151–181 [27] Narendar, S., Gopalakrishnan, S (2011) Spectral finite element formulation for nanorods via nonlocal continuum mechanics Journal of Applied Mechanics, 78(6) [28] Eltaher, M A., Alshorbagy, A E., Mahmoud, F F (2013) Determination of neutral axis position and its effect on natural frequencies of functionally graded macro/nanobeams Composite Structures, 99: 193–201 [29] Eltaher, M A., Alshorbagy, A E., Mahmoud, F F (2013) Determination of neutral axis position and its effect on natural frequencies of functionally graded macro/nanobeams Composite Structures, 99: 193–201 [30] Eltaher, M A., Abdelrahman, A A., Al-Nabawy, A., Khater, M., Mansour, A (2014) Vibration of nonlinear graduation of nano-Timoshenko beam considering the neutral axis position Applied Mathematics and Computation, 235:512–529 [31] Eltaher, M A., Khairy, A., Sadoun, A M., Omar, F.-A (2014) Static and buckling analysis of functionally graded Timoshenko nanobeams Applied Mathematics and Computation, 229:283–295 [32] Aria, A I., Friswell, M I (2019) A nonlocal finite element model for buckling and vibration of functionally graded nanobeams Composites Part B: Engineering, 166:233–246 [33] Trinh, L C., Vo, T P., Thai, H.-T., Nguyen, T.-K (2018) Size-dependent vibration of bi-directional functionally graded microbeams with arbitrary boundary conditions Composites Part B: Engineering, 134:225–245 [34] Su, H., Banerjee, J R (2015) Development of dynamic stiffness method for free vibration of functionally graded Timoshenko beams Computers & Structures, 147:107–116 [35] Lien, T V., Duc, N T., Khiem, N T (2019) Free and forced vibration analysis of multiple cracked FGM multi span continuous beams using dynamic stiffness method Latin American Journal of Solids and Structures, 16(2) [36] Lien, T V., Khiem, N T., Duc, N T (2016) Free vibration analysis of functionally graded Timoshenko beam using dynamic stiffness method Journal of Science and Technology in Civil Engineering (STCE) NUCE, 10(5):19–28 [37] Adhikari, S., Murmu, T., McCarthy, M A (2013) Dynamic finite element analysis of axially vibrating nonlocal rods Finite Elements in Analysis and Design, 63:42–50 [38] Taima, M S., El-Sayed, T A., Farghaly, S H (2020) Free vibration analysis of multistepped nonlocal Bernoulli–Euler beams using dynamic stiffness matrix method Journal of Vibration and Control, 27 (7-8):774–789 [39] Li, X.-F., Wang, B.-L (2009) Vibrational modes of Timoshenko beams at small scales Applied Physics Letters, 94(10):101903 [40] Ghavanloo, E., Rafii-Tabar, H., Fazelzadeh, S A (2019) Computational Continuum Mechanics of Nanoscopic Structures Springer International Publishing 64 ... toán dao động dầm nano FGM theo lý thuyết đàn hồi không cục Trong báo này, tác giả phát triển DSM để phân tích dao động tự dầm bậc nano FGM sở lý thuyết dầm Timoshenko lý thuyết đàn hồi không cục. .. giải tích xác cho tốn dao động dầm nano tìm chúng cho phép nghiên cứu dao động xác dầm dải tần số tùy ý Adhikari cs [4, 37] thu ma trận độ cứng động lực nano đồng chịu dao động dọc Hàm đáp ứng tần... ổn định dầm Timoshenko FGM phương pháp giải tích Ngồi ra, Mechab cs [24] nghiên cứu dao động tự do, Uymaz [25] nghiên cứu dao động cưỡng dầm nano, hai sử dụng lý thuyết biến dạng cắt bậc cao