• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.1
• /
• pp.79-87
• /
• 2013

An isogeometric topological shape optimization method is developed using the level sets and Heaviside enrichments. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set functions, which facilitates to handle complicated topological shape changes. The Heaviside enrichment improves the isogeometric analysis by adding some enrichment functions to model the internal boundaries. The proposed topological shape optimization method has several benefits: exact geometric models can be obtained using the isogeometric approach and the limitation of tensor-product patches can be overcome using the Heaviside enrichments to represent the internal voids. Even in a single patch, discontinuous displacement fields as well as smooth stress field can be obtained. Since the level sets offer the implicit moving boundary inside the domain, it is easy to represent the topological shape variations in the isogeometric analysis using Heaviside enrichments.

• Structural Engineering and Mechanics
• /
• v.59 no.3
• /
• pp.503-526
• /
• 2016

In-plane and out-of-plane free vibration analysis of Timoshenko curved beams is addressed based on the isogeometric method, and an effective scheme to avoid numerical locking in both of the two patterns is proposed in this paper. The isogeometric computational model takes into account the effects of shear deformation, rotary inertia and axis extensibility of curved beams, and is applicable for uniform circular beams, and more complicated variable curvature and cross-section beams as illustrated by numerical examples. Meanwhile, it is shown that, the $C^{p-1}$-continuous NURBS elements remarkably have higher accuracy than the finite elements with the same number of degrees of freedom. Nevertheless, for in-plane or out-of-plane vibration analysis of Timoshenko curved beams, the NURBS-based isogeometric method also exhibits locking effect to some extent. To eliminate numerical locking, the selective reduced one-point integration and $\bar{B}$ projection element based on stiffness ratio is devised to achieve locking free analysis for in-plane and out-of-plane models, respectively. The suggested integral schemes for moderately slender models obtain accurate results in both dominated and non-dominated regions of locking effect. Moreover, this strategy is effective for beam structures with different slenderness. Finally, the influence factors of structural parameters of curved beams on their natural frequency are scrutinized.

• Journal of Korean Association for Spatial Structures
• /
• v.12 no.2
• /
• pp.45-53
• /
• 2012

This paper derives Isogeometric analysis and post-processing method of surface that are generated by NURBS basis function for accurate geometric modeling and structure analysis of free-form. By deforming these parameters that are consisted of control points, knots, polynomial, variable geometric models are derived. The basis function that is used to Isogeometric analysis is same to the basis function of NURBS that is used to generate geometric models. For performing isogeometric analysis, h-p-k refinement is performed without changing of original geometry. To visualize the results of isogeometric analysis that control points' displacements, post-processing method that is the interface method between IGES format and Rhinoceros is derived.

• 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.

• Computers and Concrete
• /
• v.23 no.3
• /
• pp.171-188
• /
• 2019

This study proposed a new and efficient 2D damage-plasticity model within the framework of Isogeometric analysis (IGA) for the geometrically nonlinear damage analysis of concrete. Since concrete exhibits complicated material properties, two internal variables are introduced to measure the hardening/softening behavior of concrete in tension and compression, and an implicit gradient-enhanced formulation is adopted to restore the well-posedness of the boundary value problem. The numerical results calculated by the model is compared with the experimental data of three benchmark problems of plain concrete (three-point and four-point bending single-notched beams and four-point bending double-notched beam) to illustrate the geometrical flexibility, accuracy, and robustness of the proposed approach. In addition, the influence of the characteristic length on the numerical results of each problem is investigated.

• 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.23 no.3
• /
• pp.275-281
• /
• 2010

In this paper, the design optimization of structures with stress constraints is performed using isogeometric shape optimization method. The stress constraints have an important role in design optimization problems since stress concentration could result in structural failure. To represent exact geometry in analysis, the isogeometric analysis method uses the same basis functions as used in the CAD geometry. The geometrically exact model can be used in both stress and design sensitivity analyses so that it can yield more precise optimal design than finite element one. Through numerical examples, the isogeometric approach turns out to be effective in shape optimization problems under stress constraints.

• Structural Engineering and Mechanics
• /
• v.41 no.5
• /
• pp.617-632
• /
• 2012

An isogeometric approach is presented for static analysis of thin plate problems of various geometries. Non-Uniform Rational B-Splines (NURBS) basis function is applied for approximation of the thin plate deflection, as for description of the geometry. The governing equation based on Kirchhoff plate theory, is discretized using the standard Galerkin method. The essential boundary conditions are enforced by the Lagrange multiplier method. Several typical examples of thin plate and thin plate on elastic foundation are solved and compared with the theoretical solutions and other numerical methods. The numerical results show the robustness and efficiency of the proposed approach.

• Steel and Composite Structures
• /
• v.21 no.6
• /
• pp.1389-1410
• /
• 2016

This paper proposes a hybrid of topological derivative-based level set method (LSM) and isogeometric analysis (IGA) for structural topology optimization. In topology optimization a significant drawback of the conventional LSM is that it cannot create new holes in the design domain. In this study, the topological derivative approach is used to create new holes in appropriate places of the design domain, and alleviate the strong dependency of the optimal topology on the initial design. Furthermore, the values of the gradient vector in Hamilton-Jacobi equation in the conventional LSM are replaced with a Delta function. In the topology optimization procedure IGA based on Non-Uniform Rational B-Spline (NURBS) functions is utilized to overcome the drawbacks in the conventional finite element method (FEM) based topology optimization approaches. Several numerical examples are provided to confirm the computational efficiency and robustness of the proposed method in comparison with derivative-based LSM and FEM.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.6
• /
• pp.515-522
• /
• 2017

In this paper, an isogeometric analysis for multipatch problem is investigated, in which two or more geometries are connected at the interface in a conforming or non-conforming conditions. To express higher continuity at the patch interface, two approaches such as Nitsche based method and master-slave method are formulated for the linear elasticity problem and discretized using the isogeometric approach using NURBS basis functions. A short comparison between two approaches in formulations reveals the pros and cons of them with the applicability in the isogeometric multipatch problem. In addition, a NURBS based stress recovery is adopted to express a better stress continuity through the post-processing. Numerical examples indicate the effectiveness of Nitsche method in the non-conforming patch, following the exact solution well. For the stress concentration problem with the conforming patch, introduced two methodologies show comparative results, meanwhile the NURBS based stress recovery presents an improved smooth stress contour in the whole domain including the patch interface.