Phương pháp Proper Generalized Decomposition cho bài toán tấm mỏng chịu uốn

Bài viết trình bày phương pháp Proper Generalized Decomposition (PGD) để giải quyết bài toán tấm mỏng chịu uốn trong không gian hai chiều. Phương pháp PGD được áp dụng để đưa bài toán hai chiều thành chuỗi các bài toán một chiều.

Tạp chí Khoa học Đại học Thủ Dầu Một Số 1(32)-2017 PHƯƠNG PHÁP PROPER GENERALIZED DECOMPOSITION CHO BÀI TOÁN TẤM MỎNG CHỊU UỐN Lê Quốc Cường(1), Nguyễn Bá Duy(2) (1) Trường Đại học Sư phạm Kỹ thuật TP.HCM; (2) Trường Đại học Thủ Dầu Một Ngày nhận 29/12/2016; Chấp nhận đăng 29/01/2017; Email: lecuong2109@gmail.com Tóm tắt Trong báo này, chúng tơi trình bày phương pháp Proper Generalized Decomposition (PGD) để giải toán mỏng chịu uốn không gian hai chiều Phương pháp PGD áp dụng để đưa toán hai chiều thành chuỗi tốn chiều Sau đó, tốn chiều giải phương pháp sai phân hữu hạn Kết mơ số áp dụng cho tốn mỏng chịu uốn với điều kiện biên khác Các kết tính tốn so sánh với lời giải giải tích Từ khóa: giảm bậc mơ hình, Proper Generalized Decomposition, mỏng chịu uốn Abstract PROPER GENERALIZED DECOMPOSITION METHOD FOR THE THIN PLATE BENDING PROBLEM In this paper, we present Proper Generalized Decomposition (PGD) method to solve the problem of thin plate bending in two-dimensional space PGD method is applied to transform the two-dimensional problem into a series of one-dimensional problems Then, each onedimensional problem is solved by the finite difference method Numerical simulation results are applied to thin plate bending problem with different boundary conditions The calculation results are compared with analytical solutions Giới thiệu Nhiều mơ hình tốn thường gặp khoa học kỹ thuật thường định nghĩa không gian đa chiều, điều làm cho vấn đề chiều thứ nguyên trở nên phức tạp áp dụng kỹ thuật chia lưới rời rạc thông thường Hơn mơ hình theo tiêu chuẩn trở thành đa chiều thơng số thay đổi Vì việc phát triển phương pháp nhằm giải tốn cách nhanh chóng cần thiết Phương pháp PGD lần giới thiệu giáo sư Chinesta cộng [1] Sự đời phương pháp PGD góp phần hỗ trợ giải tốn có số chiều khơng gian lớn cách hiệu với thời gian xử lý nhanh độ xác cao Phương pháp PGD ngày mở rộng ứng dụng để giải toán đa chiều lĩnh vực lưu chất [2], truyền nhiệt [3], vật liệu composite [4] Phương pháp PGD phương pháp giảm bậc mô hình dựa sở tách biến, lời giải tốn tìm dạng tổng tích hàm số chiều không gian 63 Lê Quốc Cường Phương pháp proper generalized decomposition cho toán Giả sử trường u phụ thuộc N biến số  x1 , x2 , , xN  , giá trị u viết dạng tách biến sau: Q u( x1 , x2 , , xN )   Fi ( x1 ) Fi ( x2 ) Fi ( xN ) (1) i 1 xi  i  1, 2, , N  biến số không gian, thời gian hay tham số mà tốn cần khảo sát Bài tốn phân tích mỏng chịu uốn thực thành công với nhiều phương pháp số khác (phần tử hữu hạn, sai phân hữu hạn, phương pháp phổ, …) Trong báo này, phương pháp PGD đề xuất để giải toán mỏng chịu uốn Phương pháp PGD áp dụng để đưa phương trình vi phân đạo hàm riêng toán mỏng chịu uốn khơng gian hai chiều thành chuỗi phương trình vi phân khơng gian chiều Sau đó, phương pháp sai phân hữu hạn dựa sơ đồ sai phân trung tâm bậc hai áp dụng để giải phương trình chiều Bài báo tổ chức sau, phần trình bày phương trình vi phân chủ đạo toán mỏng chịu uốn Phần trình bày phương pháp PGD cho phương trình biharmonic không gian hai chiều Sau cùng, kết mơ trình bày phần Phương trình vi phân chủ đạo cho tốn mỏng chịu uốn Phương trình vi phân chủ đạo tốn mỏng chịu uốn khơng gian hai chiều có dạng phương trình biharmonic [7] sau 4w 4w  w p , (2) w độ võng tấm, p lực tác dụng lên bề mặt    x x 2y y D D  Eh 12 1    độ cứng chịu uốn Trong đó, E mơ đun đàn hồi, h chiều dày  hệ số poisson Phương pháp PGD cho phương trình biharmonic Xét phương trình biharmonic khơng gian hai chiều sau  4u  4u  4u    f  x, y  miền    x   y (3) x x 2y y Mục tiêu áp dụng phương pháp PGD để tìm nghiệm xấp xỉ phương trình (3) Giả sử nghiệm xấp xỉ phương trình viết dạng tách biến sau N u  x, y    X i  x   Yi  y  (4) i 1 Giả sử lời giải bước lặp thứ n biết, cần tìm lời giải bước lặp thứ n  u n 1 n  x, y    X i  x   Yi  y   R  x   S  y  (5) đây: i 1 Phương trình (3) đưa dạng yếu sau  4u  4u *  u u      x x 2y y   x  y R  x   X n1  x  S  y   Yn1  y   f  dxdy  ,  với u*  R*  S  R  S * hàm trọng số 64 (6) (7) Tạp chí Khoa học Đại học Thủ Dầu Một Số 1(32)-2017 Thay (5) (7) vào phương trình (6), ta  R  x  y *  d 4R d 2R d 2S d 4S   S  R  S *     S  2   R   dxdy    R *  S  R  S *   fdxdy dx dy dy   dx  x  y  R  *  x  y 2  R n  d 4Xi   S  R  S*     Yi dxdy dx i 1   *  x  y   R  x  y * (8) n  d X i d 2Yi   S  R  S*     dxdy dy  i 1  dx n  d 4Y   S  R  S *     X i  4i dxdy dy  i 1  Để giải phương trình (8) tìm R  x  S  y  , sử dụng giải thuật lặp cố định luân phiên gồm bước sau: Bước 1: Tìm hàm R  x  Giả sử S  y  biết, S *  , thay vào phương trình (8) ta  d 4R d 2R d 2S d 4S  * R  S   S    R  dxdy   R*  S  fdxdy  2   dx dx dy dy    x  y  x  y  n  d 4Xi  * R  S   Yi dxdy    i 1  dx   x  y 2  (9) n  d X i d 2Yi  * R  S   dxdy    dy  i 1  dx  x  y n  d 4Yi  * R  S  X    i dxdy  dy  i 1   x  y Vì tất hàm phụ thuộc y phương trình (9) biết, thực tích phân chiều  y   a y   S dy  y  b  S  d S dy y   dy y   d 4S  c y   S  dy dy  y    f y  x    S  fdy y   i  a y   S  Yi dy y   byi  S  d Yi dy  dy  y   i d 4Yi c  S  y   dy dy y   (10) 65 Lê Quốc Cường Phương pháp proper generalized decomposition cho tốn Khi phương trình (9) trở thành n  d 4R d 2R *  * * i d Xi R  a  b  c R dx  R  f x dx  R  a      y y y y y   dx   i1 dx dx dx  x x x n   R*   byi x   x i 1 d 2Xi dx (11) dx n R*   ciy X i dx i 1 Phương trình (11) dạng yếu chiều định nghĩa  x Ngoài đưa dạng mạnh sau n n n d 4R d 2R i d Xi i d Xi a y  2by  c y R  f y  x    a y  b  ciy X i (12)   y dx dx dx dx i 1 i 1 i 1 Bước 2: Tìm hàm S  y  Với R  x  vừa tính bước trên, R*  , tiến hành tương tự bước tìm hàm R  x  , ta phương trình dạng mạnh hàm S  y  sau n  d 4S d 2S *  * * i d Yi S  a  b  c S dy  S  f y dy  S  a dy   x x   x dy x dy x     dy i 1 y y y n   S *   bxi y   y i 1 n S *   cxi Yi dy i 1  a x   R dx  x  b  R  d R dx  x  dx x   d 4R c  R  dx  x  dx  x    f x  y    R  fdx x   i a x   R  X i dx x   d 2Xi bxi   R  dx dX  x  c i  R  d X i dx  x  dx x  (14) 66 d 2Yi dy dY (13) Tạp chí Khoa học Đại học Thủ Dầu Một Số 1(32)-2017 Phương trình (13) đưa dạng mạnh sau n n n d 4Yi d 4S d 2S i d Yi a x  2bx  cx S  f x  y    a xi  b   cxi Yi (15)  x dy dy dy dy i 1 i 1 i 1 Các bước giải phương trình (12) phương trình (15) để tìm R  x  S  y  lặp q q 1 kết hội tụ Nếu kí hiệu R    x  R    x  hàm R  x  tính q q 1 bước lặp bước lặp trước, tương tự với S    y  S    y  , tiêu chuẩn dừng chọn sau e  R q  x   S q  y   Rq1  x   S q1  y    RS (16)  RS số chọn đủ bé để đảm bảo độ xác Sau bước lặp tìm R  x  S  y  hội tụ, xác định X n1  x   R  x  Yn1  y   S  y  Quá trình tìm cặp hàm  X i  x  , Yi  y   phải tiếp tục đạt hội tụ tồn cục tốn bước lặp thứ N , nghiệm xấp xỉ tốn tính sau N u  x, y    X i  x   Yi  y  (17) i 1 Điều kiện dừng toàn cục tốn tính sau res E   u (18) f  x, y   u số chọn đủ nhỏ res hàm thặng dư toán  4u  4u  4u res   2   f  x, y  (19) x x y y Chúng ta thấy phương trình Biharmonic hai chiều ban đầu định nghĩa    x   y chuyển đổi thành toán chiều  x  y với phương pháp PGD Sơ đồ sai phân trung tâm bậc hai sử dụng để giải phương trình vi phân chiều có dạng phương trình (12) (15) để tìm R  x  S  y  tương ứng Sơ đồ sai phân bậc hai cho đạo hàm bậc hai bậc hàm f  x  tính sau: 2 f  x  f  x  h   f  x   f  x  h    O  h  (20) x h2 f  x  2h   f  x  h   f  x   f  x  h   f  x  2h  4 f  x   x h kích thước bước lưới h2  O  h  (21) Kết mơ số Bài tốn 1: Xét hình chữ nhật 0  x  1,  y  1 có phương trình vi phân chủ đạo phương trình (2), với vế phải cho sau: 67 Lê Quốc Cường Phương pháp proper generalized decomposition cho toán q (22) sin  x  sin  y  D Trong trường hợp với điều kiện biên gối tựa đơn bốn cạnh sau 2w 2w w  0,  x  x  (23), w  0,  y  y  (24) x y q lời giải xác toán w exact  sin  x  sin  y  (25) 4 D Phương pháp PGD áp dụng cho tốn với thơng số mô sau: p  1 , E  1092000 MPa , h  0.01 ,   0.3 Miền tính tốn rời rạc với lưới 100 cho chiều trục x y Hình trình bày lời giải PGD cho độ võng mỏng với điều kiện biên gối tựa đơn bốn cạnh Sai số lời giải PGD lời giải xác trình bày hình f  x, y   Hình Lời giải toán phương pháp PGD cho mỏng với điều kiện biên gối tựa đơn bốn cạnh lưới 100  100 Hình Sai số wexact  w lời giải giải tích lời giải PGD lưới 100  100 Bài toán 2: Trong trường hợp bị ngàm bốn cạnh, ta có điều kiện biên tương ứng sau w (26) w  0,  x  x  , x w (27) w  0,  y  y  y Các thông số mô cho tương tự trường hợp toán 1, vế phải phương trình vi phân chủ đạo cho sau: q 56400  a  10ax  15x   b  y  y D q  18800x  6a  20ax  15 x  y  6b2  20by  15 y  D q  56400  a  x  x  b2  10by  15 y  D f  x, y   (28) Với điều kiện biên bị ngàm bốn cạnh lời giải xác tốn 68 Tạp chí Khoa học Đại học Thủ Dầu Một Số 1(32)-2017 q 2 (29) 2350 x  x  a  y  y  b  D Hình trình bày lời giải PGD cho độ võng với điều kiện biên ngàm chặt bốn cạnh Sai số lời giải PGD lời giải xác trình bày hình wexact  Hình Lời giải tốn phương pháp PGD cho mỏng với điều kiện biên ngàm bốn cạnh Hình Sai số wexact  w lời giải giải tích lời giải PGD lưới 100  100 Kết luận Trong báo này, sử dụng phương pháp PGD để phân tích tốn mỏng chịu uốn với điều kiện biên khác Các kết mô cho thấy đồng thuận tốt phương pháp đề xuất với lời giải xác toán Với việc đưa toán đa chiều toán chiều, phương pháp PGD giúp làm giảm chi phí tính tốn tiết kiệm nhớ máy tính TÀI LIỆU THAM KHẢO [1] F Chinesta, A Ammar, E Cueto (2010), Proper generalized decomposition of multiscale models, Int J Numer Methods Eng, 83(8–9), pp 1114–1132 [2] A Dumon, C Allery, A Ammar (2011), Proper general decomposition (PGD) for the resolution of Navier–Stokes equations, Journal of Computational Physics, 230, pp 1387–1407 [3] E Prulière, F Chinesta, A Ammar, A Leygue, A Poitou (2013), On the solution of the heat equation in very thin tapes, International Journal of Thermal Sciences, 65, pp 148–157 [4] P Vidal, L Gallimard, O Polit (2012), Composite beam finite element based on the Proper Generalized Decomposition, Computers and Structures, 102–103, pp 76–86 [5] A.J Chorin (1968), Numerical solution of the Navier-Stokes equations, Math Comput, 22, pp 745–762 [6] U Ghia, K Ghia, C Shin (1982), High-re solutions for incompressible flow using the Navier– Stokes equations and a multigrid method, Journal of Computational Physics, 48, pp 387–411 [7] S P Timoshenko, S Woinowsky_Krieger (1970), Theory of Plates and Shells, McGraw-Hill, New York 69 ... toán mỏng chịu uốn Phương pháp PGD áp dụng để đưa phương trình vi phân đạo hàm riêng toán mỏng chịu uốn không gian hai chiều thành chuỗi phương trình vi phân khơng gian chiều Sau đó, phương pháp. .. tham số mà toán cần khảo sát Bài toán phân tích mỏng chịu uốn thực thành công với nhiều phương pháp số khác (phần tử hữu hạn, sai phân hữu hạn, phương pháp phổ, …) Trong báo này, phương pháp PGD... dụng phương pháp PGD để phân tích tốn mỏng chịu uốn với điều kiện biên khác Các kết mô cho thấy đồng thuận tốt phương pháp đề xuất với lời giải xác toán Với việc đưa toán đa chiều toán chiều, phương