Một số kỹ thuật đối ngẫu ứng dụng trong phân tích dao động và ổn định

Nội dung text: Một số kỹ thuật đối ngẫu ứng dụng trong phân tích dao động và ổn định

1. Tuyển tập Hội nghị khoa học toàn quốc lần thứ nhất về Động lực học và Điều khiển Đà Nẵng, ngày 19-20/7/2019, tr. 238-243, DOI 10.15625/vap.2019000284 Một số kỹ thuật đối ngẫu ứng dụng trong phân tích dao động và ổn định Nguyễn Đông Anh(1,2,3*), Nguyễn Cao Thắng(1,2) (1) Viện Cơ học, Viện Hàn lâm Khoa học và Công nghệ Việt Nam (2) Trường Đại học Công nghệ, Đại học Quốc gia Hà Nội (3) Học viện Khoa học và Công nghệ, Viện Hàn lâm Khoa học và Công nghệ Việt Nam tuyến tính hóa k và hệ số trở về như sau 1. Mở đầu dn dn ABAB Thiên nhiên và cuộc sống luôn chứa đựng các kppkdndndndn1; (4) khuynh hướng đối ngẫu nhau. Đó là các yếu tố có tính BA22 chất trái ngược nhau, hoặc bổ trợ cho nhau. Do vậy, việc Giải hệ phương trình (4) được nghiên cứu khoa học cũng cần phản ánh được tính chất này. Gần đây, cách tiếp cận đối ngẫu đã được đề xuất để 1 p AB 1 p k , r 2 (5) nghiên cứu đáp ứng của các hệ phi tuyến [1,2] và đã dn 2 dn 2 1 rp B2 1 rp được phát triển trong nhiều công trình, ví dụ [3-17]. Một ưu điểm quan trọng của cách tiếp cận đối ngẫu đối với trong đó r là đại lượng đặc trưng cho mức độ phụ một vấn đề khoa học là luôn xem xét hai khía cạnh khác thuộc của A và B, xác định theo công thức nhau (đối ngẫu) của vấn đề; điều này cho phép việc nghiên cứu trở nên hài hòa và phản ánh được thực chất AB hơn. Bài báo giới thiệu tóm tắt một số kỹ thuật đối ngẫu r (6) ứng dụng trong phân tích dao động và ổn định. Một số ví AB22 dụ minh họa về kỹ thuật đối ngẫu được giới thiệu và bình luận nhằm đánh giá các nghiên cứu phát triển tiếp theo. Trong bài báo [17] giá trị trọng số p được chọn theo công thức 2. Thay thế tương đương đối ngẫu 1 r 2 p (7) Trong Cơ học ta thường hay thay gần đúng đối tượng 1 r 2 A được mô tả bằng hàm A(x) bằng đối tượng B được mô Có thể nhận thấy từ (6) và (7) giá trị của trọng số p tả bằng hàm B(x). Như vậy, ta sẽ thay A(x) xấp xỉ bằng phụ thuộc vào đại lượngr 2 mà là chỉ số cho mức độ khác kttB(x) với hệ số thay thế tương đương ktt được xác định biệt giữa A và B. Thay (7) vào (5) sẽ cho giá trị của hệ số bằng tiêu chuẩn bình phương tối thiểu: thay thế tương đương 2 SA xk B( x )( )min (1) 2 AB tttt k 2r tt k (8) w 1 r 4 2 trong đó chỉ số tt viết tắt của chữ thông thường, là B phép lấy trung bình tiền định hoặc trung bình xác suất. Như vậy sử dụng tiêu chuẩn thay thế đối ngẫu (3) với Trong trường hợp B(x) là hàm tuyến tính hệ số thay thế giá trị trọng số xác định theo (7), (6) hàm A(x) được thay tương đương k thường được gọi là hệ số tuyến tính hóa tt thế tương đương bằng kdnB(x) trong đó hệ số thay thế tương đương. Với cách tiếp cận đối ngẫu cho bài toán tương đương kdn xác định theo (8). Với cách thay thế thay thế, ta xét dạng mở rộng của tiêu chuẩn (1) thành trên ta không phải làm việc với hàm A(x) nữa mà chỉ phải tiêu chuẩn đối ngẫu như sau làm việc với hàm B(x). Đặc biệt khi khi B(x) là hàm 2 SpA xk1( B )( x ) tuyến tính của x ta sẽ đưa bài toán phi tuyến về bài toán dndn tuyến tính dễ giải hơn nhiều. (2) 2 p kdndn B( xA )( x )min kdn, dn 3. Thay thế tương đương địa phương-toàn thể Năm 2012, dựa trên cách tiếp cận đối ngẫu, N. D. trong đó chỉ số dn ký hiệu đối ngẫu, dn là hệ số trở về, p là trọng số của sự thay thế trở về được chuẩn hóa Anh, L.X. Hung và L. D. Viet [5] đã phát triển tiêu chuẩn sai số bình phương địa phương – tổng thể cho các hệ 01p (3) ngẫu nhiên phi tuyến bằng cách kết hợp hai phạm trù địa Với p=0 tiêu chuẩn (3) trở thành tiêu chuẩn (2). Giả phương và tổng thể. Giả sử ta thay hàm ngẫu nhiên A(x) sử cho trước trọng số p, từ điều kiện cực tiểu (2) dẫn tới xấp xỉ bằng kB(x) với hệ số thay thế k được xác định các phương trình xác định hệ số thay thế tương đương bằng tiêu chuẩn bình phương tối thiểu:
2. Nguyễn Đông Anh, Nguyễn Cao Thắng 2 1 b SA xkBx()()min (9) [[()]][()]yxyxdr k srbs s (17) Viết chi tiết (7) ta có 11br yxdxdr() 2 bsra sa A( x ) kB ( x ) P ( x ) dx min (10) k Ví dụ với y(x) = xn và a=0, b=1 ta có r r n x n1 n 1 1 n 1 r Trong đó Px()là hàm mật độ xác suất (PDF) của . []x x dx x (18) r 0 Suy ra: r0 r( n 1) n 1 A( x ) B ( ) x Giá trị trung bình địa phương-toàn thể tại lát cắt s k (11) của hàm xn sẽ là 2 Bx() 1 Do khoảng tích phân trong (10) là ( , ), tiêu 111 r n [[)]][()]xyxdrdrn chuẩn (10) có thể được gọi là tiêu chuẩn bình phương tối sr111ssn s s (19) thiểu toàn thể. Với giả thiết cho rằng phép lấy tích phân 1 1 (1) sn 1 cần tập trung hơn để cho nghiệm chính xác hơn, Anh và r n 1 22s Di Paola vào năm 1995 đề nghị tiêu chuẩn bình phương (1)nsns (1)(1) (1) tối thiểu địa phương Ta thấy rằng giá trị trung bình địa phương-toàn thể n r x tại lát cắt s của hàm x là một hàm của s, và khi s = 1 giá [(()())A xkB ](()())()min xA xkB22 xP x dx tr ị này sẽ b ằ ng giá tr ị trung bình thông thường của xn k r x 1 [[]]xxnn (20) (12) 1 (1)n với r là một số dương nào đó, x là độ lệch chuẩn của x Như vậy giá trị trung bình địa phương-toàn thể tại lát và [.] ký hiệu giá trị trung bình xác suất địa phương cắt s là sự mở rộng của giá trị trung bình thông thường và r x chứa giá trị trung bình thông thường như trường hợp [.] (.)P ( x ) dx (13) riêng. r x Tương tự (11) từ (12) ta có 5. Giá trị trung bình trọng số đối ngẫu A()() xBx Một sự phát triển khác của cách lấy trung bình thông kr() (14) thường là cách lấy trung bình trọng số. Ta xét hàm x(t) 2 Bx() trên miền từ 0 đến vô cùng. Giá trị trung bình trọng số của x(t) được xác định như sau Ta thấy từ (12) và (13) hệ số k sẽ là hàm số của r. Sử dụng quan điểm đối ngẫu ta có thể chọn k bằng giá trị Wx(( th ))( t )(x t ) dt (21) trung bình toàn thể như sau [5]: 0 trong đó h(t) thỏa mãn điều kiện chuẩn hóa 1 s k k()() r Lim k r dr (15) d s htdt()1 (22) s 0 0 trong đó là ký hiệu giá trị trung bình thông thường Hàm số h(t) gọi là hàm trọng số. Phép lấy trung bình cho hàm số. Ta thu được một tiêu chuẩn thay thế tương trọng số được sử dụng nhiều trong toán học cũng như đương mới gọi là tiêu chuẩn thay thế địa phương – toàn trong cuộc sống thường ngày. Sự hiệu quả của nó phụ thể. Theo tiêu chuẩn này hàm ngẫu nhiên A(x) được thay thuộc rất nhiều vào cách chọn hàm trọng số. Trong bài tương đương bằng hàm kB(x) với hệ số thay thế k được báo này hàm trọng số đối ngẫu sau đây được giới thiệu xác định theo (13), (14), (15). [10] h( tstes ); ,0 t22 st (23) 4. Giá trị trung bình địa phương-toàn thể Như ta thấy hàm trọng số đối ngẫu (23) thu được Giả sử y(x) là hàm khả tích của xє[a, b]. Ta có thể bằng cách lấy tích của hàm đồng biến t với hàm nghịch đưa ra các giá trị trung bình địa phương-toàn thể của y(x) biến e st. Kết quả là hàm trọng số đối ngẫu (23) vừa như sau. Ta xác định giá trị trung bình địa phương tại lát đồng biến vừa nghịch biến (xem Hình 1). cắt r của hàm y(x) với r є [a,b] theo công thức 1 r [y ( x )]r y ( x ) dx (16) raa Giá trị trung bình địa phương-toàn thể tại lát cắt s của hàm y(x) với s є [a,b] theo công thức
3. Một số kỹ thuật đối ngẫu ứng dụng trong phân tích dao động và ổn định Hình 1. Hàm trọng số đối ngẫu (23) Hình 3. Plots of c osn versus n for s=2 Dựa trên hàm trọng số đối ngẫu (23) một giá trị trung s bình trọng số mới được đề nghị áp dụng cho hàm số ω-tuần hoàn x(ωt) như sau Wxtxtstext(())()() dt 22 st ss0 (24) sexd2 s () 0 Ta sử dụng hai ký hiệu cho giá trị trung bình trọng số (24) để dùng chúng trong các trường hợp tùy theo sự tiện lợi. Chỉ số s dùng để phân biệt với giá trị trung bình thông thường thu được khi cho s→0. Sử dụng phép biến đổi Laplace ta có Wn(cos)cos() tsten t dt 22 st Hình 4. Plots of s inn versus s for n=2 s 0 s sn22 (25) sneds22cos() s 0 ()sn22 2 W(())sin() sin n tsten t dt 22 st s 0 2sn (26) sneds22sin() s 0 ()sn22 2 Đồ thị của cosn s theo s và n được cho trên các Hình 2-3. Ta thấy rằng giá trị trung bình trọng số (25) có một giá trị cực tiểu tại 0<s<n và giá trị cực đại bằng 1 khi s → ∞. Đối với các hàm điều hòa bậc cao Hình 5. Plots of sinn s versus n for s=2 W(cos n t ) , W(sin n t ) tiến tới 0 khi n → ∞. Ta lưu s s Sử dụng (25) ta có ý rằng giá trị trung bình trọng số (25) Wnts (cos) bằng 2 1 cos2nt Wnss(cos)() t W 0 tại hai giá trị của s, cụ thể s=0 and s=n. Đồ thị của 2 1cos21 14 n tsn 22 sinn theo s và n được cho trên các Hình 4-5. WWs( )(), 2 s ss222 2 (4sn22 2) (27) 1 cos2nt Wn(sin)() t2 W ss2 1cos21 14 n tsn 22 WWs( )(). 2 ss222 2 (4sn22 2) Ta có thể khai triển hàm ω-tuần hoàn x(ωt) theo chuỗi Fourier x( txxi )( t cossin)0 xi t icis (28) i 1 và sử dụng (25) và (26) để thu được Hình 2. Plots of cosn s versus s for n=2
4. Nguyễn Đông Anh, Nguyễn Cao Thắng ta tính WxtWxxitxitssicis(())((cossin)) 0 i 1 3 Axx22622 15, xxWitx0 Wit((cos)(sin))icsiss (29) i 1 2 42 siis22 2 ABxx 3, (34) xsxx 2 () 0 icis 222222 3 i 1 ()()sisi Bxr222 , Từ (29) thu được giá trị trung bình thông thường của 5 hàm x(ωt) bằng cách cho s=0 Sử dụng (34) để tính giá trị trọng số p theo (7) và hệ số thay thế tương đương theo (8) sẽ cho xtWxtx()(()) (30) 00 p 1/ 4 (35) So sánh (29) với (30) ta thấy rằng giá trị trung bình 45 trọng số có chứa các thông tin về các thành phần điều hòa 2 kxw (36) bậc cao của hàm ω-tuần hoàn x(ωt) (cụ thể là các hệ số 17 xic và xis) trong khi giá trị trung bình thông thường không Thay (36) vào (33) được phương trình xác định chứa các thong tin này. Từ đó có thể nhận dịnh rằng giá momen đáp ứng bậc 2 của hệ Duffing (31) trị trung bình trọng số có thể cho các kết quả chính xác 2 2 45 222 hơn về hàm x(ωt). xx0 0 174 w w h 6. Một số ví dụ ứng dụng hay 1745 2 Có nhiều ví dụ khác nhau về ứng dụng của các kỹ x 224 (37) thuật đối ngẫu trình bày ở trên. Trước hết kỹ thuật thay w 9017 00 h thế tương đương đối ngẫu đã được áp dụng để phân tích dao động ngẫu nhiên phi tuyến cho hệ cơ học 1 bậc tự do Để áp dụng kỹ thuật thay thế kinh điển (thông 3 và nhiều bậc tự do trong [3,4,17]; áp dụng để tính toán thường) ta lấy cực tiểu sai số thay thế x bằng kclx các tham số tối ưu của bộ giảm chấn TMD cho hệ kết cấu ()minxkx32 (38) cl k có cản [6,8]; áp dụng để phân tích dao động tuần hoàn cl (flutter) của mô hình dao động uốn xoắn cho thiết diện 2 Suy ra kx3 (39) cánh máy bay [11]; áp dụng để phân tích bài toán nhiệt vệ tinh quỹ đạo thấp [15,16]. Kỹ thuật trung bình trọng Thay (39) vào (33) tìm được phương trình xác định số đối ngẫu được áp dụng để phân tích tần số dao động momen đáp ứng bậc 2 x 2 của hệ Duffing (31) theo kỹ phi tuyến cho nhiều hệ khác nhau [14,18,19], cũng như cl cho kết cấu dầm micro và nano [20,21]. Kỹ thuật thay thế thuật thay thế tương đương kinh điển tương đương dịa phương-toàn thể được áp dụng để phân 2 2 30xx222 (40) tích dao động ngẫu nhiên phi tuyến cho hệ cơ học 1 và GLGL o 4h nhiều bậc tự do trong [5,9]; áp dụng để tính toán các So sánh sai số của các mômen đáp ứng x 2 tính theo tham số tối ưu của bộ giảm chấn TMD cho hệ kết cấu có w cản [12]; Kỹ thuật tuyến tính hóa điều chỉnh-một dạng kỹ thuật thay thế đối ngẫu và x 2 tính theo kỹ thuật của kỹ thuật đối ngẫu-được áp dụng để phân tích dao cl động ngẫu nhiên phi tuyến cho kết cấu dầm [7,13]. Trong thay thế thông thường của hệ Duffing với bài báo này chúng ta chỉ xét minh họa hệ sau đây. h 0.5 ; 1; 2 và γ thay đổi được cho trong Dao động phi tuyến trong hệ Duffing 0 Ta xét hệ sau Bảng 1, trong đó x 2 mômen đáp ứng chính xác cx 23 x2 hx0 x x t (31) [xem 17]. Ta thấy rằng sai số nghiệm đã được giảm khá nhiều khi tính phi tuyến của hệ tăng dần. trong đó h,,, là các hằng số dương, t là quá 0 Bảng 1. So sánh các lời giải gần đúng của hệ Duffing trình ngẫu nhiên ồn trắng có trung bình bằng không và (kỹ thuật thay thế đối ngẫu) cường độ đơn vị. Sai 2 2 Sai số 2 Đặt A x3, B x , ta sẽ thay Axk Bk3 x x x x số ww cx cl (%) w (%)) và đưa phương trình phi tuyến (31) về phương trình tuyến tính 0.01 0.9721 0.9717 0.05 0.9748 0.28 0.05 0.889 0.883 0.64 0.894 0.62 xhxk2() xt 2 (32) 0w 0.1 0.818 0.805 1.49 0.821 0.47 Từ đó sẽ thu được phương trình xác định mômen đáp 0.5 0.579 0.549 5.29 0.570 1.59 ứng bậc 2 [17] 1.0 0.468 0.434 7.19 0.454 2.95 2 2 2 2 5.0 0.254 0.227 10.74 0.240 5.76 kw0 x x 0 (33) 4h 10.0 0.189 0.167 11.77 0.176 6.62 Để áp dụng kỹ thuật thay thế tương đương đối ngẫu
5. Một số kỹ thuật đối ngẫu ứng dụng trong phân tích dao động và ổn định Kỹ thuật thay thế tương đương địa phương-toàn thể momen đáp ứng bậc 2 x 2 của hệ Duffing (31) như dựa trên giả thiết rằng các giá trị của dịch chuyển và vận GL tốc của hệ phi tuyến chỉ tập trung trong một miền hữu sau theo kỹ thuật thay thế tương đương địa phương-toàn hạn. Giả thiết này có thể được minh họa bằng đồ thị của thể hàm phân bố xác suất của hệ Dufing (26) được cho trên 2 2 2.41190 xx222 (46) hình 6. GLGL o 4h So sánh sai số của các mômen đáp ứng x 2 tính GL theo kỹ thuật thay thế địa phương-toàn thể và x 2 tính cl theo kỹ thuật thay thế thông thường của hệ Duffing với o 1,0.25,1h ; hệ số đàn hồi phi tuyến 0.1 2 0.075p thay đổi được cho trong Bảng 2, trong đó x cx 0.05 2 0.025 mômen đáp ứng chính xác [xem 5]. Ta cũng thấy rằng sai 0 số nghiệm đã được giảm khá nhiều khi tính phi tuyến của 0 x -2 hệ tăng dần. 0 Bảng 2. So sánh các lời giải gần đúng của hệ Duffing -2 x (kỹ thuật thay thế địa phương-toàn thể) 2 2 2 2 x x sai số (%) x sai số (%) Hình 6. Đồ thị hàm PDF p x x, , γ =1.0 cx kd GL 0.1 0.8176 0.8054 -1.490 0.8327 1.857 Để xác định k theo kỹ thuật thay thế tương đương địa phương-toàn thể trước hết ta tính cực tiểu địa phương sau 1.0 0.4680 0.4343 -7.194 0.4692 0.263 10 0.1889 0.1667 -11.768 0.1839 -2.626 ()minxkx32 (41) k 100 0.2543 0.2270 -10.740 0.2624 3.150 trong đó [.] ký hiệu giá trị trung bình xác suất địa phương r x [.]()() Pxdx (42) 7. Kết luận r Thiên nhiên và cuộc sống luôn chứa đựng các x khuynh hướng đối ngẫu nhau. Đó là các yếu tố có tính Suy ra chất trái ngược nhau, hoặc bổ trợ cho nhau. Do vậy, việc r x nghiên cứu khoa học cũng cần phản ánh được tính chất x4 P() x dx 4 này. Gần đây, cách tiếp cận đối ngẫu đã được đề xuất để x r kr() x nghiên cứu đáp ứng của các hệ phi tuyến [1,2] và đã r x 2 x được phát triển trong nhiều công trình, ví dụ [3-17]. Một 2 (43) x P() x dx ưu điểm quan trọng của cách tiếp cận đối ngẫu đối với r x một vấn đề khoa học là luôn xem xét hai khía cạnh khác 2TxT 22 2,2,rrx 2 nhau (đối ngẫu) của vấn đề; điều này cho phép việc 2 T nghiên cứu trở nên hài hòa và phản ánh được thực chất 2Tx1,r 1,r hơn. Bài báo giới thiệu tóm tắt một số kỹ thuật đối ngẫu trong đó ký hiệu ứng dụng trong phân tích dao động và ổn định. Một ví dụ 2 r 1 x minh họa về kỹ thuật đối ngẫu được giới thiệu cho thấy Txx dxxe 2n ( ),( ) 2 (44) nr, cách tiếp cận đối ngẫu là một cách tiếp cận hiệu quả 2 0 trong nghiên cứu khoa học có thể tiếp tục nghiên cứu Tiếp theo hệ số thay thế tương đương k tính theo kỹ phát triển cho nhiều bài toán khác. thuật thay thế tương đương địa phương-toàn thể bằng giá trị trung bình cộng của k(r). Như vậy ta có Lời cám ơn. Bài báo được sự tài trợ của Quỹ Phát triển Khoa học và Công nghệ Quốc gia (NAFOSTED) 1 s k k()() r Lim k r dr mã số 107.04-2018.12. và của Viện Hàn lâm KHCNVN s cho NCVCC mã số NCVCC 03.07/19-19. s 0 s 1 T Lim2,r x2 dr (45) s sT 0 1,r Tài liệu trích dẫn 1 s T x22 Lim2,r dr2.4119 x 1. N.D. Anh, Duality in the analysis of responses to s nonlinear systems. Vietnam J. Mech. 32(4) (2010) sT0 1,r 263–266. Thay (45) vào (33) tìm được phương trình xác định
6. Nguyễn Đông Anh, Nguyễn Cao Thắng 2. ND. Anh, Dual approach to averaged values of nonlinear analysis of thermal behavior of a two-node functions. Vietnam J. Mech. 34(3), 211–214 (2012). model for small satellites in Low Earth Orbit, 3. N.D. Anh, N.N. Hieu, N.N. Linh, A dual criterion of International Journal of Mechanical Sciences, 133 equivalent linearization method for nonlinear systems (2017), pp 513-523. subjected to random excitation, Acta. Mech.(3) (2012) 17. N. D. Anh, N. N. Linh, A weighted dual criterion of 645–654. the equivalent linearization method for nonlinear 4. Anh, N.D., Zakovorotny, V.L., Hieu, N.N., Diep, systems subjected to random excitation, Acta. Mech., D.V.: A dual criterion of stochastic linearization 2018, Volume 229, Issue 3, pp 1297–1310. method for multi-degreeof-freedom systems subjected 18. D. V. Hieu and N. Q. Hai. Analyzing of Nonlinear to random excitation. Acta Mec., 2667–2684. Generalized Duffing Oscillators Using the Equivalent 5. N.D. Anh, L.X. Hung, L.D. Viet, Dual approach to Linearization Method with a Weighted Averaging. local mean square error criterion for stochastic Asian Research Journal of Mathematics, vol. 9(1): equivalent linearization. Acta Mech. 224, 241–253 1-14, 2018. (2013). 19. D. V. Hieu, N. Q. Hai, and D. T. Hung. The 6. ND Anh, NX Nguyen, Design of TMD for damped Equivalent Linearization Method with a Weighted linear structures using the dual criterion of equivalent Averaging for Solving Undamped Nonlinear linearization method, International Journal of Oscillators. Journal of Applied Mathematics, Volume Mechanical Sciences, 2013, Volume 77, December 2018, Article ID 7487851, 15 pages. 2013, Pages 164-170. 20. D. V. Hieu, N. D. Anh, L.M. Quy, N. Q. Hai, 7. N.D. Anh, I. Elishakoff, N.N. Hieu Extension of the Nonlinear vibration of microbeams based on the regulated stochastic linearization to beam vibrations, nonlinear elastic foundation using the equivalent Prob. Eng. Mech 35(2014), 2–10. linearization methodwith a weighted averaging, 8. N.D. Anh and N.X. Nguyen, Design of Archive of Applied Mechanics, 2019 (under reviewer). non-traditional dynamic vibration absorber for 21. D. V. Hieu, N. D. Anh, L.M. Quy, D. T. Hung, damped linear structures, J. Mechanical Engineering Nonlinear vibration of nanobeams under electrostatic Science, vol. 228(1), pp.45-55, 2014. force based on the nonlocal strain gradient theory, 9. N.D. Anh, L.X. Hung, L.D. Viet, N.C. Thang, International Journal of Mechanics and Materials in Global-local mean square error criterion for Design, 2019 (accepted for publication). equivalent linearization of nonlinear systems under random excitation, Acta. Mech., 226(9) (2015) 3011-3029. 10. N. D. Anh, Dual approach to averaged values of functions: A form for weighting coefficient, Vietnam J. Mech., 37(2) (2015) 145–150. 11. Triet, N.M., Extension of dual equivalent linearization technique to flutter analysis of two dimensional nonlinear airfoils. Vietnam J. Mech. 37(3), 217–230 (2015) 12. N.D.Anh, N.X.Nguyen, N.H.Quan, Global-local approach to the design of dynamic vibration absorber for damped structures, Journal of Vibration and Control, vol. 22(14), pp. 3182-3201,2016 13. N.D. Anh, I. Elishakoff, N.N. Hieu, Generalization of Seide’s problem by the regulated stochastic linearization technique, Meccanica, 2017 (52), pp1003-1016. 14. N. D. Anh, N. Q. Hai, and D. V. Hieu. The equivalent linearization method with a weighted averaging for analyzing of nonlinear vibrating systems. Latin American J. Solids and Structures, 14(9) (2017) 1723–1740. 15. N.D. Anh, NN Hieu, PN Chung, NT Anh Thermal radiation analysis for small satellites with single-node model using techniques of equivalent linearization, Applied Thermal Engineering, 2016, Volume 94, Pages 607-614. 16. P. N. Chung, N. D. Anh, N. N. Hieu, D. V. Manh, Extension of dual equivalent linearization to