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

• Steel and Composite Structures
• /
• v.30 no.6
• /
• pp.517-534
• /
• 2019

In this paper, the effects of inevitable out-of-plane defects on the postbuckling behavior of single-layered graphene sheets (SLGSs) under in-plane loadings are investigated based on nonlocal first order shear deformation theory (FSDT) and von-Karman nonlinear model. A generic imperfection function, which takes the form of the products of hyperbolic and trigonometric functions, is employed to model out-of-plane defects as initial geometrical imperfections of SLGSs. Nonlinear equilibrium equations are derived from the principle of virtual work and variational formulation. The postbuckling equilibrium paths of imperfect graphene sheets (GSs) are presented by solving the governing equations via isogeometric analysis (IGA) and Newton-Raphson iterative method. Finally, the sensitivity of the postbuckling behavior of GS to shape, amplitude, extension on the surface, and location of initial imperfection is studied. Results showed that the small scale and initial imperfection effects on the postbuckling behavior of defective SLGS are important and cannot be ignored.

• 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.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 computational fluids engineering
• /
• v.15 no.2
• /
• pp.35-40
• /
• 2010

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number, 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorrector algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.56-59
• /
• 2009

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number is 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorretor algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.2
• /
• pp.127-134
• /
• 2009

T-splines are recently proposed geometric modeling tools. A T-spline surface is a NURBS surface with T-junctions and is defined by a control grid called T-mesh. Local refinement can be performed very easily for T-splines while it is limited for B-splines or NURBS. Using T-splines, patches with unmatched boundaries can be combined easily without special technique. In this study, the analysis methodology using T-splines is proposed. In this methodology, T-splines are used both for description of geometries and for approximation of solution spaces. Two-dimensional linear elastic and dynamic problems will be solved by employing the proposed T-spline finite element method, and the effectiveness of the current analysis methodology will be verified.

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

• 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.5
• /
• pp.337-343
• /
• 2014

This paper presents a continuum-based shape design sensitivity analysis(DSA) method for crack propagation problems using a reproducing kernel method(RKM), which facilitates the remeshing problem required for finite element analysis(FEA) and provides the higher order shape functions by increasing the continuity of the kernel functions. A linear elasticity is considered to obtain the required stress field around the crack tip for the evaluation of J-integral. The sensitivity of displacement field and stress intensity factor(SIF) with respect to shape design variables are derived using a material derivative approach. For efficient computation of design sensitivity, an adjoint variable method is employed tather than the direct differentiation method. Through numerical examples, The mesh-free and the DSA methods show excellent agreement with finite difference results. The DSA results are further extended to a shape optimization of crack propagation problems to control the propagation path.