Phát triển phương pháp phần tử hữu hạn đẳng hình học để phân tích và điều khiển đáp ứng kết cấu tấm nhiều lớp

THE WORK IS COMPLETED AT HCM CITY UNIVERSITY OF TECHNOLOGY AND EDUCATION Supervisor 1: Assoc Prof Dr NGUYEN XUAN HUNG Supervisor 2: Assoc Prof Dr DANG THIEN NGON PhD thesis is protected in front of EXAMINATION COMMITTEE FOR PROTECTION OF DOCTORAL THESIS HCM CITY UNIVERSITY OF TECHNOLOGY AND EDUCATION, Date month year ORIGINALITY STATEMENT I, Nguyen Thi Bich Lieu, hereby assure that this dissertation is my own work, done under the guidance of Assoc Prof Dr Nguyen Xuan Hung and Assoc Prof Dr Dang Thien Ngon with the best of my knowledge The data and results stated in the dissertation are honest and were not been published by any works Ho Chi Minh City, October 2019 Nguyen Thi Bich Lieu i ACKNOWLEDGEMENTS This dissertation has been carried out in the Faculty of Civil Engineering, HCM City University of Technology and Education, Viet Nam The process of conducting this thesis brings excitement but has quite a few challenges and difficulties And I can say without hesitation that it has been finished thanks to the encouragement, support and help of my professors and colleagues First of all, I would like to express my deepest gratitude to Assoc Prof Dr Nguyen Xuan Hung and Assoc Prof Dr Dang Thien Ngon, especially Assoc Prof Dr Nguyen Xuan Hung from CIRTech Institute, Ho Chi Minh City University of Technology (HUTECH), Vietnam for having accepted me as their PhD student and for the enthusiastic guidance and mobilization during my research Also, I would like to sincerely thank Dr Thai Hoang Chien, a close brother, for his helpful guidance at first step of doing research and his support for my overcoming of the hardest time Secondly, I would like also to acknowledge Msc Nguyen Van Nam, Faculty of Mechanical Technology, Industrial University of Ho Chi Minh City, Vietnam for their troubleshooting and the cooperation in my study Furthermore, I am grateful to Chau Nguyen Khanh and the staffs at CIRTech Institute, HUTECH, Vietnam for their professional knowledge, interactive discussion, and immediate support Thirdly, I take this chance to thank all my nice colleagues at the Faculty of Civil Engineering, Ho Chi Minh City University of Technology and Education, for their professional advice and friendly support Finally, this dissertation is dedicated to my family, especially my beloved husband, who has always given me valuable encouragement and assistance Nguyen Thi Bich Lieu ii ABSTRACT Isogeometric analysis (IGA) was introduced in 2005 by Hughes et al [5] as a breakthrough in numerical simulation The main advantage of the IGA is to use the same basis function to describe the geometry and to approximate the problem unknowns It integrates Computer Aided Design (CAD) and Computer Aided Engineering (CAE) and so far the effectively numerical tool for the analysis of a variety of practical problems The computational cost is decreased significantly as the meshes are generated within the CAD IGA produces the results with higher accuracy because of the smoothness and the higher-order continuity between elements For the last decade of development, isogeometric analysis has surpassed the standard finite elements in terms of effectiveness and reliability for various problems, especially for the ones with complex geometry Owing to its important role in many engineering structures and modern industries, laminated plate structures are widely used in a diverse array of structures in many areas such as aviation, shipbuilding and civil engineering Laminated plates have excellent mechanical properties, including high strength to weight and stiffness to weight ratios, wear resistance, light weight and so on Besides possessing the superior material properties, the laminated composites also supply the advantageous design through the arrangement of the stacking sequence and layer thickness to obtain the desired characteristics, that’s why they have received considerable attention of many researchers worldwide In this dissertation, an isogeometric finite element formulation is developed based on Bézier extraction to solve various plate problems, using a seven-dof higher- order shear deformation theory for both analysis and control the responses of plate structures One key point in this dissertation is to exploit the distinctive advantage of Bézier extraction in analysis of plate structures In the conventional isogeometric analysis, the B-spline or Non-uniform Rational B-spline (NURBS) basis functions span over the entire domain of structures not just a local domain as Lagrangian shape iii functions in FEM The global structure induces the complex implementation in a traditional finite element context In addition, in order to compute the shape functions, the Gaussian integration points force to transform to parametric space By choosing Bernstein polynomials as the basis functions, IGA will be performed easily similar to the way of implementation in FE framework The B-spline/NURBS basis can be rewritten in form of the combination of Bernstein polynomials and Bézier extraction operator That is called Bézier extraction for B-spline/NURBS/Tspline Although IGA is suitable for the problems which have the higher-order continuity, the findings of using a higher-order shear deformation theory with the C0- continuity show the convieniences for plate analysis Furthermore, both linear and nonlinear responses for four material models including laminated composite plates, piezoelectric laminated composite plates, piezoelectric functionally graded porous plates with graphene platelets reinforcement and functionally graded piezoelectric material porous plates are investigated The control algorithms based on the constant displacement and velocity feedbacks are applied to control linear and geometrically nonlinear static and dynamic responses of the plates, where the effect of the structural damping is considered, based on a closed- loop control with piezoelectric sensors and actuators The predictions of the proposed approach agree well with analytical solutions and several other available approaches Through the analysis, numerical results indicated that the proposed method achieves high reliability as compared with other published solutions Besides, some numerical solutions for PFGPM plates and FG porous reinforced by GPLs may be considered as reference solutions for future work because there have not yet been analytical solutions so far iv TĨM TẮT Phân tích đẳng hình học (IGA) giới thiệu năm 2005 Hughes cộng [5] đột phá tính tốn mơ số Ưu điểm IGA sử dụng hàm sở để mô tả cho hình học xấp xỉ nghiệm số Nó tích hợp việc thiết kế dựa máy tính công nghệ liên quan đến việc sử dụng hệ thống máy tính để phân tích đối tượng hình học CAD (CAE) công cụ số hiệu khác nhằm giải nhiều lớp toán kỹ thuật khác Chi phí tính tốn giảm đáng kể hình học xác tạo CAD, sau đưa vào tính tốn mà khơng bị sai số hình học Hơn nữa, IGA cho kết nghiệm số với độ xác cao tính trơn tính liên tục bậc cao phần tử Trong thập kỷ phát triển gần đây, phân tích đẳng hình học vượt qua phân tích phần tử hữu hạn (FEM) tính hiệu độ tin cậy toán khác nhau, đặc biệt tốn có hình học phức tạp Bởi đóng vai trị quan trọng nhiều kết cấu kỹ thuật công nghiệp đại, kết cấu nhiều lớp sử dụng rộng rãi nhiều lĩnh vực khác chẳng hạn hàng khơng, đóng tàu, kỹ thuật dân dụng, vv Kết cấu nhiều lớp có tính chất học tuyệt vời, bao gồm độ bền độ cứng cao, khả chống mài mịn cao, trọng lượng nhẹ nhiều đặc tính khác ưu việt khác Bên cạnh việc sở hữu đặc tính tốt đó, vật liệu tổng hợp nhiều lớp cịn cung cấp thiết kế thuận lợi thơng qua việc xếp trình tự xếp chồng độ dày lớp để có đặc tính học mong muốn, lý chúng nhận quan tâm nghiên cứu đáng kể nhiều nhà nghiên cứu toàn giới Trong luận án này, cơng thức phần tử hữu hạn đẳng hình học phát triển dựa trích xuất Bézier để giải toán khác nhau, sử dụng lý thuyết biến dạng cắt bậc cao liên tục C0 cho phân tích điều khiển đáp ứng cấu trúc Một điểm luận án khai thác lợi ích vượt trội trích xuất Bézier phân tích kết cấu Trong phân tích đẳng hình học truyền thống thơng thường, hàm sở B-spline hàm NURBS phân bố tồn miền cấu trúc khơng miền cục hàm dạng v Lagrangian FEM Việc hàm dạng phân bố toàn cục làm cho việc thực tính tốn phức tạp Ngồi ra, để tính tốn hàm dạng, điểm tích phân Gauss buộc phải chuyển đổi sang không gian tham số Bằng cách chọn đa thức Bernstein làm hàm sở, IGA thực dễ dàng tương tự cách triển khai phương pháp phần tử hữ hạn Các hàm sở B-spline / NURBS viết lại dạng kết hợp đa thức Bernstein tốn tử trích xuất Bézier Đó gọi trích xuất Bézier cho B-spline / NURBS / T-spline Lý thuyết biến dạng cắt bậc cao với bậc liên tục C0 sử dụng thống cho tất chương Hơn nữa, đáp ứng tuyến tính phi tuyến cho bốn loại vật liệu bao gồm composite nhiều lớp, composite nhiều lớp có dán lớp áp điện, vật liệu chức dán lớp áp điện có lỗ rỗng gia cường graphene vật liệu áp điện chức có lỗ rỗng nghiên cứu Các thuật tốn điều khiển dựa tín hiệu phản hồi chuyển vị vận tốc không đổi áp dụng để điều khiển đáp ứng tĩnh động cho đáp ứng tuyến tính phi tuyến hình học, hiệu ứng giảm chấn cấu trúc xem xét, dựa điều khiển kín với cảm biến truyền động áp điện Thông qua phân tích phần ví dụ số, kết đạt phương pháp đề xuất đạt độ tin cậy cao so với giải pháp khác cơng bố tạp chí uy tín Ngồi ra, số lời giải số cho vật liệu chức dán lớp áp điện có lỗ rỗng gia cường graphene vật liệu áp điện chức có lỗ rỗng coi nguồn tài liệu tham khảo cho nghiên cứu khác tương lai chưa có lời giải giải tích đưa vi CONTENTS ORIGINALITY STATEMENT i ACKNOWLEDGEMENTS ii ABSTRACT iii CONTENTS vii NOMENCLATURE .xii LIST OF TABLES xvi LIST OF FIGURES xx Chapter LITERATURE REVIEW .1 1.1 Introduction 1.2 An overview of isogeometric analysis 1.3 Literature review about materials used in this dissertation 1.3.1 Laminated composite plate .5 1.3.2 Piezoelectric laminated composite plate 1.3.3 Piezoelectric Functionally Graded Porous plates reinforced by Graphene Platelets (PFGP-GPLs) .7 1.3.4 Functionally Graded Piezoelectric Material Porous plates (FGPMP) 1.4 Goal of the dissertation 11 1.5 The novelty of dissertation 12 1.6 Outline .13 1.7 Concluding remarks 15 Chapter .16 ISOGEOMETRIC ANALYSIS FRAMEWORK 16 2.1 Introduction 16 2.2 Advantages of IGA compared to FEM .16 2.3 Some disadvantages of IGA .17 2.4 B-spline geometries 17 vii 2.4.1 B-spline curves 18 2.4.2 B-spline surface .20 2.5 Refinement technique 20 2.5.1 h-refinement 21 2.5.2 p-refinement 23 2.5.3 k-refinement 25 2.6 NURBS basis function 26 2.7 Isogeometric discretization 29 2.8 Numerical integration 30 2.9 Bézier extraction 33 2.9.1 Introduction of Bézier extraction 33 2.9.2 Bézier decomposition and Bézier extraction [97-98] 34 2.10 Concluding remarks 37 Chapter .39 THEORETICAL BASIS 39 3.1 Overview 39 3.2 An overview of plate theories 39 3.2.1 The higher-order shear deformation theory 40 3.2.2 The generalized unconstrained higher-order shear deformation theory (UHSDT) 43 3.2.3 3.3 The C0-type higher-order shear deformation theory (C0-type HSDT) 45 Laminated composite plate 46 3.3.1 Definition of laminated composite plate 46 3.3.2 Constitutive equations of laminated composite plate 47 3.4 Piezoelectric material 50 3.4.1 Introduce to piezoelectric material 50 3.4.2 The basic equation of piezoelectric material 51 3.5 Piezoelectric functionally graded porous plates reinforced by graphene platelets (PFGP-GPLs) 52 viii 3.6 Functionally graded piezoelectric material porous plates (FGPMP) 56 3.7 Concluding remarks 59 Chapter .60 ANALYZE AND CONTROL THE LINEAR RESPONSES OF THE PIEZOELECTRIC LAMINATED COMPOSITE PLATES 60 4.1 Overview 60 4.2 Laminated composite plate formulation based on Bézier extraction for NURBS 60 4.2.1 The weak form for laminated composite plates 60 4.2.2 Approximated formulation based on Bézier extraction for NURBS 62 4.3 Theory and formulation of the piezoelectric laminated composite plates 65 4.3.1 Variational forms of piezoelectric composite plates .65 4.3.2 Approximated formulation of electric potential field 66 4.3.3 Governing equations of motion .67 4.4 Active control analysis 68 4.5 Results and discussions .69 4.5.1 Static analysis of the four-layer [00/900/900/00] square laminated plate 71 4.5.2 Static analysis of laminated circular plate subjected to a uniform distributed load 76 4.5.3 Free vibration of laminated composite square plate 80 4.5.4 Free vibration of laminated circular plate .81 4.5.5 Transient analysis 83 4.5.6 Static analysis of the square piezoelectric laminated composite plate 88 4.5.7 Free vibration analysis of an elliptic piezoelectric composite plate 92 4.5.8 Dynamic control of piezoelectric laminated composite plate 94 4.6 Concluding remarks 95 Chapter 5: 97 ANALYSIS AND CONTROL THE RESPONSES OF PIEZOELECTRIC FUNCTIONALLY GRADED POROUS PLATES REINFORCED BY ix