Luận án Development of novel meshless method for limit and shakedown analysis of structures and materials

Nội dung tài liệu: Luận án Development of novel meshless method for limit and shakedown analysis of structures and materials

1. MINISTRY OF EDUCATION AND TRAINING HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY AND EDUCATION HO LE HUY PHUC DEVELOPMENT OF NOVEL MESHLESS METHOD FOR LIMIT AND SHAKEDOWN ANALYSIS OF STRUCTURES & MATERIALS DOCTORAL THESIS MAJOR: ENGINEERING MECHANICS Ho Chi Minh city, 3rd May 2020
2. MINISTRY OF EDUCATION AND TRAINING HO CHI MINH CITY UNIVERSITY OF TECHNOLOGY AND EDUCATION HO LE HUY PHUC DEVELOPMENT OF NOVEL MESHLESS METHOD FOR LIMIT AND SHAKEDOWN ANALYSIS OF STRUCTURES & MATERIALS MAJOR: ENGINEERING MECHANICS Supervisors: 1. Assoc. Prof Le Van Canh 2. Assoc. Prof Phan Duc Hung Reviewer 1: Reviewer 2: Reviewer 3:
3. Declaration of Authorship I declare that this is my own research. The data and results stated in the thesis are honest and have not been published by anyone in any other works. Ho Chi Minh city, 3rd August 2020 PhD candidate HO LE HUY PHUC i
4. Acknowledgements The research presented in this thesis has been carried out in the framework of a doctorate at Faculty of Civil Engineering, Ho Chi Minh city University of Technology and Education, Vietnam. This work would have never been possible without the support and help of many people to whom I feel deeply grateful. First and foremost, I would like to express my most sincere thanks to my su- pervisors, Assoc. Prof. Le Van Canh and Assoc. Prof. Phan Duc Hung, for their guidance, valuable academic advice, mental support and constant encouragement during the course of this work. I am deeply indebted to my major supervisor, As- soc. Prof. Le Van Canh. He is one of most influential people in my life, both profes- sionally and personally. His guidance is precious, helping me develop the personal skills needed to succeed in future work. I would like to thank the co-author of my papers - Prof. Tran Cong Thanh for his encouragement, support and guidance. I would also like to express my admiration for his unsurpassed knowledge of mathematics and numerical methods. I really appreciate the financial support received from the Institute for Com- putational Science and Technology (ICST) - HCMC, the Science and Technology Incubator Youth Program - HCMC, and International University - VNU-HCMC throughout the research projects. I take this opportunity to thank my colleagues in International University - VNU-HCMC, HCMC University of Technology and Education, and HUTECH Uni- versity, especially Dr. Tran Trung Dung, PhD candidate Nguyen Hoang Phuong, PhD candidate Do Van Hien, Dr. Khong Trong Toan and Dr. Vo Minh Thien, for fruitful discussions about a range of topics and their mental support. I sincerely thank my parents and my younger sisters for their unconditional love and support. I am also definitely indebted to my wife, Nguyen My Lam, for her love, understanding and encouraging me whenever I needed motivation. ii
5. Acknowledgements Finally, I would like to dedicate this thesis to my little son - Ho Nguyen Nhat Duy. No word can describe my love for him. Ho Chi Minh city, 3rd August 2020 PhD candidate HO LE HUY PHUC iii
6. Abstract The proposed research is essentially concerning on the development of powerful numerical methods to deal with practical engineering problems. The direct methods requiring the use of a strong mathematical tool and a proper numerical discretiza- tion are considered. The current work primarily focuses on the study of limit and shakedown analysis allowing the rapid access to the requested information of structural design with- out the knowledge of whole loading history. For the mathematical treatment, the problems are formulated in form of minimizing a sum of Euclidean norms which are then cast as suitable conic programming depending on the yield criterion, e.g. second order cone programming (SOCP). In addition, a robust numerical tool also requires an excellent discretization strat- egy which is capable of providing stable and accurate solutions. In this study, the so-called integrated radial basis functions-based mesh-free method (iRBF) is em- ployed to approximate the computational fields. To eliminate numerical instability problems, the stabilized conforming nodal integration (SCNI) scheme is also intro- duced. Consequently, all constrains in resulting problems are directly enforced at scattered nodes using collocation method. That not only keeps size of the optimiza- tion problem small but also ensures the numerical procedure truly mesh-free. One more advantage of iRBF method, which is absent in almost meshless ones, is that the shape function satisfies Kronecker delta property leading the essential boundary conditions to be imposed easily. In summary, the iRBF-based mesh-free method is developed in combination with second order cone programming to provide solutions for direct analysis of structures and materials. The most advantage of proposed approach is that the highly accu- rate solutions can be obtained with low computational efforts. The performance of proposed method is justified via the comparison of obtained results and available ones in the literature. iv
7. Tóm tắt Luận án này hướng đến việc phát triển một phương pháp số mạnh để giải quyết các bài toán kỹ thuật, và phương pháp phân tích trực tiếp được sử dụng. Phương pháp này yêu cầu một thuật toán tối ưu hiệu quả và một công cụ rời rạc thích hợp. Trước tiên, nghiên cứu này tập trung vào lý thuyết phân tích giới hạn và thích nghi, phương pháp được biết đến như một công cụ hữu hiệu để xác định trực tiếp những thông tin cần thiết cho việc thiết kế kết cấu mà không cần phải thông qua toàn bộ quá trình gia tải. Về mặt toán học, các bài toán được phát biểu dưới dạng cực tiểu một chuẩn của tổng bình phương các biến trong không gian Euclide, sau đó được đưa về dạng chương trình hình nón phù hợp với tiêu chuẩn dẻo, ví dụ chương trình hình hón bậc hai (SOCP). Hơn nữa, một công cụ số mạnh còn đòi hỏi phải có kỹ thuật rời rạc tốt để đạt được kết quả tính toán chính xác với tính ổn định cao. Nghiên cứu này sử dụng phương pháp không lưới dựa trên phép tích phân hàm cơ sở hướng tâm (iRBF) để xấp xỉ các trường biến. Kỹ thuật tích phân nút ổn định (SCNI) được đề xuất nhằm loại bỏ sự thiếu ổn định của kết quả số. Nhờ đó, tất cả các ràng buộc trong bài toán được áp đặt trực tiếp tại các nút bằng phương pháp tụ điểm. Điều này không những giúp kích thước bài toán được giữ ở mức tối thiểu mà còn đảm bảo phương pháp là không lưới thực sự. Một ưu điểm nữa mà hầu hết các phương pháp không lưới khác không đáp ứng được, đó là hàm dạng iRBF thỏa mãn đặc trưng Kronecker delta. Nhờ vậy, các điều kiện biên chính có thể được áp đặt dễ dàng mà không cần đến các kỹ thuật đặc biệt. Tóm lại, nghiên cứu này phát triển phương pháp không lưới iRBF kết hợp với thuật toán tối ưu hình nón bậc hai cho bài toán phân tích trực tiếp kết cấu và vật liệu. Thế mạnh lớn nhất của phương pháp đề xuất là kết quả số với độ chính xác cao có thể thu được với chi phí tính toán thấp. Hiệu quả của phương pháp được đánh giá thông qua việc so sánh kết quả số với những phương pháp khác. v
8. Contents Declaration of Authorship i Acknowledgements iii Abstract v Contents ix List of Tables xi List of Figures xvi List of Abbreviations xvii Chapter 1: Introduction 1 1.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1 1.2 Literature review . . . . . . . . . . . . . . . . . . . . . . . . . . . .3 1.2.1 Limit and shakedown analysis . . . . . . . . . . . . . . . . .3 1.2.2 Mathematical algorithms . . . . . . . . . . . . . . . . . . . .4 1.2.3 Discretization techniques . . . . . . . . . . . . . . . . . . . .5 1.2.4 The direct analysis for microstructures . . . . . . . . . . . .7 1.2.5 Mesh-free methods - state of the art . . . . . . . . . . . . . .8 1.3 Research motivation . . . . . . . . . . . . . . . . . . . . . . . . . . 21 1.4 The objectives and scope of thesis . . . . . . . . . . . . . . . . . . . 24 1.5 Original contributions of the thesis . . . . . . . . . . . . . . . . . . 24 1.6 Thesis outline . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25 Chapter 2: Fundamentals 27 2.1 Plasticity relations in direct analysis . . . . . . . . . . . . . . . . . 27 vi
9. Contents 2.1.1 Material models . . . . . . . . . . . . . . . . . . . . . . . . . 27 2.1.2 Variational principles . . . . . . . . . . . . . . . . . . . . . . 31 2.2 Shakedown analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 2.2.1 Upper bound theorem of shakedown analysis . . . . . . . . . 35 2.2.2 The lower bound theorem of shakedown analysis . . . . . . . 36 2.2.3 Separated and unified methods . . . . . . . . . . . . . . . . 38 2.2.4 Load domain . . . . . . . . . . . . . . . . . . . . . . . . . . 38 2.3 Limit analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 2.3.1 Upper bound formulation of limit analysis . . . . . . . . . . 40 2.3.2 Lower bound formulation of limit analysis . . . . . . . . . . 41 2.4 Conic optimization programming . . . . . . . . . . . . . . . . . . . 41 2.5 Homogenization theory . . . . . . . . . . . . . . . . . . . . . . . . . 43 2.6 The iRBF-based mesh-free method . . . . . . . . . . . . . . . . . . 45 2.6.1 iRBF shape function . . . . . . . . . . . . . . . . . . . . . . 46 2.6.2 The integrating constants in iRBF approximation . . . . . . 48 2.6.3 The influence domain and integration technique . . . . . . . 49 Chapter 3: Displacement and equilibrium mesh-free formulation based on integrated radial basis functions for dual yield design 53 3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53 3.2 Kinematic and static iRBF discretizations . . . . . . . . . . . . . . 54 3.2.1 iRBF discretization for kinematic formulation . . . . . . . . 55 3.2.2 iRBF discretization for static formulation . . . . . . . . . . . 57 3.3 Numerical examples . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 3.3.1 Prandtl problem . . . . . . . . . . . . . . . . . . . . . . . . 60 3.3.2 Square plates with cutouts subjected to tension load . . . . 63 3.3.3 Notched tensile specimen . . . . . . . . . . . . . . . . . . . . 65 3.4 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67 vii
10. Contents Chapter 4: Limit state analysis of reinforced concrete slabs using an integrated radial basis function based mesh-free method 68 4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68 4.2 Kinematic formulation using the iRBF method for reinforced con- crete slab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69 4.3 Numerical examples . . . . . . . . . . . . . . . . . . . . . . . . . . . 73 4.3.1 Rectangular slabs . . . . . . . . . . . . . . . . . . . . . . . 73 4.3.2 Regular polygonal slabs . . . . . . . . . . . . . . . . . . . . 77 4.3.3 Arbitrary geometric slab with a rectangular hole . . . . . . 79 4.4 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81 Chapter 5: A stabilized iRBF mesh-free method for quasi-lower bound shakedown analysis of structures 82 5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82 5.2 iRBF discretization for static shakedown formulation . . . . . . . . 83 5.3 Numerical examples . . . . . . . . . . . . . . . . . . . . . . . . . . . 88 5.3.1 Punch problem under proportional load . . . . . . . . . . . 88 5.3.2 Thin plate with a central hole subjected to variable tension loads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 5.3.3 Grooved plate subjected to tension and in-plane bending loads 95 5.3.4 A symmetric continuous beam . . . . . . . . . . . . . . . . . 98 5.3.5 A simple frame with different boundary conditions . . . . . 101 5.4 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104 Chapter 6: Kinematic yield design computational homogenization of micro-structures using the stabilized iRBF mesh-free method 106 6.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106 6.2 Limit analysis based on homogenization theory . . . . . . . . . . . 107 6.3 Discrete formulation using iRBF method . . . . . . . . . . . . . . . 109 6.4 Numerical examples . . . . . . . . . . . . . . . . . . . . . . . . . . . 110 viii