• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.339-345
• /
• 2007

In this paper, 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 for response analysis are generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Furthermore, the solution space for the response analysis can be represented in terms of the same functions to represent the geometry, which enables to provide a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling and analyze arbitrarily shaped structures without re-meshing. In this paper, a continuum-based adjoint sensitivity analysis method using the isogeometric approach is extensively derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry In the 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. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.14 no.1
• /
• pp.94-101
• /
• 2000

In this paper, a multi-channel active noise control(ANC) system using subband hybrid control techniques is proposed. Subband techniques could reduce computational burden and improve the performance of ANC systems by dividing several frequency subband and adjusting adaptive filter coefficients. So it can effectively cancel noises at wanted frequency range and use lower order adaptive filter than the existing algorithms. The adjoint LMS algorithm, which prefilter the error signals instead of the divided reference signals in frequency band, is also used for adaptive filter algorithms to reduce the computational burden of the subband adaptive systems. To improve performance of the ANC system, a weighted hybrid control technique, which has weightily properties of feedforward control systems and feedback control systems, is applied. This algorithm shows higher stability and good noise attenuation property in broad band ANC systems. Computer simulations were performed to show the effectiveness of the proposed algorithm.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.1
• /
• pp.82-87
• /
• 1987

In this study a method of shape optimization is applied to two dimensional heat transfer system. For this the optimization problem is defined in a functional form including cost, constraints and the system governing equation. Then the material derivative concept in continuum mechanics and the adjoint variable method are employed for the shape design sensitivity analysis. With the sensitivity analysis results, an optimum is sought with the gradient projection optimization algorithm. The two dimensional isoparametric finite elements are used for accurate analysis and sensitivity calculations. The above method is employed to find the boundary shape to achieve a desired temperature distribution along a segment of the boundary subject to the maximum area constraint.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.7
• /
• pp.821-827
• /
• 2010

This paper proposes a density method based topology optimization of a perpendicular magnetic recording system design in which the saturation effect is taken into account. During the topology optimization process in magnetic fields, the magnetic reluctivity is updated in accordance with the changes in element density determined by a sensitivity analysis. The magnetic reluctivity is determined from a B-H curve and is used to represent nonlinear material property, i.e., the saturation effect. The sensitivity for a generalized response functional is formulated using the adjoint variable method in which the nonlinear property is taken into account and the objective function is set such that the magnetic energy in the media is maximized. Effects due to the nonlinear property can be observed from a numerical study in which the linear and the nonlinear topology optimization results are compared.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.20 no.8
• /
• pp.1487-1493
• /
• 2016

The susceptible, exposed, infectious, and recovered model (SEIR) is used extensively in the field of epidemiology. On the other hand, dissemination information among users through internet grows exponentially. This information spreading can be modeled as an epidemic. In this paper, we derive the mathematical model of SEIR in viral communication from the view of optimal control theory. Overall the methods based on classical calculus, In order to solve the optimal control problem, proved to be more efficient and accurate. According to Pontryagin's minimum principle (PMP) the Hamiltonian function must be optimized by the control variables at all points along the solution trajectory. We present our method based on the PMP and forward backward algorithm. In this algorithm, one should integrate forward in time for the state equations then integrate backward in time for the adjoint equations resulting from the optimality conditions. The problem is mathematically analyzed and numerically solved as well.

• 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.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.1
• /
• pp.37-43
• /
• 2010

In this paper, a reliability-based design optimization is developed for the topology design of linear structures using a performance measure approach. Spatial domain is discretized using three dimensional Reissner-Mindlin plate elements and design variable is taken as the material property of each element. A continuum based adjoint variable method is employed for the efficient computation of sensitivity with respect to the design and random variables. The performance measure approach of RBDO is employed to evaluate the probabilistic constraints. The topology optimizationproblem is formulated to have probabilistic displacement constraints. The uncertainties such as material property and external loads are considered. Numerical examples show that the developed topology optimization method could effectively yield a reliable design, comparing with the other methods such as deterministic, safety factor, and worst case approaches.

• Bulletin of the Korean Mathematical Society
• /
• v.30 no.1
• /
• pp.71-77
• /
• 1993

Every ring considered in the paper will be assumed to be commutative and have a unit element. An ideal A of a ring R will be called primal if the elements of R which are zero divisors modulo A, form an ideal of R, say pp. If A is a primal ideal of R, P is called the adjoint ideal of A. The adjoint ideal of a primal ideal is prime [2]. The definition of primal ideals may also be formulated as follows: An ideal A of a ring R is primal if in the residue class ring R/A the zero divisors form an ideal of R/A. If Q is a primary idel of a ring R then every zero divisor of R/Q is nilpotent; therefore, Q is a primal ideal of R. That a primal ideal need not be primary, is shown by an example in [2]. Let R[X], and R[[X]] denote the polynomial ring and formal power series ring in an indeterminate X over a ring R, respectively. Let S be a multiplicative system in a ring R and S$^{-1}$ R the quotient ring of R. Let Q be a P-primary ideal of a ring R. Then Q[X] is a P[X]-primary ideal of R[X], and S$^{-1}$ Q is a S$^{-1}$ P-primary ideal of a ring S$^{-1}$ R if S.cap.P=.phi., and Q[[X]] is a P[[X]]-primary ideal of R[[X]] if R is Noetherian [1]. We search for analogous results when primary ideals are replaced with primal ideals. To show an ideal A of a ring R to be primal, it sufficies to show that a-b is a zero divisor modulo A whenever a and b are zero divisors modulo A.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• v.5B no.3
• /
• pp.233-237
• /
• 2005

The reduction of the ripple of driving force is especially required for the practical utilization of linear synchronous motor for rope-less elevator. In this paper, the magnetic region of the linear motor is optimized by using topology optimization techniques (density method and ON/OFF method) in order to reduce the ripple of driving force. The optimal results of both methods are compared, and useful information for the optimal design of linear motor is obtained.