• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.4
• /
• pp.425-431
• /
• 2015

Based on a bond-based peridynamics theory for dynamic crack propagation problems, this paper presents a design sensitivity analysis and optimization method. Peridynamics has a peculiar advantage over the existing continuum theory in the mathematical modelling of problems where discontinuities arise. For the design optimization of the crack propagation problems, a non-shape design sensitivity is derived using the adjoint variable method. The obtained adjoint sensitivity of displacement and strain energy turns out to be very accurate and efficient compared to the finite different sensitivity. The obtained design sensitivities are futher utilized to optimally control the position of bifurcation point in the design optimization of crack propagation in a plate under tension. A numerical experiment demonstrates that the optimal distribution of material density could delay the position of bifurcation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.5
• /
• pp.275-282
• /
• 2018

Using isogeometric analysis(IGA) gives more accurate results for higher order mode in eigenvalue problem than using the finite element method(FEM). This is because the FEM has $C^0$ continuity between elements, whereas IGA guarantee $C^{P-1}$ between elements for p-th order basis functions. In this paper, a mode based reduced model is constructed by using IGA and dynamic behavior analysis is performed using this advantage. Craig-Bampton(CB) method is applied to construct the reduced model. Several numerical examples were performed to compare the eigenvalue analysis results for various order of element basis function by applying the IGA and FEM to simple rod analysis. We have confirmed that numerical error increases in the higher order mode as the continuity between elements decreases in the IGA by allowing internal knots multiplicity. The accuracy of the solution can be improved by using the IGA with high inter-element continuity when high-frequency external force acts on the reduced model for dynamic behavior analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.441-449
• /
• 2015

In this paper we perform a linearized buckling analysis using the Kirchhoff plate theory and the von Karman nonlinear strain-displacement relation. Design sensitivity analysis(DSA) expressions for plane elasticity and buckling problems are derived with respect to Young's modulus and thickness. Using the design sensitivity, we can formulate the topology optimization method for minimizing the compliance and maximizing eigenvalues. We develop a topology optimization method applicable to plate buckling problems using the prestress for buckling analysis. Since the prestress is needed to assemble the stress matrix for buckling problem using the von Karman nonlinear strain, we introduced out-of-plane motion. The design variables are parameterized into normalized bulk material densities. The objective functions are the minimum compliance and the maximum eigenvalues and the constraint is the allowable volume. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with the finite difference ones and the topology optimization yields physically meaningful results.

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

This study is concerned with the shape optimization of structural members frequently found in critical area in a structure system, that is, highly stressed zone. Isogeometry analysis is well known to be the very efficient way to integrate the geometric modeling(CAD) and computational analysis(CAE). This can be accomplished by directly using the geometric modeling by NURBS(Non-Uniform Rational Basis Spline). In this study, an efficient computer code adopting the isogeometry concept has been developed for the structural analysis, in which CAD information can be directly used in the finite element modeling. In order to show the validity of the present code, the present results are compared with those by using the commercial package, that is, MSC/NASTRAN. The present isogeometric analysis procedure has been integrated with the optimization procedure to deal with the optimization problem found in the context of structural mechanics. The present system has been successfully applied to the shape optimization of cantilever structure having bracket. From the present study, it can be seen the validity of the present approach and computer codes developed in this study. This paper ends with some discussions about the practical usefulness of the present approach which is based on isogeometry analysis, and extension of the present study.

• Steel and Composite Structures
• /
• v.33 no.5
• /
• pp.717-727
• /
• 2019

This paper is motivated by the lack of studies in the technical literature concerning to vibration analysis of a single-layered graphene sheet (SLGS) with corner cutout based on the nonlocal elasticity model framework of classical Kirchhoff thin plate. An isogeometric analysis (IGA) based upon non-uniform rational B-spline (NURBS) is employed for approximation of the L-shape SLGS deflection field. Trimming technique is employed to create the cutout in geometry of L-shape plate. The L-shape plate is assumed to be Free (F) in the straight edges of cutout while any arbitrary boundary conditions are applied to the other four straight edges including Simply supported (S), Clamped (C) and Free (F). The Numerical studies are carried out to express the influences of the nonlocal parameter, cutout dimensions, boundary conditions and mode numbers on the variations of the natural frequencies of SLGS. It is precisely shown that these parameters have considerable effects on the free vibration behavior of the system. In addition, numerical results are validated and compared with those achieved using other analysis, where an excellent agreement is found. The effectiveness and the accuracy of the present IGA approach have been demonstrated and it is shown that the IGA is efficient, robust and accurate in terms of nanoplate problems. This study serves as a benchmark for assessing the validity of numerical methods used to analyze the single-layered graphene sheet with corner cutout.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.4
• /
• pp.297-303
• /
• 2014

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that $C^1$-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.4
• /
• pp.199-204
• /
• 2021

In this study, an isogoemetric analysis (IGA) method that uses NURBS (Non-Uniform Rational B-Spline) basis functions in computer-aided design (CAD) systems is employed to account for the geometric exactness of curved electrodes constituting an electro-adhesive pad in electrostatic problems. The IGA is advantageous for obtaining precise normal vectors when computing the electro-adhesive forces on curved surfaces. By performing parametric studies using numerical examples, we demonstrate the superior performance of the curved electrodes, which is attributed to the increase in the normal component of the electro-adhesive forces. In addition, concave curved electrodes exhibit better performance than their convex counterparts.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.1
• /
• pp.9-18
• /
• 2015

In this paper, we verify the optimal topology design for heat conduction problems in steady stated which is obtained numerically using the adjoint design sensitivity analysis(DSA) method. In adjoint variable method(AVM), the already factorized system matrix is utilized to obtain the adjoint solution so that its computation cost is trivial for the sensitivity. For the topology optimization, the design variables are parameterized into normalized bulk material densities. The objective function and constraint are the thermal compliance of the structure and the allowable volume, respectively. For the experimental validation of the optimal topology design, we compare the results with those that have identical volume but designed intuitively using a thermal imaging camera. To manufacture the optimal design, we apply a simple numerical method to convert it into point cloud data and perform CAD modeling using commercial reverse engineering software. Based on the CAD model, we manufacture the optimal topology design by CNC.

• Structural Engineering and Mechanics
• /
• v.30 no.2
• /
• pp.225-245
• /
• 2008

T-splines are recently proposed mathematical tools for geometric modeling, which are generalizations of B-splines. Local refinement can be performed effectively using T-splines while it is not the case when B-splines or NURBS are used. Using T-splines, patches with unmatched boundaries can be combined easily without special techniques. In the present study, an analysis framework using T-splines is proposed. In this framework, T-splines are used both for description of geometries and for approximation of solution spaces. This analysis framework can be a basis of a CAD/CAE integrated approach. In this approach, CAD models are directly imported as the analysis models without additional finite element modeling. Some numerical examples are presented to illustrate the effectiveness of the current analysis framework.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2014

Using the level set and the meshfree methods, we develop a topological shape optimization method applied to linear elasticity problems. Design gradients are computed using an efficient adjoint design sensitivity analysis(DSA) method. The boundaries are represented by an implicit moving boundary(IMB) embedded in the level set function obtainable from the "Hamilton-Jacobi type" equation with the "Up-wind scheme". Then, using the implicit function, explicit boundaries are generated to obtain the response and sensitivity of the structures. Global nodal shape function derived on a basis of the reproducing kernel(RK) method is employed to discretize the displacement field in the governing continuum equation. Thus, the material points can be located everywhere in the continuum domain, which enables to generate the explicit boundaries and leads to a precise design result. The developed method defines a Lagrangian functional for the constrained optimization. It minimizes the compliance, satisfying the constraint of allowable volume through the variations of boundary. During the optimization, the velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian functional. Compared with the conventional shape optimization method, the developed one can easily represent the topological shape variations.