• 한국정밀공학회지
• /
• 제20권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.

• 한국전자파학회:학술대회논문집
• /
• 한국전자파학회 2005년도 종합학술발표회 논문집 Vol.15 No.1
• /
• pp.169-174
• /
• 2005

본 논문에서는 시간영역 유한차분법(FDTD: Finite-Difference Time-Domain Method)과 설계민감도법(Design Sensitivity Analysis)을 이용하여 유전체 산란체(Dielectric Scatterer)를 복원하기 위한 역산란문제(Inverse Scattering의 새로운 해석기법을 제안하였다. 이때 복원의 빠른 수렴을 위하여 도함수를 이용한 설계민감도법을 도입하였고 본 연구에서는 시간영역 유한차분법으로부터 직접 설계민감도 수식을 도출하였다. 계산의 효율성을 위하여 보조변수법(Adjoint Variable Method)을 도입하여 보조변수 방정식을 도출하고 최적화 알고리듬으로 최대경사도법을 이용하여 반복적인 추정을 통하여 유전체를 복원하였다. 본 연구의 타당성의 보이기 위하여 2차원 $TM^2$에서의 유전체 복원 사례를 제시한다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1053-1057
• /
• 2004

This paper presents the adjoint variable design sensitivity analysis for thermal systems considering both conduction and convection heat transfer. Both nodal temperature and total heat flow are considered to be objective functions and design sensitivity formulas are derived for each case. For the case of convection heat transfer, the adjoint analysis is carefully proceeded to obtain a precise result. A topology optimization example is examined for a simple planar square plate in order to design a heat exchanger as verification.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2004년도 봄 학술발표회 논문집
• /
• pp.127-134
• /
• 2004

In this paper, using an adjoint variable method, we develop a design sensitivity analysis (DSA) method applicable to heat conduction problems in steady state. Also, a topology design optimization method is developed using the developed DSA method. Design sensitivity expressions with respect to the thermal conductivity are derived. Since the already factorized system matrix is utilized to obtain the adjoint solution, the cost for the sensitivity computation is trivial. For the topology design optimization, the design variables are parameterized into normalized bulk material densities. The objective function and constraint are the thermal compliance of structures and allowable material volume, respectively. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with finite difference ones, requiring less than 0.3% of CPU time far the finite differencing. Also, the topology optimization yields physical meaningful results.

• 한국전산유체공학회지
• /
• 제11권1호
• /
• pp.57-66
• /
• 2006

An optimal shape design approach is presented for a subsonic S-shaped intake using aerodynamic sensitivity analysis. Two-equation turbulence model is employed to capture strong counter vortices in the S-shaped duct more precisely. Sensitivity analysis is performed for the three-dimensional Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method For code validation, the result of the flow solver is compared with experiment data and other computational results of bench marking test. To study the influence oj turbulence models and grid refinement on the duct flow analysis, the results from several turbulence models are compared with one another and the minimum number of grid points, which can yield an accurate solution is investigated The adjoint variable code is validated by comparing the complex step derivative results. To realize a sufficient and flexible design space, NURBS equations are introduced as a geometric representation and a new grid modification technique, Least Square NURBS Grid Approximation is applied With the verified flow solver, the sensitivity analysis code and the geometric modification technique, the optimization of S-shaped intake is carried out and the enhancement of overall intake performance is achieved The designed S-shaped duct is tested in several off-design conditions to confirm the robustness of the current design approach. As a result, the capability and the efficiency of the present design tools are successfully demonstrated in three-dimensional highly turbulent internal flow design and off-design conditions.

• 수산해양기술연구
• /
• 제35권4호
• /
• pp.435-447
• /
• 1999

It is main objective of this approach to present a method to analyse stochastic design sensitivity for problems of structural dynamics with randomness in design parameters. A combination of the adjoint variable approach and the second order perturbation method is used in the finite element approach. An alternative form of the constant functional that holds for all times is introduced to consider the time response of dynamic sensitivity. The terminal problem of the adjoint system is solved using equivalent homogeneous equations excited by initial velocities. The numerical procedures are shown to be much more efficient when based on the fold superposition method: the generalized co-ordinates are normalized and the correlated random variables are transformed to uncorrelated variables, whereas the secularities are eliminated by the fast Fourier transform of complex valued sequences. Numerical algorithms have been worked out and proved to be accurate and efficient : they can be readily adapted to fit into the existing finite element codes whose element derivative matrices can be explicitly generated. The numerical results of two cases -2 dimensional portal frame for the comparison with reference and 3-dimensional frame structure - for the deterministic sensitivity analysis are presented.

• 한국전산구조공학회논문집
• /
• 제27권4호
• /
• pp.297-303
• /
• 2014

페리다이나믹스 이론과 이진분해 기법의 병렬연산을 이용하여 동적 균열진전 문제에 대한 애조인 형상 설계민감도 해석법을 개발하였다. 페리다이나믹스에서는 균열의 연속적인 분기를 다룰 수 있으며, Explicit 시간적분법을 채택한다. 설계민감도 해석은 애조인 변수법은 경로의존성 문제에는 적합하지 않으나 여기서는 응답해석의 경로를 이미 알고 있으므로 채택하여 사용할 수 있었다. 얻어진 해석적 설계민감도는 유한차분과 비교하여 그 정확성을 검증하였다. 유한차분법은 설계섭동량에 민감하여 비선형성이 강한 페리다이나믹스 문제에서 부정확한 설계민감도를 제시할 수 있다. 정확한 설계민감도 해석을 위해서는 이산화과정에서 $C^1$ 연속성을 가지는 체적율이 필요함을 알 수 있었다.

• Journal of Mechanical Science and Technology
• /
• 제14권2호
• /
• pp.197-206
• /
• 2000

A boundary integral equation method in the shape design sensitivity analysis is developed for the elasticity problems with axisymmetric non-homogeneous bodies. Functionals involving displacements and tractions at the zonal interface are considered. Sensitivity formula in terms of the interface shape variation is then derived by taking derivative of the boundary integral identity. Adjoint problem is defined such that displacement and traction discontinuity is imposed at the interface. Analytic example for a compound cylinder is taken to show the validity of the derived sensitivity formula. In the numerical implementation, solutions at the interface for the primal and adjoint system are used for the sensitivity. While the BEM is a natural tool for the solution, more generalization should be made since it should handle the jump conditions at the interface. Accuracy of the sensitivity is evaluated numerically by the same compound cylinder problem. The endosseous implant-bone interface problem is considered next as a practical application, in which the stress value is of great importance for successful osseointegration at the interface. As a preliminary step, a simple model with tapered cylinder is considered in this paper. Numerical accuracy is shown to be excellent which promises that the method can be used as an efficient and reliable tool in the optimization procedure for the implant design. Though only the axisymmetric problem is considered here, the method can be applied to general elasticity problems having interface.

• 대한기계학회논문집A
• /
• 제24권3호
• /
• pp.645-652
• /
• 2000

A method for shape design sensitivity analysis is proposed utilizing commercial software ANSYS for thermal conduction problems. While the sensitivity formula is derived analytically by introduing adjoint variable concept, sensitivity calculation in practice as well as the primal and adjoint solution of thermal conduction is performed using the ANSYS very easily. Since the formula always takes boundary integral form, sensitivity evaluation in ANSYS requires a little more addition of post-processing routine which involves evaluation of boundary variable from the obtained solution. Though the BEM has been used as a better tool for this purpose, the present study shows it can also be calculated using any kind of analysis code such as ANSYS since the formula is based on analytic nature. Therefore the present study provides a new and efficient way of optimization which was not possible before using commercial software. The usefulness of the method is illustrated via a weight minimization problem of thermal diffuser.

• 한국항공우주학회지
• /
• 제33권1호
• /
• pp.1-10
• /
• 2005

비정렬 격자계에서 continuous adjoint 방정식을 사용하여 제자리 비행을 하는 헬리콥터 로터 블레이드에 대한 공력 형상 최적설계 기법을 개발하였다. 효율적인 민감도 계산을 위해 회전좌표계에서 continuous adjoint 민감도 해석 기법을 유도하였다. 설계과정의 반복적인 수치계산의 효율을 높이기 위해서 영역 분할 기법에 기반을 둔 병렬처리 기법을 도입하였다. 끝단 와류의 정확한 포착을 위해서 끝단와류를 따른 격자적응을 수행하였다. 이러한 방법은 Caradonna와 Tung의 실험형상 및 UH60 헬리콥터 로터 블레이드의 공력 최적설계에 적용되었으며, 본 연구에서 사용된 최적설계 기법을 이용하면 일정한 추력을 유지하면서 요구동력을 현저하게 줄일 수 있음을 보였다.