Bảng 3.7: Vị trí đặt giá đỡ, độ võng lớn nhất và tọa độ tương ứng của bản với 2 LPISe1/π e2/π umax (xmax1−xmax2, ymax) e1/π e2/π umax (xmax1−xmax2, ymax)
0.05 0.1 1.0136 (0.4938 - 0.5, 0) 0.1 0.2 0.5222 (0.4906 - 0.5, 0) 0.2 0.3 0.2319 (0.4844 - 0.5, 0.4063) 0.3 0.4 0.1644 (0.4375 - 0.5, 0.5) 0.1 0.3 0.2341 (0.4844 - 0.5, 0.3906) 0.2 0.4 0.1639 (0.4656 - 0.5, 0.5) 0.1 0.4 0.1638 (0.4672 - 0.5, 0.5) Hình 3.25: Đường đồng mức của độ võng trong trường hợp e1/π = 0.1, e2/π = 0.2 Hình 3.26: Đường đồng mức của độ võng trong trường hợp e1/π = 0.2, e2/π= 0.4
chuyển từ cạnh gối đơn đến tâm của bản và khi độ dài của giá đỡ tăng. (c) Trường hợp 3: Khi tải trọng phân bố không đều và là một hàm chẵn theo biến y, cụ thể q = cos(y). Tất cả các tham số lặp và các dữ liệu khác giống như Trường hợp 2 và kết quả nhận được độ chính xác 10−4
sau 16 lần lặp.
Qua các kết quả nghiên cứu về độ võng của bản có giá đỡ bên trong với các điều kiện biên ngàm hoặc gối tự do giúp cho chúng ta có thể xác định được độ võng của bản phụ thuộc vào các thông số như: số giá đỡ bên trong, vị trí đặt các giá đỡ, độ dài của giá đỡ, tải trọng tác động lên bản,...Bằng phương pháp lặp tìm nghiệm của bài toán song điều hòa với điều kiện biên hỗn hợp được đề xuất ở trên, chúng ta có thể xác định được độ võng của
bản trong từng trường hợp có thể xảy ra. Điều này rất có ý nghĩa trong thực tế.
Kết luận chương 3.
Chương 3 đã trình bày các kết quả nghiên cứu mới trên cơ sở kết hợp hai ý tưởng hạ cấp phương trình và chia miền, luận án đã đề xuất một phương pháp lặp giải bài toán song điều hòa với các điều kiện biên phức tạp, nghiên cứu sự hội tụ của phương pháp và kiểm tra sự hội tụ bằng nhiều ví dụ thử nghiệm. Áp dụng giải hai mô hình bài toán trong cơ học, đó là bài toán về vết nứt và bài toán về độ võng của bản có một hoặc hai giá đỡ bên trong. Qua các kết quả nghiên cứu cả về lý thuyết và thử nghiệm, so sánh với một số phương pháp như: Phương pháp SFBIM giải bài toán vết nứt, phương pháp phương trình chuỗi cặp giải bài toán bản với một giá đỡ, các kết quả tính toán đã chứng tỏ hiệu quả của phương pháp là: Các phép lặp hội tụ nhanh, sự tính toán là đơn giản và dễ dàng thực hiện trên máy tính điện tử, khắc phục được những khó khăn của các phương pháp khác khi trên biên xuất hiện nhiều điểm phân cách các dạng điều kiện biên. Phương pháp lặp kết hợp này còn có thể áp dụng cho các bài toán:
• Với một hoặc nhiều giá đỡ bên trong, các giá đỡ này có thể được đặt ở vị trí không đối xứng qua tâm của bản.
• Khi tải trọng được phân bố không đều và là một hàm chẵn theo biến
y.
• Khi có một cạnh được gối ngang với hạn chế quay đàn hồi không đều:
D∆u+ K(s)∂u ∂n = 0.
Phương pháp lặp cũng có thể sử dụng để giải bài toán về thanh trượt phẳng (The planar stick-slip) [26] với nhiều phần thanh trượt, trong khi với phương pháp SFBIM thì khó có thể áp dụng được.
KẾT LUẬN CHUNG
Kết hợp các ý tưởng chia miền, hạ cấp phương trình và lặp hiệu chỉnh đạo hàm, luận án đã đề xuất các phương pháp hiệu quả giải một số bài toán đối với phương trình elliptic với các điều kiện biên hỗn hợp phức tạp nảy sinh từ cơ học và vật lý. Các kết quả chính của luận án bao gồm: (1) Phát triển phương pháp chia miền kết hợp với kỹ thuật lặp hiệu chỉnh
giá trị đạo hàm trên biên phân chia giải bài toán biên elliptic với hệ số gián đoạn, đưa bài toán mặt phân cách trong môi trường không đồng nhất về các bài toán con, trong các miền ở đó tính chất của môi trường là liên tục. Các kết quả đã được chứng minh bằng lý thuyết và kiểm tra bằng các chương trình thực nghiệm tính toán trên một số miền cả đơn giản và phức tạp.
(2) Đề xuất một phương pháp lặp song song giải bài toán biên của phương trình elliptic với điều kiện biên hỗn hợp mạnh, chứng minh sự hội tụ của phương pháp, tìm khoảng tham số lặp tối ưu cho một trường hợp riêng và áp dụng giải bài toán Motz. Phương pháp này cho phép giải bài toán biên hỗn hợp trên các hệ thống tính toán song song.
song điều hòa với điều kiện biên hỗn hợp mạnh. Các kết quả đã được chứng minh chặt chẽ về mặt lý thuyết và được kiểm tra trên nhiều thực nghiệm tính toán. Đặc biệt, hiệu quả của phương pháp đã được khẳng định qua các kết quả giải các bài toán:
- Bài toán vết nứt
- Bài toán về độ uốn của bản mỏng có một hoặc hai giá đỡ bên trong. Đặc biệt, kết quả về bài toán bản với hai giá đỡ bên trong là hoàn toàn mới, và đang trong quá trình in xuất bản trên tạp chí Journal of Engineering Mathematics (SCI).
Các thử nghiệm tính toán đã được thực hiện trên nhiều ví dụ và so sánh với các phương pháp khác (Phương pháp SFBIM, phương pháp phương trình chuỗi cặp,...) chứng tỏ tính hiệu quả của các phương pháp đã được đề xuất.
(4) Xây dựng một thư viện chương trình giải số bài toán biên hỗn hợp yếu của phương trình elliptic với hệ số hằng trong miền chữ nhật với các điều kiện biên khác nhau dựa trên thuật toán thu gọn khối lượng tính toán của Samarskij-Nikolaev. Thư viện chương trình này là công cụ quan trọng để giải số các bài toán phức tạp được nghiên cứu trong luận án.
Luận án mở ra một số vấn đề có thể tiếp tục nghiên cứu:
• Nghiên cứu biện pháp làm tăng tốc độ hội tụ của các phương pháp lặp trên cơ sở phương pháp ngoại suy tham số .
• Phát triển các phương pháp lặp đã đề xuất đối với các miền hình học phức tạp và điều kiện biên phức tạp hơn.
trong cơ học và vật lý.
• Nghiên cứu phương pháp lặp dựa trên chia miền đối với các dạng phương trình hyperbolic và parabolic.
Danh mục các công trình đã công bố
(1) Dang Q. A, Trương H. H, Vu V. Q,Domain Decomposition Method for Elliptic Interface Problems, Vietnam Journal of Math. Vol. 39, No. 4, 485-499, 2011.
(2) Dang Q. A, Truong H. H, Vu V. Q, Iterative method for a biharmonic problem with crack, Applied Mathematical Sciences. Vol. 6, No. 62, 3095 - 3108, 2012.
(3) Dang Q. A, Trương H. H, Simple Iterative Method for Solving Prob- lems for Plates with Partial Internal Supports, Journal of Engineering Mathematics (SCI), DOI:10.1007/s10665-013-9652-7 (in press).
(4) Đặng Quang Á, Trương Hà Hải, Phương pháp lặp song song giải bài toán biên hỗn hợp mạnh đối với phương trình elliptic, Kỷ yếu hội thảo quốc gia "Một số vấn đề chọn lọc của CNTT và TT", NXB KH và KT, 370-382, 2010.
(5) Đặng Quang Á, Vũ Vinh Quang, Trương Hà Hải, Tiếp cận chia miền tới bài toán về mặt phân cách, Kỷ yếu hội thảo quốc gia "Một số vấn đề chọn lọc của CNTT và TT", NXB KH và KT, 311-320, 2011. (6) Trương Hà Hải, Vũ Vinh Quang, Nguyễn Thị Tuyển, Xây dựng bộ
chương trình RC2009 giải số bài toán biên elliptic với hệ số hằng số. Tạp chí Khoa học và Công nghệ, Đại học Thái nguyên, 7(69), 63-70, 2010.
(7) Vũ Vinh Quang, Trương Hà Hải. Kết quả cải tiến sơ đồ chia miền đối với bài toán song điều hòa. Tạp chí Khoa học và Công nghệ, Đại học Thái nguyên, 1(75), 38-47, 2008.
(8) Vũ Vinh Quang, Trương Hà Hải, Cao Thị Anh Thư, Mô hình tính toán song song giải bài toán biên hỗn hợp mạnh dựa trên chia miền. Tạp chí Khoa học và Công nghệ, Đại học Thái nguyên, 2 (50), 52-57, 2009.
Tài liệu tham khảo
[1] Aubin J. P.Approximation of elliptic boundary-value problems, Wiley- Interscience, 1971.
[2] Adams, R.A. Sobolev spaces. Academic Press, New York. 1975.
[3] Abramov A. A. and Ulijanova V. I., A method for solving biharmonic- type equations with a singularly occurring small parameter, J. of Comp. Math. and Math. Phys., v. 32, No 4, 481-487, 1992.
[4] Andrea Toselli, Olof Widlund, Domain Decomposition Methods: Algo- rithms and Theory, Springer-Verlag Berlin Heidelberg 2005.
[5] Bosko S. Jovanovie and Lubin G. Vulkov, Finite Difference Approx- imation of an Elliptic Interface Problem with Variable Coefficients, Numerical Analysis and Its Applications, Lecture Notes in Computer Science, Vol 3401/2005, 46-55, 2005.
[6] Bialecki B., Karageorghis A., A Legendre Special Galerkin Method for the Biharmonic Dirichlet Problem, SIAM Journal on Scientific Computing, Vol 22, Issue 5, 1549-1569, 2000.
[7] Cioranescu D. and Patrizia D., An Introduction to Homogenization, Oxford Press, 1999.
[8] Cakoni F., Hsiao G. C. and Wendland W. L., On the Boundary In- tegral Equation Method for a Mixed Boundary Value Problem of the Biharmonic Equation, Complex Variables, Theory and Application, Vol 50, 681-696, 2005.
[9] Ciarlet P. G., The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, 1978.
[10] Dang Q. A, Approximate method for solving an elliptic problem with discontinuous coefficients, Journal of Comput. and Applied Math.,51, 193-203, 1994.
[11] Dang Q. A, Parametric extrapolation as a parallel method in math- ematical physics., Tạp chí Tin học và Điều khiển học , No 1, 1-9, 2001.
[12] Dang Q. A, V. V. Quang, Domain decomposition method for solv- ing an elliptic boundary value problem, in book: Method of Complex and Clifford Analysis, SAS International Publications, Delhi, 309-319, 2006.
[13] Dang Q. A, Vu V. Q., A domain decomposition method for strongly mixed boundary value problems for the Poisson equation, In book H.G. Bock et al (eds): Modeling, Simulation and Optimization of Complex Processes, Springer, 65-76, 2012.
[14] Dang Q. A, Boundary operator method for approximate solution of biharmonic type equation, Vietnam Journal of Math. Vol. 22, No. 1-2, 114-120, 1994.
[15] Dang Q. A, Mixed boundary-domain operator method in approximate solution of biharmonic type equation, Vietnam Journal of Math. Vol. 26, No. 3, 243-252, 1998.
[16] Dang Q. A, Iterative method for solving the Neumann boundary value problem for biharmonic type equation, Journal of Comput. and Ap- plied Math., Vol. 196, No 2, 634-643, 2006.
[17] Dang Q. A and Le T. S, Iterative methods for Solving a Problem with Mixed Boundary Conditions for Biharmonic Equations, Advances in Applied Mathematics and Mechanics, Vol. 1, No. 5, 1-16, 2009.
[18] Dang Q. A, Trương H. H, Vu V. Q, Domain Decomposition Method for Elliptic Interface Problems, Vietnam Journal of Math. Vol. 39, No. 4, 485-499, 2011.
[19] Dang Q. A, Truong H. H, Vu V. Q, Iterative method for a biharmonic problem with crack, Applied Mathematical Sciences. Vol. 6, No. 62, 3095 - 3108, 2012.
[20] Dang Q. A, Trương H. H, Simple Iterative Method for Solving Prob- lems for Plates with Partial Internal Supports, Journal of Engineering Mathematics (SCI), DOI:10.1007/s10665-013-9652-7 (in press).
[21] Đặng Quang Á, Trương Hà Hải, Phương pháp lặp song song giải bài toán biên hỗn hợp mạnh đối với phương trình elliptic, Kỷ yếu hội thảo quốc gia "Một số vấn đề chọn lọc của CNTT và TT", NXB KH và KT, 370-382, 2010.
[22] Đặng Quang Á, Vũ Vinh Quang, Trương Hà Hải, Tiếp cận chia miền tới bài toán về mặt phân cách, Kỷ yếu hội thảo quốc gia "Một số vấn đề chọn lọc của CNTT và TT", NXB KH và KT, 311-320, 2011. [23] Ewing. R. E. , Z. Li, T. Lin, and Y. Lin, The immersed finite volume
element methods for the elliptic interface problems, Mathematics and Computers in Simulation 50, no. 1–4, 63–76, 1999.
[24] Elliotis M., Georgiou G. and Xenophontos C., Solution of the stick- slip problem with the singular function boundary integral method, 5th GRACM International Congress on Computational Mechanics, Limas- sol, 29 June, 2005.
[25] Elliotis M., Georgiou G. and Xenophontos C., The singular function boundary intergral method for biharmonic problems with crack singu- larities, Eng. Anal. Bound. Elem. 31, 209-215, 2007.
[26] Elliotis M., Georgiou G. and Xenophontos C.,Solution of planar New- tonian stick-slip problem with the singular function boundary integral method, Int. J. Numer. Mech. Fluids, 48: 1001-1021, 2005.
[27] Grisvard P., Elliptic problems in nonsmooth domains, Pitman Pub- lishing, Boston 1985.
[28] Glowinski R. and Pironneau O., Numerical methods for the first bi- harmonic equation and for the two-dimensional Stokes problem, SIAM, Society for Industrial and Applied Mathematics, Vol. 21, No. 2, 167- 212, 1979.
[30] Guminiak M, Sygulski R., The analysis of internally supported thin plates by the boundary element method. Part 1 Static analysis. Foun- dation of Civil and Environmental Engineering, 9: 17-41, 2007.
[31] Guminiak M, Sygulski R., The analysis of internally supported thin plates by the boundary element method. Part 2 Free vibration analysis. Foundation of Civil and Environmental Engineering, 9: 43-74, 2007. [32] Herrera I., Unified Theory of Domain Decomposition Methods, 14th
Int. Conf. on Domain Decomposition Methods, Editors: Ismael Her- rera, David E. Keyes, Olof B. Widlund, Robert Yates, www.DDM.org, 243-248 2003.
[33] He. X. M. , T. Linand Y. Lin, A Bilinear Immersed Finite Volume Element Method for the Diffusion Equation with Discontinuous Coef- ficient, Commun. Comput. Phys., 6, 185-202, 2009.
[34] Hou L.S. and Jangwoon Lee, A Robin-Robin non-overlapping domain decomposition method for an elliptic boundary control problem, In- ternational Journal of Numerical Analysis and Modeling, Volume 8, Number 3, 443–465, 2011.
[35] Jeon Y.,Scalar boundary integral equation formulas for the biharmonic equation-numerical experiments, Journal of Comp. and App. Math. (JCAM) 115, 269-282, 2000.
[36] Jeon Y., New indirect scalar boundary intergral equation formulas for the biharmonic equation, JCAM 135, 313-324, 2001.
[37] Katsikadelis J. T. and Kallivokas L. F.,Clamped Plates on Pasternak- Type Elastic Foundation by The Boundary Element Method , J. Appl. Mech. ASME, 53, 909-917, 1986.
[38] Katsikadelis J. T. and Kallivokas L. F.,Plates on Biparametric Elastic Foundation by BDIE Method, J. Engrg. Mech., 114, 847-875, 1998. [39] Li. Z, An overview of the immersed interface method and its applica-
tions, Taiwanese Journal of Mathematics, v. 7, no 1, 1-49, 2003. [40] Li. Z. and K. Ito, The Immersed Interface Method: Numerical So-
lutions of PDEs Involving Interfaces and Irregular Domains, SIAM, Philadelphia, 2006.
[41] Liang Chern, Yu-Chen Shu, A coupling interface method for elliptic interface problems, Journal of Computational Physics, 225, 2138–2174, 2007.
[42] Lions J. L. and E. Magenes,Problemes aux Limites Non Homogens et Application, Vol. 1, Dunod, Paris, 1968.
[43] Lions P. L., On the Schwarz alternating methods III: a variant for nonoverlapping subdomains, in Third International Symposium on Do- main Decomposition Methods for Partial Differential Equations, T. F. Chan, R. Glowinski, J. Periaux, and O. B. Wildlund, eds., SIAM, Philadelphia, PA, 202-231, 1990.
[44] Li Z. C. and T. T. Lu, Singularities and Treatments of Elliptic Bound- ary Value Problems, Mathematical and Computer Modeling, 31, 97-45, 2000.
[45] Li Z. C, Lu T. T, Hu H. Y, The collocation Trefftz method for bihar- monic equations with crack singularities, Eng. Anal. Bound. Elem.; 28, 79-96, 2004.
[46] Lurie A. I., Theory of Elasticity, published by Springer, 2005.
[47] Le Tallec P. et al., Domain decomposition methods for large linearly elliptic three dimensional problems, Rapprts de Recherche, No 1182, INRIA, France, 1990.
[48] Li. Z. C., Y. L. Chan, G. C. Georgiov, C. Xenophontos, Special Boundary Approcimation Methods For Laplace Equation problems with Boundary Singularities Application to Motz problem, J. Com- puters and Mathematics with Applications 51, 115-142, 2006.
[49] Lube G., M¨uller L. and M¨uller H., A new non-overlapping domain decomposition method for stabilized finite element methods applied to the non-stationary Navier-Stokes equations. Numer. Linear Algebra Appl., 7: 449–472, 2000.
[50] Meller N. A. and Dorodnisyn A. A., Some approaches to solving sta- tionary Navier–Stokes equations, Zh. Vychisl. Mat. Mat. Fiz., 8:2, 393–402, 1968.
[51] Motz. H, The treatment of singularities in relaxation methods, Quart. Appl. Math. 4, 371-377, 1946.
[52] Mandal B.N. and N. Mandal, Advances in Dual Integral Equations. Research Notes in Mathematics, Chapman and Hall-CRC, New York, 1999.
[53] Mandrik P.A.,The method of the dual integral equation for analysis of heat transfer processes, Mathematical Modeling and Analysis, 6:280- 288, 2001 .
[54] Meleshko V. V., Selected topics in the history of the two-dimensional biharmonic problem, Applied Mechanics Reviews, vol. 56, issue 1, 33- 85, 2003.
[55] Pawlak Z., Guminiak M., The application of fundamental solutions in static analysis of thin plates resting on the internal elastic support, Foundation of Civil and Environmental Engineering, No 11: 67-96, 2008.
[56] Pollikkas A., Karageorghis A., Georgiou G., Methods of fundamen- tal solutions for harmonic and biharmonic boundary value problems, Computational Mechanics 21, 416-423, 1998.
[57] Popov G.I, Rostovtsev N. A., Contact and Mixed Problems of Elastic- ity, Proceedings of Symposium on Theoretical and Applied Mechanics, Moscow, No 3, 235-252, 1966.
[58] Palsev B. V., On the expansion of the Dirichlet problem and a mixed