A Fortran implementation of isogeometric analysis for thin plate problems with the penalty method. Issue 7 (3rd October 2016)

• Record Type: Journal Article
• Authors: Chang, Feng; Wang, Weiqiang; Liu, Yan; Qu, Yanpeng
• Is Part Of: Engineering computations. Volume 33:Issue 7(2016)
• Page Start: 2149
• Page End: 2164
• Publication Date: 2016-10-03
• Subjects: Fortran -- NURBS -- Isogeometric analysis -- Kirchhoff-Love plates -- The penalty method; Computer-aided engineering -- Periodicals; Computer graphics -- Periodicals; 620.00285
• DOI: 10.1108/EC-10-2015-0306
• Languages: English
• ISSNs: 0264-4401

Abstract

Abstract : Purpose: As one of the earliest high-level programming languages, Fortran with easy accessibility and computational efficiency is widely used in the engineering field. The purpose of this paper is to present a Fortran implementation of isogeometric analysis (IGA) for thin plate problems. Design/methodology/approach: IGA based on non-uniform rational B-splines (NURBS) offers exact geometries and is more accurate than finite element analysis (FEA). Unlike the basis functions in FEA, NURBS basis functions are non-interpolated. Hence, the penalty method is used to enforce boundary conditions. Findings: Several thin plate examples based on the Kirchhoff-Love theory were illustrated to demonstrate the accuracy of the implementation in contrast with analytical solutions, and the efficiency was validated in comparison with another open method. Originality/value: A Fortran implementation of NURBS-based IGA was developed to solve Kirchhoff-Love plate problems. It easily obtained high-continuity basis functions, which are necessary for Kirchhoff formulation. In comparison with theoretical solutions, the numerical examples demonstrated higher accuracy and faster convergence of the Fortran implementation. The Fortran implementation can well solve the time-consuming problem, and it was validated by the time-consumption comparison with the Matlab implementation. Due to the non-interpolation of NURBS, the penalty method was used to impose boundary conditions. A suggestion of the selection of penalty coefficients was given.