• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.3
• /
• pp.233-238
• /
• 2008

In this paper, a shape design optimization method for linearly elastic problems is developed using isogeometric approach. In many design optimization problems for practical engineering models, initial raw data usually come from a CAD modeler. Then, designers should convert the CAD data into finite element mesh data since most of conventional design optimization tools are based on finite element analysis. During this conversion, there are some numerical errors due to geometric approximation, which causes accuracy problems in response as well as design sensitivity analyses. As a remedy for this phenomenon, the isogeometric analysis method can be one of the promising approaches for the shape design optimization. The main idea of isogeometric approach is that the basis functions used in analysis is exactly the same as the ones representing the geometry. This geometrically exact model can be used in the shape sensitivity analysis and design optimization as well. Therefore the shape design sensitivity with high accuracy can be obtained, which is very essential for a gradient-based optimization. Through numerical examples, it is verified that the shape design optimization based on an isogeometic approach works well.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.497-504
• /
• 2012

Finite element analysis is to approximate a geometry model developed in computer-aided design(CAD) to a finite element model, thus the conventional shape design sensitivity analysis and optimization using the finite element method have some difficulties in the parameterization of geometry. However, isogeometric analysis is to build a geometry model and directly use the functions describing the geometry in analysis. Therefore, the geometric properties can be embedded in the NURBS basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. In this study, the isogeometric structural analysis and shape design sensitivity analysis in the generalized curvilinear coordinate(GCC) systems are discussed for the curved geometry. Representing the higher order geometric information, such as normal, tangent and curvature, yields the isogeometric approach to be the best way for generating exact GCC systems from a given CAD geometry. The developed GCC isogeometric structural analysis and shape design sensitivity analysis are verified to show better accuracy and faster convergency by comparing with the results obtained from the conventional isogeometric method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.345-352
• /
• 2014

In this paper, the CAD data for the optimal shape design obtained by isogeometric shape optimization is directly used to fabricate the specimen by using 3D printer for the experimental validation. In a conventional finite element method, the geometric approximation inherent in the mesh leads to the accuracy issue in response analysis and design sensitivity analysis. Furthermore, in the finite element based shape optimization, subsequent communication with CAD description is required in the design optimization process, which results in the loss of optimal design information during the communication. Isogeometric analysis method employs the same NURBS basis functions and control points used in CAD systems, which enables to use exact geometrical properties like normal vector and curvature information in the response analysis and design sensitivity analysis procedure. Also, it vastly simplify the design modification of complex geometries without communicating with the CAD description of geometry during design optimization process. Therefore, the information of optimal design and material volume is exactly reflected to fabricate the specimen for experimental validation. Through the design optimization examples of elasticity problem, it is experimentally shown that the optimal design has higher stiffness than the initial design. Also, the experimental results match very well with the numerical results. Using a non-contact optical 3D deformation measuring system for strain distribution, it is shown that the stress concentration is significantly alleviated in the optimal design compared with the initial design.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.1-7
• /
• 2011

In this paper, based on an isogeometric approach, we have developed a shape design optimization method for plane elasticity problems subjected to design-dependent loads. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry. In an isogeometric analysis, however, the geometric properties are already embedded in the B-spline basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. The solution space for the response analysis can be represented in terms of the same NURBS basis functions to represent the geometry, which yields a precise analysis model that exactly represents the normal and curvature depending on the applied loads. A continuum-based isogeometric adjoint sensitivity is extensively derived for the plane elasticity problems under the design-dependent loads. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Steel and Composite Structures
• /
• v.32 no.6
• /
• pp.821-832
• /
• 2019

An analysis on thermal buckling and postbuckling of composite laminated plates reinforced with a low amount of graphene platelets is performed in the current investigation. It is assumed that graphaene platelets are randomly oriented and uniformly dispersed in each layer of the composite media. Elastic properties of the nanocomposite media are obtained by means of the modified Halpin-Tsai approach which takes into account the size effects of the graphene reinforcements. By means of the von $K{\acute{a}}rm{\acute{a}}n$ type of geometrical nonlinearity, third order shear deformation theory and nonuniform rational B-spline (NURBS) based isogeometric finite element method, the governing equations for the thermal postbuckling of nanocomposite plates in rectangular shape are established. These equations are solved by means of a direct displacement control strategy. Numerical examples are given to study the effects of boundary conditions, weight fraction of graphene platelets and distribution pattern of graphene platelets. It is shown that, with introduction of a small amount of graphene platelets into the matrix of the composite media, the critical buckling temperature of the plate may be enhanced and thermal postbuckling deflection may be alleviated.

• Steel and Composite Structures
• /
• v.34 no.3
• /
• pp.361-376
• /
• 2020

This paper presents a free vibration analysis of shell panels made of functionally graded material (FGM) in the form of the ordinary and sandwich FGM and laminated shells using the isogeometric B3-spline finite strip method (IG-SFSM). B3-spline and Lagrangian interpolation are employed along the longitudinal and transverse directions respectively in this type of finite strip. The introduced finite strip formulation is based on the degenerated shell method, which provides variable thickness, arbitrary geometries, and analysis of thin or thick shells. Validity of the obtained natural frequencies by IG-SFSM is checked by comparison with results extracted from references for similar cases in different examples. These examples incorporate several geometries, materials, boundary conditions, and continuous thickness variation. A comparison of these two kinds of results and their proximity showed that the introduced IG-SFSM is a reliable tool which can be used in analysis of shells with the aforementioned properties.

• Structural Engineering and Mechanics
• /
• v.86 no.1
• /
• pp.93-107
• /
• 2023

In this paper, a density based topology optimization is proposed for generating of supports required in additive manufacturing to maintain the overhanging regions of main structures during layer by layer fabrication process. For this purpose, isogeometric analysis method is employed to model geometry and structural analysis of main and support structures. In order to model the problem two cases are investigated. In the first case, design domain of supports can easily be separated from the main structure by using distinct isogeometric patches. The second case happens when the main structure itself is optimized by using topology optimization and the supports should be designed in the voids of optimum layout. In this case, in order to avoid boundary identification and re-meshing process for separating design domain of supports from main structure, a parameterization technique is proposed to identify the design domain of supports. To achieve this, two density functions are defined over the entire domain to describe the main structure and supporting areas. On the other hand, since supports are under gravity loads while main structure and its stiffness is not completed during manufacturing process, in the proposed method, stiffness of the main structure is considered to be trivial and the gravity loads are also naturally applied to design support structures. By doing so, the results show reasonable supports are created to protect, continuously, overhanging surfaces of the main structure. Several examples are presented to demonstrate the efficiency of the proposed method and compare the results with literature.

• Advances in nano research
• /
• v.12 no.5
• /
• pp.441-455
• /
• 2022

This study explores the linear and nonlinear solutions of sigmoid functionally graded material (S-FGM) nanoplate with porous effects. A size-dependent numerical solution is established using the strain gradient theory and isogeometric finite element formulation. The nonlinear nonlocal strain gradient is developed based on the Reissner-Mindlin plate theory and the Von-Karman strain assumption. The sigmoid function is utilized to modify the classical functionally graded material to ensure the constituent volume distribution. Two different patterns of porosity distribution are investigated, viz. pattern A and pattern B, in which the porosities are symmetric and asymmetric varied across the plate's thickness, respectively. The nonlinear finite element governing equations are established for bending analysis of S-FGM nanoplates, and the Newton-Raphson iteration technique is derived from the nonlinear responses. The isogeometric finite element method is the most suitable numerical method because it can satisfy a higher-order derivative requirement of the nonlocal strain gradient theory. Several numerical results are presented to investigate the influences of porosity distributions, power indexes, aspect ratios, nonlocal and strain gradient parameters on the porous S-FGM nanoplate's linear and nonlinear bending responses.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.49-55
• /
• 2011

Turbulent flows with wavy wall are simulated using Residual-based Variational Multiscale Method (RB-VMS) which is proposed by Bazilves et al(2007) as new Large Eddy Simulation methodology. Incompressible Navier-Stokes equations are integrated using Isogeometric analysis which adopt the basis function as NURBS. The Reynolds number is 6760 based on the bulk velocity and averaged channel height. And the amplitude (${\alpha}/{\lambda}$) of wavy wall is 0.05. The computational domain is $2{\lambda}{\times}1.05{\lambda}{\times}{\lambda}$ in the streamwise, wall normal and spanwise direction. Mean quantities and turbulent statistics near wavy wall are compared with DNS results of Cherukat et al.(1998). The predicted results show good agreement with reference data.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.5
• /
• pp.287-292
• /
• 2021

We utilized the isogeometric analysis (IGA) method that uses NURBS basis functions in CAD systems, to account for the geometric exactness of a geometrically exact beam deformation, on a new type of metamaterial, twist-translation coupled structure showing a large twist angle. A two-dimensional unit cell structure was embedded in a cylindrical wall, using free-form deformation and an appropriate interpolation scheme. A parametric study on the effects of the dimensions of the cylinder and the number of cells, on the twisting angle was performed. Furthermore, the mechanism of the twist-translation coupled metamaterial was explored through numerical examples.