• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.8
• /
• pp.175-184
• /
• 2003

A method of the shape design sensitivity analysis based on the boundary integral equation formulation is presented for two-dimensional inhomogeneous thermal conducting solids with multiple domains. Shape variation of the external and interface boundary is considered. A sensitivity formula of a general performance functional is derived by taking the material derivative to the boundary integral identity and by introducing an adjoint system. In numerical analysis, state variables of the primal and adjoint systems are solved by the boundary element method using quadratic elements. Two numerical examples of a compound cylinder and a thermal diffuser are taken to show implementation of the shape design sensitivity analysis. Accuracy of the present method is verified by comparing analyzed sensitivities with those by the finite difference. As application to the shape optimization, an optimal shape of the thermal diffuser is found by incorporating the sensitivity analysis algorithm in an optimization program.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.216-221
• /
• 2008

Shape design optimization for linear elasticity problem is performed using isogeometric analysis method. In many design optimization problems for real engineering models, initial raw data usually comes from CAD modeler. Then designer should convert this CAD data into finite element mesh data because conventional design optimization tools are generally based on finite element analysis. During this conversion there is some numerical error due to a geometry approximation, which causes accuracy problems in not only response analysis but also design sensitivity analysis. As a remedy of this phenomenon, the isogeometric analysis method is one of the promising approaches of shape design optimization. The main idea of isogeometric analysis is that the basis functions used in analysis is exactly same as ones which represent the geometry, and this geometrically exact model can be used shape sensitivity analysis and design optimization as well. In shape design sensitivity point of view, precise shape sensitivity is very essential for gradient-based optimization. In conventional finite element based optimization, higher order information such as normal vector and curvature term is inaccurate or even missing due to the use of linear interpolation functions. On the other hands, B-spline basis functions have sufficient continuity and their derivatives are smooth enough. Therefore normal vector and curvature terms can be exactly evaluated, which eventually yields precise optimal shapes. In this article, isogeometric analysis method is utilized for the shape design optimization. By virtue of B-spline basis function, an exact geometry can be handled without finite element meshes. Moreover, initial CAD data are used throughout the optimization process, including response analysis, shape sensitivity analysis, design parameterization and shape optimization, without subsequent communication with CAD description.

• Korean Journal of Mathematics
• /
• v.25 no.3
• /
• pp.405-435
• /
• 2017

We deal with a sensitivity analysis of an optimal shape control problem for the stationary Navier-Stokes system. A two-dimensional channel flow of an incompressible, viscous fluid is examined to determine the shape of a bump on a part of the boundary that minimizes the viscous drag. By using the material derivative method and adjoint variables for a shape sensitivity analysis, we derive the shape gradient of the design functional for the model problem.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.9
• /
• pp.1320-1327
• /
• 2004

An efficient boundary-based technique is developed for addressing shape design sensitivity analysis in supercavitating flow problem. An analytical sensitivity formula in the form of a boundary integral is derived based on the continuum formulation for a general functional defined in potential flow problems. The formula, which is expressed in terms of the boundary solutions and shape variation vectors, can be conveniently used for gradient computation in a variety of shape design in potential flow problems. While the sensitivity can be calculated independent of the analysis means, such as the finite element method (FEM) or the boundary element method (BEM), the FEM is used for the analysis in this study because of its popularity and easy-to-use features. The advantage of using a boundary-based method is that the shape variation vectors are needed only on the boundary, not over the whole domain. The boundary shape variation vectors are conveniently computed by using finite perturbations of the shape geometry instead of complex analytical differentiation of the geometry functions. The supercavitating flow problem is chosen to illustrate the efficiency of the proposed methodology. Implementation issues for the sensitivity analysis and optimization procedure are also addressed in this flow problem.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.05a
• /
• pp.90-95
• /
• 2001

A continuum-based shape design sensitivity analysis (DSA) method is presented for the acoustic radiation from thin body. The normal derivative integral formulation is employed as an analysis formulation and differentiated directly by using material derivative to get the acoustic shape design sensitivity. In the acoustic sensitivity formulation, derivative coefficients of the structural normal velocities on the surface are required as the input. Thus, the shape design sensitivities of structural velocities on the surface with respect to the shape change are also calculated with continuum approach. A simple disk is considered as a numerical example to validate the accuracy and efficiency of the analytical shape design sensitivity equations derived in this research. This research should be very helpful to design an application involving thin body and to change its acoustic characteristics.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.577-582
• /
• 2007

A variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions in analysis domain arc generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Thus. the solution space can be represented in terms of the same functions to represent the geometry. The coefficients of basis functions or the control variables play the role of degrees-of-freedom. Furthermore, due to h-. p-, and k-refinement schemes, the high order geometric features can be described exactly and easily without tedious re-meshing process. The isogeometric sensitivity analysis method enables us to analyze arbitrarily shaped structures without re-meshing. Also, it provides a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling. To obtain precise shape sensitivity, the normal and curvature of boundary should be taken into account in the shape sensitivity expressions. However, in conventional finite element methods, the normal information is inaccurate and the curvature is generally missing due to the use of linear interpolation functions. A continuum-based adjoint sensitivity analysis method using the isogeometric approach is derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of boundary. In isogeometric analysis, however, the geometric properties arc already embedded in the B-spline shape functions and control points. The perturbation of control points in isogeometric analysis automatically results in shape changes. Using the conventional finite clement method, the inter-element continuity of the design space is not guaranteed so that the normal vector and curvature arc not accurate enough. On tile other hand, in isogeometric analysis, these values arc continuous over the whole design space so that accurate shape sensitivity can be obtained. Through numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.335-342
• /
• 2002

A continuum-based design sensitivity analysis (DSA) method fur non-shape problems is developed for geometrically nonlinear elastic structures. The non-shape problem is characterized by the design variables that are not associated with the domain of system like sizing, material property, loading, and so on. Total Lagrangian formulation with the Green-Lagrange strain and the second Piola-Kirchhoff stress is employed to describe the geometrically nonlinear structures. The spatial domain is discretized using the 4-node isoparametric plane stress/strain elements. The resulting nonlinear system is solved using the Newton-Raphson iterative method. To take advantage of the derived analytical sensitivity In topology optimization, a fast and efficient design sensitivity analysis method, adjoint variable method, is employed and the material property of each element is selected as non-shape design variable. Combining the design sensitivity analysis method and a gradient-based design optimization algorithm, an automated design optimization method is developed. The comparison of the analytical sensitivity with the finite difference results shows excellent agreement. Also application to the topology design optimization problem suggests a very good insight for the layout design.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.5
• /
• pp.1458-1467
• /
• 1996

A direct differentiationmethod is presented for the shape design sensitivity analysis of axisymmeetric elastic solids. Based on the exisymmetric boundary integralequaiton formulation, a new boundary ntegral equatio for sensitivity analysis is derived by taking meterial derivative to the same integral identity that was used in the adjoint variable melthod. Numerical implementation is performed to show the applicaiton of the theoretical formulation. For a simple example with analytic solution, the sensitivities by present method are compared with analytic sensitivities. As an application to the shape optimization, an optimal shape of a gas turbine disc toinimize the weight under stress constraints is found by incorporating the sensitivity analysis algorithm in an optimizatio program.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.52 no.7
• /
• pp.307-314
• /
• 2003

This paper presents a 3D shape optimization algorithm for electromagnetic devices using the design sensitivity analysis with finite element method. The structural deformation analysis based on the deformation theory of the elastic body under stress is used for mesh renewing. The design sensitivity and adjoint variable formulae are derived for the 3D finite element method with edge element. The results of sensitivity analysis are used as the input data of the structural analysis to calculate the relocation of the nodal points. This method makes it possible that the new mesh of analysis region can be obtained from the initial mesh without regeneration. The proposed algorithm is applied to the shape optimization of 3D electromagnet pole to net a uniform flux density at the target region.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1047-1052
• /
• 2004

An efficient boundary-based technique is developed for addressing shape design sensitivity analysis in supercavitating flow problem. An analytical sensitivity formula in the form of a boundary integral is derived based on the continuum formulation for a general functional defined in potential flow problems. The formula, which is expressed in terms of the boundary solutions and shape variation vectors, can be conveniently used for gradient computation in a variety of shape design in potential flow problems. While the sensitivity can be calculated independent of the analysis means, such as the finite element method (FEM) or the boundary element method (BEM), the FEM is used for the analysis in this study because of its popularity and easy-touse features. The advantage of using a boundary-based method is that the shape variation vectors are needed only on the boundary, not over the whole domain. The boundary shape variation vectors are conveniently computed by using finite perturbations of the shape geometry instead of complex analytical differentiation of the geometry functions. The supercavitating flow problem is chosen to illustrate the efficiency of the proposed methodology. Implementation issues for and optimization procedure are addressed in this flow problem.